Huebsch: Created page with "A crystal is characterized by the fact that it obeys translational symmetry. In many calculations, we only consider the primitive crystal unit cell to save computational time. However in a realistic bulk material, interactions go beyond the primitive unit cell. The concept of '''crystal momentum''' is crucial in order to take into account interactions that go beyond the primitive unit cell during the simulation and it is a consequence of translational invariance. Fo..."

2024-10-21T11:32:21Z

Created page with "A crystal is characterized by the fact that it obeys translational symmetry. In many calculations, we only consider the primitive crystal unit cell to save computational time. However in a realistic bulk material, interactions go beyond the primitive unit cell. The concept of '''crystal momentum''' is crucial in order to take into account interactions that go beyond the primitive unit cell during the simulation and it is a consequence of translational invariance. Fo..."

New page

A crystal is characterized by the fact that it obeys translational [[symmetry]]. In many calculations, we only consider the primitive crystal unit cell to save computational time. However in a realistic bulk material, interactions go beyond the primitive unit cell. The concept of '''crystal momentum''' is crucial in order to take into account interactions that go beyond the primitive unit cell during the simulation and it is a consequence of translational invariance.

Formally, translational invariance can be written as a translation operator that commutes with the Hamiltonian:

<math display="block">[T_{\mathbf R},H]=0 </math>

For the [[Projector-augmented-wave formalism|KS orbitals]] this implies that each translation can only add a phase factor:

<math display="block">T_\mathbf{R} \psi_{n\mathbf{k}}(\mathbf{r}) = \text{e}^{\text{i}\mathbf{k}\cdot\mathbf{R}}\psi_{n\mathbf{k}}(\mathbf{r}),</math>

where <math>n</math> is the index that goes over all [[Projector-augmented-wave formalism|KS orbitals]].
Performing two consecutive translations should yield the sum of the individual translations:

<math display="block">T_{\mathbf{R}_1}T_{\mathbf{R}_2} \psi_{n\mathbf k}(\mathbf{r})=T_{\mathbf{R}_1+\mathbf{R}_2}\psi_{n\mathbf k}(\mathbf{r}).</math>

For a system with translational invariance, we obtain a periodic potential and, hence, it is most convenient to use periodic boundary conditions. The figure illustrates a periodic potential and the blue box is a possible unit cell.

This means, by virtue of the ''Bloch theorem'', we can separate each KS orbital into

<math display="block"> \psi_{n\mathbf{k}}(\mathbf{r})=u_{n\mathbf{k}}(\mathbf{r})\text{e}^{\text{i}\mathbf{k}\cdot\mathbf{r}} </math>

a cell periodic part <math>u</math> and a phase factor <math>\phi = \mathbf{k}\cdot\mathbf{r}</math>. The phase <math>\phi</math> is called Bloch factor and <math>\mathbf{k}</math> is the '''crystal momentum''' which lives in reciprocal space. Mind that <math>u_{n\mathbf{k}}</math> depends on <math>\mathbf{k}</math> and recall that n is the index that goes over all KS orbitals.

The reciprocal space is spanned by reciprocal lattice vectors <math>\mathbf{b}_i</math>:

<math display="block">\mathbf{b}_1 = \frac{2\pi}{\Omega} \mathbf{a}_2 \times \mathbf{a}_3</math>
<math display="block">\mathbf{b}_2 = \frac{2\pi}{\Omega} \mathbf{a}_3 \times \mathbf{a}_1</math>
<math display="block">\mathbf{b}_3 = \frac{2\pi}{\Omega} \mathbf{a}_1 \times \mathbf{a}_2</math>

These are defined in terms of the real space lattice vectors <math>a_i</math> ({{FILE|POSCAR}}). <math>\Omega = \mathbf{a}_1 \cdot \mathbf{a}_2 \times \mathbf{a}_3</math> is the volume of the unit cell. Note that a short real space direction actually yields a long direction in reciprocal space:

<math display="block">\mathbf{a}_i\cdot \mathbf{b}_j = 2\pi \delta_{ij}</math>

Based on the reciprocal lattice vectors <math>\mathbf{b}_i</math>, we can identify a ''primitive cell in reciprocal space''; This is the so-called first Brillouin zone (1. BZ). For an [[Electronic minimization|SCF calculation]], in practice the <math>\mathbf{k}</math> points are defined on a regular mesh with evenly spaced <math>\mathbf{k}</math> points in the 1. BZ. These are specified in the {{FILE|KPOINTS}} file or by using {{TAG|KSPACING}}.

To get a deeper understanding, consider an integration over all of real space (as is frequently required to compute properties): The connection between real space and reciprocal space is a Fourier transformation. That is, the integral over all of space can be expressed in terms of an integral over the 1. BZ:

<math display="block"> \int_{-\infty}^{\infty} \text{d}^3 r\, f(\mathbf{r}) =\frac{1}{\Omega_{BZ}}\int_{1.\,BZ} \text{d}^3 k\, \tilde{f}(\mathbf{k}) </math>

In principle one has to include an infinite number of <math>\mathbf{k}</math> points in the 1. BZ to describe interactions with periodic replica of the unit cell. In practice, the <math>\mathbf{k}</math> points are defined on a regular <math>\mathbf{k}</math> mesh and beyond a certain <math>\mathbf{k}</math>-points density the result (i.e. some quantity of interest) will not change (i.e. it converges w.r.t. the number of <math>\mathbf{k}</math> points). This is because the crystal momentum vectors that are close together are almost identical. Hence, we obtain

<math display="block"> \int_{-\infty}^{\infty} \text{d}^3 r\, f(\mathbf{r}) =\frac{1}{\Omega_{BZ}}\int_{1.\,BZ} \text{d}^3 k\, \tilde{f}(\mathbf{k}) = \sum_{\mathbf k \in 1.\,BZ} \tilde{f}(\mathbf{k}) \Delta^3k</math>

With a closer look at the first Brillouin zone, we notice that it also has a certain symmetry. As a consequence, some <math>\mathbf{k}</math> points on the regular mesh are equivalent. So, we only have to compute the KS orbitals on an irreducible set of <math>\mathbf{k}</math> points. VASP automatically finds these irreducible <math>\mathbf{k}</math> points and applies an appropriate weight <math> w_\mathbf{k} </math> in the sum over k points.

<math display="block"> \int_{-\infty}^{\infty} \text{d}^3 r\, f(\mathbf{r}) =\frac{1}{\Omega_{BZ}}\int_{1.\,BZ} \text{d}^3 k\, \tilde{f}(\mathbf{k}) = \sum_{\mathbf k \in 1.\,BZ} \tilde{f}(\mathbf{k}) \Delta^3k= \sum_{\mathbf{k} \in irred\, 1.\,BZ} w_\mathbf{k} \tilde{f}(\mathbf{k}) </math>

To see which [[symmetry]] VASP identified and which irreducible <math>\mathbf{k}</math> points are used during the calculation, look at the {{FILE|OUTCAR}} and {{FILE|IBZKPT}} files.

[[Category:Symmetry]]

Category:Crystal momentum - Revision history