Spatio-spectral characteristics of parametric down-conversion in waveguide arrays

Spatio-spectral characteristics of parametric down-conversion in waveguide arrays

R. Kruse, F. Katzschmann, A. Christ, A. Schreiber,
S. Wilhelm, K. Laiho, A. Gábris, C. S. Hamilton, I. Jex and C. Silberhorn
Applied Physics, University of Paderborn, Warburger Straße 100, 33098 Paderborn, Germany Max Planck Institute for the Science of Light, Günther-Scharowsky Straße 1 - Building 24, 91058 Erlangen, Germany University of Innsbruck, Institute for Experimental Physics/Photonics, Technikerstr. 25d, 6020 Innsbruck, Austria Department of Physics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, Břehová 7, 115 19 Praha, Czech Republic. Wigner Research Centre for Physics, Hungarian Academy of Sciences, H-1525 Budapest, P. O. Box 49, Hungary
August 20, 2019

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.

03.65.-w, 42.50.-p, 42.50.Dv, 42.65.Lm, 42.65.Wi

1 Introduction

The process of parametric down-conversion (PDC) [1] 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 [10], 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 [18], 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 [22]. 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 [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 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 [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 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.

Figure 1: Sketch of PDC in a WGA: a pump photon, in the illuminated waveguide channel decays via the nonlinearity, into a signal and idler photon pair which then couples from waveguide to waveguide until it exits the crystal.

2.1 Electric fields in WGAs

Following the theoretical treatment of type-I 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


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 [26] or [40].. 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 [21]. 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 [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


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


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 [42]. 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 [21]. 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 .

Figure 2: Phase-mismatch in the spectral and spatial domain. (a) In the spectral domain the phase-mismatch consists of curved lines oriented along the axis in frequency space. (b) is defined by for frequencies as depicted in (a). For degenerate frequencies implying identical coupling parameters for signal and idler, is given by circular shapes around the middle or the edges of the Brillouin zone with corresponding color coding. If the signal and idler frequencies are distinct, the patterns get stretched from their symmetric shape in (b) into the forms depicted in (c) and (d) .

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[40]. () 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.

Figure 3: In (a)-(a) we show, how different spectral pump shapes (turquoise lines) are able to select an individual frequency mismatch (dashed grey lines). This enables us to excite photon pairs with a variety of -correlations as depicted by the dashed black lines in (b)-(b). These can further be modified by adapting the spatial pump shape , by displacing the rectangular spatial pump shape in -space, as colour-coded with green, red and blue in (b)-(b). The corresponding photon-pair correlations in real space are depicted by red, green and blue squares in (c)-(c). Due to the finite width of the spatial pump shape, the real space correlations will also have a finite width in -direction. The grey shaded areas in (a)-(c) indicate regions with no PDC signal. C gives a numerical simulation without spatial pump shaping.

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 [15], 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.

Figure 4: Introducing different phase distributions in the spatial pump shape as indicated in (a)-(a) changes the real space correlations (d)-(d). Since the absolute value of the spatial pump shape does not vary, the -space correlations in (b)-(b) remain unchanged, while the phase correlations in (c)-(c) show the influence of the pump phase. Linear phases introduce a shift of the real space correlations, without changing the internal shape of the correlation function, while quadratic phases introduce a ”stretching” of the internal structure.

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, zero-phase distribution in , the resulting real space correlations in figure 4(d) show the anti-bunching 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). Post-selecting on both photons passing their individual filters, this procedure modifies the output state to

Figure 5: Filtering specific wavelengths from the PDC spectrum as indicated by the squares in (a)-(a) selects different features in -space (b)-(b). The corresponding correlations in real space are depicted in (c)-(c), where labels the waveguide channel.

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 .

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, 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.

Figure 6: Schematical picture of the experimental setup for measuring the spectral and spatial distributions of the PDC photons created in a PPLN WGA. Abbreviations: SM = spectrometer; HWP = half-wave plate; AL = aspheric lens; WG = waveguide; SF = spectral filter; InGaAs-DA = InGaAs Detector array. For more details see text.

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.

Figure 7: Spectrally and spatially resolved PDC emission from the WGA. The contour plots illustrate the PDC light in different waveguide channels with respect to the wavelength. The spectra on right hand side show the marginal distribution of the central waveguide channel. The measured data was recorded at the pump wavelengths near (a) \unit774.9\nano\meter (b) \unit774.5\nano\meter (c) \unit774.2\nano\meter and (d) \unit773.9\nano\meter. The numerical simulations in (e)-(h) correspond to the measurements in (a)-(d), respectively. The solid lines illustrate the marginals, whereas the dotted lines are fits to the data to obtain the maxima of the spectral 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.

Figure 8: The spatial marginal distributions deduced from figures 7 (d) and (h) for the upper and lower spectral branches. The measured spatial marginals are illustrated at the wavelength regions of (a) \unit1550-1750\nano\meter and (b) \unit1350-1550\nano\meter. The simulated spatial marginals are shown in (c) and (d) for the same spectral regions as in (a) and (b), respectively. The intensity in the central waveguide channel is in each case normalized to unity.
Figure 9: The phase-matching curve of the WGA (squares) extracted from the marginal spectra in figure 7(a)-(d). The error bars in the signal and idler wavelength are given by the full width at half maximum (FWHM) of the measured marginal distributions. Circles correspond to the simulated values including the data shown in figure 7(e)-(h), whereas the black solid lines and the shaded gray area provide guides for the eye and represent the maxima and FWHM of the theoretical marginal spectra, respectively. As an effect of the wavelength dependent coupling, the theory predicts that the upper wavelength branch splits in two spectral regions.

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.

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 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.

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 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 [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.

Figure 10: Form of the electric fields for a rectangular refractive index profile. The overlap of the electric fields determines the coupling parameter.

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]


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


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


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 [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 phase-matching [21], 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.

The wavelength dependent coupling parameter is fitted to the experimental results using the description in A. Its exact form is depicted in figure 11.

Figure 11: Wavelength dependence of the coupling parameter used to fit the numerical simulations to the experimental results.

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)
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
Table 1: Parameters for the simulation of the spatial correlations in figure 5 and the modelling of the experimental results in figure 7 and 9.

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.

Figure 12: Numerical simulation of the output correlation patterns in the WGA by exploiting the tunability of the correlations via the spectral pump shape. (a)-(a) show the corresponding joint spectral amplitudes of the process, illustrating the different pump frequencies. The inlays give the position of the pump shape in relation to the phase-matching function. The resulting correlation functions in -space are depicted in figures (b)-(b), according to the discussed spatio-spectral interplay in figure 2(a) and (b). In figures (c)-(c) the real space correlation function of the WGA is shown. We observe a distinct X-like shape in figure (c), while the coincidence intensity in figures (c) and (c) is localized around the pumped waveguide.

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.



  • [1] D. C. Burnham and D. L. Weinberg. Observation of simultaneity in parametric production of optical photon pairs. Phys. Rev. Lett., 25:84, 1970.
  • [2] 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.
  • [3] 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.
  • [4] 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.
  • [5] 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.
  • [6] A. Christ, K. Laiho, A. Eckstein, K. N. Cassemiro, and C. Silberhorn. Probing multimode squeezing with correlation functions. New. J. Phys., 13:033027, 2010.
  • [7] A. Mair, A. Vaziri, G. Weihs, and A. Zeilinger. Entanglement of the orbital angular momentum states of photons. Nature, 412:313, 2001.
  • [8] G. Molina-Terriza, J. P. Torres, and L. Torner. Twisted photons. Nature Phys., 3:305, 2008.
  • [9] 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.
  • [10] 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.
  • [11] 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.
  • [12] 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.
  • [13] 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.
  • [14] 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.
  • [15] 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.
  • [16] 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.
  • [17] 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.
  • [18] 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.
  • [19] 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.
  • [20] 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.
  • [21] 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.
  • [22] 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.
  • [23] 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.
  • [24] 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.
  • [25] 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.
  • [26] 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.
  • [27] 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.
  • [28] Y. Bromberg, Y. Lahini, R. Morandotti, and Y. Silberberg. Quantum and classical correlations in waveguide lattices. Physical Review Letters, 102(25):253904, 2009.
  • [29] 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.
  • [30] 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.
  • [31] 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.
  • [32] 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.
  • [33] 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.
  • [34] 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.
  • [35] M. Bellec, G. M. Nikolopoulos, and S. Tzortzakis. Faithful communication hamiltonian in photonic lattices. Optics Letters, 37(21):4504–4506, 2012.
  • [36] H. Trompeter. Discrete optics in inhomogeneous waveguide array. PhD thesis, University of Jena, 2006.
  • [37] K. G. Makris and D. N. Christodoulides. Method of images in optical discrete systems. Physical Review E, 73(3):036616, 2006.
  • [38] 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.
  • [39] 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.
  • [40] Charles Kittel. Introduction to Solid State Physics. Wiley, 8 edition, November 2004.
  • [41] 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.
  • [42] 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.
  • [43] 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.
  • [44] Agata M. Brańczyk. Non-classical states of Light. PhD thesis, University of Queensland, August 2010.
  • [45] 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.
  • [46] 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.
  • [47] 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.
  • [48] A. Yariv. Quantum Electronics 3e. John Wiley & Sons, 3. auflage edition, January 1989.
  • [49] T. Tamir, editor. Integrated Optics. Springer-Verlag, 2 edition, September 1979.
  • [50] D. H. Jundt. Temperature-dependent sellmeier equation for the index of refraction, $n_e$, in congruent lithium niobate. Optics Letters, 22:1553–1555, 1997.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description