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

Enhancing teleportation via noisy channels: effects of the induced multipartite entanglement

Victor H. T. Brauer [email protected]    Andrea Valdés-Hernández [email protected] Instituto de Física, Universidad Nacional Autónoma de México
Apartado Postal 20-364, Ciudad de México, Mexico
Abstract

Quantum teleportation in the presence of noisy channels acting on a bipartite resource state is considered. We consider a family of generalized noisy channels that continuously connect the amplitude damping and the dephasing channels, encompassing a wide family of in-between scenarios, to delve into the relation between the teleportation success and the amount of 3- and 4-partite entanglement (distributed among the qubits of the resource state and those representing local environments) generated during the evolution. Our analysis reveals that for a fixed entanglement of the resource state, the channels that better protect the teleportation fidelity against the detrimental effects of noise are those that generate higher amounts of (ghz-type) multipartite entanglement This suggests that the dynamically induced multipartite correlations may serve as an additional resource for teleportation, and throws light into the characterization of processes, and of the type of induced entanglement, according to their capability of assisting the protocol.

I Introduction

Quantum teleportation stands out as one of the most fascinating applications of quantum entanglement. It allows the transmission of quantum information between two spatially separated agents, Alice and Bob, by means of local operations, classical communication, and a key element: a shared correlated state known as resource state 111Some authors refer also to the resource state as quantum channel. Here we will keep the term ‘quantum channel’ to denote a completely positive trace-preserving map, that in the present context gives the reduced dynamics of an open system.. The original standard teleportation protocol [2] uses a (pure, maximally entangled, two-qubit) Bell state as the resource state. Subsequent investigations extended the scheme noticing that there exist resource states that are useful for teleportation yet are neither maximally entangled [3] nor pure, and their relation with violations of Bell inequalities [4, 5] and discord-like correlations [6] has been discussed.

A mixed rather than a pure resource state is more realistic, particularly when taking into consideration the interaction of the (Alice and Bob’s) entangled pair with its surroundings, resulting in mixing of the resource state. There has been extensive research on such noisy quantum teleportation schemes [7, 8, 9, 10, 11, 12, 13, 14, 15, 16], e.g., resorting to the Lindblad formalism to explore the fidelity of teleportation in terms of decoherence rates [7], or to determine the optimal Bell resource state under different local Pauli noises [8]. Recent methods for protecting teleportation against some decoherence channels have also been advanced [9, 10]. The effect on the teleportation fidelity of different noisy channels acting on the resource state —typically representing the interaction of Alices’s and Bob’s particles with additional subsystems (which in the present case are regarded as local environments)—, has been studied considering the Kraus operators corresponding to dissipative interactions via an amplitude damping channel [11, 12], together with other paradigmatic noise or decoherence channels on qubits such as bit flip, phase flip, depolarizing [13], and phase damping [14]. Also, the teleportation protocol under noisy channels in higher dimensional systems has been explored [15]. A general theoretical and experimental [16] conclusion that ensues from these investigations is that there exist appropriate channels (acting on suitable initially pure resource states) for which the detrimental effects of noise on the teleportation fidelity are minimal, compared to other noisy channels.

Further extensions of the original teleportation protocol have also been advanced that consider multi-party resource states exhibiting some type of multipartite entanglement. This has led to the development of strategies that exploit multipartite entanglement to teleport multiple qubit states, as e.g. in [17, 18, 19, 20, 21]. Multi-directional teleportation, allowing quantum information transmission between several agents, has also been explored [22, 23, 24, 25], including the effect of noisy channels [26, 27, 28, 29].

Despite the advances achieved regarding the teleportation success under noisy channels, whether acting on bipartite or multipartite resource states, the relation between the teleportation fidelity and the multipartite entanglement generated among the resource qubits and the environment, has been much less explored. Such an analysis would allow us to identify the type of processes —characterized by the type of entanglement they induce—, that favor a more successful teleportation, and to possibly explain the fidelity improvement as an effect assisted by the created multipartite entanglement. In [30, 31] some progress has been made, by relating the teleportation fidelity with the 3-partite entanglement resulting from the local interaction of one qubit of the resource state with a two-level environment. Here we contribute along these lines, by focusing on the standard teleportation protocol in the presence of noisy channels acting on the bipartite resource state. We consider a generalized noisy quantum channel that can be continuously transformed from the amplitude damping to the dephasing channel [32], and explore the correlation between the maximal average fidelity and the amount of 3- and 4-partite entanglement, distributed among the pair of qubits that conform the resource state and the qubits that represent the corresponding local environments.

We first revisit the standard teleportation protocol and the notion of maximal average fidelity FmaxF_{\max} as a quantifier of the teleportation success (Sec. II), and express the latter in terms of the Kraus operators of an arbitrary quantum channel acting on the (arbitrary yet initially pure) resource state, conformed by two qubits AA and BB (Sec. III). After these preliminary sections, we relate the maximal average fidelity above the classical threshold value, 𝔽max\mathbb{F}_{\max}, with the bipartite entanglement between AA and BB, considering A+BA+B as an ideal closed system (Sec. IV). Assuming then that BB interacts with a local, two-level environment EBE_{B} via the generalized channel, we investigate the relation between 𝔽max\mathbb{F}_{\max} and the 3-partite entanglement distributed among A,BA,B and EBE_{B} (Sec. V). The analysis is extended to the 4-partite case by considering that both AA and BB locally interact with their respective environments EAE_{A} and EBE_{B} under independent generalized channels, and a non-trivial correlation between 𝔽max\mathbb{F}_{\max} and the 4-partite entanglement is disclosed (Sec. VI). Finally, some concluding remarks are presented (Sec. VII).

II The Standard Teleportation Protocol and maximal average fidelity

The main idea behind the standard teleportation protocol is that Alice wants to send to Bob an arbitrary input state, ρin\rho_{\textrm{in}}, encoded in a qubit aa in her possession. For this task, they share a pair of qubits AA and BB (in Alice and Bob’s possession, respectively) in an entangled state ρAB\rho_{AB}, called resource state. Alice then performs a Bell measurement [33, 34, 35] on her pair of qubits aa and AA, and communicates the outcome to Bob via a classical channel. Upon receiving this information —and knowing ρAB\rho_{AB} [36]—, Bob performs a unitary operation σ(i){𝖨2,σx,σy,σz}\sigma^{(i)}\in\{\mathsf{I}_{2},\sigma^{x},\sigma^{y},\sigma^{z}\} on his qubit BB, thus putting it into the output state ρout\rho_{\textrm{out}} (σx,y,z\sigma^{x,y,z} stand for the Pauli matrices, and 𝖨n\mathsf{I}_{n} denotes the n×nn\times n identity operator).

Refer to caption
Figure 1: Scheme of the standard teleportation protocol. Vertical dashed lines divide the different stages that correspond, from left to right, to: i) Preparation of the initial resource state of the system A+BA+B; ii) Application of a Bell measurement in Alice’s qubits (aa and AA); iii) Post-processing according to Alice’s results.

Figure 1 shows a schematic generalization of the standard teleportation protocol. The qubits AA and BB are initially in the state |0\left|{0}\right\rangle, whereas aa is already in ρin\rho_{\textrm{in}}, the state to be teleported. Vertical dashed lines separate the different stages of the protocol corresponding, from left to right, to:

  • A unitary transformation U(ϕ,φ)U(\phi,\varphi) rotates the qubit AA, and a cnot gate is employed to prepare the initial resource state

    |ϕ0AB=cosϕ|00+eiφsinϕ|11,\left|{\phi_{0}}\right\rangle_{AB}=\cos\phi\left|{00}\right\rangle+e^{i\varphi}\sin\phi\left|{11}\right\rangle, (1)

    with ϕ[0,π/2]\phi\in[0,\pi/2], and φ[π/2,3π/2]\varphi\in[-\pi/2,3\pi/2]. (This generalizes the application of the Hadamard gate, resulting in the Bell state 12(|00+|11)\frac{1}{\sqrt{2}}(\left|{00}\right\rangle+\left|{11}\right\rangle)).

  • Alice applies a cnot and a Hadamard gate to the qubits in her possession, then performs a measurement in the Bell basis and communicates the outcome to Bob via classical channels.

  • Depending on Alice’s measurement outcome, Bob applies appropriate unitary operations on his qubit BB, so the output state of BB coincides with the input state of aa. The whole strategy leads to perfect teleportation if the resource state (1) is the Bell state 12(|00+|11)\frac{1}{\sqrt{2}}(\left|{00}\right\rangle+\left|{11}\right\rangle) (if |ϕ0\left|{\phi_{0}}\right\rangle was another Bell state, then the operations in Bob’s strategy must be changed).

For pure input states, ρin=|χχ|\rho_{\textrm{in}}=\left|{\chi}\right\rangle\!\left\langle{\chi}\right|, the success of the teleportation can be quantified by means of the maximal average fidelity FmaxF_{\max}, which measures the probability that ρout\rho_{\textrm{out}} coincides with the (unknown) ρin\rho_{\textrm{in}}, averaged over all input states, provided the appropriate σ(i)\sigma^{(i)} is chosen. The maximal average fidelity for an arbitrary 2-qubit resource state ρAB\rho_{AB} can be written as [37]

Fmax=13[2max(ρAB)+1],F_{\max}=\frac{1}{3}\Big{[}2\mathcal{F}_{\max}(\rho_{AB})+1\Big{]}, (2)

where max\mathcal{F}_{\max} is the maximal singlet fraction [38]

max=maxi{Φi|ρAB|Φi}\mathcal{F}_{\max}=\max_{i}\{\langle\Phi_{i}|\rho_{AB}|\Phi_{i}\rangle\} (3)

corresponding to the maximum fidelity between the resource state and any of the Bell states 222In Ref. [38] the maximal singlet fraction is defined considering the maximum over all the maximally entangled states. Here we maximize only over those states that can be obtained from Bell states by means of unitary transformations of the form 𝖨2σ(i)\mathsf{I}_{2}\otimes\sigma^{(i)}. With this restriction we adhere to the standard teleportation protocol (in which Bob’s operations are implemented via Pauli operators).

|Φi=(𝖨2σ(i))|Φ+,|Φ+=12(|00+|11).\left|{\Phi_{i}}\right\rangle=(\mathsf{I}_{2}\otimes\sigma^{(i)})\left|{\Phi^{+}}\right\rangle,\quad\left|{\Phi^{+}}\right\rangle=\frac{1}{\sqrt{2}}(\left|{00}\right\rangle+\left|{11}\right\rangle). (4)

Putting σ(0)=𝖨2\sigma^{(0)}=\mathsf{I}_{2}, and σ(1),(2),(3)=σx,y,z\sigma^{(1),(2),(3)}=\sigma^{x,y,z}, we get

|Φ0/3=|Φ±=12(|00±|11),|Φ1/2=|Ψ±=12(|01±|10).\begin{split}\left|{\Phi_{0/3}}\right\rangle=\left|{\Phi^{\pm}}\right\rangle=\frac{1}{\sqrt{2}}(\left|{00}\right\rangle\pm\left|{11}\right\rangle),\\ \left|{\Phi_{1/2}}\right\rangle=\left|{\Psi^{\pm}}\right\rangle=\frac{1}{\sqrt{2}}(\left|{01}\right\rangle\pm\left|{10}\right\rangle).\end{split} (5)

Therefore, max\mathcal{F}_{\max} picks out the optimal Bell state, i.e., the Bell state that is closest to ρAB\rho_{AB}. The above expressions show that if ρAB\rho_{AB} is a Bell state, the protocol guarantees the complete reconstruction of the input state.

For the quantum teleportation to be considered successful FmaxF_{\max} must be greater than 2/32/3, which is the value of the maximal average fidelity corresponding to the best possible reconstruction through a purely classical channel [40]. Throughout this paper we will be interested in the success of non-classical teleportation, hence we will focus on those regimes in which Fmax2/3F_{\max}\geq 2/3, and accordingly pay attention to the quantity

𝔽max=max{23,Fmax}.\mathbb{F}_{\max}=\max\Big{\{}\frac{2}{3},F_{\max}\Big{\}}. (6)

III Maximal average fidelity under quantum channels

Let us assume that AA and BB are initially prepared in the state |ϕ0AB\left|{\phi_{0}}\right\rangle_{AB}, given by Eq. (1), which then passes through a quantum channel ΛAB\Lambda_{AB} represented by a set of Kraus operators {Πμ}\{\Pi_{\mu}\}, acting on AB\mathcal{H}_{A}\otimes\mathcal{H}_{B} [41]. A scheme of this process is shown in Fig. 2.

Refer to caption
Figure 2: The pair of qubits AA and BB, once in the entangled state |ϕ0\left|{\phi_{0}}\right\rangle, passes through an arbitrary quantum channel represented by the set of Kraus operators Πμ\Pi_{\mu}, resulting in a mixed resource state ρAB\rho_{AB}.

The channel thus transforms |ϕ0ϕ0|\left|{\phi_{0}}\right\rangle\!\left\langle{\phi_{0}}\right| into an effective (typically mixed) resource state ρAB\rho_{AB} given by

ρAB=ΛAB(|ϕ0ϕ0|)=μΠμ|ϕ0ϕ0|Πμ.\rho_{AB}=\Lambda_{AB}(\left|{\phi_{0}}\right\rangle\!\left\langle{\phi_{0}}\right|)=\sum_{\mu}\Pi_{\mu}\left|{\phi_{0}}\right\rangle\!\left\langle{\phi_{0}}\right|\Pi^{\dagger}_{\mu}. (7)

Direct substitution into Eq. (3) gives, with the aid of (2),

Fmax=13+23maxi{μ|Φi|Πμ|ϕ0|2}.\begin{split}F_{\max}&=\frac{1}{3}+\frac{2}{3}\max_{i}\Big{\{}\sum_{\mu}\big{|}\langle\Phi_{i}\big{|}\Pi_{\mu}\big{|}\phi_{0}\rangle\big{|}^{2}\Big{\}}.\end{split} (8)

This simple but central expression for FmaxF_{\max} in terms of the Kraus operators allows for studying the teleportation success when the initially pure resource state is subject to an arbitrary quantum channel.

Since Alice and Bob are typically spatially separated, we will focus on cases in which AA and BB undergo independent local channels, so ΛAB=ΛAΛB\Lambda_{AB}=\Lambda_{A}\otimes\Lambda_{B} and

ΠμΠαβ=QαKβ,\Pi_{\mu}\rightarrow\Pi_{\alpha\beta}=Q_{\alpha}\otimes K_{\beta}, (9)

where {Qα}\{Q_{\alpha}\} and {Kβ}\{K_{\beta}\} stand for the sets of Kraus operators associated to ΛA\Lambda_{A} and ΛB\Lambda_{B}, respectively, each set having at most (dimA(B))2=4(\dim\mathcal{H}_{A(B)})^{2}=4 elements. Equation (8) becomes then

Fmax=13+23maxi{αβ|Φi|QαKβ|ϕ0|2}.\begin{split}F_{\max}&=\frac{1}{3}+\frac{2}{3}\max_{i}\Big{\{}\sum_{\alpha\beta}\big{|}\langle\Phi_{i}\big{|}Q_{\alpha}\otimes K_{\beta}\big{|}\phi_{0}\rangle\big{|}^{2}\Big{\}}.\end{split} (10)

The channel ΛAΛB\Lambda_{A}\otimes\Lambda_{B} acts as if AA and BB interact locally with a corresponding party EAE_{A} and EBE_{B}. Assuming that the systems (A+B)(A+B), EAE_{A} and EBE_{B} are initially uncorrelated, and that EAE_{A} and EBE_{B} are qubits in the initial state |0\left|{0}\right\rangle, we may write

|ψ0ABEAEB=|ϕ0AB|0EA|0EB\left|{\psi_{0}}\right\rangle_{ABE_{A}E_{B}}=\left|{\phi_{0}}\right\rangle_{AB}\otimes\left|{0}\right\rangle_{E_{A}}\otimes\left|{0}\right\rangle_{E_{B}} (11)

for the initial 4-partite state. Further, if the interaction between AA and EAE_{A} is represented by the unitary operator UAEAU_{AE_{A}} (and similarly for BB and EBE_{B}), then

Qα=α|UAEA|0EAEA,Kβ=β|UBEB|0EBEB,\begin{split}Q_{\alpha}={}_{E_{A}}\!\left\langle{\alpha}\right|U_{AE_{A}}\left|{0}\right\rangle_{E_{A}},\\ K_{\beta}={}_{E_{B}}\!\left\langle{\beta}\right|U_{BE_{B}}\left|{0}\right\rangle_{E_{B}},\end{split} (12)

with {|α}\{\left|{\alpha}\right\rangle\} and {|β}\{\left|{\beta}\right\rangle\} basis of EA\mathcal{H}_{E_{A}} and EB\mathcal{H}_{E_{B}}, respectively. In addition, the (unitary) evolution of the complete system can be obtained from these Kraus operators as

|ψABEAEB\displaystyle\left|{\psi}\right\rangle_{ABE_{A}E_{B}} =\displaystyle= UAEAUBEB|ψ0ABEAEB\displaystyle U_{AE_{A}}U_{BE_{B}}\left|{\psi_{0}}\right\rangle_{ABE_{A}E_{B}} (13)
=\displaystyle= αβQαKβ|ϕ0AB|αEA|βEB.\displaystyle\sum_{\alpha\beta}Q_{\alpha}K_{\beta}\left|{\phi_{0}}\right\rangle_{AB}\left|{\alpha}\right\rangle_{E_{A}}\left|{\beta}\right\rangle_{E_{B}}.

A particular family of local channels, that involves a 3-party system instead of a 4-partite one, is that in which one of the resource qubits, say AA, remains unaffected so ΛA=𝖨2\Lambda_{A}=\mathsf{I}_{2}, while BB goes through an arbitrary channel ΛB\Lambda_{B}. In this case Qα=𝖨2δα0Q_{\alpha}=\mathsf{I}_{2}\,\delta_{\alpha 0}, and from Eq. (10) we are led to

Fmax=13+23maxi{β|Φi|(𝖨2Kβ)|ϕ0|2}.F_{\max}=\frac{1}{3}+\frac{2}{3}\max_{i}\Big{\{}\sum_{\beta}\big{|}\langle\Phi_{i}\big{|}(\mathsf{I}_{2}\otimes K_{\beta})|\phi_{0}\rangle\big{|}^{2}\Big{\}}. (14)

Clearly for ΛAB=𝖨2ΛB\Lambda_{AB}=\mathsf{I}_{2}\otimes\Lambda_{B} the subsystem EAE_{A} is superfluous, and Eq. (11) reduces to the 3-qubit initial state

|ψ0ABEB=|ϕ0AB|0EB.\left|{\psi_{0}}\right\rangle_{ABE_{B}}=\left|{\phi_{0}}\right\rangle_{AB}\otimes\left|{0}\right\rangle_{E_{B}}. (15)

IV Fidelity and bipartite entanglement in the standard teleportation protocol

The standard protocol, depicted in the circuit of Fig. 1, corresponds to ΛAB=𝖨4\Lambda_{AB}=\mathsf{I}_{4}. Equation (8) thus becomes

Fmax=13+23maxi{|Φi|ϕ0|2}.F_{\max}=\frac{1}{3}+\frac{2}{3}\max_{i}\Big{\{}\big{|}\langle\Phi_{i}\big{|}\phi_{0}\rangle\big{|}^{2}\Big{\}}. (16)

With |ϕ0\left|{\phi_{0}}\right\rangle given by (1), the optimal strategy corresponds to the state |Φ+\left|{\Phi^{+}}\right\rangle for φ[π2,π2]\varphi\in[-\frac{\pi}{2},\frac{\pi}{2}], and |Φ\left|{\Phi^{-}}\right\rangle for φ[π2,3π2]\varphi\in[\frac{\pi}{2},\frac{3\pi}{2}], so (16) reduces to Fmax=23+130|cosφ|F_{\max}=\frac{2}{3}+\frac{1}{3}\,\mathcal{E}_{0}\,|\!\cos\varphi|, and consequently, in the non-interacting case, Eq. (6) gives

𝔽maxnon-int=23+130|cosφ|,\mathbb{F}^{\textrm{non-int}}_{\max}=\frac{2}{3}+\frac{1}{3}\,\mathcal{E}_{0}\,|\!\cos\varphi|, (17)

where 0\mathcal{E}_{0} stands for the entanglement of the initial resource state

0C(|ϕ0)=2cosϕsinϕ=sin2ϕ.\mathcal{E}_{0}\equiv C(|\phi_{0}\rangle)=2\cos\phi\sin\phi=\sin 2\phi. (18)

Here CC is the concurrence, quantifying the amount of qubit-qubit entanglement [42]. For an arbitrary (in general mixed) 2-qubit state ρAB\rho_{AB} it is defined as

CABC(ρAB)=max{0,λ0λ1λ2λ3},C_{AB}\equiv C(\rho_{AB})=\max\{0,\sqrt{\lambda_{0}}-\sqrt{\lambda_{1}}-\sqrt{\lambda_{2}}-\sqrt{\lambda_{3}}\}, (19)

where {λn}\{\lambda_{n}\} are the eigenvalues of the matrix ρAB(σyσy)ρAB(σyσy)\rho_{AB}(\sigma^{y}\otimes\sigma^{y})\rho^{*}_{AB}(\sigma^{y}\otimes\sigma^{y}) ordered in decreasing order, and ρAB\rho^{*}_{AB} stands for the complex conjugate of ρAB\rho_{AB} expressed in the computational basis. When the state is pure (ρAB=|χχ|\rho_{AB}=\left|{\chi}\right\rangle\!\left\langle{\chi}\right|), the expression for the concurrence simplifies and reads

CAB=C(|χ)=2(1TrρA2)=2(1TrρB2),C_{AB}=C(\left|{\chi}\right\rangle)=\sqrt{2(1-\textrm{Tr}\rho^{2}_{A})}=\sqrt{2(1-\textrm{Tr}\rho^{2}_{B})}, (20)

with ρA(B)=TrB(A)ρAB\rho_{A(B)}=\textrm{Tr}_{B(A)}\rho_{AB} the reduced density matrix of either one of the qubits.

Equation (17) makes explicit that the teleportation success enhances as the resource state’s entanglement 0\mathcal{E}_{0} increases. It also shows that as φ\varphi tends to ±π/2\pm\pi/2, the fidelity decreases up to its minimal —classically attainable— value 2/32/3, irrespective of the initial entanglement, as can be seen in Fig. 3.

Refer to caption
Figure 3: Maximal average fidelity in absence of interaction (Eq. (17)) as a function of φ\varphi and ϕ\phi, in color scale ranging from 2/32/3 to 1. At φ=±π/2\varphi=\pm\pi/2, 𝔽maxnon-int\mathbb{F}^{\textrm{non-int}}_{\max} shows no improvement with respect to the classically attainable value (2/32/3), irrespective of ϕ\phi, which determines the initial entanglement 0=sin2ϕ\mathcal{E}_{0}=\sin 2\phi.

That is, there are maximally entangled resource states for which the maximal average fidelity does not exceed its classical limit, and the relative phase φ\varphi determines the fraction of 0\mathcal{E}_{0} that ultimately improves the fidelity.

V Fidelity and three-partite entanglement generated via a noisy channel

We now focus on the case in which only BB undergoes through a quantum channel, effectively representing an interaction with an additional qubit EBE_{B}. Figure 4 illustrates this situation, leading to a mixed resource state ρAB\rho_{AB} before the measurement stage of the protocol.

Refer to caption
Figure 4: Qubits AA, BB and EBE_{B} are prepared in the initial state (15). Then BB and EBE_{B} interact via a unitary operation UBEBU_{BE_{B}}, and the effect on BB is that of a quantum channel described by the Kraus operators {Kβ}\{K_{\beta}\}.

The scenario under consideration corresponds to that in which |ϕ0ϕ0|\left|{\phi_{0}}\right\rangle\!\left\langle{\phi_{0}}\right| is subject to a channel ΛAB=𝖨2ΛB\Lambda_{AB}=\mathsf{I}_{2}\otimes\Lambda_{B}, where ΛB\Lambda_{B} encodes the interaction between BB and EBE_{B}, giving rise to the possible creation of tripartite entanglement in the system A+B+EBA+B+E_{B}. The channel thus transforms the initial state (15)

|ψ0ABEB=cosϕ|000+eiφsinϕ|110\left|{\psi_{0}}\right\rangle_{ABE_{B}}=\cos\phi\left|{000}\right\rangle+e^{i\varphi}\sin\phi\left|{110}\right\rangle (21)

into the 3-qubit state (see Eq. (13) with Qα=𝖨2δα0Q_{\alpha}=\mathsf{I}_{2}\,\delta_{\alpha 0})

|ψABEB=βKβ|ϕ0AB|βEB.\left|{\psi}\right\rangle_{ABE_{B}}=\sum_{\beta}K_{\beta}\left|{\phi_{0}}\right\rangle_{AB}\left|{\beta}\right\rangle_{E_{B}}. (22)

In [43] necessary and sufficient conditions on the Kraus operators Kβ=β|UBEB|0K_{\beta}=\left\langle{\beta}\right|U_{BE_{B}}\left|{0}\right\rangle were established, that ensure the emergence of bipartite and tripartite entanglement among the parties A,B,EBA,B,E_{B}. For the present analysis we concentrate on the dynamics of the resource state’s (bipartite) entanglement CAB=C(ρAB)C_{AB}=C(\rho_{AB}), and the tripartite entanglement, as measured by the so-called 3-tangle. The latter stands as a legitimate measure of residual entanglement in a 3-qubit pure state |ψijk\left|{\psi}\right\rangle_{ijk}, and quantifies the amount of three-way ghz-type entanglement in the state [44]. The 3-tangle is defined as [45]

τijk=τ(|ψijk)=Ci|jk2Cij2Cik2,\tau_{ijk}=\tau(\left|{\psi}\right\rangle_{ijk})=C^{2}_{i|jk}-C^{2}_{ij}-C^{2}_{ik}, (23)

where CijC_{ij} is given by Eq. (19), and Ci|jkC_{i|jk} by

Ci|jk=2(1Trρi2).C_{i|jk}=\sqrt{2(1-\textrm{Tr}\rho^{2}_{i})}. (24)

The last expression generalizes (20) for a 3-party pure state and quantifies the entanglement across the bipartition i|(j+k)i|(j+k) [46].

In [43] it is found that for ΛAB=𝖨2ΛB\Lambda_{AB}=\mathsf{I}_{2}\otimes\Lambda_{B}, CABC_{AB} and τABEB\tau_{ABE_{B}} evolve according to

CAB2\displaystyle C^{2}_{AB} =\displaystyle= 02(|detK0|+|detK1|)2\displaystyle\mathcal{E}^{2}_{0}(|\det K_{0}|+|\det K_{1}|)^{2} (25)
1202(|u||v|+|vu|),\displaystyle-\frac{1}{2}\mathcal{E}^{2}_{0}\big{(}|u|-|v|+|v-u|\big{)},

and

τABEB=02|uv|,\tau_{ABE_{B}}=\mathcal{E}^{2}_{0}|u-v|, (26)

where u=4detK0K1u=4\det K_{0}K_{1} and v=g2(K0,K1)v=g^{2}(K_{0},K_{1}), with g(M,N)=TrMTrNTr(MN)g(M,N)=\textrm{Tr}M\,\textrm{Tr}N-\textrm{Tr}(MN). The evolution parameter is encoded in the Kraus operators and is not explicitly written in the above expressions.

V.1 Generalized noisy channel

In order to study the role of the 3-tangle in the teleportation fidelity, we focus on channels ΛB\Lambda_{B} whose Kraus operators have the following structure [32] (in the basis {|0=(1,0),|1=(0,1)}\{\left|{0}\right\rangle=(1,0)^{\top},\left|{1}\right\rangle=(0,1)^{\top}\} of B\mathcal{H}_{B})

K0=(1001p),K1=p(0cosζ0sinζ),K_{0}=\left(\begin{array}[]{cc}1&0\\ 0&\sqrt{1-p}\end{array}\right),\quad K_{1}=\sqrt{p}\left(\begin{array}[]{cc}0&\cos\zeta\\ 0&\sin\zeta\end{array}\right), (27)

where ζ[0,π/2]\zeta\in[0,\pi/2] and 0p10\leq p\leq 1. From the second equation in (12) and the explicit form of K0K_{0} in (27), we get 1|K0|1=10|UBEB|10=1p\left\langle{1}\right|K_{0}\left|{1}\right\rangle=\left\langle{10}\right|U_{BE_{B}}\left|{10}\right\rangle=\sqrt{1-p}. This means that under UBEBU_{BE_{B}}, the state |10\left|{10}\right\rangle transforms into a state that writes as

UBEB|10=1p|10+peiθ|10,U_{BE_{B}}\left|{10}\right\rangle=\sqrt{1-p}\left|{10}\right\rangle+\sqrt{p}e^{i\theta}\left|{10_{\perp}}\right\rangle, (28)

with |10\left|{10_{\perp}}\right\rangle a normalized state that is orthogonal to |10\left|{10}\right\rangle. Consequently, pp can be interpreted as the probability that the state |10BEB\left|{10}\right\rangle_{BE_{B}} evolves into an orthogonal state under the transformation UBEBU_{BE_{B}}. Clearly pp is a function of the evolution parameter of UBEBU_{BE_{B}}, typically the time for Hamiltonian evolutions UBEB(t)=eiHt/U_{BE_{B}}(t)=e^{-iHt/\hbar}. Initially (for UBEB(0)=𝖨4U_{BE_{B}}(0)=\mathsf{I}_{4}) pp vanishes, and increases up to p=1p=1 when the state |10\left|{10_{\perp}}\right\rangle (completely distinguishable from |10\left|{10}\right\rangle) is reached. This allows us to identify pp as a useful parameter to track the evolution induced by UBEBU_{BE_{B}}, without making specific assumptions regarding such unitary transformation.

Each value of the parameter ζ\zeta in (27) determines a specific channel, so comparison of the dynamics under different channels can be achieved by varying the values of ζ\zeta. When ζ=0\zeta=0, (27) reduce to the Kraus operators of the amplitude damping channel (ac), whereas for ζ=π/2\zeta=\pi/2 the Kraus operators of the dephasing channel (dc) are recovered. The ac and the dc are paradigmatic decoherence channels [47] that generate, respectively, w-type and ghz-type genuine entanglement in the 3-qubit system [43, 48, 49]. By means of Eq. (27) —corresponding to what we will call the generalized (noisy) channel, or gc for short— we can analyze intermediate situations lying between the ac and the dc, and particularly extend some of the results reported in, e.g., [16, 30, 31], to a wider range of channels.

When BB is subject to the gc, the evolved state (22) reads explicitly

|ψ=cosϕ|000+eiφsinϕ(1p|110++pcosζ|101+psinζ|111).\begin{split}|\psi\rangle=&\cos\phi\left|{000}\right\rangle+e^{i\varphi}\sin\phi\big{(}\sqrt{1-p}\left|{110}\right\rangle+\\ &+\sqrt{p}\cos\zeta\left|{101}\right\rangle+\sqrt{p}\sin\zeta\left|{111}\right\rangle\big{)}.\end{split} (29)

The entanglement of the corresponding resource state and the generated amount of 3-tangle become, using Eqs. (25) and (26),

CAB=01p,τABEB=02psin2ζ,C_{AB}=\mathcal{E}_{0}\sqrt{1-p},\quad\tau_{ABE_{B}}=\mathcal{E}^{2}_{0}\,p\sin^{2}\zeta, (30)

so only the 3-tangle depends on the specific channel. Further, psin2ζp\sin^{2}\zeta determines the fraction of the initial entanglement 0\mathcal{E}_{0} that can be converted into 3-partite entanglement. For fixed 0\mathcal{E}_{0} and pp, as ζ\zeta increases from 0 to π/2\pi/2 the 3-tangle goes from its minimum (0) to its maximum (02p\mathcal{E}^{2}_{0}\,p) value. The minimum corresponds to the ac case (ζ=0\zeta=0), and the maximum to the dc case (ζ=π/2\zeta=\pi/2), which is the only channel for which all the initial entanglement can be transformed into 3-tangle (at p=1p=1).

From Eq. (14) it follows that for BB subject to the gc, the maximal average fidelity can be written as Fmax=maxi{FΦi}F_{\max}=\max_{\,i}\,\{F_{\Phi_{i}}\}, with

FΦi=13+23β|Φi|(𝖨2Kβ)|ϕ0|2.F_{\Phi_{i}}=\frac{1}{3}+\frac{2}{3}\sum_{\beta}\big{|}\langle\Phi_{i}\big{|}(\mathsf{I}_{2}\otimes K_{\beta})|\phi_{0}\rangle\big{|}^{2}. (31)

Direct calculation gives

FΦ±\displaystyle F_{\Phi^{\pm}}\!\! =\displaystyle= 23+13(±01pcosφ𝒫1pcos2ζ),\displaystyle\!\frac{2}{3}+\frac{1}{3}\Big{(}\!\pm\mathcal{E}_{0}\sqrt{1-p}\cos\varphi-\mathcal{P}_{1}p\cos^{2}\zeta\Big{)}, (32a)
FΨ+\displaystyle F_{\Psi^{+}}\!\! =\displaystyle= FΨ=2313(1𝒫1pcos2ζ),\displaystyle\!F_{\Psi^{-}}=\frac{2}{3}-\frac{1}{3}\Big{(}1-\mathcal{P}_{1}p\cos^{2}\zeta\Big{)}, (32b)

where

𝒫1=sin2ϕ\mathcal{P}_{1}=\sin^{2}\!\phi (33)

stands for the initial population of the state |11AB\left|{11}\right\rangle_{AB}. On one hand 𝒫1pcos2ζ1\mathcal{P}_{1}p\cos^{2}\zeta\leq 1 implies that FΨ±2/3F_{\Psi^{\pm}}\leq 2/3. On the other hand,

max{FΦ+,FΦ}\displaystyle\max\{F_{\Phi^{+}},F_{\Phi^{-}}\} =\displaystyle= {FΦ+φ[π/2,π/2]FΦφ[π/2,3π/2]\displaystyle\begin{cases}F_{\Phi^{+}}&\varphi\in[-\pi/2,\pi/2]\\ F_{\Phi^{-}}&\varphi\in[\pi/2,3\pi/2]\end{cases}
=\displaystyle= 23+13(01p|cosφ|𝒫1pcos2ζ).\displaystyle\frac{2}{3}+\frac{1}{3}\Big{(}\mathcal{E}_{0}\sqrt{1-p}\,|\!\cos\varphi|-\mathcal{P}_{1}p\cos^{2}\zeta\Big{)}.

Consequently, under the generalized channel, the quantity of interest (6) is given by

𝔽maxgc=23+13max{0,[CAB(p)|cosφ|𝒫1pcos2ζ]},\mathbb{F}^{\textsc{gc}}_{\max}=\frac{2}{3}+\frac{1}{3}\max\Big{\{}0,\Big{[}C_{AB}(p)\,|\!\cos\varphi|-\mathcal{P}_{1}p\cos^{2}\zeta\Big{]}\Big{\}}, (35)

where we have used Eq. (30) for CABC_{AB}.

Refer to caption
Figure 5: Evolution of the entanglement of the resource state (left panel), the maximal average fidelity above the classical threshold value (central panel), and the 3-tangle (right panel) when the qubit BB is subject to a gc channel, for different values of the channel parameter ζ[0,π/2]\zeta\in[0,\pi/2], and ϕ=π/4\phi=\pi/4, φ=0\varphi=0. The channel induces a loss in the resource entanglement CABC_{AB}, and the teleportation fidelity decreases accordingly. However, when different channels (values of ζ\zeta) are considered, it is observed that those that generate more 3-tangle favor higher fidelities at each pp. Therefore, as ζ\zeta increases from 0 (amplitude damping channel) to π/2\pi/2 (dephasing channel), the generated 3-tangle helps to enhance the maximal average fidelity (provided CAB0)C_{AB}\neq 0).

It follows from Eq. (35) that, as in the noiseless scenario (see Eq. (17)), when the maximal average fidelity exceeds the classical threshold it has a non-negative contribution proportional to the entanglement of the resource state, attenuated by |cosφ||\!\cos\varphi|. This indicates that, as in the noiseless case, the optimal φ\varphi is 0,π0,\pi. In the gc case, however, an additional negative contribution (a loss in the fidelity) appears that depends on the channel (via cos2ζ\cos^{2}\zeta), the evolution (pp), and the initial excited population (𝒫1\mathcal{P}_{1}). For fixed 𝒫1\mathcal{P}_{1} and pp, such loss decreases as τABEBsin2ζ\tau_{ABE_{B}}\sim\sin^{2}\zeta increases. This means that given an initial state, and at each stage of the evolution (determined by a fixed pp), the channels that produce higher amounts of 3-tangle lead to higher values of the maximum average fidelity.

The above conclusion can also be extracted from Figure 5, which shows CABC_{AB}, 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max}, and τABEB\tau_{ABE_{B}} for different values of ζ\zeta and varying pp, considering a maximally entangled initial state (ϕ=π/4\phi=\pi/4) and setting φ\varphi to its optimal value φ=0\varphi=0. At p=0p=0, 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} reduces to 𝔽maxnon-int\mathbb{F}^{\textrm{non-int}}_{\max}, no 3-tangle exists, and the success of the teleportation is ascribable to the (bipartite) entanglement 0=CAB(0)\mathcal{E}_{0}=C_{AB}(0) only, in line with the discussion in Sec. IV. As pp increases, the action of ΛB\Lambda_{B} degrades the entanglement of the resource state (whose dynamics is independent of ζ\zeta), and a concomitant gradual loss in the maximal average fidelity is observed for all ζ\zeta. This is an expected behavior under generic local channels ΛAΛB\Lambda_{A}\otimes\Lambda_{B} (meaning non-increasing entanglement operations).

In its turn, as pp increases the 3-tangle increases as well, as a result of the redistribution of the entanglement. Therefore, for any given generalized noisy channel, the decay of the maximal average fidelity along the evolution is accompanied by an increase of the 3-tangle. The question then arises as to under which channels 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} is more robust against the noise, that is, under which channels the loss in the teleportation success is reduced at any stage of the evolution. To answer it, we compare the values of 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} for all possible gc channels while keeping pp fixed. Comparison of the central and rightmost panels in Fig. 5 shows that the decaying path of 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} varies with ζ\zeta in such a way that for each and every pp, the channels that generate higher amounts of 3-tangle yield higher fidelities. Table 1 exemplifies this by explicitly showing the values of 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} and τABEB\tau_{ABE_{B}} extracted from Fig. 5 for p=0.8p=0.8. While the entanglement of the resource state is constant as ζ\zeta varies (see the first equation in (30)), the maximal average fidelity and the 3-tangle do change in a correlated fashion: an increase in τABEB\tau_{ABE_{B}} is accompanied by an increase in 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max}. This behavior replicates for any other value of p(0,1)p\in(0,1).

ζ\zeta CABC_{AB} 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} τABEB\tau_{ABE_{B}}
0 0.447214 0.682405 0
π/6\pi/6 0.447214 0.715738 0.2
π/4\pi/4 0.447214 0.749071 0.4
π/3\pi/3 0.447214 0.782405 0.6
π/2\pi/2 0.447214 0.815738 0.8
Table 1: Values of CABC_{AB}, 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max}, and τABEB\tau_{ABE_{B}} taken from Fig. 5 (corresponding to ϕ=π/4\phi=\pi/4 and φ=0\varphi=0) for p=0.8p=0.8 and different values of the channel parameter ζ\zeta.

Therefore, despite the adverse influence of noise on the fidelity, the detrimental effects are lessened under channels that give rise to higher amounts of tripartite entanglement. It is in this sense that the 3-tangle improves the teleportation success, and may help to maintain 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} above the classical threshold value longer (i.e., for larger values of pp, as seen in the central panel of Fig. 5).

It is also clear from Fig. 5 that the channel that gives better fidelities throughout the evolution is the dc, corresponding to ζ=π/2\zeta=\pi/2. Only for this channel 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} exceeds the classically attainable value for all 0<p<10<p<1. This holds not only for a maximally entangled initial state, but also for 00\mathcal{E}_{0}\neq 0. In fact, it follows from Eq. (35) that the condition 𝔽maxgc>2/3\mathbb{F}^{\textsc{gc}}_{\max}>2/3 amounts to

𝒫1pcos2ζ<CAB(p)|cosφ|=01p|cosφ|,\mathcal{P}_{1}p\cos^{2}\zeta<C_{AB}(p)|\!\cos\varphi|=\mathcal{E}_{0}\sqrt{1-p}\,|\!\cos\varphi|, (36)

which is trivially satisfied for p(0,1)p\in(0,1) and 00\mathcal{E}_{0}\neq 0 taking ζ=π/2\zeta=\pi/2. Further, by writing 0=2𝒫1(1𝒫1)\mathcal{E}_{0}=2\sqrt{\mathcal{P}_{1}(1-\mathcal{P}_{1})}, the condition (36) rewrites as

12cos2ζ|cosφ|p1p<1𝒫1𝒫1,\frac{1}{2}\frac{\cos^{2}\zeta}{|\!\cos\varphi|}\frac{p}{\sqrt{1-p}}<\sqrt{\frac{1-\mathcal{P}_{1}}{\mathcal{P}_{1}}}, (37)

which exhibits the role of the initial population 𝒫1\mathcal{P}_{1} in the dynamics of 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} (assuming fixed φ\varphi and ζ\zeta): the left-hand side of the inequality is an increasing function of pp, whereas the upper bound decreases with 𝒫1\mathcal{P}_{1}; consequently, as 𝒫1\mathcal{P}_{1} increases the inequality becomes more restrictive, and a point in the evolution is reached sooner (for lower values of pp) at which (37) does not longer hold, and the fidelity drops below the threshold value 2/32/3.

Refer to caption
Figure 6: Trajectories in the space (τABEB,𝔽maxgc)(\tau_{ABE_{B}},\mathbb{F}^{\textsc{gc}}_{\max}) as the channel parameter ζ\zeta runs from 0 (black circles) to π/2\pi/2 (red stars), for three initial states with ϕ=π/6\phi=\pi/6 (red dashed line), ϕ=π/4\phi=\pi/4 (black solid line), and ϕ=π/3\phi=\pi/3 (blue dotted-dashed line), at three stages of the evolution: p=0.25p=0.25 (left), p=0.5p=0.5 (center), and p=0.75p=0.75 (right). In all cases the relative phase was fixed to its optimal value φ=0\varphi=0. For each initial state, and at fixed pp, the channels that produce higher 3-tangle lead to higher maximal average fidelities. A non-trivial role of the initial entanglement 0\mathcal{E}_{0} and the initial excited population 𝒫1\mathcal{P}_{1} is observed (see text).

Finally, as p1p\rightarrow 1 we have CAB0C_{AB}\rightarrow 0, and from Eq. (35) 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} cannot exceed 2/32/3 despite τABEB\tau_{ABE_{B}} attains its maximum value. This evinces that the even though the presence of tripartite entanglement improves the teleportation success (in the sense described above), a non-zero entanglement of the resource state is key to trigger the potential of the 3-tangle to assist the teleportation.

Figure 6 shows, at different stages of the evolution, curves in the space (τABEB,𝔽maxgc)(\tau_{ABE_{B}},\mathbb{F}^{\textsc{gc}}_{\max}) sweeping as the channel parameter goes from ζ=0\zeta=0 (ac, black circles) to ζ=π/2\zeta=\pi/2 (dc, red stars). Each trajectory corresponds to an initial state determined by ϕ=π/6\phi=\pi/6 (red dashed line), ϕ=π/4\phi=\pi/4 (black solid line), and ϕ=π/3\phi=\pi/3 (blue dotted-dashed line), all with φ=0\varphi=0. In all curves the relation between the 3-tangle and the maximal average fidelity discussed above is manifest: throughout the evolution, better fidelities are attained under channels that produce higher amounts of 3-tangle.

Interestingly, Fig. 6 shows that whereas an initial resource state with maximal entanglement (ϕ=π/4\phi=\pi/4) leads to the highest value of 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} (for sufficiently large values of τABEB\tau_{ABE_{B}}), a maximally entangled state is not always the optimal one for achieving a better fidelity. This is clearly seen, for example, in the rightmost panel of Fig. 6, where the fidelity for ϕ=π/6\phi=\pi/6 (00.866\mathcal{E}_{0}\sim 0.866) is greater than that corresponding to ϕ=π/4\phi=\pi/4 (0=1\mathcal{E}_{0}=1) for some channels. Moreover, noticing that ϕ=π/3\phi=\pi/3 gives the same 0\mathcal{E}_{0} as ϕ=π/6\phi=\pi/6, the difference between the red dashed and the blue dotted-dashed curves brings out the effect of the initial population 𝒫1\mathcal{P}_{1} on the behavior of 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max}. In particular, for ϕ=π/3\phi=\pi/3 we have 𝒫1=0.75\mathcal{P}_{1}=0.75, while for ϕ=π/6\phi=\pi/6 we have 𝒫1=0.25\mathcal{P}_{1}=0.25. As seen from Eq. (35), the greater 𝒫1\mathcal{P}_{1}, the greater the loss in the fidelity (which goes as 𝒫1pcos2ζ-\mathcal{P}_{1}p\cos^{2}\zeta), which explains why the blue dotted-dashed line runs below the red dashed one, despite both curves correspond to the same initial entanglement.

In summary, from all the generalized channels acting on BB, those that improve the quantum teleportation success (at each instant) are the ones that produce higher amounts of 3-tangle among the qubits A,B,EBA,B,E_{B}. The optimal channel corresponds to ζ=π/2\zeta=\pi/2 (dc), and the worst to ζ=0\zeta=0 (ac). This can be seen graphically, and is corroborated by Eq. (35), which gives

𝔽maxac𝔽maxgc𝔽maxdc,\mathbb{F}^{\textsc{ac}}_{\max}\leq\mathbb{F}^{\textsc{gc}}_{\max}\leq\mathbb{F}^{\textsc{dc}}_{\max}, (38)

with

𝔽maxdc=23+13CAB(p)|cosφ|,\mathbb{F}^{\textsc{dc}}_{\max}=\frac{2}{3}+\frac{1}{3}\,C_{AB}(p)\,|\!\cos\varphi|, (39)

and

𝔽maxac\displaystyle\mathbb{F}^{\textsc{ac}}_{\max} =\displaystyle= 23+13max{0,CAB(p)|cosφ|𝒫1p}\displaystyle\frac{2}{3}+\frac{1}{3}\max\Big{\{}0,C_{AB}(p)\,|\!\cos\varphi|-\mathcal{P}_{1}p\Big{\}} (40)
=\displaystyle= max{23,𝔽maxdc13𝒫1p}.\displaystyle\max\Big{\{}\frac{2}{3},\mathbb{F}^{\textsc{dc}}_{\max}-\frac{1}{3}\mathcal{P}_{1}p\Big{\}}.

Below we explore in more detail the states produced by the limiting cases dc and ac.

V.2 Amplitude Damping Channel vs Dephasing Channel

The dc produces ghz-type states, whose form follows from Eq. (29) with ζ=π/2\zeta=\pi/2:

|ψghz\displaystyle\left|{\psi_{\textsc{ghz}}}\right\rangle \displaystyle\equiv |ψ(ζ=π/2)\displaystyle\left|{\psi(\zeta=\pi/2)}\right\rangle
=\displaystyle= cosϕ|000+eiφsinϕ(1p|110+\displaystyle\cos\phi\left|{000}\right\rangle+e^{i\varphi}\sin\phi\big{(}\sqrt{1-p}\left|{110}\right\rangle+
+p|111).\displaystyle+\sqrt{p}\left|{111}\right\rangle\big{)}.

If ϕ=π/4\phi=\pi/4 and φ=0\varphi=0 then, at p=1p=1, the usual ghz state |ghz=12(|000+|111)\left|{\textsc{ghz}}\right\rangle=\frac{1}{\sqrt{2}}\big{(}\!\left|{000}\right\rangle+\left|{111}\right\rangle\!\big{)} is reached, having null qubit-qubit entanglement and τABEB=1.\tau_{ABE_{B}}=1.

For ζ=0\zeta=0 the gc channel reduces to the ac, and states that are equivalent (up to local unitary transformations) to w-type states arise. This can be seen by putting ζ=0\zeta=0 in Eq. (29), obtaining

|ψw\displaystyle\left|{\psi_{\textsc{w}}}\right\rangle \displaystyle\equiv |ψ(ζ=0)\displaystyle\left|{\psi(\zeta=0)}\right\rangle
=\displaystyle= cosϕ|000+eiφsinϕ(1p|110+\displaystyle\cos\phi\left|{000}\right\rangle+e^{i\varphi}\sin\phi\big{(}\sqrt{1-p}\left|{110}\right\rangle+
+p|101).\displaystyle+\sqrt{p}\left|{101}\right\rangle\big{)}.

For ϕ=arccos(1/3)\phi=\arccos(1/\sqrt{3}) and φ=0\varphi=0 this state becomes, at p=1/2p=1/2, equivalent (up to a local unitary transformation) to the usual w state |w=13(|000+|110+|101)\left|{\textsc{w}}\right\rangle=\frac{1}{\sqrt{3}}\big{(}\left|{000}\right\rangle+\left|{110}\right\rangle+\left|{101}\right\rangle\big{)}, characterized by having all qubit-qubit entanglements equal to Cij=2/3C_{ij}=2/3, and null 3-tangle.

Refer to caption
Figure 7: Evolution of the maximal average fidelity (above the threshold value 2/32/3), employing initial states that evolve towards the states |ghz\left|{\textsc{ghz}}\right\rangle (red solid line) and |w\left|{\textsc{w}}\right\rangle (blue dashed line) , under the dc and the ac channel, respectively.

Figure 7 depicts the evolution of 𝔽max\mathbb{F}_{\max} for the initial states |ψghz|ϕ=π4,φ=0\left|{\psi_{\textsc{ghz}}}\right\rangle|_{\phi=\frac{\pi}{4},\varphi=0} (red solid line), and |ψw|ϕ=arccos13,φ=0\left|{\psi_{\textsc{w}}}\right\rangle|_{\phi=\arccos\frac{1}{\sqrt{3}},\varphi=0} (blue dashed line). It shows that the maximal average fidelity improves when the resource state involves two qubits from a 3-qubit system that evolves towards the state |ghz\left|{\textsc{ghz}}\right\rangle, rather than to |w\left|{\textsc{w}}\right\rangle.

VI Fidelity and four-partite entanglement generated via noisy channels

We now consider that both qubits AA and BB undergo local channels ΛA\Lambda_{A} and ΛB\Lambda_{B}, as a result of their separate interaction with initially uncorrelated qubits, EAE_{A} and EBE_{B}. The initial resource state is thus (see Eqs. (1) and (11))

|ψ0=cosϕ|0000+eiφsinϕ|1100,\left|{\psi_{0}}\right\rangle=\cos\phi\left|{0000}\right\rangle+e^{i\varphi}\sin\phi\left|{1100}\right\rangle, (43)

and the Kraus operators of the local channels are generically given by Eq. (12). Figure 8 shows the first two stages of the corresponding teleportation protocol.

Refer to caption
Figure 8: Qubits AA, BB, EAE_{A} and EBE_{B} are prepared in the state (43). The pairs A,EAA,E_{A}, and B,EBB,E_{B} interact via unitary operations, thus AA and BB undergo local channels whose corresponding Kraus operators are {Qα}\{Q_{\alpha}\} and {Kβ}\{K_{\beta}\}.

Different forms of multipartite entanglement may arise in this 4-qubit system [50]. Here we will focus on the multipartite entanglement as measured by the 4-tangle τ4\tau_{4}, defined for a 4-qubit pure state |ψ\left|{\psi}\right\rangle as [51]

τ4=|ψ|σy4|ψ|2.\tau_{4}=|\left\langle{\psi}\right|\sigma^{\otimes 4}_{y}\left|{\psi^{*}}\right\rangle|^{2}. (44)

This quantity becomes maximal (τ4=1\tau_{4}=1) for the 4-partite ghz state 12(|0000+|1111)\frac{1}{\sqrt{2}}(\left|{0000}\right\rangle+\left|{1111}\right\rangle), and vanishes for the 4-partite w state 14(|1000+|0100+|0010+|0001)\frac{1}{\sqrt{4}}(\left|{1000}\right\rangle+\left|{0100}\right\rangle+\left|{0010}\right\rangle+\left|{0001}\right\rangle).

VI.1 Two local generalized noisy channels

Assuming that ΛA\Lambda_{A} and ΛB\Lambda_{B} are generalized noisy channels, the corresponding Kraus operators read, following (27),

Q0=(1001pA),Q1=pA(0cosζA0sinζA),Q_{0}=\left(\begin{array}[]{cc}1&0\\ 0&\sqrt{1-p_{A}}\end{array}\right),\;Q_{1}=\sqrt{p_{A}}\left(\begin{array}[]{cc}0&\cos\zeta_{A}\\ 0&\sin\zeta_{A}\end{array}\right), (45a)
and
K0=(1001pB),K1=pB(0cosζB0sinζB),K_{0}=\left(\begin{array}[]{cc}1&0\\ 0&\sqrt{1-p_{B}}\end{array}\right),\;K_{1}=\sqrt{p_{B}}\left(\begin{array}[]{cc}0&\cos\zeta_{B}\\ 0&\sin\zeta_{B}\end{array}\right), (45b)

with ζA(B)[0,π/2]\zeta_{A(B)}\in[0,\pi/2] the parameters that determine the specific channels and 0pA(B)10\leq p_{A(B)}\leq 1, again with pA(B)p_{A(B)} appropriate parametrizations of the evolution parameter of the transformation UAEA(BEB)U_{AE_{A}(BE_{B})}.

Resorting to Eq. (13) the initial state (43) evolves into

|ψ=ABEAEB(cosϕ|00+eiφsinϕqAqB|11)|00+eiφsinϕ{qApB(cosζB|10+sinζB|11)|01+pAqB(cosζA|01+sinζA|11)|10+pApB(cosζA|0+sinζA|1)(cosζB|0+sinζB|1)|11},\begin{split}\left|{\psi}\right\rangle&{}_{ABE_{A}E_{B}}=\Bigl{(}\cos\phi\left|{00}\right\rangle+e^{i\varphi}\sin\phi\sqrt{q_{A}q_{B}}\left|{11}\right\rangle\Bigr{)}\left|{00}\right\rangle+e^{i\varphi}\sin\phi\Big{\{}\sqrt{q_{A}p_{B}}\Bigl{(}\cos\zeta_{B}\left|{10}\right\rangle+\sin\zeta_{B}\left|{11}\right\rangle\Bigr{)}\left|{01}\right\rangle\\ +&\sqrt{p_{A}q_{B}}\Bigl{(}\cos\zeta_{A}\left|{01}\right\rangle+\sin\zeta_{A}\left|{11}\right\rangle\Bigr{)}\left|{10}\right\rangle+\sqrt{p_{A}p_{B}}\Bigl{(}\cos\zeta_{A}\left|{0}\right\rangle+\sin\zeta_{A}\left|{1}\right\rangle\Bigl{)}\Bigl{(}\cos\zeta_{B}\left|{0}\right\rangle+\sin\zeta_{B}\left|{1}\right\rangle\Bigr{)}\left|{11}\right\rangle\Big{\}},\end{split} (46)

where we wrote qA(B)=1pA(B)q_{A(B)}=1-p_{A(B)}. From this state direct calculation gives

τ4\displaystyle\tau_{4} =\displaystyle= pApB|0sinζAsinζB\displaystyle p_{A}p_{B}\big{|}\mathcal{E}_{0}\sin\zeta_{A}\sin\zeta_{B} (47)
+ 4eiφ𝒫1qAqBcosζAcosζB|2.\displaystyle+\,4e^{i\varphi}\mathcal{P}_{1}\sqrt{q_{A}q_{B}}\cos\zeta_{A}\cos\zeta_{B}\big{|}^{2}.

Unlike the previous case in which only BB passes through the generalized channel, where CABC_{AB} was independent of the channel parameter ζ\zeta (see Eq. (30)), in the present scenario the entanglement of the resource state depends on both ζA\zeta_{A} and ζB\zeta_{B}. However, the general and explicit dependence of CABC_{AB} on ζA,ζB\zeta_{A},\zeta_{B} is far from trivial and will be omitted here. It suffices to recall that CABC_{AB} can be obtained from Eq. (19) using the resource state ρAB=TrEAEB|ψψ|\rho_{AB}=\textrm{Tr}_{E_{A}E_{B}}\left|{\psi}\right\rangle\!\left\langle{\psi}\right|, with |ψ\left|{\psi}\right\rangle given by (46).

From Eq. (10) the maximal average fidelity for two gc channels is Fmax=maxi{FΦi}F_{\max}=\max_{\,i}\,\{F_{\Phi_{i}}\} with

FΦi=13+23{αβ|Φi|QαKβ|ϕ0|2}.\begin{split}F_{\Phi_{i}}&=\frac{1}{3}+\frac{2}{3}\Big{\{}\sum_{\alpha\beta}\big{|}\langle\Phi_{i}\big{|}Q_{\alpha}\otimes K_{\beta}\big{|}\phi_{0}\rangle\big{|}^{2}\Big{\}}.\end{split} (48)

Substitution of the Kraus operators (45) gives

FΦ±\displaystyle F_{\Phi^{\pm}} =\displaystyle= 23+13(±0qAqBcosφ𝒫1Δ),\displaystyle\frac{2}{3}+\frac{1}{3}\Big{(}\pm\mathcal{E}_{0}\sqrt{q_{A}q_{B}}\cos{\varphi}-\mathcal{P}_{1}\Delta_{\mp}\Big{)}, (49a)
FΨ±\displaystyle F_{\Psi^{\pm}} =\displaystyle= 2313(1𝒫1Δ±),\displaystyle\frac{2}{3}-\frac{1}{3}\Big{(}1-\mathcal{P}_{1}\Delta_{\pm}\Big{)}, (49b)

where Δ±[0,1]\Delta_{\pm}\in[0,1] and is given by

Δ±\displaystyle\Delta_{\pm} =\displaystyle= qBpAcos2ζA+qApBcos2ζB\displaystyle q_{B}p_{A}\cos^{2}\zeta_{A}+q_{A}p_{B}\cos^{2}\zeta_{B} (50)
+pApBsin2(ζA±ζB).\displaystyle+\,p_{A}p_{B}\sin^{2}(\zeta_{A}\pm\zeta_{B}).

Consequently FΨ±2/3F_{\Psi^{\pm}}\leq 2/3; further, since ζA,ζB[0,π/2]\zeta_{A},\zeta_{B}\in[0,\pi/2] it holds that sin2(ζA+ζB)sin2(ζAζB)\sin^{2}(\zeta_{A}+\zeta_{B})\geq\sin^{2}(\zeta_{A}-\zeta_{B}), whence Δ+Δ\Delta_{+}\geq\Delta_{-}. This in turn leads us to

cosφ0FΦ+FΦ,\cos\varphi\geq 0\Rightarrow F_{\Phi^{+}}\geq F_{\Phi^{-}}, (51)

so in the present case 𝔽max\mathbb{F}_{\max}, given by Eq. (6), writes as

𝔽max±gc/gc=23+13max{0,[0qAqB|cosφ|𝒫1Δ]}\mathbb{F}^{\textsc{gc/gc}}_{\max\pm}=\frac{2}{3}+\frac{1}{3}\max\Big{\{}0,\Big{[}\mathcal{E}_{0}\sqrt{q_{A}q_{B}}\,|\!\cos{\varphi}|-\mathcal{P}_{1}\Delta_{\mp}\Big{]}\Big{\}} (52)

where the upper/lower sign must be chosen accordingly with

cosφ14𝒫11𝒫11qAqB(Δ+Δ).\cos\varphi\gtrless-\frac{1}{4}\frac{\sqrt{\mathcal{P}_{1}}}{\sqrt{1-\mathcal{P}_{1}}}\frac{1}{\sqrt{q_{A}q_{B}}}(\Delta_{+}-\Delta_{-}). (53)

As follows from (52), the optimal phase φ\varphi is the same as in the previous cases, namely φ=0,π\varphi=0,\pi. Also, a loss in the fidelity emerges, encoded in the negative term proportional to 𝒫1\mathcal{P}_{1}.

VI.2 Fidelity and multipartite entanglement for 4-qubits under parallel generalized noisy channels

In our forthcoming analysis, we will simplify the equations taking pA=pB=pp_{A}=p_{B}=p 333This means that the evolution parameters of UAEAU_{AE_{A}} and UBEBU_{BE_{B}} are the same (for example they are both the time tt), and the parametrizations are identical, so p(t)p(t) is the same for both pairs of Kraus operators (45).. Equation (47) thus reads

τ4\displaystyle\tau_{4} =\displaystyle= p2[0sinζAsinζB+4𝒫1(1p)cosζAcosζB]2\displaystyle p^{2}\Big{[}\mathcal{E}_{0}\sin\zeta_{A}\sin\zeta_{B}+4\mathcal{P}_{1}(1-p)\cos\zeta_{A}\cos\zeta_{B}\Big{]}^{2} (54)
20𝒫1p2(1p)sin2ζAsin2ζB(1cosφ).\displaystyle-2\mathcal{E}_{0}\mathcal{P}_{1}p^{2}(1-p)\sin 2\zeta_{A}\sin 2\zeta_{B}(1-\cos\varphi).

Recalling that ζA,B[0,π/2]\zeta_{A,B}\in[0,\pi/2], the second line in Eq. (54) is a non-positive term that reduces the amount of 4-partite entanglement. Notice, however, that for φ=0\varphi=0 such term vanishes, and consequently the relative phase that maximizes τ4\tau_{4} maximizes also 𝔽max±gc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max\pm}. By choosing this optimal phase, Eq. (53) holds with the upper inequality sign, and accordingly Eq. (52) reduces to

𝔽maxgc/gc\displaystyle\mathbb{F}^{\textsc{gc/gc}}_{\max} =\displaystyle= 23+13max{0,[0(1p)𝒫1p[psin2(ζAζB)\displaystyle\frac{2}{3}+\frac{1}{3}\max\Big{\{}0,\Big{[}\mathcal{E}_{0}(1-p)-\mathcal{P}_{1}p[p\sin^{2}(\zeta_{A}-\zeta_{B}) (55)
+(1p)(cos2ζA+cos2ζB)]]},\displaystyle+(1-p)(\cos^{2}\zeta_{A}+\cos^{2}\zeta_{B})]\Big{]}\Big{\}},

whereas τ4\tau_{4} becomes

τ4=[0psinζAsinζB+4𝒫1p(1p)cosζAcosζB]2.\tau_{4}=\Big{[}\mathcal{E}_{0}p\sin\zeta_{A}\sin\zeta_{B}+4\mathcal{P}_{1}p(1-p)\cos\zeta_{A}\cos\zeta_{B}\Big{]}^{2}. (56)

If 0p4𝒫1p(1p)\mathcal{E}_{0}p\geq 4\mathcal{P}_{1}p(1-p) then

0pcosζAcosζB4𝒫1p(1p)cosζAcosζB,\mathcal{E}_{0}p\cos\zeta_{A}\cos\zeta_{B}\geq 4\mathcal{P}_{1}p(1-p)\cos\zeta_{A}\cos\zeta_{B}, (57)

whence, by adding 0psinζAsinζB\mathcal{E}_{0}p\sin\zeta_{A}\sin\zeta_{B} on both sides of the inequality we get

0pcos(ζAζB)τ4.\mathcal{E}_{0}p\cos(\zeta_{A}-\zeta_{B})\geq\sqrt{\tau_{4}}. (58)

An analogous reasoning applies for 0p4𝒫1p(1p)\mathcal{E}_{0}p\leq 4\mathcal{P}_{1}p(1-p), so we finally get

τ4T2cos2(ζAζB),\tau_{4}\leq T^{2}\cos^{2}(\zeta_{A}-\zeta_{B}), (59)

where T=max{0p,4𝒫1p(1p)}T=\max\{\mathcal{E}_{0}p,4\mathcal{P}_{1}p(1-p)\}. This means that the maximum value of τ4\tau_{4} over all possible channels is located along the line ζA=ζB=ζ\zeta_{A}=\zeta_{B}=\zeta, and explicitly reads

τ4max={02p2=τ4|ζ=π/2if04𝒫1(1p),16𝒫12p2(1p)2=τ4|ζ=0if04𝒫1(1p).\tau^{\max}_{4}=\begin{cases}\mathcal{E}^{2}_{0}p^{2}=\tau_{4}|_{\zeta=\pi/2}\!&\textrm{if}\;\,\mathcal{E}_{0}\geq 4\mathcal{P}_{1}(1-p),\\ 16\mathcal{P}^{2}_{1}p^{2}(1-p)^{2}=\tau_{4}|_{\zeta=0}\!&\textrm{if}\;\,\mathcal{E}_{0}\leq 4\mathcal{P}_{1}(1-p).\end{cases} (60)

The maximum 4-partite entanglement is thus reached whenever both channels are either dc or ac. In the former case also 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} reaches its maximum value, over all gc/gc channels. This can be seen by noticing that the loss in the fidelity, encoded in the (non-positive) term proportional to 𝒫1p\mathcal{P}_{1}p in Eq. (55), vanishes irrespective of the initial state and the evolution parameter only for ζA=ζB=π/2\zeta_{A}=\zeta_{B}=\pi/2, resulting in

𝔽maxgc/gc𝔽maxdc/dc\displaystyle\mathbb{F}^{\textsc{gc/gc}}_{\max}\leq\mathbb{F}^{\textsc{dc/dc}}_{\max} =\displaystyle= 23+130(1p)\displaystyle\frac{2}{3}+\frac{1}{3}\mathcal{E}_{0}(1-p)
=\displaystyle= 23+13CABdc/dc(p),\displaystyle\frac{2}{3}+\frac{1}{3}C^{\textsc{dc/dc}}_{AB}(p),

where CABdc/dc=0(1p)C^{\textsc{dc/dc}}_{AB}=\mathcal{E}_{0}(1-p) stands for the entanglement of the resource state when two parallel dc channels are implemented. Consequently, whenever the condition in the first line of (60) holds, i.e., whenever

1121𝒫1𝒫1p,1-\frac{1}{2}\sqrt{\frac{1-\mathcal{P}_{1}}{\mathcal{P}_{1}}}\leq p, (62)

the channel that maximizes τ4\tau_{4} maximizes also 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max}. Notice that (62) is satisfied for all pp provided 𝒫11/5\mathcal{P}_{1}\leq 1/5. That is, for sufficiently low populations 𝒫1\mathcal{P}_{1}, the maximal value of τ4\tau_{4} and of 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} are jointly reached.

If the condition (62) is not met, τ4max\tau^{\max}_{4} is given by the second line in Eq. (60) and corresponds to ζA=ζB=0\zeta_{A}=\zeta_{B}=0, i.e., to the combination ac/ac. In this case (55) reduces to

𝔽maxac/ac\displaystyle\mathbb{F}^{\textsc{ac/ac}}_{\max} =\displaystyle= 23+13CABac/ac(p)\displaystyle\frac{2}{3}+\frac{1}{3}C^{\textsc{ac/ac}}_{AB}(p)
=\displaystyle= max{23,𝔽maxdc/dc23𝒫1p(1p)},\displaystyle\max\Big{\{}\frac{2}{3},\mathbb{F}^{\textsc{dc/dc}}_{\max}-\frac{2}{3}\mathcal{P}_{1}p(1-p)\Big{\}},

where in the first line we used that under parallel ac, CABac/ac=max{0,(1p)(02𝒫1p)}C^{\textsc{ac/ac}}_{AB}=\max\{0,(1-p)(\mathcal{E}_{0}-2\mathcal{P}_{1}p)\}.

Now, from Eq. (55) we obtain the following expression for 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} along the line ζA=ζB=ζ\zeta_{A}=\zeta_{B}=\zeta:

𝔽maxgc/gc(ζ)\displaystyle\mathbb{F}^{\textsc{gc/gc}}_{\max}(\zeta) =\displaystyle= 23+13max{0,[0(1p)\displaystyle\frac{2}{3}+\frac{1}{3}\max\Big{\{}0,\Big{[}\mathcal{E}_{0}(1-p) (64)
2𝒫1p(1p)cos2ζ]},\displaystyle-2\mathcal{P}_{1}p(1-p)\cos^{2}\zeta\Big{]}\Big{\}},

whose minimum above the threshold value 2/32/3 corresponds to the channel ac/ac. This observation, together with the fact that τ4max\tau^{\max}_{4} is attained at ζA=ζB=ζ=0,π/2\zeta_{A}=\zeta_{B}=\zeta=0,\pi/2, indicates that the channels that maximize the 4-partite entanglement are not always those that maximize 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max}. Rather, 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} is either maximal or minimal (within the family of twin channels with ζA=ζB=ζ\zeta_{A}=\zeta_{B}=\zeta) on points where the 4-tangle is maximal.

As for the vanishing value of τ4\tau_{4}, we see from Eq. (56) that it corresponds to the combination of a dc on either one of the qubits, and an ac on the other one. This null 4-tangle is accompanied by the following maximal average fidelity

𝔽maxac/dc=𝔽maxdc/ac\displaystyle\mathbb{F}^{\textsc{ac/dc}}_{\max}=\mathbb{F}^{\textsc{dc/ac}}_{\max} =\displaystyle= 23+13max{0,0(1p)𝒫1p}\displaystyle\frac{2}{3}+\frac{1}{3}\max\Big{\{}0,\mathcal{E}_{0}(1-p)-\mathcal{P}_{1}p\Big{\}} (65)
=\displaystyle= max{23,𝔽maxdc/dc13𝒫1p},\displaystyle\max\Big{\{}\frac{2}{3},\mathbb{F}^{\textsc{dc/dc}}_{\max}-\frac{1}{3}\mathcal{P}_{1}p\Big{\}},

an expression that is analogous to Eq. (40) for the case with vanishing 3-tangle (ac).

Refer to caption
Figure 9: 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} versus τ4\tau_{4} for different values of ϕ\phi (rows, from top to bottom: ϕ=π/6,π/4,π/3\phi=\pi/6,\pi/4,\pi/3) and pp (columns, from left to right: p=0.25,0.5,0.75p=0.25,0.5,0.75). The regions are generated as all the channels are considered (by varying ζA\zeta_{A}, ζB\zeta_{B}), and the color scale indicates the amount of entanglement of the corresponding resource state, CABC_{AB}. The confining curves comprise different families of channels: twin channels with ζA=ζB=ζ\zeta_{A}=\zeta_{B}=\zeta along the solid black curve; ac/gc (or gc/ac) channels along the dashed blue curve; dc/gc (or gc/dc) channels along the dashed-dotted red curve. Black circles correspond therefore to ζA=ζB=0\zeta_{A}=\zeta_{B}=0 (parallel ac), red stars to ζA=ζB=π/2\zeta_{A}=\zeta_{B}=\pi/2 (parallel dc), and blue triangles to ζA(B)=0,ζB(A)=π/2\zeta_{A(B)}=0,\zeta_{B(A)}=\pi/2 (hybrid ac/dc, dc/ac).

In order to make a complete analysis of the relation between τ4\tau_{4} and 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max}, a key element should be considered, namely the entanglement of the resource state. The evolution of CABC_{AB} and 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} for a given initial state and fixed channel parameters is qualitatively similar to that shown in Fig. 5, i.e., the quantum channel induces a loss in CABC_{AB}, which is verified by noticing that Eq. (55) is a decreasing function of pp. Yet, as occurred in the tripartite case, there exists a positive correlation between 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} and the multipartite entanglement (here measured by τ4\tau_{4}), that holds at each stage of the evolution, and for a fixed amount of CABC_{AB}. In order to verify this, it becomes crucial to recall first an important difference between the 3- and the 4-partite scenarios: as mentioned below Eq. (47), in the former case CABC_{AB} is independent of ζ\zeta (hence is the same for all channels), whereas in the 4-qubit example CABC_{AB} typically varies with both channel parameters, ζA\zeta_{A} and ζB\zeta_{B}. This observation, together with the essential role played by the entanglement of the resource state in the teleportation success, calls for an analysis that incorporates the value of CABC_{AB} when examining the dynamic interplay between τ4\tau_{4} and 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max}. To this end, we compare the effect of various channels (scanned by varying ζA\zeta_{A} and ζB\zeta_{B}) on the quantities of interest, namely τ4,𝔽maxgc/gc\tau_{4},\mathbb{F}^{\textsc{gc/gc}}_{\max} and CABC_{AB}, at the same stage of the evolution for different initial states. We do so by constructing triads (τ4,𝔽maxgc/gc,CAB)(\tau_{4},\mathbb{F}^{\textsc{gc/gc}}_{\max},C_{AB}), which for a given pp and fixed initial parameters φ\varphi and ϕ\phi, depend only on ζA,ζB\zeta_{A},\zeta_{B}. We then consider 16,384 different channels (pairs {ζA,ζB}\{\zeta_{A},\zeta_{B}\}), and display the resulting triads as points in the plane (τ4,𝔽maxgc/gc)(\tau_{4},\mathbb{F}^{\textsc{gc/gc}}_{\max}) colored according to the corresponding range of values of CABC_{AB}. Figure 9 shows the ensuing triads for φ=0\varphi=0 and different values of ϕ\phi (rows) and pp (columns). In all cases, τ4\tau_{4} and 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} were directly computed from Eqs. (56) and (55), respectively, whereas CABC_{AB} was numerically calculated as explained below Eq. (47).

All points in each panel of Fig. 9 lie within a region delimited by three curves:

  • The black solid line, encompassing the family of twin channels with ζA=ζB=ζ\zeta_{A}=\zeta_{B}=\zeta.

  • The blue dashed line, including those cases where one of the qubits (either AA or BB) undergoes an ac, meaning ζA/B=0\zeta_{A/B}=0, while the other one undergoes an arbitrary gc.

  • The red dotted-dashed line, comprising the channels where one of the qubits (either AA or BB) undergoes a dc, so ζA/B=π/2\zeta_{A/B}=\pi/2, while the other one is subject to an arbitrary gc.

Accordingly, the vertices of the regions are identified as:

  • Blue triangles represent the ac/dc (dc/ac) channels, with τ4=0\tau_{4}=0 and maximal average fidelity given by (65).

  • Red stars correspond to dc/dc, where 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} reaches its maximum value, in line with Eq. (VI.2).

  • Black circles represent the channel ac/ac, and correspond to points with the lowest 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} along the line of twin channels, in agreement with the statement below Eq. (64).

In all the graphs in Fig. 9, the colored bands, containing points whose CABC_{AB} lies within a specific range of values, reveal a positive correlation between the 4-tangle and the maximal average fidelity. Such correlation persists as the width of the range is reduced, as exemplified in Table 2 showing values of CABC_{AB}, τ4\tau_{4} and 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} taken from the central panel in Fig. 9. The values of τ4\tau_{4} and 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} are displayed in increasing order, and both increase simultaneously as the value of CABC_{AB} (numerically obtained) remains constant (up to 4 digits). The results indicate that despite the specificities present in each panel of Fig. 9, a feature common to all of them is that for a fixed entanglement of the resource state CABC_{AB}, 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} increases as τ4\tau_{4} increases.

ζA/π\zeta_{A}/\pi ζB/π\zeta_{B}/\pi CABC_{AB} τ4\tau_{4} 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max}
181/500 37/500 0.422003 0.0954376 0.760765
91/250 2/25 0.422002 0.0984991 0.761839
93/250 29/250 0.422003 0.120289 0.768959
187/500 143/1000 0.422001 0.139887 0.774978
187/500 19/125 0.422008 0.146878 0.777085
369/1000 99/500 0.422009 0.18453 0.788234
42/125 281/1000 0.422005 0.24261 0.805185
167/500 71/250 0.422008 0.243882 0.805556
Table 2: Values of CABC_{AB}, τ4\tau_{4} and 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} extracted from the panel e) of Fig. 9 (corresponding to ϕ=π/4\phi=\pi/4, φ=0\varphi=0, and p=0.5p=0.5), for different values of the channels parameters ζA\zeta_{A} and ζB\zeta_{B}.

In other words, for a given initial state (ϕ\phi fixed), and at each stage of the evolution (pp fixed), from among all the channels gc/gc that correspond to the same entanglement of the resource state, those that generate higher amounts of 4-partite entanglement lead to higher fidelities. This conclusion is analogous to that reached in the 3-partite case, indicating that the multipartite entanglement may act as a resource that protects the teleportation fidelity against the noisy channel (provided the same amount of CABC_{AB} is at disposal).

For the values of ϕ\phi and pp considered in panels a), b), c), f), and i) the condition (62) is satisfied, and accordingly the maximum of 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} is reached along with the maximum of the 4-tangle, given in this case by τ4max=02p2\tau^{\max}_{4}=\mathcal{E}^{2}_{0}p^{2}. In contrast, the cases depicted in panels d), g) and h) do not comply with the inequality (62), thus correspond to τ4max=16𝒫12p2(1p)2\tau^{\max}_{4}=16\mathcal{P}^{2}_{1}p^{2}(1-p)^{2} and 𝔽maxgc/gc=𝔽maxac/ac\mathbb{F}^{\textsc{gc/gc}}_{\max}=\mathbb{F}^{\textsc{ac/ac}}_{\max}, that is, to the minimal fidelity along the black solid curve. Panel e) corresponds to the case in which 0=4𝒫1(1p)\mathcal{E}_{0}=4\mathcal{P}_{1}(1-p), so as follows from Eq. (60) the maximum value of τ4\tau_{4} is attained simultaneously at the red star (dc/dc) and the black circle (ac/ac).

In all panels of Fig. 9 it is observed that as the black solid curve is traversed from the black circle to the red star, CABC_{AB} increases along with 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max}, so the highest fidelity is reached along with the highest CABC_{AB}. This shows that, within the family of twin channels, the resource state’s entanglement helps to improve the maximal average fidelity. However, this does not hold for other channels. For example, as the dashed blue line is traversed from the black circle to the blue triangle —that is, along the family of channels ac/gc—, there are cases in which an increase in CABC_{AB} is accompanied by a decrease in 𝔽maxac/gc\mathbb{F}^{\textsc{ac/gc}}_{\max}, as the one depicted in panel c). In such case it is only the 4-tangle what is seemingly enhancing the maximal average fidelity. Further, for the family of channels dc/gc, as the red dotted-dashed line is traversed, CABC_{AB} is kept constant (for fixed pp) and maximal —in fact, in this case we get CABdc/gc=0(1p)C^{\textsc{dc/gc}}_{AB}=\mathcal{E}_{0}(1-p)—, and the improvement of 𝔽maxdc/gc\mathbb{F}^{\textsc{dc/gc}}_{\max} is therefore ascribable to the increase of τ4\tau_{4}.

VII Closing remarks

We investigated the role of multipartite entanglement in the dynamics of the maximal average fidelity (above the classical threshold value) when the teleportation protocol includes the action of a local quantum channel, ΛAB=ΛAΛB\Lambda_{AB}=\Lambda_{A}\otimes\Lambda_{B}, acting on the qubits AA and BB that conform the resource state. To facilitate our goal we expressed the maximal average fidelity in terms of the Kraus operators corresponding to a general 2-qubit channel (Eq. (8)), and introduced the Kraus operators of a generalized noisy channel, Eq. (27), which encompasses the paradigmatic amplitude damping and dephasing channels, and connects them via a continuous parameter that also determines the amount of multipartite entanglement created along the evolution.

We first considered the case ΛA=𝖨2\Lambda_{A}=\mathsf{I}_{2}, and ΛB\Lambda_{B} representing the generalized noisy channel, rooted at the interaction of BB with an additional qubit EBE_{B}. 3-partite entanglement thus typically emerges among the qubits A,BA,B, and EBE_{B}, here quantified by the 3-tangle τABEB\tau_{ABE_{B}}. In the second scenario, both ΛA\Lambda_{A} and ΛB\Lambda_{B} correspond to generalized channels. Interpreting them as the effective result of a local interaction of AA and BB with additional qubits EAE_{A} and EBE_{B}, we focused on the ensuing 4-partite entanglement among the parties A,B,EAA,B,E_{A}, and EBE_{B}, as measured by τ4\tau_{4}.

In both cases, we found that the relative phase φ\varphi that optimizes both 𝔽max\mathbb{F}_{\max} and τ4\tau_{4} is φ=0\varphi=0. More importantly, the analytical and numerical analysis (considering identically parametrized channels in the 4-party case), revealed that for a fixed non-zero amount of the resource state’s entanglement (i.e., for fixed CAB0C_{AB}\neq 0), and at each stage of the evolution, the teleportation success improves under channels that induce higher amounts of multipartite entanglement. Here it is important to stress that both the 3-tangle and the 4-tangle quantify a specific type of multipartite entanglement —namely the 3-way and the 4-way entanglement, respectively— characteristic of ghz-type states, and thereby absent in w-type states [44, 51]. Consequently, our findings indicate that it is specifically the amount of ghz-type entanglement which favors better teleportation fidelities. It should be stressed that this conclusion does not go against the (expected) fact that, given a specific noisy channel (a fixed ζ\zeta), as the evolution takes place the teleportation fidelity decays while the 3- and 4-tangle may increase; rather, the conclusion compares the effect of different generalized noisy channels on 𝔽max\mathbb{F}_{\max} throughout the evolution, and establishes that under channels that produce more multipartite (ghz-type) entanglement, the detrimental effects on 𝔽max\mathbb{F}_{\max} are lessened. Further, when CABC_{AB} vanishes, as occurs at p=1p=1 —when the channel has suppressed all the entanglement of the resource state—, the maximal average fidelity drops below the classical threshold value, despite the multipartite entanglement may acquire relatively large values. This highlights CABC_{AB} as a necessary element that triggers the power of the nn-way entanglement to enhance the teleportation success.

In the 3-party case, CABC_{AB} is the same for all the channels considered, and the relation between 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max} and τABEB\tau_{ABE_{B}} was clearly revealed, along with the identification of the dc as the channel that produces the higher values of 𝔽maxgc\mathbb{F}^{\textsc{gc}}_{\max}. In the 4-qubit case, in contrast, CABC_{AB} depends on the specific channel, and the relation between 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} and τ4\tau_{4} only (i.e., without considering the value of CABC_{AB}) is more subtle than in the 3-party counterpart. In particular, we found that an increment in τ4\tau_{4} induces an increment in 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max}, irrespective of the initial state and stage of the evolution, only when the composite channel is dc/gc, that is, when either one of the qubits is subject to a dephasing channel. Notably, in this case, CABC_{AB} does not depend on the specific gc, and the enhancement of the maximal average fidelity is due solely to the increase in the 4-partite entanglement.

For the families of twin channels gc/gc with ζA=ζB=ζ\zeta_{A}=\zeta_{B}=\zeta, and ac/gc, a higher value of τ4\tau_{4} is not always accompanied by a higher value of 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max}, ultimately because CABC_{AB} changes with the channels’ parameters. Instead, for a fixed initial state and at a given stage of the evolution, the composite channel for which the 4-tangle is maximal is either dc/dc —in which case 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} attains its global maximal value—, or ac/ac, corresponding to a 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max} that is minimal within the family of twin channels (ζA=ζB=ζ\zeta_{A}=\zeta_{B}=\zeta).

Interestingly, guaranteeing that 𝔽max\mathbb{F}_{\max} exceeds the classically attainable value depends not only on the initial entanglement at disposal, but also on the initial population 𝒫1\mathcal{P}_{1} of the state |11AB\left|{11}\right\rangle_{AB}, which plays against the improvement of the teleportation success. Further, for sufficiently low values of 𝒫1\mathcal{P}_{1} (𝒫10.2\mathcal{P}_{1}\leq 0.2) the maximal value of τ4\tau_{4} is attained together with the maximal value of 𝔽maxgc/gc\mathbb{F}^{\textsc{gc/gc}}_{\max}.

Our analysis led us to conclude that the previously reported improvement of the teleportation fidelity under some types of noisy channels [11, 12, 13, 14, 15, 16] may be rooted at the emergence of multipartite entanglement induced by the interaction of the resource qubits with their surroundings. This offers valuable insights into the power of multipartite correlations, as well as into the characterization of the processes that better protect the teleportation fidelity in the more realistic scenario in which A+BA+B is an open system. In particular, processes that generate ghz-type states in the 3- or 4-qubit system have the potential to assist the protocol by reducing the detrimental effects of noise, as a result of the induced generation of 3- and 4-way entanglement. This highlights the ghz-type entanglement as a useful auxiliary resource in noisy quantum teleportation.

Acknowledgements.
The authors acknowledge financial support from DGAPA, UNAM through project PAPIIT IN112723. VHTB acknowledges CONAHCYT scholarship with CVU: 863195.

References