Microstructure-Dependent Rate Theory Model of Radiation-Induced Segregation in Binary Alloys

Abstract

Conventional rate theory often uses the mean field concept to describe the effect of inhomogeneous microstructures on the evolution of radiation induced defect and solute/fission product segregation. However, the spatial and temporal evolution of defects and solutes determines the formation and spatial distribution of radiation-induced second phase such as precipitates and gas bubbles/voids, especially in materials with complicated microstructures and subject to high dose radiation. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed to investigate the effect of inhomogeneous thermodynamic and kinetics properties of defects on diffusion and accumulations of solute A in AB binary alloys. Four independent concentrations: atom A, interstitial A, interstitial B, and vacancy on [A, B] sublattice are used as field variables to describe temporal and spatial distribution and evolution of defects and solute A. The independent concentrations of interstitial A and interstitial B allow to describe their different generation rates, thermodynamic and kinetic properties, and release the assumptions of interstitial generation and sink strength used in the conventional rate theory. Microstructure and concentration dependent chemical potentials of defects are used to calculate the driving forces of defect diffusions. With the model, the effects of defect chemical potentials and mobilities on the RIS in polycrystalline AB model alloys have been simulated. The results demonstrated the model capability in predicting defect evolution in materials with inhomogeneous thermodynamic and kinetic properties of defects. The model can be extended to materials with complicated microstructures such as a wide range of grain size distribution, coating structure and multiphases as well as radiation-induced precipitation subject to severe radiation damage.

Introduction

Radiation-induced segregation (RIS) and precipitation are important material property degradation mechanisms in irradiated materials (Lam et al., 1978; Odette and Lucas, 1986; Farrell et al., 1994; Akamatsu et al., 1995; Lu et al., 2008; Was et al., 2011; Hu et al., 2012; Nastar and Soisson, 2012; Field et al., 2013; Wharry and Was, 2013; Yang et al., 2016). In Ferritic/Martensitic (F/M) steels, experimental results suggested that the RIS of Cr at grain boundaries (GBs) enhances Cr precipitation, which not only led to embrittlement but also altered the corrosion resistance of grain boundaries (Bruemmer et al., 1999; Simonen and Bruemmer, 1999; Nastar and Soisson, 2012). Moreover, experimental data indicated that a number of factors such as temperature, grain boundary structure, and radiation rate might affect whether solutes segregate or deplete on grain boundaries (Damcott et al., 1995; Lu et al., 2008; Was et al., 2011). For instance, the RIS of Cr showed a bell-shape temperature dependence, changing from Cr enrichment to Cr depletion at grain boundaries with increasing temperatures (Wharry and Was, 2013). Low angle grain boundaries (LAGB) exhibited a suppressed RIS response when compared to relatively high angle grain boundaries (HAGB), and at the special Sigma coincident site lattice (CSL) boundaries, the RIS of Cr was also suppressed (Hu et al., 2012; Field et al., 2013).

In different length and time scales, theoretical models including rate theory (Woo and Singh, 1992; Grandjean et al., 1994; Allen and Was, 1998; Golubov et al., 2001; Was, 2016), cluster dynamics (CD) (Faney and Wirth, 2014; Jourdan et al., 2014; Ke et al., 2018), kinetic Monte Carlo (KMC) (Ortiz, 2007; Huang and Marian, 2016), atomistic kinetic Monte Carlo (AKMC) (Senninger et al., 2016; Soisson and Jourdan, 2016), Object kinetic Monte Carlo (OKMC) (Souidi et al., 2006; Dunn et al., 2013), and phase-field approaches (Badillo et al., 2015; Piochaud et al., 2016; Li et al., 2017) have been developed to describe defect evolution and RIS. AKMC takes into account all atoms in the system which could more accurately describe and calculate the effect of defect interaction on defect diffusion and clustering. But it is challenging to apply it in large simulation domains and high radiation doses. Coarse grained methods such as KMC and OKMC only consider defects in the system and are therefore computationally more efficient. However, to describe the complicated thermodynamic and kinetic properties of defects and defect clusters such as the elastic interaction, reaction, and anisotropic mobility requires a large set of model parameters which might be coupled and make the simulation difficult. The mean field approaches such as rate theory and CD treat structural defects such as dislocations, grain boundaries and second phase particle interfaces as sinks. With the effective sink strength of structural defects, the temporal evolution of average defect concentrations over the representative volume can be obtained efficiently by solving one dimensional rate theory (ODE). Spatially dependent rate theory has been developed to consider the effect of free surface on diffusion, trapping and detrapping kinetics of helium in BCC iron foils (Xu and Wirth, 2009; Xu and Wirth, 2010). It is well known that the thermodynamic and kinetic properties of defects on structural defects such as interface and grain boundaries are very different from those inside grains, which may result in inhomogeneous RIS and inhomogeneous precipitation of second phase. But the inhomogeneous thermodynamic and kinetics properties haven’t been considered in the existing rate theory models. Phase-field approach is good at describing spatial and temporal evolution of chemistry and microstructure in a system with inhomogeneous thermodynamic and kinetic properties. Current phase-field models of RIS (Piochaud et al., 2016; Xia et al., 2020) in bi-crystalline structure were extended from the rate theory. Like the rate theory, the sink strength of grain boundaries is derived by assuming that defects at sinks remain thermal equilibrium concentrations. However, this assumption might not be true for a system far from equilibrium. In nuclear materials, especially in materials with high irradiation dose and/or high dose rate, rich microstructure change such as recrystallization and gas bubble evolution can lead the system far from equilibrium.

In this work, we extended the rate theory to considers the effect of the inhomogeneous thermodynamic and kinetic properties of defects on RIS in polycrystalline structures. Defect diffusion is assumed to be driven by chemical potential gradient. The sink strength at grain boundaries that depends on local defect concentrations, diffusivity, and reaction rate, is replaced by local reaction. The model can release the assumptions discussed above. Because the same thermodynamic and kinetic models are used in phase-field models of phase transition the current model can be coupled with phase-field model to study radiation-induced precipitation.

Microstructure-Dependent Rate Theory Model

Description of Chemistry and Defects

We consider an AB binary model alloy and regard atom A as solute atom while atom B as solvent atom. In the conventional rate theory (Was, 2016), three independent variables, i.e., concentrations of atom A, interstitial (A and B) and vacancy (A and B) are used to describe spatial and temporal distributions of solute A and defects. Considering the fact that interstitials A and B might have 1) different generation rates during the cascade; 2) different chemical potentials on sinks, 3) different recombination rates with vacancies; 4) different binding energy with atom A and/or B, four independent variables, i.e., concentrations of atom A, interstitial A, interstitial B, and vacancy are used in the current model. They are denoted by , , , and , respectively. is the spatial coordinate and is time. It should be pointed out that is the total concentration including atom A in the lattice and the interstitial A which is also used in the rate equations (Was, 2016). Radiation generates the Frenkel pair such as interstitial A and vacancy A. When solute A has a low concentration, vacancy A has a low concentration as well. However, vacancies at atom A and atom B sites should have the similar local atomic configuration so that vacancy A and vacancy B have the same thermodynamic and kinetic properties. Therefore, it should be reasonable to use one variable describing the total vacancy concentration. Our model uses one sublattice i.e., sublattice [A, B, vacancy, Ai, or Bi] to describe the chemistry and defects, and uses one set of parameters to describe the grain orientations in the polycrystalline structure where is the total number of grain orientations in the simulation cell. The order parameter , which is obtained from phase-field modeling of grain growth, has the values of 1 inside the grain m, and 0 outside of the grain m, and continuously varies from 1 to 0 across the grain boundary. We define a shape function with the order parameters as that varies smoothly from 0 inside the grains to 1.0 at the centers of GBs. With the sets of concentration and order parameter variables, the inhomogeneous thermodynamic and kinetic properties of defects can be described.

Chemical Potentials

For a regular solution, the Gibbs free energy can be written aswhere and are the concentrations of atom A and B, are the Gibbs free energies of pure A and B, is the activity coefficient of species j, is the Boltzmann constant, and is the absolute temperature.

In irradiated materials, the vacancy and interstitial concentrations can be a few orders of magnitude higher than their thermal equilibrium concentrations although their concentrations are still lower compared with the concentrations of alloy atoms. The high defect concentrations affect not only the diffusivities of species, but also chemical potentials of species. If one assumes that 1) A and B form a regular solution; and 2) vacancy and interstitial A and B are dilute defects, the Gibbs free energy of the system can be written asThe concentration of solvent atom B is a dependent variable and can be calculated by if defect concentrations are low (i.e., , and <<1).

Chemical potentials of species at structural defects such as GBs, interfaces, and surfaces are usually different from those inside grains. The microstructure and/or spatial dependent chemical potential in the polycrystalline AB alloy cause an additional diffusion driving force that affects the fluxes of defects and species, hence, their segregation or depletion, and phase stability. To capture the effect of thermodynamic and kinetic properties of GBs on RIS, we extend the conventional rate theory to a microstructure-dependent rate theory model. The thermodynamic properties of the GBs are expected to closely correlate with the atomic density and composition. van der Waals (1979) showed that the energy of an interface can be described as a function of mass density and its variations within the interface region. Very recently, Kamachali et al. (2019) proposed a density-based thermodynamic model of GBs. The Gibbs free energy of GBs was described by the bulk free energy, atomic density, and gradient coefficient of atomic density. A Gaussian function was used to describe the atomic density change across a GB that evolves with composition. Atomic simulations (Kamachali et al., 2019) in alpha Fe showed that the atomic density of a GB with high symmetry has the Gaussian distribution. We use a shape function of to express the atomic density of a GB, which is equal to inside the grain and smoothly varies to at the center of the GB. is the density of a perfect crystal while is the difference of atomic density between those inside a grain and on its grain boundary center due to free volumes at the center of the grain boundary. Similarly, other spatially dependent thermodynamic and kinetics properties such as chemical potential and diffusivity of species j can also be described as , . () is the property of species inside a grain, and is the difference of the property of species on the GB from those inside the grain. denotes any of the defects or solute A i.e., . Therefore, the microstructure-dependent chemical potential of species j can be written aswhere is the bulk chemical potential of species j, is the difference of chemical potential of species j inside a grain and its grain boundary center. The diffusivity inhomogeneity can also be described by a similar expression as Eq. 3where is the diffusivity of species j inside grains (also called bulk diffusivity) and is the difference of diffusivity of species j between the bulk diffusivity and the diffusivity at the center of grain boundaries.

The spatial dependent chemical potential of species j on grain boundaries can be expressed as , where and are the difference of formation enthalpy and entropy due to the presence of grain boundaries, respectively. Low-angle symmetric tilt grain boundaries can be modeled by a wall of parallel edge dislocations (Hirth et al., 1968; Jiang et al., 2014; Xia et al., 2020). The elastic interaction between stress field associated with the dislocation array and point defects affect the formation (enthalpy) energy of species j near the grain boundaries. (Jiang et al., 2014) analyzed the effect of grain misorientation angles on the elastic interaction energy. Their results show that when the misorientation angle becomes greater than the stress fields of dislocations overlap and cancel each other so that the elastic interaction is negligible. Therefore, can also be written as , where and are the effect of grain boundary atomic structures and elastic interaction on the chemical potential of species j, respectively. For a given grain boundary can be calculated by atomistic simulations (Tschopp et al., 2011). can be assessed by analyzing the effect of the stress associated with the dislocation array on the formation enthalpy and entropy of species j. For symmetric tilt grain boundaries can be expressed as Eqs 3, 4 by ignoring the effect of dislocation stresses on the formation entropy (Jiang et al., 2014). is not zero within grain boundaries about few lattice constant thickness (Tschopp et al., 2011). However, reflects the elastic interaction which is a long-range interaction. For small angle grain boundaries the interaction region may reach tens of nanometers (Jiang et al., 2014). In this work, we focused on studying the effect of inhomogeneous chemical potentials on RIS in model alloys. Chemical potentials was constructed to reasonably capture the value and the length scale of elastic interaction by adjusting the shape function and the characteristic length . Similarly, could be constructed by the diffusivity data near grain boundaries from atomistic simulations and experiments.

Microstructure-Dependent Rate Theory

The rate theory model of RIS (Was, 2016) has been successfully used in predicting the effects of defect mobility, solute size, and defect-binding energies on solute segregation on surfaces and grain boundaries. However, the effects of grain morphology and inhomogeneous thermodynamic and kinetic properties of defects on RIS were ignored because most modeling was carried out in one dimension. In this work, we extend the rate theory model of RIS into polycrystalline structures. According to the rate theory, the evolution of chemical and defect concentrations is given by:where and are the fluxes of atom and defect , respectively. , , and is the generation rate, sink rate, and recombination rate of defect , respectively. The conservation principle requires the fluxes should satisfy the following equations (Piochaud et al., 2016):where for and for . In the framework of the thermodynamics of irreversible process (Onsager, 1931; Piochaud et al., 2016), the fluxes of chemical and defects driven by chemical potential gradient can be approximately calculated by,where is the Onsager kinetic coefficient. If we define the normalized Onsager coefficient asthe evolutions Eqs 5, 6 can be expressed asIn polycrystalline structures the inhomogeneous chemical potential and are expressed asAssuming that 1) the chemical potentials , satisfy the Gibbs-Duhem relationship ; 2) the defect concentration is low i.e., , then the chemical potential gradient can be calculated byWith the assumption of for the evolution Eqs 11, 12 x are simplified aswhere are functions of defect concentrations and chemical potentials on grain boundaries. They are given byIf we further assume that , the evolution Eqs 17–20 inside the grains (i.e., ) become to the rate Eqs 6–20 of ref. (Was, 2016). The normalized Onsager coefficients are the same as in Eqs 6–20 of ref. (Was, 2016) which is the diffusivity of atom diffusing via exchanging with defect on a given neighboring site. The diffusivity coefficient is given bywhere is the jump distance when atom and defect j exchange sites. is the correlation factor (Wharry and Was, 2014). is the effective exchange-jump frequency of atom - defect pair, and is the migration energy. In this work, and .

In the conventional rate theory, the generation rate of defects is assumed to be proportional to the alloy composition. Since interstitials A and B are treated as individual variables, more accurate interstitial generation rates from MD cascade simulations can be used when they satisfy the relationship of . and are calculated by the equation of where is the dose rate and is the production efficiency of defect . The recombination rate of Frenkel defects . And . is a rate constant, is the lattice constant, and is the combination factor of vacancies and interstitials. In polycrystalline structures, the main sinks of defects are dislocations and grain boundaries. The sink rate of defects on distributed dislocations can be calculated by , where is the sink strength constant of defect on dislocations. is the dislocation density. is the diffusivity of defect . is the thermal equilibrium concentration of defect at dislocations. Grain boundaries can be described by an array of dislocations such as small angle grain boundaries (Duh et al., 2001) and/or can be described as plane defects with high defect concentrations as discussed in Section 2.2. Defects (Ai, Bi, V) on grain boundaries may have different chemical potentials from those inside grains. Distributed dislocations on grain boundaries have elastic interaction with defects. Both the chemical potential and elastic interaction cause driving forces for defect diffusion. Conventional rate theory assumes that defects on grain boundaries remain their thermal equilibrium concentrations. It should be a reasonable assumption after the system reaches a steady state. However, for a system far from equilibrium, the assumption might be inappropriate. For example, if solute interstitials have much higher diffusivity than that of vacancies, the solute interstitials may accumulate on grain boundaries and form supper saturated solution. The excess interstitials might be emitted back to interior grains. But it is also possible that the super saturated solution or vacancies leads to the second phase formation or void formation on grain boundaries which are often observed in irradiated materials (Kuksenko et al., 2012; Zhao et al., 2018). In this work, the polycrystalline structure is naturally described in terms of inhomogeneous thermodynamic and kinetic properties. The sink of defects on grain boundaries is described by the local defect recombination enhanced by high defect concentrations and high defect mobility. Therefore, the assumption that defects on grain boundaries remain their thermal equilibrium concentrations in the conventional rate theory is released. A sink term on grain boundaries is added in to force defects having equilibrium concentrations on grain boundaries and study the effect of this sink term on the solute segregation.

Model Parameters

The normalized time , diffusivity , gradient , and energy are used for numerically solving Eqs 17–20. and are, respectively, the characteristic length and the largest partial diffusion coefficient of . Solving the equations, one has the spatial and temporal evolution of species concentration . To simulate the effect of inhomogeneous thermodynamic and kinetic properties on RIS, parametric studies are carried out in model alloys. Table 1 lists the model parameters and thermodynamic and kinetic properties of defects used in the simulations.

Results

We consider a model 10at % AB alloy. A denotes the solute atom while B the solvent atom. A phase field model of grain growth was used to generate a polycrystalline structure. The grain boundaries are defined by a shape function that has the value of 0 inside the grains and continuously varies to 1 at the center of the grain boundaries. Figure1A shows a two dimensional (2D) polycrystalline structure which has a dimension and is used in the simulations. The average grain size is about . From the 2D grain structure, we can see the microstructure features such as different sized grains, plate grain boundaries and triple points of grain boundaries. These features are like those in 3D grain structures. It is anticipated that the 2D simulations captured the physics of the effect of microstructures on RIS. Figure1B shows the distributions of the shape function along the lines and respectively. The inhomogeneous thermodynamic and kinetic properties of species and defects in the polycrystalline structure are defined by as described in Section 2.2.

FIGURE 1

Effect of Thermodynamic and Kinetic Properties of Grain Boundaries on Solute Segregation

In the conventional rate theory, the grain boundary is treated as a perfect sink of defects by assuming that defects have their thermal equilibrium concentrations. Current model releases this assumption and uses inhomogeneous chemical potentials () and diffusivity () of defects to describe the features of grain boundaries as sinks of defects. First, we consider the effect of solute atom’s chemical potential on RIS. In the simulations, the following model parameters are set up to be , , , and thus . The other model parameters are listed in Table 1. takes different values of 0.4, 0.0–0.4, and –0.6°[eV/atom], respectively, which means that the chemical potential of the solute A and self-interstitial on the grain boundaries is higher (positive value) or lower (negative value) than that inside grains by the value of . Figure 2 shows the distribution of the solute concentration change at time . where is the initial concentration of atom A and uniform in the simulation cell. The same color bar is used in the figures so that the color presents not only the effect of chemical potentials on the solute segregation or depletion but also the relative strength of RIS in the polycrystalline structure. To make more clear comparison, the distributions of the solute concentration change along the line are plotted in Figure 3, where the dashed lines mark the locations of grain boundaries.

FIGURE 2

FIGURE 3

From Figures 2, 3, it can be clearly seen that 1) the solute A segregates on grain boundaries when i.e., solute A has the same chemical potential on grain boundaries and inside grains. The solute segregation results from the inhomogeneous mobility of radiation defects. The larger diffusivity of defects on grain boundaries causes a faster defect recombination, hence, bigger defect concentration gradient, defect flux and RIS; 2) when solute A on grain boundaries has a higher chemical potential () than that inside grains the solute depletes on the grain boundaries. In this case, the chemical potential gradient of solute A dominates its flux; and 3) solute segregation increases with the decrease of the chemical potential (). In this case, both chemical potential gradient and defect concentration gradient drive the solute diffusing to grain boundaries. The distributions of solute concentration on grain boundaries also show that the solute segregation or depletion on grain boundaries are not uniform. Solute concentration at triple points of grain boundaries is much higher or lower than that on straight grain boundaries. In addition, solute segregation on short grain boundaries is more uniform than that on large grain boundaries. The results imply more inhomogeneous RIS in polycrystalline structure with larger grains than that in polycrystalline structures with smaller grains. It should point out that a number of factors may cause the inhomogeneous RIS on grain boundaries. For example, from Figure 1B we can find that the shape function which is defined by the order parameters obtained from phase-field modeling is not exactly equal to 1 on the grain boundaries. Consequently, the thermodynamic and kinetic properties of defects on grain boundaries are not the same. Actually, such a shape function naturally describes the inhomogeneous structures of grain boundaries at triple points i.e., curved boundaries connected different oriented grains. The inhomogeneous RIS on grain boundaries, which are often observed in experiments, confirms that the shape function can naturally describe the thermodynamic and kinetic properties of grain boundaries.

The equilibrium concentrations of interstitials () and vacancy () on grain boundaries are setup to be , respectively. The distributions of defect concentrations along the line are plotted in Figure 4. It can be seen that vacancies deplete on grain boundaries because the given defect diffusivity on grain boundaries is two order magnitude higher () than that inside grains. Therefore, a fast defect recombination reduces the vacancy concentration. By the natural recombination, concentrations of vacancies closely reach to its thermal equilibrium concentration. However, the concentration of interstitial A increases with the decrease of chemical potential because the flux of interstitial A toward grain boundaries increases. When is positive, the interstitial A depletes at grain boundaries. A negative leads to a segregation of interstitial A at grain boundaries as shown in Figure 4B. This is because interstitial A has higher solubility (or equilibrium concentration) at grain boundaries with a negative . Concentration of interstitial B is much higher than its thermal equilibrium concentration and increases with the decrease of chemical potential . Besides the defect recombination on grain boundaries, other material processes may affect the defect concentrations as well. For example, interstitials can be absorbed by grain growth, and interstitials may diffuse to free surface through grain boundaries. Both processes reduce the interstitial concentration on grain boundaries. Current model uses the sink term “” to describe the absorption mechanisms of interstitials on grain boundaries. We can turn on or off this term to study the effect of sink on RIS which releases the assumption of perfect sink in the conventional rate theory. Increasing the means more defects sinking on the grain boundaries. Figure 5A shows the solute concentration changes along at time for . The defect concentrations are plotted in Figures 5B,C. With the sink condition the results indicate that 1) the concentrations of all defects reach their thermal equilibrium concentrations on grain boundaries 2) all defect concentrations inside grains decreases, and 3) the solute segregation on grain boundaries reduces. The reduction of RIS with the increase of sink strength () is because the absorption of defects on grain boundaries lowers defect concentrations, hence the RIS kinetics. The results indicate that the assumption that all defects remain their thermal equilibrium concentrations in the rate theory is not necessarily correct unless the extra defects can be absorbed or emitted by grain boundaries. Grain boundaries acting as perfect sinks ( reduces defect concentrations inside grain and the solute segregation kinetics on grain boundaries.

FIGURE 4

FIGURE 5

Effect of Binding Energy on Solute Segregation

For a regular solution AB, the mixing entropy can be approximated by the ideal solution mixing entropy:The non-ideality of the solution is represented by the enthalpy of mixing, which, in the quasi-chemical approximation, is given bywhere is defined as the binding energy of species A and B, and is given bywhere is the number of nearest neighbors, is the bond enthalpy of the i − j bond. The free energy of mixing can be calculated by:The Gibbs free energy for a regular solution can be written aswhere and are the Gibbs free energies of pure A and B, respectively. The first three terms in Eq. 28 describe the free energy of ideal solution, and the last term describes the mixing enthalpy of the regular solution. Equation 28 is an alternative expression of the Gibbs free energy described by Eq. 1. The thermodynamic factor in Eq. 15 can be calculated by i.e., . With the same model parameters used in Section 4.1, but and , the effect of bonding energy on RIS is simulated and presented in Figure 6. From Eq. 15, a negative AB binding energy increases the thermodynamic factor . In contrast, a positive AB binding energy decreases the thermodynamic factor . Therefore, AB binding energy affects the fluxes driven by the concentration gradients of atom A and B. As a result, it affects the solute segregation. The results in Figure 6 demonstrate a clear tendency that a weaker AB binding energy () increases the RIS while a stronger AB binding energy () decreases the solute segregation.

FIGURE 6

If the binding among defects is strong, the faster diffusive defect has a drag effect on the slower diffusive defect. For instance, undersize solutes may tightly bind to interstitials forming interstitial-solute complexes that migrate as solute interstitials. In the rate theory (Was, 2016), the binding energy is used to scale the interstitial concentration of the solute which affect the partial diffusivity. The Gibbs free energy of the system with strong defect interaction can be generally written asThe binding of defects and results in an additional flux on the flux of defect associated with the concentration gradient of defect . The current model can consider the effect of defect binding on RIS once the binding energy is available.

Effect of Defect Generation Efficiency on Solute Segregation

In conventional rate theory, the generation rate of solute interstitials is assumed to be proportional to its concentration in the alloy i.e., where is the dpa rate. However, this generation rate may be only correct in the ideal alloy. MD simulations of cascades (Zhang et al., 2017) show that in Fe-Cr, the Cr interstitial fraction is much higher than the Cr solute concentration and the Cr interstitial production efficiency decreases with the increasing Cr concentration. In 10 at%CrFe, the fraction of the Cr interstitial is about 65–78% while the fraction of Fe interstitial is only about 35% or even lower. In contrast, in Fe-Cu, Cu interstitials are barely produced. In current model, the concentrations of interstitial A and B are treated as independent variables so that it can describe the effect of defect generation efficiency on RIS. With the same model parameters in Section 4.1, but and , and, the effect of production efficiency on RIS is simulated and displayed in Figures 7, 8. In the case of , the solute segregate on grain boundaries as shown by the green line in Figure 8. This indicate that solutes have a flux toward grain boundaries. The results that increasing the generation rate of solute interstitials increase the RIS is expected. Therefore, in CrFe and CuFe alloys, it is important to consider the effect of solute interstitial production efficiency on RIS because their production efficiency is largely different from the concentration of solute.

FIGURE 7

FIGURE 8

Effect of Defect Diffusivity on Solute Segregation

In the rate theory of RIS, the inhomogeneity of chemical potential in polycrystalline structures is ignored i.e., . At steady state the rate theory predicts that solute and vacancy concentration gradients have the relationship: (Was, 2016). With the migration energies and temperature listed in Table 1, we have the diffusivities of , respectively. Thus, we have . The negative value of means . The solute segregation on grain boundaries shown by the red line in Figure 8 and the vacancy depletion on grain boundaries shown in the red line in Figure 10 indicate . Two other model systems with different values of are considered to validate the model as well as the prediction of rate theory. In Model 1, the partial diffusivities () are setup, which give . In Model 3, the partial diffusivities () are setup, which give . The results in Figure 9 shows that effect of the relative partial diffusivity on RIS. For Model 2, is equal to be . It is clear that vacancy always deplete on grain boundaries shown in Figure 9B. Vacancy concentration inside grains increases with the decrease of . Figure 9A demonstrates that the relationship hold which is in agreement with the rate theory’s prediction.

FIGURE 9

FIGURE 10

Conclusion and Discussion

In this work the conventional rate theory of RIS has been extended by taking into account inhomogeneous thermodynamic and kinetics properties of defects. The developed model has the following features: 1) concentrations of solute and solvent interstitials are treated as independent variables. In conventional rate theory the total concentration of solute and solvent interstitials is treated as an independent variable and the fraction of the solute interstitial in the total interstitial concentration is assumed to be the same as the solute concentration in the alloy. The independent solute and solvent interstitial concentrations enable one to describe different thermodynamic and kinetic properties of solute and solvent interstitials including the chemical potentials and diffusivities on grain boundaries, production efficiency and binding energy among defects. 2) thermodynamic and kinetic properties of defects in a polycrystalline structure are expressed to be spatial or microstructure dependent, which enables one to examine the effect of thermodynamic and kinetic properties of defects on RIS at grain boundaries. 3) the assumption that grain boundaries are perfect sink of defects is released. Defect concentrations on grain boundaries are determined by their chemical potentials, recombination and absorption rates. The defect absorption refers to the scenario that defects are consumed by grain growth or migrate to free surface along grain boundaries. With the model, the effect of structural and concentration dependent thermodynamic and kinetic properties on RIS were simulated. For the given kinetic properties of defects inside grains and on grain boundaries, the results indicate that 1) vacancy depletion on grain boundaries is always observed and the relative kinetic property of determines the solute flux i.e., and segregation or depletion on grain boundaries, which is in agreement with the steady state solution predicted by the conventional rate theory; 2) RIS strongly depends on solute chemical potential on grain boundaries. Increasing the solute chemical potential on grain boundaries may change solute segregation to solute depletion, or vice versa; 3) solute segregation or depletion on grain boundaries is non-uniform. Solutes have strong segregation or depletion at triple points of grain boundaries. It implies that second phase may form first at the triple points which are, actually, often observed in materials. 4) the assumption that all defects remain their thermal equilibrium concentrations in the rate theory is not necessarily correct unless the extra defects can be absorbed or emitted by grain boundaries. Grain boundaries acting as perfect sinks reduces defect concentrations inside grain and the solute segregation kinetics on grain boundaries. 5) the production efficiency of solute interstitials and the binding energy among defects are important parameters which affect the RIS.

The simulations demonstrate that developed model extends and strengthens the capability of the conventional rate theory of RIS. However, the inhomogeneous thermodynamic and kinetic properties add challenges in simulations. For example, to capture the thickness of grain boundaries about few nanometers we have to use a small grid size (Piochaud et al., 2016; Xia et al., 2020). As we know that the time step for solving the diffusion equations can be estimated by , where is the largest diffusivity of defects of interest. Decreasing and increasing reduce the time step and increase the computational cost. This is the reason why we presented the RIS at very early stage , and used a high dpa rate which can increase the defect concentration and speed up the RIS in this work. We can run a one dimensional (1D) model with larger grains, lower dpa rate and longer radiation time. Figure 10 shows the 1D result. It confirms that decreasing the dpa rate slows down the RIS kinetics, but doesn’t affect the RIS feature (depletion or segregation) for given thermodynamic and kinetics properties of defects. In irradiated materials, interstitials, which are dominant defects, usually have large and strong anisotropic diffusivity (such as one dimensional diffusion). The large and anisotropic diffusivity adds more difficulty in solving the diffusion equations. The model can be extended to use the first-passage approach (Hu and Henager, 2009; Hu et al., 2016) to describe the large and anisotropic diffusivity and to increase the time step and computational efficiency.

The characteristic length is an important model parameter. A larger allows to use a larger time step in the simulations. Two ways might be used to increase the . One is to construct a shape function which could describe an arrow or sharp grain boundaries. The other is to use an adaptive mesh to capture the thermodynamic and kinetic properties of defect on structural defects (Tonks et al., 2012; Chakraborty et al., 2016). Structural defects such as grain boundaries, second phase particles in ODS, and interfaces in coated cladding materials act as sinks of defects. The local chemistry change affects the phase stability, hence, microstructure change and material property degradation. The developed model is based on fundamental properties of defects such as the properties of grain boundaries (Admal et al., 2018) and interface, partial diffusivity (Messina et al., 2014) and binding energies. With the thermodynamic and kinetics properties of materials of interest, the model can be extended and applied in studying the effect of microstructure dependent thermodynamic and kinetic properties on defect evolution and RIS in materials with complicated microstructures.

References

• 1

  AdmalN. C.PoG.MarianJ. (2018). A Unified Framework for Polycrystal Plasticity with Grain Boundary Evolution. Int. J. Plasticity106, 1–30. 10.1016/j.ijplas.2018.01.014

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AkamatsuM.Van DuysenJ. C.PareigeP.AugerP. (1995). Experimental Evidence of Several Contributions to the Radiation Damage in Ferritic Alloys. J. Nucl. Mater.225, 192–195. 10.1016/0022-3115(95)00028-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  AllenT. R.WasG. S. (1998). Modeling Radiation-Induced Segregation in Austenitic Fe-Cr-Ni Alloys. Acta Materialia46 (10), 3679–3691. 10.1016/s1359-6454(98)00019-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BadilloA.BellonP.AverbackR. S. (2015). A Phase Field Model for Segregation and Precipitation Induced by Irradiation in Alloys. Model. Simulation Mater. Sci. Eng.23 (3). 10.1088/0965-0393/23/3/035008

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BruemmerS. M.SimonenE. P.ScottP. M.AndresenP. L.WasG. S.NelsonJ. L. (1999). Radiation-induced Material Changes and Susceptibility to Intergranular Failure of Light-Water-Reactor Core Internals. J. Nucl. Mater.274 (3), 299–314. 10.1016/s0022-3115(99)00075-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  ChakrabortyP.ZhangY.TonksM. R. (2016). Multi-scale Modeling of Microstructure Dependent Intergranular Brittle Fracture Using a Quantitative Phase-Field Based Method. Comput. Mater. Sci.113, 38–52. 10.1016/j.commatsci.2015.11.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  DamcottD. L.AllenT. R.WasG. S. (1995). Dependence of Radiation-Induced Segregation on Dose, Temperature and Alloy Composition in Austenitic Alloys. J. Nucl. Mater.225, 97–107. 10.1016/0022-3115(94)00690-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  DuhT. S.KaiJ. J.ChenF. R.WangL. H. (2001). Numerical Simulation Modeling on the Effects of Grain Boundary Misorientation on Radiation-Induced Solute Segregation in 304 Austenitic Stainless Steels. J. Nucl. Mater.294 (3), 267–273. 10.1016/s0022-3115(01)00493-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  DunnA. Y.CapolungoL.MartinezE.CherkaouiM. (2013). Spatially Resolved Stochastic Cluster Dynamics for Radiation Damage Evolution in Nanostructured Metals. J. Nucl. Mater.443 (1-3), 128–139. 10.1016/j.jnucmat.2013.07.009

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  FaneyT.WirthB. D. (2014). Spatially Dependent Cluster Dynamics Modeling of Microstructure Evolution in Low Energy Helium Irradiated Tungsten. Model. Simulation Mater. Sci. Eng.22 (6). 10.1088/0965-0393/22/6/065010

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  FarrellK.MahmoodS. T.StollerR. E.MansurL. K. (1994). An Evaluation of Low Temperature Radiation Embrittlement Mechanisms in Ferritic Alloys. J. Nucl. Mater.210 (3), 268–281. 10.1016/0022-3115(94)90181-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  FieldK. G.BarnardL. M.ParishC. M.BusbyJ. T.MorganD.AllenT. R. (2013). Dependence on Grain Boundary Structure of Radiation Induced Segregation in a 9 wt.% Cr Model Ferritic/martensitic Steel. J. Nucl. Mater.435 (1-3), 172–180. 10.1016/j.jnucmat.2012.12.026

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  GolubovS. I.OvcharenkoA. M.BarashevA. V.SinghB. N. (2001). Grouping Method for the Approximate Solution of a Kinetic Equation Describing the Evolution of Point-Defect Clusters. Philos. Mag. A81 (3), 643–658. 10.1080/01418610108212164

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  GrandjeanY.BellonP.MartinG. (1994). Kinetic Model for Equilibrium and Nonequilibrium Segregation in Concentrated Alloys under Irradiation. Phys. Rev. B50 (6), 4228–4231. 10.1103/physrevb.50.4228

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  HuR.SmithG. D. W.MarquisE. A. (2012). Atom Probe Study of Radiation Induced Grain Boundary Segregation/depletion in a Fe-12%Cr Alloy. Prog. Nucl. Energ57, 14–19. 10.1016/j.pnucene.2011.10.011

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  HuS.BurkesD. E.LavenderC. A.SenorD. J.SetyawanW.XuZ. (2016). Formation Mechanism of Gas Bubble Superlattice in UMo Metal Fuels: Phase-Field Modeling Investigation. J. Nucl. Mater.479, 202–215. 10.1016/j.jnucmat.2016.07.012

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  HuS.HenagerC. H. (2009). Phase-field Modeling of Void Lattice Formation under Irradiation. J. Nucl. Mater.394 (2), 155–159. 10.1016/j.jnucmat.2009.09.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  HuangC. H.MarianJ. (2016). A Generalized Ising Model for Studying Alloy Evolution under Irradiation and its Use in Kinetic Monte Carlo Simulations. J. Phys-condens Mat28 (42). 10.1088/0953-8984/28/42/425201

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  JiangC.SwaminathanN.DengJ.MorganD.SzlufarskaI. (2014). Effect of Grain Boundary Stresses on Sink Strength. Mater. Res. Lett.2 (2), 100–106. 10.1080/21663831.2013.871588

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  JourdanT.BencteuxG.AdjanorG. (2014). Efficient Simulation of Kinetics of Radiation Induced Defects: A Cluster Dynamics Approach. J. Nucl. Mater.444 (1-3), 298–313. 10.1016/j.jnucmat.2013.10.009

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  HirthJ. P.LotheJ.AndersonP. M. (1968). Theory of Dislocation. New York: McGraw-Hill Book Company.

  • Google Scholar
  • View reference in article

• 22

  KamachaliR. D.da SilvaA. K.McEniryE.PongeD.GaultB.NeugebauerJ.et al (2019). Segregation–Assisted Spinodal and Transient Spinodal Phase Separation at Grain Boundaries. npj Comput Mater6, 191 (2020). https://doi.org/10.1038/s41524-020-00456-7

  • Google Scholar
  • View reference in article

• 23

  KeJ.-H.KeH.OdetteG. R.MorganD. (2018). Cluster Dynamics Modeling of Mn-Ni-Si Precipitates in Ferritic-Martensitic Steel under Irradiation. J. Nucl. Mater.498, 83–88. 10.1016/j.jnucmat.2017.10.008

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  KuksenkoV.PareigeC.PareigeP. (2012). Intra Granular Precipitation and Grain Boundary Segregation under Neutron Irradiation in a Low Purity Fe-Cr Based Alloy. J. Nucl. Mater.425 (1), 125–129. 10.1016/j.jnucmat.2011.10.031

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  LamN. Q.OkamotoP. R.WiedersichH.TaylorA. (1978). Radiation-Induced Solute Segregation and Precipitation in Alloys. Mta9 (12), 1707–1714. 10.1007/bf02663400

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  LiB. Y.HuS. Y.LiC. L.LiQ. L.ChenJ.ShuG. G.et al (2017). Simulations of Irradiated-Enhanced Segregation and Phase Separation in Fe-Cu-Mn Alloys. Model. Simulation Mater. Sci. Eng.25 (6). 10.1088/1361-651x/aa7197

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  LuZ.FaulknerR. G.WasG.WirthB. D. (2008). Irradiation-induced Grain Boundary Chromium Microchemistry in High Alloy Ferritic Steels. Scripta Materialia58 (10), 878–881. 10.1016/j.scriptamat.2008.01.004

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  MessinaL.NastarM.GarnierT.DomainC.OlssonP. (2014). Exact Ab Initio Transport Coefficients in Bcc Fe-X (X=Cr, Cu, Mn, Ni, P, Si) Dilute Alloys. Phys. Rev. B90 (10). 10.1103/physrevb.90.104203

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  NastarM.SoissonF. (2012). “Radiation-Induced Segregation,” in Comprehensive Nuclear Materials. Editor KoningsR. J. M. (Oxford: Elsevier), 471–496. 10.1016/b978-0-08-056033-5.00035-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  OdetteG. R.LucasG. E. (1986). Irradiation Embrittlement of Reactor Pressure Vessel Steels: Mechanisms, Models, and Data Correlations. Editor SteeleL. E. (West Conshohocken, PA: ASTM International), 206–241.

  • Google Scholar
  • View reference in article

• 31

  OnsagerL. (1931). Reciprocal Relations in Irreversible Processes. I.Phys. Rev.37 (4), 405–426. 10.1103/physrev.37.405

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  OrtizM. J. C. C. J. (2007). Simulation of Defect Evolution in Irradiated Materials: Role of Intracascade Clustering and Correlated RecombinationPhys. Rev. B75, 184101. 10.1103/PhysRevB.75.184101

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  PiochaudJ. B.NastarM.SoissonF.ThuinetL.LegrisA. (2016). Atomic-based Phase-Field Method for the Modeling of Radiation Induced Segregation in Fe-Cr. Comput. Mater. Sci.122, 249–262. 10.1016/j.commatsci.2016.05.021

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  SenningerO.SoissonF.MartínezE.NastarM.FuC.-C.BréchetY. (2016). Modeling Radiation Induced Segregation in Iron-Chromium Alloys. Acta Materialia103, 1–11. 10.1016/j.actamat.2015.09.058

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  SimonenE. P.BruemmerS. M. (1999). Radiation-induced Grain Boundary Segregation in Austenitic Stainless Steels. Mater. Sci. Forum294-2, 755–758. 10.4028/www.scientific.net/MSF.294-296.755

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  SoissonF.JourdanT. (2016). Radiation-accelerated Precipitation in Fe-Cr Alloys. Acta Materialia103, 870–881. 10.1016/j.actamat.2015.11.001

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  SouidiC. S. B. A.DomainC.TerentyevD.MalerbaL.CalderA. F.BaconD. J.et al (2006). Dependence of Radiation Damage Accumulation in Iron on Underlying Models of Displacement Cascades and Subsequent Defect Migration. J. Nucl. Mater.355 (1–3), 89–103. 10.1016/j.jnucmat.2006.04.009

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  TonksM. R.GastonD.MillettP. C.AndrsD.TalbotP. (2012). An Object-Oriented Finite Element Framework for Multiphysics Phase Field Simulations. Comput. Mater. Sci.51 (1), 20–29. 10.1016/j.commatsci.2011.07.028

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  TschoppM. A.HorstemeyerM. F.GaoF.SunX.KhaleelM. (2011). Energetic Driving Force for Preferential Binding of Self-Interstitial Atoms to Fe Grain Boundaries over Vacancies. Scripta Materialia64 (9), 908–911. 10.1016/j.scriptamat.2011.01.031

  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  van der WaalsJ. D. (1979). The Thermodynamic Theory of Capillarity under the Hypothesis of a Continuous Variation of Density. J. Stat. Phys.20 (2), 200–244. 10.1007/bf01011514

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  WasG. S. (2016). Fundamentals of Radiation Materials Science: Metals and Alloys. Springer.

  • Google Scholar
  • View reference in article

• 42

  WasG. S.WharryJ. P.FrisbieB.WirthB. D.MorganD.TuckerJ. D.et al (2011). Assessment of Radiation-Induced Segregation Mechanisms in Austenitic and Ferritic-Martensitic Alloys. J. Nucl. Mater.411 (1-3), 41–50. 10.1016/j.jnucmat.2011.01.031

  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  WharryJ. P.WasG. S. (2013). A Systematic Study of Radiation-Induced Segregation in Ferritic-Martensitic Alloys. J. Nucl. Mater.442 (1-3), 7–16. 10.1016/j.jnucmat.2013.07.071

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  WharryJ. P.WasG. S. (2014). The Mechanism of Radiation-Induced Segregation in Ferritic-Martensitic Alloys. Acta Materialia65, 42–55. 10.1016/j.actamat.2013.09.049

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  WooC. H.SinghB. N. (1992). Production Bias Due to Clustering of Point Defects in Irradiation-Induced Cascades. Philos. Mag. A65 (4), 889–912. 10.1080/01418619208205596

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  XiaL. D.JiY. Z.LiuW. B.ChenH.YangZ. G.ZhangC.et al (2020). Radiation Induced Grain Boundary Segregation in Ferritic/martensitic Steels. Nucl. Eng. Technol.52 (1), 148–154. 10.1016/j.net.2019.07.009

  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  XuD. H.WirthB. D. (2010). Modeling Spatially Dependent Kinetics of Helium Desorption in BCC Iron Following He Ion Implantation. J. Nucl. Mater.403 (1-3), 184–190. 10.1016/j.jnucmat.2010.06.025

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  XuD.WirthB. D. (2009). Spatially Dependent Rate Theory Modeling of Thermal Desorption Spectrometry of Helium-Implanted Iron. Fusion Sci. Technol.56 (2), 1064–1068. 10.13182/fst09-a9052

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  YangY.FieldK. G.AllenT. R.BusbyJ. T. (2016). Roles of Vacancy/interstitial Diffusion and Segregation in the Microchemistry at Grain Boundaries of Irradiated Fe-Cr-Ni Alloys. J. Nucl. Mater.473, 35–53. 10.1016/j.jnucmat.2016.02.007

  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  ZhangY. X.SchwenD.BaiX. M. (2017). Molecular Dynamics Simulations of Concentration- Dependent Defect Production in Fe-Cr and Fe-Cu Alloys. J. Appl. Phys.122 (22). 10.1063/1.5008757

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  ZhaoH.De GeuserF.Kwiatkowski da SilvaA.SzczepaniakA.GaultB.PongeD.et al (2018). Segregation Assisted Grain Boundary Precipitation in a Model Al-Zn-Mg-Cu Alloy. Acta Materialia156, 318–329. 10.1016/j.actamat.2018.07.003

  • CrossRef
  • Google Scholar
  • View reference in article