Texture Development and Stress–Strain Partitioning in Periclase + Halite Aggregates

Keywords: multiphase deformation; high pressure; texture; plasticity modeling


Introduction
Multiphase materials are ubiquitous in nature and are widely used in engineering due to their desirable mechanical properties. Recently, there has been a growing interest in high-pressure mechanical behavior in many fields such as materials science of super hard materials and minerals physics of the Earth's interior. Materials in the Earth's interior are composed of multiphase materials deformed under high-pressure and -temperature conditions. Thus, it is of great importance to geoscience to understand the deformation behavior of multiphase materials under extreme pressure conditions.
The deformation behavior of individual phases in multiphase aggregates can be quite different compared with a single phase. In multiphase deformation experiments at high pressure, one can utilize synchrotron X-ray diffraction to directly measure lattice strains and estimate the flow stress in each phase. Kaercher et al. (2016) [1] performed in-situ deformation experiments on NaCl + NaMgF 3 two-phase aggregates with a series of phase proportions and reported the flow stress in each phase. The flow stress in the strong NaMgF 3 phase decreases nonlinearly with increasing weak NaCl phase

Experiment Details
High-purity MgO and NaCl powders (MgO: Alfa Aesar, 99.95% purity, stock# 14684, lot# J03Q037; NaCl: Macron Fine Chemicals, 99.0% purity, product# 7581-12, lot# 0000006355) were mixed in varying phase proportions and were ground together in an agate mortar. The samples were then heated in a vacuum oven at −0.08 MPa and 120 • C for 2 weeks to eliminate moisture. Samples were loaded in a boron-kapton gasket [15] with a 50 µm diameter sample chamber. Samples were compressed using a panoramic type DAC with 200 µm culet diamond anvils. Pressure was calibrated with the equation of state of NaCl from Brown (1999) [16]. In order to avoid the NaCl phase transition to the B2 structure, which occurs at~30 GPa (e.g., [17][18][19][20]), samples were compressed to just under 30 GPa. Synchrotron X-ray diffraction in radial geometry was performed at beamline 12.2.2 of the advanced light source (ALS), Lawrence Berkeley National Lab, Berkeley, CA, USA and HPCAT (sector 16), advanced photon source (APS), Argonne National Lab, Lemont, IL, USA. Sample-detector distance and detector non-orthogonality were calibrated using a CeO 2 standard. An X-ray wavelength of 0.49592 Å Minerals 2019, 9,679 4 of 23 at ALS and 0.4066 Å at APS was used. Diffraction images were recorded using a MAR345 detector at ALS and MARCCD detector at APS.

Experiment Data Analysis
Due to axial geometry in DAC experiments, the stress state in the sample can be decomposed into hydrostatic and deviatoric components as: where σ p = 2σ 1 +σ 3 3 is the mean stress and t = σ 3 − σ 1 is the differential stress, which can be used to infer flow stress of the sample. Due to the deviatoric stress on the sample, the measured d-spacings of crystallites are not uniform but vary as a function of the angle between diffraction plane normal and compression axis, χ, as described in [21,22] where: d m is the measured d-spacing from diffraction, d p is the d-spacing under hydrostatic stress, and Q(hkl) describes the deviation of the measured d-spacing from the hydrostatic d-spacing and can be used to represent the average lattice strain of a plane (hkl) (Figure 1). Differential stress t can be used to infer the flow stress of the sample using the following estimation from [22]: where G R (hkl) and G V are shear moduli calculated under the Reuss (iso-stress) and Voigt (iso-strain) approximation. The factor α, which lies between 0 and 1, describe the relative weights of Reuss and Voigt conditions. However, an estimation of t by plasticity modeling of lattice strain and texture evolution is more reliable [23][24][25]. For multiphase deformation, t from Equation (3) can only provide estimation of the stress in each individual phase. In DAC deformation experiments, samples are subjected to both strain and pressure increase simultaneously. Thus, the effects of strain hardening and pressure strengthening are convoluted. However, since strain hardening is small compared to pressure strengthening and strains are relatively low in the DAC, the measured flow strength should be close to the yield stress at the corresponding pressure. Systematic intensity variations along diffraction rings indicate the existence of texture ( Figure 1). Diffraction images are analyzed using the Rietveld method, implemented in the materials analysis using diffraction (MAUD) software [26], to extract phase proportion, lattice strain, and texture information. Q-factors are obtained by fitting d-spacing variation with azimuth. The orientation distribution function (ODF) is computed using the E-WIMV algorithm, which is similar to the WIMV method [27]. The calculated ODF is then exported to the program BEARTEX [28] for smoothing. Refinement of DAC data in MAUD follows the procedure outlined in [29].

Plasticity Modeling
Lattice strain and texture development are modeled using the elasto-viscoplastic self-consistent method (EVPSC) [14] to deduce the deformation mechanisms of each individual phase and to calculate the aggregate yield stress. EVPSC is an effective medium method. This does not explicitly account for microstructural interactions but rather treats single grains as an inclusion in an anisotropic homogenous effective medium (HEM). The properties of the HEM are determined by the average of all inclusions. The plastic behavior of the inclusion at the local level is described by a nonlinear rate-sensitive constitutive law of varying slip systems: where is the reference shear strain rate, is the resolved shear stress of a slip system s at the reference strain rate, is the symmetric Schmid factor for the slip system s, n is an empirical stress exponent, and is the local stress tensor. The Schmid factor describes how stress is resolved onto a slip system. When the stress resolved onto a given slip system is close to the threshold value , deformation will occur on the slip system. This critical resolved shear stress (CRSS) describes the strength of the individual slip systems and to a large extent controls slip systems activity. Due to grain-grain interactions and microstructural effects, CRSS in EVPSC will not be equivalent to single crystal CRSS, but are rather an effective polycrystal CRSS. In the self-consistent scheme, the constitutive equation needs to be linearized. There are several possible schemes, namely tangent, affine, and secant [30]. Affine is used in this study since it allows the choice of different stress exponents for the various phases. At each deformation step, each inclusion interacts with the HEM and macroscopic properties get updated iteratively until the average stress and strain of all inclusions equals the macroscopic stress and strain. Compared with elastoplastic self-consistent (EPSC) [31] and viscoplastic self-consistent (VPSC) [32], both elastic and viscoplastic behavior are considered in EVPSC. Lattice strains (Q-factors) and grain rotation (texture) are modeled together, giving more constraints on slip system activities. The details of EVPSC modeling for high-pressure deformation experiment can be found in [25]. Since strain and pressure increase simultaneously in DAC experiments, it is not possible to separate the pressure and strain hardening effects on CRSS. Both are included in the pressure dependence of the CRSS. In this study, = + / • , where is the

Plasticity Modeling
Lattice strain and texture development are modeled using the elasto-viscoplastic self-consistent method (EVPSC) [14] to deduce the deformation mechanisms of each individual phase and to calculate the aggregate yield stress. EVPSC is an effective medium method. This does not explicitly account for microstructural interactions but rather treats single grains as an inclusion in an anisotropic homogenous effective medium (HEM). The properties of the HEM are determined by the average of all inclusions. The plastic behavior of the inclusion at the local level is described by a non-linear rate-sensitive constitutive law of varying slip systems: where . γ 0 is the reference shear strain rate, τ s is the resolved shear stress of a slip system s at the reference strain rate, m s kl is the symmetric Schmid factor for the slip system s, n is an empirical stress exponent, and σ kl is the local stress tensor. The Schmid factor describes how stress is resolved onto a slip system. When the stress resolved onto a given slip system is close to the threshold value τ s , deformation will occur on the slip system. This critical resolved shear stress (CRSS) describes the strength of the individual slip systems and to a large extent controls slip systems activity. Due to grain-grain interactions and microstructural effects, CRSS in EVPSC will not be equivalent to single crystal CRSS, but are rather an effective polycrystal CRSS. In the self-consistent scheme, the constitutive equation needs to be linearized. There are several possible schemes, namely tangent, affine, and secant [30]. Affine is used in this study since it allows the choice of different stress exponents for the various phases. At each deformation step, each inclusion interacts with the HEM and macroscopic properties get updated iteratively until the average stress and strain of all inclusions equals the macroscopic stress and strain. Compared with elastoplastic self-consistent (EPSC) [31] and viscoplastic self-consistent (VPSC) [32], both elastic and viscoplastic behavior are considered in EVPSC. Lattice strains (Q-factors) and grain rotation (texture) are modeled together, giving more constraints on slip system activities. The details of EVPSC modeling for high-pressure deformation experiment can be found in [25]. Since strain and pressure increase simultaneously in DAC experiments, it is not possible to separate the pressure and strain hardening effects on CRSS. Both are included in the pressure dependence of the CRSS. In this Minerals 2019, 9, 679 6 of 23 study, τ s = τ s 0 + dτ/dP·P, where τ s 0 is the initial CRSS and dτ/dP is the pressure dependence of CRSS. Sample strains in DAC experiments are typically not measured, but a reasonable axial strain value (~24%), which is based on the texture intensity of the pure NaCl sample, is assumed in the EVPSC models of the multiphase samples. Single crystal MgO/NaCl elastic constants and their pressure dependence from Zha et al. (2000) and Whitfield et al. (1976) [33,34] are used for the simulations.

Pure MgO and NaCl
Experiment and modeling results of MgO and NaCl single-phase deformation are shown first in order to demonstrate the intrinsic properties contrast between these two phases. MgO and NaCl share the same f.c.c structure and slip system sets, which are {110} 110 , {100} 011 , and {111} 110 , but with different easiest slip systems. Deformation experiments on single-crystal NaCl at low-pressure and -temperature conditions found that {110} 110 is the easiest slip system followed by {111} 110 and {100} 011 [35]. Likewise, deformation experiments on single-crystal MgO find {110} 110 to be the easiest slip system at ambient conditions but activity of {100} 011 slip increases with elevated temperature (e.g., [36][37][38][39][40][41]). Theoretical calculations indicate that {111} 110 is unlikely to be activated at room temperature (e.g., [42]). Usually, in single-crystal deformation, strain is accommodated by a single slip mode. In polycrystal deformation, back stresses induced at the grain boundary will force multiple slip systems to be activated. Auxiliary systems are also geometrically necessary to accommodate deformation of neighboring grains. Thus, CRSS values used in EVPSC can be much higher than single crystal CRSS due to grain interactions and boundary effects. Single crystal CRSS values will not be directly transferrable to EVPSC to predict polycrystal behaviors. As a result, the CRSS values used in these models are chosen such that simulated lattice strain and texture evolution matches that of the experiment.
The sample composed of pure NaCl powder was compressed up to 29 GPa. Lattice strain of the diffracted planes are in an order of Q(111) > Q(220) > Q(200) (Figure 2a). NaCl gradually develops a (110) compression texture (Figure 3), similar to that observed in compression experiments at 400 • C [43] and consistent with the complementary (111) and (100) texture observed in room temperature extension experiment [44,45]. Based on previous self-consistent modeling results, a high activity of {110} 110 slip will give a (100) compression texture [46] while {100} 011 and {111} 110 give an intense (110) texture [25]. A linear CRSS increase with pressure is given to each slip system in order to fit the Q-factors increase, which is caused mainly by pressure hardening (Table 1). Based on the modeling results, a high activity of {111} 110 slip combined with minor {110} 110 and {100} 011 give the best fit for both texture and Q-factors ( Figure 2b). Although the initial CRSS for {110} 110 used in modeling is almost half of {111} 110 (Table 1), most of the plastic strain is still accommodated by the {111} 110 system, similar to VPSC simulation results for shear deformation of halite [47]. In cubic crystals, there are only six equivalent slip modes for {110} 110 system while there are 12 for {111} 110 . Thus, slip on {110} planes is more likely to be suppressed due to an unfavorable Schmid factor, than slip on {111} planes.
The sample composed of pure MgO powder was compressed up to 47 GPa and modeled using EVPSC in a separate study [48]. Lattice strains are in an order of Q(220) > Q(111) > Q(200) at relative low pressures (<20 GPa). Above 20 GPa, Q(111) increases faster than Q(220) with pressure and gradually surpasses Q(220) (Figure 4). A reordering of Q-factors usually indicates a change of slip system activity. MgO gradually develops an intense (100) texture with increasing pressure ( Figure 5). Based on modeling results, both {110} 110 and {100} 011 are activated with a ratio of~60:40 below 20 GPa. Above 20 GPa, the reordering of Q-factors is caused by increasing activity of the {110} 110 slip system due to pressure hardening of the {100} 011 system. This is fulfilled by giving a steep increase of CRSS with pressure (Table 1). Accompanying this slip system activity change at high pressure is a Minerals 2019, 9, 679 7 of 23 rapid increase in the measured differential stress value. Based on our modeling results on pure MgO and NaCl, the rheology contrast between these two materials is~14.          [48]. Scale bar is given in multiples of random distribution (m.r.d.)

20% MgO + 80% NaCl
The sample of approximately 20% MgO plus 80% NaCl is compressed up to 24 GPa. The lattice strains of NaCl do not change significantly with addition of MgO. Q-factors are still in an order of Q(111) > Q(220) > Q(200) with values that resemble the pure NaCl sample (Figure 6a). However, NaCl gradually develops a (100) texture ( Figure 7) in contrast with the (110) texture observed for pure NaCl. Among the available NaCl slip systems, a (100) compression texture is generated by a high activity of the {110}〈11 0〉 slip system. Thus, it is likely that there is a dominant slip system transition from {111}〈11 0〉 to {110}〈11 0〉 caused by the addition of strong MgO. For MgO, the lattice strains are much lower than in the pure MgO sample. All Q-factors generally increase with pressure ( Figure  6a), which is not observed in the same pressure range in pure MgO (Figure 4a). There is no observable texture development of MgO during compression ( Figure 7). Since dislocation glide leads to grain rotation, a lack of texture development likely indicates that there is little to no plastic deformation in    [48]. Scale bar is given in multiples of random distribution (m.r.d.)

20% MgO + 80% NaCl
The sample of approximately 20% MgO plus 80% NaCl is compressed up to 24 GPa. The lattice strains of NaCl do not change significantly with addition of MgO. Q-factors are still in an order of Q(111) > Q(220) > Q(200) with values that resemble the pure NaCl sample (Figure 6a). However, NaCl gradually develops a (100) texture ( Figure 7) in contrast with the (110) texture observed for pure NaCl. Among the available NaCl slip systems, a (100) compression texture is generated by a high activity of the {110}〈11 0〉 slip system. Thus, it is likely that there is a dominant slip system transition from {111}〈11 0〉 to {110}〈11 0〉 caused by the addition of strong MgO. For MgO, the lattice strains are much lower than in the pure MgO sample. All Q-factors generally increase with pressure ( Figure  6a), which is not observed in the same pressure range in pure MgO (Figure 4a). There is no observable texture development of MgO during compression ( Figure 7). Since dislocation glide leads to grain rotation, a lack of texture development likely indicates that there is little to no plastic deformation in MgO. Given the low Q-factors, the stress in MgO is also likely below the yield stress. Low values of  [48]. Scale bar is given in multiples of random distribution (m.r.d.)

20% MgO + 80% NaCl
The sample of approximately 20% MgO plus 80% NaCl is compressed up to 24 GPa. The lattice strains of NaCl do not change significantly with addition of MgO. Q-factors are still in an order of Q(111) > Q(220) > Q(200) with values that resemble the pure NaCl sample (Figure 6a). However, NaCl gradually develops a (100) texture ( Figure 7) in contrast with the (110) texture observed for pure NaCl. Among the available NaCl slip systems, a (100) compression texture is generated by a high activity of the {110} 110 slip system. Thus, it is likely that there is a dominant slip system transition from {111} 110 to {110} 110 caused by the addition of strong MgO. For MgO, the lattice strains are much lower than in the pure MgO sample. All Q-factors generally increase with pressure ( Figure 6a), which is not observed in the same pressure range in pure MgO (Figure 4a). There is no observable texture development of MgO during compression ( Figure 7). Since dislocation glide leads to grain rotation, a lack of texture development likely indicates that there is little to no plastic deformation in MgO. Given the low Q-factors, the stress in MgO is also likely below the yield stress. Low values of experimental lattice strains and lack of texture development indicate that MgO deforms elastically. For dominant slip on {110} 110 in NaCl/MgO, the spin of symmetric variants (e.g., (110) [-110] and (-110)[110]) cancel each other's rotation, resulting in slow (100) compression texture evolution until strain is large enough to significantly flatten the grains [49]. In our EVPSC simulation, the grain shape is initially flattened in order to yield a faster (100) texture evolution within the~24% axial strain imposed. Based on modeling results, dominant {110} 110 slip combined with auxiliary {100} 011 slip for NaCl and almost no slip for MgO give a good fit to the experiment lattice strain and texture ( Figure 6b).

30% MgO + 70% NaCl
The sample composed of approximately 30% MgO plus 70% NaCl is compressed up to 27 GPa. The Q-factors of NaCl are still in an order of Q(111) > Q(220) > Q(200) (Figure 8a) and slightly higher than the pure NaCl sample. At 4 GPa, where the first diffraction image is taken, NaCl has already developed a (100) texture, corresponding to a dominant {110} 110 slip. As deformation continues, the texture gradually weakens then shifts to (110) (Figure 9), due to activation of the {111} 110 slip system. The Q-factors of MgO are still relatively low compared with pure MgO. However, MgO gradually develops a weak (100) texture ( Figure 9). This indicates that MgO starts to deform plastically in the 30% phase proportion sample. This is near the percolation threshold; thus, MgO-MgO grain contacts are likely to occur in this sample. Experiment textures of MgO/NaCl at 4 GPa are used as starting textures in these models.

50% MgO + 50% NaCl
The sample composed of approximately 50% MgO + 50% NaCl is compressed to 27 GPa. Q-factors of NaCl below 10 GPa are in an order of Q(111) > Q(220) > Q(200) (Figure 10a). Above 10 GPa, the NaCl (222) peak is weak and overlaps with the MgO (220) peak. Thus, we do not report values for Q(111) of NaCl. The Q-factor values are higher than in the 70% NaCl sample. NaCl gradually develops a (110) texture, similar to pure NaCl ( Figure 11). In MgO, the Q-factors are higher than in the 20% and 30% MgO samples but are still lower than pure MgO. MgO gradually develops an obvious (100) texture (Figure 11 for Q(111) of NaCl. The Q-factor values are higher than in the 70% NaCl sample. NaCl gradually develops a (110) texture, similar to pure NaCl ( Figure 11). In MgO, the Q-factors are higher than in the 20% and 30% MgO samples but are still lower than pure MgO. MgO gradually develops an obvious (100) texture (Figure 11

80% MgO + 20% NaCl
The sample of approximately 80% MgO + 20% NaCl is compressed up to 27 GPa. The Q-factors of NaCl are still in an order of Q(111) > Q(220) > Q(200) (Figure 12a). Q(111) is not available above 12 GPa due to weak peak intensity and overlap with the MgO (220) peak, similar to the 50% NaCl sample. The Q-factor values are lower than in the 50% NaCl sample, indicating a decreasing stress in the NaCl phase. This does not follow the trend of higher Q-factor values with increasing MgO phase for Q(111) of NaCl. The Q-factor values are higher than in the 70% NaCl sample. NaCl gradually develops a (110) texture, similar to pure NaCl ( Figure 11). In MgO, the Q-factors are higher than in the 20% and 30% MgO samples but are still lower than pure MgO. MgO gradually develops an obvious (100) texture (Figure 11), indicating a significant amount of plastic strain. Based on modeling results, NaCl has a dominant {111}〈11 0〉 slip with minor auxiliary {110}〈11 0〉 and {100}〈011〉.

80% MgO + 20% NaCl
The sample of approximately 80% MgO + 20% NaCl is compressed up to 27 GPa. The Q-factors of NaCl are still in an order of Q(111) > Q(220) > Q(200) (Figure 12a). Q(111) is not available above 12 GPa due to weak peak intensity and overlap with the MgO (220) peak, similar to the 50% NaCl sample. The Q-factor values are lower than in the 50% NaCl sample, indicating a decreasing stress in the NaCl phase. This does not follow the trend of higher Q-factor values with increasing MgO phase Figure 11. IPFs of experimental and modeled textures of 50% MgO + 50% NaCl. Scale bar is given in multiples of random distribution (m.r.d.).

80% MgO + 20% NaCl
The sample of approximately 80% MgO + 20% NaCl is compressed up to 27 GPa. The Q-factors of NaCl are still in an order of Q(111) > Q(220) > Q(200) (Figure 12a). Q(111) is not available above 12 GPa due to weak peak intensity and overlap with the MgO (220) peak, similar to the 50% NaCl sample. The Q-factor values are lower than in the 50% NaCl sample, indicating a decreasing stress in the NaCl phase. This does not follow the trend of higher Q-factor values with increasing MgO phase proportion in samples with NaCl > 20%. NaCl gradually develops a (110) texture ( Figure 13).  (Figure 12a). A rapid increase and reordering of Q-factors is also observed in pure MgO deformation, which can be explained by "pressure hardening" of {100} 011 slip system of MgO [48]. MgO gradually develops a (100) texture ( Figure 13). Based on modeling results, NaCl has a dominant {111} 110 slip with a total activity of~35%. MgO has a dominant {110} 110 slip with a total activity of~65% (Figure 12b).  (Figure 12a). A rapid increase and reordering of Q-factors is also observed in pure MgO deformation, which can be explained by "pressure hardening" of {100}〈011〉 slip system of MgO [48]. MgO gradually develops a (100) texture ( Figure 13). Based on modeling results, NaCl has a dominant {111}〈11 0〉 slip with a total activity of ~35%. MgO has a dominant {110}〈11 0〉 slip with a total activity of ~65% (Figure 12b).    (Figure 12a). A rapid increase and reordering of Q-factors is also observed in pure MgO deformation, which can be explained by "pressure hardening" of {100}〈011〉 slip system of MgO [48]. MgO gradually develops a (100) texture ( Figure 13). Based on modeling results, NaCl has a dominant {111}〈11 0〉 slip with a total activity of ~35%. MgO has a dominant {110}〈11 0〉 slip with a total activity of ~65% (Figure 12b).

Determination of Stress and Strain
Stress and strain behavior in the individual NaCl/MgO phases is quite different under varying phase proportions. Knowing stress and strain of the individual phases can provide insight into the deformation behavior of multiphase aggregates under varying microstructures. For single-phase aggregates, the measured average strain is equal to the macroscopic strain imposed on the sample. For multiphase deformation, the average strain of each phase and its contribution to bulk strain can be derived from activity results obtained from modeling. Direct stress measurements on the bulk multiphase aggregate are not feasible in high-pressure experiments. Usually, stress is evaluated based on lattice strain attained from diffraction using Equation (3) or from sample thickness and pressure gradient measured from ruby fluorescence [50]. However, measurements of differential stress from different studies show inconsistencies. Previous studies on NaCl polycrystals report different differential stresses in a wide pressure range ( Figure 14) [51][52][53][54][55][56][57]. Much of this inconsistency comes from different assumptions of mechanical behavior of polycrystals (iso-stress, iso-strain, etc.), and the choice of Q-factor to represent the average lattice strain of the sample can also bias the results. The differential stress of NaCl based on measured lattice strains in this study for different mechanical assumptions is calculated using Equation (3) (Figure 14). Given the uncertainty of differential stress measurements, it is important to adopt a unified method in order to compare different experiments in this study. The results from EVPSC are consistent and avoid issues due to lattice plane sampling bias. The self-consistent scheme accounts for strain/stress heterogeneity, which is closer to the real deformation state of aggregates. In order to study the stress relationship between aggregates and individual phases, we use EVPSC results for stresses in the individual phases and for the stresses in the multiphase.

Determination of Stress and Strain
Stress and strain behavior in the individual NaCl/MgO phases is quite different under varying phase proportions. Knowing stress and strain of the individual phases can provide insight into the deformation behavior of multiphase aggregates under varying microstructures. For single-phase aggregates, the measured average strain is equal to the macroscopic strain imposed on the sample. For multiphase deformation, the average strain of each phase and its contribution to bulk strain can be derived from activity results obtained from modeling. Direct stress measurements on the bulk multiphase aggregate are not feasible in high-pressure experiments. Usually, stress is evaluated based on lattice strain attained from diffraction using Equation (3) or from sample thickness and pressure gradient measured from ruby fluorescence [50]. However, measurements of differential stress from different studies show inconsistencies. Previous studies on NaCl polycrystals report different differential stresses in a wide pressure range ( Figure 14) [51][52][53][54][55][56][57]. Much of this inconsistency comes from different assumptions of mechanical behavior of polycrystals (iso-stress, iso-strain, etc.), and the choice of Q-factor to represent the average lattice strain of the sample can also bias the results. The differential stress of NaCl based on measured lattice strains in this study for different mechanical assumptions is calculated using Equation (3) (Figure 14). Given the uncertainty of differential stress measurements, it is important to adopt a unified method in order to compare different experiments in this study. The results from EVPSC are consistent and avoid issues due to lattice plane sampling bias. The self-consistent scheme accounts for strain/stress heterogeneity, which is closer to the real deformation state of aggregates. In order to study the stress relationship between aggregates and individual phases, we use EVPSC results for stresses in the individual phases and for the stresses in the multiphase.

Stress and Strain of Individual Phase
In EVPSC, the aggregates plastic strain is the sum of the volume-averaged plastic strain of each phase. Total activity indicates the contribution from corresponding phases to the aggregate. The ratio of activity to phase proportion indicates the average strain of the corresponding phase relative to the aggregate. Generally, activity of MgO and NaCl both decrease with increasing volume fraction of the counterpart phase ( Figure 15). This trend is not obvious when the volume fraction of NaCl is above 70% (MgO below 30%). In this range, NaCl likely forms an IWL microstructure and most of the plastic strain is accommodated by the weak NaCl matrix, leaving hard MgO undeformed. When the volume fraction of MgO is above 30%, passing the percolation threshold, MgO starts to deform plastically. Partial plastic strain is accommodated by MgO, resulting in decreasing activity in NaCl. However,

Stress and Strain of Individual Phase
In EVPSC, the aggregates plastic strain is the sum of the volume-averaged plastic strain of each phase. Total activity indicates the contribution from corresponding phases to the aggregate. The ratio of activity to phase proportion indicates the average strain of the corresponding phase relative to the aggregate. Generally, activity of MgO and NaCl both decrease with increasing volume fraction of the counterpart phase ( Figure 15). This trend is not obvious when the volume fraction of NaCl is above 70% (MgO below 30%). In this range, NaCl likely forms an IWL microstructure and most of the plastic strain is accommodated by the weak NaCl matrix, leaving hard MgO undeformed. When the volume fraction of MgO is above 30%, passing the percolation threshold, MgO starts to deform plastically. Partial plastic strain is accommodated by MgO, resulting in decreasing activity in NaCl. However, average strain of MgO and NaCl show different trends with decreasing phase proportion ( Figure 15). Average strain of NaCl is larger than the average strain of the aggregate and it increases with decreasing phase proportions. NaCl strain increases from 100% in pure NaCl to 165% in the 20% NaCl sample. NaCl grains tend to deform more with a lower phase proportion. In contrast, average strain of MgO decreases with decreasing phase proportion and is always lower than the aggregate strain. MgO strain decreases from 100% in pure MgO to nearly 0 in the 20% MgO sample.  (Figure 15). Average strain of NaCl is larger than the average strain of the aggregate and it increases with decreasing phase proportions. NaCl strain increases from 100% in pure NaCl to 165% in the 20% NaCl sample. NaCl grains tend to deform more with a lower phase proportion. In contrast, average strain of MgO decreases with decreasing phase proportion and is always lower than the aggregate strain. MgO strain decreases from 100% in pure MgO to nearly 0 in the 20% MgO sample. To better show the stress relationship between individual phases and the aggregate, the average differential stress of MgO, NaCl, and the aggregate under different volume fractions is plotted in Figure 16a-d. When NaCl is the dominant phase and an IWL microstructure is expected ( Figure  16a,b), aggregate stress is low and close to the stress of NaCl since most stress is concentrated in the weak NaCl matrix. MgO is below its yield point and the stress increase of MgO with pressure is an elastic response to small strain imposed by the NaCl matrix. When NaCl and MgO have comparable phase proportions (Figure 16c), both NaCl and MgO percolate the other phase and both deform plastically. Stress in MgO reaches ~2 GPa early in deformation and does not change much with continued compression. Aggregate stress is intermediate to stresses in NaCl and MgO. When MgO is the dominant phase and an LBF microstructure is expected, aggregate stress is much higher than NaCl and close to MgO since most of the stress is concentrated in the strong MgO framework ( Figure  16d). Aggregates stress increases abruptly starting at ~13 GPa due to pressure hardening of MgO while stress in NaCl is not affected. Similar behavior is also observed in the pure MgO sample but not in the 50% MgO sample where MgO also deforms plastically. The activity of MgO is only 25% in the 50% MgO sample, among which only 2-3% is from the {100}〈011〉 system (Figure 10b). If high pressure makes {100}〈011〉 harder to activate, the activity change will be too small to have an effect on stress in the 50% MgO sample. 'Pressure hardening' of the {100}〈011〉 system of MgO will not necessarily activate the harder {110}〈11 0〉 system of MgO. In the presence of the weak NaCl phase, strain can be accommodated by slip in NaCl rather than the harder {110}〈11 0〉 system of MgO. However, for pure MgO or the 80% MgO sample, which likely has an LBF microstructure, large amounts of strain must be accommodated by the MgO framework. An increased activation of {110}〈11 0〉 slip in MgO is required when {100}〈011〉 is hard to activate in order to accommodate geometrically required strain.
Stress of NaCl initially increases with increasing fractions of MgO but reaches a peak value around 50% MgO and then decreases with increasing MgO (Figure 16e). Stress in MgO generally decreases with increasing amounts of NaCl but not in a linear fashion (Figure 16f). This difference is caused by two things. First, due to a high yield strength contrast between NaCl and MgO, MgO does not deform plastically when its phase proportion is low but NaCl always deforms plastically. Second, MgO can obstruct plastic flow of NaCl at phase boundaries while NaCl cannot since MgO is much To better show the stress relationship between individual phases and the aggregate, the average differential stress of MgO, NaCl, and the aggregate under different volume fractions is plotted in Figure 16a-d. When NaCl is the dominant phase and an IWL microstructure is expected (Figure 16a,b), aggregate stress is low and close to the stress of NaCl since most stress is concentrated in the weak NaCl matrix. MgO is below its yield point and the stress increase of MgO with pressure is an elastic response to small strain imposed by the NaCl matrix. When NaCl and MgO have comparable phase proportions (Figure 16c), both NaCl and MgO percolate the other phase and both deform plastically. Stress in MgO reaches~2 GPa early in deformation and does not change much with continued compression. Aggregate stress is intermediate to stresses in NaCl and MgO. When MgO is the dominant phase and an LBF microstructure is expected, aggregate stress is much higher than NaCl and close to MgO since most of the stress is concentrated in the strong MgO framework (Figure 16d). Aggregates stress increases abruptly starting at~13 GPa due to pressure hardening of MgO while stress in NaCl is not affected. Similar behavior is also observed in the pure MgO sample but not in the 50% MgO sample where MgO also deforms plastically. The activity of MgO is only 25% in the 50% MgO sample, among which only 2-3% is from the {100} 011 system (Figure 10b). If high pressure makes {100} 011 harder to activate, the activity change will be too small to have an effect on stress in the 50% MgO sample. 'Pressure hardening' of the {100} 011 system of MgO will not necessarily activate the harder {110} 110 system of MgO. In the presence of the weak NaCl phase, strain can be accommodated by slip in NaCl rather than the harder {110} 110 system of MgO. However, for pure MgO or the 80% MgO sample, which likely has an LBF microstructure, large amounts of strain must be accommodated by the MgO framework. An increased activation of {110} 110 slip in MgO is required when {100} 011 is hard to activate in order to accommodate geometrically required strain. remain below its yield point and, thus, are lower compared to samples where MgO is interconnected and must deform plastically. Stress of NaCl initially increases with increasing fractions of MgO but reaches a peak value around 50% MgO and then decreases with increasing MgO (Figure 16e). Stress in MgO generally decreases with increasing amounts of NaCl but not in a linear fashion (Figure 16f). This difference is caused by two things. First, due to a high yield strength contrast between NaCl and MgO, MgO does not deform plastically when its phase proportion is low but NaCl always deforms plastically. Second, MgO can obstruct plastic flow of NaCl at phase boundaries while NaCl cannot since MgO is much stronger than NaCl. When both MgO and NaCl percolate and deform plastically, NaCl cannot impinge into MgO grains and instead must flow around harder MgO. On the other hand, MgO is hard enough to impose deformation on NaCl. Thus, stress in NaCl is related to obstruction of flow of the NaCl interconnected network by MgO. Stress of MgO is more affected by the local environment surrounding each MgO grain. Deviatoric strain (and strain rate) of NaCl and MgO under varying volume fractions is different ( Figure 15). To clarify the effect of phase proportions on stress, the stresses of pure NaCl and MgO are normalized to the strain rate of each phase proportion using the rate-sensitive power law equation (dashed line in Figure 16e,f). This can then be compared with 'true' experiment stress. For stress of NaCl, when phase proportions of MgO are below 30%, the rate normalized and experimental stresses of NaCl show little variance except for 20% MgO at high pressure. This might be due to changing deformation mechanisms or large uncertainty from fluctuation (experimental error) of Q-factors ( Figure 6a). Within this range, almost all the strain is accommodated by plastic flow of the weak NaCl matrix around MgO inclusions. MgO does not obstruct NaCl flow due to wide dispersion of MgO grains in the weak matrix and because most grain to grain interactions are NaCl-NaCl. The effective polycrystal CRSS of NaCl does not change significantly compared to pure NaCl. When the phase proportion of MgO is 50%, experimental stress is higher than the rate normalized stress. At 50% MgO, NaCl flow is obstructed at MgO-NaCl boundaries and NaCl must flow through the MgO matrix, resulting in high stress concentration in the NaCl phase. When the phase proportion of MgO is 80%, experimental stress drops and is lower than the rate normalized stress. At this phase proportion, an LBF microstructure is expected and NaCl grains are isolated from each other. NaCl grains form weak pockets protected by a hard MgO framework. Most of the strain and stress are accommodated by the MgO framework. Stress in NaCl is far lower than is required to displace the MgO framework. NaCl thus deforms passively in response to the evolution of the MgO framework to maintain compatibility at phase boundaries. Thus, there is no flow of NaCl through the MgO matrix and NaCl inclusions are passive and deform in response to the LBF, resulting in lower stresses compared to pure NaCl.
When studying stress concentration in the weak matrix as a function of the strong phase proportion, it is helpful to define an α factor to describe the average spacing between inclusions [3]: where d is the average diameter of the inclusions and z is the average center distance between two inclusions. At a given inclusion size, a higher phase proportion of the weak phase results in a larger value of α. A value of α = 1 corresponds to 100% weak phase and α approaches 0 when inclusions touch to form an LBF microstructure. Handy (1990) [3] examined the flow stress in quartz feldspar mylonites as a function of α. For quartz, dynamic recrystallization grain size is inversely related with flow stress. By examining quartz grain size between several feldspar inclusion pairs in natural rocks, grain size distribution can be plotted as a function of α ([3], Figure 7b). Quartz grain size is constant with an α value in the range of 0.4-1.0, indicating stress of the quartz matrix is not sensitive to the feldspar phase proportion in this range. Quartz grain size decreases with lower α values in the range of 0.15 to 0.4, indicating stress in the quartz matrix increases with increased feldspar phase proportion. This is consistent with results in this study that show that the stress of NaCl does not increase until MgO phase proportions reach values between 30% and 50%. When MgO reaches the percolation threshold (~30-40%), LBF type behavior is observed, similar to that observed in stress percolation model shown in [58].
For MgO, when the NaCl is 20%, stress decreases significantly compared with pure MgO and is also much lower than the rate normalized stress. Since NaCl is much weaker than MgO, NaCl-MgO boundaries have little effect on obstruction to MgO flow compared with MgO-MgO boundaries. When parts of the MgO-MgO boundaries are replaced by MgO-NaCl boundaries, the overall obstruction to MgO flow decreases, resulting in lower stress compared with pure MgO. When the phase proportion of NaCl is 50%, stress does not change significantly compared to 20% NaCl. When NaCl reaches the percolation threshold, the percentage of MgO-NaCl boundaries in a local grain environment does not change significantly with added NaCl. Thus, stress of MgO is not sensitive to phase proportions of NaCl within this range. However, when the phase proportion of NaCl is above 70%, MgO no longer percolates. Most of the strain and stress are accommodated by NaCl. NaCl controls the aggregate stress and stress levels are too low to induce yielding in MgO. Stresses in MgO remain below its yield point and, thus, are lower compared to samples where MgO is interconnected and must deform plastically.
Based on these experimental results and discussion, simple schemes for stress evolution in NaCl and MgO with respect to phase proportions can be constructed (Figure 17). Stresses in NaCl do not increase until MgO reaches the percolation threshold and begins to deform plastically. Within the phase proportion range where both MgO and NaCl percolate, stress in NaCl increases with increased MgO due to obstruction of NaCl flow by the MgO. When NaCl is below the percolation threshold and an MgO LBF framework is expected, stress in NaCl decreases with increased MgO due to shielding of NaCl grains by the MgO framework and the passive nature of the isolated weak pockets. In MgO, stress initially decreases with increased NaCl due to more freedom for MgO to flow by increasing MgO-NaCl boundaries. Once NaCl starts to percolate MgO, the local environment of MgO grains does not change with added NaCl. Thus, stress of MgO does not decrease until the volume of MgO falls below the percolation threshold. When MgO is no longer interconnected and is surrounded by NaCl, the stress in MgO decreases due to a transition from plastic to elastic deformation. Based on these experimental results and discussion, simple schemes for stress evolution in NaCl and MgO with respect to phase proportions can be constructed ( Figure 17). Stresses in NaCl do not increase until MgO reaches the percolation threshold and begins to deform plastically. Within the phase proportion range where both MgO and NaCl percolate, stress in NaCl increases with increased MgO due to obstruction of NaCl flow by the MgO. When NaCl is below the percolation threshold and an MgO LBF framework is expected, stress in NaCl decreases with increased MgO due to shielding of NaCl grains by the MgO framework and the passive nature of the isolated weak pockets. In MgO, stress initially decreases with increased NaCl due to more freedom for MgO to flow by increasing MgO-NaCl boundaries. Once NaCl starts to percolate MgO, the local environment of MgO grains does not change with added NaCl. Thus, stress of MgO does not decrease until the volume of MgO falls below the percolation threshold. When MgO is no longer interconnected and is surrounded by NaCl, the stress in MgO decreases due to a transition from plastic to elastic deformation.

Aggregates Stress
Aggregate stress generally decreases nonlinearly as the weak NaCl phase proportion increases ( Figure 18). The addition of only 20% NaCl lowers the aggregate stress by 50% compared to pure MgO. With additional NaCl, aggregates stress decreases but with a more gradual slope. When NaCl dominates the aggregate's composition, aggregate stress is not sensitive to added NaCl. Aggregate stress under the IWL and LBF microstructures are also estimated using equations from [4], which can be treated as upper and lower bounds for true aggregate stress (solid and dotted lines in Figure 18). For materials like MgO and NaCl, which have a high yield strength contrast, LBF stress is reasonable only in a very narrow range of phase proportions that lie close to pure MgO. IWL stress is close to the aggregate stress when the NaCl phase proportion is over 70%, which is where MgO no longer percolates through the NaCl matrix. At intermediate compositions, aggregate stress lies between LBF and IWL bounds. At low NaCl phase proportions, the value is closer to the LBF end member and at high NaCl phase proportions it is closer to the IWL bound. Aggregate stress for low and high NaCl

Aggregates Stress
Aggregate stress generally decreases nonlinearly as the weak NaCl phase proportion increases ( Figure 18). The addition of only 20% NaCl lowers the aggregate stress by 50% compared to pure MgO. With additional NaCl, aggregates stress decreases but with a more gradual slope. When NaCl dominates the aggregate's composition, aggregate stress is not sensitive to added NaCl. Aggregate stress under the IWL and LBF microstructures are also estimated using equations from [4], which can be treated as upper and lower bounds for true aggregate stress (solid and dotted lines in Figure 18). For materials like MgO and NaCl, which have a high yield strength contrast, LBF stress is reasonable only in a very narrow range of phase proportions that lie close to pure MgO. IWL stress is close to the aggregate stress when the NaCl phase proportion is over 70%, which is where MgO no longer percolates through the NaCl matrix. At intermediate compositions, aggregate stress lies between LBF and IWL bounds. At low NaCl phase proportions, the value is closer to the LBF end member and at high NaCl phase proportions it is closer to the IWL bound. Aggregate stress for low and high NaCl phase proportions in this study is in good agreement with predicted values under similar yield contrast, as given in [3]. At intermediate phase proportions, the aggregate stress of Handy is lower than that of this study. This may be due to foliation of the weak phase in the highly strained geologic samples studied by Handy (1990) [3].

Texture Evolution with Varying Phase Proportions
For MgO, when the phase proportion of MgO is below 30%, most of the strain is accommodated by the weak NaCl matrix. Stresses in the aggregate are too low for MgO to deform plastically, resulting in little to no texture development. When the phase proportion of MgO is above 30%, MgO starts to percolate the NaCl matrix and stress in MgO become high enough for plastic deformation. This results in a (100) texture development due to {110}〈11 0〉 slip similar to pure MgO. Texture strength of MgO increases with increased MgO phase proportions due to an increase in plastic strain in the MgO phase.
In these samples, NaCl always develops a texture regardless of the phase proportion since it is weak and aggregate stresses are always high enough to induce plastic deformation. However, deformation experiments on NaMgF3 + NaCl aggregates [1], which have a similar yield strength contrast as MgO + NaCl, shows that the weak NaCl does not develop obvious texture even if the phase proportion of NaCl is as high as 50%. Deformation experiments on CaGeO3 + MgO aggregates [10], where MgO is the weak phase, shows that MgO develops a (100) texture even if the MgO phase proportion is as low as 28% with an LBF microstructure. In contrast, deformation experiments on bridgmanite + ferropericlase ((Fe, Mg)O) aggregates [2] under similar phase proportion show that weak ferropericlase does not develop texture. Given these results and the observation in this study, plastic anisotropy may play an important role on texture development in multiphase aggregates. NaMgF3/bridgmanite have low-symmetry orthorhombic structures and have higher plastic anisotropy than MgO/CaGeO3, which have high-symmetry cubic or pseudo-cubic structures. A high plastic anisotropy in the strong matrix is likely to induce local heterogeneity of strain and stress in the weak phase, causing a disruption of coherent texture development.
The texture pattern of NaCl can also vary depending of the phase proportion of MgO. When the MgO phase proportion is 20%, NaCl develops a (100) texture in contrast to the (110) texture observed in the pure NaCl sample. When dealing with multiphase deformation, whether individual phases 'percolate' (are continuous throughout the sample) or are isolated as inclusions in a matrix greatly affects the material behavior. When the strong phase is isolated as inclusions in the weak matrix, it will cause large local heterogeneity of stress and strain, which in turn affects the flow pattern of the weak phase [59]. In this case, the volume fraction of MgO is under the percolation threshold (~30-

Texture Evolution with Varying Phase Proportions
For MgO, when the phase proportion of MgO is below 30%, most of the strain is accommodated by the weak NaCl matrix. Stresses in the aggregate are too low for MgO to deform plastically, resulting in little to no texture development. When the phase proportion of MgO is above 30%, MgO starts to percolate the NaCl matrix and stress in MgO become high enough for plastic deformation. This results in a (100) texture development due to {110} 110 slip similar to pure MgO. Texture strength of MgO increases with increased MgO phase proportions due to an increase in plastic strain in the MgO phase.
In these samples, NaCl always develops a texture regardless of the phase proportion since it is weak and aggregate stresses are always high enough to induce plastic deformation. However, deformation experiments on NaMgF 3 + NaCl aggregates [1], which have a similar yield strength contrast as MgO + NaCl, shows that the weak NaCl does not develop obvious texture even if the phase proportion of NaCl is as high as 50%. Deformation experiments on CaGeO 3 + MgO aggregates [10], where MgO is the weak phase, shows that MgO develops a (100) texture even if the MgO phase proportion is as low as 28% with an LBF microstructure. In contrast, deformation experiments on bridgmanite + ferropericlase ((Fe, Mg)O) aggregates [2] under similar phase proportion show that weak ferropericlase does not develop texture. Given these results and the observation in this study, plastic anisotropy may play an important role on texture development in multiphase aggregates. NaMgF 3 /bridgmanite have low-symmetry orthorhombic structures and have higher plastic anisotropy than MgO/CaGeO 3 , which have high-symmetry cubic or pseudo-cubic structures. A high plastic anisotropy in the strong matrix is likely to induce local heterogeneity of strain and stress in the weak phase, causing a disruption of coherent texture development.
The texture pattern of NaCl can also vary depending of the phase proportion of MgO. When the MgO phase proportion is 20%, NaCl develops a (100) texture in contrast to the (110) texture observed in the pure NaCl sample. When dealing with multiphase deformation, whether individual phases 'percolate' (are continuous throughout the sample) or are isolated as inclusions in a matrix greatly affects the material behavior. When the strong phase is isolated as inclusions in the weak matrix, it will cause large local heterogeneity of stress and strain, which in turn affects the flow pattern of the weak phase [59]. In this case, the volume fraction of MgO is under the percolation threshold (~30-40%) and little to no flow occurs in MgO. Thus, MgO behaves like strong inclusions in a flowing weak NaCl matrix. In pure NaCl, {111} 110 slip is geometrically favored since this slip system has most symmetric variants, even though {110} 110 is the so-called 'easiest' slip system. With the presence of hard MgO inclusions, plastic flow by {111} 110 is geometrically hampered, forcing activation of the softer but geometrically less favored {110} 110 slip system. This results in development of a (100) NaCl texture in the 20% MgO sample (Figure 19). It is worth noting that a dominant slip transition from {111} 110 to {110} 110 is fulfilled only by increasing the CRSS value of {111} 110 slip to suppress its activation without greatly modifying the CRSS of other slip systems from those of pure NaCl CRSS (Table 1). A similar phenomenon is also observed in Ag and Ag + Fe composite, where it is assumed that the texture of weak Ag is affected by the strong Fe [13].  (Table 1). A similar phenomenon is also observed in Ag and Ag + Fe composite, where it is assumed that the texture of weak Ag is affected by the strong Fe [13]. When the MgO phase proportion is 30%, close to the percolation threshold, a texture transition from (100) to (110) of NaCl can be observed. NaCl initially develops a (100) texture maximum, resembling the case were MgO does not percolate, but gradually transitions to a (110) maximum as deformation proceeds, similar to when MgO percolates NaCl (Figure 19). This likely corresponds to a microstructure evolution in the sample from isolated MgO inclusions towards percolation of MgO as relatively undeformed MgO grains are brought into contact with each other when the sample flattens. Local heterogeneity of strain and stress decrease in NaCl as an LBF develops and begins to impose homogeneous strain in the aggregate. When the MgO phase proportion exceeds 50%, MgO is interconnected and fully percolates the sample. The LBF framework of MgO deforms plastically and controls deformation and imposes a more homogeneous strain in the sample. This causes the dominant slip system of NaCl to switch back to {111}〈11 0〉, resulting in a (110) texture maximum ( Figure 19).

Conclusions
We performed several DAC deformation experiments on NaCl + MgO aggregates with a series of phase proportions from pure NaCl to MgO with the purpose of systematically studying texture and mechanical behavior of multiphase materials with high yield strength contrast components at high pressure. NaCl always deforms plastically regardless of the MgO phase proportion while MgO only deforms plastically when it percolates through the NaCl matrix. When deformed plastically, MgO develops a consistent (100) compression texture similar to pure MgO. NaCl develops a surprisingly different (100) compression texture in the presence of 20% MgO compared with a typical (110) texture. This may be due to changing deformation mechanisms introduced by MgO inclusions and resultant strain and stress heterogeneity. Based on experiments and modeling results in this study, we obtain a scheme for stress evolution of the aggregate and individual phases as a function of phase proportions. As shown in this study, the EVPSC approach is a valuable tool for predicting texture development and also extracting strain and stress partitioning in individual phases and in the aggregate. The latter is currently not feasible to measure directly in high-pressure multiphase deformation experiments. DAC deformation experiments have the disadvantages that phase proportion and microstructure are not controlled precisely, and samples are not easy to recover for microstructural analysis. In the future, with the help of large volume deformation (e.g., the deformation-DIA), additional deformation experiments on varying yield strength contrast aggregates with controllable phase proportions and microstructures should be done to fully constrain the mechanical behavior and texture development in multiphase materials over a range of pressures and temperature conditions. When the MgO phase proportion is 30%, close to the percolation threshold, a texture transition from (100) to (110) of NaCl can be observed. NaCl initially develops a (100) texture maximum, resembling the case were MgO does not percolate, but gradually transitions to a (110) maximum as deformation proceeds, similar to when MgO percolates NaCl ( Figure 19). This likely corresponds to a microstructure evolution in the sample from isolated MgO inclusions towards percolation of MgO as relatively undeformed MgO grains are brought into contact with each other when the sample flattens. Local heterogeneity of strain and stress decrease in NaCl as an LBF develops and begins to impose homogeneous strain in the aggregate. When the MgO phase proportion exceeds 50%, MgO is interconnected and fully percolates the sample. The LBF framework of MgO deforms plastically and controls deformation and imposes a more homogeneous strain in the sample. This causes the dominant slip system of NaCl to switch back to {111} 110 , resulting in a (110) texture maximum ( Figure 19).

Conclusions
We performed several DAC deformation experiments on NaCl + MgO aggregates with a series of phase proportions from pure NaCl to MgO with the purpose of systematically studying texture and mechanical behavior of multiphase materials with high yield strength contrast components at high pressure. NaCl always deforms plastically regardless of the MgO phase proportion while MgO only deforms plastically when it percolates through the NaCl matrix. When deformed plastically, MgO develops a consistent (100) compression texture similar to pure MgO. NaCl develops a surprisingly different (100) compression texture in the presence of 20% MgO compared with a typical (110) texture. This may be due to changing deformation mechanisms introduced by MgO inclusions and resultant strain and stress heterogeneity. Based on experiments and modeling results in this study, we obtain a Minerals 2019, 9, 679 20 of 23 scheme for stress evolution of the aggregate and individual phases as a function of phase proportions. As shown in this study, the EVPSC approach is a valuable tool for predicting texture development and also extracting strain and stress partitioning in individual phases and in the aggregate. The latter is currently not feasible to measure directly in high-pressure multiphase deformation experiments. DAC deformation experiments have the disadvantages that phase proportion and microstructure are not controlled precisely, and samples are not easy to recover for microstructural analysis. In the future, with the help of large volume deformation (e.g., the deformation-DIA), additional deformation experiments on varying yield strength contrast aggregates with controllable phase proportions and microstructures should be done to fully constrain the mechanical behavior and texture development in multiphase materials over a range of pressures and temperature conditions.