Magnetic Resonance Linewidth of Alkali-Metal Vapor in Unresolved Zeeman Resonance Regime
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 considerably exceeds the Larmor precession frequency of spin rotation (i.e., ), 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
II.1 Master equation in Liouville space
We study alkali metal spins in a static magnetic field in the direction and a radio-frequency (RF) driving field in the direction. An circularly polarized laser beam travels along the 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 of the alkali-metal atom is [8]
(1) |
The coherent evolution is determined by the Hamiltonian
(2) |
where
(3) | ||||
(4) |
represent the atomic spin Hamiltonian within the magnetic field , along with the RF driving field. In Eq. (3), and correspond to the electronic spin and nuclear spin operators of the alkali metal atom. Here, denotes the Larmor frequency ( being the gyromagnetic ratio), and is the hyperfine interaction constant. Throughout this paper, it is assumed that the magnetic field is sufficiently weak so that holds true. In Eq. (4), , and stand for the Rabi frequency, the frequency and the phase of the RF driving field, respectively. For simplicity, 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 , where 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 . 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 . The state
(5) |
is a purely nuclear spin operator of the alkali-metal atom, without electronic spin polarization [12].
In Liouville space, the density matrix is mapped to a state vector . Thee state vector is expanded as
(6) |
where
(7) |
is the basis vector of the Liouville space, and is the eigen vector of the alkali-metal Hamiltonian (3), with the total angular momentum quantum number or , and the magnetic quantum number . For the state vectors constrained in the subspace of population and Zeeman coherences with , we define
(8) |
where and [8].
The master Eq. (II.1) in Liouville space becomes
(9) |
where is the spin expectation value expressed in Liouville space, and , and 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]
(10) | ||||
(11) | ||||
(12) |
where and are the left- and right-translation superoperators, corresponding to the operator in Hilbert space, and is the commutator superoperator [13].
II.2 Rotating wave approximation
We consider the near-resonant RF driving field with frequency . The driving field only induces coherences between the Zeeman sublevels within the or subspaces. In the absence of microwave excitations, the hyperfine coherences between levels and (e.g., ) will decay within the time scale of the hyperfine relaxation time, typically in the order of 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
(13) | ||||
(14) | ||||
(15) |
and the state vector projected to the population and Zeeman coherence subspace is
(16) |
The evolution of is governed by
(17) |
where
(18) |
where , and represents the projection of the superoperator to the subspace of population and Zeeman coherences. As an example, Fig. 1 shows the matrix of the superoperator on the basis of the population and Zeeman coherences.
To further simplify the master equation, we define the rotating frame with the frequency and make the rotating wave approximation (RWA). We introduce the th order Zeeman coherence projector.
(19) |
and the free evolution superoperator
(20) |
With the rotation transform generated by , i.e.,
(21) |
the master equation of the slow-varying state vector in the rotating frame is
(22) |
The evolution superoperator is
(23) |
where is a diagonal matrix representing the frequency detuning of the th order coherences ( is an identity matrix). The driving superoperator in the rotating frame is
(24) |
and the relaxation superoperator is
(25) |
In Eqs. (23) - (II.2), is the block matrix of the th row and the th column, corresponding to a time-dependent factor in the rotating frame. The spin expectation value in the rotating frame is
(26) |
where
(27) |
is the spin expectation value on the th order Zeeman coherence .
Keeping only the zero-frequency terms of , we obtain the RWA equation
(28) |
where the evolution superoperator can be decomposed as with
(29) |
In Eq. (II.2), the spin exchange due to the transverse and longitudinal spin components are
(30) |
and
(31) |
II.3 Resolved and unresolved Zeeman resonance regimes
With the RWA, the evolution superoperator 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 atoms in different magnetic fields.
In a weak field, the frequency difference between the th order Zeeman coherences is negligible, and the detuning matrix in Eq. (II.2) is proportional to an identity matrix, i.e., . 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 and is lifted because of the slight difference of their gyromagnetic ratios. Taking atoms as an example, the difference between the gyromagnetic ratios of and levels is , which causes a splitting of in a field of [see Fig. 2(b)]. In this field, the Zeeman coherences within the or subspaces can still be regarded as degenerate, and the spectrum consists of two peaks.
II.4 Weak driving approximation
II.4.1 Weak driving approximation in UZR regime
With the operator applied on both sides of Eqs. (28), the master equation satisfied by the first order Zeeman coherence is
(32) |
where the diagonal block is
(33) |
and the off-diagonal blocks
(34) | ||||
(35) |
couple the first order Zeeman coherence to the population and the second order Zeeman coherence , due to the RF driving field and the transverse components of the spin exchange. In Eq. (33), we have replaced the detuning matrix by the detuning frequency (the identity matrix is omitted for simplicity) in the UZR regime.
In a weak RF driving field, the th order Zeeman coherence is on the order of [8], where is the resonance linewidth. Although the exact value of is unknown before solving the master equation, it is reasonable that is in the same order as the relaxation rates , or . We can always control the strength of the RF driving field, so that is well satisfied. In this case, the population is approximately constant, and we replace in Eq. (32) by the equilibrium state in the absence of RF driving (with optical pumping and spin relaxation only), i.e.,
(36) |
Furthermore, the second order Zeeman coherence in Eq. (32) is neglected. With this weak driving approximation (WDA) Eq. (32) becomes
(37) |
In Eq. (II.4.1), the superoperator due to the spin exchange of transverse components is
(38) |
where and are matrices associated with the equilibrium state as
(39) |
The inhomogeneous term of Eq. (II.4.1) arises from the RF driving field, with the state vector given by
(40) |
Notice that, with the WDA, Eq. (II.4.1) is a closed linear equation for 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 . 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
(41) |
where , and is the normalization factor. The parameter is related to the electron spin polarization by
(42) |
The probability of occupation of the spin state is a function of spin polarization
(43) |
where . As the state is determined by , the matrices defined in Eq. (39), together with the superoperator defined in Eq. (38), are also functions of .
With the spin temperature distribution, the dimensionless WDA superoperator (in the unit of ) becomes
(44) | ||||
where . In Eq. (44), we have used the fact that the spin polarization of the spin temperature state is
(45) |
or, equivalently, the dimensionless pumping rate is expressed in terms of as
(46) |
Notice that the dimensionless WDA superoperator is fully characterized by two dimensionless parameters, namely, the relative strength of spin exchange and the spin polarization .
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. , via the Faraday rotation effect. According to the discussion above, the time-dependent observable is
(47) |
and the in-phase and out-of-phase quadratures are
(48) | ||||
(49) |
Assume the right- and left-eigenstates of the superoperator , sharing the same eigenvalue , are denoted as and , respectively,
(50) | |||
(51) |
The observable is decomposed as
(52) |
where the weight factor is
(53) |
As shown in Fig. 3, the normalized weight factors are actually dominated by the one corresponding to the smallest eigenvalue of . The resonant line shape is very well approximated by a single Lorentizan function, whose linewidth is the real part of the smallest eigenvalue of . In the following, we will focus on the smallest eigenvalue of , 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 linewidth as a function of pumping rate and spin exchange rate. At low spin exchange rates (e.g., ), the linewidth monotonically increases with pumping rate. In contrast, at high spin exchange rates (), 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 ()
To analyze the linewidth in the low polarization limit, the dimensionless matrix Eq. (44) is expanded in the power series of up to the second order as
(54) |
where
(55) | ||||
(56) | ||||
(57) |
In Eqs. (55)-(57), we have expanded the matrix to the second order of , i.e.,
(58) |
The leading order matrix is Hermitian, whose eigenvalue and corresponding eigenvector are obtained by solving the eigenvalue problem
(59) |
The smallest eigenvalue is denoted as . The correction to due to the high-order matrices and is calculated using perturbation theory. Up to the second order of , the minimum eigenvalue is
(60) |
where
(61) | ||||
(62) | ||||
(63) |
For a given nuclear spin quantum number , in the low polarization limit, the leading order contribution of linewidth is
(64) |
where two slowing-down factors and are defined as
(65) | ||||
(66) |
The spin exchange slowing-down factor in Eq. (66) is identical to that obtained in the spin-exchange relaxation free (SERF) regime, while the spin destruction slowing-down factor 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 is
(67) |
Note that the first order correction is proportional to the optical pumping rate , and it is independent on the spin-exchange rate , 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 .
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 is lengthy. Here, we show the second order correction for (e.g., ) as an example
(68) |
The second order correction brings about a local maximum of the linewidth, as long as . 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 , the local maximum exists when the ratio is above a threshold value . This agrees with the analytic result in Eq. (68), where becomes positive if the spin exchange rate is too small. For weak spin exchange rate (e.g., in low density vapor) with , 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 ()
Similar analysis is performed in the high polarization limit, where the matrix is expanded as
(69) |
where
(70) | ||||
(71) | ||||
(72) |
In Eqs. (71) & (72), we have used the expansion of the transverse spin exchange matrix
(73) |
Note that the matrix is not Hermitian. With the eigenvalues and the corresponding right- and left-eigenvectors and of the leading order matrix , which satisfy the following equations
(74) | ||||
(75) |
the minimum eigenvalues expanded in a power series of is
(76) |
where
(77) | ||||
(78) | ||||
(79) |
With Eqs. (77)-(79), the linewidth in the high polarization limit is
(80) |
where is factor depending on the nuclear spin quantum number . The complete expressions of linewidth for different isotopes are listed in Table 1. Although the coefficient in front of can be small in the high polarization limit with , the third term in Eq. (80) is not negligible if the spin exchange rate is large.
isotope | ||||
---|---|---|---|---|
13.82 | ||||
, , , , | 6.22 | |||
2.77 | ||||
1.66 |
Equation (80) implies that the linewidth has a local minimum when changing the optical pumping rate . Figure 5 shows the position of the local minimum for different values of the spin exchange rate . At high spin exchange rate (i.e., ), 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 . For weak spin exchange rate with , the line width monotonically increases with increasing optical pumping rate . To observe the light-narrowing effect, the spin exchange rate must be larger than the threshold value, that is, . The threshold values 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 , while, according to Eq. (80), the linewidth at the local minimum scales as , i.e.,
(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]
(82) |
where or . Indeed, for a given Zeeman resonance labeled with , the linewidth in the RZR regime is the diagonal element of the matrix in Eq. (II.4.1). A direct comparison between the linewidth in the UZR and RZR regimes for the case is shown in Fig. 7. Here, we focus on the narrowest Zeeman resonance, i.e., .
In the low polarization limit , the linewidth in the RZR regime is
(83) |
where the slowing-down factors are defined as
(84) | ||||
(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 depends on the ratio . As increasing , the coefficient become negative and the light-narrowing effect occurs when
(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
(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 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 , as functions of the spin polarization parameter and the relative spin-exchange strength . 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 and the spin destruction rate ) 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
- Budker and Kimball [2013] D. Budker and D. F. J. Kimball, Optical Magnetometry (Cambridge University Press, 2013).
- Wright et al. [2022] M. J. Wright, L. Anastassiou, C. Mishra, J. M. Davies, A. M. Phillips, S. Maskell, and J. F. Ralph, Frontiers in Physics 10 (2022).
- Walker and Larsen [2016] T. G. Walker and M. S. Larsen, in Advances in Atomic, Molecular and Optical Physics, Vol. 65 (2016) pp. 373–401 .
- Safronova et al. [2018] M. S. Safronova, D. Budker, D. DeMille, D. F. J. Kimball, A. Derevianko, and C. W. Clark, Reviews of Modern Physics 90, 025008 (2018).
- Happer and Tam [1977] W. Happer and A. C. Tam, Physical Review A 16, 1877 (1977).
- Bhaskar et al. [1981] N. D. Bhaskar, J. Camparo, W. Happer, and A. Sharma, Physical Review A 23, 3048 (1981).
- Appelt et al. [1999] S. Appelt, A. B.-A. Baranga, A. R. Young, W. Happer, A. Ben-Amar Baranga, A. R. Young, W. Happer, A. B.-A. Baranga, A. R. Young, and W. Happer, Physical Review A 59, 2078 (1999).
- Appelt et al. [1998] S. Appelt, A. Baranga, C. Erickson, M. Romalis, A. Young, and W. Happer, Physical Review A 58, 1412 (1998).
- Schiff and Snyder [1939] L. I. Schiff and H. Snyder, Physical Review 55, 59 (1939).
- Julsgaard et al. [2004] B. Julsgaard, J. Sherson, J. L. Sørensen, and E. S. Polzik, Journal of Optics B: Quantum and Semiclassical Optics 6, 5 (2004).
- Jau et al. [2004] Y. Y. Jau, A. B. Post, N. N. Kuzma, A. M. Braun, M. V. Romalis, and W. Happer, Physical Review Letters 92, 110801 (2004).
- Happer [1972] W. Happer, Reviews of Modern Physics 44, 169 (1972).
- Happer et al. [2010] W. Happer, Y.-Y. Jau, and T. G. Walker, Optically Pumped Atoms (Wiley-VCH, Weinheim, 2010).
- Seltzer [2008] S. J. Seltzer, Developments in Alkali-Metal Atomic Magnetometry, Ph.D. thesis (2008).