Квантовый метод Монте-Карло

Квантовые методы Монте-Карло — большая семья методов, для исследования сложных квантовых систем. Одна из главных задач — обеспечить надёжное решение (или достаточно точное приближение) квантовой задачи многих тел[en]. Различные варианты этого метода имеют общую особенность: они используют метод Монте-Карло для вычисления многомерных интегралов, возникающих в различных формулировках задачи многих тел. Квантовые методы Монте-Карло позволяют описывать сложные эффекты многих частиц, зашифрованные в волновой функции, выходя за рамки теории среднего поля и предлагая в некоторых случаях точные решения задачи многих тел. В частности, существует численно точный и полиномиальный масштабируемый алгоритм точного изучения статических свойств системы бозонов без геометрической фрустрации. Для фермионов не известно таких алгоритмов, но существуют отдельно алгоритмы, которые дают очень хорошие приближения их статических свойств, и отдельно квантовые алгоритмы Монте-Карло, которые численно точны, но экспоненциально масштабируемы.

Введение править

В принципе любая физическая система описывается уравнением Шрёдингера для многих частиц, если только частицы не двигаются слишком быстро (то есть чтобы их скорость оставалась малой по сравнению со скоростью света, и релятивистскими эффектами можно было пренебречь). Это требование выполняется для широкого круга электронных задач в физике конденсированных сред, в Бозе-Эйнштейновском конденсате и в сверхтекучих жидкостях вроде жидкого гелия. Умение решать уравнения Шрёдингера для заданной системы позволяет предсказывать её поведение и имеет важные приложения во многих областях науки, начиная с материаловедения и заканчивая сложными биологическими системами. Сложность в том, что решения уравнения Шрёдингера требует знания многочастичной волновой функции в многомерном гильбертовом пространстве, размер которой, как правило, растёт экспоненциально при увеличении числа частиц.

Решение для большого числа частиц в основном невозможно за разумное время, даже для современных параллельных вычислений. Традиционно используются приближения многочастичных антисимметричных функций, составленных из одночастичных молекулярных орбиталей[1], что сводит задачу решения уравнения Шрёдингера к форме, с которой можно работать. Формулировка такого рода имеют несколько недостатков. Они либо ограничиваются учётом квантовых корреляций, например метод Хартри-Фока, или сходятся очень медленно, как в случае применения конфигурационных взаимодействий в квантовой химии.

Квантовые методы Монте-Карло открывает путь к непосредственному изучению многочастичных задач и многочастичных волновых функций без этих ограничений. Наиболее совершенные квантовые методы Монте-Карло дают точные решения многочастичных задачи системы бозонов без фрустраций, одновременно с приближенным, но как правило корректным описанием систем фермионов со взаимодействием. Большинство методов имеют целью нахождение волновой функции основного состояния системы, за исключением методов Монте-Карло для интегралов по траекториям и метода Монте-Карло для конечных температур, которые используются для вычисления матрицы плотности. Кроме стационарных задач можно решать также зависящее от времени уравнение Шрёдингера, хотя лишь приближено, ограничивая функциональную форму зависимой от времени волновой функции. Для этого разработан зависящий от времени вариационный метод Монте-Карло. С точки зрения теории вероятности вычисления ведущих собственных значений и соответствующих им волновых функций основного состояния опирается на численное решение задачи интегралов вдоль траекторий Фейнмана-Кака[2][3]. Математическая база модели поглощения частиц Фейнмана-Кака, секвенционного метода Монте-Карло и интерпретаций среднего поля заложена в работах[4][5][6][7][8].

Существует несколько квантовых методов Монте-Карло, в каждом из них Монте-Карло используется для решения задачи многих тел различными способами.

Методы править

Нулевая температура (только основное состояние) править

  • Вариационный метод Монте-Карло: неплохой исходный пункт; используется при решении широкого круга различных квантовых задач.
  • Диффузный метод Монте-Карло: самый популярный высокоточный метод для системы электронов (то есть, для химических расчётов), поскольку он сравнительно эффективно сходится к точному значению энергии основого состояния. Используется также для воспроизведения квантовой поведения атомов и тому подобное.
  • Рептационный метод Монте-Карло: современный метод вычислений при нулевой температуре, связанный с интегралами по траекториям, область применения та же что и диффузионного метода Монте-Карло, но предположения другие, поэтому преимущества и недостатки отличаются. Рептация — термин из физики полимеров, описывает переползания длинных цепочек змейкой.
  • Гауссов квантовый метод Монте-Карло
  • Нахождение основного состояния через интегралы вдоль траекторий: в основном используется для системы бозонов; для тех где физические наблюдаемые величины можно вычислить точно, то есть с произвольно малой погрешностью.

Ненулевые температуры (термодинамика) править

Динамика реального времени (замкнутые квантовые системы) править

Проекты и программные продукты править

Ссылки править

  1. Functional form of the wave function Архивная копия от 18 июля 2009 на Wayback Machine
  2. 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. — Bibcode1988JChPh..88.1088C. Архивировано 12 июня 2015 года. Архивированная копия. Дата обращения: 18 января 2018. Архивировано 12 июня 2015 года.
  3. 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. — Bibcode1992PhRvL..69..893K.
  4. 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 года.
  5. 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.
  6. Del Moral, Pierre. Mean field simulation for Monte Carlo integration (англ.). — Chapman & Hall/CRC Press, 2013. — P. 626. Архивировано 8 июня 2015 года.. — «Monographs on Statistics & Applied Probability».
  7. Del Moral, Pierre. Feynman-Kac formulae. Genealogical and interacting particle approximations (англ.). — Springer, 2004. — P. 575.. — «Series: Probability and Applications».
  8. 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.
  9. Rousseau, V. G. Stochastic Green function algorithm (англ.) // Physical Review E : journal. — 2008. — 20 May (vol. 77). — P. 056705. — doi:10.1103/physreve.77.056705. — Bibcode2008PhRvE..77e6705R. — arXiv:0711.3839. (недоступная ссылка)