<div dir="ltr">Dear all,<div><br></div><div>I am running into some problem when using the same species twice in a calculation (eg. O1, O2 for two inequivalent oxygen atoms) and then applying "standard U" on both species, e.g.</div><div><br></div><div>ATOMIC_POSITIONS (angstrom)<br>V             0.0 0.0 0.0  <br>O1            0.0 0.0 2.0   <br>O2            0.0 0.0 -2.0  <br></div><div><br></div><div><br></div><div>lda_plus_u = .true.<br>lda_plus_u_kind = 2<br>U_projection_type = 'ortho-atomic',<br clear="all"><div>Hubbard_V(1,1,1)=1<br>Hubbard_V(2,2,1)=2<br><br></div><div>will work, but adding </div><div><br></div><div>Hubbard_V(3,3,1)=3<br></div><div><br></div><div>will throw the "offset error"; </div><div><br></div><div> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>     Error in routine offset_atom_wfc (2):<br>     wrong offset: your pseudopotential file for atomic species  3 likely does not contain the needed atomic wavefunctions<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br></div><div><br></div><div>I guess that is because the offset is already set from the previous call to offset_atom which is checked at the end of that routine:</div><div><br></div><div>      IF ( ( is_hubbard(nt) .OR. is_hubbard_back(nt) ) .AND. offset(na) < 0 ) THEN<br>        CALL errore('offset_atom_wfc', 'wrong offset: your pseudopotential &<br>                &file for atomic species ' // s // new_line('a') // 'likely &<br>                &does not contain the needed atomic wavefunctions', l_back)<br>      ENDIF <br></div><div><br></div><div>Is there a way to safely detect if the wave-functions are present in the PP without running into this problem? </div><div><br></div><div>Best,</div><div>Chris</div>--<div dir="ltr" data-smartmail="gmail_signature"><div dir="ltr">IBS Center for Quantum Nanoscience<br>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>