Error estimation in molecular dynamics simulations: nitrogen impurities in iron with empirical potentials

  1. Error estimation in molecular dynamics simulations: nitrogen impurities in iron with empirical potentials

    17MODEV04 / Computational material research on the nanoscale, Model development
    Promotor(en): T. Verstraelen, S. Cottenier / Begeleider(s): S. De Waele, K. Lejaeghere, M. Cools-Ceuppens

    As electrons in matter obey the laws of quantum mechanics, one should in principle solve the Schrödinger equation when performing molecular simulations. On the nanoscale (tracking thousands of atoms during multiple nanoseconds) this is completely unfeasible. Even for approximate electronic structure methods, such as Density Functional Theory (DFT), the computational burden is too high to study phenomena such as diffusion of defects in solids. Empirical potentials can resolve this issue, because atomic forces are calculated with very simple analytical expressions. An example is the venerable embedded atom model, which was originally introduced to study hydrogen impurities in nickel [Daw1983] but was later extended for many other types of materials. The major disadvantage of empirical potentials is the loss of accuracy compared to electronic structure methods. Indeed, determining appropriate analytic expressions and corresponding parameters is a daunting task. This has led to a situation where force fields are often deemed unreliable, as one cannot anticipate to what extent a force field is able to correctly predict certain properties for the molecular systems at hand.
    A particular case where reliable force fields are required, is for the development of new steels by adding nitrogen. In conventional steel processing, nitrogen alloying is challenging due to its limited solubility during casting and solidification. Alternatively, nitriding as thermochemical treatment on the final material can be used to significantly improve surface and bulk properties. By studying the Fe-N metallurgy on the atomic scale, the aim is to gain a good understanding of the mechanisms of nitrogen diffusion and precipitation. To study the latter phenomenon, system sizes of thousands of atoms are required. Therefore, the development of an accurate empirical potential can provide a major leap forward in modeling the Fe-N system.

    Figure 1: (left) Fe16N2 Is a tetragonally distorted iron structure, where the nitrogen atoms form a sublattice on the interstitial sites. (right) Precipitation requires the diffusion of atoms to a local region with a different crystal structure and/or composition. This region is typically composed of many thousands of atoms.

    In recent years, quantifying uncertainty in simulation predictions has attracted the attention of several scientific communities. The Technology Roadmap for Computational Chemistry ( ranked as most critical the “lack of methods to estimate the intrinsic accuracy of calculations, leading to a critical barrier to the more widespread use of these methods for solving practical engineering problems.” For example, in this thesis it would be desirable to estimate the error on the simulated nitrogen diffusion constant, due to the approximations in the force field.

    Despite the fact that the necessity for uncertainty management in constructing force-field models has been expressed some years ago, the topic has received very little attention so far.


    The goal of this thesis is to construct two types of empirical potentials with controlled accuracy for the Fe-N system, and to apply them to study the diffusion and precipitation of nitrogen in iron and the geometry of nitrogen precipitates. (See Fig. 1.) Two standard empirical potential models will be used, EAM and ReaxFF. Both models have their strengths and weaknesses and are often criticized for their unclear accuracy in predictive simulations. For the calibration of parameters in these models, ample reference data is available in the frame of ongoing research at the Center for Molecular Modeling. Good starting points can be found in the literature, e.g. an EAM parameterization for Fe-P [Ackland2004] and initial parameters for ReaxFF can be found in a database curated by the company SCM in Amsterdam.

    The calibration of the parameters will in first instance be carried out with established methods, such as PotFit []. One major difficulty with the standard methods is that they try to find a single optimal set to describe the training data, which may be misleading. A good performance on the training set does not imply reliable results in applications afterwards. Bayesian inference is ideally suited to estimate errors (due to overfitting) of EAM parameters. [Federiksen2004] In this thesis, a recent extension of the Bayesian approach will be used, which was recently developed at the Center for Molecular Modeling. This extended Bayesian approach can also includes errors due to the approximations in the empirical potential and uncertainties in the reference data.

    The end result of the thesis would consist of relevant theoretical predictions with error estimates for nitrogen impurities in steel. For a fixed set of parameters in the empirical potential, diffusion constants can be computed using the Einstein relation (mean squared displacement versus time). By repeating such simulations with different "reasonable" parameter vectors, one obtains a distribution of predictions, to which standard statistics can be applied. A more challenging question is the structure of nitrogen precipitates. These domains are only a few atomic layers thick but several micrometers long. The structure of these domains at their boundaries is unclear and can be investigated with molecular dynamics simulations as well.

    Mobility: Mobility is not required for this thesis. The topic is related to a collaboration between the center for molecular modeling and OCAS, which are all Ghent-based.

    Motivation Appl. Phys.: This thesis requires physical insight because a good understanding of interatomic interactions is to construct physically sensible force field models. The final results of the thesis are used in the development of new types of steel, which is clearly an engineering problem.