Skip to content
Snippets Groups Projects

Release version

Merged Kostas Vilkelis requested to merge release_version into main
All threads resolved!
Files
8
+ 0
25
@@ -57,31 +57,6 @@ def generate_vectors(cutoff, dim):
return [*product(*([[*range(-cutoff, cutoff + 1)]] * dim))]
def compute_gap(tb, fermi_energy=0, nk=100):
"""Compute gap.
Parameters
----------
tb : dict
Tight-binding model for which to compute the gap.
fermi_energy : float
Fermi energy.
nk : int
Number of k-points to sample along each dimension.
Returns
-------
gap : float
Indirect gap.
"""
kham = tb_to_khamvector(tb, nk, ks=None)
vals = np.linalg.eigvalsh(kham)
emax = np.max(vals[vals <= fermi_energy])
emin = np.min(vals[vals > fermi_energy])
return np.abs(emin - emax)
def calculate_fermi_energy(tb, filling, nk=100):
"""Calculate the Fermi energy for a given filling."""
kham = tb_to_khamvector(tb, nk, ks=None)
Loading