Warning! Note that elemental potentials taken from alloy descriptions may not work well for the pure species. This is particularly true if the elements were fit for compounds instead of being optimized separately. As with all interatomic potentials, please check to make sure that the performance is adequate for your problem.
Citation: P. Tiwary, A. Walle, B. Jeon, and N. Grønbech-Jensen (2011), "Interatomic potentials for mixed oxide and advanced nuclear fuels", Physical Review B, 83(9), 094104. DOI: 10.1103/physrevb.83.094104.
Abstract: We extend our recently developed interatomic potentials for UO2 to the fuel system (U,Pu,Np)O2. We do so by fitting against an extensive database of ab initio results as well as to experimental measurements. The applicability of these interactions to a variety of mixed environments beyond the fitting domain is also assessed. The employed formalism makes these potentials applicable across all interatomic distances without the need for any ambiguous splining to the well-established short-range Ziegler-Biersack-Littmark universal pair potential. We therefore expect these to be reliable potentials for carrying out damage simulations (and molecular dynamics simulations in general) in nuclear fuels of varying compositions for all relevant atomic collision energies.
Citation: P. Tiwary, A. van de Walle, and N. Grønbech-Jensen (2009), "Ab initio construction of interatomic potentials for uranium dioxide across all interatomic distances", Physical Review B, 80(17), 174302. DOI: 10.1103/physrevb.80.174302.
Abstract: We provide a methodology for generating interatomic potentials for use in classical molecular-dynamics simulations of atomistic phenomena occurring at energy scales ranging from lattice vibrations to crystal defects to high-energy collisions. A rigorous method to objectively determine the shape of an interatomic potential over all length scales is introduced by building upon a charged-ion generalization of the well-known Ziegler-Biersack-Littmark universal potential that provides the short- and long-range limiting behavior of the potential. At intermediate ranges the potential is smoothly adjusted by fitting to ab initio data. Our formalism provides a complete description of the interatomic potentials that can be used at any energy scale, and thus, eliminates the inherent ambiguity of splining different potentials generated to study different kinds of atomic-materials behavior. We exemplify the method by developing rigid-ion potentials for uranium dioxide interactions under conditions ranging from thermodynamic equilibrium to very high atomic-energy collisions relevant for fission events.
GULP (2011--Tiwary-P--MOx--GULP--ipr1)
Notes: These files were posted on 17 June 2011 with the approval of Pratyush Tiwary and Axel van de Walle (California Institute of Technology). 30 Jan. 2012 Update: These files are identified as having problems with thermal expansion in UO2 and are superseded by the next GULP implementation below. File(s): superseded
Notes: 30 Jan. 2012 Update: These files (version 2.0) were provided by Pratyush Tiwary (California Institute of Technology) and posted with his permission. According to him, this version corrects the "U-U interaction term (truncated at 4 Angstroms) to resolve problems with thermal expansion in UO2." Additional information is located in readme_v2.txt File(s):