Equation of state of molten fayalite (Fe2SiO4)

We have conducted new equation of state measurements on liquid fayalite (Fe2SiO4) in a collaborative, multi-technique study. Using a shared bulk starting material, we have measured the liquid density, the bulk modulus (K), and its pressure derivative (K’) from 1 atm to 163 GPa using 1-atm double-bob Archimedean and ultrasonic, sink/float, and shock wave techniques to form a coherent, internally consistent equation of state. Previous shock studies of liquid fayalite were conducted up to pressures of 40 GPa1; we extended this data set with two additional pre-heated, molten (1573 K) fayalite shock compression experiments at 121 and 163 GPa. Linear fitting of this data in shock velocity (US)-particle velocity (up) space defines a Hugoniot with an unconstrained zero-pressure intercept that crosses within error at the bulk sound speed (Co) determined by ultrasonic techniques. Fixing the intercept at this ultrasonic value reduces the error on the linear fit and yields the relation: US =1.65(0.02)up+ 2.4377(0.006) km/s. This relationship indicates that the behavior of the liquid is relaxed during shock compression and demonstrates consistency across experimental methods. Likewise, results from new static compression sink/float experiments conducted in piston-cylinder and multi-anvil devices are in agreement with shock wave and ultrasonic data, consistent with an isothermal K=19.4 and K’=5.57 at 1500°C. In solid materials, the Grüneisen parameter (γ) generally decreases upon compression. However, preliminary calculations for γ of this liquid using additional initially solid shock data from Chen et al.(2002) indicate that γ increases upon compression. Using the functional form γ = γo(ρo/ρ)q at a density of 7.65 Mg/m3 yields a q value of -1.77 (γo = 0.41 is known from low-pressure data), which is similar to the reported q values of forsterite2, enstatite3, and anorthite-diopside liquids4. This result shows that iron-bearing mafic to ultramafic silicate liquids follow the same general behavior as iron-free liquids such that -2.0 ≤ q ≤ -1.5 for the compression range 1 ≥ ρo/ρ ≥ 0.50. We will be performing an additional shock wave experiment on initially solid (300 K) fayalite to confirm this result. We will be continuing collaborative equation of state measurements on additional iron-bearing silicate liquids, working to further clarify the properties of melts and their importance to understanding the dynamics of the early magma ocean and of melt migration within the mantle. In particular, understanding the properties of iron-rich silicates and their melts will constrain hypotheses of melting and of iron enrichment for explaining the occurrence and characteristics of ultra-low velocity zones near the CMB.