arxiv: v1 [quant-ph] 21 Aug 2012

Transcription

1 Non-classical interference in integrated 3D multiports arxiv: v1 [quant-ph] 21 Aug 212 Thomas Meany, 1 Michael Delanty, 1 Simon Gross, 1 Graham D. Marshall, M. J. Steel, 1 and Michael J. Withford 1 1 Centre for Ultra-high bandwidth Devices for Optical Systems (CUDOS), MQ Photonics Research Centre, Department of Physics and Astronomy, Macquarie University, North Ryde, 219 NSW, Australia 2 Centre for Quantum Photonics, H. H. Wills Physics Laboratory and Department of Electrical and Electronic Engineering, University of Bristol, Merchant Venturers Building, Woodland Road, Bristol, BS8 1UB, UK thomas.meany@mq.edu.au Abstract: We demonstrate three and four input multiports in a three dimensional glass platform, fabricated using the femtosecond laser directwrite technique. Hong-Ou-Mandel (HOM) interference is observed and a full quantum characterisation is performed, obtaining two photon correlation matrices for all combinations of input and output ports. For the three-port case, the quantum visibilities are accurately predicted solely from measurement of the classical coupling ratios. OCIS codes: ( ) Quantum information and processing; (25.53) Photonic integrated circuits; ( ) Glass waveguides; (14.339) Laser materials processing; (23.737) Waveguides. References and links 1. M. Reck, A. Zeilinger, H. J. Bernstein, and P. Bertani, Experimental realization of any discrete unitary operator, Phys. Rev. Lett. 73, (1994). 2. G. J. Pryde and A. G. White, Creation of maximally entangled photon-number states using optical fiber multiports, Phys. Rev. A 68, (23). 3. M. Żukowski, A. Zeilinger, and M. A. Horne, Realizable higher-dimensional two-particle entanglements via multiport beam splitters, Phys. Rev. A 55, (1997). 4. A. Perez-Leija, R. Keil, H. Moya-Cessa, A. Szameit, and D. N. Christodoulides, Quantum state transformations by multiport array beam splitters, in CLEO: QELS-Fundamental Science, (Optical Society of America, 212), p. QTh1E T. Brougham, V. Košták, I. Jex, E. Andersson, and T. Kiss, Entanglement preparation using symmetric multiports, The European Physical Journal D - Atomic, Molecular, Optical and Plasma Physics 61, (211). 6. K. Mattle, M. Michler, H. Weinfurter, A. Zeilinger, and M. Żukowski, Nonclassical statistics at multiport beamsplitters, Applied Physics B, Lasers and Optics 6, (1995). 7. G. Weihs, M. Reck, H. Weinfurter, and A. Zeilinger, Two-photon interference in optical fiber multiports, Phys. Rev. A 54, (1996). 8. A. Politi, M. J. Cryan, J. G. Rarity, S. Yu, and J. L. O Brien, Silica-on-silicon waveguide quantum circuits, Science 32, (28). 9. Y. Zhang, L. McKnight, E. Engin, I. M. Watson, M. J. Cryan, E. Gu, M. G. Thompson, S. Calvez, J. L. O Brien, and M. D. Dawson, GaN directional couplers for integrated quantum photonics, Applied Physics Letters 99, (211). 1. H. B. Perets, Y. Lahini, F. Pozzi, M. Sorel, R. Morandotti, and Y. Silberberg, Realization of quantum walks with negligible decoherence in waveguide lattices, Phys. Rev. Lett. 1, 1756 (28). 11. Y. Bromberg, Y. Lahini, R. Morandotti, and Y. Silberberg, Quantum and classical correlations in waveguide lattices, Phys. Rev. Lett. 12, (29). 12. A. Peruzzo, M. Lobino, J. C. F. Matthews, N. Matsuda, A. Politi, K. Poulios, X.-Q. Zhou, Y. Lahini, N. Ismail, K. Wörhoff, Y. Bromberg, Y. Silberberg, M. G. Thompson, and J. L. OBrien, Quantum walks of correlated photons, Science 329, (21).

2 13. A. Martin, A. Issautier, H. Herrmann, W. Sohler, D. B. Ostrowsky, O. Alibart, and S. Tanzilli, A polarization entangled photon-pair source based on a type-ii ppln waveguide emitting at a telecom wavelength, New Journal of Physics 12, 135 (21). 14. S. Tanzilli, A. Martin, F. Kaiser, M. De Micheli, O. Alibart, and D. Ostrowsky, On the genesis and evolution of integrated quantum optics, Laser and Photonics Reviews 6, (212). 15. P. J. Shadbolt, M. R. Verde, A. Peruzzo, A. Politi, A. Laing, M. Lobino, J. C. F. Matthews, M. G. Thompson, and J. L. O Brien, Generating, manipulating and measuring entanglement and mixture with a reconfigurable photonic circuit, Nature Photonics 6, (212). 16. D. Bonneau, M. Lobino, P. Jiang, C. M. Natarajan, M. G. Tanner, R. H. Hadfield, S. N. Dorenbos, V. Zwiller, M. G. Thompson, and J. L. O Brien, Fast path and polarization manipulation of telecom wavelength single photons in lithium niobate waveguide devices, Phys. Rev. Lett. 18, 5361 (212). 17. D. Bonneau, E. Engin, K. Ohira, N. Suzuki, H. Yoshida, N. Iizuka, M. Ezaki, C. M. Natarajan, M. G. Tanner, R. H. Hadfield, S. N. Dorenbos, V. Zwiller, J. L. O Brien, and M. G. Thompson, Quantum interference and manipulation of entanglement in silicon wire waveguide quantum circuits, New Journal of Physics 14, 453 (212). 18. T. Gerrits, N. Thomas-Peter, J. C. Gates, A. E. Lita, B. J. Metcalf, B. Calkins, N. A. Tomlin, A. E. Fox, A. L. Linares, J. B. Spring, N. K. Langford, R. P. Mirin, P. G. R. Smith, I. A. Walmsley, and S. W. Nam, Onchip, photon-number-resolving, telecommunication-band detectors for scalable photonic information processing, Phys. Rev. A 84, 631 (211). 19. A. S. Solntsev, A. A. Sukhorukov, D. N. Neshev, and Y. S. Kivshar, Spontaneous parametric down-conversion and quantum walks in arrays of quadratic nonlinear waveguides, Phys. Rev. Lett. 18, 2361 (212). 2. A. Peruzzo, A. Laing, A. Politi, T. Rudolph, and J. L. O Brien, Multimode quantum interference of photons in multiport integrated devices, Nature Communications 2, 224 (211). 21. A. Aspuru-Guzik and P. Walther, Photonic quantum simulators, Nature Photonics 8, (212). 22. M. Delanty, S. Rebić, and J. Twamley, Novel collective effects in integrated photonics, The European Physical Journal D - Atomic, Molecular, Optical and Plasma Physics 66, 1 7 (212). 23. G. D. Marshall, A. Politi, J. C. F. Matthews, P. Dekker, M. Ams, M. J. Withford, and J. L. O Brien, Laser written waveguide photonic quantum circuits, Opt. Express 17, (29). 24. 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, and A. G. White, Two-photon quantum walks in an elliptical direct-write waveguide array, New J. Phys. 13, 753 (211). 25. L. Sansoni, F. Sciarrino, G. Vallone, P. Mataloni, A. Crespi, R. Ramponi, and R. Osellame, Two-particle bosonicfermionic quantum walk via integrated photonics, Phys. Rev. Lett. 18, 152 (212). 26. R. Keil, A. Szameit, F. Dreisow, M. Heinrich, S. Nolte, and A. Tünnermann, Photon correlations in twodimensional waveguide arrays and their classical estimate, Phys. Rev. A 81, (21). 27. A. W. Snyder, Coupled-mode theory for optical fibers, J. Opt. Soc. Am. 62, (1972). 28. S. Eaton, H. Zhang, P. Herman, F. Yoshino, L. Shah, J. Bovatsek, and A. Arai, Heat accumulation effects in femtosecond laser-written waveguides with variable repetition rate, Opt. Express 13, (25). 29. A. Szameit, F. Dreisow, T. Pertsch, S. Nolte, and A. Tünnermann, Control of directional evanescent coupling in fs laser written waveguides, Opt. Express 15, (27). 3. I. J. Myung, Tutorial on maximum likelihood estimation, Journal of Mathematical Psychology 47, 9 1 (23). 31. C. K. Hong, Z. Y. Ou, and L. Mandel, Measurement of subpicosecond time intervals between two photons by interference, Phys. Rev. Lett. 59, (1987). 32. D. F. V. James, P. G. Kwiat, W. J. Munro, and A. G. White, Measurement of qubits, Phys. Rev. A 64, (21). 1. Introduction Optical multiports provide the ability to produce an arbitrary unitary transformation of a set of optical modes [1]. Multiports are therefore an extremely important component in optical quantum information processing (QIP) and state preparation [2, 3], where producing entangled states is a key motivator for such devices [4, 5]. Multiports and their applications in quantum optics were originally explored using bulk optic beamsplitters and phase shifters [6], for which non-classical interference was observed in both three- and four-port splitters. Following from this work a study of a three-port fiber coupler was undertaken [7]. However neither of these approaches provide the scalability and stability requirements needed for large scale optical quantum circuits. The extreme difficulty of working with large numbers bulk optical elements and fibre components has limited the development of these devices.

3 The challenges of scaling and stability in various quantum optical applications have been the main drivers for the recent explosion of interest in integrated quantum devices, such as integrated waveguide circuits [8, 9], waveguide array quantum walks [1, 11, 12] and single photon sources [13, 14]. Desirable qualities such as reconfigurability have been demonstrated [15], notably in lithium niobate-based devices which allow for telecommunications wavelength photon manipulation at GHz speeds [16]. Miniaturisation using silicon-on-insulator nanowire based photonic circuits have been shown [17] and there is also great promise to incorporate these platforms with existing on chip detectors [18]. In addition the integration of photon pair sources and circuits have been achieved using a nonlinear waveguide array fabricated in lithium niobate [19]. The advantages of integrated components have been applied to the construction of integrated quantum multiports. Peruzzo et al. have demonstrated multiport functionality with an on-chip 4-port multimode interference device [2] in the silica-on-silicon platform. The majority of these demonstrations have relied on mature planar platforms, but there are limitations associated with the inherently two dimensional (2D) nature of planar quantum circuits. For instance in the realm of optical quantum simulation [21], there are interesting Hamiltonians which can only be effectively simulated using waveguides in three dimensional (3D) arrangements [22]. Moreover, as circuits become more complicated, both the circuit area and the proportion of neutral elements, whose only purpose is to provide transparent waveguide crossings, will become significant. Hence techniques which allow the fabrication of structures in three dimensions have some specific and unique possibilities. Such a technique is the femtosecond laser direct-write technique (FLDW) which allows for fabrication of waveguides in glass subsrates in 3D. This technique has been used to produce a number of devices showing quantum operation including directional couplers [23] (evanescently-coupled waveguides which act as beamsplitters), a 6-port continuous quantum walk on an elliptical array [24], a discrete quantum walk on a line using entangled photons [25] and investigations of quantum interference in 2D arrays [26]. This has demonstrated the potential of the FLDW technique to fabricate useful structures for quantum applications. In this paper we implement a 3-port and 4-port device using a 3D FLDW technique. Each device highlights a different advantage of the 3D geometry. The 3-port is implemented as a single coupling region that replaces three 2-port devices [1]. The 4-port is composed of four 2-port directional couplers, but avoids a crossing that would require an additional coupler in 2D. Taken together these kinds of structures could reduce the complexity and sensitivity of higher order concatenated multiport based circuits. This fact is established by our observation of 2 photon interference in both structures. 2. Implementing multiports The fundamental idealised multiport device is the familiar beamsplitter which implements a desired 2 2 unitary operation. This can be generalised to a multiport operation on N input modes â N into N output modes ˆb N described by the unitary N N matrix [6] ˆb i U 1,1 U U 1,N â i ˆb j. = U 2,1 U 2,2 U 2,N â j (1) ˆb N U N,1 U N,2 U N,N â N The integrated photonics analog of the basic 2-port beamsplitter is the directional coupler, in which two waveguides are brought into close proximity so that their evanescent fields can overlap and exchange energy [8, 27]. The degree of coupling is determined by the waveguide separation and interaction length which can be adjusted to produce directional couplers of any

4 reflectivity. Using the well-known decomposition of the N N unitary into a concatenation of beamsplitters and phase shifters [1], it is possible to construct an arbitrary multiport in two dimensions in integrated optics by designing a suitable arrangement of directional couplers and phase delays. However as the number of directional couplers scales with N 2, at large N, such a device is prohibitive to build due to fabrication errors. Previous investigations have looked at cases of bulk optic beamsplitters [6], fibre optic directional couplers [7], and multimode interference devices [2], to produce multiports In the next section we apply these techniques to the specific cases of a three and four port beamsplitter based on an arrangement of coupled waveguides 2.1. Mode evolution in the three port (tritter) In three dimensions a 3 3 multiport or tritter can be constructed from a single coupling element involving three waveguides (see Fig. 1a). This single element can replace up to three directional couplers and phase shifters [1]. The coupling coefficient G(z) describes coupling between waveguides 1 and 3, and g(z) describes the coupling between waveguides 1 and 2 and waveguides 2 and 3. In general, all three of the couplings could have different values, but, as will be discussed later, the isosceles geometry considered here is appropriate as a result of the laser writing process we used. (a) (b) (c) Fig. 1: (a) Idealised coupling region of a tritter device. Writing considerations induce an effective isosceles geometry with two coupling strengths g(z) and G(z). (b) A schematic of a 3D implementation of a tritter, where three waveguides taper from a planar arrangement to a triangular interaction region and taper back to a planar array and (c) the layout of the four directional couplers and the waveguide avoided crossing in the 4-port device Such a structure is described by a Hamiltonian of the form where the coupling matrix is Ĥ = hv p 3 n,m=1 C n,m â nâ m, (2) C = β g G g β g. (3) G g β The propagation constant β = ω/v p of the three waveguide modes is assumed to be common (v p is the phase velocity and ω/2π the optical frequency). Explicitly including the common phase velocity v p in Ĥ means that operator evolution with distance z rather than time t follows naturally. After propagation through the interaction region, photons have encountered the mode transformation ˆb = Uâ (4)

5 where the input and output mode vectors are â = (â 1,â 2,â 3 ) and ˆb = (ˆb 1, ˆb 2, ˆb 3 ), and the transfer matrix U is the solution of Heisenberg s equation, dâ/dz = ic(z)â. The transfer matrix U completely determines the operation of the device, in particular the classical output intensities and quantum interference between photons in the device Mode evolution in the four-port In contrast to the single element tritter the four port (depicted in Fig. 1c) is made up of four directional couplers. Here we have used the 3D advantage to cross two waveguides without interaction, reducing the number of required directional couplers by one. As the device is made up of four discrete elements, the device can be modelled as a product of beam splitter operators. Assuming each directional coupler is identical, the operator corresponding to the four port transfer matrix is, Û = ˆB 2,4 (η) ˆB (η) ˆP 2 (φ) ˆP 3 (φ) ˆB 3,4 (η) ˆB (η) (5) where the directional couplers have the beam splitter action, [ ] ] [ ] ˆb i [âi η i 1 η ][âi = ˆB (η) = ˆb j â j i, (6) 1 η η â j and ˆP j (φ) describes the phase difference between the inner crossing arms and the outer straight-through arms (See Fig. 1c) where ˆP j (φ)â j = e iφ â j. Using these relations we find the transfer matrix corresponding to the four port operator (5) is η ie iφ (1 η)η i (1 η)η e iφ (η 1) U = i (1 η)η e iφ η η 1 ie iφ (1 η)η ie iφ (1 η)η η 1 e iφ η i (1 η)η e iφ (η 1) i (1 η)η ie. (7) iφ (1 η)η η 3. Fabrication methods and design The devices were fabricated using the femtosecond laser direct-write (FLDW) technique which employs a tightly focused femtosecond laser to inscribe a localized refractive index change in glass. By translating the glass sample in the (x, y, z) directions with respect to the incident laser, arbitrary 3D regions of net-positive refractive index change (which act as waveguides) can be produced inside the host glass. The laser used to fabricate the structures is a titanium sapphire oscillator (Femtolasers GmbH, FEMTOSOURCE XL 5, 8 nm centre wavelength, < 5 fs pulse duration) with 5.1 MHz repetition rate. A 1X oil immersion objective (Zeiss N-Achroplan, NA = 1.25) was used to focus the laser inside boro-aluminosilicate (Eagle 2). The sample translation was completed using Aerotech high precision motion control stages with 1 nm precision at writing speeds of 2 mm/min. The waveguides were written using pulse energies of 28 µj which created waveguides suitable for single mode operation at 8 nm. The combination of low pulse energy, high translation speed and high repetition rate created waveguides which were at the at the very onset of cumulative heating modification mechanism [28] (implying that pulses arrive before the characteristic thermal diffusion time of the material). This results in smooth refractive index change and short device production time, allowing the fabrication of complex photonic circuits. The refractive index change is ellipsoidal in nature and it has been shown that laser written waveguides display non-isotropic coupling due to waveguide asymmetry [29]. Since the tritter is a triangular arrangement of modes (as shown in Fig. 1a) there are two effective coupling ratios, g(z) G(z). In both cases the waveguides are initially separated, in a planar array, by 127 µm before tapering an interaction region,

6 where evanescent coupling occurs with typical spacings of 15 µm and subsequently tapering out to 127 µm planar spacing.this is best illustrated in Fig. 1, which shows the waveguides tapering towards interaction regions and returning to the original spacign at the output of the device. We set out to fabricate a tritter with equal splitting between output ports of 33/33/33 and a 4-port with equal splitting between output ports of 25/25/25/25. Hence we fabricated a range of devices with different length interaction regions targeting the ideal coupling ratios. We performed a classical characterisation to determine which structures best matched our intended constraints. 4. Classical characterisation In order to determine the coupling ratios in the tritter and splitting ratios for the four port a classical characterisation was performed. If classical light of intensity M is injected into port j, the output intensity N at port k can be modelled as, N k j = ε in j ε out k M U j,k 2, (8) where the input and output losses at port l are characterized by the pre-factors, < ε in/out l < 1. By forming the following ratios, F k,s j,r Nk j Ns r Nr k N s = U j,k 2 U r,s 2 j U r,k 2 U j,s 2, (9) it is possible to relate the measured classical intensities directly to the unknown tritter couplings, g(z) and G(z) or four-port splitting ratios and cancel the loss terms. For the tritter we fitted a model to the classical measurements by assuming uniform coupling over an interaction region of length L, i.e. g( z L)L/ν p = ḡ and G( z L)L/ν p = Ḡ with zero coupling elsewhere. Using maximum likelihood estimation [3], we varied ḡ and Ḡ in the transfer matrix U to find the best fit to the three experimentally determined fractions F and F. This involved minimizing the difference between the experimentally measured left hand side and theoretical right hand side of Eq. (9) for F, F and F, weighted by the F uncertainties in the fractions (found from the uncertainties of the measured intensities N k j ) [?]. This fit determined ḡ.81ν p /L and Ḡ.51ν p /L, indicating that waveguides one and two couple more strongly than waveguides one and three. Using this fit we obtained the matrix of classical output intensities U tritter 2 = ,. (1) The ideal symmetric tritter has U 2 = 1/3, therefore we predict our device has a slightly asymmetric power splitting [6]. Similar to the tritter, we can use a maximum likelihood approach to fit η, the beamsplitter reflectivity, from (7) using the following six fractions, F, F, F1,4 1,4, F, F2,4 2,4, and F3,4 3,4. This fit determined η =.377. Substituting this fit into Eq. (7) we find the matrix of classical output intensities for the four-port U 4port 2 = (11)

7 Our four-port device differs from the ideal symmetric four-port, U 2 =.25 and has a nonuniform power splitting [6].This is due to the fact that the value we measure for η =.377 is not the ideal value of η = Quantum characterisation We performed two photon characterisation on the devices using 84 nm photons produced using type-i spontaneous parametric down conversion (SPDC) source, as described in [23]. The setup used a Toptica 42 nm blue diode laser which was focused into a 1 mm thick BiBO crystal cut to give a 6 opening cone angle. Down converted photons were produced at 84 nm in an output cone and then coupled into polarisation maintaining single mode fibres. These were then butt-coupled to the device under test using 127 µm spaced V-groove arrays. The outputs were then monitored using silicon avalanche photodiode detectors (SAPD) from Perkin Elmer. The 2-fold coincidences across all the output modes were measured simultaneously (in a 5 ns window) using a 4 channel time tagging unit. The arrival time of the photons into the chips was varied as a means of continuously controlling the degree of indistinguishability between the two photons and thereby moving between classical propagation of light and quantum interference. We detected signatures of two photon interference in the devices through the Hong-Ou- Mandel (HOM) effect [31], whereby the probability of two photon coincidences are reduced (HOM dip) or enhanced (HOM peak) due to quantum interference in the device. The degree of quantum interference is quantified using the visibility, V k,l, defined for injecting a single photon into waveguides and detecting in waveguides k,l as, V k,l = Ck,l Qk,l C k,l where, the quantum and classical coincidence probabilities are, Q k,l =, (12) δ U i,k U j,l +U i,l U j,k 2 (13) C k,l = U i,k U j,l 2 + U i,l U j,k 2. (14) We see that the quantum two photon coincidence probability Q k,l is strongly dependent on the relative phases of the elements of U. When these phases destructively interfere, two photon coincidences are reduced compared to the classical case, leading to a large positive visibility (a HOM dip). Whereas when these phases constructively interfere, two photon coincidences are enhanced compared to the classical case, leading to a large negative visibility (a HOM peak) Visibilities in tritter As the visibilities are strongly dependent on the relative phases of the matrix elements of U, the tritter visibilities depend on the length and geometry of the coupling region as defined by the coupling functions g(z) and G(z) in Eq. (2). To illustrate the variety of tritter visibilities that can be observed in our 3D device, we plot the visibilities as a function of effective coupling ḡl in (Fig. 2) assuming uniform coupling over the interaction region. Figure 2a shows a symmetric tritter as gl is varied. Only at one coupling ( the symmetric 33/33/33 case) do we see equal visibilities of 5

8 Visibility V V Visibility (a) (b) Fig. 2: Tritter visibilities as a function of effective coupling, ḡl. We assume uniform coupling over an interaction region of length L, g( z L)L/ν p = ḡ and G( z L)L/ν p = Ḡ with zero coupling elsewhere. a) symmetric coupling region (Ḡ = ḡ). b) Isosceles coupling region (Ḡ =.6234ḡ). Vertical line corresponds to the fit to our experimental tritter with ḡl.81ν p V V V C o u n ts (1 s -1 ) 4 P o rt 1-2 P o rt 2-3 P o rt C o u n ts (1 s -1 ) (a) (b) C o u n ts (1 s -1 ) Visibility (c) (d) Fig. 3: Tritter HOM dips for (a) input ports 1 and 2, (b) for input ports 2 and 3 (c) for input ports 1 and 3. (d) Comparison of measured tritter visibilities (error bars) and theoretical predictions (columns) from the simple tritter model in section V V We measured nine HOM dips for the tritter (Fig. 3) by injecting two photons into the device at different delays and observing coincidences at the output ports. The measured visibilities (Fig. 3d) were found from the relative depth of each dip. We see a clear reduction in the visi-

9 bilities observed in the cases of V, this corresponds to Fig. 2b, we also see reductions in the visibilities for different input ports and V. In section 4 we fitted a simple model of the tritter couplings g(z) and G(z) using only classical intensity measurements of the device. As this fit uniquely determines U we can predict the visibilities for this model using Eq. (12). Therefore we are able to predict, using only classical intensity measurements measurements, the quantum correlations for the device. This is a convenient tool, since typically in order to produce a device with the correct coupling ratios a range of device must be fabricated around this indended parameter. It is much more convenient to perform a range of intensity measurements prior to any quantum measurements to confirm the device has the desired functioning. The predicted visibilities are compared to the experimental visibilities in Fig. 3d, where we see a qualitative match, issues which are causing the minor mismatch between this prediction and the experimentally observed values are the assumptions uppon which the intensity dependent measurement is based. The model assumes a constant input intensity which can fluctuate during the course of the measurement and input output facet losses which may change also. Furthermore the model assumes the interaction region is constant in z but input and out tapering regions do contribute to the overall coupling which was not accounted for. In addition since the maximum likelihood measurment is a numerical simulation it has it s own inherent errors [3, 32] port We also measured HOM interference for all input and output combinations in the 4-port device. These measurements are shown in Fig. 4a 4f. We see a variety of different magnitudes and signs of the visibilities for the various input/output combinations. This is due to the quantum term Q k,l in Eq. (12), which exhibits constructive or destructive interference, depending on the relative phases of the elements of U. The fact that all outputs had a visibility greater than 5% (or less than ( 5%)) indicates that the output states are strongly non-classical. We now compare the experimental visibilities to those predicted from the 4-port model in section 4. In contrast to the tritter the classical intensity measurements in the 4-port are insufficient to determine the phase φ due to the additional path lengths of waveguides two and three. Therefore after fitting η =.377 from the classical intensity measurements we used visibility measurements from a single input, namely inputs two and three, to fit φ. Using maximum likelihood estimation we fitted φ from the six HOM dips in Fig. 4d. This fit determined φ =.7π. After substituting the determined η and φ into the transfer matrix Eq. (7) we can use (12) to predict the other 3 visibilities in Fig. 4. The predicted visibilities are compared to the measured visibilities (see Fig. 4g) and display a qualitative match, correctly identifying the sign and magnitude of each visibility. We note that using other input combinations to determine φ produced similar predictions.

10 C o u n ts (3 s ) P o r t 1-2 P o r t 1-3 P o r t 1-4 P o r t 2-3 P o r t 2-4 P o r t 3-4 C o u n ts (3 s ) (a) (b) 4 C o u n ts (3 s ) -1 C o u n ts (3 s ) (c) (d) 2 5 C o u n ts (3 s ) C o u n ts (3 s ) (e) (f) 1. 1,4 2,4.5 Visibility ,4 V V 1,4 V V 2,4 V 3,4 (g) Fig. 4: Four-port HOM dips for (a) input ports 1 and 2, (b) for input ports 1 and 3 (c) for input ports 1 and 4 (c) for input ports 2 and 3 (d) for input ports 2 and 4 (e) for input ports 3 and 4. (f) Experimental four-port visibility matrix determined from the HOM dips. (g) Comparison of measured four port visibilities (error bars) and theoretical predictions (columns) from the

11 6. Conclusion Two photon correlation measurements were performed on integrated 3 and 4 port devices. The visibilities observed for our tritter compare extremely well with bulk optic demonstrations where visibilities, for only one pair of input ports, V 5%,V 48% and V 26% are observed [6]. In addition our tritter displays greater symmetry than a fibre optic based example where visibilities, observed only for a single pair of input ports V 5%,V 3% and V 25% [7]. Our four port has shown a uniform output very similar to a bulk optic example, where visibilities of V 1,4 2,4 75%,V 66%,V 66%,V 8%,V 75% and V 3,4 66% were observed [6] which out performs a multimode interference device described in [2]. This demonstrates the applicability of laser written photonic circuits for some quantum logic applications where a 3D capability is advantageous for circuit size or complexity requirements. Acknowledgments This research was conducted in part by the Australian Research Council Centre of Excellence for Ultrahigh bandwidth Devices for Optical Systems (project number CE11118). This work was performed in part at the Macquarie University node of the Australian National Fabrication Facility. A company established under the National Collaborative Research Infrastructure Strategy to provide nano and microfabrication facilities for Australia s researchers. GDM acknowledges support from the Australian Academy of Science s International Science Linkages scheme and the Marie Curie International Incoming Fellowship scheme. We acknowledge the generosity of Prof. Gabriel Molina-Terriza and Dr Mathieu Juan in lending equipment and Dr. Douglas J. Little for offering valuable insights and observations.