GAMESS - как задать ECP для более одного тяжелого атома?

stat52349

Делаю по мануалу - не работает.
input_
$ECP
CS-ECP GEN 46 3
5 ----- F POTENTIAL -----
-0.78916699 2 0.93849999
-8.42115784 2 2.31629992
-30.98544312 2 6.00729990
-95.03477478 2 20.37969971
-30.07960320 1 59.32889938
7 ----- S-F POTENTIAL -----
42.85466766 2 1.38530004
-138.00901794 2 1.63240004
275.99960327 2 2.20580006
-280.45663452 2 3.22149992
199.82038879 2 4.64960003
27.73096657 1 15.15250015
3.76870608 0 19.00049973
7 ----- P-F POTENTIAL -----
48.66250992 2 1.25950003
-145.70526123 2 1.44169998
264.46368408 2 1.87639999
-279.85159302 2 2.65750003
184.35585022 2 3.63870001
23.30001831 1 10.65320015
5.76792908 0 14.68060017
7 ----- D-F POTENTIAL -----
34.86072540 2 0.76810002
-106.79302979 2 0.87459999
188.23532104 2 1.10769999
-217.63992310 2 1.48230004
137.74559021 2 1.92019999
34.42418671 1 6.21829987
7.19875193 0 17.38809967
7.41062880 0 9.28610039
CS-ECP
$END
output_
WRONG ECP LIST IN $ECP GROUP!
$ECP GROUP SHOULD DEFINE THE VALID ECP PARAMETERS FOR ALL ATOMS.

Komandor

Ты чё считаешь только один атом? Если нет, то надо писать в "$ECP" остовные потенциалы для всех атомов. Если на некоторых ты считаешь в полном базисе, то там пишешь "NONE". Даже, если атомы повторяются, всё-равно надо записывать "ecp" для каждого.

stat52349

Глянь внимательнее
Считаю молекулу Cs2. Для второго атома написано "CS-ECP", как по мануалу. Вылетает с вышеуказанной ошибкой.
Кстати, пробовал для обоих атомов ECP писать - не помогает.

Komandor

<Бл* чего париться, в G03 надо такие вещи считать, нефиг GAMESS юзать.
P/S: Ты какой базис для внешних оболочек берешь?

stat52349

Базис CRENBL ECP.

Komandor

СКОКА у тебя там базисных функций?

stat52349

Не знаю. Можешь сам глянуть.

Komandor

Если тебе очень важно посчитать эту молекулу, то после праздников можешь мне написать, я тебе помогу в любом базисе посчитать или остовном потенциале (причем с учетом релятивистских эффектов, это очень важно для цезия). А сейчас я уезжаю...

stat52349

Хорошо.

Soror

А где G03 можно скачать?

vovatroff

Обратите внимание на последнее поле в вашем примере:
7 ----- D-F POTENTIAL -----
34.86072540 2 0.76810002
-106.79302979 2 0.87459999
188.23532104 2 1.10769999
-217.63992310 2 1.48230004
137.74559021 2 1.92019999
34.42418671 1 6.21829987
7.19875193 0 17.38809967
7.41062880 0 9.28610039
Число 7 означает количество слагаемых в разложении для гармоник данного типа.
Оно должно совпадать с количеством нижеследующих числовых записей. У вас
же далее следует 8 записей, вот и вся ваша ошибка.

stat52349

Да. Последняя строчка оказалась из рубидия.
Оставить комментарий
Имя или ник:
Комментарий: