Search

We present a minimal pole method for analytically continuing matrix-valued imaginary frequency correlation functions to the real axis, enabling precise access to off-diagonal elements and thus improving the interpretation of self-energies and susceptibilities in quantum simulations. Traditional methods for matrix-valued analytic continuation tend to be either noise-sensitive or make ad-hoc positivity assumptions. Our approach avoides these issues via the construction of a compact pole representation with shared poles through exponential fits, expanding upon prior work focused on scalar functions. We test our method across various scenarios, including fermionic and bosonic response functions, with and without noise, and for both continuous and discrete spectra of real materials and model systems. Our findings demonstrate that this technique addresses the shortcomings of existing methodologies, such as artificial broadening and positivity violations. The paper is supplemented with a sample implementation in Python.

We perform a fluctuation analysis of the pairing interaction in the hole-doped Hubbard model within the dynamical cluster approximation. Our analysis reveals that spin-fluctuation-mediated pairing differs qualitatively in the over- and underdoped regimes. In the underdoped regime spin fluctuations open a pseudogap. We show that in this regime the spin-fermion coupling mediates a giant attraction between antinodal fermions. As a consequence, despite the lack of coherent fermionic quasiparticles, a strong superconducting pairing is sustained down to small dopings., Comment: 6 pages, 2 figures

Quantum transport is often characterized not just by mean observables like the particle or energy current, but by their fluctuations and higher moments, which can act as detailed probes of the physical mechanisms at play. However, relatively few theoretical methods are able to access the full counting statistics (FCS) of transport processes through electronic junctions in strongly correlated regimes. While most experiments are concerned with the steady state properties, most accurate theoretical methods rely on computationally expensive propagation from a tractable initial state. Here, we propose a simple approach for computing the FCS through a junction directly at the steady state, utilizing the propagator noncrossing approximation (NCA). Compared to time propagation, our method offers reduced computational cost at the same level of approximation; but the idea can also be used within other approximations or as a basis for numerically exact techniques. We demonstrate the method's capabilities by investigating the impact of lead dimensionality on electronic transport in the nonequilibrium Anderson impurity model at the onset of Kondo physics. Our results reveal a distinct signature of one dimensional leads in the noise and Fano factor not present for other dimensionalities, showing the potential of FCS measurements as a probe of the environment surrounding a quantum dot.

By merging algorithmic Matsubara integration with discrete pole representations we present a procedure to generate fully analytic closed form results for impurity problems at fixed perturbation order. To demonstrate the utility of this approach we study the Bethe lattice and evaluate the second order self-energy for which reliable benchmarks exist. We show that, when evaluating diagrams on the Matsubara axis, the analytic sums of pole representations are extremely precise. We point out the absence of a numerical sign problem in the evaluation, and explore the application of the same procedure for real-frequency evaluation of diagrams. We find that real-frequency results are subject to noise that is controlled at low temperatures and can be mitigated at additional computational expense. We further demonstrate the utility of this approach by evaluating dynamical mean-field and bold diagrammatic self-consistency schemes at both second and fourth order and compare to benchmarks where available., Comment: 10 pages, 8 figures

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

The accurate ab initio simulation of molecules and periodic solids with diagrammatic perturbation theory is an important task in quantum chemistry, condensed matter physics, and materials science. In this article, we present the WeakCoupling module of the open-source software package Green, which implements fully self-consistent diagrammatic weak coupling simulations, capable of dealing with real materials in the finite-temperature formalism. The code is licensed under the permissive MIT license. We provide self-consistent GW (scGW) and self-consistent second-order Green's function perturbation theory (GF2) solvers, analysis tools, and post-processing methods. This paper summarizes the theoretical methods implemented and provides background, tutorials and practical instructions for running simulations., Comment: 23 pages, 2 figures

The exploitation of space group symmetries in numerical calculations of periodic crystalline solids accelerates calculations and provides physical insight. We present results for a space-group symmetry adaptation of electronic structure calculations within the finite-temperature self-consistent GW method along with an efficient parallelization scheme on accelerators. Our implementation employs the simultaneous diagonalization of the Dirac characters of the orbital representation. Results show that symmetry adaptation in self-consistent many-body codes results in substantial improvements of the runtime, and that block diagonalization on top of a restriction to the irreducible wedge results in additional speedup.

The non-equilibrium Green's function gives access to one-body observables for quantum systems. Of particular interest are quantities such as density, currents, and absorption spectra which are important for interpreting experimental results in quantum transport and spectroscopy. We present an integration scheme for the Green's function's equations of motion, the Kadanoff-Baym equations (KBE), which is both adaptive in the time integrator step size and method order as well as the history integration order. We analyze the importance of solving the KBE self-consistently and show that adapting the order of history integral evaluation is important for obtaining accurate results. To examine the efficiency of our method, we compare runtimes to a state of the art fixed time step integrator for several test systems and show an order of magnitude speedup at similar levels of accuracy.

Excitons, bound electron-hole pairs, influence the optical properties in strongly interacting solid state systems. Excitons and their associated many-body physics are typically most stable and pronounced in monolayer materials. Bulk systems with large exciton binding energies, on the other hand, are rare and the mechanisms driving their stability are still relatively unexplored. Here, we report an exceptionally large exciton binding energy in single crystals of the bulk van der Waals antiferromagnet CrSBr. Utilizing state-of-the-art angle-resolved photoemission spectroscopy and self-consistent ab-initio GW calculations, we present direct spectroscopic evidence that robust electronic and structural anisotropy can significantly amplify the exciton binding energy within bulk crystals. Furthermore, the application of a vertical electric field enables broad tunability of the optical and electronic properties. Our results indicate that CrSBr is a promising material for the study of the role of anisotropy in strongly interacting bulk systems and for the development of exciton-based optoelectronics.

Imaginary-time response functions of finite-temperature quantum systems are often obtained with methods that exhibit stochastic or systematic errors. Reducing these errors comes at a large computational cost -- in quantum Monte Carlo simulations, the reduction of noise by a factor of two incurs a simulation cost of a factor of four. In this paper, we relate certain imaginary-time response functions to an inner product on the space of linear operators on Fock space. We then show that data with noise typically does not respect the positive definiteness of its associated Gramian. The Gramian has the structure of a Hankel matrix. As a method for denoising noisy data, we introduce an alternating projection algorithm that finds the closest positive definite Hankel matrix consistent with noisy data. We test our methodology at the example of fermion Green's functions for continuous-time quantum Monte Carlo data and show remarkable improvements of the error, reducing noise by a factor of up to 20 in practical examples. We argue that Hankel projections should be used whenever finite-temperature imaginary-time data of response functions with errors is analyzed, be it in the context of quantum Monte Carlo, quantum computing, or in approximate semianalytic methodologies., Comment: 5 figures

We show that the amplitude mode in superconductors exhibits chirped oscillations under resonant excitation and that the chirping velocity increases as we approach the critical excitation strength. The chirped amplitude mode enables us to determine the local modification of the effective potential even when the system is in a long-lived pre-thermal state. We then show that this chirped amplitude mode is an experimentally observable quantity since the photo-induced (super)-current in pump-probe experiments serves as an efficient proxy for the dynamics of the order parameter, including the chirped dynamics. Our result is based on the attractive Hubbard model using dynamical mean-field theory within the symmetry-broken state after a resonant excitation across the superconducting gap. Since the collective response takes place on emergently long timescales, we extend the hierarchical low-rank compression method for nonequilibrium Green's functions to symmetry-broken states and show that it serves as an efficient representation despite long-lived memory kernels., Comment: 7 pages, 5 figures

Correlated electron systems may give rise to multiple effective interactions whose combined impact on quasiparticle properties can be difficult to disentangle. We introduce an unambiguous decomposition of the electronic self-energy which allows us to quantify the contributions of various effective interactions simultaneously. We use this tool to revisit the hole-doped Hubbard model within the dynamical cluster approximation, where commonly spin fluctuations are considered to be the origin of the pseudogap. While our fluctuation decomposition confirms that spin fluctuations indeed suppress antinodal electronic spectral weight, we show that they alone can not capture the pseudogap self-energy quantitatively. Nonlocal multi-boson Feynman diagrams yield substantial contributions and are needed for a quantitative description of the pseudogap., Comment: 5 pages, 4 figures

The many-body Green's function provides access to electronic properties beyond density functional theory level in ab inito calculations. In this manuscript, we propose a deep learning framework for predicting the finite-temperature Green's function in atomic orbital space, aiming to achieve a balance between accuracy and efficiency. By predicting the self-energy matrices in Lehmann representation using an equivariant message passing neural network, our method respects its analytical property and the $E(3)$ equivariance. The Green's function is obtained from the predicted self-energy through Dyson equation with target total number of electrons. We present proof-of-concept benchmark results for both molecules and simple periodic systems, showing that our method is able to provide accurate estimate of physical observables such as energy and density of states based on the predicted Green's function.

Analytical continuation is a central step in the simulation of finite-temperature field theories in which numerically obtained Matsubara data is continued to the real frequency axis for physical interpretation. Numerical analytic continuation is considered to be an ill-posed problem where uncertainties on the Matsubara axis are aplified exponentially. Here, we present a systematic and controlled procedure that approximates any Matsubara function by a minimal pole representation to within a predefined precision. We then show systematic convergence to the exact spectral function on the real axis as a function of our control parameter for a range of physically relevant setups. Our methodology is robust to noise and paves the way towards reliable analytic continuation in many-body theory and, by providing access to the analytic structure of the functions, direct theoretical interpretation of physical properties.

Computational models are an essential tool for the design, characterization, and discovery of novel materials. Hard computational tasks in materials science stretch the limits of existing high-performance supercomputing centers, consuming much of their simulation, analysis, and data resources. Quantum computing, on the other hand, is an emerging technology with the potential to accelerate many of the computational tasks needed for materials science. In order to do that, the quantum technology must interact with conventional high-performance computing in several ways: approximate results validation, identification of hard problems, and synergies in quantum-centric supercomputing. In this paper, we provide a perspective on how quantum-centric supercomputing can help address critical computational problems in materials science, the challenges to face in order to solve representative use cases, and new suggested directions., Comment: 65 pages, 15 figures; comments welcome

A description of long-lived photo-doped states in Mott insulators is challenging, as it needs to address exponentially separated timescales. We demonstrate how properties of such states can be computed using numerically exact steady state techniques, in particular Quantum Monte Carlo, by using a time-local ansatz for the distribution function with separate Fermi functions for the electron and hole quasiparticles. The simulations show that the Mott gap remains robust to large photo-doping, and the photo-doped state has hole and electron quasiparticles with strongly renormalized properties., Comment: 9 pages, 9 figures

We study bulk particle transport in a Fermi-Hubbard model on an infinite-dimensional Bethe lattice, driven by a constant electric field. Previous numerical studies showed that one dimensional analogs of this system exhibit a breakdown of diffusion due to Stark many-body localization (Stark-MBL) at least up to time which scales exponentially with the system size. Here, we consider systems initially in a spin density wave state using a combination of numerically exact and approximate techniques. We show that for sufficiently weak electric fields, the wave's momentum component decays exponentially with time in a way consistent with normal diffusion. By studying different wavelengths, we extract the dynamical exponent and the generalized diffusion coefficient at each field strength. Interestingly, we find a non-monotonic dependence of the dynamical exponent on the electric field. As the field increases towards a critical value proportional to the Hubbard interaction strength, transport slows down, becoming sub-diffusive. At large interaction strengths, however, transport speeds up again with increasing field, exhibiting super-diffusive characteristics when the electric field is comparable to the interaction strength. Eventually, at the large field limit, localization occurs and the current through the system is suppressed.

Response functions of quantum systems, such as electron Green's functions, magnetic, or charge susceptibilities, describe the response of a system to an external perturbation. They are the central objects of interest in field theories and quantum computing and measured directly in experiment. Response functions are intrinsically causal. In equilibrium and steady-state systems, they correspond to a positive spectral function in the frequency domain. Since response functions define an inner product on a Hilbert space and thereby induce a positive definite function, the properties of this function can be used to reduce noise in measured data and, in equilibrium and steady state, to construct positive definite extensions for data known on finite time intervals, which are then guaranteed to correspond to positive spectra.

We present the TRIQS/Nevanlinna analytic continuation package, an efficient implementation of the methods proposed by J. Fei et al in [Phys. Rev. Lett. 126, 056402 (2021)] and [Phys. Rev. B 104, 165111 (2021)]. TRIQS/Nevanlinna strives to provide a high quality open source (distributed under the GNU General Public License version 3) alternative to the more widely adopted Maximum Entropy based analytic continuation programs. With the additional Hardy functions optimization procedure, it allows for an accurate resolution of wide band and sharp features in the spectral function. Those problems can be formulated in terms of imaginary time or Matsubara frequency response functions. The application is based on the TRIQS C++/Python framework, which allows for easy interoperability with other TRIQS-based applications, electronic band structure codes and visualization tools. Similar to other TRIQS packages, it comes with a convenient Python interface., Comment: 15 pages, 5 figures

We present a temperature extrapolation technique for self-consistent many-body methods, which provides a causal starting point for converging to a solution at a target temperature. The technique employs the Carath\'eodory formalism for interpolating causal matrix-valued functions and is applicable to various many-body methods, including dynamical mean field theory, its cluster extensions, and self-consistent perturbative methods such as the self-consistent GW approximation. We show results that demonstrate that this technique can efficiently simulate heating and cooling hysteresis at a first-order phase transition, as well as accelerate convergence., Comment: 11 pages, 9 figures

We introduce a Julia implementation of the recently proposed Nevanlinna analytic continuation method. The method is based on Nevanlinna interpolants and, by construction, preserves the causality of a response function. For theoretical calculations without statistical noise, this continuation method is a powerful tool to extract real-frequency information from numerical input data on the Matsubara axis. This method has been applied to first-principles calculations of correlated materials. This paper presents its efficient and full-featured open-source implementation of the method including the Hamburger moment problem and smoothing., Comment: Submission to SciPost

The Hubbard model on an anisotropic triangular lattice in two dimensions, a fundamental model for frustrated electron physics, displays a wide variety of phases and phase transitions. This work investigates the model using the ladder dual fermion approximation which captures local correlations non-perturbatively but approximates non-local correlations. We find metallic, one-dimensional antiferromagnetic, non-collinear antiferromagnetic, square-lattice antiferromagnetic, and spiral phases but no evidence of collinear antiferromagnetic order in different parts of the phase diagram. Analyzing the spin susceptibility in detail, we see both regions of agreement and of discrepancy with previous work. The case of Cs$_2$CuCl$_4$ is discussed in detail., Comment: 13 pages, 6 figures

Green's functions of fermions are described by matrix-valued Herglotz-Nevanlinna functions. Since analytic continuation is fundamentally an ill-posed problem, the causal space described by the matrix-valued Herglotz-Nevanlinna structure can be instrumental in improving the accuracy and in enhancing the robustness with respect to noise. We demonstrate a three-pronged procedure for robust analytic continuation called PES: (1) Projection of data to the causal space. (2) Estimation of pole locations. (3) Semidefinite relaxation within the causal space. We compare the performance of PES with the recently developed Nevanlinna and Carath\'{e}odory continuation methods and find that PES is more robust in the presence of noise and does not require the usage of extended precision arithmetics. We also demonstrate that a causal projection improves the performance of the Nevanlinna and Carath\'{e}odory methods. The PES method is generalized to bosonic response functions, for which the Nevanlinna and Carath\'{e}odory continuation methods have not yet been developed. It is particularly useful for studying spectra with sharp features, as they occur in the study of molecules and band structures in solids., Comment: 14 pages, 5 figures

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with non-local correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles., Comment: 12 pages, 10 figures

We present a numerically exact steady-state inchworm Monte Carlo method for nonequilibrium quantum impurity models. Rather than propagating an initial state to long times, the method is directly formulated in the steady-state. This eliminates any need to traverse the transient dynamics and grants access to a much larger range of parameter regimes at vastly reduced computational costs. We benchmark the method on equilibrium Green's functions of quantum dots in the noninteracting limit and in the unitary limit of the Kondo regime. We then consider correlated materials described with dynamical mean field theory and driven away from equilibrium by a bias voltage. We show that the response of a correlated material to a bias voltage differs qualitatively from the splitting of the Kondo resonance observed in bias-driven quantum dots.

We present algorithmic and implementation details for the fully self-consistent finite-temperature $GW$ method in Gaussian Bloch orbitals for solids. Our implementation is based on the finite-temperature Green's function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of the self-energy are explicitly evaluated. The method is tested by evaluating the band gaps of selected semiconductors and insulators. We show agreement with other, differently formulated finite-temperature sc$GW$ implementations when finite-size corrections and basis set errors are taken into account. By migrating computationally intensive calculations to GPUs, we obtain scalable results on large supercomputers with nearly optimal performance. Our work demonstrates the applicability of Gaussian orbital based sc$GW$ for $\emph{ab initio}$ correlated materials simulations and provides a sound starting point for embedding methods built on top of $GW$., Comment: 17 pages, 10 figures, 2 tables

The real-time contour formalism for Green's functions provides time-dependent information of quantum many-body systems. In practice, the long-time simulation of systems with a wide range of energy scales is challenging due to both the storage requirements of the discretized Green's function and the computational cost of solving the Dyson equation. In this manuscript, we apply a real-time discretization based on a piece-wise high-order orthogonal-polynomial expansion to address these issues. We present a superconvergent algorithm for solving the real-time equilibrium Dyson equation using the Legendre spectral method and the recursive algorithm for Legendre convolution. We show that the compact high order discretization in combination with our Dyson solver enables long-time simulations using far fewer discretization points than needed in conventional multistep methods. As a proof of concept, we compute the molecular spectral functions of H$_2$, LiH, He$_2$ and C$_6$H$_4$O$_2$ using self-consistent second-order perturbation theory and compare the results with standard quantum chemistry methods as well as the auxiliary second-order Green's function perturbation theory method.

We study the fluctuations responsible for pairing in the $d$-wave superconducting state of the two-dimensional Hubbard model at intermediate coupling within a cluster dynamical mean-field theory with a numerically exact quantum impurity solver. By analyzing how momentum and frequency dependent fluctuations generate the $d-$wave superconducting state in different representations, we identify antiferromagnetic fluctuations as the pairing glue of superconductivity both in the underdoped and the overdoped regime. Nevertheless, in the intermediate coupling regime, the predominant magnetic fluctuations may differ significantly from those described by conventional spin-fluctuation theory.

We study the contribution of the electron-spin fluctuation coupling to the superconducting state of the two dimensional Hubbard model within dynamical cluster approximation (DCA) using a numerical exact continuous time Monte Carlo solver. By analyzing the frequency dependence of the self energy, we show that only about half of the superconductivity can be attributed to a "pairing glue" arising from treating spin fluctuations as a pairing boson in the standard one-loop theory.

We present a formulation of relativistic self-consistent $GW$ for solids based on the exact two-component formalism with one-electron approximation (X2C1e) and non-relativistic Coulomb interactions. Our theory allows us to study scalar relativistic effects, spin-orbit coupling, and the interplay of relativistic effects with electron correlation without adjustable parameters. Our all-electron implementation is fully $ab$ $initio$ and does not require a pseudopotential constructed from atomic calculations. We examine the effect of the X2C1e approximation by comparison to the established four-component formalism and reach excellent agreement. The simplicity of X2C1e enables the construction of higher order theories, such as embedding theories, on top of perturbative calculations., Comment: 18 pages, 5 figures, 7 tables

Multi-orbital quantum impurity models with general interaction and hybridization terms appear in a wide range of applications including embedding, quantum transport, and nanoscience. However, most quantum impurity solvers are restricted to a few impurity orbitals, discretized baths, diagonal hybridizations, or density-density interactions. Here, we generalize the inchworm quantum Monte Carlo method to the interaction expansion and explore its application to typical single- and multi-orbital problems encountered in investigations of impurity and lattice models. Our implementation generically outperforms bare and bold-line quantum Monte Carlo algorithms in the interaction expansion. So far, for the systems studied here, it remains inferior to the more specialized hybridization expansion and auxiliary field algorithms. The problem of convergence to unphysical fixed points, which hampers so-called bold-line methods, is not encountered in inchworm Monte Carlo.

We present a theory of modified reduced dynamics in the presence of counting fields. Reduced dynamics techniques are useful for describing open quantum systems at long emergent timescales when the memory timescales are short. However, they can be difficult to formulate for observables spanning the system and its environment, such as those characterizing transport properties. A large variety of mixed system--environment observables, as well as their statistical properties, can be evaluated by considering counting fields. Given a numerical method able to simulate the field-modified dynamics over the memory timescale, we show that the long-lived full counting statistics can be efficiently obtained from the reduced dynamics. We demonstrate the utility of the technique by computing the long-time current in the nonequilibrium Anderson impurity model from short-time Monte Carlo simulations.

The accurate determination of magnetic phase transitions in electronic systems is an important task of solid state theory. While numerically exact results are readily available for model systems such as the half-filled 3D Hubbard model, the complexity of real materials requires additional approximations, such as the restriction to certain classes of diagrams in perturbation theory, that reduce the precision with which magnetic properties are described. In this work, we examine the description of magnetic properties in second order perturbation theory, GW, FLEX, and two TMatrix approximations to numerically exact CT-QMC reference data. We assess finite-size effects and compare periodic lattice simulations to cluster embedding. We find that embedding substantially improves finite size convergence. However, by analyzing different partial summation methods we find no systematic improvement in the description of magnetic properties, with most methods considered in this work predicting first-order instead of continuous transitions, leading us to the conclusion that non-perturbative methods are necessary for the accurate determination of magnetic properties and phase transitions., Comment: 14 pages, 12 figures