Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dyson equation”

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

BPZ equations for higher degenerate fields and nonperturbative Dyson-Schwinger equations

In the two-dimensional Liouville conformal field theory, correlation functions involving a degenerate field satisfy partial differential equations due to the decoupling of the null descendant field. On the other hand, the instanton partition function of a four-dimensional N = 2 supersymmetric theory in the Ω -background at a special point of the parameter space also satisfies a partial differential equation resulting from the constraints of the gauge field configurations. This partial differential equation can be proved using the nonperturbative Dyson-Schwinger equations. We show for the next-to-simplest case that the partial differential equations obtained from two different perspectives can be identified, thereby confirming an assertion of the Bogomol’nyi-Prasad-Sommerfield/conformal field theory correspondence. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

Semiclassical Trans-Series from the Perturbative Hopf-Algebraic Dyson-Schwinger Equations: $\phi$ 3 QFT in 6 Dimensions

We analyze the asymptotically free massless scalar $\phi$ 3 quantum field theory in 6 dimensions, using resurgent asymptotic analysis to find the trans-series solutions which yield the non-perturbative completion of the divergent perturbative solutions to the Kreimer–Connes Hopf-algebraic Dyson–Schwinger equations for the anomalous dimension. This scalar conformal field theory is asymptotically free and has a real Lipatov instanton. In the Hopf-algebraic approach we find a trans-series having an intricate Borel singularity structure, with three distinct but resonant non-perturbative terms, each repeated in an infinite series. These expansions are in terms of the renormalized coupling. The resonant structure leads to powers of logarithmic terms at higher levels of the trans-series, analogous to logarithmic terms arising from interactions between instantons and anti-instantons, but arising from a purely perturbative formalism rather than from a semi-classical analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Electromagnetic wave scattering in a two-layer anisotropic random medium

For electromagnetic wave propagation and scattering in an anisotropic random medium, the Dyson equation for the mean field and the Bethe-Salpeter equation for the correlation or the covariance of the field were derived. With the random permittivity expressed in a general anisotropic form, the bilocal and the nonlinear approximations are employed to solve the Dyson equation, and the ladder approximation to solve the Bethe-Salpeter equation. The mean dyadic Green's function for a two-layer anisotropic random medium with arbitrary three-dimensional correlation functions has been investigated with the zeroth-order solutions to the Dyson equation under the nonlinear approximation. The effective propagation constants are calculated for the four characteristic waves associated with the coherent vector fields, propagating in an anisotropic random-medium layer, which are the ordinary and extraordinary waves with upward- and downward-propagating vectors.

Lee, J. K.↗
NASA NTRS1985

Remote sensing of earth terrain

Progress on the investigation of the anisotropy of the terrain media, such as vegetation canopy and sea ice, and the study of the fluctuation-dissipation theorem in conjunction with the application of strong fluctuation theory for passive remote sensing of snowpacks is reported. The Feynman diagrammatic technique is used to derive the Dyson equation for the mean field and the Bethe-Salpeter equation for the correlation or the covariance of the field for electromagnetic wave propagation and scattering in an anisotropic random medium. With the random permittivity expressed in a general form, the bilocal and the nonlinear approximations are employed to solve the Dyson equation and the ladder approximation to the Bethe-Salpeter equation. The mean dyadic Green's function for a two layer anisotropic random medium with arbitrary three dimensional correlation function was investigated with the zeroth-order solutions to the Dyson equation under the four characteristic waves associated with the coherent vector fields propagating in an anisotropic random medium layer, which are the ordinary and extraordinary waves with upward and downward propagating vectors.

Kong, J. A.↗
DOE OSTI2023

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Split representation of adaptively compressed polarizability operator

The polarizability operator plays a central role in density functional perturbation theory and other perturbative treatment of first principle electronic structure theories. The cost of computing the polarizability operator generally scales as O(Ne4) where Ne is the number of electrons in the system. The recently developed adaptively compressed polarizability operator (ACP) formulation [L. Lin, Z. Xu and L. Ying, Multiscale Model. Simul. 2017] reduces such complexity to O(Ne 3 ) in the context of phonon calculations with a large basis set for the first time, and demonstrates its effectiveness for model problems. In this paper, we improve the performance of the ACP formulation by splitting the polarizability into a near singular component that is statically compressed, and a smooth component that is adaptively compressed. The new split representation maintains the O(Ne 3 ) complexity, and accelerates nearly all components of the ACP formulation, including Chebyshev interpolation of energy levels, iterative solution of Sternheimer equations, and convergence of the Dyson equations. For simulation of real materials, we discuss how to incorporate nonlocal pseudopotentials and finite temperature effects. In this work, we demonstrate the effectiveness of our method using one-dimensional model problem in insulating and metallic regimes, as well as its accuracy for real molecules and solids.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2018

An Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. II: Approach Based on the Dyson and Bethe-Salpeter Equations

In this paper, the vector radiative transfer equation is derived by means of the vector integral Foldy equations describing the electromagnetic scattering by a group of particles. By assuming that in a discrete random medium the positions of the particles are statistically independent and by applying the Twersky approximation to the order-of-scattering expansion of the total field, we derive the Dyson equation for the coherent field and the ladder approximated Bethe–Salpeter equation for the dyadic correlation function. Then, under the far-field assumption for sparsely distributed particles, the Dyson equation is reduced to the Foldy integral equation for the coherent field, while the iterated solution of the Bethe–Salpeter equation ultimately yields the vector radiative transfer equation.

Electromagnetic scattering↗
NASA NTRS1980

Modified radiative transfer theory for a two-layer random medium

Modified radiative transfer (MRT) equations appropriate for electromagnetic wave propagation in bounded random media are derived from the Bethe-Salpeter equation with the ladder approximation and the Dyson equation with the nonlinear approximation. The MRT equations are then solved with the first-order renormalization approximation to obtain analytical results for the backscattering cross sections of a two-layer random medium with arbitrary three-dimensional correlation functions. The coherent effects of the MRT theory are illustrated and comparisons are made with backscattering cross sections obtained with the first Born approximation to the wave equation.

Zuniga, M. A.↗
NASA NTRS1988

Modified radiative transfer theory for a two-layer anisotropic random medium

The modified radiative transfer (MRT) equations which describe propagation and scattering of the electromagnetic field intensity in a layered anisotropic random medium are derived from the Bethe-Salpeter equation with the ladder approximation and the Dyson equation with the nonlinear approximation. The extinction matrices and scattering phase matrices for both mean and incoherent field intensities are obtained. Backscattering enhancement is observed due to the spatial coherent effect between upward and downward propagating waves. It also occurs for the half-space case, because of coupling between ordinary and extraordinary waves. The depolarization effect is predicted in the first-order renormalization approximation to the MRT equation.

Lee, J. K.↗
DOE OSTI2022

Replica symmetry breaking for the integrable two-site Sachdev–Ye–Kitaev model

We analyze a two-body non-Hermitian two-site Sachdev–Ye–Kitaev (SYK) model with the couplings of one site complex conjugated to the other site. This model, with no explicit coupling between the sites, shows an infinite number of second-order phase transitions, which is a consequence of the factorization of the partition function into a product over Matsubara frequencies. We calculate the quenched free energy in two different ways: first in terms of the single-particle energies and second by solving the Schwinger–Dyson equations of the two-site model. The first calculation can be done entirely in terms of a one-site model. The conjugate replica enters due to non-analyticities when Matsubara frequencies enter the spectral support of the coupling matrix. The second calculation is based on the replica trick of the two-site partition function. Both methods give the same result. The free-fermion partition function can also be rephrased as a matrix model for the coupling matrix. Up to minor details, this model is the random matrix model that describes the chiral phase transition of QCD, and the order parameter of the two-body model corresponds to the chiral condensate of QCD. Comparing to the corresponding four-body model, we are able to determine which features of the free energy are due to the chaotic nature of the four-body model. The high-temperature phase of both models is entropy dominated, and in both cases, the free energy is determined by the spectral density. The chaotic four-body SYK model has a low-temperature phase whose free energy is almost temperature-independent, signaling an effective gap of the theory even though the actual spectrum does not exhibit a gap. On the other hand, the low-temperature free energy of the two-body SYK model is not flat; in fact, it oscillates to arbitrarily low temperature. This indicates a less desirable feature that the entropy of the two-body model is not always positive in the low-temperature phase, which most likely is a consequence of the non-hermiticity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Equation-of-Motion Coupled-Cluster Cumulant Green’s Function for Excited States and X-ray Spectra

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a set of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function associated with the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then obtained as a convolution of the core spectral function and an effective one-body absorption spectrum. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. Themethod is illustrated for the XPS and XAS of NH 3 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Field redefinitions can be nonlocal

We revisit the lore establishing the allowed space of field redefinitions and show that there are essentially no restrictions. Our conclusions hold to all orders in perturbation theory and for any dispersion relation. Field redefinitions can be nonlocal, symmetry breaking, or in certain cases have explicit dependence on spacetime. We address field redefinitions that can be resummed into the propagator, which demonstrates how to perform perturbative calculations away from the minimum in field space. Field redefinitions are used to derive higher-order Schwinger-Dyson equations, which imply multiparticle soft theorems. Non-standard field redefinitions are showcased using both relativistic and nonrelativistic examples.

effective field theories↗
DOE OSTI2024

Thermal quasiparticle theory

The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Foundations of variational discrete action theory

Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Analytical continuation of matrix-valued functions: Carathéodory formalism

Finite-temperature quantum field theories are formulated in terms of Green's functions and self-energies on the Matsubara axis. In multiorbital systems, these quantities are related to positive semidefinite matrix-valued functions of the Carathéodory and Schur class. Analysis, interpretation, and evaluation of derived quantities such as real-frequency response functions requires analytic continuation of the off-diagonal elements to the real axis. Here we derive the criteria under which such functions exist for given Matsubara data and present an interpolation algorithm that intrinsically respects their mathematical properties. For small systems with precise Matsubara data, we find that the continuation exactly recovers all off-diagonal and diagonal elements. In real-materials systems, we show that the precision of the continuation is sufficient for the analytic continuation to commute with the Dyson equation, and we show that the truncation of the off-diagonal self-energy elements leads to considerable approximation artifacts.

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

Green's functions applied to the theory of spectroscopy

Green’s functions are a powerful analytical and computational tool for ab initio calculations of X-ray spectra. For example, Green’s functions provide an efficient means for calculations over broad energy ranges since many-body effects can be incorporated naturally in terms of the electron self-energy. Here, their role in the theory of X-ray absorption and related spectroscopies is discussed, with particular focus on many-body effects such as quasi-particle energy shifts and lifetimes, core-hole interactions and electron–phonon interactions. Additionally, the cumulant expansion for the Green’s function is reviewed and compared with the usual GW Dyson equation approach for including multielectron excitations that produce satellite features in X-ray spectra.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗