Dear all,<br><br>I comment out the line to impose total force = 0 in forces.f90, and perform a calculation of one atom in sawlike potential. The input file is attached below. I had expected the force exerted on nucleus by the field  would cancel the force on ion due to change in electron charge density, so the total force should vanish since the system is neutral. I also varied the cutoff energy, but failed to see a convergence behavior as shown below. Could anyone give me a clue? Thank you!<br>
<br>Ecut(Ry)   Total force(Ry/au)<br>30            0.002729<br>40            0.003052<br>50            0.003848<br>60            0.000086<br>70            0.000962<br>80            0.000137<br>90            0.003208<br>100          0.001271<br>
<br>The input file:<br>-------------------------------------------------------------------------<br> &control<br>    calculation = 'scf'<br>    restart_mode='from_scratch',<br>    prefix='c',<br>
    tstress = .true.<br>    tprnfor = .true.<br>    pseudo_dir = '.',<br>    outdir='/scratch/local'<br>    tefield=.true.<br> /<br> &system<br>    ibrav=  0, nat=  1, ntyp= 1,<br>ecutwfc=100<br>    edir=3<br>
    emaxpos=0.95<br>    eopreg=0.1<br>    eamp=0.001<br> /<br> &electrons<br>    diagonalization='cg'<br>    mixing_mode = 'plain'<br>    mixing_beta = 0.7<br>    conv_thr =  1.0d-8<br> /<br>CELL_PARAMETERS<br>
10 0 0<br>0 10 0<br>0 0 20<br>ATOMIC_SPECIES<br> C 1 C.pbe-rrkjus.UPF<br>ATOMIC_POSITIONS {bohr}<br> C 0.00 0.00 9.0<br>K_POINTS gamma<br>----------------------------------------------------------------------------<br><br>
Best regards<br><br>zah<br clear="all"><br>-- <br>------------------------------------------------------------------<br>Research Fellow,<br>Dept. of Physics, National Univ. of Singapore<br>10 Kent Ridge Crescent, Singapore, 119260<br>
Tel: +65 6516 7844<br>------------------------------------------------------------------<br><br>