We address the phenomenon of critical Kondo destruction in pseudogap Bose-Fermi Anderson and Kondo quantum impurity models. These models describe a localized level coupled both to a fermionic bath having a density of states that vanishes like ${|\ensuremath{\epsilon}|}^{r}$ at the Fermi energy ($\ensuremath{\epsilon}=0$) and, via one component of the impurity spin, to a bosonic bath having a sub-Ohmic spectral density proportional to ${|\ensuremath{\omega}|}^{s}$. Each bath is capable by itself of suppressing the Kondo effect at a continuous quantum phase transition. We study the interplay between these two mechanisms for Kondo destruction using continuous-time quantum Monte Carlo for the pseudogap Bose-Fermi Anderson model with
0<r<21 and $\frac{1}{2}\ensuremath{\le}s<1$, and applying the numerical renormalization group to the corresponding Kondo model. At particle-hole symmetry, the models exhibit a quantum-critical point between a Kondo (fermionic strong-coupling) phase and a localized (Kondo-destroyed) phase. The two solution methods, which are in good agreement in their domain of overlap, provide access to the many-body spectrum, as well as to correlation functions including, in particular, the single-particle Green's function and the static and dynamical local spin susceptibilities. The quantum-critical regime exhibits the hyperscaling of critical exponents and $\ensuremath{\omega}/T$ scaling in the dynamics that characterize an interacting critical point. The
(r,s) plane can be divided into three regions: one each in which the calculated critical properties are dominated by the bosonic bath alone or by the fermionic bath alone, and between these two regions, a third in which the bosonic bath governs the critical spin response but both baths influence the renormalization-group flow near the quantum-critical point.
In quantum critical heavy fermion systems, local moments are coupled to both collective spin fluctuations and conduction electrons. As such, the Bose-Fermi Kondo model, describing the coupling of a local moment to both a bosonic and a fermionic bath, has been of extensive interest. For the model in the presence of SU(2) spin rotational symmetry, questions have been raised about its phase diagram. Here we develop a version of continuous-time Quantum Monte Carlo (CT-QMC) method suitable for addressing this issue; this procedure can reach sufficiently low temperatures while preserving the SU(2) symmetry. Using this method for the Bose-Fermi Anderson model, we clarify the renormalization-group fixed points and the phase diagram for the case with a constant fermionic-bath density of states and a power-law bosonic-bath spectral function
ρb(ω)∝ωs (
0<s<1). We find two types of Kondo destruction QCP, depending on the power-law exponent
s in the bosonic bath spectrum. For
s∗<s<1, both types of QCPs exist and, in the parameter regime accessible by an analytical
ϵ-expansion renormalization-group calculation (here
ϵ=1−s), the CT-QMC result is fully consistent with prior predictions by the latter method. For
s<s∗, there is only one type of QCP. At both type of Kondo destruction QCPs, we find that the exponent of the local spin susceptibility
η obeys the relation
η=ϵ, which has important implications for Kondo destruction QCP in the Kondo lattice problem.
Efficient continuous-time quantum Monte Carlo (CT-QMC) algorithms that do not suffer from time discretization errors have become the state of the art for most discrete quantum models. They have not been widely used yet for fermionic quantum lattice models, such as the Hubbard model, nor other fermionic lattice systems due to a suboptimal scaling of $O({\ensuremath{\beta}}^{3})$ with inverse temperature $\ensuremath{\beta}$, compared to the linear scaling of discrete-time algorithms. Here we present a CT-QMC algorithm for fermionic lattice systems that matches the scaling of discrete-time methods but is more efficient and free of time discretization errors. This provides an efficient simulation scheme that is free from the systematic errors opening an avenue to more precise studies of large systems at low and zero temperature.
A vector bosonic field coupled to the electronic spin is treated by means of the continuous-time quantum Monte Carlo method. In the Bose-Kondo model with a sub-Ohmic density of states ${\ensuremath{\rho}}_{\mathrm{B}}(\ensuremath{\omega})\ensuremath{\propto}{\ensuremath{\omega}}^{s}$ with
s=0.2, two contributions to the spin susceptibility, the Curie term ${T}^{\ensuremath{-}1}$ and the term ${T}^{\ensuremath{-}s}$ due to bosonic fluctuations, are observed separately. This result indicates the existence of a residual moment and a hidden critical behavior. By including hybridization with itinerant electrons, a quantum critical point is identified between this local-moment state and the Kondo singlet state. It is demonstrated that the energy scale of the bosonic fluctuations is not affected by the quantum phase transition.
We present a quantum Monte Carlo (QMC) technique for calculating the exact finite-temperature properties of Bose-Fermi mixtures. The Bose-Fermi auxiliary-field quantum Monte Carlo (BFAFQMC) algorithm combines two methods, a finite-temperature AFQMC algorithm for bosons and a variant of the standard AFQMC algorithm for fermions, into one algorithm for mixtures. We demonstrate the accuracy of our method by comparing its results for the Bose-Hubbard and Bose-Fermi-Hubbard models against those produced using exact diagonalization for small systems. Comparisons are also made with mean-field theory and the worm algorithm for larger systems. As is the case with most fermion Hamiltonians, a sign or phase problem is present in the BFAFQMC algorithm. We discuss the nature of these problems in this framework and describe how they can be controlled with well-studied approximations to expand the BFAFQMC algorithm's reach. This algorithm can serve as an essential tool for answering many unresolved questions about many-body physics in mixed Bose-Fermi systems.