Equations (REF ) and (REF ) define the multiple scattering problem. The exciting field \({\bf E}_\text{exc}^p\)  depends on the positions and properties of all other particles and is generally unknown. This set of equations can be solved numerically on finite ensembles of particles using T-matrix methods [1]}, as mentioned above, however, the scattering response of a specific configuration of particles is generally of limited interest. Instead, one looks for statistical properties of the scattering response, such as the average field or intensity, intensity fluctuations, field-field or intensity-intensity correlations, etc [2]}. Theoretically, the problem is tackled in general by taking the configurational average of the relevant quantity (the field, intensity, etc.) from Eqs. (REF ) and (REF ), and making certain approximations on the interaction between particles.
