Search NASA⌕ Search

SEARCH · Search NASA

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

A workflow for automatic generation and efficient refinement of individual pressure-dependent networks

Manual chemical kinetic model construction often requires modelists to at least implicitly guess all possible decay pathways and their relative fluxes for all important species. We present a workflow and associated tool for automatic generation and re finement of pressure-dependent networks that should enable modelists to efficiently and comprehensively identify decay pathways and estimate respective parameters for chemical species. This tool combines the capabilities of the Reaction Mechanism Generator (RMG) software to generate possible reaction paths, determine thermochemistry, approximate rate coefficients and estimate frequencies with the capabilities of the Arkane software to use quantum chemical parameters to compute thermochemistry and pressure dependent rates. A flux-based algorithm is used to decide which isomers to add to the network. Isomers added to the network are reacted to form other channels. When enough of the flux is accounted for by the isomers and bimolecular product channels, the generation process terminates to yield a comprehensive network. Network sensitivity analysis is applied to the network to identify the most important wells and barriers that can be refined using quantum chemistry calculations. Iterative refinement can then be used to achieve accuracy. The comprehensive network can be easily reduced using energy and flux-based algorithms to the most important channels and isomers in the network.

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

The C 2 H 4 O isomers in the oxidation of ethylene

A combined rotational spectroscopy, thermochemistry, and kinetic modeling study explores the mechanism of acetaldehyde (ethanal), vinyl alcohol (ethenol), and ethylene oxide (oxirane) formation in the oxidation of ethylene (ethene). Multiplexed quantitative detection of oxidation intermediates and products exiting a SiO 2 / SiC microreactor at 1700 K is demonstrated with the BrightSpec W-band chirped-pulse Fourier transform millimeter-wave spectrometer. The broadband rotational spectrum contains transitions of formaldehyde (CH 2 O), methoxy (CH 3 O), methanol (CH 3 OH), ketene (CH 2 CO), acetaldehyde (CH 3 CHO), syn-vinyl alcohol (syn- CH 2 CHOH), anti-vinyl alcohol (anti-CH 2 CHOH), oxirane (c-CH 2 OCH 2 ), propyne (CH 3 CCH), and syn-propanal (syn-CH 3 CH 2 CHO). We focus on the three C 2 H 4 O species and deduce their concentration ratio [CH 3 CHO]:[CH 2 CHOH]:[c-CH 2 OCH 2 ] = 1:0.7(2):0.06(2) by comparing the observed line intensities to those simulated with the PGOPHER software. Detailed thermochemistry of the C 2 H 4 O isomers and conformers is provided by Active Thermochemical Tables. The observed excess concentrations of vinyl alcohol and oxirane relative to the more stable acetaldehyde compared to the equilibrium ratio at 1700 K (1:0.087:0.000024), point to direct chemical pathways to these higher energy isomers. The mechanism for the formation of the three C 2 H 4 O isomers is analyzed using a 0-D homogenous reactor kinetics simulation for ethylene oxidation. The ratios of the C 2 H 4 O isomers concentrations predicted by the kinetic model are compared to the experimental values.

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

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. Here, the transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Developing Scaling Relationships for Molecular Electrocatalysis through Studies of Fe-Porphyrin-Catalyzed O 2 Reduction

The oxygen reduction reaction (ORR) is a multiproton/multielectron transformation in which dioxygen (O 2 ) is reduced to water or hydrogen peroxide and serves as the cathode reaction in most fuel cells. The ORR (O 2 + 4e – + 4H + → 2H 2 O) involves up to nine substrates and thus requires navigating a complicated reaction landscape, typically with several high-energy intermediates. Many catalysts can perform this reaction, though few operate with fast rates and at low overpotentials (close to the thermodynamic potential). Attempts to optimize these parameters, both in homogeneous and heterogeneous electrocatalytic systems, have focused on modifying catalyst design and understanding kinetic/thermodynamic relationships between catalytic intermediates. One such method for analyzing and predicting catalyst reactivity and efficiency has been the development of “molecular scaling relationships”. In this work, we share our experience deriving and utilizing molecular scaling relationships for soluble, iron-porphyrin-catalyzed O 2 reduction in organic solvents. These relationships correlate turnover frequencies (TOF max ) and effective overpotentials (η eff ), properties uniquely defined for homogeneous catalysts. Following a general introduction of scaling relationships for both homogeneous and heterogeneous electrocatalysis, we describe the components of such scaling relationships: (i) the overall thermochemistry of the reaction and (ii) the rate and rate law of the catalyzed reaction. We then show how connecting these thermodynamic and kinetic parameters reveals multiple molecular scaling relationships for iron-porphyrin-catalyzed O 2 reduction. For example, the log(TOF max ) responds steeply to changes in ηeff that result from different catalyst reduction potentials (18.5 decades in TOF max /V in η eff ) but much less dramatically to changes in η eff that arise from varying the pKa of the acid buffer (5.1 decades in TOF max /V in η eff ). Thus, a single scaling relationship is not always sufficient for describing molecular electrocatalysis. This is particularly evident when the catalyst identity and reaction conditions are coupled. Using these multiple scaling relationships, we demonstrate that the metrics of turnover frequency and effective overpotential can be predictably tuned to achieve faster rates at lowered overpotentials. This Account uses a collection of related stories describing our research on soluble iron-porphyrin-catalyzed ORR to show how molecular scaling relationships can be derived and used for any electrocatalytic reaction. Such scaling relationships are powerful tools that connect the thermochemistry, mechanism, and rate law for a catalytic system. We hope that this collection shows the utility and simplicity of the molecular scaling approach for understanding catalysis, for enabling direct comparisons between catalyst systems, and for optimizing catalytic processes.

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

Uncertainty Quantification and Error Propagation in the Enthalpy and Entropy of Surface Reactions Arising from a Single DFT Functional

Accounting for parametric uncertainty in models is essential for quantifying the models’ predictive ability. Recently, approaches have been introduced to estimate parametric uncertainty in kinetic models while accounting for correlations among energy parameters. However, correlations have been estimated indirectly and correlations in entropies have not been accounted for. For surface-catalyzed microkinetic models of >C2 (more than two carbon-containing) molecules, which consist of thousands of reaction steps and intermediate surface species, first-principles density functional theory (DFT) is costly, and thus, estimation of thermochemistry and reaction barriers requires surrogate methods of DFT, such as group additivity and Brønsted–Evans–Polanyi relationships, respectively. For such parametrization, model uncertainty is unclear. This work develops a framework to overcome these gaps using group additivity and a single DFT functional. We estimate correlations in parameters of kinetic models and quantify uncertainty for thermochemistry, reaction barriers, reaction paths, and ultimately reaction rates, accounting also for the contribution of entropic uncertainty. Furthermore, the approach is illustrated on propane combustion and ethane oxidative dehydrogenation reactions.

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

Proton-Coupled Defects Impact O–H Bond Dissociation Free Energies on Metal Oxide Surfaces

Proton-coupled electron transfer (PCET) reactions on metal oxides require coupling between proton transfer at the solid-liquid interface and electron transfer involving defects at or near the band edge. Herein, hybrid functional periodic density functional theory is used to elucidate the impact of proton-coupled defects on the bond dissociation free energies (BDFEs) of O-H bonds on anatase TiO 2 surfaces. These O-H BDFEs are directly related to interfacial PCET thermochemistry. Comparison between geometrically similar O-H bonds associated with different defect types, namely conduction d-band electrons or valence p-band holes, reveals that the BDFEs differ by ~81 kcal/mol (3.50 eV), comparable to the wide TiO 2 band gap. These differences are shown to be determined primarily by differences in electron transfer driving forces, which are analyzed by using band energies and inner-sphere reorganization energies within a Marcus theory framework. Furthermore, these fundamental insights about the impact of proton-coupled defects on PCET thermochemistry at semiconductor surfaces have broad implications for electrocatalysis.

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

A graph-based network for predicting chemical reaction pathways in solid-state materials synthesis

Abstract Accelerated inorganic synthesis remains a significant challenge in the search for novel, functional materials. Many of the principles which enable “synthesis by design” in synthetic organic chemistry do not exist in solid-state chemistry, despite the availability of extensive computed/experimental thermochemistry data. In this work, we present a chemical reaction network model for solid-state synthesis constructed from available thermochemistry data and devise a computationally tractable approach for suggesting likely reaction pathways via the application of pathfinding algorithms and linear combination of lowest-cost paths in the network. We demonstrate initial success of the network in predicting complex reaction pathways comparable to those reported in the literature for YMnO 3 , Y 2 Mn 2 O 7 , Fe 2 SiS 4 , and YBa 2 Cu 3 O 6.5 . The reaction network presents opportunities for enabling reaction pathway prediction, rapid iteration between experimental/theoretical results, and ultimately, control of the synthesis of solid-state materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter ($α$) was set to recover the exact dissociation limit of H 2 in a minimal basis set. Alternatively $α$ can be viewed as a regularization parameter, where each value of $α$ represents a valid variant of BW-s2, which we denote as BW-s2($α$). In this work, we semi-empirically optimize $α$ for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal $α$ parameter ($α$ = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2($α$) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2($α$) is iterative $\mathcal{O}$($N^5$), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

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

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermochemically-informed mass transport model for interdiffusion of U and Zr in irradiated U-Pu-Zr fuel with fission products

A new formulation for redistribution of constituents in metallic U-Pu-Zr nuclear fuel is presented that can incorporate the contributions from fission products. The formulation is based on the thermodynamic driving forces derived from the generalized chemical potential that includes effects of composition and temperature. As a result, the redistribution model can readily account for the composition changes due to the generation of fission products while using only a limited set of transport coefficients. The thermodynamic model for the metallic fuel and thermochemistry solver Thermochimica were coupled with the nuclear fuel performance code BISON to implement the redistribution model. The simulations reproduce an experimentally observed Zr-depletion zone in the mid-radius region of an irradiated fuel slug. The generation of fission products during burnup is shown to have a stabilizing effect on the fuel chemistry, slowing the rate of U-Zr interdiffusion and reducing the size of the Zr-depleted zone. Here, the overall objective of the work is to advance predictive capabilities of fission product behaviour in U-Zr metallic fuel in the context of fuel performance and safety.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Association Kinetics for Perfluorinated n -Alkyl Radicals

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

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

r 2 SCAN-D4: Dispersion corrected meta-generalized gradient approximation for general chemical applications

We combine a regularized variant of the strongly constrained and appropriately normed semilocal density functional with the latest generation semi-classical London dispersion correction. The resulting density functional approximation r 2 SCAN-D4 has the speed of generalized gradient approximations while approaching the accuracy of hybrid functionals for general chemical applications. We demonstrate its numerical robustness in real-life settings and benchmark molecular geometries, general main group and organo-metallic thermochemistry, and non-covalent interactions in supramolecular complexes and molecular crystals. Main group and transition metal bond lengths have errors of just 0.8%, which is competitive with hybrid functionals for main group molecules and outperforms them for transition metal complexes. The weighted mean absolute deviation (WTMAD2) on the large GMTKN55 database of chemical properties is exceptionally small at 7.5 kcal/mol. This also holds for metal organic reactions with an MAD of 3.3 kcal/mol. The versatile applicability to organic and metal–organic systems transfers to condensed systems, where lattice energies of molecular crystals are within the chemical accuracy (errors <1 kcal/mol).

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

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗
DOE OSTI2021

Low-Cost Radiant Heater for Rapid Response, High-Temperature Heating

High-temperature processing has an irreplaceable role in many research and industrial applications. Despite remarkable development spanning over a century, the pursuit of even higher thermal flux density and more rapid thermal transients has not slowed down. As part of the ongoing energy evolution, many industrial applications are transitioning from direct combustion of fossil fuels as primary energy sources to increasing electrification, capable of adapting to renewable power grids. Thus, there is an emerging need for electrical heaters that can replace burners and supply the heat demand, especially at the highest temperatures. In this study, we report on a radiant heater design that can achieve cyclic heating/cooling rates of up to 400 K min –1 and a temperature range in excess of 1,800 K, comparable to those of commercial infrared gold image furnaces, at high surface and volumetric power densities. The heater consists of a modular unit of incandescent tungsten filament and is enclosed in an evacuated ceramic envelope, chemically inert, tolerant of thermal shock, and impervious to gasses. The material and manufacture cost of such heaters, which is estimated at ~$0.05/W, is less than 0.03% of that for infrared gold image furnaces, which is at >$2/W. Tests of more than 10,000 demanding cycles (high temperature and high heating/cooling rate) over 350 h of total operational time and in different temperature ranges confirm the robust performance of radiant heater prototypes. The design is widely applicable to high-temperature reactor and furnace designs. In thermochemistry research and practice, these radiant heaters could offer multiple benefits compared to solar simulators, lasers, infrared gold furnaces, ceramic heaters, or direct concentration of solar input.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Defying Thermodynamics: Stabilization of Alane Within Covalent Triazine Frameworks for Reversible Hydrogen Storage

The highly unfavorable thermodynamics of direct aluminum hydrogenation can be overcome by stabilizing alane within a nanoporous bipyridine-functionalized covalent triazine framework (AlH 3 @CTF-bipyridine). This material and the counterpart AlH 3 @CTF-biphenyl rapidly desorb H 2 between 95 and 154°C, with desorption complete at 250°C. Sieverts measurements, 27 Al MAS NMR and 27 Al{ 1 H} REDOR experiments, and computational spectroscopy reveal that AlH 3 @CTF-bipyridine dehydrogenation is reversible at 60°C under 700 bar hydrogen, >10 times lower pressure than that required to hydrogenate bulk aluminum. DFT calculations and EPR measurements support an unconventional mechanism whereby strong AlH 3 binding to bipyridine results in single-electron transfer to form AlH 2 (AlH 3 ) n clusters. Here the resulting size-dependent charge redistribution alters the dehydrogenation/rehydrogenation thermochemistry, suggesting a novel strategy to enable reversibility in high-capacity metal hydrides.

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

Defying Thermodynamics: Stabilization of Alane Within Covalent Triazine Frameworks for Reversible Hydrogen Storage

Metastable metal hydrides such as AlH 3 have many attractive features as hydrogen storage media, but generally require complex reaction schemes for regeneration following H 2 release. Here in this paper, we demonstrate that the highly unfavorable thermodynamics of direct aluminum hydrogenation can be overcome by stabilizing alane within a nanoporous bipyridine-functionalized Covalent Triazine Framework (AlH 3 @CTF-bipyridine). This material and the counterpart AlH 3 @CTF-biphenyl rapidly desorb H 2 between 95 and 154 °C, with desorption complete at 250 °C. Sieverts measurements, 27 Al MAS NMR and 27 Al{ 1 H} REDOR experiments, and computational spectroscopy reveal that AlH 3 @CTF-bipyridine dehydrogenation is reversible at 60 °C under 700 bar hydrogen, >10 times lower pressure than that required to hydrogenate bulk aluminum. DFT calculations and EPR measurements support an unconventional mechanism whereby strong AlH 3 binding to bipyridine results in single-electron transfer to form AlH 2 (AlH 3 ) n clusters. The resulting size-dependent charge redistribution alters the dehydrogenation/rehydrogenation thermochemistry, suggesting a novel strategy to enable reversibility in high-capacity metal hydrides.

Coordination Chemistry↗