Dear PWSCF users,<br>

      I have done a calculation of LaAlO3 recently, using 'vc-relax'
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>
      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).  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 'relax' 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  put the atom in the
ideal position (i.e. no ionic polarization), the reference ionic
polarization is 73.24 a.u. <br>
      Shall I trust the 'vc-relax' calculation in the presence of finite electric field?<br>
      Thank you very much.<br>
<br>
Input file reads:<br>
&CONTROL<br>
 calculation='vc-relax'<br>
 wf_collect=.true.<br>
 pseudo_dir = './psp'<br>
 outdir='/scratch'<br>
 wfcdir='/scratch'<br>
 prefix='sLAO-vc'<br>
 tprnfor = .true.<br>
 tstress = .true.<br>
 disk_io='low'<br>
 verbosity='high'<br>
 dt=80.D0<br>
 lelfield=.true.<br>
 gdir=3<br>
 nppstr=16<br>
 nberrycyc=1<br>
/<br>
&SYSTEM<br>
 ibrav=  6<br>
 celldm(1) = 7.27<br>
 celldm(3) = 0.9506<br>
 nat=  5<br>
 ntyp= 3<br>
 ecutwfc = 30.0<br>
 ecutrho = 180.0<br>
 occupations='fixed'<br>
 nosym=.true.<br>
/&ELECTRONS<br>
 diagonalization='david'<br>
 mixing_beta = 0.7D0<br>
 diago_david_ndim = 4<br>
 startingwfc='random'<br>
 startingpot='atomic'<br>
 efield=0.00348<br>
 conv_thr = 1.0d-9<br>
/<br>
&IONS<br>
 ion_dynamics = 'damp'<br>
 phase_space = 'full'<br>
 pot_extrapolation = 'second_order'<br>
 wfc_extrapolation = 'second_order'<br>
/<br>
&CELL<br>
 cell_dynamics = 'damp-w'<br>
 cell_dofree = 'z'<br>
/<br>
ATOMIC_SPECIES<br>
 La  138.91 057-La-ca-nsp-hanghui.uspp.format.UPF<br>
 Al  26.98  013-Al-ca--hanghui.uspp.format.UPF<br>
 O   16.00  008-O-ca--vgrp.uspp.format.UPF<br>
ATOMIC_POSITIONS crystal<br>
 La 0.00 0.00 0.00<br>
 Al 0.50 0.50 0.50<br>
 O  0.50 0.50 0.00<br>
 O  0.00 0.50 0.50<br>
 O  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>