Generalization of the Unified Analytic Melt-Shear Model to Multi-Phase Materials: Molybdenum as an Example

The unified analytic melt-shear model that we introduced a decade ago is generalized to multi-phase materials. A new scheme for calculating the values of the model parameters for both the cold ( T = 0 ) shear modulus ( G ) and the melting temperature at all densities ( ρ ) is developed. The generalized melt-shear model is applied to molybdenum, a multi-phase material with a body-centered cubic (bcc) structure at low ρ which loses its dynamical stability with increasing pressure (P) and is therefore replaced by another (dynamically stable) solid structure at high ρ . One of the candidates for the high- ρ structure of Mo is face-centered cubic (fcc). The model is compared to (i) our ab initio results on the cold shear modulus of both bcc-Mo and fcc-Mo as a function of ρ , and (ii) the available theoretical results on the melting of bcc-Mo and our own quantum molecular dynamics (QMD) simulations of one melting point of fcc-Mo. Our generalized model of G ( ρ , T ) is used to calculate the shear modulus of bcc-Mo along its principal Hugoniot. It predicts that G of bcc-Mo increases with P up to ∼240 GPa and then decreases at higher P. This behavior is intrinsic to bcc-Mo and does not require the introduction of another solid phase such as Phase II suggested by Errandonea et al. Generalized melt-shear models for Ta and W also predict an increase in G followed by a decrease along the principal Hugoniot, hence this behavior may be typical for transition metals with ambient bcc structure that dynamically destabilize at high P. Thus, we concur with the conclusion reached in several recent papers (Nguyen et al., Zhang et al., Wang et al.) that no solid-solid phase transition can be definitively inferred on the basis of sound velocity data from shock experiments on Mo. Finally, our QMD simulations support the validity of the phase diagram of Mo suggested by Zeng et al.

Experimental constraints on light elements in the Earth's outer core

Earth's outer core is liquid and dominantly composed of iron and nickel (~5-10 wt%). Its density, however, is ~8% lower than that of liquid iron, and requires the presence of a significant amount of light element(s). A good way to specify the light element(s) is a direct comparison of density and sound velocity measurements between seismological data and those of possible candidate compositions at the core conditions. We report the sound velocity measurements of a model core composition in the Fe-Ni-Si system at the outer core conditions by shock-wave experiments. Combining with the previous studies, we found that the best estimate for the outer core's light elements is ~6 wt% Si, ~2 wt% S, and possible ~1-2.5 wt% O. This composition satisfies the requirements imposed by seismology, geochemistry, and some models of the early core formation. This finding may help us to further constrain the thermal structure of the Earth and the models of Earth's core formation.

Practical uncertainty reduction and quantification in shock physics measurements

We report the development of a simple error analysis sampling method for identifying intersections and inflection points to reduce total uncertainty in experimental data. This technique was used to reduce uncertainties in sound speed measurements by 80% over conventional methods. Here, we focused on its impact on a previously published set of Mo sound speed data and possible implications for phase transition and geophysical studies. However, this technique's application can be extended to a wide range of experimental data.

Hidden carbon in Earth's inner core revealed by shear softening in dense Fe7C3

Earth's inner core is known to consist of crystalline iron alloyed with a small amount of nickel and lighter elements, but the shear wave (S wave) travels through the inner core at about half the speed expected for most iron-rich alloys under relevant pressures. The anomalously low S-wave velocity (vS) has been attributed to the presence of liquid, hence questioning the solidity of the inner core. Here we report new experimental data up to core pressures on iron carbide Fe7C3, a candidate component of the inner core, showing that its sound velocities dropped significantly near the end of a pressure-induced spin-pairing transition, which took place gradually between 10 GPa and 53 GPa. Following the transition, the sound velocities increased with density at an exceptionally low rate. Extrapolating the data to the inner core pressure and accounting for the temperature effect, we found that low-spin Fe7C3 can reproduce the observed vS of the inner core, thus eliminating the need to invoke partial melting or a postulated large temperature effect. The model of a carbon-rich inner core may be consistent with existing constraints on the Earth's carbon budget and would imply that as much as two thirds of the planet's carbon is hidden in its center sphere.

Some recent advances in understanding the mineralogy of Earth's deep mantle

Understanding planetary structure and evolution requires a detailed knowledge of the properties of geological materials under the conditions of deep planetary interiors. Experiments under the extreme pressure-temperature conditions of the deep mantle are challenging, and many fundamental properties remain poorly constrained or are inferred only through uncertain extrapolations from lower pressure-temperature states. Nevertheless, the last several years have witnessed a number of new developments in this area, and a broad overview of the current understanding of the Earth's lower mantle is presented here. Some recent experimental and theoretical advances related to the lowermost mantle are highlighted. Measurements of the equation of state and deformation behaviour of (Mg,Fe)SiO3 in the CaIrO3-type (post-perovskite) structure yield insights into the nature of the core-mantle boundary region. Theoretical studies of the behaviour of MgSiO3 liquids under high pressure-temperature conditions provide constraints on melt volumes, diffusivities and viscosities that are relevant to understanding both the early Earth (e.g. deep magma oceans) and seismic structure observed in the present Earth (e.g. ultra-low-velocity zones).

The physical and chemical composition of the lower mantle

This article reviews some of the recent advances made within the field of mineral physics. In order to link the observed seismic and density structures of the lower mantle with a particular mineral composition, knowledge of the thermodynamic properties of the candidate materials is required. Determining which compositional model best matches the observed data is difficult because of the wide variety of possible mineral structures and compositions. State-of-the-art experimental and analytical techniques have pushed forward our knowledge of mineral physics, yet certain properties, such as the elastic properties of lower mantle minerals at high pressures and temperatures, are difficult to determine experimentally and remain elusive. Fortunately, computational techniques are now sufficiently advanced to enable the prediction of these properties in a self-consistent manner, but more results are required.A fundamental question is whether or not the upper and lower mantles are mixing. Traditional models that involve chemically separate upper and lower mantles cannot yet be ruled out despite recent conflicting seismological evidence showing that subducting slabs penetrate deep into the lower mantle and that chemically distinct layers are, therefore, unlikely.Recent seismic tomography studies giving three-dimensional models of the seismic wave velocities in the Earth also base their interpretations on the thermodynamic properties of minerals. These studies reveal heterogeneous velocity and density anomalies in the lower mantle, which are difficult to reconcile with mineral physics data.