Maser oscillation threshold (source code)

= Maser oscillation threshold
{title2=$GN_0>\kappa$}

The following original model isolates amplification by inverted <electric dipoles>. Let $N=N_u-N_l$ be the excess number in the upper member of one equally degenerate resonant pair. Let $a$ be a slowly varying cavity amplitude with <photon> number $n=|a|^2$, and let $p$ describe the collective molecular dipole component in the amplifying phase convention. On exact <resonance>, stipulate the linear equations
$$
\dot a=-\frac{\kappa}{2}a+g_0p,
\qquad
\dot p=-\Gamma p+g_0Na,
\qquad \Gamma>0.
$$
The dipole relaxation rate is $\Gamma$, and $g_0=|d_{ul}\mathcal E_{\mathrm{vac}}|/\hbar$ is a single-molecule coupling in the chosen mode normalization. Here $d_{ul}$ is an <electric dipole moment> <matrix element>, $\mathcal E_{\mathrm{vac}}$ is the field's zero-point amplitude, and $\hbar=h/(2\pi)$. Field overlap and dipole orientation enter this coupling. The <characteristic polynomial> for growth rate $\lambda$ is
$$
(\lambda+\kappa/2)(\lambda+\Gamma)-g_0^2N=0.
$$
For $N>0$, one root becomes positive precisely when $g_0^2N>\kappa\Gamma/2$. If dipole relaxation is fast enough to eliminate $p$, then $p\approx g_0Na/\Gamma$ and
$$
\dot n=(GN-\kappa)n,
\qquad G=\frac{2g_0^2}{\Gamma}.
$$
Thus positive <population inversion> alone is insufficient: the gain rate $GN$ must exceed the passive decay rate. This reduction describes the onset near <resonance>; it does not predict output <phase> noise, actual transit dynamics or an exact historical threshold. <Nikolay Basov> and <Aleksandr Prokhorov>'s molecular-generator theory likewise treats self-excitation as competition between molecular gain and resonator-plus-load losses.<Theory of the Molecular Generator and Molecular Power Amplifier by Nicolay Gennadiyevich Basov and Aleksandr Mikhailovich Prokhorov>[p. 426, Introduction; pp. 427–428, dipole response and self-excitation]{ref}

To include replenishment, suppose upper-state <molecules> enter at rate $F$ and both levels leave with residence rate $\gamma$. A deliberately simple <ordinary differential equation> model is
$$
\dot N_u=F-\gamma N_u-GNn,
\qquad
\dot N_l=-\gamma N_l+GNn,
\qquad
\dot n=(GN-\kappa)n.
$$
One downward transition removes two units of the difference $N_u-N_l$, so $\dot N=F-\gamma N-2GNn$. Below onset the unsaturated inversion is $N_0=F/\gamma$, giving
$$
GN_0>\kappa,
\qquad
F>F_{\mathrm{th}}=\frac{\gamma\kappa}{G}.
$$
For fixed mode volume $V$, the equivalent <number density> condition is $N_0/V>\kappa/(GV)$. Raising selected beam intensity raises gain; increasing <quality factor> lowers the loss rate at fixed <frequency>. These statements hold with $G$, residence time and mode overlap fixed.

Above onset, the nonzero <steady state> has
$$
N_{\mathrm{ss}}=\frac{\kappa}{G},
\qquad
n_{\mathrm{ss}}=\frac{F-\gamma\kappa/G}{2\kappa},
\qquad
P_{\mathrm{out}}=\kappa_{\mathrm{out}}h\nu_0n_{\mathrm{ss}}.
$$
The inversion is clamped while excess supply supports output. A finite initial seed is understood: the noiseless equations leave $n=0$ exactly unchanged even above onset. Actual startup has <spontaneous emission> and field fluctuations.

\Image[/ammonia-maser-gain-load-and-energy-budget.png]
{title=<Maser oscillation threshold>, output loading and the <maser microwave energy budget>}
{description=Original illustrative model, not historical measurements. Left: gain $GN_0$ rises with selected beam supply, crossing two chosen total loss rates. Right: with fixed supply, output initially grows with output coupling, then falls to zero when that coupling makes the total loss too large. The lower panel accounts for incoming excitation, unused excitation leaving in the beam, wall dissipation and useful output for the worked example.}