We develop an improved stochastic formalism for the Bethe-Salpeter equation,
based on an exact separation of the effective-interaction W to two parts,
W=(WβvWβ)+vWβ where the latter is formally any translationally-invariant
interaction vWβ(rβrβ²). When optimizing the fit of vWβ exchange kernel to
W, by using a stochastic sampling of W, the difference WβvWβ becomes
quite small. Then, in the main BSE routine, this small difference is
stochastically sampled. The number of stochastic samples needed for an accurate
spectrum is then largely independent of system size. While the method is
formally cubic in scaling, the scaling prefactor is small due to the constant
number of stochastic orbitals needed for sampling W.Comment: 9 pages, 5 figures, 2 table