Development of molecular dynamics potential for uranium silicide fuels
The use of uranium-silicide (U-Si) in place of uranium dioxide (UO{sub 2}) is one of the promising concepts being proposed to increase the accident tolerance of nuclear fuels for light water reactors. First-principles based density functional theory (DFT) electronic structure calculations may provide the most reliable information. However, it is probably not possible to obtain kinetic information such as amorphization under irradiation directly from DFT simulations due to size and time limitations. Thus, a more feasible way may be to employ molecular dynamics (MD) simulation. Unfortunately, so far no MD potential is available for U-Si to discover the underlying mechanisms. Here, we present our recent progress in developing a U-Si potential. The development is based on the Tersoff type potentials for single element U and Si. The Si potential is taken from the literature and a Tersoff type U potential is developed in this project. With the primary focus on the U{sub 3}Si{sub 2} phase, some other U-Si systems such as U{sub 3}Si are also included as a test of the transferability of the potentials for binary U-Si phases. Based on the potentials for unary U and Si, two sets of parameters for the binary U-Si system are developed using the Tersoff mixing rules and the cross-term fitting, respectively. The cross-term potential is found to give better results on the enthalpy of formation, lattice constants and elastic constants than those produced by the Tersoff mixing potential, with the reference data taken from either experiments or DFT calculations. In particular, the results on the formation enthalpy and lattice constants for the U{sub 3}Si{sub 2} phase and lattice constants for the high temperature U{sub 3}Si (h-U{sub 3}Si) phase generated by the cross-term potential agree well with experimental data. Reasonable agreements are also reached on the elastic constants of U{sub 3}Si{sub 2}, on the formation enthalpy for the low temperature U{sub 3}Si (m-U{sub 3}Si) and h-U{sub 3}Si phases, and on the lattice constants of m-U{sub 3}Si phase. All these phases are predicted to be mechanically stable. The unary U potential is also tested for three metallic U phases (α, β, γ). The potential is found capable to predict the cohesive energies well against experimental data for all three phases. It matches reasonably with previous experiments on the lattice constants and elastic constants of αU. (authors)