• Photonics Research
  • Vol. 10, Issue 3, 777 (2022)
Li Zeng1、2、†, Rong Ma1、2、3、†, Hong Wen1、2, Meihong Wang1、2, Jun Liu1、2, Zhongzhong Qin1、2、4, and Xiaolong Su1、2、*
Author Affiliations
  • 1State Key Laboratory of Quantum Optics and Quantum Optics Devices, Institute of Opto-Electronics, Shanxi University, Taiyuan 030006, China
  • 2Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan 030006, China
  • 3College of Physics and Electronic Engineering, Shanxi University, Taiyuan 030006, China
  • 4e-mail: zzqin@sxu.edu.cn
  • show less
    DOI: 10.1364/PRJ.442925 Cite this Article Set citation alerts
    Li Zeng, Rong Ma, Hong Wen, Meihong Wang, Jun Liu, Zhongzhong Qin, Xiaolong Su. Deterministic distribution of orbital angular momentum multiplexed continuous-variable entanglement and quantum steering[J]. Photonics Research, 2022, 10(3): 777 Copy Citation Text show less

    Abstract

    Orbital angular momentum (OAM) multiplexing provides an efficient method to improve data-carrying capacity in various quantum communication protocols. It is a precondition to distribute OAM multiplexed quantum resources in quantum channels for implementing quantum communication. However, quantum steering of OAM multiplexed optical fields and the effect of channel noise on OAM multiplexed quantum resources remain unclear. Here, we generate OAM multiplexed continuous-variable (CV) entangled states and distribute them in lossy or noisy channels. We show that the decoherence property of entanglement and quantum steering of the OAM multiplexed states carrying topological charges l=1 and l=2 are the same as that of the Gaussian mode with l=0 in lossy and noisy channels. The sudden death of entanglement and quantum steering of high-order OAM multiplexed states is observed in the presence of excess noise. Our results demonstrate the feasibility to realize high data-carrying capacity quantum information processing by utilizing OAM multiplexed CV entangled states.

    1. INTRODUCTION

    Einstein–Podolsky–Rosen (EPR) entanglement plays a crucial role in quantum information processing, such as quantum communication, quantum computation, and quantum precision measurement [15]. Besides entanglement, quantum steering, which stands between entanglement [1] and Bell nonlocality [6] in the hierarchy of quantum correlations [7], has been identified as a useful quantum resource. Different from entanglement and Bell nonlocality, quantum steering shows unique asymmetry or even one-way characteristics [816] and, thus, allows asymmetric quantum information processing. For example, quantum steering enables one-side device-independent quantum key distribution [1719].

    Multiplexing provides an efficient method to enhance the data-carrying capability in both classical and quantum communication systems by combining multiple channels into a single channel. By utilizing different degrees of freedom (DOFs) of light, such as wavelength [20,21], polarization [22], temporal [2325] or spatial [26,27] modes, different types of multiplexing can be realized. Orbital angular momentum (OAM) of light [28] has also been found to be an attractive DOF to realize multiplexing due to its infinite range of possibly achievable topological charges [29,30]. OAM has found applications in discrete-variable quantum information processing, such as high-dimensional OAM entanglement generation [31], and 18-qubit entanglement with six photons’ three DOFs including OAM [32].

    Four-wave mixing (FWM) process in warm alkali vapor cell has found a wide range of applications [3339]. Especially, spatial-multi-mode advantage of the FWM process, attributed to its cavity-free configuration, makes it an ideal optical parametric amplifier to generate entangled images [35] and reconfigurable multipartite entanglement [36]. Quantum correlated twin beams carrying OAM were generated based on the FWM process in rubidium vapor [37]. OAM multiplexed bipartite and multipartite continuous-variable (CV) entangled states have been generated based on the FWM process [4042]. Furthermore, OAM multiplexed deterministic all-optical quantum teleportation has also been demonstrated by utilizing the OAM multiplexed bipartite CV entangled state generated from the FWM process [43]. To enhance the data-carrying capacity in quantum communication based on OAM multiplexed CV entangled states, it is essential to distribute them in lossy and noisy quantum channels towards practical applications. The distributions of weak coherent field and single photons carrying OAM in fiber, free space, and underwater have been experimentally investigated [4446]. However, it remains unclear whether the quantum entanglement and steering of OAM multiplexed CV entangled states are more sensitive to loss and noise than a commonly used Gaussian mode with l=0.

    Here, we present the deterministic distribution of OAM multiplexed CV quantum entanglement and steering in lossy and noisy channels. In the experiment, the OAM multiplexed entangled fields are generated deterministically based on the FWM process in warm cesium vapor and distributed deterministically in quantum channels. We show that the CV entangled states carrying topological charges l=1 and l=2 are as robust against loss as the Gaussian mode with l=0. Sudden death of entanglement and quantum steering of a high-order OAM multiplexed CV entangled state is observed in the presence of noise. Our results pave the way for applying OAM multiplexed CV entanglement and quantum steering in high data-carrying capacity quantum communication.

    2. PRINCIPLE AND EXPERIMENTAL SETUP

    Figure 1(a) shows the schematic of the experimental setup, and Fig. 1(b) shows the double-Λ energy-level structure used for the FWM process, which is formed from the Cs133D1 line with an excited level (6P1/2,F=4) and two ground levels (6S1/2,F=3 and F=4). The pump beam is about 1.6-GHz blue detuned from 6S1/2,F=36P1/2,F=4 transition, and the probe beam is 9.2-GHz redshifted relative to the pump beam. The pump and probe beams are combined by a GL and then cross each other in the center of the cesium vapor cell at an angle of 6 mrad [47]. The gain of the FWM process is around 3 with a pump power of 240 mW and a probe power of 3 μW. By injecting the probe beam carrying topological charge l of the OAM mode, a conjugate beam carrying topological charge l of the OAM mode is generated on the other side of the pump, which satisfies OAM conservation in the FWM process. The topological charge of the OAM mode l=1 or l=2 is added to the probe beam by passing it through a VPP. The pump beam is filtered out by using a GT with an extinction ratio of 105:1 after the vapor cell.

    (a) Experimental setup for the generation and distribution of OAM multiplexed CV quantum entanglement and steering in a lossy or noisy channel. Pr, probe beam; Conj, conjugate beam; LO−l,P and LOl,C, local oscillators of Pr and Conj fields; AM, amplitude modulator; PM, phase modulator; GL, Glan-laser polarizer; GT, Glan–Thompson polarizer; PBS, polarization beam splitter; HWP, half-wave plate; VPP, vortex phase plate; M, mirror; BS, 50:50 beam splitter; BHD1, BHD2, balanced homodyne detectors; SA, spectrum analyzer. (b) Double-Λ energy-level structure for the FWM process in a cesium vapor cell. Δ, one-photon detuning.

    Figure 1.(a) Experimental setup for the generation and distribution of OAM multiplexed CV quantum entanglement and steering in a lossy or noisy channel. Pr, probe beam; Conj, conjugate beam; LOl,P and LOl,C, local oscillators of Pr and Conj fields; AM, amplitude modulator; PM, phase modulator; GL, Glan-laser polarizer; GT, Glan–Thompson polarizer; PBS, polarization beam splitter; HWP, half-wave plate; VPP, vortex phase plate; M, mirror; BS, 50:50 beam splitter; BHD1, BHD2, balanced homodyne detectors; SA, spectrum analyzer. (b) Double-Λ energy-level structure for the FWM process in a cesium vapor cell. Δ, one-photon detuning.

    The Conj field is kept by Alice, whereas the Pr field is distributed to a remote quantum node owned by Bob through a lossy or noisy channel. The lossy channel is simulated by an HWP and a PBS. The noisy channel is modeled by combining the Pr field with an auxiliary beam at a PBS followed by an HWP and a PBS. The auxiliary beam carries the same frequency and topological charge as the Pr field and is modulated by an AM and a PM with white noise [16]. To characterize the OAM multiplexed CV entangled state, its covariance matrix is experimentally measured by utilizing two sets of BHDs. In order to extract the CV quadrature information carried by the OAM mode with a topological charge l, an LO with opposite topological charge l is required. In our experiment, the spatially mode-matched LO beams used in the BHDs are obtained from a second set of FWM processes which is around 5 mm above the first set of FWM processes in the same vapor cell [35]. More details of the experimental parameters can be found in Appendix B.

    The Hamiltonian of the OAM multiplexed FWM process can be expressed as [40] H^=liγla^l,Pa^l,C+H.c.,where γl, which is defined as the interaction strength of each OAM pair, is proportional to the pump power. a^l,P and a^l,C are the creation operators related to OAM modes of the Pr and Conj fields, respectively, and H.c. represents the Hermitian conjugate. Since the pump beam does not carry OAM (l=0), the topological charges of the Pr and Conj fields are opposite due to OAM conservation in the FWM process. The output state of the OAM multiplexed FWM process is as follows: |Ψout=|Ψl|Ψ0|Ψl,where |Ψl=|ψl,P,ψl,C presents a series of independent OAM multiplexed CV EPR entangled states generated in the FWM process. |ψl,P and |ψl,C represent Pr field carrying topological charge l and Conj field carrying topological charge l, respectively.

    All Gaussian properties of the CV entangled state |Ψl can be determined by the covariance matrix σAB with matrix element σij=ξ^iξ^j+ξ^jξ^i/2ξ^iξ^j, where ξ^(X^l,C,Y^l,C,X^l,P,Y^l,P)T, X^=a^+a^, and Y^=(a^a^)/i represent amplitude and phase quadratures of the Conj and Pr fields, respectively. The covariance matrix of the OAM multiplexed entangled state after distribution in a lossy or noisy channel is as follows: σAB,δ,η=(Va0Vc00Va0VcVc0Vb00Vc0Vb),with Va=V+V2, Vb=ηV+V2+(1η)(1+δ), and Vc=ηVV2. V and V represent the variances of squeezed and antisqueezed quadratures of the optical mode, respectively. The derivation of the covariance matrix can be found in Appendix A. For pure squeezed states, VV=1, whereas VV>1 after pure squeezed states suffer from loss or noise. η and δ represent transmission efficiency and excess noise of the quantum channel, respectively. We have δ=0 in a lossy channel, whereas we have δ>0 in a noisy channel. The submatrices σA=VaI and σB=VbI with I as the identity matrix are the covariance matrices corresponding to the states of Alice’s and Bob’s subsystems, respectively. The CV entangled state is a symmetric state and an asymmetric state when σA=σB and σAσB, respectively.

    The Peres–Horodecki criterion of positivity under the partial transpose (PPT) criterion is a sufficient and necessary criterion to characterize the entanglement of CV bipartite entanglement [48]. If the smallest symplectic eigenvalue ν of the partially transposed covariance matrix is smaller than 1, bipartite entanglement exists. Otherwise, it is a separable state. Furthermore, smaller ν represents stronger entanglement.

    Quantum steering for bipartite Gaussian states of CV systems can be quantified by [49] GAB(σAB)=max{0,12lndetσAdetσAB},where GAB(σAB)>0 represents that Alice has the ability to steer Bob’s state. Similarly, we have GBA(σAB)=max{0,12lndetσBdetσAB},which represents Bob’s ability to steer Alice’s state. From the expressions GAB(σAB) and GBA(σAB), it can be seen that Alice and Bob have the same steerability if detσA=detσB is satisfied, i.e., the bipartite Gaussian state is a symmetric state. If the state is an asymmetric state, the steerabilities of Alice and Bob will be different.

    3. RESULTS

    To verify the OAM property of the optical fields, we measure the spatial beam patterns of quantum states |Ψ1 and |Ψ2 transmitted through a lossy channel, which are shown in the top rows of Figs. 2(a) and 2(b), respectively. It is obvious that the Pr and Conj fields are both Laguerre–Gaussian beams. To infer their topological charges, they are passed through a tilted lens and imaged on a camera. As shown in the bottom rows of Figs. 2(a) and 2(b), the number of dark stripes gives the number of the topological charge, and the direction gives its sign [50]. As the transmission efficiency of the Pr field decreases, its optical intensity also decreases, whereas its topological charge remains unchanged. Additional beam patterns of the Pr and Conj fields can be found in Appendix D.

    (a) and (b) Beam patterns of the OAM multiplexed CV entanglement for l=1 and l=2 in a lossy channel and corresponding transmitted patterns through a tilted lens. (c) Dependence of PPT values of the OAM multiplexed CV entanglement on transmission efficiency η for l=0, l=1, and l=2 in lossy channels. The PPT value is 0.46±0.01 at η=1. Curves and data points show theoretical predictions and experimental results, respectively. Error bars of experimental data represent one standard deviation and are obtained based on the statistics of the measured data.

    Figure 2.(a) and (b) Beam patterns of the OAM multiplexed CV entanglement for l=1 and l=2 in a lossy channel and corresponding transmitted patterns through a tilted lens. (c) Dependence of PPT values of the OAM multiplexed CV entanglement on transmission efficiency η for l=0, l=1, and l=2 in lossy channels. The PPT value is 0.46±0.01 at η=1. Curves and data points show theoretical predictions and experimental results, respectively. Error bars of experimental data represent one standard deviation and are obtained based on the statistics of the measured data.

    The covariance matrices of the OAM multiplexed entangled states are reconstructed by measuring the noise variances of the amplitude and phase quadratures of the Conj and Pr fields Δ2X^l,C, Δ2Y^l,C, Δ2X^l,P, Δ2Y^l,P, as well as their correlation variances of amplitude and phase quadratures Δ2(X^l,PX^l,C) and Δ2(P^l,P+P^l,C), respectively. Details about the measurement of the covariance matrices can be found in Appendix C. Based on the covariance matrix of each OAM multiplexed entangled state at different loss and noise levels, its quantum entanglement and quantum steering characteristics are evaluated experimentally.

    Figure 2(c) shows the dependence of PPT values of the CV bipartite entangled state carrying different topological charges on the transmission efficiency of the Pr field. The correlation and anticorrelation levels of the initial CV entangled states carrying topological charges l=0, l=1, and l=2 are all around 3.3 and 6.1 dB, which correspond to V=0.47 and V=4.11, respectively. The entanglement between Pr and Conj fields degrades as the transmission efficiency decreases. However, the entanglement is robust against loss, i.e., it always exists until the transmission efficiency reaches 0. It is obvious that the CV bipartite entangled state carrying topological charges l=1 and l=2 is as robust to loss as its Gaussian counterpart l=0.

    Figure 3 shows the dependence of PPT values of the CV bipartite entangled state carrying different topological charges in noisy channels. Compared with the results in Fig. 2(c), the entanglement disappears at a certain transmission efficiency of the Pr field in the presence of excess noise, which demonstrates the sudden death of CV quantum entanglement. Furthermore, the higher the excess noise is, the sooner entanglement disappears. The transmission efficiencies where entanglement starts to disappear are η=0.10,0.28, and 0.44, respectively, for the excess noise δ=0.15,0.5, and 1 in the units of shot-noise level (SNL). We show that OAM multiplexed CV entangled states carrying high-order topological charges l=1 and l=2 exhibit the same decoherence tendencies as their Gaussian counterpart l=0 in noisy channels.

    Dependence of PPT values of the OAM multiplexed CV entanglement on transmission efficiency η for l=0, l=1, and l=2 in noisy channels. Three different amounts of excess noise δ=0.15 (black), δ=0.5 (red), and δ=1 (blue) are compared. The light blue plane shows the boundary for sudden death of entanglement where the PPT value equals 1. The three vertical dashed lines indicate corresponding transmission efficiencies where entanglement starts to disappear. Curves and data points show theoretical predictions and experimental results, respectively. Error bars of experimental data represent one standard deviation and are obtained based on the statistics of the measured data.

    Figure 3.Dependence of PPT values of the OAM multiplexed CV entanglement on transmission efficiency η for l=0, l=1, and l=2 in noisy channels. Three different amounts of excess noise δ=0.15 (black), δ=0.5 (red), and δ=1 (blue) are compared. The light blue plane shows the boundary for sudden death of entanglement where the PPT value equals 1. The three vertical dashed lines indicate corresponding transmission efficiencies where entanglement starts to disappear. Curves and data points show theoretical predictions and experimental results, respectively. Error bars of experimental data represent one standard deviation and are obtained based on the statistics of the measured data.

    The dependence of steerabilities GAB and GBA on the transmission efficiency η and topological charge l in lossy and noisy channels is shown in Figs. 4(a) and 4(b), respectively. In a lossy channel, the steerabilities for both directions always decrease when the transmission efficiency decreases. One-way steering is observed in the region of 0<η<0.72 for the OAM multiplexed CV entangled state carrying different topological charges l=0, l=1, and l=2. In a noisy channel where the excess noise δ=0.15 (in the units of SNL) exists, the steerabilities GAB and GBA are lower than those in the lossy channel [13,51]. Furthermore, Alice loses steerability in the region of 0<η<0.45, whereas Bob loses steerability in the region of 0<η<0.75, which confirms sudden death of quantum steering in a noisy channel [51]. It is worth noting that the CV entangled state carrying topological charges l=1 and l=2 has the same steerabilities as its counterpart l=0.

    Quantum steerabilities of OAM multiplexed CV entangled state distributed in a (a) lossy or (b) noisy channel. The excess noise shown in (b) is δ=0.15. Solid and dashed curves show theoretical predictions of GA→B and GB→A, respectively. Data points show experimental results. Error bars of experimental data represent one standard deviation and are obtained based on the statistics of the measured data.

    Figure 4.Quantum steerabilities of OAM multiplexed CV entangled state distributed in a (a) lossy or (b) noisy channel. The excess noise shown in (b) is δ=0.15. Solid and dashed curves show theoretical predictions of GAB and GBA, respectively. Data points show experimental results. Error bars of experimental data represent one standard deviation and are obtained based on the statistics of the measured data.

    4. CONCLUSION

    The distribution of OAM multiplexed CV entanglement and quantum steering in quantum channels with homogeneous loss and noise, such as fiber channels, was experimentally simulated in our paper. There were also other quantum channels with inhomogeneous loss and noise, such as atmospheric turbulence and diffraction. Recently, it was shown that other optical fields carrying OAM, such as vector beams, were turbulence resilient in atmospheric turbulence [52]. Thus, it is worthwhile to investigate the turbulence-resilient characteristics of OAM multiplexed CV quantum entanglement and steering, which have the potential to substantially improve the quantum communication distance and fidelity.

    To summarize, we experimentally demonstrated quantum steering of OAM multiplexed optical fields and investigated the distribution of OAM multiplexed CV entanglement and quantum steering in quantum channels. We showed that the decoherence property of CV entanglement and quantum steering of the OAM multiplexed optical fields carrying topological charges l=1 and l=2 were the same as that of the counterpart Gaussian mode with l=0 in lossy and noisy channels. The sudden death of entanglement and quantum steering of high-order OAM multiplexed optical fields was observed in the presence of excess noise. Our results demonstrated the feasibility to improve the quantum communication capacity in practical quantum channels by utilizing OAM multiplexed CV entanglement and quantum steering.

    APPENDIX A: THEORETICAL MODEL

    The Hamiltonian of the OAM multiplexed FWM process can be expressed as H^=likla^0,Pumpa^0,Pumpa^l,Pa^l,Ceiθ+H.c.,where a^0,Pump, a^l,P, and a^l,C are the annihilation operator of the pump field and the creation operators related to OAM modes of the Pr and Conj fields, respectively. θ is the phase of the pump field. The first term represents the process in which two pump photons are converted into a Pr photon and a Conj photon. kl is the interaction strength of each OAM pair. Since the pump beam does not carry OAM (l=0), the topological charges of the Pr and Conj fields are opposite due to OAM conservation in the FWM process. The second term H.c. is the H.c. of the first term and represents the reversed process in which a Pr photon and a Conj photon are converted to two pump photons.

    The pump field is much stronger than the Pr and Conj fields in the FWM process, so it can be regarded as a classical field. By combining the intensity of the pump field |αPump|2 with kl, i.e., |αPump|2kl=γl, and taking θ=0, the Hamiltonian can be simplified as [40] H^=liγla^l,Pa^l,C+H.c.

    The output state of the OAM multiplexed FWM process is as follows: |Ψout=|Ψl|Ψ0|Ψl,where |Ψl=S^(γl)|vac=|ψl,P,ψl,C means the EPR entangled state with Pr and Conj fields carrying topological charges l and l, respectively. |vac represents the vacuum state, and S^(γl)=eγl(a^l,Pa^l,Ca^l,Pa^l,C) represents the two-mode squeezing operator for the vacuum state carrying OAM topological charge l. It is obvious that a series of CV EPR entangled states carrying independent topological charges is generated in the FWM process; i.e., OAM multiplexing is realized.

    All Gaussian properties of the CV Gaussian entangled state |Ψl can be determined by its covariance matrix σAB with the matrix element σij=ξ^iξ^j+ξ^jξ^i/2ξ^iξ^j, where ξ^(X^l,C,Y^l,C,X^l,P,Y^l,P)T, and X^=a^+a^ and Y^=(a^a^)/i represent amplitude and phase quadratures of the Conj and Pr fields.

    The covariance matrix of the CV bipartite entangled state can be written as σAB=(V+V20VV200V+V20VV2VV20V+V200VV20V+V2)=(V+V2IVV2ZVV2ZV+V2I),where V and V represent the variances of squeezed and antisqueezed quadratures of the optical mode, respectively. I and Z are the Pauli matrices, I=(1001),Z=(1001).

    Then, we consider the distribution of CV entangled state |Ψl in a lossy or noisy channel. Let a^l,P and a^l,C represent the annihilation operators of the Pr and Conj fields, respectively. After the Pr field a^l,P is distributed in a lossy channel, it becomes a^l,P=ηa^l,P+1ημ^, where μ^ represents the vacuum state with a variance of 1. Similarly, the Pr field becomes a^l,P=ηa^l,P+1η(ϵ^+μ^) after it is distributed in a noisy channel with excess noise Δ2(X^ϵ)=Δ2(Y^ϵ)=δ [13]. δ=0 means that there is no excess noise, and only loss exists in the channel. δ>0 means that there exists excess noise in the channel. So the covariance matrix of CV entangled state |Ψl after distribution in the lossy or noisy channel is as follows: σAB,δ,η=(Va0Vc00Va0VcVc0Vb00Vc0Vb)=(VaIVcZVcZVbI),with Va=V+V2, Vb=ηV+V2+(1η)(1+δ), and Vc=ηVV2.

    APPENDIX B: DETAILS OF THE EXPERIMENT

    The Ti:sapphire laser (Coherent MBR-110) is about 1.6 GHz blue detuned from the Cs133 D1 line 6S1/2,F=36P1/2,F=4 transition with a total power of 1.2 W. As shown in Fig. 5, the laser beam is split into two beams by PBS1. The horizontally polarized beam is split into two beams by PBS2, and they serve as the pump beams for two sets of FWM processes in the same cesium vapor cell. The vertically polarized beam with a power of 30 mW passes through a resonance EOM (Qubig GmbH PM-Cs) and three successive temperature-stabilized etalons to realize a 9.2-GHz frequency redshift [47]. The OAM of topological charge l=1 or l=2 is added to the frequency redshifted beam by passing it through a VPP (RPC Photonics). Then the frequency redshifted beam is split into three beams, among which the first two beams serve as the probe beams of two FWM processes, and the third beam serves as the auxiliary beam for the noisy channel. The pump and probe beams are combined in a GL and then cross each other in the center of the cesium vapor cell (Traid Technology, Inc.) at an angle of 6 mrad. The vapor cell is 25-mm long, and its temperature is stabilized at 103°C.

    Detailed experimental schematic for distributing OAM multiplexed CV entanglement in a noisy channel. The lossy channel is realized by blocking the auxiliary beam. D-shaped mirrors (DMs) are utilized to combine or separate light beams with small distances. HWP, half-wave plate; PBS, polarization beam splitter; EOM, electro-optic modulator; VPP, vortex phase plate; GL, Glan-laser polarizer; GT, Glan–Thompson polarizer; Pr, probe beam; Conj, conjugate beam; AM, amplitude modulator; PM, phase modulator; M, mirror; DM, D-shaped mirror; PZT, piezoelectric ceramics; BS, 50:50 beam splitter; BHD, balanced homodyne detector; SA, spectrum analyzer.

    Figure 5.Detailed experimental schematic for distributing OAM multiplexed CV entanglement in a noisy channel. The lossy channel is realized by blocking the auxiliary beam. D-shaped mirrors (DMs) are utilized to combine or separate light beams with small distances. HWP, half-wave plate; PBS, polarization beam splitter; EOM, electro-optic modulator; VPP, vortex phase plate; GL, Glan-laser polarizer; GT, Glan–Thompson polarizer; Pr, probe beam; Conj, conjugate beam; AM, amplitude modulator; PM, phase modulator; M, mirror; DM, D-shaped mirror; PZT, piezoelectric ceramics; BS, 50:50 beam splitter; BHD, balanced homodyne detector; SA, spectrum analyzer.

    The two sets of FWM processes are constructed in the same cesium vapor cell with a height difference of 5 mm. The bottom FWM process is used to generate the OAM multiplexed CV entangled state, whereas the top FWM process is used to generate spatially matched LOs with Pr and Conj fields. The pump power in the bottom FWM process for generating OAM multiplexed CV entanglement is 240 mW. The probe gain of the bottom FWM process is around 3, and the degree of initially generated CV entanglement is around 3.3±0.1  dB. The pump power and seed probe power of the top FWM processes are 450 mW and 100 μW, respectively, so that the SNL is around 10 dB higher than the electronic noise of the homodyne detector. The bottom FWM process is weakly seeded with a probe power of around 3 μW for relative phase locking of the Pr/Conj fields and their LOs in the balanced homodyne detections.

    The lossy channel is simulated by HWP10 and PBS6. The noisy channel is modeled by combining the vertically polarized Pr field with a horizontally polarized auxiliary beam at PBS5 followed by HWP10 and PBS6. The auxiliary beam carries the same frequency and topological charge as the Pr field and is modulated by an AM and a PM with excess noise at 1.2 MHz. The amount of excess noise is adjusted by tuning the amplitude of the signal applied to the AM and PM and evaluated in the units of SNL. For example, excess noise δ=1 corresponds to a noise level that is 3 dB higher than the SNL. By tuning HWP10, the lower the transmission efficiency of the Pr field, the higher the excess noise coupled to the Pr field. In practical quantum communication protocols, higher excess noise is coupled to the quantum entangled state as the communication distance increases, accompanying lower transmission efficiency. Therefore, our experimental setting is similar to the realistic scenarios in the practical noisy quantum channel. To characterize the OAM multiplexed CV entangled state, its covariance matrix is experimentally obtained by utilizing two sets of BHDs (Thorlabs PDB450A). The interference visibilities for the two sets of BHDs are both around 99%. The electrical gains of these two BHDs are both 105  V/A. The original photodiodes are replaced by high quantum efficiency (QE) photodiodes with QE=98% at 895 nm (first sensor). To measure amplitude quadrature X^ or phase quadrature Y^ of the Pr/Conj fields, the relative phases between them and their LOs are locked by applying the feedback signal from proportional-integral–derivative circuits and high-voltage amplifiers to PZTs.

    APPENDIX C: MEASUREMENT OF COVARIANCE MATRIX

    To reconstruct the covariance matrix of the CV quantum entangled state, we perform six different measurements on the output optical modes. These measurements include the variances of the amplitude and phase quadratures of Conj and Pr fields Δ2X^l,C, Δ2Y^l,C, Δ2X^l,P, Δ2Y^l,P, as well as noise variances of their joint amplitude and phase quadratures Δ2(X^l,PX^l,C) and Δ2(Y^l,P+Y^l,C), respectively. Δ2X^l,C and Δ2Y^l,C (Δ2X^l,P and Δ2Y^l,P) are experimentally measured by locking the relative phase of the Conj (Pr) field, and its corresponding local oscillator of BHD1 (BHD2) at amplitude quadrature or phase quadrature. Δ2(X^l,PX^l,C) and Δ2(Y^l,P+Y^l,C) are experimentally measured by locking the relative phases of BHD1 and BHD2 at amplitude quadrature or phase quadrature and then subtracting or adding the photocurrents with a radio-frequency subtractor or adder. SNL is achieved by blocking the Pr and/or Conj fields so that only the noise of vacuum is measured. The settings of the SA (Agilent E4411B) are 30-kHz resolution bandwidth, 100-Hz video bandwidth, and zero span at 1.2 MHz.

    Figure 6 shows the measured six noise variance levels of the CV entanglement of the OAM multiplexed CV entangled state before being distributed in lossy and noisy channels (δ=0,η=1). As shown in Fig. 6, the noise levels of amplitude and phase quadratures are Δ2X^l,C=Δ2Y^l,C=Δ2X^l,P=Δ2Y^l,P=3.6±0.1  dB, and the noise levels of correlated quadratures are Δ2(X^l,PX^l,C)=Δ2(Y^l,P+Y^l,C)=3.3±0.1  dB. Our FWM process works in the amplification regime, so amplitude quadratures and phase quadratures of the Pr/Conj fields show strong correlations and anticorrelations, respectively. It is clear that the degrees of entanglement for CV entangled states carrying topological charges l=1 and l=2 are close to that of their spatially Gaussian counterpart l=0.

    Measured quantum correlation noises for initially generated OAM multiplexed CV entangled states carrying topological charges (a) l=0, (b) l=1, and (c) l=2, respectively. Brown curve at 0 dB shows the SNL. The other six curves show the noise variances of Δ2X^−l,C, Δ2Y^−l,C, Δ2X^l,P, Δ2Y^l,P, as well as the noise variances of their joint amplitude or phase quadrature Δ2(X^l,P−X^−l,C) and Δ2(Y^l,P+Y^−l,C). These six curves are all normalized to the same SNL. All the measurements are performed at 1.2 MHz. The electronic noise of the BHDs and the background noise from leaked pump fields are subtracted from the SNL and signals, respectively.

    Figure 6.Measured quantum correlation noises for initially generated OAM multiplexed CV entangled states carrying topological charges (a) l=0, (b) l=1, and (c) l=2, respectively. Brown curve at 0 dB shows the SNL. The other six curves show the noise variances of Δ2X^l,C, Δ2Y^l,C, Δ2X^l,P, Δ2Y^l,P, as well as the noise variances of their joint amplitude or phase quadrature Δ2(X^l,PX^l,C) and Δ2(Y^l,P+Y^l,C). These six curves are all normalized to the same SNL. All the measurements are performed at 1.2 MHz. The electronic noise of the BHDs and the background noise from leaked pump fields are subtracted from the SNL and signals, respectively.

    With the measured six noise variances, the cross-correlation matrix elements are calculated via Cov(ξ^i,ξ^j)=12[Δ2(ξ^i+ξ^j)Δ2ξ^iΔ2ξ^j],Cov(ξ^i,ξ^j)=12[Δ2(ξ^iξ^j)Δ2ξ^iΔ2ξ^j].

    In the experiment, we obtain all the covariance matrices of quantum states with different transmission efficiencies and excess noise and then calculate the smallest symplectic eigenvalue ν of the partially transposed covariance matrices GAB(σAB) and GBA(σAB) to verify whether quantum entanglement and quantum steering exist.

    APPENDIX D: SUPPLEMENTAL BEAM PATTERNS

    Supplemental beam patterns of the OAM modes of Pr and Conj beams generated from the FWM process are shown in Fig. 7. Pr and Conj beams carrying different topological charges l=0, l=1, and l=2 are generated by removing or switching different VPPs, and their beam patterns are taken after the GT with a camera. It is obvious that the dark hollow patterns of the Laguerre–Gaussian beams are enlarged as the topological charge increases. Furthermore, the beam patterns of Pr and Conj fields are symmetric with respect to the pump beam.

    Images of OAM modes of the Pr beam and Conj beam generated from the FWM process. (a) l=0; (b) l=1; (c) l=2. From left to right: Conj beam, leaked pump beam, and Pr beam.

    Figure 7.Images of OAM modes of the Pr beam and Conj beam generated from the FWM process. (a) l=0; (b) l=1; (c) l=2. From left to right: Conj beam, leaked pump beam, and Pr beam.

    References

    [1] R. Horodecki, P. Horodecki, M. Horodecki, K. Horodecki. Quantum entanglement. Rev. Mod. Phys., 81, 865-942(2009).

    [2] S. L. Braunstein, P. van Loock. Quantum information with continuous variables. Rev. Mod. Phys., 77, 513-577(2005).

    [3] H. J. Kimble. The quantum internet. Nature, 453, 1023-1030(2008).

    [4] C. Weedbrook, S. Pirandola, R. García-Patrón, N. J. Cerf, T. C. Ralph, J. H. Shapiro, S. Lloyd. Gaussian quantum information. Rev. Mod. Phys., 84, 621-669(2012).

    [5] X. Wang, T. Hiroshima, A. Tomita, M. Hayashi. Quantum information with Gaussian states. Phys. Rep., 448, 1-111(2007).

    [6] J. S. Bell. On the Einstein Podolsky Rosen paradox. Physics, 1, 195-200(1964).

    [7] R. Uola, A. C. S. Costa, H. C. Nguyen, O. Gühne. Quantum steering. Rev. Mod. Phys., 92, 015001(2020).

    [8] H. M. Wiseman, S. J. Jones, A. C. Doherty. Steering, entanglement, nonlocality, and the Einstein-Podolsky-Rosen paradox. Phys. Rev. Lett., 98, 140402(2007).

    [9] V. Händchen, T. Eberle, S. Steinlechner, A. Samblowski, T. Franz, R. F. Werner, R. Schnabel. Observation of one-way Einstein-Podolsky-Rosen steering. Nat. Photonics, 6, 596-599(2012).

    [10] S. Armstrong, M. Wang, R. Y. Teh, Q. Gong, Q. He, J. Janousek, H. A. Bachor, M. D. Reid, P. K. Lam. Multipartite Einstein-Podolsky-Rosen steering and genuine tripartite entanglement with optical networks. Nat. Phys., 11, 167-172(2015).

    [11] S. Wollmann, N. Walk, A. J. Bennet, H. M. Wiseman, G. J. Pryde. Observation of genuine one-way Einstein-Podolsky-Rosen steering. Phys. Rev. Lett., 116, 160403(2016).

    [12] K. Sun, X.-J. Ye, J.-S. Xu, X.-Y. Xu, J.-S. Tang, Y.-C. Wu, J.-L. Chen, C.-F. Li, G.-C. Guo. Experimental quantification of asymmetric Einstein-Podolsky-Rosen steering. Phys. Rev. Lett., 116, 160404(2016).

    [13] Z. Qin, X. Deng, C. Tian, M. Wang, X. Su, C. Xie, K. Peng. Manipulating the direction of Einstein-Podolsky-Rosen steering. Phys. Rev. A, 95, 052114(2017).

    [14] Y. Xiao, X.-J. Ye, K. Sun, J.-S. Xu, C.-F. Li, G.-C. Guo. Demonstration of multisetting one-way Einstein-Podolsky-Rosen steering in two-qubit systems. Phys. Rev. Lett., 118, 140404(2017).

    [15] A. Cavaillès, H. Le Jeannic, J. Raskop, G. Guccione, D. Markham, E. Diamanti, M. D. Shaw, V. B. Verma, S. W. Nam, J. Laurat. Demonstration of Einstein-Podolsky-Rosen steering using hybrid continuous- and discrete-variable entanglement of light. Phys. Rev. Lett., 121, 170403(2018).

    [16] M. Wang, Y. Xiang, H. Kang, D. Han, Y. Liu, Q. He, Q. Gong, X. Su, K. Peng. Deterministic distribution of multipartite entanglement and steering in a quantum network by separable states. Phys. Rev. Lett., 125, 260506(2020).

    [17] C. Branciard, E. G. Cavalcanti, S. P. Walborn, V. Scarani, H. M. Wiseman. One-sided device-independent quantum key distribution: security, feasibility, and the connection with steering. Phys. Rev. A, 85, 010301(2012).

    [18] T. Gehring, V. Händchen, J. Duhme, F. Furrer, T. Franz, C. Pacher, R. F. Werner, R. Schnabel. Implementation of continuous-variable quantum key distribution with composable and one-sided-device-independent security against coherent attacks. Nat. Commun., 6, 8795(2015).

    [19] N. Walk, S. Hosseini, J. Geng, O. Thearle, J. Y. Haw, S. Armstrong, S. M. Assad, J. Janousek, T. C. Ralph, T. Symul, H. M. Wiseman, P. K. Lam. Experimental demonstration of Gaussian protocols for one-sided device-independent quantum key distribution. Optica, 3, 634-642(2016).

    [20] M. Pysher, Y. Miwa, R. Shahrokhshahi, R. Bloomer, O. Pfister. Parallel generation of quadripartite cluster entanglement in the optical frequency comb. Phys. Rev. Lett., 107, 030505(2011).

    [21] J. Roslund, R. M. de Araújo, S. Jiang, C. Fabre, N. Treps. Wavelength-multiplexed quantum networks with ultrafast frequency combs. Nat. Photonics, 8, 109-112(2014).

    [22] H. Liu, J. Wang, H. Ma, S. Sun. Polarization-multiplexing-based measurement-device-independent quantum key distribution without phase reference calibration. Optica, 5, 902-909(2018).

    [23] S. Yokoyama, R. Ukai, S. C. Armstrong, C. Sornphiphatphong, T. Kaji, S. Suzuki, J. Yoshikawa, H. Yonezawa, N. C. Menicucci, A. Furusawa. Ultra-large-scale continuous-variable cluster states multiplexed in the time domain. Nat. Photonics, 7, 982-986(2013).

    [24] A. Marandi, Z. Wang, K. Takata, R. L. Byer, Y. Yamamoto. Network of time-multiplexed optical parametric oscillators as a coherent Ising machine. Nat. Photonics, 8, 937-942(2014).

    [25] M. Larsen, X. Guo, C. Breum, J. Neergaard-Nielsen, U. Andersen. Deterministic generation of a two-dimensional cluster state. Science, 366, 369-372(2019).

    [26] D. J. Richardson, J. M. Fini, L. E. Nelson. Space-division multiplexing in optical fibres. Nat. Photonics, 7, 354-362(2013).

    [27] Y. Pu, Y. Wu, N. Jiang, W. Chang, C. Li, S. Zhang, L. Duan. Experimental entanglement of 25 individually accessible atomic quantum interfaces. Sci. Adv., 4, eaar3931(2018).

    [28] L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, J. P. Woerdman. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes. Phys. Rev. A, 45, 8185-8189(1992).

    [29] J. Wang, J. Yang, I. M. Fazal, N. Ahmed, Y. Yan, H. Huang, Y. Ren, Y. Yue, S. Dolinar, M. Tur, A. E. Willner. Terabit free-space data transmission employing orbital angular momentum multiplexing. Nat. Photonics, 6, 488-496(2012).

    [30] N. Bozinovic, Y. Yue, Y. Ren, M. Tur, P. Kristensen, H. Huang, A. E. Willner, S. Ramachandran. Terabit-scale orbital angular momentum mode division multiplexing in fibers. Science, 340, 1545-1548(2013).

    [31] R. Fickler, G. Campbell, B. Buchler, P. K. Lam, A. Zeilinger. Quantum entanglement of angular momentum states with quantum numbers up to 10,010. Proc. Natl. Acad. Sci. USA, 113, 13642-13647(2016).

    [32] X. L. Wang, Y. H. Luo, H. L. Huang, M. C. Chen, Z. E. Su, C. Liu, C. Chen, W. Li, Y. Q. Fang, X. Jiang, J. Zhang, L. Li, N. L. Liu, C. Y. Lu, J. W. Pan. 18-qubit entanglement with six photons’ three degrees of freedom. Phys. Rev. Lett., 120, 260502(2018).

    [33] Z. Qin, A. S. Prasad, T. Brannan, A. MacRae, A. Lezama, A. I. Lvovsky. Complete temporal characterization of a single photon. Light Sci. Appl., 4, e298(2015).

    [34] Z. Qin, L. Cao, H. Wang, A. M. Marino, W. Zhang, J. Jing. Experimental generation of multiple quantum correlated beams from hot rubidium vapor. Phys. Rev. Lett., 113, 023602(2014).

    [35] V. Boyer, A. M. Marino, R. C. Pooser, P. D. Lett. Entangled images from four-wave mixing. Science, 321, 544-547(2008).

    [36] K. Zhang, W. Wang, S. Liu, X. Pan, J. Du, Y. Lou, S. Yu, S. Lv, N. Treps, C. Fabre, J. Jing. Reconfigurable hexapartite entanglement by spatially multiplexed four-wave mixing processes. Phys. Rev. Lett., 124, 090501(2020).

    [37] A. M. Marino, V. Boyer, R. C. Pooser, P. D. Lett, K. Lemons, K. M. Jones. Delocalized correlations in twin light beams with orbital angular momentum. Phys. Rev. Lett., 101, 093602(2008).

    [38] F. Hudelist, J. Kong, C. Liu, J. Jing, Z. Y. Ou, W. Zhang. Quantum metrology with parametric amplifier-based photon correlation interferometers. Nat. Commun., 5, 3049(2014).

    [39] R. C. Pooser, B. Lawrie. Ultrasensitive measurement of microcantilever displacement below the shot-noise limit. Optica, 2, 393-399(2015).

    [40] X. Pan, S. Yu, Y. Zhou, K. Zhang, K. Zhang, S. Lv, S. Li, W. Wang, J. Jing. Orbital-angular-momentum multiplexed continuous-variable entanglement from four-wave mixing in hot atomic vapor. Phys. Rev. Lett., 123, 070506(2019).

    [41] S. Li, X. Pan, Y. Ren, H. Liu, S. Yu, J. Jing. Deterministic generation of orbital-angular-momentum multiplexed tripartite entanglement. Phys. Rev. Lett., 124, 083605(2020).

    [42] W. Wang, K. Zhang, J. Jing. Large-scale quantum network over 66 orbital angular momentum optical modes. Phys. Rev. Lett., 125, 140501(2020).

    [43] S. Liu, Y. Lou, J. Jing. Orbital angular momentum multiplexed deterministic all-optical quantum teleportation. Nat. Commun., 11, 3875(2020).

    [44] D. Cozzolino, E. Polino, M. Valeri, G. Carvacho, D. Bacco, N. Spagnolo, L. Oxenløwe, F. Sciarrino. Air-core fiber distribution of hybrid vector vortex-polarization entangled states. Adv. Photon., 1, 046005(2019).

    [45] G. Vallone, V. D’Ambrosio, A. Sponselli, S. Slussarenko, L. Marrucci, F. Sciarrino, P. Villoresi. Free-space quantum key distribution by rotation-invariant twisted photons. Phys. Rev. Lett., 113, 060503(2014).

    [46] Y. Chen, W. Shen, Z. Li, C. Hu, Z. Yan, Z. Jiao, J. Gao, M. Cao, K. Sun, X. Jin. Underwater transmission of high-dimensional twisted photons over 55 meters. PhotoniX, 1, 5(2020).

    [47] R. Ma, W. Liu, Z. Qin, X. Su, X. Jia, J. Zhang, J. Gao. Compact sub-kilohertz low-frequency quantum light source based on four-wave mixing in cesium vapor. Opt. Lett., 43, 1243-1246(2018).

    [48] R. Simon. Peres-Horodecki separability criterion for continuous variable systems. Phys. Rev. Lett., 84, 2726-2729(2000).

    [49] I. Kogias, A. R. Lee, S. Ragy, G. Adesso. Quantification of Gaussian quantum steering. Phys. Rev. Lett., 114, 060403(2015).

    [50] P. Vaity, J. Banerji, R. P. Singh. Measuring the topological charge of an optical vortex by using a tilted convex lens. Phys. Lett. A, 377, 1154-1156(2013).

    [51] X. Deng, Y. Liu, M. Wang, X. Su, K. Peng. Sudden death and revival of Gaussian Einstein–Podolsky–Rosen steering in noisy channels. npj Quantum Inf., 7, 65(2021).

    [52] Z. Zhu, M. Janasik, A. Fyffe, D. Hay, Y. Zhou, B. Kantor, T. Winder, R. W. Boyd, G. Leuchs, Z. Shi. Compensation-free high-dimensional free-space optical communication using turbulence-resilient vector beam. Nat. Commun., 12, 1666(2021).

    Li Zeng, Rong Ma, Hong Wen, Meihong Wang, Jun Liu, Zhongzhong Qin, Xiaolong Su. Deterministic distribution of orbital angular momentum multiplexed continuous-variable entanglement and quantum steering[J]. Photonics Research, 2022, 10(3): 777
    Download Citation