<div dir="ltr"><div>You are correct, the code is not. Once upon a time rho(:,1)=spin up, rho(:,2)=spin down components, but then we moved to a different logic: rho(:,1)=total charge, rho(:,2)=spin polarization. Apparently that piece of code was forgotten. Thank you for reporting it<br></div><div><br></div><div>Paolo<br></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jan 16, 2021 at 8:31 AM 杨 冬 <<a href="mailto:yangdong297056853@hotmail.com">yangdong297056853@hotmail.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 style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hi, all</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
 </div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I'm trying to add DORI analysis(<span style="background-color:rgb(255,255,255);display:inline">DOI: 10.1021/ct500490b</span>) and other functions into pp.x code(<span style="background-color:rgb(255,255,255);display:inline">QE 6.7</span>), 
 but I'm not sure that rho%of_r(:,1)(<span style="background-color:rgb(255,255,255);display:inline">from PP/elf.f90</span>) is the total electron density or the spin-up electron density. In my opnions, <span style="background-color:rgb(255,255,255);display:inline">rho%of_r(:,1)
 means <span style="background-color:rgb(255,255,255);display:inline">total electron density(by reading PP/punchplof.90) and there is no need to put the total (up+down) charge density in rho%of_r(:,1) in the subroutine do_elf and do_sl2rho. I am
 not sure whether I understand it correctly.</span></span><span style="background-color:rgb(255,255,255);display:inline"></span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="background-color:rgb(255,255,255);display:inline"><span style="background-color:rgb(255,255,255);display:inline"><br>
</span></span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Best regards,
<div><br>
</div>
<div>Dong Yang</div>
<div><br>
</div>
<div>JSG</div>
</div>
</div>

_______________________________________________<br>
developers mailing list<br>
<a href="mailto:developers@lists.quantum-espresso.org" target="_blank">developers@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/developers" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/developers</a><br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>