<div>Dear all,</div>
<div> </div>
<div>I was doing an electronic minimization using the CPV code cp.x. To make sure I did</div>
<div>not make mistakes, I also did a similar calculation using the PWscf code pw.x.</div>
<div> </div>
<div>Well, when I compared the results, I found quite some difference, in particular, </div>
<div>the forces acting on atoms are quite different while the total energy are close.</div>
<div> </div>
<div>My question is whether I should expect to get same results from the two codes?</div>
<div>If so, there must be something wrong with the input files that I used (see below).</div>
<div> </div>
<div>Thanks in advance for any suggestions that you may have!<br><br>-Peter</div>
<div> </div>
<div>--------------------------pw.x input file---------<br>C2H4<br> &control<br>    calculation='scf'<br>    restart_mode='from_scratch',<br>    wf_collect=.true.,<br>    prefix='pwscf',<br>    pseudo_dir = '/espresso-
3.0/pseudo',<br>    outdir='.',<br>    tprnfor=.true.<br> /<br> &system<br>                       ibrav = 1<br>                   celldm(1) = 37.7945197733<br>    nat=6, ntyp=2,<br>    ecutwfc = 30.0,<br>    occupations='fixed',
<br> /<br> &electrons<br>    conv_thr =  1.0d-8<br>    mixing_beta = 0.7<br> /<br>ATOMIC_SPECIES<br>C   12.0107 C.pbe-van_ak.UPF<br>H   1.00794 H.pbe-van_ak.UPF<br>K_POINTS {gamma}<br>ATOMIC_POSITIONS {crystal}<br>C 0.46664509
 0.5 0.5<br>C 0.53335491 0.5 0.5<br>H 0.43793698 0.5 0.45351762<br>H 0.56206302 0.5 0.45351762<br>H 0.56206302 0.5 0.54648238<br>H 0.43793698 0.5 0.54648238<br>
<div> </div>
<div>--------------------------cp.x input file---------<br> &control<br>    calculation='cp',<br>    restart_mode='from_scratch',<br>    ndr=52, ndw=53,<br>    nstep=400, iprint=10, isave=50,<br>    dt=
5.0,<br>    etot_conv_thr  = 1.d-9,<br>    ekin_conv_thr  = 1.d-5,<br>    prefix = 'cp'<br>    pseudo_dir='/espresso-3.0/pseudo/',<br>    tprnfor = .TRUE.<br>    outdir='.',<br> /<br> &system<br>
                       ibrav = 1<br>                   celldm(1) = 37.7945197733<br>    nat=6, ntyp=2,<br>    ecutwfc = 30.0,<br> /<br> &electrons<br>    electron_dynamics='damp'<br> /<br> &ions<br>    ion_dynamics='none',
<br> /<br>ATOMIC_SPECIES<br>C   12.0107 C.pbe-van_ak.UPF<br>H   1.00794 H.pbe-van_ak.UPF<br>ATOMIC_POSITIONS {crystal}<br>C 0.46664509 0.5 0.5<br>C 0.53335491 0.5 0.5<br>H 0.43793698 0.5 0.45351762<br>H 0.56206302 0.5 0.45351762
<br>H 0.56206302 0.5 0.54648238<br>H 0.43793698 0.5 0.54648238</div></div>