Spatio-spectral characteristics of parametric down-conversion in waveguide arrays
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 down-conversion (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.Wi
The process of parametric down-conversion (PDC)  provides a versatile resource for the creation and engineering of sophisticated quantum states of light. In particular, the creation of photon-pair 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 , polarization [11, 12] and energy-time [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 , and are grouped so close to each other, that the created PDC photons are able to couple evanescently to the nearest-neighbour 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 well-known tool for state engineering in bulk crystals . PDC in WGAs, however, benefits from the bright and efficient photon-pair production in waveguides [23, 24], and offers several further advantages: these systems are not subjected to any in-coupling losses and thus maintain the process efficiency and quantum features. Further, a miniaturized, on-chip 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 .
The linear and nonlinear characteristics of WGAs have attracted enormous interest for the investigation of classical wave phenomena . In the context of quantum optics these systems enable the implementation of quantum walks [27, 28, 29]. In combination with non-classical 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 photon-pair states , displaced Fock states , and faithful quantum state transfer .
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 spatio-spectral 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 spatio-spectral 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 phase-matching 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 type-I PDC in nonlinear WGAs by Solntsev et al. in  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  altering the propagation vector in the -direction
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 waveguides111A more detailed introduction to transverse momenta of particles and bandstructures may be found in  or .. The effective refractive index characterizes the dispersion of a light field which propagates inside a single waveguide taking into account material and modal dispersion contributions . The term corresponds to the modification of dispersion due to the presence of the WGA. Hereby the impact of the WGA on the dispersion is given by the coupling parameter , which is firstly dependent on the distance between the individual waveguide channels and secondly impacted by the frequencies of the propagating fields (see A). The signal and idler fields in the th waveguide can be expressed in the form
where collects all constants . 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
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
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 first-order perturbation theory
We have extended the integration boundaries to plus and minus infinity, since we regard the PDC state long after the interaction in the crystal . We further post-select 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 two-photon PDC state emerging from the WGA as
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 spatio-spectral amplitude of the generated signal () and idler () photons. It assumes the form
in which the phase-mismatch is defined as
The four terms determining the properties of the two-photon 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 phase-matching 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 Spatio-spectral properties of PDC in WGAs
3.1 Origin of the spatio-spectral correlations between the created photon pairs
The key for understanding the spatio-spectral properties of the generated PDC state lies in the phase-matching function . Significant down-conversion rates require nearly perfect phase-matching, satisfying . This constraint induces spatio-spectral correlations between the frequencies and -components of the generated signal-idler photon pairs.
In order to understand its exact nature we split the phase-mismatch in equation (8), into its contributions from the single, isolated waveguide and the dispersion effects from the WGA
The first term describes the spectral phase-mismatch, known from PDC in a single isolated waveguide,
which stems mostly from the dispersion of the material with a contribution from the modal properties of the waveguides . The second term is the phase-mismatch induced by the dispersion of the WGA,
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 phase-mismatch 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 phase-mismatch () the contours (red lines) for different values of run below the central contour, colour-coded in black, (), while for positive spectral phase-mismatch () the contours run above (green lines).
In order to achieve phase-matching in the WGA, i.e , the spectral phase-mismatch has to be compensated by a spatial phase-mismatch of with opposite sign to satisfy . The pairing of these values, and creates the spatio-spectral correlations between the PDC photons. The different contours for which the different values of and yield perfect phase-matching are colour-coded 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 zone222The 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. () 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 circle-like 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 phase-mismatch.
- 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 phase-mismatch 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 type-I 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 spatio-spectral correlations of PDC in WGA, are located in the phase-matching condition . Every spatial correlation pattern, as depicted in figure 2(b)-(d) has a fixed spatial phase-mismatch , which must be compensated by an appropriate spectral phase-mismatch , satisfying . Next, we will illustrate, how to utilize this spatio-spectral 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 phase-mismatch have a significant influence on the emitted photon pairs. Driving the PDC process in the WGA with an adapted spatio-spectral pump shape enables us to engineer a variety of different spatial correlations between the generated photon pairs. The emerging joint spatio-spectral amplitude is given by (7)
This means we can regard as a function which, via multiplication, selects a subset of the spatio-spectral 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 phase-mismatch 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 phase-mismatch 333In this configuration the impacts of the curvature of are negligible.. Due to the spatio-spectral coupling, discussed in section 3.1, this in turn demands that all created photons satisfy in the spatial domain. The resulting -correlations in the created photon pair for a flat spatial pump shape 444This corresponds to pumping in a single waveguide. are depicted by the dashed black lines in figure 3(b)-(b), i.e. a rectangle in figure 3(b) for pumping at the central frequency, outward circles in figure 3(b) in the case of pumping below the central frequency and a circular shape in figure 3(b) for pumping above the central frequency. The corresponding correlations in real space, obtained via Fourier transformation, are depicted by the coloured squares in figure 3(c)-(c), where each square represents a different signal-idler output channel combination. Here is the correlation function for measuring two photons at positions , as defined in B. They vary from an X-shape in figure 3(c) to the square-like555For higher spectral phase-mismatch the shapes are of a circular form. shapes sketched in figure 3(c)-(c).
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 colour-coded 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 phase-mismatch , as given by colour-coded 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 effects666Due to the Fourier limitation, the finite width of the spatial pump will lead to finite widths in the real space correlations as well., as already discussed in , but it is also evident that steering the created photon pairs to the right or to the left is possible as well.
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 signal-idler 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).
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 . 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). Post-selecting on both photons passing their individual filters, this procedure modifies the output state to
A particular example is visualized in figure 5, for non-degenerate photon-pair 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: Spatio-spectral correlations
In summary, we explored the spatio-spectral structure of PDC in WGAs with special attention to spectrally dependent coupling effects. Our investigations revealed that the spatio-spectral correlations are fundamentally determined by the phase-mismatch . 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 phase-mismatch .
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, type-I, 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 InGaAs-detector array sensitive down to the few-photon level. Note that by using a two-dimensional detector array, this configuration allows us to simultaneously resolve the properties of impinging light both spectrally and spatially. Thus, we can analyze any spatio-spectral 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 continuous-wave 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 second-harmonic 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 spatio-spectral 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 phase-matching 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 wavelength-dependent 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 \unit1520-1600\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 wavelength-dependent 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 non-degenerate 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 phase-matching 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 photon-number 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 spatio-spectral characteristics, as discussed in section 3. As such, our measurements are equivalent to measuring many photon-pair 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 phase-matching 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 phase-matching curve. Small discrepancies arise from the accuracy at which the degeneracy point can be determined, limited resolution of the home-made spectrometer and fast decrease of its sensitivity at higher wavelengths. The splitting of the signal and idler spectra into two non-overlapping spectral regions caused by the curvature of the phase-matching function is already visible when tuning the pump wavelength by only a fraction of a nanometer from the degeneracy.
In regard to the spatio-spectral characteristics, we find a very good agreement between the theory for photon pairs and our experiment, even though our detector lacks true single-photon 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.
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 phase-mismatch . A phase-mismatch introduced in the spectral domain is compensated via an opposite phase-mismatch in the spatial domain , leading to spatio-spectral correlations in the generated quantum states. Thus, being able to select only specific regions of the phase-matching, by means of appropriate spectral filtering or by choosing the suitable spectral and spatial characteristics of the pump beam, we can modify the joint spatio-spectral amplitude. Consequently, we are able to engineer the emerging spatial correlations of the two-photon 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 two-photon 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 phase-matching 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 spatio-spectral correlations. Our results are important for the advancement of future applications such as quantum walks or other sophisticated optical networks on WGAs.
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 13-33906S. 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 . 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 
Here is a constant depending on the mode number and is the refractive index difference between the waveguide and the surrounding bulk material  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
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 . 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
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 Spatio-spectral correlation function
The most straightforward way to obtain information about the joint spatio-spectral 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 two-photon PDC state in (6) is given as:
By tracing over the spectral degree of freedom, the correlation function in -space can be written in the form
In the real space the correlation function is determined as
in which , and the bi-photon state in real space, is given by a two-dimensional Fourier transformation of in -space.
Appendix C Numerical simulation of the WGA
In the scope of this paper we consider extraordinarily polarized type-I PDC in lithium niobate using the dispersion relations from . 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 .
We adjust the signal and idler frequencies with the help of quasi phase-matching , where a periodic sign change of the nonlinearity in the material introduces a modification in the spectral phase-matching condition
In the simulation is an effective grating period, fitting the theoretically predicted wavelengths of the phase-matched PDC process to the experimental values.
|Properties of WGA (numerical simulation)|
|Waveguide length [\unit\metre]||0.04|
|Total number of waveguides on the array||41|
|Effective coupling parameter|
|Effective damping coefficient [m]|
|Phase engineering (Figure 4)|
|Phase-matched 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)|
|Phase-matched wavelength pair (, ) [nm]||(1400.0,1600.0)|
|Effective coupling parameter|
|Effective damping coefficient [m]|
|Experimental comparison (Figure 7(e-h))|
|Phase-matched wavelength pair (, ) [nm]||(1549.8,1549.8)|
|Spectral pump bandwidth [nm]||0.5|
|Near-degenerate correlation calculation (Figure 12)|
|Phase-matched 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 phase-matching function, serving as orientation in the spectral domain. Figures 12(b)-(b) depict the selected spatial correlation in -space in accordance to the spatio-spectral 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.
-  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 two-photon entanglement for spontaneous parametric down-conversion 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 space-time 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 down-conversion. 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 parametric-down 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. Molina-Terriza, 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. Franke-Arnold, D. G. Ireland, R. W. Boyd, S. M. Barnett, and M. J. Padgett. Quantum correlations in optical angle-orbital angular momentum variables. Science, 662:329, 2010.
-  P. J. Mosley, A. Christ, A. Eckstein, and C. Silberhorn. Direct measurement of the spatial-spectral structure of waveguided parametric down-conversion. Phys. Rev. Lett., 103:233901, 2009.
-  P. G. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, A. V. Sergienko, and Y. Shih. New high-intensity source of polarization-entangled photon pairs. Phys. Rev. Lett., 75:4337, 1995.
-  Z. H. Levine, J. Fan, J. Chen, and A. L. Migdall. Polarization-entangled 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. Two-photon interference in a mach-zehnder 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 down-conversion 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. Photon-pair 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 photon-pair 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. Opto-Electron. Rev., 13:113, 2005.
-  M. Karpinski, C. Radzewicz, and K. Banaszek. Experimental characterization of three-wave 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 guided-wave parametric down-conversion. 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 quasi-phase-matched 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 photon-pair 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. Two-photon quantum walks in an elliptical direct-write 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. Two-particle bosonic-fermionic 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. Perez-Leija, F. Dreisow, M. Heinrich, H. Moya-Cessa, S. Nolte, D.N. Christodoulides, and A. Szameit. Classical analogue of displaced fock states and quantum correlations in glauber-fock 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 fs-laser 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 two-dimensional 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 type-II down-conversion 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 down-conversion. New Journal of Physics, 12(6):063001, June 2010.
-  Agata M. Brańczyk. Non-classical 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 type-II parametric down-conversion in the high-gain regime. New Journal of Physics, 15(5):053038, May 2013.
-  R. C. Alferness, R. V. Schmidt, and E. H. Turner. Characteristics of ti-diffused lithium niobate optical directional couplers. Applied Optics, 18(23):4012–4016, December 1979.
-  J. Ctyroky, M. Hofman, J. Janta, and J. Schrofel. 3-d 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. Springer-Verlag, 2 edition, September 1979.
-  D. H. Jundt. Temperature-dependent sellmeier equation for the index of refraction, $n_e$, in congruent lithium niobate. Optics Letters, 22:1553–1555, 1997.