<div dir="ltr">Dear Saiyed,<div>it is better if next time you also include output, just to understand what is not converging.</div><div><br></div><div>In this case:</div><div>i) for convergence issues of the electronic self-consistent cycle, you should always try to reduce mixing_beta to 0.3, 0.1, 0.05 ...</div><div><br></div><div>ii) before running ANY calculations, mostly if you're dealing with a material that is novel for you, it is quite a good practice to visualize your structure.</div><div>I'm not expert about this material, however what I can notice is that, if you're running the same calculation as that you described in your first email,</div><div>that is, 3 layers of calcite, I cannot see three layers but just one and the one I see I suspect does not resemble that much to the models in the paper you mention.</div><div>You could figure out that you're likely not simulating three layers, by looking at your coordinates: the z coordinates range from 0.46234 to 0.53766, so they span an</div><div>interval 0.07532 wide. Because these are crystal units, this in Angstrom is 0.07532 *  43.2578.... Ang ~ 3.258 Ang. It is quite unlikely that three layers of any material</div><div>(even graphene. the flattest you can think about) are just 3 Ang thick.</div><div><br></div><div>Giovanni<br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(0,0,0)">-- <br><br>Giovanni Cantele, PhD<br>CNR-SPIN<br>c/o Dipartimento di Fisica<br>Universita' di Napoli "Federico II"<br>Complesso Universitario M. S. Angelo - Ed. 6<br>Via Cintia, I-80126, Napoli, Italy<br><a href="mailto:giovanni.cantele@spin.cnr.it" style="color:rgb(17,85,204)" target="_blank">e-mail: giovanni.cantele@spin.cnr.it</a><br>Phone: +39 081 676910<br>Skype contact: giocan74<br><br>ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" style="color:rgb(17,85,204)" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>Web page: </span><a href="https://sites.google.com/view/giovanni-cantele/home" style="color:rgb(17,85,204)" target="_blank">https://sites.google.com/view/giovanni-cantele/home</a><br></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno sab 16 mar 2024 alle ore 22:45 Saiyed Tasnim Md Fahim <<a href="mailto:sfahim@usc.edu">sfahim@usc.edu</a>> ha scritto:<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 dir="ltr"><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">Dear Dr. Cantele, </div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">Thank you for your response. </div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">I tried the steps that you mentioned. However, it still didn't converge.  I am following this paper as the benchmark : <a href="https://doi.org/10.1016/j.cej.2024.148940" rel="noreferrer noopener" aria-label="Persistent link using digital object identifier" title="Persistent link using digital object identifier" style="background-color:transparent;box-sizing:border-box;margin:0px;padding:0px;word-break:break-word;color:rgb(31,31,31);text-decoration-line:none;display:inline-block;outline-width:0px;font-family:ElsevierSans,Arial,Helvetica,Roboto,"Lucida Sans Unicode","Microsoft Sans Serif","Segoe UI Symbol",STIXGeneral,"Cambria Math","Arial Unicode MS",sans-serif;font-size:14px" target="_blank"><span style="box-sizing:border-box;margin:0px;padding:0px;border-bottom:2px solid rgb(235,101,0)">https://doi.org/10.1016/j.cej.2024.148940</span></a> </div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">I have used the norm conserving pseudopotentials from Pseudodojo website as USPP/PAW wasn't implemented for TS dispersion. I used the ASE (atomic simulation environment) to generate the atomic positions of the 2x2 supercell.</div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">I have attached the new input file with this mail. Could you please tell me if any parameters I have missed or need to be changed or removed? Thank you. <br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div style="color:rgb(34,34,34)"><font face="trebuchet ms, sans-serif">Sincerely,</font></div><div style="color:rgb(34,34,34)"><font face="trebuchet ms, sans-serif">Saiyed Tasnim Md Fahim,</font></div><div style="color:rgb(34,34,34)"><font face="trebuchet ms, sans-serif">Ph.D. Student,</font></div><div style="color:rgb(34,34,34)"><font face="trebuchet ms, sans-serif">Environmental Engineering,</font></div><div style="color:rgb(34,34,34)"><font face="trebuchet ms, sans-serif">University of Southern California.</font></div></div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Feb 28, 2024 at 12:27 AM Giovanni Cantele <<a href="mailto:giovanni.cantele@spin.cnr.it" target="_blank">giovanni.cantele@spin.cnr.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 dir="ltr">Dear Saiyed,<div><br></div><div>when convergence issues show up there are several tricks you might consider, one of these is certainly reducing mixing_beta,</div><div>another could be increasing degauss.</div><div><br></div><div>However, before trying any trick to cure a calculation particularly difficult to converge, as a rule of thumb the first step is to</div><div>check your input file parameters including cut-off, k-point grid, degauss, GEOMETRY, pseudopotentials, etc.</div><div><br></div><div>After a quick look to your input file I can notice:</div><div>- an incredibly huge value of degauss (1d-1 Ry = 0.1 Ry = 1.36 eV, this would smear even the band gap of silicon I think)</div><div>- a mixing beta of 0.6, actually the advise of reducing it in cases of difficult convergence is meant to go down to even 0.1 or 0.05</div><div>- but MOST IMPORTANTLY you specify ATOMIC_POSITIONS in crystal units, that are usually expected to be numbers in the range</div><div>[-1,1] because in units of the lattice vectors there should be no reason to specify the position of an atom in a unit cell far apart from the</div><div>origin. On the other hand, the numbers in your input file seem to indicate that you're actually giving the positions in Angstrom, because they are</div><div>huge. If you try to visualize your input file using, for example XCrysDen, you can easily see atoms at incredibly huge distances. The geometry</div><div>becomes much more reasonables if you replace "crystal" with angstrom. After doing that, recheck the geometry (with XCrysDen) to understand if it is</div><div>what you expect or if there are still issues, before running any calculation.</div><div><br></div><div>Giovanni</div><div><br></div><div><br></div><div><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><span style="color:rgb(0,0,0)">-- <br><br>Giovanni Cantele, PhD<br>CNR-SPIN<br>c/o Dipartimento di Fisica<br>Universita' di Napoli "Federico II"<br>Complesso Universitario M. S. Angelo - Ed. 6<br>Via Cintia, I-80126, Napoli, Italy<br><a href="mailto:giovanni.cantele@spin.cnr.it" style="color:rgb(17,85,204)" target="_blank">e-mail: giovanni.cantele@spin.cnr.it</a><br>Phone: +39 081 676910<br>Skype contact: giocan74<br><br>ResearcherID: <a href="https://urldefense.com/v3/__http://www.researcherid.com/rid/A-1951-2009__;!!LIr3w8kk_Xxm!uHeu46wuGTYOEVDT0_V9s92ksXMy-0mbaezxExpDSYHQ7RiP0fPUrWo0ZpZ_1BDUcfbWaoy76a1xhBwsJ_ARu6Qgfu-r$" style="color:rgb(17,85,204)" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>Web page: </span><a href="https://urldefense.com/v3/__https://sites.google.com/view/giovanni-cantele/home__;!!LIr3w8kk_Xxm!uHeu46wuGTYOEVDT0_V9s92ksXMy-0mbaezxExpDSYHQ7RiP0fPUrWo0ZpZ_1BDUcfbWaoy76a1xhBwsJ_ARuypOJcP_$" style="color:rgb(17,85,204)" target="_blank">https://sites.google.com/view/giovanni-cantele/home</a><br></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mer 28 feb 2024 alle ore 08:41 Saiyed Tasnim Md Fahim <<a href="mailto:sfahim@usc.edu" target="_blank">sfahim@usc.edu</a>> ha scritto:<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 dir="ltr"><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><div class="gmail_default">Dear users,</div><div class="gmail_default">I am working with calcite in Quantum Espresso 7.1. I want to perform a vc-relax calculation for 3 layers of calcite but after 100 iterations convergence is not achieved. I changed the mixing beta but it did not work.</div><div class="gmail_default"><br></div><div class="gmail_default">I have attached the input file with this mail. Kindly let me know how I should proceed for both the convergence.</div><div class="gmail_default"><br></div></div><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div style="color:rgb(34,34,34)"><font face="trebuchet ms, sans-serif">Sincerely,</font></div><div style="color:rgb(34,34,34)"><font face="trebuchet ms, sans-serif">Saiyed Tasnim Md Fahim,</font></div><div style="color:rgb(34,34,34)"><font face="trebuchet ms, sans-serif">Ph.D. Student,</font></div><div style="color:rgb(34,34,34)"><font face="trebuchet ms, sans-serif">Environmental Engineering,</font></div><div style="color:rgb(34,34,34)"><font face="trebuchet ms, sans-serif">University of Southern California.</font></div></div></div></div></div>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="https://urldefense.com/v3/__http://www.max-centre.eu__;!!LIr3w8kk_Xxm!uHeu46wuGTYOEVDT0_V9s92ksXMy-0mbaezxExpDSYHQ7RiP0fPUrWo0ZpZ_1BDUcfbWaoy76a1xhBwsJ_ARu28Aj1dA$" rel="noreferrer" target="_blank">www.max-centre.eu</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://urldefense.com/v3/__https://lists.quantum-espresso.org/mailman/listinfo/users__;!!LIr3w8kk_Xxm!uHeu46wuGTYOEVDT0_V9s92ksXMy-0mbaezxExpDSYHQ7RiP0fPUrWo0ZpZ_1BDUcfbWaoy76a1xhBwsJ_ARuyERRL04$" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="https://urldefense.com/v3/__http://www.max-centre.eu__;!!LIr3w8kk_Xxm!uHeu46wuGTYOEVDT0_V9s92ksXMy-0mbaezxExpDSYHQ7RiP0fPUrWo0ZpZ_1BDUcfbWaoy76a1xhBwsJ_ARu28Aj1dA$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__http://www.max-centre.eu__;!!LIr3w8kk_Xxm!uHeu46wuGTYOEVDT0_V9s92ksXMy-0mbaezxExpDSYHQ7RiP0fPUrWo0ZpZ_1BDUcfbWaoy76a1xhBwsJ_ARu28Aj1dA$</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://urldefense.com/v3/__https://lists.quantum-espresso.org/mailman/listinfo/users__;!!LIr3w8kk_Xxm!uHeu46wuGTYOEVDT0_V9s92ksXMy-0mbaezxExpDSYHQ7RiP0fPUrWo0ZpZ_1BDUcfbWaoy76a1xhBwsJ_ARuyERRL04$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__https://lists.quantum-espresso.org/mailman/listinfo/users__;!!LIr3w8kk_Xxm!uHeu46wuGTYOEVDT0_V9s92ksXMy-0mbaezxExpDSYHQ7RiP0fPUrWo0ZpZ_1BDUcfbWaoy76a1xhBwsJ_ARuyERRL04$</a> </blockquote></div>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</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>