Jump to content

Quantum Monte Carlo

From Wikipedia, the free encyclopedia

This is an old revision of this page, as edited by Quantum-Collapse (talk | contribs) at 10:06, 19 January 2021 (→See also). The present address (URL) is a permanent link to this revision, which may differ significantly from the current revision.

Electronic structure methods
Valence bond theory
Coulson–Fischer theory
Generalized valence bond
Modern valence bond theory
Molecular orbital theory
Hartree–Fock method
Semi-empirical quantum chemistry methods
Møller–Plesset perturbation theory
Configuration interaction
Coupled cluster
Multi-configurational self-consistent field
Quantum chemistry composite methods
Quantum Monte Carlo
Density functional theory
Time-dependent density functional theory
Thomas–Fermi model
Orbital-free density functional theory
Linearized augmented-plane-wave method
Projector augmented wave method
Electronic band structure
Nearly free electron model
Tight binding
Muffin-tin approximation
k·p perturbation theory
Empty lattice approximation
GW approximation
Korringa–Kohn–Rostoker method

Quantum Monte Carlo encompasses a large family of computational methods whose common aim is the study of complex quantum systems. One of the major goals of these approaches is to provide a reliable solution (or an accurate approximation) of the quantum many-body problem. The diverse flavor of quantum Monte Carlo approaches all share the common use of the Monte Carlo method to handle the multi-dimensional integrals that arise in the different formulations of the many-body problem. The quantum Monte Carlo methods allow for a direct treatment and description of complex many-body effects encoded in the wave function, going beyond mean-field theory and offering an exact solution of the many-body problem in some circumstances. In particular, there exist numerically exact and polynomially-scaling algorithms to exactly study static properties of boson systems without geometrical frustration. For fermions, there exist very good approximations to their static properties and numerically exact exponentially scaling quantum Monte Carlo algorithms, but none that are both.

Background

In principle, any physical system can be described by the many-body Schrödinger equation as long as the constituent particles are not moving "too" fast; that is, they are not moving at a speed comparable to that of light, and relativistic effects can be neglected. This is true for a wide range of electronic problems in condensed matter physics, in Bose–Einstein condensates and superfluids such as liquid helium. The ability to solve the Schrödinger equation for a given system allows prediction of its behavior, with important applications ranging from materials science to complex biological systems. The difficulty is however that solving the Schrödinger equation requires the knowledge of the many-body wave function in the many-body Hilbert space, which typically has an exponentially large size in the number of particles. Its solution for a reasonably large number of particles is therefore typically impossible, even for modern parallel computing technology in a reasonable amount of time. Traditionally, approximations for the many-body wave function as an antisymmetric function of one-body orbitals[1] have been used, in order to have a manageable treatment of the Schrödinger equation. However, this kind of formulation has several drawbacks, either limiting the effect of quantum many-body correlations, as in the case of the Hartree–Fock (HF) approximation, or converging very slowly, as in configuration interaction applications in quantum chemistry.

Quantum Monte Carlo is a way to directly study the many-body problem and the many-body wave function beyond these approximations. The most advanced quantum Monte Carlo approaches provide an exact solution to the many-body problem for non-frustrated interacting boson systems, while providing an approximate, yet typically very accurate, description of interacting fermion systems. Most methods aim at computing the ground state wavefunction of the system, with the exception of path integral Monte Carlo and finite-temperature auxiliary-field Monte Carlo, which calculate the density matrix. In addition to static properties, the time-dependent Schrödinger equation can also be solved, albeit only approximately, restricting the functional form of the time-evolved wave function, as done in the time-dependent variational Monte Carlo. From the probabilistic point of view, the computation of the top eigenvalues and the corresponding ground states eigenfunctions associated with the Schrödinger equation relies on the numerical solving of Feynman–Kac path integration problems.[2][3] The mathematical foundations of Feynman–Kac particle absorption models and their Sequential Monte Carlo and mean-field interpretations are developed in.[4][5][6][7][8]

There are several quantum Monte Carlo methods, each of which uses Monte Carlo in different ways to solve the many-body problem:

Quantum Monte Carlo methods

Zero-temperature (only ground state)

  • Variational Monte Carlo: A good place to start; it is commonly used in many sorts of quantum problems.
    • Diffusion Monte Carlo: The most common high-accuracy method for electrons (that is, chemical problems), since it comes quite close to the exact ground-state energy fairly efficiently. Also used for simulating the quantum behavior of atoms, etc.
    • Reptation Monte Carlo: Recent zero-temperature method related to path integral Monte Carlo, with applications similar to diffusion Monte Carlo but with some different tradeoffs.
  • Gaussian quantum Monte Carlo
  • Path integral ground state: Mainly used for boson systems; for those it allows calculation of physical observables exactly, i.e. with arbitrary accuracy

Finite-temperature (thermodynamic)

Real-time dynamics (closed quantum systems)

See also

Implementations

Notes

  1. ^ "Functional form of the wave function". Archived from the original on July 18, 2009. Retrieved April 22, 2009.
  2. ^ Caffarel, Michel; Claverie, Pierre (1988). "Development of a pure diffusion quantum Monte Carlo method using a full generalized Feynman–Kac formula. I. Formalism". The Journal of Chemical Physics. 88 (2): 1088–1099. Bibcode:1988JChPh..88.1088C. doi:10.1063/1.454227. ISSN 0021-9606.
  3. ^ Korzeniowski, A.; Fry, J. L.; Orr, D. E.; Fazleev, N. G. (August 10, 1992). "Feynman–Kac path-integral calculation of the ground-state energies of atoms". Physical Review Letters. 69 (6): 893–896. Bibcode:1992PhRvL..69..893K. doi:10.1103/PhysRevLett.69.893. PMID 10047062.
  4. ^ "EUDML | Particle approximations of Lyapunov exponents connected to Schrödinger operators and Feynman–Kac semigroups – P. Del Moral, L. Miclo". eudml.org. Retrieved June 11, 2015.
  5. ^ Del Moral, Pierre; Doucet, Arnaud (January 1, 2004). "Particle Motions in Absorbing Medium with Hard and Soft Obstacles". Stochastic Analysis and Applications. 22 (5): 1175–1207. doi:10.1081/SAP-200026444. ISSN 0736-2994. S2CID 4494495.
  6. ^ Del Moral, Pierre (2013). Mean field simulation for Monte Carlo integration. Chapman & Hall/CRC Press. p. 626. Monographs on Statistics & Applied Probability
  7. ^ Del Moral, Pierre (2004). Feynman–Kac formulae. Genealogical and interacting particle approximations. Probability and its Applications. Springer. p. 575. ISBN 9780387202686. Series: Probability and Applications
  8. ^ Del Moral, Pierre; Miclo, Laurent (2000). "Branching and Interacting Particle Systems Approximations of Feynman–Kac Formulae with Applications to Non-Linear Filtering". In Jacques Azéma; Michel Ledoux; Michel Émery; Marc Yor (eds.). Séminaire de Probabilités XXXIV (PDF). Lecture Notes in Mathematics. Vol. 1729. pp. 1–145. doi:10.1007/bfb0103798. ISBN 978-3-540-67314-9.
  9. ^ Rousseau, V. G. (May 20, 2008). "Stochastic Green function algorithm". Physical Review E. 77 (5): 056705. arXiv:0711.3839. Bibcode:2008PhRvE..77e6705R. doi:10.1103/physreve.77.056705. PMID 18643193. S2CID 2188292.

References

Retrieved from "https://en.wikipedia.org/w/index.php?title=Quantum_Monte_Carlo&oldid=1001359826"