Проблема с псевдопотенциалом в Gamess

[AlexKsen]

Уважаемые форумчане, добрый вечер.
При TDDFT (def-2TZVP) в gamess выскочила ошибка связанная с псевдопотенциалом:

*** ERROR ***
AN UNEXPECTED $END WAS ENCOUNTERED IN $ECP GROUP.
EVERY ATOM IN THE MOLECULE MUST BE DEFINED IN $ECP.

Проверил, вроде, в синтаксе все правильно. Помогите разобраться в чем дело.

inp.txt

out.txt

[Tonnie Fox]

Не уверен, но вроде во входном файле блок DATA должен быть примерно такого вида для каждого атома:

[AlexKsen]

Если Вы про лишний Enter между симметрией и координатами, то он не лишний.

[Гесс]

У вас стоит ECP=READ но при этом нет блока $ECP - типа как в http://www.ccl.net/chemistry/resources/ ... 1.008-dir/

[AlexKsen]

Блок ЕСР прописан в inp

$ECP
C NONE
O NONE
...
H NONE
I-ECP GEN 28 3
4 ----- f-ul potential -----
-21.84204000 2 19.45860900
-28.46819100 2 19.34926000
...
0.24371300 2 4.82376700
0.32080400 2 4.88431500
H NONE
...
H NONE
$END

[Гесс]

Извините, незаметил. Но гамесс его кажется невидит.

Судя по примерам надо писать

O-ecp none

а не

O none

[Tonnie Fox]

Если Вы про лишний Enter между симметрией и координатами, то он не лишний.

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

[Гесс]

Стоп, а откуда симметрия Cs?
Там C1, а ECP наверняка ждет All atoms must be given here in $ECP, not just the symmetry unique atoms.

[Гесс]

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

TOTAL NUMBER OF BASIS SET SHELLS = 404
NUMBER OF CARTESIAN GAUSSIAN BASIS FUNCTIONS = 1250
NUMBER OF ELECTRONS = 262

орбитали вчитались. Невчитался именно ECP

[AlexKsen]

Вот оптимизация этой структуры с использованием еср с симметрией Сs. Прошу прощения, что файлы с облока:
out: https://cloud.mail.ru/public/Fjfs/MYvpdgBPQ

[Гесс]

Ну так там то структура Cs
А у вас С1

И покажите пожалуйста его инпут

[AlexKsen]

Ну так там то структура Cs
А у вас С1

И покажите пожалуйста его инпут

Кажется проглядел, когда inp делал. Проверю, спасибо.

[Tonnie Fox]

Кстати, а откуда нынче ECP брать для гамеса? Попробовал на basissetexchange, выбираю базис (тот же def-2TZVP), атомы, а в файле в итоге нет блока $ECP

[AlexKsen]

Кстати, а откуда нынче ECP брать для гамеса? Попробовал на basissetexchange, выбираю базис (тот же def-2TZVP), атомы, а в файле в итоге нет блока $ECP

Странно. Брал там же.

[Гесс]

Кстати, а откуда нынче ECP брать для гамеса? Попробовал на basissetexchange, выбираю базис (тот же def-2TZVP), атомы, а в файле в итоге нет блока $ECP

Вроде все норм

ECP.jpg

Или вы про элементы до криптона?

[Tonnie Fox]

Спасибо, уже нашел проблему в своем MS Edge, тупенький он или не обновлялся давно ))
Проверил Ваш инпут на домашнем калькуляторе, правда гамесса у меня нет, запустил в FF и с ключом dfttyp=b3lyp. Действительно на симметрию ругается прога, при том что задано всё вроде верно. Мне не доводилось проводить расчеты в условиях гамесс/ECP/группа симметрия, теперь самому интересно, проблема ли это или есть какие то ключи решения без перехода к С1?

[AlexKsen]

Спасибо, уже нашел проблему в своем MS Edge, тупенький он или не обновлялся давно ))
Проверил Ваш инпут на домашнем калькуляторе, правда гамесса у меня нет, запустил в FF и с ключом dfttyp=b3lyp. Действительно на симметрию ругается прога, при том что задано всё вроде верно. Мне не доводилось проводить расчеты в условиях гамесс/ECP/группа симметрия, теперь самому интересно, проблема ли это или есть какие то ключи решения без перехода к С1?

Буду проверять.

[Гесс]

Мануал говорит лишь что нужно указывать ECP для ВСЕХ атомов а не только для symmetry unique atoms.
И в облаке был аутпут с ECP и Cs. Думаю надо лишь посмотреть на его инпут.

[AlexKsen]

Мануал говорит лишь что нужно указывать ECP для ВСЕХ атомов а не только для symmetry unique atoms.
И в облаке был аутпут с ECP и Cs. Думаю надо лишь посмотреть на его аутпут.

Какой нужен out?

[Гесс]

опечатка. инпут. Чтото я сегодня несконцентрирован.