Buoyancy and Structure of Volatile-Rich Silicate Melts

The early Earth was marked by at least one global magma ocean. Melt buoyancy played a major role for its evolution. Here we model the composition of the magma ocean using a six-component pyrolite melt, to which we add volatiles in the form of carbon as molecular CO or CO₂ and hydrogen as molecular H₂O or through substitution for magnesium. We compute the density relations from first-principles molecular dynamics simulations. We find that the addition of volatiles renders all the melts more buoyant compared to the reference volatile-free pyrolite. The effect is pressure dependent, largest at the surface, decreasing to about 20 GPa, and remaining roughly constant to 135 GPa. The increased buoyancy would have enhanced convection and turbulence, and thus promoted the chemical exchanges of the magma ocean with the early atmosphere. We determine the partial molar volume of both H₂O and CO₂ throughout the magma ocean conditions. We find a pronounced dependence with temperature at low pressures, whereas at megabar pressures the partial molar volumes are independent of temperature. At all pressures, the polymerization of the silicate melt is strongly affected by the amount of oxygen added to the system while being only weakly affected by the specific type of volatile added.