Phase-Conjugate-State Pairs in Entangled States

We consider the probability that a bipartite quantum state contains phase-conjugate-state (PCS) pairs and/or identical-state pairs as signatures of quantum entanglement. While the fraction of the PCS pairs directly indicates the property of a maximally entangled state, the fraction of the identical-state pairs negatively determines antisymmetric entangled states such as singlet states. We also consider the physical limits of these probabilities. This imposes fundamental restrictions on the pair appearance of the states with respect to the local access of the physical system. For continuous-variable system, we investigate similar relations by employing the pairs of phase-conjugate coherent states. We also address the role of the PCS pairs for quantum teleportation in both discrete-variable and continuous-variable systems.


Introduction
The spooky action induced by the entangled particles at a distance is an interesting starting point in studies of quantum entanglement [1][2][3][4][5]. For the singlet state, the spins are always antiparallel to each other independently of the measured spin angle. Another familiar example is the Einstein-Podolski-Rosen (EPR) state. For the EPR states, not only the positions of the two particles but also the momentums of the two particles have perfect correlation. When we consider the standard form of the maximally entangled state in two d-level (qudit) systems |Φ 0,0 := d−1 n=0 |n |n / √ d, an interesting property is the appearance of a phase-conjugate state |ψ * := n ψ|n |n when a subsystem is projected onto a local state |ψ , namely, for any state |ψ , we have the relation An example with two spin systems is schematically shown in Figure 1. A natural question is whether this coherent appearance of the phase-conjugate-state (PCS) pairs {|ψ |ψ * } can be a signature of entanglement. It has been known that the use of PCS pairs {|ψ |ψ * } is more efficient than the use of identical-state pairs, {|ψ |ψ }, in the transmission of unknown quantum states. Gisin and Popescu showed that the set of antiparallel spins stores the quantum information better than the set of parallel spins [6]. Zhou et al. provided a generalization of this relation by exploiting unknown PCS and identical-state pairs for twoqudit systems [7]. Similar relations have been found in the continuous variable system in which a set of phase-conjugate coherent-state pairs, say {|α |α * }, is employed [8]. Other insightful properties of the pairs of coherent states have been investigated from the aspect of the optimal cloning [9]. A scheme of continuous variable cloning of phase-conjugate coherent states is proposed in [10] and an experimental demonstration has been reported in [11].
In this paper we define the probabilities that PCS and/or identical-state pairs appear in bipartite quantum states. We investigate physical limits of these probabilities and their relations to quantum entanglement. In Section 2, we define the fraction of the PCS and/or identical-state pairs on a two-qudit system and investigate their properties. For the case of the uniform distribution, the fractions determine the symmetry of bipartite quantum states in the sense of [12]. It is shown that the fraction of the PCS pairs corresponds to the average fidelity of quantum teleportation. We consider similar relations with respect to two mutually unbiased bases in Section 2.5. We also consider the phase-conjugate pairs of coherent states for a two-mode continuous-variable system in Section 3. The results are summarized in Section 4.  Figure 1: The pair of phase-conjugate states in the Bloch sphere when the fixed basis is the Z basis with the parameterization |ψ A = cos(θ/2)|0 + e iϕ cos(θ/2)|1 . In this representation the phase-conjugation changes the sign of the y-coordinate and the phase-conjugate state vector of the system B is given by |ψ * B = cos(θ/2)|0 + e −iϕ cos(θ/2)|1 . The directions of the Bloch vectors are in the reflection relation with respect to the z-x plane. The maximally entangled state (|0, 0 +|1, 1 )/ √ 2 is a coherent superposition of these pairs and produces the phase-conjugate pairs associated with the local projection in any direction as described in (1).

Two-Qudit System
This state has an interesting property associated with a local projection (see Figure 1): if a subsystem A is projected onto a pure state ψ then the state of the other subsystem B becomes the phase conjugate of the pure state |ψ * := d−1 n=0 |n ψ|n . In equation, we have the relation A ψ||Φ 0,0 AB = |ψ * B / √ d for any pure state ψ. Here, the factor 1/ √ d comes from the probability that the subsystem A is in ψ, that is, Tr ( A ψ||Φ 0,0 Φ 0,0 ||ψ A ) = 1/d. Note that the phase conjugation is defined with respect to a fixed basis. Any maximally entangled state has the same property up to local unitary operation.
In order to find the role of entanglement in such phenomena, we may define the fraction of the PCS pairs by where dψ is the uniform (Haar) measure which satisfies dψ|ψ ψ| = 1 d , and {dψ|ψ ψ|} corresponds to the positive operator valued measure (POVM) elements of the random projective measurement of a subsystem. The integrand Tr(|ψ ψ| ⊗ |ψ * ψ * | ρ) implies the probability that a bipartite state ρ contains the state pair |ψ |ψ * . Hence the fraction represents the probability that the local states of the system are in phase-conjugate relation. For the MES Φ 0,0 := |Φ 0,0 Φ 0,0 |, the PCS pair appears for any outcome of the random measurement as in (1), and we have the unit probability Now, we would like to ask two questions. (i) How high one can attain the probability of the pair appearance without entanglement? (ii) Does the unit probability f = 1 uniquely determine the MES? Is there another state that yields the pair appearance, unconditionally? In order to answer the question (i) we estimate the maximum value of the fraction under the constraint that the state is not entangled. We call a bipartite state classically correlated or separable if the state can be written in the form ρ AB = i p i ρ (i) A ⊗ ρ (i) B with the probability distribution p i ≥ 0 and p i = 1. If the state is not separable we call the state entangled. We define the classical limit fraction of the PCS pairs as the maximum value of the fraction f achieved by classically correlated states where the maximization is taken over separable states Since any state being invariant under the unitary operation U ⊗ U * (isotropic state) can be decomposed into the [13] we can write the isotropic operator in (2) as follows: where we use (3) in order to calculate the corresponding term The second term is the maximally entangled fraction achieved by separable states, which is known to be max ρ∈Sep Tr( Φ 0,0 ρ) = 1/d [13] (see Appendix A for a formal proof). This implies From the definition of the classical limit fraction, any separable state ρ must satisfy and cl is a signature of entanglement (see Figure 2). From (5), we can verify that where the minimum is achieved by the states in the orthogonal subspace of Φ 0,0 . A notable feature is that there is no physical state that achieves the probability f smaller than 1/(d + 1). In other words, any bipartite state must contain PCS pairs with probability no smaller than 1/(d + 1). This is a sort of inequality to generally limit the capability of the physical process similar to the case of the optimal cloning [9]. It gives a physical limit on the probability that concerns the pair appearance of local states with respect to the local measurement. The physically possible regime and classically possible regime for f are summarized in Figure 2.
From (5), it is also clear that the MES Φ 0,0 maximizes the probability of the pair appearance and f = 1 uniquely determines ρ = Φ 0,0 . We thus have positive answers to the question (ii) and obtain the separable condition from the classical limit on the fraction of the PCS pairs associated with the question (i).

Phase-Conjugate-State Pairs and Average
Fidelity of Quantum Teleportation. In this section we consider the process of quantum teleportation and show that the teleportation fidelity corresponds to the fraction of the PCS pairs introduced in the previous section (see (2)).
Let us define the generalized Pauli Z and X operators: e iωd j j j (10) with |d := |0 and ω d := 2π/d. They satisfy By using Z, X, and MES |Φ 0,0 = d−1 n=0 |n A |n B / √ d, we define the Bell states as follows: Note that Φ l,m |Φ l ,m AB = δ l,l δ m,m , and the set of states {|Φ l,m AB } with (l, m) ∈ {0, 1, 2, . . . , d − 1}, forms an orthonormal basis (Bell basis) on the two-qudit system. The process of the quantum teleportation [14] which transfers the state |ψ in in system A to system B by using a pre-shared entangled state ρ AB is described by where we defined the Kraus operator which satisfies If the shared entangled state ρ AB is Φ 0,0 , we have and the state on A is perfectly transferred into the state on B so that If the shared resource is a maximally entangled state but not Φ 0,0 , we can redefine the local basis so that it becomes the standard form. Then, the teleportation procedure can be proceeded in the same manner. For a given quantum channel Λ, we define the average fidelity as the probability that the state of the channel output corresponds to the state of the channel input averaged all over the input states [13], For the teleportation channel of (14) the average fidelity becomes as follows:  Figure 2: The probability of phase-conjugate-state-pair appearance f (ρ) defined in (2) must be no smaller than 1/(1 + d) where d is the dimension of each of the two subsystems. For separable states, holds, and thus the state which satisfies f > 2/(1 + d) must be entangled. The probability of identicalstate-pair appearance g(ρ) defined in (22) must be no larger than 2/(1 + d). For separable states, 1/(1 + d) ≤ g ≤ 2/(1 + d) holds, and thus the state which satisfies g < 1/(1 + d) must be entangled. This type of diagram was already presented in [12] where the maximally entangled fraction d Φ 0,0 and the fraction of the flip operator F are used to specify the conditions.
where we used the relations such as and changed the variable of the integration as | ψ = Z −m X −l |ψ in the last line. Then, we use l,m = d 2 and (2) to reach the final expression.
Equation (20) shows that the teleportation fidelity is determined by the fraction of the PCS pairs. Hence, we can understand that the role of the entangled state in quantum teleportation is to generate the phase-conjugate state in a (remote) local system when the state of the other local system is projected onto a pure state. Note that the classical boundary of the fidelity achieved by the entanglement breaking channel F cl = 2/(d + 1) [13] corresponds to the classical limit fraction of the PCS pairs f (d) cl defined in (4).

Identical-State Pairs in Two-Qudit States.
In this section we consider the probability that a two-qudit state contains identical-state pairs {|ψ |ψ }. We show its physical limits and relation to entanglement. The probability that a state contain an identical-state pair can be written by ψ| ψ| ρ|ψ |ψ . We define the fraction of the identical-state pairs by This is the probability that an identical-state pair appears subject to the random local measurement. It is defined parallel to the fraction of the PCS pairs of (2). In order to see the property of g we give a diagonal expression of the operator dψ|ψ ψ| ⊗ |ψ ψ|. Let us write the transposition of the second subsystem with respect to the fixed basis by Γ : (|l m|⊗| j k|) → (|l m|⊗|k j|). Since Γ(|φ φ| ⊗ |ψ ψ|) = |φ φ| ⊗ |ψ * ψ * | we can write from (5) where we use (5) and F d is the flip operator defined as The operator dψ|ψ ψ| ⊗ |ψ ψ| is invariant under the unitary operations U ⊗ U [15]. The quantum state being invariant under these operations is called the Werner state, and such a state is characterized by the coefficients of the identity 1 and flip operators F as the form of (23). We can see that the so-called singlet state |Ψ (±) is the eigenstate of the flip operator belong to the eigenvalue ±1. The superscript (±) corresponds to the sign of the eigenvalue ±1, respectively. It is also clear that the product states {| j | j } are the eigenstates of F belong to the eigenvalue 1. The set of the eigenstates {| j | j } and {|Ψ (±) i j } forms an orthonormal basis in the d 2 space. Therefore, we have decompositions of the identity operator and the flip operator: From these relations and (23) we obtain a diagonal form of the Werner-type operator: This expression implies 0 ≤ Γ( f ) ≤ 2/(1 + d), and we have This relation is parallel to the relation on f in (9) as in Figure 2. Similar to the case of f , it gives another physical limitation: the probability of the identical-pair appearance cannot be higher than 2/(1 + d) for any physical state. As one Journal of Atomic, Molecular, and Optical Physics 5 can see from (27) the separable states {| j | j } can achieve the maximum of g = 2/(1 + d). Thereby, any higher value of the identical-state-pair fraction cannot be a signature of entanglement different from the case of the PCS-pair fraction.
Next, we determine the minimum value of the fraction attained by separable states. Since the expectation value of a partial transposed operator for a product state corresponds to a expectation value of the original operator for another product state [16] Using this relation we obtain min ρ∈Sep. ρ∈Sep.
where we use (9). This implies that the state ρ is entangled if From (27) we can see that the left equality g = 0 is achieved by |Ψ (−) j,k with j / = k. For such states, there is no possibility of the identical-state-pair appearance.
To show a physical meaning of the identical-pair appearance, let us consider the following quantity: Since the operator 1 d − |ψ ψ| filters out the state ψ, the fraction of the orthonormal pairs |ψ |ψ ⊥ contributes to this quantity (and reduces the value of g). Hence, the identicalpair appearance g negatively quantifies the appearance of orthogonal-state pairs. The appearance of orthogonal-state pairs is characteristics of the singlet state in which the local states are antiparallel to each other as mentioned in the introduction. Our result implies that such a phenomena cannot occur for the case of separable states and is a signature of entanglement.

Summary of the Statements for the Uniform Average over the d-Level Local States.
Here, we summarize the main statements in Sections 2.1, 2.2, and 2.3. We have defined the fraction of the PCS pairs f ( ρ) (see (2)) and the fraction of the identical-state pairs g( ρ) (see (22)) in bipartite quantum states. We have shown that, for the probability f , the following statements hold (see also Figure 2).
We have shown that, for the probability g, the following statements hold (see also Figure 2).

Two Mutually Unbiased Bases.
In the previous sections the pair appearances are considered with respect to the uniform distribution that includes arbitrary PCS/identicalstate pairs. In this section we consider the appearance of the PCS pairs with respect to the elements of two mutually unbiased bases.
The two orthonormal bases of a d-level system, say {| j } and {| j }, are said to be mutually unbiased if they satisfy the relation | k| j | = 1/ √ d for any k and j [17]. Here we use a fixed Z basis {| j } and its Fourier basis {| j } defined by | j := Z j |0 with |0 := (1/ √ d) d−1 k=0 |k . We refer to {| j } as the X basis since they are eigenstates of X defined in (11). By definition, the elements of the Z basis are unchanged under the phase conjugation | j * = | j whereas the elements of the X basis changes under the phase conjugation as | j * = |− j .
Suppose that the subsystem A is randomly measured on either the X basis or the Z basis. We are interested in the probability that the state of the subsystem B is the phase-conjugate state. Let us define the positive operator as follows: The partial trace of this operator corresponds to the POVM of the random measurement of the two mutually unbiased bases (it satisfies the condition for the POVM on a local system Tr B σ = Tr A σ = 1 d ). We define the fraction of the PCS pair with respect to two mutually unbiased bases as follows. 6 Journal of Atomic, Molecular, and Optical Physics Using the completeness of the Bell basis of (13) we have From this expression we can verify 0 ≤ f m ≤ 1, where the maximum is achieved by Φ 0,0 and the minimum is achieved by any of Φ l,m in the last summation. Hence we could not find the physical limitation on the appearance of PCS pairs with respect to the mutually unbiased bases. Since the last term in (34) satisfies d−1 l,m=1 Φ l,m ≥ 0 we have σ ≤ (1/2)(1 + Φ 0,0 ). Using this relation and max ρ∈Sep Tr( Φ 0,0 ρ) = 1/d [13] we obtain The equality is achieved by the state |0, 0 . Equality in (35) gives the classical limit of f m . Consequently, we have the inseparable condition [18] associated with the PCS-pair appearance with respect to the two mutually unbiased bases: ρ is entangled if Next, we consider the lower bound of f m for classically correlated states. Interestingly, we can find that f m = 0 is achieved by separable states when the dimension d is not a prime number. For the separable state |φ |ϕ = d−1 j,k=0 φ j ϕ k | j |k we have Suppose that d = pq with integers p ≥ 2 and q ≥ 2, that is, d is not a prime number. Then we can verify σ = 0 is achieved by the separable state |φ |ϕ with |φ = (1/ The minimum eigenvalue of M corresponds to a minimum of σ achieved by the separable state of the following form: (α|1 +β|3 )⊗ (d−1)/2 k=0 u k |2k . Figure 3 shows the minimum eigenvalues for the first 20th prime numbers. For d ≥ 5, it is observed that min ρ∈Sep. f m ( ρ) is at least smaller than d −2 . Unfortunately, we could not have found nontrivial lower bound on this minimization problem for the prime dimensions d ≥ 5. However, it is likely that there is a such lower bound to give a separable condition for each prime number of d and that a smaller value of the identical-statepair fraction becomes a signature of entanglement when we employ two mutually unbiased bases for the prime dimensions.

Two-Mode Continuous-Variable System
In the following sections we will consider continuousvariable analog of the pair appearance of phase-conjugate states. The analysis in Sections 3.1 and 3.2 is continuousvariable counterpart of the analysis on the two-qudit system described in Sections 2.1 and 2.2.

Conjugate-State Pairs on Continuous-Variable Systems.
The EPR state is defined as a simultaneous eigenstate of the relative position and total momentum of a two-particle Journal of Atomic, Molecular, and Optical Physics Figure 4: If one of the EPR particles is projected onto a coherent state |α , the other particle turns out to be in the conjugate state |α * . This is a continuous-variable analog of the phase-conjugatepair appearance of Figure 1.
system. The EPR state thus has strong correlation on the positions and strong anticorrelation on the momentums. This suggests that the complex amplitudes of the two particles maintain the complex conjugate relation as in Figure 4. To be concrete, let us consider the form of the EPR state |EPR = |x |x dx = |p | − p dp [1]. Here, |x and |p := e ipx |x dx/ √ 2π are the eigenkets of the canonical operators x and p belong to the eigenvalues x and p, respectively. Here and in what follows we assume the commutation relations for canonical operators [ x, p] = i. We can verify ( x A − x B )|EPR = 0 and ( p A + p B )|EPR = 0. We can also rewrite this state in the coherent state basis by using the over completeness relation π −1 |α α|d 2 α = 1 as |EPR = π −1 |α |α * d 2 α. This expression demonstrates the pair appearance of the complex conjugate coherent states, {|α |α * }.
As an experimental implementation of the EPR state in the quantum optics, we usually work with the two-modesqueezed states (TMSS): where |ζ| 2 < 1 and we use the number states {|n } as the fixed basis. The TMSS is a simultaneous eigenstate of where ( a, a † ) and ( b, b † ) are the annihilation and creation operators of mode A and mode B, respectively [19]. We can verify the conditions for the simultaneous eigenstate as A ζ |ψ ζ AB = 0 and B ζ |ψ ζ AB = 0. The eigenvalue zeros for A ζ and B ζ imply a − ζ b † = 0 and b − ζ a † = 0. This suggests the existence of strong correlations on the complex amplitudes of the two modes. We can see that 1 − |ζ| 2 ) in the limit |ζ| → 1. Hence, in this limit, the TMSS approaches the EPR state.
In what follows we assume that ζ is real and positive. One may define the fraction of TMSS as an analog of the maximally entangled fraction in finite dimension systems. Since the TMSS |ψ ζ of (39) has 1 − ζ 2 as its largest Schmidt coefficient, using the result of Appendix A we have max ρ∈Sep. ψ ζ = 1 − ζ 2 and obtain the separable condition: Any separable state must satisfy For the TMSS of (39), we can see the following pair appearance phenomena. When one find that the local state of A is a coherent state |α := e −|α| 2 /2 ∞ n=0 α n |n / √ n!, the local state of B turns out to be another coherent state |ζα * . Similar to (1), the following relation holds for any complex amplitude α: In order to see the role of entanglement in such pair appearance we may define the fraction of the phase-conjugate coherent-state pairs as follows: where we introduced the real parameter λ ≥ 0 and the probability density In the limit λ → 0, we have λ/ p λ (α) → 1/π, and (p λ (α)/λ)|α α|d 2 α corresponds to the POVM elements of the double homodyne measurement on mode A (It satisfies |α α|d 2 α/π = 1). For notation convention we define the normalized state which satisfies Tr Ξ λ = 1. Then, we can write the fraction as For the TMSS of (39) we have This quantity takes the maximum value of when Note that, in the limit λ → 0, we have the unit probability f max = 1 and ζ → ζ. In what follows we show that (i) the maximum value of f in (43) for physical states is actually given by f max of (48) and determine (ii) the maximum value of f for classically correlated states.
Proof. Let us define the covariance matrix of ρ where R := ( x A , p A , x B , p B ) t is the set of the canonical operators of mode A and mode B. From (45) the covariance matrix of the operator Ξ λ is given by where 1 is the identity operator of 4-by-4 matrices.
In order to diagonalize this matrix we define another matrix U(r) corresponding to the two-mode squeezing operator U := e ( a † b † − a b)r that transforms the canonical operators as follows: The covariance matrix is diagonalized as follows: where and the squeezing parameter r is determined by which is equivalent to tanh r = ζ . Therefore, the operator Ξ is diagonalized by the two-mode squeezer as the product of thermal states: Thus, we have the maximum eigenvalue corresponding to the eigenstate U(r)|0, 0 = |ψ ζ . Equation (57) proves the first statement (i) f max (λ, ζ) = max ρ f R (λ, ζ, ρ).
Note that in the limit λ → 0, we have f max → 1 for any ζ > 0. This means that with nearly unit probability, the TMSS assigns the complex-conjugate pair of coherent states as the EPR state does (see Figure 4). Note also that the maximum eigenvalue can be associated with the physical limit of the amplification task in [20]. Now, our question is the limit of the pair appearance f in (43) for classical correlated states. We can prove the following statement: (ii) When the state has no entanglement the maximum value of the operator becomes Proof. In order to show this, let us write the partial transposition of the operator Ξ in the diagonal form: where V is the beam-splitter transformation act on the coherent states as V | 1 + ζ 2 α ⊗ |0 = |α ⊗ |ζα and we used the fact that the thermal state is diagonalized on the number-state basis in the final expression. From the diagonal form the maximum eigenvalue (operator norm) of Γ( Ξ λ ) is given by Γ( Ξ) ∞ = λ/(1 + λ + ζ 2 ) = λ f r (λ, ζ). Since an expectation value of a partial transposed operator for a product state corresponds to a expectation value of the original operator for another product state [16], we have Γ( Ξ) |a |b := Tr(Γ( Ξ)|a a| ⊗ |b b|) = Tr( ΞΓ(|a a| ⊗ |b b|)) = Tr( Ξ|a a| ⊗ |b * b * |) = Ξ |a |b * . Using this relation we have max ρ∈Sep.
Journal of Atomic, Molecular, and Optical Physics Entangled Separable 1 + ζ 2 Figure 5: The fraction of the state pairs f R (0, ζ, ρ) defined in (43) cannot exceed 1/(1 + ζ 2 ) for separable states, and the observation of the fraction beyond this limit is a signature of entanglement. This diagram is a continuous-variable analog of Figure 2. A counterpart of the identical-state-pair fraction g may be defined by lim λ → 0 [Γ(Ξ)/λ)]. Its classical limit is essentially zero.
Ancilla state for teleportation Figure 6: A continuous-variable quantum teleportation process E of an input coherent state |α from mode A to mode B by using a possibly entangled state ψ on modes A and B. The outcome z of the Bell measurement on the joint mode AA is forwarded to make a displacement D † (g z * ) on mode B, where g ≥ 0 is the teleportation gain.
Consequently, it has turned out that the pair appearance of the phase-conjugate coherent states can be a signature of entanglement, and we obtain the following statement: the state ρ is entangled if there exist λ > 0 and ζ > 0 such that f (λ, ζ, ρ) > f r (λ, ζ) [21] (see Figure 5).
In Section 2.3, it is shown that the fraction of the identical-state pairs vanishes for the antisymmetric states and a smaller value of g can be a signature of entanglement on the two-qudit system. For the fraction of the identical-coherentstate pairs {|α |ζα } with the prior distribution p λ (α), it is clear that one can lower the expectation value less than any given positive number by choosing a product of number states with a sufficiently large photon number. Hence, the fraction of the identical-coherent-state pairs is not useful as a signature of entanglement in this regard.

Teleportation Fidelity and the Pair Appearance of Phase-Conjugate Coherent States.
In this section we consider the process of continuous-variable quantum teleportation and investigate the equivalence between the teleportation fidelity and the fraction of the phase-conjugate coherent-state pairs.
Let us consider the process of continuous-variable quantum teleportation [22,23] for a coherent state |α from mode A to mode B by using a possibly entangled state ψ on the joint mode AB as in Figure 6. The teleportation process consists of the following two steps: (i) a Bell measurement of the joint mode AA and (ii) a feedforward displacement operation on mode B associated with the Bellmeasurement outcome (x, p). The first step is accomplished by a half-beamsplitter transformation R and two homodyne measurements. Suppose that the state of the joint mode AA is a product of coherent states |γ A |α A . Then, the joint probability density that the Bell-measurement outcome is (x, p) can be written as where the beam-splitter transformation is given by R AA , and the complex amplitude is defined by z := √ 2(x + ip). With this measurement outcome z and a gain parameter g ≥ 0, the second step can be described by the displacement operation D † (g z * ) on mode B. By using (61) and the P representation ψ BA = d 2 βd 2 γP(β, γ)|β β| B ⊗ |γ γ| A , we can express the state after the teleportation process as where we define Note that ψ BA is not a normalized density operator while α * |ψ BA |α * A is a normalized density operator on mode B. For the case of continuous-variable quantum channels [24] we may define the average gate fidelity as The factor ζ is essential to describe the effect of loss or amplification [20,21,24]. Such an effect is inherent concept in the class of Gaussian operations [25].
If we set the gain parameter g = ζ we can write with β = α + z * . From this expression we have where we performed the integration of z with a change of the integration variable d 2 β = d 2 α. Hence, we have shown that the teleportation fidelity (64) is equivalent to the fraction of the phase-conjugate coherent-state pairs with λ = 0. Note that the final expression of (66) no more includes the parameter λ. This is because the displacement all over the phase space eliminates the information about the phasespace position of ψ (see (63)).

Summary
We have considered the probability that a bipartite quantum state contains the PCS pairs and/or identical-state pairs. We determine the physical limits and classical limits of these probability for the case of uniform distribution on qudit states. The classical limits give the separable conditions. We have also shown the equivalence between the average fidelity of quantum teleportation process and the probability of the PCS pairs on the resource state of the teleportation. A summary of the obtained statements in Sections 2.1-2.3 is given in Section 2.4. For the case of uniform distribution on two mutually unbiased bases, a part of the problems becomes highly dependent on the dimension d. We have shown that the probability of the identical-state-pair appearance is useless for the entanglement verification when d is not a prime number. We have conjectured that the probability can be a signature of entanglement when d is a prime number. It is true when d = 2 and is numerically confirmed when d = 3.
We have also considered the probability that a two-mode continuous-variable state contains the phase-conjugate coherent-state pairs. We determine its physical limit and classical limit. We have also addressed its role in the process of continuous-variable quantum teleportation.

A. Maximum Separable Eigenvalue of Pure Entangled States
In the main text, we have considered the maximum expectation value of certain observables under the constraint that the state is separable. This optimization problem is called the separable eigenvalue problem [26]. If the observable is an entangled pure state, the maximum value is the square of the largest Schmidt coefficient. Here we provide a formal proof.
Proof. Any pure entangled state can be written in the Schmidt decomposed form as follows: λ k |k |k , (A.1) with 1 ≥ λ 0 ≥ λ 1 ≥ · · · ≥ λ km ≥ 0 and k λ 2 k = 1. For any pure separable state |s = i u i |i j v j | j with j |u j | 2 = 1 and j |v j | 2 = 1, we have j |u j v j | 2 ≤ i |u i | 2 j |v j | 2 = 1 from Schwarz inequality and obtain The equality is achieved when u k = v k = δ k,0 . Since any separable state can be written as a convex combination of pure separable states, we can verify The convexity also implies that the optimization over the pure separable states is always sufficient for the separable eigenvalue problem.

B. Another Expression of the Flip Operator
Here, we give another diagonal expression of the flip operator defined in (24). In the main text we use the diagonal expression of (26). The matrix elements of F d with respect to the Bell basis of (13) is given by Φ l,m | F d |Φ l ,m = e iωdlm δ m,m δ l,d−l . From this relation we can diagonalize the flip operator:

C. Generation of the Fraction in a Symmetric Form
In the expression of f R of (43), the systems A and B are not equivalently treated as the parameter ζ is only put on the states of B. Here, we will consider a symmetric formula and present a related separable condition. Let us define the operator, in which the position of ζ is moved from B to A on (45), Then, from the same procedure done on Ξ, we obtain the separable condition (see (60)) We can verify the following relation for the convex combination of operators, with χ > 0: The bound can be achieved by ρ 0 = |0, 0 0, 0| again. Using this inequality we can obtain the following separable criterion. For any ξ > 0, and λ 1 , λ 2 , ζ 1 , ζ 2 > 0, separable states must satisfy the following: ≤ χ −1 f max (λ 1 , ζ 1 ) + χ f max (λ 2 , ζ 2 ) . (C.4)