Results for aerosol scattering (aspherical particles) can be found in the following publication:
Wauben, W. M. F. and Hovenier, J. W.: Polarized radiation of an atmosphere containing randomly-oriented spheroids, J. Quant. Spectrosc. Radiat. Transfer, 47, 491–504, doi:10.1016/0022-4073(92)90108-G, 1992.
Two homogeneous layers above a Lambertian surface. Non-zero depolarization factor is assumed for the top (Rayleigh) layer:
Wauben, W.M.F., J.F. de Haan, and J.W. Hovenier, 1994: A method for computing visible and infrared polarized monochromatic radiation in planetary atmospheres. Astron. Astrophys. 282, 277-290.
Some results for I and Q for low orders of scattering, as well as sum of all scatterings are given in:
Wauben WMF, de Haan JF, and Hovenier JW, 1993: Low orders of scattering in plane-parallel homogeneous atmosphere. Astron.Astroph., 276, pp.589-602.
Two homegeneous layers (aerosol+Rayleigh)
Haan, J.F. de, P.B. Bosma, and J.W. Hovenier, 1987: The adding method for multiple scattering calcualtions of polarized light. Astron. Astrophys. 183, 371-391.
Reflected, transmitted and internal light field: Two types of aerosols with L=13 and 60 moments are considered in
Garcia, R.D.M., C.E. Siewert, 1989: The FN method for radiative transfer models that include polarization effects, . J. Quant. Spectrosc. Radiat. Transfer 41, No 2, 117-145.
Results for similar scenarios but with smaller number of significant digits are tabulated in:
Garcia, R.D.M., C.E. Siewert, 1986: A generalized spherical harmonics solution for radiative transfer models that include polarization effects, . J. Quant. Spectrosc. Radiat. Transfer 36, No 5, 401-423.
One of the two scenarios mentioned above, the L=13 case, is tabulated in:
Siewert, C. E., 2000: A discrete-ordinates solution for radiative-transfer models that include polarization effects. J. Quant. Spectrosc. Radiat. Transfer 64, 227-254.
Those working with Fourier decomposition of the light field, might find useful the following paper where results are given for several first Fourier moments separately:
Vestrucci, P., C.E. Siewert, 1984: A numerical evaluation of analytical representation of the components in a Fourier decomposition of the phase matrix for the scattering of polarized light. J. Quant. Spectrosc. Radiat. Transfer, 31, No 2, 177-183.