Crystal Structures & Symmetry - Complete Tutorial | Dr. Nabil Khossossi

💎 Crystal Structures & Symmetry

A Complete Tutorial from Fundamentals to Computational Applications

Tutorial 3 | Fundamentals Series | Dr. Nabil Khossossi
Estimated time: 3-4 hours | Level: Beginner to Intermediate

📚 Prerequisites

  • Basic Python knowledge (Tutorial 2)
  • Familiarity with ASE (Tutorial 1)
  • Basic linear algebra (vectors, matrices)
  • High school chemistry (atoms, bonds)

🎯 Learning Objectives

By the end of this tutorial, you will be able to:

  • Understand the mathematical language of crystal symmetry
  • Identify and classify crystal structures using the 7 crystal systems and 14 Bravais lattices
  • Read and interpret space group symbols (Hermann-Mauguin notation)
  • Determine Wyckoff positions and site symmetries
  • Predict material properties from crystal symmetry
  • Use Python tools (spglib, pymatgen, ASE) for symmetry analysis

📂 Tutorial Code Repository

All Python scripts, Jupyter notebooks, and example structures are available on GitHub:

⬇️ View on GitHub

Includes: Jupyter notebooks • CIF files • Analysis scripts • Solutions

1. Why Does Symmetry Matter in Crystals?

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.

Every physical property of a crystal—electrical conductivity, thermal expansion, piezoelectricity, optical behavior—is fundamentally determined by its atomic arrangement and symmetry. Understanding symmetry is understanding materials.

Real-World Examples

Table Salt (NaCl)

Cubic symmetry → isotropic properties

Quartz (SiO₂)

Trigonal symmetry → piezoelectric

Graphite

Hexagonal layers → anisotropic

Why Symmetry Matters for DFT Calculations

Computational Efficiency

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!

Property Tensors

Symmetry determines which tensor elements are zero, equal, or independent. Cubic dielectric tensor: 1 component. Triclinic: 6 components.

2. Building Blocks: Lattices and Unit Cells

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.

a b 2D Lattice • Blue dots = lattice points • Red box = unit cell • Red arrows = lattice vectors Translation vector: T = n₁a + n₂b + n₃c

Figure 2.1: A 2D lattice showing lattice points, unit cell, and lattice vectors.

Unit Cell Parameters

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

Lattice Points vs. Atoms

Lattice points are mathematical abstractions showing where the pattern repeats.
Atoms are the physical matter. A single lattice point can have multiple atoms associated with it (the basis).

Crystal Structure = Lattice + Basis

Example: NaCl vs. Diamond

Both have FCC lattices, but different bases:

  • NaCl: Basis = Na at (0,0,0) + Cl at (½,½,½)
  • Diamond: Basis = C at (0,0,0) + C at (¼,¼,¼)

Same lattice, different basis → completely different materials!

3. Miller Indices: The Language of Crystal Planes

Miller indices (hkl) describe the orientation of crystal planes. They represent reciprocals of the fractional intercepts.

How to Determine Miller Indices

Step 1: Find the intercepts

Where does the plane intersect each axis? If parallel to an axis, intercept = ∞.

Step 2: Take reciprocals

Calculate 1/(intercept). Note: 1/∞ = 0.

Step 3: Clear fractions

Multiply to get integers.

Step 4: Write as (hkl)

Negative values use a bar: \(\bar{1}\)

(100) Intercepts: 1, ∞, ∞ (110) Intercepts: 1, 1, ∞ (111) Intercepts: 1, 1, 1 (210) Intercepts: ½, 1, ∞

Figure 3.1: Common Miller indices showing planes in cubic unit cells.

Interplanar Spacing (d-spacing)

For cubic crystals:

\[d_{hkl} = \frac{a}{\sqrt{h^2 + k^2 + l^2}}\]

Example: Silicon d-spacing

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{ Å}\]

4. The Seven Crystal Systems

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₇

Physical Example: BaTiO₃ Phase Transitions

Barium titanate changes crystal system with temperature:

  • T > 120°C: Cubic (paraelectric)
  • 5°C < T < 120°C: Tetragonal (ferroelectric)
  • -90°C < T < 5°C: Orthorhombic (ferroelectric)
  • T < -90°C: Trigonal (ferroelectric)

5. Bravais Lattices: The 14 Fundamental Patterns

Within the seven crystal systems, there are 14 ways to arrange lattice points—the Bravais lattices.

Types of Centering

P (Primitive)

Corners only

I (Body-centered)

+ center

F (Face-centered)

+ face centers

C (Base-centered)

+ top/bottom

Lattice Points per Unit Cell

P: 8 × ⅛ = 1  |  I: 1 + 1 = 2  |  F: 1 + 6×½ = 4  |  C: 1 + 2×½ = 2

Density Calculation

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³ ✓

6. Symmetry Operations: The Mathematical Tools

Symmetry operations are transformations that leave a crystal indistinguishable from its original state.

The Four Basic Operations

1. Identity (E or 1)

Do nothing. Every object has this symmetry.

2. Rotation (Cn or n)

Rotate by 360°/n. Only n = 1, 2, 3, 4, 6 allowed (crystallographic restriction).

3. Reflection (σ or m)

Mirror reflection across a plane.

4. Inversion (i or 1̄)

Point (x,y,z) maps to (−x,−y,−z).

Rotation (C₄) 90° rotation Reflection (m) Mirror plane Inversion (1̄) (x,y,z) → (−x,−y,−z) Rotoinversion (4̄) Rotate + invert

Figure 6.1: The fundamental symmetry operations in crystallography.

Matrix Representation

\[\mathbf{r'} = \mathbf{W} \cdot \mathbf{r}\]

where \(\mathbf{r} = (x, y, z)^T\) and \(\mathbf{W}\) is the operation matrix

90° Rotation about z-axis

\[\mathbf{C}_4^z = \begin{pmatrix} 0 & -1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 1 \end{pmatrix}\]

Inversion

\[\mathbf{i} = \begin{pmatrix} -1 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & -1 \end{pmatrix}\]
Symmetry operations form mathematical groups: closure, associativity, identity, and inverses. This is why group theory is the natural language of crystallography.

7. Point Groups: Local Symmetry (32 Types)

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.

Hermann-Mauguin Notation

SymbolMeaningExample
1, 2, 3, 4, 6n-fold rotation axis4 = 90° rotation
mMirror planereflection
1̄, 3̄, 4̄, 6̄Rotoinversion axis4̄ = 90° + inversion
/Perpendicular to4/m = 4-fold ⊥ mirror

The 32 Crystallographic Point Groups

Crystal SystemPoint GroupsCount
Triclinic1, 1̄2
Monoclinic2, m, 2/m3
Orthorhombic222, mm2, mmm3
Tetragonal4, 4̄, 4/m, 422, 4mm, 4̄2m, 4/mmm7
Trigonal3, 3̄, 32, 3m, 3̄m5
Hexagonal6, 6̄, 6/m, 622, 6mm, 6̄m2, 6/mmm7
Cubic23, m3̄, 432, 4̄3m, m3̄m5
Total32

Centrosymmetric vs. Non-centrosymmetric

A point group is centrosymmetric if it contains inversion (1̄). Only 11 of the 32 point groups are centrosymmetric.

Physical significance: Non-centrosymmetric crystals can exhibit piezoelectricity (20 groups) and pyroelectricity (10 polar groups).

8. Space Groups: Complete 3D Symmetry (230 Types)

Space groups combine point group symmetry with translational symmetry, including new elements:

Screw Axes (nm)

Rotation by 360°/n + translation by m/n along the axis.

Example: 2₁ = 180° rotation + ½ translation

Glide Planes

Mirror reflection + translation parallel to the mirror.

Types: a, b, c (axial), n (diagonal), d (diamond)

Space Group Notation

Example: P2₁/c (most common for organics!)

  • P: Primitive lattice
  • 2₁: 2₁ screw axis along b
  • /c: c-glide plane perpendicular to b

~35% of organic crystal structures have this space group!

Example: Fm3̄m (FCC metals)

  • F: Face-centered lattice
  • m3̄m: Full cubic point group symmetry

Space group of Cu, Al, Au, Ag, Ni, and many FCC metals.

Although there are 230 space groups, their occurrence is highly uneven:
  • P2₁/c: ~35% of organic crystals
  • P1̄: ~15% of organic crystals
  • P2₁2₁2₁: ~10% (chiral molecules)

9. Wyckoff Positions: Where Atoms Can Sit

Wyckoff positions describe allowed locations for atoms within a unit cell, classified by site symmetry.

Each Wyckoff position has:

  • Multiplicity: Number of equivalent positions
  • Wyckoff letter: Label (a, b, c, ...)
  • Site symmetry: Point symmetry at that position
  • Coordinates: Fractional positions

Example: Space Group Pm3̄m (No. 221)

Mult.LetterSite Sym.Coordinates
1am3̄m(0,0,0)
1bm3̄m(½,½,½)
3c4/mmm(0,½,½), etc.
48l1(x,y,z) general
For Structure Setup: You only need to specify ONE atom per Wyckoff position. Symmetry automatically generates the rest, reducing parameters in DFT calculations.

10. Reciprocal Space and Diffraction Basics

The reciprocal lattice is the Fourier transform of real space, essential for understanding diffraction and band structures.

\[\mathbf{a}^* = 2\pi \frac{\mathbf{b} \times \mathbf{c}}{\mathbf{a} \cdot (\mathbf{b} \times \mathbf{c})}\]

Similar expressions for \(\mathbf{b}^*\) and \(\mathbf{c}^*\)

Key Properties

  • \(\mathbf{a} \cdot \mathbf{a}^* = 2\pi\), but \(\mathbf{a} \cdot \mathbf{b}^* = 0\)
  • Reciprocal of FCC is BCC, and vice versa
  • The first Brillouin zone contains all unique k-points

Connection to Diffraction

Bragg's Law: \(n\lambda = 2d_{hkl}\sin\theta\)

Reciprocal space: Diffraction when \(\Delta\mathbf{k} = \mathbf{G}_{hkl}\)

11. Systematic Absences: Symmetry in Diffraction

Crystal symmetry causes certain diffraction reflections to have zero intensity. These are the key to determining space groups.

Structure Factor

\[F_{hkl} = \sum_j f_j \exp[2\pi i(hx_j + ky_j + lz_j)]\]

Intensity: \(I_{hkl} \propto |F_{hkl}|^2\)

Example: Body-Centered Cubic (BCC)

Atoms at (0,0,0) and (½,½,½):

\[F_{hkl} = f[1 + e^{i\pi(h+k+l)}]\]
  • h+k+l = even: F = 2f ✓ Allowed
  • h+k+l = odd: F = 0 ✗ Absent

Common Absence Rules

Symmetry ElementAffected ReflectionsCondition for Absence
Body-centered (I)All hklh + k + l = odd
Face-centered (F)All hklMixed parity
C-centeredAll hklh + k = odd
2₁ screw ∥ c00l onlyl = odd
c-glide ⊥ bh0l onlyl = odd

12. Structure-Property Relationships

Crystal symmetry directly determines which physical properties are allowed.

Neumann's Principle: The symmetry of any physical property must include all symmetry elements of the crystal's point group. Properties cannot have lower symmetry than the crystal.

Tensor Properties and Symmetry

PropertyTensor RankCubicTriclinic
Dielectric constant21 component6 components
Thermal expansion21 component6 components
Elastic stiffness43 components21 components
Piezoelectric30 (if centrosym.)18 components

Piezoelectricity

Requirement: Non-centrosymmetric

Allowed: 20 of 32 point groups

Examples: Quartz, ZnO, BaTiO₃

Ferroelectricity

Requirement: Polar point group

Allowed: 10 of 32 point groups

Examples: BaTiO₃, LiNbO₃, PZT

Property Prediction

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.

13. Computational Tools (Hands-On)

All code examples, Jupyter notebooks, and CIF files are available in the GitHub repository:

📂 Complete Code Repository

Clone or download all examples:

github.com/nabil-khossossi/crystal-symmetry-tutorial

Repository Contents

📓 Jupyter Notebooks

  • 01_spglib_basics.ipynb - Symmetry detection
  • 02_pymatgen_analysis.ipynb - Structure analysis
  • 03_wyckoff_positions.ipynb - Site symmetry
  • 04_property_prediction.ipynb - Tensors

📁 Example Files

  • structures/ - CIF files for practice
  • solutions/ - Exercise answers
  • scripts/ - Standalone Python scripts
  • data/ - Reference tables

Quick Start

# 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

Key Functions Overview

# 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)
The notebooks contain step-by-step tutorials with explanations, visualizations, and exercises. Start with 01_spglib_basics.ipynb for the fundamentals.

14. Practice Problems and Summary

Summary: What We've Learned

Core Concepts Mastered

  • ✓ Unit cells and lattice parameters
  • ✓ Miller indices for planes/directions
  • ✓ 7 crystal systems
  • ✓ 14 Bravais lattices
  • ✓ Symmetry operations
  • ✓ 32 point groups
  • ✓ 230 space groups
  • ✓ Hermann-Mauguin notation
  • ✓ Wyckoff positions
  • ✓ Reciprocal space basics
  • ✓ Systematic absences
  • ✓ Structure-property relationships
  • ✓ Python tools (spglib, pymatgen)
  • ✓ Analysis workflow

Practice Problems

Problem 1: Crystal System

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)

Problem 2: Space Group

Monoclinic crystal. Absences: (h0l) l=odd, (0k0) k=odd. What space group?

P2₁/c (No. 14)

  • Primitive (no general absences)
  • (h0l) l=odd → c-glide ⊥ b
  • (0k0) k=odd → 2₁ screw ∥ b

Problem 3: Properties

Space group P6₃mc (point group 6mm). Can it be: piezoelectric? ferroelectric?

Both YES

  • 6mm is non-centrosymmetric → piezoelectric ✓
  • 6mm is polar (unique c-axis) → ferroelectric possible ✓

Example: ZnO (wurtzite)

Essential Resources

Books

  • International Tables for Crystallography, Vol. A
  • Space Groups for Solid State Scientists — Burns & Glazer
With this foundation, you're ready for:
  • Tutorial 4: Advanced DFT calculations with symmetry
  • Tutorial 5: Electronic band structures and k-paths
  • Tutorial 6: Phonons and thermal properties