Effect of Carbon on Void Nucleation in Iron

The study reports the significance of carbon presence in affecting void nucleation in Fe. Without carbon, void nucleation rates decrease gradually at high temperatures but remain significantly high and almost saturated at low temperatures. With carbon present, even at 1 atomic parts per million, void nucleation rates show a low-temperature cutoff. With higher carbon levels, the nucleation temperature window becomes narrower, the maximum nucleation rate becomes lower, and the temperature of maximum void nucleation shifts to a higher temperature. Fundamentally, this is caused by the change in effective vacancy diffusivity due to the formation of carbon-vacancy complexes. The high sensitivity of void nucleation to carbon comes from the high sensitivity of void nucleation to the vacancy arrival rate in a void. The void nucleation is calculated by first obtaining the effective vacancy diffusivity considering the carbon effect, then calculating the defect concentration and defect flux change considering both carbon effects and pre-existing dislocations, and finally calculating the void nucleation rate based on the recently corrected homogeneous void nucleation theory. The study is important not only in the fundamental understanding of impurity effects in ion/neutron irradiation but also in alloy engineering for judiciously introducing impurities to increase swelling resistance, as well as in the development of simulation and modeling methodologies applicable to other metals.


Introduction
Subjected to ion, neutron, or electron irradiation, metals undergo the buildup of point defects in concentrations markedly surpassing their equilibrium levels.Consequently, small defect clusters nucleate and grow, evolving into complex defect structures like dislocations, dislocation loops, and voids.Understanding the response of metals to irradiation holds the utmost importance for both fission and fusion reactors, as irradiation causes material degradation [1].Typical material degradation in fission reactors includes stress corrosion cracking [2], swelling [3], fuel restructuring [4], delayed hydride cracking [5], fracture [6], embrittlement [7], and creep [8].Typical degradation expected in fusion reactors includes cracking, flaking, blistering, and embrittlement [9,10].These phenomena greatly limit reactor performance and reactor lifetime.Concerning pressured water reactors and fast neutron reactors, the void swelling raises a substantial challenge regarding reactor safety.The process of swelling in irradiated metals typically includes three phases: an initial incubation period, a transient phase, and a phase of steady-state growth.While rich data exist for the steady-state growth phase, understanding of the first phase within the incubation period is very limited.Primarily, this is due to the challenge of identifying void embryos that fall beneath the detection threshold of conventional microscopy methods.Research concerning austenitic steels has revealed that the initial incubation phase of swelling is highly sensitive to slight variations in composition and minor deviations in thermal-mechanical treatments [11].Furthermore, the duration of the incubation period is significantly influenced by environmental conditions, including the rate of displacement per atom (dpa), temperatures, and stresses.Gaining insights into void nucleation during the swelling incubation period is important for material development toward radiation-tolerant high-performance steels.
Among the numerous factors influencing void nucleation, the impact of impurities, especially carbon, holds significance.Recent studies indicate that a carbon concentration of ~100 atomic parts per million (appm) can alter the void distribution profile in irradiated iron [12].Given that carbon is a common contaminant and is often utilized as an alloying element, whether intentionally or unintentionally, it becomes imperative to investigate its effects.The impact of carbon on interstitials is less pronounced when compared to its effect on vacancies [13][14][15].Carbon can form bonds with a <110> Fe-Fe interstitial dumbbell, exhibiting a weak binding energy of approximately 0.5 eV [13].Concerning vacancies, carbon exerts its influence through two possible mechanisms.First, carbon can elevate the migration energy of vacancies without trapping vacancies [16].Second, carbon can bond with vacancies, resulting in the formation of C-V complexes [17].The first mechanism influences the diffusivity of vacancies without altering their concentration.Conversely, the second mechanism doesn't affect the diffusivity of free vacancies but reduces the number of mobile vacancies.Both effects can be represented by introducing effective vacancy diffusivity.The second effect is more significant and representative and can be modeled by considering various complex configurations beyond the interaction between a single vacancy and a single carbon atom.
The effect of carbon as a contaminant during prolonged accelerator ion irradiation has been a well-known issue.The problem is serious enough that accelerator laboratories involved in nuclear material testing have developed various methods to mitigate it.The accelerator laboratory at Los Alamos National Laboratory uses a coating technique [18], the accelerator laboratory at the University of Michigan employs plasma cleaning [19], and the accelerator laboratory at Texas A&M University uses multiple beam deflection [20].Universal contamination may come from the beam line, target chamber, or the sample itself.Most likely, it is a combination of all, but the primary contamination source may vary for specific laboratories.Contamination on the sample surface, in one situation, can come from surface preparation, shipping, or pasting on the target holder.
The effect of carbon needs a thorough investigation.Additionally, there is a lack of a quantitative method to predict the level of severity under various carbon levels.Studying the effect of carbon on void swelling requires a multiscale modeling approach.The first step is to evaluate the effect of carbon on defect kinetics, which requires atomic-scale modeling.The second step is to use defect kinetics to predict void swelling.This step involves obtaining the defect population and defect flux under a given microstructure (i.e., pre-existing dislocation density).
Here, the sensitivity of void nucleation in α-iron (Fe) to carbon background levels is studied, driven by the need to evaluate the tolerance of impurity levels in both irradiation testing and material synthesis.Fe is selected for its representation of Fe-based steels.The former pertains to the credibility of accelerated irradiation testing for emulating neutron environments, while the latter pertains to the qualification of emerging manufacturing techniques such as additive manufacturing.

Modeling Procedure
Void nucleation theory was previously established by Katz and Wiedersich [21] and Russel [22].The mathematical treatment was recently revisited by Shao, and modification was proposed to fix a problem in the original treatment [23].The original treatment uses a specific selection of energy reference points to simplify the expression.It, however, mistakenly removes an exponential modification function, leading to an underestimation of nucleation rates by orders of magnitude.Below, the original theory and recent corrections are explained as the foundation of the modeling procedure in the present work.
In homogenous void nucleation theory [21], the growth rate of voids from size x to the next size (x + 1) for an arbitrary size function f (x, t) considers the competing effects of four different interactions of voids with point defect, as described by the following: where J(x, t) is the void nucleation rate from size x to x + 1, x is the number of vacancies contained in a void, t is time, β v and β i are the arrival rates of vacancies and interstitials to a void per unit area and unit time, respectively.γ v and γ i are the emission of vacancies and interstitials by a void per unit area and unit time.s(x) is the surface area of a void of size x.
Following the theoretical framework established by Katz,Wiedersich,and Russel [21,22], the void nucleation rate under a steady state growth is given by the following: where n(x) is a specific constrained void size distribution satisfying J = 0. Additional details of the derivation are provided in Appendix A. The J value from Equation (2) explicitly depends on vacancy concentration, as reflected by the appearance of β v .Additionally, the value implicitly depends on interstitial concentrations, with their effect included in n(x).The function n(x) represents a specific size distribution of zero net growth flux.J = 0 is achieved through the cancellation of flux contributions from both vacancy-driven and interstitial-driven growth [21].n(x) has the solution, as follows: Additional details of the derivation are provided in Appendix A. The size distribution n 0 (x) represents another specific size function at which voids are in equilibrium with vacancies.The correction factor P is primarily determined by mathematical deviations, but it does possess a physical interpretation.P represents the ratio of the magnitude of the flow of void embryos when vacancies are in equilibrium with interstitials to the magnitude of the flow when actual interstitials are present.The magnitude of the flow is composed of the combined effects of vacancy-absorption-induced growth from size x to x + 1 and interstitial-absorption-induced reverse growth from size x + 1 to x, regardless of the flux direction.n 0 (x) is also called the constrained equilibrium distribution of voids [21].It means the formation of a void size x from x vacancies does not change the total free energies.n 0 (x) was derived as the following [23]: where N is the atomic density of the system, k is the Boltzmann constant, T is temperature, and S is defined as vacancy supersaturation ratio (S≡ n 1 n eq 1 ). n eq 1 is the equilibrium number of vacancies per unit volume.W(x) is work needed to create a void of size x.
Equation ( 4) is different from the original formula in reference [21], which is expressed as follows: The problem with Equation (5) in the original theory comes from the selection of a reference state.The theory was built upon the assumption that the system contains voids and vacancies reach equilibrium, in which there is no chemical potential change when x vacancies form a void of size x.In order to simplify the expression, Katz and Wiedersich set the chemical potential of a vacancy in the system containing n eq 1 as zero [21].Consequently, the effect of chemical potential disappears.In contrast, the derivation shows chemical potentials change the base that is raised to the power of x in the expression of n 0 (x), introducing a size-dependent.In the correct equation, Equation (4), the effect of chemical potential is included in n eq 1 .More details of the derivation can be found in reference [23].
In a different way to explain, the original treatment ignores the effect of entropy in forming a vacancy, which is supposed to contribute a temperature-independent but size-dependent term in n 0 (x).This effect is included in in Equation ( 4).
The energy required to form a void in pure Fe, W(x), can be obtained using various methods, including continuum mechanics models [21], first-principle quantum mechanics calculations [24], and molecular dynamics simulations (MD) [25].The present study selected MD-obtained results for the work needed to create a void of size x, as given by the following [25]: The unit of W(x) in Equation ( 6) is eV.The selection of MD-obtained results considers the fact that void swelling is frequently observed in pure Fe at temperatures around 500 • C [26].MD-obtained kinetics yield a swelling temperature range in good agreement with experimental observations, as will be shown in this study.
Ab initio calculation suggests that the vacancy migration energy in body-center cubic (BCC) Fe is 0.73 eV [27].If the vacancy migration energy is set at 0.73 eV and the vacancy formation energy at 2.59 eV [25], the activation energy for Fe self-diffusion in BCC Fe is about 3.32 eV.Note that this value matches the modeling result [27].This value is also reasonably close to the experimental values ranging from 2.95 eV to 3.10 eV [28,29].As a summary, Table 1 lists the parameters of pure Fe used in the present study.The preexponential factor of the self-diffusion coefficient, A SD , is 11.75 cm 2 /s [27].The vacancy formation entropy ∆S v is 2.17 k/vacancy [30].2.17 [30] Self-diffusion coefficient (D SD = A SD exp − H SD kT is introduced to increase the accuracy in the calculation of β v , the arrival rates of vacancies to a void per unit area and unit time.For self-diffusion dominated by the monovacancy diffusion mechanism, D SD is given by the following: where f v is the diffusion correlation factor when the sequential jumping directions are correlated.Correlation is always expected when a trace atom diffuses via the vacancy diffusion mechanisms.The most commonly used f v value is 0.72722 for the monovacancy diffusion mechanism in the BCC crystal structure [31,32].D V is the vacancy diffusivity.c eq v is the vacancy atomic fraction at equilibrium.D v is given by the following: where v is the number of successful jumps per second, and λ is the jumping distance.D v is related to β v (vacancy flux to a void surface).β v is given by the following: where c v is the vacancy atomic fraction under irradiation, and N is the substrate atomic density (8.482 × 10 22 /cm 3 for Fe) [33].The product of c v , N, and λ gives the number of vacancies in a volume defined by a unit area and the thickness of one jumping distance.
The vacancy supersaturation ratio is given by the following: where C v and C eq v are the vacancy concentration and the vacancy equilibrium concentrations, respectively.Combining Equations ( 7)- (10), one obtains the following: Using D SD to obtain β v increases the accuracy, since the attempt jumping frequency, which is the pre-exponential factor of v in Equation (9), is difficult to calculate.Another advantage is that S appears as a variable in the expression.

Void Nucleation at Arbitrary T, S, and
As a summary of the results obtained using Equations ( 4) and ( 5), Figure 1a-c shows n(x) distribution changes by varying temperatures (T ), vacancy supersaturation ratios S, and the ratios of interstitial-to-vacancy arrival rates to a void (β i /β v ).All n(x) curves feature a dip at the critical size, x c .With increasing T, the size curves move upwards, exhibiting smaller x c and higher n(x c ) (Figure 1a).Similarly, with increasing S values, the curves shift upward, reducing x c and increasing n(x c ) (Figure 1b).Increasing β i /β v increases x c while reducing void density at large sizes (Figure 1c).
Once n(x) is obtained, Equation ( 2) is used to calculate void nucleation rate J. Figure 2 plots the J values as a function of S for different T and β i /β v .The J values are very sensitive to all three parameters.For example, at 900 K and β i /β v = 0.99, a change of S from 1 × 10 4 to 1 × 10 5 leads to an increase in J by more than six orders of magnitude.Figure 2 also shows a compensating effect among parameters.For instance, to maintain the same J value, the impact of lowering T can be counterbalanced by increasing S.

Void Nucleation Considering Irradiation for Defect Production and Dislocations as Defect Sinks
After establishing the relationship of n(x) with T, S, and β i /β i , the study proceeds to rate theory calculations to establish how these parameters relate to the displacements per atom K 0 and defect sinks.Two master equations describe the time rate of change of vacancy concentration (C v ) and interstitial concentration (C i ), expressed as follows [1]: where t is time, f survive is the survival fraction of defects after the initial damage creation, N is the atomic density of Fe, K 0 is the displacements-per-atom (dpa), K th v and K th i are thermal generation rates of vacancies and interstitials, respectively, K ⊥(v) and K ⊥(i) are sink strength for vacancies and interstitials, respectively, ρ v and ρ i are sink densities for vacancies and interstitials, respectively, K iv is the interstitial-vacancy recombination rate, D v and D i are diffusivities of vacancies and interstitials, respectively.
Under equilibrium conditions, without ion irradiation and diffusion, the above two equations lead to the expression of the thermal generation rates. ) Substituting the above two equations into Equations ( 12) and ( 13), and assuming a steady state where ∂C v ∂t = 0 and ∂C i ∂t = 0, the following equation can be easily derived.The equation can also be obtained, without mathematics, but based on the physics that the net point defect fluxes of interstitials and vacancies to defect sinks are equal.Otherwise, a steady state cannot be reached [34].
Combining Equations ( 12)-( 16) for a steady state under irradiation, a quadratic equation is derived for vacancy concentration deviation from its equilibrium ∆C v = C v − C eq v , expressed as follows: The solution is as follows: Substituting C v obtained from Equation (18) into Equation ( 16) obtains C i .Then the ratio β i /β v can be obtained using The dislocation sink strength K ⊥ is described by the following: where ρ ⊥ is dislocation density, and r ⊥(i,v) is the defect-trapping radius for trapping interstitials and vacancies.The point defect recombination rate is calculated by the following [1]: where Ω is atomic volume of one lattice atom and a is the lattice constant of Fe.Table 2 summarizes the parameters used in the present study.The black square lines in Figure 3 show the calculated vacancy atomic fraction concentration as a function of temperature at various dpa rates.Under a given dpa rate, the vacancy concentration exhibits a V-shaped pattern.At very high temperatures, the concentration follows C eq v .As the temperature decreases, C v initially decreases and then begins to rise at a certain temperature point and deviates from C eq v .At sufficiently low temperatures, where C eq v and C eq i are negligible, C v becomes proportional to √ K 0 .Changing K 0 by one order of magnitude results in a nearly parallel shift of the logarithmic plot of C v at low temperatures.This shift is obvious when comparing the square lines for K 0 ranging from 1 × 10 −2 dpa/s to 1 × 10 −8 dpa/s.Lowering K 0 causes the temperature point at which C v deviates from C eq v to shift to lower temperatures accordingly.In Figure 3, the calculations assume a dislocation density of 1 × 10 10 /cm 2 .Changing the dislocation density does not cause a significant shift in the curve, unlike K 0 .Instead, the dislocation density mainly influences the concentration near the turning point.A higher dislocation density results in a lower defect concentration and shifts the turning point to a lower temperature.The effect of dislocation density becomes more evident under lower dpa rate irradiation.The color bands in Figure 3 represent the contour map of the void nucleation rate as a function of vacancy concentration and temperature.Each color or contour line corresponds to a constant void nucleation rate.It exhibits a V-shaped pattern.For lower void nucleation rates, the minimum point of the V-shaped pattern shifts to lower defect concentrations and temperatures.The V-shaped behavior primarily arises from the change in the critical void size (x c ).When x c is larger, the void density at that size, n(x c ), is lower, resulting in lower void nucleation rates.In the high-temperature region (i.e., >900 K), x c decreases with increasing temperature, leading to higher nucleation rates at higher temperatures.The effect of irradiation is less significant in this temperature range because C eq v is already high, and the irradiation-induced S changes become insignificant.At extremely high temperatures, S approaches 1.In the low-temperature region (i.e., T < 600 K), the effect of S becomes more significant as C eq v is very low.Although lowering the temperature would tend to increase x c , the high value of S counteracts this and reduces x c .The significant changes in S alter the trend, resulting in the V-shaped temperature dependence as shown in Figure 3.
For two points on the contour line with the same nucleation rate, the following observations were made: Taking the example of a nucleation rate of 1 × 10 18 voids/cm 3 per second, when the temperature decreases from 700 K to 600 K, the significantly increased vacancy concentration (and subsequently, S values) shift the entire n(x) profile upwards, resulting in higher void densities at each size.Conversely, β v is reduced when the temperature changes from 700 K to 600 K, due to large changes in vacancy diffusivity.This reduction in β v compensates for the increase in void density, thereby maintaining the same nucleation rate.

Void Nucleation Considering Irradiation, Dislocation, and Carbon Incorporation
Carbon is well-known for causing the suppression of void swelling, but no quantitative evaluation has been established to assess its impact, which motivates the present study.The effect was not significant during the short-term ion irradiation of materials in early days but has become significant for prolonged irradiation of current advanced alloys, which are highly swelling-resistant.More recent relevant studies show that carbon is the reason why the width of the void-denuded zone is governed by an activation energy that significantly deviates from the expected vacancy migration energies [12].
C atoms temporarily immobilize vacancies through the formation of V-C n complexes.The complexes are expected to be dissociated later, and the dissociation probability is determined by the V-C binding energies.The binding energies of various C-V complexes were previously calculated using ab initio calculations [17].The binding energies E b are 0.41 eV for VC, 1.18 eV for VC 2 , and 1.30 eV for VC 3 .The complex concentration, expressed as the atomic fraction concentration, ∼ C V m C n , of various vacancy-carbon complexes V m C n is approximated by the mass-action law [17]: where m is the number of vacancies and n is the number of carbon atoms in a vacancycarbon complex.E b is the binding energy of the complex.Under the approximation that (1) the amount of carbon bonded with vacancies in complexes is significantly less than the total carbon dissolved in the system, and (2) the major complexes consist of those containing one vacancy (m = 1) and multiple carbons with n = 1, 2, and 3, the effective vacancy diffusivity D v e f f can be calculated by the following: where D v is the intrinsic diffusivity of vacancies, represents the atomic fraction concentration of free/isolated vacancies, ∼ C C is the atomic fraction concentration of carbon, and E bn is the binding energy of vacancy-carbon complex V 1 C n with n = 1, 2, and 3. Figure 4 plots the effective vacancy diffusivity as a function of temperature for carbon concentrations at 1, 100, and 10,000 appm levels.The solid line represents C-free Fe, exhibiting a single activation energy (0.73 eV).With the addition of carbon, the effective diffusivity deviates from the solid line at temperatures below a critical temperature.At higher C levels, this deviation occurs at higher temperatures.Notably, in the affected-temperature region, the effective vacancy diffusivity exhibits roughly the same activation energy of 1.91 eV, regardless of the carbon level.The diffusivity curve is parallelly shifted downward with increasing carbon levels.The constant activation energy of 1.91 eV for all C levels can be approximated as the sum of the migration energy (0.73 eV) of free vacancies and the binding energy (1.18 eV) of VC 2 complexes [17].The effect of VC 2 is dominant.The VC complex plays a less significant role due to its much lower E b value, and VC 3 is also less significant due to its relatively lower concentrations [17].Figure 4 also marks the temperatures at which the effective vacancy diffusivities begin to deviate from the carbon-free case.As discussed soon, these critical temperature points also play an important role in determining the temperature dependence of void nucleation.The method of calculating effective vacancy diffusivity in the present study follows the early approach of Fu et al. [17].The diffusivity calculation in Reference [17] used a vacancy migration energy of 0.67 eV.In the plot of Figure 4, an activation energy of 0.73 eV is used [27], to be consistent with the calculations in the preceding sections.
Substituting D V e f f into rate theory calculations to determine S, β v , and β i /β v , and then substituting these parameter values into the void nucleation rate calculation, one obtains J as a function of both carbon concentration and temperatures.As shown in Figure 5, introducing carbon dramatically changes the contour map.At zero carbon concentration, the nucleation rate remains high at temperature below ~700 K.However, with a small amount of carbon addition, nucleation rates quickly evolve into a temperature-dependent peak.As the carbon concentration increases, the peak height becomes lower, the peak width is narrower, and the peak center slightly shifts to a higher temperature.In Figure 6, the void nucleation rates are plotted as a function of temperature while fixing carbon concentrations at distinct levels of 0, 1, 10, 100, and 1000 appm.At zero carbon concentration (dot line), J remains above 3 × 10 17 /cm 3 per second when temperature is at 600 K and below.Even with carbon levels as low as 1 appm, void nucleation begins to show a peak.The nucleation rates start to drop at temperatures ~530 K and below.At a higher carbon level of 10 appm, the nucleation allowable temperature window is narrower, and the temperature of the maximum void nucleation rate increases to about 600 K. Additionally, at this carbon level, the peak nucleation rate is noticeably lower than in the carbon-free case.As shown by two arrows (one at 530 K for 1 appm carbon and the other at 620 K for 10 appm carbon), the temperatures at which void nucleation begins to deviate from the carbon-free case correspond to the temperatures at which effective diffusivities of vacancies begin to deviate from the carbon-free case (as marked in Figure 4).At a higher carbon level of 100 appm, the peak shifts to about 670 K, and the peak height is about 25% of that in the carbon-free case at the same temperature.At the highest carbon level of 1000 appm, the nucleation peak is reduced by more than two orders of magnitude compared to the carbon-free curve at the same temperature.
One way to validate the carbon effect on effective vacancy diffusivity is to measure the width of the void-denuded zone, denoted by ∆x.According to the analytical solution of rate theory equations, ∆x ∝ (D V /K) 1/4 [38].Recent studies on Fe irradiated by self- ions at various temperatures, beam energies, and dpa rates have measured an effective vacancy activity energy of 1.65 eV for single crystal Fe containing a carbon background concentration of about 100 appm [12].The experimentally extracted value of 1.65 eV aligns with the 1.91 eV predicted from Figure 4 and is significantly larger than the 0.73 eV predicted for carbon-free cases.
Validation of void nucleation in carbon-free Fe is indeed very challenging due to various factors, including the difficulty of manufacturing and irradiation testing.Carbon contamination has been an issue in prolonged accelerator-based ion irradiation testing, and void disappearance has been frequently reported [39][40][41].Carbon contamination is not expected in reactor irradiation testing.However, parameters such as dpa rate, temperature, and neutron flux are all interlinked and coupled.Fixing other parameters and allowing temperature to be the single variable proves to be challenging in reactor experiments.The prediction of a non-zero void nucleation rate in carbon-free Fe at low temperatures is a subject that will be of interest for future studies, especially if low-temperature neutron irradiation can be achieved.Testing low-temperature void nucleation is challenging in accelerator-based irradiation, with complexities extending beyond beam contamination.In the case of acceleratorbased ion irradiation, heavy ions suffer from the injected interstitial effect [42][43][44], which exhibits strong temperature dependence [44].The swelling suppression by the injected interstitials becomes more significant at lower temperatures [44], leading to a suppression behavior similar to that of carbon.Proton irradiation can minimize the injected interstitial effects to a large extent.However, proton irradiation can introduce local beam heating, making it difficult to test temperatures lower than 600 K. Commercial tandem accelerators used in university laboratories typically can obtain proton beam currents at a few microamperes.For advanced ion sources, beam currents up to tens of microamperes can be achieved.With typical proton beam energy at a few MeV, the amount of energy deposition is significant.Using a weak proton beam is not realistically helpful since a strong beam is needed to achieve a reasonable damage level.
The findings presented in this study can provide insights into the significant data scattering observed in additively manufactured (AM) steels.Previous reports have indicated that AM steels (316 L) exhibit reduced void swelling compared to wrought materials [45,46].However, it is noteworthy that void swelling behavior can vary significantly among AM steels manufactured by different groups.One likely cause of such variation is the difficulty in controlling impurities during the additive manufacturing process.These impurities include elements such as carbon, oxygen, and nitrogen, as well as ambient gases used during the printing process.The presence and concentration of these impurities can have a considerable impact on the void-swelling behavior of AM steels.
The high sensitivity of void swelling to carbon levels suggests that it is possible to alloy steels with a small amount of carbon.This can be done at a level that does not significantly alter the optimized mechanical properties but is sufficient to improve the resistance to void swelling.Indeed, nitrogen is another element that can have a similar effect on void swelling in steels.For instance, nitrogen-doped full-ferritic HT-9 exhibits much less swelling compared to conventional .
The present study is useful in the development of modeling methodologies.As explained in the introduction, the effect of carbon contamination is well known but has never been quantitatively evaluated.This study provides a method to analyze the carbon effect on void nucleation.In void swelling, nucleation during the incubation period is the most significant stage for determining the radiation tolerance of alloys.The proposed method includes the sequential steps of (a) calculating the effective diffusivity of vacancies in the presence of carbon at different levels; (b) calculating defect fluxes and defect supersaturation once effective diffusivity is known from step (a) and the density of pre-existing dislocations is known; and (c) calculating void nucleation rates.This approach can be applied to other metals.The required defect kinetics and carbon-vacancy complex bonding energy can be calculated from atomic-scale simulations.

Conclusions
The nucleation of voids in alpha-iron under particle irradiation and the influence of carbon were investigated through a combination of rate theory and nucleation theory calculations.The major finding is that the steady rates of void nucleation exhibit high sensitivity to carbon.Even a seemingly negligible carbon concentration, as low as a few appm, can dramatically reduce the nucleation rates and narrow the temperature window for nucleation.Such high sensitivity to carbon impurities comes from the effectiveness of carbon in trapping vacancies and reducing their effective diffusivity, and the importance of vacancy diffusivity in determining void nucleation rates.This study highlights the challenges in certifying and qualifying nuclear materials for use in reactor environments.It emphasizes the necessity for precise control of impurity levels during material processing and stringent requirements for irradiation testing, which are susceptible to contamination.Furthermore, the study provides an explanation for the significant data scattering observed in void swelling behaviors as reported in the literature.

Figure 1 .
Figure 1.Void density n(x) as a function of void size (number of vacancies contained in a void) under various conditions in α-iron: (a) different temperatures, (b) different S values, and (c) different

Figure 2 .
Figure 2. Void nucleation rates as a function of S at different temperatures (700, 800 K, 900 K) and

) 10 10 1 Figure 3 .
Figure 3. Vacancy concentration (atomic fraction, on a logarithmic scale) at different temperatures for various steady-state void nucleation rates (color contour) and under different dpa rates (square lines) in α-iron.The point of intersection between a square line and a colored contour line gives the nucleation rate at a specific temperature and dpa rate.

Figure 4 .
Figure 4. Effective vacancy diffusivity as a function of temperatures and C concentrations in α-iron.

Figure 5 .
Figure 5.The map of void nucleation rates as a function of C concentrations and temperatures in α-iron.

Figure 6 .
Figure 6.The plot of void nucleation rates as a function of temperature for C concentrations ranging from 0 to 1000 appm in α-iron.

Table 1 .
Parameters used in the present study for calculating void nucleation rates.

Table 2 .
Parameters used in the present study for rate theory calculations.