A Molecular Dynamics Study of Tungsten’s Interstitial Dislocation Loops Formation Induced by Irradiation under Local Strain

: A molecular dynamics simulation was used to investigate the effect of applied strain on the formation of primary defects and the probability of interstitial dislocation loops (IDLs) formation of tungsten (W) during a collision cascade event. The research investigated primary knock-on atom energies of 1, 6, 10, and 14 keV, applied on a deformed W structure (form − 1.4~1.6%). The peak and surviving number of Frenkel pairs (FPs) increased with increasing tension; however, these increases were more pronounced under higher strain due to the formation of IDLs. For 10 self-interstitial atoms (SIA) lengths, the strain effect reduces the clustering energy of the IDLs by about 7 eV. In general, the current ﬁndings suggest that strain effects should be carefully considered in radiation-damaged environments, particularly in low-temperature, high-radiation-energy environments. The compressed condition may advantage materials used in high-radiation-damage devices and power systems.


Introduction
Due to its excellent mechanical properties, high melting point, high thermal conductivity, and resistance to sputtering and erosion, tungsten (W) has garnered interest for use in plasma-facing components [1][2][3]. A material suitable for use in fusion reactors is the material that functions well in high-energy ion impact, high temperature, and strain environments [4][5][6][7][8]. The interaction of irradiation sources with the construction materials of fusion reactors affects their structural and mechanical properties [9,10]. Energetic particles such as protons, neutrons, electrons, and ions are more likely to produce recoils, which are high-energy primary knock-on atoms (PKA) that initiate the collision cascade events. These collisions generate localized point or line defects, including dislocations, that promote the formation of defect clusters, voids, and stacking faults, resulting in material integrity loss [10,11].
Dislocations are one of the most critical defect structures to examine when attempting to comprehend a material's mechanical properties. Interstitial dislocation loops (IDLs) have been observed by the deformation and irradiation interactions with the structure.
IDLs can be considered mechanically as barriers that interact with dislocations based on the Orowan mechanism and result in the material's hardening. IDLs are unwanted defects that contribute to irradiation creep, swelling, and other undesirable effects that compromise the long-term reliability and safety of a system's structure. In W, there are mainly two types of IDLs observed along the following Burgers vectors (B): 1/2 <111> and <100> [12][13][14]. Based on previous investigations, both IDLs types likely start to occur at 150 keV of PKA [15].
However, W and its alloys have been subjected to many studies for their structural and mechanical properties exploration, while the strain is the topic of interest [16][17][18]. Several strain sources were observed during reactor operation, including the strain associated with void expansion generated by irradiation and the strain associated with solute segregation, which is also referred to as local strains. These local strains had a significant effect on the understanding of W's defects formation, such as a decrease in the threshold displacement energy (TDE) due to tensile strain and an increase in the TDE due to compression strain [19][20][21]. Due to the volumetric expansion in one direction, void-metal interactions have the potential to introduce strains of up to 5% [22,23]. Local strains may influence the number, type, and shape of defects generated, thereby altering the material's properties. Cascade in W and its alloys have been subjected to several studies by molecular dynamics (MD) simulations methods concerning the dynamics of cascade and defects production [24][25][26][27]. MD simulations can describe both the heat spike and the recombination of damage, and they are efficient enough to track the entire evolution of cascades in systems with millions of atoms.
Although high energetic cascade events have been proposed as a source of highly mobile IDLs, we believe that the creation probability of IDLs processed under local strain has not been observed in experiments or MD simulations for W.
We use an atomistic model to estimate the chances of IDLs development by irradiation under strain. Six strain magnitudes (−1.4, −0.8, −0.2, 0.4, 1.0, and 1.6 percent) were hydrostatically applied to each of the five distinct deposition energies of PKA (specifically, 1, 3, 6, 10, and 14 keV). Additionally, we examined the peak number of Frenkel pairs (FPs), which indicated the number of survived vacancies and self-interstitial atoms (SIAs) in the irradiated W structure. The number of FPs that survived the cascade's events, as well as the efficiency of defects production, was of interest.

Simulation Methods
For all MD simulations, the large-scale atomic/molecular massively parallel simulator (LAMMPS) code [27] was used in conjunction with the materials square platform [28]. An embedded atom method (EAM) potential [29] was used to characterize the interatomic interactions. This potential was originally parameterized by Derlet et al. [30] and refined by Björkas et al. [31] (DB) for recoil modeling, and it was used in this study.
It has been demonstrated that this EAM potential has adequate accuracy in predicting TDE and defect energies [19,[32][33][34] To begin the MD simulations, the simulation cell of a free deformation structure is first equilibrated with 30 K and 0 Pa under an NPT ensemble of the nose-hoover thermostat, which are the initial conditions that are employed before the strain is applied to the structure. Additionally, 30 K is within a typical temperature range for recoil calculations and for the TDE evaluation. The systems of 30 × 30 × 30, 45 × 45 × 45, 55 × 55 × 55, and 70 × 70 × 70 supercells of BCC, each containing 54,000, 182,250, 332,750, and 686,000 of W atoms (four intact structures), respectively, were designed utilizing the materials square platform. To calculate the dimension of each supercell's side, multiply it by the W lattice constant, which is 3.16 Å (i.e., 3.16 × 30 = 94.8 Å). Cell sizes were set to ensure that the system's average temperature after collision was less than 200 K on average, and that the displacement cascade would not cross cell boundaries. These four supercells were created to stand for the cascade energies of 1, 6, 10, and 14 keV. Next, each of the four structures was subjected into six strain values, each with an incremental hydrostatic strain of 0.6%, going from −1.6% to 1.4% (Table 1). The hydrostatic strain was applied in all directions, with negative values corresponding to compressive deformation and positive values corresponding to tensile deformation.
With a maximum displacement of 0.01 per step (x max ) and a maximum time step (t max ) of 0.02 ps, an adaptive time step was used. We determined that the settings were accurate enough to evaluate the formation of defects in comparison to numerous different x max and t max values after thorough examination and analysis. Each recoil MD simulation was initiated by imparting recoil energy to an atom located near the lattice center of the strained structures in one of four recoil directions (i.e., ), referred to as the primary knock-on atom (PKA). When recoil energy was introduced, it was converted to the velocity components of the PKA. Voronoi analysis [35] was utilized to identify whether or not a defect was generated: if the lattice site included two or more atoms, then an SIA was formed; if the site was empty, then a vacancy was formed.
The irradiation-induced displacement cascade was allowed to evolve for 20~30 ps for each PKA energy in a canonical NVT ensemble (i.e., constant atoms, volume, and temperature). The evolution time demonstrates that greater PKA energies required more evolution time to reach the stability state in comparison with the lower PKA energies, thus we confirmed that this setting was sufficient to observe the three ballistic phases.
The microcanonical/NVE ensemble was used for all cascade simulations (i.e., constant atoms, constant volume, and constant energy). Although the purpose of collision cascade simulation is to determine the number of FPs generated by the mean of MD, the standard number of atomic displacements can also be theoretically approximated using the Norgett-Robinson-Torrens (NRT) technique [33].
where E de is the nuclear deposition energy, which is about equivalent to the deposition energy collision carried by PKA for cascade, and E d is the TDE of the material, which is 90 eV according to International American Society for Testing and Materials (ASTM) standards [36] and our previous work [37]. The modeling parameters employed and the standard NRT displacements obtained are summarized in Table 2. After analyzing the final atomic configurations using OVITO postprocessing and visualization software [38], its modifiers were used to derive the cascade configuration and the clusters' sizes. The dislocation extraction algorithm (DXA) [39], which determines the burger vectors of interstitial dislocations and outputs as a line representation, was used to detect the presence of dislocations in a crystal.

Results and Discussion
The FP's production number were conducted over 16 independent MD simulations (4 different timings of each of the 4 different directions) at specific strain values of specific PKA energy at 30 K. It is noted that error bars included in the figures denote the standard error of the mean (SEM) and is determined as SEM = σ √ n , where σ is the standard deviation of the average FP over the set of n = 16 independent simulations. The error bars in the figures show this SEM. Because each strain value showed a similar thermal spike trend for all PKA energies, only the typical response of damage evolution over time for 10 keV and 14 keV are depicted in Figure 1a,b, respectively. Peak onset occurred earlier in compressed strains, whereas peaks reached greater heights and lasted longer in tensile strains. A comparison of the number of defects generated under each PKA energy and strain condition was included to clearly determine the effects of each PKA energy and strain type.

Results and Discussion
The FP's production number were conducted over 16 independent MD simulation (4 different timings of each of the 4 different directions) at specific strain values of specifi PKA energy at 30 K. It is noted that error bars included in the figures denote the standard error of the mean (SEM) and is determined as SEM = √ , where σ is the standard deviation of the average FP over the set of n = 16 independent simulations. The error bars in th figures show this SEM. Because each strain value showed a similar thermal spike trend for all PKA energies, only the typical response of damage evolution over time for 10 keV and 14 keV are depicted in Figure 1a,b, respectively. Peak onset occurred earlier in com pressed strains, whereas peaks reached greater heights and lasted longer in tensile strains A comparison of the number of defects generated under each PKA energy and strain con dition was included to clearly determine the effects of each PKA energy and strain type.  After the PKA received its kinetic energy, the number of atomic displacements in creased sharply during the ballistic phase and the highest number of FP was reached, bu most of the resulting defects moved to the original crystal lattice sites-a process known as the recombination or stabilization phase-with time. The peak time, which is the period between event initiations until the maximum number of created defects is reached, is th main feature of the ballistic stage. Different cascade core volumes were observed for var ious strain levels, resulting in a different number of defects being generated [40]. When After the PKA received its kinetic energy, the number of atomic displacements increased sharply during the ballistic phase and the highest number of FP was reached, but most of the resulting defects moved to the original crystal lattice sites-a process known as the recombination or stabilization phase-with time. The peak time, which is the period between event initiations until the maximum number of created defects is reached, is the main feature of the ballistic stage. Different cascade core volumes were observed for various strain levels, resulting in a different number of defects being generated [40]. When the W structure was compressed, the atoms became closer together, allowing the PKA energy to be dissipated over a smaller structure volume, resulting in fewer collisions without the generation of defects, as the PKA had a lower energy than the TDE [10]. Additionally, because the lattice distances between atoms are smaller in a compressed structure, the cascade atoms would move shorter distances. However, when the lattice distance between atoms increases, the atoms with more energy move freely and faster away from the cascade core, decreasing the possibility of recombination and generating more defects. The probability of secondary PKA occurrence increases with increasing PKA levels, and those secondary PKAs initiate the cascade with a lower energy level than the parent PKA.

PKA Energy and Direction Effects
Aside from anisotropic crystal structure and strain application, the chaotic nature of collision events (such as subcascade generation and channeling) [41] complicates primary radiation damage morphology and distribution. As a result, where the cascade procedure varies in stochastic nature from case to case, deterministic redistribution of collided PKA energy is impossible. As a result, we managed to run 16 different displacement cascades of each PKA recoil energy of each strained structure. The PKA energy and PKA recoil direction dependence were evaluated in terms of the number of FPs that remained after the relaxation time at the end of the cascade simulations, as shown in Figure 2.
Solids 2022, 3, FOR PEER REVIEW the W structure was compressed, the atoms became closer together, allowing the energy to be dissipated over a smaller structure volume, resulting in fewer collisions out the generation of defects, as the PKA had a lower energy than the TDE [10]. Add ally, because the lattice distances between atoms are smaller in a compressed stru the cascade atoms would move shorter distances. However, when the lattice distan tween atoms increases, the atoms with more energy move freely and faster away fro cascade core, decreasing the possibility of recombination and generating more d The probability of secondary PKA occurrence increases with increasing PKA level those secondary PKAs initiate the cascade with a lower energy level than the parent

PKA Energy and Direction Effects
Aside from anisotropic crystal structure and strain application, the chaotic nat collision events (such as subcascade generation and channeling) [41] complicates pr radiation damage morphology and distribution. As a result, where the cascade proc varies in stochastic nature from case to case, deterministic redistribution of collided energy is impossible. As a result, we managed to run 16 different displacement cas of each PKA recoil energy of each strained structure. The PKA energy and PKA direction dependence were evaluated in terms of the number of FPs that remained the relaxation time at the end of the cascade simulations, as shown in Figure 2.  Figure 2. While the results obtained for different dire varied by 8 FPs for strains less than 1%, the variance is relatively larger for larger s such as 1.6%. However, because the high-indexed direction <321> tends to avoid cha ing [41], whereas <111> has a higher probability of channeling effects, and one of the est numbers of FP generated would be with the <111> direction; consequently, <321 produce fewer defects due to its higher TDE.

Point Defects and Dislocation Loop Formation
There is a direct correlation between the volume changes of the applied hydro strain and the number of defects generated. The strained peak number of FP of al  Figure 2. While the results obtained for different directions varied by 8 FPs for strains less than 1%, the variance is relatively larger for larger strains such as 1.6%. However, because the high-indexed direction <321> tends to avoid channeling [41], whereas <111> has a higher probability of channeling effects, and one of the highest numbers of FP generated would be with the <111> direction; consequently, <321> may produce fewer defects due to its higher TDE.

Point Defects and Dislocation Loop Formation
There is a direct correlation between the volume changes of the applied hydrostatic strain and the number of defects generated. The strained peak number of FP of all PKA energies is plotted in Figure 3a, and each point in the figure is averaged over 16 different individual simulations for simplicity. The strained peak number of FPs tended to increase with tensile strain, most notably at 10 and 14 keV of PKA, and such effects on defect generation have been reported [41]. The number of defects FPs produced is directly proportional to the volume changes caused by the applied hydrostatic strain. The number of strained survived FPs after equilibration period is important for estimating radiation damage effects, particularly for understanding the subsequent processes of defect motion and clustering. Our reported FPs result nicely confirms previous MD work on tungsten by Wahyu et al. [42] at free strain conditions for 1 and 10 keV. The surviving number of FPs in Figure 3b varied significantly for the various PKA energies. The surviving number of FPs increased with increasing tension and decreased with compression strain for all PKA energies, due to the change of FP formation energy and the consequent TDE [33]. The variation in the survived number of FPs is more pronounced for higher PKA energies, which is consistent with previous findings [40,43].
Solids 2022, 3, FOR PEER REVIEW 7 individual simulations for simplicity. The strained peak number of FPs tended to increase with tensile strain, most notably at 10 and 14 keV of PKA, and such effects on defect generation have been reported [41]. The number of defects FPs produced is directly proportional to the volume changes caused by the applied hydrostatic strain. The number of strained survived FPs after equilibration period is important for estimating radiation damage effects, particularly for understanding the subsequent processes of defect motion and clustering. Our reported FPs result nicely confirms previous MD work on tungsten by Wahyu et al. [42] at free strain conditions for 1 and 10 keV. The surviving number of FPs in Figure 3b varied significantly for the various PKA energies. The surviving number of FPs increased with increasing tension and decreased with compression strain for all PKA energies, due to the change of FP formation energy and the consequent TDE [33]. The variation in the survived number of FPs is more pronounced for higher PKA energies, which is consistent with previous findings [40,43]. We found, with the DXA algorithm [39], that the relatively large shift in FPs formation at larger strains was due to the IDLs formation. The ½ <111> of <100> IDLs begin to appear at 10 keV of 1.0 % strain, whereas when the strain is increased to 1.6%, the IDLs We found, with the DXA algorithm [39], that the relatively large shift in FPs formation at larger strains was due to the IDLs formation. The 1 2 <111> of <100> IDLs begin to appear at 10 keV of 1.0 % strain, whereas when the strain is increased to 1.6%, the IDLs begin to appear at lower PKA as 8 keV, indicating that strain has a significant effect on IDL formation. At 1.6% strain, the probability of 1 2 <111> IDL formation is approximately two times that of <100> formation for the PKA from 8 keV, 9 keV, and 10 keV, while the ratio became nearly similar when the PKA is 14 keV. To confirm this, we evaluated the formation energy of both types of IDLs, and it is confirmed that as the number of SIA increases, the formation energy of the <100> clusters decrease, and this can tell why the more <100> IDLs as compared with the 1 2 <111> IDLs appeared at larger PKA and larger strains (Figure 4b). Note that this atomic potential for tungsten at high PKA stated that the <100> cluster was more stable than the <111> cluster, which is consistent with previous findings of the current potential model [15]. Solids 2022, 3, FOR PEER REVIEW 8 two times that of <100> formation for the PKA from 8 keV, 9 keV, and 10 keV, while the ratio became nearly similar when the PKA is 14 keV. To confirm this, we evaluated the formation energy of both types of IDLs, and it is confirmed that as the number of SIA increases, the formation energy of the <100> clusters decrease, and this can tell why the more <100> IDLs as compared with the ½ <111> IDLs appeared at larger PKA and larger strains (Figure 4b). Note that this atomic potential for tungsten at high PKA stated that the <100> cluster was more stable than the <111> cluster, which is consistent with previous findings of the current potential model [15].  It is also worth noting that the number and type of IDLs vary depending on vibration timing and collision displacement direction, as shown in Figure 5. The dislocation lines show that no embedded atoms interfere with the FP's point defects, and the defects formed are classified as monovacancy, mono SIA, cluster (two or more SIAs or vacancies), and dislocation loops. However, as PKA levels rise, the length of the dislocation loop It is also worth noting that the number and type of IDLs vary depending on vibration timing and collision displacement direction, as shown in Figure 5. The dislocation lines show that no embedded atoms interfere with the FP's point defects, and the defects formed are classified as monovacancy, mono SIA, cluster (two or more SIAs or vacancies), and dislocation loops. However, as PKA levels rise, the length of the dislocation loop appears to be extended (see Supplementary Materials Videos S1 and S2). Even if the kinetics differ from ab initio calculations [44] due to the absence of subcascades in tungsten, as opposed to iron [34], then the conclusion that local strain application leads to IDLs formation at lower PKA energies will not be affected by the different kinetics. Because of the decrease in TDE value at increased strain levels, IDL creation began to occur at the peak time of the ballistic phase. With more FPs, SIA clustering into a loop occurred with lower energy than monodefects ( Figure 6). The IDL of 10 SIAs has an energy that is 8 eV lower than the energy of 10 monodefects formation. The energy cutoff in this study for observing interstitial loops was 8 keV. This means that IDLs can occur even at lower PKA energies and lower PKA atomic masses, but only when the local strain is high enough, and this can be considered an addition to the study, which stated that a different types of IDLs can only be observed at very high PKA energies or larger PKA masses [41]. Because of the decrease in TDE value at increased strain levels, IDL creation began to occur at the peak time of the ballistic phase. With more FPs, SIA clustering into a loop occurred with lower energy than monodefects ( Figure 6). The IDL of 10 SIAs has an energy that is 8 eV lower than the energy of 10 monodefects formation. The energy cutoff in this study for observing interstitial loops was 8 keV. This means that IDLs can occur even at lower PKA energies and lower PKA atomic masses, but only when the local strain is high enough, and this can be considered an addition to the study, which stated that a different types of IDLs can only be observed at very high PKA energies or larger PKA masses [41].

Extra Analysis for the Point Defects and Production Efficiency
The defect production efficiency is defined as the ratio of theoretical methods for evaluating the primary radiation damage of materials to MD results in terms of FPs counts. The defect production efficiency is denoted by the abbreviation NFP/NNRT. Figure  7a,b show the deformed NRT displacements (i.e., the number of FPs predicted by the NRT theoretical model) and the value of Ed,j(V) from our previous work [45]. The defect production efficiency for strained conditions as a function of PKA energy; defect production efficiency decreased as PKA increased for all strain conditions except for 1.6% (4.8% volume change). The defect production efficiency becomes stable or slightly increased for higher PKA energies because of the formations of IDLs [46,47]. At higher PKA energies and higher tensile strain-mainly 14 keV at 1.6%, respectively-the dislocation loop formation led to further defect generation, as clarified in Section 3.2.1; hence, the efficiency increases again.

Extra Analysis for the Point Defects and Production Efficiency
The defect production efficiency is defined as the ratio of theoretical methods for evaluating the primary radiation damage of materials to MD results in terms of FPs counts. The defect production efficiency is denoted by the abbreviation N FP /N NRT . Figure 7a,b show the deformed NRT displacements (i.e., the number of FPs predicted by the NRT theoretical model) and the value of E d,j (V) from our previous work [45]. The defect production efficiency for strained conditions as a function of PKA energy; defect production efficiency decreased as PKA increased for all strain conditions except for 1.6% (4.8% volume change). The defect production efficiency becomes stable or slightly increased for higher PKA energies because of the formations of IDLs [46,47]. At higher PKA energies and higher tensile strain-mainly 14 keV at 1.6%, respectively-the dislocation loop formation led to further defect generation, as clarified in Section 3.2.1; hence, the efficiency increases again.

Extra Analysis for the Point Defects and Production Efficiency
The defect production efficiency is defined as the ratio of theoretical methods for evaluating the primary radiation damage of materials to MD results in terms of FPs counts. The defect production efficiency is denoted by the abbreviation NFP/NNRT. Figure  7a,b show the deformed NRT displacements (i.e., the number of FPs predicted by the NRT theoretical model) and the value of Ed,j(V) from our previous work [45]. The defect production efficiency for strained conditions as a function of PKA energy; defect production efficiency decreased as PKA increased for all strain conditions except for 1.6% (4.8% volume change). The defect production efficiency becomes stable or slightly increased for higher PKA energies because of the formations of IDLs [46,47]. At higher PKA energies and higher tensile strain-mainly 14 keV at 1.6%, respectively-the dislocation loop formation led to further defect generation, as clarified in Section 3.2.1; hence, the efficiency increases again.   [45]. (b) T duction efficiency which is defined as the ratio between MD and NRT result for def These findings imply that tensile strain conditions, particularly with h ergies, should be avoided for irradiated structures, and that if local defects a accumulate in structures, the structure should be treated with caution as local strain. Furthermore, our data, specifically the formation of IDLs, can comprehensive database of cascades obtained from MD for object kinetic (OKMC) simulations of primary cascade damage in bulk tungsten.  [45]. (b) The defects production efficiency which is defined as the ratio between MD and NRT result for defects calculation.
These findings imply that tensile strain conditions, particularly with higher PKA energies, should be avoided for irradiated structures, and that if local defects and impurities accumulate in structures, the structure should be treated with caution as the source of local strain. Furthermore, our data, specifically the formation of IDLs, can be used as a comprehensive database of cascades obtained from MD for object kinetic Monte Carlo (OKMC) simulations of primary cascade damage in bulk tungsten.

Conclusions
Using molecular dynamic simulation, we studied the formation of interstitial dislocation loops caused by strain in tungsten. PKA energies of 1, 6, 10, and 14 keV were applied at 30 K on a volumetric strained structures (hydrostatic) of six different strain values. Our findings suggest that the number of FPs that survive, because structure deformations increase with tensile strain and decrease with compression strain. Additionally, according to the DXA analysis, the IDLs begin to cluster at higher strains of 1.0 and 1.6%, starting from 7 keV of PKA. According to the findings, strain caused a decrease in the barrier to the occurrence of IDLs and thus material degradation. Overall, the findings of this study indicate that tensile loading should be avoided in materials used in irradiation environments, whereas compression deformation could be a promising factor in the design of structural materials used in nuclear reactor applications.