A minimal Tersoff potential for diamond silicon with improved descriptions of elastic and phonon transport properties

Abstract

Silicon is an important material and many empirical interatomic potentials have been developed for atomistic simulations of it. Among them, the Tersoff potential and its variants are the most popular ones. However, all the existing Tersoff-like potentials fail to reproduce the experimentally measured thermal conductivity of diamond silicon. Here we propose a modified Tersoff potential and develop an efficient open source code called GPUGA (graphics processing units genetic algorithm) based on the genetic algorithm and use it to fit the potential parameters against energy, virial and force data from quantum density functional theory calculations. This potential, which is implemented in the efficient open source GPUMD (graphics processing units molecular dynamics) code, gives significantly improved descriptions of the thermal conductivity and phonon dispersion of diamond silicon as compared to previous Tersoff potentials and at the same time well reproduces the elastic constants. Furthermore, we find that quantum effects on the thermal conductivity of diamond silicon at room temperature are non-negligible but small: using classical statistics underestimates the thermal conductivity by about 10\% as compared to using quantum statistics.

arXiv comment9 pages, 6 figures

References in corpus (13)

• Restoring the density-gradient expansion for exchange in solids and surfaces
• Generalized gradient approximation for solids and their surfaces
• Temperature dependent effective potential method for accurate free energy calculations of solids
• Force and heat current formulas for many-body potentials in molecular dynamics simulation with applications to thermal conductivity calculations
• Efficient molecular dynamics simulations with many-body potentials on graphics processing units
• Homogeneous nonequilibrium molecular dynamics method for heat transport and spectral decomposition with many-body potentials
• Screened empirical bond-order potentials for Si-C
• A revisit to phonon-phonon scattering in single-layer graphene
• Global optimization of parameters in the reactive force field ReaxFF for SiOH
• Equivalence of the equilibrium and the nonequilibrium molecular dynamics methods for thermal conductivity calculations: From bulk to nanowire silicon
• Classical interaction potentials for diverse materials from ab initio data: a review of potfit
• First-principles calculation of intrinsic defect formation volumes in silicon
• An optimized interatomic potential for silicon and its application to thermal stability of silicene

Cited by in corpus (5)

• Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport
• GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
• Spectral Decomposition of Thermal Conductivity: Comparing Velocity Decomposition Methods in Homogeneous Molecular Dynamics Simulations
• Interpretation of apparent thermal conductivity in finite systems from equilibrium molecular dynamics simulations
• Effects of transverse geometry on the thermal conductivity of Si and Ge nanowires