Abstract

The foundation of modern theory of lattice dynamics rests on the separation between short-range and long-range interatomic force constants, where the latter are associated with macroscopic electric fields acting in a neighborhood of the Brillouin zone center. In three-dimensional (3D) insulators, the famous dipole-dipole formula was established long ago by Cochran and Cowley; the subsequent work of Pick, Cohen and Martin provided a formal derivation in the framework of first-principles theory. Recently, we have incorporated higher-order interactions involving, e.g., dynamical quadrupoles,[1] and demonstrated their role in the accurate interpolation of phonon band structures.[2] In this talk, I shall retrace an analogous journey in the context of two-dimensional (2D) crystals, by presenting a rigorous derivation of the long-range screening and interatomic forces in 2D.[3] This enables a systematic generalization of the existing formulas [4] to an arbitrary multipolar order. In particular, I will discuss how to incorporate out-of-plane dipoles and dynamical quadrupoles in the long-range part of the dynamical matrix, and how to achieve an optimal representation of the dielectric function. Numerical tests on monolayer BN, SnS2 and BaTiO3 membranes demonstrate a drastic improvement in the description of the long-range electrostatic interactions, with comparable benefits to the quality of the interpolated phonon band structure.[3]

[1] M. Royo and M. Stengel, Physical Review X 9 , 021050 (2019).

[2] M. Royo, K. Hahn and M. Stengel, Physical Review Letters 125, 217602 (2020).

[3] M. Royo and M. Stengel, Physical Review X 11 , 041027 (2021).

[4] T. Sohier et al., Nano Letters 6, 3758 (2017).