<div dir="ltr">If I'm not wrong, for spin polarized calculations the kpoints are doubled. Each subset is assigned to a given spin. In this case, to visualize your HOMO and LUMO you should carry out<div>a run with </div><div>kpoint(1) = 1</div><div>kpoint(2) = 2<br>that will double the number of output wave functions. Among these, you should consider the ones that correspond to your HOMO and your LUMO.</div><div>Giovanni</div><div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(0,0,0)">-- <br><br>Giovanni Cantele, PhD<br>CNR-SPIN<br>c/o Dipartimento di Fisica<br>Universita' di Napoli "Federico II"<br>Complesso Universitario M. S. Angelo - Ed. 6<br>Via Cintia, I-80126, Napoli, Italy<br><a href="mailto:giovanni.cantele@spin.cnr.it" style="color:rgb(17,85,204)" target="_blank">e-mail: giovanni.cantele@spin.cnr.it</a><br>Phone: +39 081 676910<br>Skype contact: giocan74<br><br>ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" style="color:rgb(17,85,204)" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>Web page: </span><a href="https://sites.google.com/view/giovanni-cantele/home" style="color:rgb(17,85,204)" target="_blank">https://sites.google.com/view/giovanni-cantele/home</a><br></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mar 15 ott 2024 alle ore 23:18 Prem Prakash Sahu <<a href="mailto:premprakash.sahu@unifi.it">premprakash.sahu@unifi.it</a>> ha scritto:<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>Dear QE Community,</div><div><br></div><div>I am currently working 
on a transition metal-porphyrin system where vanadyl (V=O) is ligated to
 a porphyrin molecule, in an isolated state, using Quantum ESPRESSO. I 
am performing calculations at the gamma point only. After running an SCF
 calculation, I observed the following:</div><div><br></div>Spin-up: HOMO-LUMO: 178-179<br><div>Spin-down: HOMO-LUMO: 177-178.</div><div><br></div><div>I
 wanted to visualize the individual molecular orbitals to inspect how my
 alpha (spin-up) and beta (spin-down) molecular orbitals (MOs) look. For
 this, I ran a postprocessing step using `pp.x` with the following 
input:</div><div><br></div>***********************************************************************************************************************************************<br>&INPUTPP<br>    outdir = './out'<br>    prefix = 'SP_VO-por_chrom_PBE_vdW+U_111'<br>    filplot   = 'SP_VO-por_chrom_PBE_vdW+U_111_MO.dat',<br>    plot_num  = 7,<br>    kpoint    = 1,<br>    kband(1)  = 177,<br>    kband(2)  = 179,<br>    lsign     = .true.,<br>/<br>&PLOT<br>    fileout       = 'SP_VO-por_chrom_PBE_vdW+U_111_MO.cube',<br>    iflag         = 3,<br>    nfile         = 1,<br>    output_format = 6,<br>    weight(1)     = 1.0,<br>    filepp(1)     = 'SP_VO-por_chrom_PBE_vdW+U_111_MO.dat'<br>/<br>**************************************************************************************************************************************************<br><br><div>However,
 I did not obtain the expected individual orbitals for the specified 
range (177-179) in the spin-up and spin-down channels. From the 
documentation, I noticed the `spin_component` option for `plot_num = 7`,
 but it seems to deal with directional magnetization, which is not what I
 need.</div><div><br></div><div>I must plot these orbitals individually 
because I need to confirm the correct d-orbital assignment on the 
vanadium atom, specifically to crosscheck the spin-up and spin-down 
HOMO-LUMO orbitals. Having this information will help ensure that the 
electronic structure is correctly represented.</div><div><br></div><div> Could
 anyone kindly suggest how I can resolve this or point me in the right 
direction to achieve the individual plots for these orbitals?</div><div><br></div><div>Any guidance would be much appreciated.</div><div><br></div>Best regards,  <font color="#888888"><br>Prem Prakash Sahu<br>University of Florence</font><br clear="all"><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><br></div></div></div></div></div>
_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>