Nonclassical Preparation of Quantum Remote States
Abstract
Remote state preparation (RSP) enables a sender to remotely prepare the quantum state of a receiver without sending the state itself. Recently, it has been recognized that quantum discord is a necessary resource for RSP. Here, we theoretically and experimentally investigate whether RSP can outperform dynamic classical remote state preparation processes. We show that such classical processes can describe certain RSPs powered by quantum discord. Rather, we argue that a new kind of Einstein-Podolsky-Rosen steering for dynamical processes, called quantum process steering, is the resource required for performing nonclassical RSP. We show how to measure quantum process steering by experimentally realizing nonclassical RSP of photonic quantum systems. Moreover, we demonstrate the transition from classical to quantum RSP. Our results also have applications in realizing genuine quantum RSP for quantum-enabled engineering.
In the RSP protocol Pati00 ; Bennett01 , Alice and Bob initially share an Einstein-Podolsky-Rosen (EPR) pair . Because of rotational symmetry the state can be rewritten as Cabello02 , where and are arbitrary pure states that constitute an orthonormal basis, and is an arbitrary single qubit unitary operator. If the quantum state, , on the equatorial plane of the Bloch sphere, is the target state envisaged by Alice, she can measure her qubit in the basis to prepare the state for Bob. Depending on the outcome of Alice’s measurement, one classical bit (cbit) is sent from Alice to Bob in order to tell him whether or not he needs to apply a rotation about the axis for state preparation to be achieved. For example, the state can be remotely prepared by using the RSP protocol with and a rotation operator : . See Fig. 1.
The RSP confers the obvious advantage of remotely manipulating Bob’s qubit. In particular, if Bob is unable to locally implement the target operator but only the simple rotation, Alice can help. Moreover, compared to teleportation Bennett93 , which requires Bell-state measurements Pirandola15 for transmitting an unknown qubit, the RSP protocol only needs local measurements on Alice’s particle. These two features makes RSP appealing and well adapted for a range of applications in quantum information and quantum engineering — from deterministically creating single-photon states Jeffrey04 to preparing single-photon hybrid entanglement Barreiro10 , initializing atomic quantum memory Rosenfeld07 , and assisting teleportation between atomic-ensemble quantum nodes Bao12 .

Identifying truly quantum RSP, irreproducible by any classical means, is not only significant in its own right, but can be useful for realizing the power of quantum-enabled technologies in which RSP plays a role Jeffrey04 ; Barreiro10 ; Rosenfeld07 ; Bao12 . Validation of RSP naturally leads to the following fundamental question: what are the physical resources actually required to realize nonclassical RSP?
More recently, it was discovered that quantum discord Ollivier01 ; Zurek00 ; Henderson01 , which measures nonclassical correlations between two systems in terms of the difference between two prescribed expressions for mutual information, serves as a resource for RSP Dakic12 . While maximally-entangled EPR pairs are assumed to be the default states used in the ideal RSP protocol, rather remarkably, states with no entanglement can outperform entangled states in RSP once they possess non-zero geometric quantum discord Dakic10 ; Luo10 . Importantly, this kind of verification relies on assumptions deduced from quantum theory. It does not generally hold for practical implementations where experimental imperfections can be described by classical physics. The same assumption is made in the fidelity benchmarks for RSP Killoran10 , which offer a seminal paradigm for ruling out certain cheating strategies. This issue turns out to hinge on whether it is possible to simulate RSP using classical methods.
Driven by the desire to understand the resources necessary for the nonclassical preparation of quantum remote states, we investigate to what extent RSP can be performed with generic classical dynamical processes. We introduce a model to capture the classical dynamic change to a state which can occur as the result of a remote-state preparation process. The scenario considered in this model generally goes beyond the assumptions made in the existing verification methods Dakic12 ; Killoran10 . We show that certain RSPs based on quantum discord can be described by such classical processes. A new kind of EPR steering Wiseman07 ; Reid09 is found and proven to be necessary for surpassing such classical preparation of quantum remote states. This enables nonclassicality of experimental implementations to be measured in terms of quantum process steering, which we will introduce below. Experimental demonstration using polarization-correlated photon pairs support our theoretical predictions.
Theory.—The whole RSP process can be considered as a quantum operation, , implementing the following state transformation:
(1) |
where is the initial state before the RSP process, and is the output state after the process occurs. It is worth emphasizing that RSP is nonclassical since performs a unitary quantum-mechanical transformation . Moreover, comprises all the operations and resources required to implement the RSP protocol. Therefore, any experimental imperfections in these necessary elements can affect the realization of . The worst case scenario is when a real-world RSP process may be interpreted as a classical process.
To demonstrate our classical RSP model, we first note that the input state can be represented by , where is the identity matrix, is a real three-dimensional vector (called the Bloch vector Nielsen00 , see Fig. 1) with , and consists of , , and Pauli matrices. After a RSP process , the remote state of Bob’s qubit is prepared in the state . Explicitly, the observables transform as and , where and are the eigenvalues and eigenvectors of , respectively. As shown above, for a given initial state , the remote state is determined by the RSP process applied to the eigenstates, .
In our model of classical RSP, we assume that the whole process is performed in the absence of EPR pairs. Furthermore, Alice is assumed to be unable to implement any unitary transformations. Instead, the classical RSP consists of the following three steps. First, there exists a recipe which specifies how the eigenstates are described by classical pre-existing states. This step represents Alice identifying which pre-existing state she has. The second step describes how the classical states change during the RSP process. The transition of classical states mimics the unitary transformation in the RSP task. Finally, the final classical states determine the qubit states that Bob has at the end of the classical RSP process.
A classical pre-existing state is a outcome set composed of the pre-existing outcomes corresponding to the observables under the assumption of classical realism Brunner14 , and is described by , where for is a unique index determining the set. In a classical RSP, denoted as , an initial state, say , is interpreted as a pre-existing (realistic) outcome of certain pre-existing states . For example, the classical states (for ) possess the same state property of . Moreover, causes the to classically transit to another classical pre-existing state. After the state transition has occurred, the final classical state, say , indicates a possible remote state for Bob’s “qubit”.
The RSP in the above scenario can be described by the map
(2) |
where is the probability of finding conditioned on the pre-existing outcome with , denotes Bob’s possible remote states, and is the probability of transition from to . The process then classically prepares an arbitrary input state as , where
(3) |
with a Bloch vector . This process of preparing remote states is classical in the sense that the initial states and their subsequent evolution are both classical, i.e., the essential elements in a dynamical map all follow a local and realistic classical theory.
The resulting state [Eq. (12)] of the classical map [Eq. (2)] is a generalization of the local hidden state (LHS) model Wiseman07 ; Reid09 , including the temporal analogue of the LHS model Chen14 ; Chen16 . The LHS model demonstrates how Alice uses a classical strategy to affect the states of Bob’s qubit without EPR pairs; the invalidation of this model shows EPR steering of the shared pair. In contrast, Eq. (12) describes the relationship between all possible inputs from Alice and outputs for Bob of a classical dynamical process designed to simulate a RSP task. In addition to Alice’s measurement, which is only the action taken in a steering scenario, Eq. (12) emulates all the necessaties required to realize the RSP protocol, including Alice’s manipulation related to , classical communication, and Bob’s correction operations.
Compared to the conventional model describing Wiseman07 ; Reid09 and quantifying Skrzypczyk2014 ; piani2015 EPR steering, the generalized LHS model is stricter because it is a larger class consisting of LHS models and the criterion describing dynamical processes. These constraints make the generalized LHS model effective in examining all possible output states for all in a RSP task. Any with outputs that cannot be represented in the form [Eq. (12)] for any pure inputs implements a nonclassical map, which is called quantum process steering. See Supplemental Material SuppMaterial for detailed comparisons between this generalized LHS model [Eq. (12)] and the standard LHS model.


To evaluate how well a given experimental process surpasses the generalized LHS model, we compare its output, , with all possible outputs from classical processes, , for all pure input states . We illustrate four different methods to measure and (or) identify the capability of to exhibit quantum process steering, called process steerability.
(i) Quantum composition :
(4) |
where describes the minimum proportion of as being the that cannot be represented by .
(ii) Quantum robustness :
(5) |
where is the minimum amount of the noise added such that the becomes classical.
(iii) Average-state-fidelity criterion for quantum RSP:
(6) |
where , and the integral is over the uniform measure on state space and .

(iv) Complementary-state-fidelity criterion for quantum RSP:
(7) |
where for denote the state fidelities over arbitrary two orthonormal and complementary state bases, and , with for .
In measuring or , or determining the threshold , the optimization problem with respect to the dynamical maps can be solved by semi-definite programming SuppMaterial . The criterion [Eq. (19)] is based on the fact that, since Hofmann05 , beating the classical limit implies a quantum RSP. While the methods defined in (i)-(iii) can characterize the whole RSP process for all pure input states , practical experimental realizations require complete process tomography Nielsen00 . In contrast, implementing the criterion (iv), Eq. (19), is efficient and requires only two complementary sets of input states. For example, to identify experimental RSPs close to a target of , one can follow the RSP protocol and examine the experimental results with the criterion [Eq. (19)], using , for , as input states.
We now consider the process steerability of the ideal RSP process for arbitrary . In this ideal case, has a quantum composition of , a quantum robustness of , and for the fidelity criteria [Eqs. (18) and (19)], gives . These values will decrease when contaminated with noise. For example, when the Werner state , is used for RSP to represent noise in the state creation process, the process steerability decreases with increasing noise intensity and disappears at , while the underlying Werner state still possesses a non-zero quantum discord.
Process steerability serves as the resource required to yield a higher payoff for nonclassical RSP. The fidelity thresholds in the criteria [Eqs. (18) and (19)] are stricter than the fidelity threshold, , for the RSPs enabled by quantum discord Dakic12 . Furthermore, compared with existing fidelity benchmarks Killoran10 , our formalism allows, in the classical regime, for imperfections in Alice’s results, and removes the assumption that can cause false-positive results in the benchmarks proposed in Ref. Killoran10 , i.e., when Bob knows the target states in the cheating strategies. Therefore, our fidelity criteria eases the necessary requirements of high state fidelity in the benchmarks.
Experiment and results.—In the experiment, we generate photon pairs correlated in the polarization degree of freedom. A schematic diagram of our experimental set-up is shown in Fig. 2. We first implement four different kinds of RSP, denoted as , while we denote the ideal RSP process as , which are identified by the choice of single-qubit unitaries, . The experimental results for the criteria outlined in (i)-(iv), and the concrete comparisons between , and are shown in Fig. 3. The experimental RSPs exhibit quantum process steering and are seen to be close to the target RSP according to the fidelity criteria, which confirms the realizations of nonclassical RSPs.
In Fig. 4(a), we demonstrate how the photon walk-off effect affects the process steerability of the resulting RSPs. The walk-off effect describes that the ordinary and the extraordinary photons created via type-II SPDC possess different velocities inside the BBO crystal. Such effect makes generated photon pairs separable and the resulting experimental RSP classical. The walk-off effect can be compensated by a HWP and a CBBO in each mode and (see Fig. 2). By varying the compensation of the walk-off of the photons from SPDC, the generated photon pairs can contribute to different process steerabilities. Figure 4(b) shows that experimental process steerability varies with the mixedness of the Werner state . To experimentally set the noise weight of , the Werner states were generated by appropriately arranging the measurement durations Lavoie10 ; Amselem09 .
In Fig. 4(a), we show an example where a photon pair with non-zero quantum discord and EPR steerability results in an RSP that does not satisfy the fidelity criteria.The experimental results in Fig. 4(b) illustrate further that can even describe certain RSP processes which rely on quantum discord. They are consistent with the theoretical analysis showing that quantum process steering enables quantum RSP.
We have shown that quantum process steering is the necessary resource that enables RSP to surpass classical dynamical methods of preparing remote states. We experimentally demonstrated this new kind of EPR steering, via several criteria, using different types of polarization-correlated photon states. Our demonstration also shows that other quantum characteristics of polarization-correlated photon pairs, such as quantum discord and standard EPR steering, do not guarantee the existence of nonclassical RSP. Our results outline how to achieve genuine quantum RSP, and are of importance for applications where ruling out any classical mimicries of RSP is necessary. They may provide a useful insight into identifying nonclassical processes for other quantum-information tasks, such as universal one-way quantum computing raussendorf2001 ; briegel2009 ; you2007 ; tanamoto2009 ; wang2010 or quantum communication in quantum networks Chou18 ; Pirker18 .
Acknowledgement.—We acknowledge the fruitful discussions with He Lu. This work was partially supported by the Ministry of Science and Technology, Taiwan, under Grant No. MOST 107-2628-M-006-001-MY4. F.N. is supported in part by: NTT Research, Army Research Office (ARO) (Grant No. W911NF-18-1-0358), Japan Science and Technology Agency (JST) (via the CREST Grant No. JPMJCR1676), Japan Society for the Promotion of Science (JSPS) (via the KAKENHI Grant No. JP20H00134, and the grant JSPS-RFBR Grant No. JPJSBP120194828), and the Grant No. FQXi-IAF19-06 from the Foundational Questions Institute Fund (FQXi), a donor advised fund of the Silicon Valley Community Foundation. N.L. is partially supported by JST PRESTO through Grant No. JPMJPR18GC.
References
- (1) A.K. Pati, Minimum classical bit for remote preparation and measurement of a qubit, Phys. Rev. A 63, 014302 (2000).
- (2) C. H. Bennett, D. P. DiVincenzo, P. W. Shor, J. A. Smolin, B. M. Terhal, and W. K. Wootters, Remote state preparation, Phys. Rev. Lett. 87, 077902 (2001).
- (3) A. Cabello, N-particle N-level singlet states: some properties and applications, Phys. Rev. Lett. 89, 100402 (2002).
- (4) C. H. Bennett, G. Brassard, C. Crépeau, R. Jozsa, A. Peres, and W. K. Wootters, Teleporting an unknown quantum state via dual classical and Einstein-Podolsky-Rosen channels, Phys. Rev. Lett. 70, 1895 (1993).
- (5) S. Pirandola, J. Eisert, C. Weedbrook, A. Furusawa, and S. L. Braunstein, Advances in quantum teleportation, Nat. Photonics 9, 641 (2015).
- (6) E. Jeffrey, N. A. Peters, and P. G. Kwiat, Towards a periodic deterministic source of arbitrary single-photon states, New J. Phys. 6, 100 (2004).
- (7) J. T. Barreiro, T.-C. Wei, and P. G. Kwiat, Remote preparation of single-photon“hybrid” entangled and vector-polarization states, Phys. Rev. Lett. 105, 030407 (2010).
- (8) W. Rosenfeld, S. Berner, J. Volz, M. Weber, and H. Weinfurter, Remote preparation of an atomic quantum memory, Phys. Rev. Lett. 98, 050504 (2007).
- (9) X.-H. Bao, X.-F. Xu, C.-M. Li, Z.-S. Yuan, C.-Y. Lu, and J.-W. Pan, Quantum teleportation between remote atomic ensemble quantum memories, Proc. Natl. Acad. Sci. U.S.A. 109, 20347 (2012).
- (10) H. Ollivier and W. H. Zurek, Quantum discord: A measure of the quantumness of correlations, Phys. Rev. Lett. 88, 017901 (2001).
- (11) W. H. Zurek, Einselection and decoherence from an information theory perspective, Ann. Phys. (Leipzig) 9, 855 (2000).
- (12) L. Henderson and V. Vedral, Classical, quantum and total correlations, J. Phys. A 34, 6899 (2001).
- (13) B. Dakić et al., Quantum discord as resource for remote state preparation, Nat. Phys. 8, 666 (2012).
- (14) B. Dakić, V. Vedral, and Č. Brukner, Necessary and sufficient condition for nonzero quantum discord, Phys. Rev. Lett. 105, 190502 (2010).
- (15) S. Luo and S. Fu, Geometric measure of quantum discord, Phys. Rev. A 82, 034302 (2010).
- (16) N. Killoran, D. N. Biggerstaff, R. Kaltenbaek, K. J. Resch, and N. Lütkenhaus, Derivation and experimental test of fidelity benchmarks for remote preparation of arbitrary qubit states, Phys. Rev. A 81, 012334 (2010).
- (17) H. M. Wiseman, S. J. Jones, and A. C. Doherty, Steering, entanglement, nonlocality, and the Einstein-Podolsky-Rosen paradox, Phys. Rev. Lett. 98, 140402 (2007).
- (18) M. D. Reid, et al., Colloquium: the Einstein-Podolsky-Rosen paradox: from concepts to applications, Rev. Mod. Phys. 81, 1727 (2009).
- (19) M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge Univ. Press, UK, 2000).
- (20) N. Brunner, D. Cavalcanti, S. Pironio, V. Scarani, and S. Wehner, Bell nonlocality, Rev. Mod. Phys. 86, 419 (2014).
- (21) Y.-N. Chen, C.-M. Li, N. Lambert, S.-L. Chen, Y. Ota, G.-Y. Chen, and F. Nori, Temporal steering inequality, Phys. Rev. A 89, 032112 (2014).
- (22) S.-L. Chen, N. Lambert, C.-M. Li, A. Miranowicz, Y.-N. Chen, and F. Nori, Quantifying non-Markovianity with temporal steering, Phys. Rev. Lett. 116, 020503 (2016).
- (23) P. Skrzypczyk, M. Navascués, and D. Cavalcanti, Quantifying Einstein-Podolsky-Rosen steering, Phys. Rev. Lett. 112, 180404 (2014).
- (24) M. Piani and J.Watrous, Necessary and sufficient quantum information characterization of Einstein-Podolsky-Rosen steering, Phys. Rev. Lett. 114, 060404 (2015).
- (25) See Supplemental Material for the comparison between the classical RSP model and the LHS model, and the details of solving optimization problems of the , , and . Moreover, the theoretical model for the compensation of the photon walk-off effect acting on noisy created photon pairs is detailed. Finally, we describe how to measure the discord and EPR steerability of polarization-correlated photon pairs created in our experiments.
- (26) H. F. Hofmann, Complementary classical fidelities as an efficient criterion for the evaluation of experimentally realized quantum operations, Phys. Rev. Lett. 94, 160504 (2005).
- (27) P. G. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, A. V. Sergienko, and Y. Shih, New high-intensity source of polarization-entangled photon pairs, Phys. Rev. Lett. 75, 4337 (1995).
- (28) J. Lavoie, R. Kaltenbaek, M. Piani, and K. J. Resch, Experimental bound entanglement in a four-photon state, Phys. Rev. Lett. 105, 130501 (2010).
- (29) E. Amselem and M. Bourennane, Experimental four-qubit bound entanglement, Nat. Phys. 5, 748 (2009).
- (30) R. Raussendorf and H.J. Briegel, A one-way quantum computer, Phys. Rev. Lett. 86, 5188 (2001).
- (31) J. Q. You, X.-B. Wang, T. Tanamoto, and F. Nori, Efficient one-step generation of large cluster states with solid-state circuits, Phys. Rev. A 75, 052319 (2007).
- (32) H. J. Briegel, D. E. Browne, W. Dr, R. Raussendorf, and M. Van den Nest, Measurement-based quantum computation, Nat. Phys. 5, 19 (2009).
- (33) T. Tanamoto, Y. X. Liu, X. Hu, and F. Nori, Efficient quantum circuits for one-way quantum computing, Phys. Rev. Lett. 102, 100501 (2009).
- (34) H. Wang, C.-P. Yang, and F. Nori, Robust and scalable optical one-way quantum computation, Phys. Rev. A 81, 052332 (2010).
- (35) K. S. Chou et al., Deterministic teleportation of a quantum gate between two logical qubits, Nature (London) 561, 368 (2018).
- (36) A. Pirker, J. Wallnöfer, and W. Dür, Modular architectures for quantum networks, New J. Phys. 20, 053054 (2018).
Supplementary Information
I Comparisons between the generalized LHS model and the standard LHS model
The local hidden state (LHS) model Wiseman07 ; Reid09 describes the state of Bob’s qubit conditioned on Alice’s measurement:
(8) |
where denotes the th measurement outcome of Alice’s th measurement. Alternatively, this can be represented by unnormalized states and rephrased as
(9) |
where and , and is the probability of finding the outcome conditioned on the classical state . The relation of has be used to derive Eq. (9) from Eq. (8). Given Alice’s choice of measurement and outcome , any assemblage of unnormalized states which Alice’s measurement affects Bob’s particle into which cannot be represented in the form (9) is called steerable Wiseman07 ; Reid09 .
To discuss the difference between our classical model and the LHS model, we first note that, given the input state , a classical RSP, , generates an output described by Eq. (2) in the main text:
(10) |
which can be written as
(11) |
under the condition of . When comparing Eq. (8) and Eq. (11) we summarize the following two points:
(i) While they are of the same form and Eq. (11) can be considered as a LHS model, cannot be fully defined by considering only some of the output states under the inputs in the LHS model via Eq. (11). Rather, since the is a dynamical process, it can be completely defined only when output states for all the pure input states are taken into account. Note that, it is unnecessary to include mixed input states because they can be generated by incoherently mixing pure states.
(ii) From the task-oriented viewpoint, any assemblage of output states derived from an experimental RSP process which cannot described by
(12) |
as shown in Eq. (3) in the main text, is called process steerable, which corresponds to the term introduced in the main text: quantum process steering.
For the above reasons, the resulting state (12) of the classical map therefore is a generalized LHS model. Such a generalization makes identifying and quantifying quantum process steering rather different from the methods for quantitatively characterizing EPR steering. See the next section for further detailed discussions.
The above analysis, comparisons and conclusion are applicable to the temporal LHS model for temporal steering Chen14 ; Chen16 . In the temporal scenario, Alice performs a measurement on a single system in a certain initial state at time . After the measurement, the initial state becomes, say . Then this state is sent into a quantum channel for a time . At time , Bob receives the system of the state . Following the un-normalized assemblage Wiseman07 ; Reid09 in the standard LHS model, the un-normalized states in temporal scenario is defined by , where the superscript reminds one that the assemblage is for temporal steering. In the temporal LHS model, an unsteerable assemblage is defined as
(13) |
where and is created from a source which determines the possible correlations between Alice’s (at time ) and Bob’s measurement results (at the time ) under classical realism. By using the same method as shown above, Eq. (13) can be rephrased as
(14) |
Still, when comparing Eq. (11) to Eq. (14), cannot be fully defined by considering the finite output states under the inputs in the temporal LHS model via Eq. (11). The can be defined only when the output states, , for all the pure input states are considered, where
(15) |
II Quantum composition, quantum robustness, and the fidelity criteria
In this section, we will introduce how to calculate quantum composition , quantum robustness , and the fidelity bound in the fidelity criteria. The methods introduced in the main text to measure and (or) identify process steerability of an experimental process are summarized as follows:
(i) Quantum composition . For all possible pure states (i.e., ), the corresponding prepared remote states, , can be decomposed as
(16) |
where the quantum composition describes the minimum proportion of that cannot be represented by the classically prepared states . An ideal RSP shows that is fully composed of the that cannot be represented by , i.e., .
(ii) Quantum robustness . In this case the prepared states can be written as
(17) |
where the quantum robustness is the minimum amount of the noise added such that the becomes a classically prepared remote state. When , the experimental RSP is already classical. Increasing increases the difference between and .
(iii) Average state fidelity. The is identified as being faithful to show quantum process steering if the state fidelity of and : , averaged over all possible input pure states is larger than the maximum averaged state fidelity derived from classical RSP, , i.e.,
(18) |
where the integral is over the uniform measure on state space and .
(iv) Complementary state fidelity. Suppose that and are two orthonormal and complementary state bases, where for . Let us consider the state fidelities for . The is faithful and possess quantum process steerability if the complementary state fidelity Hofmann05 beating the classical limit implies quantum RSP:
(19) |
In measuring the , , and determining the threshold and , it is necessary to consider the contributions from all pure states in the respective computational tasks. To solve such questions, we first note that Eq. (16) and Eq. (17) can be rephrased as
(20) |
and
(21) |
respectively, where and . Then the optimization problems with respect to states become to the optimization of and with respect to dynamical maps Hsieh17 . Moreover, it has been shown that Gilchrist05 ; Hofmann05 , implying that the can be derived from , i.e., the process fidelity of and . Equivalently, the criterion (18) becomes , where the is obtained by considering the best similarity between and . Finally, as is experimentally measured, , and , can therefore be determined by semi-definite programming (SDP) Lofberg ; sdpsolver and the can be calculated through .
In practical experiments, the experimental process can be fully described by using the process tomography algorithm Nielsen00 ; Chuang97 and then be represented by a positive Hermitian matrix, called process matrix, . The process matrix is constructed by the specific output states which are the output states of the eigenstates of three Pauli matrices. Considering the output states from a classical RSP , we can use process tomography algorithm to get a process matrix for the classical RSP . Thus, Eq. (20) and Eq. (21) can be rephrased as
(22) |
and
(23) |
respectively, and the can be represented as , where is the process matrix of and and are the process matrices of and respectively. The process matrices and can be used to solve the optimization problems of the optimization of , , and . To calculate , , and via SDP, we let the processes in SDP satisfy the definition of process matrices or density matrices by the following constraint. A process matrix and a density matrix must be positive semi-definite, i.e., and . Since a process matrix is constructed of the output states and the output states of the for classical RSP model is determined by Bob’s possible remote states in Eq. (2) in the main text, the process matrix must satisfy .
According to Eq. (22), by representing as a unnormalized process matrix , quantum composition can be obtained by minimizing the following quantity via SDP with MATLAB Lofberg ; sdpsolver : , under the following constraints:
The first constraint ensures that the must be positive semi-definite since it is a process matrix. The second constraint ensures that all possible remote states in classical RSP model are positive semi-definite. The third constraint makes the process matrix of satisfy the definition of process matrices which need to be positive semi-definite.
As the method used to solve , to calculate the quantum robustness , we define . Thus, according to Eq. (23), the quantum robustness of can be obtained by using SDP to solve . The constraints for solving via SDP are
The first and second constraints are the same as the constraints used in calculating . The third constraint ensures that . The fourth constraint is the condition of positive semidefiniteness for in Eq. (23).
With the process matrix of , the fidelity bounds can be solved under the constraints
which make the be a normalized matrix to calculate fidelity. While we get the fidelity bound , the average-state-fidelity bound and the complementary-state-fidelity bound can be calculated through and .
III Experimental photon pairs
The experimental set-up for generating polarization-correlated photon pairs to achieve certain RSP is shown in Fig. 2 in the main text. A type-II -barium borate (BBO) crystal with 2 mm thickness is pumped by ultraviolet laser pulses with a central wavelength of 390 nm, pulse duration of 135 fs, repetition rate of 76 MHz, and power of 250 mW to generate polarization-correlated photon pairs at a central wavelength of 780 nm via spontaneous parametric down-conversion (SPDC) process Kwiat95 . The photon pairs then pass through birefringent compensators, i.e., two half-wave plates (HWP) and two correction BBO (CBBO) crystals with 1 mm thickness, to eliminate the longitudinal and spatial walk-off effects. Note that the setting angle of the HWPs in both mode and mode are under ideal compensation circumstance, i.e., . Subsequently, we utilize a phase controller consisting of two quarter-wave plates (QWP) and a HWP to obtain the desired state, . Finally, the photon pairs are filtered by narrow-band filters with 3 nm bandwidth and are detected by silicon avalanche single-photon detectors.
The resulting state of the created photon pairs, denoted as , can be experimentally characterized via state tomography using a pair of polarization analyzers; each consisting of a QWP, a HWP, and a polarizer, in mode and mode , respectively. The measured density matrix of the created photon pairs after the designed compensation () is of the following form:
(24) |
where the fidelity of and the target state, , is in our experiment. The maximum-likelihood technique James01 has been used to make the desnsity matrix reasonable.
The created entangled photon pairs of the state is used to perform nonclassical RSPs, , for (a) , (b) , (c) , and (d) shown in Fig. 3 in the main text. Moreover, is also used to demonstrate in the case of the desired compensation () in Fig. 4(a) and in the case of Werner state at in Fig. 4(b) in the main text.
IV Noise model for experimental photon pairs
As shown in Eq. (24), there exist undesired components in the created photon pairs, mainly due to the imperfection of BBO crystals, the higher-order terms of down-converted photons, the imperfection of optical components, the accidental coincidences, and the imperfect detection efficiency, which cannot be corrected by compensating the photon walk-off effect. In other words, despite the use of the HWPs and the CBBO crystals for compensating the photon walk-off effect, the polarization state of the created photon pairs is not the ideal state in our experiment. To quantitatively describe the resulting noise and imperfection that are present in our experiment, we effectively represent the state in terms of an ideal state mixed with white noise. Therefore, after the compensation of the photon walk-off effect, the created photon pairs in our noise model is described by
(25) |
where denotes the intensity of white noise in created photon pairs. The in our experiment can be estmated by maximizing the fidelity between and the experimentally created photon pairs, (24). The fidelity between them is defined as
After performing the maximization task, the fidelity between and is , and the estimated noise intensity is . Therefore, the experimental photon pairs in our noise model is
(26) |
We utilize the effective noise state to calculate the theoretical values of , , and at in Fig. 4(a) in the main text. Furthermore, as shown in Fig. 4(b) in the main text, the Werner state in our noise model is described as
(27) |
where is the noise weight in Werner state. We utilize to estimate the theoretical values of , , , and quantum discord in Fig. 4(b) in the main text. It is worth noting that the experimental results are highly consistent with the theoretical values. It can be inferred that is a good estimation of the created photon pairs in Werner state varying with the noise weight.
V Compensation of the photon walk-off effect
While a pumped laser beam enters the BBO crystal, it is partially converted into a pair of mutually orthogonally polarized photons, the ordinary and extraordinary beams which are respectively horizontally and vertically polarized. The two-photon state at the output surface of the BBO crystal can be described as Rubin94
(28) |
where the direction is assumed to be parallel to the pump beam, is the constant comparing with SPDC efficient; , and are the frequencies of the pump, ordinary, and extraordinary beams with the wave vectors , , and , respectively; is the spectral profile of the pump pulse; -function shows the energy conservation in SPDC; mm is the length of the BBO crystal; and are the creation operators for the ordinary and extraordinary beams respectively; represents the vacuum state. Since the detectors are in the modes and and narrow bandwidth filters are added in front of the detectors to select frequencies , Eq. (28) can be simplified as
where is the normalization constant, and , are the creation operators of photon with polarized state for spatial mode , respectively. After the normalization, the polarized state of the generated photon pairs is , where 0 and 1 indicate horizontal and vertical polarizations. Since the BBO crystal is a birefringent crystal and the ordinary and extraordinary beams have orthogonal polarization, a longitudinal and spatial walk-off between the ordinary and extraordinary beams is induced by the BBO crystal Poh_thesis ; Lee16 . We assume that SPDC occurs in the center of the BBO crystal, , the walk-off effect makes horizontally polarized photons be delayed by a time with respect vertically polarized photons, where and denote the velocities of the ordinary wave and extraordinary wave in the BBO crystal, respectively. When the photon pairs exit the BBO crystal, the walk-off effect makes the state becomes
(29) |
where denotes the photon with polarized state is delayed time . The delay time makes and distinguishable. Since we use entangled polarized state in our RSP experiment, we perform a partial trace over the time system, the polarized state of photon pairs is a completely mixed state
(30) |
To compensate the photon walk-off effect, a HWP and a CBBO crystal are used in each mode and . See Fig. 2 in the main text. The unitary transform of HWP for different setting angle is
By setting the angle of the HWP in each mode and at the same angle , i.e., , the output state from the BBO crystal becomes
(31) |
While the photon pairs go through the CBBO crystals which make horizontally polarized photons be delayed by a time respect vertically polarized photons, the final state after the CBBO crystals is
(32) |
Then, we use two QWPs and one HWP in mode to control the phase to become . Thus, after the partial trace over the time system, the generated polarized state of the photon pairs is
(33) |
where . We utilize to calculate , , and in Fig. 4(a) in the main text.
If the setting angle is , the unitary
makes the horizontal and vertical polarization be exchanged. Thus, Eqs. (31) and (32) for are
and
respectively. Since all photons of a pair have the same delay which means the photon walk-off effect is compensated after the CBBO crystals, after the partial trace over the time system, for , the generated state in Eq. (33) is the ideal creation of the state . By contrast, if the setting angle , the HWP do not exchanged the horizontal and vertical polarization and the photon walk-off remains after photon pairs pass through the CBBO crystals. Since the walk-off does not be compensated, the states and remain distinguishable and the generated state is a completely mixed state.
Combining the noise model given in Eq. (26) and the polarization state of the created photon pairs varying with the photon walk-off effect as shown in Eq. (33), the theoretical estimation of the created photon pairs is described by
(34) |
As shown in Fig. 4(a) in the main text and will be discussed in the following section, we utilize this theoretical model to calculate the theoretical values of , , , quantum discord and EPR steerability.
VI Implementation of RSP protocol and resulting theoretical prediction of the process steerability
As shown in Fig. 4(a) in the main text, we realize the RSP protocol with respect to a target RSP, , of . To obtain the theoretical prediction of the process steerability of the RSPs varying with the walk-off compensation (), we utilize to achieve the RSPs and obtain the resulting process matrices. In experimental RSPs, the process matrices were measured by using the process tomography technique Nielsen00 and the optical setup shown in Fig. 2 in the main text. The main steps for implementing the RSP protocol are as follows.
First, as illustrated in Fig. 2 in the main text, we set the HWP which is after the phase controller and before the polarization analyzer in mode at to realize the operation required on Alice’s side.
Second, we prepare the states as input states by the polarization analyzer consisting of a QWP, a HWP, and a polarizer in mode . Note that are eigenvectors of for and shown in the main text. The preparation basis is determined by adjusting the HWP and the QWP at different angles. For example, by adjusting the HWP and the QWP at angles of with respect to the vertical axis, the preparation basis for the polarization states is set for . Similarly, the preparation bases are set for and by adjusting the HWP and the QWP at angles of and , respectively. Then, we set the polarizer at and to prepare and , respectively.
Third, we perform state tomography to obtain the density matrix of output state by the polarization analyzer in mode . As mentioned in the previous step, the measurement basis is also determined by the components of the polarization analyzer in mode . Taking , for example, our goal is to obtain the density matrix of . We prepare the input state and by adjusting the HWP, the QWP, and the polarizer at angles of and in mode , respectively. Then, we perform the measurement of , , and on in mode and record the coincident counts of them conditioned on and prepared in mode . The density matrix of and can then be constructed through coincidence measurement results via state tomography. According to the RSP protocol, the density matrix of can be obtained by the mean of and after correction, i.e.,
Fourth, considering the theoretical measurement results of , we can get the theoretical output states for and . Therefore, the theoretical process matrix of the RSP can be obtained through process tomography Nielsen00 . According to , the process fidelity and the process steerability , can be calculated via SDP Lofberg ; sdpsolver . In addition, the theoretical average state fidelity can be obtained by .
As shown in Fig. 4(a) in the main text, the experimental results are consistent with the theoretical predictions of , , and . Therefore, we can conclude that is a good estimation of the state of the experimental photon pairs varying with the photon walk-off effect.
VII quantum discord and EPR steerability of the created photon pairs
In this section, we utilize the geometric discord Dakic10 and steerable weight skrzypczyk2014 to measure quantum discord and EPR steerability of the created photon pairs, respectively. We will show how to calculate the geometric discord and steerable weight in the following subsections.
VII.1 Geometric discord
For a given , the geometric discord is calculated by the Bloch vector with and the correlation tensor with Dakic10 . Here, is composed of , , and Pauli matrices. The geometric discord of takes the following form
(35) |
where the superscript indicates the transpose operation and is the largest eigenvalue of matrix . Thus, we can obtain the experimental geometric discord, denoted as , by .
As the results shown in Fig. 4(a) in the main text, the density matrix of the created photon pairs at in our experiment is of the following form:
(36) |
where the fidelity of and the target state, , is in our experiment. The resulting geometric discord of is . Similarly, the theoretical geometric discord, , can be obtained by . With and , the state of the created photon pairs has quantum discord in terms of geometric discord, however the resulting RSP cannot satisfy the fidelity criteria (6) in the main text.
As shown in Fig. 4(b) in the main text, the density matrix of the created photon pairs at in our experiment can be represented as
(37) |
where the fidelity of and the target state, , is in our experiment. The resulting geometric discord of is and the theoretical geometric discord of is . It can be seen that the created photon pairs possess quantum discord in terms of geometric discord but the process steerability disappears. That is, it is possible that a RSP which is based on the resource of quantum discord can be described by classical RSP, .
VII.2 Steerable weight
Given an assemblage of unnormalized states which Alice’s measurement affects Bob’s particle into, steerable weight is defined as the minimum amount of steerable resource in that cannot be represented in the form Eq. (9). Any assemblage can be decomposed as
(38) |
where and () is the unsteerable (steerable) assemblage in which the elements can (cannot) be represented in the form Eq. (9). The steerable weight is then defined as , where denotes the maximum and can be obtained via SDP; see skrzypczyk2014 for more detail. Here, to compare with quantum process steering that is measured by process matrix of the RSP protocol constructed via process tomography algorithm, we choose the measurements of the three Pauli matrices, , , and for , respectively. Through the density matrices and , we can get the experimental and theoretical steerable weight, and , respectively. As illustrated in Fig. 4(a) in the main text, the steerable weight of shown in Eq. (37) is which implies that the state of the created photon pairs has EPR steerability in terms of steerable weight, however the resulting RSP cannot satisfy the fidelity criteria (6) in the main text.
References
- (1) H. M. Wiseman, S. J. Jones, and A. C. Doherty, Steering, entanglement, nonlocality, and the Einstein-Podolsky-Rosen paradox, Phys. Rev. Lett. 98, 140402 (2007).
- (2) M. D. Reid, et al., Colloquium: the Einstein-Podolsky-Rosen paradox: from concepts to applications, Rev. Mod. Phys. 81, 1727 (2009).
- (3) Y.-N. Chen, C.-M. Li, N. Lambert, S.-L. Chen, Y. Ota, G.-Y. Chen, and F. Nori, Temporal steering inequality, Phys. Rev. A 89, 032112 (2014).
- (4) S.-L. Chen, N. Lambert, C.-M. Li, A. Miranowicz, Y.-N. Chen, and F. Nori, Quantifying non-Markovianity with temporal steering, Phys. Rev. Lett. 116, 020503 (2016).
- (5) H. F. Hofmann, Complementary classical fidelities as an efficient criterion for the evaluation of experimentally realized quantum operations, Phys. Rev. Lett. 94, 160504 (2005).
- (6) J.-H. Hsieh, S.-H. Chen, and C.-M. Li, Quantifying quantum-mechanical processes, Sci. Rep. 7, 13588 (2017).
- (7) A. Gilchrist, N. K. Langford, and M. A. Nielsen, Distance measures to compare real and ideal quantum processes, Phys. Rev. A 71, 062310 (2005).
- (8) J. Löfberg, Yalmip: A toolbox for modeling and optimization in MATLAB, in proceedings of the CACSD conference, Taipei, Taiwan, 2004; MATLAB, http://users.isy.liu.se/johanl/yalmip/.
- (9) K. C. Toh, M. J. Todd, and R. H. Ttnc, SDPT3 – a MATLAB software package for semidefinite-quadratic-linear programming, version 4.0., Available: http://www.math.nus.edu.sg/mattohkc/sdpt3.html.
- (10) M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge Univ. Press, UK, 2000).
- (11) I. L. Chuang and M. A. Nielsen, Prescription for experimental determination of the dynamics of a quantum black box, J. Mod. Opt. 44, 2455 (1997).
- (12) P. G. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, A. V. Sergienko, and Y. Shih, New high-intensity source of polarization-entangled photon pairs, Phys. Rev. Lett. 75, 4337 (1995).
- (13) D. F. V. James, P. G. Kwiat, W. J. Munro, and A. G. White, Measurements of qubits, Phys. Rev. A 64, 052312 (2001).
- (14) M. H. Rubin, D. N. Klyshko, Y. H. Shih, and A. V. Sergienko, Theory of two-photon entanglement in type-II optical parametric down-conversion, Phys. Rev. A 50, 5122 (1994).
- (15) P. H. Shun, Towards a high quality polarization-entangled multi-photon source, A Thesis Submitted for the Degree of Master of Science, Department of Physics, National University of Singapore (2009).
- (16) S. M. Lee, H. Kim, M. Cha, and H. S. Moon, Polarization-entangled photon-pair source obtained via type-II non-collinear SPDC process with PPKTP crystal, Opt. Express 24, 2941 (2016).
- (17) B. Dakić, V. Vedral, and Č. Brukner, Necessary and sufficient condition for nonzero quantum discord, Phys. Rev. Lett. 105, 190502 (2010).
- (18) P. Skrzypczyk, M. Navascués, and D. Cavalcanti, Quantifying Einstein-Podolsky-Rosen steering, Phys. Rev. Lett. 112, 180404 (2014).