В Gaussian - эта информация уже выдаётся по умолчанию в стандартном выходном файле, никаких дополнительных ключевых слов писать не надо.
Смотрите здесь (заголовок одинаковый во всех файлах):
**********************************************************************
Population analysis using the SCF density.
**********************************************************************
Alpha occ. eigenvalues -- -10.24526 -10.23438 -10.23438 -10.21601 -10.18354
Alpha occ. eigenvalues -- -10.18353 -10.17911 -10.17911 -10.16531 -10.16531
Alpha occ. eigenvalues -- -10.15957 -10.15956 -6.78162 -6.78162 -6.77260
Alpha occ. eigenvalues -- -6.77260 -6.76466 -6.76466 -6.75676 -6.75675
Alpha occ. eigenvalues -- -0.87860 -0.86543 -0.86504 -0.84445 -0.77296
Alpha occ. eigenvalues -- -0.75749 -0.75696 -0.74067 -0.58446 -0.56484
Alpha occ. eigenvalues -- -0.56293 -0.54386 -0.52109 -0.50858 -0.50356
Alpha occ. eigenvalues -- -0.49352 -0.49297 -0.46147 -0.46084 -0.44863
Alpha occ. eigenvalues -- -0.43465 -0.41233 -0.41215 -0.40574 -0.40196
Alpha occ. eigenvalues -- -0.38600 -0.38219 -0.36018 -0.34780 -0.34234
Alpha occ. eigenvalues -- -0.32401 -0.32035 -0.31268 -0.30980 -0.28379
Alpha occ. eigenvalues -- -0.28133 -0.28070 -0.25616 -0.24927 -0.20868
Alpha virt. eigenvalues -- -0.08174 -0.03664 -0.02988 -0.02225 -0.00212
Alpha virt. eigenvalues -- 0.01170 0.01514 0.02059 0.03472 0.03928
Alpha virt. eigenvalues -- 0.04769 0.06459 0.07065 0.07148 0.07608
Alpha virt. eigenvalues -- 0.07760 0.08577 0.08656 0.08869 0.08910
Alpha virt. eigenvalues -- 0.09781 0.10228 0.10690 0.11123 0.11573
Alpha virt. eigenvalues -- 0.12481 0.12640 0.13083 0.14026 0.15505
Alpha virt. eigenvalues -- 0.16236 0.16568 0.18606 0.19197 0.19957
Alpha virt. eigenvalues -- 0.20099 0.22010 0.22163 0.22979 0.24732
Alpha virt. eigenvalues -- 0.24915 0.25028 0.25488 0.25738 0.26332
Alpha virt. eigenvalues -- 0.26399 0.27039 0.27377 0.27381 0.28886
Alpha virt. eigenvalues -- 0.29202 0.29586 0.31008 0.31132 0.31476
Alpha virt. eigenvalues -- 0.31613 0.32723 0.34160 0.35172 0.37465
Это список энергии всех молекулярных орбиталей,
в электронвольтах (написал с утра, не проснувшись

конечно же,
в а.е.), отсортированный снизу вверх. Находите то место, где "occ" (occupied) превращается в "virt" (virtual) - и находите нужные вам числа.
P.S.: Только следите за тем,
к какому шагу геометрической оптимизации относятся прочитанные вами числа!