Access the full text.
Sign up today, get DeepDyve free for 14 days.
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
We present a theoretical investigation, based on ab initio calculations and the quasi-harmonic ap- proximation, on the stability properties of magnesium (MgCO ) and calcium (CaCO ) carbonates 3 3 at high temperatures and pressures. The results indicate that those carbonates should be stable in Earth’s lower mantle, instead of dissociating into other minerals, in chemical environments with excesses of SiO , MgO, or MgSiO . Therefore, considering the lower mantle chemical composition, 2 3 ﬁlled with major minerals such as MgSiO and MgO, calcium and magnesium carbonates are the primary candidates as carbon hosts in that region. For the thermodynamic conditions of the mantle, the results also indicate that carbon should be primarily hosted on MgCO . Finally, the results indicate that carbon, in the form of free CO , is unlikely in the lower mantle. keywords: carbonates, high-pressure minerals, lower mantle arXiv:1804.07252v1 [cond-mat.mtrl-sci] 19 Apr 2018 I. INTRODUCTION Carbon is a unique chemical element, mainly due to its rich bonding nature, which provides a wide range of stable and metastable structures in several hybridizations and conﬁgurations. Particularly, the role of carbon on Earth’s natural phenomena has been extensively studied over the last few decades. The carbon cycle aﬀects atmosphere, oceans and other shallow crustal phenomena, directly inﬂuencing climates and ecosystems, and consequently life on Earth. While there is considerable accumulated knowledge on the carbon cycle near Earth’s surface, there is still scarce information on the processes associated with its deep layers . In order to build consistent models on the Earth’s carbon cycle, it is important to establish a proper understanding on its chemical composition. The knowledge on the solar system composition, based on information from carbonaceous chondritic meteorites that have hit the Earth, allows to estimate the expected amount of carbon that should be present on Earth . However, the current known carbon reservoirs on its shallow layers is about two orders of magnitude smaller than the expected value, suggesting that this missing carbon should be stored on its deep layers, which has been labelled as deep carbon reservoirs . Estimates indicate that the Earth’s deep interior may contain as much as 90% of all available carbon. Such conclusions have been supported by a number of indirect evidences, such as 2 3–5 the presence of CO in magmas and mantle mineral inclusions in natural diamonds . There are several questions on the properties of deep carbon that remain open, such as determining the amount and distribution of those carbon reservoirs within major mantle minerals, and understanding the complete carbon cycle, associated with exchange of carbon between Earth’s surface and its deep interior. For example, current estimates suggest that the carbon ﬂux into Earth’s mantle through subduction substantially exceeds the carbon ﬂux emitted by volcanoes, opening the question on the minerals that host such carbon in the mantle. Those issues could only be addressed with a deep understanding on the behavior of carbon-related minerals under high pressures and temperatures, particularly on the ther- modynamic conditions of the mantle . Although some carbon may be stored in the core, the mantle is thought to be its largest reservoir. It is reasonable to assume that a certain amount of carbon could be dissolved in silicates, such as MgSiO , which is by far the most 2 abundant mineral in the lower mantle. Therefore, even if this mineral had a low to moderate solubility for carbon, it could still be the largest carbon reservoir in that region. However, it has been experimentally shown that the carbon solubility in silicates is very small, in- dicating that most of the Earth’s carbon must be stored in other minerals . Furthermore, silicates are large diﬀuse reservoirs, which contrasts with the concentrated phases of the deep 8 9 carbon identiﬁed as magma carbonates and diamonds , as well as carbon species emitted by volcanoes . In the lower mantle, it is still not known if carbon exists in a reduced state, such as in the form of diamond, or in an oxidized one, such as in the form of carbonates. To further explore the potential carriers of carbon, it is important to explore the physical properties of those minerals at lower mantle thermodynamic conditions, particularly in terms of carbon-related 11–13 14,15 stable phases. Several theoretical and experimental investigations have explored the high-pressure stability of major carbonates, such as MgCO , CaCO , and MgCa(CO ) . 3 3 3 2 16,17 However, there is still scarce information on those properties at high temperatures . Particularly for the lower mantle, such conditions would mean pressures up to 140 GPa and temperatures up to 3000 K. This investigation explores the stability of carbonates, particularly MgCO and CaCO , 3 3 at lower mantle conditions. Here, we did not explore the MgCa(CO ) mineral, since it 3 2 is well established that it dissociates into MgCO and CaCO for pressures higher than a 3 3 few GPa . The results were obtained by ﬁrst principles total energy calculations, combined with the quasi-harmonic approximation, in order to obtain the respective Gibbs free energies of those minerals in diﬀerent crystalline phases, at high temperatures and pressures. We explored the stability of those carbonates, taking into account a number of crystalline phases that have been identiﬁed by recent theoretical and experimental investigations. Our results indicate that, at high temperatures and high pressures, MgCO and CaCO should be 3 3 stable against dissociation into other minerals at several mantle conditions. Considering the chemical composition of the lower mantle, with major concentrations of MgSiO and MgO, calcium and magnesium carbonates should be the primary candidates for carbon hosts in the lower mantle. However, in the thermodynamic conditions of the mantle, along its geotherm, magnesium carbonate is more favorable than calcium carbonate. Additionally, the results suggested that carbon in the form of isolated CO is unlikely in the lower mantle. 3 II. METHODS A. Ab initio calculations The ﬁrst principles calculations were performed using the Quantum ESPRESSO compu- tational package . The electronic interactions were described within the density functional theory, considering the exchange-correlation (XC) potential based on the local density ap- proximation (LDA) functional . This functional has been widely used to obtain both static and dynamic properties of minerals of the Earth’s mantle, although many authors have investigated those properties with the generalized gradient approximation (GGA) . It is well established in the literature that static calculations with the LDA underestimate the mineral lattice parameters (generally by about 1 to 2 %) and overestimate the respective phase transition pressures and elastic constants when compared to experimental values at ﬁnite temperatures, while calculations with the GGA provide the opposite eﬀects . This investigation used only the LDA functional, since by incorporating thermal eﬀects later on, the lattice parameters increase (and elastic constants decrease), going with the appropriate trend toward the respective experimental values. The electronic wave functions were expanded using the projected augmented wave (PAW) method , with a plane-wave cutoﬀ of 1200 eV. The valence electronic conﬁgurations were 2 6 2 0 0 2 2 2 6 2 0 described with (3s 3p 4s 3d 4p ) for calcium, (2s 2p ) for carbon, (2s 2p 3s 3p ) 2 2 2 4 24 for magnesium, (3s 3p ) for silicon, and (2s 2p ) for oxygen . The Brillouin zones for electronic states of crystalline phases were sampled by a 4× 4× 4 k-mesh for the carbonates and silicates, a 8× 8× 8 k-mesh for alkaline earth oxides (MgO and CaO), and a 6× 6× 6 k- mesh for SiO and CO , in order to provide an approximately equivalent density of k-points 2 2 for all materials. Strict convergence criteria were taken into account for the simulations, with the atomic positions being considered converged when all forces acting on atoms were smaller than −1 0.01 eV A . For each pressure, from 0 to 150 GPa, the structures were optimized using the damped variable cell shape molecular dynamics method . Then, a second order Birch- Murnaghan equation of state was used to ﬁt the compression data. 4 B. Thermodynamic properties The thermodynamic properties were investigated by computing the vibrational modes (phonons) of the crystals using the Density Functional Perturbation Theory , with a 12 × 12 × 12 q-mesh to calculate the vibrational density of states (VDOS) , and the quasi- harmonic approximation (QHA) to compute the respective Gibbs free energies. The results for ﬁnite temperatures were reported within the validity range of the QHA . It is well established in the literature that the QHA results for the thermal expansion, α(P,T), start diverging beyond a certain temperature, in opposition to available experi- mental results. Such divergence results from the fact that this methodology disregards 28,29 anharmonic eﬀects . However, this quantity is also clearly quite sensitive to the choice of XC functional used , and could still be used at high temperatures, providing results consistent with experimental data in the region in which the thermal expansion coeﬃcient does not diverge. Therefore, at high temperatures, the validity region of the QHA has been 28,30 2 2 established as [∂ α(P,T)/∂T ] ≤ 0. Figure 1 shows the thermal expansion coeﬃcient at several pressures for CaCO in arago- nite and post-aragonite phases, which represents a stringent test for this methodology. The ﬁgure shows the respective validity regions of the QHA for those phases, indicating that at very high pressures, the QHA is still valid at temperatures up to 3000 K. The ﬁgure also shows that within the validity limit of the QHA, the thermal expansion coeﬃcient is in good agreement with available experimental data . We explored the validity of the QHA at the lowest pressure phases for all minerals con- sidered in this investigation. The region of low pressures and high temperatures represents a stringent test for the QHA methodology, since it is when the thermal expansion coeﬃ- cients present major divergences when compared to experimental data. The QHA provided appropriate thermal expansion coeﬃcients under the thermodynamic conditions of interest in this investigation. Our results on thermal properties of all minerals investigated here are within the validity region of the QHA, and in good agreement with available theoretical and experimental 28,31–34 data . 5 C. Crystalline phases In order to explore the stability of carbon-related minerals, within the thermodynamic conditions of the lower mantle, we initially studied several crystalline phases of CaCO and MgCO carbonates and CaSiO and MgSiO silicates in a wide pressure range, by computing 3 3 3 their respective enthalpies and Gibbs free energies. For any mineral at a certain temperature and pressure, the respective stable phase is determined as the one with the lowest Gibbs free energy value, using the theoretical model described in the previous section. Additionally, a mineral follows a phase transition at a certain pressure when the free energy of the stable phase becomes higher than the one of a diﬀerent phase. Our results on stability and phase transition pressures of several carbon- ate and silicate minerals were in good agreement with results from other theoretical and 11–13,15 experimental investigations . We considered a number of crystalline phases for CaCO , MgCO , MgSiO , CaSiO , 3 3 3 3 MgO, CaO, SiO , and CO . First of all, MgO was considered only in the Fm3m phase, since 2 2 another investigation has shown that it remains in that phase up to 227 GPa . CO was considered only in the I42d phase . We performed static calculations for several other CO phases and we obtained that the above mentioned phase is the most stable for pressures above 10 GPa, which is fully consistent with results of another investigation that found this phase for pressures from 19 to 150 GPa . Additionally, ignoring those other low pressure stable phases does not compromise our conclusions on the properties of Earth’s lower mantle. For CaSiO , we considered only the tetragonal (I/4mnm) phase for all the pressures and temperatures studied here , since it is the most important phase in the thermodynamic conditions of interest. Moreover, within such conditions, our methodology provides results which comply with the validity criteria described in section II B, i.e. the thermal expansion coeﬃcient of this phase does not diverge up to 3000 K. Figure 2 shows the stable crystalline phases for all other minerals considered in this investigation, as function of pressure at temperatures of 300 and 2000 K, in which there are available experimental data for comparison. For CaCO at 300 K between 0 and 150 GPa, the material goes from Pmcn to P2 /c-l at 15 GPa, from P2 /c-l to Pmmn at 42 GPa, and 1 1 from Pmmn to P2 /c-h at 51 GPa, results which are consistent with another theoretical investigation . Figure 2 also shows that, as the temperature increases, there is a major 6 change on the transition pressures. At a high temperature, 2000 K, the Pmmn to P2 /c-h phase transition occurs at about 80 GPa, in agreement with the theoretical and experimental 13 38 values of 76 GPa and 105 ± 5 GPa , respectively. For MgCO at 300 K, the material goes from R3c to P1 at 68 GPa, and from P1 to C2/m at 110 GPa, all results in good agreement with other theoretical investigations using static calculations . Our results indicate that at temperatures over 1850 K, this mineral follows a direct transition from R3c to C2/m, i.e. above that temperature the P1 is not stable at any pressure. Therefore, our results indicate that P1 should be of low geophysical interest for studies of the lower mantle properties. For MgSiO at 300 K, the mineral goes from Pbnm to Cmcm at 94 GPa, which is consistent with another theoretical investigation . At 2000 K, this transition occurs at 111 GPa, in good agreement with the experimental value of 120 ± 3 GPa . 3m to Pm3m. The transition between these For CaO at 300 K, the material goes from Fm two phases is found at 55 GPa, in good agreement with experimental data that identiﬁed this transition between 59.8 and 63.2 GPa . For SiO at 300 K (2000 K), the material goes from P4 /mnm to Pnnm at 48 GPa (62 2 2 GPa), and from Pnnm to Pbnc at 88 GPa (99 GPa). Those transition values are in good 42,43 44,45 agreement with experimental data at 300 K and 2000 K . III. RESULTS A. Decomposition of carbonates at high pressures and temperatures Most of the Earth’s oxidized carbon is expected to be harbored by Mg and/or Ca car- bonate forms under mantle pressures and temperatures. In order to identify the potential carbon hosts in this region, we initially explore the energetics associated with the following decomposition reactions: CaCO ⇒ CaO + CO (R1) 3 2 MgCO ⇒ MgO + CO (R2) 3 2 Figure 3 shows the relative Gibbs free energy per unit formula (u.f.) as a functions of pressure and temperature. It shows that the direct decompositions of CaCO and MgCO 3 3 7 into their respective alkaline earth oxides plus CO are unfavorable all over the lower mantle. CaCO and MgCO show similar trends with pressure, i.e., increasing the pressure reduces 3 3 the Gibbs free energy diﬀerence for decomposition, which are positive for all the pressures of interest of the lower mantle. However, the energy cost for the CaCO reaction is much higher than the one for MgCO . It should be stressed that our static results, black lines in ﬁgure 3, are in very good 12,13 agreement with the static results presented in other recent investigations , as well as with experimental data . Those investigations, without taking into account temperature eﬀects, have suggested that free CO does not occur as an independent phase within the Earth’s mantle. According to the results presented in ﬁgures 3(a) and (b), temperature eﬀects do not alter the relative stability of carbonates, when compared to their most elementary constituents. In reality, the phenomenology is quite the contrary, as the ﬁgure shows, a temperature increase further reduces the possibility of free CO to exist in that region. B. Stability of carbonates under excess of SiO We now evaluate the stability of carbonates in a condition of excess of SiO , as represented by reactions (R3) and (R4). This condition is particularly important when one takes into account that the upper and lower mantle have material mixing, resulting from the basaltic part of subducting slabs, which is rich in SiO . MgCO + SiO ⇒ MgSiO + CO (R3) 3 2 3 2 CaCO + SiO ⇒ CaSiO + CO (R4) 3 2 3 2 Figure 4 shows the relative Gibbs free energies for Mg and Ca carbonates to transform into their respective silicates. The results indicated that reactions (R3) and (R4) are mostly unfavorable, carbonates do not react to form silicates within those thermodynamic condi- tions. Therefore, there should be no CO formation in the lower mantle as result of those reactions. Static results, in ﬁgure 4(a) indicate that reaction (R3) is unfavorable up to 127 GPa, i.e. MgCO is more stable than MgSiO . This pressure is lower than the one in the 3 3 lower mantle-core boundary of 136 GPa, such that static results suggested that reaction (R3) would be favorable at the bottom of the lower mantle, which could lead to the generation 8 of free CO . The ﬁgure shows that reaction (R3) becomes less favorable with increasing temperature, such that this reaction seems unfavorable at typical temperatures of the lower mantle bottom. Our static results for reaction (R3), black line in ﬁgure 4(a), are in good agreement with 12,13 static results presented in recent theoretical investigations . However, the transition pressure found here is lower than the one found by those authors. Such diﬀerences could be explained by the functionals used to describe electron-electron interactions, since our 20 21 investigation used the LDA , while those investigations use the GGA . The choice of LDA over GGA in this investigation was discussed in detail in section II A. According to ﬁgure 4 (b), the reaction (R4) indicates that CaCO is more stable than CaSiO up to 150 GPa at any temperature. Since this pressure is much larger than the ones at the bottom of the lower mantle, this reaction would not occur in that region. Our static results on the transition pressure for reaction (R4) disagree with another theoreti- cal investigation , which predicted a reaction of calcium carbonate with SiO within the pressure range of the lower mantle. In fact, that investigation did not take into account the P2 /c-h phase of CaCO , which may explain the diﬀerences in the conclusions. On the other 1 3 hand, our static results are in good agreement with the ones of a more recent investigation . It should be pointed out that our investigation carries some uncertainties on the transition pressures of SiO , as shown in ﬁgure 2, which could play some role on the ﬁnal conclusions about the reactions (R3) and (R4), particularly on the properties near the core-mantle boundary. C. Stability of carbonates under excess of MgO or MgSiO It is reasonably well established that MgSiO and MgO are the two major minerals in the lower mantle. Therefore, it is important to study the stability of MgCO or CaCO 3 3 carbonates in a rich environment with those two major minerals. To explore those conditions, we evaluate the energetics associated to reactions (R5) and (R6). MgCO + CaO ⇒ CaCO + MgO (R5) 3 3 CaCO + MgSiO ⇒ MgCO + CaSiO (R6) 3 3 3 3 9 Figure 5(a) shows that the reaction (R5) is energetically favorable, i.e. CaCO + MgO is more stable than MgCO + CaO, at lower mantle conditions. Although by increasing the temperature, the relative Gibbs free energy is reduced, this eﬀect is not large enough to change the stability of MgCO . Therefore, when there is MgO in excess, as expected in a pyrolitic mantle, CaCO is the stable carbonate under those thermodynamic condi- tions. Additionally, our static results are in good agreement with recent results from other 12,13 investigations . We now explore the conditions of MgSiO excess, which is the main mineral in the lower mantle. Figure 5(b) shows that reaction (R6) presents a much richer phenomenology than reaction (R5) across the pressure range of interest. Our static calculations show that MgCO + CaSiO is more favorable than CaCO + MgSiO at low pressures, but at pressures higher 3 3 3 than 75 GPa, then CaCO + MgSiO becomes favorable. However, as the ﬁgure shows, there 3 3 is a dramatic change in the behavior at high temperatures. An increase in temperature increases the pressure in which such reaction could be favorable. For typical lower mantle temperatures, in the order of at least 2000 K, the reaction (R6) is not favorable anymore. Therefore, the results in ﬁgure 5(b) indicated that, at lower mantle conditions under excess of MgSiO , carbonates appear to be favorable in the form of MgCO . 3 3 This conclusion is consistent with an experimental investigation , which suggested MgCO as the main oxidized carbon host in Earth’s mantle. Figure 6 shows the phase diagram for reaction (R6), along with some Earth’s 15,46–48 geotherms . According to the ﬁgure, all geotherms lie in the region of stability of MgCO , indicating this mineral as the most likely carbon host in the lower mantle, as ex- pected in a pyrolitic mantle. However, it should be pointed that there is still controversy on the temperatures at the lower mantle bottom, therefore, the stability of CaCO over MgCO 3 3 would require a much colder mantle. IV. SUMMARY In summary, this investigation explored the stability of MgCO and CaCO carbonates 3 3 in the thermodynamic conditions of Earth’s lower mantle. The results indicated that a direct decomposition of those carbonates is unfavorable at high temperatures and pressures. Assuming an iron free pyrolitic lower-mantle composition, the investigation also explored 10 the stability of those carbonates in conditions of excess of SiO , MgO, and MgSiO . 2 3 In a MgO-rich environment, the results showed that calcium carbonate is more stable than magnesium carbonate, and the relative energy reduction with increasing temperature was not large enough to change the stability of CaCO . On the other hand, in a MgSiO - 3 3 rich environment, the static results showed that at the upper half of the lower mantle, the MgCO + CaSiO reaction is more favorable than the CaCO + MgSiO one, while at 3 3 3 3 pressures higher than 75 GPa, the latter becomes favorable. However, when the eﬀects of temperature are taken into account, this behavior changes dramatically and the increase in temperature increases the pressure at which this reaction becomes favorable, and ﬁnally the magnesium carbonate turns out to be more stable than the calcium carbonate. Since magnesium silicate is the main component of a pyrolitic mantle, it can be inferred that carbonates appear to be favorable in the form of MgCO . Therefore, the magnesium carbonate should be the main host of oxidized carbon in most of the lower mantle. Only in the bottom of the mantle or in a region with MgO in excess that calcium carbonate could become preferable. However, in the bottom of the mantle, this carbonate would be favorable only considering geotherms with very low increase in temperature close to the core-mantle boundary. The results also showed that both carbonates do not decompose into their respective alkaline oxides plus CO through the entire lower mantle, which indicates low concentration of free carbon dioxide in those regions. However, the decomposition reaction of MgCO into MgO+CO would only be possible very close to the core-mantle boundary. Furthermore, free CO could be produced in an environment with excess of SiO , as in small silica-rich 2 2 basaltic parts of the subducted slabs. All those results add new evidences for the presence of carbon on deep mantle in the 1,11 form of carbonates . However, it is still uncertain how the presence of iron, and its rich 49,50 phenomenology associated to spin transition at high pressures , would aﬀect the stability of carbonates in the lower mantle. Moreover, the decomposition of CO into diamond plus oxygen should be explored to enrich the discussion of the presence of carbon in the deep mantle in a reduced state. Acknowledgments This investigation was supported by the Brazilian agencies CNPq and CAPES. We ac- 11 knowledge resources from the Blue Gene/Q supercomputer supported by the Center for Re- search Computing (Rice University) and the Superintendˆencia de Tecnologia da Informac¸˜ao (Universidade de S˜ao Paulo). R. M. Hazen, C. M. Schiﬀries, Why deep carbon?, Rev. Mineral. Geochem. 75 (2013) 1. B. Marty, C. M. O. Alexander, S. N. Raymond, Primordial origins of Earth’s carbon, Rev. Mineral. Geochem. 75 (2013) 149. M. J. Walter, S. C. Kohn, D. Araujo, G. P. Bulanova, C. B. Smith, E. Gaillou, J. Wang, A. Steele, S. B. Shirey, Deep mantle cycling of oceanic crust: evidence from diamonds and their mineral inclusions., Science 334 (2011) 54–7. D. G. Pearson, F. E. Brenker, F. Nestola, J. McNeill, L. Nasdala, M. T. Hutchison, S. Matveev, K. Mather, G. Silversmit, S. Schmitz, B. Vekemans, L. Vincze, Hydrous mantle transition zone indicated by ringwoodite included within diamond., Nature (London) 507 (2014) 221. F. Maeda, E. Ohtani, S. Kamada, T. Sakamaki, N. Hirao, Y. Ohishi, Diamond formation in the deep lower mantle: a high-pressure reaction of MgCO and SiO ., Sci. Rep-UK 7 (2017) 40602. 3 2 A. R. Oganov, R. J. Hemley, R. M. Hazen, A. P. Jones, Structure, bonding, and mineralogy of carbon at extreme conditions, Rev. Mineralog. Geochem. 75 (2013) 47–77. F. E. Brenker, C. Vollmer, L. Vincze, B. Vekemans, A. Szymanski, K. Janssens, I. Szaloki, L. Nasdala, W. Joswig, F. Kaminsky, Carbonates from the lower part of transition zone or even the lower mantle, Earth Planet. Sc. Lett. 260 (2007) 1–9. A. P. Jones, M. Genge, L. Carmody, Carbonate melts and carbonatites, Rev. Mineralog. Geochem. 75 (2013) 289. S. B. Shirey, P. Cartigny, D. J. Frost, S. Keshav, F. Nestola, P. Nimis, D. G. Pearson, N. V. Sobolev, M. J. Walter, Diamonds and the geology of mantle carbon, Rev. Mineralog. Geochem. 75 (2013) 355. M. R. Burton, G. M. Sawyer, D. Granieri, Deep carbon emissions from volcanoes, Rev. Miner- alog. Geochem. 75 (2013) 323–354. M. L. Marcondes, J. F. Justo, L. V. C. Assali, Carbonates at high pressures: possible carriers for deep carbon reservoirs in the Earth’s lower mantle, Phys. Rev. B 94 (2016) 104112. 12 12 A. R. Oganov, S. Ono, Y. Ma, C. W. Glass, A. Garcia, Novel high-pressure structures of MgCO , CaCO and CO and their role in Earth’s lower mantle, Earth Planet. Sci. Lett. 273 (2008) 3 2 38–47. C. J. Pickard, R. J. Needs, Structures and stability of calcium and magnesium carbonates at mantle pressures, Phys. Rev. B 91 (2015) 104101. G. Fiquet, F. Guyot, M. Kunz, J. Matas, D. Andrault, M. Hanﬂand, Structural reﬁnements of magnesite at very high pressure, Am. Mineralog. 87 (2002) 1261–1265. M. Isshiki, T. Irifune, K. Hirose, S. Ono, Y. Ohishi, T. Watanuki, E. Nishibori, M. Takata, M. Sakata, Stability of magnesite and its high-pressure form in the lowermost mantle., Nature 427 (2004) 60–63. P. N. Gavryushkin, N. S. Martirosyan, T. M. Inerbaev, Z. I. Popov, S. V. Rashchenko, A. Y. Likhacheva, S. S. Lobanov, A. F. Goncharov, V. B. Prakapenka, K. D. Litasov, Aragonite-II and CaCO -VII: new high-pressure, high-temperature polymorphs of CaCO ., Cryst. Growth 3 3 Des. 17 (2017) 6291–6296. D. Smith, K. V. Lawler, M. Martinez-Canales, A. W. Daykin, Z. Fussell, G. A. Smith, C. Childs, J. S. Smith, C. J. Pickard, A. Salamat, Postaragonite phases of CaCO at lower mantle pressures, Phys. Rev. Materials 2 (2018) 013605. M. Shirasaka, E. Takahashi, Y. Nishihara, K. Matsukage, T. Kikegawa, In situ X-ray observation of the reaction dolomite = aragonite + magnesite at 900-1300 K, Am. Mineralog. 87 (2002) 922. P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococcioni, I. Dabo, A. Dal Corso, S. de Gironcoli, S. Fabris, G. Fratesi, R. Gebauer, U. Gerstmann, C. Gougoussis, A. Kokalj, M. Lazzeri, L. Martin-Samos, N. Marzari, F. Mauri, R. Mazzarello, S. Paolini, A. Pasquarello, L. Paulatto, C. Sbraccia, S. Scandolo, G. Sclauzero, A. P. Seitsonen, A. Smogunov, P. Umari, R. M. Wentzcovitch, QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materi- als., J. Phys.: Condens. Mat. 21 (2009) 395502. W. Kohn, L. J. Sham, Self-consistent equations including exchange and correlation eﬀects, Phys. Rev. 140 (1965) A1133–A1138. J. P. Perdew, K. Burke, M. Ernzerhof, Generalized gradient approximation made simple., Phys. Rev. Lett. 77 (1996) 3865–3868. M. L. Marcondes, R. M. Wentzcovitch, L. V. C. Assali, Importance of van der Waals interaction 13 on structural, vibrational, and thermodynamic properties of NaCl, Solid State Commun. 273 (2018) 11–16. P. E. Blochl, Projector augmented-wave method, Phys. Rev. B 50 (1994) 17953–17979. N. A. W. Holzwarth, A. R. Tackett, G. E. Matthews, A projector augmented wave (PAW) code for electronic structure calculations, part I: atompaw for generating atom-centered functions, Comput. Phys. Commun. 135 (2001) 329–347. R. M. Wentzcovitch, J. L. Martins, G. D. Price, Ab initio molecular dynamics with variable cell shape: application to MgSiO , Phys. Rev. Lett. 70 (1993) 3947–3951. S. Baroni, S. de Gironcoli, A. Dal Corso, P. Giannozzi, Phonons and related crystal properties from density-functional perturbation theory, Rev. Mod. Phys. 73 (2001) 515–562. Y. Wang, J. J. Wang, W. Y. Wang, Z. G. Mei, S. L. Shang, L. Q. Chen, Z. K. Liu, A mixed- space approach to ﬁrst-principles calculations of phonon frequencies for polar materials, J. Phys.: Condens. Mat. 22 (2010) 202201. R. M. Wentzcovitch, Y. G. Yu, Z. Wu, Thermodynamic properties and phase relations in mantle minerals investigated by ﬁrst principles quasiharmonic theory, Rev. Mineralog. Geochem. 71 (2010) 59. R. M. Wentzcovitch, B. B. Karki, M. Cococcioni, S. de Gironcoli, Thermoelastic properties of MgSiO -perovskite: insights on the nature of the Earth’s lower mantle, Phys. Rev. Lett. 92 (2004) 18501. P. Carrier, R. Wentzcovitch, J. Tsuchiya, First-principles prediction of crystal structures at high temperatures using the quasiharmonic approximation, Phys. Rev. B 76 (2007) 064116. K. D. Litasov, A. Shatskiy, P. N. Gavryushkin, A. E. Bekhtenova, P. I. Dorogokupets, B. S. Danilov, Y. Higo, A. T. Akilbekov, T. M. Inerbaev, P-V-T equation of state of CaCO aragonite to 29 GPa and 1673 K: In situ X-ray diﬀraction study, Phys. Earth Planet. Int. 265 (2017) 82. P. I. Dorogokupets, Equation of state of magnesite for the conditions of the earth’s lower mantle, Geochem. Int+ 45 (2007) 561–568. L. Li, D. J. Weidner, J. Brodholt, D. Alf`e, G. D. Price, R. Caracas, R. Wentzcovitch, Elasticity of CaSiO perovskite at high pressure and high temperature, Phys. Earth Planet. Int. 155 (2006) 249–259. B. B. Karki, R. M. Wentzcovitch, S. de Gironcoli, S. Baroni, Ab initio lattice dynamics of MgSiO perovskite at high pressure, Phys. Rev. B 62 (2000) 14750–14756. 14 35 T. S. Duﬀy, R. J. Hemley, H.-k. Mao, Equation of state and shear strength at multimegabar pressures: magnesium oxide to 227 GPa, Phys. Rev. Lett. 74 (1995) 1371–1374. F. Datchi, B. Mallick, A. Salamat, S. Ninet, Structure of Polymeric Carbon Dioxide CO −V, Phys. Rev. Lett. 108 (2012) 125701. S. Shim, R. Jeanloz, T. S. Duﬀy, Tetragonal structure of CaSiO perovskite above 20 GPa, Geophys. Res. Lett. 29 (2002) 2166. S. S. Lobanov, X. Dong, N. S. Martirosyan, A. I. Samtsevich, V. Stevanovic, P. N. Gavryushkin, K. D. Litasov, E. Greenberg, V. B. Prakapenka, A. R. Oganov, A. F. Goncharov, Raman spectroscopy and x-ray diﬀraction of sp CaCO at lower mantle pressures, Phys. Rev. B 96 (2017) 104101. T. Tsuchiya, J. Tsuchiya, K. Umemoto, R. M. Wentzcovitch, Phase transition in MgSiO per- ovskite in the earth’s lower mantle, Earth Planet. Sci. Lett. 224 (2004) 241 – 248. M. Murakami, K. Hirose, K. Kawamura, N. Sata, Y. Ohishi, Post-perovskite phase transition in MgSiO ., Science 304 (2004) 855. T. Yamanaka, K. Kittaka, T. Nagai, B1-B2 transition in CaO and possibility of CaSiO - perovskite decomposition under high pressure, J. Miner. Petrol. Sci. 97 (4) (2002) 144–152. R. J. Hemley, J. Shu, M. A. Carpenter, J. Hu, H. K. Mao, K. J. Kingma, Strain/order parameter coupling in the ferroelastic transition in dense SiO , Solid State Commun. 114 (2000) 527–532. D. Andrault, G. Fiquet, F. Guyot, M. Hanﬂand, Pressure-induced landau-type transition in stishovite, Science 282 (1998) 720–724. S. Ono, K. Hirose, M. Murakami, M. Isshiki, Post-stishovite phase boundary in SiO determined by in situ X-ray observations, Earth Planet. Sci. Lett. 197 (2002) 187–192. M. Murakami, K. Hirose, S. Ono, Y. Ohishi, Stability of CaCl -type and α-PbO -type SiO at 2 2 2 high pressure and temperature determined by insitu X-ray measurements, Geophys. Res. Lett. 30 (2003) 1207. J. M. Brown, T. J. Shankland, Thermodynamic parameters in the earth as determined from seismic proﬁles, Geophys. J. Int. 66 (1981) 579–596. R. Boehler, Highpressure experiments and the phase diagram of lower mantle and core materials, Rev. Geophys. 38 (2000) 221–245. O. L. Anderson, The earth’s core and the phase diagram of iron, Philos. T. R. Soc. A 306 (1982) 21–35. 15 49 R. M. Wentzcovitch, J. F. Justo, Z. Wu, C. R. S. da Silva, D. A. Yuen, D. Kohlstedt, Anomalous compressibility of ferropericlase throughout the iron spin cross-over., Proc. Nat. Acad. Sci. USA 106 (2009) 8447–8452. J. Liu, J.-F. Lin, V. B. Prakapenka, High-pressure orthorhombic ferromagnesite as a potential deep-mantle carbon carrier, Sci. Rep-UK 5 (2015) 7640. 16 4 (a) (b) Aragonite 1 Post Aragonite 0 0 0 1000 2000 3000 0 1000 2000 3000 Temperature (K) Temperature (K) FIG. 1: Thermal expansion coeﬃcient at several pressures as a function of temperature for CaCO in (a) aragonite and (b) post-aragonite phases. The QHA boundary is deﬁned by the position of the inﬂection points of α(P,T), discussed in section II B. Experimental results at 0 GPa are presented with black symbols . 60 GPa 55 GPa 15 GPa 50 GPa 45 GPa 40 GPa 12 GPa 9 GPa 6 GPa 3 GPa 0 GPa −5 −1 α (10 K ) −5 −1 α (10 K ) CaCO Pmcn Pmmn P2 /c−l P2 /c−h 1 1 MgCO − − C2/m R3c P1 MgSiO Pbnm Cmcm CaO − − Fm3m Pm3m SiO Pnnm Pbnc P4 /mnm 0 20 40 60 80 100 120 140 P (GPa) FIG. 2: Stable crystalline phases, as function of pressure, of CaCO , MgCO , MgSiO , CaO, 3 3 3 and SiO materials considered in this investigation, computed with theoretical approximations presented in section II A. The ﬁgure presents results for two temperatures: 300 K (orange line) and 2000 K (blue line), with the respective transitions represented by vertical lines. Experimental values for phase transitions, with respective experimental error bars, are presented with brown 38 40 41 42–45 symbols for: CaCO , MgSiO , CaO , and SiO . 3 3 2 18 3 (a) . CaCO CaO + CO 3 2 2.5 static 1.5 CaO + CO 2 0 K 300 K 1 1000 K 2000 K CaCO 3 3000 K 0.5 −0.5 0 20 40 60 80 100 120 140 P (GPa) 2.5 MgCO MgO + CO 3 2 (b) 2 static 0 K 300 K MgO + CO 1.5 1000 K 2000 K 3000 K 0.5 MgCO −0.5 0 20 40 60 80 100 120 140 P (GPa) FIG. 3: The relative Gibbs free energies per u.f. as function of pressure, at several temperatures, for (a) CaCO → CaO + CO and (b) MgCO → MgO + CO reactins. The dotted black lines 3 2 3 2 represent the respective relative enthalpies. The vertical gray lines indicate the pressures at the top (23 GPa) and bottom (136 GPa) of the lower mantle. The kinks in the curves arise from phase transitions that occur in the minerals described in section II C and shown in ﬁgure 2, which are used to explore the dissociation reactions. ΔG (eV) ΔG (eV) 2.5 (a) MgCO + SiO MgSiO + CO 3 2 3 2 static 1.5 0 K MgSiO + CO 3 2 300 K 1000 K 2000 K 3000 K 0.5 MgCO + SiO 3 2 −0.5 0 20 40 60 80 100 120 140 P (GPa) 2.5 (b) CaCO + SiO CaSiO + CO 3 2 3 2 static 0 K 1.5 300 K CaSiO + CO 1000 K 3 2 2000 K 3000 K 0.5 CaCO + SiO 3 2 −0.5 0 20 40 60 80 100 120 140 P (GPa) FIG. 4: The relative Gibbs free energy per u.f. as function of pressure, at several temperatures, for (a) MgCO + SiO → MgSiO + CO and (b) CaCO + SiO → CaSiO + CO reactions. 3 2 3 2 3 2 3 2 The dotted lines represent the relative enthalpies for the reactions. The vertical gray lines indicate the pressures at the top and bottom of the lower mantle. ΔG (eV) ΔG (eV) 0.5 MgCO + CaO CaCO + MgO 3 3 (a) static MgCO + CaO 0 K −0.5 300 K 1000 K 2000 K 3000 K −1 −1.5 CaCO + MgO −2 0 20 40 60 80 100 120 140 P (GPa) 0.5 (b) CaCO + MgSiO MgCO + CaSiO 3 3 3 3 CaCO + MgSiO 3 3 static 0 K −0.5 300 K MgCO + CaSiO 3 3 1000 K 2000 K 3000 K −1 0 20 40 60 80 100 120 140 P (GPa) FIG. 5: The relative Gibbs free energy per u.f. as function of pressure, at several temperatures, for (a) MgCO + CaO → CaCO + MgO and (b) CaCO + MgSiO → MgCO + CaSiO reactions. 3 3 3 3 3 3 The black dotted lines represent the relative enthalpies in (a) and (b). The vertical gray lines indicate the pressures at the top and bottom of the lower mantle. ΔG (eV) ΔG (eV) 3000 3000 2500 2500 2000 2000 1500 1500 MgCO MgCO + CaSiO + CaSiO 3 3 3 3 1000 1000 500 500 CaCO CaCO + MgSiO + MgSiO 3 3 3 3 0 0 60 60 80 80 100 100 120 120 140 140 P (GPa) P (GPa) FIG. 6: Phase diagram for the stability (brown line) of MgCO + CaSiO versus CaCO + MgSiO . 3 3 3 3 15 46 47 48 The ﬁgure also shows several geotherms: thick gray , green , blue , and red lines. Geotherm T (K) T (K)
Condensed Matter – arXiv (Cornell University)
Published: Apr 19, 2018
Access the full text.
Sign up today, get DeepDyve free for 14 days.