Return the jellium Hamiltonian as QubitOperator in the dual basis.
spinless: bool = False,
include_constant: bool = False
grid (Grid): The discretization to use.
spinless (bool): Whether to use the spinless model or not.
include_constant (bool): Whether to include the Madelung constant.
Note constant is unsupported for non-uniform, non-cubic cells with