Учет растворителя в Gamess в рамках PCM(континуальная модель
Учет растворителя в Gamess в рамках PCM(континуальная модель
Хотелось узнать подробнее об этой процедуре и данной модели. Мне нужно учесть растворитель (этанол). Я использовал дополнительный тег
$PCM IEF=0 SOLVNT=C2H5OH EPC=27 $END
на сколько это уместно здесь или нужно все таки добавить и другие опции? Я вычисляю возбужденные состояния молекулы( вертикальные переходы).
$PCM IEF=0 SOLVNT=C2H5OH EPC=27 $END
на сколько это уместно здесь или нужно все таки добавить и другие опции? Я вычисляю возбужденные состояния молекулы( вертикальные переходы).
Re: Учет растворителя в Gamess в рамках PCM(континуальная мо
Обратите внимание на вот эту опцию в группе $TDDFT
NONEQ is a flag controlling PCM's solvent behavior: .TRUE. splits the dielectric constant into a bulk value (EPS in $PCM) and a fast component (EPSINF), see Cossi and Barone, 2001. The idea is that NONEQ=.t. is appropriate for vertical excitations, and .f. for adiabatic. (the default is .TRUE.)
В Вашем случае надо ставить .f.
Опция IEF=0 включает D-PCM. В мануале пишут:
D-PCM should be considered obsolete, while C-PCM is normally a better choice than IEF-PCM.
Но у меня с D-PCM результаты получались лучше, чем с C-PCM.
EPC=27 -- наверное, все-таки EPS=27? это лишнее. Если Вы задали SOLVNT из тех, чьи параметры уже зашиты в программу, то вводить параметры растворителя имеет смысл только тогда, когда Вас зашитые чем-то не устраивают. А если растворитель в Гамесс не зашит, но зашит в Гауссиан, то можно запустить маленький расчетик в Гауссиане (лучше G03, а не G09) и вытащить параметры растворителя из выдачи, они там печатаются.
Ссылки и описания вариантов PCM есть в мануале Гамесса (input.pdf и refs.pdf) и в мануале Гауссиана.
NONEQ is a flag controlling PCM's solvent behavior: .TRUE. splits the dielectric constant into a bulk value (EPS in $PCM) and a fast component (EPSINF), see Cossi and Barone, 2001. The idea is that NONEQ=.t. is appropriate for vertical excitations, and .f. for adiabatic. (the default is .TRUE.)
В Вашем случае надо ставить .f.
Опция IEF=0 включает D-PCM. В мануале пишут:
D-PCM should be considered obsolete, while C-PCM is normally a better choice than IEF-PCM.
Но у меня с D-PCM результаты получались лучше, чем с C-PCM.
EPC=27 -- наверное, все-таки EPS=27? это лишнее. Если Вы задали SOLVNT из тех, чьи параметры уже зашиты в программу, то вводить параметры растворителя имеет смысл только тогда, когда Вас зашитые чем-то не устраивают. А если растворитель в Гамесс не зашит, но зашит в Гауссиан, то можно запустить маленький расчетик в Гауссиане (лучше G03, а не G09) и вытащить параметры растворителя из выдачи, они там печатаются.
Ссылки и описания вариантов PCM есть в мануале Гамесса (input.pdf и refs.pdf) и в мануале Гауссиана.
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)
Re: Учет растворителя в Gamess в рамках PCM(континуальная мо
У меня возник такой вопрос - если растворитель не зашит в firefly, нужно ли указывать его название, к примеру:
$SOLVNT
IEF=3
SOLVNT=1,4-DIOXANE
EPS=2,209
$END
или можно
$SOLVNT
IEF=3
EPS=2,209
$END
$SOLVNT
IEF=3
SOLVNT=1,4-DIOXANE
EPS=2,209
$END
или можно
$SOLVNT
IEF=3
EPS=2,209
$END
Re: Учет растворителя в Gamess в рамках PCM(континуальная мо
Правильно так:
$PCM SOLVNT=input
RSOLV=2.155000 EPS=36.640000 EPSINF=1.806000
$END
$PCM SOLVNT=input
RSOLV=2.155000 EPS=36.640000 EPSINF=1.806000
$END
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)
Re: Учет растворителя в Gamess в рамках PCM(континуальная мо
спасибо. а откуда можно взять эти данные, к примеру для того же диоксана или эпихлоргидрина? в Гауссиане? и что значат RSOLV и EPSINF? EPS - это понятно, диэлектр. прониц.
Re: Учет растворителя в Gamess в рамках PCM(континуальная мо
Ну, если такой растворитель зашит в Гауссиан, то оттуда. Но ни диоксана ни, тем более, эпихлоргидрина, там нет Диоксан мне в свое время самой был нужен, но тогда знаний не хватило добыть нужные параметры.
RSOLV -- радиус растворителя. Видимо, придется взять из какой-нибудь примитивной молекулярной модели (нарисовали молекулу, грубо соптимизировали мол. механикой, взяли макс. диаметр, прибавили ван-дер-ваальсовские радиусы крайних атомов, поделили пополам).
EPSINF -- диэл. проницаемость на бесконечной частоте. Не знаю, откуда ее брать. Может, в справочниках найдется. Насколько я понимаю, EPSINF нужна только для расчетов с NONEQ=.t., а вертикальные спектры считаются с NONEQ=.f.
В общем, тут непросто все
RSOLV -- радиус растворителя. Видимо, придется взять из какой-нибудь примитивной молекулярной модели (нарисовали молекулу, грубо соптимизировали мол. механикой, взяли макс. диаметр, прибавили ван-дер-ваальсовские радиусы крайних атомов, поделили пополам).
EPSINF -- диэл. проницаемость на бесконечной частоте. Не знаю, откуда ее брать. Может, в справочниках найдется. Насколько я понимаю, EPSINF нужна только для расчетов с NONEQ=.t., а вертикальные спектры считаются с NONEQ=.f.
В общем, тут непросто все
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)
Re: Учет растворителя в Gamess в рамках PCM(континуальная мо
так тут наоборот написано что для вертикального возбуждения .t. а для адиабатического .f.sanya1024 писал(а):Обратите внимание на вот эту опцию в группе $TDDFT
NONEQ is a flag controlling PCM's solvent behavior: .TRUE. splits the dielectric constant into a bulk value (EPS in $PCM) and a fast component (EPSINF), see Cossi and Barone, 2001. The idea is that NONEQ=.t. is appropriate for vertical excitations, and .f. for adiabatic. (the default is .TRUE.)
В Вашем случае надо ставить .f.
Re: Учет растворителя в Gamess в рамках PCM(континуальная мо
Sorry, перепутала. Вы правы.
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 29 гостей