Гауссиан: Spurious integrated density or basis function

[Darth Vasya]

Вообще в нормальных кв.-хим. лаб. на PC устанавливают 8Гб оперативной памяти.

В нормальных кв.-хим. лаб. считают не на PC

[PWRD]

Тогда для Me10 хороший большой потенциал, а для Sn4 подбирайте ЕСР.

пробывал штутгардские ECP использовать,
либо не сходится SCF, либо разносит молеклу, либо в кучу

[VTur]

Не заморачивайтесь. Значительная часть данных в публикациях того же качества.
У элементов 4 периода и дальше есть три проблемные области - около ядер, в третьем слое электронов и вдали от ядер. Хартри с Фоком их, по определению, не видят, потому и сходятся. Можете их и печатать.
Хотите корреляцию - базис должен быть сбалансирован, а значит для металлов очень большим.
Есть, правда, теория возмущений, но думаю, на РС ждать придется долго.

[Yurii]

У меня сложилось такое впечатление, что эту ошибку ребята из гауссиана специально не исправили в g03, хотя нареканий на нее хватает: оставили под g09 - там все ок при неизменной точности расчета. Да и под EM64T могли программу g09 значительно быстрей сделать, а не сделали (искуственно тормознули).

[o-oxhem]

Да и под EM64T могли программу g09 значительно быстрей сделать, а не сделали (искуственно тормознули).

Почему возникло такое мнение?

[Yurii]

Да и под EM64T могли программу g09 значительно быстрей сделать, а не сделали (искуственно тормознули).

Почему возникло такое мнение?

Для тех, кто занимается сборкой Гауссиана из исходников, это не является большим секретом. Из исходного кода можно выжать, смотря от задачи, еще до 50% скорости.

[Izberg]

Прошу помощи. Не сходится расчет. Пробовались DIIS, GDIIS, xqc, calcall. Болтается около одной геометрии. Хартри-фоком оптимизируется на ура.

[o-oxhem]

B3LYP - это функционал hybrid GGA approximation. А у Вас предполагается существенное разделение заряда, с которым GGA не может справиться. Попробуйте mp2 в том же базисе - всё должно "оптимизироваться на ура".

[Izberg]

GGA не может справиться. Попробуйте mp2 в том же базисе.

а кроме расчета мр2 еще что-то может помочь?

[o-oxhem]

Ну, попытайтесь немного изменить длины связей. Но я подчеркиваю, что Вы столкнулись с тем, что GGA почти не может описать. Впрочем, ещё можно попробовать функционалы с большей долей хатри-фоковского обмена, например, BHandLyp. Но ... может все-таки mp2?
Или уменьшите базис для b3lyp. Если вдруг сойдется, то стартуйте с опт. структуры в более широком базисе.
Можно попытаться с баз. наб. cc-pvdz или cc-pvtz.
Можно доюавить Maxstep=5 в группу opt.

[Yurii]

Прошу помощи. Не сходится расчет. Пробовались DIIS, GDIIS, xqc, calcall. Болтается около одной геометрии. Хартри-фоком оптимизируется на ура.

Используйте вместо b3lyp m062x и у вас тоже будет оптимизироваться на ура.

[Droog_Andrey]

Итак, Вашу проблему можно обойти.

Я добавил в скрипт: 6D 10F IOP(3/47=28)

Вот, теперь понятно. Иопы - хорошая штука.

Оказалось, дело в том, что в базисных наборах UGBS присутствуют базисные функции с очень узкими максимумами, для которых необходимо большее разрешение интегральной сетки в радиальном направлении, вплоть до Integral(DKHSO,Grid=600590) для тяжёлых атомов. Если же просто указать IOp(3/47=28), то в околоядерном пространстве получится чёрт знает что, а это существенно для того же золота и т.п.

[Yurii]

В g03 от этой ошибки для некоторых систем нельзя было избавиться и при "приличных" атомах (С, H), и при "приличных" базисах. Вроде я об этом уже писал.