Kinetics and detectability of the bridgmanite to post-perovskite transformation in the Earth's D″ layer

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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{{\prime\prime} }$$\end{document}″ 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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{0.86}$$\end{document}0.86,Fe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{0.14}$$\end{document}0.14)SiO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{3}$$\end{document}3 is fast on geological timescales. The transformation kinetics, however, affects reflection coefficients of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P$$\end{document}P and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S$$\end{document}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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{{\prime\prime} }$$\end{document}″ 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.

T he 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 postperovskite (pPv) structure [3][4][5][6] . The exact nature of the mechanism of the Pv to pPv transition remains largely unknown 7-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 00 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 % 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 ±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, ±2 GPa on relative pressure (i.e. pressure differences between the reported points) and ±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 ξðtÞ is the extend of the transformation, itself deduced from the weight proportion of pPv deduced from the x-ray diffraction data with ξ ¼ 0 and ξ ¼ 1 for the starting and final weight proportion of pPv, respectively, n the Avrami exponent, and τ 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 % 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 0 2 are constants, T the absolute temperature, Q 0 and V Ã the absolute activation energy and activation volume of the transformation, ΔP the overpressure relative to the Pv/pPv phase   (Table 1). Insets: representative ±1σ uncertainties on the experimental data.
boundary, R the gas constant, and Δ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 ± 180 kJ mol −1 and V Ã ¼ 14:4 ± 6:4 cm 3 mol −1 while the shear transformation model leads to Q 0 ¼ 437 ± 186 kJ mol −1 and V Ã ¼ 16:1 ± 6:3 cm 3 mol −1 , where the expressed uncertainties account for a ±1σ standard deviation from the non-linear least squares as well as an effect of ±5 GPa in the location of the Pv/Pv + pPv phase boundary and ±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 10 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 μ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 τ. 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, τ ranges from 10 À1 to 10 4 s (Fig. 3b, c). For mm grains and diffusion-controlled transformation, we then expect τ 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 τ ranging between 10 À1 s and 10 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  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 00 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 00 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 00 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 μm culet diameter with bevels at 7.5°up to 300 μm diameter). In all cases, sample and pressure medium were contained in a 40-60 μm diameter hole drilled in a rhenium gasket indented to %25 μ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 % 2:5 μ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 2048 pixels of 200 200 μ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 μ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 ±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 ±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 multiphase 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 ±2 GPa for each relative pressure value. This value of ±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 ±3 GPa ( ±1σ) in our pressure-temperature range. Accounting this systematic error and potential disagreements between EOS, we assign a conservative ±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 ξ 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 ξ values are listed in Supplementary Tables 1 and 2.
Values for n and τ are obtained using either a non-linear least-squares fit to Eq. (1) or a linear least-square fit of ln ln where A and B are adjustable variables and A ¼ n and B ¼ Àn ln τ.
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 ξ 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 τ is on the order of 498 s. The final result is τ ¼ 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 Δ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 Δ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 τ 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, ΔP the overpressure, ΔG r the free energy change of the reaction and R the universal gas constant. For the Pv to pPv transformation, G and Δ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 () 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 Q 0 =ðRTÞ ½ accounts for the thermal activation and sinh V Ã ΔP=ðRTÞ ½ 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 , ΔG r is estimated from ΔPΔV, where ΔP is the overpressure and ΔV is the molar volume change of the transformation. ΔV was determined directly from the X-ray diffraction data. Δ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 nonlinear 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 −17,13,27,28 , somewhat consistent with a seismological estimate of 6 MPa K −12 , with recent measurements leading to 8.5(4) MPa K −128 .
According to our estimate for absolute pressures, we test the effect of a ±5 GPa error in P e . We also explore Clapeyron slopes of 8.5 ± 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 ±5 GPa in absolute pressures and ±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 τ ranging between 10 À1 s, and 10 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.