“…The coupled differential Heisenberg equations describing the nonlinear optical response are truncated to third order in the exciting electro magnetic field O (E 3 ) [65][66][67][68][69][70]. In this limit, the relevant quantities are the optically addressable interband transitions c † ξ,s,k v ξ,s,k , building up a macroscopic polarization, which also couples to four-point correlation functions beyond a Hartree-Fock description [71,72], see figure 2: In order to find a numerically feasible approach, the interband transitions are expanded in terms of excitonic wavefunctions (section 2.1) and the four-particle correlations, equation (1), in a set of corresponding biexcitonic wavefunctions, including both bound and continuum states with respect to the energetically lowest 1s excitonic resonances (section 2.2). In particular, besides bound biexcitons, the biexcitonic continuum of exciton-exciton scattering states is required to get a quantitative understanding.…”