Role of Nanolaminated Crystal Structure on the Radiation Damage Tolerance of Ti3SiC2: Theoretical Investigation of Native Point Defects
Nanolaminated , 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 . In addition, the migration energy barriers of Si vacancy, Si interstitial, and antisite yield very low values: 0.9, 0.6, and 0.3 eV, respectively. The intercalation of Si atomic plane between nanotwinning structures dominates the formation and migration of intrinsic native point defects in . 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.
Ti3SiC2 belongs to the family of nanolaminated carbides and nitrides, the so-called (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 layers are intercalated by the A atomic plane, and the structural units are alternatively stacked along the direction. The 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 . 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 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 . Furthermore, its Weibull modulus is 29.1, ranking the highest value for monolithic ceramics . 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 . 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 . 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 . 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) . The electron-ion interactions were represented by the projector augmented wave (PAW) method . The electronic exchange correlation energy was treated as the generalized gradient approximation (GGA-PBE) . The plane wave basis set cut off was 450 eV. Calculation of the defect structure employed a supercell, which contains 48 atoms. The special -point sampling integration was used over the Brillouin zone by using the Monkhorst-Pack method with for unit cell and for supercell . The lattice constants and internal freedom of the unit cell were fully optimized until the total energy difference was smaller than 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, ) is calculated from the following equation: where the and are total energies of a defected supercell and a perfect supercell, respectively, represents the change in the number of atoms of species (, , or ) during the process of defect formation, and denotes the chemical potential of species . The chemical potentials of , , and atoms are assumed to be the bulk hcp , cubic , 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 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 (, , VSi, and VC) and 8 kinds of antisites (, , , , 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.
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.
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 or position in the Si atomic plane, on the or position in the Ti2 atomic plane, and on the , , or position between the Si and Ti2 atomic planes (totally seven configurations as seen in Figure 2(a)). These positions are labeled as I--Si, I--Ti2, and I--, 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 configurations and one kind of configuration have relatively low formation energies: 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; 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). atoms are initially put on all of the interstitial positions, and atoms are initially put on the I-CTi1 and I-CTi2 positions.
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 is not stable in Ti3C2 unit, and the possible configuration is I--SiTi2 with the DFE of 3.6 eV. will be stable in both the SiTi2 and the Ti3C2 units. The I--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 is the I--Si with the DFE of 2.1 eV. Note that 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.
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 . 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 in Ti3C2 and SiTi2 units. If we compare the formation energies, the most possible disorders are the 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 . 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 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 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, , 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, , and TiSi along the Si atomic plane are very low, only 0.9, 0.6, and 0.3 eV, respectively. This indicates that VSi, , 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.
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 antisite pair, Si Frenkel pair, and C Frenkel pair. The migration energy barriers of VSi, , 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.
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.
H. B. Zhang, Y. W. Bao, and Y. C. Zhou, “Current status in layered ternary carbide Ti3SiC2, a review,” Journal of Materials Science and Technology, vol. 25, no. 1, pp. 1–38, 2009.View at: Google Scholar
J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1996.View at: Google Scholar