Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Green's function methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 829 records · Page 46

Radiation and scattering from printed antennas on cylindrically conformal platforms

The goal was to develop suitable methods and software for the analysis of antennas on cylindrical coated and uncoated platforms. Specifically, the finite element boundary integral and finite element ABC methods were employed successfully and associated software were developed for the analysis and design of wraparound and discrete cavity-backed arrays situated on cylindrical platforms. This work led to the successful implementation of analysis software for such antennas. Developments which played a role in this respect are the efficient implementation of the 3D Green's function for a metallic cylinder, the incorporation of the fast Fourier transform in computing the matrix-vector products executed in the solver of the finite element-boundary integral system, and the development of a new absorbing boundary condition for terminating the finite element mesh on cylindrical surfaces.

Kempel, Leo C.↗

Computational Nuclear Physics and Post Hartree-Fock Methods

We present a computational approach to infinite nuclear matter employing Hartree-Fock theory, many-body perturbation theory and coupled cluster theory. These lectures are closely linked with those of chapters 9, 10 and 11 and serve as input for the correlation functions employed in Monte Carlo calculations in chapter 9, the in-medium similarity renormalization group theory of dense fermionic systems of chapter 10 and the Green's function approach in chapter 11. We provide extensive code examples and benchmark calculations, allowing thereby an eventual reader to start writing her/his own codes. We start with an object-oriented serial code and end with discussions on strategies for porting the code to present and planned high-performance computing facilities.

Lietz, Justin↗

Derivation of instanton rate theory from first principles

Instanton rate theory is used to study tunneling events in a wide range of systems including low-temperature chemical reactions. Despite many successful applications, the method has never been obtained from first principles, relying instead on the “Im F” premise. In this paper, the same expression for the rate of barrier penetration at finite temperature is rederived from quantum scattering theory [W. H. Miller, S. D. Schwartz, and J. W. Tromp, J. Chem. Phys. 79, 4889 (1983)] using a semiclassical Green’s function formalism. This justifies the instanton approach and provides a route to deriving the rate of other processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring the exact limits of the real-time equation-of-motion coupled cluster cumulant Green’s functions

In this paper, we analyze the properties of the recently proposed real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green’s function approach [Rehr et al., J. Chem. Phys. 152, 174113 (2020)]. We specifically focus on identifying the limitations of the original time-dependent coupled cluster (TDCC) ansatz and propose an enhanced double TDCC ansatz, ensuring the exactness in the expansion limit. In addition, we introduce a practical cluster-analysis-based approach for characterizing the peaks in the computed spectral function from the RT-EOM-CC cumulant Green’s function approach, which is particularly useful for the assignments of satellite peaks when many-body effects dominate the spectra. Our preliminary numerical tests focus on reproducing, approximating, and characterizing the exact impurity Green’s function of the three-site and four-site single impurity Anderson models using the RT-EOM-CC cumulant Green’s function approach. The numerical tests allow us to have a direct comparison between the RT-EOM-CC cumulant Green’s function approach and other Green’s function approaches in the numerical exact limit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stochastic bias from non-Gaussian initial conditions

In this article, we show that a stochastic form of scale-dependent halo bias arises in multi-source inflationary models, where multiple fields determine the initial curvature perturbation. We derive this effect for general non-Gaussian initial conditions and study various examples, such as curvaton models and quasi-single field inflation. Furthermore, we present a general formula for both the stochastic and the non-stochastic parts of the halo bias, in terms of the N-point cumulants of the curvature perturbation at the end of inflation. At lowest order, the stochasticity arises if the collapsed limit of the four-point function is boosted relative to the square of the three-point function in the squeezed limit. We derive all our results in two ways, using the barrier crossing formalism and the peak-background split method. In a companion paper, we prove that these two approaches are mathematically equivalent.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On the integration of image sources in exact image method of field analysis

Convergence conditions of image integration in the exact image method of field calculation have been investigated, and the method is extended to include more general media than previously considered. It is demonstrated that the integral is well behaved and the method works best when the field is calculated in a medium with less loss. If the medium has more loss, it is shown that the image line may enter the improper half-space and still produce valid results. Means of correctly selecting the integration path branch for the case of crossing branch cuts of the Green function in complex integration planes are proposed.

Lindell, I. V.↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth Observing System) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the "Independent Pixel Approximation" or IPA. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

Multi-Beam Approach for Accelerating Alignment and Calibration of HyspIRI-Like Imaging Spectrometers

A paper describes an optical stimulus that produces more consistent results, and can be automated for unattended, routine generation of data analysis products needed by the integration and testing team assembling a high-fidelity imaging spectrometer system. One key attribute of the system is an arrangement of pick-off mirrors that provides multiple input beams (five in this implementation) to simultaneously provide stimulus light to several field angles along the field of view of the sensor under test, allowing one data set to contain all the information that previously required five data sets to be separately collected. This stimulus can also be fed by quickly reconfigured sources that ultimately provide three data set types that would previously be collected separately using three different setups: Spectral Response Function (SRF), Cross-track Response Function (CRF), and Along-track Response Function (ARF), respectively. This method also lends itself to expansion of the number of field points if less interpolation across the field of view is desirable. An absolute minimum of three is required at the beginning stages of imaging spectrometer alignment.

Eastwood, Michael L.↗

Quantifying EGS Reservoir Complexity with an Integrated Geophysical Approach-Improved Resolution Ambient Seismic Noise Interferometry

The primary goal of the project was met, as improved resolution (from kilometers to tens/hundreds of meters and to a depth of at least 3 km) was added to the seismic interferometry-estimated parameters relative to the work reported in AltaRock Energy Inc. (2014) for the Dixie Valley Geothermal Wellfield. A procedure for initial fault identification, using seismic interferometry, was developed. The project team further improved, and applied a new, cost-effective, and non-invasive exploration method, developed and validated since 2010, during two previous experiments, at Dixie Valley and Soda Lake, Nevada. The method, initially designed for hydrothermal exploration, was adapted to EGS exploration conditions, which require deeper crustal information. Regions of elevated energy for both the 2-7 Hz and 5-15 Hz frequency bands are located above intersections of multiple faults. Thus, an unprecedented, innovative seismic field deployment and ambient seismic noise and signal processing were successfully used to estimate crustal structure and to identify fault locations at depths up to 3 km, at a cost per kilometer of survey that in some cases may be as much as ten times less than the cost of a controlled-source survey. Based on valid and useful seismic data collected and analyzed during Phase I, in Phase II the team has further quantified statistically the advantage of using seismic analysis, in combination with other geological and geophysical techniques, for EGS target identification. A power spectral density analysis of the passive-seismic empirical Green’s functions reveals significantly higher power recorded over the exploration area at the power plant, and in the drilling area, near the BB station 10201, especially on the horizontal components. In this study, the most important relationship observed during the prior report, namely the one between velocity and temperature was explored. The areas that show anomalously high predicted temperature in the geostatistics occur where multiple faults intersect. Not surprisingly, such areas span between the Dixie Valley range front fault and the basin-ward Piedmont faults, in the vicinity of numerous wells. These locations also fall in the area of high EGS favorability from the prior 2014 study. The anomalies also coincide with the high power spectral density values. The result of altering all favorability criteria for the new data reached a maximum of 1/8th of a point and were not large enough to change any favorability value reported in the 2014 study. The new data reinforce the results obtained from the previous 2014 study. The current high-resolution ambient seismic noise project provided excellent results in defining faults in the wellfield. As such the identification of known faults by the project demonstrated the utility of this method for identifying faults at an exploration stage in geothermal development. Such fault identification, while extremely useful in exploration investigation, has not added sufficient new data for the generation of significantly different DVGW favorability maps in the current project. We believe the cost savings of our passive-seismic technique for EGS, over common industry active-source survey techniques, may have a minimum of 30%, and a maximum of 90%. There are many significant non-cost benefits of passive seismic surveys. Our passive-source method provides an environmentally-friendly seismic survey with minimal to negligible impacts. The economic benefit is greater EGS project feasibility.

15 GEOTHERMAL ENERGY↗

3D-Printed Porous Supramolecular Sorbents for Cobalt Recycling

Electronic waste recycling is a recognized global challenge that requires new strategies to bind and release critical materials selectively, such as cobalt present in lithium-ion batteries. To address this challenge, hierarchical 3D-printed porous polymer scaffolds bearing supramolecular receptors were prepared using vat photopolymerization and their cobalt binding profiles were examined as a function of matrix polarity. By combining high-resolution digital light processing (DLP) with polymerization-induced phase separation (PIPS), functional acrylic copolymer networks with micrometer-level precision of geometry and nanometer-level pores were generated. Covalent integration of a methacrylate-functionalized bisdicyclohexyl acetamide (BDCA-MA) receptor enabled binding and release of cobalt(II) chloride (CoCl 2 ) via a solvent polarity switch mechanism involving a change in solvent from ethanol to water. The present structures proved reusable as shown by sustained high binding efficiency over five bind and release cycles. Finally, this platform represents a “green” and energy conscious method for future electronic waste recycling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deep Learning Prediction of Protein Complex Structures

Proteins interact to form protein complex to carry out biological functions such as catalytic chemical reaction. Therefore, it is important to develop computational methods to predict protein-protein interaction and the structures of protein complexes to study and enhance protein function. In this project, we successfully developed several deep learning methods to predict inter-protein contacts and the reinforcement learning and optimization methods to reconstruct protein complex structures from predicted inter-chain contacts. The methods were integrated with the MULTICOM protein complex structure prediction system and applied to predict the complex structures of biomass production-related proteins of green algae. During the two and a half years of research and development, all the specific milestones of the project were achieved successfully. 16 publications/manuscripts were produced. 10 software tools were developed. A patent application was submitted. Our MULTICOM predictors leveraging some tools developed in this project were ranked among the top predictors in the 15th Critical Assessment of Techniques for Protein Structure Prediction (CASP15) in 2022.

59 BASIC BIOLOGICAL SCIENCES↗

Strategies for designing biocatalysts with new functions

Enzymes can be optimized to accelerate chemical transformations via a range of methods. In this review, we showcase how protein engineering and computational design techniques can be interfaced to develop highly efficient and selective biocatalysts.

09 BIOMASS FUELS↗

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

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM (atmospheric radiation measurement) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth observing system) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of 'radiative smoothing' of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. 'Offbeam' Lidar returns are now being used to estimate the cloud 'radiative Green's function' (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the 'Independent Pixel Approximation' or IPA. This and other measurements of 3D (three dimensional) transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international 'Intercomparison of 3D Radiation Codes' (I3RC) program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

Transformer-based operator learning framework for self-energy in strongly correlated systems

We introduce Σ-Attention, a transformer-based operator-learning framework for approximating the self-energy operator of strongly correlated electronic systems. By creating a batched dataset that combines results from three complementary approaches, i.e., many-body perturbation theory, strong-coupling expansion, and exact diagonalization, each effective in specific parameter regimes, Σ-Attention is applied to learn an accurate approximation for the self-energy operator that is valid across a wide range of parameter regimes. This hybrid strategy leverages the strengths of existing methods while relying on the transformer's ability to generalize beyond individual limitations. More importantly, the scalability of the transformer architecture allows the learned self-energy to be extended to systems with larger sizes, leading to much improved computational scaling. Using the one-dimensional Hubbard model, we demonstrate that Σ-Attention can accurately predict the Matsubara Green's function of large systems with a wide range of coupling strength. Our framework offers a promising and scalable pathway for studying strongly correlated systems with many possible generalizations.

Zhu, Yuanran↗

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↗

The small-voxel tracking algorithm for simulating chemical reactions among diffusing molecules

Simulating the evolution of a chemically reacting system using the bimolecular propensity function, as is done by the stochastic simulation algorithm and its reaction-diffusion extension, entails making statistically inspired guesses as to where the reactant molecules are at any given time. Those guesses will be physically justified if the system is dilute and well-mixed in the reactant molecules. Otherwise, an accurate simulation will require the extra effort and expense of keeping track of the positions of the reactant molecules as the system evolves. One molecule-tracking algorithm that pays careful attention to the physics of molecular diffusion is the enhanced Green's function reaction dynamics (eGFRD) of Takahashi, Tănase-Nicola, and ten Wolde [Proc. Natl. Acad. Sci. U.S.A. 107, 2473 (2010)]. We introduce here a molecule-tracking algorithm that has the same theoretical underpinnings and strategic aims as eGFRD, but a different implementation procedure. Called the small-voxel tracking algorithm (SVTA), it combines the well known voxel-hopping method for simulating molecular diffusion with a novel procedure for rectifying the unphysical predictions of the diffusion equation on the small spatiotemporal scale of molecular collisions. Indications are that the SVTA might be more computationally efficient than eGFRD for the problematic class of non-dilute systems. A widely applicable, user-friendly software implementation of the SVTA has yet to be developed, but we exhibit some simple examples which show that the algorithm is computationally feasible and gives plausible results.

97 MATHEMATICS AND COMPUTING↗