Repository logo
 

Fitting Interatomic Potentials to Reproduce Phase Transitions


Type

Thesis

Change log

Authors

Abstract

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. 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 calculated could be upgraded in future.

Description

Date

2019-03-29

Advisors

Csányi, Gábor

Keywords

interatomic potentials, phase diagrams, phase transitions, copper, gold, titanium, nested sampling, Nelder-Mead method, fitting interatomic potentials, CuAu, Ti, melting curve, eutectic

Qualification

Doctor of Philosophy (PhD)

Awarding Institution

University of Cambridge
Sponsorship
EPSRC Centre for Doctoral Training in Computational Methods for Materials Science under grant number EP/L015552/1. Computing time on the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk) via the UKCP consortium funded by EPSRC grants ref EP/K013564/1 and ref EP/P022561/1. Additional funding from King's College (Cambridge) after the third year of the PhD.