Band gap

ase.dft.bandgap.bandgap(calc=None, direct=False, spin=None, output='-', eigenvalues=None, efermi=None, kpts=None)[source]

Calculates the band-gap.

Parameters:

calc: Calculator object
Electronic structure calculator object.
direct: bool
Calculate direct band-gap.
spin: int or None
For spin-polarized systems, you can use spin=0 or spin=1 to look only at a single spin-channel.
output: file descriptor
Use output=None for no text output or ‘-‘ for stdout (default).
eigenvalues: ndarray of shape (nspin, nkpt, nband) or (nkpt, nband)
Eigenvalues.
efermi: float
Fermi level (defaults to 0.0).
kpts: ndarray of shape (nkpt, 3)
For pretty text output only.

Returns a (gap, p1, p2) tuple where p1 and p2 are tuples of indices of the valence and conduction points (s, k, n).

Example:

>>> gap, p1, p2 = bandgap(silicon.calc)
Gap: 1.2 eV
Transition (v -> c):
    [0.000, 0.000, 0.000] -> [0.500, 0.500, 0.000]
>>> print(gap, p1, p2)
1.2 (0, 0, 3), (0, 5, 4)
>>> gap, p1, p2 = bandgap(silicon.calc, direct=True)
Direct gap: 3.4 eV
Transition at: [0.000, 0.000, 0.000]
>>> print(gap, p1, p2)
3.4 (0, 0, 3), (0, 0, 4)