Самосогласование HF Гауссиан

вопросы строения молекул и квантовой химии
Ответить
annac
Сообщения: 3
Зарегистрирован: Пн фев 20, 2017 1:53 pm

Самосогласование HF Гауссиан

Сообщение annac » Пн фев 20, 2017 2:08 pm

Добрый день.
Провожу расчет двухатомной молекулы в Гауссиане, при расстояниях более 2.85 ангстрем не согласуется Хартри Фок, а при написании SCF=QC пишет, что слишком большой градиент для этого метода (Gradient too large for Newton-Raphson -- use scaled steepest descent instead.). Подскажите возможные решения проблемы.

Аватара пользователя
amge
Сообщения: 2021
Зарегистрирован: Вт июл 31, 2007 11:42 am

Re: Самосогласование HF Гауссиан

Сообщение amge » Пн фев 20, 2017 3:00 pm

Гауссиан после этого вылетает? Если нет, то попробуйте подождать: в сети масса примеров, когда после неоднократных таких предупреждений все кончается благополучно.

Аватара пользователя
Droog_Andrey
Сообщения: 2670
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: Самосогласование HF Гауссиан

Сообщение Droog_Andrey » Пн фев 20, 2017 3:05 pm

Добавьте диффузных функций. Ну и при той же мультиплетности состояние может стать возбуждённым при сильном возрастании межъядерного расстояния.
2^74207281-1 is prime!

Аватара пользователя
Гесс
Сообщения: 13063
Зарегистрирован: Ср фев 15, 2012 11:19 pm

Re: Самосогласование HF Гауссиан

Сообщение Гесс » Пн фев 20, 2017 8:15 pm

И если это заметно длиннее "традиционной" связи то там может просто начинаться мультиреференс от которого все умирает. N2 тому хорошо обкатанный грановским пример. Тут уже хоть кьюси, хоть некьюси, пофиг.

Аватара пользователя
Vinianel
Сообщения: 49
Зарегистрирован: Пт ноя 25, 2011 12:32 pm

Re: Самосогласование HF Гауссиан

Сообщение Vinianel » Вт мар 07, 2017 3:14 pm

N2 тому хорошо обкатанный грановским пример
А можно ссылочку, где об этом почитать подробнее? :shuffle: "N2 Granovsky multireference" не гуглится.

Аватара пользователя
Гесс
Сообщения: 13063
Зарегистрирован: Ср фев 15, 2012 11:19 pm

Re: Самосогласование HF Гауссиан

Сообщение Гесс » Вт мар 07, 2017 6:26 pm

где то около 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

Аватара пользователя
Гесс
Сообщения: 13063
Зарегистрирован: Ср фев 15, 2012 11:19 pm

Re: Самосогласование HF Гауссиан

Сообщение Гесс » Пт май 12, 2017 1:39 am

диссоциация азота похоже становится "морской свинкой" http://aip.scitation.org/doi/full/10.1063/1.4944087

Аватара пользователя
Droog_Andrey
Сообщения: 2670
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: Самосогласование HF Гауссиан

Сообщение Droog_Andrey » Сб май 13, 2017 9:03 pm

Ещё интереснее диссоциация квадратного синглетного C42+ на два квадруплетных C2+.
2^74207281-1 is prime!

Аватара пользователя
Гесс
Сообщения: 13063
Зарегистрирован: Ср фев 15, 2012 11:19 pm

Re: Самосогласование HF Гауссиан

Сообщение Гесс » Сб май 13, 2017 10:51 pm

Droog_Andrey писал(а):
Сб май 13, 2017 9:03 pm
Ещё интереснее диссоциация квадратного синглетного C42+ на два квадруплетных C2+.
Нууу... оно может и интереснее, но частица больно иллюзорная. А азота - хоть залейся, какой-то умный эксперимент сделать можно.

Аватара пользователя
madschumacher
Сообщения: 883
Зарегистрирован: Ср авг 05, 2015 4:30 pm

Re: Самосогласование HF Гауссиан

Сообщение madschumacher » Сб май 13, 2017 11:40 pm

Droog_Andrey писал(а):
Сб май 13, 2017 9:03 pm
Ещё интереснее диссоциация квадратного синглетного C42+ на два квадруплетных C2+.
А чем квадруплет то интересен? :?
И да узрел Охламон, что сие есть круть несусветная!

Аватара пользователя
Droog_Andrey
Сообщения: 2670
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: Самосогласование HF Гауссиан

Сообщение Droog_Andrey » Ср май 17, 2017 12:00 pm

Сам квадруплет - ничем, это просто основное состояние для катиона C2+. Интересна диссоциация :)
2^74207281-1 is prime!

Аватара пользователя
madschumacher
Сообщения: 883
Зарегистрирован: Ср авг 05, 2015 4:30 pm

Re: Самосогласование HF Гауссиан

Сообщение madschumacher » Ср май 17, 2017 12:07 pm

Droog_Andrey писал(а):
Ср май 17, 2017 12:00 pm
Интересна диссоциация
Ну с многоатомными частицами всегда не очень хорошо изучать диссоциацию, хотя бы из-за того, что ППЭ сразу становится многомерной (даже из квадрата C4, при условии, что выбраны и зафиксированы связи C--C в двух "осколках", и мы не учитываем координату "линеаризации", есть координаты "растяжения", поворота связей C--C друг относительно друга по оси проходящей через центр). Короче, там нельзя уже так посмотреть на кривульку и что-то сказать... :( Многомерные (т.е. имеющие > 2-х координат) = зло! :evil:
И да узрел Охламон, что сие есть круть несусветная!

Аватара пользователя
Droog_Andrey
Сообщения: 2670
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: Самосогласование HF Гауссиан

Сообщение Droog_Andrey » Ср май 17, 2017 11:05 pm

Но путь через минимальный барьер всё равно один. Кстати, эта диссоциация эндотермична, чего DFT в упор не видит :)
2^74207281-1 is prime!

Аватара пользователя
madschumacher
Сообщения: 883
Зарегистрирован: Ср авг 05, 2015 4:30 pm

Re: Самосогласование HF Гауссиан

Сообщение madschumacher » Чт май 18, 2017 12:36 am

Droog_Andrey писал(а):
Ср май 17, 2017 11:05 pm
чего DFT в упор не видит :)
Ой, где в DFT начинаются игры со спиновыми состояниями, там начинается Мория неизведанных открытий и НЁХ... :lol:
И да узрел Охламон, что сие есть круть несусветная!

Ответить

Вернуться в «квантовая химия и моделирование»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 53 гостя