METHOD AND APPARATUS FOR DISTRIBUTING A QUANTUM KEY

Abstract

A method for distributing a quantum key is provided, including sending a first photon to a first receiver; sending a second photon to a second receiver, the first and second photons being a pair of time-energy entangled photons; and providing a coding scheme comprising a plurality of time bins and a plurality of frequency bins, wherein a combination of a time bin and a frequency bin corresponds to a character.

Description

FIELD

The disclosed subject matter generally relates to integrated optical systems. More particularly, the disclosed subject matter relates to a new design for a photonic crystal useful for a wide variety of applications.

BACKGROUND

Quantum key distribution (QKD) enables two parties to establish a secure key at a distance, even in the presence of one or multiple eavesdroppers. The key can then be used to secretly transmit information using unconditionally secure one-time pad encryption. Unconditional security is achieved by the laws of physics rather than by assumptions about the computational abilities of the eavesdropper. Various QKD protocols have been proposed and implemented. A first protocol due employed polarization states of photons passed between the two parties. There has been growing interest in schemes employing photons in Hilbert spaces of high dimension, resulting in a potentially very large alphabet size. Different degrees of freedom have been considered, including polarization, time, and spatial modes.

What is needed is a protocol that allows two parties to generate their secure key at the maximum rate allowed using the time-energy basis. Moreover, a protocol is needed that is compatible with fiber-based dense WDM (DWDM) systems commonly used in classical fiber communications. It is also desirable to have an extremely compact, stable, and scalable platform for the protocol's implementation.

SUMMARY

A method for distributing a quantum key is provided, including sending a first photon to a first receiver; sending a second photon to a second receiver, the first and second photons being a pair of time-energy entangled photons; and providing a coding scheme comprising a plurality of time bins and a plurality of frequency bins, wherein a combination of a time bin and a frequency bin corresponds to a character.

In some embodiments, the first photon and the second photon are sent via optical fiber. In some embodiments, the first photon and the second photon are sent via a photonic integrated chip. In some embodiments, the first photon and the second photon are sent through free space.

A method for distributing a quantum key is provided including generating a plurality of pairs of time-energy entangled photons; for each pair of time-energy entangled photons, sending one photon to the first receiver and one photon to the second receiver; and providing a coding scheme comprising a plurality of time bins and a plurality of frequency bins wherein each pair of one of the plurality of time bins and one of the plurality of frequency bins corresponds to one of a plurality of characters.

A method of receiving a quantum key is provided including receiving a first photon, the first photon being one of a pair of time-energy entangled photons; detecting a time of arrival of the first photon; detecting a frequency of the first photon; and determining a character based on the detected time and frequency of the first photon.

In some embodiments, the determining a character step further comprises: assigning the detected time to a time bin; assigning the detected frequency to a frequency bin; and determining the character from a coding scheme based on the time bin and the frequency bin.

In some embodiments, the method further includes detecting whether an eavesdropper has detected the first photon by determining the frequency distribution of the first photon.

A method of receiving a quantum key is provided including receiving a plurality of photons, each of the plurality of photons being one of a pair of time-energy entangled photons; detecting the time of arrival of each of the plurality of photons; detecting the frequency of each of the plurality of photons; and determining a character based on the detected time and frequency of each of the plurality of photons.

An apparatus for receiving a quantum key is provided including a channel configured to receive photons; a plurality of multi-channel filtering elements, each corresponding to a frequency and each configured to pass photons within a range of its corresponding frequency; and a plurality of photon detectors, each configured to receive photons passed by one of the plurality of multi-channel filtering elements and to send an indication of photon detection.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1( a) illustrates a system using the time-coding scheme in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 1( b) illustrates time bins in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 1( c) illustrates a further system using the time-coding scheme in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 1( d) illustrates time bins in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 2( a) illustrates Franson interferometers in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 2( b) illustrates the coincidence counting probability P_C displayed versus δt in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 3( a) illustrates mutual information as a function of detector jitter, σ_det, normalized to σ_bin in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 3( b) illustrates mutual information as a function of frequency channels, n_f in accordance with exemplary embodiments of the disclosed subject matter.

FIG. 4 illustrates a further Franson interferometer in accordance with exemplary embodiments of the disclosed subject matter.

DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS

It is understood that the subject matter described herein is not limited to particular embodiments described, as such may, of course, vary. It is also understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting, since the scope of the present subject matter is limited only by the appended claims. Where a range of values is provided, it is understood that each intervening value between the upper and lower limit of that range, and any other stated or intervening value in that stated range, is encompassed within the disclosed subject matter.

Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosed subject matter belongs. Although any methods and materials similar or equivalent to those described herein can also be used in the practice or testing of the present disclosed subject matter, this disclosure may specifically mention certain exemplary methods and materials.

All publications mentioned in this disclosure are, unless otherwise specified, incorporated by reference herein for all purposes, including, without limitation, to disclose and describe the methods and/or materials in connection with which the publications are cited.

The publications discussed herein are provided solely for their disclosure prior to the filing date of the present application. Nothing herein is to be construed as an admission that the present disclosed subject matter is not entitled to antedate such publication by virtue of prior invention. Further, the dates of publication provided may be different from the actual publication dates, which may need to be independently confirmed.

As used herein and in the appended claims, the singular forms “a,” “an” and “the” include plural referents unless the context clearly dictates otherwise.

Nothing contained in the Abstract or the Summary should be understood as limiting the scope of the disclosure. The Abstract and the Summary are provided for bibliographic and convenience purposes and due to their formats and purposes should not be considered comprehensive.

As will be apparent to those of skill in the art upon reading this disclosure, each of the individual embodiments described and illustrated herein has discrete components and features which may be readily separated from or combined with the features of any of the other several embodiments without departing from the scope or spirit of the present disclosed subject matter. Any recited method can be carried out in the order of events recited, or in any other order that is logically possible.

Reference to a singular item includes the possibility that there are plural of the same item present. When two or more items (for example, elements or processes) are referenced by an alternative “or,” this indicates that either could be present separately or any combination of them could be present together, except where the presence of one necessarily excludes the other or others.

As summarized above and as described in further detail below, in accordance with the various embodiments of the present invention, a large-alphabet protocol is described herein referred to as a ‘Wavelength Division Multiplexed Quantum Key Distribution’ (WDM-QKD) that employs time-energy entangled photon pairs. Given a certain photon flux n and channel bandwidth ΔΩ, WDM-QKD allows two parties, e.g., “Alice and Bob,” to generate their secure key at the maximum rate allowed using the time-energy basis. WDM-QKD forms the quantum-cryptography analog of present-day WDM systems and is compatible with fiber-based dense WDM (DWDM) systems commonly used in classical fiber communications. An implementation of WDM-QKD is described herein as both a modern fiber optic network and in a photonic integrated chip (PIC). The latter provides an extremely compact, stable, and scalable platform for the protocol's implementation.

WDM-QKD enables Alice and Bob to generate their shared key at a maximum rate of n log₂(ΔΩ/n) at log₂(ΔΩ/n) bits per photon (bpp), using present-day single photon counters, WDM equipment, and simple components such as beam splitters or directional couplers. The information per photon per bandwidth is log₂ (ΔΩ/n)/ΔΩ.

QKD protocols derive security from the fact that a measurement changes an unknown state if measured in a conjugate basis. If Alice and Bob employ two conjugate bases, they can detect a measurement by a third party eavesdropper, “Eve,” including quantum nondemolition measurements. In the protocol described herein, the conjugate bases are time and energy. Alice and Bob make time and energy measurements using an extended Franson interferometer. Assuming 400 wavelength channels with detectors having 40 ps jitter, it is expected that for a bandwidth ΔΩ=10 nm around 1550 nm, Alice and Bob can generate a key at a maximum of 10 Tera bits per second (Tbps) at 1 bit per photon (bpp); or at 10 bpp and 100 Gbps; or at 200 Mbps at 20 bpp.

Time-energy entangled photon pairs with correlation time σ_cor generated by spontaneous parametric down conversion (SPDC) were considered herein, assuming a pump field at frequency ω_p with coherence time σ_coh. The photon pair wave function can be written as

|Ψ=∫_−∞^∞∫_−∞^∞ψ(t_A,t_B)|t_A,t_B,o_A,o_Bdt_Adt_B

|t_A,t_B,o_A,o_B=â_oA^t(t_A)â_oN^t(t_N)|0

ψ(t_A,t_B)∝e^{−(tA−tB)2/4σcor2}e^{−tA+tB)2/16σco2}e^{iωp/2(tA+tB)}  (1)

The creation operator â_oi^t(t_j) denotes creation in spatial mode o_i at time t_j. If the correlation time of the photons is sufficiently smaller than the time bin duration, Alice and Bob can establish a secret key using the correlated photon arrival time. A schematic of this protocol is shown in FIG. 1. The continuous biphoton wave function can be discretized as a sum over time bins |σ_binⁱ of duration σ_bin as

$\begin{array}{cc}\stackrel{\_}{\Psi }〉=\sum _{i=-\infty }^{\infty }\sum _{j=-\infty }^{\infty }G^{i,j}{\sigma }_{\mathrm{bin}}^i,{\sigma }_{\mathrm{bin}}^j〉G^{i,j,}={\int }_{i{\sigma }_{\mathrm{bin}}}^{\left(i+1\right){\sigma }_{\mathrm{bin}}}{\int }_{j{\sigma }_{\mathrm{bin}}}^{\left(j+1\right){\sigma }_{\mathrm{bin}}}\psi \left(t_A,t_B\right)t_At_B.& \left(2\right)\end{array}$

The probability that Alice and Bob project into time bins σ_binⁱ and σ_bin^j is therefore p^i,j=σ_binⁱ,σ_bin^j| ψ|²=|G^i,j|². Detector jitter of magnitude σ_det influences the fidelity of this projection, so that roughly σ_bin>σ_det for reliable communication. Jitter therefore reduces the maximum number of characters to σ_coh/σ_det from the maximum allowed, given by the Schmidt number, K≈σ_coh/σ_cor. The fastest single photon detectors provide σ_det≈40 ps, whereas σ_cor can be on the order of 10 fs-1 ps.

FIG. 1( a) illustrates a system 100 using the time-coding scheme. A strong laser (not shown) pumps a nonlinear crystal. Photons pairs are generated by SPDC 102 and sent across channels of equal length to detectors 106 of Alice's computer 108 and Bob's computer 110 who measure their arrival times. FIG. 1(b) illustrates time bins agreed upon by Alice and Bob over a public channel. If a photon pair is detected in a given time bin, then that character is shared between Alice and Bob. FIG. 1(c) illustrates a system 200 in which photons pairs are also generated by SPDC 202. Alice's computer 208 and Bob's computer 210 include a grating, dispersive element or multi-channel filter 212 before their detectors 206 to obtain frequency information. FIG. 1(d) illustrates time bins.

For technological reasons, detector jitter is unlikely to approach the sub-ps regime in the near future. To overcome the technological mismatch between photon correlation time and detector jitter, the protocol described herein utilizes the circumstance that photons are not only entangled in time, but also in frequency. The biphoton wave function in the frequency domain, |Ψ_F=FT₂|Ψ, where FT₂ denotes the two-dimensional Fourier transform. Therefore, |Ψ_F=∫∫(ω_A,ω_B)|ω_A,ω_B,o_A,o_Bdω_Adω_B, where (ω_A,ω_B)α exp[−σ_cor²/4(ω_A−ω_B)²]exp[−σ_coh²/4(ω_A+ω_B−ω_p)²].

States of this form show non-local frequency correlations. Thus, if Alice measures a frequency ω_A on one photon, then Bob must measure a frequency ω_B≈ω_p−ω_A on the other, where ω_B=ω_p−ω_A for σ_coh→∞. If both Alice and Bob place multi-channel filtering elements before their detectors (as shown in FIG. 1(c)), then they can collect frequency information in addition to timing information. This filtering can be modeled as Gaussian projections onto discrete output spatial modes |ζⁱ with frequency bandwidth δv and center frequency v_i. Such filtering increases the correlation time of the photon pair, but this can be kept lower than the detector jitter to minimize the loss of timing information. Alternatively, one can increase the time bin duration σ_bin to allow for smaller δv. Optimization of these parameters given limitations of the source and transmission channel are described below.

The QKD protocol is secure as discussed herein. Measurements of the frequency and creation time of the biphoton packet disturb its wave function in these bases; temporal positive operator valued measurements (POVM) reduce the coherence time while frequency POVMs increase the correlation time as discussed in greater detail below.

Alice and Bob can check the security of their channel by testing for changes in the correlation and coherence times. They use an extended Franson interferometer (eFI) 300/301, depicted in FIG. 2(a), in which photons pairs are also generated by SPDC 302. eFI 300/301 include, e.g., beam splitters BS 326, partial beam splitters PBS 328, and detectors 306. Measurement in Alice's and Bob's arm of the eFI are described by the annihilation operators for the long and short paths, with respective times t_Lⁱ and t_Sⁱ, and iεA, B denoting Alice and Bob: âA(tA)=1/√{square root over (2)}[âA(t_L^A)+âA(t_S^A)] and âB(tB)=1/√{square root over (2)} [âB(t_L^B)+âA(t_S^B)].

FIG. 2( a) illustrates for Bob's eFI 301, the short arm is switched between position A and B to achieve lengths δt₁ and δt₂. This allows determination of both and σ_cor so that weak frequency and time measurements on the photon pair can be detected. FIG. 2(b) illustrates the coincidence counting probability P_C displayed versus δt. For example, at position A, δt=0 and P_c=1.0.

These times are redefined using standard notation, where Δt=(t_L^A−t_S^A) and δt=(t_L^A−t_S^A)−(t_L^B−t_S^B). A selection is made to use Δt, Δt−δt>>σ_cor in order to avoid single photon interference between long and short paths of a single arm of the eFI. Selection of two different values, δt_1,2, is allowed by placing a switch in the long path of one arm of the eFI. The parameter δt_1,2 is then scanned to provide phase-dependent coincidence counting measurements. Alice and Bob measure the coincidence counting probability P_C, which evaluates to P_c∝½+½ cos [ω(2Δt−δt)]e^{−δt2/8σcoh2} e^{−Δt2/8σcoh2} with visibility V=e^{−δt2/8σcor2}e^{−Δt2/8σcoh2}. The correlation time and coherence time can be deduced from two visibility measurements V₁ and V₂ using δt₁ and δt₂, respectively. These extrapolated values σ_coh^E′ and σ_cor^E′, are given by

$\begin{array}{cc}{\left({\sigma }_{\mathrm{cor}}^{E^{\prime }}\right)}^2=\frac{1}{8}\frac{\delta t_1^2-\delta t_2^2}{\ln V_1-\ln V_2}& \left(3\right)\\ {\left(\text{?}\right)}^2=\frac{1}{8}\frac{\Delta t^2\left(\delta t_1^2+\delta t_2^2\right)}{\delta t_1^2\ln V_2-\delta t_2^2\ln V_1}\text{?}\text{indicates text missing or illegible when filed}& \left(4\right)\end{array}$

The parameters δt₁ and δt₂ are set so that δt₁−δt₂≈σ_cor and |δt₁≠|δt₂| to sample the Franson curve at different points.

FIGS. 2( a) and 2(b) illustrate the choice of δt₁=0 and δt₂≈σ_cor. Using (σ_coh^E)²=1/[σ_coh^E′)⁻²−σ_coh⁻²] and (σ_cor^E)²=1/[(σ_cor^E′)⁻²−σ_cor⁻²] derived from this measurement, the bound on Eve's information per photon is I_E≦log₂(σ_coh/σ_coh^E)+log₂(σ_cor^E/σ_cor), which is the sum of her information obtained from temporal and frequency measurements. While all change in the wave function can be due to Eve's action, the analysis does not assume that Alice and Bob have perfect detectors or receive unaltered photon pairs. A more pessimistic calculation is required for the mutual information between Alice and Bob.

The probability of Alice and Bob measuring in frequency channels i and j and time bins k and l, is p^ijkl=ζ_Aⁱ,ζ_B^j,σ_bin,A^k,σ_bin,B¹|Ψ|²=|G^ijkl|². These coefficients form a joint probability density function which can be used to compute the Shannon entropy S=Σ_ijklp^ijkl log p^ijkl and mutual information I (A, B)=S(A)+S(B)−S(A, B).

$\begin{array}{cc}I\left(A,B\right)=-2\sum _i^{}{{\Gamma }^{\mathrm{ik}}}^2\log {{\Gamma }^{\mathrm{ik}}}^2+\sum _{i,j}^{}{\text{?}}^2\log {\text{?}}^2{\Gamma }_{\mathrm{ik}}={\int }_{k{\sigma }_{\mathrm{bin}}}^{\left(k+1\right){\sigma }_{\mathrm{bin}}}{\int }_{-\infty }^{\infty }{\mathrm{FT}}_2[\text{?}{\int }_{-\infty }^{\infty }\psi \left({\omega }_A,{\omega }_B\right){\omega }_A{\omega }_B]t_At_B\text{?}\text{indicates text missing or illegible when filed}& \left(5\right)\end{array}$

Detector jitter is included by operating a Gaussian spreading function on the biphoton state as {circumflex over (σ)}_det∫e^{−tx2/2σdet2}|tt+t_x|dt_x. FIG. 3 shows the mutual information as a function of the detector jitter normalized to the time bin duration, as discussed hereinbelow. The joint probability density function P of the new state, |Ψ={circumflex over (σ)}_det,A{circumflex over (σ)}_det,B{Ê}Ŝ_BŜ_A|Ψ_F, incorporating the set of Eve's POVMs, {Ê}, has elements p^ijkl=ζ_A^kζ_B^l,σ_bin,Aⁱ,σ_bin,B^j|Ψ′².

FIG. 3( a) illustrates mutual information as a function of detector jitter, σ_det, normalized to σ_bin. FIG. 3(b) illustrates mutual information as a function of frequency channels, n_f. The individual frequency channels have bandwidth δv=Δω/n_f and spacing 4 δv. Increasing the number of frequency channels increases the bits per photon, until the filter photon temporal bandwidth approaches σ_bin.

The number of time- and frequency-encoded bits given a photon budget and channel bandwidth can be optimized. For example, the photon budget is specific to the photon pair source and is given by the maximum emission rate, R_v. The allotted channel frequency bandwidth is Δω. The frequency bandwidth of an individual photon pair source is Δω∝1/σ_cor, so the number of sources to fill the channel bandwidth is N=ΔΩ/Δω. The communication rate, R_C, in bits per second is therefore

$\begin{array}{cc}{R}_C={\mathrm{NR}}_v{\log }_2\left(\frac{\mathrm{\Delta \omega }/\delta v}{R_v{\sigma }_{\mathrm{bin}}}\right)& \left(6\right)\end{array}$

using σ_bin and δv defined above.

The often limited photon budget for quantum communication makes high-dimensional encoding desirable. However, achieving the limit on this dimensionality in the time domain using energy-time entangled photon pairs requires detectors with sub-ps timing jitter and resolution. By invoking conjugate frequency correlations, a protocol to approach this fundamental limit has been developed using current detectors and existing telecom networks. The conjugate nature of time and energy encoding means that one can trade frequency for temporal bits (and vice versa) to minimize the effect of channel distortion such as nonlinear frequency conversion and dispersion, in addition to optimizing over transmission rate and channel bandwidth.

Gaussian filtering functions are assumed with center frequencies v_i and bandwidth δv, which project onto spatial modes |ζ_i. This is summarized in the following relation, where the total filtering operator Ŝ is taken as a sum over the individual channel filters Ŝⁱ.

$\begin{array}{cc}\begin{array}{c}\hat{S}=\sum _i^{}{\hat{S}}_i\\ =\sum _i^{}{\int }_{-\infty }^{\infty }\exp \left[\frac{-{\left(\omega -v_i\right)}^2}{2\delta v^2}\right]\omega ,{\zeta }^i〉〈\omega ,o\omega \end{array}& \left(7\right)\end{array}$

The Franson interference derived hereinabove assumes lossless propagation through the interferometer. However, this assumption is not valid in photonic integrated chips or fiber networks. Loss can be taken into account by placing a beam splitter in the long path of the Franson, which couples the waveguide mode with a vacuum mode (see FIG. 4). Working in the Heisenberg construction, evolving the annihilation operator through the virtual loss beam splitter and the two Franson beam splitters. The matrix for beam splitters 1 and 2, which leave the third mode undisturbed is given by

$\begin{array}{cc}{\hat{U}}_i=\left(\begin{array}{ccc}\sqrt{r_i}& \sqrt{1-r_i}& 0\\ \sqrt{1-r_i}& -\sqrt{r_i}& 0\\ 0& 0& 1\end{array}\right)& \left(8\right)\end{array}$

where iε1,2. The virtual loss beam splitter is given by

$\begin{array}{cc}{\hat{U}}_L=\left(\begin{array}{ccc}1& 0& 0\\ 0& \sqrt{r_L}& \sqrt{1-r_L}\\ 0& \sqrt{1-r_L}& -\sqrt{r_L}\end{array}\right)& \left(9\right)\end{array}$

The resulting annihilation operators are then â_A(t_A)=C₁â(t)+C₂â(t−Δt) and â_B(t_B)=C₁â(t)+C₂â(t−Δt−δt), disregarding the vacuum term, which will not affect coincidence counting. C₁=√{square root over (r₁)}√{square root over (r₂)} and C₂=√{square root over (1−r₁)}√{square root over (1−r₂)}√{square root over (r_L)}. For maximum visibility, C₁=C₂, so

$\begin{array}{cc}\frac{\sqrt{r_1}\sqrt{r_2}}{\sqrt{1-r_1}\sqrt{1-r_{2}}}=\sqrt{r_L}& \left(10\right)\end{array}$

This is plotted. For, r₁=r₂=½, the visibility simplifies to

$\begin{array}{cc}V_{\mathrm{PIC}}=\frac{2\text{?}}{1+\text{?}}{}^{-\delta t^2/2{\sigma }_{\mathrm{cor}}^2}\text{?}\text{?}\text{indicates text missing or illegible when filed}& \left(11\right)\end{array}$

where τ_a is the photon lifetime in the waveguide due to loss.

Detector jitter refers to the added uncertainty in the photon detection time of some stimulus, purely a result of detector electronics. Superconducting nanowire single photon detectors and InGaAs APDs both exhibit jitter of roughly 30 to 40 ps. Detector jitter is modeled as a Gaussian projection, {circumflex over (σ)}_det=∫e^{−tx2/2σdet2}|tt+t_x|dt_z. The jitter profile is not truly Gaussian and can be quite asymmetric, however this model allows for first-order analysis. If this is operated on both Alice and Bob's photons, assuming Eq. (1),

$\begin{array}{cc}\text{?}\Psi 〉\propto {\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }\exp \left[\frac{-{\left(t_A+t_B\right)}^2}{4{\sigma }_{\det }^2+16{\sigma }_{\mathrm{cab}}^2}\right]\exp \left[\frac{-{\left(t_A-t_B\right)}^2}{4{\sigma }_{\det }^2+4{\sigma }_{\mathrm{cor}}^2}\right]t_A,t_B,o_A,o_B〉t_At_B\text{?}\text{indicates text missing or illegible when filed}& \left(12\right)\end{array}$

Since σ_cohσ_det, the most important effect of jitter is to increase the observed correlation time roughly from σ_cor to σ_det. This can have a significant effect on the mutual information between Alice and Bob if ∝_det is on the order of σ_bin, as shown in FIG. 3.

The case of a single eavesdropper measuring a single photon of the photon pair is considered herein. Eve's temporal measurement is a POVM that can be written as a Gaussian filtering function:

Ê _t=∫_−∞^∞e^−t₂^{/2(σcohE′)2}|t|dt  (13)

Following, the amplitude function ψ(t_A,t_B)α exp[−(t_A−t_B)²/4σ_cor²]exp[−t_A²/4σ_coh²] for σ_coh>>σ_cor. Therefore

$\begin{array}{cc}{\Psi }_E〉=\text{?}\Psi 〉\propto {\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }\exp \left[-t_A^2\left(\frac{1}{4{\sigma }_{\mathrm{cab}}^2}+\frac{1}{4{\left(\text{?}\right)}^2}\right)\right]\exp \left[\frac{-{\left(t_A-t_B\right)}^2}{\text{?}}\right]t_A,t_B,o_A,o_B〉t_At_B\text{?}\text{indicates text missing or illegible when filed}& \left(14\right)\end{array}$

so the coherence time of the biphoton packet is strongly influenced by Eve's filtering bandwidth for σ_coh^E>>σ_coh, which gives the bound on her timing information.

Similarly, a weak frequency POVM is defined,

$\begin{array}{cc}{\hat{E}}_{\omega }={\int }_{-\infty }^{\infty }\text{?}\omega 〉〈\omega \omega \text{?}\text{indicates text missing or illegible when filed}& \left(15\right)\end{array}$

For 1/σ_cor>>1/σ_coh, and ignoring the spatial modes, |Ψ_F can be written as

|Ψ_F∫∫exp[−σ_cor²/4(2ω_A−ω_p)²]exp[−σ_coh²(ω_A+ω_B−ω_p)²]|ω_A,ω_Bdω_Adω_B.

Ê _ω|Ψ_F≈∫∫exp[−(σ_cor²/4+(σ_cor^E′)²/4)(2ω_A−ω_p)²]×exp[−σ_coh²(ω_A+ω_B−ω_p)²]|ω_A,ω_Bdω_Adω_B.  (17)

Thus, Eve projects the biphoton pair onto a narrower frequency distribution.

As will be apparent to those of skill in the art upon reading this disclosure, each of the individual embodiments described and illustrated herein has discrete components and features which may be readily separated from or combined with the features of any of the other several embodiments without departing from the scope or spirit of the present disclosed subject matter.

Claims

1.-9. (canceled)

10. An apparatus for receiving a quantum key, comprising: a channel configured to receive photons;a plurality of multi-channel filtering elements, each corresponding to a frequency and each configured to pass photons within a range of its corresponding frequency; anda plurality of photon detectors, each configured to receive photons passed by one of the plurality of multi-channel filtering elements and to send an indication of arrival time of photon detection for each frequency range.

11. The apparatus of claim 10, wherein the multi-channel filtering element comprises a grating.

12. The apparatus of claim 10, wherein the multi-channel filtering element comprises a dispersive element.

13. The apparatus of claim 10, further comprising a laser and a non linear crystal.

14. The apparatus of claim 13, wherein the non linear crystal generates pairs of photons by spontaneous parametric down conversion.

15. The apparatus of claim 14, wherein the pairs of photons are time-energy entangled photon pairs.

16. The apparatus of claim 10, wherein the photons are sent via a photonic integrated chip.

17. The apparatus of claim 10, wherein the photons are sent via optical fiber.

18. An apparatus for receiving a quantum key, comprising: a channel configured to receive photons;a plurality of multi-channel filtering elements, each corresponding to a frequency and each configured to pass photons within a range of its corresponding frequency; anda plurality of photon detectors, each configured to receive photons passed by one of the plurality of multi-channel filtering elements and to determine a time bin based on the arrival time of photon detection and frequency range.

19. The apparatus of claim 18, wherein the multi-channel filtering element comprises a grating.

20. The apparatus of claim 18, wherein the multi-channel filtering element comprises a dispersive element.

21. The apparatus of claim 18, further comprising a laser and a non linear crystal.

22. The apparatus of claim 21, wherein the non linear crystal generates pairs of photons by spontaneous parametric down conversion.

23. The apparatus of claim 22, wherein the pairs of photons are time-energy entangled photon pairs.

24. The apparatus of claim 18, wherein the photons are sent via a photonic integrated chip.

25. The apparatus of claim 18, wherein the photons are sent via optical fiber.

26. A communication system for distributing a quantum key, comprising: a laser and a non linear crystal adapted to generate pairs of photons by spontaneous parametric down conversion.a first apparatus comprising a first channel configured to receive photons;a first plurality of multi-channel filtering elements, each corresponding to a frequency and each configured to pass photons within a range of its corresponding frequency; anda first plurality of photon detectors, each configured to receive photons passed by one of the plurality of multi-channel filtering elements and to send an indication of arrival time of photon detection for each frequency range; anda second apparatus comprising a second channel configured to receive photons;a second plurality of multi-channel filtering elements, each corresponding to a frequency and each configured to pass photons within a range of its corresponding frequency; anda second plurality of photon detectors, each configured to receive photons passed by one of the plurality of multi-channel filtering elements and to send an indication of arrival time of photon detection for each frequency range; and