Abstract

Iron-technetium alloys are of relevance to the development of waste forms for disposition of radioactive technetium-99 obtained from spent nuclear fuel. Corrosion of candidate waste forms is a function of the local cohesive energy () of surface atoms. A theoretical model for calculating is developed. Density functional theory was used to construct a modified embedded atom (MEAM) potential for iron-technetium. Materials properties determined for the iron-technetium system were in good agreement with the literature. To explore the relationship between local structure and corrosion, MEAM simulations were performed on representative iron-technetium alloys and intermetallics. Technetium-rich phases have lower , suggesting that these phases will be more noble than iron-rich ones. Quantitative estimates of based on numbers of nearest neighbors alone can lead to errors up to 0.5 eV. Consequently, atomistic corrosion simulations for alloy systems should utilize physics-based models that consider not only neighbor counts, but also local compositions and atomic arrangements.

1. Introduction

The development of long-term containment strategies for spent nuclear fuel requires a combination of careful experiments and theoretical studies so that the realistic lifetimes of these strategies can be confidently predicted. Currently, a number of different containment strategies are being considered, one of which involves the storage of certain fission products within a metallic alloy waste form [1]. Chief among these fission products is technetium-99, which has a half-life of ~105 years [2]. Studies have demonstrated that technetium can be alloyed with stainless steel with significant mass-loadings and that these waste forms have considerable corrosion resistance [3].

Predictions of the long-term stability of this material can be obtained through the coupling of the results of accelerated testing studies with a rigorous, physics-based theoretical model. Corrosion is believed to be the chief process by which stored, spent nuclear fuel will degrade over time and have the potential to release radionuclides to the environment [1, 4, 5]. Modeling corrosion at long-time scales requires a thorough and fundamental understanding of the electrochemical stability of such systems, particularly the surface phenomena (Figure 1) [6]. Surface processes that are important to consider include passivation (formation of protective oxide films), depassivation (the rupture of such films by localized dissolution, mechanical damage, or presence of deleterious chemistries), mass-transport (at the metal/oxide interfaces, through the oxide, and at the oxide/solution interface), and bare-metal dissolution. The construction of a comprehensive corrosion model requires the elucidation of numerous processes acting across multiple length-scales and materials types (metals, oxides, and fluids).

One of the tools that can be applied to assist in the construction of lifetime prediction models is atomistic simulation. Contemporary atomistic modeling tools in the solid-state and molecular sciences allow the first-principles calculation of thermodynamic energies, activation barriers associated with phase-transformations, and the prediction of phase diagrams. First-principles modeling has been used by theorists at UNLV (Weck et al.) to explore the stability of Tc nanoparticles [7], and the nature of Tc-Zr alloy materials [8]. First-principles models have also been used to address problems in Tc corrosion. In a previous study, for example, we considered the reaction of water with the metallic surfaces containing iron and technetium [9]. The reaction of water to form oxygen monolayers on metal surfaces is one of the first steps in the formation of metal oxide layers. Using density functional theory, we were able to compute the energy associated with the reaction of water with the surface to form hydroxide or oxide monolayers, and using this energy, we computed the surface coverage as a function of the electrochemical pH and potential.

In the present work, we have applied atomistic modeling to illuminate a second set of corrosion reactions relevant to the iron-technetium system. The tendency for an atom to dissolve or be released from the metal by other means (such as diffusion through the passive film or ligand complexation) is related to the intrinsic stability of the atom in its local surface configuration (composition, coordination number, and local orientation, see Figure 2). Metal atom surface stability is a fundamental property of metal atoms at the surface and can be quantified by the atomic cohesive energy. The cohesive energy of an atom can be determined using an embedding functional approach, such as the modified embedded atom method developed by Baskes and Johnson [10]. By developing a modified embedded atom potential for the iron-technetium binary system, we have then computed cohesive energies for surface metal atoms in this system and, subsequently, identified which atoms and configurations of atoms will be more prone to release from the system due to corrosion effects. We plan to use this tool in future simulations of the active dissolution of candidate iron-technetium alloy waste forms.

2. Computational Details

Our strategy for computing the dissolution tendency of metal atoms in the iron-technetium binary system required the calculation of cohesive energies using the modified embedded atom method. It was necessary for us to construct a modified embedded atom potential that is accurate for the iron-technetium interactions. Modified embedded atom potentials for binary systems can be constructed by consulting accurate materials properties for the individual metals and their intermetallic phases, and these can be determined either from experiments, or higher-level theoretical calculations [10]. In the absence of experimental data and materials properties for the intermetallic phase, we have used density functional theory to compute the enthalpy of mixing, bulk modulus, and lattice parameter of a hypothetical B2 phase. We have also performed first-principles calculations of the moduli and lattice parameters of hexagonal close-packed technetium metal. These data-points were then used to construct the rest of the intermetallic iron-technetium modified embedded atom potential, which was then applied to investigate cohesive relations for some model iron-technetium surface systems.

2.1. Density Functional Theory (DFT)

Density functional theory is a powerful theoretical technique developed by Hohenberg and Kohn for the determination of the electronic structure associated with a given crystallographic and/or molecular system [11, 12]. This method is based upon the quantum mechanical physics that describe the interactions occurring between systems of electrons and nuclei. Because the physics behind this method operates on such a fundamental level, the method can very accurately describe the properties of materials that relate to the electronic structure. Such properties include the binding energy that holds a material together and the variations in this binding energy that are due to distortions in the structure or composition of a material.

Density functional theory has been implemented in numerous software packages that are commercially available. We utilized the Vienna Ab-initio Simulation Package (VASP) v.4.6.36 to solve for the electronic structures and energies for the reference structures studied in the iron-technetium binary system [13, 14]. Fully relaxed density functional theory calculations were performed for technetium in the hexagonal close-packed configuration, and materials properties such as the elastic moduli were evaluated by performing small distortions on the unit cell configuration [15]. Similar calculations were performed on a hypothetical iron-technetium (space group B2, Figure 3(b)) intermetallic to provide reference points for the intrinsic properties of iron-technetium metallic systems. Although a sigma-phase intermetallic for the iron-technetium system has been identified by Darby, the low-symmetry disordered structure makes it nontrivial for use in the development of a potential [16]. Furthermore, the materials properties, other than the lattice constant, of the Darby phase have not yet been determined.

Systems parameters within VASP that were specified to expedite the computationally intensive quantum mechanical calculations are described briefly as follows: the projector augmented wave method (PAW) was used to treat the core electrons and core-valence interactions [13]. Plane-waves were used as a basis set to expand the electronic structure, with energies taken up to 400 eV. The Perdew-Burke-Erzenrhof functional (PBE) was used to solve the exchange-correlation contribution to the electronic energy [17]. For the unit cell of hexagonal close-packed technetium, a 13 × 13 × 7 Monkhorst-Pack k-point mesh was used, and for the B2 structure of the hypothetic iron-technetium intermetallic, a 17 × 17 × 17 k-point mesh was used [18]. Electron smearing was used to accelerate convergence using the Methfessel-Paxton method [19].

2.2. Modified Embedded Atom Method (MEAM)

The modified embedded atom method is, in essence, a simplified approach to evaluating the energies and forces on atoms that is achieved by reducing and collecting the equations of density functional theory into a sum of single atom and pairwise expressions [10]. The energy is evaluated by constructing “embedding densities” contributed by all the nearby atoms, based upon simple parameterized forms. Similarly, the pairwise interactions are parameterized by comparison with experimentally or theoretically derived equations of state. The procedure for developing such potentials has been well described in the literature [10]. Because this method is much more computationally expedient than density functional theory, it can be used to simulate materials systems of greater spatial and temporal extent.

A modified embedded atom potential for iron was previously developed by Lee and Baskes [20] and was adopted for this work. Prior to this study, there were no preexisting potentials for technetium and the pairwise iron-technetium interaction. Hence, development of these potentials was an important step in this project. The materials properties evaluated using density functional theory were used to determine the necessary modified embedded atom parameters.

The large-scale atomic/molecular massively parallel simulator (LAMMPS), developed at Sandia National Laboratory, was used to perform total energy computations, conjugate gradient optimizations, and molecular dynamics simulations of iron-technetium systems using the modified embedded atom method [21]. To simulate a diversity of iron-technetium surfaces, spherical particles were simulated. Particles of hexagonal close-packed technetium, body-centered cubic iron and B2 iron-technetium were constructed. Technetium particles with iron skins, and iron particles with technetium skins, as well as technetium particles with iron adatoms and iron particles with technetium adatoms were simulated. It was important to consider such a diversity of surface configurations in order to span the entire range of possible local surface configurations that may be relevant to waste form corrosion scenarios. The energy of the surface atoms in these systems was calculated following a general optimization of the atomic coordinates on the basis of a local energy minimization procedure. Furthermore, molecular dynamics simulations at 1000 K were also performed to anneal the alloyed particles, thereby allowing even more configurations to be sampled.

The modified embedded atom method involves the computation of system energies by summation of individual and pairwise atomic cohesive energies using the formulae with being the embedding function of the density at atom and being the pair function between atoms and . The density is composed of a sum of spherical harmonic terms, as defined above, and the pairpotential is related to the Rose equation of state, also defined above. Using this formalism, it is, therefore, possible to distinguish between the local energies of surface atoms in different local configurations. Following the simulations described above, the energies of the surface atoms were tabulated and plotted as a function of the number of nearest technetium atoms versus the number of nearest iron atoms. In this way, we were able to discern the physics underlying the relative stabilities of metal atoms in different local orientations.

Attempts to estimate cohesive energies using counts of the number of nearest neighbors have been used in kinetic Monte Carlo simulations of alloy dissolution and nanoporous dealloying [22, 23]. However, in those cases, simple linear relations were assumed to construct the local energy. On the other hand, a modified embedded atom method was combined with a Metropolis Monte Carlo approach by Lillard et al. to simulate crystallographic pitting in hexagonal close-packed metals [24]. By direct computation of the cohesive energies for the alloy systems, in this work, we can directly compare the linear assumptions with the modified embedded atom method energies and determine which atomic configurations and phases will have the greatest impact on the lifetime performance of these alloys.

3. Results and Discussion

3.1. Density Functional Theory Computations of Material Properties

Density functional theory computations were performed for the hexagonal close-packed phase of technetium (Figure 3(a)). Following the procedure described by Vitos for hexagonal close-packed materials, the elastic moduli were computed via the finite distortion method [15]. The elastic moduli are reported in Table 1 and compared with the existing experimental values.

Examination of Table 1 indicates that the materials properties calculated using density functional theory for technetium are in reasonable agreement with the experimental results. The moduli are all consistently higher than the experimental values, with the calculated results being closer to the measurements made by Love et al. [25]. We attribute this to the fact that moduli determined from experiment will be sensitive to the temperature at which the experiment is performed. Guillermet reports values determined from experiments performed on Tc metal at room temperature, [26] which will be softer than the values by Love et al. obtained at 4 K, which in turn will be softer than the density functional theory values which correspond to the moduli obtained when distorting the lattice about the equilibrium positions at absolute zero.

While the B2 iron-technetium intermetallic is not a known phase for the Tc-Fe binary system, it has composition within the known range of the iron-technetium sigma-phase intermetallic (40%–60%), and its symmetry is consistent with the body-centered cubic phase of iron (Figure 3(b)) [16, 27]. Furthermore, the modified embedded atom method has been designed to use this high symmetry phase as a reference for computing the intermetallic terms in the potential energy expression. To develop this potential for the B2 phase, it is not necessary to compute all the elastic moduli (as was the case for technetium). The advantage in this case is due to the way that the pairwise intermetallic potentials are derived in the modified embedded atom method: it is only necessary to compute the lattice constant, formation energy, and bulk modulus of the intermetallic phase. We find that the energy of mixing for iron and technetium to produce the B2 phase is +0.325 eV, and that the bulk modulus for this phase is 290 GPa, similar to the bulk modulus for hexagonal close-packed technetium metal.

3.2. Development of the Modified Embedded Atom Potential

Following the procedure given by Baskes and Johnson [10], we utilize the density functional theory data in Table 1 to develop the modified embedded atom potentials for technetium and the iron-technetium intermetallic and alloy phases. In addition to the moduli and lattice parameters in Table 1, estimates of the vacancy formation energy and cohesive energy were required. The cohesive energy for technetium has been computed from first-principles previously [28]. Since no estimate of the vacancy formation energy of technetium could be found, the vacancy formation energy of rhenium was used instead [10]. Rhenium is commonly used as a surrogate material for investigations of technetium, for reasons of handling and similarity in the materials properties [29]. Vacancy energies can be calculated from first-principles, but they require supercell calculations that are not always practical. In future work, we hope to augment this value with a rigorous vacancy energy calculation. For iron, we utilized the existing potential for body-centered cubic iron created by Lee and Baskes [20]. The resulting MEAM parameters are given in Table 2.

3.3. Application of MEAM Potential to Map the Surface Stability of Tc-Fe Alloys

Spherical clusters of technetium, iron, iron-technetium, and various core-shell combinations were modeled using the LAMMPS package. Spherical clusters based on technetium were generated based on a hexagonal close-packed parent lattice. Spherical clusters based on iron were generated using a body-centered cubic lattice, and spherical clusters based on iron-technetium were generated using the hypothetical iron-technetium B2 structure. Geometry optimizations as well as molecular dynamics simulations at 1000 K for 150 ps were performed on these systems so that the atoms could migrate to low-energy sites. Due to the geometry of the spherical particles, multiple different surface planes (and hence atomic configurations for the surface atoms) were sampled within each simulation. Example clusters are shown in Figure 4.

For the purposes of probing a greater set of surface ensembles, each of the spherical particles shown in Figure 4 was subject to various additional modifications: surface atoms were replaced either completely by technetium, or completely by iron, and randomly selected atoms were switched from one type (Tc, Fe) to another (Fe, Tc). The cohesive energies calculated for all of the surface atoms within these various ensembles were collated and assembled into a master database. The database consisted of the atom type, the number of nearest neighbors of type Tc, and the number of nearest neighbors of type Fe as well as the cohesive energy for that atom, as determined by the modified embedded atom potential calculation. For each atom in the simulation, the number of nearest neighbors and their type were determined by counting the number of neighbors with a bond-length less than the sum of the metallic radii for the two atom types (i.e., + , or + , or + ) plus 10%. Metallic radii were given by half of the nearest neighbor distance in the pure metal (hexagonal close-packed technetium or body-centered cubic iron).

Following this compilation process, we plotted the cohesive energy of technetium and iron atoms that occupy surface sites in the metal versus the number of nearest neighbors of technetium (#Tc) and the number of nearest neighbors of iron (#Fe). These cohesive energies are plotted in contour form for the two variables (#Tc and #Fe) in Figure 5. Below the cohesive energies I have also plotted the standard deviation in this cohesive energies, based on the variations found in the overall database.

The cohesive energy plots shown in Figure 5 (top) indicate that surface atoms with more nearest neighbors of either iron or technetium possess a lower cohesive energy and, hence, a higher stability as regards to corrosion. This observation is in agreement with recent studies on copper nanoparticles that indicate that an increase in local coordination results in an increase in the cohesive stability of an atom [30]. Areas shown in “hot” colors (reds) correspond to surface environments that will be more likely to corrode. Areas shown in “cool” colors (blues) will be stabilized through stronger local bonding and, therefore, less likely to corrode.

What is also apparent from Figure 5 is that the cohesive energy of a surface atom is not only a function of the number of nearest neighbors, but also on the type of those nearest neighbors. That is, the degree of stabilization granted to an atom by increase coordination will vary depending upon what kind of neighbors the atom is bonding to. For both iron and technetium atoms, surface bonding to technetium provides a much more significant degree of stabilization than coordination in an iron-rich surface environment. Thus, technetium-rich surfaces would be expected to be “more noble” than well-mixed or iron-rich surface phases. Recent measurements of iron-technetium candidate alloy waste form corrosion support this observation in that iron-rich phases are found to have higher corrosion rates than the intermetallic phases [1].

In the lower portion of Figure 5, we have also plotted the standard deviations in the collated surface atom cohesive energies. This data was plotted to highlight areas in which the assumption that cohesive energies can be determined simply by considering the number and type of nearest neighbors has the potential to break down. For iron surrounded by mostly iron atoms, the standard deviation is relatively small. In this case, an analytical expression for the variation in the cohesive energy as a function of the number of nearest neighbors may still apply. However, as the number of technetium neighbors coordinating the iron atom increases, the standard deviation in the cohesive energy can be as high as 0.5 eV. This implies that the orientation, or the ordering, of these nearest neighbors, which is not captured in this two-dimensional representation, is also a significant factor in the stability of those exposed surface atoms. For technetium atoms, the standard deviations are significantly larger than for the iron atoms. This result suggests that the local cohesive energy of technetium atoms is more sensitive to the local configuration of atoms than iron. It may also be an indicator that further testing and development of the modified embedded atom method potential for technetium is required.

4. Conclusions

In this work, we have performed density functional theory computations for metallic technetium and iron-technetium intermetallic systems. We have used materials properties obtained from the density functional theory computations to develop a modified embedded atom method potential for the iron-technetium system. The materials properties calculated using density functional theory were in good agreement with the available experimental literature and were, therefore, used to parameterize the MEAM potential.

Simulations containing model surfaces of iron-technetium alloys in the iron-rich, technetium-rich, and 50-50 iron-technetium intermetallic limit were then performed. The simulations were designed to capture the variations in surface cohesive energies that occur as iron or technetium atoms are placed in different local compositional and coordination environments. We collated the local surface atom cohesive energies obtained from the various computations so that we could investigate how the local environments of iron or technetium atoms at the surface of hypothetical iron-technetium alloy waste forms may influence their corrosion behavior. The simulations indicated that the local bonding arrangements are very critical to the ultimate stability of these surface atoms. Generally speaking, increases in local coordination lead to more stable surface atoms; that is, more bonds to nearest neighbors will increase the energy required to dissolve or remove an atom from the alloy surface. Another observation that is revealed by this study is that technetium-rich environments are more noble, not only for technetium atoms on the surface, but also for iron atoms.

Our work also shows that a consideration of the local bonding environment, in terms of numbers of nearest neighbors, is not sufficient for computing or estimating the local cohesive energy and that variations in the particular way the atoms are situated about the central atom can lead to perturbations in the cohesive energy of up to 0.5 eV (as estimated by one standard deviation from the mean). Therefore, although linear relationships for dissolution and diffusion energies have been assumed in previous works, significant errors in these simulations may apply.

The modified embedded atom potential for iron-technetium alloys that we have developed in this work is currently being implemented within a kinetic Monte Carlo scheme, which we will apply in future work for the simulation of relative corrosion rates for iron and technetium in candidate metal alloy waste forms, as well as the differences in corrosion morphology and extending the potential to include other alloy components such as molybdenum and chromium. This technique will allow us to investigate the relative corrosion properties of different intermetallic and alloy phases and the roles of electrochemical potential and pH, and provide a starting point for more complex models that combine alternative surface processes such as oxidation, grain-boundary dissolution, and surface adsorption.

Acknowledgments

The authors acknowledge funding through the Fundamental Waste Form Science campaign under the auspices of the US DOE Fuel Cycle R&D under the direction of John Vienna, Waste Forms and Separations Campaign Manager. Helpful discussions with Dave Moore, Gordon Jarvinen, Scott Lillard, and Dave Kolman at Los Alamos National Laboratory are also acknowledged. The Los Alamos National Laboratory is operated by Los Alamos National Security LLC for the National Nuclear Security Administration of the U.S. Department of Energy under Contract no. DE-AC52-06NA25396. A portion of the research was performed using EMSL, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory.