Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Greens function”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Exploring connections between statistical mechanics and Green’s functions for realistic systems: Temperature dependent electronic entropy and internal energy from a self-consistent second-order Green’s function

Including finite-temperature effects from the electronic degrees of freedom in electronic structure calculations of semiconductors and metals is desired; however, in practice it remains exceedingly difficult when using zero-temperature methods, since these methods require an explicit evaluation of multiple excited states in order to account for any finite-temperature effects. Using a Matsubara Green’s function formalism remains a viable alternative, since in this formalism it is easier to include thermal effects and to connect the dynamic quantities such as the self-energy with static thermodynamic quantities such as the Helmholtz energy, entropy, and internal energy. However, despite the promising properties of this formalism, little is known about the multiple solutions of the non-linear equations present in the self-consistent Matsubara formalism and only a few cases involving a full Coulomb Hamiltonian were investigated in the past. In this study, to shed some light onto the iterative nature of the Green’s function solutions, we self-consistently evaluate the thermodynamic quantities for a one-dimensional (1D) hydrogen solid at various interatomic separations and temperatures using the self-energy approximated to second-order (GF2). At many points in the phase diagram of this system, multiple phases such as a metal and an insulator exist, and we are able to determine the most stable phase from the analysis of Helmholtz energies. Additionally, we show the evolution of the spectrum of 1D boron nitride to demonstrate that GF2 is capable of qualitatively describing the temperature effects influencing the size of the band gap.

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

Communication: Energy-dependent resonance broadening in symmetric and asymmetric molecular junctions from an ab initio non-equilibrium Green's function approach

The electronic structure of organic-inorganic interfaces often features resonances originating from discrete molecular orbitals coupled to continuum lead states. An example is molecular junction, individual molecules bridging electrodes, where the shape and peak energy of such resonances dictate junction conductance, thermopower, I-V characteristics, and related transport properties. In molecular junctions where off-resonance coherent tunneling dominates transport, resonance peaks in the transmission function are often assumed to be Lorentzian functions with an energy-independent broadening parameter Γ. Here we define a new energy-dependent resonance broadening function, Γ(E), based on diagonalization of non-Hermitian matrices, which can describe resonances of a more complex, non-Lorentzian nature and can be decomposed into components associated with the left and right leads, respectively. We compute this quantity via an ab initio non-equilibrium Green's function (NEGF) approach based on density functional theory (DFT) for both symmetric and asymmetric molecular junctions, and show that our definition of Γ(E), when combined with Breit-Wigner formula, reproduces the transmission calculated from DFT-NEGF. Through a series of examples, we illustrate how this approach can shed new light on experiments and understanding of junction transport properties in terms of molecular orbitals.

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

GFCCLib: Scalable and efficient coupled-cluster Green's function library for accurately tackling many-body electronic structure problems

Coupled-cluster Green’s function (GFCC) calculation has drawn much attention in the recent years for targeting the molecular and material electronic structure problems from a many-body perspective in a systematically improvable way. However, GFCC calculations on scientific computing clusters usually suffer from expensive higher di- mensional tensor contractions in the complex space, expensive inter-process communi- cation, and severe load imbalance, which limits it’s routine use for tackling electronic structure problems. Here we present a numerical library prototype that is specifically designed for large-scale GFCC calculations. The design of the library is focused on a systematically optimal computing strategy to improve its scalability and efficiency. The performance of the library is demonstrated by the relevant profiling analysis of running GFCC calculations on remote giant computing clusters. The capability of the library is highlighted by computing a wide near valence band of a fullerene C60 molecule for the first time at the GFCCSD level that shows excellent agreement with the experimental spectrum.

Peng, Bo↗
DOE OSTI2019

Monte Carlo Second- and Third-Order Many-Body Green’s Function Methods with Frequency-Dependent, Nondiagonal Self-Energy

We fully develop the Monte Carlo many-body Green’s function (MC-GF) method with the following enhancements: (1) The truncation order of the perturbation expansion of the Dyson self-energy is raised from the second order (MC-GF2) to the third order (MC-GF3) with the aid of a computerized procedure to enumerate and transform all 84 third-order Goldstone diagrams into Monte Carlo integrable expressions and then into central processing unit (CPU)/graphical processing unit (GPU)-parallel computer codes. (2) An efficient algorithm is proposed that computes all off-diagonal and diagonal elements of the MC-GF2 and MC-GF3 self-energy matrices by common subexpression elimination. (3) The frequency-independent approximation is lifted by introducing a method that computes frequency derivatives of the MC-GF2 and MC-GF3 self-energies up to any arbitrarily high order at nearly no additional computational cost. (4) The imaginary-time integration in the Laplace-transformed expressions of the self-energy is carried out stochastically (instead of using a quadrature in the previous implementations), resulting in a 50- to 200-fold speedup. (5) The efficiency of the redundant-walker convergence acceleration scheme is analyzed numerically, and the guidelines are established to select an optimal number of walkers for maximal efficiency. When such an optimal number is used, the cost per sample is constant of molecular size on either many CPUs or many GPUs. (6) The computational cost to obtain a binding energy within a given statistical uncertainty is observed to increase as (tentatively) O(n 4 ) and O(n 5 ) of molecular size (n is the number of electrons) with and without the diagonal approximation, respectively, as opposed to O(n 5 ) and O(n 6 ) of the corresponding deterministic algorithms. With this method applied to the electron binding energies of C 60 , we show that the third-order corrections to the self-energies are much greater in electron binding energies than in ground-state energies. As a result, they display a sign of oscillatory convergence toward experimental results, not necessarily improving the agreement with increasing perturbation order, justifying MC-GF3 and motivating even higher-order methods.

Algorithms↗
DOE OSTI2017

Tartarus: A relativistic Green's function quantum average atom code

A relativistic Green’s Function quantum average atom model is implemented in the Tartarus code for the calculation of equation of state data in dense plasmas. We first present the relativistic extension of the quantum Green’s Function average atom model described by Starrett [1]. The Green’s Function approach addresses the numerical challenges arising from resonances in the continuum density of states without the need for resonance tracking algorithms or adaptive meshes, though there are still numerical challenges inherent to this algorithm. We discuss how these challenges are addressed in the Tartarus algorithm. The outputs of the calculation are shown in comparison to PIMC/DFT-MD simulations of the Principal Shock Hugoniot in Silicon. Finally, we also present the calculation of the Hugoniot for Silver coming from both the relativistic and nonrelativistic modes of the Tartarus code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
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 OSTI2024

Green's function methods for multiphysics simulations (Final Report)

Green's functions are important tools for analyzing mathematical properties of partial differential equations (PDEs), and for numerically solving PDEs, especially when equations for the same operator but with multiple right hand sides need to be solved simultaneously. This proposal aims at developing efficient and accurate numerical methods for computing Green's functions, which can be used to tackle a challenging question in multiphysics simulation of DOE-mission science problems: how to couple quantum physics with classical physics. The key mathematical difficulty is properly formulate a "boundary condition" for the region described by quantum physics, and conventional approaches often model such boundary conditions in an empirical way. The proposed Green's function methods use the Dirichlet-to-Neumann map to formulate a boundary condition that is non-empirical and can couple the quantum and classical regions in an in principle exact way. The key ingredient of the new methods is to construct the Dirichlet-to-Neumann map in an efficient, accurate and versatile manner. The new methods have provably low complexity and are ideally suited for massively parallel and emerging many-core computational systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Theory of Green’s Functions of Nambu-Goldstone and Higgs Modes in Superconductors

We examine fundamental properties of Green’s functions of Nambu-Goldstone and Higgs modes in superconductors with multiple order parameters. Nambu-Goldstone and Higgs modes are determined once the symmetry of the system and that of the order parameters are specified. Multiple Nambu-Goldstone modes and Higgs modes exist when we have multiple order parameters. The Nambu-Goldstone Green function D(ω, q) has the form 1/(gN(0)){sup 2}⋅(2Δ){sup 2}/(ω{sup 2}−c{sub s}{sup 2}q{sup 2}) with the coupling constant g and c{sub s}=v{sub F}/√3 for small ω and q, with a pole at ω = 0 and q = 0 indicating the existence of a massless mode. It is shown, based on the Ward-Takahashi identity, that the massless mode remains massless in the presence of intraband scattering due to non-magnetic and magnetic impurities. The pole of D(ω, q), however, disappears as ω increases as large as 2Δ: ω∼2Δ. The Green function H(ω, q) of the Higgs mode is given by H(ω,q)∝(2Δ){sup 2}/((2Δ){sup 2}−(1/3)ω{sup 2}+(1/3)c{sub s}{sup 2}q{sup 2}) for small ω and q, and H(ω, q) is proportional to 1/(gN(0)){sup 2}⋅Δ/√((2Δ){sup 2}+c{sub s}{sup 2}q{sup 2}−ω{sup 2}) for ω∼2Δ with the cut in the interval − ω(q) ≤ z ≡ q{sub 0} ≤ ω(q) where ω(q)=√((2Δ){sup 2}+c{sub s}{sup 2}q{sup 2}). The constant part of the action for the Higgs modes is important since it determines the coherence length of a superconductor. There is the case that it has a large eigenvalue, indicating that the large upper critical field H{sub c2} may be realized in a superconductor with multiple order parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Gauge invariant quark Green’s functions with polygonal Wilson lines

Properties of gauge invariant two-point quark Green’s functions, defined with polygonal Wilson lines, are studied. The Green’s functions can be classified according to the number of straight line segments their polygonal lines contain. Functional relations are established between the Green’s functions with different numbers of segments on the polygonal lines. An integrodifferential equation is obtained for the Green’s function with one straight line segment, in which the kernels are represented by a series of Wilson loop vacuum averages along polygonal contours with an increasing number of segments and functional derivatives on them. The equation is exactly solved in the case of two-dimensional QCD in the large-N{sub c} limit. The spectral properties of the Green’s function are displayed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Functional renormalization-group approaches, one-particle (irreducible) reducible with respect to local Green’s functions, with dynamical mean-field theory as a starting point

We consider formulations of the functional renormalization-group (fRG) flow for correlated electronic systems with the dynamical mean-field theory as a starting point. We classify the corresponding renormalization-group schemes into those neglecting one-particle irreducible six-point vertices (with respect to the local Green’s functions) and neglecting one-particle reducible six-point vertices. The former class is represented by the recently introduced DMF{sup 2}RG approach [31], but also by the scale-dependent generalization of the one-particle irreducible representation (with respect to local Green’s functions, 1PI-LGF) of the generating functional [20]. The second class is represented by the fRG flow within the dual fermion approach [16, 32]. We compare formulations of the fRG approach in each of these cases and suggest their further application to study 2D systems within the Hubbard model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Efficient three-dimensional Poisson solvers in open rectangular conducting pipe

Three-dimensional (3D) Poisson solver plays an important role in the study of space-charge effects on charged particle beam dynamics in particle accelerators. In this paper, we propose three new 3D Poisson solvers for a charged particle beam in an open rectangular conducting pipe. These three solvers include a spectral integrated Green function (IGF) solver, a 3D spectral solver, and a 3D integrated Green function solver. These solvers effectively handle the longitudinal open boundary condition using a finite computational domain that contains the beam itself. This saves the computational cost of using an extra larger longitudinal domain in order to set up an appropriate finite boundary condition. Using an integrated Green function also avoids the need to resolve rapid variation of the Green function inside the beam. The numerical operational cost of the spectral IGF solver and the 3D IGF solver scales as , where is the number of grid points. The cost of the 3D spectral solver scales as , where is the maximum longitudinal mode number. We compare these three solvers using several numerical examples and discuss the advantageous regime of each solver in the physical application.

Computer Science↗
DOE OSTI2024

The Green's Function Model Intercomparison Project (GFMIP) Protocol

The atmospheric Green's function method is a technique for modeling the response of the atmosphere to changes in the spatial field of surface temperature. While early studies applied this method to changes in atmospheric circulation, it has also become an important tool to understand changes in radiative feedbacks due to evolving patterns of warming, a phenomenon called the “pattern effect.” To better study this method, this paper presents a protocol for creating atmospheric Green's functions to serve as the basis for a model intercomparison project, GFMIP. The protocol has been developed using a series of sensitivity tests performed with the HadAM3 atmosphere-only general circulation model, along with existing and new simulations from other models. Our preliminary results have uncovered nonlinearities in the response of the atmosphere to surface temperature changes, including an asymmetrical response to warming versus cooling patch perturbations, and a change in the dependence of the response on the magnitude and size of the patches. These nonlinearities suggest that the pattern effect may depend on the heterogeneity of warming as well as its location. These experiments have also revealed tradeoffs in experimental design between patch size, perturbation strength, and the length of control and patch simulations. The protocol chosen on the basis of these experiments balances scientific utility with the simulation time and setup required by the Green's function approach. Running these simulations will further our understanding of many aspects of atmospheric response, from the pattern effect and radiative feedbacks to changes in circulation, cloudiness, and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Multiscale Neural Networks for Approximating Green’s Functions

Neural networks (NNs) have been widely used to solve partial differential equations (PDEs) in the applications of physics, biology, and engineering. One effective approach for solving PDEs with a fixed differential operator is learning Green’s functions. However, Green’s functions are notoriously difficult to learn due to their poor regularity, which typically requires larger NNs and longer training times. In this work, we address these challenges by leveraging multiscale NNs to learn Green’s functions. Through theoretical analysis using multiscale Barron space methods and experimental validation, we show that the multiscale approach significantly reduces the necessary NN size and accelerates training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Coupled cluster Green's function: Past, Present, and Future

Coupled cluster Green's function (CCGF) approach has drawn much attention in recent years for targeting the molecular and material electronic structure problems from a many-body perspective in a systematically improvable way. Here, we will present a brief review of the history of how the Green's function method evolved with the wavefunction, early and recent development of CCGF theory, and more recently scalable CCGF software development. We will highlight some of the recent applications of CCGF approach and propose some potential applications that would emerge in the near future.

Peng, Bo↗