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