SEARCH · Engineering Papers

Results for “Enthalpy”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

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 OSTI2024

Low-temperature formation of Ti2AlN during post-deposition annealing of reactive multilayer systems

M n+1 AXn-phase Ti 2 AlN thin-films were synthesized using reactive sputtering-based methods involving the deposition of single-layer TiAlN, and Ti/AlN and TiN/TiAl multilayers of various modulation periods at ambient temperature and subsequent annealing at elevated temperatures. Ex situ and in situ x-ray diffraction measurements were used to characterize the Ti 2 AlN formation temperature and phase fraction. During annealing, Ti/AlN multilayers yielded Ti 2 AlN at a significantly lower in situ temperature of 650 °C compared to TiN/TiAl multilayers or single-layer TiAlN (750 °C). The results suggest a reactive multilayer mechanism whereby distinct Ti and AlN layers react readily to release exothermic energy resulting in lower phase transition temperatures compared to TiN and TiAl layers or mixed TiAlN. With a modulation period of 5 nm, however, Ti/AlN multilayers yielded Ti 2 AlN at a higher temperature of 750 °C, indicating a disruption of the reactive multilayer mechanism due to a higher fraction of low-enthalpy interfacial TiAlN within the film.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Polymorphism and phase transitions in layered uranium(VI) hydroxides: Ab initio lattice dynamics simulations of UO 2 (OH) 2

The phase transitions and thermodynamics of stoichiometric α-, β-, and γ-UO 2 (OH) 2 polymorphs are investigated using density functional perturbation theory. The pressure-induced β(Pbca) → α(Cmca) phase transition is reproduced by calculations, with a volume reduction of ΔV/V = -14.7% similar to experiment. Consistent with observation, a temperature-driven γ(P2 1 /c) → β(Pbca) phase transition is predicted near 533 K. At 298.15K, the computed standard molar heat capacity of α-UO 2 (OH) 2 is C p 0 = 112.1 J mol -1 K -1 , only 1.6% smaller than the value of C p 0 = 113.96 ± 0.12 J mol -1 K -1 measured by calorimetry. C p 0 = 112.4 and 104.8 J mol -1 K -1 are predicted for the β- and γ-UO 2 (OH) 2 polymorphs, respectively. The calculated molar enthalpy and Gibbs energy functions of the α-, β-, and γ-UO 2 (OH) 2 polymorphs are also reported.

74 ATOMIC AND MOLECULAR PHYSICS↗

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 OSTI2022

On the mechanisms of ion adsorption to aqueous interfaces: air-water vs. oil-water

The adsorption of ions to water-hydrophobe interfaces influences a wide range of phenomena, including chemical reaction rates, ion transport across biological membranes, and electrochemical and many catalytic processes; hence, developing a detailed understanding of the behavior of ions at water-hydrophobe interfaces is of central interest. Here, we characterize the adsorption of the chaotropic thiocyanate anion (SCN − ) to two prototypical liquid hydrophobic surfaces, water-toluene and water-decane, by surface-sensitive nonlinear spectroscopy and compare the results against our previous studies of SCN − adsorption to the air-water interface. For these systems, we observe no spectral shift in the charge transfer to solvent spectrum of SCN − , and the Gibb’s free energies of adsorption for these three different interfaces all agree within error. We employed molecular dynamics simulations to develop a molecular-level understanding of the adsorption mechanism and found that the adsorption for SCN − to both water-toluene and water-decane interfaces is driven by an increase in entropy, with very little enthalpic contribution. This is a qualitatively different mechanism than reported for SCN − adsorption to the air-water and graphene-water interfaces, wherein a favorable enthalpy change was the main driving force, against an unfavorable entropy change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Thermodynamic properties and enhancement of diamagnetism in nitrogen doped lutetium hydride synthesized at high pressure

Nitrogen doped lutetium hydride has drawn global attention in the pursuit of room-temperature superconductivity near ambient pressure and temperature. However, variable synthesis techniques and uncertainty surrounding nitrogen concentration have contributed to extensive debate within the scientific community about this material and its properties. We used a solid-state approach to synthesize nitrogen doped lutetium hydride at high pressure and temperature (HPT) and analyzed the residual starting materials to determine its nitrogen content. High temperature oxide melt solution calorimetry determined the formation enthalpy of LuH 1.96 N 0.02 (LHN) from LuH 2 and LuN to be −28.4 ± 11.4 kJ/mol. Magnetic measurements indicated diamagnetism which increased with nitrogen content. Ambient pressure conductivity measurements observed metallic behavior from 5 to 350 K, and the constant and parabolic magnetoresistance changed with increasing temperature. High pressure conductivity measurements revealed that LHN does not exhibit superconductivity up to 26.6 GPa. We compressed LHN in a diamond anvil cell to 13.7 GPa and measured the Raman signal at each step, with no evidence of any phase transition. Despite the absence of superconductivity, a color change from blue to purple to red was observed with increasing pressure. Thus, our findings confirm the thermodynamic stability of LHN, do not support superconductivity, and provide insights into the origins of its diamagnetism.

Science & Technology - Other Topics↗

DOE OSTI2024

Medium-density amorphous ice unveils shear rate as a new dimension in water’s phase diagram

Recent experiments revealed a new amorphous ice phase, medium-density amorphous ice (MDA), formed by ball-milling ice I h at 77 K [Rosu-Finsen et al., Science 379, 474–478 (2023)]. MDA has density between that of low-density amorphous (LDA) and high-density amorphous (HDA) ices, adding to the complexity of water’s phase diagram, known for its glass polyamorphism and two-state thermodynamics. The nature of MDA and its relation to other amorphous ices and liquid water remain unsolved. Here, we use molecular simulations under controlled pressure and shear rate at 77 K to produce and investigate MDA. Here. we find that MDA formed at constant shear rate is a steady-state nonequilibrium shear-driven amorphous ice (SDA), that can be produced by shearing ice I h , LDA, or HDA. Our results suggest that MDA could be obtained by ball-milling water glasses without crystallization interference. Increasing the shear rate at ambient pressure produces SDAs with densities ranging from LDA to HDA, revealing shear rate as a new thermodynamic variable in the nonequilibrium phase diagram of water. Indeed, shearing provides access to amorphous states inaccessible by controlling pressure and temperature alone. SDAs produced with shearing rates as high as 10 6 s −1 sample the same region of the potential energy landscape than hyperquenched glasses with identical density, pressure, and temperature. Intriguingly, SDAs obtained by shearing at ~10 8 s −1 have density, enthalpy, and structure indistinguishable from those of water “instantaneously” quenched from room temperature to 77 K over 10 ps, making them good approximants for the “true glass” of ambient liquid water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Thermodynamics of calcium binding to heparin: Implications of solvation and water structuring for polysaccharide biofunctions

Heparan sulfates are found in all animal tissues and have essential roles in living systems. This family of biomacromolecules modulates binding to calcium ions (Ca 2+ ) in low free energy reactions that influence biochemical processes from cell signaling and anticoagulant efficacy to biomineralization. Despite their ubiquity, the thermodynamic basis for how heparans and similarly functionalized biomolecules regulate Ca 2+ interactions is not yet established. Using heparosan (Control) and heparins with different positions of sulfate groups, we quantify how SO 3 − and COO − content and SO 3 − position modulate Ca 2+ binding by isothermal titration calorimetry. The free energy of all heparin-Ca2+ interactions (ΔG rxn ) is dominated by entropic contributions due to favorable water release from polar, hydrophilic groups. Heparin with both sulfate esters (O-SO 3 − ) and sulfamides (N-SO 3 − ) has the strongest binding to Ca 2+ compared to heparosan and to heparin with only O-SO 3 − groups (~3X). By linking Ca 2+ binding thermodynamics to measurements of the interfacial energy for calcite (CaCO 3 ) crystallization onto polysaccharides, we show molecule-specific differences in nucleation rate can be explained by differences in water structuring during Ca 2+ interactions. A large entropic term (-TΔS rxn ) upon Ca 2+ –polysaccharide binding correlates with high interfacial energy to CaCO 3 nucleation. Combining our measurements with literature values indicates many Ca 2+ –polysaccharide interactions have a shared thermodynamic signature. The resulting enthalpy–entropy compensation relationship suggests these interactions are generally dominated by water restructuring involving few conformational changes, distinct from Ca 2+ –protein binding. Our findings quantify the thermodynamic origins of heparin-specific interactions with Ca 2+ and demonstrate the contributions of solvation and functional group position during biomacromolecule-mediated ion regulation.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations

Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, the formation of ice XI is kinetically hindered, and experimentally it is obtained by doping with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping with a bias potential that enhances the formation of ionic defects in molecular dynamics simulations. The recombination of the ions thus formed proceeds through fast migration of the hydroxide along hydrogen bond loops, providing a physical and expedite way to change the proton configuration. Here, a key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modelling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter that distinguishes the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of the calculation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Effects of cyclic plasma heating on surface damage of micro-porous tungsten

The operating temperature window of solid tungsten (W) is dictated by its Ductile to Brittle Transition Temperature (DBTT) and re-crystallization temperature; roughly between 300-1300 °C. The brittleness of W at lower temperatures is exasperated when it undergoes recrystallization. We investigate here the thermal shock resistance of micro-porous W as a meta-material fabricated in a 3D open-cell network structure. We present experimental results for the effects of cyclic high-enthalpy arc-jet plasma on surface damage in three testing categories. Observed damage includes surface ablation of asperities, melting and solidification of W-fuzz on samples that have been exposed to a prior helium plasma, and micro-cracks at ligament triple junctions. Scanning electron microscope (SEM) observations show more micro-cracking on 54% and 23% foams compared to the 43% ones. In all tested samples, thermal expansion/contraction displacements were accommodated by ligament rotation and a network of micro-cracks. Furthermore, no large through-thickness crack were observed.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Improving tropical cyclone rapid intensification forecasts with satellite measurements of sea surface salinity and calibrated machine learning

Forecasting rapid intensification (RI) of tropical cyclones (TC) is a mission known for large errors. One under-researched factor that affects TC intensification is salinity, which is important for density stratification in certain ocean regions and can affect the surface enthalpy flux under a strengthening hurricane. To investigate the impact and efficacy of using salinity information in state-of-the-art forecasting, we use a statistical model consisting of a variety of machine learning (ML) methods. For salinity data, we use satellite measurements of pre-storm sea surface salinity (SSS) as a proxy for the salinity stratification. We train and test the model on various ocean basins, including the Atlantic, eastern North Pacific and western North Pacific. A calibrator is trained on top of the ML models to correct and enhance probability forecasts. The calibrator significantly improves probability forecasts relative to recent works. The ML model performance is improved with the addition of SSS in the Eastern North Pacific, western North Pacific, and the Caribbean subregion of the North Atlantic, and the overall model performance is better than previous studies. SSS decreases model skill for a model trained on the full Atlantic basin. In the Indian Ocean, SSS is also notably correlated with RI occurrence, but the TC samples are not sufficient to train ML models.

hurricane↗

DOE OSTI2025

The search for high-entropy fuel-cell catalysts using disorder descriptors

The transition to a hydrogen economy depends on efficient, affordable catalysts for fuel cells. Platinum—the industry standard for fuel-cell electrodes—is costly and scarce, highlighting the need for practical alternatives. High-entropy alloys offer vast compositional diversity and tunable properties that can mitigate these issues, yet their chemical complexity and configurational disorder have hindered rational discovery. Here, we introduce a data-driven framework that couples machine learning with first-principles disorder descriptors—including the entropy forming ability, disordered enthalpy-entropy descriptor, and electronic-structure similarity metrics to platinum—to predict alloy synthesizability and catalytic performance. These descriptors are applied for the first time in the context of fuel-cell catalyst discovery. The workflow rapidly screens more than 20 000 compositions and identifies several platinum-free candidates that are economically viable, readily scalable, and exhibit promising predicted activity. These results demonstrate that disorder descriptors are reliably predicted by machine learning models and can be effectively integrated into materials-discovery pipelines, accelerating innovation across complex compositional spaces.

fuel-cell catalysts↗

DOE OSTI2026

The origin of metallic conductivity in Pt 3 O 4 : a first principles study

The platinum oxide Pt 3 O 4 exhibits metallic conductivity even though it contains square-planar PtO4 units, which in related oxides such as PtO are usually associated with insulating behavior. To identify the electronic origin of this anomalous metallicity, we performed a comprehensive first-principles study using the PBE and r 2 SCAN functionals together with Hubbard U corrections and spin-orbit coupling (SOC). Structural benchmarks show that r 2 SCAN with SOC and a moderate U value (<4 eV) reproduces the experimental lattice constants and formation enthalpy, whereas larger U values (~8 eV) destabilize the cubic structure. Across all functionals and U values considered in this work, Pt 3 O 4 remains metallic. Analyses of the projected density of states, band structures, charge-density isosurfaces, and bonding characteristics demonstrate that the dominant contribution to the metallic character originates from delocalized Pt–O–Pt hybridized antibonding states at the Fermi level. Direct Pt–Pt interactions are present but contribute less strongly to the conductivity. Bader charge analysis reveals only weak Pt charge disproportionation, consistent with mixed Pt II /Pt III character, and a small charge-transfer energy that prevents localization of the Pt 5d electrons even at elevated U. In contrast, PtO develops a Mott or charge-transfer gap under modest U despite having the same PtO 4 coordination environment. These findings demonstrate that persistent Pt–O–Pt covalency is the primary driver of metallicity in Pt 3 O 4 and support the view that this phase can remain conductive under oxygen reduction and oxygen evolution reaction conditions in fuel cell and electrolyzer environments.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Multi-task graph neural networks for simultaneous prediction of global and atomic properties in ferromagnetic systems *

Abstract We introduce a multi-tasking graph convolutional neural network, HydraGNN, to simultaneously predict both global and atomic physical properties and demonstrate with ferromagnetic materials. We train HydraGNN on an open-source ab initio density functional theory (DFT) dataset for iron-platinum with a fixed body centered tetragonal lattice structure and fixed volume to simultaneously predict the mixing enthalpy (a global feature of the system), the atomic charge transfer, and the atomic magnetic moment across configurations that span the entire compositional range. By taking advantage of underlying physical correlations between material properties, multi-task learning (MTL) with HydraGNN provides effective training even with modest amounts of data. Moreover, this is achieved with just one architecture instead of three, as required by single-task learning (STL). The first convolutional layers of the HydraGNN architecture are shared by all learning tasks and extract features common to all material properties. The following layers discriminate the features of the different properties, the results of which are fed to the separate heads of the final layer to produce predictions. Numerical results show that HydraGNN effectively captures the relation between the configurational entropy and the material properties over the entire compositional range. Overall, the accuracy of simultaneous MTL predictions is comparable to the accuracy of the STL predictions. In addition, the computational cost of training HydraGNN for MTL is much lower than the original DFT calculations and also lower than training separate STL models for each property.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Direct measurement of electron heating in electron-only reconnection in a laboratory mini-magnetosphere

We report on the experimental observation of electron heating in electron-only magnetic reconnection in laser-driven laboratory mini-magnetosphere on the Large Plasma Device (LAPD) at the University of California, Los Angeles. In this experiment, a fast-flowing plasma impacts a pulsed magnetic dipole embedded within LAPD’s magnetized ambient plasma, creating an ion-scale magnetosphere and driving electron-only magnetic reconnection between the background and dipole field lines. The electron velocity distribution is measured across the reconnection region using noncollective Thomson scattering, enabling determination of electron temperature and density. Significant electron heating is observed in the electron diffusion region, increasing from an initial temperature of 1.8 to 9.5 eV, corresponding to a 40% local conversion of Poynting flux into electron enthalpy flux. Particle-in-cell simulations that provide insights into the heating mechanisms are also presented.

Rovige, Lucas [Department of Physics and Astronomy↗

DOE OSTI2020

Polymorphous nature of cubic halide perovskites

Many common crystal structures can be described by a single (or very few) repeated structural motif (“monomorphous structures”) such as octahedron in cubic halide perovskites. Interestingly, recent accumulated evidence suggests that electronic structure calculations based on such macroscopically averaged monomorphous cubic ($\textit{Pm} – 3m$) halide perovskites obtained from x-ray diffraction, show intriguing deviations from experiment. These include systematically too small band gaps, dielectric constants dominated by the electronic, negative mixing enthalpy of alloys, and significant deviations from the measured pair distribution function. We show here that a minimization of the systems $\textit{T}$ = 0 internal energy via density functional theory reveals a distribution of different low-symmetry local motifs, including tilting, rotations, and B-atom displacements (“polymorphous networks”). This is found only if one allows for larger-than-minimal cell size that does not geometrically exclude low symmetry motifs. As the (super) cell size increases, the energy is lowered relative to the monomorphous cell, and stabilizes after ~32 formula units (≥160 atoms) are included. Being a result of nonthermal energy minimization in the internal energy without entropy, this correlated set of displacements must represent the intrinsic geometry preferred by the underlying chemical bonding (lone pair bonding), and as such has a different origin than the normal, dynamic thermal disorder modeled by molecular dynamics. Indeed, the polymorphous network, not the monomorphous ansatz, is the kernel structure from which high temperature thermal agitation develops. The emerging physical picture is that the polymorphous network has an average structure with high symmetry, yet the local structural motifs have low symmetries. In this work we find that, compared with monomorphous counterparts, the polymorphous networks have significantly lower predicted total energies, larger band gaps, and ionic dominated dielectric constants, and agree much more closely with the observed pair distribution functions. An analogous polymorphous situation is found in the paraelectric phase of a few cubic oxide perovskites where local polarization takes the role of local displacements in halide perovskites, and in the paramagnetic phases of a few 3 d oxides where the local spin configuration takes that role.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗