Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

The tensorial form of Hooke's law for the strain-stress relationship in a crystal is (in the Voigt notation):

Hooke's law

where $\sigma$ is the strain, $\epsilon$ is the stress and C is the stiffness tensor:

stiffness tensor

For a crystalline system of the cubic symmetry class, the stiffness tensor reduces to:

enter image description here

The Born criterion for the stability of an unstrained crystal is that free energy must be represented by a positive defined quadratic form. In the case of a cubic crystal, it is known that this is equivalent to the following three conditions on the elastic constants:

$$C_{11} - C_{12} > 0$$ $$ C_{44} > 0$$ $$ C_{11} + 2 C_{12} > 0$$

But what about lower symmetry classes? What is the generic Born criterion for stability of a crystal? I have quite convinced myself that all the eigenvalues of $C$ must be positive, but I cannot find confirmation of that anywhere. Is it right? Is there a reference on that topic?

share|cite|improve this question

This derivation is correct for cubic crystals without external pressure.

A fresh review can be found in Rev. Mod. Phys. 84, 945 (2012) Lattice instabilities in metallic elements

share|cite|improve this answer

While I can't find a reference for this, I think your criterion is correct. Here's the argument I would use - if it's the same as yours, then maybe it's right! The elastic energy is ( $$E = \frac{1}{2}\int d^d r \epsilon_i C_{ij} \epsilon_j$$ where $C_{ij}$ is symmetric. Because $C_{ij}$ is symmetric and real, it can be diagonalized, and its eigenvectors are complete and orthogonal ( Then you can expand $\epsilon_j$ in terms of the eigenvectors of $C_{ij}$, $\bf{v}^{(k)}$, i.e. $C_{ij} v_j^{(k)} = \lambda_k v_i^{(k)}$. $$\epsilon_j = \sum_k (\bf{v}^{(k)} \cdot \bf{\epsilon}) v^{(k)}_j$$ (This assumes that $\bf{v}^{(k)}$ is orthonormal, which we can always choose without loss of generality.) We then find that $$E = \frac{1}{2} \int d^d r \lambda_k (\bf{v}^{(k)} \cdot \bf{\epsilon})^2$$.
Stability means that there are no modes that will lead to an unbounded decrease in energy, and no marginal modes - i.e. that $\lambda_k > 0$ for all $k$ - your positive definite quadratic form.

share|cite|improve this answer
Yes, that's how I figured it too… – F'x Apr 2 '13 at 7:43
up vote 1 down vote accepted

I've found a good analysis of the stability conditions for a crystal's elastic constants, both unstrained and under stress, in:

J. W. Morris Jr and C. R. Krenn, Philos. Mag. A 2000, 12, 2827–2840

To quote them:

In the linear elastic limit the conditions of internal stability reduce to the familiar condition that the 6 x 6 matrix $C_{ij}$ of elastic moduli have no negative eigenvalues.

So, almost two years later, after realizing a lot of people had trouble with this question (and there are mistakes made in some of the literature on the topic), we have published a short pedagogical reference on the issue: F. Mouhat and F.-X. Coudert, Phys. Rev. B, 90, 224104 (2014) (also on arXiv).

share|cite|improve this answer
R. Cowley Phys Rev B 1976 is a great article as well for this subject. He covers many of the crystallographic point groups and shows their stability criterion. – John M Oct 24 '14 at 3:01

A detailed analysis may be found also in arXiv:1104.0173 [astro-ph.SR], D. A. Baiko "Shear modulus of neutron star crust", Eqs. (29)-(33).

share|cite|improve this answer
That paper is, at best, only very marginally related to my question :( – F'x Apr 2 '13 at 7:43

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.