Search NASA⌕ Search

SEARCH · Search NASA

Results for “THERMODYNAMIC PROPERTY”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Hydrothermal solution calorimetry in acidic aqueous solutions and revisiting the standard partial molal thermodynamic properties of Nd 3+ from 25 to 300 °C

The mobility of rare earth elements (REE) can be predicted in aqueous fluids using geochemical modeling but the accuracy of these models strongly depends on the availability of robust thermodynamic properties for the REE aqueous species. The REE 3+ aqua ions are important in the derivation of the formation constants of all the major REE complexes including the chloride, sulfate, and fluoride species which predominate in many hydrothermal-magmatic systems. However, the thermodynamic properties of the REE 3+ aqua ions are still commonly derived from the Helgeson-Kirkham-Flowers (HKF) equation of state parameters tabulated several decades ago. The standard state thermodynamic properties at reference conditions (25 °C and 1bar) and their extrapolations to high temperature need to be verified, if not revised, based on hydrothermal experiments. In this study, the enthalpy of solution was measured for synthetic Nd hydroxide from 25 to 150 ºC to retrieve the standard partial molal thermodynamic properties of Nd 3+ as a function of temperature. The experiments were conducted in aqueous perchloric acid based solutions with starting pH of 2 and varying ionic strength (0.01 to 0.09 mol/kg NaClO 4 ). The standard partial molal enthalpy of formation (Δ f H°) of Nd 3+ derived from the experimental study displays differences of up to 10 kJ/mol compared to the enthalpy values derived from the HKF equation of state in the studied temperature range. These inaccuracies are resolved by adjusting the standard partial molal Gibbs energy of formation (Δ f G°) of Nd 3+ at 25 °C and 1 bar from -672.0 to -679.7 kJ/mol. The heat capacity function (C p °) derived between 25 and 150 ºC can be described by: C p ° = a 0 + a 1 ·T + a 2 ·T -2 , with a 0 = 1256, a 1 = -2.68, a 2 = -55.56·10 6 and T in Kelvin. A set of recommended thermodynamic properties is provided for the Nd 3+ aqua ions and corrections are provided for the chloride and fluoride species to remain internally consistent with the experimentally derived properties. These results allow predicting accurately the solubility of monazite between 25 and 300 ºC. Before these corrections, the properties for the Nd 3+ aqua ions derived from the HKF parameters resulted in up to ~1.5 orders of magnitude lower monazite solubility than determined experimentally. Here, a revision of the REE +3 aqua ions properties is necessary to accurately predict the mobility of REE in hydrothermal acidic solutions.

58 GEOSCIENCES↗
DOE OSTI2025

Thermodynamic properties of “near-perfect” gas

When applied thermodynamics requires simple representations of thermodynamic quantities, one often finds quantities such as pressure and the thermal energy density expressed as a monomial linear product of density and temperature power laws. This is a simple generalization of the product of density and temperature used for the perfect gas that admits a broader range of thermodynamic behavior into analytic fluid calculations or provides an analytic form that is readily fit to tabulated equation-of-state data for an arbitrary material. This paper reviews the thermodynamic properties of this generalized perfect-gas model, treating it as a class of “near-perfect-gas” models that are defined and unified here, for the first time, through a shared Helmholtz free energy. This generalization from perfect to near-perfect-gas models preserves some perfect-gas properties (e.g., constant adiabatic exponents) but not others (e.g., constant specific heats) and is constrained by thermodynamic consistency. It is important to be aware of the properties of this class of near-perfect-gas models when using them as surrogates for real nonperfect materials that may have properties that cannot be captured by near-perfect-gas models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The universal thermodynamic properties of extremely compact objects

An extremely compact object (ECO) is defined as a quantum object without horizon, whose radius is just a small distance s outside its Schwarzschild radius. We show that any ECO of mass M in d + 1 dimensions with s << (M/m p ) 2/(d-2)(d+1) lp must have (at leading order) the same thermodynamic properties—temperature, entropy and radiation rates—as the corresponding semiclassical black hole of mass M. An essential aspect of the argument involves showing that the Tolman–Oppenheimer–Volkoff equation has no consistent solution in the region just outside the ECO surface, unless this region is filled with radiation at the (appropriately blueshifted) Hawking temperature. In string theory it has been found that black hole microstates are fuzzballs—objects with no horizon—which are expected to have a radius that is only a little larger than the horizon radius. Thus the arguments of this paper provide a nice closure to the fuzzball paradigm: the absence of a horizon removes the information paradox, and the thermodynamic properties of the semiclassical hole are nonetheless recovered to an excellent approximation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Computing chemical potentials with machine-learning-accelerated simulations to accurately predict thermodynamic properties of molten salts

The successful design and deployment of next-generation nuclear technologies heavily rely on thermodynamic data for relevant molten salt systems. However, the lack of accurate force fields and efficient methods has limited the quality of thermodynamic predictions from atomistic simulations. Here we propose an efficient free energy framework for computing chemical potentials, which is the central free energy quantity behind many thermodynamic properties. We accelerate our simulations without sacrificing accuracy by using machine learning interatomic potentials trained on density functional theory (DFT) data. Using lithium chloride as our model system, we compute chemical potentials with DFT-accuracy for solid and liquid phases by transmuting ions into noninteracting particles. Notably, in the liquid phase, we demonstrate consistency whether we transmute one ion pair or the entire system into ideal gas particles. By locating the temperature where the chemical potential of solid and liquid phases cross, we predict a melting point of 880 ± 18 K for lithium chloride, which is remarkably close to the experimental value of 883 K. With this successful demonstration, we lay the foundation for high-throughput thermodynamic predictions of many properties that can be derived from the chemical potentials of the minority and majority components in molten salts.

Gibson, Luke D. [Oak Ridge National Laboratory (OR↗
DOE OSTI2025

Equation of State for the Thermodynamic Properties of Trans-1,2-dichloroethene [R-1130(E)]

We present an empirical equation of state in terms of the Helmholtz energy for trans-1,2-dichloroethene [R-1130(E)]. The range of validity is from the triple-point temperature, 223.31 K to 525 K with pressures up to 30 MPa. It may be used to calculate all thermodynamic properties in the fluid phase, including liquid, gas, and supercritical regions. Comparisons are given with existing literature data and estimated uncertainties are provided. In addition, checks were made for correct extrapolation behavior so that the equation behaves in a physically realistic manner when used outside of its range of validity, enabling its use in mixture models. The estimated uncertainties (at a k = 2 or 95 % level of confidence) are based on comparisons with critically assessed data and are 0.25 % for vapor pressure for temperatures in the range 300 K < T < 454 K, rising to 1.5 % as the temperature decreases from 300 K to 265 K. For density in the liquid phase the estimated uncertainty is 0.14 % for temperatures 270 K < T < 410 K and for pressures up to 30 MPa. For the vapor phase the estimated uncertainty in density is 3 %. The uncertainty for liquid-phase heat capacity is 1 % at atmospheric pressure over the temperature range 268 K < T < 309 K, and the uncertainty for the speed of sound in the liquid phase is 0.25 % for temperatures 230 K < T < 420 K and for pressures up to 30 MPa. The uncertainties are larger outside of these specified ranges and in the critical region.

1,2-Dichloroethene↗
DOE OSTI2025

A noncollinear density functional theory ansatz for the phononic and thermodynamic properties of α -Pu

Plutonium's phase diagram is host to complex structures and interactions that make the description of its ground state properties elusive. Using all-electron density functional theory calculations, we study the thermodynamic properties of α-Pu. To do this, we build on recent work in the literature by introducing a novel noncollinear magnetic ansatz for α-Pu's ground state. We show that the noncollinear ansatz improves the description of the experimental phonon density of states, heat capacity, and thermal expansion over previous results obtained with a collinear ansatz. Here, despite lacking a bond-equivalent picture of the bonding, these new results for α-Pu, along with recent results on δ-Pu, demonstrate the efficacy of noncollinear ansatzes for the description of plutonium.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Consistent thermodynamic properties for alicyclic components of jet fuels: Experimental data, estimation methods, and homologous series trends

Alkylcycloalkanes represent a significant fraction of jet fuel components. An evaluation of their thermodynamic properties, enthalpies of formation in liquid and gas phases and enthalpies of vaporization, was conducted. A combination of available experimental data, up-to-date group-contribution methods, high-level quantum-chemical calculations, and homologous series trends was used to identify outliers and to recommend the most reliable values. The group-contribution approach was found to work well for the enthalpies of vaporization. Its performance for the enthalpies of formation in the liquid and gas phases was found to be substantially less effective, especially considering notable differences in this property among stereoisomers. Computationally affordable high-level ab initio results and homologous series trend analysis appeared more reliable. In conclusion, the recommended property values for 212 individual compounds and their isomeric mixtures were provided.

09 BIOMASS FUELS↗
DOE OSTI2024

First-principles investigation of elastic, vibrational, and thermodynamic properties of kagome metals CsM 3 Te 5 (M = Ti, Zr, Hf)

Kagome metals are a unique class of quantum materials characterized by their distinct atomic lattice arrangement, featuring interlocking triangles and expansive hexagonal voids. These lattice structures impart exotic properties, including superconductivity, interaction-driven topological many-body phenomena, and magnetism, among others. The kagome metal CsM 3 ⁢Te 5 (where M = Ti, Zr, or Hf) exhibits both superconductivity and nontrivial topological electronic properties, offering a promising platform for exploring topological superconductivity. This study employs first-principles density functional theory calculations to systematically analyze the elastic, mechanical, vibrational, thermodynamic, and electronic properties of CsM 3 ⁢Te 5 (M = Ti, Zr, Hf). Our calculations reveal that the studied compounds—CsTi 3 ⁢Te 5 , CsZr 3 ⁢Te 5 , and CsHf 3 ⁢Te 5 —are ductile metals with elastic properties akin to the hexagonal Bi and Sb, with average elastic constants, including a bulk modulus of 27 GPa, a shear modulus of 11 GPa, and Young's modulus of 29 GPa. We observe peculiar dispersionless, flat, phonon branches in the vibrational spectra of these metals. Additionally, we thoroughly analyze the symmetries of the zone-center phonon eigenvectors and predict vibrational fingerprints of the Raman- and infrared-active phonon modes. The analysis of thermodynamic properties reveals the Einstein temperature for CsTi 3 ⁢Te 5 , CsZr 3 ⁢Te 5 , and CsHf 3 ⁢Te 5 to be 66, 54, and 53 K, respectively. Our orbital-decomposed electronic structure calculations reveal significant in-plane steric interactions and multiple Dirac band crossings near the Fermi level. We further investigate the role of spin-orbit coupling effect on the studied properties. Furthermore, this theoretical investigation sheds light on the intriguing quantum behavior of kagome metals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2025

Dual modulation of the anion-driven thermodynamic properties of aqueous choline halide-based deep eutectic solvents

Deep eutectic solvents (DESs) are considered tunable solvents because their specific properties can be achieved based on the choice of components and their relative concentrations in a mixture. In this article, we investigate the influence of the variation in halide ions (F − , Cl − , Br − , I − ) of choline salts used on the thermodynamic and physicochemical properties of choline halide-based DESs. Our findings show that the density of choline halide-based DESs decreases nonlinearly with an increasing mole fraction of water, following a trend based on the size of the halides, with choline iodide showing the highest density. Temperature-dependent density data reveal that the thermal expansion coefficient decreases slightly with increasing water content, indicating a more stable volume at a higher mole fraction of water. The excess molar volume (V E ) of the DES mixtures exhibits complex behavior depending on the choline halide used, with both negative and positive V E values observed across different water mole fractions. These variations are linked to the hydrogen bonding interactions between the DES components and water molecules. In addition, viscosities decrease with increasing water content, suggesting the disruption of hydrogen bonding networks and enhanced mobility of the ions, which contributes to the observed increase in conductivity. The excess molar Gibbs energies, enthalpies, and entropies of activation have also been determined.

Choline halide↗
DOE OSTI2025

Thermodynamic properties of ZrCl 4 with LiCl, NaCl, KCl, CsCl, MgCl 2 , and UCl 3 for molten salt reactor applications

Here, a set of self-consistent Gibbs energy functions has been developed to describe the thermochemical behavior of the major reactive fission product ZrCl 4 with key components in chloride molten salt reactors (MSRs): LiCl–ZrCl 4 , NaCl–ZrCl 4 , KCl–ZrCl 4 , CsCl–ZrCl 4 , MgCl 2 –ZrCl 4 and UCl 3 –ZrCl 4 . Low ZrCl 4 concentration phase equilibria most relevant to MSR applications have been confirmed via differential scanning calorimetry for NaCl–ZrCl 4 , and X-ray diffraction analysis performed on equilibrated samples of NaCl–ZrCl 4 , KCl–ZrCl 4 , MgCl 2 –ZrCl 4 , and UCl 3 –ZrCl 4 . Within the framework of the modified quasi-chemical model in the quadruplet approximation, extrapolations of pseudo-binary models were generated to represent KCl–MgCl 2 –ZrCl 4 , KCl–NaCl–ZrCl 4 , and MgCl 2 –NaCl–ZrCl 4 , which show agreement with phase equilibria data where available without the use of ternary interaction parameters. The optimized thermodynamic descriptions for these systems and others are available in the open-source compendium Molten Salt Thermal Properties Database – Thermochemical (MSTDB–TC).

CALPHAD↗
DOE OSTI2025

Impact of chemical ordering on thermodynamic properties of point defects and Xe substitutional in U-10Mo

The accurate knowledge of defect energetics is critical to understanding the aging and irradiation behavior of U-10Mo nuclear fuel, which is selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). Furthermore, using hybrid molecular dynamics and Monte Carlo (MDMC) simulation, we studied the impact of chemical ordering on the formation energies of vacancies, interstitials, and the solution energy of the Xe substitutionals. Instead of forming a random solid solution (RSS), substantial short-range-order (SRO) develops in U-10Mo, particularly at low temperatures. Mo atoms are found to repel each other and prefer U-rich local atomic environments within the 1st nearest neighbor (1NN) cutoff. Compared to the case of a RSS, the state with equilibrated Mo ordering shifts the distributions of vacancy and interstitial formation energies due to the dependence of defect energies on the local atomic environment, without a clear effect on Xe solution energy. In the operation temperature range (100–250 °C) of U-Mo fuels, neglecting SRO can lead to an inaccurate estimate of thermal equilibrium point defect concentrations by over an order of magnitude and incorrectly predict the preference among different types of dumbbells, highlighting the critical importance of accounting for the impact of chemical ordering for accurate atomistic calculations of defect properties.

Molecular dynamics and Monte Carlo (MDMC)↗
DOE OSTI2026

The solubility of NdPO 4 and DyPO 4 and stability of Nd and Dy chloride and hydroxyl complexes as a function of pH and salinity to 450 °C

The mobility of rare earth elements (REE) in geological systems is often controlled by the stability of monazite and xenotime, which are important hosts for the light and heavy REE. These REE phosphates constitute important resources and provide information on ore formation conditions and timing due to their uses as geothermometers and geochronometers. While the thermodynamic properties for these minerals are well-established up to high temperature and pressure, the properties for the aqueous REE complexes are not well constrained with limited information up to 300 °C. In this study, the solubility of synthetic NdPO 4 and DyPO 4 endmembers were measured in hydrothermal sub- to supercritical NaCl-bearing aqueous solutions from 350 at saturated water vapor pressure to 450 °C at 700 bar. The speciation of Nd and Dy was determined in acidic to alkaline solutions (pH 25 °C values from 2 to 10), and indicates that the hydroxyl species REE(OH) 2 + and REE(OH) 3 0 predominate at low salinity (0.01 mol/kg NaCl) with some contribution of REE chloride species REECl 2+ and REECl 2 + in acidic fluids. The REE phosphate solubility measured in these experiments is up to two orders of magnitude higher than predicted using existing thermodynamic properties from literature for aqueous species extrapolated from lower temperature data to supercritical conditions. To address this discrepancy, the thermodynamic properties of the REE aqueous species were optimized using GEMSFITS to derive the formation constants for the REE chloride (β Cl ) and hydroxyl (β OH ) complexes in the studied temperature range. This study highlights the importance of the hydroxyl species REE(OH) 2 + and REE(OH) 3 0 over a wide range of pH and temperature. As a result, the revised thermodynamic properties for Dy and Nd chloride and hydroxyl species more accurately predict the solubility of NdPO 4 and DyPO 4 and provide insight into the speciation and mobility of the light and heavy REE in hydrothermal supercritical fluids in the crust.

58 GEOSCIENCES↗
DOE OSTI2025

Profiles of Radiative Fluxes at ENA

Profiles of radiative fluxes observed at the Atmospheric Radiation Measurement (ARM)’s Eastern North Atlantic (ENA) observatory along with the ancillary measurements are reported. The below-cloud drizzle properties were derived by combining the data from the ceilometer and Ka-band ARM Zenith Radar (KAZR) following the technique explained by Ghate et al. (2021 JAMC). The cloud and drizzle water path values were derived from the brightness temperatures reported by the microwave radiometer following the technique of Cadeddu et al. (2020 AMT). The cloud water path was then scaled to the KAZR-reported radar reflectivity to calculate profiles of liquid water content (LWC). Following the analysis from Ghate et al. (2023 JGR), cloud droplet effective radius was calculated using the number concentration value of 100 cm-3. The cloud properties, along with the thermodynamic properties, served as an input to the Rapid Radiative Transfer Model (RRTM) to yield profiles of radiative fluxes at a 1-minute temporal and 50-m vertical resolution. The fluxes were then averaged to hourly temporal resolution for analysis. In Mitra et al. (2025 JClim), the calculated profiles were compared against those derived from the satellite measurements (SYN1deg). Flux profiles from the SYN1deg and the thermodynamic and cloud properties used for deriving them are also reported here. Both all-sky and clear-sky radiative flux profiles were calculated. Due to the large data volume, the surface and top-of-atmosphere (TOA) radiative fluxes for the six-year period, and the hourly profiles of the radiative fluxes for January 2018, are submitted here. Full profiles of radiative fluxes calculated from the thermodynamic and cloud properties measured at the ENA site at 1-minute temporal and 50-m vertical resolution for a six-year period are available from the authors. Six files here correspond to the following data: 1_ENARAD_CERES_with_cld_amount_timeseries.nc: Time-series of hourly values of RRTM-simulated values of upwelling and downwelling fluxes at the surface and TOA, observed boundary-layer cloud fractions, and upwelling and downwelling fluxes from the SYN1deg from July 2015 to January 2022. 2_CERES_2018_at_CERES_levels.nc: SYN1deg radiative fluxes at six levels for the year 2018. 3_ENARad_2018_at_CERES_levels.nc: RRTM calculated fluxes at the SYN1deg vertical levels for the year 2018. 4_ENARad_rrtminputs_hourly_201801.nc: Thermodynamic and cloud properties used as an input to the RRTM for January 2018. 5_CERES_inputs_hourly_201801.nc: Thermodynamic and cloud properties utilized by SYN1deg algorithm for January 2018. 6_ENARAD_hourly_201801.nc: Full profiles of hourly averaged radiative fluxes from the RRTM simulations for January 2018.

Atmosphere↗
DOE OSTI2026

The solubility and speciation of REE phosphate endmembers (CePO 4 and YPO 4 ) in Cl-rich aqueous fluids from 350 to 450 °C and implications for natural systems

The rare earth elements (REE) are important metals used increasingly in advanced technologies. Within the crust, the elements Ce and Y are commonly more abundant compared to other lanthanides and comprise important end-member constituents of REE-bearing minerals. Specifically, Ce is part of the light (L) REE which have larger ionic radii than the heavy (H) REE, which are grouped together with Y. These differences in ionic radius can lead to important physico-chemical trends within the lanthanide group. Despite a recent increase in experimental and thermodynamic data for the REE at high temperature and pressure, there is still a significant lack of these data at supercritical conditions. In this study we conducted batch-type experiments to measure the solubility of REE phosphates (CePO 4 and YPO 4 ) at varying starting pH (1.5–10), and salinity (0.01–1.4 mol/kg NaCl) at 350 and P sat , and from 400 to 450 °C at 700 bar. Results show that the solubility of Ce (33–0.14 ppb) is generally higher than Y (13–0.13 ppb) and that Ce complexes more strongly with both chloride and hydroxyl ligands compared to Y. The solubilities of both REE phosphates are highly pH-dependent and, to a lesser extent, depend on salinity at the studied conditions. The solubility data from this study were implemented into the GEMSFITS program to optimize the thermodynamic properties of Ce and Y hydroxyl and chloride species. The updated standard partial molal Gibbs energies of formation (Δ f G 0 T,P ) are used within the experimental temperature and pressure range to accurately predict the CePO 4 and YPO 4 solubility and Ce and Y speciation behavior. Based on the updated thermodynamic properties we also provide formation constants (log β n Cl,OH ) for Ce and Y hydroxyl and chloride species. Updated thermodynamic properties are applied to model REE-apatite dissolution and REE mobility based on the Pea Ridge iron oxide apatite deposit in Missouri, USA. The apatite dissolution model replicates natural observations including the replacement of monazite and xenotime after apatite and is an example of the utility of the new thermodynamic constants applied to supercritical crustal fluids. Furthermore, the findings of this study advance the predictive capabilities of geochemical models, our understanding of the behavior of individual REE, and permit modeling the overarching fractionation trends between LREE and HREE in supercritical crustal fluids.

58 GEOSCIENCES↗