Carbon dioxide in silicate melts: A molecular dynamics simulation study
  The distribution, recycling and storage of carbon in the Earth are of fundamental importance to understand the global carbon cycle between the deep Earth and near surface reservoirs. Degassing of CO2 at mid-ocean ridges may give information on the source region but the very low solubility of CO2 in tholeitic basalts has for consequence that near all Mid-Ocean Ridge Basalts glasses exsolve their CO2 rich vapor at shallow depth as they approach the ocean floor. Hence their CO2 contents mostly represent the pressure at eruption and not the source region. Recent petrological investigations have shown that the presence of carbonates at depth in the upper mantle has a large effect on the solidus of carbonated silicates by inducing incipient melting at much lower temperature. So the role of carbon-rich melts at great depth is now becoming a credible scenario to explain the extraction of CO2 from the source region to the surface. During the last three decades many studies have been devoted to measure the solubility of CO2 in silicate melts of various composition. But due to experimental difficulties these studies were generally restricted to low and moderate pressures (below ~20 kbar). By performing a series of molecular dynamics simulation where a supercritical CO2 phase is in contact with a silicate melt of various composition (from felsic to ultrabasic) at different temperatures (1473-2273K) and pressures (20-150kbar), we have been able to evaluate the solubility of CO2, the population of molecular and carbonate species, their diffusivity through the melt and the local structure. We show that this kind of molecular simulation is a useful theoretical guide to better understand the behavior of CO2 in magmas at depth. 
B. Guillot and N. Sator
cond-mat.stat-mech physics.chem-ph