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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Anharmonic Origin of the Giant Thermal Expansion of NaBr

ll phonons in a single crystal of NaBr are measured by inelastic neutron scattering at temperatures of 10, 300, and 700 K. Even at 300 K, the phonons, especially the longitudinal-optical phonons, show large shifts in frequencies and show large broadenings in energy owing to anharmonicity. Ab initio computations are first performed with the quasiharmonic approximation (QHA) in which the phonon frequencies depend only on V and on T only insofar as it alters V by thermal expansion. This QHA is an unqualified failure for predicting the temperature dependence of phonon frequencies, even 300 K, and the thermal expansion is in error by a factor of 4. Ab initio computations that include both anharmonicity and quasiharmonicity successfully predict both the temperature dependence of phonons and the large thermal expansion of NaBr. The frequencies of longitudinal-optical phonon modes decrease significantly with temperature owing to the real part of the phonon self-energy from explicit anharmonicity originating from the cubic anharmonicity of nearest-neighbor NaBr bonds. Anharmonicity is not a correction to the QHA predictions of thermal expansion and thermal phonon shifts but dominates the behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1997

The rate of the reaction between CN and C2H2 at interstellar temperatures

The rate coefficient for the important interstellar reaction between CN and C2H2 has been calculated as a function of temperature between 10 and 300 K. The potential surface for this reaction has been determined through ab initio quantum chemical techniques; the potential exhibits no barrier in the entrance channel but does show a small exit channel barrier, which lies below the energy of reactants. Phase-space calculations for the reaction dynamics, which take the exit channel barrier into account, show the same unusual temperature dependence as determined by experiment, in which the rate coefficient at first increases as the temperature is reduced below room temperature and then starts to decrease as the temperature drops below 50-100 K. The agreement between theory and experiment provides strong confirmation that the reaction occurs appreciably at cool interstellar temperatures.

NASA Discipline Exobiology↗
DOE OSTI2026

Quantum scattering of HC 5 N and para -H 2 on a new potential energy surface

In the interstellar medium (ISM), non-local thermodynamic equilibrium situations are common due to low density, and one needs to consider the effect of molecular collisions in order to interpret the observations. Among the species detected in the ISM, cyanopolyynes, with the general molecular formula HC 2n+1 N (n = 1, 2, …), are characterized by large dipole moments and small rotational constants and constitute an indispensable class of candidates for the sensitive tracers of local density and temperature. We present a study of the collisional (de-) excitation of HC 5 N by para -H 2 (p-H 2 ) in its ground rotational state, namely HC 5 N ( j 1 ) + H 2 ( j 2 = 0) → HC 5 N (j$_1^′$) + H2 (j$_2^′$ = 0), where j 1 (or j$_1^′$) and j 2 (or j$_2^′$) denote the initial (or final) rotational quantum numbers of HC 5 N and H 2 , respectively. We performed the quantum scattering calculations at low collision energy using a new four-dimensional ab initio potential energy surface. In the regime where p-H 2 remains in its rotational ground state, converged cross sections did not require including excited rotational states of p-H 2 in the rotational basis. State-to-state cross sections were computed by means of the quantum-mechanical close-coupling (CC) method and the coupled states (CS) approximation, and rate coefficients for the first 61 levels of HC 5 N were computed for the first time up to 20 K with the CC approach and up to 50 K with the CS method. CC and CS results were found to agree well at temperatures up to 20 K. Finally, these data should allow a more accurate derivation of the HC 5 N abundance in molecular clouds.

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

Ionization state and dielectric constant in cold rarefied hydrocarbon plasmas of inertial confinement fusion

A combined approach to study cold rarefied matter is introduced that includes a semianalytical method based on the free-energy minimization and ab initio calculations based on the finite-temperature density-functional theory. The approach is used to calculate the ionization state of hydrocarbon (CH) under the shock-release conditions in inertial confinement fusion. Furthermore, the dielectric constant of CH is calculated using the Kubo-Greenwood formulation and contribution from atomic polarizabilities is found to be as important as the free-electron contribution. Using the ionization state and dielectric constant, the electron density profile in the rarefaction wave of the shock-release plasma is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Metal-insulator transition in V 2 O 3 with intrinsic defects

Vanadium sesquioxide (V 2 O 3 ) is a Mott insulator exhibiting a temperature-dependent metal-insulator transition (MIT) at 165 K accompanied by both a magnetic and structural transition. Although it is expected to be a metal under conventional band theory, electron interactions at low temperature cause it to behave like an insulator, making it difficult to accurately model its electronic properties with standard ab initio methods. As such, accurate theoretical assessments of the MIT with point defects requires special attention to the type of functionals used. In this study, we conclude that the PBE + U functional provides the best compromise between accuracy and efficiency in calculating the properties related to the MIT between low-temperature and high-temperature V 2 O 3 . We use this functional to explore the various influences that intrinsic point defects will have on the MIT in V 2 O 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Discrimination of Hexane Isomers by Temperature Swing Adsorption in a Rigid Aluminum Metal–Organic Framework

The efficient separation of alkane isomers with similar physicochemical properties remains a persistent challenge for the petrochemical industry. Adsorptive separation using metal− organic frameworks (MOFs) offers an energy-efficient alternative to conventional distillation. Herein, we report temperature swing discrimination of hexane isomers with different degrees of branching using MIL-120, a rigid aluminum pyromellitate-based MOF. MIL-120 features uniform one-dimensional channels with an aperture of ∼5.5 Å. At 30 °C, it selectively adsorbs linear and monobranched hexanes while excluding the dibranched isomer. Upon heating to 120 °C, both mono- and dibranched isomers are completely excluded, whereas linear hexane remains strongly adsorbed. Breakthrough experiments validate the temperature swing separation performance. Adsorption heat analysis combined with ab initio calculations provides a quantitative measure of distinct differences in adsorption enthalpies, binding energies, and diffusion barriers responsible for the observed separation efficiency, highlighting the potential of this MOF for efficient separation of alkane isomers via temperature swing adsorption.

Adsorption↗
NASA NTRS1986

Ab initio determination of the proton affinities of small neutral and anionic molecules

The proton affinity of a molecule in the gas phase is a fundamental measure of its basicity and is the factor controlling the course of many ion-molecule reactions. In this article, ab initio molecular orbital theory at the MP4/6-311 ++ G(3df, 3pd) level of theory is demonstrated to predict proton affinities (PA's) for small neutral and anionic bases to within 2 kcal mol-1. Furthermore, the errors are random, indicating that there are likely no systematic errors in either the experimental or theoretical PA's. Also, this level of theory is used to calibrate less sophisticated theoretical models which are suitable for larger molecules; the MP4/6-311 ++ G(2d, 2p) and MP2/6-311 ++ G(d, p) theoretical models should be particularly useful. A procedure for predicting the vibrational frequencies for anion is proposed and applied to CH3-, NH2-, OH-, and CN-.

Non-NASA Center↗
DOE OSTI2020

Modeling Hydrogen Evolution Reaction Kinetics through Explicit Water–Metal Interfaces

Despite the apparent simplicity of the hydrogen evolution reaction (HER) and the decades of research into it, controversy remains in the literature regarding the identity of the active site and the competition between the Heyrovsky and Tafel steps. In this work, we use charge-extrapolated ab initio simulations with explicit water in conjunction with mean-field microkinetic modeling to explore the mechanism for HER on both close-packed (111) and stepped (211) transition metals. First, we show that atop H*, beyond a monolayer of hollow H*, is unlikely to play a role in the HER mechanism, given its very positive adsorption energies. The energetics suggests the Volmer–Heyrovsky mechanism to predominate on fcc transition metals under typical operating conditions. Here, we evaluate our theoretical results vs several experimental observations. We show that the Volmer–Heyrovsky mechanism predicts an activity volcano with its peak at a H* binding ΔG H* ≈ 0 eV, consistent with experiment. In contrast, the Volmer–Tafel volcano shows a broad rate plateau between ΔG H* ≈ 0 eV and ΔG H* ≈ – 0.4 eV. We find our theoretical Tafel slopes to be consistent with experimental ones on a range of transition metals. We show that, in line with experimental observations, the introduction of a CO(g) atmosphere shifts the strong binding metals toward the weak binding leg. Our study suggests that the simple thermodynamic approach to HER activity still holds, even when a detailed kinetic picture is considered.

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

Stabilization of Hydroxide Ions at the Interface of a Hydrophobic Monolayer on Water via Reduced Proton Transfer

We report a joint study using surface-specific sum-frequency vibrational spectroscopy and ab initio molecular dynamics simulations, respectively, on a pristine hydrophobic (sub)monolayer hexane-water interface, namely, the hexane/water interface with varied vapor pressures of hexane and different pHs in water. We show clear evidence that hexane on water revises the interfacial water structure in a way that stabilizes the hypercoordinated solvation structure and slows down the migration of hydroxide ion (OH – ) relative to that in bulk water. Furthermore, this mechanism effectively attracts the OH – to the water-hydrophobic interface with respect to its counterion. The result illustrates the striking difference of proton transfer of hydrated OH – at the interface and in the bulk, which is responsible for the intrinsic charging effect at the hydrophobic interface.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Electrocatalysis in Alkaline Media and Alkaline Membrane-Based Energy Technologies

Hydrogen energy-based electrochemical energy conversion technologies offer the promise of enabling a transition of the global energy landscape from fossil fuels to renewable energy. Here, we present a comprehensive review of the fundamentals of electrocatalysis in alkaline media and applications in alkaline-based energy technologies, particularly alkaline fuel cells and water electrolyzers. Anion exchange (alkaline) membrane fuel cells (AEMFCs) enable the use of nonprecious electrocatalysts for the sluggish oxygen reduction reaction (ORR), relative to proton exchange membrane fuel cells (PEMFCs), which require Pt-based electrocatalysts. However, the hydrogen oxidation reaction (HOR) kinetics is significantly slower in alkaline media than in acidic media. Understanding these phenomena requires applying theoretical and experimental methods to unravel molecular-level thermodynamics and kinetics of hydrogen and oxygen electrocatalysis and, particularly, the proton-coupled electron transfer (PCET) process that takes place in a proton-deficient alkaline media. Extensive electrochemical and spectroscopic studies, on single-crystal Pt and metal oxides, have contributed to the development of activity descriptors, as well as the identification of the nature of active sites, and the rate-determining steps of the HOR and ORR. Among these, the structure and reactivity of interfacial water serve as key potential and pH-dependent kinetic factors that are helping elucidate the origins of the HOR and ORR activity differences in acids and bases. Additionally, deliberately modulating and controlling catalyst-support interactions have provided valuable insights for enhancing catalyst accessibility and durability during operation. The design and synthesis of highly conductive and durable alkaline membranes/ionomers have enabled AEMFCs to reach initial performance metrics equal to or higher than those of PEMFCs. We emphasize the importance of using membrane electrode assemblies (MEAs) to integrate the often separately pursued/optimized electrocatalyst/support and membranes/ionomer components. Operando/in situ methods, at multiscales, and ab initio simulations provide a mechanistic understanding of electron, ion, and mass transport at catalyst/ionomer/membrane interfaces and the necessary guidance to achieve fuel cell operation in air over thousands of hours. We hope that this Review will serve as a roadmap for advancing the scientific understanding of the fundamental factors governing electrochemical energy conversion in alkaline media with the ultimate goal of achieving ultralow Pt or precious-metal-free high-performance and durable alkaline fuel cells and related technologies.

30 DIRECT ENERGY CONVERSION↗
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

Thermal mean-field theories

Several closely related ab initio thermal mean-field theories for fermions, both well-established and new ones, are compared with one another at the formalism level and numerically. The theories considered are Fermi–Dirac theory; thermal Hartree–Fock (HF) theory; two modifications of the thermal single-determinant and the first-order finite-temperature many-body perturbation theory based on a zero-temperature or thermal HF reference. Furthermore, thermal full-configuration-interaction theory is used as the benchmark.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Confinement effects and acid strength in zeolites

Chemical reactivity and sorption in zeolites are coupled to confinement and—to a lesser extent—to the acid strength of Brønsted acid sites (BAS). In presence of water the zeolite Brønsted acid sites eventually convert into hydronium ions. The gradual transition from zeolite Brønsted acid sites to hydronium ions in zeolites of varying pore size is examined by ab initio molecular dynamics combined with enhanced sampling based on Well-Tempered Metadynamics and a recently developed set of collective variables. While at low water content (1–2 water/BAS) the acidic protons prefer to be shared between zeolites and water, higher water contents (n > 2) invariably lead to solvation of the protons within a localized water cluster adjacent to the BAS. At low water loadings the standard free energy of the formed complexes is dominated by enthalpy and is associated with the acid strength of the BAS and the space around the site. Conversely, the entropy increases linearly with the concentration of waters in the pores, favors proton solvation and is independent of the pore size/shape.

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

Tuning aminopolycarboxylate chelators for efficient complexation of trivalent actinides

The complexation of trivalent lanthanides and minor actinides (Am 3+ , Cm 3+ , and Cf 3+ ) by the acyclic aminopolycarboxylate chelators 6,6'-((ethane-1,2-diylbis–((carboxymethyl)azanediyl))bis–(methylene))dipicolinic acid (H 4 octapa) and 6,6'-((((4-(1-(2-(2-(2-hydroxyethoxy)ethoxy)ethyl)-1H-1,2,3-triazol-4-yl)pyridine-2,6-diyl)bis–(methylene))bis–((carboxymethyl)azanediyl))bis–(methylene)) dipicolinic acid (H 4 pypa-peg) were studied using potentiometry, spectroscopy, competitive complexation liquid–liquid extraction, and ab initio molecular dynamics simulations. Two studied reagents are strong multidentate chelators, well-suited for applications seeking radiometal coordination for in-vivo delivery and f-element isolation. The previously reported H 4 octapa forms a compact coordination packet, while H 4 pypa-peg is less sterically constrained due to the presence of central pyridine ring. The solubility of H 4 octapa is limited in a non-complexing high ionic strength perchlorate media. However, the introduction of a polyethylene glycol group in H 4 pypa-peg increased the solubility without influencing its ability to complex the lanthanides and minor actinides in solution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.'s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U parameter for the 5f electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl3) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Melting curves of ice polymorphs in the vicinity of the liquid–liquid critical point

The possible existence of a liquid–liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164–168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid–liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid–liquid critical point lies between the melting curves of ice III and ice V. Herein we explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid–liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.

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

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1996

An ab initio benchmark study of the H + CO --> HCO reaction

The H + CO --> HCO reaction has been characterized with correlation consistent basis sets at five levels of theory in order to benchmark the sensitivities of the barrier height and reaction ergicity to the one-electron and n-electron expansions of the electronic wave function. Single and multireference methods are compared and contrasted. The coupled cluster method RCCSD(T) was found to be in very good agreement with Davidson-corrected internally-contracted multireference configuration interaction (MRCI+Q). Second-order Moller-Plesset perturbation theory (MP2) was also employed. The estimated complete basis set (CBS) limits for the barrier height (in kcal/mol) for the five methods, including harmonic zero-point energy corrections, are MP2, 4.66; RCCSD, 4.78; RCCSD(T), 4.15; MRCI, 5.10; and MRCI+Q, 4.07. Similarly, the estimated CBS limits for the ergicity of the reaction are: MP2, -17.99; RCCSD, -13.34; RCCSD(T), -13.79; MRCI, -11.46; and MRCI+Q, -13.70. Additional basis set explorations for the RCCSD(T) method demonstrate that aug-cc-pVTZ sets, even with some functions removed, are sufficient to reproduce the CBS limits to within 0.1-0.3 kcal/mol.

NASA Discipline Exobiology↗