Molecular dynamics calculations of the thermal expansion properties and melting points of Si and Ge

The thermal expansion properties and melting points of silicon and germanium are calculated using molecular dynamics simulations within the density functional theory framework. An isothermal-isobaric (NPT) ensemble is considered in a periodic system with a relatively small number of particles per unit cell to obtain the thermal expansion data over a range of temperatures, and it is found that the calculated thermal expansion coefficients and bond lengths agree well with experimental data. Also, the positions of discontinuities in the potential energy as a function of temperature are in good agreement with the experimental melting points