Fitting Interatomic Potentials to Reproduce Phase Transitions

Martin Schlegel, Gábor Csányi, Apollo-University Of Cambridge Repository, Apollo-University Of Cambridge Repository
2020
Due to the high computational demand of quantum mechanical simulations, researchers still rely on models for atomic interactions to simulate a large number of atoms. These so called "interatomic potentials" take a variety of forms ranging from the simple Lennard-Jones potential to machine learning potentials. They are usually fitted for specific structures as well as pressure and temperature realms and do not necessarily reproduce the full experimentally known solid and liquid phase diagram.
more » ... d phase diagram. There have been limited trials of fitting the parameters of interatomic potentials to phase transitions. In this work, I present improvements upon two already existing potentials from the literature – a CuAu potential and a Ti potential. These potentials are modified so that they reproduce experimental phase transitions better than before. For the CuAu potential, I fit the eutectic melting curve using my own implementation of the Nelder-Mead algorithm. The Ti potential is modified to predict a triple point which is in qualitative agreement with experimental predictions. To calculate phase transitions, I use the nested sampling algorithm which was already established in the literature. The Ti nested sampling results are verified with thermodynamic integration. My Nelder-Mead implementation allows for semi-automatic submission and analysis of the nested sampling runs. Both potentials are improved considerably. The new CuAu potential now accurately reproduces the melting curve for the alloy compositions used in the fitting procedure. The new Ti potential now – as opposed to before the modifications – predicts a triple point which is in qualitative accordance with experimental predictions. My work demonstrates how the transferability of potentials can be improved by just slight modifications to the original parameters using experimental data of the phase diagram. The Nelder-Mead implementation used for the CuAu potential could be used to improve other potentials. To fit more complex systems, the way the objective function is calc [...]
doi:10.17863/cam.50769 fatcat:omx5ese3ozfflf7jcxfj4fmhme