# Enhancing quantum transport in a photonic network using controllable decoherence

## Abstract

Transport phenomena on a quantum scale appear in a variety of systems, ranging from photosynthetic complexes to engineered quantum devices. It has been predicted that the efficiency of quantum transport can be enhanced through dynamic interaction between the system and a noisy environment. We report the first experimental demonstration of such environment-assisted quantum transport, using an engineered network of laser-written waveguides, with relative energies and inter-waveguide couplings tailored to yield the desired Hamiltonian. Controllable decoherence is simulated via broadening the bandwidth of the input illumination, yielding a significant increase in transport efficiency relative to the narrowband case. We show integrated optics to be suitable for simulating specific target Hamiltonians as well as open quantum systems with controllable loss and decoherence.

Recent research into photosynthetic antenna complexes has shown evidence of coherence in excitonic energy transport Engel et al. (2007); Panitchayangkoon et al. (2010); Collini et al. (2010); Wong et al. (2012), despite the noisy cellular environment in which such complexes are found. Indeed, environmental decoherence has been credited with increasing the efficiency of transport through these systems, an effect known as environment-assisted quantum transport (ENAQT) Rebentrost et al. (2009) or dephasing-assisted transport Plenio and Huelga (2008); Chin et al. (2010). While ENAQT has been the subject of many theoretical studies—whether in the photosynthetic context Mohseni et al. (2008); Cao and Silbey (2009); Caruso et al. (2009); Wu et al. (2010); Hoyer et al. (2010); Fassioli and Olaya-Castro (2010); Wu et al. (2013); Cleary and Cao (2013); León-Montiel et al. (2014) or in other nanoscale transport systems Semião et al. (2010); Scholak et al. (2011a, b); Lim et al. (2014)—and despite its potential for improving transport in artificial quantum systems, it has so far never been directly observed.

We use an integrated photonic simulator to demonstrate the first implementation of ENAQT. Our simulator was fabricated using femtosecond-laser direct writing, which allows waveguides to be drawn directly into glass using a focused pulsed laser. This permits the creation of three-dimensional waveguide arrays, as well as precision and repeatability in engineering interactions Szameit et al. (2007); Szameit and Nolte (2010); Marshall et al. (2009); Owens et al. (2011); Corrielli et al. (2013); Poulios et al. (2014); Crespi et al. (2015). We used control over the wavelength and bandwidth of the guided light to simulate effective decoherence, thereby enhancing transport efficiency.

We consider a single excitation on a network of coupled sites, governed by a tight-binding Hamiltonian Rebentrost et al. (2009)

(1) |

where denotes the excitation being localised at site , the energy of that site, and the coupling between sites and . Although ENAQT can occur on an ordered lattice where all the energies are equal Kassal and Aspuru-Guzik (2012), transport enhancement was first explained in disordered systems, which we consider here.

We are interested in the efficiency of transport from a particular initial site to a particular target site, where the excitation is trapped. In the case of a photosynthetic complex (figure 1a), trapping describes the transfer of excitons to a reaction center, where they drive charge separation. It can be modelled as irreversible coupling of the target site to a sink at rate (figure 1b). The efficiency is then the probability of finding the exciton in the sink after some particular time.

ENAQT occurs when decoherence increases the trapping probability over the fully coherent case. Decoherence results from coupling of a quantum system to inaccessible degrees of freedom. For example, in photosynthetic antenna complexes, the energies of chromophores are coupled to molecular vibrations; tracing out this environment results in decoherence in the excitonic subspace.

In the absence of decoherence, energetic disorder tends to localise the wavepacket through processes such as destructive interference or Anderson localisation Anderson (1958); Lahini et al. (2008); Segev et al. (2013), thus preventing it from reaching the target. Since these are coherent processes, they are diminished by decoherence, possibly resulting in enhancement of transport efficiency to the target site. An alternative but ultimately equivalent view of ENAQT is that eigenstates of are stationary, making it difficult to reach the target if the initial and target sites differ in energy. Incoherent processes, however, permit transitions between the eigenstates of , yielding greater mobility.

The first theoretical explorations of ENAQT focused on the case where the decoherence takes the form of site-independent, Markovian, pure dephasing Plenio and Huelga (2008); Rebentrost et al. (2009). Although the decoherence in our simulation is neither site-based nor Markovian, this is not an obstacle to ENAQT; Mohseni et al. (2014); Chen and Silbey (2011) transport efficiency can be enhanced as long as the decoherence allows population transfer between otherwise-stationary eigenstates.

We simulate ENAQT in an array of coupled single-mode optical waveguides obeying the equation Perets et al. (2008); Longhi (2009)

(2) |

where the light is propagating in the -direction, is a creation operator for a photon in waveguide at position , and and are respectively the propagation constants of the waveguides and the couplings between them. The former are determined by the waveguides’ refractive index profiles, while the latter also depend on the separations between them. Light propagation governed by this Schrödinger-like equation directly simulates evolution under , with replacing and the replacing . We can thus simulate different Hamiltonians by controlling the number, position, and refractive indices of the waveguides.

The intrinsic stability of laser-written waveguide arrays renders decoherence challenging to simulate. One approach is to stochastically modulate the index of refraction along each waveguide Levi et al. (2012). Though every realization with a particular longitudinal index profile will be fully coherent, decoherence can be simulated by averaging over the recorded optical outputs from many realizations in post-processing. This approach was recently used to simulate decoherence-enhanced navigation of a maze Caruso et al. (2015).

By contrast, we simulate decoherence by averaging over the results from a single array illuminated with many optical wavelengths. Although each individual wavelength propagates through the waveguide array coherently, decoherence can be achieved using broadband illumination and a single output intensity measurement which does not resolve wavelength. In other words, the wavelength degree of freedom is traced out yielding a partially mixed state. Similar approaches are well established in optics, where for instance thick birefringent quartz plates followed by wavelength-insensitive measurements have been used to decohere the polarization states of single photons Kwiat et al. (2000); Xu et al. (2010).

As an example of our approach to simulating decoherence, we consider two uncoupled waveguides and that have a difference in their propagation constants at a wavelength . As light propagates along the waveguides, it will accumulate a phase difference between them. Broadband illumination can then be seen to cause effective decoherence on a length scale comparable to the illumination coherence length: Any initial coherence between the waveguides will decay:

(3) |

where is the normalised first-order temporal correlation function of the light, which—being proportional to the Fourier transform of the spectrum—decays to zero faster for spectrally broader illumination. The strength of the decoherence can be quantified as the inverse of the optical coherence length, , and is usually proportional to the FWHM bandwidth . For a uniform distribution centered at ,

(4) |

In the case of coupled waveguides, wavelength dependence affects couplings in addition to the propagation constants. The resulting decoherence will therefore not only have characteristics of pure dephasing (as in equation (3)), but will also include off-diagonal terms. However, this kind of decoherence can also result in ENAQT. Indeed, the off-diagonal terms imply that our decoherence could, in principle, be used to simulate ENAQT in ordered systems Kassal and Aspuru-Guzik (2012), but the effect would be weaker because of the absence of the pure-dephasing contribution.

For our simulation we chose a network of four waveguide sites—arranged as in figure 1c—because it is one of the smallest systems in which ENAQT is possible and because it can give significant enhancements even with relatively weak decoherence. Waveguide 1 was the input, waveguide 3 the target, and the sink consisted of a long linear array of tightly coupled waveguides. The coupling between waveguides in this sink was significantly higher than between the four main waveguides, so that any light entering the sink from waveguide 3 was largely transported away Longhi (2006); Biagioni et al. (2008); Delanty et al. (2012). The sink need only be long enough to prevent light reflected from the far end from returning into the main simulation waveguides.

We need to choose propagation constants and separations among the four main waveguides in order to best approximate this Hamiltonian:

(5) |

where all parameters depend on the wavelength of the input light. Due to their wide separations, couplings between non-neighbouring sites are negligible ( of neighbouring-site couplings). Our simulator is designed so that, at our central simulation wavelength ,

(6) |

In this case, one of the eigenstates of has no support on site 3, , while the remaining three eigenstates all have substantial support on site 3. Because cannot couple to the sink—at least at —the maximum trapping efficiency at infinite time is .

Considering the wavelength dependence of provides a different way to see the decohering effects of broadband illumination. Due to the dependence of and on the wavelength, equation (6) only holds at a particular wavelength . At other wavelengths, will have some support on waveguide 3 and thus be able to couple to the sink, increasing the efficiency above , as shown in figure 1d. Unlike in other examples of ENAQT Rebentrost et al. (2009), the efficiency increases monotonically with the strength of decoherence, meaning there is no optimal level of decoherence in this model.

Based on measurements of couplings and propagation constants in isolated pairs of waveguides (see Methods), we selected the following design parameters: , , and (see figure 1c). Figure 2a shows numerical modelling of light propagation given these parameters. Although these were designed for a center wavelength of 800 nm, variations in the implementation of the waveguide parameters resulted in equation 6 being satisfied at nm.

We measured the efficiency using narrowband light (less than 1 nm bandwidth and always horizontally polarised for consistency) from a tuneable Ti:sapphire laser (Spectra-Physics Tsunami) in quasi-cw mode (figure 2b). The output was imaged using a custom-built magnifying telescope and the optical power was measured using a large-area power-meter after isolating either the system or sink waveguides using a variable slit. Examples of the output distribution are given in figure 2c, showing the significant difference between illumination at and an off-center wavelength. Figure 3a shows the measured efficiency (fraction of light output in the sink modes) for wavelengths ranging from 745 to 835 nm.

ENAQT—shown in figure 3b—is the average enhancement in efficiency over the spectral band of interest, relative to the efficiency at ,

(7) |

where the average over is taken over the top-hat distribution on . As predicted, it is an increasing function of the bandwidth, i.e., of the decoherence strength . The highest ENAQT observed was at a bandwidth of nm. We thus demonstrated that coherent transport in a coupled, statically disordered system can be enhanced through decoherence.

The theoretical prediction in figure 3b contains no free parameters. It is a simulation of the dynamics under the simulated Hamiltonian , together with trapping from site 3 at the rate discussed in Methods and pure dephasing between waveguide and the other three waveguides at rate . The disagreement between theory and experiment is small considering the number of possible contributing factors. These include the off-diagonal decoherence when the waveguides are coupled, the fact that the trapping is not perfectly exponential, errors in the measurements of the coupling constants, optical losses, and error in satisfying equation 6. As an example, the shaded band in figure 3b represents the error that would arise if deviated from by up to 10%.

In our experiment, the magnitude of ENAQT was limited by the maximum achievable decoherence, which at cm was small compared to the inverse of the propagation length. We were limited by two components of equation (4): the tunability of our laser limited , while was limited (via equation (6)) by the need to keep small enough to stay in the tight-binding approximation. These are not fundamental limitations, and future improvements that increase would result in significantly larger transport enhancement; as shown in figure 1d, the transport efficiency of this model can get arbitrarily close to 1 for sufficiently long propagation distances and decoherence strengths. Stronger decoherence would also allow ENAQT to be observed in networks that are less sensitive to decoherence than our model.

Our results demonstrate that integrated photonics are well-suited for simulating open quantum systems, capable of implementing a disordered target Hamiltonian with controllable loss. The technique of using broadband excitation to introduce tuneable levels of decoherence ensures that photonics can simulate open quantum systems; previously, low intrinsic noise in photonic devices rendered decoherence difficult to realize in integrated optics, particularly without averaging over results from many different device realizations. Our approach will not only improve the controllability of photonic quantum simulation, but also aid in the experimental optimisation of transport in other engineered quantum systems.

## Methods

Our experimental setup is shown in figure 2b. The waveguides were fabricated in high-purity fused silica (Corning 7980) using a laser direct-write technique whereby Ti:sapphire laser pulses are tightly focused into the sample, which is then translated in three dimensions to yield continuous regions of positive refractive index change which act as waveguides Szameit and Nolte (2010). To obtain single-mode waveguides with the desired propagation and coupling characteristics, laser pulses of duration 150 fs, energy 400 µJ, central wavelength 800 nm, and repetition rate 100 kHz were focused 400 µm below the surface using a 40 microscope objective and translated at 75 mm/min (for waveguides other than 4).

Permanent index changes result in the focus, yielding elliptical, vertically oriented, single-mode waveguides with modes approximately µm in size at 800 nm. Waveguide 3 and the sink waveguides are in a plane parallel to the surface, while angles of 120° between this plane and the other system waveguides minimise non-neighbour coupling. Because couplings between elliptical waveguides are dependent on angular orientation Szameit et al. (2007), we determined the couplings by writing pairs of waveguides oriented at the specified angles. Couplings as a function of separation were determined by writing the pairs at different separations and measuring the output intensities after a known propagation length when only one is optically excited.

The final waveguide separations were chosen to ensure that the tight-binding approximation is maintained, that can be matched by , that , and that the number of sink modes remains manageable. The system-sink population transfer can be modelled as a constant effective rate Longhi (2006) if . In our case, although meant that the decay was not perfectly exponential, it was nevertheless effectively irreversible, considering that sufficiently many sink modes were present to prevent reflection from the far end and back into the main waveguides.

For waveguide 4, the writing translation speed was decreased to increase the propagation constant. In a pair of waveguides with mismatched propagation constants, the maximum power transfer depends only on the ratio . We used this fact to determine the translation speed decrease necessary to set at 800 nm. Due to the measurement uncertainty (about ) and day-to-day variability in the waveguide writing process, several arrays were written with slightly different writing speeds for waveguide 4 to ensure that one set would yield in our desired wavelength range. In the best sample, the writing speed was decreased by 12%.

## Acknowledgements

DNB, AAZ, MAB, AF, AGW, and IK were supported by Australian Research Council (ARC) Centres of Excellence for Engineered Quantum Systems (CE110001013) and Quantum Computation and Communication Technology (CE110001027). RH, MG, SN, and AS thank the German Ministry of Education and Research (Center for Innovation Competence program, 03Z1HN31), the Thuringian Ministry for Education, Science and Culture (Research group Spacetime, 11027-514), the Deutsche Forschungsgemeinschaft (NO462/6-1), and the German-Israeli Foundation for Scientific Research and Development (1157-127.14/2011). AF and IK acknowledge ARC Discovery Early Career Researcher Awards (DE130100240 and DE140100433).

### References

- G. S. Engel, T. R. Calhoun, E. L. Read, T.-K. Ahn, T. Mančal, Y.-C. Cheng, R. E. Blankenship, and G. R. Fleming, Nature 446, 782 (2007).
- G. Panitchayangkoon, D. Hayes, K. A. Fransted, J. R. Caram, E. Harel, J. Wen, R. E. Blankenship, and G. S. Engel, Proc. Natl. Acad. Sci. 107, 12766 (2010).
- E. Collini, C. Y. Wong, K. E. Wilk, P. M. G. Curmi, P. Brumer, and G. D. Scholes, Nature 463, 644 (2010).
- C. Y. Wong, R. M. Alvey, D. B. Turner, K. E. Wilk, D. A. Bryant, P. M. G. Curmi, R. J. Silbey, and G. D. Scholes, Nature Chem. 4, 396 (2012).
- P. Rebentrost, M. Mohseni, I. Kassal, S. Lloyd, and A. Aspuru-Guzik, New J. Phys. 11, 033003 (2009).
- M. B. Plenio and S. F. Huelga, New J. Phys. 10, 113019 (2008).
- A. Chin, A. Datta, F. Caruso, S. Huelga, and M. Plenio, New J. Phys. 12, 065002 (2010).
- M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru-Guzik, J. Chem. Phys. 129, 174106 (2008).
- J. Cao and R. J. Silbey, J. Phys. Chem. A 113, 13825 (2009).
- F. Caruso, A. W. Chin, A. Datta, S. F. Huelga, and M. B. Plenio, J. Chem. Phys. 131, 105106 (2009).
- J. Wu, F. Liu, Y. Shen, J. Cao, and R. Silbey, New J. Phys. 12, 105012 (2010).
- S. Hoyer, M. Sarovar, and K. B. Whaley, New J. Phys. 12, 065041 (2010).
- F. Fassioli and A. Olaya-Castro, New J. Phys. 12, 085006 (2010).
- J. Wu, R. J. Silbey, and J. Cao, Phys. Rev. Lett. 110, 200402 (2013).
- L. Cleary and J. Cao, New J. Phys. 15, 125030 (2013).
- R. León-Montiel, I. Kassal, and J. P. Torres, J. Phys. Chem. B 118, 10588 (2014).
- F. L. Semião, K. Furuya, and G. J. Milburn, New J. Phys. 12, 083033 (2010).
- T. Scholak, F. De Melo, T. Wellens, F. Mintert, and A. Buchleitner, Phys. Rev. E 83, 021912 (2011a).
- T. Scholak, T. Wellens, and A. Buchleitner, J. Phys. B 44, 184012 (2011b).
- J. Lim, M. Tame, K. H. Yee, J.-S. Lee, and J. Lee, New J. Phys. 16, 053018 (2014).
- A. Szameit, F. Dreisow, T. Pertsch, S. Nolte, and A. Tünnermann, Opt. Express 15, 1579 (2007).
- A. Szameit and S. Nolte, J. Phys. B 43, 163001 (2010).
- G. D. Marshall, A. Politi, J. C. F. Matthews, P. Dekker, M. Ams, M. J. Withford, and J. L. O’Brien, Opt. Express 17, 12546 (2009).
- J. O. Owens, M. A. Broome, D. N. Biggerstaff, M. E. Goggin, A. Fedrizzi, T. Linjordet, M. Ams, G. D. Marshall, J. Twamley, M. J. Withford, et al., New J. Phys. 13, 075003 (2011).
- G. Corrielli, A. Crespi, G. Della Valle, S. Longhi, and R. Osellame, Nature Commun. 4, 1555 (2013).
- K. Poulios, R. Keil, D. Fry, J. D. A. Meinecke, J. C. F. Matthews, A. Politi, M. Lobino, M. Gräfe, M. Heinrich, S. Nolte, et al., Phys. Rev. Lett. 112, 143604 (2014).
- A. Crespi, L. Sansoni, G. Della Valle, A. Ciamei, R. Ramponi, F. Sciarrino, P. Mataloni, S. Longhi, and R. Osellame, Phys. Rev. Lett. 114, 090201 (2015).
- N. Lambert, Y.-N. Chen, Y.-C. Cheng, C.-M. Li, G.-Y. Chen, and F. Nori, Nature Phys. 9, 10 (2013).
- I. Kassal and A. Aspuru-Guzik, New J. Phys. 14, 053041 (2012).
- P. W. Anderson, Phys. Rev. 109, 1492 (1958).
- Y. Lahini, A. Avidan, F. Pozzi, M. Sorel, R. Morandotti, D. N. Christodoulides, and Y. Silberberg, Phys. Rev. Lett. 100, 013906 (2008).
- M. Segev, Y. Silberberg, and D. N. Christodoulides, Nature Photon. 7, 197 (2013).
- M. Mohseni, A. Shabani, S. Lloyd, and H. Rabitz, J. Chem. Phys. 140, 035102 (2014).
- X. Chen and R. J. Silbey, J. Phys. Chem. B 115, 5499 (2011).
- H. B. Perets, Y. Lahini, F. Pozzi, M. Sorel, R. Morandotti, and Y. Silberberg, Phys. Rev. Lett. 100, 170506 (2008).
- S. Longhi, Laser Photon. Rev. 3, 243 (2009).
- L. Levi, Y. Krivolapov, S. Fishman, and M. Segev, Nature Phys. 8, 912 (2012).
- F. Caruso, A. Crespi, A. G. Ciriolo, F. Sciarrino, and R. Osellame (2015), arXiv:1501.06438.
- P. G. Kwiat, A. J. Berglund, J. B. Altepeter, and A. G. White, Science 290, 498 (2000).
- J.-S. Xu, X.-Y. Xu, C.-F. Li, C.-J. Zhang, X.-B. Zou, and G.-C. Guo, Nature Commun. 1, 7 (2010).
- S. Longhi, Phys. Rev. Lett. 97, 110402 (2006).
- P. Biagioni, G. D. Valle, M. Ornigotti, M. Finazzi, L. Duò, P. Laporta, and S. Longhi, Opt. Express 16, 3762 (2008).
- M. Delanty, S. Rebić, and J. Twamley, Eur. Phys. J. D 66, 1 (2012).