Есть некий стандартный базисный набор (тот который зашит и может быть вызван ключевым словом), пусть как и в предыдущем вопросе def2TZVPP, который для некоторых элементов содержащий ECP (I например). И вот мы такие делаем расчеты, ну например оптимизируем геометрию в основном состоянии.Спрашивается, нужно ли ручками G09 говорить (gen pseudo=read) чтоб он для атома йода вытащил ECP и применил его.
Потому как если нужно, то это адский криминал, а если нет, не нужно (то есть когда код видит что там I и для него есть вшитый ECP, и он его использует), то я решительно ничего не понимаю
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.
И действительно, в одном и моих случаев дело могло быть как-то связано с псевдопотенциалом, потому как были таковые атомы (I, например), но в других случаях не было атомов с псевдопотенциалами...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.
PS
Гесс, если вы решите что тема недостойна быть самостоятельной, то может пришить ее к def2***