<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN">
<html><body style='font-family: Verdana,Geneva,sans-serif'>
<p>Dear Iurii, <br /><br />I've ran the calculations using JTH PAW pseudopotentials, the ones I have tested and use typically. <br /><br />I get two issues: with my QE 6.7 compiled with Intel 19 I get a "Index of the second rotated atom= 0   Error in routine symonpair (1):  Out of bounds" error, while with QE 6.8 compiled with Intel 16 I get a weird EOF error (https://gitlab.com/QEF/q-e/-/issues/415). Both builds have natx increased and should accept an 80 atom cell. <br />The in and out files are here. https://drive.google.com/file/d/1XCQPA8j_Z9-R1os5YKMbVPINkeVCOofQ/view?usp=sharing<br /><br />So, now I am waiting for the resolution of the EOF issue, which is some strange compiler dependence of file read or write (I am not sure, but it also shows up in the boron nitride HP tests, while other tests are passed). <br /><br />Regarding the occupations - I've read the THEOS article on that long ago and for some reason went with Gaussian smearing in particular. <br />Smearing is not even needed for insulators, I could have used "fixed". But, as long as I was studying systems with defects (including some where occupied defect levels were at proximity of conduction band) - I wanted to give my systems more convergence flexibility. Gaussian smearing with small degauss gives a result very similar to the one with the fixed occupations - in terms of optimized positions, forces and energies. I might have had some convergence issues with "fixed" though, which I did not have with "Gaussian". <br /><br />However, I've ran a test series, varying degauss and occupations kind. For degauss=0.001 and 0.01, "gauss" and "mv" smearing result in the same total force as "fixed", and the values is the same for k-point grids 2x2x2, 3x3x3, 4x4x4, 5x5x5, 6x6x6 (cell size is 10.4 Angstom though).<br />Wth degauss 0.1 and Gaussian smearing the force was still the same, and went up at degauss 1.0. For MV, the increase happened at degauss=0.1. <br />I found no significant run time difference between the smearing kinds - but that was without defects, pure lutetium oxide. <br />I think I will go with MV for further projects. <br /><br /></p>
<table border="0" cellspacing="0"><colgroup span="7" width="114"></colgroup>
<tbody>
<tr>
<td align="left" height="23"> </td>
<td align="left">grid ></td>
<td align="left">2x2x2</td>
<td align="left">3x3x3</td>
<td align="left">4x4x4</td>
<td align="left">5x5x5</td>
<td align="left">6x6x6</td>
</tr>
<tr>
<td align="left" height="23">occupations</td>
<td align="left">degauss</td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
</tr>
<tr>
<td align="left" height="23">fixed</td>
<td align="left"> </td>
<td align="right">0.046397</td>
<td align="right">0.046403</td>
<td align="right">0.046401</td>
<td align="right">0.046402</td>
<td align="right">0.046404</td>
</tr>
<tr>
<td align="left" height="23"> </td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
</tr>
<tr>
<td align="left" height="23">gauss</td>
<td align="right">0.001</td>
<td align="right">0.046368</td>
<td align="right">0.046379</td>
<td align="right">0.046377</td>
<td align="right">0.046375</td>
<td align="right">0.046389</td>
</tr>
<tr>
<td align="left" height="23">gauss</td>
<td align="right">0.01</td>
<td align="right">0.046368</td>
<td align="right">0.046379</td>
<td align="right">0.046377</td>
<td align="right">0.046375</td>
<td align="right">0.046389</td>
</tr>
<tr>
<td align="left" height="23">gauss</td>
<td align="right">0.1</td>
<td align="right">0.046488</td>
<td align="right">0.046499</td>
<td align="right">0.046496</td>
<td align="right">0.046494</td>
<td align="right">0.046508</td>
</tr>
<tr>
<td align="left" height="23">gauss</td>
<td align="right">1</td>
<td align="right">0.071618</td>
<td align="right">0.07521</td>
<td align="right">0.07425</td>
<td align="right">0.074528</td>
<td align="right">0.074605</td>
</tr>
<tr>
<td align="left" height="23"> </td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
<td align="left"> </td>
</tr>
<tr>
<td align="left" height="23">mv</td>
<td align="right">0.001</td>
<td align="right">0.046368</td>
<td align="right">0.046379</td>
<td align="right">0.046377</td>
<td align="right">0.046375</td>
<td align="right">0.046389</td>
</tr>
<tr>
<td align="left" height="23">mv</td>
<td align="right">0.01</td>
<td align="right">0.046368</td>
<td align="right">0.046379</td>
<td align="right">0.046377</td>
<td align="right">0.046375</td>
<td align="right">0.046389</td>
</tr>
<tr>
<td align="left" height="23">mv</td>
<td align="right">0.1</td>
<td align="right">0.064322</td>
<td align="right">0.064315</td>
<td align="right">0.064313</td>
<td align="right">0.064315</td>
<td align="right">0.064315</td>
</tr>
<tr>
<td align="left" height="23">mv</td>
<td align="right">1</td>
<td align="right">0.064302</td>
<td align="right">0.064879</td>
<td align="right">0.064763</td>
<td align="right">0.06543</td>
<td align="right">0.064942</td>
</tr>
</tbody>
</table>
<p><br />Finally, as DFT+U works fine, I am running tests to see if fixed occupations will result in any significant difference in the calculated U. <br /><br />Thank you. <br />Andrii <br /><br /><br /></p>
<p><br />W dniu 2021-10-15 12:11, Iurii TIMROV napisał:</p>
<blockquote type="cite" style="padding-left:5px; border-left:#1010ff 2px solid; margin-left:5px"><!-- html ignored --> <!-- head ignored --><!-- meta ignored -->
<div id="divtagdefaultwrapper" style="font-size: 12pt; color: #000000; font-family: Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Dear Andrii,</p>
<p> </p>
<p>I received your input files (private communication) and I have investigated them. There are several issues with your inputs.</p>
<p> </p>
<ol style="margin-bottom: 0px; margin-top: 0px;">
<li>Setting Hubbard_V to 0.0 or 1.d-10 is not the same: 0.0 means that you do not consider V while 1.d-10 (or some other very small number) means that you consider V.</li>
<li>You put <span>Hubbard_V(i,j,1) = 1d-10 in a weird way (I mean I do not understand the logic when you choose for which atoms to set Hubbard_V): you just need to put it on those atoms which you want to consider for the Hubbard V interactions (one per atomic type is enough, because the pw.x code will automatically consider all the other atoms of the same type as Hubbard atoms). In your system you have Hf, Lu, and O, and I see that you do not put V on Lu atoms: in this case you must list Lu atoms the last in the ATOMIC_POSITIONS card otherwise the HP code will stop (because you need to list first Hubbard atoms and then all the remaining non-Hubbard atoms, this is explained in the HP documentation).</span></li>
<li><span>You are using norm-conserving <span>FHI98PP</span> pseudopotentials. I recommend to use the SSSP library (PBEsol functional for solids): <a id="LPlnk799088" class="OWAAutoLink" href="https://www.materialscloud.org/discover/sssp/table/efficiency"> https://www.materialscloud.org/discover/sssp/table/efficiency</a></span></li>
<li><span>Since you are using <span>norm-conserving</span> <span>pseudopotentials</span>, you are making a mistake by setting <span>ecutwfc = 40</span> and <span>ecutrho = 400</span> : ecutrho must be exactly 4 times larger than ecutwfc (not more, not less) - in fact, there is a warning in the pw.x output file (which you should inspect). For ultrasoft <span><span>pseudopotentials</span></span> and PAW this factor is typically 8-12 (one has to perform convergence tests).  </span></li>
<li><span><span>smearing = 'gaussian'</span> : Gaussian smearing is not the most efficient one, I suggest to try MV.</span></li>
<li><span><span>degauss = 0.001</span> this is too small. Typical values are in the range 0.01-0.02 Ry but this depends on which k-mesh do you use, check here: <a id="LPlnk278123" class="OWAAutoLink" href="http://theossrv1.epfl.ch/Main/ElectronicTemperature"> http://theossrv1.epfl.ch/Main/ElectronicTemperature</a></span></li>
<li><span>When you are preparing a QE input it is good to start with the QE input generator: <a id="LPlnk274582" class="OWAAutoLink" href="https://www.materialscloud.org/work/tools/qeinputgenerator"> https://www.materialscloud.org/work/tools/qeinputgenerator</a></span></li>
<li><span><span>find_atpert = 2</span> in HP is not the default option: I would start with the default which is find_atpert = 1, it is better because the HP code checks the occupations of all atoms and perturbs non-equivalent atoms that have different occupations. Since in your system you have many O atoms it is likely that there are non-equilavent O atoms (this depends on distortions in your system).</span></li>
<li><span>You need to clarify whether you want only Hubbard V for Hf(5d)-O(2p) or also for Lu(5d)-O(2p). If you want to consider both then you need to set Hubbard_V also for Lu, otherwise put Lu the last in the <span>ATOMIC_POSITIONS</span> card.</span></li>
<li><span>Please try to make these changes and try to run again the code using QE v6.8. You system is too large to run on my workstation, I need to use a cluster, but before that try to make tests on your own and let us know if there still problems (as before, please share your input AND output files via e.g. Google Drive).</span></li>
</ol>
<div> </div>
<div>HTH</div>
<div> </div>
<div>Greetings,</div>
<div>Iurii</div>
<p> </p>
<div id="Signature">
<div id="divtagdefaultwrapper" style="font-size: 12pt; color: #000000; font-family: Calibri, Helvetica, sans-serif, 'EmojiFont', 'Apple Color Emoji', 'Segoe UI Emoji', NotoColorEmoji, 'Segoe UI Symbol', 'Android Emoji', EmojiSymbols;" dir="ltr">
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;">
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;"><span style="color: #808080; font-family: 'Times New Roman', Times, serif; font-size: medium;">--<br /> <span style="font-family: Cambria,Georgia,serif;">Dr. Iurii TIMROV</span><br /> <span style="font-family: Cambria,Georgia,serif;">Senior Research Scientist</span></span></div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;"><span style="color: #808080; font-family: 'Times New Roman', Times, serif; font-size: medium;"><span style="font-family: Cambria,Georgia,serif;">Theory and Simulation of Materials (THEOS)</span></span></div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;"><span style="font-size: 12pt; color: #808080; font-family: Cambria,Georgia,serif;">Swiss Federal Institute of Technology Lausanne (EPFL</span><span style="font-size: 12pt; color: #808080; font-family: 'Times New Roman',Times,serif;"><span style="font-family: 'Times New Roman', Times, serif;"><span style="font-family: Cambria,Georgia,serif;">)</span></span></span></div>
<span style="font-family: Cambria,Georgia,serif;"></span><span style="color: #808080;"></span><span style="font-family: Cambria,Georgia,serif;"></span>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;"><span style="font-family: Cambria,Georgia,serif;"></span><span style="color: #808080; font-family: 'Times New Roman', Times, serif; font-size: medium;"><span style="font-family: Cambria,Georgia,serif;">CH-1015 Lausanne, Switzerland</span><br /> <span style="font-family: Cambria,Georgia,serif;">+41 21 69 34 881</span></span></div>
<span style="font-family: Cambria,Georgia,serif;"></span>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;"><span style="font-family: Cambria,Georgia,serif;"></span><a id="LPNoLP" href="http://people.epfl.ch/265334"><span style="font-family: Cambria,Georgia,serif;">http://people.epfl.ch/265334</span></a></div>
</div>
</div>
</div>
</div>
<hr style="display: inline-block; width: 98%;" />
<div id="divRplyFwdMsg" dir="ltr"><span style="font-size: 11pt; color: #000000; font-family: Calibri, sans-serif;"><strong>From:</strong> users <users-bounces@lists.quantum-espresso.org> on behalf of Iurii TIMROV via users <users@lists.quantum-espresso.org><br /> <strong>Sent:</strong> Thursday, October 14, 2021 10:35:02 AM<br /> <strong>To:</strong> Andrii Shyichuk; Quantum ESPRESSO users Forum<br /> <strong>Subject:</strong> Re: [QE-users] Hubbard U+V for a large system</span>
<div> </div>
</div>
<div>
<div id="divtagdefaultwrapper" style="font-size: 12pt; color: #000000; font-family: Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Dear Andrii,</p>
<p> </p>
<p><span style="font-size: small;"><span style="font-size: 10pt;">> SCF run is fine, file hp.x run gives:<br /> > Index of the second rotated atom=           0<br /> >       Error in routine symonpair (1):<br /> >       Out of bounds</span></span></p>
<p> </p>
<p>It is strange that the pw.x codes works fine while you have some issues with hp.x. These two codes use the same routine symonpair. </p>
<p> </p>
<p>Which version of QE do you use? Can you share please your input and output files (via Google Drive) so that we can investigate your problem in more detail?</p>
<p> </p>
<p>Greetings,</p>
<p>Iurii</p>
<p> </p>
<div id="Signature">
<div id="divtagdefaultwrapper" style="font-size: 12pt; color: #000000; font-family: Calibri, Helvetica, sans-serif, 'EmojiFont', 'Apple Color Emoji', 'Segoe UI Emoji', NotoColorEmoji, 'Segoe UI Symbol', 'Android Emoji', EmojiSymbols;" dir="ltr">
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;">
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;"><span style="color: #808080; font-family: 'Times New Roman', Times, serif; font-size: medium;">--<br /> <span style="font-family: Cambria,Georgia,serif;">Dr. Iurii TIMROV</span><br /> <span style="font-family: Cambria,Georgia,serif;">Senior Research Scientist</span></span></div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;"><span style="color: #808080; font-family: 'Times New Roman', Times, serif; font-size: medium;"><span style="font-family: Cambria,Georgia,serif;">Theory and Simulation of Materials (THEOS)</span></span></div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;"><span style="font-size: 12pt; color: #808080; font-family: Cambria,Georgia,serif;">Swiss Federal Institute of Technology Lausanne (EPFL</span><span style="font-size: 12pt; color: #808080; font-family: 'Times New Roman',Times,serif;"><span style="font-family: 'Times New Roman', Times, serif;"><span style="font-family: Cambria,Georgia,serif;">)</span></span></span></div>
<span style="font-family: Cambria,Georgia,serif;"></span><span style="color: #808080;"></span><span style="font-family: Cambria,Georgia,serif;"></span>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;"><span style="font-family: Cambria,Georgia,serif;"></span><span style="color: #808080; font-family: 'Times New Roman', Times, serif; font-size: medium;"><span style="font-family: Cambria,Georgia,serif;">CH-1015 Lausanne, Switzerland</span><br /> <span style="font-family: Cambria,Georgia,serif;">+41 21 69 34 881</span></span></div>
<span style="font-family: Cambria,Georgia,serif;"></span>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; margin: 0;"><span style="font-family: Cambria,Georgia,serif;"></span><a id="LPNoLP" href="http://people.epfl.ch/265334"><span style="font-family: Cambria,Georgia,serif;">http://people.epfl.ch/265334</span></a></div>
</div>
</div>
</div>
</div>
<hr style="display: inline-block; width: 98%;" />
<div id="divRplyFwdMsg" dir="ltr"><span style="font-size: 11pt; color: #000000; font-family: Calibri, sans-serif;"><strong>From:</strong> users <users-bounces@lists.quantum-espresso.org> on behalf of Johannes Meusburger - STFC UKRI <Johannes.Meusburger@stfc.ac.uk><br /> <strong>Sent:</strong> Wednesday, October 13, 2021 10:25:41 PM<br /> <strong>To:</strong> Andrii Shyichuk; Quantum ESPRESSO users Forum<br /> <strong>Subject:</strong> Re: [QE-users] Hubbard U+V for a large system</span>
<div> </div>
</div>
<div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;">Dear Andriii,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;"> </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;">As far as I know, the default settings limit DFT + U + V to 50 atoms. However, this should be readily adjusted by changing the natx parameter in the Modules/parameters.f90 (<a href="https://gitlab.com/QEF/q-e/-/blob/develop/Modules/parameters.f90">Modules/parameters.f90 · develop · QEF - Quantum Espresso Foundation / q-e · GitLab</a>) to your system size.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;"> </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;">Also, having a look at your input file I do not think that you should <span style="color: #201f1e; font-family: 'Segoe UI', 'Segoe UI Web (West European)', 'Segoe UI', -apple-system, BlinkMacSystemFont, Roboto, 'Helvetica Neue', sans-serif; font-size: 14.6667px; background-color: #ffffff; display: inline !important;"> Hubbard_U(3) = 1d-10  AND Hubbard_V(3,3,1), since Hubbard_V(3,3,1) already refers to the on-site Hubbard U for atom 3.  To sum up, if you are using DFT + U use Hubbard_U, if you are using DFT + U + V just use the Hubbard_V keyword.</span> </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;"> </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;">The solutions to both of these problems are well described in the pw.x  input description (<a href="https://www.quantum-espresso.org/Doc/INPUT_PW.html#idm471">pw.x: input description (quantum-espresso.org)</a>, hence I thought you might find it helpful to give it a read. Also, I highly recommend working your way through the examples on the QE github page (<a href="https://github.com/QEF/q-e/blob/master/HP/examples/example10/reference/LiCoO2.scf.in">q-e/LiCoO2.scf.in at master · QEF/q-e · GitHub</a>).</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;"> </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;">HTH,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;"> </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;">Johannes Meusburger</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: #000000;">ISIS Neutron and Muon Source & Diamond Light Source & University of Exeter, UK</div>
<div id="appendonsend"> </div>
<hr style="display: inline-block; width: 98%;" />
<div id="divRplyFwdMsg" dir="ltr"><span style="font-size: 11pt; color: #000000; font-family: Calibri, sans-serif;"><strong>From:</strong> users <users-bounces@lists.quantum-espresso.org> on behalf of Andrii Shyichuk via users <users@lists.quantum-espresso.org><br /> <strong>Sent:</strong> Wednesday, October 13, 2021 8:52 PM<br /> <strong>To:</strong> users@lists.quantum-espresso.org <users@lists.quantum-espresso.org><br /> <strong>Subject:</strong> [QE-users] Hubbard U+V for a large system</span>
<div> </div>
</div>
<div class="BodyFragment">
<div class="PlainText">Dear Users,<br /> <br /> I'm trying to run a DFT+U+V Hubbard parameters calculation (hp.x) on a <br /> system with 80 atoms.<br /> <br /> SCF run is fine, file hp.x run gives:<br /> <br /> Index of the second rotated atom=           0<br />       Error in routine symonpair (1):<br />       Out of bounds<br /> <br /> How can the second rotated atom be zero?<br /> <br /> My input is based on the YouTube "Tutorial on DFT+U+V using Quantum <br /> ESPRESSO (v6.7)".<br /> The +U+V settings are:<br />                     lda_plus_u = .true.<br />                lda_plus_u_kind = 2<br />                   Hubbard_U(3) = 1d-10<br />               Hubbard_V(1,1,1) = 1d-10<br />               Hubbard_V(2,2,1) = 1d-10<br />               Hubbard_V(3,3,1) = 1d-10<br />               Hubbard_V(4,4,1) = 1d-10<br />               ...<br />               Hubbard_V(80,80,1) = 1d-10<br /> <br /> Is it a bug, or is it just the question of some parameter tweaking?<br /> I've already increased natx to 150.<br /> <br /> Thank you.<br /> Andrii Shyichuk, University of Wrocław.<br /> _______________________________________________<br /> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu">www.max-centre.eu</a>)<br /> users mailing list users@lists.quantum-espresso.org<br /> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div>
</div>
<p><span style="font-size: 6pt;">This email and any attachments are intended solely for the use of the named recipients. If you are not the intended recipient you must not use, disclose, copy or distribute this email or any of its attachments and should notify the sender immediately and delete this email from your system. UK Research and Innovation (UKRI) has taken every reasonable precaution to minimise risk of this email or any attachments containing viruses or malware but the recipient should carry out its own virus and malware checks before opening the attachments. UKRI does not accept any liability for any losses or damages which the recipient may sustain due to presence of any viruses. </span></p>
</div>
</div>
</blockquote>
<p> </p>
<div> </div>
</body></html>