<div>Dear PWSCF users,</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 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.&nbsp;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&nbsp;(turning on&nbsp;sawtooth potential) are given below. </div>

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