Нужно рассчитать колебательный спектр при возбуждении из основного электронного состояние в первое возбужденное электронное состояние. При расчете этого спектра из основного колебательного состояния v"=0 проблем не возникает, так как тут многое по умолчанию и доп. команд не требуется.
Проблема начинается, когда я хочу посчитать спектр отдельно при возбуждении из состояния v"=1, v"=2, v"=3 и так далее. Я, вроде, пишу команды как советуют в справочном файле Гауссиана (k_freq.htm), но гауссиан упорно считает одно и тоже (из v"=0) Подробнее (молекула плоская 16 атомов, 42 колебательных нормальных моды): Я использую за основу пример из test902.com и добавляю команды из мануала:
Код: Выделить всё
%nproc=12
%mem=18GB
%Chk=test902-s0
#P B3LYP/cc-pVTZ Int=UltraFine Freq=SaveNM test
Gaussian Test job (Part 1):
ground state
0 1
C -1.942902 -0.541377 0.000000
................
--Link1--
%nproc=12
%mem=18GB
%Chk=test902-t1
#P B3LYP/cc-pVTZ Int=UltraFine Freq=(SelectNormalModes, SaveNM) test
Gaussian Test job (Part 2):
excited state
0 1
C -1.923461 -0.546519 0.000000
.........
1-42
--Link1--
%nproc=12
%mem=18GB
%Chk=test902-s0
%nosave
#P Freq=(ReadFC,FC,ReadFCHT,VibRot,SelFCModes) Geom=Check NoSymm
Gaussian Test job 902 (Part 3):
gr state -> ex state
0,1
SPECHWHM=1 SPECRES=0.1 PRTMAT=12 PRTINT=0.01 SPECMIN=-4000. SPECMAX=+6000.
test902-t1.chk
2-2в результатах расчета. Заранее спасибо.
