<div dir="auto">Dear Pooja <div dir="auto">I am not expert but I am suggesting to add this line <span style="background-color:rgb(255,255,153);font-family:arial,helvetica,sans-serif;font-size:medium;font-weight:700">electron_maxstep = 500</span></div><div dir="auto"><font face="arial, helvetica, sans-serif" size="3"><span style="background-color:rgb(255,255,153)"><b>Under </b></span></font><span style="color:rgb(33,68,120);font-family:arial,helvetica,sans-serif;font-size:18px">&</span><span style="color:rgb(33,68,120);font-family:arial,helvetica,sans-serif;font-size:18px">ELECTRONS</span></div><div dir="auto"><span style="color:rgb(33,68,120);font-family:arial,helvetica,sans-serif;font-size:18px"><br></span></div><div dir="auto"><span style="color:rgb(33,68,120);font-family:arial,helvetica,sans-serif;font-size:18px">Try and see if it works </span></div><div dir="auto"><span style="color:rgb(33,68,120);font-family:arial,helvetica,sans-serif;font-size:18px">500 is just from my mind you can put any number.  Because the default is 100. </span></div><div dir="auto"><span style="color:rgb(33,68,120);font-family:arial,helvetica,sans-serif;font-size:18px"><br></span></div><div dir="auto"><span style="color:rgb(33,68,120);font-family:arial,helvetica,sans-serif;font-size:18px">Thanks </span></div><div dir="auto"><font color="#214478" face="arial, helvetica, sans-serif"><span style="font-size:18px">Malathe </span></font></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, 27 Jun 2020, 1:56 pm Pooja Vyas, <<a href="mailto:poojavyas595@gmail.com">poojavyas595@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>I tried the way you suggested, tried replacing "relax" by "scf", however , I still face the same issue: No convergence even after 100 iterations. Following is the input file:<br></div><div><br></div><div>PS: I want to calculate isolated energy of CaO using example05/PW (i.e Martyna-Tuckerman) method.<br></div><div>#!/bin/sh<br><br>for a in 30<br>do<br><br>cat > ${a}.in << EOF<br>&CONTROL<br>  calculation  = 'relax'<br>  prefix       = "${a}",<br>  pseudo_dir   = "/home/userpooja/pv/cohesive/pseudo/",<br>  outdir       = "/home/userpooja/pv/cohesive/",<br>/<br>&SYSTEM<br>  ibrav     = 1,<br>  celldm(1) = $a,<br>  nat       = 2,<br>  ntyp      = 2,<br>  ecutwfc   = 100.D0,<br>  nbnd      = 8<br>  assume_isolated = 'martyna-tuckerman'<br>/<br>&ELECTRONS<br>  conv_thr    = 1.D-8,<br>  mixing_beta = 0.7D0,<br>  startingwfc = 'random'<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</div><div><br></div><div>Any kind of help is appreciated.</div><div>Thanks and Regards.<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jun 26, 2020 at 11:58 AM Pietro Delugas <<a href="mailto:pdelugas@sissa.it" target="_blank" rel="noreferrer">pdelugas@sissa.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div lang="EN-US"><div><p class="MsoNormal"><span><span style="font-size:12pt">Hello <u></u><u></u></span></span></p><p class="MsoNormal"><span><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal"><span><span style="font-size:12pt">In ATOMIC_POSITIONS  you need to specify the units that you are  using for the coordinates, if you don’t put anything he program assumes alat units and puts the 2 atoms very far apart. <u></u><u></u></span></span></p><p class="MsoNormal"><span><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal"><span><span style="font-size:12pt">If you actually wanted  compute the two separated atoms specifying startingwfc=’random’ in the &electrons name list may help.  In this case I would also chose calculation=’scf’ in &control. <u></u><u></u></span></span></p><p class="MsoNormal"><span><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal"><span><span style="font-size:12pt">Regards - Pietro  <u></u><u></u></span></span></p><p class="MsoNormal"><span><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal">Sent from <a href="https://go.microsoft.com/fwlink/?LinkId=550986" target="_blank" rel="noreferrer">Mail</a> for Windows 10</p><p class="MsoNormal"><span><span style="font-size:12pt"><u></u> <u></u></span></span></p><div style="border-color:rgb(225,225,225) currentcolor currentcolor;border-style:solid none none;border-width:1pt medium medium;padding:3pt 0in 0in"><p class="MsoNormal" style="border:medium none;padding:0in"><b>From: </b><a href="mailto:poojavyas595@gmail.com" target="_blank" rel="noreferrer">Pooja Vyas</a><br><b>Sent: </b>Friday, June 26, 2020 8:12 AM<br><b>To: </b><a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">Quantum ESPRESSO users Forum</a><br><b>Subject: </b>[QE-users] Energy of isolated atom of CaO</p></div><p class="MsoNormal"><span><span style="font-size:12pt"><u></u> <u></u></span></span></p><div><div><p class="MsoNormal">Dear users,</p></div><div><p class="MsoNormal">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:</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">#!/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</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">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?</p></div><div><p class="MsoNormal">Any kind of help is appreciated.</p></div></div><p class="MsoNormal">Thanks and regards.</p><p class="MsoNormal"><span><span style="font-size:12pt"><u></u> <u></u></span></span></p></div></div>_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>