Exact diagonalization of Landau levels in Bernal graphite

The magneto-electronic structure of Bernal graphite is derived by the Peierls tight-binding model. Seven important atomic interactions are taken into account and the corresponding Hamiltonian matrix is rearranged as a band-like form for efficient calculation. The results derived from this exact diagonalization method demonstrate that interlayer interactions play remarkable roles in Landau level structures.