<div dir="ltr">Hi All,<div><br></div><div>I want to relax a system of 90 water molecules and 1 sodium atom, the calculation does not converge. </div><div><br></div><div>My input file looks like this: </div><div><br></div><div><br></div><div><div>&control</div><div>   calculation = 'vc-relax',</div><div>   restart_mode ='from_scratch',</div><div>   pseudo_dir = '/nfs/share/potentials/psl/pbe',</div><div>   outdir = './2work/',</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 = 70.0,</div><div>   ecutrho = 280.0,</div><div>     occupations = 'smearing',</div><div>     smearing = 'mp',</div><div>     degauss = 0.05</div><div>   </div><div>   </div><div>/</div><div>&electrons</div><div>   mixing_beta = 0.5,</div><div>/</div><div>&ions</div><div>/</div><div>&cell</div><div>cell_dofree = '2Dxy' </div><div>   </div><div>/</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> </div></div></div>