SYSTEM AND METHOD FOR DISPERSION-ENABLED QUANTUM STATE CONTROL OF PHOTONS

Information

  • Patent Application
  • 20170371105
  • Publication Number
    20170371105
  • Date Filed
    June 27, 2017
    7 years ago
  • Date Published
    December 28, 2017
    6 years ago
Abstract
Devices and methods are described for selecting a level of entanglement between two nondegenerate photons. The method may include receiving two non degenerate photons through a single input port of a directional photonic coupler; adjusting one of a first-order coupler dispersion M or a power splitting ratio η (λ00) of the directional optical coupler to select a Δη; and, emitting the photons from corresponding output ports of the directional optical coupler, wherein the emitted photons have a spectral entanglement corresponding to the selected Δη.
Description
FIELD OF THE INVENTION

The present invention generally pertains to the field of photonic couplers, and particular embodiments or aspects relate to photonic directional couplers adapted for dispersion-enabled quantum state control of photons.


BACKGROUND

Integrated optics has brought unprecedented levels of stability and performance to quantum photonic circuits. However, integrated devices are not merely micron-scale equivalents of their bulk-optics counterparts.


The quantum properties of light can unlock a variety of enhanced and novel technological capabilities. Among these are secure communications, nonclassical simulation, nonlocal imaging, and pathway-selective exciton spectroscopy. Such quantum photonic technologies have traditionally been implemented on the bench top with discrete optical components. More recently, the need for improved scalability has fuelled widespread interest in the development of on-chip quantum circuits. Much of this work has concentrated on the generation, manipulation, and detection of entangled photon pairs, often with the goal of replicating tasks previously performed using bulk optics.


There may be a need, however, for a system and method for quantum photonics that is not subject to one or more limitations of the prior art. In some embodiments, there is a need for an integrated optical component that may support quantum photonics in place of conventional bulk optics.


This background information is intended to provide information that may be of possible relevance to the present invention. No admission is necessarily intended, nor should be construed, that any of the preceding information constitutes prior art against the present invention.


SUMMARY

Integrated optical components can exhibit highly wavelength-dependent (i.e., dispersive) behavior compared to their bulk optics counterparts. Such dispersion has been shown to provide unprecedented tailorability over the properties of two-photon states generated by engineered nonlinear interactions in integrated waveguides.


By exploiting the inherently dispersive characteristics of the integrated optical component setting, such devices can play a remarkably more versatile role in quantum circuit architectures. In an embodiment, an ordinary photonic directional coupler may be adapted to support linear dispersion of photons. Dispersion unlocks several novel capabilities for the adapted photonic directional coupler including in situ control over photon spectral and polarization entanglement, tunable photon time ordering, and entanglement-sensitive two-photon coincidence generation.


In an embodiment, perfect two-photon anti-coalescence may be maintained while tuning the interference visibility.


In an embodiment, a method is provided for selecting a level of entanglement between two nondegenerate photons. The method may include receiving two non degenerate photons through a single input port of a directional photonic coupler; adjusting one of a first-order coupler dispersion M or a power splitting ratio η(λ00) of the directional optical coupler to select a Δη; and, emitting the photons from corresponding output ports of the directional optical coupler, wherein the emitted photons have a spectral entanglement corresponding to the selected Δη. In an implementation, the adjusting may comprise adjusting a waveguide core-cladding index of the directional optical coupler. In an implementation, the adjusting may comprise adjusting a separation of waveguides of the directional optical coupler. In an implementation, the method may further comprise selecting a level of time-ordering between the two nondegenerate photons by, applying a time delay τ to one of the output ports, wherein the time delay converts the selected Δη path-entanglement characteristics into time-ordering characteristics of the emitted photons. In an implementation, the method may further comprise selecting a level of interference visibility at a fixed photon pair anti-bunching (separation) rate for the two nondegenerate photons, wherein the receiving two non degenerate photons further comprises receiving two nondegenerate photons that are path-entangled across two input ports of the directional optical coupler with a relative phase of θ=π between paths; and, wherein the adjusting comprises: for a given M, adjusting η(λ00) of the directional optical coupler to select a desired interference visibility.


In an embodiment, a device is provided for selecting a level of entanglement between two nondegenerate photons. The device may include: a directional optical coupler adapted to receive two non degenerate photons through a single input port; means for adjusting a power splitting ratio η(λ) of the directional optical coupler to a selected Δη; wherein emitted photons from the directional optical coupler have a spectral entanglement corresponding to the selected Δη. In an implementation, the means for adjusting the power splitting ratio η(λ) comprises means for adjusting a waveguide core-cladding index of the directional optical coupler. In an implementation, the means for adjusting the power splitting ratio η(λ) comprises means for adjusting a separation of waveguides of the directional optical coupler. In an implementation, the means for adjusting the power splitting ratio η(λ) comprises a coupling length of the directional optical coupler longer than a minimum necessary value (L=π/(4κ(λ00))). In an implementation, the device is further operative to select a level of time-ordering between the two nondegenerate photons, wherein the device further comprises an extended path in one leg of the directional optical coupler, the extended path applying a time delay τ to one of the output ports, wherein the time delay τ converts the selected Δη path-entanglement characteristics into time-ordering characteristics of the emitted photons. In an implementation, the device is further operative to select a level of interference visibility at a fixed photon pair anti-bunching (separation) rate for the two nondegenerate photons, wherein the directional optical coupler is adapted to receive the two non degenerate photons as path-entangled photons across two input ports of the directional optical coupler with a relative phase of θ=π between paths; and, for a given M, the means for adjusting a power splitting ratio η(λ) of the directional optical coupler comprises means for adjusting η(λ00) of the directional optical coupler to select a desired interference visibility.


In an embodiment, a method is provided for selecting a level of entanglement between two nondegenerate photons comprising: receiving two non degenerate photons through a single input port of a directional coupler; adjusting M or η(λ00) of the directional coupler to select a desired Δη; emitting the photons from corresponding output ports of the directional coupler, wherein the emitted photons have a spectral entanglement corresponding to the selected Δη.


In an embodiment, a method is provided for selecting a level of time-ordering between two nondegenerate photons comprising: receiving two non degenerate photons through a single input port of a directional coupler; adjusting M or η(λ00) of the directional coupler to select a desired Δη; emitting the photons from corresponding output ports of the directional coupler, wherein the emitted photons have path-entanglement characteristics corresponding to the selected Δη; and, applying a time delay τ to one of the output paths, wherein this time delay converts the path-entanglement characteristics selected by Δη into time-ordering characteristics.


In an embodiment, a method is provided for selecting a level of interference visibility at a fixed photon pair anti-bunching (separation) rate for two nondegenerate photons comprising: receiving two nondegenerate photons that are path-entangled across two input ports of a directional coupler with a relative phase of θ=π between paths; for a given M, adjusting η(λ00) of the directional coupler to select a desired interference visibility.


In an embodiment, a method is provided for extracting the spectral entanglement of two photons of known nondegeneracy (Λ) and bandwidth (Δλ) comprising: receiving two photons that are path-entangled across two input ports of a directional coupler; measuring photon counts at the directional coupler output ports, including a rate at which the two photons exit from different ports; and, calculating the spectral entanglement of the photons using the measured count rates, and the values of M, η(λ00), θ, Δλ, and Λ.


Embodiments have been described above in conjunctions with aspects of the present invention upon which they can be implemented. Those skilled in the art will appreciate that embodiments may be implemented in conjunction with the aspect with which they are described, but may also be implemented with other embodiments of that aspect. When embodiments are mutually exclusive, or are otherwise incompatible with each other, it will be apparent to those skilled in the art. Some embodiments may be described in relation to one aspect, but may also be applicable to other aspects, as will be apparent to those of skill in the art.





BRIEF DESCRIPTION OF THE FIGURES

Further features and advantages of the present invention will become apparent from the following detailed description, taken in combination with the appended drawings, in which:



FIG. 1a illustrates an embodiment of a two-port directional coupler.



FIG. 1b illustrates a map of possible coupler responses to a two-photon input state, according to an embodiment



FIG. 2a illustrates an embodiment of a coupler.



FIG. 2b illustrates coupling strength for the coupler of FIG. 2a, according to an embodiment.



FIGS. 3a-3c, illustrate the dependence of two-photon path correlations on coupler response, according to an embodiment.



FIGS. 4a-4f illustrate the dependence of two-photon path correlations on coupler response, according to an embodiment.



FIGS. 5a-b illustrate the tunability of output state entanglement, according to an embodiment.



FIG. 6 is a plot illustrating the dependence of PS on entanglement, according to an embodiment.



FIG. 7 illustrates an embodiment of an all-integrated SN measurement.



FIG. 8 illustrates an embodiment of a multipurpose dispersive coupler for state characterization.



FIGS. 9 & 10 illustrate useful equations.





DETAILED DESCRIPTION

In some embodiments, new ways of leveraging dispersion in photonic directional couplers for the manipulation of two-photon states and their correlation properties are presented.


Photonic directional couplers (also referred to as “couplers” or “directional couplers”) are a common building block of integrated quantum circuits whose dispersion properties have yet to be fully exploited. They are typically implemented through the evanescent coupling of two identical waveguides and are characterized by a power splitting ratio η(λ). Their primary role has been to serve as on-chip beam splitters, often to mediate quantum interference. Due to the presence of dispersion in η(λ), these same couplers can also act as a wavelength demultiplexer (WD) for specific sets of nondegenerate wavelengths, without relying on waveguide modal mismatch. In fact, dispersion can cause the coupler's behavior to transition between “ideal” beam splitter operation and “ideal” WD operation in response to either the properties of the quantum state or systematic shifts to the coupling strength. The implications this has for two-photon state manipulation has yet to be studied. We show that this attribute of photonic directional couplers grants them a versatile set of new functionalities, which includes the post-selective tuning of spectral entanglement, entanglement-sensitive coincidence detection, and the ability to maintain perfect anti-coalescence while allowing full tunability over the two-photon interference visibility.


In some embodiments, symmetric 2×2 directional couplers are described as an example of quantum state engineering in integrated photonic systems without the loss of generality. As such, an essential step is to parameterize the photonic directional coupler's response for the two-photon state in terms of generic dimensionless variables that can be mapped to any combination of coupler and state properties. The details of this parametrization are described below, but some key definitions are introduced here.



FIG. 1a illustrates a generic two-port directional coupler, shown with simple implementations of thermal and/or electro-optic tuning for in situ control over η(λ00) on each of the coupler paths. FIG. 1b illustrates an illustrative map of possible coupler responses to a two-photon input state, as characterized by Δη. The coordinates labelled BS denote 50:50 beamsplitter behavior, while WD denotes perfect demultiplexing of central wavelengths λ01 and λ02.


Suppose two single-mode waveguides are coupled over a length L, such as in FIG. 1a. For symmetric rectangular waveguides, this leads to a power splitting ratio of η(λ)=cos2(κ(λ)L), where κ(λ) is the coupling strength. The wavelength dependence of the coupler is usually sufficiently described by its first-order coupler dispersion M=d(λ)L/dλ at a reference wavelength λ00, together with the value of η(λ00). Let λ01 and λ02 be the central wavelengths of a photon pair that evolves through this coupler, with Λ=|λ02−λ01| giving the nondegeneracy. Defining Δη=|(λ02)−η(λ01)| allows the coupler response to be classified as beam splitter-like for Δη→0 or WD-like for Δη→1.


The space of all possible Δη is spanned by (λ00) and the dimensionless product MΛ. This has been plotted in FIG. 1b, assuming λ00=(λ0102)/2 and negligible higher-order coupler dispersion. This plot provides a useful guide for relating the results described in the present application to transitions between beam splitter and WD behavior. A special condition, η(λ01)+η(λ02)=1, occurs along the lines η(λ00)=0 and MΛ=π/2, and corresponds to the splitting ratios η(λ01) and η(λ02) being antisymmetric about the 50:50 splitting value η=0.5.


Parameterization of Coupler Response

The analysis presented in this application assumes that the waveguides are single-mode, identical, and rectangular (i.e. non-tapered) as illustrated in FIG. 1a, although more sophisticated design geometries are possible. Spatial mode overlap between the waveguides leads to η(λ)=cos2 (κ(λ)L) in terms of a coupling strength κ(λ) over an interaction length L. As defined, η(λ) represents the probability that a photon exits from the same waveguide it enters from (i.e. η(λ)=1 means no power is transferred).


It is useful to parameterize the coupler's response to the two-photon state in terms of generic dimensionless variables that can be mapped to any combination of coupler and state properties. The dimensionless product MΛ, where Λ=|λ02−λ01| is the photon pair non-degeneracy and M=dκ(λ)L/dλ is the first-order coupler dispersion, gives the absolute difference in κ(λ)L between the photon central wavelengths. For discussing spectral dependencies, the product MΔλ, similarly gives the absolute difference in κ(λ)L across the FWEIM of the marginal spectra.


A convenient parameter space for navigating the coupler response can be created from η(λ00) and MΛ if the reference wavelength λ00 is taken to be the average of the photon central wavelengths λ01 and λ02. For photon pairs with a tunable non-degeneracy, such as those generated through spontaneous nonlinear interactions, λ00 can be set as the photon pair degeneracy wavelength, since λ01 and λ02 tend to remain approximately equidistant from the degeneracy point for Λ of up to hundreds of nanometers. All possible coupler responses to the quantum state then occur within the bounds η(λ00)∈[0, 1] and MΛ∈[0, π]. Behaviours for MΛ>π can be mapped back to the interval MΛ∈[0, π]



FIG. 1b shows how Δη varies within these bounds. There are four coordinates where the coupler responds as a 50:50 beamsplitter with η(λ01)=η(λ02)=0.5, and one central coordinate where it responds as a WD with Δη=1. These provide a reference for tracking transitions between beamsplitter and WD behaviour. The special condition η(λ01)+η(λ02)=1 occurs along the lines η(λ00)=0 and MΛ=π/2, where the splitting ratios at λ01 and λ02 are anti-symmetric about η=0.5. We note that if the assumptions of linear κ(λ) or λ00=|λ0102|/2 break down, the parameter space shown in FIG. 1b becomes skewed with respect to the horizontal axis.


Dispersive Coupler Example

An example of an embodiment of a photonic directional coupler is now presented. The example is intentionally simplistic to show an embodiment with limited deviation from conventional coupler designs. More optimal approaches will also be discussed.


We consider the manipulation of photon pairs degenerate at 1550 nm in the telecom band, having a maximum tunable nondegeneracy of at least Λ=50 nm. Such states can be generated through waveguide-based SPDC. In this example we seek a coupler that can reach the operating point MΛ=π/2 within this tunable range.


The design of the coupler is shown in FIG. 2a and is straightforward to fabricate with conventional coupler manufacturing techniques. FIG. 2b shows the coupling strength for the coupler of FIG. 2a in the vicinity of 1550 nm, which is linear and described by κ(λ)=2.2055×1010λ−2.0245×104 m−1. For 50:50 splitting at the degeneracy point, the smallest suitable interaction length is L=56.3 μm. From the definition of M, this yields MΛ=0.0621≈π/50 at the maximum non-degeneracy of Λ=50 nm, which is below our target. However, since M scales linearly with L, we can multiply the dispersion by choosing a larger value of L that still gives 50:50 splitting at degeneracy. An interaction length of L=1521 μm achieves this and gives MΛ=1.07×π/2 for our design, meeting our objective


While the approach above shows that the dispersion can be made arbitrarily large by increasing the device length, this comes at the price of increasing its footprint and insertion losses. Typical losses at 1550 nm for this silicon-nitride waveguide geometry are around 3 dB/cm, and hence roughly 10% of the photons would be lost in the design we described. However, this serves merely as an illustrative example. More sophisticated coupler designs exhibiting appreciable dispersion have been studied in the past including grating assisted couplers and couplers implemented in asymmetric vertical structures, including Bragg Reflection waveguides. These can alternative structures may provide more compact and efficient ways of achieving the necessary dispersion.


Calculation of Spectral Entanglement

The spectral entanglement of a state is completely described by its biphoton amplitude (BPA), which is defined in the Appendix of the main text. For a given BPA, the Schmidt Number is calculated from SN=1/[Σnp2n], where the pn are the eigenvalues of the matrix





ρωω′=∫dω″φ(ω,ω″)φ*(φ′,φ″),   (1)


and are normalized according to Σnpn=1. To quantify the entanglement of anti-bunched states at the coupler output, we associate the labels 1 and 2 with output paths A and B respectively, and post-select for terms containing





{circumflex over (b)}A†1){circumflex over (b)}B†2)|{circumflex over (v)}accustom-character.


The associated BPA is proportional to





ΞAB12)=ΦA→AB12)+ΦB→AB12),   (2)


which replaces φ(ω1, ω2) in Equation (1). For the non-path-entangled input state |ψcustom-characterA, we set ΦB→AB1, ω2) to zero.


Bandwidth Dependence


FIGS. 3a-3c, illustrate the dependence of two-photon path correlations on coupler response. The dependence changes with bandwidth.



FIGS. 3a 3c illustrate techniques for probing matter with tunable time ordering. FIG. 3a illustrates photons leaving a coupler from different output ports having two possible pathways: |λ01custom-characterA|λ02custom-characterB or |λ01custom-characterB|λ02custom-characterA. These coincide temporally and hence are mutually coherent. The photon in waveguide A is then temporally delayed by an interval τ relative to its twin photon in waveguide B, so that one photon always arrives at the sample before the other. The wavelength of the delayed photon depends on whether the pathway was |λ01custom-characterA|λ02custom-characterB or |λ01custom-characterB|λ02custom-characterA. FIG. 3b illustrates the case for μ=0, where only the |λ01custom-characterA|λ02custom-characterB pathway is allowed, such that the photon of wavelength λ02 is always absorbed first. FIG. 3c illustrates the case for μ=1, where the superposition permits two absorption pathways: λ02 followed by λ01, and λ01 followed by λ02. In certain systems where it is not possible to distinguish which of these pathways led to the final state of the sample, the pathways destructively interfere to suppress the two-photon absorption probability. Note that at μ=1 the pathways |λ01custom-characterB|λ02custom-characterA and |λ01custom-characterA|λ02custom-characterB are also present due to nondeterministic separation (the coupler behaves as a beam splitter rather than a WD), yielding photons with no relative delay. These are not time ordered but do support both absorption pathways and therefore compliment the path-interference effects.



FIGS. 4a-4f illustrate the dependence of two-photon path correlations on coupler response. Calculations depict FIG. 1a the “classical” separation probability, FIGS. 1b-c the contribution of quantum interference, FIG. 1d the resultant interference visibility, and FIGS. 1e-f total separation probability. Toggling the phase shift from θ=0 to θ=π leads to a sign change for PIS but leaves its magnitude |PIS| unaltered. This sign change, in turn, toggles the line of maximal PS between η(λ00)=0.5 and MΛ=π/2, respectively.


When the product MΔλ increases but spectral entanglement remains low (i.e. SN≈1), FIGS. 4a-f all begin to flatten. In FIG. 4a, the classical contribution PCS at all coordinates approaches a value of 0.5; the interference contributions PIS and visibility VS in FIGS. 4(b)-(d) all approach zero; correspondingly the total separation probability PS approaches 0.5 in FIGS. 4(e)-(f). In comparison, when the photons are highly frequency-entangled, increases to MΔλ do not flatten the surfaces uniformly in this way. Instead, for FIGS. 4(a)-(d) it causes the surfaces to ‘smear’ along the MΛ axis, with the effect of averaging the values along this axis. FIGS. 4(e)-(f) are exceptions: for large values of SN, PS at θ=0 remains relatively unchanged from its values at small bandwidths; however, PS at θ=π instead flattens to approach values of 0.5. These differ because the smearing of FIGS. 4(a)-(d) along the MΛ axis alters the symmetry in how the PCS and PIS contributions sum between the two cases.


Bunched Probabilities and Visibility

In addition to the separated (anti-bunched) probabilities PCS,PIS, and PS, there is naturally a complementary set of bunched probabilities PCB, PIB, and PB, corresponding to outcomes where the photons exit together from the same output port. For anticoalescence, these are related as follows: PS+PB=1; PCS+PCB=1; and |PIS|=|PIB|. It is likewise possible to define a bunched outcome interference visibility VB=|PIB|/PCB, which behaves differently from VS. The behaviour of these visibilities also depends on whether we are implementing coalescence (i.e. with photons beginning in different waveguides) or anti-coalescence (i.e. with photons beginning in the same waveguide). For simplicity, consider the familiar case where the coupler is non-dispersive and hence η is a fixed value. For coalescence such as in the HOM effect, VS=2η(1−η)/[η2+(1−η)2], while VB=1 and is independent of η because the classical and non-classical contributions to PB scale identically. These behaviours are reversed for anti-coalescence. We also note that without dispersion, both visibilities must be equal to unity for perfect coalescence or anti-coalescence to occur. However, with dispersion, this requirement is lifted.


Dispersion-Enabled Capabilities
Tunable Spectral Entanglement

Suppose two nondegenerate photons enter a directional coupler from a single input port so that the input state takes the form |ψcustom-characterin=|λ01custom-characterj02custom-characterj, where j∈{A,B}. The two-photon state at the output of the coupler is then post-selected for outcomes where the photons exit from different waveguides (i.e., separated). Depending on the coupler response, the output waveguide taken by a given photon can reveal information about that photon's spectral properties, which in turn alters the spectral entanglement of the post-selected output state. A WD-like response with Δη=1 predetermines which photon emerges from each output port. This leads to an output state of the form |ψcustom-characterout=|λ01custom-characterA|λ02custom-characterB (or |ψcustom-characterout=|λ01custom-characterB02custom-characterA, depending on the input port), where entanglement of the central wavelengths is lost. On the other hand, a beam splitter-like response with Δη=0 leads to the superposition |ψcustom-characterout=[|λ01custom-characterA02custom-characterB+|λ01custom-characterB02custom-characterA] √2, where the full spectral entanglement of the input state is retained. By controlling Δη through the selection of M or η(λ00) (and thus controlling, effectively, the amount of spectral information known about the output state), a directional coupler can select any level of entanglement between these extremes.



FIGS. 5a-b illustrate the tunability of output state entanglement. FIG. 5a illustrates the dependence of Schmidt number on the coupler response for post-selected outcomes where the photons are found in different waveguides. The maximum value of SN=2.31 corresponds to the input state entanglement. FIG. 5b is a plot of SN vs. κ(λ00)L taken for a slice along MΛ=π/2, plotted in terms of η(λ00). FIGS. 5a and 5b show how the choice of M and η(λ00) can tailor spectral entanglement in the post-selected output state. Spectral entanglement has been quantified using the Schmidt number (SN), which has a minimum value of unity in the absence of entanglement, and increases with greater entanglement. The input state used in this example has Λ=10 nm, SN=2.31, and equal FWHM intensity bandwidths of Δλ=1 nm for the photon marginal spectra. It is modeled after a Type-I spontaneous parametric downconversion (SPDC) process with a degeneracy wavelength of λ00=1550 nm and pump bandwidth of ΔλP=0.25 nm. As the coupler response moves away from the beam splitter-like coordinates and toward the WD-like coordinate at (MΛ=π/2, η(λ00)=0), the Schmidt number of the output state smoothly transitions from its input value of SN=2.31 down to a value of SN≃1.15. Note that some spectral entanglement remains at the WD-like coordinate even though the output paths reveal the central wavelengths. This is because the photon spectra are still inherently anticorrelated about their central wavelengths, due to energy and momentum conservation in the pair generation process. Such residual entanglement vanishes as Δλ→0.


In situ tuning of the Schmidt number becomes possible through active control of η(λ00). Effectively, this prepares states of the form |ψcustom-characterout=[|λ01custom-characterλ02custom-characterB+μ|λ01custom-characterB02custom-characterA]/√(1+μ2) with a tunable value of μ. We emphasize that this tuning occurs post-generation, without requiring changes to pump bandwidth, nonlinear interaction length, or any other parameters affecting the photon pair generation process. This makes it particularly well suited for tailoring spectral entanglement in a monolithically integrated setting, in applications where the photons remain path-distinguishable. In some embodiments control of η(λ00), and thereby the Schmidt number, can be achieved electro-optically or thermally, for instance, by modifying the waveguide core-cladding index contrast to systematically shift κ(λ00). In some embodiments, control of η(λ00) (i.e. “tuning”) may be accomplished by the quantum-confined Stark effect and, for certain fiber-based coupler assemblies, a micrometer-controlled waveguide separation. Operation along the line MΛ=π/2 offers the most precise control over entanglement at any nonzero Λ. The value of M is fixed but can be tailored through a judicious design of the coupler dimensions and material system. Note that, since M scales with L, dispersion can be enhanced by increasing the 50:50 coupling length beyond its minimum necessary value of L=π(4κ(λ00)).


This tuning approach also provides control over polarization entanglement, since correlations in the spectral and polarization degrees of freedom are coupled, except in the special case of maximal polarization entanglement. A state's polarization entanglement can be quantified using its concurrence C, with C=0 and C=1 indicating minimal and maximal entanglement, respectively. As the state Schmidt number increases, polarization entanglement tends to decrease, and vice-versa. This inverse relationship between SN and C allows for the on-chip preparation of non-maximally entangled states |ψcustom-character=(|H,Vcustom-character+r exp iφ|V,Hcustom-character)/√(1+r2) with a tunable value of r<1, with r related to the concurrence by C=2r/(1+r2). Such states offer significant advantages over maximally entangled states in certain applications such as closing the detection loophole in quantum nonlocality tests.


The tunable spectral entanglement we present may also have useful capabilities for two-photon spectroscopy and light-induced matter correlations. In these applications, the time ordering of when each photon reaches the sample can affect the two-photon absorption probability. This is because a particular two-photon transition can have pairings of absorption pathways corresponding to whether λ01 or λ02 is absorbed first. For some systems, when both time orderings are permitted by the incident light, these pathways destructively interfere to suppress the two-photon absorption probability, as is the case for two uncoupled two-level atoms. Such transitions can thus be selectively controlled by changing which time orderings (and hence absorption pathways) are allowed.


As illustrated in FIGS. 3a-c, control over the allowed time orderings may be achieved by placing a time delay in one path (e.g., path A) and tuning μ by tuning the coupler parameter η(λ00). Suppose μ=0 [FIG. 3b] so that the post-selected state at the coupler output is |ψcustom-characterout=|λ01custom-characterA02. In this case, λ01 is always delayed relative to λ02. Hence, only one set of time-ordered pathways is allowed. On the other hand, when μ=1 [FIG. 3c] so that |ψcustom-characterout=[|λ01custom-characterB02custom-characterA] √2, the delay is applied in superposition to either λ01 or λ02, and hence both sets of time-ordered pathways are allowed.


Such control over the time ordering adds to the versatility of a single on-chip light source for manipulating and probing two-photon processes, such as controlling the degree to which bi-exciton transitions may be blocked. Note that the ability to selectively excite a single absorption path (e.g., using μ=0) is only possible with quantum light sources. Classical sources have no intrinsic time ordering and hence will excite both paths equally (as with μ=1). Accordingly in some embodiments a tunable dispersive coupler may be provided that allows a sample's behavior for both the classical and nonclassical conditions to be directly compared, without the need to change the light source and with virtually no disruption to the experimental setup.


Perfect Anti-Coalescence with Tunable Visibility


Control over two-photon path correlations is another important ability for quantum photonics. In this section, we start by exploring how such correlations can be impacted by dispersion. We then describe how this enables conditions with no bulk optics equivalent; namely, perfect photon anti-coalescence that remains independent of the visibility of interference effects, even as this visibility is tuned via η(λ00) or MΛ.


Path correlations are commonly engineered using quantum interference. In the famous Hong-Ou-Mandel effect, two photons enter a 50:50 beam splitter from different input paths (antibunched), and coalesce to exit as a bunched state where they are most likely to be found in the same output path. Ideally the antibunched (i.e., separated) outcome probability becomes PS=0 under conditions of maximal interference, compared to the “classical” value of PCS=0.5 if interference were completely absent. The reverse process, called anti-coalescence, wherein PS→1, is useful for providing interference-facilitated pair separation (IFPS) to separate photons generated by integrated sources. Note that the subscript S is used to delineate these from probabilities corresponding to bunched (i.e., non-separated) outcomes. The two-photon interference can be quantified by the interference visibility VS=|PIS|/PCS, where PIS=PS−PCS represents the contribution of quantum interference toward the antibunched outcome probability.


We shall now look specifically at anti-coalescence. While perfect coalescence requires VS to be unity, coupler dispersion can lift this restriction for anti-coalescence. As we shall see, for the first time VS can be made to have any arbitrary value between 0 and 1 while the separation probability is kept constant at PS=1. Anti-coalescence requires a path-entangled input state of the form





custom-character=[|ψcustom-characterA|0custom-characterB+e−10|0custom-characterAcustom-characterB]/√2, |  custom-character(1)


where |0custom-character refers to vacuum, |ψcustom-characterj represents a photon pair in path j, and θ is a relative phase shift. Such states can be generated by coherently pumping two sources of photon pairs. This places no restrictions on the tunability of the photon pair sources. The spectral properties of |ψcustom-characterj are described by the biphoton amplitude (BPA) φ(ω1, ω2). We will assume perfect path indistinguishability such that φ(ω12)=φB12)≡φ(ω12).



FIGS. 4a-f illustrate how PS, PCS, PIS, and VS change as a function of the coupler parameters when the relative phase shift is either θ=0 or θ=π. These plots have been generated for a co-polarized input state from Type-I SPDC having Δλ=0.25 nm, ΔλP=0.1 nm and a degeneracy wavelength of λ00=780 nm. The value of |PIS| is maximal at coordinates where the coupler responds as a 50:50 beam splitter, and minimal when it responds as a WD. The “classical” probability PCS follows roughly the opposite trend, obtaining its maximal value of PCS=1 for a WD-like response, and decreasing to PCS=0.5 for beam splitter-like responses. Curiously, along the lines η(λ00)=0.5 and MΛ=π/2, changes to PCS and |PIS| are in perfect balance such that their sum always equals unity. This balancing is associated with the condition η(λ01)+η(λ02)=1, which leads to PS=1 and hence perfect anti-coalescence (i.e., deterministic separation) along either η(λ00)=0.5 or MΛ=π/2, selected through the choice of θ. Along these two lines, the interference visibility VS varies smoothly between 0 and 1. By operating at MΛ=π/2 with θ=π, and actively controlling η(λ00) through thermal or electro-optic tuners, any value of VS can be selected while maintaining a perfect separation fidelity. Note that, unlike before, this does not alter the spectral entanglement of post-selected output states, due to the presence of path entanglement at the input.


Opportunities for State Characterization
Entanglement-Sensitive Coincidence Detection

For most permutations of coupler and state attributes, PS is accurately described by the behavior in FIGS. 4a-f. However, deviations from the values of PS shown can occur when the dimensionless product MΔλ, involving coupler dispersion and photon bandwidth, becomes large. FIG. 6 indicates that the extent of these deviations depends not only on MΔλ, but also on the spectral entanglement of the input state. This opens up the possibility of discerning the Schmidt number of the input state from the anti-bunched coincidence count rate at the coupler output, which is proportional to PS.


The results in FIG. 6 were calculated for degenerate input states having Λ=0 nm, Λλ=10 nm, λ00=780 nm, and θ=0. The product MΔλwas swept by varying M, with η(λ00)=0.5 kept constant. Input state entanglement was controlled through the Type-I SPDC pump bandwidth ΔλP. In the limit of MΔλ→0, the above calculation parameters give PS=1, in agreement with FIG. 4e. Larger values of MΔλ lead to decreases in PS. However, increasing the SN of the input state has the effect of asymptotically restoring PS to unity.


State Representation: A co-polarized pair with both photons beginning in waveguide j can be represented by the pure state





|ψ>1=∫dω12φj12j†1j†2)|vaccustom-character,|  custom-character(A.1)


where âj†(ω) is the canonical mode creation operator for waveguide j. The BPA is normalized according to ∫dω12|φ(ω12)|=1. Rather than generating the BPA from device-specific mode dispersion parameters, it is more convenient to define the BPA directly in terms of the photon bandwidths and central wavelengths of interest. A BPA that mimics the output of a Type I SPDC process can be constructed from





φ(ω12)=φP12)[φ1122)+φ2112)],|   custom-character(A.2)


where φn(ω) are the marginal photon spectra and φP12) is the pump spectrum. This construction satisfies the necessary exchange symmetry and has all of the key qualitative features of a typical Type I BPA computed from SPDC theory. The marginal spectra were Gaussian and defined in terms of wavelength as φ(λ)=exp(−2 ln 2[λ−λ0n]2/Δλ2), with equal FWHM intensity bandwidths of Δλ. The pump spectrum was also Gaussian with a FWHM intensity bandwidth of ΔλP. Narrowing ΔλP below Δλ has the effect of increasing the spectral correlations, and hence Schmidt Number, of the two-photon state.


Evolution through a directional coupler: Consider the evolution of the pure state |Ψcustom-character of Eq. (A.1) through a directional coupler of length L and coupling strength (ω). It is assumed that the output remains in a pure state. Let {circumflex over (b)} if (ω) represent the mode operators at the coupler output. These are related to the input mode operators by










[






b
^


A
+




(
ω
)









b
^


B
+




(
ω
)





]

=




[




cos


(


κ


(
ω
)



L

)





i






sin


(


κ


(
ω
)



L

)








i






sin


(


κ


(
ω
)



L

)






cos


(


κ


(
ω
)



L

)





]



[






a
^


A
+




(
ω
)









a
^


B
+




(
ω
)





]


.

|







(

A

.3

)








Note that the magnitude of the matrix elements in Eq. (A.3) are related to the power-splitting ratio by |cos(κ(ω)L)|=[η(ω)]1/2 and |sin(κ(ω)L)|=[1−η(ω)]−1/2. Using this transformation, the state BPAs at the output of the coupler can be written as












Φ

j

pq




(


ω
1

,

ω
2


)


=



φ
j



(


ω
1

,

ω
2


)





G

j

p




(

ω
1

)





G

j

q




(

ω
2

)




,

|




Where







(

A

.4

)









G

j

q




(
ω
)


=

{






cos


(


κ


(
ω
)



L

)


,






if





j

=
q

,







sin


(


κ


(
ω
)



L

)


,





if





j



q
.





|








(

A

.5

)








In terms of our notation, Φj→p12) is the amplitude associated with photons 1 and 2 being coupled from input path j to output paths p and q, respectively. While the form of Eq. (A.4) is general, the Gj→q(ω) will change if a different coupler architecture is used (such as an asymmetric coupler).


Two-photon outcome probabilities: The probability of finding photons 1 and 2 in output paths p and q, respectively, is calculated from Ppq=custom-characterΨ|{circumflex over (b)}p†{circumflex over (b)}q†{circumflex over (b)}q{circumflex over (b)}pcustom-character and found to be






P
pq
=R
pq
C+cos(πδpq)RpqI(θ),|  custom-character(A.6)


where δpq is the Kronecker delta,






R
pq
C
=∫dω
1

2(|ΦA→pq12)|2+|ΦB→pq12)|2),|  custom-character(A.7)


is the “classical” probability contributed by sources A and B in the absence of interference, and






R
pq
I(θ)=∫122Re{e−ωΦB→pq1,1072)Φ′A→pq12)},|  custom-character(A.8)


is a nonclassical modifier accounting for the effects of path interference. These expressions are given in their most general form so that they can be readily applied to any arbitrary set of coupler and two-photon state attributes. Note that ΣpqPpq=1. The probability PS of obtaining an antibunched (separated) outcome is then






P
S
=P
AB
+P
BA
=P
S
C
+P
S
J,|  custom-character(A.9)


with “classical” and “interference” components given by PCS=RCAB+RCBA and PIS=RIAB+RIBA.


Obtaining VS for on-chip measurement of Λ: We refer to the configuration shown in FIG. 7. Let P(Λ,τ) represent the total antibunched outcome probability at nondegeneracy Λ and relative time delay τ. Assuming (λ00)=0.5 and θ=0, PS(Λ,0)=1 at all values of Λ. The coincidence detection rate R0 at zero delay therefore corresponds to maximum separation fidelity; thus the probability of pair separation at nonzero delay τ can be obtained from P(Λ,τ)=Rτ/R0. Provided τ is large enough that |ψcustom-characterA and |ψcustom-characterB (the possible photon-pair histories) are no longer coherent, quantum interference will not occur at that delay time; thus, PI(Λ,τ)=0 and PS(Λ,τ)=PCS(Λ,τ). It then follows from the definition of VS that










V
S

=







P
S



(

Λ
,
0

)


-


P
S



(

Λ
,
τ

)







P
S



(

Λ
,
τ

)



=





1
-


R
τ

/

R
0







R
τ

/

R
0



=







R
0


R
τ


-
1



.

|









(

A

.10

)








For η(λ00)=0.5, the visibility VS maps to a unique value of MΛ provided MΛ≦π/2[due to periodicity of VS; see FIG. 4d].


This behavior can be understood by examining Eqs. (A.3)-(A.9). The probability PS is determined from a sum over all possible combinations of frequencies ω1=2πc/λ1 and ω2=2πc/λ2 weighted by the BPA. When the state is spectrally uncorrelated (i.e., SN=1), the combinations of η(λ1) and η(λ2) contributing to this sum are not necessarily equidistant from η(λ00)=0.5 and hence can deviate from the η(λ1)+η(λ2)=1 condition required for perfect anti-coalescence. However, when the photons are spectrally anticorrelated due to entanglement, the BPA restricts all contributing λ1, λ2 combinations to be approximately equidistant from λ00, which acts to restore the splitting ratio antisymmetry. Larger products of MΔλ allow PS to be more severely degraded because a greater proportion of the nonvanishing λ1, λ2 combinations are able to violate the antisymmetry. Only in the limit of Δλ→0, where the state is entirely described by the central wavelengths λ01 and λ02, is the splitting ratio antisymmetry condition strictly enforced.


The bandwidth and entanglement sensitivity of PS grants dispersive couplers additional capabilities for state characterization. For example, dispersive couplers could empower a simple, fast, all-integrated technique for measuring the Schmidt number of an ensemble of states without needing to perform full state tomography to reconstruct the BPA. FIG. 7 illustrates an embodiment of a device for measuring the Schmidt number of an ensemble of states. In this case, we show the photons being characterized immediately after leaving the source, in the context of source calibration. However, they could also be measured after interacting with a bath or system. This could be helpful, for example, in metrological applications where the Schmidt number is monitored as an indication of state purity and hence the interaction under investigation. To obtain SN, first the marginal photon spectra are measured with a waveguide-assisted spectrograph method that uses chromatic group velocity dispersion (GVD) to map spectral components to time-of-arrival at a single-photon detector. Next, provided M is known, the values of Δλ and Λ measured in the first step are used to discern SN from standard two-photon coincidence measurements at the coupler outputs. The sensitivity of the technique diminishes as the photons are made narrowband or increasingly entangled, but can be enhanced by designing the coupler to have M as large as possible.


Obtaining SN by previous methods would require a measurement of the full BPA, which hinges on the spectral resolution of the measurement system. Measuring the BPA entirely on chip is possible using spectrographs, but its resolution can be severely limited by detector timing jitter. In comparison, precise values of Δλ and Λ for the coupler-assisted technique are more easily obtained, in part due to the straightforward use of interpolation to increase confidence in these values, but also because uncertainties from the limited spectral resolution enter only in one axis, as opposed to two. Hence, the trade-offs between the number of measurements, the total measurement time, and precision in SN scale more favorably for the coupler-assisted technique. A direct, rapid, and precise measurement of SN would be particularly useful for the real-time monitoring of sources where SN is tunable and is being used as a control parameter. Additionally, it would be advantageous for monitoring a stream of states whose properties reveal real-time information about a dynamic system or environment.


The converse functionality—estimating the photon bandwidth for a known Schmidt number—could also be useful, in the context of indistinguishable pure photons having tunable attributes. As long as the SN remains reasonably close to unity, Δλ could be measured entirely on chip using only the coupler and coincidence detectors, without the need for tunable bandpass filters, GVD fibers, or spectrometer capabilities. Presently, highly bandwidth-tunable pure photons can be generated in a free-space setup, but recent trends toward integration suggest that this capability may eventually be available in a monolithic platform, where on-chip characterization would be helpful for source calibration and monitoring drift.


Versatility of Dispersive Couplers

Referring to FIG. 7, an embodiment of an all-integrated SN measurement is presented. To apply the technique of the embodiment of FIG. 7, the photon pairs must be in the generic path-entangled state |ψcustom-character of Eq. (1). The relative phase is ideally θ=0; for other values of θ, PS is less sensitive to SN. To measure SN, the state is sampled at three locations (shown as Y-junctions for simplicity). Detectors A and B sample the two-photon statistics at the coupler output to obtain PS. Detector C obtains spectrographs, and hence Λ and Δλ, by sampling |ψcustom-character via a high-dispersion element such as a fiber or a waveguide grating operated near its band edge. It is sufficient to measure these spectrographs from only one of the source output paths, since the photon pair properties are assumed to be path-indistinguishable (i.e., |ψcustom-characterA=|ψcustom-characterB). The data obtained for Λ and Δλ (together with the dispersive coupler attributes) can then be used to map the measured PS to a corresponding value of SN (e.g. see FIG. 6).


Referring to FIG. 8, an embodiment of a multipurpose dispersive coupler for state characterization is presented. A path superposition of the form |Ψcustom-character (Eq. A.1) is created through coherent pumping of two waveguide sources of photon pairs (e.g., generated via parametric downconversion). A tunable MachZehnder interferometer (MZI) allows the relative time delay to be set to either zero (φτ=0) or τ (φ96 =π) by selecting between paths, with the longer path introducing a time delay τ when the photon is switched onto the longer path. Pump power can be adjusted between paths via φP to compensate for asymmetric losses when the delay of τ is implemented. Unconverted pump photons are removed using ring filters. MZIs at the output can be toggled (φB=π) to sample the two-photon correlations with single-photon detectors. The rate of detection coincidences for zero time delay and a delay of τ can be used to determine VS, which in turn reveals MΛ. The dispersive directional coupler must have η(λ00)=½ for this measurement. Note that adding electro-optic or thermal tuners to the dispersive coupler can enable arbitrary control over VS by tuning η(λ00). Spectral-entanglement tuning is also possible when φP is set to deliver pump power to only one of the two photon pair sources.


Since couplers are already an essential on-chip device, the state characterization capabilities granted to them by dispersion can be exploited with minimal increase to the circuit complexity or footprint. This allows dispersive couplers to provide an extremely versatile set of functionalities in a compact form factor, which the following example highlights. Consider the reconfigurable circuit in FIG. 8. The dispersive coupler in this circuit can serve several purposes. It can provide IFPS to deterministically separate the photons at the coupler output. With the addition of electro-optic or thermal tuning, it can also be utilized for other previously described state engineering functionalities, such as tunable spectral entanglement. On top of this, the circuit could easily be modified for coupler-based SN measurements by tapping photon source B with a high-dispersion element and an additional detector, as per FIG. 7. Accomplishing all of these tasks through a single dispersive coupler may help to make most efficient use of precious on-chip real estate.


Even without adding a tap to source B for a spectrograph measurement, the circuit in FIG. 8 can already access some information about the state. The relationship between VS and Λ described above provides a route for measuring the nondegeneracy Λ of an ensemble of states entirely on chip. This requires the toggling of a time delay τ between the dispersive coupler input paths. The interference visibility is obtained from VS=|R0/Rτ−1|, where R0 is the coincidence count rate at zero time delay (as measured by on-chip single photon detectors), and R96 is the coincidence rate at a time delay τ that is much larger than the two-photon coherence time. Provided M is known, this value of VS can be mapped back to the state nondegeneracy Λ, as per FIG. 4d. This technique is best applied to narrowband photons since the sensitivity of VS to Λ decreases as MΛ becomes large.


Integrated couplers are already becoming a key building block of photonic quantum circuits. This is partly because they offer greater stability and scalability than bulk-optics beam splitters and other bench-top components. It is also because the highly precise micron-scale fabrication of such couplers helps eliminate path-length mismatches and other path asymmetries, which is critical for achieving high-fidelity quantum interference. However, in addition to these known benefits, our work has revealed an as-of-yet untapped potential for integrated couplers to be utilized in a more versatile way, far beyond their traditional role as a beam splitter substitute.


We found that harnessing the full dispersion properties of an integrated directional coupler unlocks many novel capabilities for the device. These include tunable photon entanglement and time ordering, as well as bandwidth-sensitive and entanglement-sensitive two-photon effects that can be exploited for state characterization. Some of these capabilities can be achieved in bulk optics, but not with the convenience nor stability that this integrated approach provides. Yet others have no bulk-optics counterpart, such as the ability to fully tune the two-photon interference visibility (i.e., the sensitivity to time delays at the coupler input) while maintaining a constant flux of separated (i.e., anticoalesced) photon pairs. Particularly remarkable is that all of these functionalities can be provided by a single integrated coupler, making it a versatile yet compact tool for both state engineering and on-chip state characterization. This is made possible by the capacity of dispersive couplers to smoothly transition between the extremes of beam splitter and wavelength-demultiplexer behavior, in a manner without parallel in bulk optics.



FIGS. 9 & 10 illustrate some useful equations when reviewing the above description.

Claims
  • 1. A method for selecting a level of entanglement between two nondegenerate photons comprising: receiving two non degenerate photons through a single input port of a directional photonic coupler;adjusting one of a first-order coupler dispersion M or a power splitting ratio η(λ00) of the directional optical coupler to select a Δη; and,emitting the photons from corresponding output ports of the directional optical coupler, wherein the emitted photons have a spectral entanglement corresponding to the selected Δη.
  • 2. The method of claim 1, wherein the adjusting comprises adjusting a waveguide core-cladding index of the directional optical coupler.
  • 3. The method of claim 1, wherein the adjusting comprises adjusting a separation of waveguides of the directional optical coupler.
  • 4. The method of claim 1, further comprising selecting a level of time-ordering between the two nondegenerate photons by, applying a time delay τ to one of the output ports, wherein the time delay converts the selected Δη path-entanglement characteristics into time-ordering characteristics of the emitted photons.
  • 5. The method of claim 1, further comprising selecting a level of interference visibility at a fixed photon pair anti-bunching (separation) rate for the two nondegenerate photons, wherein the receiving two non degenerate photons further comprises receiving two nondegenerate photons that are path-entangled across two input ports of the directional optical coupler with a relative phase of θ=π between paths; and, wherein the adjusting comprises: for a given M, adjusting η(λ00) of the directional optical coupler to select a desired interference visibility.
  • 6. A device for selecting a level of entanglement between two nondegenerate photons comprising: a directional optical coupler adapted to receive two non degenerate photons through a single input port; and,means for adjusting a power splitting ratio η(λ) of the directional optical coupler to a selected Δη;wherein emitted photons from the directional optical coupler have a spectral entanglement corresponding to the selected Δη.
  • 7. The device of claim 6, wherein the means for adjusting the power splitting ratio η(λ) comprises means for adjusting a waveguide core-cladding index of the directional optical coupler.
  • 8. The device of claim 6, wherein the means for adjusting the power splitting ratio η(λ) comprises means for adjusting a separation of waveguides of the directional optical coupler.
  • 9. The device of claim 6, wherein the means for adjusting the power splitting ratio η(λ) comprises a coupling length of the directional optical coupler longer than a minimum necessary value (L=π/(4κ(λ00)).
  • 10. The device of claim 6, further operative to select a level of time-ordering between the two nondegenerate photons, wherein the device further comprises an extended path in one leg of the directional optical coupler, the extended path applying a time delay τ to one of the output ports, wherein the time delay τ converts the selected Δη path-entanglement characteristics into time-ordering characteristics of the emitted photons.
  • 11. The device of claim 6, further operative to select a level of interference visibility at a fixed photon pair anti-bunching (separation) rate for the two nondegenerate photons, wherein the directional optical coupler is adapted to receive the two non degenerate photons as path-entangled photons across two input ports of the directional optical coupler with a relative phase of θ=π between paths; and, for a given M, the means for adjusting a power splitting ratio η(λ) of the directional optical coupler comprises means for adjusting η(λ00) of the directional optical coupler to select a desired interference visibility.
CROSS REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of priority to U.S. Provisional Patent Application Ser. No. 62/354,906 filed on Jun. 27, 2016 and entitled System and Method for Dispersion-Enabled Quantum State Control of Photons, the contents of which are incorporated by reference.

Provisional Applications (1)
Number Date Country
62354906 Jun 2016 US