Abstract
We have formulated and experimentally demonstrated an improved algorithm for design of arbitrary two-dimensional holographic traps for ultracold atoms. Our method builds on the best previously available algorithm, MRAF and improves on it in two ways. First, it allows for creation of holographic atom traps with a well defined background potential. Second, we experimentally show that for creating trapping potentials free of fringing artifacts it is important to go beyond the Fourier approximation in modelling light propagation. To this end, we incorporate full Helmholtz propagation into our calculations.
Similar content being viewed by others
Introduction
Optical dipole traps have become ubiquitous in disciplines ranging from live cell manipulation in biophysics1 to single atom manipulation for quantum information processing2. In the field of ultracold atomic gases a variety of optical potential shapes, of ever increasing complexity, are used for fundamental studies of many-body physics3 in different optical crystal lattices4,5, reduced dimensionality6,7 and non-trivial trapping topologies8,9. Optical sculpting is also likely to offer great benefits for trapped-atom interferometry with Bose-Einstein condensates (BECs)10 and for creating âatomtronicâ11 optical circuits.
The simplest optical trap for an ultracold atomic gas is formed by a single focused Gaussian laser beam12. Natural extensions include standing-wave optical lattice potentials produced by the interference of multiple laser beams4 and ring traps produced using Laguerre-Gauss laser modes8,9. However there is so far no universal approach for creating an arbitrary optical trapping potential on the required micrometer scale. Existing techniques fall into three categories: (1) a time-averaged potential can be created by the fast scanning of a focused laser beam13, (2) an intensity mask or micromirror-device can be imaged onto the trapping plane of the atoms14,15 and (3) in the holographic method one manipulates the phase of the laser beam so as to create the desired intensity pattern after further propagation of the light. None of these methods is currently a clear overall winner; their relative merits depend in practice on various factors, such as the desired spatial resolution, efficiency of use of laser light and temporal stability of the light pattern.
In this paper, we provide a computational and experimental procedure for generating arbitrary potentials based on an improved holographic, phase only method. We illustrate our technique by producing high fidelity patterns which demonstrate suitable characteristics for atom trapping, opening a versatile toolbox for optical manipulation in exotic geometries.
As shown in Fig 1(a), we use a spatial light modulator (SLM) to imprint a custom phase pattern (kinoform) on a laser beam in order to form a diffraction pattern of the desired shape in the trapping plane. Since this is a diffractive method, the kinoform steers the light into required trap shapes without blocking out any of the incident power, making the method potentially very efficient. In some cases, the required kinoform can be calculated analytically8,9, but in general, this is a highly non-trivial numerical problem. The difficulty arises from the fact that we have incomplete a priori knowledge of the light field in both planes: In the trapping plane only the modulus of the light field is relevant for optical trapping, so we place no constraint on the phase; similarly, at the SLM, we impose the unaltered Gaussian modulus of the incident laser profile, leaving the phase as the variable to be calculated [Fig. 1(b)].
The numerical task of calculating the kinoform is approached as an optimisation procedure which can be achieved by âsteepest descentâ methods such as direct binary searches16,17,18, genetic algorithms19 or Gerchberg and Saxton's iterative Fourier transform methods20. Our method is based on the Gerchberg-Saxton algorithm (see the Methods section) since it is very computationally efficient and is easily executed on a graphics processing unit. So far, the best computational results were obtained by Pasienski and DeMarco21, using an adaptation of the Gerchberg-Saxton scheme known as the Mixed-Region-Amplitude-Freedom (MRAF) algorithm. The idea of MRAF is to enhance convergence of the iterative algorithm in one region of the trapping plane by giving up control of the remaining regions21,22,23,24. This approach produced excellent computational results for a range of generic trap shapes and some shapes have also been successfully realised with laser light25,26. However, one practical limitation of MRAF is that in order to achieve good convergence in bright areas of the trapping plane (green areas in Fig. 1), we must surrender control over almost all of the dark, background regions. In particular, the original MRAF calculations provide for only a very narrow region defining the zero background potential to which all the potentials are referenced. This could lead to experimental problems such as inefficient loading of atoms into, or percolation of atoms out of the trap; this issue is particularly important if the desired trap has sharp edges so the region of uncontrolled light intensity is immediately adjacent to the region of high atomic density. If the region of zero potential is extended, the fidelity of MRAF kinoforms is severely compromised. The first result of this paper is to show that this problem can be eliminated by properly redefining the target potential so as to include both the conventionally defined trapping region and a sufficiently large well defined background region (a âcanvasâ on which an arbitrary potential landscape is âdrawnâ). We formulate an offset-MRAF (OMRAF) algorithm through a simple set of practical rules which obey the limits imposed by Nyquist's theorem and minimise the occurrence of convergence-stalling optical vortices.
The second result of this paper is to demonstrate that the paraxial approximation often employed in literature21,27 can lead to undesirable artifacts in the traps. We illustrate the improvement that can be made by direct numerical calculation of the light propagation using the Helmholtz equation. Using our computational algorithms and a simple apparatus, we experimentally demonstrate examples of optical traps suitable for atomic experiments and compare the results to previous realisations25,26.
The paper is organised as follows: First we compare computational results produced using the existing MRAF and our OMRAF algorithms. Then we present and analyse experimental images of light fields shaped by our method. Throughout this investigation, we primarily test our technique by creating an atomtronic OR-gate trap shape11,21,28 (see Fig 1). To illustrate the universality of our method we also include the final results for a uniform square trap and an annular BEC stirrer21,25.
Results
The Gerchberg-Saxton algorithm and MRAF
Within the Gerchberg-Saxton framework, the core algorithm relies on linking the moduli in the SLM and trapping planes by simulating the light propagation back and fourth between these planes. After each propagation, we manually impose the known modulus constraints while leaving the phase to converge on the required solution. The mathematical details are presented in the Methods section and here we just briefly review the algorithms for kinoform calculation before presenting the numerical and experimental results. It is well documented that convergence of the naïve Gerchberg-Saxton algorithm is highly erratic because the modulus constraints in each plane are non-convex29,30. The MRAF algorithm improves on this by defining a âdrawingâ in the trapping plane containing all the points in which the desired potential is non-zero and an additional âcanvasâ region, , of zero potential around the drawing [Fig. 2(a)]. The modulus constraint is then manually imposed only inside and at each iterative step, thus reducing the number of constraints the algorithm aims to satisfy.
In the original MRAF paper, region is limited to a tight border around which has a typical width of around 1% of the trapping region length scales. Although is a featureless region of zero potential, it is an integral part of the trap since it defines the zero to which all the pixels in are referenced. For practical purposes such a small canvas can be a severe limitation. Our method aims to expand to approach its theoretical limit which is set by the Nyquist theorem to be 25% of the area of the simulated trapping plane (see the Methods section). In this large canvas regime, MRAF convergence stagnates and produces poor solutions even after many iterations (see Fig. 2)
Optical vortices and OMRAF
We observe the problems of the MRAF algorithm for extended to be due to the formation of large populations of optical vortices during the early iterations (approximately 1 vortex for every 10 pixels in Fig. 2(a)). Optical vortices are topological features of the light field corresponding to a phase winding around a point of zero intensity. Large vortex populations are problematic because individual vortices are difficult to eliminate in further iterations of the algorithm. The reason for this is that the Gerchberg-Saxton algorithm monotonically reduces the RMS error in each iteration31, but topological (un)winding operations cause disruption across the entire kinoform. Such global disruption is not in keeping with the monotonic error reduction in later iterations32. Therefore any erroneous vortices are âfrozen inâ early on and the algorithm gets stuck in a local RMS minimum corresponding to a particular vortex distribution.
One possibility for handling vortices is by pairwise creation and annihilation, which requires only local operations; computational algorithms which encourage such pairwise operations are given in Ref.[33]. However, we consider a much simpler solution. Our method simply offsets the trapping pattern inside and by a uniform intensity, |Î|2, to remove all points of zero light intensity. This redefines the zero of the potential but does not change the physics of the trap. With an offset intensity, we know a priori that we should aim to create no vorticies at all inside , so we feed the algorithm with an initial kinoform which contains no phase winding. Specifically, we choose a parabolic initial kinoform21,24 which defocuses the light into a patch with a characteristic size set by .
The optimal value of the offset intensity is calculated semi-empirically as follows: After just one iteration, the intensity in the trapping plane adopts approximately the correct shape, but exhibits considerable fluctuations around the desired intensity [Fig. 2(b)]. If any of these fluctuations bring the intensity locally to zero, a vortex may form at this point. Therefore, we increase |Î|2 until the trapping plane contains no vortices inside after the first iteration [Fig. 2(b) and (c)] and then trust that very few vortices will be formed in subsequent iterations. For maximum light-usage efficiency, we want to minimise the proportion of the incident light which is steered into the featureless background, i.e. choose the minimal |Î| which gives satisfactory results. We find that for all the patterns we tested a suitable offset was |Î| ~ 10 â 15% of the maximum amplitude in the trapping plane (i.e. |Î|2 ~ 1% of the maximum intensity). The computational simulation of the trapping potential produced using the OMRAF algorithm is shown in Fig. 2. Our method reduces the number of vortices seen in after 30 iterations by a factor of ~ 20 compared to MRAF. This allows excellent convergence of the algorithm, leaving only 4% RMS error after 30 iterations and 1% RMS after 1000.
One compromise associated with the OMRAF method is a reduction in the efficiency of light use. We define efficiency as the ratio of the integral of the trap intensity above the background offset to the total integral across the trapping plane. Including the non-zero Î causes a drop in efficiency by a factor of approximately 2, from 43% to 24%, for the OR-gate. Nevertheless, our method remains more efficient than intensity masking methods such as Ref.[14], which report 3% efficiency for simple patterns. We calculate that for more complex patterns such as the OR-gate, intensity masking would be less than 1% efficient.
Experimental realisation with laser light
We test our OMRAF algorithm experimentally by producing light patterns using 532â nm laser light in the arrangement shown in Fig. 1(a). We use the standard SLM-based Shack-Hartmann algorithm34 to crudely correct for the low spatial frequency phase aberrations in the optical system and characterise our input laser beam. We then use an active feedback algorithm25 to optimise the final experimental patterns. This feedback routine is essential to remove the âsinc envelopeâ caused by the pixellation of the SLM, which would otherwise globally modulate the intensity pattern. Finally, in practice, we find that when we try to produce experimental patterns with covering the maximum theoretically permitted area (25% of the trapping plane), we cannot achieve patterns with less than 20% RMS fluctuations. This is significantly improved in the results presented below by reducing to cover only 10% of the trapping plane.
In the previous section, we presented computational results generated under the paraxial approximation in which the propagator mapping the light in the SLM plane to the trapping plane is given by a scaled Fourier transform27. Using this approach in our experiments with an f = 200â mm lens, we are able to produce a trapping pattern with RMS error of â 11% [Fig. 3(a)]. As highlighted in Fig. 3(b), the most prominent form of error is semi-regular fringing. We suggest that the source of this error is the inadequacy of the paraxial approximation and use a numerical Helmholtz solver (see the Methods section) to test this hypothesis. Specifically, we use the OMRAF algorithm under the paraxial approximation to produce a âFourier-generatedâ kinoform and then input this kinoform into our numerical Helmholtz propagator to simulate the experimental apparatus. This indeed produces the fringing pattern similar to the one observed experimentally, which gives us confidence to continue with the Helmholtz method. We thus replace all instances of the Fourier propagator in the OMRAF algorithm with a Helmholtz propagator.
In this way, we are able to eliminate the erroneous fringing and reduce the RMS error of the experimental pattern to 7% [see Fig. 3(a)]. As shown in Fig. 3(c), we can produce other trap shapes with similar fidelity. Thus, for traps with sharp edges and an extended canvas region, we achieve RMS variation comparable to the 4% fluctuations previously seen only for simple smooth traps25. In Ref.[35], it has already been shown that 5% RMS errors are sufficiently low for atomtronic applications. We therefore believe that our methods show great promise for future applications.
Discussion
We have presented a simple method for producing optical traps of arbitrary 2D profile suitable for use in cold atom experiments. A useful optical potential must be referenced to a well defined background and we showed that in order to do this, we must modify the existing algorithms to include a non-zero background light intensity. Specifically, we developed the OMRAF algorithm guided by the principle that points of zero intensity seed convergence-stalling vortices and therefore should be avoided. In computer simulations we produce traps covering 25% of the trapping plane (the Nyquist limit) with 4% RMS error after only 30 iterations. We also took the next step of realising these 2D profiles in laser light. Here, we remove systematic aberrations using a Shack-Hartmann method and active feedback. In addition, we found that using a beyond-paraxial Helmholtz solver improves the fidelity of the experimentally produced light patterns.
Methods
Modified Gerchberg-Saxton algorithms
Both the MRAF and the OMRAF approach have the Gerchberg-Saxton algorithm at their core. This algorithm attempts to obtain the kinoform, Ï(r), which must be applied to the input Gaussian field, E0(r), in the SLM plane (spanned by r) to obtain the desired trap shape, T(R), in the trapping plane (spanned by R) at z = 2f [see Fig. 1(a)]. These two planes are related by a projection operator, , so we may write the implicit equation for Ï(r) as:
The Gerchberg-Saxton algorithm solves this equation iteratively as illustrated in Fig. 4.
MRAF relaxes the modulus constraint outside the trapping region by defining a drawing, (defined by T(R) > 0) together with a narrow canvas, , of zero intensity around and then applying the following algorithm:
We find that a mixing parameter m = 0.4 gives the minimum RMS deviation for all the patterns trialled (see also Ref.[21]).
OMRAF further modifies this algorithm in the following ways: First we expand the canvas so that the trapping region covers the maximum theoretical area set by Nyquist's theorem (25% of the whole trapping plane - see the next section). Secondly, throughout and , we offset the desired light intensity according to:
This shifts both the trapping potential and its zero reference level equally and does not change the physics of the trap.
Maximum canvas size
In order to best use the degrees of freedom offered by an N Ã N pixel SLM, we should allow region to cover as many pixels as possible in the discretized trapping plane. This theoretical limit is found by the following Nyquist argument: The N Ã N SLM plane can be fast-Fourier-transformed to an N Ã N trapping plane with a diffraction limited spot size of exactly 1 pixel. The Nyquist limit states that the diffraction limit should be reduced to half a pixel (the smallest feature in the trap shape). This can be done by padding the computer-simulated SLM plane with zeros to increase its size to 2N Ã 2N. This means that only 25% of the simulated SLM plane is actually covered by the physical SLM, so we can at best control only 25% of the trapping plane.
A computationally efficient helmholtz solver
We see improvement in the experimentally realised trapping intensities when the following method is used to model the propagation operator in the OMRAF iterations:
The Helmholtz equation for light propagation in free space gives that the Fourier transform of a light field, E(r; z), is modified by a phase factor Î(ξ; z) as it propagates in the +z direction27
where denotes the Fourier transform (FT) of g, making the Fourier variable pairing r ξ explicit and k = 2Ï/λ denotes the wave vector of the light. The routine for propagation of E0 via the apparatus in Fig. 1(a) can then be conceptually decomposed into 3 stages:
where represents convolution and represents the FT of the phase pattern imprinted by the lens, which, after aberration correction, we assume to be:
References
Ashkin, A. & Dziedzic, J. M. Optical trapping and manipulation of viruses and bacteria. Science 235, 1517â1520 (1987).
Beugnon, J. et al. Two-dimensional transport and transfer of a single atomic qubit in optical tweezers. Nature Phys. 3, 696â699 (2007).
Bloch, I., Dalibard, J. & Zwerger, W. Many-body physics with ultracold gases. Rev. Mod. Phys. 80, 885â964 (2008).
Greiner, M., Mandel, O., Esslinger, T., Hänsch, T. W. & Bloch, I. Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms. Nature 415, 39â44 (2002).
Struck, J. et al. Quantum simulation of frustrated classical magnetism in triangular optical lattices. Science 333, 996â999 (2011).
Kinoshita, T., Wenger, T. & Weiss, D. S. Observation of a one-dimensional Tonks-Girardeau gas. Science 305, 1125â1128 (2004).
Hadzibabic, Z., Krüger, P., Cheneau, M., Battelier, B. & Dalibard, J. Berezinskii-Kosterlitz-Thouless crossover in a trapped atomic gas. Nature 441, 1118â1121 (2006).
Ramanathan, A. et al. Superflow in a toroidal Bose-Einstein condensate: an atom circuit with a tunable weak link. Phys. Rev. Lett. 106, 130401 (2011).
Moulder, S., Beattie, S., Smith, R. P., Tammuz, N. & Hadzibabic, Z. Quantized supercurrent decay in an annular Bose-Einstein condensate. Phys. Rev. A 86, 013629 (2012).
Shin, Y. et al. Atom interferometry with Bose-Einstein condensates in a double-well potential. Phys. Rev. Lett. 92, 050405 (2004).
Seaman, B. T., Krämer, M., Anderson, D. Z. & Holland, M. J. Atomtronics: ultracold-atom analogs of electronic devices. Phys. Rev. A 75, 023615 (2007).
Stamper-Kurn, D. M., Andrews, M. R., Chikkatur, A. P., Inouye, S., Miesner, H.-J. Stenger, J. & Ketterle, W. Optical confinement of a BoseâEinstein condensate. Phys. Rev. Lett. 80, 2027â2030 (1998).
Henderson, K., Ryu, C., MacCormick, C. & Boshier, M. G. Experimental demonstration of painting arbitrary and dynamic potentials for Bose-Einstein condensates. New J. Phys. 11, 043030 (2009).
Scherer, D. R., Weiler, C. N., Neely, T. W. & Anderson, B. P. Vortex formation by merging of multiple trapped Bose-Einstein condensates. Phys. Rev. Lett. 98, 110402 (2007).
Liang, J., Kohn, R. N., Becker, M. F. & Heinzen, D. 1.5% root-mean-square flat-intensity laser beam formed using a binary-amplitude spatial light modulator. Appl. Optics 48, 1955â1962 (2009).
Seldowitz, M. A., Allebach, J. P. & Sweeney, D. W. Synthesis of digital holograms by direct binary search. Appl. Optics 26, 2788â2798 (1987).
Boyer, V., Chandrashekar, C. M. & Foot, C. J. Dynamic optical trap generation using FLC SLMs for the manipulation of cold atoms. J. Mod. Optic. 51, 2235â2240 (2004).
Boyer, V. et al. Dynamic manipulation of Bose-Einstein condensates with a spatial light modulator. Phys. Rev. A 73, 031402 (2006).
Zhou, G., Chen, Y., Wang, Z. & Song, H. Genetic local search algorithm for optimization design of diffractive optical elements. Appl. Opt. 38, 4281â4290 (1999).
Gerchberg, R. W. & Saxton, W. O. A practical algorithm for the determination of phase from image and diffraction plane pictures. Optik 35, 237â246 (1972).
Pasienski, M. & DeMarco, B. A high-accuracy algorithm for designing arbitrary holographic atom traps. Opt. Express 16, 2176â2190 (2008).
Wyrowski, F. Diffractive optical elements: iterative calculation of quantized, blazed phase structures. J. Opt. Soc. Am. A 7, 961â969 (1990).
Akahori, H. Spectrum leveling by an iterative algorithm with a dummy area for synthesizing the kinoform. Appl. Optics 25, 802â811 (1986).
Aagedal, H., Schmid, M., Beth, T., Teiwes, S. & Wyrowski, F. Theory of speckles in diffractive optics and its application to beam shaping. J. Mod. Optic. 43, 1409â1421 (1996).
Bruce, G. D., Mayoh, J., Smirne, G., Torralbo-Campo, L. & Cassettari, D. A smooth, holographically generated ring trap for the investigation of superfluidity in ultracold atoms. Phys. Scr. T143, 014008 (2011).
Bruce, G. D., Bromley, S. L., Smirne, G., Torralbo-Campo, L. & Cassettari, D. Holographic power-law traps for the efficient production of Bose-Einstein condensates. Phys. Rev. A 84, 053410 (2011).
Ersoy, O. K. Diffraction, Fourier Optics and Imaging. Wiley Series in Pure and Applied Optics (Wiley, 2006).
Pepino, R. A., Cooper, J., Meiser, D., Anderson, D. Z. & Holland, M. J. Open quantum systems approach to atomtronics. Phys. Rev. A 82, 013640 (2010).
Levi, A. & Stark, H. Image restoration by the method of generalized projections with application to restoration from magnitude. J. Opt. Soc. Am. A 1, 932â943 (1984).
Bauschke, H. H., Combettes, P. L. & Luke, D. R. Phase retrieval, error reduction algorithm and Fienup variants: a view from convex optimization. J. Opt. Soc. Am. A 19, 1334â1345 (2002).
Fienup, J. R. Phase retrieval algorithms: a comparison. Appl. Opt. 21, 2758â2769 (1982).
Wyrowski, F. & Bryngdahl, O. Iterative Fourier-transform algorithm applied to computer holography. J. Opt. Soc. Am. A 5, 1058â1065 (1988).
Senthilkumaran, P., Wyrowski, F. & Schimmel, H. Vortex stagnation problem in iterative Fourier transform algorithms. Opt. Laser Eng. 43, 43â56 (2005).
Bowman, R. W., Wright, A. J. & Padgett, M. J. An SLM-based Shack-Hartmann wavefront sensor for aberration correction in optical tweezers. J. Opt. 12, 124004 (2010).
Wright, K. C., Blakestad, R. B., Lobb, C. J., Phillips, W. D. & Campbell, G. K. Driving phase slips in a superfluid atom circuit with a rotating weak link. eprint arXiv: 1208.3608 (2012).
Acknowledgements
We thank Naaman Tammuz, Stuart Moulder, Scott Beattie and Richard Bowman for helpful discussions and Robert Smith for comments on the manuscript. This work was supported by EPSRC (Grant No. EP/G026823/1).
Author information
Authors and Affiliations
Contributions
AG performed the simulations and experiments. Both authors contributed equally to the conception of the project, discussion of the results and writing of the manuscript.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
This work is licensed under a Creative Commons Attribution 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by/3.0/
About this article
Cite this article
Gaunt, A., Hadzibabic, Z. Robust Digital Holography For Ultracold Atom Trapping. Sci Rep 2, 721 (2012). https://doi.org/10.1038/srep00721
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep00721