Bold diagrammatic Monte Carlo method applied to fermionized frustrated spins

Pseudo-fermion functional renormalization group for spin models

For decades, frustrated quantum magnets have been a seed for scientific progress and innovation in condensed matter. As much as the numerical tools for low-dimensional quantum magnetism have thrived and improved in recent years due to breakthroughs inspired by quantum information and quantum computation, higher-dimensional quantum magnetism can be considered as the final frontier, where strong quantum entanglement, multiple ordering channels, and manifold ways of paramagnetism culminate. At the same time, efforts in crystal synthesis have induced a significant increase in the number of tangible frustrated magnets which are generically three-dimensional in nature, creating an urgent need for quantitative theoretical modeling. We review the pseudo-fermion (PF) and pseudo-Majorana (PM) functional renormalization group (FRG) and their specific ability to address higher-dimensional frustrated quantum magnetism. First developed more than a decade ago, the PFFRG interprets a Heisenberg model Hamiltonian in terms of Abrikosov pseudofermions, which is then treated in a diagrammatic resummation scheme formulated as a renormalization group flow ofm-particle pseudofermion vertices. The article reviews the state of the art of PFFRG and PMFRG and discusses their application to exemplary domains of frustrated magnetism, but most importantly, it makes the algorithmic and implementation details of these methods accessible to everyone. By thus lowering the entry barrier to their application, we hope that this review will contribute towards establishing PFFRG and PMFRG as the numerical methods for addressing frustrated quantum magnetism in higher spatial dimensions.

On the magnetization of the 120° order of the spin-1/2 triangular lattice Heisenberg model: a DMRG revisited

We revisit the issue about the magnetization of the 120° order in the spin-1/2 triangular lattice Heisenberg model with density matrix renormalization group (DMRG). The accurate determination of the magnetization of this model is challenging for numerical methods and its value exhibits substantial disparities across various methods. We perform a large-scale DMRG calculation of this model by employing bond dimension as large asD=24000and by studying the system with width as large asLy=12. With careful extrapolation with truncation error and suitable finite size scaling, we give a conservative estimation of the magnetization asM0=0.208(8). The ground state energy per site we obtain isEg=-0.5503(8). Our results provide valuable benchmark values for the development of new methods in the future.

Diagrammatic Monte Carlo Approach to Angular Momentum in Quantum Many-Particle Systems

We introduce a diagrammatic Monte Carlo approach to angular momentum properties of quantum many-particle systems possessing a macroscopic number of degrees of freedom. The treatment is based on a diagrammatic expansion that merges the usual Feynman diagrams with the angular momentum diagrams known from atomic and nuclear structure theory, thereby incorporating the non-Abelian algebra inherent to quantum rotations. Our approach is applicable at arbitrary coupling, is free of systematic errors and of finite-size effects, and naturally provides access to the impurity Green function. We exemplify the technique by obtaining an all-coupling solution of the angulon model; however, the method is quite general and can be applied to a broad variety of systems in which particles exchange quantum angular momentum with their many-body environment.

Resummation of Diagrammatic Series with Zero Convergence Radius for Strongly Correlated Fermions

We demonstrate that a summing up series of Feynman diagrams can yield unbiased accurate results for strongly correlated fermions even when the convergence radius vanishes. We consider the unitary Fermi gas, a model of nonrelativistic fermions in three-dimensional continuous space. Diagrams are built from partially dressed or fully dressed propagators of single particles and pairs. The series is resummed by a conformal-Borel transformation that incorporates the large-order behavior and the analytic structure in the Borel plane, which are found by the instanton approach. We report highly accurate numerical results for the equation of state in the normal unpolarized regime, and reconcile experimental data with the theoretically conjectured fourth virial coefficient.

Spin-charge transformation of lattice fermion models: duality approach for diagrammatic simulation of strongly correlated systems

I derive a dual description of lattice fermions, specifically focusing on the t-J and Hubbard models, that allow diagrammatic techniques to be employed efficiently in the strongly correlated regime, as well as for systems with a restricted Hilbert space. These constructions are based on spin-charge transformation, where the lattice fermions of the original model are mapped onto spins and spin-less fermions. This mapping can then be combined with Popov-Fedotov fermionisation, where the spins are mapped onto lattice fermions with imaginary chemical potential. The resulting models do not contain any large expansion parameters, even for strongly correlated systems. Also, they exhibit dramatically smaller corrections to the density matrix from nonlinear terms in the Hamiltonian. The combination of these two properties means that they can be addressed with diagrammatic methods, including simulation techniques based on stochastic sampling of diagrammatic expansions.

Auxiliary-Field Monte Carlo Method to Tackle Strong Interactions and Frustration in Lattice Bosons

We introduce a new numerical technique, the bosonic auxiliary-field Monte Carlo method, which allows us to calculate the thermal properties of large lattice-boson systems within a systematically improvable semiclassical approach, and which is virtually applicable to any bosonic model. Our method amounts to a decomposition of the lattice into clusters, and to an ansatz for the density matrix of the system in the form of a cluster-separable state-with nonentangled, yet classically correlated clusters. This approximation eliminates any sign problem, and can be systematically improved upon by using clusters of growing size. Extrapolation in the cluster size allows us to reproduce numerically exact results for the superfluid transition of hard-core bosons on the square lattice, and to provide a solid quantitative prediction for the superfluid and chiral transition of hardcore bosons on the frustrated triangular lattice.

Determinant Diagrammatic Monte Carlo Algorithm in the Thermodynamic Limit

We present a simple trick that allows us to consider the sum of all connected Feynman diagrams at fixed position of interaction vertices for general fermionic models, such that the thermodynamic limit can be taken analytically. With our approach one can achieve superior performance compared to conventional diagrammatic Monte Carlo algorithm, while rendering the algorithmic part dramatically simpler. By considering the sum of all connected diagrams at once, we allow for massive cancellations between different diagrams, greatly reducing the sign problem. In the end, the computational effort increases only exponentially with the order of the expansion, which should be contrasted with the factorial growth of the standard diagrammatic technique. We illustrate the efficiency of the technique for the two-dimensional Fermi-Hubbard model.

Stability of Dirac Liquids with Strong Coulomb Interaction

We develop and apply the diagrammatic Monte Carlo technique to address the problem of the stability of the Dirac liquid state (in a graphene-type system) against the strong long-range part of the Coulomb interaction. So far, all attempts to deal with this problem in the field-theoretical framework were limited either to perturbative or random phase approximation and functional renormalization group treatments, with diametrically opposite conclusions. Our calculations aim at the approximation-free solution with controlled accuracy by computing vertex corrections from higher-order skeleton diagrams and establishing the renormalization group flow of the effective Coulomb coupling constant. We unambiguously show that with increasing the system size L (up to ln(L)∼40), the coupling constant always flows towards zero; i.e., the two-dimensional Dirac liquid is an asymptotically free T=0 state with divergent Fermi velocity.

Spin-Ice State of the Quantum Heisenberg Antiferromagnet on the Pyrochlore Lattice

We study the low-temperature physics of the SU(2)-symmetric spin-1/2 Heisenberg antiferromagnet on a pyrochlore lattice and find "fingerprint" evidence for the thermal spin-ice state in this frustrated quantum magnet. Our conclusions are based on the results of bold diagrammatic Monte Carlo simulations, with good convergence of the skeleton series down to the temperature T/J=1/6. The identification of the spin-ice state is done through a remarkably accurate microscopic correspondence for the static structure factor between the quantum Heisenberg, classical Heisenberg, and Ising models at all accessible temperatures, and the characteristic bowtie pattern with pinch points observed at T/J=1/6. The dynamic structure factor at real frequencies (obtained by the analytic continuation of numerical data) is consistent with diffusive spinon dynamics at the pinch points.

Nonexistence of the Luttinger-Ward functional and misleading convergence of skeleton diagrammatic series for hubbard-like models

The Luttinger-Ward functional Φ[G], which expresses the thermodynamic grand potential in terms of the interacting single-particle Green's function G, is found to be ill defined for fermionic models with the Hubbard on-site interaction. In particular, we show that the self-energy Σ[G]∝δΦ[G]/δG is not a single-valued functional of G: in addition to the physical solution for Σ[G], there exists at least one qualitatively distinct unphysical branch. This result is demonstrated for several models: the Hubbard atom, the Anderson impurity model, and the full two-dimensional Hubbard model. Despite this pathology, the skeleton Feynman diagrammatic series for Σ in terms of G is found to converge at least for moderately low temperatures. However, at strong interactions, its convergence is to the unphysical branch. This reveals a new scenario of breaking down of diagrammatic expansions. In contrast, the bare series in terms of the noninteracting Green's function G0 converges to the correct physical branch of Σ in all cases currently accessible by diagrammatic Monte Carlo calculations. In addition to their conceptual importance, these observations have important implications for techniques based on the explicit summation of the diagrammatic series.

Diagrammatic Monte Carlo method for many-polaron problems

We introduce the first bold diagrammatic Monte Carlo approach to deal with polaron problems at a finite electron density nonperturbatively, i.e., by including vertex corrections to high orders. Using the Holstein model on a square lattice as a prototypical example, we demonstrate that our method is capable of providing accurate results in the thermodynamic limit in all regimes from a renormalized Fermi liquid to a single polaron, across the nonadiabatic region where Fermi and Debye energies are of the same order of magnitude. By accounting for vertex corrections, the accuracy of the theoretical description is increased by orders of magnitude relative to the lowest-order self-consistent Born approximation employed in most studies. We also find that for the electron-phonon coupling typical for real materials, the quasiparticle effective mass increases and the quasiparticle residue decreases with increasing the electron density at constant electron-phonon coupling strength.

Flat histogram diagrammatic Monte Carlo method: calculation of the Green's function in imaginary time

The diagrammatic Monte Carlo (DiagMC) method is a numerical technique which samples the entire diagrammatic series of the Green's function in quantum many-body systems. In this work, we incorporate the flat histogram principle in the diagrammatic Monte Carlo method, and we term the improved version the "flat histogram diagrammatic Monte Carlo" method. We demonstrate the superiority of this method over the standard DiagMC in extracting the long-imaginary-time behavior of the Green's function, without incorporating any a priori knowledge about this function, by applying the technique to the polaron problem.