Самосогласование HF Гауссиан
Самосогласование HF Гауссиан
Добрый день.
Провожу расчет двухатомной молекулы в Гауссиане, при расстояниях более 2.85 ангстрем не согласуется Хартри Фок, а при написании SCF=QC пишет, что слишком большой градиент для этого метода (Gradient too large for Newton-Raphson -- use scaled steepest descent instead.). Подскажите возможные решения проблемы.
Провожу расчет двухатомной молекулы в Гауссиане, при расстояниях более 2.85 ангстрем не согласуется Хартри Фок, а при написании SCF=QC пишет, что слишком большой градиент для этого метода (Gradient too large for Newton-Raphson -- use scaled steepest descent instead.). Подскажите возможные решения проблемы.
Re: Самосогласование HF Гауссиан
Гауссиан после этого вылетает? Если нет, то попробуйте подождать: в сети масса примеров, когда после неоднократных таких предупреждений все кончается благополучно.
- Droog_Andrey
- Сообщения: 2670
- Зарегистрирован: Сб сен 29, 2007 8:29 pm
- Контактная информация:
Re: Самосогласование HF Гауссиан
Добавьте диффузных функций. Ну и при той же мультиплетности состояние может стать возбуждённым при сильном возрастании межъядерного расстояния.
2^74207281-1 is prime!
Re: Самосогласование HF Гауссиан
И если это заметно длиннее "традиционной" связи то там может просто начинаться мультиреференс от которого все умирает. N2 тому хорошо обкатанный грановским пример. Тут уже хоть кьюси, хоть некьюси, пофиг.
Re: Самосогласование HF Гауссиан
А можно ссылочку, где об этом почитать подробнее? "N2 Granovsky multireference" не гуглится.N2 тому хорошо обкатанный грановским пример
Re: Самосогласование HF Гауссиан
где то около 24 и далее страниц тут http://www.msg.ameslab.gov/tutorials/mcscf1.pdf и по мелочи чтото есть во второй лекции (заменить 1 на 2 в имени файла)
Upd from 2022: new links https://www.msg.chem.iastate.edu/tutorials/mcscf1.pdf https://www.msg.chem.iastate.edu/tutorials/mcscf2.pdf
Upd from 2022: new links https://www.msg.chem.iastate.edu/tutorials/mcscf1.pdf https://www.msg.chem.iastate.edu/tutorials/mcscf2.pdf
Re: Самосогласование HF Гауссиан
диссоциация азота похоже становится "морской свинкой" http://aip.scitation.org/doi/full/10.1063/1.4944087
- Droog_Andrey
- Сообщения: 2670
- Зарегистрирован: Сб сен 29, 2007 8:29 pm
- Контактная информация:
Re: Самосогласование HF Гауссиан
Ещё интереснее диссоциация квадратного синглетного C42+ на два квадруплетных C2+.
2^74207281-1 is prime!
Re: Самосогласование HF Гауссиан
Нууу... оно может и интереснее, но частица больно иллюзорная. А азота - хоть залейся, какой-то умный эксперимент сделать можно.Droog_Andrey писал(а): ↑Сб май 13, 2017 9:03 pmЕщё интереснее диссоциация квадратного синглетного C42+ на два квадруплетных C2+.
- madschumacher
- Сообщения: 883
- Зарегистрирован: Ср авг 05, 2015 4:30 pm
Re: Самосогласование HF Гауссиан
А чем квадруплет то интересен?Droog_Andrey писал(а): ↑Сб май 13, 2017 9:03 pmЕщё интереснее диссоциация квадратного синглетного C42+ на два квадруплетных C2+.
И да узрел Охламон, что сие есть круть несусветная!
- Droog_Andrey
- Сообщения: 2670
- Зарегистрирован: Сб сен 29, 2007 8:29 pm
- Контактная информация:
Re: Самосогласование HF Гауссиан
Сам квадруплет - ничем, это просто основное состояние для катиона C2+. Интересна диссоциация
2^74207281-1 is prime!
- madschumacher
- Сообщения: 883
- Зарегистрирован: Ср авг 05, 2015 4:30 pm
Re: Самосогласование HF Гауссиан
Ну с многоатомными частицами всегда не очень хорошо изучать диссоциацию, хотя бы из-за того, что ППЭ сразу становится многомерной (даже из квадрата C4, при условии, что выбраны и зафиксированы связи C--C в двух "осколках", и мы не учитываем координату "линеаризации", есть координаты "растяжения", поворота связей C--C друг относительно друга по оси проходящей через центр). Короче, там нельзя уже так посмотреть на кривульку и что-то сказать... Многомерные (т.е. имеющие > 2-х координат) = зло!
И да узрел Охламон, что сие есть круть несусветная!
- Droog_Andrey
- Сообщения: 2670
- Зарегистрирован: Сб сен 29, 2007 8:29 pm
- Контактная информация:
Re: Самосогласование HF Гауссиан
Но путь через минимальный барьер всё равно один. Кстати, эта диссоциация эндотермична, чего DFT в упор не видит
2^74207281-1 is prime!
- madschumacher
- Сообщения: 883
- Зарегистрирован: Ср авг 05, 2015 4:30 pm
Re: Самосогласование HF Гауссиан
Ой, где в DFT начинаются игры со спиновыми состояниями, там начинается Мория неизведанных открытий и НЁХ...
И да узрел Охламон, что сие есть круть несусветная!
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 53 гостя