## Capacities of quantum Markovian noise for large times

Omar Fawzi , Mizanur Rahaman , and Mostafa Taheri 1 1 1

1 Univ Lyon, Inria, ENS Lyon, UCBL, LIP, Lyon, France

February 27, 2025

## Abstract

Given a quantum Markovian noise model, we study the maximum dimension of a classical or quantum system that can be stored for arbitrarily large time. We show that, unlike the fixed time setting, in the limit of infinite time, the classical and quantum capacities are characterized by efficiently computable properties of the peripheral spectrum of the quantum channel. In addition, the capacities are additive under tensor product, which implies in the language of Shannon theory that the one-shot and the asymptotic i.i.d. capacities are the same. We also provide an improved algorithm for computing the structure of the peripheral subspace of a quantum channel, which might be of independent interest.

## 1 Introduction

Consider a quantum system that we would like to use for storing quantum or classical information. This system is affected by noise that we assume is Markovian. It is natural to ask what is the minimum error that can be achieved for storing log D (qu)bits of information for some fixed time t . This is a typical question studied in Shannon theory, but here we focus on the limit t → ∞ , i.e., the information should remain in the system for arbitrarily long times.

Building a quantum system able to store quantum information for large times is one of the important goals of quantum information theory and it has been studied from different angles. Quantum error correction gives a mechanism to actively preserve the quantum information in a system undergoing local noise [26]. Such methods can achieve much more than a memory and can be used for fault-tolerant quantum computation [10]. A related important area of research is the study of self-correcting (or passive) quantum memories [6], i.e., physical systems that are robust to different forms of imperfections including thermal noise.

In this paper, we study the question of quantum

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00116v2-figure-8f63dc17-e051-4352-b2d5-91d32be8489a.jpg)

time

Figure 1: (a) Passive error correction : The system evolves under noise ( T ) over time, starting with encoding ( E ) and ending with decoding ( D ). (b) Active error correction : Noise ( N ) and recovery ( R ) maps are applied periodically to maintain system integrity, from encoding ( E ) to decoding ( D ).

memory from an abstract perspective where given a dynamical quantum system, the objective is to characterize the maximum amount of information that can be stored in this system without placing restrictions on the encoding/decoding operations.

As an example, in the passive model, the noise T is applied at each time step and we would like to characterize how many qubits can be stored reliably for time t as a function of T and t ? Our framework can also model the active setting where a fixed recovery operation R is applied at each time step; see Fig 1 for an illustration. Note that in our modeling, T is not limited to representing undesirable noise; it can also be a model for an engineered system such as cat qubits (see Examples 3.2 and 3.4 for a discussion). Such questions were studied in the recent work [22] and they obtained, among other results about scrambling, conditions for the classical capacity to be zero as well as lower bounds on the classical capacity of ergodic channels.

Our results In this work, we focus on the setting of arbitrarily large time, i.e., t → ∞ and we characterize both the classical and quantum capacities for a fixed error δ in terms of the peripheral spectrum of the noise model (Theorem 2.3). In addition, we show that such capacities are additive for the tensor product of channels (Theorem 2.4) and

can be computed in polynomial-time in the dimension (Theorem 2.6). This algorithm for computing the structure of the peripheral subspace (or the fixed point subspace) of a quantum channel improves on previous works and might of independent interest. We note that the fact that we can efficiently characterize the infinite-time capacities of noisy channels is in contrast with other settings for which capacities or optimal success probabilities correspond to hard problems, such as maximum independent set [21] or maximum coverage problems [4].

Our analysis extends to infinite-dimensional Hilbert spaces, where we identify specific conditions under which our findings remain valid (see Proposition 2.10). This generalization broadens the applicability of our results to Markovian noise on continuous variable systems. We then illustrate our framework through examples in Section 3.

## 2 Results

For a Hilbert space H , we use the notation Tr( H ) for the trace class operators on H . When we do not need to make the Hilbert space explicit and when it has dimension d , we denote the algebra of linear operators (or d × d complex matrices) by M d ( C ) or simply M d for short. We recall that a quantum channel E : M →M D d is a completely positive and trace preserving linear map [28].

We start with the standard definition of an error correcting code for a noisy quantum channel T . Note that throughout the paper d ≥ 1 is an integer and T is a quantum channel from M d to M d .

A quantum channel E : M D ( C ) → M d ( C ) is a ( D,δ ) classical code for a channel T : M d ( C ) → M d ( C ) if there exists a recovery channel R : M d ( C ) → M D ( C ) such that the average fidelity of the channel R◦ T ◦ E over all diagonal density matrices in M D is at least 1 -δ , i.e.,

$$\frac { 1 } { D } \sum _ { i = 1 } ^ { D } \langle i | \mathcal { R } \circ \mathcal { T } \circ \mathcal { E } ( | i \rangle \langle i | ) | i \rangle \geq 1 - \delta, \quad ( 1 ) \begin{array} { c } \text{for} \\ \text{space} \\ \text{such} \end{array}.$$

where {| 〉} i D i =1 is a fixed orthonormal basis of C D . We say that E is a ( D,δ ) quantum code if there exist channel R such that the entanglement fidelity of R◦T ◦ E is at least 1 -δ , i.e.,

$$\langle \Phi ^ { + } | \left ( \mathcal { I } \otimes \left ( \mathcal { R } \circ \mathcal { T } \circ \mathcal { E } \right ) \right ) \left ( | \Phi ^ { + } \rangle \langle \Phi ^ { + } | \right ) | \Phi ^ { + } \rangle \geq 1 - \delta, \begin{array} { c } d _ { k } & = \\ ( 2 ) & H _ { k, 2 } & [ \frac { \epsilon } { 2 } \\ \text{where } | \Phi ^ { + } \rangle & = & \frac { 1 } { \pi } \sum _ { i = 1 } ^ { D } | i \rangle \otimes | i \rangle \, \in \, ( \mathbb { C } ^ { D } ) ^ { \otimes 2 } \text{ and } \mathcal { I } & \text{Theor} \end{array}$$

where | Φ + 〉 = 1 √ D ∑ D i =1 | i 〉 ⊗ | 〉 i ∈ ( C D ) ⊗ 2 and I is the identity quantum channel on the reference system of dimension D .

Definition 2.1 (Capacity) . The (one-shot) classical capacity for a channel T is defined as

$$C _ { \delta } ( \mathcal { T } ) = \sup _ { \varepsilon } \log D,$$

where the supremum is over all possible ( D,δ ) classical codes E for T . Similarly, the (one-shot) quantum capacity for the channel T is defined as

$$Q _ { \delta } ( \mathcal { T } ) = \sup _ { \varepsilon } \log D,$$

where the supremum is taken over all possible ( D,δ ) quantum codes.

Remark 2.2 . This definition is concerned about passive error correction where no recovery is allowed at regular intervals. One could define an active capacity as follows:

$$Q _ { \delta } ^ { \text{active,t} } ( \mathcal { T } ) & = \log \{ \max D \, \colon \exists \mathcal { E }, \mathcal { R } _ { 1 }, \dots, \mathcal { R } _ { t } \text{ s.t.} \\ & \quad F _ { E } ( \mathcal { R } _ { t } \circ \mathcal { T } \circ \dots \circ \mathcal { T } \circ \mathcal { R } _ { 1 } \circ \mathcal { T } \circ \mathcal { E } ) \geq 1 - \delta \}, \\ \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

where E : M → M D d and R i : M → M d d for 1 ≤ i ≤ t -1 and R t : M → M d D are quantum channels. Note that the active capacity can be larger than the passive one; in particular, it is simple to see that for any t ≥ 1, Q active ,t 0 ( T ) = Q 0 ( T ) by choosing the recovery maps R i to decode and re-encode information. The works [17, 8, 19] have studied some Shannon-theoretic aspects of active error correction.

The infinite-time classical or quantum capacity of a quantum channel T is defined as C ∞ δ ( T ) = lim t →∞ C δ ( T t ), and Q ∞ δ ( T ) = lim t →∞ Q δ ( T t ) .

As we will demonstrate in the following, the peripheral subspace of a channel is directly related to its infinite-time capacity. The peripheral subspace of a quantum channel T is defined as follows:

$$\text{density} \ \chi _ { \mathcal { T } } = \text{span} \{ X \in \mathcal { M } _ { d } \, | \ \mathcal { T } ( X ) = \lambda X, \, | \lambda | = 1, \, \lambda \in \mathbb { C } \}.$$

For peripheral subspace χ T there exists Hilbert space decomposition C d = H 0 ⊕ ⊕ K k =1 ( H k, 1 ⊗ H k, 2 ) such that

$$\chi _ { \mathcal { T } } = 0 \oplus \bigoplus _ { k = 1 } ^ { K } \mathcal { M } _ { d _ { k } } \otimes \omega _ { k }, \quad \ \ ( 3 ) \\ \dots \cdots \cdots e \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{array}$$

where M d k is the full matrix algebra on H k, 1 , with d k = dim H k, 1 , and ω k is a density operator on H k, 2 [28].

Theorem 2.3. Let T : M → M d d be a quantum channel and let K and { d k } K k =1 be the integers from the peripheral subspace decomposition in (9) .

For any δ ∈ [0 , 1) , the infinite time classical capacity of T is given by

$$C _ { \delta } ^ { \infty } ( \mathcal { T } ) = \log \left ( \left \lfloor \frac { \sum _ { k } d _ { k } } { 1 - \delta } \right \rfloor \right ), \quad \begin{array} { c } \text{the } \alpha \\ ( 4 ) & \text{scrib} \\ \text{form} \\ \text{the $I$} \end{array} \end{array}$$

and the quantum capacity can be bounded as

$$\log \left ( \left \lfloor \frac { \max _ { k } d _ { k } } { \sqrt { 1 - \delta } } \right \rfloor \right ) \leq Q _ { \delta } ^ { \infty } ( \mathcal { T } ) \leq \log \left ( \frac { \max _ { k } d _ { k } } { 1 - \delta } \right ) \cdot \sum _ { ( 5 ) } ^ { m a c r o } \, R e m a r$$

We note that, in independent work, the special case δ = 0 of (4) and (5) were proved in [23].

Continuous-time Markovian noise T t are known as quantum Markov semigroups (QMS) generated by a Lindblad operator L , i.e., T t = e L t [20, 9, 16]. Similarly, the results presented in Theorem 2.3 hold for a QMS. Furthermore, the peripheral subspace of a QMS for any time can be expressed in terms of the spectrum of the generator L as follows:

$$\chi _ { e } \varsigma _ { t } = \text{span} \{ X \in \mathcal { M } _ { d } \, | \, \exists \theta \in \mathbb { R }, \, \mathcal { L } ( X ) = i \theta X \}. \quad \text{enc}$$

For more details check Proposition A.4.

One of the important questions regarding any capacity is additivity under tensor product. Since the infinite-time capacity is determined by the peripheral subspace, it is essentially additive under the tensor product.

Theorem 2.4. For any two quantum channels T and S , the infinite-time zero-error classical and quantum capacities are additive under the tensor product. Specifically, we have:

$$C _ { 0 } ^ { \infty } \left ( \mathcal { T } \otimes \mathcal { S } \right ) & = C _ { 0 } ^ { \infty } \left ( \mathcal { T } \right ) + C _ { 0 } ^ { \infty } \left ( \mathcal { S } \right ), \\ Q _ { 0 } ^ { \infty } \left ( \mathcal { T } \otimes \mathcal { S } \right ) & = Q _ { 0 } ^ { \infty } \left ( \mathcal { T } \right ) + Q _ { 0 } ^ { \infty } \left ( \mathcal { S } \right ).$$

$$Q _ { 0 } ^ { \infty } \left ( \mathcal { T } \otimes \mathcal { S } \right ) = Q _ { 0 } ^ { \infty } \left ( \mathcal { T } \right ) + Q _ { 0 } ^ { \infty } \left ( \mathcal { S } \right ).$$

For the asymptotic scenario, considering the limit of the infinite tensor product of a channel, the following results hold:

Proposition 2.5. Let δ ∈ [0 , 1) and T : M →M d d be a quantum channel. Then and

$$\lim _ { m \to \infty } \frac { 1 } { m } C _ { \delta } ^ { \infty } \left ( \mathcal { T } ^ { \otimes m } \right ) = C _ { 0 } ^ { \infty } ( \mathcal { T } ),$$

$$\lim _ { m \to \infty } \frac { 1 } { m } Q _ { \delta } ^ { \infty } \left ( \mathcal { T } ^ { \otimes m } \right ) = Q _ { 0 } ^ { \infty } ( \mathcal { T } ). \\ \text{compute the infinite-time capacities of a.}$$

To compute the infinite-time capacities of a given channel T , it is crucial to determine the structure of its peripheral subspace. An algorithm is proposed in [5] that, in polynomial time, maps the structure of the peripheral subspace to the structure of a finitedimensional von Neumann algebra. Consequently, our algorithm can also be utilized for determining the structure of finite-dimensional von Neumann algebras.

Theorem 2.6. Let T : M → M d d be a quantum channel, and let χ T = 0 ⊕ ⊕ K k =1 M ⊗ d k ω k be the decomposition of its peripheral subspace as described in (3) .Algorithm 1, takes the super-operator form of T as input and returns a representation of the Hilbert space decomposition and the fixed density matrices ω k in time O d ( 6 log d ) .

Remark 2.7 . An implementation of this algorithm with examples can be found in [24].

## 2.1 Infinite-dimensional Hilbert spaces

So far, we have focused on finite-dimensional Hilbert spaces. Now, let us shift our attention to the case of infinite-dimensional separable Hilbert spaces. Some concepts from finite-dimensional spaces do not extend directly to infinite dimensions, introducing new challenges and subtleties. To illustrate the differences in infinite-dimensional spaces, let us consider an example:

Example 2.8. Let {| 〉} i i ∈ Z be an orthonormal basis for the Hilbert space H , and let Tr( H ) denote the trace class space, i.e., the space of compact operators on H with finite trace norm. Now, consider the quantum channel T : Tr( H ) → Tr( H ) defined by

$$\mathcal { T } ( X ) = U X U ^ { \dagger }, \\ \dots$$

where U = ∑ i ∈ Z | i +1 〉〈 | i is known as the bilateral shift operator. This operator shifts each basis state | i 〉 to | i +1 , effectively 'shifting the indices up by 〉 one.'

It is easy to see that the channel has an empty peripheral subspace [18]. Despite the absence of a peripheral subspace, this channel preserves the distinguishability of input states, meaning that tr( ρσ ) = tr( T ( ρ ) T ( σ )), and thus its capacity is infinite despite the lack of a peripheral subspace.

The example above makes it clear that the peripheral subspace is not, in general, the relevant subspace for determining long-time capacity in the infinite-dimensional case.

Instead, for an infinite-dimensional system, the isometric subspace of the channel T , denoted Λ T , provides suitable properties [18]. This isometric subspace is defined as

$$\Lambda _ { \mathcal { T } } \coloneqq \{ & x \in \text{Tr} ( H ) \, \colon \forall t \in \mathbb { N }, \\ & \| \mathcal { T } ^ { t } ( x ) \| _ { 2 } = \| \mathcal { T } ^ { * t } ( x ) \| _ { 2 } = \| x \| _ { 2 } \}.$$

Remark 2.9 . If T is a quantum channel on a finitedimensional Hilbert space, then its isometric subspace coincides with its peripheral subspace [18].

Proposition 2.10. Let T : Tr( H ) → Tr( H ) be a completely positive, trace-preserving (CPTP) map, where H is an infinite-dimensional separable Hilbert space. Assume

- · T be contractive in the operator norm.
- · Λ T is a finite-dimensional subspace (spanned by a finite number of generators).

Then, Λ T has the form in (3) .

As a consequence, the classical and quantum capacities of channel T satisfy:

$$\log \left ( \left \lfloor \frac { \sum _ { k } d _ { k } } { 1 - \delta } \right \rfloor \right ) \leq C _ { \delta } ^ { \infty } ( \mathcal { T } ), \quad \ \ ( 6 ) \ \min _ { \substack { T t \\ \text{sipat} } }$$

and

$$\log \left ( \left \lfloor \frac { \max _ { k } d _ { k } } { \sqrt { 1 - \delta } } \right \rfloor \right ) \leq Q _ { \delta } ^ { \infty } ( \mathcal { T } ). \quad \quad ( 7 ) \\.$$

## 3 Examples

In this section, we present three examples to illustrate our theorem, each highlighting a different type of quantum dynamic. Example 3.1 considers T as noise affecting the system. Example 3.2 examines an engineered dynamic designed to control the evolution. Finally, Example 3.4 explores a scenario where noise, engineered dynamics, and a recovery process interact to restore the system.

Example 3.1. Let T be a collective noise that simultaneously affects multiple qubits due to uniform interactions with the environment. For an n -qubit system, T is given by the Kraus operators { 1 √ 3 exp( iσ ( n ) j ) } j ∈{ x,y,z } , where

$$\sigma _ { j } ^ { ( n ) } & = \sum _ { k = 1 } ^ { n } I _ { 2 } ^ { \otimes k - 1 } \otimes \sigma _ { j } \otimes I _ { 2 } ^ { \otimes n - k }, \\ \sim \quad. \quad. \cdot \cdot \cdot \cdot \cdot$$

and σ j are Pauli matrices [13].

For a 4-qubit system, the only eigenvalue with unit modulus is 1 + 0 i (easily checked numerically), thus the peripheral subspace of T is same as its fixed point subspace, and has the structure M ⊕ M ⊗ 2 ( 3 I 3 ) ⊕ C I 5 [13]. Therefore, by Theorem 2.3 its longtime classical and quantum capacity are C ∞ δ ( T ) = log ⌊ 6 1 -δ ⌋ , and Q ∞ δ ( T ) /similarequal log ( 3 1 -δ ) .

Example 3.2 (Cat code) . In quantum information, cat codes are an approach to encode qubits in coherent superpositions of photon-number states, known as cat states . These states are particularly valuable for creating robust qubits that resist certain types of noise, which is crucial for quantum error correction [12, 1].

To construct an n -component cat code, the system must be engineered so that the system and environment exchange photons while preserving a fixed photon-number parity modulo n [12]. This evolution can be modeled by a QMS with a generator L , where the jump operator is given by L = a n -α I n , with a as the photon annihilation operator and α ∈ C as a parameter that controls the steady-state properties.

This dissipative process gradually drives the system toward the n -component cat code subspace [1]. As a result, cat codes exhibit resilience against common errors, such as dephasing, allowing them to maintain coherence over extended timescales.

The peripheral subspace of n -photon driven dissipation evolution is identical to its fixed-point subspace and is given by

$$\chi _ { e } c _ { t } = \text{span} \{ | \beta \rangle \langle \beta ^ { \prime } | \, \colon \, \beta ^ { n } = \beta ^ { \prime n } = \alpha ^ { n } \},$$

where | β 〉 is a coherent state with parameter β [1, 12]. Although the coherent states | β 〉 are not orthogonal to each other, the set of orthogonal states {| ψ i 〉} n -1 i =0 , defined as superposition of coherent states as

$$| \psi _ { i } \rangle = \sum _ { j = 0 } ^ { n - 1 } \omega ^ { i j } | \omega ^ { j } | \alpha | \rangle, \quad \quad ( 8 )$$

where ω = exp(2 πi/n ) , satisfies

$$\chi _ { e ^ { c } } = \text{span} \{ | \psi _ { i } \rangle \langle \psi _ { j } | \, \colon i, j = 0, \dots, n - 1 \}.$$

Thus, the peripheral subspace has the structure 0 ⊕ M ⊗ n 1.

Remark 3.3 . Although the Hilbert space of the cat code in Example 3.2 is infinite-dimensional, the coefficients of coherent states in the Fock basis decay exponentially. Thus, we can apply the algorithm in Theorem 2.6 by truncating the number of photons, achieving a good accuracy (see [24]).

Although Proposition 2.10 and Theorem 2.3 are not applicable for this family of Markov processes-due to the infinite dimensionality of the Hilbert space and the non-contractivity of the processes in the operator norm-the fact that this evolution drives system into χ L [1] allows us to leverage the results of Theorem 2.3. Consequently, we have C ∞ δ ( e L t ) = log (⌊ n 1 -δ ⌋) , and Q ∞ δ ( e L t ) /similarequal log ( n 1 -δ ) .

Example 3.4. Photon loss, a common noise in bosonic systems, is modeled as a Lindblad process with the jump operator a (the photon annihilation operator) [1, 12]. Photon loss disrupts the structure of the cat code by altering the photon-number

parity, which is a key property for the code's errorcorrecting capabilities.

To approximate this noise over a finite time interval, we consider a simplified model where, with probability p , no photon is lost, and with probability 1 -p , exactly one photon is lost. The noise channel N is described by the Kraus operators: N 1 = √ pI , N 2 = √ 1 -p ∑ ∞ m =1 | m -1 〉〈 m | , and N 3 = √ 1 - | 〉〈 | p 0 0 .

∑ Let the noise channel N and the recovery channel R act on the system at each time interval ˜ . t Then, the evolution of the system in each time interval ˜ t is given by:

A recovery channel R can be defined to counteract the effects of photon loss. The Kraus operators for R , in a simplified version of those given in [15], are: R 1 = ∑ ∞ m =0 | 2 m + 1 〉〈 2 m | and R 2 = ∞ m =0 | 2 m +1 〉〈 2 m +1 . |

$$\mathcal { T } = \mathcal { R } \circ \mathcal { N } \circ e ^ { \mathcal { L } \tilde { t } }, \\. \quad \dots \quad \sim \mathcal { N } \circ$$

where L is the generator of the QMS describing the system's dynamics.

Let us focus on the 4-component cat code setup ( n = 4). One can verify that

$$\chi _ { \mathcal { T } } = \text{span} \{ | \psi ^ { i } \rangle \langle \psi ^ { j } | \, \colon i, j \in \{ 1, 3 \} \},$$

where | ψ i 〉 are defined as in (8). As a result, instead of having zero capacity in the long-time limit due to noise, the system achieves a classical capacity of

$$C _ { \delta } ^ { \infty } ( \mathcal { T } ) & = \log \left \lfloor \frac { 2 } { 1 - \delta } \right \rfloor, \\ \cdots \cdots$$

and a quantum capacity satisfying

$$Q _ { \delta } ^ { \infty } ( \mathcal { T } ) \geq \log \left \lfloor \frac { 2 } { \sqrt { 1 - \delta } } \right \rfloor.$$

## 4 Conclusion

In summary, we have shown that in the setting of infinite time, the classical and quantum capacities are essentially given by the zero-error capacities and can be efficiently computed, unlike the fixed time setting. We see this result as a first step towards understanding channel capacities for quantum evolutions. It would be interesting to study the behavior of capacities for finite t both for active and passive error correction.

## Acknowledgements

We acknowledge support from the European Research Council (ERC Grant AlgoQIP, Agreement No. 851716). MR is supported by the Marie Sk/suppressodowskal Curie Fellowship from the European

Union's Horizon Research and Innovation programme, grant Agreement No. HORIZON-MSCA2022-PF-01 (Project number: 101108117). We would like to thank Nilanjana Datta for a talk at the 'Quantum information' conference at the SwissMAP Research Station in Les Diablerets on [22] which stimulated this work.

## References

- [1] Victor V. Albert. Bosonic coding: introduction and use cases. arXiv preprint arXiv:2211.05714 , 2022. Available at https://arxiv.org/abs/2211.05714 .
- [2] Ravi B Bapat and Tirukkannamangai ES Raghavan. Nonnegative matrices and applications . Number 64. Cambridge university press, 1997.
- [3] Ivan Bardet, Marius Junge, Nicholas Laracuente, Cambyse Rouz´ e, and Daniel Stilck Fran¸ ca. Group transference techniques for the estimation of the decoherence times and capacities of quantum markov semigroups. IEEE Transactions on Information Theory , 67(5):2878-2909, 2021.
- [4] Siddharth Barman and Omar Fawzi. Algorithmic aspects of optimal channel coding. IEEE Transactions on Information Theory , 64(2):1038-1045, 2017.
- [5] Robin Blume-Kohout, Hui Khoon Ng, David Poulin, and Lorenza Viola. Informationpreserving structures: A general framework for quantum zero-error information. Phys. Rev. A , 82:062306, Dec 2010.
- [6] Benjamin J. Brown, Daniel Loss, Jiannis K. Pachos, Chris N. Self, and James R. Wootton. Quantum memories at finite temperature. Rev. Mod. Phys. , 88:045005, Nov 2016.
- [7] Eric V Denardo. Periods of connected networks and powers of nonnegative matrices. Mathematics of Operations Research , 2(1):20-24, 1977.
- [8] Omar Fawzi, Alexander M¨ uller-Hermes, and Ala Shayeghi. A Lower Bound on the Space Overhead of Fault-Tolerant Quantum Computation. In Mark Braverman, editor, 13th Innovations in Theoretical Computer Science Conference (ITCS 2022) , volume 215 of Leibniz International Proceedings in Informatics (LIPIcs) , pages 68:1-68:20, Dagstuhl, Germany, 2022. Schloss Dagstuhl -LeibnizZentrum f¨r Informatik. u

- [9] Vittorio Gorini, Andrzej Kossakowski, and E. C. G. Sudarshan. Completely positive dynamical semigroups of N-level systems. Journal of Mathematical Physics , 17(5):821-825, 05 1976.
- [10] Daniel Gottesman. Surviving as a quantum computer in a classical world. Textbook manuscript preprint , 2016.
- [11] Ji Guan, Yuan Feng, and Mingsheng Ying. Decomposition of quantum markov chains and its applications. Journal of Computer and System Sciences , 95:55-68, 2018.
- [12] J´ er´mie Guillaud, Joachim Cohen, and Mazyar e Mirrahimi. Quantum computation with cat qubits. SciPost Phys. Lect. Notes , page 72, 2023.
- [13] John A Holbrook, David W Kribs, and Raymond Laflamme. Noiseless subsystems and the structure of the commutant in quantum error correction. Quantum Information Processing , 2:381-419, 2003.
- [14] R. A. Horn and C. R. Johnson. Matrix Analysis . Cambridge University Press, 2013.
- [15] Zaki Leghtas, Gerhard Kirchmair, Brian Vlastakis, Robert J. Schoelkopf, Michel H. Devoret, and Mazyar Mirrahimi. Hardwareefficient autonomous quantum memory protection. Phys. Rev. Lett. , 111:120501, Sep 2013.
- [16] Goran Lindblad. On the generators of quantum dynamical semigroups. Communications in mathematical physics , 48:119-130, 1976.
- [17] Alexander M¨ uller-Hermes, David Reeb, and Michael M Wolf. Quantum subdivision capacities and continuous-time quantum coding. IEEE Transactions on Information Theory , 61(1):565-581, 2014.
- [18] R. Olkiewicz. Environment-Induced Superselection Rules in Markovian Regime. Communications in Mathematical Physics , 208(1):245265, December 1999.
- [19] Stefano Pirandola. End-to-end capacities of a quantum communication network. Communications Physics , 2(1):51, 2019.
- [20] Angel Rivas and Susana F Huelga. Open quantum systems , volume 10. Springer, 2012.
- [21] Claude Shannon. The zero error capacity of a noisy channel. IRE Transactions on Information Theory , 2(3):8-19, 1956.
- [22] Satvik Singh, Mizanur Rahaman, and Nilanjana Datta. Zero-error communication, scrambling, and ergodicity. arXiv preprint arXiv:2402.18703v1 , 2024. Available at https://arxiv.org/abs/2402.18703v1 .
- [23] Satvik Singh, Mizanur Rahaman, and Nilanjana Datta. Zero-error communication under discrete-time markovian dynamics. arXiv preprint arXiv:2402.18703v2 , 2024. Available at https://arxiv.org/abs/2402.18703v2 .
- [24] Mostafa Taheri. Code for computing the peripheral structure of a quantum channel. Available at https://github.com/mostafataheri75 /structureof-the-peripheral-subspace.
- [25] Robert Tarjan. Depth-first search and linear graph algorithms. SIAM journal on computing , 1(2):146-160, 1972.
- [26] Barbara M. Terhal. Quantum error correction for quantum memories. Rev. Mod. Phys. , 87:307-346, Apr 2015.
- [27] John Watrous. The Theory of Quantum Information . Cambridge University Press, 2018.
- [28] Michael M. Wolf. Quantum channels and operations - guided tour, 2012. Graue Literatur.

## A Infinite-time capacities

As previously mentioned, the peripheral subspace of the channel plays a crucial role in determining the infinite-time capacity. Before diving into the proof, let us first highlight some of its key properties.

Proposition A.1. [28] The peripheral subspace of a quantum channel T is defined by

$$\chi _ { \mathcal { T } } = \text{span} \{ X \in \mathcal { M } _ { d } \, \colon \mathcal { T } ( X ) = \lambda X, \ s. t. \, | \lambda | = 1 \}. \quad \text{a time}$$

It satisfies the following:

- · There exists a direct sum decomposition of the Hilbert space C d = H 0 ⊕ ⊕ K k =1 ( H k, 1 ⊗ H k, 2 ) for some nonnegative integer K such that

$$\chi _ { \mathcal { T } } & = 0 \oplus \bigoplus _ { k = 1 } ^ { K } \mathcal { M } _ { d _ { k } } \otimes \omega _ { k }, \quad \ \ ( 9 ) \\ \dots \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad \ \text{rm} \quad \text{whe}$$

where M d k is the full matrix algebra on H k, 1 , with d k = dim H k, 1 , and ω k is a density operator on H k, 2 .

- · There exists unitaries U k acting on H k, 1 and a permutation π which permutes systems { 1 , . . . , K } having the same dimension such that for every X ∈ χ T of the form X = 0 ⊕ ⊕ K k =1 x k ⊗ ω k , we have

$$\mathcal { T } ( X ) = 0 \oplus \bigoplus _ { k = 1 } ^ { K } U _ { k } x _ { \pi ^ { - 1 } ( k ) } U _ { k } ^ { \dagger } \otimes \omega _ { k }.$$

## A.1 Proof of Theorem 2.3

Proof. We start with the achievability statement, i.e., lower bounds on the capacities. Using the decomposition of the space in Proposition A.1, C d = H 0 ⊕ ⊕ K k =1 ( H k, 1 ⊗ H k, 2 ), let {| e k,j 〉} d k j =1 be orthonormal bases of H k, 1 . For different H k, 1 having the same dimension, we will identify the orthonormal bases.

We define the classical code E as follows. Let σ : { 1 , . . . , ∑ K k =1 d k } → { ( k, j ) : k ∈ { 1 , . . . , K } , j ∈ { 1 , . . . , d k }} be an arbitrary bijection, then

$$\mathcal { E } ( | i \rangle \langle i | ) & = \begin{cases} \ | e _ { \sigma ( i ) } \rangle \langle e _ { \sigma ( i ) } | \otimes \omega _ { k } & \text{if $1\leq i\leq \sum_{k=1}^{K}d_{k}$} \\ | e _ { 1 } \rangle \langle e _ { 1 } | \otimes \omega _ { k } & \text{if $\sum_{k=1}^{K}d_{k}<i\leq D$} \cdot$ and $\mathcal{E}(|(k)$} \end{cases} \mathcal { E } ( | ( k ) \\ \text{For } \ 1 \ \leq \ i \ \leq \ \sum _ { k=1}^{K} d_{k}$, \ we have $\mathcal{ }E(|i\rangle \langle i | ) \ = } ^ { K } \quad \text{we get}$$

For 1 ≤ i ≤ ∑ K k =1 d k , we have E | 〉 〈 ( i i | ) = | e k,j 〉 〈 e k,j | ⊗ ω k ∈ B H ( k, 1 ⊗ H k, 2 ) with ( k, j ) = σ i ( ). As a result, applying T , we get

$$\mathcal { T } ( \mathcal { E } ( | i \rangle \langle i | ) ) & = U _ { \pi ( k ) } | e _ { k, j } \rangle \langle e _ { k, j } | \otimes \omega _ { \pi ( k ) } U _ { \pi ( k ) } ^ { \dagger } \\ & \quad \in B ( H _ { \pi ( k ), 1 } \otimes H _ { \pi ( k ), 2 } ),$$

using Proposition A.1. Composing t times the map T , we get

$$\mathbf n g \, \text{the} \\ \text{proof}, \quad & \mathcal { T } ^ { t } ( \mathcal { E } ( | i \rangle \langle i | ) ) = \\ \mathfrak s. \quad & U _ { \pi ^ { t } ( k ) } \cdots U _ { \pi ( k ) } | e _ { k, j } \rangle \langle e _ { k, j } | U _ { \pi ( k ) } ^ { \dagger } \cdots U _ { \pi ^ { t } ( k ) } ^ { \dagger } \otimes \omega _ { \pi ^ { t } ( k ) } \\ \intertext { a c e of } \quad & \in B ( H _ { \pi ^ { t } ( k ), 1 } \otimes H _ { \pi ^ { t } ( k ), 2 } ). \\ & \quad \vdots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \vdots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \end{array}$$

Note that we always have T t ( E | 〉 〈 ( i i | )) ∈ χ T . For a time t , we choose a recovery map R defined by

$$\mathcal { R } ( X ) = \sum _ { k = 1 } ^ { K } & \left [ V _ { k } U _ { \pi ( k ) } ^ { \dagger } \cdots U _ { \pi ^ { t } ( k ) } ^ { \dagger } \\ & \quad \text{tr} _ { H _ { \pi ^ { t } ( k ), 2 } } ( P _ { \pi ^ { t } ( k ) } X P _ { \pi ^ { t } ( k ) } ) \\ & U _ { \pi ^ { t } ( k ) } \cdots U _ { \pi ( k ) } V _ { k } ^ { \dagger } \right ], \\ \text{where } P _ { k } \text{ is the orthogonal projection onto } H _ { k. 1 } \otimes$$

where P k is the orthogonal projection onto H k, 1 ⊗ H k, 2 and V k = ∑ d k j =1 | σ -1 ( k, j ) 〉 〈 e k,j | . Clearly, an error occurs only if i &gt; ∑ K i =1 d K and as such the success probability is given by

$$\stackrel { \Phi _ { \Phi }, \Phi _ { \Phi } } { \stackrel { \Phi \Phi } { \stackrel { \Phi \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Psi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \left | \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \stackrel { \Phi } { \text. } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ) } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } { } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } _ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }. } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } { _ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ & \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \. \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ | \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ _ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ \\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\. \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \...........................................................................................................................................$$

As a result, if we choose D = ⌊ ∑ K k =1 d k 1 -δ ⌋ , we obtain the desired achievability.

Now let us construct a quantum code of dimension D satisfying the condition max ( γ k k ) { ∑ K k =1 γ 2 k D 2 } ≥ 1 -δ .

∑ In order to define our encoding map E , we first define the Kraus operators E k = ∑ γ k j =1 | e k,j 〉 〈 ( k, j ) | for k = 1 to K , and F i = | e (1 1) , 〉 〈 i | for i ∈ S ′ . Then for X ∈ M D we define E ( X ) = ∑ K k =1 ( E XE k † k ) ⊗ ω k + ∑ i ∈ S ′ F XF i † i ⊗ ω 1 . As a result, we have

Let γ k achieve the maximum. It is convenient to label the basis of C D by elements in S ∪ S ′ , where S = { ( k, j ) : k ∈ { 1 , . . . , K } , j = 1 ∈ { 1 , . . . , γ k }} and S ′ = { 1 , . . . , D -∑ K k =1 γ k } . By the condition K k =1 γ k ≤ D , the set S is of size at most D .

$$\mathcal { E } ( | ( k, j ) \rangle \langle ( k, j ^ { \prime } ) | ) = | e _ { k, j } \rangle \langle e _ { k, j ^ { \prime } } | \otimes \omega _ { k },$$

/negationslash and E | ( ( k, j ) 〉 〈 ( k , j ′ ′ ) ) = 0 for | k = k ′ . For i, i ′ ∈ S ′ we get

$$\text{ge} \\ \mathcal { E } ( | ( k, j ) \rangle \langle i | ) = 0, \\ \mathcal { E } ( | i \rangle \langle i | ) = | e _ { 1, 1 } \rangle \langle e _ { 1, 1 } | \otimes \omega _ { 1 } \quad, \text{ and} \\ \mathcal { E } ( | i \rangle \langle i ^ { \prime } | ) = 0 \text{ for } i \neq i ^ { \prime }.$$

$$\mathcal { E } ( | i \rangle \langle i ^ { \prime } | ) = 0 \text{ for } i \neq i ^ { \prime }.$$

/negationslash

Thus, applying T t , we get

$$( \mathcal { T } ^ { t } \circ \mathcal { E } ) ( | ( k, j ) \rangle \langle ( k, j ^ { \prime } ) | ) & = \underset { \begin{array} { c } U _ { \pi ^ { t } ( k ) } \cdots U _ { \pi ( k ) } | e _ { k, j } \rangle \langle e _ { k, j ^ { \prime } } | U _ { \pi ( k ) } ^ { \dagger } \cdots U _ { \pi ^ { t } ( k ) } ^ { \dagger } & \text{capa} \\ \otimes \omega _ { \pi ^ { t } ( k ) }. & \text{turn} \\ \mathbb { N } \circ \mathbb { T } ^ { t } \circ \mathcal { T } ^ { t } \sim \mathcal { E } / | / L _ { \dot { i } } \vee \vee = \mathbb { W } \circ \text{choose same} \end{array}$$

Note that T t ◦ E ( ( | k, j ) 〉 〈 ( k , j ′ ′ ) ) | ∈ χ T . We choose R as

$$\mathcal { R } ( X ) = \sum _ { k = 1 } ^ { K } & \left [ V _ { k } U _ { \pi ( k ) } ^ { \dagger } \cdots U _ { \pi ^ { t } ( k ) } ^ { \dagger } \\ & \quad t r _ { H _ { \pi ^ { t } ( k ), 2 } } ( P _ { \pi ^ { t } ( k ) } X P _ { \pi ^ { t } ( k ) } ) \\ & U _ { \pi ^ { t } ( k ) } \cdots U _ { \pi ( k ) } V _ { k } ^ { \dagger } \right ], \\ \text{where } P _ { k } \text{ is the orthogonal projection onto } H _ { k, 1 } \otimes$$

where P k is the orthogonal projection onto H k, 1 ⊗ H k, 2 and V k = γ k j =1 | ( k, j ) 〉 〈 e k,j | .

∑ Then we can compute the entanglement fidelity as

$$& \text{$I$hen we can compute the entanglement hdelity} \ \mathcal { B } _ { k }$} \\ & \text{$\left(\Phi^{+}|\left(\mathcal{I} \otimes \mathcal{R} \circ \mathcal{T}^{t} \circ \mathcal{E} \right) \left(\|Phi^{+})\left(\Phi^{+}|\right)|\Phi^{+} \right \rangle} \\ & \geq \frac { 1 } { D ^ { 2 } } \sum _ { ( k, j ), ( k ^ { \prime }, j ^ { \prime } ) \in S } \left [ \quad & \text{wl} \ \text{in} \ \text{di} \ \right \rangle \\ & \left \langle ( k, j ) | \mathcal { R } \circ \mathcal { T } ^ { t } \circ \mathcal { E } ( | ( k, j ) \right \rangle ( ( k ^ { \prime }, j ^ { \prime } ) | ) | ( k ^ { \prime }, j ^ { \prime } ) \right \rangle \right ] \ \text{$\left(\|$in} \ \text{di} \ \right \rangle \\ & = \frac { 1 } { D ^ { 2 } } \sum _ { k = 1 } ^ { K } \sum _ { j, j ^ { \prime } = 1 } ^ { \gamma _ { k } } 1 \ \text{$\left(\|$a$} \ \text{$\alpha$} \ \text{$\alpha$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rho$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \\ & \geq 1 - \gamma, \quad & \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \ \text{$\rfloor$} \\ & \text{which proves the desired result.} \quad & \text{$\mathcal{T}_{P}$} \\ & \text{We now move to the converse bound. For this, it} \\. \quad & \text{$\mathcal{ }$} \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdots$$

which proves the desired result.

We now move to the converse bound. For this, it is convenient to consider the peripheral projection channel T P [28, Proposition 6.3] which satisfies the following properties: [28]

- · T P ( M d ) = χ T
- · there exist an increasing sequence { t i } such that lim i →∞ T t i = T P
- · T P ( X ) = X for any X ∈ χ T .

Note that if E is a ( D,δ ) code for T , and ‖T T ‖ ′ /diamondmath ≤ η , where the diamond norm is defined as ‖T - T ′ ‖ /diamondmath = sup ρ ‖ I ⊗ ( ( T - T ′ ))( ρ ) ‖ 1 , then E is also a ( D,δ + ) code for η T ′ .

Consider a sequence t i such that lim i →∞ T t i = T P and let η i = ‖T t i - T P ‖ /diamondmath . Let E i be a ( D,δ ) code for the channel T t i . Then E i is also a ( D,δ + η i ) code for the channel T P . Thus C δ ( T t i ) ≤ C δ + η i ( T P ). Taking the limit i → ∞ , we have that sup δ &lt;δ ′ C δ ′ ( T P ) ≤ C ∞ δ ( T ) ≤ inf δ &gt;δ ′ C δ ′ ( T P ). The exact same result holds for the quantum capacity as well. It now suffices to find upper bounds on the capacities of the channel T P .

First, let us show that the classical and quantum capacities of T P and ∑ K k =0 A k ◦ T P are the same, where we define A k as follows. Define P 0 to be the orthogonal projector onto H 0 and A 0 ( X ) = P XP 0 0 . Then define P k to be the orthogonal projector onto H k, 1 ⊗ H k, 2 and A k : B ( C d ) → B H ( k, 1 ) by A k ( X ) = tr H k, 2 ( P XP k k ). Then ∑ K k =0 A ◦ T k P is clearly a quantum channel. As T P ( M d ) = χ T , ∑ K k =0 A ◦T k P = ∑ K k =1 A ◦T k P so ∑ K k =1 A ◦T k P is a quantum channel mapping B ( C d ) to ⊕ K k =1 B H ( k, 1 ). The inequality C δ ( ∑ K k =1 A k ◦ T P ) ≤ C δ ( T P ) is clear. For the other inequality, let E be a code for T P and R be a corresponding recovery map. We define the recovery map R ′ = ∑ K k =1 R ◦ B k where B k : B H ( k, 1 ) → B H ( k, 1 ⊗ H k, 2 ) and maps B k ( x k ) = x k ⊗ ω k . It is easy to see that

$$\mathcal { R } ^ { \prime } \circ \sum _ { k = 1 } ^ { K } \mathcal { A } _ { k } \circ \mathcal { T } _ { P } \circ \mathcal { E } = \mathcal { R } \circ \mathcal { T } _ { P } \circ \mathcal { E }, \\ \vdots \cdots \dots \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

which proves the desired statement. For this reason, in what follows, we assume that dim H 0 = 0 and dim H k, 2 = 1 and the space decomposes as C d = K k =1 H k, 1 .

$$\text{$iP\Psi\rho\rho_{k=1}^{\prime}_{H_{k,1}}$}, \quad \text{$iMus$,} \\ \text{this, it} \\ \text{jection} \\ \text{files the} \quad \frac { 1 } { D } \sum _ { i = 1 } ^ { D } \langle i | \mathcal { R } \circ \mathcal { T } _ { P } \circ \mathcal { E } ( | i \rangle \langle i | ) | i \rangle & \leq \frac { 1 } { D } \sum _ { i = 1 } ^ { D } \langle i | \mathcal { R } \left ( \bigoplus _ { k = 1 } ^ { K } I _ { H _ { k, 1 } } \right ) | i \rangle \\ & = \frac { 1 } { D } \text{ tr} \left ( \bigoplus _ { k = 1 } ^ { K } I _ { H _ { k, 1 } } \right ) \\ \text{$i$} \} \text{ such} \\ & = \frac { \sum _ { k = 1 } ^ { K } d _ { k } } { D }. \\ \\. \quad \cdots \widetilde { \longrightarrow } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\array} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ }
\text{
\text{$\quad \ } \, \Phi$}$$

⊕ Let E be a ( D,δ ) classical code for T P and R a corresponding recovery channel. Note that because T P ( M d ) = χ T , for any density operator ρ ∈ B ( C d ), there exists positive operators ρ k on H k, 1 , T P ( ρ ) = ⊕ K k =1 ρ k . In addition, as T P is trace-preserving ∑ k tr( ρ k ) = 1 and ρ k ≤ I H k, 1 , where I H k, 1 is the identity on H k, 1 . As a result, T P ( ρ ) ≤ ⊕ K k =1 I H k, 1 . Thus,

Thus, for any ( D,δ ) code, we should have ∑ K k =1 d k D ≥ 1 -δ which implies D ≤ ⌊ ∑ K k =1 d k 1 -δ ⌋ . This implies that C δ ( T P ) ≤ ⌊ ∑ K k =1 d k 1 -δ ⌋ . As a result C ∞ δ ( T ) ≤ ⌊ ∑ K k =1 d k 1 -δ ′ ⌋ for any δ ′ &gt; δ which gives the desired result.

Let us now move to the quantum capacity. Let E be a ( D,δ ) classical code for T P and R a corresponding recovery channel.

In order to compute the entanglement fidelity, recall that the entanglement fidelity of channel E with Kraus operators { E i } is given by F E ( E ) = ∑ i | tr( E i ) | 2 . Let us introduce Kraus operators for { E j } j for T P ◦E , { R i } i for R . As T P ( M d ) = χ T , the operator { P E k j } k,j are also Kraus operators for the map T P ◦E , where we recall that P k is the projector onto H k, 1 .

Then we have

$$F _ { E } ( \mathcal { R } \circ \mathcal { T } _ { P } \circ \mathcal { E } ) = \frac { 1 } { D ^ { 2 } } \sum _ { i, j, k } \left | \text{tr} ( R _ { i } P _ { k } E _ { j } ) \right | ^ { 2 }. \\.$$

Let us denote α 2 k = ∑ i,j | tr( R P E i k j ) | 2 . We show two properties on α k : α 2 k ≤ d β k k with ∑ k β k = D and α 2 k ≤ d 2 k . We start with the first property

$$\alpha _ { k } ^ { 2 } & \leq \sum _ { i, j } \text{tr} ( R _ { i } P _ { k } P _ { k } ^ { \dagger } R _ { i } ^ { \dagger } ) \text{ tr} ( E _ { j } ^ { \dagger } P _ { k } ^ { \dagger } P _ { k } E _ { j } ) \\ & = \text{tr} ( P _ { k } ) \sum _ { j } \text{tr} \left ( E _ { j } ^ { \dagger } P _ { k } E _ { j } \right ). \\ \text{For the first incoming} \text{ } we need the Cambry-Software}$$

For the first inequality, we used the Cauchy-Schwarz inequality | tr( A B † ) | 2 ≤ tr( A A † ) tr( B B † ). Then we used the fact that R is trace preserving. Now note that because T P ◦ E is a quantum channel, we have ∑ k,j tr ( E P P E † j † k k j ) = D . Calling β k = ∑ j tr ( E P P E † j † k k j ) , we proved the first claimed inequality.

To show the second inequality α 2 k ≤ d 2 k , we write

$$\alpha _ { k } ^ { 2 } & \leq \sum _ { i, j } \text{tr} ( E _ { j } R _ { i } P _ { k } P _ { k } ^ { \dagger } R _ { i } ^ { \dagger } E _ { j } ^ { \dagger } ) \text{tr} ( P _ { k } ^ { \dagger } P _ { k } ) \\ & = \text{tr} ( P _ { k } ) d _ { k } \\ & = d _ { k } ^ { 2 }.$$

where we used again the Cauchy-Schwartz inequality and the fact that T P ◦ E ◦ R is trace-preserving. Defining γ k = min( d , β k k ), we have that ∑ K k =1 γ k ≤ D and γ k ≤ d k and the entanglement fidelity can be bounded as

$$\stackrel { \dots } { \dots } F _ { E } ( \mathcal { R } \circ \mathcal { T } _ { P } \circ \mathcal { E } ) & \leq \frac { 1 } { D ^ { 2 } } \sum _ { k = 1 } ^ { K } \alpha _ { k } ^ { 2 } \\ & \leq \frac { 1 } { D ^ { 2 } } \sum _ { k = 1 } ^ { K } d _ { k } \gamma _ { k }, \\ \text{h concludes the proof of } (10).$$

which concludes the proof of (10).

Remark A.2 . For the quantum capacity, we do not have an exact expression for the capacity but upper and lower bounds that differ by roughly 1 2 log(1 -δ ). The upper and lower bounds we establish in the proof is slightly stronger and it is simpler to express in terms of the optimal error δ for a fixed code size D :

$$\left | \text{el} \ \mathcal { E } \quad D \colon \\ \text{s for} \quad \underset { ( \gamma _ { k } ) _ { k } } { \max } \frac { \sum _ { k = 1 } ^ { K } \gamma _ { k } ^ { 2 } } { D ^ { 2 } } \\ \text{r the} \quad \leq \sup \{ 1 - \delta \colon \exists ( D, \delta ) \, \text{quantum code for} \ \mathcal { T } \} \ \text{(10)} \\ \text{ctor} \quad \leq \max \frac { \sum _ { k = 1 } ^ { K } d _ { k } \gamma _ { k } } { D ^ { 2 } } \\ \quad \text{where the supremum is taken over integers $\gamma_{k}$ sat} \quad \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot$$

where the supremum is taken over integers γ k satisfying γ k ≤ d k and ∑ K k =1 γ k ≤ D . Note that it is simple to see that the optimal choice for γ k is simply to take γ 1 = d , . . . , γ 1 s = d s and γ s +1 = D -∑ s k =1 γ k where s = argmax { s : ∑ s k =1 d k &lt; D } . When D is the sum of the largest r elements of ( d k ) k , then the upper and lower bounds match and we obtain an exact characterization. Observe that (5) follows from (10) by observing that max ( γ k k ) ∑ K k =1 γ 2 k ≥ min ( D , 2 max k d 2 k ) and ∑ K k =1 d γ k k ≤ (max k d k ) D . We leave the problem of computing the exact optimal fidelity as an open problem.

/negationslash

Remark A.3 (Classical special case) . Note that a classical stochastic d × d matrix M can be seen as a special case of a quantum channel T ( | i 〉 〈 i ′ | ) = 0 if i = i ′ and T ( | i 〉 〈 i | ) = ∑ j M j,i | j 〉 〈 j | . It is straightforward to verify that the eigenvalues of M are the eigenvalues of T and that the quantity ∑ K k =1 d k for such a channel is the number of eigenvalues of M of modulus 1 counted with multiplicity or the dimension of the peripheral subspace. This quantity even has a combinatorial interpretation as the sum of the periods of the bottom strongly connected components of the directed graph associated with the Markov chain M , as shown in [11]. Note that this combinatorial interpretation holds more generally for nonnegative matrices and such a decomposition into bottom strongly connected components is sometimes called the Frobenius normal form of the matrix [2, Section 1.7], see also e.g., in [14, Theorem 8.5.3 and Remark 8.5.4] for the period of each component. Using the combinatorial interpretation, one can find an algorithm running in linear time in the size of the graph for computing the capacity by using e.g., Tarjan's algorithm [25] to find the bottom strongly connected components and then find the periods of each component using [7].

## A.2 Quantum Markov semi-group's infinite-time capacities

Proposition A.4. Let ( T t ) t ≥ 0 be a quantum Markov semigroup. Then we have

$$\lim _ { t \to \infty } C _ { \delta } ( \mathcal { T } _ { t } ) = C _ { \delta } ^ { \infty } ( \mathcal { T } _ { 1 } ) \\ \lim _ { t \to \infty } Q _ { \delta } ( \mathcal { T } _ { t } ) = Q _ { \delta } ^ { \infty } ( \mathcal { T } _ { 1 } ).$$

$$\lim _ { t \to \infty } Q _ { \delta } ( \mathcal { T } _ { t } ) = Q _ { \delta } ^ { \infty } ( \mathcal { T } _ { 1 } ).$$

In addition, the peripheral subspace χ T 1 can be expressed in terms of the spectrum of generator L as follows: is the size of the block. Thus the super-operator of T 1 = e L has the form

$$\chi _ { \mathcal { T } _ { 1 } } = s p a n \{ X \in \mathcal { M } _ { d } | \exists \theta \in \mathbb { R }, \ \mathcal { L } ( X ) = i \theta X \}.$$

Remark A.5 . In Proposition A.4, the choice of t = 1 in T 1 is arbitrary and made for simplicity. The result holds for any fixed t 0 &gt; 0, with T 1 replaced by T t 0 .

We note that the behavior capacity of a special family of QMS (transferred QMS) over time has been thoroughly analyzed in [3], providing an asymptotic capacity evolution over time.

Proof. Let E be a ( D,δ ) classical (or quantum) code for T t , and R be its corresponding recovery channel. If t ′ ≤ t , then E is a ( D,δ ) classical (or quantum) code for T t ′ with recovery operator R◦T t -t ′ . Therefore, we have the following inequalities:

$$C _ { \delta } ( \mathcal { T } _ { \lfloor t \$$

and

$$Q _ { \delta } ( \mathcal { T } _ { [ t ] } ) \$$

In the limit as t → ∞ , both C δ ( T /floorleft t /floorright ) and C δ ( T /ceilingleft t /ceilingright ) converge to C ∞ δ ( T 1 ). Therefore, we have

$$e ^ { L } = V \left ( \bigoplus _ { \ell = 1 } ^$$

with

$$\begin{array} {$$

Thus the eigenvalues of T 1 are { e λ /lscript } /lscript . Now let X be an eigenvector of T 1 with eigenvalue e λ /lscript with | e λ /lscript | = 1. We know from [28, Proposition 6.2] that for such eigenvalues d /lscript = 1. Therefore, X is also an eigenvector of L with eigenvalue λ /lscript and λ /lscript is pure imaginary.

A simple corollary is that when taking tensor powers of a channel, as is common in Shannon theory, the classical and quantum capacities are given by the zero-error capacities at infinite time.

$$\lim _ { t \to \infty } C _ { \delta } ( \mathcal$$

Next, we examine the peripheral subspace of T 1 . It is straightforward to see that if L ( X ) = iθX for some operator X and real number θ , then T t ( X ) = e t L ( X ) = e iθt X , which implies X ∈ χ T t . Consequently, if X is spanned by eigenoperators of L corresponding to imaginary eigenvalues, then X belongs to χ T 1 .

To show the converse, we use the 'super-operator form' which represents the linear map E as a matrix ˆ acting on the vector space of operators. E Specifically, for a quantum channel E and a density matrix ρ , the action of E can be written as vec( E ( ρ )) = ˆ vec( E ρ ), where vec( ) denotes the column-stacking · vectorization of a matrix. Let L be the superoperator form of L , with a Jordan decomposition given by

$$L = V \left ( \bigoplus _ { \ell = 1 } ^ { s } J _ { \ell } ( \lambda _ { \ell } ) \right ) V ^ { - 1 }, \\ \therefore \cdots \cdots \cdots \cdots \cdots$$

where J /lscript ( λ /lscript ) is a Jordan block corresponding to an eigenvalue λ /lscript of L , i.e., J /lscript ( λ /lscript ) = λ I /lscript + J /lscript (0) with

$$& \text{eigenvalue} \ \lambda _ { \ell } \text{ of } L, \text{ i.e., } J _ { \ell } ( \lambda _ { \ell } ) = \lambda _ { \ell } I + J _ { \ell } ( 0 ) \text{ with } \\ & J _ { \ell } ( 0 ) = \begin{pmatrix} 0 & 1 & 0 & \dots & 0 \\ 0 & 0 & 1 & \dots & 0 \\ \vdots & \vdots & \ddots & \vdots & \vdots \\ 0 & \dots & \dots & \dots & 0 \end{pmatrix} \in \mathcal { M } _ { d _ { \ell } }, \text{ where } d _ { \ell }$$

. It has been demonstrated that the zero-error capacities of quantum channel T : M → M d d attain their infinite-time capacities after d 2 time concatenations[22]. Consequently, for any QMS or infinitely divisible quantum channel on finite Hilber space, we have:

Proposition A.6. Let ( T t ) t&gt; 0 be a quantum Markov semigroup. Then, for any ˜ t &gt; 0 , we have:

$$C _ { 0 } ( \mathcal { T } _ { \tilde { t } } ) & = C _ { 0 } ^ { \infty } ( \mathcal { T } _ { 1 } ), \\ Q _ { 0 } ( \mathcal { T } _ { \tilde { t } } ) & = Q _ { 0 } ^ { \infty } ( \mathcal { T } _ { 1 } ).$$

$$Q _ { 0 } ( \mathcal { T } _ { \tilde { t } } ) = Q _ { 0 } ^ { \infty } ( \mathcal { T } _ { 1 } ). \\ \cdot \quad t \in \in \mathfrak { m } \ \dots$$

Proof. It is shown in [22][Proposition 5.2] that any quantum channel T : M d → M d reaches its infinite-time classical and quantum zero-error capacities after t ≥ d 2 . Specifically:

$$\underset { \text{eq} - } { \longleftarrow } \sum C _ { 0 } ( \mathcal { T } ^ { d ^ { 2 } } ) = C _ { 0 } ^ { \infty } ( \mathcal { T } ) \quad \text{and} \quad Q _ { 0 } ( \mathcal { T } ^ { d ^ { 2 } } ) = Q _ { 0 } ^ { \infty } ( \mathcal { T } ).$$

Now, for any ˜ t &gt; 0, we can express T ˜ t as:

$$\mathcal { T } _ { \tilde { t } } = \left ( \mathcal { T } _ { \tilde { t } / d ^ { 2 } } \right ) ^ { d ^ { 2 } }. \\ \text{omposition, the zero}$$

Using this decomposition, the zero-error classical and quantum capacities of T ˜ t are equivalent to the infinite-time zero-error capacities of T ˜ t/d 2 . By applying Proposition A.4, we have:

$$C _ { 0 } ( \mathcal { T } _ { \tilde { t } } ) = C _ { 0 } ^ { \infty } ( \mathcal { T } _ { \tilde { t } / d ^ { 2 } } ) = C _ { 0 } ^ { \infty } ( \mathcal { T } _ { 1 } ),$$

and

$$Q _ { 0 } ( \mathcal { T } _ { \tilde { t } } ) = Q _ { 0 } ^ { \infty } ( \mathcal { T } _ { \tilde { t } / d ^ { 2 } } ) = Q _ { 0 } ^ { \infty } ( \mathcal { T } _ { 1 } ).$$

This completes the proof.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00116v2-figure-d5142047-7ad5-46c0-a9a2-a449002ed911.jpg)

## B Additivity of infinite-time capacities

Lemma B.1. Consider two quantum channels, T : M →M d d and S : M →M d ′ d ′ , with their respective peripheral projections denoted by T P and S P , and their respective peripheral subspaces denoted by χ T and χ S . For the tensor product of these channels, T ⊗ S , the peripheral projection is given by T P ⊗S P , and the corresponding peripheral subspace is χ T ⊗ χ S .

Proof. By definition there exist increasing sequences { m i } , { n i } , and { k i } such that lim i →∞ T m i = T P , lim i →∞ S n i = S P , and lim i →∞ ( T ⊗ S ) k i = ( T ⊗ S ) P . Since any power of a peripheral projection is itself, we have

$$\tilde { \ } m _ { i } n _ { i } k _ { i }$$

$$\text{und} \quad \lim _ { i \to \infty } ( \mathcal { T } \otimes \mathcal { S } ) ^ { m _ { i } n _ { i } k _ { i } } = ( \mathcal { T } \otimes \mathcal { S } ) _ { P } \,.$$

$$\lim _ { i \to \infty } \mathcal { T } ^ { m _ { i } n _ { i } k _ { i } } & = \mathcal { T } _ { P }, \\ \lim _ { i \to \infty } \mathcal { S } ^ { m _ { i } n _ { i } k _ { i } } & = \mathcal { S } _ { P }, \\ \text{and} \quad \lim _ { i \to \infty } ( \mathcal { T } \otimes \mathcal { S } ) ^ { m _ { i } n _ { i } k _ { i } } & = ( \mathcal { T } \otimes \mathcal { S } ) _ { P } \,.$$

Therefore, ( T ⊗ S ) P = T P ⊗S P .

Next, we consider the peripheral subspaces. The peripheral subspace is the fixed-point subspace of the peripheral projection. If X ∈ χ T ⊗ χ S , then it can be written as X = ∑ i Y i ⊗ Z i where Y i ∈ χ T and Z i ∈ χ S . So we have

$$\mathcal { T } _ { P } \otimes \mathcal { S } _ { P } ( X ) & = \sum _ { i } \mathcal { T } _ { P } ( Y _ { i } ) \otimes \mathcal { S } _ { P } ( Z _ { i } ) = \sum _ { i } Y _ { i } \otimes Z _ { i } = X, \\ \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \ \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

showing that X is a fixed point of T P ⊗S P . Therefore χ T ⊗ χ S ⊆ χ T ⊗S .

For the other direction, let X ∈ M ⊗ M d d ′ belongs to peripheral subspace of T ⊗ S , and decompose X as X = ∑ i Y i ⊗ Z i where Y i ∈ M d and Z i ∈ M d ′ . Since X is a fixed point of T P ⊗S P , we have

$$X = \mathcal { T } _ { P } \otimes \mathcal { S } _ { P } ( X ) = \sum _ { i } \mathcal { T } _ { P } ( Y _ { i } ) \otimes \mathcal { S } _ { P } ( Z _ { i } ). \\ \quad \cdots \cdots \cdots \cdots \cdots.$$

Because T P ( Y i ) ∈ χ T and S P ( Z i ) ∈ χ S , it follows that X ∈ χ T ⊗ χ S . Thus χ T ⊗S ⊆ χ T ⊗ χ S . Therefore, we conclude that χ T ⊗S = χ T ⊗ χ S .

## B.1 Proof of Theorem 2.4

Proof. By Lemma B.1, the peripheral subspace of T ⊗ S is the tensor product of the peripheral subspace of T and S . So if χ T = 0 ⊕ ⊕ K k =1 M ⊗ d k ω k for the decomposition C d = H 0 ⊕ ⊕ K k =1 H k, 1 ⊗ H k, 2 and χ S = 0 ⊕ ⊕ K ′ k ′ =1 M d ′ k ′ ⊗ ω ′ k ′ for the decomposition C d ′ = H ′ 0 ⊕ ⊕ K ′ k ′ =1 H ′ k , ′ 1 ⊗ H ′ k , ′ 2 , then we can decompose

$$\mathbb { C } ^ { d } \otimes \mathbb { C } ^ { d ^ { \prime } } & = \overline { H } _ { 0 } \oplus \\! s, \mathcal { T } \colon \quad & \bigoplus _ { k \in \{ 1, \dots, K \}, k ^ { \prime } \in \{ 1, \dots, K ^ { \prime } \} } H _ { k, 1 } \otimes H ^ { \prime } _ { k ^ { \prime }, 1 } \otimes H _ { k, 2 } \otimes H ^ { \prime } _ { k ^ { \prime }, 2 }, \\ \dots \quad & \mathcal { S } _ { P }, \\ \dots$$

where

$$\cdot \quad \text{where} \\ \cdot \quad \begin{array} { c } \overline { H } _ { 0 } = ( H _ { 0 } \otimes H ^ { \prime } _ { 0 } ) \\ \oplus \begin{pmatrix} H _ { 0 } \otimes ( \bigoplus _ { k ^ { \prime } = 1 } ^ { K ^ { \prime } } H ^ { \prime } _ { k ^ { \prime }, 1 } \otimes H ^ { \prime } _ { k ^ { \prime }, 2 } ) \\ \oplus \\ \oplus \begin{pmatrix} \overset { K } { ( \bigoplus _ { k = 1 } ^ { K } H _ { k, 1 } \otimes H _ { k, 2 } ) \otimes H ^ { \prime } _ { 0 } \end{pmatrix} \\ \end{array} \\ \quad \text{and get for this decomposition}$$

and get for this decomposition

$$\chi \mathcal { T } \otimes \mathcal { S } = 0 \oplus \bigoplus _ { k, k ^ { \prime } } \mathcal { M } _ { d _ { k } \times d _ { k ^ { \prime } } ^ { \prime } } \otimes \omega _ { k } \otimes \omega _ { k ^ { \prime } } ^ { \prime }.$$

By Theorem 2.3, C ∞ 0 ( T ⊗ S ) = log( ∑ k,k ′ d d k ′ k ′ ) = C ∞ 0 ( T ) + C ∞ 0 ( S ) and Q ∞ 0 ( T ⊗ S ) = log(max k,k ′ d d k ′ k ′ ) = Q ∞ 0 ( T ) + Q ∞ 0 ( S ).

## B.2 Proof of Proposition 2.5

Proof. By Theorem 2.3 and Theorem 2.4, we have

$$C _ { \delta } ^ { \infty } ( \mathcal { T } ^ { \otimes m } ) & = \log \left ( \left \lfloor \frac { 2 ^ { C _ { 0 } ^ { \infty } ( \mathcal { T } ^ { \otimes m } ) } } { 1 - \delta } \right \rfloor \right ) \\ & = \log \left ( \left \lfloor \frac { 2 ^ { m C _ { 0 } ^ { \infty } ( \mathcal { T } ) } } { 1 - \delta } \right \rfloor \right ).$$

But

$$\stackrel { \bullet } { \mathcal { S } } _ { P }, \, \text{we} \\ & \quad m C _ { 0 } ^ { \infty } ( \mathcal { T } ) + \log \left ( \frac { 1 } { 1 - \delta } \right ) - 1 \leq \log \left ( \left \lfloor \frac { 2 ^ { m C _ { 0 } ^ { \infty } } ( \mathcal { T } ) } { 1 - \delta } \right \rfloor \right ) \leq \\ _ { i } \right ). & \quad m C _ { 0 } ^ { \infty } ( \mathcal { T } ) + \log \left ( \frac { 1 } { 1 - \delta } \right ), \\ \text{follows} \\ & \quad \text{which browse the declared record} }$$

which proves the desired result.

For the quantum capacity, the same argument gives

$$\text{of} \quad m Q _ { 0 } ^ { \infty } ( \mathcal { T } ) + \frac { 1 } { 2 } \log \left ( \frac { 1 } { 1 - \delta } \right ) - 1 \leq Q _ { \delta } ^ { \infty } \left ( \mathcal { T } ^ { \otimes m } \right ) \\ \omega _ { k } \quad & \leq m Q _ { 0 } ^ { \infty } ( \mathcal { T } ) + \log \left ( \frac { 1 } { 1 - \delta } \right ), \\ \omega _ { k, 2 } \quad \text{with} \quad \text{the row one to avoid reading} \quad.$$

which proves the desired result.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00116v2-figure-d5142047-7ad5-46c0-a9a2-a449002ed911.jpg)

## C Algorithm

Proof of Theorem 2.6. The first step is to compute the peripheral projection channel T P . Let ˆ T ∈ M d 2 be the super-operator form of T . Using the Jordan normal form, we can express ˆ as T

- ˆ = T ∑ s /lscript =1 λ P /lscript /lscript + N /lscript , where λ /lscript are the eigenvalues of ˆ , T P /lscript are projections and N /lscript are nilpotent. The super-operator form of T P is given by [28, Proposition 6.3]

$$\hat { T } _ { P } = \sum _ { \ell \colon | \lambda _ { \ell } | = 1 } P _ { \ell }. \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots$$

Thus, ˆ T P can be computed using the Jordan normal form of ˆ which can be obtained in time T O (( d 2 ) 3 ) = O d ( 6 ).

Recall that χ T = Fix( T P ), where Fix( S ) is the fixed point space of a map S , i.e., Fix( S ) = span { X ∈ M d : S ( X ) = X } . The rest of the algorithm computes the structure of the fixed point space of the quantum channel T P .

$$\mathbb { C } ^ { d } = H _ { 0 } \oplus \bigoplus _ { k = 1 } ^ { K } H _ { k, 1 } \otimes H _ { k, 2 }. \quad ( 1 2 ) \quad \text{$\check{N}_{\text{c}$} \\ \text{ad of working} \sigma \text{ direct} l v \text{ with } \nu _ { \tau } \text{ it is convenient} \quad \underset { \text{$\check{d}$} \colon$}$$

Let χ T = Fix( T P ) = 0 ⊕ ⊕ K k =1 M ⊗ d k ω k for the Hilbert space decomposition

Instead of working directly with χ T , it is convenient to work with the matrix algebra A ⊆ B ( C d ) which is defined as

$$\mathcal { A } = \bigoplus _ { k = 1 } ^ { K } \mathcal { M } _ { d _ { k } } \otimes I _ { d _ { k } ^ { \prime } }. \quad \text{(13)} \quad \text{are} \ ] \\ \text{ilhert some decommission (12)} \text{ with } d \iota \, = \ \{ \mathcal { P } _ { k } \ \text{ }$$

for the Hilbert space decompostion (12) with d k = dim H k, 1 and d ′ k = dim H k, 2 . Then, χ T is called a distortion of the matrix algebra A = ⊕ K k =1 M ⊗ d k I d ′ k , i.e., there is a completely positive map D such that D A ( ) = χ T .

In [5, Theorem 5 and Section V], an algorithm is given to find operators A , . . . , A 1 N ∈ M d such that A = span( A , . . . , A 1 N ). This is done by computing left and right eigenvectors of ˆ and thus can be done T in time O d ( 6 ).

What remains is to determine the structure of A . To do this, it is sufficient to find a complete set of basis vectors | e k,i,j 〉 , where k = 1 , . . . , K , i = 1 , . . . , d k , and j = 1 , . . . , d k ′ , such that any element of A has a block diagonal matrix form 0 ⊕ ⊕ K k =1 A k ⊗ I d ′ k . In other words, we need to satisfy the following condition:

$$\langle e _ { k ^ { \prime }, i ^ { \prime }, j ^ { \prime } } | A | e _ { k, i, j } \rangle = \delta _ { k, k ^ { \prime } } \delta _ { j, j ^ { \prime } } \langle i ^ { \prime } | A _ { k } | i \rangle. \\ \ast \dots \dots \dots \dots \dots \dots \dots \dots.$$

Note that the choice of this basis is not unique because we can choose any arbitrary basis for M d k and M d k ′ , and form a basis for the algebra by taking their tensor product. To construct such a basis, the process involves three steps:

- 1. Find Minimal Central Projections for A : First, identify orthogonal projectors P k on the spaces H k, 1 ⊗ H k, 2 .
- 2. Find Minimal Projection for A : Next, decompose each projector P k as a sum of projections P k,i , i.e., P k = ∑ d k i =1 P k,i , where P k,i is a projection onto a subspace of the form | ψ i 〉 ⊗ H k, 2 , with | ψ i 〉 ∈ H k, 1 .
- 3. Construct the basis: Finally, use these projections to construct the required set of basis vectors.

This approach will yield the desired structure of the algebra A .

Note that for the matrix algebra A in (13), we have A ′ = ⊕ K k =1 I H k, 1 ⊗ M d ′ k and Z ( A ) = ⊕ K k =1 C I H k, 1 ⊗ I H k, 2 . Given an algebra B , a minimal projection in B is an orthogonal projection P such that P B P = C P [13]. The minimal projections of Z ( A ), that are also called the minimal central projections, are exactly the projectors P k we are looking for.

We start with the first step, i.e., computing P k . For that we use that the projectors P k are the minimal projections in the center of A . Recall that the center Z ( A ) = { X ∈ A : XY = Y X ∀ Y ∈ A} and that Z ( A ) = A∩A ′ where A ′ is the commutant of A defined by A ′ = { X ∈ B ( C d ) : XY = Y X ∀ Y ∈ A} .

Thus, in order to compute the projectors {P } k K k =1 , we first compute a representation of Z ( A ) as a linear span. We do this by computing a representation of the commutant A ′ as the linear span of some operators B , . . . , B 1 N ′ , which can be done by computing the kernel of the matrix Γ described in Lemma C.1. Then, the center can be computed by taking the intersection of subspaces A and A ′ ; see Algorithm 3. Then Algorithm 6 described a general algorithm to find all minimal projections of an algebra in M d in time O Nd ( 3 log d ), where N is the number of operators describing the algebra as a linear span. Note that the center Z ( A ) has dimension at most d and so we may assume N ≤ d . As such, we have computed all the {P } k K k =1 in time O d ( 5 log d ).

Now, the second step is within each block k , we want to compute minimal projections P k,i such that ∑ d k i =1 P k,i = P k where d k = dim H k, 1 . For that we now compute minimal projections P k,j in A satisfying P k,j P k = P k,j in algebra A by using Algorithm 6 with inputs P k and { A , . . . , A 1 N } . For each k , the runtime for finding the minimal projections { P k,j } k j =1 O d ( k × d 5 log d ) (as the dimension N of the algebra A can be up to d 2 ). As ∑ K k =1 d k ≤ d , the runtime of this step is O d ( 6 log d ).

For the last step we will construct the basis by using P k and P k,i . We have that d ′ k = rank( P k,j ) and d k is the number of minimal projections found { P k,j } j . From the structure of the algebra A (as given in (13)), we know that each minimal projection P k,i takes the form | α k,i 〉〈 α k,i | ⊗ I H k, 2 , where {| α k,i 〉} d k i =1 forms a complete basis for H k, 1 . Thus, the support of P k,i (eigenvectors with unite eigenvalues) is of form of {| α k,i 〉 ⊗ | β k,i,j 〉} , where {| β k,i,j 〉} d ′ k j =1 is a complete basis for H k, 2 . One can find these set of vectors by computing the eigenvalue and eigenvectors of P k,i . The last challenge is that although {| β k,i,j 〉} is a complete basis of H k, 2 for any i , they are not necessarily same for any two i and i ′ . To overcome this problem, we should find the unitary maps U k,m,n such | β k,m,j 〉 = U k,m,n | β k,n,j 〉 for all j .

For finding unitary map U k,m,n , we define matrix V k,m,n for A ∈ A as

$$V _ { i, j } ^ { k, m, n } \coloneqq \langle \alpha _ { k, n } | \langle \beta _ { k, n, j } | \mathcal { P } _ { k } A \mathcal { P } _ { k } | \alpha _ { k, m } \rangle | \beta _ { k, n, i } \rangle. \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

As A has the structure in form of ⊕ K k =1 A k ⊗ I H k, 2 , we have

$$V _ { i, j } ^ { k, m, n } = \langle \alpha _ { k, n } | A _ { k } | \alpha _ { k, m } \rangle \times \langle \beta _ { k, n, j } | \beta _ { k, m, i } \rangle. \quad \begin{smallmatrix} 1 \\ 1 \end{smallmatrix}$$

If 〈 α k,n | A k | α k,m 〉 is non-zero, then

$$U ^ { k, m, n } = \frac { V ^ { k, m, n } } { \text{tr} ( V ^ { k, m, n ^ { \dagger } } V ^ { k, m, n } ) }.$$

Since A is spanned by { A i } , for any k , m , and n , there exists at least one A ∈ { A i } such that V k,m,n is a non-zero matrix . This ensures that we can construct all of the U k,m,n .By using U k,m,n , we can construct the basis in the form of

$$| e _ { k, i, j } \rangle & \colon = | \alpha _ { k, i } \rangle \otimes | \beta _ { k, 1, j } \rangle \\ & = \sum _ { m } \langle \beta _ { k, i, m } | \beta _ { k, 1, j } \rangle | \alpha _ { k, i, m } \rangle | \beta _ { k, i, m } \rangle \\ & = \sum _ { m } U _ { j, m } ^ { k, 1, i } | \alpha _ { k, i } \rangle | \beta _ { k, i, m } \rangle. \\ & \quad \Pi _ { c }$$

Below we provide a lemma that was used in the above proof. For this lemma we need the concept of operator-vector correspondence: given a matrix X ∈ M d , it represents a vector, | X 〉〉 ∈ C n ⊗ C n .

Lemma C.1. Let A be a matrix algebra generated by { A , 1 · · · , A N } , and A ′ be the commutant of A . Then X belongs to A ′ if and only if | X 〉〉 belongs to the kernel of Γ , where

$$\Gamma = \sum _ { i = 1 } ^ { N } \left ( A _ { i } \otimes I - I \otimes A _ { i } ^ { T } \right ) ^ { \dagger } \left ( A _ { i } \otimes I - I \otimes A _ { i } ^ { T } \right ). \quad \text{Lem} _ { \substack { \text{Tr} (E)} \\ \text{in $t$} } } ^ { \text{oum} }$$

(14)

## Algorithm 1 Find the structure of χ T

- 2: Output: The structure of χ T as in (9)
- 1: Input: ˆ , T the super-operator form of T
- 3: procedure Peripheral-Subspace-Structure ( T )

/triangleright Jordan decomposition

- 5: ˆ T P ← ∑ /lscript : | λ /lscript | =1 P /lscript
- 4: ˆ = T ∑ s /lscript =1 λ P /lscript /lscript + N /lscript
- 6: return Fixed-Subspace-Structure ( ˆ T P ) /triangleright Algorithm 2 7: end procedure

## Algorithm 2 Find the fixed point structure Fix( S )

- 2: Output: The structure of Fix( S )
- 1: Input: ˆ , S the super-operator form of S
- 3: procedure Fixed-Subspace-Structure ( S )
- 5: { C , . . . , C 1 M } ← Center-of-algebra ( A , . . . , A 1 N ) /triangleright Algorithm 3
- 4: { A , . . . , A 1 N } ← Per-algebra-as-linear-span ( S ) /triangleright [5] repr. of A as a linear span
- 6: Ensure { A j } j and { C j } j are Hermitian via A → A + A † and A → i ( A -A † )
- 8: minimalCentralProj ← Minimal-projections ( P, { C i } ) /triangleright Algorithm 6
- 7: P ← projector on support of the center of A /triangleright P projector onto ⊕ K k =1 H k, 1 ⊗ H k, 2
- 9: for k ← 1 to K do

11:

- 10: P k ← minimalCentralProj[ k ]

minimalProj[

k

]

projections

←

(

P ,

k

A

i

)

Find-minimal-

{

}

/triangleright

Algorithm 6

- 12: end for 13: basisSet ← Construct-Basis ( { A i } , minimalCentralProj, minimalProj ) /triangleright Algorithm 7
- 14: Compute { ω k } as in [5, Lemma 5.4]
- 16: end procedure
- 15: return basisSet, {P } k

Proof. By the fact that A ⊗ B C | 〉〉 = | ACB T 〉〉 (see [27, Proposition 2.20]) applies if the dimensions of A B , , and C indicate ACB T is a valid matrix. So X ⊗ | I Y 〉〉 = | XY 〉〉 and I ⊗ X T | Y 〉〉 = | Y X 〉〉 . Therefore if we define

$$\text{can} \quad \Gamma _ { A _ { i } } \coloneqq \left ( A _ { i } \otimes I - I \otimes A _ { i } ^ { T } \right ) ^ { \dagger } \left ( A _ { i } \otimes I - I \otimes A _ { i } ^ { T } \right ), \\ \text{then an operator } X \text{ commutes with } A _ { i } \text{ if and only} \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdots$$

then an operator X commutes with A i if and only if Γ A i | X 〉〉 = 0. As Γ A i is positive, the kernel of Γ is the intersection of kernels of all Γ A i . Thus, X ∈ A ′ if and only if | X 〉〉 ∈ Γ.

Remark C.2 . The structure of the peripheral (fixed) subspace of collective noise for different numbers of qubits is summarized in Table 1. Additionally, the time required to compute this structure using the code provided in [24] on an Apple M2 Pro processor is included.

## D Infinite dimension

Before commencing the proof of Proposition 2.10, we present a proposition, and then we proceed to outline the proof.

Lemma D.1 (Theorem 19 of [18]) . Let T : Tr( H ) → Tr( H ) be a CPTP map and contractive in the operator norm. Then the isometric subspace

|   n | Dimension   | χ T                                  | Time       |
|-----|-------------|--------------------------------------|------------|
|   3 | 8 × 8       | CI 4 ⊕ ( M 2 ⊗ I 2 )                 | 18.576 ms  |
|   4 | 16 × 16     | M 2 ⊕ CI 5 ⊕ ( M 3 ⊗ I 3 )           | 647.805 ms |
|   5 | 32 × 32     | CI 6 ⊕ ( M 5 ⊗ I 2 ) ⊕ ( M 4 ⊗ I 4 ) | 4.763 s    |
|   6 | 64 64       | 5 CI 7 ( M 5 I 5 ) ( 9 I 3 )         | 223 s      |

×

M ⊕

⊕

M ⊗

⊕ M ⊗

Table 1: Structure of χ T , matrix sizes, and computation times for various number of qubit.

## Algorithm 3 Compute center of algebra A = span { A , . . . , A 1 N }

## Algorithm 5 Find one minimal projection in the range of projection P

- 1: Input: A , A , . . . , A 1 2 N d /triangleright = span( A , . . . , A 1 N )
- 2: Output: { C , . . . , C 1 M } that span the center of A
- ∈ M A
- 3: procedure Center-of-algebra ( A , A , . . . , A 1 2 N )
- 4: Γ ← ∑ i ( A A † i i ) ⊗ I -A † i ⊗ A T i -A i ⊗ A i + I ⊗ ( A A i T i ) /triangleright ¯ is the complex conjugate of A A
- 5: Compute the kernel of Γ as span { B , . . . , B 1 N ′ } /triangleright See Lemma C.1, span { B , . . . , B 1 N ′ } = A ′
- 7: Compute the kernel of Γ + Γ ′ as span { C , . . . , C 1 M } /triangleright Γ Γ , ′ ≥ 0, so ker Γ + Γ ′ = ker Γ ∩ ker Γ ′
- 6: Γ ′ ← ∑ i ( B B † i i ) ⊗ I -B † i ⊗ B T i -B i ⊗ B i + I ⊗ ( B B i T i ) /triangleright ker Γ ′ corresponds to A ′′ = A
- 8: return { C , . . . , C 1 M }
- 9: end procedure O Nd ( 4 + d 6 )

/triangleright The runtime of this procedure is

## Algorithm 4 Find projection smaller than P

- N
- 2: Output: Projection Q ∈ A such that QP = Q , and tr( Q ) ≤ tr( P / ) 2 if P is not minimal
- 1: Input: Orthogonal projection P and { A i } i =1 spanning A
- 3: procedure Reduce-Projection ( P, { A i } )
- 5: Write spectral decomposition PA P i = ∑ j λ P j j /triangleright At least 2 nonzero eigenvalues
- 4: if ∃ i is such that PA P i ∈ C P then /triangleright P is not minimal

/negationslash

6:

7:

- /triangleright [13] shows that P j ∈ A for all j

return

Q

= arg min tr(

Q

)

- ≤ 8: else

tr(

P /

)

2 as tr(

∑

{

j

P

j

- 9: Return P

10:

- 11: end procedure /triangleright The runtime of this procedure is O Nd ( 3 )

end if

Λ T is a T -invariant subspace and Λ T decomposes as

$$\Lambda _ { \mathcal { T } } = \bigoplus _ { k } \Lambda _ { \mathcal { T } } ^ { ( k ) }, \\ \Lambda _ { \mathcal { T } } \left ( k \right ) \, \underset { 1 \, \mathcal { H } } { \dots } \, \left ( k \right )$$

/negationslash where, for all φ ∈ Λ ( k ) T and ψ ∈ Λ ( ) l T , we have tr( ψφ ) = 0 for k = l . For each k , there exists a Banach space isomorphism α k : Λ ( k ) T → Tr( ˜ H k ) , such that the action of T on Λ ( k ) T corresponds to U † k · U k , where U k is a unitary operator on ˜ H k , i.e.

$$\alpha _ { k } \circ \mathcal { T } \circ \alpha _ { k } ^ { - 1 } (. ) = U _ { k } \,. \, U _ { k } ^ { \dagger } \\ \sim \wedge \ w \, \ w \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot$$

Remark D.2 . If Λ T is not finite, then the infinitetime classical capacity of T will be infinite. However, the infinite-time quantum capacity is not necessarily infinite. For example, consider the quantum channel T : Tr( H ) → Tr( H ) defined as:

$$\mathcal { T } ( X ) & = \sum _ { i \in \mathbb { Z } } K _ { i } X K _ { i } ^ { \dagger }, \\. & \models / \psi _ { i } \quad \mathbb { T } \quad \dots \quad \mathbb { T } \quad \colon = \mathbb { T } \quad \dots$$

where K i = | i 〉〈 i - | 1 , and {| 〉} i is the orthonormal basis of H . In this case, the infinite-time classical

P

j

)

/triangleright

We have

≤

P

/triangleright P is minimal

}

1: Input: Orthogonal projection P and { A i } spanning A

2: Output: A minimal projection Q ∈ A such that Q ≤ P

3: procedure Find-one-minimal-projection ( P, { A i } )

4: Q ← Reduce-Projection ( P, { A i } )

/triangleright

Algorithm 4

5: while P = Q do /triangleright As the trace is divided by 2 at each step, at most log d steps

6: P ← Q

7:

Q ← Reduce-Projection ( P, { A i } )

8: end while

9: return P

10: end procedure /triangleright Reduce-Projection is called at most

log d times

/negationslash

## Algorithm 6 Decomposing P to minimal projection in algebra A

/negationslash

1: Input: Orthogonal projection P and A i spanning

∈ A { } A such

2: Output: A set of minimal projections Q ,. . . , Q 1 s in A that Q 1 + · · · + Q s = P

3: procedure Find-minimal-projections ( P, { A i } )

4: minimalProjections ←{}

5:

while P = 0 do /triangleright Iterate until P

becomes zero

6: Q ← Find-one-minimal-projection ( P ) /triangleright Algorithm 5 7: Add Q to minimalProjections

8: P ← P - Q /triangleright Update P by removing the minimal projector that was found

9: /triangleright Note that P - Q is also an orthogonal projector end while

10:

11: return minimalProjections

12: end procedure

capacity is infinite ( C ∞ 0 ( T ) = ∞ ), but the infinitetime quantum capacity is zero ( Q ∞ 0 ( T ) = 0).

Proof of Proposition 2.10. We begin by showing that the isometric and peripheral subspaces are the same if Λ T is finite-dimensional. Afterward, we define a restricted channel ¯ T on a finite-dimensional Hilbert space that acts similarly to T on the isometric subspace. Finally, we use the result from Theorem 2.3 to determine the achievable capacity of T t in the limit as t →∞ .

For any CPTP map T , the peripheral subspace is always a subset of the isometric subspace (see Proposition 9 of [18]). To prove Λ T = χ T , it suffices to show that Λ T ⊆ χ T when Λ T is spanned by a finite number of generators.

By Lemma D.1, the isometric subspace Λ T can be decomposed into K subspaces { Λ ( k ) T } K k =1 , where each Λ ( k ) T is isomorphic to Tr( ˜ H k ), with ˜ H k a finitedimensional Hilbert space. Let {| i k 〉} d k i =1 be a complete basis of ˜ H k , where d k = dim( ˜ H k ), such that the unitary operator U k acts as U k | i k 〉 = e iθ i | i k 〉 ,

## Algorithm 7 Construct the Basis e k,i,j for A

- 1: Input : { A i } spanning A , Orthogonal Minimal Central Projections { P k } , Orthogonal Minimal Projections { P k,i } for each P k
- 2: Output : The basis vectors e k,i,j 3: procedure Construct-Basis ( { A i } { , P k } { , P k,i }
- )
- 4: K ← number of minimal central projections { P k }

5:

d k

[

]

number of minimal projections

P

k,i

6:

←

{

}

)

/triangleright Step 1: Compute Eigenvectors for Minimal d-prime[ k ] ← tr( P k, 1 Projections

7: for each P k,i where k = 1 , . . . , K and i = 1 , . . . , d [ k ] do 8: eigenVectors[ k ][ i ] ← of eigenvectors of P k,i with unit eigenvalues

9:

end for

/triangleright and

k

for

= 1

, . . . , K

Initialize

V

while

←

V

Pick

= 0

A

for

Step 2:

Construct matrices

n

0

do from

i, j

/triangleright

= 1

, . . . , d

U

[

k

]

do

Exit loop once a valid

{

A

V

i

= 1

, . . . ,

V i, j

[

]

eigenVectors[

k

][1][

i

]

end for end while

U k

[

][

n

]

V

V

tr(

†

←

V

)

}

d-prime[

k

]

do eigenVectors[

k

][

n j

][

]

†

·

A

·

/triangleright Normalize U k [ ][ n ]

19:

10:

11:

12:

13:

14:

15:

16:

17:

18:

←

end for

/triangleright Step 3: Construct the basis vectors e k,i,j 20: for k = 1 , . . . , K , i = 1 , . . . , d [ k ], and j = 1 , . . . , d-prime[ k ] do

21:

Initialize basisVector

0

22:

←

for

m

= 1

, . . . , d-prime[

k

]

do

23: basisVector ← basisVector + U k [ ][ i ][ j, m ] × eigenVectors[ k ][ i ][ m ]

24:

end for

25:

26:

27:

basisSet[

basisSet[

k

][

i

k

][

i

][

j

]

basisVector

][

j

]

end for

Return

The set basisSet

28:

end procedure for some phases θ i . Considering the action of U k on the operators | i k 〉〈 j k | , we have:

$$U _ { k } | i _ { k } \rangle \langle j _ { k } | U _ { k } ^ { \dagger } = e ^ { i ( \theta _ { i } - \theta _ { j } ) } | i _ { k } \rangle \langle j _ { k } |.$$

Using the isomorphism α -1 k that maps Tr( ˜ H k ) to Λ ( k ) T , we find:

$$\mathcal { T } \circ \alpha _ { k } ^ { - 1 } ( | i _ { k } \rangle \langle j _ { k } | ) & = \alpha _ { k } ^ { - 1 } ( U _ { k } | i _ { k } \rangle \langle j _ { k } | U _ { k } ^ { \dagger } ) \\ & = e ^ { i ( \theta _ { i } - \theta _ { j } ) } \alpha _ { k } ^ { - 1 } ( | i _ { k } \rangle \langle j _ { k } | ). \\ \dots \. \. \quad. \quad - 1 \dots \dots \dots \. \. \quad.$$

This shows that α -1 k ( | i k 〉〈 j k | ) belongs to the peripheral subspace χ T for all i, j . Since Tr( ˜ H k ) is spanned by the operators {| i k 〉〈 j k |} , it follows that Λ ( k ) T is spanned by { α -1 k ( | i k 〉〈 j k | ) } . Therefore, Λ ( k ) T ⊆ χ T . Summing over all k , we conclude that Λ T ⊆ χ T . Combined with the fact that χ T ⊆ Λ , T we conclude that Λ T = χ T .

Since we assumed Λ T is spanned by a finite number of generators, and we know that Λ T is generated by finite-dimensional orthogonal projections [18, Proposition 5],therefore, Λ T is generated by a finite number of finite-dimensional projections, i.e

$$\Lambda _ { \mathcal { T } } = \text{span} \{ P _ { i } \in \text{Tr} ( H ) ; i = 1, \cdots d \}. \quad ( 1 5 )$$

←

/triangleright Store in is

k,

1

,n found

Let us define S = ∑ d i =1 P i . We let H ′ be the support of S . As H ′ is finite-dimensional, there exists a subspace H ′′ such that H = H ′ ⊕ H ′′ . So any X ∈ Λ T can be written as x ⊕ 0 in the decomposition H ′ ⊕ H ′′ .

Now let us show that the map T keeps Tr( H ′ ) ⊕ 0 invariant. Let x ∈ Tr( H ′ ) be positive semi-definite operator, and S = s ⊕ 0. Recall that x and s are both positive semi-definite operator on finite-dimensional Hilbert space and supp( x ) ⊆ supp( ), s thus there exists /epsilon1 &gt; 0 such that s -/epsilon1x &gt; 0 [5, Lemma 1.1]. From Lemma D.1 Λ T is invariant under action of T , so T ( s ⊕ 0) = s ′ ⊕ 0. As T is linear we have

$$\cdot \ A \ \cdot$$

Since T is complete positive, both s ′ ⊕ 0 and T ( x ⊕ 0) are positive semi-definite. Consequently, the support of T ( x ⊕ 0) must be a subset of the support of s ′ ⊕ 0, thus supp( T ( x ⊕ 0)) ⊆ H ′ . It can be easily generalized to any x ∈ Tr( H ′ ) by using the fact that any operator can be written as linear combination of four positive semi-definite operator. Therefore T keeps Tr( H ′ ) invariant.

Let us now define the map ¯ : Tr( T H ′ ) → Tr( H ′ ) as follows: For any x ∈ Tr( H ′ ), let ¯ ( T x ) = T ( x ⊕ 0), where x ⊕ 0 is written in the decomposition H ′ ⊕ H ′′ . Since Tr( H ′ ) is invariant under the action of T , it is easy to verify that ¯ T is a CPTP map on operators on the finite-dimensional Hilbert space H ′ .

Now we can apply Eq. (3) and Theorem 2.3 to ¯ . T By definition, the peripheral subspace of ¯ ( T X ) is also Λ T (up to 0 acting on H ′′ ).

This gives the achievable infinite-time capacity of T . The only remaining task is to clarify how the encoder and decoder on Tr( H ′ ) can be extended to Tr( H ). For an encoder E : Tr( C d ) → Tr( H ′ ), we define the extension to Tr( H ) as ˜ ( E X ) = E ( X ) ⊕ 0. For a decoder D : Tr( H ′ ) → Tr( C d ), we define the extension to Tr( H ) as ˜ ( D X ) = D ( P H ′ XP H ′ ) + (1 -tr( P H ′ X )) | ψ ψ 〉 〈 | where | ψ 〉 is an arbitrary state in H ′′ and P H ′ is the orthogonal projector on H ′ .
Omar Fawzi Mizanur Rahaman Mostafa Taheri
Capacities of quantum Markovian noise for large times
Given a quantum Markovian noise model, we study the maximum dimension of a
classical or quantum system that can be stored for arbitrarily large time. We
show that, unlike the fixed time setting, in the limit of infinite time, the
classical and quantum capacities are characterized by efficiently computable
properties of the peripheral spectrum of the quantum channel. In addition, the
capacities are additive under tensor product, which implies in the language of
Shannon theory that the one-shot and the asymptotic i.i.d. capacities are the
same. We also provide an improved algorithm for computing the structure of the
peripheral subspace of a quantum channel, which might be of independent
interest.