Search NASA⌕ Search

SEARCH · Search NASA

Results for “Formal analysis”

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

Mean Time to Failure Analysis for the Los Alamos National Laboratory Weapons Engineering Tritium Facility Tritium Containment Vessels

The Weapons Engineering Tritium Facility (WETF) is a Hazard Category 2 nonreactor nuclear facility located in Technical Area 16 (TA-16) of Los Alamos National Laboratory. Tritium related research, development, and processing in support of stockpile stewardship program activities is conducted at WETF. Within WETF, various devices can be used to contain tritium and certain devices are designated as Tritium Containment Vessels (TCVs). Several different device types are identified as TCVs within WETF such as LP-50s and AL-M1s. Historically, a threshold tritium mass value has been used as the mechanism for designating a device as a TCV. When the tritium content of the device was below the tritium mass threshold it was not considered to be a TCV in the formal context. When the tritium content was above the administrative threshold, the device was considered to be a TCV. An upper tritium mass limit was also established for TCVs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Selenium interaction with iron minerals: Quantitative comparison of sorption and coprecipitation impacts on mobility

Given the significance of selenium (Se) as a micronutrient, the radioactive nature of some of its isotopes, and its affinity to iron (Fe) minerals, extensive research has been conducted on the sorption mechanisms between Se and these minerals. Here, in this study, we employ sorption data sourced from the L-SCIE database and coprecipitation data from available literature to achieve the following objectives: i) establish coherence between adsorption and coprecipitation processes, ii) quantitatively evaluate the importance of these processes in nuclear waste repository science, and iii) propose a forward-looking approach for integrating coprecipitation into reactive transport models. Our findings indicate that a correlation between Se adsorption and coprecipitation can be established using the λ formalism. The comparable log(λ Se(IV) /λ Se(VI) ) ratios derived from adsorption and coprecipitation experiments suggest that these processes can be quantitatively compared and evaluated using our numerical approach. Across all iron oxide phases examined, coprecipitation leads to significantly greater immobilization of Se compared to adsorption. Specifically, for hydrous ferric oxide, hematite, and goethite, coprecipitation is predicted to result in 100–1000 times more Se immobilization compared to adsorption, irrespective of the Se oxidation state (Se(IV) or Se(VI)); notably stronger immobilization potential via coprecipitation was observed for magnetite. The modeling approach and quantitative analysis presented herein clearly highlight the importance of including coprecipitation processes when simulating Se (and other elements) transport, particularly under conditions where mineral compositions are transient or evolving with time. Neglecting coprecipitation in models is likely to lead to significant overestimates of migration.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Benders Decomposition Using Graph Modeling and Multi-Parametric Programming

Benders decomposition is a widely used method for solving large and structured optimization problems, but its performance is affected by the repeated solution of subproblems. We propose a flexible and modular algorithmic framework for accelerating Benders decomposition. Specifically, we express the problem structure by using a graph-theoretic modeling abstraction in which nodes represent optimization subproblems and edges represent connectivity between subproblems. A key innovation of our approach is that we embed multiparametric programming (mp) surrogates for node subproblems, which maps the exact analytical map of the subproblem solution space. The use of mp surrogates allows us to replace subproblem solves with fast look-ups and function evaluations for primal and dual variables during the iterative Benders process. We formally show the equivalence between classical Benders cuts and those derived from the mp solution. We implement our framework in the open-source PlasmoBenders.jl software package. To demonstrate the capabilities of the proposed framework, we apply it to a two-stage stochastic programming problem, which aims to make optimal capacity expansion decisions under market uncertainty. We evaluate both single-cut and multicut variants of Benders decomposition and show that the use of mp surrogates achieves substantial speedups in subproblem solve time, while preserving the convergence guarantees of Benders decomposition. We highlight advantages in solution analysis and interpretability that is enabled by mp critical region tracking; specifically, we show that these reveal how decisions evolve geometrically across the Benders search. Our results aim to demonstrate that combining surrogate modeling with graph modeling offers a promising and extensible foundation for structure-exploiting decomposition. In addition, by decomposing the problem into more tractable subproblems, the proposed approach also aims to overcome scalability issues of mp. Finally, the use of mp surrogates provides a unifying and modular optimization framework that enables the representation of heterogeneous node subproblems as modeling objects with a homogeneous structure.

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

Competition between roughness and strength for scale-dependent surfaces

Rocks famously have scale-dependent strength, yet the actual dependence is notoriously hard to measure or incorporate into any theoretical framework. Natural rough surfaces present an opportunity to solve the problem. Surfaces sliding in shear evolve as protrusions collide. These asperities can deform or break, thus creating a new surface shape. In particular, natural surfaces have roughness at all scales as well as scale-dependent strength. Based on a scaling analysis, we have previously suggested that the scale-dependent aspect ratio of steady-state surfaces should be proportional to the scale-dependent shear strain at yield. If true, scale-dependent strength could easily be inferred from natural surfaces. Thus, moving beyond the scaling argument to a rigorous treatment of scale-dependent strength for multiscale rough surfaces in shear is important. However, analytic frameworks for analyzing multiscale problems are challenging, as conventional continuum mechanics typically involves a single value for a material property across scales. Here, in this work, we build on the formalism of Persson (2001) that presents a method to compute contact area for rough surfaces with a prescribed topographic spectrum using a stochastic differential equation. The Persson formalism allows for plastic yield under normal loading of otherwise elastic materials and leaves open the possibility of scale-dependent yield stress. In this study, we pursue this route to develop a theory and numerical results for the yielding of a rough, elastoplastic surface with scale-dependent yield stress. Here, we examine surfaces for which the power spectrum of the topography 𝐶 and yield stress 𝑌 follow power laws as a function of scale 𝜆, such that 𝐶∼𝜆 −𝑚 and 𝑌∼𝜆 −𝑛 , respectively. In this formal treatment of the problem, we focus on surfaces in contact and the resulting yield and do not impose shear. Numerical solutions show that the deviation from the elastic scaling solution is bounded as expected by the prior 1D heuristic scaling argument that anticipates the Hurst exponent as 1−𝑛. We also show that the plasticity is expected to erode the contacts if 𝑚 is lower than 𝑛−3, which corresponds to a Hurst exponent lower than 1−𝑛/2. This result is rigorously sound for 2D, i.e., realistic surfaces, and quantitatively different than the prior scaling argument. The theory now permits a correspondingly quantitative approach to interpreting natural surfaces.

elasticity↗
DOE OSTI2024

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

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

Laboratory Constraints on the Neutron-Spin Coupling of feV-Scale Axions

Ultralight axionlike particles can contribute to the dark matter near the Sun, leading to a distinct, stochastic signature in terrestrial experiments. We search for such particles through their neutron-spin coupling by reanalyzing approximately 40 days of data from a K- 3 He comagnetometer with a new frequency-domain likelihood-based formalism that properly accounts for stochastic effects over all axion coherence times relative to the experimental time span. Assuming that axions make up all of the dark matter in the Sun’s vicinity, we find a median 95% upper limit on the neutron-spin coupling of 2.4 × 10 −10 GeV −1 for most axion masses from 0.4 to 4 feV, which is about 5 orders of magnitude more stringent than previous laboratory bounds in that mass range. Although several peaks in the experiment’s magnetic power spectrum suggest the rejection of a white-noise null hypothesis, further analysis of their line shapes yields no positive evidence for a dark-matter axion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Volume extrapolation via eigenvector continuation

Here, we develop an extension of eigenvector continuation (EC) that makes it possible to extrapolate simulations of quantum systems in finite periodic boxes across large ranges of box sizes. The formal justification for this approach, which we call finite-volume eigenvector continuation (FVEC), is provided by matching periodic functions at different box sizes. As concrete FVEC implementation we use a discrete variable representation based on plane-wave states and present several applications calculated within this framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

DOE Deep Energy Retrofit Cost Survey

A survey was conducted by the Lawrence Berkeley National Laboratory on deep energy retrofit (DER) market drivers, opportunities, and challenges. The survey was part a research study sponsored by the U.S. Department of Energy to gather information on the costs of DER from home performance contractors and stakeholders. Cost data was gathered from DER projects that use a comprehensive, whole-home approach to drastically reduce energy use and improve performance. DER projects often aim at reducing energy use by 50% or more. In addition, these projects can improve home comfort and potentially benefiting occupant health. Yet, market adoption of DER has been limited. Major limiting factors include complex projects, high costs, perceived risks, extensive disruption, and unfamiliar work scopes to some contractors. In order to better understand what motivates and deters DER projects in today’s market, a survey was conducted to gather this information, and to learn about promising approaches and technologies from the industry perspective. Past surveys on homeowners and home energy performance professionals have studied the motivations and barriers of energy efficiency retrofits. Two surveys of home energy performance professionals were conducted in recent years. The Resources for the Future (RFF) Home Energy Audit and Retrofit Survey was conducted in 2011 by recruiting energy auditors and retrofit installers through members of Efficiency First and Building Performance Institute (BPI) accredited contractors. The survey asked about the business and services that respondents provide, how often homeowners follow their recommendations to retrofit their homes, and the respondents’ opinion on barriers faced by the industry. The survey found that not enough homeowners know about energy audits, but more importantly, it is the high cost of retrofits compared to low energy prices that is responsible for few energy audits and retrofits being completed. The RESNET Deep Retrofit Industry Stakeholder Survey was conducted shortly after launching of the EnergySmart Home Performance Team program. The EnergySmart Team program involves a formal agreement among allied contractors who are engaged in high performance retrofits. Using this allied team approach, teams can pool their expertise and provide each other with customer referrals. The survey asked EnergySmart Team members and outside stakeholders on questions about market and technical barriers in performing home energy retrofits. Survey respondents identified lack of consumer awareness and lack of affordable financing for consumers as the leading market barriers to home energy retrofits. In their written comments, many survey respondents also echoed that the high costs of retrofit compared to low energy prices is a market barrier. Respondents found “certain housing characteristics that prevent effective retrofit” and “energy analysis software inaccuracy or limitations” are the two leading technical barriers. Their choices for technical barriers were reflective of the energy rater/auditor role played by the majority (78%) of the survey respondents. In comparison to past surveys, this work aimed to gather inputs from a broader segment of the home performance industry to identify the opportunities and barriers faced by DERs from all perspectives. The survey asks for project costs to help breakdown the high costs of DERs. This survey is also motivated by a need to better understand the role of DERs in reducing energy use by the residential sectors and meeting climate goals. The survey is designed towards obtaining more substantive inputs from survey respondents by encouraging written comments, rather than setting the goal to reach a large number of respondents. We took this approach because DER is currently still a niche market, so it is more valuable to gather in-depth inputs from individuals who are performing this work rather than getting to the masses.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Epitaxial growth of hexagonal BaNi O 3–δ thin films on SrTi O 3 (111) substrates

Transition metal oxides containing nickel species in oxidation states higher than 3+ often exhibit high catalytic activity, which makes them promising for applications as advanced electrocatalysts for water splitting and fuel cells. Here, we examine the structure and properties of BaNiO 3 (BNO) thin films, containing formally Ni 4+ , grown on SrTiO 3 (111) substrates using oxygen-plasma–assisted molecular beam epitaxy. X-ray diffraction and scanning transmission electron microscopy measurements reveal that BNO films have a hexagonal structure with $c$- and $a$-axes of mixed textures showing epitaxial relationships BNO (0001) ∥ SrTiO 3 (111) and BNO (10$\overline{1}$0) ∥ SrTi O 3 (111), respectively. The formation of the $a$-axis texture is dominant due to the smaller lattice mismatch with the substrate. Here, density functional theory calculations confirm that the hexagonal BNO film with the $a$-axis texture is energetically more favorable than the competing $c$-axis texture. Detailed spectroscopy data analysis indicates that hexagonal BNO films contain mixtures of Ni 2+ , Ni 3+ , and Ni 4+ species, with Ni 4+ being dominant. Our study provides insights into stabilizing Ni 4+ in complex oxides, which is important for further exploration of the potential of materials containing Ni 4+ .

36 MATERIALS SCIENCE↗
DOE OSTI2020

The hidden geometry of particle collisions

We establish that many fundamental concepts and techniques in quantum field theory and collider physics can be naturally understood and unified through a simple new geometric language. The idea is to equip the space of collider events with a metric, from which other geometric objects can be rigorously defined. Our analysis is based on the energy mover’s distance, which quantifies the “work” required to rearrange one event into another. This metric, which operates purely at the level of observable energy flow information, allows for a clarified definition of infrared and collinear safety and related concepts. A number of well-known collider observables can be exactly cast as the minimum distance between an event and various manifolds in this space. Jet definitions, such as exclusive cone and sequential recombination algorithms, can be directly derived by finding the closest few-particle approximation to the event. Several area- and constituent-based pileup mitigation strategies are naturally expressed in this formalism as well. Finally, we lift our reasoning to develop a precise distance between theories, which are treated as collections of events weighted by cross sections. In all of these various cases, a better understanding of existing methods in our geometric language suggests interesting new ideas and generalizations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Impact of impurities on peeling–ballooning modes and turbulence in tokamak plasmas

This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.

BOUT++ simulation↗
DOE OSTI2024

Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics

Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 10 1 –10 4 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction–Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. Here, the model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).

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

Assessment of Survey Results from Advanced Reactor Industry Domain Experts on Nonlinear Soil-Structure Interaction Analysis Software Verification and Validation

The seismic load case has a significant impact on the design and cost of nuclear power plants. Given the need to substantially reduce cost, advanced reactor designers are looking to leverage numerical tools that enable seismic analysis of an integrated vessel/support/structure/soil system. Modern nonlinear analysis tools provide a solution, capturing dynamic coupling between components (including soil-structure interaction (SSI)) concurrent with nonlinear behavior in one or more parts of the system. Although these methods have a rich history of technical development and implementation, software quality assurance (SQA) following nuclear industry standards remains a significant burden for those wishing to adopt such methods for advanced reactor design and licensing. To reduce the SQA burden, this project will develop the guidance for SQA verification and validation (V&V) of coupled nonlinear SSI analysis tools, to include a test matrix of software features and test problems to support commercial grade dedication (CGD). The guidance is intended to be technology-neutral: supporting designers of all advanced reactors. To maximize the value of the project to reactor designers, the project team performed focused surveys of and interviews with advanced reactor designers. Additionally, the project surveyed members of the broader industry involved in reactor design and licensing. The aggregated industry feedback consists of written survey responses, live polling responses, focused interviews, and informal feedback provided after outreach presentations, collectively referred to as “survey results”. Herein, the survey results are reported, reviewed, and assessed to inform follow-on project activities. The survey results confirmed the familiarity of the industry with the coupled nonlinear SSI analysis. They also affirm the project need based upon the expressed intent to (1) incorporate nonlinear features and (2) pursue the coupled nonlinear SSI analysis to reduce seismic demands and construction costs. All reactor designers indicated their intent to incorporate two or more nonlinear features, and all reactor designer respondents opined that the coupled nonlinear SSI analysis would allow for the optimization of the reactor design and construction. However, the programmatic challenges, whether real or perceived, present a significant barrier to reactor designers. The barrier most commonly identified by the reactor designer population was regulatory risk, with 70% citing this as a reason not to pursue the approach. The perceived regulatory risk identified by the reactor designers underscores the importance of regulator engagement and dialogue in this project. Additionally, half of the reactor designers identified cost and lack of guidance as a deterrent. The survey results also provide insights to tailor specific aspects of the guidance document and test matrix. The project plans to prepare both a formal referenceable guidance document and a collaborative, web-based test matrix and problem set. The project will focus on more complete test problem definitions for the prioritized nonlinear features over shallower problem descriptions for a larger set of features. The nonlinearities prioritized as high based upon survey feedback include fluid-structure interaction and seismic isolation and energy dissipation devices. The nonlinearities prioritized as intermediate include nonlinear geomaterials, nonlinear concrete, nonlinear steel, and interface/contact nonlinearity. Deep embedment and the associated nonlinear phenomena are assigned the lowest priority based upon survey feedback.

42 ENGINEERING↗
DOE OSTI2021

Linear polarization–direction correlations in $\gamma $-ray scattering experiments

Scattering measurements with incident linearly polarized γ rays provide information on spins, parities, and γ-ray multipolarity mixing coefficients, and, therefore, on the nuclear matrix elements involved in the transitions. Here, we present the general formalism for analyzing the observed angular correlations. The expressions are used to compute three-dimensional radiation patterns, which are important tools for optimizing experimental setups. Frequently, γ-ray transitions can proceed via two multipolarities that mix coherently. In such cases, the relative phases of the nuclear matrix elements are important when comparing results from different measurements. We discuss different phase conventions that have been used in the literature and present their relationships. Finally, we propose a basic experimental geometry consisting of detectors located at four different spatial locations. For this geometry, we present the measured anisotropies of the emitted γ rays in graphical format as an aid in the data analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Dynamical coupled-channel models for hadron dynamics

Dynamical coupled-channel (DCC) approaches parametrize the interactions and dynamics of two and more hadrons and their response to different electroweak probes. The inclusion of unitarity, three-body channels, and other properties from scattering theory allows for a reliable extraction of resonance spectra and their properties from data. Here, we review the formalism and application of the ANL-Osaka, the Juelich-Bonn-Washington, and other DCC approaches in the context of light baryon resonances from meson, (virtual) photon, and neutrino-induced reactions, as well as production reactions, strange baryons, light mesons, heavy meson systems, exotics, and baryon-baryon interactions. Finally, we also provide a connection of the formalism to study finite-volume spectra obtained in Lattice QCD, and review applications involving modern statistical and machine learning tools.

Amplitude analysis↗
DOE OSTI2021

A Bayesian analysis of classical shadows

Abstract The method of classical shadows proposed by Huang, Kueng, and Preskill heralds remarkable opportunities for quantum estimation with limited measurements. Yet its relationship to established quantum tomographic approaches, particularly those based on likelihood models, remains unclear. In this article, we investigate classical shadows through the lens of Bayesian mean estimation (BME). In direct tests on numerical data, BME is found to attain significantly lower error on average, but classical shadows prove remarkably more accurate in specific situations—such as high-fidelity ground truth states—which are improbable in a fully uniform Hilbert space. We then introduce an observable-oriented pseudo-likelihood that successfully emulates the dimension-independence and state-specific optimality of classical shadows, but within a Bayesian framework that ensures only physical states. Our research reveals how classical shadows effect important departures from conventional thinking in quantum state estimation, as well as the utility of Bayesian methods for uncovering and formalizing statistical assumptions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Concepts Relevant for the Kinetic Analysis of Reversible Reaction Systems

The net rate of a reversible chemical reaction is the difference between unidirectional rates of traversal along forward and reverse reaction paths. In a multistep reaction sequence, the forward and reverse trajectories, in general, are not the microscopic reverse of one another; rather, each unidirectional route is comprised of distinct rate-controlling steps, intermediates, and transition states. Consequently, traditional descriptors of rate (e.g., reaction orders) do not reflect intrinsic kinetic information but instead conflate unidirectional contributions determined by (i) the microscopic occurrence of forward/reverse reactions (i.e., unidirectional kinetics) and (ii) the reversibility of reaction (i.e., nonequilibrium thermodynamics). This review aims to provide a comprehensive resource of analytical and conceptual tools which deconvolute the contributions of reaction kinetics and thermodynamics to disambiguate unidirectional reaction trajectories and precisely identify rate- and reversibility-controlling molecular species and steps in reversible reaction systems. The extrication of mechanistic and kinetic information from bidirectional reactions is accomplished through equation-based formalisms (e.g., De Donder relations) grounded in principles of thermodynamics and interpreted in the context of theories of chemical kinetics developed in the past 25 years. Here, the aggregate of mathematical formalisms detailed herein is general to thermochemical and electrochemical reactions and encapsulates a diverse body of scientific literature encompassing chemical physics, thermodynamics, chemical kinetics, catalysis, and kinetic modeling.

42 ENGINEERING↗
DOE OSTI2023

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗