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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Adaptive Variational Quantum Computing Approaches for Green’s Functions and Nonlinear Susceptibilities

Here, we present and benchmark quantum computing approaches for calculating real-time single-particle Green’s functions and nonlinear susceptibilities of Hamiltonian systems. The approaches leverage adaptive variational quantum algorithms for state preparation and propagation. Using automatically generated compact circuits, the dynamical evolution is performed over sufficiently long times to achieve adequate frequency resolution of the response functions. We showcase accurate Green’s function calculations using a statevector simulator on classical hardware for Fermi-Hubbard chains of 4 and 6 sites, with maximal ansatz circuit depths of 65 and 424 layers, respectively, and for the molecule LiH with a maximal ansatz circuit depth of 81 layers. Additionally, we consider an antiferromagnetic quantum spin-1 model that incorporates the Dzyaloshinskii-Moriya interaction to illustrate calculations of the third-order nonlinear susceptibilities, which can be measured in two-dimensional coherent spectroscopy experiments. These results demonstrate that real-time approaches using adaptive parametrized circuits to evaluate linear and nonlinear response functions can be feasible with near-term quantum processors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Lattice Green’s Functions for High-Order Finite Difference Stencils

Lattice Green's Functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Elastic Green’s Function in Anisotropic Bimaterials Considering Interfacial Elasticity

Here, the two-dimensional elastic Green’s function is calculated for a general anisotropic elastic bimaterial containing a line dislocation and a concentrated force while accounting for the interfacial structure by means of a generalized interfacial elasticity paradigm. The introduction of the interface elasticity model gives rise to boundary conditions that are effectively equivalent to those of a weakly bounded interface. The equations of elastic equilibrium are solved by complex variable techniques and the method of analytical continuation. The solution is decomposed into the sum of the Green’s function corresponding to the perfectly bonded interface and a perturbation term corresponding to the complex coupling nature between the interface structure and a line dislocation/concentrated force. Such construct can be implemented into the boundary integral equations and the boundary element method for analysis of nano-layered structures and epitaxial systems where the interface structure plays an important role.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Exchange-correlation contributions to thermodynamic properties of the homogeneous electron gas from a cumulant Green's function approach

Thermodynamic properties of the interacting homogeneous electron gas are calculated using a finite-temperature cumulant Green's function approach over a broad range of densities and temperatures up to the warm dense matter regime $T ~ T_F$, where $T_F$ is the Fermi degeneracy temperature. These properties can be separated into independent particle and exchange-correlation contributions, and our focus here is on the latter. In this work our approach is based on the Galitskii-Migdal-Koltun and electron number sum rules from the finite temperature many-body Green's function formalism, together with an extension of the cumulant Green's function to finite temperature. Previously this approach yielded exchange-correlation energies and potentials in good agreement with quantum Monte-Carlo calculations. Here the method is extended for various thermodynamic quantities including the chemical potential, total energy, Helmholtz free-energy, electronic equation of state, specific heat, and isothermal compressibility, which optionally include spin dependence. We find that the exchange-correlation contributions are weakly varying at low temperature but exhibit significant temperature dependence in the WDM regime, as well as a crossover from exchange- to correlation-dominated behavior. In contrast to the $T = 0^+$ limit, we also find that renormalization effects are largely but not completely suppressed at finite temperature. Comparisons with other approaches at various levels of approximation are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Relation between two-point Green’s functions of N=1 SQED with N{sub f} flavors, regularized by higher derivatives, in the three-loop approximation

We verify the identity which relates the two-point Green’s functions of N=1 SQED with N{sub f} flavors, regularized by higher derivatives, by explicit calculations in the three-loop approximation. This identity explains why, in the limit of the vanishing external momentum, the two-point Green’s function of the gauge superfield is given by integrals of double total derivatives in the momentum space. It also makes it possible to exactly derive the NSVZ β-function in all loops if the renormalization group functions are defined in terms of the bare coupling constant. In order to verify the considered identity, we use it to construct integrals giving the three-loop β-function starting from the two-point Green’s functions of matter superfields in the two-loop approximation. Then we demonstrate that the results for these integrals coincide with the sums of the corresponding three-loop supergraphs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Properties of advanced coupled-cluster Green's function

We present an extension of the analysis previously applied to the retarded part of the coupled cluster (CC) Green's function to its advanced part. In analogy to our earlier studies for the retarded part, we demonstrate that the advanced CC Green's function is expressed in terms of connected diagrams only, which is a direct consequence of algebraic form of equations satisfied by CC amplitudes. We also demonstrate that omega-derivatives of the advanced CC Green's function can be calculated analytically and can be expressed in terms of connected diagrams only. We analyze the structure of connected diagrams and the role of intermediate operators which satisfy electron affinity equation-of-motion CC type conditions.

electronic structure theory↗
DOE OSTI2014

Non-equilibrium Green's functions method: Non-trivial and disordered leads

The non-equilibrium Green's function algorithm requires contact self-energies to model charge injection and extraction. All existing approaches assume infinitely periodic leads attached to a possibly quite complex device. This contradicts today's realistic devices in which contacts are spatially inhomogeneous, chemically disordered, and impacting the overall device characteristics. This work extends the complex absorbing potentials method for arbitrary, ideal, or non-ideal leads in atomistic tight binding representation. The algorithm is demonstrated on a Si nanowire with periodic leads, a graphene nanoribbon with trumpet shape leads, and devices with leads of randomly alloyed Si{sub 0.5}Ge{sub 0.5}. It is found that alloy randomness in the leads can reduce the predicted ON-state current of Si{sub 0.5}Ge{sub 0.5} transistors by 45% compared to conventional lead methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

One-particle many-body Green’s function theory: Algebraic recursive definitions, linked-diagram theorem, irreducible-diagram theorem, and general-order algorithms

A thorough analytical and numerical characterization of the whole perturbation series of one-particle many-body Green’s function (MBGF) theory is presented in a pedagogical manner. Three distinct but equivalent algebraic (first-quantized) recursive definitions of the perturbation series of the Green’s function are derived, which can be combined with the well-known recursion for the self-energy. Six general-order algorithms of MBGF are developed, each implementing one of the three recursions, the ΔMPn method (where n is the perturbation order) [S. Hirata et al., J. Chem. Theory Comput. 11, 1595 (2015)], the automatic generation and interpretation of diagrams, or the numerical differentiation of the exact Green’s function with a perturbation-scaled Hamiltonian. They all display the identical, nondivergent perturbation series except ΔMPn, which agrees with MBGF in the diagonal and frequency-independent approximations at 1 ≤ n ≤ 3 but converges at the full-configuration-interaction (FCI) limit at n=∞ (unless it diverges). Numerical data of the perturbation series are presented for Koopmans and non-Koopmans states to quantify the rate of convergence towards the FCI limit and the impact of the diagonal, frequency-independent, or ΔMPn approximation. The diagrammatic linkedness and thus size-consistency of the one-particle Green’s function and self-energy are demonstrated at any perturbation order on the basis of the algebraic recursions in an entirely time-independent (frequency-domain) framework. The trimming of external lines in a one-particle Green’s function to expose a self-energy diagram and the removal of reducible diagrams are also justified mathematically using the factorization theorem of Frantz and Mills. Equivalence of ΔMPn and MBGF in the diagonal and frequency-independent approximations at 1 ≤ n ≤ 3 is algebraically proven, also ascribing the differences at n = 4 to the so-called semi-reducible and linked-disconnected diagrams

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

The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

42 ENGINEERING↗
DOE OSTI2016

Description of quasiparticle and satellite properties via cumulant expansions of the retarded one-particle Green's function

A cumulant-based GW approximation for the retarded one-particle Green's function is proposed, motivated by an exact relation between the improper Dyson self-energy and the cumulant generating function. We explore qualitative aspects of this method within a simple one-electron independent phonon model, where it is seen that the method preserves the energy moment of the spectral weight while also reproducing the exact Green's function in the weak-coupling limit. For the three-dimensional electron gas, this method predicts multiple satellites at the bottom of the band, albeit with inaccurate peak spacing. But, its quasiparticle properties and correlation energies are more accurate than both previous cumulant methods and standard G0W0. These results point to features that may be exploited within the framework of cumulant-based methods and suggest promising directions for future exploration and improvements of cumulant-based GW approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Efficient numerical method for evaluating normal and anomalous time-domain equilibrium Green's functions in inhomogeneous systems

In this work we develop EPOCH (equilibrium propagator by orthogonal polynomial chain), a computationally efficient method to calculate the time-dependent equilibrium Green's functions, including the anomalous Green's functions of superconductors, to capture the time evolution in large inhomogeneous systems. The EPOCH method generalizes the Chebyshev wave-packet propagation method from quantum chemistry and efficiently incorporates the Fermi-Dirac statistics that is needed for equilibrium quantum condensed matter systems. The computational cost of EPOCH scales only linearly in the system degrees of freedom, generating an extremely efficient algorithm also for very large systems. We demonstrate the power of the EPOCH method by calculating the time evolution of an excitation near a superconductor–normal metal interface in two and three dimensions, capturing transmission as well as normal and Andreev reflections.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Detecting Multipartite Entanglement Patterns Using Single-Particle Green’s Functions

Here, we present a protocol for detecting multipartite entanglement in itinerant many-body electronic systems using single-particle Green’s functions. To achieve this, we first establish a connection between the quantum Fisher information and single-particle Green’s functions by constructing a set of witness operators built out of single electron creation and destruction operators in a doubled system. This set of witness operators is indexed by a momentum k. We compute the quantum Fisher information for these witness operators and show that for thermal ensembles it can be expressed as an autoconvolution of the single-particle spectral function. We then apply our framework to a one-dimensional fermionic system to showcase its effectiveness in detecting entanglement in itinerant electron models. We observe that the detected entanglement level is sensitive to the wave vector associated with witness operator. Our protocol will permit detecting entanglement in many-body systems using scanning tunneling microscopy and angle-resolved photoemission spectroscopy, two spectroscopies that measure the single-particle Green’s function. It offers the prospect of the experimental detection of entanglement through spectroscopies beyond the established route of measuring the dynamical spin response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

A Green's function quantum average atom model

A quantum average atom model is reformulated using Green's functions. This allows integrals along the real energy axis to be deformed into the complex plane. The advantage being that sharp features such as resonances and bound states are broadened by a Lorentzian with a half-width chosen for numerical convenience. An implementation of this method therefore avoids numerically challenging resonance tracking and the search for weakly bound states, without changing the physical content or results of the model. A straightforward implementation results in up to a factor of 5 speed-up relative to an optimized orbital based code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Green’s Function of the Schrödinger Equation in the Potential Quantization Method

The problem for determining Green’s function G(r, r') for the time-independent Schrödinger equation is considered using the potential quantization method. The expression for G(r, r') is obtained in terms of the eigenfunctions and eigenvalues of this method. The equivalence of the used and traditional methods for solving the given problem has been proven.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Electronic Excitations in Complex Molecular Environments: Many-Body Green’s Functions Theory in VOTCA-XTP

Many-body Green’s functions theory within the GW approximation and the Bethe-Salpeter Equation (BSE) is implemented in the open-source VOTCA-XTP software, aiming at the calculation of electronically excited states in complex molecular environments. Based on Gaussian-type atomic orbitals and making use of resolution of identity techniques, the code is designed specifically for nonperiodic systems. Application to a small molecule reference set successfully validates the methodology and its implementation for a variety of excitation types covering an energy range from 2 to 8 eV in single molecules. Further, embedding each GW-BSE calculation into an atomistically resolved surrounding, typically obtained from Molecular Dynamics, accounts for effects originating from local fields and polarization. Using aqueous DNA as a prototypical system, different levels of electrostatic coupling between the regions in this GW-BSE/MM setup are demonstrated. Particular attention is paid to charge-transfer (CT) excitations in adenine base pairs. It is found that their energy is extremely sensitive to the specific environment and to polarization effects. The calculated redshift of the CT excitation energy compared to a nucelobase dimer treated in vacuum is of the order of 1 eV, which matches expectations from experimental data. Predicted lowest CT energies are below that of a single nucleobase excitation, indicating the possibility of an initial (fast) decay of such an UV excited state into a binucleobase CT exciton. Furthermore, the results show that VOTCA-XTP’s GW-BSE/MM is a powerful tool to study a wide range of types of electronic excitations in complex molecular environments.

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

The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation

We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.

42 ENGINEERING↗
DOE OSTI2025

A Green’s function fast multipole method for computation of micromechanical fields in heterogeneous materials

Computation of micromechanical fields in heterogeneous materials is usually performed using either the finite element method or the Green’s function method based on FFTs. The finite element method allows for accurate discretization and for non-periodic boundary conditions but is computationally expensive. On the other hand, the FFT-based method is computationally efficient but requires discretization on a regular grid of hexahedral voxels. In this paper, a Green’s function method allowing for accurate discretization using tetrahedral elements and for non-periodic boundary conditions is proposed. The convolution is computed using the fast multipole method, which provides good accuracy even for low-order expansion due to the fast decay of interactions between elements. The proposed Green’s function fast multipole method is verified by comparison with analytical and FFT-based solutions. Furthermore, the computational time is analyzed and compared to the FFT-based method for non-periodic convolution. Finally, effective properties of an elastic polycrystalline microstructure containing thin intergranular cracks are computed and analyzed.

36 MATERIALS SCIENCE↗