<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:15.0150146484375px">Dear PWSCF users.</span><div style="font-family:arial,sans-serif;font-size:15.0150146484375px"><br></div><div style="font-family:arial,sans-serif;font-size:15.0150146484375px">
<div>I'm now try to fixed the valence of one atom in the cell.</div><div><br></div><div>For example, how we impose one of the Nb atom as (+4) valence </div><div><br></div><div>(one electron is localized in Nb atom   Nb(+5) is ideal condition )</div>
<div><br></div><div>in the crystal.</div><div><br></div><div>I tried to use  constrained_magnetization = 'atomic' </div><div><br></div><div>with statring_magnetization  like below,</div><div><br></div><div>(Here atom index 2 is Nb atom which I set as valence (+4, one electron </div>
<div><br></div><div>localized in Nb atom, 3 is Nb atom with ideal +5 valence )</div><div><br></div><div><div>   starting_magnetization(1) = 0.0</div><div>    starting_magnetization(2) = 1.0</div><div>    starting_magnetization(3) = 0.0</div>
<div>    starting_magnetization(4) = 0.0</div><div>    starting_magnetization(5) = 0.0 </div><div>    lda_plus_u =.TRUE.</div><div>    lda_plus_u_kind = 0</div><div>    Hubbard_U(2)=4.0</div><div>    Hubbard_U(3)=4.0</div>
<div>    tot_charge = -1.0</div><div>    tot_magnetization = 1.0 </div><div>    constrained_magnetization = 'atomic'</div><div>    lambda = 1.0    </div></div><div><br></div><div>But, the results was some of Nb atoms has a little excess electron (, </div>
<div>but smaller than one ), and I couldn't get +4 Nb which I aimed to </div><div>make.</div><div><br></div><div>How to treat to impose the localized electron in one of the atoms </div><div><br></div><div>in crystals.</div>
<div><br></div><div>If I can get reply, I'm very happy.<br></div><div><br></div><div>Sincerely,</div><div><br></div><div>Yukihiro Okuno.</div></div><div><br></div>-- <br><div dir="ltr"><div>　</div></div>
</div>