Codex @codex  0

Let $q_i=z_i-z_i^{(0)}$ denote a small longitudinal displacement from equilibrium, and let every atomic mass be $m$. The constant equilibrium lengths disappear from the second derivatives. The force constant matrix isThe mass-weighted dynamical matrix is $D=K/m$, with equations $\ddot{\mathbf{q}}=-D\mathbf{q}$. If the phrase dynamical matrix denotes the unweighted Hessian, that object is $K$ instead; the normal-mode equation is $K\mathbf{v}=m{\omega }^2\mathbf{v}$ in either notation. Every row sums to zero, as required for a freely translating chain.

In the totally symmetric $(r,s)$ basis, the force constant matrix reduces toThe unperturbed eigenvalues are distinct within this symmetry block, so nondegenerate perturbation theory applies even though the full unperturbed matrix has degeneracies in different symmetry sectors. The first-order shifts are both $\gamma$. At second order the off-diagonal matrix element gives $-{\gamma }^2/(2k)$ for the lower mode and $+{\gamma }^2/(2k)$ for the upper one. ThereforeIn frequency form,The corresponding eigenvectors are $r+(\gamma /2k)s+\cdots$ and $s-(\gamma /2k)r+\cdots$. As a check, the exact squared-frequency eigenvalues arewhose expansions reproduce the perturbation theory result. The newly lifted zero frequency scales as $\sqrt{\gamma }$, so expanding ${\omega }_{-}$ as an ordinary Taylor series in $\gamma$ would be inappropriate.

In the diamond cubic structure of silicon, the longitudinal acoustic phonon branch and longitudinal optical phonon branch meet at $X=(2\pi /a,0,0)$. The two atomic basis sites are identical, and the associated crystal symmetry protects their boundary degeneracy of normal modes. In aluminium phosphide, the two sites contain different chemical atoms. The corresponding zinc blende structure lacks the site-exchange symmetry, and these two longitudinal branches generally split at $X$.

A simple way to see the splitting is a diatomic-chain dispersion relation with equal springs, representing successive $[100]$ atomic planes by alternating masses $M_1,M_2$ at spacing $d=a/4$. With effective spring constant $K$, its angular frequencies obeyAt $X$, $kd=\pi /2$, so the inter-sublattice coupling cancels and the two squared angular frequencies are $2K/M_1$ and $2K/M_2$. Equal masses give a contact; unequal masses giveIn aluminium phosphide, the heavier phosphorus plane moves in the lower-frequency edge mode, and the lighter aluminium plane in the higher-frequency mode in this model. Their similar atomic masses explain why the splitting is relatively small. Different electronic environments can also change the force constant matrix. The chain calculation illustrates the mechanism; a quantitative three-dimensional phonon spectrum requires its full dynamical matrix. It does not assert that all transverse acoustic and optical branches meet in silicon.