Search NASA⌕ Search

SEARCH · Search NASA

Results for “reparameterization”

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.

12 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Machine learning for reparameterization of multi-scale closures

Scientific machine learning (ML) is becoming increasingly useful in learning closure models for multi-scale physics problems; however, many ML approaches require a vast array of training data and can struggle with generalization and interpretability. Here, rather than learning an entire closure operator, we adopt an existing reduced-dimension model of the microphysics and learn an optimal re-parameterization of the solver. We demonstrate two approaches for training the reduced dimension closure model (1) an a priori method that optimizes the closure parameterization and the neural network parameters separately and (2) an a posteriori method that simultaneously optimizes both. Using the simulation of biomass pyrolysis as a motivating example, we show that the a posteriori method achieves better target losses and is less dependent on training dataset size for generalizability. We then demonstrate the impact that implementing this reparameterization has at the macroscale, showing improved predictive performance with no modification to the underlying macroscale solvers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Alleviating prior dependencies for DESI DR1 clustering fits through reparameterization

Bayesian analyses of the full-shape clustering of Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) exhibit prior-volume projection effects, whereby weakly constrained nuisance parameters of the Effective Field Theory of Large Scale Structure (EFTofLSS) shift marginalized cosmological posteriors away from the posterior maximum. We reanalyze DESI DR1 power spectrum multipoles using two complementary mitigation strategies: (i) nonlinear orthogonalization to decorrelate nuisance and cosmological parameter priors, and (ii) a fully reparameterization-invariant Jeffreys prior over all EFTofLSS coefficients, evaluated on-the-fly via closed-form Jacobians. Including data from DESI, Big-Bang Nuclesynthesis and a constraint on $n_{\mathrm{s}}$, baseline priors lead to multi-$σ$ projection in the Hubble parameter $H_{0}$ and dark energy equation of state parameters $w_{0}$ and $w_{a}$; the Jeffreys prior successfully recenters these posteriors to enclose the maximum a posteriori estimate within the 68% credible regions, demonstrating clear mitigation of projection effects for these late-time expansion parameters. A hybrid Jeffreys+baseline-Gaussian configuration controls residual over-broad tails in the physical cold dark matter density $ω_{\mathrm{c}}$ while preserving the volume correction, and is our favoured approach. We compare the credible intervals derived using our methodology to those obtained using Halo Occupation Distribution (HOD)-informed priors and to confidence intervals derived using frequentist profile likelihood analyses, finding agreement in both central values and degeneracy directions in the $w_{0}$--$w_{a}$ plane. This demonstrates that, once projection effects are properly controlled, we can make robust inferences about the late-time cosmological expansion independent of the statistical framework adopted.

Bonici, M. [Waterloo U.; Perimeter Inst. Theor. Ph↗
DOE OSTI2026

Reparameterizing a Lipid Force Field Using Small-Angle X-ray Scattering to Improve Predictions of Multicomponent Membranes under Organic Solvent Stress

Understanding and predicting properties of lipid-bilayer membranes are essential to elucidating organismal physiology and pathophysiology. Therefore, substantial efforts have been undertaken to derive accurate molecular-mechanics force fields (FF) to allow simulation of their properties; however, much of these past efforts focused on tuning force fields to accurately reproduce the properties of model single-component membranes and not multicomponent or stressed membranes. Here, we tuned the CHARMM36 FF for a simple 2-component model of a Gram-positive bacterium. This updated force field is parametrized against a single condensed-phase property, the small-angle X-ray scattering (SAXS) intensities of the 2-component lipid system, using a version of ForceBalance implemented previously (named ForceBalance-SAS) with SAXS and small-angle neutron scattering (SANS) intensities as optimization targets. After tuning, we observe improved agreement with the experimental SAXS of the membrane under 1-butanol and tetrahydrofuran solvent stresses, with a factor of 10.0 and 7.5 reduction, respectively, in χ 2 , the measure of the discrepancy between experimental and computed SAXS intensities. Furthermore, this reparametrized force field yielded improved agreement between experimental and simulated membrane thicknesses for the pure lipid systems. However, we note limitations in transferability and diagnose the source of such limitations, particularly the need for SANS in addition to SAXS.

Demerdash, Omar N. A. [Oak Ridge National Laborato↗
DOE OSTI2024

Mixed gauge-global symmetries, elliptic modes, and black hole thermodynamics in Hořava-Lifshitz gravity

In Hořava-Lifshitz gravity, a putative consistent theory of quantum gravity for which there is evidence for both black hole thermodynamics and a holographic construction, spacetime is endowed with a preferred dynamical spacelike foliation. The theory has a leaf reparameterization symmetry that is neither global nor local gauge, hyperbolic and elliptic equations of motion, a lack of splittability, and universal horizon black hole solutions. The reparameterization symmetry is “mixed”: it is a local symmetry in one coordinate yet global on each leaf. More broadly it is an example of both unfree and projectable gauge symmetries. The mixed symmetry and associated charge has not yet been accounted for in calculations of universal horizon thermodynamics in Hořava-Lifshitz gravity. This has led to problems, in particular the failure of the first law in a class of asymptotically AdS solutions where the normal to the leaves of the foliation is not aligned with the time translation Killing vector at infinity. We show how the dynamics of the charge corresponding to this symmetry coupled with the other features above resolves this issue. We then briefly comment how this mixed symmetry, the corresponding charge, and the elliptic equations of motion also conspire to evade recent holographic arguments for only local gauge fields in consistent theories of quantum gravity due to the lack of splittability of the elliptic equation and associated mode.

Global Symmetries↗
DOE OSTI2025

Off-shell vertices in heavy particle effective theories and B → Dπℓν

We study the modifications to decay amplitudes in heavy to heavy semileptonic decays with multiple hadrons in the final state due to intermediate heavy hadrons being off-shell or having a finite width. Combining Heavy Hadron Chiral Perturbation Theory (HHχPT) with a BCFW on-shell factorization formula, we show that these effects induce O (1/M) corrections to the standard results computed in the narrow-width approximation and therefore are important in extracting form factors from data. A combination of perturbative unitarity, analyticity, and reparameterization invariance fully determine these corrections in terms of known Isgur-Wise functions without the need to introduce new form factors. In doing so, we develop a novel technique to compute the boundary term at complex infinity in the BCFW formula for theories with derivatively coupled scalars. While we have used the $\overline{B}$ → Dπℓν decay as an example, these techniques can generally be applied to effective field theories with (multiple) distinct reference vectors. Article PDF

Chiral Lagrangian↗
DOE OSTI2025

Atomistic Insights into the Reactive Diffusion of CO 2 in Guanidine-Based Facilitated Transport Membranes

The pressing need to address climate change has led to significant advancements in carbon dioxide (CO 2 ) capture technologies. Notably, facilitated transport membranes (FTMs) are distinguished by their exceptional selectivity and permeance, attributed to their reversible chemical reactions with CO 2 . This study, for the first time, sheds light on the reactive diffusion mechanism of CO 2 in FTMs, utilizing 1,1,3,3-tetramethylguanidine (TMG) as a mobile carrier. Specifically, state-ofthe- art molecular dynamics (MD) simulations, augmented by a reparameterized reactive force field (ReaxFF) capable of describing atomistic interactions and reaction pathways, are conducted to investigate the transport of CO 2 in TMG. The analysis of mean squared displacement (MSD) and diffusion coefficients reveals a clear hierarchy in the mobility of reaction components. Our findings highlight a unique hopping diffusion mechanism between bicarbonate ions and TMG molecules, increasing the diffusivity of reacted CO 2 by 1.4 times. The hopping events observed not only enhance our understanding of molecular mobility but also offer a means to boost the performance of FTMs in CO 2 capture applications. Overall, this research lays the groundwork for the future design of FTMs with optimal carrier properties.

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

Prediction and Experimental Verification of Electrolyte Solvation Structure from an OMol25-Trained Interatomic Potential

A molecular-level understanding of electrolyte solvation structure and ion–ion correlations is critical to developing next-generation battery chemistries. Atomistic simulation capabilities with sufficient accuracy, speed, and transferability to deliver reliable structural insights while avoiding arduous system-specific reparameterization are thus highly desirable. Machine learning interatomic potentials (MLIPs) trained on large, chemically diverse data sets are revolutionizing computational chemistry, enabling molecular dynamics simulations of battery electrolytes with near-DFT accuracy over 10,000× faster than DFT. While previous MLIP training data sets with suitable elemental coverage for electrolytes have been based on inorganic materials, the Open Molecules 2025 (OMol25) data set provides large-scale molecular DFT MLIP training data with broad elemental coverage and specifically samples tens of millions of electrolyte configurations. Here, we integrate computational modeling with experimental validation to systematically assess the ability of large-scale MLIPs pretrained on materials data or on OMol25 to accurately resolve nanoscale structural organization and ion-solvation characteristics in Na-ion battery electrolytes across diverse physicochemical conditions and compositional regimes. We find that the OMol25-trained Universal Model of Atoms (UMA-OMol) predicts experimentally measured densities and X-ray structure factors in substantially better agreement compared to state-of-the-art models trained only on inorganic materials data. Using UMA-OMol, we further analyze systematic trends in solvation structure as a function of cation identity, anion chemistry, salt concentration, and solvent topology. We observe that increasing system temperature amplifies the heterogeneity within the solvation environment, perturbing cation–solvent interactions and promoting the formation of contact ion pairs (CIPs). Moreover, subtle variations in the solvent topology of glyme-based electrolytes cause pronounced changes in ion correlations and solvation structure. The experimental agreement and microscopic insights shown here position OMol25-trained MLIPs as a practical route to predictive, high-throughput electrolyte simulations beyond the limits of classical force fields and direct DFT molecular dynamics, serving as a powerful tool for accelerating the design of next-generation Na-ion battery electrolytes and beyond.

MLIPs↗
DOE OSTI2025

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗
DOE OSTI2024

A modular and extensible CHARMM-compatible model for all-atom simulation of polypeptoids

Peptoids (N-substituted glycines) are a class of sequence-defined synthetic peptidomimetic polymers with applications including drug delivery, catalysis, and biomimicry. Classical molecular simulations have been used to predict and understand the conformational dynamics of single chains and their self-assembly into morphologies including sheets, tubes, spheres, and fibrils. The CGenFF-NTOID model based on the CHARMM General Force Field has demonstrated success in accurate all-atom molecular modeling of peptoid structure and thermodynamics. Extension of this force field to new peptoid side chains has historically required reparameterization of side chain bonded interactions against ab initio data. This fitting protocol improves the accuracy of the force field but is also burdensome and precludes modular extensibility of the model to arbitrary peptoid sequences. In this work, we develop and demonstrate a Modular Side Chain CGenFF-NTOID (MoSiC-CGenFF-NTOID) as an extension of CGenFF-NTOID employing a modular decomposition of the peptoid backbone and side chain parameterizations, wherein arbitrary side chains within the large family of substituted methyl groups (i.e., –CH 3 , –CH 2 R, –CHRR', and –CRR'R") are directly ported from CGenFF. We validate this approach against ab initio calculations and experimental data to develop a MoSiC-CGenFF-NTOID model for all 20 natural amino acid side chains along with 13 commonly used synthetic side chains and present an extensible paradigm to efficiently determine whether a novel side chain can be directly incorporated into the model or whether refitting of the CGenFF parameters is warranted. We make the model freely available to the community along with a tool to perform automated initial structure generation.

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

Switchgrass ( Panicum Virgatum ) and Miscanthus ( Miscanthus × Giganteus ) Long-Term Yield Patterns Reveal Consistent Productivity Declines

Perennial grasses like switchgrass ( Panicum virgatum ) and miscanthus ( Miscanthus × giganteus ) are expected to supply a substantial amount of the United States bioeconomy's feedstock demand. However, uncertainties around their long-term yields challenge the viability of their potential and limit their wider adoption. To resolve their long-term yield patterns, we analyzed over 200 plantings of switchgrass and miscanthus across Michigan and Wisconsin, USA, measured over 5–15 years. We found a consistent two-phase long-term yield dynamic; during a yield-building phase , peak yields occurred within 4–5 years after planting, followed by a yield-decline phase in which switchgrass and miscanthus lost 30%–47% and 14%–40% of peak yields, respectively. Among the potential drivers of this dynamic and the yield decline, we found that weather conditions had little impact, as the variation across years was not large enough to drive the observed yield differences. Added nitrogen increased peak yields by 10%–20% and attenuated the yield decline by 20%–50%. However, since fertilized stands still showed a yield decline, other factors became limiting as stands aged. This conserved long-term yield dynamic has direct implications on management. A farm-to-gate economic analysis suggests replanting switchgrass and miscanthus 5 and 9 years following their peak yields maximizes profit over a 30-year time horizon. Results call for further management and breeding strategies to mitigate the yield-decline phase, and for reparameterization of global bioenergy models with carbon capture and storage, which may overestimate yields and the economic and environmental benefits of crops grown for bioenergy feedstocks.

bioenergy↗
NASA NTRS1981

Software models: A Bayesian approach to parameter estimation in the Jelenski-Moranda software reliability model

Maximum likelihood estimation procedures for the Jelinski-Moranda software reliability model often give misleading answers. A reparameterization and a Bayesian analysis eliminate some of the problems incurred by MLE methods and often give better predictions on sets of real and simulated data. Practical difficulties in estimating the initial number of errors N and the failure rate of each error phi by the method of maximum likelihood are: N, the MLE of N, is occasionally infinite (i.e., the routines for calculating N and phi do not converge). It is shown that N is finite sub i only if the regression line of the interevent times t sub i vs. i has positive slope. A serious problem is that often N approximates n, the sample size, and sometimes N = n. Thus the MLE predicts that the program is perfect even when it is far from being so. Only when almost all failures have been removed can N and phi be trusted near the end of debugging.

Littlewood, B.↗
DOE OSTI2025

Energy efficient integrated photonic systems based on inverse design

The energy footprint of modern information processing and communications systems is immense and utilizes a significant fraction of total global energy usage. Data centers alone consume over 70 billion kilowatt-hours per year. Much of this energy usage is intrinsic to the use of electronic wiring, making optical-based technologies a necessary and promising route to mitigating energy consumption in short- to medium-distance communication links. In this project, we developed and implemented a framework for the design of optical components, based on machine learning, which enables optical components relevant to optical information processing to be realized at their physical performance limits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗