Research Article | Open Access
Lei Liu, Longxing Yang, Chunqiang Zhuang, Guangshu Yang, Li Yi, Hong Liu, Fengxia Sun, Xiaoyu Gu, Hanyu Wang, "Diffusion of CO2 in Magnesite under High Pressure and High Temperature from Molecular Dynamics Simulations", Geofluids, vol. 2021, Article ID 6621425, 9 pages, 2021. https://doi.org/10.1155/2021/6621425
Diffusion of CO2 in Magnesite under High Pressure and High Temperature from Molecular Dynamics Simulations
CO2 transports in the Earth’s interior play a crucial role in understanding the deep carbon cycle and the global climate changes. Currently, CO2 transports inside of the Earth under extreme condition of pressure and temperature have not been understood well. In this study, the molecular dynamics (MD) calculations were performed to study CO2 transports under different CO2 pressures in slit-like magnesite pores with different pore sizes at 350~2500 K and 3~50 GPa are presented. Diffusion of CO2 in magnesite was improved as the temperature increases but showed the different features as a function of pressure. The diffusion coefficients of CO2 in magnesite were found in the range of . Magnesite with the pore size of 20~25 Å corresponds to the highest transports. Anisotropic diffusion of CO2 in magnesite may help to understand the inhomogeneous distribution of carbon in the upper mantle. The time of CO2 diffusion from the mantle to Earth surface was estimated to be around several tens of Ma and has an important effect on deep carbon cycle. The simulation of CO2 transports based on the Earth condition provides new insights to revealing the deep carbon cycle in the Earth’s interiors.
Prediction of the changes in the Earth’s climate due to the increase of carbon dioxide (CO2) concentrations needs to deeply understand the chemical and physical processes that control the carbon cycle . Generally, the carbon cycle involves the deep and shallow subcycles of carbon. So far, intensive attention has been paid to the shallow carbon cycle that is mainly related to the atmosphere, oceans, and shallow crustal. However, there is relatively little known about the deep carbon cycle. Due to the Earth’s deep interior that contains more than 90% of the total Earth’s carbon [2, 3], the carbon transport and storage in the deep Earth play a key role in determining the global carbon cycle [1, 4–9]. Therefore, it is very necessary to understand the carbon cycle in the Earth’s deep interior.
Carbon in the mantle exists in the different oxidized forms (e.g., CO, CO2, and CO32−) and the different states (solid, molten carbonate, fluid, etc.) [10–13]. Magnesite (MgCO3) is the most common carbonate phase in the subducted slabs . High-pressure experiments suggest that magnesite with a rhombohedral structure can be stable up to 115 GPa at 2000~3000 K which may be the major carbon host in the mid-lower mantle [14–17].
CO2 is an important carbon-bearing form in the Earth’s interior. About 40 megatons of CO2 emits from the mantle to the surface every year by diffuse outgassing . Therefore, CO2 plays a fundamental role in determining the global carbon cycle. In situ Raman spectroscopy in laser-heated diamond anvil cells revealed that CO2 molecular was stable up to 34 GPa at ~1770 K [6, 19]. Therefore, the transport properties of CO2 in the deep Earth under high pressure and high temperature are a vital link to understand the deep carbon cycle.
Except for high-pressure and high-temperature (high P-T) experiments, molecular dynamics (MD) simulations are powerful tools to investigate the molecular-scale dynamic properties of materials ([20–26]; Hu et al. 018). There are some achievements on researching diffusion of CO2 in different materials by MD simulation . For example, the CO2 diffusion coefficients of crystalline polystyrene at 1 atm. are (at 333.15 K), (at 343.15 K), and (at 353.15 K) . Diffusion coefficients of supercritical CO2 and its mixtures monotonously decreased as its increasing molar density at 317.5 K . The diffusion coefficients of H2, CH4, CO, O2, and CO2 in water near the critical point (600~670 K, 250 atm) were investigated by MD simulation , and CO2 can greatly influence the diffusion of mobile species in clay minerals . However, most of those works focused on model structures, which are different from the potential minerals in the mantle. In addition, previous works were mainly performed under the static conditions without concerning high temperature and high pressure. Accordingly, the transport properties of CO2 under high temperature-pressure in the Earth’s interior conditions have not yet been identified. In this work, the diffusion of CO2 in magnesite was investigated under the high temperature-pressure conditions mimicking the extreme environment in the mantle.
2.1. Crystalline Models
Magnesite has a space group of R3c. Its structure is composed of the alternative layers of Mg and CO3 groups. Among them, one C atom and three neighbor O atoms lie in the same plane that is parallel to the - plane. Under the ambient condition, the cell parameters of magnesite are as follows: , , , . A supercell of magnesite was established. A vacuum layer was added between two CO3 planes to build slit-like pores. The pore size (H, Å) was determined by the height of a vacuum layer. The two-layer slit-like pore structures of magnesite with electrically neutral were built and simulated (Figure 1).
2.2. Equilibrium Adsorption Configurations: The Grand Canonical Monte Carlo Method (GCMC) Simulation
The GCMC employed the μVT ensemble. The chemical potential , the volume (), and the temperature () were kept constant, while the number of molecules was allowed to fluctuate . GCMC simulations were carried out by sorption code as implemented in Materials Studio 2017. The structures of CO2 in magnesite pores with energy minimized adsorption were obtained under the different conditions. During the GCMC simulation, CO2 molecules within the framework were randomly rotated and translated. Only were the positions and orientations of CO2 molecules sampled. Each configuration was treated as a rigid body. The COMPASS force field was used to describe the hypersurface of the potential energy on which the atomic nuclei moved . The Ewald method was adopted to calculate the Coulomb force interaction. The atom interaction-based method was used to calculate van der Waals force with the Lennard-Jones potential. The cut-off distance was 0.15 Å. The number of load step in each calculated run is , in which the balance step number is and the process step number is . This GCMC method was successfully used to calculate the absorption of CO2 in magnesite in our previous works , and its validity also was proved by different works [31, 36–40].
2.3. Molecular Dynamics (MD) Simulations
MD simulations were carried out by the Forcite code as implemented in Materials Studio 2017. The equilibrium configurations from GCMC calculations were further used in MD simulations. The MD simulations were performed with a time step of 1 fs (up to 1000 ps) using an NPT ensemble at 350~2500 K and 3~50 GPa. Interlayer molecular configurations were sampled every 1000 steps for further analysis. The temperature was kept constant using a velocity scale thermostat algorithm, and the convergence criterion is 1 K. The Ewald summation method was employed with an Ewald accuracy of kcal/mol. Like the GCMC simulations, the COMPASS force field was used for describing the potential energy hypersurface .
2.4. Calculations for Diffusivities
The mean square displacement (MSD) is a mode of the displacement of particles followed over time. It can be defined aswhere is the position of the atom at time , and the angular bracket denotes that the quantity is an ensemble average property. In an equilibrium ensemble, the average squared displacement is independent of the time , which can be averaged out, leaving the mean square displacement over an interval :
If there are equivalent particles, the MSD can be further averaged as
If the particle is diffusing, the MSD becomes linear as a function of time (diffusive regime). The self-diffusion coefficient () can be described by the slope according to Einstein equation:where is the number of diffusion molecules, is the time, and is the displacement vector of the th molecule at time . The angular bracket denotes that the quantity is an ensemble average property.
3.1. The Diffusivities at High Pressure and High Temperature
Diffusions of CO2 in magnesite pore were calculated with the pore size of 25 Å and the CO2 pressure of 50 MPa at high pressures (3~50 GPa) and high temperatures (350~2500 K); the results are shown in Figure 2. Diffusion of CO2 in magnesite increased as temperature increases but decreased as pressure increases. Specifically, the ( in 10-12 m2/s) increased from 1.67 at 350 K to 4.45 at 2500 K (3 GPa), while it decreased from 4.45 at 3 GPa to 3.29 at 50 GPa (2500 K). However, the change trend of the diffusion with temperature is different under different pressures. At low pressure (<9 GPa), the diffusion coefficient consistently increased with increasing the temperature at the whole range explored in this work. However, as the pressure increased, the diffusion coefficients showed the different trends when the temperature is higher or lower than 900 K.
For the temperature higher than 900 K, the values showed a good linear relation with the , which corresponds well with the Arrhenius law:where is the preexponential factor, is the gas constant, and is the activation energy for diffusion. Fitting parameters at 900~2500 K and 3~50 GPa are shown in Figure 3 and Table 1. The decreased from 4.929 kJ/mol at 3 GPa to 4.501 kJ/mol at 50 GPa. The increased from 25.2 kJ/mol at 3 GPa to 60.7 kJ/mol at 50 GPa, indicating that diffusion of CO2 in magnesite becomes more difficult under high pressure.
The different diffusion properties of CO2 in magnesite under different pressure and temperature (9 GPa and 900 K) might be related to the phase transition of CO2 and stability of magnesite under different pressure and temperature. The CO2 might transform from fluid phase to Pa3 molecular structure, Pbcn, P42/mnm, and I42d with increasing pressure [41–44]; however, the stability of CO2 under high temperature is not yet well understood. For example, at ambient temperature, Pa3 CO2 undergoes a pressure-induced transformation to the Cmca symmetry at about 9-12 GPa; when pressure is larger than 12 GPa, there are two molecular phases: P42/mnm phase for K and Pbcn phase for K [43, 44]. At the same time, the presence of CO2 might affect the stability of magnesitite [45, 46].
3.2. Diffusion Anisotropy
Figure 4 shows an anisotropic feature of CO2 diffusion in magnesite. For the magnesite structure, it has the same -axis and -axis; therefore, the diffusion in the - plane of magnesite is the same. Accordingly, we only calculated diffusion coefficient along -axis and -axis direction under high T-P. The diffusion coefficients () of CO2 along the -axis and -axis varied differently with temperature and pressure (Figure 4). At temperature of 350~500 K, the ( along -axis direction) was larger than the ( along -axis direction) at a low pressure (<9 GPa) and then the became smaller than the with increasing pressure. When temperature was higher than 700 K, the was larger than the , but the differences became very small with increasing pressure. Anisotropy of the CO2 diffusion in magnesite decreased with increasing pressure and temperature. At 25 GPa, anisotropy increased from 3% at 2500 K to 60% at 350 K. For comparison, at 3 GPa, the anisotropy increased from 3.5% at 2500 K to 280% at 350 K. As the pressure was higher than 25 GPa, the difference between and was very small. The results indicated diffusion of CO2 in magnesite is anisotropic at a low pressure and is isotropic under high pressure.
Anisotropy of CO2 diffusion in magnesite is important to understand the anisotropic properties of deep carbon cycle. A study on olivine-hosted melt inclusions from the Mid-Atlantic Ridge indicated that carbon content in the upper mantle is highly heterogeneous, varying by almost two orders of magnitude globally . Volatiles in the mantle can be hosted either in independent phases or as dissolved components in minerals. The content of magnesite in the upper mantle is 1000~8000 ppm [17, 48]. Under the upper mantle ( and ), diffusion of CO2 in magnesite shows an obvious anisotropy. Therefore, this diffusion difference of CO2 probably leads to the heterogeneous transport and storage of carbon in the mantle.
3.3. The Diffusion under Different CO2 Pressures
The diffusion of CO2 in magnesite is closely related to its content in the mantle. Therefore, values of CO2 in magnesite under different CO2 pressures were calculated in order to understand the effect of CO2 content on its diffusion behavior. Figure 5 shows the values of CO2 in magnesite with the CO2 pressures of 10, 30, 50, 70, and 100 MPa under the 900~2500 K and 3~50 GPa. The diffusibility increased as the CO2 pressure increased. The averages of values increased by 1.5%-4.8% as the CO2 pressure increased from 10 MPa to 100 MPa under the different pressures. As detected by Raman spectroscopy, the diffusion coefficient of CO2 in aqueous solutions decreases significantly at high CO2 concentrations . Our results here show different diffusion behavior from the diffusivity of CO2 in aqueous solutions under geologic carbon sequestration conditions. In addition, the effects of CO2 pressure on its diffusion show a slight decrease with increasing the temperature. Compared with the change caused by temperature or pressure (orders of magnitude), the effect of CO2 pressure on its diffusion is limited.
3.4. The Diffusion under Different Magnesite Pore Sizes
Slit-like magnesite pores of 7.5, 10, 12.5, 15, 20, 25, and 30 Å were constructed. Figure 6 shows diffusions of CO2 in these models. It is clear that firstly increases and then decreases with increasing pore sizes. Under the same simulated condition of temperature and pressure, the diffusion has the highest ability in magnesite with a pore size of 20~25 Å. The differences of the diffusion caused by the change of pore size become smaller with increasing pressure.
Plenty of CO2 is emitted from the deep Earth to the surface by diffusion [18, 47] and made significant influence on the Earth’s carbon cycle. Therefore, it is critical to understand the diffusion process of CO2 in the Earth’s interior. The diffusion coefficient of CO2 in magnesite is at 350~2500 K and 3~50 GPa. CO2 diffusion in crystalline polystyrene was at 1 atm and 333.15~353.15 K . The diffusion coefficients of CO2 in silicate melt (Jadeite+albite+quartz) are at 1323~1673 K and 500 MPa . Under the similar condition of temperature and pressure, our calculated diffusion coefficient of CO2 is smaller than that in crystalline polystyrene, however is larger than that of CO2 in silicate melt. The results indicate the faster diffusion ability of CO2 in magnesite crystal than in silicate melt.
Another critical issue on deep carbon cycle is the distribution of the diffusion ability of CO2 along the Earth’s geothermal profile. From Figure 7, we can estimate diffusion distance of CO2 in magnesite is 0.017~0.933 m each year. Accordingly, CO2 diffusion from the mantle to the surface takes around dozens of Ma. Carbon hosted in the mantle ultimately affected the surface carbon budgets over billion-year time scales [50, 51]. Similarly, isotopic results indicated three time-scale carbon cycles in the entire upper mantle: 1: >120 km, 5~10 Ma; 2: 120~410 km, 48~53 Ma; and 3: 410~660 km, >60 Ma [52–56]. The diffusion of CO32- would be three to four orders of magnitude slower than that of molecular CO2 . Therefore, CO2 diffusion plays the important role in determining the deep carbon cycles. For the first time, we present the CO2 diffusion behavior in real minerals (magnesite) under the Earth’s interior conditions. The results here provide new data to understand the deep carbon cycle.
The above results indicate that temperature, pressure, and pore size of magnesite have important effect on CO2 diffusion. The excess adsorption amounts of CO2 in the pores of magnesite decrease with increasing temperature but increase with increasing pressure and pore size . So interior circumstance changes from deep to shallow of the Earth might drive the transport of CO2 in its interior, such as the changes of stress and porosity caused by the movement of fault or tectonic slab. However, the grain or mineral boundaries are always the important channel for fluid migration in the Earth’s interior, which also needs to be paid more attention.
For understanding the Earth’s deep carbon cycle, diffusion of CO2 in magnesite pore with different CO2 pressure and pore sizes was calculated under 350~2500 K and 3~50 GPa. The diffusion coefficient of CO2 in magnesite increases with increasing temperature but decreases with increasing pressure. The increases with increasing pressure, indicating the diffusion of CO2 in magnesite becomes more difficult under high pressure. CO2 has the strongest diffusion in the magnesite when the pore size is 20~25 Å, and the diffusion slightly increases with increasing CO2 pressure. Diffusion of CO2 in magnesite shows obvious anisotropy and might be responsible for the inhomogeneous distribution of carbon in the upper mantle. The diffusion coefficient of CO2 in magnesite is in the range of under the pressure and temperature investigated here. Our calculated results indicate that diffusion of CO2 is faster in magnesite than in silicate melt. The results of CO2 diffusion in magnesite under high pressure and high temperature provide new insights to understand the carbon cycle in the Earth’s interiors.
All data included in this study are available upon request by contact with the corresponding author.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This work was supported by the Foundation of the Key Laboratory of Earthquake Forecasting and United Laboratory of High-Pressure Physics and Earthquake Science, the Institute of Earthquake Forecasting, CEA (Grant No. 2019IEF0101-1, 2019IEF0502, 2017KLEP03, and 2019HPPES06), and the National Natural Science Foundation of China (Grant No. 41762008).
- B. Orcutt, I. Daniel, and R. Dasgupta, Deep Carbon: Past to Present, Cambridge University Press, Cambridge, UK, 2019.
- R. M. Hazen, A. P. Jones, and J. A. Baross, Carbon in Earth, Mineralogical Society of America, 2013.
- M. Javoy, “The major volatile elements of the Earth: their origin, behavior, and fate,” Geophysical Research Letters, vol. 24, no. 2, pp. 177–180, 1997.
- M. R. Burton, G. M. Sawyer, and D. Granieri, “Deep carbon emissions from volcanoes,” Reviews in Mineralogyand Geochemistry, vol. 75, no. 1, pp. 323–354, 2013.
- R. Dasgupta and D. Walker, “Carbon solubility in core melts in a shallow magma ocean environment and distribution of carbon between the Earth’s core and the mantle,” Geochimica et Cosmochimica Acta, vol. 72, no. 18, pp. 4627–4641, 2008.
- K. Spickenbom, M. Sierralta, and M. Nowak, “Carbon dioxide and argon diffusion in silicate melts: insights into the CO2 speciation in magmas,” Geochimica et Cosmochimica Acta, vol. 74, no. 22, pp. 6541–6564, 2010.
- R. Tao, L. Zhang, M. Tian et al., “Formation of abiotic hydrocarbon from reduction of carbonate in subduction zones: constraints from petrological observation and experimental simulation,” Geochimica et Cosmochimica Acta, vol. 239, pp. 390–408, 2018.
- T. Yoshioka, D. Nakashima, T. Nakamura, S. Shcheka, and H. Keppler, “Carbon solubility in silicate melts in equilibrium with a CO- CO2 gas phase and graphite,” Geochimica et Cosmochimica Acta, vol. 259, pp. 129–143, 2019.
- Y. Zhang, C. Yuan, M. Sun et al., “Molybdenum and boron isotopic evidence for carbon-recycling via carbonate dissolution in subduction zones,” Geochimica et Cosmochimica Acta, vol. 278, pp. 340–352, 2020.
- R. Dasgupta and M. M. Hirschmann, “The deep carbon cycle and melting in Earth’s interior,” Earth and Planetary Science Letters, vol. 298, no. 1-2, pp. 1–13, 2010.
- R. Dasgupta, A. Mallik, K. Tsuno, A. C. Withers, G. Hirth, and M. M. Hirschmann, “Carbon-dioxide-rich silicate melt in the Earth’s upper mantle,” Nature, vol. 493, no. 7431, pp. 211–215, 2013.
- F. V. Kaminsky and R. Wirth, “Iron carbide inclusions in lower-mantle diamond from Juina, Brazil,” The Canadian Mineralogist, vol. 49, no. 2, pp. 555–572, 2011.
- V. Stagno, “Carbon, carbides, carbonates and carbonatitic melts in the Earth’s interior,” Journal of the Geological Society, vol. 176, no. 2, pp. 375–387, 2019.
- I. Martinez, E. M. C. Peréz, J. Matas, P. Gillet, and G. Vidal, “Experimental investigation of silicate-carbonate system at high pressure and high temperature,” Journal of Geophysical Research: Solid Earth, vol. 103, no. B3, pp. 5143–5163, 1998.
- M. Isshiki, T. Irifune, K. Hirose et al., “Stability of magnesite and its high-pressure form in the lowermost mantle,” Nature, vol. 427, no. 6969, pp. 60–63, 2004.
- J. F. Lin, J. Liu, C. Jacobs, and V. B. Prakapenka, “Vibrational and elastic properties of ferromagnesite across the electronic spin-pairing transition of iron,” American Mineralogist, vol. 97, no. 4, pp. 583–591, 2012.
- V. Stagno, Y. Tange, N. Miyajima, C. A. McCammon, T. Irifune, and D. J. Frost, “The stability of magnesite in the transition zone and the lower mantle as function of oxygen fugacity,” Geophysical Research Letters, vol. 38, no. 19, article L19309, 2011.
- T. Plank and C. E. Manning, “Subducting carbon,” Nature, vol. 574, no. 7778, pp. 343–352, 2019.
- K. D. Litasov, A. F. Goncharov, and R. J. Hemley, “Crossover from melting to dissociation of CO2 under pressure: implications for the lower mantle,” Earth and Planetary Science Letters, vol. 309, no. 3-4, pp. 318–323, 2011.
- S. G. Charati and S. A. Stern, “Diffusion of gases in silicone polymers molecular dynamics simulations,” Macromolecules, vol. 31, no. 16, pp. 5529–5535, 1998.
- E. Chiavazzo, M. Fasano, P. Asinari, and P. Decuzzi, “Scaling behaviour for the water transport in nanoconfined geometries,” Nature Communications, vol. 5, no. 1, p. 3565, 2014.
- I. Hussain and J. O. Titiloye, “Molecular dynamics simulations of the adsorption and diffusion behavior of pure and mixed alkanes in silicalite,” Microporous and Mesoporous Materials, vol. 85, no. 1-2, pp. 143–156, 2005.
- J. Karger, D. M. Ruthven, and D. N. Theodorou, Diffusion in Nanoporous Materials, Wiley-VCH, Weinheim, Germany, 2012.
- M. Meunier, “Diffusion coefficients of small gas molecules in amorphous cis-1, 4-polybutadiene estimated by molecular dynamics simulations,” The Journal of Chemical Physics, vol. 123, no. 13, p. 134906, 2005.
- J. M. Stubbs, “Molecular simulations of supercritical fluid systems,” The Journal of Supercritical Fluids, vol. 108, pp. 104–122, 2016.
- J. Wagner, O. Adjaoud, K. Marquardt, and S. Jahn, “Anisotropy of self-diffusion in forsterite grain boundaries derived from molecular dynamics simulations,” Contribition to Mineralogy and Petrology, vol. 171, no. 12, p. 98, 2016.
- A. S. Mohammad, J. J. Biernacki, S. Northrup, and M. Adenson, “Diffusion of CO2 and fractional free volume in crystalline and amorphous cellulose,” Journal of Analytical and Applied Pyrolysis, vol. 134, pp. 43–51, 2018.
- G. Milano, G. Guerra, and F. Müller-Plathe, “Anisotropic diffusion of small penetrants in the δ crystalline phase of syndiotactic polystyrene: a molecular dynamics simulation study,” Chemistry of Materials, vol. 14, no. 7, pp. 2977–2982, 2002.
- J.-H. Yoo, A. Breitholz, Y. Iwai, and K.-P. Yoo, “Diffusion coefficients of supercritical carbon dioxide and its mixtures using molecular dynamic simulations,” Korean Journal of Chemical Engineering, vol. 29, no. 7, pp. 935–940, 2012.
- X. Zhao, H. Jin, Y. Chen, and Z. Ge, “Numerical study of H2, CH4, CO, O2 and CO2 diffusion in water near the critical point with molecular dynamics simulation,” Computers & Mathematics with Applications, vol. 81, pp. 759–771, 2021.
- A. Botan, B. Rotenberg, V. Marry, P. Turq, and B. Noetinger, “Carbon dioxide in montmorillonite clay hydrates: thermodynamics structure and transport from molecular simulation,” The Journal of Physical Chemistry C, vol. 114, no. 35, pp. 14962–14969, 2010.
- K. D. Oh, H. Morikawa, S. Iwai, and H. Aoki, “The crystal structure of magnesite,” American Mineralogist, vol. 58, pp. 1029–1033, 1973.
- D. Frenkel and B. Smit, Understanding Molecular Simulation: From Algorithms to Applications, Academic Press, 2002.
- H. Sun, “COMPASS: an ab Initio forcefield optimized for condensed-phase applications - overview with details on alkane and benzene compounds,” The Journal of Physical Chemistry B, vol. 102, no. 38, pp. 7338–7364, 1998.
- L. Yang, L. Liu, H. Liu et al., “Adsorption behavior of CO2 in magnesite micro-pores at high temperature and pressure,” Geoscience Frontiers, vol. 12, no. 2, pp. 991–999, 2021.
- Z. Jin and A. Firoozabadi, “Methane and carbon dioxide adsorption in clay-like slit pores by Monte Carlo simulations,” Fluid Phase Equilibria, vol. 360, pp. 456–465, 2013.
- A. Kadoura, A. K. N. Nair, and S. Sun, “Adsorption of carbon dioxide, methane, and their mixture by montmorillonite in the presence of water,” Microporous and Mesoporous Materials, vol. 225, pp. 331–341, 2016.
- M. Krishnan, M. Saharay, and R. J. Kirkpatrick, “Molecular dynamics modeling of CO2 and poly (ethylene glycol) in montmorillonite: the structure of clay–polymer composites and the incorporation of CO2,” The Journal of Physical Chemistry C, vol. 117, no. 40, pp. 20592–20609, 2013.
- N. Yang, S. Liu, and X. Yang, “Molecular simulation of preferential adsorption of CO2 over CH4 in Na-montmorillonite clay material,” Applied Surface Science, vol. 356, pp. 1262–1271, 2015.
- C. H. Yu, S. Q. Newton, M. A. Norman, L. Schäfer, and D. M. Miller, “Molecular dynamics simulations of adsorption of organic compounds at the clay mineral/aqueous solution Interface,” Structural Chemistry, vol. 14, no. 2, pp. 175–185, 2003.
- K. Aoki, H. Yamawaki, and M. Sakashita, “Phase study of solid CO2 to 20 GPa by infrared-absorption spectroscopy,” Physical Review B, vol. 48, no. 13, pp. 9231–9234, 1993.
- S. A. Bonev, F. Gygi, T. Ogitsu, and G. Galli, “High-pressure molecular phases of solid carbon dioxide,” Physcial Review Letters, vol. 91, no. 6, article 065501, 2003.
- F. A. Gorelli, V. M. Giordano, P. R. Salvi, and R. Bini, “Linear carbon dioxide in the high-pressure high-temperature crystalline phase IV,” Physical Review Letters, vol. 93, no. 20, 2004.
- C. S. Yoo, V. Iota, and H. Cynn, “Nonlinear carbon dioxide at high pressures and temperatures,” Physical Review Letters, vol. 86, no. 3, pp. 444–447, 2001.
- R. Dasgupta, M. M. Hirschmann, and N. Dellas, “The effect of bulk composition on the solidus of carbonated eclogite from partial melting experiments at 3 GPa,” Contributions to Mineralogy and Petrology, vol. 149, no. 3, pp. 288–305, 2005.
- K. D. Litasov, A. Shatskiy, E. Ohtani, and G. M. Yaxley, “Solidus of alkaline carbonatite in the deep mantle,” Geology, vol. 41, no. 1, pp. 79–82, 2013.
- M. Le Voyer, K. A. Kelley, E. Cottrell, and E. H. Hauri, “Heterogeneity in mantle carbon content from CO2-undersaturated basalts,” Nature Communications, vol. 8, no. 1, article 14062, 2017.
- A. Rohrbach and M. W. Schmidt, “Redox freezing and melting in the Earth’s deep mantle resulting from carbon- iron redox coupling,” Nature, vol. 472, no. 7342, pp. 209–212, 2011.
- P. N. Perera, H. Deng, P. J. Schuck, and B. Gilbert, “Diffusivity of carbon dioxide in aqueous solutions under geologic carbon sequestration conditions,” Journal of Physical Chemistry B, vol. 122, no. 16, pp. 4566–4572, 2018.
- R. Dasgupta and M. M. Hirschmann, “Melting in the Earth’s deep upper mantle caused by carbon dioxide,” Nature, vol. 440, no. 7084, pp. 659–662, 2006.
- N. H. Sleep and K. Zahnle, “Carbon dioxide cycling and implications for climate on ancient Earth,” Journal of Geophysical Research: Planets, vol. 106, no. E1, pp. 1373–1399, 2001.
- J. Huang, S.-G. Li, Y. Xiao, S. Ke, W.-Y. Li, and Y. Tian, “Origin of low δ26Mg Cenozoic basalts from South China Block and their geodynamic implications,” Geochimica et Cosmochimica Acta, vol. 164, pp. 298–317, 2015.
- P. B. Kelemen and C. E. Manning, “Reevaluating carbon fluxes in subduction zones, what goes down, mostly comes up,” Proceedings of the National Academy of Sciences of the United States of America, vol. 112, no. 30, pp. E3997–E4006, 2015.
- S.-G. Li, W. Yang, S. Ke et al., “Deep carbon cycles constrained by a large-scale mantle Mg isotope anomaly in eastern China,” National Science Review, vol. 4, no. 1, pp. 111–120, 2017.
- H. C. Tian, W. Yang, S. G. Li, S. Ke, and Z. Y. Chu, “Origin of low δ26Mg basalts with EM-I component: evidence for interaction between enriched lithosphere and carbonated asthenosphere,” Geochimica et Cosmochimica Acta, vol. 188, pp. 93–105, 2016.
- W. Yang, F. Z. Teng, H. F. Zhang, and S. G. Li, “Magnesium isotopic systematics of continental basalts from the north China craton: implications for tracing subducted carbonate in the mantle,” Chemical Geology, vol. 328, pp. 185–194, 2012.
Copyright © 2021 Lei Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.