Drawing Wigner-Seitz Cells

To draw a Wigner-Seitz cell, first the Bravais lattice vectors must be constructed from the primitive lattice vectors. The Bravais lattice vectors are, $\vec{R}_{hkl}=h\vec{a}_1+k\vec{a}_2+l\vec{a}_3$, where $\vec{a}_1,\,\vec{a}_2,\,\vec{a}_3$ are the primitive lattice vectors and $h,\,k$ and $l$ are integers.

The Bravais lattice vectors are,

$$\vec{R}_{hkl}=h\vec{a}_1+k\vec{a}_2+l\vec{a}_3=\left(ha_{1x}+ka_{2x}+la_{3x}\right)\hat{x}+\left(ha_{1y}+ka_{2y}+la_{3y}\right)\hat{y}+\left(ha_{1z}+ka_{2z}+la_{3z}\right)\hat{z}.$$

A plane normal to each Bravais lattice vector is drawn that passes through $\vec{R}_{hkl}/2$. All of the points that can be reached from the origin without crossing any of these planes is in the Wigner-Seitz cell.

A brief review of vectors and planes
The set of planes perpendicular to a vector $A_x\hat{x}+A_y\hat{y}+A_z\hat{z}$ is,

$$A_xx+A_yy+A_zz = C,$$

where $C$ is any constant. If a point $ (x_0,y_0,z_0)$ on the plane is known, $C$ can be calculated,

$$C=A_xx_0+A_yy_0+A_zz_0.$$

To find a point where three planes intersect, write down the equations that describe these three planes. These will be three linear equations with three unknowns. The app below will solve for the three unknowns and determine the point of intersection of the planes.

The form below takes the primitive lattice vectors in real space as input and calculates the Bravais lattice vectors $\vec{R}_{hkl} $, the planes $(hkl)$ that form the Wigner-Seitz cell boundaries, and the corners of the Wigner-Seitz cell.

A boundary of the Wigner-Seitz cell is a plane normal to $\vec{R}_{hkl}$, that passes through the point $\frac{\vec{R}_{hkl}}{2}$. For the planes that make up the Wigner-Seitz cell boundary, the distance from $\frac{\vec{R}_{hkl}}{2}$ to the origin is smaller than the distance from $\frac{\vec{R}_{hkl}}{2}$ to any of the other Bravais lattice vectors. By computing these distances, the planes that make up the Wigner-Seitz cell can be determined.

Once the planes are known, the points at the corners of the Wigner-Seitz cell can be determined by considering the intersections of the planes. The formula for the $(hkl)$ plane is,

$$R_{hkl,x}k_x+R_{hkl,y}k_y+R_{hkl,z}k_z=\frac{R_{hkl,x}^2}{2}+\frac{R_{hkl,y}^2}{2}+\frac{R_{hkl,z}^2}{2}.$$

By solving the sets of linear equations, the corners can be determined.

The edges of the Wigner-Seitz cell can be determined by considering all pairs of corners. If a pair of corners share two planes, there is an edge between these two corners.

Rotatable model of the Wigner-Seitz cell

//set up the three.js variables

var canvas1 = document.getElementById("canvas_bz");
var width = canvas1.clientWidth;
var height = canvas1.clientHeight;

var scene = new THREE.Scene();
scene.background = new THREE.Color(0xffffff);
var camera = new THREE.PerspectiveCamera(3, width/height, 1, 1000);

var renderer = new THREE.WebGLRenderer();
renderer.setSize(width, height);
canvas1.appendChild(renderer.domElement);

var controls = new THREE.OrbitControls(camera);

var material = new THREE.MeshBasicMaterial({
    color: 0x000000
});

var geometry;

//this while loop erases the previous BZ
  //while(scene.children.length > 0){ scene.remove(scene.children[0]); }

for (i = 0; i < Ne; i++) {
    geometry = new THREE.Geometry();
    geometry.vertices.push(new THREE.Vector3(edges[i][0], edges[i][1], edges[i][2]));
    geometry.vertices.push(new THREE.Vector3(edges[i][3], edges[i][4], edges[i][5]));

var line = new THREE.Line(geometry, material);

scene.add(line);

}

camera.position.z = 50;
  var animate = function() {
    
    requestAnimationFrame(animate);

//controls.update() must be called after any manual changes to the camera's transform
    controls.update();
    renderer.render(scene, camera);

}

animate();

</script>

Right-click and select 'view source' to see the code that generates this page. For some primitive lattice vectors, numerical errors occur and not all corners and edges are found.