<div>Dear PWSCF users,&nbsp;&nbsp;&nbsp;</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; In my previous email, I reported that if I calculated a symmetric<br>&nbsp;&nbsp; system (both the material and the external sawtooth potential are<br>&nbsp;&nbsp; symmetric), I found that the forces are not symmetric. I think I may find<br>
&nbsp;&nbsp; the reason.<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The total force in the presence of sawtooth potential has two parts:<br>&nbsp;&nbsp; one is from the variation of Kohn-Sham functional; the other is from<br>&nbsp;&nbsp; sawtooth potential which in the code is &quot;forcefield&quot; (see the routine<br>
&nbsp;&nbsp; forces.f90). The additional force is calculated by the routine<br>&nbsp;&nbsp; add_efield.f90. However in routine add_efield.f90, it only calculates the<br>&nbsp;&nbsp; magnitude of the force but does not specify the direction (the additional<br>
&nbsp;&nbsp; force only has one direction).&nbsp; Therefore all the forces in the &#39;reversed<br>&nbsp;&nbsp; region&#39; are not correctly calculated because the direction of &#39;forcefield&#39;<br>&nbsp;&nbsp; is NOT reversed as is the external electric field.<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I sincerely hope that the people who wrote the routine add_efield.f90<br>&nbsp;&nbsp; can reexamine the code and add the directions to the &quot;forcefield&quot;.<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; No offense intended.<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Thank you very much.<br>
&nbsp;&nbsp; Hanghui Chen<br>&nbsp;&nbsp; Department of Physics, Yale University</div>