Double lattice potential for molecular dynamics simulation of silicon with demonstrated validity

Hui Zhang,Chongyang Wei,Zhongwu Liu,Xichun Zhong,Dongling Jiao,Wanqi Qiu,Hongya Yu
DOI: https://doi.org/10.48550/arXiv.2002.00669
2020-08-11
Abstract:To reproduce the diamond structure of silicon, double lattice (DL) potential constructed from two interatomic potentials for face centered cubic (fcc) lattice, is proposed for molecular dynamics (MD) simulations. For the validity test of MD simulation, the Tersoff potential, the Stillinger and Weber (SW) potential, the environment-dependent interatomic (EDI) potential, the charge optimized many-body (COMB) potential, and the modified embedded-atom (MEAM) potential have been also employed for comparison. The crystal lattice of simulated silicon system is identified by calculating the distribution functions of the distances between the atoms and the angles between the lines linking an atom with its nearest neighbors. The results are also compared with the perfect silicon crystal. The crystal lattice, the crystallization temperature, and elastic constants have been calculated from MD simulations using above potentials. The results show that the systems with modified Tersoff, SW, EDI, COMB, and MEAM potentials could not exhibit the diamond structure and only the DL potential gives diamond lattice. The ground state for DL potential is the wurtzite structure, and the metastable state formed during rapid cooling is the cubic diamond structure. The physical parameters obtained from the simulation with DL potential are in agreement with the experiment results. This work indicated that only DL potential is valid for MD simulation of silicon crystal among above various potentials.
Materials Science,Chemical Physics
What problem does this paper attempt to address?
The problem this paper attempts to address is: how to reproduce the diamond structure of silicon through molecular dynamics (MD) simulations. Although various existing interatomic potential functions (such as Tersoff, Stillinger and Weber (SW), Environment-Dependent Interatomic (EDI), Charge Optimized Many-Body (COMB), and Modified Embedded Atom Method (MEAM) potentials) match well with experimental data in terms of physical parameters, they fail to form the diamond structure of silicon in MD simulations. Therefore, the authors propose a new Dual-Lattice (DL) potential function, aiming to reproduce the diamond structure of silicon through MD simulations and to verify its effectiveness and accuracy. Specifically, the goals of the paper include: 1. Constructing a new DL potential function that can reproduce the diamond structure of silicon through MD simulations. 2. Comparing the performance of the DL potential function with other existing potential functions (such as Tersoff, SW, EDI, COMB, and MEAM) in MD simulations to verify the effectiveness of the DL potential function. 3. Calculating and comparing the crystal structure, crystallization temperature, lattice constant, and elastic constants under different potential functions to evaluate the accuracy and reliability of the DL potential function in simulating silicon crystals. Through these objectives, the paper aims to provide a more accurate and effective potential function for atomic-scale simulations of silicon, thereby promoting the application and development of silicon materials in semiconductor technology.