Excessive mixing of frozen core and valence orbitals...

вопросы строения молекул и квантовой химии
Ответить
Аватара пользователя
surius
Сообщения: 1489
Зарегистрирован: Пт сен 21, 2007 11:20 am

Excessive mixing of frozen core and valence orbitals...

Сообщение surius » Пт окт 30, 2015 2:38 pm

Господа, знатоки и им сочувствующие, очередной вопрос от плохо образованной робочей молодежи:

Есть некий стандартный базисный набор (тот который зашит и может быть вызван ключевым словом), пусть как и в предыдущем вопросе def2TZVPP, который для некоторых элементов содержащий ECP (I например). И вот мы такие делаем расчеты, ну например оптимизируем геометрию в основном состоянии.Спрашивается, нужно ли ручками G09 говорить (gen pseudo=read) чтоб он для атома йода вытащил ECP и применил его.
Потому как если нужно, то это адский криминал, а если нет, не нужно (то есть когда код видит что там I и для него есть вшитый ECP, и он его использует), то я решительно ничего не понимаю :issue:

UPD:
Все работает в авто режиме

И в связи с этим продолжение вопроса: в основном состоянии все прекрасно работает. А вот когда начинается vertical excitation расчет в некоторых, избранных случаях расчет вываливается с ошибкой:

Код: Выделить всё

 
Largest valence mixing into a core orbital is  3.44D-01
Largest core mixing into a valence orbital is  4.69D-01
Excessive mixing of frozen core and valence orbitals.
Error termination via Lnk1e in /***/***/****/***/gaussian_intel//g09/l801.exe at Fri Oct 30 12:47:39 2015.
Гугл показывает что такая проблема случается с post-SCF расчетами, и зачастую связана с
This typically means you are using a custom effective core potential and you forgot to use genECP or gen pseudo = read as basis set specification. If you just use Gen it still runs, I'm unsure whether it uses a default ECP or just tries to make do with the specified basis set, either way it ends bad.
И действительно, в одном и моих случаев дело могло быть как-то связано с псевдопотенциалом, потому как были таковые атомы (I, например), но в других случаях не было атомов с псевдопотенциалами...

PS
Гесс, если вы решите что тема недостойна быть самостоятельной, то может пришить ее к def2***
Последний раз редактировалось surius Пт окт 30, 2015 5:06 pm, всего редактировалось 4 раза.
"Bite my shiny metal ass"
Bender

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

Re: ECP в базисном наборе

Сообщение Гесс » Пт окт 30, 2015 2:57 pm

surius писал(а):Спрашивается, нужно ли ручками G09 говорить (gen pseudo=read) чтоб он для атома йода вытащил ECP и применил его.
Pseudo without any options defaults to Pseudo=Read.
Gen имхо тоже ненужен. А вот сказать гауссиану что в базисе есть псевдопотенциал - имхо невредно. Поэтому пустое Pseudo я бы поставил. Орка кстати не требует. Собственно - а почему не попробовать загнать HI в def2-SVP с и без Pseudo и посмотреть на разницу?

Аватара пользователя
surius
Сообщения: 1489
Зарегистрирован: Пт сен 21, 2007 11:20 am

Re: ECP в базисном наборе

Сообщение surius » Пт окт 30, 2015 3:45 pm

Гесс писал(а):а почему не попробовать загнать HI в def2-SVP с и без Pseudo и посмотреть на разницу?
светлая мысль. Сделано. Все работает в автомате, что не удивительно. колличество базисных функции, примитивных гауссиан и т.п. в точности одно и тожею Ну и конечно энергии в обоих случаях эквивалентны.
"Bite my shiny metal ass"
Bender

Аватара пользователя
surius
Сообщения: 1489
Зарегистрирован: Пт сен 21, 2007 11:20 am

Re: ECP в базисном наборе

Сообщение surius » Пт окт 30, 2015 4:37 pm

Товарищи, а кто-нибудь вкурсе как в гауссиане рулить frozen core для DFT?
"Bite my shiny metal ass"
Bender

Ответить

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

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

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