# Masers

### Short Topical Videos

### Reference Materials

- Population Inversion (Wikipedia)
- Astrophysical Masers (Wikipedia)
- "Mega-Masers and Galaxies"(K.Y. Lo, Annual Review)
- "Astrophysical Masers" (Gray, Malcolm, RAS review)

### Need a Review?

<latex>
\documentclass[11pt]{article}
\def\inv#1Template:1 \over
\def\ddtTemplate:D \over dt
\def\mean#1{\left\langle {#1}\right\rangle}
\def\sigot{\sigma_{12}}
\def\sigto{\sigma_{21}}
\def\eval#1{\big|_{#1}}
\def\tr{\nabla}
\def\dce{\vec\tr\times\vec E}
\def\dcb{\vec\tr\times\vec B}
\def\wz{\omega_0}
\def\ef{\vec E}
\def\ato{{A_{21}}}
\def\bto{{B_{21}}}
\def\bot{{B_{12}}}
\def\bfieldTemplate:\vec B
\def\apTemplate:A^\prime
\def\xp{{x^{\prime}}}
\def\yp{{y^{\prime}}}
\def\zp{{z^{\prime}}}
\def\tp{{t^{\prime}}}
\def\upxTemplate:U x^\prime
\def\upyTemplate:U y^\prime
\def\e#1{\cdot10^{#1}}
\def\hf{\frac12}
\def\^{\hat }
\def\.{\dot }
\usepackage{fullpage}
\usepackage{amsmath}
\usepackage{eufrak}
\begin{document}
\subsection*{ Cosmic Masers }

$OH$ and $H_20$ masers can occur in dusty, star-forming regions which are cold enough for these molecules to form. The dust's black-body radiation in the infrared band is absorbed by these molecules and a population inversion is established. When maser emission is caused by via stimulated emission, these clouds can get very bright (brightness temperatures $\sim 10^{14}K$). Temperatures this high cannot be thermal, so we know a maser when we see one. \par Generally, masers are useful for tracing the galaxy's magnetic field (emission lines are Zeeman split), and for following disks of gas and dust around stars in star-forming regions. In order to detect them, we need clouds which are moving uniformly together, and have velocity coherence both $\perp, \|$ to our line-of-sight through a disk of rotation gas. In general, the intensity we observe depends on the path length through the masing cloud, so we like long path lengths with the same velocity.

\subsection*{ How masers Work } We will disucss population inversion a bit later, for now lets just assume that there are many more atoms in state 2 then state 1, where state 2 has higher energy than state 1. If a photon with an energy equal to E2-E1 comes along and interacts with an atom, it will thus be significantly more likely to produce stimulated emission (ie cause the atom to emit an identical photon with energy E2-E1 and drop from state 2 to 1) than to radiatively excite the atom (ie absorb the photon and transition from state 1 to 2). We will ignore the possibility of radiative excitation from state 1 to 2 entirely. The duplicate pair of photons produced by stimulated emission then go on to interact with other atoms leading to an exponential increase in the number of photons (and thus the intensity of radiation). Now onto the math. Consider a molecule with two rotational energy levels. Observing a homogeneous slab of this molecule, the intensity we receive is given by the familiar (see Radiative Transfer Equation): $$I_\nu=S_\nu(1-e^{-\tau_\nu})$$ where $S_\nu={j_\nu\over\alpha_\nu}$. $j_\nu$ and $\alpha_\nu$ are given by (see Einstein Coefficients): $$\begin{aligned}j_\nu&=\overbrace{n_2\ato}^{per time}\overbrace{h\nu}^{per E} \overbrace{\phi(\nu)}^{per\ Hz}\overbrace{\inv{4\pi}}^{per steradian}\\ \alpha_\nu&={h\nu\over4\pi}\phi(\nu)[\overbrace{n_1\bot}^{abs}-\overbrace{ n_2\bto}^{stim\atop emis}]\\ \end{aligned}$$ Thus, our source function looks like: $$S_\nu={n_2\ato\phi(\nu)h\nu\inv{4\pi}\over{h\nu\over4\pi}\phi(\nu)[n_1\bot- n_2\bto]}$$ Then since $g_1\bot=g_2\bto$ and ${\ato\over\bto}={2h\nu^3\over c^2}$, we have: $$S_\nu={2h\nu^3\over c^2}\inv{{n_1g_2\over g_1n_2}-1}$$ A population is said to be inverted when $n_1g_2<g_1n_2$ (not $n_1<n_2$). $n_1\over g_1$ is an expression for the population per degenerate sub-level in energy level 1. If ${n_1\over g_1}<{n_2\over g_2}$, then $S_\nu<0$, and we have a maser. Expressed in terms of the excitation temperature (${n_2\over n_1} ={g_2\over g_1}e^{-{h\nu\over kT_{ex}}}$), we have: $$S_\nu={2h\nu^3\over c^2}\inv{e^{h\nu\over kT_{ex}}-1}$$ which is less than 0 when $T_{ex}<0$. We can express the optical depth of this slab to maser radiation as: $$\tau_\nu=\alpha_\nu L ={h\nu\over4\pi}\phi(\nu)\bto\left[{n_1g_2\over g_1}-n_2\right] L ={h\nu\over4\pi}\phi(\nu)\bto n_2\left[{n_1g_2\over g_1n_2}-1\right] L$$ If ${n_2\over g_2}>{n_1\over g_1}$, then $\tau_\nu<0$. Now you might think we're talking nonsense with a negative source function and a negative optical depth, but we're not. Only observable quantities need to be positive, since optical depth and the source function are just mathematical record keeping objects, they can take on any values.

Recall that the intensity is:

$$I_\nu=S_\nu(1-e^{-\tau})$$ If $\tau\to-\infty$, then $I_\nu\to-S_\nu e^{\tau_\nu}$, so $I_\nu\gg1$. On the other hand, if $\tau<0$ and $|\tau|\ll1$, then $I_\nu\to S_\nu\tau_\nu$, which is the product of two negatives = positive. This should be convincing you that $I_\nu$ is always positive, and therefore, actually manifested.

\subsection*{ Maser Species }

$OH$ mases at around 18 cm, and is found around AGB (asymptotic-giant-branch) stars in star-forming regions and around the galactic nucleus. AGB's are important because they have lots of dust. The two masing transitions are from $1667\to1612MHz$ and $1720\to1665MHz$.\par $H_2O$ mases at $1.35 cm$ in a transition to its ground rotational state. Note that an order-of-magnitude calculation of $\Delta E={\hbar\over 2I}$ for $H_2O$ gives us an estimate of $\sim 1 mm$, which is incorrect. The correct transitional energy is caused by a slight degeneracy in water molecules.\par $SiO$ mases at $3.4 mm$ in its ground vibrational state, and is typically found in star-forming regions and around AGB stars.\par Other molecules mase, and there is even potential for detecting atomic lasers around massive stars ($L\sim10^5 L_\odot$). Hydrogen transitions from $10\to9$ have been observed (at $55\mu m$), and Stielnitski 1996 claims to have observed a population inversion in atomic H.

\subsection*{ Saturated vs. Unsaturated masers }

There are two modes of operation for masers. For unsaturated masers, the gain is exponential with the path length, and for saturated ones, the gain only grows linearly with path length. The masers we find in the cosmos are typically saturated. The following is working toward understanding why there are two modes in masers. Before we start with the math lets imagine a sequence of events that might lead to saturation. Initially we have lots of atoms in state 2. If a photon that has energy E2-E1 interacts with one of these atoms we get two identical photons. Both of these photons can then interact with other atoms in state 2 giving 4 photons. This leads to an exponential amplification of the number of photons. If this increasing number of photons continues to propogate through the atoms, a point will eventually be reached where there are more photons with energy E2-E1 then there are atoms in state 2 for them to interact with. This means that we transition from an exponential amplification (from photon doubling) to a linear amplification (from continuing to travel through a medium with some density of atoms in state 2). This is called saturation. The graph given below shows intensity $I_\nu$ as a function of optical depth $\tau$. The discontinuity in the derivative on the graph is the saturation point.

Now let us prove this intuitive solution mathematically.

We begin with the expression for specific intensity from

emissivity (see Radiative Transfer Equation and Einstein Coefficients):
$$\begin{aligned}{dI_\nu\over dz}&=j_\nu-\alpha_\nu I_\nu\\
&={h\nu\over4\pi}\phi(\nu)n_2A_{21}-{h\nu\over4\pi}
\phi(\nu)[n_1\bot-n_2\bto]I_\nu\\ \end{aligned}$$
In principle, the Line Profile Functions governing spontaneous emission and the
line-profile governing stimulated emission might not have to be the same, but
evidence seems to suggest they are. Let's assume they are, and rewrite this:
$$\begin{aligned}{dI_\nu\over dz}&={h\nu\over4\pi}\phi(\nu){n_2\ato\over g_2}\cdot g_2
-{h\nu\over4\pi}\phi(\nu)g_2\bto\left[{n_1\over g_1}-{n_2\over g_2}\right]
\end{aligned}$$
Defining $N_1\equiv{n_1\over g_1}$, $N_2\equiv{n_2\over g_2}$, $A\equiv
\ato g_2$, and $B=\bto g_2$, our equation looks like:
$${dI_\nu\over dz}={h\nu\over4\pi}\phi(\nu)[N_2A+B(N_2-N_1)I_\nu]$$
Now we'll integrate over frequency. $\phi(\nu)$ is a sharply peaked function,
so we'll treat it as a $\delta$-function. Then substituting
$I\equiv\int{\phi(\nu)I_\nu d\nu}$, we have:
$${\int{I_\nu d\nu}\over\Delta\nu}\equiv\int{I_\nu\phi(\nu)d\nu}$$
which is true by definition of $\Delta\nu$: $\int{I_\nu d\nu}=I\Delta\nu$.
So now we have:
$${dI\over dz}={h\nu\over4\pi\Delta\nu}[(N_2-N_1)BI+N_2A]$$
This is an equation with three unknowns. To close this system, we'll use the
equations of statistical equilibrium, which say:
$${dN_2\over dt}=-N_2BJ+N_1BJ-N_2A+R_2(N-N_1-N_2)-\Gamma_2N_2$$
where $J\equiv\int{I_\nu\phi(\nu)d\nu}$, and $J_\nu=\inv{4\pi}\int{I_\nu
d\Omega}$. $R_2$ is the ``pumping rate*. It describes the rate at which*
molecules not in state 1 or 2 (counted by $N-N_1-N_2$) are radiatively or
collisionally
knocked into state 2. $\Gamma_2$ is the loss rate of molecules in state 2 into
any state other than state 1. In the above equation, we've neglected two
terms: collisional excitation $1\to2$ and $2\to1$. These are crucial terms,
being tightly related to local thermal equilibrium. However, we will neglect
them to simplify analysis. We also have a similar equation for state 1:
$${dN_1\over dt}=N_2BJ-N_1BJ+N_2A+R_1(N-N_1-N_2)-\Gamma_1N_1$$
We'll further simply matters by setting the loss-rates for the two populations
equal to each other ($\Gamma_1=\Gamma_2=\Gamma$). Now let's solve for the
steady-state solution (${dN\over dt}=0$).
$${d(N_2-N_1)\over dt}=-(N_2-N_1)2BJ-2N_2A+(R_2-R_1)(N-N_{12})-\Gamma(N_2-N_1)$$
where $N_{12}=N_1+N_2$.
Two terms in this equation are acting to reduce the population inversion of
state 2 with respect to state 1: $(N_1-N_1)2BJ$ and $\Gamma(N_2-N_1)$. In the
{\it unsaturated regime} where $BJ\ll\Gamma$, stimulated emission is a minor
perturbation to the inverse. In this case we can solve the system of equations:
$$\begin{aligned}{d\Delta N\over dt}={d(N_2-N_1)\over dt}&=-2N_2A+(R_2-R_1)(N-N_{12})
-\Gamma\Delta N=0\\
{dN_{12}\over dt}&=(R_2+R_1)(N-N_{12})-\Gamma(N_1+N_2)=0\\ \end{aligned}$$
Then solving for $N-N_{12}$:
$$N-N_{12}={\Gamma(N_1+N_2)\over R_1+R_2}={\Gamma N_{12}\over R_1+R_2}$$
Substituting this into the difference equation, which reads:
$$-2N_2A={R_2-R_1\over R_1+R_2}\Gamma N_{12}-\Gamma\Delta N=0$$
the using $N_2={N_{12}+\Delta N\over2}$ (this is just an identity), we get:
$$\underbrace{\left(1+{\Gamma\over A}\right)}_{\equiv2\beta}\Delta N
=N_{12}\left[\underbrace{{R_2-R_1\over R_1+R_2}{\Gamma\over A}}_{\equiv
\inv{\alpha}}-1\right]$$
Thus:
$$\Delta N={N_{12}\over2\beta}\left[\inv{\alpha}-1\right]={N_{12}(1-\alpha)\over
1\alpha\beta}$$
Now we'll define $S\equiv{\alpha\beta\over1-\alpha}$. $S$ is meant to connote
the source function here. The reason we might expect $S$ to be related to
the source function is that $S_\nu\propto{N_2\over N_2-N_1}\propto{N_2\over
\Delta N}$, so $\Delta N\propto {N_{12}\over2S}$. This is why we use an $S$
here. Getting back to our original equation, we have:
$${4\pi\Delta\nu\over h\nu A}{dI\over dz}=\Delta N\underbrace{BI\over A}_{
\equiv\mathfrak{I}}+N_2$$
Using our newly defined $S$, this becomes:
$${4\pi\Delta\nu\over h\nu B}{d\mathfrak{I}\over dz}={N_{12}\over2S}\mathfrak{I}+
\hf(N_{12}+\Delta N)$$
Finally, defining an unsaturated gain length $L\equiv{4\pi\Delta\nu\over Bh\nu}
{2S\over N_{12}}$, we have:
$$L=\mathfrak{I}+S+\hf$$
Then using $ds={dz\over L}$, we have:
$${d\mathfrak{I}\over ds}=\mathfrak{I}+\hf+S$$
To order of magnitude:
$$\begin{aligned}\mathfrak{I}&\equiv{B\over A}I={B\over A}B_\nu(T_{bright})\\
&={c^2\over2h\nu^3}{2h\nu^3\over c^2}\inv{e^{h\nu\over kT_{bright}}-1}\\
&=\inv{e^{h\nu\over kT_{bright}}-1}\\ \end{aligned}$$
Now since $1\ll{kT_{bright}\over h\nu}$ (typical maser wavelengths are in
mm and $T_{bright}$ is typically several K), we can throw away our $\hf$ in
our equation for ${d\mathfrak{I}\over ds}$:
$${d\mathfrak{I}\over ds}=\mathfrak{I}+S\Rightarrow
\int{d\mathfrak{I}\over\mathfrak{I}+S}=\int{ds}\Rightarrow
\mathfrak{I}+S=De^s$$
Choosing $\mathfrak{I}=0, s=0$ (we're assuming there is no background source), then
$D=S$, so:
$$\boxed{\mathfrak{I}=S(e^s-1)}$$
If $s\ll1$, we have $\mathfrak{I}=Ss$ from spontaneous emission, and this is the
saturated case. If $s\gg1$, the $\mathfrak{I}=Se^s$ from stimulated emission, and
this is the unsaturated case. Earlier we threw away the $BJ$ term, but we
could have done all of this including that term, and the algebra would have
been the same. If we do this, we get the answer:
$$\boxed{{d\mathfrak{I}\over ds}={\beta(\mathfrak{I}+\hf)\over\beta+\mathfrak{J}}+S}$$
where $\mathfrak{J}$ is the non-dimensionalized, integrated flux,
$\mathfrak{J}={BJ\over A}$.

\subsection*{ The Saturated Mode }

Recall our masing equation: $${d\mathfrak{I}\over ds}={\beta(\mathfrak{I}+\hf)\over\beta+\mathfrak{J}}+S$$ In the unsaturated case, ${\beta\over\beta+\mathfrak{J}}\to1\iff\beta\gg\mathfrak{J}$. This is equivalent to saying: $$\hf\left(1+{\Gamma\over A}\right)\gg{BJ\over A}$$ And if $\Gamma\gg A$, then: $$A+\Gamma\gg BJ\iff\Gamma\gg BJ$$ In the saturated case, then $\beta\ll\mathfrak{J}$. In this case: $${d\mathfrak{I}\over ds}={\beta\mathfrak{I}\over\mathfrak{J}}+S={4\pi\beta\mathfrak{I}\over \int{\mathfrak{I}ds}}+S$$ where we used that $\mathfrak{J}\sim\inv{4\pi}\int{\mathfrak{I}d\Omega}$. If we consider a single beam of photons through a cloud, then $\int{\mathfrak{I}d\Omega}\approx \mathfrak{I}\Delta\Omega$, so: $$\boxed{ \begin{aligned}{d\mathfrak{I}\over ds}&={4\pi\beta\mathfrak{I}\over\mathfrak{I}\Delta\Omega}+S\\ &={4\pi\beta\over\Delta\Omega}+S\\ \end{aligned}}$$

\subsection*{ How We Get Population Inversions } \n There is another dichotomy in masers: ones which are excited radiatively and those which are excited collisionally. We'll discuss a cloud of molecules which have 3 energy states which are populated by {\it simple collisional pumping}. Before we dive into the math lets consider conceptually what is going on. The diagram shows a simple 3 state atom. We can only go between state E1 and E3 collisionally (ie they do not communicate radiatively). State E3 spontaneously decays radiatively to state E2 with a short lifetime. State E2 radiatively decays to state E1, but with a very long lifetime. This means that there can be a large number of atoms that are in state E2. Thus we have the requisite population inversion for masing to occur. Now onto the math.

The rate of change of the population of energy state

1 in this case is given by (sources - sinks):
$${dn_1\over dt}=n_2\ato+n_2C_{21}+n_3C_{31}+n_2\bto J_{21}-n_1\bot J_{12}
-n_1C_{12}-n_1C_{13}$$
Then $J_{12}=J_{21}$, and because of local collisional thermal equilibrium,
$C_{12}=C_{21}{g_2\over g_1}e^{-{E_{21}\over kT}}$. Similarly,
$C_{13}=C_{31}{g_3\over g_1}e^{-{E_{21}\over kT}}$. So dividing by $g$, and
defining $N_1\equiv{n_1\over g_1}$, we have:
$$\begin{aligned}{dN_1\over dt}=&\underbrace{n_2\over g_2}_{N_2}\underbrace{
{\ato\over g_1}g_2}_{\ato}+\underbrace{n_2\over g_2}_{N_2}\underbrace{
{C_{21}\over g_1}g_2}_{\equiv C}+\underbrace{n_2\over g_3}_{N_3}\underbrace{
{C_{31}\over g_1}g_3}_{\equiv C}+\underbrace{n_2\over g_2}_{N_2}\underbrace{
{\bto\over g_1}g_2}_{\bto}J_{21}\\
&-N_1\underbrace{\bto}_{\bto{g_2\over g_1}\equiv
\bto}-N_1\underbrace{C_{21}{g_2\over g_1}}_{\equiv C}e^{-{E_{21}\over kt}}-
N_1\underbrace{C_{31}{g_3\over g_1}}_{\equiv C}e^{-{E_{31}\over kT}}\\
\end{aligned}$$
Phew. Notice that we set $C_{21}=C_{31}$. This is just to make our lives
easier. We can do the same for $dN_2\over dt$, but omitting $B_{23}$, because
we're deciding not to have absorptions from $2\to3$ and no stimulated emission
from $3\to2$. Then our total population is $N=N_1+N_2+N_3$. Without being
careful, instinct tells us that in steady state, we'll have a population
inversion ${N_2\over N_1}>1$ if $C\ll A_{32}$. This
instinct is correct, but let's do this carefully. First we'll make some
assumptions:
\begin{itemize}\item $\ato\ll C$. For $H_2O$:
$$\ato\sim10^8s^{-1}\left({1216\AA\over1.35cm}\right)^3
\sim10^8\e{-15}s^{-1}
\sim10^{-7}s^{-1}\left({\mu\over ea_0}\right)^2$$
where $\mu$ is our way of accommodating the fact that the dipole moment for
$H_2O$ might not be the same as for the fine structure of hydrogen. It turns
out the answer is $\ato\sim2\e{-9}s^{-1}$.\par
Estimating C:
$$C\sim n_{H_2}\sigma v_{rel}\sim n_{H_2}\e{-15}\cdot(3{km\over s})
\sim n_{H_2}\cdot3\e{-10}$$
which is $\gg10^{-7}s^{-1}\left({\mu\over ea_0}\right)^2$ when
$n_{H_2}\gg10^3cm^{-3}\left({\mu\over ea_0}\right)^2$.
\item Next we'll assume $E_{21}\ll kT$. Define ${E_{21}\over kT}\equiv\delta
\ll1$.\end{itemize}
Now we have a 2-step:
\begin{itemize}\item Step 1: Since $1\to3$ are not linked by radiation,
$${N_3\over N_1}\approx e^{-{E_{31}\over kT}}\equiv\theta\le1$$
\item Step 2: Radiative decays from $3\to2$. To get an inversion, we'll
argue that the sources into 2 are larger than the sinks out of 2 (this is
a little weird because we're solving our steady-state equations, but whatever):
$$N_3A_{32}+N_3C+N_1Ce^{-{E_{21}\over kT}}>N_2Ce^{-{E_{32}\over kT}}+N_2C$$
We can rewrite this as:
$$\begin{aligned}
\underbrace{N_2\over N_1}_{=\theta}{A_{32}\over C}&>\left[{N_2\over N_1}
\underbrace{e^{-{E_{32}\over kT}}}_{E_{32}=E_{32}-E_{12}}-{N_3\over N_2}\right]
+\left[{N_2\over N_1}-e^{-{E_{21}\over kT}}\right]\\
{A_{32}\theta\over C}&>\left[\underbrace{N_2\over N_1}_{1}\theta(1+\delta)-
\theta\right]+\left[\underbrace{N_2\over N_1}_{1}-(1-\delta)\right]\\
&=\theta\delta+\delta=\delta(\theta+1)\\ \end{aligned}$$
Thus, the $\delta$ really helps get masing started.
\end{itemize}
Now one last thing: we'd
chosen to ignore stimulated radiative transfer between energy states 3 and 2.
In general, this process
will tend to reduce the population inversion. However, for optically thick
clouds $\tau\gg1$, photons only have a $P\sim\inv{\tau}$ probability of escaping,
so we can describe this by ``diluting* the $A_{32}$ term by $\inv{\tau}$.*
\subsubsection{The 4 state atom}
While 3 state masers are possible, most actual masers involve 4 states. Lets discuss population inversion in a 4 state maser. We will relax the condition that collisions
be responsible for the transition to the highest state. As shown in the diagram we require that there is some strong transition from state 1 to 4. This provides atoms in a state that can dexcite to our population-inverted state. Imagine that we want states 3 and 2 to be the population inverted states. To achieve this, we need a relatively strong transition from state 4 to 3 to populate state 3. Similarly we need a strong 2 to 1 transition to empty out state 2. We are thus left with lots of atoms in state 3, and almost none in 2. This provides the inversion necessary for masing. Our masing transition will be 3 to 2 (ie it will correspond to our population inverted states).
\end{document}
<\latex>