<div>Dear PWSCF users,</div>
<div>      I calculate a symmetric STO slab (separated by vacuum) and turn on a symmetric sawtooth potential (the turning point is right at the center of the STO slab). I fix all the atoms. Without the sawtooth potential, the forces are symmetric and there are no forces in the central layer. However, when I turn on the sawtooth potential which I set to be symmetric, the forces on the atoms are NO longer symmetric and the forces in the layer are NO longer zero. The input file and the forces (turning on sawtooth potential) are given below. </div>

<div>     Why are the forces NOT symmetric but the configuration (including the sawtooth potential) is symmetric? Does the turning point in the material cause any anomaly?</div>
<div>     Thank you very much.</div>
<div> </div>
<div>Hanghui Chen</div>
<div>Department of Physics, Yale University</div>
<div> </div>
<div>Input file:</div>
<div>&CONTROL<br>  calculation='scf'<br>  wf_collect=.true.<br>  pseudo_dir = './psp'<br>  outdir='.'<br>  wfcdir='/scratch'<br>  prefix='test1'<br>  disk_io='low'<br>  verbosity='high'  <br>
  tprnfor = .true.<br>  tstress = .true.<br>  dt=80.D0<br>  tefield=.true.<br>/<br>&SYSTEM   <br>  ibrav=  6<br>  celldm(1) = 7.27<br>  celldm(3) = 6.0<br>  nat=  12<br>  ntyp= 3<br>  ecutwfc = 30.0<br>  ecutrho = 180.0<br>
  occupations='smearing'<br>  smearing='gauss'<br>  degauss=0.005<br>  edir=3<br>  emaxpos=0.5<br>  eopreg=0.5<br>  eamp=0.01<br>/<br>&ELECTRONS<br>  diagonalization='david'<br>  mixing_beta = 0.7D0<br>
  diago_david_ndim = 4<br>  mixing_mode= 'local-TF'<br>  electron_maxstep = 300<br>  startingpot = 'file'<br>  conv_thr= 1.d-8<br>/<br>&IONS<br>  ion_dynamics = 'bfgs'<br>  phase_space = 'full'<br>
  pot_extrapolation = 'first_order'<br>  wfc_extrapolation = 'first_order'<br>/<br>&CELL<br>  cell_dynamics = 'damp-w'<br>/<br>ATOMIC_SPECIES<br> Sr  87.62  038-Sr-ca-sp-vgrp.uspp.format.UPF<br>
 Ti  47.90  022-Ti-ca-sp-vgrp.uspp.format.UPF<br> O   16.00  008-O-ca--vgrp.uspp.format.UPF<br>ATOMIC_POSITIONS alat<br> Sr 0.00 0.00 2.00 <br> O  0.50 0.50 2.00 <br> Ti 0.50 0.50 2.50 <br> O  0.50 0.00 2.50 <br> O  0.00 0.50 2.50 <br>
 Sr 0.00 0.00 3.00 <br> O  0.50 0.50 3.00 <br> Ti 0.50 0.50 3.50 <br> O  0.50 0.00 3.50<br> O  0.00 0.50 3.50<br> Sr 0.00 0.00 4.00<br> O  0.50 0.50 4.00<br>K_POINTS {automatic}<br>8 8 1 0 0 0<br></div>
<div>Forces:</div>
<div> atom   1 type  1   force =     0.00000001   -0.00000001   -0.09670746<br> atom   2 type  3   force =     0.00000000    0.00000000   -0.21698818<br> atom   3 type  2   force =     0.00000000    0.00000001   -0.17490285<br>
 atom   4 type  3   force =     0.00000000    0.00000000   -0.17341486<br> atom   5 type  3   force =     0.00000001    0.00000000   -0.17341491<br> atom   6 type  1   force =    -0.00000002    0.00000000    0.06540405<br>
 atom   7 type  3   force =    -0.00000001   -0.00000001   -0.01253063<br> atom   8 type  2   force =     0.00000001    0.00000001    0.32344891<br> atom   9 type  3   force =     0.00000000    0.00000000    0.08311003<br>
 atom  10 type  3   force =     0.00000000    0.00000000    0.08311005<br> atom  11 type  1   force =     0.00000000    0.00000000    0.16600421<br> atom  12 type  3   force =     0.00000001    0.00000000    0.12688164<br>
</div>
<div> </div>