Квантовый метод Монте-Карло
Квантовые методы Монте-Карло — большая семья методов, для исследования сложных квантовых систем. Одна из главных задач — обеспечить надёжное решение (или достаточно точное приближение) квантовой задачи многих тел[англ.]. Различные варианты этого метода имеют общую особенность: они используют метод Монте-Карло для вычисления многомерных интегралов, возникающих в различных формулировках задачи многих тел. Квантовые методы Монте-Карло позволяют описывать сложные эффекты многих частиц, зашифрованные в волновой функции, выходя за рамки теории среднего поля и предлагая в некоторых случаях точные решения задачи многих тел. В частности, существует численно точный и полиномиальный масштабируемый алгоритм точного изучения статических свойств системы бозонов без геометрической фрустрации. Для фермионов не известно таких алгоритмов, но существуют отдельно алгоритмы, которые дают очень хорошие приближения их статических свойств, и отдельно квантовые алгоритмы Монте-Карло, которые численно точны, но экспоненциально масштабируемы.
Введение
[править | править код]В принципе любая физическая система описывается уравнением Шрёдингера для многих частиц, если только частицы не двигаются слишком быстро (то есть чтобы их скорость оставалась малой по сравнению со скоростью света, и релятивистскими эффектами можно было пренебречь). Это требование выполняется для широкого круга электронных задач в физике конденсированных сред, в Бозе-Эйнштейновском конденсате и в сверхтекучих жидкостях вроде жидкого гелия. Умение решать уравнения Шрёдингера для заданной системы позволяет предсказывать её поведение и имеет важные приложения во многих областях науки, начиная с материаловедения и заканчивая сложными биологическими системами. Сложность в том, что решения уравнения Шрёдингера требует знания многочастичной волновой функции в многомерном гильбертовом пространстве, размер которой, как правило, растёт экспоненциально при увеличении числа частиц.
Решение для большого числа частиц в основном невозможно за разумное время, даже для современных параллельных вычислений. Традиционно используются приближения многочастичных антисимметричных функций, составленных из одночастичных молекулярных орбиталей[1], что сводит задачу решения уравнения Шрёдингера к форме, с которой можно работать. Формулировка такого рода имеют несколько недостатков. Они либо ограничиваются учётом квантовых корреляций, например метод Хартри-Фока, или сходятся очень медленно, как в случае применения конфигурационных взаимодействий в квантовой химии.
Квантовые методы Монте-Карло открывает путь к непосредственному изучению многочастичных задач и многочастичных волновых функций без этих ограничений. Наиболее совершенные квантовые методы Монте-Карло дают точные решения многочастичных задачи системы бозонов без фрустраций, одновременно с приближенным, но как правило корректным описанием систем фермионов со взаимодействием. Большинство методов имеют целью нахождение волновой функции основного состояния системы, за исключением методов Монте-Карло для интегралов по траекториям и метода Монте-Карло для конечных температур, которые используются для вычисления матрицы плотности. Кроме стационарных задач можно решать также зависящее от времени уравнение Шрёдингера, хотя лишь приближено, ограничивая функциональную форму зависимой от времени волновой функции. Для этого разработан зависящий от времени вариационный метод Монте-Карло. С точки зрения теории вероятности вычисления ведущих собственных значений и соответствующих им волновых функций основного состояния опирается на численное решение задачи интегралов вдоль траекторий Фейнмана-Кака[2][3]. Математическая база модели поглощения частиц Фейнмана-Кака, секвенционного метода Монте-Карло и интерпретаций среднего поля заложена в работах[4][5][6][7][8].
Существует несколько квантовых методов Монте-Карло, в каждом из них Монте-Карло используется для решения задачи многих тел различными способами.
Методы
[править | править код]Нулевая температура (только основное состояние)
[править | править код]- Вариационный метод Монте-Карло: неплохой исходный пункт; используется при решении широкого круга различных квантовых задач.
- Диффузный метод Монте-Карло: самый популярный высокоточный метод для системы электронов (то есть, для химических расчётов), поскольку он сравнительно эффективно сходится к точному значению энергии основого состояния. Используется также для воспроизведения квантовой поведения атомов и тому подобное.
- Рептационный метод Монте-Карло: современный метод вычислений при нулевой температуре, связанный с интегралами по траекториям, область применения та же что и диффузионного метода Монте-Карло, но предположения другие, поэтому преимущества и недостатки отличаются. Рептация — термин из физики полимеров, описывает переползания длинных цепочек змейкой.
- Гауссов квантовый метод Монте-Карло
- Нахождение основного состояния через интегралы вдоль траекторий: в основном используется для системы бозонов; для тех где физические наблюдаемые величины можно вычислить точно, то есть с произвольно малой погрешностью.
Ненулевые температуры (термодинамика)
[править | править код]- Метод Монте-Карло вспомогательного поля: в основном применяется для задач, определённых на решётке, хотя существуют новые работы, которые применяют этот метод к электронам в химических системах.
- Квантовый метод Монте-Карло с непрерывным временем.
- Детерминантный квантовый метод Монте-Карло или Квантовый метод Монте-Карло Хирша-Фая
- Гибридный квантовый метод Монте-Карло
- Квантовый метод Монте-Карло через интегралы вдоль траекторий: методика вычислений при ненулевых температурах, которую в основном используют для систем, где температурные эффекты имеют большое значение, в частности для сверхтекучего гелия.
- Стохастический алгоритм для функции Грина[9]: алгоритм, сконструированный для бозонов, моделирует определённый на решетке гамильтониан любой сложности, если только в нём нет проблемы со знаком.
- Квантовый метод Монте-Карло мировых линий.
Динамика реального времени (замкнутые квантовые системы)
[править | править код]- Времязависимый вариационный квантовый метод Монте-Карло: расширение вариационного метода Монте-Карло на динамику чистых квантовых состояний.
Проекты и программные продукты
[править | править код]- ALF
- ALPS Архивировано 18 мая 2016 года.
- CASINO Архивная копия от 8 июня 2015 на Wayback Machine
- CHAMP Архивная копия от 6 августа 2017 на Wayback Machine
- Monte Python Архивная копия от 25 октября 2012 на Wayback Machine
- PIMC++
- pi-qmc Архивная копия от 14 июля 2012 на Wayback Machine
- QMcBeaver Архивная копия от 11 июля 2017 на Wayback Machine
- QmcMol
- QMCPACK Архивная копия от 28 ноября 2019 на Wayback Machine
- Qumax
- Qwalk Архивная копия от 4 февраля 2020 на Wayback Machine
- TurboRVB Архивная копия от 18 апреля 2017 на Wayback Machine
- Zori
Ссылки
[править | править код]- ↑ Functional form of the wave function Архивная копия от 18 июля 2009 на Wayback Machine
- ↑ Caffarel, Michel; Claverie, Pierre. Development of a pure diffusion quantum Monte Carlo method using a full generalized Feynman–Kac formula. I. Formalism (англ.) // Journal of Chemical Physics : journal. — 1988. — Vol. 88, no. 2. — P. 1088—1099. — ISSN 0021-9606. — doi:10.1063/1.454227. — . Архивировано 12 июня 2015 года. Архивированная копия . Дата обращения: 18 января 2018. Архивировано 12 июня 2015 года.
- ↑ Korzeniowski, A.; Fry, J. L.; Orr, D. E.; Fazleev, N. G. Feynman-Kac path-integral calculation of the ground-state energies of atoms (англ.) // Physical Review Letters : journal. — 1992. — 10 August (vol. 69, no. 6). — P. 893—896. — doi:10.1103/PhysRevLett.69.893. — .
- ↑ EUDML | Particle approximations of Lyapunov exponents connected to Schrödinger operators and Feynman–Kac semigroups - P. Del Moral, L. Miclo. eudml.org. Дата обращения: 11 июня 2015. Архивировано 4 февраля 2017 года.
- ↑ Del Moral, Pierre; Doucet, Arnaud. Particle Motions in Absorbing Medium with Hard and Soft Obstacles (англ.) // Stochastic Analysis and Applications : journal. — 2004. — 1 January (vol. 22, no. 5). — P. 1175—1207. — ISSN 0736-2994. — doi:10.1081/SAP-200026444.
- ↑ Del Moral, Pierre. Mean field simulation for Monte Carlo integration (англ.). — Chapman & Hall/CRC Press, 2013. — P. 626. Архивировано 8 июня 2015 года.. — «Monographs on Statistics & Applied Probability».
- ↑ Del Moral, Pierre. Feynman-Kac formulae. Genealogical and interacting particle approximations (англ.). — Springer, 2004. — P. 575.. — «Series: Probability and Applications».
- ↑ Del Moral, Pierre; Miclo, Laurent. Branching and Interacting Particle Systems Approximations of Feynman-Kac Formulae with Applications to Non-Linear Filtering (англ.). — 2000. — Vol. 1729. — P. 1—145. — doi:10.1007/bfb0103798.
- ↑ Rousseau, V. G. Stochastic Green function algorithm (англ.) // Physical Review E : journal. — 2008. — 20 May (vol. 77). — P. 056705. — doi:10.1103/physreve.77.056705. — . — arXiv:0711.3839. (недоступная ссылка)