<div dir="ltr">Deal All,<div><br></div><div>I want to relax a system of water molecules with sodium ion, I have tried to run it before and it did not converge. Please tell me what i need to improve on my input file, this is what I have done for my input file:</div><div><br></div><div>&control</div><div>   calculation = 'relax',</div><div>   restart_mode ='from_scratch',</div><div>   pseudo_dir = '/nfs/share/potentials/psl/pbe',</div><div>   outdir = './H2ONa/',</div><div>   tprnfor = .true.,</div><div><br></div><div>/</div><div>&system</div><div>   ibrav= 0,</div><div>   nat = 91,</div><div>   ntyp = 3,</div><div>   ecutwfc = 90.0,</div><div>   ecutrho = 400.0,</div><div>   tot_charge = +1,</div><div>     occupations = 'smearing',</div><div>     smearing = 'mp',</div><div>     degauss = 0.02,</div><div>     nosym = .TRUE.                                           </div><div>/</div><div>&electrons</div><div>mixing_beta = 0.015,</div><div><br></div><div>/</div><div>&ions</div><div>ion_dynamics = 'bfgs',</div><div>  </div><div>/</div><div><br></div><div>ATOMIC_SPECIES</div><div>O   15.9994    O.pbe-n-kjpaw_psl.1.0.0.UPF</div><div>H   1.00794    H.pbe-kjpaw_psl.1.0.0.UPF</div><div>Na  22.9898   Na.pbe-spn-kjpaw_psl.1.0.0.UPF  </div></div>