Search NASA⌕ Search

SEARCH · Search NASA

Results for “diagrammatic Monte Carlo”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

18 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Diagrammatic Monte Carlo Approach to Real‐frequency Response Functions and the Spin‐Fermion Model of Hot Spots (Final Report)

This three‐year research project was focused on two outstanding condensed matter problems (i) Real‐frequency response functions in Coulomb systems at finite temperature, and Fermi surface (FS) reconstruction in the Spin-Fermion (SF) model. The uniform electron gas (jellium) model, describing electrons interacting via long‐range Coulomb forces on positive neutralizing background, is fundamentally important both for understanding the physics of correlated electrons and for formulation of the time‐dependent density functional theory (TDDFT). The spin-fermion model has found a wide range of applications in the physics of cuprates and iron‐based superconductors to explain the “strange metal" behavior and to suggest a possible pairing mechanism for high‐temperature superconductivity. The project’s goals were development and application of the Diagrammatic Monte Carlo (DiagMC) techniques for solving the above problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

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); A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlancPhys. Rev. B101, 125109 (2020); A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlancPhys. Rev. B102, 045115 (2020), B. Holm and U. von Barth, Phys. Rev. B 57, 2108 (1998), J. Vičičevic and M. Ferrero, Phys. Rev. B 101, 075113 (2020)], 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. Finally, for illustrative purposes we consider the problem of the plasmon linewidth in a homogeneous electron gas (jellium).

3-dimensional systems↗
DOE OSTI2022

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. In this work, 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 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.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

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. Here, 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.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Interaction-enhanced nesting in spin-fermion and Fermi-Hubbard models

The spin-fermion (SF) model postulates that the dominant coupling between low-energy fermions in near critical metals is mediated by collective spin fluctuations (paramagnons) peaked at the Néel wave vector, Q N , connecting hot spots on opposite sides of the Fermi surface. It has been argued that strong correlations at hot spots lead to a Fermi surface deformation (FSD) featuring flat regions and increased nesting. This conjecture was confirmed in the perturbative self-consistent calculations when the paramagnon propagator dependence on momentum deviation from Q N is given by χ − 1 ∝ | Δ q | . Using diagrammatic Monte Carlo (diagMC) technique we show that such a dependence holds only at temperatures orders of magnitude smaller than any other energy scale in the problem, indicating that a different mechanism may be at play. Instead, we find that a χ − 1 ∝ | Δ q | 2 dependence yields a robust finite- T scenario for achieving FSD. To link phenomenological and microscopic descriptions, we applied the connected determinant diagMC method to the ( t − t ′ ) Hubbard model and found that at large U / t > 5.5 before the formation of electron and hole pockets (i) the FSD defined as a maximum of the spectral function is not very pronounced; instead, it is the lines of zeros of the renormalized dispersion relation that deforms toward nesting, and (ii) the static spin susceptibility is well described by χ − 1 ∝ | Δ q | 2 . Flat FS regions yield a nontrivial scenario for realizing a non-Fermi liquid state. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Many-body Green's function approaches to the doped Fröhlich solid: Exact solutions and anomalous mass enhancement

In polar semiconductors and insulators, the Fröhlich interaction between electrons and long-wavelength longitudinal optical phonons induces a many-body renormalization of the carrier effective masses and the appearance of characteristic phonon sidebands in the spectral function, commonly dubbed “polaron satellites.” The simplest model that captures these effects is the Fröhlich model, whereby electrons in a parabolic band interact with a dispersionless longitudinal optical phonon. The Fröhlich model has been employed in a number of seminal papers, from early perturbation-theory approaches to modern diagrammatic Monte Carlo calculations. One limitation of this model is that it focuses on undoped systems, thus ignoring carrier screening and Pauli blocking effects that are present in real experiments on doped samples. Here, to overcome this limitation, we extend the Fröhlich model to the case of doped systems, and we provide exact solutions for the electron spectral function, mass enhancement, and polaron satellites. We perform the analysis using two approaches, namely, Dyson’s equation with the Fan-Migdal self-energy, and the second-order cumulant expansion. We find that these two approaches provide qualitatively different results. In particular, Dyson’s approach yields better quasiparticle masses and worse satellites, while the cumulant approach provides better satellite structures, at the price of worse quasiparticle masses. Both approaches yield an anomalous enhancement of the electron effective mass at finite doping levels, which in turn leads to a breakdown of the quasiparticle picture in a significant portion of the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Ab initio self-consistent many-body theory of polarons at all couplings

We present a theoretical framework to describe polarons from first principles within a many-body Green’s function formalism. Starting from a general electron-phonon Hamiltonian, we derive a self-consistent Dyson equation in which the phonon-mediated self-energy is composed by two distinct terms. One term is the Fan-Migdal self-energy and describes dynamic electron-phonon processes, the other term is a contribution to the self-energy originating from the static displacements of the atomic nuclei in the polaronic ground state. The lowest-order approximation to the present theory yields the standard many-body perturbation theory approach to electron-phonon interactions in the limit of large polarons, and the ab initio polaron equations introduced in the limit of small polarons. Here, a practical recipe to implement the present unifying formalism in first-principles calculations is outlined. We apply our method to the Fröhlich model, and obtain remarkably accurate polaron energies at all couplings, in line with Feynman’s polaron theory and diagrammatic Monte Carlo calculations. We also recover the well-known results of Fröhlich and Pekar at weak and strong coupling, respectively. The present approach enables predictive many-body calculations of polarons in real materials at all couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Unified Approach to Polarons and Phonon-Induced Band Structure Renormalization

Ab initio calculations of the phonon-induced band structure renormalization are currently based on the perturbative Allen-Heine theory and its many-body generalizations. These approaches are unsuitable to describe materials where electrons form localized polarons. Here, we develop a self-consistent, many-body Green’s function theory of band structure renormalization that incorporates localization and self-trapping. Here, we show that the present approach reduces to the Allen-Heine theory in the weak-coupling limit, and to total energy calculations of self-trapped polarons in the strong-coupling limit. To demonstrate this methodology, we reproduce the path-integral results of Feynman and diagrammatic Monte Carlo calculations for the Fröhlich model at all couplings, and we calculate the zero point renormalization of the band gap of an ionic insulator including polaronic effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Understanding polaronic transport in complex oxides by combining precise synthesis and first-principles many-body theory

In complex oxides, charge carriers often couple strongly with lattice vibrations to form polarons–entangled electron–phonon quasiparticles whose transport properties remain difficult to characterize. Experimental access to intrinsic polaronic transport requires ultraclean samples, while theoretical description demands methods beyond low-order perturbation theory. Here, we show a predictive theory–experiment workflow to study polaron transport in complex oxides. Focusing on a prototypical polaronic oxide, anatase TiO 2 , we combine growth of high-quality oxygen-vacancy-doped films using hybrid molecular beam epitaxy with a first-principles electron–phonon diagrammatic Monte-Carlo (FEP-DMC) framework recently developed for accurate polaron predictions. Our films exhibit record-high electron mobility for anatase TiO 2 , in excellent agreement with FEP-DMC calculations conducted prior to experiment, which predict a room-temperature mobility of 45 ± 15 cm −2 V −1 s −1 and a mobility-temperature scaling of μ ∝ T −1.9 ± 0.077 . Microscopic analysis using scanning transmission electron microscopy and x-ray photoelectron spectroscopy reveals the role of oxygen vacancies in modulating transport at lower temperatures. FEP-DMC further provides quantitative insight into polaron formation energy, phonon cloud distribution, lattice distortion around the polaron, and the polaronic contribution to mobility. Together, these results provide a deeper microscopic understanding of large-polaron transport in a complex oxide and provide the blueprint to characterize other polaronic materials.

anatase TiO2↗
DOE OSTI2023

Tensor train continuous time solver for quantum impurity models

The simulation of strongly correlated quantum impurity models is a significant challenge in modern condensed matter physics that has multiple important applications. Thus far, the most successful methods for approaching this challenge involve Monte Carlo techniques that accurately and reliably sample perturbative expansions to any order. However, the cost of obtaining high precision through these methods is high. Recently, tensor train decomposition techniques have been developed as an alternative to Monte Carlo integration. In this study, we apply these techniques to the single-impurity Anderson model at equilibrium by calculating the systematic expansion in power of the hybridization of the impurity with the bath. Furthermore, we demonstrate the performance of the method in a paradigmatic application, examining the first-order phase transition on the infinite-dimensional Bethe lattice, which can be mapped to an impurity model through dynamical mean field theory. Our results indicate that using tensor train decomposition schemes allows the calculation of finite-temperature Green's functions and thermodynamic observables with unprecedented accuracy. The methodology holds promise for future applications to frustrated multiorbital systems, using a combination of partially summed series with other techniques pioneered in diagrammatic and continuous time quantum Monte Carlo.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Krylov spaces for truncated spectrum methodologies

We propose herein an extension of truncated spectrum methodologies, a nonperturbative numerical approach able to elucidate the low energy properties of quantum field theories. TSMs, in their various flavors, involve a division of a computational Hilbert space, H , into two parts, one part, H 1 that is “kept” for the numerical computations, and one part, H 2 , that is discarded or “truncated.” Even though H 2 is discarded, truncated spectrum methodologies will often try to incorporate the effects of H 2 in some effective way. In these terms, we propose to keep the dimension of H 1 small. We pair this choice of H 1 with a Krylov subspace iterative approach able to take into account the effects of H 2 . This iterative approach can be taken to arbitrarily high order and so offers the ability to compute quantities to arbitrary precision. In many cases it also offers the advantage of not needing an explicit UV cutoff. To compute the matrix elements that arise in the Krylov iterations, we employ a Feynman diagrammatic representation that is then evaluated with Monte Carlo techniques. Each order of the Krylov iteration is variational and is guaranteed to improve upon the previous iteration. The first Krylov iteration is akin to the next-to-leading order approach of Elias-Miró [NLO renormalization in the Hamiltonian truncation, ]. To demonstrate this approach, we focus on the ( 1 + 1 d )-dimensional ϕ 4 model and compute the bulk energy and mass gaps in both the Z 2 -broken and unbroken sectors. We estimate the critical ϕ 4 coupling in the broken phase to be g c = 0.2645 ± 0.002 . Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Interaction-expansion inchworm Monte Carlo solver for lattice and impurity models

Multiorbital 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. We generalize the inchworm quantum Monte Carlo method to the interaction expansion, and we explore its application to typical single- and multiorbital 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. For the systems studied here, our implementation remains inferior to the more specialized hybridization expansion and auxiliary field algorithms. So, the problem of convergence to unphysical fixed points, which hampers so-called bold-line methods, is not encountered in inchworm Monte Carlo.

36 MATERIALS SCIENCE↗
DOE OSTI2021

How correlations change the magnetic structure factor of the kagome Hubbard model

The kagome Hubbard model (KHM) is a paradigmatic example of a frustrated two-dimensional model. While its strongly correlated regime, described by a Heisenberg model, is of topical interest due to its enigmatic prospective spin-liquid ground state, the weakly and moderately correlated regimes remain largely unexplored. Motivated by the rapidly growing number of metallic kagome materials (e.g., Mn 3 Sn, Fe 3 Sn 2 , FeSn, Co 3 Sn 2 S 2 , Gd 3 Ru 4 Al 12 , and AV 3 Sb 5 with A = K, Rb, Cs), here we study the respective regimes of the KHM by means of three complementary numerical methods: the dynamical mean-field theory, the dynamical vertex approximation, and determinant quantum Monte Carlo. In contrast to the archetypal square lattice, we find no tendencies toward magnetic ordering, as magnetic correlations remain short-range. Nevertheless, the magnetic correlations undergo a remarkable crossover as the system approaches the metal-to-insulator transition. The Mott transition itself does not affect the magnetic correlations. Our equal-time and dynamical structure factors can be used as a reference for inelastic neutron scattering experiments on the growing family of metallic kagome materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Accelerated impurity solver for DMFT and its diagrammatic extensions

Here, we present ComCTQMC, a GPU accelerated quantum impurity solver. It uses the continuous-time quantum Monte Carlo (CTQMC) algorithm wherein the partition function is expanded in terms of the hybridisation function (CT-HYB). ComCTQMC supports both partition and worm-space measurements, and it uses improved estimators and the reduced density matrix to improve observable measurements whenever possible. ComCTQMC efficiently measures all one and two-particle Green's functions, all static observables which commute with the local Hamiltonian, and the occupation of each impurity orbital. ComCTQMC can solve complex-valued impurities with crystal fields that are hybridized to both fermionic and bosonic baths. Most importantly, ComCTQMC utilizes graphical processing units (GPUs), if available, to dramatically accelerate the CTQMC algorithm when the Hilbert space is sufficiently large. We demonstrate acceleration by a factor of over 600 (100) in a simulation of δ-Pu at 600 K with (without) crystal fields. In easier problems, the GPU offers less impressive acceleration or even decelerates the CTQMC. Here we describe the theory, algorithms, and structure used by ComCTQMC in order to achieve this set of features and level of acceleration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Charge correlations suppress unconventional pairing in the Holstein model

In a recent work by Schrodi et al. [Phys. Rev. B. 104, L140506 (2021)], the authors find an unconventional superconducting state with a sign-changing order parameter using the Migdal-Eliashberg theory, including the first vertex correction. This unconventional solution arises despite using an isotropic bare electron-phonon coupling in the Hamiltonian. We examine this claim using hybrid quantum Monte Carlo for a single-band Holstein model with a cuprate-like noninteracting band structure and identical parameters to Schrodi et al.. Our Monte Carlo results for these parameters suggest that unconventional pairing correlations do not exceed their noninteracting values at any carrier concentration we have checked. Instead, strong charge-density-wave correlations persist at the lowest accessible temperatures for dilute and nearly half-filled bands. As a result, we present arguments for how vertex-corrected Migdal-Eliashberg calculation schemes can lead to uncontrolled results in the presence of Fermi surface nesting.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quantum Monte Carlo Method in the Steady State

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.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Efficient phase-space generation for hadron collider event simulation

We present a simple yet efficient algorithm for phase-space integration at hadron colliders. Individual mappings consist of a single t-channel combined with any number of s-channel decays, and are constructed using diagrammatic information. The factorial growth in the number of channels is tamed by providing an option to limit the number of s-channel topologies. We provide a publicly available, parallelized code in C++ and test its performance in typical LHC scenarios.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Advancing the understanding of energy-energy correlators in heavy-ion collisions

We investigate the collinear limit of the energy-energy correlator (EEC) in a heavy-ion context. First, we revisit the leading-logarithmic (LL) resummation of this observable in vacuum following a diagrammatic approach. We argue that this route allows to naturally incorporate medium-induced effects into the all-orders structure systematically. As an example, we show how the phase-space constraints imposed by the static medium on vacuum-like emissions can be incorporated into the LL result by modifying the anomalous dimensions. On the fixed-order side, we calculate the $\mathcal{O}$(α s ) expansion of the in-medium EEC for a γ → $q\overline{q}$ splitting with arbitrary kinematics including, for the first time, subleading colour corrections. When comparing this result to previously used approximations in the literature, we find up to $\mathcal{O}$(1) deviations in the regime of interest for jet quenching signatures. Energy loss effects are also quantified and further suppress the EEC at large angles. These semi-analytic studies are complemented with a phenomenological study using the jet quenching Monte Carlo JetMed. Finally, we argue that the imprint of medium-induced effects in energy-energy correlators can be enhanced by using an alternative definition that takes as input Lund primary declusterings instead of particles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗