Diffusion Monte Carlo


Diffusion Monte Carlo

Diffusion Monte Carlo (DMC) is a quantum Monte Carlo method that uses a Green's function to solve the Schrödinger equation. DMC is potentially numerically exact, meaning that it can find the exact ground state energy within a given error for any quantum system. When actually attempting the calculation, one finds that for bosons, the algorithm scales as a polynomial with the system size, but for fermions, DMC is exponentially scaling with the system size. This makes exact large-scale DMC simulations for fermions impossible; however, with a clever approximation known as fixed-node, very accurate results can be obtained. What follows is an explanation of the basic algorithm, how it works, why fermions cause a problem, and how the fixed-node approximation resolves this problem.

The Projector Method

To motivate the algorithm, let's look at the Schrödinger equation for a particle in some potential in one dimension:

i\frac{d\Psi(x,t)}{dt}=-\frac{1}{2}\frac{d^2 \Psi(x,t)}{dx^2} + V(x)\Psi(x,t).

We can condense the notation a bit by writing it in terms of an operator equation, with

H=-\frac{1}{2}\frac{d^2 }{dx^2} + V(x).

So then we have

i\frac{d\Psi(x,t)}{dt}=H\Psi(x,t),

where we have to keep in mind that H is an operator, not a simple number or function. There are special functions, called eigenfunctions, for which HΨ(x) = EΨ(x), where E is a number. These functions are special because no matter where we evaluate the action of the H operator on the wave function, we always get the same number E. These functions are called stationary states, because the time derivative at any point x is always the same, so the amplitude of the wave function never changes in time. Since the overall phase of a wave function is not measurable, the system does not change in time.

We are usually interested in the wave function with the lowest energy eigenvalue, the ground state. We're going to write a slightly different version of the Schrödinger equation that will have the same energy eigenvalue, but, instead of being oscillatory, it will be convergent. Here it is:

-\frac{d\Psi(x,t)}{dt}=(H-E_0)\Psi(x,t).

We've removed the imaginary number from the time derivative and added in a constant offset of E0, which is the ground state energy. We don't actually know the ground state energy, but there will be a way to determine it self-consistently which we'll introduce later. Our modified equation(some people call it the imaginary-time Schrödinger equation) has some nice properties. The first thing to notice is that if we happen to guess the ground state wave function, then HΦ0(x) = E0Φ0(x) and the time derivative is zero. Now suppose that we start with another wave function(Ψ), which is not the ground state but is not orthogonal to it. Then we can write it as a linear sum of eigenfunctions:

\Psi=c_0\Phi_0+\sum_{i=1}^\infty c_i\Phi_i

Since this is a linear differential equation, we can look at the action of each part separately. We already determined that Φ0 is stationary. Suppose we take Φ1. Since Φ0 is the lowest-energy eigenfunction, the associate eigenvalue of Φ1 satisfies the property E1 > E0. Thus the time derivative of c1 is negative, and will eventually go to zero, leaving us with only the ground state. This observation also gives us a way to determine E0. We watch the amplitude of the wave function as we propagate through time. If it increases, then decrease the estimation of the offset energy. If the amplitude decreases, then increase the estimate of the offset energy.

Stochastic Implementation

Now we have an equation that, as we propagate it forward in time and adjust E0 appropriately, we find the ground state of any given Hamiltonian. This is still a harder problem than classical mechanics, though, because instead of propagating single positions of particles, we must propagate entire functions. In classical mechanics, we could simulate the motion of the particles by setting x(t + τ) = x(t) + τv(t) + 0.5F(t2, if we assume that the force is constant over the time span of τ. For the imaginary time Schrödinger equation, instead, we propagate forward in time using a convolution integral with a special function called a Green's function. So we get  \Psi(x,t+\tau)=\int G(x,x',\tau) \Psi(x',t) dx' . Similarly to classical mechanics, we can only propagate for small slices of time; otherwise the Green's function is inaccurate. As the number of particles increases, the dimensionality of the integral increases as well, since we have to integrate over all coordinates of all particles. We can do these integrals by Monte Carlo integration.

References

  • [1] R.C. Grimm and R.G. Storer, J. Comput. Phys. 7, 134 (1971)
  • [2] J. Anderson, J. Chem. Phys. 63, 1499 (1975)
  • [3] B.L. Hammond, W.A Lester, Jr. & P.J. Reynolds "Monte Carlo Methods in Ab Initio Quantum Chemistry" (World Scientific, 1994)s by Monte Carlo.

Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Monte Carlo TMC — TMC Monte Carlo Création 19 novembre 1954 Slogan « Mais c est bien Sud ! » (1988 1991) « Votre famille du Sud » (1993 1996) « Le meilleur pour vous ! » Langue Français Pays d origine …   Wikipédia en Français

  • Monte Carlo methods for electron transport — The Monte Carlo method for electron transport is a semiclassical Monte Carlo(MC) approach of modeling semiconductor transport. Assuming the carrier motion consists of free flights interrupted by scattering mechanisms, a computer is utilized to… …   Wikipedia

  • Monte Carlo method — Not to be confused with Monte Carlo algorithm. Computational physics …   Wikipedia

  • Monte Carlo method for photon transport — Modeling photon propagation with Monte Carlo methods is a flexible yet rigorous approach to simulate photon transport. In the method, local rules of photon transport are expressed as probability distributions which describe the step size of… …   Wikipedia

  • Monte Carlo methods for option pricing — In mathematical finance, a Monte Carlo option model uses Monte Carlo methods to calculate the value of an option with multiple sources of uncertainty or with complicated features. [1] The term Monte Carlo method was coined by Stanislaw Ulam in… …   Wikipedia

  • Monte Carlo Doualiya — Création 1972 Propriétaire RFI Slogan « La radio qui va changer la radio » Langue Arabe Ancien nom Radio Monte Carlo Moyen Orient RMC Moyen Orient Siège social …   Wikipédia en Français

  • Monte carlo cinétique — Méthode de Monte Carlo cinétique La méthode de Monte Carlo cinétique, kinetic Monte Carlo (KMC) en anglais, est une méthode de Monte Carlo de simulation informatique permettant de simuler des processus se produisant à des taux connus. En cela… …   Wikipédia en Français

  • Quantum Monte Carlo — is a large class of computer algorithms that simulate quantum systems with the idea of solving the many body problem. They use, in one way or another, the Monte Carlo method to handle the many dimensional integrals that arise. Quantum Monte Carlo …   Wikipedia

  • Variational Monte Carlo — Variational Monte Carlo(VMC) is a quantum Monte Carlo method that applies the variational method to approximate the ground state of the system. The expectation value necessary can be written in the x representation as frac{langle Psi(a) | H |… …   Wikipedia

  • TMC MONTE CARLO — Création 19 novembre 1954 Slogan « Mais c est bien Sud ! » (1988 1991) « Votre famille du Sud » (1993 1996) « Le meilleur pour vous ! » Langue Français Pays d origine …   Wikipédia en Français