We investigate the relationship between the Rabi frequency (ΩRF, related to the applied electric field) and Autler-Townes (AT) splitting, when performing atom-based radio-frequency (RF) electric (E) field strength measurements using Rydberg states and electromagnetically induced transparency (EIT) in an atomic vapor. The AT splitting satisfies, under certain conditions, a well-defined linear relationship with the applied RF field amplitude. The EIT/AT-based E-field measurement approach derived from these principles is currently being investigated by several groups around the world as a means to develop a new SI-traceable RF E-field measurement technique. We establish conditions under which the measured AT-splitting is an approximately linear function of the RF electric field. A quantitative description of systematic deviations from the linear relationship is key to exploiting EIT/AT-based atomic-vapor spectroscopy for SI-traceable field measurement. We show that the linear relationship is valid and can be used to determine the E-field strength, with minimal error, as long as the EIT linewidth is small compared to the AT-splitting. We also discuss interesting aspects of the thermal dependence (i.e., hot- versus cold-atom) of this EIT-AT technique. An analysis of the transition from cold- to hot-atom EIT in a Doppler-mismatched cascade system reveals a significant change of the dependence of the EIT linewidth on the optical Rabi frequencies and of the AT-splitting on ΩRF.

Recently, there has been an interest in the development of a novel spectroscopic approach for electric (E) field measurements1–9 that can lead to a direct International System of Units (SI) traceable, self-calibrated measurement, that has the capability to perform measurements on a fine spatial resolution. The approach utilizes the phenomena of electromagnetically induced transparency (EIT) and Autler-Townes (AT) splitting.1–3,10 This technique has the capability of becoming a new international standard for E-field measurements and calibrations. For a new measurement method to be accepted by National Metrology Institutes, including the National Institute of Standards and Technology (NIST), the accuracy and uncertainties of the approach must be assessed. Knowing the uncertainties of a technique is an important step when establishing a new international measurement standard for an E-field strength and is a necessary step for this method to be accepted as a standard calibration technique. Here, we look at one aspect of the measurement uncertainties/errors of this E-field measurement approach. We further investigate the effects of the atom temperature on the radio frequency (RF) E-field measurement by investigating the transitional behavior of EIT spectra from the cold-atom ( T = 1 μ K) to the hot-atom (T = 300 K) case. Our analysis sheds light on several important aspects, including the emergence of the Doppler correction factor in hot-atom, vapor-cell EIT, and a change in the dependence of EIT linewidth on the optical Rabi frequencies.

This measurement approach can be represented by the four-level atomic system shown in Fig. 1. In effect, the “probe” laser is used to probe the response of the ground-state transition of the atoms, and a second laser (“coupling” laser) is used to excite the atoms to a high energy state. In the presence of the coupling laser, a destructive quantum interference occurs and the atoms become transparent to the resonant probe laser (this is the concept of EIT). A transparency window is opened for the probe laser light: probe light transmission is increased. The coupling laser wavelength is chosen such that the atom is in a sufficiently high state (a Rydberg state) such that a radio frequency (RF) field couples two Rydberg states (levels 3 and 4 in Fig. 1). The RF field in the four-level atomic system causes constructive interference of excitation pathways within the EIT transmission window, resulting in a decreased transmission of the probe laser and splitting (AT splitting) of the EIT peak. A typical measured spectrum with and without an RF source is shown in Fig. 2. In this figure, Δp is the detuning of the probe laser from the ground state transition of the atom (where Δ p = ω p ω o ; ωo is the on-resonance angular frequency of the ground state transition and ωp is the angular frequency of the probe laser).

FIG. 1.

Illustration of a four-level system and the vapor cell setup for measuring EIT, with counter-propagating probe and coupling beams. The RF field propagates transversely to the optical beams in the vapor cell.

FIG. 1.

Illustration of a four-level system and the vapor cell setup for measuring EIT, with counter-propagating probe and coupling beams. The RF field propagates transversely to the optical beams in the vapor cell.

Close modal
FIG. 2.

Experimental data for the EIT signal (i.e., probe laser transmission through the cell) as a function of probe laser detuning Δ p . This dataset is for a RF of 15.095 GHz and corresponds to the following 85Rb 4-level atomic system: 5 S 1 / 2 5 P 3 / 2 52 D 5 / 2 53 P 3 / 2 .

FIG. 2.

Experimental data for the EIT signal (i.e., probe laser transmission through the cell) as a function of probe laser detuning Δ p . This dataset is for a RF of 15.095 GHz and corresponds to the following 85Rb 4-level atomic system: 5 S 1 / 2 5 P 3 / 2 52 D 5 / 2 53 P 3 / 2 .

Close modal

The AT splitting (defined as 2 π Δ f o ) of the probe laser spectrum is easily measured and under certain conditions is equal to the Rabi frequency of the RF transition, or

AT plitting = 2 π Δ f o = Ω R F ,
(1)

where Ω R F = | E | / is the Rabi frequency of the RF transition, is Planck's constant, and is the dipole moment of the RF atomic transition. This relationship between the AT splitting and the Rabi frequency is obtained under the weak probe limit and for no Doppler averaging. By measuring this splitting ( Δ f m ), we get a direct measurement of the RF E-field strength. In this approach, either the probe or the coupling laser can be scanned or detuned. For either case, the E-field strength is given by Refs. 2 and 3 

| E | = 2 π D Δ f m = 2 π Δ f o ,
(2)

where Δ f m is the measured splitting, Δ f o = D Δ f m , and D is a parameter whose value depends on which of the two lasers is scanned during the measurement. If the probe laser is scanned, D = λ p λ c ; λp and λc are the wavelengths of the probe and coupling laser, respectively. The λ p / λ c ratio is needed to account for the Doppler mismatch of the probe and coupling lasers10 when the probe laser is scanned. If the coupling laser is scanned, it is not required to correct for the Doppler mismatch, and as a result D = 1. We consider this type of measurement of the E-field strength a direct SI-traceable, self-calibrated measurement in that it is related to Planck's constant (which will become an SI-defined quantity by standard bodies in the near future) and only requires a frequency measurement ( Δ f m , which can be measured very accurately). The atomic dipole moment is a parameter which can be calculated accurately.2,11

Various aspects of the uncertainties of this type of measurement are currently being investigated.1–3,7,14,15 In this paper, we investigate the validity of Eq. (2) and determine conditions when it is valid. Under certain conditions, the linear relationship between Δ f o and | E | as predicted from Eq. (2) can break down and the relationship becomes non-linear and, hence, calls into question the validity of Eq. (2). One such case is shown in Fig. 3, where we plot Δ f o versus | E | . (The quantity “ W ” on the x-axis is proportional to the applied RF E-field, where W represents the RF generator power.) The results in Fig. 3 are for measurements performed in 133Cs at 10.7 GHz. We see that when | E | is small (or for small Δ f o ), the data exhibit a nonlinear behavior. Linearity only applies for large values of | E | ; i.e., the linear relationship fails for small | E | . This same nonlinear behavior has been observed by others.16–18 The expected linear behavior predicted in Eq. (2) is based on the assumption that the AT splitting of the EIT signal is proportional to the RF Rabi frequency (ΩRF). The validity of Eq. (2) is a key aspect of the uncertainty of this measurement approach. When the data in the non-linear region are used with Eq. (2), then the E-field would be either underestimated or overestimated. The reason this non-linear behavior is observed is due to the linewidth of the EIT signal. If the EIT linewidth is comparable to or larger than Δ f o , then the linear relationship fails. In most of these cases, including the non-linear regime in Fig. 3, the AT splitting is artificially pulled closer together for simultaneously large EIT linewidths and small E-field values (or small Δ f o ). We investigate these effects by performing numerical studies of the four-level EIT signal for a range of conditions. We also explore how the temperature of the atomic vapor affects the linewidth of the EIT signal and thus the range of the non-linear behavior. We note that the splitting observed in the case of off-resonant AT effect is intrinsically non-linear;19 in the present work, we focus on non-linear AT splittings observed in the case of on-resonant AT effect.

FIG. 3.

Example of non-linear behavior. This dataset is for a RF of 10.7 GHz and corresponds to the following 133Cs 4-level atomic system: 6 S 1 / 2 6 P 3 / 2 41 D 5 / 2 42 P 3 / 2 . The three different symbols correspond to three different datasets. The numbers (square root of power, i.e., W ) stated on the x-axis are the values from the signal generator (SG) that feeds the waveguide which, in turn, feeds the horn antenna. As such, the quantity “ W ” on the x-axis is proportional to the applied RF E-field.

FIG. 3.

Example of non-linear behavior. This dataset is for a RF of 10.7 GHz and corresponds to the following 133Cs 4-level atomic system: 6 S 1 / 2 6 P 3 / 2 41 D 5 / 2 42 P 3 / 2 . The three different symbols correspond to three different datasets. The numbers (square root of power, i.e., W ) stated on the x-axis are the values from the signal generator (SG) that feeds the waveguide which, in turn, feeds the horn antenna. As such, the quantity “ W ” on the x-axis is proportional to the applied RF E-field.

Close modal

In the model and the results presented in the following, we use parameters for Rb 85 as a typical, specific example. Similar results apply to other typical cases, such as Rb 87 (which merely differs in partial vapor pressure) and Cs 133 . We start by noting that the power of the probe beam measured on the detector (the EIT signal, i.e., the probe transmission through the vapor cell) is given by13 

P = P 0 exp ( 2 π L Im [ χ ] λ p ) = P 0 exp ( α L ) ,
(3)

where P0 is the power of the probe beam at the input of the cell, L is the length of the cell, λp is the wavelength of the probe laser, χ is the susceptibility of the medium seen by the probe laser, and α = 2 π Im [ χ ] / λ p is the Beer's absorption coefficient for the probe laser. The susceptibility for the probe laser is related to the density matrix component (ρ21) that is associated with the | 1 | 2 transition by the following:13 

χ = 2 N 0 12 E p ϵ 0 ρ 21 D ,
(4)

where the subscript D on ρ12 presents a Doppler averaged value, N 0 is the total density of atoms in the cell and is given by

N 0 = 0.7217 p k B T ,
(5)

where kB is the Boltzmann constant, T is the temperature in Kelvin, and the pressure p (in units of Pa) is given by12 

p = 10 5.006 + 4.857 4215 T .
(6)

The factor 0.7217 in Eq. (5) reflects the natural abundance of 85Rb. Note that this density is for vapor above solid Rb. A different expression applies for vapor above liquid Rb.12 In Eq. (4), 12 is the dipole moment for the | 1 | 2 transition, ϵ0 is the vacuum permittivity, and Ep is the amplitude of the probe laser E-field.

The density matrix component (ρ21) is obtained from the master equation13 

ρ ̇ = ρ t = i [ H , ρ ] + L ,
(7)

where H is the Hamiltonian of the atomic system under consideration and L is the Lindblad operator that accounts for the decay processes in the atom.

For the four level system, the Hamiltonian can be expressed as

H = 2 [ 0 Ω p 0 0 Ω p 2 Δ p Ω c 0 0 Ω c 2 ( Δ p + Δ c ) Ω R F 0 0 Ω R F 2 ( Δ p + Δ c + Δ R F ) ] ,
(8)

where Δp, Δc, and ΔRF are the detunings of the probe laser, couple laser, and the RF source, respectively; and Ωp, Ωc, and ΩRF are the Rabi frequencies associated with the probe laser, coupling laser, and the RF field, respectively. The detuning for each field is defined as

Δ p , c , R F = ω p , c , R F ω o p , c , R F ,
(9)

where ω o p , c , R F are the on-resonance angular frequencies of transitions | 1 | 2 , | 2 | 3 , and | 3 | 4 , respectively; and ω p , c , R F are the angular frequencies of the probe laser, coupling laser, and the RF source, respectively. The Rabi frequencies are defined as Ω p , c , R F = | E p , c , R F | p , c , R F , where | E p , c , R F | are the magnitudes of the E-field of the probe laser, the coupling laser, and the RF source, respectively. Finally, p , c , and R F are the atomic dipole moments corresponding to the probe, coupling, and RF transitions, respectively.

For the four-level system, the L matrix is given by

L = [ Γ 2 ρ 22 γ 12 ρ 12 γ 13 ρ 13 γ 14 ρ 14 γ 21 ρ 21 Γ 3 ρ 33 Γ 2 ρ 22 γ 23 ρ 23 γ 24 ρ 24 γ 31 ρ 31 γ 32 ρ 32 Γ 4 ρ 44 Γ 3 ρ 33 γ 34 ρ 34 γ 41 ρ 41 γ 42 ρ 42 γ 43 ρ 43 Γ 4 ρ 44 ] ,
(10)

where γ i j = ( Γ i + Γ j ) / 2 and Γ i , j are the transition decay rates. Since the purpose of the present study is to explore the intrinsic limitations of Rydberg-EIT field sensing in vapor cells, no collision terms or dephasing terms are added. While Rydberg-atom collisions, Penning ionization, and ion electric fields can, in principle, cause dephasing, such effects can, for instance, be alleviated by reducing the beam intensities, lowering the vapor pressure, or limiting the atom-field interaction time. In this analysis, we set Γ 1 = 0 , Γ 2 = 2 π × (6 MHz), Γ 3 = 2 π × (3 kHz), and Γ 4 = 2 π × (2 kHz). Γ 2 is for the D2 line in 85Rb,12 and Γ 3 and Γ 4 are typical Rydberg decay rates.

Substituting H and L into Eq. (7), the diagonal density matrix components are given by

ρ ̇ 11 = i Ω p 2 ( ρ 12 ρ 21 ) + Γ 2 ρ 22 ρ ̇ 22 = i Ω p 2 ( ρ 12 ρ 21 ) + i Ω c 2 ( ρ 23 ρ 32 ) Γ 2 ρ 22 + Γ 3 ρ 33 ρ ̇ 33 = i Ω c 2 ( ρ 23 ρ 32 ) + i Ω R F 2 ( ρ 34 ρ 43 ) Γ 3 ρ 33 + Γ 4 ρ 44 ρ ̇ 44 = i Ω R F 2 ( ρ 34 ρ 43 ) Γ 4 ρ 44
(11)

and the off-diagonal density matrix components are given by

ρ ̇ 21 = [ i Δ p γ 21 ] ρ 21 + i Ω p 2 ( ρ 22 ρ 11 ) i Ω c 2 ρ 31 ρ ̇ 31 = [ i ( Δ p + Δ c ) γ 31 ] ρ 31 + i Ω p 2 ρ 32 i Ω c 2 ρ 21 i Ω R F 2 ρ 41 ρ ̇ 41 = [ i ( Δ p + Δ c + Δ R F ) γ 41 ] ρ 41 + i Ω p 2 ρ 42 i Ω R F 2 ρ 31 ρ ̇ 32 = [ i Δ c γ 32 ] ρ 32 + i Ω c 2 ( ρ 33 ρ 22 ) + i Ω p 2 ρ 31 i Ω R F 2 ρ 42 ρ ̇ 42 = [ i ( Δ c + Δ R F ) γ 42 ] ρ 42 + i Ω p 2 ρ 41 + i Ω c 2 ρ 43 i Ω R F 2 ρ 32 ρ ̇ 43 = [ i Δ R F γ 43 ] ρ 43 + i Ω R F 2 ( ρ 44 ρ 33 ) + i Ω c 2 ρ 42 ,
(12)

where ρ ̇ i j * = ρ ̇ j i and * represents the complex conjugate.

We numerically solve these equations to find the steady-state solution for ρ21 for various values of Ωc, Ωp, and ΩRF. This is done by forming a matrix with the system of equations for ρ ̇ i j = 0 . The null-space of this resulting system matrix is the steady-state solution. The steady-state solution for ρ21 is Doppler averaged in the usual way

ρ 21 D = 1 π u 3 u 3 u ρ 21 ( Δ p , Δ c ) e v 2 u 2 d v ,
(13)

where u = 2 k B T / m and m is the mass of the atom. We use the case where the probe and coupling laser are counter-propagating. Thus, the frequency seen by the atom moving toward the probe beam is upshifted by 2 π v / λ p (where v is the velocity of the atoms), while the coupling beam is downshifted by 2 π v / λ c . Thus, the probe and coupling beam detuning are modified by the following:

Δ p = Δ p 2 π λ p v and Δ c = Δ c + 2 π λ c v .
(14)

In this section, we investigate the AT-splitting (i.e., Δ f o ) dependence on Ωc and Ωp. Figure 4 shows the EIT signal obtained from Eq. (3) for Ω R F = 2 π × (9.2 MHz) and different combinations of Ωp and Ωc. In our discussion, we use the combined quantity Ω = Ω c 2 + Ω p 2 . The results are for T = 300 K and for scanning the probe laser (indicated by Δp). Once the Doppler mismatch is taken into account, the separation between the two peaks for Ω R F = 2 π × (9.2 MHz) and probe detuning should be 5.66 MHz (i.e., Δ f m = Δ f o λ c / λ p , where Δ f o = Ω R F , assuming the validity of the linear relationship). This figure shows that the AT-splitting (or Δ f o ) is not constant for a given value of ΩRF for the different values for Ω. We see that the AT-splitting is either smaller or larger than ΩRF depending on the values of Ωp and Ωc. Recall that Eqs. (1) and (2) state that the AT-splitting should be constant for a given value of ΩRF. In the figure, we also indicate the percent difference between the measured splitting and ΩRF defined by

Δ % = 100 × ( 2 π × Δ f o ) Ω R F Ω R F and Δ f o = D Δ f m .
(15)
FIG. 4.

Numerical data for the EIT signal (i.e., probe laser transmission through the cell) as a function of probe laser detuning Δp. This dataset is for Ω R F = 2 π × (9.2 MHz) and for different values Ω = Ω c 2 + Ω p 2 and with on resonant RF and coupling fields ( Δ c = Δ R F = 0 ).

FIG. 4.

Numerical data for the EIT signal (i.e., probe laser transmission through the cell) as a function of probe laser detuning Δp. This dataset is for Ω R F = 2 π × (9.2 MHz) and for different values Ω = Ω c 2 + Ω p 2 and with on resonant RF and coupling fields ( Δ c = Δ R F = 0 ).

Close modal

Notice that Δ % ranges from 30 % to 19 % for the different values of Ω, where a negative value indicates an underestimate and a positive value indicates an overestimate.

This non-linear behavior is further illustrated with the family of curves shown in Figs. 5(a) for Ω p = Ω c . The dashed reference line in this figure corresponds to Δ f o = Ω R F . For a given Ωc and Ωp, the linear behavior occurs only as ΩRF exceeds some minimum value. Below this value the curves become non-linear. The percentage difference ( Δ % ) for these data is shown in Fig. 6(a). We see Δ % can be large in the non-linear region. Using Δ f o in these non-linear regions could result in large errors in the calculated E-field strength.

FIG. 5.

Numerical calculations for Δ f o versus Ω R F / 2 π for T = 300 K: (a) Δ f o versus Ω R F / 2 π for Ω p = Ω c . (b) Δ f o versus Ω R F / 2 π for Ω p = 2 π × (1.0 MHz). (c) Δ f o versus Ω R F / 2 π for Ω c = 2 π × (1.0 MHz).

FIG. 5.

Numerical calculations for Δ f o versus Ω R F / 2 π for T = 300 K: (a) Δ f o versus Ω R F / 2 π for Ω p = Ω c . (b) Δ f o versus Ω R F / 2 π for Ω p = 2 π × (1.0 MHz). (c) Δ f o versus Ω R F / 2 π for Ω c = 2 π × (1.0 MHz).

Close modal
FIG. 6.

Numerical calculations for Δ % versus Ω R F / 2 π for T = 300 K: The gray boxes in these curves indicate the region where | Δ % | < 1 % . (a) Δ % versus Ω R F / 2 π for Ω p = Ω c . (b) Δ % versus Ω R F / 2 π for Ω p = 2 π × (1.0 MHz). (c) Δ % versus Ω R F / 2 π for Ω c = 2 π × (1.0 MHz).

FIG. 6.

Numerical calculations for Δ % versus Ω R F / 2 π for T = 300 K: The gray boxes in these curves indicate the region where | Δ % | < 1 % . (a) Δ % versus Ω R F / 2 π for Ω p = Ω c . (b) Δ % versus Ω R F / 2 π for Ω p = 2 π × (1.0 MHz). (c) Δ % versus Ω R F / 2 π for Ω c = 2 π × (1.0 MHz).

Close modal

Figure 5(b) shows a set of curves for Ωp held fixed at 2 π × (1.0 MHz) and various Ωc, and Fig. 5(c) shows a set of curves for Ωc held fixed at 2 π × (1.0 MHz) and various Ωp. The percentage difference ( Δ % ) for these data are shown in Figs. 6(b) and 6(c). The results in Fig. 5(b) [i.e., for Ω p = 2 π × (1.0 MHz)] have an additional interesting feature in that the curves first overshoot the Δ f o = Ω R F reference line before dropping to Δ f o 0 . The overshoot increases with Ωc. The overshooting is caused by the enhanced absorption dips in the EIT signal. These enhanced absorption dips were first observed in Ref. 10 and are seen by the dashed curve in Fig. 4 at the Δ p = 2 π × (7.5 MHz) location. An enhanced absorption dip from one EIT line interferes with the peak of another EIT line, causing the peaks to be pushed further apart and resulting in a larger AT-splitting. As indicated in Fig. 6(b), the enhanced absorption dips cause much larger overestimates (i.e., positive Δ % ) when compared to the other cases [i.e., Figs. 6(a) and 6(c)].

The full-width half-maxima (FWHM) of the EIT linewidths (ΓEIT) for each curve are indicated in Figs. 5 and 6. (These values for ΓEIT were obtained for Ω R F = 0 .) We see that values for ΩRF where the non-linear behavior begins (i.e., the curves drop below the Δ f o = Ω R F line) decrease as ΓEIT decreases (or as Ω decreases). The highly non-linear behavior occurs when Ω R F Γ E I T . The AT peaks are artificially pulled closer together for simultaneously large ΓEIT and small E-field values (or small Δ f o ). We see that for the cases in Fig. 5(c), as long as Ω R F 2 Γ E I T , the curves fall within the linear regions to within 1 % . For the cases with enhanced absorption where there are large overshoots [see Fig. 5(b)], the curves are more complicated and ΩRF needs to be larger than 3 Γ E I T to 4 Γ E I T in order for the curves to start approaching the linear region to within 1 % . These deviations from the linear behavior for all the cases in Fig. 5 result in a measurement uncertainty when determining the E-field strength with this EIT method. In Fig. 6, we have highlighted where these deviations are below 1 % . This is indicated by the gray box region, which corresponds to | Δ % | < 1 % . For cases in Fig. 6(c), we see as long as Ω R F 2 Γ E I T , the deviation or the error in using Eq. (2) is less than 1% (i.e., Δ % < 1 % ), and if Ω R F 2 Γ E I T the deviation is decreased further to well below 0.5 % . However, the cases in Figs. 6(b) and 6(a) are more problematic and can result in larger deviations. For the cases in Fig. 6(a), Δ % > 1 % unless Ω R F > 3 Γ E I T ; for the cases in Fig. 6(b), Δ % > 1 % unless ΩRF is larger than 3 Γ E I T or 4 Γ E I T . These results indicate that the most useful region to work within is shown in the cases in Fig. 6(c). An additional benefit of these cases is that they have the smallest value of Ωc. In fact, powers for the coupling laser to reach Rydberg states can be a critical experimental limitation, so working in regions where Ωc is small is convenient.

From the results in Fig. 5, it is seen that this non-linear behavior also leads to a critical value for ΩRF (defined as Ωcrit) before AT-splitting occurs. These critical values correspond to the minimum splitting that can be detected and is given by the locations where the curves touch the x-axis. The E-field cannot be measured with the AT-splitting method when Ω R F < Ω crit . No matter how small the value of Ω, the value of ΩRF, where non-linear behavior begins, never approaches zero. This is due to the Doppler effect as a result of the temperature of the atomic vapor. The reason Ωcrit never approaches zero can be seen in Fig. 7, where we show the calculated EIT linewidth versus Ω for both our room temperature atomic vapor and a cold-atom scenario. We find that, as Ω is reduced, ΓEIT approaches a minimum value for the hot-atom case, whereas it reaches near-zero for cold atoms. The Ωcrit has a lower limit of about 1.8 MHz for the hot-atom 85Rb case studies in this work. In Sec. IV, we explain how the temperature of the atomic vapor changes the measurement of an RF E-field by affecting the characteristics of the AT-splitting.

FIG. 7.

Numerical data for ΓEIT versus Ω. We show results for both hot-atoms (T =300 K) and cold-atoms (T =1 μK). These values were obtained for Ω R F = 0 .

FIG. 7.

Numerical data for ΓEIT versus Ω. We show results for both hot-atoms (T =300 K) and cold-atoms (T =1 μK). These values were obtained for Ω R F = 0 .

Close modal

The region of the non-linear behavior is related to the temperature of the atomic vapor. From the results shown in Fig. 5, we see that as Ω gets smaller, the location where the non-linear region starts is pulled to smaller values of ΩRF. However, at some point, no matter how small Ω becomes, we cannot force the start of the non-linear region to zero. This is related to the Doppler effect of the room temperature atomic vapor and the natural linewidth of the probe transition. A question of interest is how vapor-cell (hot-atom) EIT emerges from the behavior of cold-atom EIT. In cold-atom EIT, there is no Doppler factor [i.e., D = 1 in Eq. (2)], and for small Ω, the EIT linewidth follows Γ E I T = ( Ω c 2 + Ω p 2 ) / Γ 2 (see Fig. 7). This allows, in principle, the measurement of very small RF-induced AT splittings, limited only by laser and Rydberg-level linewidths. The system gradually transitions from the cold- into the hot-atom case when the Doppler shifts in the sample become on the order of the linewidth of the probe transition. In the case of Rb Rydberg-EIT on the | 5 S 1 / 2 | 5 P 3 / 2 | R y d cascade, the transition occurs between 1 mK and 1 K. In the hot-atom limit, the velocity averaging in Eq. (13) produces the Doppler factor in Eq. (2) and limits the EIT linewidth to values above several MHz as shown in Fig. 7. This, in turn, limits the ability to measure RF fields via the AT-splitting method to fields with ΩRF larger than about 2 MHz, as evident in Fig. 5. Deviations from a linear relationship between the observed AT-splitting and ΩRF, as well as the loss in splitting contrast at low ΩRF and Ωc also follow from the velocity averaging.

Figure 8 explains the transition from cold- to hot-atom EIT by considering the absorption coefficient per velocity class [which is proportional to the integrand in Eq. (13)] as a function of atom velocity and the probe detuning Δp for representative Rabi-frequency choices. In the cold-atom case (a), the velocity is practically zero, and the overall scattering profile is a Lorentzian with an FWHM of Γ 2 along the Δp-axis. With a small laser linewidth, the Rydberg-EIT linewidth, Γ E I T = ( Ω c 2 + Ω p 2 ) / Γ 2 , can be made very narrow; in Fig. 8(a) Γ E I T = 2 π × 300 kHz. According to Section II, at a fixed velocity, the AT-splitting is given by ΩRF. The observable EIT spectrum follows from the integral along the velocity axis in Fig. 8(a), leading to two narrow EIT dips within a 2 π × 6 MHz wide Lorenztian. It is thus seen that at low temperature, atom-based RF electric-field measurement via determination of AT splittings in EIT spectra, and using Eq. (2) with a Doppler factor D = 1, can be performed down to very small values of ΩRF. While there is a lower limit set by laser linewidth, interaction time, level dephasing, etc., lower limits for ΩRF in the range 100  kHz seem quite feasible. For a (typical) transition dipole moment of 1000 ea0 (where e is the elementary charge and a0 is the Bohr radius), this translates into an RF-field sensitivity in the range of 10  mV/m. When scanning the probe laser, there is an upper limit to the strongest RF fields that can be measured; here, the cold-atom Rydberg-EIT-AT method requires Ω R F 2 π × 10  MHz, due to the finite probe-transition linewidth Γ 2 . This results in an upper limit for the RF field that can be measured. However, that limitation can be easily eliminated by scanning the coupler-laser detuning, Δc, and using a fixed Δ p = 0 .

FIG. 8.

Calculated absorption per velocity, in arb. un., vs atom velocity and probe detuning Δp: (a) in a cold-atom case with T = 10 μ K, Ω p / ( 2 π ) = 1 MHz, Ω c / ( 2 π ) = 1 MHz and Ω R F / ( 2 π ) = 5 MHz, (b) in a cold-atom case with T = 10 mK and with the same parameters as in case (a), and (c) in a hot-atom (vapor-cell) case with T = 300 K, Ω p / ( 2 π ) = 5 MHz, Ω c / ( 2 π ) = 20 MHz, and Ω R F / ( 2 π ) = 50 MHz. The dashed lines in (b) and (c) show the functions in Eq. (16).

FIG. 8.

Calculated absorption per velocity, in arb. un., vs atom velocity and probe detuning Δp: (a) in a cold-atom case with T = 10 μ K, Ω p / ( 2 π ) = 1 MHz, Ω c / ( 2 π ) = 1 MHz and Ω R F / ( 2 π ) = 5 MHz, (b) in a cold-atom case with T = 10 mK and with the same parameters as in case (a), and (c) in a hot-atom (vapor-cell) case with T = 300 K, Ω p / ( 2 π ) = 5 MHz, Ω c / ( 2 π ) = 20 MHz, and Ω R F / ( 2 π ) = 50 MHz. The dashed lines in (b) and (c) show the functions in Eq. (16).

Close modal

Atom-based RF electric-field measurement via determination of AT splittings in cold-atom EIT spectra benefits from the absence of significant Doppler effects. A typical case of hot-atom EIT is shown in Fig. 8(c). The dressed-atom energy levels that arise from Eq. (14) correspond to the following lines in the velocity vs Δp map:

Ground state : v 1 ( Δ p ) = 0 Intermediate state : v 2 ( Δ p ) = λ p Δ p 2 π AT - split Rydberg pairs : v R y d ( Δ p ) = 1 2 π ( Δ p ± Ω R F 2 ) λ p λ c λ p λ c .
(16)

The lines for ground and Rydberg states are largely hidden in the plots in Figs. 8(b) and 8(c) because those levels do not scatter photons, except near the dressed-state avoided crossings. The coupler field leads to avoided crossings between intermediate and Rydberg states; in the vicinity of the crossings, the intermediate-state character is partially shared with the Rydberg levels, making them “visible” in the absorption map. At fixed velocity, the AT-split Rydberg levels are still separated by ΩRF (same as in the cold-atom case). Also, away from the crossings, the width of the absorption profile associated with the intermediate level is Γ 2 .

The measurable EIT spectrum arises from the velocity integral of the results given in Fig. 8. The observed splitting between the EIT-AT dips is given by the separation of the projections of the red dots in Fig. 8(c) onto the x-axis. Straightforward analysis leads to Eq. (2) with the Doppler factor D. It also becomes obvious that the velocity integration leads to a minimum observable ΩRF that is given by Γ 2 , Ωc, and Ωp. Basically, Ωc cannot be reduced to arbitrarily small values without losing the EIT absorption dips after performing the velocity integral. This is visualized in Fig. 8(b), where ΓEIT is so low that the EIT signatures get lost in the velocity averaging process. The lowest Ωc at which the EIT lines can still be observed scales with Γ 2 . This, in turn, sets a lower limit of ΩRF below which the RF field will not cause an AT splitting of the EIT line. From Eq. (16), it is further noted that in a Doppler matched case (i.e., λ p = λ c ), the velocities of the Rydberg dressed states versus Δp in Fig. 8 would be vertical lines. In this case, the EIT features would not get washed out in the velocity averaging process. Hence, the lower limit of ΩRF can be avoided by selecting a Doppler-matched case.

The EIT spectra in Fig. 9 show the Beer's absorption coefficient [as defined in Eq. (3)] in a Rb sample of fixed density for a range of temperatures. The top four traces all use the same combination of Ω p / 2 π = Ω c / 2 π = 1 MHz and Ω R F / 2 π = 5 MHz. For traces 1, 2, 3, 4, the temperature is varied from T = 10 μ K, 1 mK, 0.1 K, to 10 K, and the measured AT-splitting for each is 5.0 MHz, 4.8 MHz, 3.8 MHz, and 3.1 MHz. At 10 μ K, the measured AT-splitting is equal to the incident ΩRF, for a Doppler factor of D = 1, while beyond 10 K the measured AT-splitting is a fraction of ΩRF, leading to the Doppler correction of D = λ p / λ c as in Eq. (2). Trace 5 shows that this factor remains as Ωp, Ωc, and ΩRF are varied, and trace 6 shows that D does not change as the temperature is increased from 10 K to 300 K.

FIG. 9.

Calculated EIT spectra vs. probe detuning for a range of temperatures. To show the effect of temperature on the probe absorption coefficient, α, the atom density is held fixed at 7.3 × 10 9 cm– 3. The parameters T, Ω p / ( 2 π ) , Ω c / ( 2 π ) , and Ω R F / ( 2 π ) are, for the plots labeled 1–6, (1) 10 μ K, 1, 1, and 5 MHz; (2) 1 mK, 1, 1, and 5 MHz; (3) 0.1 K, 1, 1, and 5 MHz; (4) 10 K, 1, 1, and 5 MHz; (5) 10 K, 5, 5, and 10 MHz; and (6) 300 K, 5, 20, and 50 MHz. The frequencies of the AT-splitting for each case are indicated on the figure.

FIG. 9.

Calculated EIT spectra vs. probe detuning for a range of temperatures. To show the effect of temperature on the probe absorption coefficient, α, the atom density is held fixed at 7.3 × 10 9 cm– 3. The parameters T, Ω p / ( 2 π ) , Ω c / ( 2 π ) , and Ω R F / ( 2 π ) are, for the plots labeled 1–6, (1) 10 μ K, 1, 1, and 5 MHz; (2) 1 mK, 1, 1, and 5 MHz; (3) 0.1 K, 1, 1, and 5 MHz; (4) 10 K, 1, 1, and 5 MHz; (5) 10 K, 5, 5, and 10 MHz; and (6) 300 K, 5, 20, and 50 MHz. The frequencies of the AT-splitting for each case are indicated on the figure.

Close modal

Another feature that can be seen in Fig. 9 is that the visibility of the AT-splitting varies with Ω. As Ω is reduced, both the minimum observable RF Rabi frequency Ωcrit and the visibility of the AT-splitting are reduced. This introduces another trade-off for measuring a weak RF E-field. As Ω is lowered to detect a smaller minimum ΩRF, the ability to resolve the AT-splitting is also reduced. The visibility V for a range of values of Ω is plotted in Fig. 10. Here we define visibility as V = ( α m a x α m i n ) / α m a x , where αmin and αmax are the minimum absorption coefficients at the bottom of the EIT dips and the maximum coefficient in between the two dips, respectively. For Rb, requiring a splitting visibility V 0.1 leads to a lowest value of Ω R F that can be measured on the order of 3 MHz. This limit applies to RF E-field measurements via the Rydberg-EIT-AT splitting method in hot vapor cells, which is the most common and robust method for the atom-based RF field measurement. Lower limits can be achieved in less robust methods, such as cold-atom EIT, a detailed signal analysis of Rydberg-EIT lines in RF fields that are too weak to split the EIT line into a resolved AT pair,3 or a 3-photon EIT method.20,21

FIG. 10.

Visibility V of the Rydberg-EIT-AT splitting, as defined in the text, vs RF Rabi frequency ΩRF for the indicated values of Ωp and Ωc in a 300-K vapor cell. It is seen that high coupler Rabi frequencies lead to better visibility of the AT-split EIT lines at high ΩRF, but greatly reduce the applicability of the Rydberg-EIT-AT splitting method at low ΩRF.

FIG. 10.

Visibility V of the Rydberg-EIT-AT splitting, as defined in the text, vs RF Rabi frequency ΩRF for the indicated values of Ωp and Ωc in a 300-K vapor cell. It is seen that high coupler Rabi frequencies lead to better visibility of the AT-split EIT lines at high ΩRF, but greatly reduce the applicability of the Rydberg-EIT-AT splitting method at low ΩRF.

Close modal

As an atomic vapor transitions from the cold- to the hot-atom case, the Doppler factor D in Eq. (2) changes from D = 1 to D = λ p / λ c ; for Rb, this transition occurs between 100 μ K and 1 K. For small values of Ω, the EIT linewidth scaling with the pump and probe Rabi frequencies changes from Γ E I T = ( Ω p 2 + Ω c 2 ) / Γ 2 to a more complicated relationship that has a lower limit for ΓEIT of 1.5 MHz (as shown in Fig. 7). Consequently, there is a minimum RF Rabi frequency (Ωcrit) above which an AT-splitting in the hot-atom case can be observed. The visibility of the AT-splitting is likewise affected, setting an effective minimum resolvable ΩRF for this particular method.

We explored the conditions under which the linear relationship between the measured Autler-Townes splitting ( Δ f o ) and the RF Rabi frequency (ΩRF) is valid. This is key to assessing the uncertainty of this measurement approach. We show that the linear relationship of these two quantities is valid as long as ΩRF is greater than twice the EIT linewidth ( 2 Γ E I T ). In this region, we determined that the error in the assumed linear relationship can be limited to < 1 % , and if Ω R F 2 Γ E I T , the deviation is decreased further to well below 1 % . While the AT split is resolvable when the RF Rabi frequency is below 2 Γ E I T , the pump and probe Rabi frequencies can no longer be neglected, causing the relationship between Δ f o and ΩRF to become non-linear. The non-linear behavior also leads to a critical value for ΩRF, which corresponds to the minimum splitting that can be detected, and as a result, an E-field cannot be measured with the AT-splitting method when Ω R F < Ω crit . If we used the linear relationship [i.e., Eq. (2)] to estimate the E-field strength when ΩRF approaches Ωcrit, the uncertainties in the field strength measurements can be very large (over 30 % to 50 % ). With the knowledge of the probe and coupling Rabi frequencies, the relationship between Δ f o and ΩRF can, in principle, be calculated with the model presented here, reducing the measurement uncertainty when ΩRF is below the 2 Γ E I T limit, but a closed-form expression cannot be used (because it is not available for these cases).

Another strategy to reduce this limit is to reduce the EIT linewidth. We find that reducing the probe and coupling Rabi frequencies reduces the EIT linewidth only to a minimum that is limited by the Doppler broadening. To reduce this further requires lowering the temperature of the atomic vapor. In our discussion of the effect of temperature on the relationship between Δ f o and ΩRF, we show that the dependence of ΓEIT on Ωp and Ωc is different above and below a temperature threshold. Below this threshold, ΓEIT can be made arbitrarily small, allowing measurements of very weak ΩRF. Above this threshold, there is a lower limit on ΓEIT, which limits the range in which the relationship between Δ f o and ΩRF is linear.

Finally, the 3-photon scheme discussed in Ref. 20 allows, in principle, for narrower EIT linewidths (ΓEIT) in vapor-cell Rydberg EIT and for measurements of weaker E-field strengths (i.e., smaller ΩRF). The technique in Ref. 20 is based on a coincidental wavelength combination and has shown theoretically to produce linewidths of about 10 kHz. To observe a kHz linewidth, one must decrease the power to avoid power broadening. This will increase the shot noise, reducing the sensitivity.20,21 The technique has not been realized experimentally as of yet, but does show promise.

This work was partially supported by the Defense Advanced Research Projects Agency (DARPA) under the QuASAR Program and by NIST through the Embedded Standards program, as well as by the NSF (PHY-1506093) and Rydberg Technologies, LLC. We thank Dr. J. P. Shaffer from The University of Oklahoma for his useful technical discussions on 3-photon Rydberg excitation.

1.
C. L.
Holloway
,
M. T.
Simons
,
J. A.
Gordon
,
P. F.
Wilson
,
C. M.
Cooke
,
D. A.
Anderson
, and
G.
Raithel
,
IEEE Trans. Electromagn. Compat.
59
(
2
),
717
728
(
2017
).
2.
C. L.
Holloway
,
J. A.
Gordon
,
A.
Schwarzkopf
,
D. A.
Anderson
,
S. A.
Miller
,
N.
Thaicharoen
, and
G.
Raithel
,
IEEE Trans. Antenna Propag.
62
(
12
),
6169
6182
(
2014
).
3.
J. A.
Sedlacek
,
A.
Schwettmann
,
H.
Kübler
,
R.
Löw
,
T.
Pfau
, and
J. P.
Shaffer
,
Nat. Phys.
8
,
819
(
2012
).
4.
C. L.
Holloway
,
J. A.
Gordon
,
A.
Schwarzkopf
,
D. A.
Anderson
,
S. A.
Miller
,
N.
Thaicharoen
, and
G.
Raithel
,
Appl. Phys. Lett.
104
,
244102
(
2014
).
5.
J. A.
Gordon
,
C. L.
Holloway
,
A.
Schwarzkopf
,
D. A.
Anderson
,
S. A.
Miller
,
N.
Thaicharoen
, and
G.
Raithel
,
Appl. Phys. Lett.
105
,
024104
(
2014
).
6.
J. A.
Sedlacek
,
A.
Schwettmann
,
H.
Kübler
, and
J. P.
Shaffer
,
Phys. Rev. Lett.
111
,
063001
(
2013
).
7.
H.
Fan
,
S.
Kumar
,
J.
Sedlacek
,
H.
Kübler
,
S.
Karimkashi
, and
J. P.
Shaffer
,
J. Phys. B: At. Mol. Opt. Phys.
48
,
202001
(
2015
).
8.
M.
Tanasittikosol
,
J. D.
Pritchard
,
D.
Maxwell
,
A.
Gauguet
,
K. J.
Weatherill
,
R. M.
Potvliege
, and
C. S.
Adams
,
J. Phys. B
44
,
184020
(
2011
).
9.
C. G.
Wade
,
N.
Sibalic
,
N. R.
de Melo
,
J. M.
Kondo
,
C. S.
Adams
, and
K. J.
Weatherill
,
Nat. Photonics
11
,
40
43
(
2017
).
10.
A. K.
Mohapatra
,
T. R.
Jackson
, and
C. S.
Adams
,
Phys. Rev. Lett.
98
,
113003
(
2007
).
11.
M. T.
Simons
,
J. A.
Gordon
, and
C. L.
Holloway
,
J. Appl. Phys.
120
,
123103
(
2016
).
12.
D. A.
Steck
, see http://steck.us/alkalidata for Rubidium 85 D line data (revision 2.1.6, 20 September
2013
).
13.
P. R.
Berman
and
V. S.
Malinovsky
,
Priciples of Laser Spectroscopy and Quantum Optics
(
Princeton Univerity Press
,
2011
).
14.
C. L.
Holloway
,
J. A.
Gordon
,
M. T.
Simons
,
H.
Fan
,
S.
Kumar
,
J. P.
Shaffer
,
D. A.
Anderson
,
A.
Schwarzkopf
,
S. A.
Miller
,
N.
Thaicharoen
, and
G.
Raithel
, in
EMC 2015: Joint IEEE International Symposium on Electromagnetic Compatibility and EMC Europe, Dresden, Germany, August 16–22, 2015
, pp.
467
472
.
15.
H.
Fan
,
S.
Kumar
,
J.
Sheng
,
J. P.
Shaffer
,
C. L.
Holloway
, and
J. A.
Gordon
,
Phys. Rev. Appl.
4
,
044015
(
2015
).
16.
H.
Zhang
,
L.
Zhang
,
L.
Wang
,
S.
Bao
,
J.
Zhao
, and
G.
Raithel
,
Phy. Rev. A
90
,
043849
(
2014
).
17.
P. M.
Anisimov
,
J. P.
Dowling
, and
B. C.
Sanders
,
Phy Rev. Lett.
107
,
163604
(
2011
).
18.
T. Y.
Abi-Salloum
,
J. Mod. Opt.
57
(
14
),
1366
1376
(
2010
).
19.
M. T.
Simons
,
J. A.
Gordon
,
C. L.
Holloway
,
D. A.
Anderson
,
S. A.
Miller
, and
G.
Raithel
,
Appl. Phys. Lett.
108
,
174101
(
2016
).
20.
S.
Kumar
,
H.
Fan
,
J.
Sheng
, and
J. P.
Shaffer
,
Sci. Rep.
7
,
42981
(
2017
).
21.
S.
Kumar
,
H.
Fan
,
H.
Kübler
,
A.
Jozani
, and
J. P.
Shaffer
,
Opt. Express
25
,
8625
(
2017
).