× Notice! This site is currently being redesigned. Please let us know any feedback on the new design or if you find something incorrect/not working.
× Cool! Click on a potential's identifier for computed properties.
 
Citation: J.R. Fernández, and M.I. Pascuet (2014), "On the accurate description of uranium metallic phases: a MEAM interatomic potential approach", Modelling and Simulation in Materials Science and Engineering, 22(5), 55019. DOI: 10.1088/0965-0393/22/5/055019.
Abstract: A new interatomic potential in the framework of the modified embedded atom method (MEAM) to model U metal is presented. The potential acceptably reproduces the lattice parameters and cohesive energy of the orthorhombic αU. The relative stability of the experimentally observed phase at low temperatures with respect to several other structures (bct, bcc, simple cubic, tetragonal β Np, fcc and hcp) is also taken into account. Intrinsic point defect properties compare reasonably well with data from the literature. To determine the quality of the interaction, the potential is used to study a number of properties for the pure metal at finite temperatures and the results are compared with the available data. The obtained allotropic αU ↔ γU transformation and melting temperatures are in good agreement with experimental values. Based on the simulations, a new αU ↔ γU transformation mechanism is proposed.

LAMMPS pair_style meam (2014--Fernandez-J-R--U--LAMMPS--ipr1)
Notes: These files were sent by M.I. Pascuet (CONICET) on 25 Apr. 2016 and posted with her permission.
File(s):
Citation: D.E. Smirnova, S.V. Starikov, and V.V. Stegailov (2011), "Interatomic potential for uranium in a wide range of pressures and temperatures", Journal of Physics: Condensed Matter, 24(1), 15702. DOI: 10.1088/0953-8984/24/1/015702.
Abstract: Using the force-matching method we develop an interatomic potential that allows us to study the structure and properties of α-U, γ-U and liquid uranium. The potential is fitted to the forces, energies and stresses obtained from ab initio calculations. The model gives a good comparison with the experimental and ab initio data for the lattice constants of α-U and γ-U, the elastic constants, the room-temperature isotherm, the normal density isochore, the bond-angle distribution functions and the vacancy formation energies. The calculated melting line of uranium at pressures up to 80 GPa and the temperature of the α-γ transition at 3 GPa agree well with the experimental phase diagram of uranium.
Citation: D.E. Smirnova, S.V. Starikov, and V.V. Stegailov (2012), "Interatomic potential for uranium in a wide range of pressures and temperatures", Journal of Physics: Condensed Matter, 24(14), 149501. DOI: 10.1088/0953-8984/24/14/149501.

LAMMPS pair_style eam/alloy (2012--Smirnova-D-E--U--LAMMPS--ipr1)
Notes: This file was sent by Daria Smirnova (Joint Institute for High Temperatures, Russian Academy of Sciences) and posted on 18 Jan. 2013 with her permission.
File(s):
 
Citation: M.I. Pascuet, and J.R. Fernández (2015), "Atomic interaction of the MEAM type for the study of intermetallics in the Al-U alloy", Journal of Nuclear Materials, 467, 229-239. DOI: 10.1016/j.jnucmat.2015.09.030.
Abstract: Interaction for both pure Al and Al–U alloys of the MEAM type are developed. The obtained Al interatomic potential assures its compatibility with the details of the framework presently adopted. The Al–U interaction fits various properties of the Al2U, Al3U and Al4U intermetallics. The potential verifies the stability of the intermetallic structures in a temperature range compatible with that observed in the phase diagram, and also takes into account the greater stability of these structures relative to others that are competitive in energy. The intermetallics are characterized by calculating elastic and thermal properties and point defect parameters. Molecular dynamics simulations show a growth of the Al3U intermetallic in the Al/U interface in agreement with experimental evidence.

LAMMPS pair_style meam (2015--Pascuet-M-I--Al-U--LAMMPS--ipr2)
Notes: These files were sent by M.I. Pascuet (CONICET) on 22 Apr. 2016 and posted with her permission.
File(s):
 
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), 94104. 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
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


GULP
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):
 
Citation: S.V. Starikov, L.N. Kolotova, A.Y. Kuksin, D.E. Smirnova, and V.I. Tseplyaev (2018), "Atomistic simulation of cubic and tetragonal phases of U-Mo alloy: Structure and thermodynamic properties", Journal of Nuclear Materials, 499, 451-463. DOI: 10.1016/j.jnucmat.2017.11.047.
Abstract: We studied structure and thermodynamic properties of cubic and tetragonal phases of pure uranium and U-Mo alloys using atomistic simulations: molecular dynamics and density functional theory. The main attention was paid to the metastable γ0-phase that is formed in U-Mo alloys at low temperature. Structure of γ0-phase is similar to body-centered tetragonal (bct) lattice with displacement of a central atom in the basic cell along [001] direction. Such displacements have opposite orientations for part of the neighbouring basic cells. In this case, such ordering of the displacements can be designated as antiferro-displacement. Formation of such complex structure may be interpreted through forming of short U-U bonds. At heating, the tetragonal structure transforms into cubic γs-phase, still showing ordering of central atom displacements. With rise in temperature, γs-phase transforms to γ-phase with a quasi body-centered cubic (q-bcc) lattice. The local positions of uranium atoms in γ-phase correspond to γs-phase, however, orientations of the central atom displacements become disordered. Transition from γ0 to γ can be considered as antiferro-to paraelastic transition of order-disorder type. \n\nThis approach to the structure description of uranium alloy allows to explain a number of unusual features found in the experiments: anisotropy of lattice at low temperature; remarkably high self-diffusion mobility in γ-phase; decreasing of electrical resistivity at heating for some alloys. In addition, important part of this work is the development of new interatomic potential for U-Mo system made with taking into account details of studied structures.

LAMMPS pair_style adp (2018--Starikov-S-V--U-Mo--LAMMPS--ipr1)
Notes: These files were sent by S.V. Starikov (Joint Institute for High Temperatures, Russian Academy of Sciences) on 3 Dec. 2017 and posted with his permission.
File(s): superseded


LAMMPS pair_style adp (2018--Starikov-S-V--U-Mo--LAMMPS--ipr2)
Notes: This file was recieved by Sergei Starikov on August 2, 2018. He noted: "In the previous version, one function for Mo-Mo interaction had poor smoothing with r=r_cut. This "bug" led to the sake of the energy conservation during high-temperature (or long) calculations for pure Mo or U-Mo alloy. I made an additional check and found a small potential compilation issue leading to this energy drift. I fixed the file, so now it can be safely used with the time-step = 0.5 fs, even for long runs at high temperatures (the test was made for 1500K)."
File(s):
 
Citation: D.E. Smirnova, A.Y. Kuksin, S.V. Starikov, V.V. Stegailov, Z. Insepov, J. Rest, and A.M. Yacout (2013), "A ternary EAM interatomic potential for U-Mo alloys with xenon", Modelling and Simulation in Materials Science and Engineering, 21(3), 35011. DOI: 10.1088/0965-0393/21/3/035011.
Abstract: A new interatomic potential for a uranium–molybdenum system with xenon is developed in the framework of an embedded atom model using a force-matching technique and a dataset of ab initio atomic forces. The verification of the potential proves that it is suitable for the investigation of various compounds existing in the system as well as for simulation of pure elements: U, Mo and Xe. Computed lattice constants, thermal expansion coefficients, elastic properties and melting temperatures of U, Mo and Xe are consistent with the experimentally measured values. The energies of the point defect formation in pure U and Mo are proved to be comparable to the density-functional theory calculations. We compare this new U–Mo–Xe potential with the previously developed U and Mo–Xe potentials. A comparative study between the different potential functions is provided. The key purpose of the new model is to study the atomistic processes of defect evolution taking place in the U–Mo nuclear fuel. Here we use the potential to simulate bcc alloys containing 10 wt% of intermetallic Mo and U2Mo.

LAMMPS pair_style eam/alloy (2013--Smirnova-D-E--U-Mo-Xe--LAMMPS--ipr1)
Notes: This file was sent by Daria Smirnova (Joint Institute for High Temperatures, Russian Academy of Sciences) and posted on 14 March 2013 with her permission.
File(s):
 
Citation: V.I. Tseplyaev, and S.V. Starikov (2016), "The atomistic simulation of pressure-induced phase transition in uranium mononitride", Journal of Nuclear Materials, 480, 7-14. DOI: 10.1016/j.jnucmat.2016.07.048.
Abstract: In this work we studied the pressure-induced phase transition between different structures of uranium mononitride: cubic Fm-3m-structure and rhombohedral R-3m-structure. We used molecular dynamics together with a new interatomic potential developed for this purpose. We estimated phase diagram of uranium mononitrde in a wide range of temperature and pressure using thermodynamic and mechanical criteria of stability. From simulations we see that at zero temperature the phase transition Fm-3m -> R-3m takes place at pressure about 35 GPa, which agrees well with the available experimental and theoretical data. Results of the calculations show that the lattice of rhombohedral phase becomes close to cubic structure with increase in temperature.

LAMMPS pair_style adp (2016--Tseplyaev-V-I--U-N--LAMMPS--ipr1)
Notes: These files were sent by S.V. Starikov (Joint Institute for High Temperatures, Russian Academy of Sciences) on 3 Nov. 2016 and posted with his permission.
File(s):
 
Citation: B. Beeler, M. Baskes, D. Andersson, M.W.D. Cooper, and Y. Zhang (2017), "A modified Embedded-Atom Method interatomic potential for uranium-silicide", Journal of Nuclear Materials, 495, 267-276. DOI: 10.1016/j.jnucmat.2017.08.025.
Abstract: Uranium-silicide (U-Si) fuels are being pursued as a possible accident tolerant fuel (ATF). This uranium alloy fuel benefits from higher thermal conductivity and higher fissile density compared to uranium dioxide (UO2). In order to perform engineering scale nuclear fuel performance simulations, the material properties of the fuel must be known. Currently, the experimental data available for U-Si fuels is rather limited. Thus, multiscale modeling efforts are underway to address this gap in knowledge. In this study, a semi-empirical modified Embedded-Atom Method (MEAM) potential is presented for the description of the U-Si system. The potential is fitted to the formation energy, defect energies and structural properties of U3Si2. The primary phase of interest (U3Si2) is accurately described over a wide temperature range and displays good behavior under irradiation and with free surfaces. The potential can also describe a variety of U-Si phases across the composition spectrum.

LAMMPS pair_style meam (modified)
Notes: These files were sent by B. Beeler (Idaho National Laboratory) on 21 Mar. 2018 and posted with his permission. Dr. Beeler noted that the provided MEAM parameter files also require the use of a MEAM modification file to be compiled with LAMMPS.
File(s):
 
Citation: A.P. Moore, B. Beeler, C. Deo, M.I. Baskes, and M.A. Okuniewski (2015), "Atomistic modeling of high temperature uranium\textendashzirconium alloy structure and thermodynamics", Journal of Nuclear Materials, 467, 802-819. DOI: 10.1016/j.jnucmat.2015.10.016.

Dynamo MEAM
Notes: These files were sent by Alexander Moore (Georgia Institute of Technology) on 13 Aug. 2015 and posted with his permission. He noted that "This is a MEAM potential for U, Zr, and U-Zr alloys. The files attached are the potential files for DYNAMO. It should be noted that use of this potential in LAMMPS requires LAMMPS to have a modified cut-off function before it is compiled." Update 27 April 2018: Publication information was added.
File(s):
 
Citation: A.E. Thompson, B. Meredig, M. Stan, and C. Wolverton (2014), "Interatomic potential for accurate phonons and defects in UO2", Journal of Nuclear Materials, 446(1-3), 155-162. DOI: 10.1016/j.jnucmat.2013.11.040.
Abstract: We have developed an improved uranium dioxide interatomic potential by fitting to forces, energies, and stresses of first principles molecular dynamics calculations via a genetic algorithm approach called Iterative Potential Refinement (IPR). We compare the defect energetics and vibrational properties of the IPR-fit potential with other interatomic potentials, density functional theory calculations, and experimental phonon dispersions. We find that among previously published potentials examined, there is no potential that simultaneously yields accurate defect energetics and accurate vibrational properties. In contrast, our IPR-fit potential produces both accurate defects and the best agreement with the experimental phonon dispersion and phonon density of states. This combination of accurate properties makes this IPR-fit potential useful for simulating UO2 in high temperature, defect-rich environments typical for nuclear fuel. Additionally, we verify that density functional theory with a Hubbard U correction accurately reproduces the experimentally derived UO2 phonon density of states.

GULP
Notes: These parameters were sent by Dr. Alexander Thompson and posted with his approval on 25 Nov. 2014.
Parameters:
O core 1.772816
O shell -3.737358
U core -2.069827
U shell 5.998911
spring
U 115.5906
O 261.0604
buck
U shell O shell 1062.316 0.398196 0.0 0.0 15.0
U shell U shell 183.6 0.420031 0.0 0.0 15.0
buck4
O shell O shell 10402.58 0.238539 81.75059 0.0 1.150671 2.22 2.460169 15.0

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.

Notes: 31 Jan. 2012 Update: This potential has been superseded by the (U,Pu,Np)O2 interatomic potential.

GULP
Notes: These were supplied by Pratyush Tiwary and Axel van de Walle (California Institute of Technology) and posted with their approval on 22 June 2011.
File(s): superseded


Date Created: October 5, 2010 | Last updated: October 05, 2018