<div dir="ltr"><font face="comic sans ms,sans-serif"><font size="4">The other way I tried computing isolated energy of Ca and O, is using the following separate inputs for Ca and O. Are the input scripts valid? The energies got converged for the following scripts.</font></font><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jun 27, 2020 at 4:45 PM Pietro Delugas <<a href="mailto:pdelugas@sissa.it">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 class="gmail-m_-3191345924098248726WordSection1"><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><span style="font-size:12pt">Dear Pooja <u></u><u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><span style="font-size:12pt">I might not have understooo what you want to do. <u></u><u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><span style="font-size:12pt">Looking at  your input it looks like you are not actually computing the energy of an isolated system, but instead the energy of <u></u><u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><span style="font-size:12pt">rocksalt CaO. In that case it is not isolated system, remove assume_isolated from the input, use sensible values for alat and the computation should converge. <u></u><u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><span style="font-size:12pt">Pietro <u></u><u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><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">Mail</a> for Windows 10</p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><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">Pooja Vyas</a><br><b>Sent: </b>Saturday, June 27, 2020 11:56 AM<br><b>To: </b><a href="mailto:users@lists.quantum-espresso.org" target="_blank">Quantum ESPRESSO users Forum</a><br><b>Subject: </b>Re: [QE-users] Energy of isolated atom of CaO</p></div><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><div><div><p class="MsoNormal">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:</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">PS: I want to calculate isolated energy of CaO using example05/PW (i.e Martyna-Tuckerman) method.</p></div><div><p class="MsoNormal">#!/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</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Any kind of help is appreciated.</p></div><div><p class="MsoNormal">Thanks and Regards.</p></div></div><p class="MsoNormal"><u></u> <u></u></p><div><div><p class="MsoNormal">On Fri, Jun 26, 2020 at 11:58 AM Pietro Delugas <<a href="mailto:pdelugas@sissa.it" target="_blank">pdelugas@sissa.it</a>> wrote:</p></div><blockquote style="border-color:currentcolor currentcolor currentcolor rgb(204,204,204);border-style:none none none solid;border-width:medium medium medium 1pt;padding:0in 0in 0in 6pt;margin-left:4.8pt;margin-right:0in"><div><div><p class="MsoNormal"><span class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><span style="font-size:12pt">Hello </span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><span style="font-size:12pt"> </span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><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. </span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><span style="font-size:12pt"> </span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><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. </span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><span style="font-size:12pt"> </span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><span style="font-size:12pt">Regards - Pietro  </span></span></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><span style="font-size:12pt"> </span></span></p><p class="MsoNormal">Sent from <a href="https://go.microsoft.com/fwlink/?LinkId=550986" target="_blank">Mail</a> for Windows 10</p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><span style="font-size:12pt"> </span></span></p><div style="border-style:solid none none;border-width:1pt medium medium;padding:3pt 0in 0in;border-color:currentcolor"><p class="MsoNormal"><b>From: </b><a href="mailto:poojavyas595@gmail.com" target="_blank">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">Quantum ESPRESSO users Forum</a><br><b>Subject: </b>[QE-users] Energy of isolated atom of CaO</p></div><p class="MsoNormal"><span class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><span style="font-size:12pt"> </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"> </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"> </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 class="gmail-m_-3191345924098248726gmail-m-704748173674001904defaultfonthxmailstyle"><span style="font-size:12pt"> </span></span></p></div></div></blockquote></div><p class="MsoNormal" style="margin-left:4.8pt">_______________________________________________<br>Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></p><p class="MsoNormal"><span class="gmail-m_-3191345924098248726DefaultFontHxMailStyle"><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" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>