Reptation Monte Carlo

Last updated

Reptation Monte Carlo is a quantum Monte Carlo method.

It is similar to Diffusion Monte Carlo, except that it works with paths rather than points. This has some advantages relating to calculating certain properties of the system under study that diffusion Monte Carlo has difficulty with.

In both diffusion Monte Carlo and reptation Monte Carlo, the method first aims to solve the time-dependent Schrödinger equation in the imaginary time direction. When you propagate the Schrödinger equation in time, you get the dynamics of the system under study. When you propagate it in imaginary time, you get a system that tends towards the ground state of the system.

When substituting in place of , the Schrodinger equation becomes identical with a diffusion equation. Diffusion equations can be solved by imagining a huge population of particles (sometimes called "walkers"), each diffusing in a way that solves the original equation. This is how diffusion Monte Carlo works.

Reptation Monte Carlo works in a very similar way, but is focused on the paths that the walkers take, rather than the density of walkers.

In particular, a path may be mutated using a Metropolis algorithm which tries a change (normally at one end of the path) and then accepts or rejects the change based on a probability calculation.

The update step in diffusion Monte Carlo would be moving the walkers slightly, and then duplicating and removing some of them. By contrast, the update step in reptation Monte Carlo mutates a path, and then accepts or rejects the mutation.

Related Research Articles

Lattice gauge theory

In physics, lattice gauge theory is the study of gauge theories on a spacetime that has been discretized into a lattice.

Lattice QCD

Lattice QCD is a well-established non-perturbative approach to solving the quantum chromodynamics (QCD) theory of quarks and gluons. It is a lattice gauge theory formulated on a grid or lattice of points in space and time. When the size of the lattice is taken infinitely large and its sites infinitesimally close to each other, the continuum QCD is recovered.

Numerical relativity is one of the branches of general relativity that uses numerical methods and algorithms to solve and analyze problems. To this end, supercomputers are often employed to study black holes, gravitational waves, neutron stars and many other phenomena governed by Einstein's theory of general relativity. A currently active field of research in numerical relativity is the simulation of relativistic binaries and their associated gravitational waves.

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 of the quantum many-body problem. The diverse flavors 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.

Car–Parrinello molecular dynamics or CPMD refers to either a method used in molecular dynamics or the computational chemistry software package used to implement this method.

In computational physics, variational Monte Carlo (VMC) is a quantum Monte Carlo method that applies the variational method to approximate the ground state of a quantum system.

Diffusion Monte Carlo (DMC) or diffusion quantum Monte Carlo 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 scales exponentially with the system size. This makes exact large-scale DMC simulations for fermions impossible; however, DMC employing a clever approximation known as the fixed-node approximation can still yield very accurate results.

In physics, a quantum amplifier is an amplifier that uses quantum mechanical methods to amplify a signal; examples include the active elements of lasers and optical amplifiers.

In applied mathematics, the numerical sign problem is the problem of numerically evaluating the integral of a highly oscillatory function of a large number of variables. Numerical methods fail because of the near-cancellation of the positive and negative contributions to the integral. Each has to be integrated to very high precision in order for their difference to be obtained with useful accuracy.

The dihydrogen cation or hydrogen molecular ion is a cation with formula H+
2
. It consists of two hydrogen nuclei (protons) sharing a single electron. It is the simplest molecular ion.

In condensed matter physics, biexcitons are created from two free excitons.

A trion is a localized excitation which consists of three charged particles. A negative trion consists of two electrons and one hole and a positive trion consists of two holes and one electron. The trion itself is a quasiparticle and is somewhat similar to an exciton, which is a complex of one electron and one hole. The trion has a ground singlet state and an excited triplet state. Here singlet and triplet degeneracies originate not from the whole system but from the two identical particles in it. The half-integer spin value distinguishes trions from excitons in many phenomena; for example, energy states of trions, but not excitons, are split in an applied magnetic field. Trion states were predicted theoretically in 1958; they were observed experimentally in 1993 in CdTe/Cd1−xZnxTe quantum wells, and later in various other optically excited semiconductor structures. There are experimental proofs of their existence in nanotubes supported by theoretical studies. Despite numerous reports of experimental trion observations in different semiconductor heterostructures, there are serious concerns on the exact physical nature of the detected complexes. The originally foreseen 'true' trion particle has a delocalized wavefunction while recent studies reveal significant binding from charged impurities in real semiconductor quantum wells.

The following timeline starts with the invention of the modern computer in the late interwar period.

The quantum jump method, also known as the Monte Carlo wave function (MCWF) is a technique in computational physics used for simulating open quantum systems and quantum dissipation. The quantum jump method was developed by Dalibard, Castin and Mølmer at a similar time to the similar method known as Quantum Trajectory Theory developed by Carmichael. Other contemporaneous works on wave-function-based Monte Carlo approaches to open quantum systems include those of Dum, Zoller and Ritsch and Hegerfeldt and Wilser.

Ali Alavi

Ali Alavi FRS is a professor of theoretical chemistry in the Department of Chemistry at the University of Cambridge and a Director of the Max Planck Institute for Solid State Research in Stuttgart.

In computational solid state physics, Continuous-time quantum Monte Carlo (CT-QMC) is a family of stochastic algorithms for solving the Anderson impurity model at finite temperature. These methods first expand the full partition function as a series of Feynman diagrams, employ Wick's theorem to group diagrams into determinants, and finally use Markov chain Monte Carlo to stochastically sum up the resulting series.

Ab initio methods (nuclear physics)

In nuclear physics, ab initio methods seek to describe the atomic nucleus from the bottom up by solving the non-relativistic Schrödinger equation for all constituent nucleons and the forces between them. This is done either exactly for very light nuclei or by employing certain well-controlled approximations for heavier nuclei. Ab initio methods constitute a more fundamental approach compared to e.g. the nuclear shell model. Recent progress has enabled ab initio treatment of heavier nuclei such as nickel.

Applying classical methods of machine learning to the study of quantum systems is the focus of an emergent area of physics research. A basic example of this is quantum state tomography, where a quantum state is learned from measurement. Other examples include learning Hamiltonians, learning quantum phase transitions, and automatically generating new quantum experiments. Classical machine learning is effective at processing large amounts of experimental or calculated data in order to characterize an unknown quantum system, making its application useful in contexts including quantum information theory, quantum technologies development, and computational materials design. In this context, it can be used for example as a tool to interpolate pre-calculated interatomic potentials or directly solving the Schrödinger equation with a variational method.

Quantum Trajectory Theory (QTT) is a formulation of quantum mechanics used for simulating open quantum systems, quantum dissipation and single quantum systems. It was developed by Howard Carmichael in the early 1990s around the same time as the similar formulation, known as the quantum jump method or Monte Carlo wave function (MCWF) method, developed by Dalibard, Castin and Mølmer. Other contemporaneous works on wave-function-based Monte Carlo approaches to open quantum systems include those of Dum, Zoller and Ritsch, and Hegerfeldt and Wilser.

In many-body physics, the problem of analytic continuation is that of numerically extracting the spectral density of a Green function given its values on the imaginary axis. It is a necessary post-processing step for calculating dynamical properties of physical systems from quantum Monte Carlo simulations, which often compute Green function values only at imaginary-times or Matsubara frequencies.

References

S. Baroni & S. Moroni (1999). "Reptation Quantum Monte Carlo: A Method for Unbiased Ground-State Averages and Imaginary-Time Correlations". Phys. Rev. Lett. 82 (24): 4745–4748. Bibcode:1999PhRvL..82.4745B. doi:10.1103/PhysRevLett.82.4745.

S. Baroni & S. Moroni (1998). "Reptation Quantum Monte Carlo". arXiv: cond-mat/9808213 .

G. Carleo; F. Becca; S. Moroni & S. Baroni (2010). "Reptation quantum Monte Carlo algorithm for lattice Hamiltonians with a directed-update scheme". Phys. Rev. E. 82 (4): 046710. arXiv: 1003.3696 . Bibcode:2010PhRvE..82d6710C. doi:10.1103/PhysRevE.82.046710. PMID   21230415. S2CID   23090095.