Physics:Laue equations

From HandWiki
Short description: Equations describing diffraction in a crystal lattice
Laue equation

In crystallography and solid state physics, the Laue equations relate incoming waves to outgoing waves in the process of elastic scattering, where the photon energy or light temporal frequency does not change upon scattering by a crystal lattice. They are named after physicist Max von Laue (1879–1960).

The Laue equations can be written as [math]\displaystyle{ \mathbf{\Delta k}= \mathbf{k}_{\mathrm{out}} - \mathbf{k}_{\mathrm{in}} = \mathbf{G} }[/math] as the condition of elastic wave scattering by a crystal lattice, where [math]\displaystyle{ \mathbf{\Delta k} }[/math] is the scattering vector, [math]\displaystyle{ \mathbf{k}_{\mathrm{in}} }[/math], [math]\displaystyle{ \mathbf k_{\mathrm{out}} }[/math] are an incoming and outgoing wavevectors (to the crystal and from the crystal, by scattering), and [math]\displaystyle{ \mathbf{G} }[/math] is a crystal reciprocal lattice vector. Due to elastic scattering [math]\displaystyle{ |\mathbf{k}_{\mathrm{out}}|^2=|\mathbf{k}_{\mathrm{in}}|^2 }[/math], three vectors. [math]\displaystyle{ \mathbf{G} }[/math], [math]\displaystyle{ \mathbf{k}_{\mathrm{out}} }[/math], and [math]\displaystyle{ -\mathbf{k}_{\mathrm{in}} }[/math] , form a rhombus if the equation is satisfied. If the scattering satisfies this equation, all the crystal lattice points scatter the incoming wave toward the scattering direction (the direction along [math]\displaystyle{ \mathbf k_{\mathrm{out}} }[/math]). If the equation is not satisfied, then for any scattering direction, only some lattice points scatter the incoming wave. (This physical interpretation of the equation is based on the assumption that scattering at a lattice point is made in a way that the scattering wave and the incoming wave have the same phase at the point.) It also can be seen as the conservation of momentum as [math]\displaystyle{ \hbar \mathbf{k}_{\mathrm{out}} = \hbar \mathbf{k}_{\mathrm{in}} + \hbar \mathbf{G} }[/math] since [math]\displaystyle{ \mathbf{G} }[/math] is the wavevector for a plane wave associated with parallel crystal lattice planes. (Wavefronts of the plane wave are coincident with these lattice planes.)

The equations are equivalent to Bragg's law; the Laue equations are vector equations while Bragg's law is in a form that is easier to solve, but these tell the same content.

The Laue equations

Let [math]\displaystyle{ \mathbf{a}\, ,\mathbf{b}\, ,\mathbf{c} }[/math] be primitive translation vectors (shortly called primitive vectors) of a crystal lattice [math]\displaystyle{ L }[/math], where atoms are located at lattice points described by [math]\displaystyle{ \mathbf x = p\,\mathbf a+q\,\mathbf b+r\,\mathbf c }[/math] with [math]\displaystyle{ p }[/math], [math]\displaystyle{ q }[/math], and [math]\displaystyle{ r }[/math] as any integers. (So [math]\displaystyle{ \mathbf x }[/math] indicating each lattice point is an integer linear combination of the primitive vectors.)

Let [math]\displaystyle{ \mathbf{k}_{\mathrm{in}} }[/math] be the wavevector of an incoming (incident) beam or wave toward the crystal lattice [math]\displaystyle{ L }[/math], and let [math]\displaystyle{ \mathbf k_{\mathrm{out}} }[/math] be the wavevector of an outgoing (diffracted) beam or wave from [math]\displaystyle{ L }[/math]. Then the vector [math]\displaystyle{ \mathbf k_{\mathrm{out}} - \mathbf k_{\mathrm{in}} = \mathbf{\Delta k} }[/math], called the scattering vector or transferred wavevector, measures the difference between the incoming and outgoing wavevectors.

The three conditions that the scattering vector [math]\displaystyle{ \mathbf{\Delta k} }[/math] must satisfy, called the Laue equations, are the following:

[math]\displaystyle{ \mathbf{\Delta k}\cdot \mathbf{a} =2\pi h }[/math]
[math]\displaystyle{ \mathbf{\Delta k}\cdot \mathbf{b} =2\pi k }[/math]
[math]\displaystyle{ \mathbf{\Delta k}\cdot \mathbf{c} =2\pi l }[/math]

where numbers [math]\displaystyle{ h, k, l }[/math] are integer numbers. Each choice of integers [math]\displaystyle{ (h,k,l) }[/math], called Miller indices, determines a scattering vector [math]\displaystyle{ \mathbf{\Delta k} }[/math]. Hence there are infinitely many scattering vectors that satisfy the Laue equations as there are infinitely many choices of Miller indices [math]\displaystyle{ (h,k,l) }[/math]. Allowed scattering vectors [math]\displaystyle{ \mathbf{\Delta k} }[/math] form a lattice [math]\displaystyle{ L^* }[/math], called the reciprocal lattice of the crystal lattice [math]\displaystyle{ L }[/math], as each [math]\displaystyle{ \mathbf{\Delta k} }[/math] indicates a point of [math]\displaystyle{ L^* }[/math]. (This is the meaning of the Laue equations as shown below.) This condition allows a single incident beam to be diffracted in infinitely many directions. However, the beams corresponding to high Miller indices are very weak and can't be observed. These equations are enough to find a basis of the reciprocal lattice (since each observed [math]\displaystyle{ \mathbf{\Delta k} }[/math] indicates a point of the reciprocal lattice of the crystal under the measurement), from which the crystal lattice can be determined. This is the principle of x-ray crystallography.

Mathematical derivation

For an incident plane wave at a single frequency [math]\displaystyle{ \displaystyle f }[/math] (and the angular frequency [math]\displaystyle{ \displaystyle \omega =2\pi f }[/math]) on a crystal, the diffracted waves from the crystal can be thought as the sum of outgoing plane waves from the crystal. (In fact, any wave can be represented as the sum of plane waves, see Fourier Optics.) The incident wave and one of plane waves of the diffracted wave are represented as

[math]\displaystyle{ \displaystyle f_{\mathrm {in} }(t,\mathbf {x} )=A_{\mathrm {in} }\cos(\omega \,t-\mathbf {k} _{\mathrm {in} }\cdot \mathbf {x} +\varphi _{\mathrm {in} }), }[/math]
[math]\displaystyle{ \displaystyle f_{\mathrm {out} }(t,\mathbf {x} )=A_{\mathrm {out} }\cos(\omega \,t-\mathbf {k} _{\mathrm {out} }\cdot \mathbf {x}+\varphi _{\mathrm {out} }), }[/math]

where [math]\displaystyle{ \displaystyle \mathbf {k} _{\mathrm {in} } }[/math] and [math]\displaystyle{ \displaystyle \mathbf {k} _{\mathrm {out} } }[/math] are wave vectors for the incident and outgoing plane waves, [math]\displaystyle{ \displaystyle \mathbf {x} }[/math] is the position vector, and [math]\displaystyle{ \displaystyle t }[/math] is a scalar representing time, and [math]\displaystyle{ \varphi _{\mathrm {in} } }[/math] and [math]\displaystyle{ \varphi _{\mathrm {out} } }[/math] are initial phases for the waves. For simplicity we take waves as scalars here, even though the main case of interest is an electromagnetic field, which is a vector. We can think these scalar waves as components of vector waves along a certain axis (x, y, or z axis) of the Cartesian coordinate system.

The incident and diffracted waves propagate through space independently, except at points of the lattice [math]\displaystyle{ L }[/math] of the crystal, where they resonate with the oscillators, so the phases of these waves must coincide.[1] At each point [math]\displaystyle{ \mathbf x = p\,\mathbf a+q\,\mathbf b+r\,\mathbf c }[/math] of the lattice [math]\displaystyle{ L }[/math], we have

[math]\displaystyle{ \cos (\omega\,t-\mathbf k_{\mathrm{in}}\cdot\mathbf x+\varphi _{\mathrm {in} }) = \cos (\omega\,t-\mathbf k_{\mathrm{out}}\cdot\mathbf x+\varphi _{\mathrm {out} }), }[/math]

or equivalently, we must have

[math]\displaystyle{ \omega\,t-\mathbf k_{\mathrm{in}}\cdot\mathbf x + \varphi _{\mathrm {in} }= \omega\,t-\mathbf k_{\mathrm{out}}\cdot\mathbf x + \varphi _{\mathrm {out} } + 2\pi n, }[/math]

for some integer [math]\displaystyle{ n }[/math], that depends on the point [math]\displaystyle{ \mathbf x }[/math]. Since this equation holds at [math]\displaystyle{ \mathbf x=0 }[/math], [math]\displaystyle{ \varphi _{\mathrm {in} } = \varphi _{\mathrm {out} }+2\pi n' }[/math] at some integer [math]\displaystyle{ n' }[/math]. So

[math]\displaystyle{ \omega\,t-\mathbf k_{\mathrm{in}}\cdot\mathbf x = \omega\,t-\mathbf k_{\mathrm{out}}\cdot\mathbf x + 2\pi n. }[/math]

(We still use [math]\displaystyle{ n }[/math] instead of [math]\displaystyle{ (n-n') }[/math] since both the notations essentially indicate some integer.) By rearranging terms, we get

[math]\displaystyle{ \mathbf{\Delta k}\cdot \mathbf x = (\mathbf k_{\mathrm{out}}-\mathbf k_{\mathrm{in}})\cdot \mathbf x = 2\pi n. }[/math]

Now, it is enough to check that this condition is satisfied at the primitive vectors [math]\displaystyle{ \mathbf a,\mathbf b,\mathbf c }[/math] (which is exactly what the Laue equations say), because, at any lattice point [math]\displaystyle{ \mathbf x = p\,\mathbf a+q\,\mathbf b+r\,\mathbf c }[/math], we have

[math]\displaystyle{ \mathbf{\Delta k} \cdot \mathbf x = \mathbf{\Delta k}\cdot (p\,\mathbf a+q\,\mathbf b+r\,\mathbf c) = p(\mathbf{\Delta k}\cdot \mathbf a) + q(\mathbf{\Delta k}\cdot \mathbf b) + r(\mathbf{\Delta k}\cdot \mathbf c) = p\,(2\pi h) + q\,(2\pi k) + r\,(2\pi l) = 2\pi(ph+qk+rl)=2\pi n, }[/math]

where [math]\displaystyle{ n }[/math] is the integer [math]\displaystyle{ ph+qk+rl }[/math]. The claim that each parenthesis, e.g. [math]\displaystyle{ (\mathbf{\Delta k} \cdot \mathbf{a}) }[/math], is to be a multiple of [math]\displaystyle{ 2\pi }[/math] (that is each Laue equation) is justified since otherwise [math]\displaystyle{ p(\mathbf{\Delta k}\cdot \mathbf a) + q(\mathbf{\Delta k}\cdot \mathbf b) + r(\mathbf{\Delta k}\cdot \mathbf c) = 2\pi n }[/math] does not hold for any arbitrary integers [math]\displaystyle{ p, q, r }[/math].

This ensures that if the Laue equations are satisfied, then the incoming and outgoing (diffracted) wave have the same phase at each point of the crystal lattice, so the oscillations of atoms of the crystal, that follows the incoming wave, can at the same time generate the outgoing wave at the same phase of the incoming wave.

Relation to reciprocal lattices and Bragg's Law

If [math]\displaystyle{ \mathbf{G}=h \mathbf{A}+k \mathbf{B}+l \mathbf{C} }[/math] with [math]\displaystyle{ h }[/math], [math]\displaystyle{ k }[/math], [math]\displaystyle{ l }[/math] as integers represents the reciprocal lattice for a crystal lattice [math]\displaystyle{ L }[/math] (defined by [math]\displaystyle{ \mathbf x = p\,\mathbf a+q\,\mathbf b+r\,\mathbf c }[/math]) in real space, we know that [math]\displaystyle{ \mathbf{G}\cdot \mathbf{x} = \mathbf{G} \cdot (p \mathbf{a}+q \mathbf{b}+r \mathbf{c})= 2\pi(hp+kq+lr) = 2\pi n }[/math] with an integer [math]\displaystyle{ n }[/math] due to the known orthogonality between primitive vectors for the reciprocal lattice and those for the crystal lattice. (We use the physical, not crystallographer's, definition for reciprocal lattice vectors which gives the factor of [math]\displaystyle{ 2\pi }[/math].) But notice that this is nothing but the Laue equations. Hence we identify [math]\displaystyle{ \mathbf{\Delta k}= \mathbf{k}_{\mathrm{out}} - \mathbf{k}_{\mathrm{in}} = \mathbf{G} }[/math], means that allowed scattering vectors [math]\displaystyle{ \mathbf{\Delta k}= \mathbf{k}_{\mathrm{out}} - \mathbf{k}_{\mathrm{in}} }[/math] are those equal to reciprocal lattice vectors [math]\displaystyle{ \mathbf{G} }[/math] for a crystal in diffraction, and this is the meaning of the Laue equations. This fact is sometimes called the Laue condition. In this sense, diffraction patterns are a way to experimentally measure the reciprocal lattice for a crystal lattice.

The Laue condition can be rewritten as the following.[2]

[math]\displaystyle{ \begin{align}
  \mathbf{G} &= \mathbf{k}_{\mathrm{out}} - \mathbf{k}_{\mathrm{in}}\\ \rightarrow
 |\mathbf{k}_{\mathrm{in}}|^2 &= |\mathbf{k}_{\mathrm{out}} - \mathbf{G}|^2\\ \rightarrow
 |\mathbf{k}_{\mathrm{in}}|^2 &= |\mathbf{k}_{\mathrm{out}}|^2 - 2\mathbf{k}_{\mathrm{out}}\cdot\mathbf{G} + |\mathbf{G}|^2.
 \end{align} }[/math]

Applying the elastic scattering condition [math]\displaystyle{ |\mathbf{k}_{\mathrm{out}}|^2=|\mathbf{k}_{\mathrm{in}}|^2 }[/math] (In other words, the incoming and diffracted waves are at the same (temporal) frequency. We can also say that the energy per photon does not change.) to the above equation, we obtain

[math]\displaystyle{ 2\mathbf{k}_{\mathrm{out}}\cdot\mathbf{G}=|\mathbf{G}|^2, }[/math]
[math]\displaystyle{ 2{{\mathbf{k}}_{\text{in}}}\cdot (-\mathbf{G})=|\mathbf{G}{{|}^{2}}. }[/math]

The second equation is obtained from the first equation by using [math]\displaystyle{ \mathbf{k}_{\mathrm{out}} - \mathbf{k}_{\mathrm{in}} = \mathbf{G} }[/math].

The result [math]\displaystyle{ 2\mathbf{k}_{\mathrm{out}}\cdot\mathbf{G}=|\mathbf{G}|^2 }[/math] (also [math]\displaystyle{ 2{{\mathbf{k}}_{\text{in}}}\cdot (-\mathbf{G})=|\mathbf{G}{{|}^{2}} }[/math]) is an equation for a plane (geometry) (as the set of all points indicated by [math]\displaystyle{ \mathbf{k}_{\mathrm{out}} }[/math] satisfying this equation) as its equivalent equation [math]\displaystyle{ \mathbf{G}\cdot (2{{\mathbf{k}}_{\text{out}}}-\mathbf{G})=0 }[/math] is a plane equation in geometry. Another equivalent equation, that may be more easier to understand, is [math]\displaystyle{ {{\mathbf{k}}_{\text{out}}}\cdot \widehat{\mathbf{G}}=\frac{1}{2}\left| \mathbf{G} \right| }[/math] (also [math]\displaystyle{ (-{{\mathbf{k}}_{\text{in}}})\cdot \widehat{\mathbf{G}}=\frac{1}{2}\left| \mathbf{G} \right| }[/math]). This indicates the plane that is perpendicular to the straight line between the reciprocal lattice origin [math]\displaystyle{ \mathbf{G}=0 }[/math] and [math]\displaystyle{ \mathbf{G} }[/math] and located at the middle of the line. Such a plane is called Bragg plane.[3] This plane can be understood since [math]\displaystyle{ \mathbf{G} = \mathbf{k}_{\mathrm{out}} - \mathbf{k}_{\mathrm{in}} }[/math] for scattering to occur (It is the Laue condition, equivalent to the Laue equations.) and the elastic scattering [math]\displaystyle{ |\mathbf{k}_{\mathrm{out}}|^2=|\mathbf{k}_{\mathrm{in}}|^2 }[/math] has been assumed so [math]\displaystyle{ \mathbf{G} }[/math], [math]\displaystyle{ \mathbf{k}_{\mathrm{out}} }[/math], and [math]\displaystyle{ -\mathbf{k}_{\mathrm{in}} }[/math] form a rhombus. Each [math]\displaystyle{ \mathbf{G} }[/math] is by definition the wavevector of a plane wave in the Fourier series of a spatial function which periodicity follows the crystal lattice (e.g., the function representing the electronic density of the crystal), wavefronts of each plane wave in the Fourier series is perpendicular to the plane wave's wavevector [math]\displaystyle{ \mathbf{G} }[/math], and these wavefronts are coincident with parallel crystal lattice planes. This means that X-rays are seemingly "reflected" off parallel crystal lattice planes perpendicular [math]\displaystyle{ \mathbf{G} }[/math] at the same angle as their angle of approach to the crystal [math]\displaystyle{ \theta }[/math] with respect to the lattice planes; in the elastic light (typically X-ray)-crystal scattering, parallel crystal lattice planes perpendicular to a reciprocal lattice vector [math]\displaystyle{ \mathbf{G} }[/math] for the crystal lattice play as parallel mirrors for light which, together with [math]\displaystyle{ \mathbf{G} }[/math], incoming (to the crystal) and outgoing (from the crystal by scattering) wavevectors forms a rhombus.

Since the angle between [math]\displaystyle{ \mathbf{k}_{\mathrm{out}} }[/math] and [math]\displaystyle{ \mathbf{G} }[/math] is [math]\displaystyle{ \pi/2 - \theta }[/math], (Due to the mirror-like scattering, the angle between [math]\displaystyle{ \mathbf{k}_{\mathrm{in}} }[/math] and [math]\displaystyle{ \mathbf{G} }[/math] is also [math]\displaystyle{ \pi/2 - \theta }[/math].) [math]\displaystyle{ \mathbf{k}_{\mathrm{out}}\cdot\mathbf{G} = |\mathbf{k}_{\mathrm{out}}||\mathbf{G}|\sin\theta }[/math]. Recall, [math]\displaystyle{ |\mathbf{k}_{\mathrm{out}}| = 2\pi/\lambda }[/math] with [math]\displaystyle{ \lambda }[/math] as the light (typically X-ray) wavelength, and [math]\displaystyle{ |\mathbf{G}| = \frac{2\pi}{d}n }[/math] with [math]\displaystyle{ d }[/math] as the distance between adjacent parallel crystal lattice planes and [math]\displaystyle{ n }[/math] as an integer. With these, we now derive Bragg's law that is equivalent to the Laue equations (also called the Laue condition):

[math]\displaystyle{ \begin{align} 2\mathbf{k}_{\mathrm{out}}\cdot\mathbf{G}=|\mathbf{G}|^2 \\ 2|\mathbf{k}_{\mathrm{out}}||\mathbf{G}|\sin\theta =|\mathbf{G}|^2 \\ 2 (2\pi/\lambda) (2\pi n/d) \sin\theta =(2\pi n/d )^2 \\ 2d\sin\theta=n\lambda. \end{align} }[/math]

References

Notes
  1. More realistically, the oscillators of the lattice should lag behind the incoming wave, and the outcoming wave should lag behind the oscillator. But since the lag is the same at all point of the lattice, the only effect of this correction would be global shift of phase of the outcoming wave, which we are not taking into consideration.
  2. Chaikin, P. M.; Lubensky, T. C.. Principles of condensed matter physics. p. 47. ISBN 0521794501. 
  3. Ashcroft, Neil; Mermin, Nathaniel (1976). Solid State Physics. Saunders College Publishing. pp. 99. ISBN 0030839939.