Physical State of an Early Magma Ocean Constrained by the Thermodynamics and Viscosity of Iron Silicate Liquid

Yicheng Sun,Huiqun Zhou,Xiandong Liu,Kun Yin,Xiancai Lu
DOI: https://doi.org/10.1016/j.epsl.2020.116556
IF: 5.3
2020-01-01
Earth and Planetary Science Letters
Abstract:The initial cooling stage and crystallization of a magma ocean were significantly affected by the thermodynamics and viscosity of iron silicate liquid. Here, we use first-principles molecular dynamics based on density functional theory plus the Hubbard U method to study the structures, thermodynamics and viscosities of Fe2SiO4 liquid at 3000-6000 K under pressure conditions spanning the entire mantle. Our calculations show that the compressibility of Fe2SiO4 liquid is weakened at pressures above 50 GPa. The densification of Fe2SiO4 liquid under compression is mainly manifested by a decrease in the distance among Si-O polyhedra and an increase in the cation-anion coordination number. We find that the viscosity of the Fe2SiO4 liquid has a positive pressure dependence, and the effect of pressure on viscosity decreases with increasing temperature. The predicted behaviour of viscosity coefficients of the Fe2SiO4 liquid at different P-T conditions can be described well as: eta(P, T) = exp[-7.89 + 0.0042P - 0.0000632P(2) + (3800 + 61.6P + 0.18P(2))/(T - 1000)]. Compared to Mg2SiO4 liquid, the viscosity of Fe2SiO4 liquid is lower under the pressure of the lower mantle along isotherms. We calculate the adiabat of the Fe2SiO4 liquid in a magma ocean, and the results show that the iron content has little effect on the adiabats of (Mg,Fe)(2)SiO4 liquid in the magma ocean. Combining our calculated thermodynamic properties and viscosities with those of previous studies, the physical properties of an early magma ocean were constrained. The upper bound of the potential temperature for magma ocean crystallization is 3250 K. We constructed the range of the viscosity profile of a magma ocean and found that the lower viscosity of iron-rich silicate liquid would maintain the viscosity of an early magma ocean in the range of several mPas. As the magma ocean starts to crystallize, the lower bound of the viscosity at the surface and the core-mantle boundary are 0.0020 Pas and 0.0068 Pas, respectively. The weakening of the viscosity of the iron silicate liquids would have a significant impact on the evolution of the magma ocean. (C) 2020 Elsevier B.V. All rights reserved.
What problem does this paper attempt to address?