Rapid calculation method for Frenkel-type two-exciton states in one to three dimensions

Biexciton and two-exciton dissociated states of Frenkel-type excitons are well described by a tight-binding model with a nearest-neighbor approximation. Such two-exciton states in a finite-size lattice are usually calculated by numerical diagonalization of the Hamiltonian, which requires an increasing amount of computational time and memory as the lattice size increases. I develop here a rapid, memory-saving method to calculate the energies and wave functions of two-exciton states by employing a bisection method. In addition, an attractive interaction between two excitons in the tight-binding model can be obtained directly so that the biexciton energy agrees with the observed energy, without the need for the trial-and-error procedure implemented in the numerical diagonalization method.