Supermode-based second harmonic generation in a nonlinear interferometer

We demonstrate supermode-based second harmonic generation in an integrated nonlinear interferometer made of linear and nonlinear directional couplers. We use a fully-fibered pump shaper to demonstrate second harmonic generation pumped by the symmetric or anti-symmetric fundamental spatial modes. The selection of the pumping mode and thus of a specific SHG spectral profile is achieved through the selection of the fundamental wavelength and via a robust phase setting scheme. We use two methods: either post-selecting or actively setting the pumping mode. Such a modal phase matching paves the way for classical and quantum applications of coupled nonlinear photonic circuits, where multimode excitation, encoding and detection are a route for multiplexing and scaling up light-processing.


Introduction
The nonlinear directional coupler (NDC) is a core device in integrated optics without analog in bulk optics: evanescent coupling and nonlinear interaction effects are intertwined instead of being sequential as it is the case for daisy-chained beamsplitters and active nonlinear elements.The NDC has been widely studied in the classical regime for ultrafast all-optical switching [1][2][3][4][5][6][7][8][9].Recently, the  (2) NDC has found a flourishing field of application: quantum optics [10].
Its key strengths in quantum information processing as a source of entangled photons and entangled field quadratures are actively explored [11][12][13][14][15][16].Particularly, strong bipartite and four-partite continuous-variable entanglement are predicted in a second harmonic generation (SHG) configuration when pump light at the fundamental frequency is in propagation eigenmodes -i.e. in the spatial supermodes (see Figure 1) [17]-of the NDC [18].Such SHG schemes are especially interesting as they can harness widely available telecom domain sources to generate spatial entanglement at visible wavelengths compatible with atomic memories [19].In addititon, they offer the possibility to generate in a controlled way multipartite entangled states, that are key resources for quantum technologies such as measurement-based quantum computing, multiuser quantum communications and quantum metrology [20][21][22][23][24].In this paper we experimentally study SHG pumped by supermodes for the first time.
SHG in waveguides is commonly produced by quasi-phasematching [25].This method is based on periodical reset of wavevector phase mismatch to maintain a coherent buildup of the nonlinear interaction and it is conventionally obtained by periodic modulation of the nonlinear arXiv:2106.13204v2[physics.optics]27 Nov 2021 coefficient (periodic poling) [26,27].Ultra-high SHG conversion efficiencies in lithium niobate waveguides have been recently reported [28][29][30].In arrays of evanescently coupled waveguides, the poling period is usually fixed in order to fullfil the phasematching condition at a given wavelength for the mode propagating in each waveguide.However, considering instead optical pump in spatial supermodes is a fruitful approach for multimode light processing in coupled waveguides and makes new and different phasematching conditions arise [31,32].The SHG is then generated efficiently at different wavelengths depending on the propagating supermode [12,33].Supermode-based SHG has however remained elusive up to now due to the difficulty in exciting specific supermodes at the input of the array.In this paper we overcome those experimental issues for the NDC, an array of two nonlinear waveguides, through a fully-fibered pump shaper and demonstrate supermode-based SHG through a specifically-designed integrated nonlinear interferometer made of linear (LDC) and nonlinear (NDC) directional couplers.Note that the term nonlinear interferometer corresponds here to a linear interferometer with a phase changed due to the nonlinear effect [34], in contrast with SU(1,1) nonlinear interferometers or similar schemes [35,36].The paper is organized as follows: in section 2 we introduce the theoretical formalism that governs the SHG light in the nonlinear interferometer including the phase mismatch and the different pump configurations.In section 3 we detail the experimental setup.We exhibit experimental SHG generation in NDC via engineered pumping to access and harness the supermodes of the device in section 4. Finally, the main results of this work are summarized in section 5.

Theoretical model
Our device is composed by three sections of coupled waveguides produced by High Vacuum Proton Exchange (HiVacPE) (Figure 2): a periodically poled lithium niobate waveguides-based NDC embedded between two LDCs.In the three sections, the energy of fundamental fields propagating in each waveguide is exchanged between the coupled waveguides through evanescent waves, whereas the interplay between second harmonic fields is negligible for the considered propagation lengths due to their high confinement into the waveguides.We consider a single polarization mode in this model, since only TM fields can propagate in Z-cut HiVacPE lithum niobate waveguides.The input-output LDCs are engineered to act as 3-dB beam splitters at the fundamental wavelength.The fundamental fields after passing a 3-dB LDC with two identical waveguides are where  and  are the slowly-varying amplitudes of fundamental (f) modes propagating in the upper (a) and lower (b) waveguides, respectively, and  is the coordinate along the direction of propagation.The input LDC allows to generate fields along specific axes of a Bloch sphere at the input of the NDC section at  = 15: i) pumping single waveguides i.e. (1, 0)  or (0, 1)  2, for respectively Δ = 0 and  (see Figure 1) [17].
The NDC section is made of two identical waveguides with  (2) nonlinear response due to periodic poling (PPLN).In each waveguide, the input fundamental field at frequency   is up-converted into a second-harmonic field at frequency  ℎ = 2  (SHG).The physical processes involved in the NDC, i.e. fundamental-fields evanescent coupling and nonlinear generation, are described in the perturbative regime by the following system of equations [5] where the subindex (h) refers to second harmonic fields,  is the nonlinear constant, proportional to  (2) and to the spatial overlap of the fundamental and harmonic fields in each waveguide,  the linear coupling constant, and Δ ≡ ( ℎ ) − 2(  ) − 2/Λ is the wavevector phase mismatch with () the propagation constant at wavelength  and Λ the poling period.This system of equations is non integrable and thus does not present analytical solutions in general [37], although analytical solutions in terms of Jacobi elliptic functions arise under symmetric pumping conditions [38].
Rewriting the Eqs (2) for a NDC in terms of fundamental slowly-varying supermodes of a LDC,   = (   +   ) and   are fundamental supermodes, symmetric and antisymmetric, of the LDC region of the device, and can be also used as a suitable basis for the NDC region (Figure 1).In contrast, and as the fields  ℎ and  ℎ are not coupled via evanescent fields,  ℎ and  ℎ are just a convenient change of variable, they are not supermodes of either the NDC or LDC parts of the device.Note that in the case of non-negligible evanescent coupling of harmonic modes -e.g. in the case of nanowaveguide arrays [30]-harmonic supermodes would appear modifying the dynamics of the system.
The nonlinear and coupled Eqs (3) do not have analytical solutions in general, but show two phasematching conditions for the fundamental supermodes.Indeed, at the fundamental wavelength   =   , fulfilling the quasi-phase matching condition Δ(  ) − 2 = 0, the nonlinear interaction between the supermode field   and the field  ℎ builds up coherently and Eqs (3) can be approximated by for wavelengths around   .A similar condition is obtained for the fields   and  ℎ at the fundamental wavelength   =   such that Δ(  ) + 2 = 0 and for wavelengths around   .The above approximations are correct as long as the coupling constant  and the length of the NDC are large enough to avoid overlapping with the phasematching condition for the fundamental single modes   (  ), i.e.Δ(  ) = 0, around the wavelength   =  single .Note that  ℎ is coupled to both fundamental supermodes   and   in Eqs (3) with phasematching at Δ(  ) = 0.However, due to the spatial orthogonality of the fundamental supermodes, it is not possible to excite both at the same time, so the generation of an antisymmetric linear combination of harmonic modes  ℎ is not physically feasible.Eqs (4) and ( 5) are formally equivalent to the well-known equations of single mode SHG and have analytical solutions [26].In the low conversion regime the fundamental fields remain essentially constant over the interaction length, and the solutions for the harmonic supermodes when pumping the even and odd fundamental supermodes are respectively Then it is splitted in two by a fiber-beam splitter (fbs).Two fiber polarizers (P1 and P2) control the relative power injected in each waveguide that guide only TM-polarized modes.A fiber phase shifter (PS) controls the relative phase of the injected fields.The injection into the nonlinear interferometer is carried out via an AR-coated V-groove array of fibers.SHG fields (red) are produced on chip.The two spatial modes at fundamental and harmonic frequencies are collected with a lens (L) and frequency-demultiplexed with a dichroic mirror (DM).The fundamental (harmonic) fields, drawn in blue (red) are then directed to two photodiodes (Pd) sensitive to telecom (visible) frequencies with the help of a square mirror (M).Lenses in front of each photodiode are used to collect all the light (not shown).The four detectors (Pd1-Pd4) characterize completely the physical process in the nonlinear interferometer.
Thus, equal pump powers with phase difference Δ = 0 () [39] at the input of the NDC PPLN region at wavelengths   (  ) , i.e. excitation of the   (  ) supermode, generates efficiently a symmetric linear combination of harmonic fields  ℎ around wavelengths  ℎ =   () /2.The mode  ℎ builds up coherently when pumping at wavelengths   =   () and the individual modes  ℎ and  ℎ of the individual single waveguides strive at   =  single .

Experiment
In practice the nonlinear interferometer of Figure 2 has been implemented in high vaccum proton exchanged periodically poled Lithium Niobate (HiVacPE PPLN) 6 µm wide-waveguides [40].
As described in section 2, such a system is governed by the nonlinearity  of the periodically poled region and the evanescent coupling  between the coupled waveguides.To experimentally determine the linear coupling constant , we used arrays of 50-100 waveguides fabricated under the same conditions for different separations  between waveguides.We obtained the following empiric dependance of the coupling constant as a function of the distance between the waveguides  () = 16.4  −/3.8 −1 for  ≥ 13  [41].We designed the device to allow the use of commercial V-groove arrays of fibers with 127  spacing between fibers for input and output.
Our device has thus a -shaped adiabatic region where the distance between the waveguides goes from 127  to 16  within a total length of 15 mm.The coupling strength for  = 16  is  = 0.25  −1 deduced from the empirical law.Taking into account the coupling strength of the region where  > 16 , the LDC could be approximated by an effective coupler with constant separation between waveguides and an effective length  eff = 9.3 .The corresponding coupling is such that   eff ≈ 3/4 and it is correctly described by the transfer matrix of Eq. (1).Our integrated nonlinear interferometer thus presents a NDC region embedded between two 3-dB LDCs (Figure 2) where we obtain a 3-dB performance in the LDCs thanks to the coupling constant and bends engineering.
The NDC section is a 20 mm periodically poled region.The poling period is Λ = 16.4 .It has been chosen for SHG of  single = 1.56  radiation to produce 0.78 m field at a temperature phase matching of 59.1 • C, stabilized within 0.1 • C with a resistive heater controlled through a PID controller (Thorlabs).The typical propagation losses of our waveguides are 0.2 dB cm −1 at fundamental wavelength and 1 dB cm −1 at SH wavelength.The device (daisy-chained LDC+NDC+LDC) total length is 5 cm.To fully characterize our device, we need to determine the nonlinear strength .This is done by measuring the SHG efficiency in a 6-m-wide single waveguide of similar PPLN length, obtaining    = 42% W −1 cm −2 , which corresponds to  ≈ 0.07 mm −1  −1/2 .The refractive indices for fundamentals and harmonics modes are calculated with a finite element method simulation (COMSOL) using the Sellmeier coefficients for bulk lithium niobate and a model of the refractive index of the HiVacPE waveguides [42].
With these indices we calculate the Δ we use in our theoretical model.
The experimental arrangement used to study supermode-based SHG in a nonlinear interferometer is depicted in Figure 3.It consists of cw tunable laser (TUNICS T100R) emitting 10 mW between 1490 nm and 1620 nm.The optical radiation is amplified with an Er-doped amplifier (KEOPSYS), generating up to 2 watts in the C-band.For the purpose of our investigation, we set the nominal laser power at 500 W (-3 dBm) and amplify it up to 200 mW (23 dBm).The telecom field (fundamental) is split in two paths by a 3dB fiber beam splitter and coupled into the chip by an AR-coated V-groove array of fibers with a pitch that matches the waveguides spacing.Typical coupling efficiencies are estimated to be about 60%.The amplitude of each input is adjusted via 3-paddle polarizers that efficiently decrease the amount of light coupled into the TM guided mode of the HiVacPE waveguides without additional effect.The relative phase between each input is adjusted when necessary via a fiber phase shifter (IDIL FS20).The stability of the phase between inputs is better than /20 over few minutes, which is the typical duration of a full wavelength scan as depicted in Figures 4 to 6.Both the polarization controllers and the phase shifter allow to fully describe the Bloch sphere and excite the NDC in all possible configurations.We collect the output signal from each waveguide with a collimating lens (THORLABS C-220 TMD-C) having a transmission > 85% at the SHG wavelength and > 99% at the fundamental wavelength.We separate the output fundamental fields from the harmonic fields using a dichroic mirror (THORLABS DMSP950).The two reflected fundamental beams are further separated by means of a gold mirror (M2 in Figure 3) and are focused on two InGaAs photodiodes (Newport 1811-FS) using two  = 100 mm lenses.The two transmitted harmonic beams are also separated by means of a gold mirror (M1) and focused on two Si photodiodes (Newport 1801-FS) with two  = 150 mm lenses.To improve the detection sensitivity of the generated SH beams, a 30 % depth intensity modulation is applied to the fundamental laser beam by driving its current with an RF signal.The modulated fundamental and harmonic intensities are then detected using SRS 830 lock-in amplifiers.The frequency of the modulation is fixed at 80 kHz to be compatible with the bandwidth of the lock-in amplifiers.

Pumping a single waveguide
We start our analysis by exciting only one of the inputs of the LDC.As described by Equation (1), due the linear coupling effect, the pump will be present at both inputs of the NDC.The relative powers and the relative phases will be defined by both the propagation length in the LDC and the effective coupling constant.Accordingly, whatever the injected waveguide, the pumps arriving at the NDC inputs have equal amplitudes and a /2 relative phase.Figure 4 shows the SH spectra as the pump wavelength is tuned from 1.558 m to 1.562 m, exciting either waveguide  or waveguide .The SHG spectra are measured simultaneously at outputs  and  in both configurations.The two top plots correspond to the injection into one waveguide and the measurement of the SH spectra at the two output waveguides and the bottom plots are obtained when pumping the other waveguide.The apparent noise in the spectra is a Fabry-Perot periodic oscillation effect, not resolved in this scale, due to the Fresnel reflection at the input and output facets of the sample (∼14% of reflectivity).The dashed black lines are theoretical results obtained from Eq.( 2) for  0 = 0.255 ± 0.005 mm −1 and  = 0.07 mm −1  −1/2 , in excellent agreement with the measurements.This experimental value of  0 is close to the one obtained from the empirical expression deduced from other experimental measurements, and in excellent agreement with the value of design thus demonstrating our ability to engineer the coupling constant.We note that our theoretical model fits with a high degree of accuracy the direct pumping configuration input  → output  (or  → ).However, there is a small wavelength mismatch for the indirect pumping configuration input  → output  (or  → ).A small change in the phase mismatch Δ used in our simulations can fit the indirect pumping traces at the cost of degrading those related to direct pumping.Our theoretical model considers the temperature as homogeneous along the NDC.We hypothesize that our experimental result shows the impact of temperature inhomogeneity in the sample on the SHG spectra and the sensitivity of SHG spectra measurement to the temperature stability [43][44][45].
We emphasize here the fact that the theoretical lines in Fig.( 4) describing well the behaviour of the NDC are obtained using Eq.( 2), exciting the two inputs of the NDC with the same pump powers,   =   (|   | = |  |) at  = 15 mm, and with a /2 relative phase.Any deviation from these conditions change drastically the shape of the fits.Thus, this demonstrates our control on the input fields at the NDC.Moreover, this shows our ability to control the LDC design and in particular to achieve a perfect 3 dB integrated beamsplitter, with equal amplitude outputs and relative phase of /2.
The spectra depicted in Fig.( 4) exhibit three peaks.The central peak at the fundamental wavelength   ∼ 1.56 m corresponds to the phasematching of single fundamental modes (i.e.Δ = 0).The two side peaks appear as a consequence of the interplay between the phase mismatch at wavelengths   ≠   and the beating of the fundamental mode amplitudes.Thus, the wavelengths of the side peaks are strongly dependent on the coupling coefficient.Note that in this configuration there is no modal phasematching: the fundamental supermodes   (  ) are never excited due to the /2 relative phase between the fundamental single modes   and   .Indeed, the fundamental field amplitudes oscillate between (1, )  / √ 2, (0, 1)  , (1, −)  / √ 2, and (1, 0)  along propagation in the NDC due to the coupling.This effect produces a -dependent phase mismatch between the fundamental modes and the generated harmonic modes.For wavelengths   ≠   the energy flows forth and back from the fundamental modes to the harmonic modes at a rate given by the phase mismatch.Thus, for wavelengths very far from   the second harmonic generated is negligible.For the wavelengths where this beating is slow enough and that for the length of our NDC are in a cycle of harmonic generation -and not of parametric downconversion-appear side harmonic peaks.

Pumping both waveguides
Now, we explore SHG of the NDC when exciting simultaneous both inputs of the LDC.The two fundamental fields are obtained after the amplifier using a fibered beamsplitter.With the LDC, they constitute a Mach-Zehner interferometer whose two outputs feed the inputs of NDC (see Figure 3).Ideally, one would like to equalize the NDC input intensities and fix the relative phases to 0 or  in order to excite either the symmetric supermode or the antisymmetric one.However, the LDC that is necessary to on-chip coupling mixes the two inputs prior to the injection into the NDC according to Eq. ( 1).In the following, we introduce the two methods we explored and implemented to circumvent the issue of the LDC that demonstrate that this is an asset rather than a limitation.We first present how to extract the SHG spectra associated to a specific pumping mode via post-selection.For this scheme, we harness the simultaneous monitoring of the relative phase at the input of the NDC provided by the LDC.We further combine the SHG with a control of the phase at the input of the LDC which maps into a power balance control and enables supermode based generation at any fundamental wavelength.

Length-difference phase
We take advantage of the fact that the two arms of the Mach-Zehnder interferometer have, by construction, different lengths.This is due to the small difference in length of the similar input fibers (Figure 3).The fundamental fields experience a wavelength-dependent interference depending on the phase Δ (  ) = 2  Δ/  , where   is the effective refractive index of the fundamental mode propagating in the fiber and Δ is the length difference.Our LDC has been specifically designed such that   eff = 3/4.For such  eff and equal input powers   (0) =   (0), from Eq. (1) the fundamental powers and relative phase at the input of the PPLN NDC region are Thus, according to Eq. ( 8), when tuning the wavelength   , the relative phase of the fundamental fields reaching the NDC inputs is either 0 or , exciting respectively the even   or odd   fundamental supermodes when   =   at  = 15 mm.Conversely for wavelengths where Δ (  ) = /2 (3/2), only one of the two waveguides of the NDC,  or , is injected.In this way we can select the input configuration by only adjusting the fundamental wavelength.Depending on the choice of   the symmetric and antisymmetric fundamental supermodes are excited (Figure 1) and act as pumps for SHG.The excellent control on fabrication parameters allows us to target and fabricate a device functioning at this convenient working point.In order to demonstrate this performance, we first confirm that we pump equally both waveguides checking the ratio of fundamental powers at the output (Pd3 and Pd4 in Fig. 3) when pumping first waveguide  and then .We then sweep the fundamental wavelength around 1.560 m across a 4 -bandwidth with 1 m resolution and collect the SHG at each output.Figure 5 shows the theoretical predictions (upper) and the experimental results (lower).The discrepancy in the amplitudes between theory and experiment can be explained by a wavelength-dependent collection efficiency e.g.caused by an imperfect alignment of the output modes in the collecting lens.The number of oscillations showed is related to the length difference Δ and it has been fitted manually in Figure 5 top-left and top-right to have a similar shape to the experimental spectra.
Using the transmitted fundamental output in the upper waveguide  as a reference, we follow the evolution of the relative phase as the wavelength is tuned.According to Eq.( 8), the maxima indeed correspond to Δ = 0 and thus Δ = 0 (even supermode), and the minima correspond to Δ =  and thus Δ =  (odd supermode).We can thus select and separate the experimental data in the SHG spectra corresponding to the maxima of the transmitted fundamental field (Δ = 0) from those corresponding to the minima (Δ = ).This is reported in Figure 6 in agreement with the theoretical expectations (dashed blue and red lines).The supermode-based SHGs are obtained for two fundamental pump phases with a  difference.Sinc 2 -shaped spectra appear centered at wavelengths equivalent to Δ − 2 = 0 (even) and Δ + 2 = 0 (odd) as expected from Eqs (6-7), thus demonstrating supermode-based SHG.The spurious points in Figure 6 are due to phase jumps inherent to experimental instabilities in the interferometer.This effect could be avoided by actively locking the phases of the fundamental fields.

Phase controlled
In the previous section, the relative phase Δ is tuned by scanning the laser wavelength.Hence, for each wavelength, we only have access to a single supermode related to the value of Δ (  ).We have adopted an alternative strategy to access both supermodes and gain a new tuning parameter to balance pump power inputs.We have included in one arm of the interferometer a voltage-driven phase shifter achieving 2 phase shift with a 10 V triangular signal oscillating at 100 Hz.The full phase range becomes thus accessible for every wavelength.In the following we pump equally both waveguides and tune the phase Δ from −/2 to 3/2 thanks to the phase shifter.
The input phase Δ is thus transformed into an adjustable pumping profile at the NDC region with Δ as  , ( ]) as relative phase Δ at the input of the NDC. Figure 7 shows the experimental results (lower) and the related simulations (upper) for a suitably chosen set of wavelengths.i) At   = 1.5592  (left) we excite the even supermode at the NDC when Δ = 0, producing supermode-based SHG (yellow).This wavelength is highlighted by the solid vertical lines in Fig. 6.At Δ = , however, no SHG is present as the odd supermode is not phase matched at   = 1.5592 .ii) At the phase-matching wavelength of a single waveguide of 1.5601  (dashed vertical lines in Fig. 6), the supermodes are not phase matched, thus almost no SHG is produced both at 0 and  phases.SHG is produced at /2 and 3/2 since then individual modes are excited at the NDC.The physical process is here similar to that discussed in section 4.1.The different height of the SH peaks is related to the effective interaction length seen by each individual fundamental mode when pumping one or the other input of the NDC.iii) Finally, at   = 1.5610 (dotted vertical lines in Fig. 6) we excite the odd supermode at the NDC when Δ = , producing supermode-based SHG.At Δ = 2, however, no SHG is present as the even supermode is not phase matched at this wavelength.Note that the second LDC does not affect the generated second-harmonic fields as they are not evanescently coupled.This experiment demonstrates our ability to produce supermode-based SHG at specific fundamental wavelength working points.This phase analysis at specific wavelength also sheds some light on the SHG observed in the case of the   -phase dependence.Indeed Figure 7 also accounts for the double beating in the central peak of Figure 5: around  single , the two peaks of SHG seen in the middle pannel of Figure 7 contribute whereas only a single peak contributes in the case of the side peaks around   and   .

Conclusions
We have experimentally demonstrated supermode-based second harmonic generation in an integrated nonlinear interferometer, made of linear and nonlinear directional couplers.We used a fibered and integrated pump shaper to engineer the pumping profile and thus exhibit SHG related to the symmetric and antisymmetric fundamental modes.More precisely, the selection of the pumping mode and thus of a specific SHG spectral profile is achieved through the selection of the fundamental wavelength that sets the pumping profile at the input of the non directional coupler.The experimental results match the theoretical predictions with an excellent agreement and without adjustable parameters.Indeed the governing parameters of the device -nonlinearity and coupling -have been characterized independently.
Full characterization is done either by sampling second harmonic emission at the appropriate fundamental frequency working point or by changing the pump amplitude profile at a fixed fundamental wavelength.This latter scheme that actively sets the pumping mode incorporates a fibered phase shifter and a specifically designed linear directional coupler.This scheme maps the driving voltage of the phase shifter to the power balance at the input of the nonlinear direction coupler while enabling a switch between even-supermode or odd-supermode excitation.
In addition to its interest in spatial multimode quantum optics, we stress that this single shot supermode-based approach is fruitful to characterize a coupled nonlinear waveguide array device [46].We have combined on a single chip easy input and output coupling via V-groove fiber arrays and directional couplers and direct addressing of supermodes through a robust scheme that does produce reliably 0 and  dephasing at the input of the nonlinear directional coupler.The generated SH signal results from modal phasematching and is governed by coupling and nonlinearity.Such a multimode and modal approach is thus useful for integrated-optics characterization, especially as the widths, positions and height of the peaks in the SH spectra are degraded and displaced by inhomogenous poling, lack of accuracy for the coupling or inefficient nonlinearity [43][44][45].From a quantum technologies point of view, these results are directly related to the harnessing of nonlinearities [28,47] and supermodes [32,48,49] for quantum information and more specifically the generation and manipulation of advanced entangled states.Indeed, the NDC is the simplest of important devices for quantum information to generate robust entangled multicolor states in the discrete and continuous variable domain [31].Our fiber and integrated excitation scheme to adress the NDC is thus an interesting addition to explore with versatility the supermode-based and individual mode-based entanglement.

Fig. 3 .
Fig.3.Sketch of our experimental setup to measure supermode-based SHG.A Telecom fundamental field (blue) from a wavelength-tunable laser is amplified by an erbium amplifier.Then it is splitted in two by a fiber-beam splitter (fbs).Two fiber polarizers (P1 and P2) control the relative power injected in each waveguide that guide only TM-polarized modes.A fiber phase shifter (PS) controls the relative phase of the injected fields.The injection into the nonlinear interferometer is carried out via an AR-coated V-groove array of fibers.SHG fields (red) are produced on chip.The two spatial modes at fundamental and harmonic frequencies are collected with a lens (L) and frequency-demultiplexed with a dichroic mirror (DM).The fundamental (harmonic) fields, drawn in blue (red) are then directed to two photodiodes (Pd) sensitive to telecom (visible) frequencies with the help of a square mirror (M).Lenses in front of each photodiode are used to collect all the light (not shown).The four detectors (Pd1-Pd4) characterize completely the physical process in the nonlinear interferometer.

Fig. 5 .
Fig. 5. Two-pump-input second harmonic generation in a nonlinear directional coupler with a   -dependent relative phase: theoretical SHG (upper) in waveguides a (left) and b (right), and experimental counterpart (lower) for waveguides a (left) and b (right).

Fig. 6 .
Fig. 6.Experimental supermode-based SHG sampled from the SHG of Figure 5 lower.SHG generated by the even fundamental supermode   (upper) and by the odd fundamental supermode   (lower).The generated harmonic modes  ℎ and  ℎ in blue and red, respectively.Dashed lines correspond to our theoretical model.Vertical lines stand for the phasematching fundamental wavelengths for even supermode (  =   = 1.5592 , solid), individual modes (  =   = 1.5601 , dashed) and odd supermode (  =   = 1.5610 , dotted).The spurious points are due to phase jumps inherent to experimental instabilities in the interferometer.

Fig. 7 .
Fig.7.Theoretical (upper) versus experimental (lower) SHG (yellow) for equal pump power in both waveguides and varying input phase Δ at three phase-matching pump wavelengths   (see vertical lines in Fig.6).Left: even supermode-based SHG for   =   = 1.5592  (Δ − 2 = 0).Center: individual-mode based SHG for   =  single = 1.5601  (Δ = 0).Right: odd supermode-based SHG for   =   = 1.5610  (Δ + 2 = 0).The fundamental field measured at the upper output  which was modulated with a triangular signal at 480 Hz is used as a phase reference (blue).The windows on the x axis have been chosen to display conveniently the peaks as central.

16 µm 127 µm 20 mm 15 mm Leff= 9.3 mm Leff= 9.3 mm 15 mm
Fig.2.Sketch of our nonlinear interferometer in lithium niobate.The input and output ports are 127 m apart to match the interfiber separation in a commercial V-groove array of fibers.Waveguides are adiabatically coupled in 4.5 mm into an 8.5 mm linear directional coupler (LDC) with 16 m center to center separation.The evanescent coupling constant for the straight section is  = 0.25  −1 .We compute the evanescently coupled region of the first 15  of the chip are equivalent to a total effective coupling length of  eff = 9.3 mm (see text).  eff ≈ 3/4 and the LDC acts as a 3 dB beam splitter.The light then encounters a 20-mm long periodically poled region where SHG is produced in a NDC.The poling period is 16.4  and the phase matching temperature 59.1 ºC.The device is symmetric with respect to the center.