Electro-optic analyzer of angular momentum hyperentanglement.

Characterizing a high-dimensional entanglement is fundamental in quantum information applications. Here, we propose a theoretical scheme to analyze and characterize the angular momentum hyperentanglement that two photons are entangled simultaneously in spin and orbital angular momentum. Based on the electro-optic sampling with a proposed hyper-entanglement analyzer and the simple matrix operation using Cramer rule, our simulations show that it is possible to retrieve effectively both the information about the degree of polarization entanglement and the spiral spectrum of high-dimensional orbital angular momentum entanglement.

A twisted light beam can carry spin and orbital angular momentum (OAM), even down to the single photon limit12. Photon spin is associated with circular polarizations, while OAM arises from helical phase structures described by , where is integer and is the azimuthal angle3. Recent years have witnessed a growing interest in angular momentum hyperentanglement, where two photons are entangled simultaneously in spin and OAM degrees of freedom4. Hyperentanglement has been proved an effective tool to enlarge the Hilbert space5. Benefitting from hyperentanglement, the violation of Einstein-Podolsky-Rosen (EPR) local realism was enhanced6, and the possibility to beat the channel capacity limit for linear photonic superdense coding was shown7. Besides, increase of Shannon dimensionality8, remote preparation of single-photon hybrid state9, and super-resolving phase measurement10 were reported.

As is well known, characterizing the OAM spectrum (or spiral spectrum) of entangled photons generated by spontaneous parametric down-conversion (SPDC) is crucial, since most quantum applications are based on the availability of specific quantum states. In theory, a full characterization of quantum spiral bandwidth in terms of Laguerre-Gaussian modes has been given analytically11. But the experimental implementation is not a trivial task. Based on half-integer spiral phase plate or multisector phase masks, the Shannon dimensionality that characterize the quantum channel capacity was measured1213. As an overall quantifier of entanglement, Shannon measurement does not touch the full knowledge of spiral spectrum. By combining an interferometer with image rotator, a clever experimental scheme was reported for the first time to measure the complete spiral spectrum14. Besides, compressive sensing was employed as an effective technique to reconstruct the high-dimensional OAM entanglement15. Spatial light modulators are also used to measure the full spiral bandwidth, and angular version of strong EPR correlation16 as well as violation of generalized Bell inequalities17 was reported. Here, we report theoretically another feasible scheme to retrieve both the information about the degree of polarization entanglement and the high-dimensional spiral spectrum involved in an angular momentum hyperentanglement. Our approach is mainly based on the electro-optic sampling with a proposed hyperentanglement analyzer and the simple matrix operations using Cramer’s rule.

Results

The proposed setup includes two stages, namely, hyperentanglement preparation and analysis, shown in Fig. 1. An ultraviolet (UV) beam (frequency ) is incident on a thin nonlinear crystal cut for collinear type-II phase-matching, where degenerate photon pairs (frequency ) with orthogonal polarizations are produced via SPDC18. A dispersion prism is used to separate down-converted photons from the pump ones. Then pairs of down-converted photons are further separated probabilistically by a non-polarizing beam splitter (50:50 BS) with 50–50% reflection and transmission coefficients. Therefore post-selecting the cases where photons exit in different ports by means of a coincidence measurement (&) yields the entanglement of polarization, . To show the ability of our scheme in characterizing polarization entanglement, we prepare more general states with variable degree of entanglement. We employ a partial polarizer, consisting of a series of coated glass slabs tilted in approximately their Brewster’s angle, to make such nonmaximally entangled states19. Assume that the transmission probability for horizontal and vertical polarizations of the partial polarizer are denoted by T and T, respectively. If we place such a partial polarizer in one path, then we are able to transfer into nonmaximally entangled states after polarization filtering, , where is the tunable degree of entanglement. Meanwhile, quantum entanglement is also established well in the spatial degree of freedom. In general, if we assume that the UV pump beam is Gaussian, then the wave function in terms of OAM can be written in a more general form, , where describes the probability amplitude of finding one photon in OAM mode and the other in 20. An additional mirror (M) compensates the reflection occurring in BS, such that the OAM signs are not flipped in optical path 1. In this scenario, two photons have been made in a hyperentangled state,

Figure 1

Our proposed experimental setup for (a) preparation and (b) analysis of two-photon spin-orbit hyperentanglement.

In the analysis stage, the key elements are our proposed hyperentanglement analyzers. Let us first briefly introduce its function of spin-dependent and electrically controllable OAM generation, see ref 21 for more details. Such an analyzer consists of two nominally conjugated spiral phase plate (SPP) made of z-cut electro-optic ZnTe crystal, and both incident and exist faces are coated with transparent electrodes for applying an electric field to induce the Pockels effect. As the crystalline coordinates of two SPP are designed to have a relative orientation of 90°, so the horizontal and vertical polarization after traversing the analyzer will acquire optical vortices of opposite helix. And the topological charge Q is linearly dependent on the applied voltage U, namely, Q = αU, where (n0 and γ63 are refractive index and electro-optic coefficient at wavelength λ). At the single photon level, the function of our electro-optic (EO) analyzer can be described by an operator,

where Q can take in principle any integer or non-integer values by adjusting a suitable voltage conveniently. Of particular interest is that fractional Q can be decomposed into a coherent superposition of theoretically infinite OAM eigenstates , namely, , where . Recently, we used such an EO analyzer to demonstrate a quantum protocol of hybrid teleportation22. Here we further exploit its possibility for engineering and characterizing high-dimensional angular momentum hyperentanglement. To this end, we incorporate two EO analyzer biased with the same voltage U, each in one path. After interaction with two EO analyzers as described by equation (2), the hyperentangled state of equation (1) becomes,

As mentioned above, the horizontal and vertical polarizations, H and V, yield opposite optical vortices, +Q and −Q, respectively. Reasonably, we can infer that the information of polarization entanglement can be transferred and contained in the OAM entanglement. Thus the labels of H and V in equation (3) can be discarded, and two-photon state is simplified to,

In contrast, if the vortex generation is not spin-dependent (e.g. using traditional spiral phase plates2324), namely, and , then equation (4) only becomes . Obviously, the information about the polarization entanglement has been lost.

The emerging photons from our EO analyzers are collected by a lens and focused into the single-mode fiber (SMF) which supports exclusively the fundamental Gaussian mode with zero OAM. Inversely, SMF in connection with single-photon detector (D) projects the incoming photons into the OAM state of. The output of both single-photon detectors are fed to a coincidence counting circuit (&), which records a joint detection probability,

Considering 14, equation (5) is naturally simplified into,

which denotes the joint probability amplitude.

The equation (6) forms the key result in present work. One of its characteristic features is the nontrivial role played by fractional Q. If Q is integer, and , then only those OAM modes of can be detected, as can be inferred from equation (6). In contrast, if Q is fractional, and , then each mode will contribute towards the coincidence counts via . The ability of fractional Q to access many OAM modes can be understood from its multidimensional nature. As mentioned above, a fractional vortex corresponds to a multidimensional vector state residing in the infinite OAM Hilbert space. In this regard, our work can be connected to those reporting the so-called infinite dimensional entanglement exploited by half-integer spiral phase plates122324. Besides, if we further make Q vary continually in a certain range (by adjusting U), then we can image that the corresponding vector state traces a multidimensional curve in that space.

Another characteristic feature of equation (6) is the dependence of joint probability P, equivalently , on the parameters and , which suggests inversely that it is possible for us to retrieve and by measuring and analyzing the coincidence rates. Based on equation (6), we plot a group of coincidence curves for different degrees of polarization entanglement in Fig. 2, where can be easily obtained by adjusting the partial polarizer introduced in Fig. 1 (a). Besides, under the thin crystal approximation, the expression of spiral spectrum can be analytically given by11,

Figure 2

The joint probability amplitude vs. applied voltage (in Q) for different . Inset (a) shows the original entangled OAM spectrum based on equation (7).

where , is the Gauss hypergeometric function, γ (γ) denotes the ratio of the pump beam waist w to the signal (idler) beam waist w (w), p and p are the radial mode indices of signal and idler beam, respectively. Thus, the spectrum, , as a sum over all the p and p indices for each mode, is written as, .

As illustrated in Fig. 2, each coincidence curve is symmetrical, as a result of the symmetry of spiral spectrum given by equation (7). Besides, from a careful investigation on equation (6), we find that parameter can be directly retrieved from the following equation,

where , denotes the coincidence rate normalized by that with . The above equation can be confirmed by the plotted coincidence curves at Q = 0, see Fig. 2: , , , namely, , , . These results are in good agreement with predictions given by equation (8), as shown in Fig. 3, where the red circles and the blue line are calculated based on equations (6) and (8), respectively. A little surprisingly, regardless of no use of any polarizers in the measurement stage, our EO analyzers are capable to characterize the polarization entanglement, owing to quantum information transferred from polarization to OAM.

Figure 3

The normalized joint probability amplitude vs. the tunable degree of entanglement for Q = 0.

Furthermore, we can retrieve based on equation (6), and illustrate a typical example of spiral spectrum in Fig. 4, which is retrieved by using an 41-point electro-optic sampling, where (upper row) and (lower row) are trivially considered. Fig. 4 (a,d) indicate the sampling points of applied voltages (in terms of Q, as Q = αU). By means of electro-optic sampling and the matrix operations using Cramer’s rule25, we show the retrieved (solid bars) and original (empty bars) spiral spectra () in Fig. 4 (b,e). After making comparisons between the retrieved and original spectra, we illustrate the relative errors of retrieved spiral spectrum in Fig. 4 (c,f). And we find that of the lower OAM modes () are well retrieved, but with an increasing , will have a larger deviation from the expected, which is obviously due to the very limited sampling points selected in Fig. 4 (a). As is well known, the entangled spiral spectrum generated by SPDC is generally of limited bandwidth; photon pairs with a smaller OAM are more frequent than the higher-order ones, namely, if i < j. So it can be expected mathematically that a more precise spiral spectrum containing more OAM modes could be retrieved by our EO analyzer, only if the number of electro-optic sampling points increases enough.

Figure 4

Results for (upper row) and (lower row). (a,d) The joint probability amplitude vs. applied voltage (in Q), where 41 sampling points are marked by circles. (b,e) The retrieved spiral spectrum (solid bars), in comparison with the original one (empty bars). (c,f) The relative errors for OAM modes of .

Discussion

In an actual experiment, there are several reasons that can reduce the signal-to-noise ratio, such as the imperfections in the SPPs orientation and the misalignment of the SPPs. For the SPPs with non-integer topological charges, there is a radial phase discontinuity that characterizes its orientation. In our simulations, we have assumed that the SPPs in two arms are both placed radically along the x-axis direction. However, it is noted that a slight mismatch of the SPP orientations, β, will cause the change of non-integer OAM states. This can be understood well by decomposing the non-integer vortex into a superposition of integer OAM states, namely, . One can see that this imperfection will bring a phase shift of to each OAM mode26. Thus the joint probability of Eq. (6), and therefore, the performance of retrieve algorithm of Eqs. (9) and (10) will be both affected. Besides, as the OAM state measurements are highly symmetric, the good alignment of the SPPs in the optical paths is highly desirable. It was found that both lateral displacement and angular deflection will cause the coupling between different OAM and therefore broaden the OAM spectrum of a nonlocal vortex27. Accordingly, the measurement state with non-integer SPPs will be changed also in our case. In other words, both the imperfect orientation and misalignment of SPPs will affect the performance of our scheme, such as a reduction of the signal-to-noise ratio. These practical factors certainly deserve our very careful attention in the further experimental study.

In summary, we have demonstrated a theoretical scheme to analyze and retrieve the information of two-photon spin-orbit hyperentanglement. Our method is based on electro-optic sampling using a proposed hyperentanglement analyzer and simple matrix operation using Cramer’s rule. It is noted that an arbitrary pump beam will give birth to a more complicated OAM spectrum that is characterized by its complex elements, , where and are arbitrary integers28. This case also deserves our further investigation. From the experimental point of view, it is noted that a variety of high-quality SPPs have been made with the state-of-the-art micro-machining and a molding process29. Besides, spiral phase mirrors were also produced by direct machining with a diamond tool30, and this technique can be suitably adopted to make such an OAM generator. So we think that our proposed scheme of electro-optic hyper-entanglement analyzer is experimentally feasible with current lab technology, and may find potential application in designing a fast electro-optic switching for OAM-based quantum information processing.

Methods

For spiral spectrum analysis, we can retrieve in the following procedure: 1) Pick the joint probability amplitude for each sampling point , where (K + 1 is total number of sampling points). 2) Using these collected data of vs. Q, we can construct a system of K + 1 linear equations based on equation (6), as follows:

3) The above linear system can be expressed by a matrix equation, AX = b, where x is the column vector of unknown variables of , b is the column vector of , and A is the matrix of the coefficients of . The element of matrix A can be specified by,

where and . According to Cramer rule, we know if K = 2L and A is nonsingular, then the linear system of equation (9) has a unique solution, x = A−1b, where A−1 denotes the inverse of A. In other words, a desired spiral spectrum with ranging from −L to L can be retrieved using electro-optic sampling by our EO analyzers.

Additional Information

How to cite this article: Wu, Z. and Chen, L. Electro-optic analyzer of angular momentum hyperentanglement. Sci. Rep. 6, 21856; doi: 10.1038/srep21856 (2016).