Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized functions”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The Generalized Green’s function Cluster Expansion: A Python package for simulating polarons

We present an efficient implementation of the Generalized Green’s function Cluster Expansion (GGCE), which is a new method for computing the ground-state properties and dynamics of polarons (single electrons coupled to lattice vibrations) in model electron-phonon systems. The GGCE works at arbitrary temperature and is well suited for a variety of electron-phonon couplings, including, but not limited to, site and bond Holstein and Peierls (Su-SchriefferHeeger) couplings, and couplings to multiple phonon modes with different energy scales and coupling strengths. Quick calculations can be performed efficiently on a laptop using solvers from NumPy and SciPy, or in parallel at scale using the PETSc sparse linear solver engine.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spectroscopy and crystallography define carotenoid oxygenases as a new subclass of mononuclear non-heme Fe II enzymes

Carotenoid cleavage dioxygenases (CCDs) are non-heme Fe II enzymes that catalyze the oxidative cleavage of alkene bonds in carotenoids, stilbenoids, and related compounds. How these enzymes control the reaction of dioxygen (O 2 ) with their alkene substrates is unclear. Here, we apply spectroscopy in conjunction with X-ray crystallography to define the iron coordination geometry of a model CCD, CAO1 (Neurospora crassa carotenoid oxygenase 1), in its resting state and following substrate binding and coordination sphere substitutions. Resting CAO1 exhibits a five-coordinate (5C), square pyramidal Fe II center that undergoes steric distortion toward a trigonal bipyramidal geometry in the presence of piceatannol. Titrations with the O 2 -analog, nitric oxide, show a >100-fold increase in iron–nitric oxide affinity upon substrate binding, defining a crucial role for the substrate in activating the Fe II site for O 2 reactivity. The importance of the 5C Fe II structure for reactivity was probed through mutagenesis of the second-sphere Thr151 residue of CAO1, which occludes ligand binding at the sixth coordination position. A T151G substitution resulted in the conversion of the iron center to a six-coordinate state and a 135-fold reduction in apparent catalytic efficiency toward piceatannol compared with the wildtype enzyme. Substrate complexation resulted in partial six-coordinate to 5C conversion, indicating solvent dissociation from the iron center. Additional substitutions at this site demonstrated a general functional importance of the occluding residue within the CCD superfamily. Taken together, these data suggest an ordered mechanism of CCD catalysis occurring via substrate-promoted solvent replacement by O 2 . CCDs thus represent a new class of mononuclear non-heme Fe II enzymes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Reaction losses of charged particles in CsI(Tl) crystals

To efficiently detect energetic light charged particles, it is common to use arrays of energy-loss telescopes involving two or more layers of detection media. As the energy of the particles increases, thicker layers are usually needed. However, carrying out measurements with thick-telescopes may require corrections for the losses due to nuclear reactions induced by the incident particles on nuclei within the detector and for the scattering of incident particles out of the detector, without depositing their full energy in the active material. In this paper, we develop a method for measuring such corrections and determine the reaction and out-scattering losses for data measured with the silicon-CsI(Tl) telescopes of the newly developed HiRA10 array. Here, the extracted efficiencies are in good agreement with model predictions using the GEANT4 reaction loss algorithm for Z = 1 and Z = 2 isotopes. After correcting for the HiRA10 geometry, we obtain a general function that describes the loss of efficiency due to reaction losses in CsI(Tl) crystals as a function of range.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Coupling High-Throughput Experiments and Regression Algorithms to Optimize PGM-Free ORR Electrocatalyst Synthesis

Over the past decades, significant improvement has been achieved in the performance of platinum group metal-free (PGM-free) materials as an alternative to Pt-based electrocatalysts for oxygen reduction reaction (ORR). However, further progress in ORR activity requires evaluation of precursors and synthesis approaches. In response to this challenge, we generated a first of its kind experimental data set of 36 samples using high-throughput synthesis and activity measurements. Several control parameters (e.g., Fe precursor identity, the precursor content, and pyrolysis temperature) were varied. We then developed several state-of-the-art machine learning (ML) based regression models to predict ORR activity, dependent on selected synthesis variables. Through an iterative algorithm, higher prediction accuracy (smaller root-mean-square error) was achieved. We identified that gradient boosting regression (GBR) and support vector regression (SVR), among several methods, work best for this data set. Aided by our ML-based surrogate models, we decided to alter catalyst synthesis conditions, which resulted in a 36% increase in measured ORR activity in comparison to the maximum ORR mass activity value of 21.9 A/g catalyst in the original data set. Overall, this combined experiment and machine learning approach represents a promising path forward toward developing highly efficient next-generation ORR electrocatalysts and, more generally, functional materials.

25 ENERGY STORAGE↗
DOE OSTI2023

Marine energy converters: Potential acoustic effects on fishes and aquatic invertebrates

The potential effects of underwater anthropogenic sound and substrate vibration from offshore renewable energy development on the behavior, fitness, and health of aquatic animals is a continuing concern with increased deployments and installation of these devices. Initial focus of related studies concerned offshore wind. However, over the past decade, marine energy devices, such as a tidal turbines and wave energy converters, have begun to emerge as additional, scalable renewable energy sources. Because marine energy converters (MECs) are not as well-known as other anthropogenic sources of potential disturbance, their general function and what is known about the sounds and substrate vibrations that they produce are introduced. Furthermore, while most previous studies focused on MECs and marine mammals, this paper considers the potential of MECs to cause acoustic disturbances affecting nearshore and tidal fishes and invertebrates. In particular, the focus is on particle motion and substrate vibration from MECs because these effects are the most likely to be detected by these animals. Finally, an analysis of major data gaps in understanding the acoustics of MECs and their potential impacts on fishes and aquatic invertebrates and recommendations for research needed over the next several years to improve understanding of these potential impacts are provided.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

BCLink User Documentation

BCLink is a shared programming library providing general functionality to read temporally- and spatially-varying boundary condition data (contained in Exodus files) into MDG codes (ParaDyn and Diablo). This document is intended for prospective users of BCLink, presenting a summary of the available features of the library, including details regarding the input syntax with accompanying examples. The library’s capabilities are demonstrated through example problems run in ParaDyn, and using the BCRemap command line utility – a stand-alone preprocessing tool which leverages the native functionality provided by the BCLink library to remap surface data between dissimilar meshes.

42 ENGINEERING↗
DOE OSTI2026

A tale of two extremes: Temperature sensitivity of carbon loss from cool and hot soils

Soils represent the largest terrestrial carbon (C) pool, and the flux of carbon dioxide (CO 2 ) from soils to the atmosphere is ~ 6-10 times more than anthropogenic emissions. Understanding responses of soil CO 2 emissions to warming is crucial for evaluating feedback to ongoing environmental changes. The relationship between microbial respiration and temperature is typically modeled using a Q 10 function. Generally, observations of the apparent Q 10 of soil respiration are higher for cold vs. warm ecosystems, reflecting expected biophysical controls of Arrhenius kinetics. However, results from two field warming experiments in the tropics contradict this expectation, both observing extraordinarily high soil respiration responses to in situ warming. Our overall objective for the proposed work is to reduce uncertainty in temperature sensitivity of soil C loss by systematically synthesizing underlying mechanisms related to soil C turnover and stabilization. We are evaluating the temperature sensitivity of soil respiration in ecosystems across temperature extremes (e.g., arctic/boreal and tropical systems) by integrating data collected from field warming experiments with machine learning and biogeochemical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Exploring the Nucleon Structure via Deep Electroproduction Processes

Understanding the internal structure of the nucleon is a fundamental goal of modern physics, which aims at a comprehensive framework describing the internal dynamics of quarks and gluons. Among other structure functions, Generalized Parton Distributions (GPDs) offer a powerful framework for describing the nucleon dynamics by correlating the longitudinal momentum and the transverse position of its internal partons. Such a correlation provides a three-dimensional picture of the nucleon and enables access to fundamental properties, including the internal pressure distributions and the parton’s angular momentum contribution to the nucleon's total spin, thereby playing a central role in resolving the nucleon spin puzzle. At the Thomas Jefferson National Facility (JLab), polarized electron beam experiments allow for probing GPDs through the measurement of hard exclusive processes. Among the cleanest experimental channels, we find the electro-production of a real photon through the Deeply Virtual Compton Scattering (DVCS) mechanism. The first data-taking period of the CLAS12 program, taking place in 2018, allowed for unique DVCS Beam Spin Asymmetry (BSA) measurements in the phase space covered by a 10.6 GeV polarized electron beam impinging on an unpolarized liquid hydrogen target. Although detecting all final-state particles ensures exclusivity of the process, conservation laws indicate that it is not mandatory. I adopt an approach omitting the direct detection of the recoil proton, providing a simplified yet effective event selection strategy that boosts statistics and gives access to a larger phase space sensitive to the underlying GPD dynamics through BSA and cross section measurements. The Double DVCS (DDVCS) process promises a dedicated mapping of GPDs. Contributing to the electro-production of a lepton pair cross-section, the DDVCS reaction extends DVCS by allowing the final-state photon to be virtual, enriching the kinematic phase space and providing unique access to the internal correlations encoded by GPDs. A feasibility study is conducted to assess the potential of future DDVCS measurements at Jefferson Lab and the future Electron-Ion Collider (EIC). While Jefferson Lab will provide DDVCS measurements in the valence region through the SoLID$\mu$ and $\mu$CLAS12 experimental projects, in the long term, the EIC will provide complementary measurements in the sea region, both accessing unprecedented information about GPDs in a phase space region otherwise inaccessible. Taken together, these investigations demonstrate both the current capabilities and future opportunities for probing GPDs through exclusive processes. The experimental analysis of DVCS at CLAS12 provides precise measurements within an established framework, while the phenomenological study of DDVCS opens the door to richer and more comprehensive explorations with future detectors and facilities.

Alvarado, Juan [Université Paris-Saclay: Gif-sur-Y↗
DOE OSTI2024

A generalized Selberg zeta function for flat space cosmologies

Flat space cosmologies (FSCs) are time dependent solutions of three-dimensional (3D) gravity with a vanishing cosmological constant. They can be constructed from a discrete quotient of empty 3D flat spacetime and are also called shifted-boost orbifolds. Using this quotient structure, we build a new and generalized Selberg zeta function for FSCs, and show that it is directly related to the scalar 1-loop partition function. We then propose an extension of this formalism applicable to more general quotient manifolds $\mathcal{M}$/ℤ, based on representation theory of fields propagating on this background. Our prescription constitutes a novel and expedient method for calculating regularized 1-loop determinants, without resorting to the heat kernel. We compute quasinormal modes in the FSC using the zeroes of a Selberg zeta function, and match them to known results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Identification of crystal plasticity model parameters by multi-objective optimization integrating microstructural evolution and mechanical data

Crystal plasticity models evolve a polycrystalline yield surface using meso-scale descriptions of deformation mechanisms. The activation of deformation mechanisms is governed by crystallography and a set of model parameters, which are typically calibrated through the fitting of mechanical data such as stress–strain curves and elastic lattice strains. Microstructural data such as phase fractions and texture evolution are used for verifying crystal plasticity parameters. In this study, we use a multi-objective genetic algorithm to identify hardening parameters from flow stress curves with an option to incorporate texture into the optimization approach. Robust, generalized objective functions are developed and used to identify sets of parameters pertaining to dislocation density-based hardening laws in visco-plastic and elasto-plastic self-consistent (VPSC and EPSC) homogenization models. First, the parameters are identified for pure Nb directly from texture using an objective function based on generalized spherical harmonics. Since texture evolution is driven by the relative contribution of active slip systems, the parameters governing the evolution of slip resistance ratios can be recovered from fitting discrete textures at a series of strains. Next, a comprehensive set of load reversal data for dual phase (DP) 780 steel is used to fit a hardening law and a back-stress law in EPSC. Finally, parameters pertaining to a complex hardening law for the evolution of slip and twinning in pure α-Ti are identified. Remarkably, using texture as an objective in combination with stress–strain objectives constrains the model of Ti to fully reproduce not only stress–strain and texture evolution but also hierarchical twinning measurements as a function of initial grain size and texture. Furthermore, given an appropriate model fit to representative experimental texture evolution, underlying twin volume fractions contributing to texture evolution can be predicted.

42 ENGINEERING↗
DOE OSTI2022

Correlators of four light-ray operators in CCFT

Building on results in arXiv:2203.04255, we compute the correlator of four gluon light-ray operators in celestial CFT. We find that it is described by Fox H-functions and generalized I-functions of multiple variables. We also analyze light-ray correlators for a scalar amplitude involving the exchange of a massive scalar.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Beyond leaf habit: generalities in plant function across 97 tropical dry forest tree species

Leaf habit has been hypothesized to define a linkage between the slow-fast plant economic spectrum and the drought resistance-avoidance trade-off in tropical forests (‘slow-safe vs fast-risky’). However, variation in hydraulic traits as a function of leaf habit has rarely been explored for a large number of species. We sampled leaf and branch functional traits of 97 tropical dry forest tree species from four sites to investigate whether patterns of trait variation varied consistently in relation to leaf habit along the ‘slow-safe vs fast-risky’ trade-off. Leaf habit explained from 0% to 43.69% of individual trait variation. We found that evergreen and semi-deciduous species differed in their location along the multivariate trait ordination when compared to deciduous species. While deciduous species showed consistent trait values, evergreen species trait values varied as a function of the site. Last, trait values varied in relation to the proportion of deciduous species in the plant community. Here We found that leaf habit describes the strategies that define drought avoidance and plant economics in tropical trees. However, leaf habit alone does not explain patterns of trait variation, which suggests quantifying site-specific or species-specific uncertainty in trait variation as the way forward.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Some aspects of the five-dimensional Lovelock black hole spacetime: Strong homotopy retract, perihelion precession and quasistationary levels

Highlights: • Some topological aspects are explored. • The perihelion precession is calculated. • The general Heun functions are applied to obtain analytical solutions. • The quasibound states and the radial wave eigenfunctions are obtained. In this work we explore some mathematical physics aspects of the spherically symmetric Lovelock black hole in high dimensions. Intended for this aim, we thoroughly consider the metric corresponding to the five-dimensional Lovelock black hole spacetime. We construct the strong retractions by the geodesic equations on the background under consideration. As a result, from the topological point of view, we construct the theory of strong homotopy retract, which will allow us, in principle, to better understand some of its suitable applications on astrophysics and cosmology, in particular, in the analysis of the spacetime singularities. We find the solutions of the equation of motion for both radial and angular coordinates, and then we describe the outer (“exterior”) and lower (“interior”) apparent horizons. Indeed, the outer apparent horizon is the last surface from which the light waves could still escape from the black hole. Thus, it is meaningful to analyze some physical phenomena related to quantum particles propagating outside the exterior apparent horizon, in particular, we discuss the quasistationary levels of scalar fields and their radial wave functions, which are given in terms of the general Heun functions. We also calculate the perihelion precession in this background.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Molecular dipole moment learning via rotationally equivariant derivative kernels in molecular-orbital-based machine learning

This study extends the accurate and transferable molecular-orbital-based machine learning (MOB-ML) approach to modeling the contribution of electron correlation to dipole moments at the cost of Hartree–Fock computations. A MOB pairwise decomposition of the correlation part of the dipole moment is applied, and these pair dipole moments could be further regressed as a universal function of MOs. The dipole MOB features consist of the energy MOB features and their responses to electric fields. An interpretable and rotationally equivariant derivative kernel for Gaussian process regression (GPR) is introduced to learn the dipole moment more efficiently. The proposed problem setup, feature design, and ML algorithm are shown to provide highly accurate models for both dipole moments and energies on water and 14 small molecules. To demonstrate the ability of MOB-ML to function as generalized density-matrix functionals for molecular dipole moments and energies of organic molecules, we further apply the proposed MOB-ML approach to train and test the molecules from the QM9 dataset. The application of local scalable GPR with Gaussian mixture model unsupervised clustering GPR scales up MOB-ML to a large-data regime while retaining the prediction accuracy. In addition, compared with the literature results, MOB-ML provides the best test mean absolute errors of 4.21 mD and 0.045 kcal/mol for dipole moment and energy models, respectively, when training on 110 000 QM9 molecules. The excellent transferability of the resulting QM9 models is also illustrated by the accurate predictions for four different series of peptides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗
DOE OSTI2023

Resource frugal optimizer for quantum machine learning

Quantum-enhanced data science, also known as quantum machine learning (QML), is of growing interest as an application of near-term quantum computers. Variational QML algorithms have the potential to solve practical problems on real hardware, particularly when involving quantum data. However, training these algorithms can be challenging and calls for tailored optimization procedures. Specifically, QML applications can require a large shot-count overhead due to the large datasets involved. In this work, we advocate for simultaneous random sampling over both the dataset as well as the measurement operators that define the loss function. We consider a highly general loss function that encompasses many QML applications, and we show how to construct an unbiased estimator of its gradient. This allows us to propose a shot-frugal gradient descent optimizer called Refoqus (REsource Frugal Optimizer for QUantum Stochastic gradient descent). Our numerics indicate that Refoqus can save several orders of magnitude in shot cost, even relative to optimizers that sample over measurement operators alone.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Resource selection functions based on hierarchical generalized additive models provide new insights into individual animal variation and species distributions

Habitat selection studies are designed to generate predictions of species distributions or inference regarding general habitat associations and individual variation in habitat use. Such studies frequently involve either individually indexed locations gathered across limited spatial extents and analyzed using resource selection functions (RSFs) or spatially extensive locational data without individual resolution typically analyzed using species distribution models. Both analytical methodologies have certain desirable features, but analyses that combine individual- and population-level inference with flexible non-linear functions may provide improved predictions while accounting for individual variation. Here, we describe how RSFs can be fit using hierarchical generalized additive models (HGAMs) using widely available software, providing a means to explore individual variation in habitat associations and to generate species distribution maps. We used GPS tracking data from golden eagles Aquila chrysaetos from across eastern North America with four environmental predictors to generate monthly distribution models. We considered three model structures that assumed different amounts of individual variation in the functional relationship between predictors and habitat use and used k-fold cross-validation to compare model performance. Models accounting for individual variability in shape and smoothness of functional responses performed best. Eagles exhibited the least amount of individual variation in response to land cover variables during winter months, with most individuals more closely adhering to the population-level trend. During the summer months, eagles exhibited more substantial individual variation in shape and smoothness of the functional relationships, suggesting some need to account for individual variation in eagle habitat use for both inferential and predictive purposes, during this time of year. Because they allow users to blend flexible functions with random effects structures and are well-supported by a variety of software platforms, we believe that HGAMs provide a useful addition to the suite of analyses used for modeling habitat associations or predicting species distributions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Exact reaction coordinates for flap opening in HIV-1 protease

The primary goal of protein science is to understand how proteins function, which requires understanding the functional dynamics responsible for transitions between different functional structures of a protein. A central concept is the exact reaction coordinates that can determine the value of committor for any protein configuration, which provide the optimal description of functional dynamics. Despite intensive efforts, identifying the exact reaction coordinates (RCs) in complex molecules remains a formidable challenge. Using the recently developed generalized work functional, we report the discovery of the exact RCs for an important functional process—the flap opening of HIV-1 protease. Our results show that this process has six RCs, each one is a linear combination of ~240 backbone dihedrals, providing the precise definition of collectivity and cooperativity in the functional dynamics of a protein. Applying bias potentials along each RC can accelerate flap opening by 10 3 to 10 4 folds. The success in identifying the RCs of a protein with 198 residues represents a significant progress beyond that of the alanine dipeptide, currently the only other complex molecule for which the exact RCs for its conformational changes are known. Our results suggest that the generalized work functional (GWF) might be the fundamental operator of mechanics that controls protein dynamics.

Science & Technology - Other Topics↗