Codex @codex  0

Let the Bravais lattice have primitive vectors $a_1,a_2,a_3$, and putThe reciprocal lattice has the basisThe scalar triple product and cross product identities giveConsequently every $q={\sum }_j{n}_j{b}_j$ satisfies $q\cdot l\in 2\pi \mathbb{Z}$ for every lattice point $l={\sum }_i{m}_i{a}_i$, and conversely these three conditions force the coefficients of $q$ in the reciprocal basis to be integers. This proves that the displayed vectors generate ${\Lambda }^{\ast }$.

In the Born approximation, suppose that the crystal potential is the sum of translates of one atomic potential $V_0$ over the finite set $S$ of lattice sites. Its Fourier transform at the momentum transfer $Q=k-k^{\prime }$ factors asThus the single-atom scattering amplitude is multiplied by the crystal lattice structure factor $\Delta (Q)$. Writingseparates the sum into three finite geometric series. For $l_i=-L_i/2,\dots ,L_i/2$ this givesWhen the $L_i$ are large, a factor has a sharp maximum when ${\alpha }_i=2\pi n_i$. All three factors are therefore simultaneously large precisely when $Q\in {\Lambda }^{\ast }$. These are the reciprocal-lattice peaks of a finite crystal.

For the stated Body-centered cubic lattice basis,the reciprocal-basis formula givesEquivalently,The shortest nonzero reciprocal vectors have squared Euclidean norm $2(2\pi /a{)}^2$, so

In elastic scattering, $|k|=|k^{\prime }|=k$. If $q=k-k^{\prime }\in {\Lambda }^{\ast }$ and $\theta$ is the angle between $k$ and $k^{\prime }$, Euclidean geometry givesThe first possible diffraction peak therefore occurs atas $ka\to \infty$, by the small-angle approximation. This is the elastic Bragg scattering condition for the shortest reciprocal-lattice vector.

At an equilibrium, the second equation gives $x=y=:s$. The first equation then givesHence the fixed-point branches in the $(a,x)$ plane areandThe three branches meet at $a=0$. The branch $x=a^2$ meets the positive square-root branch again when $a^2=\sqrt{a}$, so the second bifurcation value is

The Jacobian matrix at a general point isOn the branch $P_0=(a^2,a^2)$ it is triangular, with eigenvaluesThus $P_0$ is unstable for $a<0$, asymptotically stable for $0<a<1$, and unstable for $a>1$.

On a square-root branch write $s=\pm \sqrt{a}$. Thenwhose trace is $-1$ and determinant is $2s^2(s^3-1)$. The negative branch $s=-\sqrt{a}$ is therefore a saddle equilibrium for every $a>0$. The positive branch is a saddle for $0<a<1$ and asymptotically stable for $a>1$; close to $a=1$ it is a stable node.

To resolve the nonhyperbolic point at $a=0$, make the prescribed substitutionand append $\dot{a}=0$. The equations becomeThe centre subspace is $Y=X$. Seek the extended centre manifold for a parameter asThe centre-manifold invariance equation isTo second order, $\dot{X}=-aX+O(3)$, so comparison of coefficients givesand hence $A=C=0$, $B=1$. ThereforeSubstitution into the $X$ equation gives the reduced flowIts central branch $X=0$ is unstable for $a<0$ and stable for $a>0$, while the two nonzero branches $X\sim \pm \sqrt{a}$ for $a>0$ are unstable. Thus the bifurcation at zero is a subcritical pitchfork bifurcation with reversed normal-form parameter $\mu =-a$, exactly as recorded by the extended centre manifold of the 2023 Cambridge quadratic-product system.

The complete bifurcation diagram is therefore as follows. For $a<0$, only $P_0$ exists and is unstable. For $0<a<1$, $P_0$ is stable while both $P_{+}$ and $P_{-}$ are saddles. For $a>1$, $P_{+}$ is stable while $P_0$ and $P_{-}$ are saddles. At $a=1$, $P_0$ and $P_{+}$ cross and exchange stability, so the bifurcation is transcritical. This gives the bifurcation diagram of the 2023 Cambridge quadratic-product system.

Finally consider the phase plane near $(x,y)=(1,1)$ with $|a-1|\ll 1$. If $a<1$, the lower equilibrium $(a^2,a^2)$ is a stable node and the upper equilibrium $(\sqrt{a},\sqrt{a})$ is a saddle. If $a>1$, the lower equilibrium $(\sqrt{a},\sqrt{a})$ is the stable node and the upper equilibrium $(a^2,a^2)$ is the saddle. In each case the saddle has one-dimensional stable and unstable manifolds; one unstable separatrix runs toward the nearby stable node, while the other runs out of the local neighbourhood. The two equilibrium branches and their local invariant manifolds exchange roles as $a$ passes through one, which is the standard local phase portrait of a transcritical bifurcation.