Engineering Papers⌕ Search

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.

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

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ab initio and force field molecular dynamics study of bulk organophosphorus and organochlorine liquid structures

Here, we performed ab initio molecular dynamics (AIMD) simulations to benchmark bulk liquid structures and to evaluate results from all-atom force field molecular dynamics (FFMD) simulations with the generalized Amber force field (GAFF) for organophosphorus (OP) and organochlorine (OC) compounds. Our work also addresses the current and important topic of force field validation, applied here to a set of nonaqueous organic liquids. Our approach differs from standard treatments, which validate force fields based on thermodynamic data. Utilizing radial distribution functions (RDFs), our results show that GAFF reproduces the AIMD-predicted asymmetric liquid structures moderately well for OP compounds that contain bulky alkyl groups. Among the OCs, RDFs obtained from FFMD overlap well with AIMD results, with some offsets in position and peak structuring. However, re-parameterization of GAFF for some OCs is needed to reproduce fully the liquid structures predicted by AIMD. The offsets between AIMD and FFMD peak positions suggest inconsistencies in the developed force fields, but, in general, GAFF is able to capture short-ranged and long-ranged interactions of OPs and OCs observed in AIMD. Along with the local coordination structure, we also compared enthalpies of vaporization. Overall, calculated bulk properties from FFMD compared reasonably well with experimental values, suggesting that small improvements within the FF should focus on parameters that adjust the bulk liquid structures of these compounds.

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

The Li–F–H ternary system at high pressures

Evolutionary crystal structure prediction searches have been employed to explore the ternary Li-F-H system at 300 GPa. Metastable phases were uncovered within the static lattice approximation, with LiF 3 H 2 , LiF 2 H, Li 3 F 4 H, LiF 4 H 4 , Li 2 F 3 H and LiF 3 H lying within 50 meV/atom of the 0 K convex hull. All of these phases contain H n F¯ n+1 (n = 1; 2) anions, and Li + cations. Other structural motifs such as LiF slabs, $H$ $^{+}_{3}$ molecules and F δ- ions are present in some of the low enthalpy Li-F-H structures. The bonding within the H n F¯ n+1 molecules, which may be bent or linear, symmetric or asymmetric, is analyzed. The five phases closest to the hull are insulators, while LiF 3 H is metallic and predicted to have a vanishingly small superconducting critical temperature. Li 3 F 4 H is predicted to be stable at zero pressure. Furthermore, this study lays the foundation for future investigations of the role of temperature and anharmonicity on the stability and properties of compounds and alloys in the Li-F-H ternary system.

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

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

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

A chemical kinetics perspective on thermoelectric transport

There is a growing community of chemists and materials scientists exploring thermoelectric properties of organic, polymeric, hybrid, and composite materials. Many of these materials are nondegenerate semiconductors, meaning that the Fermi and transport charge transport energy levels are significantly offset. Others are ionic conductors. While the meaning of the essential parameters electrical conductivity (σ) and Seebeck coefficient (S, the voltage difference per degree of temperature difference) are accessible, the origins of S are not readily apparent to one trained mainly in chemistry. The purpose of this manuscript is to illustrate S as being the result of a steady state difference between concentrations of two designations of charge carriers, those on hotter and those on colder sides of a material sample, analogous to the equilibrium among interconverting chemical species that differ in enthalpy and entropy. The ion-based Seebeck coefficient, also known as the Soret effect, can be explained using principles similar to those applicable to electrons and holes. We hope that this analysis leads to wider understanding of the origins of S through an explanation using the language of chemical kinetics and thermodynamics, and appreciation of ways that thermoelectric efficiency can be enhanced in emerging materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronic states and transitions of PrO and PrO + probed by threshold ionization spectroscopy and spin–orbit multiconfiguration perturbation theory

The precise ionization energy of praseodymium oxide (PrO) seeded in supersonic molecular beams is measured with mass-analyzed threshold ionization (MATI) spectroscopy. Here, a total of 33 spin–orbit (SO) states of PrO and 23 SO states of PrO + are predicted by second-order multiconfigurational quasi-degenerate perturbation (MCQDPT2) theory. Electronic transitions from four low-energy SO levels of the neutral molecule to the ground state of the singly charged cation are identified by combining the MATI spectroscopic measurements with the MCQDPT2 calculations. The precise ionization energy is used to reassess the ionization energies and the reaction enthalpies of the Pr + O → PrO + + e – chemi-ionization reaction reported in the literature. An empirical formula that uses atomic electronic parameters is proposed to predict the ionization energies of lanthanide monoxides, and the empirical calculations match well with available precise experimental measurements.

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

Influence of mass ablation on ignition and burn propagation in layered fusion capsules

After decades of research, recent laser-driven inertial fusion experiments have demonstrated rapid progress toward achieving thermonuclear ignition using capsule designs with cryogenic fuel layers. The ignition physics for these layered capsules involves a complex interplay between the dynamically forming hot spot and the dense surrounding fuel. Using analytic theory and numerical simulations, we demonstrate that the mass ablation rate into the hot spot depends sensitively upon the temperature of the dense fuel, resulting in ablative inflows up to [Formula: see text] faster than previous estimates. This produces an enthalpy flux into the hot spot that plays a critical role in controlling the hot spot temperature, the ignition threshold, and the subsequent burn propagation. The net influence of mass ablation on the ignition threshold is regulated by a dimensionless parameter that depends upon the temperature of the dense fuel. As a consequence, the ignition threshold is sensitive to any mechanism that heats the dense fuel, such as neutrons or radiation emitted from the hot spot. These predictions are confirmed using radiation-hydrodynamic simulations for a series of capsules near ignition conditions. This analysis may have relevance for understanding the variable performance of recent experiments and for guiding new capsule designs toward higher fusion yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

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

Boron coordination change in barium borate melts and glasses and its contribution to configurational heat capacity, entropy, and fragility

High-energy x-ray diffraction from molten and glassy BaB 2 O 4 and BaB 4 O 7 has been performed using aerodynamic levitation and laser heating over a wide range of temperatures. Remarkably, even in the presence of a heavy metal modifier dominating x-ray scattering, it was possible to extract accurate values for the tetrahedral, sp 3 , boron fraction, N 4 , which declines with increasing temperature, using bond valence-based mapping from the measured mean B–O bond lengths while accounting for vibrational thermal expansion. These are used within a boron-coordination-change model to extract enthalpies, ΔH, and entropies, ΔS, of isomerization between sp 2 and sp 3 boron. The results for BaB 4 O 7 , ΔH = 22(3) kJ mol -1 boron, ΔS = 19(2) J mol -1 boron K -1 , agree quantitatively with those found previously for Na 2 B 4 O 7 . Analytical expressions for N 4 (J, T) and associated configurational heat capacity, $C_{P}^{conf}$(J, T), and entropy, S conf (J, T), contributions are extended to cover a wide composition range 0 ≤ J = BaO/B 2 O 3 ≤ 3 using a model for ΔH(J) and ΔS(J) derived empirically for lithium borates. Maxima in the $C_{P}^{conf}$(J, T g ) and fragility index contributions are thereby predicted for J ≲ 1, higher than the maximum observed and predicted in N 4 (J, T g ) at J ≃ 0.6. We discuss the utility of the boron-coordination-change isomerization model in the context of borate liquids containing other modifiers and the prospect of neutron diffraction to aid in empirical determinations of modifier-dependent effects, illustrated by new neutron diffraction data on Ba 11 B 4 O 7 glass, its well-known α-polymorph, and lesser-known δ-phase.

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

Predicted thermophysical properties of UN, PuN, and (U,Pu)N

Molecular dynamics and density functional theory simulations are used to predict the lattice and electronic contributions of thermophysical properties for UN, PuN, and mixed (U,Pu)N systems. The properties predicted include the lattice parameter, linear thermal expansion, enthalpy, and specific heat capacity, as a function of temperature. The simulation predictions for high temperature specific heat capacity are compared against experimental measurements to understand the behavior, and why differences in the experimental measurements are observed. The influence of adding U vacancies, N interstitials, and Pu to UN is also examined. For this, a new PuN potential parameter set is developed and used with the Kocevski UN potential, enabling the dynamics of mixed (U,Pu)N systems to be studied. How defects impact the thermophysical properties is important for understanding fuel behavior under different reactor conditions, and these mechanistic predictions can be used to support fuel performance codes where data is scarce.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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↗