Curve of Growth The curve of growth describes how the line strength increases as the optical depth increases. An appropriate measure of line strength is the equivalent width: the area of a rectangle with the same area as the line profile and height equal to the continuum. Webegin with apurelyabsorbing, simple slab modelwith radiation incident on one side. Then from the solution to the transfer equation, we have the intensity emerging: I ν = I e τn where the optical depth is τ n = n x κ φ v dl ' n x Lκ φ v where n x is the density of atoms of element x, L is the thickness of the slab, κ is the absorption coecent per atom at line center, and φ ν is the line profile function. The factor R n x dl N x ' n x L is called the column density The area under the curve I ν (ν) is R I ν dν and the area of the line is R (I I ν ) dν The equivalent width of the line is defined to be µ W = 1 I ν dν = 1 e N xκ φ v dν I Of course this is an idealization. In practice the limits are frequencies far enough from the line center that we are clearly in the continuum. In the classical oscillator model, ³ hν/kt κ = πe2 mc f 1 e where the oscillator strength f takes quantum e ects into account and ν is the frequency at line center. The line profile function depends on the broadening mechanisms in e ect. 1. We always have Doppler broadening due to the velocity distribution of the atoms. In LTE, " µ # 2 1 ν ν φ ν,doppler = p exp πνd ν D 1
and the line breadth factor is ν D = ν c 2kT M + (turbulent velocity)2 1/2 2. Natural broadening is due to radiation reaction (my notes, R&L pg 287, Shore pg 197-8) γ φ ν = 4π 2 (ν ν ) 2 + (γ/2) 2 This is sometimes called the Lorentz profile. We can include the e ect of collisions by adding a term γ C to the natural broadening damping factor γ..1 Weak lines N x κ φ 1 Then W = 1 1 e N xκ φ v dν ' N x κ φ v dν = N x κ φ v dν = N x κ The result is independent of the shape of φ ν and thus of the type of broadening. The equivalent width is linearly proportional to the column density..2 Strong lines τ = N x κ φ À 1 Notice that N x κ φ ν cannot be À 1 everywhere, since φ ν! far from the line center. Let the frequency where τ ν ' 1 be ν. Then W ' ν ν + 1 1 e N xκ φ v dν+ 1 e N xκ φ v dν+ ν ν + In the middle factor, e τν is 1. The two other terms are each of order N x κ ν + φ v dν N x κ 1 e N xκ φ v dν and we neglect them. Thus W ' 2 2
.2.1 Moderately strong lines In these lines N x κ φ ν becomes <1 while still in the Doppler core. So N x κ φ v (ν + ) = 1 " µ 2 # 1 N x κ p exp = 1 πνd ν D µ Nx κ = ν D sln p πνd where the square root is. 1 for self consistency, and it follows that µ Nx κ W ' 2ν D sln p πνd Thus the equivalent width increases logarithmically with N. Check that the terms we neglected are small: W neglected W kept So the approximations are consistent..2.2 Very strong lines N x κ ³ 2ν D rln pπνd. e Nx κ 2» 1 Here N x κ φ ν does not become <1 until the wings of the line. Thus the transition occurs at where N x κ γ 4π 2 2 + (γ/2) 2 = 1 Uusally γ 2, so we may approximate: p Nx κ γ = 2π and then p Nx κ γ W ' π and increases as the square root of N. Check as in the previous case. The neglected term is: N x κ γ 4π 2 x 2 + (γ/2) 2 dx = N xκ γ 4π 2 = N xκ γ 4π 4π 2 γ» N xκ γ 4π 2» 3 dx x 2 + γ 2 /16π 2 µ π 2 4π tan 1 γ p Nx κ γ 2π
Thus we get the curve of growth: let x = N xκ ν D and y = W/ν D then The transition to the saturated regime occurs sooner (lower W) for lower temperatures because ν D is smaller. The damped (square root) portion is higher is the damping paramater γ is larger. The line profile becomes squarer in the saturated portion abd devlops broad wings in the damping portion. (See figure on next page.) Touse the curve of growth one compares the theoretical curve with an emipiricalcurve contructed from the data. The column density and temperature can then be deduced. One can use line multiplets of the same atom, or di erent trancitions (but beware here if γ is dominated by naturalrather than collisional broadening.) Here s some data for hydrogen from http://nedwww.ipac.caltech.edu/level5/charlton/charlton1_1.html b is the Doppler parameter ν D in velocity units: b = p 2kT/m and m is the atomic mass. 4
Here s another one from http://www-int.stsci.edu/~leonidas/voss3/ 5
6
Q42-388 z=3.12 Lya Forest + metals Ly-limit Damped Lya QSO Lya em. 7