Search NASA⌕ Search

SEARCH · Search NASA

Results for “Green's function methods”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
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 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 OSTI2020

Evaluating the exfoliation of two-dimensional materials with a Green's function surface model

Previous methods for the evaluation of the exfoliation of two-dimensional (2D) layered materials have drawbacks in computational efficiency and are unable to describe cases with semi-infinite substrates. Based on a Green's function surface (GFS) model, here we develop an approach to efficiently determine the tendency of exfoliation of 2D materials from their bulk crystals or semi-infinite substrates. By constructing appropriate surface configurations, we may calculate the exfoliation energy more precisely and quickly than the traditional way with the slab model. Furthermore, the GFS approach can provide angle-resolved photoemission spectroscopy of surface systems for direct comparison with experimental data. In this work, our findings indicate that the GFS approach is powerful for studies of 2D materials and various surface problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1987

Full CI benchmark calculations for several states of the same symmetry

Full CI (FCI) wave functions are used to compute energies for several electronic states of the same symmetry for SiH2, CH2, and CH2(+). It is found that CASSCF/multireference CI wave functions yield results very similar to FCI, irrespective of whether the CASSCF MOs are optimized independently for each state or using an average of the CASSCF energies for all desired states. The ionization potentials and excitation energies obtained from the FCI calculations should help calibrate methods (such as Green's function approaches, equations of motion and propagator methods, and cluster expansions) in which energy differences are computed directly.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2022

Analysis of spin frustration in an Fe III 7 cluster using a combination of computational, experimental, and magnetostructural correlation methods

The synthesis, structure, and magnetic properties are reported for [Fe 7 O 3 (O 2 C t Bu) 9 (mda) 3 (H 2 O) 3 ] ( 1 ), where mdaH 2 is N -methyldiethanolamine. 1 was prepared from the reaction of [Fe 3 O(O 2 C t Bu) 6 (H 2 O) 3 ](NO 3 ) with mdaH 2 in a 1:~3 ratio in MeCN. The core of 1 consists of a central octahedral Fe III ion held within a non-planar Fe 6 loop by three μ 3 -O 2- and three μ 2 -RO - arms from the three mda 2- chelates. Variable-temperature dc and ac magnetic susceptibility studies revealed dominant antiferromagnetic coupling, leading to a ground state spin of S = 5 / 2 . The ground state was confirmed by a fit of magnetization data collected in the 0.1–7.0 T and 1.8–10.0 K ranges. The four Fe 2 pairwise exchange parameters ( J 1 - J 4 ) were estimated by independent methods: theoretical calculations using either broken symmetry energy differences (-46.3, -16.2, -3.9, and - 28.1 cm -1 , respectively) or Green’s function approximation methods (-41.4, -14.8, -13.2, and - 24.7 cm -1 ), and a magnetostructural correlation (MSC) previously developed for high nuclearity Fe III /O complexes (-39.5, -13.8, -6.7, and - 23.5 cm -1 ). Additionally, the J 1 - J 4 obtained from the MSC and theoretical methods were used with the program PHI to both simulate χ M T vs T as well as to serve as reasonable input values to fit the experimental data (-41.0, -11.4, -5.0, and - 27.3 cm -1 ). Analysis of the J ij led to identification of the spin frustration effects operative and the resultant spin vector alignments at each Fe III ion, thus allowing for the rationalization of the experimental ground state.

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

Electric fields near undulating dielectric membranes

Dielectric interfaces are crucial to the behavior of charged membranes, from graphene to synthetic and biological lipid bilayers. Understanding electrolyte behavior near these interfaces remains a challenge, especially in the case of rough dielectric surfaces. A lack of analytical solutions consigns this problem to numerical treatments. Here, we report an analytic method for determining electrostatic potentials near curved dielectric membranes in a two-dimensional periodic “slab” geometry using a periodic summation of Green’s functions. This method is amenable to simulating arbitrary groups of charges near surfaces with two-dimensional deformations. We concentrate on one-dimensional undulations. We show that increasing membrane undulation increases the asymmetry of interfacial charge distributions due to preferential ionic repulsion from troughs. In the limit of thick membranes, we recover results mimicking those for electrolytes near a single interface. Our work demonstrates that rough surfaces generate charge patterns in electrolytes of charged molecules or mixed-valence ions.

Chemistry↗
NASA NTRS1980

Note on the eigensolution of a homogeneous equation with semi-infinite domain

The 'variation-iteration' method using Green's functions to find the eigenvalues and the corresponding eigenfunctions of a homogeneous Fredholm integral equation is employed for the stability analysis of fluid hydromechanics problems with a semiinfinite (infinite) domain of application. The objective of the study is to develop a suitable numerical approach to the solution of such equations in order to better understand the full set of equations for 'real-world' flow models. The study involves a search for a suitable value of the length of the domain which is a fair finite approximation to infinity, which makes the eigensolution an approximation dependent on the length of the interval chosen. In the examples investigated y = 1 = a seems to be the best approximation of infinity; for y greater than unity this method fails due to the polynomial nature of Green's functions.

Wadia, A. R.↗
DOE OSTI2021

Fast inversion, preconditioned quantum linear system solvers, fast Green's-function computation, and fast evaluation of matrix functions

Preconditioning is the most widely used and effective way for treating ill-conditioned linear systems in the context of classical iterative linear system solvers. We introduce a quantum primitive called fast inversion, which can be used as a preconditioner for solving quantum linear systems. The key idea of fast inversion is to directly block encode a matrix inverse through a quantum circuit implementing the inversion of eigenvalues via classical arithmetics. We demonstrate the application of preconditioned linear system solvers for computing single-particle Green's functions of quantum many-body systems, which are widely used in quantum physics, chemistry, and materials science. We analyze the complexities in three scenarios: the Hubbard model, the quantum many-body Hamiltonian in the plane-wave-dual basis, and the Schwinger model. We also provide a method for performing Green's function calculation in second quantization within a fixed-particle manifold and note that this approach may be valuable for simulation more broadly. Aside from solving linear systems, fast inversion also allows us to develop fast algorithms for computing matrix functions, such as the efficient preparation of Gibbs states. Furthermore, we introduce two efficient approaches for such a task, based on the contour-integral formulation and the inverse transform, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Dynamic, symmetry-preserving, and hardware-adaptable circuits for quantum computing many-body states and correlators of the Anderson impurity model

We present a hardware-reconfigurable ansatz on N q -qubits for the variational preparation of many-body states of the Anderson impurity model (AIM) with N imp + N bath = N q /2 sites, which conserves total charge and spin z component within each variational search subspace. The many-body ground state of the AIM is determined as the minimum over all minima of O(N$^2_ q$) distinct charge-spin sectors. Hamiltonian expectation values are shown to require ω(N q ) < N meas. $\leqslant$ O(N imp N bath ) symmetry-preserving, parallelizable measurement circuits, each amenable to postselection. To obtain the one-particle impurity Green’s function we show how initial Krylov vectors can be computed via midcircuit measurement and how Lanczos iterations can be computed using the symmetry-preserving ansatz. For a single-impurity Anderson model with a number of bath sites increasing from one to seven, we show using numerical emulation that the ease of variational ground-state preparation is suggestive of linear scaling in circuit depth and subquartic scaling in optimizer complexity. We therefore expect that, combined with time-dependent methods for Green’s function computation, our ansatz provides a useful tool to account for electronic correlations on early fault-tolerant processors. Finally, with a view towards computing real materials properties of interest like magnetic susceptibilities and electron-hole propagators, we provide a straightforward method to compute many-body, time-dependent correlation functions using a combination of time evolution, midcircuit measurement-conditioned operations, and the Hadamard test.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomistic modeling of energy band alignment in CdSeTe surfaces

Atomistic modeling based on Desnity Functional Theory (DFT)- 1/2 method coupled with surface Green’s function has been employed to investigate the energy band alignment results in cadmium selenium telluride (CdSeTe) surfaces. The structural and electronic properties of the bulk ternary alloy were established before cleaving the CdSeTe low index surface facets. The dependency of atomic-scale electronic features on different plane orientations was explored by studying the energy band alignment diagrams in the unreconstructed and reconstructed surfaces. While the low index CdSeTe surface geometry reconstructions replicate the CdTe surface geometries, the energy band alignment features of the unreconstructed and reconstructed CdSeTe surfaces differ from those observed in CdTe low index surfaces. The structurally relaxed reconstructed CdSeTe surface results in purely flat bands as opposed to an internal cusp feature observed in CdTe surfaces. The presence of an internal energy cusp in CdTe(111) surfaces is expected to play a key role in enhancing the hole charge transport towards the back of CdTe-only solar cell device and the absence of such feature in CdSeTe surface may explain one of the reason for the lower performance of CdSeTe-only solar cell seen experimentally.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Scattering from perfectly conducting and resistive strips on a grounded dielectric slab

The scattering properties of perfectly conducting and resistive strips are predicted for strips which are located on a dielectric slab backed by a perfectly conducting ground plane. The spectral domain Green's function is used to relate the currents and fields on the strip, and the resulting integral equation is solved using the method of moments. Both TE and TM strips are examined using piecewise linear and pulse subdomain basis functions, respectively, to model the current on the strip. Calculated results are compared with results measured at the NASA Langley Research Center.

PULSE SUBDOMAIN BASIS FUNCTION↗
DOE OSTI2023

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

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

Space-Charge Aberrations in Single-Shot Time-Resolved Transmission Electron Microscopy

In this paper, we discuss effects of space-charge fields on imaging performance in a single-shot time-resolved transmission electron microscope. Using a Green’s function perturbation method, we derive analytical estimates for the effects of space-charge nonlinearity on the image formation process and the associated aberration coefficients. The dependence of these coefficients on the initial beam phase space distribution is elucidated. The results are validated by particle tracking simulations and provide fundamental scaling laws for the trade-off between temporal and spatial resolutions in single-shot time-resolved transmission electron microscopy.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Effect of Time Window and Spectral Measurement Options on Empirical Green’s Function Analysis Using DAS Array and Seismic Stations

The recorded seismic waveform is a convolution of event source term, path term, and station term. Removing high-frequency attenuation due to path effect is a challenging problem. Empirical Green’s function (EGF) method uses nearly collocated small earthquakes to correct the path and station terms for larger events recorded at the same station. However, this method is subject to variability due to many factors. Here, we focus on three events that were well recorded by the seismic network and a rapid response distributed acoustic sensing (DAS) array. Using a suite of high-quality EGF events, we assess the influence of time window, spectral measurement options, and types of data on the spectral ratio and relative source time function (RSTF) results. Increased number of tapers (from 2 to 16) tends to increase the measured corner frequency and reduce the source complexity. Extended long time window (e.g., 30 s) tends to produce larger variability of corner frequency. The multitaper algorithm that simultaneously optimizes both target and EGF spectra produces the most stable corner-frequency measurements. The stacked spectral ratio and RSTF from the DAS array are more stable than two nearby seismic stations, and are comparable to stacked results from the seismic network, suggesting that DAS array has strong potential in source characterization.

58 GEOSCIENCES↗
NASA NTRS1983

Acoustic loading in planar networks

The acoustic loading in a complex planar network of ducts is determined by a method in which Green function surface elements are used. The network consists of straight ducts, elbows and branched ducts. A transfer matrix technique is developed in which each duct is treated separately and the matrix of the influence coefficients is transformed to tri-diagonal form allowing efficient inversion.

El-Raheb, M.↗
NASA NTRS1984

Sensitivity analysis of discrete structural systems: A survey

Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications addressed to structural analysis, but also includes a number of methods developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.

Adelman, H. M.↗
NASA NTRS1984

Measurements of large scale disturbances in rectangular wall jets

Detailed measurements of the velocity field and wall surface pressure distribution of a tone-generating wall jet were obtained in the Ames 7- by 10-foot Wind Tunnel and anechoeic chamber. The jet issued from a 0.508 cm by 10.16 cm nozzle adjacent to a 3.81 cm wall, with an initially parabolic velocity profile at the exit. The jet was operated at Reynolds numbers, R, of 4910, 7740, and 9410. At R = 4910 and 7740, the flow appeared as a double-rowed vortex array which formed upstream of the edge and convected downstream without change in structure. At R = 9410, apparent vortex pairing occurred at the edge location. These flow characteristics were observed in both velocity field measurements and flow visualization records. The unsteady vorticity field of the jet was computed from phase-averaged hot-wire anemometer measurements of the velocity field. The associated acoustic field was estimated by the method of Green's functions, and was found to be similar to measurements of the sound field. These observations tend to support theoretical predictions that vortex shedding at the edge tends to cancel the effect of acoustic sources in other parts of the flow field.

Horne, W. C.↗
NASA NTRS1986

Sensitivity analysis of discrete structural systems

Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.

Adelman, H. M.↗