Deformation of (Mg,Fe)SiO 3 Post-Perovskite and D'' Anisotropy

An improved setup for radial diffraction experiments at high pressures and high temperatures in a resistive graphite-heated diamond anvil cell

We present an improved setup for the experimental study of deformation of solids at simultaneous high pressures and temperatures by radial x-ray diffraction. This technique employs a graphite resistive heated Mao-Bell type diamond anvil cell for radial x-ray diffraction in combination with a water-cooled vacuum chamber. The new chamber has been developed by the sample environment group at PETRA III and implemented at the Extreme Conditions Beamline P02.2 at PETRA III, DESY (Hamburg, Germany). We discuss applications of the new setup to study deformation of a variety of materials, including ferropericlase, calcium perovskite, bridgmanite, and tantalum carbide, at high-pressure/temperature.

Anisotropic diffusion creep in postperovskite provides a new model for deformation at the core-mantle boundary

The lowermost portion of Earth's mantle (D″) above the core-mantle boundary shows anomalous seismic features, such as strong seismic anisotropy, related to the properties of the main mineral MgSiO3 postperovskite. But, after over a decade of investigations, the seismic observations still cannot be explained simply by flow models which assume dislocation creep in postperovskite. We have investigated the chemical diffusivity of perovskite and postperovskite phases by experiment and ab initio simulation, and derive equations for the observed anisotropic diffusion creep. There is excellent agreement between experiments and simulations for both phases in all of the chemical systems studied. Single-crystal diffusivity in postperovskite displays at least 3 orders of magnitude of anisotropy by experiment and simulation (Da = 1,000 Db; Db ≈ Dc) in zinc fluoride, and an even more extreme anisotropy is predicted (Da = 10,000 Dc; Dc = 10,000 Db) in the natural MgSiO3 system. Anisotropic chemical diffusivity results in anisotropic diffusion creep, texture generation, and a strain-weakening rheology. The results for MgSiO3 postperovskite strongly imply that regions within the D″ region of Earth dominated by postperovskite will 1) be substantially weaker than regions dominated by perovskite and 2) develop a strain-induced crystallographic-preferred orientation with strain-weakening rheology. This leads to strain localization and the possibility to bring regions with significantly varying textures into close proximity by strain on narrow shear zones. Anisotropic diffusion creep therefore provides an attractive alternative explanation for the complexity in observed seismic anisotropy and the rapid lateral changes in seismic velocities in D″.

Competing Deformation Mechanisms in Periclase: Implications for Lower Mantle Anisotropy

Seismic anisotropy is observed above the core-mantle boundary in regions of slab subduction and near the margins of Large Low Shear Velocity Provinces (LLSVPs). Ferropericlase is believed to be the second most abundant phase in the lower mantle. As it is rheologically weak, it may be a dominant source for anisotropy in the lowermost mantle. Understanding deformation mechanisms in ferropericlase over a range of pressure and temperature conditions is crucial to interpret seismic anisotropy. The effect of temperature on deformation mechanisms of ferropericlase has been established, but the effects of pressure are still controversial. With the aim to clarify and quantify the effect of pressure on deformation mechanisms, we perform room temperature compression experiments on polycrystalline periclase to 50 GPa. Lattice strains and texture development are modeled using the Elasto-ViscoPlastic Self Consistent method (EVPSC). Based on modeling results, we find that { 110 } ⟨ 1 1 ¯ 0 ⟩ slip is increasingly activated with higher pressure and is fully activated at ~50 GPa. Pressure and temperature have a competing effect on activities of dominant slip systems. An increasing { 100 } ⟨ 011 ⟩ : { 110 } ⟨ 1 1 ¯ 0 ⟩ ratio of slip activity is expected as material moves from cold subduction regions towards hot upwelling region adjacent to LLSVPs. This could explain observed seismic anisotropy in the circum-Pacific region that appears to weaken near margins of LLVSPs.

A simple variant selection in stress-driven martensitic transformation

The study of orientation variant selection helps to reveal the mechanism and dynamic process of martensitic transformations driven by temperature or pressure/stress. This is challenging due to the multiple variants which may coexist. While effects of temperature and microstructure in many martensitic transformations have been studied in detail, effects of stress and pressure are much less understood. Here, an in situ variant selection study of Mn₂O₃ across the cubic-to-orthorhombic martensitic transformation explores orientation variants at pressures up to 51.5 GPa and stresses up to 5.5 GPa, using diamond anvil cells in radial geometry with synchrotron X-ray diffraction. The diamonds not only exert pressure but also impose stress and cause plastic deformation and texture development. The crystal orientation changes were followed in situ and a {110} _c 〈001〉 _c // (100) _o 〈010〉 _o relationship was observed. Only the {110} _c plane perpendicular to the stress direction was selected to become (100) _o , resulting in a very strong texture of the orthorhombic phase. Contrary to most other martensitic transformations, this study reveals a clear and simple variant selection that is attributed to structural distortions under pressure and stress.

Prediction of Mechanical Twinning in Magnesium Silicate Post-Perovskite

The plastic properties of MgSiO3 post-perovskite are considered to be one of the key issues necessary for understanding the seismic anisotropy at the bottom of the mantle in the so-called D" layer. Although plastic slip in MgSiO3 post-perovskite has attracted considerable attention, the twinning mechanism has not been addressed, despite some experimental evidence from low-pressure analogues. On the basis of a numerical mechanical model, we present a twin nucleation model for post-perovskite involving the emission of 1/6 <110> partial dislocations. Relying on first-principles calculations with no adjustable parameters, we show that {110} twin wall formation resulting from the interaction of multiple twin dislocations occurs at a twinning stress comparable in magnitude to the most readily occurring slip system in post-perovskite. Because dislocation activities and twinning are competitive strain-producing mechanisms, twinning should be considered in future models of crystallographic preferred orientations in post-perovskite to better interpret seismic anisotropy in the lowermost lower mantle.

Note: Modified anvil design for improved reliability in DT-Cup experiments

The Deformation T-Cup (DT-Cup) is a modified 6-8 multi-anvil apparatus capable of controlled strain-rate deformation experiments at pressures greater than 18 GPa. Controlled strain-rate deformation was enabled by replacing two of the eight cubic "second-stage" anvils with hexagonal cross section deformation anvils and modifying the "first-stage" wedges. However, with these modifications approximately two-thirds of experiments end with rupture of the hexagonal anvils. By replacing the hexagonal anvils with cubic anvils and, split, deformation wedge extensions, we restore the massive support to the deformation anvils that were inherent in the original multi-anvil design and prevent deformation anvil failure. With the modified parts, the DT-Cup has an experimental success rate that is similar to that of a standard hydrostatic 6-8 multi-anvil apparatus.

Seismic anisotropy of the D″ layer induced by (001) deformation of post-perovskite

Crystallographic preferred orientation (CPO) of post-perovskite (Mg,Fe)SiO₃ (pPv) has been believed to be one potential source of the seismic anisotropic layer at the bottom of the lower mantle (D″ layer). However, the natural CPO of pPv remains ambiguous in the D″ layer. Here we have carried out the deformation experiments of pPv-(Mg_0.75,Fe_0.25)SiO₃ using synchrotron radial X-ray diffraction in a membrane-driven laser-heated diamond anvil cell from 135 GPa and 2,500 K to 154 GPa and 3,000 K. Our results show that the intrinsic texture of pPv-(Mg_0.75,Fe_0.25)SiO₃ should be (001) at realistic P–T conditions of the D″ layer, which can produce a shear wave splitting anisotropy of ∼3.7% with V_SH>V_SV. Considering the combined effect of both pPv and ferropericlase, we suggest that 50% or less of deformation is sufficient to explain the origin of the shear wave anisotropy observed seismically in the D″ layer beneath the circum-Pacific rim.

The source of the anisotropic layer (D'' layer) at the bottom of the lower mantle remains unclear. Here, using high pressure and temperature experiments, the authors find that seismic anisotropy observed at the D'' layer is caused by 50% deformation of the minerals post-perovskite and ferropericlase.

Modeling defects and plasticity in MgSiO3 post-perovskite: Part 3-Screw and edge [001] dislocations

In this study, we investigate the complex structure of [001] screw and edge dislocation cores in MgSiO₃ post-perovskite at the atomic scale. Both [001] screw and edge dislocations exhibit spontaneous dissociation in (010) into two symmetric partials characterized by the presence of <100> component. In case of edge dislocations, dissociation occurs into ½<101> partials, while for the screw dislocations the <100> component reaches only 15%. Under applied stress, both [001](010) screw and edge dislocations behave similarly. Above the Peierls stress, the two partials glide together while keeping their stacking-fault widths (~11 and ~42 Å for the screw and edge dislocations, respectively) constant. The Peierls stress opposed to the glide of [001](010) screw dislocations is 3 GPa, while that of edge dislocations is 33% lower. Relying on the observed characteristics of the dislocation cores, we estimate the efficiency of [001](010) dislocation glide under the P-T conditions relevant to the lowermost mantle and demonstrate that dislocation creep for this slip system would occur in the so-called athermal regime where lattice friction for the considered slip system vanishes when the temperature rises above the critical T _a value of ~2,000 K.

High-pressure studies with x-rays using diamond anvil cells

Pressure profoundly alters all states of matter. The symbiotic development of ultrahigh-pressure diamond anvil cells, to compress samples to sustainable multi-megabar pressures; and synchrotron x-ray techniques, to probe materials' properties in situ, has enabled the exploration of rich high-pressure (HP) science. In this article, we first introduce the essential concept of diamond anvil cell technology, together with recent developments and its integration with other extreme environments. We then provide an overview of the latest developments in HP synchrotron techniques, their applications, and current problems, followed by a discussion of HP scientific studies using x-rays in the key multidisciplinary fields. These HP studies include: HP x-ray emission spectroscopy, which provides information on the filled electronic states of HP samples; HP x-ray Raman spectroscopy, which probes the HP chemical bonding changes of light elements; HP electronic inelastic x-ray scattering spectroscopy, which accesses high energy electronic phenomena, including electronic band structure, Fermi surface, excitons, plasmons, and their dispersions; HP resonant inelastic x-ray scattering spectroscopy, which probes shallow core excitations, multiplet structures, and spin-resolved electronic structure; HP nuclear resonant x-ray spectroscopy, which provides phonon densities of state and time-resolved Mössbauer information; HP x-ray imaging, which provides information on hierarchical structures, dynamic processes, and internal strains; HP x-ray diffraction, which determines the fundamental structures and densities of single-crystal, polycrystalline, nanocrystalline, and non-crystalline materials; and HP radial x-ray diffraction, which yields deviatoric, elastic and rheological information. Integrating these tools with hydrostatic or uniaxial pressure media, laser and resistive heating, and cryogenic cooling, has enabled investigations of the structural, vibrational, electronic, and magnetic properties of materials over a wide range of pressure-temperature conditions.

Low viscosity and high attenuation in MgSiO3 post-perovskite inferred from atomic-scale calculations

This work represents a numerical study of the thermal activation for dislocation glide of the [100](010) slip system in MgSiO₃ post-perovskite (Mg-ppv) at 120 GPa. We propose an approach based on a one-dimensional line tension model in conjunction with atomic-scale calculations. In this model, the key parameters, namely, the line tension and the Peierls barrier, are obtained from density functional theory calculations. We find a Peierls stress σ_p = 2.1 GPa and a line tension Γ = 9.2 eV/Å, which lead to a kink-pair enthalpy (under zero stress) of 2.69 eV. These values confirm that this slip system bears a very low lattice friction because it vanishes for temperatures above approximately 500 K under mantle conditions. In the Earth’s mantle, high-pressure Mg-ppv silicate is thus expected to become as ductile as ferropericlase. These results confirm the hypothesis of a weak layer in the D″ layer where Mg-ppv is present. Easy glide along [100](010) suggests strong preferred orientations with (010) planes aligned. Highly mobile [100] dislocations are also likely to respond to stresses related to seismic waves, leading to energy dissipation and strong attenuation.

Deformation Behavior across the Zircon-Scheelite Phase Transition

The pressure effects on plastic deformation and phase transformation mechanisms of materials are of great importance to both Earth science and technological applications. Zircon-type materials are abundant in both nature and the industrial field; however, there is still no in situ study of their deformation behavior. Here, by employing radial x-ray diffraction in a diamond anvil cell, we investigate the dislocation-induced texture evolution of zircon-type gadolinium vanadate (GdVO_{4}) in situ under pressure and across its phase transitions to its high-pressure polymorphs. Zircon-type GdVO_{4} develops a (001) compression texture associated with dominant slip along ⟨100⟩{001} starting from 5 GPa. This (001) texture transforms into a (110) texture during the zircon-scheelite phase transition. Our observation demonstrates a martensitic mechanism for the zircon-scheelite transformation. This work will help us understand the local deformation history in the upper mantle and transition zone and provides fundamental guidance on material design and processing for zircon-type materials.

Modeling defects and plasticity in MgSiO3 post-perovskite: Part 1-generalized stacking faults

In this work, we examine the transferability of a pairwise potential model (derived for MgSiO₃ perovskite) to accurately compute the excess energies of the generalized stacking faults (GSF, also called γ-surfaces) in MgSiO₃ post-perovskite. All calculations have been performed at 120 GPa, a pressure relevant to the D″ layer. Taking into account an important aspect of crystal chemistry for complex materials, we consider in detail all possible locations of slip planes in the post-perovskite structure. The γ-surface calculations emphasize the easiness of glide of slip systems with the smallest shear vector [100] and of the [001](010) slip system. Our results are in agreement with previous ab initio calculations. This validates the use the chosen potential model for further full atomistic modeling of dislocations in MgSiO₃ post-perovskite.

Simple metal under tensile stress: layer-dependent herringbone reconstruction of thin potassium films on graphite

While understanding the properties of materials under stress is fundamentally important, designing experiments to probe the effects of large tensile stress is difficult. Here tensile stress is created in thin films of potassium (up to 4 atomic layers) by epitaxial growth on a rigid support, graphite. We find that this "simple" metal shows a long-range, periodic "herringbone" reconstruction, observed in 2- and 3- (but not 1- and 4-) layer films by low-temperature scanning tunneling microscopy (STM). Such a pattern has never been observed in a simple metal. Density functional theory (DFT)simulations indicate that the reconstruction consists of self-aligned stripes of enhanced atom density formed to relieve the tensile strain. At the same time marked layer-dependent charging effects lead to substantial variation in the apparent STM layer heights.

Deformation T-Cup: a new multi-anvil apparatus for controlled strain-rate deformation experiments at pressures above 18 GPa

A new multi-anvil deformation apparatus, based on the widely used 6-8 split-cylinder, geometry, has been developed which is capable of deformation experiments at pressures in excess of 18 GPa at room temperature. In 6-8 (Kawai-type) devices eight cubic anvils are used to compress the sample assembly. In our new apparatus two of the eight cubes which sit along the split-cylinder axis have been replaced by hexagonal cross section anvils. Combining these anvils hexagonal-anvils with secondary differential actuators incorporated into the load frame, for the first time, enables the 6-8 multi-anvil apparatus to be used for controlled strain-rate deformation experiments to high strains. Testing of the design, both with and without synchrotron-X-rays, has demonstrated the Deformation T-Cup (DT-Cup) is capable of deforming 1-2 mm long samples to over 55% strain at high temperatures and pressures. To date the apparatus has been calibrated to, and deformed at, 18.8 GPa and deformation experiments performed in conjunction with synchrotron X-rays at confining pressures up to 10 GPa at 800 °C .

Combined resistive and laser heating technique for in situ radial X-ray diffraction in the diamond anvil cell at high pressure and temperature

To extend the range of high-temperature, high-pressure studies within the diamond anvil cell, a Liermann-type diamond anvil cell with radial diffraction geometry (rDAC) was redesigned and developed for synchrotron X-ray diffraction experiments at beamline 12.2.2 of the Advanced Light Source. The rDAC, equipped with graphite heating arrays, allows simultaneous resistive and laser heating while the material is subjected to high pressure. The goals are both to extend the temperature range of external (resistive) heating and to produce environments with lower temperature gradients in a simultaneously resistive- and laser-heated rDAC. Three different geomaterials were used as pilot samples to calibrate and optimize conditions for combined resistive and laser heating. For example, in Run#1, FeO was loaded in a boron-mica gasket and compressed to 11 GPa then gradually resistively heated to 1007 K (1073 K at the diamond side). The laser heating was further applied to FeO to raise temperature to 2273 K. In Run#2, Fe-Ni alloy was compressed to 18 GPa and resistively heated to 1785 K (1973 K at the diamond side). The combined resistive and laser heating was successfully performed again on (Mg0.9Fe0.1)O in Run#3. In this instance, the sample was loaded in a boron-kapton gasket, compressed to 29 GPa, resistive-heated up to 1007 K (1073 K at the diamond side), and further simultaneously laser-heated to achieve a temperature in excess of 2273 K at the sample position. Diffraction patterns obtained from the experiments were deconvoluted using the Rietveld method and quantified for lattice preferred orientation of each material under extreme conditions and during phase transformation.

A method for intragranular orientation and lattice strain distribution determination

A novel approach to quantifying intragranular distributions is developed and applied to the α → ∊ phase transition in iron. The approach captures both the distribution of lattice orientation within a grain and the orientation dependence of the lattice strain. Use of a finite element discretization over a ball in Rodrigues space allows for the efficient use of degrees of freedom in the numerical approach and provides a convenient framework for gradient-based regularization of the inverse problem. Application to the α → ∊ phase transition in iron demonstrates the utility of the method in that intragranular orientation and lattice strain distributions in the α phase are related to the observed ∊ orientations. Measurement of the lattice strain distribution enables quantitative analysis of the driving forces for ∊ variant selection. The measurement and analysis together indicate quantitatively that the Burgers mechanism is operative under the experimental conditions examined here.

High-pressure experimental geosciences: state of the art and prospects

This paper aims at reviewing the current advancements of high pressure experimental geosciences. The angle chosen is that of in situ measurements at the high pressure (P) and high temperature (T) conditions relevant of the deep Earth and planets, measurements that are often carried out at large facilities (X-ray synchrotrons and neutron sources). Rather than giving an exhaustive catalogue, four main active areas of research are chosen: the latest advancements on deep Earth mineralogy, how to probe the properties of melts, how to probe Earth dynamics, and chemical reactivity induced by increased P-T conditions. For each area, techniques are briefly presented and selected examples illustrate their potentials, and what that tell us about the structure and dynamics of the planet.

Crystal structures of (Mg1-x,Fe(x))SiO3 postperovskite at high pressures

X-ray diffraction experiments on postperovskite (ppv) with compositions (Mg(0.9)Fe(0.1))SiO(3) and (Mg(0.6)Fe(0.4))SiO(3) at Earth core-mantle boundary pressures reveal different crystal structures. The former adopts the CaIrO(3)-type structure with space group Cmcm, whereas the latter crystallizes in a structure with the Pmcm (Pmma) space group. The latter has a significantly higher density (ρ = 6.119(1) g/cm(3)) than the former (ρ = 5.694(8) g/cm(3)) due to both the larger amount of iron and the smaller ionic radius of Fe(2+) as a result of an electronic spin transition observed by X-ray emission spectroscopy (XES). The smaller ionic radius for low-spin compared to high-spin Fe(2+) also leads to an ordered cation distribution in the M1 and M2 crystallographic sites of the higher density ppv structure. Rietveld structure refinement indicates that approximately 70% of the total Fe(2+) in that phase occupies the M2 site. XES results indicate a loss of 70% of the unpaired electronic spins consistent with a low spin M2 site and high spin M1 site. First-principles calculations of the magnetic ordering confirm that Pmcm with a two-site model is energetically more favorable at high pressure, and predict that the ordered structure is anisotropic in its electrical and elastic properties. These results suggest that interpretations of seismic structure in the deep mantle need to treat a broader range of mineral structures than previously considered.

Experimental method for in situ determination of material textures at simultaneous high pressure and high temperature by means of radial diffraction in the diamond anvil cell

We introduce the design and capabilities of a resistive heated diamond anvil cell that can be used for side diffraction at simultaneous high pressure and high temperature. The device can be used to study lattice-preferred orientations in polycrystalline samples up to temperatures of 1100 K and pressures of 36 GPa. Capabilities of the instrument are demonstrated with preliminary results on the development of textures in the bcc, fcc, and hcp polymorphs of iron during a nonhydrostatic compression experiment at simultaneous high pressure and high temperature.

Between a rock and a hot place: the core-mantle boundary

The boundary between the rocky mantle and iron core, almost 2900 km below the surface, is physically the most significant in the Earth's interior. It may be the terminus for subducted surface material, the source of mantle plumes and a control on the Earth's magnetic field. Its properties also have profound significance for the thermochemical and dynamic evolution of the solid Earth. Evidence from seismology shows that D'' (the lowermost few hundred kilometres of the mantle) has a variety of anomalous features. Understanding the origin of these observations requires an understanding of the elastic and deformation properties of the deep Earth minerals. Core-mantle boundary pressures and temperatures are achievable in the laboratory using diamond anvil cell (DAC) apparatus. Such experiments have led to the recent discovery of a new phase, 'post-perovskite', which may explain many hitherto poorly understood properties of D''. Experimental work is also done using analogue minerals at lower pressures and temperatures; these circumvent some of the limits imposed by the small sample size allowed by the DAC. A considerable contribution also comes from theoretical methods that provide a wealth of otherwise unavailable information, as well as verification and refinement of experimental results. The future of the study of the lowermost mantle will involve the linking of the ever-improving seismic observations with predictions of material properties from theoretical and experimental mineral physics in a quantitative fashion, including simulations of the dynamics of the deep Earth. This has the potential to dispel much of the mystery that still surrounds this remote but important region.

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).

Structure and dynamics of Earth's lower mantle

Processes within the lowest several hundred kilometers of Earth's rocky mantle play a critical role in the evolution of the planet. Understanding Earth's lower mantle requires putting recent seismic and mineral physics discoveries into a self-consistent, geodynamically feasible context. Two nearly antipodal large low-shear-velocity provinces in the deep mantle likely represent chemically distinct and denser material. High-resolution seismological studies have revealed laterally varying seismic velocity discontinuities in the deepest few hundred kilometers, consistent with a phase transition from perovskite to post-perovskite. In the deepest tens of kilometers of the mantle, isolated pockets of ultralow seismic velocities may denote Earth's deepest magma chamber.