# Experimental detection of polarization-frequency quantum correlations in a photonic quantum channel by local operations

###### Abstract

The measurement of correlations between different degrees of freedom is an important, but in general extremely difficult task in many applications of quantum mechanics. Here, we report an all-optical experimental detection and quantification of quantum correlations between the polarization and the frequency degrees of freedom of single photons by means of local operations acting only on the polarization degree of freedom. These operations only require experimental control over an easily accessible two-dimensional subsystem, despite handling strongly mixed quantum states comprised of a continuum of orthogonal frequency states. Our experiment thus represents a photonic realization of a scheme for the local detection of quantum correlations in a truly infinite-dimensional continuous-variable system, which excludes an efficient finite-dimensional truncation.

## I Introduction

Photons represent natural carriers of quantum information and play a crucial role in quantum communication and quantum computation Zoller (). As flying qubits they provide the messengers among different nodes in quantum communication protocols. During transmission, however, a photonic qubit will unavoidably couple to an environment, i.e., it will become an open quantum system. On the one hand, correlations between the qubit state and its environment can significantly influence the qubit’s time evolution ringbauer2015 () and thereby lead to detrimental effects on the capacity of the quantum channel. System-environment correlations are responsible for non-Markovian effects TheWork (); lofranco2013 (); breuer2015 (), and, when present in the initial state, can even lead to a complete breakdown of any standard master equation description Gorini (); Pechukas (); Lindblad (); Royer (); Buzek (). On the other hand, quantum correlations form the foundation for quantum information technologies. Quantum entanglement, for instance, provides the quantum advantage for quantum key distribution ekert1991 () and quantum teleportation bennett1993 (), while quantum discord modi2012 () plays an important role for entanglement activation piani2011 (); streltsov2011 (); Adesso2014 (); Orieux2015 () and distribution streltsov2012 (), as well as for quantum metrology modi2011 (); girolami2013 (). It is therefore of both fundamental and practical importance to develop efficient ways to characterize quantum correlations in photonic systems.

In optical experiments, quantum correlations of small systems, e.g., a photon pair, can be detected by performing tomography on the total system xujs2010 (). However, for high-dimensional systems, resource constraints and lack of control soon render such measurements impossible. A particularly complex, yet common situation is encountered when the photon’s polarization state interacts with its own frequency modes, for instance in a birefringent crystal licf2011 (). In this case, correlations emerge between the polarization state, which encodes the qubit, and a continuum of frequency states. This type of quantum channel is frequently used in quantum information applications even though the environmental state describing the frequency modes is not fully controllable and hard to determine.

In contrast, the -dimensional subsystem, i.e., the polarization qubit, can be fully controlled and detected with local operations even in such a complex system. The question thus arises, does this allow us to extract information about the quantum correlation of the entire system? One strategy to answer this question is to use the time evolution of reduced states of the -dimensional subsystem. Detecting, for instance, an increase of the trace distance between pairs of states in the subsystem constitutes a witness for total correlations in one of the initial states laine2010 (). This scheme was applied experimentally in optical architectures licf2011 (); smirne2011 (). Further theoretical extensions allow for the detection of quantum correlations of an unknown system-environment state using only local operations on the system gessner2011 (), as demonstrated in a trapped-ion experiment gessner2014 (). The method can also be used to discriminate between quantum and classical correlations Cialdi2014 () and, quite recently, an application to the local detection of a quantum phase transition has been developed Gessner2014epl ().

In this letter we focus on the open quantum system formed by the polarization degree of freedom of a single photon which is coupled to its environmental frequency degrees of freedom in a birefringent quantum channel licf2011 (). Single photons are generated from a self-assembled quantum dot (QD) and initial system-environment correlations are created by sending the photon through a calcite crystal. We show that it is possible to experimentally detect and quantify the quantum correlations in this system by carrying out only local operations on the polarization degree of freedom of the photon.

## Ii Theoretical framework

First, we briefly recall the theoretical framework of the local detection method, introduced in Ref. gessner2011 (). Given a photon state , where “” represents the polarization subsystem and “” represents the frequency modes, how can we measure the quantum correlation of this state? One method is to map this state to an associated state , which contains only classical correlations, and then calculate the distance between and . This map can be realized through a controlled local dephasing process, i.e.,

(1) |

where and the are the eigenstates of the reduced state . The quantum discord of is zero gessner2011 (). Therefore, the quantum correlations in can be quantified as luo2008 (); bellomo2012 (); aaronson2013 ()

(2) |

where denotes the trace norm. We note that the trace distance (2) can be interpreted in terms of the probability for a successful discrimination of the two quantum states by a single measurement Hayashi (). Suppose now that both and evolve in time through the same unitary evolution operator leading to the reduced open system state and , respectively. Since the partial trace over the environment is a positive operation, the local trace distance

(3) |

represents a lower bound for the quantum correlations in the state ruskai1994 (), i.e. we have

(4) |

As this relation holds for all we find

(5) |

According to this general inequality the quantity on the left-hand side – a local quantity accessible through only local operations on the open system – represents a witness for the quantum correlations in the total system-environment state gessner2013 ().

## Iii Experimental framework and results

In our experiment, we suppose Alice prepares the initial state and sends it to Bob. The state is unknown to Bob, who now can perform local operations on the polarization state to determine the quantum correlations in the initial state. The experimental setup is shown in Fig. 1.

On Alice’s side, the job is divided into two parts. The first part is the generation of single photons (not shown). The photons are emitted from a single InAs/GaAs QD, which is isolated by a -m-diameter pinhole in a gold mask. The exciting power intensity is about , which makes the QD oversaturated. This power intensity, together with the gold mask on the surface of the QD sample, induces a significant broadening of the QD emission line because of noise effects, such as spectral diffusion sallen2010 (); wolters2013 (). The frequency state then becomes a mixed state because the inhomogeneous broadening is far larger than the homogeneous broadening. The resulting emission line is then filtered by a Fabry-Perot (FP) cavity with a finesse of 28. The line shape of the single photons is measured by a Michelson interferometer and the result is shown in Fig. 2(a), where is the delay between the two arms of the interferometer. Assuming a Lorentzian line shape

(6) |

with a center frequency corresponding to nm, we can extract the line width ps, by fitting the measured data to with .

The second part of Alice’s job is the preparation of the initial state . A polarizer prepares the photon into the pre-initial state , where , are the horizontal and vertical polarization states, respectively. Then, the photon moves through a birefringent quantum channel, generated by a calcite of length with its optical axis along the vertical direction. Finally a half-wave plate (HWP0) driven by a random number generator (RNG) rotates the basis of the state to a random direction in order to make it impossible for Bob to know the local eigenbasis on the side of Alice. At this point, Alice’s preparation of the initial state is completed and the state is passed on to Bob.

Knowing the parameters used for state preparation, Alice can theoretically determine the quantum correlations in the state based on Eq. (2), which in the described situation leads to the following expression (see Supplement 1),

(7) |

with corresponding to the length of the calcite. To measure the birefringence of the calcite we use again the Michelson interferometer. The result for a calcite of length mm is shown in Fig. 2(b). From the shift of the zero delay, we obtain the birefringence , which coincides with the parameter provided by the calcite’s manufacturer.

The setup on Bob’s side consists of four sections. The first section is made up of three removable mirrors (RM1, RM2 and RM3). By placing and removing these mirrors, the photons are led to different functional modules in the experimental setup. One of them is the second section, which implements the controlled dephasing map, Eq. (1), using a computer-controlled half-wave plate (HWP1) and a very long polarization-maintaining (PM) optical fiber ( m). The third section is a Michelson delay setup, which represents the unitary time evolution. Here, the angle of the half-wave plate (HWP2) allows us to generate a set of time-evolution operators . This section further includes a polarizing beam splitter (PBS1), two -placed quarter-wave plates (QWP1 and QWP2), a fixed mirror (M1), and a delayable mirror (M2), whose coordinate corresponds to the time variable. The last section is the tomography setup (T1 and T2). T1 includes QWP3, HWP3, PBS2 and two avalanche photodiodes (APD1 and APD2). The setup of T2 is the same as T1.

The measurement of the quantum correlations occurs in two steps. The first step is the preparation of . With RM1 in place, the photons are directed to T2, where the tomography of the reduced state is performed on the polarization subsystem. Then the computer calculates the eigenstates and of , and controls HWP1 to rotate to the angle . At this point, we remove RM1 and place RM2, causing the photons to go through the controlled dephasing map. The birefringence of the PM fiber is , so we have . This means that the PM fiber is long enough to completely remove coherences between the orthogonal subspaces spanned by and , which are responsible for the quantum correlations in terms of the quantum discord. This way the state is prepared at the end of the PM fiber.

In the second step, we first remove all of the removable mirrors and send the photons with state to go through the unitary evolution of . Considering only the polarization state, we obtain the evolved reduced state . This state can be measured by the tomography setup in T1. We then place RM2 and RM3. Now, instead of , we subject to the same evolution, and by local tomography we obtain another reduced state . We can then calculate the trace distance between the two reduced states. Here, we choose the parameters , as and , with and . Taking mm as an example, we obtain the results shown in Fig. 3. The theoretical analysis yields the following expression for the local trace distance (see Supplement 1),

(8) |

Since this expression is independent of , as clearly confirmed by the experimental data (see Fig. 3), we write for simplicity. Note that according to expressions (7) and (8) the general inequality (4) is, in the present special case, a simple consequence of the triangular inequality. For the Lorentzian spectrum one finds

(9) |

which is the theoretical prediction shown in Fig. 3.

Now Alice changes the length of the calcite crystal, thus creating a set of different initial states which correspond to the photon states at different positions in a quantum channel. The crystal length is chosen as mm with . For each crystal length, Alice knows the total state and determines the theoretical quantum correlation employing Eq. (7). The results are shown in Fig. 4 by the blue line. Bob receives the initial states, and uses local operations to experimentally detect the quantum correlations. Using Eq. (9) we find that the theoretical prediction for Bob’s measurement results is given by

(10) |

which is shown by the red line in Fig. 4, and compared to the measurement result shown as black dots. This experimental result is derived by fitting the experimental data of (Fig. 3 shows an example of these data) using Eq. (9), and the uncertainty in the result is determined by the fittings. As is evident from the figure, the local detection scheme of Bob provides a good lower bound for the quantum correlations in the initial state prepared by Alice: The bound is rather close to the actual value of the quantum correlations and, moreover, exhibits the same monotonicity as the initial correlations.

## Iv Conclusions

To conclude, we have prepared quantum correlated initial states in a photonic system composed of a single photon’s polarization state and its frequency modes using a birefringent quantum channel. Our all-optical experiment demonstrates that these quantum correlations can be detected and quantified by performing only local operations on the polarization subsystem. The theoretical analysis shows that the local detection scheme leads to a close lower bound for the quantum correlations of the total initial state. The ability to detect and estimate these bipartite correlations with access to only one of the subsystems distinguishes our work from previously developed methods Girolami2012 (), which had been realized experimentally with two-qubit systems based on nuclear magnetic resonance Silva2013 (); Paula2013 (). The methods presented in Refs. Girolami2012 (); Silva2013 () allow for estimation of the correlations of unknown states by local operations on both subsystems if one of the subsystems is a qubit, and both subsystems are finite-dimensional and experimentally under control. In contrast, the local detection method employed here only requires control over one of the two systems (of arbitrary dimensions) and the presence of a suitable interaction between them.

Finally, we emphasize that our work represents the first experimental realization of this scheme for mixed correlated initial states in an infinite-dimensional tensor product Hilbert space, where a finite-dimensional truncation is prevented by the contributions of a truly infinite number of orthogonal continuous-variable states. The developed experimental techniques can be naturally extended to the study of correlations in other photonic open systems and quantum channels for which only a local part can be controlled and measured.

## Funding Information

This work is supported by the National Natural Science Foundation of China (Grants Nos. 61490711, 11274289, 11304305, 11325419, 61327901 and 91321313), the National Basic Research Program of China (Grants Nos. 2011CB921200 and 2013CB933304), the Strategic Priority Research Program(B) of the Chinese Academy of Sciences (Grant No. XDB01030300), the Fundamental Research Funds for the Central Universities (No. WK2470000011), the China Postdoctoral Science Foundation funded project (Grant No. 2012M521229), and the German Academic Exchange Service (DAAD). M.G. thanks the German National Academic Foundation for support. C.-F.L. and H.-P.B. acknowledge support from the European Union (EU) through the Collaborative Project QuProCS (Grant Agreement No. 641277).

See Supplement 1 for supporting content.

## References

- (1) L.-M. Duan, M. D. Lukin, J. I. Cirac and P. Zoller, Nature 414, 413 (2001).
- (2) M. Ringbauer et al., Phys. Rev. Lett. 114, 090402 (2015).
- (3) H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press, Oxford, 2007).
- (4) R. Lo Franco, B. Bellomo, S. Maniscalco, and G. Compagno, Int. J. Mod. Phys. B 27, 1345053 (2013).
- (5) H.-P. Breuer, E.-M. Laine, J. Piilo, and B. Vacchini, arXiv:1505.01385 (2015).
- (6) V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, J. Math. Phys. (N.Y.) 17, 821 (1976).
- (7) P. Pechukas, Phys. Rev. Lett. 73, 1060(1994).
- (8) G. Lindblad, J. Phys. A 29, 4197 (1996).
- (9) A. Royer, Phys. Rev. Lett. 77, 3272 (1996).
- (10) P. Štelmachovič and V. Bužek, Phys. Rev. A 64, 062106 (2001); Phys. Rev. A 67, 029902(E) (2003).
- (11) A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991).
- (12) C. H. Bennett et al., Phys. Rev. Lett. 70, 1895 (1993).
- (13) K. Modi, A. Brodutch, H. Cable, T. Paterek, and V. Vedral, Rev. Mod. Phys. 84, 1655 (2012).
- (14) M. Piani et al., Phys. Rev. Lett. 106, 220403 (2011).
- (15) A. Streltsov, H. Kampermann, and D. Bruß, Phys. Rev. Lett. 106, 160401 (2011).
- (16) G. Adesso, V. D’Ambrosio, E. Nagali, M. Piani, and F. Sciarrino, Phys. Rev. Lett. 112, 140501 (2014).
- (17) A. Orieux et al., Phys. Rev. Lett. 115, 160503 (2015).
- (18) A. Streltsov, H. Kampermann, and D. Bruß, Phys. Rev. Lett. 108, 250501 (2012).
- (19) K. Modi, H. Cable, M. Williamson, and V. Vedral, Phys. Rev. X 1, 021022 (2011).
- (20) D. Girolami, T. Tufarelli, and G. Adesso, Phys. Rev. Lett. 110, 240402 (2013).
- (21) J.-S. Xu et al., Nat. Commun. 1, 7 (2010).
- (22) C.-F. Li, J.-S. Tang, Y.-L. Li, and G.-C. Guo, Phys. Rev. A 83, 064102 (2011).
- (23) E.-M. Laine, J. Piilo, and H.-P. Breuer, Europhys. Lett. 92, 60010 (2010).
- (24) A. Smirne, D. Brivio, S. Cialdi, B. Vacchini, and M. G. A. Paris, Phys. Rev. A 84, 032112 (2011).
- (25) M. Gessner and H.-P. Breuer, Phys. Rev. Lett. 107, 180402 (2011).
- (26) M. Gessner et al., Nature Phys. 10, 105 (2014).
- (27) S. Cialdi, A. Smirne, M. G. A. Paris, S. Olivares, and B. Vacchini, Phys. Rev. A 90, 050301(R) (2014).
- (28) M. Gessner, M. Ramm, H. Häffner, A. Buchleitner, and H.-P. Breuer, Europhysics Letters 107, 40005 (2014).
- (29) S. Luo, Phys. Rev. A 77, 022301 (2008).
- (30) B. Bellomo, R. Lo Franco, and G. Compagno, Phys. Rev. A 86, 012312 (2012).
- (31) B. Aaronson, R. Lo Franco, G. Compagno, and G. Adesso, New J. Phys. 15, 093022 (2013).
- (32) M. Hayashi, Quantum Information (Springer, Berlin, 2006).
- (33) M. B. Ruskai, Rev. Math. Phys. 6, 1147 (1994).
- (34) M. Gessner and H.-P. Breuer, Phys. Rev. A 87, 042107 (2013).
- (35) G. Sallen et al., Nature Photonics 4, 696 (2010).
- (36) J. Wolters, N. Sadzak, A. W. Schell, T. Schröder, and O. Benson, Phys. Rev. Lett. 110, 027401 (2013).
- (37) D. Girolami, and G. Adesso, Phys. Rev. Lett. 108, 150403 (2012).
- (38) I. A. Silva et al., Phys. Rev. Lett. 110, 140501 (2013).
- (39) F. M. Paula et al., Phys. Rev. Lett. 111, 250401 (2013).

SUPPLEMENTARY MATERIAL

### .1 Preparation of the initial state

#### .1.1 The pre-initial state

The pre-initial state is a product state, , where the density matrix of the two-level system represented in the basis reads

(S1) |

and the environmental state is a mixed state given by a sum over closely spaced frequencies,

(S2) |

In order to avoid problems with non-normalizable continuum states, we formally represent the state by a discrete sum of frequencies, and perform the continuum limit () after obtaining the relevant quantities, such as the trace distance. The frequency states form a complete orthonormal set, i.e., and . The conditions and then ensure that indeed describes a positive and normalized quantum state. In the following we assume for simplicity that is a unimodal function with a width of order which is symmetric about the center frequency .

#### .1.2 The correlated initial state

In order to prepare a quantum state with polarization-frequency correlations, we send the pre-initial state

(S3) |

through a birefringent crystal, which generates an interaction of the form

(S4) |

where corresponds to the time spent inside the crystal. The state evolves into the initial state

(S5) |

Its reduced state for the two-level system is given by

(S6) |

with and denotes the center frequency of the distribution . For fixed , the initial phase can be chosen such that , i.e., . Since is symmetric about , we obtain a real-valued function

(S7) |

For the implementation of the local detection protocol, the first step consists in obtaining the eigenbasis of . Diagonalization of this matrix yields the eigenvalues,

(S8) |

and corresponding eigenvectors,

(S9) |

Before passing the state on to Bob, Alice implements a random local unitary rotation to hide her local eigenbasis. This is implemented by a half-wave plate which is controlled by a random number generator. As a consequence, the eigenstates are rotated to and . This local unitary operation has no effect on the correlation properties of the state. Redefining the basis vectors and accordingly always allows us to map the state back onto the original non-rotated state. Thus, for the following theoretical analysis, we can restrict to the case , i.e., effectively disregard this additional local rotation.

### .2 The local dephasing operation

The local dephasing must be carried out in the eigenbasis of the reduced state (S6). Hence, the desired operation is

(S10) |

Local application to the state yields

(S11) | ||||

where we have used Eq. (.1.2).

#### .2.1 Experimental implementation

We implement the controlled dephasing operation by sending the quantum state through a very long polarization-maintaining fiber after rotating the eigenstates and onto the principal axes of the fiber with a computer-controlled half-wave plate. Formally, the generated evolution is then described by

(S12) |

where corresponds to the time the photon spends inside the fiber. This can be reformulated as

(S13) |

with

(S14) |

We thus obtain

(S15) | ||||

This state realizes the desired dephased state to a good approximation. To see this, consider for instance,

(S16) |

Due to the condition , the function oscillates very rapidly compared to variations of and the remaining exponential functions. Such rapid oscillations cannot be resolved and, hence, must be replaced by its mean value over one period (in frequency space):

(S17) |

We obtain correspondingly

(S18) |

With these replacements reduces to .

#### .2.2 Reduced states

The reduced states of and coincide by construction, i.e., and . This can be confirmed explicitly as follows. First note that

(S19) |

where

(S20) |

is purely imaginary and anti-symmetric about . One immediately confirms that

(S21) |

Tracing over the environmental degrees of freedom thus yields

(S22) |

We now perform the continuum limit , leading to

(S23) |

In the final step, we have used that is a symmetric function about , while is anti-symmetric. Hence, we confirm

(S24) |

#### .2.3 Total trace distance

The difference between the two states and is given by

(S25) |

Diagonalization yields eigenvalues

(S26) |

and eigenstates

(S27) |

The trace norm of , which corresponds to the trace distance of and can thus be written as

(S28) |

Performing the continuum limit finally leads to

(S29) |

which for the parameter used in the experiment () yields the result given in Eq. (7) of the main manuscript.

### .3 Open-system evolution depending on initial correlations

For the dynamical local detection of the correlations in , we subject both and to the Michelson delay setup (see main manuscript), generating the dynamics

(S30) |

where correponds to the coordinate of the delayable mirror (M2) and the half-wave plate (HWP2) is rotated to the angle , generating the basis vectors

(S31) |

The original polarization vectors and can be expressed in terms of the new basis vectors and as

(S32) |

Hence, the evolution of the original states is given by

(S33) |

#### .3.1 Local trace distance

The operator evolves under this dynamics as

(S34) |

Since we only observe the dynamics of the polarization states, we take the partial trace, which yields

(S35) |

Represented as a matrix in the basis , this reads

(S36) |

where

(S37) |

Note that this result is independent of . The trace norm of the operator (S36) corresponds to the trace distance of the locally accessible polarization states and and can be determined as

(S38) | ||||

The continuum limit then yields

(S39) |

By comparing this result to Eq. (S29), one readily confirms that the general inequality

(S40) |

here is satisfied as a consequence of the triangular inequality.

In the present experiment, the frequency distribution of is characterized by a Lorentzian spectrum

(S41) |

We have , and ps. The function is normalized and has its support to a good approximation in the positive domain since ,

(S42) |

Its Fourier transform is given by

(S43) |

with

(S44) |

With this, the local trace distance, Eq. (S39), reduces to

(S45) |

which for yields Eq. (9) of the main manuscript.