$$\require{cancel}$$

11.5: Curve of Growth for Lorentzian Profiles

[ "article:topic", "authorname:tatumj" ]

The optical depth of a line broadened by radiation damping is given, as a function of wavelength, by

$\label{11.5.1}\tau (x)=\tau(0)\frac{l^2}{x^2+l^2},$

where the HWHM is

$\label{11.5.2} l=\frac{\Gamma \lambda_0^2}{4\pi c}$

and the optical thickness at the line center is

$\tau(0) = \frac{e^2\mathcal{N}_1 f_{12}}{m\epsilon_0 c\Gamma}.$

Here $$x = \lambda − \lambda_0$$, and the damping constant $$\Gamma$$ may include a contribution from pressure broadening. As in Section 11.4, in figure XI.4 I draw line profiles for optical thicknesses at the line center $$\tau(0)=\frac{1}{16},\frac{1}{8},\frac{1}{4},\frac{1}{2}, 1 , 2 , 4 , 8$$.

We see that the wings continue to add to the equivalent width as soon as, and indeed before, the central depth has reached unity. On combining equations 11.3.4, \ref{11.5.1},2 and 3, we obtain the following expression for the equivalent width:

$\label{11.5.4}W=\frac{\Gamma \lambda_0^2}{2\pi c}\int_0^\infty \left ( 1-\text{exp}\left \{-\frac{\tau(0)}{y^2+1}\right \}\right )\,dy,$

in which $$y = x/l$$ - i.e. distance from the line center in units of the HWHM.

If we now substitute $$y = \tan \theta$$, the expression for the equivalent width becomes

$\label{11.5.5}W=\frac{\Gamma \lambda_0^2}{2\pi c} \int_0^{\pi/2} \frac{1-\text{exp}\left [-\tau(0)\cos^2 \theta \right ]}{\cos^2 \theta }\,d\theta .$

Now that we have a finite upper limit, the expression can be integrated numerically without artificial and unjustified truncation. As described in the Appendix to Chapter 10, calculation of the trigonometric function cos can be avoided, and hence the integration much speeded up, by the substitution of $$t=\tan (\frac{1}{2}\theta)$$. Although the denominator of the integrand is obviously zero at the upper limit, so is the numerator, and the value of the integrand at the upper limit is finite and equal to $$\tau (0)$$. Figure XI.5 shows the equivalent width, in units of $$\frac{\Gamma \lambda_0^2}{2\pi c}$$ as a function of $$\tau(0)$$.

$$\text{FIGURE XI.4}$$    $$\text{FIGURE XI.5}$$