Fitting electron density as a physically sound basis for the development of interatomic potentials of complex alloys.

The development of new interatomic potentials to model metallic systems is a difficult task, due in part to the dependence between the parameters that describe the electron density and the short-range interactions. Parameter search methods are prone to false convergence. To solve this problem, we have developed a methodology for obtaining the electron density parameters independently of the short-range interactions, so that physically sound parameters can be obtained to describe the electron density, after which the short-range parameters can be fitted, thus reducing the complexity of the process and yielding better interatomic potentials. With the new method we can develop self-consistent, accurate force fields, using solely calculations, without the need to fit to experimental data. Density functional theory calculations are used to compute the observables with which the potential is fit. We applied the method to a Ni-based Inconel 625 superalloy (IN625), modelled here as Ni, Cr, Mo and Fe solid solution alloys. The capability of the force fields developed using this new method is validated, by comparing the structural and thermo-elastic properties predicted with the force fields, with the corresponding experimental data, both for single crystals and polycrystalline alloys.

Thermostructural behaviour of Ni-Cr materials: modelling of bulk and nanoparticle systems.

The thermostructural properties of Ni-Cr materials, as bulk and nanoparticle (NP) systems, have been predicted with a newly developed interatomic potential, for Ni/Cr ratios from 100/0 to 60/40. The potential, which has been fitted using experimental data and further validated using Density Functional Theory (DFT), describes correctly the variation with temperature of lattice parameters and the coefficient of thermal expansion, from 100 K to 1000 K. Using this potential, we have performed Molecular Dynamics (MD) simulations on bulk Ni-Cr alloys of various compositions, for which no experimental data are available. Similarly, NPs with diameters of 3, 5, 7, and 10 nm were studied. We found a very rapid convergence of NP properties with the size of the systems, showing already the 5 nm NPs with a thermostructural behaviour similar to the bulk. MD simulations of two 5 nm NPs show very little sintering and thermally induced damage, for temperatures between 300 K and 1000 K, suggesting that materials formed by agglomeration of Ni-Cr NPs meet the thermostructural stability requirements for catalysis applications.