<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Dear Kosuke Nakano,<br>
  the total charge option acts on deciding the number of states
used and how to fix the fermi energy/chemical potential.<br>
  the occupation matrix return the projection of the occupied
manifold (determined by the previous recipe) on the localized
states so it should be what one expects it to be.<br>
  so yes it looks like it's Ni ~ +2 (d8) but with partial
occupation of the eg orbitals and NOT spin polarized... Is this
what you expect/desire ?.<br>
Did you suggest a starting magnetization to see it the system does
not prefer being d5_up d3_dw ?<br>
<br>
stefano<br>
<br>
On 13/01/2017 05:26, <a class="moz-txt-link-abbreviated" href="mailto:kosuke.nakano@agc.com">kosuke.nakano@agc.com</a> wrote:<br>
</div>
<blockquote
cite="mid:89a85045f3084b75a5e3fc93fc1dd50e@T180MLZPEXCMB3.agc.jp"
type="cite">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<meta name="Generator" content="Microsoft Word 15 (filtered
medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"ï¼ï¼³ ゴシック";
panose-1:2 11 6 9 7 2 5 8 2 4;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:"\@ï¼ï¼³ ゴシック";
panose-1:2 11 6 9 7 2 5 8 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0mm;
margin-bottom:.0001pt;
text-align:justify;
text-justify:inter-ideograph;
font-size:10.5pt;
font-family:"Arial","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.17
{mso-style-type:personal-compose;
font-family:"Arial","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Arial","sans-serif";}
/* Page Definitions */
@page WordSection1
{size:612.0pt 792.0pt;
margin:99.25pt 30.0mm 30.0mm 30.0mm;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026">
<v:textbox inset="5.85pt,.7pt,5.85pt,.7pt" />
</o:shapedefaults></xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Hi, all:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p>Â </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I have a question about
tot_charge flag with DFT+U calculations<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p>Â </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I would like to
calculate d state of transition metal with tot_charge flag.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p>Â </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I know that tot_charge
flag with DFT+U calculations is valid as written in
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><a class="moz-txt-link-rfc2396E" href="http://qe-forge.org/pipermail/pw_forum/2015-August/107486.html"><http://qe-forge.org/pipermail/pw_forum/2015-August/107486.html></a><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p>Â </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I am wondering if
occupation matrices are also valid in tot_charge
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">flag with DFT+U
calculations.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p>Â </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">For example, I obtained
the following occupation matrices with a finite
tot_charge(+4)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The system has 24 atoms
and 4 of them are Ni atoms.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p>Â </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">--The final result--<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">atom  1  Tr[ns(na)]
(up, down, total) =Â Â 4.26150Â 4.08319Â 8.34469 <- Ni
(+2, d8?)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  spin 1<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â Â Â eigenvalues: <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â Â 0.658Â 0.658Â 0.981Â
0.982Â 0.982<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â Â Â eigenvectors:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.000Â 0.000Â 0.996Â
0.003Â 0.002<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.588Â 0.074Â 0.000Â
0.291Â 0.046<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.074Â 0.588Â 0.001Â
0.046Â 0.291<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.038Â 0.300Â 0.002Â
0.089Â 0.571<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.300Â 0.038Â 0.001Â
0.571Â 0.090<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â Â Â occupations:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.981 -0.000 -0.000
-0.000 -0.000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">-0.000Â 0.768Â 0.000
-0.000Â 0.153<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">-0.000Â 0.000Â 0.768Â
0.153 -0.000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">-0.000 -0.000Â 0.153Â
0.873 -0.000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">-0.000Â 0.153 -0.000
-0.000Â 0.873<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  spin 2<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â Â Â eigenvalues: <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â Â 0.649Â 0.649Â 0.923Â
0.931Â 0.931<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â Â Â eigenvectors:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.000Â 0.000Â 1.000Â
0.000Â 0.000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.649Â 0.018Â 0.000Â
0.078Â 0.255<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.018Â 0.649Â 0.000Â
0.255 Â 0.078<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.009Â 0.324Â 0.000Â
0.510Â 0.157<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.324Â 0.009Â 0.000Â
0.157Â 0.510<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â Â Â occupations:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.923Â 0.000Â 0.000Â
0.000Â 0.000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.000Â 0.743Â 0.000Â
0.000Â 0.133<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.000Â 0.000Â 0.743Â
0.133Â 0.000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.000Â 0.000Â 0.133Â
0.837Â 0.000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Â 0.000Â 0.133Â 0.000Â
0.000Â 0.837<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">atomic mag. moment =Â Â Â Â
0.178309<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p>Â </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I am wondering if it
means Ni(+2, d8) or not.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p>Â </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Best regards<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p>Â </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Kosuke Nakano<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Asahi Glass Co., Japan.<o:p></o:p></span></p>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Pw_forum mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>
<a class="moz-txt-link-freetext" href="http://pwscf.org/mailman/listinfo/pw_forum">http://pwscf.org/mailman/listinfo/pw_forum</a></pre>
</blockquote>
<p><br>
</p>
</body>
</html>