Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Accelerating Instanton Theory with the Line Integral Nudged Elastic Band Method and Gaussian Process Regression

Quantum tunneling plays a fundamental role in many chemical reactions, particularly proton transfer processes. Ring polymer instanton theory offers a practical framework for computing tunneling rates in complex molecular systems. However, applying the ring polymer instanton method with a potential energy surface generated on-the-fly using electronic structure calculations can be computationally demanding. Here, in this work, we present a new efficient implementation of the ring polymer instanton method by combining the Line Integral Nudged Elastic Band (LI-NEB) approach with Gaussian Process Regression (GPR). We benchmarked this method on prototypical ground-state proton transfer systems, including the benchmark gas-phase hydrogen abstraction reaction H + CH 4 → H 2 + CH 3 , malonaldehyde, and Z-3-amino-propenal (aminopropenal). Our results show that this approach is an order of magnitude faster than traditional instanton algorithms while maintaining excellent agreement with their tunneling rates. This development opens the door to studying proton transfer in larger systems with improved efficiency.

chemical physics↗
DOE OSTI2025

Many-Body Basis Set Amelioration Method for Incremental Full Configuration Interaction

Incremental full configuration interaction (iFCI) is a polynomial-cost electronic structure method that systematically approaches the FCI limit by employing the method of increments to solve the Schrödinger equation through a many-body expansion. This article introduces the many-body basis set amelioration (MBBSA) method, which is designed to allow iFCI to be applicable to larger atomic orbital basis sets. MBBSA uses a series of inexpensive iFCI calculations to approximate the correlation energy that would be found using a more expensive, highly accurate iFCI calculation. Here, when compared to standard iFCI computations on smaller molecules in triple-zeta and larger basis sets, MBBSA provides approximations to the total and relative energies within chemical accuracy. MBBSA exhibits a reduced cost of between 60-92% when compared to standard iFCI calculations, with larger systems experiencing the largest benefit. Tests of MBBSA on two reactions that involve highly correlated systems, the automerization of cyclobutadiene and a Criegee intermediate reaction, show that MBBSA has practical utility for studying realistic chemistries.

Basis sets↗
DOE OSTI2025

Lowering the Scaling of Self-Consistent Field Methods by Combining Tensor Hypercontraction and a Density Difference Ansatz

We present the tensor hypercontraction difference self-consistent field (SCF) method, an approach that reduces the formal computational scaling of traditional naive self-consistent field methods from 𝑂(𝑁 4 ) to 𝑂(𝑁 3 ) with system size 𝑁. The scaling reduction is achieved by developing a new technique for constructing the tensor hypercontraction decomposition based on the fundamental approximation made in density fitting. Combining this scheme with the difference self-consistent field methodology, we achieve a method that enables 𝑂(𝑁 3 ) scaling SCF calculations with only 𝑂(𝑁 2 ) storage requirements. In conclusion, our proof-of-concept numerical tests demonstrate robust performance with errors in total energies below 8 × 10 –4 E h and with sub 1 kcal/mol errors for relative energies.

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

Aggregation Methods for Quantifying PTM and Structural Changes in Bottom-Up Proteomics

Bottom-up proteomic workflows rely on sequential preprocessing steps, commonly including peptide-to-protein aggregation (“roll-up”), to enhance data reliability and interpretability. While roll-up is effective for protein-centered analyses, it may be suboptimal for applications focused on post-translational modifications (PTMs) or protein structural changes, such as limited proteolysis–mass spectrometry (LiP-MS). Here, we investigate how different roll-up strategies influence site-level quantification in PTM differential analysis. Moreover, we introduce a novel site-centric roll-up approach tailored for LiP-MS, which quantifies proteolytic fragments rather than solely tryptic peptides. We benchmark these methods through simulation studies, comparing their sensitivity and specificity in detecting structural and PTM-driven changes. We found that the median and mean roll-up methods outperform the sum method in both PTM and LiP proteomics, and site-level quantification in LiP outperforms peptide-level quantification. Our findings offer the first systematic, data-driven guidance for selecting roll-up techniques in site-level proteomic analyses, with implications for both PTM-focused and structural proteomics studies.

aggregation↗
DOE OSTI2024

A Novel Method to Train Classification Models for Structure Detection in In Situ Spacecraft Data

We present a method for creating spacecraft-like data which can be used to train Machine Learning (ML) models to detect and classify structures in in situ spacecraft data. First, we use the Grad-Shafranov equation to numerically solve for several magnetohydrostatic equilibria which are variations on a known analytic equilibrium. These equilibria are then used as the initial conditions for Particle-In-Cell simulations in which the structures of interest are observed and labeled. We then take one-dimensional slices through the simulations to replicate what a spacecraft collecting data from the simulation would observe. This sliced data then can be used as training data for the initial training of ML models intended for use on spacecraft data. We demonstrate the method applied to the problem of detecting small-scale plasmoids in the magnetotail, which is important for understanding complex magnetotail reconnection dynamics. The simple 1D classifier we train is able to detect more than 70% of the plasmoid points in the data set but also produces a large number of false positives. Our further work on this example problem is detailed, and further potential uses of the method are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Single Compartment Relaxed Eddy Accumulation Method

The relaxed eddy accumulation (REA) method is a widely-known technique that measures turbulent fluxes of scalar quantities. The REA technique has been used to measure turbulent fluxes of various compounds, such as methane, ethene, propene, butene, isoprene, nitrous oxides, ozone, and others. The REA method requires the accumulation of scalar concentrations in two separate compartments that conditionally sample updrafts and downdraft events. It is demonstrated here that the assumptions behind the conventional or two-compartment REA approach allow for one-compartment sampling, therefore called a one compartment or 1-C-REA approach, thereby expanding its operational utility. The one-compartment sampling method is tested across various land cover types and atmospheric stability conditions, and it is found that the one-compartment REA can provide results comparable to those determined from conventional two-compartment REA. This finding enables rapid expansion and practical utility of REA in studies of surface-atmosphere exchanges, interactions, and feedbacks.

Banerjee, T.↗
DOE OSTI2025

Hot or Not? An Evaluation of Methods for Identifying Hot Moments of Nitrous Oxide Emissions From Soils

Abstract Effectively quantifying hot moments of nitrous oxide (N 2 O) emissions from agricultural soils is critical for managing this potent greenhouse gas. However, we are challenged by a lack of standard approaches for identifying hot moments, including (a) determining thresholds above which emissions are considered hot moments, and (b) considering seasonal variation in the magnitude and frequency distribution of net N 2 O fluxes. We used one year of hourly N 2 O flux measurements from 16 autochambers that varied in flux magnitude and frequency distribution in a conventionally tilled maize field in central Illinois, USA, to compare three approaches to identify hot moment thresholds: standard deviations (SD) above the mean, 1.5x the interquartile range (IQR), and isolation forest (IF) identification of anomalous values. We also compared these approaches on seasonally subdivided data (early, late, and non‐growing seasons) versus the whole year. Our analyses revealed that 1.5x IQR method best identified N 2 O hot moments. In contrast, using 2 or 4 SD both yielded hot moment threshold values too high, and IF yielded threshold values too low, leading to missed N 2 O hot moments or low net N 2 O fluxes mischaracterized as hot moments, respectively. Furthermore, seasonally subdividing the data set not only facilitated identification of smaller hot moments in the late‐ and non‐growing seasons when N 2 O hot moments were generally smaller but it also increased hot moment threshold values in the early growing season when N 2 O hot moments were larger. Consequently, of the methods evaluated here, we recommend using the 1.5x IQR method on whole year data sets to identify N 2 O hot moments.

Stuchiner, Emily R. [Institute for Sustainability,↗
DOE OSTI2024

A Diffusion‐Based Uncertainty Quantification Method to Advance E3SM Land Model Calibration

Abstract Calibrating land surface models and accurately quantifying their uncertainty are crucial for improving the reliability of simulations of complex environmental processes. This, in turn, advances our predictive understanding of ecosystems and supports climate‐resilient decision‐making. Traditional calibration methods, however, face challenges of high computational costs and difficulties in accurately quantifying parameter uncertainties. To address these issues, we develop a diffusion‐based uncertainty quantification (DBUQ) method. Unlike conventional generative diffusion methods, which are computationally expensive and memory‐intensive, DBUQ innovates by formulating a parameterized generative model and approximates this model through supervised learning, which enables quick generation of parameter posterior samples to quantify its uncertainty. DBUQ is effective, efficient, and general‐purpose, making it suitable for site‐specific ecosystem model calibration and broadly applicable for parameter uncertainty quantification across various earth system models. In this study, we applied DBUQ to calibrate the Energy Exascale Earth System Model land model at the Missouri Ozark AmeriFlux forest site. Results indicated that DBUQ produced accurate parameter posterior distributions similar to those from Markov Chain Monte Carlo sampling but with 30 times less computing time. This significant improvement in efficiency suggests that DBUQ can enable rapid, site‐level model calibration at a global scale, enhancing our predictive understanding of climate impacts on terrestrial ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Exponential concentration in quantum kernel methods

Kernel methods in Quantum Machine Learning (QML) have recently gained significant attention as a potential candidate for achieving a quantum advantage in data analysis. Among other attractive properties, when training a kernel-based model one is guaranteed to find the optimal model’s parameters due to the convexity of the training landscape. However, this is based on the assumption that the quantum kernel can be efficiently obtained from quantum hardware. In this work we study the performance of quantum kernel models from the perspective of the resources needed to accurately estimate kernel values. We show that, under certain conditions, values of quantum kernels over different input data can be exponentially concentrated (in the number of qubits) towards some fixed value. Thus on training with a polynomial number of measurements, one ends up with a trivial model where the predictions on unseen inputs are independent of the input data. We identify four sources that can lead to concentration including expressivity of data embedding, global measurements, entanglement and noise. For each source, an associated concentration bound of quantum kernels is analytically derived. Lastly, we show that when dealing with classical data, training a parametrized data embedding with a kernel alignment method is also susceptible to exponential concentration. Our results are verified through numerical simulations for several QML tasks. Altogether, we provide guidelines indicating that certain features should be avoided to ensure the efficient evaluation of quantum kernels and so the performance of quantum kernel methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Passive high explosive neutron inspection (PHENIX): a new method to confirm the presence or absence of high explosives for nuclear treaty verification

Advanced instruments and methods need to be developed now to create a technical basis to support the negotiation of future nuclear arms control treaties. One new capability that is anticipated is the ability to confirm either the declared presence or declared absence of high explosive (HE) material in the presence of special nuclear material (SNM). Towards this goal, Passive HE Neutron Inspection (PHENIX) has been developed and demonstrated as a method for confirming the presence or absence of HE in the presence of plutonium. The method exploits the inherent presence of neutrons associated with the decay of plutonium as an internal probe source for performing prompt gamma-ray neutron activation analysis (PGNAA), searching for the presence of HE as revealed by the emission of characteristic gamma rays following neutron absorption in hydrogen and nitrogen which are building blocks of present-day, military-grade HE. Tests using stoichiometrically-correct hemishells of mock HE with plutonium show that a system can be expected to positively confirm the presence or absence of these signatures, supporting determination of HE presence or absence with Pu, in a few hours. To protect other potentially sensitive gamma-ray signatures from a treaty accountable item, an analog information barrier has been conceptualized and tested which physically prevents the collection of gamma-ray spectral data outside of user selected energy windows strategically chosen to view only narrow spectral regions corresponding to the hydrogen (2223.2 keV) and nitrogen (9807.2 keV, 10,318.2 keV, and 10,829.2 keV) PGNAA signatures.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

A simple and practical wax-encapsulation method for air-sensitive XAS samples

To facilitate X-ray absorption spectroscopy (XAS) measurements of air-sensitive samples, we present a simple method in which materials are encased in common paraffin wax to protect them from air and moisture. We demonstrate the efficacy of this approach using a highly reducing, air- and moisture-sensitive uranium(III) complex, the tris(amide) U[N(SiMe 3 ) 2 ] 3 (1). When finely dispersed in a boron nitride matrix and subsequently encased in inert paraffin wax, samples of 1 remain stable with no visible or spectroscopic degradation after several days under ambient conditions. The viability of this method for XAS measurements was further evaluated across a series of uranium compounds, ranging from uranyl species to highly air- and moisture-sensitive molecular complexes, at the uranium L 3 -edge. Edge energy determinations were highly reproducible (±0.1 eV between replicates) and, where available, showed excellent agreement with literature values. In conclusion, this low-cost, effective, and versatile method offers a viable solution for XAS studies of air-sensitive compounds and materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A novel method may reveal bulk metallic glass compressive ductility trends in high data rate nanoindentation

Recent methods allow novel amorphous alloy compositions to be rapidly manufactured at small scale; however, obtaining materials properties such as compressive ductility from these smaller specimens has remained a challenge. Here, we suggest a potential high-throughput nanoindentation method that may be able to rapidly characterize the relative compressive ductility between these alloys based on their serration characteristics. The properties of emergent serrations, when interpreted in a simple micromechanical stress relaxation model, may order these materials by their compressive plastic strain to failure. These results are consistent with the ordering obtained from compressed specimens as well as with model simulations, suggesting that this model may be broadly useful for interpreting compressive ductility from nanoindentation serrations. After it is validated on more materials, this new method will match the rapid pace of amorphous alloy development, thus allowing metallic glass properties to be fine-tuned for each application prior to scale prototyping.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Jacobian-free pseudo-arclength continuation method for phase transitions in inhomogeneous thermodynamic systems

Developing phase diagrams for inhomogeneous systems in thermodynamics is difficult, in part, due to the large phase space and the possibility of unstable and metastable solutions arising from first-order phase transitions. Pseudo-arclength continuation (PAC) is a method that allows one to trace out stable and unstable solutions of nonlinear systems. Typically, PAC utilizes the Jacobian in order to implement Newton (or quasi-Newton) steps. In this work, we present a Jacobian-free PAC method that is amenable to the usual workflows in inhomogeneous thermodynamics. We demonstrate our method in systems that have first-order phase transitions, including a novel example of polyelectrolyte complex coacervation in confinement, where multiple surface phase transitions occur and can overlap with one another.

Chemistry↗
DOE OSTI2024

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗
DOE OSTI2024

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING↗
DOE OSTI2025

Measuring the thermal conductivity of hydrogels with a bidirectional 3w method

Hydrogels are soft, water-absorbing polymer materials with diverse applications in biomedicine and agriculture. Recently, hydrogels have been proposed to encapsulate water-soluble phase change materials which store energy in their latent heat of solidification. In these applications, the thermal conductivity of these materials affects their performance. Few methods exist for measuring the thermal conductivity of small quantities of hydrogels. Here, we describe an implementation of the bidirectional 3w technique to measure the thermal conductivity of hydrogels with particular attention to their moisture content. Our implementation of the technique can probe sample volumes as little as ~20 mL and yields the thermal conductivity without requiring fitting of additional thermal parameters. We numerically simulate 3w sensor designs with frequency-domain 3-D models to quantify and reduce errors introduced by the choice of substrate and insulation layer thickness. Frequencies in the ~1−20 Hz range yield less error for the materials considered here. We verify our setup with measurements on water and report values for polyacrylamide and poly(2-acrylamido-2-methylpropane sulfonic acid) (PAMPS) hydrogels. Our swollen hydrogels exhibited thermal conductivities nearly equivalent to water, 0.6 W m-1 K-1, and we estimate thermal conductivities of 0.43 and 0.42 W m-1 K-1 for neat polyacrylamide and PAMPS, respectively. Finally, we estimate an error of ±7%, consistent with other 3ω methods, with the largest error coming from the sensor calibration. We find our implementation of the bidirectional 3w method gives reasonable results and can be employed for prototyping soft materials relevant for thermal storage.

3-omega, thermal conductivity, hydrogel, moisture ↗
DOE OSTI2025

Symmetry-projected spin-AGP methods applied to spin systems

Symmetry-projected wave function methods capture static correlation by breaking and restoring the symmetries of a system. In this article, we present the symmetry-projected spin antisymmetrized geminal power (spin-AGP) state projected onto space group symmetry as well as complex conjugation, spin-flip, and time-reversal symmetries. The method is benchmarked on the 1D XXZ model and the 2D J 1 − J 2 model with square and triangular lattices. Our results indicate that symmetry projection methods provide a powerful tool for frustrated spin systems.

Antisymmetrized geminal power↗
DOE OSTI2025

Performance of the spin-component-scaled methods for energy bands

The performance of various spin-component-scaled parameterisations is examined for the second-order many-body Green's-function [MBGF(2)] calculations of valence energy bands, taking three of the experimentally well-characterised polymers as examples: polyethylene, polytetrafluoroethylene, and polyacetylene. The parameterisations considered are Grimme's original SCS parameter set, Jung et al.'s original SOS set (retaining the opposite-spin component only), Śmiga et al.'s SCS(IP) set (calibrated specifically for ionization energies), and Śmiga et al.'s SOS(IP) set (calibrated for ionization energies with the opposite-spin component only; implicit in the os-D2 model of Opoku et al.). The SCS(IP) and SOS(IP) parameterisations are found to shift both outer and inner valence bands by up to a few electronvolts away from the experimental data. The original SCS and SOS parameter sets do not improve upon, but largely maintain the accuracy of the unscaled MBGF(2) methods. Given that the SOS-MBGF(2) method can be implemented in a quartic-scaling algorithm (for all roots), it is most promising for solid-state applications. Furthermore this observation is consistent with the success of the quartic-scaling GW methods without the vertex correction based on a density-functional theory reference.

Green’s-function theory↗