Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “green 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 505 records · Page 28

Effect of Flight on the Noise from Turbulent Jets in the Generalized Acoustic Analogy

This presents an extension of a noise prediction method for static jets based on the Generalized Acoustic Analogy to include a non-zero ambient stream to simulate the effects of flight. Extended formulae for the propagator functions, source terms and the far-field acoustic spectrum are derived. The formulations of the equation and the boundary conditions for numerical solution of the scalar adjoint Green’s function in the presence of a non-zero ambient stream are also given. Wherever possible, the corresponding equations in the original papers for the static case and their extensions to include a flight stream are identified. The extended formulation is implemented into an existing code and used to make predictions of the effect of the flight stream on the observed far-field noise for a series of subsonic, unheated, axisymmetric jets over a range of jet-exit and free-stream velocities using Reynolds-averaged Navier-Stokes flow solutions as input. The predictions are compared with data from experiments conducted at NASA Glenn Research Center which extend a fundamental jet noise database to include a flight stream of various speeds. Descriptions of the experimental setup, data acquisition, post-processing and corrections applied to account for the external shear-layer are given. It is shown that the acoustic analogy-based method can provide reasonably good predictions of the impact of an external flight stream on the spectral characteristics of noise from unheated round jets.

Jet noise↗

The coordination chemistry of oxide and nanocarbon materials

Understanding how a ligand affects the steric and electronic properties of a metal is the cornerstone of the inorganic chemistry enterprise. What happens when the ligand is an extended surface? This question is central to the design and implementation of state-of-the-art functional materials containing transition metals. Here, this perspective will describe how these two very different sets of extended surfaces can form well-defined coordination complexes with metals. In the Green formalism, functionalities on oxide surfaces react with inorganics to form species that contain X-type or LX-type interactions between the metal and the oxide. Carbon surfaces are neutral L-type ligands; this perspective focuses on carbons that donate six electrons to a metal. The nature of this interaction depends on the curvature, and thereby orbital overlap, between the metal and the extended π-system from the nanocarbon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Equation of state for 304L stainless steel

Here, a new equation of state for 304L stainless steel has been constructed and entered in the Los Alamos National Laboratory SESAME database with material ID 4273. The equation of state is based on a multiphase approach. The solid and liquid phases were fitted to available experimental data for thermal expansion, heat capacity, and shock Hugoniot, as well as to data from density functional theory calculations of the cold curve and quantum molecular dynamics of the liquid regime. The electronic thermal contribution over the entire equation of state range is based on Tartarus (Green’s functions) average atom calculations. For this low-carbon stainless steel both the density functional and electronic thermal calculations are comprised of 19 wt% Cr, 10 wt% Ni, and the remainder Fe.

36 MATERIALS SCIENCE↗

Validating first-principles phonon lifetimes via inelastic neutron scattering

Phonon lifetimes are a key component of quasiparticle theories of transport; yet first-principles lifetimes are rarely directly compared with inelastic neutron scattering (INS) results. Existing comparisons show discrepancies even at temperatures where perturbation theory is expected to be reliable. In this paper, we demonstrate that the reciprocal space voxel (q voxel), which is the finite region in reciprocal space required in INS data analysis, must be explicitly accounted for within theory in order to draw a meaningful comparison. Here, we demonstrate accurate predictions of peak widths of the scattering function when accounting for the q voxel in CaF 2 and ThO 2 . Passing this test implies high fidelity of the phonon interactions and the approximations used to compute the Green's function, serving as a critical benchmark of theory and indicating that other material properties should be accurately predicted, which we demonstrate for thermal conductivity.

36 MATERIALS SCIENCE↗

Explosive Yield Estimation Using Regional Seismic Moment Tensors

Here, we use the Pasyanos and Chiang (2022) data set to calculate the seismic moment M 0 for each explosion and use the measured explosive yield W to validate the W~M 0 relationship in Denny and Johnson (1991; hereafter, DJ91). The M 0 is corrected by transforming to a potency tensor and applying more appropriate near-source geophysical parameter values in the moment estimate. The mean residual between observed and predicted yield is near zero; however, the standard deviation of the residuals results in an F-value (a 95% confidence factor) of about 5. We re-estimate the coefficients in the DJ91 model and find similar values and only a slight improvement in the F-value. Next, we embark on a similar model selection process as DJ91, allowing for non-cube-root yield scaling and other plausible near-source elastic moduli. As was found by DJ91, the yield dependence is not significantly different from unity, and a cube root assumption is valid. Therefore, we yield scale the seismic moment and test the significance of all plausible explanatory variables. Isotropic moment performs better in the response variable than total moment. The preference for isotropic moment could be due to its relationship to volume change, which would be more directly affected by explosive yield. Surprisingly, we find that the overburden pressure, which is a function of depth, is not a significant parameter in the model. We hypothesize that this is due to the competing depth effects on source asymmetry and the incorporation of depth in the Green’s functions used to calculate the seismic moment tensors. Importantly, this emphasizes that only seismic moment tensor-derived moments should be used in these models. After removing insignificant model parameters, we are left with a simple model to predict explosive yield $\widehat{W}$ in kt from isotropic moment M I in N·m, $\widehat{W}$=κ –1.4132 10 0.035626GP M I , in which κ and GP are the near-source bulk modulus and gas porosity in Pa and %, respectively. The F-value for this model is approximately 3.

58 GEOSCIENCES↗

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↗

Nanoscale functionalized superconducting transport channels as photon detectors

Herein, single-photon detectors have typically consisted of macroscopic materials where both the photon absorption and transduction to an electrical signal happen. Newly proposed designs suggest that large arrays of nanoscale detectors could provide improved performance in addition to decoupling the absorption and transduction processes. Here we study the properties of such a detector consisting of a nanoscale superconducting (SC) transport channel functionalized by a photon absorber. We explore two detection mechanisms based on photoinduced electrostatic gating and magnetic effects. To this end we model the narrow channel as a one-dimensional atomic chain and use a self-consistent Keldysh-Nambu Green's function formalism to describe nonequilibrium effects and SC phenomena. We consider cases where the photon creates electrostatic and magnetic changes in the absorber, as well as devices with strong and weak coupling to the metal leads. Our results indicate that the most promising case is when the SC channel is weakly coupled to the leads and in the presence of a background magnetic field, where photoexcitation of a magnetic molecule can trigger a SC-to-normal transition in the channel that leads to a change in the device current several times larger than in the case of a normal-phase channel device.

1-dimensional systems↗

Load-Displacement Relation and Gap Distribution Between Rough Surfaces: Partial Differential Equations Approach

We develop a theoretical model to predict the load-displacement relation and probability density function for gaps between contacting rough surfaces. We derive partial differential equations from the previous model (Joe et al., 2018), and extend them to the non-adhesive contact problem. The predictions of the present theory are compared with numerical results using the Green's function molecular dynamics algorithm, and a good agreement is obtained.

adhesion↗

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Linear response to nonstationary random excitation.

Development of a method for computing the mean-square response of linear systems to nonstationary random excitation of the form given by y(t) = f(t) x(t), in which x(t) = a stationary process and f(t) is deterministic. The method is suitable for application to multidegree-of-freedom systems when the mean-square response at a point due to excitation applied at another point is desired. Both the stationary process, x(t), and the modulating function, f(t), may be arbitrary. The method utilizes a fundamental component of transient response dependent only on x(t) and the system, and independent of f(t) to synthesize the total response. The role played by this component is analogous to that played by the Green's function or impulse response function in the convolution integral.

Hasselman, T.↗

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.↗

Static shape determination and control for large space structures. I - The flexible beam. II - A large space antenna

A method for determining and controlling the shape of large, continuous space structures by means of discrete or pointwise observations and control devices is presented. The general linear boundary value problem satisfied by a one-dimensional shape function is defined, and the existence of solutions is studied. The static shape control problems for one-dimensional systems with and without rigid body modes and the static shape estimation problem are presented and solved. Eigenfunction expansions are presented which provide approximations to the algorithm terms when the associated Green's function is not known. An integral operator approach is applied to the multidimensional static problem, and the results are illustrated with a finite element model of the disk of a large space antenna which assumes no rigid body modes. It is shown that the shape control algorithm must be modified for systems with rigid body modes.

Weeks, C. J.↗

Design Rules for Open-Shell Molecular Wires: Insights from Correlated Many-Body Transport

Open-shell π-conjugated systems have emerged as promising molecular wires that can sustain unusually high low-bias conductance over tens of nanometers. However, predictive design rules remain limited because near-degeneracy and spin polarization give rise to correlated electronic states that are not reliably captured by standard density functional theory and Landauer transport descriptions. Here, we establish physically transparent design rules for long-range, weakly length-dependent charge transport based on the interplay between bond-length alternation and electron–electron interactions. Using a fully correlated transport framework that combines density matrix renormalization group calculations with nonequilibrium Green’s function embedding, we show that weak dimerization maintains spatially extended edge states, while moderate on-site interaction stabilizes open-shell character without excessively separating transport-relevant resonances. This balance yields zero-bias transmission that is largely insensitive to molecular length, providing a route to connect chemically tunable structure to correlated transport beyond mean-field descriptions.

Charge transport↗

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↗

Reduced Artifact Approach for Determining Diffusion Coefficients in Time-Resolved Microscopy

Ultrafast microscopy methods traditionally assume a Gaussian profile to extract excited state diffusivities from transport measurements. Although this fitting method recovers accurate diffusion coefficients when the point spread function is well-represented by a Gaussian, even minor spatial aberrations introduced by the imaging system cause significant errors in the determined value. To provide a more accurate measure of excited state transport in nano- and microscale materials systems, in this work an alternative analysis protocol is proposed that numerically convolves the Green’s function solution to the diffusion equation with the experimentally measured point spread function. In contrast to the Gaussian fitting approach, the numerical convolution is shown to be robust against artifacts caused by nonideal point spread functions. Furthermore, the numerical convolution approach is highly effective at resolving anisotropic diffusion in modeled data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Crucial role of thermal fluctuations and vertex corrections for the magnetic pseudogap

It is generally believed that in a two-dimensional metal, whose ground state is antiferromagnetically ordered with Q = (π,π), thermal (static) magnetic fluctuations give rise to precursor behavior above TN in which the spectral function of a hot fermion (the one for which k and k + Q are Fermi surface points) contains two peaks, separated by roughly the same energy as in the antiferromagnetically ordered state. The two peaks persist in some range of T > T N and eventually merge into a single peak at zero frequency. This behavior is obtained theoretically by departing from free fermions in a paramagnet and evaluating the dressed fermionic Green's function by summing up infinite series of diagrams with contributions from thermal magnetic fluctuations. We show, following that keeping vertex renormalization diagrams in these series is crucial as other terms only broaden the spectral function of a hot fermion but do not shift its maximum away from zero frequency. As the consequence, the magnetic pseudogap should be treated as an input for theories that neglect vertex corrections, such as, e.g., Eliashberg theory for magnetically mediated superconductivity. We also analyze the potential pseudogap behavior at T = 0. Here, we argue that it may exist, but only at a finite correlation length, and not as a precursor to antiferromagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗