## Abstract

Bridgmanite, the dominant mineral in the Earth’s lower mantle, crystallizes in the perovskite structure and transforms into post-perovskite at conditions relevant for the D\({}^{{\prime\prime} }\) layer. This transformation affects the dynamics of the Earth’s lowermost mantle and can explain a range of seismic observations. The thickness over which the two phases coexist, however, can extend over 100 km, casting doubt on the assignment of the observed seismic boundaries. Here, experiments show that the bridgmanite to post-perovskite transition in (Mg\({}_{0.86}\),Fe\({}_{0.14}\))SiO\({}_{3}\) is fast on geological timescales. The transformation kinetics, however, affects reflection coefficients of \(P\) and \(S\) waves by more than one order of magnitude. Thick layers of coexisting bridgmanite and post-perovskite can hence be detected using seismic reflections. Morever, the detection and wave period dependence of D\({}^{{\prime\prime} }\) reflections can be used to constrain significant features of the Earth’s lowermost mantle, such as the thickness of the coexistence layer, and obtain information on temperature and grain sizes.

## Introduction

The region \(200\) km above the core-mantle boundary, the D″ layer, is heterogeneous, anisotropic, and acts as a boundary layer for mantle convection^{1}. Seismic reflections in D″ have been interpreted as the result of a solid-state phase change^{2} and, most frequently, attributed to the transition of the main lower mantle mineral bridgmanite from a perovskite (Pv) to a post-perovskite (pPv) structure^{3,4,5,6}. The exact nature of the mechanism of the Pv to pPv transition remains largely unknown^{7,8,9,10} and there is scarce information regarding its kinetics in the literature. Phase transformation kinetics, however, has been shown to strongly affect dynamics in other mantle regions through buoyancy, thermal, and rheological effects^{11,12}. The presence and seismic detectability of the Pv-pPv transition in the Earth’s mantle also remains a matter of debates as chemical composition strongly affect the thickness and pressure-temperature (P/T) range of the Pv-pPv coexistence layer^{13,14,15}. The interpretation of seismic travel times, however, tends to favor the existence of pPv in D″^{16,17,18}. Some discrepancies still exist; the observed velocity contrast, for instance, is larger than that predicted by a Pv-pPv phase transition^{16,19}. Seismic reflection amplitudes and polarities could provide strong additional constrains for deciphering D\({}^{{\prime\prime} }\)composition and processes^{20,21,22} but their interpretation requires a good understanding of the underlying phase transformations, which is lacking for the case of the Pv to pPv conversion.

Here, we report a high pressure, high temperature (HP/HT) study of Pv to pPv transformations in (Mg\({}_{0.86}\),Fe\({}_{0.14}\))SiO\({}_{3}\) conducted in the double-sided laser heated diamond anvil cell. The experimental data is used to constrain the kinetics of the Pv to pPv transformation at D″ P/T conditions. We then evaluate the combined effect of both kinetics and thickness of the Pv-pPv coexistence region on the detectability of the Pv-pPv transition using seismic wave reflections in the Earth’s mantle. Stresses induced by seismic waves can lead to a re-equilibration process^{23,24} and strongly affect the amplitude of the reflected waves. Finally, we suggest new seismic measurements that could help constrain the physical state of the Earth’s lowermost mantle.

## Results

### In situ experiments

A sample of natural enstatite of composition (Mg\({}_{0.86}\),Fe\({}_{0.14}\))SiO\({}_{3}\) was loaded in a diamond anvil cell inside an argon pressure medium (see Methods). Enstatite was first compressed at ambient temperature to \(\approx\! 90\) GPa and converted to bridgmanite by laser heating to above 2000 K for 30 min while scanning around the sample. Bridgmanite was then compressed to the target pressure (between 116 and 130 GPa, Table 1) at ambient temperature. We then activated the conversion from the Pv structure into pPv by further heating the sample and keeping it at a constant temperature in the 1600–2400 K range. To ensure kinetics data quality, all heating runs with temperature gradients over 100 K between both sample sides or changes in temperature over 100 K over time were discarded. Conversion of Pv to pPv was continuously monitored by collecting X-ray diffraction images every 10 s at the P02.2 beamline of the PETRA III synchrotron (Fig. 1). When no change in X-ray diffraction pattern could be observed, the sample was quenched to room temperature. If possible, the experiment was repeated in another location of the same sample which had not yet been fully converted to pPv. Otherwise, a new sample was used. Out of nine attempts, we successfully studied four different samples, with data at 12 P/T points (Table 1), 2 of which with workable Pv to pPv kinetics data. Measurements at higher temperature were not attempted to allow for reliable time measurement with the beamline setup we used. Sample pressures were determined from the unit cell volumes and the pressure-volume-temperature equations of state (PVT-EOS) of Pv^{25} and pPv^{26}, with a relative uncertainty of \(\pm\!\)2 GPa. There are, however, additional uncertainties on absolute pressure due to the inconsistencies in PVT-EOS calibrations^{27,28}. This study will hence report two uncertainties for pressure, \(\pm\!\)2 GPa on relative pressure (i.e. pressure differences between the reported points) and \(\pm\!\)5 GPa on absolute pressures (see Methods). Temperature uncertainty is below 100 K. Resolution on times during data collection is 10 s.

### Kinetics of the Pv to pPv transformation

A typical sequence of diffraction patterns with the Pv to pPv conversion is shown in Fig. 1. As diffraction peaks of Pv decrease in intensity, the pPv diffraction line intensities increase, indicative of the ongoing transformation (Fig. 2). In some cases, we observe a transformation from Pv to pPv. In others, no transformation to pPv was observed during the time of the experiment (Fig. 3a, Table 1), indicating either that kinetics is slower than the experiment duration or that the P/T conditions are not in the pPv stability domain. For P/T points with partial or full conversion to pPv, the resulting proportions vs. time plots (Fig. 2) were analyzed in the framework of the Avrami model^{29} with

where \(\xi (t)\) is the extend of the transformation, itself deduced from the weight proportion of pPv deduced from the x-ray diffraction data with \(\xi =0\) and \(\xi =1\) for the starting and final weight proportion of pPv, respectively, \(n\) the Avrami exponent, and \(\tau\) the characteristic transformation time (Table 1, Methods). The Avrami exponent \(n\) ranges between \(0.87\) and \(1.74\) and, in most cases, is close to 1.5. Such a value for \(n\) suggests a mechanism of instantaneous homogeneous nucleation with a diffusion-controlled growth^{30}. Previous publications have suggested mechanisms involving homogeneous shear for the Pv-pPv transformation^{7,9,10}, with contributions of nucleation and growth^{8}, and this can seem contradictory. It has been shown, however, that transformation mechanisms in silicates can involve both shear and diffusion stages^{31,32}. The present result of \(n\approx 1.5\) hence indicates that the controlling factor for kinetics is a diffusion-controlled stage in the transformation. This suggested mechanism will have to be reinforced by additional future observations. In the meantime, and in the rest of the analysis, we will consider both mechanisms involving shear and nucleation and growth.

The kinetics data of Table 1 are extrapolated to all D″ P/T conditions using either a nucleation and growth model

or a shear transformation model

\({k}_{2}\) and \({k}_{2}^{\prime}\) are constants, \(T\) the absolute temperature, \({Q}_{0}\) and \({V}^{* }\) the absolute activation energy and activation volume of the transformation, \(\Delta P\) the overpressure relative to the Pv/pPv phase boundary, \(R\) the gas constant, and \(\Delta G\) the free energy change of the transformation (Fig. 3b, c). Both were adapted from models in the literature^{33,34,35} (see Methods).

The parameters of Eqs. (2) and (3) can be obtained using unconstrained non-linear least squares curve fitting (Methods). The nucleation and growth model leads to \({Q}_{0}=426\,\pm\,180\) kJ mol^{−1} and \({V}^{* }=14.4\,\pm\,6.4\) cm\({}^{3}\) mol^{−1} while the shear transformation model leads to \({Q}_{0}=437\,\pm\,186\) kJ mol^{−1} and \({V}^{* }=16.1\,\pm\,6.3\) cm\({}^{3}\) mol^{−1}, where the expressed uncertainties account for a \(\pm\! 1\sigma\) standard deviation from the non-linear least squares as well as an effect of \(\pm\! 5\) GPa in the location of the Pv/Pv + pPv phase boundary and \(\pm\! 2\) MPa K^{−1} for the corresponding Clapeyron slope. Both models yield very similar results regarding characteristic transformation times at D“ conditions (Fig. 3b, c). The goodness of the fit of kinetics parameters to Eqs. (2) and (3) does not allow to favor one model over another.

Interestingly, for P/T conditions at which we did not observe the Pv/pPv phase transition, both models predict characteristic transformation times above \(1{0}^{3}\) s, which is longer than the duration of our experiments and is, hence, in agreement with our reported non-observations (Table 1). The slow kinetics of the Pv to pPv transformation for temperatures below 2000 K and pressures below 120 GPa is also consistent with reports of “sluggish” Pv-pPv transformation in the literature and the difficulty of establishing consistent phase diagrams in this P/T region^{13,14,15}. Above \(3000\) K and \(125\) GPa, both kinetics models suggest transformation characteristic times below \(1\) s. With most current HP-HT experimental facilities, such fast kinetics can not be measured and appears instantaneous.

## Discussion

### Upscaling to the D″ layer

Grain size becomes a critical parameter for scaling experimental kinetics results to the Earth’s D″ layer. Diamond anvil cell experiments are performed with grain sizes in the \(\upmu\)m range. Grain sizes in D″ are unknown but samples of the upper mantle typically show grain sizes in the mm range.

Shear transformation are driven by the cooperative and homogeneous movement of many atoms. Such transformation (i.e. Eq. (3)) are somewhat independent of grain size. Transformations controlled by diffusion (i.e. Eq. (2)) typically assume a constant product phase growth rate \(G\). A dimensional analysis indicates that growth rate \(G\), is proportional to the ratio of a characteristic length \(l\) over a characteristic time \(t\). In our problem, the relevant length and time units are the average sample grain size \(d\) and transformation characteristic time \(\tau\). In a first approximation, we can hence define

where \({k}_{0}\) is a constant. We can hence expect that, for transformations controlled by diffusion, transformation times scale with grain sizes (Eq. (4)). For experimental grain sizes and D″ P/T conditions, \(\tau\) ranges from \(1{0}^{-1}\) to 10\({}^{4}\) s (Fig. 3b, c). For mm grains and diffusion-controlled transformation, we then expect \(\tau\) ranging from 10\({}^{2}\) to 10\({}^{7}\) s (Fig. 3d), or 10\({}^{5}\) to 10\({}^{10}\) s (i.e. 300 years) with meter-size grains.

Whatever the pressure, temperature, or grain size, the transformation from Pv to pPv is fast on geological timescales and will not affect long term mantle dynamics. For smaller grain sizes or shear-driven transformations, however, the characteristic times fall in the range of possible interaction with seismic waves. We hence speculate that the Pv to pPv transformation kinetics could be studied with seismic data.

### Detectability of the phase transition

In regions where two phases coexist, stresses induced by seismic waves propagating through the transition loop can locally disrupt the thermodynamically stable system and lead to a re-equilibration process that significantly softens the elastic moduli of the aggregate^{23,24}. If the kinetics of the transformation is not instantaneous, it will affect attenuation and reflection coefficients of the seismic waves. We thus propose to use the kinetics timescales determined experimentally to investigate the effect of the Pv-pPv transition on seismic wave reflection amplitudes.

Following the approach used in previous work^{22,36}, we explore the scenario of vertically incident \(SH\) and \(P\) waves, with a 1–40 s period, and a D″ interface caused by the Pv to pPv transition using a micro-mechanical model of a phase coexistence loop^{24} and the elastic parameters and densities previously computed at 128 GPa and 2800 K^{37} (Methods). Following experimental results^{15}, the Pv/pPv equilibrium coexistence region is assumed to vary over 0–100 km (0–8 GPa). The characteristic transformation times are set to \(\tau\) ranging between \(1{0}^{-1}\) s and \(1{0}^{4}\) s. In such scenarios, we observe significant variations of the reflection coefficients, ranging from a minimum value below 1% up to 33% (Fig. 4).

The detectability of the Pv-pPv transition using seismic waves can be strongly enhanced, particularly if the Pv-pPv coexistence layer extends over several kilometers and transition times are short, which is plausible for small grain sizes or grain size independent kinetics (Fig. 3). Without the effect of kinetics, a thick Pv-pPv coexistence region would not be detected by seismic waves because of the small amplitudes of reflection, and especially for short period waves, but, with the inclusion of fast kinetics, a gradual transition from Pv to pPv can produce large \(P\) and \(SH\) reflections. For a 50 km thick interface, for instance, reflection coefficients increase strongly for periods above 5 and 10 s for \(P\) and \(SH\) waves, respectively, potentially leading to clearly visible seismic reflections.

The experimental results for kinetics were measured in (Mg\({}_{0.86}\),Fe\({}_{0.14}\))SiO\({}_{3}\). The addition of other elements such as Ca or Al can affect the timescales of the transformation. The reflection coefficients are weakly affected by a change of density and elastic properties. The results of Fig. 4, however, strongly depend on kinetics and the thickness of the Pv-pPv coexistence region and both are thoroughly tested in the figures. Figure 4 will hence remain applicable when further experiments quantifying the influence of incorporating other elements on kinetics will become available.

Our calculations only consider vertically incident waves. A variation of the incidence angle would change the reflection coefficients^{16}. Other parameters also affect the amplitude of reflected seismic waves such as anisotropy, that induces azimuthal-dependent reflection coefficients^{38}, or the topography of the reflector, through focusing and defocusing effects^{39}. We suggest that the respective effects of kinetics, anisotropy, and topography on D\({}^{{\prime\prime} }\)reflections could be resolved in a region with a dense distribution of \(P\) and \(S\)-waves reflections, using a range of periods and azimuths. The effect of kinetics is independent of azimuth, unlike that of anisotropy. Topography would generate strongly varying amplitudes of \(P\) and \(S\)-waves within one region, which is not expected for the effect of kinetics.

Our study shows that reflection coefficients can vary over orders of magnitude due to the effect of phase transformation kinetics which could explain previously observed strong variations of D\({}^{{\prime\prime} }\) reflected waves^{16,19}. Our results also show a dependence of the reflection coefficients with the period of the seismic wave. Reflection coefficients for \(SH\) waves, for instance, increase strongly above 2, 10, and 20 s for Pv/pPv coexistence region of 10, 50, and 100 km, respectively. As such, the frequency dependence of the detection and amplitude of D\({}^{{\prime\prime} }\) reflections could be used as a proxy for the existence of a Pv/pPv interface in the deep mantle and constraining the thickness of a potential Pv + pPv coexistence region. Moreover, a measure of the absolute value of the reflection coefficient would constrain the transformation kinetics between both phases. Combined with the results of Fig. 3, such measurement will allow constraining the temperature and grain size in the Earth’s D″ layer.

## Methods

### Experimental details

The starting material was a pure natural enstatite of composition (Mg\({}_{0.86}\),Fe\({}_{0.14}\))SiO\({}_{3}\), as determined from electron microprobe analysis. It was loaded with argon pressure medium in diamond anvil cells (DAC) equipped with two beveled diamond anvils with conical support (100–150 \(\upmu\)m culet diameter with bevels at 7.5° up to 300 \(\upmu\)m diameter). In all cases, sample and pressure medium were contained in a 40–60 \(\upmu\)m diameter hole drilled in a rhenium gasket indented to \(\approx\)25 \(\upmu\)m thickness.

The experiments were performed using monochromatic synchrotron X-ray radiation on beamline P02.2 of the PETRA III synchrotron in Hamburg, Germany^{40}. The X-ray beam was focused to \(\approx\! 2.5\ \upmu\)m both vertically and horizontally, with a wavelength fixed at 0.2918 Å (42.49 keV). Diffraction images were acquired on a Perkin Elmer detector with \(2048\times 2048\) pixels of \(200\,\times\,200\;\upmu {\text{m}}^{2}\) size. The detector to sample distance (551.61 mm), beam center position and detector tilt were determined using a CeO\({}_{2}\) standard.

Sample heating was performed with off-axis Yb-Fiber laser heating system, as provided by the beamline (10–30 \(\upmu\)m heating spot on sample).

### Data processing

The collected X-ray diffraction data were integrated and processed using the Rietveld software MAUD^{41} in order to extract the cell parameters of Pv, pPv, and the Ar pressure medium along with the weight proportions of all phases (Supplementary Fig. 1). The background was corrected for by using polynomial functions. The fitting did not include the effect of stress nor texture, which is not directly observable in such experiment.

### Pressure and temperature

Sample temperatures were estimated on both sample sides by pyrometry of the emitted light with a Czerny–Turner spectrograph equipped with a CCD. Temperature uncertainty on each individual measurement is below 30 K but the sample temperature can evolve during the time of a kinetic experiment. As all data points with more than 100 K dispersion in temperature were discarded, the uncertainty on temperature can hence be safely set to \(\pm\! 100\) K.

Sample pressure at high temperature was determined from the unit cell volumes and the equations of state of Pv^{25} and pPv^{26}. We did notice some fluctuations in the exact values of pressure deduced over time during transformation (varying over \(\pm\!\)2 GPa at maximum). The fluctuations in cell parameters (leading to an apparent change of pressure) could be due to Fe-Mg chemical exchange between Pv and pPv, for instance, but could also be an artifact of the Rietveld refinement on a multi-phase recrystallizing sample.

Accounting for those pressure changes in our kinetic model is a priori feasible but introduces numerical instabilities in the rest of the analysis. For each sample, we hence decided to use an average value of pressure determined from either Pv or pPv over the whole duration of the heating, assigning an uncertainty of \(\pm\!\)2 GPa for each relative pressure value. This value of \(\pm\!\)2 GPa includes the elusive effect of potential Fe-content exchanges between Pv and pPv.

There are, however, additional uncertainties on absolute pressure due to experimental uncertainties and equation of state (EOS) calibrations. Recent publications and reviews^{27,28} discuss the details of determining the Pv-pPv absolute transition pressure. At this point, there is no consensus and conflicting results cannot be solely addressed using self-consistent pressure scales nor accounting for effect of the exact sample chemical composition, with up to 10–15 GPa disagreement between studies. The EOS data we used for bridgmanite^{25} and their reported uncertainties on parameters leads to an uncertainty of \(\pm\!\)3 GPa (\(\pm\! 1\sigma\)) in our pressure-temperature range. Accounting this systematic error and potential disagreements between EOS, we assign a conservative \(\pm\!\)5 GPa on our absolute pressures.

### Fitting of the Avrami model

Data for pressures and temperatures with successful conversion from Pv to pPv are analyzed in the framework of the Avrami model^{29,42,43} (Eq. (1)). The starting, current, and final proportions of pPv define the values for \(\xi\) using:

where \(s\), \(f\), and \(w\) are the starting, final, and current weight proportion of pPv. Experimental results for weight proportions of pPv vs. time and the corresponding \(\xi\) values are listed in Supplementary Tables 1 and 2.

Values for \(n\) and \(\tau\) are obtained using either a non-linear least-squares fit to Eq. (1) or a linear least-square fit of

where \(A\) and \(B\) are adjustable variables and \(A=n\) and \(B=-n\,\mathrm{ln}\, \tau\).

During the experiment, we sometimes noticed a counter-intuitive variation of the diffraction peaks intensity with time, resulting is an apparent reversion of the transformation (Supplementary Fig. 2). These are due to heating of the mechanical parts holding the diamond anvil cell in place inducing unavoidable movements of the sample relative to the heating and X-ray diffraction spot. In this case, the sample was brought back to the correct position by adjusting motors below the diamond anvil cell. The corresponding \(\xi\) vs \(t\) curves were corrected to account for a time offset between both sections of the curves. In the example of Supplementary Fig. 2 the first estimate of \(\tau\) is on the order of 498 s. The final result is \(\tau =502\) s. The effect of such sample drifting corrections is minor on the final results.

### Kinetics model

Characteristic transformation times at constant pressure with successful conversions to pPv are consistent when plotted in an Arrhenius plot, leading to an apparent activation energy of 237 kJ mol^{−1} (Supplementary Fig. 3a). The data, however, does not fall on the same trend for the different hydrostatic pressures. Such effect of pressure is not accounted for when expressing results using a simple Arrhenius law.

In fact, phase transition kinetics increases both with temperature and distance in pressure from the phase boundary. We hence need to introduce a term of overpressure \(\Delta P\), i.e. the distance in pressure from the equilibrium boundary, in addition of temperature in our analysis. Do note that confusion can arise as some kinetics study rely on absolute pressures \(P\) rather than overpressure \(\Delta P\). Such assumption is correct for phase transitions with a low value of Clapeyron slope but not for other cases.

We tested multiple kinetic laws available in the literature for both nucleation and growth and shear transformation models^{30,44,45,46,47,48}. Out of these, two were successful at reproducing all experimental data: the interface-controlled model of refs. ^{33,34} and the shear model of ref. ^{35} (Supplementary Fig. 3b, c).

Kinetics equations are typically fitted in terms a growth rate of a product phase vs. physical parameters. Extracting growth rates from experimental data requires an estimation of grain-sizes. Extracting precise information on grain size is not straightforward in in situ diamond anvil cell experiments. We hence rely on the dimensional analysis of Eq. (4) and express our results in terms of characteristic transformation times \(\tau\) which can be readily extracted from and compared to the true experimental data.

A law for interface-controlled transformation is typically expressed as^{33,34,49}

where \({k}_{1}\) is a constant, \({Q}_{0}\) and \({V}^{* }\) the activation energy and volume for growth, \(\Delta P\) the overpressure, \(\Delta {G}_{r}\) the free energy change of the reaction and \(R\) the universal gas constant. For the Pv to pPv transformation, \(G\) and \(\Delta P\) increase as pressure increases away from the Pv-pPv coexistence boundary, hence the minus sign for the definition of a positive activation volume. This expression becomes Eq. (2) when expressed in terms of characteristic transformation times. \({k}_{2}\) in Eq. (2) is a constant that depends on grain size with \({k}_{2}={k}_{0}d/{k}_{1}\), where \({k}_{1}\) is the constant in Eq. (7), \({k}_{0}\) the constant in Eq. (4), and \(d\) the grain size.

Multiple expressions have been suggested for shear transformation models^{35}. Most depend on parameters which are unknown for the Pv-pPv transformation, i.e. because they depend on the exact transformation mechanism. Models derived for the calcite \(\iff\) aragonite transition, for instance^{35}, assume a mechanism similar to that proposed for the Pv to pPv transformation^{7,8,10}, with a transformation controlled by the development of stacking faults of the daughter phase into the parent phase and kinetics controlled by the motion of partial dislocations associated with these fault. The expression in Eq. (3) is simplified from this model and was found to fit the data reasonably well. In the expression of Eq. (3), \(\exp \left[{Q}_{0}/(RT)\right]\) accounts for the thermal activation and \(\sinh \left[{V}^{* }\Delta P/(RT)\right]\) models the effect of distance from the phase boundary.

### Fitting of kinetics models

For both the nucleation and growth and shear transformation mechanisms, \({k}_{2}\), \({Q}_{0}\), and \({V}^{* }\) were adjusted to the data of Table 1.

Due to the lack of a proper thermodynamic database for Pv and pPv, and as in previous studies^{49}, \(\Delta {G}_{r}\) is estimated from \(\Delta P\Delta V\), where \(\Delta P\) is the overpressure and \(\Delta V\) is the molar volume change of the transformation. \(\Delta V\) was determined directly from the X-ray diffraction data. \(\Delta P\) was calculated by comparing the pressure difference between the measured P/T point and the Pv/Pv+pPv equilibrium boundary.

We obtain the fitting parameters of Eqs. (2) and (3) using unconstrained non-linear least squares curve fitting, relying on python script and a Levenberg-Marquardt optimization as implemented in the scipy.optimize.curve_fit routine in scipy^{50}. The script has been made open-source and can be found online at [https://github.com/smerkel/kinetics-py]. Standard deviation errors on the fit parameters are calculated from the diagonal coefficients of the covariance matrix. Experimental uncertainties in pressure and temperature are not directly accounted for in the Levenberg-Marquardt optimization. They are, however, accounted for indirectly by changing the location of the Pv/Pv + pPv equilibrium boundary (see below).

Calculations of overpressures are sensitive to the location of the Pv/Pv+pPv equilibrium boundary, which is ill-defined for (Mg,Fe)SiO\({}_{3}\). Here, we define the Pv/Pv+pPv equilibrium boundary using a reference point and a Clapeyron slope. The reference point, belonging to the Pv/Pv+pPv equilibrium boundary, was extracted from the extrapolation of ref. ^{14} and is equal to \([{P}_{e};\;{T}_{e}]\) = [128 GPa; 3300 K]. Estimation of the Clapeyron slope of the Pv/pPv equilibrium boundary range between 5.0 and 11.5 MPa K^{−1}^{7,13,27,28}, somewhat consistent with a seismological estimate of 6 MPa K^{−1}^{2}, with recent measurements leading to 8.5(4) MPa K^{−1}^{28}.

According to our estimate for absolute pressures, we test the effect of a \(\pm\! 5\) GPa error in \({P}_{e}\). We also explore Clapeyron slopes of 8.5 \(\pm\!\) 2.0 MPa K^{−1}, which encompases most publish values in the litterature. Extrapolated transformation characteristic times away from the boundary do not change significantly (Supplementary Figs. 4 and 5), as for the numerical values of \({k}_{2}\), \({Q}_{0}\), and \({V}^{* }\) (Supplementary Table 3). Errors reported for \({k}_{2}\), \({Q}_{0}\), and \({V}^{* }\) in the Results section are the sum of the statistical errors of Supplementary Table 3 and that induced by an error of \(\pm\! 5\) GPa in absolute pressures and \(\pm\! 2\) MPa K^{−1} for the Clapeyron slope.

### Relevance of activation volume and energy values

The apparent activation energy obtained from the Arrhenius plot for the Pv/(Pv + pPv) transition, 237 kJ mol^{−1}, is consistent with other mantle phase transformations^{36,49,51,52,53,54,55}. It is, however, significantly lower than that predicted from numerical modeling^{9} indicating, probably, differences in the transformation mechanism. The absolute activation energies for the nucleation and growth and shear models refined from Eqs. (5) and (6) are 426 and 437 KJ mol^{−1}, respectively, two times higher than the apparent activation energy. The absolute activation energy accounts for the fact experimental measurements were performed away from the equilibrium boundary, which is not taken into account in the apparent activation energies.

The absolute activation energy accounts for the effect of overpressure relative to the phase boundary (e.g. Supplementary Fig. 3b, c) and are hence more appropriate for extrapolations to D“ conditions. The activation volumes for the nucleation and growth and shear models are 14.4 and 16.1 cm\({}^{3}\) mol^{−1}, respectively. This value is significantly larger than that measured for phase transformations in (Mg,Fe)\({}_{2}\)SiO\({}_{4}\)^{36,49}. It is, however, lower than those reported for other transitions such the calcite to aragonite transformation^{35}.

### Reflection coefficients of vertically incident SH and P waves

The propagation of the seismic front within a phase loop induces a small perturbation of pressure that can significantly soften the elastic moduli of the aggregate^{23}. Hence, body waves crossing a phase loop will experience a strong attenuation that will affect their propagation. This attenuation can be predicted using a mechanical model^{24} and depends on the period of the wave relative to the kinetics of the phase transition. The effect is maximum when both are of the same order of magnitude. As shown previously^{56}, this softening of elastic moduli in the wave front will also affect the reflectivity of the interface.

Reflection coefficients of body waves are computed assuming a phase loop separating two elastic half-spaces divided into multiple sublayers with a linearly increasing proportion of the high pressure phase with depth. Applying Eqs. (13) and (25) of ref. ^{24} and using the elastic parameters and densities previously computed at 128 GPa and 2800 K in ref. ^{37}, we obtain the depth-dependence of the elastic moduli through the loop presented in Supplementary Fig. 6. Density is assumed to vary linearly through the loop. Reflection and transmission coefficients are evaluated by propagating elastic waves using the Thomson–Haskell matricial propagator method^{57,58} solving for continuity of displacement and stress at each interface. The obtained coefficients depend on the elastic parameters inside the loop as well as the wave frequency and the phase transformation kinetics.

Here, we compute reflection coefficients for vertically incident \(SH\) and \(P\) waves of period 1 to 40 s. The Pv/pPv coexistence region is set to 0, 10, 50, or 100 km (corresponding to 0–8 GPa) and the characteristic transformation times are set to \(\tau\) ranging between \(1{0}^{-1}\) s, and \(1{0}^{4}\) s. For a sharp interface (no phase loop between Pv and pPv), reflections coefficients do not depend on kinetics nor on wave period with \({R}_{P}\,=\,0.5 \%\) and \({R}_{SH}\,=\,2.5 \%\) and can also be compared with results of the Zoeppritz equations^{59}. Using a gradient layer and no effect of kinetics, the amplitude of reflected wave should be reduced for thicker Pv-pPv coexistence layers. Due to phase transformation kinetics, however, the reflection coefficients depend on wave period and can be either reduced for slow transition times or enhanced for fast transition times. Results of the calculations are shown in Fig. 4.

Reflection coefficients depend on the angle of incidence. Using the parameters used in this work and the Zoeppritz equations one can show that reflection coefficient will undergo a polarity change, with small reflection coefficients in its vicinity. Modeling the full effects of both kinetics and incidence angle is beyond the scope of this work. As a starting point, we hence assume that the kinetics will modulate the reflection coefficient as shown in Fig. 4. Full calculations of the reflection coefficient, accounting for both kinetics and incidence angle of the \(P\) and \(SH\) wave will have to be performed in future work to confirm our assumption.

## Data availability

The raw diffraction data used to produce Fig. 1 and the results presented in this paper will be made available by the corresponding author upon reasonable request. The data used to produce Fig. 2 is available in Supplementary Tables 1 and 2. The data used for producing Fig. 3 as well as Supplementary Figs. 3, 4, 5 is in Table 1.

## Code availability

The python code used to fit the kinetics parameters of Eqs. (2) and (3) and generate Fig. 3 and Supplementary Figs. 4 and 5 is available online at [https://github.com/smerkel/kinetics-py]. The program used to generate Fig. 4 will be made available by the corresponding author upon reasonable request.

## References

- 1.
Garnero, E. J. & McNamara, A. K. Structure and dynamics of Earthas lower mantle.

*Science***320**, 626–628 (2008). - 2.
Sidorin, I., Gurnis, M. & Helmberger, D. V. Evidence for a ubiquitous seismic discontinuity at the base of the mantle.

*Science***286**, 1326–1331 (1999). - 3.
Murakami, M., Hirose, K., Kawamura, K., Sata, N. & Ohishi, Y. Post-perovskite phase transition in MgSiO\({}_{3}\).

*Science***304**, 855–858 (2004). - 4.
Oganov, A. R. & Ono, S. Theoretical and experimental evidence for a post-perovskite phase of MgSiO

_{3}in Earthas D″ layer.*Nature***430**, 445–448 (2004). - 5.
Hernlund, J. W., Thomas, C. & Tackley, P. J. A doubling of the post-perovskite phase boundary and structure of the Earthas lowermost mantle.

*Nature***434**, 882–886 (2005). - 6.
Nowacki, A., Walker, A., Wookey, J. & Kendall, J.-M. Evaluating post-perovskite as a cause of D” anisotropy in regions of palaeosubduction.

*Geophys. J. Int.***192**, 1085–1090 (2013). - 7.
Tsuchiya, T., Tsuchiya, J., Umemoto, K. & Wentzcovitch, R. M. Phase transition in MgSiO\({}_{3}\) perovskite in the Earth’s lower mantle.

*Earth Planet. Sci. Lett.***224**, 241–248 (2004). - 8.
Oganov, A. R., Martoňák, R., Laio, A., Raiteri, P. & Parrinello, M. Anisotropy of Earthas D\({}^{{\prime\prime} }\) layer and stacking faults in the MgSiO

_{3}post-perovskite phase.*Nature***438**, 1142–1144 (2005). - 9.
Zahn, D. Nucleation mechanism and kinetics of the perovskite to post-perovskite transition of MgSiO\({}_{3}\) under extreme conditions.

*Chem. Phys. Lett.***573**, 5–7 (2013). - 10.
Dobson, D. P. et al. Strong inheritance of texture between perovskite and post-perovskite in the D\({}^{{\prime\prime} }\) layer.

*Nat. Geosci.***6**, 575–578 (2013). - 11.
Bina, C. R., Stein, S., Marton, F. C. & Ark, E. M. V. Implications of slab mineralogy for subduction dynamics.

*Phys. Earth Planet. Inter.***127**, 51–66 (2001). - 12.
Faccenda, M. & Zilio, L. D. The role of solid-solid phase transitions in mantle convection.

*Lithos***268–271**, 198–224 (2017). - 13.
Catalli, K., Shim, S. & Prakapenka, V. Thickness and clapeyron slope of the post-perovskite boundary.

*Nature***462**, 782–785 (2009). - 14.
Andrault, D. et al. Experimental evidence for perovskite and post-perovskite coexistence throughout the whole D\({}^{{\prime\prime} }\) region.

*Earth Planet. Sci. Lett.***293**, 90–96 (2010). - 15.
Grocholski, B., Catalli, K., Shim, S.-H. & Prakapenka, V. Mineralogical effects on the detectability of the post-perovskite boundary.

*PNAS***109**, 2275–2279 (2012). - 16.
Cobden, L., Thomas, C. & Trampert, J.

*Seismic Detection of Post-perovskite Inside the Earth*pp. 391–440. (Springer, Cham, 2015). - 17.
Kito, T., Rost, S., Thomas, C. & Garnero, E. J. New insights into the P- and S-wave velocity structure of the D\({}^{{\prime\prime} }\) discontinuity beneath the cocos plate.

*Geophys. J. Int.***169**, 631–645 (2007). - 18.
Koelemeijer, P., Schuberth, B., Davies, D., Deuss, A. & Ritsema, J. Constraints on the presence of post-perovskite in Earth’s lowermost mantle from tomographic-geodynamic model comparisons.

*Earth Planet. Sci. Lett.***494**, 226–238 (2018). - 19.
Wysession, M.et al. The D\({}^{{\prime\prime} }\) discontinuity and its implications. in

*The Core-Mantle Boundary Region*(eds Gurnis, M., Wysession, M., Knittle, E. and Buffet, B.) 273–297 (American Geophysical Union, Washington, DC, 1998). - 20.
Cobden, L. & Thomas, C. The origin of D\({}^{{\prime\prime} }\) reflections: a systematic study of seismic array data sets.

*Geophys. J. Int.***194**, 1091–1118 (2013). - 21.
Thomas, C., Wookey, J., Brodholt, J. & Fieseler, T. Anisotropy as cause for polarity reversals of D\({}^{{\prime\prime} }\) reflections.

*Earth Planet. Sci. Lett.***307**, 369–376 (2011). - 22.
Durand, S., Chambat, F., Matas, J. & Ricard, Y. Constraining the kinetics of mantle phase changes with seismic data.

*Geophys. J. Int.***189**, 1557–1564 (2012). - 23.
Li, L. & Weidner, D. J. Effect of phase transitions on compressional-wave velocities in the Earthas mantle.

*Nature***454**, 984–986 (2008). - 24.
Ricard, Y., Matas, J. & Chambat, F. Seismic attenuation in a phase change coexistence loop.

*Phys. Earth Planet. Inter.***176**, 124–131 (2009). - 25.
Wolf, A. S., Jackson, J. M., Dera, P. & Prakapenka, V. B. The thermal equation of state of (Mg,Fe)SiO\({}_{3}\) bridgmanite (perovskite) and implications for lower mantle structures.

*J. Geophys. Res. Solid Earth***120**, 7460–7489 (2015). - 26.
Guignot, N., Andrault, D., Morard, G., Bolfan-Casanova, N. & Mezouar, M. Thermoelastic properties of post-perovskite phase MgSiO\({}_{3}\) determined experimentally at core-mantle boundary P-T conditions.

*Earth Planet. Sci. Lett.***256**, 162–168 (2007). - 27.
Hirose, K., Wentzcovitch, R., Yuen, D. and Lay, T. Mineralogy of the deep mantle—the post-perovskite phase and its geophysical significance. in

*Treatise on Geophysics*(ed. Schubert, G.) Vol. 2, 85–115 (Elsevier, 2015). - 28.
Sun, N. et al. Phase transition and thermal equations of state of (Fe,Al)-bridgmanite and post-perovskite: implication for the chemical heterogeneity at the lowermost mantle.

*Earth Planet. Sci. Lett.***490**, 161–169 (2018). - 29.
Avrami, M. Granulation, phase change, and microstructure kinetics of phase change.

*III. J. Chem. Phys.***9**, 177–184 (1941). - 30.
Christian, J.

*Transformations in metals and alloys*. (Pergamon, New York, 1975). - 31.
Kerschhofer, L., Rubie, D., Sharp, T., McConnell, J. & Dupas-Bruzek, C. Kinetics of intracrystalline olivine-ringwoodite transformation.

*Phys. Earth Planet. Int.***121**, 59–76 (2000). - 32.
Chen, J., Weidner, D. J., Parise, J. B., Vaughan, M. T. & Raterron, P. Observation of cation reordering during the olivine-spinel transition in fayalite by In Situ synchrotron X-ray diffraction at high pressure and temperature.

*Phys. Rev. Lett.***86**, 4072–4075 (2001). - 33.
Turnbull, D. Phase changes.

*Solid state physics***3**, 225–306 (1956). - 34.
Carlson, W. D. & Rosenfeld, J. L. Optical determination of topotactic aragonite-calcite growth kinetics: metamorphic implications.

*J. Geol.***89**, 615–638 (1981). - 35.
Sotin, C. & Madon, M. Generalized nonlinear inversion of kinetics data: application to the calcite-aragonite transformation.

*Phys. Earth Planet. Inter.***52**, 159–171 (1988). - 36.
Perrillat, J. et al. Kinetics of the olivine-ringwoodite transformation and seismic attenuation in the Earth’s mantle transition zone.

*Earth Planet. Sci. Lett.***433**, 360–369 (2016). - 37.
Wookey, J., Stackhouse, S., Kendall, J., Brodholt, J. & Price, G. D. Efficacy of the post-perovskite phase as an explanation for lowermost-mantle seismic properties.

*Nature***438**, 1004–1007 (2005). - 38.
Pisconti, A., Thomas, C. and Wookey, J. Discriminating between causes of D\({}^{{\prime\prime} }\) anisotropy using reflections and splitting measurements for a single path.

*J. Geophys. Res. Solid Earth***124**, 4811–4830 (2019). https://doi.org/10.1029/2018JB01699 - 39.
Thomas, C. & Weber, M. P velocity heterogeneities in the lower mantle determined with the German Regional Seismic Network: improvement of previous models and results of 2D modelling.

*Phys. Earth Planet. Inter.***101**, 105–117 (1997). - 40.
Liermann, H.-P. et al. The extreme conditions beamline P02.2 and the extreme conditions science infrastructure at PETRA III.

*J. Synchrotron Rad.***22**, 908–924 (2015). - 41.
Lutterotti, L., Bortolotti, M., Ischia, G., Lonardelli, I. & Wenk, H.-R. Rietveld texture analysis from diffraction images.

*Z. Kristallogr., Suppl.***26**, 125–130 (2007). - 42.
Avrami, M. Kinetics of phase change. I: General theory.

*J. Chem. Phys.***7**, 1103–1112 (1939). - 43.
Avrami, M. Kinetics of phase change. II transformation-time relations for random distribution of nuclei.

*J. Chem. Phys.***8**, 212–224 (1940). - 44.
Poirier, J. On the kinetics of olivine-spinel transition.

*Phys. Earth Planet. Inter.***26**, 179–187 (1981). - 45.
Carlson, W. D. Aragonite-calcite nucleation kinetics: an application and extension of Avrami transformation theory.

*J. Geol.***91**, 57–71 (1983). - 46.
Madon, M. & Gillet, P. A theoretical approach to the kinetics of calcite \(\rightleftharpoons\) aragonite transition: application to laboratory experiments.

*Earth Planet. Sci. Lett.***67**, 400–414 (1984). - 47.
Rubie, D. and Thompson, A. Kinetics of metamorphic reactions at elevated temperatures and pressures: an appraisal of available experimental data. in

*Metamorphic Reactions*, 27–79 (Springer, 1985). - 48.
Liu, F., Sommer, F., Bos, C. & Mittemeijer, E. Analysis of solid state phase transformation kinetics: models and recipes.

*Int. Mater. Rev.***52**, 193–212 (2007). - 49.
Kubo, T., Ohtani, E. & Funakoshi, K -i. Nucleation and growth kinetics of the α-β transformation in Mg

_{2}SiO_{4}determined by in situ synchrotron powder X-ray diffraction.*Am. Mineral.***89**, 285–293 (2004). - 50.
Jones, E., et al. SciPy: Open source scientific tools for Python (2001). http://www.scipy.org/.

- 51.
Perrillat, J. P., Daniel, I., Lardeaux, J.-M. & Cardon, H. Kinetics of the coesite-quartz transition: application to the exhumation of ultrahigh-pressure rocks.

*J. Petrol.***44**, 773–788 (2003). - 52.
Perrillat, J.-P. et al. Kinetics of antigorite dehydration: a real-time X-ray diffraction study.

*Earth Planet. Sci. Lett.***236**, 899–913 (2005). - 53.
Perrillat, J. P. Kinetics of high-pressure mineral phase transformations using in situ time-resolved x-ray diffraction in the Paris-Edinburgh cell: a practical guide for data acquisition and treatment.

*Mineral. Mag.***72**, 683–695 (2008). - 54.
Diedrich, T., Sharp, T., Leinenweber, K. & Holloway, J. The effect of small amounts of H

_{2}O on olivine to ringwoodite transformation growth rates and implications for subduction of metastable olivine.*Chem. Geol.***262**, 87–99 (2009). - 55.
Perrillat, J. et al. Mechanism and kinetics of the α-β transition in San Carlos olivine Mg

_{1.8}Fe_{0.2}SiO_{4}.*J. Geophys. Res. Solid Earth***118**, 110–119 (2013). - 56.
Durand, S., Chambat, F., Matas, J. & Ricard, Y. Constraining the kinetics of mantle phase changes with seismic data.

*Geophys. J. Int.***189**, 1557–1564 (2012). - 57.
Haskell, N. A. The dispersion of surface waves on multilayered media.

*Bull. Seismol. Soc. Am.***43**, 17–34 (1953). - 58.
Kennett, B. L.

*Seismic Wave Propagation in Stratified Media*. (Australian National University Electronic Press, Canberra, 2009). - 59.
Zoeppritz, K. Über Reflexion und Durchgang seismischer Wellen durch Unstetigkeitsflächen.

*Nachrichten von der Königlichen Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-physikalische Klasse***1919**, 66–84 (1919). - 60.
Lutterotti, L., Matthies, S. & Wenk, H.-R. Maud: a friendly java program for materials analysis using diffraction.

*IUCr: Newsletter of the CPD***21**, 14–15 (1999).

## Acknowledgements

We wish to thank S. Bellayer for assistance with the electron probe analysis and J-P. Perrillat and M. Roskosz for useful discussions. The study was financed from a grant by the Université de Lille and the Région Hauts-de-France for C.L., with support from the Programme National de Planétologie (PNP) of the CNRS, the bilateral ANR-DFG TIMEleSS project (ANR-17-CE31-0025; TH 1530/18-1; SA 2585/3-1), the bilateral PROCOPE Grant n. 40555PC, and the ARCHI-CM project. The project ARCHI-CM of Chevreul Institute (FR 2638) is a grant from the Ministère de l’Enseignement Supérieur et de la Recherche, Région Nord-Pas de Calais and European Regional Development Fund (ERDF). The electron probe microanalysis (EPMA) facility in Lille (France) is supported by the European Regional Development Fund (ERDF). Parts of this research were carried out at the P02.2 beamline of the PETRA III synchrotron at DESY, a member of the Helmholtz Association (HGF). S.D. was supported by DFG grant HAADES DU1634/1-1. We would like to thank H-P. Liermann for assistance during the experiment. The research leading to this result has been supported by the project CALIPSOplus under the Grant Agreement 730872 from the EU Framework Programme for Research and Innovation HORIZON 2020.

## Author information

### Affiliations

### Contributions

The study was originally designed by S.M. and N.H., with input from D.A. C.L., S.M., D.A., and Z.K. conducted the synchrotron experiment. C.L. processed the experimental data and built the kinetics model with the assistance of S.M., N.H, and D.A. S.D. and C.T. provided the seismic reflection model and input for the relevance of the results for seismic observations. S.M., C.L., and N.H. wrote the paper with input and discussion from all co-authors.

### Corresponding author

## Ethics declarations

### Competing interests

The authors declare no competing interests

## Additional information

**Peer review information** *Nature Communications* thanks Jiuhua Chen, Artem Oganov and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Publisher’s note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Supplementary information

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Langrand, C., Andrault, D., Durand, S. *et al.* Kinetics and detectability of the bridgmanite to post-perovskite transformation in the Earth's D″ layer.
*Nat Commun* **10, **5680 (2019). https://doi.org/10.1038/s41467-019-13482-x

Received:

Accepted:

Published:

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.