The skylight that we observe on a clear day is the result of the scattering of sunlight by the gas of the atmosphere. In the literature, two approaches have been used to account for this scattering. The first approach involves the direct summation of the sunlight scattered by the individual molecules of the gas and is the method originally used by Rayleigh. The second approach involves the scattering of sunlight by fluctuations in the relative permittivity of the gas and is the method developed by Smoluchowski and Einstein. In discussions of the origin of skylight, only one of these two approaches is generally followed, and it is sometimes stated or implied that the other approach is inaccurate or not applicable to the scattering of light by the atmosphere. In this paper, a simple model is adopted for the gas in the atmosphere. Then this model is used to obtain the irradiance for the sunlight scattered by the molecules contained in a coherence volume. The irradiance involves the statistical evaluation of a sum, and this sum is taken as the common starting point for investigating the two approaches mentioned before. Simple numerical calculations, based on a random number generator, are used to perform the statistical calculations, and in the end, the observed irradiances are shown to be the same for both approaches.
REFERENCES
The superscript sr is used to indicate that this is the “scattered radiated” field: the part of the scattered field that behaves as , where is the radial distance from the molecule.
The calculations reported here were performed using MATLAB with the random number generator rand.
The numbers used for , , and in Figs. 4 and 6 were chosen to make the graphs clearly illustrate the points we wished to make with regard to the distributions.
To construct the histogram, the range for is divided into small increments , and the ratio of the number of elements within bins of width to is plotted versus .
We are dealing with the fluctuations in the number of molecules, , in the slices of our equivalent volume, shown in Figs. 2(c) and 5. This volume was obtained from the coherence volume, see Fig. 2(b), by combining all regions with phase that differ by integer multiples of . Thus, the fluctuations are actually spread throughout the coherence volume.