## The Effect of Chemical Disorder on Defect Formation and Migration in Disordered MAX Phases

Prashant Singh, a* Daniel Sauceda,  and Raymundo Arroyave a a a Department of Materials Science &amp; Engineering, Texas A&amp;M University, College Station, TX 77843, USA

## Abstract

MAX phases have attracted increased attention due to their unique combination of ceramic and metallic properties. Point-defects are known to play a vital role in the structural, electronic and transport properties of alloys in general and this system in particular. As some MAX phases have been shown to be stable in non-stoichiometric compositions, it is likely that such alloying effects will  affect  the  behavior  of  lattice  point  defects.  This  problem,  however,  remains  relatively unexplored. In this work, we investigate the alloying effect on the structural-stability, energystability, electronic-structure, and diffusion barrier of point defects in MAX phase alloys within a first-principles density functional theory framework. The vacancy (VM, VA, VX) and antisite (M-A; M-X) defects are considered with M and A site disorder in (Zr-M)2 (AA')C, where M=Cr,Nb,Ti and AA'=Al,  Al -Sn,  Pb-Bi.  Our  calculations  suggest  that  the  chemical  disorder  helps  lower  the  VA formation energies compared to VM and VX. The VA diffusion barrier is also significantly reduced for M-site disorder compared to their ordered counterpart. This is very important finding because reduced barrier height will ease the Al diffusion at high-operating temperatures, which will help the formation of passivating oxide layer (i.e., Al2O3 in aluminum-based MAX phases) and will slow down or stop the material degradation. We believe that our study will provide a fundamental understanding and an approach to tailor the key properties that can lead to the discovery of new MAX phases.

Keywords: Disorder, MAX phase, Point-defects, Vacancy-diffusion, Density-functional theory

## Introduction

Research  in  high-temperature  structural  materials  is  gaining  momentum  due  to  the  ongoing quest for ever increasing operating efficiencies in power generation technologies [1]. An example is the nuclear power industry, which seeks new reactor designs that in turn demand structural and  coating  materials  capable  of  withstanding  higher  operating  temperatures  at  higher irradiation  conditions.  MAX  phases  are  a  new  class  of  high-temperature  materials  with  the general formula M n+1 AX n  [n=1,3; M=transition element; A=A-group, and X=C or N]. Their hightemperature thermodynamic and structural stability, damage-tolerance and oxidation-resistance [2-5]  make MAX-phases a potential candidate structural material for nuclear applications [6]. Beyond structural and coating applications, research on MAX phases has also led to the discovery of MXenes based on the selective leaching of the A element [7,8]. Many MAX phase (and MXene) properties of interest in these materials depend on the nature and behavior of point defects within the lattice.

In  the  past,  irradiation  has  been  shown  to  significantly  alter  the  microstructures  of  systems through the creation of point defects [9-12]. For example, in the context of nuclear applications, irradiation degrades the material by damage accumulation.  These defects interact with other defects and atoms within the solid-solution while they migrate through the lattice [13]. It is thus to be expected that understanding the behavior of (extrinsic or intrinsic) point defects in MAX phases  is  important  if  one  is  to  understand  the  factors  controlling  their  performance  under irradiation environments. In many materials systems, chemical disorder-induced alloy complexity is  another  factor  that  can  significantly  impact  the  defect  evolution,  energy  dissipation,  and radiation resistance [14-17] that arises from the interaction of (energetic) ions with solids in a radiation environment.

Following the trends observed in the development of other materials systems, recent efforts on the synthesis, characterization and modeling of MAX phases has focused on the expansion of the composition palette through mixing in the different sublattices of the MAX phase crystal system. While some work has been carried out already in understanding the energetics and behavior of point defects in stoichiometric MAX phases [12], the investigation of defect physics of disordered MAX  phases  remains  relatively  unexplored  relative  to  other  materials  systems.  It  becomes desirable  to  understand  the  interplay  of  disorder  and  the  point  defects  to  develop  an understanding  towards  the  microstructural  changes  due  to  irradiation  for  the  future  nuclear structural  materials  [18-27].  Recent  studies  show  that  disorder  can  be  used  to  tailor  alloy properties  to  achieve  desirable  physical  and  chemical  properties,  e.g.,  thermoelectric-effects [28],  charge-transport  [29],  and  magnetoresistance  [30].  The  tailoring  of  properties  in  these applications ultimately arises from the strong connection between chemical disorder and the stability of point defects in those materials systems. It is thus to be expected that, in a similar fashion, chemical disorder in MAX phase crystal systems plays an important role in determining the local configuration, formation energetics and migration behavior of point defects in these systems.

In this work, we discuss the impact of chemical disorder on the behavior of point defects of MAX phases. Five MAX phase solid solution systems [based on the stoichiometry M2AlC, M=Ti, Zr, Nb, Cr,  Ta  (for  Ta  see  supplement)],  of  the  211  stoichiometry,  are  chosen  to  investigate  the connection between alloying and point defect behavior. We choose to investigate the effect of single-site (M-, or A-), and two-site (M- &amp; A-) disorder on the energetics and migration behavior of point defects. Our calculations show, for example, that modifying chemical disorder leads to notable  changes  on  the  defect  formation  and  migration  energies.  Our  focus  remains  on  Alcontaining systems, which arises from their superior oxidation resistance that results from the formation (in some of these compounds) of stable alumina-based passivating layers that protect them from further oxidation [31]. Specifically, we apply first-principles density functional theory and climbing nudge elastic band (cNEB) schemes to study VM, VA and VC vacancymediated diffusion in (Zr-M)2AlC and (Zr-M)2 (AA')C , where M=Cr,Nb,Ti and AA'=Al, Al -Sn, Pb-Bi. Zr-containing MAX phases were selected because they have already been identified as candidate materials for nuclear structural applications. The detailed discussion about the dependence of the migration energy on the alloying elements and defects is provided. We also investigate the responsible electronic structure features that control the properties of defects and its connection to elemental properties, such as, atomic-radii and electronegativity.

## Computational details

While  in  the  context  of  Density  Functional  Theory  (DFT)  there  are  many  approaches  used nowadays to account for chemical disorder in a crystal structure, the use of periodic supercells that  mimic  disordered  configurations  constitutes  the  most  common  method  to  compute  the electronic and structural properties of disorder materials with and without defects [32]. In this work,  we  use  the  of  special  quasi-random  supercell  structures  (SQS)  to  model  chemicallydisordered MAX phase supercells as implemented within the Alloy Theoretic Automated Toolkit (ATAT) [33]. SQS are periodic unit cells in which atomic occupancy is engineered to match as close as possible the statistics of truly random alloys at the same chemical composition. The quality of the generated SQS (periodically repeated) supercell is quantified by the correlation function (we choose two and three body interactions to achieve best correlation function).

In this work, an SQS of 200 atoms has been generated to mimic disorder on M-, A- and/or X-site sublattices.  Both  cationic  and  anionic  vacancy  defects  are  considered  in  the  disordered configurations (MM') 2(AA')(X X'). -Antisite pairs are created by interchanging the M (A) or A (M) atoms within neighboring layers. We note here that while smaller cells could have been used, the selected supercell size of O(10 2 ) is necessary in order to accurately compute the energetics of periodic structures without the pathologies arising from their periodic images [34,35].

We use first-principles density functional theory as implemented in Vienna Ab-initio Simulation Package (VASP) [36,37] to study structural properties, electronic properties and vacancy diffusion in disorder MAX phases. The gamma-centered Monkhorst-Pack [38] k-mesh of 1x1x1 and 3x3x3 is used for Brillouin zone integration during geometry-optimization and charge self-consistency, respectively.  Total  energies  and  forces  are  converged  to  10 -5 eV/cell  and  -0.001  eV/Å.  The

Perdew, Burke and Ernzerhof (PBE) generalized gradient approximation is used with a planewave cut-off energy of 533 eV [39].

The vacancy formation energy ( 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 ) in MAX phases is defined as:

$$E _ { f o r m } ^ { V a c } \left ( V _ { M / A / X } \right ) = E _ { \text{ tot} } \left ( V _ { M / A / X } \right ) - E _ { \text{ tot} } ( \text{MAX} ) + \mu _ { M / A / X }$$

where M/A/X, M (=Ti, Ta, Cr, Nb, Zr), A=Al, X=C. E tot (VM/A/X) is the calculated total energy of a cell with defect and Etot(MAX) is the total energy of a pure MAX phase without defects, and µM/A/X is the chemical potential of VM/A/X. For the calculation of defect energies, we use the chemical potential of an isolated gas-phase atom and consider the vacancy defects as charge-neutral due to metallic nature of MAX phases. The chemical potential of isolated atoms is calculated in the same MAX phase SQS supercell by keeping M/A/X atoms at the center of the unit-cell. For this work, we decided not to include the entropic effect because the energy scales (for the formation of defects and migration energies) tend to be in the order of 100 meV and above, with some processes requiring energies well above 4 eV. Therefore, the major conclusions are unlikely to change.

The defect migration (or diffusion barrier) energies are calculated using the cNEB method [40] as implemented in VASP. The cNEB method can be used to identify the minimum energy path of migrating defects from one site to another site within the same unit cell by essentially identifying the set of configurations along the path with vanishingly small atomic forces normal to the path itself. Thus, the cNEB method allows us to understand plausible mechanisms for defect migration as  well  as  the  energy  barriers  associated  to  these  atomic  long-range  displacements.  This information is important not only to understand the kinetics of defect evolution but can also contribute to better understanding and identification of plausible atomic mechanisms leading to defect recovery [41].   Here, the vacancy migration energies in ordered and disordered MAX phase crystals is determined by a transition state search, where the initial and final states are relaxed structures with one vacancy, respectively. We use 6 intermediate images to accurately calculate the  migration  energy  barrier  of  MAX  phases  [42]  from  lowest  energy  position  to  the  other symmetry equivalent position.

## Results and discussion

Structural  properties  and  formation  energies  of  ordered  MAX  phases  with  and  without vacancies :  M2AlC  phase  forms  a  hexagonal structure  (space-group=P63/mmc) that consists of edge-sharing  M6X  (M=cation,  X=anion)  octahedra  interleaved  with  cation  (A)-layer.  The octahedra is similar to the rock salt binary carbides. While our interest is on disordered MAX phases, we first present structural and energy stability of ordered MAX phase and vacancy cases. Table I shows  the  ground-state  structural  properties  and  effect  of  vacancy  on  M2AlC,  where M=Cr,Ta,Nb,Ti,Zr. We compute the equilibrium lattice configurations of each case and find good agreement with experimentally known MAX phases [43-46], ensuring the reliability and accuracy

of our calculations. We allow the vacancy supercell to relax homogeneously to the ground state energies before evaluating the complete set of the elastic constants. In Table. S1 ,  we  list  the elastic constants (Cij), bulk modulus (K), shear modulus (G), Young's modulus I, Poisson ratio (  ) and Pugh ratio k= G/K as obtained by Voigt-Reuss-Hill approximation [47]. The lower Pugh (k) ratio  for  some  MAX  phases  indicates  relative  brittleness  compared  to  other  chemistries. Comparing the bulk moduli of Zr2AlC with Nb2AlC from Table S1 shows that Nb2AlC is elastically much stiffer than usually known allotropes. This is because Nb2AlC has larger elastic moduli all across: a bulk modulus of 170 GPa, a shear modulus of 119 GPa and a Young's modulus of 290 GPa compared to Zr2AlC.

Figure 1 . The  formation  enthalpy  of  (a) pure  and  (b) vacancy  (VM,VA,VX) M2AlC  MAX  phases, M=Cr,Ta,Nb,Ti,Zr. (c) Volume change of vacancy MAX phase with respect no-vacancy cases.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-c3ebda27-d638-411d-84dd-20b12086d267.jpg)

In Fig . 1 , we show (a) the formation enthalpy of ordered MAX phases, (b) defect stability ( 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 ), and (c) relative change in cell volume [  V=Vfinal(M2AlC-vacancy) - Vinitial(M2AlC)] between pure and defect-containing structures. The formation energy (Eform) in Fig. 1a suggests that Ti- and Zrbased MAX alloys are more stable --relative to their elemental constituents --than their other transition metal counterparts. The large VM and VC vacancy formation energies in Fig. 1b makes it  harder  to  create  M-&amp;  C-vacancies  for  Ti2AlC  and  Zr2AlC.  On  other  hand,  the  A-site  vacancy formation is less energetically favorable in Nb2AlC and Ta2AlC. The volume of the MAX phase, in Fig. 1c , is sensitive to the vacancy with positive  V. The Cr2AlC with VC is the only exception in Fig. 1c that shows either very small or no change in volume (attributed to magnetic character of Cr). If we make a comparison, the stable phases in Fig. 1b have smaller shear and Young's moduli (see Table. S1) . The Poisson's ratio of most of the ordered and vacancy MAX phases lie within the usual MAX phase range, i.e., from 0.200 to 0.260. The Pugh ratio (k) for crystalline Ti2AlC (1.29) is smaller than other MAX phases, while some vacancy-containing alloys in Table S1 are fairly large. This indicates that those phases are relatively close to the ductile regime --here we note that  the  ability  of  MAX  phases  to  undergo  significant  deformation  is  the  result  of  complex mechanisms that require analysis well beyond the relatively simplistic metrics of bulk and shear modulus and that remain outside the scope of this work.

A major focus of our work is to look at the effect of chemical alloying on the formation energy, vacancy stability and vacancy migration of disordered MAX phases. As mentioned above, the fact that MAX phases can be synthesized with solid solutions on the M, A, or X sites greatly increases their chemical versatility and significantly expands the MAX phase composition space. In order to emulate disordered solid solutions, we designed an SQS supercell with 50:50 disorder on Msite,  and  M-&amp;A-site  without vacancy to test the  Eform compared to their ordered MAX phase counterparts. The calculated Eform of (Zr-M)2AlC and (Zr-M)2(AA')C are shown in Fig. 2 where M=Cr, Nb, Ti and AA' =Sn -Al, Pb-Bi. The horizonal lines in Fig. 2 represents Eform of the order phase. For example, alloying M-site in Nb2AlC with Zr stabilizes the (ZrNb)2AlC by -0.2 eV/atom, however, if we think of alloying the M-site in Zr2AlC with Nb this slightly reduces the stability by + 0.05 eV/atom. This indicates the tendency of C-Nb-Al to order while C-ZrAl tends to cluster compared to their ordered counterpart. The overall picture of negative formation energy shows the  thermodynamic  stability  of  each  of  disordered  MAX  phase,  relative  to  their  elemental constituents,  considered  in  this  study.  Our  simulations  also  suggest  that  the  compositional dependence of the energy of mixing is complex and non-ideal [48]. The minimum in the energy of mixing obtained for ordered ( Fig. 1 ) and disordered ( Fig. 2 ) MAX phases is consistent with the existence of a solubility limit.

Figure 2 . Formation enthalpy (Eform) of MAX phase for M-site disorder and horizonal dashed (black) lines represent order MAX-phase Eform. Right-panel shows (top) ordered, M2AlC, and (bottom) disordered (MM') 2 AlC supercell. Note that comparison with formation energies of pure MAX phases is just qualitative as we have not computed the convex hull against all possible phases in these systems.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-91eeca71-33a0-4b2d-aded-32dc76fbf201.jpg)

Vacancy defects: In this section, we discuss the effect of vacancies on volume and 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 of pure and disordered max phases. We investigate three types of vacancies, i.e., VM, VA and VX. Fig.

3 (a,c,e) shows the trends of positive volume change with respect to no-vacancy disorder max phase. Point-defecst can lead to overall change in volume with respect to ideal cell, e.g., the defects in nuclear graphite induced by irradiation increases the clat and decreases (alat &amp; blat) in [0001]  plane  in  response  to  interstitial  carbons.  The  introduction  of  vacancies  increases  the volume as shown in Fig. 1c &amp; 3 . The present study shows that M-elements from group-VI of the periodic table, i.e. Cr or Nb based MAX phase are the only exceptions in that they exhibit a drop in relative volume with respect to other cases. Two site-disordered (ZrNb)2(SnAl)C MAX phase shows negative  V. The change in volume is very small but a significant drop in 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 is found compared to one-site disorder. The calculated 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 in Fig. 3(b,d,f ) shows that VM in (ZrTi)2AlC are harder to be formed compared to (ZrNb)2AlC and (ZrCr)2AlC, i.e., 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 [(ZrTi)2AlC] &gt; 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 [(ZrNb)2  AlC]  &gt; 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 [(ZrCr)2  AlC].  However,  the 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 is lower  for  two  site-disordered (ZrM)2[ AA' ]C, where M=Cr,Nb,Ti and AA' = Sn-Al, Pb-Bi, i.e., vacancy creation is easier compared to  the  single-site  disordered  case.  The  (ZrCr)2(SnAl)C  is  the  only  exception,  where  two  site disorder increases VM and VA 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 compared to (ZrNb)2(SnAl)C and (ZrTi)2(SnAl)C. On the other hand, [(ZrCr)2(SnAl)C] and [(ZrNb)2(SnAl)C] have competing 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 as shown Fig. 3d .  The 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 trend  in Fig. 3f for [(ZrM)2(PbBi)C]  is  significantly  reduced  compared  to Fig. 3b : 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 [(ZrTi)2(PbBi)C]  &gt; 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 [(ZrNb)2(PbBi)C]  &gt; 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 [(ZrCr)2(PbBi)C],  i.e.,  it  is  easier  to  create  V A vacancies compared VX in (ZrCr)2AlC.

Figure 3 .  The  effect  of  vacancy  on volume and vacancy formation energy of disorder MAX phase. We consider M-and A-site disorder along with three possible vacancy scenarios: VM, VA and VX. We consider (Zr-M)2AlC and (Zr-M)2(AA')C cases where M=Cr, Nb, &amp; Ti, and A A'=Sn --Al &amp; Pb-Bi, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-64519279-03b4-42fc-bf07-c71f355b7e1d.jpg)

For M=Cr, Nb and Ti in (Zr-M)2AlC, the 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 approximately follows the trend [E(VM) &gt; E(VA) &gt; E(VX)],  [E(VM)  &gt;  E(VX)  &gt;  E(VA)]  and  E(VM)  &gt;  E(VX)  &gt;  E(VA)  for  M=Ti,  respectively.  Notably, 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 depends strongly on choice of atomic µM/A/X. This is the reason, 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 for VA, and/or VC is lower in (Zr-M)2AlC than (Zr-M)2(SnAl)C or [(Zr-M)2(PbBi)C as shown Fig. 3d &amp; 3f . The VZr is the least  affected  while  VC  is  the  most  affected  defect  in  disordered  MAX  phases  due  to  the respective atomic µZr/C. Therefore, irrespective of the choice of µ in ordered MAX phase, it is difficult to create VZr as shown in Fig. 1b . Depending on choice of disorder, VZr becomes relatively easier  to  form,  while  it  is  harder  to  create  VA  compared  to  VX  in  Zr2AlC.  Consequently,  the synthesis  condition,  i.e.,  choice  of  µ,  decides  the  relative  stability  of    VAl  or  VC  in  Zr2AlC. Contrastingly, choice of µ on other MAX phases has a negligible effect on the 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 , where VC remain the most stable defect configurations.

Antisite defect : The 𝐸 𝑓𝑜𝑟𝑚 𝑎𝑛𝑡𝑖𝑠𝑖𝑡𝑒 indicates the difficulty of antisite defect formation and irradiation incited recovery mechanism, which depends on the vacancy types and the target sublattices. The effect of antisite defects on volume (  V), energy stability and charge density difference of pure and defected MAX phases are shown in Fig.  4a &amp; 4b for (Zr-M)2AlC. The 𝐸 𝑓𝑜𝑟𝑚 𝑎𝑛𝑡𝑖𝑠𝑖𝑡𝑒 do not have significant dependence on the atomic configuration as the atomic species are only interchanged, i.e., no atoms are moved in-or-out of the chemical reservoir. We show in Fig. 4a that M-Al antisite pair does not affect  V, but the exchange between M and X increases the volume -with Ti-C as only exception that shows negative  V. Antisite defects leads to positive  V for (ZrCr)2AlC while negative  V  for  (ZrCr)2AlC  and  (ZrTi)2AlC.  Based  on 𝐸 𝑓𝑜𝑟𝑚 𝐴𝑛𝑡𝑖𝑠𝑖𝑡𝑒 in Fig.  4b ,  it  is  harder  to  create antisite pairs at Zr (Zr-C or Zr-Al). Comparison of 𝐸 𝑓𝑜𝑟𝑚 𝐴𝑛𝑡𝑖𝑠𝑖𝑡𝑒 ( Fig. 4b ) with 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 ( Fig. 3b ) suggests that antisite defects for M=Cr are favorable as they have lower formation energies.

We can understand this  in  terms  of  empirical  parameters,  such  as,  electronegativity (  )  and atomi-radii  (R).  Here,  the  electronegativity  difference  of  two  atoms  can  tell  us  about  the displacement  of  charge  when  two  chemical  species  are  mixed  together.  For  example,  in (ZrCr)2AlC,  Cr-Al is less than  Zr-Al or  Zr-Al for  Zr = 1.33,  Cr = 1.66,  Al = 1.61, and  C = 2.5 and  R of Cr-Al is less than Zr-C or Al-C for RZr = 0.86 Å , RCr = 0.76 Å , RAl = 0.53 Å , and RC = 0.29 Å , respectively. The large difference in  R (62% vs 43%) and  (21% vs 3%)) while comparing Zr-Al with Cr-Al, we see  R  and  in  (ZrCr) 2AlC explains large formation energy of the Zr-Al pair. However,  and  R suggest that Cr-Al is the preferred pair of antisite defects where electronic effects play a pivotal role. Clearly, the C-layer is not energetically favorable for Zr/Cr or Al for antisite  defect  formation,  therefore,  it  creates  a  vacancy  by  migrating  into  the  next  available cation M- (or A) layer. This shows that interchanging M- and A-cation position with X-anion costs more energy than anion-anion interchange.

Figure  4 .  Effect  of  antisite  defects  on  (a)  cell  volume  and  (b)  formation  energies  of  (ZrM) AlC,  where 2 M=Cr,Nb,TI. (c-f) Charge density difference between (ZrM)2AlC antisite defect and disorder cells.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-1caf91f6-bb85-4fcd-bb80-a8fe58387899.jpg)

The charge density difference between (Zr-M)2AlC antisite defect and disordered cells in Fig. 4e shows that the in-plane Zrd /Crd interactions shows large charge transfer for M= Cr. The Zr-Cr/Zr-Zr/CrCr bonds in (ZrCr)2AlC are populated with in-plane Zr-Cr/Zr-Zr/Cr-Cr bonds, respectively, whereas the in-plane Zr-Cr/Zr-Zr/Cr-Cr bonds are emptied in favor of the Zr-Cr/Zr-Zr/Cr-Cr basal bonds. This can help us to understand the reason behind the increased stability of the Zr/Cr-C layers in the (ZrCr)2AlC system. To exemplify, we plot Ti2AlC and Ta2AlC band structure in Fig. S7 that shows metallic behavior in the basal plane with multiple bands crossing at the Fermi-level (E Fermi ) while no bands cross the E Fermi in c-direction along (  -A).

Electronic properties : In Fig. 5 , we show the total density of state of disordered (Zr-M)2AlC MAX phases for M=Cr, Nb and Ti. The transition elements Zr, Cr, Nb and Ti have partially filled bonding d-states (also see Fig. S8 ), which means that the electrons at E Fermi are mainly from M elements and form conducting bands. If we look at the partial DOS of the Zr2AlC and Cr2AlC in Fig. S8, the Zrd and  Crd states  in  (-5.5  eV  to  -2.5  eV) and  (-8.5  eV  to  -5  eV)  energy  range,  respectively, overlap significantly with the Cp states. Addition of Cr to Zr2AlC at the Zr site add more electrons that increase the DOS near the E Fermi ( Fig. 5 ), i.e., Cr addition increases the localized d-states. This also signifies that Zrd /Crd &amp; Cp in the alloy form strong covalent bonds, while, the Alp with Zrd and Crd states in the energy range -3eV to 0 eV and -5 eV to -2 eV show a weaker hybridization, respectively. This is suggesting metallic bonding between Alp and Md states.

Figure 5 . The total DOS of the pure (top), antisite (middle) disorder and monovacancy (bottom) defects in (MM')2AlC (M=Zr; M'=Cr,Nb,Ti).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-31521e15-5365-4328-b718-b7aad89c4688.jpg)

We note on comparing TDOS of (Zr-M)2AlC in Fig. 5 with antisite (middle) and vacancy (bottom) defects that the hybridization energy peaks near the E Fermi for vacancy cases changes greatly. This also depends on type of vacancy, i.e., VM, VA, or VX, which alters the electronic density and the strength of hybridization near E Fermi because VM and VC reduces the number of M-C bonds. On the contrary, for the A vacancy, the neighboring atoms readjusts itself to maintain the energy stability by charge transfer. Comparing total DOS in the presence of the vacancy (bottom-panel) and no vacancies (top-panel), we could clearly see the hybridization peak in the same energy range near the E Fermi for vacancy cases. We attribute this to the reduced numbers of M-A bonds. Notably, no such changes are observed for the case of antisite defects as system s electronic ' density (only atomic positions are interchanged) remains conserved

To  further  investigate  this,  we  plot  the  charge  density  difference  between  vacancy  and  novacancy cases in the inset Fig. 5 . The DOS above the E Fermi is the result of weaker Zrd /Crd and Alp bonding with a little contribution form the Cp states (also see supplementary material). On the other hand, the top of the valance band is formed by Zrd /Crd , Cp and Alp hybridized states. For the ordered MAX phase, in Fig. S8 , the states at the EFermi are mainly from the Zrd /Crd , Cp and Alp with a small contribution of Cp .

Figure 6 . Charge difference of M-site disorder (ZrCr)2AlC MAX phase with (a-d) antisite disorder and (e) VM, VA and Vx vacancies. Clearly, antisite and vacancy defect changes the local chemical activity by stronger charge transfer.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-7d0a3c8f-5359-46cb-a02d-e520e8bea1c9.jpg)

Charge density difference of antisite and vacancy defects with disorder (ZrCr)2AlC :  The MAX phases  form  M-X,  M-A,  and  A-A  bonds,  where  M-X  is  the  covalent  and  A-A  is  the  metallic respectively. This is the reason M-X bonds are stronger compared to A-A. If we compare two ordered MAX phases, e.g., Zr2AlC and Cr2AlC, the Zr-based compounds are weaker than the Crbased compounds due to significant difference in  between Zr and Cr (  Zr =  1.33,  Cr = 1.66). Large  of Cr attracts more charges and makes the bonds in Cr2AlC stronger. This also suggests that  Zr2AlC  will  have  more  agility  towards  re-organizing  defects  compared  to  Cr-based  MAX phases. The local bonds near the vacancy sites are disturbed due to vacancy formation as shown in Fig.  6 by  charge  density  difference  between  pure  and  defect  (ZrCr)2AlC.  For  example,  the cationic defects lead to much stronger bonding between M(=Cr/Zr), viz. Zr-Zr/Cr-Cr and Al-C.

Therefore, we believe that alloying and defect (antisite/vacancy) can be used as an effective approach in modifying materials properties. The smaller atomic radius of Cr (0.76 Å ) compared to Zr (0.86 Å ) increases the C-Cr interaction compared C-Zr, which redistributed the electronicdensity  in  and  around  the  A-layer.  We  believe  that  this  electronic  rearrangement  plays  an important role in the increased stability of (ZCr)2AlC. Electronegativity is another factor important to understand the chemical bonding. Following the Pauling electronegativity scale (  Cr =1.66,  Zr =1.33,  Al = 1.61 and  C = 2.55), we find a significant charge transfer in Fig. 6 from Zr/Cr to C, where  the  Zr/Cr-C  bonding  exhibits  directionality  and  therefore  the  covalent  character.  The driving force behind the displacement of the bonding charge is the greater ability of C to attract electrons as a result of the difference between atomic electronegativities (  =  M - C).

Vacancy migrations energies : Vacancies are the most common defects at equilibrium and they greatly control the kinetics of diffusion-mediated processes. For high-temperature applications, vacancy diffusion is an important aspect to explore, particularly because the exchange between vacancies and atoms in the lattice is often times the limiting step when considering possible chemistry-driven transformations, including the selective oxidation of Al in Al-containing MAX phases,  for  example.  To  date,  very  little  progress  has  been  made  towards  increasing  the understanding  of  vacancy  diffusion  mechanism  in disordered  MAX  phases from  a  theoretical framework --we  note  that  a  companion  paper  from  the  present  group  explicitly  looks  at migration energies of vacancies for a large number of pure 211 MAX phase systems [49].

Figure 7. Schematic of VM, VA and VC vacancy migration in order and disorder MAX phase alloy.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-458839bd-49b7-4cb7-896f-8fbe2773241d.jpg)

We study VM, VA, and VX migration in disordered MAX phases as shown by schematic in Fig. 7 . We  use  cNEB  to  generate  a  succession  of  configurations  (images)  along  the  initial  band connecting the initial-to-final state and calculate vacancy formation energies. cNEB relaxes each state on the band and provides information for deforming that band towards a lower energy embedding in the potential energy landscape. The band location in the energy landscape after sufficient number of iterations corresponds to the minimum path connecting the reactants and products.

In Fig.  8(a-c) ,  we  plot  the  VM,  VAl,  and  VC  migration  energies  for  five  ordered  MAX  phase compounds: M2AlC, M=Ti, Ta, Cr, Nb, Zr. The Zr2AlC shows lowest energy barrier for VM, VAl, and VC in comparison to all other case with samaller VAl diffusion barrier height. The overall order of vacancy migration is VM &gt; VAl &gt; VC. The Zr2AlC is the only exception where VM has lower migration energy than VC. Comparing 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 in Fig. 1b with barrier energy in Fig. 8a-c of ordered MAX phase, we can make the inference that stable VAl is easier to move within the basal plane than other vacancies.  In  spite  of  competing  migration  energies  of  VM  and  VC  in  Zr2AlC, 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 is  in  sharp

contrast in term of energy stability in Fig. 1b . Higher 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 for Zr2AlC suggest that vacancies are hard to create, while (relatively) low diffusion barriers suggest facile vacancy migration.

Figure 8 . Calculated energy barrier of VM, VA and VC vacancy migration energies in (a-c) ordered and (d-f) disorder MAX phases using cNEB method. We consider Msite disorder (MM') AlC, M, M'=Ti, Zr, Nb, Ta, 2 Cr. The vacancy at A-sublattice, VA, shows lower vacancy diffusion barrier compared to VM and VC. The plot describes the energy variation along the minimum path. The energy threshold (y-axis), i.e., the activation energy, must be overcome for vacancy migration.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-6d89a7da-117a-499e-a109-8892b600854e.jpg)

The vacancy migration barrier for disordered (ZrM)2AlC in Fig. 8(d-f) follows the trend of 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 in Fig. 3 , where A-site migration is easier in (ZrNb)2AlC compared to (ZrCr)2AlC and (ZrTi)2AlC, i.e., Emig [(ZrCr)2AlC] &gt; Emig [(ZrTi)2AlC] &gt; Emig [(ZrNb)2AlC]. The barrier energies in Fig. 8(d-f) for  the disordered alloys seem to depend greatly on the local environment. This is further established by zero end-point energies in ordered MAX phases, whereas the non-zero end-point energies in disorder  MAX  phases  are  configuration  dependent.  Creating  two  vacancies  at  symmetrically related Wycoff-positions changes the neighbor distributions, which affects the total energy of the vacancy supercell. The vacancy formation energies and vacancy concentration can directly be  connected.  If  we  compare  vacancy  formation  energies  from Fig.  2 ,  e.g.,  the  lower  and comparable formation energies of (ZrNb)2AlC and (ZrTi)2AlC will allow larger vacancy concentration compared to (ZrCr)2AlC. Similarly, for (ZrM)2 (AA')C ,  the  V A  vacancy migration in (ZrM)2AlC within the basal plane is easier and proceeds with an energy barrier of 0.5 to 1.0 eV. This is much lower than the VM or VC vacancy migration barrier. However, the high-migration energy  makes  VM  diffusion  unfavorable  compared  to  VA  and  VX.  Therefore,  the  significant population of the vacancies created in M-layers may remain intact even after a significant long time.

Figure 9 . (a-f) The M-site (VM) vacancy migration barrier for A-site disorder for Zr2(AA')C, A=Al, A'=As, Bi, Pb, S, Sb, Sn. The energy threshold (activation energy) must be overcome for vacancy migration. With Asite disorder, vacancy migration at Zr-site becomes harder.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-76d9c933-f7e8-4af1-9217-fd720de75fd2.jpg)

In Fig. 9 and Fig. 10 , we present barrier energy for A-site disorder [Zr2(AA')C, A=Al, A'= As, Bi, Pb, S,  Sb,  Sn]  and two (M and A)-site disorder ((Zr-M)2(AA')C, where M=Cr,Nb,Ti and A A'=SnAl), -respectively. The high vacancy migration energies in Fig. 9 shows very high barrier of 9-12 eV for A-site disorder with VM. The high barrier makes vacancy diffusion almost impossible, however, once the vacancies are created, the significant population of the vacancies will be intact for a considerable time under the irradiation process. On the other hand, for two-site disorder with VM, VA and VC vacancies in Fig. 10 (also see Fig. S6 ), the diffusion barrier show competing energies for VM and VC for M=Nb and Ti. The A-site vacancies for M=Cr show almost barrierless diffusion. If compared to single M-site disorder migration energies as presented in Fig. 8d-f , adding more disorder, i.e., to both to M- and A-site, further improves the vacancy diffusion by reducing the barrier  height.  The  VM  and  VX  site  vacancies  show  competing 𝐸 𝑓𝑜𝑟𝑚 𝑉𝑎𝑐 in Fig.  3b ,  but  vacancy migration becomes easier for VA with small barrier size of 0.40 eV due to increased chemical disorder.

Figure 10. The VM, VA and VC vacancy migration energies of MAX alloys with M- and A-site disorder using cNEB method: (M-M')2(AA')C, M=Zr, M'=Cr,Nb,Ti, A=Sn and A'=Al.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-91b176ad-0deb-44a5-8d1d-2c442feec118.jpg)

## Conclusion

MAX  phases  are  very  exciting  new  class  of  materials  because  of  their  high-temperature applications and the vacancies are most common defects at high operating conditions. Therefore, we set out the main objective of this work to find out: (1) How the role of vacancies changes in going  from  order  to  disorder  MAX  phases;  and  (2)  Can  we  use  alloying  in  combination  with vacancies as an effective tool to tailor MAX phase properties? We perform first-principles DFT calculations to explore the alloying effect on the structural-stability, energy-stability, electronicstructure,  and  vacancy diffusion of MAX phases. Both as a test-bed and useful application in nuclear engineering, Zr2AlC was an interesting example to be explored for the effect of alloying and  defect  formation  because  the  lower  formation  energy  of  ordered  Zr2AlC  suggests  the possibility of higher vacancy concentration which may in turn be beneficial for the production of protective passive oxide layers.

We considered single-site (M-/A-) and two-site (M-A) disorder and explored the vacancy and antisite defect stability. We found that disordering has a significant effect on energy-stability, electronic-structure and vacancy migration behavior of (ZrM') 2(AA')X (M = Cr, Nb, Ti, AA' = Al, SnAl, X = C), whereas in some cases antisite defects are much easier to form compared to point vacancies.  For  example,  Al  preferentially  goes  to  Cr  compared  to  Zr  in  Cr-doped  Zr2AlC,  i.e., (ZrCr)2AlC. Other focus of this work was to understand the effect of chemical alloying and defects on vacancy diffusion profile of MAX phases. While the diffusion barrier in strongly bonded metal alloys is usually large, our DFT+cNEB calculations show that alloying helps in lowering down both vacancy formation energy and diffusion-barrier compared to ordered MAX phases. The reduced barrier  height  for  Al-diffusion  as  found  in  our  study  through  alloying  is  an  important  finding because Al is known to promote the formation of protective oxide layer (i.e., Al2O3 in aluminumbased MAX phases) at high-operating temperatures. Therefore, the reduced vacancy formation energy  and  barrier  height  will  help  the  formation  of  protective  oxide  layer  in  Al-based  MAX phases. Based on our study, we believe that the chemical alloying route will surely help the MAX phase  community  to  understand  the  defect  formation  and  migration  mechanism  as  well  as provide ways to manipulate the electronic and mechanical behavior of MAX phases.

## Acknowledgement

Valuable discussion with Dr. M. Radovic is acknowledged. The financial support was provided from National Science Foundation through grants no. (DMREF) CMMI-1729350. First-principles calculations  were  carried  out  at  the  Supercomputing  Facility  at  Texas  A&amp;M  University.    DS acknowledges the support of NSF through grant no. NSF-DGE-1545403.

## *prashant40179@gmail.com

## References :

- 1. Barsoum, M.W. MAX Phases: Properties of Machinable Carbides and Nitrides. Wiley VCH, 2013.
- 2. Radovic, M., &amp; Barsoum, M.W. MAX phases: Bridging the gap between metals and ceramics. American Ceramic Society Bulletin 92 , 20-27 (2013).
- 3. Barsoum, M.W., &amp; Radovic, M. Elastic and Mechanical Properties of the MAX phases. Annul. Rev. Mater. Res. 41 , 195-227 (2011).
- 4. Wang, J., &amp; Zhou, Y. Recent progresses in theoretical prediction, preparation and characterization of layered ternary transition metal carbides. Annul. Rev. Mater. Res. 39 10 (129) (2009).
- 5. Lin, Z., Li, M., &amp; Zhou, Y. TEM investigation on layered ternary ceramics. J. Mater. Sci. Technol. 23 , 145 (2007).
- 6. Hoffman, E.N., Vinson, D.W., Sindelar, R.L., Tallman, D.J., Kohse, G. and Barsoum, M.W. MAX phase carbides and nitrides: Properties for future nuclear power plant in-core applications and neutron transmutation analysis. Nuclear Engineering and Design 244 , 17-24 (2012).
- 7. Naguib, M., Mochalin, V.N., Barsoum, M.W., Gogotsi, &amp; Y. MXenes: A New Family of TwoDimensional Materials. Adv. Mater. 26 , 992-1005 (2014).
- 8. Naguib, M., Kurtoglu, M., Presser, V., Lu, J., Niu, J., Heon, M., Hultman, L., Gogotsi, Y., &amp; Barsoum, M.W. Two- Dimensional Nanocrystals Produced by Exfoliation of Ti3AlC2. Adv. Mater. 23 , 4248-4253 (2011).
- 9. Ito, T. &amp; Okazaki, S. Pushing the limits of lithography. Nature 406 , 1027 (2000).
- 10. Sickafus, K. et al. Radiation tolerance of complex oxides. Science 289 , 748 -751 (2000).
- 11. Sickafus, K. E. et al. Radiation-induced amorphization resistance and radiation tolerance in structurally related oxides. Nat. Mater. 6 , 217 -223 (2007).
- 12. Wang, C., Yang, T., Tracy, C.L., Lu, C., Zhang, H., Hu, Y.-J. Wang, L. Qi, L. Gu, L. Huang, Q., Zhang, J., Wang, J., Xue, J., Ewing, R.C., &amp; Wang, Y. Disorder in Mn+1AXn phases at the atomic scale. Nat. Commun. 10 , 622 (2019).
- 13. Eklund, P., Beckers, M., Jansson, U., H ö gberg, H., &amp; Hultman, L. The Mn+1AXn phases: Materials science and thin film processing. Thin Solid Films 518 , 1851 (2010).
- 14. Nordlund, K., Keinonen, J., Ghaly, M. &amp; Averback, R. S. Coherent displacement of atoms during ion irradiation. Nature 398 , 49-51 (1999).
- 15. Zhang, Y. et al. Ionization-induced annealing of pre-existing defects in silicon carbide. Nat. Commun . 6 , 8049 (2015).
- 16. Zhang, Y. et al. Influence of chemical disorder on energy dissipation and defect evolution in concentrated solid solution alloys. Nat. Commun . 6 , 8736 (2015).
- 17. Was, G. S. &amp; Averback, R. S. in Comprehensive Nuclear Materials Vol. 1 (ed. Konings, R. J. M.) 293-332 (Elsevier Ltd., 2012).
- 18. Jager, W., &amp; Wilkens, M. Phys. Formation of vacancy-type dislocation loops in tungsten bombarded by 60 keV Au ions. Status Solidi (a) 32 , 89 (1975).
- 19. Liao, T., Wang, J.Y., &amp; Zhou, YC. First-principles investigation of intrinsic defects and (N, O) impurity atom stimulated Al vacancy in Ti2AlC. Appl. Phys. Lett . 93 , 261911 (2008).

- 20. Zhao, S.J., Xue, J.M., Wang, Y.G., &amp; Huang, Q. Ab initio study of irradiation tolerance for different Mn+1AXn phases: Ti3SiC2 and Ti3AlC2. J. Appl. Phys . 115 , 023503 (2014).
- 21. Marion, L.F., &amp; Monnet, I. Saturation of irradiation damage in (Ti,Zr)3(Si,Al)C2 compounds. J. Nucl. Mater . 433 , 534-537 (2013).
- 22. Xu, Y.G., Ou, X.D., &amp; Rong, X.M. Vacancy trapping behaviors of hydrogen in Ti3SiC2: A firstprinciples study. Materials Letters 116 , 322-327 (2014).
- 23. Ou, X.D., Wang, Y.X., Shi, l.Q., Ding, W., Wang, M., &amp; Zhu, Y.S. Effect of hydrogen-doping on bonding properties of Ti3SiC2. Physica B 406 , 4460-4465 (2011).
- 24. Xiao, J.R., Wang, C.X., Yang, T.F., Kong, S.Y., Xue, J.M., &amp; Wang, Y.G. Theoretical investigation on helium incorporation in Ti3AlC2. Nucl. Instr. Meth. Phys. Res. B 304 , 27-31 (2013).
- 25. Middleburgh, S.C., Lumpkin, G.R., &amp; Riley, D. Accommodation, Accumulation, and Migration of Defects in Ti3SiC2 and Ti3AlC2 MAX Phases. J. Am. Ceram. Soc. 96 , 3196-3201 (2013).
- 26. Nordlund, K., Keinonen, J., Ghaly, M. &amp; Averback, R. Coherent displacement of atoms during ion irradiation. Nature 398 , 49 (1999).
- 27. Wirth, B. D. How does radiation damage materials? Science 318 , 923 -924 (2007).
- 28. Balandin, A. A. Thermal properties of graphene and nanostructured carbon materials. Nat. Mater. 10 , 569 -581 (2011).
- 29. Noriega, R. et al. A general relationship between disorder, aggregation and charge transport in conjugated polymers. Nat. Mater. 12 , 1038 (2013).
- 30. Meneghini, C. et al. Nature of 'disorder' in the ordered double perovskite Sr 2 FeMoO6. Phys. Rev. Lett. 103 , 046403 (2009).
- 31. Tallman, D.J., Anasori, B., &amp; Barsoum, M.W. A Critical Review of the Oxidation of Ti2AlC, Ti3AlC2 and Cr2AlC in Air. Mater. Res. Lett . 1 , 115-125 (2013).
- 32. Zunger, A., Wei, S.H., Ferreira, L.G., &amp; Bernard J.E. Special quasirandom structures. Phys. Rev. Lett . 65 , 353-356 (1990).
- 33. van de Walle, A., Tiwary, P., de Jong, M.M., Olmsted, D.L., Asta, M., Dick, A., Shin, D., Wang, Y., Chen, L.-Q., &amp; Liu, Z.-K. Efficient stochastic generation of Special Quasirandom Structures. Calphad 42 , 13-18 (2013).
- 34. Makov, G., &amp; Payne, M.C. Periodic boundary-conditions in ab-initio calculations. Phys. Rev. B 51 , 4014-4022 (1995).
- 35. Lany, S., &amp; Zunger, A. Assessment of correction methods for the band-gap problem and for finite-size effects in supercell defect calculations: case studies for ZnO and GaAs. Phys. Rev. B 78 , 235104 (2008) .
- 36. Kresse, G. &amp; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 47 , 558-561 (1993).
- 37. Kresse, G., &amp; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 59 , 1758-1775 (1999).
- 38. Monkhorst, H.J., &amp; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 13, 51885192 (1976).
- 39. Perdew, J.P., Burke, K., &amp; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett . 77 3865-3868 (1996).
- 40. Henkelman, G., Uberuaga, B.P., &amp; Johnsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys 113, 9901-9904 (2000).
- 41. Middleburgh, S.C., &amp; Grimes, R.W. Defects and Transport Processes in Beryllium. Acta Mater. 59 , 7095-7103 (2011).

- 42. Tahini, H., Chroneos, A., &amp; Grimes, R.W. Diffusion of E Centers in Germanium Predicted Using GGA+U Approach. Applied Physics 99 , 072112 (1-3) (2011).
- 43. Kun, L., Yuan, Q., &amp; Ji-Zheng, D. First-principles investigation of the vacancy effect on the electronic properties in M2AlC (M = V and Nb). AIP Advances 4 , 107137 (1-7) (2014).
- 44. Wang, J.Y., Zhou, Y.C., Liao, T., Zhang, J., &amp; Lin, Z.J. A first-principles investigation of the phase stability of Ti AlC with Al vacancies. 2 Scripta Mater . 58 , 227-230 (2008).
- 45. Liao, T., Wang, J.Y., &amp; Zhou, Y.C. A first-principles investigation of the phase stability of Ti2AlC with Al vacancies. Scripta Mater . 59 , 854-857 (2008).
- 46. Meyer, B., &amp; F ä hnle, M. Ab initio calculation of the formation energy and the formation volume of monovacancies in Mo. Phys. Rev. B 56 , 13595-13598 (1997).
- 47. Yao, H., Ouyang, L., &amp; Ching, W.Y. Ab initio Calculation of the elastic constants of ceramic crystals. J. Am. Ceram. Soc. 90, 3194-3204 (2007).
- 48. Talapatra, A., Duong, T., Son, W., Gao, H. Radovic, M., &amp; Arroyave, R. High-throughput combinatorial study of the effect of M site alloying on the solid solution behavior of M2AlC MAX phases. Phys. Rev. B 94 , 104106 (1-15) (2016).
- 49. Talapatra, A. Radovic, M., &amp; Arroyave, R. High-throughput Investigation of Vacancy Formation and Migration Energies in 211 MAX Phase Systems . unpublished .

## Supplemental material

## The Effect of Chemical Disorder on Defect Formation and Migration in Disordered MAX Phases

Prashant Singh,  Daniel Sauceda,  and Raymundo Arroyave a a a a Department of Materials Science &amp; Engineering, Texas A&amp;M University, College Station, TX 77843, USA

Figure S1 . Vacancy dependence change in (a) volume, (b) formation energy in (TiTa)2AlC, Ti2(AlSi)C, and TiAl(CN).  (c)  Energy  barrier  of  V M,  V A  and  VC  vacancy  migration  in  (TiTa)2AlC.  We  calculate  vacancy migration minimum energy path at the equilibrium volume in a 200-lattice supercell.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-746524ae-99ff-415c-9b05-b87fe66d1919.jpg)

Figure S2 . The total DOS of the pure (top), antisite (middle) disorder and monovacancy (bottom) defects in (TiTa) AlC. 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-a4eac5cf-c9f4-4d5d-b6b3-9c05c1e5a0dd.jpg)

Figure S3 .  The partial DOS of the pure (left-panel) and vacancy (right-panel) in disorder (TiTa)2AlC. 2D charge density plot showing charge depletion at the vacancy site.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-bfb4274e-23bc-4f5c-af41-8637624b8e8c.jpg)

Figure  S4 .  The  partial  charge-densities  of  the  pure  (left-panel)  and  vacancy  (right-panel)  in  disorder (TiTa)2AlC.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-f438d740-8197-4f8c-877a-1de2a4c4461f.jpg)

Figure S5 . The volume and average bond-length change of +0.6% and 2.3% for M-site disorder with M-vacancy with respect to (TiTa)2AlC with no vacancies.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-0fb4f7c3-4efb-4a26-b170-06c2fed38465.jpg)

Figure S6 . The VM, VA and VC vacancy migration energies of MAX alloys with M- and A-site disorder using cNEB method: (MM')2 (AA')C, M=Zr, M'=Cr,Nb,Ti, A=Sn and A'=Al. We use 200 atom supercell to calculate minimum energy path for vacancy migration at the equilibrium volume.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-52ae581a-c548-4f91-8211-994c1980012d.jpg)

Electronic structure : The local feature of the TDOS near the Fermi-level (E Fermi ) is indicative of the stability. For example, a valley in DOS at E Fermi signifies a lower energy barrier for electrons (E &lt; 0 eV) to move to the unoccupied states, whereas peak usually implies higher energy barrier (E&gt;0 eV). This criterion works reasonably well for Fig. S7 &amp; S8 . Ti 2AlC and Ta2AlC have a valley at E Fermi ,  suggests increased stability. Out of 9-MAX phase ordered alloys in Fig. S7 ,  only Zr2AlC shows a peak at E Fermi in the TDOS. Our qualitative assumption explains the observation of the early transition metal based MAX phases very well [S 1 ].

Figure S7 . Electronic band-structure and partial density of states of Ti2AlC and Ta2AlC ordered MAX phase alloys.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-cd7df2ac-16fe-4d1e-87c5-f459f89617de.jpg)

Figure S8 . Partial density of states (PDOS) of M2AlC ordered MAX  phase alloys for M=Ti, V,Cr,Zr,Nb,Mo,Hf,Ta,W. Only Zr2AlC shows a peak exactly on Fermi-level, which indicates alloy instability.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00038v3-figure-81282649-5368-423d-a5fe-35a71afa8ee5.jpg)

Table S1 - Structural properties of pure ordered MAX phases M2AlC with no-vacancy and withvacancy (VM, VA, VX).

| System       | a    | b    | c     | V     | C 11   | C 12   | C 44   | K     | G     | E     |     | k    |
|--------------|------|------|-------|-------|--------|--------|--------|-------|-------|-------|------|------|
|              | Å    | Å    | Å     | Å 3   | GPa    | GPa    | GPa    | GPa   | GPa   | GPa   |      |      |
| Cr 2 AlC     | 2.86 | 2.48 | 12.70 | 90.2  | 334.8  | 74.1   | 131.1  | 175.1 | 126.0 | 304.9 | 0.21 | 1.39 |
| Ta 2 AlC     | 3.10 | 2.69 | 13.92 | 115.9 | 332.7  | 128.6  | 159.0  | 194.1 | 121.5 | 301.6 | 0.24 | 1.60 |
| Nb 2 AlC     | 3.14 | 2.72 | 13.91 | 118.7 | 323.9  | 83.4   | 141.7  | 170.1 | 119.6 | 290.7 | 0.22 | 1.42 |
| Ti 2 AlC     | 3.05 | 2.64 | 13.78 | 111.2 | 293.1  | 65.2   | 107.7  | 139.2 | 107.5 | 256.6 | 0.19 | 1.29 |
| Zr 2 AlC     | 3.31 | 2.87 | 14.63 | 138.9 | 253.5  | 72.2   | 90.5   | 130.1 | 87.6  | 214.7 | 0.22 | 1.48 |
| Zr 2 AlC-V M | 6.65 | 5.76 | 14.68 | 561.5 | 214.5  | 70.9   | 51.4   | 112.8 | 63.2  | 159.9 | 0.26 | 1.78 |
| Ti 2 AlC-V M | 6.12 | 5.30 | 13.81 | 448.3 | 257.4  | 73.3   | 75.8   | 125.4 | 85.2  | 208.5 | 0.22 | 1.47 |
| Ta 2 AlC-V M | 6.18 | 5.35 | 13.90 | 460.1 | 328.5  | 111.6  | 133.4  | 184.6 | 114.7 | 285.1 | 0.24 | 1.61 |
| Nb 2 AlC-V M | 6.25 | 5.42 | 13.90 | 470.7 | 293.5  | 88.6   | 116.7  | 160.8 | 102.5 | 253.8 | 0.24 | 1.57 |
| Cr 2 AlC-V M | 5.69 | 4.92 | 12.62 | 353.3 | 322.3  | 64.2   | 122.6  | 163.2 | 122.2 | 293.4 | 0.20 | 1.34 |
| Zr 2 AlC-V X | 6.64 | 5.75 | 14.63 | 559.2 | 228.9  | 72.4   | 82.7   | 121.9 | 79.2  | 195.5 | 0.23 | 1.54 |
| Ti 2 AlC-V X | 6.11 | 5.29 | 13.75 | 445.2 | 267.7  | 68.5   | 101.4  | 134.0 | 99.6  | 239.5 | 0.20 | 1.34 |
| Ta 2 AlC-V X | 6.16 | 5.33 | 13.99 | 459.5 | 318.8  | 138.0  | 141.8  | 197.0 | 111.8 | 282.2 | 0.26 | 1.76 |
| Nb 2 AlC-V X | 6.22 | 5.39 | 14.02 | 469.7 | 283.7  | 120.6  | 124.7  | 173.9 | 99.4  | 250.5 | 0.26 | 1.75 |
| Cr 2 AlC-V X | 5.65 | 4.89 | 12.65 | 349.8 | 359.3  | 86.1   | 138.7  | 188.6 | 133.7 | 324.5 | 0.21 | 1.41 |
| Zr 2 AlC-V A | 6.69 | 5.79 | 14.45 | 559.3 | 243.4  | 62.1   | 58.7   | 118.5 | 74.1  | 184.1 | 0.24 | 1.60 |
| Ti 2 AlC-V A | 6.16 | 5.33 | 13.61 | 446.7 | 277.0  | 62.5   | 70.5   | 128.7 | 88.7  | 216.5 | 0.22 | 1.45 |
| Ta 2 AlC-V A | 6.20 | 5.37 | 13.85 | 460.5 | 324.6  | 129.8  | 126.0  | 186.2 | 106.7 | 268.9 | 0.26 | 1.74 |
| Nb 2 AlC-V A | 6.29 | 5.45 | 13.72 | 471.3 | 301.4  | 90.0   | 105.7  | 161.8 | 99.2  | 247.2 | 0.25 | 1.63 |
| Cr 2 AlC-V A | 5.69 | 4.92 | 12.62 | 353.3 | 332.3  | 67.4   | 121.4  | 166.9 | 123.3 | 297.0 | 0.20 | 1.35 |

Table S2 - Structural properties of disorder (ZrCr)2AlC MAX phases with no-vacancy and VM, VA, and VX vacancies.

| System   | a     | b      | c      | V         | Zr-C       | Cr-C       | Al-Al      |
|----------|-------|--------|--------|-----------|------------|------------|------------|
|          | Å     | Å      | Å      | Å 3 /atom | Avg BL (Å) | Avg BL (Å) | Avg BL (Å) |
| No-vac   | 8.174 | 13.979 | 25.355 | 14.077    | 2.23       | 2.02       | 2.89       |
| V M      | 8.142 | 14.002 | 25.330 | 14.118    | 2.24       | 1.96       | 2.78       |
| V A      | 8.159 | 14.001 | 25.372 | 14.147    | 2.22       | 2.02       | 2.79       |
| V X      | 8.158 | 14.000 | 25.317 | 14.136    | 2.27       | 1.98       | 2.96       |

## References :

S1. Mo, Y. Rulis, P., &amp; Ching, W.Y. Electronic Structure and Optical Conductivities Of 20 Max-Phase Compounds. Phys. Rev. B 86, 165122 (1-10) (2012).
Prashant Singh Daniel Sauceda Raymundo Arroyave
The Effect of Chemical Disorder on Defect Formation and Migration in Disordered MAX Phases
MAX phases have attracted increased attention due to their unique combination
of ceramic and metallic properties. Point-defects are known to play a vital
role in the structural, electronic and transport properties of alloys in
general and this system in particular. As some MAX phases have been shown to be
stable in non-stoichiometric compositions, it is likely that such alloying
effects will affect the behavior of lattice point defects. This problem,
however, remains relatively unexplored. In this work, we investigate the
alloying effect on the structural-stability, energy-stability,
electronic-structure, and diffusion barrier of point defects in MAX phase
alloys within a first-principles density functional theory framework. The
vacancy (V$_{M}$, V$_{A}$, V$_{X}$) and antisite (M-A; M-X) defects are
considered with M and A site disorder in (Zr-M)$_{2}$(AA${'}$)C, where
M=Cr,Nb,Ti and AA${'}$=Al, Al-Sn, Pb-Bi. Our calculations suggest that the
chemical disorder helps lower the V$_{A}$ formation energies compared to
V$_{M}$ and V$_{X}$. The V$_{A}$ diffusion barrier is also significantly
reduced for M-site disorder compared to their ordered counterpart. This is very
important finding because reduced barrier height will ease the Al diffusion at
high-operating temperatures, which will help the formation of passivating oxide
layer (i.e., Al$_{2}$O$_{3}$ in aluminum-based MAX phases) and will slow down
or stop the material degradation. We believe that our study will provide a
fundamental understanding and an approach to tailor the key properties that can
lead to the discovery of new MAX phases.