Hugues Raimbourg
email: hugues.raimbourg@univ-orleans.fr

Toshihiro Kogure

Tsuyoshi Toyoshima

Crystal bending, subgrain boundaries development and recrystallization in orthopyroxene during granulite-facies deformation

Keywords: Pyroxene, grain boundary, dislocations, recrystallization, TEM, EBSD

1-Abstract

A prominent feature of a granulite-facies shear zone from the Hidaka Main Zone (Japan) is the folding of orthopyroxene (opx) porphyroclasts. Dislocation density estimated by Transmission Electron Microscope (TEM) and chemical etching in homogeneously folded domains is too low to account for the amplitude of crystallographic bending, leading us to propose a model similar to "flexural slip" folding, where folded layers are μm-wide opx layers between thin planar clinopyroxene (cpx) exsolutions. Extension (compression) in the extrados (intrados) of the folded layer is accommodated by dislocations at the cpx-opx interfaces.

Alternatively to distributed deformation, crystal bending also localizes in grain boundaries (GB), mostly oriented close to the (001) plane and with various misorientation angles but misorientation axes consistently close to the b-axis. For misorientation up to a few degrees, GBs were imaged as tilt walls composed of regularly spaced (100)[001] dislocations. For misorientation angles of 7°, individual dislocations are no longer visible, but high-resolution TEM (HRTEM) observation showed the partial continuity of opx tetrahedral chains through the boundary. For 21° misorientation, the two adjacent crystals are completely separated by an incoherent boundary. In spite of these atomic-scale variations, all GBs share orientation and rotation axis, suggesting a continuous process of misorientation by symmetric incorporation of (100) [001] dislocations.

In addition to the dominant GBs perpendicular to the (100) plane, boundaries at low angle with (100) planes are also present, incorporating dislocations with a component of Burgers vector along the a-axis. The two kinds of boundaries combine to delimit subgrains, which progressively rotate with respect to host grains around the b-axis, eventually leading to recrystallization of large porphyroclasts.

2-Introduction

The interplay between deformation, grain size decrease and mechanical strength is a major process to account for the localization of the ductile deformation; deformation leads to dynamic recrystallization to a finer grain-size, which can cause a transition in deformation mechanisms and rock weakening, which, in turn, enhances deformation (e.g. [START_REF] Boullier | SP-Mylonites: Origin of some mylonites by superplastic flow[END_REF], [START_REF] Precigout | Strain localisation in the subcontinental mantle -a ductile alternative to the brittle mantle[END_REF], [START_REF] Raimbourg | Grain size reduction mechanisms and rheological consequences in high-temperature gabbro mylonites of Hidaka[END_REF]).

Though the microscopic mechanisms operative are not constant depending on strain rate and temperature [START_REF] Hirth | Dislocation creep regimes in quartz aggregates[END_REF], grain size reduction involves the creation of a large surface of grain boundaries and is therefore ultimately controlled by GB structure and properties.

One possible mechanism for recrystallization is the evolution of subgrain boundaries into "true" GBs with progressive misorientation. The transition between subgrain -regular arrays of dislocations within a single crystal-and GBs -incoherent zone separating two distinct crystals-is often associated with a misorientation angle of 10-15° (e.g. [START_REF] White | Geological significance of recovery and recrystallization processes in quartz[END_REF]), on the idea that for larger angles dislocations should be too closely spaced for coherent crystal to be present between them [START_REF] Read | Dislocations in crystals[END_REF]. On the other hand, high-resolution transmission electronic microscopy observations of the boundary core, for example in ZnO [START_REF] Sato | Arrangement of multiple structural units in a [0001] 49 tilt grain boundary in ZnO[END_REF] or in olivine [START_REF] Heinemann | Synthetic [100] tilt grain boundaries in forsterite: 9.9° to 21.5°[END_REF], revealed a structure made of dislocation-like units separating good crystal domains, with misorientations up to ~20°. This "low-angle" model of the boundary in terms of alternating good crystal and dislocation-like structural units was extended theoretically up to 53.1° of misorientation for [001] symmetrical boundaries in silicon [START_REF] Levi | Calculated structures for [001] symmetrical grain boundaries in silicon[END_REF]. Leaving apart these peculiar cases, where the boundary and the misorientation axis have particular crystallographic orientations, generally, high-angle GBs cannot be described like low-angle GBs only in terms of dislocations and other approaches need to be developed (see the review in [START_REF] Smith | On the general grain boundary[END_REF]). Independently of the atomic-scale descriptions considered, high-and low-angle GB have different properties in terms of energy [START_REF] Duyster | Grain boundary energies in olivine derived from natural microstructures[END_REF], element segregation [START_REF] Hiraga | Grain boundaries as reservoirs of incompatible elements in the Earth's mantle[END_REF]) and diffusivity. Another key property with respect to recrystallization is GB mobility: There is an increase in the mobility of the GB with misorientation, although there is no clear-cut transition between low-and high-angle GBs(see the review in [START_REF] Drury | Subgrain rotation recrystallization in minerals[END_REF]).

As most of the concepts and models on these topics are derived from works on metals, their applicability to silicates is not trivial, as special problems arise due to silicates' low symmetry and low number of slip systems with respect to metals. Deformation of opx is primarily achieved by the glide on (100)[001] slip system [START_REF] Etheridge | Deformation and recrystallisation of orthopyroxene from the Giles Complex, central Australia[END_REF][START_REF] Poirier | Plastic rheology of crystals[END_REF][START_REF] Steuten | An optical and electron microscopy study of defect microstructures in naturally deformed orthopyroxene[END_REF], far short of the five independent systems required for an arbitrary deformation geometry, but also insufficient to account for the development of subgrains, involved in models of recrystallization based on subgrain rotation (e.g. [START_REF] Guillopé | Dynamic recrystallization during creep of single-crystalline halite: an experimental study[END_REF]).

Deformation microstructure, such as crystal bending, subgrains and recrystallized grains, are extensively present in the mylonitic textures presented here and their boundaries have been exceptionally preserved from late-stage retrogression. Based on electron back-scattered diffraction (EBSD) and TEM observations, we propose in this study (i) to describe how the deformation of orthopyroxene prophyroclasts is accommodated within the crystal, including the development and progressive misorientation of GBs and (ii) to relate the properties of such intra-grain deformation structures with the whole process of recrystallization operative during mylonitization.

3-Geological settings

The mylonite studied here is located in the western border of the Hidaka Main Zone, a unit within Hidaka Metamorphic Belt interpreted as continental gabbros exhumed from the middle crust [START_REF] Arita | The metamorphic rocks and migmatites of the southern part of the Hidaka Metamorphic Belt[END_REF][START_REF] Jolivet | The Hidaka Shear Zone (Hokkaido, Japan): Genesis during a right-lateral strike slip movement[END_REF][START_REF] Miyashita | The basic plutonic and metamorphic rocks from the northern Hidaka Metamorphic Belts[END_REF][START_REF] Toyoshima | Tectonic evolution of lower crustal rocks in an exposed magmatic arc section in the Hidaka metamorphic belt, Hokkaido, northern Japan[END_REF]. The mylonite, constituted principally of elongated (up to a few mm) orthopyroxene (opx) grains embedded in a matrix of small (~50μm) plagioclase grains (details of the microstructure in [START_REF] Raimbourg | Grain size reduction mechanisms and rheological consequences in high-temperature gabbro mylonites of Hidaka[END_REF]), is characterized by a strong dextral deformation at granulite-facies conditions (T~800 -900℃ and P~6.1 -7.2 kbar) [START_REF] Osanai | Metamorphism and granite genesis In the Hidaka Metamorphic Belt, Hokkaido, Japan[END_REF][START_REF] Toyoshima | Granulite-grade gabbroic mylonites and crustal-scale decollement in a deep crust exposed in the Hidaka metamorphic belt[END_REF], which was completely preserved from lower-grade retrogression described in [START_REF] Kanagawa | Mylonitic deformation of gabbro in the lower crust: A case study from the Pankenushi gabbro in the Hidaka metamorphic belt of central Hokkaido, Japan[END_REF].

4-Internal structure of deformed opx grains

Grain-scale distribution of deformation

As described in detail in [START_REF] Raimbourg | Grain size reduction mechanisms and rheological consequences in high-temperature gabbro mylonites of Hidaka[END_REF] and [START_REF] Toyoshima | Tectonic evolution of the Hidaka metamorphic belt and its implication in late Cretaceous-Middle Tertiary Tectonics of Hokkaido, Japan[END_REF], the major feature of deformed opx grains is their elongation in the simple shear strain field, accompanied by the acquisition of a strong crystallographic preferred orientation (CPO) with the (001) planes concentrated around the foliation plane and the c-axes around the lineation.

Definitions

On the basis of their geometry at atomic-scale, which is defined later (See section 4.2), we define as low-, medium-and high-angle GB (LAGB, MAGB, HAGB), a boundary where the misorientation is in the range 0-5°, 5-15°, >15°, respectively. We differentiate recrystallized grains from subgrains as the former are completely independent grains while the latter preserve some connection, i.e. some elastic strain, with surrounding porphyroclast. In accordance with this definition, recrystallized grains (subgrains) are bordered by HAGB (LAGB/MAGB). Note that with this terminology a subgrain boundary is simply a particular type of GB.

Exsolution of clinopyroxene lamellae during deformation

Opx deformation is accompanied by the exsolution of a dense network of parallel lamellae, a feature absent in opx collected far from the mylonite. An energy dispersive X-ray spectrometer (EDS: JEOL JED-2200) in a TEM (JEOL JEM-2010UHR) unraveled a significant concentration of Ca within the lamellae. Crystallographic structure (Fig.

Supp

. Mat. 1) as well as this calcium content leads to interpret these lamellae as calcic clinopyroxene (cpx) exsolutions, similar to the larger, semi-coherent exsolutions we analyzed with electronic microprobe [START_REF] Raimbourg | Grain size reduction mechanisms and rheological consequences in high-temperature gabbro mylonites of Hidaka[END_REF]. These lamellae share the (100) plane with opx (Fig. Supp. Mat. 1) and constitute therefore an easy way to track orientation and potential distorsion of opx lattice.

Grain folding

In addition to their elongation, some opx grains are deformed with a "fold"-like geometry, i.e. with opposite hinges within a single grain, in contrast to simple crystal bending. Concomitantly with the deformation of the grain shape, opx crystallographic lattice itself is folded, as evidenced by the curvature of the densely distributed lamellae.

We carried out the characterization (axis and rotation angle) of such internal rotations using an electron-backscattered diffraction (EBSD) system (HKL, Channel 5) attached to an FE-SEM microscope (JEOL JSM-7000F).

Measurements at the points regularly spaced along the length of an opx with a curved shape (Fig. 1) show that the large internal rotations of the crystal structure are concentrated around the b-axes, themselves near the Y axis of the deformation ellipsoid, while the a-and c-axes are much more scattered within girdles. At most, between points 8 and 31, the internal rotation angle reaches as much as , while the rotation axis lies very close to the b-axis.

A fraction of the crystal rotation in Fig. 1 can be accounted for by several GBs, underlined by sharp orientation contrast in the fore-scatterred image (Fig. 1A), while elsewhere in the grain the rotation is more gradual, without visible contrast. To further unravel the grain internal microstructures, we carried out the chemical etching of polished surfaces, following the method described in [START_REF] Wegner | Chemical etching of amphiboles and pyroxenes[END_REF], which resulted in the preferential dissolution of (i) cpx lamellae parallel to the opx (100)

planes and (ii) GBs, which appear as thin black lines (Fig. 2). Please note that the etched surfaces, as well as thin sections described throughout this manuscript, are prepared in the XZ deformation plane, i.e. opx are cut in average close to their (010) plane.

Chemically etched sections confirmed the two concurrent modes of crystal bending, either homogeneous lattice distorsion, without visible boundary (Fig. 2C), or heterogeneous lattice distorsion (Fig. 2A andB), where crystal distorsion is concentrated in boundaries.

Homogeneous folding

In homogeneous folding, the progressive rotation of the lattice is homogeneously distributed through the folded domain. Discontinuities in crystal orientation, i.e. GBs, are absent, as evidenced by the fore-scattered image (Fig. 3A) and EBSD map (Fig. 3B).

This smooth bending is confirmed by the profile XY (Fig. 3D), where the crystal, with respect to reference point X, is progressively misorientated along the fold, by rotation about the b-axis (Fig. 3C).

Heterogeneous folding and formation of GBs

In contrast to homogeneous folding, many bent areas show sharp discontinuities in crystal orientation, revealed either by fore-scattered images and EBSD (Fig. 1 and4) or by chemical etching (Fig. 2A andB). The precision to be expected in the automatic mapping and detection of boundary by EBSD, which can be estimated in Fig. 3D from the scatter of the points around a smooth trend of increasing misorientation, is of the order of 1°. In contrast, etching made boundaries with misorientation as low as 0.2° degrees apparent (e.g. Fig. 2A) and constitutes therefore a finer tool to analyze grain internal structures.

Most of the GBs are straight and at large angle to the traces of opx (100) planes (Fig. 2A). The other kind of GB, less abundant, appear in etched surfaces as tortuous lines, at low angle to the traces of opx (100) planes (Fig. 2B). Finally, GBs with a structure intermediate between these 2 endmember geometries are rare.

Recrystallization

In some opx porphyroclasts, we observed also aggregates of subgrains/grains, localized on porphyroclasts rims or in cross-cutting deformation bands. The grains have a size from a few to a few tens of microns, misoriented around b-axis with respect to surrounding opx domains (Fig. 4A), with misorientation from a few degrees up to 30-40°. Their boundaries can be interpreted as an array of the two kinds of boundaries described above. In rare examples, the grain aggregates within an opx porphyroclast do not follow the same simple geometrical pattern (Fig. 4B) and seem to be randomly misoriented with respect to host grain. In such cases, grains are larger and devoid of cpx exsolutions.

In a final stage of deformation, these small grains get separated by on-going deformation from parent opx porphyroclasts to form asymetric tails [START_REF] Raimbourg | Grain size reduction mechanisms and rheological consequences in high-temperature gabbro mylonites of Hidaka[END_REF]. The crystallographic fabrics of these grains is random, as a result of deformation by a grain-size-sensitive deformation mechanism.

Microstructure and low-angle grain boundaries at atomic-scale

We carried out the atomic-scale observations of the deformation microstructures described above using a TEM (JEOL JEM-2010UHR) operated at 200kV on ultra-thin sections prepared from etched samples by focused ion beam (FIB: HITACHI FB-2100).

To image structures misoriented about the opx b-axis and as high-resolution TEM (HRTEM) imaging requires to align the incident electron beam along the crystallographic zone axis of low index, most sections were prepared perpendicular to the b-axis.

Lamellae

In all these ultra-thin sections, lamellae with an orientation parallel to (100) planes are present (Fig. Supp. Mat. 1). Their thickness is variable, with typical values of the order of tens to hundred Å, while their length is in most cases larger than 10μm. They are irregularly distributed, with an average spacing of the order of 1μm. Some of these lamellae concentrate dislocations on their surfaces, with a density often larger than surrounding opx domains.

Distributed dislocations

In spite of the large crystal distorsion unveiled by EBSD, in most ultrathin sections dislocations ditributed in the bulk of opx are surprisingly rare, which is in agreement with the relative lack of etch pits in opx compared to hosted inclusions (ilmenite or cpx).

They were found to be most abundant in homogeneously folded portions of the porphyroclasts (See 5.3), reaching concentrations of ~10 12 m -2 .

Straight LAGBs

Straight LAGBs are systematically oriented, in etched sections, close to the perpendicular to the traces of the opx (100) planes. Misorientation axis between adjacent domains is close to opx b-axis, though for very low misorientation angle (<0.5°) it was difficult to be determined precisely. As seen in TEM images, these boundaries, for low misorientation angle (<2-3°), are composed of regularly spaced dislocations, which appear strongly when h0l (small h index and l = 2,4,..) reflections are excited (Fig. 5A) and almost disappear when h00 reflections are excited (fig. 5B).

The contrast conditions, the geometry of the boundary as well as the misorientation axis are all fully consistent with the interpretation of these LAGB as tilt walls made of (100)[001] dislocations.

Dislocation spacing (d) decreases with the misorientation angle θ (Fig. 6), following a relationship close to the theoretical value for dislocation tilt walls [START_REF] Hirth | Theory of dislocations[END_REF]:

d ≈ B θ (Eq. 1)
where B, the length of the Burgers vector, was taken as unit cell parameter c ≈5.24 Å [START_REF] Van Duysen | Transmission electron microscope study of dislocations in orthopyroxene (Mg, Fe)2Si206[END_REF].

At higher resolution, these dislocations are slightly dissociated into two partials along the (001) plane (Fig. 7). In a thin domain, parallel to (100) plane, which includes and slightly extends over the partials, opx lattice periodicity is divided by two, as the result of ortho-to clino-enstatite transformation [START_REF] Buseck | High-resolution electron microscopy of silicates[END_REF].

Curved low angle grain boundaries

Another type of LAGB is constituted of tortuous boundaries at low angle to the traces of opx (100) planes. TEM imaging of these curved GB, with a misorientation of ~2°

principally around b-axis, showed that they are composed of dislocations, similarly to the straight LAGB (Fig. 8). The orientation of the GB is not constant; even at this very small scale, the boundary has a zigzag structure made of a succession of segments either perpendicular or at low angle to (100) planes. Portions at low angle to (100) planes have irregular features: the orientations of the dislocation lines as well as of the boundary itself vary (Fig. 8 (D) and (E)), strain patterns are often irregular, in contrast with the portions perpendicular to (100) planes (C). Our attempts (extinction of the dislocations or HRTEM imaging) to rigorously constrain the Burgers vector of the dislocations involved in boundaries at low angle to (100) planes were unsuccessful.

Evolution of GB structure with increasing misorientation

Although we could observe several kinds of GB in opx porphyroclasts, the straight ones perpendicular to opx (100) planes are by large the most abundant and have the simplest structure, so that they constitute an ideal target to study the evolution of the boundary with increasing misorientation.

Presence of dislocations and elastic strain

While in Fig. 5, where misorientation is ~0.5°, dislocations appear as black strain dots separating domains of the boundary where the lattice is less strained, in Fig. 8

(C),
where misorientation is ~1.5°, i.e. three times larger, strain fringes cover the whole boundary.

Higher misorientations, of the order of 5-15°, yield, for (100)[001] dislocations organized into tilt walls, theoretical spacing of the order of 3-1 a (Eq. 1). In such case, where calculated dislocation spacing is roughly equivalent to lattice cell parameters, the description of the boundary in terms of dislocations does probably not hold [START_REF] Read | Dislocations in crystals[END_REF]).

In contrast with LAGB (fig. 5), MAGB (Fig. 9) and HAGB (Fig. 10 and 11) do not appear an alignement of strain dots between continuous opx, but as a line separating two distinct opx domains. Some of the cpx exsolution lamellae end on the boundary, while some of them are continuous through it. When the h0l reflections with low h indices are excited, regularly spaced strain zones are apparent on GB with 7° (Fig. 9) and 11° (Fig. 10) of misorientation, while we could not observe any strain contrast near the GB with 21° of misorientation, whatever the diffraction conditions.

Atomic structures

In the 7°-misoriented LAGB (Fig. 9), the atomic scale structure of the boundary is well imaged especially with 400 lattice fringes (Fig. 9B), where every fourth fringe has a brighter contrast (probably because of a slight misorientation of the incident beam from the exact [010] direction). These lattice fringes correspond to the repetition of four units (a plane for single tetrahedral chains and a plane for octahedral slabs) within each lattice cell in the a-direction. The boundary varies between domains where the crystal on both sides is not connected, and domains where the fringes are continuous through the boundary. The nano-scale strain patterns associated (not to be confused with micron-scale ones seen with TEM in Fig. 9A), have a pseudo-periodic spacing roughly equivalent to dislocation spacing in the dislocation tilt wall model. At the sub-lattice scale, some 400 lattice fringes are continuous and some are not, which must correspond to the fact that some tetrahedral chains and octahedral slabs are continuous and some are broken, a structure that cannot be expressed by conventional dislocation scheme.

For 11° of misorientation, i.e. at the MAGB/ HAGB transition, although micron-scale strain is visible when h0l reflections are excited and disappear for h00 excitation(Fig. 10), HRTEM images are equivocal as to the continuity of (100) planes through the boundary. Finally, for HAGB (misorientation ~21°), no micron-scale strain is visible and the boundary appears in HRTEM as a band of width 1nm or less made of bright, incoherent material (Fig. 11).

Three-dimensional geometry of the boundary

The macroscopic, geometrical features of the GB oriented close to the perpendicular to the traces of opx (100) planes seem to be independent of the misorientation: They are relatively straight and the misorientation axis is close to opx b-axis (e.g. Fig. 11C).

As regards the boundary geometry, all our observations, restricted to a plane of section, In conclusion, the geometry of the GB is constant and independent of its atomic-scale structure, which varies from coherent to incoherent with increasing misorientation.

5-Discussion

Active slip systems in opx

Although the clino-to orthoenstatite transformation has been proposed as a deformation process in opx [START_REF] Coe | The orthoenstatite to clinoenstatite transformation by shearing and reversion by annealing: Mechanism and potential applications[END_REF][START_REF] Coe | Crystallographic orientation of clinoenstatite produced by deformation of orthoenstatite[END_REF][START_REF] Mclaren | A transmission electron microscope study of naturally deformed orthopyroxene. I. Slip mechanisms[END_REF][START_REF] Raleigh | Slip and the clinoenstatite transformation as competing rate processes in enstatite[END_REF], the paucity of clino-enstatite lamellae observed in our samples (mostly related to the dissociation of (100)[001] dislocations into two partials, see Fig. 7) as well as the very large elongation (shape ratio correspond to finite shear strain locally up to more than 20), imply that most of the deformation results from dislocation activity.

The dislocation system involved in plastic deformation is principally (100)[001], as evidenced by (i) the porphyroclast CPO, (ii) the GB geometry and misorientation axis and (iii) the direct characterization by TEM of the nature of the dislocations involved in GB. This slip system is the most common one in opx, both from natural and experimental samples [START_REF] Kohlstedt | Transmission electron microscopy investigation of the defect microstructure of four natural orthopyroxenes[END_REF][START_REF] Naze | A TEM study of lattice defects in naturally and experimentally deformed orthopyroxenes[END_REF][START_REF] Van Duysen | Transmission electron microscope study of dislocations in orthopyroxene (Mg, Fe)2Si206[END_REF]. Several other minor systems are described by these authors, which have as Burgers vector either b or c. [START_REF] Steuten | An optical and electron microscopy study of defect microstructures in naturally deformed orthopyroxene[END_REF] noted that none of these slip systems can account for the development of GBs parallel to (100) planes. We face in our study the same question, as the dislocations present in the boundary parallel to (100) planes and which accommodate a rotation around b (Fig. 8) have by construction, as Burgers vector, the a-axis [START_REF] Hirth | Theory of dislocations[END_REF]). The spacing between the elastic strain concentrations observed in these boundaries (Fig. 6 and8E) is also in agreement with the a-axis as Burgers vector.

On the other hand, the detailed structure of the boundary parallel to (100) planes appears much more complex than this simple model (Fig. 8C to E), as the boundary is not straight and perfectly symetric with respect to adjacent opx domains and the strain concentrations are irregularly distributed along the boundary. The lattice parameter a is much larger (~18.2 Å (Van Duysen, Doukhan and Doukhan 1985)) than b or c, so that the energy of dislocations with the a-axis as Burgers vector is probably quite high and such dislocations are unfavoured. Furthermore, as four silicate chains are superposed in the a direction within a unit cell, one may suppose that these chains terminate at different points within the boundary and that the associated dislocation is dissociated.

In spite of the complexity of the boundary parallel or at low angle to a-planes, it is unambiguously made of dislocations, which implies the activation, in addition to the dominant (100)[001], of another slip system, whose Burgers vector, if not simply a, has some component along the a-axis. Such dislocations, which are not observed elsewhere in the crystal nor reported elsewhere, may be activated in very specific cases where they are geometrically necessary, for example in segments of GB joining tilt walls made of (100)[001] dislocations. In materials such as silicates, the low number of easy slip system may lead, because of locally high stresses, to the activation of less-favoured systems.

Homogeneous lattice folding

Geometrically necessary dislocations

From the theory of bending of plates [START_REF] Timoshenko | Theory of elasticity[END_REF], if the opx lattice curvature was purely elastic, the magnitude of stresses acting to bend the crystal would be of the order of σ~E Considering the example of the opx homogeneously bent in Fig. 13, lattice curvature radius is of the order of R=700μm, yielding, for a Burgers vector of length c ≈ 5.24Å, (Van Duysen, Doukhan and Doukhan 1985)), a theoretical dislocation density of 2.7/μm 2 . Actual dislocation density, calculated on the TEM micrograph by counting dislocations away from cpx lamellae (Fig. 13C) or on the SEM image by counting triangular etch pits away from lamellae (Fig. 13B), is in the range 1-1.2/μm 2 , i.e.

roughly half of the theoretical density.

One can wonder whether our low dislocation density is not the result of repeated unfortunate sampling, where we would have missed dislocation tilt walls concentrating the curvature outside of the zones we analyzed by TEM or SEM. Our "limit of detection" of tilt walls corresponds to misorientation of ~0.2°, above which they are unraveled by etch pit as thin black lines. Assuming that all the curvature is accommodated by 0.2° misorientation tilt walls, the radius of curvature R=700μm

implies that these walls should be spaced by 2.4μm, i.e. there should be several of them visible in the TEM ultrathin section, whose width is ~10μm. As this is not the case, lattice curvature has to be accommodated homogeneously, i.e. by dislocations distributed more or less evenly in the volume.

Model of "flexural slip" folding

While the actual density of dislocations in opx bulk is too low to explain the lattice curvature, HRTEM images unveiled the presence, at the interface between folded cpx lamellae and opx layers, of densely distributed dislocations (Fig. 14).

To understand the role of these dislocations in folding, let us assume that opx is not folded as a thick monolayer but as a stack of μm-thick layers, separated by cpx exsolution lamellae (Fig. 15A). Let us further suppose that these lamellae can accommodate differential displacement, similarly to macroscopic models of "flexural slip folding" [START_REF] Ramsay | Folding and fracturing of rocks[END_REF]). If the folding was purely elastic, then for a layer of thickness t, curved with a radius R, a given length l (in unstressed state) parallel to c-direction along the extrados (intrados) of the layer is stretched (shortened) to a length

l(1 + t 2R ) (l(1 - t 2R
)). In reality, this deformation is not elastic and we suppose that it is accommodated by regularly spaced dislocations on opx layers interfaces (i.e. along cpx lamellae), which account for the difference in (001) plane density between two adjacent layers (Fig. 15C).

The difference in displacement between upper and lower surfaces is:

l (1 + t 2R ) -l (1 - t 2R ) = l t R
If this difference is taken equal to the distance between (001) planes (i.e. d 001 ), mismatch between upper surface of the lower layer and lower surface of the upper layer is equivalent to inserting an extra c-plane on one side of the boundary, i.e. an edge dislocation on the boundary (Fig. 15C). On layer interface, the spacing l between dislocations accommodating difference in c-plane density is thus:

l t R = c → l = Rc t (Eq.
3)

The dislocation density associated with this distribution is finally given by:

ρ = 1 lt = 1 Rc
i.e. the density is the same as Eq.2, which was calculated for homogeneously distributed dislocations.

"Flexural slip" folding and concentration of dislocations on opx layer interfaces (i.e.

parallel to (100) planes) constitutes therefore another possible distribution of dislocations to accommodate lattice curvature, alternative to either (i) dislocations homogeneously distributed in the bulk or (ii) dislocations organized into tilt walls (i.e.

parallel to (001) planes). On the example of Fig. 13, it accounts for about half of the curvature. The occurrence of flexural "flexural slip" folding requires the presence of decoupling interfaces. This role is played by exsolutions, but not in the sense that there is actually slip along the lamellae, but rather in the sense that exsolution process itself, through the growth of the lamellae, enables the opx on both sides of the lamellae to move with respect to each other, similarly to the way synkinematic veins accommodate differential displacement of their rims [START_REF] Passchier | Microtectonics[END_REF].

Nature and role of the dislocations at cpx-opx interfaces

Van Duysen, Doukhan and Doukhan (1985) interpreted these dislocations (Fig. 14) as thermal mismatch dislocations, emplaced as a result of the difference in d 001 between cpx and opx at temperature of the deformation (~800°C) and subsequently preserved throughout cooling. The model of flexural slip we propose above (Fig. 15A) leads to an alternate explanation, where these dislocations accommodate the lattice curvature (Fig.

Supp. Mat. 3).

Dislocation spacing cannot be directly estimated using [START_REF] Van Duysen | Transmission electron microscope study of dislocations in orthopyroxene (Mg, Fe)2Si206[END_REF]'s results because of a mistake in opx thermal expansion coefficient.

After correction and integration of the effect of pressure (See Appendix), mismatch dislocation spacing is of the order of a few hundreds of nanometers, with further precision impossible because of uncertainty in crystallographic parameters.

In the "flexural slip" folding model, on the other hand, with R~1mm, t~1-5μm (Fig. 13), the dislocation spacing predicted by Eq.3 is in the range 100-500nm. Actual average spacing of dislocations at cpx-opx interfaces (Fig. 14) is in good agreement with values from both models.

One can nevertheless note that the dislocation density is highly variable along a given lamellae and even varies between its upper and lower surface (Fig. 14). Such features

Energy budget in flexural slip folding

The total energy stored in the crystal when dislocations are arranged according to "flexural slip" model is difficult to compute. On the other hand, one can qualitatively estimate the variations in energy when, within the model framework, the density of decoupling interfaces diminishes/layer thickness increases.

To estimate energy budget, let us take as example the case where the thickness of opx layer is 5μm, the curvature radius 1mm, Eq. 3 predicts that dislocation spacing is ≈ 100nm along the cpx lamellae. In contrast, the distance between dislocations in two adjacent lamellae is of the order of opx layer thickness, 5μm, i.e. 50 times farther than neighboring dislocations on the same lamellae. As elastic stresses associated with dislocations scale as r -1 (Hirth and Lothe 1968), where r is the distance to the dislocation core, one can neglect, in the energy of the dislocation distribution, the interaction energy between adjacent layers. The interaction energy is thus principally the result of parallel (100)[001] dislocations located on a common (100) plane. Let us consider two neighboring dislocations on (100) plane. Their energy of interaction W per unit length L is given by:

W L = - μ 2π(1 -ν) 𝑐 2 ln r r α
Where μ, ν, 𝑐 , r α are respectively the shear modulus, Poisson ration, unit cell parameter and a constant for the logarithmic argument to be dimensionless (Hirth and Lothe, Eq-5.17 ( 1968)). This energy increases as r, in other words dislocation spacing l in Eq.3, decreases. The total energy of a unit volume of folded pyroxene is, when considering only the interaction energy of neighboring dislocations on (100) plane, porportional to Wρ, where ρ is the dislocation density of Eq.2.

Let us finally compare two configurations of flexural folding, the first one where there are few, thick opx lamellae and the other one where there are many, thin opx lamellae.

Dislocation spacing l is inversely proportional to the thickness of opx lamellae, hence W increases as thickness increases. On the other hand, dislocation density is constant, irrespective of the opx layer thickness. As a result, the volumic energy of folded opx decreases as the thickness of the folded layers decreases. The case is much similar to the elastic case, where the stresses necessary to bend a layer scale as σ~E t R , i.e. σ decreases as the layer thickness decreases. In terms of energy, it is easier to bend plastically or elastically a multilayer assemblage made of numerous, thin layers than few, thick layers.

As thickness is inversely proportional to the density of decoupling surfaces, hence on density of cpx exsolutions, "flexural slip" folding energy is controlled by the process of exsolution.

Subgrain/grain transition

While there is some regularly spaced elastic strain on the GB with 11° misorientation shown in Fig. 10, other GB we analyzed had a similar misorientation but no elastic strain. The transition between MAGB and HAGB, i.e. between subgrain and recrystallized grains, occurs therefore probably for misorientation in the range 10-15°.

10° corresponds to a critical dislocation spacing estimated as d crit ≈ 𝑐 θ crit ≈ 30nm ≈ 1.6 |𝒂|, i.e. for dislocations spacing smaller than 1.5 times the lattice parameter in the direction of the boundary where dislocations pile up, lattice distortions in the boundary become too large for the two adjacent crystals to remain connected.

In olivine, Heinemann, Wirth, Gottschalk and Dresen (2005) observed GB made of regularly spaced dislocations for tilt angles up to 21.5°, which correspond to a spacing

d crit = 𝑐 θ crit ≈ 16nm ≈ 1.6 |𝒃|,
where 𝒃 is the cell dimension along the dislocation array (equivalent to 𝒂 in the opx tilt walls). Therefore, the large difference of transition between low-and high-angle grain boundary structure between olivine and pyroxene might be explained in terms of critical spacing of dislocations in the boundary.

In the transitional domain of MAGB, where lattices are connected through the boundary but dislocations are too densely distributed to be individualized, we observed periodic strain along the boundary, with a spacing ~150 nm (Fig. 9), i.e. much larger than spacing calculated from theoretical law. Such long-range periodic strain pattern was also observed in olivine GB with large misorientation by Heinemann, Wirth, Gottschalk and Dresen ( 2005), who related its periodicity to the coincidence, on the GB, of lattice sites from both crystals.

Models of recrystallization

Any model invoked to account for recrystallization involves at some point the increase in GB surface with respect to original, undeformed grain. The two mechanisms of crystal folding described here have contrasted effects on recrystallization:

Homogeneous folding (4.1.4) does not lead to GB growth, on the contrary to the localized accommodation of bending, resulting in organization of dislocations into tilt walls and further evolution into GBs (4.1.5). Consequently, bending is a factor of recrystallization only if deformation is heterogeneous at the grain scale. Our observations do not provide a clear pattern as to which factors controls whether bending is homogeneous or heterogeneous. Nevertheless, on the basis of the energy analysis in 5.4, homogeneous folding is favoured when the thickness of folded opx layers is small, i.e. when exsolutions are densely distributed. The process of cpx exsolution is therefore potentially controlling the style of folding in elongated opx porphyroclasts, hence its potential for recrystallization.

Considering the cases where GB develop in response to folding, two distinct mechanisms may lead to recrystallization: Migration of grain boundaries from low-to high-density of dislocations domains [START_REF] Sellars | Recrystallization of metals during hot deformation[END_REF][START_REF] Tullis | Dynamic recrystallization of feldspar: a mechanism for ductile shear zone formation[END_REF] or dislocation climb and organization into tilt walls and formation of subgrains [START_REF] Guillopé | Dynamic recrystallization during creep of single-crystalline halite: an experimental study[END_REF]. [START_REF] Hirth | Dislocation creep regimes in quartz aggregates[END_REF] describe for example in quartz several regime of dislocation creep, based on different combination of these two possible mecanisms for restoration. For opx, data is limited to [START_REF] Etheridge | Deformation and recrystallisation of orthopyroxene from the Giles Complex, central Australia[END_REF]'s analysis of recrystallization as the result of the migration of high-angle grain boundaries, generated in "kink bands" domains, which are similar to the GB at large angle to (100) planes traces that are described here.

Subgrain rotation-recrystallization

As already proposed, on the basis of CPO, in [START_REF] Raimbourg | Grain size reduction mechanisms and rheological consequences in high-temperature gabbro mylonites of Hidaka[END_REF], all the characteristics of the GBs we described here (geometries and atomic scale structure of the LAGB, misorientation axes) point to the major role played by (100)[001] slip system and associated rotation in subgrain formation. The subgraingrain evolution is a bit more problematic: As described in halite by [START_REF] Drury | Subgrain rotation recrystallization in minerals[END_REF] and [START_REF] Pennock | The development of subgrain misorientations with strain in dry synthetic NaCl measured using EBSD[END_REF], the model of recrystallization where all subgrains rotate progressively with strain, eventually reaching a misorientation large enough for them to become independent grains may be oversimplified, as these authors described several kind of subgrains, of which only certain ones evolve into recrystallized grains with increasing strain. In this analysis, among the low-angle GB formed during deformation, [START_REF] Drury | Subgrain rotation recrystallization in minerals[END_REF] makes the distinction between incidental and geometrically necessary boundaries (Kuhlmann-Wilsdorf and Hansen 1991), only the latter evolving potentially into high-angle GB separating recrystallized grains.

As shown by their CPO [START_REF] Raimbourg | Grain size reduction mechanisms and rheological consequences in high-temperature gabbro mylonites of Hidaka[END_REF], elongated opx porphyroclasts are in average oriented so that the c axes are close to elongation axis X and (100) planes close to the foliation plane XY. As a result, slip on the dominant (100)[001] slip system allows only further elongation parallel to X. The GBs parallel to (001) planes, also made of (100)[001] dislocation, combine to fold the whole grain, hence they enable shortening parallel to c-axis, a mode of deformation which is not possible simply by slip. In that respect, these GBs can be considered as geometrically necessary, in the sense that they are controlled at large scale, by the deformation field of the grain and its surroundings. Similarly, as they accommodate the crystal rotations around the b-axis induced by the GBs parallel to (001) planes, GBs parallel to (100)

planes can be also considered as geometrically necessary, which can explain the activation of slip systems whose Burgers vector has a component along a-axis, not reported so far in deformed opx and probably not an "easy" slip system. Furthermore, EBSD patterns of recrystallized grains (Fig. 5 of [START_REF] Raimbourg | Grain size reduction mechanisms and rheological consequences in high-temperature gabbro mylonites of Hidaka[END_REF]), show that all grains are derived from parent grain by rotation only in a single direction (counterclockwise in the present case), which means that these rotations are not incidental but are kinematically necessary, i.e. they accomodate much larger-scale (e.g. porphyroclast scale) deformation. The corollary of this property is that rotations are expected to increase with strain, resulting in the evolution of subgrains into recrystallized grains. A large fraction of recrystallization in opx seems therefore to be related to localized bending within large porphyroclasts.

Grain boundary migration

Extensive migration of HAGB is not supported by their geometry: Although they are incoherent, most of them preserve the geometry of low angle GB, i.e. boundary plane is bissector of (100) planes of opx domains and rotation is about opx b-axis (e.g. Fig. 11, 12 and 4 (A)). This implies that they grow from the incorporation of (100) [001] dislocations and that their mobility, if they have some, is restricted to the c direction.

The major obstacle to GB movement, is the presence of coherent cpx exsolutions along the boundary, which pin it and hinder its mobility (e.g. Herwegh and Kunze ( 2002 is the result of deformation-assisted diffusion, which may greatly enhance Ca mobility.

As a result, grain boundary migration is possible, but only in a late stage of porphyroclast disruption. Recrystallization in opx is therefore a multi-staged process, dominated by subgrain rotation recrystallization but where GB migration is eventually activated after expulsion of cpx lamellae.

6-Conclusion

The plastic deformation of opx porphyroclasts of the Hidaka mylonite zone led, in addition to strong grain elongation and CPO acquisition, to extensive intra-grain deformation, which manifests as lattice folding about b-axis as a result of the presence of (100)[001] dislocations. This deformation is accommodated either homogeneously through the grain or by the formation of boundaries concentrating the misorientation. In the former case, dislocations distributed in the bulk opx are too few for the lattice curvature. To account for the discepancy, we developed a model similar to flexural slip folding, where opx layers, folded between cpx exsolutions, are able to slip with respect to each other, which results in the dense distribution of dislocations along the slip interfaces, i.e. along cpx lamellae surfaces. In the alternative mode of deformation, where curvature is accommodated by the formation of boundaries, we observed a transition between subgrain (connected crystals through the boundary) and recrystallized grain (disconnected crystals) around 10-15°. In addition, while for LAGB, individual dislocations are apparent, for MAGB, though elastic strain patterns are visible, the boundary appears as a continuous, thin frontier and individual dislocations are no longer visible. The connection of the crystals on the two sides of the boundary is only apparent at atomic-scale, as some lattice fringes are continuous through the boundary. Finally, as recrystallization is the result of the formation of subgrains and their progressive misorientation, the two modes of grain folding (homogeneous/heterogeneous) have a great impact on the potential to recrystallize, as only heterogeneous recrystallization favours the formation of GBs.

7-Appendix

We estimated the length of opx and cpx unit cell c parameter for T=800°C and P=0.7

GPa as:

𝑐 = 𝑐 0 (1 + 800α T + 0.7α P )
where 𝑐 0 is the value for ambient P and T and α T and α P expresses the thermal dilatation and compressibility coefficients. For opx, α T ≈ 0.16 10 -4 /℃, irrespective of its composition [START_REF] Frisillo | Linear expansion coefficients of orthopyroxene to 1000°C[END_REF][START_REF] Sueno | Orthoferrosilite: high temperature crystal chemistry[END_REF]) and α P ≈ -3.5 10 -3 / GPa [START_REF] Angel | Equations of state and thermodynamic properties of enstatite pyroxenes[END_REF] for enstatite). For cpx, α T ≈ 0.06 10 -4 /℃ [START_REF] Cameron | High temperature crystal chemistry of acmite, diopside, hedenbergite, jadeite, spodumene and ureyite[END_REF] for either hedenbergite of diopside) and α P ≈ -2.8 10 -3 /GPa [START_REF] Levien | High-pressure structural study of diopside[END_REF] for diopside).

In other words, pressure-dependent variations are similar for cpx and opx, while thermal expansion of opx is much larger than in cpx. In contrast, unit cell parameter 𝑐 0 is consistently larger in cpx [START_REF] Levien | High-pressure structural study of diopside[END_REF][START_REF] Salviulo | C2/c pyroxene from two alkaline sodic suites (Western Ross Embayment-Antarctica): crystal chemical characterization and its petrologic significance[END_REF] than in opx [START_REF] Angel | Equations of state and thermodynamic properties of enstatite pyroxenes[END_REF][START_REF] Smyth | An orthopyroxene structure up to 850°C[END_REF][START_REF] Sueno | Orthoferrosilite: high temperature crystal chemistry[END_REF]. To take into account the large compositional dependence of 𝑐 0 , we considered compositions as close as possible to ours [START_REF] Raimbourg | Grain size reduction mechanisms and rheological consequences in high-temperature gabbro mylonites of Hidaka[END_REF]: For opx, [START_REF] Smyth | An orthopyroxene structure up to 850°C[END_REF] reported 𝑐 0 opx = 5.232Å, while for cpx, [START_REF] Salviulo | C2/c pyroxene from two alkaline sodic suites (Western Ross Embayment-Antarctica): crystal chemical characterization and its petrologic significance[END_REF] reported 𝑐 0 cpx1 = 5.26Å and 𝑐 0 cpx2 = 5.28Å for two samples of similar composition (VG3187 and VG3185Pc).

When considering the first value 𝑐 0 cpx1 , parameter c has roughly the same length in opx and in cpx for T=800°C and P=0.7 GPa, i.e. no mismatch dislocations are required. For the same P-T conditions, when considering the second value 𝑐 0 cpx2 , 𝑐 cpx2 is larger than 𝑐 opx by ~0.2%, which implies dislocations spaced by 250nm.

These calculations tends to show that 𝑐 cpx > 𝑐 opx at high P and T, requiring emplacement of mismatch dislocations along the interface. Our own HRTEM observations of the polarity of regularly spaced dislocations along the interface of a lamellae in a domain devoid of bending show similarly that 𝑐 cpx > 𝑐 opx . On the other hand, the large uncertainties in unit cell parameter at ambient conditions, of the order of the pressure and temperature dependent variations, preclude to precisely assess the density of such mismatch dislocations and we retain therefore only the possible magnitude of dislocation spacing, as of the order of a few hundreds of nanometers. where d 100 is half that of ortho-enstatite [START_REF] Buseck | High-resolution electron microscopy of silicates[END_REF]. One can also note that the clino-enstatite domain (enclosed in dashed white line in Fig. 7B) is not restricted to the interval between the two partials but extends further outside of it, as apparent in the half-spaced fringes with three dark lines, interpreted as clino-enstatite with a unit cell width. EPMA mapping of Ca distribution within a deformed opx porphyroclast. Recrystallized opx, either within the porphyroclast ("intra-opx deformed domain") or on its rim ("small opx domain") have a lower Ca content than the porphyroclast ("coherent opx domain"). Furthermore, within the recrystallized grains on the upper right of (A), grain boundaries are visible in (C) because of their low Ca concentration compared with the grains interiors. Within recrystallized zones, abundant cpx grains probably result from the depletion in Ca of the opx. Note that when considering a single boundary cutting across the opx, the Ca-depletion is visible only when the misorientation is very large ("GB1"), while no chemical heterogeneity is apparent for a lower misorientation ("GB2").

  are insufficient to fully constrain it. Nevertheless, knowing the crystallographic orientation of both sides of the boundary, we can determine what orientation, in the plane of section, would have a theoretical boundary bissector (100) planes (similar to a tilt wall boundary made of (100)[001] dislocations) and compare it to the actual boundary. Irespective of the misorientation angle, GB trace are systematically oriented close to the trace of the plane bissector to (100) planes from the two adjacent opx domains (Fig. 12).

R

  the Young's modulus (~100 GPa[START_REF] Duffy | Elasticity of enstatite and its relationship to crystal structure[END_REF])), t the thickness of the folded layer and R the curvature of the folded section. Taking the example in Fig.1, t~100 μm, R~2000 μm, therefore σ~E 20 , i.e. of the order of a few GPa, which is by orders larger than what is recorded in ductilely deformed rocks. This implies that lattice curvature is accommodated plastically by the presence of edge dislocations of a given sign[START_REF] Hirth | Theory of dislocations[END_REF], namely the (100)[001] dislocations in our case. The density ρ of these geometrically necessary dislocations can be related to their Burgers vector (in this case unit cell c) and to the curvature 1 of the lattice along the fold as[START_REF] Nicolas | Crystalline plasticity and solid state flow in metamorphic rocks[END_REF]):

  are not accounted for by the simple crystallographic thermal mismatch model by[START_REF] Van Duysen | Transmission electron microscope study of dislocations in orthopyroxene (Mg, Fe)2Si206[END_REF]. On the contrary, in the flexural slip model, the density of dislocation may vary along a lamellae, as a result of variations in curvature.Furthermore, in Fig.15, we considered the whole cpx lamellae as concentrating the relative slip and the dislocation, irrespective of its lower/upper surface. In reality, depending on how the exsolution proceeds and on which side the lamellae grows, the density of the dislocations on the lamellae upper and lower surfaces may be different (Fig. Supp. Mat. 3). As a consequence, the heterogeneous dislocation distribution at cpx-opx interfaces can only be accounted for when considering our folding model, where dislocations are geometrically necessary to accommodate lattice curvature. On the other hand, one should note that the two models (thermal mismatch and mechanical folding) are not mutually incompatible and the actual distribution of dislocations probably reflects the combination of both processes (Fig. Supp. Mat. 3).

  )). A contrario, in recrystallized zones, where most grains are exsolution-free, GB are likely to be more mobile, which could explain the transition between recrystallization patterns shown in Fig.4: When exsolution lamellae still densely distributed (A), boundaries grow symmetrically by the addition of (100)[001] dislocations and misorientation axes are consistently close to b-axis. In wider recrystallized domains (B), where exsolution lamellae have been expelled, the misorientation axes seem to be randomly distributed, which could be interpreted as the result of combined (i) migration of GB and (ii) potential activation of additional slip systems.It seems therefore that the GB mobility, irrespectively of the misorientation angle, is primarily controlled by the density of exsolution present within opx lattice. The process of expulsion of cpx lamellae from opx is not clear, although it can be unequivocally associated with deformation (Fig. Supp. Mat. 2), as recrystallized grains have a much fewer exsolutions[START_REF] Etheridge | Deformation and recrystallisation of orthopyroxene from the Giles Complex, central Australia[END_REF][START_REF] Raimbourg | Grain size reduction mechanisms and rheological consequences in high-temperature gabbro mylonites of Hidaka[END_REF].The reequilibration of opx, over several microns, by Ca loss, near boundaries of recrystallized grains also supports a relatively high mobility of Ca in opx. This contrasts sharply with[START_REF] Skemer | Effects of solute segregation on the grain-growth kinetics of orthopyroxene with implications for the deformation of the upper mantle[END_REF]'s observation of grain growth inhibition: Calcium, because of its low mobility, is segregated near the GBs and prevents their migration. One may tentatively relate this difference to the static vs dynamic character of the microstructure: In Skemer and Karato (2007)'s static grain growth experiment, only element diffusion is involved, while the Ca-distribution shown in Fig. Supp. Mat. 2

Figure 1 :

 1 Figure 1: Internal folding of an opx porphyroclast (XZ section, with horizontal dextral

Figure 2 :

 2 Figure 2: Optical microscope image in reflected light of the etched surface of the opx

Figure 3 :

 3 Figure 3: (A) Fore-scatter image and (B) EBSD map of homogeneously folded opx.

Figure 4 :

 4 Figure 4: Domains of incipient (A) and extensive (B) recrystallization within opx

Figure 5 :

 5 Figure 5: TEM bright-field images and corresponding diffraction conditions

Figure 6 :

 6 Figure 6: Comparison between measured dislocation spacing and theoretical value for

Figure 7 :

 7 Figure 7: Detail of the dislocations within a straight LAGB with the misorientation

Figure 8 :

 8 Figure 8: (A) TEM bright field image and (B) sketch of a LAGB (thick black line) at

Figure 9 :

 9 Figure 9: (A-B) TEM bright field inages and corresponding diffraction conditions

Figure 10 :

 10 Figure 10: (A) TEM bright field images and corresponding diffraction conditions

Figure 11 :

 11 Figure 11: Structure of a GB with 21° misorientation. (A) Optical reflected light: the

Figure 12 :

 12 Figure 12: Angle, in the section plane, between the trace of the GB and the trace of the

Figure 14 :

 14 Figure 14: Bright-field TEM image showing (100)[001] dislocations (white triangles)

Figure 15 :Figure

 15 Figure 15: (A) Model of opx flexural slip folding, where the equivalent of the folded

Figure

  Figure Supplementary material 2: Gradients in composition associated with

Figure

  Figure Supplementary material 3: Mixed origin of dislocations at cpx-opx interfaces. In

  Fig.3

  Fig.4

  Fig.8

  Fig.11

  Fig.12

  Fig.13

Fig. 14

 14 Fig.14

  Fig.15

  Figure Supplementary material 1

  Figure Supplementary material 2

  Figure Supplementary material 3