Search NASA⌕ Search

SEARCH · Search NASA

Results for “extrapolation method”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Recommendations for an Applicant to Calculate Activity Data for Greenhouse Gases Estimates

In 2009, the U.S. Nuclear Regulatory Commission (NRC) directed the NRC staff to address climate change issues and consider the impacts of the emissions of carbon dioxide (CO 2 ) and other greenhouse gases (GHGs) in its environmental reviews for major licensing actions (NRC 2009b). To implement this direction from the Commission, the staff issued guidance in 2011 and updated guidance in 2014 in Attachment 1 to Interim Staff Guidance COL/ESP-ISG-026 (NRC 2011; NRC 2014). This guidance provides a simpler method than the method described in RG 4.2 Rev. 3, that an applicant can use to meet the plant parameter envelope (PPE) value from the Generic Environmental Impact Statement for Licensing of New Nuclear Reactors (NR GEIS). NRC staff estimated the 97-year lifecycle GHG emissions from a reference 1000 megawatt electrical (MWe) light-water reactor (LWR) for various activities associated with construction, operation (including uranium fuel cycle), and decommissioning of nuclear power plants and presented the results in Appendix H of the NR GEIS. Appendix H of the NR GEIS includes estimates of direct emissions from construction equipment and emergency diesel engines in a nuclear facility and indirect emissions from workforce vehicular traffic, fuel transportation and the uranium fuel cycle. The NR GEIS Section 3.3 extended the estimates in Appendix H for the installation of two 1000 MWe nuclear reactors on the same site. Scaling factors were used to extrapolate the GHG emissions of a reference 1000 MWe reactor to a two-unit nuclear reactor plant (each reactor unit generating 1000 MWe). GHG emission estimates for building, operation, decommissioning and safe storage (SAFSTOR) for a two-unit nuclear reactor plant would be based on the plant’s physical size, and therefore estimates for these source categories were assumed to be twice the value of the reference 1000 MWe reactor. However, GHG emissions from the fuel cycle (including fuel transportation) were scaled upward by a factor of 3, based on plant efficiencies greater than the 80 percent assumption in Appendix H. Table 1 below shows the PPE emissions for two 1000 MWe nuclear reactors as provided in NR GEIS. The total GHG emissions for two 1000 MWe reactors were calculated as 2,534,000 metric tons (MT) of CO 2 equivalent (CO 2 (e)) based on a 97 year GHG life cycle period. The GHG emissions lifetime of 97 years for a reference nuclear reactor includes a 7-year building phase, 40 years of operation, 10 years of active decommissioning, and 40 years of SAFSTOR operations (NRC 2024). Construction equipment and vehicular traffic from workers commute would contribute to the GHG emissions during a 7-year building phase. Uranium fuel cycle, vehicular traffic, fuel and waste transportation, and testing of standby diesel generators would contribute to GHG emissions during the 40-year operations phase. While NRC’s regulations allow up to 60 years of reactor facility decommissioning, Appendix H estimated that most of the GHGs would occur over an estimated 10-year period during which to the licensee would engage in significant demolition and earth-moving activities, as discussed in Supplement 1 to NUREG-0586 (NRC 2002). Vehicular traffic by the workforce during a 40-year SAFSTOR period would additionally contribute GHG emissions. The carbon footprint for a 40-year SAFSTOR period was separately analyzed from the decommissioning activities as provided in Table YYYY-2 of the staff issued guidance in 2011 (NRC 2011).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Learning two-phase microstructure evolution using neural operators and autoencoder architectures

Abstract Phase-field modeling is an effective but computationally expensive method for capturing the mesoscale morphological and microstructure evolution in materials. Hence, fast and generalizable surrogate models are needed to alleviate the cost of computationally taxing processes such as in optimization and design of materials. The intrinsic discontinuous nature of the physical phenomena incurred by the presence of sharp phase boundaries makes the training of the surrogate model cumbersome. We develop a framework that integrates a convolutional autoencoder architecture with a deep neural operator (DeepONet) to learn the dynamic evolution of a two-phase mixture and accelerate time-to-solution in predicting the microstructure evolution. We utilize the convolutional autoencoder to provide a compact representation of the microstructure data in a low-dimensional latent space. After DeepONet is trained in the latent space, it can be used to replace the high-fidelity phase-field numerical solver in interpolation tasks or to accelerate the numerical solver in extrapolation tasks.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An analysis of carrier dynamics in methylammonium lead triiodide perovskite solar cells using cross correlation noise spectroscopy

Using cross correlation current noise spectroscopy, we have investigated carrier dynamics in methylammonium lead triiodide solar cells. This method provides space selectivity for devices with a planar multi-layered structure, effectively amplifying current noise contributions coming from the most resistive element of the stack. In the studied solar cells, we observe near full-scale shot noise, indicating the dominance of noise generation by a single source, likely the interface between the perovskite and the spiro-organic 2,2',7,7'-tetrakis(N,N-di-p-methoxyphenyl-amine)9,9'-spirobifluorene hole-transport layer. In this paper, we argue that the strong 1/f noise term has contributions from both the perovskite layer and interfaces. It displays a non-ideal dependence on photocurrent, S $\propto$ I 1.4 (instead of usual S $\propto$ I 2 ), which is likely due to current-induced halide migration. Finally, we observe generation–recombination noise. We argue that this contribution is due to bimolecular recombination in the perovskite bulk absorption layer. Extrapolating our results, we estimate that at standard 1 sun illumination, the electron–hole recombination time is 5 microseconds.

14 SOLAR ENERGY↗
DOE OSTI2023

Data reduction considerations for the burning velocity of spherical constant volume flames of R32 (CH 2 F 2 ) with air

Here, the present work explores data reduction techniques for the measurement of the laminar burning velocities of R32(CH 2 F 2 )-air mixtures using a constant volume combustion device, in which the pressure-time history is the only measured parameter. To allow clear assessment of the accuracy of the data reduction methods, the pressure-time histories used for analysis are synthetically generated via a detailed numerical simulation employing full kinetics and with and without an optically-thin radiation model. Various data reduction models are employed, including a two-zone model and two multi-zone models, and these are compared with the results from the burning velocity obtained from the output of the numerical simulation. The data reduction schemes are shown to be accurate if the same radiation model is employed in the data reduction as was used in the flame simulation to generate the pressure trace used for post-processing. If the incorrect radiation model is employed, however, the errors can be quite large. The effects of stretch, radiation, and different data post-processing methodologies are explored and the errors quantified. Stretch is shown to be important for the early stages and the selected data range that is used for extrapolation has a significant effect on the extrapolated burning velocity. However, with an appropriate choice of data considered for extrapolation, the prediction of the unstretched burning velocity can be quite accurate.

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

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

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

SHαDE: survey description and mass–kinematics scaling relations for dwarf galaxies

ABSTRACT The Study of H α from Dwarf Emissions (SHαDE) is a high spectral resolution (R = 13 500) H α integral field survey of 69 dwarf galaxies with stellar masses 106 < M⋆ < 109 M⊙. The survey used FLAMES on the ESO Very Large Telescope. SHαDE is designed to study the kinematics and stellar populations of dwarf galaxies using consistent methods applied to massive galaxies and at matching level of detail, connecting these mass ranges in an unbiased way. In this paper, we set out the science goals of SHαDE, describe the sample properties, outline the data reduction, and analysis processes. We investigate the log M⋆–log S0.5 mass–kinematics scaling relation, which has previously shown potential for combining galaxies of all morphologies in a single scaling relation. We extend the scaling relation from massive galaxies to dwarf galaxies, demonstrating this relation is linear down to a stellar mass of M⋆ ∼ 108.6 M⊙. Below this limit, the kinematics of galaxies inside one effective radius appears to be dominated by the internal velocity dispersion limit of the H α-emitting gas, giving a bend in the log M⋆–log S0.5 relation. Replacing stellar mass with total baryonic mass using gas mass estimate reduces the severity but does not remove the linearity limit of the scaling relation. An extrapolation to estimate the galaxies’ dark matter halo masses, yields a log Mh–log S0.5 scaling relation that is free of any bend, has reduced curvature over the whole mass range, and brings galaxies of all masses and morphologies on to the virial relation.

Barat, Dilyar↗
DOE OSTI2021

Lessons Learned in Employing Data Analytics to Predict Oxidation Kinetics and Spallation Behavior of High-Temperature NiCr-Based Alloys

Machine learning (ML) can offer many advantages in predicting material properties over traditional materials development methods based solely on limited experimental investigations or physical-based simulations with the capability to reduce development cost, risk, and time. However, so far, limited efforts have been made to predict alloy oxidation kinetics and spallation behavior via ML due to the lack of consistently measured and sufficient experimental data and the inherent complexity in oxidation behavior of multicomponent high-temperature alloys. A previous study reported the ability of ML to predict oxidation kinetics of NiCr-based alloys as a function of alloy composition and operating conditions. Here, the performance of a ML model in predicting rate constants and spallation probability was evaluated in light of the roles of the data distribution of the experimental dataset (data analytics), the alloy composition, the exposure environment and the chosen oxidation approach to extracting kinetic values from the measured mass changes (but using either a simple parabolic law or a statistical cyclic oxidation model). Potential strategies to improve the predictions and enhance the extrapolative capability of the previously trained model will be discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Rate Coefficients for OH + NO (+N 2 ) in the Fall-off Regime and the Impact of Water Vapor

The termolecular, association reaction between OH and NO is a source of nitrous acid (HONO), an important atmospheric trace gas. Rate coefficients for the title reaction as recommended by evaluation panels differ substantially at the temperatures and pressures that prevail in the Earth’s boundary layer where the reaction is in the fall-off regime between low- and high-pressure limiting rate coefficients. Using pulsed laser methods for generation and detection of OH, we have reinvestigated the kinetics of the title reaction at pressures of 22–743 Torr (1 Torr = 1.333 hPa) and temperatures (273, 298, and 333 K) in pure N 2 and in N 2 –H 2 O bath gases. In situ optical absorption measurements were used to rule out any bias due to NO 2 or HONO impurities. Our rate coefficients (k 1 ) in N 2 bath gas are parametrized in terms of low-pressure (k o ) and high-pressure (k ∞ ) rate coefficients and a fall-off parameter (F C ) with k 1,0 N 2 = 7.24 × 10 –31 (T/300 K) –2.17 cm 6 molecule –2 s –1 , k 1,∞ = 3.3 × 10 –12 (T/300 K) –3 cm 3 molecule –1 s –1 , and F C = 0.53. Used with the “Troe” expression for termolecular reactions, these parameters accurately reproduce the current data in the fall-off regime and also capture literature rate coefficients at extrapolated temperatures. The presence of water vapor was found to enhance the rate coefficients of the title reaction significantly. The low-pressure limiting rate coefficient in H 2 O bath gas is a factor 5–6 larger than in N 2 , at room temperature (k 1,o H 2 O = 4.55 × 10 –30 (T/300 K) –4.85 cm 6 molecule –2 s –1 ) indicating that H 2 O is much more efficient in quenching the association complex HONO* through collisional energy transfer. Based on measurements in N 2 –H 2 O mixtures, a parametrization of k 1 including both N 2 and H 2 O as third-body quenchers was derived. Neglecting the effect of H 2 O results, e.g., in an underestimation of k 1 by >10% in the tropical boundary layer.

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

A consistent, Bayesian, approach to the cross section probability distribution in the unresolved resonance region

The cross sections of neutron-induced reactions can be divided into three energy ranges: the resolved resonance region (RRR), the unresolved resonance region (URR), and the fast region. In general, the cross sections in the URR show significant fluctuations that cannot be predicted and cannot be experimentally resolved, thus, it is commonly assumed that the cross section at a specific energy is given by a probability distribution function (PDF) over a range of values that can span several orders of magnitude. The current methodology used to describe such behavior is to construct the PDF by stochastically generating resonance ladders and numerically measuring the PDF. The resonance ladders are sampled using known resonance statistical properties and average resonance widths and spacings extrapolated from the RRR. Although this is a standard and widely used technique, it is computationally very expensive, therefore, an alternative, analytical, approach would be preferable due to the considerable speed up of the computational time in real life applications. Moreover, the current methodology does not take into account existing experimental data, such for total and capture cross sections, that are available for many nuclei. Finally, this approach was developed to be used in reactor-scale applications and it is not suited for use in single-event applications. In this work we will rethink the entire approach to the PDF construction using a Bayesian mindset. This will allow us to provide a different definition of the PDF that allows a much faster calculation of the higher-temperature PDFs and a proper combination of theoretical and experimental PDFs following the probability theory. We will also show that our definition is well suited for single-event applications and we will make an explicit connection between our method and the standard approach. We do this by showing that the central limit theorem applies and our method leads to the same PDF obtained with the standard methodology, for a large number of events per history.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes

Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.

Chemical engineering↗
DOE OSTI2022

Accelerating Multiscale Materials Modeling with Machine Learning

The focus of this project is to accelerate and transform the workflow of multiscale materials modeling by developing an integrated toolchain seamlessly combining DFT, SNAP, LAMMPS, (shown in Figure 1-1) and a machine-learning (ML) model that will more efficiently extract information from a smaller set of first-principles calculations. Our ML model enables us to accelerate first-principles data generation by interpolating existing high fidelity data, and extend the simulation scale by extrapolating high fidelity data (10 2 atoms) to the mesoscale (10 4 atoms). It encodes the underlying physics of atomic interactions on the microscopic scale by adapting a variety of ML techniques such as deep neural networks (DNNs), and graph neural networks (GNNs). We developed a new surrogate model for density functional theory using deep neural networks. The developed ML surrogate is demonstrated in a workflow to generate accurate band energies, total energies, and density of the 298K and 933K Aluminum systems. Furthermore, the models can be used to predict the quantities of interest for systems with more number of atoms than the training data set. We have demonstrated that the ML model can be used to compute the quantities of interest for systems with 100,000 Al atoms. When compared with 2000 Al system the new surrogate model is as accurate as DFT, but three orders of magnitude faster. We also explored optimal experimental design techniques to choose the training data and novel Graph Neural Networks to train on smaller data sets. These are promising methods that need to be explored in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Medial axis and local thickness computation using the Fast Sweeping Method

This report describes an efficient and robust voxel-based methodology for computing the medial axis, local thickness, and distance-to-skeleton of arbitrary three-dimensional geometries. It is assumed that the object can be represented by an exact or approximate signed distance function on a discrete grid. The gradient of such function is used to formulate a hyperbolic partial differential equation (PDE) that models the collapse of the position vector in space. By exploiting the causality property of the PDE, the Fast Sweeping Method is able to obtain the solution in a finite number of sweeps independent of the mesh resolution. The intersection of characteristic lines leads to the formation of shocks and a discrete bisector function is used to identify the medial axis. The same PDE approach is used to compute the local thickness inside the object and obtain the distance-to-skeleton field. Multiple examples are given in two and three dimensions along with a resolution study. The methodology has optimal complexity and yields subsecond computational times for geometries with over a million zones on a single core. The methodology is also capable of parallelization across shared and distributed memory architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Identifying Abandoned Well Sites Using Database Records and Aeromagnetic Surveys

Oil and natural gas are primary sources of energy in the United States. Improved drilling and extracting techniques have led to a renewed interest in historic oil and gas fields, but limited records of legacy wells make new drilling efforts more difficult, as abandoned wells may provide conduits for liquids and gases to migrate into groundwater reservoirs or the atmosphere. Well finding using aeromagnetic surveys pinpoints the location of steel-cased wells, detecting both active and abandoned wells, including buried casings lacking aboveground markers. Here, we present six aeromagnetic surveys conducted in Pennsylvania and Wyoming as case studies, comparing the magnetic points to locations known in databases. In all study sites, more magnetic points were detected than recorded in databases. Based on differences between theoretical database well counts and the actual number of wells detected in surveys, we estimated the total number of wells in Pennsylvania to be 395 000–466 000 and 181 000–182 000 in Wyoming. Extrapolating to the national level, in this work, we estimate the average number of wells in the continental United States is 6.04 ± 19.97 million wells with 1.16 ± 3.84 million of those designated as abandoned wells, within the range of previous abandoned well count estimations. Although aeromagnetic surveys are limited to detecting steel-cased wells and do not differentiate sites based on well status, this study nevertheless demonstrates the utility of aeromagnetic surveys in well finding efforts across the country and shows limitations in database records of oil and natural gas wells.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Embedded training of neural-network subgrid-scale turbulence models

We report that the weights of a deep neural-network model are optimized in conjunction with the governing flow equations to provide a model for subgrid-scale stresses in a temporally developing plane turbulent jet at Reynolds number Re 0 = 6000 . The objective function for training is first based on the instantaneous filtered velocity fields from a corresponding direct numerical simulation, and the training is by a stochastic gradient descent method, which uses the adjoint Navier-Stokes equations to provide the end-to-end sensitivities of the model weights to the velocity fields. In-sample and out-of-sample testing on multiple dual-jet configurations show that its required mesh density in each coordinate direction for prediction of mean flow, Reynolds stresses, and spectra is half that needed by the dynamic Smagorinsky model for comparable accuracy. The same neural-network model trained directly to match filtered subgrid-scale stresses, without the constraint of being embedded within the flow equations during the training, fails to provide a qualitatively correct prediction. The coupled formulation is generalized to train based only on mean-flow and Reynolds stresses, which are more readily available in experiments. The mean-flow training provides a robust model, which is important, though a somewhat less accurate prediction for the same coarse meshes, as might be anticipated due to the reduced information available for training in this case. The anticipated advantage of the formulation is that the inclusion of resolved physics in the training increases its capacity to extrapolate. This is assessed for the case of passive scalar transport, for which it outperforms established models due to improved mixing predictions.

42 ENGINEERING↗
DOE OSTI2024

Nuclear-matter saturation and symmetry energy within Δ -full chiral effective field theory

Nuclear saturation and the symmetry energy are key properties of low-energy nuclear physics that depend on fine details of the nuclear interaction. The equation of state around saturation is also an important anchor for extrapolations to higher densities and studies of neutron stars. Here we develop a unified statistical framework that uses realistic nuclear forces to link the theoretical modeling of finite nuclei and infinite nuclear matter. We construct fast and accurate emulators for nuclear-matter observables and employ an iterative history-matching approach to explore and reduce the enormous parameter domain of Δ -full chiral interactions. We perform rigorous uncertainty quantification and find that model calibration including O 16 observables gives saturation predictions that are more precise than those that only use few-body data. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Nonlinear X-ray Optics

Novel X-ray free-electron lasers (XFELs) are capable of producing radiation with unprecedented properties. Specifically, their ultrashort coherent X-ray pulses can generate extreme peak intensities where ordinary rules of light-matter interaction may no longer apply, and nonlinear processes start to become important. XFELs have made it possible to observe some of the most important coherent non-sequential nonlinearities, including X-ray-optical sum frequency generation (SFG), X-ray second harmonic generation (XSHG), two-photon absorption (TPA) and nonlinear two-photon Compton scattering (2PCS). In particular, the 2PCS experiment has led to unexpected results, namely the observation of a substantial anomalous red shift in the energy of photons generated by nonlinear two-photon X-ray Compton scattering in beryllium. The energy shift is in addition to the predicted nonlinear Compton shift and was neither expected from extrapolations from linear X-ray interactions nor from nonlinear effects at optical wavelengths. Our results suggest a novel nonlinear scattering mechanism that can only be observed at high fields at X-ray wavelengths, where the photon energy is a significant fraction of the electron rest mass (and thus a significant momentum transfer to the electron occurs during the scattering) and where the interaction is not dominated by the quiver motion of the electron in the field (ponderomotive potential) as is usually the case at longer (optical) wavelengths. The goal of this research project was the investigation of fundamental nonlinear X-ray-matter interactions. As these processes are typically weak, part of this research is the exploration of more efficient nonlinear X-ray processes by studying novel mechanisms including solid-state effects, the assistance of optical radiation and novel phase-matching techniques. These experiments require extremely large X-ray field strengths, such as produced by XFELs. Nonlinear X-ray effects can lead to important applications in many research fields. For example, it could lead to instantaneous plasma diagnostics for materials in extreme conditions, as a method for combining atomic-scale structural sensitivity with chemical specificity, or as a probe of electron dynamics in solids on the attosecond (1 as = 10 -18 s) time- and Ångstrom (10 -10 m) length scales. In addition, a detailed understanding of nonlinear X-ray matter interactions is crucial for understanding and extending the limits of structural determination at high intensity such as in serial femto-crystallography. The results could have a profound impact on advanced XFEL light sources such as the LCLS-II.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

$\overline{TKE}$ Parameterization and $\bar{v}$ Uncertainty Analysis for CGMF

Previous work was performed on tuning CGMF parameters for 235 U, 238 U, and Plutonium isotopes. Now work is being done to tune minor uranium isotopes. However, uranium isotopes like 232 U and 236 U have almost no experimental data. We are applying cross-isotope models to extrapolate and tune CGMF on isotopes that lack experimental data. There exist several internal CGMF physics quantities that affect the output of CGMF—multi-chance fission probability, excitation energy sharing, spin-cutoff factor, spin scaling, and fragment total kinetic energy to name a few. The mean fragment total kinetic energy, $\overline{TKE}$, is particularly interesting because of its strong anti-correlation with $\bar{v}$. We are most interested in the mean fragment total kinetic energy before neutron emissions. $\overline{TKE}$ is assumed to be pre-neutron emission unless otherwise stated. Currently in CGMF, the $\overline{TKE}$ model for 233,234,235,238 U are tuned independently to reproduce ν for the associated isotopes. In this report, we will tune a cross-isotope $\overline{TKE}$ model to experimental $\overline{TKE}$ data for 232,233,234,235,236,238 U. Because of the unreliable and sparse nature of $\overline{TKE}$ experimental data, future work will use more reliable experimental $\bar{v}$ data to infer the $\overline{TKE}$ model (and likely other internal CGMF parameters) for uranium isotopes. Such work has been performed previously using a sensitivity analysis and Kalman filter methods.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

From chiral effective field theory to perturbative QCD: A Bayesian model mixing approach to symmetric nuclear matter

Constraining the equation of state (EOS) of strongly interacting, dense matter is the focus of intense experimental, observational, and theoretical effort. Chiral effective field theory (𝜒⁢EFT ) can describe the EOS between the typical densities of nuclei and those in the outer cores of neutron stars, while perturbative QCD (pQCD) can be applied to properties of deconfined quark matter, both with quantified theoretical uncertainties. However, describing the full range of densities in between with a single EOS that has well-quantified uncertainties is a challenging problem. Bayesian multimodel inference from 𝜒⁢EFT and pQCD can help bridge the gap between the two theories. In this work, we introduce a correlated Bayesian model mixing framework that uses a Gaussian process (GP) to assimilate different information into a single QCD EOS for symmetric nuclear matter. The present implementation uses a stationary GP to infer this mixed EOS solely from the EOSs of 𝜒⁢EFT and pQCD while accounting for the truncation errors of each theory. The GP is trained on the pressure as a function of number density in the low- and high-density regions where 𝜒⁢EFT and pQCD are, respectively, valid. We impose priors on the GP kernel hyperparameters to suppress unphysical correlations between these regimes. This, together with the assumption of stationarity, results in smooth 𝜒⁢EFT-to-pQCD curves for both the pressure and the speed of sound. We show that using uncorrelated mixing requires uncontrolled extrapolation of at least one of 𝜒⁢EFT or pQCD into regions where the perturbative series breaks down and leads to an acausal EOS. Here, we also discuss extensions of this framework to nonstationary and less differentiable GP kernels, its future application to neutron-star matter, and the incorporation of additional constraints from nuclear theory, experiment, and multimessenger astronomy.

Bayesian methods↗