Elastic Constants and Bending Rigidities from Long-Wavelength Perturbation Expansions

Lin, Changpeng;Poncé, Samuel;Macheda, Francesco;Mauri, Francesco;Marzari, Nicola
(2026) PRX Energy — Vol. 5, n° 1, p. 13012 (2026)

Files

Lin2026.pdf
  • Open Access
  • Adobe PDF
  • 1.13 MB

Details

Authors
  • Lin, Changpengorcid-logoÉcole Polytechnique Fédérale de Lausanne
    Author
  • Poncé, Samuelorcid-logoUniversité catholique de Louvain
    Author
  • Macheda, Francescoorcid-logoUniversità di Roma La Sapienza
    Author
  • Mauri, FrancescoUniversità di Roma La Sapienza
    Author
  • Marzari, Nicolaorcid-logoÉcole Polytechnique Fédérale de Lausanne
    Author
Abstract
Mechanical and elastic properties of materials are among the most fundamental quantities for many engineering and industrial applications. Here, we present a formulation that is efficient and accurate for calculating the elastic and bending rigidity tensors of crystalline solids, leveraging interatomic force constants and long-wavelength perturbation theory. Crucially, in the long-wavelength limit, lattice vibrations induce macroscopic electric fields, which further couple with the propagation of elastic waves, and a separate treatment on the long-range electrostatic interactions is thereby required to obtain elastic properties under the appropriate electrical boundary conditions. A cluster expansion of the charge-density response and dielectric screening function in the long-wavelength limit has been developed to efficiently extract multipole and dielectric tensors of arbitrarily high order. We implement the proposed method in a first-principles framework and perform extensive validations on silicon, NaCl, GaAs and rhombohe-dral BaTiO3 as well as monolayer graphene, hexagonal BN, MoS2 , and InSe, obtaining good to excellent agreement with other theoretical approaches and experimental measurements. Notably, we establish that multipolar interactions up to at least octupoles are necessary to obtain the accurate short-circuit elastic tensor of bulk materials, while higher orders beyond octupole interactions are required to converge the bending rigidity tensor of two-dimensional crystals. The present approach greatly simplifies the calculations of bending rigidities and will enable the automated characterization of the mechanical properties of novel functional materials.
Affiliations

Citations

Lin, C., Poncé, S., Macheda, F., Mauri, F., & Marzari, N. (2026). Elastic Constants and Bending Rigidities from Long-Wavelength Perturbation Expansions. PRX Energy, 5(1), 13012. https://doi.org/10.1103/hc53-g1p3 (Original work published 2026)