Lingyun
Qian
a,
Honggang
Bao
a,
Rui
Li
*a and
Qing
Peng
*bc
aSchool of Mechanical Engineering, University of Science and Technology Beijing, Beijing, 100083, China. E-mail: lirui@ustb.edu.cn
bPhysics Department, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia. E-mail: qing.peng@kfupm.edu.sa
cK.A.CARE Energy Research & Innovation Center at Dhahran, Dhahran, 31261, Saudi Arabia
First published on 21st December 2021
High irradiation tolerance is a key feature required for designing structural materials for next-generation nuclear reactors, for which high entropy alloys (HEAs) and equiatomic multicomponent single-phase alloys are good candidates. In this paper, the effect of chemical complexity on irradiation resistance is investigated using molecular dynamics simulations. Taking the classical Cantor alloy CoCrFeMnNi as a paradigm, we have studied the number of defects, the size of defect clusters and dislocation of six FCC materials, pure Ni, FeNi, CrFeNi, CoCrFeNi, CoCrCuFeNi and CoCrFeMnNi, under consecutive bombardment. The fewest defects can be observed in two HEAs, and so can the size of defect clusters. Multicomponent alloys also show better performance in terms of irradiation resistance than pure Ni. The reason can be attributed to the complexity of high entropy alloys and multicomponent alloys, which leads to lattice distortion and hysteresis diffusion effects. The longest and shortest dislocation lengths are observed in CoCrCuFeNi and CoCrFeMnNi. A larger number of stacking faults is observed in CoCrCuFeNi. The fewer defects and dislocations of CoCrFeMnNi indicate its promising applications in nuclear structural materials.
HEAs are a type of alloys formed by five or more equal or approximately equal amounts of metallic elements.9–11 The definition of EAMCs is similar, but they usually have a smaller number of types of constituent elements than HEAs. It has been confirmed that many multicomponent alloys and high entropy alloys have excellent mechanical properties, good wear resistance,12,13 structural stability14,15 and corrosion resistance.16,17 Alloys consisting of Co, Cr, Fe, Ni, Mn, Cu elements are one kind. Gludovatz et al.18 showed that due to strain hardening, CoCrNi alloy could achieve excellent damage resistance at low temperatures, which implied a good combination of strength, ductility and fracture toughness. Wu et al.19 indicated that some alloys (e.g. CoCrNi, CoCrMnNi and CoCrFeNi) exhibited yields and ultimate strength that increased strongly with decreasing temperature, while others (e.g. CoNi and Ni) exhibited very weak temperature dependency. Buluc et al.20 measured the hardness and friction properties of non-equivalent AlCrFeMnNi high entropy alloys through experiments. Ye et al.21 found that CoCrFeMnNi had excellent corrosion resistance, which was higher than that of A36 steel substrate and similar to that of 304 stainless steel. Chen Zhen et al.22 reported that compared to Al0.75CrCoFeNi, the removal of the Co element decreased the strength, hardness and ductility. With the addition of Ti, Al0.75CrCoFeNiTi0.25 had higher hardness and strength. The results showed that the addition of different elements influenced the alloy properties, helping to design high entropy alloys.
EAMCs and HEAs in general have better radiation resistance than pure metals or conventional alloys.23,24 This could be attributed to their lattice distortion and sluggish diffusion effects,25–28 which reduce the migration and diffusion of defects. Wang et al.29 studied the irradiation of CoCrCuFeNi HEA under different doses of He+ ion beams. Kumar et al.30 investigated the radiation resistance of CrFeMnNi and FeCrNi. The results showed that the introduction of Mn enhanced the radiation resistance of the material. Do et al.31 indicated that the origin of high resistance to irradiation in the CoCrFeMnNi HEA was lattice distortion, which was caused by the complexity of the alloy. Lu et al.32 found that the enhanced swelling resistance of FeNi, CoFeNi, CoCrFeNi and CoCrFeMnNi could be attributed to a change in tailored interstitial defect cluster motion, which enhanced the recombination of point defects. Li et al.33 indicated that CoCrFeMnNi had stronger resistance to surface bombardment. Xia et al.34 pointed out that AlxCoCrFeNi (x = 0.1, 0.75 and 1.5) had a lower volume swelling rate under high-dose irradiation than other commonly used radiation-resistant materials, such as M316 stainless steel or pure zirconium. Granberg et al.35,36 compared the defects in pure Ni, FeNi, CoFeNi and CoCrNi under high-energy ion irradiation. They indicated that defect accumulation was related to the difference in edge dislocation mobility, and the irradiation damage did not necessarily decrease with an increase in the amount of complexity in the alloys.
Despite these studies, the influence of elements in alloys on irradiation resistance is still elusive, especially in the range from a medium entropy alloy to a high entropy alloy. In this work, we have investigated the irradiation resistance of Ni, FeNi, CrFeNi, CoCrFeNi, CoCrCuFeNi and CoCrFeMnNi, in increasing order of chemical complexity, in a single-phase FCC structure under multiple cascading collisions by molecular dynamics simulation. The results are helpful for the design of irradiation-resistant materials.
![]() | ||
Fig. 1 Illustration of six FCC materials. (a–f) are Ni, FeNi, CrFeNi, CoCrFeNi, CoCrCuFeNi and CoCrFeMnNi, respectively. |
The interatomic interactions of CoCrCuFeNi are described by the EAM potential.37 The short-range interaction of atomic collisions in the cascade process is described by the ZBL potential.38 The interatomic interactions in CoCrFeMnNi are described by the 2NN MEAM potential.39 Meanwhile, two parameters are adjusted in order to describe the short-range interatomic interaction in high-energy atomic collisions in cascades correctly. This change has been confirmed by Do et al.31 to be sufficient, and so the effect on the fitted physical properties could be ignored.
The MD simulations use the LAMMPS package.40 The systems are annealed before conducting bombardment simulation. The ensemble is NPT with a timestep of 0.001 ps. The following are the details of the annealing process. First, the simulation systems relax at room temperature 300 K for 100 ps. Then a temperate gradient of 3.4 K ps−1 is applied to increase the temperature to 2000 K. The temperature of 2000 K is held for 1000 ps for relaxation and equilibrium. Finally, the same temperature gradient is applied to decrease the temperature to 300 K. The NVE ensemble is applied in the simulation of irradiation.
An adaptive time step is applied because the movement of atoms in a step should be controlled to within 0.002 nm. The energy dissipation during the bombardment process is achieved by applying a Berendsen thermostat41 to the boundary atoms within 0.5 nm. The temperature is 300 K. The simulation systems relax for 2000 ps at room temperature to make sure that the systems reach equilibrium before the bombardment process starts. The PKA energy of Ni atoms is 3000 eV,33 which is selected because it can cause a reasonable ratio of defects in the matrix.42 The six models are all bombarded 500 times. After each bombardment, 50000 steps of relaxation are conducted before the next bombardment, which is sufficient to cause the energy of the system to stabilize before the next bombardment starts. The visualization and analysis of the simulation process are done using the OVITO visualization tool.43
After 300 periods of bombardment, the number of defects in CrCoCuFeNi and CoCrFeMnNi HEAs is relatively stable at about 500. The number of defects in CoCrFeNi decreases significantly, but is still more than the defects in the two HEAs. The number of defects in pure Ni and FeNi rise rapidly. The number of defects in CrFeNi increases first, followed by a decrease at the end. At the end of the bombardment, the largest number of defects occurs in pure Ni, followed by FeNi, CrFeNi, CoCrFeNi, CoCrCuFeNi, and CoCrFeMnNi. The annihilation of defects in the two HEAs is relatively slight. However, there are multiple annihilation processes in the other four materials. The different behaviors can be attributed to the more complex elements in HEAs, which lead to lattice distortion and a hysteresis diffusion effect as other researchers have reported. Therefore, the movement of atoms and defects is suppressed, so fewer defects and annihilation processes occur.
Defect clusters of different sizes are generated during the bombardment process, as shown as Fig. 2b. Large clusters in which the number of defects is more than 200 occur in FeNi and CrFeNi. Small clusters dominate in the other materials, especially for the CoCrFeMnNi HEA, in which the defect clusters are all smaller than 100. In addition, the growth of defect clusters is different in the six materials, as shown as in Fig. S1 in the ESI.† Clusters where the number of defects is more than 100 occur in Ni, FeNi, CrFeNi and CoCrFeNi after 100 periods of bombardment. However, the growth of defect clusters in the two HEAs is relatively slow. The results show that the defects in the HEAs disperse and will be helpful for better irradiation resistance performance.
To further validate the above speculation, the dislocation and atomic distribution of stacking faults during the whole bombardment in CoCrCuFeNi and CoCrFeMnNi are examined, as shown as in Fig. 4 and 5, respectively. The results of the other four materials can be found in Fig. S3 in the ESI.† As shown as in Fig. 4, during the irradiation process, multiple cascade collisions cause stacking faults in CoCrCuFeNi. A large number of dislocations are formed between the region including stacking faults and other regions, in which 1/6〈112〉 Shockley dislocations dominate. Most stacking fault atoms are in dislocation loops. However, the dislocations in CoCrFeMnNi are much fewer, and tend to self-heal during the bombardment process, as shown as in Fig. 5.
![]() | ||
Fig. 4 Dislocation and atomic distribution of stacking faults in CoCrCuFeNi HEA. The dots and lines have the same meaning as in Fig. 3. (a)–(e) are the cases after 100, 200, 300, 400, and 500 periods of bombardment, respectively. |
![]() | ||
Fig. 5 Dislocation and atomic distribution of stacking faults in CoCrFeMnNi HEA. The dots and lines have the same meaning as in Fig. 3. (a)–(e) are also the cases after 100, 200, 300, 400, and 500 periods of bombardment, respectively. |
The difference in components between the two HEAs are the elements Cu and Mn. The stacking fault energies of CoCrCuFeNi and CoCrFeMnNi were calculated, as shown in Fig. 6, and are 67.52 mJ m−2 and −31.29 mJ m−2, respectively. Zhao et al.45 investigated the stacking fault energies of a series of multiple alloys using first principles theory and indicated that the presentation of negative stacking fault energies in multiple alloys was important for vibrational entropy in stabilizing their FCC structures. In our results, the negative stacking fault in CoCrFeMnNi suggests that the FCC structure of CoCrFeMnNi is more stable; therefore, few dislocations or stacking faults occur.
The longest and shortest dislocation lengths are observed in CoCrCuFeNi and CoCrFeMnNi, respectively, which show different trends in defects. A longer length of dislocation in an alloy leads to a large amount of stacking faults in CoCrCuFeNi. The lower number of defects and dislocations indicates that the CoCrFeMnNi high entropy alloy has excellent irradiation resistance. Our atomistic insights into the effect of chemical complexity on irradiation resistance might be helpful in the design of structural materials for next-generation nuclear reactors.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ma01184g |
This journal is © The Royal Society of Chemistry 2022 |