Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermochemistry”

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

Active Thermochemical Tables: Should the enthalpy of formation of gas phase boron atom be revised?

The Active Thermochemical Tables approach produces the enthalpy of formation of gas phase boron atom: Δ f H° 298 (B (g) ) = 570.48±0.61 kJ/mol and Δ f H° 0 (B (g) ) = 565.38±0.61 kJ/mol. This is about 5 kJ/mol higher and nearly an order of magnitude more accurate than the CODATA value. Here, while the ATcT value is in excellent agreement with the revisions proposed by Bauschlicher, Martin, and Taylor [J. Phys. Chem. A 103 (1999) 7715] and by Karton and Martin [J. Phys. Chem. A 111 (2007) 5936], it invalidates several earlier theoretical revisions that are too high by up to 5 kJ/mol.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Exploring diversion-pathway analysis of a generic molten-salt fast reactor using multiphysics informed signatures

Molten salt reactors are being explored by multiple commercial ventures due to their inherent safety features, flexibility in fuel sources, and high fuel utilization and thermal efficiency. The continual flow of fuel salt, large fissile quantities present, and ability to add or divert material due to the liquid nature introduces new challenges for international safeguards. To understand how international safeguards should be applied, it is important to capture the inherent multi-physics nature of a molten salt reactor. This work examines a generic molten salt fast reactor to understand how potential diversion scenarios would affect the concentration of radionuclides in the primary and auxiliary systems. Three types of diversion were examined: a slow drip of fuel salt, gaseous plutonium extraction, and uranium metal plating. The analysis determined that several key isotopes become statistically significant once diversion begins, indicating that detection of such diversion cases would be possible through measuring specific signatures such as gamma spectra.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstraction by $CH_3O\dot{O}$ radicals

Hydrogen atom abstraction by methyl peroxy ($CH_3O\dot{O}$) radicals can play an important role in gasoline/ethanol interacting chemistry for fuels that produce high concentrations of methyl radicals. Detailed kinetic reactions for hydrogen atom abstraction by $CH_3O\dot{O}$ radicals from the components of FGF-LLNL (a gasoline surrogate) including cyclopentane, toluene, 1-hexene, n -heptane, and isooctane have been systematically studied in this work. Here, the M06-2X/6-311 ++ G(d,p) level of theory was used to obtain the optimized structure and vibrational frequency for all stationary points and the low-frequency torsional modes. The 1-D hindered rotor treatment for low-frequency torsional modes was treated at M06-2X/6-31G level of theory. The UCCSD(T)-F12a/cc-pVDZ-F12 and QCISD(T)/CBS level of theory were used to calculate single point energies for all species. High pressure limiting rate constants for all hydrogen atom abstraction channels were performed using conventional transition state theory with unsymmetric tunneling corrections. Individual rate constants are reported in the temperature range from 298.15 to 2000 K. Our computed results show that the abstraction of allylic hydrogen atoms from 1-hexene is the fastest at low temperatures. When the temperature increases, the hydrogen atom abstraction reaction channel at the primary alkyl site gradually becomes dominant. Thermodynamics properties for all stable species and high-pressure limiting rate constants for each reaction pathway obtained in this work were incorporated into the latest gasoline surrogate/ethanol model to investigate the influence of the rate constants calculated here on model predicted ignition delay times.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Status and Challenges of Density Functional Theory

Here, we discuss some of the challenges facing density functional theory (DFT) and recent progress in DFT for both ground and excited electronic states. We discuss key aspects of the results we have been able to obtain with the strategy of designing density functionals to have various ingredients and functional forms that are then optimized to accurately predict various types of properties and systems with as much universality as possible. Finally, we make specific recommendations of approximate density functionals that are well suited for particular kinds of applications.

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

Thermochemical Transformation of Calcium during Biomass Burning and the Effects on Postfire Aqueous Dissolution of Macronutrients

Calcium is commonly the most abundant element in fire residues and its speciation largely determines the geochemical properties of fire residues and their effects on postfire soil chemistry. To explore the effects of biomass composition and fire conditions on ash Ca speciation, this study characterizes the speciation of Ca in charcoal and ash samples that were derived from different plant compartments and thermal conditions, using Ca K-edge X-ray absorption near edge spectroscopy. Results showed that biomass contains abundant organic Ca complexes, which were mineralized into fairchildite and calcite after heating at 450 to 600 °C and then CaO, as temperature increased to 750 °C. Apatite could be an abundant Ca species in fire residues if the Ca/P molar ratio of the biomass is small (<2). The mineralization of organic Ca to the identified Ca minerals during burning was negligibly affected by the oxygen level. Calcium speciation in prescribed fire residues resembled that of biomass ash burned at 550 °C with similar Ca/P molar ratios. Batch experiments showed that macronutrients (Ca, Mg, K, and P) were differentially released, as a result of different solubility of minerals in ashes and reprecipitation of minerals. The aqueous solubility of Ca, Mg, and P decreased as pH increased from 5 to 9, while K showed no pH dependency and was almost completely soluble. Results from this study improve our understanding of the chemistry of fire residues and their geochemical behaviors, which can help evaluate the impact of fire on postfire soil properties and macronutrient cycling.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Formation of 6H-Ba 3 Ce 0.75 Mn 2.25 O 9 during Thermochemical Reduction of 12R-Ba 4 CeMn 3 O 12 : Identification of a Polytype in the Ba(Ce,Mn)O 3 Family

The resurgence of interest in a hydrogen economy and the development of hydrogen-related technologies has initiated numerous research and development efforts aimed at making the generation, storage, and transportation of hydrogen more efficient and affordable. Solar thermochemical hydrogen production (STCH) is a process that potentially exhibits numerous benefits such as high reaction efficiencies, tunable thermodynamics, and continued performance over extended cycling. Although CeO 2 has been the de facto standard STCH material for many years, more recently 12R-Ba 4 CeMn 3 O 12 (BCM) has demonstrated enhanced hydrogen production at intermediate H 2 /H 2 O conditions compared to CeO 2 , making it a contender for large-scale hydrogen production. However, the thermo-reduction stability of 12R-BCM dictates the oxygen partial pressure (pO 2 ) and temperature conditions optimal for cycling. In this study, we identify the formation of a 6H-BCM polytype at high temperature and reducing conditions, experimentally and computationally, as a mechanism and pathway for 12R-BCM decomposition. 12R-BCM was synthesized with high purity and then controllably reduced using thermogravimetric analysis (TGA). Synchrotron X-ray diffraction (XRD) data is used to identify the formation of a 6H-Ba 3 Ce 0.75 Mn 2.25 O 9 (6H-BCM) polytype that is formed at 1350 degrees C under strongly reducing pO 2 . Density functional theory (DFT) total energy and defect calculations show a window of thermodynamic stability for the 6H-polytype consistent with the XRD results. These data provide the first evidence of the 6H-BCM polytype and could provide a mechanistic explanation for the superior water-splitting behaviors of 12R-BCM.

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

Thermochemical Properties of a Water-Soluble Zr(IV)-Flavonolate Complex for Fluoride Sensing

Detection of fluoride in water samples is critical for monitoring many chemical, environmental, and biological processes, however, these measurements often rely on costly instrumentation. Here, this situation has motivated the design of simple to use, fieldable fluoride sensors. Previously, a fluorometric sensor based on (EDTA)­Zr­(H 2 O) 2 (EDTA = ethylenediamine-N,N,N′,N′-tetraacetate) and 3-hydroxyflavone (FlvH) was reported for fluoride detection in ethanol-H 2 O mixtures. Herein, we report the synthesis and characterization of the previously proposed Zr­(IV)-flavonolate complex as the dimethylammonium salt, [(EDTA)­Zr­(Flv)]­[NMe 2 H 2 ] (Flv = 3-hydroxyflavonolate), and detail its utility for aqueous fluoride sensing. This species exhibits an intense blue fluorescence upon excitation with near-UV light (λ ex = 390 nm; λ em = 460 nm) that is quenched in the presence of NaF (100 mM acetate buffer, pH = 5) with good sensitivity (LOD = 20.9 nM; [F – ] from 10 –8 to 10 –4 M) under neat aqueous conditions. The equilibrium constant (K eq = 1.26 ± 0.00184 × 10 8 M –2 ) and free energy (Δ$G^°_{\textrm{rxn}}$ = –11.0 kcal mol –1 ) were determined under the conditions used for fluoride sensing, indicating a highly favorable reaction with F – in water. Finally, we developed paper-based sensing strips for qualitative fluoride sensing by immobilizing the sensor in a deposited gel enabling good selectivity and reduced interference due to slower ion diffusion.

aqueous chemistry↗
DOE OSTI2021

Reaction Mechanism Generator v3.0: Advances in Automatic Mechanism Generation

In chemical kinetics research, kinetic models containing hundreds of species and tens of thousands of elementary reactions are commonly used to understand and predict the behavior of reactive chemical systems. Reaction Mechanism Generator (RMG) is a software suite developed to automatically generate such models by incorporating and extrapolating from a database of known thermochemical and kinetic parameters. Here, we present the recent version 3 release of RMG and highlight improvements since the previously published description of RMG v1.0. Most notably, RMG can now generate heterogeneous catalysis models in addition to the previously available gas- and liquid-phase capabilities. For model analysis, new methods for local and global uncertainty analysis have been implemented to supplement first-order sensitivity analysis. The RMG database of thermochemical and kinetic parameters has been significantly expanded to cover more types of chemistry. The present release includes parallelization for faster model generation and a new molecule isomorphism approach to improve computational performance. RMG has also been updated to use Python 3, ensuring compatibility with the latest cheminformatics and machine learning packages. Overall, RMG v3.0 includes many changes which improve the accuracy of the generated chemical mechanisms and allow for exploration of a wider range of chemical systems.

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

Energetic and Electronic Properties of UX +/0/– for X = Li and Be and Comparison of the Properties of the Uranium Atom Binding to 2nd Row Elements Li–F

The bonding and spectroscopic properties of ULi +/0/– and UBe +/0/– to complete the series for UX +/0/– for X = Li to F were investigated by high-level ab initio SO-CASPT2 and CCSD(T) electronic structure calculations. The low-lying spin–orbit states were obtained at the SA-CASPT2/aQ-PP level; bond dissociation energies (BDEs), ionization energies (IEs), adiabatic electronic affinities (AEAs), and vertical detachment energies (VDEs) were calculated at the Feller-Peterson-Dixon (FPD) level. A dense manifold of low-lying states was predicted for ULi +/0/– and UBe +/0/– . Here, the calculated BDEs for ULi (37.7 kJ/mol) and UBe (8.0 kJ/mol) show that UBe is weakly bound. For redox processes, the BDEs increased for ULi + (109.3 kJ/mol), ULi – (47.4 kJ/mol), UBe + (35.6 kJ/mol), and UBe – (72.3 kJ/mol). The IE(ULi) = 4.650 eV is lower than IE(Li); the IE(UBe) = 5.901 eV is close to the IE(U) and to the IEs of UB, UC, UN, UO, and UF. The AEAs of ULi (0.708 eV) and UBe (0.989 eV) are lower than those for UB, UC, UN, and UO but higher than that for EA(UF). Natural bond orbital (NBO) calculations show that ULi has the 5f 3 6d 1 7s 2 configuration for U and 2s 1 for Li, with a small partial negative charge slightly delocalized on U. UBe arises from the U(5f 3 6d 1 7s 2 ) and Be(2s 2 ) electron configurations with no charge separation. The same calculations were made for WX (X = Li, Be, C–F) to enable detailed comparisons of the properties for UX with WX (X = Li–F). For WX, BDE(WX) is higher than that for UX for X = Li to N and lower than BDE(UX) for X = O and F, mostly due to the higher IE of W than U as ionic character becomes more important going from Li to F.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Discrepancies between Theory and Experiment in Determining the Ionization Energy of NF 3

High-accuracy ab initio thermochemical predictions for the ionization energy of NF 3 , the barrier height (to inversion) of NF 3 + , and the dissociative ionization threshold of NF 3 to NF 2 + + F are presented and incorporated into Active Thermochemical Tables. The adiabatic ionization energy of the first ionization band of NF 3 , calculated at 12.647 ± 0.010 eV, is at odds with previous experimental interpretations by nearly 0.36 eV due to unfavorable Franck-Condon factors associated with this transition. The barrier (to inversion) height is calculated to be about 0.6 eV lower in energy than the prior interpretation, which instigates a discussion of the supposed vibrational structure of the first ionization band of NF 3 . Updated assignments of the photoelectron spectrum are proposed, and the loss in vibrational spacing on the high-energy side of the experimental ionization band is discussed. Rudimentary anharmonic Franck-Condon simulations qualitatively reproduce the broad spectral features observed in experiment.

Active Thermochemical Tables↗
DOE OSTI2025

The Electronic Structure of Zirconium and Hafnium Monochalcogenides

High-level ab initio CCSD(T) and spin–orbit icMRCI+Q calculations were used to predict potential energy curves (PECs) for the lowest-lying states of ZrO, ZrS, HfO, and HfS. The prediction of the ground state is basis set dependent at the icMRCI+Q level for ZrO and ZrS due to the small singlet–triplet splitting between the lowest 1 Σ + and 3 Δ states. CCSD(T) with a spin orbit correction predicted the 1 Σ + ground state in agreement with experiment. New all-electron basis sets were developed for Hf to improve the results over those predicted by use of effective core potentials (ECPs) that subsume the 4f electrons into the definition of the core. The use of the new DK-4f basis sets rather than ECPs became more important for HfO and HfS where there is a lack of a good core–valence separation. icMRCI+Q, CCSD(T), and DFT calculations for the spectroscopic parameters of ZrO, ZrS, HfO, and HfS were benchmarked with available experimental data. Bond dissociation energies (BDEs) of these four systems were calculated at the Feller–Peterson–Dixon (FPD) level to be 762.1 (ZrO), 543.5 (ZrS), 803.8 (HfO), and 575.1 kJ/mol (HfS), in excellent agreement with experiment. The HfS BDE was remeasured using the R3PI method, providing an updated experimental measurement of D 0 (HfS) = 5.978 ± 0.002 eV = 576.8 ± 0.2 kJ/mol. This experimental value, combined with experimental measurements of the ionization energies of Hf and HfS, gives the cationic BDE of D 0 (Hf + -S) = 5.124 ± 0.002 eV = 494.4 ± 0.2 kJ/mol.

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

Uranium Hexafluoride Hydrolysis Reaction Dynamics from Cryogenic Layering, FTIR Spectroscopy, and Isotopic Substitution

The first direct evidence that the hydrolysis reaction of uranium hexafluoride (UF 6 ) follows multiple competing pathways which are driven by the ratio of water to UF 6 , temperature, and isotopic composition is presented. Using temperature dependent infrared spectroscopy, it is shown the hydrolysis can be prevented at temperatures below 150 K, and that water-rich environments promote the formation of uranium oxyfluoride intermediates. Spectral shifts reveal isomeric transitions and the growth of polymeric species, with reaction reversibility observed at high water concentrations. Additionally, controlled heating rates affect the emergence of intermediates. The final particulate product consistently forms as uranyl fluoride hydrate, though its morphology and spectral signature vary with reaction conditions and annealing. These findings help clarify long-standing uncertainties surrounding UF 6 hydrolysis.

Chemical reactions↗
DOE OSTI2021

Adsorbate Partition Functions via Phase Space Integration: Quantifying the Effect of Translational Anharmonicity on Thermodynamic Properties

A new method for computing anharmonic thermophysical properties for adsorbates on metal surfaces is presented. Classical Monte Carlo phase space integration is performed to calculate the partition function for the motion of a hydrogen atom on Cu(111). Here, a minima-preserving neural network potential energy surface is used within the integration routine. Two different sampling schema for generating the training data are presented, and two different density functionals are used. The results are benchmarked against direct state counting results by using discrete variable representation. The phase space integration results are in excellent quantitative agreement with the benchmark results. Additionally, both the discrete variable representation and the phase space integration results confirm that the motion of H on Cu(111) is highly anharmonic. The results were applied to calculate the free energy of dissociative adsorption of H 2 and the resulting Langmuir isotherms at 400, 800, and 1200 K in a partial pressure range of 0–1 bar. It shows that the anharmonic effects lead to significantly higher predicted surface site fractions of hydrogen.

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

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

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

Bond Dissociation Energy, Ionization Energy, and Electronic Structure of Thorium Dimer

Diatomic thorium, Th 2 , has been investigated using a laser ablation, supersonic expansion source to produce the molecule and resonant two-photon ionization spectroscopy to measure its bond dissociation energy (BDE) and ionization energy (IE). The molecule has a high density of states in the vicinity of its bond dissociation energy, leading to rapid predissociation as soon as this energy is exceeded. The BDE is identified from this predissociation threshold as D 0 (Th 2 ) = 2.857(7) eV, where the assigned error limit is provided in parentheses in units of the last quoted digit. Similarly, the one-photon ionization threshold has been measured, providing the ionization energy IE(Th 2 ) = 5.042(4) eV. Together with a thermochemical cycle and the atomic ionization energy, these values provide the BDE of the cation, giving D 0 (Th 2 + ) = 4.122(8) eV. Computations show that Th 2 has three nearly degenerate low-lying electronic states (1 3 Σ u + , 1 1 Σ g + , and 1 3 Δ g ) with bonding dominated by 7s and 6d orbitals, indicating predominantly transition-metal-like behavior. The 1 3 Σ u + state exhibits a triple bond, whereas the 1 1 Σ g + and 1 3 Δ g states possess quadruple-bond character and correspondingly shorter bonds. Although 1 3 Σ u + is predicted to be the lowest state without spin–orbit coupling, the large spin–orbit stabilization of the 1 3 Δ g state makes its Ω = 1 g component the ground state. Furthermore, the calculated dissociation energy (2.840 eV) and ionization energy of Th 2 (5.098 eV) are both in excellent agreement with experiment.

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

Consequences of sp 2 –sp 3 boron isomerization in supercooled liquid borates

In this study, time-resolved high-energy synchrotron x-ray total scattering measurements on supercooled molten lithium metaborate (LiBO 2 ) reveal an isomerization reaction involving conversion of trigonal sp 2 boron to tetrahedral sp 3 boron during quenching and glass formation. Van't Hoff analysis yields an accurate enthalpy change, ΔH = 21(1) kJ mol –1 boron, from which we develop an analytical model for the sp 3 isomer fraction and its contribution to configurational heat capacity (C p conf ) and entropy as a function of temperature and composition. Isomerization constitutes 40% of the total calorimetric C p conf at the glass transition for LiBO 2 and directly contributes to the observed rise in liquid fragility with the lithium content.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

A guided ion beam investigation of UO 2 + thermodynamics and f orbital participation: Reactions of U + + CO 2 , UO + + O 2 , and UO + + CO

A guided ion beam tandem mass spectrometer was employed to study the reactions of U + + CO 2 , UO + + O 2 , and the reverse of the former, UO + + CO. Reaction cross sections as a function of kinetic energy over about a three order of magnitude range were studied for all systems. The reaction of U + + CO 2 proceeds to form UO + + CO with an efficiency of 118% ± 24% as well as generating UO 2 + + C and UCO + + O. The reaction of UO + + O 2 forms UO 2 + in an exothermic, barrierless process and also results in the collision-induced dissociation of UO + to yield U + . In the UO + + CO reaction, the formation of UO 2 + in an endothermic process is the dominant reaction, but minor products of UCO + + O and U + + (O + CO) are also observed. Analysis of the kinetic energy dependences observed provides the bond energies, D 0 (U + –O) = 7.98 ± 0.22 and 8.05 ± 0.14 eV, D 0 (U + –CO) = 0.73 ± 0.13 eV, and D 0 (OU + –O) = 7.56 ± 0.12 eV. The values obtained for D 0 (U + –O) and D 0 (OU + –O) agree well with the previously reported literature values. To our knowledge, this is the first experimental measurement of D 0 (U + –CO). Furthermore, an analysis of the oxide bond energies shows that participation of 5f orbitals leads to a substantial increase in the thermodynamic stability of UO 2 + relative to ThO 2 + and especially transition metal dioxide cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗