<div dir="ltr">Typically the 'frontier orbitals' - those close to the HOMO or LUMO are most important in deciding the system's properties. So those are what you want to examine.<div><br></div><div>But you should choose the orbitals that you want to look at. You need to decide what information you need and change NHOMO and NLUMO appropriately.</div><div><br></div><div>Matt<br><br>On Monday, March 25, 2019 at 10:33:13 AM UTC, Akshay Malik wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div dir="ltr"><div>Dear cp2k users,</div><div>I have done AIMD of liquid water. I want to know what the following section is doing- <br></div><div>    &PRINT<br>       ! at the end of the SCF procedure generate cube files of the density<br>       &E_DENSITY_CUBE OFF<br>       &END E_DENSITY_CUBE<br>       ! compute eigenvalues and homo-lumo gap each 10nd MD step<br>       &MO_CUBES<br>          NLUMO 4<br>          NHOMO 4<br>          WRITE_CUBE .FALSE.<br>          &EACH<br>            MD 10<br>          &END<br>       &END<br>    &END<br></div><div>How do we choose this NLUMO and NHOMO?    <br></div><div><br></div></div></div>
</blockquote></div></div>