<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear all,<div><br></div><div>I am giving the linear response method for Hubbard U and Hund J another go. A recent publication has obtained quite resonable values for U and J based on a hexahydrated transition metal model and I thought it might be worthwhile trying (<a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.98.235157">https://journals.aps.org/prb/abstract/10.1103/PhysRevB.98.235157</a>) ; however, either I am too stupid or there is a bug in the Hubbard_beta routine that I cannot figure out. Anyways, if someone else has managed to do this please do comment!</div><div><br></div><div>The procedure should be easy, I make a series of calculation for values of Hubbard_beta keeping Hubbard_U and Hubbard_alpha=1D-40; the on-site magnetization is define by m=n(up)-n(down) which I can get from the scf output <br clear="all"><div><br></div><div>bare response:</div><div><div> --- enter write_ns ---</div><div> LDA+U parameters:</div><div>U( 3)     =  0.00000000</div><div>J0( 3)     =  0.00000000</div><div>alpha( 3) =  0.00000000</div><div>beta( 3) =  0.05000000</div><div>atom    1   Tr[ns(na)] (up, down, total) =   <b>4.95101</b>  <b>1.38740</b>  6.33841</div></div><div><br></div><div>converged response</div><div>atom    1   Tr[ns(na)] (up, down, total) =   <b>4.95032</b>  <b>1.39090</b>  6.34122<br></div><div><br></div><div>If I do this for a series of calculations I obtain the following:</div><div><br></div><div><div># beta bare conv</div><div>-0.08 3.57155 3.57023</div><div>-0.05 3.56860 3.56776</div><div>-0.02 3.56561 3.56528</div><div>0.02  3.56160 3.56194</div><div>0.05  3.55856 3.55942</div><div>0.08  3.55548 3.55686</div><div># slopes: bare= -0.100419 conv= -0.0835161 <br></div></div><div><br></div><div>the respective \partial m/\partial \beta look linear in the region of the chosen Hubbard_beta.</div><div><br></div><div>but the resulting J=-2*(bare^-1-conv^-1) makes no sense and even the sign is wrong.</div><div><br></div><div>The value for Hubbard_U obtained by this method using \partial n/\partial \alpha is around 5.4, which seems perfectly suitable by the way.</div><div><br></div><div>Input for a single calculation is attached below.</div><div><br></div><div>Thanks for your help in advance!</div><div><br>Chris </div><div><br></div><div><br></div><div><div>&CONTROL</div><div>calculation = "scf",<br></div><div>...</div><div>/</div><div>&SYSTEM</div><div>ibrav=1,A=15.000000,ecutwfc=30, ecutrho=400,</div><div>starting_magnetization(3)=0.7,</div><div> nspin=2,occupations="fixed", tot_magnetization=4.000000<br></div><div>nat=19,ntyp=3,tot_charge=2.000000</div><div>lda_plus_U=.true.,U_projection_type="ortho-atomic"</div><div>Hubbard_U(3)=1D-40,Hubbard_alpha(3)=1D-40, Hubbard_beta(3)=-0.08</div><div>/</div><div>&ELECTRONS</div><div>startingwfc="file", startingpot="file"</div><div>mixing_beta=0.6,conv_thr=1D-12,diago_thr_init=1D-14,electron_maxstep=150, mixing_mode="local-TF"</div><div>/</div><div>ATOMIC_SPECIES</div><div>H 1.008 H.pbe-kjpaw_psl.1.0.0.UPF</div><div>O 15.999 O.pbe-n-kjpaw_psl.1.0.0.UPF</div><div>Fe1 55.845 Fe.pbe-spn-kjpaw_psl.0.2.1.UPF</div><div>ATOMIC_POSITIONS (angstrom)</div><div>Fe1 7.500000000 7.500000000 7.500000000</div><div>O 7.500000000 7.500000000 9.491812807</div><div>O 7.500000000 9.491814995 7.500000000</div><div>O 9.491776410 7.500000000 7.500000000</div><div>O 5.508226144 7.500000000 7.500000000</div><div>O 7.500000000 5.508190791 7.500000000</div><div>O 7.500000000 7.500000000 5.508184515</div><div>H 6.674174667 10.066164672 7.500000000</div><div>H 8.326002124 10.065917296 7.500000000</div><div>H 10.065724550 7.500000000 8.326098474</div><div>H 10.066137163 7.500000000 6.674195836</div><div>H 7.500000000 6.674283010 4.933686902</div><div>H 7.500000000 8.326129049 4.934261135</div><div>H 6.674236823 4.933772022 7.500000000</div><div>H 8.326063477 4.934192685 7.500000000</div><div>H 4.934274010 7.500000000 6.673909845</div><div>H 4.933876539 7.500000000 8.325806755</div><div>H 7.500000000 6.674482948 10.066592251</div><div>H 7.500000000 8.326364525 10.065407565</div><div>K_POINTS (automatic)</div><div>1 1 1 0 0 0 </div></div><div><br></div><div><br></div><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr">Postdoctoral Researcher<br>Center for Quantum Nanoscience, Institute for Basic Science<br>Ewha Womans University, Seoul, South Korea<blockquote type="cite" style="font-size:12.8px"><div dir="ltr"><div><div dir="ltr"></div></div></div></blockquote></div></div></div></div></div></div></div></div>