Abstract:We have performed a series of molecular dynamics simulations aimed at the evaluation of the solubility of CO2 in silicate melts of natural composition (from felsic to ultramafic). In making in contact within the simulation cell a supercritical CO2 phase with a silicate melt of a given composition, we have been able to evaluate (i) the solubility of CO2 in the P–T range 1473–2273K and 20–150kbar, (ii) the density change experienced by the CO2-bearing melt, (iii) the respective concentrations of CO2 and CO 3 2 - species in the melt, (iv) the lifetime and the diffusivity of these species and (v) the structure of the melt around the carbonate groups. The main results are the following: (1) The solubility of CO2 increases markedly with the pressure in the three investigated melts (a rhyolite, a mid-ocean ridge basalt and a kimberlite) from about ∼2wt% CO2 at 20kbar to ∼25wt% at 100kbar and 2273K. The solubility is found to be weakly dependent on the melt composition (as far as the present compositions are concerned) and it is only at very high pressure (above ∼100kbar) that a clear hierarchy between solubilities occurs (rhyolite<MORB<kimberlite). Furthermore at a given pressure the calculated solubility is negatively correlated with the temperature. (2) In CO2-saturated melts, the proportion of carbonate ions ( CO 3 2 - ) is positively correlated with the pressure at isothermal condition and is negatively correlated with the temperature at isobaric condition (and vice versa for molecular CO2). Furthermore, at fixed (P, T) conditions the proportion of carbonate ions is higher in CO2-undersaturated melts than in the CO2-saturated melt. Although the proportion of molecular CO2 decreases when the degree of depolymerization of the melt increases, it is still significant in CO2-saturated basic and ultrabasic compositions at high temperatures. This finding is at variance with experimental data on CO2-bearing glasses which show no evidence of molecular CO2 as soon as the degree of depolymerization of the melt is high (e.g. basalt). These conflicting results can be reconciled with each other by noticing that a simple low temperature extrapolation of the simulation data predicts that the proportion of molecular CO2 in basaltic melts might be negligible in the glass at room temperature. (3) The carbonate ions are found to be transient species in the liquid phase, with a lifetime increasing exponentially with the inverse of the temperature. Contrarily to a usual assumption, the diffusivity of carbonate ions into the liquid silicate is not vanishingly small with respect to that of CO2 molecules: in MORB they differ from each other by a factor of ∼6 at 1473K and only a factor of ∼2 at 2273K. Although the bulk diffusivity of CO2 is governed primarily by the diffusivity of CO2 molecules, the carbonate ions contribute significantly to the diffusivity of CO2 in depolymerized melts. (4) Concerning the structure of the CO2-bearing silicate melt, the carbonate ions are found to be preferentially associated with NBO’s of the melt, with an affinity for NBOs which exceeds that for BOs by almost one order of magnitude. This result explains why the concentration in carbonate ions is positively correlated with the degree of depolymerization of the melt and diminishes drastically in fully polymerized melts where the number of NBO’s is close to zero. Furthermore, the network modifier cations are not randomly distributed in the close vicinity of carbonate groups but exhibit a preferential ordering which depends at once on the nature of the cation and on the melt composition. However at the high temperatures investigated here, there is no evidence of long lived complexes between carbonate groups and metal cations.

Effect of pressure on the short-range structure and speciation of carbon in alkali silicate and aluminosilicate glasses and melts at high pressure up to 8 GPa: 13 C, 27 Al, 17 O and 29 Si solid-state NMR study

The Effect of Na/Si on the Structure of Sodium Silicate and Aluminosilicate Glasses Quenched from Melts at High Pressure : A Multi-Nuclear (Al-27, Na-23, O-17) 1D and 2D Solid-State NMR Study

Effect of Network Polymerization on the Pressure-Induced Structural Changes in Sodium Aluminosilicate Glasses and Melts: 27al and 17O Solid-State NMR Study

Effect of Network Polymerization on the Pressure-Induced Structural Changes in Sodium Aluminosilicate Glasses and Melts: 27 Al and 17 O Solid-State NMR Study

Nature of polymerization and properties of silicate melts and glasses at high pressure

Effect of Al<sub>2</sub>O<sub>3</sub> on Sound Velocity of MgSiO<sub>3</sub> Glass at High Pressure

Oxygen-17 nuclear magnetic resonance study of the structure of mixed cation calcium-sodium silicate glasses at high pressure: implications for molecular link to element partitioning between silicate liquids and crystals.

Phase Relations in the System MgSiO3‐Al2O3 Up to 2300 K at Lower Mantle Pressures

Probing Carbon-Bearing Species and CO2inclusions in Amorphous Carbon-Mgsio3enstatite Reaction Products at 1.5 GPa: Insights From13c High-Resolution Solid-State NMR

Transformations to the aluminum coordination environment and network polymerization in amorphous aluminosilicates under pressure

Development of chemical and topological structure in aluminosilicate liquids and glasses at high pressure

Densification of Sodium and Magnesium Aluminosilicate Glasses at Ambient Temperature: Structural Investigations by Solid-state Nuclear Magnetic Resonance and Molecular Dynamics Simulations

Order and Disorder in Sodium Silicate Glasses and Melts at 10 Gpa

Pressure dependent structure of amorphous magnesium aluminosilicates: The effect of replacing magnesia by alumina at the enstatite composition

Carbon dioxide in silicate melts: A molecular dynamics simulation study

First-principles molecular dynamics simulations of anorthite (CaAl2Si2O8) glass at high pressure

X-ray Raman Scattering Study of MgSiO3 Glass at High Pressure: Implication for Triclustered MgSiO3 Melt in Earth's Mantle

Coordination Changes in Densified Aluminate Glass upon Compression up to 65 GPa: A View from Solid-State Nuclear Magnetic Resonance

Quantitative studies on the microstructure of ternary CaO–Al2O3–SiO2 glasses by Raman spectroscopy, 27Al MAS NMR and quantum chemistry ab initio calculation

Structure of Alkali Borate Glasses at High Pressure: B and Li K-edge Inelastic X-ray Scattering Study.

Raman Spectroscopy of SiO2 Glass at High Pressure.