# Topologically Robust Transport of Photons in a Synthetic Gauge Field

###### Abstract

Electronic transport in low dimensions through a disordered medium leads to localization. The addition of gauge fields to disordered media leads to fundamental changes in the transport properties. For example, chiral edge states can emerge in two-dimensional systems with a perpendicular magnetic field. Here, we implement a “synthetic” gauge field for photons using silicon-on-insulator technology. By determining the distribution of transport properties, we confirm the localized transport in the bulk and the suppression of localization in edge states, using the “gold standard” for localization studies. Our system provides a new platform to investigate transport properties in the presence of synthetic gauge fields, which is important both from the fundamental perspective of studying photonic transport and for applications in classical and quantum information processing.

1.2 Photons provide a natural and convenient medium to investigate fundamental quantum transport properties Segev et al. (2013); Lagendijk et al. (2009). Using photons, one can selectively excite states, and observe both spectral and spatial response throughout the material, which are challenging tasks in electronic systems. However, experimental efforts studying gauge fields with photons, have been limited to the microwave domain Raghu and Haldane (2008); Wang et al. (2009), while investigations in the optical domain have remained elusive. This is due to the fact that magneto-optic effects — the simplest source of coupling between gauge fields and photons — are extremely weak at optical frequencies. Recently though, there have been a significant number of proposals to synthesize gauge fields for optical photons Hafezi et al. (2011); Umucalilar and Carusotto (2011); Fang et al. (2012); Kraus et al. (2012); Rechtsman et al. (2013a); Khanikaev et al. (2013). In particular, two concurrent experiments showed exemplary signatures of topological edge states through direct imaging Hafezi et al. (2013); Rechtsman et al. (2013b). Here we report the first observation of the robust nature of topologically-protected edge states using an analysis of the statistics of transport properties (transmission and delay). We use a 2D lattice of coupled ring resonators with a synthetic magnetic field, implemented using silicon-on-insulator technology. By considering the distribution of Wigner delay times Chabanov et al. (2000); Chabanov and Genack (2001), we can unambiguously distinguish non-localized diffusive transport in lossy edge states from tunneling through localized bulk states. Finally, we compare the transmission of topologically ordered edge states against the transmission in a topologically trivial one-dimensional system.

Our experiments are performed on a two dimensional lattice of coupled ring resonators Cooper et al. (2010a) (Fig. 1a). The ring resonators are coupled using another set of link rings which are designed to be anti-resonant to the main ring resonators, i.e., the length of the connecting rings is slightly longer than main rings so as to acquire an extra phase shift. The link resonators are spatially shifted, along the y-axis, with respect to the main lattice-site resonators such that transiting photons acquire a phase when hopping along the x-axis at a lattice site with row index Hafezi et al. (2011). Therefore a round trip along any plaquette (consisting of 4 ring and 4 link resonators, see Fig. 1a) results in a total accumulated phase of magnitude with a sign corresponding to the direction (clockwise or counter-clockwise) of travel along the plaquette. Here we only excite and measure the counter-clockwise mode in the main ring resonators, with the input port as indicated in Fig. 1a. This system is equivalent to a uniform synthetic magnetic field with flux penetrating each plaquette of a 2D photon gas, with the tight-binding Hamiltonian

(1) |

where is the coupling rate between the on-site rings, and and are the photon annihilation and creation operators at a main resonator site with indices .

For an infinite lattice, the energy eigenvalues of this Hamiltonian constitute the famous Hofstadter butterfly spectrum. The eigenvalues group into allowed energy bands separated by band-gaps, forming a topological insulator. For a finite lattice, the band-gaps are populated with so-called edge states. The edge states are unidirectional, clockwise (or counter-clockwise) propagating states, with their wavefunction confined to the perimeter of the lattice. We call these long-edge- and short-edge-states (Fig. 1a), respectively, because of the length they travel along the lattice edge from input to output port. These states are in sharp distinction with the eigenstates in the allowed energy bands, which are called bulk states because their eigenfunctions occupy resonators in the bulk of the lattice. In the presence of lattice disorders, such as resonance frequency mismatch and coupling variations, bulk states become localized as the lattice size exceeds the localization length Wen (2007). For our system, bulk states in even the smallest lattice (44) are localized. Edge states, on the other hand, are topologically protected and their wavefunctions are robust against disorder in the lattice. As a result, edge state wavefunctions propagate along the entire edge of the lattice irrespective of the lattice size, although their intensity fall due to absorption and scattering loss. Our goal here is to leverage quantitative measurements of transmission and delay-time to unambiguously demonstrate the robust nature of edge states and distinguish them from bulk states. We show that the edge state transport is diffusive and the delay distribution is gaussian and centered at the average, while for bulk states the delay distribution is asymmetric with the peak value being well below the average, similar to the localization typical in one-dimension Chabanov and Genack (2001); Texier and Comtet (1999); Xu and Wang (2011).

We implemented this system using silicon-on-insulator (SOI) technology, as described in the Supplementary. Fig. 2(a,b) shows the observed transmission and delay spectra at the drop port for eight different 88 lattice size devices. While the spectra differ significantly because of intrinsic fabrication variations in waveguide dimensions, we can already see the first manifestation of robust edge states in the form of two regions with suppressed variance across devices, in both the transmission and delay spectra (red and green shaded). Since edges states are topologically constrained to travel along the lattice edge, device-to-device fabrication variations in system parameters do not affect the edge state wavefunctions as much as they do for bulk states. Edge states therefore show reduced variation. Using numerical modeling including our measured values for disorder (see Supplementary), as shown in Fig. 2(c,d), we can identify these regions as the long edge and the short edge.

We next analyze the delay distribution to distinguish edge state transport from bulk states. This approach provides an unequivocal signature of localization Texier and Comtet (1999); Chabanov and Genack (2001). Fig. 2(e,f) show the measured and simulated delay distributions for the edge and bulk states in 88 lattice sized devices and highlight the remarkable difference between edge and bulk states. For edge states, the delay distribution normalized to its average (in actuality, we used root mean square to allow for negative delay values, see Supplementary) is essentially gaussian with a gaussian width independent of system size. This behavior is characteristic of diffusive transport as seen previously in one-dimensional systems Genack et al. (1999). The bulk state distribution is however, asymmetric with the most probable value being less than the average. This feature is reminiscent of transport governed by localization which also has been observed earlier in the microwave regime for one-dimensional systems Chabanov and Genack (2001). For localized transport, the delay spectrum exhibits spikes (see 2(b)) which manifest in the asymmetric delay statistics. These spikes appear due to resonant tunneling through (delocalized) necklace states which are common to finite-size open systems Pendry (1987). Therefore, even in the presence of loss, delay distribution can clearly differentiate two different regimes of transport in the same photonic system. Our measured results show a good match with numerical modeling. We observe similar behavior for other lattice sizes as well (see Supplementary).

A test to further establish the topologically-protected nature of edge states would be a comparison of the transmission scaling with system size for an edge state against that of a topologically trivial 1D system Cooper et al. (2010a); Notomi et al. (2008), both with similar degrees of disorder. Fabrication-induced disorder in a 1D ring resonator array leads to a spread in the resonance wavelengths of the resonators. This impedes the forward propagation of light, increases back-reflection, i.e., less light is coupled into the array, and hence the transmission at output is reduced Ferrari et al. (2009). Ultimately, as the array length increases, Anderson localization halts transmission of the light Schwartz et al. (2007); Mookherjea et al. (2008); Topolancik et al. (2007); Atlasov et al. (2010). Edge states, on the other hand, are unidirectional and immune to reflection caused by disorder. Therefore, transmission through edge states is expected to be less affected.

Fig. 3a shows the measured average transmission and its standard deviation across a number of chips (95 in total) for the long edge state band in 2D lattice and the mid-band of the 1D array as a function of system size, i.e. the number of resonators travelled from input to output (excluding the link resonators). Transmission in both the long-edge state as well as the 1D system decays exponentially with system size. A linear fit to measured transmission (in dB) in the long edge band gives the decay slope as -0.75(20) dB per ring but for 1D transmission the slope is -0.93(16) dB per ring, where uncertainties represent one sigma standard deviation. Transmission along long edge state can be seen to decay at a slower rate when compared to 1D transport. Simulation results using the experimentally estimated parameters are also presented in the figure. The simulated transmission decay slopes are -0.66(2) and -1.06(5) dB per ring respectively for long edge state and 1D. The experimental and simulation results are seen to agree, given that the number of devices measured for each lattice size is only 8 versus the 5000 realizations for each simulation. In order to differentiate the decay of transmission with system size resulting from resonator losses characterized by , from losses due to disorder - both resulting in exponential attenuation - we plot the simulated result for transmission without disorder (presented as a dashed line). In that zero disorder limit, both the 2D and 1D systems are similarly attenuated by loss with a decay rate of 0.46dB per ring. We therefore observe that disorder affects both 2D and 1D systems, but transport in edge states is less susceptible to disorder.

Fig. 3b shows the measured and simulated average delay and its standard deviation for short-edge and long-edge-state-bands. The measured delay, when plotted against the number of rings on the short and long edges of the lattice, increases linearly with a slope 3.9(9)ps and 5.4(1.0)ps per ring respectively for long and short edge states. The simulated delay slopes are 3.2(2)ps and 4.4(1)ps per ring respectively. Again, the experimental results are in agreement with the simulation. Also shown in the figure, for comparison, is the measured delay in 1D devices. That delay follows the same scaling as the edge states. However, it can be seen that the standard deviation in delay for 1D devices is less than that for edge states. This is contrary to the case of a uniform magnetic field where the standard deviation of delay in edge states remains smaller than in a 1D system Hafezi et al. (2011). Using simulations we have verified that this is due to the fabrication disorder of in the 2D lattice. We further compare the delay distribution for edge states and a 1D system. The normalized delay distribution for long-edge states in a 1515 lattice, and for mid-band and band-ends of a 30 ring 1D array are shown in Fig. 3c. We see that for both, edge states and the mid-band of a 1D array, the transport is diffusive, the distribution is gaussian and the width of the distribution is independent of system size. However, the band ends of the 1D array are localized. Using simulations, we also find that as the array length increases beyond 70 rings, even the mid-band of the 1D system shows localization.

The silicon-on-insulator technology provides a suitable platform to investigate the statistical effects of synthetic gauge fields on various transport properties and to demonstrate the localization of bulk states and robustness of edge states. Such a system could pave the way to investigate the effects of other gauge fields, including those with magnetic monopoles and floquet properties, with or without abelian features on bosonic transport. Moreover, the addition of enhanced optical nonlinearity in these ring resonator structures opens the door to intriguing questions on the nature of solutions for nonlinear transport in systems with topological order.

## References

- Segev et al. (2013) M. Segev, Y. Silberberg, and D. N. Christodoulides, Nature Photon. 7, 197 (2013).
- Lagendijk et al. (2009) A. Lagendijk, v. B. Tiggelen, and D. Wiersma, Physics Today pp. 24–29 (2009).
- Raghu and Haldane (2008) S. Raghu and F. D. M. Haldane, Phys. Rev. A 78, 033834 (2008).
- Wang et al. (2009) Z. Wang, Y. Chong, J. D. Joannopoulos, and M. Soljačić, Nature 461 (2009).
- Hafezi et al. (2011) M. Hafezi, E. A. Demler, M. D. Lukin, and J. M. Taylor, Nature Phys. 7, 907 (2011).
- Umucalilar and Carusotto (2011) R. O. Umucalilar and I. Carusotto, Phys. Rev. A 84, 043804 (2011).
- Fang et al. (2012) K. Fang, Z. Yu, and S. Fan, Nature Photon. 6, 782 (2012).
- Kraus et al. (2012) Y. Kraus, Y. Lahini, Z. Ringel, M. Verbin, and O. Zilberberg, Phys. Rev. Lett. 109, 106402 (2012).
- Rechtsman et al. (2013a) M. C. Rechtsman, J. M. Zeuner, A. Tünnermann, S. Nolte, M. Segev, and A. Szameit, Nature Photon. 7, 153 (2013a).
- Khanikaev et al. (2013) A. B. Khanikaev, S. H. Mousavi, W.-K. Tse, M. Kargarian, A. H. MacDonald, and G. Shvets, Nature Mater. 12, 233 (2013).
- Hafezi et al. (2013) M. Hafezi, S. Mittal, J. Fan, A. Migdall, and J. M. Taylor, Nature Photon. 7, 1001 (2013).
- Rechtsman et al. (2013b) M. C. Rechtsman, J. M. Zeuner, Y. Plotnik, Y. Lumer, D. Podolsky, F. Dreisow, S. Nolte, M. Segev, and A. Szameit, Nature 496, 196 (2013b).
- Chabanov et al. (2000) A. A. Chabanov, M. Stoytchev, and A. Z. Genack, Nature 404, 850 (2000).
- Chabanov and Genack (2001) A. A. Chabanov and A. Genack, Phys. Rev. Lett. 87, 233903 (2001).
- Cooper et al. (2010a) M. L. Cooper, G. Gupta, M. A. Schneider, W. M. J. Green, S. Assefa, F. Xia, Y. A. Vlasov, and S. Mookherjea, Opt. Express 18, 26505 (2010a).
- Wen (2007) X.-G. Wen, Quantum Field Theory of Many-body Systems (Oxford University Press, USA, 2007).
- Texier and Comtet (1999) C. Texier and A. Comtet, Phys. Rev. Lett. 82, 4220 (1999).
- Xu and Wang (2011) F. Xu and J. Wang, Phys. Rev. B. 84, 024205 (2011).
- Genack et al. (1999) A. Genack, P. Sebbah, M. Stoytchev, and B. van Tiggelen, Phys. Rev. Lett. 82, 715 (1999).
- Pendry (1987) J. B. Pendry, J. Phys. C: Solid State Phys. 20, 733 (1987).
- Notomi et al. (2008) M. Notomi, E. Kuramochi, and T. Tanabe, Nature Photon. 2, 741 (2008).
- Ferrari et al. (2009) C. Ferrari, F. Morichetti, and A. Melloni, J. Opt. Soc. Am. B 26, 858 (2009).
- Schwartz et al. (2007) T. Schwartz, G. Bartal, S. Fishman, and M. Segev, Nature 446, 52 (2007).
- Mookherjea et al. (2008) S. Mookherjea, J. S. Park, S.-H. Yang, and P. R. Bandaru, Nature Photon. 2, 90 (2008).
- Topolancik et al. (2007) J. Topolancik, B. Ilic, and F. Vollmer, Phys. Rev. Lett. 99, 253901 (2007).
- Atlasov et al. (2010) K. A. Atlasov, M. Felici, K. F. Karlsson, P. Gallo, A. Rudra, B. Dwir, and E. Kapon, Opt. Express 18, 117 (2010).
- Gifford et al. (2005) D. K. Gifford, B. J. Soller, M. S. Wolfe, and M. E. Froggatt, Appl. Opt. 44, 7282 (2005).
- Cooper et al. (2010b) M. L. Cooper, G. Gupta, M. A. Schneider, W. M. J. Green, S. Assefa, F. Xia, D. K. Gifford, and S. Mookherjea, Opt. Lett. 35, 3030 (2010b).
- van Tiggelen et al. (1999) B. A. van Tiggelen, P. Sebbah, M. Stoytchev, and A. Z. Genack, Phys. Rev. E 59, 7166 (1999).

Acknowledgments: We thank A. Lobos, S. Ganeshan and J. Keutchayan for fruitful discussions. This research was supported by the U.S. Army Research Office Multidisciplinary University Research Initiative award W911NF0910406, and the NSF through the Physics Frontier Center at the Joint Quantum Institute. SF acknowledges financial support from the European Research Council (Advanced Grant SiMoSoMa).

Disclaimer: Certain commercial equipment, instruments or materials are identified in this paper to foster understanding. Such identification does not imply recommendation or endorsement by the National Institute of Standards and Technology, nor does it imply that the materials or equipment are necessarily the best available for the purpose.

Supplementary Material

## .1 Device Fabrication and Measurements

The devices were fabricated at the IMEC foundry using deep-UV projection photolithography Hafezi et al. (2013). The resonator waveguide cross-section is 510 nm in width and 220 nm in height, which allows only a single TE mode to exist in waveguides. The coupling region between all rings consists of a linear waveguide section of 7 m with coupling gap of 180 nm, resulting in a uniform coupling rate . The system is probed using input/output waveguides coupled to the lattice, with a coupling rate (Fig. 1a). Light coupled to the lattice at the input port travels through the lattice and appears at the drop port. The fraction of light which does not couple to the lattice travels to the through port. The light backscattered, due to waveguide surface roughness and reflections in the coupling region, is directed to the back-scattering port. The back-scattered light intensity is about 30dB lower than that observed on the drop port indicating negligible spin-flip disorder in the lattice. For transmission and delay-time measurements, we use an optical vector analyzer (LUNA OVA5000) based on swept wavelength interferometry Gifford et al. (2005). The delay-time (Wigner time) for propagation is then calculated as a derivative of phase with respect to angular frequency. Unlike transport time, Wigner time can be negative for anomalous dispersion regions around a phase jump. The negative delay values also appear in simulations where the delay is calculated similarly.

## .2 Characterization

To characterize the system parameters, we use one ring resonator with two coupling waveguides (add-drop filter - ADF) and one ring resonator with a single coupling waveguide (all-pass filter - APF), designed with the same parameters as those of the rings in the lattice. Transmission measurements were made on 26 single ring ADFs on different chips. Using these measurements , and were measured to be 2.35 GHz, 37.8 GHz and 32.0 GHz with a relative standard deviation of 20%, 4% and 4% respectively. The magnetic flux is designed to be . For 1D devices, the measured bandwidth was less than that expected from simulations with =32 GHz. Therefore, for 1D simulations, a corrected value of =25 GHz was used. A similar procedure was used earlier to estimate for a 1D array of coupled resonators in Cooper et al. (2010b). Fabrication errors also result in a variation of the resonance frequency of the rings in a given lattice. The standard deviation was estimated using transmission measurements on three chips with five all-pass filters each, with a physical separation of the APFs commensurate to rings in a lattice. can then be used to calculate the deviation of optical path length and hence in the link rings. was estimated to be 27.5 GHz and to be 0.1.

## .3 Calibration of spectra

Each transmission and delay spectrum shown in this work is normalized to the corresponding measurement made away from the resonance band, at the through port. The measurements thus normalized give the actual transmission and delay incurred only through the lattice and excludes those in the coupling waveguides and connecting fibers. Because of the intrinsic spread in resonance frequencies resulting from fabrication disorders, the measured and simulated spectra have been shifted along the frequency axis to superpose them. Since the spectra are expected to be disparate in the bulk region Hafezi et al. (2011), we can rely only on the edge state regions to superpose them. For measured spectra, we therefore first do a manual coarse shift to align similar looking features in the expected edge state regions of the spectra. This accounts for across various chips (which is much greater than for a given lattice). Then we analyze the standard deviation of transmission and delay across devices as a function of frequency and find that the edge states are evident as regions with reduced noise. To verify this evidence for edge states and also to align them further, we require an algorithm based on quantitative measurements of transmission (T) and delay-time (). Weighted delay time is one parameter that accounts for both our measurements and has been used extensively to study transport properties in random media Chabanov and Genack (2001); van Tiggelen et al. (1999). For completely random transport, as is the case for bulk states, we expect increased variations in and and hence also in as a function of frequency. On the contrary, transport through the edge states band follows a definite path and should therefore display regions with reduced variance in . We accordingly use the standard deviation of to look for edge state regions and align the measured spectra. Each spectrum is shifted such that given by

(S1) |

where refers to device index, is minimized in the designated edge state bands. The bandwidth of the long edge is found to be 10 GHz, independent of the device size while the short edge is wider (12.5 GHz - 19 GHz). For our analysis, we fix the bandwidth of the short- and long-edge regions to be 10 GHz for all devices. For simulated spectra, we follow exactly the same protocol except for the course shift which is not required. The spectra have not been shifted along y axis and there is no re-scaling of the spectrum.

## .4 Numerical simulations with noise

For numerical simulations of the transmission and delay spectrum we use coupled mode analysis treating the ring resonators as lumped elements Hafezi et al. (2011). In the tight binding approximation, the resonators are coupled only to their nearest neighbors with a coupling rate . As shown in the supplementary section of Hafezi et al. (2013), includes the response of connecting rings. The Hamiltonian of the system is characterized by the resonance frequency of each ring resonator , its coupling rate to its nearest neighbors , the magnetic phase acquired when hopping along nearest neighbors along x-axis and , the coupling rate to probe waveguide.

To include lattice disorder into this Hamiltonian, we impose random variations on each of the parameters with a gaussian probability distribution around the mean. For each numerical realization of the lattice, each resonator ring has noise added to it resonance frequency, coupling rate to neighboring resonators and also to the magnetic phase acquired in hopping along x-axis. To get mean transmission and delay, we then average the results over 5000 realizations for each device type. Required simulation parameters and their deviations, characterizing the system and its disorder, have been measured using multiple add-drop filters as described above.

The fact that the measured transmission scales exponentially and the delay scales linearly with the number of rings on the edge of the lattice, and both match well with the simulated results, reinforce our claim that the low noise areas are in fact the short and long edge regions. From simulation, we found that the main disorder terms affecting transmission in a 2D lattice are and variations in the otherwise uniform magnetic field i.e. , whereas for a 1D array only the first term is applicable since there is no magnetic field. In the absence of , the transmission in the edge state would be even closer to the dashed line with no disorder. The short-edge transmission in our system was however consistently found to be much lower (8 dB for 6x6 devices) than expected using simulations, but it tends to match simulation results for bigger sized devices. We expect this to be the result of some systematic problem with our fabrication process which couples less light to the short edge at the input port and hence produces a non-zero intercept on the transmission axis. Using through port data, we verified the reduction in coupling efficiency for the short edge band, for all devices.

## .5 Spectrum and delay statistics for 1515-lattice-sized devices

Fig. S1 shows the measured and simulated transmission and delay statistics for 1515-lattice-sized devices. As was seen in Fig. 2, the transmission and delay spectra show two regions with reduced noise. The delay statistics is also similar to what is observed for 88-lattice-sized devices. Edge state transport is diffusive whereas the bulk state is localized.

## .6 Gaussian fits to delay distribution

Fig. S2 plots Gaussian fits to measured and simulated delay distributions for 88-lattice-sized devices. The bulk states deviate systematically from the Gaussian fit, towards longer delay times.

## .7 Localization in 2D Lattice

Fig. S3 (a,b) shows the simulated long edge and bulk state intensity for an 88 lattice, with disorder, averaged over 50 realizations. In the presence of disorder, the long edge state spans the complete edge of the lattice, whereas the bulk states are localized near the input port. As the system size increases, as shown in Fig. S3(c,d), the edge state wavefunction still extends across the lattice edge, but the bulk state is again localized with localization extent independent of lattice size.

## .8 Localization in a 1D array

Fig. S4 shows the simulated delay statistics for a 1D array with 70 rings. The probability distribution is asymmetric indicating the onset of localization.

## .9 Measured spectrum for 1D devices with 10 resonator rings

Fig. S5 shows the measured and simulated transmission and delay spectrum for eleven, 1D devices with 10 resonator rings. The device to device variations in transmission and delay are more or less independent of frequency.

## .10 Device yield

Table I shows the detailed device yield. For 2D devices, the device yield was found to be 23%. For 1D devices (other than 50 rings), the device yield was 100%. For 50 ring devices, the yield was 50%. The devices with a very noisy or attenuated spectrum were considered as bad.

Device size | Worked | Bad | Did not scan | Total |
---|---|---|---|---|

66 | 7 | 28 | 0 | 35 |

88 | 8 | 27 | 0 | 35 |

1010 | 9 | 26 | 0 | 35 |

1515 | 8 | 27 | 0 | 35 |

1818 | 8 | 27 | 0 | 35 |

21 | 11 | 0 | 21 | 35 |

101 | 15 | 0 | 20 | 35 |

201 | 11 | 0 | 24 | 35 |

301 | 12 | 0 | 23 | 35 |

501 | 6 | 7 | 22 | 35 |