Вычислительная химия — это раздел химии , который использует компьютерное моделирование для решения химических задач. [2] Она использует методы теоретической химии, встроенные в компьютерные программы, для расчета структур и свойств молекул , групп молекул и твердых тел. [3] Важность этого предмета проистекает из того факта, что, за исключением некоторых относительно недавних открытий, связанных с молекулярным ионом водорода ( дигидрогенным катионом ), достижение точного квантово-механического описания химических систем аналитически или в замкнутой форме не представляется возможным. [4] Сложность, присущая проблеме многих тел, усугубляет проблему предоставления подробных описаний квантово-механических систем. [5] Хотя результаты вычислений обычно дополняют информацию, полученную в ходе химических экспериментов , они иногда могут предсказывать ненаблюдаемые химические явления . [6]
Вычислительная химия отличается от теоретической химии , которая включает математическое описание химии. Однако вычислительная химия включает использование компьютерных программ и дополнительных математических навыков для точного моделирования различных химических проблем. В теоретической химии химики, физики и математики разрабатывают алгоритмы и компьютерные программы для прогнозирования атомных и молекулярных свойств и путей химических реакций. Вычислительные химики, напротив, могут просто применять существующие компьютерные программы и методологии к конкретным химическим вопросам. [7]
Исторически вычислительная химия имела два различных аспекта:
Эти аспекты, наряду с целью вычислительной химии, привели к появлению целого ряда алгоритмов.
Основываясь на основополагающих открытиях и теориях в истории квантовой механики , первые теоретические расчеты в химии были сделаны Уолтером Гайтлером и Фритцем Лондоном в 1927 году с использованием теории валентных связей . [9] Книги, оказавшие влияние на раннее развитие вычислительной квантовой химии, включают « Введение в квантовую механику – с приложениями к химии» Лайнуса Полинга и Э. Брайта Уилсона 1935 года , [10] «Квантовая химия» Айринга , Уолтера и Кимбалла 1944 года , [11] «Элементарная волновая механика – с приложениями к квантовой химии» Гайтлера 1945 года , [12] и позднее учебник Коулсона 1952 года «Валентность» , каждая из которых служила основным источником информации для химиков в последующие десятилетия. [13]
С развитием эффективной компьютерной технологии в 1940-х годах решения сложных волновых уравнений для сложных атомных систем стали достижимой целью. В начале 1950-х годов были выполнены первые полуэмпирические расчеты атомных орбиталей. Теоретические химики стали активными пользователями первых цифровых компьютеров. Одно из значительных достижений было отмечено статьей Клеменса К. Дж. Рутана 1951 года в Reviews of Modern Physics. [14] [15] Эта статья в основном была сосредоточена на подходе «LCAO MO» (линейная комбинация атомных орбиталей молекулярных орбиталей). В течение многих лет она была второй по частоте цитирования статьей в этом журнале. [14] [15] Очень подробный отчет о таком использовании в Соединенном Королевстве дан Смитом и Сатклиффом. [16] Первые ab initio расчеты по методу Хартри-Фока для двухатомных молекул были выполнены в 1956 году в Массачусетском технологическом институте с использованием базисного набора орбиталей Слейтера . [17] Для двухатомных молекул систематическое исследование с использованием минимального базисного набора и первый расчет с большим базисным набором были опубликованы Рансилом и Несбетом соответственно в 1960 году. [18] Первые полиатомные расчеты с использованием гауссовых орбиталей были выполнены в конце 1950-х годов. Первые расчеты конфигурационного взаимодействия были выполнены в Кембридже на компьютере EDSAC в 1950-х годах с использованием гауссовых орбиталей Бойзом и его коллегами. [19] К 1971 году, когда была опубликована библиография ab initio расчетов, [20] самыми большими включенными молекулами были нафталин и азулен . [21] [22] Рефераты многих более ранних разработок в области ab initio теории были опубликованы Шефером. [23]
В 1964 году на компьютерах в Беркли и Оксфорде были проведены расчеты по методу Хюккеля (с использованием метода простой линейной комбинации атомных орбиталей (ЛКАО) для определения электронных энергий молекулярных орбиталей π-электронов в сопряженных углеводородных системах) молекул, начиная от бутадиена и бензола и заканчивая оваленом . [24] Эти эмпирические методы были заменены в 1960-х годах полуэмпирическими методами, такими как CNDO . [25]
В начале 1970-х годов эффективные компьютерные программы ab initio, такие как ATMOL, Gaussian , IBMOL и POLYAYTOM, начали использоваться для ускорения ab initio расчетов молекулярных орбиталей. [26] Из этих четырех программ только Gaussian, теперь значительно расширенная, все еще используется, но сейчас используются и многие другие программы. [26] В то же время были разработаны методы молекулярной механики , такие как силовое поле MM2 , в первую очередь Норманом Аллингером . [27]
Одно из первых упоминаний термина «вычислительная химия» можно найти в книге 1970 года « Компьютеры и их роль в физических науках» Сиднея Фернбаха и Абрахама Хаскелла Тауба, где они утверждают: «Поэтому кажется, что «вычислительная химия» может наконец стать все более и более реальностью». [28] В 1970-х годах самые разные методы стали рассматриваться как часть новой зарождающейся дисциплины вычислительной химии . [29] Журнал «Вычислительная химия» был впервые опубликован в 1980 году.
Вычислительная химия была отмечена несколькими Нобелевскими премиями, наиболее заметными из которых были премии 1998 и 2013 годов. Вальтер Кон , «за разработку теории функционала плотности», и Джон Попл , «за разработку вычислительных методов в квантовой химии», получили Нобелевскую премию по химии 1998 года. [30] Мартин Карплус , Майкл Левитт и Арье Варшел получили Нобелевскую премию по химии 2013 года за «разработку многомасштабных моделей для сложных химических систем». [31]
В вычислительной химии существует несколько направлений.
Эти поля могут иметь несколько применений, как показано ниже.
Вычислительная химия — это инструмент для анализа каталитических систем без проведения экспериментов. Современная электронная теория структуры и теория функционала плотности позволили исследователям обнаружить и понять катализаторы . [37] Вычислительные исследования применяют теоретическую химию к исследованию катализа. Методы теории функционала плотности вычисляют энергии и орбитали молекул, чтобы получить модели этих структур. [38] Используя эти методы, исследователи могут предсказывать такие значения, как энергия активации , реакционная способность сайта [39] и другие термодинамические свойства. [38]
Данные, которые трудно получить экспериментально, можно найти, используя вычислительные методы для моделирования механизмов каталитических циклов. [39] Опытные вычислительные химики дают прогнозы, близкие к экспериментальным данным с надлежащим учетом методов и базисных наборов. С хорошими вычислительными данными исследователи могут предсказать, как можно улучшить катализаторы, чтобы снизить стоимость и повысить эффективность этих реакций. [38]
Вычислительная химия используется в разработке лекарств для моделирования потенциально полезных молекул лекарств и помогает компаниям экономить время и деньги при разработке лекарств. Процесс открытия лекарств включает анализ данных, поиск способов улучшения текущих молекул, поиск синтетических путей и тестирование этих молекул. [36] Вычислительная химия помогает в этом процессе, давая прогнозы о том, какие эксперименты лучше всего провести, не проводя другие эксперименты. Вычислительные методы также могут находить значения, которые трудно найти экспериментально, такие как pKa соединений. [40] Такие методы, как теория функционала плотности, могут использоваться для моделирования молекул лекарств и нахождения их свойств, таких как их энергии HOMO и LUMO и молекулярные орбитали. Вычислительные химики также помогают компаниям в разработке информатики, инфраструктуры и дизайна лекарств. [41]
Помимо синтеза лекарств, носители лекарств также исследуются вычислительными химиками для наноматериалов . Это позволяет исследователям моделировать среды для проверки эффективности и стабильности носителей лекарств. Понимание того, как вода взаимодействует с этими наноматериалами, обеспечивает стабильность материала в организме человека. Эти вычислительные симуляции помогают исследователям оптимизировать материал, найти лучший способ структурировать эти наноматериалы перед их изготовлением. [42]
Базы данных полезны как для вычислительных, так и для невычислительных химиков в исследованиях и проверке достоверности вычислительных методов. Эмпирические данные используются для анализа погрешности вычислительных методов по сравнению с экспериментальными данными. Эмпирические данные помогают исследователям с их методами и базисными наборами иметь большую уверенность в результатах исследователей. Базы данных вычислительной химии также используются при тестировании программного обеспечения или оборудования для вычислительной химии. [43]
Базы данных также могут использовать чисто расчетные данные. Чисто расчетные данные используют расчетные значения вместо экспериментальных значений для баз данных. Чисто расчетные данные избегают работы с этими корректировками для различных экспериментальных условий, таких как энергия нулевой точки. Эти расчеты также могут избежать экспериментальных ошибок для сложных для тестирования молекул. Хотя чисто расчетные данные часто не идеальны, выявление проблем часто проще для расчетных данных, чем для экспериментальных. [43]
Базы данных также предоставляют публичный доступ к информации для использования исследователями. Они содержат данные, которые другие исследователи нашли и загрузили в эти базы данных, чтобы любой мог их искать. Исследователи используют эти базы данных, чтобы найти информацию о молекулах, представляющих интерес, и узнать, что можно сделать с этими молекулами. Некоторые общедоступные химические базы данных включают следующее. [43]
Программы, используемые в вычислительной химии, основаны на многих различных квантово-химических методах, которые решают молекулярное уравнение Шредингера, связанное с молекулярным гамильтонианом . [46] Методы, которые не включают никаких эмпирических или полуэмпирических параметров в свои уравнения, будучи выведенными непосредственно из теории, без включения экспериментальных данных, называются методами ab initio . [47] Теоретическое приближение строго определяется на первых принципах, а затем решается в пределах погрешности, которая качественно известна заранее. Если необходимо использовать численные итерационные методы, цель состоит в том, чтобы повторять до тех пор, пока не будет получена полная машинная точность (наилучшее, что возможно при конечной длине слова на компьютере и в пределах сделанных математических и/или физических приближений). [48]
Методы ab initio должны определять уровень теории (метод) и базисный набор. [49] Базисный набор состоит из функций, центрированных на атомах молекулы. Затем эти наборы используются для описания молекулярных орбиталей с помощью линейной комбинации атомных орбиталей (LCAO) молекулярный орбитальный метод ansatz . [50]
Распространенным типом расчета электронной структуры ab initio является метод Хартри-Фока (HF), расширение теории молекулярных орбиталей , где отталкивание электронов в молекуле специально не учитывается; в расчет включается только средний эффект электронов. По мере увеличения размера базисного набора энергия и волновая функция стремятся к пределу, называемому пределом Хартри-Фока. [50]
Многие типы вычислений начинаются с вычисления Хартри-Фока и впоследствии корректируются для электрон-электронного отталкивания, называемого также электронной корреляцией . [51] Эти типы вычислений называются пост-Хартри-Фоковскими методами. Постоянно совершенствуя эти методы, ученые могут все больше приближаться к идеальному предсказанию поведения атомных и молекулярных систем в рамках квантовой механики, как определено уравнением Шредингера. [52] Чтобы получить точное согласие с экспериментом, необходимо включить определенные термины, некоторые из которых гораздо важнее для тяжелых атомов, чем для легких. [53]
В большинстве случаев волновая функция Хартри-Фока занимает одну конфигурацию или определитель. [54] В некоторых случаях, особенно для процессов разрыва связей, этого недостаточно, и необходимо использовать несколько конфигураций . [55]
Полная молекулярная энергия может быть оценена как функция молекулярной геометрии ; другими словами, поверхность потенциальной энергии . [56] Такая поверхность может быть использована для динамики реакции. Стационарные точки поверхности приводят к предсказаниям различных изомеров и переходных структур для конверсии между изомерами, но их можно определить без полного знания полной поверхности. [53]
Особенно важной задачей, называемой вычислительной термохимией , является расчет термохимических величин, таких как энтальпия образования, с химической точностью. Химическая точность — это точность, необходимая для создания реалистичных химических прогнозов, и обычно считается равной 1 ккал/моль или 4 кДж/моль. Чтобы достичь этой точности экономичным способом, необходимо использовать ряд методов пост-Хартри–Фока и объединить результаты. Эти методы называются композитными методами квантовой химии . [57]
После разделения электронных и ядерных переменных в представлении Борна-Оппенгеймера) волновой пакет, соответствующий ядерным степеням свободы, распространяется посредством оператора эволюции времени (физики), связанного с зависящим от времени уравнением Шредингера (для полного молекулярного гамильтониана ). [58] В подходе, зависящем от дополнительной энергии, независимое от времени уравнение Шредингера решается с использованием формализма теории рассеяния . Потенциал, представляющий межатомное взаимодействие, задается поверхностями потенциальной энергии . В общем случае поверхности потенциальной энергии связаны посредством членов вибронной связи . [59]
Наиболее популярными методами распространения волнового пакета , связанными с молекулярной геометрией, являются:
То, как вычислительный метод решает квантовые уравнения, влияет на точность и эффективность метода. Метод разделения операторов является одним из таких методов решения дифференциальных уравнений. В вычислительной химии метод разделения операторов снижает вычислительные затраты на моделирование химических систем. Вычислительные затраты — это то, сколько времени требуется компьютерам для расчета этих химических систем, поскольку для более сложных систем это может занять несколько дней. Квантовые системы сложны и требуют много времени для решения людьми. Методы разделения операторов помогают компьютерам быстро рассчитывать эти системы, решая подзадачи в квантовом дифференциальном уравнении . Метод делает это путем разделения дифференциального уравнения на два разных уравнения, например, когда имеется более двух операторов. После решения разделенные уравнения снова объединяются в одно уравнение, чтобы дать легко вычисляемое решение. [62]
Этот метод используется во многих областях, где требуется решение дифференциальных уравнений, например, в биологии . Однако, этот метод сопровождается ошибкой расщепления. Например, со следующим решением для дифференциального уравнения. [62]
Уравнение можно разделить, но решения не будут точными, только похожими. Это пример расщепления первого порядка. [62]
Существуют способы уменьшить эту ошибку, в том числе взять среднее значение двух разделенных уравнений. [62]
Другой способ повышения точности — использовать расщепление более высокого порядка. Обычно расщепление второго порядка — это самое большее, что делается, потому что расщепление более высокого порядка требует гораздо больше времени для расчета и не стоит затрат. Методы более высокого порядка становятся слишком сложными для реализации и бесполезны для решения дифференциальных уравнений, несмотря на более высокую точность. [62]
Химики-вычислители тратят много времени на то, чтобы сделать системы, рассчитанные с помощью техники оператора разделения, более точными, минимизируя при этом вычислительные затраты. Методы расчета представляют собой сложную задачу для многих химиков, пытающихся моделировать молекулы или химические среды. [62]
Методы теории функционала плотности (DFT) часто считаются методами ab initio для определения молекулярной электронной структуры, хотя многие из наиболее распространенных функционалов используют параметры, полученные из эмпирических данных или из более сложных вычислений. В DFT полная энергия выражается в терминах полной одноэлектронной плотности, а не волновой функции. В этом типе вычислений есть приближенный гамильтониан и приближенное выражение для полной электронной плотности. Методы DFT могут быть очень точными при небольших вычислительных затратах. Некоторые методы объединяют функционал обмена функционала плотности с обменным членом Хартри-Фока и называются гибридными функциональными методами. [63]
Полуэмпирические методы квантовой химии основаны на формализме метода Хартри-Фока , но делают много приближений и получают некоторые параметры из эмпирических данных. Они были очень важны в вычислительной химии с 60-х по 90-е годы, особенно для обработки больших молекул, где полный метод Хартри-Фока без приближений был слишком затратным. Использование эмпирических параметров, по-видимому, позволяет включить в методы некоторое количество эффектов корреляции. [64]
Примитивные полуэмпирические методы были разработаны еще раньше, где двухэлектронная часть гамильтониана явно не включена. Для π-электронных систем это был метод Хюккеля , предложенный Эрихом Хюккелем , а для всех валентных электронных систем — расширенный метод Хюккеля , предложенный Роальдом Хоффманном . Иногда методы Хюккеля называют «полностью эмпирическими», поскольку они не выводятся из гамильтониана. [65] Тем не менее, термин «эмпирические методы» или «эмпирические силовые поля» обычно используется для описания молекулярной механики. [66]
Во многих случаях большие молекулярные системы можно успешно моделировать, полностью избегая квантово-механических расчетов. Например, молекулярно-механические симуляции используют одно классическое выражение для энергии соединения, например, гармонический осциллятор . Все константы, появляющиеся в уравнениях, должны быть получены заранее из экспериментальных данных или расчетов ab initio . [64]
База данных соединений, используемых для параметризации, т. е. результирующий набор параметров и функций называется силовым полем , имеет решающее значение для успеха расчетов молекулярной механики. Силовое поле, параметризованное против определенного класса молекул, например, белков, как ожидается, будет иметь какое-либо значение только при описании других молекул того же класса. [64] Эти методы могут быть применены к белкам и другим крупным биологическим молекулам и позволяют изучать подход и взаимодействие (стыковку) потенциальных молекул лекарств. [67] [68]
Молекулярная динамика (МД) использует либо квантовую механику , либо молекулярную механику , либо их смесь для расчета сил, которые затем используются для решения законов движения Ньютона для изучения поведения систем, зависящего от времени. Результатом моделирования молекулярной динамики является траектория, которая описывает, как положение и скорость частиц изменяются со временем. Фазовая точка системы, описываемая положениями и импульсами всех ее частиц в предыдущей точке времени, определит следующую фазовую точку во времени путем интегрирования по законам движения Ньютона. [69]
Монте-Карло (МК) генерирует конфигурации системы, внося случайные изменения в положения ее частиц, вместе с их ориентацией и конформацией, где это уместно. [70] Это метод случайной выборки, который использует так называемую выборку по важности . Методы выборки по важности способны генерировать состояния с низкой энергией, поскольку это позволяет точно рассчитывать свойства. Потенциальная энергия каждой конфигурации системы может быть рассчитана вместе со значениями других свойств из положений атомов. [71] [72]
QM/MM — это гибридный метод, который пытается объединить точность квантовой механики со скоростью молекулярной механики. Он полезен для моделирования очень больших молекул, таких как ферменты . [73]
Квантовая вычислительная химия нацелена на использование квантовых вычислений для моделирования химических систем, отличаясь от подхода QM/MM (квантовая механика/молекулярная механика). [74] В то время как QM/MM использует гибридный подход, объединяя квантовую механику для части системы с классической механикой для остальной части, квантовая вычислительная химия использует исключительно квантовые методы вычислений для представления и обработки информации, такие как гамильтоновы операторы. [75]
Обычные методы вычислительной химии часто сталкиваются со сложными квантово-механическими уравнениями, особенно из-за экспоненциального роста волновой функции квантовой системы. Квантовая вычислительная химия решает эти проблемы, используя квантовые вычислительные методы , такие как кубитизация и квантовая фазовая оценка , которые, как считается, предлагают масштабируемые решения. [76]
Квантовая фазовая оценка, с другой стороны, помогает точно определить собственные энергетические состояния, которые имеют решающее значение для понимания поведения квантовой системы. [77]
Хотя эти методы продвинули область вычислительной химии, особенно в моделировании химических систем, их практическое применение в настоящее время ограничено в основном более мелкими системами из-за технологических ограничений. Тем не менее, эти разработки могут привести к значительному прогрессу в достижении более точных и ресурсоэффективных симуляций квантовой химии. [76]
Вычислительная стоимость и алгоритмическая сложность в химии используются для понимания и прогнозирования химических явлений. Они помогают определить, какие алгоритмы/вычислительные методы использовать при решении химических проблем. В этом разделе основное внимание уделяется масштабированию вычислительной сложности в зависимости от размера молекулы и подробно описываются алгоритмы, обычно используемые в обеих областях. [78]
В квантовой химии, в частности, сложность может расти экспоненциально с числом электронов, вовлеченных в систему. Этот экспоненциальный рост является существенным препятствием для точного моделирования больших или сложных систем. [79]
Продвинутые алгоритмы в обеих областях стремятся сбалансировать точность с вычислительной эффективностью. Например, в МД такие методы, как интеграция Верле или алгоритм Бимана, используются для их вычислительной эффективности. В квантовой химии гибридные методы, объединяющие различные вычислительные подходы (например, КМ/ММ), все чаще используются для решения больших биомолекулярных систем. [80]
Следующий список иллюстрирует влияние вычислительной сложности на алгоритмы, используемые в химических вычислениях. Важно отметить, что хотя этот список и содержит ключевые примеры, он не является исчерпывающим и служит лишь руководством для понимания того, как вычислительные требования влияют на выбор конкретных вычислительных методов в химии.
Решает уравнения движения Ньютона для атомов и молекул. [81]
Стандартный расчет парного взаимодействия в МД приводит к сложности для частиц. Это происходит потому, что каждая частица взаимодействует с каждой другой частицей, что приводит к взаимодействиям. [82] Продвинутые алгоритмы, такие как суммирование Эвальда или быстрый мультипольный метод, сводят это к или даже группируют отдаленные частицы и рассматривают их как единое целое или используют умные математические аппроксимации. [83] [84]
Объединяет квантово-механические расчеты для небольшой области с молекулярной механикой для более крупной среды. [85]
Сложность методов QM/MM зависит как от размера квантовой области, так и от метода, используемого для квантовых вычислений. Например, если для квантовой части используется метод Хартри-Фока, сложность может быть аппроксимирована как , где — число базисных функций в квантовой области. Эта сложность возникает из-за необходимости решать набор связанных уравнений итеративно до тех пор, пока не будет достигнута самосогласованность. [86]
Находит единственное состояние Фока, которое минимизирует энергию. [87]
NP-трудный или NP-полный, как показано путем встраивания примеров модели Изинга в вычисления Хартри-Фока. Метод Хартри-Фока включает решение уравнений Рутана-Холла, которые масштабируются в зависимости от реализации, причем является числом базисных функций. Вычислительные затраты в основном связаны с оценкой и преобразованием двухэлектронных интегралов. Это доказательство NP-трудности или NP-полноты происходит от встраивания таких проблем, как модель Изинга, в формализм Хартри-Фока. [87]
Исследует электронную структуру или ядерную структуру многочастичных систем, таких как атомы, молекулы и конденсированные фазы . [89]
Традиционные реализации DFT обычно масштабируются как , в основном из-за необходимости диагонализации матрицы Кона-Шэма . [90] Шаг диагонализации, который находит собственные значения и собственные векторы матрицы, вносит наибольший вклад в это масштабирование. Последние достижения в DFT направлены на уменьшение этой сложности с помощью различных приближений и алгоритмических улучшений. [91]
Методы CCSD и CCSD(T) представляют собой передовые методы электронной структуры, включающие одиночные, двойные и, в случае CCSD(T), тройные пертурбативные возбуждения для расчета эффектов электронной корреляции. [92]
Масштабируется как , где - число базисных функций. Эта интенсивная вычислительная потребность возникает из-за включения одиночных и двойных возбуждений в расчет электронной корреляции. [92]
С добавлением пертурбативных троек сложность увеличивается до . Эта повышенная сложность ограничивает практическое использование более мелкими системами, обычно до 20-25 атомов в обычных реализациях. [92]
Адаптация стандартного метода CCSD(T) с использованием локальных естественных орбиталей (NO) для значительного снижения вычислительной нагрузки и возможности применения в более крупных системах. [92]
Достигает линейного масштабирования с размером системы, что является значительным улучшением по сравнению с традиционным масштабированием пятой степени CCSD. Это достижение позволяет применять практические приложения к молекулам до 100 атомов с разумными базисными наборами, что является значительным шагом вперед в способности вычислительной химии обрабатывать более крупные системы с высокой точностью. [92]
Доказательство классов сложности для алгоритмов включает комбинацию математического доказательства и вычислительных экспериментов. Например, в случае метода Хартри-Фока доказательство NP-трудности является теоретическим результатом, полученным из теории сложности, в частности, посредством редукций из известных NP-трудных задач. [93]
Для других методов, таких как MD или DFT, вычислительная сложность часто наблюдается эмпирически и поддерживается анализом алгоритма. В этих случаях доказательство корректности в меньшей степени касается формальных математических доказательств, а в большей — последовательного наблюдения за вычислительным поведением в различных системах и реализациях. [93]
Вычислительная химия не является точным описанием реальной химии, поскольку математические и физические модели природы могут дать только приближение. Однако большинство химических явлений можно описать в определенной степени в качественной или приблизительной количественной вычислительной схеме. [94]
Молекулы состоят из ядер и электронов, поэтому применяются методы квантовой механики . Химики-вычислители часто пытаются решить нерелятивистское уравнение Шредингера с добавлением релятивистских поправок, хотя некоторый прогресс был достигнут в решении полностью релятивистского уравнения Дирака . В принципе, можно решить уравнение Шредингера либо в его зависящей от времени, либо в независимой от времени форме, в зависимости от решаемой задачи; на практике это невозможно, за исключением очень малых систем. Поэтому большое количество приближенных методов стремятся достичь наилучшего компромисса между точностью и вычислительными затратами. [95]
Точность всегда можно улучшить за счет больших вычислительных затрат. Значительные ошибки могут возникнуть в моделях ab initio , включающих много электронов, из-за вычислительных затрат полных релятивистских инклюзивных методов. [92] Это усложняет изучение молекул, взаимодействующих с атомами с высокой атомной массой, такими как переходные металлы и их каталитические свойства. Современные алгоритмы в вычислительной химии могут в обычном режиме рассчитывать свойства малых молекул, которые содержат до 40 электронов с ошибками для энергий менее нескольких кДж/моль. Для геометрий длины связей могут быть предсказаны в пределах нескольких пикометров, а углы связей в пределах 0,5 градуса. Обработка более крупных молекул, которые содержат несколько десятков атомов, вычислительно поддается более приближенным методам, таким как теория функционала плотности (DFT). [96]
В этой области ведутся споры о том, достаточны ли последние методы для описания сложных химических реакций, например, в биохимии. Большие молекулы можно изучать полуэмпирическими приближенными методами. Еще более крупные молекулы обрабатываются методами классической механики , которые используют то, что называется молекулярной механикой (ММ). В методах КМ-ММ небольшие части больших комплексов обрабатываются квантово-механически (КМ), а остальная часть обрабатывается приближенно (ММ). [97]
Существует множество самодостаточных пакетов программного обеспечения для вычислительной химии . Некоторые включают множество методов, охватывающих широкий диапазон, в то время как другие концентрируются на очень узком диапазоне или даже на одном методе. Подробную информацию о большинстве из них можно найти в: