<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">N.B., I should of course have said Hermitian, not symmetric.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">N.N.B., the extension to include spin-orbit is simple.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jul 31, 2019 at 11:24 PM Laurence Marks <<a href="mailto:laurence.marks@gmail.com">laurence.marks@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">A little math/analysis may be helpful -- with the caveat that this is my interpretation and I am currently on travel.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">The density matrix (ignoring spin-orbit for simplicity) "D" is a conjugate symmetric MxM matrix, M=5 for d electrons. From standard linear algebra it can be decomposed as</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">    D = U^T L U</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">where "U" is the MxM unitary matrix of eigenvectors, and "L" is the diagonal eigenvalue matrix. You can interpret the different eigenvectors as pseudo orbitals of the different p, d, f, and the eigenvalues as the occupancies. (Think about this please.) The total d (or p, f) occupancy is the sum of the eigenvalues which is equal to the trace of D.<br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">A small fortran, C or some other code such as matlab etc can be used to generate U & L.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">One can setup different density matrices by choosing different L values, and recreating D. As Fabien said before one would use this in a +U method to setup the values then run some cycles with fixed +U (i.e. -orbc) to establish the particular spin state. Then let it run normally to see if the spin state is stable. (One can also do this with -eece but it is trickier.) You can either do this choosing your L values by hand, or even using a random number generator to Monte Carlo explore the spin space. (I did this once for 4f electrons.)</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">Warning: this is not completely trivial and you need to know how to manipulate Wien2k, not just run it blackbox.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">N.B., I suggested some time ago adding the decomposition to (I think) lapwdm, but did not push. Any scientist who uses Wien2k should be capable of doing this analysis themselves.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jul 31, 2019 at 5:29 PM prasad jayasena <<a href="mailto:prasad.j86@yahoo.com" target="_blank">prasad.j86@yahoo.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear experts<br>
<br>
I am trying a test case to learn the correct procedure to check local minima.<br>
Dear experts<br>
<br>
I am sending this followup questions related to my previous post on " electron occupancy in dmatup/dn" in the mailing list. <br>
<br>
<br>
In the firt step I finished a "runsp_lapw" and saved. I checked case.dmatup/dn and they were empty.<br>
Then "runsp_lapw -orb -dm" is employed.<br>
<br>
the dmatup/dn both have values similar to the one in my previous email. I believe this work flow is corresponding to a random localization.<br>
<br>
I appreciate your support on following questions:<br>
<br>
(1)<br>
If I want to fix initial density matrix, is this the correct place to do it? Is the only file I have to edit case.dmatup/dn?<br>
<br>
Or,  do I have to do that without trying a random PBE+U as I have done? But in that case since the case.dmatup/dn from PBE is empty, I wonder how to create that file correctly.<br>
<br>
(2)<br>
If I need to include Spin orbital coupling, how I should proceed after this?<br>
<br>
(3)<br>
In my case I expected to have d-orbital electrons. But I do not see them in the case.dmatup/dn. Are there any special flow to find d-matrix in the output?<br>
<br>
Thank you in advance<br>
<br>
Prasad<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=JsiPLhhnGY5vWjIHSTaErXCmOUvmmcdRuaraV6JWK24&s=Ja--Kl05WF99HyMK2Quj7Vyc9Du2dPeJS3oF5nd85vA&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=JsiPLhhnGY5vWjIHSTaErXCmOUvmmcdRuaraV6JWK24&s=Ja--Kl05WF99HyMK2Quj7Vyc9Du2dPeJS3oF5nd85vA&e=</a><br>
SEARCH the MAILING-LIST at:  <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=JsiPLhhnGY5vWjIHSTaErXCmOUvmmcdRuaraV6JWK24&s=t-_GN5vjxe2Z9BCic3AZSLz1Rz8lTy19azdsfr7982E&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=JsiPLhhnGY5vWjIHSTaErXCmOUvmmcdRuaraV6JWK24&s=t-_GN5vjxe2Z9BCic3AZSLz1Rz8lTy19azdsfr7982E&e=</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail-m_8670004259231055703gmail_signature"><div dir="ltr">Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu/" target="_blank">www.numis.northwestern.edu</a><div>Corrosion in 4D: <a href="http://www.numis.northwestern.edu/MURI" target="_blank">www.numis.northwestern.edu/MURI</a><br>Co-Editor, Acta Cryst A<br>"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>Albert Szent-Gyorgi</div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr">Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu/" target="_blank">www.numis.northwestern.edu</a><div>Corrosion in 4D: <a href="http://www.numis.northwestern.edu/MURI" target="_blank">www.numis.northwestern.edu/MURI</a><br>Co-Editor, Acta Cryst A<br>"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>Albert Szent-Gyorgi</div></div></div>