<div dir="ltr"><font color="#000000">Hi all,</font><div><font color="#000000"><br></font></div><div><font color="#000000">I'm interested in doing DFT+U calculations on the rare earth mono-nitride GdN. The most successful calculations in the past have used multiple Hubbard U terms (U_f and U_d) on the 4f and 5d electron channels on the Gd ion. This has been done before using QE (<font face="arial, sans-serif"><span style="white-space:nowrap">Comput. Mater. Sci 95 (2014) 263-270</span>)</font>. </font></div><div><font color="#000000"><br></font></div><div><font color="#000000">I suspect that in v6.6 this may now be possible in the stable branch via an onsite Hubbard V term on both the standard and background orbitals?</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">The pw.x input guide gives that Hubbard_V(na,na,k) will be the onsite Hubbard interaction. After setting the appropriate values for the background n, l and occupations of Gd in:</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">Modules/set_hubbard_l.f90</font></div><div><font color="#000000">Modules/set_hubbard_n.f90</font></div><div><font color="#000000">PW/src/tabd.f90</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">and increasing sc_size = 3 in Modules/parameters.f90</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">and compiling the code, I've used the following in the &SYSTEM block:</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">lda_plus_u=.true.</font></div><div><font color="#000000">lda_plus_u_kind=2</font></div><div><font color="#000000">Hubbard_V(1,1,1)=a</font></div><div><font color="#000000">Hubbard_V(1,1,3)=b<br></font></div><div><font color="#000000"><br></font></div><div><font color="#000000">I guess there will now be two Hubbard terms on that atom at side na, on the standard and background orbitals.</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">However, I find the following error:</font></div><div><font color="#000000"><br></font></div><div>Error in routine init_lda_plus_u (1):<br>     lda_plus_u calculation but Hubbard_l not set</div><div><font color="#000000"><br></font></div><div><font color="#000000">Any help is appreciated, </font><span style="color:rgb(0,0,0)">I'm running QE 6.6 on Ubuntu 20.04.</span></div><div><font color="#000000"><br></font></div><div><font color="#000000">Thanks</font></div><div><br></div><div><font color="#000000">Will Hewett</font></div><div><font color="#000000">Postdoctoral Researcher</font></div><div><font color="#000000">Victoria University of Wellington</font></div><div><font color="#000000">New Zealand</font></div></div>