Reciprocal lattice

From HandWiki
Short description: Fourier transform of real-space lattices, important in solid-state physics
The computer-generated reciprocal lattice of a fictional monoclinic 3D crystal.
A two-dimensional crystal and its reciprocal lattice

In physics, the reciprocal lattice represents the Fourier transform of another lattice (usually a Bravais lattice). In normal usage, the initial lattice (whose transform is represented by the reciprocal lattice) is usually a periodic spatial function in real-space and is also known as the direct lattice. While the direct lattice exists in real-space and is what one would commonly understand as a physical lattice (e.g., a lattice of a crystal, for example), the reciprocal lattice exists in reciprocal space (also known as momentum space or less commonly as K-space, due to the relationship between the Pontryagin duals momentum and position). The reciprocal lattice of a reciprocal lattice is equivalent to the original direct lattice, because the defining equations are symmetrical with respect to the vectors in real and reciprocal space. Mathematically, direct and reciprocal lattice vectors represent covariant and contravariant vectors, respectively.

The reciprocal lattice plays a very fundamental role in most analytic studies of periodic structures, particularly in the theory of diffraction. In neutron and X-ray diffraction, due to the Laue conditions, the momentum difference between incoming and diffracted X-rays of a crystal is a reciprocal lattice vector. The diffraction pattern of a crystal can be used to determine the reciprocal vectors of the lattice. Using this process, one can infer the atomic arrangement of a crystal.

The Brillouin zone is a Wigner-Seitz cell of the reciprocal lattice.

Wave-based description

Reciprocal space

Reciprocal space (also called k-space) provides a way to visualize the results of the Fourier transform of a spatial function. It is similar in role to the frequency domain arising from the Fourier transform of a time dependent function. The domain of the spatial function itself is often referred to as real space. In physical applications, such as crystallography, both real and reciprocal space will often each be two or three dimensional. Whereas spatial dimensions of these two associated spaces will be the same, the spaces will differ in their units of length, so that when the real space has units of length L, its reciprocal space will have units of one divided by the length L−1 (the reciprocal of length).

Reciprocal space comes into play regarding waves, both classical and quantum mechanical. Because a sinusoidal plane wave with unit amplitude can be written as an oscillatory term [math]\displaystyle{ \cos{(kx {-} \omega t {+} \phi_0)} }[/math], with initial phase [math]\displaystyle{ \phi_0 }[/math], angular wavenumber [math]\displaystyle{ k }[/math] and angular frequency [math]\displaystyle{ \omega }[/math], it can be regarded as a function of both [math]\displaystyle{ k }[/math] and [math]\displaystyle{ x }[/math] (and the time-varying part as a function of both [math]\displaystyle{ \omega }[/math] and [math]\displaystyle{ t }[/math]). This complementary role of [math]\displaystyle{ k }[/math] and [math]\displaystyle{ x }[/math] leads to their visualization within complementary spaces. The spatial periodicity of this wave is defined by its wavelength [math]\displaystyle{ \lambda }[/math], where [math]\displaystyle{ k \lambda = 2\pi }[/math]; hence the corresponding wavenumber in reciprocal space will be [math]\displaystyle{ k = 2\pi / \lambda }[/math].

In three dimensions, the corresponding plane wave term becomes [math]\displaystyle{ \cos{(\mathbf{k} {\cdot} \mathbf{r} {-} \omega t {+} \phi_0)} }[/math], which simplifies to [math]\displaystyle{ \cos{(\mathbf{k} {\cdot} \mathbf{r} {+} \phi)} }[/math] at a fixed time [math]\displaystyle{ t }[/math], where [math]\displaystyle{ \mathbf{r} }[/math] is the position vector of a point in real space and now [math]\displaystyle{ \mathbf{k}=2\pi \mathbf{e} / \lambda }[/math] is the wavevector in the three dimensional reciprocal space. The constant [math]\displaystyle{ \phi }[/math] is the phase of the wavefront (a plane of a constant phase) through the origin [math]\displaystyle{ \mathbf{r}=0 }[/math] at time [math]\displaystyle{ t }[/math], and [math]\displaystyle{ \mathbf{e} }[/math] is a unit vector perpendicular to this wavefront. The wavefronts with phase [math]\displaystyle{ \phi + (2\pi)n }[/math], where [math]\displaystyle{ n }[/math] represents any integer number, comprise a set of parallel planes, equally spaced by the wavelength [math]\displaystyle{ \lambda }[/math].

Reciprocal lattice

In general, a geometric lattice is an infinite, regular array of vertices (points) in space, which can be modelled vectorially as a Bravais lattice. Some lattices may be skew, which means that their primary lines may not necessarily be at right angles. A reciprocal lattice is the periodic set of the wave vectors [math]\displaystyle{ \mathbf{k} }[/math] in reciprocal space that make up the Fourier series of any function [math]\displaystyle{ f(\mathbf{r}) }[/math] whose periodicity is compatible with that of an initial direct lattice in real space. Equivalently, a wavevector is a vertex of the reciprocal lattice if it corresponds to a plane wave in real space whose phase at any given time is the same (actually differs by [math]\displaystyle{ (2\pi)n }[/math] with an integer [math]\displaystyle{ n }[/math]) at every direct lattice vertex.

One heuristic approach to constructing the reciprocal lattice in three dimensions is to write the position vector of a vertex of the direct lattice as [math]\displaystyle{ \mathbf{R} = n_{1}\mathbf{a}_{1} {+} n_{2}\mathbf{a}_{2} {+} n_{3}\mathbf{a}_{3} }[/math], where the [math]\displaystyle{ n_i }[/math] are integers defining the vertex and the [math]\displaystyle{ \mathbf{a}_i }[/math] are linearly independent primitive translation vectors (or shortly called primitive vectors) that are characteristic of the lattice. There is then a unique plane wave (up to a factor of negative one), whose wavefront through the origin [math]\displaystyle{ \mathbf{R} = 0 }[/math] contains the direct lattice points at [math]\displaystyle{ \mathbf{a}_2 }[/math] and [math]\displaystyle{ \mathbf{a}_3 }[/math], and with its adjacent wavefront (whose phase differs by [math]\displaystyle{ 2\pi }[/math] or [math]\displaystyle{ -2\pi }[/math] from the former wavefront passing the origin) passing through [math]\displaystyle{ \mathbf{a}_1 }[/math]. Its angular wavevector takes the form [math]\displaystyle{ \mathbf{b}_1 = 2\pi \mathbf{e}_1 / \lambda_{1} }[/math], where [math]\displaystyle{ \mathbf{e}_1 }[/math] is the unit vector perpendicular to these two adjacent wavefronts and the wavelength [math]\displaystyle{ \lambda_1 }[/math] must satisfy [math]\displaystyle{ \lambda_1 = \mathbf{a}_1 \cdot \mathbf{e}_1 }[/math], means that [math]\displaystyle{ \lambda_1 }[/math] is equal to the distance between the two wavefronts. Hence by construction [math]\displaystyle{ \mathbf{a}_1 \cdot \mathbf{b}_1 = 2\pi }[/math] and [math]\displaystyle{ \mathbf{a}_2 \cdot \mathbf{b}_1 = \mathbf{a}_3 \cdot \mathbf{b}_1 = 0 }[/math].

Cycling through the indices in turn, the same method yields three wavevectors [math]\displaystyle{ \mathbf{b}_j }[/math] with [math]\displaystyle{ \mathbf{a}_i \cdot \mathbf{b}_j = 2\pi \, \delta_{ij} }[/math], where the Kronecker delta [math]\displaystyle{ \delta_{ij} }[/math] equals one when [math]\displaystyle{ i=j }[/math] and is zero otherwise. The [math]\displaystyle{ \mathbf{b}_j }[/math] comprise a set of three primitive wavevectors for the reciprocal lattice, each of whose vertices takes the form [math]\displaystyle{ \mathbf{G} = m_{1}\mathbf{b}_{1} {+} m_{2}\mathbf{b}_{2} {+} m_{3}\mathbf{b}_{3} }[/math], where the [math]\displaystyle{ m_j }[/math] are integers. Simple algebra then shows that, for any plane wave with a wavevector [math]\displaystyle{ \mathbf{G} }[/math] on the reciprocal lattice, the total phase shift [math]\displaystyle{ \mathbf{G} \cdot \mathbf{R} }[/math] between the origin and any point [math]\displaystyle{ \mathbf{R} }[/math] on the direct lattice is a (possibly zero) multiple of [math]\displaystyle{ 2\pi }[/math], so the phase will indeed be equal for every direct lattice vertex, in conformity with the reciprocal lattice definition above. (Although any wavevector [math]\displaystyle{ \mathbf{G} }[/math] on the reciprocal lattice does always take this form, this derivation is motivational, rather than rigorous, because it has omitted the proof that no other possibilities exist.)

The Brillouin zone is a primitive cell (more specifically a Wigner-Seitz cell) of the reciprocal lattice, which plays an important role in solid state physics due to Bloch's theorem. In pure mathematics, the dual space of linear forms and the dual lattice provide more abstract generalizations of reciprocal space and the reciprocal lattice.

Mathematical description

Demonstration of relation between real and reciprocal lattice. A real space 2D lattice (red dots) with primitive vectors [math]\displaystyle{ a_1 }[/math] and [math]\displaystyle{ a_2 }[/math] are shown by blue and green arrows respectively. Atop, plane waves of the form [math]\displaystyle{ e^{iG\cdot r} }[/math] are plotted. From this we see that when [math]\displaystyle{ G }[/math] is any integer combination of reciprocal lattice vector basis [math]\displaystyle{ b_1 }[/math] and [math]\displaystyle{ b_2 }[/math] (i.e. any reciprocal lattice vector), the resulting plane waves have the same periodicity of the lattice - that is that any translation from point [math]\displaystyle{ r }[/math] (shown orange) to a point [math]\displaystyle{ R+r }[/math] ([math]\displaystyle{ R }[/math] shown red), the value of the plane wave is the same. These plane waves can be added together and the above relation will still apply.

Assuming a two-dimensional Bravais lattice and labelling each lattice vector by the subscript [math]\displaystyle{ n = (n_1, n_2), }[/math]

[math]\displaystyle{ \mathbf{R}_n = n_1 \mathbf{a}_1 + n_2 \mathbf{a}_2 }[/math] where [math]\displaystyle{ n_1, n_2 \in \mathbb{Z} }[/math].

Taking a function [math]\displaystyle{ f(\mathbf{r}) }[/math] where [math]\displaystyle{ \mathbf{r} }[/math] is a vector from the origin to any position, if [math]\displaystyle{ f(\mathbf{r}) }[/math] follows the periodicity of the lattice, e.g. the electronic density in an atomic crystal, it is useful to write [math]\displaystyle{ f(\mathbf{r}) }[/math] as a multi-dimensional Fourier series

[math]\displaystyle{ \sum_m {f_m e^{i \mathbf{G}_m \cdot \mathbf{r}}} = f\left(\mathbf{r}\right) }[/math]

where now the subscript [math]\displaystyle{ m = (m_1, m_2), }[/math] so this is a double sum.

As [math]\displaystyle{ f(\mathbf{r}) }[/math] follows the periodicity of the lattice, translating [math]\displaystyle{ \mathbf{r} }[/math] by any lattice vector [math]\displaystyle{ \mathbf{R}_n }[/math] we get the same value, hence

[math]\displaystyle{ f(\mathbf{r} + \mathbf{R}_n) = f(\mathbf{r}) }[/math]

Expressing the above instead in terms of their Fourier series we have

[math]\displaystyle{ \sum_m {f_m e^{i \mathbf{G}_m \cdot \mathbf{r}}} = \sum_m {f_m e^{i \mathbf{G}_m \cdot (\mathbf{r} + \mathbf{R}_n)}} = \sum_m {f_m e^{i \mathbf{G}_m \cdot \mathbf{R}_n} \, e^{i \mathbf{G}_m \cdot \mathbf{r}}} }[/math]

Because equality of two Fourier series implies equality of their coefficients, [math]\displaystyle{ e^{i \mathbf{G}_m \cdot \mathbf{R}_n} = 1 }[/math], which only holds when

[math]\displaystyle{ \mathbf{G}_m \cdot \mathbf{R}_n = 2\pi N }[/math] where [math]\displaystyle{ N \in \mathbb{Z}. }[/math]

This criterion restricts the values of [math]\displaystyle{ \mathbf{G}_m }[/math] to vectors that satisfy this relation. Mathematically, the reciprocal lattice is the set of all vectors [math]\displaystyle{ \mathbf{G}_m }[/math] that satisfy the above identity for all lattice point position vectors [math]\displaystyle{ \mathbf{R}_n }[/math]. As such, any function which exhibits the same periodicity of the lattice can be expressed as a Fourier series with angular frequencies taken from the reciprocal lattice.

This reciprocal lattice is itself a Bravais lattice, and the reciprocal of the reciprocal lattice is the original lattice, which reveals the Pontryagin duality of their respective vector spaces.

Two dimensions

For an infinite two-dimensional lattice, defined by its primitive vectors [math]\displaystyle{ \left(\mathbf{a}_1, \mathbf{a}_2\right) }[/math], its reciprocal lattice can be determined by generating its two reciprocal primitive vectors, through the following formulae,

[math]\displaystyle{ \mathbf{G}_m = m_1 \mathbf{b}_1 + m_2 \mathbf{b}_2 }[/math]

Where,

[math]\displaystyle{ \begin{align} \mathbf{b}_1 &= 2\pi \frac{-\mathbf{Q} \, \mathbf{a}_2}{-\mathbf{a}_1 \cdot \mathbf{Q} \, \mathbf{a}_2} = 2\pi \frac{ \mathbf{Q} \, \mathbf{a}_2}{ \mathbf{a}_1 \cdot \mathbf{Q} \, \mathbf{a}_2} \\ \mathbf{b}_2 &= 2\pi \frac{ \mathbf{Q} \, \mathbf{a}_1}{ \mathbf{a}_2 \cdot \mathbf{Q} \, \mathbf{a}_1} \end{align} }[/math]

Here [math]\displaystyle{ \mathbf{Q} }[/math] represents a 90 degree rotation matrix, i.e. a quarter turn. The anti-clockwise rotation and the clockwise rotation can both be used to determine the reciprocal lattice: If [math]\displaystyle{ \mathbf{Q} }[/math] is the anti-clockwise rotation and [math]\displaystyle{ \mathbf{Q'} }[/math] is the clockwise rotation, [math]\displaystyle{ \mathbf{Q}\,\mathbf{v}=-\mathbf{Q'}\,\mathbf{v} }[/math] for all vectors [math]\displaystyle{ \mathbf{v} }[/math]. Thus, using the permutation

[math]\displaystyle{ \sigma = \begin{pmatrix} 1 & 2 \\ 2 & 1 \end{pmatrix} }[/math]

we obtain

[math]\displaystyle{ \begin{align} \mathbf{b}_n = 2\pi \frac{ \mathbf{Q} \, \mathbf{a}_{\sigma(n)}}{ \mathbf{a}_n \cdot \mathbf{Q} \, \mathbf{a}_{\sigma(n)}}=2\pi \frac{ \mathbf{Q}' \, \mathbf{a}_{\sigma(n)}}{ \mathbf{a}_n \cdot \mathbf{Q}' \, \mathbf{a}_{\sigma(n)}}. \end{align} }[/math]

Three dimensions

For an infinite three-dimensional lattice, defined by its primitive vectors [math]\displaystyle{ \left(\mathbf{a_{1}}, \mathbf{a}_2, \mathbf{a}_3\right) }[/math], its reciprocal lattice can be determined by generating its three reciprocal primitive vectors, through the formulae

[math]\displaystyle{ \mathbf{G}_m = m_1 \mathbf{b}_1 + m_2 \mathbf{b}_2 + m_3 \mathbf{b}_3 }[/math]

where the subscript [math]\displaystyle{ m = (m_1, m_2, m_3) }[/math] in three dimensions, and for the scalar triple product [math]\displaystyle{ V = \mathbf{a}_1 \cdot \left(\mathbf{a}_2 \times \mathbf{a}_3\right) }[/math]:

[math]\displaystyle{ \begin{align} \mathbf{b}_1 &= \frac{2\pi}{V} \ \mathbf{a}_2 \times \mathbf{a}_3 \\ \mathbf{b}_2 &= \frac{2\pi}{V} \ \mathbf{a}_3 \times \mathbf{a}_1 \\ \mathbf{b}_3 &= \frac{2\pi}{V} \ \mathbf{a}_1 \times \mathbf{a}_2 \end{align} }[/math]

Using column vector representation of (reciprocal) primitive vectors, the formulae above can be rewritten using matrix inversion:

[math]\displaystyle{ \left[\mathbf{b}_1\mathbf{b}_2\mathbf{b}_3\right]^\mathsf{T} = 2\pi\left[\mathbf{a}_1\mathbf{a}_2\mathbf{a}_3\right]^{-1}. }[/math]

This method appeals to the definition, and allows generalization to arbitrary dimensions. The cross product formula dominates introductory materials on crystallography.

The above definition is called the "physics" definition, as the factor of [math]\displaystyle{ 2 \pi }[/math] comes naturally from the study of periodic structures. An equivalent definition, the "crystallographer's" definition, comes from defining the reciprocal lattice to be [math]\displaystyle{ e^{2 \pi i\mathbf{K} \cdot \mathbf{R}} = 1 }[/math] which changes the definitions of the reciprocal lattice vectors to be

[math]\displaystyle{ \mathbf{b}_1 = \frac{\mathbf{a}_2 \times \mathbf{a}_3}{\mathbf{a}_1 \cdot \left(\mathbf{a}_2 \times \mathbf{a}_3\right)} }[/math]

and so on for the other vectors. The crystallographer's definition has the advantage that the definition of [math]\displaystyle{ \mathbf{b}_1 }[/math] is just the reciprocal magnitude of [math]\displaystyle{ \mathbf{a}_1 }[/math] in the direction of [math]\displaystyle{ \mathbf{a}_2 \times \mathbf{a}_3 }[/math], dropping the factor of [math]\displaystyle{ 2 \pi }[/math]. This can simplify certain mathematical manipulations, and expresses reciprocal lattice dimensions in units of spatial frequency. It is a matter of taste which definition of the lattice is used, as long as the two are not mixed.

Each point [math]\displaystyle{ (hkl) }[/math] in the reciprocal lattice corresponds to a set of lattice planes [math]\displaystyle{ (hkl) }[/math] in the real space lattice. The direction of the reciprocal lattice vector corresponds to the normal to the real space planes. The magnitude of the reciprocal lattice vector is given in reciprocal length and is equal to the reciprocal of the interplanar spacing of the real space planes.

[math]\displaystyle{ n }[/math] dimensions

The formula for [math]\displaystyle{ n }[/math] dimensions can be derived assuming an [math]\displaystyle{ n }[/math]-dimensional real vector space [math]\displaystyle{ V }[/math] with a basis [math]\displaystyle{ (\mathbf{a}_1,\ldots,\mathbf{a}_n) }[/math] and an inner product [math]\displaystyle{ g\colon V\times V\to\mathbf{R} }[/math]. The reciprocal lattice vectors are uniquely determined by the formula [math]\displaystyle{ g(\mathbf{a}_i,\mathbf{b}_j)=2\pi\delta_{ij} }[/math]. Using the permutation

[math]\displaystyle{ \sigma = \begin{pmatrix} 1 & 2 & \cdots &n\\ 2 & 3 & \cdots &1 \end{pmatrix}, }[/math]

they can be determined with the following formula:

[math]\displaystyle{ \mathbf{b}_i = 2\pi\frac{\epsilon_{\sigma^1i\ldots\sigma^ni}}{\omega(\mathbf{a}_1,\ldots,\mathbf{a}_n)}g^{-1}(\mathbf{a}_{\sigma^{n-1}i}\,\lrcorner\ldots\mathbf{a}_{\sigma^1i}\,\lrcorner\,\omega)\in V }[/math]

Here, [math]\displaystyle{ \omega\colon V^n\to\mathbf{R} }[/math] is the volume form, [math]\displaystyle{ g^{-1} }[/math] is the inverse of the vector space isomorphism [math]\displaystyle{ \hat{g}\colon V\to V^* }[/math] defined by [math]\displaystyle{ \hat{g}(v)(w)=g(v,w) }[/math] and [math]\displaystyle{ \lrcorner }[/math] denotes the inner multiplication.

One can verify that this formula is equivalent to the known formulas for the two- and three-dimensional case by using the following facts: In three dimensions, [math]\displaystyle{ \omega(u,v,w)=g(u\times v,w) }[/math] and in two dimensions, [math]\displaystyle{ \omega(v,w)=g(Rv,w) }[/math], where [math]\displaystyle{ R\in\text{SO}(2)\subset L(V,V) }[/math] is the rotation by 90 degrees (just like the volume form, the angle assigned to a rotation depends on the choice of orientation[1]).

Reciprocal lattices of various crystals

Reciprocal lattices for the cubic crystal system are as follows.

Simple cubic lattice

The simple cubic Bravais lattice, with cubic primitive cell of side [math]\displaystyle{ a }[/math], has for its reciprocal a simple cubic lattice with a cubic primitive cell of side [math]\displaystyle{ \textstyle \frac{2\pi}{a} }[/math] ([math]\displaystyle{ \textstyle \frac{1}{a} }[/math] in the crystallographer's definition). The cubic lattice is therefore said to be self-dual, having the same symmetry in reciprocal space as in real space.

Face-centered cubic (FCC) lattice

The reciprocal lattice to an FCC lattice is the body-centered cubic (BCC) lattice, with a cube side of [math]\displaystyle{ \textstyle \frac{4\pi}{a} }[/math].

Consider an FCC compound unit cell. Locate a primitive unit cell of the FCC; i.e., a unit cell with one lattice point. Now take one of the vertices of the primitive unit cell as the origin. Give the basis vectors of the real lattice. Then from the known formulae, you can calculate the basis vectors of the reciprocal lattice. These reciprocal lattice vectors of the FCC represent the basis vectors of a BCC real lattice. Note that the basis vectors of a real BCC lattice and the reciprocal lattice of an FCC resemble each other in direction but not in magnitude.

Body-centered cubic (BCC) lattice

The reciprocal lattice to a BCC lattice is the FCC lattice, with a cube side of [math]\displaystyle{ \textstyle \frac{4\pi}{a} }[/math].

It can be easily proven that only the Bravais lattices which have 90 degrees between [math]\displaystyle{ \left(\mathbf{a}_1, \mathbf{a}_2, \mathbf{a}_3\right) }[/math] (cubic, tetragonal, orthorhombic) have [math]\displaystyle{ \left(\mathbf{b}_1, \mathbf{b}_2, \mathbf{b}_3\right) }[/math] parallel to their real-space vectors.

Simple hexagonal lattice

The reciprocal to a simple hexagonal Bravais lattice with lattice constants c and a is another simple hexagonal lattice with lattice constants [math]\displaystyle{ \textstyle \frac{2\pi}{c} }[/math] and [math]\displaystyle{ \textstyle \frac{4\pi}{a\sqrt{3}} }[/math] rotated through 30° about the c axis with respect to the direct lattice. The simple hexagonal lattice is therefore said to be self-dual, having the same symmetry in reciprocal space as in real space. vectors a 1 = (a(3) 1/2 /2) i + (a/2) j ; a 2 = -(a(3) 1/2 /2) i + (a/2) j ve a 3 = a k

Arbitrary collection of atoms

Shadow of a 118-atom faceted carbon-pentacone's intensity reciprocal-lattice lighting up red in diffraction when intersecting the Ewald sphere.

One path to the reciprocal lattice of an arbitrary collection of atoms comes from the idea of scattered waves in the Fraunhofer (long-distance or lens back-focal-plane) limit as a Huygens-style sum of amplitudes from all points of scattering (in this case from each individual atom).[2] This sum is denoted by the complex amplitude F in the equation below, because it is also the Fourier transform (as a function of spatial frequency or reciprocal distance) of an effective scattering potential in direct space:

[math]\displaystyle{ F[\vec{g}] = \sum_{j=1}^N f_j \left[\vec{g}\right] e^{2 \pi i \vec{g} \cdot \vec{r}_j}. }[/math]

Here g = q/(2π) is the scattering vector q in crystallographer units, N is the number of atoms, fj[g] is the atomic scattering factor for atom j and scattering vector g, while rj is the vector position of atom j. Note that the Fourier phase depends on one's choice of coordinate origin.

For the special case of an infinite periodic crystal, the scattered amplitude F = M Fhkl from M unit cells (as in the cases above) turns out to be non-zero only for integer values of [math]\displaystyle{ (hkl) }[/math], where

[math]\displaystyle{ F_{hkl} = \sum_{j=1}^m f_j\left[g_{hkl}\right] e^{2\pi i \left(h u_j + k v_j + l w_j\right)} }[/math]

when there are j=1,m atoms inside the unit cell whose fractional lattice indices are respectively {uj, vj, wj}. To consider effects due to finite crystal size, of course, a shape convolution for each point or the equation above for a finite lattice must be used instead.

Whether the array of atoms is finite or infinite, one can also imagine an "intensity reciprocal lattice" I[g], which relates to the amplitude lattice F via the usual relation I = F*F where F* is the complex conjugate of F. Since Fourier transformation is reversible, of course, this act of conversion to intensity tosses out "all except 2nd moment" (i.e. the phase) information. For the case of an arbitrary collection of atoms, the intensity reciprocal lattice is therefore:

[math]\displaystyle{ I[\vec{g}] = \sum_{j=1}^N \sum_{k=1}^N f_j \left[\vec{g}\right] f_k \left[\vec{g}\right] e^{2 \pi i \vec{g} \cdot \vec{r}_{jk}}. }[/math]

Here rjk is the vector separation between atom j and atom k. One can also use this to predict the effect of nano-crystallite shape, and subtle changes in beam orientation, on detected diffraction peaks even if in some directions the cluster is only one atom thick. On the down side, scattering calculations using the reciprocal lattice basically consider an incident plane wave. Thus after a first look at reciprocal lattice (kinematic scattering) effects, beam broadening and multiple scattering (i.e. dynamical) effects may be important to consider as well.

Generalization of a dual lattice

Main page: Dual lattice

There are actually two versions in mathematics of the abstract dual lattice concept, for a given lattice L in a real vector space V, of finite dimension.

The first, which generalises directly the reciprocal lattice construction, uses Fourier analysis. It may be stated simply in terms of Pontryagin duality. The dual group V^ to V is again a real vector space, and its closed subgroup L^ dual to L turns out to be a lattice in V^. Therefore, L^ is the natural candidate for dual lattice, in a different vector space (of the same dimension).

The other aspect is seen in the presence of a quadratic form Q on V; if it is non-degenerate it allows an identification of the dual space V* of V with V. The relation of V* to V is not intrinsic; it depends on a choice of Haar measure (volume element) on V. But given an identification of the two, which is in any case well-defined up to a scalar, the presence of Q allows one to speak to the dual lattice to L while staying within V.

In mathematics, the dual lattice of a given lattice L in an abelian locally compact topological group G is the subgroup L of the dual group of G consisting of all continuous characters that are equal to one at each point of L.

In discrete mathematics, a lattice is a locally discrete set of points described by all integral linear combinations of dim = n linearly independent vectors in Rn. The dual lattice is then defined by all points in the linear span of the original lattice (typically all of R^n) with the property that an integer results from the inner product with all elements of the original lattice. It follows that the dual of the dual lattice is the original lattice.

Furthermore, if we allow the matrix B to have columns as the linearly independent vectors that describe the lattice, then the matrix

[math]\displaystyle{ A = B\left(B^\mathsf{T} B\right)^{-1} }[/math] has columns of vectors that describe the dual lattice.

See also

References

  1. Audin, Michèle (2003). Geometry. Springer. p. 69. 
  2. B. E. Warren (1969/1990) X-ray diffraction (Addison-Wesley, Reading MA/Dover, Mineola NY).

External links