Formation, preservation and extinction of high-pressure minerals in meteorites: temperature effects in shock metamorphism and shock classification

The goal of classifying shock metamorphic features in meteorites is to estimate the corresponding shock pressure conditions. However, the temperature variability of shock metamorphism is equally important and can result in a diverse and heterogeneous set of shock features in samples with a common overall shock pressure. In particular, high-pressure (HP) minerals, which were previously used as a solid indicator of high shock pressure in meteorites, require complex pressure-temperature-time (P-T-t) histories to form and survive. First, parts of the sample must be heated to melting temperatures, at high pressure, to enable rapid formation of HP minerals before pressure release. Second, the HP minerals must be rapidly cooled to below a critical temperature, before the pressure returns to ambient conditions, to avoid retrograde transformation to their low-pressure polymorphs. These two constraints require the sample to contain large temperature heterogeneities, e.g. melt veins in a cooler groundmass, during shock. In this study, we calculated shock temperatures and possible P-T paths of chondritic and differentiated mafic-ultramafic rocks for various shock pressures. These P-T conditions and paths, combined with observations from shocked meteorites, are used to constrain shock conditions and P-T-t histories of HP-mineral bearing samples. The need for rapid thermal quench of HP phases requires a relatively low bulk-shock temperature and therefore moderate shock pressures below ~ 30 GPa, which matches the stabilities of these HP minerals. The low-temperature moderate-pressure host rock generally shows moderate shock-deformation features consistent with S4 and, less commonly, S5 shock stages. Shock pressures in excess of 50 GPa in meteorites result in melt breccias with high overall post-shock temperatures that anneal out HP-mineral signatures. The presence of ringwoodite, which is commonly considered an indicator of the S6 shock stage, is inconsistent with pressures in excess of 30 GPa and does not represent shock conditions different from S4 shock conditions. Indeed, ringwoodite and coexisting HP minerals should be considered as robust evidence for moderate shock pressures (S4) rather than extreme shock (S6) near whole-rock melting.

Seismically determined elastic parameters for Earth's outer core

Turbulent convection of the liquid iron alloy outer core generates Earth's magnetic field and supplies heat to the mantle. The exact composition of the iron alloy is fundamentally linked to the processes powering the convection and can be constrained by its seismic properties. Discrepancies between seismic models determined using body waves and normal modes show that these properties are not yet fully agreed upon. In addition, technical challenges in experimentally measuring the equation-of-state (EoS) parameters of liquid iron alloys at high pressures and temperatures further complicate compositional inferences. We directly infer EoS parameters describing Earth's outer core from normal mode center frequency observations and present the resulting Elastic Parameters of the Outer Core (EPOC) seismic model. Unlike alternative seismic models, ours requires only three parameters and guarantees physically realistic behavior with increasing pressure for a well-mixed homogeneous material along an isentrope, consistent with the outer core's condition. We show that EPOC predicts available normal mode frequencies better than the Preliminary Reference Earth Model (PREM) while also being more consistent with body wave-derived models, eliminating a long-standing discrepancy. The velocity at the top of the outer core is lower, and increases with depth more steeply, in EPOC than in PREM, while the density in EPOC is higher than that in PREM across the outer core. The steeper profiles and higher density imply that the outer core comprises a lighter but more compressible alloy than that inferred for PREM. Furthermore, EPOC's steeper velocity gradient explains differential SmKS body wave travel times better than previous one-dimensional global models, without requiring an anomalously slow ~90- to 450-km-thick layer at the top of the outer core.

Continuous Sound Velocity Measurements along the Shock Hugoniot Curve of Quartz

We report continuous measurements of the sound velocity along the principal Hugoniot curve of α quartz between 0.25 and 1.45 TPa, as determined from lateral release waves intersecting the shock front as a function of time in decaying-shock experiments. The measured sound velocities are lower than predicted by prior models, based on the properties of stishovite at densities below ∼7  g/cm^{3}, but agree with density functional theory molecular dynamics calculations and an empirical wide-regime equation of state presented here. The Grüneisen parameter calculated from the sound velocity decreases from γ∼1.3 at 0.25 TPa to 0.66 at 1.45 TPa. In combination with evidence for increased (configurational) specific heat and decreased bulk modulus, the values of γ suggest a high thermal expansion coefficient at ∼0.25-0.65  TPa, where SiO_{2} is thought to be a bonded liquid. From our measurements, dissociation of the molecular bonds persists to ∼0.65-1.0  TPa, consistent with estimates by other methods. At higher densities, the sound velocity is close to predictions from previous models, and the Grüneisen parameter approaches the ideal gas value.

Three-phase junction for modulating electron-hole migration in anatase-rutile photocatalysts

The heterophase solid-solid junction as an important type of structure unit has wide applications for its special mechanics and electronic properties. Here we present a first three-phase atomic model for the anatase-rutile TiO2 heterophase junction and determine its optical and electronic properties, which leads to resolution of the long-standing puzzles on the enhanced photocatalytic activity of anatase-rutile photocatalysts. By using a set of novel theoretical methods, including crystal phase transition pathway sampling, interfacial strain analysis and first principles thermodynamics evaluation of holes and electrons, we identify an unusual structurally ordered three-phase junction, a layer-by-layer "T-shaped" anatase/TiO2-II/rutile junction, for linking anatase with rutile. The intermediate TiO2-II phase, although predicted to be only a few atomic layers thick in contact with anatase, is critical to alleviate the interfacial strain and to modulate photoactivity. We demonstrate that the three-phase junction acts as a single-way valve allowing the photogenerated hole transfer from anatase to rutile but frustrating the photoelectron flow in the opposite direction, which otherwise cannot be achieved by an anatase-rutile direct junction. This new model clarifies the roles of anatase, rutile and the phase junction in achieving high photoactivity synergistically and provides the theoretical basis for the design of better photocatalysts by exploiting multi-phase junctions.

Equations of state for solids under strong compression

Various approaches for the representation of equation of state data for solids under strong compression are discussed. The theoretical background for reasonable extrapolations to higher pressures and higher as well as lower temperatures is described. The distinction between ideal, regular, and anomalous behaviour allows to gain deeper insight into the electronic changes occurring in various solids under strong compression. The discussion of experimental data for various regular solids leads finally to an estimate of the accuracy obtained in the present realisation of a practical pressure scale based on equation of state measurements.

Achieving high-density states through shock-wave loading of precompressed samples

Materials can be experimentally characterized to terapascal pressures by sending a laser-induced shock wave through a sample that is precompressed inside a diamond-anvil cell. This combination of static and dynamic compression methods has been experimentally demonstrated and ultimately provides access to the 10- to 100-TPa (0.1-1 Gbar) pressure range that is relevant to planetary science, testing first-principles theories of condensed matter, and experimentally studying a new regime of chemical bonding.

Applicability of isothermal unrealistic two-parameter equations of state for solids

The aim of the present study, an extension of a recent one (Bose Roy and Bose Roy 2005 J. Phys.: Condens. Matter 17 6193), is to assess and compare the curve-fitting utility of the isothermal unrealistic two-parameter equations of state for solids (EOS), proposed at different stages in the development of the EOS field, for the purposes of smoothing and interpolation of pressure-volume data, and extraction of accurate values of the isothermal bulk modulus and its pressure derivative. To this end, 21 such EOSs are considered, formulated by/labelled as Born-Mie (1920), Born-Mayer (1932), Bardeen (1938), Slater-Morse (1939), Birch-Murnaghan (1947), Pack-Evans-James (1948), Lagrangian (1951), Davydov (1956), Davis and Gordon (1967), Onat and Vaisnys (1967), Grover-Getting-Kennedy (1973), Brennan-Stacey (1979), Walzer-Ullmann-Pan'kov (1979), Rydberg (1981), Dodson (1987), Holzapfel (1991), Parsafar-Mason (1994), Shanker-Kushwah-Kumar (1997), Poirier-Tarantola (1998), Deng-Yan (2002) and Kun-Loa-Syassen (2003). Furthermore, all these EOSs are compared with our three-parameter EOS, as well as its two-parameter counterpart proposed in this work. We have applied all the EOS models, with no constraint on the parameters, to the accurate and model-independent isotherms of nine solids. The applicability has been assessed in terms of an unbiased composite test, comprising fitting accuracy, agreement of the fit parameters with experiment, stability of the fit parameters with variation in the compression/pressure ranges and on the basis of the number of wiggles of the data deviation curves about the fit parameters. Furthermore, a rigorous method is devised to scale the relative adequacy of the EOSs with respect to the test parameters. A number of remarkable findings emerge from the present study. Surprisingly, both the old EOSs, the Born-Mie and the Pack-Evans-James, are significantly better in their curve-fitting capability than the Birch-Murnaghan EOS which has been widely used and continues to be used for curve-fitting purposes as a standard EOS in the literature. The Born-Mayer as well as the Walzer-Ullmann-Pan'kov models also fit isotherms better than the Birch. The performance of the EOS based on the Rydberg potential-that has been rediscovered by Rose et al (1984 Phys. Rev. B 29 2963), and strongly promoted by Vinet et al (1989 J. Phys.: Condens. Matter 1 1941) as the so-called universal equation of state, and is currently used as a standard EOS along with that of the Birch-is very poor, on a comparative scale. Furthermore, the curve-fitting capability of our original three-parameter EOS, and more importantly its two-parameter counterpart, is superior to all the isothermal unrealistic two-parameter EOSs so far proposed in the literature.

Density and composition of the lower mantle

The observed density distribution of the lower mantle is compared with density measurements of the (M g,Fe)SiO 3 perovskite and (Mg,Fe)O magnesiowtistite highpressure phases as functions of pressure, tem perature and composition. We find that for plausible bounds on the composition of the upper mantle (ratio of magnesium to iron + magnesium components x M g ^ 0.88) and the temperature in the lower mantle ( T ^ 2000 K), the high-pressure mineral assemblage of upper-mantle composition is at least 2 .6 ( ± 1 ) % less dense than the lower m antle over the depth range 1000-2000 km. Thus, we find that a model of uniform m antle composition is incompatible with the existing mineralogical and geophysical data. Instead, we expect that the mantle is stratified, with the upper and lower m antle convecting separately, and we estimate that the compositional density difference between these regions is about 5 ( + 2) %. The stratification may not be perfect (‘leaky layering’), but significant intermixing and homogenization of the upper and lower m antle over geological timescales are precluded.

Condensed-Matter Energetics from Diatomic Molecular Spectra

Analyses of molecular spectra and compression data from crystals show that a single function successfully describes the dependence on interatomic separation of both the potential energy of diatomic molecules and the cohesive binding energy of condensed matter. The empirical finding that one function describes interatomic energies for such diverse forms of matter and over a wide range of conditions can be used to extend condensed-matter equations of state but warrants further theoretical study.