Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для чайни

вопросы строения молекул и квантовой химии
Ответить
lmihael
Сообщения: 227
Зарегистрирован: Вт ноя 21, 2006 3:33 pm

Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для чайни

Сообщение lmihael » Пт май 07, 2010 1:38 pm

Здравствуйте. Надеюсь не очень надоел вам своими вопросами про возбужденное состояние. Хочу попробовать что-то более серьезное чем TDDFT/b3lyp и CIS.
Слышал про CASSCF/MCQDPT, caspt2 и eom-ccsd.Что это такое, чем они отличаются и не одно ли это и тоже :)? А также что из них реализовано в Firefly? Каким методом лучше оптимизировать геометрию, а каким определять энергию и дипольный момент первого возбужденого состояния? Какие базисы следует использовать? Был бы очень благодарен за алгоритм расчета на примере молекулы этаналя.

VTur
Сообщения: 7357
Зарегистрирован: Пт авг 31, 2007 1:36 pm

Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч

Сообщение VTur » Пт май 07, 2010 2:06 pm

В малых базисах эти методы работают плохо, в больших базисах не пробовал.
После отстоя требуйте долива

lmihael
Сообщения: 227
Зарегистрирован: Вт ноя 21, 2006 3:33 pm

Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч

Сообщение lmihael » Пт май 07, 2010 2:53 pm

VTur писал(а):В малых базисах эти методы работают плохо, в больших базисах не пробовал.
А в малых это каких? И можете ли вы посоветовать методы, которые работают хорошо?

Аватара пользователя
Lantano
Сообщения: 189
Зарегистрирован: Чт окт 09, 2008 10:35 pm

Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч

Сообщение Lantano » Пт май 07, 2010 9:13 pm

В Firefly есть CASSCF с возможностью расчета аналитических градиентов. Т.е. можно "естественным образом" оптимизировать геометрические параметры. Есть теория возмущений на его основе MCQDPT, а также какие-то варианты связанных кластеров. Но в них градиенты аналитически не считаются.
Какое у Вас первое возбужденное состояние? Для ацетальдегида по энергии это триплет n-->pi*. В этом случае выбор методов расчета побогаче.
P.S. Можете написать мне на e-mail. При более предметном разговоре может быть что-нибудь дельное посоветую.
"Если миру не нужен свет - Не беда, проживёт и так.
Если стал не нужен поэт - Не беда, он наймётся в кабак..."

lmihael
Сообщения: 227
Зарегистрирован: Вт ноя 21, 2006 3:33 pm

Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч

Сообщение lmihael » Пн май 10, 2010 9:13 am

Lantano писал(а):В Firefly есть CASSCF с возможностью расчета аналитических градиентов. Т.е. можно "естественным образом" оптимизировать геометрические параметры. Есть теория возмущений на его основе MCQDPT, а также какие-то варианты связанных кластеров. Но в них градиенты аналитически не считаются.
Какое у Вас первое возбужденное состояние? Для ацетальдегида по энергии это триплет n-->pi*. В этом случае выбор методов расчета побогаче.
P.S. Можете написать мне на e-mail. При более предметном разговоре может быть что-нибудь дельное посоветую.
Очень бы хотелось на примере ацетальдегида или любой другой несимметричной молекулы оптимизировать геометрию первого возбужденного синглетного состояния (S1) методом CASSCF, а затем рассчитать дипольный момент и полную энергию этого состояния в MCQDPT. Затем хотел бы применить полученные навыки в расчете систем с 1. водородными связями и 2. внутримолекулярным переносом протона, например, гидроксибензалдегида и т.п. Сейчас мне нужны рекомендации:
1. по выбору базиса
2. кол-ва орбиталей
3. общие советы по расчету, чтобы все нормально сходилось

Аватара пользователя
Lantano
Сообщения: 189
Зарегистрирован: Чт окт 09, 2008 10:35 pm

Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч

Сообщение Lantano » Пн май 10, 2010 10:36 am

"Для любой другой несимметирчной молекулы" что-то конкретное посоветовать сложно. Начнем с карбонильных соединений.
Как правило, низшие синглетные и триплетные состояния у них одного типа (n-->pi*). Равновесные геометрические параметры тоже могут быть близки. Поэтому рекомендую начать с низшего триплета Т1. Посчитайте для него параметры минимума на ППЭ. Также рекомендую использовать стартовые орбитали триплета при расчете CASSCF синглета.
Но учтите
1) при расчете синглетного состояния вариационный коллапс в CASSCF никто не отменял. Поэтому либо усреднение по состояниям, что в Firefly появилось не очень давно и может быть пока недоработано, либо надеяться на удачу, правильно выбирая стартовое приближение. На основе богатого личного опыта советую для начала выбрать второй путь. Для этого полезно стартовать с симметрично точки на ППЭ (там возможно S0 и S1 состояния отделены по симметрии и можно сделать правильные орбитали для S1, которые потом использовать как стартовые для несимметирчных структур).
2) если нет сопряжения карбонильной группы с ненасыщенным фрагментом, весьма вероятна пирамидализация карбонильного фрагмента. Отсюда увеличение числа конформеров и прочие проблемы.
В общем случае для алифатических систем достаточно орбиталей, локализованных на карбонильном фрагменте (две сигмы CO, две pi CO и орбиталь неподеленной пары на кислороде). Это самая устойчивая и прилично ведущая себя комбинация. Но возможны и другие варианты. По ацетальдегиду в низших возбужденных электронных состояниях работ опубликовано много.
Салициловый альдегид, боюсь, - большая "бяка". Возможно, в активное пространство придется включать все орбитали pi-типа + еще какие-то сигма. Это надо пробовать. Все зависит от Ваших задач и от требуемой точности к расчету.
Может быть сумбурно получилось, но для начала такие общие советы.
"Если миру не нужен свет - Не беда, проживёт и так.
Если стал не нужен поэт - Не беда, он наймётся в кабак..."

Аватара пользователя
Lantano
Сообщения: 189
Зарегистрирован: Чт окт 09, 2008 10:35 pm

Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч

Сообщение Lantano » Пн май 10, 2010 10:56 am

Вдогонку старый обзор наших работ по расчетам алифатических карбонильных соединений в возбужденных электронных состояниях. Несмотря на то, что это Журнал Физической химии, на русском языке статья никогде не публиковалась.
У вас нет необходимых прав для просмотра вложений в этом сообщении.
"Если миру не нужен свет - Не беда, проживёт и так.
Если стал не нужен поэт - Не беда, он наймётся в кабак..."

lmihael
Сообщения: 227
Зарегистрирован: Вт ноя 21, 2006 3:33 pm

Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч

Сообщение lmihael » Сб май 22, 2010 3:41 pm

Lantano писал(а):"Для любой другой несимметирчной молекулы" что-то конкретное посоветовать сложно. Начнем с карбонильных соединений.
Спасибо большое за ответы, но я даже не могу создать инпут файл :(. В справочной информации CASSCF упоминается только три раза. Не могу найти ни одного примера. Не могли бы вы мне показать такой пример - буду его изучать.

lmihael
Сообщения: 227
Зарегистрирован: Вт ноя 21, 2006 3:33 pm

Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч

Сообщение lmihael » Вт май 25, 2010 1:13 pm

Нашел в сети инпуты Грановского с комментариями - пытаюсь разобраться. Самое интересное, что на свежем FF не сходится тот расчет который нормально сходился на старом.

Ответить

Вернуться в «квантовая химия и моделирование»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 73 гостя