Расчеты в gaussian с учетом растворителя... Проблемы.....

[Vladchimic]

Входной файл
%chk=kation.chk
%mem=2100Mb
%nproc=2
#T B3LYP/6-31+G(d,p) SCRF(SCIPCM) Opt

и т.д.
Проводится оптимизация геометрии и расчет энергии положительно заряженной частицы (C8H15N2) в ацетонитриле (e=36.64). После долгий раздумий получаю

The electronic state of the initial guess is 1-A.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
>>>>>>>>>> Convergence criterion not met.
SCF Done: E(RB+HF-LYP) = -423.288866741 A.U. after 129 cycles
Convg = 0.2125D-04 -V/T = 2.0101
S**2 = 0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in /install/g03/l502.exe at Thu Nov 16 16:44:30 2006.

Может ли кто-нибудь посказать мне, в чем дело и что может мне помочь ? Спасибо....

[Nord]

Процедура ССП у Вас не сошлась

[Vladchimic]

Очень лаконично. Спасибо.
Что же можно или нужно предпринять ?

[Nord]

1 . Можно попробовать отключить разностное формирование матрицы плотности
2. Можно попробовать использовать другое начальное приближение
3. Можно взять другой базис и/или увеличить точность вычисления интегралов в исходном.
4. Можно взять другой ускоритель сходимости.
5. Можно просто увеличить максимальное число итераций.
6. Можно считать в другом приближении.

Но всегда справедливо то, что проблемы сходимости означают проблемы с описанием физики задачи.

[Vladchimic]

Спасибо
1 . Можно попробовать отключить разностное формирование матрицы плотности
2. Можно попробовать использовать другое начальное приближение
3. Можно взять другой базис и/или увеличить точность вычисления интегралов в исходном.
4. Можно взять другой ускоритель сходимости.
5. Можно просто увеличить максимальное число итераций.
6. Можно считать в другом приближении.

1. Относитеьно первого - без понятия. Есть ли это в помощи к G98w и где ?
2. В качестве начального приближения использовал уже оптимизированную структуру в газовой фазе. Возможны иные, более лучшие варианты ?
3. Другой базис... Как всегда все упирается в мощности.. Это сейчас я считаю катион, а когда ионную пару ?
4. также к стыду своему не совсем понимаю.. можно поподроднее ?
5. Можно... попробую
6. да, тоже можно попробовать. спасибо еще раз.

[Nord]

По поводу конкретно G98 сказать ничего не могу - у меня его нету.
Ускорители (по англ. convergers) бывают разные. Один из наиболее популярных - DIIS

[Vladchimic]

DIIS используется по умолчанию......

[Raiden]

Попробуйте поменять начальное предположение. Т.е. к вашей строке добавить Guess=Hukkel либо Guess=Core. Подробнее смотрите в справке. 2003 гаусиан поддерживает команду Guess=AM1. Теоретически если сначала расчитать структуру в AM1, а затем написать Guess=Read, то может сработать и в 98.

[Vladchimic]

Спасибо. Попробую..

[duke2003]

попробуйте использовать другую модель для рассчета.
Я юзал так
#T opt rb3lyp/6-31g SCRF=(Dipole,a0=6.82,solvent=thf) SCF=Tight Test

a0 - вроде радиус сферы для модели. Его нужно тоже считать исходя из геометрии молекулы, реккомендуется делать это так

#T rb3lyp/6-31g volume geom=connectivity

Удачи.

[Vladchimic]

Использовал все четыре возможных модели... ни одна не считает.. видимо молекула очень велика....

[duke2003]

А действительно, попробуй сначало для более простой молекулки посчитать, если получиться, то переходи на более сложную.

[Vladchimic]

Уже попробовал... Обычные бензол, анилин, толуол считаются без проблем, причем с использованием любой из чеырех предложенных моделей расчета в Gaussian.... Видимо объем молекулы критичен..

[Rt19]

К стыду своему не знаю, что такое SCRF (пойду вечером документацию к Гауссиану читать). Но если это что-то вроде моделей поляризационного континуума, то тогда напрашивается увеличить размер полости.

[duke2003]

Верно говорят, попробуй увечить размер полости процентов на 10, может пройдет.
ЗЫ. У меня тоже было, что для больших молекул не удавалось оптимизирвать геометрию в растворе. Не судьба может..