We introduce a new computational framework that incorporates multiple scattering for large-scale 3D particle localization using single-shot in-line holography. Traditional holographic techniques rely on single-scattering models that become inaccurate under high particle densities and large refractive index contrasts. Existing multiple scattering solvers become computationally prohibitive for large-scale problems, which comprise of millions of voxels within the scattering volume. Our approach overcomes the computational bottleneck by slice-wise computation of multiple scattering under an efficient recursive framework. In the forward model, each recursion estimates the next higher-order multiple scattered field among the object slices. In th...