A Complete Tutorial from Fundamentals to Computational Applications
By the end of this tutorial, you will be able to:
All Python scripts, Jupyter notebooks, and example structures are available on GitHub:
⬇️ View on GitHubIncludes: Jupyter notebooks • CIF files • Analysis scripts • Solutions
Before diving into mathematics, let's understand why symmetry is so fundamental to materials science. When you examine a salt crystal under a microscope, you see perfect cubic faces. Quartz crystals display hexagonal prisms. These shapes directly reflect the internal arrangement of atoms.
Cubic symmetry → isotropic properties
Trigonal symmetry → piezoelectric
Hexagonal layers → anisotropic
A cubic crystal like silicon has 48 symmetry operations. You only need to calculate in 1/48th of the Brillouin zone. For triclinic? The entire zone—48× more k-points!
Symmetry determines which tensor elements are zero, equal, or independent. Cubic dielectric tensor: 1 component. Triclinic: 6 components.
Every crystal is built from a repeating pattern. The mathematical framework describing this repetition is called a lattice. The smallest repeating unit is the unit cell.
Figure 2.1: A 2D lattice showing lattice points, unit cell, and lattice vectors.
A 3D unit cell is described by six parameters:
Three edge lengths: \(a\), \(b\), \(c\) (in Ångströms)
Three angles:
\(\alpha\) = angle between b and c | \(\beta\) = angle between a and c | \(\gamma\) = angle between a and b
Crystal Structure = Lattice + Basis
Both have FCC lattices, but different bases:
Same lattice, different basis → completely different materials!
Miller indices (hkl) describe the orientation of crystal planes. They represent reciprocals of the fractional intercepts.
Where does the plane intersect each axis? If parallel to an axis, intercept = ∞.
Calculate 1/(intercept). Note: 1/∞ = 0.
Multiply to get integers.
Negative values use a bar: \(\bar{1}\)
Figure 3.1: Common Miller indices showing planes in cubic unit cells.
For cubic crystals:
\[d_{hkl} = \frac{a}{\sqrt{h^2 + k^2 + l^2}}\]Silicon: a = 5.431 Å. Calculate d₁₁₁:
\[d_{111} = \frac{5.431}{\sqrt{1^2 + 1^2 + 1^2}} = \frac{5.431}{\sqrt{3}} = 3.136 \text{ Å}\]All crystals belong to one of seven crystal systems, defined by unit cell parameter constraints:
| Crystal System | Unit Cell Constraints | Minimum Symmetry | Examples |
|---|---|---|---|
| Cubic | a = b = c; α = β = γ = 90° | Four 3-fold axes | NaCl, Diamond, Cu |
| Tetragonal | a = b ≠ c; α = β = γ = 90° | One 4-fold axis | TiO₂, BaTiO₃ |
| Orthorhombic | a ≠ b ≠ c; α = β = γ = 90° | Three 2-fold axes | Sulfur, BaSO₄ |
| Hexagonal | a = b ≠ c; α = β = 90°, γ = 120° | One 6-fold axis | Graphite, ZnO |
| Trigonal | a = b = c; α = β = γ ≠ 90° | One 3-fold axis | Quartz, Calcite |
| Monoclinic | a ≠ b ≠ c; α = γ = 90° ≠ β | One 2-fold axis | Gypsum, organics |
| Triclinic | a ≠ b ≠ c; α ≠ β ≠ γ ≠ 90° | None | K₂Cr₂O₇ |
Barium titanate changes crystal system with temperature:
Within the seven crystal systems, there are 14 ways to arrange lattice points—the Bravais lattices.
Corners only
+ center
+ face centers
+ top/bottom
P: 8 × ⅛ = 1 | I: 1 + 1 = 2 | F: 1 + 6×½ = 4 | C: 1 + 2×½ = 2
Question: Copper has FCC structure, a = 3.615 Å, M = 63.55 g/mol. Calculate density.
FCC has 4 atoms per cell.
\[\rho = \frac{4 \times 63.55}{(3.615 \times 10^{-8})^3 \times 6.022 \times 10^{23}} = 8.94 \text{ g/cm}^3\]Experimental: 8.96 g/cm³ ✓
Symmetry operations are transformations that leave a crystal indistinguishable from its original state.
Do nothing. Every object has this symmetry.
Rotate by 360°/n. Only n = 1, 2, 3, 4, 6 allowed (crystallographic restriction).
Mirror reflection across a plane.
Point (x,y,z) maps to (−x,−y,−z).
Figure 6.1: The fundamental symmetry operations in crystallography.
where \(\mathbf{r} = (x, y, z)^T\) and \(\mathbf{W}\) is the operation matrix
A point group is a collection of symmetry operations that leave at least one point unmoved. In crystallography, only 32 point groups are compatible with 3D periodicity.
| Symbol | Meaning | Example |
|---|---|---|
| 1, 2, 3, 4, 6 | n-fold rotation axis | 4 = 90° rotation |
| m | Mirror plane | reflection |
| 1̄, 3̄, 4̄, 6̄ | Rotoinversion axis | 4̄ = 90° + inversion |
| / | Perpendicular to | 4/m = 4-fold ⊥ mirror |
| Crystal System | Point Groups | Count |
|---|---|---|
| Triclinic | 1, 1̄ | 2 |
| Monoclinic | 2, m, 2/m | 3 |
| Orthorhombic | 222, mm2, mmm | 3 |
| Tetragonal | 4, 4̄, 4/m, 422, 4mm, 4̄2m, 4/mmm | 7 |
| Trigonal | 3, 3̄, 32, 3m, 3̄m | 5 |
| Hexagonal | 6, 6̄, 6/m, 622, 6mm, 6̄m2, 6/mmm | 7 |
| Cubic | 23, m3̄, 432, 4̄3m, m3̄m | 5 |
| Total | 32 | |
Space groups combine point group symmetry with translational symmetry, including new elements:
Rotation by 360°/n + translation by m/n along the axis.
Example: 2₁ = 180° rotation + ½ translation
Mirror reflection + translation parallel to the mirror.
Types: a, b, c (axial), n (diagonal), d (diamond)
~35% of organic crystal structures have this space group!
Space group of Cu, Al, Au, Ag, Ni, and many FCC metals.
Wyckoff positions describe allowed locations for atoms within a unit cell, classified by site symmetry.
Each Wyckoff position has:
| Mult. | Letter | Site Sym. | Coordinates |
|---|---|---|---|
| 1 | a | m3̄m | (0,0,0) |
| 1 | b | m3̄m | (½,½,½) |
| 3 | c | 4/mmm | (0,½,½), etc. |
| 48 | l | 1 | (x,y,z) general |
The reciprocal lattice is the Fourier transform of real space, essential for understanding diffraction and band structures.
Similar expressions for \(\mathbf{b}^*\) and \(\mathbf{c}^*\)
Bragg's Law: \(n\lambda = 2d_{hkl}\sin\theta\)
Reciprocal space: Diffraction when \(\Delta\mathbf{k} = \mathbf{G}_{hkl}\)
Crystal symmetry causes certain diffraction reflections to have zero intensity. These are the key to determining space groups.
Intensity: \(I_{hkl} \propto |F_{hkl}|^2\)
Atoms at (0,0,0) and (½,½,½):
\[F_{hkl} = f[1 + e^{i\pi(h+k+l)}]\]| Symmetry Element | Affected Reflections | Condition for Absence |
|---|---|---|
| Body-centered (I) | All hkl | h + k + l = odd |
| Face-centered (F) | All hkl | Mixed parity |
| C-centered | All hkl | h + k = odd |
| 2₁ screw ∥ c | 00l only | l = odd |
| c-glide ⊥ b | h0l only | l = odd |
Crystal symmetry directly determines which physical properties are allowed.
| Property | Tensor Rank | Cubic | Triclinic |
|---|---|---|---|
| Dielectric constant | 2 | 1 component | 6 components |
| Thermal expansion | 2 | 1 component | 6 components |
| Elastic stiffness | 4 | 3 components | 21 components |
| Piezoelectric | 3 | 0 (if centrosym.) | 18 components |
Requirement: Non-centrosymmetric
Allowed: 20 of 32 point groups
Examples: Quartz, ZnO, BaTiO₃
Requirement: Polar point group
Allowed: 10 of 32 point groups
Examples: BaTiO₃, LiNbO₃, PZT
Question: Can pure silicon (space group Fd3̄m, point group m3̄m) be piezoelectric?
No. m3̄m contains inversion (3̄), making it centrosymmetric. Piezoelectricity requires non-centrosymmetric crystals.
All code examples, Jupyter notebooks, and CIF files are available in the GitHub repository:
Clone or download all examples:
github.com/nabil-khossossi/crystal-symmetry-tutorial01_spglib_basics.ipynb - Symmetry detection02_pymatgen_analysis.ipynb - Structure analysis03_wyckoff_positions.ipynb - Site symmetry04_property_prediction.ipynb - Tensorsstructures/ - CIF files for practicesolutions/ - Exercise answersscripts/ - Standalone Python scriptsdata/ - Reference tables# Install required packages pip install spglib pymatgen ase # Clone the repository git clone https://github.com/nabil-khossossi/crystal-symmetry-tutorial.git cd crystal-symmetry-tutorial # Launch Jupyter jupyter notebook
# spglib - Symmetry detection
import spglib
symmetry = spglib.get_symmetry_dataset(cell)
print(f"Space group: {symmetry['international']}")
# pymatgen - Structure analysis
from pymatgen.symmetry.analyzer import SpacegroupAnalyzer
analyzer = SpacegroupAnalyzer(structure)
print(f"Point group: {analyzer.get_point_group_symbol()}")
# ASE - Structure building
from ase.build import bulk
cu = bulk('Cu', 'fcc', a=3.615)
01_spglib_basics.ipynb for the fundamentals.
A crystal has: a = b = c = 4.0 Å, α = β = γ = 90°. All hkl with h+k+l = odd are absent.
Questions: (a) Crystal system? (b) Bravais lattice?
(a) Cubic (a=b=c, all 90°)
(b) Body-centered (I) (h+k+l odd absent)
Monoclinic crystal. Absences: (h0l) l=odd, (0k0) k=odd. What space group?
P2₁/c (No. 14)
Space group P6₃mc (point group 6mm). Can it be: piezoelectric? ferroelectric?
Both YES
Example: ZnO (wurtzite)
You now have the tools to understand and analyze crystal structures.
Symmetry is the key to understanding materials!
Dr. Nabil Khossossi | DIFFER | sustai-nabil.com