Continuous-time quantum Monte Carlo explained

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.

The attribute continuous-time was introduced to distinguish the method from the then-predominant Hirsch–Fye quantum Monte Carlo method, which relies on a Suzuki–Trotter discretisation of the imaginary time axis.

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).

Partition function expansion

In second quantisation, the Hamiltonian of the Anderson impurity model reads:

H = \underbrace{\sum _{ij} E _{ij} c ^{\dagger} _{i} c _{j}} ^{H _{loc0}} + \underbrace{ (12)/(\sum) _{ijkl} U _{ikjl} c ^{\dagger} _{i} c ^{\dagger} _{j} c _{l} c _{k}} ^{H _{int}} + \underbrace{\sum _{p,i} ( V _{pi} f ^{\dagger} _{p} c _{i} + V ^{*} _{pi} c ^{\dagger} _{i} f _{p} ) } _{H _{hyb}} + \underbrace{\sum _{p} \epsilon _{p} f ^{\dagger} _{p} f _{p}} ^{H _{bath}} ,where c ^{\dagger} _{i} 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 ^{\dagger} _{p} 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 _{I} . Different choices correspond to different expansions and thus different algorithmic descriptions. Common choices are:

• Interaction expansion (CT-INT):

H _{I} = H _{int}

• Hybridization expansion (CT-HYB):

H _{I} = H _{hyb}

• Auxiliary field expansion (CT-AUX): 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( e ^{-\beta} \right) = \sum ^{inf ty} _{n =0} (( -1 ) ^{n})/(n!) \int ^{\beta} _{0} d ^{n\tau} \operatorname{tr}\left[ e ^{- \beta H _{0}} T _{\tau} H _{I(\tau} _{1)} H _{I(\tau} _{2)} … H _{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.

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:

• Interaction expansion (CT-INT):

Z = \sum ^{inf ty} _{n =0} \prod ^{n} _{\alpha =1} \sum _{i _{\alpha} j _{\alpha} k _{\alpha} l _{\alpha}} \int d\tau _{\alpha} \left( - (12)/(U) ^{i _{\alpha} k _{\alpha} j _{\alpha} l _{\alpha}} \right) \det \left[\begin{array}{cc} \langle T _{\tau} c ^{\dagger} _{i _{\alpha}} (\tau _{\alpha) c} ^{k _{\beta}} (\tau _{\beta)} \rangle _{0} & \langle T _{\tau} c ^{\dagger} _{i _{\alpha}} (\tau _{\alpha) c} ^{l _{\beta}} (\tau _{\beta)} \rangle _{0} \\ \langle T _{\tau} c ^{\dagger} _{j _{\alpha}} (\tau _{\alpha) c} ^{k _{\beta}} (\tau _{\beta)} \rangle _{0} & \langle T _{\tau} c ^{\dagger} _{j _{\alpha}} (\tau _{\alpha) c} ^{l _{\beta}} (\tau _{\beta)} \rangle _{0} \end{array} \right] _{\alpha\beta}

• Hybridisation expansion (CT-HYB):

Z = \sum ^{inf ty} _{n =0} \prod ^{n} _{\alpha =1} \sum _{i _{\alpha} j _{\alpha}} \int d\tau _{\alpha} d\tau ^{\prime} _{\alpha} \operatorname{tr} [ e ^{- \beta (H _{loc0} + H _{int)}} T _{\tau} \prod ^{n} _{\alpha =1} c ^{\dagger} _{i _{\alpha}} (\tau _{\alpha)} c _{i _{\alpha}} (\tau ^{\prime} _{\alpha)} ] \det [ \Delta _{i _{\alpha} j _{\beta}} (\tau _{\alpha} - \tau ^{\prime} _{\beta) ]} _{\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.

See also

• Metropolis–Hastings algorithm
• Quantum Monte Carlo
• Dynamical mean field theory

References

^{[1] [2] [3] [4] [5] [6]}

Notes and References

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