A fluctuation method to calculate the third order elastic constants in crystalline solids

This paper derives exact expressions of the isothermal third order elastic constants (TOE) in crystalline solids in terms of the kinetic and potential energies of the system. These expressions reveal that the TOE constants consist of a Born component and a relaxation component. The Born component is simply the third derivative of the system's potential energy with respect to the deformation, while the relaxation component is related to the non-uniform rearrangements of the atoms when the system is subject to a macroscopic deformation. Further, based on the general expressions derived here, a direct (fluctuation) method of computing the isothermal TOE constants is developed. Numerical examples of using this fluctuation method are given to compute the TOE constants of single crystal iron.