<div dir="ltr">Dear CP2K experts,<div><br></div><div>I am running a system with MO print out and would like to check on the changes of HOMO and LUMO energy of my system.</div><div>However, i only managed to get the LUMO energy as per below.</div><div><br></div><div><div>Â Lowest Eigenvalues of the unoccupied subspace spin Â Â Â Â Â Â 1</div><div>Â -----------------------------------------------------</div><div>Â Reached convergence in Â Â Â Â Â 119 Â iterations</div><div>Â Â Â Â 0.02037228 Â Â Â 0.05911493</div><div><br></div><div>Â Lowest Eigenvalues of the unoccupied subspace spin Â Â Â Â Â Â 2</div><div>Â -----------------------------------------------------</div><div>Â Reached convergence in Â Â Â Â Â 89 Â iterations</div><div>Â Â Â Â 0.05501903 Â Â Â 0.08163833</div><div><br></div><div>Â HOMO - LUMO gap [eV] : Â Â 0.554358</div><div>Â HOMO - LUMO gap [eV] : Â Â 1.497144</div></div><div><br></div><div>is the HOMO-LUMO gap calculated by using the eigenvalue? for instance (eigenvalue of HOMO) - (eigenvalue of LUMO 0.02037228)Â </div><div><b>May i know is it possible to print out the value of HOMO and LUMO for the gap as calculated?</b></div><div><br></div><div>below is the part of my input file.</div><div><br></div><div><div>Â Â Â Â &DFT</div><div>Â Â Â Â Â Â Â Â BASIS_SET_FILE_NAME Â /share/apps/cp2k/cp2k-2.6.1/data/BASIS_MOLOPT</div><div>Â Â Â Â Â Â Â Â POTENTIAL_FILE_NAME Â /share/apps/cp2k/cp2k-2.6.1/data/GTH_POTENTIALS</div><div>Â Â Â Â Â Â Â Â MULTIPLICITY 2</div><div>Â Â Â Â Â Â Â Â CHARGE -1</div><div>Â Â Â Â Â Â Â Â ROKS ON</div><div>Â Â Â Â Â Â Â Â &SIC</div><div>Â Â Â Â Â Â Â Â Â Â Â Â ORBITAL_SET UNPAIRED</div><div>Â Â Â Â Â Â Â Â Â Â Â Â SIC_METHOD MAURI_SPZ</div><div>Â Â Â Â Â Â Â Â Â Â Â Â SIC_SCALING_A 0.3</div><div>Â Â Â Â Â Â Â Â Â Â Â Â SIC_SCALING_B 0.2</div><div>Â Â Â Â Â Â Â Â &END SIC</div><div>Â Â Â Â Â Â Â Â &QS</div><div>Â Â Â Â Â Â Â Â Â Â Â Â EPS_DEFAULT 1.0E-10</div><div>Â Â Â Â Â Â Â Â &END QS</div><div>Â Â Â Â Â Â Â Â &MGRID</div><div>Â Â Â Â Â Â Â Â Â Â Â Â NGRIDS 5</div><div>Â Â Â Â Â Â Â Â Â Â Â Â CUTOFF 280</div><div>Â Â Â Â Â Â Â Â Â Â Â Â REL_CUTOFF 40</div><div>Â Â Â Â Â Â Â Â &END MGRID</div><div>Â Â Â Â Â Â Â Â &XC</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &XC_FUNCTIONAL PBE</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â &PBE</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â Â PARAMETRIZATION REVPBE</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â &END PBE</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &END XC_FUNCTIONAL</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &VDW_POTENTIAL</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â &PAIR_POTENTIAL</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â Â TYPE DFTD3(BJ)</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â Â PARAMETER_FILE_NAME /share/apps/cp2k/cp2k-2.6.1/data/dftd3.dat</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â Â D3BJ_SCALING 1.000 0.5238 2.3550 3.5016</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â Â CALCULATE_C9_TERM OFF</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â &END PAIR_POTENTIAL</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &END VDW_POTENTIAL</div><div>Â Â Â Â Â Â Â Â &END XC</div><div>Â Â Â Â Â Â Â Â &SCF</div><div>Â Â Â Â Â Â Â Â Â Â Â Â SCF_GUESS RESTART</div><div>Â Â Â Â Â Â Â Â Â Â Â Â EPS_SCF 5.0E-6</div><div>Â Â Â Â Â Â Â Â Â Â Â Â MAX_SCF 100</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &OT ON</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â MINIMIZER DIIS</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â ROTATION ON</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &END OT</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &MIXING T</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â METHOD BROYDEN_MIXING</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â ALPHA 0.4</div><div>Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â NBROYDEN 8</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &END MIXING</div><div>Â Â Â Â Â Â Â Â &END SCF</div><div>Â Â Â Â Â Â Â Â &POISSON</div><div>Â Â Â Â Â Â Â Â Â Â Â Â PERIODIC NONE</div><div>Â Â Â Â Â Â Â Â Â Â Â Â POISSON_SOLVER MT</div><div>Â Â Â Â Â Â Â Â &END POISSON</div><div>Â Â Â Â Â Â Â Â &PRINT</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &MO</div><div>Â Â Â Â Â Â Â Â Â Â Â Â EIGENVALUES</div><div>Â Â Â Â Â Â Â Â Â Â Â Â OCCUPATION_NUMBERS</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &END MO</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &MO_CUBES</div><div>Â Â Â Â Â Â Â Â Â Â Â Â NHOMO 2</div><div>Â Â Â Â Â Â Â Â Â Â Â Â NLUMO 2</div><div>Â Â Â Â Â Â Â Â Â Â Â Â &END MO_CUBES</div><div>Â Â Â Â Â Â Â Â &END PRINT</div><div>Â Â Â Â &END DFT</div><div>&END FORCE_EVAL</div></div><div><br></div><div>Rgds</div><div>Alfred</div><div><br></div></div>