First-Principles Calculation of Third-Order Elastic Constants via Numerical Differentiation of the Second
Tengfei Cao1,2, David Cuffari1,3, Angelo Bongiorno1,2,3,4
1Department of Chemistry, College of Staten Island, Staten Island, New York 10314, USA.
Abstract:
A general method is presented to calculate from first principles the full set of third-order elastic constants of a material of arbitrary symmetry. The method here illustrated relies on a plane-wave density functional theory scheme to calculate the Cauchy stress and the numerical differentiation of the second Piola-Kirchhoff stress tensor to evaluate the elastic constants. It is shown that finite difference formulas lead to a cancellation of the finite basis set errors, whereas simple solutions are proposed to eliminate numerical errors arising from the use of Fourier interpolation techniques. Applications to diamond, silicon, aluminum, magnesium, graphene, and a graphane conformer give results in excellent agreement with both experiments and previous calculations based on fitting energy density curves, demonstrating both the accuracy and generality of our new methodology to investigate nonlinear elastic behaviors of materials.
Related Concept Videos
01:24Numerical Calculations
The solution to a problem is obtained using different methods. While manually solving algebraic symbols is one of the most common methods, the graphical method is often preferred. Computers...
02:46Calculating the Equilibrium Constant
For example, gaseous nitrogen dioxide forms dinitrogen tetroxide according to this equation:
01:22Elastic Strain Energy for Normal Stresses
If...
01:20Elastic Strain Energy for Shearing Stresses
01:21Kirchhoff's Rules
Kirchhoff's first rule is called the junction rule. A junction, also known as a node, is a connection of three or more wires. The rule states that the sum of all currents entering a junction must equal the sum of all currents leaving the junction.
02:27Le Chatelier's Principle: Changing Concentration