The Combined Strength of Thermodynamics and Comparative Planetology: Application of Activity Models to Core Formation in Terrestrial Bodies
Recent models for accretion of terrestrial bodies involve metal-silicate equilibrium as the metallic core formed during growth. Most elements considered are either refractory or well studied elements for which effects of pressure, temperature, oxygen fugacity, and metallic liquid composition are well known. There are a large number of elements that are both siderophile and volatile, whose fate in such models is unknown, largely due to a lack of data at comparable conditions and com-positions (FeNi core with light elements such as S, C, Si, and O). We have focused on Ge, In, As, Sb and determined the effect of Si and C on metal-silicate partitioning, and developed a thermo-dynamic model that allows application of these new data to a wide range of planetary bodies. New experiments: We have previously carried out experiments with FeSi metallic liquid at C-saturated conditions at 1600 and 1800 C [4]. In a new series of experiments we investigate the effect of Si in carbon-free systems at 1600 C for comparison. Experiments were carried out at 1 GPa in MgO capsules using the same basaltic starting composition as in previous studies. The MgO capsule reacts with the silicate melt to form more MgO-rich liquids that have 22-26 wt% MgO. Experimental met-als and silicates were analyzed using a combination of electron microprobe analysis and laser ablation ICP-MS. Results: The new results can be interpreted by considering Ge as an example, in the simple exchange equilibrium Fe + GeO = FeO + Ge, where the equilibrium constant Kd can be examined as a function of Si content of the metal. The slope of lnKd vs. (1-XSi) for this new series allows derivation of the epsilon interaction parameter for each of these four elements and Si (both C-saturated and C-free).All four elements have positive epsilon values, indicating that Si causes a decrease in the partition coefficients; values are 6.6, 6.5, 27.8 and 25.2 for In, Ge, As, and Sb, respectively, at 1 GPa and 1600 C. As an example of how large the effect of Si can be, these epsilon values correspond to activity coefficients (gamma) for As of 0.01 when XSi = 0, and up to gamma = 23 when XSi = 0.2. Combining these new results with previous determinations [5,6] of epsilon parameters for S and C for these elements allows us calculate activity of Ge, In, As, and Sb in Fe-Ni-Si-S-C-O metallic liquids. We apply this new model to sever-al terrestrial bodies such as Earth (Si-rich core), Mars (S-rich core), Moon (S-, C-, and Si-poor core), and Vesta, and examine the resulting core and mantle concentrations of these elements. Mantle concentrations of these four elements are well explained for Earth and Mars in models that call for mid-mantle equilibration between Si-bearing and S-bearing FeNi cores, respectively. Modeling results for the Moon and Vesta will also be presented.