Можем ли мы моделировать химические реакции с помощью квантовой механики? Если да, то какую самую сложную реакцию мы можем смоделировать?

Не физик и не химик, просто интересуюсь QM и его приложениями.

Недавно я читал о квантовой химии, и мне пришло в голову, что, поскольку мы можем моделировать электронные орбитали в квантовой механике и определять форму получаемых орбиталей (S, P, D, F и т. д.), это объясняет, почему нейтральные атомы могут образовывать связи вместе (волновая функция валентного электрона распространяется в более низкое энергетическое состояние между атомами, а разница переходит в кинетическую энергию).

Итак, мне приходит в голову, что мы должны уметь представлять химические реакции как некое взаимодействие между группами атомов. Итак, возможно ли это реально (вычислительно)? И если да, то какой комплекс реакций мы можем представить в КМ? Пытаемся ли мы смоделировать H2 или реакцию 4H + O2 => 2H2O, или компьютеры находятся на этапе, когда мы можем моделировать метилирование бензола или синтез АТФ?

Бонусные баллы, если кто-то может дать ссылку на видеоролики с визуализацией реакций, демонстрирующих их квантовую природу!

Спасибо!

К 2009 году квантовая механика может смоделировать 1 миллион атомов sciencedirect.com/science/article/pii/S0009261409005041 , я думаю, что для реакции это будет включать более сложные детали, тем не менее диапазон может быть аналогичным
Ну, это во многом зависит от того, какая точность вам требуется. Если вы хотите, чтобы уровни энергии и полные волновые функции соответствовали машинной точности, мне сказали, что текущий практический предел составляет около 6 или 7 частиц, что очень трудно превзойти из-за экспоненциальности гильбертова пространства. Группа Джима Митроя в ХДС в этом профессионалы. Методы для большего числа частиц основаны на определенных схемах приближения. Квантовый химик должен знать, какие схемы полезны для каких практических целей, например, DFT, как упоминал @user26143.
Эй, что! Я пошел в ХДС! Я видел Джима Митроя, но я не занимался инженерией, поэтому он никогда не читал мне лекции. Маленький мир! До сих пор я понятия не имел, что он был исследователем QM, не говоря уже о мировом уровне. Удивительный.
Я бы с осторожностью относился к сравнению производительности симуляции металла (как в ссылке на симуляцию «1 mil атом» выше) с химией. Я не квантовый химик, но я думаю, что в металлическом случае расположение атомов фиксировано и оценивается электронная плотность, в то время как то, что спрашивает OP, должно включать как моделирование электронной орбиты, так и динамику молекулярных / атомов. Эта комбинация называется QM/MM (или QMMM) и используется для имитации, например, активных сайтов ферментов вместе с их движением, но у меня нет никаких данных о производительности.
Другая проблема заключается в том, принадлежит ли MM QM, как спросил OP.
Относительно Джима Митроя: я сотрудник JCU , поэтому я встретил его только один раз, когда посетил Дарвин (именно тогда он рассказал мне о пределе в 6-7 частиц), но он сотрудничает с несколькими другими нашими аспирантами по переносу позитронов в мягкое уплотненное вещество. Я думаю, что он занимается действительно квантовыми вещами, а они занимаются транспортной теорией.
@MichaelBrown, интересно, возможно, мне придется попытаться попасть в его рабочие часы.

Ответы (1)

Да, с помощью КМ можно моделировать химические реакции. В квантовой химии часто решают уравнение Шредингера молекулярного гамильтониана, предполагая приближение Борна-Оппенгеймера. В этом приближении энергия является функцией координат атомных ядер. Таким образом, для моделирования реакции вы можете следовать пути с наименьшей энергией от координат атомов реагентов к координатам продуктов.

Теперь, насколько сложными могут быть реакции, и при этом их можно смоделировать с вычислительной точки зрения? Это зависит от двух вещей: (1) какая точность вам нужна? (2) какие эффекты участвуют в реакции? Например, с современными компьютерами и с помощью теории функционала плотности (DFT) вполне реально смоделировать метилирование бензола (или даже гораздо более крупных систем), но я бы не ожидал точности лучше, чем около 4 ккал/моль в модели. энергии. Однако я не знаю ни одного метода, способного дать точное описание диссоциации Cr. 2 димер, хотя более крупные системы могут быть точно смоделированы различными методами. Причина этого в том, что эффекты электронной корреляции очень сильны в Cr 2 , и эти эффекты по своей природе трудно и вычислительно дорого учитывать.

Если вам нужен наилучший возможный результат (т. е. учет всех эффектов корреляции), вам необходимо использовать метод, называемый полным взаимодействием конфигурации (FCI). В FCI вы минимизируете энергию волновой функции, которая представляет собой линейную комбинацию всех возможных конфигураций (детерминантов Слейтера), которые может принять система. Это самая общая волновая функция из всех возможных, и, таким образом, согласно вариационному принципу, когда вы минимизируете ее энергию, вы получаете точный численный результат для уравнения Шредингера для заданного базисного набора (то есть функций, которые вы используете для представления своих орбиталей). Однако временная сложность этого метода составляет примерно О ( М ! ) , где М есть число базисных функций, и результат действительно точен только тогда, когда М стремится к бесконечности. На практике, М должны быть очень большими, чтобы сходиться, но это все еще настолько дорого, что мы можем позволить себе делать FCI только для систем с не более чем 10 электронами. Вот почему в квантовой химии так много методов приближения; так что можно использовать метод, соответствующий размеру системы и требуемой точности. ТПФ чрезвычайно популярен, потому что он относительно дешев и точен (формально О ( М 4 ) и средние ошибки около 4 ккал/моль в энергиях, как указано выше). Другими хорошими методами, которые находятся между FCI и DFT по точности и стоимости, являются методы связанных кластеров ( О ( М 7 ) ) и CASPT2 (комбинаторная стоимость), который является своего рода FCI в подмножестве электронов и орбиталей и добавляет поверх этого поправку теории возмущений второго порядка.

Что касается примеров реакций, вам, вероятно, следует проверить правила Вудворда-Хоффмана ( http://en.wikipedia.org/wiki/Woodward%E2%80%93Hoffmann_rules ). Это несколько простых правил, основанных на квантовой механике и теории молекулярных орбиталей, которые могут предсказать исход множества реакций, которые нельзя объяснить никаким другим способом.

Может быть, вы могли бы взглянуть на этот вопрос, который меня интересует. Есть ли у этих методов что-нибудь сказать?
@WetSavannaAnimalakaRodVance Я разместил ответ на вопрос по вашей ссылке. Вы можете строго моделировать любые спектры поглощения/излучения любым методом, который дает волновую функцию. Все методы, которые я упомянул, дают волновые функции, за исключением DFT, который дает плотности. Однако вы можете использовать TDDFT ( en.wikipedia.org/wiki/Time-dependent_density_functional_theory ) или Δ SCF ( pubs.acs.org/doi/abs/10.1021/ct400356k ) для расчета спектров с помощью DFT.