Расчёты возбуждённых состояний полуэмпирикой

[Vit Nhoc]

Я пробую в Gaussian рассчитывать возбуждённые состояния бодипаев полуэмпирикой. Нам нужны не только спектры поглощения, но и свойства возбуждённых состояний - геометрия, ППЭ (мы пытается считать флюоресценцию).
Я запустил на G09A на AM1 расчёт оптимизации геометрии возбуждённого состояния (#P AM1 TD(NSTATES=6, ROOT=1)), и он оборвался: первую точку рассчитал, напечатал спектр поглощения, и прервался с ошибкой. Ниже в этом сообщении – входной файл и что выдал Gaussian. Подскажите, как запустить этот расчёт.
G03 и G98 не рассчитали даже спектр.
Может быть, Gaussian не может рассчитать первые производные? Я что-то не нашёл в помощи, как сделать численный расчёт первых производных.
И ещё вопрос общего плана. Чем отличается возбуждённое состояние от состояния, например, с мультиплетностью 3 (у меня основное состояние – обычный синглет)? Расчёт TD даёт список состояний с другой мультиплетностью? А в спектре поглощения может быть переход из синглета в триплет?

Входной файл:

#P AM1 OPT(LOOSE) TD(NSTATES=6, ROOT=1)

BODIPY

0 1
6 0.000000000 0.000000000 -5.213443000
6 0.799543000 0.905233000 -4.515222000
6 0.806775000 0.904269000 -3.121838000
6 0.000000000 0.000000000 -2.415208000
6 0.000000000 0.000000000 -0.953499000
6 -2.569243000 0.149779000 -0.726258000
6 -1.207479000 0.057751000 -0.243367000
7 -1.257454000 0.055850000 1.174782000
5 0.000000000 0.000000000 2.129864000
9 0.046612000 1.058017000 2.903742000
9 -0.046612000 -1.058017000 2.903742000
7 1.257454000 -0.055850000 1.174782000
6 -3.403060000 0.194465000 0.390581000
6 -2.549515000 0.135211000 1.554527000
6 1.207479000 -0.057751000 -0.243367000
6 2.569243000 -0.149779000 -0.726258000
6 3.403060000 -0.194465000 0.390581000
6 2.549515000 -0.135211000 1.554527000
6 -0.799543000 -0.905233000 -4.515222000
6 -0.806775000 -0.904269000 -3.121838000
1 0.000000000 0.000000000 -6.313754000
1 1.427928000 1.622676000 -5.063980000
1 1.436719000 1.619724000 -2.569983000
1 -1.427928000 -1.622676000 -5.063980000
1 -1.436719000 -1.619724000 -2.569983000
1 2.849087000 -0.183926000 -1.776740000
1 -2.849087000 0.183926000 -1.776740000
1 -4.486094000 0.263212000 0.417497000
1 4.486094000 -0.263212000 0.417497000
1 -2.848229000 0.151298000 2.607420000
1 2.848229000 -0.151298000 2.607420000


В выходном файле:

dumping /fiocom/, unit = 3 NFiles = 1 SizExt = 524288 WInBlk = 512
defal = T LstWrd = 67072 FType=2 FMxFil=10000

Number 0
Base 20480
End 67072
End1 67072
Wr Pntr 20480
Rd Pntr 20480
Length 46592
Error termination in NtrErr:
NtrErr Called from FileIO.

[Гесс]

Я запустил инпут у себя.

(Enter /usr/local/g09/l607.exe)
******************************Gaussian NBO Version 3.1******************************
N A T U R A L A T O M I C O R B I T A L A N D
N A T U R A L B O N D O R B I T A L A N A L Y S I S
******************************Gaussian NBO Version 3.1******************************
/RESON / : Allow strongly delocalized NBO set

Analyzing the SCF density

Job title: BODIPY

Storage needed: 27848 in NPA, 34362 in NBO (2097146885 available)

Subroutine NAOANL could not find a S-type valence orbital on atom C 1.
IVAL : 1 1 0 0 M : 1 LA : 1
Error
Error termination via Lnk1e in /usr/local/g09/l607.exe at Fri Nov 6 16:14:11 2015.
Job cpu time: 0 days 0 hours 2 minutes 56.9 seconds.
File lengths (MBytes): RWF= 57 Int= 0 D2E= 0 Chk= 6 Scr= 1
170.959u 6.040s 0:18.86 938.4% 0+0k 205256+188312io 1523pf+0w

Я считал спектры поглощения в AM1 в VAMP, там переходы синглет-триплет нулевой интенсивности.
У вас довольно небольшая молекула - вы не думали делать ее в TDDFT?

[Гесс]

ему nbo с полуэмпирикой не нравится.
дайте ему

#P AM1 OPT TD=(NSTATES=6,ROOT=1) Pop=none

расчет заканчивается корректно, делает ли он то что от него хотят - не проверял

[Vit Nhoc]

На моём Gaussian-е пока это не проходит, ну да ладно. Я всё-таки хотел бы понять: чем отличается возбуждённое состояние от состояний с мультиплетностью 3, 5, 7 и т.д.?

[Гесс]

На моём Gaussian-е пока это не проходит, ну да ладно. Я всё-таки хотел бы понять: чем отличается возбуждённое состояние от состояний с мультиплетностью 3, 5, 7 и т.д.?

файлы дошли?
Заселенность орбиталей
Основное состояние 2 2 2 0 0 0. Все остальные состояния являются возбужденными.
Триплет 2 2 1 1 0 0 причем 1 и 1 параллельны. ЕМНИП диаграмма яблонского запрещает его образование при поглощении фотона из основного синглетного состояния. Поэтому в спектре поглощения переходов синглет-триплет быть не должно. Тут пусть лучше спектроспециалисты разложат на пальцах.
Возбужденный синглет либо 2 2 1 1 0 0 (причем 1 и 1 антипараллельны) либо 2 2 0 2 0 0 если вверх подняли сразу пару.
Квинтет 2 1 1 1 1 0