Path-polarization hyperentangled and cluster states of photons on a chip.

Encoding many qubits in different degrees of freedom (DOFs) of single photons is one of the routes toward enlarging the Hilbert space spanned by a photonic quantum state. Hyperentangled photon states (that is, states showing entanglement in multiple DOFs) have demonstrated significant implications for both fundamental physics tests and quantum communication and computation. Increasing the number of qubits of photonic experiments requires miniaturization and integration of the basic elements, and functions to guarantee the setup stability, which motivates the development of technologies allowing the precise control of different photonic DOFs on a chip. We demonstrate the contextual use of path and polarization qubits propagating within an integrated quantum circuit. We tested the properties of four-qubit linear cluster states built on both DOFs, and we exploited them to perform the Grover's search algorithm according to the one-way quantum computation model. Our results pave the way toward the full integration on a chip of hybrid multi-qubit multiphoton states.

Introduction

Novel integrated photonic circuits built on a single chip have been recently introduced within the realm of quantum information[1], disclosing new perspectives toward quantum communication[2], quantum computation[3] and the quantum simulation of physical phenomena[4, 5, 6, 7, 8, 9, 10]. The miniaturization of integrated photonic devices represents a necessary step toward the implementation of state-of-the-art quantum information protocols, which require cascading a high number of elements with an excellent stability that is impossible to achieve with standard bulk optical setup.

The next generation of integrated quantum circuits (IQCs), which incorporates highly efficient photon sources[11, 12, 13, 14] and detectors[15, 16, 17, 18, 19], is expected to have a large impact on future photonic quantum technologies and are essential to achieve a level of complexity and stability of the operations higher than what has been previously demonstrated. Nevertheless, the generation of a large number of photons, in particular on-chip, is still a very challenging task; therefore, it is appealing to encode more than one qubit on the same photon by exploiting different degrees of freedom (DOFs).

Inherently stable interferometer networks, composed of waveguides, beam splitters and phase shifters, built in two dimensions on different material platforms, such as silicon, silicon nitride and others, are achieved by lithography, a well-established technique already developed for telecom wavelengths. This approach allows the fabrication of a large number of replicas of the same circuit using a single mask and represents the strongest candidate for a large-scale production of IQCs. Although it has been demonstrated that the operation complexity performed by such systems may be very high, only path-encoded qubits are allowed in such systems because polarization qubits are degraded by the intrinsic large birefringence of the material substrate or of the waveguide itself. On the other hand, several applications in the quantum domain, such as quantum computation and quantum communications, may greatly benefit from the possibility of manipulating and controlling polarization qubits.

Femtosecond laser writing, recently introduced for IQC applications[20, 21, 22], allows writing circular transverse waveguide profiles in three dimensions that are able to support the propagation of nearly Gaussian modes with any polarization state while simultaneously keeping the phase of path-encoded qubits highly stable. Moreover, this technique makes it possible to perform arbitrary transformations of the polarization state by suitable integrated devices, such as polarization beam splitters[23] and waveguide-based optical waveplates[22, 24].

Currently, the building blocks necessary to perform the basic operations with path and polarization qubits are available. The time is right to demonstrate the simultaneous control of two different DOFs of the photons within the same chip. In this work, we manipulate path- and polarization-encoded qubits belonging to a two-photon four-qubit hyperentangled (HE)/cluster state[25, 26, 27, 28, 29] and propagating through an IQC. This enables us to demonstrate the Grover's algorithm on a four-element database in a one-way quantum-computing approach.

Materials and methods

Our experimental setup is depicted in Figures 1 and 2, and consists of a HE-photon source, a manipulation stage, which includes the integrated photonic chip, and a detection stage.

Figure 1

Experimental setup. The state is generated by the hyperentangled source in the blue area. The four-hole screen in (a) selects the spatial modes , rA and rB. In the green area, the four modes are independently addressed using half-mirrors, half-waveplates and half-lenses to compensate their polarization and to couple them into the fiber array. A translation stage is used to control the path indistinguishability between left and right modes. A glass plate in mode is used to control the path phase. The chip is connected to the fiber array using a 6-axis stage (NanoMax, Thorlabs Inc., Newton, NJ, USA). In the orange area, measurements are performed: two of the outputs of the chip are addressed into multi-mode fibers connected to SPADs in coincidence mode. Polarization analysis is performed using QWP, HWP and PBS. (b) Schematic representation of the integrated device. (c) Hong-Ou-Mandel (HOM) dip of BSA and BSB. BBO, beta barium borate; BS, beam splitter; HWP, half waveplate; MMF, multi-mode fiber; PBS, polarising beam splitter; QWP, quarter waveplate; SMF, single-mode fiber; SPAD, single-photon avalanche diode.

Figure 2

Chip mount assembly. GRIN, gradient index.

The source generates pairs of photons that are HE in the path and polarization DOFs via spontaneous-parametric down-conversion at 710 nm by a BBO type I nonlinear crystal[25]. Polarization entanglement is produced by optical superposition of two cones of perpendicular polarization created by double passage of a ultraviolet laser pump through the BBO crystal from both sides[30]. Path entanglement is generated by selecting two pairs of correlated spatial modes with a four-hole screen, namely , rB and rA, .

The generated HE state is

where A and B are the two photons; rA, , rB and identify the four spatial modes; and H and V represent the two possible polarization states for each photon. The phases θ and φ, corresponding, respectively, to the polarization and path DOFs, can be independently controlled.

The chip consists of balanced waveguide beam splitters (BSA and BSB), yielding polarization insensitive behavior[31] and fabricated by femtosecond laser waveguide writing[32] using the second harmonic (λ=515 nm) of a Yb:KYW cavity-dumped laser oscillator (300 fs pulse duration, 1 MHz repetition rate). Femtosecond laser pulses are focused by a 0.6 numerical aperture microscope objective into the volume of the glass substrate (EAGLE 2000, Corning Inc., Corning, NY, USA), where nonlinear energy absorption creates a permanent and localized refractive index increase. Waveguides are produced by smoothly translating the sample under the laser beam, using 3-axis air-bearing stages (FiberGLIDE 3D, Aerotech Inc., Pittsburgh, PA, USA). Under proper irradiation conditions (100 nJ pulse energy and 10 mm s−1 translation speed), single-mode waveguides at 710 nm are produced at 170 μm depth below the glass surface and are characterized by a mode diameter of ~8 μm, propagation loss of 1.5 dB cm−1 and coupling loss to single-mode fibers <1 dB per facet. Integrated beam splitters are realized following a particular three-dimensional directional-coupler design[30] that ensures insensitivity to polarization. To achieve a balanced splitting ratio, the waveguides, initially 127 μm apart, are brought closer (with a bending radius of 90 mm) to 10 μm for an interaction length of 1.8 mm. The overall chip length is 34 mm.

The four modes rA, , rB and are coupled to the input ports of BSA and BSB through an 8-cm-long, single-mode fiber array, terminated at the input side by a set of four Gradient Index lenses. One of the main technical issues is the independent injection of each mode into the corresponding Gradient Index lens, which is performed using a set of custom half-mirrors and one half-lens for each mode. Moreover, polarization compensation is individually performed for each mode through a properly chosen set of half-waveplates (HWPs) and quarter waveplates. Finally, the fiber array is thermally insulated from the environment to guarantee path stability.

Output light from the chip is collected by a 10 × objective. The coupling ratios of the four modes in each optical component and the overall transmission efficiency are provided in the Supplementary Information.

Two interference filters centered at 710 nm select 10-nm bandwidth and ensure photon indistinguishability, and two avalanche photodiode detectors measure coincidences between output modes and , which are coupled to the detectors through multi-mode fibers (we label , and , as the output modes of BSA and BSB, respectively).

Figure 1c shows the Hong-Ou-Mandel dips obtained when the two photons are injected within BSA and BSB. We obtained the following visibilities: VA=(0.976±0.003) for BSA and VB=(0.982±0.003) for BSB, confirming the correct operation of the two systems. After removing Nacc=12 accidental coincidences every 30 s, we obtain VA–net=(0.985±0.003) and VB–net=(0.991±0.003).

Results and discussion

In the first experiment, we injected the path-polarization HE state into the chip:

where and . This state was obtained from Equation (1) by introducing the polarization transformation using the compensation waveplates shown in Figure 1 such that H→H and V→V on modes and rA and H→V and V→H on and rB. In this way, we were able to guarantee the polarization compensation over the entire system, which includes the fiber array and the chip, on the computational polarization basis. To compensate the two pairs of correlated modes on the diagonal basis, we introduced an additional half-HWP at zero degrees before the compensation plates on mode . Because the addition of the plate preserves the compensation in the computational basis, we tilted it along its vertical axis to keep the phase difference between the two couples of modes. This ensures the simultaneous compensation on the polarization DOF constant. Arbitrary values for the parameters θ and φ in Equation (1) can be set by translating the spherical mirror within the HE source and by tilting an additional glass plate on mode , respectively. The plus signs appearing in Equation (2) were achieved by setting φ, θ=0, and the minus signs were obtained by setting φ, θ=π. The overall symmetry of the state determines the behavior of the two photons: if the wavefunction of the HE state is symmetric, they emerge from the same output port of the beam splitters, corresponding to a coincidence dip, whereas the expected result in the case of an antisymmetric wavefunction is a coincidence peak. This behavior can be analyzed by recovering the dips and peaks of path entanglement by varying both φ and θ so that a dip can be obtained with (φ, θ)=(0, 0), (π, π) and a peak with (φ, θ)=(0, π), (π, 0). The results are shown in Figure 3. The average peak/dip visibilities are Vpeak=(0.93±0.20) and Vdip=(0.860±0.005). These results are comparable with those of Ref. 25 and prove the achievement of path-polarization hyperentanglement on a chip with good fidelity.

Figure 3

HE state interference. Interference pattern between modes and when varying both θ and φ. Peaks and dips derive from the symmetry of the entire HE wave function. Experimental data of each curve are fitted with the convolution of a Gaussian with a sinc function; error bars are derived from the Poisson statistical distribution of counts. Errors on the peak/dip visibilities are estimated using a Monte Carlo statistical analysis of the experimental data.

In a second experiment, the same hyperentanglement source was used to engineer a four-qubit cluster state encoded in the path and polarization DOFs of the two photons:

where . At variance with the standard HE states, four-qubit cluster states are not biseparable and present genuine multipartite entanglement[33]. The states can be thought of as graphs in which the vertices are the physical qubits initially in the state , and each edge represents a controlled-phase gate entangling the two connected nodes. As an example, in this graphical representation, an HE state is equivalent, up to single-qubit transformations, to the graph formed by two disjointed couples of vertices, as shown in Figure 4a. A cluster state can then be obtained by connecting these pairs of qubits. In our case, this was easily performed by inserting a zero-order HWP in mode rA oriented along the optical axis, as explained later. In the one-way quantum computation model[34], cluster states provide the entire resource for the computation because the information is written, processed and read out by single-qubit measurements on the physical qubits of the cluster. This approach to quantum computation transfers the main complexity of the process from the ability to implement multi-qubit gates to the capability to create the initial cluster state. It is therefore well-suited for quantum optical schemes where states can be produced with high fidelity and photon–photon interactions are difficult to achieve, thus precluding two-qubit gates from linear optical circuits.

Figure 4

Cluster state. (a) Graphical representation of four-qubit hyperentangled and cluster states. (b) Real part of the two-qubit polarization tomography of pair of the cluster state, which is compensated to be a state. (c) Real part of the two-qubit polarization tomography of pair of the cluster state, compensated to be a state. The imaginary parts of the two tomographies are negligible.

Starting from the HE state and inserting a zero-order, zero-degree HWP in mode rA, we were able to change the polarization of from to while keeping unchanged, thus creating the cluster state . The quantum-state tomographies corresponding to the two mode pairs are reported in Figure 4b and 4c and correspond to the following parameters:

Here, F is the fidelity and C is the concurrence of the experimental state, and errors are calculated from a Monte Carlo analysis of the experimental data.

The obtained results are comparable with those of the tomographies of the hyperentanglement source reported in the Supplementary Information, proving that the cluster state is correctly generated.

The stabilizer formalism, explained in Ref. 35 can be adopted to measure a genuine multipartite entanglement witness:

where upper case X and Z define the Pauli operators for the polarization of the state, and lower case x and z define the Pauli operators for the momentum. The state is entangled when −1≤W<0 and is purely entangled for . The polarization-dependent stabilizers are measured by rotating the analysis waveplates (Figure 1). The two beam splitters perform the transformation in the path of each photon and , where is achieved by changing the phase φ, that is, tilting the corresponding glass plate on mode r. This manipulation allows us to measure the momentum-dependent stabilizers involving x.

We report in Table 1 the measured outcomes for the stabilizers. The overall value of demonstrates that the state presents genuine multipartite entanglement, and we can derive a lower bound[35] for the fidelity of the created cluster state This result is comparable with that of Ref. 28.

Table 1

Measured outcomes of the stabilizer used for calculating the genuine multipartite entanglement witness W

ZAZB	+0.940±0.028
XAXBzA	−0.860±0.030
XAXBzB	+0.860±0.030
zAzB	−0.990±0.007
ZAxAxB	+0.809±0.036
ZBxAxB	+0.808±0.035

Finally, the quality of the created cluster state was tested by performing the Grover's search algorithm for a four-element database. A linear optics implementation of the Grover's algorithm using multiple DOFs was originally proposed by Ref. 36. Here we present, to the best of our knowledge, the first achievement of a one-way quantum computation basic operation based on multiple DOFs in integrated photonics. The quantum circuit is represented in Figure 5a. It consists of two qubits initially prepared in the state , a black box tagging one item of the database and an operation that allows us to identify the tagged item in the readout. The algorithm can be implemented[34, 37] using the four-qubit box cluster state defined in Refs. 38 and 39. Following Figure 5a, the black box tags the item by choosing the bases α and β for the measurements on qubits 1 and 4. The information is processed and read on qubits 2 and 3. Labeling the physical qubits in the order (1,2,3,4)=(kB, πA, kA, πB), with k (π) standing for the path (polarization) qubits, the cluster state is equivalent to up to the single-qubit unitaries

Figure 5

Grover's algorithm scheme. (a) Quantum circuit for Grover's algorithm. The black box tags the item through the measurements of qubits 1 and 4. The information is then encoded and processed in qubits 2 and 3 of the cluster state . The single-qubit operations are implemented by the choice of the measurement basis. (b) Outcome probability for different tagged items for the probabilistic Grover's algorithm. The average success rate of the algorithm is s=(0.960±0.007).

can be implemented by simply rotating the measurement basis because is a single-qubit transformation.

First, we performed a probabilistic computation where we post-selected the cases with no errors occurring in the one-way computation model[39, 40]. The results are reported in Figure 5b, in which we show that the average success rate in identifying the correct item in the database is s=(0.960±0.007) at an average protocol rate of 17 Hz. This result is probabilistic and depends on the post selection of the measurement outcomes. We then applied a feed-forward protocol in which the outcomes are relabeled depending on the results of the measurements performed by the black box[38, 39]. In this case, the computation is deterministic with a success rate of s=0.964±0.003 and a protocol rate of 68 Hz. This procedure of passive feed-forward corresponds to corrections made in the post-selection process by relabeling the outputs.

In this work, four-qubit HE states built on the path and the polarization of two photons have been adopted to manipulate qubits based on the two DOFs. Furthermore, they propagated through an integrated photonic circuit fabricated by the femtosecond laser writing technique.

The device was used to test the presence of path and polarization entanglement. We also engineered cluster states, and we measured a multipartite genuine entanglement witness to estimate the purity of correlations within the entangled state. For both experiments, we achieved fidelities comparable with previous bulk experiments[25, 28], thus showing the good quality of the experimental results. We exploited the cluster state to perform Grover's search algorithm as an experimental realization on a chip of one-way quantum computation using different DOFs. Our experiment provides the first demonstration that proves it is possible to simultaneously use different DOFs of photons within an integrated photonic circuit, thus increasing the number of qubits. A future step forward in this approach will include the use of more complex circuits for the active manipulation of both phase and polarization inside the same device.

Notes

During the review process of this work, an implementation of measurement-based, on-chip quantum operations has been reported[41].