# Frequency Redistribution in Non-Coherent Scattering

### Reference Materials

- "Non-coherent scattering: I. The redistribution function with Doppler broadening (Hummer, 1962)"
- "Non-coherent scattering and the absence of polarization in Fraunhofer lines (Zanstra, 1946)"
- "Notes on the Theory of Noncoherent Scattering (Spitzer, 1944)"
- "Noncoherent Formation of Absorption Lines (Woolley, 1938)"
- "The Time-Relaxation of Resonance Line Profiles (Field, 1959)"
- "Noncoherent Scattering due to Collisions. I. Zanstra's Ratio of Coherent to Uncorrelated Noncoherent Scattering. (Edmonds, 1955)"
- "Note on the Zanstra Redistribution in Planetary Nebulae (Unno, 1952)"

### Need to Review

## Introduction

When a photon is scattered by an atom, it may undergo a change in frequency. Scattering events in which such **frequency redistribution** occurs are said to be "**non-coherent**." Hummer, Field, Spitzer, Mihalas, Avrett, and others have explored the physical causes for frequency noncoherence in the scattering of light by atoms and the effects of noncoherence arising from both the natural width of atomic levels and the thermal motion of particles in stellar atmospheres. Here we develop a mathematical formalism to generally describe the process of non-coherent scattering. We apply it to four physical cases of noncoherence.

## Formulation of the Redistribution Functions

The process of non-coherent **line** scattering can be described very generally by a redistribution function, . describes scattering due to a single electronic transition – it does not describe continuum scattering or relativistic scattering. Consider the picture of non-coherent scattering presented in Figure 1. The probability that an incident photon with frequency between and , traveling in direction into an element of solid angle would be scattered to a different frequency, between , in direction , into , is:

**This is the definition of a redistribution function.**

This probability is normalized to unity:

For convenience, we take all integrals over frequency and /or solid angle to be complete, unless otherwise specified.

Now, integrating over and reduces to . This is the normalized probability that the atom absorbs a photon with frequency between and traveling in direction into an element of solid angle , and re-emits it at *some* frequency in *some* direction. This defines the **Line Profile Functions** of the transition, which could have an angular dependence. For spatially symmetric absorption profiles, this quantity is given by the familiar function, :

By the definition of ,

The line profile function gives a transition’s normalized absorption probability as a function of . We can therefore define a transition’s frequency-dependent absorption coefficient, , in terms of the line profile function:

where is the line-center absorption coefficient of the transition.

## Redistribution and the Radiative Transfer Equation

How does non-coherent scattering add energy to or remove energy from a beam of radiation? We now consider the effect of a single, idealized non-coherent scatterer on the spectrum of a radiation field, measured along a line of sight. Following Figure 2, define the specific intensity of the field in direction before scattering to be . Let the detector lie along direction . The energy per volume per time per solid angle per frequency removed from the beam is given by . This follows explicitly from the definition of the absorption coefficient, .
What about light scattered non-coherently *into* the beam? To visualize that case, we follow Figure 3, where our detector lies along . Now, the energy per volume per time per solid angle per frequency **added** into the incident beam at frequency is related to the **amount of incident radiation at frequency and direction scattered into frequency and direction **.

In the formalism of redistribution functions,

where we have convolved the redistribution function and the specific intensity of the radiation field over frequency and solid angle to account for the continuity of the field.

This defines a "non-coherent scattering emissivity"

The constant of proportionality is related to the Einstein ’s (relative likelihoods) of other atomic transitions that may occur when an atom is excited to the level that produces the transition corresponding to . For example, in the case of Ly- (3 -> 1) scattering, an H atom in the ground state that is excited into the level may either undergo Ly- decay or H- decay. The relative likelihood of H- is given by

Therefore the relative likelihood that a transition from the state is Ly- is given by . Therefore, the non-coherent scattering emissivity is given by

The total emissivity of the medium is therefore the sum of the scattering emissivity and the purely "emissive" (non-scattering) emissivity . So the entire Radiative Transfer Equation for a medium undergoing non-coherent scattering is given by

where we have summed over all possible upper level destinations, .

## Functional Forms of The Redistribution Functions I: Two Simple Cases

We now seek specific functional forms of . **Often, but not always, we can decompose into the product of an angular function , and a frequency dependent function .** First, we consider two very simple cases that permit this simplification: **coherent scattering **, and **complete redistribution ( completely uncorrelated)**:

Further, we consider two subcases, corresponding to two distinct functional forms of .

**Isotropic Scattering**(e.g., the scattering of light by slow-moving dust.) In this case,

since incident radiation is equally likely to be reemitted in any direction.

**Dipole Scattering**(e.g., Thomson Scattering.) In this case,

since light is more likely to be scattered into the dipole field.

### Case I: Coherent Scattering

We consider the case of coherent scattering, where the frequency of incident radiation is equal to the frequency of scattered radiation . In this case, is given by the Dirac Delta Function, centered on . This is because there is a 100% probability that incident radiation is reemitted at its original frequency. We assume corresponds to either of the isotropic or dipole cases. Therefore, we can write down:

The scattering emissivity is therefore given by

If the scattering is isotropic, then , and becomes

### Case II: Complete Redistribution

Complete redistribution occurs when and are completely uncorrelated. In this case, is determined exclusively by :

If the scattering is isotropic, , so

### Was that realistic?

No. In reality, redistribution occurs **in the frame of the atom**. Additionally, the phenomenon of **Doppler redistribution** ensures that laboratory-observed redistribution is only partial. To account for these effects, we henceforth place ourselves into the rest-frame of the atom, making the following changes of variables:

## Non-Coherent Scattering with Doppler Redistribution

### Formalism

Let define the absorption profile (line profile) in the atom’s frame. We therefore have

From our earlier discussion of redistribution functions, it must be the case that

defines the probability that a photon with frequency between and is absorbed in solid angle . Now, let

define the frequency-dependent component of the linearly-separated redistribution function, in the atom’s rest-frame. Further, let

define the corresponding angular function. Bringing everything together, we have that

gives the probability that *if* a photon is absorbed at , it is re-emitted at . Therefore,

for a single particle with velocity .

## Three Physical Cases of Non-Coherent Scattering with Doppler Redistribution

We now consider three specific cases of non-coherent scattering in the rest-frames of individual atoms. Later, we will generalize our results to entire stellar atmospheres.

### Case I: The Infinitely Sharp Atom

The first case we examine is doppler redistribution in the rest-frame of a simple two-level atom with infinitely sharp (i.e. 0 width) energy states. In this case, there is **no frequency redistribution in the frame of the atom**. Instead, we have:

### Case II: Resonance Lines

Here we treat the case of an atom with a perfectly sharp ground state and a broadened upper state, whose finite lifetime against radiative decay leads to a Lorentzian line profile for the atom, consistent with natural broadening:

where is the total radiative decay rate (or "radiative damping width") of the atom, that is,

We further assume that there are no additional perturbations of the atom while it is in its upper state. Then there will be no reshuffling of electrons among substates of the upper state, and the decay back down to the lower state will produce a photon of exactly the same frequency as the one originally absorbed. Thus we have

This case applies to resonance lines in media of such low densities that collisional broadening of the upper state is completely negligible, such as the Ly- line of Hydrogen in the interstellar medium.

### Case III: Complete Redistribution

The basic physical picture here is of an atom with a perfectly sharp lower state, and a broadened upper state, in a medium where collisions are so frequent that *all* excited electrons are randomly reshuffled over the substates of the upper state before emission occurs. The absorption profile is again the Lorentz profile, where now represents the full width (radiative plus collisional) of the upper state. In this extreme limit, the frequency of the emitted photon will have *no* correlation with the frequency of the absorbed photon; the probability for emission at any particular frequency is then proportional to the number of substates present at that frequency, and hence to the absorption profile itself. When complete redistribution in the atom’s frame occurs we thus have

## Solving the Three Cases: Doppler Shift Redistribution in the Laboratory Frame

In this section, we shall consider the effects of the Doppler shifts introduced by the motion of the scattering atoms relative to the laboratory frame by deriving expressions that describe the full angular and frequency dependence of redistribution in the scattering process.

### A Thermal Atmosphere

Suppose an atom moving with velocity , which remains fixed during the scattering process, absorbs a photon () and emits a photon (), as measured in the laboratory frame. We may therefore express the photon frequencies in the atom’s frame according to:

and

The redistribution function for this scattering event may be rewritten as

This is the *specific* redistribution function for a single species of atoms with speed . To find the *net* result for the entire population of atoms, we must take the average over the velocity distribution, assumed to be Maxwellian. To perform this average, we make a change of coordinates.

Let lie in the plane. bisects . Now,

Write:

Then

Introduce a quantity called the reduced velocity, .

Call the doppler width .

Now, the Maxwell-Boltzmann distribution can be re-expressed:

The integral is trivial by our choice of coordinate system. It is .

### Cases 1 and 2: Coherent Scattering in the Rest Frame of the Atom

To do the integral, let . Then

Let .

So, finally, for coherent scattering:

#### Case I: An Atmosphere of Infinitely Sharp Atoms in the Laboratory Frame

Recall that for Case I, we established that

We therefore seek to evaluate:

Make the following change of variables:

Then the above integral becomes:

Therefore

Therefore, although the scattering is perfectly coherent in the frame of the atom, it is non-coherent in the laboratory frame, due to the doppler effect. Now, people often rewrite this result in terms of a new set of variables, and , which define a frequency shift relative to the doppler width:

**Failed to parse (syntax error): {\displaystyle \mathbb {R}_ I(x’, \hat{n}^\prime ; x, \hat{n}) = \mathbb {R}_ I(\nu ^{\prime }, \hat{n}^{\prime };\nu , \hat{n})\underbrace{\frac{d\nu ^\prime }{dx^\prime }\frac{d\nu }{dx}}_{\Delta ^2} \,\!}**

**Failed to parse (syntax error): {\displaystyle \mathbb {R}_ I(x’, \hat{n}^\prime ; x, \hat{n}) = \frac{g(\hat{n}^\prime , \hat{n})}{\pi \sin (\theta )}e^{-x^2-(x^\prime -x\cos (\theta ))^2\csc ^2(\theta )} \,\!}**

#### Case II: An Atmosphere of Resonance Line Atoms in the Laboratory Frame

Here we have

Therefore, the integral over the thermal velocities becomes:

This integrand is similar to the Voigt function, . Therefore,

**Failed to parse (syntax error): {\displaystyle R_{II}(x’, \hat{n}^\prime ; x, \hat{n}) = \frac{g(\hat{n}^\prime , \hat{n})}{\pi \sin (\theta )}e^{-\frac{1}{2}(x-x^\prime )^2\csc ^2(\frac{\theta }{2})}H\left(a\sec (\theta /2), \frac{1}{2}(x+x^\prime )\sin (\theta /2)\right) \,\!}**

where

and

We leave case III to the concerned reader!