Carbon dioxide in silicate melts at upper mantle conditions: Insights from atomistic simulations

• Published in: Chemical geology 2015-12, Vol.418, p.77-88
• Main Authors: Vuilleumier, Rodolphe, Seitsonen, Ari P., Sator, Nicolas, Guillot, Bertrand
• Format: Article
• Language: English
• Subjects: Basalt; Carbon dioxide; Carbonates; CO2; Dynamical systems; Earth Sciences; First-principle and classical molecular dynamics simulations; Geochemistry; Kimberlite; Mantle; Melts; Molecular dynamics; Sciences of the Universe; Silicates; Simulation; Speciation
• ISSN: 0009-2541; 1872-6836
• DOI: 10.1016/j.chemgeo.2015.02.027

The detail of the incorporation of carbon dioxide in silicate melts at upper mantle conditions is still badly known. To give some theoretical guidance, we have performed first-principle molecular dynamics simulations (FPMD) to quantify the speciation and the incorporation of carbon dioxide in two CO2-rich silicate melts (~20wt.% CO2 at 2073K and 12GPa), a basaltic and a kimberlitic composition chosen in the CaO–MgO–Al2O3–SiO2 system. In the basaltic composition, carbon dioxide is incorporated under the form of a minority population of CO2 molecules and a prevailing population of carbonate ions (CO32−). In contrast, the amount of CO2 molecules is found to be very small in the kimberlitic melt. Moreover, a new (transient) species has been identified, the pyrocarbonate ion C2O52− issued from the reaction between CO2 and CO32−. With regard to the structure of the CO2-bearing melts, it is shown that the carbonate ions modify the silicate network by transforming some of the oxygen atoms into bridging carbonates, non-bridging carbonates, and free carbonates, with a distribution depending on the melt composition. In the basaltic melt a majority of carbonate ions are non-bridging or free, whereas in the kimberlitic melt, most of the carbonate ions are under the form of free carbonates linked to alkaline earth cations. Surprisingly, the addition of CO2 only has a weak influence on the diffusion coefficients of the elements of the melt. The consequence is that the strong enhancement of the electrical conductivity reported recently for carbonated basalts (Sifré et al., 2014, Nature 509, 81), can be reproduced by simulation only if one assumes that the ionic charges assigned to the elements of the melt depend, in a non-trivial way, on the CO2 content. Finally, a comparison of the FPMD calculations with classical molecular dynamics simulations using an empirical force field of the literature (Guillot and Sator, 2011, GCA 75, 1829) shows that the latter one needs some improvement. •Incorporation of CO2 in basalt and kimberlite is investigated by atomistic simulation.•First-principle and classical molecular dynamics simulations were performed.•The speciation of CO2 and the structure of CO2-bearing melts are evaluated.•CO2 has a weak influence on the diffusion coefficients of the elements in the melt.•The electrical conductivity of carbonated melts is evaluated and discussed.