Dear PWSCF users,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I have done a calculation of LaAlO3 recently, using &#39;vc-relax&#39;
with the contraints that only z-direction is relaxed and x,y directions
are fixed, and in the presence of finite electric field (based on Berry
Phase method). The input file is attached at the end.<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; However I am very confused with the
results. The equilibrium
lattice constant in z direction is c=0.95*a (a is the lattice constant
in x and y direction).&nbsp; The final ionic polarization is 72.03 a.u.
and
electronic polarization is 4.50 a.u. Then I use the equilibrium unit
cell (c=0.95*a) and just &#39;relax&#39; the atoms in the presence of the same
electric field. I obtain: ionic polarization is 75.13 a.u. and
electronic polarization is 1.36 a.u. Even more confusing is that if I
use the equilibrium unit cell (c=0.95*a) and&nbsp; put the atom in the
ideal position (i.e. no ionic polarization), the reference ionic
polarization is 73.24 a.u. <br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Shall I trust the &#39;vc-relax&#39; calculation in the presence of finite electric field?<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Thank you very much.<br>
<br>
Input file reads:<br>
&amp;CONTROL<br>
&nbsp;calculation=&#39;vc-relax&#39;<br>
&nbsp;wf_collect=.true.<br>
&nbsp;pseudo_dir = &#39;./psp&#39;<br>
&nbsp;outdir=&#39;/scratch&#39;<br>
&nbsp;wfcdir=&#39;/scratch&#39;<br>
&nbsp;prefix=&#39;sLAO-vc&#39;<br>
&nbsp;tprnfor = .true.<br>
&nbsp;tstress = .true.<br>
&nbsp;disk_io=&#39;low&#39;<br>
&nbsp;verbosity=&#39;high&#39;<br>
&nbsp;dt=80.D0<br>
&nbsp;lelfield=.true.<br>
&nbsp;gdir=3<br>
&nbsp;nppstr=16<br>
&nbsp;nberrycyc=1<br>
/<br>
&amp;SYSTEM<br>
&nbsp;ibrav=&nbsp; 6<br>
&nbsp;celldm(1) = 7.27<br>
&nbsp;celldm(3) = 0.9506<br>
&nbsp;nat=&nbsp; 5<br>
&nbsp;ntyp= 3<br>
&nbsp;ecutwfc = 30.0<br>
&nbsp;ecutrho = 180.0<br>
&nbsp;occupations=&#39;fixed&#39;<br>
&nbsp;nosym=.true.<br>
/&amp;ELECTRONS<br>
&nbsp;diagonalization=&#39;david&#39;<br>
&nbsp;mixing_beta = 0.7D0<br>
&nbsp;diago_david_ndim = 4<br>
&nbsp;startingwfc=&#39;random&#39;<br>
&nbsp;startingpot=&#39;atomic&#39;<br>
&nbsp;efield=0.00348<br>
&nbsp;conv_thr = 1.0d-9<br>
/<br>
&amp;IONS<br>
&nbsp;ion_dynamics = &#39;damp&#39;<br>
&nbsp;phase_space = &#39;full&#39;<br>
&nbsp;pot_extrapolation = &#39;second_order&#39;<br>
&nbsp;wfc_extrapolation = &#39;second_order&#39;<br>
/<br>
&amp;CELL<br>
&nbsp;cell_dynamics = &#39;damp-w&#39;<br>
&nbsp;cell_dofree = &#39;z&#39;<br>
/<br>
ATOMIC_SPECIES<br>
&nbsp;La&nbsp; 138.91 057-La-ca-nsp-hanghui.uspp.format.UPF<br>
&nbsp;Al&nbsp; 26.98&nbsp; 013-Al-ca--hanghui.uspp.format.UPF<br>
&nbsp;O&nbsp;&nbsp; 16.00&nbsp; 008-O-ca--vgrp.uspp.format.UPF<br>
ATOMIC_POSITIONS crystal<br>
&nbsp;La 0.00 0.00 0.00<br>
&nbsp;Al 0.50 0.50 0.50<br>
&nbsp;O&nbsp; 0.50 0.50 0.00<br>
&nbsp;O&nbsp; 0.00 0.50 0.50<br>
&nbsp;O&nbsp; 0.50 0.00 0.50<br>
K_POINTS crystal<br>
256<br>
(The k points are omitted here)<br>
<br>
Hanghui Chen<br>
Department of Physics <br>
Yale University<br>
<br>
<br>