This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Nonclassical Preparation of Quantum Remote States

Shih-Hsuan Chen Department of Engineering Science, National Cheng Kung University, Tainan 70101, Taiwan Center for Quantum Frontiers of Research &\& Technology, National Cheng Kung University, Tainan 70101, Taiwan    Yu-Chien Kao Department of Engineering Science, National Cheng Kung University, Tainan 70101, Taiwan Center for Quantum Frontiers of Research &\& Technology, National Cheng Kung University, Tainan 70101, Taiwan    Neill Lambert Theoretical Quantum Physics Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Saitama 351-0198, Japan    Franco Nori Theoretical Quantum Physics Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Saitama 351-0198, Japan Department of Physics, University of Michigan, Ann Arbor, Michigan 48109-1040, USA    Che-Ming Li [email protected] Department of Engineering Science, National Cheng Kung University, Tainan 70101, Taiwan Center for Quantum Frontiers of Research &\& Technology, National Cheng Kung University, Tainan 70101, Taiwan Center for Quantum Technology, Hsinchu 30013, Taiwan
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 |Ψ=(|01|10)/2\ket{\Psi^{-}}=(\ket{01}-\ket{10})/\sqrt{2}. Because of rotational symmetry the state |Ψ\ket{\Psi^{-}} can be rewritten as |Ψ=(U|𝕤0U|𝕤0U|𝕤0U|𝕤0)/2\ket{\Psi^{-}}=(U\ket{\mathbb{s}_{0}}\otimes U\ket{\mathbb{s}_{0}^{\bot}}-U\ket{\mathbb{s}_{0}^{\bot}}\otimes U\ket{\mathbb{s}_{0}})/\sqrt{2} Cabello02 , where |𝕤0\ket{\mathbb{s}_{0}} and |𝕤0\ket{\mathbb{s}_{0}^{\bot}} are arbitrary pure states that constitute an orthonormal basis, and UU is an arbitrary single qubit unitary operator. If the quantum state, |𝕤=(|0+eiϕ|1)/2\ket{\mathbb{s}}=(\ket{0}+e^{i\phi}\ket{1})/\sqrt{2}, on the equatorial plane of the Bloch sphere, is the target state envisaged by Alice, she can measure her qubit in the basis {U|𝕤0,U|𝕤0}\{U\ket{\mathbb{s}_{0}},U\ket{\mathbb{s}_{0}^{\bot}}\} to prepare the state |𝕤=U|𝕤0\ket{\mathbb{s}}=U\ket{\mathbb{s}_{0}} 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 π\pi rotation about the zz axis for state preparation to be achieved. For example, the state |𝕤\ket{\mathbb{s}} can be remotely prepared by using the RSP protocol with |𝕤0=|0\ket{\mathbb{s}_{0}}=\ket{0} and a rotation operator UU: R(ϕ)=(|00|+eiϕ|10|+|01|eiϕ|11|)/2R(\phi)\!=\!(\ket{0}\!\!\bra{0}\!+\!e^{i\phi}\ket{1}\!\!\bra{0}\!+\!\ket{0}\!\!\bra{1}\!-\!e^{i\phi}\ket{1}\!\!\bra{1})\!/\!\sqrt{2}. 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 UU but only the simple π\pi 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 .

Refer to caption
Figure 1: Remote state preparation (RSP). To prepare Bob’s remote state |𝕤\ket{\mathbb{s}} with a desired operation UU, Alice first performs the operation UU^{{\dagger}} and then measures her qubit (part of the state |Ψ\ket{\Psi^{-}}) in the basis {|𝕤0,|𝕤0}\{\ket{\mathbb{s}_{0}},\ket{\mathbb{s}_{0}^{\bot}}\}. Suppose |𝕤0=|0\ket{\mathbb{s}_{0}}=\ket{0}, |𝕤0=|1\ket{\mathbb{s}_{0}^{\bot}}=\ket{1}, as shown on the Bloch sphere by the Bloch vectors, and U=R(ϕ)U=R(\phi) is a rotation operator. If Alice obtains the result |0\ket{0} (|1\ket{1}), she sends Bob the message +1+1 (1-1) via a classical channel to show that his qubit is in the prepared state: R(ϕ)|0=(|0+eiϕ|1)/2R(\phi)\ket{0}=(\ket{0}+e^{i\phi}\ket{1})/\sqrt{2}, without correction, denoted as the identity II in the figure (with a correction of the Pauli-ZZ operation). It is worth noting that, while Bob may know the UU, the prepared states |𝕤\ket{\mathbb{s}} are unknown to him as the |𝕤0\ket{\mathbb{s}_{0}} is unknown as well. The above procedure equivalently results in a quantum operation RSP\mathcal{E}_{\text{RSP}} which transforms |𝕤0\ket{\mathbb{s}_{0}} to |𝕤\ket{\mathbb{s}}, as described in Eq. (1).

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, RSP\mathcal{E}_{\text{RSP}}, implementing the following state transformation:

ρ𝕤=RSP(ρ𝕤0),\rho_{\mathbb{s}}=\mathcal{E}_{\text{RSP}}(\rho_{\mathbb{s}_{0}}), (1)

where ρ𝕤0=|𝕤0𝕤0|\rho_{\mathbb{s}_{0}}=\ket{\mathbb{s}_{0}}\!\!\bra{\mathbb{s}_{0}} is the initial state before the RSP process, and ρ𝕤=|𝕤𝕤|\rho_{\mathbb{s}}=\ket{\mathbb{s}}\!\!\bra{\mathbb{s}} is the output state after the process occurs. It is worth emphasizing that RSP is nonclassical since RSP\mathcal{E}_{\text{RSP}} performs a unitary quantum-mechanical transformation UU. Moreover, RSP\mathcal{E}_{\text{RSP}} 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 RSP\mathcal{E}_{\text{RSP}}. 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 |𝕤0\ket{\mathbb{s}_{0}} can be represented by ρ𝕤0=(I+𝕤0σ)/2\rho_{\mathbb{s}_{0}}=(I+\vec{\mathbb{s}}_{0}\cdot\vec{\sigma})/2, where II is the identity matrix, 𝕤0=(s01,s02,s03)\vec{\mathbb{s}}_{0}=(s_{01},s_{02},s_{03}) is a real three-dimensional vector (called the Bloch vector Nielsen00 , see Fig. 1) with |𝕤0|=1|\vec{\mathbb{s}}_{0}|=1, and σ=(σ1,σ2,σ3)\vec{\sigma}=(\sigma_{1},\sigma_{2},\sigma_{3}) consists of σ1=X\sigma_{1}=X, σ2=Y\sigma_{2}=Y, and σ3=Z\sigma_{3}=Z Pauli matrices. After a RSP process \mathcal{E}, the remote state of Bob’s qubit is prepared in the state (ρ𝕤0)=[(I)+m=13s0m(σm)]/2\mathcal{E}(\rho_{\mathbb{s}_{0}})=[\mathcal{E}(I)+\sum_{m=1}^{3}s_{0m}\mathcal{E}(\sigma_{m})]/2. Explicitly, the observables transform as (I)=n=01(|nmmn|)\mathcal{E}(I)=\sum_{n=0}^{1}\mathcal{E}(\ket{n}_{mm}\!\!\bra{n}) and (σm)=n=01vnm(|nmmn|)\mathcal{E}(\sigma_{m})=\sum_{n=0}^{1}v_{nm}\mathcal{E}(\ket{n}_{mm}\!\!\bra{n}), where vnm=(1)nv_{nm}=(-1)^{n} and |nm\ket{n}_{m} are the eigenvalues and eigenvectors of σm\sigma_{m}, respectively. As shown above, for a given initial state ρ𝕤0\rho_{\mathbb{s}_{0}}, the remote state (ρ𝕤0)\mathcal{E}(\rho_{\mathbb{s}_{0}}) is determined by the RSP process applied to the eigenstates, (|nmmn|)\mathcal{E}(\ket{n}_{mm}\!\!\bra{n}).

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 |nmmn|\ket{n}_{mm}\!\!\bra{n} 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 σm\sigma_{m} under the assumption of classical realism Brunner14 , and is described by vλ(vn11,vn22,vn33)v_{\lambda}\equiv(v_{n_{1}1},v_{n_{2}2},v_{n_{3}3}), where λ=4n1+2n2+n3+1\lambda=4n_{1}+2n_{2}+n_{3}+1 for n1,n2,n3{0,1}n_{1},n_{2},n_{3}\in\{0,1\} is a unique index determining the set. In a classical RSP, denoted as c\mathcal{E}_{c}, an initial state, say |nmmn|\ket{n}_{mm}\!\!\bra{n}, is interpreted as a pre-existing (realistic) outcome vnmv_{nm} of certain pre-existing states vλv_{\lambda}. For example, the classical states vλv_{\lambda} (for λ=1,2,3,4\lambda=1,2,3,4) possess the same state property of vn11=v01v_{n_{1}1}=v_{01}. Moreover, c\mathcal{E}_{c} causes the vλv_{\lambda} to classically transit to another classical pre-existing state. After the state transition has occurred, the final classical state, say vμv_{\mu}, indicates a possible remote state ρμ\rho_{\mu} for Bob’s “qubit”.

The RSP in the above scenario can be described by the map

c(|nmmn|)=λp(λ|vnm)ρ~λn,m,\mathcal{E}_{c}(\ket{n}_{mm}\!\!\bra{n})=\sum_{\lambda}p(\lambda|v_{nm})\;\tilde{\rho}_{\lambda}\ \ \ \ \forall\ n,m, (2)

where p(λ|vnm)p(\lambda|v_{nm}) is the probability of finding vλv_{\lambda} conditioned on the pre-existing outcome vnmv_{nm} with λp(λ|vnm)=1\sum_{\lambda}p(\lambda|v_{nm})=1, ρ~λ=μΩμλρμ\tilde{\rho}_{\lambda}=\sum_{\mu}\Omega_{\mu\lambda}\rho_{\mu} denotes Bob’s possible remote states, and Ωμλ\Omega_{\mu\lambda} is the probability of transition from vλv_{\lambda} to vμv_{\mu}. The process c\mathcal{E}_{c} then classically prepares an arbitrary input state ρ𝕤0\rho_{\mathbb{s}_{0}} as c(ρ𝕤0)=ρ𝕣c|𝕤0\mathcal{E}_{c}(\rho_{\mathbb{s}_{0}})=\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}, where

ρ𝕣c|𝕤0=12[c(I)+m,ns0mvnmc(|nmmn|)]|𝕤0,\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}\!\!=\!\frac{1}{2}[\mathcal{E}_{c}(I)+\!\sum_{m,n}s_{0m}v_{nm}\mathcal{E}_{c}(\ket{n}_{mm}\!\!\bra{n})]\ \ \forall\ket{\mathbb{s}_{0}}, (3)

with a Bloch vector 𝕣c\vec{\mathbb{r}}_{c}. 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 UU, 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 |𝕤0\ket{\mathbb{s}_{0}} in a RSP task. Any \mathcal{E} 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.

Refer to caption
Figure 2: Experimental set-up for RSP. Polarization-correlated photon pairs are created via type-II spontaneous parametric down-conversion (SPDC) Kwiat95 from a β\beta-barium borate (BBO) crystal with 2 mm thickness at a wavelength of 780 nm, where the BBO crystal is pumped with a laser beam at 390 nm with a repetition rate of 76 MHz. The longitudinal and spatial walk-off of the photons in modes aa and bb are compensated by a half-wave plate (HWP) at 9090^{\circ} and a correction BBO (CBBO) of 1 mm thickness. A combination of two quarter-wave plates (QWP) and one HWP is used to control the phase introduced by optical elements such that the entangled photons are prepared in the state |Ψ\ket{\Psi^{-}}. All photons are filtered by narrow bandwidth filters (Δλ3\Delta\lambda\sim 3 nm) and are monitored by silicon avalanche single-photon detectors. Coincidences are recorded by a field-programmable gate-array-based coincidence unit.
Refer to caption
Figure 3: Visualization of nonclassical RSPs. We experimentally demonstrated nonclassical RSPs, ϕ\mathcal{E}_{\phi}, for (a) ϕ=0\phi=0, (b) ϕ=π/2\phi=\pi/2, (c) ϕ=π\phi=\pi, and (d) ϕ=3π/2\phi=3\pi/2. The errors from Poissonian noise in all fidelities are below 0.5%0.5\%. The quantumness of ϕ\mathcal{E}_{\phi} can be visualized from the output states ρ𝕣|𝕤0\rho_{\mathbb{r}|\mathbb{s}_{0}} shown on the Bloch sphere. The ρ𝕣|𝕤0\rho_{\mathbb{r}|\mathbb{s}_{0}} generated by the ideal preparation RSP(ϕ)\mathcal{E}^{(\phi)}_{\text{RSP}} are on the surface of the Bloch sphere (gray mesh). Conditioned on 9898 randomly chosen pure input states ρ𝕤0\rho_{\mathbb{s}_{0}}, we observed all the corresponding output states (blue points) of ϕ\mathcal{E}_{\phi} having larger Bloch vectors, |𝕣|>|𝕣c||\vec{\mathbb{r}}|>|\vec{\mathbb{r}}_{c}|, compared to the states ρ𝕣c|𝕤0\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}} (red points) derived from the classical RSP c\mathcal{E}_{c} with the best process fidelity.

To evaluate how well a given experimental process \mathcal{E} surpasses the generalized LHS model, we compare its output, (ρ𝕤0)=ρ𝕣|𝕤0\mathcal{E}(\rho_{\mathbb{s}_{0}})=\rho_{\mathbb{r}|\mathbb{s}_{0}}, with all possible outputs from classical processes, ρ𝕣c|𝕤0\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}, for all pure input states ρ𝕤0\rho_{\mathbb{s}_{0}}. We illustrate four different methods to measure and (or) identify the capability of \mathcal{E} to exhibit quantum process steering, called process steerability.

(i) Quantum composition α\alpha:

ρ𝕣|𝕤0=αρ𝕣Q|𝕤0+(1α)ρ𝕣c|𝕤0|𝕤0,\rho_{\mathbb{r}|\mathbb{s}_{0}}=\alpha\rho_{\mathbb{r}_{Q}|\mathbb{s}_{0}}+(1-\alpha)\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}\ \ \ \ \forall\ket{\mathbb{s}_{0}}, (4)

where α\alpha describes the minimum proportion of ρ𝕣|𝕤0\rho_{\mathbb{r}|\mathbb{s}_{0}} as being the ρ𝕣Q|𝕤0\rho_{\mathbb{r}_{Q}|\mathbb{s}_{0}} that cannot be represented by ρ𝕣c|𝕤0\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}.

(ii) Quantum robustness β\beta:

ρ𝕣|𝕤0+βρnoise|𝕤01+β=ρ𝕣c|𝕤0|𝕤0,\frac{\rho_{\mathbb{r}|\mathbb{s}_{0}}+\beta\rho_{\text{noise}|\mathbb{s}_{0}}}{1+\beta}=\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}\ \ \ \ \forall\ket{\mathbb{s}_{0}}, (5)

where β\beta is the minimum amount of the noise ρnoise|𝕤0\rho_{\text{noise}|\mathbb{s}_{0}} added such that the ρ𝕣|𝕤0\rho_{\mathbb{r}|\mathbb{s}_{0}} becomes classical.

(iii) Average-state-fidelity criterion for quantum RSP:

F¯s()𝑑𝕤0𝕤0|U(ρ𝕤0)U|𝕤0>F¯sc 0.79,\bar{F}_{s}(\mathcal{E})\!\equiv\!\int\!\!d\mathbb{s}_{0}\!\bra{\mathbb{s}_{0}}\!U^{{\dagger}}\!\mathcal{E}(\rho_{\mathbb{s}_{0}})U\!\ket{\mathbb{s}_{0}}\;>\;\!\bar{F}_{sc}\;\sim\;0.79, (6)

where F¯scmaxcF¯s(c)\bar{F}_{sc}\!\equiv\!\max_{\mathcal{E}_{c}}\bar{F}_{s}(\mathcal{E}_{c}), and the integral is over the uniform measure d𝕤0d\mathbb{s}_{0} on state space and 𝑑𝕤0=1\int d\mathbb{s}_{0}=1.

Refer to caption
Figure 4: Experimental transition from classical to quantum RSP. (a) Experimental RSPs varying with the photon walk-off effect. Given a target RSP, RSP(0)\mathcal{E}^{(0)}_{\text{RSP}}, the process steering of experimental RSPs changes with the setting angle (θHWP\theta_{\text{HWP}}) of the half-wave plate used for compensating the photon walk-off effect (Fig. 2). The walk-off compensation does not exist at θHWP=0\theta_{\text{HWP}}=0^{\circ}, which results in a completely mixed state of photon pair. The experimental average state fidelity, F¯s,expt()\bar{F}_{s,\text{expt}}(\mathcal{E}), and process steerability, αexpt\alpha_{\text{expt}} and βexpt\beta_{\text{expt}}, monotonically increase with θHWP\theta_{\text{HWP}} and are highly consistent with the theoretical predictions of F¯s,theory()\bar{F}_{s,\text{theory}}(\mathcal{E}), αtheory\alpha_{\text{theory}}, and βtheory\beta_{\text{theory}}, where the effect of noise in the photon pairs created has been taken into account SuppMaterial . Compared to the theoretical predictions that are calculated using effective noise states, the ideal predictions of F¯s,ideal()\bar{F}_{s,\text{ideal}}(\mathcal{E}), αideal\alpha_{\text{ideal}}, and βideal\beta_{\text{ideal}} are obtained from perfect photon pairs, |ΨΨ|\ket{\Psi^{-}}\!\!\bra{\Psi^{-}} SuppMaterial . A RSP is nonclassical and close to RSP(0)\mathcal{E}^{(0)}_{\text{RSP}} only when satisfying the criterion [Eq. (18)]. For example, the created state at θHWP=40\theta_{\text{HWP}}=40^{\circ} possesses quantum discord Dakic10 (𝒟theory0.07\mathcal{D}_{\text{theory}}\sim 0.07, 𝒟expt0.08\mathcal{D}_{\text{expt}}\sim 0.08) SuppMaterial and steerability in terms of steerable weight Skrzypczyk2014 (SWtheory0.15SW_{\text{theory}}\sim 0.15, SWexpt0.02SW_{\text{expt}}\sim 0.02) SuppMaterial , but the resulting RSP is not superior to c\mathcal{E}_{c}. (b) Experimental RSPs by the Werner states. A variety of polarization-correlated photon pairs in the Werner states ρW(pnoise)\rho_{W}(p_{\text{noise}}) are created to realize the RSP protocol. As theoretically shown and experimentally demonstrated, a RSP is nonclassical only when the noise is lower than certain thresholds. Moreover, a ρW(pnoise)\rho_{W}(p_{\text{noise}}) which possesses quantum discord can even result in a classical RSP, e.g., (𝒟theory0.11\mathcal{D}_{\text{theory}}\sim 0.11, 𝒟expt0.10\mathcal{D}_{\text{expt}}\sim 0.10) for pnoise=0.5p_{\text{noise}}=0.5, where the effect of noise in created photon pairs is also taken into account SuppMaterial . The ideal values of F¯s,ideal()\bar{F}_{s,\text{ideal}}(\mathcal{E}), αideal\alpha_{\text{ideal}}, and βideal\beta_{\text{ideal}} are calculated using ρW(pnoise)\rho_{W}(p_{\text{noise}}).

(iv) Complementary-state-fidelity criterion for quantum RSP:

F¯comp()14q=12n=01Fsn(q)()>14(3F¯sc+1)0.84,\bar{F}_{\text{comp}}(\mathcal{E})\!\!\equiv\!\!\frac{1}{4}\!\sum_{q=1}^{2}\sum_{n=0}^{1}\!F^{(q)}_{sn}(\mathcal{E})\!\;>\;\!\frac{1}{4}(3\bar{F}_{sc}+1)\sim 0.84, (7)

where Fsn(q)()𝕤0n(q)|U(ρ𝕤0n(q))U|𝕤0n(q)F^{(q)}_{sn}(\mathcal{E})\equiv\bra{\mathbb{s}^{(q)}_{0n}}\!U^{{\dagger}}\!\mathcal{E}(\rho_{\mathbb{s}^{(q)}_{0n}})U\!\ket{\mathbb{s}^{(q)}_{0n}} for q=1,2q=1,2 denote the state fidelities over arbitrary two orthonormal and complementary state bases, S1={|𝕤00(1),|𝕤01(1)}S_{1}=\{\ket{\mathbb{s}^{(1)}_{00}},\ket{\mathbb{s}^{(1)}_{01}}\} and S2={|𝕤00(2),{|𝕤01(2)}S_{2}=\{\ket{\mathbb{s}^{(2)}_{00}},\{\ket{\mathbb{s}^{(2)}_{01}}\}, with |𝕤0n(1)|𝕤0n(2)|=1/2|\langle\mathbb{s}^{(1)}_{0n}|\mathbb{s}^{(2)}_{0n^{\prime}}\rangle|=1/\sqrt{2} for n,n=0,1n,n^{\prime}=0,1.

In measuring α\alpha or β\beta, or determining the threshold F¯sc\bar{F}_{sc}, the optimization problem with respect to the dynamical maps c\mathcal{E}_{c} can be solved by semi-definite programming SuppMaterial . The criterion [Eq. (19)] is based on the fact that, since F¯comp()(3F¯s()+1)/4\bar{F}_{\text{comp}}(\mathcal{E})\leq(3\bar{F}_{s}(\mathcal{E})+1)/4 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 |𝕤0\ket{\mathbb{s}_{0}}, 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 RSP\mathcal{E}_{\text{RSP}} of U=R(ϕ)U=R(\phi), one can follow the RSP protocol and examine the experimental results with the criterion [Eq. (19)], using |𝕤0Sn={|0n+1,|1n+1}\ket{\mathbb{s}_{0}}\in S_{n}\!=\!\{\ket{0}_{n+1},\ket{1}_{n+1}\}, for n=1,2n=1,2, as input states.

We now consider the process steerability of the ideal RSP process for arbitrary UU. In this ideal case, RSP\mathcal{E}_{\text{RSP}} has a quantum composition of α=1\alpha=1, a quantum robustness of β=0.464\beta=0.464, and for the fidelity criteria [Eqs. (18) and (19)], gives F¯s(RSP)=F¯comp(RSP)=1\bar{F}_{s}(\mathcal{E}_{\text{RSP}})=\bar{F}_{\text{comp}}(\mathcal{E}_{\text{RSP}})=1. These values will decrease when contaminated with noise. For example, when the Werner state ρW(pnoise)=(1pnoise)|ΨΨ|+pnoiseII/4\rho_{W}(p_{\text{noise}})=(1-p_{\text{noise}})\ket{\Psi^{-}}\!\!\bra{\Psi^{-}}+p_{\text{noise}}I\otimes I/4, is used for RSP to represent noise in the state creation process, the process steerability decreases with increasing noise intensity pnoisep_{\text{noise}} and disappears at pnoise=11/30.42p_{\text{noise}}=1-1/\sqrt{3}\sim 0.42, 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, 𝕤0|U(ρ𝕤0)U|𝕤0=1/2\bra{\mathbb{s}_{0}}\!U^{{\dagger}}\!\mathcal{E}(\rho_{\mathbb{s}_{0}})U\!\ket{\mathbb{s}_{0}}=1/2, 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 ϕ\mathcal{E}_{\phi}, while we denote the ideal RSP process as RSP(ϕ)\mathcal{E}^{(\phi)}_{\text{RSP}}, which are identified by the choice of single-qubit unitaries, U=R(ϕ)U=R(\phi). The experimental results for the criteria outlined in (i)-(iv), and the concrete comparisons between ϕ\mathcal{E}_{\phi}, RSP(ϕ)\mathcal{E}^{(\phi)}_{\text{RSP}} and c\mathcal{E}_{c} 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 aa and bb (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 ρW(pnoise)\rho_{W}(p_{\text{noise}}). To experimentally set the noise weight of pnoisep_{\text{noise}}, 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 c\mathcal{E}_{c} 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 α\alpha, β\beta, and F¯sc\bar{F}_{sc}. 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. Du¨\ddot{\text{u}}r, 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:

ρnm=λp(λ|vnm)ρλn,m,\rho_{nm}=\sum_{\lambda}p(\lambda|v_{nm})\rho_{\lambda}\ \ \ \ \ \ \ \forall n,m, (8)

where vnmv_{nm} denotes the nnth measurement outcome of Alice’s mmth measurement. Alternatively, this can be represented by unnormalized states and rephrased as

ρnm=λp(vnm|λ)ρλn,m,\rho^{\prime}_{nm}=\sum_{\lambda}p(v_{nm}|\lambda)\rho^{\prime}_{\lambda}\ \ \ \ \ \ \ \forall n,m, (9)

where ρnm=p(vnm)ρnm\rho^{\prime}_{nm}=p(v_{nm})\rho_{nm} and ρλ=p(λ)ρλ\rho^{\prime}_{\lambda}=p(\lambda)\rho_{\lambda}, and p(vnm|λ)p(v_{nm}|\lambda) is the probability of finding the outcome vnmv_{nm} conditioned on the classical state vλv_{\lambda}. The relation of p(vnm)p(λ|vnm)=p(λ)p(vnm|λ)p(v_{nm})p(\lambda|v_{nm})=p(\lambda)p(v_{nm}|\lambda) has be used to derive Eq. (9) from Eq. (8). Given Alice’s choice of measurement mm and outcome vnmv_{nm}, any assemblage {ρnm}nm\{\rho^{\prime}_{nm}\}_{nm} 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 |nmmn|\ket{n}_{mm}\!\!\bra{n}, a classical RSP, c\mathcal{E}_{c}, generates an output described by Eq. (2) in the main text:

c(|nmmn|)=λμp(λ|vnm)Ωμλρμn,m,\mathcal{E}_{c}(\ket{n}_{mm}\!\!\bra{n})=\sum_{\lambda}\sum_{\mu}p(\lambda|v_{nm})\Omega_{\mu\lambda}\rho_{\mu}\ \ \ \ \forall\ n,m, (10)

which can be written as

c(|nmmn|)=λp(λ|vnm)ρλn,m,\mathcal{E}_{c}(\ket{n}_{mm}\!\!\bra{n})=\sum_{\lambda}p(\lambda|v_{nm})\rho_{\lambda}\ \ \ \ \forall\ n,m, (11)

under the condition of Ωμλ=δμλ\Omega_{\mu\lambda}=\delta_{\mu\lambda}. 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, c\mathcal{E}_{c} cannot be fully defined by considering only some of the output states c(|nmmn|)\mathcal{E}_{c}(\ket{n}_{mm}\!\!\bra{n}) under the inputs |nm\ket{n}_{m} in the LHS model via Eq. (11). Rather, since the c\mathcal{E}_{c} 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 {ρ𝕣|𝕤0}\{\rho_{\mathbb{r}|\mathbb{s}_{0}}\} of output states derived from an experimental RSP process which cannot described by

ρ𝕣c|𝕤0=12[c(I)+m,ns0mvnmc(|nmmn|)]|𝕤0,\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}\!\!=\!\frac{1}{2}[\mathcal{E}_{c}(I)+\!\sum_{m,n}s_{0m}v_{nm}\mathcal{E}_{c}(\ket{n}_{mm}\!\!\bra{n})]\ \ \forall\ket{\mathbb{s}_{0}}, (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 c\mathcal{E}_{c} 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 t=0t=0. After the measurement, the initial state becomes, say |nmmn|\ket{n}_{mm}\!\!\bra{n}. Then this state is sent into a quantum channel \mathcal{E} for a time tt. At time tt, Bob receives the system of the state (|nmmn|)=ρnm(t)\mathcal{E}(\ket{n}_{mm}\!\!\bra{n})=\rho_{nm}(t). Following the un-normalized assemblage Wiseman07 ; Reid09 in the standard LHS model, the un-normalized states in temporal scenario is defined by ρnmT=p(vnm)ρnm(t)\rho^{T}_{nm}=p(v_{nm})\rho_{nm}(t), where the superscript TT reminds one that the assemblage {ρnmT}nm\{\rho^{T}_{nm}\}_{nm} is for temporal steering. In the temporal LHS model, an unsteerable assemblage is defined as

ρnmT=λp(vnm|λ)ρλTn,m,\rho^{T}_{nm}=\sum_{\lambda}p(v_{nm}|\lambda)\rho^{T}_{\lambda}\ \ \ \ \ \ \ \forall n,m, (13)

where ρλT=p(λ)ρλ(t)\rho^{T}_{\lambda}=p(\lambda)\rho_{\lambda}(t) and ρλ(t)\rho_{\lambda}(t) is created from a source which determines the possible correlations between Alice’s (at time t=0t=0) and Bob’s measurement results (at the time tt) under classical realism. By using the same method as shown above, Eq. (13) can be rephrased as

ρnm(t)=λp(λ|vnm)ρλ(t)n,m.\rho_{nm}(t)=\sum_{\lambda}p(\lambda|v_{nm})\rho_{\lambda}(t)\ \ \ \ \ \ \ \forall n,m. (14)

Still, when comparing Eq. (11) to Eq. (14), c\mathcal{E}_{c} cannot be fully defined by considering the finite output states c(|nmmn|)\mathcal{E}_{c}(\ket{n}_{mm}\!\!\bra{n}) under the inputs |nmm\ket{n}_{mm} in the temporal LHS model via Eq. (11). The c\mathcal{E}_{c} can be defined only when the output states, c(ρ𝕤0)=ρ𝕣c|𝕤0(t)\mathcal{E}_{c}(\rho_{\mathbb{s}_{0}})=\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}(t), for all the pure input states are considered, where

ρ𝕣c|𝕤0(t)=12[c(I)+m,n,λs0mvnmc(|nmmn|)]|𝕤0.\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}(t)\!\!=\!\frac{1}{2}[\mathcal{E}_{c}(I)+\!\!\!\sum_{m,n,\lambda}\!\!\!s_{0m}v_{nm}\mathcal{E}_{c}(\ket{n}_{mm}\!\!\bra{n})]\ \ \forall\ket{\mathbb{s}_{0}}. (15)

II Quantum composition, quantum robustness, and the fidelity criteria

In this section, we will introduce how to calculate quantum composition α\alpha, quantum robustness β\beta, and the fidelity bound F¯sc\bar{F}_{sc} in the fidelity criteria. The methods introduced in the main text to measure and (or) identify process steerability of an experimental process \mathcal{E} are summarized as follows:

(i) Quantum composition α\alpha. For all possible pure states |𝕤0\ket{\mathbb{s}_{0}} (i.e., |𝕤0|=1|\vec{\mathbb{s}}_{0}|=1), the corresponding prepared remote states, (ρ𝕤0)=ρ𝕣|𝕤0\mathcal{E}(\rho_{\mathbb{s}_{0}})=\rho_{\mathbb{r}|\mathbb{s}_{0}}, can be decomposed as

ρ𝕣|𝕤0=αρ𝕣Q|𝕤0+(1α)ρ𝕣c|𝕤0|𝕤0,\rho_{\mathbb{r}|\mathbb{s}_{0}}=\alpha\rho_{\mathbb{r}_{Q}|\mathbb{s}_{0}}+(1-\alpha)\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}\ \ \ \ \forall\ket{\mathbb{s}_{0}}, (16)

where the quantum composition α\alpha describes the minimum proportion of ρ𝕣|𝕤0\rho_{\mathbb{r}|\mathbb{s}_{0}} that cannot be represented by the classically prepared states ρ𝕣c|𝕤0\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}. An ideal RSP shows that ρ𝕣|𝕤0\rho_{\mathbb{r}|\mathbb{s}_{0}} is fully composed of the ρ𝕣Q|𝕤0\rho_{\mathbb{r}_{Q}|\mathbb{s}_{0}} that cannot be represented by ρ𝕣c|𝕤0\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}, i.e., α=1\alpha=1.

(ii) Quantum robustness β\beta. In this case the prepared states can be written as

ρ𝕣|𝕤0+βρnoise|𝕤01+β=ρ𝕣c|𝕤0|𝕤0,\frac{\rho_{\mathbb{r}|\mathbb{s}_{0}}+\beta\rho_{\text{noise}|\mathbb{s}_{0}}}{1+\beta}=\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}\ \ \ \ \forall\ket{\mathbb{s}_{0}}, (17)

where the quantum robustness β\beta is the minimum amount of the noise ρnoise|𝕤0\rho_{\text{noise}|\mathbb{s}_{0}} added such that the ρ𝕣|𝕤0\rho_{\mathbb{r}|\mathbb{s}_{0}} becomes a classically prepared remote state. When β=0\beta=0, the experimental RSP is already classical. Increasing β\beta increases the difference between ρ𝕣|𝕤0\rho_{\mathbb{r}|\mathbb{s}_{0}} and ρ𝕣c|𝕤0\rho_{\mathbb{r}_{c}|\mathbb{s}_{0}}.

(iii) Average state fidelity. The \mathcal{E} is identified as being faithful to show quantum process steering if the state fidelity of (ρ𝕤0)\mathcal{E}(\rho_{\mathbb{s}_{0}}) and RSP(ρ𝕤0)\mathcal{E}_{\text{RSP}}(\rho_{\mathbb{s}_{0}}): Fs,𝕤0()𝕤0|U(ρ𝕤0)U|𝕤0F_{s,\mathbb{s}_{0}}(\mathcal{E})\equiv\bra{\mathbb{s}_{0}}\!U^{{\dagger}}\!\mathcal{E}(\rho_{\mathbb{s}_{0}})U\!\ket{\mathbb{s}_{0}}, averaged over all possible input pure states is larger than the maximum averaged state fidelity derived from classical RSP, F¯scmaxcF¯s(c)\bar{F}_{sc}\!\equiv\!\max_{\mathcal{E}_{c}}\bar{F}_{s}(\mathcal{E}_{c}), i.e.,

F¯s()𝑑𝕤0Fs,𝕤0()>F¯sc0.79,\bar{F}_{s}(\mathcal{E})\!\equiv\!\int\!d\mathbb{s}_{0}F_{s,\mathbb{s}_{0}}(\mathcal{E})>\!\bar{F}_{sc}\sim 0.79, (18)

where the integral is over the uniform measure d𝕤0d\mathbb{s}_{0} on state space and 𝑑𝕤0=1\int d\mathbb{s}_{0}=1.

(iv) Complementary state fidelity. Suppose that S1={|𝕤00(1),|𝕤01(1)}S_{1}=\{\ket{\mathbb{s}^{(1)}_{00}},\ket{\mathbb{s}^{(1)}_{01}}\} and S2={|𝕤00(2),{|𝕤01(2)}S_{2}=\{\ket{\mathbb{s}^{(2)}_{00}},\{\ket{\mathbb{s}^{(2)}_{01}}\} are two orthonormal and complementary state bases, where |𝕤0n(1)|𝕤0n(2)|=1/2|\langle\mathbb{s}^{(1)}_{0n}|\mathbb{s}^{(2)}_{0n^{\prime}}\rangle|=1/\sqrt{2} for n,n=0,1n,n^{\prime}=0,1. Let us consider the state fidelities Fsn(q)()𝕤0n(q)|U(ρ𝕤0n(q))U|𝕤0n(q)F^{(q)}_{sn}(\mathcal{E})\equiv\bra{\mathbb{s}^{(q)}_{0n}}\!U^{{\dagger}}\!\mathcal{E}(\rho_{\mathbb{s}^{(q)}_{0n}})U\!\ket{\mathbb{s}^{(q)}_{0n}} for q=1,2q=1,2. The \mathcal{E} is faithful and possess quantum process steerability if the complementary state fidelity F¯comp()\bar{F}_{\text{comp}}(\mathcal{E}) Hofmann05 beating the classical limit implies quantum RSP:

F¯comp()14q=12n=01Fsn(q)()>F¯comp,c14(3F¯sc+1)0.84.\bar{F}_{\text{comp}}(\mathcal{E})\!\!\equiv\!\!\frac{1}{4}\!\sum_{q=1}^{2}\sum_{n=0}^{1}\!F^{(q)}_{sn}(\mathcal{E})\!>\!\bar{F}_{\text{comp},c}\equiv\frac{1}{4}(3\bar{F}_{sc}+1)\sim 0.84. (19)

In measuring the α\alpha, β\beta, F¯s\bar{F}_{s} and determining the threshold F¯sc\bar{F}_{sc} and F¯comp,c\bar{F}_{\text{comp},c}, it is necessary to consider the contributions from all pure states |𝕤0\ket{\mathbb{s}_{0}} in the respective computational tasks. To solve such questions, we first note that Eq. (16) and Eq. (17) can be rephrased as

(ρ𝕤0)=αQ(ρ𝕤0)+(1α)c(ρ𝕤0)\mathcal{E}(\rho_{\mathbb{s}_{0}})=\alpha\mathcal{E}_{Q}(\rho_{\mathbb{s}_{0}})+(1-\alpha)\mathcal{E}_{c}(\rho_{\mathbb{s}_{0}}) (20)

and

(ρ𝕤0)+βnoise(ρ𝕤0)1+β=c(ρ𝕤0),\frac{\mathcal{E}(\rho_{\mathbb{s}_{0}})+\beta\mathcal{E}_{\text{noise}}(\rho_{\mathbb{s}_{0}})}{1+\beta}=\mathcal{E}_{c}(\rho_{\mathbb{s}_{0}}), (21)

respectively, where Q(ρ𝕤0)=ρ𝕣Q|𝕤0\mathcal{E}_{Q}(\rho_{\mathbb{s}_{0}})=\rho_{\mathbb{r}_{Q}|\mathbb{s}_{0}} and noise(ρ𝕤0)=ρnoise|𝕤0\mathcal{E}_{\text{noise}}(\rho_{\mathbb{s}_{0}})=\rho_{\text{noise}|\mathbb{s}_{0}}. Then the optimization problems with respect to states become to the optimization of α\alpha and β\beta with respect to dynamical maps Hsieh17 . Moreover, it has been shown that F¯s()=(2F+1)/3\bar{F}_{s}(\mathcal{E})=(2F_{\mathcal{E}}+1)/3 Gilchrist05 ; Hofmann05 , implying that the F¯s\bar{F}_{s} can be derived from FF_{\mathcal{E}}, i.e., the process fidelity of \mathcal{E} and RSP\mathcal{E}_{\text{RSP}}. Equivalently, the criterion (18) becomes F>FcF_{\mathcal{E}}>F_{\mathcal{E}_{c}}, where the FcF_{\mathcal{E}_{c}} is obtained by considering the best similarity between c\mathcal{E}_{c} and RSP\mathcal{E}_{\text{RSP}}. Finally, as \mathcal{E} is experimentally measured, α\alpha, and β\beta, F¯sc\bar{F}_{sc} can therefore be determined by semi-definite programming (SDP) Lofberg ; sdpsolver and the F¯comp,c\bar{F}_{\text{comp},c} can be calculated through F¯sc\bar{F}_{sc}.

In practical experiments, the experimental process \mathcal{E} can be fully described by using the process tomography algorithm Nielsen00 ; Chuang97 and then be represented by a positive Hermitian matrix, called process matrix, χexpt\chi_{\text{expt}}. The process matrix χexpt\chi_{\text{expt}} is constructed by the specific output states (|nmmn|)\mathcal{E}(\ket{n}_{mm}\!\!\bra{n}) which are the output states of the eigenstates of three Pauli matrices. Considering the output states (|nmmn|)\mathcal{E}(\ket{n}_{mm}\!\!\bra{n}) from a classical RSP c\mathcal{E}_{c}, we can use process tomography algorithm to get a process matrix χc\chi_{c} for the classical RSP c\mathcal{E}_{c}. Thus, Eq. (20) and Eq. (21) can be rephrased as

χexpt=αχQ+(1α)χc\chi_{\text{expt}}=\alpha\chi_{Q}+(1-\alpha)\chi_{c} (22)

and

χexpt+βχnoise1+β=χc,\frac{\chi_{\text{expt}}+\beta\chi_{\text{noise}}}{1+\beta}=\chi_{c}, (23)

respectively, and the FcF_{\mathcal{E}_{c}} can be represented as Fc=maxχc[tr(χcχRSP)]F_{\mathcal{E}_{c}}=\max_{\chi_{c}}[\text{tr}(\chi_{c}\chi_{\text{RSP}})], where χRSP\chi_{\text{RSP}} is the process matrix of RSP\mathcal{E}_{\text{RSP}} and χQ\chi_{Q} and χnoise\chi_{\text{noise}} are the process matrices of Q\mathcal{E}_{Q} and noise\mathcal{E}_{\text{noise}} respectively. The process matrices χexpt\chi_{\text{expt}} and χc\chi_{c} can be used to solve the optimization problems of the optimization of α\alpha, β\beta, and F¯sc\bar{F}_{sc}. To calculate α\alpha, β\beta, and F¯sc\bar{F}_{sc} via SDP, we let the processes in SDP satisfy the definition of process matrices or density matrices by the following constraint. A process matrix χ\chi and a density matrix ρ\rho must be positive semi-definite, i.e., χ0\chi\geq 0 and ρ0\rho\geq 0. Since a process matrix is constructed of the output states (|nmmn|)\mathcal{E}(\ket{n}_{mm}\!\!\bra{n}) and the output states of the c\mathcal{E}_{c} for classical RSP model is determined by Bob’s possible remote states ρ~λ\tilde{\rho}_{\lambda} in Eq. (2) in the main text, the process matrix χc\chi_{c} must satisfy ρ~λ0,λ\tilde{\rho}_{\lambda}\geq 0,\forall\lambda.

According to Eq. (22), by representing (1α)χc(1-\alpha)\chi_{c} as a unnormalized process matrix χ~c=(1α)χc\tilde{\chi}_{c}=(1-\alpha)\chi_{c}, quantum composition α\alpha can be obtained by minimizing the following quantity via SDP with MATLAB Lofberg ; sdpsolver : α=minχ~c[1tr(χ~c)]\alpha=\min_{\tilde{\chi}_{c}}[1-\text{tr}(\tilde{\chi}_{c})], under the following constraints:

χ~c0;ρ~λ0,λ;χexptχ~c0.\begin{split}\tilde{\chi}_{c}\geq 0;\\ \tilde{\rho}_{\lambda}\geq 0,\forall\lambda;\\ \chi_{\text{expt}}-\tilde{\chi}_{c}\geq 0.\\ \end{split}

The first constraint ensures that the χ~c\tilde{\chi}_{c} 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 Q\mathcal{E}_{Q} satisfy the definition of process matrices which need to be positive semi-definite.

As the method used to solve α\alpha, to calculate the quantum robustness β\beta, we define χ~c=(1+β)χc\tilde{\chi}_{c}=(1+\beta)\chi_{c}. Thus, according to Eq. (23), the quantum robustness β\beta of χexpt\chi_{\text{expt}} can be obtained by using SDP to solve β=minχ~c[tr(χ~c)1]\beta=\min_{\tilde{\chi}_{c}}[\text{tr}(\tilde{\chi}_{c})-1]. The constraints for solving β\beta via SDP are

χ~c0;ρ~λ0,λ;χ~cχexpt0;tr(χ~c)1.\begin{split}\tilde{\chi}_{c}\geq 0;\\ \tilde{\rho}_{\lambda}\geq 0,\forall\lambda;\\ \tilde{\chi}_{c}-\chi_{\text{expt}}\geq 0;\\ \text{tr}(\tilde{\chi}_{c})\geq 1.\end{split}

The first and second constraints are the same as the constraints used in calculating α\alpha. The third constraint ensures that β0\beta\geq 0. The fourth constraint is the condition of positive semidefiniteness for βχnoise=(1+β)χcχexpt\beta\chi_{\text{noise}}=(1+\beta)\chi_{c}-\chi_{\text{expt}} in Eq. (23).

With the process matrix χRSP\chi_{\text{RSP}} of RSP\mathcal{E}_{\text{RSP}}, the fidelity bounds Fcmaxχ~c[tr(χ~cχRSP)]F_{\mathcal{E}_{c}}\equiv\max_{\tilde{\chi}_{c}}[\text{tr}(\tilde{\chi}_{c}\chi_{\text{RSP}})] can be solved under the constraints

χ~c0;ρ~λ0,λ;tr(χ~c)=1,\begin{split}\tilde{\chi}_{c}\geq 0;\\ \tilde{\rho}_{\lambda}\geq 0,\forall\lambda;\\ \text{tr}(\tilde{\chi}_{c})=1,\end{split}

which make the χ~c\tilde{\chi}_{c} be a normalized matrix to calculate fidelity. While we get the fidelity bound FcF_{\mathcal{E}_{c}}, the average-state-fidelity bound F¯sc\bar{F}_{sc} and the complementary-state-fidelity bound F¯comp,c\bar{F}_{\text{comp},c} can be calculated through F¯sc=(2Fc+1)/3\bar{F}_{sc}=(2F_{\mathcal{E}_{c}}+1)/3 and F¯comp,c=(3F¯sc+1)/4\bar{F}_{\text{comp},c}=(3\bar{F}_{sc}+1)/4.

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 β\beta-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 aa and mode bb are 9090^{\circ} under ideal compensation circumstance, i.e., θHWP,a=θHWP,b=θHWP=90\theta_{\text{HWP},a}=\theta_{\text{HWP},b}=\theta_{\text{HWP}}=90^{\circ}. Subsequently, we utilize a phase controller consisting of two quarter-wave plates (QWP) and a HWP to obtain the desired state, |Ψ\ket{\Psi^{-}}. 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 ρexpt\rho_{\text{expt}}, 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 aa and mode bb, respectively. The measured density matrix of the created photon pairs after the designed compensation (θHWP=90\theta_{\text{HWP}}=90^{\circ}) is of the following form:

ρexpt(90)=[0.020.01+0.02i0.010.01i0.01i0.010.02i0.490.45+0.03i0.01+0.01i0.01+0.01i0.450.03i0.480.010.02i0.01i0.010.01i0.01+0.02i0.01],\rho_{\text{expt}}(90^{\circ})=\left[\begin{matrix}0.02&-0.01+0.02i&-0.01-0.01i&0.01i\\ -0.01-0.02i&0.49&-0.45+0.03i&0.01+0.01i\\ -0.01+0.01i&-0.45-0.03i&0.48&0.01-0.02i\\ -0.01i&0.01-0.01i&0.01+0.02i&0.01\\ \end{matrix}\right], (24)

where the fidelity of ρexpt(90)\rho_{\text{expt}}(90^{\circ}) and the target state, |ΨΨ|\ket{\Psi^{-}}\!\bra{\Psi^{-}}, is F=tr[ρexpt(90)|ΨΨ|]0.94F=\text{tr}\left[\rho_{\text{expt}}(90^{\circ})\ket{\Psi^{-}}\!\bra{\Psi^{-}}\right]\sim 0.94 in our experiment. The maximum-likelihood technique James01 has been used to make the desnsity matrix ρexpt(90)\rho_{\text{expt}}(90^{\circ}) reasonable.

The created entangled photon pairs of the state ρexpt(90)\rho_{\text{expt}}(90^{\circ}) is used to perform nonclassical RSPs, ϕ\mathcal{E}_{\phi}, for (a) ϕ=0\phi=0, (b) ϕ=π/2\phi=\pi/2, (c) ϕ=π\phi=\pi, and (d) ϕ=3π/2\phi=3\pi/2 shown in Fig. 3 in the main text. Moreover, ρexpt(90)\rho_{\text{expt}}(90^{\circ}) is also used to demonstrate 0\mathcal{E}_{0} in the case of the desired compensation (θHWP=90\theta_{\text{HWP}}=90^{\circ}) in Fig. 4(a) and in the case of Werner state at pnoise=0p_{\text{noise}}=0 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 |Ψ\ket{\Psi^{-}} 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 |ΨΨ|\ket{\Psi^{-}}\!\bra{\Psi^{-}} 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

ρnoise(90)=(1p)|ΨΨ|+pII4, 0p1,\rho_{\text{noise}}(90^{\circ})=(1-p)\ket{\Psi^{-}}\!\bra{\Psi^{-}}+p\frac{I\otimes I}{4},\ \ \ \ \ \ \ 0\leq p\leq 1, (25)

where pp denotes the intensity of white noise in created photon pairs. The pp in our experiment can be estmated by maximizing the fidelity between ρnoise(90)\rho_{\text{noise}}(90^{\circ}) and the experimentally created photon pairs, ρexpt(90)\rho_{\text{expt}}(90^{\circ}) (24). The fidelity between them is defined as

=trρnoise(90)12ρexpt(90)ρnoise(90)12.\mathcal{F}=\text{tr}\sqrt{\rho_{\text{noise}}(90^{\circ})^{\frac{1}{2}}\rho_{\text{expt}}(90^{\circ})\rho_{\text{noise}}(90^{\circ})^{\frac{1}{2}}}.

After performing the maximization task, the fidelity between ρnoise(90)\rho_{\text{noise}}(90^{\circ}) and ρexpt(90)\rho_{\text{expt}}(90^{\circ}) is 0.99\mathcal{F}\sim 0.99, and the estimated noise intensity is p=0.06p=0.06. Therefore, the experimental photon pairs in our noise model is

ρnoise(90)=0.94|ΨΨ|+0.06II4.\rho_{\text{noise}}(90^{\circ})=0.94\ket{\Psi^{-}}\!\bra{\Psi^{-}}+0.06\frac{I\otimes I}{4}. (26)

We utilize the effective noise state ρnoise(90)\rho_{\text{noise}}(90^{\circ}) to calculate the theoretical values of αtheory\alpha_{\text{theory}}, βtheory\beta_{\text{theory}}, and F¯s,theory()\bar{F}_{s,\text{theory}}(\mathcal{E}) at θHWP=90\theta_{\text{HWP}}=90^{\circ} 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

ρW,noise(pnoise)=(1pnoise)ρnoise(90)+pnoiseII4,\rho_{W,\text{noise}}(p_{\text{noise}})=(1-p_{\text{noise}})\rho_{\text{noise}}(90^{\circ})+p_{\text{noise}}\frac{I\otimes I}{4}, (27)

where 0pnoise10\leq p_{\text{noise}}\leq 1 is the noise weight in Werner state. We utilize ρW,noise(pnoise)\rho_{W,\text{noise}}(p_{\text{noise}}) to estimate the theoretical values of αtheory\alpha_{\text{theory}}, βtheory\beta_{\text{theory}}, F¯s,theory()\bar{F}_{s,\text{theory}}(\mathcal{E}), 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 ρW,noise(pnoise)\rho_{W,\text{noise}}(p_{\text{noise}}) 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

|Ψ=ηA(ωp)dωpd3kod3keδ(ωpωoωe)0Ldzei(kpkoke)zakoake|vac,\ket{\Psi}=\eta\int\text{A}(\omega_{\text{p}})\text{d}\omega_{\text{p}}\int\text{d}^{3}\textbf{{k}}_{\text{o}}\int\text{d}^{3}\textbf{{k}}_{\text{e}}\delta(\omega_{\text{p}}-\omega_{\text{o}}-\omega_{\text{e}})\int_{0}^{L}\text{d}ze^{i(\textbf{{k}}_{\text{p}}-\textbf{{k}}_{\text{o}}-\textbf{{k}}_{\text{e}})z}a_{\textbf{{k}}_{\text{o}}}^{\dagger}a_{\textbf{{k}}_{\text{e}}}^{\dagger}\ket{\text{vac}}, (28)

where the zz direction is assumed to be parallel to the pump beam, η\eta is the constant comparing with SPDC efficient; ωp\omega_{\text{p}}, ωo\omega_{\text{o}} and ωe\omega_{\text{e}} are the frequencies of the pump, ordinary, and extraordinary beams with the wave vectors kp\textbf{{k}}_{\text{p}}, ko\textbf{{k}}_{\text{o}}, and ke\textbf{{k}}_{\text{e}}, respectively; A(ωp)\text{A}(\omega_{\text{p}}) is the spectral profile of the pump pulse; δ\delta-function δ(ωpωoωe)\delta(\omega_{\text{p}}-\omega_{\text{o}}-\omega_{\text{e}}) shows the energy conservation in SPDC; L=2L=2 mm is the length of the BBO crystal; akoa_{\textbf{{k}}_{\text{o}}}^{\dagger} and akea_{\textbf{{k}}_{\text{e}}}^{\dagger} are the creation operators for the ordinary and extraordinary beams respectively; |vac\ket{\text{vac}} represents the vacuum state. Since the detectors are in the modes aa and bb and narrow bandwidth filters are added in front of the detectors to select frequencies ωo=ωe=ωp/2\omega_{\text{o}}=\omega_{\text{e}}=\omega_{\text{p}}/2, Eq. (28) can be simplified as

|Ψ=η(a0b1+eiϕa1b0)|vac,\ket{\Psi^{\prime}}=\eta^{\prime}(a_{\text{0}}^{\dagger}b_{\text{1}}^{\dagger}+e^{i\phi}a_{\text{1}}^{\dagger}b_{\text{0}}^{\dagger})\ket{\text{vac}},

where η\eta^{\prime} is the normalization constant, and aia_{i}^{\dagger}, bib_{i}^{\dagger} are the creation operators of photon with polarized state |i\ket{i} for spatial mode aa, bb respectively. After the normalization, the polarized state of the generated photon pairs is |ψ=(|01+eiϕ|10)/2\ket{\psi}=(\ket{01}+e^{i\phi}\ket{10})/\sqrt{2}, 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, z=L/2z=L/2, the walk-off effect makes horizontally polarized photons be delayed by a time t=|1/u01/u1|L/2t=\left|1/u_{0}-1/u_{1}\right|L/2 with respect vertically polarized photons, where u0u_{0} and u1u_{1} 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 |ψ\ket{\psi} becomes

|ψ=12(|0t10+eiϕ|100t),\ket{\psi^{\prime}}=\frac{1}{\sqrt{2}}(\ket{0_{t}1_{0}}+e^{i\phi}\ket{1_{0}0_{t}}), (29)

where |ij\ket{i_{j}} denotes the photon with polarized state |i\ket{i} is delayed time jj. The delay time tt makes |0t10\ket{0_{t}1_{0}} and |100t\ket{1_{0}0_{t}} 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

ρwalk-off=12(|0101|+|1010|).\rho_{\text{walk-off}}=\frac{1}{2}(\ket{01}\!\!\bra{01}+\ket{10}\!\!\bra{10}). (30)

To compensate the photon walk-off effect, a HWP and a CBBO crystal are used in each mode aa and bb. See Fig. 2 in the main text. The unitary transform of HWP for different setting angle θHWP\theta_{\text{HWP}} is

UHWP(θHWP)=[cos(θHWP)sin(θHWP)sin(θHWP)cos(θHWP)].U_{\text{HWP}}(\theta_{\text{HWP}})=\left[\begin{matrix}\text{cos}(\theta_{\text{HWP}})&\text{sin}(\theta_{\text{HWP}})\\ \text{sin}(\theta_{\text{HWP}})&-\text{cos}(\theta_{\text{HWP}})\end{matrix}\right].

By setting the angle of the HWP in each mode aa and bb at the same angle θHWP\theta_{\text{HWP}}, i.e., θHWP,a=θHWP,b=θHWP\theta_{\text{HWP,a}}=\theta_{\text{HWP,b}}=\theta_{\text{HWP}}, the output state from the BBO crystal |ψ\ket{\psi^{\prime}} becomes

|ψ′′(θHWP)=12[(cos(θHWP)|0t+sin(θHWP)|1t)(sin(θHWP)|00cos(θHWP)|10)+eiϕ(sin(θHWP)|00cos(θHWP)|10)(cos(θHWP)|0t+sin(θHWP)|1t)].\begin{split}\ket{\psi^{\prime\prime}(\theta_{\text{HWP}})}=&\frac{1}{\sqrt{2}}[(\text{cos}(\theta_{\text{HWP}})\ket{0_{t}}+\text{sin}(\theta_{\text{HWP}})\ket{1_{t}})(\text{sin}(\theta_{\text{HWP}})\ket{0_{0}}-\text{cos}(\theta_{\text{HWP}})\ket{1_{0}})\\ &+e^{i\phi}(\text{sin}(\theta_{\text{HWP}})\ket{0_{0}}-\text{cos}(\theta_{\text{HWP}})\ket{1_{0}})(\text{cos}(\theta_{\text{HWP}})\ket{0_{t}}+\text{sin}(\theta_{\text{HWP}})\ket{1_{t}})].\end{split} (31)

While the photon pairs go through the CBBO crystals which make horizontally polarized photons be delayed by a time tt respect vertically polarized photons, the final state after the CBBO crystals is

|ψout(θHWP)=12[(cos(θHWP)|02t+sin(θHWP)|1t)(sin(θHWP)|0tcos(θHWP)|10)+eiϕ(sin(θHWP)|0tcos(θHWP)|10)(cos(θHWP)|02t+sin(θHWP)|1t)].\begin{split}\ket{\psi_{\text{out}}(\theta_{\text{HWP}})}=&\frac{1}{\sqrt{2}}[(\text{cos}(\theta_{\text{HWP}})\ket{0_{2t}}+\text{sin}(\theta_{\text{HWP}})\ket{1_{t}})(\text{sin}(\theta_{\text{HWP}})\ket{0_{t}}-\text{cos}(\theta_{\text{HWP}})\ket{1_{0}})\\ &+e^{i\phi}(\text{sin}(\theta_{\text{HWP}})\ket{0_{t}}-\text{cos}(\theta_{\text{HWP}})\ket{1_{0}})(\text{cos}(\theta_{\text{HWP}})\ket{0_{2t}}+\text{sin}(\theta_{\text{HWP}})\ket{1_{t}})].\end{split} (32)

Then, we use two QWPs and one HWP in mode aa to control the phase to become eiϕ=1e^{i\phi}=-1. Thus, after the partial trace over the time system, the generated polarized state of the photon pairs is

ρent(θHWP)=12sin2(θHWP)(|10|01)(10|01|)+12cos2(θHWP)(|0101|+|1010|)=sin2(θHWP)|ΨΨ|+cos2(θHWP)ρwalk-off,\begin{split}\rho_{\text{ent}}(\theta_{\text{HWP}})&=\frac{1}{2}\text{sin}^{2}(\theta_{\text{HWP}})(\ket{10}-\ket{01})(\bra{10}-\bra{01})+\frac{1}{2}\text{cos}^{2}(\theta_{\text{HWP}})(\ket{01}\!\!\bra{01}+\ket{10}\!\!\bra{10})\\ &=\text{sin}^{2}(\theta_{\text{HWP}})\ket{\Psi^{-}}\!\!\bra{\Psi^{-}}+\text{cos}^{2}(\theta_{\text{HWP}})\rho_{\text{walk-off}},\end{split} (33)

where |Ψ=(|10|01)/2\ket{\Psi^{-}}=(\ket{10}-\ket{01})/\sqrt{2}. We utilize ρent(θHWP)\rho_{\text{ent}}(\theta_{\text{HWP}}) to calculate αideal\alpha_{\text{ideal}}, βideal\beta_{\text{ideal}}, and F¯s,ideal()\bar{F}_{s,\text{ideal}}(\mathcal{E}) in Fig. 4(a) in the main text.

If the setting angle is 9090^{\circ}, the unitary

UHWP(90)=[0110]U_{\text{HWP}}(90^{\circ})=\left[\begin{matrix}0&1\\ 1&0\end{matrix}\right]

makes the horizontal and vertical polarization be exchanged. Thus, Eqs. (31) and (32) for θHWP=90\theta_{\text{HWP}}=90^{\circ} are

|ψ′′(90)=12(|1t00+eiϕ|001t)\ket{\psi^{\prime\prime}(90^{\circ})}=\frac{1}{\sqrt{2}}(\ket{1_{t}0_{0}}+e^{i\phi}\ket{0_{0}1_{t}})

and

|ψout(90)=12(|1t0t+eiϕ|0t1t),\ket{\psi_{\text{out}}(90^{\circ})}=\frac{1}{\sqrt{2}}(\ket{1_{t}0_{t}}+e^{i\phi}\ket{0_{t}1_{t}}),

respectively. Since all photons of a pair have the same delay pairpair which means the photon walk-off effect is compensated after the CBBO crystals, after the partial trace over the time system, for eiϕ=1e^{i\phi}=-1, the generated state ρent(90)\rho_{\text{ent}}(90^{\circ}) in Eq. (33) is the ideal creation of the state ρent(90)=|ΨΨ|\rho_{\text{ent}}(90^{\circ})=\ket{\Psi^{-}}\!\!\bra{\Psi^{-}}. By contrast, if the setting angle θHWP=0\theta_{\text{HWP}}=0, 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 |0t10\ket{0_{t}1_{0}} and |100t\ket{1_{0}0_{t}} remain distinguishable and the generated state ρent(0)=ρwalk-off\rho_{\text{ent}}(0)=\rho_{\text{walk-off}} 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

ρent,noise(θHWP)=0.94ρent(θHWP)+0.06II4.\rho_{\text{ent,noise}}(\theta_{\text{HWP}})=0.94\rho_{\text{ent}}(\theta_{\text{HWP}})+0.06\frac{I\otimes I}{4}. (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 αtheory\alpha_{\text{theory}}, βtheory\beta_{\text{theory}}, F¯s,theory()\bar{F}_{s,\text{theory}}(\mathcal{E}), 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, RSP(0)\mathcal{E}^{(0)}_{\text{RSP}}, of U:R(0)=(|00|+|10|+|01||11|)/2U\!:~{}\!\!\!R(0)\!=\!(\ket{0}\!\!\bra{0}\!+\!\ket{1}\!\!\bra{0}\!+\!\ket{0}\!\!\bra{1}\!-\!\ket{1}\!\!\bra{1})\!/\!\sqrt{2}. To obtain the theoretical prediction of the process steerability of the RSPs varying with the walk-off compensation (θHWP\theta_{\text{HWP}}), we utilize ρent,noise(θHWP)\rho_{\text{ent,noise}}(\theta_{\text{HWP}}) 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 aa at 22.522.5^{\circ} to realize the operation UU^{{\dagger}} required on Alice’s side.

Second, we prepare the states |nmmn|\ket{n}_{mm}\!\!\bra{n} as input states by the polarization analyzer consisting of a QWP, a HWP, and a polarizer in mode aa. Note that |nm\ket{n}_{m} are eigenvectors of σm\sigma_{m} for n=0,1n=0,1 and m=1,2,3m=1,2,3 shown in the main text. The preparation basis σm\sigma_{m} is determined by adjusting the HWP and the QWP at different angles. For example, by adjusting the HWP and the QWP at angles of (22.5,0)\left(22.5^{\circ},0^{\circ}\right) with respect to the vertical axis, the preparation basis for the polarization states is set for σ1\sigma_{1}. Similarly, the preparation bases are set for σ2\sigma_{2} and σ3\sigma_{3} by adjusting the HWP and the QWP at angles of (0,45)\left(0^{\circ},45^{\circ}\right) and (0,0)\left(0^{\circ},0^{\circ}\right), respectively. Then, we set the polarizer at 00^{\circ} and 9090^{\circ} to prepare |0mm0|\ket{0}_{mm}\!\!\bra{0} and |1mm1|\ket{1}_{mm}\!\!\bra{1}, respectively.

Third, we perform state tomography to obtain the density matrix of output state (|nmmn|)\mathcal{E}(\ket{n}_{mm}\!\!\bra{n}) by the polarization analyzer in mode bb. As mentioned in the previous step, the measurement basis is also determined by the components of the polarization analyzer in mode bb. Taking |𝕤0=|0\ket{\mathbb{s}_{0}}=\ket{0}, |𝕤0=|1\ket{\mathbb{s}_{0}^{\bot}}=\ket{1} for example, our goal is to obtain the density matrix of (|0330|)\mathcal{E}(\ket{0}_{33}\!\bra{0}). We prepare the input state |0330|\ket{0}_{33}\!\bra{0} and |1331|\ket{1}_{33}\!\bra{1} by adjusting the HWP, the QWP, and the polarizer at angles of (0,0,0)\left(0^{\circ},0^{\circ},0^{\circ}\right) and (0,0,90)\left(0^{\circ},0^{\circ},90^{\circ}\right) in mode aa, respectively. Then, we perform the measurement of σ1\sigma_{1}, σ2\sigma_{2}, and σ3\sigma_{3} on (|0330|)\mathcal{E}(\ket{0}_{33}\!\bra{0}) in mode bb and record the coincident counts of them conditioned on |0330|\ket{0}_{33}\!\bra{0} and |1331|\ket{1}_{33}\!\bra{1} prepared in mode aa. The density matrix of (|𝕤0𝕤0|)\mathcal{E}\left(\ket{\mathbb{s}_{0}}\bra{\mathbb{s}_{0}}\right) and (|𝕤0𝕤0|)\mathcal{E}\left(\ket{\mathbb{s}_{0}^{\bot}}\bra{\mathbb{s}_{0}^{\bot}}\right) can then be constructed through coincidence measurement results via state tomography. According to the RSP protocol, the density matrix of (|0330|)\mathcal{E}(\ket{0}_{33}\!\bra{0}) can be obtained by the mean of (|𝕤0𝕤0|))\mathcal{E}\left(\ket{\mathbb{s}_{0}}\bra{\mathbb{s}_{0}}\right)) and (|𝕤0𝕤0|)\mathcal{E}\left(\ket{\mathbb{s}_{0}^{\bot}}\bra{\mathbb{s}_{0}^{\bot}}\right) after correction, i.e.,

(|0330|)=12(|𝕤0𝕤0|))+12Z(|𝕤0𝕤0|)Z.\mathcal{E}(\ket{0}_{33}\!\bra{0})=\frac{1}{2}\mathcal{E}\left(\ket{\mathbb{s}_{0}}\bra{\mathbb{s}_{0}}\right))+\frac{1}{2}Z\mathcal{E}\left(\ket{\mathbb{s}_{0}^{\bot}}\bra{\mathbb{s}_{0}^{\bot}}\right)Z^{\dagger}.

Fourth, considering the theoretical measurement results of ρent,noise(θHWP)\rho_{\text{ent,noise}}(\theta_{\text{HWP}}), we can get the theoretical output states (|nmmn|)\mathcal{E}(\ket{n}_{mm}\!\!\bra{n}) for n=0,1n=0,1 and m=1,2,3m=1,2,3. Therefore, the theoretical process matrix χtheory\chi_{\text{theory}} of the RSP can be obtained through process tomography Nielsen00 . According to χtheory\chi_{\text{theory}}, the process fidelity F,theoryF_{\mathcal{E},\text{theory}} and the process steerability αtheory\alpha_{\text{theory}}, βtheory\beta_{\text{theory}} can be calculated via SDP Lofberg ; sdpsolver . In addition, the theoretical average state fidelity F¯s,theory\bar{F}_{s,\text{theory}} can be obtained by F¯s,theory()=(2F,theory+1)/3\bar{F}_{s,\text{theory}}(\mathcal{E})=(2F_{\mathcal{E},\text{theory}}+1)/3.

As shown in Fig. 4(a) in the main text, the experimental results are consistent with the theoretical predictions of αtheory\alpha_{\text{theory}}, βtheory\beta_{\text{theory}}, and F¯s,theory()\bar{F}_{s,\text{theory}}(\mathcal{E}). Therefore, we can conclude that ρent,noise(θHWP)\rho_{\text{ent,noise}}(\theta_{\text{HWP}}) 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 𝒟\mathcal{D} Dakic10 and steerable weight SWSW 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 ρexpt\rho_{\text{expt}}, the geometric discord is calculated by the Bloch vector x(x1,x2,x3)\vec{x}\equiv\left(x_{1},x_{2},x_{3}\right) with xi=Tr[ρexpt(σiI)]x_{i}=\text{Tr}\left[\rho_{\text{expt}}\left(\sigma_{i}\otimes I\right)\right] and the correlation tensor AA with Aij=Tr[ρexpt(σiσj)]A_{ij}=\text{Tr}\left[\rho_{\text{expt}}\left(\sigma_{i}\otimes\sigma_{j}\right)\right] Dakic10 . Here, {σi(j)}i(j)=13\{\sigma_{i(j)}\}^{3}_{i(j)=1} is composed of σ1=X\sigma_{1}=X, σ2=Y\sigma_{2}=Y, and σ3=Z\sigma_{3}=Z Pauli matrices. The geometric discord of ρexpt\rho_{\text{expt}} takes the following form

𝒟(ρexpt)=14[Tr(xTx)+Tr(ATA)kmax],\mathcal{D}(\rho_{\text{expt}})=\frac{1}{4}\left[\text{Tr}\left(\vec{x}^{T}\vec{x}\right)+\text{Tr}\left(A^{T}A\right)-k_{\text{max}}\right], (35)

where the superscript TT indicates the transpose operation and kmaxk_{\text{max}} is the largest eigenvalue of matrix K=xxT+AATK=\vec{x}\vec{x}^{T}+AA^{T}. Thus, we can obtain the experimental geometric discord, denoted as 𝒟expt\mathcal{D}_{\text{expt}}, by ρexpt\rho_{\text{expt}}.

As the results shown in Fig. 4(a) in the main text, the density matrix of the created photon pairs at θHWP=40\theta_{\text{HWP}}=40^{\circ} in our experiment is of the following form:

ρexpt(40)=[0.180.130.04i0.04+0.02i0.01+0.09i0.13+0.04i0.400.28+0.04i0.010.01i0.040.02i0.280.04i0.300.08+0.05i0.010.09i0.01+0.01i0.080.05i0.12],\rho_{\text{expt}}(40^{\circ})=\left[\begin{matrix}0.18&0.13-0.04i&-0.04+0.02i&-0.01+0.09i\\ 0.13+0.04i&0.40&-0.28+0.04i&0.01-0.01i\\ -0.04-0.02i&-0.28-0.04i&0.30&-0.08+0.05i\\ -0.01-0.09i&0.01+0.01i&-0.08-0.05i&0.12\\ \end{matrix}\right], (36)

where the fidelity of ρexpt(40)\rho_{\text{expt}}(40^{\circ}) and the target state, |ΨΨ|\ket{\Psi^{-}}\!\bra{\Psi^{-}}, is F0.63F\sim 0.63 in our experiment. The resulting geometric discord of ρexpt(40)\rho_{\text{expt}}(40^{\circ}) is 𝒟expt0.08\mathcal{D}_{\text{expt}}\sim 0.08. Similarly, the theoretical geometric discord, 𝒟theory0.07\mathcal{D}_{\text{theory}}\sim 0.07, can be obtained by ρent,noise(40)\rho_{\text{ent,noise}}(40^{\circ}). With 𝒟theory0.07\mathcal{D}_{\text{theory}}\sim 0.07 and 𝒟expt0.08\mathcal{D}_{\text{expt}}\sim 0.08, 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 pnoise=0.5p_{\text{noise}}=0.5 in our experiment can be represented as

ρW,expt(0.5)=[0.130.01+0.01i0.01i0.01i0.010.01i0.380.23+0.01i0.01i0.01i0.230.01i0.360.010.01i0.01i0.01i0.01+0.01i0.13],\rho_{W,\text{expt}}(0.5)=\left[\begin{matrix}0.13&-0.01+0.01i&-0.01i&0.01i\\ -0.01-0.01i&0.38&-0.23+0.01i&0.01i\\ 0.01i&-0.23-0.01i&0.36&0.01-0.01i\\ -0.01i&-0.01i&0.01+0.01i&0.13\\ \end{matrix}\right], (37)

where the fidelity of ρW,expt(0.5)\rho_{W,\text{expt}}(0.5) and the target state, |ΨΨ|\ket{\Psi^{-}}\!\bra{\Psi^{-}}, is F0.59F\sim 0.59 in our experiment. The resulting geometric discord of ρW,expt(0.5)\rho_{W,\text{expt}}(0.5) is 𝒟expt0.10\mathcal{D}_{\text{expt}}\sim 0.10 and the theoretical geometric discord of ρW,noise(0.5)\rho_{W,\text{noise}}(0.5) is 𝒟theory0.11\mathcal{D}_{\text{theory}}\sim 0.11. 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, c\mathcal{E}_{c}.

VII.2 Steerable weight

Given an assemblage {ρnm}nm\{\rho^{\prime}_{nm}\}_{nm} of unnormalized states which Alice’s measurement affects Bob’s particle into, steerable weight is defined as the minimum amount of steerable resource in {ρnm}nm\{\rho^{\prime}_{nm}\}_{nm} that cannot be represented in the form Eq. (9). Any assemblage {ρnm}\{\rho^{\prime}_{nm}\} can be decomposed as

ρnm=μρnm,US+(1μ)ρnm,Sn,m,\rho^{\prime}_{nm}=\mu\rho^{\prime}_{nm,\text{US}}+(1-\mu)\rho^{\prime}_{nm,\text{S}}\ \ \ \ \forall n,m~{}, (38)

where 0μ10\leq\mu\leq 1 and {ρnm,US}nm\{\rho^{\prime}_{nm,\text{US}}\}_{nm} ({ρnm,S}nm\{\rho^{\prime}_{nm,\text{S}}\}_{nm}) is the unsteerable (steerable) assemblage in which the elements can (cannot) be represented in the form Eq. (9). The steerable weight SW(ρexpt)SW(\rho_{\text{expt}}) is then defined as SW(ρexpt)=1μSW(\rho_{\text{expt}})=1-\mu^{*}, where μ\mu^{*} denotes the maximum μ\mu 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, XX, YY, and ZZ for m=1,2,3m=1,2,3, respectively. Through the density matrices ρexpt\rho_{\text{expt}} and ρent,noise\rho_{\text{ent,noise}}, we can get the experimental and theoretical steerable weight, SWexptSW_{\text{expt}} and SWtheorySW_{\text{theory}}, respectively. As illustrated in Fig. 4(a) in the main text, the steerable weight of ρexpt(40)\rho_{\text{expt}}(40^{\circ}) shown in Eq. (37) is SWexpt0.02SW_{\text{expt}}\sim 0.02 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. Tu¨\ddot{\text{u}}tu¨\ddot{\text{u}}ncu¨\ddot{\text{u}}, 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).