<div dir='auto'><div>Dear Omar A. Ashour,</div><div dir="auto"><br></div><div dir="auto">Thank you so much for your detailed guidance to tackle the issue. I will take a step back and re-implement the input parameters as per your suggestions and try to run the calculation again.</div><div dir="auto"><br></div><div dir="auto">Regards,</div><div dir="auto">Bhargab</div><div><br></div><div data-smartmail="gmail_signature"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 8 Sept 2024 1:20 am, "Omar A. Ashour via users" <users@lists.quantum-espresso.org> wrote:<br type="attribution" /><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div dir="ltr">
Be careful translating VASP input files directly into QE because the implementation details can differ, so just converting the units isn’t enough. </div><div dir="ltr"><br /></div><div dir="ltr">electron_maxtep = 3000 —> if you can’t converge in a couple hundred steps, you should really take a step back and assess what your system is doing. Your calculation did 628 steps, and the SCF accuracy was all over the place--the last step is ~2 Ry. Something is clearly wrong.<br /></div><div dir="ltr"><br /></div><div dir="ltr">Increasing your cutoff to at least 60 Ry will possibly fix your convergence problems. VASP’s PAW sets are incredibly soft; you can’t just convert their cutoff to Ry and expect it to work with a different pseudopotential. Read the UPF files for your pseudos: Ga has a suggested cutoff of 60 Ry, and that’s the minimum you should use.</div><div dir="ltr"><br /></div><div dir="ltr">Other notes:</div><div dir="ltr"><ol start="1" dir="ltr" style="margin:0px"><li>Your pseudos are wrong. QE implements SOC differently than VASP. You need to use fully relativistic pseudos.</li><li>If this is metallic, use cold smearing with a larger degauss; your k-grid is somewhat coarse. It will probably help convergence, too.</li><li>conv_thr is very large; there is no good reason to go above the default 1e-6 for a cell of this size, and you may need to go lower depending on what you’re trying to do. Note that QE’s implementation is different from VASP’s, see equation A7 in this paper, and I believe it’s missing a factor of V^2/2 with V the cell volume. </li><li>for_conv_thr is far too large; again, don’t go above the default unless you have a good reason. If I recall correctly, QE applies the force convergence criteria to each component, whereas VASP applies it to the magnitude of the force, which is something to be aware of.</li><li>It won’t make a difference here for this system, but if you use ibrav=0 you should have more decimal places in your CELL_PARAMETERS card (10 or so) or the symmetry analysis may not work properly. </li></ol><div><br /></div><div dir="ltr">You’re relaxing a noncolinear magnet, so due to the complicated energy landscape, you should expect more difficult convergence, even with a proper input file.</div><div dir="ltr"><br /></div><div dir="ltr">Best,</div><div dir="ltr"><br /></div><div dir="ltr">Omar A. Ashour</div><div dir="ltr">UC Berkeley </div><div dir="ltr"><br /></div></div>
<div class="elided-text">
<div dir="ltr">On Sep 7, 2024 at 1:53:16 AM, Bhargab Kakati <<a href="mailto:bhrgabkakati@iasst.res.in">bhrgabkakati@iasst.res.in</a>> wrote:<br /></div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb( 204 , 204 , 204 );padding-left:1ex">
<div>
<div>
On 2024-09-05 23:58, Bhargab Kakati wrote:<br /><blockquote> Thank you so for your suggestion, I'll try and get back yo you if I<br /></blockquote><blockquote> encounter some issues.<br /></blockquote><blockquote> <br /></blockquote><blockquote> Regards,<br /></blockquote><blockquote> Bhargab<br /></blockquote><blockquote> <br /></blockquote><blockquote> On 5 Sept 2024 5:13 pm, <a href="mailto:sogenyi@student.eaifr.org">sogenyi@student.eaifr.org</a> wrote:<br /></blockquote><blockquote> <br /></blockquote><blockquote>> Dear Kakati,<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> I think you need to apply same condition for ecut and ecutrho when<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> using ultrasoft potential to PAW pseudopotential i.e ecutrho = 8 or<br /></blockquote><blockquote>> 10<br /></blockquote><blockquote>> times ecut.<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> Sunday.<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> Quoting Bhargab Kakati <<a href="mailto:bhrgabkakati@iasst.res.in">bhrgabkakati@iasst.res.in</a>>:<br /></blockquote><blockquote>> <br /></blockquote><blockquote>>> Dear experts,<br /></blockquote><blockquote>>> <br /></blockquote><blockquote>>> I am trying to relax my alloy system but it is not converging even<br /></blockquote><blockquote>> <br /></blockquote><blockquote>>> after running for a very long time. I hava attached here the input<br /></blockquote><blockquote>> <br /></blockquote><blockquote>>> and outputs of the run. Any help in resloving the issue would be<br /></blockquote><blockquote>> <br /></blockquote><blockquote>>> highly appreciated.<br /></blockquote><blockquote>>> <br /></blockquote><blockquote>>> Regards,<br /></blockquote><blockquote>>> Bhargab Kakati<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> _______________________________________________<br /></blockquote><blockquote>> The Quantum ESPRESSO community stands by the Ukrainian<br /></blockquote><blockquote>> people and expresses its concerns about the devastating<br /></blockquote><blockquote>> effects that the Russian military offensive has on their<br /></blockquote><blockquote>> country and on the free and peaceful scientific, cultural,<br /></blockquote><blockquote>> and economic cooperation amongst peoples<br /></blockquote><blockquote>> _______________________________________________<br /></blockquote><blockquote>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu">www.max-centre.eu</a>)<br /></blockquote><blockquote>> users mailing list <a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br /></blockquote><blockquote>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br /></blockquote><blockquote> Thank you so for your suggestion, I'll try and get back yo you if I <br /></blockquote><blockquote> encounter some issues.<br /></blockquote><blockquote> <br /></blockquote><blockquote> Regards,<br /></blockquote><blockquote> Bhargab<br /></blockquote><blockquote> <br /></blockquote><blockquote> <br /></blockquote><blockquote> On 5 Sept 2024 5:13 pm, <a href="mailto:sogenyi@student.eaifr.org">sogenyi@student.eaifr.org</a> wrote:<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> Dear Kakati,<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> I think you need to apply same condition for ecut and ecutrho when <br /></blockquote><blockquote>> using ultrasoft potential to PAW pseudopotential i.e ecutrho = 8 or <br /></blockquote><blockquote>> 10 <br /></blockquote><blockquote>> times ecut.<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> Sunday.<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> <br /></blockquote><blockquote>> Quoting Bhargab Kakati <<a href="mailto:bhrgabkakati@iasst.res.in">bhrgabkakati@iasst.res.in</a>>:<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> > Dear experts,<br /></blockquote><blockquote>> ><br /></blockquote><blockquote>> > I am trying to relax my alloy system but it is not converging even <br /></blockquote><blockquote>> > after running for a very long time. I hava attached here the input <br /></blockquote><blockquote>> > and outputs of the run. Any help in resloving the issue would be <br /></blockquote><blockquote>> > highly appreciated.<br /></blockquote><blockquote>> ><br /></blockquote><blockquote>> > Regards,<br /></blockquote><blockquote>> > Bhargab Kakati<br /></blockquote><blockquote>> <br /></blockquote><blockquote>> <br /></blockquote><blockquote>> _______________________________________________<br /></blockquote><blockquote>> The Quantum ESPRESSO community stands by the Ukrainian<br /></blockquote><blockquote>> people and expresses its concerns about the devastating<br /></blockquote><blockquote>> effects that the Russian military offensive has on their<br /></blockquote><blockquote>> country and on the free and peaceful scientific, cultural,<br /></blockquote><blockquote>> and economic cooperation amongst peoples<br /></blockquote><blockquote>> _______________________________________________<br /></blockquote><blockquote>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu">www.max-centre.eu</a>)<br /></blockquote><blockquote>> users mailing list <a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br /></blockquote><blockquote>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br /></blockquote><br />Hello,<br /><br />I tried to relax my system bu increasing ecutrho = 8 times of ecut and <br />also changed the pseudopotentials too. It is still not converging. I <br />have attached the new input and output files here.<br /><br />Regards,<br />Bhargab
</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="http://www.max-centre.eu">www.max-centre.eu</a>)<br />users mailing list <a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br /><a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a>
</div>
</div>
</blockquote>
</div></div>
</blockquote></div><br></div>