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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Role of water model on ion dissociation at ambient conditions

In this work, we study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion–ion interactions.

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

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗
DOE OSTI2024

Depth and microtopography influence microbial biogeochemical processes in a forested peatland

Background and aims: Peat-accumulating wetlands have undulating surfaces of raised areas (hummocks) and depressions (hollows). Hummock-hollow microtopography in relation to the water table influences the distribution of plant species, root density, and microbial community composition, which could in turn alter carbon (C) and nitrogen (N) cycling within peatlands. We used paired hummock and hollow cores from a boreal, forested peatland to assess how microtopography influences peatland microbial function and, in turn, ecosystem C and N cycling. Methods: The peat was analyzed for microbial biomass and potential enzyme activity in 10 cm depth increments relative to the water table, resulting in two increments for hollows and three for hummocks, which has a raised increment above the water table. Results: Across hummocks and hollows, microbial C and N and fungal biomass generally decreased with depth from the peat surface. In contrast, potential enzyme activity often increased with depth, but this varied within enzyme functional groups according to topography, depth, or both. The potential enzyme activity of C-N degrading peptidases, for example, differed across the five topography × depth increments with the lowest rate in the aerated hummocks. Hummocks compose approximately 66% of the land area at our study site and would therefore underestimate C turnover by an average of 25% if solely used to extrapolate patterns across a forested bog. Conclusion: In conclusion, our results suggest that asynchrony in C and N cycling across the undulating surface of forested peatlands impacts our ability to accurately predict biogeochemical cycling across this important ecosystem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Enhanced intersystem crossing of boron dipyrromethene by TEMPO radical

Radical enhanced intersystem crossing (EISC) of organic chromophores is an important approach to generate a long-lived triplet state for various electronic and optoelectronic applications. However, structural factors and design rules to promote EISC are not entirely clear. In this work, we report a series of boron dipyrromethene (BODIPY) derivatives covalently linked with a 2,2,6,6-tetramethyl-1-piperidinyloxyl (TEMPO) radical with varying distances and topologies. We show that the incorporation of the TEMPO radical to BODIPY results in strong fluorescence quenching by up to 85% as a result of EISC and enhanced internal conversion. In BDP-2AR [2-(4-methyleneamino-TEMPO) BODIPY], a dyad with the shortest BODIPY–TEMPO through-bond distance, we observe the fastest EISC rate (τ isc = 1.4 ns) and the longest triplet excited state lifetime (τ T = 32 μs) compared to other distance and geometry variations. Contrary to previous reports and a general presumption, the BODIPY–TEMPO through-bond distance in this system does not play a significant role on the triplet formation rate and yield. Finally, density functional theory suggests a folding of the TEMPO radical to form a sandwich-like structure with a BODIPY ring that leads to a decrease in the through-space distance, providing a new and an interesting insight for the radical enhanced intersystem.

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

Coaction and double-copy properties of configuration-space integrals at genus zero

We investigate configuration-space integrals over punctured Riemann spheres from the viewpoint of the motivic Galois coaction and double-copy structures generalizing the Kawai-Lewellen-Tye (KLT) relations in string theory. For this purpose, explicit bases of twisted cycles and cocycles are worked out whose orthonormality simplifies the coaction. We present methods to efficiently perform and organize the expansions of configuration-space integrals in the inverse string tension α' or the dimensional-regularization parameter ϵ of Feynman integrals. Generating-function techniques open up a new perspective on the coaction of multiple polylogarithms in any number of variables and analytic continuations in the unintegrated punctures. We present a compact recursion for a generalized KLT kernel and discuss its origin from intersection numbers of Stasheff polytopes and its implications for correlation functions of two-dimensional conformal field theories. We find a non-trivial example of correlation functions in (p, 2) minimal models, which can be normalized to become uniformly transcendental in the p → ∞ limit.

Conformal Field Theory↗
DOE OSTI2021

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

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

Nearest neighbour distributions: New statistical measures for cosmological clustering

ABSTRACT The use of summary statistics beyond the two-point correlation function to analyse the non-Gaussian clustering on small scales, and thereby, increasing the sensitivity to the underlying cosmological parameters, is an active field of research in cosmology. In this paper, we explore a set of new summary statistics – the k-Nearest Neighbour Cumulative Distribution Functions (kNN-CDF). This is the empirical cumulative distribution function of distances from a set of volume-filling, Poisson distributed random points to the k-nearest data points, and is sensitive to all connected N-point correlations in the data. The kNN-CDF can be used to measure counts in cell, void probability distributions, and higher N-point correlation functions, all using the same formalism exploiting fast searches with spatial tree data structures. We demonstrate how it can be computed efficiently from various data sets – both discrete points, and the generalization for continuous fields. We use data from a large suite of N-body simulations to explore the sensitivity of this new statistic to various cosmological parameters, compared to the two-point correlation function, while using the same range of scales. We demonstrate that the use of kNN-CDF improves the constraints on the cosmological parameters by more than a factor of 2 when applied to the clustering of dark matter in the range of scales between 10 and $40\, h^{-1}\, {\rm Mpc}$. We also show that relative improvement is even greater when applied on the same scales to the clustering of haloes in the simulations at a fixed number density, both in real space, as well as in redshift space. Since the kNN-CDF are sensitive to all higher order connected correlation functions in the data, the gains over traditional two-point analyses are expected to grow as progressively smaller scales are included in the analysis of cosmological data, provided the higher order correlation functions are sensitive to cosmology on the scales of interest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Composition, Growth, Succession, and Function in the Cladophora Microbiome: Insights From Quantitative Stable Isotope Probing and NanoSIMS Imaging

The branching green macroalga Cladophora glomerata and its epiphytic microbiome dominate summer biomass in the Eel River, a Northern California river under Mediterranean (summer drought, winter rain) seasonality. Green Cladophora streamers proliferate in early summer, then change to yellow and then red-brown as epiphyte loads increase. Here, we characterised successional changes in epiphytic bacteria on Cladophora, examining both community composition and growth rates, using quantitative Stable Isotope Probing (qSIP) and 16S rRNA gene amplicon sequencing. The number of bacterial taxa increased with succession while growth rates peaked in the middle stage. NanoSIMS imaging confirmed high sulphur (S) concentrations in Cladophora cell walls relative to surrounding biomass, coinciding with a bloom of sulphur bacteria (bacteria that reduce or oxidise sulphur/sulphates). In general, relative abundances and growth rates were independent, indicating that either metric alone is insufficient for understanding how taxonomy and functional groups affect ecosystem processes. For instance, the relative abundance of nitrogen fixers peaked in the late summer when their relative growth rates were slowest. Such patterns may be driven by space competition limiting growth. Together, changes in abundance and relative growth rates suggest different limiting factors for different functional groups in the Cladophora microbiome at multiple successional stages.

Biological and medical sciences↗
DOE OSTI2023

Impact of solvation on the GW quasiparticle spectra of molecules

First-principles calculations for electrochemistry require accurate treatment of both electronic structure and solvation. The perturbative GW approximation starting from density functional theory (DFT) calculations accurately models materials systems with varying dimensionality. Continuum solvation models enable efficient treatment of solvation effects in DFT calculations, but their applications with beyond-DFT electronic structure methods such as GW have been limited. Here, we introduce the frequency-dependent liquid polarizability from a nonlocal continuum solvation model in the screened Coulomb interaction of full-frequency GW calculations with a solvated DFT starting point. In this study, we show that the liquid screening contributions substantially reduce the HOMO–LUMO gap of molecules by 3–5 eV, while solvent effects on the DFT starting point negligibly impact the GW gap. The resulting framework facilitates the simultaneous electronic and solvation accuracy needed for first-principles electrochemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Nickel/gallium modified HZSM-5 for ethane aromatization: Influence of metal function on reactivity and stability

Bimetallic catalysts often outperform their monometallic counterparts due to the synergistic effect. Here we reported that the Ni/Ga co-functionalized HZSM-5 catalyst, specifically with Ni 1 /Ga 1 stoichiometric ratio, shows significantly improved selectivity (above 80 %) and stability than the Ni/HZSM-5 catalyst and higher activity (with identical selectivity) than the Ga/HZSM-5 catalyst in ethane aromatization. According to the transient experiment, the Ga and Ni 1 Ga 1 /HZSM-5 catalysts show an induction period of a few minutes. However, for Ni/HZSM-5 and Ni 1 Ga 1/3 /HZSM-5, the induction period increased to >10 min and >30 min, respectively, where up to 100 % of ethane was converted into methane and coke through cracking/hydrogenolysis. Furthermore, based on the generalized power-law equation model, a modified model for deactivation kinetic was proposed. Consequently, the influence of the metal function on the kinetic parameters was discussed. We anticipate the formation of Ni 3 Ga intermetallic alloy and its synergy with the exchanged Ga δ+ to be responsible for optimal BTX formation.

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

An SYK-inspired model with density–density interactions: Spectral & wave function statistics, Green’s function and phase diagram

The Sachdev–Ye–Kitaev (SYK) model is a rare example of a strongly-interacting system that is analytically tractable. Tractability arises because the model is largely structureless by design and therefore artificial: while the interaction is restricted to two-body terms, interaction matrix elements are “randomized” and therefore the corresponding interaction operator does not commute with the local density. Unlike conventional density–density-type interactions, the SYK-interaction is, in this sense, not integrable. We here investigate a variant of the (complex) SYK model, which restores this integrability. It features a randomized single-body term and a density–density-type interaction. We present numerical investigations suggesting that the model exhibits two integrable phases separated by several intermediate phases including a chaotic one. The chaotic phase carries several characteristic SYK-signatures including in the spectral statistics and the frequency scaling of the Green’s function and therefore should be adiabatically connected to the non-Fermi liquid phase of the original SYK model. Thus, our model Hamiltonian provides a bridge from the SYK-model towards microscopic realism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

AGGREGATE: dAta-driven modelinG preservinG contRollable dEr for outaGe mAnagemenT and rEsiliency (Report for Task 13: Functional Specification Report(Deliverable D9))

This report provides functional specification for the AGGREGATE project and associated modules/ tools being developed for outage management and restoration with the Distributed Energy Resources (DERs). This report also provide initial validation of the developed modules and performance metrics. The AGGREGATE team at this stage has been collaborating with General Electric (GE) and Seattle City Light (SCL) to move forward from offline validation to online validation using commercial advanced distribution management system (ADMS). Hence, the functional specifications documents, which outline the requirement of various modules, is provided to address the integration of the modules, including 1) real-time model update and aggregation tools module, 2) observability and controllability metric module, 3) fault location solution service restoration module, 4) transmission-distribution co-simulation module, and 5) advanced distribution management system module. integration of different AGGREGATE modules and functional specification provides the foundation for online multi-scenario testing using SCL system model.

42 ENGINEERING↗
DOE OSTI2021

Formation of a Ti–Cu(111) single atom alloy: Structure and CO binding

A single atom Ti–Cu(111) surface alloy can be generated by depositing small amounts of Ti onto Cu(111) at slightly elevated surface temperatures (~500 to 600 K). Here, scanning tunneling microscopy shows that small Ti-rich islands covered by a Cu single layer form preferentially on ascending step edges of Cu(111) during Ti deposition below about 400 K but that a Ti–Cu(111) alloy replaces these small islands during deposition between 500 and 600 K, producing an alloy in the brims of the steps. Larger partially Cu-covered Ti-containing islands also form on the Cu(111) terraces at temperatures between 300 and 700 K. After surface exposure to CO at low temperatures, reflection absorption infrared spectroscopy (RAIRS) reveals distinct C–O stretch bands at 2102 and 2050 cm –1 attributed to CO adsorbed on Cu-covered Ti-containing domains vs sites in the Ti–Cu(111) surface alloy. Calculations using density functional theory (DFT) suggest that the lower frequency C–O stretch band originates specifically from CO adsorbed on isolated Ti atoms in the Ti–Cu(111) surface alloy and predicts a higher C–O stretch frequency for CO adsorbed on Cu above subsurface Ti ensembles. DFT further predicts that CO preferentially adsorbs in flat-lying configurations on contiguous Ti surface structures with more than one Ti atom and thus that CO adsorbed on such structures should not be observed with RAIRS. The ability to generate a single atom Ti–Cu(111) alloy will provide future opportunities to investigate the surface chemistry promoted by a representative early transition metal dopant on a Cu(111) host surface.

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

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Plant diversity and functional identity drive grassland rhizobacterial community responses after 15 years of CO 2 and nitrogen enrichment

Abstract Improved understanding of bacterial community responses to multiple environmental filters over long time periods is a fundamental step to develop mechanistic explanations of plant–bacterial interactions as environmental change progresses. This is the first study to examine responses of grassland root‐associated bacterial communities to 15 years of experimental manipulations of plant species richness, functional group and factorial enrichment of atmospheric CO 2 (eCO 2 ) and soil nitrogen (+N). Across the experiment, plant species richness was the strongest predictor of rhizobacterial community composition, followed by +N, with no observed effect of eCO 2 . Monocultures of C 3 and C 4 grasses and legumes all exhibited dissimilar rhizobacterial communities within and among those groups. Functional responses were also dependent on plant functional group, where N 2 ‐fixation genes, NO 3− ‐reducing genes and P‐solubilizing predicted gene abundances increased under resource‐enriched conditions for grasses, but generally declined for legumes. In diverse plots with 16 plant species, the interaction of eCO 2 +N altered rhizobacterial composition, while +N increased the predicted abundance of nitrogenase‐encoding genes, and eCO 2 +N increased the predicted abundance of bacterial P‐solubilizing genes. Synthesis : Our findings suggest that rhizobacterial community structure and function will be affected by important global environmental change factors such as eCO 2 , but these responses are primarily contingent on plant species richness and the selective influence of different plant functional groups.

Revillini, Daniel↗
DOE OSTI2022

Spectral Filtering Induced by Non-Hermitian Evolution with Balanced Gain and Loss: Enhancing Quantum Chaos

Here, the dynamical signatures of quantum chaos in an isolated system are captured by the spectral form factor, which exhibits as a function of time a dip, a ramp, and a plateau, with the ramp being governed by the correlations in the level spacing distribution. While decoherence generally suppresses these dynamical signatures, the nonlinear non-Hermitian evolution with balanced gain and loss (BGL) in an energy-dephasing scenario can enhance manifestations of quantum chaos. In the Sachdev-Ye-Kitaev model and random matrix Hamiltonians, BGL increases the span of the ramp, lowering the dip as well as the value of the plateau, providing an experimentally realizable physical mechanism for spectral filtering. The chaos enhancement due to BGL is optimal over a family of filter functions that can be engineered with fluctuating Hamiltonians.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

On the Rapid Calculation of Binding Affinities for Antigen and Antibody Design and Affinity Maturation Simulations

The accurate and efficient calculation of protein-protein binding affinities is an essential component in antibody and antigen design and optimization, and in computer modeling of antibody affinity maturation. Such calculations remain challenging despite advances in computer hardware and algorithms, primarily because proteins are flexible molecules, and thus, require explicit or implicit incorporation of multiple conformational states into the computational procedure. The astronomical size of the amino acid sequence space further compounds the challenge by requiring predictions to be computed within a short time so that many sequence variants can be tested. In this study, we compare three classes of methods for antibody/antigen (Ab/Ag) binding affinity calculations: (i) a method that relies on the physical separation of the Ab/Ag complex in equilibrium molecular dynamics (MD) simulations, (ii) a collection of 18 scoring functions that act on an ensemble of structures created using homology modeling software, and (iii) methods based on the molecular mechanics-generalized Born surface area (MM-GBSA) energy decomposition, in which the individual contributions of the energy terms are scaled to optimize agreement with the experiment. When applied to a set of 49 antibody mutations in two Ab/HIV gp120 complexes, all of the methods are found to have modest accuracy, with the highest Pearson correlations reaching about 0.6. In particular, the most computationally intensive method, i.e., MD simulation, did not outperform several scoring functions. The optimized energy decomposition methods provided marginally higher accuracy, but at the expense of requiring experimental data for parametrization. Within each method class, we examined the effect of the number of independent computational replicates, i.e., modeled structures or reinitialized MD simulations, on the prediction accuracy. We suggest using about ten modeled structures for scoring methods, and about five simulation replicates for MD simulations as a rule of thumb for obtaining reasonable convergence. We anticipate that our study will be a useful resource for practitioners working to incorporate binding affinity calculations within their protein design and optimization process.

59 BASIC BIOLOGICAL SCIENCES↗