Search NASA⌕ Search

SEARCH · Search NASA

Results for “sensitivity theory”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Implementation of Perturbation Theory and Sensitivity Capabilities in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor Multiphysics analysis application, jointly developed by Argonne and Idaho National Laboratories under the DOE-NE NEAMS program. This fiscal year, capabilities for reactivity and sensitivity evaluation using perturbation methods were implemented and verified. The First Order Perturbation Method (FOPT) was employed to compute reactivity worth resulting from small perturbations in input parameters, while the Generalized Perturbation Theory (GPT) was used to evaluate sensitivities of a range of response types, including reaction rate ratio, k-eigenvalue, neutron generation time, and effective delayed neutron fraction. These perturbation methods enable users to quantify how response quantities change due to a perturbation in a input parameter without explicitly performing an additional transport simulation for each perturbed state. In particular, the GPT formulation accounts for indirect effects arising from flux changes by solving generalized inhomogeneous equations, for which a Neumann series-based iterative solution method was developed and implemented in Griffin. The implemented reactivity and sensitivity evaluation capabilities were verified using two test problems: an infinite homogeneous system and a two-dimensional hexagonal core. The results showed excellent agreement with reference solutions obtained by a direct method based on finite difference approximation as well as GPT-based results from the PERSENT code, confirming the accuracy of both reactivity and sensitivity evaluations. Additionally, preliminary uncertainty quantification (UQ) results were obtained by combining the sensitivity values computed using GPT and external covariance data, demonstrating that the implemented sensitivity results can be reliably used for uncertainty calculations. To further demonstrate the generality and practical strength of the implementation, the sensitivity evaluation capability was successfully applied to the Empire microreactor with a geometrically complex design that poses significant modeling challenges. The results confirm that Griffin enables sensitivity evaluations even for irregular and highly heterogeneous reactor configurations, thereby establishing a foundation for UQ applications in advanced reactor designs and analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

A Theory for the Balance between Warm Rain and Ice Crystal Processes of Precipitation in Mixed-Phase Clouds

Abstract Mixed-phase clouds contain both supercooled cloud liquid and ice crystals. In principle, precipitation may be initiated either by the liquid phase or by the ice phase. Ice crystals may grow by vapor diffusion to become snow (“ice crystal process”), forming “cold” precipitation. Equally, cloud droplets, when large enough, coalesce to form “warm” precipitation by the “warm rain process.” Warm rain could be supercooled and freeze as “warm” graupel. In the present paper, a new simplified theoretical analysis is provided to examine the microphysical system consisting of three species of hydrometeor, namely, cloud liquid, “cold ice” (crystals, snow), and “warm rain” (frozen or supercooled). This is obtained by nondimensionalizing and simplifying the evolution equations for the mass of each species. Analytical formulas are given for equilibria. Feedback analysis shows that the sign of the feedback is linked to the abundance of precipitation, with a neutral surface in the 3D phase space. The system’s precipitation amount explodes while in the initial unstable regime, crossing the neutral surface and approaching the equilibrium point that is a stable attractor. Positive and negative feedbacks are elucidated. In a standard case, the cold ice mass is about 1000 times larger than the warm rain mass. To illustrate the physical behavior of the theory, sensitivity tests are performed with respect to environmental conditions (e.g., aerosol, updraft speed) and microphysical parameters (e.g., riming and sedimentation rates for cold ice). Cold ice prevails, especially in fast ascent, due to its low bulk density, favoring slow sedimentation and a wide cross-sectional area for riming. Significance Statement The theory elucidates how the ice phase can prevail in the precipitation from any mixed-phase clouds with supercooled cloud liquid and crystals. The ice phase radically suppresses cloud liquid by riming when active and “wins” the competition against coalescence. This prevalence of ice is shown to arise from the low bulk density of snow. The cloud is viewed as a system of negative and positive feedbacks that prevail in realms of stability and instability in a 3D phase space.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Tunable magneto-optical properties in MoS 2 via defect-induced exciton transitions

The presence of chalcogen vacancies in monolayer transition metal dichalcogenides (TMDs) leads to excitons with mixed localized-delocalized character and to reduced valley selectivity. Recent experimental advances in defect design in TMDs allow for a close examination of such mixed exciton states as a function of their degree of circular polarization under external magnetic fields, revealing strongly varying defect-induced magnetic properties. A theoretical understanding of these observations and their physical origins demands a predictive, structure-sensitive theory. In this work, we study the effect of chalcogen vacancies on the exciton magnetic properties in monolayer MoS 2 . Using many-body perturbation theory, we show how the complex excitonic picture associated with the presence of defects—with reduced valley and spin selectivity due to hybridized electron-hole transitions—leads to a structurally controllable exciton magnetic response. We find a variety of $\textit{g}$-factors with changing magnitudes and sign depending on the exciton energy and character. Further, our findings suggest a pathway to tune the nature of the excitons—and by that their magneto-optical properties—through defect architecture.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Fractonic coset construction for spontaneously broken translations

We study the homogeneous breaking of spatial translation symmetry concomitantly with the spontaneous breaking of other internal and spacetime symmetries, including dilations. We use the symmetry-breaking pattern as the only input to derive, via the coset construction, general effective field theories for the symmetry-originated modes associated with Goldstone’s theorem, namely the Nambu-Goldstone candidates. Through explicit computations, we show that integrating out the explicit massive Nambu-Goldstone candidates or imposing symmetric constraints, namely the inverse Higgs constraints, to express massive modes in terms of the massless ones leads to physically distinct effective field theories. This sensitivity to the chosen method can be traced back to the homogeneous breaking of translations, the homogeneous aspect of the breaking induces a mixing between internal and spacetime symmetries at the level of the Lie algebra. This, in turn, leads to subtle discussions about the inverse Higgs constraints, in particular that they lead to a loss of generality in our specific examples. The derived general effective field theories also give rise to a broad class of theories exhibiting emergent enhanced shift symmetries, which constrain the mobility of the modes. The latter are referred to as fractonic modes. Published by the American Physical Society 2025

Chavda, Ameya (ORCID:0000000211731605)↗
DOE OSTI2022

An experimental and modeling study on auto-ignition kinetics of ammonia/methanol mixtures at intermediate temperature and high pressure

A rapid compression machine (RCM) has been applied to measure the ignition delay times of NH 3 /CH 3 OH mixtures covering pressures of 20 and 40 bar, equivalence ratios of 0.5, 1.0 and 2.0, and temperatures between 845 and 1100 K. Here the measurements show that the NH 3 /CH 3 OH mixtures become more reactive with increasing methanol addition. Addition of merely 1% (molar basis) of CH 3 OH to NH 3 lowers the ignition temperature around 100 K at 40 bar in comparison to pure NH 3 . The ignition delay is a complex function of fuel mixture and stoichiometry. For the 1% CH 3 OH mixture, the leaner mixtures are more reactive, while the reverse trend is found for mixtures with 5%, 20% and pure CH 3 OH. Analysis of the pressure profiles shows three distinct ignition modes for NH 3/ CH 3 OH mixtures, facilitated by the pre-ignition heat release from NH 3 consumption. A detailed mechanism for ignition of NH 3 /CH 3 OH fuel blends has been developed, capable of reproducing the ignition behavior of mixtures with reasonable accuracy. A subset for amine / methanol interactions was established, with rate constants for the key reaction between NH 2 and CH 3 OH calculated from ab initio theory. A sensitivity analysis indicates that the critical reactions during the auto-ignition process vary with the CH 3 OH mole fraction in the fuel. The ammonia chemistry, namely NH 2 + NO, NH2 + NO 2 and NH 2 + HO 2 , is dominant for the mixture with 1% CH 3 OH, while the reactions related to CH 3 OH and H 2 O 2 are more important for the 20% CH 3 OH mixture. The interaction between ammonia and methanol shows a more prominent effect on auto-ignition for mixtures with 5% CH 3 OH in fuel as compared to those with 1% and 20% CH 3 OH. According to the modeling results, methanol addition is found to enrich the O/H radical pool, consuming ammonia and promoting auto-ignition through different reaction pathways.

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

Revisiting superconductivity in the extended one-band Hubbard model: Pairing via spin and charge fluctuations

The leading superconducting instabilities of the two-dimensional extended repulsive one-band Hubbard model within spin-fluctuation pairing theory depend sensitively on electron density, band, and interaction parameters. We map out the phase diagrams within a random-phase-approximation spin- and charge-fluctuation approach, and find that while B 1g (d x 2 –y 2 ) and B 2g (d xy ) pairing dominates in the absence of repulsive longer-range Coulomb interactions V NN , the latter induces pairing in other symmetry channels, including, e.g., A 2g (g-wave), nodal A 1g (extended s-wave), or nodal E u (p-wave) spin-triplet superconductivity. At the lowest temperatures, transition boundaries in the phase diagrams between symmetry-distinct spin-singlet orders generate complex time-reversal symmetry broken superpositions. By contrast, we find that boundaries between singlet and triplet regions are characterized by first-order transitions. Finally, motivated by recent photoemission experiments, we have determined the influence of an additional explicitly attractive nearest-neighbor interaction, V NN < 0, on the superconducting gap structure. Furthermore, depending on the electronic filling, such an attraction boosts E u (p-wave) spin-triplet or B 1g (d x 2 –y 2 ) spin-singlet ordering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Extreme Risk Mitigation in Reinforcement Learning using Extreme Value Theory

Risk-sensitive reinforcement learning (RL) has garnered significant attention in recent years due to the growing interest in deploying RL agents in real-world scenarios. A critical aspect of risk awareness involves modelling highly rare risk events (rewards) that could potentially lead to catastrophic outcomes. These infrequent occurrences present a formidable challenge for data-driven methods aiming to capture such risky events accurately. While risk-aware RL techniques do exist, they suffer from high variance estimation due to the inherent data scarcity. Our work proposes to enhance the resilience of RL agents when faced with very rare and risky events by focusing on refining the predictions of the extreme values predicted by the state-action value distribution. To achieve this, we formulate the extreme values of the state-action value function distribution as parameterized distributions, drawing inspiration from the principles of extreme value theory (EVT). We propose an extreme value theory based actor-critic approach, namely, Extreme Valued Actor-Critic (EVAC) which effectively addresses the issue of infrequent occurrence by leveraging EVT-based parameterization. Importantly, we theoretically demonstrate the advantages of employing these parameterized distributions in contrast to other risk-averse algorithms. Our evaluations show that the proposed method outperforms other risk averse RL algorithms on a diverse range of benchmark tasks, each encompassing distinct risk scenarios.

Wang, Yu↗
DOE OSTI2022

Sensitivity of time-dependent density functional theory to initial conditions

Time-dependent density-functional theory is mathematically formulated through nonlinear coupled time-dependent three-dimensional partial differential equations, and it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrödinger equation for an interacting many-body system is, however, linear and mathematically the exact equations of the density-functional theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a density-functional theory framework. Whether for realistic implementations of the time-dependent density-functional theory the question of the absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density-functional theory is a unique tool allowing us to study the nonequilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Here our analysis, for a number of quantum superfluid many-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom O(10 10 ) and larger, suggests that its maximum Lyapunov exponents are negligible for all practical purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

42 ENGINEERING↗
DOE OSTI2022

Computational Response Theory for Dynamics

Quantifying the sensitivity - how a quantity of interest (QoI) varies with respect to a parameter – and response – the representation of a QoI as a function of a parameter - of a computer model of a parametric dynamical system is an important and challenging problem. Traditional methods fail in this context since sensitive dependence on initial conditions implies that the sensitivity and response of a QoI may be ill-conditioned or not well-defined. If a chaotic model has an ergodic attractor, then ergodic averages of QoIs are well-defined quantities and their sensitivity can be used to characterize model sensitivity. The response theorem gives sufficient conditions such that the local forward sensitivity – the derivative with respect to a given parameter - of an ergodic average of a QoI is well-defined. We describe a method based on ergodic and response theory for computing the sensitivity and response of a given QoI with respect to a given parameter in a chaotic model with an ergodic and hyperbolic attractor. This method does not require computation of ensembles of the model with perturbed parameter values. The method is demonstrated and some of the computations are validated on the Lorenz 63 and Lorenz 96 models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Unraveling the Structural Sensitivity of Metal Catalysts in Ethylene Hydroformylation: Insights from Theory and Experiments

Here, in this study, we combined experimental and theoretical methods to investigate the structural sensitivity of metal catalysts in the ethylene hydroformylation reaction. Among Rh, Pt, Ir, Ni, Au, Ag, Pd, and Cu catalysts studied using experimental and theoretical methods, Rh showed the highest selectivity toward the C-C coupling product from CO and C 2 H 5 (i.e., C 2 H 5 CHO). The results from DFT and microkinetic simulations revealed that the activation energy barrier for C-C coupling is lowest on the Rh nanocluster, which explains the experimentally observed highest C 2 H 5 CHO selectivity on the Rh catalyst. Furthermore, DFT results demonstrated that the sites located on the flat surfaces of nanoparticles primarily promote the hydrogenation reaction, leading to the formation of undesired C 2 H 6 . In contrast, undercoordinated edge and corner sites of the nanocluster promote the C-C coupling reaction. Thus, our results illustrate that the selectivity toward C 3 oxygenates in ethylene hydroformylation reaction can be steered by tuning the size of Rh nanoparticles (the best-performing catalyst) to optimize the active (edge and corner) sites that preferentially promote the C-C coupling reaction.

58 GEOSCIENCES↗
DOE OSTI2025

A stress-sensitive precipitate nucleation model beyond classical nucleation theory

The dynamic evolution of precipitates and second phases dictates the strength and stability of most engineering alloys. By design, or as a consequence of thermo-mechanical aging, engineering metals and alloys often form precipitates of second phases when subjecting to diverse thermal and mechanical loads. Precipitation is governed by several factors, including the alloy’s composition, processing/operating temperature, and stresses — either as a result of external loads or from residual stresses. However, state-of-the-art models for precipitate nucleation (i.e., classical nucleation theory) typically lacks consistent method to capture the effects of externally applied and/or internal stresses on nucleation; thereby severely limiting the applicability of these models to complex materials systems and to representative loading scenarios. Here, in this work, we extend upon classical nucleation theory to account for the effect of stresses on precipitation kinetics and thermodynamics. This is achieved via the use of an Eshelbian micromechanics framework keeping track of (i) the stress build up resulting from second phase formation as a function of mechanical load and, (ii) the effects of dislocations on precipitate formation. This new model is applied to σ precipitate in Fe–Cr binary alloys and M 23 C 6 precipitate in 316H stainless steel (SS). Simulations demonstrate the important role of both the remotely applied loads and dislocation pile ups on precipitate nucleation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dark Matter-Induced Nuclear De-Excitation at SBND with Ab Initio Nuclear Theory

We explore the sensitivity of the Short-Baseline Near Detector (SBND) experiment to light dark matter using MeV-scale electromagnetic activity. Inelastic scattering of dark matter with argon nuclei can excite nuclear states that subsequently de-excite via the emission of MeV-scale photons, producing localized low-energy "blip" signatures in a liquid argon time projection chamber. We perform state-of-the-art ab initio nuclear calculations, including all relevant argon excited states with energies up to 18 MeV, to provide reliable predictions for these signals. After accounting for relevant backgrounds, we find that SBND can probe previously unexplored regions of parameter space for light dark matter.

Dutta, Bhaskar [Bucharest, IFIN-HH; Texas A-M] (OR↗
DOE OSTI2024

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Statistical uncertainty of fission matrix eigenvalues using perturbation theory

Eigenvalue search of high dominance ratio systems may be slow to converge. The fission matrix element is defined by its element (FM){sub ij}, which are the probability for a neutron born in cell i to create a fission in cell j for a spatial mesh of n{sub i}*n{sub j} elements. Fission matrices are used in Monte Carlo criticality simulations to enhance computing speed, but also to find higher order eigenvalues. However, few studies have been made on the link between statistical uncertainties of fission matrix elements and eigenvalues uncertainties. Thus, dominance ratio statistical uncertainties remain unknown. This paper uses a new generalized perturbation theory (GPT) method to estimate sensitivities of eigenvalues to fission matrix elements and then to calculate dominance ratio uncertainties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

High-Level Reverse Intersystem Crossing and Molecular Rigidity Improve Spin Statistics for Triplet–Triplet Annihilation Upconversion

The structural factors affecting triplet−triplet annihilation (TTA) at the molecular level are not well-understood. Here, our steady-state photoluminescence and transient absorption results demonstrate that the spin statistical factor, η, decreases from 0.60 to 0.46 and 0.14 going from 9,10-diphenylanthracene (DPA) to the 1,5-DPA and 2,6-DPA isomers, respectively, during photon upconversion with a platinum octaethylporphyrin sensitizer. Density functional theory (DFT) shows that η depends on the energetics of hot triplet states and molecular rigidity. The significantly high conical intersection energy between the S 0 and T 1 states for 9,10-DPA gives its longer triplet lifetime. Time-dependent DFT calculations show that 9,10-DPA and 1,5-DPA can undergo high-level reverse intersystem crossing from their T 2 and T 3 states, respectively, to the bright S 1 state, increasing the limit of the spin statistical factor. Both factors ultimately serve to enhance the TTA efficiency. Furthermore, this work provides insight into designing molecules for efficient light-emitting and photon upconversion applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗