arxiv: v4 [physics.atom-ph] 24 Mar 2014

Transcription

1 Magneto-Optical Trapping of Holmium Atoms J. Miao, J. Hostetter, G. Stratis, and M. Saffman Department of Physics, University of Wisconsin, 115 University Avenue, Madison, Wisconsin 5376 (Dated: July 11, 18) arxiv:11.156v [physics.atom-ph] Mar 1 We demonstrate sub-doppler laser cooling and magneto-optical trapping of the rare earth element Holmium. Atoms are loaded from an atomic beam source andcapturedin six-beamσ + σ molasses using a strong J = 15/ J = 17/ cycling transition at λ = 1.5 nm. Due to the small difference in hyperfine splittings and Landé g-factors in the lower and upper levels of the cooling transition the MOT is self-repumped without additional repump light, and deep sub-doppler cooling is achieved with the magnetic trap turned on. We measure the leakage out of the cycling transition to metastable states and find a branching ratio < 1 5 which is adequate for state resolved measurements on hyperfine encoded qubits. The magneto-optical trap (MOT) is a standard and very widely used tool in cold atom physics. To date MOT operation has been demonstrated for about 3 different neutral elements. Operation of a MOT on an open transition which allows for pumping into metastable states which are not cooled requires repump light to return atomic population to the levels participating in laser cooling. The rare earth lanthanides did not originally appear amenable to laser cooling due to the presence of an open f shell (except for Yb which was laser cooled in 1999[1]), and a correspondingly complex atomic structure. It was recognized by McClelland that the lanthanides are amenable to laser cooling due to the presence of cycling transitions between odd(even) parity ground states with angular momentum J and even(odd) parity excited states with angularmomentum J = J+1. MOT operation was first demonstrated with an open f shell lanthanide in 6 using Er[]. That demonstration was followed in 1 by demonstration of MOTs for Dy[3] and Tm[]. In this work we report on the operation of a Ho MOT loaded from an atomic beam with and without laser slowing[5]. The Doyle group has also recently demonstrated lanthanide MOTs including Ho using buffer gas precooling of an atomic beam [6]. Interest in laser cooling and trapping of lanthanide atoms is motivated by several topics in current research. Quantum degenerate bosonic and fermionic gases of Dy[7, 8], and Er[9, 1] open up new research directions due to the large magnetic moments Er (7 µ B ), Dy (1 µ B ) which provide for much stronger dipolar magnetic interactions between neutral atoms than are present in the more widely studied alkali gases. Ho has one stable isotope 165 Ho, which is bosonic with nuclear spin of I nuc = 7/ and a ground electronic configuration [Xe]f 11 6s. The ground state term is I o,j = 15 giving eight hyperfine levels F g =...11, and the magnetic moment is 9 µ B. Ho is distinguished by having 18 hyperfine-zeeman states, the largest number of any stable atomic isotope. This large number of states is of interest for collective encoding of multi-qubit quantum registers[11, 1]. Implementation of collective encoding will rely on Rydberg blockade interactions[13] in a dense atomic sample. The demonstration of sub-doppler laser Energy (cm -1 ) γ=. s -1 5 γ=17 s -1 γ=53 s nm 9/ 11/ 13/ 15/ 17/ 19/ 1/ J 8, 9 F 7 e=1 s f 11 6s6p 6 [Xe]f J=17/ 6s6p [Xe]f 1 5d6s 1 [Xe]f 11 6s even odd f 11 6s J=15/ F g=11 slower repumper FIG. 1. (color online) Energy levels of Ho. The upper levels of the cycling and cooling transitions are labeled 1-7. The diagram on the right shows the cooling, repumper, and slowing light used on the 1.5 nm transition to level e. TABLE I. Ho cycling and cooling transitions. The columns list the vacuum wavelength, natural linewidth, and Doppler cooling limit (T D = γ/k B). Transition λ(nm) γ/π(mhz) Doppler limit(µk) unknown unknown unknown a b,c unknown b,c d b,c b,c 78. Wavelengths are calculated from [1]. a Linewidth derived from oscillator strength value reported in [15], b Ref. [16], c Ref. [17], d Two electron jump transition. cooling, MOT operation, and optical cycling transitions with fractional leakage rates at the 1 5 level reported here are first steps towards Rydberg spectroscopy and qubit encoding in Ho atoms. The level diagram of Ho together with cycling transitions suitable for cooling is shown in Fig. 1. The ground 1 9 c cooling 7.67 GHz GHz 7.77

2 N (1 ) t (ms) N(t)/N() τ=1 ms τ=31 ms µm density (1 1 m -3 ) FIG.. Vacuum and laser cooling setup. The length of the apparatus from end to end is about 1 m. The slowing beam enters the vacuum chamber through a vertical window and is reflected from a Cu mirror to propagate towards the atomic source. state has J = 15/ with odd parity. There are eight dipole allowed transitions at wavelengths longer than 1 nm to levels with even parity and J = 17/. Since there are no odd parity levels with 15/ J 19/ in between the ground state and the upper levels of the six transitions labeled 1-6 these transitions are cycling. Relevant transition parameters are given in Table I. The 1.5 nm transition from the ground state to [Xe]f 11 6s6p( 1 P),J = 17/ at 361 cm 1 (labeled 9 in Fig. 1) is not cycling since there are odd parity levels above the ground state accessible by electric dipole decay from e. Nevertheless a closer examination of the odd parity levels below e with 15/ J 19/ reveals that almost all possible transitions have j > 1 on a single electron or are forbidden due to intercombination spin changes. The only transition which is allowed for single electron jumps is e f 11 5d6s( 1 D),J = 17/ at a transition energy of 75 cm 1. This transition is very weak due to the ω 3 factor in the expression for the radiative linewidth, with ω the transition frequency. We can roughly estimate the decay rate using hydrogenic orbitals which gives γ. s 1. The actual radial matrix elements are unknown so this value is not quantitatively correct but provides some guidance. To the extent that LS coupling is a good description for the level structure of Ho we expect that the cooling transition with g as an upper level will have small leakage since γ /γ 1 5. The experimental apparatus is shown in Fig.. A water cooled effusion cell with Ta crucible operated at T = 115 C provides a beam of Ho atoms with a mean velocity of 51 m/s. The atomic beam passes through a.5 cm diameter tube to prevent outflow of any melted Ho from the horizontally oriented effusion cell, and a.5 cm diameter aperture for differential pumping before entering the MOT chamber. Two ion pumps provide a base pressure of 1 9 mbar in the MOT chamber. A pair of electric coils provide a quadrupole mag c /π (MHz) FIG. 3. (color online) Number of trapped atoms (circles) and peak density (squares) as a function of cooling light detuning with the slower light turned on and the repump light turned off. Themagnetic fieldgradientwas.3t/m. Theinsetsshow an averaged fluorescence image of the trapped atoms from 1 exposures at -36 MHz detuning and population decay curves with exponential fits at -15 MHz (fast decay) and -5 MHz (slow decay). netic field with a gradient of up to. T/m (vertical axis) and. T/m (horizontal axis). The cooling beams were arranged in a standard 6 beam σ + σ configuration. The beams had Gaussian waists (1/e intensity radius) of. mm and total incident power of mw which was doubled by retroreflecting the beams. The cooling light was detuned by c 1.5γ from the F g = 11 F e = 1 cycling transition. Repump light was tuned to F g = 1 F e = 11 and overlapped with all MOT beams. In addition a circularly polarized slowing beam counterpropagating to the atomic beam was detuned by s = π 3 MHz from the F g = 11 F e = 1 transition. The slowing beam had a power of 1 mw and was focused to a waist of 1. mm. Measurements of the MOT atom number and density were made with an electron multiplying charge coupled device (EMCCD) camera using either fluorescence imaging or absorption imaging with an additional beam tuned to be resonant with the cooling transition. All laser beams were derived from a frequency doubled Ti:Sa laser (M Solstis with ECD-X) providing up to 1.5 W of 1.5 nm light. The laser was locked to the cooling transition using saturation spectroscopy in a hollow cathode lamp. The frequency and power level of the cooling, repump, and slowing light was controlled by acousto-optic modulators. The hyperfine energies shown in Fig. 1 were calculated from known values of the A and B constants for the ground state[18] and measured values for the excited state. We measured the hyperfine constants of the upper level of the cooling transition using modulation transfer spectroscopy in the hollow cathode lamp. Fits to our data gave A = 65.9±.3, B = 6± MHz. These values agree well with earlier measurements reported in [19]. Withtheslowingbeamturnedon, butnorepumplight, -3

3 3 we achieved a typical MOT population of N and atomic density of n a m 3 as shown in Fig. 3. Additional data taken with 9 mw of MOT light, beam waists of 1.1 cm, and 38 mw of slower light gave larger MOTs with up to N = 1 5 atoms. The atom number measurements were calibrated by integrating the detected EMCCD MOT image using the measured camera sensitivity to 1.5 nm light and an integration time of 5 ms. The measurement relies on knowing the rate of scattered photons per atom which we estimated by the two-level expression[] r = γρ ee = γ I T /I s 1+ c/γ +I T /I s. (1) Here ρ ee is the excited state fraction, I T is the total intensity of the six MOT beams and for the saturation intensity we use I s =.76 I sc. Here I sc = 61. W/m is the saturation intensity for the cycling transition F g = 11,M = 11 F e = 1,M = 1 and the factor of.76 = 3( 11+1)/( 11+3) accounts for averaging over Zeeman substates and the random light polarization in the MOT region. The scattered light was collected with a lens of numerical aperture.5 and imaged onto a calibrated electron multiplying charge coupled device camera which allowed us to deduce the atom number on the basis of the camera photoeletron counts. When the repumper was turned on tuned to the F g = 1 F e = 11 transition the atom number values in Fig. 3 increased by less than 1%. The negligible influence of the repump light is due to the fact that the cooling light also repumps population in F g = 1 much faster than the depumping rate out of F g = 11. The depumping rate due to Raman transitions via F e = 1 or 11 is calculated by averaging over M levels and light polarization, and accounting for the branching ratios of the fluorescence decay. We find r R = γ 3 [ ] F e=f g F e=f g 1 I T I sc.() In this expression F g = 11 is the F value for the lower levelofthecyclingtransition, Fe=F g = (F+1)e,F e + c, Fe=F g 1 = (F+1)e,(F 1) e + c are the detunings of the MOT light from F e = 11 and F e = 1, and we have assumed the light is far detuned so that we may replace factors of (1 + /γ + I T /I s ) 1 by γ /( ). The excited state hyperfine splittings shown in Fig. 1 are 1,11 = π 7.67 GHz and 1,1 = π 1.79 GHz. We find that for c = γ, the largest detuning we have used, r R = 8 s 1 and r R /r = There is also depumping due to leakage to metastable states. The rate for this process we estimate below to be negligible compared to the Raman rate. The total depumping rate is therefore < 1 s 1. This rate is balanced by the repumping rate for the Raman process F g = 1 F e = 11 F g = 11. The cooling light is detuned from the F g = 1 F e = 11 transition by 61 MHz+ c /π 65 MHz. The resulting repumping rate is about s 1. The ratio of the depumping to repumping rates implies that 99.9% of the population is in F g = 11. We verified that turning on a separate repumper increased the atom number by less than 1%, with the resolution limited by the stability of the MOT number measurement. Thus, even without a separate repumper, the atoms are pumped with better than 99% purity into the F g = 11 level. Further pumping into a specific Zeeman sublevel as a starting point for quantum state control experiments has not been demonstrated, but should be straightforward using an additional π polarized beam tuned to F g = 11 F e = 11 for pumping into M =, or a σ + polarized beam tuned to F g = 11 F e = 1 for pumping into M = 11. The effect ofthe repumper ismoresignificantwhen the MOT is operated without a slowing beam. In this case the MOT captures only a very small fraction of the low velocity tail of the atomic beam and the loading rate is reduced by almost three orders of magnitude, resulting in a very small MOT with N = atoms, again without a repumper. Turning on the repumper increases N by a factor of up to.5, depending on the detuning of the repumper, as shown in Fig.. The increase in atom number cannot be due to additional repumping since, as explained above, the cooling light acts as its own repumper. We instead interpret the data as being due to the repump light cooling and trapping additional atoms on the F g = 1 F e = 11 transition. The high temperature atomic source creates a beam which we assume to be uniformly distributed among ground Zeeman states. Since there are 3 Zeeman states in F g = 11 and 1 in F g = 1 we expect roughlyafactor ofincrease, which is in reasonable agreement with the data. No such increase is seen when the slower is on since the slower light is only effective for atoms in F g = 11 and the additional loading of unslowed F g = 1 atoms is negligible. As shown in Fig. 5 we observe deep sub-doppler cooling well below T D = 78 µk using σ + σ MOT beams. Data were acquired by observing MOT expansion and fitting to Gaussian density profiles. The theoretical dependence of the molasses temperature with respect to detuning is given by[1] T = T + γ I k B I s ( ) b a+ 1+ /γ (3) with a,b constants and T the low intensity temperature limit. The parametersa,b are known for low angularmomentum transitions[1], but have not been calculated for the F g = 11 F e = 1 Ho cooling transition. Using separate fit parameters for the horizontal and vertical temperatures Eq. (3) reproduces the observed dependence on detuning quite well. Cooling anisotropy has also been observed in Dy MOTs[] and may be due to the large magnetic moments present in rare earth atoms. Due to the near equality of the ground state F g = 11 and ex-

4 N/N(repump off) mw 1 mw 5 mw Fg=1 - Fe=11 1. P= r /π (MHz) FIG.. MOT population as a function of repumper detuning relative to the F g = 11 F e = 1 transition, and repumper power, with the slowing beam off. The population is normalized to the value with no repumper. The data stops at 6 MHz which was the limit of the AOM used for frequency shifting. T (µk) 6 vertical horizontal /π (MHz) FIG. 5. (Color online) MOT temperature with quadrupole magnetic field on versus detuning at saturation parameter I/I s = 1.1. The temperatures were extracted from fits to time of flight data for free expansion times up to 15 ms. The Doppler temperature is 78 µk. citedstatef e = 1g-factors,.8and.83,temperatures reaching ten times below the Doppler limit are achieved with relatively small detunings and with the quadrupole MOT field on during the entire cooling phase. In order to further study the internal state dynamics we loaded the MOT to steady state with the slowing beam on and then measured the decay lifetime with the slowing beam turned off, which greatly reduces the loading into the MOT, but the magnetic quadrupole field still on. The MOT population decay can be modeled with the rate equation[3] dn dt = R off ΓN βn. () Here R off is the loading rate with the slowing beam off, and the loss rate is Γ = (1 ρ ee )γ bg + ρ ee (γ la + γ ms ) where γ bg accounts for losses due to background collisions with ground state atoms, γ la is the loss rate due to light assisted collisions between optically excited and backgroundatoms[], γ ms is the loss rate due to scattering into metastable states, and β accounts for losses from collisions of trapped atoms. Leakage into metastable states may lead to atom loss since the metastable states do not interact with the cooling light and the fraction of atoms which are not magnetically trapped will drift away. The lowest energy odd parity states with J = 15/,13/,11/,9/ are fine structure levels of the ground term. They decay via magnetic dipole transitions which have been calculated in the LS coupling approximation. The resulting lifetimes indicated in Fig. 1 range from. to. s. Also the lowest energy even parity states with J = 17/,19/,1/ are expected to be long lived, but with lifetimes that are unknown. Measurements of the MOT number decay curve with the slowing beam turned off showed no indication of fast two-body loss (β ). The MOT lifetime was found to range from 1-3 ms as the detuning was scanned from 5 to 5 MHz. Representative decay curves are shown in Fig. 3. Fitting Eq. () to the data with ρ ee determined from Eq. (1) we found γ bg 1 ± 1 s 1 and γ ms +γ la 3.±1. s 1. The value for γ la is unknown so the value of 3. s 1 should be taken as an upper limit on the scattering rate into metastable loss states. This value is much larger than those reported for Yb[5] and Tm[] but is much smaller than the value reported for 165 Ho[6] (γ ms = 151 s 1 ). The Ho measurements in [6] were taken with a much higher background pressure, mainly He, which might contribute to the discrepancy. The branching ratio to metastable states is in any case very small γ ms /(γ +γ ms ) 1 6 so that loss out of the cycling transition is dominated by Raman transitions to other ground hyperfine levels with a branching ratio we estimatedabovetobe 1 5. Onthebasisofthesemeasurements we conclude that it should be possible to make high fidelity hyperfine state resolved measurements of Ho atoms, which will be important for qubit experiments. In conclusion, we demonstrated a magneto-optical trap of Ho atoms. Single frequency molasses light together with a slowing beam is sufficient to load a MOT with as many as 1 5 atoms at temperatures below 1 µk from an atomic beam source. The atoms are prepared in the F g = 11 level with probability 99% and fractional leakage rates out of the cooling transition are 1 5. Future work will explore the use of this cold sample as the starting point for quantum control experiments which take advantage of the large manifold of hyperfine states. This work was supported by NSF grant PHY and the University of Wisconsin graduate school. We are grateful to Jake Covey, Johannes Nipper, and Hannes Gorniaczyk for contributions at early stages of this experiment and we thank Boerge Hemmerling from the Doyle group for helpful discussions concerning measurements of the metastable decay rate.

5 5 [1] K. Honda, Y. Takahashi, T. Kuwamoto, M. Fujimoto, K. Toyoda, K. Ishikawa, and T. Yabuzaki, Phys. Rev. A 59, R93 (1999). [] J. J. McClelland and J. L. Hanssen, Phys. Rev. Lett. 96, 135 (6). [3] M. Lu, S. H. Youn, and B. L. Lev, Phys. Rev. Lett. 1, 631 (1). [] D. Sukachev, A. Sokolov, K. Chebakov, A. Akimov, S. Kanorsky, N. Kolachevsky, and V. Sorokin, Phys. Rev. A 8, 115(R) (1). [5] J. Miao, J. Hostetter, G. Stratis, and M. Saffman, Bull. Am. Phys. Soc. 58, U.5 (13). [6] B. Hemmerling, G. K. Drayna, E. Chae, A. Ravi, and J. M. Doyle, arxiv: (13). [7] M. Lu, N. Q. Burdick, S. H. Youn, and B. L. Lev, Phys. Rev. Lett. 17, 191 (11). [8] M. Lu, N. Q. Burdick, and B. L. Lev, Phys. Rev. Lett. 18, 1531 (1). [9] K. Aikawa, A. Frisch, M. Mark, S. Baier, A. Rietzler, R. Grimm, and F. Ferlaino, Phys. Rev. Lett. 18, 11 (1). [1] K. Aikawa, A. Frisch, M. Mark, S. Baier, R. Grimm, and F. Ferlaino, Phys. Rev. Lett. 11, 1 (1). [11] E. Brion, K. Mølmer, and M. Saffman, Phys. Rev. Lett. 99, 651 (7). [1] M. Saffman and K. Mølmer, Phys. Rev. A 78, 1336 (8). [13] M. Saffman, T. G. Walker, and K. Mølmer, Rev. Mod. Phys. 8, 313 (1). [1] A. Kramida, Y. Ralchenko, J. Reader, and NIST ASD Team, version 5.1 (13). [15] V. N. Gorshkov and V. A. Komarovskii, Opt. Spectrosk. 7, 631 (1979), [Opt. Spectrosc. (USSR) 7, 35 (1979)]. [16] E. A. D. Hartog, L. M. Wiese, and J. E. Lawler, J. Opt. Soc. Am. B 16, 78 (1999). [17] G. Nave, J. Opt. Soc. Am. B, 193 (3). [18] L. S. Goodman, H. Kopfermann, and K. Schlüpmann, Naturwiss. 9, 11 (196). [19] B. K. Newman, N. Brahms, Y. S. Au, C. Johnson, C. B. Connolly, J. M. Doyle, D. Kleppner, and T. J. Greytak, Phys. Rev. A 83, 1713 (11). [] H. J. Metcalf and P. van der Straten, Laser cooling and trapping (Springer-Verlag, New York, 1999). [1] J. Dalibard and C. Cohen-Tannoudji, J. Opt. Soc. Am. B 6, 3 (1989). [] S. H. Youn, M. Lu, and B. L. Lev, Phys. Rev. A 8, 33 (1). [3] D. Sesko, T. Walker, C. Monroe, A. Gallagher, and C. Wieman, Phys. Rev. Lett. 63, 961 (1989). [] J. E. Bjorkholm, Phys. Rev. A 38, 1599 (1988). [5] T. Loftus, J. R. Bochinski, R. Shivitz, and T. W. Mossberg, Phys. Rev. A 61, 511 ().