<div dir="ltr">I first relax system and then do vc-relax calculation but show this in the force and strange for me?<div><br></div><div>Because of the before do vc-relax calculation I didn't show this behavior in this force?</div><div><br></div><div>Total force of system converged to   0.000335 but last scf  total force is different and larger than this force (0.004711). in the most calculation this two amount must be same.</div><div><br></div><div><br></div><div>all of total force in vc-relax calculation:</div><div><br></div><div><div>     Total force =     0.004513     Total SCF correction =     0.000180</div><div>     Total force =     0.018979     Total SCF correction =     0.000142</div><div>     Total force =     0.017385     Total SCF correction =     0.000143</div><div>     Total force =     0.039273     Total SCF correction =     0.000107</div><div>     Total force =     0.013076     Total SCF correction =     0.000222</div><div>     Total force =     0.006658     Total SCF correction =     0.000104</div><div>     Total force =     0.005537     Total SCF correction =     0.000050</div><div>     Total force =     0.005065     Total SCF correction =     0.000046</div><div>     Total force =     0.002880     Total SCF correction =     0.000044</div><div>     Total force =     0.000900     Total SCF correction =     0.000005</div><div>     Total force =     0.000570     Total SCF correction =     0.000011</div><div><span style="background-color:rgb(182,215,168)">     Total force =     0.000335     Total SCF correction =     0.000015</span></div><div><span style="background-color:rgb(224,102,102)">     Total force =     0.004711     Total SCF correction =     0.000023</span></div></div><div><br></div><div>my input:</div><div><br></div><div><div>&CONTROL</div><div>  calculation  = "vc-relax",</div><div>  prefix       = "221-50GrporFe",</div><div>  pseudo_dir = ./pseudo',</div><div>  outdir='./tmp3',</div><div>  disk_io='none' ,</div><div>  tprnfor=.TRUE. ,</div><div>  tstress=.TRUE. ,</div><div><span style="background-color:rgb(234,153,153)">  forc_conv_thr=1.0D-4 ,</span></div><div>  etot_conv_thr=1.0D-5,</div><div>/</div><div>&SYSTEM</div><div>  ibrav     = 8,</div><div>  celldm(1) =23.255385645,</div><div>  celldm(2) =1.0392304835091129925263003159134,</div><div>  celldm(3) =1.3001554254551417531019473809222,</div><div>  nat       = 59,</div><div>  ntyp      = 3,</div><div>  ecutwfc   = 50.0,</div><div>  ecutrho   = 300.0,</div><div>  occupations='smearing', </div><div>  smearing='gaussian', </div><div>  degauss=0.01,</div><div>  nspin =2 ,</div><div>  starting_magnetization(1)=  0.0 ,</div><div>  starting_magnetization(2)=  0.5 ,</div><div>  starting_magnetization(3)=  0.0 ,</div><div>  vdw_corr='DFT-D'</div><div>/</div><div>&ELECTRONS</div><div>  conv_thr    = 1.D-8,</div><div>  mixing_beta = 0.7D0,</div><div>/</div><div>&IONS</div><div>  ion_dynamics = 'bfgs'</div><div>/</div><div>&CELL</div><div>cell_factor=1.2 ,</div><div>cell_dofree='xy'</div><div>/</div><div>ATOMIC_SPECIES</div><div>C  12.0107     C.pbe-n-kjpaw_psl.0.1.UPF</div><div>Fe 55.8450     Fe.pbe-spn-kjpaw_psl.0.2.1.UPF</div><div>N  14.0067     N.pbe-n-kjpaw_psl.0.1.UPF</div><div>ATOMIC_POSITIONS (angstrom)</div><div><br></div><div>K_POINTS  {automatic}</div><div>2 2 1 0 0 0</div></div><div><br></div></div>