skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Ab initio simulations of liquid carbon monoxide at high pressure

Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
SC0010517; SC0016248
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
High Energy Density Physics
Additional Journal Information:
Journal Volume: 22; Journal Issue: C; Related Information: CHORUS Timestamp: 2017-10-04 21:05:13; Journal ID: ISSN 1574-1818
Country of Publication:

Citation Formats

Leonhardi, Tanis C., and Militzer, Burkhard. Ab initio simulations of liquid carbon monoxide at high pressure. Netherlands: N. p., 2017. Web. doi:10.1016/j.hedp.2017.02.005.
Leonhardi, Tanis C., & Militzer, Burkhard. Ab initio simulations of liquid carbon monoxide at high pressure. Netherlands. doi:10.1016/j.hedp.2017.02.005.
Leonhardi, Tanis C., and Militzer, Burkhard. Wed . "Ab initio simulations of liquid carbon monoxide at high pressure". Netherlands. doi:10.1016/j.hedp.2017.02.005.
title = {Ab initio simulations of liquid carbon monoxide at high pressure},
author = {Leonhardi, Tanis C. and Militzer, Burkhard},
abstractNote = {},
doi = {10.1016/j.hedp.2017.02.005},
journal = {High Energy Density Physics},
number = C,
volume = 22,
place = {Netherlands},
year = {Wed Mar 01 00:00:00 EST 2017},
month = {Wed Mar 01 00:00:00 EST 2017}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.hedp.2017.02.005

Save / Share:
  • We have performed a detailed study of molecular dissociation in liquid hydrogen using both Born-Oppenheimer molecular dynamics with Density Functional Theory and Coupled Electron-Ion Monte Carlo simulations. We observe a range of densities where (dP/d{rho}){sub T} = 0 that coincides with sharp discontinuities in the electronic conductivity, which is clear evidence of the plasma phase transition for temperatures 600K {le} T {le} 1500K. Both levels of theory exhibit the transition, although Quantum Monte Carlo predicts higher transition pressures. Based on the temperature dependence of the discontinuity in the electronic conductivity, we estimate the critical point of the transition at temperaturesmore » slightly below 2000 K. We examine the influence of proton zero point motion by using Path Integral Molecular Dynamics with Density Functional Theory; the main effect is to shift the transition to lower pressures. Furthermore, we calculate the melting curve of molecular hydrogen up to pressures of 200 GPa, finding a reentrant melting line in good agreement with previous calculations. The melting line crosses the metalization line at 700 K and 220 GPa using density functional energetics and at 550 K and 290 GPa using Quantum Monte Carlo energetics.« less
  • Crystalline polyethylene was investigated under pressure between 0 and 40 GPa, up to 280 deg. C, by means of synchrotron x-ray powder diffraction and ab initio calculations. A rich polymorphism was unveiled, consisting of two new high-pressure monoclinic phases, in addition to the well-known orthorhombic one, which appear reversibly, although with strong hysteresis, upon increasing pressure above 6 GPa (P2{sub 1}/m, Z{sub chain}=1) and 14-16 GPa (A2/m, Z{sub chain}=2), respectively. The equation of state was determined for the three solid phases. We find that polyethylene is characterized by a sharp separation between strong covalent intrachain and weaker interchain interactions upmore » to the maximum investigated pressure, which, in turn, places the ultimate chemical stability limit of polyethylene far beyond these thermodynamic conditions.« less
  • A computational strategy is devised for the accurate ab initio simulation of elastic properties of crystalline materials under pressure. The proposed scheme, based on the evaluation of the analytical stress tensor and on the automated computation of pressure-dependent elastic stiffness constants, is implemented in the CRYSTAL solid state quantum-chemical program. Elastic constants and related properties (bulk, shear and Young moduli, directional seismic wave velocities, elastic anisotropy index, Poisson's ratio, etc.) can be computed for crystals of any space group of symmetry. We apply such a technique to the study of high-pressure elastic properties of three silicate garnet end-members (namely, pyrope,more » grossular, and andradite) which are of great geophysical interest, being among the most important rock-forming minerals. The reliability of this theoretical approach is proved by comparing with available experimental measurements. The description of high-pressure properties provided by several equations of state is also critically discussed.« less
  • We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm –3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreementmore » with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.« less
  • Production of C{sub 1} to C{sub 6} alcohols from syngas is one of the most promising developments in C{sub 1} chemistry. The reaction mechanism, which obeys a Schultz--Flory polymerization law, leads to simultaneous production of alcohols and hydrocarbons. The rate of heavier alcohols and hydrocarbons decreases with their increasing carbon number. A water-gas shift reaction occurs as a side reaction. Therefore, it is necessary to know the vapor-liquid equilibria of the H{sub 2}, CO, CO{sub 2}, CH{sub 4}, C{sub 2}H{sub 6}, C{sub 3}H{sub 8}, H{sub 2}O, CH{sub 3}OH, C{sub 2}H{sub 5}OH, and {ital n}-C{sub 3}H{sub 7}OH system for separation processmore » design. In a previous publication, the authors reported vapor-liquid equilibrium data on CO{sub 2}-alcohols, CH{sub 4}-alcohols, and C{sub 2}H{sub 6}-alcohols. In this work the authors report isothermal high-pressure vapor-liquid equilibrium data for the H{sub 2}, CO, CO{sub 2}, CH{sub 4}, C{sub 2}H{sub 6}, C{sub 3}H{sub 8}, H{sub 2}O, CH{sub 3}OH, C{sub 2}H{sub 5}OH, and {ital n}-C{sub 3}H{sub 7}OH system in the 313--333 K temperature range.« less