Gamess-Us. Использование псевдопотенциалов и базисов для Lu.

[etzell]

Доброе время. на сайте http://www.theochem.uni-stuttgart.de/ps ... se.en.html можно выбрать псевдопотенциал и базис для лантанидов. Gamess - формат там недоступен, но есть формат для Gaussian, который можно переделать в Gamess, кам не сказали. Какие ключи нужно использовать в инпут файле и как вообще подключать внешний базис (или его можно сразу вставить в input-файл Gamess) для рачсета с ECP (молекула содержит лютеций)?

[alxyppv]

это опять я. в Орке вшиты базисы из турбрмола. Если вы делаете расчет в Орке с исползованием опции PrintBasis, то в выдаче будет напечатан базис в фомате гамесс. Можно просто копировать
А так - мне, например, неизвестны опции автоматического пернеделывания форматов, но здесь и вручную один раз можно переделать просто сравнивая форматы. Или посмотрет на сайте https://bse.pnl.gov/bse/portal - там есть и гамессовский формат.
Пр подключении внешнего базиса я обычно вставляю его непосредственно в группу $DATA, но тогда нужно ставить базис после каждого атома, а группу $BASIS убить.

[etzell]

Хорошо, попробую. Спасибо!

[etzell]

Кстати а как быть с Лютецием? Дл лютеция Gamess-Us- формат выдает на сайте следующее:
$ECP
LU-ECP GEN 28 5
1 ----- h-ul potential -----
0.000000000 2 1.000000000
1 ----- s-ul potential -----
989.995584000 2 35.162097000
1 ----- p-ul potential -----
278.865652000 2 19.464402000
1 ----- d-ul potential -----
71.009178000 2 10.006865000
1 ----- f-ul potential -----
-47.405890000 2 23.517932000
1 ----- g-ul potential -----
-35.557146000 2 29.412238000
$END
------------------------
Если я буду вставлять это все в $DATA, то у меня выходит ошибка, так как у меня ключ $ECP открывается. Это нужно вставлять после $DATA, тогда как быть с атомом лютецием, для которого задаются координаты?

[sanya1024]

check.zip

Здесь в файле используется внешний базис,

BASIS.ZIP

но можно вставлять базисы прямо в $DATA, как обычно.

[VTur]

Эти базисы
а) "CRENBL ECP"
La: R.B. ROSS, W.C. ERMLER, P.A. CHRISTIANSEN ET AL. J. CHEM.
! PHYS. 93,(1990)6654.
!Ce - Lu: R.B. Ross, W.C. Ermler, S. Das, To be published
!Hf - Rn: R.B. ROSS, W.C. ERMLER, P.A. CHRISTIANSEN ET AL. J. CHEM.
! PHYS. 93,(1990)6654.
!Fr - Pu: W.C. ERMLER, R.B. ROSS, P.A. CHRISTIANSEN, INT. J. QUANT.
! CHEM. 40,(1991)829.
!Am - Uuo:C.S. NASH, B.E. BURSTEN, W.C. ERMLER, J. CHEM. PHYS. 1997

б) "Stuttgart RSC 1997 ECP"
Rf - Ha: M. Dolg, H. Stoll, H. Preuss, R.M. Pitzer, J. Phys. Chem. 97, 5852
! (1993).
есть в ЧемКрафте. Вы их можете автоматически использовать для генерации входного файла ГАМЕСС,
если недостаточно, подгрузите свой и автоматом создавайте свой файл




.

[etzell]

Спасибо. не знал про эту возможность Чемкрафта) Надо будет посмотреть еще и другие возможности этой программы.
Кстати, когда я выбираю базисы (Штуттгарт), то в $DATA нет базисных функции для атомов водорода. Для них в ручную нужно выбирать базисные функции из другого базиса? Или как в этом случае быть?

[VTur]

Либо прописать для него свой базис руками,
либо подгрузить другой базис с водородами,
либо взять понравившийся базис, поставить водороды первыми, создать новый файл, вырезать из него водороды с базисом, вставить в свой.

[sanya1024]

Я посмотрела, как в ChemCraft это работает... что-то странновато получается. На мой взгляд, руками проще Базисы тут: https://bse.pnl.gov/bse/portal. Можно создать в ChemCraft некую "рыбу", а потом доработать ее ручками.

[VTur]

Если на портале нужный базис есть, то ничего дорабатывать не надо.
etzell, вышлите мне вашу структуру и укажите базис на портале.

[sanya1024]

Вот именно: взять базис с портала и если где чего не вставилось ChemCraft-ом, вставить ручками.
Но по-моему гораздо удобнее создать внешний файл с базисом (взять нужные базисы для элементов с портала) и подгружать его. Тогда, например, при рестарте оптимизации геометрии не придется заново вписывать базисы в инпут (ручками или ChemCraft-ом), побеспокоиться придется только о том, чтобы не потерялась группа $ECP.