Abstract
The probabilistic nature of singlephoton sources and photon–photon interactions encourages encoding as much quantum information as possible in every photon for the purpose of photonic quantum information processing. Here, by encoding highdimensional units of information (qudits) in time and frequency degrees of freedom using onchip sources, we report deterministic twoqudit gates in a single photon with fidelities exceeding 0.90 in the computational basis. Constructing a twoqudit modulo SUM gate, we generate and measure a singlephoton state with nonseparability between time and frequency qudits. We then employ this SUM operation on two frequencybin entangled photons—each carrying two 32dimensional qudits—to realize a fourparty highdimensional Greenberger–Horne–Zeilinger state, occupying a Hilbert space equivalent to that of 20 qubits. Although highdimensional coding alone is ultimately not scalable for universal quantum computing, our design shows the potential of deterministic optical quantum operations in large encoding spaces for practical and compact quantum information processing protocols.
Introduction
Quantum information processing has drawn massive attention due to its power in solving some crucial algorithms exponentially faster than their classical counterparts,^{1} as well as its ability to transmit information in a fully secure fashion, two capabilities looked to be combined in the emerging quantum internet.^{2} Among the platforms that can exhibit quantum behavior, optical states have the advantages of low decoherence and suitability for longdistance communications, yet the weak coupling of photons to their surroundings also makes it extremely difficult to manipulate the state of one photon based on the state of another. This operation, needed for a twoqubit gate, is probabilistic with standard linear optics and photon counting.^{3} Quantum gates have been demonstrated in a number of different photonic degrees of freedom, such as polarization,^{4} orbital angular momentum,^{5} time,^{6} and frequency,^{7,8} and to sidestep the challenges of probabilistic multiphoton interactions, encoding qubits in different degrees of freedom (DoFs) in a single photon has been demonstrated, where each DoF carries one qubit and, now, operations between different qubits can be made deterministic.^{9,10} This scheme allows encoding more quantum information in single photons, and can find use in standalone processing tasks or be subsequently incorporated into larger systems built on true photon–photon interactions, thus offering a potentially more efficient method for photonic quantum information processing. Even though in this case two and threequbit operations can be executed with unity success probability, each DoF contains only one qubit, and the number of a photon’s DoFs are limited; thus the size of the Hilbert space in which these deterministic transformations can happen is fairly moderate (e.g., an eightdimensional Hilbert space has been demonstrated by encoding three qubits in three different DoFs of a single photon^{10}).
In this article, we take advantage of the high dimensionality in two particular DoFs of a single photon—namely, time and frequency, which are both compatible with fiber optical transmission—to encode one qudit in each DoF. We consider multiple time bins and frequency bins; as long as the frequency spacing between different modes (Δf) and the timebin spacing (Δt) are chosen such that they exceed the Fourier transform limit (i.e., Δf Δt > 1), we are able to manipulate the time and frequency DoFs independently in a hyperencoding fashion, using concepts developed in timedivision and wavelengthdivision multiplexing, respectively.^{11,12} In other words, each time–frequency mode pair constitutes a welldefined entity, or plaquette,^{11,12} which is sufficiently separated from its neighbors to provide stable encoding (Fig. 1). Alternatively, this can be understood by considering bandwidthlimited plaquettes with individual spectral linewidth δf (corresponding to temporal duration ~1/δf). These will not overlap in time–frequency space as long as the chosen bin separations satisfy Δf > δf and Δt > 1/δf. Combined, then, these two equations yield the aforementioned condition Δf Δt > 1. An analogous process is at work in the advanced optical modulation formats gaining adoption in modern digital communications, where many bits are encoded in a single symbol via modulation of canonically conjugate quadratures.^{13} Since our single photons can potentially be generated in a superposition of many time and frequency bins, multiple qubits can be encoded in each DoF, making our proposed scheme a favorable platform for deterministic optical quantum information processing on Hilbert spaces dramatically larger than previously demonstrated deterministic qubitbased gates. Ultimately, the total number of DoFs carried by a single photon is limited, so one cannot increase the Hilbert space indefinitely by encoding in increasingly more properties within individual photons. The Hilbert space can be increased, though, by expanding the dimensionality within each DoF. While enabling only linear scaling of the Hilbert space with the number of modes,^{14} and thereby not facilitating the exponential scaling required for faulttolerant quantum computing, qudit encoding promises significant potential in the current generation of quantum circuits. It has been shown, for example, that twoqudit optical gates are useful in transmitting quantum states with higher information content per photon by means of qudit teleportation,^{15} a task that requires twoqudit gates which can operate on the different degrees of freedom of a single photon^{16,17}—precisely the functionality we demonstrate here.
Results
To enable the realization of all singlequdit unitaries, it is sufficient to demonstrate the generalized Pauli gates X (cyclic shift) and Z (statedependent phase), which are universal for singlequdit operations,^{5} and from which all ddimensional Weyl operators can be constructed.^{18} The Z gate applies a unique phase shift to each of the d basis states, which can be easily executed with a phase modulator and a pulse shaper in the time domain and frequency domain, respectively. Specifically, for the basis state n〉 (n = 0, 1, …, d − 1), we have \({\mathrm{Z}}\left n \right\rangle = \exp (2{\mathrm{\pi i}}n{\mathrm{/}}d)\left n \right\rangle\) and Xn〉 = n ⊕ 1〉, where ⊕ denotes addition modulo d. An X gate in the frequency domain can be realized using a Z gate sandwiched between two highdimensional DFT gates. Such a DFT operation has been recently demonstrated,^{7} completing in principle the universal gate set for singlequdit frequencydomain operations. To complete the gate set in the time domain, we demonstrate the timebin X gate presented in Fig. 2a, operating on time bins in three dimensions, a process which corresponds to statedependent delay. Because the gate operates on each photon individually, we can simulate its performance with coherent states; the statistics of the input field have no impact on the principle of operation. Of course, to apply this gate in multiphoton quantum information processing, true single photons would need to be tested as well, the preparation or heralding of which is technically demanding and could introduce additional noise. However, as this noise is extrinsic to the gate itself, we focus on weak coherent states for initial characterization here. To test for the correct modal operation of this gate, we use a continuouswave (CW) laser and prepare the desired weak coherent state by carving out three time bins {0〉_{t}, 1〉_{t}, 2〉_{t}} using an intensity modulator and manipulating their relative phases with a phase modulator. The time bins are 3ns wide with Δt = 6ns centertocenter spacing. To perform the X operation, we need to separate the time bins 0〉_{t} and 1〉_{t} from 2〉_{t} and delay the route for time bins 0〉_{t} and 1〉_{t} by 3 bins (18 ns). We realize the necessary spatial separation between the time bins with a Mach–Zehnder modulator (MZM) switch. We emphasize that while most MZM designs are oneport devices, with one of the two output paths terminated, this 1 × 2 version permits access to both interferometer outputs, and accordingly it is in principle lossless—as required for a unitary operation. (In practice, of course, insertion loss reduces throughput, but it should be possible in the future to significantly reduce this loss through, e.g., onchip integration.) After the pathdependent delay, another 1 × 2 MZM, but operated in reverse, can be used to recombine the time bins deterministically as well. However, due to lack of equipment availability, in this proofofprinciple experiment we employ a 2 × 2 fiber coupler for recombination, which introduces an additional 3dB power penalty. To measure the gate output, we synchronize a singlephoton detector and time interval analyzer with the generated time bins. The transformation matrix performed by the X gate when probed by single time bins yields a computational basis fidelity \({\cal{F}}_{\mathrm{C}}\) of 0.996 ± 0.001, shown in Fig. 2b (see the Methods section). As such computationalbasisonly measurements do not reflect the phase coherence of the operation, we next prepare superposition states as input and interfere the transformed time bins after the gate with a cascade of 1bin and 2bin delay unbalanced interferometers. In order to combat environmentally induced phase fluctuations, we stabilize both these interferometers and the X gate by sending a CW laser in the backwards direction and using a feedback phase control loop. We apply a phase of 0, ϕ, and 2ϕ to the timebins 0〉_{t}, 1〉_{t}, and 2〉_{t}, respectively, with the phase modulator in the state preparation stage and sweep ϕ from 0 to 2π, obtaining the interference pattern shown in Fig. 2c. After subtraction of the background, we calculate a visibility of 0.94 ± 0.01 from the maximum and minimum points, showing strong phase coherence (the ability to preserve and utilize coherent superpositions) between the time bins after the gate. If for concreteness we assume a channel model consisting of pure depolarizing (white) noise,^{18} we can use this visibility to estimate the process fidelity \({\cal{F}}_{\mathrm{P}}\), finding \({\cal{F}}_{\mathrm{P}} = 0.92 \pm 0.01\) for the X gate (see the Methods section). Given the ability to perform arbitrary onequdit operations using combinations of X and Z gates, it follows that it is in principle possible to generate and measure photons in all mutually unbiased bases^{19}—an essential capability for highdimensional quantum key distribution (QKD),^{20} which has been proven to offer greater robustness to noise compared with qubitbased QKD^{21} and can enable significantly higher secret key rates over metropolitanscale distances.^{22}
With this highperformance timebin X gate in hand, we are then in a position to incorporate it into a frequency network to realize deterministic twoqudit gates, where the frequency DoF acts as the control and the time DoF is the target qudit. For this demonstration, instead of a weak coherent state, we utilize true single photons, heralded by detecting the partner photon of a frequencybin entangled pair generated through spontaneous fourwave mixing in an onchip silicon nitride microresonator. The time bins, defined by intensity modulation of the pump, couple into a microring resonator with a free spectral range (FSR) Δf = 380 GHz and resonance linewidths δf≃250 MHz, generating a biphoton frequency comb. The timebin and frequencybin entanglement of such sources have been proven recently.^{23,24,25,26} As our time and frequencybins exceed the Fourier limit (Δf Δt = 2280, δf Δt = 1.5), our time–frequency entangled photons can be considered hyperentangled—that is, entangled in two fully separable DoFs. The signal and idler photons from the first three combline pairs are then selected and separated with a commercial pulse shaper, as shown in Fig. 3a. Now that the time bins and frequency bins are all generated in the state preparation stage, the idler photons are sent to a singlephoton detector for heralding, and the signal photons are what carry the two qudits in the three time bins {0〉_{t}, 1〉_{t}, 2〉_{t}} and frequency bins {0〉_{f}, 1〉_{f}, 2〉_{f}}. This procedure lets us prepare any timebin/frequencybin product state m〉_{t}n〉_{f} (m, n = 0, 1, 2) of the full computational basis set. In principle, we could also herald arbitrary time–frequency superposition states in this setup, by first sending the idler photon through a combination of time or frequencybin interferometers prior to detection in the temporal and spectral eigenbases. This more general case would permit the preparation of any twoqudit state and is an important area for further research.
As the first twoqudit gate, we demonstrate the controlledincrement (CINC) operation, where an X gate is applied to the timebin qudit only when the frequency qudit is in the state 2〉_{f}. This twoqudit gate along with arbitrary singlequdit gates [which, as noted above, can be formed from singlequdit X and Z operations^{5}] complete a universal set for any quantum operation.^{27} To implement this gate, we separate 2〉_{f} from the other two frequency bins with a dense wavelengthdivision multiplexing (DWDM) filter and route it to a timebin X gate (Fig. 3a); no operation happens on the route of the other two frequency bins. The frequency bins are then brought back together with another DWDM with zero relative delay to complete the twoqudit gate operation. To measure the transformation matrix of this gate in the computational basis, we prepare the input state in each of the nine combinations of single time bins and frequency bins, using the first intensity modulator and the pulse shaper, respectively. We then record the signal counts in all possible output timebin/frequencybin pairs, conditioned on detection of a particular idler time–frequency mode, by inserting three different DWDMs in the path of the signal photons to pick different frequency bins. The measured transformation matrix is shown in Fig. 3b, with accidentalsubtracted fidelity \({\cal{F}}_{\mathrm{C}} = 0.90 \pm 0.01\) (see the Methods section).
For the next step, we implement an even more complex operation, the SUM gate—a generalized controlledNOT gate^{28}—which adds the value of the control qudit to the value of the target qudit, modulo 3. In this gate, the time bins associated with 0〉_{f} are not delayed, the time bins associated with 1〉_{f} experience a cyclic shift by 1 slot, and the time bins corresponding to 2〉_{f} go through a cyclic shift of 2 slots. To delay the time bins dependent on their frequencies, we induce a dispersion of −2 ns.nm^{−1} on the photons using a chirped fiber Bragg grating (CFBG); this imparts 6 ns (1bin) and 12 ns (2bin) delays for the temporal modes of 1〉_{f} and 2〉_{f}, respectively, as required for the SUM operation. However, this delay is linear—not cyclic—so that some of the time bins are pushed outside of the computational space, to modes 3〉_{t} and 4〉_{t}. Returning these bins to overlap with the necessary 0〉_{t} and 1〉_{t} slots can be achieved using principles identical to the timebin X gate with a relative delay of three bins. The experimental setup is shown in Fig. 3a, where we use the same techniques as for the CINC gate to measure the transfer matrix shown in Fig. 3c, with \({\cal{F}}_{\mathrm{C}} = 0.92 \pm 0.01\). The fact that this SUM gate is implemented with qudits in a single step potentially reduces the complexity and depth of quantum circuits in algorithms that require an addition operation.^{29} Lack of frequency shifting components in these gates can be confirmed by the small offdiagonal terms in Fig. 3b, c for which the input and output frequency bins differ.
To show the ability of our design to operate on large Hilbert spaces, we extend the dimensions of our qudits and encode two 16dimensional quantum states in the time and frequency DoFs of a single photon. For this demonstration, as we want to use more time bins and set a smaller frequency spacing between modes, we use a periodically poled lithium niobate (PPLN) crystal as a broadband source of time–frequency entangled photons followed by a programmable pulse shaper to set the frequency spacing and linewidth, instead of a microring with fixed frequency spacing. (We note that, in principle, one could still use an integrated source for these experiments by appropriately engineering a microring’s FSR, bandwidth, and resonance linewidth to realize spectral and temporal spacings tighter than the integrated photon sources currently available to us.) In this experiment, we first shine a 773nm CW laser on the PPLN crystal, generating entangled photons with a bandwidth of ~5 THz.^{30} We then carve 16 time bins with a full width at half maximum of ~200 ps and 1.2ns spacing between them, to generate the timebin qudits. Then, a pulse shaper is used to carve out the frequency of these entangled photons to generate sixteen 22GHzwide frequency bins on both the signal and idler side of the spectrum, each spaced by 75 GHz from each other. Now that we have 16dimensional qudits in both time and frequency, we send a heralded signal photon into the same SUMgate structure. We note that after the CFBG, the individual time bins will spread to ~350 ps due to their now larger (22 GHz) linewidth. While not necessary in this proofofprinciple experiment, such spreading could be reduced either by using a smaller linewidth for our frequency modes (e.g., with a Fabry–Perot etalon), or by using a dispersive element with a stepwise frequencydependent delay profile.^{31,32,33} To verify the operation, we send in different input twoqudit states, chosen from one of 256 basis states, and measure the output after the gate. While this yields a total of 256 × 256 (2^{16}) computational input/output combinations to test, we have no active frequencyshifting elements in the SUM gate to shift 75 GHzspaced frequencies into each other, so we make the reasonable assumption that the frequency qudit remains unchanged through the operation. This is also enforced by the high extinction ratio of the pulse shaper (~40 dB), which blocks unwanted frequency bins. This allows us to focus on results in the sixteen 16 × 16 transfer matrices measured in Fig. 4a–p (a subset with a total of 2^{12} input/output combinations). In each matrix, 16 different inputs with the same frequency and different time bins are sent into the SUM gate, and the output time bins are measured. For this experiment, we use superconducting nanowire singlephoton detectors (SNSPDs), which allow us to report our data without accidental subtraction. The average computational space fidelity for the whole process, with the assumption that frequencies do not leak into each other, can be calculated as \(\bar {\cal{F}}_{\mathrm{C}} = 0.9589 \pm 0.0005\), which shows the high performance of our operation. This high fidelity benefits greatly from the high extinction ratio of the intensity modulator used to carve the time bins (~25 dB). To show the coherence of our SUM gate, we use this setup to perform a SUM operation on a threedimensional input state, \(\left \psi \right\rangle _{{\mathrm{in}}} = \frac{1}{{\sqrt 3 }}\left( {\left 0 \right\rangle _{\mathrm{f}} + \left 1 \right\rangle _{\mathrm{f}} + \left 2 \right\rangle _{\mathrm{f}}} \right)\left 0 \right\rangle _{\mathrm{t}}\), which results in a maximally nonseparable state^{34} between time and frequency DoFs: \(\left \psi \right\rangle _{{\mathrm{out}}} = \frac{1}{{\sqrt 3 }}\left( {\left {00} \right\rangle _{{\mathrm{ft}}} + \left {11} \right\rangle _{{\mathrm{ft}}} + \left {22} \right\rangle _{{\mathrm{ft}}}} \right)\). To quantify the dimensionality of this state, we use an entanglement certification measure called entanglement of formation (E_{of}).^{35,36} We experimentally obtain E_{of} ≥ 1.19 ± 0.12 ebits, where 1 ebit corresponds to a maximally nonseparable pair of qubits, while 1.585 ebits represents the maximum for two threedimensional parties (see Methods); in exceeding the qubit limit, our state thus possesses true highdimensional nonseparability.
One of the most crucial challenges toward optical quantum operations is the lack of ondemand photon sources. Therefore, it is interesting to consider our scheme for application to quantum communication and networking, for which operations with just a few qudits have potential impact. A gate very similar to the SUM gate is the XOR gate, which subtracts the control qudit from the target and is a requirement for qudit teleportation protocols.^{15,37,38} Since teleportation of quantum states is possible using different degrees of freedom of an entangled photon pair,^{165} a singlephoton twoqudit gate in our time–frequency paradigm could be applied directly for teleporting highdimensional states. Specifically, the XOR gate can be demonstrated by using positive dispersion and reconfiguring the switching in the SUM gate, or in the threedimensional case, by simply relabeling the frequency bins 0〉_{f} → 2〉_{f} and 2〉_{f} → 0〉_{f} and performing the same process as the SUM operation. In addition, these twoqudit gates can be used for the purpose of beating the channel capacity limit for standard superdense coding for highdimensional entangled states.^{39} In such quantum communications applications for the twoqudit gates, a modest number of state manipulations brings significant value.
The demonstrated SUM gate can also be used to produce highdimensional Greenberger–Horne–Zeilinger (GHZ) states.^{40} GHZ states consist of more than two parties, entangled with each other in a way that measurement of one party in the computational basis determines the state of all the other parties.^{41} It has been only recently that these states were demonstrated in more than two dimensions, where a threedimensional threeparty GHZ state was realized using the orbital angular momentum of optical states.^{40} Here, we take advantage of our SUM gate and the large dimensionality of time–frequency states to generate a fourparty GHZ state with 32 dimensions in each DoF. We start from the state \(\left \psi \right\rangle _{{\mathrm{in}}} = \frac{1}{{\sqrt {32} }}\left {0,0} \right\rangle _{{\mathrm{t}}_{\mathrm{s}}{\mathrm{t}}_{\mathrm{i}}}\mathop {\sum}\nolimits_{m = 0}^{31} {\left {m,m} \right\rangle _{{\mathrm{f}}_{\mathrm{s}}{\mathrm{f}}_{\mathrm{i}}}}\), which means both signal and idler photons are initialized in the zero timebin state and are maximally entangled in the frequency domain. Then, we operate deterministic SUM gates separately on both signal and idler photons, resulting in a fourparty GHZ state of the form \(\left \psi \right\rangle _{{\mathrm{out}}} = \frac{1}{{\sqrt {32} }}\mathop {\sum}\nolimits_{m = 0}^{31} {\left {m,m,m,m} \right\rangle _{{\mathrm{f}}_{\mathrm{s}}{\mathrm{t}}_{\mathrm{s}}{\mathrm{f}}_{\mathrm{i}}{\mathrm{t}}_{\mathrm{i}}}}\). Since the initial state only consists of the zeroth time bins, the dispersion module does not shift any of the bins outside of the computational space; hence the interferometric structure used in the full SUM gate is not required when operating within this subspace. The GHZ state is measured in the computational basis (Fig. 5); the plot contains coincidences for all basis states in the set \(\{ \left {m,n,k,l} \right\rangle _{{\mathrm{f}}_{\mathrm{s}}{\mathrm{t}}_{\mathrm{s}}{\mathrm{f}}_{\mathrm{i}}{\mathrm{t}}_{\mathrm{i}}};0 \le m,n,k,l \le 31\}\). Only states whose four qudits match (i.e., \(\left {m,m,m,m} \right\rangle _{{\mathrm{f}}_{\mathrm{s}}{\mathrm{t}}_{\mathrm{s}}{\mathrm{f}}_{\mathrm{i}}{\mathrm{t}}_{\mathrm{i}}}\)) have high counts, as expected for a GHZ state. Of course, full characterization of the state requires measurements in superposition bases as well,^{42} but due to the additional insertion loss associated with superposition measurements in time and frequency using interferometers and phase modulators, respectively, we were unable to measure such projections. Remarkably, the demonstrated GHZ state resides in a Hilbert space equivalent to that of 20 qubits, an impressive 1,048,576 (32^{4}) dimensions. We emphasize that the four parties of the demonstrated GHZ state are carried by only two photons, and hence cannot be used for genuine multipartite GHZ applications such as demonstration of Bell’s theorem without inequalities,^{41} quantum secret sharing,^{43} or opendestination teleportation.^{44} However, the realization of such highdimensional GHZ states indicates the potential of our time–frequency platform for quantum technologies such as nearterm quantum computation and clusterstate generation.^{33,45}
Discussion
Hyperentangled time–frequency states, as opposed to other highdimensional optical degrees of freedom like orbital angular momentum, can be generated in integrated onchip sources, which have gained tremendous attention in recent years due to their low cost, room temperature operation, compatibility with CMOS foundries and the ability to be integrated with other optical components. Pulse shapers,^{46} phase modulators,^{47} and switches^{48} can all be demonstrated on a chip, and a series of DWDMs and delay lines can be used to realize the equivalent functionality of onchip CFBG. In addition, demonstration of balanced and unbalanced interferometers onchip eliminates the need for active stabilization, which is of considerable profit for the scalability of the scheme.^{49} These contributions can potentially lead to combining these sources with onchip components designed for manipulation of these states, realizing the whole process on an integrated circuit.
Highdimensional optical states^{25,26,49,50,51} can open the door to deterministically carry out various quantum operations in relatively large Hilbert spaces,^{52} as well as enable higher encoding efficiency in quantum communication protocols, such as quantum key distribution^{22} and quantum teleportation.^{16,53} We have demonstrated deterministic single and twoqudit gates using the time and frequency degrees of freedom of a single photon for encoding—operating on up to 256 (2^{8})dimensional Hilbert spaces—and carried out these gates with a high computationalspace fidelity. We have shown the application of such twoqudit gates in nearterm quantum computation by using them to realize a GHZ state of four parties with 32 dimensions each, corresponding to a Hilbert space of more than one million modes. Such deterministic quantum gates add significant value to the photonic platform for quantum information processing and have direct application in, e.g., simulation of quantum manybody physics.^{54,55,56}
Methods
For the timebin singlequdit X gate shown in Fig. 2, we split the experimental setup in three stages: state preparation, X gate operation, and state measurement. For the state preparation, we use an Agilent 81645A CW laser tuned to 1553.9 nm and send it into an intensity modulator (~4dB insertion loss) and phase modulator (~3dB insertion loss), both manufactured by EOSpace, which are used to create the time bins and control their relative phases, respectively. To implement the X gate operation, we used an MZM with two complementary outputs (~4dB insertion loss), also manufactured by EOSpace. We use a piezobased fiberphase shifter (General Photonics FPS001) to control the phase difference between the two paths following the MZM. Then a 2 × 2 fiber coupler is used to merge the branches. For the state measurement, we used 1bin and 2bin delay interferometers implemented with 2 × 2 fiber couplers and additional piezobased fiber phase shifters. For the timebin X gate and computationalbasis measurements of threedimensional twoqudit gates, gated InGaAs singlephoton detectors (Aurea Technologies SPD_AT_M2) were used. For the rest of the measurements, we used superconducting nanowire singlephoton detectors (Quantum Opus). To measure the arrival times of the photons on the singlephoton detectors, a timeinterval analyzer (PicoQuant HydraHarp 400) is used. The stabilization of the interferometers is done by sending a CW laser at 1550.9 nm in the backward direction and feeding the output power into a computerbased feedback loop to maintain the phase. To stabilize the X gate, we use a similar scheme with an additional circulator at the input of the gate (not shown in the figures) to retrieve the optical power in the backward direction. The signals applied to the intensity modulators and phase modulator, as well as the trigger and synchronization signal of the singlephoton detector and time interval analyzer, are generated by an electronic arbitrary waveform generator (Tektronix AWG7122B) and adjusted to the proper level by linear amplifiers.
To assess the performance of our one and twoqudit quantum gates, we first focus on the computationalbasis fidelity \({\cal{F}}_{\mathrm{C}}\)—one example of a socalled “classical” fidelity in the literature.^{57} Defining n〉 (n = 0, 1, …, N − 1) as the set of all computational basis states and u_{n}〉 as the corresponding output states for a perfect operation, we have the fidelity
where p(u_{n}n) is the probability of measuring the output state u_{n}〉 given an input of n〉. In the operations considered here, the ideal output states u_{n}〉 are members of the computational basis as well, so there is no need to measure temporal or spectral superpositions in determination of \({\cal{F}}_C\). Given the measured counts, we retrieve the N conditional probability distributions via Bayesian mean estimation (BME),^{58,59} where our model assumes that each set of count outcomes (after accidentals subtraction) follows a multinomial distribution with tobedetermined probabilities; for simplicity, we take the prior distributions as uniform (equal weights for all outcomes). We then compute the mean and standard deviation of each value p(u_{n}n) and sum them to arrive at \({\cal{F}}_{\mathrm{C}}\). Specifically, if \(C_{u_nn}\) signifies the counts measured for outcome u_{n}, and C_{totn}, the total counts over all outcomes (both for a given input state n〉), BME predicts:
where the standard deviation in the estimate is used for the error. Since the probabilities here each actually come from N different distributions, we estimate the total error in \({\cal{F}}_C\) by adding these constituent errors in quadrature. Explicitly, we find \({\cal{F}}_{\mathrm{C}} = 0.996 \pm 0.001\) for the X gate, 0.90 ± 0.01 for the CINC operation, 0.92 ± 0.01 for the 3 × 3 SUM gate, and \(\bar {\cal{F}}_{\mathrm{C}} = 0.9589 \pm 0.0005\) for the 16 × 16 SUM gate. The reduction in \({\cal{F}}_{\mathrm{C}}\) for the twoqudit gates is due in large part to the fewer total counts in these cases, from our use of heralded single photons rather than a weak coherent state. As seen by the presence of N in the denominator of Eq. (2), even when \(C_{u_nn} = C_{{\mathrm{tot}}n}\), the estimate p(u_{n}n) is not unity unless C_{totn} ≫ N. In our experiments, the twoqudit tests have only ~100–300 total counts per input computational basis state for the 9 × 9 matrices (with N = 9) and ~500–800 counts per input state for the 16 × 16 matrices (with N = 16), thereby effectively bounding the maximum p(u_{n}n) and, by extension, fidelity \({\cal{F}}_{\mathrm{C}}\). This behavior is actually a strength of BME, though, as it ensures that we have a conservative estimate of the fidelity that is justified by the total amount of data acquired.^{58}
While extremely useful for initial characterization, however, the computationalbasis fidelity above provides no information on phase coherence. On the other hand, process tomography would offer a complete quantification of the quantum gate. Yet due to the challenging experimental complexity involved in quantum process tomography, here we choose a much simpler test which—while limited—nonetheless offers strong evidence for the coherence of our timebin X gate. To begin with, note that all threedimensional quantum processes can be expressed in terms of the nine Weyl operations^{60}:
The quantum process itself is a completely positive map \({\cal{E}}\),^{61} which for a given input density matrix ρ_{in} outputs the state
The process matrix with elements χ_{mn} uniquely describes the operation. The ideal threebin X gate with process matrix χ_{X} has only one nonzero value, [χ_{X}]_{11} = 1. To compare to this ideal, we assume the actual operation consists of a perfect X gate plus depolarizing (white) noise.^{18} In this case we have a total operation modeled as
whose process matrix we take to be \(\chi _N = \lambda \chi _{\mathrm{X}} + \frac{{1  \lambda }}{9}{\Bbb I}_9\), which can be calculated by using \({\Bbb I}_3 = \frac{1}{3}\mathop {\sum}\nolimits_{n = 0}^8 {{\mathrm{U}}_n\rho _{{\mathrm{in}}}{\mathrm{U}}_n^\dagger }\).^{18} if we then assume a pure input superposition state ρ_{in} = ψ_{in}〉〈ψ_{in}, where ψ_{in}〉 ∝ 0〉_{t} + e^{iϕ}1〉_{t} + e^{2iϕ}2〉_{t}, and measure the projection onto the output ψ_{out}〉 ∝ 0〉_{t} + 1〉_{t} + 2〉_{t} (as in Fig. 2c), λ can be estimated from the interference visibility V as^{62}:
and the process fidelity is then given by:
as discussed in the article.
To show the coherence of our SUM gate, we generate an input state in the signal photon which is in timebin 0_{t}〉 and an equiamplitude superposition in frequency \(\left \psi \right\rangle _{{\mathrm{in}}} = \frac{1}{{\sqrt 3 }}\left( {\left 0 \right\rangle _{\mathrm{f}} + \left 1 \right\rangle _{\mathrm{f}} + \left 2 \right\rangle _{\mathrm{f}}} \right)\left 0 \right\rangle _{\mathrm{t}}\). After passing this state through the SUM gate, the timebin state of the photon is shifted based on the frequency, leaving us with a maximally nonseparable state \(\left \psi \right\rangle _{{\mathrm{out}}} = \frac{1}{{\sqrt 3 }}\left( {\left {00} \right\rangle _{{\mathrm{ft}}} + \left {11} \right\rangle _{{\mathrm{ft}}} + \left {22} \right\rangle _{{\mathrm{ft}}}} \right)\). We note that since we are starting with timebin zero, the time bins will not fall out of the computational space; therefore, the interferometric structure is not needed for the SUM gate and a dispersion module alone can do the operation. This saves us the extra insertion loss of the interferometer, which is an important parameter due to the low photon pair rate on the detectors in this particular experiment. To measure the threedimensional nonseparability in ψ_{out}〉, we must vary the phases of different signal frequency bins and time bins with a pulse shaper and phase modulator, respectively. To observe the effect of this phase sweep with our relatively slow singlephoton detectors (with ~100 ps jitter), an indistinguishable projection of all three time bins and frequency bins should be created. In general, the time bins can be projected on an indistinguishable state by using a cascade of interferometers, as illustrated in Fig. 2a. However, in our specific experiment, it is simpler to use a dispersion module with opposite dispersion to that of the module used in the SUM gate to perform the same projection. The frequency bins are then projected on an indistinguishable state using a phase modulator and pulse shaper to mix the frequencies (Fig. 6a)—a technique used previously in ref. ^{26}. We note that our measurements on the signal photons are conditioned on heralding by idler frequency superposition states. To measure the interference between different signal frequency bins, the idler photons too have to be projected on an indistinguishable frequency bin using a phase modulator and pulse shaper (Fig. 6a). This projection guarantees that detection of an idler photon does not give us any information on the frequency of the signal photon. Here the phases of the idler frequency bins are held constant; only the phases of the signal frequency and time bins are varied. This is in contrast to experiments in ref. ^{26}, where the phases of both signal and idler frequency bins were varied.
In our experiment, threedimensional interference measurements were not possible since mixing all three frequencies together adds extra projection loss, which we cannot afford. Therefore, we vary the phases of different time bins and frequency bins to measure twodimensional interference patterns between all three time bins and frequency bins (Fig. 6c). Using the visibilities of these interference patterns along with a joint spectral intensity (JSI) measurement (Fig. 6b) can give us a lower bound on the amount of nonseparability present in our system by measuring entanglement of formation.^{35,36} The JSI denotes the correlations between the time bins and frequency bins of a signal photon heralded by an idler photon in its computational basis. This measurement was done using the same experimental setup used in Fig. 6a, without the equipment used for sweeping the phase of different signal time bins and projection measurements. For this measurement, the idler photons were detected after PS1, and the signal photons were detected right after the SUM gate.
Having the JSI measurement and the twodimensional interference visibilities in hand, we have all the data needed to calculate the entanglement of formation with the assumption of having only white noise in our system, which can be expressed as:
where
Here, C is the number of indices (j, k) used in the sum. This measurement is useful when we do not have access to all the elements of the density matrix. 〈j, jρk, k〉 (j ≠ k) elements indicate the coherence between modes j and k, and can be lowerbounded using the twodimensional visibilities. The terms 〈j, kρj, k〉 can be calculated using the elements of the JSI. Using these values, we measure E_{of} ≥ 1.19 ± 0.12 ebits, which indicates greater than twodimensional nonseparability in our twoparty system, more than one standard deviation away from the threshold.
To generate the 32dimensional fourparty GHZ state, the signal and idler go through the same dispersion module (−2 ns.nm^{−1}). After dispersion, the signal frequency bins farther away from the center of the spectrum are delayed more, but the idler frequency bins are delayed less as we move farther away from the center. In order to write the GHZ state in the form \(\left \psi \right\rangle _{{\mathrm{out}}} = \frac{1}{{\sqrt {32} }}\mathop {\sum}\nolimits_{m = 0}^{31} {\left {m,m,m,m} \right\rangle _{{\mathrm{f}}_{\mathrm{s}}{\mathrm{t}}_{\mathrm{s}}{\mathrm{f}}_{\mathrm{i}}{\mathrm{t}}_{\mathrm{i}}}}\), we label the signal time bins after dispersion 0 to 31 starting from earlier time bins (time bin 0 the earliest, time bin 31 the latest), while on the idlers, we label the time bins such that the earliest time bin is 31 and the latest time bin is 0. Another choice would be to send signal and idler through separate modules with equal but opposite dispersion, in which case we would use identical time labeling. To measure the state illustrated in Fig. 5, we individually measured coincidences for the 32 different settings of both signal and idler frequency bins (32 × 32 measurements). For each of these measurements, we used our event timer to assign signal and idler time bins for each coincidence, which results in a 32 × 32 submatrix for each signalidler frequency setting. Therefore, we have 32^{4} measurements in total. Two of the 32 × 32 timebin submatrices are shown in Fig. 5b, c.
We use bulk switches, dispersion modules, pulse shapers, and phase modulators in out experiments, which have high insertion loss (switch: 3 dB, dispersion module: 3 dB, pulse shaper: 5 dB, phase modulator: 3 dB). Therefore, we use very bright entangled photons at the input in order to have reasonable coincidence counts on our detectors in our acquisition time. Using bright biphotons gives rise to multipair generation which leads to the relatively high accidental rate here.
Data availability
The data and analysis codes used in this study are available from the corresponding author on request.
References
 1.
Shor, P. W. Polynomialtime algorithms for prime factorization and discrete logarithms on a quantum computer. SIAM Review 41, 303–332 (1999).
 2.
Kimble, H. J. The quantum internet. Nature 453, 1023–1030 (2008).
 3.
Knill, E., Laflamme, R. & Milburn, G. J. A scheme for efficient quantum computation with linear optics. Nature 409, 46–52 (2001).
 4.
Crespi, A. et al. Integrated photonic quantum gates for polarization qubits. Nat. Commun. 2, 566 (2011).
 5.
Babazadeh, A. et al. Highdimensional singlephoton quantum gates: concepts and experiments. Phys. Rev. Lett. 119, 180510 (2017).
 6.
Humphreys, P. C. et al. Linear optical quantum computing in a single spatial mode. Phys. Rev. Lett. 111, 150501 (2013).
 7.
Lu, H.H. et al. Electrooptic frequency beam splitters and tritters for highfidelity photonic quantum information processing. Phys. Rev. Lett. 120, 030502 (2018).
 8.
Lu, H.H. et al. A controlledNOT gate for frequencybin qubits. npj Quantum Information 5, 24 (2019).
 9.
Fiorentino, M. & Wong, F. Deterministic controlledNOT gate for singlephoton twoqubit quantum logic. Phys. Rev. Lett. 93, 070502 (2004).
 10.
Kagalwala, K. H., Di Giuseppe, G., Abouraddy, A. F. & Saleh, B. E. A. Singlephoton threequbit quantum logic using spatial light modulators. Nat. Commun. 8, 739 (2017).
 11.
Fang, W.T. et al. Onchip generation of timeand wavelengthdivision multiplexed multiple timebin entanglement. Opt. Express 26, 12912–12921 (2018).
 12.
Humphreys, P. C. et al. Continuousvariable quantum computing in optical timefrequency modes using quantum memories. Phys. Rev. Lett. 113, 130502 (2014).
 13.
MarinPalomo, P. et al. Microresonatorbased solitons for massively parallel coherent optical communications. Nature 546, 274–279 (2017).
 14.
Cerf, N. J., Adami, C. & Kwiat, P. G. Optical simulation of quantum logic. Phys. Rev. A 57, R1477 (1998).
 15.
Roa, L., Delgado, A. & FuentesGuridi, I. Optimal conclusive teleportation of quantum states. Phys. Rev. A 68, 022310 (2003).
 16.
Boschi, D., Branca, S., De Martini, F., Hardy, L. & Popescu, S. Experimental realization of teleporting an unknown pure quantum state via dual classical and einsteinpodolskyrosen channels. Phys. Rev. Lett. 80, 1121 (1998).
 17.
Sheng, Y.B., Deng, F.G. & Long, G. L. Complete hyperentangledBellstate analysis for quantum communication. Phys. Rev. A 82, 032318 (2010).
 18.
Wilde, M. M. Quantum information theory. (Cambridge University Press, Cambridge, UK, 2013).
 19.
Wootters, W. K. & Fields, B. D. Optimal statedetermination by mutually unbiased measurements. Ann. Phys. 191, 363–381 (1989).
 20.
Lukens, J. M., Islam, N. T., Lim, C. C. W. & Gauthier, D. J. Reconfigurable generation and measurement of mutually unbiased bases for timebin qudits. Appl. Phys. Lett. 112, 111102 (2018).
 21.
Sheridan, L. & Scarani, V. Security proof for quantum key distribution using qudit systems. Phys. Rev. A 82, 030301 (2010).
 22.
Islam, N. T., Lim, C. C. W., Cahall, C., Kim, J. & Gauthier, D. J. Provably secure and highrate quantum key distribution with timebin qudits. Sci. Adv. 3, e1701491 (2017).
 23.
JaramilloVillegas, J. A. et al. Persistent energytime entanglement covering multiple resonances of an onchip biphoton frequency comb. Optica 4, 655–658 (2017).
 24.
Reimer, C. et al. Generation of multiphoton entangled quantum states by means of integrated frequency combs. Science 351, 1176–1180 (2016).
 25.
Kues, M. et al. Onchip generation of highdimensional entangled quantum states and their coherent control. Nature 546, 622–626 (2017).
 26.
Imany, P. et al. 50GHzspaced comb of highdimensional frequencybin entangled photons from an onchip silicon nitride microresonator. Opt. Express 26, 1825–1840 (2018).
 27.
Brennen, G. K., Bullock, S. S. & O’Leary, D. P. Efficient circuits for exactuniversal computations with qudits. Quantum Info Comput. 6, 436–454 (2005).
 28.
Wang, X., Sanders, B. C. & Berry, D. W. Entangling power and operator entanglement in qudit systems. Phys. Rev. A 67, 042323 (2003).
 29.
Draper, T. G., Kutin, S. A., Rains, E. M. & Svore, K. M. A logarithmicdepth quantum carrylookahead adder. Quantum Info Comput. 6, 351–369 (2006).
 30.
Imany, P., Odele, O. D., JaramilloVillegas, J. A., Leaird, D. E. & Weiner, A. M. Characterization of coherent quantum frequency combs using electrooptic phase modulation. Phys. Rev. A 97, 013813 (2018).
 31.
Li, L., Yi, X., Huang, T. X. H. & Minasian, R. A. Microwave photonic filter based on dispersion controlled spectrum slicing technique. Electron. Lett. 47, 511 (2011).
 32.
Li, L., Yi, X., Huang, T. X. H. & Minasian, R. A. Distortionfree spectrum sliced microwave photonic signal processor: analysis, design and implementation. Opt. Express 20, 11517–11528 (2012).
 33.
Reimer, C. et al. Highdimensional oneway quantum processing implemented on dlevel cluster states. Nat. Phys. 15, 148 (2019).
 34.
Karimi, E. & Boyd, R. W. Classical entanglement? Science 350, 1172–1173 (2015).
 35.
Martin, A. et al. Quantifying Photonic HighDimensional Entanglement. Phys. Rev. Lett. 118, 110501 (2017).
 36.
Tiranov, A. et al. Quantification of multidimensional entanglement stored in a crystal. Phys. Rev. A 96, 040303 (2017).
 37.
Hu, X.M. et al. Experimental multilevel quantum teleportation. arXiv preprint arXiv:1904.12249 (2019).
 38.
Luo, Y.H. et al. Quantum teleportation in high dimensions. arXiv preprint arXiv:1906.09697 (2019).
 39.
Barreiro, J. T., Wei, T. C. & Kwiat, P. G. Beating the channel capacity limit for linear photonic superdense coding. Nat. Phys. 4, 282 (2008).
 40.
Erhard, M., Malik, M., Krenn, M. & Zeilinger, A. Experimental Greenberger–Horne–Zeilinger entanglement beyond qubits. Nat. Photonics 12, 759–764 (2018).
 41.
Pan, J. W., Bouwmeester, D., Daniell, M., Weinfurter, H. & Zellinger, A. Experimental test of quantum nonlocality in threephoton GreenbergerHorneZeilinger entanglement. Nature 403, 515–519 (2000).
 42.
Wang, X.L. et al. 18Qubit entanglement with six photons’ three degrees of freedom. Phys. Rev. Lett. 120, 260502 (2018).
 43.
Hillery, M., Bužek, V. & Berthiaume, A. Quantum secret sharing. Phys. Rev. A 59, 1829 (1999).
 44.
Zhao, Z. et al. Experimental demonstration of fivephoton entanglement and opendestination teleportation. Nature 430, 54–58 (2004).
 45.
Pant, M., Towsley, D., Englund, D. & Guha, S. Percolation thresholds for photonic quantum computing. Nat. Commun. 10, 1070 (2019).
 46.
Wang, J. et al. Reconfigurable radiofrequency arbitrary waveforms synthesized in a silicon photonic chip. Nat. Commun. 6, 5957 (2015).
 47.
Wang, C. et al. Integrated lithium niobate electrooptic modulators operating at CMOScompatible voltages. Nature 562, 101–104 (2018).
 48.
Wang, C. et al. Monolithic lithium niobate photonic circuits for Kerr frequency comb generation and modulation. Nat. Commun. 10, 978 (2019).
 49.
Wang, J. et al. Multidimensional quantum entanglement with largescale integrated optics. Science 360, 285–291 (2018).
 50.
Ikuta, T. & Takesue, H. Fourdimensional entanglement distribution over 100 km. Sci. Rep. 8, 817 (2018).
 51.
Malik, M. et al. Multiphoton entanglement in high dimensions. Nat. Photonics 10, 248–252 (2016).
 52.
Erhard, M., Fickler, R., Krenn, M. & Zeilinger, A. Twisted photons: new quantum perspectives in high dimensions. Light Sci. Appl. 7, 17146 (2018).
 53.
Bouwmeester, D. et al. Experimental quantum teleportation. Nature 390, 575–579 (1997).
 54.
Lu, H.H. et al. Simulations of subatomic manybody physics on a quantum frequency processor. arXiv Preprint arXiv:1810.03959 (2018).
 55.
Sparrow, C. et al. Simulating the vibrational quantum dynamics of molecules using photonics. Nature 557, 660–667 (2018).
 56.
Bernien, H. et al. Probing manybody dynamics on a 51atom quantum simulator. Nature 551, 579–584 (2017).
 57.
De Greve, K. et al. Complete tomography of a highfidelity solidstate entangled spinphoton qubit pair. Nat. Commun. 4, 2228 (2013).
 58.
BlumeKohout, R. Optimal, reliable estimation of quantum states. New J. Phys. 12, 043034 (2010).
 59.
Williams, B. P. & Lougovski, P. Quantum state estimation when qubits are lost: a nodataleftbehind approach. New J. Phys. 19, 043003 (2017).
 60.
Bertlmann, R. A. & Krammer, P. Bloch vectors for qudits. J. Phys. A Math. Theor. 41, 235303 (2008).
 61.
O’Brien, J. L. et al. Quantum process tomography of a controlledNOT gate. Phys. Rev. Lett. 93, 80502 (2004).
 62.
Thew, R. T., Acín, A., Zbinden, H. & Gisin, N. Belltype test of energytime entangled qutrits. Phys. Rev. Lett. 93, 010503 (2004).
Acknowledgements
This work was funded in part by the National Science Foundation under award number 1839191ECCS. We thank B.P. Williams for discussions on Bayesian estimation. We thank M. Hosseini, N. Lingaraju, and Nathan O’Malley for discussions. J.M.L. acknowledges support from a Wigner Fellowship at ORNL. A portion of this work was performed at Oak Ridge National Laboratory, operated by UTBattelle for the U.S. Department of Energy under Contract No. DEAC0500OR22725.
Author information
Affiliations
Contributions
P.I. initiated and developed the idea. P.I., J.A.J. and M.S.A. performed the experiments with help from D.E.L. P.I., J.A.J., J.M.L., M.S.A., O.D.O. and A.M.W. performed the theoretical work and analyzed the results. P.I. and A.J.M. analyzed the data to quantify the nonseparablity in the system. M.Q. supervised the fabrication of the microring resonator. P.I., J.A.J., J.M.L., M.S.A., O.D.O. and A.M.W. prepared the manuscript. A.M.W. managed and supervised the project.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Imany, P., JaramilloVillegas, J.A., Alshaykh, M.S. et al. Highdimensional optical quantum logic in large operational spaces. npj Quantum Inf 5, 59 (2019). https://doi.org/10.1038/s4153401901738
Received:
Accepted:
Published:
Further reading

648 Hilbertspace dimensionality in a biphoton frequency comb: entanglement of formation and Schmidt mode decomposition
npj Quantum Information (2021)

Separability of heterogeneous quantum systems using multipartite concurrence and tangle
Quantum Information Processing (2021)

Nonbinary entanglementassisted stabilizer codes
Quantum Information Processing (2021)

Separability of heterogeneous multipartite quantum systems using Bloch Vectors
Quantum Information Processing (2020)