ORIGINAL RESEARCH article
Polymerized 4-Fold Coordinated Carbonate Melts in the Deep Mantle
- 1Muséum National d'Histoire Naturelle, UMR CNRS 7590, IRD, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, IMPMC, Sorbonne Université, Paris, France
- 2Institut des Sciences de la Terre de Paris, Sorbonne Université, CNRS, Paris, France
- 3SUPA, Centre for Science at Extreme Conditions and School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom
- 4DESY Photon Science, Hamburg, Germany
- 5European Synchrotron Radiation Facility, ESRF, Grenoble, France
Our understanding of the deep carbon cycle has witnessed amazing advances in the last decade, including the discovery of tetrahedrally coordinated high pressure (P) carbonate phases. However, little is known about the physical properties of their molten counterpart at moderate depths, while their properties at lower mantle conditions remain unexplored. Here, we report the structure and density of FeCO3 melts and glasses from 44 to 110 GPa by means of in situ x-ray synchrotron diffraction, and ex situ Raman and x-ray Raman spectroscopies. Carbon is fully transformed to 4-fold coordination, a bond change recoverable at ambient P. While low P melts react with silica, resulting in the formation of silico-carbonate glasses, high P melts are not contaminated but still quench as glasses. Carbonate melts are therefore polymerized, highly viscous and poorly reacting with silicates in the lower mantle, in stark opposition with their low P properties.
Although the lower mantle is mostly a reducing environment with the presence of reduced Fe (Frost et al., 2004; Smith et al., 2016), significant amount of subducted carbonates are estimated to be preserved (Litasov and Shatskiy, 2018). Transition to 4-fold carbon was first predicted for crystalline CaCO3 (Oganov et al., 2006; Arapan et al., 2007). This transition strongly depends on the carbonate composition, occurring for CaCO3 above 105 GPa (Lobanov et al., 2017), 80 GPa for MgCO3 (Oganov et al., 2008; Boulard et al., 2011), and 50 GPa for FeCO3 (Liu et al., 2015), while intermediate CaCO3-MgCO3-FeCO3 compositions form a single tetrahedral carbonate phase (Merlini et al., 2017) unlike silicates. This transition induces polymerization such as sheets or 3-membered rings for MgCO3 (Oganov et al., 2008), and chains for CaCO3 (Oganov et al., 2006). In contrast, our knowledge of carbonate melts structure at depth is scarce and limited to upper mantle pressures. The melting curves of CaCO3, Na2CO3, and FeCO3 have been measured over most of the upper mantle regime (Kang et al., 2015; Li et al., 2017), and viscosity measurements up to 6 GPa span several compositions [K2Ca(CO3)2 and K2Mg(CO3)2 by (Dobson et al., 1996), CaCO3 and natural dolomite by (Kono et al., 2014), Na2CO3 by (Stagno, 2018)]. Structural data instead have only been collected on molten CaCO3 below 10 GPa (Hudspeth et al., 2018) while theoretical investigations of the properties of carbonate melts cover a larger P-range but are also limited to the carbon 3-fold stability field (Vuilleumier et al., 2014; Zhang and Liu, 2015; Du et al., 2018; Desmaele et al., 2019). One main question is therefore how this 3 to 4-fold transition translates in the molten state, and what are the consequences on the physical and chemical properties of carbonate melts? Of particular interest is the mobility and reactivity of carbonate melts in the lower mantle, knowing that these properties underpin the key role played by carbonate melts in mantle geodynamics through lubrication of plate tectonics, cratonic roots (Foley, 2008) and ascending plumes (Litasov et al., 2013).
The role of Fe in the deep carbon cycle is emphasized by the predominance of Fe-rich ferropericlase in diamond inclusions from the lower mantle (Kaminsky, 2012). The lowest transition P from 3 to 4-fold C in FeCO3 amongst carbonates justifies its choice as the first composition to investigate. Not only this transition occurs at less challenging experimental conditions, but it might be driven by Fe high spin to low spin transition at 40.4 GPa (Weis et al., 2017), a consequence of which being the large enrichment in Fe of (Mg,Fe)-carbonates coexisting with bridgmanite to almost pure FeCO3 (Lobanov et al., 2015). Besides, high Fe concentration stabilizes (Ca,Mg,Fe)IVCO3 with respect to single cation 3-fold carbonates at mid mantle conditions (30–50 GPa) (Solomatova and Asimow, 2018). Formation of Fe-carbonates in the lower mantle might also result from carbonation of Fe-oxides [(Mg,Fe)O, FeOOH] with CO2 (Boulard et al., 2018, 2012). Last but not least, FeCO3 is a technical choice as it can be laser heated, which is required to reach lower mantle conditions without the need for additional laser coupler.
2. Materials and Methods
2.1. Materials and Chemical Analyses
The starting natural crystalline siderite sample (mineralogical collection at Sorbonne Université) was loaded in the sample chamber laser-drilled in a rhenium gasket as 20 μm-thick platelet between two equally thick platelets of compressed SiO2 powder. The SiO2 platelets act as thermal insulators and P-transmitting medium. Only one sample was used per P point (Figure 1) to avoid repeated laser-heatings, and preserve the chemical integrity of the sample. Six samples could be recovered after the experiments, embedded in epoxy and polished for analysis. Samples 8, 9, and 15 were carbon-coated for SEM imaging (Figure 2), samples 8 and 15 were then repolished and gold-coated along with samples 13, 14, and 20 for electron microprobe analysis using a CAMECA SX-FIVE analyzer (EMPA) at the Camparis center of Sorbonne Université (Table 1), using the following operating conditions: 15 keV, 10 nA. We used a defocused beam size of 10 μm to get an average composition at the laser-heated spot.
Figure 1. Microphotograph of the sample after laser heating at 110 GPa. Single shot laser heating resulted in the formation of a quasi-spherical pure carbonate glass that was removed from the gasket for EPMA and/or SEM analyses.
Figure 2. SEM images of recovered samples. Low P sample 8 (a) shows pervasive contamination of carbonate sample with SiO2 P-transmitting medium. High P samples 9 (b) and 15 (c) show that chemical integrity of carbonate melt (homogeneous light gray zone) was preserved.
Table 1. Run conditions, quenched products and their chemical composition in wt% obtained from EMPA. One standard deviations are given in parentheses. Starting natural siderite sample also contained <0.1 wt% CaO and MnO.
2.2. P-T Conditions
We used diamond-anvil cells and a double-sided infra-red laser focussed down to 20 μm to generate high T and P. For each P point, targeted power was increased in 2 W increments from 20 to 50 W of power on each laser depending on P until complete melting of the sample. Melting was identified by disappearance of diffraction peaks apart from SiO2 peaks, and by the appearance of diffuse scattering. As we used the off-axis heating system to avoid using carbon mirrors that would add to the x-ray background signal and compromise processing of the scattered signal, T could not be measured by pyrometric techniques. FeCO3 melting curve has only been measured up to 20 GPa (Kang et al., 2015), where it reaches 1,865 K. The stishovite to CaCl2 SiO2 transition has been investigated up to 90 GPa (Fischer et al., 2018), this constrains T to a maximum of 2,300 K at 79 GPa and 2,500 K at 83 GPa as CaCl2 is the observed SiO2 structure for the three highest P runs, while stishovite is observed below. We therefore consider that x-ray diffraction patterns were collected on molten FeCO3 within the 2,000 K–2,500 K interval except for the highest P point that is only constrained to below 3,500 K from extrapolation of the stishovite-CaCl2 Clapeyron slope (Fischer et al., 2018). P is measured at room T using fluorescence of a ruby sphere added in the sample chamber (Mao et al., 1986) and SiO2 equations of state (Andrault et al., 1998; Nishihara et al., 2005) for quenched samples, and using only SiO2 equations of state for molten samples with error bars on P including the effect of a 2,000 K–2,500 K T-range, and up to 3,500 K for the 110 GPa data point.
2.3. X-ray Diffraction Methods
We collected in situ high P-T x-ray diffraction data in laser-heated diamond anvil cells at the extreme conditions beamline P02.2 at the PETRAIII synchrotron. We used symmetric diamond-anvil cells equipped with 70° opening Boehler-Almax seats in order to access a wider q-range up to 10 Å−1, and reduce the diamond Compton contribution as Boehler-Almax anvils are only 1.5 mm thick. The x-ray monochromatic beam (42.7 keV) was focussed down to a size of 4 × 6 μm2, allowing high spatial resolution in direct space. To limit iron migration away from the laser heating spot due to Soret effect, the laser shutters were opened only once the targeted power was reached, and held open for 10 s during which 10 x-ray diffraction patterns of 1 s acquisition time were recorded on a Perkin-Elmer 2-D detector. 2-D patterns were integrated using the Fit2D software (Hammersley et al., 1996). In order to isolate the scattered intensity from the molten FeCO3 only, each sample was removed from the gasket, and the gasket put back in place to collect x-ray data on the empty cell. Obtained patterns were then scaled vertically to match the baseline of x-ray patterns collected on the starting crystalline sample under P (Sanloup and de Grouchy, 2018). This last step ensures that any P effect on the background is corrected for. Amongst eight successful runs (Table 1) for which full melting was observed, intensity from molten FeCO3 could only be processed for the highest P run for which the sample vs. SiO2 platelets thickness ratio was slightly higher, the scattered intensity being too weak for the lower P points. All glass patterns could be processed. The x-ray diffracted intensity data are converted into the structure factor, S(q) (Figures 3A, 4), using the Ashcroft-Langreth formalism. The radial distribution function g(r) (Figure 3B), that describes ion-ion contributions in real space, is obtained by Fourier transforming of S(q),
where , is the Avogadro number, M the mean atomic molar mass, and ρ the density.
Figure 3. Structure of non-crystalline FeCO3 at high pressures. (A) Structure factor, S(q), for all quenched glasses (black) and the highest P melt (red). (B) Corresponding radial distribution functions, g(r). Dashed lines are fits to the C-O contribution at 1.2-1.3 Å where there is no overlap with farther contributions.
Figure 4. Structure factor, S(q), for low P reacted FeCO3+SiO2 glass (black), compared to SiO2 glass at 20 GPa (Sato and Funamori, 2008) (brown) and FeCO3 glass at 44 GPa (red). Low P sample 8 (15 GPa) shows intermediate structure between SiO2 glass and high P FeCO3 glasses.
2.4. Density Measurements
The method to derive density from x-ray diffraction data on melts compressed in diamond-anvil cell experiments (Eggert et al., 2002; Sanloup et al., 2013) consists in minimizing the oscillations in g(r) where there should not be any signal, i.e., below the minimum interatomic distance (r < 0.95 Å here). This method requires that the background, essentially the Compton signal from the diamond anvils that dominates the total diffracted intensity, is perfectly subtracted.
As the C-O contribution is distinct on g(r) of quenched glasses up to 83 GPa, we also ran consistency checks by fixing the C-O coordination number to 4 as indicated by x-ray Raman spectra (cf Results section), and simulating the C-O contribution using the obtained density values against a gaussian with the following equation:
where cC and cO are the atomic proportions of carbon and oxygen, KC and KO are defined as the average effective atomic number over the experimental q-range (Eggert et al., 2002) and calculated using form factors from Hajdu (1972).
with Kp, the effective atomic number (Eggert et al., 2002), Ztot the total atomic number of the compositional unit (e.g., FeCO3), CN the C-O coordination number fixed to 4, d the C-O inter-atomic distance, and σ a parameter depending on structural disorder, where k is an adjustable parameter (Hosemann and Bagchi, 1962) with a value of 0.11 here. The C-O contribution to g(r) thus calculated adequately fits the experimental ion-ion contribution (dashed lines on Figure 2b), hence comforting the obtained density values.
2.5. X-ray Raman and Raman Methods
Raman and x-ray Raman spectra were collected at ambient conditions on glassy FeCO3 recovered from x-ray diffraction experiments and from additional laser-heated diamond anvil cell synthesis, respectively.
X-ray Raman data were collected at an incident energy of 9.7 keV at the C K-edge on beamline ID20 of the European Synchrotron Radiation Facility (ESRF), beamsize was 15 × 15 μm2. The large-solid-angle x-ray scattering spectrometer (Huotari et al., 2017) was set up with 24 Si(660) analyzer crystals for an average momentum transfer of 7.3 ± 0.2 Å−1 and an overall energy resolution of 0.7 eV. All experimental data were analyzed using the XRStools software package (Sahle et al., 2015). The integrated intensity of each spectrum was normalized over a 35 eV energy range. Glassy FeCO3 spheres had been previously synthesized at 59 GPa using the same P02.2 laser heating system in PetraIII as for x-ray diffraction experiments. LiF was used instead of SiO2 as a P-transmitting medium to avoid any contamination of the x-ray Raman signal by oxygen from SiO2 as measurements at the O K-edge were initially planned but signal was to weak for data to be processed. Despite its higher melting curve than siderite (Boehler et al., 1997), LiF salt could not be used for the x-ray diffraction experiments due to its continuous powder diffracted signal that prevents a qualitative analysis of the diffuse scattering signal from molten FeCO3.
Raman spectra were collected on glassy FeCO3 recovered from x-ray diffraction experiments using 633 nm wavelength in order to preserve the samples, using more energetic lower wavelengths resulted in dissociation of the sample and detection of hematite signal.
All samples are systematically quenched as a glass. Chemical integrity of FeCO3 molten spheres is observed for runs conducted above 40 GPa, apart from a marginal fraction at the glass-SiO2 interface in one sample showing enrichment of the P-transmitting medium in Fe and C. Instead, the lowest P samples, i.e., 11 and 15 GPa, have reacted with the SiO2 P-transmitting medium. This is shown by SEM imaging (Figure 2) and EMPA analysis on sample 8 (Table 1). High P carbonate melts are thus much less reactive than low P melts. This might not contradict the observed reactivity of high P crystalline MgCO3 with SiO2 (Seto et al., 2008; Maeda et al., 2017) due to the much longer heating durations (20–240 min against 10 s heating duration in this work); alternatively, Fe stabilizing effect on high P carbonates could be at stake. We observe no disproportionation of Fe as was reported in the crystalline state in some studies (Boulard et al., 2011; Cerantola et al., 2017) but not in others (Liu et al., 2015). This might be due to different P-T paths followed, i.e., flash heating here instead of continuous T increase (Boulard et al., 2011; Cerantola et al., 2017).
A striking characteristic of glassy FeCO3 is its strong first sharp diffraction peak (FSDP) that persists in the structure factor up to the highest P investigated (Figure 3A), indicative of a strong medium-range order. This is in stark contrast to silicate glasses that lose their medium-range order with increased P (Sato and Funamori, 2008), but consistent with ab initio calculations on carbon-bearing silicate melts reporting P-induced polymerization of carbonate species into dimers and with the silicate network (Ghosh et al., 2017; Solomatova and Asimow, 2019). A second noticeable feature is the decrease of the contribution at 4 Å−1 attributed in molten carbonates to the O-O bond (Wilding et al., 2016). On radial distribution functions, g(r) (Figure 3B), the C-O contribution is clearly visible at 1.2-1.3 Å with none or little overlap with the second contribution (Fe-O and O-O) at ~2 Å in the glass, and with some overlap in the melt. No significant structural changes are observed between molten and quenched glassy state at 110 GPa, apart from a generally lower intensity in the melt due to the high T and consequent higher degree of disorder. For g(r), this weaker intensity translates into broader C-O and Fe-O/O-O contributions in the molten state. For glasses quenched at 11 GPa and 15 GPa, the x-ray structure factor, S(q), is intermediate between that of pure SiO2 glass (Sato and Funamori, 2008) and high-P FeCO3 glasses (Figure 4). SEM image of sample 8 (15 GPa, Figure 2) shows heterogeneities in the quenched glass, which indicates that the x-ray structure factor likely averages at least two types of glass structure and therefore data cannot be interpreted quantitatively.
The x-ray Raman C K-edge spectrum of quenched FeCO3 glass shows no presence of sp2 3-fold carbon characterized by an intense π* peak at 290 eV (Figure 5, π* peak). Only the σ* peak of tetrahedrally coordinated carbon (Shieh et al., 2013) is visible (Figure 5, σ* peak). The totally missing π* peak is indicative of a fully sp3 state of carbon atoms in the siderite glass. P-induced coordination changes of major cations in silicate melts (e.g., Si, Al) were first reported from the study of glasses quenched from high P (Yarger et al., 1995; Meade, Hemley and Mao, 1992), and later confirmed by in situ studies in the molten phase (Sanloup et al., 2013; Drewitt, 2015). However, the opposite, i.e., coordination change occurring only in the quenched glass, not in the high P melt, have not been reported nor been theoretically predicted. The 3 to 4-fold transition therefore occurs in molten Fe-carbonates at P less or equal to 51 GPa. This transition is preserved upon quenching to the glassy state, and is recoverable at ambient conditions, opening the way to the synthesis of a new class of glassy materials. Two broad bands are observed in the Raman spectra (Figure 6), very different from those of the only two carbonate systems that quench as glasses at room P, MgCO3-K2CO3 and La(OH)3-Ca(OH)2-CaCO3-CaF2BaSO4 (Sharma and Simons, 1979), that are essential dominated by the strong CO stretching mode at ~1,080 cm−1. Instead, present Raman spectra are reminiscent of those reported for calcium silicate glasses (Figure 6) (Mysen et al., 1982) albeit at higher Raman shift values for the broadest band (1,200–1,600 cm−1 for glassy FeCO3 vs. 850–1,100 cm−1 for calcium silicate glasses).
Figure 5. X-ray Raman spectra collected at the carbon K-edge on crystalline siderite and high P-quenched FeCO3 glasses at ambient conditions. The disappearance of the π* feature, which is solely related to the three-fold coordinated carbon, is a spectroscopic evidence of a full four-fold coordination state in the glassy structure of FeCO3.
Figure 6. Raman spectra collected on high P-quenched FeCO3 glasses (runs 9 and 14) at ambient conditions.
Density values are reported in Figure 7 along with predictions for lower P melt properties (Kang et al., 2015), P-evolution of crystalline siderite, and with the Earth's seismological PREM model (Dziewonski and Anderson, 1981). Density profile below 40 GPa is calculated using KT, 0 value of 80.23 GPa (Kang et al., 2015), consistent with that reported for molten calcite (Hudspeth et al., 2018), and density at room P of 2,500 kg·m−3 by assuming a similar density jump upon melting as for other carbonates for which room P density is known. Comparison with PREM model shows that Fe-carbonate melts are buoyant at all depths. Density contrast between the high P polymerized melt or glass and extrapolated equation of state for low P melt is 15%, i.e., similar to volume collapse reported upon transition from crystalline high spin siderite I to low spin siderite II (Liu et al., 2015). The volume collapse is smoothed out over a ~ 30 GPa range in the molten state with, as a direct consequence, a steepening of the melting curve from 55 GPa on (Cerantola et al., 2017).
Figure 7. Density evolution of glassy, molten, and cristalline siderite with pressure. Molten low P siderite (plain curve), high P data on glass (black points) and the highest P melt (red point), crystalline equation of state (dashed curve) includes the transition from high spin siderite I to low spin siderite II at 50 GPa (Liu et al., 2015).
The 3 to 4-fold transition occurs in molten Fe-carbonates at P less or equal to 51 GPa, compared to 50 GPa for crystalline FeCO3 (Liu et al., 2015), 80 GPa for MgCO3 (Oganov et al., 2008; Boulard et al., 2011), and 130 GPa for CaCO3 (Oganov et al., 2006; Arapan et al., 2007). A consequence of the effect of Fe on the 3-fold C to 4-fold C transition P is that crystalline Fe-poor (Ca,Mg,Fe)IIICO3 and Fe-enriched (Ca,Mg,Fe)IVCO3 melts could co-exist at depth. In the case of Si isotopes, fractionation between VISi bridgmanite and IVSi olivine structures is theoretically estimated to ~−1 ‰28Si at 2000 K (Huang et al., 2014). If this effect can be scaled to C simply using mass difference considerations, then a few ‰13C fractionation is expected, and could potentially explain isotopic differences between calcite inclusions from super-deep diamonds (Kaminsky et al., 2016). This effect might be sufficient to confer a mantle-like signature to deep diamonds grown from slab-derived carbonate melts while co-existing tetrahedral crystalline carbonate are expected to get lighter.
Density of non-crystalline FeCO3 remains considerably lower than that of its crystalline counter parts, even at the highest investigated P, by 15%. The situation is thus very different from that of molten and crystalline silicates which density converge at deep mantle conditions (Petitgirard et al., 2015; Sanloup, 2016), and such difference could be attributed to the very strong medium-range order preserved in tetrahedral high P carbonate melts while it is mostly collapsed by 5 GPa in silicate melts. That high P FeCO3 melts quench as glasses contrasts with the behavior observed at lower P, and suggests an important increase of carbonate melt viscosity consistent with the observation of a very strong medium-range order. It is also opposite to the behavior of molten basalt that systematically quenches as crystalline phases above 11 GPa (Sanloup et al., 2013) and as a glass below. The strongly reduced chemical reactivity of high P FeCO3 melts with silica along with their glass-forming ability suggest that unlike at lower P, tetrahedral carbonate melts are not pervasive, which could contribute to the longevity of carbonates in the deep mantle where allowed by oxidizing conditions or slow reduction kinetics (Litasov and Shatskiy, 2018).
JH and CS: devised the project; CS: processed x-ray diffraction data and wrote the paper with input from GL and LC; JH, CS, BC, VA, and ZK: participated in x-ray diffraction data acquisition; JH, GL, LC, CS, and CC: participated in x-ray Raman data acquisition; JH: collected Raman data.
The research leading to these results has received funding from the European Community's Seventh Framework Programme (FP7/2007-2013) under grant agreements no. 312284 and 259649 (European Research Council starting grant to CS). Portions of this research were carried out at the light source PETRAIII at DESY, a member of the Helmholtz Association (HGF). The laser heating system on beamline P02.2 is funded by the German BMBF (project number 05K10RFA).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We acknowledge E. Boulard for providing the starting siderite sample, K. Glazyrin for his help with ex situ diamond-anvil cell laser-heating synthesis in PETRAIII, L. Rémusat at Museum National d'Histoire Naturelle (Paris, France) for gold coating of the recovered samples, the ESRF (Grenoble, France) and PETRAIII (Hamburg, Germany) for provision of synchrotron radiation facilities.
Boulard, E., Gloter, A., Corgne, A., Antonangeli, D., Auzende, A.-L., Perrillat, J.-P., et al. (2011). New host for carbon in the deep Earth. Proc. Natl Acad. Sci. U.S.A. 108, 5184–5187. doi: 10.1073/pnas.1016934108
Boulard, E., Guyot, F., and Fiquet, G. (2012). The influence on Fe content on Raman spectra and unit cell parameters of magnesite-siderite solid solutions. Phys. Chem. Miner. 39, 239–246. doi: 10.1007/s00269-011-0479-3
Boulard, E., Guyot, F., Menguy, N., Corgne, A., Auzende, A.-L., Perrillat, J.-P., et al. (2018). CO2-induced destabilization of pyrite-structured FeO2Hx in the lower mantle. Natl. Sci. Rev. 5, 870–877. doi: 10.1093/nsr/nwy032
Cerantola, V., Bykova, E., Kupenko, I., Merlini, M., Ismailova, L., McCammon, C., et al. (2017). Stability of iron-bearing carbonates in the deep Earth's interior. Nat. Commun. 8:15960. doi: 10.1038/ncomms15960
Desmaele, E., Sator, N., and Vuilleumier, R. Guillot, B. (2019). Atomistic simulations of molten carbonates: thermodynamic and transport properties of the Li2CO3-Na2CO3-K2CO3 system. J. Chem. Phys. 150:094504. doi: 10.1063/1.5082731
Dobson, D. P., Jones, A. P., Rabe, R., Sekine, T., Kurita, K., Taniguchi, T., et al. (1996). In-situ measurement of viscosity and density of carbonate melts at high pressure. Earth Planet. Sci. Lett. 143, 207–215. doi: 10.1016/0012-821X(96)00139-2
Drewitt, J. W. E., Jahn, S., Sanloup, C., de Grouchy, C., and Garbarino, G. Hennet L. (2015). Development of chemical and topological structure in aluminosilicate liquids and glasses at high pressure. J. Phys.: Cond. Matt. 27:105103. doi: 10.1088/0953-8984/27/10/105103
Eggert, J. H., Weck, G., Loubeyre, P., and Mezouar, M. (2002). Quantitative structure factor and density measurements of high-pressure in diamond anvil cells by x-ray diffraction: Argon and water. Phys. Rev. B 65:174105. doi: 10.1103/PhysRevB.65.174105
Fischer, R. A., Campbell, A. J., Chidester, B. A., Reaman, D. M., Thompson, E. C., Pigott, J. S., et al. (2018). Equations of state and phase boundary for stishovite and CaCl2-type SiO2. Am. Mineral. 103, 792–802. doi: 10.2138/am-2018-6267
Frost, D., Liebske, C., Langenhorst, F., McCammon, C., Tronnes, R., and Rubie, D. (2004). Experimental evidence for the existence of iron-rich metal in the Earth's lower mantle. Nature 428, 409–412. doi: 10.1038/nature02413
Hammersley, A. P., Svensson, S. O., Hanfland, M., Fitch, A. N., and Hausermann, D. (1996). Two-dimensional detector software: From real detector to idealised image or two-theta scan. High Press. Res. 14, 235–248. doi: 10.1080/08957959608201408
Huang, F., Wu, Z., Huang, S., and Wu, F. (2014). First-principles calculations of equilibrium silicon isotope fractionation among mantle minerals. Geochim. Cosmochim. Acta 140, 509–520. doi: 10.1016/j.gca.2014.05.035
Huotari, S., Sahle, C. J., Henriquet, C., Al-Zein, A., Martel, K., Simonelli, L., et al. (2017). A large-solid-angle X-ray Raman scattering spectrometer at ID20 of the European Synchrotron Radiation Facility. J. Sync. Rad. 24, 521–530. doi: 10.1107/S1600577516020579
Kaminsky, F., Matzel, J., Jacobsen, B., Hutcheon, I., and Wirth, R. (2016). Isotopic fractionation of oxygen and carbon in decomposed lower-mantle inclusions in diamond. Mineral. Petrol. 110, 379–385. doi: 10.1007/s00710-015-0401-7
Kang, N., Schmidt, M. W., Poli, S., Franzolin, E., and Connolly, A. D. (2015). Melting of siderite to 20 GPa and thermodynamic properties of FeCO3-melt. Chem. Geol. 400, 34–43. doi: 10.1016/j.chemgeo.2015.02.005
Li, Z., Li, J., Lange, R., Liu, J., and Militzer, B. (2017). Determination of calcium carbonate and sodium carbonate melting curves up to Earth's transition zone pressures with implications for the deep carbon cycle. Earth Planet. Sci. Lett. 457, 395–402. doi: 10.1016/j.epsl.2016.10.027
Lobanov, S. S., Goncharov, A. F., and Litasov, K. D. (2015). Optical properties of siderite (FeCO3) across the spin transition: Crossover to iron-rich carbonates in the lower mantle. Am. Mineral. 100, 1059–1064. doi: 10.2138/am-2015-5053
Lobanov, S. S., Stevanovic, V., Gavryushkin, P. N., Litasov, K. D., Greenberg, E., Prakapenka, V. B., and Oganov, A. R. Goncharov, Al. F. (2017). Raman spectroscopy and x-ray diffraction of sp(3) CaCO3 at lower mantle pressures. Phys. Rev. B 96:104101. doi: 10.1103/PhysRevB.96.104101
Maeda, F., Ohtani, E., Kamada, S., Sakamaki, T., Hirao, N., and Ohishi, Y. (2017). Diamond formation in the deep lower mantle: a high-pressure reaction of MgCO3 and SiO2. Sci. Rep. 7:40602. doi: 10.1038/srep40602
Merlini, M., Cerantola, V., Gatta, G. D., Gemmi, M., Hanfland, M., Kupenko, I., et al. (2017). Dolomite-IV: candidate structure for a carbonate in the Earth's lower mantle. Am. Mineral. 102, 1763–1766. doi: 10.2138/am-2017-6161
Nishihara, Y., Nakayama, K., Takahashi, E., Iguchi, T., and i. Funakoshi, K. (2005). P-V-T equation of state of stishovite to the mantle transition zone conditions. Phys. Chem. Min. 31, 660–670. doi: 10.1007/s00269-004-0426-7
Oganov, A. R., Ono, S., Ma, Y., Glass, C. W., and Garcia, A. (2008). Novel high-pressure structures of MgCO3, CaCO3 and CO2 and their role in Earth's lower mantle. Earth Planet. Sci. Lett. 273, 38–47. doi: 10.1016/j.epsl.2008.06.005
Petitgirard, S., Malfait, W. J., Sinmyo, R., Kupenko, I., Hennet, L., Harries, D., et al. (2015). Fate of MgSiO3 melts at core-mantle boundary conditions. P. Natl. Acad. Sci. U.S.A. 112, 14186–14190. doi: 10.1073/pnas.1512386112
Sahle, C. J., Mirone, A., Niskanen, J., Inkinen, J., Krisch, M., and Huotari, S. (2015). Planning, performing and analyzing X-ray Raman scattering experiments. J. Sync. Rad. 22, 400–409. doi: 10.1107/S1600577514027581
Sanloup, C., Drewitt, J. W. E., Konôpková, Z., Dalladay-Simpson, P., Morton, D. M., Rai, N., et al. (2013). Structural change in molten basalt at deep mantle conditions. Nature 503, 104–107. doi: 10.1038/nature12668
Seto, Y., Hamane, D., Nagai, T., and Fujino, K. (2008). Fate of carbonates within oceanic plates subducted to the lower mantle, and a possible mechanism of diamond formation. Phys. Chem. Miner. 35, 223–229. doi: 10.1007/s00269-008-0215-9
Shieh, S. R., Jarrige, I., Wu, M., Hiraoka, N., Tse, J. S., Mi, Z., et al. (2013). Electronic structure of carbon dioxide under pressure and insights into the molecular-to-nonmolecular transition. Proc. Natl. Acad. Sci. U. S. A. 110, 18402–18406. doi: 10.1073/pnas.1305116110
Smith, E. M., Shirey, S. B., Nestola, F., Bullock, E. S., Wang, J., Richardson, S. H., et al. (2016). Large gem diamonds from metallic liquid in earth's deep mantle. Science 354, 1403–1405. doi: 10.1126/science.aal1303
Solomatova, N. V. and Asimow, P. D. (2018). First-principles calculations of high-pressure iron-bearing monoclinic dolomite and single-cation carbonates with internally consistent Hubbard U. Phys. Chem. Miner. 45, 293–302. doi: 10.1007/s00269-017-0918-x
Stagno, V., Stopponi, V., Kono, Y., and Manning, C. E. Tetsuo, I. (2018). Experimental determination of the viscosity of Na2CO3 melt between 1.7 and 4.6 GPa at 1200-1700 degrees C: implications for the rheology of carbonatite magmas in the Earth's upper mantle. Chem. Geol. 501,19–25. doi: 10.1016/j.chemgeo.2018.09.036
Vuilleumier, R., Seitsonen, A., Sator, N., and Guillot, B. (2014). Structure, equation of state and transport properties of molten calcium carbonate (CaCO3) by atomistic simulations. Geochim. Cosmchim. Acta 141, 547–566. doi: 10.1016/j.gca.2014.06.037
Weis, C., Sternemann, C., Cerantola, V., Sahle, C. J., Spiekermann, G., Harder, M., et al. (2017). Pressure driven spin transition in siderite and magnesiosiderite single crystals. Sci. Rep. 7:16526. doi: 10.1038/s41598-017-16733-3
Wilding, M. C., Wilson, M., Alderman, O. L. G., Benmore, C., Weber, J. K. R., Parise, J. B., et al. (2016). Low-dimensional network formation in molten sodium carbonate. Sci. Reports 6:24415. doi: 10.1038/srep24415
Yarger, J. L., Smith, K. H., Nieman, R. A., Diefenbacher, J., Wolf, G. H., Poe, B. T., and McMillan, P. F. (1995). Al coordination changes in high-pressure aluminosilicate liquids. Science 270, 1964–1967.
Keywords: carbonate melts, high pressure, x-ray diffraction, Raman, x-ray Raman, polymerization, deep mantle
Citation: Sanloup C, Hudspeth JM, Afonina V, Cochain B, Konôpková Z, Lelong G, Cormier L and Cavallari C (2019) Polymerized 4-Fold Coordinated Carbonate Melts in the Deep Mantle. Front. Earth Sci. 7:72. doi: 10.3389/feart.2019.00072
Received: 15 January 2019; Accepted: 19 March 2019;
Published: 10 April 2019.
Edited by:Dan J. Bower, University of Bern, Switzerland
Reviewed by:Zhicheng Jing, Southern University of Science and Technology, China
Alexander F. Goncharov, Carnegie Institution for Science (CIS), United States
Copyright © 2019 Sanloup, Hudspeth, Afonina, Cochain, Konôpková, Lelong, Cormier and Cavallari. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chrystèle Sanloup, firstname.lastname@example.org