How did Lord Rayleigh derive/determine the phase function for his scattering model? I've been researching the question for quite some time, as I understand it the phase function is actually an approximation due to the particle-wave duality inherent in participating media such as the atmosphere or anything else.
As they are reemitted from particles, some EM waves mingle with their neighbours and amplify or kill off each other, therefore general approximate lobes are constructed from empirical data.
The Rayleigh scattering phase function is symmetrical as defined:
$$\Phi_R(\theta) ~=~ \frac{1}{4\pi}\frac{3}{4}(1 + \cos^2\theta).$$
However, I have no idea where this came from. I tried researching different phase functions, even those more modern from Henyey and Greenstein in 1941. But they just stated it is an approximation and gave the final form with no details.
Could someone show me, please?
 A: If the incident radiation is unpolarized, it can be seen as the sum of two mutually incoherent terms: the first one due to a polarization perpendicular to the scattering plane and the second one contained in this plane. For each term, the scattered radiation preserves the polarization.

*

*The first term does not vary with $\theta$, since the induced dipole is also perpendicular to the scattering plane, and its emission is constant along the equator


*The amplitude of the second one goes as $\cos (\theta)$ (and the intensity as $\cos ^2(\theta)$), because the dipole does not emit along its direction. See the last slide of this for an illustration.
Notes: a) the formula also holds for circularly polarized light. b) for other polarization states, the two terms must be weighted according to their contribution to the incident field.
A: It looks like there is a derivation at (http://irina.eas.gatech.edu/ATOC5235_2003/Lec9.pdf ) (the coefficient differs, but that may be due to a somewhat different definition of the  phase function).
