Based on first-principles density-functional calculations, we have developed and tested a force-field for silicon, which can be used for molecular dynamics simulations and the calculation of its thermal properties. This force field uses the exact Taylor expansion of the total energy about the equilibrium positions up to 4th order. In this sense, it becomes systematically exact for small enough displacements, and can reproduce the thermodynamic properties of Si with high fidelity. Having the harmonic force constants, one can easily calculate the phonon spectrum of this system. The cubic force constants, on the other hand, will allow us to compute phonon lifetimes and scattering rates. Results on equilibrium Green-Kubo molecular dynamics simulations of thermal conductivity as well as an alternative calculation of the latter based on the relaxation-time approximation will be reported. The accuracy and ease of computation of the lattice thermal conductivity using these methods will be compared. This approach paves the way for the construction of accurate bulk interatomic potentials database, from which lattice dynamics and thermal properties can be calculated and used in larger scale simulation methods such as Monte Carlo.
Predicting CaO-(MgO)-Al2O3-SiO2 glass reactivity in alkaline environments from force field molecular dynamics simulations
