Расчет DFT [Ti(H2O)6]3+

[StudentQC]

Здравствуйте, я хотел бы узнать правильно ли задан инпут для рассчета оптимизации геометрии и энергии комплекса [Ti(H2O)6]3+
( Насколько данный расчет верен/правдоподобен? ).
( До этого методом DFT считал простые небольшие молекулы).
Рассчитываю в GAMESS (US) , открываю ChemCraft .
Инпут:

$CONTRL UNITS=ANGS COORD=CART RUNTYP=OPTIMIZE SCFTYP=ROHF
ICHARG= 3 MULT= 2 maxit=200 $END
$GUESS GUESS= HUCKEL $END
$DFT DFTTYP=B3LYP $END
$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 $END
$SYSTEM MWORDS=2 MEMDDI=2 $END
$Statpt NSTEP=10000 $END

$DATA
Untitled-2
C1
Ti 22.0 0.061480 -0.001395 0.005673
O 8.0 -1.324443 1.384529 0.005673
O 8.0 0.754442 0.691567 1.703075
O 8.0 1.261725 1.198850 -0.974323
O 8.0 -1.138765 -1.201639 0.985669
O 8.0 -0.619445 -0.682319 -1.701444
O 8.0 1.447403 -1.387318 0.005673
H 1.0 -0.897380 2.279884 0.005673
H 1.0 -1.889898 1.287366 0.814908
H 1.0 1.036818 1.635151 1.584939
H 1.0 1.532018 0.146726 1.990421
H 1.0 1.420210 2.018542 -0.438572
H 1.0 0.848088 1.445010 -1.841705
H 1.0 -0.680333 -2.065624 1.151235
H 1.0 -1.958959 -1.359214 0.450417
H 1.0 -1.563834 -0.965499 -1.591982
H 1.0 -0.072568 -1.457858 -1.990421
H 1.0 1.023129 -2.279884 -0.080124
H 1.0 1.958959 -1.344055 0.854485
$END

[Fanfate]

$SCF DIRSCF=.TRUE. $END
Думаю эту карту следует добавить в любом случае, иначе расчет становится более долгим (значительно), а ещё сбережет ваш диск, если вы считаете дома.

[StudentQC]

Спасибо

[amge]

Еще для ускорения оптимизации геометрии добавить
$control nzvar=51 $end
$zmat dlc=.t. auto=.t. $end

[himdim2012]

( Насколько данный расчет верен/правдоподобен? ).

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

[StudentQC]

Еще для ускорения оптимизации геометрии добавить
$control nzvar=51 $end
$zmat dlc=.t. auto=.t. $end

Ускорение расчетов очень пригодится, спасибо.

[Shorku]

( Насколько данный расчет верен/правдоподобен? ).

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

Прямой учет растворителя, тем более целой сольватной оболочки - это не так легко и просто, как Вам могло бы показаться.

[ Post made via Android ]

[himdim2012]

Прямой учет растворителя, тем более целой сольватной оболочки - это не так легко и просто, как Вам могло бы показаться.

естественно это многократно усложняет расчетную задачу, но зато приближает её к реальным объектам. Вопрос ТС в том числе касался правдоподобности данной модели. Вот я и высказал своё мнение.

[sanya1024]

Подумаешь, бином Ньютона всего лишь достроить сольватную оболочку с помощью EFP При расчете это почти ничего не будет стоить.

Да, неплохо было бы на водороды добавить поляризационных функций NPFUNC=1.
И для Ti я бы лучше взяла, например, LANL с квазирелятивистским псевдопотенциалом -- как-никак, переходный металл.

[StudentQC]

Подумаешь, бином Ньютона всего лишь достроить сольватную оболочку с помощью EFP При расчете это почти ничего не будет стоить.

Да, неплохо было бы на водороды добавить поляризационных функций NPFUNC=1.
И для Ti я бы лучше взяла, например, LANL с квазирелятивистским псевдопотенциалом -- как-никак, переходный металл.

Я бы с удовольствие осуществил расчет с использованием EFP и LANL ( читал о них в мануале GAMESS ).
Но не представляю как это должно выглядеть в инпут файле. Есть ли где-нибудь образцы таких инпутов?
Да и как навешать поляризационные функции на водороды ?

[comcon1]

на водороды:
$basis npfunc=1 $end

LANL ECP скачать на EMSL. можно выбрать LANL2 и базис DZVP например.

[sanya1024]

По EFP -- погодите, разгребусь немного и напишу главу методички (в ЖЖ).