License: CC BY 4.0
arXiv:2402.15114v1 [physics.atom-ph] 23 Feb 2024

Cavity based non-destructive detection of photoassociation in a dark MOT

V. I. Gokul    Arun Bahuleyan    S. P. Dinesh    V. R. Thakar    S. A. Rangwala Raman Research Institute, C.V. Raman Avenue, Sadashivanagar, Bangalore 560080, India
(February 23, 2024)
Abstract

The photoassociation (PA) of rubidium dimer (Rb22{}_{2}start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT) in a dark magneto-optic trap (MOT) is studied using atom-cavity collective strong coupling. This allows non-destructive detection of the molecule formation process as well as rapid and repeated interrogation of the atom-molecule system. The vacuum Rabi splitting (VRS) measurements from the bright MOT are carefully calibrated against equivalent measurements with fluorescence. Further loading rates in dark MOT are determined using VRS. This method provides a reliable, fast, and non-destructive detection scheme for ultracold molecules when the atoms are non-fluorescing using the free atoms coupled to a cavity.

preprint: APS/123-QED

I Introduction

The detection of multiparticle interaction is a challenging experimental problem. In an ensemble, this can be done by detecting the post-collision product or by interrogating individual participants. This methods only allow partial information about the process of interest. For atomic systems, cavity-based measurements are very precise for studying atom-field interactions. In this work, we adapt cavity QED technique to detect the photoassociation (PA) of atoms. We aim to demonstrate the power of cavity based technique with PA as a representative example, though this technique has wider application, for instance, collision rate measurement between the two overlapping ensembles in the cavity mode volume. PA results when three particles, two free atoms in a scattering state absorb one photon and form a bound, excited molecule [1, 2](Fig. 1). The absorbed photon is resonant with the free-to-bound transition, and the excited state molecule can subsequently decay into free atoms by emitting a photon (Γf)subscriptΓ𝑓(\Gamma_{f})( roman_Γ start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ) or into a bound, ground-state molecule in one of many possible ro-vibrational states (Γb)subscriptΓ𝑏(\Gamma_{b})( roman_Γ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ). PA is efficient in an ensemble of ultracold atoms, where the energy of the scattering state is well-defined, favoring resonant excitation by a single, narrow-band laser. Suitable conditions for PA, therefore, require trapped atoms in magneto-optical traps (MOT’s) [3, 4, 5, 6, 7, 8], magnetic traps [9, 10], or dipole traps [11, 12, 13].

Detection of PA in a bright MOT (where the MOT is observable in fluorescence) is done by measuring the resonant loss of the fluorescing atoms. Fluorescence detection is not an option in a dark MOT [14] and for other dark traps. Here, resonant multi-photon ionization (ReMPI) of the resulting ground state molecule and the subsequent detection of the molecular ion is the main detection tool available [15, 16, 9, 17, 18]. ReMPI has limited efficiency principally because the ground-state molecules produced by PA are created in various states due to spontaneous emission from the excited state.

In this paper, we present cavity-based, non-destructive detection of the atoms, which are collectively strongly coupled to a low finesse cavity, during PA. In the experiment, we probe the vacuum Rabi splitting (VRS) of the atoms coupled to the cavity [19, 20]. As will be shown below, the technique works very well for a dark MOT, where atomic fluorescence is not detectable. In addition, this measurement is non-destructive, as it probes the atoms and not the molecules, with very low light intensity detuned with respect to the atomic resonance, resulting in no additional loss of atom from the trap. This combination of properties allows for continuous monitoring of PA, which we demonstrate so that its temporal evolution can also be measured and studied. Below, we first explain our experimental setup, calibrate our technique with a bright MOT, demonstrate results with the dark MOT, discuss the different ratios of bright and dark MOT resonance intensities, and conclude with some applications of this method. Our technique complements the recent direct strong coupling signature of PA [21] observed in a high finesse (F5×104)similar-to𝐹5superscript104(F\sim 5\times 10^{4})( italic_F ∼ 5 × 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ) cavity. Even though the finesse in our experiment doesn’t allow for a direct strong coupling signature of PA, detecting atoms instead of molecules allows us to detect PA and calculate rates. A relatively low finesse (330)similar-toabsent330(\sim 330)( ∼ 330 ) cavity used in the experiment implies that one can perform these measurements even by having a cavity outside the vacuum.

Refer to caption
Figure 1: Potential energy curves for Rb22{}_{2}start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT molecule relevant in the experiment and various transitions. b) Trap-loss spectrum for different PA transitions in bright MOT. Here, PA frequency is referenced to F=3𝐹3F=3italic_F = 3 to F=4superscript𝐹4F^{\prime}=4italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 4 D2222 transition of Rb.

II Experimental setup

The experimental setup consists of a Fabry-Perot (FP) cavity co-centered with a dilute gas of ultracold 8585{}^{85}start_FLOATSUPERSCRIPT 85 end_FLOATSUPERSCRIPTRb atoms that are trapped in an MOT. A detailed overview of the experimental setup is discussed in previous work [22, 23]. The FP cavity has length L45.7approximately-equals-or-equals𝐿45.7L\approxeq 45.7italic_L ≊ 45.7 mm, finesse F330similar-to𝐹330F\sim 330italic_F ∼ 330, and linewidth κ4.4similar-to𝜅4.4\kappa\sim 4.4italic_κ ∼ 4.4 MHz. A cooling beam detuned 1212-12- 12 MHz from F=3𝐹3F=3italic_F = 3 to F=4superscript𝐹4F^{\prime}=4italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 4 and repumper on F=2𝐹2F=2italic_F = 2 to F=3superscript𝐹3F^{\prime}=3italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 3 forms the bright MOT. The dark MOT is implemented by obstructing the center of an independent repumper beam with 2absent2\approx 2≈ 2 mm disk shelving the atoms in the F=2𝐹2F=2italic_F = 2 non-fluorescing state. A double pass acousto-optic modulator (AOM) is used to switch both repumping lasers, allowing us to change between bright and dark MOT in our system. Additionally, a weak depumping beam (1μ(\sim 1\mu( ∼ 1 italic_μW)))) from F=3𝐹3F=3italic_F = 3 to F=3superscript𝐹3F^{\prime}=3italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 3 is used to maximize dark MOT density. For a cooling power of 25252525 mW and repumping power of 3.53.53.53.5 mW measured peak density of dark MOT, ρ01.7×1010subscript𝜌01.7superscript1010\rho_{0}\approx 1.7\times 10^{10}italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ≈ 1.7 × 10 start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPT cm33{}^{-3}start_FLOATSUPERSCRIPT - 3 end_FLOATSUPERSCRIPT. For all dark MOT measurements, the cavity is referenced to F=2𝐹2F=2italic_F = 2 to F=3superscript𝐹3F^{\prime}=3italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 3 transition [24], and a weak probe beam is scanned across the same transition. The cavity output is monitored using a photomultiplier tube (PMT) and a charge-coupled device (CCD) camera. A separate free-space setup is used to collect MOT fluorescence [23]. The PA beam has 2similar-toabsent2\sim 2∼ 2 mm diameter and is incident perpendicular to the cavity axis. The wavelength of this PA laser is continuously monitored using a wavemeter.

In the experiment, κ𝜅\kappaitalic_κ and excited state decay (Γ)Γ(\Gamma)( roman_Γ ) are much greater than single atom-cavity coupling constant g0=μ232ω23/2ϵ0Vc0.141subscript𝑔0superscriptsubscript𝜇232subscript𝜔232Planck-constant-over-2-pisubscriptitalic-ϵ0subscript𝑉𝑐0.141g_{0}=\sqrt{\mu_{23}^{2}\omega_{23}/{2\hbar\epsilon_{0}V_{c}}}\approx 0.141italic_g start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = square-root start_ARG italic_μ start_POSTSUBSCRIPT 23 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ω start_POSTSUBSCRIPT 23 end_POSTSUBSCRIPT / 2 roman_ℏ italic_ϵ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG ≈ 0.141 MHz for F=2F=3𝐹2superscript𝐹3F=2\rightarrow F^{\prime}=3italic_F = 2 → italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 3 transition [25], where μ23subscript𝜇23\mu_{23}italic_μ start_POSTSUBSCRIPT 23 end_POSTSUBSCRIPT is the transition dipole matrix element for probe transition and Vcsubscript𝑉𝑐V_{c}italic_V start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT is the cavity mode volume. So, while a single atom cannot couple strongly to the cavity [26, 19] when Ncsubscript𝑁𝑐N_{c}italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT atoms are present in Vcsubscript𝑉𝑐V_{c}italic_V start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, each of these atoms can couple to a single cavity photon, resulting in an effective coupling strength of gN=g0Ncsubscript𝑔𝑁subscript𝑔0subscript𝑁𝑐g_{N}=g_{0}\sqrt{N_{c}}italic_g start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT = italic_g start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT square-root start_ARG italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG [27, 28, 19, 29]. In this collective strong coupling regime, the condition for observing strong coupling effects becomes gN=g0NcΓ,κformulae-sequencesubscript𝑔𝑁subscript𝑔0subscript𝑁𝑐much-greater-thanΓ𝜅g_{N}=g_{0}\sqrt{N_{c}}\gg\Gamma,\kappaitalic_g start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT = italic_g start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT square-root start_ARG italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG ≫ roman_Γ , italic_κ. For a MOT with density distribution ρ(x,y,z)𝜌𝑥𝑦𝑧\rho(x,y,z)italic_ρ ( italic_x , italic_y , italic_z ) co-centered with a cavity mode, Ncsubscript𝑁𝑐N_{c}italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT can be calculated using Nc=ρ(x,y,z)|ψ(x,y,z)|2𝑑Vsubscript𝑁𝑐𝜌𝑥𝑦𝑧superscript𝜓𝑥𝑦𝑧2differential-d𝑉N_{c}=\int\rho(x,y,z)|\psi(x,y,z)|^{2}dVitalic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = ∫ italic_ρ ( italic_x , italic_y , italic_z ) | italic_ψ ( italic_x , italic_y , italic_z ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_V, where ψ𝜓\psiitalic_ψ is the cavity mode function [23].

For our experimental parameters, atom numbers as low as Nc1000greater-than-or-equivalent-tosubscript𝑁𝑐1000N_{c}\gtrsim 1000italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ≳ 1000 can be measured via VRS in a collective strong coupling regime. Since the separation between VRS peaks is given by 2g0Nc2subscript𝑔0subscript𝑁𝑐2g_{0}\sqrt{N_{c}}2 italic_g start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT square-root start_ARG italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG, any change in VRS corresponds to a change in Ncsubscript𝑁𝑐N_{c}italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, which is a direct measure of change in MOT atom number [29]. This allows detection of PA transitions in dark MOT using VRS.

III Results

III.1 PA resonances and loading rates in bright MOT

The potential energy curves (PEC) for long-range photoassociated Rb22{}_{2}start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT  [30, 31] molecule in D2 transition are shown in fig. 1. From the selection rule only transitions relevant and observed in our experiments are 0gsuperscriptsubscript0𝑔0_{g}^{-}0 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT, 0g+superscriptsubscript0𝑔0_{g}^{+}0 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT and 1gsubscript1𝑔1_{g}1 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT  [32]. In the presence of a PA laser, two ground state Rb atoms can form a weakly bound excited state Rb22{}_{2}start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT molecule (0gsuperscriptsubscript0𝑔0_{g}^{-}0 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT, 0g+superscriptsubscript0𝑔0_{g}^{+}0 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT or 1gsubscript1𝑔1_{g}1 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT). This excited state molecule decays back into a high vibrational ground state molecule (Σg+1superscriptsuperscriptsubscriptΣ𝑔1{}^{1}\Sigma_{g}^{+}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT roman_Σ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT or Σu+3superscriptsuperscriptsubscriptΣ𝑢3{}^{3}\Sigma_{u}^{+}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT roman_Σ start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT) or two free atoms by photoemission.

In the experiment, bright MOT has a peak density of 8×1010similar-toabsent8superscript1010\sim 8\times 10^{10}∼ 8 × 10 start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPT cm33{}^{-3}start_FLOATSUPERSCRIPT - 3 end_FLOATSUPERSCRIPT and FWHM of 180similar-toabsent180\sim 180∼ 180 μ𝜇\muitalic_μm. The PA laser is tuned from 12814128141281412814 cm11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT to 12816128161281612816 cm11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT and six prominent PA resonances are addressed. These transitions are labeled as pa1-pa6. The peak intensity of the PA laser is 15similar-toabsent15\sim 15∼ 15 W///cm22{}^{2}start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT, which is sufficient to saturate the PA transition. Fig. 1 (b) shows the trap-loss measurement for these PA transitions. The dashed lines in the figure represent the closest calculated vibrational level for the excited state ab initio molecular potentials [31]. As this method is non-destructive, one can perform continuous measurements to study the evolution and dynamics of the system. As proof of principle, we measure the loading rates in MOT using VRS and calculate loss rates in the system during PA.

Refer to caption
Figure 2: Loading curve for bright MOT a) time evolution of VRS as the MOT loads b) loading curve from MOT fluorescence. Here, the green curve is when PA is off-resonant, and the blue curve is when PA is on-resonance. The loading rate ΓΓ\Gammaroman_Γ and γ𝛾\gammaitalic_γ are calculated from the fit. c) loading curve measured from VRS. Each point is an average of 5 repeated measurements and the error bar is the standard deviation. All the measurements are done on pa5 transition.
Refer to caption
Figure 3: Direct detection of PA resonances in dark MOT using VRS. Here,(a)-(f) shows the change in VRS for various PA resonances. The blue line is when the PA laser is kept on-resonance, and the green plot is when the PA laser is kept off-resonant. The figure shows the VRS for a single sweep and the theoretical fit of the data to determine the peak separation. (g) Loading curve using atoms coupled to cavity mode(Nc𝑐{}_{c}start_FLOATSUBSCRIPT italic_c end_FLOATSUBSCRIPT). Here, the green circle is when the PA laser is kept off-resonant to the transition, and the blue circle is when the PA laser is kept on-resonant. (h) The ratio of Ncsubscript𝑁𝑐N_{c}italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT when PA laser is on-transition(Nc;1subscript𝑁𝑐1N_{c;1}italic_N start_POSTSUBSCRIPT italic_c ; 1 end_POSTSUBSCRIPT) to off-transition(Nc;0subscript𝑁𝑐0N_{c;0}italic_N start_POSTSUBSCRIPT italic_c ; 0 end_POSTSUBSCRIPT) as a function of PA intensity. The arrow represents the peak intensity used for dark MOT loading measurements. Each point is an average of 5 measurements. All the measurements are done on the pa5 transition.

The MOT fluorescence with a bright MOT provides a direct calibration of the VRS detection. As bright MOT atoms are predominantly in F=3𝐹3F=3italic_F = 3 state, the cavity is locked to F=3𝐹3F=3italic_F = 3 to F=3superscript𝐹3F^{\prime}=3italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 3 transition, and the probe is scanned across the same transition. For this measurement, value of g0~=μ332ω33/2ϵ0Vc0.13~subscript𝑔0superscriptsubscript𝜇332subscript𝜔332Planck-constant-over-2-pisubscriptitalic-ϵ0subscript𝑉𝑐0.13\tilde{g_{0}}=\sqrt{\mu_{33}^{2}\omega_{33}/{2\hbar\epsilon_{0}V_{c}}}\approx 0% .13over~ start_ARG italic_g start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG = square-root start_ARG italic_μ start_POSTSUBSCRIPT 33 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ω start_POSTSUBSCRIPT 33 end_POSTSUBSCRIPT / 2 roman_ℏ italic_ϵ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG ≈ 0.13 MHz. The frequency of PA transition is monitored by a wavemeter. Change in the VRS signal is observed in addition to the trap-loss fluorescence signal on PA resonances. The ground state molecules and majority of the free atoms formed from excited state molecules during a PA are not trapped in a MOT [33], resulting in a reduction of atom number in MOT. This results in a smaller VRS when the PA laser is resonant to any of the transitions. This reduction in VRS measures the formation of excited state Rb22{}_{2}start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT molecules by PA.

As shown in fig 2, a direct measure of the loading curve using VRS is also performed for bright MOT with PA laser kept off and on transition. Here, the probe laser is scanned at a rate of 8 Hz for 10 seconds and frequency is referenced with respect to saturated absorption spectroscopy (SAS). The MOT coil is turned off for a short time, and the probe output is continuously monitored. As a result, once the MOT starts to load rapid VRS measurements are possible. The loading rate is calculated from the VRS signal. Fig 2 a) shows the evolution of VRS as the MOT loads. Initially, when there are no atoms in the MOT, the cavity transmission corresponds to empty cavity transmission. As MOT builds up, the number of atoms in the cavity mode increases, resulting in an increase in VRS.

Fig. 2 (b) shows the MOT loading from fluorescence when the PA laser is kept on-resonance (blue) and off-resonance (green). This is fitted with N(t)=N0(1exp(Γt))𝑁𝑡subscript𝑁01Γ𝑡N(t)=N_{0}\left(1-\exp\left({-\Gamma t}\right)\right)italic_N ( italic_t ) = italic_N start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 1 - roman_exp ( - roman_Γ italic_t ) ) (see appendix) to extract loading rate (ΓΓ\Gammaroman_Γ). When the PA laser is kept off-resonance, the loading rate is γ=0.41±0.03𝛾plus-or-minus0.410.03\gamma=0.41\pm 0.03italic_γ = 0.41 ± 0.03 s11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT. Similarly, loading rate with on-resonant PA is Γ=γ+γPA=0.83±0.04Γ𝛾subscript𝛾𝑃𝐴plus-or-minus0.830.04\Gamma=\gamma+\gamma_{PA}=0.83\pm 0.04roman_Γ = italic_γ + italic_γ start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT = 0.83 ± 0.04 s11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT. The loss rate in MOT due to PA is therefore as γPA=0.42±0.07subscript𝛾𝑃𝐴plus-or-minus0.420.07\gamma_{PA}=0.42\pm 0.07italic_γ start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT = 0.42 ± 0.07 s11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT. Fig. 2 (c) shows the same loading measured using VRS with on-resonant (blue circle) and off-resonant (green circle) PA. From the fit, the loading rate with off-resonant PA is γ=0.43±0.02𝛾plus-or-minus0.430.02\gamma=0.43\pm 0.02italic_γ = 0.43 ± 0.02 s11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT and with on-resonant PA is Γ=0.80±0.11Γplus-or-minus0.800.11\Gamma=0.80\pm 0.11roman_Γ = 0.80 ± 0.11 s11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT. Hence, the loss rate in MOT calculated from VRS is γPA=0.37±0.13subscript𝛾𝑃𝐴plus-or-minus0.370.13\gamma_{PA}=0.37\pm 0.13italic_γ start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT = 0.37 ± 0.13 s11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT. These values agree very well with MOT fluorescence measurement, validating VRS-based measurements.

III.2 PA resonances and loading rates in dark MOT

In systems like dark MOT, most of the atoms are in non-fluorescing state (F=2𝐹2F=2italic_F = 2). As a result, a direct fluorescence trap-loss signal for PA resonances is not possible. The VRS measurement makes such detection possible even when the system is in a dark state. For this detection, the cavity is locked to F=2𝐹2F=2italic_F = 2 to F=3superscript𝐹3F^{\prime}=3italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 3 transition. A weak probe laser is scanned across the same transition and referenced to SAS. Since the atoms are in F=2𝐹2F=2italic_F = 2 state, all the PA resonances (pa1-pa6) are shifted by 3.03similar-toabsent3.03\sim 3.03∼ 3.03 GHz, which is the difference in frequency between F=2𝐹2F=2italic_F = 2 and F=3𝐹3F=3italic_F = 3 hyperfine levels. The PA laser with peak intensity 15similar-toabsent15\sim 15∼ 15 W///cm22{}^{2}start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT is incident on the dark MOT. The cavity transmission is monitored with PA kept on and off-resonant for all six transitions.

Figure 3 (a)-(f) shows the change in VRS signal when the PA laser is on-resonance (blue) and off-resonance (green) for all the transitions. In all cases, the blue curve is shifted for better visibility. The values of all VRS and corresponding molecular transitions are given in the appendix (Table 1). As seen from fig 3 the value of VRS decreases when a PA transition is addressed, which is a direct measurement of the formation of Rb*2superscriptsubscriptabsent2{}_{2}^{*}start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT * end_POSTSUPERSCRIPT molecule due to PA in a dark MOT.

To demonstrate the time evolution capabilities, the loading rates for dark MOT are measured in the presence of the PA. For this measurement, the probe laser is locked to F=2𝐹2F=2italic_F = 2 to F=3superscript𝐹3F^{\prime}=3italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 3 transition using SAS. This locked laser is passed through two double-pass AOMs. One of the AOMs is scanned continuously at a rate of 4444 Hz for 20202020 seconds. For this loading measurements peak intensity of the PA laser was kept 10similar-toabsent10\sim 10∼ 10 W///cm22{}^{2}start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT. Figure 3 (g) shows the Ncsubscript𝑁𝑐N_{c}italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT as a function of time when the PA laser is kept off-resonance (green circle) and on-resonance (blue circle). Loading rate of dark MOT on resonance, ΓD=γD+γPA;D=0.33±0.06subscriptΓ𝐷subscript𝛾𝐷subscript𝛾𝑃𝐴𝐷plus-or-minus0.330.06\Gamma_{D}=\gamma_{D}+\gamma_{PA;D}=0.33\pm 0.06roman_Γ start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT = italic_γ start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT + italic_γ start_POSTSUBSCRIPT italic_P italic_A ; italic_D end_POSTSUBSCRIPT = 0.33 ± 0.06 s11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT, where γD=0.13±0.03subscript𝛾𝐷plus-or-minus0.130.03\gamma_{D}=0.13\pm 0.03italic_γ start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT = 0.13 ± 0.03 s11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT, for pa5, as shown in fig 3(g). From this, the loss rate in dark MOT due to PA is calculated as γPA;D=0.20±0.09subscript𝛾𝑃𝐴𝐷plus-or-minus0.200.09\gamma_{PA;D}=0.20\pm 0.09italic_γ start_POSTSUBSCRIPT italic_P italic_A ; italic_D end_POSTSUBSCRIPT = 0.20 ± 0.09 s11{}^{-1}start_FLOATSUPERSCRIPT - 1 end_FLOATSUPERSCRIPT. Figure 3 (h) shows the ratio of Ncsubscript𝑁𝑐N_{c}italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT when PA is on and off-resonant to the transition as a function of PA laser intensity. As the PA intensity increases the fraction of atoms getting converted to PA molecule increases, showing a reduction in the ratio of Ncsubscript𝑁𝑐N_{c}italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT value. The dashed line, which is a linear fit shows good agreement for lower PA intensity to theoretical expectation. For higher intensity saturation effect results in deviation from the linear fit. This saturation in PA rate for higher PA intensity is well known both theoretically [34] and experimentally [35, 10, 36].

III.3 Temperature measurement during PA

Refer to caption
Figure 4: Temperature measurement during PA. Here, the green circle represents measurement with PA, and the blue circle without PA. Each data point is an average of 5 measurements.

A small fraction of atoms that are dissociated from excited molecules can have kinetic energy less than the trap depth. As a result, these atoms don’t contribute to the trap-loss signal. However, if they have sufficient energy to increase the average temperature of the MOT [37], an extension of this detection technique is possible. To demonstrate this, temperature measurement using cavity VRS [22] is done for bright MOT with and without PA. The cavity is locked to F=3F=4𝐹3superscript𝐹4F=3\rightarrow F^{\prime}=4italic_F = 3 → italic_F start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 4 transition. The probe laser is locked to the same transition, and AOM is scanned at a rate of 2222 kHz. MOT beams are switched off using an AOM-based switch for 5555 ms. During this period, MOT expands ballistically, and from the measured value of Ncsubscript𝑁𝑐N_{c}italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT the standard deviation (σ)𝜎(\sigma)( italic_σ ) of MOT can be calculated. This gives a direct measure of temperature using the expression [38], σ(τ)2=σ02+kBTmτ2𝜎superscript𝜏2superscriptsubscript𝜎02subscript𝑘𝐵𝑇𝑚superscript𝜏2\sigma(\tau)^{2}=\sigma_{0}^{2}+\frac{k_{B}T}{m}\tau^{2}italic_σ ( italic_τ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG italic_k start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT italic_T end_ARG start_ARG italic_m end_ARG italic_τ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. The same procedure is repeated with an on-resonant PA, which is turned off just before the MOT expansion. This ensures that the atoms that are recaptured in the trap contribute to the temperature of the MOT. As shown in fig 4, the measured value of temperature of the MOT is TMOT=20.10±0.66subscript𝑇𝑀𝑂𝑇plus-or-minus20.100.66T_{MOT}=20.10\pm 0.66italic_T start_POSTSUBSCRIPT italic_M italic_O italic_T end_POSTSUBSCRIPT = 20.10 ± 0.66 μ𝜇\muitalic_μK and with on-resonant PA is TPA=21.86±1.01subscript𝑇𝑃𝐴plus-or-minus21.861.01T_{PA}=21.86\pm 1.01italic_T start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT = 21.86 ± 1.01 μ𝜇\muitalic_μK, which shows that fraction of atoms recaptured by the decay of the excited molecules is negligible.

IV Discussion & Conclusion

The advantages of PA from an optically dark ensemble are numerous, the most prominent is the larger fraction of atoms converted to molecules (fig.3), since the dark MOT atoms are all in the same ground state. Further in systems like dark MOT, the PA rate (γPA;D)subscript𝛾𝑃𝐴𝐷\left(\gamma_{PA;D}\right)( italic_γ start_POSTSUBSCRIPT italic_P italic_A ; italic_D end_POSTSUBSCRIPT ) is larger compared to the loading rate (γD)subscript𝛾𝐷\left(\gamma_{D}\right)( italic_γ start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ), giving a higher fraction of PA (see appendix). The cavity-based technique of detection which probes the free atoms instead of the molecules directly [39, 40], can be used very effectively for detection, in scenarios where detection is difficult and in cases like 3-particle PA [41] or PA of Rydberg molecules [42] with accuracy.

This technique is expected to be particularly useful where direct fluorescence detection is not possible. As the method is non-destructive by nature, atom-molecule population dynamics in the system can be studied continuously. Since ultracold molecules are not trapped in our experiment, a direct detection of atom-molecule collision was not possible. However, the extension of this method for such experiments is straightforward, making this way of studying interactions more universal. In summary, the present work demonstrates how cavity techniques can be applied to detect complex processes and yield accurate results in very challenging systems and spatially compact geometries.

Acknowledgements.
The authors would like to thank the Department of Science and Technology and Ministry of Electronics and Information Technology (MeitY), Government of India, under a Centre for Excellence in Quantum Technologies grant with Ref. No. 4(7)/2020-ITEA.

Appendix A Loading rates in MOT

The rate equation for the number of atoms (N𝑁Nitalic_N) in the trap with PA laser is given by [43, 35]

dNdt=LγN(β+βPA)n2(r)d3r𝑑𝑁𝑑𝑡𝐿𝛾𝑁𝛽subscript𝛽𝑃𝐴superscript𝑛2𝑟superscript𝑑3𝑟\frac{dN}{dt}=L-\gamma N-\left(\beta+\beta_{PA}\right)\int n^{2}(r)d^{3}rdivide start_ARG italic_d italic_N end_ARG start_ARG italic_d italic_t end_ARG = italic_L - italic_γ italic_N - ( italic_β + italic_β start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT ) ∫ italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r ) italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_r (1)

where, L is the rate at which atoms are loaded into the trap, γN𝛾𝑁\gamma Nitalic_γ italic_N is the rate at which atoms are lost from the trap due to background collisions, β𝛽\betaitalic_β is the loss rate due to collision with trapped atoms, and βPAsubscript𝛽𝑃𝐴\beta_{PA}italic_β start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT is the loss rate due to photoassociation. Dark MOT has a uniform density [23] (nmax)subscript𝑛𝑚𝑎𝑥\left(n_{max}\right)( italic_n start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT ) inside the trap. As a result, eq 1 can be rewritten as

dNdt=L(γ+(β+βPA)nmax)NN(t)=N0(1exp(Γt))𝑑𝑁𝑑𝑡𝐿𝛾𝛽subscript𝛽𝑃𝐴subscript𝑛𝑚𝑎𝑥𝑁𝑁𝑡subscript𝑁01Γ𝑡\begin{split}\frac{dN}{dt}&=L-\left(\gamma+\left(\beta+\beta_{PA}\right)n_{max% }\right)N\\ N(t)&=N_{0}\left(1-\exp\left({-\Gamma t}\right)\right)\end{split}start_ROW start_CELL divide start_ARG italic_d italic_N end_ARG start_ARG italic_d italic_t end_ARG end_CELL start_CELL = italic_L - ( italic_γ + ( italic_β + italic_β start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT ) italic_n start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT ) italic_N end_CELL end_ROW start_ROW start_CELL italic_N ( italic_t ) end_CELL start_CELL = italic_N start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 1 - roman_exp ( - roman_Γ italic_t ) ) end_CELL end_ROW (2)

where, Γ=γ+(β+βPA)nmaxΓ𝛾𝛽subscript𝛽𝑃𝐴subscript𝑛𝑚𝑎𝑥\Gamma=\gamma+\left(\beta+\beta_{PA}\right)n_{max}roman_Γ = italic_γ + ( italic_β + italic_β start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT ) italic_n start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT. In steady state,

NPANat=1γ+βnmaxγ+(β+βPA)nmaxsubscript𝑁𝑃𝐴subscript𝑁𝑎𝑡1𝛾𝛽subscript𝑛𝑚𝑎𝑥𝛾𝛽subscript𝛽𝑃𝐴subscript𝑛𝑚𝑎𝑥\frac{N_{PA}}{N_{at}}=1-\frac{\gamma+\beta n_{max}}{\gamma+\left(\beta+\beta_{% PA}\right)n_{max}}divide start_ARG italic_N start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT end_ARG start_ARG italic_N start_POSTSUBSCRIPT italic_a italic_t end_POSTSUBSCRIPT end_ARG = 1 - divide start_ARG italic_γ + italic_β italic_n start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT end_ARG start_ARG italic_γ + ( italic_β + italic_β start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT ) italic_n start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT end_ARG (3)

where, Natsubscript𝑁𝑎𝑡N_{at}italic_N start_POSTSUBSCRIPT italic_a italic_t end_POSTSUBSCRIPT is the steady state atom number without PA and NPAsubscript𝑁𝑃𝐴N_{PA}italic_N start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT is the fraction of atoms converted into PA. Even for a bright MOT with peak density, ρ08×1010similar-tosubscript𝜌08superscript1010\rho_{0}\sim 8\times 10^{10}italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∼ 8 × 10 start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPT cm33{}^{-3}start_FLOATSUPERSCRIPT - 3 end_FLOATSUPERSCRIPT this can be used to find the approximate value of βPAsubscript𝛽𝑃𝐴\beta_{PA}italic_β start_POSTSUBSCRIPT italic_P italic_A end_POSTSUBSCRIPT assuming a radiation trap model for MOT density [43].

Appendix B PA resonances in dark MOT

Table 1: VRS values for different PA resonances. VRS(off) is when the PA laser is off-resonant, and VRS(on) is when the PA laser is on-resonant.

Nearest PA

Resonance

VRS(off)

MHz

VRS(on)

MHz

pa1

0g,v=50superscriptsubscript0𝑔𝑣500_{g}^{-},v=500 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_v = 50

1g,v=164subscript1𝑔𝑣1641_{g},v=1641 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_v = 164

0u+,v=251superscriptsubscript0𝑢𝑣2510_{u}^{+},v=2510 start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , italic_v = 251

25.4125.4125.4125.41

±0.8plus-or-minus0.8\pm 0.8± 0.8

22.5322.5322.5322.53

±1.01plus-or-minus1.01\pm 1.01± 1.01

pa2

0g,v=51superscriptsubscript0𝑔𝑣510_{g}^{-},v=510 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_v = 51

1g,v=165subscript1𝑔𝑣1651_{g},v=1651 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_v = 165

0u+,v=252superscriptsubscript0𝑢𝑣2520_{u}^{+},v=2520 start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , italic_v = 252

24.4324.4324.4324.43

±0.54plus-or-minus0.54\pm 0.54± 0.54

20.9520.9520.9520.95

±1.44plus-or-minus1.44\pm 1.44± 1.44

pa3

0g,v=52superscriptsubscript0𝑔𝑣520_{g}^{-},v=520 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_v = 52

1g,v=166subscript1𝑔𝑣1661_{g},v=1661 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_v = 166

0u+,v=253superscriptsubscript0𝑢𝑣2530_{u}^{+},v=2530 start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , italic_v = 253

24.8224.8224.8224.82

±1.141plus-or-minus1.141\pm 1.141± 1.141

19.3319.3319.3319.33

±1.9plus-or-minus1.9\pm 1.9± 1.9

pa4

0g,v=53superscriptsubscript0𝑔𝑣530_{g}^{-},v=530 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_v = 53

1g,v=168subscript1𝑔𝑣1681_{g},v=1681 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_v = 168

0u+,v=255superscriptsubscript0𝑢𝑣2550_{u}^{+},v=2550 start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , italic_v = 255

22.1422.1422.1422.14

±1.19plus-or-minus1.19\pm 1.19± 1.19

15.1815.1815.1815.18

±1.66plus-or-minus1.66\pm 1.66± 1.66

pa5

0g,v=54superscriptsubscript0𝑔𝑣540_{g}^{-},v=540 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_v = 54

1g,v=169subscript1𝑔𝑣1691_{g},v=1691 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_v = 169

0u+,v=256superscriptsubscript0𝑢𝑣2560_{u}^{+},v=2560 start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , italic_v = 256

25.3825.3825.3825.38

±1.3plus-or-minus1.3\pm 1.3± 1.3

12.1312.1312.1312.13

±0.58plus-or-minus0.58\pm 0.58± 0.58

pa6

0g,v=58superscriptsubscript0𝑔𝑣580_{g}^{-},v=580 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_v = 58

1g,v=173subscript1𝑔𝑣1731_{g},v=1731 start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_v = 173

0u+,v=260superscriptsubscript0𝑢𝑣2600_{u}^{+},v=2600 start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , italic_v = 260

24.2824.2824.2824.28

±1.32plus-or-minus1.32\pm 1.32± 1.32

12.712.712.712.7

±0.46plus-or-minus0.46\pm 0.46± 0.46

The table B shows the values of VRS measured for different PA transitions in dark MOT. Here, VRS (off) represents the vacuum Rabi split when the PA laser is kept off-resonance, and VRS (on) is when the PA laser is kept on-resonance. The data is fitted with a theoretical expression for intra-cavity intensity derived from a two-level atom coupled to a cavity to extract VRS values. Each value given in the table is an average of 10 measurements, and the error bar is the standard deviation. As seen from the table B the average value of VRS is always smaller when PA resonance is addressed. The fraction of excited Rb22{}_{2}start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT molecules created can be calculated from 1(Nc;1/Nc;0)=1(VRS1/2VRS0)21-\left(N_{c;1}/N_{c;0}\right)=1-\left(VRS_{1}{}^{2}/VRS_{0}{}^{2}\right)1 - ( italic_N start_POSTSUBSCRIPT italic_c ; 1 end_POSTSUBSCRIPT / italic_N start_POSTSUBSCRIPT italic_c ; 0 end_POSTSUBSCRIPT ) = 1 - ( italic_V italic_R italic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT / italic_V italic_R italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT ), where VRS1𝑉𝑅subscript𝑆1VRS_{1}italic_V italic_R italic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is the split when PA is on resonance and VRS0𝑉𝑅subscript𝑆0VRS_{0}italic_V italic_R italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is corresponding off-resonance value. All the PA transitions are identified from ab initio calculation done in ref [31].

Table 1: VRS values for different PA resonances. VRS(off) is when the PA laser is off-resonant, and VRS(on) is when the PA laser is on-resonant.