## High-dimensional quantum teleportation under noisy environments

Alejandro Fonseca

Departamento de F´ ısica, CCEN, Universidade Federal de Pernambuco, 50670-901, Recife, PE, Brazil (Dated: December 6, 2019)

We study the protocol of qudit teleportation using quantum systems subjected to several kinds of noise for arbitrary dimensionality d . We consider four classes of noise: dit-flip, d -phase-flip, ditphase-flip and depolarizing, each of them corresponding to a family of Weyl operators, introduced via Kraus formalism. We derive a general expression for the average fidelity of teleportation in arbitrary dimension d for any combination of noise on the involved qudits. Under a different approach we derive the average fidelity of teleportation for a more general scenario involving the d -dimensional generalization of amplitude damping noise as well. We show that all possible scenarios may be classified in four different behaviours and discuss the cases in which it is possible to improve the fidelity by increasing the associated noise fractions. All our results are in agreement with previous analysis by Fortes and Rigolin for the case of qubits (Phys. Rev. A, 92 012338, 2015).

## I. INTRODUCTION

Since its proposal in 1993 [1] until nowadays, the teleportation protocol represents one of the most known and widely studied applications of quantum entanglement [2]. According to it, if one has a source of maximally entangled qudits and a measurement apparatus capable of discriminating the d 2 elements of the generalized Bell basis, then it is possible to send an arbitrary qudit state between two locations even without prior knowledge of it. Nevertheless, in day by day experiments some additional features arise from the unavoidable interaction of the involved parts with the environment and/or imperfections in the preparation of the system, leading to losses of the resources responsible for the improvement in the execution of the task over its classical analogue. Thus, in addition to adopting strategies aiming to diminish the action of noise, one may also modify the scheme of measurements and operations on the parts and consequently to improve the performance of the protocol [3, 4]. For this reason, in order to effectively optimize the strategies to be implemented it is important to have a general picture of the scenario. In this respect, the problem of characterizing the protocol of teleportation in the presence of noisy environments have been addressed from several perspectives: In [5], Oh and collaborators employ the Lindblad operators formalism, obtaining the fidelity of teleportation for several classes of noise in the quantum channel. More recently, Fortes and Rigolin have presented a set of results within the frame of Kraus operators for some of the most known instances of noise in the literature [6]. By using the same approach, the authors of [7] contrast theoretical predictions with experimental results. In addition, there are even some approaches to the problem of multipartite noisy teleportation [8-10].

Apart from a few exceptions, most of the effort in describing quantum information protocols under noisy scenarios has been focused on systems involving qubits, nevertheless it has been shown that the performance of several tasks is enhanced when high dimensional systems are used instead [11-15]. In fact, two independent groups have been able to successfully implement the teleporta- tion protocol of qutrits recently [16, 17]. Furthermore, because of the limited efficiency in usual detectors, it is more convenient the usage of single qudit systems in the process of generation and manipulation of high dimensional entanglement, rather than multi-qubit based structures [18]. In this way, the exploration of quantum information protocols employing qudits under realistic conditions becomes an important aspect of quantum information science.

In this paper we present a characterization of the standard protocol of qudit teleportation considering several sources of imperfection, including non-maximally entangled channel and/or joint measurements and qudits susceptible to modifications of their states due to the interaction with their environment. The paper is organized as follows: first we present the protocol, derive general expressions for the fidelity of teleportation and study the noiseless case. In section III we introduce the Kraus operators used throughout this work to model the effect of noise in the system. Then we consider the simplest case in which noise acts on a single qudit and extend results to more general scenarios involving more than one part affected. The last part is devoted to discuss our main results and present some conclusions.

## II. TELEPORTATION PROTOCOL

The standard teleportation protocol involves two parts: Alice and Bob, as usual (see Fig. 1), sharing a pair of entangled qudits whose state may be described by a density operator ˆ ρ ch . Alice has an additional qudit prepared in an arbitrary, not necessarily known qudit state | φ 〉 = ∑ d -1 j =0 α j | j 〉 , her task is to send it to Bob. For this, she carries out a d 2 -outcome joint projective measurement on her pair of qudits in a generalized Bell-like basis {| Φ mn 〉} , with elements given by

$$\left | \Phi _ { m n } \right \rangle = \sum _ { k = 0 } ^ { d - 1 } \beta _ { k m } \left | k, k \oplus n \right \rangle, \text{ \quad \ \ } ( 1 )$$

FIG. 1. Teleportation scheme: Alice and Bob share a channel composed by two entangled qudits prepared in a state ˆ ρ ch (red dashed). Alice performs a joint measurement on qudits I and A (blue dotted). Under appropriate conditions (maximal entanglement) and after LOCC (green dot dashed), corresponding to transmission of a pair of dits ( m,n ) using a classical channel and subsequent application of a local unitary operation ˆ U mn on qudit B , it finally holds in the initial state of the input qudit I .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01097v4-figure-ab196d0d-8e7e-49fa-bbb1-3c6c3199d6da.jpg)

where the symbol ' ⊕ ' denotes sum modulo d , the β km coefficients account for the extent of entanglement and satisfy the relation ∑ d -1 k =0 β km β ∗ km ′ = δ mm ′ . In this way, the usual maximally entangled joint measurements are recovered whenever β km = ω k m · d / √ d , where ω d = exp(2 πi/d ) is the primitive d -th root of unity.

By employing a classical channel, Alice sends the information about her measurement outcome ( m,n ) to Bob who applies a local unitary operation on his qudit, given by one out of the d 2 Weyl operators ˆ U mn , defined as [19]

$$\hat { U } _ { m n } = \sum _ { j = 0 } ^ { d - 1 } \omega _ { d } ^ { j m } \left | j \right \rangle \langle j \oplus n \right |. \quad \ \ ( 2 ) \quad \ \max _ { ( \gamma _ { k } ) } ^ { \ddots }$$

After each run, the state of Bob's qudit (up to normalization) holds

$$\hat { \rho } _ { m n } = \hat { U } _ { m n } \ t r _ { A } \left \{ \left ( | \Phi _ { m n } \rangle \langle \Phi _ { m n } | \otimes \hat { 1 } _ { B } \right ) | \phi \rangle \langle \phi | \otimes \hat { \rho } _ { c h } \right \} \hat { U } _ { m n } ^ { \dagger }, \ \text{                                                                                                                                                                                                        
                                                                                                                                                                                                        }$$

where tr A denotes the partial trace on the subsystem associated to Alice.

The reliability of the protocol is usually assesed by calculating the fidelity of teleportation i.e. how close the state in Bob is after the process to that initially possessed by Alice: F mn = tr (ˆ ρ mn | φ 〉 〈 φ | ) tr (ˆ ρ mn ) -1 , where the normalization factor tr (ˆ ρ mn ) -1 is equal to the probability of occurrence of the ( m,n ) output. By considering the whole possible measurement results, the mean fidelity reads

$$F = \sum _ { m n } \text{tr} \left \{ | \phi \rangle \langle \phi | \, \hat { \rho } _ { m n } \right \}. \quad \quad ( 3 ) \quad \text{ent}$$

Given that our goal is to assess the quality of the protocol, independent of the teleported state and due to the fact that F typically depends on the coefficients α j , then it is more convenient to calculate the average fidelity over the set of input states

$$\langle F \rangle = \frac { 1 } { V _ { d } } \int \mathrm d \Gamma _ { d } \ F.$$

General expressions for dΓ d and V d are given in appendix A. For the sake of simplicity in what follows we will refer to 〈 F 〉 as the fidelity instead of average fidelity of teleportation.

## A. Noise-free environment

Let us assume the qudits composing the channel are initially prepared in a pure partially entangled state, i.e. ˆ ρ ch = | ψ ψ 〉 〈 | , with | ψ 〉 = ∑ d -1 k =0 γ k | kk 〉 . After some calculations the fidelity of teleportation is reduced to

$$\int _ { \substack { \stackrel { m \lambda } { l } \stackrel { \lambda } { \text{state} } } } ^ { \stackrel { m \lambda } { \text{state} } } \left \langle F \right \rangle = f _ { C } \left \{ 1 + \frac { 1 } { d } \sum _ { \substack { m, n = 0 \\ j > k = 0 } } ^ { d - 1 } \text{Re} \left [ \omega _ { d } ^ { m ( k - j ) } \beta _ { j m } \beta _ { k m } ^ { * } \gamma _ { k \oplus n } \gamma _ { j \oplus n } ^ { * } \right ] \right \},$$

where f C = 2 d +1 is known as classical fidelity and corresponds to the maximal value attained by the fidelity of teleportation after the usage of any classical strategy [20, 21], e.g. by using a classical channel, Alice sends to Bob the result of a projective measurement carried out on the input qudit. Then, based on this information he guesses the state his qudit has to be prepared.

Remarkably, in the expression above the classical and quantum contributions to the fidelity are made clear. Furthermore, it is easy to see that whenever we have maximal entanglement in the measurements and channel ( γ k = 1 / √ d ), the second term is reduced to ( d -1) / d ( +1), and in this way the fidelity reaches its maximum value, 〈 F 〉 = 1, as expected.

In order to have a qualitative picture of how the amount of entanglement in the channel is related to the quantum contribution to the fidelity of teleportation f Q = 2 d d ( +1) ∑ d -1 m,n =0 j&gt;k =0 Re ( ω m k ( -j ) d β jm β ∗ km γ k ⊕ n γ ∗ j ⊕ n ) , we produced a sample of random entangled states uniformly distributed in the space of the Schmidt basis [22]. Results of the quantum contribution to the fidelity of teleportation (normalised to one), f ′ Q [= ( d +1) f Q / d ( -1)], as a function of the amount of entanglement in the channel for the case in which the measurements basis is maximally entangled are presented in figure 2. Note that d = 2 is the only case in which there is a one to one relation between f ′ Q and the amount entanglement. In contrast, it is only possible to infer bounds in the extent of channel's entanglement for a given value of the fidelity within the high dimensional case. It is not difficult to show that for

FIG. 2. Normalized quantum contribution to the fidelity f ′ Q [= ( d +1) f Q / d ( -1)] in function of the amount of entanglement in the channel (normalized to log d ), for a maximally entangled measurement basis and a set of N random pure states of two entangled qudits, for several values of d . (a) d = 2, N = 10 . 4 (b) d = 3, N = 10 . 5 (c) d = 4, N = 10 . 5 (d) d = 5, N = 10 . 6 Each line corresponds to a particular family of 'boundary' states ϕ µ , given by Eq. 5: red dotted : µ = 1, green dashed : µ = 2, blue solid : µ = 3 and cyan dash-dotted : µ = 4. Note also that any intersection point between two family lines corresponds to a maximally symmetric ν -rank state φ ν (Eq. 6).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01097v4-figure-ee2529ca-5bd2-4283-ad0d-b89b609a2989.jpg)

each dimensionality d , such a boundary is fully described by d -1 families of states ϕ µ :

$$| \varphi _ { \mu } \rangle = a _ { \mu } \, | 0 0 \rangle + \sqrt { \frac { 1 - a _ { \mu } ^ { 2 } } { \mu } } \sum _ { k = 1 } ^ { \mu } | k k \rangle \,, \quad \ ( 5 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where µ = 1 , · · · , d -1, with a µ ∈ [ 0 1 , / √ µ +1 ] for µ = 1 , · · · , d -2 and a µ ∈ [0 , 1] for µ = d -1. It is possible to see that maximally symmetric ν -rank states,

$$| \phi _ { \nu } \rangle = \frac { 1 } { \sqrt { \nu } } \sum _ { k = 0 } ^ { \nu - 1 } | k k \rangle \,, \quad \quad ( 6 ) \, \begin{array} { c } \Delta \sigma \\ \text{ell} \sigma \\ \text{wh} \\ \text{isf} \end{array}.$$

with ν = 1 , · · · , d -1, are extremal cases of the families ϕ µ .

It is worth to remark that the relation between fidelity of teleportation and bounds of entanglement presented above holds whenever the source is capable of providing entangled pairs in pure states. Nevertheless, the previous analysis does not hold when the channel is prepared in a mixed state, given that as shown in the next section, the classical contribution is also modified in this case. In this respect, by employing the recently introduced concept of nonclassical teleportation witnesses [23, 24], ˇ upi´, S c Skrzypczyk and Cavalcanti have found that it is possible to infer lower bounds of channel entanglement from teleportation data only [25], even for scenarios in which Alice and Bob share pairs prepared in some special families of mixed entangled states not capable of attaining fidelity of teleportation values above its classical limit, also known as bound entangled states [26].

FIG. 3. Weyl operators ˆ U mn and their relation with Kraus operators for several kinds of noise on d -dimensional systems. The blue row represents dit-flip like operators, the yellow column d -phase-flip like operators and the pink squares are related to matrices corresponding to dit-phase-flip like noise. Note that the three classes mentioned before are employed to define depolarizing noise (blue dashed line). In addition, the set of Pauli matrices corresponds to the nontrivial operators for d = 2 (red dashed line).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01097v4-figure-b88e3395-729c-4fb8-a522-b985e7c21f79.jpg)

## III. NOISE AND KRAUS OPERATORS

In addition to the technical limitations in the preparation of the system and realization of measurements in maximally entangled states, noise is an unavoidable feature of real experiments, for this reason it is very important to establish strategies which lead to the improvement of the final results. In the following sections we explore the influence of protecting one or more qudits from noise on the fidelity of teleportation.

To date there are several methods to study the evolution of the state ˆ, ρ associated to an open quantum system [27]. One of the most widely used is the Kraus operators formalism, in which the evolution may be modelled by a trace preserving map ˆ ρ → ˆ ρ ′ = ∑ k ˆ E ρE k ˆ ˆ † k , where the ˆ E k 's are known as Kraus operators and satisfy the completeness relation ∑ k ˆ E E † k ˆ k = ˆ 1 [27]. Under this approach the ˆ E k operators contain information about the effects of the system-environment interaction, without the necessity of deepening into the involved physical processes behind. Let us introduce the relation between high-dimensional generalizations of the most known classes of error in quantum systems and Weyl operators.

## A. Noise and Weyl operators

Kraus operators corresponding to some of the most known instances of quantum noise, namely: bit-flip, phase-flip, bit-phase-flip and depolarizing for d = 2 and d = 3 are presented in [27] and [28] respectively. In addition, it is not difficult to see that arbitrary d -dimensional generalizations are proportional to families of Weyl oper-

ators ˆ U mn (Eq. 2), such a correspondence is illustrated in figure 3. This class of operators represent a very important tool when dealing with quantum information tasks based on high dimensional systems. Among them it is possible to highlight the quantum teleportation protocol [1, 21], representation of quantum states [19], random unitary evolutions [29], quantum computing [30, 31] and quantum error correction [32]. In the latter case, the elements of set the of Weyl operators (also known as single-qudit Pauli group ) are often found in literature as products of the form ω X k d ˆ m ˆ Z n , where ˆ X and ˆ Z correspond to the operators ˆ U 01 and ˆ U 10 , respectively and 0 ≤ k, m, n ≤ d -1 [31-33]. In addition, note that the set { ˆ U mn } constitute a natural basis for the d × d HilbertSchmidt space and that Weyl operators are proportional to the set of Pauli matrices for d = 2 (see Fig. 3).

Let us briefly describe of each particular noise and the corresponding expressions for the Kraus operators used in this work:

- · dit-flip noise: In analogy to bit-flip for d = 2, this kind of noise considers perturbations that flip | j 〉 either to the state | j ⊕ 〉 1 , | j ⊕ 〉 2 , . . . , or | j ⊕ d - 〉 1 , with probability p . The associated Kraus operators are: ˆ E 00 = √ 1 -p U ˆ 00 , ˆ E 01 = √ p d -1 ˆ U 01 , . . . , ˆ E 0 ,d -1 = √ p d -1 ˆ U 0 ,d -1 .
- · d-phase-flip noise: A qudit | j 〉 subjected to d -phase-flip noise may with probability p suffer one out of d -1 phase shifts of the form: ω d | j 〉 , ω 2 d | j 〉 , . . . , ω d -1 d | j 〉 . The corresponding Kraus operators are given by: ˆ E 00 = √ 1 -p U ˆ 00 , ˆ E 10 = √ p d -1 ˆ U 10 , . . . , ˆ E d -1 0 , = √ p d -1 ˆ U d -1 0 , .
- · dit-phase-flip noise: This is a special case in which a combination of both former kinds of noise may take place, e.g. a qudit suffer a flip and a phase shift at the same time. The related Kraus operators are: ˆ E 00 = √ 1 -p U ˆ 00 and ˆ E mn = √ p d -1 ˆ U mn , with 1 ≤ m,n ≤ d -1.

glyph[negationslash]

- · Depolarizing noise: Under this, a system initially prepared in an arbitrary state evolves to a maximally mixed state, ˆ 1 /d with probability p . The Kraus operators for this are given by [30, 31]: ˆ E 00 = √ 1 -d 2 -1 d 2 p U ˆ 00 and ˆ E mn = √ p d ˆ U mn , with 0 ≤ m,n ≤ d -1, for ( m,n ) = (0 , 0).

It is important to mention that the present work does not refer to any particular implementation, for this reason the noise models we have just described consider that possible errors can occur with equal probability. For instance, a qutrit under interactions with the environment leading to a dit-flip noise may jump either one or two levels above, each with probability p/ 2.

The Kraus operators for the classes of noise mentioned above may be written as: ˆ E mn = a mn ˆ U mn , with coefficients a mn ∈ R , satisfying ∑ mn a 2 mn = 1. Given an arbitrary system initially prepared in a state ˆ = ρ ∑ glyph[vector]glyph[vector] kl ρ glyph[vector]glyph[vector] kl | glyph[vector] k 〉〈 | glyph[vector] l , where N is the number of subsystems, glyph[vector] k = ( k , ..., k 1 N ) and 0 ≤ k j ≤ d -1, the action of a set of Kraus operators ˆ E glyph[vector]glyph[vector] kl = ˆ E k l 1 1 ⊗ ... ⊗ ˆ E k N N l = ∏ N j =1 a k l j j ˆ U k l 1 1 ⊗ ... ⊗ ˆ U k N N l , transforms ˆ into ρ ˆ ρ ′ = ∑ glyph[vector] mglyph[vector] n ρ ′ glyph[vector] mglyph[vector] n | glyph[vector] m glyph[vector] 〉 〈 n | , with ρ ′ glyph[vector] mglyph[vector] n given by

$$\rho ^ { \prime } _ { \vec { m } \vec { n } } = \sum _ { \vec { k } \vec { l } } \omega _ { d } ^ { \vec { k } \cdot ( \vec { m } - \vec { n } ) } \rho _ { \vec { m } \oplus \vec { l }, \vec { n } \oplus \vec { l } } \prod _ { j = 1 } ^ { N } a _ { k _ { j } l _ { j } } ^ { 2 }. \quad ( 7 )$$

This expression contains the information about noise acting on the whole system and it is very useful in the calculation of a general expression for the fidelity of teleportation (see App. C).

An important class of noise related to losses in quantum systems is known as amplitude damping noise . However, in contrast to the previous cases, the related Kraus operators are not proportional to single Weyl operators. For this reason, calculations of fidelity involving amplitude damping noise were carried out using the standard computational basis. Now let us briefly present its corresponding Kraus operators and the general expression for modified density operator coefficients.

## B. Amplitude damping noise

Amplitude damping noise has been used to model a large amount of phenomena including idle errors in quantum computing [31], energy dissipation, spontaneous photon emission, attenuation, among others in two level systems [27]. A d -dimensional generalization was recently introduced in [34], the corresponding Kraus operators read

$$\hat { E } _ { 0 } = | 0 \rangle \langle 0 | + \sqrt { 1 - p } \sum _ { j = 1 } ^ { d - 1 } | j \rangle \langle j | \,, \quad \quad ( 8 )$$

and

$$\hat { E } _ { j } = \sqrt { p } \left | 0 \right \rangle \langle j \right |,$$

with j = 1 , · · · , d -1. This kind of noise may be interpreted in the following way: A d -level system interacting with its environment may with probability p lose population from the excited levels, leading the system to the ground state | 0 . 〉

In general terms, for a system composed by one part only, we can write any Kraus operator as: ˆ E k = ∑ mn a ( k ) mn | m n 〉 〈 | , with coefficients a ( k ) nm ∈ C , satisfying ∑ kn a ( k ) mn a ( k ) ∗ ln = δ ml (due to the completeness relation). The N -party case is a straightforward generalization ˆ ρ → ˆ ρ ′ = ∑ glyph[vector] k ˆ E glyph[vector] k ˆ ρ E ˆ † glyph[vector] k , where glyph[vector] k = { k , . . . , k 1 N } , as previously and Kraus operators given by ˆ E glyph[vector] k = ˆ E ( a 1 ) k 1 ⊗ ˆ E ( a 2 ) k 2 · · · ⊗ ˆ E ( a N ) k N = ∑ glyph[vector] mglyph[vector] n ( ∏ N j =1 a ( k j ) m n j j ) | glyph[vector] m glyph[vector] 〉 〈 n | . In this

case the modified density operator coefficients ρ ′ glyph[vector]glyph[vector] pq are reduced to

$$\rho _ { \vec { p } \vec { q } } ^ { \prime } = \sum _ { \vec { k } \vec { m } \vec { n } } \left ( \prod _ { j = 1 } ^ { N } a _ { p _ { j } m _ { j } } ^ { ( k _ { j } ) } a _ { q _ { j } n _ { j } } ^ { ( k _ { j } ) * } \right ) \rho _ { \vec { m } \vec { n } }. \quad ( 1 0 )$$

The coefficients a ( k j ) mn may be easily calculated for each particular case as usual: a ( k j ) mn = 〈 m E | ˆ k j | n 〉 .

We have performed calculations of fidelity of teleportation for the cases in which the three qudits involved in the process may be affected by arbitrary combinations of the classes of noise described previously. For details, we refer the reader to Appendix C. The following sections are devoted to present results concerning some particular cases.

## IV. NOISE ACTING ON A SINGLE QUDIT

## A. Weyl-like noises

In this part we consider the case of two qudits fully protected from noise, e.g. an experiment in which the production of pairs of entangled qudits is carried out in Alice's location and the Bob's qudit is affected by interacting with the environment during the transportation process.

By direct substitution into the expression for fidelity (Eq. C3), it is easy to see that when noise is acting on one qudit only, the fidelity does not depend on the qudit affected. Then, giving continuity to the example given above, we consider that the affected qudit is that on Bob's location. In this case the general expression for the fidelity of teleportation is reduced to:

$$\langle F \rangle = \frac { 1 } { d + 1 } \left [ 1 + d c _ { p } ^ { 2 } + c _ { 0 } ^ { 2 } - c _ { p } ^ { 2 } + ( d + 1 ) ( c _ { 0 } ^ { 2 } - c _ { p } ^ { 2 } ) f _ { Q } \right ], \quad \text{nois}$$

which does not depend on the coefficient c f and for this reason the fidelity corresponding to dit-flip 〈 F F 〉 and ditphase-flip 〈 F FP 〉 noises are both equal to:

$$\langle F _ { F } \rangle = \langle F _ { F P } \rangle = \frac { 2 } { d + 1 } \left ( 1 - \frac { p } { 2 } \right ) + f _ { Q } ( 1 - p ). \quad ( 1 1 )$$

For d -phase-flip noise, the fidelity is reduced to:

$$\langle F _ { P } \rangle = \frac { 2 } { d + 1 } + f _ { Q } \left ( 1 - \frac { d } { d - 1 } p \right ), \quad ( 1 2 ) \quad \overset { \ddot { \text{d} } } { \text{By} } \, \iota _ { \text{b-t} }$$

the classical fidelity is not affected because phase shifts are exclusive elements of quantum systems. This feature will be explored in more detail in the following subsection.

The corresponding fidelity for depolarizing noise is:

$$\langle F _ { D } \rangle = \frac { 2 } { d + 1 } \left ( 1 - \frac { d - 1 } { 2 d } p \right ) + f _ { Q } \left ( 1 - p \right ). \quad \left ( 1 3 \right ) \quad \text{entar} \\ \beta _ { j m } \colon$$

Note that when we have a maximally entangled channel and measurements, the fidelities 〈 F F 〉 , 〈 F P 〉 , and 〈 F FP 〉 are all equal to:

$$\langle F \rangle = 1 - \frac { d } { d + 1 } p, \text{ \quad \ \ } ( 1 4 )$$

and the corresponding to depolarizing reduces to:

$$\langle F _ { D } \rangle = 1 - \frac { d - 1 } { d } p. \text{ \quad \ \ } ( 1 5 )$$

The noise thresholds [critical noise fractions above which the fidelity of teleportation acquire values below the classical limit, f C = 2 ( / d +1)] are given by p ∗ = ( d -1) /d and p ∗ D = d/ d ( +1).

## B. Amplitude-Damping noise

In this case it is not possible to write a closed expression in terms of the quantum contribution to the fidelity of teleportation f Q , in the same way as in the previous section, for arbitrary dimension, d . For a maximally entangled measurement basis and channel, the fidelity reads

$$\text{$\text{cut in}$} _ { \text{ten}$} _ { \text{:ation} } \quad \langle F _ { A D } \rangle = \frac { 2 } { d + 1 } \left [ \frac { d ^ { 2 } - d + 2 } { 2 d } - \frac { ( d - 1 ) ^ { 2 } } { 2 d } p + \frac { d - 1 } { d } \sqrt { 1 - p } \right ].$$

The noise threshold p ∗ AD in this case is given by

$$\Lambda _ { A } \nu \\ p _ { A D } ^ { * } = \frac { d + 2 \sqrt { d } } { \left ( \sqrt { d } + 1 \right ) ^ { 2 } }.$$

Note that the above expression for the fidelity of teleportation is valid whenever two parts are fully protected and only one qudit is being affected by amplitude damping noise, no matter which one. Nevertheless, this symmetry ceases to appear either when entanglement is not maximal and/or another qudit suffers the action of any kind of noise.

## C. Optimization of fidelity under d -phase-flip noise in one qudit

Besides the fact that the classical fidelity is not affected by the presence of d -phase-flip noise acting on a single qudit, it is possible to find some other interesting features. By analysing the expression for fidelity (Eq. 12), we see that above a noise threshold p ∗ = ( d -1) /d , the coefficient accompanying the quantum contribution becomes negative. This situation may be overcome if we make a phase addition in the measurement basis, as pointed out in [6] for d = 2. Without loss of generality and in order to simplify calculations, assume a channel initially prepared in a maximally entangled state and a maximally entangled measurement basis with arbitrary phases φ j : β jm = e iφ j / √ d , with φ 0 = 0. The fidelity then holds:

$$\langle F _ { P } \rangle = \frac { 2 } { d + 1 } \left \{ 1 + \frac { 1 } { d } \left ( 1 - \frac { p d } { d - 1 } \right ) \left ( \sum _ { k = 1 } ^ { d - 1 } \cos \phi _ { k } + \sum _ { k > l = 1 } ^ { d - 1 } \cos ( \phi _ { l } - \phi _ { k } ) \right ) \right \}.$$

FIG. 4. (solid): Optimal fidelity of teleportation for the case in which only one of the qudits may suffer d -phase-flip noise. (dashed): Classical fidelity. (a) Calculations of 〈 F p 〉 for 2 ≤ d ≤ 5. (b) Optimal fidelity of teleportation for arbitrary dimension d .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01097v4-figure-1b7ee81c-4b4a-49bc-a3fa-5079988dc7aa.jpg)

The problem is thus reduced to a optimization procedure in which we search for extremal values (maximum when p &lt; p ∗ and minimum for p &gt; p ∗ ) of the quantum contribution to the fidelity. We carried out analytical calculations up to d = 3, obtaining following results: For noise fractions below the treshold p ∗ , the whole set of phases are null, as expected. For p &gt; p ∗ , we got φ 1 = π for d = 2 and ( φ , φ 1 2 ) = (2 π/ , 3 4 π/ 3) for d = 3. The resulting fidelities are plotted in figure 4(a). Further numerical calculations were performed from which we were able to infer the following expressions for the optimized fidelity:

$$\langle F _ { P } \rangle = \begin{cases} 1 - \frac { d p } { d + 1 } & \text{for $p<p^{*}$} \\ \frac { d p + d - 1 } { d ^ { 2 } - 1 } & \text{for $p>p^{*}$.} \\ \quad & \text{to $t$} \\ \text{deli} \\ \quad & \text{theo} \end{cases}$$

In conclusion, if somehow either Alice or Bob are capable of estimating the associated d -phase-flip noise fraction on the affected qudit, then she can improve the fidelity of the teleported state by choosing one out of two measurement basis. The results are summarized in Fig. 4(b) for arbitrary d . As it can be seen, the best improvement is attained by systems composed by qubits. As Fortes and Rigolin have shown [6], this feature can be exploited if we permit a part of the system to be strongly affected by phase-flip noise. Nevertheless it is worth to remark that such a recovery in the fidelity of qubit teleportation reported in [6] may be explained by the fact that an increase in the noise fraction p leads to an effective suppression of the mixedness in the final state as p → 1. Unfortunately this is not the case for arbitrary dimension, given that there are more than one phases involved, thus such a recovery becomes lower as we increase d , as it can be seen in Fig. 4(a).

## V. NOISE IN MORE THAN ONE QUDIT

In this section we explore the case in which protection may be applied in at most one of the qudits [35]. In order to have a best insight from the results, we assume maximal entanglement in the channel and measurements.

Before examining several cases in detail let summarize some general results: When entanglement is maximal, given either one, two or three classes of Weyl-like noises acting on the system, the fidelity does not depend on how those are distributed on the qudits. In this way we have 〈 F X,Y,Z 〉 ( p X , p Y , p Z ) = 〈 F Y,X,Z 〉 ( p Y , p X , p Z ) = 〈 F Z,X,Y 〉 ( p Z , p X , p Y ) = . . . , where 〈 F X,Y,Z 〉 ( p X , p Y , p Z ) indicates the teleportation fidelity given that X Y , and Z noises are acting on the input, Alice's and Bob's part of the channel, with noise fractions p X , p Y and p Z respectively. For instance, any situation in which two qudits may be affected is equivalent to that of having the input protected only, i.e. 〈 F X, ∅ ,Y 〉 = 〈 F X,Y, ∅ 〉 = 〈 F ∅ ,X,Y 〉 , where the symbol ' ∅ ' stands for a noise-free qudit.

glyph[negationslash]

glyph[negationslash]

In addition, some interesting results arise when we consider dit-flip, d -phase-flip and dit-phase-flip noises: 〈 F F,F,X 〉 = 〈 F P,P,X 〉 with X = { ∅ , F P, D } , 〈 F X ,Y ′ ′ ,F 〉 = 〈 F X ,Y ′ ′ ,P 〉 = 〈 F X ,Y ′ ′ ,FP 〉 for X ′ = Y ′ = { ∅ , D } and 〈 F X,Y, ∅ 〉 = 〈 F X,Z, ∅ 〉 for X = Y = Z = { F, P, F P } explicit expressions are not presented here, however all may be obtained by direct substitution in equation C4.

glyph[negationslash]

With the exception of some specific cases listed below, we observed that whenever the system is subjected to the instances of noise considered in this work, the fidelity of teleportation exhibits quite the same behaviour: the larger the amount of noise, the lower the values attained by 〈 F 〉 . As an illustration, we have plotted results of teleportation fidelity in function of the corresponding noise fractions for the scenario ( FP,F,P ) and d = 3 in Figure 5. Furthermore it was observed that as the dimension increases, the set of noise configurations leading to a fidelity of teleportation above the classical limit becomes larger. Nevertheless this is something we expected, given that in most of the cases despite the teleportation fidelity shows a sharper fall as d increases, the value of the classical limit decreases even more abruptly as shown in Fig. 4(a).

Although in principle it is expected that the fidelity of teleportation tends to decrease for high noise fractions, Fortes and Rigolin [6] found a set of scenarios in which it is not the case, even with no change in the measurement basis it is possible to get fidelities above the classical limit for high noise fractions. The following part is devoted to explore these cases in detail.



F

FIG. 5. Fidelity of teleportation under a scenario in which the qudits are affected by dit-phase-flip, dit-flip and d -phase-flip noises respectively ( FP,F,P ), and d = 3, for configurations leading to values above the classical limit, 〈 F 〉 &gt; f C = 1 2. / It is possible to notice that the larger the noise fraction, the lower the fidelity of teleportation, as expected.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01097v4-figure-cc82026c-2ccf-485c-b660-29ae0c4869be.jpg)

## A. Fighting noise with noise in qudit teleportation

We have performed an exhaustive search of scenarios in which the addition of noise to the system leads to an enhancement in the execution of the protocol without carrying out a basis change, obtaining the following extremal scenarios: ( ∅ , F, F ), ( ∅ , P, P ), [( ∅ , F P, F P ), for d = 2 only] and commutations. Under these, whenever both associated noise fractions get values either lower or higher than the noise threshold p ∗ = 1 -1 /d , the fidelity attains values above the classical limit, with 〈 F 〉 → (2 d -1) / d ( 2 -1) as the noise fractions approach to 1. Notice though that given the dependency of the classical fidelity on d [ f C ∼ O ( d -1 )], thus the larger the dimensionality, the more negligible becomes the gap between 〈 F 〉 and f C in the regime of high noise fractions. Moreover, note that a perfect restoration in the fidelity ( 〈 F 〉 = 1), is achieved by d = 2 only. It happens because the error in one of the qubits is globally corrected by the action of the same kind of error in another part of the system, while such a correction can be well succeeded for high dimensional systems only probabilistically, due to the fact that in this case there are several possible final configurations (e.g. a qutrit | j 〉 may be flipped in two different ways: either | j ⊕ 〉 1 or | j ⊕ 〉 2 ). Furthermore in the limit p → 1, d = 2 is the only case in which there is just one Kraus operator acting on the state and in this way the purity of the subsystem is not altered. Along with the scenarios we have just described, there are some other instances for which even adding noise to the third qudit, the system is still able to exhibit fidelities beyond its classical value, as exposed in figure 6 for a system under the action of dit-flip noise independently on its three components. For this, the teleportation fidelity reaches a value above f C whenever either the noise



FIG. 6. Fidelity of teleportation for the scenario in which all the parts are affected by dit-flip noise, and d = 3. The axes represent the noise fraction associated to each qudit involved in the teleportation process. Note that fidelity of teleportation values above the classical limit are obtained either if all the noise fractions are lower or two of them are higher than the noise threshold p ∗ = 1 -1 /d .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01097v4-figure-9f0593a1-7279-4393-b715-4f4ad8d5f225.jpg)

fractions are all below or two of them are above a noise threshold p ∗ = 1 -1 /d . It is important to remark that the scenario ( P, P, P ) exhibits the same behaviour. The other relevant situations correspond to adding an arbitrary class of noise to the input's qudit and dit-flip or d -phase-flip noise on the qudits in the channel. For all these instances the fidelity of teleportation behaves in the same way. In particular, results corresponding to the scenario ( D,F,F ) are presented in figure 7. Note that there are two separated regions satisfying 〈 F 〉 &gt; f C : the usual corresponding to all noise fractions below the threshold, and on the other hand, that in which both noise fractions associated to dit-flip errors are above the threshold p ∗ = 1 -1 /d , extended throughout the p D axis up to a maximal depolarizing noise fraction p ∗ D equal to d/ d ( 2 -d +1). The whole scenarios are summarized in Table I, as well as the maximal noise fractions on the input's qudit and possible commutations among the involved parts.

Another interesting case observed in [6] (also in [7], using the singlet fraction instead [20]), is the scenario ( ∅ , AD, AD ), for d = 2. In this, the fidelity of teleportation has an asymptotic tendency to the classical limit when the noise fractions approach to their maximal value, as it can be observed in Fig. 8 for d = 3. Another feature we were able to infer is that the regions of noise parameters corresponding to 〈 F 〉 &lt; f C are relatively small when compared to other scenarios involving noise on two parts. As an illustration, for d = 2 it covers ∼ 24 44% out of the , whole configurations and exhibits a tendency to decrease with d , for instance it is only about 15 , 4% for d = 5.



FIG. 7. Fidelity of teleportation for the scenario of one qutrit affected by depolarizing and the other two by trit-flip-noise, for configurations leading to values above the classical limit. Note that there is a region in which the fidelity attains values above f C , even for large noise fractions. See table I for all possible noise configurations presenting this behaviour.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01097v4-figure-0b508f94-a85d-4718-bc42-6915197f1c1a.jpg)

## VI. CONCLUSION

We have carried out a characterization of the qudit teleportation protocol in realistic scenarios, generalizing the results recently obtained for the case of qubits [6]. Under this approach, errors may be introduced to the system either by imperfections in the preparation and measurements and/or by the unavoidable interaction with the environment. We have been able to establish relations among the fidelity of teleportation and the amount of entanglement in the channel for the noiseless scenario. We performed an exhaustive search of all possible kinds of behaviour in the fidelity of teleportation, finding four predominant sets: i) The most typical and intuitive may be described as a decay in the fidelity with the noise frac-

TABLE I. Scenarios in which the fidelity of teleportation 〈 F X,Y,Z 〉 attains values above the classical limit for p Y , p Z &gt; p ∗ = 1 -1 /d , and corresponding maximal noise fractions p ∗ X . The last column indicates the cases where the fidelity does not change upon commutation between affected qudits and noise, i.e. 〈 F X,Y,Z 〉 = 〈 F Z,Y,X 〉 = 〈 F Y,X,Z 〉 = · · · .

| X   | Y   | Z   | p ∗ X            | Commutations?   |
|-----|-----|-----|------------------|-----------------|
| P   | F   | F   | 1 /d             | Yes             |
| FP  | F   | F   | 1 /d             | Yes             |
| D   | F   | F   | d/ ( d 2 - d +1) | Yes             |
| AD  | F   | F   | - a              | No              |
| F   | P   | P   | 1 /d             | Yes             |
| FP  | P   | P   | 1 /d             | Yes             |
| D   | P   | P   | d/ ( d 2 - d +1) | Yes             |
| AD  | P   | P   | - a              | No              |

a p ∗ = 2 √ 2 -2, for d = 2.



FIG. 8. Fidelity of teleportation for a scenario in which all parts are affected by amplitude damping noise ( AD,AD,AD ) and d = 3, in function of the noise fractions associated to each qudit in the channel. Note that in the limit p I → 0, only a small region of noise parameters lead to fidelities below its classical value, f C = 1 2. / Furthermore, as the amount of noise in the channel increases, the fidelity approaches to its classical limit.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01097v4-figure-9058369c-6010-47f5-af2a-df26a71e263b.jpg)

tions, depicted in figure 5. ii) The scenario ( F, F, F ) or ( P, P, P ) (Figure 6), in which the presence of noise in one qudit leads to an effective partial correction of the same kind of error on a different part. iii) An intermediate behaviour between the previous scenarios, on the one hand a decay in the fidelity and on the other one a partial error correction, summarized in Table I and Figure 7. And iv) the corresponding to any kind of noise in the input and Amplitude damping in the channel, depicted in Figure 8. Furthermore it was possible to note that it is possible to partially correct errors by a basis change, for a very specific case ( P, ∅ ∅ , ). It is left as an open question whether considering special scenarios in which the nature of the interactions favours some sub-classes of noise, can lead to an improvement of the performance of the protocol of high dimensional quantum teleportation. As a final remark, it would be very interesting to find further strategies leading to the improvement of the protocol of high-dimensional quantum teleportation, such as the presence of non-local memory effects [36].

## Appendix A: Parametrization for pure qudit states

Any arbitrary pure state of a qudit | φ 〉 = ∑ d -1 j =0 α j | j 〉 may be parametrized as:

$$\alpha _ { j } = \begin{cases} \cos \theta _ { 0 } & \text{for $j=0$} \\ \sin \theta _ { 0 } & \dots & \sin \theta _ { j - 1 } \cos \theta _ { j } \ e ^ { i \phi _ { j } } & \text{for $1\leq j\leq d-2$} \\ \sin \theta _ { 0 } & \dots & \sin \theta _ { d - 2 } \ e ^ { i \phi _ { d - 1 } } & \text{for $j=d-1$,} \end{cases}$$

with 0

&lt; θ j ≤ π/ 2 and 0 &lt; φ j ≤ 2 π .

Under this parametrization the invariant volume ele-

ment dΓ d is given by [37-39]:

$$d \Gamma _ { d } = \sin ^ { 2 d - 3 } \theta _ { 0 } \dots \sin \theta _ { d - 2 } \cos \theta _ { 0 } \dots \\ \dots \cos \theta _ { d - 2 } d \theta _ { 0 } \dots d \theta _ { d - 2 } d \phi _ { 1 } \dots d \phi _ { d - 1 },$$

in a compact form:

$$d \Gamma _ { d } = \prod _ { j = 0 } ^ { d - 2 } \sin ^ { 2 d - 2 j - 3 } \theta _ { j } \cos \theta _ { j } \, d \theta _ { j } \, d \phi _ { j + 1 }. \quad ( A 1 ) \quad \text{(A)}$$

The total volume V d = ∫ dΓ d may be easily calculated and is equal to:

$$V _ { d } = \frac { \pi ^ { d - 1 } } { ( d - 1 )! }. \quad \quad \quad ( A 2 ) \quad \text{(s}$$

## Appendix B: Calculation of 〈 α α α α j ∗ k l ∗ m 〉

The calculation of 〈 α α α α j ∗ k l ∗ m 〉 = 1 V d ∫ dΓ d α α α α j ∗ k l ∗ m may be facilitated if we take into account some symmetries. First of all, note that the volume element dΓ d does not depend explicitly on the phases φ j . Moreover, the state coefficients α j are proportional to exp( iφ j ), then the only way in which the integration does not vanish is having both: the coefficient and its conjugate inside the argument in order to cancel the corresponding phases. In this way we must have: 〈 α α α α j ∗ k l ∗ m 〉 ∝ ( δ jk δ lm + δ jm kl δ ) . Let us determine the proportionality constant. For simplicity we only show calculations for 〈 | α 0 | 4 〉 , nevertheless as the generated states are uniformly distributed, then any choice is equivalent. The integration reads:

$$\left \langle \left | \alpha _ { 0 } \right | ^ { 4 } \right \rangle = \frac { ( d - 1 )! } { \pi ^ { d - 1 } } \int \prod _ { j = 0 } ^ { d - 2 } \sin ^ { 2 d - 2 j - 3 } \theta _ { j } \times$$

× cos θ j d θ j d φ j +1 cos 4 θ . 0

It is not hard to show that integrations of the kind above have the following solutions:

$$\text{$v$ numbers $v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{\end{v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{ $\v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{ $v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{\end{v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \end{v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \vec{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \vec{v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{$v$} \colon \text{\ end}$$

for n = 0 2 4 , , , · · · and m &gt; 0. Thus 〈 | α 0 | 4 〉 is reduced to

$$\left \langle | \alpha _ { 0 } | ^ { 4 } \right \rangle = 2 ^ { d - 1 } ( d - 1 )! \, I _ { 2 d - 3 } ^ { 4 } \, \prod _ { j = 1 } ^ { d - 2 } I _ { 2 d - 2 j - 3 } ^ { 0 }. \quad ( B 1 ) \quad \sum _ { m n } ^ { \i }$$

It is straightforward to see that I 0 2 d -2 j -3 = 1 2( d - -j 1) and I 4 2 d -3 = 1 ( d +1) ( d d -1) , then

$$\left \langle \left | \alpha _ { 0 } \right | ^ { 4 } \right \rangle = 2 ^ { d - 1 } ( d - 1 )! \frac { 1 } { ( d + 1 ) d ( d - 1 ) } \, \prod _ { j = 1 } ^ { d - 2 } \frac { 1 } { 2 ( d - j - 1 ) }. \quad \text{$2d_{telepo}$}$$

This expression reduces to

$$\left \langle \left | \alpha _ { 0 } \right | ^ { 4 } \right \rangle = \frac { 2 } { d ( d + 1 ) }.$$

Back to the general case, it is possible infer that the proportionality factor must be equal to 1 d d ( +1) . In this way we have

$$\langle \alpha _ { j } \alpha _ { k } ^ { * } \alpha _ { l } \alpha _ { m } ^ { * } \rangle = \frac { 1 } { d ( d + 1 ) } ( \delta _ { j k } \delta _ { l m } + \delta _ { j m } \delta _ { k l } ). \quad ( B 3 )$$

The result above is very useful in the calculation of reduced expressions for the average fidelity of teleportation (see next Appendix).

## Appendix C: General expressions for fidelity of teleportation

This appendix is devoted to present the derivation of general expressions for the fidelity of teleportation within noisy environments under two approaches. The first contemplates the cases in which the Kraus operators associated to the classes of noise involved are proportional to Weyl operators. In the second part we consider Kraus operators written in the standard computational basis in order to consider the cases in which amplitude damping noise may take place in any part of the system.

## 1. Weyl-like noises

The noise coefficient associated to the input qudit a jk may be expressed as a superposition of the contributions of each region in figure 3: a 0 , noiseless region (green); a f , flip region (blue); a p , phase flip region (yellow) and a c for the region of combination of flip and phase-flip (red). In this way, the squared noise coefficient reads

$$\max \, \omega, \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \   \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \text{$\varphi$} \, \colon \\ a _ { j k } ^ { 2 } & = a _ { 0 } ^ { 2 } \delta _ { j, 0 } \delta _ { k, 0 } + a _ { f } ^ { 2 } \delta _ { j, 0 } \sum _ { n = 1 } ^ { d - 1 } \delta _ { k, n } + \\ & \quad + a _ { p } ^ { 2 } \delta _ { k, 0 } \sum _ { m = 1 } ^ { d - 1 } \delta _ { j, m } + a _ { c } ^ { 2 } \sum _ { m, n = 1 } ^ { d - 1 } \delta _ { j, m } \delta _ { k, n }. \ ( C 1 ) \\ \text{Thrues $\,two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$} \, \text{$two$.}$$

Thus we have the following correspondences between noise and reduced coefficients: dit-flip: a p = a c = 0, phase-flip: a f = a c = 0, dit-phase-flip: a p = a f = 0 and a f = a p = a c for depolarizing.

After some steps, the fidelity of teleportation F = ∑ mn Tr {| φ 〉 〈 φ ρ | ˆ mn } takes the form:

$$F = \sum _ { \substack { j k m n \\ \mu \nu = 0 } } ^ { d - 1 } \alpha _ { m } \alpha _ { n } ^ { * } \beta _ { j \mu } \beta _ { k \mu } ^ { * } \omega _ { d } ^ { \mu ( n - m ) } \rho _ { k, k \oplus \nu, n \oplus \nu, j, j \oplus \nu, m \oplus \nu },$$

using equation (7) and assuming a channel initially prepared in a pure state | ψ 〉 = ∑ d -1 k =0 γ k | kk 〉 , the fidelity of teleportation holds:

$$F = \sum _ { \substack { j k \mu \nu p _ { 1 } p _ { 2 } p _ { 3 } \\ q _ { 1 } q _ { 2 } q _ { 3 } = 0 } } ^ { d - 1 } \alpha _ { j \oplus \phi _ { q } \circ \Omega _ { 3 } } \alpha _ { k \oplus \phi _ { q } \circ \Omega _ { 3 } } ^ { * } \alpha _ { j \oplus \phi _ { 1 } } \alpha _ { k \oplus q _ { 1 } } \beta _ { j } \beta _ { k \mu } ^ { * } \omega _ { d } ^ { ( k - j ) ( \mu + p _ { 1 } + p _ { 2 } + p _ { 3 } ) } \gamma _ { k \oplus \vee \oplus q _ { 2 } } \gamma _ { j \oplus \vee \oplus \phi _ { q } } ^ { * } a _ { p _ { 1 } q _ { 1 } } ^ { 2 } b _ { p _ { 2 } q _ { 2 } } ^ { 2 } c _ { p _ { 3 } q _ { 3 } } ^ { 2 },$$

where b p q 2 2 and c p q 3 3 are the noise coefficients corresponding to the channel qudits respectively. By using the result of Appendix B (eq. B3) and after calculations the average fidelity takes the form

$$\cdots \quad & \langle F \rangle = \frac { 1 } { d + 1 } \left \{ 1 + \frac { 1 } { d } \sum _ { \substack { j k \mu \nu p _ { 1 } p _ { 2 } p _ { 3 } \\ q _ { 1 } q _ { 2 } q _ { 3 } = 0 } } ^ { d - 1 } \beta _ { j \mu } \beta _ { k \mu } ^ { * } \omega _ { d } ^ { ( k - j ) ( \mu + p _ { 1 } + p _ { 2 } + p _ { 3 } ) } \gamma _ { k \oplus \nu \oplus q _ { 2 } } \gamma _ { j \oplus \nu \oplus q _ { 2 } } ^ { * } q _ { 2 } ^ { 2 } \beta _ { 1 } ^ { 2 } \beta _ { 2 q _ { 2 } } c _ { p _ { 3 } q _ { 3 } } ^ { 2 } \delta _ { q _ { 2 }, q _ { 1 } \oplus q _ { 3 } } \right \}. \\ - \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \cdots \quad & \end{$$

Note that for the noiseless case the noise coefficients read a p q j j = δ p , j 0 δ q j , 0 (the same for b p q j j and c p q j j ), then the fidelity reduces to the equation obtained in II A for a noise free environment, as expected. By using analogous expressions for the noise coefficients of the channel qudits b jk and c jk (Eq. C1), substituting into equation C3, and after some calculations the fidelity of teleportation becomes:

$$\text{expressions for the noise coefficients of the channel quidits } b _ { j k } \text{ and } c _ { j k } \text{ (Eq. C1), substituting into equation C3, and} \\ \text{for some calculations the fidelity of teleportation becomes:} \\ \left ( F \right ) & = \frac { 1 } { d + 1 } \left ( 1 + d \right ) \left ( b _ { p } ^ { 2 } \left [ a _ { 0 } ^ { 2 } c _ { 0 } ^ { 2 } + ( d - 1 ) a _ { j } ^ { 2 } c _ { j } ^ { 2 } \right ] + \left [ \left [ a _ { 0 } ^ { 2 } + ( d - 2 ) b _ { p } ^ { 2 } \right ] \left [ a _ { p } ^ { 2 } c _ { 0 } ^ { 2 } + a _ { 0 } ^ { 2 } c _ { p } ^ { 2 } + ( d - 1 ) \left ( a _ { j } ^ { 2 } c _ { c } ^ { 2 } + a _ { 2 } ^ { 2 } c _ { j } ^ { 2 } \right ) \right ] + \\ & + \left [ ( d - 2 ) b _ { 0 } ^ { 2 } + ( d ^ { 2 } - 3 d + 3 ) b _ { p } ^ { 2 } \right ] \left [ a _ { p } ^ { 2 } c _ { p } ^ { 2 } + ( d - 1 ) a _ { 2 } ^ { 2 } c _ { p } ^ { 2 } \right ] \right \} + \\ + d ( d - 1 ) \left \{ \left [ ( d - 2 ) b _ { j } ^ { 2 } + ( d ^ { 2 } - 3 d + 3 ) b _ { c } ^ { 2 } \right ] \left [ a _ { p } ^ { 2 } c _ { c } ^ { 2 } + a _ { 2 } ^ { 2 } c _ { p } ^ { 2 } + ( d - 2 ) a _ { 2 } ^ { 2 } c _ { c } ^ { 2 } \right ] + b _ { 2 } ^ { 2 } \left [ a _ { j } ^ { 2 } c _ { 0 } ^ { 2 } + a _ { 0 } ^ { 2 } c _ { j } ^ { 2 } + ( d - 2 ) a _ { 2 } ^ { 2 } c _ { j } ^ { 2 } \right ] + \\ & + \left [ b _ { j } ^ { 2 } + ( d - 2 ) b _ { c } ^ { 2 } \right ] \left [ a _ { c } ^ { 2 } c _ { 0 } ^ { 2 } + a _ { 0 } ^ { 2 } c _ { c } ^ { 2 } + a _ { j } ^ { 2 } c _ { p } ^ { 2 } + a _ { p } ^ { 2 } c _ { j } ^ { 2 } + ( d - 2 ) \left ( a _ { j } ^ { 2 } c _ { c } ^ { 2 } + a _ { 2 } ^ { 2 } c _ { j } ^ { 2 } \right ) \right ] \right \} + \\ & + ( b _ { 0 } ^ { 2 } - b _ { p } ^ { 2 } ) \left [ \left ( a _ { 0 } ^ { 2 } - a _ { p } ^ { 2 } \right ) \left ( c _ { 0 } ^ { 2 } - c _ { p } ^ { 2 } \right ) + ( d - 1 ) \left ( a _ { j } ^ { 2 } - a _ { c } ^ { 2 } \right ) \left ( c _ { j } ^ { 2 } - c _ { c } ^ { 2 } \right ) \right ] \left ( 1 + ( d + 1 ) f _ { Q } \right ) + \\ & + \left ( b _ { 2 } ^ { 2 } - b _ { 2 } ^ { 2 } \right ) \left [ \left ( a _ { 0 } ^ { 2 } - a _ { p } ^ { 2 } \right ) \left ( c _ { j } ^ { 2 } - c _ { c } ^ { 2 } \right ) + \left ( a _ { j } ^ { 2 } - a _ { 2 } ^ { 2 } \right ) \left ( c _ { 0 } ^ { 2 } - c _ { j } ^ { 2 } \right ) + ( d - 2 ) \left ( a _ { j } ^ { 2 } - a _ { c } ^ { 2 } \right ) \left ( c _ { j } ^ { 2 } - c _ { c } ^ { 2 } \right ) \right ] \right ] \right ), \\ \text{where } f _ { Q } \text{ is the quantum contribution to the fidelity of teleportation in the absence of noise and } \tilde { f } \text{ is related to the}$$

where f Q is the quantum contribution to the fidelity of teleportation in the absence of noise and ˜ is related to the f channel and measurement coefficients as:

$$\tilde { f } = \frac { 1 } { d } \sum _ { \substack { j k \mu \nu = 0 \\ q = 1 } } ^ { d - 1 } \beta _ { j \mu } \beta _ { k \mu } ^ { * } \omega _ { d } ^ { \mu ( k - j ) } \gamma _ { k \oplus \nu \oplus q } \gamma _ { j \oplus \nu \oplus q } ^ { * },$$

attaining its highest value d d ( -1) when the entanglement in the channel and measurements is maximal.

## 2. Kraus operators in the standard computational basis

Let us calculate the fidelity of teleportation. Substituting Eq. 10 in Eq. 3, we have

$$F = \sum _ { \substack { j k m n \mu \nu \\ n _ { 1 } n _ { 2 } p _ { 1 } p _ { 2 } \\ k _ { 1 } k _ { 2 } k _ { 3 } } } \alpha _ { m } \alpha _ { n } ^ { * } \alpha _ { n _ { 1 } } \alpha _ { p _ { 1 } } ^ { * } \beta _ { j \mu } \beta _ { k \mu } ^ { * } \gamma _ { n _ { 2 } } \gamma _ { p _ { 2 } } ^ { * } \omega _ { d } ^ { \mu ( n - m ) } a _ { k, n _ { 1 } } ^ { ( k _ { 1 } ) } b _ { k \oplus \nu, n _ { 2 } } ^ { ( k _ { 2 } ) } c _ { n \oplus \nu, n _ { 2 } } ^ { ( k _ { 1 } ) } a _ { j, p _ { 1 } } ^ { ( k _ { 1 } ) * } b _ { j \oplus \nu, p _ { 2 } } ^ { ( k _ { 2 } ) * } c _ { m \oplus \nu, p _ { 2 } } ^ { ( k _ { 3 } ) * }.$$

Analogously to the previous treatment, using the results of Appendix B and after some calculations, the average fidelity of teleportation holds

$$\langle F \rangle = \frac { 1 } { d ( d + 1 ) } \sum _ { \substack { j k m _ { 1 } \mu \nu \\ n _ { 1 } n _ { 2 } p _ { 2 } \\ k _ { 1 } k _ { 2 } k _ { 3 } } } \beta _ { j \mu } \beta _ { k _ { \mu } } ^ { * } \gamma _ { n _ { 2 } } \gamma _ { p _ { 2 } } a _ { k _ { n _ { 1 } } } ^ { ( k _ { 1 } ) } b _ { k _ { \oplus } \nu, n _ { 2 } } ^ { ( k _ { 2 } ) } b _ { k _ { \oplus } \nu, p _ { 2 } } ^ { ( k _ { 3 } ) * } c _ { m \oplus \nu, p _ { 2 } } ^ { ( k _ { 3 } ) * } \left [ a _ { j, n _ { 1 } } ^ { ( k _ { 1 } ) * } c _ { m \oplus \nu, n _ { 2 } } + \mu _ { d } ^ { ( n _ { 1 } - m ) } a _ { j, m } ^ { ( k _ { 1 } ) * } c _ { n _ { 1 } \oplus \nu, n _ { 2 } } \right ].$$

It is worth to mention that the expression for fidelity above (Eq. C6) reproduces the whole results of Fortes

and Rigolin [6] for the case of qubits ( d = 2).

## ACKNOWLEDGMENTS

It is a pleasure to thank M´rcio M. Cunha for enlighta ening discussions, valuable remarks and collaboration at

- [1] C. H. Bennett, G. Brassard, C. Cr´ epeau, R. Jozsa, A. Peres, and W. K. Wootters, Phys. Rev. Lett. 70 , 1895 (1993).
- [2] S. Pirandola, J. Eisert, C. Weedbrook, A. Furusawa, and S. L. Braunstein, Nature photonics 9 , 641 (2015).
- [3] K. Banaszek, Phys. Rev. A 62 , 024301 (2000).
- [4] B. G. Taketani, F. de Melo, and R. L. de Matos Filho, Phys. Rev. A 85 , 020301 (2012).
- [5] S. Oh, S. Lee, and H.-w. Lee, Physical Review A 66 , 022316 (2002).
- [6] R. Fortes and G. Rigolin, Phys. Rev. A 92 , 012338 (2015).
- [7] L. T. Knoll, C. T. Schmiegelow, and M. A. Larotonda, Phys. Rev. A 90 , 042332 (2014).
- [8] X. Hu, Y. Gu, Q. Gong, and G. Guo, Phys. Rev. A 81 , 054302 (2010).
- [9] M. M. Cunha, E. A. Fonseca, M. Moreno, and F. Parisio, Quantum Information Processing 16 , 254 (2017).
- [10] M. G. M. Moreno, A. Fonseca, and M. M. Cunha, Quantum Information Processing 17 , 191 (2018), arXiv:arXiv:1805.06517v1.
- [11] N. J. Cerf, M. Bourennane, A. Karlsson, and N. Gisin, Phys. Rev. Lett. 88 , 127902 (2002).
- [12] T. Durt, D. Kaszlikowski, J.-L. Chen, and L. C. Kwek, Phys. Rev. A 69 , 032313 (2004).
- [13] T. V´ ertesi, S. Pironio, and N. Brunner, Phys. Rev. Lett. 104 , 060401 (2010).
- [14] N. T. Islam, C. C. W. Lim, C. Cahall, J. Kim, and D. J. Gauthier, Science advances 3 , e1701491 (2017).
- [15] P. Skrzypczyk and D. Cavalcanti, Physical Review Letters 120 , 260401 (2018), arXiv:1803.05199.
- [16] X.-M. Hu, C. Zhang, B.-H. Liu, Y.-F. Huang, C.-F. Li, and G.-C. Guo, arXiv preprint arXiv:1904.12249 (2019).
- [17] Y.-H. Luo, H.-S. Zhong, M. Erhard, X.-L. Wang, L.C. Peng, M. Krenn, X. Jiang, L. Li, N.-L. Liu, C.-Y. Lu, A. Zeilinger, and J.-W. Pan, Phys. Rev. Lett. 123 , 070505 (2019).
- [18] M. Kues, C. Reimer, P. Roztocki, L. R. Cort´s, S. Sciara, e B. Wetzel, Y. Zhang, A. Cino, S. T. Chu, B. E. Little, et al. , Nature 546 , 622 (2017).
- [19] R. A. Bertlmann and P. Krammer, Journal of Physics A: Mathematical and Theoretical 41 , 235303 (2008).
- [20] M. Horodecki, P. Horodecki, and R. Horodecki, Physical Review A - Atomic, Molecular, and Optical Physics 60 , 1888 (1999), arXiv:9807091 [quant-ph].
- [21] R. Weinar, W. Laskowski, and M. Pawowski, Journal of Physics A: Mathematical and Theoretical 46 , 435301 (2013).
- [22] We produced random states uniformly distributed in the basis {| 00 〉 , | 11 〉 , · · · | d -1 , d - 〉} 1 . Note that this parametrization does not cover uniformly the space of

early stages of this project. Useful suggestions by Nadja Bernardes, Fernando de Melo, Rafael Chaves, Carlos Batista, Marcelo Terra Cunha and Fernando Parisio are also acknowledged. This work was supported by a postgraduate grant from Conselho Nacional de Desenvolvimento Cient´ ıfico e Tecnol´gico (CNPq). o Support from Instituto Nacional de Ciˆ encia e Tecnologia-Informa¸ c˜o a Quˆ antica (INCT-IQ) is also acknowledged.

pure states for two-qudits. Nevertheless for our purposes it is enough.

- [23] D. Cavalcanti, P. Skrzypczyk, and I. Supi´, Phys. Rev. ˇ c Lett. 119 , 110501 (2017).
- [24] G. Carvacho, F. Andreoli, L. Santodonato, M. Bentivegna, V. D'Ambrosio, P. Skrzypczyk, I. ˇ upi´, D. CavS c alcanti, and F. Sciarrino, Phys. Rev. Lett. 121 , 140501 (2018).
- [25] I. ˇ upi´, S c P. Skrzypczyk, and D. Cavalcanti, Physical Review A 99 , 032334 (2019), arXiv:1804.10612.
- [26] M. Horodecki, P. Horodecki, and R. Horodecki, Physical Review Letters 80 , 5239 (1998).
- [27] M. A. Nielsen and I. L. Chuang, Quantum computation and quantum information (Cambridge university press, 2010).
- [28] M. Ramzan, Quantum Information Processing 12 , 577 (2013).
- [29] D. Chrusci´ski and F. A. Wudarski, Physical Review A n -Atomic, Molecular, and Optical Physics 91 , 012104 (2015), arXiv:1408.1792.
- [30] P. Imany, J. A. Jaramillo-Villegas, M. S. Alshaykh, J. M. Lukens, O. D. Odele, A. J. Moore, D. E. Leaird, M. Qi, and A. M. Weiner, npj Quantum Information 5 , 59 (2019).
- [31] P. Gokhale, J. M. Baker, C. Duckering, N. C. Brown, K. R. Brown, and F. T. Chong, in Proceedings of the 46th International Symposium on Computer Architecture -ISCA '19 (ACM Press, New York, New York, USA, 2019) pp. 554-566, arXiv:1905.10481.
- [32] D. Miller, T. Holz, H. Kampermann, and D. Bruß, Physical Review A 98 , 052316 (2018), arXiv:1807.06030.
- [33] Daniel Gottesman, Chaos, Solitons &amp; Fractals 10 , 1749 (1999).
- [34] A. Dutta, J. Ryu, W. Laskowski, and M. ukowski, Physics Letters A 380 , 2191 (2016).
- [35] Here we refer to a protected qudit as one whose probability of being affected by any kind of noise is negligible in comparison with that corresponding to other components in the system.
- [36] E. M. Laine, H. P. Breuer, and J. Piilo, Scientific Reports 4 , 4620 (2014).
- [37] C. M. Caves, Measures and volumes for spheres, the probability simplex, projective Hilbert space, and density operators , Tech. Rep. (University of New Mexico, 2001).
- [38] I. Bengtsson and K. Zyczkowski, Geometry of quantum states: an introduction to quantum entanglement (Cambridge University Press, 2007).
- [39] K. ˙ yczkowski and H.-J. Sommers, Journal of Physics A: Z Mathematical and General 34 , 7111 (2001).
Alejandro Fonseca
High-dimensional quantum teleportation under noisy environments
We study the protocol of qudit teleportation using quantum systems subjected
to several kinds of noise for arbitrary dimensionality $d$. We consider four
classes of noise: dit-flip, $d$-phase-flip, dit-phase-flip and depolarizing,
each of them corresponding to a family of Weyl operators, introduced via Kraus
formalism. We derive a general expression for the average fidelity of
teleportation in arbitrary dimension $d$ for any combination of noise on the
involved qudits. Under a different approach we derive the average fidelity of
teleportation for a more general scenario involving the $d$-dimensional
generalization of amplitude damping noise as well. We show that all possible
scenarios may be classified in four different behaviours and discuss the cases
in which it is possible to improve the fidelity by increasing the associated
noise fractions. All our results are in agreement with previous analysis by
Fortes and Rigolin for the case of qubits (Phys. Rev. A, 92 012338, 2015).