Research ArticleAPPLIED SCIENCES AND ENGINEERING

Quantum interference enables constant-time quantum information processing

See allHide authors and affiliations

Science Advances  19 Jul 2019:
Vol. 5, no. 7, eaau9674
DOI: 10.1126/sciadv.aau9674

Abstract

It is an open question how fast information processing can be performed and whether quantum effects can speed up the best existing solutions. Signal extraction, analysis, and compression in diagnostics, astronomy, chemistry, and broadcasting build on the discrete Fourier transform. It is implemented with the fast Fourier transform (FFT) algorithm that assumes a periodic input of specific lengths, which rarely holds true. A lesser-known transform, the Kravchuk-Fourier (KT), allows one to operate on finite strings of arbitrary length. It is of high demand in digital image processing and computer vision but features a prohibitive runtime. Here, we report a one-step computation of a fractional quantum KT. The quantum d-nary (qudit) architecture we use comprises only one gate and offers processing time independent of the input size. The gate may use a multiphoton Hong-Ou-Mandel effect. Existing quantum technologies may scale it up toward diverse applications.

INTRODUCTION

Science, medicine, and engineering demand efficient information processing. It is a long-standing goal to use quantum mechanics to substantially improve these computations (1). The processing routinely involves examining data as a function of complementary variables, e.g., time and frequency. This is done by the Fourier transform (FT) approximations, which accurately compute inputs of 2n samples in O(n 2n) steps (2). In the quantum domain, an analogous process exists, namely an FT of quantum amplitudes (3), which requires exponentially fewer O(n log n) quantum gates. Here, we report a quantum fractional Kravchuk-Fourier transform (KT), a related process suited to finite string processing (4). Unlike previous demonstrations (5, 6), our architecture involves only one gate, resulting in constant-time processing of quantum information. The gate exploits a generalized Hong-Ou-Mandel (HOM) effect (7), the basis for quantum-photonic information applications (8). We perform a proof-of-concept experiment by the creation of large photon number states, interfering them on a beam splitter (BS) and using photon-counting detection. The discrete FT (DFT) is an efficient approximation to the FT. The signal (x0, x1,…, xS) is taken to be samples of one period of a continuous function and is turned into a new sequence (X0, X1, …, XS) whereXk=1S+1∑l=0Se−i2πklS+1·xl,k=0,…,S(1)The DFT does not, however, reproduce all essential features of the FT. In some cases, a transform that is a fractional power of the FT, the α-fractional FT where 0 ≤ α ≤ 1, yields advantages (9). For α = 0, this transform is the identity, while for α = 1, this is the FT. If α=1n, where n = 2, 3, 4,…, then a composition of n α-fractional FTs amounts to the FT. This intuitive property does not hold true for the α-fractional DFT (Supplementary Materials), which generalizes the DFT, but for α = 1, it reduces to Eq. 1. This is because the α-fractional FT can be seen as a circular rotation of the signal in the time-frequency plane by angle πα2, while the α-fractional DFT is an elliptical rotation in this plane, which requires additional computation steps to properly approximate the α-fractional FT (9).

The DFT is powerful because of the fast FT (FFT) algorithm (2). Using an FFT lowers the number of operations from O(22n) to O(n 2n), which nevertheless remains a bottleneck in signal processing (10). The FFT uses a “divide and conquer” method to recursively split Eq. 1 into 2n sums, which can be processed quickly, and therefore is applicable to signals of period 2n. Notably, the minimal number of operations required to implement the DFT is unknown (11). The quantum FT, the cornerstone of quantum algorithms (12, 13), enables implementation of the DFT on quantum amplitudes with O(n log n) operations by processing n qubits (n quantum bits encode 2n amplitudes) (14).

In many applications, e.g., bioimaging, the signals are typically not periodic and are random in length. For these cases, the KT is a useful alternative to the FFT because it can be applied to finite signal processing (15, 16). The KT computes orthogonal moments corresponding to the Kravchuk polynomials, which are discrete and orthogonal with respect to a binomial distribution in the data space (4). By varying a parameter of the binomial distribution, one is able to set the fractionality α of the KT (Supplementary Materials). This feature allows the exploration of a specific region of interest of an image. To illustrate the action of a KT, the numerical study in fig. S5 in the Supplementary Materials demonstrates the advantages of KT over FFT in reconstructing test images.

The KT’s computational time is equal to the DFT’s runtime (17) (Supplementary Materials), and implementations with a lower number of operations are of high demand. Recently, quantum KTs (QKTs) have been realized in waveguides with two photons, but they are difficult to scale up and their fractionality is fixed by waveguide length (5, 6).

The α-fractional KT uses the weighted Kravchuk polynomials ϕk(p)(q,S) (4), which are real-valued and correspond to wave functions of finite harmonic oscillatorsXk=∑l=0Se−iπα2S2eiπ2(l−k)ϕk(p)(l−Sp,S)·xl,k=0,…,S(2)where p=sin2(πα4). Unlike plane waves, e−i2πklS+1, the polynomials are defined and orthogonal on a set of S + 1 points. This enables one to transform the signal as a finite sequence rather than as an infinite periodic one. In the limit of S → ∞, ϕk(p)(q,S) tend to eigenfunctions of quantum harmonic oscillators and the α-fractional KT reproduces the α-fractional FT. Equation 2 can be viewed in terms of overlaps of two spin-S2 states, in which they are prepared as eigenstates of S3 and one undergoes a rotation by angle πα2 generated by S1, eiπ2(l−k)ϕk(p)(l−Sp,S)=〈S2;S2−k∣eiπα2S1∣S2;S2−l〉.

RESULTS

Here, we demonstrate a single-step QKT with tunable fractionality using quantum effects based on multiparticle bosonic interference resulting from an exchange interaction. To this end, we interfere photon number states (light pulses with definite particle number) on a BS with an adjustable splitting ratio. This leads to a multiparticle HOM effect (18), which we observe using states with up to five photons. This QKT implementation enables constant-time quantum information processing for quantum d-nary (qudit) data encoding, which is set by the total number of interfering particles S, allowing up to d = S + 1 signal samples.

Photon number (Fock) states ∣l〉=(a†)ll∣0〉 and ∣S−l〉=(b†)S−lS−l∣0〉 impinging on a BS exhibit a generalized HOM effect (Fig. 1A). A BS interaction between two such inputs described by annihilation operators a and b is UBS=exp{θ2(a†be−iφ−ab†eiφ)}, where r=sin2θ2 is the BS reflectivity (defined as the probability of reflection of a single photon) and φ is the phase difference between the reflected and transmitted fields (19). Since φ does not influence our experiments, we assume that φ=π2 for convenience. If the BS is balanced (r = 0.5), then two photons at the input ports will leave through the same exit port. This is known as photon bunching (7). Similar effects hold for multiphoton number states (18). This is reflected in the probability amplitudes of detecting ∣k〉 and ∣S − k〉 behind the BS, AS(r)(k,l)=e−iθS2〈k,S−k∣UBS∣l,S−l〉. This is important for implementing the KT, since AS(r)(k,l)=e−iθS2eiπ2(l−k)·ϕk(r)(l−Sr,S); thus, if we send a quantum state ∣Ψ〉=∑l=0Sxl∣l,S−l〉 into the BS, then the probability of measuring k and S − k photons behind is the absolute square of a fractional QKT of the input probability amplitudes, ∣Xk∣2=∣∑l=0SAS(r)(k,l)·xl∣2 (compare Eq. 2). The reflectivity r determines the QKT fractionality, α=2θπ=4πarcsinr. Since two-mode optical interference can be achieved in a single step, regardless of the number of photons involved, this process implements a constant time QKT. For full derivations, see the Supplementary Materials.

Fig. 1 Photonic implementation of a fractional QKT.

(A) HOM interference of photon number states on a variable BS, followed by two photon-counting detectors, (B) Setup: Ti:Sa, titanium-sapphire laser pump (blue); BS, 50:50 BS; τ, optical phase delay; SPDC, periodically poled potassium titanyl phosphate nonlinear spontaneous parametric down-conversion waveguide chip that produces photon number–correlated states (red); VC, variable coupler; DAQ, data acquisition unit.

A deeper understanding of the result may be gained from the Schwinger representation of the spin algebra (Supplementary Materials), which links multiphoton interference to spin systems and allows the quantum states to be visualized on a Bloch sphere. In this picture, a total of S photons correspond to a spin-S2 system. The Hamiltonian generating UBS = exp { −iθHBS} corresponds to an Sx operator for a spin-S2. The two-mode Fock state ∣l, S − l〉 corresponds to an Sz=S2−l eigenstate, known as a Dicke state. Hence, HOM interference may be considered a rotation Rθ, φ = exp { −iθSx} of Sz around the Sx axis on the sphere. It transfers the eigenstate ∣S2;S2−l〉 to a superposition of Dicke states (Fig 2, A to D). The Q-function in Fig. 2 (E to H) shows that the initial and final states are eigenstates of two complementary observables, Sz and Sy, respectively. Thus, one may identify the former with a position, and the latter with a momentum eigenstate.

Fig. 2 HOM interference and QKT on a Bloch sphere.

(A to D) Two-mode Fock states (blue) correspond to Dicke states (black)—the basis of spin-S2 states. HOM interference turns Dicke states into a superposition of them. This coincides with a rotation Rθ,ϕ in the Dicke state basis. The two most distinct cases are shown: the rotation Rπ2,π2 of the pole ∣S2;S2⟩ and of the great circle state ∣S2;0〉. (E to H) Q-function representation of (A to D). HOM interference implements a rotation on the Bloch sphere by θ=π2 around Sx of input Sz-eigenbasis Dicke states and thus the full QKT (compare Eq. 2). The sequence (x0, x1,…, xS) is (1, 0, 0, …, 0) in (A) and (0, …, 1, …, 0) in (C). The QKT transfers the input—a position eigenstate—into the same state but in Sy basis—a momentum eigenstate.

We depicted experimental setup for multiphoton HOM interference in Fig. 1B. Two pulsed spontaneous parametric down-conversion (SPDC) sources each generated two-mode photon number–correlated states (Supplementary Materials). The signal and idler were separated with a polarization BS (PBS) into four spatial modes. The modes A and D were used for heralding and creation of Fock states ∣l〉 in B and ∣S − l〉 in C, which interfered in a variable ratio fiber coupler (the BS). An optical phase delay τ in one of the pump beams ensured optimal temporal overlap at interference. Photon number–resolved measurements were achieved using transition edge sensors (TESs) that we previously estimated to achieve over 90% efficiency (20).

We interfered the vacuum ∣0〉 (l = 0) with multiphoton Fock states ∣S〉 (S − l = S) on a coupler with splitting ratios r = 0.05 (green), 0.2 (red), 0.5 (blue), and 0.95 (gray) and measured photon number statistics. They are depicted in Fig. 3 (A to C) for S = 3, 4, 5. The input states encode sequences (x0 = 1, x1 = 0,…, xS = 0), while the measured probabilities set their QKTs to (∣X0∣2, ∣X1∣2, …, ∣XS∣2), where ∣Xk∣2=∣AS(r)(k,0)∣2. The reflectivities used correspond to fractionalities α = 0.28, 0.60, 1.00, and 1.72. Errors were estimated as a square root inverse of the number of measurements (Supplementary Materials). The second-order interferometric visibility reached values between 71.4 and 98.6% for S = 5 (Supplementary Materials).

Fig. 3 Photon number statistics resulting from Fock state |l, S − l〉 interference.

The probabilities of detecting ∣k〉 and ∣S − k〉 photons behind the BS for input (A) ∣0,3〉, (B) ∣0,4〉, (C) ∣0,5〉, (D) ∣1,2〉, (E) ∣2,2〉, and (F) ∣2,3〉. The BS reflectivities are r = 0.05 (green), 0.2 (red), 0.5 (blue), and 0.95 (gray). Vertical bars represent theoretical values for an ideal system, while dots are values determined in experiment. The states in (A) to (C) encode sequences (x0 = 1, x1 = 0, …, xS = 0), and states in (D) to (F) encode (0, 1, 0, 0), (0, 0, 1, 0, 0), and (0, 0, 1, 0, 0, 0), respectively. The measured probabilities set their QKTs (∣X0∣2, ∣X1∣2, …, ∣XS∣2), ∣Xk∣2=∣∑l=0SAS(r)(k,l)·xl∣2 of fractionality α = 0.28 (green), 0.60 (red), 1.00 (blue), and 1.72 (gray).

For the same values of r, we measured photon number distribution resulting from interference of ∣1,2〉, ∣2,2〉, and ∣2,3〉. They are shown in Fig. 3 (D to F). The inputs encode (0, 1, 0, 0), (0, 0, 1, 0, 0), and (0, 0, 1, 0, 0, 0), while ∣Xk∣2=∣A3(r)(k,1)∣2, ∣A4(r)(k,2)∣2, and ∣A5(r)(k,2)∣2. The visibility was between 54.8 and 99.5% (S = 5) (Supplementary Materials). Figure 3 shows that the theoretical values computed for an ideal system (the bars) match the experimental results (the dots) well.

DISCUSSION

Realization of the fractional QKT with qudit systems opens a previously unidentified prospect for transformation of large data sequences in O(1) time. This is not possible with the implementations based on waveguides. Both cases are examples of a nonuniversal quantum computer optimized for one task, which is the basis for a variety of important applications (15). The photonic proof of concept is currently limited by the range of input states that can be prepared. However, deterministic creation of an arbitrary superposition of Fock states has been demonstrated for trapped ions and superconducting resonators (21). Since a BS sees orthogonal spectral or polarization modes independently, one can extend the transform to higher dimensions (22, 23). We note that the QKT could also be implemented on existing quantum annealing processors (24), which operate on a chain of interacting spin-12 systems (Supplementary Materials), and using HOM interference of fermions with a symmetric wave function of the interfering degrees of freedom.

Large-scale realizations of the QKT may use an increasing number of measurements as an error minimization strategy (Supplementary Materials). This is akin to the common approach in classical data processing where the accuracy can be improved by an enhanced precision of numeric data types and number of iterations, without altering the proper algorithm. Errors resulting from losses in the system, e.g., at the BS, are easily corrected by applying a post-selection scheme and filtering out the cases when the total number of photons behind the BS is lower than in the input state.

O(1) computation of the fractional FT for continuous variable systems can be implemented with a shallow system, realizing a phase shift operation (25). This is a quantum counterpart of the operation of a single focusing lens in classical optics, which produces the fractional FT of the image placed at its focal length (9). The QKT operates on discrete variables, but when the sampling rate and the sequence length of input data increase, the α-fractional QKT tends to the α-fractional FT (4). This relation is also reflected by the fact that symmetric Kravchuk functions ϕk(p=1/2) (eigenfunctions of QKT) tend to Hermite-Gauss polynomials (eigenfunctions of the FT) in this limit.

Our result, along with the fact that qudit-based algorithms exhibit much lower number of operations than qubit-based ones (26), motivates the further development of highly controllable quantum harmonic oscillator platforms with implications for quantum signal processing in a whole range of applications. Provided efficient input state preparation and detection of larger Fock states, the O(1) QKT demonstrated here, in principle, may find practical applications in imaging of unprecedented quality, fostering early diagnostics and advances in neuroscience (27).

MATERIALS AND METHODS

A light pulse from a Ti:Sapphire laser at 775 nm [full width at half maximum (FWHM), 2 nm; repetition rate, 75 kHz] pumped collinear type-II phase-matched 8-mm-long SPDC waveguides written in a periodically poled potassium titanyl phosphate crystal sample. They generate two independent photon number–correlated states, the two-mode squeezed vacua ∣Ψ〉=∑n=0∞λn∣n,n〉, where λn=tanhngcoshg is a probability amplitude for creation of a pair of n photons and g is the parametric gain. The average photon number in the signal and idler mode equals 〈n^〉=sinh2g. For small g, cosh g ≈ 1, and thus, λn2≈sinh2ng=〈n^〉n. In the experiment, the average photon number was 〈n^〉≈0.2. This value is sufficient to ensure the emission of multiphoton pairs but, at the same time, to diminish the interferometric visibility of two-photon events. In both output states, the signal and idler pulses were split with a PBS to four spatial modes A to D. Subsequently, they were filtered by bandpass filters with an FWHM of 3 nm angle-tuned to the central wavelength of their respective spectra to reduce the broadband background typically generated in dielectric nonlinear waveguides (28). The pump beam was discarded with an edge filter. The modes A and D were used for heralding and conditional creation of Fock states in modes B and C, which interfered in a variable ratio polarization-maintaining (PM) fiber coupler. The coupling ratio can be set in the range of 0 to 100% with an error of ±1.5%. The heralding signal modes (H-pol.) were centered at 1554 nm, while the interfering idler modes (V-pol.) were at 1546 nm. We used TESs running at 70 mK, which allowed for photon number–resolved measurements in all modes (29). Their voltage output was captured with an analog to digital converter (ADC) card.

Before demonstrating the HOM interference, we characterized the setup. A high–photon number resolution and single-mode input states are pivotal for this experiment. The resolution of TES detectors (the confidence that the detector gives a correct information about the number of photons) was previously confirmed to exceed 95% (20). The depth of the HOM dip of 85.9 ± 0.3% for a two-photon interference indicated an effective Schmidt-mode number K = 1.16. For the measured 4-tuples of photon numbers, losses were computed by assuming perfect setup components, each followed by a BS with a reflection coefficient introducing the loss. We estimated that the total transmission in each mode to be approximately 50%. For the details, see the Supplementary Materials.

Measurements for individual settings of the splitting ratio were taken over approximately 400 s, giving 109 data samples for each r ranging from 0 to 1 with a step of approximately 3%. Small error bars for low photon numbers and larger bars for the higher ones result from keeping the pump power fixed and near-single modeness of the interfering beams.

SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/5/7/eaau9674/DC1

Fig. S1. Symmetric Kravchuk polynomials kn(1/2)(x, N) and functions ϕn(1/2)(x, N).

Fig. S2. Basis states for a 16-point KT.

Fig. S3. Basis states for a 16-point discrete FT.

Fig. S4. KT versus DFT.

Fig. S5. Example of FFT and KT image processing.

Fig. S6. HOM dip.

Fig. S7. Photon number statistics resulting from Fock state ∣l, S − l〉 interference.

Table S1. Second-order interferometric visibilities in HOM interference.

References (30–40)

This is an open-access article distributed under the terms of the Creative Commons Attribution license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.

REFERENCES AND NOTES

  1. ↵
    1. P. Zoller,
    2. T. Beth,
    3. D. Binosi,
    4. R. Blatt,
    5. H. Briegel,
    6. D. Bruss,
    7. T. Calarco,
    8. J. I. Cirac,
    9. D. Deutsch,
    10. J. Eisert,
    11. A. Ekert,
    12. C. Fabre,
    13. N. Gisin,
    14. P. Grangiere,
    15. M. Grassl,
    16. S. Haroche,
    17. A. Imamoglu,
    18. A. Karlson,
    19. J. Kempe,
    20. L. Kouwenhoven,
    21. S. Kröll,
    22. G. Leuchs,
    23. M. Lewenstein,
    24. D. Loss,
    25. N. Lütkenhaus,
    26. S. Massar,
    27. J. E. Mooij,
    28. M. B. Plenio,
    29. E. Polzik,
    30. S. Popescu,
    31. G. Rempe,
    32. A. Sergienko,
    33. D. Suter,
    34. J. Twamley,
    35. G. Wendin,
    36. R. Werner,
    37. A. Winter,
    38. J. Wrachtrup,
    39. A. Zeilinger
    , Quantum information processing and communication: Strategic report on current status, visions and goals for research in Europe. Eur. Phys. J. D 36, 203–228 (2005).
    OpenUrl
  2. ↵
    E. Brigham, Fast Fourier Transform and Its Applications (Prentice Hall, 1988).
  3. ↵
    M. A. Nielsen, I. L. Chuang, in Quantum Computation and Quantum Information (Cambridge University Press, 2000), pp. 216–242.
  4. ↵
    1. N. M. Atakishiyev,
    2. K. B. Wolf
    , Fractional Fourier–Kravchuk transform. J. Opt. Soc. Am. A 14, 1467–1477 (1997).
    OpenUrl
  5. ↵
    1. S. Weimann,
    2. A. Perez-Leija,
    3. M. Lebugle,
    4. R. Keil,
    5. M. Tichy,
    6. M. Gräfe,
    7. R. Heilmann,
    8. S. Nolte,
    9. H. Moya-Cessa,
    10. G. Weihs,
    11. D. N. Christodoulides,
    12. A. Szameit
    , Implementation of quantum and classical discrete fractional Fourier transforms. Nat. Commun. 7, 11027 (2016).
    OpenUrl
  6. ↵
    1. A. Crespi,
    2. R. Osellame,
    3. R. Ramponi,
    4. M. Bentivegna,
    5. F. Flamini,
    6. N. Spagnolo,
    7. N. Viggianiello,
    8. L. Innocenti,
    9. P. Mataloni,
    10. F. Sciarrino
    , Suppression law of quantum states in a 3D photonic fast Fourier transform chip. Nat. Commun. 7, 10469 (2016).
    OpenUrl
  7. ↵
    1. C. K. Hong,
    2. Z. Y. Ou,
    3. L. Mandel
    , Measurement of subpicosecond time intervals between two photons by interference. Phys. Rev. Lett. 59, 2044–2046 (1987).
    OpenUrlCrossRefPubMedWeb of Science
  8. ↵
    1. K. Makino,
    2. Y. Hashimoto,
    3. J.-i. Yoshikawa,
    4. H. Ohdan,
    5. T. Toyama,
    6. P. van Loock,
    7. A. Furusawa
    , Synchronization of optical photons for quantum information processing. Sci. Adv. 2, e1501772 (2016).
    OpenUrlFREE Full Text
  9. ↵
    1. E. Sejdić,
    2. I. Djurović,
    3. L. Stanković
    , Fractional Fourier transform as a signal processing tool: An overview of recent developments. J. Signal Process. 91, 1351–1369 (2010).
    OpenUrl
  10. ↵
    1. J. Morgenstern
    , Note on a lower bound on the linear complexity of the fast Fourier transform. J. ACM 20, 305–306 (1973).
    OpenUrl
  11. ↵
    N. Ailon, Tighter fourier transform lower bounds, in International Colloquium on Automata, Languages, and Programming (Springer, 2015), pp. 14–25.
  12. ↵
    1. C. M. Bowden,
    2. G. Chen,
    3. Z. Diao,
    4. A. Klappenecker
    , The universality of the quantum Fourier transform in forming the basis of quantum computing algorithms. J. Math. Anal. Appl. 274, 69–80 (2002).
    OpenUrl
  13. ↵
    1. D. Gottesman
    , Fault-tolerant quantum computation with higher-dimensional systems. Chaos Solitons Fractals 10, 1749–1758 (1999).
    OpenUrl
  14. ↵
    L. Hales, S. Hallgren, An improved quantum Fourier transform algorithm and applications, in Proceedings of the 41st IEEE Annual Symposium on Foundations of Computer Science (IEEE, 2000), pp. 515–525.
  15. ↵
    1. P.-T. Yap,
    2. R. Paramesran,
    3. S.-H. Ong
    , Image analysis by Krawtchouk moments. IEEE Trans. Image Process. 12, 1367–1377 (2003).
    OpenUrlPubMed
  16. ↵
    1. A. Kumar
    , Nonlocal means image denoising using orthogonal moments. Appl. Opt. 54, 8156–8165 (2015).
    OpenUrl
  17. ↵
    A. Venkataramana, P. A. Raj, Recursive computation of forward Krawtchouk moment transform using Clenshaw's recurrence formula, in 2011 Third IEEE National Conference on Computer Vision, Pattern Recognition, Image Processing and Graphics (NCVPRIPG) (IEEE, 2011), pp. 200–203.
  18. ↵
    1. R. A. Campos,
    2. B. E. A. Saleh,
    3. M. C. Teich
    , Quantum-mechanical lossless beam splitter: SU(2) symmetry and photon statistics. Phys. Rev. A 40, 1371–1384 (1989).
    OpenUrlCrossRefPubMed
  19. ↵
    1. M. S. Kim,
    2. W. Son,
    3. V. Bužek,
    4. P. L. Knight
    , Entanglement by a beam splitter: Nonclassicality as a prerequisite for entanglement. Phys. Rev. A 65, 032323 (2002).
    OpenUrl
  20. ↵
    1. P. C. Humphreys,
    2. B. J. Metcalf,
    3. T. Gerrits,
    4. T. Hiemstra,
    5. A. E. Lita,
    6. J. Nunn,
    7. S. W. Nam,
    8. A. Datta,
    9. W. S. Kolthammer,
    10. I. A. Walmsley
    , Tomography of photon-number resolving continuous-output detectors. New J. Phys. 17, 103044 (2015).
    OpenUrl
  21. ↵
    1. M. Hofheinz,
    2. H. Wang,
    3. M. Ansmann,
    4. R. C. Bialczak,
    5. E. Lucero,
    6. M. Neeley,
    7. A. D. O’Connell,
    8. D. Sank,
    9. J. Wenner,
    10. J. M. Martinis,
    11. A. N. Cleland
    , Synthesizing arbitrary quantum states in a superconducting resonator. Nature 459, 546–549 (2009).
    OpenUrlCrossRefPubMedWeb of Science
  22. ↵
    1. M. G. Raymer,
    2. S. J. van Enk,
    3. C. J. McKinstrie,
    4. H. J. McGuinness
    , Interference of two photons of different color. Opt. Commun. 283, 747–752 (2009).
    OpenUrl
  23. ↵
    1. T. Kobayashi,
    2. R. Ikuta,
    3. S. Yasui,
    4. S. Miki,
    5. T. Yamashita,
    6. H. Terai,
    7. T. Yamamoto,
    8. M. Koashi,
    9. N. Imoto
    , Frequency-domain Hong–Ou–Mandel interference. Nat. Photonics 10, 441–444 (2016).
    OpenUrl
  24. ↵
    D. Dahl, Qubits, Couplers & Quantum Computing in 2017, in ISC High Performance 2017 Conference (2017).
  25. ↵
    1. H.-y. Fan,
    2. L.-y. Hu
    , Correspondence between quantum-optical transform and classical-optical transform explored by developing Dirac’s symbolic method. Front. Phys. 7, 261–310 (2012).
    OpenUrl
  26. ↵
    A. Pavlidis, E. Floratos, Arithmetic circuits for multilevel qudits based on quantum fourier transform. arXiv:1707.08834 [quant-ph] (27 July 2017).
  27. ↵
    1. D. Cyranoski
    , China launches brain-imaging factory. Nature 548, 268 (2017).
    OpenUrl
  28. ↵
    1. A. Eckstein,
    2. A. Christ,
    3. P. J. Mosley,
    4. C. Silberhorn
    , Highly efficient single-pass source of pulsed single-mode twin beams of light. Phys. Rev. Lett. 106, 013603 (2011).
    OpenUrlPubMed
  29. ↵
    1. T. Gerrits,
    2. N. Thomas-Peter,
    3. J. C. Gates,
    4. A. E. Lita,
    5. B. J. Metcalf,
    6. B. Calkins,
    7. N. A. Tomlin,
    8. A. E. Fox,
    9. A. Lamas Linares,
    10. J. B. Spring,
    11. N. K. Langford,
    12. R. P. Mirin,
    13. P. G. R. Smith,
    14. I. A. Walmsley,
    15. S. W. Nam
    , On-chip, photon-number-resolving, telecommunication-band detectors for scalable photonic information processing. Phys. Rev. A 84, 060301 (2011).
    OpenUrlCrossRef
  30. ↵
    P. Feinsilver, Philip, J. Kocik, Krawtchouk polynomials and Krawtchouk matrices, in Recent Advances in Applied Probability (Springer, 2005), pp. 115–141.
  31. G. A. Papakostas, D. E. Koulouriotis, E. G. Karakasis, V. D. Tourassis, A General Framework for Computation of Biomedical Image Moments, in Biomedical Engineering, Trends in Electronics, Communications and Software (INTECH, 2011).
  32. G. Gautam, K. Choudhary, S. Chatterjee, M. H. Kolekar, Facial expression recognition using Krawtchouk moments and support vector machine classifier, in 2017 Fourth IEEE International Conference on Image Information Processing (ICIIP) (IEEE, 2017) pp. 1–6.
    1. X. Liu,
    2. G. Han,
    3. J. Wu,
    4. Z. Shao,
    5. G. Coatrieux,
    6. H. Shu
    , Fractional Krawtchouk transform with an application to image watermarking. IEEE Trans. Signal Process. 65, 1894–1908 (2017).
    OpenUrl
  34. A. Mademlis, A. Axenopoulos, P. Daras, D. Tzovaras, M. G. Strintzis, 3D content-based search based on 3D Krawtchouk moments, in Third IEEE International Symposium on 3D Data Processing, Visualization, and Transmission (3DPVT'06) (IEEE, 2006), pp. 743–749.
  35. G. A. Papakostas, B. G. Mertzios, D. A. Karras, Performance of the orthogonal moments in reconstructing biomedical images, in 2009 16th International Conference on Systems, Signals and Image Processing (IEEE, 2009), pp. 1–4.
    1. N. M. Atakishiyev,
    2. K. B. Wolf
    , Approximation on a finite set of points through Kravchuk functions. Rev. Mex. Fis. 40, 366–377 (1994).
    OpenUrl
    1. D. N. Klyshko
    , Use of two-photon light for absolute calibration of photoelectric detectors. Sov. J. Quantum Electron. 10, 1112–1117 (1980).
    OpenUrl
    1. A. Christ,
    2. K. Laiho,
    3. A. Eckstein,
    4. K. N. Cassemiro,
    5. C. Silberhorn
    , Probing multimode squeezing with correlation functions. New J. Phys. 13, 033027 (2011).
    OpenUrl
    1. R. Filip
    , Overlap and entanglement-witness measurements. Phys. Rev. A 65, 062320 (2002).
    OpenUrl
  40. ↵
    1. M. Christandl,
    2. N. Datta,
    3. A. Ekert,
    4. A. J. Landahl
    , Perfect state transfer in quantum spin networks. Phys. Rev. Lett. 92, 187902 (2004).
    OpenUrlPubMed
Acknowledgments: Funding: M.S. and A.B. were supported by the Foundation for Polish Science “First Team” project no. POIR.04.04.00-00-220E/16-00 (originally, FIRST TEAM/2016-2/17), the National Science Centre (NCN) grant no. 2012/04/M/ST2/00789, and MNiSW Iuventus Plus project no. IP 2014 044873. A.E., I.A.W., and M.S. were supported by the Engineering and Physical Sciences Research Council project no. EP/K034480/1. Numerical computations were performed with a Zeus cluster in the ACK “Cyfronet” AGH computer center. This work was a contribution of NIST, an agency of the U.S. government, not subject to copyright. Author contributions: M.S. and A.B. developed the theory, while A.E., M.M., W.R.C., W.S.K., and I.A.W. were responsible for the realization of the experiment. J.J.R., S.W.N., T.G., and A.L. delivered and maintained the TES detection system. A.B and A.E. developed the software and performed the numerical computations. A.B. prepared the plots. All of the co-authors wrote the manuscript. Competing interests: M.S., A.B., A.E., and I.A.W. are inventors on a patent application related to this work filed by the University of Warsaw (no. P.426228, filed on 6 July 2018). The authors declare that they no other competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data available related to this paper may be requested from the authors.
View this article with LENS

Article Tools

My saved folders

Stay Connected to Science Advances

Related Content

Similar Articles in:

Citing Articles in:

Navigate This Article