Optics

The cloud light has to walk through

A beam crossing thirty scattering lengths of anything should keep e⁻³⁰ of itself — a ten-millionth of a millionth. A cloud thirty scattering lengths thick lets through nearly a third of the sunlight falling on it. The light has not crossed; it has walked, one scattering at a time, and a walk through a slab obeys a law with the shape of Ohm's rather than of an exponential.

Assumes: Everything a scatterer removes, from one direction · When the particle is the size of the wave

Everything a scatterer removes found the total light a single particle takes out of a beam in one direction, straight ahead, and it ended by naming the thing that theorem cannot do. It is a statement about one particle in a plane wave. In a medium many scattering lengths thick, the light arriving at each particle has already been scattered by others, and nothing about a single particle’s forward amplitude says what comes out of the far side.

The question is not exotic. Why a litre of water is not blue observed that a cloud is white because its optical depth is tens or hundreds, so that every colour has certainly been scattered, and when the particle is the size of the wave noted that so many scatterings wash out a droplet’s own angular pattern. Both are correct statements about colour and pattern. Neither says how much light gets through a cloud, or why the answer is so much larger than an exponential would allow. That needs a different kind of calculation, and the kind of calculation it needs is the most useful fact in the subject.

A walk instead of a beam

The direct approach is to follow light through a slab one photon at a time. A photon travels a random distance, drawn from the exponential distribution whose mean is the scattering mean free path. It then scatters into a new direction drawn from the particle’s angular pattern, travels another random distance, scatters again, and so on until it leaves through the top or the bottom.

Light that walks through a cloud. Seven photons entering a slab 8 scattering mean free paths thick from above, straight down, and followed until they leave, with every scattering equally likely to send them in any direction; their paths are projected onto the page. Of 20000 photons followed the same way, 82.5 per cent come back out of the top and 17.5 per cent out of the bottom. Diffusion theory predicts 18.2 per cent through the bottom. Only 10 photons cross without scattering at all, where e⁻⁸ predicts 6.7 — within the counting error of so few: nearly all of what is transmitted has walked, and the direction it arrived from is forgotten on the way.
Fig. 1 Seven photons entering a slab eight scattering mean free paths thick from above, followed until they leave, with every scattering equally likely to send them in any direction; their paths are projected onto the page. Of 20,000 followed the same way, 82.5 per cent come back out of the top and 17.5 per cent out of the bottom, against 18.2 per cent from diffusion theory. Only 10 cross without scattering, where e8e^{-8} predicts 6.7.

The paths are tangles. A photon entering the top goes a mean free path or so, turns, turns again, and most of the time finds its way back out of the surface it came in through. The ones that reach the bottom do so after wandering sideways as far as they have travelled down. Of twenty thousand photons sent into a slab eight mean free paths thick, 17.5 per cent emerge from the bottom.

The number to compare that with is the fraction a straight beam would keep, e8e^{-8}, which is 0.034 per cent, and the photons that crossed without scattering at all are indeed about that fraction — ten of them, against 6.7 expected, which is inside the counting error of so few. Five hundred times more light comes out of the bottom than crossed without scattering. Nearly everything transmitted has walked, and has forgotten the direction it arrived from.

Ohm’s law for light

A photon’s walk through a thick slab is a random walk, and a crowd of random walkers obeys the equation that only runs forwards: the density of photons diffuses, with a diffusion coefficient set by how far a photon goes before forgetting its direction. Light entering the top of a slab is a source just inside the top surface; light leaving through either face is a sink. The steady state is a density that falls linearly from one face to the other.

A linear density profile has a definite consequence. The flux of photons through the slab is proportional to the slope of the profile, and the slope is inversely proportional to the thickness. So the transmitted fraction falls as one over the thickness, not exponentially — the relation between current and length in a wire, with the mean free path in the role of the conductivity. The correct form, from the diffusion equation with the surfaces handled properly, is

T=1+zeL/+2ze,T = \frac{1 + z_e}{L/\ell^* + 2z_e},

where \ell^* is the transport mean free path, introduced below, and ze=0.7104z_e = 0.7104 is the extrapolation length, measured in transport paths: the distance past each face at which the diffusing density would fall to zero if the medium continued. That 0.7104 comes from solving the exact transport problem at a surface — Milne’s problem — and it is the one number in the formula that a diffusion argument alone cannot supply.

The formula has a second reading that needs no flux at all. A walker diffusing between two absorbing walls leaves through the far one with a probability equal to its distance from the near one divided by the distance between them — the one-dimensional case of walkers stopping on a boundary to solve Laplace’s equation, where the answer is always the harmonic function and in one dimension the harmonic function is a straight line. A photon entering the top makes its first effective scattering about one transport length in. The walls it diffuses between are the two extrapolated faces, zez_e transport lengths outside the slab on either side. So it starts 1+ze1 + z_e from the top wall, the walls are L/+2zeL/\ell^* + 2z_e apart, and the chance of leaving through the bottom is the ratio: the transmission formula, term for term.

Transmission that falls as one over the thickness. The fraction of light getting through a non-absorbing slab, against its thickness in scattering mean free paths, on logarithmic axes, for particles that scatter equally in all directions and for cloud droplets, which scatter strongly forward with a mean cosine of 0.85. The dots follow photons; the solid curves are diffusion theory with the step lengthened to ℓ/(1 − g); the dashed curve is the light that crosses without scattering, e to the minus the thickness. With g = 0, 1 mean free path transmits 0.662, 8 mean free paths transmit 0.180, 30 mean free paths transmit 0.060; with g = 0.85, 1 mean free path transmits 0.957, 8 mean free paths transmit 0.637, 30 mean free paths transmit 0.290, 100 mean free paths transmit 0.107. Thirty mean free paths of cloud let through 29 per cent of the light, where a beam crossing unscattered would keep 9.4 × 10⁻¹⁴ of it. Diffuse transmission falls as the inverse of thickness, not exponentially: Ohm's law, with the scattering length playing the part of a conductivity.
Fig. 2 Fraction transmitted through a non-absorbing slab against its thickness in scattering mean free paths, on logarithmic axes, for particles that scatter evenly and for cloud droplets with a mean scattering cosine of 0.85. Dots follow photons; solid curves are diffusion theory; the dashed curve is light that crosses unscattered. Thirty mean free paths of cloud let through 29 per cent of the light, where an unscattered beam keeps 9.4 × 10⁻¹⁴.

On logarithmic axes the difference between the two laws is the difference between a straight line of slope minus one and a cliff. For particles that scatter evenly, eight mean free paths transmit 0.180 and thirty transmit 0.060: nearly four times the thickness, a third of the light. Diffusion theory, which adds the extrapolation lengths to both thicknesses, gives 0.182 and 0.054, a factor of 3.3, and the photons at thirty paths sit within their counting error of it. The unscattered beam over the same range falls by nearly ten orders of magnitude. Even a slab one mean free path thick is not described by the beam: it transmits 0.662 where e1e^{-1} is 0.368, so nearly half of what crosses it has already scattered at least once.

The correspondence with Ohm’s law is closer than an analogy. Electrons in a disordered metal also scatter elastically off impurities, forget their direction after a transport mean free path, and diffuse, and the conductance of a wire is proportional to that mean free path divided by the wire’s length for exactly the reason the transmission of a cloud is. A cloud is a resistor for sunlight, and the diffusion equation does not care whether what walks through it is a photon or an electron.

Forward scattering only stretches the step

Cloud droplets are not even scatterers. When the particle is the size of the wave computed their angular pattern and found it tipped strongly forward, with a mean scattering cosine — the asymmetry parameter gg — near 0.8 for droplets many wavelengths across. A photon scattered by such a droplet is usually deflected by only a few degrees, and needs several scatterings before its direction has been randomised.

That suggests a single correction, and it is the right one. If each scattering keeps a fraction gg of the photon’s forward momentum on average, the direction is forgotten after about 1/(1g)1/(1 - g) scatterings, so the distance over which it is forgotten — the transport mean free path — is

=1g.\ell^* = \frac{\ell}{1 - g}.

For g=0.85g = 0.85 that is 6.7 scattering lengths. The claim is then that a slab of forward scatterers behaves, once it is thick, exactly like a slab of even scatterers whose mean free path is \ell^*.

A forward-scattering cloud is an ordinary one with longer steps. Transmission through a non-absorbing slab against its thickness measured in transport mean free paths, ℓ/(1 − g), for phase functions with mean cosines of 0, 0.5, 0.85, 0.95, each point from 2,500 photons. A slab of droplets with g = 0.95 is twenty times as many scattering lengths thick as one with g = 0 at the same point on this axis. Once the slab is a few transport paths thick the curves lie on one another — the largest spread among them from four transport paths onward is 9.2 per cent — and on the diffusion curve, drawn solid. Thin slabs do not collapse: at half a transport path the transmissions are 0.80, 0.83, 0.85, 0.85, because a few forward scatterings carry light through before any walk begins. Far from the surfaces, only the distance over which a photon forgets its direction matters.
Fig. 3 Transmission against thickness measured in transport mean free paths, ℓ/(1 − g), for mean scattering cosines of 0, 0.5, 0.85 and 0.95. From four transport paths onward the curves lie on one another, within 9.2 per cent, and on the diffusion curve. At half a transport path they do not — 0.80, 0.83, 0.85, 0.85 — because a few forward scatterings carry light through before any walk begins.

The photons agree. Measured in transport lengths, slabs of four very different kinds of scatterer — from particles that scatter evenly to droplets whose mean cosine is 0.95, which need about twenty scatterings to forget a direction — transmit within nine per cent of one another from four transport lengths onward. The collapse is the similarity principle of radiative transfer, and it is why an atmospheric scientist summarises a cloud’s droplets by two numbers, an optical depth and an asymmetry parameter, and gets the brightness of the cloud right without knowing anything else about how the droplets scatter.

The collapse fails for thin slabs, and the way it fails is instructive. A slab half a transport length thick lets through more light when its scatterers are forward-peaked, because the few scatterings a photon undergoes all deflect it only slightly and it emerges still heading roughly downward. The walk has not had room to start, and the angular pattern still matters. Far from the surfaces only the distance over which a direction is forgotten counts; near them, everything about the pattern does.

White on top, grey underneath

The same calculation, applied to a cloud of water droplets with g=0.85g = 0.85 and no absorption, answers a question anyone who has flown through weather has asked.

Bright on top and grey underneath, from the same droplets. The fraction of sunlight a cloud of water droplets sends back up — which is how bright its top looks from above — and the fraction it lets through to the ground — how bright its base looks from below — for clouds of optical thickness 5, 10, 30, 100, with the droplets' forward-peaked scattering, g = 0.85, and no absorption. At 5: 24 per cent up, 76 per cent down; at 10: 42 per cent up, 58 per cent down; at 30: 71 per cent up, 29 per cent down; at 100: 90 per cent up, 10 per cent down. Nothing is absorbed, so every percentage point missing from below is reflected above: a thick cloud is dark underneath because it is white on top. With a tenth of a per cent of absorption per scattering, the thickest cloud absorbs 18 per cent of the light, because a photon crossing it scatters hundreds of times on the way.
Fig. 4 The fraction of sunlight a cloud of water droplets sends back up — how bright its top looks from above — and lets through to the ground — how bright its base looks from below — at optical thicknesses of 5, 10, 30 and 100, with no absorption. At 5, 24 per cent up and 76 down; at 10, 42 and 58; at 30, 71 and 29; at 100, 90 and 10. With a tenth of a per cent of absorption per scattering, the thickest cloud absorbs 18 per cent in all.

A thin cloud, optical depth five, lets three quarters of the sunlight through and reflects a quarter, which is why a veil of cloud dims the Sun without darkening the day much. At an optical depth of thirty — an ordinary overcast deck a few hundred metres thick — the numbers have reversed, 71 per cent reflected and 29 transmitted. At a hundred, a deep storm cloud, the base receives a tenth of the light falling on the top.

Nothing in that calculation absorbs, so every per cent missing from the base has gone upward. The dark underside of a thick cloud is the same light as the dazzling white of its top seen from a satellite, and a cloud’s darkness from below is a measure of how bright it is from above. That arithmetic is one of the largest terms in the Earth’s energy balance: clouds cover about two thirds of the planet, and the fraction of sunlight they send back to space before it reaches the surface is set by the transmission law in the figure. Clouds supply roughly half of the Earth’s albedo of about 0.3, the number that sets the temperature a planet radiates at, and they send back something like 47 watts on every square metre of the planet averaged over day and night; a change of under ten per cent in that is a push on the climate as large as doubling carbon dioxide.

The base is grey rather than simply dim for a second reason the walk explains. The light emerging from the bottom of a thick slab has forgotten where the Sun is. It comes out nearly equally in all downward directions, so the underside of an overcast sky is evenly lit with no hint of the Sun’s position, and shadows on the ground below are soft or absent. The same forgetting is why the light’s polarisation is washed out after passing through fog: a photon that has lost its direction has lost its plane of polarisation along with it.

Why a trace of absorption matters

The last number in the cloud figure is the most surprising, and it follows from counting scatterings.

A photon that crosses a cloud scatters as the square of its thickness. The average number of times a photon is scattered on its way through a non-absorbing slab, counting only the photons that get through, against the slab's effective thickness in transport mean free paths — its thickness plus the extrapolation length of 0.7104 diffusion theory adds at each face — on logarithmic axes. With g = 0 the count grows as the 2.08 power of that thickness from four transport paths onward; with g = 0.85 the count grows as the 2.06 power of that thickness from four transport paths onward: the square a random walk needs to cover a distance. Fitted against the bare thickness instead, the same counts give 1.76 and 1.74, because the faces add a fixed length that matters less the thicker the slab. The forward-scattering droplets take about 1/(1 − g) ≈ 6.7 times as many scatterings for each transport step. A photon getting through 16 transport paths of cloud droplets scatters 1021 times on average, which is why a material almost perfectly transparent per scattering can still darken a thick cloud: its absorption is applied at every one of those events.
Fig. 5 Scatterings per transmitted photon against the slab’s effective thickness — its thickness in transport paths plus the 0.71 diffusion adds at each face — on logarithmic axes. The count grows as the 2.08 power for even scatterers and the 2.06 power for cloud droplets, where the bare thickness would give 1.76 and 1.74. A photon getting through 16 transport paths of droplets scatters 1,021 times.

A random walker needs a number of steps proportional to the square of the distance it has to cover — the law that makes light take so long to leave the Sun, whose interior is a ball of plasma some tens of billions of free paths in radius — and the photons confirm it: the scatterings per transmitted photon grow as the square of the slab’s effective thickness, with a measured exponent of 2.08 and 2.06. Measured against the bare thickness the same counts give 1.76 and 1.74, because each face adds a fixed extrapolation length that matters less the thicker the slab — a reminder that the thickness diffusion sees is not quite the thickness a ruler measures. Forward-scattering droplets take 1/(1g)1/(1 - g) times as many scatterings for each transport step, so a photon crossing sixteen transport lengths of cloud scatters about a thousand times.

A thousand scatterings multiply any absorption a thousandfold. Water absorbs almost nothing in the visible, but in the near infrared, at wavelengths where a droplet absorbs a tenth of a per cent of the light it scatters, a cloud with an optical depth of a hundred absorbs eighteen per cent of what falls on it. Weak absorption in a strongly scattering medium is amplified by the length of the walk, and satellite instruments exploit exactly this: comparing a cloud’s brightness at a wavelength where water absorbs slightly with one where it does not measures the droplets’ size, because larger droplets absorb more per scattering.

The amplification runs through every milky medium. A white paint contains pigment particles that scatter strongly and absorb weakly, and a small trace of absorbing contaminant greys it far more than its concentration suggests. Near-infrared light can penetrate several centimetres of living tissue because tissue scatters strongly and, between the absorption bands of water and haemoglobin, absorbs weakly; the light that emerges has walked, and clinical instruments that measure blood oxygen through the skin read how much of that walk was absorbed.

A cone that interference adds back

The diffusion picture treats light as particles and discards phase, and there is one place where that is measurably wrong — and wrong in a way that connects this cloud to the electrons in a wire once more.

Consider light scattered back out of the top of a thick slab along a closed-looking path: in at one point, a sequence of scatterings, out at another. The same sequence of scatterings traversed in the opposite order is a second path with exactly the same length. In exactly the backward direction the two paths leave in phase and interfere constructively, doubling the intensity of that pair. Away from exact backscattering they drift out of phase. So a white slab lit by a laser sends back a narrow cone around the exact backward direction, up to twice as bright at its peak as the diffuse background, with an angular width of about the wavelength divided by the transport mean free path.

That enhanced backscattering was observed in suspensions of particles in 1984 and 1985, and it is the optical version of weak localisation in metals, where the same pair of time-reversed paths makes a disordered conductor’s resistance slightly higher than the classical Drude value. The cone is interference surviving a random walk, and it is the first correction diffusion misses. Its width measures the transport mean free path directly, which is how \ell^* is measured in a material too dense to follow photons through.

Where the walk stops being a diffusion

The slab must be several transport lengths thick. Diffusion describes the walk, and a walk needs room. Below a few transport lengths the transmission depends on the angular pattern and the diffusion curve overestimates it, which the photons show for every phase function at half a transport length.

The scatterers are independent. The mean free path is computed as though each particle scattered alone. In a densely packed medium — a powder, a white pigment at high concentration — the particles’ positions are correlated, their scattered waves interfere, and the transport length departs from what one particle would give. That is why a paint does not simply get whiter with more pigment.

The geometry is a plane-parallel slab. A real cloud has sides, gaps and towers, and light leaks out of its edges and finds its way down through its thin places. Averaging a broken cloud field as a uniform slab overestimates its reflection, a correction that matters in climate models and that the one-dimensional calculation cannot supply.

The phase function is Henyey and Greenstein’s. The photons here scatter according to a one-parameter pattern that has the right asymmetry and the wrong details; a real droplet’s pattern has a diffraction spike, a rainbow and a glory. In a thick cloud the similarity collapse makes the details irrelevant to the brightness, which is the point of the figure; for anything that depends on the angle of the emerging light, such as the halo of a thin cloud around the Sun, they are everything.

And phase has been discarded. Every photon here is a particle. The coherent backscattering cone, the speckle a laser makes on a white wall, and the whole question of whether disorder can trap light are wave effects the walk contains none of.

What the photon paths hide

The paths figure draws seven walks projected onto the page, and a projection flattens three dimensions of wandering into two; the sideways excursions a photon makes perpendicular to the page are invisible, which makes every path look more tangled across and less tangled in depth than it is. And seven photons say nothing statistical: the reflected and transmitted fractions come from twenty thousand others that are not drawn.

The transmission and cloud figures report fractions and not where the light goes once it leaves. A cloud sends its reflected light back in all upward directions, brighter near the edges of its top for some geometries and near the centre for others, and it is that angular distribution — not the total — that a satellite looking straight down actually measures. The total is the number the energy balance needs; the angular distribution is the number an instrument sees, and the figures give only the first.

Still open: whether disorder can stop light

For electrons, enough disorder does more than raise the resistance. Past a critical strength the interference of multiply scattered paths traps them entirely, and a disordered metal becomes an insulator — Anderson localisation. The backscattering cone is the first hint of that interference; localisation is the same effect grown until diffusion stops.

Whether light can be localised the same way in three dimensions is an open question with an uncomfortable history. Claims to have observed it in strongly scattering powders were later traced, at least in part, to absorption, which also makes transmission fall faster than one over the thickness and is hard to distinguish from localisation. Calculations for dense collections of point-like scatterers suggest that the vector nature of light and the near-field coupling between closely packed scatterers may prevent three-dimensional localisation of light altogether, while localisation of sound and of matter waves in three dimensions has been reported. What a disordered material would have to be, if anything, to stop light by interference alone — and how an experiment could tell that from absorption — has not been settled.

The habit worth carrying away is to ask whether a quantity is crossing a medium or walking through it. An exponential describes what has not yet been scattered; an inverse law describes what has forgotten where it was going, and in a thick enough medium the second is almost all of what arrives.

Part 5 of 6

This essay is one argument about Scattering. The others:

The objects named here

The third axis, after the field and the reading path: the things themselves, and every essay that touches each one.

AlbedoAsymmetry parameterCoherent backscatteringDiffusionMean free pathMultiple scatteringOptical depthRadiative transferRandom walk