Nanolaminated Ti3SiC2, a representative MAX phase,
shows excellent tolerance to radiation damage. In this paper, first-principles
calculations were used to investigate the mechanism of intrinsic point defects
in order to explain this outstanding property. Formation energies of intrinsic point
defects are calculated and compared; and the results establish a low-energy disorder
mechanism in Ti3SiC2. In addition, the migration energy
barriers of Si vacancy, Si interstitial, and TiSi antisite yield very low values: 0.9, 0.6, and 0.3 eV, respectively.
The intercalation of Si atomic plane between Ti3C2 nanotwinning
structures dominates the formation and migration of intrinsic native point defects
in Ti3SiC2. The present study also highlights a novel method
to improve radiation damage tolerance by developing nanoscale-layered structure which
can serve as a sink or rapid recovery channel for point defects.
1. Introduction
Ti3SiC2 belongs to the family of nanolaminated carbides and nitrides, the so-called Mn+1AXn (MAX) phases, where M is an early transition metal, A is an IIIA or IVA group element, and X is carbon or nitrogen. In the crystal structure of MAX phase, the nanotwinned Mn+1Xn layers are intercalated by the A atomic plane, and the structural units are alternatively stacked along the c direction. The Mn+1Xn layer has the rock-salt-type structure and consists of strong covalent M–X bonds, while the M–A bonds are relatively weaker. Due to the nanolaminated structure, MAX phases exhibit good damage tolerance and machinability, which are the most salient properties compared with typical brittle binary carbides and nitrides [1]. Considering the high modulus, high strength, and thermal shock resistance, MAX phases show a great potential applied as high-temperature structure materials. Ti3SiC2 shares the common properties of MAX phase: its fracture toughness (K1C) is as high as 7 MPa·m1/2; Vicker’s hardness is about 5 GPa; Young’s modulus is within 322–333 GPa; and flexural strength is about 450 MPa [2]. Furthermore, its Weibull modulus is 29.1, ranking the highest value for monolithic ceramics [3]. In addition, Ti3SiC2 attracts extensive attentions due to its good high-temperature oxidation resistance, good thermal shock resistance, and good radiation damage tolerance, which highlights the applications in harsh environment, under high temperature, or under nuclear radiation.
Recently, experimental results showed good radiation damage tolerance of Ti3SiC2 [4–9]. This ternary carbide illustrated structure disorder, instead of amorphization, under heavy ion irradiation up to a damage of ~25 dpa [6]. Ti3SiC2 also demonstrated small swelling of the (0001) atomic planes after irradiation. The swelling is only 2.2% for an average irradiation dose of 4.3 dpa at room temperature [7]. At a temperature of above 300°C, the radiation damage can be easily recovered [8, 9]. It is interesting to compare the radiation damage tolerance of Ti3SiC2 with the binary carbide SiC. Amorphization appears in SiC at low level of damage, usually <1 dpa; and its swelling is about seven times higher than that of Ti3SiC2 [7]. The main origin of good radiation damage tolerance of Ti3SiC2 might be traced back to the nanolaminated crystal structure in which the nanotwinning structure could serve as a sink or recombination channel of point defects.
The behavior of point defect is one of the key factors to control the response of a material to radiation damage. For Ti3SiC2 under radiation, structure disorder instead of amorphization, small swelling, and absence of extended defects, should mainly relate to the formation and recovery of point defects. Therefore, the formation and migration of intrinsic point defects are studied using the first-principles calculations in this work. The results were compared with those of SiC to illustrate the influence of the nanolaminated crystal structure on the point defect behaviors. The aim of this work is twofold: firstly, to explain the unusual radiation damage tolerance by exploring the mechanisms of native point defects in Ti3SiC2; secondly, to elucidate the key role of the nanolaminated crystal structure on the generation and recombination of native point defects in Ti3SiC2. More important, the present study could highlight a novel method to improve radiation damage tolerance by developing nanoscale-layered structure which can serve as a sink or rapid recovery channel for point defects.
2. Computational Method
This calculation was performed by using the VASP code, based on the density functional theory (DFT) [10]. The electron-ion interactions were represented by the projector augmented wave (PAW) method [11]. The electronic exchange correlation energy was treated as the generalized gradient approximation (GGA-PBE) [12]. The plane wave basis set cut off was 450 eV. Calculation of the defect structure employed a 2×2×1 supercell, which contains 48 atoms. The special k-point sampling integration was used over the Brillouin zone by using the Monkhorst-Pack method with 12×12×2 for unit cell and 6×6×2 for supercell [13]. The lattice constants and internal freedom of the unit cell were fully optimized until the total energy difference was smaller than 1×10-6 eV. According to our previous studies on defects in MAX phases, the present calculation method has been proved to be accurate enough to reproduce the defect structure and formation energies [14–16]. To study the migration of defects, the diffusion energy barrier was calculated by searching the transition state linking the defect configurations before and after the migration process. A generalized synchronous transit method, LST/QST, employed in the CASTEP code was used for locating the transition state [17, 18]. It combines the linear (LST) or quadratic synchronous transit (QST) methods with conjugate gradient refinements. The LST/optimization and QST/maximization calculations were firstly performed to search the transition state, and then CG minimization was carried out to refine the saddle point geometry. The cycle was repeated until a stationary saddle point (transition state) was located.
The defect formation energy (DFE, Ef) is calculated from the following equation:
(1)Ef=Edef-Eperf+niμi,
where the Edef and Eperf are total energies of a defected supercell and a perfect supercell, respectively, ni represents the change in the number of atoms of species i (i=Ti, Si, or C) during the process of defect formation, and μi denotes the chemical potential of species i. The chemical potentials of Ti, Si, and C atoms are assumed to be the bulk hcp Ti, cubic Si, and graphite, respectively. The chemical potentials of these bulk materials are obtained from the total energies of first-principles calculation.
3. Results and Discussion
The crystal structure of Ti3SiC2 is shown in Figure 1. There are two sites (Ti1 and Ti2) of Ti atoms: one Si site and one C site. According to the crystal structure, the structures of on-lattice defects in Ti3SiC2 can be easily constructed. There are 4 possible on-lattice vacancies (VTi1, VTi2, VSi, and VC) and 8 kinds of antisites (SiTi1, CTi1, SiTi2, CTi2, TiSi, CSi, TiC, and SiC). Calculation results show that the vacancy formation energy on the Ti2 site is 1.8 eV lower than that on the Ti1 site, indicating that the concentration of Ti2 vacancy is much higher than that of Ti1 vacancy. Therefore, the vacancy on Ti2 site is used for further discussion of vacancy on Ti lattice site.
Crystal structure of Ti3SiC2. Letters on the left side are the stacking sequence of Ti and Si atoms.
Calculated formation energies of on-lattice defects in Ti3SiC2 are listed in Table 1. VTi has a relatively high formation energy ~5.5 eV. The values of VSi and VC are both 2.1 eV, only half of the DFE of VTi in Ti3SiC2. From Table 1, it can be found that the DFE of antisites between Ti and C atoms is the highest and the DFE of antisites between Ti and Si atoms is the lowest. The DFE of antisites between C and Si atoms locates between the two end values.
Defect formation energies (DFEs) (in eV) of on-lattice defects in Ti3SiC2.
VTi
VSi
VC
SiTi
TiSi
CTi
TiC
CSi
SiC
DFE (eV)
5.5
2.1
2.1
3.2
1.8
5.9
8.0
3.2
3.4
The configurations of interstitial defects are complicated compared with the on-lattice defects. For simplification, the unit cell of Ti3SiC2 is divided into two regions, namely, the SiTi2 and Ti3C2 units. In the SiTi2 unit, interstitial atoms are initially placed on the b or c position in the Si atomic plane, on the a or c position in the Ti2 atomic plane, and on the a, b, or c position between the Si and Ti2 atomic planes (totally seven configurations as seen in Figure 2(a)). These positions are labeled as I-b/c-Si, I-a/c-Ti2, and I-a/b/c-SiTi2, respectively. The Ti3C2 unit in Ti3SiC2 is similar to the rock-salt structure. Since all of the octahedral interstitials are occupied by C atoms in the Ti3C2 unit, it is very hard for the Si atom migrating inside and locating at the interstitial positions. Consequently, only the Ti and C interstitials are considered in Ti3C2 unit. The interstitials in rock-salt carbide have been discussed before [19, 20]. Two kinds of Ci configurations and one kind of Tii configuration have relatively low formation energies: Ci atom is inserted into the center of two neighboring C atoms to form a linear C-C-C trimmer or to occupy the site near the centre of C or Ti tetrahedron; Tii atom locates at the centre of C or Ti tetrahedron. As a result, four interstitial configurations are adopted in the Ti3C2 unit: interstitial atoms locate between C and Ti1 atomic planes (labeled as I-CTi1), between C and Ti2 atomic planes (labeled as I-CTi2), inside the C atomic plane (labeled as I-C), and inside the Ti1 atomic plane (labeled as I-Ti1). The configurations are related to the sites on the center of C or Ti tetrahedron as seen in Figures 2(b) and 2(c) and the C-C-C trimmer as seen in Figures 2(d) and 2(e). Ci atoms are initially put on all of the interstitial positions, and Tii atoms are initially put on the I-CTi1 and I-CTi2 positions.
Interstitial configurations in SiTi2 unit (a) and in Ti3C2 unit (b)–(e). Letters above the Si plane in (a) are the stacking positions. Small red balls represent the positions of interstitial atoms.
After the geometry optimization, some configurations are not stable and would change to other configurations. The DFEs of all of the stable interstitial configurations are listed in Table 2. It is found that Tii is not stable in Ti3C2 unit, and the possible configuration is I-c-SiTi2 with the DFE of 3.6 eV. Ci will be stable in both the SiTi2 and the Ti3C2 units. The I-b-Si and I-CTi2 configurations are the most possible defects configurations with the DFEs of 0.9 eV and 3.6 eV in SiTi2 and Ti3C2 units, respectively. The most possible configuration of Sii is the I-c-Si with the DFE of 2.1 eV. Note that Sii initially located in Ti2 atomic plane or between the Si and Ti2 atomic planes converges to the interstitial sites between the C and Ti2 atomic planes after the geometry optimization, while the corresponding DFE has a very high value of 7.0 eV.
Calculated defect formation energies (in eV) of stable interstitials in Ti3SiC2.
Blocks
Sites
C
Ti
Si
Ti3C2 unit
I-C
I-Ti1
4.8
I-CTi1
3.9
I-CTi2
3.6
7.0
SiTi2 unit
I-a-Ti2
8.2
I-c-Ti2
2.8
I-a-SiTi2
3.1
5.3
I-b-SiTi2
5.0
I-c-SiTi2
3.6
I-b-Si
0.9
5.1
4.2
I-c-Si
1.8
4.1
2.1
The main consequence of displacive radiation damage is the structural disorder caused by the accumulation of point defects. The possibility of accommodating structure disorder is the key factor to prevent amorphization. Sickafus and collaborators did experimental and theoretical investigations on the radiation damage tolerance of A2B2O7 pyrochlores [21]. In their work, the calculated formation energies of cation antisite pair have very low values and were used to evaluate the resistance to radiation-induced amorphization. The related method has been proved to be a reliable way to predict radiation damage tolerance of a material. Therefore, the structure disorders in Ti3SiC2 are analyzed and compared with those in SiC to illustrate the mechanisms of accommodating radiation damage.
From the calculated point defect formation energies, the formation energies of Frenkel pairs and the antisite pairs in Ti3SiC2 can be obtained accordingly. The following equations are the possible Frenkel and antisite pairs in Ti3SiC2 and their formation energies.
Two formation energies of C Frenkel pairs correspond to the Ci in Ti3C2 and SiTi2 units. If we compare the formation energies, the most possible disorders are the TiSi+SiTi antisite pair, Si Frenkel pair, and C Frenkel pair for their low formation energies.
The formation energies of point defects in 3C-SiC had been reported by Lucas and Pizzagalli [22]. The lowest formation energies of C and Si Frenkel pairs are 6.73 eV and 13.46 eV, respectively. The formation energy of antisite pair of CSi and SiC is 7.5 eV in SiC. So the mechanism of disorder in SiC should be the C Frenkel pair and CSi+SiC antisite pair. It is noted that the formation energies of SiC are obviously much higher than those of the Ti3SiC2.
Based on previous analysis, it can be found that the intercalated Si atomic plane, as well as the nanolaminated crystal structure, contributes to the mechanisms of disorders in Ti3SiC2 by two aspects: firstly, to provide a possible process of Si Frenkel pair to accommodate disorder with very low formation energy and, secondly, to present a possible mechanism of TiSi+SiTi antisite pair with low formation energy. Therefore, Ti3SiC2 is more resistant to radiation damage than SiC, which had been proved in experiments.
The migrations of point defects dominate the recombination process of defects under radiation damage. Therefore, lower migration energy of a point defect means an easier recovery of Ti3SiC2 from radiation damage. As discussed above, the weakly bonded Si atomic plane provides two kinds of additional disorder mechanisms: Si Frenkel pair and Si-Ti-related antisite pair. As a consequence, the migrations of VSi, Sii, and TiSi have the close relationships to the recovery of radiation damages. Since TiSi locates in the Si atomic plane, it is expected that it can migrate via on-lattice VSi site along the Si atomic plane. The calculated migration energy barriers of VSi, Sii, and TiSi along the Si atomic plane are very low, only 0.9, 0.6, and 0.3 eV, respectively. This indicates that VSi, Sii, and TiSi defects can easily migrate along the region between the neighbouring nanotwined Ti3C2 blocks. These processes provide quick recovery mechanisms for Ti3SiC2 under radiation damage. The results also illustrate the predominant role of nanolaminated crystal structure on the excellent radiation tolerance of Ti3SiC2.
4. Conclusion
In summary, the stable defect structures, defect formation energies, and migration barriers of intrinsic point defects in Ti3SiC2 are studied by using first-principles calculations. The results clearly show that the intercalated Si atomic plane, as well as the nanolaminated crystal structure of Ti3SiC2, not only provides new disorder mechanisms with low formation energies but also presents the quick migration paths of the related point defects. The reported disorder mechanisms with low formation energies are the TiSi+Si Ti antisite pair, Si Frenkel pair, and C Frenkel pair. The migration energy barriers of VSi, Sii, and TiSi along the Si atomic plane are only 0.9, 0.6, and 0.3 eV, respectively. The low formation energies of the defect pairs ensure the high ability of Ti3SiC2 to accommodate structure disorders under displacive radiation damage. At the same time, the low migration energy barriers of point defects in Ti3SiC2 indicate easy recovery of radiation damage. These theoretical results are helpful to explain the good radiation damage tolerance of Ti3SiC2 and can shed a light on the theoretical predication of ceramics with good radiation damage tolerance. This further highlights a novel method to improve the radiation damage tolerance by developing nanoscale-layered structure which can serve as the sink or rapid recovery channel for radiation-induced point defects.
Conflict of Interests
The authors declare no conflict of interests.
Acknowledgments
Haibin Zhang is grateful to the foundation by the Recruitment Program of Global Youth Experts and the Youth Hundred Talents Project of Sichuan Province. This work was also supported by the Natural Sciences Foundation of China under Grant Nos. 50832008, 51372252, and 91326102, the Science and Technology Development Foundation of China Academy of Engineering Physics (Grants nos. 2012B0302035 and 2013A0301012), and Science and Technology Innovation Research Foundation of Institute of Nuclear Physics and Chemistry.
BarsoumM. W.MN+1AXN phases: a new class of solids; thermodynamically stable nanolaminates2000281–42012812-s2.0-003451263610.1016/S0079-6786(00)00006-6ZhangH. B.BaoY. W.ZhouY. C.Current status in layered ternary carbide Ti3SiC2, a review20092511382-s2.0-60949090321BaoY. W.ZhouY. C.ZhangH. B.Investigation on reliability of nanolayer-grained Ti3SiC2 via Weibull statistics20074212447044752-s2.0-3454730884310.1007/s10853-006-0427-7NappéJ. C.GrosseauP.AudubertF.GuilhotB.BeauvyM.BenabdesselamM.MonnetI.Damages induced by heavy ions in titanium silicon carbide: effects of nuclear and electronic interactions at room temperature200938523043072-s2.0-6184916456710.1016/j.jnucmat.2008.12.018Le FlemM.LiuX.DoriotS.CozzikaT.MonnetI.Irradiation damage in Ti3(Si,Al)C2: a TEM investigation2010767667752-s2.0-7844927537610.1111/j.1744-7402.2010.02523.xWhittleK. R.BlackfordM. G.AughtersonR. D.MoriccaS.LumpkinG. R.RileyD. P.ZaluzecN. J.Radiation tolerance of Mn+1AXn phases, Ti3AlC2 and Ti3SiC220105813436243682-s2.0-7995331583210.1016/j.actamat.2010.04.029NappéJ. C.MauriceC.GrosseauP.AudubertF.ThoméL.GuilhotB.BeauvyM.BenabdesselamM.Microstructural changes induced by low energy heavy ion irradiation in titanium silicon carbide2011318150315112-s2.0-7995331243710.1016/j.jeurceramsoc.2011.01.002LiuX. M.Le FlemM.BéchadeJ.OnimusF.CozzikaT.MonnetI.XRD investigation of ion irradiated Ti3Si0.90Al0.10C2201026855065122-s2.0-7674916221410.1016/j.nimb.2009.11.017LiuX. M.Le FlemM.BéchadeJ. L.MonnetI.Nanoindentation investigation of heavy ion irradiated Ti3(Si,Al)C220114011–31491532-s2.0-7875158876810.1016/j.jnucmat.2010.04.015KresseG.FurthmüllerJ.Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set199654161116911186
10.1103/PhysRevB.54.11169KresseG.JoubertJ.From ultrasoft pseudopotentials to the projector augmented-wave method19995931758177510.1103/PhysRevB.59.1758PerdewJ. P.BurkeK.ErnzerhofM.Generalized gradient approximation made simple19967718386538682-s2.0-4243943295PackJ. D.MonkhorstH. J.‘Special points for Brillouin-zone integrations’—a reply1977164174817492-s2.0-424345641210.1103/PhysRevB.16.1748WangJ. Y.ZhouY. C.LiaoT.ZhangJ.LinZ. J.A first-principles investigation of the phase stability of Ti2AlC with Al vacancies200858322723010.1016/j.scriptamat.2007.09.048LiaoT.WangJ. Y.ZhouY. C.Ab initio modeling of the formation and migration of monovacancies in Ti2AlC20085988548572-s2.0-4874910115510.1016/j.scriptamat.2008.06.044LiaoT.WangJ. Y.ZhouY. C.First-principles investigation of intrinsic defects and (N, O) impurity atom stimulated Al vacancy in Ti2AlC20089326326191110.1063/1.3058718SegallM. D.LindanP. J. D.ProbertM. J.PickardC. J.HasnipP. J.ClarkS. J.PayneM. C.First-principles simulation: ideas, illustrations and the CASTEP code20021411271727432-s2.0-003717100510.1088/0953-8984/14/11/301GovindN.PetersenM.FitzgeraldG.King-SmithD.AndzelmJ.A generalized synchronous transit method for transition state location20032822502582-s2.0-014200838210.1016/S0927-0256(03)00111-3TsetserisL.PantelidesS. T.Vacancies, interstitials and their complexes in titanium carbide20085612286428712-s2.0-4464909886110.1016/j.actamat.2008.02.020KimS.SzlufarskaI.MorganD.Ab initio study of point defect structures and energetics in ZrC20101075805352110.1063/1.3309765SickafusK. E.MinerviniL.GrimesR. W.ValdezJ. A.IshimaruM.LiF.McClellanK. J.HartmannT.Radiation tolerance of complex oxides200028954807487512-s2.0-003460448610.1126/science.289.5480.748LucasG.PizzagalliL.Structure and stability of irradiation-induced Frenkel pairs in 3C-SiC using first principles calculations200725511241292-s2.0-3384694793510.1016/j.nimb.2006.11.047