To verify the validity of the ISA and EFA models for the problem of interest, we use an in-house multiple-scattering code developed by Jean-Paul Hugonin at Laboratoire Charles Fabry (Institut d'Optique Graduate School, CNRS, Université Paris Saclay). The numerical method belongs to the broad family of T-matrix methods [1]}, that is probably the most adapted to solve wave scattering problems by discrete media. T-matrix methods benefit from the high degree of analyticity of the functions describing the wave propagation between particles, typically dyadic Green functions or vector spherical wave functions, depending on the problem of interest and the specific implementation. With the recent emergence of disordered media in photonics [2]}, the modelling of multiple scattering by particles in layered media has gained considerable attention, leading to the development of new methods [3]} and powerful publicly available tools [4]}.
