<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body>
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Dear <span>Quang</span>,</p>
<p><br>
</p>
<p><span>> starting_magnetization(1) = 0.05, starting_magnetization(2) = -0.05,</span><br>
</p>
<p><br>
</p>
<p>You have only one type in your input, while here you specify the starting magnetization for two types, which is wrong. Moreover, bulk Si is nonmagnetic so you should not specify starting_magnetization. But if you do, then you need to uncomment the lines
 below</p>
<p><span>!   occupations = ’smearing’, smearing = ’mv’, degauss = 0.01,</span></p>
<p><br>
</p>
<p><span>> Since qe-7.1, the flag lda_plus_u=.true. is not needed. Input files are below.</span><br>
</p>
<p><br>
</p>
<p>You are right</p>
<p><br>
</p>
<p><span>> The HP code can be used only when lda_plus_u=.true.</span><br>
</p>
<p><br>
</p>
<p>This will be changed, the message is confusing, we need to update it. <br>
</p>
<p><br>
</p>
<p>In the documentation Doc/Hubbard_input.pdf on page 15 it is written:</p>
<p>"<i>Hubbard parameters U and V can be computed using the hp.x code of Quantum ESPRESSO.<br>
</i></p>
<div><i>However, the hp.x currently supports the calculations of U and V for one Hubbard channel per<br>
atomic type. In other words, the advanced features presented above (i.e. cross-manifold inter-<br>
actions) are currently not implemented in hp.x.</i>"<br>
</div>
<p></p>
<p><br>
</p>
<p>So what you want to do cannot be currently done with HP.</p>
<p><br>
</p>
<p>HTH</p>
<p><br>
</p>
<p>Iurii</p>
<p><br>
</p>
<div id="Signature">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, "EmojiFont", "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font size="3" face="'Times New Roman', Times, serif" color="808080">--<br>
<span style="font-family:Cambria,Georgia,serif">Dr. Iurii TIMROV</span><br>
<span style="font-family:Cambria,Georgia,serif">Senior Research Scientist</span></font></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font size="3" face="'Times New Roman', Times, serif" color="808080"><span style="font-family:Cambria,Georgia,serif">Theory and Simulation of Materials (THEOS)</span></font></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<span style="font-size:12pt; color:rgb(128,128,128); font-family:Cambria,Georgia,serif">Swiss Federal Institute of Technology Lausanne (EPFL</span><font style="font-size:12pt; color:rgb(128,128,128); font-family:"Times New Roman",Times,serif" color="808080"><font face="'Times New Roman', Times, serif"><span style="font-family:Cambria,Georgia,serif">)</span></font></font><br>
</div>
<span style="font-family:Cambria,Georgia,serif"></span><font color="808080"></font><span style="font-family:Cambria,Georgia,serif"></span>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<span style="font-family:Cambria,Georgia,serif"></span><font size="3" face="'Times New Roman', Times, serif" color="808080"><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></font></div>
<span style="font-family:Cambria,Georgia,serif"></span>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<span style="font-family:Cambria,Georgia,serif"></span><a href="http://people.epfl.ch/265334" tabindex="0" id="LPNoLP"><span style="font-family:Cambria,Georgia,serif">http://people.epfl.ch/265334</span></a><br>
</div>
</div>
</div>
</div>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> users <users-bounces@lists.quantum-espresso.org> on behalf of VNU Quang Tran <tranquang@vnu.edu.vn><br>
<b>Sent:</b> Friday, July 1, 2022 10:48:56 AM<br>
<b>To:</b> users@lists.quantum-espresso.org<br>
<b>Subject:</b> [QE-users] DFT+U+V: The HP code can be used only when lda_plus_u=.true. in qe-7.1</font>
<div> </div>
</div>
<div>
<div dir="ltr">Dear All,
<div>I just want to check DFT+U+V for a simple case (before going beyond), Si with qe-7.1.</div>
<div>PWSCF runs well. However, the self-consistent U V calculation with hp.x stops with the following message</div>
<div>
<div>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
     Error in routine hp_readin (1):<br>
     The HP code can be used only when lda_plus_u=.true.<br>
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
<br>
     stopping ...<br>
</div>
<div>Since qe-7.1, the flag lda_plus_u=.true. is not needed. Input files are below.</div>
</div>
<div>Can you tell me how to overcome this, please?</div>
<div>Many thanks,</div>
<div>Quang</div>
<div><font face="monospace">PhD, </font><span style="font-family:monospace">VNU University of Engineering and Technology, Hanoi, Vietnam</span></div>
<div><span style="font-family:monospace"><br>
</span></div>
<div>---PW input--------<br>
</div>
<div> &control<br>
</div>
<div>   prefix       = 'Si',<br>
   calculation  = 'scf',<br>
   verbosity    = 'high',<br>
   pseudo_dir = '~/pseudo/PBE'<br>
   outdir       = './tmp',<br>
 /<br>
 &system<br>
   ibrav = 2, celldm(1) = 10.261212854,<br>
   nat = 2,<br>
   ntyp = 1,<br>
   ecutwfc = 60,<br>
!   occupations = ’smearing’, smearing = ’mv’, degauss = 0.01,<br>
   starting_magnetization(1) = 0.05, starting_magnetization(2) = -0.05,<br>
 /<br>
 &electrons<br>
   conv_thr        =  1.0D-10,<br>
   mixing_beta     =  0.7,<br>
 /<br>
 ATOMIC_SPECIES<br>
   Si 28.086  Si.pbe-rrkj.UPF<br>
 ATOMIC_POSITIONS crystal<br>
   Si   0.000   0.000   0.000<br>
   Si   0.250   0.250   0.250<br>
 K_POINTS automatic<br>
   6  6  6  0  0  0<br>
HUBBARD {ortho-atomic}<br>
V Si-3p Si-3p 1 1 0.01<br>
V Si-3p Si-3s 1 1 0.01<br>
V Si-3s Si-3s 1 1 0.01<br>
V Si-3s Si-3p 1 1 0.01<br>
V Si-3p Si-3p 1 2 0.01<br>
V Si-3p Si-3s 1 2 0.01<br>
V Si-3s Si-3s 1 2 0.01<br>
V Si-3s Si-3p 1 2 0.01<br>
</div>
<div>--------</div>
<div>HP input</div>
<div>&inputhp<br>
   prefix = 'Si'<br>
   outdir = './tmp'<br>
   lmin = 1<br>
   nq1 = 2, nq2 = 2, nq3 = 2<br>
   conv_thr_chi = 1.d-6<br>
   iverbosity = 2<br>
/</div>
</div>
</div>
</body>
</html>