<div dir="ltr">Dear researchers and developers,<div>I am a beginner in CP2K, and I major in computational mathematics. So the concepts about chemistry are new to me. I want to know how to get the HOMO-LUMO gap in a QMMM simulation. The setting in my input file is as follows</div><div><br></div><div><div>&PRINT</div><div>         &MO_CUBES</div><div>             NHOMO -1</div><div>             NLUMO -1</div><div>             WRITE_CUBE F</div><div>         &END MO_CUBES</div></div><div><br></div><div>which locates at the &DFT subsection. </div><div><br></div><div>But the results seems unreasonable, the result of 2-QM-atom is even larger than the situation if I use full QM calculation. The whole system consists of 48 atoms. But if I specify 12 atoms as QM-atom, the result is less than the result of full QM calculation. So, is there any error in the setting ? </div><div><br></div></div>