<div dir="ltr"><div>Dear users,</div><div>I want to calculate cohesive energy of CaO. I'm using cluster_example/PW which uses assume_isolated = 'martyna_tuckerman' . Following is my input script:</div><div><br></div><div>#!/bin/sh<br>for a in 30<br>do<br>cat > ${a}.in << EOF<br>&CONTROL<br>  calculation  = 'relax'<br>  prefix       = "${a}",<br>  pseudo_dir   = "/home/user/pv/cohesive/pseudo/",<br>  outdir       = "/home/user/pv/cohesive/",<br>/<br>&SYSTEM<br>  ibrav     = 1,<br>  celldm(1) = $a,<br>  nat       = 2,<br>  ntyp      = 2,<br>  ecutwfc   = 100.D0,<br>  assume_isolated = 'martyna-tuckerman'<br>/<br>&ELECTRONS<br>  conv_thr    = 1.D-8,<br>  mixing_beta = 0.7D0,<br>/<br>&IONS<br>/<br>ATOMIC_SPECIES<br>Ca  40.078  Ca.pbe-nsp-van.UPF<br>O  15.999  O.pbe-van_ak.UPF<br>ATOMIC_POSITIONS <br>Ca 0.0 0.0 0.0<br>O  0.5 0.5 0.5<br>K_POINTS Gamma<br><br>EOF<br><br>mpirun -np 40 -machinefile x1 /apps/codes/qe/6.4/bin/pw.x < ${a}.in > ${a}.out<br>done<br>done</div><div><br></div><div>My lattice constant is around 8 a.u , so I assumed the box size to be 3 times of it. But when I run the script, the convergence was not achieved even after 100 iterations. Then I tried to take small box size of about 16 bohr. That too didn't work. Then I even tried reducing ecutwfc and mixing_beta..but same error was obtained. Can anyone suggest what can I do to resolve the issue?</div><div>Any kind of help is appreciated.</div><div>Thanks and regards.<br></div></div>