Слышал про CASSCF/MCQDPT, caspt2 и eom-ccsd.Что это такое, чем они отличаются и не одно ли это и тоже
Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для чайни
Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для чайни
Здравствуйте. Надеюсь не очень надоел вам своими вопросами про возбужденное состояние. Хочу попробовать что-то более серьезное чем TDDFT/b3lyp и CIS.
Слышал про CASSCF/MCQDPT, caspt2 и eom-ccsd.Что это такое, чем они отличаются и не одно ли это и тоже
? А также что из них реализовано в Firefly? Каким методом лучше оптимизировать геометрию, а каким определять энергию и дипольный момент первого возбужденого состояния? Какие базисы следует использовать? Был бы очень благодарен за алгоритм расчета на примере молекулы этаналя.
Слышал про CASSCF/MCQDPT, caspt2 и eom-ccsd.Что это такое, чем они отличаются и не одно ли это и тоже
Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч
В малых базисах эти методы работают плохо, в больших базисах не пробовал.
После отстоя требуйте долива
Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч
А в малых это каких? И можете ли вы посоветовать методы, которые работают хорошо?VTur писал(а):В малых базисах эти методы работают плохо, в больших базисах не пробовал.
Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч
В Firefly есть CASSCF с возможностью расчета аналитических градиентов. Т.е. можно "естественным образом" оптимизировать геометрические параметры. Есть теория возмущений на его основе MCQDPT, а также какие-то варианты связанных кластеров. Но в них градиенты аналитически не считаются.
Какое у Вас первое возбужденное состояние? Для ацетальдегида по энергии это триплет n-->pi*. В этом случае выбор методов расчета побогаче.
P.S. Можете написать мне на e-mail. При более предметном разговоре может быть что-нибудь дельное посоветую.
Какое у Вас первое возбужденное состояние? Для ацетальдегида по энергии это триплет n-->pi*. В этом случае выбор методов расчета побогаче.
P.S. Можете написать мне на e-mail. При более предметном разговоре может быть что-нибудь дельное посоветую.
"Если миру не нужен свет - Не беда, проживёт и так.
Если стал не нужен поэт - Не беда, он наймётся в кабак..."
Если стал не нужен поэт - Не беда, он наймётся в кабак..."
Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч
Очень бы хотелось на примере ацетальдегида или любой другой несимметричной молекулы оптимизировать геометрию первого возбужденного синглетного состояния (S1) методом CASSCF, а затем рассчитать дипольный момент и полную энергию этого состояния в MCQDPT. Затем хотел бы применить полученные навыки в расчете систем с 1. водородными связями и 2. внутримолекулярным переносом протона, например, гидроксибензалдегида и т.п. Сейчас мне нужны рекомендации:Lantano писал(а):В Firefly есть CASSCF с возможностью расчета аналитических градиентов. Т.е. можно "естественным образом" оптимизировать геометрические параметры. Есть теория возмущений на его основе MCQDPT, а также какие-то варианты связанных кластеров. Но в них градиенты аналитически не считаются.
Какое у Вас первое возбужденное состояние? Для ацетальдегида по энергии это триплет n-->pi*. В этом случае выбор методов расчета побогаче.
P.S. Можете написать мне на e-mail. При более предметном разговоре может быть что-нибудь дельное посоветую.
1. по выбору базиса
2. кол-ва орбиталей
3. общие советы по расчету, чтобы все нормально сходилось
Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч
"Для любой другой несимметирчной молекулы" что-то конкретное посоветовать сложно. Начнем с карбонильных соединений.
Как правило, низшие синглетные и триплетные состояния у них одного типа (n-->pi*). Равновесные геометрические параметры тоже могут быть близки. Поэтому рекомендую начать с низшего триплета Т1. Посчитайте для него параметры минимума на ППЭ. Также рекомендую использовать стартовые орбитали триплета при расчете CASSCF синглета.
Но учтите
1) при расчете синглетного состояния вариационный коллапс в CASSCF никто не отменял. Поэтому либо усреднение по состояниям, что в Firefly появилось не очень давно и может быть пока недоработано, либо надеяться на удачу, правильно выбирая стартовое приближение. На основе богатого личного опыта советую для начала выбрать второй путь. Для этого полезно стартовать с симметрично точки на ППЭ (там возможно S0 и S1 состояния отделены по симметрии и можно сделать правильные орбитали для S1, которые потом использовать как стартовые для несимметирчных структур).
2) если нет сопряжения карбонильной группы с ненасыщенным фрагментом, весьма вероятна пирамидализация карбонильного фрагмента. Отсюда увеличение числа конформеров и прочие проблемы.
В общем случае для алифатических систем достаточно орбиталей, локализованных на карбонильном фрагменте (две сигмы CO, две pi CO и орбиталь неподеленной пары на кислороде). Это самая устойчивая и прилично ведущая себя комбинация. Но возможны и другие варианты. По ацетальдегиду в низших возбужденных электронных состояниях работ опубликовано много.
Салициловый альдегид, боюсь, - большая "бяка". Возможно, в активное пространство придется включать все орбитали pi-типа + еще какие-то сигма. Это надо пробовать. Все зависит от Ваших задач и от требуемой точности к расчету.
Может быть сумбурно получилось, но для начала такие общие советы.
Как правило, низшие синглетные и триплетные состояния у них одного типа (n-->pi*). Равновесные геометрические параметры тоже могут быть близки. Поэтому рекомендую начать с низшего триплета Т1. Посчитайте для него параметры минимума на ППЭ. Также рекомендую использовать стартовые орбитали триплета при расчете CASSCF синглета.
Но учтите
1) при расчете синглетного состояния вариационный коллапс в CASSCF никто не отменял. Поэтому либо усреднение по состояниям, что в Firefly появилось не очень давно и может быть пока недоработано, либо надеяться на удачу, правильно выбирая стартовое приближение. На основе богатого личного опыта советую для начала выбрать второй путь. Для этого полезно стартовать с симметрично точки на ППЭ (там возможно S0 и S1 состояния отделены по симметрии и можно сделать правильные орбитали для S1, которые потом использовать как стартовые для несимметирчных структур).
2) если нет сопряжения карбонильной группы с ненасыщенным фрагментом, весьма вероятна пирамидализация карбонильного фрагмента. Отсюда увеличение числа конформеров и прочие проблемы.
В общем случае для алифатических систем достаточно орбиталей, локализованных на карбонильном фрагменте (две сигмы CO, две pi CO и орбиталь неподеленной пары на кислороде). Это самая устойчивая и прилично ведущая себя комбинация. Но возможны и другие варианты. По ацетальдегиду в низших возбужденных электронных состояниях работ опубликовано много.
Салициловый альдегид, боюсь, - большая "бяка". Возможно, в активное пространство придется включать все орбитали pi-типа + еще какие-то сигма. Это надо пробовать. Все зависит от Ваших задач и от требуемой точности к расчету.
Может быть сумбурно получилось, но для начала такие общие советы.
"Если миру не нужен свет - Не беда, проживёт и так.
Если стал не нужен поэт - Не беда, он наймётся в кабак..."
Если стал не нужен поэт - Не беда, он наймётся в кабак..."
Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч
Вдогонку старый обзор наших работ по расчетам алифатических карбонильных соединений в возбужденных электронных состояниях. Несмотря на то, что это Журнал Физической химии, на русском языке статья никогде не публиковалась.
У вас нет необходимых прав для просмотра вложений в этом сообщении.
"Если миру не нужен свет - Не беда, проживёт и так.
Если стал не нужен поэт - Не беда, он наймётся в кабак..."
Если стал не нужен поэт - Не беда, он наймётся в кабак..."
Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч
Спасибо большое за ответы, но я даже не могу создать инпут файлLantano писал(а):"Для любой другой несимметирчной молекулы" что-то конкретное посоветовать сложно. Начнем с карбонильных соединений.
Re: Расчет возб.сост.(не TDDFT) в Firefly.Нужен ликбез для ч
Нашел в сети инпуты Грановского с комментариями - пытаюсь разобраться. Самое интересное, что на свежем FF не сходится тот расчет который нормально сходился на старом.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 73 гостя