The negative sign problem in quantum Monte Carlo (QMC) simulations of cluster impurity problems is the major bottleneck in cluster dynamical mean field calculations. In this paper we systematically investigate the dependence of the sign problem on the single-particle basis. We explore both the hybridization-expansion and the interaction-expansion variants of continuous-time QMC for three-site and four-site impurity models with baths that are diagonal in the orbital degrees of freedom. We find that the sign problem in these models can be substantially reduced by using a non-trivial single-particle basis. Such bases can be generated by diagonalizing a subset of the intracluster hoppings.
We present a continuous-time Monte Carlo method for quantum impurity models, which combines a weak-coupling expansion with an auxiliary-field decomposition. The method is considerably more efficient than Hirsch-Fye and free of time discretization err
ors, and is particularly useful as impurity solver in large cluster dynamical mean field theory (DMFT) calculations.
We describe an open-source implementation of the continuous-time interaction-expansion quantum Monte Carlo method for cluster-type impurity models with onsite Coulomb interactions and complex Weiss functions. The code is based on the ALPS libraries.
We discuss a projector Monte Carlo method for quantum spin models formulated in the valence bond basis, using the S=1/2 Heisenberg antiferromagnet as an example. Its singlet ground state can be projected out of an arbitrary basis state as the trial s
tate, but a more rapid convergence can be obtained using a good variational state. As an alternative to first carrying out a time consuming variational Monte Carlo calculation, we show that a very good trial state can be generated in an iterative fashion in the course of the simulation itself. We also show how the properties of the valence bond basis enable calculations of quantities that are difficult to obtain with the standard basis of Sz eigenstates. In particular, we discuss quantities involving finite-momentum states in the triplet sector, such as the dispersion relation and the spectral weight of the lowest triplet.
In the present paper, we present an efficient continuous-time quantum Monte Carlo impurity solver with high acceptance rate at low temperature for multi-orbital quantum impurity models with general interaction. In this hybridization expansion impurit
y solver, the imaginary time evolution operator for the high energy multiplets, which decays very rapidly with the imaginary time, is approximated by a probability normalized $delta$-function. As the result, the virtual charge fluctuations of $f^{n}rightarrow f^{npm1}$ are well included on the same footing without applying Schrieffer-Wolff transformation explicitly. As benchmarks, our algorithm perfectly reproduces the results for both Coqblin-Schriffeer and Kondo lattice models obtained by CT-J method developed by Otsuki {it et al}. Furthermore, it allows capturing low energy physics of heavy-fermion materials directly without fitting the exchange coupling $J$ in the Kondo model.
We show how efficient loop updates, originally developed for Monte Carlo simulations of quantum spin systems at finite temperature, can be combined with a ground-state projector scheme and variational calculations in the valence bond basis. The metho
ds are formulated in a combined space of spin z-components and valence bonds. Compared to schemes formulated purely in the valence bond basis, the computational effort is reduced from up to O(N^2) to O(N) for variational calculations, where N is the system size, and from O(m^2) to O(m) for projector simulations, where m>> N is the projection power. These improvements enable access to ground states of significantly larger lattices than previously. We demonstrate the efficiency of the approach by calculating the sublattice magnetization M_s of the two-dimensional Heisenberg model to high precision, using systems with up to 256*256 spins. Extrapolating the results to the thermodynamic limit gives M_s=0.30743(1). We also discuss optimized variational amplitude-product states, which were used as trial states in the projector simulations, and compare results of projecting different types of trial states.