Research overview

Introduction

The term «quantum Monte Carlo» (QMC) refers to a class of numerical algorithms aimed at solving the Schrödinger equation for a quantum-mechanical system in a stochastic manner. Its best-known flavours are the variational and diffusion Monte Carlo methods (VMC and DMC):

• VMC is a simple application of the Monte Carlo integration technique to obtain the expectation value of the Hamiltonian of the system with respect to a given trial wave function. This expectation value is guaranteed to be an upper bound of the ground-state energy of the system, as per the variational principle. In addition, one can optimize the parameters of the wave function in order to minimize the variational energy, thus improving the accurary of the description of the system.
• DMC is a more sophisticated method which uses the imaginary-time-dependent Schrödinger equation to project out the higher-energy components of the trial wave function to give, in principle, the exact ground state energy of the system.

Advantages

The advantages of QMC over other methods are:

• Because of using wave functions, QMC is a truly many-body method. This is in contrast with Density Functional Theory, whose density-based description of the system removes many-body effects. These must be re-introduced approximately by choosing an appropriate exchange-correlation functional, which is arbitrary in any case.
• Because of being a numerical integration technique, QMC can use any wave function to describe the system. In particular, physically-sensible constraints like the Kato cusp conditions can be enforced directly without prohibitively increasing the cost of the calculation. This is in contrast with some quantum Chemistry methods which rely on being able to calculate integrals analytically, and therefore require large expansions of simple terms to achieve a good accuracy.

Thus QMC provides more accuracy than faster methods (e.g. DFT), while it features a more favourable scaling with system size than the more accurate ones (e.g. Coupled Cluster).

Limitations

VMC is strictly driven by the trial wave function, and it is therefore limited by how good a trial wave function we can construct. DMC should in principle be free from this constraint, but due to the fermion sign problem the method needs to be reformulated into a fixed-node (FN) version, which is only exact when the nodes (zero-valued regions) of the trial wave function coincide with those of the exact ground-state wave function.

There are two main lines of research aimed at solving (or alleviating the effects of) the fermion sign problem:

• Understanding the nodal properties of trial wave functions and designing better functional forms.
• Designing new QMC methods.

The former has seen promising, gradual developments in the past couple of years, whereas the second is still to yield a satisfactory solution to the problem.

Introduction

The standard QMC wave function for fermionic systems is the Slater-Jastrow wave function, consisting of a Slater determinant times a positive Jastrow factor,

Ψ(R)=e ^J(R) D(R)

The Slater determinant is constructed using one-particle orbitals, usually obtained from a cheaper method, and accounts for exchange. The Jastrow factor, which usually contains optimizable parameters, depends explicitly on the inter-particle distances so it is capable of describing correlation effects. It is also used to enforce the Kato cusp conditions, which cancel out the divergencies of the Coulomb potential and give a more stable and accurate calculation.

The Jastrow factor gives the functional form a great deal of flexibility, and it can typically recover about 60%–80% of the correlation energy of the system at the VMC level. However it cannot modify the nodes of the wave function, since e ^J ≥0, therefore it does not affect the DMC energy. Note that this statement can be reversed to give an useful conclusion: DMC is equivalent to VMC with a perfect Jastrow.

(empty)

• Optimal sampling
• The zero variance limit

• Phase diagram of the electron-hole system