DOS

sisl.physics.electron.DOS(E, eig, distribution='gaussian')[source]

Calculate the density of states (DOS) for a set of energies, E, with a distribution function

The DOS(E) is calculated as:

DOS(E)=iD(Eϵi)δ(Eϵi)

where D(ΔE) is the distribution function used. Note that the distribution function used may be a user-defined function. Alternatively a distribution function may be retrieved from Distribution functions.

Parameters:
  • E (ndarray) – energies to calculate the DOS at

  • eig (ndarray) – electronic eigenvalues

  • distribution (Literal['gaussian', 'lorentzian', 'fermi', 'bose-einstein', 'cold', 'step-function', 'heaviside'] | ~typing.Callable[[~collections.abc.Buffer | ~numpy._typing._array_like._SupportsArray[~numpy.dtype[~typing.Any]] | ~numpy._typing._nested_sequence._NestedSequence[~numpy._typing._array_like._SupportsArray[~numpy.dtype[~typing.Any]]] | bool | int | float | complex | str | bytes | ~numpy._typing._nested_sequence._NestedSequence[bool | int | float | complex | str | bytes]], ~numpy.ndarray]) – a function that accepts ΔE as argument and calculates the distribution function.

See also

Distribution functions

a selected set of implemented distribution functions

COP

calculate COOP or COHP curves

PDOS

projected DOS (same as this, but projected onto each orbital)

spin_moment

spin moment

Returns:

DOS calculated at energies, has same length as E

Return type:

numpy.ndarray

Parameters:

distribution (Literal['gaussian', 'lorentzian', 'fermi', 'bose-einstein', 'cold', 'step-function', 'heaviside'] | ~typing.Callable[[~collections.abc.Buffer | ~numpy._typing._array_like._SupportsArray[~numpy.dtype[~typing.Any]] | ~numpy._typing._nested_sequence._NestedSequence[~numpy._typing._array_like._SupportsArray[~numpy.dtype[~typing.Any]]] | bool | int | float | complex | str | bytes | ~numpy._typing._nested_sequence._NestedSequence[bool | int | float | complex | str | bytes]], ~numpy.ndarray])