# Periclase deforms more slowly than bridgmanite under mantle conditions

Patrick Cordier, Karine Gouriet, Timmo Weidner, James van Orman, Olivier

Castelnau, Jennifer M. Jackson, Philippe Carrez

► To cite this version:

Patrick Cordier, Karine Gouriet, Timmo Weidner, James van Orman, Olivier Castelnau, et al. Periclase deforms more slowly than bridgmanite under mantle conditions. Nature, 2023, Nature, 613 (7943), pp.303-307. 10.1038/s41586-022-05410-9. hal-03944762

HAL Id: hal-03944762

https://hal.univ-lille.fr/hal-03944762

Submitted on 18 Jan 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Periclase deforms more slowly than bridgmanite under mantle conditions

https://doi.org/10.1038/s41586-022-05410-9

Received: 10 February 2022

Accepted: 3 October 2022

Published online: 11 January 2023

Open access

<figure><img src="image_2.png" /></figure>

Check for updates

Patrick Cordier<fnref n="1" />,²,³, Karine Gouriet<fnref n="1" />, Timmo Weidner<fnref n="1" />, James Van Orman³, Olivier Castelnau⁴, Jennifer M. Jackson<fnref n="5" /> & Philippe Carrez<fnref n="1" />

Transport of heat from the interior of the Earth drives convection in the mantle, which involves the deformation of solid rocks over billions of years. The lower mantle of the Earth is mostly composed of iron-bearing bridgmanite MgSiO₃ and approximately 25% volume periclase MgO (also with some iron). It is commonly accepted that ferropericlase is weaker than bridgmanite<fnref n="1" />. Considerable progress has been made in recent years to study assemblages representative of the lower mantle under the relevant pressure and temperature conditions²,³. However, the natural strain rates are 8 to 10 orders of magnitude lower than in the laboratory, and are still inaccessible to us. Once the deformation mechanisms of rocks and their constituent minerals have been identified, it is possible to overcome this limitation thanks to multiscale numerical modelling, and to determine rheological properties for inaccessible strain rates. In this work we use 2.5-dimensional dislocation dynamics to model the low-stress creep of MgO periclase at lower mantle pressures and temperatures. We show that periclase deforms very slowly under these conditions, in particular, much more slowly than bridgmanite deforming by pure climb creep. This is due to slow diffusion of oxygen in periclase under pressure. In the assemblage, this secondary phase hardly participates in the deformation, so that the rheology of the lower mantle is very well described by that of bridgmanite. Our results show that drastic changes in deformation mechanisms can occur as a function of the strain rate.

This work aims to model steady-state creep of periclase (MgO) under lower mantle conditions, with an emphasis on the role of extremely low strain rates. It is based on the extensive knowledge we have of the deformation of this oxide following (1) the large number of experimental studies carried out on this ceramic since the 1970s, including the more recent mineral physics studies taking into account the effect of pressure, and (2) the important modelling effort carried out in recent years. At room pressure, periclase is easily deformed by glide of $1/2 < 110°$ dislocations on {110} planes. The thermal activation of this mechanism can now be modelled from the atomic scale⁴, even at low strain rates<fnref n="5" />. Pressure strongly affects this mechanism and induces a transition of glide towards the {100} planes above 50 GPa⁶⁻⁹. At high temperature, dislocations glide without friction (athermal regime) and the flow stress results from dislocation interactions, leading, at large strains and strain rates, to a pronounced hardening¹⁰. At low strain rates, ionic diffusion can enhance recovery mechanisms and allow steady-state creep¹¹. In this study, we model the influence of mantle pressures on steady-state creep of periclase by 2.5-dimensional (2.5D) dislocation dynamics (DD). We focus on four pressures: 30, 60, 90 and 120 GPa, which, along the geotherm¹², correspond to 2,000, 2,300, 2,500 and 2,800 K respectively, at around depths of 800, 1,500, 2,100 and 2,700 km respectively.

## Diffusion coefficients

Dislocation climb in MgO is enabled by diffusion of Mg and O atoms, and the climb rate is controlled by the slower diffusing species. Diffusion in MgO has been studied for decades, and the mechanisms and rates of cation and anion diffusion are well established¹³. Mg and O diffuse using vacant sites on the cation and anion sublattices, respectively. Because the concentration of cation vacancies is much higher than that of anion vacancies, in both nominally pure synthetic and natural periclase samples, O diffusion is orders of magnitude slower than Mg diffusion, and is the rate-limiting step in dislocation climb. Whereas Mg diffusion is controlled by extrinsic vacancies, produced to charge-balance trivalent (and tetravalent) cation impurities, O diffusion is intrinsic, enabled by thermally produced Mg and O vacancy pairs (Schottky defects).

To calculate the oxygen self-diffusion coefficient, $D_{Ox}^{sd}$, over the full range of temperatures and pressures relevant to Earth's mantle, we use the expression presented by Ita and Cohen¹⁴, which is derived from their theoretical results on vacancy formation and migration,

$$ \ln D_{Ox}^{sd}(P, T) = \ln(a^2 v_a) + \frac{S_0 + P S'_0}{k_B} - \frac{E_0 + P V_0 + P^2 V'_0}{k_B T} \quad (1) $$

<footnote n="1">¹Univ. Lille, CNRS, INRAE, Centrale Lille, UMR 8207 - UMET - Unité Matériaux et Transformations, Lille, France. ²Institut Universitaire de France, Paris, France. ³Department of Earth, Environmental and Planetary Sciences, Case Western Reserve University, Cleveland, OH, USA. ⁴Laboratoire PIMM, Arts et Metiers Institute of Technology, CNRS, CNAM, HESAM University, Paris, France.</footnote>

<footnote n="5">⁵Seismological Laboratory, Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, USA. ⁶e-mail: patrick.cordier@univ-lille.fr</footnote>
<table><caption>**Table 1 | Oxygen self-diffusion parameters for MgO global fit as a function of pressure and temperature (from Ita and Cohen¹⁴)**</caption><thead><tr><td>Parameter</td><td>Symbol</td><td>Value</td></tr></thead><tbody><tr><td>Attempt frequency</td><td>v<sub>a</sub></td><td>5.4 THz</td></tr><tr><td>Activation entropy at P=0</td><td>S<sub>0</sub></td><td>4 k<sub>B</sub></td></tr><tr><td>S<sub>0</sub> pressure derivative</td><td>S'<sub>0</sub></td><td>0.02 k<sub>B</sub> GPa<sup>-1</sup></td></tr><tr><td>Activation energy</td><td>E<sub>0</sub></td><td>0.4 × 10<sup>-19</sup> J</td></tr><tr><td>Activation volume at P=0</td><td>V<sub>0</sub></td><td>16.7 Å<sup>3</sup></td></tr><tr><td>V<sub>0</sub> pressure derivative</td><td>V'<sub>0</sub></td><td>-0.038 Å<sup>3</sup> GPa<sup>-1</sup></td></tr></tbody></table>

where *P* and *T* are in GPa and K respectively, *a* is the lattice constant and *k*<sub>B</sub> is the Boltzmann constant. The other parameters are listed in Table 1. The enthalpies are obtained by inverting the Arrhenius equation of diffusion coefficient:

$$D_{\text{Ox}}^{\text{sd}} = D_0 \exp\left(\frac{-\Delta H^{\text{sd}}}{k_B T}\right) \quad (2)$$

where $D_0 = 7.2 \times 10^{-4} \text{ m}^2 \text{ s}^{-1}$ and $\Delta H^{\text{sd}}$ are, respectively, the exponential prefactor and the activation enthalpy¹¹. The activation enthalpies for each set of pressure/temperature and the diffusion coefficients are given in Table 2.

This expression produces values of $D_{\text{Ox}}^{\text{sd}}$ that are in agreement with the most reliable experimental results on intrinsic O diffusion from Yang and Flynn¹⁵ (Fig. 1). Higher values of $D_{\text{Ox}}^{\text{sd}}$ determined by Van Orman et al.¹⁶ at high pressure are considered to reflect the contribution of dislocations¹⁷ and/or other extended defects, as in most other studies of oxygen diffusion in MgO¹³. The results of Ita and Cohen¹⁴ are also similar to those presented in several later theoretical studies¹⁸⁻²⁰.

## Modelling dislocation creep of periclase

To model plasticity of periclase under mantle conditions we use a 2.5D DD simulation approach already benchmarked for periclase at ambient

<table><caption>**Table 2 | Oxygen self-diffusion coefficients (m²s⁻¹) from equation (1) and activation enthalpy (eV) for MgO obtained from equation (2)**</caption><thead><tr><td>P (GPa)</td><td>T (K)</td><td>D<sup>sd</sup><sub>Ox</sub></td><td>ΔH</td></tr></thead><tbody><tr><td>30</td><td>2,000</td><td>6.61×10<sup>-27</sup></td><td>9.15</td></tr><tr><td>60</td><td>2,300</td><td>3.05×10<sup>-29</sup></td><td>11.59</td></tr><tr><td>90</td><td>2,500</td><td>3.46×10<sup>-31</sup></td><td>13.56</td></tr><tr><td>120</td><td>2,800</td><td>5.26×10<sup>-31</sup></td><td>15.09</td></tr></tbody></table>

pressure¹¹. DD simulation is a modelling tool that describes the collective motion and interactions of dislocations at the mesoscale. This technique is based on continuum elasticity theory, which provides the description of the elastic field induced by dislocations in a crystal, their interactions with each other and with respect to the stress field resulting from external loading. Within our DD code, dislocation glide is coupled with climb to investigate high-temperature creep (Methods). Our creep results are presented in Fig. 2a, which shows the creep strain rates obtained from the 2.5D DD simulations as a function of the creep stress for the four (*P*, *T*) conditions selected. A constant slope is found for all these strain rate-stress curves, which gives the value of the power law stress exponent, *n* = 3.1. A strong pressure dependence is found between (30 GPa, 2,000 K) and (90 GPa, 2,500 K), whereas creep rates at (90 GPa, 2,500 K) and (120 GPa, 2,800 K) are very close due to very similar values for the diffusion coefficients (Fig. 1 and Table 2). On Fig. 2a we compare the creep rates of periclase with those of bridgmanite at the same depths. For bridgmanite, we use the results of Reali et al.²¹ for pure climb creep (PCC) with a vacancy concentration $X_v = 1 \times 10^{-5}$; these results were shown to be in good agreement with empirical geophysical constraints on the viscosity profile of Earth's lower mantle. Bridgmanite and periclase exhibit similar stress exponents. However, we see that these calculations lead to bridgmanite having higher creep rates than periclase at any depth. To better illustrate this difference, we report, on Fig. 2b, the viscosity contrast between those two phases, taken as the ratio $\eta_{\text{per}}/\eta_{\text{Bdm}}$. We see that under lower mantle conditions bridgmanite creep rates are $10^3$ to $10^7$.

<figure><img src="image_3.png" /><figcaption>**Fig. 1 | Oxygen diffusion coefficients in periclase.** The solid black line represents the diffusion coefficients of Yang and Flynn¹⁵. The black dashed line corresponds to the diffusion coefficients of Ita and Cohen¹⁴ at 0 GPa. The other lines (blue, green, red and magenta) correspond to equation (1) at 30, 60, 90</figcaption></figure>

and 120 GPa respectively. The range of temperatures considered in this study [2,000 K:2,800 K] is shown with the grey shading. The symbols (filled circles) correspond to the (*P*, *T*) conditions along the geotherm considered in the present study¹² (Table 2).
<figure><img src="image_4.png" /><figcaption>**Fig. 2** | Creep rate modelling. **a**, Creep strain rates as a function of applied stress of bridgmanite (Bdm; solid line) and periclase (Per; symbols) at 30 GPa (blue), 60 GPa (green), 90 GPa (red) and 120 GPa (magenta). **b**, Viscosity</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>contrast $\eta_{Per}/\eta_{Bdm}$ (calculated at 10 MPa) between the two pure phases, periclase and bridgmanite, as a function of depth.</figcaption></figure>

greater than periclase creep rates. This result reflects the fact that in the lower mantle, anionic diffusion in periclase is much slower than cationic diffusion in bridgmanite (Extended Data Fig. 2). In the online Methods section and Extended Data Figs. 1–4 we consider other diffusion data from the literature for periclase and other vacancy concentrations in bridgmanite and show that the conclusions stated above remain valid in any case.

## Laboratory versus nature

In this paragraph we address the influence of strain rate on deformation mechanisms. Our calculations showing that bridgmanite deforms faster than periclase under mantle conditions seem to contradict the

<figure><img src="image_6.png" /><figcaption>**Fig. 3** | Influence of strain rate on deformation mechanisms. This figure corresponds to approximately 30 GPa, 2,000 K, where the experimental (exp.) data of Girard et al.³ (symbols) can be compared to numerical models (lines). At laboratory stresses and strain rates numerical models based on dislocation glide are in fair agreement with experimental data and show that bridgmanite (Bdm) is stronger than periclase (Per). At mantle strain rates, creep (Weertman creep modelled by 2.5D DD in periclase (see text) and pure climb creep (PCC) in bridgmanite) is controlled by ionic diffusion, which leads to an inversion of rheology, with periclase deforming more slowly than bridgmanite.</figcaption></figure>

common intuition that the silicate is stronger than the oxide. Beyond general considerations¹ or studies of analogous compounds²,²³, Girard et al.³ conducted experiments under mantle conditions of pressure and temperature that show very clearly that ferropericlase accommodates most of the strain compared to bridgmanite. We show here that there is no fundamental opposition between modelling and experiments and that our models are able to account equally for natural and experimental conditions. The experiments of Girard et al.³ are conducted at a strain rate of $3 \times 10^{-5} \text{ s}^{-1}$ and under very high stresses. Under these laboratory conditions, dislocation glide is expected as the dominant deformation mode. For bridgmanite, dislocation glide is thermally activated under these conditions and a model based on the Orowan equation can be built²⁴. The stress values predicted by this model are reported on Fig. 3. They are close to the experimental values (slightly lower as this model was constructed for single crystals and does not involve grain boundaries). At 30 GPa and 2,000 K, dislocations in periclase glide in the athermal regime⁶. Several deformation experiments of periclase report that in the athermal regime, no steady state is reached and that strain hardening occurs²⁵⁻²⁷ under the influence of dislocation–dislocation interactions. The origin of this strain hardening has been studied and modelled by Amodeo et al.¹⁰. It is well described by the Taylor mechanism, and the athermal stress $\tau_f$ required to overcome the forest dislocation network is given by:

$$ \tau_f = \mu b \frac{\ln(1/b \sqrt{\beta \rho_f})}{\ln(1/b \sqrt{\beta \rho_{\text{ref}}})} \sqrt{\beta \rho_f} \quad (3) $$

where $\mu$ is the shear modulus, $b$ the modulus of the Burgers vector, $\beta$ is a forest strengthening coefficient calibrated by Amodeo et al.¹⁰, $\rho_f$ is the forest dislocation density and $\rho_{\text{ref}}$ is the reference dislocation density. Transmission electron microscopy observations²⁸ on the specimens of Girard et al.³ confirm this view and show a very large dislocation density estimated at over $10^{15} \text{ m}^{-2}$. Using this dislocation density and equation (3) leads to an excellent description of experimental values for the deformation of ferropericlase under shallow lower mantle $P, T$ conditions at laboratory strain rates (Fig. 3). Our modelling approach allows us to reproduce very well the experimental data of Girard et al.³ and confirm that, under laboratory conditions, bridgmanite is indeed expected to be stronger than ferropericlase. It should be noted that the flow stress value for ferropericlase corresponds to a very high hardening and that under more ‘gentle’ laboratory
<figure><img src="image_7.png" /><figcaption>**Fig. 4 | Rheology of bridgmanite-periclase aggregates.** Pressure dependent flow stress $\hat{\sigma}_0$ for a bridgmanite (Bdm)–periclase (Per) aggregate deformed in creep at $10^{-15} \text{ s}^{-1}$. Calculations are performed within the POSO-SC scheme, with a volume fraction of periclase of 0.25 and a vacancy concentration in bridgmanite of $X_v = 1 \times 10^{-5}$. The ratio of strain rate (periclase/bridgmanite) is also indicated. Flow stresses $\sigma_0$ of pure bridgmanite (in red) and pure periclase (in black) are also shown for comparison.</figcaption></figure>

deformation conditions (that is, smaller strains or more recovery), the rheology contrast under laboratory conditions should be even larger. At 30 GPa, oxygen diffusion is already slow, and strain rates much lower than those obtainable in laboratory experiments would be required for climb to become active and promote the recovery necessary for steady-state creep.

## Deformation of a bridgmanite-periclase aggregate

Our model places first-order constraints on the intrinsic properties of the major phases of the lower mantle. The rheological behaviour of rocks in the lower mantle certainly involves more complexity. A first step towards this direction is to use dedicated micromechanical models, to connect the viscoplastic behaviour of a polycrystalline aggregate representative of a lower mantle assemblage with those of the individual phases, bridgmanite and periclase, and, at the same time, quantify the contribution of both phases to the effective strain rate, in relation with the aggregate microstructure.

We have therefore applied a state-of-the-art scale transition homogenization method²⁹ (partially optimized second-order self-consistent procedure, denoted POSO-SC) adapted for multiphase aggregates with nonlinear rheologies, such as those considered here. The model considers an isotropic microstructure, that is, a random mixture of bridgmanite and periclase with no preferred elongated shape nor crystallographic texture for each phase. Details are given in the Methods. The effective viscoplastic behaviour can be expressed with a nonlinear creep law with a stress sensitivity $n = 3.1$ (see Methods). Figure 4 shows the results obtained for a vacancy concentration $X_v = 1 \times 10^{-5}$ in bridgmanite and a volume fraction of periclase of 25% (alternative figures considering different combinations of diffusion coefficients are presented in Extended Data Fig. 5). For pressures in the range 30–120 GPa, the effective rheology (blue line) lies very close to that of pure bridgmanite (red line), with the flow stress values differing by a factor of only around 1.65, whereas the flow stress of pure periclase (black line) is one to two orders of magnitude larger. Periclase deforms at least 100 times slower than bridgmanite in the aggregate and its contribution to the effective strain rate is negligible. In the online Methods section, we show that the volume fraction of periclase would have to exceed 40% to significantly affect the rheology of the

assemblage. The presence of iron in MgO may affect anionic diffusivity but does not change our conclusions for the bulk lower mantle (see online Methods section).

Following the widespread acceptance that ferropericlase is the weakest phase¹, and inspired by the experimental results of Girard et al.³, several studies³⁰⁻³² have discussed the implication of the morphology of ferropericlase (either forming a load-bearing framework or an interconnected weak layer³³) on the rheology of the lower mantle. Our study shows that under the conditions (mainly strain rate) of the lower mantle, the behaviour of periclase is radically different and that this fundamental assumption is now challenged. Periclase, minor in volume, hardly deforms and probably induces little evolution of the microstructure. Consequently, the rheology of the lower mantle is very well described by that of bridgmanite deforming by pure climb creep.

## Online content

Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41586-022-05410-9.

1. Yamazaki, D. & Karato, S.-I. Some mineral physics constraints on the rheology and geothermal structure of Earth's lower mantle. *Am. Mineral.* **86**, 385-391 (2001).

2. Miyagi, L. & Wenk, H.-R. Texture development and slip systems in bridgmanite and bridgmanite + ferropericlase aggregates. *Phys. Chem. Minerals* **43**, 597-613 (2016).

3. Girard, J., Amulele, G., Farla, R., Mohiuddin, A. & Karato, S. I. Shear deformation of bridgmanite and magnesiowüstite aggregates at lower mantle conditions. *Science* **351**, 144-147 (2016).

4. Amodeo, J., Carrez, P. H., Devincre, B. & Cordier, P. Multiscale modelling of MgO plasticity. *Acta Mat.* **59**, 2291-2301 (2011).

5. Cordier, P., Amodeo, J. & Carrez, P. H. Modelling the rheology of MgO under Earth's mantle pressure, temperature and strain-rates. *Nature* **481**, 177-180 (2012).

6. Amodeo, J., Carrez, P. H. & Cordier, P. Modelling the effect of pressure on the critical shear stress of MgO single crystals. *Philos. Mag.* **92**, 1523-1541 (2012).

7. Girard, J., Chen, J. & Raterron, P. Deformation of periclase single crystals at high pressure and temperature: quantification of the effect of pressure on slip-system activities. *J. Appl. Phys.* **111**, 112607 (2012).

8. Marquardt, H. & Miyagi, L. Slab stagnation in the shallow lower mantle linked to an increase in mantle viscosity. *Nature Geosci.* **8**, 311-314 (2015).

9. Immoor, J. et al. Evidence for (100)\<011\> slip in ferropericlase in Earth's lower mantle from high-pressure/high-temperature experiments. *Earth Planet. Sci. Lett.* **489**, 251-257 (2018).

10. Amodeo, J., Devincre, B., Carrez, P. H. & Cordier, P. Dislocation reactions, plastic anisotropy and forest strengthening in MgO at high temperature. *Mech. Materials* **71**, 62-73 (2014).

11. Reali, R. et al. Modeling plasticity of MgO by 2.5D dislocation dynamics simulations. *Mat. Science Eng.* **A 690**, 52-61 (2017).

12. Stacey, F. & Davis, P. High pressure equations of state with applications to the lower mantle and core. *Phys. Earth Planet. Int.* **142**, 137-184 (2004).

13. Van Orman, J. & Crispin, K. Diffusion in oxides. *Rev. Mineral. Geochem.* **72**, 757-825 (2010).

14. Ita, J. & Cohen, R. E. Effects of pressure on diffusion and vacancy formation in MgO from non-empirical free-energy integrations. *Phys. Rev. Lett.* **79**, 3198 (1997).

15. Yang, M. & Flynn, C. Intrinsic diffusion properties of an oxide: MgO. *Phys. Rev. Lett.* **73**, 1809-1812 (1994).

16. Van Orman, J. A., Fei, Y., Hauri, E. H. & Wang, J. Diffusion in MgO at high pressures: constraints on deformation mechanisms and chemical transport at the core-mantle boundary. *Geophys. Res. Lett.* **30**, 1056-1059 (2003).

17. Landeiro Dos Reis, M., Giret, Y., Carrez, P. & Cordier, P. Efficiency of the vacancy pipe diffusion along an edge dislocation in MgO. *Comput. Mat. Sci.* **211**, 111490 (2022).

18. Ammann, M. W., Brodholt, J. P. & Dobson, D. P. Simulating diffusion. *Rev. Mineral. Geochem.* **71**, 201-224 (2010).

19. Karki, B. B. & Khanduja, G. Vacancy defects in MgO at high pressure. *Am. Mineral.* **91**, 511-516 (2006).

20. Runevall, O. & Sandberg, N. Self-diffusion in MgO—a density functional study. *J. Phys. Condens. Matter* **23**, 345402 (2011).

21. Reali, R. et al. The role of diffusion-driven pure climb creep on the rheology of bridgmanite under lower mantle conditions. *Sci Rep.* **9**, 2053 (2019).

22. Wang, Y. et al. High pressure, high-temperature deformation of CaGeO₃ (perovskite) ±MgO aggregates: implications for multiphase rheology of the lower mantle. *Geochem. Geophys. Geosyst.* **14**, 3389-3408 (2013).

23. Kaercher, P. et al. Two-phase deformation of lower mantle mineral analogs. *Earth Planet. Sci. Lett.* **456**, 134-145 (2016).

24. Kraych, A., Carrez P. H. & Cordier, P. On dislocation glide in MgSiO₃ bridgmanite at high-pressure and high-temperature. *Earth Planet. Sci. Lett.* **452**, 60-68 (2016).

25. Copley, S. & Pask, J. Plastic deformation of MgO single crystals up to 1600 C. *J. Am. Ceram. Soc.* **48**, 139-146 (1965).
26. Hulse, C. & Pask, J. Mechanical properties of magnesia single crystals in compression. *J. Am. Ceram. Soc.* **43**, 373–378 (1960).

27. Routbort, J. L. Work-hardening and creep of MgO. *Acta Metall.* **27**, 649–661 (1979).

28. Nzogang, B. et al. Characterization by scanning precession electron diffraction of an aggregate of bridgmanite and ferropericlase deformed at HP-HT. *Geochem. Geophys. Geosyst.* **B** **19**, 582–594 (2018).

29. Ponte Castañeda, P. Second-order homogenization estimates for nonlinear composites incorporating field fluctuations: I—theory. *J. Mech. Phys. Solids* **50**, 737–757 (2002).

30. Madi, K., Forest, S., Cordier, P. & Boussuge, M. Numerical study of creep in two-phase aggregates with a large viscosity contrast: implications for the rheology of the lower mantle. *Earth Planet. Sci. Lett.* **237**, 223–238 (2005).

31. Kasemer, M., Zepeda-Alarcon, E., Carson, R., Dawson, P. & Wenk, H.-R. Deformation heterogeneity and intragrain lattice misorientation in high strength contrast, dual-phase bridgmanite/periclase. *Acta Mat.* **189**, 284–298 (2020).

32. Thielmann, M., Golabek, G. J. & Marquardt, H. Ferropericlase control of lower mantle rheology: impact of phase morphology. *Geochem. Geophys. Geosyst.* **21**, e2019GC008688 (2020).

33. Handy, M. R. Flow laws for rocks containing two nonlinear viscous phases: a phenomenological approach. *J. Struct. Geol.* **16**, 287–301 (1994).

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

<figure><img src="image_8.png" /></figure>

**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/.

© The Author(s) 2023
## Dislocation dynamics calculations

The 2.5D DD simulation method is used throughout this study to model the creep properties of periclase under lower mantle conditions under a constant creep stress ranging from 10 to 150 MPa. We consider the *P*, *T* conditions on the geotherm of Stacey and Davis¹² at 30, 60, 90 and 120 GPa, and use the periclase elastic properties from Karki et al.³⁴ (Extended Data Table 2).

To capture the statistically representative volume element of the microstructure, dislocations of Burgers vector **b** = 1/2<110> are introduced edge-on in two intersecting slip systems within a square element of typical cell size $L_x = L_y$ between 8 to 60 µm (function of the desired dislocation densities) and periodic boundary conditions are applied. As extensively described in some of our previous works¹¹, specific characteristics of a realistic MgO microstructure are introduced as local rules, including junction formation strengthening and a dislocation multiplication rate to allow random introduction of fresh dislocations of opposite signs on all slip systems. All local rules are parameterized according to experimental observations or three-dimensional DD simulations as described in Reali et al.¹¹.

For a realistic simulation of the creep behaviour, we use a combination of glide and climb events to control the displacement of dislocations in the simulation volume. The Burgers vector contained in the reference plane defines the slip direction, whereas the climb direction is defined as the one in the reference plane, orthogonal to the Burgers vector. The positive orientation of the climb direction is taken along the vacancy emission direction. During a given timestep of the simulation, dislocations are displaced along the glide or climb direction according to the velocity formulation presented below.

In our simulations, dislocations glide in the athermal regime according to a «free-flight» velocity $v_g$, which is a function of the effective resolved shear stress $\tau = \tau_{app} + \tau_{int}$ (that is, including the applied resolved shear stress and the elastic interaction stress):

$$v_g = \frac{b\tau}{B}$$

As, in the range of *P* and *T* conditions investigated here, the lattice friction of MgO vanishes⁵,⁶, the glide velocity is thus expressed with a viscous drag coefficient *B*. Such a coefficient is known to be temperature and material dependent. Thus, we consider a linear scaling of *B* with temperature³⁵ and use as a starting point the experimental value recorded at 300 K³⁶.

Climb events occur through the self-diffusion of vacancies adsorbed or emitted by the dislocation line. To describe the climb rate in steady-state creep equilibrium conditions, we use a dislocation climb velocity $v_c$ expressed as follows³⁵:

$$v_c = \alpha \frac{D^{\text{sd}}}{b} \left[ \exp\left(\frac{\tau_c \Omega}{k_B T}\right) - \frac{c_\infty}{c_0} \right]$$

where $\tau_c$ is the climb stress, calculated similarly to $\tau$ but resolved here along the climb direction. $\alpha = 2\pi/\ln(R/r_c)$ is a geometrical factor that describes the cylindrical geometry of the vacancy flux field around the dislocation line, where *R* and $r_c$ represent the two radii of the cylindrical surfaces through which the vacancy flux is calculated. Again, in the range of temperature investigated, the dislocation line must contain a high density of jogs, and thus be saturated with jogs, enabling vacancies to be absorbed or emitted instantaneously by the dislocation. Here, $r_c$ is the dislocation core capture radius and *R* is taken as a fraction of the average dislocation distance. Being within the logarithmic term, the $R/r_c$ ratio does not significantly affect the climb velocity values and here is taken as constant and equal to 100³⁷. $\Omega$ is the formation volume of vacancies in the Mg and O sites and is calculated from the unit cell

volume of MgO—that is, $\Omega = a^3/Z$ where $Z = 4$ is the number of formula units per unit cell. $c_m$ and $c_0$ are the vacancy concentrations far from the dislocation and the equilibrium vacancy concentration in the bulk volume, respectively. Following Boioli et al.³⁷, we assume that far away from the dislocations (that is, on the external limit *R* of the cylindrical flux) the vacancy concentration is constant and equal to the equilibrium concentration in the bulk volume ($c_m = c_0$). Finally, $D_{sd}$ is the vacancy self-diffusion coefficient, which controls the flow of atomic species from and to the climbing dislocation.

## Micromechanical model of the rheology of the aggregate

**Local behaviour.** We apply state-of-the-art formulations of nonlinear homogenization methods to estimate and bound the effective viscosity of a two-phase material constituted by bridgmanite and periclase with various volume fractions, at pressures and temperatures relevant for the lowermost part of Earth's mantle. These phases are supposed to be randomly mixed so that the effective mechanical behaviour of the whole aggregate can be considered isotropic. Bridgmanite and periclase are considered to be two isotropic viscous phases too. In reality, they are both made of grains, but we consider that these grains are equiaxed and randomly oriented so that, on average, the polycrystalline aggregates can be replaced by two homogenous and isotropic phases. The role of grain or phase boundaries is not explicitly taken into account. We are interested in the effective rheological behaviour in the permanent (that is, secondary) creep regime. Thus, the elastic behaviour of both phases does not come into play. At the scale of each phase, the viscoplastic behaviour can be described by the following nonlinear creep law:

$$\dot{\epsilon}_{eq} = \dot{\epsilon}_0 \left( \frac{\sigma_{eq}}{\sigma_0} \right)^n$$

with $\dot{\epsilon}_0$ a reference strain rate, $\sigma_0$ a reference stress, *n* the stress sensitivity, $\dot{\epsilon}_{eq}$ the equivalent strain rate and $\sigma_{eq}$ the equivalent stress, defined as:

$$\dot{\epsilon}_{eq} = \sqrt{\frac{2}{3} \dot{\epsilon}_{ij} \dot{\epsilon}_{ij}} , \dot{\sigma}_{eq} = \sqrt{\frac{3}{2} \sigma'_{ij} \sigma'_{ij}}$$

with summation over repeated indices (Einstein convention), $\dot{\epsilon}_{ij}$ and $\sigma'_{ij}$ components of the strain rate tensor and deviatoric stress tensor, respectively. We chose here to take $\dot{\epsilon}_0 = 1 \times 10^{-15} \text{s}^{-1}$, that is, a value typical for lower mantle convection. Doing so, the $\sigma_0$ values are of the order of magnitude of the equivalent stress encountered in situ during mantle flow. Values of $\sigma_0$ for bridgmanite and periclase are fitted from the creep laws of Fig. 2 and available from the data deposit. The stress sensitivity *n* is the same for both phases, *n* = 3.1.

**Effective behaviour.** The effective behaviour reads in a similar way as the local behaviour:

$$\dot{\epsilon}_{eq} = \dot{\epsilon}_0 \left( \frac{\bar{\sigma}_{eq}}{\bar{\sigma}_0} \right)^{\tilde{n}}$$

where $\dot{\epsilon}$ and $\bar{\sigma}$ are the effective strain rate and stress tensors, respectively, given by the volume average (denoted $\langle \cdot \rangle$) of local strain rate and stress fields: $\dot{\epsilon} = \langle \dot{\epsilon} \rangle$ and $\bar{\sigma} = \langle \sigma \rangle$. $\tilde{n}$ is the effective stress sensitivity, which can, in general, be inferred from the homogenization procedure. In the case considered here, things are simpler as bridgmanite and periclase exhibit the same *n* value. Therefore $\tilde{n} = n = 3.1$. Finally, $\tilde{\sigma}_0$ is the effective flow stress associated with the (nonlinear) viscosity of the aggregate in situ. It is determined by the homogenization procedure detailed below.

**Homogenization procedure.** To describe the possible rheology of an assemblage representative of the lower mantle, one can reasonably
assume that bridgmanite and periclase phases are randomly mixed. In that case, a good estimation of the effective behaviour is provided by the Self-Consistent scheme (denoted SC below), based on a fully disordered microstructure³⁸,³⁹. To extend the SC scheme to nonlinear rheologies, as here, we have made use of the very efficient and accurate linearization procedure (so-called partially optimized second-order procedure (POSO)) of Ponte Castaneda²⁹.

## Sensitivity to calculation parameters

The calculations presented in the main text are made based on the diffusion coefficients for pure MgO periclase and for MgSiO₃ bridgmanite with a vacancy concentration of 10⁻⁵. For the calculation of the aggregate rheology, we considered a periclase proportion of 25% by volume. In the following, reasonable alternatives to these choices are discussed to assess the robustness of our conclusions.

### Diffusion coefficients.
As the creep rate critically depends on the diffusion coefficients, we have re-calculated the creep rates of periclase using another dataset of diffusion coefficients for oxygen in periclase, those of Yoo et al.⁴⁰, which are faster than those of Ita and Cohen¹⁴ and those of Yang and Flynn¹⁵. Fitting equation (1) on the data of Yoo et al.⁴⁰ (Extended Data Fig. 1), we determine two new constants for equation (1): $E = 7200 \times 10^{-21}$ J and $S_0 = 2k_B$. The activation enthalpies obtained with these new constants are given in Extended Data Table 1 and the evolution of the two sets of diffusion coefficients are compared in Extended Data Fig. 2. The creep rates of periclase with the diffusion coefficients of Yoo et al.⁴⁰ are then slightly faster than with the diffusion coefficients of Ita and Cohen¹⁴, but still significantly lower than those of bridgmanite (Extended Data Fig. 3). For this comparison, we have also taken into account the uncertainty on the rheology of bridgmanite owing to the lack of constraints on diffusion coefficients (this is discussed in detail by Reali et al.²¹). For this purpose, we vary the concentration of vacancies in bridgmanite $\chi_v$ between $10^{-3}$ and $10^{-6}$. Extended Data Fig. 2 shows that no matter which combination of data we may choose, Mg (or Si) diffusion in bridgmanite is always significantly faster than oxygen diffusion in periclase. The contrast is particularly marked at high pressures. Extended Data Figs. 3–7 show that even taking into account these uncertainties in the most conservative way does not change any of our conclusions.

### Influence of iron.
Our study is based on the end-member compositions MgO and MgSiO₃, whereas in the mantle these minerals contain iron (and aluminium for bridgmanite). These differences in composition are likely to affect the values of the diffusion coefficients on which our calculations are based. In bridgmanite, this is taken into account by the range of vacancy concentrations considered²¹, see discussion above and Extended Data Figs. 3, 4, 6 and 7. In ferropericlase, it has been shown that oxygen diffusion coefficients are unaffected by the presence of trivalent cations¹³,⁴¹, because these do not alter the concentration of neutral cation–anion vacancy pairs that are responsible for oxygen transport. The addition of significant amounts of FeO, however, is likely to alter oxygen diffusion coefficients. This effect was considered by Reali et al.⁴², who showed that anion diffusion in MgO and in other oxides and halides with rock salt structures is well described by a homologous temperature scaling, wherein the diffusion coefficients as a function of temperature collapse onto a common curve when normalized to the melting temperature. The addition of FeO lowers the melting temperature of periclase, and hence is expected to enhance oxygen diffusivity. We have thus calculated oxygen diffusion coefficients under lower mantle conditions for ferropericlase as a function of its Fe content, under the same assumptions made by Reali et al.⁴². Considering the addition of up to 20 mol% FeO in periclase⁴³, we do find that the enhancement of oxygen diffusion is up by a factor of approximately 2000 relative to pure MgO. However, this enhancement does not change the conclusions of this study for the bulk lower

mantle, as oxygen diffusion in ferropericlase remains significantly slower than Mg or Si in bridgmanite.

### Influence of periclase volume fraction.
We have computed the rheology of bridgmanite–periclase aggregates for various volume fractions of periclase based on four homogenization procedures (Extended Data Figs. 6 and 7), the difference between them being related to the knowledge on the specimen microstructure that is taken into account.

First, Reuss and Voigt bounds are valid for any microstructures. They only depend on the volume fraction of both phases, not on their geometrical arrangement. They provide rigorous lower and upper bounds for $\tilde{\sigma}_0$, respectively. In other words, $\tilde{\sigma}_0$ in the aggregate cannot be smaller than $\tilde{\sigma}_0^{\text{Reuss}}$ nor larger than $\tilde{\sigma}_0^{\text{Voigt}}$. The Reuss bound is obtained when assuming that both phases are submitted to the same stress, whereas the Voigt bound is obtained when assuming a uniform strain rate in both phases.

The Hashin–Strikman variational upper bound²⁹, denoted HS+ in Extended Data Figs. 6 and 7, is another rigorous upper bound for $\tilde{\sigma}_0$. It is more stringent than the Voigt upper bound as is based on an extra assumption concerning the specimen microstructure, that is, that the microstructure is isotropic. It is obtained when assuming that the softer phase (here, bridgmanite) is spread as inclusions in a stiffer matrix (periclase) while the associated nonlinear homogenization problem is addressed using the variational procedure introduced by Ponte Castaneda, also known as ‘modified secant’ method⁴⁴.

To stick more closely to the possible rheology of an assemblage representative of the lower mantle, one can reasonably assume that bridgmanite and periclase phases are randomly mixed, that is, that contrary to what is assumed within the HS+ bound, both phases are on the same footing, with neither of them playing the role of a matrix and the other one being inclusions. In that case, a good estimation of the effective behaviour is provided by the POSO extension of the Self-Consistent scheme (SC in Extended Data Figs. 6 and 7), introduced in the previous section. Here, unlike previous applications to polycrystalline aggregates, where individual crystal orientations were taken into account⁴⁵,⁴⁶, a two-phase aggregate was considered. A fully optimized linearization procedure for the SC scheme has been proposed recently⁴⁷,⁴⁸ and applied to minerals of the mantle transition zone⁴⁹; results are close to the POSO version when the mechanical contrast between the phases and the nonlinearity of their behaviour is not too large, as here.

Typical results are shown in Extended Data Figs. 6 and 7, calculated here for pressures between 30 and 120 GPa and $\chi_v = 10^{-3}-10^{-6}$. Due to the large mechanical contrast between periclase and bridgmanite, the effective rheology of the aggregate, predicted by the POSO-SC model, lies very close to the Reuss lower bound for volume fractions lower than 30%, whereas it is close to the Voigt upper bound for a volume fraction larger than 70%. In between, the effective stress $\tilde{\sigma}_0$ evolves significantly, by more than one order of magnitude. This would correspond to a change in viscosity of more than three orders of magnitude, the viscosity being proportional to $(\tilde{\sigma}_0)^n$. One other specific feature of these results is illustrated in the right panels of Extended Data Figs. 6 and 7, which show the mean equivalent strain rate in the two phases within the aggregate, normalized by the effective strain rate, as predicted by the SC model. It is obtained that the periclase strain rate is really small at volume fractions lower than 40%. This volume fraction corresponds to a ‘mechanical percolation threshold’ observed in previous work⁵⁰,⁵¹. At small volume fractions of the highly viscous periclase phase, the aggregate can deform with almost no plastic strain in that phase.

## Data availability

The data of this manuscript are available at https://doi.org/10.5281/zenodo.5970733. Source data are provided with this paper.
Code availability

The source code of the 2.5D DD is available at https://doi.org/10.5281/zenodo.6306546.

34. Karki, B. B. et al. Structure and elasticity of MgO at high pressure. Am. Mineral. 82, 51–60 (1997).

35. Hirth, J. P. & Lothe, J. *Theory of Dislocations* (John Wiley & Sons, 1992).

36. Singh, R. N. & Coble, R. L. Dynamic dislocation behavior in “pure” magnesium oxide single crystals. J. Appl. Phys. 45, 981–989 (1974).

37. Boioli, F., Carrez, P., Cordier, P., Devincre, B. & Marquille, M. Modeling the creep properties of olivine by 2.5-dimensional dislocation dynamics simulations. Phys. Rev. B 92, 014115 (2015).

38. Hershey, A. V. The elasticity of an isotropic aggregate of anisotropic cubic crystals. J. Appl. Mech. 21, 236–240 (1954).

39. Kröner, E. Self-consistent scheme and graded disorder in polycrystal elasticity. J. Phys. F: Metal Phys. 8, 2261–2267 (1978).

40. Yoo, H., Wuensch, B. J. & Petuskey, W. T. Oxygen self-diffusion in single crystal MgO: secondary-ion mass spectrometric analysis with comparison of results from gas-solid and solid-solid exchange. Solid State Ion. 150, 207–221 (2002).

41. Ando, K., Kurokawa, Y. & Oishi, Y. Oxygen self-diffusion in Fe-doped MgO single crystals. J. Chem. Phys. 78, 6890–6892 (1983).

42. Reali, R. et al. Modeling viscosity of (Mg,Fe)O at lowermost mantle conditions. Phys. Earth Planet. Int. 287, 65–75 (2019).

43. Xu, S., Lin, J.-F. & Morgan, D. Iron partitioning between ferropericlase and bridgmanite in the Earth's lower mantle. J. Geophys. Res. Solid Earth 122, 1074–1087 (2017).

44. Suquet, P. Overall properties of nonlinear composites: a modified secant moduli theory and its link with Ponte Castañeda nonlinear variational procedure. C. R. Acad. Sci., Paris—Series IIB 320, 563–571 (1995).

45. Castelnau, O., Blackman, D. K., Lebensohn, R. A. & Ponte Castañeda, P. Micromechanical modeling of the viscoplastic behavior of olivine. J. Geophys. Res. 113, B09202 (2008).

46. Raterron, P. et al. Multiscale modeling of upper mantle plasticity: from single-crystal rheology to multiphase aggregate deformation. Phys. Earth Planet. Int. 228, 232–243 (2014).

47. Ponte Castañeda, P. Fully optimized second-order variational estimates for the macroscopic response and field statistics in viscoplastic crystalline composites. Proc. Royal Soc. A. 471, 20150665 (2015).

48. Song, D. & Ponte Castañeda, P. Fully optimized second-order homogenization estimates for the macroscopic response and texture evolution of low-symmetry viscoplastic polycrystals. Int. J. Plasticity 110, 272–293 (2018).

49. Castelnau, O. et al. Multiscale modeling of the effective viscoplastic behavior of Mg₃SiO₄ wadsleyite: bridging atomic and polycrystal scales. Comptes Rendus Mécanique. 348, 827–846 (2020).

50. Budiansky, B. On the elastic moduli of some heterogeneous materials. J. Mech. Phys. Solids 13, 223–227 (1965).

51. Gilormini, P., Brenner, R. & Castelnau, O. A similarity between the classical and modified secant extensions of the self-consistent model. C. R. Acad. Sci., Paris—Series IIB 329, 523–527 (2001).

**Acknowledgements** We acknowledge funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme under grant agreement no. 787198—TimeMan and from the National Science Foundation's (NSF) Collaborative Study of Earth's Deep Interior under EAR-1161046 and EAR-2009935. We are thankful to the CIDER programme set at the Kavli Institute for Theoretical Physics, University of California, Santa Barbara (NSF EAR-1135452, funded under the FESD Program).

**Author contributions** The study was designed by P. Cordier, K.G. and T.W. performed the dislocation dynamics calculations. O.C. performed the micromechanical (self-consistent) calculations. P. Cordier, K.G., T.W., J.V.O., O.C., J.M.J. and P. Carrez discussed and analysed the data. P. Cordier wrote the paper with contributions from all authors.

**Competing interests** The authors declare no competing interests.

**Additional information**

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41586-022-05410-9.

**Correspondence and requests for materials** should be addressed to Patrick Cordier.

**Peer review information** Nature thanks Ian Jackson and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

**Reprints and permissions information** is available at http://www.nature.com/reprints.
<figure><img src="image_9.png" /></figure>

**Extended Data Fig. 1 |** Oxygen diffusion coefficients in periclase. The black dotted line corresponds to the fit of Ita and Cohen¹⁴ diffusion law at 0 GPa on the experimental data of Yoo et al.⁴⁰. Comparison with results of Yang and Flynn¹⁵ results is provided. The other lines (blue, green, red and magenta)

correspond to this new set at 30, 60, 90 and 120 GPa respectively. The square symbols correspond to the diffusion coefficients along the geotherm deduced from this data. Comparison with data from Ita and Cohen¹⁴ used in the main text (filled circles symbols) is provided.
<figure><img src="image_10.png" /></figure>

**Extended Data Fig. 2 | Pressure dependence of ionic diffusion.** The lines represent oxygen diffusion coefficients in periclase (Per). We compare the data of Ita and Cohen¹⁴ (solid lines) to those obtained from the pressure dependence

of Ita and Cohen¹⁴ but anchored on data from Yoo et al.⁴⁰ (dotted lines). The symbols correspond to magnesium diffusion in bridgmanite (Bdm) considering a range of vacancy concentrations from 10⁻³ to 10⁻⁶.
<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

**Extended Data Fig. 3** | Creep laws in bridgmanite (Bdm) and periclase (Per).
In this figure we illustrate the influence of possible uncertainties on the
diffusion coefficients in bridgmanite and periclase. For bridgmanite (Bdm)

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

vacancy concentrations of 10⁻³ and 10⁻⁶ are considered. Creep in periclase is
based on diffusion data of Ita and Cohen¹⁴ and Yoo et al.⁴⁰. Whatever the set of
data considered, bridgmanite always deforms faster than periclase.
<figure><img src="image_15.png" /></figure>

**Extended Data Fig. 4 | Viscosity contrast between bridgmanite (Bdm) and periclase (Per).** In this figure we illustrate the influence of possible uncertainties on the diffusion coefficients in bridgmanite and periclase. For bridgmanite (Bdm) vacancy concentrations in the range 10⁻³ to 10⁻⁶ are

considered. Creep in periclase is based on diffusion data of Ita and Cohen¹⁴ (Solid lines) and Yoo et al.⁴⁰ (dotted lines). Whatever the set of data considered, bridgmanite always deforms faster than periclase.
<figure><img src="image_16.png" /></figure>

**Extended Data Fig. 5** | Pressure dependence of the aggregate rheology. The flow stress $\delta_0$ (blue line) for a bridgmanite (Bdm) – periclase (Per) aggregate (volume fraction of MgO: 0.25), is obtained with the POSO Self-Consistent scheme, for vacancy concentrations of 10$^{-3}$ and 10$^{-6}$ in bridgmanite and two

oxygen diffusion coefficients for periclase (from Ita and Cohen¹⁴ and Yoo et al.⁴⁰). Flow stresses of pure bridgmanite (red line) and pure periclase (black line) are also shown for comparison.
<figure><img src="image_17.png" /><figcaption>**Extended Data Fig. 6** | Influence of periclase volume fraction ($X_v = 10^{-3}$). Behavior for a vacancy concentration in bridgmanite of $X_v = 10^{-3}$ and oxygen diffusion coefficient in periclase from Ita & Cohen¹⁴ at various pressures.</figcaption></figure>

Left: Effective flow stress $\dot{\sigma}_f$ based on four homogenization procedures: Reuss, Voigt, Hashin-Strikman (HS+) and the POSO extension of the Self-Consistent scheme (SC). Right: Relative strain-rates in bridgmanite and periclase.
<figure><img src="image_18.png" /><figcaption>**Extended Data Fig. 7** | Influence of periclase volume fraction ($X_v = 10^{-6}$). Behavior for a vacancy concentration in bridgmanite of $X_v = 10^{-6}$ and oxygen diffusion coefficient in periclase from Ita & Cohen¹⁴ at various pressures.</figcaption></figure>

Left: Effective flow stress $\dot{\sigma}_0$ based on four homogenization procedures: Reuss, Voigt, Hashin-Strikman (HS+) and the POSO extension of the Self-Consistent scheme (SC). Right: Relative strain-rates in bridgmanite and periclase.
Extended Data Table 1 | Additional diffusion coefficients

<table><thead><tr><th><i>P</i> (GPa)</th><th><i>T</i> (K)</th><th><i>D</i><sub>Ox</sub><sup>sd</sup></th><th>Δ<i>H</i></th></tr></thead><tbody><tr><td>30</td><td>2000</td><td>1.25×10<sup>-24</sup></td><td>8.25</td></tr><tr><td>60</td><td>2300</td><td>2.50×10<sup>-27</sup></td><td>10.74</td></tr><tr><td>90</td><td>2500</td><td>1.54×10<sup>-29</sup></td><td>12.74</td></tr><tr><td>120</td><td>2800</td><td>1.26×10<sup>-29</sup></td><td>14.32</td></tr></tbody></table>

Diffusion coefficients (m²/s) from eq. 1 fit on Yoo et al.⁴⁰ and activation enthalpy (eV) for MgO obtained from eq. 2.
Extended Data Table 2 | Periclase parameters

<table><thead><tr><th></th><th>30 GPa,<br>2,000 K</th><th>60 GPa,<br>2,300 K</th><th>90 GPa,<br>2,500 K</th><th>120 GPa,<br>2,800 K</th></tr></thead><tbody><tr><td>Lattice constant a (Å)</td><td>4.10</td><td>3.95</td><td>3.82</td><td>3.75</td></tr><tr><td>Burgers vector b (Å)</td><td>2.90</td><td>2.80</td><td>2.70</td><td>2.65</td></tr><tr><td>Poisson ratio</td><td>0.23</td><td>0.26</td><td>0.28</td><td>0.30</td></tr><tr><td>Shear Modulus μ (GPa)</td><td>176</td><td>218</td><td>252</td><td>280</td></tr><tr><td>Bulk Modulus (GPa)</td><td>275</td><td>387</td><td>498</td><td>600</td></tr></tbody></table>

Lattice parameters and elastic properties of periclase used in 2.5D Dislocation Dynamics.