OSA's Digital Library

Optics Express

Optics Express

  • Editor: Andrew M. Weiner
  • Vol. 21, Iss. 8 — Apr. 22, 2013
  • pp: 10188–10198
« Show journal navigation

Sensitive absorption imaging of single atoms in front of a mirror

Atreju Tauschinsky and Robert J. C. Spreeuw  »View Author Affiliations


Optics Express, Vol. 21, Issue 8, pp. 10188-10198 (2013)
http://dx.doi.org/10.1364/OE.21.010188


View Full Text Article

Acrobat PDF (1488 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

In this paper we show that the sensitivity of absorption imaging of ultracold atoms can be significantly improved by imaging in a standing-wave configuration. We present simulations of single-atom absorption imaging both for a travelling-wave and a standing-wave imaging setup, based on a scattering approach to calculate the optical density of a single atom. We find that the optical density of a single atom is determined only by the numerical aperture of the imaging system. We determine optimum imaging parameters, taking all relevant sources of noise into account. For reflective imaging we find an improvement of 1.7 in the maximum signal-to-noise ratio can be achieved. This is particularly useful for imaging in the vicinity of an atom chip, where a reflective surface is naturally present.

© 2013 OSA

1. Introduction

The detection of single atoms with high fidelity is an important requirement for many modern atomic physics experiments, and in particular for quantum information science. Single particle imaging has been achieved using fluorescence imaging of ions [1

1. W. Neuhauser, M. Hohenstatt, P. Toschek, and H. Dehmelt, “Localized visible Bâ{+} mono-ion oscillator,” Phys. Rev. A. 22, 1137–1140 (1980) [CrossRef] .

] and neutral atoms [2

2. N. Schlosser, G. Reymond, I. Protsenko, and P. Grangier, “Sub-poissonian loading of single atoms in a microscopic dipole trap.” Nature 411, 1024–7 (2001) [CrossRef] [PubMed] .

, 3

3. K. D. Nelson, X. Li, and D. S. Weiss, “Imaging single atoms in a three-dimensional array,” Nat. Phys. 3, 556–560 (2007) [CrossRef] .

, 4

4. M. Wilzbach, D. Heine, S. Groth, X. Liu, T. Raub, B. Hessmo, and J. Schmiedmayer, “Simple integrated single-atom detector,” Opt. Lett. 34, 259–261 (2009) [CrossRef] [PubMed] .

, 5

5. R. Bücker, A. Perrin, S. Manz, T. Betz, C. Koller, T. Plisson, J. Rottmann, T. Schumm, and J. Schmiedmayer, “Single-particle-sensitive imaging of freely propagating ultracold atoms,” New J. Phys. 11, 103039 (2009) [CrossRef] .

, 6

6. D. Heine, W. Rohringer, D. Fischer, M. Wilzbach, T. Raub, S. Loziczky, X. Liu, S. Groth, B. Hessmo, and J. Schmiedmayer, “A single-atom detector integrated on an atom chip: fabrication, characterization and application,” New J. Phys. 12, 095005 (2010) [CrossRef] .

], but is difficult in atom-chip experiments due to a scattering background from the chip [4

4. M. Wilzbach, D. Heine, S. Groth, X. Liu, T. Raub, B. Hessmo, and J. Schmiedmayer, “Simple integrated single-atom detector,” Opt. Lett. 34, 259–261 (2009) [CrossRef] [PubMed] .

]. A different approach is to enhance the atom-photon interaction using resonant cavities [7

7. J. Goldwin, M. Trupke, J. Kenner, a. Ratnapala, and E. a. Hinds, “Fast cavity-enhanced atom detection with low noise and high fidelity.” Nature Commun. 2, 418 (2011) [CrossRef] .

, 8

8. J. Bochmann, M. Mücke, C. Guhl, S. Ritter, G. Rempe, and D. L. Moehring, “Lossless State Detection of Single Neutral Atoms,” Phys. Rev. Lett. 104, 203601 (2010) [CrossRef] [PubMed] .

, 9

9. R. Gehr, J. Volz, G. Dubois, T. Steinmetz, Y. Colombe, B. L. Lev, R. Long, J. Estève, and J. Reichel, “Cavity-Based Single Atom Preparation and High-Fidelity Hyperfine State Readout,” Phys. Rev. Lett. 104, 203602 (2010) [CrossRef] [PubMed] .

]. However, free-space absorption imaging of single particles [10

10. D. J. Wineland, W. M. Itano, and J. C. Bergquist, “Absorption spectroscopy at the limit: detection of a single atom.” Opt. Lett. 12, 389–391 (1987) [CrossRef] [PubMed] .

, 11

11. E. W. Streed, A. Jechow, B. G. Norton, and D. Kielpinski, “Absorption imaging of a single atom.” Nature Commun. 3, 933 (2012) [CrossRef] .

, 12

12. M. K. Tey, G. Maslennikov, T. C H Liew, S. A. Aljunid, F. Huber, B. Chng, Z. Chen, V. Scarani, and C. Kurtsiefer, “Interfacing light and single atoms with a lens,” New J. Phys. 11, 043011 (2009) [CrossRef] .

] is still a daunting task.

In this paper we propose absorption imaging in a standing wave as a simple and effective way to improve the signal-to-noise ratio (SNR) of atom detection. This is particularly useful in the context of atom-chip experiments [13

13. D. A. Smith, S. Aigner, S. Hofferberth, M. Gring, M. Andersson, S. Wildermuth, P. Krüger, S. Schneider, T. Schumm, and J. Schmiedmayer, “Absorption imaging of ultracold atoms on atom chips,” Opt. Express 19, 8471 (2011) [CrossRef] [PubMed] .

] where a standing wave naturally forms if the imaging is done perpendicular to the chip surface. We will take the setup currently in use in our group as a reference [14

14. V. Y. F. Leung, A. Tauschinsky, N. J. Druten, and R. J. C. Spreeuw, “Microtrap arrays on magnetic film atom chips for quantum information science,” Quantum Inf. Proc. 10, 955–974 (2011) [CrossRef] .

, 15

15. R. Gerritsma, S. Whitlock, T. Fernholz, H. Schlatter, J. A. Luigjes, J.-U. Thiele, J. B. Goedkoop, and R. J. C. Spreeuw, “Lattice of microtraps for ultracold atoms based on patterned magnetic films,” Phys. Rev. A. 76, 033408 (2007) [CrossRef] .

, 16

16. S. Whitlock, R. Gerritsma, T. Fernholz, and R. J. Spreeuw, “Two-dimensional array of microtraps with atomic shift register on a chip,” New J. Phys. 11, 023021 (2009) [CrossRef] .

] for realistic experimental parameters.

In the first section we will calculate the absorption signal of a single atom, both for traditional one-pass imaging and for reflection imaging as employed in our experiment. In the second half of the paper we present numerical simulations of the expected absorption images and signal-to-noise ratio (SNR) under realistic conditions, taking into account all common sources of noise, for in-trap atoms as well as for untrapped atoms. To do so we use a hybrid approach, first calculating the movement of the atom in the imaging beam based on scattering individual photons, and subsequently investigating the absorption signal due to this atom, now taking the imaging beam as a classical wave. We show that using reflection imaging drastically increases the expected SNR and thus the single-atom detection accuracy.

2. Single-atom absorption imaging

Our starting point is the absorption signal of a single atom at rest, illuminated by a resonant laser beam with homogeneous intensity distribution (see Fig. 1). As is shown in the appendix the intensity in the imaging plane is then given by
Iabs(ρ)Iin(ρ)=σ1+s[Re(p(ρ))aχ|p(ρ)|2]
(1)
where Iabs(ρ) is the absorbed intensity, i.e. the missing intensity in the image plane with image-plane coordinates ρ, Iin is the incoming intensity, σ is the absorption cross section of the atom, s the saturation parameter in the object plane, p(ρ) the point-spread function of the imaging system (see below), a = [∫ |p(ρ)|2d2ρ]−1 the effective area over which the signal is distributed in the detection plane and χ is the fraction of the scattered light collected by the imaging system.

Fig. 1 Sketch of the imaging system considered for this paper. On the left the situation of a travelling wave without reflection is depicted. On the right the situation of reflective imaging (with the atom chip on top) is shown.

In an experiment the above quantities a and Iabs(0) are not directly accessible due to the finite size of any detector and the finite time resolution. Instead one would observe the quantities
Ndet(x)=Adρ0τdt(IinIabs)Nref(x)=Adρ0τdtIin
(2)
where x is an index over camera pixels, A is the size of a pixel and τ is the exposure time, i.e. the duration of the imaging pulse. Both σ and s now depend on time, as the atom acquires a finite velocity during the imaging pulse, leading to Doppler shifts. In addition the local saturation parameter in the object plane varies in the case of standing wave imaging, as will be discussed below. We can then define the apparent column density of atoms per pixel in the image plane as
n(x)=1+s0σ0Ndet(x)Nref(x)
(3)
where s0 is the saturation parameter for the intensity of the incoming beam and σ0 is the absorption cross section of an atom at rest for the incoming intensity. The apparent atom number per pixel would be given by Napp(x)= An(x).

As the atoms will move out of the focal plane during the imaging pulse, simulating this signal will require a point-spread function that is a function not only of the image plane coordinates, but also takes defocussing into account. This PSF can be numerically determined as
p(ρ,f;t)=2πρ0201rexp(ifr2)×J0(2πρrρ0)dr
(4)
for any given defocus f=2πλZ(11NA2) where Z is a real-space coordinate in the imaging direction [17

17. A. J. E. M. Janssen, “Extended Nijboer–Zernike approach for the computation of optical point-spread functions,” J. Opt. Soc. Am. A. 19, 849 (2002) [CrossRef] .

], NA is the numerical aperture of the imaging system and both ρ and Z might depend on time. ρ0 = λ/NA is proportional to the resolution of the imaging system, determined by the wavelength and the numerical aperture. Aberrations in the imaging system can be taken into account in a similar manner. We assume these to be negligible near the optical axis. In the simulations the numerical point-spread function given above is used to calculate the intensity in the image plane.

2.1. Analytical estimates

For a better intuition of the above results we consider some simple analytic estimates. Assuming a stationary atom in perfect focus of a two-lens imaging system with unit magnification (f −2ff) we can use a simple point-spread function of the form p(ρ)= 1/(ρρ0)J1(2πρ/ρ0) where p(0)=π/ρ02=a1. Using χ=(38)NA2 and assuming a 2-level system where σ0 = 3λ2/2π the maximum signal at the center of the image simplifies to
Iabs(0)Iin(0)=11+s0(32NA2916NA4)
(5)
which is thus determined purely by the numerical aperture of the imaging system and the saturation parameter s0. For a numerical aperture NA = 0.4 as is the case in our experiment we then find an absorption in the centre of the image of Iabs(0)/Iin(0)= 0.23 in the low saturation limit s0 ≪ 1 for a single atom.

We can further consider the total apparent number of atoms one would extract from this signal as the peak signal times the area over which the signal is distributed,
Napp=a1+s0σ0Iabs(0)Iin(0)=138NA2
(6)
This shows that the expected number of atoms is reduced by the second, NA-dependent term. This term is due to light which is scattered by the atom, but scattered into the solid angle of the lens, thus reducing the apparent amount of light that is absorbed. For a numerical aperture of 0.4 we therefore expect to find an apparent atom number Napp of 0.94 rather than one.

The commonly used quantity optical density (OD) can thus be defined for a single particle by the point-spread function of the imaging system, i.e. OD=ln(1IabsIin). At the center of the image (for vanishing pixel size) we find for the peak optical density a value of OD(0)32NA2+916NA4 in the low-saturation limit (where the sign change is due to a series expansion of the logarithm). For our numerical aperture this equals a peak optical density of 0.26. It is worth noting here that the atom number extracted in the usual way from optical density, in this case
NOD=a1+s0σ0ln(1Iabs(0)Iin(0))
(7)
would yield 1.07 atoms rather than 0.94. The atom number is slightly overestimated, as Lambert-Beers law is valid only for a continuously absorbing medium, not for a single absorber. In the following we always use the spatially dependent equivalent of (6) rather than (7) for the atom number.

2.2. Reflection imaging

The above discussion is for an atom imaged by a travelling wave, as depicted in Fig. 1(a). In reflection imaging, such as in Fig. 1(b) the situation is slightly different. In this case the probe light forms a standing wave at the position of the atoms. For Rb atoms in our magnetic lattice [14

14. V. Y. F. Leung, A. Tauschinsky, N. J. Druten, and R. J. C. Spreeuw, “Microtrap arrays on magnetic film atom chips for quantum information science,” Quantum Inf. Proc. 10, 955–974 (2011) [CrossRef] .

, 15

15. R. Gerritsma, S. Whitlock, T. Fernholz, H. Schlatter, J. A. Luigjes, J.-U. Thiele, J. B. Goedkoop, and R. J. C. Spreeuw, “Lattice of microtraps for ultracold atoms based on patterned magnetic films,” Phys. Rev. A. 76, 033408 (2007) [CrossRef] .

, 16

16. S. Whitlock, R. Gerritsma, T. Fernholz, and R. J. Spreeuw, “Two-dimensional array of microtraps with atomic shift register on a chip,” New J. Phys. 11, 023021 (2009) [CrossRef] .

] we expect a ground-state size of the atomic wavefunction of about 40nm, about a factor of 10 smaller than the period of the standing wave. We therefore approximate the atom as a point-particle which can initially be positioned in an anti-node of the standing wave. In our atom chip experiment the trap position can be adjusted by means of an externally controlled magnetic field to achieve this.

The intensity of the probe beam then varies sinusoidally along the probe direction and the maximum intensity in the anti-nodes is a factor four greater than the intensity of the incoming travelling wave. The coherent scattering amplitude 𝒜 scales with the the local field amplitude (at low saturation). Thus the scattered wave is two times as strong as for a travelling wave, leading to an accordingly higher observed signal. More precisely, the saturation parameter will also change.

It should be noted here that while in the limit of low saturation and short exposure times the observed signal will indeed be two times larger than for a travelling wave, for realistic imaging with finite exposure time the atom starts to probe the spatial varation of the light field, leading to an observed signal which depends on the exact imaging parameters. The expected signal can then only be predicted by simulations such as described below.

A further important difference to travelling wave imaging is that the atom is not pushed out of the focus of the imaging system by the probe beam, but performs a random walk in all dimensions. This point will also be further discussed in the simulations below.

Another factor of two could in principle be gained because the scattered wave is emitted towards and reflected by the mirror. This again doubles the amplitude of the scattered wave in the detection plane. A practical consideration here is that for single-atom imaging one will typically use high-NA imaging, with a Rayleigh length zR = λ / (πNA2) ≲ 10μm. For the typical situation in our experiment the atom is at about 10μm from the surface. This implies that the atom and its mirror image are not simultaneously in focus in the image plane. In practice this reflected wave hardly contributes to the optical density for our experimental parameters. If the imaging is done extremely close to a reflecting surface, or if the Rayleigh length of the imaging system is significantly larger this can however be a significant effect.

3. Simulations

In our simulations we distinguish four cases:
  1. free atom in travelling wave
  2. trapped atom in travelling wave
  3. free atom in standing wave
  4. trapped atom in standing wave
where the cases a) and b) are for traditional one-pass imaging, and the cases c) and d) are for our reflective atom-chip setup (see Fig. 1). While case a) and c) assume a freely moving atom, case b) and d) assume the atom to be held in a magnetic trap. For this trap we assume the potential calculated from the magnetic pattern of microtraps present on our atom chip [14

14. V. Y. F. Leung, A. Tauschinsky, N. J. Druten, and R. J. C. Spreeuw, “Microtrap arrays on magnetic film atom chips for quantum information science,” Quantum Inf. Proc. 10, 955–974 (2011) [CrossRef] .

] which is approximately equivalent to trap frequencies of (ωx, ωy, ωz)= 2π×(38.1, 36.5, 14.0)kHz and a trap depth of 8.8G. In the simulations the calculated potential is used without harmonic approximations. The atom is treated as free as soon as it acquires sufficient energy to escape the trap. We assume optical pumping to non-trapped states is negligible for an appropriate choice of polarisation for the imaging beam.

3.1. Recoil blurring

As the atom is illuminated by the imaging beam it scatters photons at a rate determined by the local intensity at the position of the atom as well as any detuning due to e.g. Doppler shifts as the atom acquires momentum with each of these scattering events. In the case of a travelling wave the photon scattering drives the atom strongly along the direction of the imaging beam, as the absorption of a photon always happens in this direction, performing a biased random walk. This quickly leads to the atom moving out of the focal plane of the imaging lens. The atom further performs an unbiased random walk in momentum space in the plane perpendicular to the imaging direction due to the spontaneous emission of a scattered photon, leading to a blurring of the image. In the case of a standing wave the atom can absorb a photon from both directions along the imaging axis, giving rise to an unbiased random walk in the imaging direction as well as the plane perpendicular to it. The atom therefore leaves the focus of the lens much more slowly than in the case of a travelling wave.

We simulate the movement of the atoms by calculating stochastic trajectories where random scattering events occur at the local scattering rate, taking the effect of the trap into account as a classical potential. In the case of an untrapped travelling wave we recover the expected behaviour, i.e. the root-mean-square (RMS) position along the imaging beam growing as t2 and the RMS position in a plane perpendicular to the imaging direction growing as t3/2, as is expected for uniform acceleration and a random walk in momentum space respectively. The effect of the trap is almost negligible in the case of a travelling wave, as the force of the light field is much greater than the restoring force of the trap potential. In the case of reflective imaging the behaviour proportional to t3/2 is also expected in the probe direction. This is however modified by the periodic intensity of the standing wave. Since the net force on the atoms is much lower than in the case of a travelling wave the influence of the trapping field is also more important, particularly at low intensities.

These calculated RMS trajectories are then used in the integrations of Eq. (2) to determine the signal after a given exposure time for a certain probe intensity.

3.2. Detection noise

To determine the feasibility of absorption imaging for single atoms one needs to consider various sources of noise in addition to the atomic trajectories to determine the signal-to-noise ratio. In the following we assume near shot-noise limited imaging, with further contributions from camera read-out noise and dark counts, as well as a finite quantum efficiency of the imaging system. Equations (2) determine the number of photons per pixel in the signal and light fields, and can be trivially modified to take a finite quantum-efficiency of the imaging system into account. With these values we can calculate the expected signal n(x) as atom density per pixel (Eq. 3) for any combination of exposure time and intensity. The variance of this signal is given by
σn(x)2=(1+s0σ0)2(σdet2Nref2+Ndet2σref2Nref4)
(8)
where Ndet and Nref are the number of electrons per pixel in the signal and reference images respectively, also integrated over pixel size and exposure time. σdet2 and σref2 are the variances of these quantities, mostly determined by photon shot noise (∝ N) and camera readout noise.

The signal to noise ratio per pixel is then given by n(x)σn(x); one can further improve the signal to noise ratio by not only evaluating the central pixel, but using an optimal estimator achieving the Cramér-Rao bound, such as
N=n(x)q(x)q(x)2
(9)
where q(x) is the spatial mode function of the signal [18

18. C. Ockeloen, A. Tauschinsky, R. Spreeuw, and S. Whitlock, “Detection of small atom numbers through image processing,” Phys. Rev. A. 82, 061606 (2010) [CrossRef] .

]. In our results we list the SNR both for evaluating only the central pixel and for an estimator achieving the Cramér-Rao bound.

We define the accuracy of the measurement as the probability of finding true positives p1,1, i.e. finding one atom if there is one atom present. Here we assume all measurement outcomes for which 0.5 < N < 1.5 to indicate the presence of exactly one atom. We have determined this value by simulations of 106 individual absorption images for each case. The results are in excellent agreement with those expected for a normal distribution, for which the accuracy is given by F=Erf(z2) with z=N2σ.

We can also determine the probability of finding false positives in the absence of any atoms p1,0, i.e. finding one atom if zero atoms are present. False positives due to the presence of more than one atom can not easily be determined, as an accurate calculating of the absorption signal of two atoms is non-trivial. We expect these to be of similar magnitude as those due to the absence of any atoms.

3.3. Simulation results

We simulate the expected absorption image and SNR assuming a numerical aperture of the objective lens of 0.4 (Edmund Optics NT47-727 with custom AR coating), a total magnification of the imaging system of 10, a pixel size of the camera of (13 × 13)μm2 (Andor iKon-M 934), a quantum efficiency of 0.9 and a readout noise of 13 counts per pixel [18

18. C. Ockeloen, A. Tauschinsky, R. Spreeuw, and S. Whitlock, “Detection of small atom numbers through image processing,” Phys. Rev. A. 82, 061606 (2010) [CrossRef] .

]. We consider resonant σ+-transitions on the D2-line of Rubidium 87, where σ0 = 2.907 × 10−9 cm2 and the saturation intensity Isat = 1.669mW/cm2[19

19. D. A. Steck, “Rubidium 87 D Line Data,” http://steck.us/alkalidata2 (2010).

]. We have confirmed that in the limit of low saturation our simulations yield the expected number of atoms in all four cases.

We find optimum imaging parameters by numerically maximizing the signal-to-noise ratio of the central pixel as a function of exposure time and intensity, i.e. maximizing SNR(0) where
SNR(x)=n(x)σn(x)
(10)
using the precalculated trajectories of an atom in the integrations of Eqs. (2). To avoid confusion in the case of reflective imaging where the intensity is spatially varying, the intensity is here always given in units of saturation intensity for the incoming beam (cf. Fig. 1). The local intensity at the initial position of the atoms is therefore higher by a factor of four for the standing wave. We also determine the SNR of an estimator achieving the Cramér-Rao bound by weighting the signal by the mode function of the atomic distribution as described above. SNR as a function of exposure time and intensity is depicted in Fig. 2. Table 1 lists the signal-to noise ratio for optimum imaging parameters. Note that this can be further improved by a factor 2 by using the fringe-removal algorithm[18

18. C. Ockeloen, A. Tauschinsky, R. Spreeuw, and S. Whitlock, “Detection of small atom numbers through image processing,” Phys. Rev. A. 82, 061606 (2010) [CrossRef] .

] which reduces the shot noise in the light field by optimal averaging of reference images even in the absence of any fringes. Fig. 3 finally shows simulated single-shot absorption images of a single atom for the parameters given above, including the use of fringe removal. The number of atoms extracted from these images using the optimum estimator is given in the figure caption.

Fig. 2 Simulated SNR as function of intensity of the incoming beam (in units of saturation intensity) and exposure time. The four subplots correspond to the four cases described in the text as indicated near the top and right axes, the color scale is normalized to the maximum SNR of case d). The black dot marks the position of the optimum.

Table 1. Optimum exposure parameters, and resulting SNR for the four cases described in the text. SNRpx is the SNR obtainable from evaluating a single pixel, while SNRCRB is the SNR for an estimator achieving the Cramér-Rao bound (without taking fringe-removal into account). The normalization factor Napp for case a) and b) is as expected from the analytic treatment of section 2.1. The normalization factor in cases c) and d) can only be obtained from the simulations. Finally we give the probability of true positive measurements p1,1 and false positive measurements due to the presence of zero atoms p1,0 (see 3.2).

table-icon
View This Table
Fig. 3 Simulated single-shot absorption images for optimum exposure parameters, including fringe removal. The four subplots correspond to the four cases described in section 3. Bright colors correspond to high optical density. Each subplot shows 5 × 5 pixels, with the atom initially located at the center. At larger distances from the initial position the PSF becomes negligibly small. Using Eq. 9 we extract N=(1.32, 1.39, 0.56, 1.26) atoms for subplots a)–d) respectively. These atom numbers are normalized by Napp for each imaging setup.

Depending on the imaging method the effect of recoil blurring on the optical resolution of the imaging system can be significant. For our optimum exposure parameters we find the strongest recoil blurring in case c), where the RMS position of the atoms in the object plane increases by 0.9μm during the imaging pulse. For all other cases considered here recoil blurring is smaller by at least a factor of 3. This blurring should be compared to the base resolution of the imaging system of 1.2μm (Rayleigh criterion). Furthermore the images are discretized by the effective pixel size of 1.3μm in object space.

4. Discussion & conclusion

In this paper we have outlined the theoretical prospects of detecting single rubidium atoms using absorption imaging. We find the optical density of a single atom in single-pass imaging to be ultimately determined only by the numerical aperture of the imaging system with a peak optical density of OD=32NA2(1+38NA2) in the low saturation limit. We have calculated optimum imaging parameters for a number of different cases and could show that the use of a reflective imaging setup can significantly improve the signal-to-noise ratio of absorption imaging. We expect an accuracy of about 82% for the detection of single atoms for the optimum parameters given above in case d), which can be improved to 94% (with a false-positive rate due to zero atoms of 4%) using the reduction in shot noise achieved by the fringe removal algorithm [20

20. S. Whitlock, C. F. Ockeloen, and R. J. C. Spreeuw, “Sub-Poissonian Atom-Number Fluctuations by Three-Body Loss in Mesoscopic Ensembles,” Phys. Rev. Lett. 104, 120402 (2010) [CrossRef] [PubMed] .

]. This compares to a detection accuracy of only 74% (with a false-positive rate due to zero atoms of 15%) for single-pass imaging with the same imaging system. For measuring statistical averages in a regular array of traps further improvements can be made using correlation analysis [20

20. S. Whitlock, C. F. Ockeloen, and R. J. C. Spreeuw, “Sub-Poissonian Atom-Number Fluctuations by Three-Body Loss in Mesoscopic Ensembles,” Phys. Rev. Lett. 104, 120402 (2010) [CrossRef] [PubMed] .

].

A. Absorption signal of a single atom

We start from the simple configuration of a single atom at rest, placed in the centre of the minimum waist of a travelling-wave laser beam, see Fig. 1(a). The laser is taken as monochromatic, Ein(xo,yo,z,t)=12in(o)(xo,yo)uei(kzωt)+c.c., where c.c. denotes complex conjugate, u is a complex polarization unit vector, in(o)(xo,yo) is the amplitude in the object plane with coordinates xo, yo. We assume that the laser is collimated in the object plane, so we can choose in(o)(xo,yo) to be real. The incident intensity is Iin(o)(xo,yo)=[in(o)(xo,yo)]2/2Z0, with Z0=(ε0c)1=μ0/ε0.

We describe the total light field as a sum of this incident plane wave plus a scattered wave. Assuming they have the same polarisation we ignore from now on the polarization vector. The total light field is imaged onto a detector plane where we measure the total intensity Idet(ρ = (x, y)). The scattered wave has both a coherent and an incoherent contribution, which we add separately as follows,
Idet(ρ)=12Z0|in(ρ)+sc(coh)(ρ)|2+Isc(incoh)(ρ)
(11)
=Iin(ρ)+Isc(ρ)+Iif(ρ)
(12)
The total scattered intensity Isc(ρ) is the sum of a coherent and an incoherent contribution, Isc(incoh)(ρ)=sIsc(coh)(ρ) where s is the saturation parameter [21

21. C. Cohen-Tannoudij, J. Dupont-Roc, and G. Grynberg, Atom—Photon Interactions (Wiley-VCH Verlag GmbH, Weinheim, Germany, 1992).

]. The third term is the interference term,
Iif(ρ)=in(ρ)Z0Re(sc(coh)(ρ))
(13)
where sc(coh)(ρ) is complex-valued. The interference intensity is typically negative at the center of the image, as it reduces the total intensity in the image plane.

Next, we introduce the point spread function to describe how the point-like atom is imaged in the detection plane:
sc(coh)(ρ)=𝒜p(ρ)
with complex amplitude 𝒜. We choose the (complex valued) point spread function (PSF) to be normalized as ∫ p(ρ)d2ρ = 1 so that p(ρ) has the dimension 1/area and a = [∫ |p(ρ)|2d2ρ]−1 can be interpreted as the effective area over which the scattered light is distributed in the detection plane.

The absorption signal, Eq. (14), can now be written as
Iabs(ρ)=in(ρ)2Z0[𝒜p(ρ)+c.c.]1+s2Z0|𝒜|2|p(ρ)|2
(18)
In most practical cases we can approximate the probe beam as a nearly plane wave, and p(ρ) as sharply peaked compared to in(ρ). Taking the latter outside the interference integral, Eq. (17), replacing in(ρ) ≈ in(0), we obtain, using Eq. (13),
Iif(ρ)d2ρ=in(0)Re(𝒜)Z0=Psc
(19)

Using the integrals defined above, we can then express the interference term as
Iif(ρ)Psc[Re(p(ρ))Im(𝒜)Re(𝒜)Im(p(ρ))]
(20)
If the atom is in the focal plane of the imaging system we can take p(ρ) to be real-valued so that the Im(p(ρ)) term vanishes. The term also vanishes if the probe beam is exactly on resonance with the atomic transition. In this case Im(p(ρ))/Re(p(ρ)) = tan(arg𝒜)= 0. The approximation of a resonant probe beam is valid as all Doppler shifts are less than 2% of the natural linewidth of the transition for the optimum exposure parameters found in the main part of the paper. In both cases the absorption signal then reduces to
Iabs(ρ)Psc[Re(p(ρ))aχ|p(ρ)|2]
(21)
The leading term in the absorption is simply given by the scattered power. The correction expresses the partial “refilling” of the absorption dip by scattered light. Within the above approximation the absorption signal can also be expressed in terms of the scattering cross section, using Psc = σIin(0)/(1 + s),
Iabs(ρ)Iin(0)=σ1+s[Re(p(ρ))aχ|p(ρ)|2]
(22)

Acknowledgments

We would like to thank Jasper Reinders for help with the simulations and N.J. van Druten for helpful discussions of the manuscript. This work is part of the research programme of the Foundation for Fundamental Research on Matter (FOM), which is part of the Netherlands Organisation for Scientific Research (NWO).

References and links

1.

W. Neuhauser, M. Hohenstatt, P. Toschek, and H. Dehmelt, “Localized visible Bâ{+} mono-ion oscillator,” Phys. Rev. A. 22, 1137–1140 (1980) [CrossRef] .

2.

N. Schlosser, G. Reymond, I. Protsenko, and P. Grangier, “Sub-poissonian loading of single atoms in a microscopic dipole trap.” Nature 411, 1024–7 (2001) [CrossRef] [PubMed] .

3.

K. D. Nelson, X. Li, and D. S. Weiss, “Imaging single atoms in a three-dimensional array,” Nat. Phys. 3, 556–560 (2007) [CrossRef] .

4.

M. Wilzbach, D. Heine, S. Groth, X. Liu, T. Raub, B. Hessmo, and J. Schmiedmayer, “Simple integrated single-atom detector,” Opt. Lett. 34, 259–261 (2009) [CrossRef] [PubMed] .

5.

R. Bücker, A. Perrin, S. Manz, T. Betz, C. Koller, T. Plisson, J. Rottmann, T. Schumm, and J. Schmiedmayer, “Single-particle-sensitive imaging of freely propagating ultracold atoms,” New J. Phys. 11, 103039 (2009) [CrossRef] .

6.

D. Heine, W. Rohringer, D. Fischer, M. Wilzbach, T. Raub, S. Loziczky, X. Liu, S. Groth, B. Hessmo, and J. Schmiedmayer, “A single-atom detector integrated on an atom chip: fabrication, characterization and application,” New J. Phys. 12, 095005 (2010) [CrossRef] .

7.

J. Goldwin, M. Trupke, J. Kenner, a. Ratnapala, and E. a. Hinds, “Fast cavity-enhanced atom detection with low noise and high fidelity.” Nature Commun. 2, 418 (2011) [CrossRef] .

8.

J. Bochmann, M. Mücke, C. Guhl, S. Ritter, G. Rempe, and D. L. Moehring, “Lossless State Detection of Single Neutral Atoms,” Phys. Rev. Lett. 104, 203601 (2010) [CrossRef] [PubMed] .

9.

R. Gehr, J. Volz, G. Dubois, T. Steinmetz, Y. Colombe, B. L. Lev, R. Long, J. Estève, and J. Reichel, “Cavity-Based Single Atom Preparation and High-Fidelity Hyperfine State Readout,” Phys. Rev. Lett. 104, 203602 (2010) [CrossRef] [PubMed] .

10.

D. J. Wineland, W. M. Itano, and J. C. Bergquist, “Absorption spectroscopy at the limit: detection of a single atom.” Opt. Lett. 12, 389–391 (1987) [CrossRef] [PubMed] .

11.

E. W. Streed, A. Jechow, B. G. Norton, and D. Kielpinski, “Absorption imaging of a single atom.” Nature Commun. 3, 933 (2012) [CrossRef] .

12.

M. K. Tey, G. Maslennikov, T. C H Liew, S. A. Aljunid, F. Huber, B. Chng, Z. Chen, V. Scarani, and C. Kurtsiefer, “Interfacing light and single atoms with a lens,” New J. Phys. 11, 043011 (2009) [CrossRef] .

13.

D. A. Smith, S. Aigner, S. Hofferberth, M. Gring, M. Andersson, S. Wildermuth, P. Krüger, S. Schneider, T. Schumm, and J. Schmiedmayer, “Absorption imaging of ultracold atoms on atom chips,” Opt. Express 19, 8471 (2011) [CrossRef] [PubMed] .

14.

V. Y. F. Leung, A. Tauschinsky, N. J. Druten, and R. J. C. Spreeuw, “Microtrap arrays on magnetic film atom chips for quantum information science,” Quantum Inf. Proc. 10, 955–974 (2011) [CrossRef] .

15.

R. Gerritsma, S. Whitlock, T. Fernholz, H. Schlatter, J. A. Luigjes, J.-U. Thiele, J. B. Goedkoop, and R. J. C. Spreeuw, “Lattice of microtraps for ultracold atoms based on patterned magnetic films,” Phys. Rev. A. 76, 033408 (2007) [CrossRef] .

16.

S. Whitlock, R. Gerritsma, T. Fernholz, and R. J. Spreeuw, “Two-dimensional array of microtraps with atomic shift register on a chip,” New J. Phys. 11, 023021 (2009) [CrossRef] .

17.

A. J. E. M. Janssen, “Extended Nijboer–Zernike approach for the computation of optical point-spread functions,” J. Opt. Soc. Am. A. 19, 849 (2002) [CrossRef] .

18.

C. Ockeloen, A. Tauschinsky, R. Spreeuw, and S. Whitlock, “Detection of small atom numbers through image processing,” Phys. Rev. A. 82, 061606 (2010) [CrossRef] .

19.

D. A. Steck, “Rubidium 87 D Line Data,” http://steck.us/alkalidata2 (2010).

20.

S. Whitlock, C. F. Ockeloen, and R. J. C. Spreeuw, “Sub-Poissonian Atom-Number Fluctuations by Three-Body Loss in Mesoscopic Ensembles,” Phys. Rev. Lett. 104, 120402 (2010) [CrossRef] [PubMed] .

21.

C. Cohen-Tannoudij, J. Dupont-Roc, and G. Grynberg, Atom—Photon Interactions (Wiley-VCH Verlag GmbH, Weinheim, Germany, 1992).

OCIS Codes
(020.0020) Atomic and molecular physics : Atomic and molecular physics
(030.4280) Coherence and statistical optics : Noise in imaging systems

ToC Category:
Atomic and Molecular Physics

History
Original Manuscript: February 22, 2013
Revised Manuscript: April 8, 2013
Manuscript Accepted: April 8, 2013
Published: April 16, 2013

Citation
Atreju Tauschinsky and Robert J. C. Spreeuw, "Sensitive absorption imaging of single atoms in front of a mirror," Opt. Express 21, 10188-10198 (2013)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-21-8-10188


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. W. Neuhauser, M. Hohenstatt, P. Toschek, and H. Dehmelt, “Localized visible Bâ{+} mono-ion oscillator,” Phys. Rev. A.22, 1137–1140 (1980). [CrossRef]
  2. N. Schlosser, G. Reymond, I. Protsenko, and P. Grangier, “Sub-poissonian loading of single atoms in a microscopic dipole trap.” Nature411, 1024–7 (2001). [CrossRef] [PubMed]
  3. K. D. Nelson, X. Li, and D. S. Weiss, “Imaging single atoms in a three-dimensional array,” Nat. Phys.3, 556–560 (2007). [CrossRef]
  4. M. Wilzbach, D. Heine, S. Groth, X. Liu, T. Raub, B. Hessmo, and J. Schmiedmayer, “Simple integrated single-atom detector,” Opt. Lett.34, 259–261 (2009). [CrossRef] [PubMed]
  5. R. Bücker, A. Perrin, S. Manz, T. Betz, C. Koller, T. Plisson, J. Rottmann, T. Schumm, and J. Schmiedmayer, “Single-particle-sensitive imaging of freely propagating ultracold atoms,” New J. Phys.11, 103039 (2009). [CrossRef]
  6. D. Heine, W. Rohringer, D. Fischer, M. Wilzbach, T. Raub, S. Loziczky, X. Liu, S. Groth, B. Hessmo, and J. Schmiedmayer, “A single-atom detector integrated on an atom chip: fabrication, characterization and application,” New J. Phys.12, 095005 (2010). [CrossRef]
  7. J. Goldwin, M. Trupke, J. Kenner, a. Ratnapala, and E. a. Hinds, “Fast cavity-enhanced atom detection with low noise and high fidelity.” Nature Commun.2, 418 (2011). [CrossRef]
  8. J. Bochmann, M. Mücke, C. Guhl, S. Ritter, G. Rempe, and D. L. Moehring, “Lossless State Detection of Single Neutral Atoms,” Phys. Rev. Lett.104, 203601 (2010). [CrossRef] [PubMed]
  9. R. Gehr, J. Volz, G. Dubois, T. Steinmetz, Y. Colombe, B. L. Lev, R. Long, J. Estève, and J. Reichel, “Cavity-Based Single Atom Preparation and High-Fidelity Hyperfine State Readout,” Phys. Rev. Lett.104, 203602 (2010). [CrossRef] [PubMed]
  10. D. J. Wineland, W. M. Itano, and J. C. Bergquist, “Absorption spectroscopy at the limit: detection of a single atom.” Opt. Lett.12, 389–391 (1987). [CrossRef] [PubMed]
  11. E. W. Streed, A. Jechow, B. G. Norton, and D. Kielpinski, “Absorption imaging of a single atom.” Nature Commun.3, 933 (2012). [CrossRef]
  12. M. K. Tey, G. Maslennikov, T. C H Liew, S. A. Aljunid, F. Huber, B. Chng, Z. Chen, V. Scarani, and C. Kurtsiefer, “Interfacing light and single atoms with a lens,” New J. Phys.11, 043011 (2009). [CrossRef]
  13. D. A. Smith, S. Aigner, S. Hofferberth, M. Gring, M. Andersson, S. Wildermuth, P. Krüger, S. Schneider, T. Schumm, and J. Schmiedmayer, “Absorption imaging of ultracold atoms on atom chips,” Opt. Express19, 8471 (2011). [CrossRef] [PubMed]
  14. V. Y. F. Leung, A. Tauschinsky, N. J. Druten, and R. J. C. Spreeuw, “Microtrap arrays on magnetic film atom chips for quantum information science,” Quantum Inf. Proc.10, 955–974 (2011). [CrossRef]
  15. R. Gerritsma, S. Whitlock, T. Fernholz, H. Schlatter, J. A. Luigjes, J.-U. Thiele, J. B. Goedkoop, and R. J. C. Spreeuw, “Lattice of microtraps for ultracold atoms based on patterned magnetic films,” Phys. Rev. A.76, 033408 (2007). [CrossRef]
  16. S. Whitlock, R. Gerritsma, T. Fernholz, and R. J. Spreeuw, “Two-dimensional array of microtraps with atomic shift register on a chip,” New J. Phys.11, 023021 (2009). [CrossRef]
  17. A. J. E. M. Janssen, “Extended Nijboer–Zernike approach for the computation of optical point-spread functions,” J. Opt. Soc. Am. A.19, 849 (2002). [CrossRef]
  18. C. Ockeloen, A. Tauschinsky, R. Spreeuw, and S. Whitlock, “Detection of small atom numbers through image processing,” Phys. Rev. A.82, 061606 (2010). [CrossRef]
  19. D. A. Steck, “Rubidium 87 D Line Data,” http://steck.us/alkalidata 2 (2010).
  20. S. Whitlock, C. F. Ockeloen, and R. J. C. Spreeuw, “Sub-Poissonian Atom-Number Fluctuations by Three-Body Loss in Mesoscopic Ensembles,” Phys. Rev. Lett.104, 120402 (2010). [CrossRef] [PubMed]
  21. C. Cohen-Tannoudij, J. Dupont-Roc, and G. Grynberg, Atom—Photon Interactions (Wiley-VCH Verlag GmbH, Weinheim, Germany, 1992).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Figures

Fig. 1 Fig. 2 Fig. 3
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited