Hi,
I would like to consult how to relax the structure at finite electric field?
My main purpose is to calculate thermal properties of GaN under electric field. So, I want to relax the structure at finite electric field. I modified the example for structure optimization calculation of wurtzite structure to do my calculations, as shown in attached input file. Is that correct?
I ran this file. It works without error, but some warnings especially the following one besides berryphase_new warning:
chkinp : WARNING 
Constant unreduced E calculation with relaxation of cell parameters is allowed.
But we strongly recommend users to use reduced ebar calculation (berryopt=14)
with the relaxation of cell parameters, for internal consistency purpose.
For more information, please refer to "M. Stengel, N.A. Spaldin and D.Vanderbilt,
Nat. Phys., 5, 304,(2009)" and its supplementary notes.
I didn’t understand this paper well since I am not familiar with basic theory relevant to electric field in solids. Could you give me some advice on this warning in calculations?
Thanks,
Daosheng
Relax structure at finite electric field [SOLVED]
Relax structure at finite electric field [SOLVED]
 Attachments

 run_relax_under_Efield.in
 (1.46 KiB) Downloaded 48 times
Re: Relax structure at finite electric field
Hi Daosheng,
this gets a bit technical. The electrical displacement D you probably know about (berryopt 16), and you'll have to read the supplementary of Stengel's paper to grasp the difference for the reduced stuff. Or just trust them and use berryopt 14 (which is what I would do)!
The issue is with cell relaxation rather than positions, and the change in spatial metric when you modify the cell (so you really have a quantity which is held fixed during the relaxation). The comments in their SI after eq 31 are quite clear, actually, have a look. Do be careful the units will be different from the E field, but if you are looking to do finite differences anyway it should be fine.
this gets a bit technical. The electrical displacement D you probably know about (berryopt 16), and you'll have to read the supplementary of Stengel's paper to grasp the difference for the reduced stuff. Or just trust them and use berryopt 14 (which is what I would do)!
The issue is with cell relaxation rather than positions, and the change in spatial metric when you modify the cell (so you really have a quantity which is held fixed during the relaxation). The comments in their SI after eq 31 are quite clear, actually, have a look. Do be careful the units will be different from the E field, but if you are looking to do finite differences anyway it should be fine.
Matthieu Verstraete
University of Liege, Belgium
University of Liege, Belgium
Re: Relax structure at finite electric field
Dear Professor,
Actually, I want to obtain the optimized structure at finite electric field which can be used for further phonon calculations. The reduced E only supports cell relaxation without atomic position. So, how should we do to relax both cells and atomic positions at finite electric field?
Thanks,
Daosheng
Actually, I want to obtain the optimized structure at finite electric field which can be used for further phonon calculations. The reduced E only supports cell relaxation without atomic position. So, how should we do to relax both cells and atomic positions at finite electric field?
Thanks,
Daosheng
Re: Relax structure at finite electric field
Dear Daosheng,
A good advise is to start to understand the theory related to the calculations you want to do, otherwise you'll not understand the output data you get and you'll not be able to criticize your results if they go wrong, this is quite important.
Regarding the status of the implementation with applied electric field.
 Relaxing the cell and atomic position with "regular" applied electric field works for both normconcerving and PAW. It is only reduced to parallelization on kpoints only.
 Getting the phonons under applied electric field works only in sequential and might have a bug in it (see viewtopic.php?f=10&t=3987&p=12184&hilit=Antonio#p12184 and viewtopic.php?f=10&t=4346&p=12989&hilit ... eld#p12989). I works ONLY with normconserving.
 Applied Dfield also works only in sequential and normconserving (see above links).
 Reduced efield, I'm not sure if it has been extended to PAW but normconserving is OK and probably only in sequential too.
The reduced efield is better than simple applied Efield when you relax the cell parameter because in E.P term P depends on the volume.
In any case, relaxing the internal atomic positions always work in all the cases, why do you say that only cell parameter can be relaxed?
Best wishes,
Eric
A good advise is to start to understand the theory related to the calculations you want to do, otherwise you'll not understand the output data you get and you'll not be able to criticize your results if they go wrong, this is quite important.
Regarding the status of the implementation with applied electric field.
 Relaxing the cell and atomic position with "regular" applied electric field works for both normconcerving and PAW. It is only reduced to parallelization on kpoints only.
 Getting the phonons under applied electric field works only in sequential and might have a bug in it (see viewtopic.php?f=10&t=3987&p=12184&hilit=Antonio#p12184 and viewtopic.php?f=10&t=4346&p=12989&hilit ... eld#p12989). I works ONLY with normconserving.
 Applied Dfield also works only in sequential and normconserving (see above links).
 Reduced efield, I'm not sure if it has been extended to PAW but normconserving is OK and probably only in sequential too.
The reduced efield is better than simple applied Efield when you relax the cell parameter because in E.P term P depends on the volume.
In any case, relaxing the internal atomic positions always work in all the cases, why do you say that only cell parameter can be relaxed?
Best wishes,
Eric
Re: Relax structure at finite electric field
Dear Eric,
Thank you very much for your detailed reply. I got it how to relax structure at finite electric field now.
It is OK to relax both cells and atomic position, which I made a mistake.
Best,
Daosheng
Thank you very much for your detailed reply. I got it how to relax structure at finite electric field now.
It is OK to relax both cells and atomic position, which I made a mistake.
Best,
Daosheng