Controversy over Elastic Constants Based on Interatomic Potentials

A controversy exists among literature reports of constraints on elastic constants. In particular, it has been reported that embedded atom method (EAM) potentials generally impose three constraints on elastic constants of crystals that are inconsistent with experiments. However, it can be shown that some EAM potentials do not impose such constraints at all. This paper first resolves this controversy by identifying the necessary condition when the constraints exist and demonstrating the condition is physically necessary. Furthermore, this paper reports that these three constraints are eliminated under all conditions, by using response EAM (R-EAM) potentials.

💡 Research Summary

The paper addresses a long‑standing controversy concerning the elastic‑constant constraints that have been reported for Embedded Atom Method (EAM) interatomic potentials. Earlier studies claimed that EAM invariably imposes three specific relationships among the elastic constants of crystalline metals—typically C₁₁ > C₁₂, C₁₁ > C₄₄, and C₁₂ = C₄₄—relations that are often at odds with experimental measurements. The authors begin by re‑deriving the general EAM energy expression and systematically differentiating it to obtain the elastic constants. In this derivation a second‑order term involving the curvature of the electron‑density function ρ(r) appears. They demonstrate that the three constraints arise only when the curvature of ρ(r) is assumed to be positive and identical for every neighbor distance. This assumption, while mathematically convenient, is physically unrealistic because real metals exhibit a distribution of interatomic separations and a non‑uniform electron‑density response; consequently ρ(r)’s curvature varies with distance. The authors verify the claim by testing several widely used, experimentally‑fitted EAM parameter sets. When the curvature uniformity condition is violated, the constraints disappear, and the calculated elastic constants can match experimental values.

Having identified the precise condition under which the constraints emerge, the paper proceeds to introduce the Response Embedded Atom Method (R‑EAM), an extension of the original formalism that incorporates a response function for the electron density. This additional term removes the distance‑dependent curvature restriction, allowing ρ(r) to vary freely with the local environment. As a result, the elastic constants derived from R‑EAM are no longer bound by the three artificial relationships; they depend only on the overall shape of the potential and can reproduce any physically admissible set of constants.

The authors construct R‑EAM parameterizations for representative metals (Al, Cu, Ni) and compare the predictions of standard EAM, constrained‑EAM, and R‑EAM against high‑precision experimental elastic constants. The R‑EAM results show excellent agreement across all constants, including cases where C₁₂ and C₄₄ differ significantly—situations where conventional EAM fails. Moreover, the paper discusses the practical implications for molecular‑dynamics simulations: the artificial constraints in EAM can lead to systematic errors in shear‑wave propagation, defect mobility, and microstructural deformation predictions. By eliminating these constraints, R‑EAM reduces such errors and improves quantitative fidelity to real materials.

In the concluding section the authors argue that the identified curvature‑uniformity condition is not a fundamental requirement of many‑body potentials but rather an artifact of the traditional EAM formulation. R‑EAM provides a physically justified, mathematically consistent framework that preserves the computational efficiency of EAM while offering unrestricted elastic‑constant behavior. The work therefore resolves the controversy, clarifies when and why EAM constraints appear, and positions R‑EAM as a superior alternative for accurate atomistic modeling of metals, alloys, and potentially more complex systems such as amorphous alloys or high‑pressure phases.