Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio thermodynamics”

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

High-Temperature Thermodynamics of Uranium from Ab Initio Modeling

We present high-temperature thermodynamic properties for uranium in its γ phase (γ-U) from first-principles, relativistic, and anharmonic theory. The results are compared to CALPHAD modeling. The ab initio electronic structure is obtained from density-functional theory (DFT) that includes spin–orbit coupling and an added self-consistent orbital-polarization (OP) mechanism for more accurate treatment of magnetism. The first-principles method is coupled to a lattice dynamics scheme that is used to model anharmonic lattice vibrations, namely, Self-Consistent Ab Initio Lattice Dynamics (SCAILD). The methodology can be summarized in the acronym DFT + OP + SCAILD. Upon thermal expansion, γ-U develops non-negligible magnetic moments that are included for the first time in thermodynamic theory. The all-electron DFT approach is shown to model γ-U better than the commonly used pseudopotential method. In addition to CALPHAD, DFT + OP + SCAILD thermodynamic properties are compared with other ab initio and semiempirical modeling and experiments. Our first-principles approach produces Gibbs free energy that is essentially identical to CALPHAD. The DFT + OP + SCAILD heat capacity is close to CALPHAD and most experimental data and is predicted to have a significant thermal dependence due to the electronic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ab initio Molecular Dynamics Assessment of Thermodynamic and Transport Properties in (K,Li)Cl and (K, Na)Cl Molten Salt Mixtures

Molten salt mixtures are integral part of highly important technological applications such as nuclear reactors. However, due to inherent difficulties associated with experiment at high temperatures and the intrinsic complexity of liquid-phase multi-component systems, understanding their properties at a molecular level remains a challenge. Here, we report on an ab initio molecular dynamics investigation on structural, electronic, transport, and thermal properties of two common molten salt mixtures, (K, Li)Cl and (K,Na)Cl, at five different compositions and three temperatures. Most of the properties were found to depend on both composition and temperature. While properties, like atomic charges, show additive behaviors, other properties, such as electrical conductivity, show considerable deviations from additivity. We shall show that the mixing of the molten salt mixtures is mainly driven by entropy, and that the KCl and LiCl mix better than KCl and NaCl. Our computational results are in general consistent with available experimental data. Comparison with available theoretical data is also provided.

Nguyen, Manh Thuong↗
DOE OSTI2025

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ab Initio Modeling on The Thermodynamic and Temperature-Dependent Elastic Properties of Subsystems of The FCC FeNiCoCr Medium Entropy Alloys (MEAs)

The stability, phonon spectra, thermodynamic, and temperature-dependent elasticity of subsystems of the FCC FeNiCoCr MEAs are systematically studied by the ab initio approach. Especially, the quasi-harmonic approximation (QHA) and the innovative Zentropy theory were utilized to predict the thermodynamic properties and elastic properties of FeNi, NiCo, FeNiCo, and FeNiCoCr MEAs with the consideration of magnetic transition. With the ensemble of the partition function based on the multiplicity of each magnetic configuration, the current work successfully reproduced the Curie temperature and the Schottky anomaly of heat capacity of these four MEAs purely based on the ab initio input. Meanwhile, the elastic properties of these alloys at finite temperatures are also successfully predicted with the consideration of magnetic transition. The overall results agree well with the available experimental data and CALPHAD prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tuning the Electronic Properties of Graphane via Hydroxylation: An Ab Initio Study

The thermodynamic stability of hydroxylated graphane, that is, fully sp 3 graphene derivatives coordinated with -H and -OH groups, has been recently demonstrated by ab initio calculations. Within the density functional theory approach, we investigate the electronic property modifications of graphane by progressive hydroxylation, that is, by progressively substituting -H with -OH groups. When 50% of graphane is hydroxylated, the energy bandgap reaches its largest value of 6.68 eV. The electronic affinity of 0.8 eV for graphane can widely change in the 0.28–1.60 eV range depending on the geometric configuration. Hydroxylated graphane has two interfaces with vacuum, hence its electron affinity can be different on each interface with the formation of an intrinsic dipole perpendicular to the monolayer. We envisage the possibility of using hydroxylated graphane allotropes with tunable electronic affinity to serve as interfacial layers in 2D material-based heterojunctions.

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

Ab initio studies on structural and thermodynamic properties of magnetic Fe

The present work systematically investigates the total energy, phonon spectra, and thermodynamic properties of different polymorphs of pure Fe, i.e., FCC, BCC, and HCP, with the ab initio approach, considering various magnetic configurations. In general, the calculated energy vs. volume curves and phonon spectra agree well with previous calculations and the experimental data. In addition, their thermodynamic properties are estimated by the quasiharmonic approximation (QHA). Specifically, a superposition approach based on the latest Zentropy theory was utilized to predict magnetic transition temperatures and thermodynamic properties of pure Fe. With the ensemble of the partition function considering the multiplicity of each magnetic microstate, the current work successfully reproduced the Curie/Néel temperature and the Schottky anomaly of heat capacity in FCC, BCC, and HCP Fe purely based on the ab initio input, which exhibits good agreement with the experimental data and CALPHAD modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chemistry Informed Machine Learning-Based Heat Capacity Prediction of Solid Mixed Oxides

Knowing heat capacity is crucial for modeling temperature changes with the absorption and release of heat and for calculating the thermal energy storage capacity of oxide mixtures with energy applications. The current prediction methods (ab initio simulations, computational thermodynamics, and the Neumann–Kopp rule) are computationally expensive, not fully generalizable, or inaccurate. Machine learning has the potential of being fast, accurate, and generalizable, but it has been scarcely used to predict mixture properties, particularly for mixed oxides. Here, we demonstrate a method for the generalizable prediction of heat capacity of solid oxide pseudobinary mixtures using heat capacity data obtained from computational thermodynamics and descriptors from ab initio databases. Further, models trained through this workflow achieved an error (mean absolute error of 0.43 J mol –1 K –1 ) lower than the uncertainty in differential scanning calorimetry measurements, and the workflow can be extended to predict other properties derived from the Gibbs free energy and for higher-order oxide mixtures.

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

Theoretical study of the electronic, thermodynamic, and thermo-conductive properties of y-LiAlO2 with 6Li isotope substitutions for tritium production

Lithium aluminate has attracted researchers’ interests due to its wide applications. The structural, electronic, optical, lattice phonon thermodynamic, and thermo-conductive properties of ?-LiAlO2 with 6Li-isotope substitution are investigated by density functional theory and lattice phonon dynamics. The calculated results show that ?-LiAlO2 possesses a wide band-gap of 4.63 eV. The valence band is mainly formed by p orbitals of Li, O, and Al. The calculated lattice thermal conductivity of ?-LiAlO2 is found to be in good agreement over a wide range of temperatures with the available experimental data. The calculated dielectric matrix shows anisotropic behavior along z and x (or y) axes. Due to the wide band-gap, the optical conductivity of ?-LiAlO2 is mainly contributed from the electron hopping between valence-conduction bands. Substitution of 7Li with 6Li isotope in ?-LiAlO2 leads to observable differences in the lattice phonon frequencies at higher frequencies (>8THz) and to slight changes in corresponding thermal conductivity and the infrared and Raman spectra which are in good agreement with the measured data. The relevance of these properties for tritium production through 6Li absorbing neutron in nuclear reactors is discussed.

ab initio thermodynamics, ?-LiAlO2 and 6Li-isotope↗
DOE OSTI2023

Electronic structural and lattice thermodynamic properties of MAlO2 and M5AlO4 (M = Li, Na, K) sorbents for CO2 capture applications

Abstract The electronic properties and thermal stabilities of MAlO 2 and M 5 AlO 4 (M = Li, Na, K) are investigated by density functional theory and lattice phonon dynamics. Based on the calculated electronic and lattice thermodynamic properties, their abilities to capture CO 2 as solid sorbents are analyzed. The calculated electronic structural properties of MAlO 2 and M 5 AlO 4 indicate that all these alkali aluminates are semiconductors with a bandgap range of 2.4 ~ 6.4 eV. The 1st valence bands of these alkali aluminates are located 0 ~ − 6 eV under Fermi levels and are mainly contributed by p orbitals of O, s and p orbitals of Al and M. The phonon vibrational frequencies of M 5 AlO 4 spread at a lower frequency range compared to their MAlO 2 phases. With increasing temperature, the calculated phonon free energies of M 5 AlO 4 decrease faster than their corresponding MAlO 2 while their entropies have opposite trends. The reaction 2MAlO 2 + CO 2 = M 2 CO 3 + Al 2 O 3 has higher reaction heat and Gibbs free energy change than those of corresponding reaction 2 / 5 M 5 AlO 4 + CO 2 = M 2 CO 3 + 1 / 5 Al 2 O 3 , which shows the former reaction possesses lower turnover temperature. Among the alkali aluminates studied, the β-NaAlO 2 , lt-KAlO 2 , and γ-LiAlO 2 are better candidates that could be applied for CO 2 capture technologies. Graphical Abstract

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling Pb(II) Adsorption on Mineral Surfaces: Bridging Density Functional Theory and Experiment with Thermodynamic Insights

Despite decades of work on aqueous lead (Pb) adsorption on a-Fe2O3 (hematite) and a-Al2O3 (alumina), gaps between measurements and modeling obscure molecular-level understanding. Achieving well-matched geometries between theory and experimental for mineral-water interfaces is a hurdle, as surface functional group type and distribution must be accounted for in determining mechanisms. Additionally, computational methods that can describe the substrate are often not appropriate to capture aqueous effects. Progress requires focusing on well-studied and relevant systems, such as key facets (001),(012), and (110) of hematite and alumina, and ubiquitous contaminants such as aqueous Pb. In the past, bulk-parametrized bond-valence principles were used to rationalize Pb(II) adsorption trends. These approaches can break down at surfaces, where flexible bonding environments and adsorption-induced surface relaxations play a critical role. Here, we adapt and apply a density functional theory (DFT) and thermodynamics framework, integrating DFT-calculated energies with experimental data and electrochemical principles, to predict Pb(II) adsorption. Our model results capture trends across the full set of surfaces and predict that inner-sphere Pb(II) sorption on (001) alumina varies from unfavorable to weakly favorable across a range of pH conditions. This aligns with experiment insights that Pb(II) interacts at that surface through outer-sphere interactions. Extending to Fe(II) adsorption, we demonstrate a coverage-dependent site preference, potentially explaining disorder in overlayers grown by the oxidative adsorption of Fe(II) on hematite (001).

lead contamination↗
DOE OSTI2022

Thermodynamics of the Metal Carbonates and Bicarbonates of Mn, Co, Ni, Cu, and Zn Relevant to Mineral Energetics

The gas phase heats of formation of ground-state MCO 3 , M(HCO 3 ) 2 , and M(HCO 3 )(OH), where M = Mn, Co, Ni, Cu, and Zn, have been predicted using the correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit using the Feller–Peterson–Dixon (FPD) approach. Cohesive energies of the carbonates were predicted based on the calculated gas phase and experimental solid heats of formation. Coulombic dissociation energies (CDEs) between metal cations and anions show a near-linear correlation with Shannon metal cation atomic radii, yet no correlation is found with the hardness of these cations. The total reaction dissociation energies (TRDEs) of transition metals are higher than their CDEs for the di-bicarbonates, in contrast to those for Mg and Ca based on our prior work. In addition to differences in the energies needed to prepare the transition metal dications, electron donation from the ligands to the 3d orbitals of open-shell transition metal dications from lone pairs of adjacent O atoms also plays a role. No electron donation from the ligands to the fully occupied 3d orbitals of Zn and Cd was found. Decomposition energies for generating MO, CO 2 , and/or H 2 O were calculated. Gas phase metal exchange energies only partially correlate with the electrochemical series for M(s) → M 2+( aq). The FPD heats of formation were used to benchmark a range of density functional theory exchange–correlation functionals, including those commonly used in solid-state mineral calculations. None of the functionals provided chemical accuracy agreement (±1 kcal/mol) with the FPD results. So, the best agreement with the FPD results is predicted for the $\tau$-HCTH functional with an average unsigned error of 8.3 kcal/mol.

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

The α-WC(0001) Surface as a Hydrogen Sponge: A First Principle Study of H 2 Dissociation and Formation of Low and High Coverages

Tungsten carbide (WC) displays a Pt-like behavior in catalysis, applied in hydrogenation processes. Numerous theoretical studies have modeled the behavior and use of adsorbed hydrogen without obtaining a general picture, missing basic links between H 2 dissociation and generation of high surface coverage (Θ H >0.5 ML). Here, in this study, the capability of C- and W-terminations of the α-WC(0001) surface is analyzed to dissociate several H 2 molecules to produce coverages, Θ H , ranging from low to very high values (0.13<Θ H <2.00 ML). Density functional theory and an ab initio atomistic thermodynamic were used to achieve the conditions for H 2 dissociation. The WC–C surface has higher capacity to dissociate H 2 molecules than WC–W. However, both surfaces can reach full surface coverage, Θ H =1 ML, at mild ambient conditions, T=300 K and P=1 atm, and even up to 500 K at low and high pressures. The H-adatoms on WC–W are more labile than on WC–C. The binding of adsorbates is hindered at high Θ H , implying a need to modulate Θ H according to the application. The results give the basis to understand the capabilities of WC-based catalysts in hydrogenation-related reactions, with the advantage of WC being a hydrogen reservoir at mild practical catalytic conditions.

03 NATURAL GAS↗
DOE OSTI2022

Combustion chemistry of alkenes and alkadienes

Alkenes formed, during the refining of crude oil, by cracking the heavier fractions are present in transportation fuels in significant amounts, up to as much as 15–20% in gasoline. Moreover, alkenes are also the major intermediate products of the oxidation of alkanes, which play a significant role in autoignition chemistry. This review has assessed the recent progress in gas-phase detailed kinetic model development for species with C=C double bond, mostly C2–C10 alkenes and 1,3-butadiene. The compiled knowledge on alkene combustion chemistry enabled a better understanding of the influence of the number and the position of the C=C double bond on the chemical kinetics and hence combustion behavior of alkenes in engines. At first, the article gives an extensive overview of fundamental combustion experiments by considering studies of C2–C10 alkenes and 1,3-butadiene in shock tubes, rapid compression machines, laminar flames, and jet-stirred and flow reactors. The value of the data from such experiments is critically discussed. Secondly, this article highlights the important reaction classes involved in alkene oxidation over low-, intermediate- and high-temperature ranges. Combustion chemistry covering C2 to C10 alkenes, with a special emphasis on C2 to C7 isomers is discussed by presenting a large body of experimental and modeling investigations. Detailed chemistry differences between alkene isomers and also between alkenes and alkanes are also addressed. Thirdly, the article presents important reaction pathways for PAH precursor formation in different alkenes. Finally, a summary of the distinguishing features of alkene combustion chemistry and an outlook towards future research in this area are presented. In conclusion, this review is focused on linear and branched chain alkenes, and the chemistry of cyclo-alkenes is not included.

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

Understanding CO 2 adsorption on the surfaces of SrO and its hydroxylated variants Sr(OH) 2 · n H 2 O ( n = 0, 1, 8)

Strontium oxide (SrO) is a promising material for CO 2 capture through a reversible cycle of carbonation and calcination, where SrO reacts with CO 2 to form SrCO 3 and can be regenerated by calcination. In the presence of moisture, SrO forms strontium hydroxide and its hydrates (Sr(OH) 2 ·nH 2 O). Here, this study, which employs density functional theory, investigates the CO 2 adsorption mechanism of these processes on various crystal surfaces, including SrO, Sr(OH) 2 , Sr(OH) 2 ·1H 2 O, and Sr(OH) 2 ·8H 2 O. A significant finding is that the interaction of CO 2 with these surfaces leads to carbonate/bicarbonate formation via electron transfer, with notable differences in CO 2 orientation and bond characteristics between SrO surfaces and its hydroxylated surfaces. To explore the effects of moisture on CO 2 adsorption, H 2 O adsorption on these stable surfaces was investigated. The results showed that H 2 O reacts with the SrO (100) surface to form hydroxyl (OH) groups while it bonds with the surfaces of Sr(OH) 2 ·nH 2 O (n = 0, 1, 8) with hydrogen bonding. A small amount of H 2 O can enhance CO2 adsorption while a large amount of H 2 O could decrease the capability of CO 2 adsorption.

CO2 capture↗
DOE OSTI2022

Synthesis and Thermophysical Property Determination of NaCl-PuCl3 Salts

Currently, a knowledge gap exists in the available data and understanding of thermophysical properties relating to fresh fuel salts, especially those containing plutonium. Determination of this data is necessary for the design and construction of test reactors, as well as the licensing of future commercial molten salt reactors. Thermophysical properties such as melting temperature, salt stability, density, and heat capacity were determined on synthesized eutectic NaCl-PuCl3 and a more sodium rich composition containing 25 mol% PuCl3. These measurements document the baseline properties of the salt as a function of temperature for future experiments on irradiated fuel salt which will provide a holistic perspective on the change of thermophysical properties during reactor operations. It was determined that the NaCl-PuCl3 ingot synthesized for this study contained 63.4 mol% NaCl, 36.3 mol% PuCl3, and was 99.7% pure. Upon heating it was shown that the NaCl- PuCl3 eutectic was stable at temperatures of 800°C. The onset of melting occurred at 541°C and the enthalpy of fusion was determined to be 140.7 ± 8.4 J/g. Specific heat capacity measurements showed a slightly decreasing trend with respect to temperature in the liquid phase ranging from 0.67 to 0.57 J/g·K, with an average value of 0.637 ± 0.03 J/g·K (104 ± 5 J/mol·K) between 500 to 720°C. Three independent trials of the molten NaCl-PuCl3 eutectic salt found the density to be ?(T) = 3.8589 – 9.5342·10-4 T(°C), validated between 500 to 800°C. In addition to salt synthesis and experimentally determining thermodynamic properties, Ab Initio molecular dynamic (AIMD) simulations were used to calculate density and heat capacity values.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗