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

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

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

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

Denoising of imaginary time response functions with Hankel projections

Imaginary-time response functions of finite-temperature quantum systems are often obtained with methods that exhibit stochastic or systematic errors. Reducing these errors comes at a large computational cost—in quantum Monte Carlo simulations, the reduction of noise by a factor of two incurs a simulation cost of a factor of four. In this paper, we relate certain imaginary-time response functions to an inner product on the space of linear operators on Fock space. We then show that data with noise typically does not respect the positive definiteness of its associated Gramian. The Gramian has the structure of a Hankel matrix. As a method for denoising noisy data, we introduce an alternating projection algorithm that finds the closest positive definite Hankel matrix consistent with noisy data. We test our methodology at the example of fermion Green's functions for continuous-time quantum Monte Carlo data and show remarkable improvements of the error, reducing noise by a factor of up to 20 in practical examples. We argue that Hankel projections should be used whenever finite-temperature imaginary-time data of response functions with errors is analyzed, be it in the context of quantum Monte Carlo, quantum computing, or in approximate semianalytic methodologies. Published by the American Physical Society 2024

Yu, Yang (ORCID:0000000186178878)↗
DOE OSTI2024

Neural Networks to Find the Optimal Forcing for Offsetting the Anthropogenic Climate Change Effects

Abstract Of great relevance to climate engineering is the systematic relationship between the radiative forcing to the climate system and the response of the system, a relationship often represented by the linear response function (LRF) of the system. However, estimating the LRF often becomes an ill-posed inverse problem due to high-dimensionality and nonunique relationships between the forcing and response. Recent advances in machine learning make it possible to address the ill-posed inverse problem through regularization and sparse system fitting. Here, we develop a convolutional neural network (CNN) for regularized inversion. The CNN is trained using the surface temperature responses from a set of Green’s function perturbation experiments as imagery input data together with data sample densification. The resulting CNN model can infer the forcing pattern responsible for the temperature response from out-of-sample forcing scenarios. This promising proof of concept suggests a possible strategy for estimating the optimal forcing to negate certain undesirable effects of climate change. The limited success of this effort underscores the challenges of solving an inverse problem for a climate system with inherent nonlinearity. Significance Statement Predicting the climate response for a given climate forcing is a direct problem, while inferring the forcing for a given desired climate response is often an inverse, ill-posed, problem, posing a new challenge to the climate community. This study makes the first attempt to infer the radiative forcing for a given target pattern of global surface temperature response using a deep learning approach. The resulting deeply trained convolutional neural network inversion model shows promise in capturing the forcing pattern corresponding to a given surface temperature response, with a significant implication on the design of an optimal solar radiation management strategy for curbing global warming. This study also highlights the technical challenges that future research should prioritize in seeking feasible solutions to the inverse climate problem.

Ren, Huiying↗
DOE OSTI2022

Transverse Enhancement, Longitudinal Quenching and Coulomb Sum Rule in e-$^{12}$C and e-$^{16}$O Scattering

We present a short summary of a phenomenological analysis of all available electron scattering data on 12 C 12 C (about 6600 differential cross section measurements) and on 16 O 16 O (about 250 measurements) within the framework of the quasielastic superscaling model (including Pauli blocking). All quasielastic and inelastic cross section measurements are included down to the lowest momentum transfer ? q (including photo-production data). We find that there is enhancement of the transverse quasielastic response function and quenching of the longitudinal response function at low ? q (in addition to Pauli blocking). We extract parameterizations of a ? ? ? ? ? ? ? ? ? ? ? ? ? ? multiplicative low ? q "Longitudinal Quenching Factor" and an ? ? ? ? ? ? ? ? additive "Transverse Enhancement" contribution. Additionally, we find that the excitation of nuclear states contribute significantly (up to 30\%) to the Coulomb Sum Rule. We extract the most accurate determination of the sum rule to date and find it to be in disagreement with random phase approximation based calculations but in reasonable agreement with recent theoretical calculations such as "First Principle Green's Function Monte Carlo".

Bodek, Arie↗
DOE OSTI2024

Robust Measurements of n -Point Correlation Functions of Driven-Dissipative Quantum Systems on a Digital Quantum Computer

We propose and demonstrate a unified hierarchical method to measure n -point correlation functions that can be applied to driven, dissipative, or otherwise open or nonequilibrium quantum systems. In this method, the time evolution of the system is repeatedly interrupted by interacting an ancilla qubit with the system through a controlled operation, and measuring the ancilla immediately afterward. We discuss the robustness of this method as compared to other ancilla-based interferometric techniques (such as the Hadamard test), and highlight its advantages for near-term quantum simulations of open quantum systems. We implement the method on a quantum computer in order to measure single-particle Green’s functions of a driven-dissipative fermionic system. This Letter shows that dynamical correlation functions for driven-dissipative systems can be robustly measured with near-term quantum computers. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Broken Luttinger theorem in the two-dimensional Fermi-Hubbard model

One of the fundamental questions about high-temperature cuprate superconductors is the size of the Fermi surface underlying the superconducting state. By analyzing the single-particle spectral function for the Fermi-Hubbard model as a function of repulsion U and chemical potential μ, we find that the Fermi surface in the normal state undergoes a transition from a large Fermi surface matching the Luttinger volume as expected in a Fermi liquid, to a Fermi surface that encloses fewer electrons that we dub the “Luttinger breaking” phase, as the Mott insulator is approached. This transition into a non-Fermi-liquid phase that violates the Luttinger count occurs at a critical density in the absence of any other broken symmetry. We obtain the Fermi-surface contour from the spectral weight A k (ω = 0) and from an analysis of the singularities of the Green's function Re G k (E = 0), calculated using determinantal quantum Monte Carlo and analytic continuation methods. We discuss our numerical results in connection with experiments on Hall measurements, scanning tunneling spectroscopy, and angle-resolved photoemission spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Single-Molecule Electron Transport in Peptoids

Peptoids are structural analogs of peptides in which side chains are appended to the backbone nitrogen rather than the α-carbon. The sequence-defined modularity of peptoids enables precise control over structure−function relationships, enabling applications in energy storage and biomedical materials. Despite recent progress, the role of sequence and conformation on electron transport in peptoid molecules is not fully understood. Here, we synthesize a library of peptoid oligomers and characterize their molecular electronic properties using the scanning tunneling microscope-break junction (STM-BJ) technique. Our results show well-defined electron transport behavior for peptoid sequences containing aromatic side groups lacking hydrogen bonds (H-bonds) and without chemical substitutions at the N−C α position. This behavior fundamentally differs from electron transport in peptides, where H-bond interactions give rise to higher conductance states. All-atom molecular dynamics (MD) simulations are used to understand the conformational heterogeneity of peptoids, and molecular conformations obtained from MD simulations are used in quantum mechanical calculations based on the nonequilibrium Green’s function−density functional theory (NEGF-DFT) formalism. In all cases, computational results are in reasonable qualitative agreement with experiments. Our work demonstrates that the conductance behavior of peptoids depends on monomer identity, including side-chain aromaticity and substitution at the N−C α position. Overall, this work provides new insights into the structure−function relationships governing electron transport in peptoid-based materials and establishes design rules for peptoid-based molecular junctions.

Charge transport↗
DOE OSTI2022

Ab initio estimation of E 2 strengths in Li 8 and its neighbors by normalization to the measured quadrupole moment

For electric quadrupole (E2) observables, which depend on the large-distance tails of the nuclear wave function, ab initio no-core configuration interaction (NCCI) calculations converge slowly, making meaningful predictions challenging to obtain. Nonetheless, the calculated values for different E2 matrix elements, particularly those involving levels with closely-related structure (e.g., within the same rotational band) are found to be robustly proportional. This observation suggests that a known value for one observable may be used to determine the overall scale of E2 strengths, and thereby provide predictions for others. In particular, we demonstrate that meaningful predictions for E2 transitions may be obtained by calibration to the ground-state quadrupole moment. Here, we test this approach for well-measured low-lying E2 transitions in 7 Li and 9 Be, then provide predictions for transitions in 8 Li and 9 Li. In particular, we address the 2 + → 1 + transition in 8 Li, for which the reported measured strength exceeds ab initio Green's function Monte Carlo (GFMC) predictions by over an order of magnitude.

6 ≤ A ≤ 19↗
DOE OSTI2024

Journey to Time-Variable Moment Tensors through Inversion of Acoustic and Seismoacoustic Data

We explore the capability of acoustic and seismoacoustic datasets to directly resolve a complex, time-variable source consisting of a buried mechanism, represented as a moment tensor, and a spall mechanism, represented as a vertical force at the surface. Traditionally, each component of a resolved moment tensor assumes one underlying source time function, which likely fails to capture the full evolution of a dynamic source, such as an explosion followed by slip on near-source joints or development of spallation. Specifically, we expand previous work to resolve a time-variable moment tensor using single-modality and joint-modality inversion frameworks through analysis of infrasound and seismoacoustic data recorded as part of the Source Physics Experiment Phase II: Dry Alluvium Geology (DAG). We investigate the impact of including signals from seismic-to-air coupling that are local to each infrasound sensor in comparison to mainly atmosphere-propagating acoustic signals, which occur from coupling of the wavefield from the subsurface to the atmosphere directly above the source. Additionally, we assess the ability of our inversion algorithm to fit observed infrasound data using a variety of time-variable source mechanisms. First, we consider the buried moment tensor source alone, which assumes that the determined Green’s functions incorporate effects from spallation or that the impact from spallation is minimal. Second, we examine the estimated buried moment tensor and vertical surface spallation as terms that must both be resolved in the inversion. Third, we assess the ability for an estimated vertical surface spallation source to fit the acoustic data on its own. Finally, we compare results from the joint inversion of both seismic geophone and infrasound acoustic data for the buried-only source compared to buried and spallation sources. Our results are a preliminary investigation into the applications of the inversion technique to recorded datasets and show the technique has limited capabilities using acoustic data alone. Instead, this method shows promise for seismic and seismoacoustic datasets to resolve the time-variable mechanisms of a buried source.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Metastable antiphase boundary ordering in CaFe 2 O 4

CaFe 2 O 4 is an $\textit{S}$ = 5/2 antiferromagnet exhibiting two magnetic orders that shows regions of coexistence at some temperatures. Using a Green’s function formalism, we model neutron scattering data of the spin wave excitations in this material, elucidating the microscopic spin Hamiltonian. In doing so, we suggest that the low-temperature A phase order (↑↑↓↓) finds its origins in the freezing of antiphase boundaries created by thermal fluctuations in a parent B phase order (↑↓↑↓). In this work, the low-temperature magnetic order observed in CaFe 2 O 4 is thus the result of a competition between the exchange coupling along $\textit{c}$, which favors the B phase, and the single-ion anisotropy, which stabilizes thermally generated antiphase boundaries, leading to static metastable A phase order at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Breakdown of the drift-diffusion model for transverse spin transport in a disordered Pt film

Spin-accumulation and spin-current profiles are calculated for a disordered Pt film subjected to an in-plane electric current within the nonequilibrium Green's function approach. In the bulklike region of the sample, this approach captures the intrinsic spin Hall effect found in other calculations. Near the surfaces, the results reveal qualitative differences with the results of the widely used spin-diffusion model, even when the boundary conditions are modified to try to account for them. One difference is that the effective spin-diffusion length for transverse spin transport is significantly different from its longitudinal counterpart and is instead similar to the mean-free path. Furthermore, this feature may be generic for spin currents generated via the intrinsic spin Hall mechanism because of the differences in transport mechanisms compared to longitudinal spin transport. Orbital accumulation in the Pt film is only significant in the immediate vicinity of the surfaces and has a small component penetrating into the bulk only in the presence of spin-orbit coupling, as a secondary effect induced by the spin accumulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Machine Learning Approach to Modeling of Neutral Particles Transport in Plasma

A propagator‐based approach is investigated for Monte‐Carlo (MC) modeling of neutral particle transport in fusion boundary plasmas. The propagator is based on a Green's function for the neutral kinetic equation, which depends on the plasma profiles. A neural network (NN)‐based model for the propagator provides a fast and accurate solution for the neutral distribution function in plasma. Preliminary results from a small 1D test problem look encouraging. The proposed approach, a propagator‐based NN model for neutral transport in plasma, has potential for generalization to higher dimensions and efficient coupling with plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Renormalization group improved photon impact factors and the high energy virtual photon scattering

We perform the renormalization group improved collinear resummation of the photon-gluon impact factors. We construct the resummed cross section for virtual photon-photon (γ ∗ γ ∗ ) scattering which incorporates the impact factors and BFKL gluon Green’s function up to the next-to-leading logarithmic accuracy in energy. The impact factors include important kinematical effects which are responsible for the most singular poles in Mellin space at next-to-leading order. Further conditions on the resummed cross section are obtained by requiring the consistency with the collinear limits. Our analysis is consistent with previous impact factor calculations at NLO, apart from a new term proportional to C F that we find for the longitudinal polarization. Finally, we use the resummed cross section to compare with the LEP data on the γ ∗ γ ∗ cross section and with previous calculations. The resummed result is lower than the leading logarithmic approximation but higher than the pure next-to-leading one, and is consistent with the experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Boundary correlators and the Schwarzian mode

The effective low temperature dynamics of near-extremal black holes is governed by the quantum fluctuations of the Schwarzian mode of JT gravity. Utilizing as a proxy a planar charged black hole in asymptotically Anti-de-Sitter spacetime, we investigate the effects of these fluctuations on a probe scalar field. The corresponding holographic real-time boundary correlators are computed following a holographic renormalization procedure, using the dubbed gravitational Schwinger-Keldysh geometry (grSK) and known exact results of boundary correlators from the near-horizon region. This analysis gives rise to a retarded Green’s function that decays as a power law for late Lorentzian times. Its analytic structure indicates the presence of a branch cut in the complex frequency domain at finite temperature. These features are a non-perturbative hallmark that prevails as long as the planar transverse space is kept compact.

2D Gravity↗
DOE OSTI2021

Constraints on quasinormal modes and bounds for critical points from pole-skipping

We consider a holographic thermal state and perturb it by a scalar operator whose associated real-time Green’s function has only gapped poles. These gapped poles correspond to the non-hydrodynamic quasinormal modes of a massive scalar perturbation around a Schwarzschild black brane. Relations between pole-skipping points, critical points and quasinormal modes in general emerge when the mass of the scalar and hence the dual operator dimension is varied. First, this novel analysis reveals a relation between the location of a mode in the infinite tower of quasinormal modes and the number of pole-skipping points constraining its dispersion relation at imaginary momenta. Second, for the first time, we consider the radii of convergence of the derivative expansions about the gapped quasinormal modes. These convergence radii turn out to be bounded from above by the set of all pole-skipping points. Furthermore, a transition between two distinct classes of critical points occurs at a particular value for the conformal dimension, implying close relations between critical points and pole-skipping points in one of those two classes. We show numerically that all of our results are also true for gapped modes of vector and tensor operators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗