Continuous-time quantum Monte Carlo

Last updated March 07, 2023

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.^[1]^[2]^[3]^[4]^[5] 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.^[1]

If the sign problem is absent, the method can also be used to solve lattice models such as the Hubbard model at half filling. To distinguish it from other Monte Carlo methods for such systems that also work in continuous time, the method is then usually referred to as Diagrammatic determinantal quantum Monte Carlo (DDQMC or DDMC).^[6]

Partition function expansion

In second quantisation, the Hamiltonian of the Anderson impurity model reads:^[1]

H=\underbrace {\sum _{ij}E_{ij}c_{i}^{\dagger }c_{j}} _{H_{\mathrm {loc0} }}+\underbrace {{\frac {1}{2}}\sum _{ijkl}U_{ikjl}c_{i}^{\dagger }c_{j}^{\dagger }c_{l}c_{k}} _{H_{\mathrm {int} }}+\underbrace {\sum _{p,i}{\big (}V_{pi}f_{p}^{\dagger }c_{i}+V_{pi}^{*}c_{i}^{\dagger }f_{p}{\big )}} _{H_{\mathrm {hyb} }}+\underbrace {\sum _{p}\epsilon _{p}f_{p}^{\dagger }f_{p}} _{H_{\mathrm {bath} }}

,

where $c_{i}^{\dagger }$ and $c_{i}$ are the creation and annihilation operators, respectively, of a fermion on the impurity. The index $i$ collects the spin index and possibly other quantum numbers such as orbital (in the case of a multi-orbital impurity) and cluster site (in the case of multi-site impurity). $f_{p}^{\dagger }$ and $f_{p}$ are the corresponding fermion operators on the non-interacting bath, where the bath quantum number $p$ will typically be continuous.

Step 1 of CT-QMC is to split the Hamiltonian into an exactly solvable term, $H_{0}$ , and the rest, $H_{\mathrm {I} }$ . Different choices correspond to different expansions and thus different algorithmic descriptions. Common choices are:

Interaction expansion (CT-INT):^[2] $H_{\mathrm {I} }=H_{\mathrm {int} }$
Hybridization expansion (CT-HYB):^[3]^[4] $H_{\mathrm {I} }=H_{\mathrm {hyb} }$
Auxiliary field expansion (CT-AUX):^[5] like CT-INT, but the interaction term is first decoupled using a discrete Hubbard-Stratonovich transformation

Step 2 is to switch to the interaction picture and expand the partition function in terms of a Dyson series:

Z=\operatorname {tr} \left(\mathrm {e} ^{-\beta H}\right)=\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{n!}}\int _{0}^{\beta }\mathrm {d} ^{n}\tau \;\operatorname {tr} \left[\mathrm {e} ^{-\beta H_{0}}T_{\tau }H_{\mathrm {I} }(\tau _{1})H_{\mathrm {I} }(\tau _{2})\cdots H_{\mathrm {I} }(\tau _{n})\right]

,

where $\beta$ is the inverse temperature and $T_{\tau }$ denotes imaginary time ordering. The presence of a (zero-dimensional) lattice regularises the series and the finite size and temperature of the system makes renormalisation unnecessary.^[2]

The Dyson series generates a factorial number of identical diagrams per order, which makes sampling more difficult and possibly worsen the sign problem. Thus, as step 3, one uses Wick's theorem to group identical diagrams into determinants. This leads to the expressions:^[1]

Interaction expansion (CT-INT):

Z=\sum _{n=0}^{\infty }\prod _{\alpha =1}^{n}\sum _{i_{\alpha }j_{\alpha }k_{\alpha }l_{\alpha }}\int \mathrm {d} \tau _{\alpha }\left(-{\frac {1}{2}}U_{i_{\alpha }k_{\alpha }j_{\alpha }l_{\alpha }}\right)\det \left[{\begin{array}{cc}\langle T_{\tau }c_{i_{\alpha }}^{\dagger }\!(\tau _{\alpha })\ c_{k_{\beta }}\!(\tau _{\beta })\rangle _{0}&\langle T_{\tau }c_{i_{\alpha }}^{\dagger }\!(\tau _{\alpha })\ c_{l_{\beta }}\!(\tau _{\beta })\rangle _{0}\\\langle T_{\tau }c_{j_{\alpha }}^{\dagger }\!(\tau _{\alpha })\ c_{k_{\beta }}\!(\tau _{\beta })\rangle _{0}&\langle T_{\tau }c_{j_{\alpha }}^{\dagger }\!(\tau _{\alpha })\ c_{l_{\beta }}\!(\tau _{\beta })\rangle _{0}\end{array}}\right]_{\alpha \beta }

Hybridisation expansion (CT-HYB):

Z=\sum _{n=0}^{\infty }\prod _{\alpha =1}^{n}\sum _{i_{\alpha }j_{\alpha }}\int \mathrm {d} \tau _{\alpha }\mathrm {d} \tau _{\alpha }^{\prime }\operatorname {tr} {\Big [}\mathrm {e} ^{-\beta (H_{\mathrm {loc0} }+H_{\mathrm {int} })}T_{\tau }\prod _{\alpha =1}^{n}c_{i_{\alpha }}^{\dagger }(\tau _{\alpha })c_{i_{\alpha }}(\tau _{\alpha }^{\prime }){\Big ]}\;\det {\big [}\Delta _{i_{\alpha }j_{\beta }}(\tau _{\alpha }-\tau _{\beta }^{\prime }){\big ]}_{\alpha \beta }

In a final step, one notes that this is nothing but an integral over a large domain and performs it using a Monte Carlo method, usually the Metropolis–Hastings algorithm.

Related Research Articles

In mechanics, the virial theorem provides a general equation that relates the average over time of the total kinetic energy of a stable system of discrete particles, bound by potential forces, with that of the total potential energy of the system. Mathematically, the theorem states

Fractional calculus is a branch of mathematical analysis that studies the several different possibilities of defining real number powers or complex number powers of the differentiation operator

In quantum mechanics, the particle in a one-dimensional lattice is a problem that occurs in the model of a periodic crystal lattice. The potential is caused by ions in the periodic structure of the crystal creating an electromagnetic field so electrons are subject to a regular potential inside the lattice. It is a generalization of the free electron model, which assumes zero potential inside the lattice.

In physics, the S-matrix or scattering matrix relates the initial state and the final state of a physical system undergoing a scattering process. It is used in quantum mechanics, scattering theory and quantum field theory (QFT).

In physics, the Polyakov action is an action of the two-dimensional conformal field theory describing the worldsheet of a string in string theory. It was introduced by Stanley Deser and Bruno Zumino and independently by L. Brink, P. Di Vecchia and P. S. Howe in 1976, and has become associated with Alexander Polyakov after he made use of it in quantizing the string in 1981. The action reads

In thermodynamics, an activity coefficient is a factor used to account for deviation of a mixture of chemical substances from ideal behaviour. In an ideal mixture, the microscopic interactions between each pair of chemical species are the same and, as a result, properties of the mixtures can be expressed directly in terms of simple concentrations or partial pressures of the substances present e.g. Raoult's law. Deviations from ideality are accommodated by modifying the concentration by an activity coefficient. Analogously, expressions involving gases can be adjusted for non-ideality by scaling partial pressures by a fugacity coefficient.

Critical exponents describe the behavior of physical quantities near continuous phase transitions. It is believed, though not proven, that they are universal, i.e. they do not depend on the details of the physical system, but only on some of its general features. For instance, for ferromagnetic systems, the critical exponents depend only on:

In the statistical mechanics of quantum mechanical systems and quantum field theory, the properties of a system in thermal equilibrium can be described by a mathematical object called a Kubo–Martin–Schwinger state or, more commonly, a KMS state: a state satisfying the KMS condition. Kubo (1957) introduced the condition, Martin & Schwinger (1959) used it to define thermodynamic Green's functions, and Rudolf Haag, M. Winnink, and N. M. Hugenholtz (1967) used the condition to define equilibrium states and called it the KMS condition.

Local-density approximations (LDA) are a class of approximations to the exchange–correlation (XC) energy functional in density functional theory (DFT) that depend solely upon the value of the electronic density at each point in space. Many approaches can yield local approximations to the XC energy. However, overwhelmingly successful local approximations are those that have been derived from the homogeneous electron gas (HEG) model. In this regard, LDA is generally synonymous with functionals based on the HEG approximation, which are then applied to realistic systems.

The Sudarshan-Glauber P representation is a suggested way of writing down the phase space distribution of a quantum system in the phase space formulation of quantum mechanics. The P representation is the quasiprobability distribution in which observables are expressed in normal order. In quantum optics, this representation, formally equivalent to several other representations, is sometimes preferred over such alternative representations to describe light in optical phase space, because typical optical observables, such as the particle number operator, are naturally expressed in normal order. It is named after George Sudarshan and Roy J. Glauber, who worked on the topic in 1963. Despite many useful applications in laser theory and coherence theory, the Glauber–Sudarshan P representation has the peculiarity that it is not always positive, and is not a bona-fide probability function.

In many-body theory, the term Green's function is sometimes used interchangeably with correlation function, but refers specifically to correlators of field operators or creation and annihilation operators.

In the context of the physical and mathematical theory of percolation, a percolation transition is characterized by a set of universal critical exponents, which describe the fractal properties of the percolating medium at large scales and sufficiently close to the transition. The exponents are universal in the sense that they only depend on the type of percolation model and on the space dimension. They are expected to not depend on microscopic details such as the lattice structure, or whether site or bond percolation is considered. This article deals with the critical exponents of random percolation.

Dynamical mean-field theory (DMFT) is a method to determine the electronic structure of strongly correlated materials. In such materials, the approximation of independent electrons, which is used in density functional theory and usual band structure calculations, breaks down. Dynamical mean-field theory, a non-perturbative treatment of local interactions between electrons, bridges the gap between the nearly free electron gas limit and the atomic limit of condensed-matter physics.

In statistical mechanics, the Griffiths inequality, sometimes also called Griffiths–Kelly–Sherman inequality or GKS inequality, named after Robert B. Griffiths, is a correlation inequality for ferromagnetic spin systems. Informally, it says that in ferromagnetic spin systems, if the 'a-priori distribution' of the spin is invariant under spin flipping, the correlation of any monomial of the spins is non-negative; and the two point correlation of two monomial of the spins is non-negative.

Heat transfer physics describes the kinetics of energy storage, transport, and energy transformation by principal energy carriers: phonons, electrons, fluid particles, and photons. Heat is energy stored in temperature-dependent motion of particles including electrons, atomic nuclei, individual atoms, and molecules. Heat is transferred to and from matter by the principal energy carriers. The state of energy stored within matter, or transported by the carriers, is described by a combination of classical and quantum statistical mechanics. The energy is different made (converted) among various carriers. The heat transfer processes are governed by the rates at which various related physical phenomena occur, such as the rate of particle collisions in classical mechanics. These various states and kinetics determine the heat transfer, i.e., the net rate of energy storage or transport. Governing these process from the atomic level to macroscale are the laws of thermodynamics, including conservation of energy.

The eigenstate thermalization hypothesis is a set of ideas which purports to explain when and why an isolated quantum mechanical system can be accurately described using equilibrium statistical mechanics. In particular, it is devoted to understanding how systems which are initially prepared in far-from-equilibrium states can evolve in time to a state which appears to be in thermal equilibrium. The phrase "eigenstate thermalization" was first coined by Mark Srednicki in 1994, after similar ideas had been introduced by Josh Deutsch in 1991. The principal philosophy underlying the eigenstate thermalization hypothesis is that instead of explaining the ergodicity of a thermodynamic system through the mechanism of dynamical chaos, as is done in classical mechanics, one should instead examine the properties of matrix elements of observable quantities in individual energy eigenstates of the system.

Tau functions are an important ingredient in the modern theory of integrable systems, and have numerous applications in a variety of other domains. They were originally introduced by Ryogo Hirota in his direct method approach to soliton equations, based on expressing them in an equivalent bilinear form. The term Tau function, or $-function$ , was first used systematically by Mikio Sato and his students in the specific context of the Kadomtsev–Petviashvili equation, and related integrable hierarchies. It is a central ingredient in the theory of solitons. Tau functions also appear as matrix model partition functions in the spectral theory of Random Matrices, and may also serve as generating functions, in the sense of combinatorics and enumerative geometry, especially in relation to moduli spaces of Riemann surfaces, and enumeration of branched coverings, or so-called Hurwitz numbers.

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.

In quantum physics, the "monogamy" of quantum entanglement refers to the fundamental property that it cannot be freely shared between arbitrarily many parties.

Phase space crystal is the state of a physical system that displays discrete symmetry in phase space instead of real space. For a single-particle system, the phase space crystal state refers to the eigenstate of the Hamiltonian for a closed quantum system or the eigenoperator of the Liouvillian for an open quantum system. For a many-body system, phase space crystal is the solid-like crystalline state in phase space. The general framework of phase space crystals is to extend the study of solid state physics and condensed matter physics into phase space of dynamical systems. While real space has Euclidean geometry, phase space is embedded with classical symplectic geometry or quantum noncommutative geometry.

References

1 2 3 4 Gull, E.; Millis, A.J.; Lichtenstein, A.I.; Rubtsov, A.N.; Troyer, M.; Werner, P. (2011). "Continuous-time Monte Carlo methods for quantum impurity models". Rev. Mod. Phys. 83 (2): 349–404. arXiv: 1012.4474 . Bibcode:2011RvMP...83..349G. doi:10.1103/RevModPhys.83.349.
1 2 3 4 Rubtsov, A.N.; Savkin, V.V.; Lichtenstein, A.I. (2005). "Continuous-time quantum Monte Carlo method for fermions". Phys. Rev. B . 72 (3): 035122. arXiv: cond-mat/0411344 . Bibcode:2005PhRvB..72c5122R. doi:10.1103/PhysRevB.72.035122.
1 2 Werner, P.; Comanac, A.; de’ Medici, L.; Troyer, M.; Millis, A.J. (2006). "Continuous-Time Solver for Quantum Impurity Models". Phys. Rev. Lett. 97 (7): 076405. arXiv: cond-mat/0512727 . Bibcode:2006PhRvL..97g6405W. doi:10.1103/PhysRevLett.97.076405.
1 2 Werner, P.; Millis, A.J. (2006). "Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models". Phys. Rev. B. 74 (15): 155107. arXiv: cond-mat/0607136 . Bibcode:2006PhRvB..74o5107W. doi:10.1103/PhysRevB.74.155107.
1 2 Gull, E.; Werner, P.; Parcollet, O.; Troyer, M. (2008). "Continuous-time auxiliary-field Monte Carlo for quantum impurity models". EPL. 82 (5): 57003. arXiv: 0802.3222 . Bibcode:2008EL.....8257003G. doi:10.1209/0295-5075/82/57003.
↑ Assaad, F.F.; Lang, T.C. (2007). "Diagrammatic determinantal quantum Monte Carlo methods: Projective schemes and applications to the Hubbard-Holstein model". Phys. Rev. B . 76 (3): 035116. arXiv: cond-mat/0702455 . Bibcode:2007PhRvB..76c5116A. doi:10.1103/PhysRevB.76.035116.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[gull11-1] 1 2 3 4 Gull, E.; Millis, A.J.; Lichtenstein, A.I.; Rubtsov, A.N.; Troyer, M.; Werner, P. (2011). "Continuous-time Monte Carlo methods for quantum impurity models". Rev. Mod. Phys. 83 (2): 349–404. arXiv: 1012.4474 . Bibcode:2011RvMP...83..349G. doi:10.1103/RevModPhys.83.349.

[rubtsov05-2] 1 2 3 4 Rubtsov, A.N.; Savkin, V.V.; Lichtenstein, A.I. (2005). "Continuous-time quantum Monte Carlo method for fermions". Phys. Rev. B . 72 (3): 035122. arXiv: cond-mat/0411344 . Bibcode:2005PhRvB..72c5122R. doi:10.1103/PhysRevB.72.035122.

[werner06-3] 1 2 Werner, P.; Comanac, A.; de’ Medici, L.; Troyer, M.; Millis, A.J. (2006). "Continuous-Time Solver for Quantum Impurity Models". Phys. Rev. Lett. 97 (7): 076405. arXiv: cond-mat/0512727 . Bibcode:2006PhRvL..97g6405W. doi:10.1103/PhysRevLett.97.076405.

[werner06b-4] 1 2 Werner, P.; Millis, A.J. (2006). "Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models". Phys. Rev. B. 74 (15): 155107. arXiv: cond-mat/0607136 . Bibcode:2006PhRvB..74o5107W. doi:10.1103/PhysRevB.74.155107.

[gull08-5] 1 2 Gull, E.; Werner, P.; Parcollet, O.; Troyer, M. (2008). "Continuous-time auxiliary-field Monte Carlo for quantum impurity models". EPL. 82 (5): 57003. arXiv: 0802.3222 . Bibcode:2008EL.....8257003G. doi:10.1209/0295-5075/82/57003.

[assaad07-6] Assaad, F.F.; Lang, T.C. (2007). "Diagrammatic determinantal quantum Monte Carlo methods: Projective schemes and applications to the Hubbard-Holstein model". Phys. Rev. B . 76 (3): 035116. arXiv: cond-mat/0702455 . Bibcode:2007PhRvB..76c5116A. doi:10.1103/PhysRevB.76.035116.

[1]

[2]

[3]

[4]

[5]

[6]

Continuous-time quantum Monte Carlo

Contents

Partition function expansion

See also

Related Research Articles

References