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

Toward ab-initio nuclear theory calculations of 𝛿 𝐶

We propose a theory framework to study the isospin-symmetry breaking correction 𝛿 C in superallowed nuclear 𝛽 decays, crucial for the precise determination of |𝑉 𝑢⁢𝑑 |. Based on a general assumptions of the isovector dominance in isospin-symmetry breaking interactions, we construct a set of functions 𝐹 𝑇 𝑧 which involve nuclear matrix elements of isovector monopole operators and the nuclear Green's function. Via the functions 𝐹 𝑇 𝑧 , a connection of 𝛿 C to measurable electroweak nuclear radii is established, providing an experimental gauge of the theory accuracy of 𝛿 C . Here, we outline a strategy to perform ab initio calculations of 𝐹 𝑇 𝑧 based on the Lanczos algorithm, and discuss its similarity with other nuclear-structure-dependent inputs in nuclear 𝛽 decays.

Ab initio calculations↗
DOE OSTI2023

New Dimensions in the Theory of Excited States and X-ray Spectra (Final Report)

This Final Technical Report briefly summarizes the achievements during the lifetime of our DOE BES grant DE-FG02-97ER45623. The long-term goal of this project has been the development of quantitative theories of the interaction between radiation and matter, with a focus on x-ray spectroscopies. X-ray spectra have long been among the most important probes of atomic-scale structure and properties of matter, ranging from atoms and molecular systems to condensed matter and exotic states. These spectroscopies are widely used in investigations at the major DOE synchrotron x-ray facilities and related centers world-wide. In addition to fundamental theory, a major goal of our project has been the development of computational software that implements the theory for calculations of x-ray spectra of various materials throughout the periodic table. Due to the complex nature of x-ray spectra, quantitative theory is essential for its interpretation. The theory is challenging since it involves excited state electronic structure and many-body correlation effects that go beyond independent particle approximations like DFT or Hartree-Fock. Moreover, the experimental investigations typically involve a broad range of energy, time, and temperature scales, from the UV-Vis to hard x-ray energies of order 10 4 eV, and temperatures T from ambient up to the warm-dense-matter regime where the Fermi energy kBTF is of order a few eV, i.e., temperatures of order 105 Kelvin. This broad range of experimental conditions has fostered many novel theoretical approaches and computational techniques, many of which we have developed systematically over the duration of the grant. In contrast to the traditional wave-function approach of quantum theory and electronic structure methods, our theoretical approach is based on modern Green's function techniques. This approach is better suited for aperiodic structures, excited states, and broad spectral ranges, since it avoids the computational bottlenecks of sum-over-states approaches, as in the Fermi golden rule. This theoretical framework has been incorporated into efficient, user-friendly x-ray spectroscopy software which is now used routinely worldwide to simulate and analyze spectra. These theoretical tools provide an essential complement to synchrotron and next-generation light sources, which are used to investigate complex materials with ever increasing precision. Moreover, the synergism between theory, computation and experiment contributed by our research enhances scientific understanding and creates opportunities for innovations in materials and energy science and in many fields. As documented in this Report, this research grant has been remarkably successful in achieving these goals. In particular, this grant has supported the development of the x-ray spectroscopy software suite known as FEFF (named for an effective scattering amplitude f eff in the theory). The FEFF codes have become one of the premier tools for quantitative simulations of x-ray spectra as documented by many thousands of citations in the Web of Science and Google-Scholar.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The effects of earth model uncertainty on the inversion of seismic data for seismic source functions

SUMMARY We use Monte Carlo simulations to explore the effects of earth model uncertainty on the estimation of the seismic source time functions that correspond to the six independent components of the point source seismic moment tensor. Specifically, we invert synthetic data using Green’s functions estimated from a suite of earth models that contain stochastic density and seismic wave-speed heterogeneities. We find that the primary effect of earth model uncertainty on the data is that the amplitude of the first-arriving seismic energy is reduced, and that this amplitude reduction is proportional to the magnitude of the stochastic heterogeneities. Also, we find that the amplitude of the estimated seismic source functions can be under- or overestimated, depending on the stochastic earth model used to create the data. This effect is totally unpredictable, meaning that uncertainty in the earth model can lead to unpredictable biases in the amplitude of the estimated seismic source functions.

58 GEOSCIENCES↗
DOE OSTI2023

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗
DOE OSTI2021

Examination of the sensitivity of quasifree reactions to details of the bound-state overlap functions

It is often stated that heavy-ion nucleon knockout reactions are mostly sensitive to the tails of the bound-state wave functions. In contrast, (p,2p) and (p,pn) reactions are known to access information on the full overlap functions within the nucleus. We analyze the oxygen isotopic chain and explore the differences between single-particle wave functions generated with potential models, used in the experimental analysis of knockout reactions, and ab initio computations from self-consistent Green's function theory. Contrary to common belief, we find that not only the tail of the overlap functions, but also their internal part is assessed in both reaction mechanisms, which are crucial to yield accurately determined spectroscopic information.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Equation of motion coupled-cluster cumulant approach for intrinsic losses in x-ray spectra

We present an equation of motion coupled cluster approach for calculating and understanding intrinsic inelastic losses in core level x-ray absorption spectra (XAS). The method is based on a factorization of the transition amplitude in the time-domain, which leads to a convolution of an effective one-body spectrum and the core-hole spectral function. The spectral function characterizes these losses in terms of shake-up excitations and satellites, and is calculated using a cumulant representation of the core-hole Green’s function that includes non-linear corrections. The one-body spectrum also includes orthogonality corrections that enhance the XAS at the edge.

Rehr, John J.↗
DOE OSTI2021

Electronic structure of chromium trihalides beyond density functional theory

In this work, we explore the electronic band structure of freestanding monolayers of chromium trihalides Cr X 3 , X = Cl, Br, I, within an advanced ab initio theoretical approach based on the use of Green's function functionals. We compare the local density approximation with the quasiparticle self-consistent GW (QS GW ) approximation and its self-consistent extension ( QS G W ^ ) by solving the particle-hole ladder Bethe-Salpeter equations to improve the effective interaction W . We show that, at all levels of theory, the valence band consistently changes shape in the sequence Cl → Br → I , and the valence band maximum shifts from the M point to the Γ point. By analyzing the dynamic and momentum-dependent self-energy, we show that QS G W ^ adds to the localization of the systems in comparison with QS GW , thereby leading to a narrower band and reduced amount of halogens in the valence band manifold. Further analysis shows that X = Cl is most strongly correlated, and X = I is least correlated (most bandlike) as the hybridization between Cr d and X p enhances in the direction Cl → Br → I . For CrBr 3 and CrI 3 , we observe remarkable differences between the QS GW and QS G W ^ valence band structures, while their eigenfunctions are very similar. We show that weak perturbations, like moderate strain, weak changes to the d - p hybridization, and adding small U , can flip the valence band structures between these two solutions in these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

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↗
DOE OSTI2023

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↗
DOE OSTI2022

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↗
DOE OSTI2023

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↗
DOE OSTI2021

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↗
DOE OSTI2021

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↗
DOE OSTI2023

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↗
DOE OSTI2025

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↗
DOE OSTI2026

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