DFT расчет для молекулярного кислорода
DFT расчет для молекулярного кислорода
Провожу расчет для синглетного и триплетного кислорода. Длина связи получается одинаковой, энергия S-T перехода получается порядка 2 эВ, в то время как должна быть, насколько я понимаю, около 1 эВ.
Считаю в B3LYP. Помогите разобраться.
Интересует именно DFT. Просьба, MCSCF и прочие более адекватные методы не предлагать)
Пример для синглетного кислорода прилагается.
Считаю в B3LYP. Помогите разобраться.
Интересует именно DFT. Просьба, MCSCF и прочие более адекватные методы не предлагать)
Пример для синглетного кислорода прилагается.
У вас нет необходимых прав для просмотра вложений в этом сообщении.
Re: DFT расчет для молекулярного кислорода
А в чем проблема?
А я вот паровоз поднимал... Но не поднял.
Re: DFT расчет для молекулярного кислорода
У Вас получился O2 (1-Sigma). Напрасно Вы так отсекаете более адекватные методы. Хотя, может, именно это состояние Вам и надо. Нужно нарушить симметрию для получения O2 (1-delta).
Re: DFT расчет для молекулярного кислорода
И как же нарушить симметрию?o-oxhem писал(а):У Вас получился O2 (1-Sigma). Напрасно Вы так отсекаете более адекватные методы. Хотя, может, именно это состояние Вам и надо. Нужно нарушить симметрию для получения O2 (1-delta).
Re: DFT расчет для молекулярного кислорода
Для расчёта синглетных бирадикалов гибридными DFT-методами нужно просто обязательно ставить guess mix (нарушение симметрии альфа- и бета- электронов). Хотя и это не гарантирует, что волновая функция не свалится в симметричную.
Re: DFT расчет для молекулярного кислорода
огромное спасибо за совет.EvgeniX писал(а):Для расчёта синглетных бирадикалов гибридными DFT-методами нужно просто обязательно ставить guess mix (нарушение симметрии альфа- и бета- электронов). Хотя и это не гарантирует, что волновая функция не свалится в симметричную.
но теперь энергия перехода получается 0,44 эВ)) что тоже как-то не очень.
Re: DFT расчет для молекулярного кислорода
Вместе с approximate spincorrection procedure proposed by Yamaguchi получается 0,89эВ. Однако E синглетного кислорода кажется заниженной.
Последний раз редактировалось ruderal Чт дек 08, 2011 10:43 pm, всего редактировалось 1 раз.
Re: DFT расчет для молекулярного кислорода
И еще такой вопрос: насколько адекватно B3LYP позволяет рассчитать активные формы кислорода, например, супероксид-анион радикал?
Re: DFT расчет для молекулярного кислорода
Ответ: скорее нет, чем да. Такую мелочь как активные формы кислорода сам бог велел изучать в CASSCF в приличном базисе. В радикалах и тем более бирадикалах работает статическая корреляция -- та самая, что учитывается в CASSCF. Добавите учет динамической корреляции (пока самый надежный способ это сделать -- XMCQDPT в FireFly) -- и будет Вам счастье. Альтернатива -- Full CI (для кислорода или супероксида это реально).
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)
Re: DFT расчет для молекулярного кислорода
Все это замечательно, но проблема в том, что меня интересует расчет ∆Е реакций молекулярного кислорода и АФК с более крупными молекулами (30 атомов). АФК не самоцель.sanya1024 писал(а):Ответ: скорее нет, чем да. Такую мелочь как активные формы кислорода сам бог велел изучать в CASSCF в приличном базисе. В радикалах и тем более бирадикалах работает статическая корреляция -- та самая, что учитывается в CASSCF. Добавите учет динамической корреляции (пока самый надежный способ это сделать -- XMCQDPT в FireFly) -- и будет Вам счастье. Альтернатива -- Full CI (для кислорода или супероксида это реально).
Re: DFT расчет для молекулярного кислорода
30 атомов -- фигня для CASSCF (задачи до 1000 баз. функций у меня нормально считались на кластере). В FireFly CASSCF отлично распараллелен, в FF8 Грановский ввел новый эффективный параллельный алгоритм для детерминантного CASSCF (ALDET), что позволяет работать с довольно большими (до (17,17)) акт. пространствами.
Конечно, Full CI для 30 атомов отпадает.
Конечно, Full CI для 30 атомов отпадает.
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)
Re: DFT расчет для молекулярного кислорода
Хорошо, большое спасибо. Я пока не очень владею этими методами. Так что это все будет задача на будущее.sanya1024 писал(а):30 атомов -- фигня для CASSCF (задачи до 1000 баз. функций у меня нормально считались на кластере). В FireFly CASSCF отлично распараллелен, в FF8 Грановский ввел новый эффективный параллельный алгоритм для детерминантного CASSCF (ALDET), что позволяет работать с довольно большими (до (17,17)) акт. пространствами.
Конечно, Full CI для 30 атомов отпадает.
Вопрос: оптимизацию геометрии крупных молекул тоже проводить в CASSCF?
Re: DFT расчет для молекулярного кислорода
Ну, я оптимизировала комплекс европия с четырьмя органическими лигандами. Потому и говорю, что 1000 баз. функций вполне влезают.
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 11 гостей