 Three-sublevel optical pumping model

Consider an original hypothetical atom with a lower angular momentum multiplet $F=1$, $m=-1,0,+1$, and one opposite-parity excited level $F^{\prime }=0,m^{\prime }=0$. Stipulate that all spontaneous emission returns to these three lower quantum states, with no other decay channel. This is an isolated four-state model, not the actual sodium or rubidium D-line structure. In particular, the physical rubidium cell above contains multiple hyperfine structure levels.

The electric-dipole selection rule permits $\Delta F=-1$ here; ${\sigma }^{+}$ absorption with $\Delta m=+1$ excites only $m=-1$. Decay from $F^{\prime }=0$ has three equal integrated branching probabilities $b_{-}=b_0=b_{+}=1/3$: rotational symmetry gives equal squared Clebsch-Gordan coefficients. The emitted photons carry the three permitted angular momentum projections; equal integrated branching does not imply equal intensity in every observation direction. Thus $m=0$ and $m=+1$ are both dark. The chosen beam cannot prepare a unique pure $m=+1$ state.

Let $R\ge 0$ be the absorption rate of an atom in $m=-1$, and let $\Gamma$ be the excited-state spontaneous emission rate, both in inverse seconds. Assume $R,\gamma \ll \Gamma$, negligible stimulated emission, and no maintained quantum coherence between sublevels. Eliminating the short-lived excited population gives normalized lower populations $p_{-}+p_0+p_{+}=1$ to this approximation. At each scattering event the atom leaves $m=-1$ with probability $2/3$; the remaining $1/3$ simply returns it there.

Model depolarizing collisions and other population relaxation by resets to the unpolarized distribution at rate $\gamma \ge 0$. The rate equations areThis relaxation prescription assumes that Zeeman splitting is small compared with thermal energy, so the unpolarized distribution approximates thermal equilibrium. It is not an energy-independent relaxation law at arbitrarily low temperature.

Every transfer has a nonnegative rate: optical transfers $-\to 0$ and $-\to +$ each have rate $R/3$, and relaxation transfers between each distinct pair have rate $\gamma /3$. These define a continuous-time Markov chain. The three derivatives sum to zero, and at $p_i=0$ the corresponding derivative cannot be negative. Hence the model conserves population and preserves nonnegative populations. Optical cycles that return to $m=-1$ contribute to scattering even though they cause no net population transfer.