Quantum discord of thermal two-photon orbital angular momentum state: mimicking teleportation to transmit an image

We formulate a density matrix to fully describe two-photon state within a thermal light source in the photon orbital angular momentum (OAM) Hilbert space. We prove the separability, i.e., zero entanglement of the thermal two-photon state. Still, we reveal the hidden quantum correlations in terms of geometric measures of discord. By mimicking the original protocol of quantum teleportation, we demonstrate that the non-zero quantum discord can be utilized to transmit a high-dimensional OAM state at the single-photon level. It is found that albeit the low fidelity of teleportation due to the inherent component of maximally mixed state, the information of all parameters that characterize the original state can still be extracted from the teleported one. Besides, we demonstrate that the multiple repetitions of the protocol, enable the transmission of a complex-amplitude light field, e.g., an optical image, regardless of being accompanied with a featureless background. We also distinguish our scheme of optical image transmission from that of ghost imaging. The formulation of two-photon orbital angular momentum state within a thermal light enables the revelation of its hidden quantumness with non-zero quantum discord and a proposal of teleporting an image.


Introduction
In a seminal paper in 1993 1 , Bennett and coworkers described the quantum protocol of teleportation to transfer any unknown quantum state from one location to another distant one, without physical traveling of the object itself. Its first experimental realization was achieved by Bouwmeester and his colleagues in 1997 using photonic polarization entanglement 2 . Since then, quantum teleportation has been playing an active role in the progress of both quantum information science and technology 3,4 . A lot of research efforts have been triggered to demonstrate numerous variants of teleportation schemes, such as entanglement swapping and quantum repeaters 5 , open-destination teleportation 6 , quantum gate teleportation and quantum computing 7,8 , port-based teleportation 9 , and teleportation of a single photon's multiple degrees of freedom 10 . The photons are an optimal choice for carrying information in the form of flying qubits and for showcasing nonlocality at a distance such that longhaul experimental teleportation was generally demonstrated with photons 11,12 . Recent progress has been made to realize the teleportation from a transmitter on Earth to a receiver on a satellite, towards a global scale 13 .
In the original scheme of quantum teleportation 1 , two communicating parties, Alice and Bob, must initially share an entangled pair of quantum particles that serves as the quantum channel. For two-level quantum systems, e.g., polarization qubits as in the 1997 experiment 2 , the entangled state generated by parametric down-conversion can be chosen as, Ψ j i ðÀÞ , where H and V represent the horizontal and vertical polarizations, respectively. In the meantime, Alice is given another third photon C whose polarization state ψ j i C will be teleported from Alice to Bob. Alice performs a joint Bell measurement on her photons A and C and obtains one possible outcome of four Bell states, . Then Alice communicates her measurement outcome to Bob via a classical channel. Accordingly, Bob applies a unitary operation, namely the identity operator I or Pauli operators σ x , σ y , σ z on his photon B and recovers the original state. Quantum teleportation has been already generalized to high-dimensional discrete-variable 14,15 or infinite-dimensional continuous-variable systems 16,17 . Besides, it was demonstrated that some mixed states that contain noisy entanglement, rather than a pure entanglement, can also be exploited as quantum channels for teleportation 18,19 . Moreover, teleportation using only classical channels has been explored and quantified in terms of fidelity bounds 20 . Here, by mimicking the standard protocol of quantum teleportation, we propose a new scheme of using a pseudo-thermal light source to realize the transmission of a high-dimensional orbital angular momentum (OAM) state and even an optical image with quantum discord. Our scheme is based on the formulation of a density matrix on the basis of OAM eigenstates to fully describe two-photon state, which reveals the inherently high-dimensional spatial correlation of non-zero quantum discord. At the singlephoton level, the fidelity of the transferred sate is quantified, and for multiple repetitions of our protocol, the distinction of our scheme for optical image transmission is also discussed from ghost imaging.

Results
Thermal two-photon state in the OAM space Our theory is based on the adoption of OAM eigenstates 21,22 , in terms of Laguerre−Gaussian (LG) modes, to formulate the density matrix describing two-photon states within a thermal light source. The quantum theory of the second-and higher-order coherence of light was introduced in 1963 by Glauber 23,24 . Here, we adopt ref. 25 to model the thermal radiation as a collection of independent atoms emitting radiation randomly, and its density operator is written as, ρ 0 ¼ P fn k g P fn k g fn k g j i fn k g h j, where P fn k g ¼ Q k Pðn k Þ with P(n k ) being the probability for n k photons in the mode k, and the symbol {n k } denotes a set of numbers n k1 , n k2 , n k3 ,…, etc, of photons excited in very mode. For twophoton case, n ¼ P k n k ¼ 2, then we can specify denoting only one excited photon in each mode k j i and k 0 j i, respectively. As the LG modes form an orthogonal and complete computation basis 26 , the light field of any mode k j i can be expressed in terms of LG modes '; p j i, where ' and p are the azimuthal and radial mode indices, respectively. Accordingly, we have which indicates that the chaotic source, e.g., thermal light, can be thought of as incoherent statistical mixtures of photon pairs, '; p j i ' 0 ; p 0 j i. In our proposed scheme sketched in Fig. 1, the light derived from a thermal source is directed to a non-polarizing beam splitter (BS) and then is divided into two paths. The interaction with BS changes the two-photon state to, Subsequently, we use coincident measurement between photons in two paths, which means that it is just this joint detection that post selects the following two-photon state, where Ψð'; p; . It is found that Eq. (4) can be equivalently written as the following density matrix (see Supplementary Information for details), with subscript, a and b denote the photon a and b in two paths, respectively. We show the physically intuitive quantum−mechanical structure of Eq. (5), which describes the two-photon state within a thermal light source: ρ 0 ab of Eq. (6) is merely a diagonal fully separable state, corresponding to the maximally mixed state. In contrast, Ψ j i ab of Eq. (7) is inherently a pure entangled state, which behaves in a very analogous way to the OAM entanglement 27 .
To further derive the specific expression of P ';p , we consider the chaotic thermal light source's crossspectrum density function. Because of its generality and validity, the Gaussian−Schell model was extensively employed, which can be written as 28 , where G 0 is a constant, σ S is the transverse size, and σ g is the transverse coherence width of the source. In biphoton case, two-photon amplitude, Φðr 1 ; r 2 Þ, can be expressed as the single-sum expansion based on the Schmidt decomposition, Φðr 1 ; r 2 Þ ¼ P n ffiffiffiffi ffi λ n p u n ðr 1 Þv n ðr 2 Þ, where u n (r 1 ) and v n (r 2 ) are Schmidt modes for signal and idler photons, respectively, and λ n is the eigenvalue 29 . We perform a similar mode decomposition on W(r 1 , r 2 ) in the full set of normalized LG modes, i.e., W(r 1 , r 2 ) = ∑ ℓ,p,ℓ′,p ′ ƒ ℓ , ℓ′,p,p′ LG ℓ,p (r 1 )LG ℓ′,p′ (r 2 ). After some algebra (see Supplementary Information for details), we can find that, where tan β = 2σ S /σ g . The main incentive of our present work is to explore the thermal two-photon OAM state of Eq. (5) for the possibility of mimicking teleportation to realize the transmission of a high-dimensional OAM state and even a two-dimensional complex-amplitude optical image, regardless of no prior entanglement.

Proof of the separability of thermal two-photon OAM state
Based on the concept of robustness of entanglement, we first prove the separability of two-photon state of Eq. (5). Following the procedure presented in ref. 30 , we know the robustness of entanglement of Eq.
which is a separable state by construction. In other words, all entanglement in ρ Q can be washed out after being mixed with ρ À S , and the required minimal amount of ρ À S is R. Thus the mixture, separable. Besides, by considering the fact that ρ À S ¼ 1 Eq. (5) as, which is obviously separable. Thus we have proved that the thermal two-photon source is not entangled per se. However, our theoretical formulation also provides an intuitive picture to reveal the quantum correlations hidden in thermal two-photon OAM state.
The revelation of the quantumness of thermal two-photon OAM state When quantum correlation is mentioned, we might be immediately tempted to think of entanglement. But it deserves our attention that entanglement is not the only aspect of quantum correlations, namely, some quantum systems even without entanglement are capable of possessing correlations that cannot be simulated by classical physics 31 . In the OAM Hilbert space, we have formulated the density operator of Eq. (5) to fully describe the thermal two-photon source. It is crucial to distinguish quantum correlations beyond entanglement. In parallel with the entanglement measure used for entanglement vs separability criteria, discord was proposed as a figure of merit for characterizing the nonclassical resource in the quantum vs classical scenario 32,33 It is defined as the discrepancy between quantum mutual information (total correlation) and classical correlation in a bipartite system. However, the evaluation of quantum discord requires considerable numerical minimization. From a geometric perspective, a distance-based notion of discord 34 , was introduced instead. For an arbitrary state ρ, this geometric discord is defined as D ρ ð Þ ¼ min denotes the set of zero-discord states, and ρ À χ k k 2 ¼ tr ρ À χ ð Þ 2 is the square norm in the Hilbert-Schmidt space. Progress in computing D(ρ) has been made by showing that 35 where the minimization is taken over all local bases η B k on photons in path B. After a lengthy algebra based on Eq. (11), we can obtain the analytic geometric discord for Eq. (5) as We have plotted D(ρ) in Fig. 2 in different OAM subspaces, whose dimension is specified by with ' ranging from −L to L and p from 0 to P. One can see that, regardless of zero entanglement, the quantum discord is surprisingly non-zero, and therefore revealing the quantumness of correlations hidden in the thermal two-photon state. We also plot the geometric discord D(ρ Q ) for the pure entangled state of Eq. (7) within the same subspaces. Two extreme cases, σ g = ∞ and σ g = 0, deserve our special attention. For σ g = ∞, the source is completely coherent, and none of the quantum correlations can be extracted from either the thermal state or the entangled state, as DðρÞ ¼ Dðρ Q Þ ¼ 0. For L = P = ∞, we can also obtain, D 1 ρ ð Þ ¼ 1= σ g =σ s þ 2σ s =σ g À Á 4 , which is shown by the purple curves in Fig. 2. The other extreme is the case of the completely incoherent light source, where ρ Q is maximally entangled and P ';p becomes a constant. For comparison, if we also apply the geometric discord to the initial state of Eq. (1), a straightforward analysis finds D(ρ 0 ) = 0. Thus we are allowed to conclude that it is just the function of BS and the post-selection of joint detection that extract the nonclassical correlations. As we will demonstrate below, we can still utilize such a non-entangled yet nonclassical thermal two-photon state to realize the transmission of a high-dimensional OAM state at the singlephoton level and that of an optical image with multiple repetitions of our protocol.
Mimicking teleportation to transmit the high-dimensional OAM state of a single photon As shown in Fig. 1, we use an additional mirror (M) in the path of photon b to flip the OAM value. For simplicity and without loss of generality, we consider a completely incoherent source with σ g = 0, which leads to tan β 2 ¼ 1 and indicates that P ';p becomes a constant, as it is independent of ' or p. If we further consider a subspace with the radial index p ranging from 0 to P and the azimuthal index ' from −L to +L, Besides, for convenience, we reindex all the eigenstates of LG modes by assigning j As a result, we can rewrite Eqs. (6) and (7) as, Here we explore the above quantum correlations hidden in the thermal two-photon state to serve as a quantum bridge to accomplish the transmission of a high-dimensional OAM state for a single photon, akin to the standard teleportation protocol. In Fig. 1, the photons a and b in two paths are spatially separated. Assume Alice holds photon a while Bob possesses photon b. To mimic the teleportation, Alice has a third photon c at hand, which is encoded with the highdimensional OAM state, mathematically corresponding to the complex light field of an optical image, ψ(r). This is realized by passing the photon c in an initial fundamental Gaussian mode through an encoded holographic mask. According to the concept of digital spiral imaging 26,36 , we can rewrite the complex light field in terms of LG modes as, ψ LG ';p ðrÞ Â Ã Ã ψðrÞdr. In our proposed experiment, we use a series of imaging lenses to image the plane of a chaotic thermal light source to the plane of an ICCD camera (via 4f 1 for photon b) and to the plane of a high-dimensional Bell state measurement device (d-BSM) (via 4f 2 for photon a), respectively. Meanwhile, another single imaging lens (4f 3 ) is used to image the mask onto the d-BSM as well in the path of photon c. Then the wavefunction of the whole system can be written as, Here the d-BSM realizes the function of projecting photons a and c onto the generalized high-dimensional Bell states in the Hilbert space spanned by the LG modes. Then, we are allowed to rewrite the whole state ρ abc , in terms of the generalized Bell states of photons a and c, as where ψ m;n j i ac ¼ 1 ffiffi d p P j e i2πjm=d j j i a j þ n ð Þmode d j i c are two-photon generalized Bell states, and φ j being a unitary matrix. Equation (15) forms the key result of the present work, which suggests the possibility of using a thermal light source to mimic the teleportation protocol for realizing the transmission of a high-dimensional OAM state of a single photon and even a complex-amplitude optical image. After performing the generalized Bell state measurement on photons a and c, Alice gets one possible result out of d 2 outcomes, k = (m,n), where m; n ¼ 0; 1; 2; :::; d À 1. Then Alice communicates classically her measurement result k = (mn) with Bob. Accordingly, Bob applies the conditional unitary matrix U mn to the photon b, namely This suggests that Bob can retrieve an output state of photon b, which is just a mixture of an exact replica of the

Discussion
Our protocol at the single-photon level We assume that the state encoded by photon c is represented by a high-dimensional vector, jψi c ¼ P j A j j j i c ¼ P ';p A ';p '; p j i c , which can be equivalently characterized by two spectra, a ';p ¼ A ';p and ϕ ';p ¼ argðA ';p Þ, respectively. It is noted that in classical optics, such a high-dimensional state just corresponds to a complex-amplitude optical image, namely, ψ r ð Þ ¼ P ';p A ';p LG ';p r ð Þ. The calculation of a ';p and ϕ ';p for a Clover image is shown in the section of "Materials and methods". From Eq. (15), without loss of generality, assume the d-BSM results at Alice's side are k ¼ ðm; nÞ ¼ ð0; 0Þ, (2,0) and (10,10), namely, ψ 0;0 j i ac , ψ 2;0 j i ac , and ψ 10;10 j i ac , respectively. Accordingly, photon b will be projected onto the state of U þ 0;0 ψ j i b , U þ 2;0 ψ j i b , and U þ 10;10 ψ j i b . For an easy visualization, we illustrate in Fig. 3 the spectra a ';p and ϕ ';p of these states, in which the trivial background of a maximally mixed state 1 d P j j j i b j h j is omitted. In Fig. 3a, b, as U þ 00 ¼ I, photon b just acquires the original state of photon c, accompanying with the trivial background. In contrast, conditional on ψ 2;0 j i ac and ψ 10;10 j i ac and without performing the desired unitary operations, we can only obtain the incorrect highdimensional state, U þ 2;0 ψ j i b , and U þ 10;10 ψ j i b , respectively. According to U mn ¼ P k e i2πkm=d k j i ðk þ nÞ mode d h j , we have U þ mn ¼ P k e Ài2πkm=d ðk þ nÞ mode d j ik h j. Therefore, we have U þ 2;0 ¼ P k e Ài4πk=d k j i k h j such that the spectrum of a ';p remains unchanged while each phase ϕ ';p is individually shifted by an amount of Δϕ ¼ À4πk=d, see Fig. 3c, d. While for U þ 10;10 ¼ P k e Ài20πk=d k þ 10 j i k h j, as shown by Fig. 3e, f, we can see that both spectra of a ';p and ϕ ';p are changed. In other words, the two latter cases without desired unitary operations can only lead to the incorrect transferred states.
As well known, quantum teleportation is distinct in that it allows a single copy of a state to be teleported with unity fidelity. Therefore, it is crucial to investigate the fidelity of the transferred state in our protocol, which can be written as, The particular case of d = 2 leads to F = 3/4, which agrees to the optimal fidelity achievable with a two spin-1/ 2 Werner stat 37 . However, high dimensions are essential for encoding an optical image. From Eq. (17), we find that the fidelity decreases inversely as the dimensionality increases, approaching to the limit of F = 1/2 for d ! 1, due to the inherent existence of the background maximally mixed state of 1 d P j j j i j h j in Eq. (16). In an actual experiment, the background can only be subtracted by sending multiple copies of the same state, which changes the resources used in the scheme since Alice performing the Bell-state measurement needs to communicate d bits to Bob for each repetition of the protocol, corresponding to one of the d possible outcomes of the Bell-state measurement. This classical communication itself may exceed the number of bits used to describe the image. However, if the density matrix of Eq. (16) could be reconstructed, e.g., via quantum state tomography, one can see that the information about all the parameters of the original state, a ';p and ϕ ';p , could still be fully known, Albeit a low fidelity of teleportation.

Our protocol with multiple repetitions
Reconstruction of an optical image requires multiple single-photon events. Thus it is also interesting to consider the multiple repetitions of our protocol. As mentioned above, the high-dimensional state just corresponds to a complex-amplitude optical image of a Clover, see the inset of Fig. 1. If we post-select the generalized d-BSM result of ψ 0;0 j i ac , we can record the single photons to reconstruct the Clover image with the ICCD camera. Then we can plot the intensity and phase profiles of the Clover image in Fig. 4a, b, respectively, which corresponds to the high-dimensional state characterized by Fig. 3a, b. Note that the background is trivially omitted. It can be seen that the reconstruction of the pure Clover image is ensured by the second term of Eqs. (15) or (16), which when viewed separately is due to the spatial correlations established by the entangled state of Eq. (7) or Eq. (14). This is just physically equivalent to using a quantum source of biphoton pure entangled state. In contrast, when the ICCD camera is triggered by another d-BSM results, e.g, ψ 2;0 j i ac and ψ 10;10 j i ac , we could only record a fuzzy image, see Fig. 4c-f, respectively, as a result of the incorrect transmissions of the high-dimensional states, see Fig. 3c-f, respectively. In principle, one should perform a desired unitary operations for each repetition in order to recover the correct image.

Comparison with ghost imaging
As it is known, another well-established application of thermal light source, is the technique of ghost imaging. The ghost imaging could be traced back to the experiment conducted in 1995 that exploited the nonlocal correlations of entangled photon pairs created by parametric down-conversion 38 . However, it was demonstrated that a coincidence image can also be obtained with thermal light source [39][40][41][42] , thus raising a question whether entanglement was truly necessary for ghost imaging [43][44][45] . Several efforts have contributed to providing some insightful understanding on this historic argument [46][47][48] . Here our scheme uses thermal light, and the recorded images in both two cases are embedded in a featureless background. However, we distinguish our scheme of image teleportation from that of ghost imaging. In ghost imaging, only a pair of photons are employed; one photon in the pair illuminates the object and is collected by a bucket detector while the other photon is utilized to record the images with a spatially resolved detector, e.g., an intensive CCD camera. In contrast, here we employ three photons; two photons (a and b) in a pair derived from the thermal light source serve as the channel for teleportation, while the third photon (c) is directed to illuminate the object and acquires the image information to be teleported. We then follow the standard procedure of teleportation to perform the generalized Bell state measurements (d-BSM) on photons a and c, and use photon b to record the teleported images with the ICCD camera nonlocally. Besides, it is shown that the correct image could be obtained only after the desired unitary transformation, conditional on the d-BSM results of photons a and c.
Besides, the formulation of thermal two-photon OAM state, ρ ¼ ρ Q þ ρ 0 ab , offers an intuitive picture of twophoton quantum interference. If we look ρ Q separately, it just accounts for the existence two-photon OAM interference. In contrast, if we look ρ wholly, the interference of ρ Q has been already canceled out by ρ 0 ab , thus leading to the inexistence of two-photon interference. In this regard, our work can also shed new light on the historical debate as to whether the origin of ghost imaging is quantummechanical or classical.
In summary, we have formulated a density matrix on the basis of OAM eigenmodes to fully describe the thermal twophoton state and reveal the quantumness with non-zero discord. We also demonstrate that albeit a low teleportation fidelity, such a thermal two-photon state can be explored for realizing the transmission of a high-dimensional OAM state at the single-photon level and that of an optical image with multiple repetitions of our protocol. This is obviously distinguished from the traditional method of using a prior entangled source. However, as in the case of using entanglement, the practical implementation of our scheme is hindered by the high-dimensional Bell state measurement (d-BSM). For our demonstration, it requires the manipulation of two-photon d-level states encoded by LG modes with both azimuthal and radial indices, which remains more challenging with the current experimental technique. In the future, our work can also call for further studies of using thermal multi-photon state to demonstrate some new quantum information tasks, such as remote state preparation 49 and novel imaging with undetected photons 50 .

Materials and methods
The calculation of both intensity and phase spiral spectra for the Clover image In Fig. 1, the Clover object together with the single-mode fiber and the single-photon detector (not shown) inversely defines the state of photon c, ψ j i c ¼ P ';p A ';p '; p j i c , where A ';p ¼ R LG ';p ðrÞ À Á Ã ψðrÞdr denotes the overlap probability amplitude and is related to digital spiral imaging of the Clover object. Accordingly, we present in Fig. 3a, b the intensity spectrum, a ';p , and the phase spectrum, ϕ ';p . Due to the three-fold rotational symmetry of the Clover considered here, one can see that only OAM components of ' ¼ 0; ± 3; ± 6, dominate in the spiral spectrum. Besides, as the adopted beam waist of the LG modes w = 0.4469 mm is comparable to the radius of the Clover object R = 1 mm, only those lower-order LG modes with a small p index contribute together to form the Clover object.