Pressure Induced sp2 sp3 Transitions in Carbon Bearing Phases - GFZpublic
←
→
Page content transcription
If your browser does not render page correctly, please read the page content below
1 Pressure‐Induced sp2‐sp3 Transitions in Carbon‐Bearing Phases Sergey S. Lobanov1 and Alexander F. Goncharov2 ABSTRACT Carbon‐bearing phases show a rich variety of structural transitions as an adaptation to pressure. Of particular interest is the crossover from sp2 carbon to sp3 carbon, as physical and chemical properties of carbon in these distinct electronic configurations are very different. In this chapter we review pressure‐induced sp2‐sp3 transi- tions in elemental carbon, carbonates, and hydrocarbons. 1.1. INTRODUCTION Carbon prefers the sp2‐configuration at low pressure, although there are exceptions (e.g. CH4), because of the Understanding of pressure‐induced phase transitions in small carbon atomic radius (70 pm) that is not sufficient carbon‐bearing phases is fundamental for physics and to reduce steric repulsion among the atoms in its chemistry of dense condensed matter and has important coordination. The increase of coordination number upon implications for Earth and planetary sciences. Carbon is the sp2‐sp3 transition always entails longer carbon‐to‐ unique in that it can form compounds with three distinct neighbor bonds in the sp3‐phase to compensate for the electronic configurations (we will focus on sp2 and sp3 con- increase in steric repulsion (Prewitt & Downs, 1998). As a figurations here) depending on the way s and p atomic result, significant atomic rearrangements are necessary for orbitals of carbon interact with each other. The different sp2‐sp3 phase transitions with large energy barriers typical geometrical arrangement of chemical bonds typical of of such crossovers (Powles et al., 2013). High‐temperature compounds with sp2‐bonded and sp3‐bonded carbon (T) conditions are often required to reach the thermody- (Figure 1.1, inset) results in threefold (lower steric repul- namic ground state, such as in the case of graphite and sion) and fourfold (higher steric repulsion) coordination, diamond. Consequently, many carbon‐bearing phases can respectively. High pressure is an efficient driver for sp2‐sp3 be preserved outside their thermodynamic stability range. transitions as it promotes denser crystal structures and This intrinsic metastability of sp3 carbons finds many use- allows negating the steric repulsion of neighboring atoms ful scientific and technological applications. Natural dia- in sp3 phases. More often than not, material properties monds, despite being metastable at ambient conditions, over such transitions change in a radical manner. It is most often secure unique samples of the Earth’s mantle due to obvious in the case of pure carbon where graphite and dia- their chemical inertness (Kopylova et al., 2010; Wirth mond represent archetypal examples of sp2‐bonded and et al., 2014; Pearson et al., 2014; Smith et al., 2016; Nestola sp3‐bonded crystals, which are extremely d issimilar in the et al., 2018; Tschauner et al., 2018). The sluggish kinetics crystal and electronic structure, and hence in the properties of the graphite‐to‐diamond transition also helps to deci- such as, for example, electrical conductivity, elasticity, and pher the complex P‐T history of high‐pressure metamor- vibrational anisotropy (Oganov et al., 2013). phic complexes and mantle xenoliths (De Corte et al., 2000; Hwang et al., 2001; Massonne, 2003; Korsakov 1 GFZ German Research Center for Geosciences, Potsdam, et al., 2010; Mikhailenko et al., 2016; Shchepetova et al., Germany 2017). The rich polymorphism of carbonates driven by 2 Geophysical Laboratory, Carnegie Institution for Science, the diverse carbon bonding patterns presents another Washington, DC, USA example with relevance to geosciences as it creates Carbon in Earth’s Interior, Geophysical Monograph 249, First Edition. Edited by Craig E. Manning, Jung-Fu Lin, and Wendy L. Mao. © 2020 American Geophysical Union. Published 2020 by John Wiley & Sons, Inc. 1
2 CARBON IN EARTH’S INTERIOR 10 has been used since to pinpoint the P‐T conditions of the sp2-carbon sp3-carbon p-orbital sp3-orbital impact. Lonsdaleite has a hexagonal crystal structure fcc that has a slightly higher energy than cubic diamond. As 8 sp2-orbital such, the formation of this phase (e.g. from graphite) sh must follow a certain kinetic route (Erskine & Nellis, Temperature (kK) 6 Melt 1991; Xie et al., 2017; Bundy & Kasper, 1967), making its Hugoniot BC8 occurrence difficult to uniquely identify in both natural Metals 4 Graphite Diamond sp3 and synthetic samples. This resulted in conflicting reports sp3 concerning the microstructure of samples synthesized at sp2 static pressure conditions, questioning the uniqueness of 2 lonsdaleite as a structurally discrete material (Németh SC Postgraphite et al., 2014; Shiell et al., 2016). Further, the existence of a 0 nanometer (nm) scale structural complexity in lonsda- 1 10 100 1000 10000 leite makes its characterization complex even at ambient Pressure (GPa) conditions. As a result of its structural complexities, the mechanism of lonsdaleite formation in early shock wave Figure 1.1 Phase diagram of carbon at extreme P‐T conditions. experiments (e.g. Erskine & Nellis, 1991) was difficult to BC8, sc (simple cubic), sh (simple hexagonal), and fcc (face‐ resolve. In recent years, however, the development of in centered cubic) are theoretically predicted carbon structures at situ X‐ray diffraction (XRD) characterization of dynamic P > 1 TPa that have never been observed in experiment. Inset: processes on a very fast time scale (ps to ns) allowed for The electronic structure of sp2 and sp3‐hybridized carbon. better understanding of the lonsdaleite formation and See electronic version for color representation of the figures in this book. the carbon phase diagram (Rygg et al., 2012; Gupta et al., 2012; Gauthier et al., 2014). Nevertheless, ambiguities still exist. For example, a direct transition from graphite athways for transporting carbon into the deep mantle. p to lonsdaleite at 50 GPa was inferred by Turneaure et al. Aside from the geological importance, the sp2‐sp3 poly- (2017), while Kraus et al. (2016) reported on a different morphism of carbon‐bearing phases is of great techno- sequence that includes a transition to cubic diamond at logical interest as it allows accessing a rich range of diverse 50 GPa and subsequently to lonsdaleite at 170 GPa. physical and chemical properties in an isochemical setting These works stimulate the development of theoretical (Field, 1992). This chapter reviews sp2‐sp3 transitions in models, which also remain contradictory concerning the carbon‐bearing phases that include pure carbon, carbon- mechanisms of graphite‐to‐diamond transformations ates, and hydrocarbons. (Xie et al., 2017; Pineau, 2013; Mundy et al., 2008; Khaliullin et al., 2011). 1.2. ELEMENTAL CARBON As with lonsdaleite, other metastable carbon phases may serve as indicators of shock metamorphism P‐T Graphite and diamond are the only two stable allotropes conditions. Several kinetically accessible carbon mate- of carbon at P < ~1000 GPa (Figure 1.1). However, a variety rials are formed by compression of graphite (and other of metastable phases are known with a broad range of metastable sp2‐bonded carbons) at room temperature. kinetic stability, which is a unique feature of carbon, making Without sufficient thermal energy to activate the its phase diagram very complex, especially at near‐ambient transition to diamond‐like phases, the system relaxes in conditions (e.g. Bundy et al., 1996). In addition, the melting structures that represent a compromise between the curves of the carbon polymorphs remain uncertain in the thermodynamic stimulus and kinetic hindrance. Cold limit of high pressure and high temperature, where obtain- compression (at 300 K) to above 15 GPa, well into the ing direct experimental data is challenging. Here we update stability field of diamond, results in an sp3‐bonded form the previous reviews (Oganov et al., 2013; Bundy et al., 1996; of carbon (termed postgraphite here) as suggested by a Goncharov, 1987; Hazen et al., 2013), mainly concentrating decrease in the electrical conductivity, appearance of on the phase relations at extreme conditions and focusing on optical transparency, and changes in inelastic X‐ray materials of interest to geosciences. scattering spectra (Goncharov et al., 1989; Utsumi & Yagi, 1991; Mao et al., 2003). XRD measurements on 1.2.1. Metastable Phases this sp3‐bearing carbon are not conclusive in the struc- tural determination (Mao et al., 2003; Wang et al., 2012) Lonsdaleite is a naturally occurring diamond polytype because only powder‐like data was examined and the that was first described in Canyon Diablo meteorite peaks of the high‐pressure phase were weak and broad. (Frondel & Marvin, 1967; Hanneman et al., 1967) and Raman spectroscopy, which is generally sensitive to the
PRESSURE‐INDUCED SP 2‐SP 3 TRANSITIONS IN CARBON‐BEARING PHASES 3 type of chemical bonding, is also not decisive as sp3 trons in the core, allowing the valence p‐electrons to be carbon has a two orders of magnitude smaller scattering closer to the nucleus, making stronger directional cross‐section compared to sp2 carbon when the spectra chemical bonds (Yin & Cohen, 1983). As a result, the dia- are excited in the visible spectral range (Ferrari, 2002). mond structure has a much larger P‐T stability range The Raman spectra of cold‐compressed graphite (Wang compared to its Si and Ge sp3 analogs. Unlike Si and Ge, et al., 2012; Goncharov et al., 1990; Xu et al., 2002) show the melting curve of diamond is positive (above the a large broadening of the C‐C stretching mode of the sp2‐ graphite‐diamond‐melt triple point) (Bundy et al., 1996), bonded carbon in the high‐pressure phase. However, and the sp3‐bonded BC8 phase stable at P > ~1 TPa is not positive spectroscopic identification of sp3‐bonded metallic (cf. β‐tin Si) but semiconducting (Correa et al., carbon has not been achieved because the first‐order dia- 2006; Martinez‐Canales et al., 2012). According to theo- mond peak of the diamond anvils and the disorder‐ retical calculations (Wang et al., 2005; Benedict et al., induced D‐band of sp2‐bonded carbon both obscure the 2014), the slope of this diamond melting curve changes spectral range of interest. Overall, the XRD and Raman from positive to negative above 450 GPa due to the soft- data suggest the presence of disorder in postgraphite ening of the transverse acoustic branches near the L and phase(s). The persistence of the sp2 fingerprints in its X symmetry points (Correa et al. 2006). Dynamic com- Raman spectra upon cold compression indicates that pression experiments detected the melting of diamond postgraphite still contains sp2‐bonded carbon, which is close to the expected range of P‐T conditions, inferred a inconsistent with the theoretically predicted M‐carbon crossover in the slope of its melting line, and provided evi- phase (Li et al., 2009) usually considered the best match dence of the second diamond‐BC8‐liquid triple point to the experiment. Interestingly, the transition to post- (Brygoo et al., 2007; Knudson et al., 2008; Eggert et al., graphite is reversible at 300 K (Utsumi & Yagi, 1991; 2009). At higher pressures, theory predicts that the BC8 Mao et al., 2003; Wang et al., 2012), but this phase can be structure undergoes metallization (Correa et al., 2006) or quenched at low temperatures 1000 K measurements did not record discontinuities that could be and P ~20 GPa results in a pressure‐quenchable phase attributed to the predicted phase transition (Smith et al., that is structurally different from postgraphite and shows 2014). In the absence of direct structural characterization, similarities to lonsdaleite (Utsumi & Yagi, 1991). Glassy however, any robust conclusion seems premature. carbon cold‐compressed to 50 GPa and heated to ~1800 K produces an amorphous and pressure‐quenchable sp3‐ 1.3. CARBONATES bonded phase (Zeng et al., 2017; Yin et al., 2011). These results highlight the importance of the initial carbon 1.3.1. (Mg,Fe) Carbonates microstructure on the dense high‐pressure synthesis products, which has the potential for creating new ultra- To the best of our knowledge, Skorodumova et al. light and ultrastrong materials (Hu et al., 2017). (2005) provided the first computational evidence of ther- Theoretical description of all these metastable carbon modynamic stability of sp3‐MgCO3 with a pyroxene phases (see Shi et al., 2018, and references therein) structure (C2/c) at P > 113 GPa with CO42‐ tetrahedral remains largely unsatisfactory. groups arranged in chains. This study, however, explored only the relative stabilities of several possible sp3‐car- 1.2.2. Phase Diagram at High Temperature bonate structures as suggested by crystallographic and and in the TPa Pressure Range chemical similarities with Mg‐bearing silicates. The first unconstrained search for the stable structures in the In the limit of high temperature, the graphite‐to‐diamond MgCO3 system by evolutionary crystal structure predic- transition is well established and sharp (Bundy et al., tion algorithms revealed an sp2‐sp3 transition at 82 GPa 1996). In contrast, the theoretically predicted sp2‐sp3 (Oganov et al., 2008) with a C2/m structure being ~0.2 eV transition in molten carbon is sluggish and the results more favorable than the pyroxene structure initially pro- appear sensitive to the level of theory used in the compu- posed by Skorodumova et al. (2005). This sp3‐MgCO3 tation (e.g. first‐principles vs. classical) (Ghiringhelli phase was termed magnesite‐II (phase II) with a wide sta- et al., 2004). With regard to solid–solid transformations bility field up to 138 GPa where another sp3 modification and melting at very high pressures, carbon is unique in of MgCO3 (P21, phase III) was predicted to take over the group IV of the periodic table in that it has no p‐elec- (Oganov et al., 2008). More structures, however, were
4 CARBON IN EARTH’S INTERIOR proposed in the 90–120 GPa range as energetically com- (Boulard et al., 2011). A crystal structure with P21/c petitive by Panero and Kabbes (2008) based on a first symmetry allowed for improved Le Bail refinement over principles evaluation of several potential structures. the predicted C2/m MgCO3 structure. Pickard and Needs (2015) revised the phase diagram of Similar experimental evidence was put forward in MgCO3 using ab initio random structure searching tech- support of the sp3‐(Fe, Mg)CO3 synthesis from nique to find new stable P‐1 (85‐101 GPa) and P212121 (P (Fe0.75Mg0.25) CO3 or (Mg0.6Fe0.4) + CO2 after laser heating >144 GPa) MgCO3 structures with sp3 carbon in addition to ~2000 K at 80 GPa and to ~2850 K at 105 GPa, respec- to the C2/m structure proposed earlier (Oganov et al., tively (Boulard et al., 2011). Regardless of the starting 2008). Most recently, another polymorph of sp3‐MgCO3 components, the measured XRD patterns showed simi- (P21) was predicted at P > 143 GPa (Yao et al., 2018). larities to that observed in pure MgCO3, and the authors Figure 1.2 summarizes theoretical predictions for the concluded that the same crystal structure formed in these MgCO3 (and other) system(s). two systems. In the iron‐containing system, however, the These first‐principles predictions stimulated the exper- synthesized phase was not only temperature‐quenchable imental search for pressure‐induced sp2‐sp3 transitions in but also pressure‐quenchable, and the recovered samples carbonates. Boulard et al. (2011) attempted to synthesize were analyzed by electron energy loss spectroscopy, which sp3‐MgCO3 in a diamond anvil cell (DAC) using syn- showed a peak at 290.7 eV (assigned to C3O96‐) as opposed chrotron XRD as a probe to detect the new phases. At 80 to 290.3 eV in reference sp2‐siderite (Boulard et al., 2011). GPa and ~2400 K, the authors observed new XRD peaks However, it is unclear whether the authors’ interpretation that could not be attributed to the starting materials is unequivocal as the reported energy resolution (width (MgCO3 or MgO + CO2). This new phase was not tem- of the zero‐loss peak at its half height) was 0.6 eV. perature‐quenchable and only the low‐pressure magne- In a subsequent work, Boulard et al. (2015) reported on site structure could be observed after cooling (Boulard the room‐temperature infrared spectrum of (Fe0.75Mg0.25) et al., 2011). Nonetheless, on the basis of Le Bail–type CO3 that was presynthesized at ~2100 K at 103 GPa. The crystallographic refinement of high‐temperature XRD synthesis of C2/m (Fe,Mg)CO3, as well as of a high‐ patterns, the authors concluded that they successfully pressure polymorph of CaMn2O4‐type Fe3O4, was synthesized magnesite II with a crystal structure based reported based on synchrotron XRD. One should bear in on three CO4‐tetrahedra linked into (C3O9)6‐ rings mind, however, that reliable identification of these low symmetry structures at P ~100 GPa based solely on a Le Bail–type refinement is challenging. For example, the 3.0 identification of the Fe3O4 structure could not have been C2/c Mg2Fe2C4O13 C2/c Fe2Fe2C4O13 possible because iron oxides show a complex high P‐T Pnma behavior (Bykova et al., 2016), the details of which were P21/c CaCO3 Ca(Mg,Fe)(CO ) 3 2 not known at the time of the Boulard et al. (2015) work. Temperature (1000 K) 2.0 The stable structure of Fe3O4 is different from the one Fe2Fe2C4O13 used for Le Bail refinement in Boulard et al. (2015), which + Fe4C3O12 undermines the robustness and the overall conclusions of that work. Merlini et al. (2015) provided a more reliable insight 1.0 into the crystal structures of sp3‐(Fe,Mg) carbonates by C2/c CaMg(CO3)2 P212121 MgCO3 means of single‐crystal structure solution methods find- P21/c CaCO3 C2/m MgCO3 P-1 MgCO3 P21 MgCO3 ing that sp2‐(Fe0.75Mg0.25) CO3 transforms to a C2/c Mg2Fe2(C4O13) together with a C2/m Fe13O19 at P ~135 GPa and T ~2650 K (Figure 1.2). The new carbonate 0.0 50 75 100 125 150 phase contained tetrahedrally coordinated sp3‐hybridized Pressure (GPa) carbon linked into corner‐shared truncated C4O13 chains and apparently remained stable on decompression down Figure 1.2 Pressure‐temperature conditions of sp2‐sp3 phase to 40 GPa. Importantly, this study provided the first evi- transitions in carbonates. Theoretical predictions (calculations dence that the composition of the (Mg,Fe) carbonate at 0 K) are shown by dashed lines. Only the most recent reports may change over the sp2‐sp3 transition, which had not are shown: MgCO3 (Pickard & Needs, 2015; Yao et al., 2018), been considered by first‐principles computations. A sim- CaMg(CO3)2 (Solomatova & Asimow, 2017), and CaCO3 (Pickard & Needs, 2015). Circles depict the synthesis condi- ilar single‐crystal XRD technique was subsequently tions of sp3‐phases in (Mg,Fe)‐ (Merlini et al., 2015), (Ca,Mg,Fe)‐ applied to the FeCO3 endmember uncovering two distinct (Merlini et al., 2017), and Ca‐carbonate (Lobanov et al., 2017) sp3 carbonates at P > ~70 GPa (Cerantola et al., 2017). systems. See electronic version for color representation of the According to this study, Fe4(C3O12) and Fe2Fe2(C4O13) figures in this book. coexist in the temperature limit of ~1500–2200 K, while
PRESSURE‐INDUCED SP 2‐SP 3 TRANSITIONS IN CARBON‐BEARING PHASES 5 only Fe2Fe2(C4O13) is present at T > ~2200 K. The crystal 2015; Yao et al., 2018). Both C2221 and P21/c CaCO3 are structure of Fe4(C3O12) involves isolated CO4 tetrahedra, based on CO4‐tetrahedra linked into 1D pyroxene‐like while Fe2Fe2(C4O13) contains truncated chains of CO4 chains with only slight differences in the arrangement of tetrahedra (Figure 1.3) and is isostructural with CO4‐groups. All vertex‐sharing helices in the C2221 struc- Mg2Fe2(C4O13) found earlier by Merlini et al. (2015). ture are right‐handed, while half helices in the P21/c Overall, the case of (Mg,Fe) carbonates shows that CaCO3 are left‐handed. they exhibit a remarkably complex crystallographic To test these theoretical predictions, Lobanov et al. behavior over the sp2‐sp3 transition. Initial theoretical (2017) used synchrotron XRD and Raman spectroscopy predictions of the stable sp3‐MgCO3 did find some exper- as probes for the sp2‐sp3 transition. Crystallographic data imental support in both Mg carbonates and (Mg,Fe) collected after the heating of CaCO3 to ~2000 K at 105 carbonates. However, subsequent advanced crystallo- GPa provided evidence of the sp2‐sp3 transition but were graphic studies of (Mg,Fe) carbonates and Fe carbonates not sufficient to discriminate between the proposed C2221 unraveled a complex behavior that always involves a and P21/c CaCO3 models. Raman spectra measured change in the carbonate chemical composition, which before and after the heating, however, were clearly dis- was not anticipated from the earlier ab initio computa- tinct. A characteristic C‐O stretching vibration appears tions. For these reasons, earlier experimental reports on after the heating with a frequency that is ~20% lower than sp3‐(Mg,Fe)CO3 (Boulard et al., 2011, 2015, 2012) must that in sp2 CaCO3 prior to the heating (Figure 1.4a), con- be deemed unreliable. sistent with the longer C‐O bond length in CO4 tetra- hedra. This spectral feature was reproduced independently 1.3.2. Ca Carbonates and Ca(Mg,Fe) Carbonates via ab initio computations of the vibrational spectrum of sp3 CaCO3 and provided strong evidence in favor of the The theoretical search for an sp2‐sp3 transition in P21/c model of CaCO3 (Figure 1.4b). Altogether, the CaCO3 using evolutionary crystal structure prediction combined crystallographic and spectroscopic approach algorithms yielded an orthorhombic pyroxene‐type proved successful in detecting the sp2‐sp3 transition in (C2221) structure at P > 137 GPa (Oganov et al., 2006). CaCO3 and provided a reliable reference for the frequency Ono et al. (2007) showed that laser‐heating of CaCO3 at of C‐O stretching in CO4 tetrahedra, which can be used P > 130 GPa results in a new phase that is structurally to identify sp2‐sp3 transitions in other carbonates consistent with the predicted C2221 CaCO3. To the best (Lobanov et al., 2017). of our knowledge, the synthesis of sp3‐CaCO3 by Ono Because of the chemical reactions with mantle min- et al. (2007) is the first experimental report of an sp3 car- erals, CaMg(CO3)2 is a more realistic composition for the bonate. The crystal structure of this phase, however, was mantle carbonate. Merlini et al. (2017) reported on the challenged by more recent computations that found that synthesis of sp3 Ca(Mg0.6,Fe0.4)(CO3)2 (dolomite‐IV) at P21/c CaCO3 is ~0.2 eV/f.u. more energetically favorable 115 GPa and 2500 K. The crystal structure of dolomite‐ than the C2221 model and is stable at P > 76 GPa IV was solved from single‐crystal XRD data and con- (Figure 1.2) and up to at least 160 GPa (Pickard & Needs, tains CO4 tetrahedra linked into threefold C3O9 rings (a) (b) Fe4(C3O12) Fe2Fe2(C4O13) C1 C2 C Fe2 Fe1 Fe1 Fe2 x z y x Figure 1.3 Crystal structures of Fe4(C3O12) and Fe2Fe2(C4O13) after Cerantola, et al. (2017). Oxygen atoms are not shown for clarity. Unit cells are shown by thin black lines. See electronic version for color representation of the figures in this book.
6 CARBON IN EARTH’S INTERIOR (a) (b) CaCO3 at 105 GPa CaCO3 at 105 GPa Raman intensity (arb. units) Raman intensity (arb. units) unheated C2221 – CaCO3 660 nm 532 nm P21 /c – CaCO3 488 nm 488 nm 100 300 500 700 900 1100 1300 100 300 500 700 900 1100 1300 Raman shift (cm–1) Raman shift (cm–1) Figure 1.4 (a) Raman spectra of CaCO3 at 105 GPa collected with 488, 532, and 660 nm excitations. Vertical bars are computed Raman modes of P21/c‐CaCO3 (bottom) and post‐aragonite CaCO3 (top). (b) Experimental spectrum of CaCO3 laser‐heated at 105 GPa in comparison with the theoretical spectra of P21/c and C2221‐CaCO3 at 105 GPa as computed by LDA‐DFT. Shaded areas are guides to compare the computed spectra with experiment. Adapted from Lobanov et al., 2017. See electronic version for color representation of the figures in this book. (a) (b) sp3 Ca(Mg,Fe)(CO3)2 (experiment) sp3 Ca(Mg,Fe)(CO3)2 (theory) C Ca z C x Ca x Mg Figure 1.5 (a) Crystal structure of sp3‐Ca(Mg0.6,Fe0.4)(CO3)2 (dolomite‐IV) after Merlini, et al. (2017). Mg and Fe atoms are not shown as they occupy identical crystallographic positions to Ca. (b) Theoretically predicted crystal structure of sp3‐CaMg(CO3)2 (Solomatova & Asimow, 2017). In a and b, oxygen atoms are not shown for clarity. Unit cells are shown by thin lines. See electronic version for color representation of the figures in this book. (Figure 1.5a). Interestingly, Ca, Mg, and Fe cations agnesite‐II phase (Oganov et al., 2008). However, an m occupy identical positions in the lattice, suggesting that evolutionary search for sp3‐CaMg(CO3)2 yielded a differ- moderate addition of iron does not bear important effects ent motif with CO4‐tetrahedra linked into pyroxene‐like on the sp3‐CaMg(CO3)2 crystal structure. The presence chains (Solomatova & Asimow, 2017), similarly to what of C3O9 rings in the dolomite‐IV structure (Merlini has been predicted and experimentally confirmed for sp3‐ et al., 2017) is an important similarity to the predicted CaCO3. Interestingly, Raman spectra of dolomite‐III at
PRESSURE‐INDUCED SP 2‐SP 3 TRANSITIONS IN CARBON‐BEARING PHASES 7 P > 63 GPa, another polymorph of CaMg(CO3)2 stable at examples of pressure‐induced polymerization in unsatu- 36–115 GPa (Merlini et al., 2017), show a new band appear- rated hydrocarbons have been reported in the literature ing as a low‐frequency shoulder of the C‐O symmetric (see Schettino & Bini, 2007, for a more thorough review). stretching in CO3 groups (Vennari & Williams, 2018). This Below we provide only a handful of these. High‐density shoulder was interpreted to be due to an increase in carbon polyethylene with a large fraction of sp3 carbon can be coordination (3+1); thus, it may indicate a sluggish sp2‐sp3 obtained by compressing ethylene (C2H4) in a DAC at transition in CaMg(CO3)2 at room temperature. room temperature to ~3.6 GPa (Chelazzi et al., 2004). In summary, CaCO3 is the only sp3 carbonate for which Solid acetylene polymerizes at 3.5 GPa with the trans‐ theoretical and experimental data on the structure of the opening of the triple bond (Aoki et al., 1988; Santoro sp3‐bonded phase are in agreement. The crystal structures et al., 2003). Benzene, a prototype of aromatic com- of more complex sp3 carbonates, such as in the (Mg,Fe) pounds, polymerizes at P > 23 GPa (Prizan et al., 1990) and (Ca,Mg) systems, have been solved directly from with the formation of a one‐dimensional sp3 carbon single crystal XRD data, but the agreement between these nanomaterial (Fitzgibbons et al., 2015). Although the structures and those from theoretically predicted models mechanism of benzene polymerization remains contro- is yet to be seen. Since the experimentally observed struc- versial (Wen et al., 2011), it is generally agreed that the tures are chemically complex and contain a large number process is governed by the pressure‐activated π‐π electron of atoms per cell, it is possible that earlier theoretical sim- density overlap of the neighboring molecules, as initially ulations on chemically simple carbonate systems did not suggested by Drickamer (1967). Ciabini et al. (2007) sample unit cells with a sufficient number of atoms. On showed that benzene polymerization is triggered at a crit- the other hand, it is possible that the experimentally ical C‐C distance of ~2.6 Å between adjacent C6H6 mole- observed phases are not representative of the ground cules. These experimentally observed polymerization state sp3 structures because of severe chemical and reactions have also been reproduced by ab initio compu- thermal gradients that are typical of laser‐heated DACs, tations and confirm the importance of intermolecular especially when iron is present in the system. Subsequent separation and molecular orientation (Ciabini et al., studies employing multiple characterization techniques 2007; Bernasconi et al., 1997). will be of significant value to the topic of sp2‐sp3 transi- The conjugated nature of the recovered products may tions in carbonates. give rise to peculiar physical properties. For example, a very high electrical conductivity (~100 Ω‐1 cm‐1) was 1.4. HYDROCARBONS observed for polyacetylene (Clark & Hester, 1991). Hydrocarbons are expected to form polymeric states not High‐pressure behavior of the carbon‐hydrogen sys- only at relatively low pressures and room temperature but tems is relevant to planetary sciences as hydrocarbons are also at P‐T conditions representative of planetary inte- an important mantle component of icy giant planets riors (Lee & Scandolo, 2011; Chau et al., 2011; Sherman (Hubbard, 1981). Application of high external pressure et al., 2012; Lobanov et al., 2013) where polymerized to hydrocarbons, which are typically soft molecular crys- hydrocarbons may promote the high electrical conduc- tals, forces the molecules to shorten intermolecular tivity expected of the icy mantles of Uranus and Neptune distances, making intermolecular and intramolecular (Stanley & Bloxham, 2004). As such, knowledge of the forces comparable (Hemley & Dera, 2000). At short conditions of hydrocarbon polymerization at extreme enough intermolecular distances, molecules of unsatu- pressure and temperature is critical for our understanding rated hydrocarbons (i.e. with π‐bonds in the chemical of the icy giant planet interiors. structure) are forced to react with each other (polymerize) Many theoretical and experimental studies have with the formation of conjugated molecular structures. explored the high P‐T stability of hydrocarbons. Polymerization reactions in unsaturated hydrocarbons Thermodynamic and molecular dynamics modeling of involve breaking of unsaturated π‐bonds with the hydrocarbon systems showed that conjugated hydrocar- formation of σ‐bonds between the adjacent molecules. As bons (molecular and polymeric) are more stable at high such, pressure‐induced polymerization of unsaturated pressure than at 1 atm (Sherman et al., 2012; Kenney hydrocarbons can be viewed as an sp2‐sp3 phase transition et al., 2002; Zhang & Duan, 2009; Spanu et al., 2011; Gao (Drickamer, 1967). This process is typically irreversible et al., 2010). This is consistent with experimental reports and hydrocarbon polymers are available for ex situ char- on hydrocarbon synthesis from methane in laser‐heated acterization after decompression. DACs (Lobanov et al., 2013; Hirai et al., 2009; Kolesnikov Pressure‐induced polymerization of unsaturated et al., 2009). Interestingly, unsaturated hydrocarbons can hydrocarbons has been extensively characterized by be formed from methane at P > 24 GPa and T > 1500 K Raman and IR spectroscopy, which detected signatures (Lobanov et al., 2013), as observed by Raman spectroscopy of sp3‐bonded chemical structures at high pressure. Many at 300 K after laser heating. This, however, is inconsistent
8 CARBON IN EARTH’S INTERIOR with the apparent instability of unsaturated carbon‐carbon Drickamer, H. G. (1967). Science, 156, 1183. double and triple bonds at high pressure (as discussed Eggert, J. H., Hicks, D. G., Celliers, P. M., Bradley, D. K., above). It is likely that Raman spectroscopy at 300 K pro- McWilliams, R. S., Jeanloz, R., Miller, J. E., Boehly, T. R., & vides a distorted snapshot of the hydrocarbon fluid at high Collins, G. W. (2009). Nat. Phys., 6, 40. Erskine, D. J., & Nellis, W. J. (1991). Nature, 349, 317. temperature and that unsaturated carbon‐carbon bonds Ferrari, A. C. (2002). Diamond Relat. Mater., 11, 1053. form upon quenching of a highly dissociated fluid. Recent Field, J. E. (1992). The properties of natural and synthetic dia- laser‐shock experiments on hydrocarbons provided an in mond. San Diego: Academic Press. situ characterization of hydrocarbon fluids at high P‐T and Fitzgibbons, T. C., Guthrie, M., Xu, E. S., Crespi, V. H., confirmed that such fluids are nonmolecular at P < 100– Davidowski, S. K., Cody, G. D., Alem, N., & Badding, J. V. 150 GPa and T < ~3000 K (Kraus et al., 2017). Further, this (2015). Nat. Mater., 14, 43. nonmolecular hydrocarbon fluid undergoes dissociation Frondel, C., & Marvin, U. B. (1967). Nature, 214, 587. into diamond and hydrogen at P > 150 GPa and T > 3000 Gao, G. Y., Oganov, A. R., Ma, Y. M., Wang, H., Li, P. F., Li, K. As such, methane dissociation to diamond and hydrogen Y. W., Iitaka, T., & Zou, G. T. (2010). J. Chem. Phys., 133, at high P‐T can be viewed as an sp3‐sp3 transition. 144508. Ghiringhelli, L. M., Los, J. H., Meijer, E. J., Fasolino, A., & Frenkel, D. (2004). Phys. Rev. B, 69, 100101. ACKNOWLEDGMENTS Gloahec, M. (2007). Nat. Mater., 6, 274. Goncharov, A. F. (1987). Sov. Phys. Uspekhi, 30, 525. This work was supported by the Deep Carbon Goncharov, A. F., Makarenko, I. N., & Stishov, S. M. (1989). Observatory. SSL acknowledges the support of the Zh. Eksper. Teor. Fiziki, 96, 670. Helmholtz Young Investigators Group CLEAR (VH‐ Gauthier, M. et al. (2014). Rev. Sci. Instrum., 85, 11E616. NG‐1325) and Deutsche Forschungsgemeinschaft within Goncharov, A. F., Makarenko, I. N., & Stishov, S. M. (1990). the research group FOR2125 CarboPaT. High Pressure Res., 4, 345. Gupta, Y. M., Turneaure, S. J., Perkins, K., Zimmerman, K., Arganbright, N., Shen, G., and Chow, P. (2012). Rev. Sci. REFERENCES Instrum., 83, 123905. Hanneman, R. E., Strong, H. M., & Bundy, F. P. (1967). Science, Aoki, K., Usuba, S., Yoshida, M., Kakudate, Y., Tanaka, K., & 155, 995. Fujiwara, S. (1988). J. Chem. Phys., 89, 529. Hazen, R. M., Downs, R. T., Jones, A. P., & Kah, L. (2013). Benedict, L. X., Driver, K. P., Hamel, S., Militzer, B., Qi, T. T., Correa, Rev. Mineral. Geochem., 75, 7. A. A., Saul, A., & Schwegler, E. (2014). Phys. Rev. B, 89, 224109. Hemley, R. J., & Dera, P. (2000). Rev. Mineral. Geochem., 41, 335. Bernasconi, M., Chiarotti, G. L., Focher, P., Parrinello, M., & Hirai, H., Konagai, K., Kawamura, T., Yamamoto, Y., & Yagi, Tosatti, E. (1997). Phys. Rev. Lett., 78, 2008. T. (2009). Phys. Earth Planet. Inter., 174, 242. Boulard, E., et al. (2012). Journal of Geophysical Research, 117, Hu, M. et al. (2017). Sci. Adv., 3, e1603213. B02208. Hubbard, W. B. (1981). Science, 214, 145. Boulard, E., Gloter, A., Corgne, A., Antonangeli, D., Auzende, Hwang, S. L., Shen, P., Chu, H. T., Yui, T. F., & Lin, C. C. A. L., Perrillat, J. P., Guyot, F., & Fiquet, G. (2011). Proc. (2001). Earth Planet. Sci. Lett., 188, 9. Natl. Acad. Sci. U.S.A., 108, 5184. Kenney, J. F., Kutcherov, V. A., Bendeliani, N. A., & Alekseev, Boulard, E., Pan, D., Galli, G., Liu, Z. X., & Mao, W. L. (2015). V. A. (2002). Proc. Natl. Acad. Sci. U.S.A., 99, 10976. Nature Comm., 6, 6311. Khaliullin, R. Z., Eshet, H., Kühne, T. D., Behler, J., & Brygoo, S., Henry, E., Loubeyre, P., Eggert, J., Koenig, M., Parrinello, M. (2011). Nat. Mater., 10, 693. Loupias, B., Benuzzi‐Mounaix, A., & Rabec Le Bundy, F. P., Knudson, M. D., Desjarlais, M. P., & Dolan, D. H. (2008). Bassett, W. A., Weathers, M. S., Hemley, R. J., Mao, H. U., & Science, 322, 1822. Goncharov, A. F. (1996). Carbon, 34, 141. Kolesnikov, A., Kutcherov, V. G., & Goncharov, A. F. (2009). Bundy, F. P., & Kasper, J. S. (1967). J. Chem. Phys. 46, 3437. Nature Geosci., 2, 566. Bykova, E., et al. (2016). Nat. Commun., 7, 10661. Kopylova, M., Navon, O., Dubrovinsky, L., & Khachatryan, Cerantola, V., et al. (2017). Nat. Commun., 8, 15960. (2010). G. Earth Planet. Sci. Lett., 291, 126. Chau, R., Hamel, S., & Nellis, W. J. (2011). Nat. Commun., 2, 1198. Korsakov, A. V., Perraki, M., Zedgenizov, D. A., Bindi, L., Chelazzi, D., Ceppatelli, M., Santoro, M., Bini, R., & Schettino, Vandenabeele, P., Suzuki, A., & Kagi, H. (2010). J. Petrol., 51, 763. V. (2004). Nat. Mater., 3, 470. Kraus, D., et al. (2016). Nat. Commun., 7, 10970. Ciabini, L., Santoro, M., Gorelli, F. A., Bini, R., Schettino, V., Kraus, D., et al. (2017). Nat. Astron., 1, 606. & Raugei, S. (2007). Nat. Mater., 6, 39. Lee, M. S., & Scandolo, S. (2011). Nat. Commun., 2, 1184. Clark, R. J. H., & Hester, R. E. Spectroscopy of advanced mate- Li, Q., Ma, Y., Oganov, A. R., Wang, H., Wang, H., Xu, Y., Cui, rials. Chichester: Wiley, Advances in Spectroscopy, 19. T., Mao, H.‐K., & Zou, G. T. (2009). Phys. Rev. Lett., 102, Correa, A. A., Bonev, S. A., & Galli, G. (2006). Proc. Natl. 175506. Acad. Sci. U.S.A., 103, 1204. Lin, Y., Zhang, L., Mao, H.‐K., Chow, P., Xiao, Y., Baldini, M., De Corte, K., Korsakov, A., Taylor, W. R., Cartigny, P., Ader, Shu, J., & Mao, W. L. (2011). Phys. Rev. Lett., 107, 175504. M., & De Paepe, P. (2000). Isl. Arc, 9, 428. Lobanov, S. S., et al. (2017). Phys. Rev. B, 96, 104101.
PRESSURE‐INDUCED SP 2‐SP 3 TRANSITIONS IN CARBON‐BEARING PHASES 9 Lobanov, S. S., Chen, P. N., Chen, X. J., Zha, C. S., Litasov, K. Schettino, V., & Bini, R. (2007). Chem. Soc. Rev., 36, 869. D., Mao, H. K., & Goncharov A. F., (2013). Nature Comm., Shchepetova, O. V., Korsakov, A., Mikhailenko, D., Zelenovskiy, 4, 2446. P., Shur, V., & Ohfuji, H. (2017). J. Raman Spectrosc., 48, 1606. Mao W. L., et al. (2003). Science, 302, 425. Sherman, B. L., Wilson, H. F., Weeraratne, D., & Militzer, B. Martinez‐Canales, M., Pickard, C. J., & Needs, R. J. (2012). (2012). Phys. Rev. B, 86, 224113. Phys. Rev. Lett., 108, 045704. Shi, X., He, C. R., Pickard, C. J., Tang, C. C., & Zhong, J. Massonne, H. J. (2003). Earth Planet. Sci. Lett., 216, 347. (2018). Phys. Rev. B, 97, 014104. Merlini, M., Cerantola, V., Gatta, G. D., Gemmi, M., Hanfland, Shiell, T. B., McCulloch, D. G., Bradby, J. E., Haberl, B., M., Kupenko, I., Lotti, P., Muller, H., & Zhang, L. (2017). Boehler, R., & McKenzie, D. R. (2016). Sci. Rep., 6, 37232. Am. Mineral., 102, 1763. Skorodumova, N. V., Belonoshko, A. B., Huang, L., Ahuja, R., Merlini, M., Hanfland, M., Salamat, A., Petitgirard, S., & & Johansson, B. (2005). Am. Mineral., 90, 1008. Muller, H. (2015). Am. Mineral., 100, 2001. Smith, E. M., Shirey, S. B., Nestola, F., Bullock, E. S., Wang, J. Mikhailenko, D. S., Korsakov, A. V., Zelenovskiy, P. S., & H., Richardson, S. H., & Wang, W. Y. (2016). Science, 354, Golovin, A. V. (2016). Am. Mineral., 101, 2155. 1403. Miller, E. D., Nesting, D. C., & Badding, J. V. (1997). Chem. Smith, R. F., et al. (2014). Nature, 511, 330. Mater., 9, 18. Solomatova, N. V., & Asimow, P. D. (2017). Am. Mineral., 102, 210. Mundy, C. J., Curioni, A., Goldman, N., Will Kuo, I. F., Reed, Spanu, L., Donadio, D., Hohl, D., Schwegler, E., & Galli, G. E. J., Fried, L. E., & Ianuzzi, M. (2008). J. Chem. Phys., 128, (2011). Proc. Natl. Acad. Sci. U.S.A., 108, 6843. 184701. Stanley, S., & Bloxham, J., (2004). Nature, 428, 151. Németh, P., Garvie, L. A. J., Aoki, T., Dubrovinskaia, N., Tschauner, O., et al. (2018). Science, 359, 1136. Dubrovinsky, L., & Buseck, P. R. (2014). Nat. Commun., 5, 5447. Turneaure, S. J., Sharma, S. M., Volz, T. J., Winey, J. M., & Nestola, F., et al. (2018). Nature, 555, 237. Gupta, Y. M. (2017). Sci. Adv., 3, eaao3561. Oganov, A. R., Glass, C. W., & Ono, S. (2006). Earth Planet. Utsumi, W., & Yagi, T. (1991). Science, 252, 1542. Sci. Lett., 241, 95. Vennari, C. E., & Williams, Q. (2018). Am. Mineral., 103, 171. Oganov, A. R., Hemley, R. J., Hazen, R. M., & Jones, A. P. Wang, X., Scandolo, S., & Car, R. (2005). Phys. Rev. Lett., 95, (2013). Rev. Mineral. Geochem., 75, 47. 185701. Oganov, A. R., Ono, S., Ma, Y. M., Glass, C. W., & Garcia, A. Wang, Y., Panzik, J. E., Kiefer, B., & Lee, K. K. M. (2012). Sci. (2008). Earth Planet. Sci. Lett., 273, 38. Rep., 2, 520. Ono, S., Kikegawa, T., & Ohishi, Y. (2007). Am. Mineral., 92, Wang, Z., et al. (2004). Proc. Natl. Acad. Sci. U.S.A., 101, 1246. 13699. Panero, W. R., & Kabbes, J. E. (2008). Geophys. Res. Lett., 35, Wen, X. D., Hoffmann, R., & Ashcroft, N. W. (2011). J. Am. L14307. Chem. Soc., 133, 9023. Pearson, D. G., et al. (2014). Nature, 507, 221. Wirth, R., Dobrzhinetskaya, L., Harte, B., Schreiber, A., & Pickard, C. J., & Needs, R. J. (2015). Phys. Rev. B, 91, 104101. Green, H. W. (2014). Earth Planet. Sci. Lett., 404, 365. Pineau, N. (2013). J. Chem. Phys. C, 117, 12778. Xie, Y.‐P., Zhang, X.‐J., & Liu, Z.‐P. (2017). J. Am. Chem. Soc., Powles, R. C., Marks, N. A., Lau, D. W. M., McCulloch, D. G., 139, 2545. & McKenzie, D. R. (2013). Carbon, 63, 416. Xu, J., Mao, H.‐K., & Hemley, R. J. (2002). J. Phys.: Cond. Prewitt, C. T., & Downs, R. T. (1998). Rev. Mineral. Geochem., Mater., 14, 11549. 37, 283 (1998). Yao, X., Xie, C., Dong, X., Oganov, A. R., & Zeng, Q. F. (2018). Pruzan, P., Chervin, J. C., Thiery, M. M., Itie, J. P., Besson, J. M., Phys. Rev. B, 98, 014108. Forgerit, J. P., & Revault, M. (1990). J. Chem. Phys., 92, 6910. Yin, M. T., & Cohen, M. L. (1983). Phys. Rev. Lett., 50, 2006. Rygg, J. R., et al. (2012). Rev. Sci. Instrum., 83, 113904. Zeng, Z., et al. (2017). Nat. Commun., 8, 322. Santoro, M., Ciabini, L., Bini, R., & Schettino, V. (2003). J. Raman Zhang, C., & Duan, Z. H. (2009). Geochim. Cosmochim. Acta, Spectrosc., 34, 557. 73, 2089.
You can also read