# Creation of matter wave Bessel beams

###### Abstract

Bessel beams are plane waves with amplitude profiles described by Bessel functions. They are important because of their property of limited diffraction and their capacity to carry orbital angular momentum. Here we report the creation of a Bessel beam of de Broglie matter waves. The Bessel beam is produced by the free evolution of a thin toroidal atomic Bose-Einstein condensate (BEC) which has been set into rotational motion. By attempting to stir it at different rotation rates, we show that the toroidal BEC can only be made to rotate at discrete, equally-spaced frequencies, demonstrating that circulation is quantized in atomic BECs. The method used here to generate matter wave Bessel beams with a âPainted Potentialâ can be viewed as a form of wavefunction engineering which might be extended to implement arbitrary cold atom matter wave holography.

^{†}

^{†}preprint: LA-UR-13-26796

Bessel beams - plane waves with amplitude profiles described by Bessel functions â are important because of their property of limited diffraction Durnin1987 () and their capacity to carry orbital angular momentum Volke-Sepulveda2002 (). These characteristics have already been exploited to improve the performance of interferometry Fortin2004 (), microlithography Erdelyi1997 (), and imaging Ok2013 (). Bessel beams have also been used for micromanipulation Arlt2001 (), where it has been proposed that they might create âtractor beamsâ Chen2011 (). To date, Bessel beams have been realized with light Durnin1987 (); McGloin2005 () and with sound Lu1990 (), and related Gauss-Laguerre beams with orbital angular momentum have also been realized with electrons Uchida2010 (); Verbeeck2010 (); McMorran2011 (). Here we report the first creation of a Bessel beam of de Broglie matter waves. The Bessel beam is produced by the free evolution of a thin toroidal atomic Bose-Einstein condensate (BEC) which has been set into rotational motion. By attempting to stir it at different rotation rates, we show that the toroidal BEC can only be made to rotate at discrete, equally-spaced frequencies. This is a direct demonstration that circulation is quantized in atomic BECs. It also dictates that atoms in the Bessel beam must carry quantized orbital angular momentum, and it means that the preparation of beams with this method does not require precise control of the initial BEC rotation rate. We show that the results of the experiment are consistent with a simple theory. The method used here to generate matter wave Bessel beams with a âPainted Potentialâ Henderson2009 () can be viewed as a form of wavefunction engineering which might be extended to implement arbitrary cold atom matter wave holography.

The useful properties of Bessel beams follow from the result that if is the th Bessel function of the first kind, then . Consequently satisfies the Helmholtz equation and so a plane wave with this amplitude profile and phase winding propagates without any change in form. While this perfectly non-diffracting plane wave has infinite transverse extent, practical implementations of these waves are still useful because the effect of finite extent is to limit the non-diffracting region to finite length Durnin1987 (). In the particular case of quantum mechanics, (with being a constant) is a diffraction-free solution of the free Schrödinger equation for a particle of mass . To make experimentally an approximation to this wavefunction, we note that at long times the solution of the 2D free Schrödinger equation has a “Fraunhofer diffraction limit” solution

(1) |

where is a constant, is the 2D momentum space wavefunction at time , and the position vector . Since a spatial wavefunction in the form of a thin torus of radius with phase winding number , , has , free expansion of such a wavefunction will create an approximation to the desired Bessel wave.

Our experiment follows this recipe with a toroidal atomic Bose-Einstein condensate (Fig. 1). In this system the requirement for a single integer phase winding number is guaranteed by the quantized rotational response that follows from the superfluid being described by a complex-valued order parameter . The associated superfluid velocity is , where is the phase of and is the mass of the particles forming the superfluid current. The line integral known as the circulation then takes on quantized values,

(2) |

because the requirement that the order parameter be single-valued restricts the phase change around the path to , where is the (integer) winding number. For a thin toroid of radius these allowed states have quantized angular velocity , where , and quantized orbital angular momentum per atom. In addition, these states are also metastable persistent current states with lifetimes of several seconds Ryu2007 ().

The discussion so far has assumed a thin toroidal BEC obeying the linear Schrödinger equation. However a real system is 3D, and the strong interatomic interactions lead to dynamics described by the nonlinear Gross-Pitaevskii equation (GPE). Fortunately we find from numerical solutions of the GPE that these factors should not cause a serious deviation from ideal Bessel evolution. Fig. 1, computed for the parameters of our experiment, shows the free expansion of toroidal BECs with and . The toroidal condensate initially expands in all directions (ms) until the wavefunction reaches the center of the toroid (ms) where the interference is either constructive (), forming a central peak, or destructive (), forming a central hole. The BEC then evolves into a long thin condensate with concentric layers structured in accordance with the relevant Bessel function (ms). Fig. 2 shows that the column density along the axial direction has fringes that match up well with the corresponding Bessel function, the main difference being that the fringe amplitude falls off with radius faster for the GPE solution. It is clear from the difference between the black and green curves that interactions only slightly expand the fringe pattern. The main effect of interactions is seen in the axial direction, where they drive the rapid expansion that leads to the elongated beam-like object depicted in Fig. 1.

Our experimental procedure to create rotating thin toroidal BECs is conceptually straightforward: we make an Rb condensate in the ground state by evaporating into a rotating toroidal trap, and we then measure the resulting superfluid velocity by releasing the condensate and imaging it after a period of free expansion. The rotating toroidal traps are realized using the ”Painted Potential” technique Henderson2009 () for creating condensates in arbitrary time-averaged optical dipole potentials formed by painting with a rapidly moving tweezer beam on top of a tight horizontal sheet trap. The technique can realize and animate any 2D potential that can be drawn on a sheet of paper. The trap diameter is 19m, for which rad/s, and the full width at half-maximums (FWHMs) of the BEC density in the axial and radial directions are 2m and 5m respectively for the typical atom number of 40,000.

For the experiments reported here, the system of Ref. Henderson2009 () has been improved in three ways. First, the horizontal light sheet is now more uniform because it is created by scanning a spherically focused beam instead of the static cylindrically-focused beam used previously. Second, the intensity of the tweezer beam which paints the potential is now modulated synchronously with the painting scan to improve the flatness of the potential. Third, separate lenses are used for focusing the tweezer beam and for imaging the trapped condensate, delivering better performance in both areas. The tweezer beam (waist =m and wavelength nm) paints at kHz scan frequency a m diameter torus with trap depth nK and radial trapping frequency Hz. The tweezer beam intensity is modulated so as to create a rotating barrier of height 45 nK and FWHM thickness m (the scheme was inspired by the proposal of Brand2001 (), although our torus is too thin to easily support defects inside the ring). The horizontal sheet optical dipole trap providing the vertical confinement has trapping frequency 400 Hz. The experimental timing is as follows. The optical dipole potential, superimposed on the cold atoms in the magnetic trap, is ramped up over 100 ms, at which point the magnetic trap is turned off. After a 200 ms equilibration time the trapped atoms are evaporated to condensation by lowering the sheet intensity over 3 seconds, forming a condensate which is rotating with the trap. A repulsive barrier with height greater than the chemical potential is initially painted across the toroid to force the cold gas to rotate with the trap as it condenses. After the BEC has formed in the rotating trap, the barrier is lowered adiabatically over 200 ms, transforming the potential into toroidal form. After 100 ms of equilibration the trapping potential is turned off and the condensate is imaged in absorption along its long axis following 34 ms of free expansion.

Fig. 3(a) and (b) show time of flight (TOF) images of the column density obtained in this way, exhibiting concentric ring patterns that agree qualitatively with the Bessel function intensity distributions predicted by Eq. (1) and shown in Fig. 3(c). Imaging these beam-like optically thick condensates is non-trivial because diffraction and refraction alter the intensity profile of the probe beam as it propagates through the condensate. We have developed a model of this process which predicts the simulated absorption images shown in Fig. 3(d) for our conditions. Modeling of the experimental images must take into account absorption, lensing, and diffraction of the probe beam as it propagates through the condensate. Our current model for this process proceeds as follows: First, the simulation of the free expansion dynamics is done in two stages: the GPE is solved numerically in polar coordinates until a time where the interaction energy of the condensate has become negligible, after which the wavefunction is projected onto a Cartesian mesh and propagated directly to any desired time using a decomposition into plane waves. Second, an appropriate two-level atom susceptibility is applied to the corresponding density distribution to model the medium through which the probe beam will propagate. A plane wave propagated through that medium using the paraxial wave equation yields the electric field distribution of the imaging beam at the exit of the condensate. That field is then propagated numerically through free space to the focal plane of the imaging system to obtain a simulation of the image seen by the camera. Full details of this procedure and the complete results will be published elsewhere, but the important result is that there exists a focus plane for the imaging system where absorption images appear very similar to the condensate column density. One notable deviation from the exact column density is that diffraction following absorption on entry to the condensate causes the intensity of the probe beam to actually increase above its input value Boshier1982 (), explaining the negative absorption seen at the center of many of the images in Figs. 3(a) and (b) with . The agreement between these experimental images and the corresponding predictions of the imaging model [Fig. 3(d)] is qualitatively good. The imaging model used here makes the paraxial approximation and it also assumes perfect imaging optics. We expect that a more sophisticated imaging model now under development will yield quantitative agreement between model and data, and also explain the apparent correlation between fringe radius and atom number seen in Fig. 4(b).

The simple classical argument that a zero-temperature toroidal cloud of gas rotating with speed would evolve asymptotically into a torus of radius after expansion time implies that central hole radius is a measure of rotational speed (the argument carries over to the quantum domain because the position of the first maximum of varies nearly linearly with ). Fig. 3(a) shows six successive TOF images taken with the barrier rotating at 2 Hz. It is clear that these images show two different hole sizes - corresponding to two different rotational speeds - illustrating directly that while the outcome of the experiment is not always the same the rotational speed of the toroidal BEC is quantized and the system always has integer phase winding, as predicted by Eq. (2). Analysis shows that images 2 and 3 correspond to winding number , while the other four images correspond to . Repeating the experiment for a range of stirring frequencies allows generation of states from up to , as seen in Fig. 3(b). We note that discrete rotational states in a toroidal BEC have been seen by us Ryu2010 () and by others Moulder2012 (); Wright2013 (), but the present work is the first to compare experimental results with theoretical predictions. Further, the quantized rotation seen here results from creating the condensate in a (non-quantized) rotating trap, so the experiment is a BEC analog of the famous Hess-Fairbank rotating bucket experiment Hess1967 ().

We fit the central fringe of each TOF image to a circular function with Gaussian radial profile to estimate the fringe diameter and the peak absorption on that fringe. Fig. 4(a) is a histogram of the fringe diameter measurements for TOF images exhibiting a central hole. The clear peaks in the distribution confirm that circulation in the torus is quantized, with the peaks corresponding to successive winding numbers from to . While some of the finite peak width in the histogram is a result of small imperfections in the potential, a scatter plot of peak absorption versus fringe radius [Fig. 4(b)] shows that most of the apparent peak width is associated with a correlation between fitted peak absorption (proportional to atom number for each n) and fitted fringe diameter. The red lines added to Fig. 4(b) to guide the eye show that an extrapolation of fringe diameter to the limit of zero absorption is in good qualitative agreement with the simple non-interacting Bessel function prediction for the radius of the first fringe, shown as green lines corresponding to the first maximum in of . We believe that the correlation between fringe radius and atom number, which is larger than the interaction effects seen in Fig. 2, is an imaging artifact (see text).

For the same data set Fig. 4(c) shows the average winding number as a function of the rotation frequency of the trap. The procedure of evaporating into the rotating trap followed by a slow lowering of the barrier should find the thermodynamic ground state of the system. For a trap rotating at angular frequency about the -axis, that is the state which minimizes the energy associated with the effective Hamiltonian in the rotating frame, where is the hamiltonian for the non-rotating trap and is the component of the angular momentum. In our case of a thin toroid, the energy is therefore minimized when the winding number takes the integer value closest to , giving rise to the equally spaced steps in the theory prediction shown as a red line on Fig. 4(c) (preliminary results from an improved apparatus exhibit these steps and suggest that the method can set the winding number deterministically with fidelity in excess of ).

Fig. 1 shows that the freely-expanding toroidal condensate eventually forms a very elongated object analogous to a pulsed beam: after 34ms of free evolution the BEC has fallen mm under gravity and grown from 2m to 300m in length. At the same time, the GPE simulations of the experiment in Fig. 5(a) and the experimental data in Fig. 3(b) show that the diameter of the central peak at this point is only m, which is still substantially smaller than the original torus. Because of the correspondingly large aspect ratio and the structure of concentric fringes seen in Fig. 1, the six images in Fig. 3(b) can be regarded as the first experimental demonstration of cold atom matter wave Bessel beams, covering all quantized orbital angular momentum states between and .

The GPE simulation in Fig. 5(a) also shows that the evolution is not completely diffraction-free. The small divergence is a result of the quadratic phase term in Eq. 1. Fig. 5(b) shows the diffraction-free evolution that would be obtained if the condensate had the negative of this phase imprinted on it to remove the curvature (equivalent in optics to collimation by a thin lens). Because the wavefront in the simulation does not have the infinite Bessel form it propagates diffraction-free for only a finite distance before blowing up, as is well-known from optical implementations Durnin1987 (). In principle such phase imprinting is easily accomplished experimentally by subjecting the BEC to a short pulse of far-detuned light Denschlag2000 (); Burger1999 () with a quadratic intensity profile. However, it was not possible to implement this procedure in our experiment because the poor optical quality of the cell windows prevented delivery of a beam with sufficient intensity smoothness to remove the phase variation of radians developed across the condensate after 20 ms of expansion. It should be easily demonstrated in future work with an improved system. Fig. 5(c) shows the evolution of just the central part of the beam after imprinting combined with masking to suppress the outer fringes. It propagates almost Heisenberg limited, with calculated beam quality factor Impens2008 () . It is interesting to note that the presence of the relatively faint outer fringes in the non-imprinted Bessel beam realized in our experiment [Fig. 5(a)] reduces the divergence of the central beam by a factor compared to the masked beam, increasing the intensity in the central peak (which contains of the atoms) by a factor . So the non-imprinted Bessel beam is, in a sense, better than Heisenberg-limited, and such beams may therefore be preferable to Heisenberg-limited atom lasers Jeppesen2008 () in applications where the additional fringes can be tolerated.

Our experiment can be regarded as a simple form of matter wave holography in which the initial atomic wavefunction is engineered so that the subsequent free evolution (possibly in conjunction with phase imprinting acting as a lens) transforms it into a desired form at some remote location. Since the Painted Potential can create arbitrary and dynamic potentials there is considerable scope for developing this idea further: the current implementation has about addressable spots, so it should already be possible to create a gray scale hologram with this resolution. While atom holography has previously been demonstrated with atomic beams passing through a mask serving as a binary transmission hologram Fujita1996 (); Morinaga1996 (), this form of holography is inefficient because of absorption by the mask, and the hologram cannot be changed once fabricated and installed inside the vacuum system. Holography with the Painted Potential would consist of forming a suitably shaped condensate by evaporation into a potential landscape corresponding to the desired amplitude distribution in the hologram, or of phase-imprinting a uniform BEC created in a flat box potential. Painted Potential holography would not waste any atoms through absorption at a mask and it would allow for dynamic and arbitrary holograms. It could, for example, lead to a new form of microlithography in which surfaces are patterned with atoms in arrangements that can be changed in real time.

This work was supported by the U.S. Department of Energy through the LANL/LDRD Program. We acknowledge inspiring conversations with Eddy Timmermans.

## References

- (1) J. Durnin, J. J. Miceli, and J. H. Eberly, Phys. Rev. Lett. 58, 1499 (1987).
- (2) K. Volke-Sepulveda, V. GarcÃ©s-ChÃ¡vez, S. ChÃ¡vez-Cerda, J. Arlt, and K. Dholakia, J. Opt. B-Quantum S. O. 4, S82 (2002).
- (3) M. Fortin, M. PichÃ©, and E. Borra, Opt. Express 12, 5887 (2004).
- (4) M. Erdelyi, Z. L. Horvath, G. Szabo, Z. Bor, F. K. Tittel, J. R. Cavallaro, and M. C. Smayling, J. Vac. Sci. Technol. B 15, 287 (1997).
- (5) G. Ok, S. W. Choi, K. H. Park, and H. S. Chun, Sensors 13, 71 (2013).
- (6) J. Arlt, V. Garces-Chavez, W. Sibbett, and K. Dholakia, Opt. Commun. 197, 239 (2001).
- (7) J. Chen, J. Ng, Z. Lin, and C. T. Chan, Nat. Photonics 5, 531 (2011).
- (8) D. McGloin, and K. Dholakia, Contemp. Phys. 46, 15 (2005).
- (9) J. y. Lu, and J. F. Greenleaf, IEEE T. Ultrason. Ferr. 37, 438 (1990).
- (10) M. Uchida, and A. Tonomura, Nature 464, 737 (2010).
- (11) J. Verbeeck, H. Tian, and P. Schattschneider, Nature 467, 301 (2010).
- (12) B. J. McMorran, A. Agrawal, I. M. Anderson, A. A. Herzing, H. J. Lezec, J. J. McClelland, and J. Unguris, Science 331, 192 (2011).
- (13) K. Henderson, C. Ryu, C. MacCormick, and M. G. Boshier, New J. Phys. 11, 043030 (2009).
- (14) C. Ryu, M. F. Andersen, P. Clade, V. Natarajan, K. Helmerson, and W. D. Phillips, Phys. Rev. Lett. 99, 260401 (2007).
- (15) J. Brand, and W. P. Reinhardt, J. Phys. B 34, L113 (2001).
- (16) M. G. Boshier, and W. J. Sandle, Opt. Commun. 42, 371 (1982).
- (17) C. Ryu, K. C. Henderson, and M. G. Boshier, Bull. Am. Phys. Soc., BAPS.2010.DAMOP.M1.31 (2010).
- (18) S. Moulder, S. Beattie, R. P. Smith, N. Tammuz, and Z. Hadzibabic, Phys. Rev. A 86, 013629 (2012).
- (19) K. C. Wright, R. B. Blakestad, C. J. Lobb, W. D. Phillips, and G. K. Campbell, Phys. Rev. Lett. 110, 025302 (2013).
- (20) G. B. Hess, and W. M. Fairbank, Phys. Rev. Lett. 19, 216 (1967).
- (21) J. Denschlag et al., Science 287, 97 (2000).
- (22) S. Burger, K. Bongs, S. Dettmer, W. Ertmer, K. Sengstock, A. Sanpera, G. V. Shlyapnikov, and M. Lewenstein, Phys. Rev. Lett. 83, 5198 (1999).
- (23) F. Impens, Phys. Rev. A 77, 013619 (2008).
- (24) M. Jeppesen, J. DuguÃ©, G. R. Dennis, M. T. Johnsson, C. Figl, N. P. Robins, and J. D. Close, Phys. Rev. A 77, 063618 (2008).
- (25) J. Fujita, M. Morinaga, T. Kishimoto, M. Yasuda, S. Matsui, and F. Shimizu, Nature 380, 691 (1996).
- (26) M. Morinaga, M. Yasuda, T. Kishimoto, F. Shimizu, J.-i. Fujita, and S. Matsui, Phys. Rev. Lett. 77, 802 (1996).