Abstract: As interest in Cerium containing alloys and Al-Ce alloys in particular grows, the need to generate predictions of mechanical, thermodynamic and transport properties across a compositional space becomes more pronounced. The absence of a reliable interaction potential to be used in classical molecular dynamics (MD) simulation of γ-Ce is a fundamental bottleneck to subsequent alloy simulation. In this work, a Modified Embedded Atom Method (MEAM) potential to be used in MD simulation for γ-Ce has been generated. The parameterization process involved two steps. First, the iterative Latin hypercube sampling (LHS) method was used to generate MEAM parameter sets based on an automated optimization of energies and forces relative to a training set of small configurations evaluated through Density Functional Theory (DFT). Second, a human-guided optimization in that local parameter space was performed to refine the parameters to satisfy an array of mechanical, thermodynamic and transport properties available from the experimental literature. When used in an MD simulation, the resulting potential provides excellent estimates of all tested material properties across a broad temperature range. This γ-Ce potential, when combined with existing MEAM potential for Al, will form the necessary foundation for the subsequent development of a mixture potential, enabling the simulation of Al-Ce and Al-Ce-X alloys.

Abstract: We have developed a reactive force-field of the ReaxFF type for stoichiometric ceria (CeO2) and partially reduced ceria (CeO2–x). We describe the parametrization procedure and provide results validating the parameters in terms of their ability to accurately describe the oxygen chemistry of the bulk, extended surfaces, surface steps, and nanoparticles of the material. By comparison with our reference electronic structure method (PBE+U), we find that the stoichiometric bulk and surface systems are well reproduced in terms of bulk modulus, lattice parameters, and surface energies. For the surfaces, step energies on the (111) surface are also well described. Upon reduction, the force-field is able to capture the bulk and surface vacancy formation energies (Evac), and in particular, it reproduces the Evac variation with depth from the (110) and (111) surfaces. The force-field is also able to capture the energy hierarchy of differently shaped stoichiometric nanoparticles (tetrahedra, octahedra, and cubes), and of partially reduced octahedra. For these reasons, we believe that this force-field provides a significant addition to the method repertoire available for simulating redox properties at ceria surfaces.