[Pw_forum] lelfield and lberry
chenhanghuipwscf at gmail.com
Tue Aug 26 00:23:52 CEST 2008
Dear PWSCF users,
I am now running a calculation of strained STO and try to calculate
its spontaneous polarization based on Berry phase method. I try two
1) run a self-consistent calculation turning on lelfield=.true. and set
2) first run a self-consistent calculation and then run a
non-self-consistent calculation turning on lberry=.true.
Are these two method equivalent? When I try these two methods on the same
system, method 1) works well but method 2) fails because "probably because
G_par is NOT a reciprocal lattice vector".
I notice that (in the version 4.0) in the two subroutines bp_c_phase.f90 and
c_phase_field.f90 the definition of gpar are different by a factor
DBLE(nppstr)/DBLE(nppstr-1). But I think these two gpar should be the same
I just want to ask whether the option lberry is fully tested and compatible
with 'lelfield = .true.'.
Thank you very much for your help.
Department of Physics
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the users