<div dir="ltr"><span style="font-size:12.8px;white-space:nowrap">Guido Fratesi,</span><div><span style="font-size:12.8px;white-space:nowrap"><br></span></div><div><span style="font-size:12.8px;white-space:nowrap">Your answer is too hard even for me. Remember what beginners do not spend their lives not only for pressing the buttons, usually people try to engage in much more useful things, like experiment)</span><br><div><span style="font-size:12.8px;font-weight:bold;white-space:nowrap"><br></span></div><div><span style="font-size:12.8px;font-weight:bold;white-space:nowrap"><br></span></div><div class="gmail_extra"><br><div class="gmail_quote">2016-08-10 15:25 GMT+04:00 Guido Fratesi <span dir="ltr"><<a href="mailto:fratesi@mater.unimib.it" target="_blank">fratesi@mater.unimib.it</a>></span>:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Rolly,<br>
<br>
let me try to expand a bit my previous concise answers to your questions<br>
(Q1-Q2):<br>
<br>
As it is implemented, molecularpdos.x needs the two systems (full &<br>
part) to be described by the very same parameters as of:<br>
(1) unit cell dimensions<br>
(2) k-point sampling<br>
(3) number spin populations considered<br>
(4) dimension of the localized basis set (atomic functions)<br>
<br>
In principle, such close correspondence could be avoided in some cases<br>
like where your reference (i.e., part) system is an isolated molecule.<br>
Its states could be computed at the gamma point for a cell of your<br>
choice, then the Bloch states (as would be obtained by computing them<br>
for the same cell as for the full system) could be reconstructed<br>
analytically. This is not implemented so far, so an explicit calculation<br>
with the very same k-grid is needed (I generally copy the k-points as<br>
given in output by pw.x for the full system).<br>
<br>
Similarly, the spin-up states of the full system will be projected on<br>
the spin-up states of the part (same for spin down); this requires you<br>
to make a calculation with nspin=2 (but no magnetization) also when the<br>
molecule by itself is not spin-polarized if it is embedded in a<br>
spin-polarized system.<br>
<br>
Smearing / not smearing is not relevant to molecularpdos.x : just use<br>
the one adequate to each of the two calculations (full & part).<br>
<br>
As for the dimension of the localized basis set, molecularpdos.x assumes<br>
that the atomic orbitals of the full & part calculation correspond 1-1<br>
to each other, keeping the same order. The range of the orbitals to be<br>
considered is specified by input variables<br>
i_atmwfc_{beg,end}_{full,part} so the atoms to be considered have to be<br>
consecutive and listed in the same order in the full & part systems.<br>
BEWARE (common mistake) this range specifies the atomic wavefunctions<br>
(see the output of projwfc.x for a list), not the atoms. You may want to<br>
add atoms to the "part" calculation, e.g., add H to saturate a dangling<br>
bond: in that case, just do not include the corresponding atomic states<br>
in the range for molecularpdos.x.<br>
<br>
Hope this helps,<br>
Guido<br>
<div><div class="h5"><br>
On 09/08/2016 18:46, Rolly Ng wrote:<br>
> Dear QE experts,<br>
><br>
> I am trying to refine computation of 3-NT adsorption on metal doped<br>
> graphene using QE. This is according to our previous work<br>
> <a href="http://dx.doi.org/10.1016/j.commatsci.2011.07.045" rel="noreferrer" target="_blank">http://dx.doi.org/10.1016/j.<wbr>commatsci.2011.07.045</a><br>
><br>
> With the help of Dr. Guido Fratesi, I am exploring molecularpdos.x to<br>
> find the change of HOMO and LUMO of the nitrated tyrosine molecule on<br>
> adsorption to metal-doped graphene.<br>
><br>
> There are my questions and answers, your comments are welcome.<br>
><br>
> Q1) Spin polarization. The full system contains a metal doped graphene<br>
> with single metallic atom of Au and Ni at the center and a 3-NT molecule<br>
> adsorbed onto the graphene sheet, so I included nspin=2 for the full<br>
> system. But the adsorbate (3-NT) molecule is likely to be nonmagnetic,<br>
> so I did nspin=1 for its gas phase. Can molecularpdos.x cope with<br>
> different nspin for the full system and adsorbate?<br>
> A1) Yes, since nspin=2 can handle a nonmagnetic case<br>
> R1) So, I am adding nspin=2 to the 3-NT molecules<br>
><br>
> Q2) K-points. I used smearing for the full system since it is<br>
> semi-metallic. Should I use smearing for the molecule in gas phase? I<br>
> believe this is not a good idea but can molecularpdos.x works with<br>
> k-point and non k-point?<br>
> A2) Smearing: no problem, K-points: should be the same<br>
> R2) So, I am using the same K-points for both full system and single<br>
> molecule.<br>
><br>
> Q3) I would like to evaluate the change of the HOMO-LUMO gap of the<br>
> adsorbate (3-NT) on adsorption to the metal-doped graphene. Can<br>
> molecularpdos.x do that?<br>
> A3) Dr. Fratesi and his team used that code also to study<br>
> molecules/graphene...<br>
> <a href="http://dx.doi.org/10.1038/srep24603" rel="noreferrer" target="_blank">http://dx.doi.org/10.1038/<wbr>srep24603</a><br>
> R3) Very useful and thank you!<br>
><br>
> Q4) I would like to determine if physisorption or chemisorption occurs as<br>
> with DMol3, can molecularpdos.x do the same job? For the physisorption<br>
> vs. chemisorption, the adsorption energy may be a good indicator.<br>
> However, I believe the electronic structure of the bond may also<br>
> indicate which type of adsorption it suppose to be?<br>
><br>
> With regards,<br>
> Rolly<br>
><br>
<br>
--<br>
</div></div>Guido Fratesi<br>
<br>
Dipartimento di Fisica<br>
Universita` degli Studi di Milano Italy<br>
<div class=""><div class="h5"><br>
______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Best wishes,</div><div>Maxim Arsent'ev, Ph.D. (Chemistry)</div><div>Laboratory of research of nanostructures</div><div>Institute of Silicate Chemistry of RAS</div></div></div></div></div>
</div></div></div>