### Author Topic: Negative energy in phonon band structure calculation  (Read 153 times)

0 Members and 1 Guest are viewing this topic.

#### Utpal Sarkar

• New ATK user
• Posts: 3
• Country:
• Reputation: 0
##### Negative energy in phonon band structure calculation
« on: September 1, 2020, 20:23 »
I am a beginner in quantum ATK and facing a problem during phonon band structure calculation I am getting negative energy states. Could you please suggest me how to overcome this?
I have re optimized the structure with high force and stress tolerance but unable to get the correct phonon band structure. I am attaching the inputs and outputs for your reference.

#### NW

• Heavy ATK user
• Posts: 82
• Country:
• Reputation: 1
##### Re: Negative energy in phonon band structure calculation
« Reply #1 on: September 3, 2020, 15:20 »
Hi,

I have one suggestion that may help you. Consider bigger cell for phonon calculation.

#### Utpal Sarkar

• New ATK user
• Posts: 3
• Country:
• Reputation: 0
##### Re: Negative energy in phonon band structure calculation
« Reply #2 on: September 3, 2020, 18:47 »
What  do you mean by bigger cell? whether you want to mean that I have to design a  super cell and re-optimize it and followed by phonon band structure calculation or you want to mean that I have to increase the size of the lattice vector of the unit cell and re-optimize it and followed by phonon band structure calculation. Kindly enlighten me on this topic.

#### NW

• Heavy ATK user
• Posts: 82
• Country:
• Reputation: 1
##### Re: Negative energy in phonon band structure calculation
« Reply #3 on: September 3, 2020, 19:03 »
I mean repeat the structure in x and y directions (2*2) and then re-optimize the cell and structure. Then calculate phonon bandstructure.

#### Utpal Sarkar

• New ATK user
• Posts: 3
• Country:
• Reputation: 0
##### Re: Negative energy in phonon band structure calculation
« Reply #4 on: September 13, 2020, 08:48 »
Sir as per your suggestion I have performed phonon band structure calculation for 2*2*1 super cell. But the system still have negative energy of -6.18 meV which is earlier -6.27 for unit cell. To further increase the super cell size is very difficult for computation. Could you please suggest any other way to solve this?
I am attaching here input and output for your reference.

#### Troels-Markussen

• Regular ATK user
• Posts: 30
• Country:
• Reputation: 1
• QuantumATK staff
##### Re: Negative energy in phonon band structure calculation
« Reply #5 on: September 14, 2020, 16:40 »
Hi,

I noticed that in your geometry optimization, you fixed the strain in the x- and y- directions, while not in the z-direction. I think this is a mistake. For such a 2D material you should rather fix the strain in the z-direction and relax it in the x- and y-directions like:

Code: [Select]
``constraints = [FixStrain(x=False, y=False, z=True)]``

Try to re-run the geometry optimization with this constraint and the recalculate the DynamicalMatrix.

Additionally, phonon calculations might require more accurate setting than 'normal' calculations. For this you might consider decreasing the tolerance as:

Code: [Select]
``iteration_control_parameters = IterationControlParameters(    tolerance=1e-05,    )``

Regards,
Troels