if dk is None:
            dk = np.linspace(self.dk_bohr[0], self.dk_bohr[-1], 100)