• Photonics Research
  • Vol. 10, Issue 7, 1763 (2022)
Pawan Kumar1、2、*, Sina Saravi1, Thomas Pertsch1、2, Frank Setzpfandt1, and Andrey A. Sukhorukov3、4
Author Affiliations
  • 1Institute of Applied Physics, Abbe Center of Photonics, Friedrich Schiller University Jena, 07745 Jena, Germany
  • 2Fraunhofer Institute for Applied Optics and Precision Engineering, 07745 Jena, Germany
  • 3Research School of Physics, Australian National University, Canberra, ACT 2601, Australia
  • 4ARC Centre of Excellence for Transformative Meta-Optical Systems (TMOS), Australia
  • show less
    DOI: 10.1364/PRJ.450410 Cite this Article Set citation alerts
    Pawan Kumar, Sina Saravi, Thomas Pertsch, Frank Setzpfandt, Andrey A. Sukhorukov. Nonlinear quantum spectroscopy with parity–time-symmetric integrated circuits[J]. Photonics Research, 2022, 10(7): 1763 Copy Citation Text show less

    Abstract

    We propose a novel quantum nonlinear interferometer design that incorporates a passive parity–time (PT)-symmetric coupler sandwiched between two nonlinear sections where signal–idler photon pairs are generated. The PT symmetry enables efficient coupling of the longer-wavelength idler photons and facilitates the sensing of losses in the second waveguide exposed to analyte under investigation, whose absorption can be inferred by measuring only the signal intensity at a shorter wavelength where efficient detectors are readily available. Remarkably, we identify a new phenomenon of sharp signal intensity fringe shift at critical idler loss values, which is distinct from the previously studied PT symmetry breaking. We discuss how such unconventional properties arising from quantum interference can provide a route to enhancing the sensing of analytes and facilitate broadband spectroscopy applications in integrated photonic platforms.

    1. INTRODUCTION

    The generation of photon pairs through spontaneous parametric downconversion (SPDC) of a pump photon into signal (s) and idler (i) photons inside a quadratic nonlinear medium is inherently probabilistic and characterized by a quantum mechanical pair-generation probability amplitude [1]. When two such sources of photon pairs constitute a quantum optical system, the resulting signal (and idler) photon amplitude after superposition from the sources does not, in general, show any first-order interference. However, if these sources are pumped by a common coherent pump laser and their idler modes are properly aligned, so that it is impossible to ascertain in which source the photon pair was created, the final signal photon intensity after superposition does show interference [2,3]. Indistinguishability of the two sources with regard to pair generation lies at the heart of this quantum interference effect [35]. The idler mode from the first nonlinear source must pass through the second source to ensure this indistinguishability and “induce” the coherence between their signal modes necessary for first-order interference. Such a configuration of two aligned nonlinear sources of photon pairs is commonly referred to as a quantum nonlinear interferometer [69].

    Recently, quantum nonlinear interferometers have been employed to perform spectroscopic and imaging applications in technologically challenging spectral ranges such as infrared and terahertz [1017]. The exciting aspect of these applications is that detection is only required in the visible spectral range on the shorter wavelength signal photon of correlated pairs emitted from nondegenerate SPDC sources. Most applications thus far have employed bulk optical platforms with nonlinear crystals as the source of photon pairs. An alternative approach is to use photonic integrated circuits to realize on-chip nonlinear interferometers, where nonlinear waveguides serve as the source of photon pairs [1821]. This allows combination of the inherent advantages of integrated platforms, such as higher nonlinear conversion efficiency, smaller device footprint, long-term stability, and scalability, with the application prospects of quantum nonlinear interferometry. On the other hand, integrated systems implementing advanced physical concepts can also enhance the capabilities of optical interferometers. One particularly promising concept explored in integrated optical devices is parity-time (PT) symmetry [22,23]. In the context of coupled waveguide systems, such as directional couplers and coupled resonators, PT symmetry is usually realized by judiciously incorporating balanced optical losses and/or gain in distinct parts of the system in association with a symmetric refractive index distribution [2426]. The intriguing aspect of such a PT-symmetric system is that it is characterized by a phase transition phenomenon symbolized by the existence of a PT symmetry breaking loss strength. This means that its response below and above this exceptional point is qualitatively different and forms the basis for a slew of novel effects [22,23,27,28]. Importantly, in recent works PT symmetry has been put to use to enhance the sensitivity of optical structures to external perturbations [29,30], paving the way for integrated spectroscopic and sensing applications with increased responsivity.

    Here, we propose a conceptually new quantum nonlinear interferometer for spectroscopic applications using an integrated waveguide platform with PT symmetry. Different to the standard spectroscopic implementations, the analyte to be probed forms part of a passive PT-symmetric coupler [22,31]. This coupler is positioned between two sources of photon pairs constituting the interferometer [see Fig. 1(a)]. The coupler offers exceptional control over linear light dynamics such that the output signal intensity can exhibit strong dependence on the idler loss introduced by the analyte in the second waveguide. We show that this opens up new possibilities to tailor the response of the nonlinear interferometer. We identify a new phenomenon of sharply shifting signal interference fringes at specific critical loss magnitudes, which is unrelated to the PT-symmetry breaking effect. This feature appears due to modulation of quantum indistinguishability between the sources of photon pairs and can provide a new route for enhancing the sensing of analytes. The scheme outlined here could also be beneficial in efforts to engineer biphoton quantum states using integrated waveguides and optical-fiber-based nonlinear interferometers [21,3234].

    (a) Sketch of hybrid nonlinear interferometer incorporating a PT coupler for sensing of analyte-induced absorption (γi) in waveguide Wg-2. Dependence of (b) real and (c) imaginary parts of normalized eigenvalues, λ˜±/Ci, which account for the loss and wavenumber respectively, on the normalized loss γi/Ci for the two idler supermodes of symmetric PT coupler composed of Wg-1 and Wg-2.

    Figure 1.(a) Sketch of hybrid nonlinear interferometer incorporating a PT coupler for sensing of analyte-induced absorption (γi) in waveguide Wg-2. Dependence of (b) real and (c) imaginary parts of normalized eigenvalues, λ˜±/Ci, which account for the loss and wavenumber respectively, on the normalized loss γi/Ci for the two idler supermodes of symmetric PT coupler composed of Wg-1 and Wg-2.

    2. NONLINEAR INTERFEROMETER INTEGRATING A PT-SYMMETRIC COUPLER

    A. Design of Interferometer

    We consider an integrated nonlinear interferometer based on waveguides (Wgs) as shown in Fig. 1(a). It consists of a pair of evanescently coupled Wgs where Wg-2 interacts with an analyte that can introduce losses to form a passive PT coupler. Two short sections of Wg-1, shown in yellow and denoted by NL1 and NL2, possess a quadratic nonlinear susceptibility and function as SPDC sources of signal–idler photon pairs when pumped by a continuous wave laser; here shown as the pump. On the other hand, the section of Wg-1 between NL1 and NL2 acts as a purely linear element. In practice, such a configuration of linear and nonlinear elements can be easily realized within the same physical Wg made out of quadratic nonlinear materials by engineering the spatial distribution of their effective nonlinearity through quasi-phase matching (QPM) [35]. For instance, sections NL1 and NL2 could be periodically poled to modulate their nonlinearity at period Λ, while the portion between these could be left unpoled. Assisted by the additional grating vector, K=2π/Λ, arising from the QPM poling, the SPDC process can occur efficiently in the two poled sections while the central section merely acts as a dispersive linear element.

    Whereas a simple device geometry employing two spatially separated sources of photon pairs in a single Wg (here Wg-1) can already function as a nonlinear interferometer [20], the key suggestion of this work is to introduce a second linear waveguide (Wg-2) adjacent to Wg-1. We show that this enriches the response of the nonlinear interferometer by making use of the exceptional dispersive properties of the resulting PT coupler.

    Throughout this paper, we consider the operation of SPDC in the nonlinear sections in the nondegenerate wavelength regime such that the idler photon wavelength is λi>2  μm while the pump and signal wavelengths lie in the λp,s<1  μm range [10,36,37]. Although the exact wavelengths of operation can be suitably chosen as needed by the application, the ranges we indicate here are meant to convey the intended utility of the scheme in facilitating spectroscopic applications in the technologically important mid-infrared spectral range with detection and excitation wavelengths being in the visible or near-infrared range.

    A salient feature of our proposed design is that the separation distance between Wg-1 and Wg-2 can be suitably chosen so that the longer wavelength idler photon generated in NL1 can tunnel to Wg-2 with coupling rate Ci and evanescently interact with the analyte deposited in its vicinity. At the same time, the shorter wavelength pump and signal photons remain confined to Wg-1 and do not interact with the analyte. However, the idler photon can tunnel back into Wg-1 ensuring that the biphoton amplitude emanating from NL1 can interfere with that from NL2. Due to the nonlocal nature of this quantum interference, the effective absorption strength in Wg-2, denoted as γi, can be determined by counting just the signal photons at the output of Wg-1. Note that neither the pump nor the signal photons ever propagate through Wg-2 in this scheme.

    B. PT Coupler

    Before we discuss the generation and evolution of photon pairs in the complete interferometer structure shown in Fig. 1(a), we provide a brief description of the evolution of classical optical fields within the linear PT coupler between the photon-pair sources. Since the pump and signal fields remain confined to Wg-1, the evolution of their field amplitudes, ap,s(z), is described by ap,s(z)/z=iβp,s·ap,s(z) with βp,s being their respective propagation constants. On the other hand, the idler field can couple to Wg-2 and, therefore, its field amplitude is described by a two-dimensional vector, Ai(z)=[ai(1)(z)ai(2)(z)]T whose evolution is given by Ai(z)z=MiAi(z),with Mi=[i(β¯i+Δβi)iCiiCii(β¯iΔβi)γi].

    Here, β¯i(βi(1)+βi(2))/2 and Δβi(βi(1)βi(2))/2 are defined in terms of the idler propagation constants βi(1,2) of the two uncoupled Wgs. The loss coefficient γi is responsible for idler loss in the coupler and is present exclusively in Wg-2. Finally, Ci denotes the coupling constant between Wg-1 and Wg-2. Solutions of Eq. (1) can be obtained analytically by making use of the eigenvalues λ±=iβ¯i+λ˜± and corresponding eigenvectors V± of the coupler propagation matrix Mi, which are λ˜±=±ξiγi2,V+=[γi+ξi2iΔβi2iCi]T,V=[2iCiγi+ξi2iΔβi]T,where ξi=(γi2iΔβi)2(2Ci)2 and Ai(z)=α+eλ+zV++αeλzV.

    The coefficients α± are determined by the input fields to the coupler. The eigenvectors V± describe the two supermodes of the coupler in terms of the individual modes of the two Wgs. The corresponding eigenvalues define their complex propagation constants, which depend on the idler loss magnitude γi in a nontrivial manner as is evident from the expression for λ˜±.

    We display the characteristic dependence of Re(λ˜±) and Im(λ˜±), which describe the loss and wavenumber of the supermodes for a symmetric coupler (Δβi=0), in Figs. 1(b) and 1(c). The effect of PT symmetry breaking [22] occurs at γi/Ci=2. For idler losses below this threshold value, the two supermodes of the PT coupler have the same effective losses equal to γi/2 [see Fig. 1(b)]. This symmetry between supermodes V± is due to the nature of their composition in terms of individual modes of the two Wgs below the threshold loss, with V±(1)=V±(2). Above the PT symmetry breaking point, V+ starts to localize in Wg-1 with V+(1)>V+(2) while V localizes in Wg-2. Consequently, the loss coefficient for the V+ supermode decreases with increasing γi, while the opposite is true for the V supermode. As we demonstrate in the following, this has a characteristically defining effect on the behavior of the proposed nonlinear interferometer. Another important feature to note here is the dependence of Im(λ˜±) on γi which is shown in Fig. 1(c). Both supermodes of the coupler have loss-dependent wavenumbers only below the PT symmetry breaking point. In the broken symmetry regime, they become independent of γi.

    3. METHODS AND RESULTS

    After reviewing the linear properties of the PT coupler, we now show its profound influence on photon-pair generation in a nonlinear interferometer and discuss applications of the ensuing interference phenomena for quantum sensing. To mathematically describe the behavior of the proposed nonlinear interferometer, we draw upon the theoretical formulation developed in Refs. [38,39] for describing the process of photon-pair generation through SPDC and its propagation in lossy media. We simplify the analysis of pair generation in the nonlinear sections NL1 and NL2 by assuming that their lengths l are much smaller than the coupling length for the idler photon (lπ/2Ci). Hence, in the following we neglect the effect of coupling between Wg-1 and Wg-2 when describing the evolution of signal–idler biphoton amplitudes in the nonlinear sections. We must note here that our assumption about small nonlinear section lengths is not fundamental to the effects we demonstrate in this work. In fact, for practical implementation, one can choose a coupler design where Wg-1 bends away at the beginning and end of the structure to decouple the waveguides and nonlinear sections of arbitrary lengths can then be placed in these uncoupled regions in Wg-1. In such a design, pair generation in the first and second nonlinear sections will be unaffected by the loss in the coupler section.

    A. Biphoton State Amplitude Evolution

    We calculate biphoton amplitudes for the quantum state in the Wg number basis. This is primarily because this basis provides a clear understanding of the spatial dynamics of the generation and propagation of the photon pairs. The state is described by three probability amplitudes [38] corresponding to the three possible localizations of the signal and idler photons: ϕ11: signal–idler pair is jointly present in Wg-1,ϕ12: signal is present in Wg-1 while idler is in Wg-2,ϕ12(s): signal is in Wg-1 and idler photon is absorbed in Wg-2.

    In this section, we first focus on the evolution of the pair amplitudes ϕ=[ϕ11ϕ12]T; the evolution of ϕ12(s) is discussed in the subsequent section. The Schrödinger-type equations for evolution of ϕ11(z) and ϕ12(z) in each of the nonlinear sections are given by [40]ϕ1,ni(z)z=i(βs+βi(ni))ϕ1,ni(z)+δ1,niAp(z),where the subscript ni refers to the number of the Wg, 1 or 2, in which the idler photon could be present while its partner signal photon is confined to Wg-1, and the Kronecker delta function δ1,ni represents the fact that pair generation can occur only in Wg-1 where the classical pump field given by Ap(z)=Ap(0)×cos(Kz)×exp(iβpz) is present. Here, Ap(0) subsumes the strength of the quadratic nonlinearity and the factor cos(Kz) represents the dominant spatial modulation of the nonlinearity that realizes QPM in the two nonlinear sections with K=2π/Λ. Poling period Λ is chosen such that β¯pβ¯sβ¯i(1)K=0 for a specific set of pump, signal, and idler modes. By neglecting the coupling between the two waveguides in Eq. (5), we implicitly also have neglected the effects of idler loss in Wg-2 on the biphoton evolution in the nonlinear sections. The solution of Eq. (5) for the first nonlinear section NL1 with 0zl gives the photon-pair amplitude originating from the first source. At z=l this is ϕ(l)=[10]Tϕ0(l)with ϕ0(l)=[Ap(0)l/2]·ei(βs+βi(1)+ΔβNL2)lsinc(ΔβNLl2),where ΔβNL=βpβsβi(1)K denotes the phase mismatch for SPDC in NL1 and NL2.

    Within the PT coupler section of length L, the evolution of amplitude ϕ(z)=[ϕ11(z)ϕ12(z)]T is governed by ϕ(z)z=(iβs1+Mi)ϕ(z),where we have used the propagation matrix Mi from Eq. (1). We propagate the photon-pair amplitudes generated in NL1 through the coupler section of the interferometer with lzl+L by solving Eq. (8) making use of the eigenvalues and eigenvectors of matrix (iβs1+Mi). At z=l+L, this results in ϕ(l+L)=ei(βs+β¯i)LeγiL/2×[cosh2ΘeξiL/2sinh2ΘeξiL/2isinhΘcoshΘ(eξiL/2eξiL/2)]ϕ0(l),where sinhΘ(2Ci)/(γi+ξi2iΔβi)2(2Ci)2 and cosh2Θsinh2Θ=1 (see Appendix A for details of derivation). Finally, to obtain the complete biphoton state amplitude at the end of the interferometer, z=2l+L, we superpose the photon-pair amplitudes in Eq. (9) with the pair generation amplitude arising from the second source, NL2. This results in the amplitude ϕ(2l+L)=[eiβi(1)l00eiβi(2)l]eiβslϕ(l+L)+[10]ϕ0eiβp(l+L),where we note that the second nonlinear source can contribute only to ϕ11 since both photons are generated in and remain confined to Wg-1 due to this source.

    B. Measurable Signal Photon Intensity

    Quantum spectroscopy based on nonlinear interference and induced coherence requires only detection of the signal photon intensity. In the implementation discussed here, depicted in Fig. 1(a), the intensity measurement is performed only on the signal mode in Wg-1, as the signal photon generated in either source will not couple to Wg-2. The detection of a signal photon in Wg-1 could result from one of three possibilities for the state of the idler photon shown schematically in Fig. 2(a), corresponding to the three probability amplitudes ϕ11, ϕ12, and ϕ12(s). These possibilities are distinguishable, and hence the total signal intensity is a sum of the intensities of the three contributions as I1s=I11+I12+I12(s).

    (a) Top panel shows waveguides in which signal and idler photons are present for the biphoton (I11 & I12) and single photon I12(s) intensity contributions (see text for their description). (b) For lossless PT coupler, biphoton amplitudes from sources NL1 and NL2 interfere destructively. (c) At the critical loss, γi=γcr, pair generation in NL1 and NL2 becomes totally distinguishable and their amplitudes add up incoherently. (d) In the broken PT regime, the nonlinear interferometer exhibits increased indistinguishability and coherence between pair generation amplitudes from the two sources. In (b)–(d), a special coupler design with CiL=π is considered.

    Figure 2.(a) Top panel shows waveguides in which signal and idler photons are present for the biphoton (I11 & I12) and single photon I12(s) intensity contributions (see text for their description). (b) For lossless PT coupler, biphoton amplitudes from sources NL1 and NL2 interfere destructively. (c) At the critical loss, γi=γcr, pair generation in NL1 and NL2 becomes totally distinguishable and their amplitudes add up incoherently. (d) In the broken PT regime, the nonlinear interferometer exhibits increased indistinguishability and coherence between pair generation amplitudes from the two sources. In (b)–(d), a special coupler design with CiL=π is considered.

    Here, the intensity contributions stemming from the two-photon amplitudes can be directly calculated from Eq. (10) as I11(z)=|ϕ11(z)|2, I12(z)=|ϕ12(z)|2. Next, we derive the I12(s) contribution to the signal intensity. The probability amplitude for detecting the signal photon at position z with its partner idler photon being absorbed at position z is given by ϕ12(s)(z,z). Following Ref. [38], its evolution is given by ϕ12(s)(z,z)z={iβsϕ12(s)(z,z)for  zz0otherwise,with the initial condition ϕ12(s)(z,z)=i2γiϕ12(z). The signal intensity I12(s)(z) is accumulated over the length of the coupler from the absorbed idler amplitudes ϕ12(s)(z,z) and is equal to I12(s)(z)=lz|ϕ12(s)(z,z)|2dz=2γilz|ϕ12(z)|2dz,for  z>l,with I12(s)(l)=0. The second nonlinear source does not contribute to this amplitude, since it is assumed to be lossless and not coupled to Wg-2.

    Overall, we find that the second nonlinear section can only affect I11. Only the first source NL1 placed before the coupler can contribute to all three intensities I11, I12, and I12(s). Hence, the two sources only interfere, if NL1 contributed to I11 at the end of the structure, since both sources are indistinguishable in this intensity component. This means that induced coherence between the sources can only appear if there is a nonzero probability for the idler photon from NL1 to be located in Wg-1 in the region NL2. The dynamics of induced coherence is thus dependent on the state of the idler photon generated from NL1, which itself is controlled by the properties of the coupler, i.e., its coupling strength, length, and the loss in Wg-2.

    C. Constructive and Destructive Photon Interference and Its Dependence on Idler Loss

    To fully understand the dependence of induced coherence between the sources and the resulting interference in the signal intensity on the idler loss, we explicitly calculate the intensity stemming from the case where both photons are in Wg-1 using Eqs. (10) and (9), as this is the only contribution to the signal intensity that can show interference. It is I11(2l+L)={1+|Vγi|2+2Re[Vγi×ei(ΔβNL+K)(L+l)]}|ϕ0|2,with Vγi=eiΔβiLγiL/2[coshξiL2+(1+2sinh2Θ)sinhξiL2].

    In I11, coherence between pair generation amplitudes from the two sources is clearly manifested in the last term of Eq. (14), where the visibility of the resulting interference fringes as a function of ΔβNL is determined by parameter Vγi. The other contributions to the signal intensity can only result from photon-pair generation in the first nonlinear section NL1, and the explicit expressions for intensities I12(2l+L) and I12(s)(2l+L) do not show any signature of interference like I11(2l+L) (see Appendix B). As such, the magnitude of |Vγi| could be regarded as the measure of the indistinguishability of the two photon-pair sources in the interferometer enabling interference in the resulting signal intensity.

    We see from Eq. (14) that whether the interference is constructive or destructive depends on the additional phase acquired by the idler photon as it propagates through the PT coupler. This is encapsulated in Vγi, which is defined through the coupler asymmetry, the coupling strength, the length of the coupler, and the idler loss.

    For a symmetric coupler with equal real parts of the idler propagation constants in both waveguides, Δβi=0, the expression in Eq. (15) simplifies to Vγi=eγiL2(cosσiL2+γiσisinσiL2)  withσi=(2Ci)2γi2below the PT threshold loss, γi2Ci. Above PT threshold loss γi>2Ci, Vγi can be similarly shown to be Vγi=eγiL2(coshηiL2+γiηisinhηiL2)withηi=γi2(2Ci)2.

    We note from Eqs. (16) and (17) that Vγi is a real function and depends on idler loss γi. Further, only constructive interference is possible in the PT-broken regime, since Vγi>0 according to Eq. (17). Importantly, the sign of Vγi is unrestricted in the PT-symmetric regime. The transition between constructive and destructive interference can therefore occur at a critical idler loss γcr, when Vγcr=0.

    These features can be employed to tailor the interferometer response to loss variations as we discuss below.

    D. Evolution of Photon Intensities

    We show the characteristic evolution of the different intensity contributions along the interferometer in Figs. 2(b)–2(d) for three different values of loss in Wg-2. We consider the case of phase-matched SPDC with ΔβNL=0 and a symmetric PT coupler where Δβi=0. The relative phase between pair amplitudes from the two sources stems from the constant phase factor K(L+l) and the loss-dependent phase of Vγi. For simplicity, here we assume that the QPM grating vector is such that K(L+l)=2nπ with nN, resulting in constructive biphoton interference if no coupling was present between the two Wgs. However, since coupling is present, the phase of Vγi additionally modifies the interference. To maximize the photon-pair interference in the second nonlinear section, we choose the length of the PT coupler equal to twice the coupling length in absence of loss, corresponding to CiL=π. The nonlinear section length is taken to be Cil=0.05, much less than the coupling length. In practice, nonlinear section length l can be always suitably chosen to ensure that K(L+l)=2nπ if the coupler length L has to be kept constant. In the following numerical simulations, we keep the normalized value of the coupling constant fixed as Ci=1 with no loss of generality. Other parameters such as lengths, propagation constants, and losses are assumed to be appropriately scaled by Ci and rendered dimensionless. As an example for a coupler with the physical coupling of Ci=0.5  cm1 [40], the normalized length of L=π corresponds to a physical length L=π/Ci=2π  cm and the normalized idler loss of γi=2 for PT threshold loss corresponds to an actual loss of γi=2Ci=1  cm1.

    We present in Fig. 2(b) the Hermitian case in which the PT coupler is lossless with γi=0. In this case, the idler photon of the pair generated in NL1 completely couples to Wg-2 and back such that |ϕ11(l+L)|=|ϕ11(l)| while I12(l+L)=0. The photon-pair amplitude generated in NL2 is coherently added to the amplitude stemming from NL1, which in the considered case leads to complete destructive interference and a signal intensity of zero at the end of the interferometer. This is the manifestation of the fact that in this case the contributions of the two nonlinear sections to the photon-pair amplitude are indistinguishable, enabling interference. Indeed, for the considered parameters according to Eq. (16) the visibility parameter is Vγi=cos(CiL), leading to Vγi=1 for the chosen coupler length, such that the final biphoton intensity I11 becomes zero. Since no loss is present in Wg-2, the signal photon amplitude associated with the absorbed idler photon is also zero in this case (see I12(s)) and thus, the total final signal intensity I1s vanishes. The intensity values plotted in these figures are normalized by the maximum possible signal intensity which, as we discuss later, is attained for infinitely large idler loss in the considered coupler design.

    Introducing loss in Wg-2 for the idler photon changes the response of the interferometer. We show this in Fig. 2(c) for a loss γi=γcr=1.35 according to Eq. (18), which we refer to as the critical loss. At γcr, at the end of the PT coupler, the idler photon of a pair generated in NL1 is either present in Wg-2 or is absorbed in Wg-2 during propagation in the coupler, such that I11(l+L)=0. Therefore, in this case, no interference takes place between the two sources, because they are completely distinguishable and I11(2l+L) only contains contributions from the second source NL2. Consequently, the final output signal intensity I1s(2l+L) from the interferometer results from a completely incoherent addition of pair amplitudes originating from the two sources.

    We note that in standard implementations of quantum nonlinear interferometers [6], biphoton amplitudes from individual sources approach complete distinguishability, leading to their fully incoherent addition and the absence of interference, only at infinitely high losses. Due to the interplay of coupling and loss in our structure, complete distinguishability can be reached for the finite loss γcr. As we show in the following, the value of γcr can be controlled by appropriately choosing the coupling constant and length of the PT coupler. We also find that even several such critical loss values can exist for a particular structure, where for a fixed coupling strength the coupler length has to be made longer to increase the number of critical points corresponding to complete distinguishability (see Appendix B).

    The behavior of the interferometer is quite different for losses above the PT symmetry breaking threshold γth=2Ci, since the points of critical loss exist only in the PT-symmetric regime as discussed in Section 3.C. This aspect is reflected in Fig. 2(d), where we plot the various intensity contributions at a loss magnitude of γi=5. Here, the final signal intensity I1s(2l+L) as well as I11(2l+L) is larger compared to Fig. 2(c) despite stronger loss. This apparently counter-intuitive result can be explained by considering the propagation dynamics of the idler photon in the PT coupler for different losses. The nature of the two supermodes V± and their corresponding eigenvalues dictate the idler photon dynamics as expressed by Eq. (4). As discussed before, these eigenmodes of the PT coupler remain symmetric in the two Wgs below γth and exhibit equal losses as is shown in Fig. 1(b). For losses above γth, this symmetry is broken and the two Wgs start to effectively decouple from each other. The V+ eigenmode, with its corresponding eigenvalue shown by the red solid curve in Figs. 1(b) and 1(c) tends to be localized in the lossless Wg-1, while the other eigenmode localizes in the lossy Wg-2. As a result, the V+ mode experiences progressively lower effective losses while the V mode becomes increasing lossy in the broken PT regime. Since the first source NL1 is positioned in Wg-1 in the interferometer, its biphoton amplitude favorably excites the V+ idler supermode in the broken PT regime and hence, the idler may experience lower effective losses for γ>γth. As this supermode features a high probability for the idler photon to be localized in Wg-1, this means that also I11(l+L) is increased for γ>γth, which can interfere with the contribution from NL2. Hence, the indistinguishability of pair generation between the two sources is improved as the idler loss strength increases in Wg-2 in the broken PT regime.

    E. Signal Fringes and Sharp Transition at Critical Loss

    After discussing the signal intensity evolution within the structure for several specific values of idler loss and a fixed phase mismatch, we now explicitly demonstrate the sharp change in the behavior of the interferometer around the critical loss magnitude for a varying mismatch ΔβNL. In experimental realizations, this mismatch usually corresponds to the wavelength of the detected signal mode for a fixed pump wavelength.

    We show in Fig. 3(a) the dependence of the total signal intensity I1s on the idler loss strength γi and the phase mismatch parameter ΔβNL. We clearly observe the fringes as a function of ΔβNL. The phase mismatch is normalized by Ci, so the period of intensity fringes given by ΔβNL=2 here, corresponds to an actual phase shift of ΔβNLL=2π. For zero loss, the visibility of interference is perfect, being equal to unity, and the interference maxima appear for specific ΔβNL. As γi increases, the visibility is reduced, and it completely vanishes exactly at the critical loss γcr=1.35.

    Normalized total signal intensity I1s versus the idler absorption strength (γi) and phase mismatch parameter (ΔβNL). (a) In a symmetric coupler, βi(2)=βi(1), the signal intensity fringe position exhibits a sharp shift by half a period at the critical loss γi=γcr. (b) For the asymmetric coupler with Δβi=0.2, the signal fringe position shifts gradually with the absorption strength. For all the plots, CiL=π and Ci=1.

    Figure 3.Normalized total signal intensity I1s versus the idler absorption strength (γi) and phase mismatch parameter (ΔβNL). (a) In a symmetric coupler, βi(2)=βi(1), the signal intensity fringe position exhibits a sharp shift by half a period at the critical loss γi=γcr. (b) For the asymmetric coupler with Δβi=0.2, the signal fringe position shifts gradually with the absorption strength. For all the plots, CiL=π and Ci=1.

    Strikingly, the signal intensity fringe exhibits a sharp shift by half a period as the idler loss strength changes across this critical loss magnitude. This is illustrated by the blue dashed line overlaid on the signal intensity contour plot that tracks the position of the central intensity maximum. We also see that such behavior does not occur with increasing idler losses beyond the PT threshold loss γth=2. In fact, the signal intensity fringe acquires a specific position in ΔβNL and its visibility improves with increasing idler loss.

    This behavior of the signal intensity fringes and the associated visibility with increasing idler loss are due to the interference of photon-pair amplitudes from the two sources. For zero loss, γi=0, with the chosen PT coupler parameters, CiL=π, the idler photon probability amplitude from the first source, tunnels to Wg-2 and comes back completely to Wg-1 at the second nonlinear source. This leads to perfect indistinguishability between the sources and the resulting signal interference fringe shows maximum visibility. In contrast, at the critical idler loss, γi=γcr, the two sources are completely distinguishable since the idler photon, if generated in the first source, is either present in the Wg-2 at the end of the interferometer or is absorbed during propagation in the coupler. This leads to the vanishing of signal interference visibility. Further increasing the idler loss reduces the effective coupling between waveguides and increases the probability that the idler photon amplitude from the first source remains in Wg-1 when it reaches the end of the coupler section of the interferometer. This leads to the recovery of the interference visibility since the sources become partially indistinguishable again.

    We have, thus far, only visualized the case of the symmetric coupler where the two Wgs have equal propagation constants for the idler mode, i.e., βi(2)=βi(1) and thus Δβi=0. We now analyze the effect of breaking this symmetry when Δβi0. The resulting signal intensity behavior is displayed in Fig. 3(b) for a detuning magnitude of Δβi=0.2. Unlike the PT-symmetric case, the interferometer employing the asymmetric coupler shows a gradual shift in the position of the signal intensity fringe with a change in the idler loss strength. Also, the visibility of the interference does not reduce to zero for any specific loss value.

    F. Sensing of Idler Losses

    Based on the nontrivial dependence of the signal intensity on the idler loss, next we suggest a specific sensing application of our scheme. We have shown the dependence of the interference fringes on the loss in Fig. 3(a), and we explicitly show the associated variation of their visibility in Fig. 4(a) with the red curve. Here, the visibility V of the signal interference fringe for loss γi is calculated as V=[maxΔβNL(I1s)minΔβNL(I1s)]/[maxΔβNL(I1s)+minΔβNL(I1s)], where I1s(L+2l) is the signal intensity as a function of the phase mismatch ΔβNL at constant loss magnitude. We note that V is different from the visibility parameter Vγi discussed before for the I11 intensity due to the contributions of I12 and I12(s) to the total signal intensity.

    (a) Visibility V of the interference in the signal intensity at the end of the nonlinear interferometer containing the PT coupler versus the idler loss for different coupler lengths. (b) Corresponding normalized total signal intensity at ΔβNL=0. The dots denote the boundaries of regions where the visibility is below 0.05 and the intensity is close to 0.5. (c) Sensitivity |dVdγi| for visibility and (d) sensitivity |dI1sdγi| for narrowband signal intensity measurements.

    Figure 4.(a) Visibility V of the interference in the signal intensity at the end of the nonlinear interferometer containing the PT coupler versus the idler loss for different coupler lengths. (b) Corresponding normalized total signal intensity at ΔβNL=0. The dots denote the boundaries of regions where the visibility is below 0.05 and the intensity is close to 0.5. (c) Sensitivity |dVdγi| for visibility and (d) sensitivity |dI1sdγi| for narrowband signal intensity measurements.

    In Fig. 4(a), we mark a region in gray where the visibility V is below 0.05. This value can serve as a specific threshold magnitude of visibility below which interference fringes cannot be detected reliably. We have highlighted the corresponding two points on the red curve between which the visibility is below this threshold value. The critical loss γcr, for which the visibility vanishes completely, lies between these two points. In Fig. 4(b) we show the corresponding signal intensity I1s(ΔβNL=0) for the phase-matched case along with the two threshold visibility points. In the intermediate region between these points, we have I1s0.5×I1s(γi) [see Fig. 4(b)], signifying that the two sources of photon pairs are almost distinguishable. Hence, no notable interference pattern in the signal intensity can be observed.

    The range of loss magnitudes where the two photon-pair sources are nearly completely distinguishable and the visibility is very close to zero can be controlled by adapting the design parameters of the nonlinear interferometer. We show this by considering designs with couplers of lengths L=2π and 3π, shown in Figs. 4(a) and 4(b) with the blue and black curves, respectively. The visibility plots clearly show that a larger loss region with visibilities below 0.05 can be realized for longer coupler lengths with several critical loss points lying within this region. It is important to realize that all such critical points lie below the PT-symmetry-breaking threshold γth=2. Furthermore, for loss magnitudes within this region, the signal intensity stays almost constant, regardless of the precise value of the loss (see Appendix B). The variation of the signal intensity with changing idler loss measured through its slope |dI1s/dγi| can be regarded as a figure-of-merit to gauge the sensitivity of the measurement scheme for sensing idler loss when measurements are made at a given signal wavelength. Alternatively, for broadband signal intensity fringe measurements sensitivity can be defined through the fringe visibility in terms of its slope |dV/dγi| for various idler losses. We show the variation of these sensitivity parameters in Figs. 4(c) and 4(d). Both of these plots show similar overall behavior. The cases with L=2π and L=3π are characteristically different from L=π in that the sensitivities reduce to zero at certain idler losses below the PT threshold loss unlike the L=π case and become nonzero again for increasing losses. In fact, there is a range of loss within which both |dV/dγi| and |dI1s/dγi| have very small magnitudes.

    The peculiar features discussed above could be applicable to sensing scenarios that require determination of whether the concentration of the analyte responsible for the loss is within a specific target range. Importantly, the boundaries of the range with approximately constant signal intensity can be fixed by appropriately choosing the length and coupling strength of the PT coupler. For such target loss-range sensing it is important that the signal intensity changes appreciably with change in loss when the background loss goes outside the specified range. For the two design cases shown in Figs. 4(a) and 4(b) with L=2π and 3π, the slope of signal intensity is larger than with L=π for small losses. At the same time, the design with L=2π shows a larger variation of signal intensity and a higher fringe visibility for losses γi>2 compared to the L=3π case. Hence, it is important to choose an appropriate length of the coupler that can provide sufficiently large sensitivities outside of the target range. We find that for longer couplers the sensitivity at larger losses is reduced, while at very small losses the sensitivity is increased, as shown by the plots in Figs. 4(c) and 4(d). Importantly, our suggested device can be preferably used to sense very high absorption. Here, contrary to standard nonlinear interferometers, the visibility of interference is revived and increasing for larger losses, as shown by Fig. 4(a), which renders their detection easier.

    It is instructive to compare the proposed nonlinear interferometric sensing scheme with the conventional linear absorption sensing based on a PT coupler, in which laser light at the idler wavelength is transmitted through the coupler and detected at its end. As was shown in Ref. [22], the transmitted intensity in such a scheme shows a nontrivial dependence on the absorption loss such that it first decreases with increasing losses up to a certain loss value and then starts to increase with increasing losses. In our work, we employ the unconventional transmission property of the PT coupler for the idler photon. However, the most important difference of the presented scheme is that we utilize the spectral correlation of the signal–idler pairs generated through SPDC in the nonlinear sections of the interferometer to effectively transfer the information about the losses at the longer idler wavelength to the short-wavelength signal-photon spectrum. This makes it possible to employ highly sensitive and low-noise detectors to measure the signal intensity and infer the losses at the idler wavelength.

    G. Broadband Absorption Spectroscopy

    Spectroscopic applications employing nonlinear interferometers typically aim to achieve operation in a broad spectral range such that the frequency-dependent absorption profile of an analyte can be measured around the idler photon frequency by performing detection at the corresponding frequencies of the signal photon. Such operation is essentially achieved by employing broadband SPDC in the interferometer where two short nonlinear sources of photon pairs are pumped by a continuous-wave pump laser at a fixed frequency. We now perform a similar frequency-dependent loss analysis for the interferometric scheme that we propose in this work.

    To illustrate the spectral response, we start by considering the frequency dependence of the phase mismatch parameter ΔβNL due to the dispersion of signal and idler modes. Under a first-order approximation, the propagation constants of the signal and idler modes can be described as βs,i(1)(ω¯s,i+Δωs,i)=βs,i(1)(ω¯s,i)+(1/vg)s,i×(Δωs,i), where ω¯s,i are the respective central frequencies and vg denotes the group velocity at these frequencies. Assuming a continuous wave pump at frequency ωp=ω¯s+ω¯i, the phase mismatch parameter then becomes ΔβNL=Δ(1vg)×Δωsωp,with Δ(1vg)=[(1/vg)s(1/vg)i]×ωp being the normalized inverse group velocity difference. Here, we have made use of the SPDC phase matching relation and the energy conservation condition Δωs+Δωi=0.

    The bandwidth of the interferometer is limited by the bandwidth of the biphoton amplitude ϕ0 from SPDC sources and is given by ΔβNL=4π/l according to Eqs. (14) and (7). This determines the width of the envelope containing the signal intensity fringes. On the other hand, the periodicity of these fringes is determined by the length of the central linear section according to Eq. (14) and is given by ΔβNL=2π/(L+l). To detect change in the spectral position and nature of fringes in response to idler loss, densely packed spectral fringes are advantageous, and this regime is realized for a long coupler length L and larger group velocity mismatch Δ(1vg) according to Eq. (19). At the same time though, a large group velocity mismatch also reduces the SPDC bandwidth and therefore, judiciously chosen small lengths of the SPDC sources ensure that the overall fringe envelope is still broad enough to not alter the relative peak intensities of the signal fringes significantly near the center of the envelope.

    We demonstrate a typical frequency response of the proposed interferometer by considering a spectrally localized loss profile for the idler mode centered around ω¯i, as shown by the black solid curve in Fig. 5. Due to the energy conservation, Δωi=Δωs in the plot. The signal intensity fringes in presence of this loss are calculated by assuming inverse group velocity difference of Δ(1vg)=103, and are shown by the red solid curve. We also display the signal intensity for the lossless idler case by the green dotted curve, which serves as a reference.

    Signal spectral intensity (red solid curve) from the nonlinear interferometer in the presence of a spectrally localized idler absorption profile (black solid curve) with the corresponding frequency given by Δωi=−Δωs. Comparison with the reference case of lossless idler, shown by green dotted curve, reveals that the nature of interference, constructive or destructive, switches at the critical idler loss γi=γcr represented by the two vertical dashed lines.

    Figure 5.Signal spectral intensity (red solid curve) from the nonlinear interferometer in the presence of a spectrally localized idler absorption profile (black solid curve) with the corresponding frequency given by Δωi=Δωs. Comparison with the reference case of lossless idler, shown by green dotted curve, reveals that the nature of interference, constructive or destructive, switches at the critical idler loss γi=γcr represented by the two vertical dashed lines.

    The broadband nature of SPDC in the interferometer ensures that the envelope of the signal spectral fringes is much wider than the spectral width of the absorption profile under investigation. Thus, the peaks in the normalized signal intensity shown by the green dotted curve in Fig. 5 are close to unity over the width of the black solid curve representing the idler loss profile. Meanwhile, the period of the signal fringes is much smaller than the absorption width such that several fringes encompass the spectral range of loss.

    We see from the red curve in Fig. 5 that the nature of interference, either constructive or destructive, and accordingly the final signal intensity, depend on the strength of the idler loss at the corresponding idler frequency. To demonstrate this, we mark the central region of the idler loss profile by two vertical black dashed lines, where the magnitude of loss is greater than the critical loss γcr. In this frequency range, maxima (and minima) in signal intensity get reversed as compared to the reference signal intensity due to reversal of the nature of biphoton interference. On the other hand, for frequency ranges in which the idler loss is less than the critical loss, the spectral positions of the signal intensity maxima (and minima) follow that of the case with no loss. In both these cases, the magnitude of fringe contrast which determines the prominence of the peaks in signal intensity depends on the magnitude of the idler loss. In both the central spectral region of higher losses and the low-loss region outside it, fringe contrast is significant as long as the loss magnitude is appreciably different from the critical loss γcr.

    4. CONCLUSIONS

    In our work, we investigated sensing with an integrated nonlinear quantum interferometer consisting of a PT coupler between two identical sources of signal and idler photon pairs. We established that a peculiar structure of the eigenmodes in the PT coupler enables the tailoring of the quantum interference effects to facilitate the sensing of losses in the second waveguide exposed to an analyte under investigation. In particular, the first-order interference between the signal photons generated in the two nonlinear sources strongly depends on the idler loss in the coupled waveguide, which the signal photons never enter.

    We identify a new effect of a sharp shift of the signal interference fringes in the vicinity of critical loss magnitudes, a phenomenon that is strikingly different from the previously studied nonlinear interferometers or PT-symmetric structures. Such a peculiar dependence of the signal intensity on the idler loss can benefit the sensing applications. For instance, if the loss depends on the analyte concentration, the interferometer can be engineered such that the signal intensity is constant for a particular allowed range of concentrations and only changes its value when the concentration leaves this target corridor. Furthermore, broadband spectroscopic information about the analyte absorption can be obtained from the dependencies of the signal interference fringes on the wavelengths of the signal and idler photons, where a sudden half-period shift of fringes can directly indicate the change of absorption beyond a critical magnitude.

    We anticipate that these results will stimulate advances in the fundamentals and applications of quantum enhanced sensing with nonlinear interferometers. The proposed integrated scheme using nonlinear waveguides could be particularly beneficial for on-chip sensing applications in the detection of gases and biological analytes with characteristically strong absorption lines in the infrared wavelength ranges [4143]. This will complement the growing field of on-chip linear absorption spectroscopy applications by expanding the operation in the infrared spectral range.

    APPENDIX A: EVOLUTION OF BIPHOTON AMPLITUDES IN THE INTERFEROMETER

    Here we provide calculations for the biphoton state amplitudes. The nonlinear interferometer is composed of two waveguides (Wg-1 and Wg-2) that are selectively coupled at the longer wavelength idler radiation through coupling constant Ci, as shown in Fig. 1(a) of the main text. Throughout these calculations we assume that the classical pump radiation and signal photon remain confined to Wg-1 due to their shorter wavelength. Two short sections each of length l at either ends of Wg-1 act as nonlinear elements in the interferometer where pair generation could take place through SPDC of the pump photon. The nonlinearity of these sections could be thought of as being “turned on” by appropriate choice of QPM poling period. The central section of length L between these two nonlinear sections in Wg-1 is assumed to behave as a purely linear element. Further, Wg-2 also acts as a linear element.

    We describe evolution of biphoton amplitudes in the interferometer in stages. These are as follows.System.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElement

    The quantum state in the interferometer is defined by the probability amplitudes ϕ11(z), ϕ12(z), and ϕ12(s)(z,z) as described in Section 3.

    Photon-Pair Amplitude from the First Source

    In the nonlinear source section NL1 (0zl), ϕ11(z) evolves as ϕ11(z)z=i(βs+βi(1))ϕ11(z)+Ap(z),withϕ11(0)=0.

    Ap(z) denotes the classical pump field under the undepleted pump approximation and is given by Ap(z)=Ap(0)×cosKz×eiβpz=Ap(0)2(eiKz+eiKz)eiβpz=Ap(0)2[ei(βpK)z+ei(βp+K)z].

    Here Ap(0) subsumes the strength of quadratic nonlinearity while its effective spatial profile is captured by the factor cosKz. The grating vector K=2π/Λ reflects the QPM that facilitates SPDC in NL1.

    By substituting ϕ11(z)=ϕ˜11(z)ei(βs+βi(1))z in Eq. (A1), we get ϕ˜1,1(z)z=Ap(0)2[ei(βpK)z+ei(βp+K)z]ei(βs+βi(1))z=Ap(0)2[ei(βpβsβi(1)K)z+ei(βpβsβi(1)+K)z].

    The two terms in the above expression denote SPDC processes which are mediated by the +K and K Fourier components of the nonlinearity. Usually only the term corresponding to ΔβNL=βpβsβi(1)K can be phase-matched such that ΔβNL=0 for a specific set of pump, signal, and idler modes. Therefore, if we neglect the non-phase-matched term corresponding to ei(βpβsβi(1)+K)z, ϕ˜1,1(z)z=Ap(0)2eiΔβNLz.

    The solution of Eq. (A3) is ϕ˜1,1(z)ϕ˜1,1(0)=Ap(0)2eiΔβNLz1iΔβNL=Ap(0)×z2eiΔβNLz2sinc(ΔβNLz2),and using ϕ˜1,1(0)=0, ϕ1,1(z)=Ap(0)×z2ei(βs+βi(1)+ΔβNL2)zsinc(ΔβNLz2).

    Since waveguides are treated as effectively uncoupled in nonlinear section, 0zl, and the pump field is present only in Wg-1, ϕ12(z) and ϕ12(s)(z,z) amplitudes are zero in this section of the interferometer. At z=l, biphoton amplitudes ϕ11 and ϕ12 can be denoted as a column vector, ϕ(l)=[ϕ11(l)ϕ12(l)]=[10]ϕ0(l),with ϕ0(l) being the biphoton amplitude given by Eq. (7) in the main text.

    Photon-Pair Amplitude within the PT Coupler

    In the coupler section, lzl+L, the evolution of ϕ(z) is given by ϕ(z)z=(iβs1+Mi)ϕ(z),with iβs1+Mi=[iβsi(β¯i+Δβi)iCiiCiiβsi(β¯iΔβi)γi].

    The eigenvalues λ¯ and eigenvectors V of the propagation matrix (iβs1+Mi) are obtained by first solving the characteristic polynomial det(iβs1+Miλ¯1)=0 to get the two eigenvalues λ¯± and then solving the corresponding eigenvector equations (iβs1+Miλ¯±1)V±=0. Explicitly, λ¯±=i(βs+β¯i)+λ˜±,with λ˜± and V± given for the PT coupler in Eq. (3). Evolution of a general state ϕ(z) in the coupler section can be expressed using λ¯± and V± as ϕ(z)=α+eλ¯+(zl)V++αeλ¯(zl)V.

    The coefficients α± are calculated by considering the initial state ϕ(l) at the input of the coupler given in Eq. (A6). These are the solution of α+V++αV=[10]Tϕ0(l).

    A straightforward calculation gives α+=(γi+ξi2iΔβi)ϕ0(l)(γi+ξi2iΔβi)2+(2iCi)2,and α=(2iCi)ϕ0(l)(γi+ξi2iΔβi)2+(2iCi)2.

    Putting these back into Eq. (A8) and after some simplification we get ϕ(z)=ei(βs+β¯i)(zl)eγi(zl)/2×[cosh2Θeξi(zl)/2sinh2Θeξi(zl)/2isinhΘcoshΘ[eξi(zl)/2eξi(zl)/2]]ϕ0(l),where sinhΘ(2Ci)/(γi+ξi2iΔβi)2(2Ci)2 and cosh2Θsinh2Θ=1.

    Due to presence of loss in the coupler section, the amplitude ϕ12(s)(z,z), which describes the contribution from pairs with lost idler photon, becomes nonzero. Its evolution is given in Eq. (12) with the solution ϕ12(s)(z,z)=[i2γi×ϕ12(z)]eiβs(zl),lzl+L,where ϕ12(z) is given by Eq. (A11) in the coupler.

    Photon-Pair Amplitude in the Second Nonlinear Section

    The amplitude for photon-pair generation in NL2 is easily obtained by using Eq. (A6). It is given by ϕ(l)×eiβp(l+L). Here, the phase factor eiβp(l+L) accounts for the evolution of the pump amplitude in the first nonlinear source and the coupler section of Wg-1. Now, to obtain the final photon-pair amplitudes from the interferometer as a whole, we superpose this biphoton amplitude from NL2 with that stemming from the NL1, ϕ(l+L), which is given by Eq. (A11). The resulting expressions for photon-pair amplitudes are given by Eq. (10) in the main text.

    APPENDIX B: SIGNAL INTENSITY FROM THE INTERFEROMETER

    At the output of the interferometer, total signal intensity I1s(2l+L) can be calculated by summing up the three different contributions I11(2l+L)=|ϕ11(2l+L)|2, I12(2l+L)=|ϕ12(2l+L)|2, and I12(s)(2l+L)=ll+L|ϕ12(s)(z,z)|2dz=2γill+L|ϕ12(z)|2dz=2γill+LI12(z)dz. Using Eq. (10), we obtain I11(2l+L)={1+|Vγi|2+2Re[Vγi×ei(ΔβNL+K)(L+l)]}|ϕ0|2,with Vγi=eiΔβiLγiL/2[coshξiL2+(1+2sinh2Θ)sinhξiL2],and I12(2l+L)=eγiL|sinh2Θ×sinhξiL2|2|ϕ0(l)|2.

    When the coupler in the interferometer is symmetric (Δβi=0), ξi=γi2(2Ci)2. We see that ξi is a purely imaginary quantity below the PT threshold and a positive real quantity above the PT threshold. This in turn affects the nature of Vγi.

    Below PT Threshold

    Since γi<2Ci in the PT-symmetric regime, ξi can be written as ξi=iσi, where σi=(2Ci)2γi2>0. For substituting ξi by iσi in Eq. (B2), we note that coshiσiL2=cosσiL2 and sinhiσiL2=isinσiL2. Further, 1+2sinh2Θ=γiiσi. Putting these into Eq. (B2) leads to Eq. (16) for Vγi in the PT-symmetric regime.

    Above PT Threshold

    For the broken PT regime, we denote ηi=γi2(2Ci)2>0. It can be shown that 1+2sinh2Θ=γi/ηi above the PT threshold. Thus, using these definitions in Eq. (B2) readily leads to Eq. (17) for Vγi.

    We show the signal intensity I1s for the interferometer design presented in the main text (l=0.05, L=π, and Ci=1) as a function of phase mismatch ΔβNL and idler loss γi in Fig. 6(a). The three contributions I11, I12, and I12(s) are presented in Figs. 6(b), 6(c), and 6(d), respectively. We see that the I11 signal intensity component shows the interference fringes whereas I12 and I12(s) add a background which is almost independent of ΔβNL for any given loss magnitude γi. This means that these two components reduce the overall visibility of fringes in I1s. Similar to the case considered in the main text, we assume K(L+l)=2nπ for the relative phase between biphoton amplitudes from the two sources.

    (a) Signal intensity I1s as a function of phase mismatch ΔβNL and idler loss γi. (b), (c), and (d) show the three individual contributions to the total signal intensity. I12 and I12(s) are devoid of any interference fringes unlike I11.

    Figure 6.(a) Signal intensity I1s as a function of phase mismatch ΔβNL and idler loss γi. (b), (c), and (d) show the three individual contributions to the total signal intensity. I12 and I12(s) are devoid of any interference fringes unlike I11.

    Signal intensity fringes for different relative phases K(L+l) in the interferometer. The feature of fringe shift occurs in each case at the same idler loss value irrespective of K(L+l) phases.

    Figure 7.Signal intensity fringes for different relative phases K(L+l) in the interferometer. The feature of fringe shift occurs in each case at the same idler loss value irrespective of K(L+l) phases.

    Illustration of the condition for existence of critical idler loss, γcr is CiL>π/2, for different L values as indicated by labels. (a)–(d) When CiL≤π/2, no critical loss points are observed and hence no shift in the fringes. (e), (f) For CiL=π, a fringe shift is observed. (a), (c), (e) Fringe visibility; (b), (d), (f) signal intensity fringes. For all the plots Ci=1.

    Figure 8.Illustration of the condition for existence of critical idler loss, γcr is CiL>π/2, for different L values as indicated by labels. (a)–(d) When CiLπ/2, no critical loss points are observed and hence no shift in the fringes. (e), (f) For CiL=π, a fringe shift is observed. (a), (c), (e) Fringe visibility; (b), (d), (f) signal intensity fringes. For all the plots Ci=1.

    Signal intensity I1s for nonlinear interferometers with PT couplers of length (a) L=π, (b) L=3π, and (c) L=5π. (d) Variation of the corresponding visibility with loss γi. The grayed area shows region with visibility less than or equal to 0.05.

    Figure 9.Signal intensity I1s for nonlinear interferometers with PT couplers of length (a) L=π, (b) L=3π, and (c) L=5π. (d) Variation of the corresponding visibility with loss γi. The grayed area shows region with visibility less than or equal to 0.05.

    APPENDIX C: SINGLE WAVEGUIDE NONLINEAR INTERFEROMETER

    We have mentioned that a nonlinear interferometer can already be implemented using a single waveguide with two spatially separated nonlinear sections [20]. In such a device the long wavelength idler mode interacts with the surrounding analyte through its evanescent fields while the pump and the signal mode are tightly confined within the core of the waveguide [20]. This allows sensing of idler loss by measurement of only the signal intensity.

    To contrast the response of a single waveguide nonlinear interferometer with the one proposed in this work that incorporates a passive PT coupler, we show the variation of their signal intensities for phase-matched SPDC (ΔβNL=0) with idler loss in Fig. 10(a). For a fair comparison, we assume broadband sources of photon pairs in both the cases and take the length of the linear section separating the sources to be the same as the length of the coupler, which is equal to L=π. In Fig. 10(b) we show the dependence of the signal fringe visibility for the two cases. The remarkable feature of the PT coupler-based scheme is that unlike the single waveguide design, where the visibility decreases exponentially with idler loss (blue curve), here the visibility first decreases to zero with increasing idler loss and then revives again at larger losses (red curve). Such a dependence improves the sensitivity of the PT coupler-based scheme at larger idler losses in comparison to a single waveguide nonlinear interferometer. To illustrate the improvement in sensitivity, we show the slope of signal intensity |dI1s/dγi| and the slope of visibility |dV/dγi| in Figs. 10(c) and 10(d), respectively, for the two cases. It can be seen that the single waveguide design can provide slightly higher sensitivity at small loss magnitudes. However, the design with the PT coupler is much more sensitive at larger losses, where the exact amount of improvement depends on the background idler loss at which the sensor is operating.

    Comparison of the performance of the proposed interferometer (incorporating the PT coupler) with that of the conventional integrated nonlinear interferometer with two separated sources in a single waveguide (single Wg). (a) Signal intensity, (b) visibility, (c) slope of signal intensity, and (d) slope of visibility for the two cases.

    Figure 10.Comparison of the performance of the proposed interferometer (incorporating the PT coupler) with that of the conventional integrated nonlinear interferometer with two separated sources in a single waveguide (single Wg). (a) Signal intensity, (b) visibility, (c) slope of signal intensity, and (d) slope of visibility for the two cases.

    APPENDIX D: EFFECT OF ANALYTE INDUCED ASYMMETRY IN THE COUPLER

    We have shown in Section 3.E that asymmetry between waveguides of the coupler, Δβi0, modifies the response of the nonlinear interferometer. In particular, the effect of sharp fringe shift at critical loss γcr does not occur for the asymmetric couplers. Since the presence of analyte can induce changes in the real part of the propagation constant of the second waveguide and make the coupler asymmetric, we now consider the effect of this induced asymmetry on the resulting signal intensity fringes.

    In Fig. 11(a) we show the signal spectral intensity in the presence of spectrally localized idler absorption profile by the solid red curve. Here, the coupler is taken to be asymmetric with idler propagation constant mismatch Δβi=0.2. We see that the influence of idler loss is still evident in the form of reduced fringe visibility. However, if the change in the propagation constant of the second waveguide becomes too large such that the resulting mismatch Δβi is greater than the coupling constant Ci=1, then the effect of idler loss on signal fringe visibility is drastically reduced due to the reduced effective coupling between the waveguides. We show this in Fig. 11(b) for Δβi=2. We see that the influence of idler loss on signal spectral fringe visibility is greatly reduced in this case. Unlike the case for small mismatch, here the fringe visibility does not show strong dependence on the idler loss magnitude. In fact, in the limit of infinitely large Δβi, the system effectively behaves like a nonlinear interferometer with two spatially separated sources in the first waveguide alone and no idler loss. Therefore, for detection of the idler loss through fringe visibility measurements, the proposed scheme can tolerate some asymmetry between the waveguides as long as the magnitude of induced Δβi change is smaller than the coupling Ci between the waveguides.

    Signal intensity fringes in the presence of waveguide asymmetry in the coupler with (a) Δβi=0.2 and (b) Δβi=2. The coupling constant in both cases is Ci=1.

    Figure 11.Signal intensity fringes in the presence of waveguide asymmetry in the coupler with (a) Δβi=0.2 and (b) Δβi=2. The coupling constant in both cases is Ci=1.

    References

    [1] D. Klyshko. Photons and Nonlinear Optics(1988).

    [2] L. J. Wang, X. Y. Zou, L. Mandel. Induced coherence without induced emission. Phys. Rev. A, 44, 4614-4622(1991).

    [3] X. Y. Zou, L. J. Wang, L. Mandel. Induced coherence and indistinguishability in optical interference. Phys. Rev. Lett., 67, 318-321(1991).

    [4] H. M. Wiseman, K. Molmer. Induced coherence with and without induced emission. Phys. Lett. A, 270, 245-248(2000).

    [5] M. Lahiri, A. Hochrainer, R. Lapkiewicz, G. B. Lemos, A. Zeilinger. Partial polarization by quantum distinguishability. Phys. Rev. A, 95, 033816(2017).

    [6] M. V. Chekhova, Z. Y. Ou. Nonlinear interferometers in quantum optics. Adv. Opt. Photonics, 8, 104-155(2016).

    [7] Z. Y. Ou, X. Y. Li. Quantum su(1,1) interferometers: basic principles and applications. APL Photonics, 5, 080902(2020).

    [8] C. M. Caves. Reframing su(1,1) interferometry. Adv. Quantum Technol., 3, 1900138(2020).

    [9] A. Ferreri, M. Santandrea, M. Stefszky, K. H. Luo, H. Herrmann, C. Silberhorn, P. R. Sharapova. Spectrally multimode integrated su(1,1) interferometer. Quantum, 5, 461(2021).

    [10] D. A. Kalashnikov, A. V. Paterova, S. P. Kulik, L. A. Krivitsky. Infrared spectroscopy with visible light. Nat. Photonics, 10, 98-102(2016).

    [11] G. B. Lemos, V. Borish, G. D. Cole, S. Ramelow, R. Lapkiewicz, A. Zeilinger. Quantum imaging with undetected photons. Nature, 512, 409-412(2014).

    [12] A. V. Paterova, H. Z. Yang, C. W. An, D. A. Kalashnikov, L. A. Krivitsky. Tunable optical coherence tomography in the infrared range using visible photons. Quantum Sci. Technol., 3, 025008(2018).

    [13] A. V. Paterova, L. A. Krivitsky. Nonlinear interference in crystal superlattices. Light Sci. Appl., 9, 82(2020).

    [14] A. Paterova, S. Lung, D. A. Kalashnikov, L. A. Krivitsky. Nonlinear infrared spectroscopy free from spectral selection. Sci. Rep., 7, 42608(2017).

    [15] A. Valles, G. Jimenez, L. J. Salazar-Serrano, J. P. Torres. Optical sectioning in induced coherence tomography with frequency-entangled photons. Phys. Rev. A, 97, 023824(2018).

    [16] C. Lindner, S. Wolf, J. Kiessling, F. Kuhnemann. Fourier transform infrared spectroscopy with visible light. Opt. Express, 28, 4426-4432(2020).

    [17] M. Kutas, B. Haase, P. Bickert, F. Riexinger, D. Molter, G. von Freymann. Terahertz quantum sensing. Sci. Adv., 6, eaaz8065(2020).

    [18] M. Ravaro, E. Guillotel, M. Le Du, C. Manquest, X. Marcadet, S. Ducci, V. Berger, G. Leo. Nonlinear measurement of mid-infrared absorption in AlOx waveguides. Appl. Phys. Lett., 92, 151111(2008).

    [19] A. S. Solntsev, P. Kumar, T. Pertsch, A. A. Sukhorukov, F. Setzpfandt. LiNbO3 waveguides for integrated SPDC spectroscopy. APL Photonics, 3, 021301(2018).

    [20] P. Kumar, S. Saravi, T. Pertsch, F. Setzpfandt. Integrated induced-coherence spectroscopy in a single nonlinear waveguide. Phys. Rev. A, 101, 053860(2020).

    [21] T. Ono, G. F. Sinclair, D. Bonneau, M. G. Thompson, J. C. F. Matthews, J. G. Rarity. Observation of nonlinear interference on a silicon photonic chip. Opt. Lett., 44, 1277-1280(2019).

    [22] A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. Volatier-Ravat, V. Aimez, G. A. Siviloglou, D. N. Christodoulides. Observation of PT-symmetry breaking in complex optical potentials. Phys. Rev. Lett., 103, 093902(2009).

    [23] C. E. Ruter, K. G. Makris, R. El-Ganainy, D. N. Christodoulides, M. Segev, D. Kip. Observation of parity-time symmetry in optics. Nat. Phys., 6, 192-195(2010).

    [24] R. El-Ganainy, K. G. Makris, D. N. Christodoulides, Z. H. Musslimani. Theory of coupled optical PT-symmetric structures. Opt. Lett., 32, 2632-2634(2007).

    [25] K. G. Makris, R. El-Ganainy, D. N. Christodoulides, Z. H. Musslimani. Beam dynamics in PT symmetric optical lattices. Phys. Rev. Lett., 100, 103904(2008).

    [26] R. El-Ganainy, K. G. Makris, M. Khajavikhan, Z. H. Musslimani, S. Rotter, D. N. Christodoulides. Non-Hermitian physics and PT symmetry. Nat. Phys., 14, 11-19(2018).

    [27] Z. Lin, H. Ramezani, T. Eichelkraut, T. Kottos, H. Cao, D. N. Christodoulides. Unidirectional invisibility induced by PT-symmetric periodic structures. Phys. Rev. Lett., 106, 213901(2011).

    [28] H. Hodaei, M. A. Miri, M. Heinrich, D. N. Christodoulides, M. Khajavikhan. Parity-time-symmetric microring lasers. Science, 346, 975-978(2014).

    [29] J. Wiersig. Enhancing the sensitivity of frequency and energy splitting detection by using exceptional points: application to microcavity sensors for single-particle detection. Phys. Rev. Lett., 112, 203901(2014).

    [30] H. Hodaei, A. U. Hassan, S. Wittek, H. Garcia-Gracia, R. El-Ganainy, D. N. Christodoulides, M. Khajavikhan. Enhanced sensitivity at higher-order exceptional points. Nature, 548, 187-200(2017).

    [31] M. Ornigotti, A. Szameit. Quasi PT-symmetry in passive photonic lattices. J. Opt., 16, 065501(2014).

    [32] P. B. Main, P. J. Mosley, A. V. Gorbach. Spontaneous parametric down-conversion in asymmetric couplers: photon purity enhancement and intrinsic spectral filtering. Phys. Rev. A, 100, 053815(2019).

    [33] J. Su, L. Cui, J. M. Li, Y. H. Liu, X. Y. Li, Z. Y. Ou. Versatile and precise quantum state engineering by using nonlinear interferometers. Opt. Express, 27, 20479-20492(2019).

    [34] J. M. Li, J. Su, L. Cui, T. Q. Xie, Z. Y. Ou, X. Y. Li. Generation of pure-state single photons with high heralding efficiency by using a three-stage nonlinear interferometer. Appl. Phys. Lett., 116, 204002(2020).

    [35] D. S. Hum, M. M. Fejer. Quasi-phasematching. C. R. Phys., 8, 180-198(2007).

    [36] M. Kumar, P. Kumar, A. Vega, M. A. Weissflog, T. Pertsch, F. Setzpfandt. Mid-infrared photon pair generation in AgGaS2. Appl. Phys. Lett., 119, 244001(2021).

    [37] Y. M. Sua, H. Fan, A. Shahverdi, J. Y. Chen, Y. P. Huang. Direct generation and detection of quantum correlated photons with 3.2 μm wavelength spacing. Sci. Rep., 7, 17494(2017).

    [38] D. A. Antonosyan, A. S. Solntsev, A. A. Sukhorukov. Effect of loss on photon-pair generation in nonlinear waveguide arrays. Phys. Rev. A, 90, 043845(2014).

    [39] D. A. Antonosyan, A. S. Solntsev, A. A. Sukhorukov. Photon-pair generation in a quadratically nonlinear parity-time symmetric coupler. Photonics Res., 6, A6-A9(2018).

    [40] A. Belsley, T. Pertsch, F. Setzpfandt. Generating path entangled states in waveguide systems with second-order nonlinearity. Opt. Express, 28, 28792-28809(2020).

    [41] P. T. Lin, V. Singh, J. J. Hu, K. Richardson, J. D. Musgraves, I. Luzinov, J. Hensley, L. C. Kimerling, A. Agarwal. Chip-scale mid-infrared chemical sensors using air-clad pedestal silicon waveguides. Lab Chip, 13, 2161-2166(2013).

    [42] M. C. Estevez, M. Alvarez, L. M. Lechuga. Integrated optical devices for lab-on-a-chip biosensing applications. Laser Photonics Rev., 6, 463-487(2012).

    [43] E. Ryckeboer, R. Bockstaele, M. Vanslembrouck, R. Baets. Glucose sensing by waveguide-based absorption spectroscopy on a silicon chip. Biomed. Opt. Express, 5, 1636-1648(2014).

    Pawan Kumar, Sina Saravi, Thomas Pertsch, Frank Setzpfandt, Andrey A. Sukhorukov. Nonlinear quantum spectroscopy with parity–time-symmetric integrated circuits[J]. Photonics Research, 2022, 10(7): 1763
    Download Citation