K-Space

The K-Space sampling (i.e., the k-points used to sample the Brillouin Zone) is an important technical aspect of Band, as it influences heavily the accuracy, the CPU time and the memory usage of the calculation (see section Recommendations for k-space).

KSpace input block

The K-Space can be controlled via the KSpace input block. Two different k-space integration methods are available: the Regular Grid (default) and the Symmetric Grid.

KSpace
   Type [Regular | Symmetric]
   Quality [Auto | GammaOnly | Basic | Normal | Good | VeryGood | Excellent]
End
KSpace
Type

Block

Description

Options for the k-space integration (i.e. the grid used to sample the Brillouin zone)

Type
Type

Multiple Choice

Default value

Regular

Options

[Regular, Symmetric]

GUI name

K-space grid type

Description

The type of k-space integration grid used to sample the Brillouin zone (BZ) used. ‘Regular’: simple regular grid. ‘Symmetric’: symmetric grid for the irreducible wedge of the first BZ (useful when high-symmetry points in the BZ are needed to capture the correct physics of the system, graphene being a notable example).

Quality
Type

Multiple Choice

Default value

Auto

Options

[Auto, GammaOnly, Basic, Normal, Good, VeryGood, Excellent]

GUI name

K-space

Description

Select the quality of the K-space grid used to sample the Brillouin Zone. If ‘Auto’, the quality defined in the ‘NumericalQuality’ will be used. If ‘GammaOnly’, only one point (the gamma point) will be used. The actual number of K points generated depends on this option and on the size of the unit cell. The larger the real space cell, the fewer K points will be generated. The CPU-time and accuracy strongly depend on this option.

Regular K-Space grid

By default, Band will look at the size of a lattice vectors and the KSpace quality to determine the number of k-points. The larger the lattice vector in real space, the smaller the reciprocal space vectors are, and as a result fewer k-points are needed. The following intervals will be distinguished: 0-5 Bohr, 5-10 Bohr, 10-20 Bohr, 20-50 Bohr, and beyond. Here is the table explaining how many k-points will be used along a lattice vector.

Lattice vector length

Basic

Normal

Good

VeryGood

Excellent

0-5 Bohr

5

9

13

17

21

5-10 Bohr

3

5

9

13

17

10-20 Bohr

1

3

5

9

13

20-50 Bohr

1

1

3

5

9

50- Bohr…

1

1

1

3

5

By preferring odd-numbered values we can use a quadratic interpolation method, and have the \(\Gamma\) point in the grid. It is then reasonable to assume a decaying error when going to a better quality setting.

It is also possible to manually specify the number of k-space points along each reciprocal lattice vector

KSpace
   Regular
      NumberOfPoints integer_list
   End
End
KSpace
Type

Block

Description

Options for the k-space integration (i.e. the grid used to sample the Brillouin zone)

Regular
Type

Block

Description

Options for the regular k-space integration grid.

NumberOfPoints
Type

Integer List

Description

Use a regular grid with the specified number of k-points along each reciprocal lattice vector. For 1D periodic systems you should specify only one number, for 2D systems two numbers, and for 3D systems three numbers.

Symmetric K-Space grid (tetrahedron method)

The tetrahedron method can be useful when high symmetry points in the BZ are needed to capture the correct physics of the system, graphene being a notable example.

The number of k-points in the symmetric grid depends on the KSpace quality and on the length of the shortest lattice vector.

It is also possible to manually specify the symmetric k-space integration parameter:

KSpace
   Symmetric
      KInteg integer
   End
End
KSpace
Type

Block

Description

Options for the k-space integration (i.e. the grid used to sample the Brillouin zone)

Symmetric
Type

Block

Description

Options for the symmetric k-space integration grid.

KInteg
Type

Integer

GUI name

Accuracy

Description

Specify the accuracy for the Symmetric method. 1: absolutely minimal (only the G-point is used) 2: linear tetrahedron method, coarsest spacing 3: quadratic tetrahedron method, coarsest spacing 4,6,… (even): linear tetrahedron method 5,7…. (odd): quadratic method The tetrahedron method is usually by far inferior.

General Remark: The tetrahedron method samples the irreducible wedge of the first BZ, whereas the regular grid samples the whole, first BZ. As a rule of thumb you need to choose roughly twice the value for the regular grid. For example kspace 2 compares to grid 4 4 4, kspace 3 to grid 5 5 5, etc.. Sticking to this rule the number of unique k-points will be roughly similar.

Recommendations for k-space

Which K-Space quality to use depends very much a) the system you are studying and b) the property you are interested in. We strongly recommend you to test the effect of different K-Space qualities on your system and properties of interest.

As an example, in the following table we list the errors on formation energy and band gap for diamond using regular k-space grids of different qualities (using Excellent kSpace quality as reference).

Table 6 Accuracy of formation energy for diamond (primitive unit cell) using different KSpace grids

KSpace quality

Energy error / atom [eV]

CPU time ratio

Gamma-Only

3.3

1

Basic

0.6

2

Normal

0.03

6

Good

0.002

16

VeryGood

0.0001

35

Excellent

reference

64

It is worthwhile noting that the errors due to finite k-space sampling in formation energies are to some extend systematic, and they partially cancel each other out when taking energy differences.

In general, metals (or narrow-gap semiconductor) require higher K-Space sampling than insulators. For insulators and wide-gap semiconductors, Normal K-Space quality often suffices. For Narrow-gap semiconductor and metals, Good K-Space quality is highly recommended. For geometry optimizations under pressure, Good K-Space quality is recommended.

Furthermore for certain properties, such as band gaps, Normal K-Space quality might not be enough to obtain reliable results. For example, in following figure we see how Normal K-Space quality is often not enough for computing band gaps (especially for the narrow-gap semiconductor of the top panel). For band gap prediction, it is recommended to use Good K-Space quality.

/scm-uploads/doc/BAND/_images/band_gaps_k-space_jetsa.png

Fig. 4 Convergence of band gaps WRT the k-space quality set for various system. (XC:PBE, Basis:TZP)

High symmetry points and the regular grid

Using the symmetric k-grid it is ensured that points with high symmetry are included. However, for the default regular grid this is not the case. An important example is graphene. Its band structure has a conical intersection at the point labeled “K”, in all other k-points there is a gap.

Table 7 The magical values for which the regular k-grid includes the high symmetry point “K” (with fractional coordinates 1/3,1/3) for a honeycomb lattice.

grid

point “K” included?

K-grid quality for Graphene

5x5

No

Normal

7x7

Yes

9x9

No

Good

11x11

No

13x13

Yes

VeryGood

15x15

No

But in general it makes more sense to use the symmetric k-grid if points of high symmetry are deemed important.