High-pressure polymeric phases of carbon dioxide

Phase Transition in Silicon from Machine Learning Informed Metadynamics

Investigating reconstructive phase transitions in large-sized systems requires a highly efficient computational framework with computational cost proportional to the system size. Traditionally, widely used frameworks such as density functional theory (DFT) have been prohibitively expensive for extensive simulations on large systems that require long-time scales. To address this challenge, this study employed well-trained machine learning potential to simulate phase transitions in a large-size system. This work integrates the metadynamics simulation approach with machine learning potential, specifically deep potential, to enhance computational efficiency and accelerate the study of phase transition and consequent development of grains and dislocation defects in a system. The new method is demonstrated using the phase transitions of bulk silicon under high pressure. This approach has revealed the transition path and formation of polycrystalline silicon systems under specific stress conditions, demonstrating the effectiveness of deep potential-driven metadynamics simulations in gaining insights into complex material behaviors in large-sized systems.

Pressure-Driven Polymorphic Transition, Emergent Insulator-To-Metal Transition, and Photoconductivity Switching in Violet Phosphorus

Polymorphic phase transition is an essential phenomenon in condensed matter that the physical properties of materials may undergo significant changes due to the structural transformation. Phase transition has thus become an important means and dimension for regulating material properties. Herein, this study demonstrates the pressure-induced multi-transition of both structure and physical properties in violet phosphorus, a novel phosphorus allotrope. Under compression, violet phosphorus undergoes sequential polymorphic phase transitions. Concomitant with the first phase transition, violet phosphorus exhibits emergent insulator-metal transition, superconductivity, and dramatic switching from positive to negative photoconductivity. Remarkably, the resistance of violet phosphorus shows a sudden drop of around 107 along with the phase transition. In addition, piezochromism from translucent red to opaque black and suppression of photoluminescence are observed upon compression. Of particular interest is that the sample irreversibly transforms into black phosphorus with a pronounced discrepancy in physical properties from the pristine violet phosphorus after decompression. The abundant polymorphic transitions and property changes in violet phosphorus have significant implications for designing novel pressure-responsive electronic/optoelectronic devices and exploring concealed polymorphic transition materials.

Size-Dependent Nucleation in Crystal Phase Transition from Machine Learning Metadynamics

In this Letter, we present a framework that combines machine learning potential (MLP) and metadynamics to investigate solid-solid phase transition. Based on the spectral descriptors and neural networks regression, we develop a scalable MLP model to warrant an accurate interpolation of the energy surface where two phases coexist. Applying it to the simulation of B4-B1 phase transition of GaN under 50 GPa with different model sizes, we observe sequential change of the phase transition mechanism from collective modes to nucleation and growths. When the size is at or below 128 000 atoms, the nucleation and growth appear to follow a preferred direction. At larger sizes, the nuclei occur at multiple sites simultaneously and grow to microstructures by passing the critical size. The observed change of the atomistic mechanism manifests the importance of statistical sampling with large system size in phase transition modeling.

Theoretical methods for structural phase transitions in elemental solids at extreme conditions: statics and dynamics

In recent years, theoretical studies have moved from a traditionally supporting role to a more proactive role in the research of phase transitions at high pressures. In many cases, theoretical prediction leads the experimental exploration. This is largely owing to the rapid progress of computer power and theoretical methods, particularly the structure prediction methods tailored for high-pressure applications. This review introduces commonly used structure searching techniques based on static and dynamic approaches, their applicability in studying phase transitions at high pressure, and new developments made toward predicting complex crystalline phases. Successful landmark studies for each method are discussed, with an emphasis on elemental solids and their behaviors under high pressure. The review concludes with a perspective on outstanding challenges and opportunities in the field.

Nucleating a Different Coordination in a Crystal under Pressure: A Study of the B1-B2 Transition in NaCl by Metadynamics

Here we propose an NPT metadynamics simulation scheme for pressure-induced structural phase transitions, using coordination number and volume as collective variables, and apply it to the reconstructive structural transformation B1-B2 in NaCl. By studying systems with size up to 64 000 atoms we reach a regime beyond collective mechanism and observe transformations proceeding via nucleation and growth. We also reveal the crossover of the transition mechanism from Buerger-like for smaller systems to Watanabe-Tolédano for larger ones. The scheme is likely to be applicable to a broader class of pressure-induced structural transitions, allowing study of complex nucleation effects and bringing simulations closer to realistic conditions.

Carbon Dioxide under Earth Mantle Conditions: From a Molecular Liquid through a Reactive Fluid to Polymeric Regimes

In both its gaseous and condensed forms, carbon dioxide has an ever-increasing impact on Earth's chemistry and human life and activities. However, many aspects of its high-pressure phase diagram remain unclear. In this work, we present a complete structural characterization of carbon dioxide fluids under geological conditions using extensive ab initio molecular dynamics simulations throughout a wide pressure and temperature range, corresponding to Earth's lower mantle. We identify and describe four different disordered regimes, including two polymeric forms and two molecular ones, all within the geothermal conditions of the lower mantle. At pressures below 40 GPa, we find that the molecular liquid becomes very reactive above 2000 K: the C-O double bond routinely breaks, resulting in small and transient chains composed of CO₂ units and frequently leading to an exchange of oxygen atoms between molecules. At higher pressures, in addition to the polymeric fluid previously reported at 3000 K, we find a polymeric system with glass-like behavior at lower temperatures, suggesting a complex interplay between kinetics and stability.

The phase diagram of carbon dioxide from correlation functions and a many-body potential

The phase stability and equilibria of carbon dioxide are investigated from 125-325 K and 1-10 000 atm using extensive molecular dynamics (MD) simulations and the Two-Phase Thermodynamics (2PT) method. We devise a direct approach for calculating phase diagrams, in general, by considering the separate chemical potentials of the isolated phase at specific points on the P-T diagram. The unique ability of 2PT to accurately and efficiently approximate the entropy and Gibbs energy of liquids allows for assignment of phase boundaries from relatively short (∼100 ps) MD simulations. We validate our approach by calculating the critical properties of the flexible elementary physical model 2, showing good agreement with previous results. We show, however, that the incorrect description of the short-range Pauli force and the lack of molecular charge polarization lead to deviations from experiments at high pressures. We, thus, develop a many-body, fluctuating charge model for CO₂, termed CO₂-Fq, from high level quantum mechanics (QM) calculations that accurately capture the condensed phase vibrational properties of the solid (including the Fermi resonance at 1378 cm^-1) as well as the diffusional properties of the liquid, leading to overall excellent agreement with experiments over the entire phase diagram. This work provides an efficient computational approach for determining phase diagrams of arbitrary systems and underscores the critical role of QM charge reorganization physics in molecular phase stability.

Water-Gas Shift Reaction Produces Formate at Extreme Pressures and Temperatures in Deep Earth Fluids

The water-gas shift reaction is one of the most important reactions in industrial hydrogen production and plays a key role in Fischer-Tropsch-type synthesis, which is widely believed to generate hydrocarbons in the deep carbon cycle but is little known at extreme pressure-temperature conditions found in the Earth's upper mantle. Here, we performed extensive ab initio molecular dynamics simulations and free energy calculations to study the water-gas shift reaction. We found the direct formation of formic acid from CO and supercritical water at 10-13 GPa and 1400 K without any catalyst. Contrary to the common assumption that formic acid or formate is an intermediate product, we found that HCOOH is thermodynamically more stable than the products of the water-gas shift reaction above 3 GPa and at 1000-1400 K. Our study suggests that the water-gas shift reaction may not happen in the Earth's upper mantle, and formic acid or formate may be an important carbon carrier in reducing environments, participating in many geochemical processes in deep Earth.

Equation of state for a chemically dissociative, polyatomic system: Carbon dioxide

A notorious challenge in high-pressure science is to develop an equation of state (EOS) that explicitly treats chemical reactions. For instance, many materials tend to dissociate at high pressures and temperatures where the chemical bonds that hold them together break down. We present an EOS for carbon dioxide (CO₂) that allows for dissociation and captures the key material behavior in a wide range of pressure-temperature conditions. Carbon dioxide is an ideal prototype for the development of a wide-ranging EOS that allows for chemical-dissociation equilibria since it is one of the simplest polyatomic systems and because it is of great interest in planetary science and in the study of detonations. Here, we show that taking dissociation into account significantly improves the accuracy of the resulting EOS compared to other EOSs that either neglect chemistry completely or treat CO₂ dissociation in a more rudimentary way.

Two-dimensional dry ices with rich polymorphic and polyamorphic phase behavior

Both carbon dioxide (CO2) and water (H2O) are triatomic molecules that are ubiquitous in nature, and both are among the five most abundant gases in the Earth's atmosphere. At low temperature and ambient pressure, both CO2 and H2O form molecular crystals--dry ice I and ice I h Because water possesses distinctive hydrogen bonds, it exhibits intricate and highly pressure-dependent phase behavior, including at least 17 crystalline ice phases and three amorphous ice phases. In contrast, due to its weak van der Waals intermolecular interactions, CO2 exhibits fewer crystalline phases except at extremely high pressures, where nonmolecular ordered structures arise. Herein, we show the molecular dynamics simulation results of numerous 2D polymorphs of CO2 molecules in slit nanopores. Unlike bulk polymorphs of CO2, 2D CO2 polymorphs exhibit myriad crystalline and amorphous structures, showing remarkable polymorphism and polyamorphism. We also show that depending on the thermodynamic path, 2D solid-to-solid phase transitions can give rise to previously unreported structures, e.g., wave-like amorphous CO2 structures. Our simulation also suggests intriguing structural connections between 2D and 3D dry ice phases (e.g., Cmca and PA-3) and offers insights into CO2 polyamorphic transitions through intermediate liquid or amorphous phases.

Creating new layered structures at high pressures: SiS2

Old and novel layered structures are attracting increasing attention for their physical, electronic, and frictional properties. SiS₂, isoelectronic to SiO₂, CO₂ and CS₂, is a material whose phases known experimentally up to 6 GPa exhibit 1D chain-like, 2D layered and 3D tetrahedral structures. We present highly predictive ab initio calculations combined with evolutionary structure search and molecular dynamics simulations of the structural and electronic evolution of SiS₂ up to 100 GPa. A highly stable CdI2-type layered structure, which is octahedrally coordinated with space group surprisingly appears between 4 and up to at least 100 GPa. The tetrahedral-octahedral switch is naturally expected upon compression, unlike the layered character realized here by edge-sharing SiS₆ octahedral units connecting within but not among sheets. The predicted phase is semiconducting with an indirect band gap of about 2 eV at 10 GPa, decreasing under pressure until metallization around 40 GPa. The robustness of the layered phase suggests possible recovery at ambient pressure, where calculated phonon spectra indicate dynamical stability. Even a single monolayer is found to be dynamically stable in isolation, suggesting that it could possibly be sheared or exfoliated from bulk -SiS₂.

Crystal structures and dynamical properties of dense CO2

Structural polymorphism in dense carbon dioxide (CO₂) has attracted significant attention in high-pressure physics and chemistry for the past two decades. Here, we have performed high-pressure experiments and first-principles theoretical calculations to investigate the stability, structure, and dynamical properties of dense CO₂ We found evidence that CO₂-V with the 4-coordinated extended structure can be quenched to ambient pressure below 200 K-the melting temperature of CO₂-I. CO₂-V is a fully coordinated structure formed from a molecular solid at high pressure and recovered at ambient pressure. Apart from confirming the metastability of CO₂-V (I-42d) at ambient pressure at low temperature, results of ab initio molecular dynamics and metadynamics (MD) simulations provided insights into the transformation processes and structural relationship from the molecular to the extended phases. In addition, the simulation also predicted a phase V'(Pna2₁) in the stability region of CO₂-V with a diffraction pattern similar to that previously assigned to the CO₂-V (P2₁2₁2₁) structure. Both CO₂-V and -V' are predicted to be recoverable and hard with a Vicker hardness of ∼20 GPa. Significantly, MD simulations found that the CO₂ in phase IV exhibits large-amplitude bending motions at finite temperatures and high pressures. This finding helps to explain the discrepancy between earlier predicted static structures and experiments. MD simulations clearly indicate temperature effects are critical to understanding the high-pressure behaviors of dense CO₂ structures-highlighting the significance of chemical kinetics associated with the transformations.

Pressure-induced Transformations of Dense Carbonyl Sulfide to Singly Bonded Amorphous Metallic Solid

The application of pressure, internal or external, transforms molecular solids into non-molecular extended network solids with diverse crystal structures and electronic properties. These transformations can be understood in terms of pressure-induced electron delocalization; however, the governing mechanisms are complex because of strong lattice strains, phase metastability and path dependent phase behaviors. Here, we present the pressure-induced transformations of linear OCS (R3m, Phase I) to bent OCS (Cm, Phase II) at 9 GPa; an amorphous, one-dimensional (1D) polymer at 20 GPa (Phase III); and an extended 3D network above ~35 GPa (Phase IV) that metallizes at ~105 GPa. These results underscore the significance of long-range dipole interactions in dense OCS, leading to an extended molecular alloy that can be considered a chemical intermediate of its two end members, CO₂ and CS₂.

Different evolutionary pathways from B4 to B1 phase in AlN and InN: metadynamics investigations

Pressure-induced B4-B1 phase transitions of AlN and InN at ambient temperature are systematically investigated using density functional-based metadynamics simulations. A homogeneous deformation path, which is energetically favorable, is through a hexagonal structure for AlN, and through a tetragonal structure for InN. Furthermore, the dynamical stability, instead of the mechanical stability, is crucial to determining the phase-transition paths: the intermediate hexagonal structure can remain stable, whereas the tetragonal structure is always unstable. The B4 phase always shows the direct band gap before the occurrence of structure transition, while the band gap of stable intermediate hexagonal phase is indirect for AlN. Finally, the band gap of the ultimate cubic phase is direct for AlN and indirect for InN, due to the strong p-d repulsion at the R point.

Modeling Polymorphic Molecular Crystals with Electronic Structure Theory

Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.

Predicting finite-temperature properties of crystalline carbon dioxide from first principles with quantitative accuracy

The temperature-dependence of the crystalline carbon dioxide (phase I) structure, thermodynamics, and mechanical properties are predicted in excellent agreement with experiment over a 200 K temperature range using high-level electronic structure calculations.

Molecular crystal structures, thermodynamics, and mechanical properties can vary substantially with temperature, and predicting these temperature-dependencies correctly is important for many practical applications in the pharmaceutical industry and other fields. However, most electronic structure predictions of molecular crystal properties neglect temperature and/or thermal expansion, leading to potentially erroneous results. Here, we demonstrate that by combining large basis set second-order Møller–Plesset (MP2) or even coupled cluster singles, doubles, and perturbative triples (CCSD(T)) electronic structure calculations with a quasiharmonic treatment of thermal expansion, experimentally observable properties such as the unit cell volume, heat capacity, enthalpy, entropy, sublimation point and bulk modulus of phase I crystalline carbon dioxide can be predicted in excellent agreement with experiment over a broad range of temperatures. These results point toward a promising future for ab initio prediction of molecular crystal properties at real-world temperatures and pressures.

High-Energy Density and Superhard Nitrogen-Rich B-N Compounds

The pressure-induced transformation of diatomic nitrogen into nonmolecular polymeric phases may produce potentially useful high-energy-density materials. We combine first-principles calculations with structure searching to predict a new class of nitrogen-rich boron nitrides with a stoichiometry of B(3)N(5) that are stable or metastable relative to solid N(2) and h-BN at ambient pressure. The most stable phase at ambient pressure has a layered structure (h-B(3)N(5)) containing hexagonal B(3)N(3) layers sandwiched with intercalated freely rotating N(2) molecules. At 15 GPa, a three-dimensional C222(1) structure with single N-N bonds becomes the most stable. This pressure is much lower than that required for triple-to-single bond transformation in pure solid nitrogen (110 GPa). More importantly, C222(1)-B(3)N(5) is metastable, and can be recovered under ambient conditions. Its energy density of ∼3.44  kJ/g makes it a potential high-energy-density material. In addition, stress-strain calculations estimate a Vicker's hardness of ∼44  GPa. Structure searching reveals a new clathrate sodalitelike BN structure that is metastable under ambient conditions.

Nitrogen Backbone Oligomers

We found that nitrogen and hydrogen directly react at room temperature and pressures of ~35 GPa forming chains of single-bonded nitrogen atom with the rest of the bonds terminated with hydrogen atoms - as identified by IR absorption, Raman, X-ray diffraction experiments and theoretical calculations. At releasing pressures below ~10 GPa, the product transforms into hydrazine. Our findings might open a way for the practical synthesis of these extremely high energetic materials as the formation of nitrogen-hydrogen compounds is favorable already at pressures above 2 GPa according to the calculations.

Perspective: crystal structure prediction at high pressures

Crystal structure prediction at high pressures unbiased by any prior known structure information has recently become a topic of considerable interest. We here present a short overview of recently developed structure prediction methods and propose current challenges for crystal structure prediction. We focus on first-principles crystal structure prediction at high pressures, paying particular attention to novel high pressure structures uncovered by efficient structure prediction methods. Finally, a brief perspective on the outstanding issues that remain to be solved and some directions for future structure prediction researches at high pressure are presented and discussed.

Predicting crystal structures and properties of matter under extreme conditions via quantum mechanics: the pressure is on

The role of quantum mechanical calculations in understanding and predicting the behavior of matter at extreme pressures is discussed in this feature contribution.

Exploring the metallic phase of N2O under high pressure

Using the CALYPSO method, we proposed a new metallic structure of N₂O under high pressure.

Novel structural phases and superconductivity of iridium telluride under high pressures

Transition metal selenide and telluride have recently receive considerable attention due to their possible structural relation to ferropnictide. Pressure is often used as an efficient way to modify the crystal or electronic structure that in many cases lead to new material states of interest. Here we search the structures of IrTe₂ up to 150 GPa using crystal structure prediction techniques combining with ab initio calculations. Three new stable phases under high pressures are predicted, and their electronic structure properties, phonon spectra, and electron-phonon couplings are also investigated. Significant reconstructions of band structures and Fermi surfaces are found in these new phases. Calculated results show that while the C2/m-2 phase has bad metal behavior and very weak electron-phonon coupling, the and I4/mmm phases have relatively higher electron-phonon coupling up to ~ 1.5 and 0.7, respectively. The variable-composition searching have been performed, newly compounds with different stoichiometries, such as IrTe₃, IrTe, and Ir₃Te, are predicted to be thermodynamically and dynamically stable at high pressures. The pressure range investigated here is accessible in the diamond anvil cell experiments, thus our results might stimulate further experimental studies.

Well-tempered metadynamics converges asymptotically

Metadynamics is a versatile and capable enhanced sampling method for the computational study of soft matter materials and biomolecular systems. However, over a decade of application and several attempts to give this adaptive umbrella sampling method a firm theoretical grounding prove that a rigorous convergence analysis is elusive. This Letter describes such an analysis, demonstrating that well-tempered metadynamics converges to the final state it was designed to reach and, therefore, that the simple formulas currently used to interpret the final converged state of tempered metadynamics are correct and exact. The results do not rely on any assumption that the collective variable dynamics are effectively Brownian or any idealizations of the hill deposition function; instead, they suggest new, more permissive criteria for the method to be well behaved. The results apply to tempered metadynamics with or without adaptive Gaussians or boundary corrections and whether the bias is stored approximately on a grid or exactly.

Stable all-nitrogen metallic salt at terapascal pressures

The phase diagram and equation of state of dense nitrogen are of interest in understanding the fundamental physics and chemistry under extreme conditions, including planetary processes, and in discovering new materials. We predict several stable phases of nitrogen at multi-TPa pressures, including a P4/nbm structure consisting of partially charged N(2)(δ+) pairs and N(5)(δ-) tetrahedra, which is stable in the range 2.5-6.8 TPa. This is followed by a modulated layered structure between 6.8 and 12.6 TPa, which also exhibits a significant charge transfer. The P4/nbm metallic nitrogen salt and the modulated structure are stable at high pressures and temperatures, and they exhibit strongly ionic features and charge density distortions, which is unexpected in an element under such extreme conditions and could represent a new class of nitrogen materials. The P-T phase diagram of nitrogen at TPa pressures is investigated using quasiharmonic phonon calculations and ab initio molecular dynamics simulations.

Well-tempered metadynamics as a tool for characterizing multi-component, crystalline molecular machines

The well-tempered, smoothly converging form of the metadynamics algorithm has been implemented in classical molecular dynamics simulations and used to obtain an estimate of the free energy surface explored by the molecular rotations in the plastic crystal, octafluoronaphthalene. The biased simulations explore the full energy surface extremely efficiently, more than 4 orders of magnitude faster than unbiased molecular dynamics runs. The metadynamics collective variables used have also been expanded to include the simultaneous orientations of three neighboring octafluoronaphthalene molecules. Analysis of the resultant three-dimensional free energy surface, which is sampled to a very high degree despite its significant complexity, demonstrates that there are strong correlations between the molecular orientations. Although this correlated motion is of limited applicability in terms of exploiting dynamical motion in octafluoronaphthalene, the approach used is extremely well suited to the investigation of the function of crystalline molecular machines.

Structural evolution of carbon dioxide under high pressure

Using an efficient structure search method based on a particle swarm optimization algorithm, we study the structural evolution of solid carbon dioxide (CO2) under high pressure. Our results show that, although it undertakes many structural transitions under pressure, CO2 is quite resistive to structures with C beyond 4-fold coordination. For the first time, we are able to identify two 6-fold structures of solid CO2 with Pbcn and Pa3 symmetries that become stable at pressures close to 1 TPa. Both structures consist of a network of C-O octahedra, showing hypervalence of the central C atoms. The C-O bond length varies from 1.30 to 1.34 Å at the 4-fold to 6-fold transition, close to the C-O distance in the transition state of a corresponding S(N)2 reaction. It has been a longstanding and challenging objective to stabilize C in a hypervalent state, particularly when it is bonded with nonmetallic elements. Most of the work so far has focused on synthesizing organic molecules with a high coordination number of C. Our results provide a good measure of the resistivity of C toward forming hypervalent compounds with nonmetallic elements and of the barrier of reaction involving C-O bonds.

Physical and chemical transformations of highly compressed carbon dioxide at bond energies

Carbon dioxide exhibits a richness of high-pressure polymorphs with a great diversity in intermolecular interaction, chemical bonding, and crystal structures. It ranges from typical molecular solids to fully extended covalent solids with crystal structures similar to those of SiO2. These extended solids of carbon dioxide are fundamentally new materials exhibiting interesting optical nonlinearity, low compressibility and high energy density. Furthermore, the large disparity in chemical bonding between the extended network and molecular structures results in a broad metastability domain for these phases to room temperature and almost to ambient pressure and thereby offers enhanced opportunities for novel materials developments. Broadly speaking, these molecular-to-non-molecular transitions occur due to electron delocalization manifested as a rapid increase in electron kinetic energy at high density. The detailed mechanisms, however, are more complex with phase metastabilities, path-dependent phases and phase boundaries, and large lattice strains and structural distortions - all of which are controlled by well beyond thermodynamic constraints to chemical kinetics associated with the governing phases and transitions. As a result, the equilibrium phase boundary is difficult to locate precisely (experimentally or theoretically) and is often obscured by the presence of metastable phases (ordered or disordered). This paper will review the pressure-induced transformations observed in highly compressed carbon dioxide and present chemistry perspectives on those molecular-to-non-molecular transformations that can be applied to other low-Z molecular solids at Mbar pressures where the compression energy rivals the chemical bond energies.

Fluid-solid equilibrium of carbon dioxide as obtained from computer simulations of several popular potential models: the role of the quadrupole

In this work the solid-fluid equilibrium for carbon dioxide (CO2) has been evaluated using Monte Carlo simulations. In particular the melting curve of the solid phase denoted as I, or dry ice, was computed for pressures up to 1000 MPa. Four different models, widely used in computer simulations of CO2 were considered in the calculations. All of them are rigid non-polarizable models consisting of three Lennard-Jones interaction sites located on the positions of the atoms of the molecule, plus three partial charges. It will be shown that although these models predict similar vapor-liquid equilibria their predictions for the fluid-solid equilibria are quite different. Thus the prediction of the entire phase diagram is a severe test for any potential model. It has been found that the Transferable Potentials for Phase Equilibria (TraPPE) model yields the best description of the triple point properties and melting curve of carbon dioxide. It is shown that the ability of a certain model to predict the melting curve of carbon dioxide is related to the value of the quadrupole moment of the model. Models with low quadrupole moment tend to yield melting temperatures too low, whereas the model with the highest quadrupole moment yields the best predictions. That reinforces the idea that not only is the quadrupole needed to provide a reasonable description of the properties in the fluid phase, but also it is absolutely necessary to describe the properties of the solid phase.

Room-temperature structures of solid hydrogen at high pressures

By employing first-principles metadynamics simulations, we explore the 300 K structures of solid hydrogen over the pressure range 150-300 GPa. At 200 GPa, we find the ambient-pressure disordered hexagonal close-packed (hcp) phase transited into an insulating partially ordered hcp phase (po-hcp), a mixture of ordered graphene-like H(2) layers and the other layers of weakly coupled, disordered H(2) molecules. Within this phase, hydrogen remains in paired states with creation of shorter intra-molecular bonds, which are responsible for the very high experimental Raman peak above 4000 cm(-1). At 275 GPa, our simulations predicted a transformation from po-hcp into the ordered molecular metallic Cmca phase (4 molecules/cell) that was previously proposed to be stable only above 400 GPa. Gibbs free energy calculations at 300 K confirmed the energetic stabilities of the po-hcp and metallic Cmca phases over all known structures at 220-242 GPa and >242 GPa, respectively. Our simulations highlighted the major role played by temperature in tuning the phase stabilities and provided theoretical support for claimed metallization of solid hydrogen below 300 GPa at 300 K.

Structure of polymeric carbon dioxide CO2-V

The structure of polymeric carbon dioxide (CO2-V) has been solved using synchrotron x-ray powder diffraction, and its evolution followed from 8 to 65 GPa. We compare the experimental results obtained for a 100% CO2 sample and a 1 mol % CO2/He sample. The latter allows us to produce the polymer in a pure form and study its compressibility under hydrostatic conditions. The high quality of the x-ray data enables us to solve the structure directly from experiments. The latter is isomorphic to the β-cristobalite phase of SiO2 with the space group I42d. Carbon and oxygen atoms are arranged in CO4 tetrahedral units linked by oxygen atoms at the corners. The bulk modulus determined under hydrostatic conditions, B0=136(10)  GPa, is much smaller than previously reported. The comparison of our experimental findings with theoretical calculations performed in the present and previous studies shows that density functional theory very well describes polymeric CO2.

Persistence and eventual demise of oxygen molecules at terapascal pressures

Computational searches for structures of solid oxygen under high pressures in the multi-TPa range are carried out using density-functional-theory methods. We find that molecular oxygen persists to about 1.9 TPa at which it transforms into a semiconducting square-spiral-like polymeric structure (I4(1)/acd) with a band gap of ~3.0 eV. Solid oxygen forms a metallic zigzag chainlike structure (Cmcm) at about 3.0 TPa, but the chains in each layer gradually merge as the pressure is increased and a structure of Fmmm symmetry forms at about 9.3 TPa in which each atom has four nearest neighbors. The superconducting properties of molecular oxygen do not vary much with compression, although the structure becomes more symmetric. The electronic properties of oxygen have a complex evolution with pressure, swapping between insulating, semiconducting, and metallic.

High-pressure phases of hydrogen cyanide: formation of hydrogenated carbon nitride polymers and layers and their electronic properties

We have performed a set of first-principles simulations to consider the possible phase transitions in molecular crystals of HCN under high pressure. Our calculations reveal several transition paths from the orthorhombic phase to tetragonal and then to triclinic phases. The transitions from the orthorhombic to the tetragonal phases are of the second order, whereas those from the tetragonal to the triclinic phases turn out to be of the first-order type and characterized by an abrupt decrease in volume. Our calculations show that, by adjustment of the temperature and pressure of the HCN molecular crystal, novel layered and polymeric crystals with insulating, semiconducting or metallic properties can be found. Based on our simulation results, two different crystal formation mechanisms are deduced. The stabilities of the predicted structures at ambient pressure are further assessed by performing phonon or MD simulations. In addition, the electron transport properties of the predicted polymers are obtained using the non-equilibrium Green's function technique combined with density functional theory. The results show that the polymers have metallic-like I-V characteristics.

Three dimensional carbon-nanotube polymers

Eight fascinating sp(2)- and sp(3)-hybridized carbon allotropes have been uncovered using a newly developed ab initio particle-swarm optimization methodology for crystal structure prediction. These crystalline allotropes can be viewed respectively as three-dimensional (3D) polymers of (4,0), (5,0), (7,0), (8,0), (9,0), (3,3), (4,4), and (6,6) carbon nanotubes, termed 3D-(n, 0) or 3D-(n, n) carbons. The ground-state energy calculations show that the carbons all have lower energies than C(60) fullerene, and some are energetically more stable than the van der Waals packing configurations of their nanotube parents. Owing to their unique configurations, they have distinctive electronic properties, high Young's moduli, high tensile strength, ultrahigh hardness, good ductility, and low density, and may be potentially applied to a variety of needs.

Using first-principles metadynamics simulation to predict new phases and probe the phase transition of NaAlH4

A first-principles-based metadynamics technique was used to study the solid-solid phase transition of NaAlH4 under normal pressure and at temperatures close to the melting point in order to understand the mechanism of hydrogen release from NaAlH4. Two orthorhombic phases with space symmetry groups of Cmcm and Pbcm were determined based on the structures isolated from the metadynamics trajectories starting from the α phase (I4(1/a)) of NaAlH(4). The I4(1/a)→Cmcm→Pbcm transitions can be correlated with the symmetry-lowering process (D(2h) → C(2v) → C(s)) of the AlH(4)(-) complex ions. Thermodynamics analysis revealed that the Cmcm phase forms from the I4(1/a) phase at around the melting point (548 K). Our results are in qualitative agreement with a previous Raman scattering measurement, which showed a transition at ∼ 503 K. The Cmcm phase is expected to be a key intermediate in hydrogen desorption prior to Na(3)AlH(6) formation.

Silicon carbonate phase formed from carbon dioxide and silica under pressure

The discovery of nonmolecular carbon dioxide under high-pressure conditions shows that there are remarkable analogies between this important substance and other group IV oxides. A natural and long-standing question is whether compounds between CO(2) and SiO(2) are possible. Under ambient conditions, CO(2) and SiO(2) are thermodynamically stable and do not react with each other. We show that reactions occur at high pressures indicating that silica can behave in a manner similar to ionic metal oxides that form carbonates at room pressure. A silicon carbonate phase was synthesized by reacting silicalite, a microporous SiO(2) zeolite, and molecular CO(2) that fills the pores, in diamond anvil cells at 18-26 GPa and 600-980 K; the compound was then temperature quenched. The material was characterized by Raman and IR spectroscopy, and synchrotron X-ray diffraction. The experiments reveal unique oxide chemistry at high pressures and the potential for synthesis of a class of previously uncharacterized materials. There are also potential implications for CO(2) segregation in planetary interiors and for CO(2) storage.