Efficacious symmetry-adapted atomic displacement method for lattice dynamical studies

Abstract

Small displacement methods have been successfully used to calculate the lattice dynamical properties of crystals. It involves displacing atoms by a small amount in order to calculate the induced forces on all atoms in a supercell for the computation of force constants. Even though these methods are widely in use, to our knowledge, there is no systematic discussion of optimal displacement directions from the crystal's symmetry point of view nor a rigorous error analysis of such methods. Based on the group theory and point group symmetry of a crystal, we propose displacement directions, with an equivalent concept of the group of $k$, deduced directly in the Cartesian coordinates rather than the usual fractional coordinates, that maintain the theoretical maximum for the triple product $V$ spanned by the three displacements to avoid possible severe roundoff errors. The proposed displacement directions are generated from a minimal set of irreducible atomic displacements that keep the required independent force calculations to a minimum. We find the error in the calculated force constant explicitly depends on the inverse of $V$ and inaccuracy of the forces. Test systems such as Si, graphene, and orthorhombic Sb2S3 are used to illustrate the method. Our displacement method is shown to be very robust in treating low-symmetry cells with a large `aspect ratio' due to huge differences in lattice parameters, use of a large vacuum height, or a very oblique unit cell due to unconventional choice of primitive lattice vectors. It is expected that our displacement strategy can be used to address higher-order interatomic interactions to achieve good accuracy and efficiency.

arXiv comment10 pages, 2 tables, 6 figures

References in corpus (20)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Phonons and related properties of extended systems from density-functional perturbation theory
• High-throughput electronic band structure calculations: challenges and tools
• Electron-phonon interactions from first principles
• High-accuracy first-principles determination of the structural, vibrational and thermodynamical properties of diamond, graphite, and derivatives
• The Effect of Lattice Vibrations on Substitutional Alloy Thermodynamics
• Interlayer breathing and shear modes in few-trilayer MoS2 and WSe2
• AiiDA: Automated Interactive Infrastructure and Database for Computational Science
• High-Throughput Computational Screening of thermal conductivity, Debye temperature and Grüneisen parameter using a quasi-harmonic Debye Model
• Lattice dynamics and electron-phonon coupling calculations using non-diagonal supercells
• First-principles study of graphene edge properties and flake shapes
• First-principles study of the lattice dynamics of Sb2S3
• Direct calculation of the linear thermal expansion coefficients of MoS2 via symmetry-preserving deformations
• Theory of thermal expansion: Quasi-harmonic approximation and corrections from quasiparticle renormalization
• Anharmonic interatomic force constants and thermal conductivity from Grüneisen parameters: an application to graphene
• Large anharmonic effect and thermal expansion anisotropy of metal chalcogenides: The case of antimony sulfide
• Origins of anisotropic thermal expansion in flexible materials
• Anisotropic thermal expansion of bismuth from first principles
• A group theoretical approach to computing phonons and their interactions
• Anharmonic phonon effects on linear thermal expansion of trigonal bismuth selenide and antimony telluride crystals

Cited by in corpus (3)

• A first-principles investigation of the linear thermal expansion coefficients of BeF${}_2$: Giant thermal expansion
• Complementary local-global approach for phonon mode connectivities
• A size-consistent Grüneisen-quasiharmonic approach for lattice thermal conductivity