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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Quantum solver for single-impurity Anderson models with particle-hole symmetry

Quantum embedding methods, such as dynamical mean-field theory (DMFT), provide a powerful framework for investigating strongly correlated materials. A central computational bottleneck in DMFT is in solving the Anderson impurity model (AIM), whose exact solution is classically intractable for large bath sizes. In this work, we benchmark a quantum-classical hybrid solver tailored for particle-hole symmetric AIMs, using the variational quantum eigensolver to prepare the ground state of the model with shallow quantum circuits. The solver uses shallow quantum ansätze and one set of variational parameters to prepare the ground state and its particle and hole excitations, enabling the construction of the impurity Green’s function through a continued-fraction expansion. We evaluate the performance of this approach across a few bath sizes and interaction strengths under noisy, shot-limited conditions. We compare three optimization routines (COBYLA, Adam, and L-BFGS-B) in terms of convergence and fidelity, assess the benefits of estimating a quantum-computed moment correction to the variational energies, and benchmark the approach by comparing the density of states computed from the impurity Green’s function against that obtained using a classical pipeline. Our results demonstrate the feasibility of Green’s function construction on near-term devices and establish practical benchmarks for quantum impurity solvers embedded within self-consistent DMFT loops.

Karabin, Mariia [ORNL]↗

DOE OSTI2025

$\gamma$-ray angular correlations, distributions and linear polarization in tracking arrays

With $\gamma $-ray tracking arrays, the angular correlation and, especially, the linear polarization of $\gamma $ rays can be measured with high precision. However, because of the complexity of these spectrometers, measuring angular distributions is more involved. The techniques and associated formalisms were developed and established in the 1950s. The experimental setups of that era were relatively simple when compared to the advanced $\gamma $-ray tracking arrays available today. This paper aims to review and refine these techniques, unify the nomenclatures and, specifically, adapt them for use with the newest tracking arrays. In addition, new methods for producing the required response functions for these arrays are introduced.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Extracting spinon self-energies from two-dimensional coherent spectroscopy

Two-dimensional coherent spectroscopy (2DCS) is a nonlinear spectroscopy technique capable of identifying whether apparent continua in linear response are made out of multiplets of sharp deconfined quasiparticles. This makes it a potent tool for experimental identification of fractionalized phases. Previous discussions have focused on limits where the quasiparticles in question are infinitely long lived. In this paper we discuss 2DCS in the regime where the fractionalized quasiparticles can themselves decay. We introduce a powerful path-integral-based approach, whereby the computation of nonlinear susceptibilities reduces to an efficient exercise in diagrammatic perturbation theory. Here, we apply this method to compute the 2DCS response of the one-dimensional transverse field Ising model, in the presence of integrability-breaking perturbations. We discuss aspects of the self-energy of the fractionalized quasiparticles that may be extracted via 2DCS, such as the momentum-dependent decay rate.

1-dimensional spin chains↗

NASA NTRS1989

Computational Aspects of Sensitivity Calculations in Linear Transient Structural Analysis

A study has been performed focusing on the calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, transient response problems. One significant goal was to develop and evaluate sensitivity calculation techniques suitable for large-order finite element analyses. Accordingly, approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite element model. Much of the research focused on the accuracy of both response quantities and sensitivities as a function of number of vectors used. Two types of sensitivity calculation techniques were developed and evaluated. The first type of technique is an overall finite difference method where the analysis is repeated for perturbed designs. The second type of technique is termed semianalytical because it involves direct, analytical differentiation of the equations of motion with finite difference approximation of the coefficient matrices. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models.

Greene, William H.↗

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

Developing Advanced Genetic and Synthetic Biology Tools for Improved Algae Productivity

Under the UCSD-led Productivity Enhanced Algae and ToolKits (PEAK) project we sought to develop tools to enable the production of valuable co-products, as well as strategies to decrease the cost of cultivation, with an intended outcome of enabling the economic production of algae-based biofuel. We have developed new genetic tools and high throughput selection methods for eukaryotic green algae and used them to express a high value recombinant protein; a growth factor known as Osteopontin (OPN). Recombinant OPN protein expression levels were improved through a newly developed rapid breeding and mutagenesis strain improvement method to generate a robust strain capable of growing in outdoors raceway conditions in brackish water that tolerated high pH and large temperature swings. We further applied the new genetic tools to express recombinant protein in a recently bioprospected strain with even more robust growth characteristics. In response to the COVID-19 pandemic we used our improved workflow to rapidly express, purify, and characterize a recombinant SARS-CoV-2 spike protein Receptor Binding Domain (RBD) and demonstrate that it functionally interacts with its cognate human host receptor ACE2. The newly identified extremophile strain Chlamydomonas sp (402), both mating type (mt+) and mating type (mt-), have been deposited in the Chlamydomonas Resource Center (https://www.chlamycollection.org/), making them available world-wide with no restrictions. In addition, all of genetic tools are also deposited at this site and are also available world-wide with no restrictions.

09 BIOMASS FUELS↗

NASA NTRS2007

Simple and Efficient Numerical Evaluation of Near-Hypersingular Integrals

Recently, significant progress has been made in the handling of singular and nearly-singular potential integrals that commonly arise in the Boundary Element Method (BEM). To facilitate object-oriented programming and handling of higher order basis functions, cancellation techniques are favored over techniques involving singularity subtraction. However, gradients of the Newton-type potentials, which produce hypersingular kernels, are also frequently required in BEM formulations. As is the case with the potentials, treatment of the near-hypersingular integrals has proven more challenging than treating the limiting case in which the observation point approaches the surface. Historically, numerical evaluation of these near-hypersingularities has often involved a two-step procedure: a singularity subtraction to reduce the order of the singularity, followed by a boundary contour integral evaluation of the extracted part. Since this evaluation necessarily links basis function, Green s function, and the integration domain (element shape), the approach ill fits object-oriented programming concepts. Thus, there is a need for cancellation-type techniques for efficient numerical evaluation of the gradient of the potential. Progress in the development of efficient cancellation-type procedures for the gradient potentials was recently presented. To the extent possible, a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. However, since the gradient kernel involves singularities of different orders, we also require that the transformation leaves remaining terms that are analytic. The terms "normal" and "tangential" are used herein with reference to the source element. Also, since computational formulations often involve the numerical evaluation of both potentials and their gradients, it is highly desirable that a single integration procedure efficiently handles both.

Fink, Patrick W.↗

NASA NTRS1993

Double layer propagation in experiments with electron beam injection

Electron beam injection into a plasma is investigated using the analytical inverted Bernstein-Green-Kruskal method. Particle number and momentum conservation laws are applied to evaluate the propagation velocity and potential drop on the leading edge of the beam. Electric potential is supposed to be monotonic, thus the leading front has a double-layer-like structure. For the case of cold particles, analytical expressions for the double layer velocity and potential drop are obtained. It is pointed out that double layer velocity differs from the initial electron speed: even for weak beams a noticeable deceleration takes place. Strong beams are found incapable of penetrating into plasma - their propagation velocity is very small. Ambient electrons undergo a considerable acceleration forming a return current which neutralizes the injector. Possible instability of the distribution functions is discussed.

Bruskin, L. G.↗

NASA NTRS1995

Electromagnetic Scattering From a Rectangular Cavity Recessed in a 3-D Conducting Surface

The problem of electromagnetic (EM) scattering from an aperture backed by a rectangular cavity recessed in a three-dimensional conducting body is analyzed using the coupled field integral equation approach. Using the free-space Green's function, EM fields scattered outside the cavity are determined in terms of (1) an equivalent electric surface current density flowing on the three-dimensional conducting surface of the object including the cavity aperture and (2) an equivalent magnetic surface current density flowing over the aperture only. The EM fields inside the cavity are determined using the waveguide modal expansion functions. Making the total tangential electric and magnetic fields across the aperture continuous and subjecting the total tangential electric field on the outer conducting three-dimensional surface of the object to zero, a set of coupled integral equations is obtained. The equivalent electric and magnetic surface currents are then obtained by solving the coupled integral equation using the Method of Moments. The numerical results on scattering from rectangular cavities embedded in various three-dimensional objects are compared with the results obtained by other numerical techniques.

Deshpande, M. D.↗

NASA NTRS2017

Cyclotron Resonant Scattering Feature Simulations II. Description of the CRSF Simulation Process

Context. Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons o_ quantized plasma electrons in the strong magnetic field (of the order 1012 G) close to the surface of an accreting X-ray pulsar. Due to the complex scattering cross-sections, the line profiles of CRSFs cannot be described by an analytic expression. Numerical methods, such as Monte Carlo (MC) simulations of the scattering processes, are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems.Aims. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time.Methods. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time-consuming simulations possible on convenient time scales. Furthermore, it can generate responses to monoenergetic photon injections, producing Green's functions, which can be used later to generate spectra for arbitrary continua.Results. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green's function tables are available online (see link in footnote, page 1).

Schwarm, F.-W.↗

DOE OSTI2021

High-Resolution Spatio-Temporal Estimation of Net Ecosystem Exchange in Ice-Wedge Polygon Tundra Using In Situ Sensors and Remote Sensing Data

Land-atmosphere carbon exchange is known to be extremely heterogeneous in arctic ice-wedge polygonal tundra regions. In this study, a Kalman filter-based method was developed to estimate the spatio-temporal dynamics of daytime average net ecosystem exchange (NEEday) at 0.5-m resolution over a 550 m by 700 m study site. We integrated multi-scale, multi-type datasets, including normalized difference vegetation indices (NDVIs) obtained from a novel automated mobile sensor system (or tram system) and a greenness index map obtained from airborne imagery. We took advantage of the significant correlations between NDVI and NEEday identified based on flux chamber measurements. The weighted average of the estimated NEEday within the flux-tower footprint agreed with the flux tower data in term of its seasonal dynamics. We then evaluated the spatial variability of the growing season average NEEday, as a function of polygon geomorphic classes; i.e., the combination of polygon types—which are known to present different degradation stages associated with permafrost thaw—and microtopographic features (i.e., troughs, centers and rims). Our study suggests the importance of considering microtopographic features and their spatial coverage in computing spatially aggregated carbon exchange.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Three-dimensional analysis of microbunched electron cooling

Microbunched electron cooling (MBEC, Ratner, ), a special category of coherent electron cooling, has been under consideration as a candidate method for the strong hadron cooling component of the Electron Ion Collider (EIC) at BNL. Most of the theoretical treatments of this scheme so far have been based on simplified models for the key space charge effect in which electrons and/or hadrons are approximated by charged disks. Although such disk-based models provide a flexible tool that captures a good part of the underlying physics, they lack the fidelity and robustness of a rigorous analysis that explicitly treats the hadron and cooler electron beams as collections of point charges. In this paper, we present such a rigorous, three-dimensional (3D) model, both from a theory-based and a simulation perspective. Specifically, we calculate the generalized wakefield (or Green’s function) of the cooling system—a crucial quantity that determines the overall performance—paying special attention to the dependence of the wake on the hadron transverse positions. Additionally, our treatment takes into account previously neglected 3D effects that can cause damping of the plasma oscillations in the amplification section of the cooler, such as focusing and the angular spread of the electron beam. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Electrode-molecule energy level offsets in a gold-benzene diamine-gold single molecule tunnel junction

One means for describing electron transport across single molecule tunnel junctions (MTJs) is to use density functional theory (DFT) in conjunction with a nonequilibrium Green’s function formalism. This description relies on interpreting solutions to the Kohn–Sham (KS) equations used to solve the DFT problem as quasiparticle (QP) states. Many practical DFT implementations suffer from electron self-interaction errors and an inability to treat charge image potentials for molecules near metal surfaces. For MTJs, the overall effect of these errors is typically manifested as an overestimation of electronic currents. Correcting KS energies for self-interaction and image potential errors results in MTJ current–voltage characteristics in close agreement with measured currents. An alternative transport approach foregoes a QP picture and solves for a many-electron wavefunction on the MTJ subject to open system boundary conditions. It is demonstrated that this many-electron method provides similar results to the corrected QP picture for electronic current. The analysis of these two distinct approaches is related through corrections to a junction’s electronic structure beyond the KS energies for the case of a benzene diamine molecule bonded between two gold electrodes.

Szepieniec, Mark S.↗

DOE OSTI2025

Archaeal protein containing domain of unknown function 2193 undergoes oligomeric reconfiguration upon iron–sulfur cluster binding

Methanogenic archaea are particularly rich in iron–sulfur proteins, yet their roles remain largely enigmatic. Here, we characterized aMethanococcus voltae(Mvo) protein from the domain of unknown function (DUF) 2193 family, a group of proteins present primarily in archaea and characterized by a conserved cysteine‐rich C‐terminal motif.MvoDUF2193 was heterologously expressed and characterized by a range of spectroscopic and analytical methods. The results demonstrate thatMvoDUF2193 binds a single [4Fe–4S] cluster per subunit and that cluster occupancy regulates the transition from an apo tetramer to a [4Fe–4S] monomeric form. We hypothesize thatMvoDUF2193 serves a regulatory role in the cell, mediated by [Fe–S] cluster binding and changes in oligomeric state.

Biochemistry & Molecular Biology↗

DOE OSTI2024

Millimeter-Scale Exfoliation of hBN with Tunable Flake Thickness for Scalable Encapsulation

As a two-dimensional (2D) dielectric material, hexagonal boron nitride (hBN) is in high demand for applications in photonics, nonlinear optics, and nanoelectronics. Unfortunately, the high-throughput preparation of macroscopic-scale, high-quality hBN flakes with controlled thickness is an ongoing challenge, limiting device fabrication and technological integration. In this report we present a metal thin-film exfoliation method to prepare hBN flakes with millimeter-scale dimension, near-unity yields, and tunable flake thickness distribution from 1-7 layers, a substantial improvement over scotch tape exfoliation. The single crystallinity and high quality of the exfoliated hBN are demonstrated with optical microscopy, atomic force microscopy, Raman spectroscopy, and second harmonic generation. We further explore a possible mechanism for the effectiveness and selectivity based on thin-film residual stress measurements, density functional theory calculations, and transmission electron microscopy imaging of the deposited metal films. We find that the magnitude of the residual tensile stress induced by thin film deposition plays a key role in determining exfoliated flake thickness in a manner which closely resembles 3D semiconductor spalling. Lastly, we demonstrate that our exfoliated, large-area hBN flakes can be readily incorporated as encapsulating layers for other 2D monolayers. Altogether, this method brings us one step closer to the high throughput, mass production of hBN-based 2D photonic, optoelectronic, and quantum devices.

2D materials↗

DOE OSTI2026

Engineering Escherichia coli for Urease-Driven Synthesis of Metal Oxide Nanomaterials

The development of functional nanomaterials with controlled morphologies is essential for advancements in medicine, electronics and computing, energy, catalysis, and environmental applications. However, conventional synthesis methods often demand high energy input and pose significant environmental challenges. Urease-based biomineralization presents an efficient, eco-friendly alternative for nanomaterial production under mild conditions. In this study, we engineered Escherichia coli ( E. coli ) to express a urease gene cluster from Sporosarcina pasteurii using CRAGE-Duet technology. The engineered strain successfully synthesized calcium carbonate and calcium phosphate crystals. Expanding the approach, we synthesized metal oxide nanoparticles, including hematite (Fe 2 O 3 ), and nanocrystalline anatase titanium dioxide (TiO 2 ). These nanomaterials were characterized by electron microscopy, demonstrating the potential of E. coli as a sustainable and versatile platform for green nanomaterial synthesis.

bacteria↗

DOE OSTI2023

Defects go green: using defects in nanomaterials for renewable energy and environmental sustainability

Induction of point defects in nanomaterials can bestow upon them entirely new physics or augment their pre-existing physical properties, thereby expanding their potential use in green energy technology. Predicting structure-property relationships for defects a priori is challenging, and developing methods for precise control of defect type, density, or structural distribution during synthesis is an even more formidable task. Hence, tuning the defect structure to tailor nanomaterials for enhanced device performance remains an underutilized tool in materials design. We review here the state of nanomaterial design through the lens of computational prediction of defect properties for green energy technology, and synthesis methods to control defect formation for optimal performance. We illustrate the efficacy of defect-focused approaches for refining nanomaterial physics by describing several specific applications where these techniques hold potential. Most notably, we focus on quantum dots for reabsorption-free solar windows and net-zero emission buildings, oxide cathodes for high energy density lithium-ion batteries and electric vehicles, and transition metal dichalcogenides for electrocatalytic green hydrogen production and carbon-free fuels.

14 SOLAR ENERGY↗

DOE OSTI2021

Chemometric optimisation of pressurised liquid extraction for the determination of alliin and S‐allyl‐cysteine in giant garlic ( Allium ampeloprasum L.) by liquid chromatography tandem mass spectrometry

Abstract Introduction Giant garlic is a functional food that contains different kinds of bioactive molecules with beneficial effects on chronic noncommunicable diseases like diabetes and cardiovascular conditions. Considering biosynthesis pathways, abundance, and biological activity, alliin and S‐allyl‐cysteine were used as chemical markers of organosulphur compounds present in giant garlic. Objective To establish a chemometric optimisation of pressurised liquid extraction for the determination of alliin and S‐allyl‐cysteine in giant garlic by liquid chromatography tandem mass spectrometry (LC–MS/MS). Methodology Samples were blanched (ca. 90°C for 10 min) to inactivate alliinase and γ‐glutamyl transpeptidase enzymes and then freeze‐dried. Chemometric optimisation was performed via response surface methodology based on central composite design (CCD). Organosulphur compound yields were determined applying a validated LC–MS/MS method in multiple reaction monitoring (MRM) mode using the following transitions: for alliin m/z 178 → 74 and for S‐allyl‐cysteine m/z 162 → 41. Results According to CCD results, under constant conditions of pressure (1500 psi) and time (20 min), the optimal conditions for pressurised liquid extraction of alliin and S‐allyl‐cysteine were 70.75 and 68.97% v/v of ethanol in extraction solvent and 76.45 and 98.88°C as extraction temperature, respectively. Multiple response optimisation for the simultaneous extraction of both organosulphur compounds was established via desirability function. Under these conditions, 2.70 ± 0.27 mg g −1 dry weight (DW) of alliin and 2.79 ± 0.22 mg g −1 DW of S‐allyl‐cysteine were extracted. Conclusions These results clearly demonstrated that pressurised liquid extraction is an efficient green technique to extract bioactive organosulphur compounds from giant garlic. Extraction yields were significantly ( p < 0.05) higher than those obtained with conventional ultra‐turrax extraction.

Peterssen‐Fonseca, Darlene↗