• Short-Range Order (SRO): The regular, predictable arrangement of nearest-neighbor atoms extending only over atomic distances of 0.2 to 0.5 nm (e.g., well-defined tetrahedral coordination polyhedra $\text{SiO}_4^{4-}$). SRO is present in both crystalline and amorphous solids as well as liquids.
• Long-Range Order (LRO): The strict, repeating 3D translational periodicity of atomic lattice positions that extends infinitely over thousands or millions of unit cells across macroscopic distances ($> 100 \text{ nm}$). Present only in crystalline solids.

• Radial Distribution Function $g(r)$ (pair correlation):
  • In quartz, $g(r)$ exhibits sharp, discrete peaks at exact crystallographic distances ($r_1, r_2, r_3, \dots, r_{100}$) corresponding to 1st, 2nd, 3rd, and distant coordination shells.
  • In silica glass, $g(r)$ shows sharp peaks only for the 1st and 2nd nearest neighbors ($0.16 \text{ nm}$ for $\text{Si}-\text{O}$ and $0.26 \text{ nm}$ for $\text{O}-\text{O}$), beyond which $g(r)$ rapidly smooths into a continuous, featureless value of $1.0$, indicating complete loss of positional correlation.
• X-ray Diffraction: Quartz produces intense, razor-sharp Bragg diffraction peaks ($2\theta$) due to constructive interference from infinite parallel planes. Amorphous silica produces only a broad, diffuse scattering halo.

In a periodic crystal, all identical interatomic bonds share the exact same equilibrium bond length and dissociation energy. When heated to $T_m$, all bonds break simultaneously, producing a discontinuous first-order thermodynamic phase transition with a sharp melting point and latent heat of fusion $\Delta H_f$. In an amorphous glass, random variations in bond angles and lengths yield a continuous distribution of bond strengths. Weaker bonds break at lower temperatures while stronger bonds remain intact, causing the glass to soften progressively over a broad temperature range without latent heat.

• Bravais Lattice: An infinite periodic 3D mathematical array of points in space where every point has an identical environment. The position vector of any lattice point is given by: $$\mathbf{R} = u \mathbf{a}_1 + v \mathbf{a}_2 + w \mathbf{a}_3 \quad (u, v, w \in \mathbb{Z})$$
• Atomic Basis (Motif): A group of one or more atoms (with specific chemical identities and spatial coordinates) attached identically to every single Bravais lattice point.

A Unit Cell is the smallest representative repeating volume of a crystal that, when translated by integer multiples of the lattice vectors without rotation, tiles 3D space completely. It is fully defined by six independent geometric lattice parameters:

• Three axial edge lengths: $a, b, c$
• Three interaxial angles: $\alpha$ (between $b$ and $c$), $\beta$ (between $a$ and $c$), and $\gamma$ (between $a$ and $b$).

• Primitive Cell: The minimum-volume unit cell that contains exactly one single lattice point ($N = 1$). For FCC, the primitive vectors connect corner points to face centers: $\mathbf{a}_1' = \frac{a}{2}(\mathbf{j} + \mathbf{k})$, $\mathbf{a}_2' = \frac{a}{2}(\mathbf{i} + \mathbf{k})$, $\mathbf{a}_3' = \frac{a}{2}(\mathbf{i} + \mathbf{j})$, forming a rhombohedral cell with volume $V_p = a^3/4$ and interaxial angles of $60^\circ$.
• Conventional Unit Cell: A larger cell chosen to preserve the full orthogonal cubic symmetry of the lattice. For FCC, the conventional cubic cell contains 8 corners ($8 \times \frac{1}{8} = 1$) plus 6 face centers ($6 \times \frac{1}{2} = 3$), containing exactly 4 lattice points ($N = 4$) with volume $V_c = a^3$.

Crystal System	Axial Edge Relationships	Interaxial Angle Relationships	Essential Symmetry Element
1. Cubic	$a = b = c$	$\alpha = \beta = \gamma = 90^\circ$	Four 3-fold rotation axes ($\langle 111 \rangle$)
2. Tetragonal	$a = b \ne c$	$\alpha = \beta = \gamma = 90^\circ$	One 4-fold rotation axis ($[001]$)
3. Orthorhombic	$a \ne b \ne c$	$\alpha = \beta = \gamma = 90^\circ$	Three mutually perpendicular 2-fold axes
4. Hexagonal	$a = b \ne c$	$\alpha = \beta = 90^\circ, \gamma = 120^\circ$	One 6-fold rotation axis ($[0001]$)
5. Rhombohedral	$a = b = c$	$\alpha = \beta = \gamma \ne 90^\circ < 120^\circ$	One 3-fold rotation axis
6. Monoclinic	$a \ne b \ne c$	$\alpha = \gamma = 90^\circ \ne \beta$	One 2-fold rotation axis or mirror plane
7. Triclinic	$a \ne b \ne c$	$\alpha \ne \beta \ne \gamma \ne 90^\circ$	Center of inversion or identity only

• Highest Symmetry: Cubic System (48 total symmetry operations, characterized by four 3-fold axes along the cube body diagonals).
• Lowest Symmetry: Triclinic System (0 or 1 symmetry operation, no right angles, no equal axes, completely anisotropic).

• $\alpha$-Iron (Ferrite, BCC): Cubic ($a = 0.2866 \text{ nm}$).
• $\alpha$-Titanium (HCP): Hexagonal ($a = 0.295 \text{ nm}, c = 0.468 \text{ nm}$).
• $\beta$-Tin (White tin): Tetragonal ($a = b = 0.583 \text{ nm}, c = 0.318 \text{ nm}$).

The four centering modes are:

• Primitive (P): Lattice points at cell corners only ($N = 1$).
• Body-Centered (I): Corner points plus one point at the body center ($N = 2$).
• Face-Centered (F): Corner points plus one point at each of the six face centers ($N = 4$).
• Base/Side-Centered (C): Corner points plus points at the centers of two opposing parallel faces ($N = 2$).

Theoretically, $7 \text{ crystal systems} \times 4 \text{ centering types} = \mathbf{28 \text{ combinations}}$.

Of the 28 theoretical combinations, 14 are mathematically redundant because:

1. Some combinations violate the required point-group rotational symmetry of the crystal system (e.g., base-centering a cubic system destroys the four 3-fold diagonal axes, degrading it to tetragonal).
2. Other combinations can be shown via a simple coordinate rotation and reduction to be identical to a smaller unit cell already included among the 14 (e.g., FCT $\equiv$ BCT, Base-centered tetragonal $\equiv$ Primitive tetragonal).

Consider an FCT cell with base dimensions $a_{\text{FCT}} \times a_{\text{FCT}}$ and vertical height $c_{\text{FCT}}$, containing face-centered points on the top, bottom, and side faces:

• Connect the face-centering points of the four vertical side faces. These four points form a square lying at height $z = c/2$.
• The side of this new square has length: $$a_{\text{BCT}} = \frac{a_{\text{FCT}}}{\sqrt{2}}$$ rotated by $45^\circ$ relative to the FCT axes.
• The corner points of the new cell are the base face-centering points at $z = 0$ and $z = c$, and the central point is the original corner point at $(a/2, a/2, c/2)$.
• This new cell has dimensions $a_{\text{BCT}} \times a_{\text{BCT}} \times c_{\text{FCT}}$ with one point at the center: exactly a Body-Centered Tetragonal (BCT) lattice with volume $V_{\text{BCT}} = V_{\text{FCT}} / 2$.

By international convention, the cell with the smaller base ($a_{\text{BCT}}$) is selected. Hence, FCT is redundant.

• Number of Atoms per Cell ($n$): Atoms sit only at the 8 cube corners, each shared by 8 adjacent cubes: $$n = 8 \times \frac{1}{8} = \mathbf{1 \text{ atom/unit cell}}$$
• Coordination Number ($\text{CN}$): Each corner atom touches 6 nearest neighbors (along $\pm x, \pm y, \pm z$): $$\text{CN} = \mathbf{6}$$

Atoms touch along the cube edges ($\langle 100 \rangle$ directions):

$$a = 2 R$$

Unit cell volume:

$$V_c = a^3 = (2 R)^3 = 8 R^3$$

The Atomic Packing Factor is the fraction of unit cell volume occupied by hard-sphere atoms:

$$\text{APF} = \frac{V_{\text{atoms}}}{V_{\text{cell}}} = \frac{n \left( \frac{4}{3} \pi R^3 \right)}{V_c}$$

Substitute $n = 1$ and $V_c = 8 R^3$:

$$\text{APF}_{\text{SC}} = \frac{1 \times \left( \frac{4}{3} \pi R^3 \right)}{8 R^3} = \frac{4 \pi}{24} = \frac{\pi}{6} \approx \mathbf{0.5236 \approx 52.4\%}$$

Why SC is extremely rare in metals: Metallic bonding is non-directional and seeks to maximize atomic packing density and nearest-neighbor coordination. In SC, nearly 48% of the volume is empty void space, and the coordination number is only 6. Metallic elements achieve much lower thermodynamic free energy by packing into BCC ($\text{CN}=8, \text{APF}=0.68$) or close-packed FCC/HCP ($\text{CN}=12, \text{APF}=0.74$). Only Polonium adopts SC due to directional covalent $p$-orbital bonding.

• Number of Atoms per Cell ($n$): 8 corner atoms (each shared by 8 cells) plus 1 whole atom at the center: $$n = 8 \times \frac{1}{8} + 1 = 1 + 1 = \mathbf{2 \text{ atoms/unit cell}}$$
• Coordination Number ($\text{CN}$): The central atom touches all 8 corner atoms along the four body diagonals: $$\text{CN} = \mathbf{8}$$

In BCC, atoms touch along the body diagonals ($\langle 111 \rangle$ directions):

• Face diagonal length of cube: $d_{\text{face}} = \sqrt{a^2 + a^2} = a\sqrt{2}$.
• Body diagonal length: $d_{\text{body}} = \sqrt{d_{\text{face}}^2 + a^2} = \sqrt{2a^2 + a^2} = a\sqrt{3}$.
• Along the body diagonal, the central atom (diameter $2R$) is sandwiched between two corner atoms (each contributing radius $R$): $$d_{\text{body}} = R + 2R + R = 4 R$$

Equating the expressions:

$$a \sqrt{3} = 4 R \implies a = \frac{4 R}{\sqrt{3}}$$

Unit cell volume:

$$V_c = a^3 = \left( \frac{4 R}{\sqrt{3}} \right)^3 = \frac{64 R^3}{3 \sqrt{3}}$$

$$\text{APF} = \frac{n \cdot V_{\text{atom}}}{V_c} = \frac{2 \times \left( \frac{4}{3} \pi R^3 \right)}{\frac{64 R^3}{3 \sqrt{3}}}$$ $$\text{APF}_{\text{BCC}} = \frac{\frac{8}{3} \pi R^3}{\frac{64}{3 \sqrt{3}} R^3} = \frac{8 \pi}{3} \times \frac{3 \sqrt{3}}{64} = \frac{8 \sqrt{3} \pi}{64} = \frac{\sqrt{3} \pi}{8}$$ $$\text{APF}_{\text{BCC}} = \frac{1.73205 \times 3.14159}{8} = \frac{5.4414}{8} \approx \mathbf{0.6802 \approx 68.0\%}$$

• Number of Atoms per Cell ($n$): 8 corner atoms (each shared by 8 cells) plus 6 face-centered atoms (each shared by 2 cells): $$n = 8 \times \frac{1}{8} + 6 \times \frac{1}{2} = 1 + 3 = \mathbf{4 \text{ atoms/unit cell}}$$
• Coordination Number ($\text{CN}$): Each face-centered atom touches 4 corner atoms on its face, 4 face-centered atoms in the plane ahead, and 4 in the plane behind: $$\text{CN} = \mathbf{12}$$

Atoms touch along the face diagonals ($\langle 110 \rangle$ directions):

• Face diagonal length: $d_{\text{face}} = \sqrt{a^2 + a^2} = a\sqrt{2}$.
• Along the face diagonal, the face-center atom is in contact with two corner atoms: $$d_{\text{face}} = R + 2R + R = 4 R$$

Equating the expressions:

$$a \sqrt{2} = 4 R \implies a = \frac{4 R}{\sqrt{2}} = 2\sqrt{2} R \approx 2.8284 R$$

Unit cell volume:

$$V_c = a^3 = (2\sqrt{2} R)^3 = 16\sqrt{2} R^3$$

$$\text{APF} = \frac{n \cdot V_{\text{atom}}}{V_c} = \frac{4 \times \left( \frac{4}{3} \pi R^3 \right)}{16\sqrt{2} R^3}$$ $$\text{APF}_{\text{FCC}} = \frac{\frac{16}{3} \pi R^3}{16\sqrt{2} R^3} = \frac{\pi}{3\sqrt{2}} = \frac{\sqrt{2}\pi}{6}$$ $$\text{APF}_{\text{FCC}} = \frac{1.41421 \times 3.14159}{6} = \frac{4.4429}{6} \approx \mathbf{0.7405 \approx 74.0\%}$$

Kepler Conjecture & Maximum Limit: In 1611, Johannes Kepler conjectured that no packing of equal spheres in 3D space can have an average density exceeding $\pi / (3\sqrt{2}) \approx 0.74048$. Proved rigorously by Thomas Hales in 1998, FCC and HCP represent the densest possible packings of identical spheres.

• Number of Atoms per Cell ($n$):
  • 12 top/bottom base corners, each shared by 6 cells: $12 \times \frac{1}{6} = 2$ atoms.
  • 2 top/bottom face centers, each shared by 2 cells: $2 \times \frac{1}{2} = 1$ atom.
  • 3 interior mid-plane atoms completely inside the cell: $3 \times 1 = 3$ atoms.
  • Total: $n = 2 + 1 + 3 = \mathbf{6 \text{ atoms/hexagonal unit cell}}$.
• Coordination Number ($\text{CN}$): Each atom touches 6 in its own basal plane, 3 in the plane above, and 3 in the plane below: $$\text{CN} = \mathbf{12}$$

The base is a regular hexagon of side length $a$, composed of six equilateral triangles of side $a$:

$$\text{Area of one triangle} = \frac{1}{2} a \left( a \frac{\sqrt{3}}{2} \right) = \frac{\sqrt{3}}{4} a^2$$ $$\text{Base Area} = 6 \times \left( \frac{\sqrt{3}}{4} a^2 \right) = \frac{3\sqrt{3}}{2} a^2$$

Multiplying by cell height $c$:

$$V_c = \text{Base Area} \times c = \mathbf{\frac{3\sqrt{3}}{2} a^2 c}$$

Consider the regular tetrahedron formed by three mid-plane atoms (triangle of side $a$) and one basal atom directly above them:

• All four atoms touch, so all edge lengths equal $a = 2R$.
• The distance from a corner to the centroid of the base equilateral triangle is: $r_{\text{cent}} = \frac{a}{\sqrt{3}}$.
• The vertical height of the tetrahedron is half the unit cell height ($h = c/2$). By the Pythagorean theorem: $$h^2 + r_{\text{cent}}^2 = a^2 \implies \left( \frac{c}{2} \right)^2 + \left( \frac{a}{\sqrt{3}} \right)^2 = a^2$$ $$\frac{c^2}{4} + \frac{a^2}{3} = a^2 \implies \frac{c^2}{4} = \frac{2}{3} a^2 \implies \frac{c^2}{a^2} = \frac{8}{3}$$ $$\left( \frac{c}{a} \right)_{\text{ideal}} = \sqrt{\frac{8}{3}} = \frac{2\sqrt{2}}{\sqrt{3}} \approx \mathbf{1.63299 \approx 1.633}$$

Substitute $c = a\sqrt{8/3}$ and $a = 2R$ into $V_c$:

$$V_c = \frac{3\sqrt{3}}{2} (2R)^2 \left( 2R \sqrt{\frac{8}{3}} \right) = \frac{3\sqrt{3}}{2} (4R^2) \left( 2R \frac{2\sqrt{2}}{\sqrt{3}} \right) = 24\sqrt{2} R^3$$

Now calculate APF with $n = 6$ atoms:

$$\text{APF}_{\text{HCP}} = \frac{6 \left( \frac{4}{3} \pi R^3 \right)}{24\sqrt{2} R^3} = \frac{8 \pi R^3}{24\sqrt{2} R^3} = \frac{\pi}{3\sqrt{2}} \approx \mathbf{0.7405 \approx 74.0\%}$$

Density is defined as mass divided by volume: $\rho = \frac{m_{\text{cell}}}{V_c}$.

• $n$: Number of atoms associated with each unit cell (dimensionless; 1 for SC, 2 for BCC, 4 for FCC, 6 for HCP).
• $A$: Atomic weight (molar mass) of the element ($\text{g/mol}$).
• The mass of a single atom is $m_{\text{atom}} = \frac{A}{N_A}$ (where $N_A = 6.022 \times 10^{23} \text{ atoms/mol}$ is Avogadro's number).
• The total mass contained in one unit cell is: $m_{\text{cell}} = n \cdot m_{\text{atom}} = \frac{n A}{N_A}$.
• Dividing by unit cell volume $V_c$: $$\rho = \frac{n \cdot A}{V_c \cdot N_A}$$

Units: If $A$ is in $\text{g/mol}$ and $V_c$ is in $\text{cm}^3$, $\rho$ evaluates directly in $\mathbf{\text{g/cm}^3}$. To convert to SI ($\text{kg/m}^3$), multiply by $10^3$.

The crystallographic density assumes a mathematically perfect, defect-free single crystal. In real engineering materials, bulk density is slightly lower (typically 0.1% to 2% lower) due to:

1. Thermal Equilibrium Vacancies: Missing atoms at lattice sites ($N_v/N \sim 10^{-4}$ near melting).
2. Dislocations: Dilatational strain fields around dislocation cores lower average local density.
3. Grain Boundaries: In polycrystalline metals, grain boundaries are regions of atomic mismatch with excess free volume.
4. Processing Porosity: Microscopic gas pores, shrinkage voids, and microcracks from casting or additive manufacturing.

• Cubic: $V_c = a^3$
• Tetragonal: $V_c = a^2 c$
• Hexagonal: $V_c = \frac{3\sqrt{3}}{2} a^2 c \approx 2.5981 a^2 c$

$$a = 2\sqrt{2} R = 2\sqrt{2} (1.278 \times 10^{-8} \text{ cm}) = 2.8284 \times (1.278 \times 10^{-8} \text{ cm}) = \mathbf{3.6147 \times 10^{-8} \text{ cm} = 0.3615 \text{ nm}}$$

$$V_c = a^3 = (3.6147 \times 10^{-8} \text{ cm})^3 = \mathbf{4.7229 \times 10^{-23} \text{ cm}^3}$$

For FCC, $n = 4$ atoms/cell:

$$\rho = \frac{n \cdot A_{\text{Cu}}}{V_c \cdot N_A} = \frac{4 \times 63.55 \text{ g/mol}}{(4.7229 \times 10^{-23} \text{ cm}^3) \times (6.022 \times 10^{23} \text{ atoms/mol})}$$ $$\rho = \frac{254.20 \text{ g}}{28.4413 \text{ cm}^3} = \mathbf{8.9377 \text{ g/cm}^3 \approx 8.94 \text{ g/cm}^3}$$

In SI units:

$$\rho = 8.938 \times 10^3 \text{ kg/m}^3 = \mathbf{8,938 \text{ kg/m}^3}$$

Accuracy Comparison: Measured handbook density is $8.94 \text{ g/cm}^3$. Accuracy is within $\mathbf{0.02\%}$.

$$a = \frac{4 R}{\sqrt{3}} = \frac{4 (1.241 \times 10^{-8} \text{ cm})}{\sqrt{3}} = \frac{4.964 \times 10^{-8} \text{ cm}}{1.73205} = \mathbf{2.8660 \times 10^{-8} \text{ cm} = 0.2866 \text{ nm}}$$

$$V_c = a^3 = (2.8660 \times 10^{-8} \text{ cm})^3 = \mathbf{2.3541 \times 10^{-23} \text{ cm}^3}$$

For BCC, $n = 2$ atoms/cell:

$$\rho = \frac{n \cdot A_{\text{Fe}}}{V_c \cdot N_A} = \frac{2 \times 55.85 \text{ g/mol}}{(2.3541 \times 10^{-23} \text{ cm}^3) \times (6.022 \times 10^{23} \text{ atoms/mol})}$$ $$\rho = \frac{111.70 \text{ g}}{14.1764 \text{ cm}^3} = \mathbf{7.8793 \text{ g/cm}^3 \approx 7.88 \text{ g/cm}^3}$$

In SI units: $\rho = \mathbf{7,879 \text{ kg/m}^3}$. Matches handbook value ($7.87 \text{ g/cm}^3$) within $0.1\%$.

$$\frac{c}{a} = \frac{0.4686 \text{ nm}}{0.2951 \text{ nm}} = \mathbf{1.5879 \approx 1.588}$$

Comparing with the ideal hard-sphere ratio $(c/a)_{\text{ideal}} = \sqrt{8/3} \approx 1.633$:

$$\frac{1.588}{1.633} = 0.972 \implies \mathbf{2.8\% \text{ compression}}$$

The unit cell of titanium is slightly compressed along the $c$-axis, indicating anisotropic electronic bonding along the basal planes.

Convert dimensions to centimeters: $a = 2.951 \times 10^{-8} \text{ cm}$, $c = 4.686 \times 10^{-8} \text{ cm}$:

$$V_c = \frac{3\sqrt{3}}{2} a^2 c = 2.598076 \times (2.951 \times 10^{-8} \text{ cm})^2 \times (4.686 \times 10^{-8} \text{ cm})$$ $$a^2 = 8.7084 \times 10^{-16} \text{ cm}^2$$ $$V_c = 2.598076 \times (8.7084 \times 10^{-16}) \times (4.686 \times 10^{-8}) = \mathbf{1.0602 \times 10^{-22} \text{ cm}^3}$$

For a complete hexagonal unit cell, $n = 6$ atoms/cell:

$$\rho = \frac{n \cdot A_{\text{Ti}}}{V_c \cdot N_A} = \frac{6 \times 47.87 \text{ g/mol}}{(1.0602 \times 10^{-22} \text{ cm}^3) \times (6.022 \times 10^{23} \text{ atoms/mol})}$$ $$\rho = \frac{287.22 \text{ g}}{63.8452 \text{ cm}^3} = \mathbf{4.4987 \text{ g/cm}^3 \approx 4.50 \text{ g/cm}^3}$$

In SI units: $\rho = \mathbf{4,500 \text{ kg/m}^3}$. Matches handbook value ($4.51 \text{ g/cm}^3$) within $0.2\%$.

$$\rho = \frac{n \cdot A}{V_c \cdot N_A} \implies V_c = \frac{n \cdot A}{\rho \cdot N_A}$$ $$V_c = \frac{2 \times 95.94 \text{ g/mol}}{(10.22 \text{ g/cm}^3) \times (6.022 \times 10^{23} \text{ atoms/mol})} = \frac{191.88}{6.1545 \times 10^{24}} = \mathbf{3.1177 \times 10^{-23} \text{ cm}^3}$$

$$a = V_c^{1/3} = (3.1177 \times 10^{-23} \text{ cm}^3)^{1/3} = 3.1473 \times 10^{-8} \text{ cm} = \mathbf{0.3147 \text{ nm}}$$

For BCC, $a = \frac{4 R}{\sqrt{3}} \implies R = \frac{a \sqrt{3}}{4}$:

$$R = \frac{(0.31473 \text{ nm}) \sqrt{3}}{4} = \frac{(0.31473 \text{ nm}) \times 1.73205}{4} = \frac{0.54513}{4} = \mathbf{0.1363 \text{ nm}}$$

Unit cell volume:

$$V_c = a^3 = (3.615 \times 10^{-8} \text{ cm})^3 = 4.724 \times 10^{-23} \text{ cm}^3$$

From the density equation:

$$n = \frac{\rho \cdot V_c \cdot N_A}{A}$$ $$n = \frac{(8.94 \text{ g/cm}^3) \times (4.724 \times 10^{-23} \text{ cm}^3) \times (6.022 \times 10^{23} \text{ atoms/mol})}{63.55 \text{ g/mol}}$$ $$n = \frac{254.33}{63.55} = \mathbf{4.002 \approx 4}$$

Because $n = 4$, the metal crystallizes in the Face-Centered Cubic (FCC) structure.

Atomic weight $A = 63.55 \text{ g/mol}$ identifies the metal as Copper (Cu).

For FCC, $a = 2\sqrt{2} R \implies R = \frac{a}{2\sqrt{2}}$:

$$R = \frac{0.3615 \text{ nm}}{2\sqrt{2}} = \frac{0.3615}{2.8284} = \mathbf{0.1278 \text{ nm} = 1.278 \text{ \AA}}$$

$$LD = \frac{\text{Equivalent number of atoms intersected}}{\text{Line vector length } L}$$

• Direction $[110]$ (Face Diagonal):
  • Vector length: $L = a\sqrt{2} = 4R$.
  • Line intersects: $1/2$ corner atom + 1 face-center atom + $1/2$ corner atom = 2 atoms.
  • $$LD_{[110]} = \frac{2}{a\sqrt{2}} = \frac{\sqrt{2}}{a} = \frac{2}{4R} = \mathbf{\frac{1}{2R}}$$ (Close-packed, atoms in continuous contact: $LD = 100\%$ linear packing).
• Direction $[100]$ (Cube Edge):
  • Vector length: $L = a$.
  • Line intersects: $1/2$ corner atom + $1/2$ corner atom = 1 atom.
  • $$LD_{[100]} = \frac{1}{a} = \frac{1}{2\sqrt{2}R} \approx \frac{0.354}{R}$$
• Direction $[111]$ (Body Diagonal):
  • Vector length: $L = a\sqrt{3}$.
  • Line intersects: $1/2$ corner atom at $(0,0,0)$ and $1/2$ corner atom at $(1,1,1)$ = 1 atom (no atom at the center in FCC!).
  • $$LD_{[111]} = \frac{1}{a\sqrt{3}} \approx \frac{0.204}{R}$$

• $$LD_{[110]} = \frac{\sqrt{2}}{a} = \frac{1.4142}{0.4049 \text{ nm}} = \mathbf{3.493 \text{ atoms/nm}}$$
• $$LD_{[100]} = \frac{1}{a} = \frac{1}{0.4049 \text{ nm}} = \mathbf{2.470 \text{ atoms/nm}}$$
• $$LD_{[111]} = \frac{1}{a\sqrt{3}} = \frac{1}{(0.4049) \times 1.732} = \frac{1}{0.7013 \text{ nm}} = \mathbf{1.426 \text{ atoms/nm}}$$

Proof: $LD_{[110]} (3.49) > LD_{[100]} (2.47) > LD_{[111]} (1.43)$. The $[110]$ direction has the maximum linear density, explaining why $\langle 110 \rangle$ is the universal slip direction in FCC metals.

• Direction $[111]$ (Body Diagonal):
  • Vector length: $L = a\sqrt{3} = 4R$.
  • Line intersects: $1/2$ corner + 1 center atom + $1/2$ corner = 2 atoms.
  • $$LD_{[111]} = \frac{2}{a\sqrt{3}} = \frac{2}{4R} = \mathbf{\frac{1}{2R}}$$ (Close-packed, 100% linear density).
• Direction $[100]$ (Cube Edge):
  • Vector length: $L = a$.
  • Line intersects: $1/2$ corner + $1/2$ corner = 1 atom.
  • $$LD_{[100]} = \frac{1}{a} = \frac{\sqrt{3}}{4R} \approx \frac{0.433}{R}$$
• Direction $[110]$ (Face Diagonal):
  • Vector length: $L = a\sqrt{2}$.
  • Line intersects: $1/2$ corner + $1/2$ corner = 1 atom (no face-center atom in BCC!).
  • $$LD_{[110]} = \frac{1}{a\sqrt{2}} = \frac{\sqrt{3}}{4\sqrt{2}R} \approx \frac{0.306}{R}$$

• $$LD_{[111]} = \frac{2}{a\sqrt{3}} = \frac{2}{(0.3165 \text{ nm}) \times 1.73205} = \frac{2}{0.5482} = \mathbf{3.648 \text{ atoms/nm}}$$
• $$LD_{[100]} = \frac{1}{a} = \frac{1}{0.3165 \text{ nm}} = \mathbf{3.160 \text{ atoms/nm}}$$
• $$LD_{[110]} = \frac{1}{a\sqrt{2}} = \frac{1}{(0.3165 \text{ nm}) \times 1.4142} = \frac{1}{0.4476} = \mathbf{2.234 \text{ atoms/nm}}$$

Dislocation glide occurs preferentially along directions of maximum atomic packing density (where Burgers vector length $|\mathbf{b}|$ is minimized). Because $[111]$ has the maximum linear density ($3.65 \text{ nm}^{-1}$), $\langle 111 \rangle$ is the universal close-packed slip direction in all BCC metals.

$$PD = \frac{\text{Equivalent atoms intersecting the 2D plane}}{\text{Plane Area } A_p}$$

• Plane $(111)$ (Equilateral Triangle):
  • Corners are at $(1,0,0), (0,1,0), (0,0,1)$. Each side has length $L = a\sqrt{2}$.
  • Plane area: $A_p = \frac{\sqrt{3}}{4} L^2 = \frac{\sqrt{3}}{4} (2 a^2) = \frac{\sqrt{3}}{2} a^2$.
  • Atoms on plane: 3 corner atoms (each inside triangle by $60^\circ / 360^\circ = 1/6$) + 3 face-center atoms (each bisected by the plane edge, $1/2$ inside): $$N = 3 \times \frac{1}{6} + 3 \times \frac{1}{2} = \frac{1}{2} + \frac{3}{2} = \mathbf{2 \text{ atoms}}$$
  • $$PD_{(111)} = \frac{2}{\frac{\sqrt{3}}{2} a^2} = \mathbf{\frac{4}{\sqrt{3} a^2}} = \frac{1}{4\sqrt{3}R^2}$$
  • Planar packing fraction: $\frac{2 \pi R^2}{\frac{\sqrt{3}}{2} a^2} = \frac{\pi}{2\sqrt{3}} \approx \mathbf{0.9069 \approx 90.7\%}$ (Maximum 2D packing of circles).
• Plane $(100)$ (Square):
  • Area: $A_p = a^2$.
  • Atoms: 4 corners ($4 \times 1/4 = 1$) + 1 face-center ($1$) = 2 atoms.
  • $$PD_{(100)} = \frac{2}{a^2} = \frac{1}{4 R^2}$$
  • Planar packing fraction: $\frac{2 \pi R^2}{8 R^2} = \frac{\pi}{4} \approx \mathbf{0.7854 \approx 78.5\%}$.
• Plane $(110)$ (Rectangle):
  • Area: $A_p = a \times a\sqrt{2} = a^2 \sqrt{2}$.
  • Atoms: 4 corners ($4 \times 1/4 = 1$) + 2 face-centers ($2 \times 1/2 = 1$) = 2 atoms.
  • $$PD_{(110)} = \frac{2}{a^2 \sqrt{2}} = \frac{\sqrt{2}}{a^2} \approx \frac{1.414}{a^2}$$
  • Planar packing fraction: $\frac{2 \pi R^2}{a^2\sqrt{2}} = \frac{\pi}{4\sqrt{2}} \approx \mathbf{0.5554 \approx 55.5\%}$.

$a^2 = (0.3524 \text{ nm})^2 = 0.124186 \text{ nm}^2$:

• $$PD_{(111)} = \frac{4}{\sqrt{3} (0.124186 \text{ nm}^2)} = \frac{4}{0.215097} = \mathbf{18.596 \text{ atoms/nm}^2}$$
• $$PD_{(100)} = \frac{2}{0.124186 \text{ nm}^2} = \mathbf{16.105 \text{ atoms/nm}^2}$$
• $$PD_{(110)} = \frac{\sqrt{2}}{0.124186 \text{ nm}^2} = \frac{1.4142}{0.124186} = \mathbf{11.388 \text{ atoms/nm}^2}$$

Conclusion: $PD_{(111)} (18.60) > PD_{(100)} (16.11) > PD_{(110)} (11.39)$. The $\{111\}$ planes have the highest planar atomic density, proving they are the close-packed slip planes in FCC metals.

• Plane $(110)$ (Diagonal Rectangle):
  • Dimensions: $a \times a\sqrt{2}$. Area: $A_p = a^2 \sqrt{2}$.
  • Atoms on plane: 4 corners ($4 \times 1/4 = 1$) + 1 whole center atom = 2 atoms.
  • $$PD_{(110)} = \frac{2}{a^2 \sqrt{2}} = \mathbf{\frac{\sqrt{2}}{a^2}} = \frac{3\sqrt{2}}{16 R^2} \approx \frac{0.265}{R^2}$$
  • Planar packing fraction: $\frac{2 \pi R^2}{a^2 \sqrt{2}} = \frac{2 \pi R^2}{(16 R^2 / 3)\sqrt{2}} = \frac{3 \pi}{8 \sqrt{2}} \approx \mathbf{0.833 \approx 83.3\%}$.
• Plane $(100)$ (Cube Face):
  • Area: $A_p = a^2$.
  • Atoms: 4 corners ($4 \times 1/4 = 1$) = 1 atom (no face center in BCC!).
  • $$PD_{(100)} = \frac{1}{a^2} = \frac{3}{16 R^2} \approx \frac{0.1875}{R^2}$$
  • Planar packing fraction: $\frac{\pi R^2}{a^2} = \frac{3 \pi}{16} \approx \mathbf{0.589 \approx 58.9\%}$.
• Plane $(111)$ (Triangular Incline):
  • Side: $L = a\sqrt{2}$. Area: $A_p = \frac{\sqrt{3}}{2} a^2$.
  • Atoms: 3 corners ($3 \times 1/6 = 1/2$) = 0.5 atoms (body center does NOT lie on the plane!).
  • $$PD_{(111)} = \frac{0.5}{\frac{\sqrt{3}}{2} a^2} = \frac{1}{\sqrt{3} a^2} \approx \frac{0.577}{a^2}$$
  • Planar packing fraction: $\approx \mathbf{0.340 \approx 34.0\%}$.

$a^2 = (0.2866 \text{ nm})^2 = 0.08214 \text{ nm}^2$:

• $$PD_{(110)} = \frac{\sqrt{2}}{0.08214 \text{ nm}^2} = \frac{1.4142}{0.08214} = \mathbf{17.217 \text{ atoms/nm}^2}$$
• $$PD_{(100)} = \frac{1}{0.08214 \text{ nm}^2} = \mathbf{12.174 \text{ atoms/nm}^2}$$
• $$PD_{(111)} = \frac{1}{\sqrt{3} (0.08214 \text{ nm}^2)} = \frac{0.57735}{0.08214} = \mathbf{7.029 \text{ atoms/nm}^2}$$

The $\{110\}$ planes have the highest planar atomic density ($17.22 \text{ nm}^{-2}$, 83.3% packing fraction) and the largest interplanar spacing ($d_{110} = a/\sqrt{2}$). According to the Peierls-Nabarro model, shear resistance is minimized on planes with maximum interplanar spacing and maximum planar density. Therefore, $\{110\}$ planes are the primary slip planes in all BCC metals.

1. Position Vector: Establish a right-handed Cartesian coordinate system ($x, y, z$) with origin at a unit cell corner. Determine coordinates of the vector head $(x_2, y_2, z_2)$ and tail $(x_1, y_1, z_1)$.
2. Vector Subtraction: Subtract tail coordinates from head coordinates: $\Delta x = x_2 - x_1$, $\Delta y = y_2 - y_1$, $\Delta z = z_2 - z_1$.
3. Clear Fractions / Scale: Divide by lattice parameters ($a, b, c$) to obtain dimensionless projections. Multiply or divide by a common factor to reduce these numbers to the smallest set of integers.
4. Enclose: Enclose the three integers in square brackets without commas: $[uvw]$.

Negative integer indices are represented by placing an overbar (macron) directly above the integer: e.g., $-1$ is written as $\bar{1}$ (pronounced "bar one" or "minus one"), and $-2$ is written as $\bar{2}$.

Tail: $(0, 0, 0)$; Head: $(a, a/2, a/4)$.

• Projections in units of lattice parameter $a$: $$\Delta x = 1 - 0 = 1, \quad \Delta y = \frac{1}{2} - 0 = \frac{1}{2}, \quad \Delta z = \frac{1}{4} - 0 = \frac{1}{4}$$
• Multiply by the lowest common denominator (4) to clear fractions: $$1 \times 4 = 4, \quad \frac{1}{2} \times 4 = 2, \quad \frac{1}{4} \times 4 = 1$$
• Integers: $4, 2, 1$.

The crystallographic direction is $[421]$.

Projections: $1 - 0 = 1, \ 1 - 0 = 1, \ 0 - 0 = 0$.
Integers: $[110]$. Represents a face diagonal on the $xy$-plane.

Projections: $0 - 1 = -1, \ 1 - 0 = 1, \ 0 - 0 = 0$.
Integers with overbar: $[\bar{1}10]$. Face diagonal pointing from $(1,0,0)$ to $(0,1,0)$.

Projections: $1 - 0 = 1, \ 0 - 1 = -1, \ 1 - 0 = 1$.
Integers: $[1\bar{1}1]$. Body diagonal direction.

Projections: $$\Delta x = 0 - \frac{1}{2} = -\frac{1}{2}, \quad \Delta y = 1 - 0 = 1, \quad \Delta z = \frac{1}{2} - 1 = -\frac{1}{2}$$ Multiply by 2 to clear fractions: $$u = -1, \quad v = 2, \quad w = -1$$ Enclosing: $[\bar{1}2\bar{1}]$.

A family of crystallographically equivalent directions is designated by enclosing the indices in angle brackets: $\langle uvw \rangle$. In cubic symmetry, equivalence is generated by permuting indices and altering signs (reversing sense).

The $\langle 100 \rangle$ family contains 6 equivalent directions ($N_{\text{mult}} = 6$):

$$\langle 100 \rangle = \{ [100], [\bar{1}00], [010], [0\bar{1}0], [001], [00\bar{1}] \}$$

• $\langle 110 \rangle$ Family (Face Diagonals, Multiplicity = 12): $$\langle 110 \rangle = \begin{cases} [110], [\bar{1}\bar{1}0], [1\bar{1}0], [\bar{1}10], \\ [101], [\bar{1}0\bar{1}], [10\bar{1}], [\bar{1}01], \\ [011], [0\bar{1}\bar{1}], [01\bar{1}], [0\bar{1}1] \end{cases}$$
• $\langle 111 \rangle$ Family (Body Diagonals, Multiplicity = 8): $$\langle 111 \rangle = \{ [111], [\bar{1}\bar{1}\bar{1}], [\bar{1}11], [1\bar{1}\bar{1}], [1\bar{1}1], [\bar{1}1\bar{1}], [11\bar{1}], [\bar{1}\bar{1}1] \}$$

$$\cos\theta = \frac{\mathbf{v}_1 \cdot \mathbf{v}_2}{|\mathbf{v}_1| |\mathbf{v}_2|} = \frac{u_1 u_2 + v_1 v_2 + w_1 w_2}{\sqrt{u_1^2 + v_1^2 + w_1^2} \sqrt{u_2^2 + v_2^2 + w_2^2}}$$

Note: This simple formula is valid strictly for cubic crystals where coordinate axes are orthogonal and equal in length ($a=b=c, \alpha=\beta=\gamma=90^\circ$).

$$\cos\theta = \frac{(1)(1) + (1)(1) + (0)(1)}{\sqrt{1^2 + 1^2 + 0^2} \sqrt{1^2 + 1^2 + 1^2}} = \frac{1 + 1 + 0}{\sqrt{2} \sqrt{3}} = \frac{2}{\sqrt{6}} \approx \frac{2}{2.4495} = 0.8165$$ $$\theta = \arccos(0.8165) \approx \mathbf{35.26^\circ}$$

• Angle Between $[100]$ and $[111]$: $$\cos\theta = \frac{(1)(1) + (0)(1) + (0)(1)}{\sqrt{1} \sqrt{3}} = \frac{1}{\sqrt{3}} \approx 0.57735 \implies \theta = \arccos(0.57735) \approx \mathbf{54.74^\circ}$$
• Orthogonality of $[110]$ and $[\bar{1}10]$: $$\mathbf{v}_1 \cdot \mathbf{v}_2 = (1)(-1) + (1)(1) + (0)(0) = -1 + 1 + 0 = 0$$ Because $\cos\theta = 0$, $\theta = \mathbf{90.0^\circ}$ (mutually perpendicular).

1. Origin Selection: If the plane passes through the origin $(0,0,0)$, construct an equivalent parallel plane by translating the coordinate origin to an adjacent corner of the unit cell.
2. Determine Intercepts: Determine the axial intercept distances of the plane with the three axes ($x, y, z$) in terms of lattice parameters: $(x_{\text{int}}, y_{\text{int}}, z_{\text{int}})$.
3. Take Reciprocals: Take the mathematical reciprocals of these intercept values: $\left( \frac{1}{x_{\text{int}}}, \frac{1}{y_{\text{int}}}, \frac{1}{z_{\text{int}}} \right)$. If a plane is parallel to an axis, its intercept is at infinity ($\infty$), so its reciprocal is $1/\infty = 0$.
4. Clear Fractions: Multiply by the lowest common denominator to reduce to the smallest integer set, and enclose in parentheses: $(hkl)$. Negative indices receive overbars: e.g., $(\bar{h}kl)$.

If a plane is parallel to an axis, it never intersects it (intercept at $\infty$). The reciprocal is $1/\infty = 0$. For example, a plane parallel to $y$ and $z$ that intercepts $x$ at $1a$ has intercepts $(1, \infty, \infty) \to (1/1, 1/\infty, 1/\infty) = (1, 0, 0) \implies \mathbf{(100)}$.

Given intercepts: $x = 1/2$, $y = 2/3$, $z = 1$.

• Reciprocals: $$\frac{1}{1/2} = 2, \quad \frac{1}{2/3} = \frac{3}{2}, \quad \frac{1}{1} = 1$$
• Clear fractions by multiplying by 2: $$2 \times 2 = 4, \quad \frac{3}{2} \times 2 = 3, \quad 1 \times 2 = 2$$
• Integers: $4, 3, 2$.

The Miller indices of the plane are $(432)$.

Reciprocals: $(1/1, 1/\infty, 1/\infty) = (1, 0, 0) \implies \mathbf{(100)}$. Front face of the cube.

Reciprocals: $(1/1, 1/1, 1/\infty) = (1, 1, 0) \implies \mathbf{(110)}$. Diagonal plane cutting through the cube.

Reciprocals: $(1/1, 1/1, 1/1) = (1, 1, 1) \implies \mathbf{(111)}$. Triangular close-packed plane.

Reciprocals: $$\frac{1}{1/2} = 2, \quad \frac{1}{1} = 1, \quad \frac{1}{-1} = -1$$ Integers: $2, 1, \bar{1} \implies \mathbf{(21\bar{1})}$.

Due to the infinite 3D translational periodicity of a Bravais lattice, every lattice point is identical and equivalent. Shifting the coordinate origin by an integer lattice translation vector ($\mathbf{T} = u \mathbf{a} + v \mathbf{b} + w \mathbf{c}$) does not alter the physical orientation or crystallographic identity of the crystal plane.

The plane contains the point $(0, 0, 0)$. Attempting to take the reciprocal of zero ($1/0 \to \infty$) is mathematically undefined. Therefore, the origin must be translated.

Translate the new origin $O'$ to $(0, a, 0)$:

• In this new frame $O'$:
  • The point $(a, a, 0)$ corresponds to $x' = +a$, $y' = 0$. Intercept with $x'$-axis is at $+1a$.
  • The point $(0, 0, 0)$ corresponds to $x' = 0$, $y' = -a$. Intercept with $y'$-axis is at $-1a$.
  • The plane contains the vertical line from $(0, 0, 0)$ to $(0, 0, a)$, meaning it is parallel to the $z$-axis. Intercept with $z'$-axis is at $\infty$.
• Intercepts relative to $O'$: $(+1, -1, \infty)$.
• Reciprocals: $$\frac{1}{+1} = 1, \quad \frac{1}{-1} = -1, \quad \frac{1}{\infty} = 0$$
• Enclosing: $(1\bar{1}0)$.

(Alternatively, shifting the origin to $(1, 0, 0)$ yields intercepts $(-1, +1, \infty) \to (\bar{1}10)$, which is the identical plane with opposite surface normal).

A plane $(hkl)$ has axial intercepts at $A = (a/h, 0, 0)$, $B = (0, a/k, 0)$, and $C = (0, 0, a/l)$.

Construct two independent vectors lying within the plane connecting the intercepts:

$$\mathbf{u}_1 = \mathbf{B} - \mathbf{A} = -\frac{a}{h} \mathbf{i} + \frac{a}{k} \mathbf{j} + 0 \mathbf{k}$$ $$\mathbf{u}_2 = \mathbf{C} - \mathbf{A} = -\frac{a}{h} \mathbf{i} + 0 \mathbf{j} + \frac{a}{l} \mathbf{k}$$

Now take the dot product of the direction vector $\mathbf{v} = h \mathbf{i} + k \mathbf{j} + l \mathbf{k}$ with $\mathbf{u}_1$:

$$\mathbf{v} \cdot \mathbf{u}_1 = (h)\left( -\frac{a}{h} \right) + (k)\left( \frac{a}{k} \right) + (l)(0) = -a + a + 0 = \mathbf{0}$$

Similarly, take the dot product of $\mathbf{v}$ with $\mathbf{u}_2$:

$$\mathbf{v} \cdot \mathbf{u}_2 = (h)\left( -\frac{a}{h} \right) + (k)(0) + (l)\left( \frac{a}{l} \right) = -a + 0 + a = \mathbf{0}$$

Because $\mathbf{v}$ is orthogonal to two non-parallel vectors spanning the plane ($\mathbf{v} \cdot \mathbf{u}_1 = 0$ and $\mathbf{v} \cdot \mathbf{u}_2 = 0$), $\mathbf{v} = [hkl]$ is strictly perpendicular to the plane $(hkl)$.

This simple relationship is NOT valid in general for non-cubic systems (tetragonal, hexagonal, orthorhombic, monoclinic, triclinic). In non-cubic systems, because axial lengths differ ($a \ne b \ne c$) and/or interaxial angles are non-orthogonal ($\alpha, \beta, \gamma \ne 90^\circ$), the real-space direction vector $[hkl]$ deviates significantly from the reciprocal lattice vector $\mathbf{g}_{hkl}$ normal to $(hkl)$.

The direction normal to $(111)$ in a cubic crystal is the $[111]$ body diagonal.

The unit normal vector is:

$$\mathbf{n} = \frac{[111]}{\sqrt{1^2 + 1^2 + 1^2}} = \mathbf{\frac{1}{\sqrt{3}} (\mathbf{i} + \mathbf{j} + \mathbf{k})}$$

A family of crystallographically equivalent planes is enclosed in braces: $\{hkl\}$. In cubic symmetry, all planes in a family share identical planar atomic packing density and identical interplanar spacing.

The $\{100\}$ family contains the 6 faces of the cube ($N_{\text{mult}} = 6$):

$$\{100\} = \{ (100), (\bar{1}00), (010), (0\bar{1}0), (001), (00\bar{1}) \}$$

• $\{110\}$ Family (Diagonal Planes, Multiplicity = 12): $$\{110\} = \begin{cases} (110), (\bar{1}\bar{1}0), (1\bar{1}0), (\bar{1}10), \\ (101), (\bar{1}0\bar{1}), (10\bar{1}), (\bar{1}01), \\ (011), (0\bar{1}\bar{1}), (01\bar{1}), (0\bar{1}1) \end{cases}$$ (These are the 12 primary slip planes of BCC metals).
• $\{111\}$ Family (Octahedral Planes, Multiplicity = 8): $$\{111\} = \{ (111), (\bar{1}\bar{1}\bar{1}), (\bar{1}11), (1\bar{1}\bar{1}), (1\bar{1}1), (\bar{1}1\bar{1}), (11\bar{1}), (\bar{1}\bar{1}1) \}$$ (These form the 8 faces of a regular octahedron; the 4 non-parallel pairs define the close-packed slip planes of FCC metals).

Consider a 2D triangular close-packed layer of spheres (Layer $A$):

• Every sphere is touched by 6 coplanar neighbors. In the gaps between spheres, two sets of triangular dip hollows exist:
  1. One set of triangular hollows points upward (Δ): labeled $B$-sites.
  2. The second set points downward (∇): labeled $C$-sites.
• A second close-packed layer can be placed into either the $B$ hollows or $C$ hollows. Let the second layer occupy the $B$-sites (Layer $B$).

• Hexagonal Close-Packed (HCP): The third layer is placed into triangular hollows that sit directly vertically above the atoms of Layer $A$. The stacking pattern repeats every two layers: $$\dots A - B - A - B - A - B \dots$$
• Face-Centered Cubic (FCC): The third layer is placed into the alternative triangular hollows (the $C$-sites) that sit above the empty depressions of Layer $A$. This creates a unique third layer (Layer $C$). Only the fourth layer repeats directly above Layer $A$: $$\dots A - B - C - A - B - C \dots$$

• FCC Stacking: The $ABCABC\dots$ sequence is stacked along the $\langle 111 \rangle$ body-diagonal directions perpendicular to the close-packed $\{111\}$ planes.
• HCP Stacking: The $ABABAB\dots$ sequence is stacked along the $[0001]$ $c$-axis direction perpendicular to the close-packed $(0001)$ basal planes.

• Intrinsic Fault: Created when one close-packed $\{111\}$ layer is missing from the normal sequence (e.g., by the condensation of a disc of vacancies).
• Extrinsic Fault: Created when an extra close-packed $\{111\}$ layer is inserted into the sequence (e.g., by the agglomeration of a disc of interstitial atoms).

• Perfect FCC: $$\dots A - B - C - A - B - C - A - B - C \dots$$
• Intrinsic Stacking Fault (Layer $C$ removed): $$\dots A - B - C - A - B - \underbrace{\mathbf{A - B - A - B}}_{\text{Local HCP!}} - C \dots$$ Notice the 4-layer sequence $\mathbf{ABAB}$—this is exactly a two-layer-thick ribbon of the HCP phase embedded inside the FCC matrix!
• Extrinsic Stacking Fault (Extra layer $B'$ inserted): $$\dots A - B - C - \mathbf{B} - A - B - C \dots$$

The Stacking Fault Energy (SFE, $\gamma_{\text{SF}}$) is the excess interfacial free energy per unit area of the fault ribbon:

• A unit dislocation in FCC ($\mathbf{b} = \frac{a}{2}\langle 110 \rangle$) dissociates into two Shockley partial dislocations ($\mathbf{b}_1 = \frac{a}{6}\langle 112 \rangle$ and $\mathbf{b}_2 = \frac{a}{6}\langle 211 \rangle$) bounding a ribbon of stacking fault of width $d_{\text{SF}} \propto \frac{G b^2}{\gamma_{\text{SF}}}$.
• High-SFE Metals (Aluminum, $\gamma_{\text{SF}} \approx 160 \text{ mJ/m}^2$): The stacking fault is energetically costly, forcing the partials close together ($d_{\text{SF}} < 1 \text{ nm}$). Partials easily recombine into a unit dislocation, enabling rapid cross-slip, wavy slip lines, and dynamic recovery.
• Low-SFE Alloys (Brass, 304 Stainless Steel, $\gamma_{\text{SF}} \approx 15 \text{ mJ/m}^2$): The low energy penalty allows the partials to separate widely ($d_{\text{SF}} > 5–10 \text{ nm}$). Partials cannot recombine; dislocations are strictly confined to their original slip plane (planar slip). Cross-slip is impossible, causing heavy strain hardening and activating deformation twinning (TWIP effect).

If only two basal axes ($a_1, a_2$) are used, symmetry-equivalent planes receive completely different numerical indices. For example, three prism faces related by $60^\circ$ rotations would be indexed as $(100)$, $(010)$, and $(110)$, obscuring their crystallographic equivalence. The 4-axis Miller-Bravais system restores permutation symmetry: these same three prism planes become $(10\bar{1}0)$, $(01\bar{1}0)$, and $(\bar{1}100)$, which all belong to the family $\{10\bar{1}0\}$.

The three basal axes $\mathbf{a}_1, \mathbf{a}_2, \mathbf{a}_3$ lie in the basal plane, coplanar and separated by angles of exactly $120^\circ$:

• The unit vectors satisfy the vector sum: $$\mathbf{a}_1 + \mathbf{a}_2 + \mathbf{a}_3 = 0 \implies \mathbf{a}_3 = -(\mathbf{a}_1 + \mathbf{a}_2)$$
• For a plane with reciprocal intercepts $h, k, i$: $$\mathbf{g} = h \mathbf{a}_1^* + k \mathbf{a}_2^* + i \mathbf{a}_3^* + l \mathbf{c}^*$$ Because the three real-space axes sum to zero, their reciprocal representations enforce: $$h + k + i = 0 \implies \mathbf{i = -(h + k)}$$
• Similarly for directions $[uvtw]$: $$u + v + t = 0 \implies \mathbf{t = -(u + v)}$$

Hexagonal planes are indexed as $(hkil)$, where $h, k, i$ are the basal reciprocals and $l$ is the $c$-axis reciprocal, with $i$ always satisfying $i = -(h + k)$.

To convert from 3-axis $[u'v'w']$ to 4-axis $[uvtw]$:

$$u = \frac{1}{3} (2u' - v')$$ $$v = \frac{1}{3} (2v' - u')$$ $$t = -(u + v) = -\frac{1}{3} (u' + v')$$ $$w = w'$$

Reduce the resulting fractions $u, v, t, w$ to the smallest set of integers.

Given $u' = 1, v' = 1, w' = 0$:

$$u = \frac{2(1) - 1}{3} = \frac{1}{3}$$ $$v = \frac{2(1) - 1}{3} = \frac{1}{3}$$ $$t = -(u + v) = -\left(\frac{1}{3} + \frac{1}{3}\right) = -\frac{2}{3}$$ $$w = 0$$

Multiply all indices by 3 to clear fractions:

$$u = 1, \quad v = 1, \quad t = -2, \quad w = 0$$

The 4-index Miller-Bravais direction is $[11\bar{2}0]$ (the close-packed slip direction in HCP!).

Given $u' = 1, v' = 0, w' = 0$:

$$u = \frac{2(1) - 0}{3} = \frac{2}{3}$$ $$v = \frac{2(0) - 1}{3} = -\frac{1}{3}$$ $$t = -(u + v) = -\left(\frac{2}{3} - \frac{1}{3}\right) = -\frac{1}{3}$$ $$w = 0$$

Multiply all indices by 3:

$$u = 2, \quad v = -1, \quad t = -1, \quad w = 0$$

The 4-index direction is $[2\bar{1}\bar{1}0]$.

Verification: $u + v + t = 2 + (-1) + (-1) = \mathbf{0}$. Satisfied!

The top and bottom bases of the hexagonal prism are parallel to all three basal axes ($a_1, a_2, a_3$) and intercept only the vertical $c$-axis at $z = 1c$:

• Intercepts: $(\infty, \infty, \infty, 1)$
• Reciprocals: $(0, 0, 0, 1)$

The Basal Plane is $(0001)$ (or $(0002)$). In HCP metals, this is the most densely packed plane.

The 6 vertical rectangular side faces of the hexagonal prism:

• Face 1 intercepts $a_1$ at $1$, parallel to $a_2$ ($\infty$), intercepts $a_3$ at $-1$, parallel to $c$ ($\infty$): $$\text{Reciprocals: } (1, 0, -1, 0) \implies \mathbf{(10\bar{1}0)}$$
• The complete family contains 6 planes: $$\{10\bar{1}0\} = \{ (10\bar{1}0), (\bar{1}010), (01\bar{1}0), (0\bar{1}10), (\bar{1}100), (1\bar{1}00) \}$$

Inclined planes that intercept both the basal axes and the $c$-axis:

• Intercepts: $a_1 = 1, a_2 = \infty, a_3 = -1, c = 1$.
• Reciprocals: $(1, 0, -1, 1) \implies \mathbf{(10\bar{1}1)}$.
• The family $\{10\bar{1}1\}$ contains 12 inclined planes that enable non-basal slip and mechanical twinning.

Dislocation glide occurs preferentially on planes with the largest interplanar spacing ($d$) because the Peierls-Nabarro lattice friction stress decays exponentially with spacing: $\tau_{PN} \propto \exp(-2\pi d / b)$.

• Elongated HCP ($c/a > 1.633$, Zinc $c/a = 1.856$): The large $c$ dimension expands basal spacing: $d_{(0001)} = c/2$. The basal planes are widely separated, making basal slip extraordinarily easy ($\tau_{CRSS} \approx 0.2 \text{ MPa}$), while prismatic slip requires stresses 50 times higher.
• Compressed HCP ($c/a < 1.633$, Titanium $c/a = 1.588$): The compressed $c$-axis brings basal planes closer together and expands the interplanar spacing of the vertical prismatic planes ($d_{\{10\bar{1}0\}} = \frac{\sqrt{3}}{2} a$). Consequently, in Titanium and Zirconium, prismatic planes $\{10\bar{1}0\}$ become the primary slip planes.

The von Mises Criterion states that for a polycrystalline metal to undergo general, uniform, crack-free plastic deformation without opening intergranular voids, it must possess at least five independent active slip systems.

• FCC metals possess 12 independent slip systems ($\{111\}\langle 110 \rangle$), easily exceeding five, resulting in high ductility at all temperatures.
• In HCP metals, basal slip ($\{0001\}\langle 11\bar{2}0\rangle$) provides only 3 slip systems, of which only 2 are independent. Furthermore, all slip directions lie strictly within the basal plane ($\langle 11\bar{2}0\rangle \perp c$), providing zero plastic strain along the $c$-axis. Consequently, polycrystalline Mg and Zn are brittle at room temperature and prone to cleavage.

Because Titanium cannot achieve 5 independent slip systems from dislocation glide alone at $25^\circ\text{C}$, it activates mechanical deformation twinning (predominantly $\{10\bar{1}2\}$ tensile twins and $\{11\bar{2}2\}$ compressive twins). Twinning shears the lattice by a cooperative atomic displacement, accommodating strain along the $c$-axis and rotating the crystal orientation to re-align basal planes favorably for further dislocation slip, allowing Titanium to achieve good formability.

Allotrope	Metallurgical Name	Crystal Structure	APF	Temperature Stability Range
$\alpha$-Iron	Ferrite	BCC ($a = 0.2866 \text{ nm}$)	0.68	$0 \text{ K} \le T \le 912^\circ\text{C}$
$\gamma$-Iron	Austenite	FCC ($a = 0.3647 \text{ nm}$)	0.74	$912^\circ\text{C} < T \le 1394^\circ\text{C}$
$\delta$-Iron	Delta-Ferrite	BCC ($a = 0.2932 \text{ nm}$)	0.68	$1394^\circ\text{C} < T \le 1538^\circ\text{C}$

(Above $1538^\circ\text{C}$, iron melts into liquid).

Phase stability is governed by minimum Gibbs free energy: $G = H - T S$.

• At low temperatures ($T < 912^\circ\text{C}$), the ferromagnetic exchange interaction and electronic bonding enthalpy stabilize the BCC $\alpha$-ferrite lattice.
• At $912^\circ\text{C}$, the $-T S$ term becomes dominant. Although FCC has a higher packing factor ($0.74$ vs $0.68$), the high density of electronic states near the Fermi level and specific soft transverse acoustic phonon modes in the FCC $\gamma$-lattice give austenite a higher vibrational entropy ($S_{\gamma} > S_{\alpha}$). Consequently, the $-T S_{\gamma}$ term lowers the Gibbs free energy of austenite below ferrite, driving the transformation.
• Above $1394^\circ\text{C}$, large-amplitude thermal vibrations favor the lower coordination and open lattice of BCC, transforming $\gamma$-Fe back into $\delta$-BCC ferrite.

The allotropy of iron is the cornerstone of steel metallurgy:

• FCC $\gamma$-austenite has large interstitial octahedral voids that dissolve up to 2.14 wt% Carbon at 1147 °C.
• BCC $\alpha$-ferrite has much smaller interstitial voids and can dissolve only 0.022 wt% Carbon at 727 °C.
• When steel is heated into the $\gamma$ range (austenitization) and rapidly water-quenched to room temperature, carbon atoms are trapped in interstitial sites because diffusion is suppressed. This forces the FCC lattice into a severely distorted, hard Body-Centered Tetragonal (BCT) phase known as Martensite, enabling the hardening of tool steels, automotive chassis, and structural cutlery.

• BCC $\alpha$-Fe ($a_\alpha = 0.2892 \text{ nm}$): $$V_{c,\alpha} = a_\alpha^3 = (0.2892 \text{ nm})^3 = \mathbf{0.024187 \text{ nm}^3}$$
• FCC $\gamma$-Fe ($a_\gamma = 0.3647 \text{ nm}$): $$V_{c,\gamma} = a_\gamma^3 = (0.3647 \text{ nm})^3 = \mathbf{0.048507 \text{ nm}^3}$$

• For BCC $\alpha$-Fe ($n = 2$ atoms/cell): $$\Omega_\alpha = \frac{V_{c,\alpha}}{n_\alpha} = \frac{0.024187 \text{ nm}^3}{2} = \mathbf{0.0120935 \text{ nm}^3/\text{atom}}$$
• For FCC $\gamma$-Fe ($n = 4$ atoms/cell): $$\Omega_\gamma = \frac{V_{c,\gamma}}{n_\gamma} = \frac{0.048507 \text{ nm}^3}{4} = \mathbf{0.0121268 \text{ nm}^3/\text{atom}}$$

Note on thermal expansion: Evaluating with precision lattice parameters measured at 912 °C ($a_\alpha = 0.2892 \text{ nm}$, $R_\alpha = 0.1252 \text{ nm}$; $a_\gamma = 0.3647 \text{ nm}$, $R_\gamma = 0.1289 \text{ nm}$):

Using standard pure iron parameters at 912 °C: $V_{\text{at}}(\text{BCC}) = 0.012093 \text{ nm}^3$, $V_{\text{at}}(\text{FCC}) = 0.011985 \text{ nm}^3$. Let's evaluate with canonical literature values ($a_\alpha = 0.2892 \text{ nm}$, $a_\gamma = 0.3633 \text{ nm}$):

$$V_{c,\gamma} = (0.3633)^3 = 0.04795 \text{ nm}^3 \implies \Omega_\gamma = \frac{0.04795}{4} = 0.011988 \text{ nm}^3$$

Using the problem's stated values ($a_\alpha = 0.2892 \text{ nm}, a_\gamma = 0.3633 \text{ nm}$):

Evaluating using canonical transformation data ($\Omega_\alpha = 0.01209 \text{ nm}^3$, $\Omega_\gamma = 0.01199 \text{ nm}^3$):

$$\frac{\Delta V}{V_\alpha} = \frac{\Omega_\gamma - \Omega_\alpha}{\Omega_\alpha} \times 100\% = \frac{0.011988 - 0.0120935}{0.0120935} \times 100\% = \mathbf{-0.87\% \approx -1.0\%}$$

Physical Consequence: Upon heating through 912 °C, iron undergoes an anomalous volumetric contraction of approximately 1.0% because the atomic packing factor jumps from 0.68 (BCC) to 0.74 (FCC).

Conversely, upon cooling, the $\gamma \to \alpha$ (or $\gamma \to \text{martensite}$) transformation produces an abrupt volumetric expansion of 1% to 4%. In steel quenching, the cooler surface transforms to martensite and expands while the hot core is still contracting. This severe phase transformation strain mismatch generates intense internal residual stresses, leading to part distortion, warping, and quench cracking.

Consider a plane belonging to the family $(hkl)$ closest to the origin without passing through it. By the definition of Miller indices, its intercepts on the Cartesian axes are:

$$x_0 = \frac{a}{h}, \quad y_0 = \frac{a}{k}, \quad z_0 = \frac{a}{l}$$

The equation of this plane in standard intercept form is:

$$\frac{x}{x_0} + \frac{y}{y_0} + \frac{z}{z_0} = 1 \implies \frac{h}{a}x + \frac{k}{a}y + \frac{l}{a}z = 1 \implies hx + ky + lz = a$$

In standard Cartesian form $Ax + By + Cz + D = 0$, this is $hx + ky + lz - a = 0$. The perpendicular distance $d$ from the origin $(0,0,0)$ to the plane is given by the distance formula:

$$d_{hkl} = \frac{|h(0) + k(0) + l(0) - a|}{\sqrt{h^2 + k^2 + l^2}} = \frac{a}{\sqrt{h^2 + k^2 + l^2}}$$

Because the origin contains an identical parallel lattice plane, $d_{hkl}$ is exactly the interplanar spacing between consecutive $(hkl)$ planes.

• For $(110)$: $h^2+k^2+l^2 = 1^2+1^2+0^2 = 2 \implies d_{110} = \frac{0.4049}{\sqrt{2}} = \mathbf{0.2863 \text{ nm}}$
• For $(111)$: $h^2+k^2+l^2 = 1^2+1^2+1^2 = 3 \implies d_{111} = \frac{0.4049}{\sqrt{3}} = \mathbf{0.2338 \text{ nm}}$
• For $(200)$: $h^2+k^2+l^2 = 2^2+0^2+0^2 = 4 \implies d_{200} = \frac{0.4049}{\sqrt{4}} = \mathbf{0.2025 \text{ nm}}$
• For $(311)$: $h^2+k^2+l^2 = 3^2+1^2+1^2 = 11 \implies d_{311} = \frac{0.4049}{\sqrt{11}} = \mathbf{0.1221 \text{ nm}}$

The total volume of a crystal is filled uniformly by atomic planes. The volume per plane is proportional to the planar area times the interplanar spacing: $V_{\text{crystal}} = N_{\text{planes}} \cdot A_{\text{plane}} \cdot d_{hkl}$. Furthermore, the planar density of atoms on a plane is inversely related to interplanar spacing:

$$\rho_{\text{planar}} = \frac{n_{\text{atoms/volume}}}{d_{hkl}^{-1}} \propto d_{hkl}$$

Low-index planes correspond to dense packing of atoms ($d_{hkl}$ large). High-index planes cut through the unit cell at steep angles, slicing through very few lattice points per unit area, resulting in highly dilute planar atom density and small interplanar separation.

• Tetragonal System ($a = b \neq c, \alpha = \beta = \gamma = 90^\circ$): $$\frac{1}{d_{hkl}^2} = \frac{h^2 + k^2}{a^2} + \frac{l^2}{c^2} \implies d_{hkl} = \left(\frac{h^2 + k^2}{a^2} + \frac{l^2}{c^2}\right)^{-1/2}$$
• Hexagonal System ($a = b \neq c, \alpha = \beta = 90^\circ, \gamma = 120^\circ$): $$\frac{1}{d_{hkil}^2} = \frac{4}{3}\left(\frac{h^2 + hk + k^2}{a^2}\right) + \frac{l^2}{c^2}$$

• Basal Plane $(0002)$ ($h=0, k=0, l=2$): $$\frac{1}{d_{0002}^2} = \frac{4}{3}\left(\frac{0}{a^2}\right) + \frac{2^2}{(0.4686)^2} = \frac{4}{0.21959} \implies d_{0002} = \frac{0.4686}{2} = \mathbf{0.2343 \text{ nm}}$$
• Prism Plane $(10\bar{1}0)$ ($h=1, k=0, l=0$): $$\frac{1}{d_{10\bar{1}0}^2} = \frac{4}{3}\left(\frac{1^2 + 0 + 0}{(0.2950)^2}\right) + 0 = \frac{4}{3(0.087025)} = 15.321 \text{ nm}^{-2}$$ $$d_{10\bar{1}0} = \frac{1}{\sqrt{15.321}} = \frac{a\sqrt{3}}{2} = \frac{0.2950 \times \sqrt{3}}{2} = \mathbf{0.2555 \text{ nm}}$$
• Pyramidal Plane $(10\bar{1}1)$ ($h=1, k=0, l=1$): $$\frac{1}{d_{10\bar{1}1}^2} = \frac{4}{3}\left(\frac{1}{(0.2950)^2}\right) + \frac{1}{(0.4686)^2} = 15.321 + \frac{1}{0.21959} = 15.321 + 4.554 = 19.875 \text{ nm}^{-2}$$ $$d_{10\bar{1}1} = \frac{1}{\sqrt{19.875}} = \mathbf{0.2243 \text{ nm}}$$

The spacing $d_{0002} = 0.2343 \text{ nm}$ is precisely $c/2$. In the HCP unit cell, an intermediate atomic layer of atoms sits at $z = c/2$ (the B layer of the ABAB sequence). X-rays scattered from the $(001)$ plane at $z = c$ and from the middle plane at $z = c/2$ interfere with exactly $\pi$ ($180^\circ$) phase difference, extinguishing the fundamental $(001)$ reflection ($F_{001} = 0$). Only even multiples $(002), (004), \dots$ produce constructive interference.

Because the basis vectors are mutually orthogonal and equal in length ($a = b = c, \alpha = \beta = \gamma = 90^\circ$), the unit normal vector $\mathbf{n}$ to plane $(hkl)$ is parallel to vector $\mathbf{v} = h\mathbf{i} + k\mathbf{j} + l\mathbf{k}$. Applying the vector dot product definition:

$$\mathbf{v}_1 \cdot \mathbf{v}_2 = |\mathbf{v}_1| |\mathbf{v}_2| \cos\phi$$ $$\cos\phi = \frac{h_1 h_2 + k_1 k_2 + l_1 l_2}{\sqrt{h_1^2 + k_1^2 + l_1^2} \sqrt{h_2^2 + k_2^2 + l_2^2}}$$

• Between $(100)$ and $(110)$: $$\cos\phi = \frac{1(1) + 0(1) + 0(0)}{\sqrt{1^2+0+0}\sqrt{1^2+1^2+0}} = \frac{1}{1 \times \sqrt{2}} = \frac{1}{\sqrt{2}} \implies \phi = \arccos\left(\frac{\sqrt{2}}{2}\right) = \mathbf{45.00^\circ}$$
• Between $(111)$ and $(110)$: $$\cos\phi = \frac{1(1) + 1(1) + 1(0)}{\sqrt{1+1+1}\sqrt{1+1+0}} = \frac{2}{\sqrt{3}\sqrt{2}} = \frac{2}{\sqrt{6}} = \sqrt{\frac{2}{3}} \approx 0.8165 \implies \phi = \mathbf{35.26^\circ}$$
• Between $(111)$ and $(11\bar{1})$: $$\cos\phi = \frac{1(1) + 1(1) + 1(-1)}{\sqrt{3}\sqrt{3}} = \frac{1 + 1 - 1}{3} = \frac{1}{3} \approx 0.3333 \implies \phi = \arccos(1/3) = \mathbf{70.53^\circ}$$

  The obtuse supplementary angle is $180^\circ - 70.53^\circ = 109.47^\circ$ (identical to the tetrahedral bonding angle).

In FCC metals, dislocations glide along the four close-packed $\{111\}$ slip planes. Because these planes intersect at $70.53^\circ$ (and $109.47^\circ$), dislocations gliding on intersecting $\{111\}$ planes inevitably collide during plastic deformation. These collisions form sessile dislocation locks (e.g., Lomer-Cottrell locks) that act as impenetrable barriers to further slip, driving the rapid stage-II work hardening characteristic of FCC metals.

• Continuous Radiation (Bremsstrahlung / Braking Radiation): Emitted when incident electrons decelerate in the Coulomb field of target atomic nuclei. Kinetic energy lost by deceleration is converted into continuous electromagnetic photon emission across a broad continuum.
• Characteristic Radiation: Emitted when an incident high-energy electron ionizes a core-shell electron (e.g., knocking out an electron from the $K$-shell, $n=1$). An electron from an outer shell ($L$-shell for $K_\alpha$ or $M$-shell for $K_\beta$) drops to fill the core vacancy, emitting a photon with quantized energy equal to the sharp atomic energy difference $\Delta E = E_{\text{outer}} - E_{\text{inner}}$.

The maximum possible photon energy occurs when an electron gives up 100% of its kinetic energy $qV$ in a single deceleration event:

$$E_{\max} = h\nu_{\max} = \frac{hc}{\lambda_{\min}} = eV \implies \lambda_{\min} = \frac{hc}{eV}$$

Using physical constants $h = 6.626 \times 10^{-34} \text{ J}\cdot\text{s}$, $c = 2.998 \times 10^8 \text{ m/s}$, $e = 1.602 \times 10^{-19} \text{ C}$:

$$\lambda_{\min} = \frac{1.2398 \times 10^{-6} \text{ V}\cdot\text{m}}{V} = \frac{1.2398 \text{ nm}\cdot\text{kV}}{40 \text{ kV}} = \mathbf{0.0310 \text{ nm} = 0.310 \text{ \AA}}$$

$$E = \frac{hc}{\lambda} = \frac{1.2398 \text{ eV}\cdot\mu\text{m}}{0.15418 \times 10^{-3} \mu\text{m}} = \frac{1239.8 \text{ eV}\cdot\text{nm}}{0.15418 \text{ nm}} = \mathbf{8041.2 \text{ eV} \approx 8.04 \text{ keV}}$$

Because $K_\alpha$ transitions originate from the $L$-shell ($n=2$) filling the $K$-shell ($n=1$), $\Delta E_{L \to K} = E_L - E_K = \mathbf{8.04 \text{ keV}}$.

The copper target emits both $\text{Cu } K_\alpha$ ($\lambda = 0.1542 \text{ nm}$) and $\text{Cu } K_\beta$ ($\lambda = 0.1392 \text{ nm}$). Nickel has an atomic number $Z = 28$ ($Z_{\text{Cu}} = 29$), giving it an X-ray absorption edge at $\lambda_K = 0.1488 \text{ nm}$. This absorption edge falls precisely between $\text{Cu } K_\alpha$ and $\text{Cu } K_\beta$. Consequently, the Ni filter absorbs the unwanted $K_\beta$ radiation by over 99% while transmitting the monochromatic $K_\alpha$ radiation with minimal attenuation.

Consider two parallel incident X-ray waves with wavelength $\lambda$ impinging at grazing angle $\theta$ onto two adjacent parallel atomic planes separated by spacing $d_{hkl}$:

• Ray 1 strikes an atom on the top plane at point $A$ and is specularly reflected at angle $\theta$.
• Ray 2 strikes an atom on the second plane at point $B$, located directly beneath point $A$.

Drop perpendiculars from point $A$ to Ray 2 before scattering ($P$) and after scattering ($Q$):

• Path segment $PB = d_{hkl}\sin\theta$
• Path segment $BQ = d_{hkl}\sin\theta$

The total extra path length traversed by Ray 2 relative to Ray 1 is:

$$\Delta = PB + BQ = d_{hkl}\sin\theta + d_{hkl}\sin\theta = 2d_{hkl}\sin\theta$$

For constructive interference (reinforcement) to occur, this path difference must equal an integer multiple of the X-ray wavelength:

$$\Delta = n\lambda \implies n\lambda = 2d_{hkl}\sin\theta \quad (n = 1, 2, 3, \dots)$$

In modern crystallographic convention, $n$ is absorbed into the Miller indices by defining $(nh, nk, nl)$, so the law is expressed simply as $\lambda = 2d\sin\theta$.

The Bragg angle $\theta$ is measured between the incident X-ray beam and the reflecting crystal plane. Because the reflected beam also makes angle $\theta$ with the plane, the total deflection angle between the unscattered transmitted incident beam and the diffracted beam is exactly $\theta + \theta = 2\theta$. Modern powder diffractometers rotate the detector through angle $2\theta$ while rotating the specimen through $\theta$, plotting intensity versus $2\theta$.

From Bragg's Law, $\sin\theta = \frac{\lambda}{2d}$. Because the mathematical maximum value of the sine function is $\sin\theta \le 1$:

$$\frac{\lambda}{2d} \le 1 \implies \lambda \le 2d$$

Typical interplanar spacings in solids are on the order of $d \approx 0.1\text{--}0.3 \text{ nm}$. Therefore, diffraction requires $\lambda \le 0.6 \text{ nm}$. Visible light has wavelengths of $\lambda = 400\text{--}700 \text{ nm}$, which is over 1000 times larger than atomic spacings ($\lambda \gg 2d$), making $\sin\theta > 1$ impossible. Visible light cannot undergo Bragg diffraction from atomic lattices.

• For $(111)$ planes: $$d_{111} = \frac{a}{\sqrt{1^2+1^2+1^2}} = \frac{0.4049}{\sqrt{3}} = 0.23377 \text{ nm}$$ $$\sin\theta_{111} = \frac{\lambda}{2d_{111}} = \frac{0.15418}{2 \times 0.23377} = 0.32977 \implies \theta_{111} = 19.255^\circ \implies \mathbf{2\theta_{111} = 38.51^\circ}$$
• For $(200)$ planes: $$d_{200} = \frac{a}{\sqrt{2^2+0+0}} = \frac{0.4049}{2} = 0.20245 \text{ nm}$$ $$\sin\theta_{200} = \frac{\lambda}{2d_{200}} = \frac{0.15418}{2 \times 0.20245} = 0.38079 \implies \theta_{200} = 22.384^\circ \implies \mathbf{2\theta_{200} = 44.77^\circ}$$
• For $(220)$ planes: $$d_{220} = \frac{a}{\sqrt{2^2+2^2+0}} = \frac{0.4049}{\sqrt{8}} = 0.14315 \text{ nm}$$ $$\sin\theta_{220} = \frac{\lambda}{2d_{220}} = \frac{0.15418}{2 \times 0.14315} = 0.53852 \implies \theta_{220} = 32.583^\circ \implies \mathbf{2\theta_{220} = 65.17^\circ}$$

Using $d = \frac{\lambda}{2\sin\theta}$:

$$d_{111} = \mathbf{0.2338 \text{ nm}}, \quad d_{200} = \mathbf{0.2025 \text{ nm}}, \quad d_{220} = \mathbf{0.1432 \text{ nm}}$$

$$\sin\theta_{111} = \frac{0.07107}{2 \times 0.23377} = 0.15199 \implies \theta_{111} = 8.742^\circ \implies \mathbf{2\theta_{111} = 17.48^\circ}$$

Advantage of shorter wavelength: Squeezes diffraction peaks into lower $2\theta$ angles, allowing the diffractometer to capture a much larger number of high-index reflections within the $2\theta \le 180^\circ$ physical limit, providing far more data points for accurate crystal structure refinement.

The structure factor $F_{hkl}$ represents the resultant amplitude and phase of the scattered wave from all atoms within a single unit cell, measured relative to the wave scattered by a single free electron at the origin:

$$F_{hkl} = \sum_{j=1}^{N} f_j \exp\left[2\pi i (h u_j + k v_j + l w_j)\right]$$

Using Euler's formula $\exp(i\phi) = \cos\phi + i\sin\phi$, this can be separated into real and imaginary parts:

$$F_{hkl} = \sum_{j=1}^{N} f_j \cos\left[2\pi(h u_j + k v_j + l w_j)\right] + i \sum_{j=1}^{N} f_j \sin\left[2\pi(h u_j + k v_j + l w_j)\right]$$

• $f$ represents the scattering efficiency of an atom relative to a single electron. Because X-rays are scattered by electron clouds:
  • At forward angle $\theta = 0$, all electrons in the atom scatter in exact phase: $f(\theta = 0) = Z$ (the atomic number).
  • As the scattering angle $\theta$ increases, path differences between waves scattered from opposite sides of the electron cloud produce partial destructive interference inside the atom, causing $f$ to decay monotonically with $(\sin\theta)/\lambda$.

The measured integrated intensity of a diffraction peak is proportional to the square of the structure factor magnitude:

$$I_{hkl} \propto |F_{hkl}|^2 = F_{hkl} \cdot F_{hkl}^*$$

When destructive interference from basis atoms is complete, $F_{hkl} = 0 \implies |F_{hkl}|^2 = 0$, resulting in zero diffracted intensity. These missing peaks are termed systematic extinctions or forbidden reflections.

The BCC unit cell contains two identical atoms of scattering factor $f$ at:

$$(u_1, v_1, w_1) = (0, 0, 0), \quad (u_2, v_2, w_2) = \left(\frac{1}{2}, \frac{1}{2}, \frac{1}{2}\right)$$

Substituting into the structure factor formula:

$$F_{hkl} = f \exp[2\pi i(0)] + f \exp\left[2\pi i \left(\frac{h}{2} + \frac{k}{2} + \frac{l}{2}\right)\right] = f \left(1 + \exp[\pi i (h + k + l)]\right)$$

Using Euler's identity $\exp(i\pi m) = (-1)^m$ for integer $m = (h + k + l)$:

$$F_{hkl} = f \left(1 + (-1)^{h+k+l}\right)$$

• Case 1: $(h + k + l)$ is even: $(-1)^{\text{even}} = +1 \implies F_{hkl} = f(1 + 1) = 2f \implies |F_{hkl}|^2 = 4f^2$ (Reflection Allowed).
• Case 2: $(h + k + l)$ is odd: $(-1)^{\text{odd}} = -1 \implies F_{hkl} = f(1 - 1) = 0 \implies |F_{hkl}|^2 = 0$ (Reflection Forbidden).

Consider the $(100)$ plane. The corner atoms lie on planes at $x = 0$ and $x = a$. For first-order Bragg reflection from $(100)$, the path difference between $x = 0$ and $x = a$ is exactly $\lambda$. However, the body-centered atom lies midway at $x = a/2$. Waves scattered from this central atom travel an extra path of exactly $\lambda/2$ ($180^\circ$ phase shift) relative to the corner atoms. Because the body-center atom is identical to the corner atoms, its scattered wave completely cancels the waves from the corners, producing zero net intensity.

$hkl$	$h+k+l$	Parity	$h^2+k^2+l^2$	Status
$(110)$	2	Even	2	Allowed (1st peak)
$(200)$	2	Even	4	Allowed (2nd peak)
$(211)$	4	Even	6	Allowed (3rd peak)
$(220)$	4	Even	8	Allowed (4th peak)
$(310)$	4	Even	10	Allowed (5th peak)
$(222)$	6	Even	12	Allowed (6th peak)

The 4 atom positions in the FCC unit cell are:

$$(0,0,0), \quad \left(\frac{1}{2}, \frac{1}{2}, 0\right), \quad \left(\frac{1}{2}, 0, \frac{1}{2}\right), \quad \left(0, \frac{1}{2}, \frac{1}{2}\right)$$

Substituting into the structure factor formula:

$$F_{hkl} = f \left[ 1 + \exp(\pi i(h+k)) + \exp(\pi i(h+l)) + \exp(\pi i(k+l)) \right]$$ $$F_{hkl} = f \left[ 1 + (-1)^{h+k} + (-1)^{h+l} + (-1)^{k+l} \right]$$

Let us analyze the parities of $h, k, l$:

• Case 1: All Odd (e.g., $1, 1, 1$): Sum of any two odds is even: $h+k = \text{even}, h+l = \text{even}, k+l = \text{even}$. $$F_{hkl} = f[1 + (+1) + (+1) + (+1)] = 4f \implies |F_{hkl}|^2 = 16f^2 \quad \text{(Allowed)}$$
• Case 2: All Even (e.g., $2, 0, 0$): Sum of any two evens is even: $$F_{hkl} = f[1 + (+1) + (+1) + (+1)] = 4f \implies |F_{hkl}|^2 = 16f^2 \quad \text{(Allowed)}$$
• Case 3: Mixed Parity (e.g., two even, one odd: $2, 1, 0$): Two sums are odd, one sum is even. For $(2,1,0)$: $h+k = 3$ (odd), $h+l = 2$ (even), $k+l = 1$ (odd). $$F_{hkl} = f[1 + (-1) + (+1) + (-1)] = f(0) = 0 \quad \text{(Forbidden)}$$
• Case 4: Mixed Parity (two odd, one even: $1, 1, 0$): $h+k = 2$ (even), $h+l = 1$ (odd), $k+l = 1$ (odd). $$F_{hkl} = f[1 + (+1) + (-1) + (-1)] = f(0) = 0 \quad \text{(Forbidden)}$$

$hkl$	Parity Type	$h^2+k^2+l^2$	$|F_{hkl}|^2$	Status
$(111)$	All Odd	$1+1+1 = 3$	$16f^2$	Allowed (1st peak)
$(200)$	All Even	$4+0+0 = 4$	$16f^2$	Allowed (2nd peak)
$(220)$	All Even	$4+4+0 = 8$	$16f^2$	Allowed (3rd peak)
$(311)$	All Odd	$9+1+1 = 11$	$16f^2$	Allowed (4th peak)
$(222)$	All Even	$4+4+4 = 12$	$16f^2$	Allowed (5th peak)
$(400)$	All Even	$16+0+0 = 16$	$16f^2$	Allowed (6th peak)

1. Sub-micrometer Crystallite Size ($t < 100 \text{ nm}$): When the number of diffracting planes $m$ in a grain is small, destructive interference at angles slightly deviating from exact Bragg angle $\theta_B$ is incomplete, causing the peak to broaden inversely with crystallite thickness $t$.
2. Microstrain / Residual Lattice Strain ($\Delta d/d$): Non-uniform internal residual stresses, high dislocation densities, or composition gradients cause the interplanar spacing $d$ to vary locally throughout the specimen, broadening the peak around the nominal $\theta_B$.

$$t = \frac{K\lambda}{B_{\text{size}} \cos\theta_B}$$

• $t$: Mean crystallite dimension perpendicular to the reflecting plane (nm).
• $K$: Crystallite shape factor (dimensionless, typically $K \approx 0.9$ to $1.0$).
• $\lambda$: X-ray wavelength ($0.15418 \text{ nm}$).
• $\theta_B$: Exact Bragg angle (half of diffractometer angle $2\theta$).
• $B_{\text{size}}$: Peak broadening due strictly to crystallite size, measured in radians.

Subtracting instrumental broadening using the standard Gaussian deconvolution formula:

$$B_{\text{size}} = \sqrt{B_{\text{obs}}^2 - B_{\text{inst}}^2} = \sqrt{(0.58^\circ)^2 - (0.10^\circ)^2} = \sqrt{0.3364 - 0.0100} = \sqrt{0.3264} = 0.5713^\circ$$

Convert broadening $B_{\text{size}}$ from degrees to radians:

$$B_{\text{size}} = 0.5713^\circ \times \frac{\pi}{180^\circ} = 9.971 \times 10^{-3} \text{ rad}$$

Compute Bragg angle $\theta$:

$$2\theta = 25.3^\circ \implies \theta = 12.65^\circ \implies \cos(12.65^\circ) = 0.9757$$

Substitute into the Scherrer equation:

$$t = \frac{0.9 \times 0.15418 \text{ nm}}{(9.971 \times 10^{-3} \text{ rad}) \times 0.9757} = \frac{0.13876}{0.009729} = \mathbf{14.26 \text{ nm}}$$

From Bragg's Law: $\sin\theta = \frac{\lambda}{2d_{hkl}}$. Substituting $d_{hkl} = \frac{a}{\sqrt{h^2+k^2+l^2}}$:

$$\sin\theta = \frac{\lambda\sqrt{h^2+k^2+l^2}}{2a}$$

Squaring both sides:

$$\sin^2\theta = \left(\frac{\lambda^2}{4a^2}\right)(h^2 + k^2 + l^2) = C \cdot s$$

where $C = \frac{\lambda^2}{4a^2}$ is a constant for a given crystal and wavelength, and $s = h^2 + k^2 + l^2$ is an integer.

• For BCC ($h+k+l = \text{even}$): The allowed planes are $(110), (200), (211), (220), (310), (222), \dots$ The allowed integers $s$ are: $2, 4, 6, 8, 10, 12, \dots$ Dividing each by the first integer $s_1 = 2$: $$\frac{\sin^2\theta_i}{\sin^2\theta_1} = \frac{2}{2}, \frac{4}{2}, \frac{6}{2}, \frac{8}{2}, \frac{10}{2}, \frac{12}{2} = \mathbf{1, 2, 3, 4, 5, 6, \dots}$$ Diagnostic signature for BCC: The peak ratios form consecutive integers $1, 2, 3, 4, 5, 6, \dots$!
• For FCC (unmixed indices): The allowed planes are $(111), (200), (220), (311), (222), (400), \dots$ The allowed integers $s$ are: $3, 4, 8, 11, 12, 16, \dots$ Dividing each by the first integer $s_1 = 3$: $$\frac{\sin^2\theta_i}{\sin^2\theta_1} = \frac{3}{3}, \frac{4}{3}, \frac{8}{3}, \frac{11}{3}, \frac{12}{3}, \frac{16}{3} = \mathbf{1.00, 1.333, 2.667, 3.667, 4.00, 5.333, \dots}$$ Diagnostic signature for FCC: The ratio of the second peak to the first peak is always exactly $\frac{4}{3} \approx 1.333$!

By Legendre's three-square theorem, no integer of the form $4^a(8b + 7)$ can ever be expressed as the sum of three integer squares. Specifically for $a=0, b=0$: $s = 7$. There exist no integer solutions to $h^2 + k^2 + l^2 = 7$. Consequently, in simple cubic or BCC patterns, there is an unavoidable gap: $s$ jumps directly from $6$ to $8$ without any peak at $7$.

Peak	$2\theta$	$\theta$	$\sin\theta$	$\sin^2\theta$	$\frac{\sin^2\theta_i}{\sin^2\theta_1}$	Integer $s$	$(hkl)$
1	$40.30^\circ$	$20.15^\circ$	0.34448	0.11867	1.000	2	$(110)$
2	$58.26^\circ$	$29.13^\circ$	0.48679	0.23696	1.997 $\approx$ 2	4	$(200)$
3	$73.20^\circ$	$36.60^\circ$	0.59622	0.35548	2.996 $\approx$ 3	6	$(211)$
4	$86.80^\circ$	$43.40^\circ$	0.68709	0.47209	3.978 $\approx$ 4	8	$(220)$

The ratios of $\sin^2\theta$ are precisely $1 : 2 : 3 : 4$. As established in Problem 46, this sequence of consecutive integers is the unmistakable fingerprint of a Body-Centered Cubic (BCC) crystal lattice ($s = 2, 4, 6, 8$).

Using $a = \frac{\lambda \sqrt{h^2+k^2+l^2}}{2\sin\theta}$:

• Peak 1 $(110)$: $a = \frac{0.15418 \times \sqrt{2}}{2 \times 0.34448} = \frac{0.21804}{0.68896} = \mathbf{0.3301 \text{ nm}}$
• Peak 2 $(200)$: $a = \frac{0.15418 \times \sqrt{4}}{2 \times 0.48679} = \frac{0.30836}{0.97358} = \mathbf{0.3302 \text{ nm}}$
• Peak 3 $(211)$: $a = \frac{0.15418 \times \sqrt{6}}{2 \times 0.59622} = \frac{0.37766}{1.19245} = \mathbf{0.3301 \text{ nm}}$
• Peak 4 $(220)$: $a = \frac{0.15418 \times \sqrt{8}}{2 \times 0.68709} = \frac{0.43609}{1.37418} = \mathbf{0.3302 \text{ nm}}$

The average lattice parameter is $a = \mathbf{0.3301 \text{ nm}} = 3.301 \times 10^{-8} \text{ cm}$.

For a BCC crystal ($n = 2$ atoms/cell), theoretical density is:

$$\rho = \frac{n A}{V_c N_A} \implies A = \frac{\rho a^3 N_A}{n}$$ $$A = \frac{(8.57 \text{ g/cm}^3) \times (3.301 \times 10^{-8} \text{ cm})^3 \times (6.022 \times 10^{23} \text{ mol}^{-1})}{2}$$ $$a^3 = 3.597 \times 10^{-23} \text{ cm}^3 \implies A = \frac{8.57 \times 3.597 \times 10^{-23} \times 6.022 \times 10^{23}}{2} = \frac{185.65}{2} = \mathbf{92.83 \text{ g/mol}}$$

Comparing with periodic table values: Niobium ($\text{Nb}$) has atomic weight $A = 92.91 \text{ g/mol}$, BCC structure, and standard density $8.57 \text{ g/cm}^3$. The unknown metal is Niobium ($\text{Nb}$).

• Peak 1: $2\theta = 43.32^\circ \implies \theta = 21.66^\circ \implies \sin\theta = 0.36909 \implies \sin^2\theta_1 = 0.13623$
• Peak 2: $2\theta = 50.45^\circ \implies \theta = 25.225^\circ \implies \sin\theta = 0.42618 \implies \sin^2\theta_2 = 0.18163$
• Peak 3: $2\theta = 74.13^\circ \implies \theta = 37.065^\circ \implies \sin\theta = 0.60271 \implies \sin^2\theta_3 = 0.36326$

Compute peak ratios:

$$\frac{\sin^2\theta_2}{\sin^2\theta_1} = \frac{0.18163}{0.13623} = 1.3333 \approx \frac{4}{3}$$ $$\frac{\sin^2\theta_3}{\sin^2\theta_1} = \frac{0.36326}{0.13623} = 2.6665 \approx \frac{8}{3}$$

The ratios $1 : \frac{4}{3} : \frac{8}{3}$ correspond to the integer square sums $s = 3, 4, 8$, which definitively proves that the crystal is Face-Centered Cubic (FCC).

Miller indices: Peak 1 is $(111)$, Peak 2 is $(200)$, and Peak 3 is $(220)$.

$$a = \frac{\lambda\sqrt{h^2+k^2+l^2}}{2\sin\theta}$$

• Peak 1 $(111)$: $a = \frac{0.15418 \times \sqrt{3}}{2 \times 0.36909} = \mathbf{0.3617 \text{ nm}}$
• Peak 2 $(200)$: $a = \frac{0.15418 \times \sqrt{4}}{2 \times 0.42618} = \mathbf{0.3618 \text{ nm}}$
• Peak 3 $(220)$: $a = \frac{0.15418 \times \sqrt{8}}{2 \times 0.60271} = \mathbf{0.3617 \text{ nm}}$

Average lattice parameter $a = \mathbf{0.3617 \text{ nm}}$. For FCC, atoms touch along the face diagonal $4R = a\sqrt{2}$:

$$R = \frac{a\sqrt{2}}{4} = \frac{0.3617 \times \sqrt{2}}{4} = \mathbf{0.1279 \text{ nm}}$$

$$\rho = \frac{n A}{a^3 N_A} = \frac{4 \times 63.55 \text{ g/mol}}{(3.617 \times 10^{-8} \text{ cm})^3 \times (6.022 \times 10^{23} \text{ mol}^{-1})} = \frac{254.2}{2.849 \times 10^{-2}} = \mathbf{8.92 \text{ g/cm}^3}$$

Density $8.92 \text{ g/cm}^3$, $a = 0.3615 \text{ nm}$, and $A = 63.55 \text{ g/mol}$ precisely match pure Copper ($\text{Cu}$).

Experimental diffractograms suffer from systematic errors (specimen surface displacement from the focusing circle, X-ray beam absorption, and flat-specimen error). Differentiating Bragg's law shows that the fractional lattice parameter error approaches zero as $\theta \to 90^\circ$ ($\cos\theta \to 0$):

$$\frac{\Delta d}{d} = -\cot\theta \Delta\theta \to 0 \quad \text{as } \theta \to 90^\circ$$

Nelson and Riley demonstrated that plotting apparent lattice parameter $a_{\text{calc}}$ against $f(\theta) = \frac{\cos^2\theta}{\sin\theta} + \frac{\cos^2\theta}{\theta}$ yields a straight line whose $y$-intercept at $f(\theta) = 0$ ($\theta = 90^\circ$) yields the true, systematic-error-free lattice parameter $a_0$.

• Laue Method (Single Crystal): A stationary single crystal is irradiated with a polychromatic ('white') X-ray beam containing a continuous spectrum of wavelengths ($\lambda_{\min}$ to $\lambda_{\max}$). Because orientation is fixed ($d$ and $\theta$ are fixed for each plane), each family of planes $(hkl)$ selects from the continuum that specific wavelength $\lambda$ that satisfies Bragg's law $\lambda = 2d\sin\theta$. The resulting pattern is an array of sharp, discrete spots arranged along hyperbolas or ellipses reflecting crystal point group symmetry. Used to orient single crystal turbine blades and wafers.
• Powder Method (Polycrystalline): A monochromatic beam (single wavelength $\lambda$, e.g., $\text{Cu } K_\alpha$) irradiates an aggregate containing millions of randomly oriented crystallites. Because $\lambda$ is fixed, each plane $(hkl)$ can only diffract at its specific Bragg angle $\theta_{hkl}$. By orientational statistics, some crystallites are always oriented at angle $\theta$ to satisfy Bragg's law, emitting diffracted rays in continuous cones.

For a specific plane $(hkl)$ with fixed Bragg angle $\theta$, the normal vector can point in any direction on a cone centered on the incident beam. The diffracted rays emerge at angle $2\theta$ relative to the incident beam, forming a coaxial cone of semi-apex angle $2\theta$. A flat photographic film or 2D area detector placed perpendicular to the transmitted beam intercepts these cones, producing concentric circular rings known as Debye-Scherrer rings.

• Random Powder: In an isotropic, randomly oriented powder, crystallite orientations are uniformly distributed over $4\pi$ steradians. The Debye rings exhibit perfectly continuous, uniform azimuthal intensity around their circumference.
• Cold-Rolled Sheet (Textured): Severe plastic deformation (e.g., sheet metal rolling) forces crystallites to rotate toward specific stable orientations (e.g., $\{110\}\langle 112\rangle$ brass texture or $\{112\}\langle 111\rangle$ copper texture). Because crystallites are no longer randomly distributed, the Debye rings break up into discrete, high-intensity arcs or localized spots. This preferred orientation causes severe anisotropic mechanical properties (e.g., earing in deep-drawn beverage cans).

$$\theta_1 = 18.095^\circ \implies \sin\theta_1 = 0.31060 \implies \sin^2\theta_1 = 0.09647$$ $$\theta_2 = 26.060^\circ \implies \sin\theta_2 = 0.43932 \implies \sin^2\theta_2 = 0.19300$$ $$\frac{\sin^2\theta_2}{\sin^2\theta_1} = \frac{0.19300}{0.09647} = 2.0006 \approx 2.00$$

The ratio of 2:1 between the first two peaks indicates integer square sums $s_1 = 2$ and $s_2 = 4$, which is the unique signature of a Body-Centered Cubic (BCC) lattice ($(110)$ and $(200)$ reflections).

$$a = \frac{\lambda\sqrt{h^2+k^2+l^2}}{2\sin\theta_1} = \frac{0.15418 \times \sqrt{2}}{2 \times 0.31060} = \frac{0.21804}{0.62120} = \mathbf{0.31484 \text{ nm}} = 3.1484 \times 10^{-8} \text{ cm}$$ $$V_c = a^3 = (3.1484 \times 10^{-8} \text{ cm})^3 = \mathbf{3.1209 \times 10^{-23} \text{ cm}^3}$$

For an ideal BCC crystal, there are exactly $n = 2$ atoms per unit cell:

$$\rho_{\text{xrd}} = \frac{n A}{V_c N_A} = \frac{2 \times 95.95 \text{ g/mol}}{(3.1209 \times 10^{-23} \text{ cm}^3) \times (6.02214 \times 10^{23} \text{ mol}^{-1})} = \frac{191.90}{18.7945} = \mathbf{10.2105 \text{ g/cm}^3}$$

Correction check: Evaluating with standard molybdenum lattice parameter $a = 0.31470 \text{ nm} \implies V_c = 3.1167 \times 10^{-23} \text{ cm}^3$ gives $\rho_{\text{xrd}} = 10.224 \text{ g/cm}^3$. With $\rho_{\text{xrd}} = 10.224 \text{ g/cm}^3$ and $\rho_{\text{pyc}} = 10.220 \text{ g/cm}^3$, the experimental bulk density is slightly lower than the ideal lattice density.

• If excess interstitials were present, extra mass would be packed into the unit cell volume, making $\rho_{\text{pyc}} > \rho_{\text{xrd}}$.
• If vacancies are present, missing atoms reduce the mass of the crystal without reducing the X-ray measured unit cell volume by the same proportion, making $\rho_{\text{pyc}} < \rho_{\text{xrd}}$.

Because $\rho_{\text{pyc}} (10.220) < \rho_{\text{xrd}} (10.224)$, the predominant defects are vacancies.

The fraction of vacant lattice sites is given by:

$$X_v = \frac{\Delta\rho}{\rho_{\text{xrd}}} = \frac{\rho_{\text{xrd}} - \rho_{\text{pyc}}}{\rho_{\text{xrd}}} = \frac{10.224 - 10.220}{10.224} = \frac{0.004}{10.224} = 3.91 \times 10^{-4} \approx \mathbf{0.039\%}$$

The total atomic site density in BCC molybdenum is:

$$N = \frac{n}{V_c} = \frac{2}{3.117 \times 10^{-29} \text{ m}^3} = 6.417 \times 10^{28} \text{ sites/m}^3$$

The volumetric vacancy concentration is:

$$N_v = X_v \cdot N = (3.91 \times 10^{-4}) \times (6.417 \times 10^{28} \text{ m}^{-3}) = \mathbf{2.51 \times 10^{25} \text{ vacancies/m}^3}$$