Dear PWSCF users,<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I am now running a calculation of
strained STO&nbsp; 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&nbsp; 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 &quot;probably
because G_par is NOT a reciprocal lattice vector&quot;.<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 &#39;lelfield = .true.&#39;.<br>
Thank you very much for your help.<br>
<br>
Hanghui Chen<br>
Department of Physics <br>
Yale University<br>