Search NASA⌕ Search

SEARCH · Search NASA

Results for “structured Gaussian process”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Model-independent constraints on Ωm and H ( z ) from the link between geometry and growth

ABSTRACT We constrain the expansion history of the Universe and the cosmological matter density fraction in a model-independent way by exclusively making use of the relationship between background and perturbations under a minimal set of assumptions. We do so by employing a Gaussian process to model the expansion history of the Universe from present time to the recombination era. The expansion history and the cosmological matter density are then constrained using recent measurements from cosmic chronometers, Type-Ia supernovae, baryon acoustic oscillations, and redshift-space distortion data. Our results show that the evolution in the reconstructed expansion history is compatible with the Planck 2018 prediction at all redshifts. The current data considered in this study can constrain a Gaussian process on H(z) to an average $9.4 {{\ \rm per\ cent}}$ precision across redshift. We find Ωm = 0.224 ± 0.066, lower but statistically compatible with the Planck 2018 cosmology. Finally, the combination of future DESI measurements with the CMB measurement considered in this work holds the promise of $8 {{\ \rm per\ cent}}$ average constraints on a model-independent expansion history as well as a five-fold tighter Ωm constraint using the methodology developed in this work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Map-level baryonification: unified treatment of weak lensing two-point and higher-order statistics

Precision cosmology benefits from extracting maximal information from cosmic structures, motivating the use of higher-order statistics (HOS) at small spatial scales. However, predicting how baryonic processes modify matter statistics at these scales has been challenging. The baryonic correction model (BCM) addresses this by modifying dark-matter-only simulations to mimic baryonic effects, providing a flexible, simulation-based framework for predicting both two-point and HOS. We show that a 3-parameter version of the BCM can jointly fit weak lensing maps' two-point statistics, wavelet phase harmonics coefficients, scattering coefficients, and the third and fourth moments to within 2% accuracy across all scales ℓ < 2000 and tomographic bins for a DES-Y3-like redshift distribution ( z ≲ 2), using the FLAMINGO simulations. These results demonstrate the viability of BCM-assisted, simulation-based weak lensing inference of two-point and HOS, paving the way for robust cosmological constraints that fully exploit non-Gaussian information on small spatial scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Planck 2018 results

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E -mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and optimal modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following final results: f NL local = -0.9 ± 5.1; f NL equil = -26 ± 47; and f NL ortho = -38 ± 24 (68% CL, statistical). These results include low-multipole (4 ≤ ℓ < 40) polarization data that are not included in our previous analysis. The results also pass an extensive battery of tests (with additional tests regarding foreground residuals compared to 2015), and they are stable with respect to our 2015 measurements (with small fluctuations, at the level of a fraction of a standard deviation, which is consistent with changes in data processing). Polarization-only bispectra display a significant improvement in robustness; they can now be used independently to set primordial NG constraints with a sensitivity comparable to WMAP temperature-based results and they give excellent agreement. In addition to the analysis of the standard local, equilateral, and orthogonal bispectrum shapes, we consider a large number of additional cases, such as scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is, however, detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 σ . Beyond estimates of individual shape amplitudes, we also present model-independent reconstructions and analyses of the Planck CMB bispectrum. Our final constraint on the local primordial trispectrum shape is g NL local = (-5.8 ± 6.5) × 10 4 (68% CL, statistical), while constraints for other trispectrum shapes are also determined. Exploiting the tight limits on various bispectrum and trispectrum shapes, we constrain the parameter space of different early-Universe scenarios that generate primordial NG, including general single-field models of inflation, multi-field models (e.g. curvaton models), models of inflation with axion fields producing parity-violation bispectra in the tensor sector, and inflationary models involving vector-like fields with directionally-dependent bispectra. Our results provide a high-precision test for structure-formation scenarios, showing complete agreement with the basic picture of the ΛCDM cosmology regarding the statistics of the initial conditions, with cosmic structures arising from adiabatic, passive, Gaussian, and primordial seed perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Spatiotemporally dynamic implicit large eddy simulation using machine learning classifiers

In this article, we utilize machine learning to dynamically determine if a point on the computational grid requires implicit numerical dissipation for large eddy simulation (LES). The decision making process is learnt through a priori training on quantities derived from direct numerical simulation (DNS) data. In particular, we compute eddy-viscosities obtained through the coarse-graining of DNS quantities and utilize their projection onto a Gaussian distribution to categorize areas that may require dissipation. If our learning determines that closure is necessary, an upwinded scheme is utilized for computing the non-linear Jacobian. In contrast, if it is determined that closure is unnecessary, a symmetric and second-order accurate energy and enstrophy preserving Arakawa scheme is utilized instead. This results in a closure framework that precludes the specification of any model-form for the small scale contributions of turbulence but deploys an appropriate numerical dissipation from explicit closure driven hypotheses. Here, this methodology is deployed for the Kraichnan turbulence test-case and assessed through various statistical quantities such as angle-averaged kinetic energy spectra and vorticity structure functions. Our framework thus establishes a link between the use of explicit LES ideologies for closure and numerical dissipation-based modeling of turbulence leading to improved statistical fidelity of a posteriori simulations.

42 ENGINEERING↗
DOE OSTI2021

Graphical Gaussian Process Regression Model for Aqueous Solvation Free Energy Prediction of Organic Molecules in Redox Flow Battery

The solvation free energy of organic molecules is a critical parameter in determining emergent properties such as solubility, liquid-phase equilibrium constants, and pKa and redox potentials in an organic redox flow battery. In this work, we present a machine learning (ML) model that can learn and predict the aqueous solvation free energy of an organic molecule using Gaussian process regression method based on a new molecular graph kernel. To investigate the performance of the ML model on electrostatic interaction, the nonpolar interaction contribution of solvent and the conformational entropy of solute in solvation free energy, three data sets with implicit or explicit water solvent models, and contribution of conformational entropy of solute are tested. We demonstrate that our ML model can predict the solvation free energy of molecules at chemical accuracy with a mean absolute error of less than 1 kcal/mol for subsets of the QM9 dataset and the Freesolv database. To solve the general data scarcity problem for a graph-based ML model, we propose a dimension reduction algorithm based on the distance between molecular graphs, which can be used to examine the diversity of the molecular data set. It provides a promising way to build a minimum training set to improve prediction for certain test sets where the space of molecular structures is predetermined.

25 ENERGY STORAGE↗
DOE OSTI2020

A Bayesian framework for adsorption energy prediction on bimetallic alloy catalysts

Abstract For high-throughput screening of materials for heterogeneous catalysis, scaling relations provides an efficient scheme to estimate the chemisorption energies of hydrogenated species. However, conditioning on a single descriptor ignores the model uncertainty and leads to suboptimal prediction of the chemisorption energy. In this article, we extend the single descriptor linear scaling relation to a multi-descriptor linear regression models to leverage the correlation between adsorption energy of any two pair of adsorbates. With a large dataset, we use Bayesian Information Criteria (BIC) as the model evidence to select the best linear regression model. Furthermore, Gaussian Process Regression (GPR) based on the meaningful convolution of physical properties of the metal-adsorbate complex can be used to predict the baseline residual of the selected model. This integrated Bayesian model selection and Gaussian process regression, dubbed as residual learning, can achieve performance comparable to standard DFT error (0.1 eV) for most adsorbate system. For sparse and small datasets, we propose an ad hoc Bayesian Model Averaging (BMA) approach to make a robust prediction. With this Bayesian framework, we significantly reduce the model uncertainty and improve the prediction accuracy. The possibilities of the framework for high-throughput catalytic materials exploration in a realistic setting is illustrated using large and small sets of both dense and sparse simulated dataset generated from a public database of bimetallic alloys available in Catalysis-Hub.org.

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

The Mira-Titan Universe. III. Emulation of the Halo Mass Function

We construct an emulator for the halo mass function over group and cluster mass scales for a range of cosmologies, including the effects of dynamical dark energy and massive neutrinos. The emulator is based on the recently completed Mira-Titan Universe suite of cosmological N-body simulations. The main set of simulations spans 111 cosmological models with 2.1 Gpc boxes. We extract halo catalogs in the redshift range z = [0.0, 2.0] and for masses M-200c >= 10(13)M(circle dot)/h. The emulator covers an eight-dimensional hypercube spanned by {Omega(m)h(2), Omega(b)h(2), Omega(nu)h(2), sigma(8), h, n(s), w(0), w(a)}; spatial flatness is assumed. We obtain smooth halo mass functions by fitting piecewise second-order polynomials to the halo catalogs and employ Gaussian process regression to construct the emulator while keeping track of the statistical noise in the input halo catalogs and uncertainties in the regression process. For redshifts z less than or similar to 1, the typical emulator precision is better than 2% for 10(13)-10(14)M(circle dot)/h and <10<^> M similar or equal to 101(circle dot)(54)/h. For comparison, fitting functions using the traditional universal form for the halo mass function can be biased at up to 30% at M similar or equal to 10(144)M(circle dot)/h for z = 0. Our emulator is publicly available at https://github.com/SebastianBocquet/MiraTitanHMFemulator.

cosmology: theory↗
DOE OSTI2025

Bayesian Calibration of Nuclear Graphite Property Models Accounting for Model Inadequacy and Impacts on Component Performance

Nuclear-grade structural graphite is extensively utilized in the core designs of various advanced nuclear reactors. In the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, including nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity requires accurate models of graphite's thermomechanical response. This report documents the Bayesian calibration of thermomechanical properties for nuclear-grade graphite and their application to graphite component modeling and simulation using the Grizzly code. As part of this work, uncertainty-quantified models were developed for the elastic modulus, coefficient of thermal expansion, irradiation-induced dimensional change, and irradiation-induced creep for graphite grades IG-110, NBG-18, NBG-17, PCEA, and 2114. Using a hierarchical Bayesian approach, multiple experimental data sources were combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due to inadequacies in the model and the inherent spread in the experimental data were quantified for three different models. These uncertainty-quantified models, with a model-form correction, were subsequently applied to a coupled-physics simulation of representative graphite components, revealing that the uncertainties have a large impact on the components' deformation.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Measurement of Charged Current Inclusive $\bar{V}_\mu$ Interactions and 2p2h Contribution Using the NOvA Near Detector

More precise measurements of the scattering of antineutrinos on nuclei in the few GeV region of incident energy, $E_{\nu}$, can reveal new aspects of $v/\bar{v}$ interactions and of nuclear structure, while providing information crucial to continued progress by neutrino oscillation experiments. This Thesis reports an investigation of charged current $\bar{v}_\mu$ interactions in the predominantly hydrocarbon medium of the NOvA Near Detector. The data are obtained with the detector exposed to a narrow-band $\bar{v}$ µ beam with spectral peak at E ν = 2.0 GeV, for a total accumulated exposure of 12.5×10$^{20}$ protons-on-target. The analysis expresses results in terms of three-momentum transfer,|$\vec{q}$| and available hadronic energy, $E$ avail . These variables approximate (q 0 , |$\vec{q}$|), the variables used in theoretical treatments of lepton scattering on multi-nucleon systems such as in 2-particle-2-hole (2p2h) processes. This Thesis delivers two core measurements: Firstly, the double-dierential cross section $d^2σ/d|\vec{q} |dE_{avail}$ for inclusive CC $\bar{v}_\mu$ interactions in the NOvA detector medium is obtained. The cross section has a roughly gaussian shape in its projection onto |$\vec{q}$|, while being largely concentrated at values of $E$ avail below 100 MeV. The single diferential cross sections dσ/d|$\vec{q}$| and dσ/dE avail are also obtained, and are compared to predictions based on the GENIE v3.0.6 neutrino generator augmented with four different implementations of 2p2h reactions. Secondly, the contribution of 2p2h reactions to the CC inclusive scattering sample is determined. Its presence is characterized in terms of double-differential and single-differential cross sections. The latter 2p2h distributions are compared to predictions from two experimental data-tune models and from two theory-based models, namely the SuSAv2 model and the Valencia model. Both of the theoretical models underpredict the measured cross-section strength and predict a shape for dσ/d|$\vec{q}$| that is flatter than the measured shape. Thus, this Thesis provides new information about CC antineutrino-nucleus scattering in the vicinity of Eν ~ 2 GeV, and it highlights improvements that are needed in theoretical treatments of 2p2h reactions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Adaptive Discovery and Mixed-Variable Optimization of Next Generation Synthesizable Microelectronic Materials

Design of new microelectronic materials is characterized by several challenges such as high-dimensionality of the atomic structure-composition variable space, formidable cost of directly using high-fidelity simulations for design optimization, dispersity in literature-reported similar materials and synthesis methods, complex physical mechanisms, and mixed qualitative and quantitative design variables that lead to a disjointed design space. Even though machine learning (ML) techniques have been employed to expedite materials innovation, existing methods treat ML and design optimization as two separate processes, failing to resolve the fundamental challenges associated with high dimensionality and mixed-variable complexity. We have developed a ML enhanced mixed-variable material design optimization framework to efficiently extract useful information from existing data in literature and physics-based simulations to guide the autonomous search for optimal materials. Our proposed framework is composed of four computational modules: (1) a natural language processing (NLP) based virtual screening module, (2) classification based concept exploration module, (3) a density functional theory (DFT)-based high-fidelity evaluation model, and (4) a novel latent-variable Gaussian process (LVGP) ML model for mixed-variable problems with uncertainty quantification, which seamlessly integrates with Bayesian Optimization (BO) and achieves superb efficiency through embedded physics-based dimension reduction. Our approach is demonstrated and validated using the testbed of functional materials exhibiting metal-insulation transitions (MITs), with the targeted reversible resistivity changes (∼10^5) near room temperature. At the end of the 30-month project, we have developed a series of new ML techniques using NLP, conditional variational autoencoders, active learning, latent-variable Gaussian processes, integrated with Bayesian optimization. Our project has resulted in new predicted MITs compounds and improved understanding of MITs microscopic mechanisms, which in turn will revolutionize microelectronics science to provide energy-saving solutions. Our research has improved both creativity and efficiency in transforming rare-event discoveries of new functional materials to persistent innovations. In addition to open-sourcing the online MIT database and the classification model, the LVGP open source code has been downloaded more than 15,000 times within two years. More than 40 MIT compounds have been identified and many have been pursued experimentally via collaborators. The research results are published in close to 20 collaborative papers in high-impact journals, such as Chem. Mater., Appl. Phys. Rev., Sci. Rep., among others of design space.

36 MATERIALS SCIENCE↗
DOE OSTI2016

The order of three lowest-energy states of the six-electron harmonium at small force constant

The order of low-energy states of six-electron harmonium is uncertain in the case of strong correlation, which is not a desired situation for the model system being considered for future testing of approximate methods of quantum chemistry. The computational study of these states has been carried out at the frequency parameter ω = 0.01, using the variational method with the basis of symmetry-projected, explicitly correlated Gaussian (ECG) lobe functions. It has revealed that the six-electron harmonium at this confinement strength is an octahedral Wigner molecule, whose order of states is different than in the strong confinement regime and does not agree with the earlier predictions. The results obtained for ω = 0.5 and 10 are consistent with the findings based on the Hund’s rules for the s{sup 2}p{sup 4} electron configuration. Substantial part of the computations has been carried out on the graphical processing units and the efficiency of these devices in calculation of the integrals over ECG functions has been compared with traditional processors.

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

Measuring the speed of light with ultra-compact radio quasars

In this paper, based on a 2.29 GHz VLBI all-sky survey of 613 milliarcsecond ultra-compact radio sources with 0.0035< z <3.787, we describe a method of identifying the sub-sample which can serve as individual standard rulers in cosmology. If the linear size of the compact structure is assumed to depend on source luminosity and redshift as l {sub m} = l L {sup β} (1+ z ) {sup n} , only intermediate-luminosity quasars (10{sup 27} W/Hz< L < 10{sup 28} W/Hz) show negligible dependence (| n |≅ 10{sup −3}, |β|≅ 10{sup −4}), and thus represent a population of such rulers with fixed characteristic length l =11.42 pc. With a sample of 120 such sources covering the redshift range 00.46< z <2.8, we confirm the existence of dark energy in the Universe with high significance under the assumption of a flat universe, and obtain stringent constraints on both the matter density Ω {sub m} =0.323{sup +0.245}{sub −0.145} and the Hubble constant H {sub 0}=66.30{sup +7.00}{sub −8.50} km sec{sup −1} Mpc{sup −1}. Finally, with the angular diameter distances D {sub A} measured for quasars extending to high redshifts (0 z ∼ 3.), we reconstruct the D {sub A} ( z ) function using the technique of Gaussian processes. This allows us to identify the redshift corresponding to the maximum of the D {sub A} ( z ) function: 0 z {sub m} =1.7 and the corresponding angular diameter distance D {sub A} ( z {sub m} )=1719.01±43.46 Mpc. Similar reconstruction of the expansion rate function H ( z ) based on the data from cosmic chronometers and BAO gives us H ( z {sub m} )=176.77±6.11 km sec{sup −1} Mpc{sup −1}. These measurements are used to estimate the speed of light: c =3.039(±0.180)× 10{sup 5} km/s. This is the first measurement of the speed of light in a cosmological setting referring to the distant past.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Active deep kernel learning of molecular properties from structural embeddings

As vast databases of chemical identities become increasingly available, the challenge shifts to how we effectively explore and leverage these resources to study molecular properties. This paper presents an active learning approach for molecular discovery using deep kernel learning (DKL), demonstrated on the QM9 dataset. DKL links structural embeddings directly to properties, creating organized latent spaces that prioritize relevant property information. By iteratively recalculating embedding vectors in alignment with target properties, DKL uncovers concentrated maxima representing key molecular properties and reveals unexplored regions with potential for innovation. This approach underscores DKL’s potential in advancing molecular research and discovery.

Artificial neural networks↗
DOE OSTI2025

Mapping causal patterns in crystalline solids

The evolution of the atomic structures of the combinatorial library of Sm-substituted thin film BiFeO 3 along the phase transition boundary from the ferroelectric rhombohedral phase to the non-ferroelectric orthorhombic phase is explored using scanning transmission electron microscopy. Localized properties, including polarization, lattice parameter, and chemical composition, are parameterized from atomic-scale imaging, and their causal relationships are reconstructed using a linear non-Gaussian acyclic model. This approach is further extended to explore the spatial variability of the causal coupling using the sliding window transform method, which revealed that new causal relationships emerged at both the expected locations, such as domain walls and interfaces, and at additional regions forming clusters in the vicinity of the walls or spatially distributed features. While the exact physical origins of these relationships are unclear, they likely represent nanophase-separated regions in the morphotropic phase boundaries. Overall, we posit that an in-depth understanding of complex disordered materials away from thermodynamic equilibrium necessitates understanding not only the generative processes that can lead to observed microscopic states but also the causal links between multiple interacting subsystems.

Causal inference↗
DOE OSTI2025

2025 ASMS Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS

Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2017

Group Additivity Determination for Oxygenates, Oxonium Ions, and Oxygen-Containing Carbenium Ions

Bio-oil produced from biomass fast pyrolysis often requires catalytic upgrading to remove oxygen and acidic species over zeolite catalysts. The elementary reactions in the mechanism for this process involve carbenium and oxonium ions. In order to develop a detailed kinetic model for the catalytic upgrading of biomass, rate constants are required for these elementary reactions. The parameters in the Arrhenius equation can be related to thermodynamic properties through structure-reactivity relationships, such as the Evans-Polanyi relationship. For this relationship, enthalpies of formation of each species are required, which can be reasonably estimated using group additivity. However, the literature previously lacked group additivity values for oxygenates, oxonium ions, and oxygen-containing carbenium ions. In this work, 71 group additivity values for these types of groups were regressed, 65 of which had not been reported previously and six of which were newly estimated based on regression in the context of the 65 new groups. Heats of formation based on atomization enthalpy calculations for a set of reference molecules and isodesmic reactions for a small set of larger species for which experimental data was available were used to demonstrate the accuracy of the Gaussian-4 quantum mechanical method in estimating enthalpies of formation for species involving the moieties of interest. Isodesmic reactions for a total of 195 species were constructed from the reference molecules to calculate enthalpies of formation that were used to regress the group additivity values. The results showed an average deviation of 1.95 kcal/mol between the values calculated from Gaussian-4 and isodesmic reactions versus those calculated from the group additivity values that were newly regressed. Importantly, the new groups enhance the database for group additivity values, especially those involving oxonium ions.

09 BIOMASS FUELS↗
DOE OSTI2024

Earth System Reanalysis in Support of Climate Model Improvements

Recent climate model developments, established through increased model resolution, have led to substantial improvements in model simulations of the time-evolving, coupled Earth system and its subcomponents. However, regardless of resolution, climate models will always produce climate features and variability that differ from the real world and will be prone to biases. This is due to many remaining uncertainties, such as in parametric and structural model uncertainty, in the initial conditions prescribed, and in the prescribed (scenario) forcing which varies on decadal to centennial timescales. Further model improvements are expected to arise specifically from improved representation of physical processes realized through model-data fusion. This will create an unprecedented opportunity to better exploit a large array of Earth observations, from in situ measurements to weather radars and satellite observations, as the resolved scales of the models approach those of the observations. For this, climate DA will be the central tool to bring models and observations into consistency, by improving initial conditions, inferring uncertain model parameters and structure, and quantifying uncertainty. Generally, there will be advantages and complementarities of adjoint-based smoother approaches, ensemble-based filter approaches, or new ML-inspired approaches. Yet, the ever-increasing model resolution will present growing challenges arising from computational cost, calling for new ways of performing data assimilation and model optimization. Using the complementarity in a hybrid approach, blending tools and concepts from variational, ensemble and ML methods might be what is required in the future. In this context ML could be important to handle non-linear responses, and to better approximate non-Gaussian distributions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Cooling rate and dendrite spacing control in direct metal deposition printed Cu-Fe alloys

Over the last several decades, additive manufacturing (AM) has been primarily used for rapid prototyping or to create novel geometries that would be difficult or impossible to create by normal manufacturing methods. More recently, research has been focused on expanding the list of materials that can be made through additive manufacturing, opening a greater range of material properties for this manufacturing method. Due to the unusual conditions during AM, including the high cooling rates and voxel by voxel method of production, AM parts often have anisotropic material microstructures and properties. In this investigation, the laser power, composition, and nozzle head speed during direct metal deposition of copper-iron alloys was varied to understand how the grain structure within the printed parts could be changed and controlled. The resulting dendrite spacing was measured and compared to calculated cooling rates from Gaussian beams on flat plates under similar material and laser properties, which resulted in a cooling rate to dendrite spacing relationship following an inverse square root, as is found in other dendritic systems [Young and Kerkwood, Metall. Trans. A 6, 197–205 (1975)]. Thus, it is demonstrated that in the Cu-Fe system, dendrite spacing can be controlled through manipulation of printing parameters.

36 MATERIALS SCIENCE↗