Questions about printing out unoccupied molecular orbitals (LUMO, LUMO+x)

[zh]

In my test calculations, I also found another problem:
If the keyword "ADDED_MOS" in the section &SCF is set as 0 or 1, the
LUMO and LUMO+1 orbitals of C6H6 molecule obtained with cp2k will be
similar to those obtained with other codes.

In the log file, there are two values printed out for the LUMO-HOMO
gap.

The first value of LUMO-HOMO gap is printed out due to the following
keyword:
"
    &MO
        EIGENVALUES
        OCCUPATION_NUMBERS
      &END MO
"
This one is reasonable. That is why I say that the unoccupied
molecular orbitals are already obtained before calling "make_lumo()".

The second value of LUMO-gap is printed out due to the following
keyword:
"
&MO_CUBES
         NHOMO  3
         NLUMO  10
         WRITE_CUBE .TRUE.
      &END MO_CUBES
"
While the second one is quite larger than the first one.