Spatiospectral characteristics of parametric downconversion in waveguide arrays
Abstract
High dimensional quantum states are of fundamental interest for quantum information processing. They give access to large Hilbert spaces and, in turn, enable the encoding of quantum information on multiple modes. One method to create such quantum states is parametric downconversion (PDC) in waveguide arrays (WGAs) which allows for the creation of highly entangled photon pairs in controlled, easily accessible spatial modes, with unique spectral properties.
In this paper we examine both theoretically and experimentally the PDC process in a lithium niobate WGA. We measure the spatial and spectral properties of the emitted photon pairs, revealing correlations between spectral and spatial degrees of freedom of the created photons. Our measurements show that, in contrast to prior theoretical approaches, spectrally dependent coupling effects have to be taken into account in the theory of PDC in WGAs. To interpret the results, we developed a theoretical model specifically taking into account spectrally dependent coupling effects, which further enables us to explore the capabilities and limitations for engineering the spatial correlations of the generated quantum states.
pacs:
03.65.w, 42.50.p, 42.50.Dv, 42.65.Lm, 42.65.Wi1 Introduction
The process of parametric downconversion (PDC) [1] provides a versatile resource for the creation and engineering of sophisticated quantum states of light. In particular, the creation of photonpair states, entangled in various degrees of freedom, is of special interest for quantum information and quantum communication applications. To date several experiments have investigated the entanglement of these photon pairs in various degrees of freedom, for example using their spectral properties [2, 3, 4, 5, 6], the orbital angular momentum of the photons [7, 8, 9], transverse modes of a single waveguide [10], polarization [11, 12] and energytime [13, 14].
Recently, PDC in a waveguide array (WGA) with a nonlinear optical response has been proposed as a resource for preparing highly entangled photon pairs into controllable and easily accessible spatial modes [15, 16, 17]. In a nonlinear WGA several waveguides are embedded into a nonlinear material, such as lithium niobate [18], and are grouped so close to each other, that the created PDC photons are able to couple evanescently to the nearestneighbour waveguides. This introduces, aside from the spectral degree of freedom present in the PDC process [19, 20, 21], a new spatial degree of freedom, where each waveguide channel represents an individual, easily accessible spatial mode. In general, the use of the spatial degree of freedom is a wellknown tool for state engineering in bulk crystals [22]. PDC in WGAs, however, benefits from the bright and efficient photonpair production in waveguides [23, 24], and offers several further advantages: these systems are not subjected to any incoupling losses and thus maintain the process efficiency and quantum features. Further, a miniaturized, onchip realization is also highly beneficial concerning the scalability, stability, adjustability and coherence of the system. Moreover, active optical devices can be implemented in the nonlinear substrate of lithium niobate, thus allowing the implementation of flexible, fast switching integrated optical components on the chip for controlling the properties of quantum light [25].
The linear and nonlinear characteristics of WGAs have attracted enormous interest for the investigation of classical wave phenomena [26]. In the context of quantum optics these systems enable the implementation of quantum walks [27, 28, 29]. In combination with nonclassical light sources, WGAs have further been utilized for the study of bosonic and fermionic behavior [30, 31, 32]. WGAs have also become a resource for simulating the properties of other quantum systems. Utilizing classical light sources, WGAs with different types of architectures have been used to mimic the properties of photonpair states [33], displaced Fock states [34], and faithful quantum state transfer [35].
Consequently, integrating WGAs in nonlinear materials offers a variety of possibilities for quantum optical experiments and a high potential for integrated quantum state generation and manipulation. For this purpose a detailed and accurate understanding of the process properties is necessary to enable precise quantum state engineering.
In this paper, we theoretically as well as experimentally, investigate PDC in a nonlinear WGA with special attention to the correlations between the spatial and spectral degrees of freedom. We explicitly take into account frequency dependent coupling effects and show how the frequencies of the created photon pairs are connected with their spatial properties. With this knowledge we then explore the preparation of specific spatial correlations between the created photon pairs.
The paper is structured into three main parts. In section 2, we introduce our theoretical PDC model including spectrally dependent coupling effects. In section 3 we investigate the occurring spatiospectral correlations and introduce different ways to utilize them in order to create highly sophisticated quantum states. Finally, in section 4, we perform a spectrally and spatially resolved measurement of the PDC emission in a periodically poled lithium niobate (PPLN) WGA. Our measurement results demonstrate the predicted spatiospectral correlations and confirm the necessity to include spectrally dependent coupling effects into the theoretical treatment of PDC in a WGA. Furthermore our experimental investigations enable us to extract the phasematching curve of the WGA from the acquired data, which is in very good agreement with the theoretical prediction.
2 PDC in nonlinear WGAs
During the PDC process a pump photon is converted inside a nonlinear crystal into a photon pair usually labelled signal and idler. In this paper we consider PDC in a WGA, where the generated photons can couple from waveguide to waveguide before they exit the crystal, while the pump is constrained to one waveguide, as schematically depicted in figure 1.
2.1 Electric fields in WGAs
Following the theoretical treatment of typeI PDC in nonlinear WGAs by Solntsev et al. in [15] we start with the mathematical description of the quantized electric fields propagating through the WGA and use coupled mode theory to solve the corresponding Maxwell’s equations [36, 37, 38, 39]. The nearest neighbour coupling in the WGA induces a modified, discretized dispersion [26] altering the propagation vector in the direction
(1) 
where describes the propagation vector in a single waveguide at frequency , given by , and is its normalized transverse momentum, with as the distance between the waveguides
(2) 
where collects all constants [41]. The pump field driving the PDC process is considered as a bright, undepleted optical beam, which allows us to treat it classically. Moreover, the pump does not couple to neighbouring channels, because it resides at wavelengths far below signal and idler. In this regime the waveguide mode size is smaller, and thus, the mode overlap governing the coupling parameter can be approximated as zero. Therefore the influence of the WGA on the propagation vector can be neglected (see A). Thus, the pump () field in the th waveguide is described as
(3)  
where is the spectral shape of the pump beam. The spatial illumination pattern of the pump in the waveguide is connected to its Bloch mode distribution via a Fourier transformation .
2.2 The PDC state
Using the electric field definitions in (2) and (3), we are able to express the effective Hamiltonian of the PDC process in a nonlinear WGA by
(4) 
where denotes the electric constant, is the nonlinearity of the material and the length of the WGA with an infinite number of channels .
The nonlinear interaction inside the medium is weak, which enables us to calculate the output state using firstorder perturbation theory
(5) 
We have extended the integration boundaries to plus and minus infinity, since we regard the PDC state long after the interaction in the crystal [42]. We further postselect on the detection of photon pairs, which enables us to drop the vacuum contribution appearing in (5) and renormalize the state accordingly. A straightforward calculation delivers the expression for the twophoton PDC state emerging from the WGA as
(6) 
i.e., two photons are created into a superposition of spectral and spatial modes, where is the normalization constant. The spatial and spectral structure of the created photon pair is determined by the form of the joint spatiospectral amplitude of the generated signal () and idler () photons. It assumes the form
(7)  
in which the phasemismatch is defined as
(8) 
The four terms determining the properties of the twophoton state in (7) can be attributed to the following sources. The spectral and spatial properties of the pump field determine the functions and respectively. The phasematching function and the phase factor are given by the nonlinear and linear properties of the photonic structure, i.e. the dispersion in the medium, periodic poling and waveguide spacing.
3 Spatiospectral properties of PDC in WGAs
3.1 Origin of the spatiospectral correlations between the created photon pairs
The key for understanding the spatiospectral properties of the generated PDC state lies in the phasematching function . Significant downconversion rates require nearly perfect phasematching, satisfying . This constraint induces spatiospectral correlations between the frequencies and components of the generated signalidler photon pairs.
In order to understand its exact nature we split the phasemismatch in equation (8), into its contributions from the single, isolated waveguide and the dispersion effects from the WGA
(9) 
The first term describes the spectral phasemismatch, known from PDC in a single isolated waveguide,
(10) 
which stems mostly from the dispersion of the material with a contribution from the modal properties of the waveguides [21]. The second term is the phasemismatch induced by the dispersion of the WGA,
(11) 
which depends on the transverse components of the generated photons and the frequency dependent coupling parameter .
In figure 2(a) we sketch the contours of the spectral phasemismatch defined in (10) for various values. Due to dispersion effects of the material, this yields curved lines aligned along the axis, for constant in space. In the case of negative spectral phasemismatch () the contours (red lines) for different values of run below the central contour, colourcoded in black, (), while for positive spectral phasemismatch () the contours run above (green lines).
In order to achieve phasematching in the WGA, i.e , the spectral phasemismatch has to be compensated by a spatial phasemismatch of with opposite sign to satisfy . The pairing of these values, and creates the spatiospectral correlations between the PDC photons. The different contours for which the different values of and yield perfect phasematching are colourcoded in figure 2(a)(d).
To investigate the induced correlations we choose three different scenarios for : , and , as shown in figures 2(b)(d).
In figure 2(b) we show the combinations yielding a constant at the degeneracy point (). A numerical simulation illustrating this scenario can be found in C. In this regime both coupling parameters are identical and, in the considered frequency range, constant: . In order to understand the depicted shapes we regard three different cases, where (i) , (ii) and (iii) :
 Scenario (i):

If the contributions from and in (11) cancel each other. This means that if the transverse momentum of the signal photon is in the center of the Brillouin zone
^{2} () the momentum of the idler has to be at the edge () or vice versa. This results in a rectangular shape in ()space, depicted in figure 2(b).  Scenario (ii):

If , the correlations between and form a circlelike shape in the ()space, as depicted by the green line in figure 2(b). If either signal or idler is at the center of the Brillouin zone ( or ) than the other one has to have a component with a value of in order to reach the required phasemismatch.
 Scenario (iii):

If , the solutions are located at the corners of the ()space, as visualized by the red line in figure 2(b). In order to reach the required negative phasemismatch the transverse momentum of scenario (ii) has to be shifted from the center to the edges of the Brillouin zone.
In figure 2(c) we depict the case in a frequency range where both coupling values are constant (). The different coupling values lead to stretched contours, in comparison to figure 2(b). Finally, in figure 2(d) we depict , which is a flipped version of figure 2(c), due to the fact that we regard typeI PDC, where signal and idler are interchangeable. Here the different coupling constants are selected such, that the difference is strong enough to impact the spatial correlations.
In conclusion our investigations reveal that the origins of the spatiospectral correlations of PDC in WGA, are located in the phasematching condition . Every spatial correlation pattern, as depicted in figure 2(b)(d) has a fixed spatial phasemismatch , which must be compensated by an appropriate spectral phasemismatch , satisfying . Next, we will illustrate, how to utilize this spatiospectral correlation in order to generate sophisticated spatial correlations between the generated signal and idler photons.
3.2 Engineering spatial correlations via pump shaping
In section 3.1 we have shown that the spatial and spectral correlations in the phasemismatch have a significant influence on the emitted photon pairs. Driving the PDC process in the WGA with an adapted spatiospectral pump shape enables us to engineer a variety of different spatial correlations between the generated photon pairs. The emerging joint spatiospectral amplitude is given by (7)
(12) 
This means we can regard as a function which, via multiplication, selects a subset of the spatiospectral correlations inherent in the WGA. The phase is neglected in the discussion, since it does not influence the space correlations, while only refining the correlation patterns in real space.
For degenerate signal and idler frequencies, the whole process is sketched, from top to bottom, in figure 3. This corresponds to the configuration presented in figure 2(b). In this special case the spectral phasemismatch is aligned along the axis in frequency space (dashed black lines in figure 3 (a)(a)).
Choosing a specific pump frequency () defines a line in frequency space, aligned along the axis (torquoise shaded area in figure 3 (a)(a)) and consequently selects a specific phasemismatch
We can gain an enhanced control by adapting the spatial illumination of the WGA to modify , i.e. we change from pumping in a single channel () to pumping in different channels simultaneously. Here, we choose a rectangularly formed spatial pump shape in space, which offers the advantage that we can finely select the spatial correlations. The spatial illumination pattern of the WGA to reach this shape in space is highly intricate, such that an experimental realization becomes more involved.
Three different scenarios are colourcoded in figure 3(b)(b) by the red, green and blue shaded areas, where we have displaced the rectangular spatial pump shape in space. Each consists of a different range of components specifically chosen to excite distinct correlations in real space. This simultaneous pumping in different channels enables us to individually select various parts of the phasemismatch , as given by colourcoded areas, only constrained by the fact that the function has to be aligned along the axis in space. This enables us, for example, to engineer photon bunching and antibunching effects
3.3 Phase engineering of
In the previous section we have concentrated on the influence of real valued pump shapes on the correlation functions in both  and real space. However, we are also able to engineer the signalidler correlations by tuning the phase of the spatial pump distribution . For this discussion, we restrict ourselves to the scenario discussed in section 3.2 figure 3(a)(c) with a rectangular pump corresponding to an selecting the red shaded area.
The applied pump shapes are given in figures 4(a)(a), where red denotes the absolute value and green the phase of the function . We vary from a constant, over a linear to a quadratic phase in the spatial pump distribution. The exact parameters for the simulation are given in C.
Note, that the absolute value of the correlations , as depicted in figure 4(b)(b) remains unchanged, as there is no variation in the absolute value of the pump function. The phases of the  correlations however, vary with the induced phase, as depicted in figures 4(c)(c).
In the case of a constant, zerophase distribution in , the resulting real space correlations in figure 4(d) show the antibunching contribution, which has already been discussed in section 3.2.
Next, we introduce a linear phase in the pump shape, as depicted in figure 4(a). This is a trivial phase distribution, since a linear phase in space, can be realized by simply centering the spatial pump pattern around a waveguide different than zero. It leads to an asymmetry in the corresponding phase correlation as depicted in figure 4(c). Under Fourier transformation the phase asymmetry results, as expected, in a shift of the real space correlations, as given in figure 4(d), leaving the internal shape of the correlation function unchanged.
Quadratic phases, as shown in figure 4(a) introduce curvatures in the phases of the correlation (figure 4(c)). The result is a ”stretching” of the real space correlations in figure 4(d). This is due to the linear slope of the quadratic functions, leading to different displacements in the real space correlations.
In general the shapes shown in figures 3 and 4 highlight both the flexibility and the limitations of PDC in WGA structures. More complicated pump shapes, like sinusodial shapes in the absolute value of the pump functions or different phase distributions only lead to various superpositions, displacements or stretching of the already discussed real space correlations in section 3.2.
3.4 Enhanced engineering of spatial correlations via frequency filtering
We are able to further engineer the spatial correlations between the emitted photon pairs via frequency filtering [43]. In this approach we use two distinct rectangular filters, with upper and lower bounds , to select photon pairs within a narrow frequency range. This enables us to select a single point in the space depicted in figure 2(a) and consequently the corresponding contour in figure 2(b)(d). Postselecting on both photons passing their individual filters, this procedure modifies the output state to
(13) 
A particular example is visualized in figure 5, for nondegenerate photonpair filtering, which corresponds to figure 2(c). In this scenario we assume signal photons with a wavelength of about 1400 nm and idler photons at about 1600 nm. The exact simulation parameters are given in C, while the corresponding coupling parameters are depicted in figure 11. Figures 5(a)(a), show the spectral shape of , where the rectangle depicts the individual filter placement. Below, in figures 5(b)(b) we depict the selected correlations. Finally figures 5(c)(c) show the corresponding spatial correlations (as defined in B), obtained via Fourier transformation, where label the individual waveguide channels.
3.5 Summary: Spatiospectral correlations
In summary, we explored the spatiospectral structure of PDC in WGAs with special attention to spectrally dependent coupling effects. Our investigations revealed that the spatiospectral correlations are fundamentally determined by the phasemismatch . Furthermore the explicit frequency dependence of the coupling parameter leads to a spectral dependence of the emerging space correlations. Via pump shaping and adapted spectral filtering this system allows for the flexible preparation of photon pairs featuring a variety of spatial correlations. It is, however, fundamentally limited to correlations already present in the phasemismatch .
4 Experiment
In this section we present a first preliminary experiment to test our theory and investigate the PDC process in a WGA. We examine the spatial and spectral properties of the PDC emission to reveal spectrally dependent coupling effects.
In our experiment, schematically depicted in figure 6, we employed tunable picosecond Ti:Sapphire laser pulses (\unit76\mega\hertz repetition rate, \unit775\nano\meter central wavelength and \unit0.8\nano\meter bandwidth) as a pump for the PDC process in a \unit40\milli\meter long, typeI, PPLN WGA. The array includes 101 waveguides packed into a region with a total width of \unit1.6\milli\meter and it is held during the experiment in an oven at the temperature of . A small portion of the pump beam is directed to a bright light spectrometer in order to control the pump wavelength. After passing through the power and polarization control (not shown), the pump beam illuminates a single waveguide. The PDC light generated in the pumped waveguide then couples to the neighbouring channels. The light launched out of the WGA was sent via a periscope (not shown) to a spectral filter, which blocks the residual pump. Thereafter, the spatially spread PDC emission passed through a dispersive prism and was imaged with a single lens (not shown) on an InGaAsdetector array sensitive down to the fewphoton level. Note that by using a twodimensional detector array, this configuration allows us to simultaneously resolve the properties of impinging light both spectrally and spatially. Thus, we can analyze any spatiospectral correlations of the generated photons and compare it with our theoretical predictions. The applied detector has a nearly flat response in the near infrared wavelengths. However, its sensitivity drops fast in the vicinity of \unit1.7\micro\meter. The spectral resolution of our detector is approximately \unit10\nano\meter, which was calibrated with a tunable telecommunication continuouswave laser coupled through a single isolated test waveguide on the same chip. This finite resolution introduces a smoothing of the experimental data.
After setting the pump wavelength close to the degeneracy point, which was found via a measurement of the secondharmonic response in the WGA, we recorded images of the spatially and spectrally resolved PDC emission. Our results in figure 7(a)(d) show the measured signal and idler wavelengths in each waveguide channel, illustrating the spatiospectral structure of the PDC state.
Close to the degeneracy point (figure 7(a)), the spatial spread spans over five waveguide channels and the central channel has a spectral bandwidth of approximately \unit100\nano\meter. We then tuned the pump wavelength in small steps from the degeneracy point towards lower wavelengths. As shown in figures 7(b)(d), we observe that the PDC emission is separated into two spectral regions due to the curvature of the phasematching function. Additionally, one can recognize that the spatial spread of the PDC photons clearly behaves differently at the upper and lower spectral branches. The coupling of the PDC photons into the neighbouring channels increases with the growing wavelength and causes an obvious and detectable imbalance to the measured distributions.
In order to evaluate the quality of our theoretical model in section 2, we compare the experimental results with numerical simulations. Since our model explicitly takes into account the wavelengthdependent coupling, it explains the differences in the observed spatial spreads at different wavelengths. The parameters governing the wavelength dependent coupling (see A) are determined by comparing the spatial spread given by the theoretical model to the measured ones. Additionally, the order of magnitude of the coupling parameter was confirmed by illuminating a single waveguide with laser featuring a fine linewidth, tuning the light in the range \unit15201600\nano\meter and measuring the linear spread of the light in the WGA.
For the comparison between the experimental and the numerical results, we have adapted the wavelength dependent coupling parameter together with the degeneracy wavelength of the process to reproduce the measured spatial spreads and central wavelengths of the PDC photons as well as their spectral marginals. The obtained numerical simulations depicted in figures 7(e)(h) show a very good quantitative agreement with the experimental data. For the details of the simulation see C. Besides the wavelengthdependent spatial spread, reflecting clear correlations between the spatial and spectral degree of freedom, the theoretical model also predicts a spectral double peak structure for the nondegenerate PDC emission, which is clearly visible at the higher wavelength branch in figures 7(f)(h). This effect is caused by the splitting of the phasematching function into two regions due to the different coupling parameters for signal and idler outside the degeneracy. We, however, are not able to observe this effect in the experiment due to the limited spectral resolution of the applied spectrometer. For the same reason the measured spectral marginals are slightly broader than the ones predicted with our theoretical model.
Due to the high pump powers utilized in our experiment, we cannot eliminate the possibility of higher order photonnumber contributions. However, in our system, the generated photons are widely spread in the spatial and spectral domain, preventing stimulation or other high gain regime effects [44, 45]. In this regime, the generated higher order photon pairs are independent and possess the spatiospectral characteristics, as discussed in section 3. As such, our measurements are equivalent to measuring many photonpair states subsequently.
For a more direct test of our theoretical model with the experiment, we firstly compare in figure 8 the spatial marginal distributions from figures 7 (d) and (h) both at the upper and the lower spectral branches. Although some differences appear in the shape of the spatial spreads between the experiment and theory, our model reproduces their widths correctly. Further, the measured spatial distributions also illustrate an excellent symmetry with respect to the central waveguide channel, which indicates a homogenous WGA quality. Secondly, we deduce the phasematching curve from the maxima of the measured spectral distributions in figure 7 and compare the experimental curve to the one gained from our theoretical model. As seen in figure 9 our simulation accurately predicts the observed shape of the phasematching curve. Small discrepancies arise from the accuracy at which the degeneracy point can be determined, limited resolution of the homemade spectrometer and fast decrease of its sensitivity at higher wavelengths. The splitting of the signal and idler spectra into two nonoverlapping spectral regions caused by the curvature of the phasematching function is already visible when tuning the pump wavelength by only a fraction of a nanometer from the degeneracy.
In regard to the spatiospectral characteristics, we find a very good agreement between the theory for photon pairs and our experiment, even though our detector lacks true singlephoton resolution and several photon pairs can be created by the pump at high power levels. Our results show that a detailed model, taking into account wavelength dependent coupling effects, is necessary in order to understand the experimental results. With our results we are able to determine the relevant physical parameters required for measuring the desired correlation properties between the signal and idler photons.
5 Conclusion
In this paper we have theoretically and experimentally investigated the PDC in nonlinear WGAs, explicitly taking into account spectrally dependent coupling effects. Our analysis revealed that the spatial and spectral degrees of freedom of the PDC process in the WGA are connected via the phasemismatch . A phasemismatch introduced in the spectral domain is compensated via an opposite phasemismatch in the spatial domain , leading to spatiospectral correlations in the generated quantum states. Thus, being able to select only specific regions of the phasematching, by means of appropriate spectral filtering or by choosing the suitable spectral and spatial characteristics of the pump beam, we can modify the joint spatiospectral amplitude. Consequently, we are able to engineer the emerging spatial correlations of the twophoton state. This leads to a variety of well defined spatial correlation patterns between the signal and idler photons limited only by the dispersion relation of the transverse momentum.
Our experimental investigations show that a wavelength dependent coupling parameter, which quantifies the strength of the evanescent overlap within the WGA, must be taken into account to understand the properties of the generated twophoton states. With our model, which takes into account the wavelength dependency of the coupling constant, we were able to accurately predict the spatial spread of the PDC photons into the adjacent waveguide channels at different wavelength regimes. Furthermore, we found that the phasematching curve of the PDC process, extracted from the measured data is in very good agreement with the theoretical prediction.
In conclusion, our analysis provides a framework for a detailed understanding of the possibilities and limitations of quantum state generation in nonlinear WGAs and enables an accurate engineering of a variety of spatiospectral correlations. Our results are important for the advancement of future applications such as quantum walks or other sophisticated optical networks on WGAs.
6 Acknowledgments
The authors thank Harald Herrmann and Hubertus Suche for useful discussions and helpful comments. The authors further thank Frank Setzpfandt and Thomas Pertsch for discussions concerning the experimental setup.
A. G., C. S. H. and I. J. received funding from MSM 6840770039, RVO 68407700 and GACR 1333906S. A. G. acknowledges partial support from the National Research Fund of Hungary under contract No. T83858.
Appendix A Wavelength dependency of the coupling parameter
Only over a narrow frequency range the coupling parameter can be approximated as a constant [26, 15]. In this paper, however, we consider spectral widths spanning over 100 nm in range. We therefore include wavelength dependent coupling effects into our analysis [38]. Much work on the frequency dependence of in LiNbO has already been performed [46, 47], but it has mainly relied on numerical calculations. Here we consider a simple model for the wavelength dependence of the coupling parameter , as shown in figure 10.
In coupled mode theory, the coupling parameter is given by the overlap integral from the mode of waveguide (I) with the mode in the neighbouring waveguide (II) as described in [48]
(14) 
Here is a constant depending on the mode number and is the refractive index difference between the waveguide and the surrounding bulk material [47] and is zero outside waveguides (I) and (II). Assuming a very simple model for the waveguide modes, which only considers the exponential part of the evanescent field, we arrive at
(15)  
(16) 
Here are the amplitudes of the electric field, and is the damping factor determined by the difference between the propagation vector in the waveguide and in bulk material [49]. We model the wavelength dependence of this parameter via , where describes the refractive index in terms of wavelength . We further neglect the possible influence of the waveguide dopant on the wavelength dependent refractive index. Solving (14) with these assumptions and integrating over the region of waveguide (II) yields
(17) 
where is a constant depending on the width of the waveguides and their mode profiles, as well as on the refractive index difference between the waveguide and substrate and is a constant that is determined by the distance between the waveguides and the waveguide mode profiles.
Appendix B Spatiospectral correlation function
The most straightforward way to obtain information about the joint spatiospectral distribution of signal and idler are spectrally and/or spatially resolved correlation function measurements. The correlation function between one photon at and the other photon at from the twophoton PDC state in (6) is given as:
(18)  
By tracing over the spectral degree of freedom, the correlation function in space can be written in the form
(19) 
In the real space the correlation function is determined as
(20) 
in which , and the biphoton state in real space, is given by a twodimensional Fourier transformation of in space.
Appendix C Numerical simulation of the WGA
In the scope of this paper we consider extraordinarily polarized typeI PDC in lithium niobate using the dispersion relations from [50]. We run our experiments far above room temperatures. Consequently our theoretical model takes into account the temperature dependence of the refractive indices, also given in [50].
We adjust the signal and idler frequencies with the help of quasi phasematching [21], where a periodic sign change of the nonlinearity in the material introduces a modification in the spectral phasematching condition
(21) 
In the simulation is an effective grating period, fitting the theoretically predicted wavelengths of the phasematched PDC process to the experimental values.
The wavelength dependent coupling parameter is fitted to the experimental results using the description in A. Its exact form is depicted in figure 11.
The process properties for the numerical modelling are given in table 1, while the experimental parameters can be found in section 4.
Properties of WGA (numerical simulation)  

Waveguide length [\unit\metre]  0.04 
Total number of waveguides on the array  41 
Temperature [\unit\celsius]  185.0 
Effective coupling parameter  
Effective damping coefficient [m]  
Phase engineering (Figure 4)  
Phasematched wavelength pair (, ) [nm]  (1550.0,1550.0) 
Narrowband filtering (, ) [nm]  (1549.9, 1550.1) 
Real valued spatial pump width  
Number of waveguides in the array (only this section)  49 
Effective coupling parameter  
Effective damping coefficient [m]  
Correlation calculation (Figure 5)  
Phasematched wavelength pair (, ) [nm]  (1400.0,1600.0) 
Effective coupling parameter  
Effective damping coefficient [m]  
Experimental comparison (Figure 7(eh))  
Phasematched wavelength pair (, ) [nm]  (1549.8,1549.8) 
Spectral pump bandwidth [nm]  0.5 
Neardegenerate correlation calculation (Figure 12)  
Phasematched wavelength pair (, ) [nm]  (1550.1,1550.1) 
Spectral pump bandwidth [nm]  0.5 
Coupling parameter [1/m]  400 
Number of waveguides in the array (only this section)  51 
Figure 12 shows the output correlations available when pumping near the degeneracy wavelength, corresponding to figure 2(b) in section 3. Here, we exploit the tunability of the spatial output correlation patterns via the spectral pump shape, as discussed in section 3.1.
The selection process is shown from top to bottom. In figures 12(a)(a) we select the spatial correlations by tuning the pump frequency according to the selection rules discussed in section 3.1. The small inlays depict the position of the pump shape in relation to the phasematching function, serving as orientation in the spectral domain. Figures 12(b)(b) depict the selected spatial correlation in space in accordance to the spatiospectral interplay shown in figure 2(a) and (b). The corresponding real space correlation functions in figure 12(c)(c) are then determined via Fourier transformation.
References
Footnotes
 A more detailed introduction to transverse momenta of particles and bandstructures may be found in [26] or [40].
 The Brillouin zone is the elemental cell of a lattice in Fourier space. Due to its periodic nature all momenta can be expressed in this part of the Fourier space[40].
 In this configuration the impacts of the curvature of are negligible.
 This corresponds to pumping in a single waveguide.
 For higher spectral phasemismatch the shapes are of a circular form.
 Due to the Fourier limitation, the finite width of the spatial pump will lead to finite widths in the real space correlations as well.
References
 D. C. Burnham and D. L. Weinberg. Observation of simultaneity in parametric production of optical photon pairs. Phys. Rev. Lett., 25:84, 1970.
 T. E. Keller and M. H. Rubin. Theory of twophoton entanglement for spontaneous parametric downconversion driven by a narrow pump pulse. Phys. Rev. A, 56:1534, 1997.
 W. P. Grice, A. B. U’Ren, and I. A. Walmsley. Eliminating frequency and spacetime correlations in multiphoton states. Phys. Rev. A, 64:063815, 2001.
 M. V. Fedorov, M. A. Efremov, P. A. Volkov, E. V. Moreva, S. S. Straupe, and S. P. Kulik. Anisotropically and high entanglement of biphoton states generated in spontaneous parametric downconversion. Phys. Rev. Lett., 99:063901, 2007.
 G. Brida, V. Caricato, M. Genovese, M. Gramegna, M.V.Fedorov, and S.P.Kulik. Characterization of spectral entanglement of spontaneous parametricdown conversion biphotons. EPL, 87:64003, 2009.
 A. Christ, K. Laiho, A. Eckstein, K. N. Cassemiro, and C. Silberhorn. Probing multimode squeezing with correlation functions. New. J. Phys., 13:033027, 2010.
 A. Mair, A. Vaziri, G. Weihs, and A. Zeilinger. Entanglement of the orbital angular momentum states of photons. Nature, 412:313, 2001.
 G. MolinaTerriza, J. P. Torres, and L. Torner. Twisted photons. Nature Phys., 3:305, 2008.
 J. Leach, B. Jack, J. Romero, A. K. Jha, A. M. Yao, S. FrankeArnold, D. G. Ireland, R. W. Boyd, S. M. Barnett, and M. J. Padgett. Quantum correlations in optical angleorbital angular momentum variables. Science, 662:329, 2010.
 P. J. Mosley, A. Christ, A. Eckstein, and C. Silberhorn. Direct measurement of the spatialspectral structure of waveguided parametric downconversion. Phys. Rev. Lett., 103:233901, 2009.
 P. G. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, A. V. Sergienko, and Y. Shih. New highintensity source of polarizationentangled photon pairs. Phys. Rev. Lett., 75:4337, 1995.
 Z. H. Levine, J. Fan, J. Chen, and A. L. Migdall. Polarizationentangled photon pairs from a periodically poled crystalline waveguide. Opt. Express, 19:6724, 2011.
 J. G. Rarity, P. R. Tapster, E. Jakeman, T. Larchuk, R. A. Campos, M. C. Teich, and B. E. A. Saleh. Twophoton interference in a machzehnder interferometer. Physical Review Letters, 65(11):1348â1351, September 1990.
 Z. Y. Ou, X. Y. Zou, L. J. Wang, and L. Mandel. Observation of nonlocal interference in separated photon channels. Physical Review Letters, 65(3):321â324, July 1990.
 A. S. Solntsev, A. A. Sukhorukov, D. N. Neshev, and Y. S. Kivshar. Spontaneous parametric downconversion and quantum walks in arrays of quadratic nonlinear waveguides. Physical Review Letters, 108:023601, 2012.
 A. S. Solntsev, A. A. Sukhorukov, D. N. Neshev, and Y. S. Kivshar. Photonpair generation in arrays of cubic nonlinear waveguides. Optics Express, 20(24):27441â27446, 2012.
 A. S. Solntsev, A. Clark, F. Setzpfandt, M. J. Collins, C. Xiong, A. Wu, F. Eilenberger, A. Schreiber, F. Katzschmann, R. Schiek, W. Sohler, A. Mitchell, C. Silberhorn, B. J. Eggleton, T. Pertsch, A. Sukhorukov, D. Neshev, and Y. Kivshar. Simultaneous photonpair generation and quantum walks in a waveguide array. In Frontiers in Optics, 2012.
 R. Iwanow, R. Schiek, G. Stegeman, T. Pertsch, F. Lederer, Y. Min, and W. Sohler. Arrays of weakly coupled, periodically poled lithium niobate waveguides: beam propagation and discrete spatial quadratic solitons. OptoElectron. Rev., 13:113, 2005.
 M. Karpinski, C. Radzewicz, and K. Banaszek. Experimental characterization of threewave mixing in a multimode nonlinear KTiOPOâ waveguide. Appl. Phys. Lett., 94:181105, 2009.
 M. F. Saleh, B. E. A. Saleh, and M. C. Teich. Modal, spectral and polarization entanglement in guidedwave parametric downconversion. Phys. Rev. A, 79:053842, 2009.
 A. Christ, K. Laiho, A. Eckstein, T. Lauckner, P. J. Mosley, and C. Silberhorn. Spatial modes in waveguided parametric downconversion. Phys. Rev. A, 80:033829, 2009.
 A. B U’Ren, K. Banaszek, and I. A Walmsley. Photon engineering for quantum information processing. Quantum Info. Comput., 3(7):480–502, October 2003. ACM ID: 2011567.
 D. K. Serkland, M. M. Fejer, R. L. Byer, and Y. Yamamoto. Squeezing in a quasiphasematched LiNbOâ waveguide. Opt. Lett., 20:1649, 1995.
 S. Tanzilli, H. De Riedmatten, H. Tittel, H. Zbinden, P. Baldi, M. De Micheli, D. B. Ostrowsky, and N. Gisin. Highly efficient photonpair source using periodically poled lithium niobate waveguide. Electronics Letters, 37:26, 2001.
 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., 108:053601, 2012.
 F. Lederer, G. I. Stegeman, D. N. Christodoulides, G. Assanto, M. Segev, and Y. Silberberg. Discrete solitons in optics. Physics Reports, 463:1–126, 2008.
 H. B. Perets, Y. Lahini, F. Pozzi, M. Sorel, R. Morandotti, and Y. Silberberg. Realization of quantum walks with negligible decoherence in waveguide lattices. Physical Review Letters, 100(17):170506, 2008.
 Y. Bromberg, Y. Lahini, R. Morandotti, and Y. Silberberg. Quantum and classical correlations in waveguide lattices. Physical Review Letters, 102(25):253904, 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, and A G White. Twophoton quantum walks in an elliptical directwrite waveguide array. New Journal of Physics, 13(7):075003, July 2011.
 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. O’Brien. Quantum walks of correlated photons. Science, 329:1500, 2010.
 L. Sansoni, F. Sciarrino, G. Vallone, P. Mataloni, A. Crespi, R. Ramponi, and R. Osellame. Twoparticle bosonicfermionic quantum walk via integrated photonics. Physical Review Letters, 108(1):010502, 2012.
 Andrea Crespi, Roberto Osellame, Roberta Ramponi, Vittorio Giovannetti, Rosario Fazio, Linda Sansoni, Francesco De Nicola, Fabio Sciarrino, and Paolo Mataloni. Anderson localization of entangled photons in an integrated quantum walk. Nature Photonics, 7(4):322–328, April 2013.
 M. Gräfe, A. S. Solntsev, R. Keil, A. A Sukhorukov, M. Heinrich, A. Tünnermann, S. Nolte, A. Szameit, and Y. S Kivshar. Biphoton generation in quadratic waveguide arrays: A classical optical simulation. Scientific Reports, 2:562, 2012.
 R. Keil, A. PerezLeija, F. Dreisow, M. Heinrich, H. MoyaCessa, S. Nolte, D.N. Christodoulides, and A. Szameit. Classical analogue of displaced fock states and quantum correlations in glauberfock photonic lattices. Physical Review Letters, 107(10):103601, 2011.
 M. Bellec, G. M. Nikolopoulos, and S. Tzortzakis. Faithful communication hamiltonian in photonic lattices. Optics Letters, 37(21):4504â4506, 2012.
 H. Trompeter. Discrete optics in inhomogeneous waveguide array. PhD thesis, University of Jena, 2006.
 K. G. Makris and D. N. Christodoulides. Method of images in optical discrete systems. Physical Review E, 73(3):036616, 2006.
 A. Szameit, D. Blömer, J. Burghoff, T. Pertsch, S. Nolte, and A. Tünnermann. Hexagonal waveguide arrays written with fslaser pulses. Applied Physics B, 82(4):507â512, March 2006.
 A. Szameit, T. Pertsch, F. Dreisow, S. Nolte, A. Tünnermann, U. Peschel, and F. Lederer. Light evolution in arbitrary twodimensional waveguide arrays. Physical Review A, 75(5):053814, 2007.
 Charles Kittel. Introduction to Solid State Physics. Wiley, 8 edition, November 2004.
 K. J. Blow, Rodney Loudon, Simon J. D. Phoenix, and T. J. Shepherd. Continuum fields in quantum optics. Physical Review A, 42(7):4102, October 1990.
 W. P. Grice and I. A. Walmsley. Spectral information and distinguishability in typeII downconversion with broadband pump. Phys. Rev. A, 56:1627, 1997.
 A. M BraÅczyk, T C Ralph, W. Helwig, and C. Silberhorn. Optimized generation of heralded fock states using parametric downconversion. New Journal of Physics, 12(6):063001, June 2010.
 Agata M. BraÅczyk. Nonclassical states of Light. PhD thesis, University of Queensland, August 2010.
 Andreas Christ, Benjamin Brecht, Wolfgang Mauerer, and Christine Silberhorn. Theory of quantum frequency conversion and typeII parametric downconversion in the highgain regime. New Journal of Physics, 15(5):053038, May 2013.
 R. C. Alferness, R. V. Schmidt, and E. H. Turner. Characteristics of tidiffused lithium niobate optical directional couplers. Applied Optics, 18(23):4012â4016, December 1979.
 J. Ctyroky, M. Hofman, J. Janta, and J. Schrofel. 3d analysis of LiNbOâ: ti channel waveguides and directional couplers. IEEE Journal of Quantum Electronics, 20(4):400 â409, April 1984.
 A. Yariv. Quantum Electronics 3e. John Wiley & Sons, 3. auflage edition, January 1989.
 T. Tamir, editor. Integrated Optics. SpringerVerlag, 2 edition, September 1979.
 D. H. Jundt. Temperaturedependent sellmeier equation for the index of refraction, $n_e$, in congruent lithium niobate. Optics Letters, 22:1553–1555, 1997.