Dear PWSCF users,<br>
      I am now running a calculation of
strained STO  and try to calculate its spontaneous polarization
based on Berry phase method. I try two different ways: <br>
1) run a self-consistent calculation turning on lelfield=.true. and set efield=0.0 <br>
2) first run a self-consistent calculation and  then run a non-self-consistent calculation turning on lberry=.true.<br>
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".<br>
<br>
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 thing?<br>
<br>
I just want to ask whether the option lberry is fully tested and compatible with 'lelfield = .true.'.<br>
Thank you very much for your help.<br>
<br>
Hanghui Chen<br>
Department of Physics <br>
Yale University<br>