Search NASA⌕ Search

SEARCH · Search NASA

Results for “pressure dependence”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. Here, the parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

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

Formation of Phenanthrene via Recombination of Indenyl and Cyclopentadienyl Radicals: A Theoretical Study

This work presents quantum chemical G3(MP2,CC)//B2PLYPD3/6-311G(d,p) calculations of the potential energy surface for the indenyl (C 9 H 7 ) + cyclopentadienyl (C 5 H 5 ) reaction followed by unimolecular decomposition of the C 14 H 11 radicals formed as the primary products, as well as the Rice–Ramsperger–Kassel–Marcus master equation (RRKM-ME) calculations to predict temperature- and pressure-dependent reaction rate constants and product branching ratios. The reaction begins with the barrierless recombination of indenyl and cyclopentadienyl forming a C 14 H 12 molecule with a new C–C bond connecting two five-membered rings, which subsequently dissociates to C 14 H 11 radicals by H losses. The primary products of the C 9 H 7 + C 5 H 5 → C 14 H 11 + H reaction can directly decompose by another H loss to benzofulvalene, and this pathway is most favorable in terms of the entropy factor and hence is preferable at higher temperatures. Otherwise, the initial C 14 H 11 isomers can undergo significant structural rearrangements before eliminating an H atom and producing phenanthrene, anthracene, or benzoazulenes, among which the formation of phenanthrene via the “spiran” pathway is clearly preferred. The calculated barriers along the computed favorable dissociation pathways are relatively low, in the ~30–40 kcal/mol range, making the C 14 H 11 radicals unstable at temperatures above 1000–1250 K at 1 atm. The results of RRKM-ME calculations show that, under typical combustion conditions, the decomposition of the C 14 H 11 radicals predominantly leads to benzofulvalene. However, the latter can be rapidly converted to phenanthrene via H-assisted isomerization with the rate constant for the benzofulvalene + H → phenanthrene + H reaction being close to 10 –11 cm 3 molecule –1 s –1 at 1000–1500 K and 1 atm. Here, the results provide further support for the hypothesis that recombination of two π radicals containing five-membered rings can lead to a growth of PAH with the formation of two fused six-membered rings, but the reaction mechanism may not be direct and is likely to involve two consecutive H atom losses leading to a fulvalene-like product, with subsequent H-assisted isomerization of the latter to a benzenoid PAH.

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

Theoretical Kinetics Predictions for Reactions on the NH 2 O Potential Energy Surface

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH 2 O potential energy surface, including the NH 2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH 3 oxidation with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations, including large basis sets and high-level corrections to yield high-accuracy (~0.2 kcal/mol 2σ uncertainty) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature- and pressure-dependent kinetics. For radical–radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constants as well as the branching ratios. We find that the NO + H 2 channel contributes 10% of the total NH 2 + O flux at combustion temperatures, although this channel is not included in modern NH 3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH 3 oxidation and NO x formation. Finally, the present results for NH 2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH 3 and thermal DeNO x .

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

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

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

Realization of Organocerium-Based Fullerene Molecular Materials Showing Mott Insulator-Type Behavior

Abstract Electron-rich organocerium complexes (C 5 Me 4 H) 3 Ce and [(C 5 Me 5 ) 2 Ce(ortho-oxa)], with redox potentials E 1/2 = –0.82 V and E 1/2 = –0.86 V versus Fc/Fc + respectively, were reacted with fullerene (C 60 ) in different stoichiometries to obtain molecular materials. Structurally characterized co-crystals: [(C 5 Me 4 H) 3 Ce] 2 •C 60 (1) and [(C 5 Me 5 ) 2 Ce(ortho-oxa)] 3 •C 60 (2) of C 60 with cerium-based rare earth molecular precursors are reported for the first time. The extent of charge transfer in 1 and 2 was evaluated using a series of physical measurements: FT-IR, Raman, solidstate UV-vis-NIR spectroscopy, X-ray absorption near edge structure (XANES) spectroscopy, and magnetic susceptibility measurements. The physical measurements indicate that 1 and 2 comprise the cerium(III) oxidation state with formally neutral C 60 as a co-crystal in both cases. Pressure-dependent periodic density functional theory calculations were performed to study the electronic structure of 1. Inclusion of a Hubbard-U parameter removes Ce f states from the Fermi level, opens up a band gap, and stabilizes FM/AFM magnetic solutions that are isoenergetic because of the large distances between the Ce(III) cations. Furthermore, the electronic structure of this strongly correlated Mott insulator-type system is reminiscent of the well-studied Ce 2 O 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dehydrogenation vs Apparent Hydrogenation: Unraveling the Mechanisms of He and O 2 Plasma Etching on Colloidal Nanocrystal Films

Removing organic ligands from colloidal nanoparticles is critical for fabricating solid-state devices, yet accurately quantifying this removal remains a significant analytical challenge. Here, we establish a robust and accessible method for this quantification by calibrating Raman spectroscopy against precise ion beam analysis (IBA) for nanoparticle assemblies (CNAs) processed by helium (He) and oxygen (O 2 ) plasmas. We demonstrate that the calibration curves are remarkably independent of plasma power and pressure, depending critically only on the choice of feed gas. He plasma induces rapid dehydrogenation and cross-linking, evidenced by a much faster decrease in the C–H Raman signal relative to the actual carbon loss. Conversely, O 2 plasma leads to a surprising “apparent hydrogenation”, where the carbon backbone is removed significantly faster than the C–H signal diminishes. This counterintuitive effect is explained by a serial mechanism of oxidative fragmentation; β-scission cleaves the alkyl chains, and subsequent stabilization steps enrich the remaining film with hydrogen-rich methyl-terminated fragments, while carbon is efficiently removed as volatile CO. This work provides calibrated functions that enable the rapid determination of absolute carbon content in processed CNAs using simple Raman spectroscopy with uncertainties of ∼8% for O 2 and ∼12% for He plasma, offering a vital tool for both process diagnostics and fundamental studies of plasma–matter interactions in colloidal nanocrystal films.

Animal feed↗
DOE OSTI2025

Vibrational Dynamics and Phase Transitions of Hydrazine to 50 GPa

The high pressure behavior of hydrazine, N 2 H 4 , has been investigated to 50 GPa at room temperature using infrared and Raman spectroscopy to explore pressure induced phase transitions and changes in hydrogen bonding. Three solid–solid phase transitions were detected at 11, 21, and 32 GPa on room temperature compression through dramatic changes in the lattice vibration and N–H stretching regions with increasing pressure in both measurement techniques. The transition to phase IV, which appears at 32 GPa, exhibits increased hydrogen bonding with significant hysteresis, persisting to 9 GPa on decompression. This work presents a detailed analysis of the pressure dependence of mode shifts and calculations of mode Grüneisen parameters as well as a determination of an approximate thermodynamic Grüneisen parameter. We compare these results to the behavior of other small molecular materials such as ammonia and water and explore the evolution of hydrogen bonding in hydrazine toward the symmetrically hydrogen bonded state, which has previously been suggested by theoretical computations.

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

A Quantitative Examination of Catalyst Component Impacts on Hydrogen Adsorption and Spillover

Hydrogen spillover is a widely recognized but poorly understood surface phenomenon. We coupled volumetric chemisorption with Fourier transform infrared spectroscopy, thermogravimetric analysis, and acid–base titrations to quantify the roles of metal (Pt, Au) and support chemistry (anatase, rutile, P25 titania) on spillover. We show that metal identity has little to no impact on the amount of spillover hydrogen. Spillover pressure dependence is essentially the same on Au/TiO 2 and Pt/TiO 2 , indicating H 2 adsorption is equilibrated across both the metal and support. Conversely, spillover is highly sensitive to support surface chemistry, as rutile TiO 2 stabilizes ∼5× more surface hydrogen than anatase TiO 2 . This change is due to a combination of factors, including the rutile surface’s higher proton affinity, stabilized surface electronic states, and larger surface entropy. This work highlights the utility of Au/MO x catalysts as control materials for spillover, as they enable quantitative evaluation of H 2 adsorption on active metals and spillover onto the support.

Adsorption↗
DOE OSTI2021

Effects of Hydrogen on the Phase Relations in Fe-FeS at Pressures of Mars-Sized Bodies

The large radius, and therefore low density, of the Martian core found in the InSight mission data analysis highlights the importance of considering other light elements besides sulfur (S), which has been considered as the main light element for Mars for decades. Hydrogen (H) is abundant in the solar system and becomes siderophile at high pressures. Although Fe-S and Fe-H systems have been studied individually, the Fe-S-H ternary system has only been investigated up to 16 GPa and 1723 K. We have investigated the Fe-S-H system at pressures and temperatures (P-T) relevant to the cores of Mars-sized planets (up to 45 GPa and well above the melting temperature of FeS) in the laser-heated diamond anvil cell combined with in situ synchrotron X-ray diffraction. Here we found that sufficient hydrogen leads to the disappearance of Fe 3 S at high P-T. Instead, separate Fe-H and Fe-S phases appear at 23–35 GPa. At pressures above 35 GPa, we found a new phase appearing while Fe-S phases disappear and Fe-H phases remain. Our analysis indicates that the new phase likely contains both S and H in the crystal structure (tentatively FeS x H y where x ≈ 1 and y ≈ 1). The observed pressure-dependent changes in the phase relation may be important for understanding the structure and dynamics of the Martian core and the cores of Mars-sized exoplanets.

58 GEOSCIENCES↗
DOE OSTI2021

Manifestations of metastable criticality in the long-range structure of model water glasses

Abstract Much attention has been devoted to water’s metastable phase behavior, including polyamorphism (multiple amorphous solid phases), and the hypothesized liquid-liquid transition and associated critical point. However, the possible relationship between these phenomena remains incompletely understood. Using molecular dynamics simulations of the realistic TIP4P/2005 model, we found a striking signature of the liquid-liquid critical point in the structure of water glasses, manifested as a pronounced increase in long-range density fluctuations at pressures proximate to the critical pressure. By contrast, these signatures were absent in glasses of two model systems that lack a critical point. We also characterized the departure from equilibrium upon vitrification via the non-equilibrium index; water-like systems exhibited a strong pressure dependence in this metric, whereas simple liquids did not. These results reflect a surprising relationship between the metastable equilibrium phenomenon of liquid-liquid criticality and the non-equilibrium structure of glassy water, with implications for our understanding of water phase behavior and glass physics. Our calculations suggest a possible experimental route to probing the existence of the liquid-liquid transition in water and other fluids.

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

Effects of Oxygen Modification on the Structural and Magnetic Properties of Highly Epitaxial La0.7Sr0.3MnO3 (LSMO) thin films

Abstract La 0.7 Sr 0.3 MnO 3 , a strong semi-metallic ferromagnet having robust spin polarization and magnetic transition temperature ( T C ) well above 300 K, has attracted significant attention as a possible candidate for a wide range of memory, spintronic, and multifunctional devices. Since varying the oxygen partial pressure during growth is likely to change the structural and other physical functionalities of La 0.7 Sr 0.3 MnO 3 (LSMO) films, here we report detailed investigations on structure, along with magnetic behavior of LSMO films with same thickness (~30 nm) but synthesized at various oxygen partial pressures: 10, 30, 50, 100, 150, 200 and 250 mTorr. The observation of only ( 00 l ) reflections without any secondary peaks in the XRD patterns confirms the high-quality synthesis of the above-mentioned films. Surface morphology of the films reveals that these films are very smooth with low roughness, the thin films synthesized at 150 mTorr having the lowest average roughness. The increasing of magnetic T C and sharpness of the magnetic phase transitions with increasing oxygen growth pressure suggests that by decreasing the oxygen growth pressure leads to oxygen deficiencies in grown films which induce oxygen inhomogeneity. Thin films grown at 150 mTorr exhibits the highest magnetization with T C = 340 K as these thin films possess the lowest roughness and might exhibit lowest oxygen vacancies and defects. Interpretation and significance of these results in the 30 nm LSMO thin films prepared at different oxygen growth pressures are also presented, along with the existence and growth pressure dependence of negative remanent magnetization (NRM) of the above-mentioned thin films.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

An Experimental and Theoretical Study of the High Temperature Reactions of the Four Butyl Radical Isomers

The high temperature gas phase chemistry of the four butyl radical isomers (n-butyl,sec-butyl, iso-butyl, andtert-butyl) was investigated in a combined experimental and theoretical study. Organic nitrites were used as convenient and clean sources of each of the butyl radical isomers. Rate coefficients for dissociation of each nitrite were obtained experimentally and are at, or close to, the high pressure limit. Low pressure experiments were performed in a diaphragmless shock tube with laser schlieren densitometry at post-shock pressures of 65, 130, and 260 Torr and post-shock temperatures of 700-1000 K. Additional experiments were conducted with iso-butyl radicals at 805 K and 8.7 bar to elucidate changes in mechanism at higher pressures. These experiments were performed in a miniature shock tube with synchrotron-based photoionization mass spectrometry. The mass spectra confirmed that scission of the O-NO bond is the primary channel by which the precursors dissociate, but they also provided evidence of a minor channel (<7.7%) through HNO loss and formation of an aldehyde. These high pressure experiments were also used to determine the disproportionation/recombination ratio for iso-butyl radicals as 0.3. Reanalysis of the lower-temperature literature and the present data yielded rate constants for the disproportionation reaction, iso-butyl + iso-butyl = iso-butene + iso-butane. A chemical kinetics model was developed for the reactions of the butyl isomers that included new paths for highly energized adducts. These adducts are formed by the addition of H, CH 3 or C 2 H 5 to the butyl radicals. Accompanying theoretical investigations show that chemically activated pathways are competitive with stabilization of the adduct by collision under the conditions of the laser schlieren experiments. These calculations also show that at 10 bar and T< 1000 K stabilization is the only important reaction, but at higher temperatures, even at 10 bar, chemically activated product channels should also be considered. Branching fractions and rate coefficients are presented for these reactions. Furthermore, this study also highlights the importance of the radical structure for determining branching ratios for disproportionation and recombination of alkyl radicals, and these were facilitated by theoretical calculations of recombination rate coefficients for the four butyl radical isomers. Overall, the results reveal previously unknown features of butyl radical chemistry under conditions that are relevant to a wide range of applications and reaction mechanisms are presented that incorporate pressure dependent rate coefficients for the key steps.

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

Accurate prediction of carbon dioxide capture by deep eutectic solvents using quantum chemistry and a neural network

Carbon dioxide (CO 2 ) emissions from fossil fuel combustion are a significant source of greenhouse gas, contributing in a major way to global warming and climate change. Carbon dioxide capture and sequestration is gaining much attention as a potential method for controlling these greenhouse gas emissions. Among the environmentally friendly solvents, deep eutectic solvents (DESs) have demonstrated the potential capability for carbon capture. To establish a theoretical framework for DES activity, thermodynamics modeling and solubility predictions are significant factors to anticipate and understand the system behavior. Here, in this study, we combine the COSMO-RS model with machine learning techniques to predict the solubility of CO 2 in various deep eutectic solvents. A comprehensive data set was established comprising 1973 CO 2 solubility data points in 132 different DESs at a variety of temperatures, pressures, and DES molar ratios. This data set was then utilized for the further verification and development of the COSMO-RS model. The CO 2 solubility (ln(x CO 2 )) in DESs calculated with the COSMO-RS model differs significantly from the experiment with an average absolute relative deviation (AARD) of 23.4%. A multilinear regression model was developed using the COSMO-RS predicted solubility and a temperature-pressure dependent parameter, which improved the AARD to 12%. Finally, a machine learning model using COSMO-RS-derived features was developed based on an artificial neural network algorithm. The results are in excellent agreement with the experimental CO 2 solubilities, with an AARD of only 2.72%. The ML model will be a potentially useful tool for the design and selection of DESs for CO 2 capture and utilization.

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

Permanently porous oligoMOFs and polyMOFs of MIL-53

Enhancing the porosity of metal–organic frameworks (MOFs) typically relies on substituting their organic ligand components. Furthermore, we show that cross-linked terephthalate ligands yield oligoMOF and polyMOF versions of MIL-53(In) in which pressure-dependent phase transformations no longer occur. These materials show enhanced uptakes and CO 2 /N 2 selectivity, while retaining the overall MIL-53 framework.

MOFs↗
DOE OSTI2026

Conformal hetero-electrolyte interface between soft oxyhalides and garnet enables low-pressure lithium-reservoir-free solid-state batteries

The operation of solid-state batteries with a lithium metal anode and a high voltage cathode requires solid electrolytes (SEs) that are chemically stable with lithium, have a wide electrochemical window, and accommodate volume changes in the electrodes. Unfortunately, no SE has exhibited satisfactory mechanical and electrochemical properties that fit these requirements to date. Dual solid electrolyte systems that use a different SE for the anolyte and catholyte present a viable solution. Here, we focus on oxyhalide SEs that demonstrate superior ionic conductivity and cathode compatibility, where their lithium metal reactivity and poor reduction stability can be resolved using a lithium garnet (Li 6.5 La 3 Zr 1.5 Ta 0.5 O 12 , LLZTO) separator. Nonetheless, this imposes a new hetero-electrolyte (H-E) interface at the anolyte|catholyte contact that defines the ion transport across the boundary. We report its promising properties, which are deconvoluted from the electrical measurements of bilayer symmetric cells, for three representative oxyhalide catholytes, LiNbOCl 4 , LiTaOCl 4 , and Li 3 Al 3 O 2 Cl 8 . Pressure-dependent measurements reveal that the relative softness of the oxyhalides (hardness ≤0.4 GPa) enables H-E resistances lower than 150 Ω cm 2 at 2–3 MPa. Mesoscale modelling reveals that the transfer-active contact area of oxyhalides with LLZTO is about 2–3-fold higher than that of argyrodite, Li 6 PS 5 Cl. The low H-E resistance of the LiNbOCl 4 |Li 6.5 La 3 Zr 1.5 Ta 0.5 O 12 dual electrolyte enables the cycling of a Li|LiNi 0.82 Mn 0.07 Co 0.11 O 2 full cell with a high discharge capacity (200 mA h g −1 ) at 60 °C and ∼7 MPa. Importantly, we demonstrate a Li-reservoir-free full cell with high Coulombic efficiency (>99.5%) and capacity at 1 MPa using this approach coupled with a garnet-silver interlayer.

Palmer, Max [University of California, Santa Barba↗
DOE OSTI2020

Heterostructure design to achieve high quality, high density GaAs 2D electron system with g -factor tending to zero

Hydrostatic pressure is a useful tool that can tune several key parameters in solid state materials. For example, the Landé g-factor in GaAs two-dimensional electron systems (2DESs) is expected to change from its bulk value g ≃ -0.44 to zero and even to positive values under sufficiently large hydrostatic pressure. Although this presents an intriguing platform to investigate electron-electron interaction in a system with g = 0, studies are quite limited because the GaAs 2DES density decreases significantly with increasing hydrostatic pressure. Here, we show that a simple model, based on pressure-dependent changes in the conduction band alignment, quantitatively explains this commonly observed trend. Furthermore, we demonstrate that the decrease in the 2DES density can be suppressed by more than a factor of 3 through an innovative heterostructure design.

Physics↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗