<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear all,</div><div><br></div><div>Some time ago I read about the implementation of a constrained DFT approach here: <a href="http://theossrv1.epfl.ch/Main/OxidationStates">http://theossrv1.epfl.ch/Main/OxidationStates</a></div><div><br></div><div>I was naturally curious on how to simply test this and decided to use H-H dissociation as outlined in the text. However I am quite skeptical if my approach works and I was hoping someone would give it a look.</div><div><br></div><div>I repeat briefly the "algorithm" given in above link:</div><div><br></div><div>To the total energy calculated by QE add a penalty term -alpha*N_target</div><div>then the resulting "CDFT" energy term reads</div><div><br></div><div>E=E_KS+alpha[Tr[n]-N_target]</div><div>with "Tr[n]" the output of the lda_plus_u implementation, i.e. the occupation of the target manifold.</div><div><br></div><div>find the "correct" alpha by maximizing partial E/partial alpha -> alpha [Tr[n]-N_target]=0 (for alpha is not zero)</div><div><br></div><div>I attempted to do this as follows:</div><div><br></div><div>for alpha in `seq 1E-40 1.0 30` ;do</div><div>...</div><div>  lda_plus_u=.true.<br>  Hubbard_U(1)=1D-40<br>  Hubbard_alpha(1)=-$alpha<br>  Hubbard_U(2)=1D-40<br>  Hubbard_alpha(2)=$alpha</div><div>....</div><div>  H1   1.0  H.pbe-kjpaw_psl.1.0.0.UPF<br>  H2   1.0  H.pbe-kjpaw_psl.1.0.0.UPF</div><div>[...]</div><div><br></div><div>and I then collect for atom 1 and 2 Tr[n].</div><div><br></div><div>The main problem I encounter is "double counting" of electrons, i.e. the total number of electrons is not 2 but over 3 and the "most integer" occupation I can get for H1 is ~1.96 (ok, close enough to 2) but only for fairly large distance H-H. I attached a plot and the run-script, if anyone could give it a glance and tell me if this is totally wrong or a little bit right I would be very grateful!</div><div><br></div><div>Thanks in advance for your help!</div><div><br></div><div>Chris <br><div><img src="cid:ii_jn1lodbs1" alt="ECDFT.png" width="544" height="425"><br></div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br>-- <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>