Current address: ]Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
Rapid generation of metastable helium Bose-Einstein condensates
Abstract
We report the realisation of Bose-Einstein condensation (BEC) of metastable helium atoms using an in-vacuum coil magnetic trap and a crossed beam optical dipole trap. A novel quadrupole-Ioffe configuration (QUIC) magnetic trap made from in-vacuum hollow copper tubes provides fast switching times while generating traps with a 10G bias, without compromising optical access. The bias enables in-trap 1D doppler cooling to be used, which is the only cooling stage between the magneto-optic trap (MOT) and the optical dipole trap. This allows direct transfer to the dipole trap without the need for any additional evaporative cooling in the magnetic trap. The entire experimental sequence takes 3.3 seconds, with essentially pure BECs observed with atoms after evaporative cooling in the dipole trap.
The experimental creation of Bose-Einstein condensates (BECs) of dilute weakly interacting gases of atoms Anderson et al. (1995); Bradley et al. (1995); Davis et al. (1995) has opened the possibility of exploring interesting phenomena of the quantum world on a macroscopic scale. Over the subsequent years the field has exploded, with BECs now used in diverse fields of quantum science, including quantum many-body systems Bloch et al. (2008), topological physics Cooper et al. (2019) and precision inertial measurements Geiger et al. (2020). While BEC experiments were initially limited to observing collective properties of the ensemble, a number of more recent detection techniques allow for the detection of individual atoms Ott (2016), opening up a broad new range of experimental possibilities. In experiments involving alkali or rare earth atoms, such single atom detection is usually performed via high resolution fluorescence imaging, either in-situ in optical lattices Bakr et al. (2009); Sherson et al. (2010) or after expansion from a trap Bücker et al. (2009); Bergschneider et al. (2018). However, such techniques usually have limitations on their spatial extent and the atom number able to be imaged.
An alternative method of single atom detection exploits the high internal energy of helium atoms trapped in the first atomic excited metastable state (He∗), which has 19.8eV of internal energy. The high internal energy allows direct detection of individual atoms with full 3D resolution using electronic detectors such as multi-channel plate and delay-line detectors (MCP-DLD) Vassen et al. (2012). This has opened up a wide range of exciting experimental possibilities, allowing quantum-optics equivalent demonstrations with atoms of iconic experiments such as the Hanbury Brown-Twiss (HBT) effect Schellekens et al. (2005), Wheeler’s delayed choice Manning et al. (2015), the Hong-Ou-Mandel effect Lopes et al. (2015), ghost imaging Khakimov R. I. et al. (2016) and the measurement of Bell correlations Shin et al. (2019). More than just replicating quantum optics though, the atomic interactions allow effects to be seen that are not possible with photons, such as Fermionic anti-bunching Jeltes et al. (2007), quantum depletion Chang et al. (2016) and strongly interacting lattice physics Cayla et al. (2020).
Crucial to many of these experiments are the measurement of HBT-style correlation functions - the key observable enabled by single atom detection. However, such correlation functions require large amounts of data, especially if higher order correlations are being measured Dall R. G. et al. (2013); Hodgman et al. (2017), often needing 10,000-100,000 individual experimental runs for a single experiment. This has led to constant improvements across generations of experimental He∗ apparatuses aiming for ever shorter experimental sequences. Most He∗ experiments Robert et al. (2001); Pereira Dos Santos et al. (2001); Tychkov et al. (2006); Dall and Truscott (2007); Keller et al. (2014) use magnetic traps, where the sequence length is limited by the slow thermalisation rates in magnetic traps that have relatively weak trap frequencies. While using liquid helium in the source stage can substantially reduce the sequence length Dall and Truscott (2007), it is less practical for everyday operation. The tight traps provided by dipole traps can overcome this, either in a hybrid combination with a magnetic trap Flores et al. (2015) or in a stand-alone cross beam configuration Bouton et al. (2015). However, due to the limited depth of the dipole trap, to enable an efficient transfer the atoms need to first be cooled, either optically via complex additional cooling schemes such as grey molasses and/or evaporatively in a magnetic trap Bouton et al. (2015).
In this work, we report on the construction of a He∗ BEC machine capable of creating condensates in a simplified, rapid sequence that takes only 3.3s in total. An initial stage of 1D Doppler cooling Tychkov et al. (2006) is performed in a biased magnetic trap, constructed using a simplified in-vacuum water-cooled coil geometry to produce a biased quadrupole trap. This lowers the temperature enough to permit direct loading of a crossed dipole trap, which allows fast and efficient evaporative cooling. The final BECs are produced with atoms. As well as the relatively small size and minimal number of coil turns, which allow fast switch off times, a major advantage of this magnetic trap design is that the small size allows good optical access to the atomic cloud. This will make the apparatus ideal for future optical lattice experiments Cayla et al. (2020).

A number of magnetic trap designs are used in BEC experiments, with some commonly used ones including the cloverleaf design Mewes et al. (1996), the Quadrupole-Ioffe configuration (QUIC) trap Esslinger et al. (1998) and atomic chip traps Fortágh and Zimmermann (2007). Any design has strengths and weaknesses, and inevitably involves compromises and trade-offs between factors such as confinement, depth, stability, switch off time, physical size, current required and heat dissipation. In our magnetic trap setup, no evaporation is implemented; we only perform 1D Doppler cooling in the trap. Hence we require a large enough bias to split the atomic energy levels (as well as prevent Majorana and Penning losses), but are not so concerned about the usual consideration of having a tight confinement, as all evaporation is conducted in the dipole trap. The trap is in a simple QUIC configuration, consisting of two 5 turn anti-Helmhotz (AH) coils of 40 mm diameter centred on the axis and an 8 turn bias coil of 20 mm diameter centred on the axis, as shown in Fig. 1. These coils are made from 2 mm hollow copper tubing (internal diameter 0.4 mm) mounted in-vacuum to ensure they are close to the atoms. Cooling is provided by pumping chilled water through the tubes, with each tube attached to hollowed out inch feedthroughs sealed with vacuum compatible solder sol . By driving A through the AH coils and A though the bias, we were able to produce a trap of mK depth and G bias field (see Fig. 2) in the trap centre, which is offset 6mm vertically along the axis from the centre of the AH coils. The radial and axial trap frequencies are Hz in the radial axes and Hz respectively.
A cold source of He∗ atoms suitable for loading into a MOT is generated from our He∗ beamline, shown schematically in Fig. 1. Helium atoms are excited into the metastable state via a hollow cathode DC discharge source Swansson et al. (2004), which is cryogenically cooled using liquid nitrogen. The expanding atomic beam is then collimated using an optical collimation stage featuring 4 beams propagating perpendicular to the He∗ atoms with an intensity ( mW/cm2) and detuned by (MHz) from resonance. Note that all detunings for laser cooling and imaging in this work are relative to the transition, with a negative sign indicating red detuning. All laser cooling and imaging beams are sourced from a home built external gain-chip laser Shin et al. (2016) with linewidth 100kHz, which is used to seed a 5W fibre amplifier. Individual frequency shifts are provided by acousto-optic modulators in each cooling beam. A Zeeman slower is then employed to reduce the velocity of the atoms to m/s, featuring a beam of intensity and detuning, counter-propagating with the direction of the atomic beam. A spatially varying field to keep the atoms on resonance while they slow down is generated by two separate coils with variable windings along their length and a zero field in the centre. These atoms are then trapped and cooled by the 1st MOT, with the 3 counter-propagating MOT beams having intensity (horizontal beam), and intensity (vertical beams), all with detuning. However, the background pressure from the source is too high for this chamber to be used to form a BEC. The horizontal MOT retro mirror located inside the vacuum chamber has a mm diameter hole in the centre, which enables atoms to leak out of the MOT and into the UHV MOT chamber, where the pressure is Torr. To assist this process, an additional ‘push’ beam is added ( intensity and blue detuned), forming a low velocity intense source scheme (LVIS) Dall and Truscott (2007). The flux from the LVIS, as measured 10 cm beyond the location of the MOT, is atoms/s, measured using a Faraday cup on a rotation stage and a picoammeter.

The second MOT consists of a quadrupole magnetic field generated by our in-vacuum AH coils with a gradient along the tight axis of G/cm, along with three pairs of counter-propagating laser beams with detunings and intensities of , and . In 1s we load a cloud with atoms at a temperature of mK. The atom number is measured using saturated fluorescence on an InGaAs photodiode Dall and Truscott (2007). To measure the temperature we switch off the MOT and allow the atoms to fall mm onto a detector comprising an mm diameter stacked pair of multi-channel plates (MCP). By using pick-off electronics we extract a signal pulse from the charge depletion of the plates that each atom causes. After fast amplification and discrimination of the pulses, they are recorded using a digital counter, which provides an integrated 1D time-of-flight (TOF) distribution of the atoms. The TOF distribution is fitted to a Maxwell-Boltzmann distribution, given by Yavin et al. (2002); Dall and Truscott (2007)
(1) |
where , , is the mass of a 4He atom, is Boltzmann’s constant and is the falling distance. From this distribution we extract the temperature , with the only other free fit parameter being the amplitude .
Following loading, the MOT is then compressed by ramping the magnetic field down to G/cm in ms, while simultaneously the laser detuning is ramped up to and the intensity of the MOT beams are reduced by two orders of magnitude. The MOT beams are then switched off and the atoms in the state remain trapped in the quadrupole magnetic field of the AH coils. To improve the trap depth and atom density we tighten the quadrupole field to G/cm in s. The bias coil current is then ramped up in ms to form a QUIC trap with a G field offset, removing the magnetic field zero in the quadrupole trap. This configuration provides weak trapping frequencies of Hz in the radial (x-y) direction and Hz in the axial (z) direction. At this stage we have atoms at mK, with the lower temperature mostly due to the reduced magnetic trap depth compared to the MOT. The atom number and temperature in the magnetic trap are measured via absorption imaging, with the imaging beam along the x axis being imaged on an InGaAs CCD camera. Mechanical flipper mirrors allow us to swap between imaging and MOT beams on the same axis, although this prevents absorption imaging during the operation of the MOT.
To cool the cloud further it is illuminated with a 1D Doppler cooling beam Tychkov et al. (2006) of intensity , aligned vertically along the bias field axis, circularly polarised to drive the transition and red detuned by . This Doppler cooling stage cools the cloud in ms, after which we have atoms at K. This is sufficiently cold to directly transfer the atoms from the magnetic trap to the crossed optical dipole trap.

The final trapping stage is a crossed optical dipole trap (CODT), formed from two intersecting laser beams far detuned from the 1083nm helium resonance at nm wavelength, from a 100kHz linewidth seed laser that is amplified to W of power. The power is distributed between the two beams such that after feedback loops to regulate the intensity and losses due to AOM and fibre coupling the first beam (aligned along the axis) has up to 6W and the second beam (aligned from the axis in the plane) has up to 1W of power. These beams are focused down to Gaussian waists of m and m, respectively, at the location of the atoms. The two beams intersect at the centre of the G biased QUIC trap, forming a CODT with trapping frequencies kHz and a trap depth of K at full power.
Since the 1D Doppler cooling process is more efficient at higher densities Tychkov et al. (2006), we observe that the Doppler cooling is more effective with the dipole beams on, in addition to the magnetic trap. Here the dipole beams mostly serve to increase the density of the atomic cloud, making the cooling process more efficient and reducing the final temperature after 1D Doppler cooling. In our experiment, this improves the final atom number by . Hence the dipole beams are linearly increased to full power over 100 ms prior to the start of the 1D Doppler cooling stage. To transfer the atoms from the QUIC to the dipole trap, after the 1D Doppler cooling stage the trap is linearly decreased in ms to a non-biased trap with 1.1G/cm gradient along the tight axis, before being abruptly turned off with a FET switch. To preserve the quantisation axis of the atoms and prevent any spin flips during this switch off, an additional large (50 cm diameter square) coil mounted above the vacuum chamber was used to generate a uniform magnetic field of 1.6 G aligned along the direction. This field is switched on before the current through the QUIC trap starts ramping down. However, once in the dipole trap, it was found that the background magnetic field (dominated by the component of the Earth’s magnetic field) was sufficient to preserve the quantisation, so the extra quantisation coil was switched off 700 ms after transfer to the CODT. Immediately after transfer we measure atoms at K with a phase space density of in the CODT, an excellent starting point for runaway evaporative cooling.
Evaporation is performed in the CODT by lowering the trap depth through reducing the power of the dipole beams exponentially over 1 second. As the hotter atoms escape and the remaining atoms rethermalise, this reduces the temperature and enhances the phase-space density. The evaporation duration is set by the relatively high trap frequencies (and thus high rethermalisation rates). Unlike in a previous similar He∗ experiment Bouton et al. (2015), we did not observe any improvement in the evaporation process by adding a gradient to cancel the weak components of the dipole potentials along the beam propagation directions.
To determine the important parameters of the gas in the CODT, we switch off the dipole trap and allow the atoms to fall under gravity onto the MCP to measure the 1D TOF profile along the z axis, integrated in the x-y dimensions. Examples of the resulting integrated 1D TOF profiles are shown in Fig. 3. Above the condensation critical temperature (note that refers to the non-interacting critical temperature, while is the critical temperature corrected for interactions Dalfovo et al. (1999)), the profile is fitted to the same Maxwell-Boltzmann distribution as in Eqn. 1 to yield the temperature . Below we fit with a bi-modal distribution consisting of the same Maxwell-Boltzmann distribution plus an inverted parabola to represent the far-field Thomas-Fermi density profile of the BEC Dalfovo et al. (1999). The half width of this parabola is the far field Thomas-Fermi radius , from which we can extract the chemical potential from TOF . From the chemical potential we can extract the number of atoms in the condensate , where is the He∗ wave scattering length. The number of atoms in the thermal cloud is given by , for . To extract the atom number above , we sum the integrated counts, correct for the fraction of the cloud at temperature that will hit the detector and then scale by an effective quantum efficiency of the detector . We vary until the thermal atom number matches for the clouds just above and below . Note that non-linearities in at high fluxes mean this is not accurate for .

The above information, along with the total atom number and , is then combined to produce Fig. 4. This shows the condensate fraction vs the reduced temperature for our sequence as we evaporate through the transition temperature. We cross at with atoms, and by evaporating further we are able to produce BECs with no discernible thermal fraction (not shown on the graph) with atoms.
We demonstrate a comparatively simple and short method for producing Bose Einstein condensates of He∗ atoms. By using a magnetic trap with a 10G bias generated from an in vacuum set of coils, we are able to cool the cloud below in the magnetic trap via 1D Doppler cooling. The degenerate state of metastable helium atoms is accomplished via direct evaporative cooling in a crossed optical dipole trap. We cross the condensation critical temperature with atoms at K and by evaporating further are able to generate essentially pure BECs with atoms. The entire sequence to produce a BEC only takes 3.3 sec. This will provide an excellent starting point for a range of experiments with He∗ BECs, such as many-body correlation experiments Hodgman et al. (2017), probes of quantum non-locality Shin et al. (2019) or strongly interacting optical lattice physics Cayla et al. (2020).
Acknowledgements.
The authors would like to thank David Clement and Robert Dall for helpful discussions, and Colin Dedman, Ross Tranter and Liangwei Wang for technical assistance. This work was supported through Australian Research Council (ARC) Discovery Project grants DP160102337 and DP190103021. SSH was supported by ARC Discovery Early Career Researcher Award DE150100315.References
- Anderson et al. (1995) M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Science 269, 198 (1995).
- Bradley et al. (1995) C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. Hulet, Phys. Rev. Lett. 75, 1687 (1995).
- Davis et al. (1995) K. B. Davis, M. O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995).
- Bloch et al. (2008) I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008).
- Cooper et al. (2019) N. R. Cooper, J. Dalibard, and I. B. Spielman, Rev. Mod. Phys. 91, 015005 (2019).
- Geiger et al. (2020) R. Geiger, A. Landragin, S. Merlet, and F. Pereira Dos Santos, AVS Quantum Science 2, 024702 (2020).
- Ott (2016) H. Ott, Rep. Prog. Phys. 79, 054401 (2016).
- Bakr et al. (2009) W. S. Bakr, J. I. Gillen, A. Peng, S. Fölling, and M. Greiner, Nature 462, 74 (2009).
- Sherson et al. (2010) J. F. Sherson, C. Weitenberg, M. Endres, M. Cheneau, I. Bloch, and S. Kuhr, Nature 467, 68 (2010).
- Bücker et al. (2009) R. Bücker, A. Perrin, S. Manz, T. Betz, C. Koller, T. Plisson, J. Rottmann, T. Schumm, and J. Schmiedmayer, New Journal of Physics 11, 103039 (2009).
- Bergschneider et al. (2018) A. Bergschneider, V. M. Klinkhamer, J. H. Becher, R. Klemt, G. Zürn, P. M. Preiss, and S. Jochim, Phys. Rev. A 97, 063613 (2018).
- Vassen et al. (2012) W. Vassen, C. Cohen-Tannoudji, M. Leduc, D. Boiron, C. I. Westbrook, A. Truscott, K. Baldwin, G. Birkl, P. Cancio, and M. Trippenbach, Rev. Mod. Phys. 84, 175 (2012).
- Schellekens et al. (2005) M. Schellekens, R. Hoppeler, A. Perrin, J. V. Gomes, D. Boiron, A. Aspect, and C. I. Westbrook, Science 310, 648 (2005).
- Manning et al. (2015) A. G. Manning, R. I. Khakimov, R. G. Dall, and A. G. Truscott, Nature Phys. 11, 539 (2015).
- Lopes et al. (2015) R. Lopes, A. Imanaliev, A. Aspect, M. Cheneau, D. Boiron, and C. I. Westbrook, Nature 520, 66 (2015).
- Khakimov R. I. et al. (2016) Khakimov R. I., Henson B. M., Shin D. K., Hodgman S. S., Dall R. G., Baldwin K. G. H., and Truscott A. G., Nature 540, 100 (2016).
- Shin et al. (2019) D. K. Shin, B. M. Henson, S. S. Hodgman, T. Wasak, J. Chwedenczuk, and A. G. Truscott, Nature Communications 10, 4447 (2019).
- Jeltes et al. (2007) T. Jeltes, J. M. McNamara, W. Hogervorst, W. Vassen, V. Krachmalnicoff, M. Schellekens, A. Perrin, H. Chang, D. Boiron, A. Aspect, and C. I. Westbrook, Nature 445, 402 (2007).
- Chang et al. (2016) R. Chang, Q. Bouton, H. Cayla, C. Qu, A. Aspect, C. I. Westbrook, and D. Clément, Phys. Rev. Lett. 117, 235303 (2016).
- Cayla et al. (2020) H. Cayla, S. Butera, C. Carcy, A. Tenart, G. Hercé, M. Mancini, A. Aspect, I. Carusotto, and D. Clément, Phys. Rev. Lett. 125, 165301 (2020).
- Dall R. G. et al. (2013) Dall R. G., Manning A. G., Hodgman S. S., RuGway Wu, Kheruntsyan K. V., and Truscott A. G., Nature Phys. 9, 341 (2013).
- Hodgman et al. (2017) S. S. Hodgman, R. I. Khakimov, R. J. Lewis-Swan, A. G. Truscott, and K. V. Kheruntsyan, Phys Rev Lett 118, 240402 (2017).
- Robert et al. (2001) A. Robert, O. Sirjean, A. Browaeys, J. Poupard, S. Nowak, D. Boiron, C. I. Westbrook, and A. Aspect, Science 292, 461 (2001).
- Pereira Dos Santos et al. (2001) F. Pereira Dos Santos, J. Léonard, J. Wang, C. J. Barrelet, F. Perales, E. Rasel, C. S. Unnikrishnan, M. Leduc, and C. Cohen-Tannoudji, Phys. Rev. Lett. 86, 3459 (2001).
- Tychkov et al. (2006) A. S. Tychkov, T. Jeltes, J. M. McNamara, P. J. J. Tol, N. Herschbach, W. Hogervorst, and W. Vassen, Phys. Rev. A 73, 031603 (2006).
- Dall and Truscott (2007) R. Dall and A. Truscott, Opt. Commun. 270, 255 (2007).
- Keller et al. (2014) M. Keller, M. Kotyrba, F. Leupold, M. Singh, M. Ebner, and A. Zeilinger, Phys. Rev. A 90, 063607 (2014).
- Flores et al. (2015) A. S. Flores, H. P. Mishra, W. Vassen, and S. Knoop, Appl. Phys. B 121, 391 (2015).
- Bouton et al. (2015) Q. Bouton, R. Chang, A. L. Hoendervanger, F. Nogrette, A. Aspect, C. I. Westbrook, and D. Clément, Phys. Rev. A 91, 061402 (2015).
- Mewes et al. (1996) M.-O. Mewes, M. R. Andrews, N. J. van Druten, D. M. Kurn, D. S. Durfee, and W. Ketterle, Phys. Rev. Lett. 77, 416 (1996).
- Esslinger et al. (1998) T. Esslinger, I. Bloch, and T. W. Hänsch, Phys. Rev. A 58, R2664 (1998).
- Fortágh and Zimmermann (2007) J. Fortágh and C. Zimmermann, Rev. Mod. Phys. 79, 235 (2007).
- (33) Note that standard solder and flux was also found to work, although required an initial period under vacuum to remove the residual flux before being installed in the final UHV chamber. This prevented out-gassing from the non-vacuum compatible flux, was sufficient to still be able to reach UHV pressures.
- Swansson et al. (2004) J. Swansson, K. Baldwin, M. Hoogerland, A. Truscott, and S. Buckman, Appl. Phys. B 79, 485 (2004).
- Shin et al. (2016) D. K. Shin, B. M. Henson, R. I. Khakimov, J. A. Ross, C. J. Dedman, S. S. Hodgman, K. G. H. Baldwin, and A. G. Truscott, Opt. Express 24, 27403 (2016).
- Yavin et al. (2002) I. Yavin, M. Weel, A. Andreyuk, and A. Kumarakrishnan, American Journal of Physics 70, 149 (2002).
- Dalfovo et al. (1999) F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999).
- (38) While this approximation for the chemical potential from the far field is only valid for condensates in cigar shaped traps, numerical solutions of the equations of motion following Castin and Dum (1996) show that it gives at most an error in for the trap geometries in our system, which is within our experimental uncertainties.
- Castin and Dum (1996) Y. Castin and R. Dum, Phys. Rev. Lett. 77, 5315 (1996).