## Talbot Effect of orbital angular momentum lattices with single photons

S. Schwarz, 1, 2, ∗ C. Kapahi, 1, 2 R. Xu, 1, 3 A. R. Cameron, 1, 2 D. Sarenac, 1, 2 J. P. W. MacLean, 1, 2 K. B. Kuntz, 1, 2 D. G. Cory, 1, 4, 5, 6 T. Jennewein, 1, 2 K.J. Resch, 1, 2 and D.A. Pushin 1, 2, †

1

Institute for

Quantum Computing, University of Waterloo, Waterloo, ON, Canada, N2L 3G1

2

Department of Physics &amp; Astronomy, University of Waterloo, Waterloo, ON, Canada, N2L 3G1

3

Department of

Electrical

&amp; Computer Engineering,

University of Waterloo, Waterloo, ON, Canada, N2L 3G1

4 Department of Chemistry, University of Waterloo, Waterloo, ON, Canada, N2L3G1

5 Perimeter Institute for Theoretical Physics, Waterloo, ON, Canada, N2L2Y5

6 Canadian Institute for Advanced Research, Toronto, ON, Canada, M5G 1Z8

The self-imaging, or Talbot Effect, that occurs with the propagation of periodically structured waves has enabled several unique applications in optical metrology, image processing, data transmission, and matter-wave interferometry. In this work, we report on the first demonstration of a Talbot Effect with single photons prepared in a lattice of orbital angular momentum (OAM) states. We observe that upon propagation, the wavefronts of the single photons manifest self-imaging whereby the OAM lattice intensity profile is recovered. Furthermore, we show that the intensity at fractional Talbot distances is indicative of a periodic helical phase structure corresponding to a lattice of OAM states. This phenomenon is a powerful addition to the toolbox of orbital angular momentum and spin-orbit techniques that have already enabled many recent developments in quantum optics.

The Talbot Effect [1] is a near-field diffraction phenomenon whereby periodic phase and amplitude modulations are self-imaged due to free-space propagation. In accordance with Fresnel diffraction [2], replicas of a periodic transverse intensity profile reappear after a specific propagation distance known as the Talbot length. The Talbot Effect has been demonstrated in numerous areas of research involving linear and nonlinear optical waves [3-5], single photons [6, 7], x-rays [8], matter-waves [912], exciton polaritons [13], and Bose-Einstein condensates [14]. The Talbot Effect has a diverse array of applications in optical metrology [15], imaging processing [16], and lithography [17-19], with potential in data transmission [20].

Here we consider the Talbot Effect manifested by lattices of orbital angular momentum (OAM) states. The OAM degree of freedom of light has garnered significant interest in various fields ranging from optical manipulation and high-bandwidth communication [21-24] to quantum information processing [25, 26]. In addition to the photonic applications, OAM beams have been extended to neutrons [27-29] and electrons [30, 31].

The Talbot Effect has been considered with classical light as well as OAM lattices [32-38]. In this Letter, we report the first demonstration of the Talbot Effect with single photons prepared in a lattice of OAM states. The extension of the Talbot Effect to single photons and OAM techniques offers the possibility of utilizing quantum information processing protocols, such as remote state preparation, to leverage quantum communication advantages [39]. Furthermore, self-imaging has potential applications in implementing quantum logic operations

∗ sacha.schwarz@uwaterloo.ca

† dmitry.pushin@uwaterloo.ca

as qudits may be encoded in the transverse spatial profile of single photons [40, 41].

We consider spin-orbit states described by the wavefunction

$$\text{ific} \quad | \Psi \rangle = A ( r, \phi ) \left [ \cos \left ( \frac { \pi r } { d } \right ) | R \rangle + i e ^ { i \ell \phi } \sin \left ( \frac { \pi r } { d } \right ) | L \rangle \right ], \ ( 1 )$$

where ( r, φ ) are the cylindrical coordinates, glyph[lscript] specifies the OAM number, d is the distance in which the polarization state performs a full rotation on the Poincar´ e sphere, | R 〉 and | L 〉 denote the right and left circular polarization states, and A r, φ ( ) denotes the envelope. A lattice of spin-orbit states can be obtained by passing circularly polarized light through pairs of birefringent linear gradients whose optical axes are perpendicular to each other [35, 42]. This may be derived by considering the Suzuki-Trotter expansion of the operator ˆ U which generates the spin-orbit state of Eq. (1) when acting on an input state | R 〉 , where

$$\hat { U } = e ^ { i \frac { \pi } { d } [ x \hat { \sigma } _ { x } + y \hat { \sigma } _ { y } ] } = \lim _ { N \to \infty } ( e ^ { i \frac { \pi x } { N d } \hat { \sigma } _ { x } } e ^ { i \frac { \pi y } { N d } \hat { \sigma } _ { y } } ) ^ { N }. \quad ( 2 )$$

Truncating the expansion to N terms, the operators in Eq. (2) can be realized by sets of perpendicular birefringent gradients with the general form

$$\hat { U } _ { x } = e ^ { i \frac { \pi } { a } ( x - x _ { 0 } ) \hat { \sigma } _ { x } }, \ \hat { U } _ { y } = e ^ { i \frac { \pi } { a } ( y - y _ { 0 } ) \hat { \sigma } _ { y } }, \quad ( 3 )$$

where the origin of the gradients is given by ( x , y 0 0 ), ˆ σ x,y are Pauli matrices, and where a = Nd . It was shown in Ref. [35] that linear gradients of Eq. (3) may be implemented via 'Lattice of Optical Vorticies' (LOV) prism pairs. A LOV prism pair consists of two wedge-shaped birefringent prisms where the optical axis of the first prism is along the wedge incline direction and that of the second is offset by 45 ◦ [35]. By sending a photon in

FIG. 1. Simulated intensity distributions in both the xy and yz planes, where the beam propagates along z . (a) Right-handed circularly polarized light is sent through N = 2 sets of LOV prism pairs, which yields a lattice of ring-shaped intensity structures when filtered with respect to the left-handed circular polarization, i.e., I ( x, y ) = |〈 L | | Ψ N =2 LOV 〉 | 2 at propagation distance z = 0. Note that here the Gaussian intensity envelope α x, y ( ) with beam waist w 0 = 3 a is added. (b) By plotting the yz intensity at x = a/ 4 (indicated in (a) with the dash dotted white line) we recover what is known as the Talbot carpet. (c) xy intensity patterns at specific propagation distances z . (d-e) The Talbot carpet and the xy intensity cross sections when the phase structure of the initial beam is removed. This demonstrates the effect of the OAM lattice phase structure on the intensity profile at the fractional Talbot distances.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01016v2-figure-b64c48e4-363f-4b14-924c-3d256258ac03.jpg)

the right circular polarization state | R 〉 through N sets of LOV prism pairs, we prepare the state

$$\left | \Psi _ { L O V } ^ { N } \right \rangle = \alpha ( x, y ) ( \hat { U } _ { x } \hat { U } _ { y } ) ^ { N } \left | R \right \rangle, \quad \quad ( 4 ) \quad \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$$$

tance results in the Fourier transform of the initial periodic pattern. The Fraunhofer distance (far-field) is given by z F = 8 w /λ 2 0 , where w 0 is the beam waist. In our setup z F ≈ 166 m.

where α x, y ( ) describes the incoming Gaussian beam envelope with beam waist w 0 . The periodic nature of polarization rotations enables the linear gradients to prepare a two-dimensional lattice of spin-orbit states.

Filtering on one circular polarization state prepares a periodically structured intensity distribution with a lattice spacing of a = λ (∆ n tan( θ )) -1 , where ∆ n and θ are the birefringence and the incline angle of the LOV prism pairs respectively, and λ is the wavelength. In our experiment we use N = 2 LOV prism pairs and we filter on | L 〉 to obtain an initial intensity distribution of the form

$$I ( x, y ) & = | \langle L | \, | \Psi _ { L O V } ^ { N = 2 } \rangle \, | ^ { 2 } \\ & = | \alpha ( x, y ) | ^ { 2 } \cos ^ { 2 } \left [ \frac { \pi x } { a } \right ] \cos ^ { 2 } \left [ \frac { \pi y } { a } \right ] \, \begin{array} { c } E ( x, \\ \text{but} \\ \text{to} \\ \end{array} \, \\ \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \end{array}$$

which is depicted in Fig. 1(a). This periodic beam structure imprinted by the LOV prism pairs sets up conditions required for the Talbot Effect. The transmitted light interferes in such a way that after a distance z T = 2 a /λ 2 , the initial periodic intensity pattern reappears. The same intensity distribution also appears at half the distance, z T / 2, but with spatial shifts ∆ a = a/ 2 along the x -and y -directions. Furthermore, the intensity distribution at propagation distances much larger than the Talbot dis-

Theory predicts the same self-imaging phenomenon for single photons as well. We describe the free-space propagation of single photons by a complex-valued transverse field distribution E x,y ( ) convoluted with the Fresnel propagator

$$K _ { F } ( x, y, z ) = \frac { e ^ { i k z } } { i \lambda z } \exp \left [ \frac { i k } { 2 z } ( x ^ { 2 } + y ^ { 2 } ) \right ], \quad ( 6 )$$

where k is the wavevector. The field E x,y ( ) at position z is evaluated via

$$\text{form} \quad & E ( x, y, z ) = \frac { e ^ { i k z } } { i \lambda z } \iint d x ^ { \prime } d y ^ { \prime } \, E ( x ^ { \prime }, y ^ { \prime }, 0 ) \, e ^ { \frac { i k } { 2 z } [ ( x - x ^ { \prime } ) ^ { 2 } + ( y - y ^ { \prime } ) ^ { 2 } ] }. \\ & \text{row} \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdots$$

Fig. 1(a) shows the simulated transverse intensity distribution, I ( x, y ) = |〈 L ψ | N =2 LOV 〉| 2 , before beam propagation. Fig. 1(b) and Fig. 1(d) depict the intensity distribution in the yz -planes at x = a/ 4 for the initial states 〈 L ψ | N =2 LOV 〉 and |〈 L ψ | N =2 LOV 〉| , respectively. Fig. 1(c) and Fig. 1(e) illustrate the intensity distribution in the xy -planes for specific propagation distances. We observe that the initial phase profile defines the transverse intensity pattern at fractional Talbot distances. Furthermore, it can be observed that the OAM phase structure induces an asymmetry between the intensity distributions at propagation distances { z T / 8, z T / 4,3 z T / 8 } and { 7 z T / 8, 3 z T / 4,5 z T / 8 } .

The experimental setup is schematically depicted in Fig. 2. Degenerate photon pairs are prepared using typeII spontaneous parametric down-conversion in a Sagnac

FIG. 2. Schematic of the experimental setup. Correlated photon pairs are generated via type-II spontaneous parametric downconversion in a Sagnac interferometer and coupled into single mode fibers (SMF). A singles rate of 18 kHz and a coincidence rate of 1 5 kHz are measured after the SMF. After propagating through a 30 m long fiber, the signal photon is sent through a . telescope with 8 3x magnification, . N = 2 sets of LOV prism pairs and a polarization filter. The free-space propagation z , can be varied via different flip mirror combinations. The signal photons are then imaged onto an intensified electron-multiplying CCD (emICCD), triggered by the detection of the corresponding idler. The imaging arrangement in the detection unit consists of a telescope with 4x demagnification ( f 3 and f 4 lenses) followed by a single-lens ( f 5 ) that images the beam onto the detection plane of the emICCD.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01016v2-figure-e937ee49-ebe6-4f45-8c24-23a4cd581bce.jpg)

interferometer. We pump a 10 mm long periodicallypoled KTP crystal (PPKTP) with a continuous wave diode laser (404 8 nm) to produce correlated photon pairs . centered at λ SP = 810 8 nm with a spectral bandwidth of . 0 4 nm. With the pump horizontally polarized, we mea-. sure the second-order correlation function at zero time delay to be g (2) (0) = 0 0251 . ± 0 0011, . implying that two-photon events in any coincidence time window are around 1% [43]. Note that a diagonal polarized pump would offer the ability to generate a polarization entangled target state, however, here we herald the signal by means of measuring the idler (see Ref. [44] for further details). The outputs of the Sagnac interferometer are coupled into two single-mode fibers, which allow for a distinct separation of signal and idler. The signal photons are sent through a telescope to magnify the beam by a factor of 8 3, followed by . N = 2 sets of LOV prism pairs. This configuration prepares a lattice of spin-orbit states where one of the polarization states is coupled to glyph[lscript] = 1. Higher values of glyph[lscript] may be achieved by employing a setup with more LOV prism pairs, while negative values of glyph[lscript] may be achieved by changing the input polarization state [35].

The polarization state of the signal photon is prepared using a manual fiber polarization controller, polarizing beamsplitter (PBS), half wave plate (HWP) and quarter wave plate (QWP). After transmission through the LOV prism pairs, the signal is filtered with respect to left-handed or right-handed circulary polarized light using a QWP. The free-space propagation of the OAM lattice is then analyzed via an arrangement of flip mirrors that effectively change the propagation distance z before measurement. The single photon detection unit consists of a telescope to demagnify the beam by a factor of 4 ( f 3 and f 4 lenses in Fig. 2) and a gated intensified electronmultiplying CCD (emICCD PI-Max4: 1024 EMB). The telescope is followed by a single lens ( f 5 lens in Fig. 2) which images the plane immediately following the telescope.

The idler is detected by an avalanche photodiode with no polarization filtering, which acts as a trigger for the emICCD to herald the single photon state. We use a 30 m spool of single-mode fiber to delay the detection of the signal with respect to the idler to accommodate for the delays from the triggering electronics. We set the delay time between the idler and signal photon for each propagation distance z and use the emICCD camera to align the coincidence window of 3 ns.

In addition to the single photon setup, we couple light from a linearly polarized laser diode (central wavelength λ LD = 813 4 nm) into the signal . channel in order to compare images generated by single photons versus laser diode light. We measure the intensity profile using a conventional CCD camera (Coherent LaserCam-HR II) at the same positions as the single photon images captured by the emICCD.

In Fig. 3 we present simulated and measured beam profiles at fractional Talbot distances. Although the theoretical Talbot length is z T = 16 m, the propagation distances in the experimental setup were increased by a constant offset of 0 85 m in order to account for the three . lens system in the detection unit [45]. Tab. I lists the experimental distances, Z exp , which effectively correspond to the theoretical distances, Z theo . The diode images

FIG. 3. Simulated and experimental self-images at different fractional Talbot lengths. We measure the two-dimensional intensity profile I ( x, y ) = |〈 L | Ψ N =2 LOV 〉 | 2 at positions z ∈ Z exp . In the simulation, we multiply a Gaussian beam envelope with the same beam waist w 0 as in the experiment (i.e., w 0 = (4 1 . ± 0 05) mm) to account for features occurring due to finite . beam sizes when propagating along the z -axis. For comparison, we couple light from a laser diode into the signal channel, and measure corresponding self-images at the same positions. Good qualitative agreement is found between the simulated and observed profiles. We measure a lattice spacing of a exp = (0 573 . ± 0 012) mm which corresponds to (2 229 . . ± 0 037) mm . before the demagnification by the three lens system in the detection unit. From the simulation, we extract a lattice spacing of a sim = (0 577 . ± 0 010) mm and (2 301 . . ± 0 031) mm, respectively. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01016v2-figure-b37347a6-2279-4680-a5a3-d179698cfc5d.jpg)

were also measured at distances z ∈ Z exp . The central wavelength difference of | λ LD -λ SP | = 2 6 nm corre-. sponds to a change in Talbot length z T of only ∼ 5 cm. The LOV prisms were realigned in the transverse plane to obtain the most pronounced doughnut structures with the diode laser.

TABLE I. Experimental propagation distances Z exp which correspond to the fractional Talbot distances Z theo , and single photon signal-to-noise ratio (SNR). The SNR is given by the ratio of the average signal to the standard deviation of the background. In the third (fourth) column, we list the SNR calculated from raw (post-processed) images.

| Z theo    | Z exp    | Measured SNR   | Post-processed SNR   |
|-----------|----------|----------------|----------------------|
| 0         | 0 . 71 m | 0 . 584        | 240 . 377            |
| z T / 8   | 2 . 86 m | 0 . 547        | 181 . 988            |
| z T / 4   | 4 . 85 m | 0 . 113        | 102 . 514            |
| 3 z T / 8 | 6 . 87 m | 0 . 159        | 126 . 298            |
| z T / 2   | 8 . 86 m | 0 . 259        | 264 . 755            |

The observed intensity profiles are measured with a total exposure time of about 1 hour and are processed using background subtraction and an adaptive two-dimensional Gaussian image filter. Including the quadratic phase profiles of the imaging lenses in the simulation yields good agreement between theoretical and observed lattice spacing. For instance, in the case of single photons, we extract from the transverse intensity distribution at z = 0 071 m a separation between two nearest-neighbor . lattice sites of a exp = (0 573 . ± 0 012) mm from exper-. imental data and a sim = (0 577 . ± 0 010) mm from the . simulation. Additionally, at half Talbot distance z T / 2, the expected half period shift ∆ a can be evaluated by comparing the effective pixel positions of the lattice sites at propagation distance z = 0 071 m with the pixel posi-. tions at z = z T / 2 yielding ∆ a exp = (0 273 . ± 0 015) mm . and ∆ a sim = (0 279 . ± 0 014) mm, respectively. .

The robustness of the Talbot Effect with a lattice of OAM states is demonstrated by the good qualitative agreement between simulation, single photon, and diode laser images in Fig. 3. The SNR decreases with larger

distances, but is increased again depending on the intensity pattern complexity. In Tab. I, we present the SNR before and after the imaging post-processing for different propagation distances. However, it can be noted that the self-imaging property of this beam can be seen clearly in the similarity between images taken at distances z = { 0 , z T / 2 } , with the correct spatial shift. Images at z = { z T / 8 3 , z T / 8 } show an orientation about each lattice site that appears counter-clockwise in z T / 8 and clockwise in 3 z T / 8. These features are indicative of the OAM state in each lattice site, as shown in Fig. 1(c). Furthermore, gaps in the outermost rings of the z T / 2 image can be mitigated by using a beam containing more lattice sites.

In this work we demonstrated and analyzed the Talbot Effect with single photons prepared in a lattice of OAM states. Heralded single photons are sent through N = 2 sets of LOV prism pairs and their transverse two-dimensional intensity distribution are measured at various fractional Talbot lengths. The propagation of

- [1] H. T. E. F.R.S., The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 9 , 401 (1836).
- [2] J. Wen, Y. Zhang, and M. Xiao, Adv. Opt. Photon. 5 , 83 (2013).
- [3] R. Iwanow, D. A. May-Arrioja, D. N. Christodoulides, G. I. Stegeman, Y. Min, and W. Sohler, Phys. Rev. Lett. 95 , 053902 (2005).
- [4] W. B. Case, M. Tomandl, S. Deachapunya, and M. Arndt, Optics Express 17 , 20966 (2009).
- [5] Y. Zhang, J. Wen, S. N. Zhu, and M. Xiao, Phys. Rev. Lett. 104 , 183901 (2010).
- [6] X.-B. Song, H.-B. Wang, J. Xiong, K. Wang, X. Zhang, K.-H. Luo, and L.-A. Wu, Phys. Rev. Lett. 107 , 033902 (2011).
- [7] S. Deachapunya, S. Srisuphaphon, P. Panthong, T. Photia, K. Boonkham, and S. Chiangga, Opt. Express 24 , 20029 (2016).
- [8] F. Pfeiffer, M. Bech, O. Bunk, P. Kraft, E. F. Eikenberry, C. Br¨ onnimann, C. Gr¨nzweig, u and C. David, Nature Materials 7 , 134 (2008).
- [9] M. S. Chapman, C. R. Ekstrom, T. D. Hammond, J. Schmiedmayer, B. E. Tannian, S. Wehinger, and D. E. Pritchard, Phys. Rev. A 51 , R14 (1995).
- [10] L. Deng, E. W. Hagley, J. Denschlag, J. E. Simsarian, M. Edwards, C. W. Clark, K. Helmerson, S. L. Rolston, and W. D. Phillips, Phys. Rev. Lett. 83 , 5407 (1999).
- [11] S. Wu, E. Su, and M. Prentiss, Phys. Rev. Lett. 99 , 173201 (2007).
- [12] F. Pfeiffer, C. Gr¨nzweig, O. Bunk, G. Frei, E. Lehmann, u and C. David, Phys. Rev. Lett. 96 , 215505 (2006).
- [13] T. Gao, E. Estrecho, G. Li, O. A. Egorov, X. Ma, K. Winkler, M. Kamp, C. Schneider, S. H¨fling, A. G. Truscott, o and E. A. Ostrovskaya, Phys. Rev. Lett. 117 , 097403 (2016).
- [14] V. Makhalov and A. Turlapov, Phys. Rev. Lett. 122 , 090403 (2019).

structured wavefronts is calculated in the near-field and shows good agreement with experimental results. We observe that the initial phase profile defines the transverse intensity pattern at various propagation distances, and thus the Talbot carpet. Future work will scrutinize the connection between OAM and Talbot physics as a new characterization tool. Another avenue of exploration includes the addition of path entangled OAM lattices using the Talbot Effect and the OAM degree of freedom. Other avenues of exploration include the addition of path entangled OAM lattices and the implementation of quantum logic using the Talbot Effect and the OAM degree of freedom.

This research was supported in part by the Canadian Excellence Research Chairs (CERC) program, the Natural Sciences and Engineering Research Council of Canada (NSERC), the NSERC Discovery program, Canada Research Chairs, Industry Canada and the Canada Foundation for Innovation (CFI), Ontario Research Fund (ORF), and the Canada First Research Excellence Fund (CFREF).

- [15] E. T. F. Rogers, J. Lindberg, T. Roy, S. Savo, J. E. Chad, M. R. Dennis, and N. I. Zheludev, Nature Materials 11 , 432 (2012).
- [16] A. Wang, P. Gill, and A. Molnar, Appl. Opt. 48 , 5897 (2009).
- [17] A. Isoyan, F. Jiang, Y. C. Cheng, F. Cerrina, P. Wachulak, L. Urbanski, J. Rocca, C. Menoni, and M. Marconi, Journal of Vacuum Science &amp; Technology B: Microelectronics and Nanometer Structures Processing, Measurement, and Phenomena 27 , 2931 (2009).
- [18] H. H. Solak, C. Dais, and F. Clube, Opt. Express 19 , 10686 (2011).
- [19] Y. Liu, X. Zhang, Y. Huang, J. Zhang, W. Hofmann, Y. Ning, and L. Wang, Optik 183 , 579 (2019).
- [20] J. Azana and M. A. Muriel, Appl. Opt. 38 , 6700 (1999).
- [21] H. He, M. E. J. Friese, N. R. Heckenberg, and H. Rubinsztein-Dunlop, Phys. Rev. Lett. 75 , 826 (1995).
- [22] M. E. J. Friese, J. Enger, H. Rubinsztein-Dunlop, and N. R. Heckenberg, Phys. Rev. A 54 , 1593 (1996).
- [23] J. Wang, J.-Y. Yang, I. M. Fazal, N. Ahmed, Y. Yan, H. Huang, Y. Ren, Y. Yue, S. Dolinar, M. Tur, and A. E. Willner, Nature Photonics 6 , 488 (2012).
- [24] R. Zhang, L. Li, Z. Zhao, G. Xie, G. Milione, H. Song, P. Liao, C. Liu, H. Song, K. Pang, A. N. Willner, B. Lynn, R. Bock, M. Tur, and A. E. Willner, Opt. Lett. 44 , 691 (2019).
- [25] M. F. Andersen, C. Ryu, P. Clade, V. Natarajan, A. Vaziri, K. Helmerson, and W. D. Phillips, Phys. Rev. Lett. 97 , 170406 (2006).
- [26] M. Erhard, R. Fickler, M. Krenn, and A. Zeilinger, Light: Science &amp;Amp; Applications 7 , 17146 (2018).
- [27] C. W. Clark, R. Barankov, M. G. Huber, M. Arif, D. G. Cory, and D. A. Pushin, Nature 525 , 504 (2015).
- [28] D. Sarenac, M. G. Huber, B. Heacock, M. Arif, C. W. Clark, D. G. Cory, C. B. Shahi, and D. A. Pushin, Opt. Express 24 , 22528 (2016).
- [29] D. Sarenac, C. Kapahi, W. Chen, C. W. Clark, D. G.

Cory, M. G. Huber, I. Taminiau, K. Zhernenkov, and D. A. Pushin, Proceedings of the National Academy of Sciences (2019), 10.1073/pnas.1906861116.

- [30] M. Uchida and A. Tonomura, Nature 464 , 737 (2010).
- [31] B. J. McMorran, A. Agrawal, I. M. Anderson, A. A. Herzing, H. J. Lezec, J. J. McClelland, and J. Unguris, Science 331 , 192 (2011).
- [32] J. Courtial, R. Zambrini, D. M. R, and M. Vasnetsov, Opt. Express 14 , 938 (2006).
- [33] S. Vyas and P. Senthilkumaran, Applied optics 46 , 2893 (2007).
- [34] G.-X. Wei, L.-L. Lu, and C.-S. Guo, Optics Communications 282 , 2665 (2009).
- [35] D. Sarenac, D. G. Cory, J. Nsofini, I. Hincks, P. Miguel, M. Arif, C. W. Clark, M. G. Huber, and D. A. Pushin, Phys. Rev. Lett. 121 , 183602 (2018).
- [36] H. Gao, Y. Li, L. Chen, J. Jin, M. Pu, X. Li, P. Gao, C. Wang, X. Luo, and M. Hong, Nanoscale 10 , 666 (2018).
- [37] J. Hu, C.-S. Br`s, e and C.-B. Huang, Opt. Lett. 43 , 4033 (2018).
- [38] S. Rasouli and D. Hebri, J. Opt. Soc. Am. A 36 , 800 (2019).
- [39] J. T. Barreiro, T.-C. Wei, and P. G. Kwiat, Phys. Rev. Lett. 105 , 030407 (2010).
- [40] O. J. Far´ ıas, F. de Melo, P. Milman, and S. P. Walborn, Phys. Rev. A 91 , 062328 (2015).
- [41] K. Sawada and S. P. Walborn, Journal of Optics 20 , 075201 (2018).
- [42] D. Sarenac, J. Nsofini, I. Hincks, M. Arif, C. W. Clark, D. G. Cory, M. G. Huber, and D. A. Pushin, New Journal of Physics 20 , 103012 (2018).
- [43] M. J. Stevens, S. Glancy, S. W. Nam, and R. P. Mirin, Opt. Express 22 , 3244 (2014).
- [44] L. Vermeyden, X. Ma, J. Lavoie, M. Bonsma, U. Sinha, R. Laflamme, and K. J. Resch, Phys. Rev. A 91 , 012120 (2015).
- [45] J. W. Goodman, Roberts &amp; Co. Publishers, (2005).
Sacha Schwarz Connor Kapahi Ruoxuan Xu Andrew R. Cameron Dusan Sarenac Jean-Philippe W. MacLean Katanya B. Kuntz David G. Cory Thomas Jennewein Kevin J. Resch Dmitry A. Pushin
Talbot Effect of orbital angular momentum lattices with single photons
The self-imaging, or Talbot Effect, that occurs with the propagation of
periodically structured waves has enabled several unique applications in
optical metrology, image processing, data transmission, and matter-wave
interferometry. In this work, we report on the first demonstration of a Talbot
Effect with single photons prepared in a lattice of orbital angular momentum
(OAM) states. We observe that upon propagation, the wavefronts of the single
photons manifest self-imaging whereby the OAM lattice intensity profile is
recovered. Furthermore, we show that the intensity at fractional Talbot
distances is indicative of a periodic helical phase structure corresponding to
a lattice of OAM states. This phenomenon is a powerful addition to the toolbox
of orbital angular momentum and spin-orbit techniques that have already enabled
many recent developments in quantum optics.