-
During the last decade, the growing interest in quantum simulation has fostered the development of several techniques for implementing effective electromagnetic fields in systems of neutral particles1, 2. In this vein, artificial gauge fields (AGFs) have been widely used in photonics to control light dynamics3-5, emulating the effect of electromagnetic fields on charged particles. Moreover, AGFs have also allowed the exploration of a plethora of phenomena stemming from their close connection to topological phases of matter6-9 (see Ozawa et al.10 for a recent review). Typically, these AGFs are introduced either by geometric manipulation4, 5 or by time-dependent modulation11-13. While in Wu et al.14, wavepackets carrying orbital angular momentum (OAM) were used to create edge states in crystalline topological insulators, here, we experimentally demonstrate that an AGF in the form of an effective magnetic flux can be induced using Laguerre–Gauss light beams carrying OAM15. Specifically, to prove the existence of this flux, we show how Aharonov-Bohm (AB) caging naturally appears when OAM modes with a specific topological charge are injected into cylindrical optical waveguides arranged in a diamond chain configuration16, 17.
AB caging, which was originally studied in the context of two-dimensional electronic systems, is a single-particle localisation effect arising from the interplay between the lattice geometry and a magnetic flux. More specifically, a constant magnetic flux modifies the phase relations of wavepackets, resulting in a destructive interference effect that binds the modes. Thus, it enables one to halt all propagation by controlling the flux. This phenomenon, which can be interpreted in terms of quantum interference18, 19, has been predicted to occur20-22 and experimentally verified23, 24 in photonic structures implementing AGFs. Unlike the previous photonic proposals based on geometric manipulation20-22, we show in this work how non-zero-energy flat bands, which are responsible for the caging effect, can be naturally and deliberately achieved by injecting light carrying OAM instead of fabricating a new sample. Therefore, our proposal enables the study of the effect of AGFs in photonic lattices just by selecting the topological charge of the input beam. In this context, our proposal differs from related works where the intrinsic angular momentum, i.e., the polarisation of the input beam, instead of the extrinsic one, i.e., the OAM, was used as the AGF switching mechanism25. Moreover, this method also allows access to different topological regimes without the need to fabricate different structures or employ high intensities, as is the case for topological phase transitions realised via non-linear optics26.
To experimentally visualise the AB caging effect induced by OAM modes, we fabricate photonic lattices composed of direct laser written optical waveguides27 arranged in a diamond chain configuration, as displayed in Fig. 1a. The unit cell j is composed of three waveguides $ \left( {S_j = A_j,B_j,C_j} \right) $ forming a triangle with a central angle θ. Each cylindrical waveguide sustains OAM modes of the form28
Fig. 1 Lattice structure and optical waveguides.
a Schematic representation of the structure composed of identical cylindrical waveguides arranged in a diamond chain configuration. Each unit cell j hosts three waveguides forming a triangle with central angle θ. The distances between waveguide centres are$ s_j \equiv A_j,\;B_j,\;C_j $ ,$ d_{A_j - B_j} = d_{A_j - C_j} \equiv d $ and$ d_{B_j - C_j} = 2d\sin (\theta /2) $ . The blue arrows indicate the couplings. b Refractive index profile of the waveguides, defined by ncore = 1.548, nclad = 1.540 and waveguide radius R = 1.9 μm. Field intensity of the$ d_{A_j - A_{j + 1}} = 2d\cos (\theta /2) $ (green) and$ \ell = 0 $ (red) modes, where$ \ell = 1 $ is the propagation constant of mode$ \beta _\ell $ , k0 = 2π/λ0 is the vacuum wavenumber and λ0 is the light wavelength in vacuum. c Numerically calculated coupling strengths for separation distances d = 5 μm, 5.5 μm, 6 μm, 6.5 μm, 7 μm and 7.5 μm using λ0 = 700 nm. In particular, c0 (crosses) accounts for the coupling between the$ \ell $ modes, and c1 (circles) and c2 (squares) account for the coupling between the$ \ell = 0 $ modes with the same or opposite circulation directions, respectively. The dashed and solid lines correspond to the exponential fittings of$ \ell = 1 $ ,$ c_0\left( d \right) \approx K_0\exp ( - \kappa _0d) $ and$ c_1(d) \approx K_1\exp ( - \kappa _1d) $ , where K0 = 387 mm−1, κ0 = 1.17 μm−1, K1 = 19.39 mm−1, κ1 = 0.52 μm−1, K2 = 56.25 mm−1 and κ2 = 0.59 μm−1. The inset in c shows c2/c1 with respect to the separation distance d$ c_2(d) \approx K_2\exp ( - \kappa _2d) $ $$ {\mathrm{\Psi }}_{S_j}^{ \pm \ell }\left( {r_{S_j},\phi _{S_j},z} \right) = \psi _{S_j}^\ell \left( {r_{S_j}} \right)e^{ \pm i\ell \left( {\phi _{S_j} - \phi _0} \right)}e^{ - i\beta _\ell z} $$ (1) where $ \ell = 0,1,2, \ldots $ is the topological charge, ± accounts for positive and negative circulation of the phase front, $ \psi _{S_j}^\ell \left( {r_{S_j}} \right) $ is the radial mode profile given by the Bessel functions15, $ (r_{S_j},\phi _{S_j}) $ are the polar coordinates with respect to the centre of each waveguide Sj in the transverse plane, z is the propagation direction, ϕ0 is an arbitrary phase origin, and $ \beta _\ell $ is the propagation constant of mode $ \ell $. Moreover, while between fundamental modes $ (\ell = 0) $, there is only one coupling amplitude $ c_{0,0} \equiv c_0 $, between OAM modes $ (\ell \, \ne \, 0) $ with the same or opposite circulation directions, there are two coupling amplitudes $ c_{\ell ,\ell } \equiv c_1 $ and $ c_{\ell , - \ell } \equiv c_2e^{i2\ell \phi _0} $29. In particular, as a proof of concept, we restrict our implementation to the $ \ell = 0 $ and $ \ell = 1 $ modes by properly engineering the refractive index contrast and the width of the step-index profile presented in Fig. 1b. In this case, between the $ \ell = 1 $ modes with the same or opposite circulation directions, there are two coupling amplitudes $ c_{1,1} \equiv c_1 $ and $ c_{1, - 1} \equiv c_2e^{i2\ell \phi _0} $29. Therefore, when dealing with OAM modes, complex coupling amplitudes between modes with different circulation directions appear naturally. The different coupling strengths $ c_0,\;c_1 $ and c2 are presented in Fig. 1c (see Supplementary I for details on the calculations). Specifically, we set the phase origin ϕ0 along the $ A_j \leftrightarrow C_j $ direction such that $ c_{1, - 1} = c_2 $ is real in this direction, while $ c_{1, - 1} = c_2e^{ - i2\ell {\mathtt θ}} $ is complex along the $ A_j \leftrightarrow B_j $ direction. In particular, we fix θ = π/2, which allows the coupling between modes propagating in next-nearest neighbour waveguides to be neglected30 (see Supplementary II for a detailed discussion). Moreover, for this specific angle, a relative phase difference of π between the $ c_{1, - 1} $ couplings in the $ A_j \leftrightarrow C_j $ and $ A_j \leftrightarrow B_j $ directions appears. This phase difference introduces a π flux into the plaquettes that opens an energy gap between the dispersive bands, as discussed in detail in the following.
Assuming periodic boundary conditions, the bulk band structure for the $ \ell = 0 $ modes consists of one flat and two dispersive bands (Fig. 2a), with energies given by20
Fig. 2 Energy band structure.
a Band structure of the considered diamond chain lattice for consisting of two dispersive bands$ \ell = 0 $ and$ E_ - ^0(k) $ (dotted green lines) and one zero-energy flat band$ E_ + ^0(k) $ (solid red line). Band structure of the considered diamond chain lattice for$ E_0^0\left( k \right) $ when b c2/c1 = 2 and c c2/c1 = 1. In b and c, each band has a two-fold degeneracy, i.e.,$ \ell = 1 $ (dashed line),$ E_ - ^1(k) \equiv E_1\left( k \right) = E_2\left( k \right) $ (solid line) and$ E_0^1(k) \equiv E_3\left( k \right) = E_4\left( k \right) $ (dashed line)$ E_ + ^1(k) \equiv E_5\left( k \right) = E_6\left( k \right) $ $$ E_0^0\left( k \right) = 0,\;E_ \pm ^0\left( k \right) = \pm 2c_0\sqrt {1 + {\mathrm{cos}}(k\sqrt 2 d)} $$ (2) where k is the quasi-momentum and $ \sqrt 2 d $ is the lattice constant. On the other hand, as presented in Fig. 2b, the band structure for $ \ell = 1 $ is composed of six energy bands, i.e., three bands with a twofold degeneracy (positive and negative circulation)16
$$ E_0^1\left( k \right) = 0,E_ \pm ^1\left( k \right) \\= \pm 2\sqrt {\left( {c_1^2 + c_2^2} \right) + \left( {c_1^2 - c_2^2} \right){\mathrm{cos}}(k\sqrt 2 d)} $$ (3) The main difference between the energy bands in the two cases is the existence of an energy gap for $ \ell = 1 $, which is absent for $ \ell = 0 $, indicating the presence of an AGF. By performing a basis rotation (Supplementary III), the original diamond chain can be decoupled into two identical chains with three energy bands and a π flux through the plaquettes that opens the energy gap16. Moreover, as illustrated in Fig. 2c, in the $ c_2/c_1 \to 1 $ limit, the dispersive bands $ E_ \pm ^1 \to \pm 2\sqrt 2 c_1 $ become flat, and the associated supermodes are localised in the $ A_j,B_j,B_{j + 1},C_j $ and $ C_{j + 1} $ waveguides. Therefore, if one excites Aj with a $ \ell = 1 $ mode, then the injected intensity will oscillate between the central and four surrounding waveguides, as predicted by the AB caging effect (Supplementary III).
To experimentally demonstrate AB caging using OAM modes, we excite a central waveguide Aj using modes with and without OAM and compare the resulting dynamics. We fabricate several samples with seven unit cells with different total lengths (ranging from z = 250 μm to z = 1000 μm) and extract the output pattern intensities. A scheme of the samples is depicted in Fig. 3. First, as displayed in Fig. 3a, we inject a mode with $ \ell = 1 $ and negative circulation into A4 (see Supplementary IV for complementary results). The injected mode spreads to the four surrounding waveguides at z = 250 μm (Fig. 3b) and recombines in the central waveguide at z = 500 μm (Fig. 3c). This spreading and recombination effect can be observed a second time at 750 μm (Fig. 3d) and 1000 μm (Fig. 3e). Even though we implement the model with $ c_2/c_1 \approx 2 $ due to experimental restrictions on the total size of the samples, we measure two full oscillations of the AB caging effect. Since the dispersive bands are not totally flat, light propagates into waveguides A3 and A5 during the second oscillation, and part of the intensity escapes from the cage (Fig. 3d, e). Additionally, although we try to excite the donut mode with negative circulation (see the input beam in Fig. 3a), the propagating mode has a lobe-shaped intensity (Fig. 3b–e) corresponding to a superposition of donut modes with positive and negative circulation. This lobe-shaped mode appears due to a slight ellipticity of the fabricated waveguides and the influence of the surrounding waveguides (see Supplementary IV for a complementary discussion). Nevertheless, since the propagation of the $ \ell = 1 $ modes with positive and negative circulation results in the same flux, the observed AB caging is the same for any superposition of both types of circulation, i.e., a lobe-shaped mode (Supplementary III). In contrast, the $ \ell = 0 $ mode injected into A4 only spreads transversally as it evolves along the propagation direction, and no caging is observed in Fig. 3f–j.
Fig. 3 Aharonov-Bohm caging effect.
Experimentally observed input and output intensities obtained by exciting the A4 waveguide using the OAM mode with and negative circulation at a z = 0 μm, b z = 250 μm, c z = 500 μm, d z = 750 μm and e z = 1000 μm and with$ \ell = 1 $ at f z = 0 μm, g z = 250 μm, h z = 500 μm, i z = 750 μm and j z = 1000 μm. Note that the image in a is taken before entering the sample. The diamond chain lattice is composed of seven unit cells, i.e., 21 waveguides with radius R = 1.9 μm and nearest-neighbour separation d = 5.5 μm. The wavelength used is λ0 = 700 nm. The intensity distribution in each figure is normalised to the maximum intensity value of the corresponding figure$ \ell = 0 $ Finally, we compare the experimental observations of the light dynamics with numerical calculations. Figure 4 shows the intensity extracted at the output port from the A4 waveguide and its associated cage formed by $ A_4,\;B_4,C_4,B_5,\;C_5 $. In Fig. 4a, we can observe how the experimentally measured intensity maxima in A4 associated with the caging phenomenon occur around z = 500 μm and 1000 μm, in agreement with finite-difference method (FDM) simulations. On the other hand, in Fig. 4b, one can observe the standard decay of the intensity in A4 when the $ \ell = 0 $ mode is injected. Moreover, we also compute the light dynamics for longer distances using coupled-mode equations (Supplementary I and II). In Fig. 4c, one can observe how for $ \ell = 1 $, the first and second intensity maxima in A4 have ~60% and 10% of the injected intensity, respectively, which can be increased by reducing the difference between c1 and c2 (see inset of Fig. 1c). For example, for $ c_2/c_1 \approx 1.25, $ i.e., d = 15 μm, the first and second maxima increase up to 97% and 80%, respectively, achieving 100% in the flat-band limit. However, larger separations between waveguides require longer samples, which were not feasible in our experiments. Alternatively, for $ \ell = 0 $, the intensity in A4 exponentially decays independent of the waveguide separation, confirming the different origins of the oscillations. Finally, note that the agreement between the experimental results obtained with different samples and waveguides and injecting modes with both types of circulation shown in Supplementary IV confirms the robustness of the AB caging effect since each measurement includes slight parameter variations (see Supplementary V for more details).
Fig. 4 Light dynamics along the propagation direction.
Intensity extracted from waveguide A4 (circles) and from the cage formed by (squares) normalised to the intensity extracted from the entire lattice as a function of the propagation distance z when the a$ A_4,\;B_4,C_4,B_5,\;C_5 $ mode and b$ \ell = 1 $ mode is injected into waveguide A4. The results shown in a are an average of the intensities extracted for$ \ell = 0 $ with positive and negative circulation. The circles and squares correspond to the experimentally extracted intensities, while the lines correspond to the best-fitting curve of the simulated results obtained using FDM numerical techniques. The error bars associated with the experimental data are estimated taking into account a refractive index error of Δn = ±0.001 in the fabrication process. Intensity propagating in waveguide A4 numerically calculated using coupled-mode equations as a function of z when the c$ \ell = 1 $ mode and d$ \ell = 1 $ mode is injected into waveguide A4. The solid lines correspond to the case with d = 5.3 μm, i.e.,$ \ell = 0 $ , while the dashed lines correspond to d = 15 μm, i.e.,$ c_2/c_1 \approx 2 $ . Note that the simulations were performed considering N = 7 unit cells and λ0 = 700 nm with a correction of Δd = –0.2 μm with respect to the expected experimental distance d = 5.5 μm. This difference may originate from slight variations in the position during the writing process (±0.05 μm) and small changes in the refractive index contrast$ c_2/c_1 \approx 1.25 $ In summary, we demonstrated that an artificial gauge field in the form of an effective magnetic flux could be induced in a photonic lattice by exploiting the orbital angular momentum carried by light beams. Specifically, we demonstrated the appearance of this synthetic flux by experimentally measuring the photonic analogue of the Aharonov-Bohm caging effect for an arrangement of direct laser written cylindrical waveguides in a diamond chain configuration. Using this structure, we showed how an energy gap is opened between the dispersive bands of the system when light carrying OAM is injected, analogous to the effect produced by an artificial gauge field20. Moreover, we proved how non-zero-energy flat bands, which yield the AB caging effect, can be achieved by properly tuning the geometry of the unit cells and the separation between waveguides. The agreement between the dynamics revealed by the coupled-mode equations, the FDM simulations and the experiments confirms the validity of the presented model, which constitutes a step towards accessing different topological regimes in an active way by controlling the input states. Moreover, the inherently infinite dimensionality of OAM modes15 can be potentially exploited to increase the transmission capacity by using mode multiplexing31, paving the way towards combining integrated spatial multiplexing32 with topological protection10.
Artificial gauge field switching using orbital angular momentum modes in optical waveguides
- Light: Science & Applications 9, Article number: (2020)
- Received: 14 April 2020
- Revised: 26 July 2020
- Accepted: 11 August 2020 Published online: 28 August 2020
doi: https://doi.org/10.1038/s41377-020-00385-6
Abstract: The discovery of artificial gauge fields controlling the dynamics of uncharged particles that otherwise elude the influence of standard electromagnetic fields has revolutionised the field of quantum simulation. Hence, developing new techniques to induce these fields is essential to boost quantum simulation of photonic structures. Here, we experimentally demonstrate the generation of an artificial gauge field in a photonic lattice by modifying the topological charge of a light beam, overcoming the need to modify the geometry along the evolution or impose external fields. In particular, we show that an effective magnetic flux naturally appears when a light beam carrying orbital angular momentum is injected into a waveguide lattice with a diamond chain configuration. To demonstrate the existence of this flux, we measure an effect that derives solely from the presence of a magnetic flux, the Aharonov-Bohm caging effect, which is a localisation phenomenon of wavepackets due to destructive interference. Therefore, we prove the possibility of switching on and off artificial gauge fields just by changing the topological charge of the input state, paving the way to accessing different topological regimes in a single structure, which represents an important step forward for optical quantum simulation.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article′s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article′s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.