Benzene under high pressure: a story of molecular crystals transforming to saturated networks, with a possible intermediate metallic phase

In a theoretical study, benzene is compressed up to 300 GPa. The transformations found between molecular phases generally match the experimental findings in the moderate pressure regime (<20 GPa): phase I (Pbca) is found to be stable up to 4 GPa, while phase II (P4(3)2(1)2) is preferred in a narrow pressure range of 4-7 GPa. Phase III (P2(1)/c) is at lowest enthalpy at higher pressures. Above 50 GPa, phase V (P2(1) at 0 GPa; P2(1)/c at high pressure) comes into play, slightly more stable than phase III in the range of 50-80 GP, but unstable to rearrangement to a saturated, four-coordinate (at C), one-dimensional polymer. Actually, throughout the entire pressure range, crystals of graphane possess lower enthalpy than molecular benzene structures; a simple thermochemical argument is given for why this is so. In several of the benzene phases there nevertheless are substantial barriers to rearranging the molecules to a saturated polymer, especially at low temperatures. Even at room temperature these barriers should allow one to study the effect of pressure on the metastable molecular phases. Molecular phase III (P2(1)/c) is one such; it remains metastable to higher pressures up to ～200 GPa, at which point it too rearranges spontaneously to a saturated, tetracoordinate CH polymer. At 300 K the isomerization transition occurs at a lower pressure. Nevertheless, there may be a narrow region of pressure, between P = 180 and 200 GPa, where one could find a metallic, molecular benzene state. We explore several lower dimensional models for such a metallic benzene. We also probe the possible first steps in a localized, nucleated benzene polymerization by studying the dimerization of benzene molecules. Several new (C(6)H(6))(2) dimers are predicted.     

Ab initio study of phase transition and bulk modulus of NaH

The phase transition of NaH from NaCl- to CsCl-type structure is investigated by an ab initio plane-wave pseudopotential density functional theory method with the norm-conserving pseudopotential scheme in the frame of the generalized gradient approximation correction; the isothermal bulk modulus and its first and second pressure derivatives of the NaCl- and CsCl-type structures under high pressure and temperature are obtained through the quasi-harmonic Debye model. The phase transition obtained from the usual condition of equal enthalpies occurs at the pressure of 32 GPa, which is consistent with the experimental and other calculated values. Through the quasi-harmonic Debye model, in which the phononic effects are considered, the dependences of cell volume V and lattice constant a on temperature T at zero pressure, the isothermal bulk modulus B₀ and its pressure derivatives B₀′and B₀″ on pressure P along isotherms 0, 300, and 600 K, are also successfully obtained.     

Pressure-induced phase transformations in LiAlH4

The pressure-induced phase transformations in pure LiAlH4 have been studied using in situ Raman spectroscopy up to 7 GPa. The analyses of Raman spectra reveal a phase transition at approximately 3 GPa from the ambient pressure monoclinic alpha-LiAlH4 phase (P2(1)/c) to a high pressure phase (beta-LiAlH4, reported recently to be monoclinic with space group I4(1)/b) having a distorted [AlH4]- tetrahedron. The Al-H stretching mode softens and shifts dramatically to lower frequencies beyond the phase transformation pressure. The high pressure beta-LiAlH4 phase was pressure quenchable and can be recovered at lower pressures ( approximately 1.2 GPa). The Al-H stretching mode in the quenched state further shifts to lower frequencies, suggesting a weakening of the Al-H bond.     

Structural, electronic, and dynamical properties of methane under high pressure

The electronic structure and lattice dynamical properties of solid methane under high pressure have been studied based on density functional theory. We identify a cubic structure with space group of I43m below 14 GPa, the Pmn2(1) structure in the range of 14-21 GPa, and the P2(1)/c structure from 21 to 65 GPa. Our obtained Raman spectra of the P2(1)/c structure agree well with the typical Raman active modes in the available experimental data. At 65 GPa, methane undergoes a phase transition from P2(1)/c to Pnma. The structures with P2(1)/c and Pnma symmetries are insulating, and under any pressure studied methane always remains in molecular form. For Pnma phase, the orientational ordering of CH(4) molecules varies significantly at 79, 88, and 92 GPa, and by further increasing pressure the rotation of the molecules freezes and orientational ordering remains unchanged.     

Polar symmetry in new high-pressure phases of chloroform and bromoform

Polar ordering has been induced by pressure in solid chloroform (trichloromethane), CHCl3, and bromoform (tribromomethane), CHBr3, obtained by isochoric and isothermal freezing in a diamond anvil cell. Structures of these new polymorphs have been determined by single-crystal X-ray diffraction, CHCl3 at 0.62 and 0.75 GPa and CHBr3 at 0.20 and 0.35 GPa. Despite different centrosymmetric structures of all low-temperature phases of CHCl3 (space group Pbcn) and CHBr3 (P6(3)/m, P1, and P3), the high-pressure phases are isostructural in space group P6(3). The polar phase of CHBr3 is formed at 295 K, already at the freezing pressure of approximately 0.1 GPa, while CHCl3 transforms from the Pbcn phase into the P6(3) phase between 0.62 and 0.75 GPa. It has been demonstrated that the electrostatic contribution to halogen...halogen and H...halogen interactions in the CHCl3 and CHBr3 molecular crystals is favorable for the polar aggregation and that this effect intensifies with increasing pressure.     

Theoretical calculations of high-pressure phases of NiF<Subscript>2</Subscript>: An ab initio constant-pressure study

We have studied the structural properties of the antiferromagnetic NiF₂ tetragonal structure with P4₂/mnm symmetry using density functional theory (DFT) under rapid hydrostatic pressure up to 400 GPa. For the exchange correlation energy we used the local density approximation (LDA) of Ceperley and Alder (CA). Two phase transformations are successfully observed through the simulations. The structures of XF₂-type compounds crystallize in rutile-type structure. NiF₂ undergoes phase transformations from the tetragonal rutile-type structure with space group P4₂/mnm to orthorhombic CaCl₂-type structure with space group Pnnm and from this orthorhombic phase to monoclinic structure with space group C2/m at 152 GPa and 360 GPa, respectively. These phase changes are also studied by total energy and enthalpy calculations. According to these calculations, we perdict these phase transformations at about 1.85 and 30 GPa.     

First Principles Investigation of Electronic Property and High Pressure Phase Stability of SmN

An investigation of electronic property and high pressure phase stability of SmN has been conducted using first principles calculations based on density functional theory. The electronic properties of SmN show a striking feature of a half metal, the majority-spin electrons are metallic and the minority-spin electrons are semiconducting. It was found that SmN undergoes a pressure-induced phase transition from NaCl-type (B1) to CsCl-type structure (B2) at 117 GPa. The elastic constants of SmN satisfy Born conditions at ambient pressure, indicating that B1 phase of SmN is mechanically stable at 0 GPa. The result of phonon spectra shows that B1 structure is dynamically stable at ambient pressure, which agrees with the conclusion derived from the elastic constants.     

Pressure-induced formation of molecular B2X4(μ-X)2 (X = Cl, Br, I) species

Boron(III) halides (BX(3), where X = F, Cl, Br, I) at ambient pressure conditions exist as strictly monomeric, trigonal-planar molecules. Using correlated ab initio calculations, the three heavier halides (X = Cl, Br, I) are shown to possess B(2)X(4)(μ-X)(2) local minima, isostructural with the diborane molecule. The calculated dissociation barrier of the B(2)I(4)(μ-I)(2) species [≈14 kJ/mol with CCSD(T)/cc-pVTZ] may be high enough to allow cryogenic isolation. The remaining dimer structures are more labile, with dissociation barriers of less than 6 kJ/mol. All three dimer species may be stabilized by application of external pressure. Periodic density functional theory calculations predict a new dimer-based P1 solid, which becomes more stable than the P6(3)/m monomer-derived solids at 5 (X = I) to 15 (X = Cl) GPa. Metadynamics simulations suggest that B(2)X(4)(μ-X)(2)-based solids are the kinetically preferred product of pressurization of the P6(3)/m solid.     

Crystal Structure and Non-Hydrostatic Stress-Induced Phase Transition of Urotropine Under High Pressure

High-pressure behavior of hexamethylenetetramine (urotropine) was studied in situ using angle-dispersive single-crystal synchrotron X-ray diffraction (XRD) and Fourier-transform infrared absorption (FTIR) spectroscopy. Experiments were conducted in various pressure-transmitting media to study the effect of deviatoric stress on phase transformations. Up to 4 GPa significant damping of molecular librations and atomic thermal motion was observed. A first-order phase transition to a tetragonal structure was observed with an onset at approximately 12.5 GPa and characterized by sluggish kinetics and considerable hysteresis upon decompression. However, it occurs only in non-hydrostatic conditions, induced by deviatoric or uniaxial stress in the sample. This behavior finds analogies in similar cubic crystals built of highly symmetric cage-like molecules and may be considered a common feature of such systems. DFT computations were performed to model urotropine equation of state and pressure dependence of vibrational modes. The first successful Hirshfeld atom refinements carried out for high-pressure diffraction data are reported. The refinements yielded more realistic C−H bond lengths than the independent atom model even though the high-pressure diffraction data are incomplete.     

High-pressure phase transitions of Cu2O

In the present study, we extensively explored the crystal structures of Cu₂O on increasing the pressure from 0 GPa to 24 GPa using the first-principles density functional calculations. A series of pressure-induced structure phase transitions of Cu₂O are examined. The calculated results show that the phase transitions (Pn-3m phase → R-3m phase → P-3m1 phase) occur at 5 GPa and 12 GPa, respectively. The P-3m1 phase is found to be the metallic phase via band-gap closure under high pressure.     

Pressure-volume relationships for CeSe and CeBe13 up to 25 GPa

The variations in the volumes of CeSe and CeBe₁₃ have been determined by powder X-ray diffraction in a diamond anvil cell up to 25 GPa at room temperature. In each case, the bulk modulus and its first pressure derivative have been determined. They do not indicate the presence of any negative contribution due to a Kondo volume-dependent interaction or a pressureinduced valence change. A crystallographic phase transformation to a cubic CsCl-type structure has been evidenced for CeSe above 18 GPa.     

Pressure-induced phase transitions in PbTiO3: a query for the polarization rotation theory

Our first-principles computations show that the ground state of PbTiO3 under hydrostatic pressure transforms discontinuously from P4mm to R3c at 9 GPa. Spontaneous polarization decreases with increasing pressure so that the R3c phase transforms to the centrosymmetric Rc phase at around 27 GPa. The first-order phase transition between the tetragonal and rhombohedral phases is exceptional since there is no evidence for a bridging phase. The essential feature of the R3c and Rc phases is that they allow the oxygen octahedron to increase its volume VB at the expense of the cuboctahedral volume VA around a Pb ion. This is further supported by the fact that neither the R3m nor Cm phase, which keep the VA/VB ratio constant, is a ground state within the pressure range between 0 and 40 GPa. Thus, tetragonal strain is dominant up to 9 GPa, whereas at higher pressures, efficient compression through oxygen octahedra tilting plays the central role for PbTiO3. Previously predicted pressure induced colossal enhancement of piezoelectricity in PbTiO3 corresponds to unstable Cm and R3m phases. This suggests that the phase instability, in contrast to the polarization rotation, is responsible for the large piezoelectric properties observed in systems like Pb(Zr,Ti)O3 in the vicinity of the morphotropic phase boundary.     

Conformational and phase transformations of chlorocyclohexane at high pressures by Raman spectroscopy

Pressure induced conformational and phase transformations of chlorocyclohexane (CCH) were investigated in a diamond anvil cell by Raman spectroscopy at room temperature. Pure CCH was compressed up to 20 GPa and then decompressed to ambient pressure. The conformational equilibrium was shifted by pressure from equatorial to axial conformers in the fluid phase below 0.7 GPa, consistent with previous observations. Upon further compression, several solid-to-solid phase transitions were identified by the observation of markedly different Raman patterns as well as different pressure dependences of characteristic Raman modes. The possible structures of these phases were analyzed in correlation with previously observed solid phases at low temperatures. Finally, CCH exhibits pressure hysteresis and partial reversibility upon decompression which result in the formation of the phases with different Raman patterns from those obtained upon compression. The difference can be interpreted as conformational contribution as well as the intrinsic plasticity of CCH crystals.     

Compression of silver sulfide: X-ray diffraction measurements and total-energy calculations

Angle-dispersive X-ray diffraction measurements have been performed in acanthite, Ag(2)S, up to 18 GPa in order to investigate its high-pressure structural behavior. They have been complemented by ab initio electronic structure calculations. From our experimental data, we have determined that two different high-pressure phase transitions take place at 5 and 10.5 GPa. The first pressure-induced transition is from the initial anti-PbCl(2)-like monoclinic structure (space group P2(1)/n) to an orthorhombic Ag(2)Se-type structure (space group P2(1)2(1)2(1)). The compressibility of the lattice parameters and the equation of state of both phases have been determined. A second phase transition to a P2(1)/n phase has been found, which is a slight modification of the low-pressure structure (Co(2)Si-related structure). The initial monoclinic phase was fully recovered after decompression. Density functional and, in particular, GGA+U calculations present an overall good agreement with the experimental results in terms of the high-pressure sequence, cell parameters, and their evolution with pressure.     

High-pressure studies of cyclohexane to 40 GPa

We present data from two room temperature synchrotron X-ray powder diffraction studies of cyclohexane up to approximately 40 and approximately 20 GPa. In the first experiment, pressure cycling was employed wherein pressure was varied up to approximately 16 GPa, reduced to 3.5 GPa, and then raised again to 40 GPa. Initially, the sample was found to be in the monoclinic phase (P12(1)/n1) at approximately 8.4 GPa. Beyond this pressure, the sample adopted triclinic unit cell symmetry (P1) which remained so even when the pressure was reduced to 3.5 GPa, indicating significant hysteresis and metastability. In the second experiment, pressure was more slowly varied, and the monoclinic unit cell structure (P12(1)/n1) was observed at lower pressures up to approximately 7 GPa, above which a phase transformation into the P1 triclinic unit cell symmetry occurred. Thus, the pressure onset of the triclinic phase may be dependent upon the pressurizing conditions. High-pressure Raman data that further emphasize a phase transition (probably into phase VI) around 10 GPa are also presented. We also have further evidence for a phase VII, which is probably triclinic.     

High-pressure study of adamantane: variable shape simulations up to 26 GPa

We report simulations of adamantane by carefully combining ab initio and empirical approaches to enable simulations with internal degrees of freedom on crystalline adamantane up to a pressure of 26 GPa. Two sets of simulations, assuming the adamantane molecule as a rigid (RB) and flexible body (FB), have been carried out as a function of pressure up to 26 GPa to understand changes in the crystal structure as well as molecular structure. The flexible body simulations have been performed by including 6 lowest frequency internal modes (obtained from DFT calculations performed with Gaussian98) out of the total of 72. The calculated variation in c/a and V/V(0) from the RB and FB calculations as a function of pressure have been compared with the experimental curve. Other relevant thermodynamic and structural properties reported are the radial distribution functions, deviation in the position of a given type of atom with respect to its position at standard pressure, delta(s), cell parameters, volume, and energy. With an increase in pressure, three additional peaks are seen to develop gradually at three different pressures in the center of mass (com)-com radial distribution function (rdf). We attribute these changes to structural transformations (probably second-order phase transitions) which is consistent with the three phase transitions reported by Vijayakumar et al. for adamantane in the pressure range of 1 atm-15 GPa. Our simulations also show that these additional peaks in the rdf's are associated with the differences between opposite and parallel spin neighbors of Greig and Pawley as well as the crystallographic directional dependence of intermolecular distances in the first three shells of the neighbors. Also, the structural quantities from the RB calculation show considerable deviation from the FB calculation for pressures greater than 5 GPa, which suggests that the rigid body assumption for molecules may not be valid above this pressure.     

Structural phase transition and 5f-electrons localization of PuSe explored by ab initio calculations

An investigation into the structural phase transformation, electronic and optical properties of PuSe under high pressure was conducted by using the full potential linearized augmented plane wave plus local orbitals (FP-LAPW+lo) method, in the presence and in the absence of spin-orbit coupling (SOC). Our results demonstrate that there exists a structural phase transition from rocksalt (B 1) structure to CsCl-type (B 2) structure at the transition pressure of 36.3 GPa (without SOC) and 51.3 GPa (with SOC). The electronic density of states (DOS) for PuSe show that the f-electrons of Pu are more localized and concentrated in a narrow peak near the Fermi level, which is consistent with the experimental studies. The band structure shows that B 1-PuSe is metallic. A pseudogap appears around the Fermi level of the total density of states of B 1 phase PuSe, which may contribute to its stability. The calculated reflectivity R(ω) shows agreement with the available experimental results. Furthermore, the absorption spectrum, refractive index, extinction coefficient, energy-loss spectrum and dielectric function were calculated. The origin of the spectral peaks was interpreted based on the electronic structures.     

Mechanistic aspects of the solid-state transformation of ammonium cyanate to urea at high pressure

The chemical transformation of ammonium cyanate into urea has been of interest to many generations of scientists since its discovery by Friedrich W?hler in 1828. Although widely studied both experimentally and theoretically, several mechanistic aspects of this reaction remain to be understood. In this paper, we apply computational methods to investigate the behavior of ammonium cyanate in the solid state under high pressure, employing a theoretical approach based on the self-consistent-charges density-functional tight-binding method (SCC-DFTB). The ammonium cyanate crystal structure was relaxed under external pressure ranging from 0 to 700 GPa, leading to the identification of five structural phases. Significantly, the phase at highest pressure (above 535 GPa) corresponds to the formation of urea molecules. At ca. 25 GPa, there is a phase transition of ammonium cyanate (from tetragonal P4/nmm to monoclinic P21/m) involving a rearrangement of the ammonium cyanate molecules. This transformation is critical for the subsequent transformation to urea. The crystalline phase of urea obtained above 535 GPa also has P21/m symmetry (Z = 2). This polymorph of urea has never been reported previously. Comparisons to the known (tetragonal) polymorph of urea found experimentally at ambient pressure suggests that the new polymorph is more stable above ca. 8 GPa. Our computational studies show that the transformation of ammonium cyanate into urea is strongly exothermic (enthalpy change -170 kJ mol-1 per formula unit between 530 and 535 GPa). The proposed mechanism for this transformation involves the transfer of two hydrogen atoms of the ammonium cation toward nitrogen atoms of neighboring cyanate anions, and the remaining NH2 group creates a C-NH2 bond with the cyanate unit.     

High-pressure induced conformational and phase transformations of 1,2-dichloroethane probed by Raman spectroscopy

1,2-Dichloroethane (DCE) was loaded into diamond anvil cells and compressed up to 30 GPa at room temperature. Pressure-induced transformations were probed using Raman spectroscopy. At pressures below 0.6 GPa, fluid DCE exists in two conformations, gauche and trans in equilibrium, which is shifted to gauche on compression. DCE transforms to a solid phase with exclusive trans conformation upon further compression. All the characteristic Raman shifts remain constant in fluid phase and move to higher frequencies in the solid phase with increasing pressure. At about 4-5 GPa, DCE transforms from a possible disordered phase into a crystalline phase as evidenced by the observation of several lattice modes and peak narrowing. At 8-9 GPa, dramatic changes in Raman patterns of DCE were observed. The splitting of the C-C-Cl bending mode at 325 cm-1, together with the observation of inactive internal mode at 684 cm-1 as well as new lattice modes indicates another pressure-induced phase transformation. All Raman modes exhibit significant changes in pressure dependence at the transformation pressure. The new phase remains crystalline, but likely with a lower symmetry. The observed transformations are reversible in the entire pressure region upon decompression.     

高压下MgO的热弹性研究

The thermoelastic properties of MgO over a wide range of pressure and temperature are studied using the first-principles plane wave pseudopotential method within the generalized gradient approximation. It is shown that MgO remains in the B1 (NaCl) structure at all pressures existing within the Earth, and transforms into the CsCl-type structure at 397 GPa. The athermal elastic moduli of MgO are calculated, as a function of pressure up to 150 GPa. The calculated results are in excellent agreement with experimental data at zero pressure and compare favorably with other pseudopotential predictions over the pressure regime studied. MgO is found to be highly anisotropic in its elastic properties, with the magnitude of the anisotropy first decreasing between 0 and 20 GPa and then increasing from 20 GPa to 150 GPa. The Cauchy condition is found to be strongly violated in MgO, reflecting the importance of noncentral many-body forces. The thermodynamic properties of MgO are consistent with the experimental data at ambient condition.