How to choose efficiently the size of the Bethe-Salpeter equation Hamiltonian for accurate exciton calculations on supercells

This is the authors' abstract. We don't add key points for this paper.

The Bethe-Salpeter Equation (BSE) is the workhorse method to study excitons in materials. The BSE Hamiltonian size, which depends on how many valence-to-conduction band transitions are considered, needs to be chosen to be sufficiently large to converge excitons' energies and wavefunctions but should be minimized to make calculations tractable, as BSE calculations are expensive and scale with the number of atoms as O(Nₐₜₒₘₛ⁶). In particular, in the case of supercell (SC) calculations composed of Nᵣₑₚ replicas of a primitive cell (PC), a natural choice to build this BSE Hamiltonian is to include all transitions derived from PC calculations by zone folding. However, this leads to a very large BSE Hamiltonian, as the number of matrix elements in it is (Nₖ Nc Nᵥ)², where Nₖ is the number of k-points and Nc(v) is the number of conduction (valence) states used. When creating a SC, the number of k-points decreases by a factor Nᵣₑₚ but both the number of conduction and valence states increase by the same factor, therefore the number of matrix elements in the BSE Hamiltonian increases by a factor Nᵣₑₚ², making exactly corresponding calculations prohibitive. Here, we provide a workflow to decide how many transitions are necessary to achieve comparable results, based on only PC results. With our method, we show that to converge the first exciton binding energy of a LiF SC composed of 64 PCs, to an energy tolerance of 0.15 eV, we only need 12% of the valence-to-conduction matrix elements that result from zone folding with a minimal set of bands. As an example, we use the number of bands from our method to obtain the absorption spectrum of LiF with a Vₖ-like defect. The procedure in our work helps in evaluating excitonic properties in large SC calculations.