First-principles calculation of higher-order elastic constants from divided differences

Abstract

A method is presented to calculate from first principles the higher-order elastic constants of a solid material. The method relies on finite strain deformations, a density functional theory approach to calculate the Cauchy stress tensor, and a recursive numerical differentiation technique homologous to the divided differences polynomial interpolation algorithm. The method is applicable as is to any material, regardless its symmetry, to calculate elastic constants of, in principle, any order. Here, we introduce conceptual framework and technical details of our method, we discuss sources of errors, we assess convergence trends, and we present selected applications. In particular, our method is used to calculate elastic constants up to the 6${}^{th}$ order of two crystalline materials with the cubic symmetry, silicon and gold. To demonstrate general applicability, our method is also used to calculate the elastic constants up to the 5${}^{th}$ order of α-quartz, a crystalline material belonging to the trigonal crystal system, and the second- and third-order elastic constants of kevlar, a material with an anisotropic bonding network. Higher order elastic constants computed with our method are validated against density functional theory calculations by comparing stress responses to large deformations derived within the continuum approximation.