Magnetic Resonance Linewidth of Alkali-Metal Vapor in Unresolved Zeeman Resonance Regime

Feng Tang Beijing Computational Science Research Center, Beijing 100193, PR China    Nan Zhao [email protected] Beijing Computational Science Research Center, Beijing 100193, PR China
(July 23, 2024)
Abstract

The study of magnetic resonance linewidth is crucial in magnetic resonance physics and its applications. Previous studies focused on the linewidth of alkali metal atoms within the spin-exchange relaxation-free regime near zero magnetic field and in strong magnetic fields where Zeeman resonances are well resolved due to the quadratic Zeeman effect. However, the linewidth in the unresolved Zeeman resonance regime, which is prevalent in various magnetometer and comagnetometer applications, is not well understood. To address this, we developed a theoretical framework based on the master equation for alkali metal atoms and solved it under the rotating wave approximation and weak driving conditions. Our numerical calculations and analytical expressions reveal that the light-narrowing effect occurs only when the ratio of the spin exchange rate to the spin destruction rate exceeds a critical value. Additionally, we show that the linewidth in the unresolved Zeeman resonance regime is significantly influenced by the mutual coupling of quantum coherence between different Zeeman sublevels. These findings provide a theoretical tool for understanding spin relaxation in alkali-metal atoms and optimizing the performance of atomic magnetometers and comagnetometers operating in this regime.

I Introduction

The resonant linewidth is a crucial physical parameter in both fundamental and applied research on magnetic resonance. In fundamental studies, the resonant linewidth reveals the interaction between spins and their microscopic environments. This linewidth can be utilized to derive various physical insights about the surroundings of the spins, such as the intensity of magnetic field fluctuations and correlation times. In practical applications, including magnetic field sensing and magnetic resonance imaging, the magnetic resonance linewidth defines essential performance metrics like magnetic detection sensitivity and resolution. Hence, understanding the mechanisms underlying the resonant linewidth is highly important.

Magnetic resonance of spin-polarized alkali metal atoms has been extensively utilized in high-precision measurement domains, such as high-sensitivity magnetometers[1], inertial sensors[2, 3], and fundamental physical model testing[4]. The resonant linewidth of alkali-metal atoms primarily depends on the spin-related collision processes involving alkali metal atoms and other atoms or molecules. These processes include spin depolarization of alkali metal atoms through spin destruction events and spin exchange processes that lead to spin phase relaxation. The spin exchange collision process among alkali metal atoms is crucial in determining the magnetic resonance linewidth. For high density alkali metal vapor in low magnetic fields, when the spin-exchange rate ΓSEsubscriptΓSE\Gamma_{\mathrm{SE}}roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT considerably exceeds the Larmor precession frequency of spin rotation ΩLsubscriptΩL\Omega_{\mathrm{L}}roman_Ω start_POSTSUBSCRIPT roman_L end_POSTSUBSCRIPT (i.e., ΓSEΩLmuch-greater-thansubscriptΓSEsubscriptΩL\Gamma_{\mathrm{SE}}\gg\Omega_{\mathrm{L}}roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT ≫ roman_Ω start_POSTSUBSCRIPT roman_L end_POSTSUBSCRIPT), the system enters a spin-exchange relaxation-free regime (SERF) with a notable reduction in the resonant linewidth [5]. Increasing the magnetic field violates the SERF conditions, causing the spin resonance linewidth of high density alkali metal vapor to be primarily determined by the spin exchange rate. In this situation, the resonant linewidth also depends on the level of atomic spin polarization, and increasing atomic spin polarization via optical pumping results in the light-narrowing effect [6, 7]. A quantitative understanding of the various collision-induced physical effects is essential for the use of alkali metal vapor in precision measurements.

The master equation is used to describe various microscopic physical processes that influence the spin resonance linewidth [8]. The microscopic processes can be quantitatively characterized by their collision cross-sections, and magnetic resonance linewidth in obtained by solving the master equation quantitatively under different conditions. In Ref. [8], Appelt et al. extensively analysed the spin relaxation mechanisms of alkali metal atoms and calculated the linewidth within the resolved Zeeman resonance (RZR) regime, where the frequency degeneracy of coherences between adjacent Zeeman levels is lifted by the quadratic Zeeman effect in strong fields [9, 10]. Unfortunately, because of the complexity of alkali metal energy levels and the nonlinear nature of the spin-exchange process, calculating the magnetic resonance linewidth under more general conditions remains a challenging task.

Magnetic resonance linewidths show varying characteristics depending on the strength of the magnetic field and the spin polarization. Extensive studies have been conducted on resonance linewidths under both weak and strong magnetic field conditions. In the SERF regime at near-zero magnetic fields, the resonance linewidth is inversely proportional to the spin exchange rate[5]. In strong magnetic fields, the behavior of the magnetic resonance linewidth in the RZR regime has been explored in Refs.[8, 7, 11]. Conversely, in the intermediate magnetic field region, the interaction between quantum coherences of different Zeeman sublevels plays a significant role in influencing the magnetic resonance linewidth. A precise description of the linewidth in the unresolved Zeeman resonance (UZR) regime remains unavailable. The advance in various atomic magnetometers operating in the UZR regime necessitates a more thorough examination of the linewidth.

Based on the master equation in Ref. [8], we investigated the linewidth of the atom resonance under various conditions. We analyze the mathematical structure of the evolution matrix of the master equation, simplifying it under the rotating wave approximation and weak driving scenarios. This simplification not only makes the master equation more tractable, but also offers a clearer understanding of the physical phenomena under differing conditions. Our findings indicate that, in the UZR regime, the interplay of quantum coherence among various Zeeman sublevels has a significant impact on the magnetic resonance linewidth. We derived analytical expressions for the linewidth in both high- and low-spin polarization limits, which will benefit the optimization of parameters for atomic magnetometers and gyroscopes operating in the UZR regime.

II Theoretical Treatment

Refer to caption
Figure 1: Heat map plot of the matrix [𝒜SD]delimited-[]subscript𝒜SD\left[\mathcal{A}_{\rm SD}\right][ caligraphic_A start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT ]. The bases are denoted by |F,mF;F,mF)=|mF;mF)|F,m_{F};F^{\prime},m^{\prime}_{F^{\prime}})=\left|m_{F};m^{\prime}_{F^{\prime% }}\right)| italic_F , italic_m start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ; italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ) = | italic_m start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ; italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ) [see Eq. (7)], and are grouped according to their frequencies in the rotating frame. The horizontal and vertical dashed lines separate the matrix blocks according to their frequencies in the rotating frame. The block marked by the solid lines is [𝒜SD]1,1subscriptdelimited-[]subscript𝒜SD11\left[\mathcal{A}_{\rm SD}\right]_{1,1}[ caligraphic_A start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT, the diagonal block of the first-order Zeeman coherences. The block marked by the dotted lines is [𝒜SD]1,1subscriptdelimited-[]subscript𝒜SD11\left[\mathcal{A}_{\rm SD}\right]_{1,-1}[ caligraphic_A start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , - 1 end_POSTSUBSCRIPT, which is a double frequency block in the rotating frame and is neglected with the RWA.

II.1 Master equation in Liouville space

We study alkali metal spins in a static magnetic field B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT in the z^^𝑧\hat{z}over^ start_ARG italic_z end_ARG direction and a radio-frequency (RF) driving field in the x^^𝑥\hat{x}over^ start_ARG italic_x end_ARG direction. An circularly polarized laser beam travels along the z^^𝑧\hat{z}over^ start_ARG italic_z end_ARG axis to pump the alkali metal atoms. The primary incoherent interactions that affect the spin dynamics of alkali metals are spin destruction collisions between alkali metal atoms and other atoms and molecules and spin exchange collisions among alkali metal atoms. The master equation governing the density matrix ρ𝜌\rhoitalic_ρ of the alkali-metal atom is [8]

dρdt=𝑑𝜌𝑑𝑡absent\displaystyle\frac{d\rho}{dt}=divide start_ARG italic_d italic_ρ end_ARG start_ARG italic_d italic_t end_ARG = i[H,ρ]+ΓSE[φ(1+4𝐒𝐒)ρ]𝑖𝐻𝜌subscriptΓSEdelimited-[]𝜑14delimited-⟨⟩𝐒𝐒𝜌\displaystyle-i[H,\rho]+\Gamma_{\text{SE}}[\varphi(1+4\langle\mathbf{S}\rangle% \cdot\mathbf{S})-\rho]- italic_i [ italic_H , italic_ρ ] + roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT [ italic_φ ( 1 + 4 ⟨ bold_S ⟩ ⋅ bold_S ) - italic_ρ ]
+ΓSD(φρ)+ROP[φ(1+2Sz)ρ].subscriptΓSD𝜑𝜌subscript𝑅OPdelimited-[]𝜑12subscript𝑆𝑧𝜌\displaystyle+\Gamma_{\text{SD}}(\varphi-\rho)+R_{\text{OP}}[\varphi(1+2S_{z})% -\rho].+ roman_Γ start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT ( italic_φ - italic_ρ ) + italic_R start_POSTSUBSCRIPT OP end_POSTSUBSCRIPT [ italic_φ ( 1 + 2 italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ) - italic_ρ ] . (1)

The coherent evolution is determined by the Hamiltonian

H=H0+Hdrv,𝐻subscript𝐻0subscript𝐻drvH=H_{0}+H_{\text{drv}},italic_H = italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_H start_POSTSUBSCRIPT drv end_POSTSUBSCRIPT , (2)

where

H0subscript𝐻0\displaystyle H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT =Ω0Sz,+Ωhf𝐈𝐒absentsubscriptΩ0subscript𝑆𝑧subscriptΩhf𝐈𝐒\displaystyle=\Omega_{0}S_{z},+\Omega_{\rm hf}\mathbf{I}\cdot\mathbf{S}= roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , + roman_Ω start_POSTSUBSCRIPT roman_hf end_POSTSUBSCRIPT bold_I ⋅ bold_S (3)
Hdrvsubscript𝐻drv\displaystyle H_{\text{drv}}italic_H start_POSTSUBSCRIPT drv end_POSTSUBSCRIPT =ΩRSxcos(ωt+ϕ0),absentsubscriptΩ𝑅subscript𝑆𝑥𝜔𝑡subscriptitalic-ϕ0\displaystyle=\Omega_{R}S_{x}\cos(\omega t+\phi_{0}),= roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT roman_cos ( italic_ω italic_t + italic_ϕ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , (4)

represent the atomic spin Hamiltonian within the magnetic field B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, along with the RF driving field. In Eq. (3), 𝐒𝐒\mathbf{S}bold_S and 𝐈𝐈\mathbf{I}bold_I correspond to the electronic spin and nuclear spin operators of the alkali metal atom. Here, Ω0=|γeB0|>0subscriptΩ0subscript𝛾𝑒subscript𝐵00\Omega_{0}=|\gamma_{e}B_{0}|>0roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = | italic_γ start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | > 0 denotes the Larmor frequency (γesubscript𝛾𝑒\gamma_{e}italic_γ start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT being the gyromagnetic ratio), and ΩhfsubscriptΩhf\Omega_{\rm hf}roman_Ω start_POSTSUBSCRIPT roman_hf end_POSTSUBSCRIPT is the hyperfine interaction constant. Throughout this paper, it is assumed that the magnetic field B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is sufficiently weak so that Ω0Ωhfmuch-less-thansubscriptΩ0subscriptΩhf\Omega_{0}\ll\Omega_{\rm hf}roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ≪ roman_Ω start_POSTSUBSCRIPT roman_hf end_POSTSUBSCRIPT holds true. In Eq. (4), ΩRsubscriptΩ𝑅\Omega_{R}roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT, ω𝜔\omegaitalic_ω and ϕ0subscriptitalic-ϕ0\phi_{0}italic_ϕ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT stand for the Rabi frequency, the frequency and the phase of the RF driving field, respectively. For simplicity, ϕ0=0subscriptitalic-ϕ00\phi_{0}=0italic_ϕ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0 is assumed in the following text.

The incoherent part of the master Eq. (II.1) comprises spin exchange (SE) collisions, spin destruction (SD) collisions, and optical pumping (OP) mechanisms. The second term in Eq. (II.1) characterizes the spin exchange collisions among alkali metal spins at a rate of ΓSEsubscriptΓSE\Gamma_{\text{SE}}roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT, where 𝐒=tr[ρ𝐒]delimited-⟨⟩𝐒trdelimited-[]𝜌𝐒\langle\mathbf{S}\rangle=\text{tr}[\rho\mathbf{S}]⟨ bold_S ⟩ = tr [ italic_ρ bold_S ] represents the spin expectation value. The third term results from spin destruction collisions between alkali metal atoms and other atoms or molecules, with a total spin destruction rate denoted as ΓSDsubscriptΓSD\Gamma_{\text{SD}}roman_Γ start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT. In Eq. (II.1), we have assumed that the optical pumping occurs in the high broadening limit of the optical transition[12, 13], and the pumping effect is described by the final term of Eq. (II.1) with a pumping rate ROPsubscript𝑅OPR_{\text{OP}}italic_R start_POSTSUBSCRIPT OP end_POSTSUBSCRIPT. The state

φ=14ρ+𝐒ρ𝐒𝜑14𝜌𝐒𝜌𝐒\varphi=\frac{1}{4}\rho+\mathbf{S}\cdot\rho\mathbf{S}italic_φ = divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_ρ + bold_S ⋅ italic_ρ bold_S (5)

is a purely nuclear spin operator of the alkali-metal atom, without electronic spin polarization [12].

In Liouville space, the density matrix ρ𝜌\rhoitalic_ρ is mapped to a state vector |ρ)|\rho)| italic_ρ ). Thee state vector is expanded as

|ρ)=m,m,F,F|mF;mF)(mF;mF|ρ),|\rho)=\sum_{m,m^{\prime},F,F^{\prime}}|m_{F};m^{\prime}_{F^{\prime}})(m_{F};m% ^{\prime}_{F^{\prime}}|\rho),| italic_ρ ) = ∑ start_POSTSUBSCRIPT italic_m , italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_F , italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_m start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ; italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ) ( italic_m start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ; italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ρ ) , (6)

where

|mF;mF)=|F,mF,m||m_{F};m^{\prime}_{F^{\prime}})=|F,m\rangle\langle F^{\prime},m^{\prime}|| italic_m start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ; italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ) = | italic_F , italic_m ⟩ ⟨ italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT | (7)

is the basis vector of the Liouville space, and |F,mket𝐹𝑚|F,m\rangle| italic_F , italic_m ⟩ is the eigen vector of the alkali-metal Hamiltonian (3), with the total angular momentum quantum number F=I+1/2a𝐹𝐼12𝑎F=I+1/2\equiv aitalic_F = italic_I + 1 / 2 ≡ italic_a or F=I1/2b𝐹𝐼12𝑏F=I-1/2\equiv bitalic_F = italic_I - 1 / 2 ≡ italic_b, and the magnetic quantum number m=F,,F𝑚𝐹𝐹m=-F,\dots,Fitalic_m = - italic_F , … , italic_F. For the state vectors constrained in the subspace of population and Zeeman coherences with F=F𝐹superscript𝐹F=F^{\prime}italic_F = italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT, we define

|F,m¯)Δm=|mF;mF),|F,\bar{m})_{\Delta m}=|m_{F};m^{\prime}_{F}),| italic_F , over¯ start_ARG italic_m end_ARG ) start_POSTSUBSCRIPT roman_Δ italic_m end_POSTSUBSCRIPT = | italic_m start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ; italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ) , (8)

where m¯=(m+m)/2¯𝑚𝑚superscript𝑚2\bar{m}=(m+m^{\prime})/2over¯ start_ARG italic_m end_ARG = ( italic_m + italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) / 2 and Δm=mmΔ𝑚𝑚superscript𝑚\Delta m=m-m^{\prime}roman_Δ italic_m = italic_m - italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT [8].

The master Eq. (II.1) in Liouville space becomes

ddt|ρ)=\displaystyle\frac{d}{dt}|\rho)=divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG | italic_ρ ) = i|ρ)(ΓSE+ΓSD+ROP)𝒜SD|ρ)\displaystyle-i\mathcal{H}|\rho)-(\Gamma_{\text{SE}}+\Gamma_{\text{SD}}+R_{% \text{OP}})\mathcal{A}_{\text{SD}}|\rho)- italic_i caligraphic_H | italic_ρ ) - ( roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT + roman_Γ start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT + italic_R start_POSTSUBSCRIPT OP end_POSTSUBSCRIPT ) caligraphic_A start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT | italic_ρ )
+(ΓSE(𝐒|ρ)+12Ropz^)𝓐SE|ρ)𝒢|ρ),\displaystyle+\left(\Gamma_{\text{SE}}(\mathbf{S}|\rho)+\frac{1}{2}R_{\rm op}% \hat{z}\right)\cdot\bm{\mathcal{A}}_{\text{SE}}|\rho)\equiv-\mathcal{G}|\rho),+ ( roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT ( bold_S | italic_ρ ) + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT over^ start_ARG italic_z end_ARG ) ⋅ bold_caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT | italic_ρ ) ≡ - caligraphic_G | italic_ρ ) , (9)

where 𝐒=(𝐒|ρ)delimited-⟨⟩𝐒conditional𝐒𝜌\langle\mathbf{S}\rangle=(\mathbf{S}|\rho)⟨ bold_S ⟩ = ( bold_S | italic_ρ ) is the spin expectation value expressed in Liouville space, and \mathcal{H}caligraphic_H, 𝒜SDsubscript𝒜SD\mathcal{A}_{\rm SD}caligraphic_A start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT and 𝓐SEsubscript𝓐SE\bm{\mathcal{A}}_{\rm SE}bold_caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT are the superoperator in the Liouville space, corresponding to the coherent evolution, spin destruction and spin exchange processes, respectively. The superoperators are defined as [13]

\displaystyle\mathcal{H}caligraphic_H =0+drv=H0©+Hdrv©,absentsubscript0subscriptdrvsubscriptsuperscript𝐻©0superscriptsubscript𝐻drv©\displaystyle=\mathcal{H}_{0}+\mathcal{H}_{\rm drv}=H^{\copyright}_{0}+H_{\rm drv% }^{\copyright},= caligraphic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + caligraphic_H start_POSTSUBSCRIPT roman_drv end_POSTSUBSCRIPT = italic_H start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_H start_POSTSUBSCRIPT roman_drv end_POSTSUBSCRIPT start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT , (10)
𝒜SDsubscript𝒜SD\displaystyle\mathcal{A}_{\text{SD}}caligraphic_A start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT =12𝐒©𝐒©,absent12superscript𝐒©superscript𝐒©\displaystyle=\frac{1}{2}\mathbf{S}^{\copyright}\cdot\mathbf{S}^{\copyright},= divide start_ARG 1 end_ARG start_ARG 2 end_ARG bold_S start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT ⋅ bold_S start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT , (11)
𝓐SEsubscript𝓐SE\displaystyle\bm{\mathcal{A}}_{\text{SE}}bold_caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT =𝐒+𝐒2i𝐒×𝐒,absentsuperscript𝐒superscript𝐒2𝑖superscript𝐒superscript𝐒\displaystyle=\mathbf{S}^{\flat}+\mathbf{S}^{\sharp}-2i\mathbf{S}^{\flat}% \times\mathbf{S}^{\sharp},= bold_S start_POSTSUPERSCRIPT ♭ end_POSTSUPERSCRIPT + bold_S start_POSTSUPERSCRIPT ♯ end_POSTSUPERSCRIPT - 2 italic_i bold_S start_POSTSUPERSCRIPT ♭ end_POSTSUPERSCRIPT × bold_S start_POSTSUPERSCRIPT ♯ end_POSTSUPERSCRIPT , (12)

where Xsuperscript𝑋X^{\flat}italic_X start_POSTSUPERSCRIPT ♭ end_POSTSUPERSCRIPT and Xsuperscript𝑋X^{\sharp}italic_X start_POSTSUPERSCRIPT ♯ end_POSTSUPERSCRIPT are the left- and right-translation superoperators, corresponding to the operator X𝑋Xitalic_X in Hilbert space, and X©=XXsuperscript𝑋©superscript𝑋superscript𝑋X^{\copyright}=X^{\flat}-X^{\sharp}italic_X start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT = italic_X start_POSTSUPERSCRIPT ♭ end_POSTSUPERSCRIPT - italic_X start_POSTSUPERSCRIPT ♯ end_POSTSUPERSCRIPT is the commutator superoperator [13].

II.2 Rotating wave approximation

We consider the near-resonant RF driving field with frequency ωΩ0𝜔subscriptΩ0\omega\approx\Omega_{0}italic_ω ≈ roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The driving field only induces coherences between the Zeeman sublevels within the F=a𝐹𝑎F=aitalic_F = italic_a or b𝑏bitalic_b subspaces. In the absence of microwave excitations, the hyperfine coherences between levels F=a𝐹𝑎F=aitalic_F = italic_a and F=b𝐹𝑏F=bitalic_F = italic_b (e.g., |ma;mb)|m_{a};m^{\prime}_{b})| italic_m start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ; italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT )) will decay within the time scale of the hyperfine relaxation time, typically in the order of mssimilar-toabsentms\sim{\rm ms}∼ roman_ms or even shorter. In this case, we can neglect the hyperfine coherences and focus on the dynamical evolution of the population and Zeeman coherences. To this end, we define the projection operators

𝒫a(hf)subscriptsuperscript𝒫hf𝑎\displaystyle\mathcal{P}^{(\text{hf})}_{a}caligraphic_P start_POSTSUPERSCRIPT ( hf ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT =m,n|ma;na)(ma;na|,\displaystyle=\sum_{m,n}|m_{a};n_{a})(m_{a};n_{a}|,= ∑ start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT | italic_m start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ; italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) ( italic_m start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ; italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT | , (13)
𝒫b(hf)subscriptsuperscript𝒫hf𝑏\displaystyle\mathcal{P}^{(\text{hf})}_{b}caligraphic_P start_POSTSUPERSCRIPT ( hf ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT =m,n|mb;nb)(mb;nb|,\displaystyle=\sum_{m,n}|m_{b};n_{b})(m_{b};n_{b}|,= ∑ start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT | italic_m start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ; italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) ( italic_m start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ; italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT | , (14)
𝒫Z(hf)subscriptsuperscript𝒫hfZ\displaystyle\mathcal{P}^{(\text{hf})}_{\text{Z}}caligraphic_P start_POSTSUPERSCRIPT ( hf ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT =𝒫a(hf)+𝒫b(hf),absentsubscriptsuperscript𝒫hf𝑎subscriptsuperscript𝒫hf𝑏\displaystyle=\mathcal{P}^{(\text{hf})}_{a}+\mathcal{P}^{(\text{hf})}_{b},= caligraphic_P start_POSTSUPERSCRIPT ( hf ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT + caligraphic_P start_POSTSUPERSCRIPT ( hf ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , (15)

and the state vector projected to the population and Zeeman coherence subspace is

|ρZ)=𝒫Z(hf)|ρ).|\rho_{\text{Z}})=\mathcal{P}_{\text{Z}}^{(\text{hf})}|\rho).| italic_ρ start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT ) = caligraphic_P start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( hf ) end_POSTSUPERSCRIPT | italic_ρ ) . (16)

The evolution of |ρZ)|\rho_{\text{Z}})| italic_ρ start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT ) is governed by

ddt|ρZ)=𝒫Z(hf)𝒢𝒫Z(hf)|ρZ)𝒢Z|ρZ)\displaystyle\frac{d}{dt}|\rho_{\text{Z}})=-\mathcal{P}_{\text{Z}}^{(\text{hf}% )}\mathcal{G}\mathcal{P}_{\text{Z}}^{(\text{hf})}|\rho_{\text{Z}})\equiv-% \mathcal{G}_{{\text{Z}}}|\rho_{\text{Z}})divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG | italic_ρ start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT ) = - caligraphic_P start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( hf ) end_POSTSUPERSCRIPT caligraphic_G caligraphic_P start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( hf ) end_POSTSUPERSCRIPT | italic_ρ start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT ) ≡ - caligraphic_G start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT | italic_ρ start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT ) (17)

where

𝒢Z=i[]+Γtot[𝒜SD](ΓSE(𝐒|ρZ)+12ROPz^)[𝓐SE],subscript𝒢Z𝑖delimited-[]subscriptΓtotdelimited-[]subscript𝒜SDsubscriptΓSEconditional𝐒subscript𝜌𝑍12subscript𝑅OP^𝑧delimited-[]subscript𝓐SE\displaystyle\mathcal{G}_{{\text{Z}}}=i\left[\mathcal{H}\right]+\Gamma_{\text{% tot}}\left[\mathcal{A}_{\text{SD}}\right]-\left(\Gamma_{\text{SE}}(\mathbf{S}|% \rho_{Z})+\frac{1}{2}R_{\text{OP}}\hat{z}\right)\cdot\left[\bm{\mathcal{A}}_{% \text{SE}}\right],caligraphic_G start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT = italic_i [ caligraphic_H ] + roman_Γ start_POSTSUBSCRIPT tot end_POSTSUBSCRIPT [ caligraphic_A start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT ] - ( roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT ( bold_S | italic_ρ start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT ) + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_R start_POSTSUBSCRIPT OP end_POSTSUBSCRIPT over^ start_ARG italic_z end_ARG ) ⋅ [ bold_caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT ] , (18)

where Γtot=ΓSE+ΓSD+ROPsubscriptΓtotsubscriptΓSEsubscriptΓSDsubscript𝑅OP\Gamma_{\text{tot}}=\Gamma_{\text{SE}}+\Gamma_{\text{SD}}+R_{\text{OP}}roman_Γ start_POSTSUBSCRIPT tot end_POSTSUBSCRIPT = roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT + roman_Γ start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT + italic_R start_POSTSUBSCRIPT OP end_POSTSUBSCRIPT, and [𝒳]=𝒫Z(hf)𝒳𝒫Z(hf)delimited-[]𝒳subscriptsuperscript𝒫hfZ𝒳subscriptsuperscript𝒫hfZ\left[\mathcal{X}\right]=\mathcal{P}^{(\text{hf})}_{\text{Z}}\mathcal{X}% \mathcal{P}^{(\text{hf})}_{\text{Z}}[ caligraphic_X ] = caligraphic_P start_POSTSUPERSCRIPT ( hf ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT caligraphic_X caligraphic_P start_POSTSUPERSCRIPT ( hf ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT represents the projection of the superoperator 𝒳𝒳\mathcal{X}caligraphic_X to the subspace of population and Zeeman coherences. As an example, Fig. 1 shows the matrix of the superoperator 𝒜SDsubscript𝒜SD\mathcal{A}_{\rm SD}caligraphic_A start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT on the basis of the population and Zeeman coherences.

To further simplify the master equation, we define the rotating frame with the frequency ω𝜔\omegaitalic_ω and make the rotating wave approximation (RWA). We introduce the k𝑘kitalic_kth order Zeeman coherence projector.

𝒫k=m¯|a,m¯)kk(a,m¯|+m¯|b,m¯)kk(b,m¯|,\displaystyle\mathcal{P}_{k}=\sum_{\bar{m}}|a,\bar{m})_{k}~{}_{k}(a,\bar{m}|+% \sum_{\bar{m}}|b,\bar{m})_{-k}~{}_{-k}(b,\bar{m}|,caligraphic_P start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = ∑ start_POSTSUBSCRIPT over¯ start_ARG italic_m end_ARG end_POSTSUBSCRIPT | italic_a , over¯ start_ARG italic_m end_ARG ) start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_FLOATSUBSCRIPT italic_k end_FLOATSUBSCRIPT ( italic_a , over¯ start_ARG italic_m end_ARG | + ∑ start_POSTSUBSCRIPT over¯ start_ARG italic_m end_ARG end_POSTSUBSCRIPT | italic_b , over¯ start_ARG italic_m end_ARG ) start_POSTSUBSCRIPT - italic_k end_POSTSUBSCRIPT start_FLOATSUBSCRIPT - italic_k end_FLOATSUBSCRIPT ( italic_b , over¯ start_ARG italic_m end_ARG | , (19)

and the free evolution superoperator

𝒢0=ikkω𝒫k.subscript𝒢0𝑖subscript𝑘𝑘𝜔subscript𝒫𝑘\displaystyle\mathcal{G}_{0}=i\sum_{k}k\omega\mathcal{P}_{k}.caligraphic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_i ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT italic_k italic_ω caligraphic_P start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT . (20)

With the rotation transform generated by 𝒢0subscript𝒢0\mathcal{G}_{0}caligraphic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, i.e.,

|ρZ)\displaystyle|\rho_{\text{Z}})| italic_ρ start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT ) =exp(𝒢0t)|ρ~Z)𝒰(t)|ρ~),\displaystyle=\exp\left(-\mathcal{G}_{0}t\right)|\tilde{\rho}_{\text{Z}})% \equiv\mathcal{U}(t)|\tilde{\rho}),= roman_exp ( - caligraphic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_t ) | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT ) ≡ caligraphic_U ( italic_t ) | over~ start_ARG italic_ρ end_ARG ) , (21)

the master equation of the slow-varying state vector |ρ~)|\tilde{\rho})| over~ start_ARG italic_ρ end_ARG ) in the rotating frame is

ddt|ρ~)=(𝒰𝒢Z𝒰𝒢0)|ρ~)𝒢rot(t)|ρ~).\displaystyle\frac{d}{dt}|\tilde{\rho})=-\left(\mathcal{U}^{\dagger}\mathcal{G% }_{\text{Z}}\mathcal{U}-\mathcal{G}_{0}\right)|\tilde{\rho})\equiv-\mathcal{G}% _{\text{rot}}(t)|\tilde{\rho}).divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG | over~ start_ARG italic_ρ end_ARG ) = - ( caligraphic_U start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT caligraphic_G start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT caligraphic_U - caligraphic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) | over~ start_ARG italic_ρ end_ARG ) ≡ - caligraphic_G start_POSTSUBSCRIPT rot end_POSTSUBSCRIPT ( italic_t ) | over~ start_ARG italic_ρ end_ARG ) . (22)

The evolution superoperator 𝒢rot(t)subscript𝒢rot𝑡\mathcal{G}_{\text{rot}}(t)caligraphic_G start_POSTSUBSCRIPT rot end_POSTSUBSCRIPT ( italic_t ) is

𝒢rot(t)=ik[Δ]k,k+𝒢drv(t)+𝒢relax(t),subscript𝒢rot𝑡𝑖subscript𝑘subscriptdelimited-[]Δ𝑘𝑘subscript𝒢drv𝑡subscript𝒢relax𝑡\mathcal{G}_{\text{rot}}(t)=i\sum_{k}\left[\Delta\right]_{k,k}+\mathcal{G}_{% \text{drv}}(t)+\mathcal{G}_{\text{relax}}(t),caligraphic_G start_POSTSUBSCRIPT rot end_POSTSUBSCRIPT ( italic_t ) = italic_i ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT [ roman_Δ ] start_POSTSUBSCRIPT italic_k , italic_k end_POSTSUBSCRIPT + caligraphic_G start_POSTSUBSCRIPT drv end_POSTSUBSCRIPT ( italic_t ) + caligraphic_G start_POSTSUBSCRIPT relax end_POSTSUBSCRIPT ( italic_t ) , (23)

where [Δ]k,k𝒫k0𝒫kkωksubscriptdelimited-[]Δ𝑘𝑘subscript𝒫𝑘subscript0subscript𝒫𝑘𝑘𝜔subscript𝑘\left[\Delta\right]_{k,k}\equiv\mathcal{P}_{k}\mathcal{H}_{0}\mathcal{P}_{k}-k% \omega\mathcal{I}_{k}[ roman_Δ ] start_POSTSUBSCRIPT italic_k , italic_k end_POSTSUBSCRIPT ≡ caligraphic_P start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT caligraphic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT caligraphic_P start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT - italic_k italic_ω caligraphic_I start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT is a diagonal matrix representing the frequency detuning of the k𝑘kitalic_kth order coherences (ksubscript𝑘\mathcal{I}_{k}caligraphic_I start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT is an identity matrix). The driving superoperator 𝒢drv(t)subscript𝒢drv𝑡\mathcal{G}_{\text{drv}}(t)caligraphic_G start_POSTSUBSCRIPT drv end_POSTSUBSCRIPT ( italic_t ) in the rotating frame is

𝒢drv(t)=iΩRcos(ωt)p,qei(pq)ωt[Sx©]p,q,subscript𝒢drv𝑡𝑖subscriptΩ𝑅𝜔𝑡subscript𝑝𝑞superscript𝑒𝑖𝑝𝑞𝜔𝑡subscriptdelimited-[]superscriptsubscript𝑆𝑥©𝑝𝑞\mathcal{G}_{\text{drv}}(t)=i\Omega_{R}\cos(\omega t)\sum_{p,q}e^{i(p-q)\omega t% }\left[S_{x}^{\copyright}\right]_{p,q},caligraphic_G start_POSTSUBSCRIPT drv end_POSTSUBSCRIPT ( italic_t ) = italic_i roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT roman_cos ( italic_ω italic_t ) ∑ start_POSTSUBSCRIPT italic_p , italic_q end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_p - italic_q ) italic_ω italic_t end_POSTSUPERSCRIPT [ italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT italic_p , italic_q end_POSTSUBSCRIPT , (24)

and the relaxation superoperator 𝒢relax(t)subscript𝒢relax𝑡\mathcal{G}_{\text{relax}}(t)caligraphic_G start_POSTSUBSCRIPT relax end_POSTSUBSCRIPT ( italic_t ) is

𝒢relax(t)=Γtotp,qei(pq)ωt[𝒜SD]p,qsubscript𝒢relax𝑡subscriptΓtotsubscript𝑝𝑞superscript𝑒𝑖𝑝𝑞𝜔𝑡subscriptdelimited-[]subscript𝒜SD𝑝𝑞\displaystyle\mathcal{G}_{\text{relax}}(t)=\Gamma_{\text{tot}}\sum_{p,q}e^{i(p% -q)\omega t}\left[\mathcal{A}_{\text{SD}}\right]_{p,q}caligraphic_G start_POSTSUBSCRIPT relax end_POSTSUBSCRIPT ( italic_t ) = roman_Γ start_POSTSUBSCRIPT tot end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_p , italic_q end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_p - italic_q ) italic_ω italic_t end_POSTSUPERSCRIPT [ caligraphic_A start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT italic_p , italic_q end_POSTSUBSCRIPT
(ΓSEk𝐒~(k)eikωt+12ROPz^)p,qei(pq)ωt[𝓐SE]p,q.subscriptΓSEsubscript𝑘superscript~𝐒𝑘superscript𝑒𝑖𝑘𝜔𝑡12subscript𝑅OP^𝑧subscript𝑝𝑞superscript𝑒𝑖𝑝𝑞𝜔𝑡subscriptdelimited-[]subscript𝓐SE𝑝𝑞\displaystyle-\left(\Gamma_{\text{SE}}\sum_{k}\tilde{\mathbf{S}}^{(k)}e^{-ik% \omega t}+\frac{1}{2}R_{\text{OP}}\hat{z}\right)\cdot\sum_{p,q}e^{i(p-q)\omega t% }\left[\bm{\mathcal{A}}_{\text{SE}}\right]_{p,q}.- ( roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT over~ start_ARG bold_S end_ARG start_POSTSUPERSCRIPT ( italic_k ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_k italic_ω italic_t end_POSTSUPERSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_R start_POSTSUBSCRIPT OP end_POSTSUBSCRIPT over^ start_ARG italic_z end_ARG ) ⋅ ∑ start_POSTSUBSCRIPT italic_p , italic_q end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_p - italic_q ) italic_ω italic_t end_POSTSUPERSCRIPT [ bold_caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT italic_p , italic_q end_POSTSUBSCRIPT . (25)

In Eqs. (23) - (II.2), [𝒳]p,q=𝒫p[𝒳]𝒫qsubscriptdelimited-[]𝒳𝑝𝑞subscript𝒫𝑝delimited-[]𝒳subscript𝒫𝑞\left[\mathcal{X}\right]_{p,q}=\mathcal{P}_{p}\left[\mathcal{X}\right]\mathcal% {P}_{q}[ caligraphic_X ] start_POSTSUBSCRIPT italic_p , italic_q end_POSTSUBSCRIPT = caligraphic_P start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT [ caligraphic_X ] caligraphic_P start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT is the block matrix of the p𝑝pitalic_pth row and the q𝑞qitalic_qth column, corresponding to a time-dependent factor exp(i(pq)ωt)𝑖𝑝𝑞𝜔𝑡\exp(i(p-q)\omega t)roman_exp ( italic_i ( italic_p - italic_q ) italic_ω italic_t ) in the rotating frame. The spin expectation value in the rotating frame is

𝐒delimited-⟨⟩𝐒\displaystyle\langle\mathbf{S}\rangle⟨ bold_S ⟩ =(𝐒|ρZ)=(𝐒|e𝒢0t|ρ~)absentconditional𝐒subscript𝜌Z𝐒superscript𝑒subscript𝒢0𝑡~𝜌\displaystyle=(\mathbf{S}|\rho_{\text{Z}})=(\mathbf{S}|e^{-\mathcal{G}_{0}t}|% \tilde{\rho})= ( bold_S | italic_ρ start_POSTSUBSCRIPT Z end_POSTSUBSCRIPT ) = ( bold_S | italic_e start_POSTSUPERSCRIPT - caligraphic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_t end_POSTSUPERSCRIPT | over~ start_ARG italic_ρ end_ARG )
=k(𝐒|𝒫k|ρ~)eikωtk𝐒~(k)eikωt,absentsubscript𝑘𝐒subscript𝒫𝑘~𝜌superscript𝑒𝑖𝑘𝜔𝑡subscript𝑘superscript~𝐒𝑘superscript𝑒𝑖𝑘𝜔𝑡\displaystyle=\sum_{k}(\mathbf{S}|\mathcal{P}_{k}|\tilde{\rho})e^{-ik\omega t}% \equiv\sum_{k}\tilde{\mathbf{S}}^{(k)}e^{-ik\omega t},= ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ( bold_S | caligraphic_P start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG ) italic_e start_POSTSUPERSCRIPT - italic_i italic_k italic_ω italic_t end_POSTSUPERSCRIPT ≡ ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT over~ start_ARG bold_S end_ARG start_POSTSUPERSCRIPT ( italic_k ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_k italic_ω italic_t end_POSTSUPERSCRIPT , (26)

where

𝐒~(k)=(𝐒|𝒫k|ρ~)=(𝐒|ρ~k)superscript~𝐒𝑘𝐒subscript𝒫𝑘~𝜌conditional𝐒subscript~𝜌𝑘\tilde{\mathbf{S}}^{(k)}=(\mathbf{S}|\mathcal{P}_{k}|\tilde{\rho})=(\mathbf{S}% |\tilde{\rho}_{k})over~ start_ARG bold_S end_ARG start_POSTSUPERSCRIPT ( italic_k ) end_POSTSUPERSCRIPT = ( bold_S | caligraphic_P start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG ) = ( bold_S | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) (27)

is the spin expectation value on the k𝑘kitalic_kth order Zeeman coherence |ρ~k)=𝒫k|ρ~)|\tilde{\rho}_{k})=\mathcal{P}_{k}|\tilde{\rho})| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) = caligraphic_P start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG ).

Keeping only the zero-frequency terms of 𝒢rot(t)subscript𝒢rot𝑡\mathcal{G}_{\text{rot}}(t)caligraphic_G start_POSTSUBSCRIPT rot end_POSTSUBSCRIPT ( italic_t ), we obtain the RWA equation

ddt|ρ~)=𝒢RWA|ρ~),\frac{d}{dt}|\tilde{\rho})=-\mathcal{G}_{\text{RWA}}|\tilde{\rho}),divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG | over~ start_ARG italic_ρ end_ARG ) = - caligraphic_G start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG ) , (28)

where the evolution superoperator 𝒢RWAsubscript𝒢RWA\mathcal{G}_{\text{RWA}}caligraphic_G start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT can be decomposed as 𝒢RWA=k𝒢RWA(k)subscript𝒢RWAsubscript𝑘subscriptsuperscript𝒢𝑘RWA\mathcal{G}_{\text{RWA}}=\sum_{k}\mathcal{G}^{(k)}_{\text{RWA}}caligraphic_G start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT = ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT caligraphic_G start_POSTSUPERSCRIPT ( italic_k ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT with

𝒢RWA(k)subscriptsuperscript𝒢𝑘RWA\displaystyle\mathcal{G}^{(k)}_{\text{RWA}}caligraphic_G start_POSTSUPERSCRIPT ( italic_k ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT =i[Δ]k,k+iΩR2([Sx©]k,k+1+[Sx©]k,k1)absent𝑖subscriptdelimited-[]Δ𝑘𝑘𝑖subscriptΩ𝑅2subscriptdelimited-[]superscriptsubscript𝑆𝑥©𝑘𝑘1subscriptdelimited-[]superscriptsubscript𝑆𝑥©𝑘𝑘1\displaystyle=i\left[\Delta\right]_{k,k}+i\frac{\Omega_{R}}{2}\left([S_{x}^{% \copyright}]_{k,k+1}+[S_{x}^{\copyright}]_{k,k-1}\right)= italic_i [ roman_Δ ] start_POSTSUBSCRIPT italic_k , italic_k end_POSTSUBSCRIPT + italic_i divide start_ARG roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ( [ italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT italic_k , italic_k + 1 end_POSTSUBSCRIPT + [ italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT italic_k , italic_k - 1 end_POSTSUBSCRIPT )
+(ΓSE+ΓSD+ROP)[𝒜SD]k,kROP2[𝒜SE(z)]k,ksubscriptΓSEsubscriptΓSDsubscript𝑅OPsubscriptdelimited-[]subscript𝒜SD𝑘𝑘subscript𝑅OP2subscriptdelimited-[]superscriptsubscript𝒜SE𝑧𝑘𝑘\displaystyle+(\Gamma_{\text{SE}}+\Gamma_{\text{SD}}+R_{\text{OP}})\left[% \mathcal{A}_{\text{SD}}\right]_{k,k}-\frac{R_{\text{OP}}}{2}\left[\mathcal{A}_% {\text{SE}}^{(z)}\right]_{k,k}+ ( roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT + roman_Γ start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT + italic_R start_POSTSUBSCRIPT OP end_POSTSUBSCRIPT ) [ caligraphic_A start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT italic_k , italic_k end_POSTSUBSCRIPT - divide start_ARG italic_R start_POSTSUBSCRIPT OP end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG [ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_z ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT italic_k , italic_k end_POSTSUBSCRIPT
ΓSE(𝒜SE(k,z)+𝒜SE(k,+)+𝒜SE(k,)).subscriptΓSEsuperscriptsubscript𝒜SE𝑘𝑧superscriptsubscript𝒜SE𝑘superscriptsubscript𝒜SE𝑘\displaystyle-\Gamma_{\text{SE}}\left(\mathcal{A}_{\text{SE}}^{(k,z)}+\mathcal% {A}_{\text{SE}}^{(k,+)}+\mathcal{A}_{\text{SE}}^{(k,-)}\right).- roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT ( caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_k , italic_z ) end_POSTSUPERSCRIPT + caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_k , + ) end_POSTSUPERSCRIPT + caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_k , - ) end_POSTSUPERSCRIPT ) . (29)

In Eq. (II.2), the spin exchange due to the transverse and longitudinal spin components are

𝒜SE(k,±)superscriptsubscript𝒜SE𝑘plus-or-minus\displaystyle\mathcal{A}_{\text{SE}}^{(k,\pm)}caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_k , ± ) end_POSTSUPERSCRIPT =S~x(1)[𝒜SE(x)]k,k±1+S~y(1)[𝒜SE(y)]k,k±1absentsuperscriptsubscript~𝑆𝑥minus-or-plus1subscriptdelimited-[]superscriptsubscript𝒜SE𝑥𝑘plus-or-minus𝑘1superscriptsubscript~𝑆𝑦minus-or-plus1subscriptdelimited-[]superscriptsubscript𝒜SE𝑦𝑘plus-or-minus𝑘1\displaystyle=\tilde{S}_{x}^{(\mp 1)}\left[\mathcal{A}_{\text{SE}}^{(x)}\right% ]_{k,k\pm 1}+\tilde{S}_{y}^{(\mp 1)}\left[\mathcal{A}_{\text{SE}}^{(y)}\right]% _{k,k\pm 1}= over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ∓ 1 ) end_POSTSUPERSCRIPT [ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_x ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT italic_k , italic_k ± 1 end_POSTSUBSCRIPT + over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ∓ 1 ) end_POSTSUPERSCRIPT [ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_y ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT italic_k , italic_k ± 1 end_POSTSUBSCRIPT (30)

and

𝒜SE(k,z)superscriptsubscript𝒜SE𝑘𝑧\displaystyle\mathcal{A}_{\text{SE}}^{(k,z)}caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_k , italic_z ) end_POSTSUPERSCRIPT =S~z(0)[𝒜SE(z)]k,k.absentsuperscriptsubscript~𝑆𝑧0subscriptdelimited-[]superscriptsubscript𝒜SE𝑧𝑘𝑘\displaystyle=\tilde{S}_{z}^{(0)}\left[\mathcal{A}_{\text{SE}}^{(z)}\right]_{k% ,k}.= over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT [ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_z ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT italic_k , italic_k end_POSTSUBSCRIPT . (31)
Refer to caption
Figure 2: (a) Magnetic resonance spectrum of weakly polarized (i.e., P1much-less-than𝑃1P\ll 1italic_P ≪ 1) Rb87superscriptRb87{}^{87}{\rm Rb}start_FLOATSUPERSCRIPT 87 end_FLOATSUPERSCRIPT roman_Rb atoms in a weak field of 5μT5𝜇T5~{}{\rm\mu T}5 italic_μ roman_T. (b) The same as (a), but for a stronger field of 100μT100𝜇T100~{}{\rm\mu T}100 italic_μ roman_T. Two resonance peaks corresponds to the resonances of F=2𝐹2F=2italic_F = 2 and F=1𝐹1F=1italic_F = 1 levels, respectively. (c) The same as (a), but for a field of 500μT500𝜇T500~{}{\rm\mu T}500 italic_μ roman_T, which drives the atom in the RZR regime. Six resonance peaks are well resolved, and the peaks are broader than the weak field cases. The parameters used in the calculation are ΓSD=10HzsubscriptΓSD10Hz\Gamma_{\rm SD}=10~{}{\rm Hz}roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT = 10 roman_Hz and ΓSE=775HzsubscriptΓSE775Hz\Gamma_{\rm SE}=775~{}{\rm Hz}roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT = 775 roman_Hz. A very weak optical pumping rate Rop=1Hzsubscript𝑅op1HzR_{\rm op}=1~{}{\rm Hz}italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT = 1 roman_Hz is assumed to ensure the weak spin polarization.

II.3 Resolved and unresolved Zeeman resonance regimes

With the RWA, the evolution superoperator 𝒢RWAsubscript𝒢RWA\mathcal{G}_{\rm RWA}caligraphic_G start_POSTSUBSCRIPT roman_RWA end_POSTSUBSCRIPT is time-independent. The magnetic resonance spectrum is obtained from the steady-state solution of Eq. (28). Figure 2 shows the numerical results of the spectrum of weakly polarized Rb87superscriptRb87{}^{87}{\rm Rb}start_FLOATSUPERSCRIPT 87 end_FLOATSUPERSCRIPT roman_Rb atoms in different magnetic fields.

In a weak field, the frequency difference between the k𝑘kitalic_kth order Zeeman coherences is negligible, and the detuning matrix [Δ]k,ksubscriptdelimited-[]Δ𝑘𝑘\left[\Delta\right]_{k,k}[ roman_Δ ] start_POSTSUBSCRIPT italic_k , italic_k end_POSTSUBSCRIPT in Eq. (II.2) is proportional to an identity matrix, i.e., [Δ]k,k=Δkksubscriptdelimited-[]Δ𝑘𝑘subscriptΔ𝑘subscript𝑘\left[\Delta\right]_{k,k}=\Delta_{k}\mathcal{I}_{k}[ roman_Δ ] start_POSTSUBSCRIPT italic_k , italic_k end_POSTSUBSCRIPT = roman_Δ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT caligraphic_I start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT. In this case, all Zeeman resonances are degenerate, and only one resonance peak appears in the spectrum [see Fig. 2(a)].

As the magnetic field increases, the frequency degeneracy of Zeeman coherences in subspaces with F=a𝐹𝑎F=aitalic_F = italic_a and F=b𝐹𝑏F=bitalic_F = italic_b is lifted because of the slight difference of their gyromagnetic ratios. Taking Rb87superscriptRb87{}^{87}{\rm Rb}start_FLOATSUPERSCRIPT 87 end_FLOATSUPERSCRIPT roman_Rb atoms as an example, the difference between the gyromagnetic ratios of F=2𝐹2F=2italic_F = 2 and F=1𝐹1F=1italic_F = 1 levels is Δγ=|γb||γa|=2π×28Hz/μTΔ𝛾subscript𝛾𝑏subscript𝛾𝑎2𝜋28Hz𝜇T\Delta\gamma=|\gamma_{b}|-|\gamma_{a}|=2\pi\times 28~{}{\rm Hz/\mu T}roman_Δ italic_γ = | italic_γ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT | - | italic_γ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT | = 2 italic_π × 28 roman_Hz / italic_μ roman_T, which causes a splitting of 2π×2.8kHz2𝜋2.8kHz2\pi\times 2.8~{}{\rm kHz}2 italic_π × 2.8 roman_kHz in a field of 100μT100𝜇T100~{}{\rm\mu T}100 italic_μ roman_T [see Fig. 2(b)]. In this field, the Zeeman coherences within the F=1𝐹1F=1italic_F = 1 or F=2𝐹2F=2italic_F = 2 subspaces can still be regarded as degenerate, and the spectrum consists of two peaks.

When the magnetic field is so strong that the quadratic splitting ωQ=2(γeB)2/Ωhfsubscript𝜔𝑄2superscriptsubscript𝛾𝑒𝐵2subscriptΩhf\omega_{Q}=2(\gamma_{e}B)^{2}/\Omega_{\rm hf}italic_ω start_POSTSUBSCRIPT italic_Q end_POSTSUBSCRIPT = 2 ( italic_γ start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_B ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / roman_Ω start_POSTSUBSCRIPT roman_hf end_POSTSUBSCRIPT well exceeds the typical resonance linewidth [9, 10], the system enters the RZR regime. Figure 2(c) shows the example of Rb87superscriptRb87{}^{87}{\rm Rb}start_FLOATSUPERSCRIPT 87 end_FLOATSUPERSCRIPT roman_Rb atoms in a magnetic field B0=500μTsubscript𝐵0500𝜇TB_{0}=500~{}{\rm\mu T}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 500 italic_μ roman_T, where the spectrum consists of six peaks with quadratic splitting ωQ=2π×3.6kHzsubscript𝜔𝑄2𝜋3.6kHz\omega_{Q}=2\pi\times 3.6~{}{\rm kHz}italic_ω start_POSTSUBSCRIPT italic_Q end_POSTSUBSCRIPT = 2 italic_π × 3.6 roman_kHz.

As shown in Fig. 2, the resonance linewidth in the RZR regime is much broader than that in the UZR regime. The linewidth in the RZR regime [Fig. 2(c)] has been well studied theoretically [8] and experimentally [7]. In the following, we will focus on the linewidth in the UZR regime.

II.4 Weak driving approximation

II.4.1 Weak driving approximation in UZR regime

With the operator 𝒫1subscript𝒫1\mathcal{P}_{1}caligraphic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT applied on both sides of Eqs. (28), the master equation satisfied by the first order Zeeman coherence |ρ~1)|\tilde{\rho}_{1})| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) is

ddt|ρ~1)=\displaystyle\frac{d}{dt}|\tilde{\rho}_{1})=divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = [𝒢RWA]1,1|ρ~1)[𝒢RWA]1,0|ρ~0)[𝒢RWA]1,2|ρ~2),\displaystyle-\left[\mathcal{G}_{\text{RWA}}\right]_{1,1}|\tilde{\rho}_{1})-% \left[\mathcal{G}_{\text{RWA}}\right]_{1,0}|\tilde{\rho}_{0})-\left[\mathcal{G% }_{\text{RWA}}\right]_{1,2}|\tilde{\rho}_{2}),- [ caligraphic_G start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) - [ caligraphic_G start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 0 end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) - [ caligraphic_G start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) , (32)

where the diagonal block is

[𝒢RWA]1,1=subscriptdelimited-[]subscript𝒢RWA11absent\displaystyle\left[\mathcal{G}_{\text{RWA}}\right]_{1,1}=[ caligraphic_G start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT = iΔ1+Γtot[𝒜SD]1,1(ROP2+ΓSES~z(0))[𝒜SE(z)]1,1,𝑖subscriptΔ1subscriptΓtotsubscriptdelimited-[]subscript𝒜SD11subscript𝑅OP2subscriptΓSEsuperscriptsubscript~𝑆𝑧0subscriptdelimited-[]subscriptsuperscript𝒜𝑧SE11\displaystyle i\Delta_{1}+\Gamma_{\text{tot}}\left[\mathcal{A}_{\text{SD}}% \right]_{1,1}-\left(\frac{R_{\text{OP}}}{2}+\Gamma_{\text{SE}}\tilde{S}_{z}^{(% 0)}\right)\left[\mathcal{A}^{(z)}_{\text{SE}}\right]_{1,1},italic_i roman_Δ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + roman_Γ start_POSTSUBSCRIPT tot end_POSTSUBSCRIPT [ caligraphic_A start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - ( divide start_ARG italic_R start_POSTSUBSCRIPT OP end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG + roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ) [ caligraphic_A start_POSTSUPERSCRIPT ( italic_z ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT , (33)

and the off-diagonal blocks

[𝒢RWA]1,0=subscriptdelimited-[]subscript𝒢RWA10absent\displaystyle\left[\mathcal{G}_{\text{RWA}}\right]_{1,0}=[ caligraphic_G start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 0 end_POSTSUBSCRIPT = iΩR2[Sx©]1,0ΓSE𝒜SE(k,),𝑖subscriptΩ𝑅2subscriptdelimited-[]superscriptsubscript𝑆𝑥©10subscriptΓSEsuperscriptsubscript𝒜SE𝑘\displaystyle\frac{i\Omega_{R}}{2}\left[S_{x}^{\copyright}\right]_{1,0}-\Gamma% _{\text{SE}}\mathcal{A}_{\text{SE}}^{(k,-)},divide start_ARG italic_i roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG [ italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 0 end_POSTSUBSCRIPT - roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_k , - ) end_POSTSUPERSCRIPT , (34)
[𝒢RWA]1,2=subscriptdelimited-[]subscript𝒢RWA12absent\displaystyle\left[\mathcal{G}_{\text{RWA}}\right]_{1,2}=[ caligraphic_G start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT = iΩR2[Sx©]1,2ΓSE𝒜SE(k,+),𝑖subscriptΩ𝑅2subscriptdelimited-[]superscriptsubscript𝑆𝑥©12subscriptΓSEsuperscriptsubscript𝒜SE𝑘\displaystyle\frac{i\Omega_{R}}{2}\left[S_{x}^{\copyright}\right]_{1,2}-\Gamma% _{\text{SE}}\mathcal{A}_{\text{SE}}^{(k,+)},divide start_ARG italic_i roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG [ italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT - roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_k , + ) end_POSTSUPERSCRIPT , (35)

couple the first order Zeeman coherence |ρ~1)|\tilde{\rho}_{1})| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) to the population |ρ~0)|\tilde{\rho}_{0})| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) and the second order Zeeman coherence |ρ~2)|\tilde{\rho}_{2})| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ), due to the RF driving field and the transverse components of the spin exchange. In Eq. (33), we have replaced the detuning matrix [Δ]1,1subscriptdelimited-[]Δ11\left[\Delta\right]_{1,1}[ roman_Δ ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT by the detuning frequency Δ1subscriptΔ1\Delta_{1}roman_Δ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (the identity matrix is omitted for simplicity) in the UZR regime.

In a weak RF driving field, the k𝑘kitalic_kth order Zeeman coherence |ρ~k)|\tilde{\rho}_{k})| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is on the order of (ΩR/Γ2)ksuperscriptsubscriptΩ𝑅subscriptΓ2𝑘(\Omega_{R}/\Gamma_{2})^{k}( roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT / roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT [8], where Γ2subscriptΓ2\Gamma_{2}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is the resonance linewidth. Although the exact value of Γ2subscriptΓ2\Gamma_{2}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is unknown before solving the master equation, it is reasonable that Γ2subscriptΓ2\Gamma_{2}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is in the same order as the relaxation rates Γ2ΓSEsimilar-tosubscriptΓ2subscriptΓSE\Gamma_{2}\sim\Gamma_{\text{SE}}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ∼ roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT, ΓSDsubscriptΓSD\Gamma_{\text{SD}}roman_Γ start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT or ROPsubscript𝑅OPR_{\text{OP}}italic_R start_POSTSUBSCRIPT OP end_POSTSUBSCRIPT. We can always control the strength of the RF driving field, so that ΩRΓ2much-less-thansubscriptΩ𝑅subscriptΓ2\Omega_{R}\ll\Gamma_{2}roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ≪ roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is well satisfied. In this case, the population |ρ~0)|\tilde{\rho}_{0})| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) is approximately constant, and we replace |ρ~0)|\tilde{\rho}_{0})| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) in Eq. (32) by the equilibrium state |ρeq)|\rho_{\text{eq}})| italic_ρ start_POSTSUBSCRIPT eq end_POSTSUBSCRIPT ) in the absence of RF driving (with optical pumping and spin relaxation only), i.e.,

|ρ~0)=𝒫0|ρeq)|ρ¯0).|\tilde{\rho}_{0})=\mathcal{P}_{0}|\rho_{\text{eq}})\equiv|\bar{\rho}_{0}).| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = caligraphic_P start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | italic_ρ start_POSTSUBSCRIPT eq end_POSTSUBSCRIPT ) ≡ | over¯ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) . (36)

Furthermore, the second order Zeeman coherence |ρ~2)|\tilde{\rho}_{2})| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) in Eq. (32) is neglected. With this weak driving approximation (WDA) Eq. (32) becomes

ddt|ρ~1)\displaystyle\frac{d}{dt}|\tilde{\rho}_{1})divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) =([𝒢RWA]1,1ΓSE[𝒜SE()]1,1)|ρ~1)ΩR|ν¯1)\displaystyle=-\left(\left[\mathcal{G}_{\text{RWA}}\right]_{1,1}-\Gamma_{\rm SE% }\left[\mathcal{A}_{\text{SE}}^{(\perp)}\right]_{1,1}\right)|\tilde{\rho}_{1})% -\Omega_{R}|\bar{\nu}_{1})= - ( [ caligraphic_G start_POSTSUBSCRIPT RWA end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT [ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT ) | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) - roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT | over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT )
𝒢WDA|ρ~1)ΩR|ν¯1)\displaystyle\equiv-\mathcal{G}_{\text{WDA}}|\tilde{\rho}_{1})-\Omega_{R}|\bar% {\nu}_{1})≡ - caligraphic_G start_POSTSUBSCRIPT WDA end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) - roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT | over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) (37)

In Eq. (II.4.1), the superoperator due to the spin exchange of transverse components is

[𝒜SE()]1,1subscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to11\displaystyle\left[\mathcal{A}_{\text{SE}}^{(\perp)}\right]_{1,1}[ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT =[𝒜SE(x)]1,0[𝒬x]0,1+[𝒜SE(y)]1,0[𝒬y]0,1,absentsubscriptdelimited-[]superscriptsubscript𝒜SE𝑥10subscriptdelimited-[]subscript𝒬𝑥01subscriptdelimited-[]superscriptsubscript𝒜SE𝑦10subscriptdelimited-[]subscript𝒬𝑦01\displaystyle=\left[\mathcal{A}_{\text{SE}}^{(x)}\right]_{1,0}\left[\mathcal{Q% }_{x}\right]_{0,1}+\left[\mathcal{A}_{\text{SE}}^{(y)}\right]_{1,0}\left[% \mathcal{Q}_{y}\right]_{0,1},= [ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_x ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 0 end_POSTSUBSCRIPT [ caligraphic_Q start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 0 , 1 end_POSTSUBSCRIPT + [ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_y ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 0 end_POSTSUBSCRIPT [ caligraphic_Q start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 0 , 1 end_POSTSUBSCRIPT , (38)

where [𝒬x]0,1subscriptdelimited-[]subscript𝒬𝑥01\left[\mathcal{Q}_{x}\right]_{0,1}[ caligraphic_Q start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 0 , 1 end_POSTSUBSCRIPT and [𝒬y]0,1subscriptdelimited-[]subscript𝒬𝑦01\left[\mathcal{Q}_{y}\right]_{0,1}[ caligraphic_Q start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 0 , 1 end_POSTSUBSCRIPT are matrices associated with the equilibrium state |ρeq)|\rho_{\text{eq}})| italic_ρ start_POSTSUBSCRIPT eq end_POSTSUBSCRIPT ) as

[𝒬x,y]0,1=𝒫0|ρeq)(Sx,y|𝒫1.\left[\mathcal{Q}_{x,y}\right]_{0,1}=\mathcal{P}_{0}|\rho_{\text{eq}})(S_{x,y}% |\mathcal{P}_{1}.[ caligraphic_Q start_POSTSUBSCRIPT italic_x , italic_y end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 0 , 1 end_POSTSUBSCRIPT = caligraphic_P start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | italic_ρ start_POSTSUBSCRIPT eq end_POSTSUBSCRIPT ) ( italic_S start_POSTSUBSCRIPT italic_x , italic_y end_POSTSUBSCRIPT | caligraphic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT . (39)

The inhomogeneous term of Eq. (II.4.1) arises from the RF driving field, with the state vector |ν¯1)|\bar{\nu}_{1})| over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) given by

|ν¯1)=i2[Sx©]1,0|ρ¯0).|\bar{\nu}_{1})=\frac{i}{2}\left[S_{x}^{\copyright}\right]_{1,0}|\bar{\rho}_{0% }).| over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = divide start_ARG italic_i end_ARG start_ARG 2 end_ARG [ italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 0 end_POSTSUBSCRIPT | over¯ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) . (40)

Notice that, with the WDA, Eq. (II.4.1) is a closed linear equation for |ρ~1)|\tilde{\rho}_{1})| over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) and is used to analyze the resonance linewidth.

II.4.2 Spin temperature distribution

In general, the WDA does not have requirement of the specific form of the population |ρ¯0)|\bar{\rho}_{0})| over¯ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ). Here we consider the spin temperature distribution, which is widely used in the study of alkali-metal spin dynamics [8]. The spin temperature distribution density matrix is

ρeq=eβFzZ=eβIzeβSzZIZS,subscript𝜌eqsuperscript𝑒𝛽subscript𝐹𝑧𝑍superscript𝑒𝛽subscript𝐼𝑧superscript𝑒𝛽subscript𝑆𝑧subscript𝑍𝐼subscript𝑍𝑆\displaystyle\rho_{\rm eq}=\frac{e^{\beta F_{z}}}{Z}=\frac{e^{\beta I_{z}}e^{% \beta S_{z}}}{Z_{I}Z_{S}},italic_ρ start_POSTSUBSCRIPT roman_eq end_POSTSUBSCRIPT = divide start_ARG italic_e start_POSTSUPERSCRIPT italic_β italic_F start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_POSTSUPERSCRIPT end_ARG start_ARG italic_Z end_ARG = divide start_ARG italic_e start_POSTSUPERSCRIPT italic_β italic_I start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_β italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_POSTSUPERSCRIPT end_ARG start_ARG italic_Z start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT italic_Z start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT end_ARG , (41)

where Fz=Sz+Izsubscript𝐹𝑧subscript𝑆𝑧subscript𝐼𝑧F_{z}=S_{z}+I_{z}italic_F start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT + italic_I start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT, ZJ=m=JJexp(βm)subscript𝑍𝐽superscriptsubscript𝑚𝐽𝐽𝛽𝑚Z_{J}=\sum_{m=-J}^{J}\exp(\beta m)italic_Z start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = ∑ start_POSTSUBSCRIPT italic_m = - italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_J end_POSTSUPERSCRIPT roman_exp ( italic_β italic_m ) and Z=ZIZS=Tr[exp(βFz)]𝑍subscript𝑍𝐼subscript𝑍𝑆Trdelimited-[]𝛽subscript𝐹𝑧Z=Z_{I}Z_{S}={\rm Tr}[\exp(\beta F_{z})]italic_Z = italic_Z start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT italic_Z start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = roman_Tr [ roman_exp ( italic_β italic_F start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ) ] is the normalization factor. The parameter β𝛽\betaitalic_β is related to the electron spin polarization P𝑃Pitalic_P by

P=2Sz=2S~z(0)=tanh(β2).𝑃2delimited-⟨⟩subscript𝑆𝑧2superscriptsubscript~𝑆𝑧0𝛽2\displaystyle P=2\langle S_{z}\rangle=2\tilde{S}_{z}^{(0)}=\tanh\left(\frac{% \beta}{2}\right).italic_P = 2 ⟨ italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ⟩ = 2 over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT = roman_tanh ( divide start_ARG italic_β end_ARG start_ARG 2 end_ARG ) . (42)

The probability of occupation pmsubscript𝑝𝑚p_{m}italic_p start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT of the spin state |F,mket𝐹𝑚|F,m\rangle| italic_F , italic_m ⟩ is a function of spin polarization P𝑃Pitalic_P

pm(P)=P(1+P)[I]/2+m(1P)[I]/2m(1+P)[I](1P)[I],subscript𝑝𝑚𝑃𝑃superscript1𝑃delimited-[]𝐼2𝑚superscript1𝑃delimited-[]𝐼2𝑚superscript1𝑃delimited-[]𝐼superscript1𝑃delimited-[]𝐼\displaystyle p_{m}(P)=\frac{P(1+P)^{[I]/2+m}(1-P)^{[I]/2-m}}{(1+P)^{[I]}-(1-P% )^{[I]}},italic_p start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ( italic_P ) = divide start_ARG italic_P ( 1 + italic_P ) start_POSTSUPERSCRIPT [ italic_I ] / 2 + italic_m end_POSTSUPERSCRIPT ( 1 - italic_P ) start_POSTSUPERSCRIPT [ italic_I ] / 2 - italic_m end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 + italic_P ) start_POSTSUPERSCRIPT [ italic_I ] end_POSTSUPERSCRIPT - ( 1 - italic_P ) start_POSTSUPERSCRIPT [ italic_I ] end_POSTSUPERSCRIPT end_ARG , (43)

where [I]=2I+1delimited-[]𝐼2𝐼1[I]=2I+1[ italic_I ] = 2 italic_I + 1. As the state ρeq=ρeq(P)subscript𝜌eqsubscript𝜌eq𝑃\rho_{\rm eq}=\rho_{\rm eq}(P)italic_ρ start_POSTSUBSCRIPT roman_eq end_POSTSUBSCRIPT = italic_ρ start_POSTSUBSCRIPT roman_eq end_POSTSUBSCRIPT ( italic_P ) is determined by P𝑃Pitalic_P, the matrices 𝒬x,y=𝒬x,y(P)subscript𝒬𝑥𝑦subscript𝒬𝑥𝑦𝑃\mathcal{Q}_{x,y}=\mathcal{Q}_{x,y}(P)caligraphic_Q start_POSTSUBSCRIPT italic_x , italic_y end_POSTSUBSCRIPT = caligraphic_Q start_POSTSUBSCRIPT italic_x , italic_y end_POSTSUBSCRIPT ( italic_P ) defined in Eq. (39), together with the superoperator [𝒜SE()]1,1=[𝒜SE()(P)]1,1subscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to11subscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to𝑃11\left[\mathcal{A}_{\text{SE}}^{(\perp)}\right]_{1,1}=\left[\mathcal{A}_{\text{% SE}}^{(\perp)}(P)\right]_{1,1}[ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT = [ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ( italic_P ) ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT defined in Eq. (38), are also functions of P𝑃Pitalic_P.

With the spin temperature distribution, the dimensionless WDA superoperator 𝔾WDAsubscript𝔾WDA\mathbb{G}_{\text{WDA}}blackboard_G start_POSTSUBSCRIPT WDA end_POSTSUBSCRIPT (in the unit of ΓSDsubscriptΓSD\Gamma_{\rm SD}roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT) becomes

𝔾WDA(η,P)𝒢WDAΓSDsubscript𝔾WDA𝜂𝑃subscript𝒢WDAsubscriptΓSD\displaystyle\mathbb{G}_{\rm WDA}(\eta,P)\equiv\frac{\mathcal{G}_{\text{WDA}}}% {\Gamma_{\rm SD}}blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT ( italic_η , italic_P ) ≡ divide start_ARG caligraphic_G start_POSTSUBSCRIPT WDA end_POSTSUBSCRIPT end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG (44)
=\displaystyle== (η+11P)([𝒜SD]1,1P2[𝒜SE(z)]1,1)η[𝒜SE()(P)]1,1,𝜂11𝑃subscriptdelimited-[]subscript𝒜SD11𝑃2subscriptdelimited-[]superscriptsubscript𝒜SE𝑧11𝜂subscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to𝑃11\displaystyle\left(\eta+\frac{1}{1-P}\right)\left(\left[\mathcal{A}_{\text{SD}% }\right]_{1,1}-\frac{P}{2}\left[\mathcal{A}_{\text{SE}}^{(z)}\right]_{1,1}% \right)-\eta\left[\mathcal{A}_{\text{SE}}^{(\perp)}(P)\right]_{1,1},( italic_η + divide start_ARG 1 end_ARG start_ARG 1 - italic_P end_ARG ) ( [ caligraphic_A start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - divide start_ARG italic_P end_ARG start_ARG 2 end_ARG [ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_z ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT ) - italic_η [ caligraphic_A start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ( italic_P ) ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT ,

where η=ΓSE/ΓSD𝜂subscriptΓSEsubscriptΓSD\eta=\Gamma_{\rm SE}/\Gamma_{\rm SD}italic_η = roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT / roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT. In Eq. (44), we have used the fact that the spin polarization P𝑃Pitalic_P of the spin temperature state is

P=RopRop+ΓSD,𝑃subscript𝑅opsubscript𝑅opsubscriptΓSD\displaystyle P=\frac{R_{\rm op}}{R_{\rm op}+\Gamma_{\rm SD}},italic_P = divide start_ARG italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT end_ARG start_ARG italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT + roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG , (45)

or, equivalently, the dimensionless pumping rate is expressed in terms of P𝑃Pitalic_P as

RopΓSD=P1P.subscript𝑅opsubscriptΓSD𝑃1𝑃\displaystyle\frac{R_{\rm op}}{\Gamma_{\rm SD}}=\frac{P}{1-P}.divide start_ARG italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG = divide start_ARG italic_P end_ARG start_ARG 1 - italic_P end_ARG . (46)

Notice that the dimensionless WDA superoperator 𝔾WDA(η,P)subscript𝔾WDA𝜂𝑃\mathbb{G}_{\text{WDA}}(\eta,P)blackboard_G start_POSTSUBSCRIPT WDA end_POSTSUBSCRIPT ( italic_η , italic_P ) is fully characterized by two dimensionless parameters, namely, the relative strength η𝜂\etaitalic_η of spin exchange and the spin polarization P𝑃Pitalic_P.

II.4.3 Observable and lineshape

In typical magnetometer experiments, a linearly polarized laser beam is used to measure the transverse spin component, e.g. Sxsubscript𝑆𝑥S_{x}italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT, via the Faraday rotation effect. According to the discussion above, the time-dependent observable Sx(t)delimited-⟨⟩subscript𝑆𝑥𝑡\langle S_{x}(t)\rangle⟨ italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ( italic_t ) ⟩ is

Sx(t)=(Sx|ρZ(t))=S~x(1)eiωt+c.c,delimited-⟨⟩subscript𝑆𝑥𝑡conditionalsubscript𝑆𝑥subscript𝜌Z𝑡superscriptsubscript~𝑆𝑥1superscript𝑒𝑖𝜔𝑡c.c\displaystyle\langle S_{x}(t)\rangle=(S_{x}|\rho_{\rm Z}(t))=\tilde{S}_{x}^{(1% )}e^{-i\omega t}+\text{c.c},⟨ italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ( italic_t ) ⟩ = ( italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT | italic_ρ start_POSTSUBSCRIPT roman_Z end_POSTSUBSCRIPT ( italic_t ) ) = over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ω italic_t end_POSTSUPERSCRIPT + c.c , (47)

and the in-phase and out-of-phase quadratures are

X𝑋\displaystyle Xitalic_X =Re[S~x(1)]=Re[(Sx|ρ~1)],absentRedelimited-[]superscriptsubscript~𝑆𝑥1Redelimited-[]conditionalsubscript𝑆𝑥subscript~𝜌1\displaystyle=\text{Re}\left[\tilde{S}_{x}^{(1)}\right]=\text{Re}\left[(S_{x}|% \tilde{\rho}_{1})\right],= Re [ over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ] = Re [ ( italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ] , (48)
Y𝑌\displaystyle Yitalic_Y =Im[S~x(1)]=Im[(Sx|ρ~1)].absentImdelimited-[]superscriptsubscript~𝑆𝑥1Imdelimited-[]conditionalsubscript𝑆𝑥subscript~𝜌1\displaystyle=-\text{Im}\left[\tilde{S}_{x}^{(1)}\right]=-\text{Im}\left[(S_{x% }|\tilde{\rho}_{1})\right].= - Im [ over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ] = - Im [ ( italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ] . (49)
Refer to caption
Figure 3: (a) The eigenvalues of 𝔾WDAsubscript𝔾WDA\mathbb{G}_{\rm{WDA}}blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT as functions of P𝑃Pitalic_P with η=100𝜂100\eta=100italic_η = 100. The lowest curve represents the smallest eigenvalue λ1subscript𝜆1\lambda_{1}italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, while the dashed curves are λ2,λ3,,λ6subscript𝜆2subscript𝜆3subscript𝜆6\lambda_{2},\lambda_{3},\dots,\lambda_{6}italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_λ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT , … , italic_λ start_POSTSUBSCRIPT 6 end_POSTSUBSCRIPT. (b) The normalized weight w~k=|wk|/i|wi|subscript~𝑤𝑘subscript𝑤𝑘subscript𝑖subscript𝑤𝑖\tilde{w}_{k}=|w_{k}|/\sum_{i}|w_{i}|over~ start_ARG italic_w end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = | italic_w start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT | / ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | italic_w start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT |. The red solid curve is the weight corresponding to λ1subscript𝜆1\lambda_{1}italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT. Dashed curves are much smaller weights corresponding to other eigenvalues. (c) and (d) The real and imaginary part of the observable S𝑆Sitalic_S, normalized to the maximum value. The solid curve is the sum of contributions of all eigenvalues [see Eq. (52)], while the dashed curve is the only contribution corresponding to the smallest eigenvalue λ1subscript𝜆1\lambda_{1}italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT.
Refer to caption
Figure 4: (a) The linewidth of Rb87superscriptRb87{}^{87}{\rm Rb}start_FLOATSUPERSCRIPT 87 end_FLOATSUPERSCRIPT roman_Rb atom as a function of spin polarization for different values of η𝜂\etaitalic_η. (b) Zoom-in of the low-polarization part of the curve corresponding to η=50𝜂50\eta=50italic_η = 50. The solid curve is the exact numerical result, and the dashed curve is the analytic result up to the second order correction of P𝑃Pitalic_P [see Eqs. (64)-(68)]. (c) Zoom-in of the high-polarization part. The solid curve is the exact numerical result, and the dashed curve is the analytic result of Eq. (80).

Assume the right- and left-eigenstates of the superoperator 𝒢WDAsubscript𝒢WDA\mathcal{G}_{\text{WDA}}caligraphic_G start_POSTSUBSCRIPT WDA end_POSTSUBSCRIPT, sharing the same eigenvalue iΔ1+λj𝑖subscriptΔ1subscript𝜆𝑗i\Delta_{1}+\lambda_{j}italic_i roman_Δ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT, are denoted as |λj)\left|\lambda_{j}\right)| italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) and {λj|\left\{\lambda_{j}\right|{ italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT |, respectively,

𝒢WDA|λj)=(iΔ1+λj)|λj),\displaystyle\mathcal{G}_{\text{WDA}}\left|\lambda_{j}\right)=\left(i\Delta_{1% }+\lambda_{j}\right)\left|\lambda_{j}\right),caligraphic_G start_POSTSUBSCRIPT WDA end_POSTSUBSCRIPT | italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) = ( italic_i roman_Δ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) | italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) , (50)
{λj|𝒢WDA=(iΔ1+λj){λj|.\displaystyle\left\{\lambda_{j}\right|\mathcal{G}_{\text{WDA}}=\left(i\Delta_{% 1}+\lambda_{j}\right)\left\{\lambda_{j}\right|.{ italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | caligraphic_G start_POSTSUBSCRIPT WDA end_POSTSUBSCRIPT = ( italic_i roman_Δ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) { italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | . (51)

The observable is decomposed as

S=(Sx|ρ~1)=j(Sx|λj){λj|ρ~1)=iΩR2jwjiΔj+λj.\displaystyle S=\left(S_{x}|\tilde{\rho}_{1}\right)=\sum_{j}\left(S_{x}\big{|}% \lambda_{j}\right)\left\{\lambda_{j}\big{|}\tilde{\rho}_{1}\right)=-\frac{i% \Omega_{R}}{2}\sum_{j}\frac{w_{j}}{i\Delta_{j}+\lambda_{j}}.italic_S = ( italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT | italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) { italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | over~ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = - divide start_ARG italic_i roman_Ω start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ∑ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT divide start_ARG italic_w start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_ARG start_ARG italic_i roman_Δ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_ARG . (52)

where the weight factor wjsubscript𝑤𝑗w_{j}italic_w start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT is

wj={λj|[Sx©]1,0[𝒬x]0,1|λj).w_{j}=\left\{\lambda_{j}\big{|}[S_{x}^{\copyright}]_{1,0}\left[\mathcal{Q}_{x}% \right]_{0,1}\big{|}\lambda_{j}\right).italic_w start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = { italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | [ italic_S start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT © end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 0 end_POSTSUBSCRIPT [ caligraphic_Q start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 0 , 1 end_POSTSUBSCRIPT | italic_λ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) . (53)

As shown in Fig. 3, the normalized weight factors are actually dominated by the one corresponding to the smallest eigenvalue λminsubscript𝜆min\lambda_{\rm min}italic_λ start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT of 𝒢WDAsubscript𝒢WDA\mathcal{G}_{\text{WDA}}caligraphic_G start_POSTSUBSCRIPT WDA end_POSTSUBSCRIPT. The resonant line shape is very well approximated by a single Lorentizan function, whose linewidth is the real part of the smallest eigenvalue of 𝒢WDAsubscript𝒢WDA\mathcal{G}_{\text{WDA}}caligraphic_G start_POSTSUBSCRIPT WDA end_POSTSUBSCRIPT. In the following, we will focus on the smallest eigenvalue of 𝒢WDAsubscript𝒢WDA\mathcal{G}_{\text{WDA}}caligraphic_G start_POSTSUBSCRIPT WDA end_POSTSUBSCRIPT, and study its behavior under different spin polarization and spin exchange rate.

III Results and Discussion

Figure 4(a) illustrates the overall behavior of the Rb87superscriptRb87{}^{87}{\rm Rb}start_FLOATSUPERSCRIPT 87 end_FLOATSUPERSCRIPT roman_Rb linewidth as a function of pumping rate and spin exchange rate. At low spin exchange rates (e.g., η=1𝜂1\eta=1italic_η = 1), the linewidth monotonically increases with pumping rate. In contrast, at high spin exchange rates (η>5𝜂5\eta>5italic_η > 5), the linewidth initially decreases with increasing pumping rate before reaching a minimum value and then increases linearly with an increase in pumping rate. This nonmonotonic behavior is the consequence of the competition between the narrowing effect due to the spin exchange at high spin polarization and the broadening induced by strong optical pumping. Indeed, the light narrowing effect at high spin exchange rate was observed in the strong magnetic field cases, where the Zeeman resonances are well resolved [6, 7]. To reveal the difference between the Zeeman resonances resolved and unresolved cases, we analyze the linewidth of alkali atoms in the low and high polarization limits, and compare with the result of previous studies.

III.1 Low polarization limit (P1much-less-than𝑃1P\ll 1italic_P ≪ 1)

To analyze the linewidth in the low polarization limit, the dimensionless matrix Eq. (44) is expanded in the power series of P𝑃Pitalic_P up to the second order as

𝔾WDA(η,P)=𝔾WDA(0)+P𝔾WDA(1)+P2𝔾WDA(2),subscript𝔾WDA𝜂𝑃superscriptsubscript𝔾WDA0𝑃superscriptsubscript𝔾WDA1superscript𝑃2superscriptsubscript𝔾WDA2\mathbb{G}_{\rm WDA}(\eta,P)=\mathbb{G}_{\rm WDA}^{(0)}+P\mathbb{G}_{\rm WDA}^% {(1)}+P^{2}\mathbb{G}_{\rm WDA}^{(2)},blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT ( italic_η , italic_P ) = blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + italic_P blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT + italic_P start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT , (54)

where

𝔾WDA(0)=superscriptsubscript𝔾WDA0absent\displaystyle\mathbb{G}_{\mathrm{WDA}}^{(0)}=blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT = (η+1)[𝒜SD]1,1η[𝒜SE()]1,1(0),𝜂1subscriptdelimited-[]subscript𝒜SD11𝜂superscriptsubscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to110\displaystyle\left(\eta+1\right)\left[\mathcal{A}_{\mathrm{SD}}\right]_{1,1}-% \eta\left[\mathcal{A}_{\mathrm{SE}}^{(\perp)}\right]_{1,1}^{(0)},( italic_η + 1 ) [ caligraphic_A start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - italic_η [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT , (55)
𝔾WDA(1)=superscriptsubscript𝔾WDA1absent\displaystyle\mathbb{G}_{\mathrm{WDA}}^{(1)}=blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = [𝒜SD]1,1η+12[𝒜SE(z)]1,1η[𝒜SE()]1,1(1),subscriptdelimited-[]subscript𝒜SD11𝜂12subscriptdelimited-[]superscriptsubscript𝒜SE𝑧11𝜂superscriptsubscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to111\displaystyle\left[\mathcal{A}_{\mathrm{SD}}\right]_{1,1}-\frac{\eta+1}{2}% \left[\mathcal{A}_{\mathrm{SE}}^{(z)}\right]_{1,1}-\eta\left[\mathcal{A}_{% \mathrm{SE}}^{(\perp)}\right]_{1,1}^{(1)},[ caligraphic_A start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - divide start_ARG italic_η + 1 end_ARG start_ARG 2 end_ARG [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_z ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - italic_η [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT , (56)
𝔾WDA(2)=superscriptsubscript𝔾WDA2absent\displaystyle\mathbb{G}_{\mathrm{WDA}}^{(2)}=blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT = [𝒜SD]1,112[𝒜SE(z)]1,1η[𝒜SE()]1,1(2).subscriptdelimited-[]subscript𝒜SD1112subscriptdelimited-[]superscriptsubscript𝒜SE𝑧11𝜂superscriptsubscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to112\displaystyle\left[\mathcal{A}_{\mathrm{SD}}\right]_{1,1}-\frac{1}{2}\left[% \mathcal{A}_{\mathrm{SE}}^{(z)}\right]_{1,1}-\eta\left[\mathcal{A}_{\mathrm{SE% }}^{(\perp)}\right]_{1,1}^{(2)}.[ caligraphic_A start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_z ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - italic_η [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT . (57)

In Eqs. (55)-(57), we have expanded the matrix [𝒜SE()(P)]1,1subscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to𝑃11\left[\mathcal{A}_{\rm SE}^{(\perp)}(P)\right]_{1,1}[ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ( italic_P ) ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT to the second order of P𝑃Pitalic_P, i.e.,

[𝒜SE()(P)]1,1=[𝒜SE()]1,1(0)+P[𝒜SE()]1,1(1)+P2[𝒜SE()]1,1(2).subscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to𝑃11superscriptsubscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to110𝑃superscriptsubscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to111superscript𝑃2superscriptsubscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to112\left[\mathcal{A}_{\mathrm{SE}}^{(\perp)}(P)\right]_{1,1}=\left[\mathcal{A}_{% \mathrm{SE}}^{(\perp)}\right]_{1,1}^{(0)}+P\cdot\left[\mathcal{A}_{\mathrm{SE}% }^{(\perp)}\right]_{1,1}^{(1)}+P^{2}\cdot\left[\mathcal{A}_{\mathrm{SE}}^{(% \perp)}\right]_{1,1}^{(2)}.[ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ( italic_P ) ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT = [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + italic_P ⋅ [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT + italic_P start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ⋅ [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT . (58)

The leading order matrix 𝔾WDA(0)superscriptsubscript𝔾WDA0\mathbb{G}_{\rm WDA}^{(0)}blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT is Hermitian, whose eigenvalue λk(0)superscriptsubscript𝜆𝑘0\lambda_{k}^{(0)}italic_λ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT and corresponding eigenvector |vk(0))\left|v_{k}^{(0)}\right)| italic_v start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ) are obtained by solving the eigenvalue problem

𝔾WDA(0)|vk(0))=λk(0)|vk(0)).\mathbb{G}_{\mathrm{WDA}}^{(0)}\left|v_{k}^{(0)}\right)=\lambda_{k}^{(0)}\left% |v_{k}^{(0)}\right).blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT | italic_v start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ) = italic_λ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT | italic_v start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ) . (59)

The smallest eigenvalue is denoted as λ0(0)superscriptsubscript𝜆00\lambda_{0}^{(0)}italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT. The correction to λ0(0)superscriptsubscript𝜆00\lambda_{0}^{(0)}italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT due to the high-order matrices 𝔾WDA(1)superscriptsubscript𝔾WDA1\mathbb{G}_{\rm WDA}^{(1)}blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT and 𝔾WDA(2)superscriptsubscript𝔾WDA2\mathbb{G}_{\rm WDA}^{(2)}blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT is calculated using perturbation theory. Up to the second order of P𝑃Pitalic_P, the minimum eigenvalue is

λ0=λ0(0)+λ0(1)P+(λ0(2,a)+λ0(2,b))P2,subscript𝜆0superscriptsubscript𝜆00superscriptsubscript𝜆01𝑃superscriptsubscript𝜆02𝑎superscriptsubscript𝜆02𝑏superscript𝑃2\lambda_{0}=\lambda_{0}^{(0)}+\lambda_{0}^{(1)}P+\left(\lambda_{0}^{(2,a)}+% \lambda_{0}^{(2,b)}\right)P^{2},italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT italic_P + ( italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 , italic_a ) end_POSTSUPERSCRIPT + italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 , italic_b ) end_POSTSUPERSCRIPT ) italic_P start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (60)

where

λ0(1)superscriptsubscript𝜆01\displaystyle\lambda_{0}^{(1)}italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT =(v0(0)|𝔾WDA(1)|v0(0)),absentsuperscriptsubscript𝑣00superscriptsubscript𝔾WDA1superscriptsubscript𝑣00\displaystyle=\left(v_{0}^{(0)}\right|\mathbb{G}_{\mathrm{WDA}}^{(1)}\left|v_{% 0}^{(0)}\right),= ( italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT | blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT | italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ) , (61)
λ0(2a)superscriptsubscript𝜆02𝑎\displaystyle\lambda_{0}^{(2a)}italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 italic_a ) end_POSTSUPERSCRIPT =(v0(0)|𝔾WDA(2)|v0(0)),absentsuperscriptsubscript𝑣00superscriptsubscript𝔾WDA2superscriptsubscript𝑣00\displaystyle=\left(v_{0}^{(0)}\right|\mathbb{G}_{\mathrm{WDA}}^{(2)}\left|v_{% 0}^{(0)}\right),= ( italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT | blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT | italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ) , (62)
λ0(2b)superscriptsubscript𝜆02𝑏\displaystyle\lambda_{0}^{(2b)}italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 italic_b ) end_POSTSUPERSCRIPT =k0(v0(0)|𝔾WDA(1)|vk(0))(vk(0)|𝔾WDA(1)|v0(0))λ0(0)λk(0).absentsubscript𝑘0superscriptsubscript𝑣00superscriptsubscript𝔾WDA1superscriptsubscript𝑣𝑘0superscriptsubscript𝑣𝑘0superscriptsubscript𝔾WDA1superscriptsubscript𝑣00superscriptsubscript𝜆00superscriptsubscript𝜆𝑘0\displaystyle=\sum_{k\neq 0}\frac{\left(v_{0}^{(0)}\right|\mathbb{G}_{\mathrm{% WDA}}^{(1)}\left|v_{k}^{(0)}\right)\left(v_{k}^{(0)}\right|\mathbb{G}_{\mathrm% {WDA}}^{(1)}\left|v_{0}^{(0)}\right)}{\lambda_{0}^{(0)}-\lambda_{k}^{(0)}}.= ∑ start_POSTSUBSCRIPT italic_k ≠ 0 end_POSTSUBSCRIPT divide start_ARG ( italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT | blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT | italic_v start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ) ( italic_v start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT | blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT | italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT - italic_λ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT end_ARG . (63)

For a given nuclear spin quantum number I𝐼Iitalic_I, in the low polarization limit, the leading order contribution of linewidth is

Γ2(0)=λ0(0)ΓSD=1qSDΓSD+1qSEΓSE,superscriptsubscriptΓ20superscriptsubscript𝜆00subscriptΓSD1subscript𝑞SDsubscriptΓSD1subscript𝑞SEsubscriptΓSE\displaystyle\Gamma_{2}^{(0)}=\lambda_{0}^{(0)}\Gamma_{\rm SD}=\frac{1}{q_{\rm SD% }}\Gamma_{\rm SD}+\frac{1}{q_{\rm SE}}\Gamma_{\rm SE},roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT = italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG italic_q start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG italic_q start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT , (64)

where two slowing-down factors qSDsubscript𝑞SDq_{\rm SD}italic_q start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT and qSEsubscript𝑞SEq_{\rm SE}italic_q start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT are defined as

qSDsubscript𝑞SD\displaystyle q_{\rm SD}italic_q start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT =(2I+1)22I2+I+1,absentsuperscript2𝐼122superscript𝐼2𝐼1\displaystyle=\frac{(2I+1)^{2}}{2I^{2}+I+1},= divide start_ARG ( 2 italic_I + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_I start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_I + 1 end_ARG , (65)
qSEsubscript𝑞SE\displaystyle q_{\rm SE}italic_q start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT =3(2I+1)22I(2I1).absent3superscript2𝐼122𝐼2𝐼1\displaystyle=\frac{3(2I+1)^{2}}{2I(2I-1)}.= divide start_ARG 3 ( 2 italic_I + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_I ( 2 italic_I - 1 ) end_ARG . (66)

The spin exchange slowing-down factor qSEsubscript𝑞SEq_{\rm SE}italic_q start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT in Eq. (66) is identical to that obtained in the spin-exchange relaxation free (SERF) regime, while the spin destruction slowing-down factor qSDsubscript𝑞SDq_{\rm SD}italic_q start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT in Eq. (65) is different from the longitudinal slowing-down factor at low polarization limit [14] and that in the resolved Zeeman resonance regime (see Sect. III.3).

The first order correction of the linewidth with P1much-less-than𝑃1P\ll 1italic_P ≪ 1 is

Γ2(1)=λ0(1)PΓSD=PqSDΓSD1qSDRop.superscriptsubscriptΓ21superscriptsubscript𝜆01𝑃subscriptΓSD𝑃subscript𝑞SDsubscriptΓSD1subscript𝑞SDsubscript𝑅op\Gamma_{2}^{(1)}=\lambda_{0}^{(1)}P\Gamma_{\rm SD}=\frac{P}{q_{\rm SD}}\Gamma_% {\mathrm{SD}}\approx\frac{1}{q_{\rm SD}}R_{\rm op}.roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT italic_P roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT = divide start_ARG italic_P end_ARG start_ARG italic_q start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT ≈ divide start_ARG 1 end_ARG start_ARG italic_q start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT . (67)

Note that the first order correction is proportional to the optical pumping rate Ropsubscript𝑅opR_{\rm op}italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT, and it is independent on the spin-exchange rate ΓSEsubscriptΓSE\Gamma_{\rm SE}roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT, which is different from the case of the RZR regime (see III.3). Furthermore, the correction in Eq. (67) is always positive, which implies that the linewidth is always broadened by optical pumping when RopΓSDmuch-less-thansubscript𝑅opsubscriptΓSDR_{\rm op}\ll\Gamma_{\rm SD}italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT ≪ roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT.

The light-narrowing effect manifests itself in the second order correction. Unfortunately, the analytic expression of the second order correction for the general nuclear spin quantum number I𝐼Iitalic_I is lengthy. Here, we show the second order correction for I=3/2𝐼32I=3/2italic_I = 3 / 2 (e.g., Rb87superscriptRb87{}^{87}{\rm Rb}start_FLOATSUPERSCRIPT 87 end_FLOATSUPERSCRIPT roman_Rb) as an example

Γ2(2)=λ0(2)P2ΓSD=7ΓSD234ΓSEΓSD98ΓSE2560ΓSEΓSD+160ΓSD2P2.superscriptsubscriptΓ22superscriptsubscript𝜆02superscript𝑃2subscriptΓSD7superscriptsubscriptΓSD234subscriptΓSEsubscriptΓSD98superscriptsubscriptΓSE2560subscriptΓSEsubscriptΓSD160superscriptsubscriptΓSD2superscript𝑃2\Gamma_{2}^{(2)}=\lambda_{0}^{(2)}P^{2}\Gamma_{\rm SD}=\frac{7\Gamma_{\mathrm{% SD}}^{2}-34\Gamma_{\mathrm{SE}}\Gamma_{\mathrm{SD}}-98\Gamma_{\mathrm{SE}}^{2}% }{560\Gamma_{\mathrm{SE}}\Gamma_{\mathrm{SD}}+160\Gamma_{\mathrm{SD}}^{2}}P^{2}.roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT = italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT italic_P start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT = divide start_ARG 7 roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 34 roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT - 98 roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 560 roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT + 160 roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_P start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (68)

The second order correction brings about a local maximum of the linewidth, as long as λ0(2)<0superscriptsubscript𝜆020\lambda_{0}^{(2)}<0italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT < 0. Figure 4(b) shows the linewidth in the low polarization regime, with a direct comparison of the the numerical result and the perturbation analysis in Eqs. (64)-(68). Figure 5(a) presents the position of the local maximum at various spin exchange rate, obtained from the numerical calculation. For the case of I=3/2𝐼32I=3/2italic_I = 3 / 2, the local maximum exists when the ratio η=ΓSE/ΓSD𝜂subscriptΓSEsubscriptΓSD\eta=\Gamma_{\rm SE}/\Gamma_{\rm SD}italic_η = roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT / roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT is above a threshold value η0=6.22subscript𝜂06.22\eta_{0}=6.22italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 6.22. This agrees with the analytic result in Eq. (68), where λ0(2)superscriptsubscript𝜆02\lambda_{0}^{(2)}italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT becomes positive if the spin exchange rate ΓSEsubscriptΓSE\Gamma_{\rm SE}roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT is too small. For weak spin exchange rate (e.g., in low density vapor) with η<η0𝜂subscript𝜂0\eta<\eta_{0}italic_η < italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, the local maximum does not exist, and the linewidth monotonically increases with the polarization as shown in Fig. 4(a).

III.2 High polarization limit (Q=1P1𝑄1𝑃much-less-than1Q=1-P\ll 1italic_Q = 1 - italic_P ≪ 1)

Refer to caption
Figure 5: Local minimum and local maximum of the linewidth of Rb87superscriptRb87{}^{87}{\rm Rb}start_FLOATSUPERSCRIPT 87 end_FLOATSUPERSCRIPT roman_Rb as a function of spin polarization P𝑃Pitalic_P. As increasing the spin-exchange rate, the position of the local maximum shifts to higher polarization, while the local minimum shifts in the opposite direction. The local minimum and local maximum merged together when the parameter η=η0𝜂subscript𝜂0\eta=\eta_{0}italic_η = italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. For η<η0𝜂subscript𝜂0\eta<\eta_{0}italic_η < italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, the linewidth increases monotonically as increasing the polarization.
Refer to caption
Figure 6: Light narrowing effect. Blue dots are the numerical results of the local minimum of the linewidth as a function of η𝜂\etaitalic_η. For η>η0𝜂subscript𝜂0\eta>\eta_{0}italic_η > italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, the linewidth of at the local minimum is smaller than that at P=0𝑃0P=0italic_P = 0 (the dashed curve) due to the light-narrowing effect. The red dashed-dotted curve is the analytic result of Eq. (81).

Similar analysis is performed in the high polarization limit, where the matrix 𝔾WDAsubscript𝔾WDA\mathbb{G}_{\rm{WDA}}blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT is expanded as

𝔾WDA(η,1Q)=Q1𝔾~WDA(1)+𝔾~WDA(0)+Q𝔾~WDA(1),subscript𝔾WDA𝜂1𝑄superscript𝑄1superscriptsubscript~𝔾WDA1superscriptsubscript~𝔾WDA0𝑄superscriptsubscript~𝔾WDA1\mathbb{G}_{\mathrm{WDA}}(\eta,1-Q)=Q^{-1}\cdot\tilde{\mathbb{G}}_{\mathrm{WDA% }}^{(-1)}+\tilde{\mathbb{G}}_{\mathrm{WDA}}^{(0)}+Q\cdot\tilde{\mathbb{G}}_{% \mathrm{WDA}}^{(1)},blackboard_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT ( italic_η , 1 - italic_Q ) = italic_Q start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ⋅ over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT + over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + italic_Q ⋅ over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT , (69)

where

𝔾~WDA(1)=superscriptsubscript~𝔾WDA1absent\displaystyle\tilde{\mathbb{G}}_{\mathrm{WDA}}^{(-1)}=over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT = [𝒜SD]1,112[𝒜SE(z)]1,1,subscriptdelimited-[]subscript𝒜SD1112subscriptdelimited-[]superscriptsubscript𝒜SE𝑧11\displaystyle\left[\mathcal{A}_{\mathrm{SD}}\right]_{1,1}-\frac{1}{2}\left[% \mathcal{A}_{\mathrm{SE}}^{(z)}\right]_{1,1},[ caligraphic_A start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_z ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT , (70)
𝔾~WDA(0)=superscriptsubscript~𝔾WDA0absent\displaystyle\tilde{\mathbb{G}}_{\mathrm{WDA}}^{(0)}=over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT = η[𝒜SD]1,1η12[𝒜SE(z)]1,1η[𝒜~SE()]1,1(0),𝜂subscriptdelimited-[]subscript𝒜SD11𝜂12subscriptdelimited-[]superscriptsubscript𝒜SE𝑧11𝜂superscriptsubscriptdelimited-[]superscriptsubscript~𝒜SEperpendicular-to110\displaystyle\eta\left[\mathcal{A}_{\mathrm{SD}}\right]_{1,1}-\frac{\eta-1}{2}% \left[\mathcal{A}_{\mathrm{SE}}^{(z)}\right]_{1,1}-\eta\left[\tilde{\mathcal{A% }}_{\mathrm{SE}}^{(\perp)}\right]_{1,1}^{(0)},italic_η [ caligraphic_A start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - divide start_ARG italic_η - 1 end_ARG start_ARG 2 end_ARG [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_z ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - italic_η [ over~ start_ARG caligraphic_A end_ARG start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT , (71)
𝔾~WDA(1)=superscriptsubscript~𝔾WDA1absent\displaystyle\tilde{\mathbb{G}}_{\mathrm{WDA}}^{(1)}=over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = η2[𝒜SE(z)]1,1η[𝒜~SE()]1,1(1).𝜂2subscriptdelimited-[]superscriptsubscript𝒜SE𝑧11𝜂superscriptsubscriptdelimited-[]superscriptsubscript~𝒜SEperpendicular-to111\displaystyle\frac{\eta}{2}\left[\mathcal{A}_{\mathrm{SE}}^{(z)}\right]_{1,1}-% \eta\left[\tilde{\mathcal{A}}_{\mathrm{SE}}^{(\perp)}\right]_{1,1}^{(1)}.divide start_ARG italic_η end_ARG start_ARG 2 end_ARG [ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_z ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT - italic_η [ over~ start_ARG caligraphic_A end_ARG start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT . (72)

In Eqs. (71) & (72), we have used the expansion of the transverse spin exchange matrix

[𝒜SE()(1Q)]1,1=[𝒜~SE()]1,1(0)+Q[𝒜~SE()]1,1(1).subscriptdelimited-[]superscriptsubscript𝒜SEperpendicular-to1𝑄11superscriptsubscriptdelimited-[]superscriptsubscript~𝒜SEperpendicular-to110𝑄superscriptsubscriptdelimited-[]superscriptsubscript~𝒜SEperpendicular-to111\left[\mathcal{A}_{\mathrm{SE}}^{(\perp)}(1-Q)\right]_{1,1}=\left[\tilde{% \mathcal{A}}_{\mathrm{SE}}^{(\perp)}\right]_{1,1}^{(0)}+Q\cdot\left[\tilde{% \mathcal{A}}_{\mathrm{SE}}^{(\perp)}\right]_{1,1}^{(1)}.[ caligraphic_A start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ( 1 - italic_Q ) ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT = [ over~ start_ARG caligraphic_A end_ARG start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + italic_Q ⋅ [ over~ start_ARG caligraphic_A end_ARG start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( ⟂ ) end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT 1 , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT . (73)

Note that the matrix 𝔾~WDA(1)superscriptsubscript~𝔾WDA1\tilde{\mathbb{G}}_{\mathrm{WDA}}^{(-1)}over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT is not Hermitian. With the eigenvalues λ~k(1)superscriptsubscript~𝜆𝑘1\tilde{\lambda}_{k}^{(-1)}over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT and the corresponding right- and left-eigenvectors |v~k(1))\left|\tilde{v}_{k}^{(-1)}\right)| over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT ) and {v~k(1)|\left\{\tilde{v}_{k}^{(-1)}\right|{ over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT | of the leading order matrix 𝔾~WDA(1)superscriptsubscript~𝔾WDA1\tilde{\mathbb{G}}_{\mathrm{WDA}}^{(-1)}over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT, which satisfy the following equations

𝔾~WDA(1)|v~k(1))\displaystyle\tilde{\mathbb{G}}_{\mathrm{WDA}}^{(-1)}\left|\tilde{v}_{k}^{(-1)% }\right)over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT | over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT ) =λ~k(1)|v~k(1)),\displaystyle=\tilde{\lambda}_{k}^{(-1)}\left|\tilde{v}_{k}^{(-1)}\right),= over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT | over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT ) , (74)
{v~k(1)|𝔾~WDA(1)\displaystyle\left\{\tilde{v}_{k}^{(-1)}\right|\tilde{\mathbb{G}}_{\mathrm{WDA% }}^{(-1)}{ over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT | over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT =λ~k(1){v~k(1)|,\displaystyle=\tilde{\lambda}_{k}^{(-1)}\left\{\tilde{v}_{k}^{(-1)}\right|,= over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT { over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT | , (75)

the minimum eigenvalues expanded in a power series of Q𝑄Qitalic_Q is

λ~0=Q1λ~0(1)+λ~0(0)+(λ~0(1,a)+λ~0(1,b))Q,subscript~𝜆0superscript𝑄1superscriptsubscript~𝜆01superscriptsubscript~𝜆00superscriptsubscript~𝜆01𝑎superscriptsubscript~𝜆01𝑏𝑄\tilde{\lambda}_{0}=Q^{-1}\tilde{\lambda}_{0}^{(-1)}+\tilde{\lambda}_{0}^{(0)}% +\left(\tilde{\lambda}_{0}^{(1,a)}+\tilde{\lambda}_{0}^{(1,b)}\right)Q,over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_Q start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT + over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + ( over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 , italic_a ) end_POSTSUPERSCRIPT + over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 , italic_b ) end_POSTSUPERSCRIPT ) italic_Q , (76)

where

λ~0(0)superscriptsubscript~𝜆00\displaystyle\tilde{\lambda}_{0}^{(0)}over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ={v~0(1)|𝔾~WDA(0)|v~0(1)),\displaystyle=\left\{\tilde{v}_{0}^{(-1)}\right|\tilde{\mathbb{G}}_{\mathrm{% WDA}}^{(0)}\left|\tilde{v}_{0}^{(-1)}\right),= { over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT | over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT | over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT ) , (77)
λ~0(1,a)superscriptsubscript~𝜆01𝑎\displaystyle\tilde{\lambda}_{0}^{(1,a)}over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 , italic_a ) end_POSTSUPERSCRIPT ={v~0(1)|𝔾~WDA(1)|v~0(1)),\displaystyle=\left\{\tilde{v}_{0}^{(-1)}\right|\tilde{\mathbb{G}}_{\mathrm{% WDA}}^{(1)}\left|\tilde{v}_{0}^{(-1)}\right),= { over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT | over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT | over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT ) , (78)
λ~0(1,b)superscriptsubscript~𝜆01𝑏\displaystyle\tilde{\lambda}_{0}^{(1,b)}over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 , italic_b ) end_POSTSUPERSCRIPT =k0{v~0(1)|𝔾~WDA(0)|v~k(1)){v~k(1)|𝔾~WDA(0)|v~0(1))λ~0(1)λ~k(1)\displaystyle=\sum_{k\neq 0}\frac{\left\{\tilde{v}_{0}^{(-1)}\right|\tilde{% \mathbb{G}}_{\mathrm{WDA}}^{(0)}\left|\tilde{v}_{k}^{(-1)}\right)\left\{\tilde% {v}_{k}^{(-1)}\right|\tilde{\mathbb{G}}_{\mathrm{WDA}}^{(0)}\left|\tilde{v}_{0% }^{(-1)}\right)}{\tilde{\lambda}_{0}^{(-1)}-\tilde{\lambda}_{k}^{(-1)}}= ∑ start_POSTSUBSCRIPT italic_k ≠ 0 end_POSTSUBSCRIPT divide start_ARG { over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT | over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT | over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT ) { over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT | over~ start_ARG blackboard_G end_ARG start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT | over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT ) end_ARG start_ARG over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT - over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( - 1 ) end_POSTSUPERSCRIPT end_ARG (79)

With Eqs. (77)-(79), the linewidth in the high polarization limit is

Γ2=λ~0ΓSD=12I+1Rop+12ΓSD+1ξIΓSDΓSERop.subscriptΓ2subscript~𝜆0subscriptΓSD12𝐼1subscript𝑅op12subscriptΓSD1subscript𝜉𝐼subscriptΓSDsubscriptΓSEsubscript𝑅op\Gamma_{2}=\tilde{\lambda}_{0}\Gamma_{\rm SD}=\frac{1}{2I+1}R_{\rm op}+\frac{1% }{2}\Gamma_{\rm SD}+\frac{1}{\xi_{I}}\frac{\Gamma_{\rm SD}\Gamma_{\rm SE}}{R_{% \rm op}}.roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 italic_I + 1 end_ARG italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG italic_ξ start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT end_ARG divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG start_ARG italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT end_ARG . (80)

where ξIsubscript𝜉𝐼\xi_{I}italic_ξ start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT is factor depending on the nuclear spin quantum number I𝐼Iitalic_I. The complete expressions of linewidth for different isotopes are listed in Table 1. Although the coefficient in front of ΓSEsubscriptΓSE\Gamma_{\rm SE}roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT can be small in the high polarization limit with RopΓSDmuch-greater-thansubscript𝑅opsubscriptΓSDR_{\rm op}\gg\Gamma_{\rm SD}italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT ≫ roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT, the third term in Eq. (80) is not negligible if the spin exchange rate ΓSEsubscriptΓSE\Gamma_{\rm SE}roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT is large.

I𝐼\displaystyle Iitalic_I isotope Γ2(P0)subscriptΓ2𝑃0\displaystyle\Gamma_{2}(P\to 0)roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_P → 0 ) Γ2(P1)subscriptΓ2𝑃1\displaystyle\Gamma_{2}(P\to 1)roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_P → 1 ) η0subscript𝜂0\eta_{0}italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT
1111 Li6superscriptLi6{}^{6}{\rm Li}start_FLOATSUPERSCRIPT 6 end_FLOATSUPERSCRIPT roman_Li Γ2=227ΓSE+49ΓSDsubscriptΓ2227subscriptΓSE49subscriptΓSD\displaystyle\Gamma_{2}=\frac{2}{27}\Gamma_{\text{SE}}+\frac{4}{9}\Gamma_{% \text{SD}}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG 2 end_ARG start_ARG 27 end_ARG roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT + divide start_ARG 4 end_ARG start_ARG 9 end_ARG roman_Γ start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT Γ2=ROP3+ΓSD2+ΓSDΓSE8RopsubscriptΓ2subscript𝑅OP3subscriptΓSD2subscriptΓSDsubscriptΓSE8subscript𝑅op\displaystyle\Gamma_{2}=\frac{R_{\mathrm{OP}}}{3}+\frac{\Gamma_{\mathrm{SD}}}{% 2}+\frac{\Gamma_{\rm SD}\Gamma_{\rm SE}}{8R_{\rm op}}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG italic_R start_POSTSUBSCRIPT roman_OP end_POSTSUBSCRIPT end_ARG start_ARG 3 end_ARG + divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG + divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG start_ARG 8 italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT end_ARG 13.82
3232\displaystyle\frac{3}{2}divide start_ARG 3 end_ARG start_ARG 2 end_ARG Li7superscriptLi7{}^{7}{\rm Li}start_FLOATSUPERSCRIPT 7 end_FLOATSUPERSCRIPT roman_Li, Na23superscriptNa23{}^{23}{\rm Na}start_FLOATSUPERSCRIPT 23 end_FLOATSUPERSCRIPT roman_Na, K39superscriptK39{}^{39}{\rm K}start_FLOATSUPERSCRIPT 39 end_FLOATSUPERSCRIPT roman_K, K41superscriptK41{}^{41}{\rm K}start_FLOATSUPERSCRIPT 41 end_FLOATSUPERSCRIPT roman_K, Rb87superscriptRb87{}^{87}{\rm Rb}start_FLOATSUPERSCRIPT 87 end_FLOATSUPERSCRIPT roman_Rb Γ2=18ΓSE+716ΓSDsubscriptΓ218subscriptΓSE716subscriptΓSD\displaystyle\Gamma_{2}=\frac{1}{8}\Gamma_{\text{SE}}+\frac{7}{16}\Gamma_{% \text{SD}}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 8 end_ARG roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT + divide start_ARG 7 end_ARG start_ARG 16 end_ARG roman_Γ start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT Γ2=ROP4+ΓSD2+5ΓSDΓSE24RopsubscriptΓ2subscript𝑅OP4subscriptΓSD25subscriptΓSDsubscriptΓSE24subscript𝑅op\displaystyle\Gamma_{2}=\frac{R_{\mathrm{OP}}}{4}+\frac{\Gamma_{\mathrm{SD}}}{% 2}+\frac{5\Gamma_{\rm SD}\Gamma_{\rm SE}}{24R_{\rm op}}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG italic_R start_POSTSUBSCRIPT roman_OP end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG + divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG + divide start_ARG 5 roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG start_ARG 24 italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT end_ARG 6.22
5252\displaystyle\frac{5}{2}divide start_ARG 5 end_ARG start_ARG 2 end_ARG Rb85superscriptRb85{}^{85}{\rm Rb}start_FLOATSUPERSCRIPT 85 end_FLOATSUPERSCRIPT roman_Rb Γ2=527ΓSE+49ΓSDsubscriptΓ2527subscriptΓSE49subscriptΓSD\displaystyle\Gamma_{2}=\frac{5}{27}\Gamma_{\text{SE}}+\frac{4}{9}\Gamma_{% \text{SD}}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG 5 end_ARG start_ARG 27 end_ARG roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT + divide start_ARG 4 end_ARG start_ARG 9 end_ARG roman_Γ start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT Γ2=ROP6+ΓSD2+3ΓSDΓSE10RopsubscriptΓ2subscript𝑅OP6subscriptΓSD23subscriptΓSDsubscriptΓSE10subscript𝑅op\displaystyle\Gamma_{2}=\frac{R_{\mathrm{OP}}}{6}+\frac{\Gamma_{\mathrm{SD}}}{% 2}+\frac{3\Gamma_{\rm SD}\Gamma_{\rm SE}}{10R_{\rm op}}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG italic_R start_POSTSUBSCRIPT roman_OP end_POSTSUBSCRIPT end_ARG start_ARG 6 end_ARG + divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG + divide start_ARG 3 roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG start_ARG 10 italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT end_ARG 2.77
7272\displaystyle\frac{7}{2}divide start_ARG 7 end_ARG start_ARG 2 end_ARG Cs133superscriptCs133{}^{133}{\rm Cs}start_FLOATSUPERSCRIPT 133 end_FLOATSUPERSCRIPT roman_Cs Γ2=732ΓSE+2964ΓSDsubscriptΓ2732subscriptΓSE2964subscriptΓSD\displaystyle\Gamma_{2}=\frac{7}{32}\Gamma_{\text{SE}}+\frac{29}{64}\Gamma_{% \text{SD}}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG 7 end_ARG start_ARG 32 end_ARG roman_Γ start_POSTSUBSCRIPT SE end_POSTSUBSCRIPT + divide start_ARG 29 end_ARG start_ARG 64 end_ARG roman_Γ start_POSTSUBSCRIPT SD end_POSTSUBSCRIPT Γ2=ROP8+ΓSD2+39ΓSDΓSE112RopsubscriptΓ2subscript𝑅OP8subscriptΓSD239subscriptΓSDsubscriptΓSE112subscript𝑅op\displaystyle\Gamma_{2}=\frac{R_{\mathrm{OP}}}{8}+\frac{\Gamma_{\mathrm{SD}}}{% 2}+\frac{39\Gamma_{\rm SD}\Gamma_{\rm SE}}{112R_{\rm op}}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG italic_R start_POSTSUBSCRIPT roman_OP end_POSTSUBSCRIPT end_ARG start_ARG 8 end_ARG + divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG + divide start_ARG 39 roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG start_ARG 112 italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT end_ARG 1.66
Table 1: Linewidth of different isotopes in low and high polarization limits. The analytic results of the linewidth Γ2subscriptΓ2\Gamma_{2}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are kept to the lowest order of polarization P𝑃Pitalic_P or 1P1𝑃1-P1 - italic_P. See the text for higher order corrections.

Equation (80) implies that the linewidth has a local minimum when changing the optical pumping rate Ropsubscript𝑅opR_{\rm op}italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT. Figure 5 shows the position of the local minimum for different values of the spin exchange rate ΓSEsubscriptΓSE\Gamma_{\rm SE}roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT. At high spin exchange rate (i.e., η1much-greater-than𝜂1\eta\gg 1italic_η ≫ 1), the local minimum and local maximum are well separated. As the spin exchange rate decreases, the local minimum and maximum of the linewidth are merged into a single point, denoted as (η0,P0)subscript𝜂0subscript𝑃0\left(\eta_{0},P_{0}\right)( italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_P start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ). For weak spin exchange rate with η<η0𝜂subscript𝜂0\eta<\eta_{0}italic_η < italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, the line width monotonically increases with increasing optical pumping rate Ropsubscript𝑅opR_{\rm op}italic_R start_POSTSUBSCRIPT roman_op end_POSTSUBSCRIPT. To observe the light-narrowing effect, the spin exchange rate must be larger than the threshold value, that is, η>η0𝜂subscript𝜂0\eta>\eta_{0}italic_η > italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The threshold values η0subscript𝜂0\eta_{0}italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT for different isotopes are listed in Table 1.

Figure 6 compares the linewidth in the low polarization limit with the local minimum. Equation (64) shows that the linewidth scales linearly with the spin-exchange rate at P=0𝑃0P=0italic_P = 0, while, according to Eq. (80), the linewidth at the local minimum scales as ΓSE1/2similar-toabsentsuperscriptsubscriptΓSE12\sim\Gamma_{\rm SE}^{1/2}∼ roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT, i.e.,

Γ2,min=2ΓSDΓSE(2I+1)ξI,subscriptΓ2min2subscriptΓSDsubscriptΓSE2𝐼1subscript𝜉𝐼\Gamma_{2,{\rm min}}=2\sqrt{\frac{\Gamma_{\rm SD}\Gamma_{\rm SE}}{(2I+1)\xi_{I% }}},roman_Γ start_POSTSUBSCRIPT 2 , roman_min end_POSTSUBSCRIPT = 2 square-root start_ARG divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG start_ARG ( 2 italic_I + 1 ) italic_ξ start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT end_ARG end_ARG , (81)

which agrees well with the numerical results (see Fig. 6).

III.3 Comparison to the RZR regime

The linewidth in the RZR regime is [8, 7, 11]

Γ2,m¯RZRsuperscriptsubscriptΓ2¯𝑚RZR\displaystyle\Gamma_{2,\bar{m}}^{\text{RZR}}roman_Γ start_POSTSUBSCRIPT 2 , over¯ start_ARG italic_m end_ARG end_POSTSUBSCRIPT start_POSTSUPERSCRIPT RZR end_POSTSUPERSCRIPT =(ΓSE+ΓSD+ROP)3(2I+1)24m¯2+14(2I+1)2absentsubscriptΓSEsubscriptΓSDsubscript𝑅OP3superscript2𝐼124superscript¯𝑚214superscript2𝐼12\displaystyle=\left(\Gamma_{\mathrm{SE}}+\Gamma_{\mathrm{SD}}+R_{\mathrm{OP}}% \right)\frac{3(2I+1)^{2}-4\bar{m}^{2}+1}{4(2I+1)^{2}}= ( roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT + roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT + italic_R start_POSTSUBSCRIPT roman_OP end_POSTSUBSCRIPT ) divide start_ARG 3 ( 2 italic_I + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 over¯ start_ARG italic_m end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 end_ARG start_ARG 4 ( 2 italic_I + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG
(PΓSE+ROP)m¯2(FI)(2I+1)ΓSEeβm¯ZI(2F+1)24m¯24(2I+1)2,𝑃subscriptΓSEsubscript𝑅OP¯𝑚2𝐹𝐼2𝐼1subscriptΓSEsuperscript𝑒𝛽¯𝑚subscript𝑍𝐼superscript2𝐹124superscript¯𝑚24superscript2𝐼12\displaystyle-\frac{(P\Gamma_{\mathrm{SE}}+R_{\mathrm{OP}})\bar{m}}{2(F-I)(2I+% 1)}-\Gamma_{\mathrm{SE}}\frac{e^{\beta\bar{m}}}{Z_{I}}\frac{(2F+1)^{2}-4\bar{m% }^{2}}{4(2I+1)^{2}},- divide start_ARG ( italic_P roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT + italic_R start_POSTSUBSCRIPT roman_OP end_POSTSUBSCRIPT ) over¯ start_ARG italic_m end_ARG end_ARG start_ARG 2 ( italic_F - italic_I ) ( 2 italic_I + 1 ) end_ARG - roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT divide start_ARG italic_e start_POSTSUPERSCRIPT italic_β over¯ start_ARG italic_m end_ARG end_POSTSUPERSCRIPT end_ARG start_ARG italic_Z start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT end_ARG divide start_ARG ( 2 italic_F + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 over¯ start_ARG italic_m end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 ( 2 italic_I + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (82)

where F=a𝐹𝑎F=aitalic_F = italic_a or F=b𝐹𝑏F=bitalic_F = italic_b. Indeed, for a given Zeeman resonance labeled with m¯¯𝑚\bar{m}over¯ start_ARG italic_m end_ARG, the linewidth in the RZR regime is the diagonal element of the matrix 𝒢WDAsubscript𝒢WDA\mathcal{G}_{\rm WDA}caligraphic_G start_POSTSUBSCRIPT roman_WDA end_POSTSUBSCRIPT in Eq. (II.4.1). A direct comparison between the linewidth in the UZR and RZR regimes for the I=3/2𝐼32I=3/2italic_I = 3 / 2 case is shown in Fig. 7. Here, we focus on the narrowest Zeeman resonance, i.e., m¯=I¯𝑚𝐼\bar{m}=Iover¯ start_ARG italic_m end_ARG = italic_I.

Refer to caption
Figure 7: Comparison the linewidth in the UZR and RZR regimes for the I=3/2𝐼32I=3/2italic_I = 3 / 2 case. In the low polarization limit, the linewidth in the RZR regime is larger than that in the UZR regime. In the high polarization limit, the two linewidths converge to approximately the same value.

In the low polarization limit P1much-less-than𝑃1P\ll 1italic_P ≪ 1, the linewidth in the RZR regime is

Γ2,IRZR=ΓSDqSD+ΓSEqSE+I(2I+1)ΓSD(4I+1)ΓSE(2I+1)2P,superscriptsubscriptΓ2𝐼RZRsubscriptΓSDsuperscriptsubscript𝑞SDsubscriptΓSEsuperscriptsubscript𝑞SE𝐼2𝐼1subscriptΓSD4𝐼1subscriptΓSEsuperscript2𝐼12𝑃\Gamma_{2,I}^{\text{RZR}}=\frac{\Gamma_{\mathrm{SD}}}{q_{\rm SD}^{\prime}}+% \frac{\Gamma_{\mathrm{SE}}}{q_{\rm SE}^{\prime}}+\frac{I(2I+1)\Gamma_{\rm SD}-% (4I+1)\Gamma_{\rm SE}}{(2I+1)^{2}}P,roman_Γ start_POSTSUBSCRIPT 2 , italic_I end_POSTSUBSCRIPT start_POSTSUPERSCRIPT RZR end_POSTSUPERSCRIPT = divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG start_ARG italic_q start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG + divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG start_ARG italic_q start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG + divide start_ARG italic_I ( 2 italic_I + 1 ) roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT - ( 4 italic_I + 1 ) roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG start_ARG ( 2 italic_I + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_P , (83)

where the slowing-down factors are defined as

qSDsuperscriptsubscript𝑞SD\displaystyle q_{\rm SD}^{\prime}italic_q start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT =2I+1I+1,absent2𝐼1𝐼1\displaystyle=\frac{2I+1}{I+1},= divide start_ARG 2 italic_I + 1 end_ARG start_ARG italic_I + 1 end_ARG , (84)
qSEsuperscriptsubscript𝑞SE\displaystyle q_{\rm SE}^{\prime}italic_q start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT =(2I+1)2I(2I+3).absentsuperscript2𝐼12𝐼2𝐼3\displaystyle=\frac{(2I+1)^{2}}{I(2I+3)}.= divide start_ARG ( 2 italic_I + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_I ( 2 italic_I + 3 ) end_ARG . (85)

The smaller slowing-down factors in Eq. (84) & (85) compared to Eq. (65) & (66) causes larger linewidth in the RZR regime than in the UZR regime. Furthermore, the sign of the coefficient in front of P𝑃Pitalic_P depends on the ratio η=ΓSE/ΓSD𝜂subscriptΓSEsubscriptΓSD\eta=\Gamma_{\rm SE}/\Gamma_{\rm SD}italic_η = roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT / roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT. As increasing ΓSEsubscriptΓSE\Gamma_{\rm SE}roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT, the coefficient become negative and the light-narrowing effect occurs when

ΓSEΓSD>I(2I+1)4I+1.subscriptΓSEsubscriptΓSD𝐼2𝐼14𝐼1\frac{\Gamma_{\rm SE}}{\Gamma_{\rm SD}}>\frac{I(2I+1)}{4I+1}.divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT end_ARG > divide start_ARG italic_I ( 2 italic_I + 1 ) end_ARG start_ARG 4 italic_I + 1 end_ARG . (86)

This is different from the UZR regime, where the first order correction in Eq. (67) is always positive, and the second order correction in Eq. (68) drives the light-narrowing effect.

In the high polarization limit, the linewidth in the RZR regime is

Γ2,IRZR=12I+1ROP+I+12I+1ΓSD+I2I+1ΓSDΓSEROP.superscriptsubscriptΓ2𝐼RZR12𝐼1subscript𝑅OP𝐼12𝐼1subscriptΓSD𝐼2𝐼1subscriptΓSDsubscriptΓSEsubscript𝑅OP\Gamma_{2,I}^{\rm RZR}=\frac{1}{2I+1}R_{\mathrm{OP}}+\frac{I+1}{2I+1}\Gamma_{% \mathrm{SD}}+\frac{I}{2I+1}\frac{\Gamma_{\mathrm{SD}}\Gamma_{\mathrm{SE}}}{R_{% \mathrm{OP}}}.roman_Γ start_POSTSUBSCRIPT 2 , italic_I end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_RZR end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 italic_I + 1 end_ARG italic_R start_POSTSUBSCRIPT roman_OP end_POSTSUBSCRIPT + divide start_ARG italic_I + 1 end_ARG start_ARG 2 italic_I + 1 end_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT + divide start_ARG italic_I end_ARG start_ARG 2 italic_I + 1 end_ARG divide start_ARG roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT end_ARG start_ARG italic_R start_POSTSUBSCRIPT roman_OP end_POSTSUBSCRIPT end_ARG . (87)

As the linewidth is dominated by strong optical pumping, the linewidth in the RZR regime [see Eq. (87)] is approximately the same as that in the UZR regime [see Eq. (80)], as shown in Fig. 7. Physically, this result is reasonable because, in the high polarization limit, the Zeeman resonances with m¯<I¯𝑚𝐼\bar{m}<Iover¯ start_ARG italic_m end_ARG < italic_I are less populated, and the interplay between different Zeeman resonances becomes less important.

IV Conclusions

We studied the linewidth of the magnetic resonance of alkali-metal atoms in the UZR regime. A theoretical framework is developed to describe the dynamics of the Zeeman coherences and the populations, in presence of the optical pumping, spin exchange, and spin destruction processes. Under the RWA and WDA, the master equation is solved, and important phenomena such as the light-narrowing effect and the light-broadening effect are discussed. We obtain the linewidth of alkali metal atoms with different nuclear spin I𝐼Iitalic_I, as functions of the spin polarization parameter P𝑃Pitalic_P and the relative spin-exchange strength η𝜂\etaitalic_η. Analytic analysis of the linewidth in the UZR regime is provided in the low- and high-polarization limits. Based on the analytic results, a direct connection between the microscopic rates (e.g. the spin exchange rate ΓSEsubscriptΓSE\Gamma_{\rm SE}roman_Γ start_POSTSUBSCRIPT roman_SE end_POSTSUBSCRIPT and the spin destruction rate ΓSDsubscriptΓSD\Gamma_{\rm SD}roman_Γ start_POSTSUBSCRIPT roman_SD end_POSTSUBSCRIPT) and the experimentally observable resonance linewidth is established. We show that the linewidth in the UZR regime is different from previous results in the RZR regime, due to the interplay between different Zeeman resonances. Our study provides a deep understanding and useful theoretical tools for developing various atomic sensors (e.g. atomic magnetometers or comagnetometers) working in the UZR regime.

References