# Prediction of mechanical twinning in magnesium silicate post-perovskite

Philippe Carrez, Alexandra M. Goryaeva, Patrick Cordier

► To cite this version:

Philippe Carrez, Alexandra M. Goryaeva, Patrick Cordier. Prediction of mechanical twinning in magnesium silicate post-perovskite. Scientific Reports, 2017, 7, pp.1-9. 10.1038/s41598-017-18018-1. hal-02628924

HAL Id: hal-02628924

https://hal.inrae.fr/hal-02628924

Submitted on 27 May 2020

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
SCIENTIFIC REPORTS

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

OPEN

# Prediction of Mechanical Twinning in Magnesium Silicate Post-Perovskite

Philippe Carrez, Alexandra M. Goryaeva & Patrick Cordier

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

Received: 18 September 2017

Accepted: 5 December 2017

Published online: 15 December 2017

Seismic anisotropy is one of our major sources of information about the dynamic processes and flow in the Earth's mantle. In contrast to the bulk of the lower mantle, which appears to be mostly isotropic, the lowermost lower mantle exhibits strong seismic anisotropy and major heterogeneities. In particular, distinct anisotropy signatures are found in regions thought to be associated with cold downwelling<sup>1,2</sup>. The discovery in 2004 that bridgmanite, the magnesium silicate with perovskite structure which is the main constituent of the lower mantle, is not stable at pressures comparable to those of the "D" layer and transforms at c.a. 120 GPa into a distinct structurally anisotropic phase, named post-perovskite<sup>3</sup>, has attracted considerable attention. The potential link between this new phase and the seismic anisotropy in the lowermost lower mantle received more support when it was established that the post-perovskite phase can form only in a relatively cold mantle<sup>4</sup>. The high-pressure magnesium silicate post-perovskite phase exhibits an orthorhombic layered crystal structure of space group *Cmcm* with strongly different lattice parameters<sup>3</sup>. Exhibiting layers of SiO<sub>6</sub> octahedrons parallel to {010}, the structure is thus highly anisotropic with such a structural characteristic being potentially related to the strong seismic anisotropy of "D". To further establish the role of post-perovskite and to ultimately decipher the flow patterns at the base of the mantle, it is necessary to understand how crystal preferred orientations (CPOs) develop in this phase during plastic flow<sup>5</sup>. Given the very high-pressure required to stabilize the magnesium silicate post-perovskite, only a few set of experiments have been conducted on this phase<sup>6-8</sup> and most experiments have been performed on analogue materials with the same crystal structure<sup>9-11</sup>, but stable at lower pressures. This includes calcium iridate (CaIrO<sub>3</sub>), which is stable at ambient pressure<sup>12-15</sup>. Unfortunately, all these experiments have led to conflicting results, possibly because of textures inherited from phase transformations<sup>11,12,16</sup> and differences in the crystal chemistry of the analogue materials<sup>8,17,18</sup>. Given the formidable difficulties of deformation experiments under very high pressures, numerical modelling currently represents a very attractive alternative. Using atomic-scale modelling of dislocations<sup>19-21</sup>, we have shown that shearing the post-perovskite structure occurs easily along the shortest [100] lattice repeat in the Mg-O layer (010) plane, with a lattice friction of 2 GPa<sup>20</sup>. The other dense direction in this plane, [001], is the second easiest<sup>21</sup> (ca. 3 GPa). However, shearing the Si-bearing layers appears to be much more difficult, because of the breaking of the strong Si-O bonds. Indeed, the lattice friction opposed to [100](001) is on the order of 17 GPa<sup>19</sup>. On the basis of these results, strong CPO along (010) is thus expected in post-perovskite. However, this cannot be the end of the story because a crystalline aggregate must sustain some strain components along the three directions of space to satisfy strain compatibility. Therefore, to provide reliable models of crystal
<figure><img src="image_3.png" /><figcaption>**Figure 1.** ½[110](1 ¯10) edge dislocation core in MgSiO₃ post-perovskite computed at 120 GPa. (a) Burgers vector density computed from the differential displacement of atoms in the dislocation core region. (b) Atomic structure of the edge dislocation core with a perovskite-like stacking fault bounded by two partial dislocations of 1/6[110] and 1/3[110] Burgers vectors.</figcaption></figure>

preferred orientations and hence of seismic properties, it is necessary to understand which deformation mechanisms are active in this structure.

Mechanical twinning is a deformation mechanism that has received little attention despite microscopic observations of its occurrence in deformed CaIrO₃ post-perovskite¹³,¹⁴. In this paper, we show that [010] dislocations are not stable in MgSiO₃ post-perovskite, leading to partial dislocations that may be linked to mechanical twinning. Hence we present a hierarchical numerical model of the mechanical twinning in MgSiO₃ post-perovskite at 120 GPa³, which is compared with the dislocation activity to assess its possible relevance in plastic flow and CPO development in post-perovskite in the lowermost lower mantle.

## Results

MgSiO₃ post-perovskite phase exhibits an orthorhombic structure. The computed lattice parameters at a pressure of 120 GPa are *a* = 2.47 Å, *b* = 8.11 Å and *c* = 6.14 Å. The very strong anisotropy of the unit cell and the characteristic layering of the structure have led to a focus on plastic shear along [100]. Not surprisingly, this shear is very easy to produce along the magnesium layer in the (010) plane¹⁹,²⁰. The second easiest slip system also corresponds to shear in (010) but along the [001] direction²¹. Shear along [010] raises some questions, because it would require activating dislocations with a very large Burgers vector. This is usually unfavourable since the elastic energy of a dislocation scales to the square of the modulus of the Burgers vector. The classical Frank criteria show that [010] dislocations are not stable and decompose into two ½<110> dislocations because $b^2_{[010]} > 2xb^2_{½<110>}$. Indeed, in the orthorhombic C-lattice, ½<110>, with a length 4.2 Å, is one of the shortest lattice repeats of the structure and consequently represents another potential Burger vectors. Therefore, in this work, we start by addressing the properties of ½<110> dislocation cores.

### ½<110> dislocation cores in MgSiO₃ post-perovskite.
Atomistic calculations of the core structures of screw and edge dislocations with ½[110] Burgers vectors (performed at a pressure of 120 GPa to account for the D" layer conditions) indicate their strong tendency to spread, and hence to glide, in {110} planes. Indeed, regardless of the dislocation character, a full ½[110] Burgers vector spontaneously dissociates in {110} into two partial dislocations bounding a stacking fault with the typical perovskite-like octahedral interconnections by corners (Fig. 1b). The measurements of the partial dislocation Burgers vectors show that the two partial dislocations are asymmetric with $b_p = 1/6[110]$ and 1/3[110] (Fig. 1a). The equilibrium distance *R* between the partial dislocations (taken as the distance between the two maximum peaks of the Burgers vector density) is on the order of a few nanometres (Fig. 1a). Resulting from a balance between a repulsive elastic force and an attractive force associated with the fault formation energy, this large equilibrium distance suggests a very low stacking fault energy $\gamma_{isf}$ associated with the 1/6<110> {110} fault configuration. The partial dislocations are characterised by different Burgers vectors ($b_p$) and also respond to an applied stress with significantly different behaviours. Here, we describe the case of a dislocation with an initial edge character (the study of an initial screw dislocation leads to the same conclusions and is detailed in the supplementary materials). For a typical simulation cell containing approximately 40,000 atoms (350 Å along the direction of the dislocation core dissociation), we calculate that a shear stress in the range of 2.5–2.8 GPa triggers the expansion or the closure of the stacking fault, but this process is due to the sole displacement of the 1/6[110] partial dislocation (i.e. the second partial dislocation remains immobile). The onset of motion of the 1/6[110] partial dislocation is fairly insensitive to the initial stacking fault width or the investigated atomic system size. The 1/3[110] partial dislocation does not glide unless a high stress is applied. However, the 1/3[110] partial dislocation never glides as a partial dislocation, indeed we compute that the full dislocation first recombines into the compact ½[110] Burgers vector dislocation core (requiring an applied stress of 6.5 GPa for the 57 Å core depicted in Fig. 1) before it can actually further glide at an applied stress of 19 GPa.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Deformation twinning in (1 ¯10) with a 1/6[110] partial dislocation of the post-perovskite structure. Starting from the perfect lattice (a) viewed along [001], panels (b), (c) and (d) show the lattice with one-, two-, and three-layer faults after the shearing of successive 1/6[110] twinning dislocations. The four atomic structures correspond to the minimum energy configurations in the GPFE landscape of MgSiO₃ as depicted in Fig. 3.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** GPFE landscape for <110>{110} deformation twinning in MgSiO₃ and CaIrO₃ post-perovskites characterized by $b_p = 1/6<110>$. The calculated fault energies are given in Table 1.</figcaption></figure>

The conclusions of this study are that neither the [010] nor the ½[110] dislocations are stable or lead to active deformation mechanisms in MgSiO₃ post-perovskite. Instead, the only dislocation that is mobile is the 1/6[110] partial dislocation, which expands a stacking fault in its wake. Activation of partial dislocation is a deformation mechanism, known to occur in some materials (e.g., Si²² and SiC²³), but that may also lead to another mechanism: mechanical twinning. Indeed, because dislocation glide and deformation twinning are complementary mechanisms potentially contributing to plastic deformation, the occurrence of the highly mobile 1/6<110> partial dislocation raises the possibility of a mechanical twinning mechanism associated with the ½<110> {110} slip system in MgSiO₃ post-perovskite at 120 GPa. Figure 2 shows that if 1/6[110] partial dislocations are emitted and glide into successive planes, the characteristic structure of the stacking fault leads to the formation of a domain, which appears to be a twin of the orthorhombic structure. To the best of our knowledge, the occurrence of twinning has never been established in MgSiO₃ post-perovskite, although, interestingly, {110} twin domains have been detected in deformed CaIrO₃ post-perovskite samples quenched to ambient pressure¹³,¹⁴. These experimental observations in this analogue material, which is stable at ambient pressure, provide additional motivation for investigating ½<110> {110} deformation twinning in MgSiO₃ post-perovskite under pressure conditions corresponding to the D" layer. Because great care must be taken with plasticity mechanism interpretations based on analogues²⁴, we conducted calculations for both phases (CaIrO₃ and MgSiO₃ post-perovskites) to assess a possible bias (for instance, whether mechanical twinning might readily occur in CaIrO₃ post-perovskite but not in MgSiO₃ post-perovskite).

**Twinning mechanism and its associated energy landscape.** Deformation twinning is a deformation mechanism that is much more complex to describe than dislocation slip. Let us first consider the nucleation of a 1/6[110] partial dislocation (later called twinning dislocation) in the (1 ¯10) plane. The $b_p = 1/6[110]$ shear
<table><caption>**Table 1.** Characteristic parameters of ½<110> {110} deformation twinning in MgSiO₃ post-perovskite at 120 GPa and CaIrO₃ at 0 GPa. For the interpretation of the different parameters and energies, the reader is invited to refer to the Results section.</caption><thead><tr><th></th><th>MgSiO<sub>3</sub></th><th>CaIrO<sub>3</sub></th></tr></thead><tbody><tr><td>μ (GPa)</td><td>327</td><td>60</td></tr><tr><td>b<sub>p</sub> (Å)</td><td>1.41</td><td>1.73</td></tr><tr><td>ζ (Å)</td><td>3.2</td><td>3.9</td></tr><tr><td>s = b<sub>p</sub>/Δh</td><td>0.60</td><td>0.57</td></tr><tr><td>γ<sub>us</sub> (J/m<sup>2</sup>)</td><td>1.23</td><td>0.96</td></tr><tr><td>γ<sub>isf</sub> (J/m<sup>2</sup>)</td><td>0.69</td><td>0.62</td></tr><tr><td>γ<sub>ut</sub> (J/m<sup>2</sup>)</td><td>1.95</td><td>0.84</td></tr><tr><td>2γ<sub>isf</sub> (J/m<sup>2</sup>)</td><td>0.86</td><td>0.69</td></tr><tr><td>γ<sub>TM</sub> (J/m<sup>2</sup>)</td><td>1.09</td><td>0.15</td></tr></tbody></table>

displacement in (110) creates the intrinsic stacking fault (labelled *isf*), forming a one-octahedron-thick perovskite-like lamella (Fig. 2b) similar to that observed in ½<110> {110} dislocation cores (Fig. 1). If this event is followed by further successive nucleation and gliding of similar twinning dislocations in parallel successive {110} planes (Fig. 2c and d), a micro twin is formed. Once the {110} deformation twin lamella is nucleated, it can grow thicker by further activation of twinning dislocations on top of the twinning walls. The geometry of <110> {110} deformation twinning in post-perovskites can be described as the rotation of the parent lattice along the [001] axis by 34.5° in MgSiO₃ and by 31.1° in CaIrO₃.

In this work, we compute the twinning energy landscape (Fig. 3) corresponding to the twin formation mechanism described in Fig. 2 for the two post-perovskite compounds. This energy, also called generalized planar fault energy (GPFE)²⁵, corresponds to the cost per unit area required to form a N-layer twin by shearing N consecutive atomic layers along the [110] direction in the (110) plane. The GPFE first involves the γ<sub>us</sub> barrier against a one-layer partial fault becoming a one-layer full fault. This barrier is followed by the one-layer intrinsic stacking fault energy γ<sub>isf</sub> (Fig. 2b). Nucleation of the second, third and subsequent 1/6[110] dislocations creates the two-, three- and further N-layer stacking faults (Fig. 2c and d for instance). By analogy, the energy barrier opposed to the formation of each subsequent N-layer fault (N > 1), i.e., the barrier preventing a N-layer partial fault from becoming a N + 1-layer partial fault, is denoted γ<sub>ut</sub>, and the energy minimum, 2γ<sub>isf</sub>, is twice the energy of the twin stacking fault (i.e., it accounts for the upper and lower twin boundary). The energy difference between γ<sub>ut</sub> and 2γ<sub>isf</sub> defines the so-called twin migration energy γ<sub>TM</sub>²⁵,²⁶. As previously discussed by Kibey et al.²⁷ and Wang et al.²⁸, the pathway barriers γ<sub>us</sub>, γ<sub>ut</sub> and γ<sub>TM</sub> (Table 1) cannot be measured experimentally but nonetheless represent important parameters that strongly affect the critical twin nucleation stress σ<sub>cr</sub> as later described. For the investigated <110> {110} twinning system in MgSiO₃ and CaIrO₃ post-perovskites, the convergence in γ<sub>TM</sub> energy is reached after nucleation of the third twinning partial dislocation, thus resulting in total shear displacement by a full ½[110] lattice repeat. Hence, further nucleation and propagation of successive 1/6[110] dislocations enable twin growth on the developed three-layer twin lamella.

**Twinning stress mechanical model.** Considering that deformation twins usually form as individual thin plates embedded in the matrix or in contact with a free surface or a grain boundary²⁹, a deformation twin can be described as a series of loops or half-loops of twinning dislocations (Fig. 4a) belonging to twinning planes equidistant from each other by Δh. After it is fully nucleated, such a twin lens grows thicker (in the direction normal to the twinning plane) through successive nucleation of new dislocation loops. This thickening (growth) stage is generally easier than nucleation³⁰. The total energy associated with the process of twin nucleation³¹ can be defined as follows:

$$E_{tot} = E_{int} + E_{GPFE} + E_{l} - W \quad (1)$$

where E<sub>GPFE</sub> is the twin boundary energy²⁸; E<sub>int</sub> is the energy term resulting from the interaction of twinning dislocations; E<sub>l</sub> is the twinning dislocation line energy; and W is the work of the applied stress. In this equation, the energy term E<sub>l</sub> is independent of the spacing d between the twinning dislocation, and it consequently does not contribute to the critical twin nucleation stress σ<sub>cr</sub> and will be therefore disregarded in the following.

Being characterized by large twinning shear values (s = b<sub>p</sub>/Δh) close to 0.6, twin domains in MgSiO₃ and CaIrO₃ post-perovskites (Table 1) are very thin relative to their length. As shown by Cooper³¹, the total energy of an extended twin lamella of thickness h and half-length l, such that l ≫ h, can be considered as the energy of two flat surfaces containing twinning dislocations, i.e., the contribution of the {110} interplanar spacing Δh (2.41 Å and 3.06 Å in MgSiO₃ and CaIrO₃, respectively) versus the distance d between twinning dislocations can be neglected. In such a twin lamella, the twinning dislocations are assumed to belong to two flat surfaces and the total interaction energy of all twinning dislocation²⁸,³¹ is given by:

$$E_{int} = \frac{\mu b_p^2}{2\pi(1-\nu)} \left\{ N^2 \ln \frac{L}{d} - \ln(N-2)! + \sum_{i=2}^{N-1} [\ln(N-i)! + \ln(i-1)!] \right\} \quad (2)$$
<figure><img src="image_6.png" /><figcaption>**Figure 4.** (a) Schematic illustration of a semi-lenticular twin morphology characterized by thickness h, length l and invariant spacing d between the neighbouring twinning dislocations. The twinning plane is (110), and the twinning direction is [110]. (b) Evolution of the twin energy as a function of the spacing d between twinning dislocations. Under a typical applied stress of 100 MPa, the equilibrium distance of twin nuclei is approximately 40 Å.</figcaption></figure>

where N is the number of twinning dislocations of Burgers vector $b_p$; L is the arbitrary size of a crystal grain (Fig. 4a); v is the Poisson ratio; μ is the anisotropic shear modulus; and d is the distance between the i'th and the (i + 1)'th twinning dislocations (Fig. 4a), which is assumed to be invariant in this study.

The twin boundary energy $E_{GPFE}$ consists of two contributions: $\gamma_{SF}$, the energy required to create the intrinsic stacking fault, and $\gamma_{twin}$, the energy required to nucleate a twin. Relying on the computed characteristics of the GPFE landscape (Fig. 3) and on the fact that twinning dislocations may have classical properties of ordinary dislocations, these two terms can be defined as follows:

$$ 
\begin{align*}
\gamma_{SF} &= \gamma_{isf} + \left(\frac{\gamma_{us} - \gamma_{isf}}{2}\right) \left[1 - \cos\left(2\pi \frac{f(y)}{b_p}\right)\right] \\
\gamma_{twin} &= \left(\frac{2\gamma_{isf} + \gamma_{isf}}{2}\right) + \frac{1}{2}\left(\gamma_{ut} - \frac{2\gamma_{isf} + \gamma_{isf}}{2}\right) \left[1 - \cos\left(2\pi \frac{f(y)}{b_p}\right)\right] \\
\text{with } f(y) &= \frac{b_p}{2} + \frac{b_p}{N\pi} \left[\arctan\left(\frac{y}{\zeta}\right) + \arctan\left(\frac{y-d}{\zeta}\right) + \dots + \arctan\left(\frac{y-(N-1)d}{\zeta}\right)\right] \tag{3}
\end{align*}
$$

where $f(y)$ describes the disregistry function of N twinning dislocations characterized by the uniform half-width $\zeta$ and separated from each other by the invariant distance d (Fig. 4a). For MgSiO$_3$, the half-width $\zeta = 3.2$ Å can be directly extracted from the geometric characteristics of the 1/6[110] partial dislocation (Fig. 1), whereas for CaIrO$_3$, we use $\zeta = 3.9$ Å, accounting for the scaling factor $b_p(\text{CaIrO}_3)/b_p(\text{MgSiO}_3)$. For $\gamma_{SF}$, $f(y)$ is considered to be in the range $0 \le f(y) \le b_p$, corresponding to the intrinsic stacking fault region, while $\gamma_{twin}$ is computed for $b_p < f(y) \le Nb_p$ with $N > 1$, describing the disregistry function of the twin domain nucleated on the existing stacking fault configuration.

Relying on the defined analytic expressions for $\gamma_{SF}$ and $\gamma_{twin}$ (Eq. 3), the twin boundary energy can be computed as:

$$ E_{GPFE}(d) = \int_{0}^{d} \gamma_{SF} dy + (N - 1) \int_{0}^{d} \gamma_{twin} dy \tag{4} $$

Assuming the applied shear stress $\sigma_a$ to be uniform within the twin domain, the work of the applied stress can be expressed as follows:

$$ W = s\sigma_a h l = 2N^2\sigma_a b_p d \tag{5} $$

Relying on the energy terms described above (and reported in Table 1), the total energy of a twin lamella can be represented as the function of the spacing d between the twinning dislocations (Fig. 4b). The minimum energy configuration corresponds to the equilibrium distance $d_{eq}$, which increases slightly with applied stress. Thus, for the applied stress in the range of 20–200 MPa, $d_{eq}$ varies from 39 Å to 50 Å and from 12 Å to 18 Å in MgSiO$_3$ and CaIrO$_3$, respectively. According to these geometric parameters, the ratio $h/l$ of a twin lamella in the considered post-perovskites is on the order of $10^{-2}$.
<table><caption>**Table 2.** Calculated barrier energy for the three different plastic strain-producing mechanisms: individual partial dislocation slip, twinning and perfect dislocation slip.</caption><thead><tr><th></th><th>MgSiO<sub>3</sub></th><th>CaIrO<sub>3</sub></th></tr></thead><tbody><tr><td>Emission of 1/6[110] partial dislocation, γ<sub>ut</sub> (J/m<sup>2</sup>)</td><td>1.23</td><td>0.96</td></tr><tr><td>Emission of twinning dislocation, γ<sub>ut</sub> - γ<sub>isf</sub> (J/m<sup>2</sup>)</td><td>1.26</td><td>0.22</td></tr><tr><td>Emission of full ½&lt;110&gt; dislocation, γ<sub>ut</sub><sup>1/3[110]</sup> - γ<sub>isf</sub> (J/m<sup>2</sup>)</td><td>&gt;3</td><td>&gt;1</td></tr></tbody></table>

Finally, the critical twin nucleation stress $\sigma_{cr}$ can be further defined while minimizing $E_{tot}$ (Eq. 1) with respect to the distance $d$: $\partial E_{tot}/\partial d = 0$. Using $N=3$ defined from the GPFE calculations, we obtain $\sigma_{cr}$ values of ~880 MPa and ~620 MPa (see Supplementary Fig. 3) for the <110> {110} twinning system in MgSiO$_3$ and CaIrO$_3$ post-perovskites.

Discussion

At the microscopic scale, plastic deformation can result from several mechanisms. In this study we focus on intracrystalline mechanisms. Owing to the very peculiar crystal structure of post-perovskite, dislocation glide appears easy only in the plane of structural layering$^{19-21}$. However, which deformation mechanism effectively accounts for strain components out of the (010) plane remains an open question which is fundamental to address in order to model CPO fabrics and the resulting seismic properties. The fact that [010] and ½[110] dislocations spontaneously decompose into partial dislocations and that only one of these is able to move raises a fundamental question: which mechanism mediates plastic strain out of the (010) plane? Is it plastic slip of either perfect or partial dislocations, or is it mechanical twinning? These three mechanisms involve the emission of a leading 1/6[110] partial dislocation as the first step. In case of full dislocation slip, this first event (i.e., the nucleation of a leading partial dislocation) must be followed by the nucleation of the associated trailing 1/3[110] partial dislocation. In case of twinning, the nucleation of the first 1/6[110] partial dislocation is followed by the nucleation of a second 1/6[110] adjacent twinning dislocation. This process may or may not occur, because it may be more favourable to nucleate a second partial dislocation at a completely different location in the crystal rather than next to the first one, thus suggesting that extended partial slip may be a third possible strain-producing mechanism.

According to the work of Rice$^{32}$, the emission of a partial dislocation can be associated with the energy barrier $\gamma_{us}$ (Fig. 3 or Supplementary Fig. 2), which a partial dislocation must overcome in order to nucleate. Because twin formation requires at least the nucleation of a second twinning dislocation on a plane where a stacking fault already exists, similarly to Rice approach, Tadmor & Hai$^{33}$ have defined the energy barrier that a twinning partial dislocation must overcome as $(\gamma_{ut} - \gamma_{isf})$. For full dislocation slip, by analogy, an energy barrier can also be associated with the nucleation of the trailing partial dislocation by considering the difference between $\gamma_{isf}$ and the unstable energy corresponding to the trailing partial dislocation. Classically, in fcc metals, the emission of the trailing partial dislocation (which is identical to the leading one) involves thus the sole term $(\gamma_{us} - \gamma_{isf})$. The situation is different here for post-perovskite compounds, in which the trailing partial dislocation corresponds to a different Burgers vector. The emission of the trailing partial dislocation of partial Burgers vector 1/3<110>, which is basically not mobile, involves a second unstable energy $\gamma_{us}^{1/3[110]}$ that is higher than $\gamma_{us}$ (ca. greater than 4 J/m$^2$ in MgSiO$_3$, see Supplementary Fig. 2).

The competition between the three distinct mechanisms can thus be inferred by considering the different energy barriers associated with the onset of each mechanism (Table 2). Consequently, for the two post-perovskite phases, full dislocation slip is rather unlikely, because it requires overcoming the highest barrier ($\gamma_{us}^{1/3[110]} - \gamma_{isf}$) (Table 2), which is higher than the single $\gamma_{us}$. This conclusion is also consistent with our calculations of the lattice friction in MgSiO$_3$ post-perovskite presented in the first part of the results section. Indeed, calculation of the minimum stress for full dislocation glide yields a stress of 19 GPa. The glide of the full ½<110> dislocation is therefore very unlikely, and the occurrence of such a deformation mechanism in MgSiO$_3$ can be ruled out from both an energetic nucleation aspect and from consideration of the lattice friction.

The competition between the slip of individual partial dislocations and twinning can be assessed by comparing $(\gamma_{ut} - \gamma_{isf})$ and $\gamma_{us}$ (Table 2). For CaIrO$_3$ post-perovskite, the energy barrier that a twinning partial dislocation must overcome is clearly smaller than $\gamma_{us}$. Consequently, deformation twinning is expected to be highly favoured in this compound, in agreement with the observation of {110} twin domains in CaIrO$_3$ post-perovskite deformed experimentally$^{13,14}$. Additionally, long stacking fault ribbons (resulting from slip of partial dislocations) have never been reported in TEM observations of deformed post-perovskite analogues. The situation is not so straightforward in MgSiO$_3$ post-perovskite, because $(\gamma_{ut} - \gamma_{isf})$ is comparable to $\gamma_{us}$. In that case, the simple consideration of $\gamma_{ut}/\gamma_{isf}$ and $\gamma_{us}/\gamma_{isf}$ ratios is not sufficient for robust conclusions$^{34}$. The twinning mechanism efficiency is known to result from a complex interplay between the various energies associated with the GPFE energy landscape, as pointed out by Kibey et al.$^{27}$ in fcc Cu-Al alloys. Going beyond these simple energetic considerations is the goal of the hierarchical model described above, which computes the twinning stress on the basis of the entire GPFE landscape. For the MgSiO$_3$ post-perovskite phase, the critical twinning stress $\sigma_{cr}$ is approximately 880 MPa, computed according to the GPFE calculated at a confining pressure of 120 GPa. By itself, $\sigma_{cr}$ already appears to be a rather low value. This value can be further compared to the lattice friction opposed to the motion of a 1/6<110> partial dislocation (for the matter of comparison, we define the lattice friction as the computed minimum stress that must be applied to observe an expansion or a reduction in the stacking fault ribbon if the 1/6<110> partial core is forced to move under an applied shear stress). By computing the ½<110> dislocation core structure and the evolution of the core under applied shear stress, we find that the motion of the leading partial dislocation can occur only at a stress level of 2.7 GPa, which is 3 times higher than the critical twinning stress. In a deforming
grain of MgSiO₃ post-perovskite, the twinning stress will be reached well before the lattice friction of partial slip is overcome. This suggests a high ability for MgSiO₃ post-perovskite (as in CaIrO₃) to readily exhibit mechanical twinning, according to the present mechanism involving <110>{110} slip system.

Finally, the twinning stress, computed here at 0 K and supposed to be mostly athermal³⁰, can be compared to the lattice friction of dislocation glide in post-perovskite materials²⁰,²¹,²⁴. In CaIrO₃, the twinning stress reported here is at least 30% below any lattice friction level²⁴. For the MgSiO₃ post-perovskite, it turns out that the twinning stress is also below the stress level of lattice friction for the easiest slip systems since the Peierls stresses for [100] (010)²⁰ and [001](010)²¹ are 2 and 3 GPa, respectively. Thus, twinning appears to be the best candidate to account for the, so far, missing mechanism responsible for shear along the [010] direction. Up to now, such twinning mechanism has never been considered to explain CPO development in post-perovskite deformation, although the contribution of deformation twinning is especially important in low-symmetry crystals in which information on five independent slip systems is needed to describe a general deformation of the material.

## Methods

**First-principles calculations of GPFE.** Generalized planar fault energy (GPFE) calculations were performed based on the density-functional theory (DFT) within the generalized gradient approximation (GGA), with the PW91 parameterisation³⁵ and the all-electron projector augment-wave (PAW) method as implemented in VASP³⁶,³⁷. To achieve computational convergence, a plane-wave cut-off $E_{cut}$ of 600 eV was applied. The first Brillouin zone was sampled using the Monkhorst-Pack scheme³⁸ with an $8 \times 6 \times 1$ k-point grid. The convergence energy is $10^{-3}$ meV/atom. Simulations were performed at constant volume, corresponding to a bulk volume under confining pressure of 120 GPa and 0 GPa for MgSiO₃ and CaIrO₃, respectively. All calculations were performed using fully periodic atomic arrays containing 180 atoms. The simulation cells were built on vectors $a_1 = \frac{1}{2}[110]$, $a_2 = [001]$ and $a_3 = [1 \bar{1} 10]$ and oriented in such a way that the twinning direction $\frac{1}{2}[110]$ was aligned with the Cartesian $x$ axis and that the (1 ¯10) twinning plane was normal to $z$. Indeed, $a_3$ is nearly normal (89.1°) to the twinning plane. Along the $z$ axis, the simulation cells were as thick as 18 octahedral SiO₆/IrO₆ layers, i.e., 43.34 Å for MgSiO₃ and 54.99 Å for CaIrO₃. The interplanar distance $\Delta h$ between the subsequent {110} stacking planes was equal the size of one SiO₆/IrO₆ octahedron along $z$.

To compute the energy landscape, we create the first layer fault by incrementally shift half of the supercell by a twinning Burgers vector $b_p$. For the second layer fault, a new shear level is defined one octahedron above the already existing first layer fault. Atoms above this new level are rigidly shifted by the same incremental amount to reach a total displacement of $2b_p$. A similar procedure is repeated to create the third, fourth, etc layers. In order to keep the periodicity of the atomic array, the $a_3$ vector is tilted by the displacement vector. Practically, the path was sampled using eleven shear steps per twinning Burgers vector $b_p$. The full GPFE curve was initially computed using the classical molecular static simulations described below using pairwise potential of Buckingham form. After this first minimization, we selected four shear steps per twinning Burgers vector with increment $0.25b_p$. These configurations were relaxed using the VASP code. For the minimum energy configurations on GPFE curves (resulting from the shear by $N \times b_p$, where $N$ is a positive integer), full atomic relaxation, including the directions normal and parallel to the shear plane, was allowed. For the high-energy configurations, the degrees of freedom along the (110) plane were restricted only for the Si/Ir sublattice. To ensure optimum geometry of the faulted structure, additional shuffling across the boundary plane was allowed, i.e., all the atoms within one octahedral layer right above and below each newly introduced stacking plane (indicated with small black arrows in Supplementary Fig. 2b) were allowed to relax fully.

**Molecular static calculations of the $\frac{1}{2}<110>{110}$ dislocation core structure in MgSiO₃ post-perovskite.** For MgSiO₃ post-perovskite, atomistic simulations of dislocation cores at 120 GPa were performed using classical molecular statics (MS) simulations. The force field used corresponds to the Buckingham form of a pairwise potential parameterized by Oganov *et al.*³⁹. Calculations were performed using the LAMMPS package, relying on Ewald summation methods for Coulombic interactions⁴⁰. Optimization of the dislocation core configurations was performed using a conjugate gradient minimizer on a so-called periodic slab geometry containing a single dislocation core, as described by Hirel *et al.*⁴¹. The simulation cells were designed to be fully periodic along the dislocation line direction (Cartesian $x$ aligned with $\frac{1}{2}[110]$ for screw dislocations and with [001] for edge dislocations) and the dislocation glide direction (Cartesian $y$ aligned with [001] for screw dislocations and with $\frac{1}{2}[110]$ for edge dislocations). Along the direction normal to the {110} glide plane (Cartesian $z$), atoms at the bottom and at the top of the supercell were kept fixed to their regular positions given by the dislocation long-range elastic field. The total width (along $z$) of the layer with fixed atoms was equal to the potential short range cut-off distance of 12 Å. For screw dislocations, the supercell parameter aligned with $y$ was tilted by $\frac{1}{2}b$ along $x$ to maintain periodicity along the dislocation glide direction. The simulation cell size was equal to one unit cell repeat along $x$, 160 Å along $z$ and varies from 300 to 450 Å along $y$. Such simulation systems contained approximately 36,000–70,000 atoms. To initiate a dislocation glide, $\epsilon_{xz}$ and $\epsilon_{xy}$ strain components were gradually increased for screw and edge dislocations, respectively. To ensure quasi-static loading, the sheared atomic configuration was optimized every shear strain increment of $10^{-5}–10^{-4}$.

## References

1. Lay, T., Williams, Q., Garnero, E. J., Kellogg, L., & Wysession, M. E. Seismic wave anisotropy in the D" region and its implications. In *The Core–Mantle Boundary Region, Geodynamics Series* 28, 299–318, eds Gurnis, M., Wysession, M. E., Knittle, E., & Buffett, B. A., American Geophysical Union, Washington, D.C., USA (1998).

2. Grand, S. P. Mantle shear-wave tomography and the fate of subducted slabs. *Phil. Trans. Of the Royal Soc. A* **360**, 2475–2491 (2002).

3. Murakami, M., Hirose, K., Kawamura, K., Sata, N. & Ohishi, Y. Post-perovskite phase transition in MgSiO₃. *Science* **304**, 855–858 (2004).

4. Hernlund, J. W., Thomas, C. & Tackley, P. J. A doubling of the post-perovskite phase boundary and structure of the Earth's lowermost mantle. *Nature* **434**, 882–886 (2005).
5. Nowacki, A., Wookey, J. & Kendall, J. M. Deformation of the lowermost mantle from seismic anisotropy. *Nature* **467**, 1091–1096 (2010).

6. Merkel, S. et al. Deformation of (Mg,Fe)SiO₃ post-perovskite and D" anisotropy. *Science* **316**, 1729–1732 (2007).

7. Miyagi, L., Kanitpanyacharoen, W., Kaercher, P., Lee, K. K. M. & Wenk, H.-R. Slip systems in MgSiO₃ post-perovskite: Implications for D" anisotropy. *Science* **329**, 1639–1642 (2010).

8. Wu, X. et al. Seismic anisotropy of the D" layer induced by (001) deformation of post-perovskite. *Nat. Commun.* **8**, 14669 (2017).

9. Merkel, S. et al. Plastic deformation of MgGeO₃ post-perovskite at lower mantle pressures. *Science* **311**, 644–646 (2006).

10. Hirose, K., Nagaya, Y., Merkel, S. & Ohishi, Y. Deformation of MnGeO₃ post-perovskite at lower mantle pressure and temperature. *Geophysical Research Letters* **37**, 20302 (2010).

11. Miyagi, L., Kanitpanyacharoen, W., Kaercher, P., Stackhouse, S., Militzer, B. & Wenk, H.-R. The enigma of post-perovskite anisotropy: deformation versus transformation textures. *Phys. Chem. Min.* **38**, 665–678 (2011).

12. Walte, N. P. et al. Transformation textures in post-perovskite: Understanding mantle flow in the D" layer of the earth. *Geophys. Res. Lett.* **36**, 3–7 (2009).

13. Miyajima, N., Niwa, K., Heidelbach, F., Yagi, T. & Ohgushi, K. Deformation micro textures in CaIrO₃ post-perovskite under high stress conditions using a laser-heated diamond anvil cell. *Journal of Physics: Conference Series* **215**, 012097 (2010).

14. Niwa, K., Miyajima, N., Seto, Y., Ohgushi, K., Gotou, H. & Yagi, T. *In situ* observation of shear stress-induced perovskite to post-perovskite phase transition in CaIrO₃ and the development of its deformation texture in a diamond-anvil cell up to 30 GPa. *Phys. Earth Plan. Int.* **194–195**, 10–17 (2012).

15. Hunt, S. A., Walker, A. M. & Mariani, E. *In-situ* measurement of texture development in CaIrO₃ post-perovskite. *Phys. Earth Plan. Int.* **257**, 91–104 (2016).

16. Dobson, D. P. et al. Strong inheritance of texture between perovskite and post-perovskite in the D" layer. *Nat. Geosci.* **6**, 575–578 (2013).

17. Tsuchiya, T. & Tsuchiya, J. Structure and elasticity of Cmcm CaIrO₃ and their pressure dependences: Ab initio calculations. *Phys. Rev. B* **76**, 144119 (2007).

18. Hustoft, J., Shim, S., Kubo, A. & Nishiyama, N. Raman spectroscopy of CaIrO₃ post-perovskite up to 30 GPa. *Am. Mineral.* **93**, 1654 (2008).

19. Goryaeva, A. M., Carrez, P. & Cordier, P. Modeling defects and plasticity in MgSiO₃ post-perovskite: Part 2-Screw and edge [100] dislocations. *Phys. Chem. Min.* **42**, 793–803 (2015).

20. Goryaeva, A. M., Carrez, P. & Cordier, P. *Low viscosity and high attenuation in MgSiO3 post-perovskite inferred from atomic-scale calculations*. *Sci. Reports* **6**, 34771 (2016).

21. Goryaeva, A. M., Carrez, P. & Cordier, P. Modeling defects and plasticity in MgSiO₃ post-perovskite: Part 3-Screw and edge [001] dislocations. *Phys. Chem. Min.* **44**, 521–533 (2017).

22. Castaing, J., Veyssière, P., Kubin, L. P. & Rabier, J. Plastic deformation of silicon between 300 °C and 600 °C. *Phil. Mag. A* **44**, 1407–1413 (1981).

23. Pirouz, P. & Yang, J. W. Polytypic transformations in silicon carbide: the role of TEM. *Ultramicroscopy* **51**, 189–214 (1993).

24. Metsue, A., Carrez, P., Mainprice, D. & Cordier, P. Numerical modeling of dislocations and deformation mechanisms in CaIrO₃ and MgGeO₃ post-perovskites – comparison with MgSiO₃ post-perovskite. *Phys. Earth Planet. Int.* **174**, 165–173 (2009).

25. Ogata, S., Ju, L. & Yip, S. Energy landscape of deformation twinning in bcc and fcc metals. *Phys. Rev. B* **71**, 224102 (2005).

26. Kibey, S., Liu, J. B., Johnson, D. D. & Sehitoglu, H. Generalized planar fault energies and twinning in Cu-Al alloys. *Appl. Phys. Let.* **89**, 191911 (2006).

27. Kibey, S. A., Wang, L. L., Liu, J. B., Johnson, H. T., Sehitoglu, H. & Johnson, D. D. Quantitative prediction of twinning stress in fcc alloys: Application to Cu-Al. *Phys. Rev. B* **79**, 214202 (2009).

28. Wang, J. & Sehitoglu, H. Modeling of pseudotwinning in Fe₃Ga. *Mod. Sim. Mat. Sci. Eng.* **22**, 055008 (2014).

29. Cahn, R. Twinned crystals. *Adv. Phys.* **3**, 363–445 (1954).

30. Christian, J. W. & Mahajan, S. Deformation twinning. *Progress in Mat. Sci.* **39**, 1–157 (1995).

31. Cooper, R. The equilibrium shape of deformation twins. *Acta Metall.* **13**, 46–48 (1965).

32. Rice, J. R. Dislocation nucleation from a crack tip: An analysis based on the Peierls concept. *J. of the Mechanics and Physics of Solids* **40**, 239–271 (1992).

33. Tadmor, E. B. & Hai, S. A Peierls criterion for the onset of deformation twinning at a crack tip. *J. of the Mechanics and Physics of Solids* **51**, 765–793 (2003).

34. Van Swygenhoven, H., Derlet, P. M. & Froseth, A. G. Stacking fault energies and slip in nanocrystalline metals. *Nat. Mat.* **3**, 399–403 (2004).

35. Perdew, J. P. & Wang, Y. Accurate and simple analytic representation of the electron-gas correlation energy. *Phys. Rev. B* **42**, 13244 (1992).

36. Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. *Phys. Rev. B* **54**, 11169 (1996).

37. Blöchl, P. E. Projector augmented-wave method. *Phys. Rev. B* **50**, 17953 (1994).

38. Monkhorst, H. J. & Pack, J. D. Special points for Brillouin-zone integrations. *Phys. Rev. B* **13**, 5188–5192 (1976).

39. Oganov, A. R., Brodholt, J. P. & Price, G. D. Comparative study of quasiharmonic lattice dynamics, molecular dynamics and Debye model applied to MgSiO₃ perovskite. *Phys. Earth Planet. Inter.* **122**, 277–288 (2000).

40. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. *J. Comput. Phys.* **117**, 1–19 (1995).

41. Hirel, P., Kraych, A., Carrez, P. & Cordier, P. Atomic core structure and mobility of [100](010) and [010](100) dislocations in MgSiO₃ perovskite. *Acta Mater.* **79**, 117–125 (2014).

## Acknowledgements

This work was supported by funding from the European Research Council under the Seventh Framework Programme (FP7), ERC Grant No 290424—RheoMan. Computational resources were provided by the CRI-Université de Lille 1.

## Author Contributions

Ph.C. and P.C. designed the study. A.G. performed the atomistic calculations. All authors discussed, interpreted the results and contributed to the paper.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-18018-1.

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_7.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 Cre-
ative 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 per-
mitted 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) 2017