Расчеты в gaussian с учетом растворителя... Проблемы.....
-
- Сообщения: 97
- Зарегистрирован: Чт июл 08, 2004 6:39 pm
- Контактная информация:
Расчеты в gaussian с учетом растворителя... Проблемы.....
Входной файл
%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.
Может ли кто-нибудь посказать мне, в чем дело и что может мне помочь ? Спасибо....
%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.
Может ли кто-нибудь посказать мне, в чем дело и что может мне помочь ? Спасибо....
-
- Сообщения: 97
- Зарегистрирован: Чт июл 08, 2004 6:39 pm
- Контактная информация:
1 . Можно попробовать отключить разностное формирование матрицы плотности
2. Можно попробовать использовать другое начальное приближение
3. Можно взять другой базис и/или увеличить точность вычисления интегралов в исходном.
4. Можно взять другой ускоритель сходимости.
5. Можно просто увеличить максимальное число итераций.
6. Можно считать в другом приближении.
Но всегда справедливо то, что проблемы сходимости означают проблемы с описанием физики задачи.
2. Можно попробовать использовать другое начальное приближение
3. Можно взять другой базис и/или увеличить точность вычисления интегралов в исходном.
4. Можно взять другой ускоритель сходимости.
5. Можно просто увеличить максимальное число итераций.
6. Можно считать в другом приближении.
Но всегда справедливо то, что проблемы сходимости означают проблемы с описанием физики задачи.
Не важно, что о вас говорят современники, важно что о вас скажут потомки
-
- Сообщения: 97
- Зарегистрирован: Чт июл 08, 2004 6:39 pm
- Контактная информация:
Спасибо
1 . Можно попробовать отключить разностное формирование матрицы плотности
2. Можно попробовать использовать другое начальное приближение
3. Можно взять другой базис и/или увеличить точность вычисления интегралов в исходном.
4. Можно взять другой ускоритель сходимости.
5. Можно просто увеличить максимальное число итераций.
6. Можно считать в другом приближении.
1. Относитеьно первого - без понятия. Есть ли это в помощи к G98w и где ?
2. В качестве начального приближения использовал уже оптимизированную структуру в газовой фазе. Возможны иные, более лучшие варианты ?
3. Другой базис... Как всегда все упирается в мощности.. Это сейчас я считаю катион, а когда ионную пару ?
4. также к стыду своему не совсем понимаю.. можно поподроднее ?
5. Можно... попробую
6. да, тоже можно попробовать. спасибо еще раз.
1 . Можно попробовать отключить разностное формирование матрицы плотности
2. Можно попробовать использовать другое начальное приближение
3. Можно взять другой базис и/или увеличить точность вычисления интегралов в исходном.
4. Можно взять другой ускоритель сходимости.
5. Можно просто увеличить максимальное число итераций.
6. Можно считать в другом приближении.
1. Относитеьно первого - без понятия. Есть ли это в помощи к G98w и где ?
2. В качестве начального приближения использовал уже оптимизированную структуру в газовой фазе. Возможны иные, более лучшие варианты ?
3. Другой базис... Как всегда все упирается в мощности.. Это сейчас я считаю катион, а когда ионную пару ?
4. также к стыду своему не совсем понимаю.. можно поподроднее ?
5. Можно... попробую
6. да, тоже можно попробовать. спасибо еще раз.
-
- Сообщения: 97
- Зарегистрирован: Чт июл 08, 2004 6:39 pm
- Контактная информация:
Попробуйте поменять начальное предположение. Т.е. к вашей строке добавить Guess=Hukkel либо Guess=Core. Подробнее смотрите в справке. 2003 гаусиан поддерживает команду Guess=AM1. Теоретически если сначала расчитать структуру в AM1, а затем написать Guess=Read, то может сработать и в 98.
Кошка, которая гуляет сама по себе - это кошка Мебиуса.
-
- Сообщения: 97
- Зарегистрирован: Чт июл 08, 2004 6:39 pm
- Контактная информация:
-
- Сообщения: 97
- Зарегистрирован: Чт июл 08, 2004 6:39 pm
- Контактная информация:
-
- Сообщения: 97
- Зарегистрирован: Чт июл 08, 2004 6:39 pm
- Контактная информация:
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 3 гостя