A combined variational and diagrammatic quantum Monte Carlo approach to the many-electron problem

Symmetry-Broken Perturbation Theory to Large Orders in Antiferromagnetic Phases

We introduce a spin-symmetry-broken extension of the connected determinant algorithm [Riccardo Rossi, Determinant diagrammatic Monte Carlo algorithm in the thermodynamic limit, Phys. Rev. Lett. 119, 045701 (2017).PRLTAO0031-900710.1103/PhysRevLett.119.045701]. The resulting systematic perturbative expansions around an antiferromagnetic state allow for numerically exact calculations directly inside a magnetically ordered phase. We show new precise results for the magnetic phase diagram and thermodynamics of the three-dimensional cubic Hubbard model at half-filling. With detailed computations of the order parameter in the low to intermediate-coupling regime, we establish the Néel phase boundary. The critical behavior in its vicinity is shown to be compatible with the O(3) Heisenberg universality class. By determining the evolution of the entropy with decreasing temperature through the phase transition we identify the different physical regimes at U/t=4. We provide quantitative results for several thermodynamic quantities deep inside the antiferromagnetic dome up to large interaction strengths and investigate the crossover between the Slater and Heisenberg regimes.

Static Self-Energy and Effective Mass of the Homogeneous Electron Gas from Quantum Monte Carlo Calculations

We discuss the methodology of quantum Monte Carlo calculations of the effective mass based on the static self-energy Σ(k,0). We then use variational Monte Carlo calculations of Σ(k,0) of the homogeneous electron gas at various densities to obtain results very close to perturbative G_{0}W_{0} calculations for values of the density parameter 1≤r_{s}≤10. The obtained values for the effective mass are close to diagrammatic Monte Carlo results and disagree with previous quantum Monte Carlo calculations based on a heuristic mapping of excitation energies to those of an ideal gas.

Fermionic physics from ab initio path integral Monte Carlo simulations of fictitious identical particles

The ab initio path integral Monte Carlo (PIMC) method is one of the most successful methods in statistical physics, quantum chemistry and related fields, but its application to quantum degenerate Fermi systems is severely hampered by an exponential computational bottleneck: the notorious fermion sign problem. Very recently, Xiong and Xiong [J. Chem. Phys. 157, 094112 (2022)] have suggested to partially circumvent the sign problem by carrying out simulations of fictitious systems guided by an interpolating continuous variable ξ ∈ [-1, 1], with the physical Fermi- and Bose-statistics corresponding to ξ = -1 and ξ = 1. It has been proposed that information about the fermionic limit might be obtained by calculations within the bosonic sector ξ > 0 combined with an extrapolation throughout the fermionic sector ξ < 0, essentially bypassing the sign problem. Here, we show how the inclusion of the artificial parameter ξ can be interpreted as an effective penalty on the formation of permutation cycles in the PIMC simulation. We demonstrate that the proposed extrapolation method breaks down for moderate to high quantum degeneracy. Instead, the method constitutes a valuable tool for the description of large Fermi-systems of weak quantum degeneracy. This is demonstrated for electrons in a 2D harmonic trap and for the uniform electron gas (UEG), where we find excellent agreement (∼0.5%) with exact configuration PIMC results in the high-density regime while attaining a speed-up exceeding 11 orders of magnitude. Finally, we extend the idea beyond the energy and analyze the radial density distribution (2D trap), as well as the static structure factor and imaginary-time density-density correlation function (UEG).

Electronic density response of warm dense hydrogen on the nanoscale

The properties of hydrogen at warm dense matter (WDM) conditions are of high importance for the understanding of astrophysical objects and technological applications such as inertial confinement fusion. In this work, we present extensive ab initio path integral Monte Carlo results for the electronic properties in the Coulomb potential of a fixed ionic configuration. This gives us unique insights into the complex interplay between the electronic localization around the protons with their density response to an external harmonic perturbation. We find qualitative agreement between our simulation data and a heuristic model based on the assumption of a local uniform electron gas model, but important trends are not captured by this simplification. In addition to being interesting in their own right, we are convinced that our results will be of high value for future projects, such as the rigorous benchmarking of approximate theories for the simulation of WDM, most notably density functional theory.

Assessing the accuracy of hybrid exchange-correlation functionals for the density response of warm dense electrons

We assess the accuracy of common hybrid exchange-correlation (XC) functionals (PBE0, PBE0-1/3, HSE06, HSE03, and B3LYP) within the Kohn-Sham density functional theory for the harmonically perturbed electron gas at parameters relevant for the challenging conditions of the warm dense matter. Generated by laser-induced compression and heating in the laboratory, the warm dense matter is a state of matter that also occurs in white dwarfs and planetary interiors. We consider both weak and strong degrees of density inhomogeneity induced by the external field at various wavenumbers. We perform an error analysis by comparing with the exact quantum Monte Carlo results. In the case of a weak perturbation, we report the static linear density response function and the static XC kernel at a metallic density for both the degenerate ground-state limit and for partial degeneracy at the electronic Fermi temperature. Overall, we observe an improvement in the density response when the PBE0, PBE0-1/3, HSE06, and HSE03 functionals are used, compared with the previously reported results for the PBE, PBEsol, local-density approximation, and AM05 functionals; B3LYP, on the other hand, does not perform well for the considered system. Additionally, the PBE0, PBE0-1/3, HSE06, and HSE03 functionals are more accurate for the density response properties than SCAN in the regime of partial degeneracy.

Dynamic Response of an Electron Gas: Towards the Exact Exchange-Correlation Kernel

Precise calculations of dynamics in the homogeneous electron gas (jellium model) are of fundamental importance for design and characterization of new materials. We introduce a diagrammatic Monte Carlo technique based on algorithmic Matsubara integration that allows us to compute frequency and momentum resolved finite temperature response directly in the real frequency domain using a series of connected Feynman diagrams. The data for charge response at moderate electron density are used to extract the frequency dependence of the exchange-correlation kernel at finite momenta and temperature. These results are as important for development of the time-dependent density functional theory for materials dynamics as ground state energies are for the density functional theory.

Evaluating Second-Order Phase Transitions with Diagrammatic Monte Carlo: Néel Transition in the Doped Three-Dimensional Hubbard Model

Diagrammatic Monte Carlo-the technique for the numerically exact summation of all Feynman diagrams to high orders-offers a unique unbiased probe of continuous phase transitions. Being formulated directly in the thermodynamic limit, the diagrammatic series is bound to diverge and is not resummable at the transition due to the nonanalyticity of physical observables. This enables the detection of the transition with controlled error bars from an analysis of the series coefficients alone, avoiding the challenge of evaluating physical observables near the transition. We demonstrate this technique by the example of the Néel transition in the 3D Hubbard model. At half filling and higher temperatures, the method matches the accuracy of state-of-the-art finite-size techniques, but surpasses it at low temperatures and allows us to map the phase diagram in the doped regime, where finite-size techniques struggle from the fermion sign problem. At low temperatures and sufficient doping, the transition to an incommensurate spin density wave state is observed.

Single-particle excitations in the uniform electron gas by diagrammatic Monte Carlo

We calculate the single-particle excitation spectrum and the Landau liquid parameters for the archetypal model of solids, the three-dimensional uniform electron gas, with the variational diagrammatic Monte Carlo method, which gives numerically controlled results without systematic error. In the metallic range of density, we establish benchmark values for the wave-function renormalization factor Z, the effective mass \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m^*/m$$\end{document}m∗/m, and the Landau parameters \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_0^s$$\end{document}F0s and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_0^a$$\end{document}F0a with unprecedented accuracy, and we resolve the long-standing puzzle of non-monotonic dependence of mass on density. We also exclude the possibility that experimentally measured large reduction of bandwidth in Na metal can originate from the charge and spin fluctuations contained in the model of the uniform electron gas.

Quasiparticle Effective Mass of the Three-Dimensional Fermi Liquid by Quantum Monte Carlo

According to Landau's Fermi liquid theory, the main properties of the quasiparticle excitations of an electron gas are embodied in the effective mass m^{*}, which determines the energy of a single quasiparticle, and the Landau interaction function, which indicates how the energy of a quasiparticle is modified by the presence of other quasiparticles. This simple paradigm underlies most of our current understanding of the physical and chemical behavior of metallic systems. The quasiparticle effective mass of the three-dimensional homogeneous electron gas has been the subject of theoretical controversy, and there is a lack of experimental data. In this Letter, we deploy diffusion Monte Carlo (DMC) methods to calculate m^{*} as a function of density for paramagnetic and ferromagnetic three-dimensional homogeneous electron gases. The DMC results indicate that m^{*} decreases when the density is reduced, especially in the ferromagnetic case. The DMC quasiparticle energy bands exclude the possibility of a reduction in the occupied bandwidth relative to that of the free-electron model at density parameter r_{s}=4, which corresponds to Na metal.

Real-Frequency Response Functions at Finite Temperature

Building on previous developments [A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlanc, Phys. Rev. B 99, 035120 (2019); PRBMDO2469-995010.1103/PhysRevB.99.035120A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlancPhys. Rev. B101, 125109 (2020); PRBMDO2469-995010.1103/PhysRevB.101.125109A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlancPhys. Rev. B102, 045115 (2020)PRBMDO2469-995010.1103/PhysRevB.102.045115, B. Holm and U. von Barth, Phys. Rev. B 57, 2108 (1998)PRBMDO0163-182910.1103/PhysRevB.57.2108, J. Vičičević and M. Ferrero, Phys. Rev. B 101, 075113 (2020)PRBMDO2469-995010.1103/PhysRevB.101.075113], we show that the diagrammatic Monte Carlo technique allows us to compute finite-temperature response functions directly on the real-frequency axis within any field-theoretical formulation of the interacting fermion problem. There are no limitations on the type and nature of the system's action or whether partial summation and self-consistent treatment of certain diagram classes are used. In particular, by eliminating the need for numerical analytic continuation from a Matsubara representation, our scheme allows us to study spectral densities of arbitrary complexity with controlled accuracy in models with frequency-dependent effective interactions. For illustrative purposes we consider the problem of the plasmon linewidth in a homogeneous electron gas (jellium).

Homotopic Action: A Pathway to Convergent Diagrammatic Theories

The major obstacle preventing Feynman diagrammatic expansions from accurately solving many-fermion systems in strongly correlated regimes is the series slow convergence or divergence problem. Several techniques have been proposed to address this issue: series resummation by conformal mapping, changing the nature of the starting point of the expansion by shifted action tools, and applying the homotopy analysis method to the Dyson-Schwinger equation. They emerge as dissimilar mathematical procedures aimed at different aspects of the problem. The proposed homotopic action offers a universal and systematic framework for unifying the existing-and generating new-methods and ideas to formulate a physical system in terms of a convergent diagrammatic series. It eliminates the need for resummation, allows one to introduce effective interactions, enables a controlled ultraviolet regularization of continuous-space theories, and reduces the intrinsic polynomial complexity of the diagrammatic Monte Carlo method. We illustrate this approach by an application to the Hubbard model.

Quantum Quasi-Monte Carlo Technique for Many-Body Perturbative Expansions

High order perturbation theory has seen an unexpected recent revival for controlled calculations of quantum many-body systems, even at strong coupling. We adapt integration methods using low-discrepancy sequences to this problem. They greatly outperform state-of-the-art diagrammatic Monte Carlo simulations. In practical applications, we show speed-ups of several orders of magnitude with scaling as fast as 1/N in sample number N; parametrically faster than 1/sqrt[N] in Monte Carlo simulations. We illustrate our technique with a solution of the Kondo ridge in quantum dots, where it allows large parameter sweeps.