中国物理B ›› 2021, Vol. 30 ›› Issue (11): 114301-114301.doi: 10.1088/1674-1056/abf4c0

• • 上一篇    下一篇

Phonon dispersion relations of crystalline solids based on LAMMPS package

Zhiyong Wei(魏志勇)1,†, Tianhang Qi(戚天航)1, Weiyu Chen(陈伟宇)2, and Yunfei Chen(陈云飞)1   

  1. 1 Jiangsu Key Laboratory for Design & Manufacture of Micro/Nano Biomedical Instruments and School of Mechanical Engineering, Southeast University, Nanjing 211189, China;
    2 College of Mechanical and Electronic Engineering, Nanjing Forestry University, Nanjing 210037, China
  • 收稿日期:2021-02-23 修回日期:2021-03-20 接受日期:2021-04-05 出版日期:2021-10-13 发布日期:2021-10-13
  • 通讯作者: Zhiyong Wei E-mail:zywei@seu.edu.cn
  • 基金资助:
    Project supported by the National Key R&D Program of China (Grant No. 2017YFB0406000) and the Southeast University “Zhongying Young Scholars” Project.

Phonon dispersion relations of crystalline solids based on LAMMPS package

Zhiyong Wei(魏志勇)1,†, Tianhang Qi(戚天航)1, Weiyu Chen(陈伟宇)2, and Yunfei Chen(陈云飞)1   

  1. 1 Jiangsu Key Laboratory for Design & Manufacture of Micro/Nano Biomedical Instruments and School of Mechanical Engineering, Southeast University, Nanjing 211189, China;
    2 College of Mechanical and Electronic Engineering, Nanjing Forestry University, Nanjing 210037, China
  • Received:2021-02-23 Revised:2021-03-20 Accepted:2021-04-05 Online:2021-10-13 Published:2021-10-13
  • Contact: Zhiyong Wei E-mail:zywei@seu.edu.cn
  • Supported by:
    Project supported by the National Key R&D Program of China (Grant No. 2017YFB0406000) and the Southeast University “Zhongying Young Scholars” Project.

摘要: The phonon dispersion relations of crystalline solids play an important role in determining the mechanical and thermal properties of materials. The phonon dispersion relation, as well as the vibrational density of states, is also often used as an indicator of variation of lattice thermal conductivity with the external stress, defects, etc. In this study, a simple and fast tool is proposed to acquire the phonon dispersion relation of crystalline solids based on the LAMMPS package. The theoretical details for the calculation of the phonon dispersion relation are derived mathematically and the computational flow chart is present. The tool is first used to calculate the phonon dispersion relation of graphene with two atoms in the unit cell. Then, the phonon dispersions corresponding to several potentials or force fields, which are commonly used in the LAMMPS package to modeling the graphene, are obtained to compare with that from the DFT calculation. They are further extended to evaluate the accuracy of the used potentials before the molecular dynamics simulation. The tool is also used to calculate the phonon dispersion relation of superlattice structures that contains more than one hundred of atoms in the unit cell, which predicts the phonon band gaps along the cross-plane direction. Since the phonon dispersion relation plays an important role in the physical properties of condensed matter, the proposed tool for the calculation of the phonon dispersion relation is of great significance for predicting and explaining the mechanical and thermal properties of crystalline solids.

关键词: phonon dispersion relation, molecular dynamics, force constant, potential function

Abstract: The phonon dispersion relations of crystalline solids play an important role in determining the mechanical and thermal properties of materials. The phonon dispersion relation, as well as the vibrational density of states, is also often used as an indicator of variation of lattice thermal conductivity with the external stress, defects, etc. In this study, a simple and fast tool is proposed to acquire the phonon dispersion relation of crystalline solids based on the LAMMPS package. The theoretical details for the calculation of the phonon dispersion relation are derived mathematically and the computational flow chart is present. The tool is first used to calculate the phonon dispersion relation of graphene with two atoms in the unit cell. Then, the phonon dispersions corresponding to several potentials or force fields, which are commonly used in the LAMMPS package to modeling the graphene, are obtained to compare with that from the DFT calculation. They are further extended to evaluate the accuracy of the used potentials before the molecular dynamics simulation. The tool is also used to calculate the phonon dispersion relation of superlattice structures that contains more than one hundred of atoms in the unit cell, which predicts the phonon band gaps along the cross-plane direction. Since the phonon dispersion relation plays an important role in the physical properties of condensed matter, the proposed tool for the calculation of the phonon dispersion relation is of great significance for predicting and explaining the mechanical and thermal properties of crystalline solids.

Key words: phonon dispersion relation, molecular dynamics, force constant, potential function

中图分类号:  (Ultrasonics, quantum acoustics, and physical effects of sound)

  • 43.35.+d
62.60.+v (Acoustical properties of liquids) 63.20.-e (Phonons in crystal lattices)