Communication Enhancement Through Quantum Coherent Control of $N$ Channels in an Indefinite Causal-order Scenario

TL;DR

This paper explores how quantum coherent control of multiple channels in indefinite causal order can significantly enhance information transmission, surpassing classical limits, with explicit formulas derived for small networks.

Contribution

It extends quantum control to networks of N channels with arbitrary capacities, providing a formalism and explicit results for the quantum switch output and Holevo limit.

Findings

Transmission for three channels is twice that of two channels with full superposition of causal orders.

Explicit formulas derived for N=2 and N=3 channels.

Quantum indefinite causal order enhances information capacity beyond classical limits.

Abstract

In quantum Shannon theory, transmission of information is enhanced by quantum features. Up to very recently, the trajectories of transmission remained fully classical. Recently, a new paradigm was proposed by playing quantum tricks on two completely depolarizing quantum channels i.e. using coherent control in space or time of the two quantum channels. We extend here this control to the transmission of information through a network of an arbitrary number $N$ of channels with arbitrary individual capacity i.e. information preservation characteristics in the case of indefinite causal order. We propose a formalism to assess information transmission in the most general case of $N$ channels in an indefinite causal order scenario yielding the output of such transmission. Then we explicitly derive the quantum switch output and the associated Holevo limit of the information transmission for $N=2$, $N=3$ as a function of all involved parameters. We find in the case $N=3$ that the transmission of information for three channels is twice of transmission of the two channel case when a full superposition of all possible causal orders is used.

Figures

Tables

Table 1. Table 1: Values of the Holevo information ratio χ Q ​ 3 ​ S / χ Q ​ 2 ​ S subscript 𝜒 Q 3 𝑆 subscript 𝜒 Q 2 𝑆 \chi_{\text{Q}3S}/\chi_{\text{Q}2S} . The mean value of the ratio is 1.9328 ± plus-or-minus \pm 0.0617.

$d$	${\chi }_{\text{Q}2S}$	${\chi }_{\text{Q}3S}$	${\chi }_{\text{Q}2S}/{\chi }_{\text{Q}3S}$
2	0.0487	0.0980	2.0123
3	0.0183	0.0339	1.8524
4	0.0085	0.0159	1.8705
5	0.0046	0.0087	1.8913
6	0.0027	0.0053	1.9629
7	0.0018	0.0034	1.8888
8	0.0012	0.0023	1.9166
9	0.0008	0.0016	2
10	0.0006	0.0012	2

Equations

$${\mathcal{N}}_{\pi }:=\pi ({\mathcal{N}}_1\circ \cdots \circ {\mathcal{N}}_N)$$

$${\mathcal{N}}_{q_j}^D(\rho )$$

$$=\frac{1-q_j}{d^2}\sum _{i_j=0}^{d^2}{U}_{i_j}^j\rho U_{i_j}^{j†}$$

$$\mathcal{S}({\mathcal{N}}_1,{\mathcal{N}}_2,\dots ,{\mathcal{N}}_N)\left(\rho \otimes {\rho }_c\right)=\sum _{\{i_j\}{|}_{j=1}^N}{W}_{\mathbf{i}}(\rho \otimes {\rho }_c)W_{\mathbf{i}}^{†}$$

$$\sum _{\{i_j\}{|}_{j=1}^N}\to \sum _{z=0}^N\sum _{A_z\in {\mathbf{A}}_z^N}\sum _{b\in B_z},$$

$$\mathcal{S}({\mathcal{N}}_1,{\mathcal{N}}_2,\dots ,{\mathcal{N}}_N)\left(\rho \otimes {\rho }_c\right)=\sum _{z=0}^N{\mathcal{S}}_z,$$

$${\mathcal{S}}_z=\sum _{k,k^{\prime }=1}^{N!}\sqrt{P_k{P}_{k^{\prime }}}\sum _{A_z\in {\mathbf{A}}_z^N}{f}_{A_z}\cdot Q_{A_z}^{k,k^{\prime }}\otimes |k⟩⟨k^{\prime }|$$

$$f_{A_z}=d^{2(z-N)}\prod _{j=1}^N(1-q_j)\prod _{a\in A_z}\frac{q_a}{1-q_a}$$

$$Q_{A_z}^{k,k^{\prime }}=\sum _{\{i_b|b\in B_z\}}{\pi }_k\left(U_{i_1}\cdots U_{i_N}\right)\rho {\left[{\pi }_{k^{\prime }}\left(U_{i_1}\cdots U_{i_N}\right)\right]}^{†}.$$

$$\mathcal{S}({\mathcal{N}}_1,{\mathcal{N}}_2)(\rho \otimes {\rho }_c)={\mathcal{S}}_0+{\mathcal{S}}_1+{\mathcal{S}}_2.$$

$$\displaystyle\begin{array}[]{ll}Q^{1,1}_{\emptyset}=\sum_{i_{1},i_{2}}\pi_{1}(U_{i_{1}}U_{i_{2}})\rho\pi_{1}(U_{i_{1}}U_{i_{2}})^{\dagger}\\ \hskip 21.33955pt=\sum_{i_{1},i_{2}}(U_{i_{1}}U_{i_{2}})\rho(U_{i_{2}}^{\dagger}U_{i_{1}}^{\dagger})\\ \hskip 21.33955pt=d\sum_{i_{1},i_{2}}U_{i_{1}}U_{i_{1}}^{\dagger}=d^{3}\mathbb{I}.\\ Q^{1,2}_{\emptyset}=\sum_{i_{1},i_{2}}\pi_{1}(U_{i_{1}}U_{i_{2}})\rho\pi_{2}(U_{i_{1}}U_{i_{2}})^{\dagger}\\ \hskip 21.33955pt=\sum_{i_{1},i_{2}}(U_{i_{1}}U_{i_{2}})\rho(U_{i_{1}}^{\dagger}U_{i_{2}}^{\dagger})\\ \hskip 21.33955pt=d\sum_{i_{1},}U_{i_{1}}{\rm tr}(\rho U_{i_{1}}^{\dagger})=d^{2}\rho.\end{array}$$

$${\mathcal{S}}_0=\sum _{(\alpha ,{\alpha }^{\prime })\in \mathfrak{A}}\frac{r_0\mathbb{I}}{d}\sqrt{P_{\alpha }{P}_{{\alpha }^{\prime }}}\otimes |\alpha ⟩⟨{\alpha }^{\prime }|+\sum _{(\beta ,{\beta }^{\prime })\in \mathfrak{B}}\frac{r_0\rho }{d^2}\sqrt{P_{\beta }{P}_{{\beta }^{\prime }}}\otimes |\beta ⟩⟨{\beta }^{\prime }|,$$

$${\mathcal{S}}_1=\sum _{k,k^{\prime }}\frac{r_1}{d}\sqrt{P_k{P}_{k^{\prime }}}\mathbb{I}\otimes |k⟩⟨k^{\prime }|=\frac{r_1}{d}\mathbb{I}\otimes {\rho }_c,$$

$${\mathcal{S}}_2=\sum _{k,k^{\prime }}{r}_2\rho \sqrt{P_k{P}_{k^{\prime }}}\otimes |k⟩⟨k^{\prime }|=r_2\rho \otimes {\rho }_c,$$

$$\begin{array}[]{lll}{\mathcal{S}}_{0}=\left(\begin{array}[]{cc}\frac{r_{0}}{d}\mathbb{I}P_{1}&\frac{r_{0}\rho}{d^{2}}\sqrt{P_{1}P_{2}}\\ \frac{r_{0}\rho}{d^{2}}\sqrt{P_{2}P_{1}}&\frac{r_{0}}{d}\mathbb{I}P_{2}\end{array}\right),\\[10.00002pt] {\mathcal{S}}_{1}=\left(\begin{array}[]{cc}\frac{r_{1}}{d}\mathbb{I}P_{1}&\frac{r_{1}}{d}\mathbb{I}\sqrt{P_{1}P_{2}}\\ \frac{r_{1}}{d}\mathbb{I}\sqrt{P_{2}P_{1}}&\frac{r_{1}}{d}\mathbb{I}P_{2}\end{array}\right),\\[10.00002pt] {\mathcal{S}}_{2}=\left(\begin{array}[]{cc}r_{2}\rho P_{1}&r_{2}\rho\sqrt{P_{1}P_{2}}\\ r_{2}\rho\sqrt{P_{2}P_{1}}&r_{2}\rho P_{2}\end{array}\right).\end{array}$$

$$\begin{array}[]{ll}{\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2})(\rho\otimes\rho_{c})=\left(\begin{array}[]{cc}a_{1}&b\\ b&a_{2}\end{array}\right),\end{array}$$

$$\mathcal{S}({\mathcal{N}}_1,{\mathcal{N}}_2,{\mathcal{N}}_3)\left(\rho \otimes {\rho }_c\right)={\mathcal{S}}_0+{\mathcal{S}}_1+{\mathcal{S}}_2+{\mathcal{S}}_3.$$

$$Q_{\varnothing }^{1,1}=\sum _{i_1,i_2,i_3}{\pi }_1(U_{i_1}{U}_{i_2}{U}_{i_3})\rho {\pi }_1(U_{i_1}{U}_{i_2}{U}_{i_3}{)}^{†}=d^5\mathbb{I}.$$

$$Q_{\varnothing }^{1,4}=\sum _{i_1,i_2,i_3}{\pi }_1(U_{i_1}{U}_{i_2}{U}_{i_3})\rho {\pi }_4(U_{i_1}{U}_{i_2}{U}_{i_3}{)}^{†}=d^4\rho .$$

$$\begin{array}[]{ll}Q^{1,2}_{\emptyset}=\sum_{i_{1},i_{2},i_{3}}\pi_{1}(U_{i_{1}}U_{i_{2}}U_{i_{3}})\rho\pi_{2}(U_{i_{1}}U_{i_{2}}U_{i_{3}})^{\dagger}\\ \hskip 21.33955pt=\sum_{i_{1},i_{2},i_{3}}(U_{i_{1}}U_{i_{2}}U_{i_{3}})\rho(U_{i_{2}}^{\dagger}U_{i_{3}}^{\dagger}U_{i_{1}}^{\dagger})\\ \hskip 21.33955pt=d\sum_{i_{1},i_{3}}U_{i_{1}}\operatorname{Tr}(U_{i_{3}}\rho)U_{i_{3}}^{\dagger}U_{i_{1}}^{\dagger}=d^{2}\sum_{i_{1}}U_{i_{1}}\rho U_{i_{1}}^{\dagger}\\ \hskip 21.33955pt=d^{3}\mathbb{I},\\[0.0pt] Q^{1,3}_{\emptyset}=\sum_{i_{1},i_{2},i_{3}}\pi_{1}(U_{i_{1}}U_{i_{2}}U_{i_{3}})\rho\pi_{3}(U_{i_{1}}U_{i_{2}}U_{i_{3}})^{\dagger}\\ \hskip 21.33955pt=\sum_{i_{1},i_{2},i_{3}}(U_{i_{1}}U_{i_{2}}U_{i_{3}})\rho(U_{i_{3}}^{\dagger}U_{i_{1}}^{\dagger}U_{i_{2}}^{\dagger})\\ \hskip 19.91684pt=d\sum_{i_{1},i_{2}}U_{i_{1}}U_{i_{2}}\mathbb{I}U_{i_{1}}^{\dagger}U_{i_{2}}^{\dagger}=d^{2}\sum_{i_{1}}\operatorname{Tr}(U_{i_{2}}\mathbb{I})U_{i_{2}}^{\dagger}\\ \hskip 21.33955pt=d^{3}\mathbb{I},\\[0.0pt] Q^{1,5}_{\emptyset}=\sum_{i_{1},i_{2},i_{3}}\pi_{1}(U_{i_{1}}U_{i_{2}}U_{i_{3}})\rho\pi_{5}(U_{i_{1}}U_{i_{2}}U_{i_{3}})^{\dagger}\\ \hskip 21.33955pt=\sum_{i_{1},i_{2},i_{3}}(U_{i_{1}}U_{i_{2}}U_{i_{3}})\rho(U_{i_{2}}^{\dagger}U_{i_{1}}^{\dagger}U_{i_{3}}^{\dagger})\\ \hskip 21.33955pt=d\sum_{i_{1},i_{3}}U_{i_{1}}\operatorname{Tr}(U_{i_{3}}\rho)U_{i_{1}}^{\dagger}U_{i_{3}}^{\dagger}=d^{3}\sum_{i_{3}}\operatorname{Tr}(U_{i_{3}}\rho)U_{i_{3}}^{\dagger}\\ \hskip 21.33955pt=d^{4}\rho,\\[0.0pt] Q^{1,6}_{\emptyset}=\sum_{i_{1},i_{2},i_{3}}\pi_{1}(U_{i_{1}}U_{i_{2}}U_{i_{3}})\rho\pi_{6}(U_{i_{1}}U_{i_{2}}U_{i_{3}})^{\dagger}\\ \hskip 21.33955pt=\sum_{i_{1},i_{2},i_{3}}(U_{i_{1}}U_{i_{2}}U_{i_{3}})\rho(U_{i_{1}}^{\dagger}U_{i_{2}}^{\dagger}U_{i_{3}}^{\dagger})\\ \hskip 21.33955pt=d\sum_{i_{1},i_{3}}U_{i_{1}}\operatorname{Tr}(U_{i_{3}}\rho U_{i_{1}}^{\dagger})U_{i_{3}}^{\dagger}=d^{2}\sum_{i_{1}}U_{i_{1}}\rho U_{i_{1}}^{\dagger}\\ \hskip 21.33955pt=d^{3}\mathbb{I}.\end{array}$$

$$Q_{\varnothing }^{2,6}=\sum _{i_1,i_2,i_3}{\pi }_2(U_{i_1}{U}_{i_2}{U}_{i_3})\rho {\pi }_6(U_{i_1}{U}_{i_2}{U}_{i_3}{)}^{†}=\sum _{i_1,i_2,i_3}(U_{i_1}{U}_{i_3}{U}_{i_2})\rho (U_{i_1}^{†}{U}_{i_2}^{†}{U}_{i_3}^{†}),$$

$$\begin{array}[]{ll}Q^{i,i^{\prime}}_{\emptyset}=d^{3}\mathbb{I},\forall\hskip 2.84544pt(i,i^{\prime})\in{\mathfrak{I}}\equiv\{(1,6),(2,4),(3,5),(4,2),\\ \hskip 99.58464pt(1,2),(2,1),(3,4),(4,3),(5,6),\\ \hskip 99.58464pt(6,5),(5,3),(6,1),(1,3),(2,5),\\ \hskip 99.58464pt(3,1),(4,6),(5,2),(6,4)\},\\ Q^{j,j^{\prime}}_{\emptyset}=d^{4}\rho,\forall\hskip 2.84544pt(j,j^{\prime})\in{\mathfrak{J}}\equiv\{(1,4),(2,6),(3,2),(4,5),\\ \hskip 99.58464pt(5,1),(6,3),(1,5),(2,3),(3,6),\\ \hskip 99.58464pt(4,1),(5,4),(6,2)\},\\ Q^{k,k^{\prime}}_{\emptyset}=d^{5}\mathbb{I},\forall\hskip 2.84544pt(k,k^{\prime})\in{\mathfrak{K}}\equiv\{(1,1),(2,2),(3,3),\\ \hskip 99.58464pt(4,4),(5,5),(6,6)\}.\\ \end{array}$$

$$\begin{array}[]{ll}{\cal S}_{0}=\displaystyle\sum_{(i,i^{\prime})\in{\mathfrak{I}}}\frac{s_{0}}{d^{3}}\mathbb{I}\sqrt{P_{i}P_{i^{\prime}}}\otimes|i\rangle\langle i^{\prime}|+\sum_{(j,j^{\prime})\in{\mathfrak{J}}}\frac{s_{0}\rho}{d^{2}}\sqrt{P_{j}P_{j^{\prime}}}\otimes|j\rangle\langle j^{\prime}|\\[15.00002pt] \hskip 170.71652pt\displaystyle+\sum_{(k,k^{\prime})\in{\mathfrak{K}}}\frac{s_{0}}{d}\mathbb{I}\sqrt{P_{k}P_{k^{\prime}}}\otimes|k\rangle\langle k^{\prime}|,\end{array}$$

$$\begin{array}[]{ll}Q^{\ell,\ell^{\prime}}_{\{1\}}=d^{2}\rho,\forall\hskip 2.84544pt(\ell,\ell^{\prime})\in{\mathfrak{L}}_{1}\equiv\{(2,3),(3,2),(2,4),(4,2),\\ \hskip 113.81102pt(3,5),(5,3),(3,6),(6,3),(4,5),\\ \hskip 113.81102pt(5,4),(4,6),(6,4),(5,1),(1,5),\\ \hskip 113.81102pt(1,2),(2,1),(1,6),(6,1)\},\\ Q^{m,m^{\prime}}_{\{1\}}=d^{3}\mathbb{I},\forall\hskip 2.84544pt(m,m^{\prime})\in{\mathfrak{M}}_{1}\equiv\{(1,1),(2,2),(3,3),(4,4),\\ \hskip 113.81102pt(5,5),(6,6),(1,3),(1,4),(4,1)\\ \hskip 113.81102pt(3,1),(2,5),(5,2),(2,6),(6,2),\\ \hskip 113.81102pt(3,4),(4,3),(5,6),(6,5)\},\\ \end{array}$$

$$\begin{array}[]{ll}Q^{\ell,\ell^{\prime}}_{\{2\}}=d^{2}\rho,\forall\hskip 2.84544pt(\ell,\ell^{\prime})\in{\mathfrak{L}}_{2}\equiv\{(1,4),(1,5),(1,6),(2,4),\\ \hskip 113.81102pt(2,5),(2,6),(3,4),(3,5),(3,6)\\ \hskip 113.81102pt(4,1),(4,2),(4,3),(5,1),(5,2),\\ \hskip 113.81102pt(5,3),(6,1),(6,2),(6,3)\},\\ Q^{m,m^{\prime}}_{\{2\}}=d^{3}\mathbb{I},\forall\hskip 2.84544pt(m,m^{\prime})\in{\mathfrak{M}}_{2}\equiv\{(1,1),(1,2),(1,3),(2,1),\\ \hskip 113.81102pt(2,2),(2,3),(3,1),(3,2),(3,3)\\ \hskip 113.81102pt(4,4),(4,5),(4,6),(5,4),(5,5),\\ \hskip 113.81102pt(5,6),(6,4),(6,5),(6,6)\},\\ \end{array}$$

$$\begin{array}[]{ll}Q^{\ell,\ell^{\prime}}_{\{3\}}=d^{2}\rho,\forall\hskip 2.84544pt(\ell,\ell^{\prime})\in{\mathfrak{L}}_{3}\equiv\{(1,3),(1,4),(1,6),(2,3),\\ \hskip 113.81102pt(2,4),(2,6),(3,1),(3,2),(3,5)\\ \hskip 113.81102pt(4,1),(4,2),(4,5),(5,3),(5,4),\\ \hskip 113.81102pt(5,6),(6,1),(6,2),(6,5)\},\\ Q^{m,m^{\prime}}_{\{3\}}=d^{3}\mathbb{I},\forall\hskip 2.84544pt(m,m^{\prime})\in{\mathfrak{M}}_{3}\equiv\{(1,1),(1,2),(1,5),(2,1),\\ \hskip 113.81102pt(2,2),(2,5),(3,3),(3,4),(3,6)\\ \hskip 113.81102pt(4,3),(4,4),(4,6),(5,1),(5,2),\\ \hskip 113.81102pt(5,5),(6,3),(6,4),(6,6)\}.\\ \end{array}$$

$$\begin{array}[]{ll}{\cal S}_{1}=\displaystyle\frac{1}{d^{2}}\sum_{s=1}^{3}\left(t_{s}d\sum_{(\ell,\ell^{\prime})\in{\mathfrak{M}}_{s}}\sqrt{P_{\ell}P_{\ell^{\prime}}}\,\mathbb{I}\otimes|\ell\rangle\langle\ell^{\prime}|\right.\left.+t_{s}\displaystyle\sum_{(m,m^{\prime})\in{\mathfrak{L}}_{s}}\sqrt{P_{m}P_{m^{\prime}}}\,\rho\otimes|m\rangle\langle m^{\prime}|\right),\end{array}$$

$$Q_{\{2,3\}}^{k,k^{\prime }}=\sum _{i_1}{\pi }_k(U_{i_1}\cdot \mathbb{I}\cdot \mathbb{I})\rho {\pi }_{k^{\prime }}(U_{i_1}\cdot \mathbb{I}\cdot \mathbb{I}{)}^{†}=d\mathbb{I},$$

$$\begin{array}[]{ll}{\cal S}_{2}=\frac{p_{1}p_{2}p_{3}}{d^{2}}\sum_{k,k^{\prime}}\sqrt{P_{k}P_{k^{\prime}}}\left(\frac{q_{2}q_{3}}{p_{2}p_{3}}Q^{k,k^{\prime}}_{\{2,3\}}\right.\\[10.00002pt] \hskip 28.45274pt+\left.\frac{q_{1}q_{3}}{p_{1}p_{3}}Q^{k,k^{\prime}}_{\{1,3\}}+\frac{q_{1}q_{2}}{p_{1}p_{2}}Q^{k,k^{\prime}}_{\{1,2\}}\right)\otimes|k\rangle\langle k^{\prime}|=\frac{s_{2}}{d}\mathbb{I}\otimes\rho_{c}\end{array}$$

$${\mathcal{S}}_3=s_3\sum _{k,k^{\prime }}\sqrt{P_k{P}_{k^{\prime }}}\rho \otimes |k⟩⟨k^{\prime }|=s_3\rho \otimes {\rho }_c,$$

$$\begin{array}[]{lll}{\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},{\mathcal{N}}_{3})(\rho\otimes\rho_{c})=\left(\begin{array}[]{cccccc}\mathcal{A}_{1}&\mathcal{B}&\mathcal{C}&\mathcal{D}&\mathcal{E}&\mathcal{F}\\ \mathcal{B}&\mathcal{A}_{2}&\mathcal{G}&\mathcal{H}&\mathcal{I}&\mathcal{J}\\ \mathcal{C}&\mathcal{G}&\mathcal{A}_{3}&\mathcal{K}&\mathcal{L}&\mathcal{M}\\ \mathcal{D}&\mathcal{H}&\mathcal{K}&\mathcal{A}_{4}&\mathcal{N}&\mathcal{P}\\ \mathcal{E}&\mathcal{I}&\mathcal{L}&\mathcal{N}&\mathcal{A}_{5}&\mathcal{Q}\\ \mathcal{F}&\mathcal{J}&\mathcal{M}&\mathcal{P}&\mathcal{Q}&\mathcal{A}_{6}\end{array}\right),\end{array}$$

Full text

I Introduction

In information theory, the main tasks to perform are the transmission, codification, and compression of information shannon1948mathematical . Incorporating quantum phenomena, such as quantum superposition and quantum entanglement, into classical information theory gives rise to a new paradigm known as quantum Shannon theory nielsen2002quantum . In this paradigm, each figure of merit can be enhanced: the capacity to transmit information in a channel is increased holevo1998capacity , the security to share a message is improved bennett2014quantum and the storing and compressing of information is optimized schumacher1995quantum . In all these enhancements, only the carriers and the channels of information are considered as quantum entities. On the other hand, connections between channels are still classical, that is, quantum channels are connected setting a definite causal order in space or time. However, principles of quantum mechanics and specifically the quantum superposition principle can be applied to the connections of channels Chiribella2013 , i.e. the trajectories either in space abbott2018communication or time chiribella2018indefinite .

Recently, it has been theoretically ebler2018enhanced and experimentally goswami2018communicating ; guo2018experimental shown that two completely depolarizing channels can surprisingly transmit classical information when combined under an indefinite causal order (i.e., when the order of application of the two channels is not one after another instead of a quantum superposition of the two possibilities). In this paper, we tackle the general situation of an arbitrary number $N$ of channels with arbitrary parameters associated to the control and depolarizing strength. As $N$ is greater than two, the number of different causal orders increases as $N!$

The indefiniteness of causal order has been recently theoretically proposed as a novel resource for applications to quantum information theory chiribella2012perfect ; Araujo2014 and quantum communication salek2018quantum ; guerin2016exponential . Initially, indefinite causal orders have been studied and implemented using two parties with the proposal of a quantum switch by Chiribella et al. Chiribella2013 followed by experimental demonstrations guo2018experimental ; procopio2015experimental ; goswami2018indefinite ; wei2018experimental ; rubino2017experimental The quantum switch is an example of quantum control where a switch can, like its classical counterpart, routes a target system to undergo through two operators in series following one causal order ($1$ then $2$) or the other ($2$ then $1$). But this quantum switch can also trigger a whole new quantum trajectory where the ordering of the two operators is indefinite. Efforts to describe the quantum switch in a multipartite scenario of more than two quantum operations have recently started wechs2018definition ; oreshkov2016causal with an application to reduce the number of queries for quantum computationAraujo2014 .

Specifically, in a quantum $N$-switch used in a second-quantized Shannon theory context chiribella2018second , the order of application of $N$ channels ${\mathcal{N}}_{j}$ on a target system $\rho$ is coherently controlled by a control system $\rho_{c}$. The state of $\rho_{c}$ encodes for the temporal combination of the $N$ channels applied to $\rho$. There are $N!$ different possibilities of definite causal orders using each channel once and only once, as sketched in Figures 1 and 2 for $N=2$ and $N=3$ respectively. In those figures, when the wiring passes through the channel, there is a single channel use, i.e. the target system passes once through one physical channel procopio2015experimental . We discard all wirings with multiple use of the same channel and missing channels abbott2018communication . For each causal order of channels, the overall operator is

[TABLE]

where $\pi$ is a permutation element of the symmetric group $S_{N}=\{\pi_{k}|k\in\{1,2,\ldots,N!\}\}$, and $k$ is associated to a specific definite causal order (equivalent to a single element of $S_{N}$) to combine the $N$ channels where each channel is used once and only once.

In a quantum $N$-switch, the control state $\rho_{c}$ in the state $\left|1\right>\left<1\right|$ for instance fixes the order of application of the channels to be ${\mathcal{N}}_{\text{Id}}={\mathcal{N}}_{1}\circ{\mathcal{N}}_{2}\circ\cdots\circ{\mathcal{N}}_{N}$. Whereas, choosing $\rho_{c}=\left|k\right>\left<k\right|$, $k\leq N!$ would assign another ordering ${\mathcal{N}}_{\pi_{k}}={\mathcal{N}}_{\pi_{k}(1)}\circ{\mathcal{N}}_{\pi_{k}(2)}\circ\cdots\circ{\mathcal{N}}_{\pi_{k}(N)}$ (defined by the effect of the permutation element $\pi_{k}\in S_{N}$ on the order of channels). The key to accessing indefinite causal order of the channels is thus to put $\rho_{c}$ in a superposition of the $\left|k\right>\left<k\right|$ states (e.g. $\rho_{c}=\left|+\right>\left<+\right|$ where $\left|+\right>\equiv\frac{1}{\sqrt{N!}}\sum\left|k\right>$).

The paper is organized as follows. Section II is devoted to the general theoretical framework for the investigation of the transmission of classical information over $N$ noisy channels with arbitrary degree of depolarization, i.e. arbitrary level of noise. Section II also gives the channels representation in terms of Kraus operators performed from those operators for a single depolarizing channel. In Section III, following the previous formalism, we explicitly analyze the case $N=2$, generalizing the outcomes in the literature ebler2018enhanced to any degree of depolarization and level of coherent control. Similarly, the case $N=3$ is developed in the same section. Finally, conclusions and perspectives are given in Section V.

II Transmission over multiple channels in quantum superposition of causal order

In the current development, the sender prepares the target system in the state $\rho$, where the information to transmit is encoded. A control system $\rho_{c}$ is associated to the target system to coherently control the causal order for the application of $N$ quantum communication channels. We relate the basis for the quantum state $\rho_{c}$ mapping their elements on those of the symmetric group of permutations $S_{N}:\frac{1}{N!}\sum_{k,k^{\prime}}|k\rangle\langle k^{\prime}|$. Then, the sender introduces as input $\rho\otimes\rho_{c}$ to a network of $N$ partially depolarizing channels ${\mathcal{N}}_{i}=\mathcal{N}_{q_{i}}^{D},1\leq i\leq N$ applied in series (i.e. the output of one channel becomes the input of the next channel). Throughout this work the $N$ depolarizing channels ${\mathcal{N}}_{1}$, ${\mathcal{N}}_{2},\ldots$, ${\mathcal{N}}_{N}$ can have different depolarization strengths $1-q_{j}$, (thus, ${\mathcal{N}}_{q_{j}}^{D}$ is sometimes used for ${\mathcal{N}}_{j}$ to improve the readability).

After the network, the receiver gets the output state ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},\ldots,{\mathcal{N}}_{N})(\rho\otimes\rho_{c})$, where ${\cal S}$ is the quantum $N$-switch channel. No information is encoded by the sender into the control system controlling the way information is transmitted. Eventually, the receiver retrieves the information decoded in $\rho$.

Communication quantum channels in a network are mathematically described with completely positive trace preserving maps (CPTP). Here, we adopt the Kraus decomposition nielsen2002quantum ${\mathcal{N}}_{(}\rho)=\sum_{i}K_{i}\rho K_{i}^{\dagger}$ to describe the action of a total depolarizing channel ${\mathcal{N}}$ on the quantum state $\rho$ ($i\in\{1,2,\ldots,d\}$): ${\mathcal{N}}_{\rho}={\operatorname{Tr}}[\rho]\frac{\mathbb{I}}{d}$. The set of $d^{2}$ non-unique and generally non-unitary Kraus operators $\{K_{i}\}$ satisfies the completeness condition $\sum_{i=1}^{d^{2}}K_{i}K_{i}^{\dagger}=\mathbb{I}$. Thus, to describe the action of the $j$-th partially depolarizing channel ${\mathcal{N}}_{j}$ on a $d$-dimensional quantum system $\rho$, we write as in ebler2018enhanced

[TABLE]

where each ${\mathcal{N}}_{j}={\mathcal{N}}_{q_{j}}^{D}$ is thus decomposed on an orthonormal basis $\{U_{i_{j}}^{j}\}|_{i_{j}=1}^{d^{2}}$. Then, we define $K_{i_{j}}^{j}=\frac{\sqrt{1-q_{j}}}{d}U_{i_{j}}^{j}$ for $i_{j}\neq 0$, where the added non-unitary operator $U_{0}^{j}=\frac{d\sqrt{q_{j}}}{\sqrt{1-q_{j}}}\mathbb{I}_{t}$, for $i_{j}=0$. Besides ${\mathcal{N}}_{j}$ has no noise when $q_{j}=1$. On the other hand, ${\mathcal{N}}_{j}$ is completely depolarizing when $q_{j}=0$. The results reported in ebler2018enhanced ; abbott2018communication are mainly related to two completely depolarizing ($q_{1}=q_{2}=0$) channels ${\mathcal{N}}_{1}$ and ${\mathcal{N}}_{2}$, despite the generalization is outlined. Below we extend the results from ebler2018enhanced to the case of a quantum switch with $N$ channels ${\mathcal{N}}_{j}$ with arbitrary individual depolarization strengths $q_{j}$.

II.1 The formalism for a quantum $N$-switch channel ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},\ldots,{\mathcal{N}}_{N})$

We define the control state $\rho_{c}$ as $\rho_{c}=\left|\psi_{c}\right>\left<\psi_{c}\right|=\sum_{k,k^{\prime}=1}^{N!}\sqrt{P_{k}P_{k^{\prime}}}\left|k\right>\left<k^{\prime}\right|$ where $P_{k}$ is the probability to apply the causal order $k$ (corresponding to the permutation $\pi_{k}$ as it was previously stated) to the channels such that $\sum_{k=1}^{N!}P_{k}=1$.

The action of the quantum $N$-switch channel ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},\ldots,{\mathcal{N}}_{N})$ can be expressed through generalized Kraus operators $W_{i_{1}i_{2}\ldots i_{N}}$ for the full quantum channel resulting from the switching of $N$ channels as

[TABLE]

where $W_{\mathbf{i}}:=W_{i_{1}i_{2}\ldots i_{N}}=\sum_{k=1}^{N!}K_{\pi_{k}}\otimes\left|k\right>\left<k\right|$ and $K_{\pi_{k}}$ has been defined similarly to equation (1) : $K_{\pi_{k}}:=\pi_{k}(K_{i_{1}}^{1}\cdots K_{i_{N}}^{N})$ where $\pi_{k}$ acts on the index $j$, and the sum over $\{i_{j}\}|_{j=1}^{N}$ means all $i_{j}$ associated to each channel ${\mathcal{N}}_{j}$ vary from [math] to $d^{2}$. We verify (see Appendix A) that these generalized Kraus operators satisfy the completeness property $\sum_{\{i_{j}\}|_{j=1}^{N}}W_{\mathbf{i}}W_{\mathbf{i}}^{\dagger}=\mathbb{I}_{t}\otimes\mathbb{I}_{c}$, where identity operators in the target and control systems spaces are denoted $\mathbb{I}_{t}$ and $\mathbb{I}_{c}$, respectively. This check of completeness suggests how the $i_{j}$ indices allow the systematic reordering of the sums by isolating and grouping the $i_{j}=0$ cases. To distinguish those terms, we introduce the number $z$ of indices $i_{j}$ equal to zero. The sums over the indices $i_{j}$ can then be rearranged as

[TABLE]

where $A_{z}$ is the set of $z$ indices equal to zero ($i_{a}=0$, $\forall a\in A_{z}$) and $B_{z}$ is the complementary set of indices in $\{1,2\ldots,N\}$ : $i_{b}\neq 0$, $i_{b}\in\{1,2,\ldots,d^{2}\}$ for all $b\in B_{z}$. Then, $U_{\pi_{k}}U_{\pi_{k}}^{\dagger}=d^{2z}h_{A_{z}}\mathbb{I}_{t}$, where $h_{A_{z}}=\prod_{a\in A_{z}}\frac{q_{a}}{1-q_{a}}$ and $h_{A_{0}}=1$.

Introducing the Kraus operators $W_{\mathbf{i}}$ into $\mathcal{S}({\mathcal{N}}_{1},{\mathcal{N}}_{2},\ldots,{\mathcal{N}}_{N})$, equation (3) can be written as a sum of $N+1$ matrices ${\cal S}_{z}$ whose $N!\times N!$ elements are matrices of dimension $d\times d$ involving exactly $z$ factors $U_{i_{j}}$ equal to the identity operator. The overall dimension ${\cal S}_{z}$ is thus $dN!\times dN!$

[TABLE]

and (see Appendices A and B)

[TABLE]

with

[TABLE]

where ${\bf A}_{z}^{N}$ is the collection of all possible subsets ${A}_{z}$ of $z$ subscripts in $\{1,2,\ldots,N\}$ corresponding to the $z$ indices equal to zero (i.e. $i_{a}=0$ $\forall a\in A_{z}$). The following subsections detail examples with $N=2$ and $N=3$. The coefficients $Q^{k,k^{\prime}}_{A_{z}}$ are given by

[TABLE]

The $U_{i_{j}}^{j}$ of equation (II) have been simplified in $U_{i_{j}}$. We can see from equation (7) that the elements of the matrix ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},\ldots,{\mathcal{N}}_{N})\left(\rho\otimes\rho_{c}\right)$ will always be linear combination of $\rho$ and $\mathbb{I}_{1}$, whatever $N$ channels. Note also the operators $U_{i_{j}}$ for ${j\in A_{z}}$ are identity operators $\mathbb{I}_{t}$ by construction. Thus, arguments $U_{i_{1}}\cdots U_{i_{N}}$ under $\pi_{k}$ in (7) involves $N$ elements, $z$ of them in $A_{z}$ and $N-z$ in $B_{z}$. The matrix $\mathcal{S}$ and the pivotal equations (5-7) contain all information about the correlations between precise causal orders coherently controlled by $\rho_{c}$ and the output of the quantum switch. $\mathcal{S}$ is a function of several parameters: the involved causal orders $\pi_{k}$ via the probabilities $P_{k}$, the depolarization strengths $q_{i}$’s of each individual channel ${\mathcal{N}}_{i}$, the dimension $d$ of the target system undergoing the operations of those channels and the number of channels $N$. Notably the sum over $k$ and $k^{\prime}$ in equation (6) can be restricted to a subset of definite causal orders via the probabilities $P_{k}$, i.e. a subset of superposition of $m$ causal orders among the $N!$ existing ones for advanced quantum control. This handle had remained unexplored up to now. It was not accessible to former explorations limited to two channels. In the current work we consider only superpositions of all causal orders. The control of causal orders will be presented elsewhere.

In the following subsections we will give the explicit expressions of the quantum switch matrices for the quantum $N$-switch channel for $N=2$ and $N=3$. We access these matrices of the quantum $N$-switch channel via the systematic ordering of the terms in equations (3) as settled in equations (5-7).

The explicit calculation of the quantum $N$-switch channel gives important insights on the transmission of information coherently controlled by $\rho_{c}$ in a fascinating multi-parameter space. We briefly review below some of the intriguing behaviors associated to the parameters exploration in the $N=2$ and the $N=3$ cases. We show indeed in those cases how the nature and number of the causal orders in the control state superposition, the dimension of the target system, the level of noise all play a role. We underline that the $N=3$ case is still untouched experimentally.

To derive equation (6) for particular cases of $N$, we first introduce the definitions of $W_{\mathbf{i}}$ and $\rho_{c}$ into equation (3). Introducing the definitions of the Kraus operators in terms of $U_{i_{j}}^{j}$ operators and applying the same reordering on the sums as in equation (4) leads to equation (6). In the following subsection, specific developments for $N=2$ and $N=3$ to evaluate the $Q^{k,k^{\prime}}_{A_{z}}$ are given simplifying $Q^{k,k^{\prime}}_{A_{z}}$ in (7) by following the relations presented in Appendix B (equations (B.1)-(B.3)).

III The quantum switch matrices for $N=2$ and $N=3$

To show the usefulness of equations (6), we derive general expressions to investigate the transmission of information through two and three channels in an indefinite causal order. Our method can be easily applied to any number of depolarizing channels provided that $\{U_{i}\}_{i=1}^{d^{2}}$ are unitary operators setting an orthonormal basis for the space of $d\times d$ matrices.

III.1 Evaluation of $\mathcal{S}$ for $N=2$

To explicitly evaluate equation (5) with two channels, we identify the two permutations in $S_{2}$ : $\pi_{1}=\left(\begin{smallmatrix}1&2\\ 1&2\end{smallmatrix}\right)$ and $\pi_{2}=\left(\begin{smallmatrix}1&2\\ 2&1\end{smallmatrix}\right)$. Equation (5) for the quantum 2-switch channel matrix acting on the input state $\rho\otimes\rho_{c}$ writes

[TABLE]

The collection of all subsets of subscripts in $\{1,2\}$ are ${\bf A}_{0}^{2}=\{\emptyset\},{\bf A}_{1}^{2}=\{\{1\},\{2\}\}$ and ${\bf A}_{2}^{2}=\{\{1,2\}\}$. Then, the corresponding complementary collections are

${\bf B}_{0}^{2}=\{\{1,2\}\},{\bf B}_{1}^{2}=\{\{2\},\{1\}\}$ and ${\bf B}_{2}^{2}=\{\emptyset\}$.

Coefficients for $\mathcal{S}_{0}$. In this case, we use ${\bf A}_{0}^{2}=\{\emptyset\}$ to calculate the coefficients $Q^{k,k^{\prime}}_{\emptyset}$, $k,k^{\prime}\in\{1,2\}$. The $Q^{k,k^{\prime}}_{\emptyset}$ then reads

[TABLE]

where we have used equations (B.1) and (B.3) for $Q^{1,1}_{\emptyset}$, equation (B.1) with $X=U_{i_{2}}\rho$ and equation (B.2) for $Q^{1,2}_{\emptyset}$. Likewise, we have $Q^{\alpha,\alpha^{\prime}}_{\emptyset}=d^{3}\mathbb{I},\quad\text{for}\quad(\alpha,\alpha^{\prime})\in{\mathfrak{A}}\equiv\{(1,1),(2,2)\}$ and $Q^{\beta,\beta^{\prime}}_{\emptyset}=d^{2}\rho,\quad\text{for}\quad(\beta,\beta^{\prime})\in{\mathfrak{B}}\equiv\{(1,2),(2,1)\}$. Then, we may write

[TABLE]

where $r_{0}=p_{1}p_{2}$ with $p_{i}=1-q_{i}$.

Coefficients for $\mathcal{S}_{1}$. In this case ${\bf A}_{1}^{2}=\{\{1\},\{2\}\}$ and ${\bf B}_{1}^{2}=\{\{2\},\{1\}\}$. Let us first consider the coefficient $Q^{\gamma,\gamma^{\prime}}_{\{1\}}=\sum_{i_{2}}\pi_{\gamma}(\mathbb{I}\cdot U_{i_{2}})\rho\pi_{\gamma^{\prime}}(\mathbb{I}\cdot U_{i_{2}})^{\dagger}=d\mathbb{I}$, using the general relations (B.1)-(B.3), for $(\gamma,\gamma^{\prime})\in{\mathfrak{G}}\equiv\{(1,1),(1,2),(2,1),(2,2)\}$. Since indices are dumb it can be shown that $Q^{\gamma,\gamma^{\prime}}_{\{2\}}=Q^{\gamma,\gamma^{\prime}}_{\{1\}}$ for all $(\gamma,\gamma^{\prime})$. Then the term ${\cal S}_{1}$ can be written as

[TABLE]

where $r_{1}=q_{1}p_{2}+q_{2}p_{1}.$

*Coefficients for $\mathcal{S}_{2}$. *Finally, let us consider the term ${\cal S}_{2}$. In this case ${\bf A}_{2}^{2}=\{\{1,2\}\}$ and hence ${{\bf B}_{2}^{2}}=\{\emptyset\}$. Note that $Q^{k,k^{\prime}}_{\{1,2\}}=\rho$ for all $k$ and $k^{\prime}$. Thus, the term with $z=2$ reads

[TABLE]

with $r_{2}=q_{1}q_{2}$. By expanding the matrices $\mathcal{S}_{0}$, $\mathcal{S}_{1}$ and $\mathcal{S}_{2}$ in the control qubit basis, $\{\left|1\right>,\left|2\right>\}$, we are able to write

[TABLE]

where $\mathbb{I}=\mathbb{I}_{t}$. Summing those matrices according to equation (5), we find that the quantum 2-switch channel matrix ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2})$ has diagonal elements $a_{k}=P_{k}[(r_{0}+r_{1})\mathbb{I}/d+r_{2}\rho],$ for $k=1,2$ and off-diagonal elements $b=\sqrt{P_{1}P_{2}}[(r_{0}+d^{2}r_{2})\rho/d^{2}+\frac{r_{1}}{d}\mathbb{I}]$, with $r_{0}=p_{1}p_{2}$, $r_{1}=q_{1}p_{2}+q_{2}p_{1}$ and $r_{2}=q_{1}q_{2}$. Thus,

[TABLE]

note that the diagonal and off-diagonal elements $a_{1}$, $a_{2}$ and $b$ are matrices and are linear combinations of matrices $\rho$ and $\mathbb{I}_{t}$. This property is non-unique for case $N=2$, instead is general for $N$ channels, an advisable aspect from equation (7) and equations (B.1)-(B.2). Indeed, (20) gives as particular outputs the predicted Holevo capacity of Figure 3 in goswami2018communicating and expressions of Holevo information in ebler2018enhanced .

We end up this subsection stressing that Figure. 1 sketches different ways to connect channels ${\mathcal{N}}_{1}$ and ${\mathcal{N}}_{2}$ in either (a) and (b) a definite causal order and (c) for an indefinite causal order combining the $2!$ possible orders.

III.2 Evaluation of $\mathcal{S}$ for $N=3$

In this section, we explicitly evaluate expression (5) considering three channels. Let us label the 6 elements of $S_{3}$ according to the following set of permutations $\pi_{1}=\left(\begin{smallmatrix}1&2&3\\ 1&2&3\end{smallmatrix}\right)$, $\pi_{2}=\left(\begin{smallmatrix}1&2&3\\ 1&3&2\end{smallmatrix}\right)$, $\pi_{3}=\left(\begin{smallmatrix}1&2&3\\ 2&1&3\end{smallmatrix}\right)$, $\pi_{4}=\left(\begin{smallmatrix}1&2&3\\ 2&3&1\end{smallmatrix}\right)$, $\pi_{5}=\left(\begin{smallmatrix}1&2&3\\ 3&1&2\end{smallmatrix}\right)$ and $\pi_{6}=\left(\begin{smallmatrix}1&2&3\\ 3&2&1\end{smallmatrix}\right)$. Equation (5) for the quantum 3-switch channel matrix acting on input state $\rho\otimes\rho_{c}$ reads

[TABLE]

Coefficients for $\mathcal{S}_{0}$. In this case note that ${\bf A}_{0}^{3}=\{\emptyset\}$, hence ${{\bf B}_{0}^{3}}=\{\{1,2,3\}\}$. Besides, the sum in $Q^{1,k^{\prime}}_{\emptyset}$ is over the indices $\{i_{1},i_{2},i_{3}\}$. These can be computed explicitly

[TABLE]

Likewise,

[TABLE]

The remaining coefficients for ${\cal S}_{0}$ are

[TABLE]

The coefficients $Q^{k,k^{\prime}}_{\emptyset}$ with $k\geq 2$ can be computed using these expressions from equations (24). For instance, consider the following

[TABLE]

which is equivalent to expression $Q^{1,4}_{\emptyset}$ because the indices $i$’s are dumb. Thus one can calculate explicitly the remaining coefficients. Results are thus summarized in the following list

[TABLE]

After calculating all these coefficients, we obtain

[TABLE]

where $s_{0}=p_{1}p_{2}p_{3}$.

Coefficients for $\mathcal{S}_{1}$. In this case ${\bf A}_{1}^{3}=\{\{1\},\{2\},\{3\}\}$ and ${{\bf B}_{1}^{3}}=\{\{2,3\},\{1,3\},\{1,2\}\}$. Let us first consider the coefficient $Q^{k,k^{\prime}}_{\{1\}}$, so that sum must be accomplished over the indices $\{i_{2},i_{3}\}$, hence $Q^{k,k^{\prime}}_{\{1\}}=\sum_{i_{2},i_{3}}\pi_{k}(\mathbb{I}\cdot U_{i_{2}}\cdot U_{i_{3}})\rho\pi_{k^{\prime}}(\mathbb{I}\cdot U_{i_{2}}\cdot U_{i_{3}})^{\dagger}$. Using the relations (B.1)-(B.3) we obtain

[TABLE]

[TABLE]

[TABLE]

Hence, the matrix ${\cal S}_{1}$ can be computed

[TABLE]

where $t_{1}=p_{2}p_{3}q_{1}$, $t_{2}=p_{1}p_{3}q_{2}$ and $t_{3}=p_{1}p_{2}q_{3}$.

Coefficients for $\mathcal{S}_{2}$. In this case ${\bf A}_{2}^{3}=\{\{1,2\},\{1,3\},\{2,3\}\}$ and hence ${{\bf B}_{2}^{3}}=\{\{3\},\{2\},\{1\}\}$. Let us consider

[TABLE]

where the operators $\mathbb{I}$ have been written for the sake of clarity as the permutations $\pi_{k}$ act on sets of three elements. In a similar way $Q^{k,k^{\prime}}_{\{1,3\}}=Q^{k,k^{\prime}}_{\{1,2\}}=d\mathbb{I}.$ Thus, we obtain

[TABLE]

where $s_{2}=q_{1}q_{2}p_{3}+q_{1}q_{3}p_{2}+q_{2}q_{3}p_{1}$.

Coefficients for $\mathcal{S}_{3}$. Finally, note that $Q^{k,k^{\prime}}_{\{1,2,3\}}=\rho$ for all $k$ and $k^{\prime}$. Thus, the term with $z=3$ reads

[TABLE]

where $s_{3}=q_{1}q_{2}q_{3}$ and using the definition of the control qudit.

For three channels, Figure 2 shows different ways to connect channels ${\mathcal{N}}_{1}$, ${\mathcal{N}}_{2}$ and ${\mathcal{N}}_{3}$ in either (a)-(f) a definite causal order, or (g) in an indefinite causal order taking into account all 3! causal orders. The quantum 3-switch matrix is again calculated with equation (5) (see Appendix C)

[TABLE]

where the diagonal and the off-diagonal elements whose expressions are given in Appendix C are also linear combinations of matrices $\rho$ and $\mathbb{I}_{t}$. From the definition of symmetric matrices horn1990matrix , we can see that the quantum switch matrices (20) and (34) are block-symmetric matrices with respect to the main diagonal. This could be seen as general from the fact $Q^{k,k^{\prime}}_{A_{z}}=Q^{k^{\prime},k}_{A_{z}}$ due to equations (6) and (7), because indices in the sums are dumb. Thus, as the number of channels increases, the number of different $d\times d$ matrices involved in the quantum $N$-switch matrix $\mathcal{S}$ scales as $N!(N!+1)/2$. Notice that those matrices also characterize information transmission of any definite causal ordering $\pi_{k}$ of channels $\mathcal{N}_{\pi_{k}}$ when setting $P_{k}=1$ and $P_{s}=0$ for all $s\neq k$.

Matrices in equation (20) or (34) are written in the basis of the control system $\rho_{c}$ which maps and weights the chosen causal orders. To know the best rate to communicate classical information with two and three channels, in the following Section we diagonalize matrices (20) and (34) to compute the Holevo information limit $\chi$, which quantifies how much classical information can be transmitted through a channel in a single use. $\chi$ gives a lower bound on the classical capacity holevo1998capacity ; abbott2018communication ; schumacher1997sending .

IV Holevo information limit for two and three channels

We compute the Holevo information limit (Holevo information for shortness in the following) $\chi(\mathcal{S})$ for $N=2$ and $N=3$ channels through a generalization of the mutual information (see for example wilde2013quantum ) and supplementary information of ebler2018enhanced . The Holevo information $\chi(\mathcal{S})$ is found by maximizing mutual information, and it can be shown that maximization over the $\rho$ pure states is sufficient wilde2013quantum . The Holevo information is then given by

[TABLE]

where $d$ is the dimension of the target system $\rho$, $H({\tilde{\rho}}_{c}^{(N)})$ is the von-Neumann entropy of the output control system $\tilde{\rho}_{c}^{(N)}$ for $N$ channels and $H^{\text{min}}(\mathcal{S})$ is the minimum of the entropy at the output of the channel $\mathcal{S}$. The minimization of $H^{\text{min}}(\mathcal{S})\equiv\underset{\rho}{\text{min}}\hskip 4.30554ptH^{\text{min}}(\mathcal{S}(\rho))$ is over all input states $\rho$ going on the channel $\mathcal{S}$ wilde2013quantum . To evaluate equation (35):

The diagonalization and minimization of $H^{\text{min}}({\mathcal{S}})$ is performed on all possible states given by $\rho$. It is done analytically for $N=2$ channels and arbitrary $q_{i}$. For $N=3$ channels we compute the eigenvalues of the full quantum 3-switch matrix $\mathcal{S}({\mathcal{N}}_{1},{\mathcal{N}}_{2},{\mathcal{N}}_{3})\left(\rho\otimes\rho_{c}\right)$ numerically. 2. 2.

${\tilde{\rho}}_{c}^{(N)}$ was analytically calculated following ebler2018enhanced . 3. 3.

We deduce $H({\tilde{\rho}}_{c}^{(N)})$ from the analytical expressions of ${\tilde{\rho}}_{c}^{(N)}$.

IV.1 Holevo information limit for $N=2$ channels

IV.1.1 Calculation of $H^{\text{min}}$

We calculate the minimum output entropy $H^{\text{min}}(\mathcal{S})$ of the channel $\mathcal{S}\equiv{\mathcal{S}}({\mathcal{N}}_{1},\ldots,{\mathcal{N}}_{N})$

[TABLE]

where the minimization is a priori over all input states $\rho$ and $\{\lambda_{{\cal S}(\rho),i}\}_{i=1}^{d}$ are the eigenvalues of ${\cal S}(\rho)$. In fact it is sufficient to minimize over the states wilde2013quantum and the eigenvalues $\{\lambda_{\rho,i}\}_{i=1}^{d}$ sum up to $1$. As $H^{\text{min}}(\mathcal{S}(\rho))$ is concave, the minimization is done as in Ref. ebler2018enhanced : the eigenvalues $\{\lambda_{{\cal S}(\rho),i}\}_{i=1}^{d}$ are taken at the border of the interval $[0,1]^{\times d}$ and as they sum up to one, the minimization is simplified to the cases where all $\lambda$ but one are set to zero and the last one is equal to 1.

In this situation, ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2})\left(\rho\otimes\rho_{c}\right)$ has only four non-zero matrix elements, (see equation (20)), which can be rewritten as $2\times 2$ matrices

[TABLE]

where $a_{0}$ and $b$ are $d\times d$ matrices and linear combinations of $\rho$ and $\mathbb{I}_{t}$:

[TABLE]

with $p\equiv P_{1}$, $q\equiv P_{2}$ are the control probabilities with $p+q=1$.

Using the commutativity of $\rho$ and $\mathbb{I}_{t}$ (so they have the same eigenvectors), we then retrieve analytically $a_{\pm}$, the matrix-eigenvalues of ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},\ldots,{\mathcal{N}}_{N})\left(\rho\otimes\rho_{c}\right)$

[TABLE]

The existence of this last expression is warranted by the positivity of the discriminant bhatia2009positive , considering the positivity of $\rho$ and the structure of $a_{0}$ and $b$, which are linear combinations of $\mathbb{I}$ and $\rho$.

The commutativity properties of $\rho$ and $\mathbb{I}_{t}$ are inherited to $a_{\pm}$. Then, the eigenvalues of ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2})\left(\rho\otimes\rho_{c}\right)$ for two causal orders are the eigenvalues of $a_{\pm}$.Thus, to diagonalize $a_{\pm}$ we just replace $\rho$ by its eigenvalues, labeled as $\lambda_{\rho,i}$, in equation (38), which generalizes the procedure obtained in ebler2018enhanced . Our procedure gives access to the transmission of information in a more general situation, where the depolarization strengths $q_{i}$ can be different for each channel and it can take any value between 0 and 1. Equation (38) gives the eigenvalues of the matrix ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2})(\rho\otimes\rho_{c})$ ($s=\pm 1$)

[TABLE]

The eigenvalues of $\lambda_{s,i}$ are well defined because of the positivity of discriminant bhatia2009positive . Finally, using the concavity of the entropy, the minimum of the entropy $H^{\rm min}$ for a state is reached by setting just one $\lambda_{\rho,i}$ to one and all the others to zero, with this we obtain

[TABLE]

[TABLE]

[TABLE]

[TABLE]

It is easy to show $\beta_{k}\leq{\alpha_{0,k}}$, then $\lambda_{\pm,i}\geq 0$ $(\lambda_{s,k}\geq 0)$ as expected. Also, $0\leq\lambda_{s,k}\leq 1$ and then $-H^{\rm min}(\mathcal{S}({\mathcal{N}}_{1},{\mathcal{N}}_{2}))\leq 0$.

If one of $q_{1}=1$ , i.e. channel $1$ is free of depolarization, then $\alpha_{0,k}=\frac{p_{2}}{d}+kq_{2}=\beta_{k}$ and

[TABLE]

depends only on the probability of depolarization for channel $2$. Thus, $-H^{\rm min}(\mathcal{S}({\mathcal{N}}_{1},{\mathcal{N}}_{2}))$ reaches its maximum value of zero only if $q_{1}=q_{2}=1$. Alternatively, it is direct to show that the discriminant reaches its maximum value when $p=\frac{1}{2}$, which is the case studied by Ebler et al. ebler2018enhanced . In addition, if $q_{1}=q_{2}=0$, i.e. both channels are fully depolarizing, then $\alpha_{0,k}=\frac{1}{d},\beta_{k}=\frac{k}{d^{2}}$, so $-H^{\rm min}(\mathcal{S}({\mathcal{N}}_{1},{\mathcal{N}}_{2}))$ reaches the minimum value

[TABLE]

For sake of shortness the entropy $H^{\rm min}(\mathcal{S}({\mathcal{N}}_{1},{\mathcal{N}}_{2}))$ will be denoted simply as $H^{\rm min}(\mathcal{S}_{N})$. To illustrate the range of parameters of equation 44, we plot the entropy $H^{\rm min}(\mathcal{S}_{2})$ map for two noisy channels. Figure 3 shows the entropy $H^{\rm min}(\mathcal{S}_{2})$. The plots are contour surfaces of $H^{\rm min}(\mathcal{S}_{2})$ when $q_{1},q_{2}$ vary from 0 to 1. Each plot contains thirty surfaces distributed in their complete range shown in the color-chart. We plot several cases of $H^{\rm min}(\mathcal{S}_{2})$ when the dimension of the target is $d=2,3,10$ and 100.

IV.1.2 Derivation of ${\tilde{\rho}}_{c}^{(2)}$

To obtain the output state of the control system ${\tilde{\rho}}_{c}^{(N)}$ after $N$ channels, we calculate

[TABLE]

where $\omega_{XIJAC}$ is an extended input state with pure conditional state as described in ebler2018enhanced . A direct calculation shows:

[TABLE]

here $\rho^{\prime}=X(i)Z(j)\rho Z(j)^{\dagger}X(i)^{\dagger}$ and $X(i)\left|l\right>=\left|i\oplus l\right>,Z(j)\left|l\right>=e^{2\pi ijl}\left|l\right>$ are the known Heisenberg$-$Weyl operators wilde2013quantum . To isolate the term ${\tilde{\rho}}_{c}^{(N)}$ we apply the following relations

[TABLE]

which are valid for $N\geq 2$ and they are obtained by direct calculation following the former definitions. Then, for $N=2$ we find that the output control state is

[TABLE]

where $p_{i}=1-q_{i}$.

Using the two previous results for $\chi_{{\rm Q}N{\rm S}}$, Figure 4 shows the transmission map of information for two noisy channels. The plots are contour surfaces of $\chi_{\rm Q2S}$ when $q_{1},q_{2}$ and $P_{1}=P_{2}=p$ vary from 0 to 1. The maximum capacity is trivially reached when $q_{1}=q_{2}=1$ simultaneously reaching the value $\chi_{\rm Q2S}=\log d$. The minimum capacity is zero, reached in the boundary of the front sides with $(q_{1}=0,p=0,1)$, $(q_{2}=0,p=0,1)$, $(q_{1}=0,q_{2}=1)$, and $(q_{1}=1,q_{2}=0)$. Notably, for $q_{1}=q_{2}=0$ there are values higher than the minimum. This phenomenon is observed in the protuberance of plots near $\chi_{\rm Q2S}=0$. For larger values of $d$, the protuberance occurs sharply near $q_{1}=0$ and $q_{2}=0$ faces. Note the nearest surface to those faces are for $\chi_{\rm Q2S}=10^{-3},10^{-3},10^{-4},10^{-7}$ respectively for each plot $d=2,3,10,100$.

IV.2 Holevo information for $N=3$ channels

We numerically calculate the eigenvalues of the entropy $H^{\rm min}$ for $N=3$ channels from equation (34). Then, using relations (47) from (27), (31), (32) and (33) we find that the output state is

[TABLE]

As before, putting those outcomes together in (35), Figures 5 (a) and (b) give the Holevo informations $\chi_{\text{Q2S}}$ and $\chi_{\text{Q3S}}$ for two and three channels respectively, as a function of the depolarization strengths $q_{i}$ and the dimension $d$ of the target system. Our model enables us to exhibit a wealth of different behaviors as a function of $d$,$q_{i}$, and $N$ from the fully noisy situation to the identity channel transmission. For the sake of simplicity, we restrict our graphical analysis to equal depolarization strengths, i.e., $q_{1}=q_{2}=q_{3}$, with a balanced superposition of $N!$ causal orders, that is, with equally weighted probabilities $P_{k}=1/N!$ for each case $N=2,3$. The analysis of these results allows us to draw the following conclusions for those particular cases:

• •

For a fixed dimension $d$, the Holevo information for indefinite causal order is always higher than the one obtained using one of the definite causal order shown in Fig. 2. This is especially the case for totally depolarized channels i.e. $q_{i}=0,\forall i$. For completely clean channels ($q=1$), the Holevo information for indefinite and definite causal order converges to the same value depending on $d$ (not shown).

• •

Two regions can be distinguished. In the strongly depolarized region (roughly $q<0.3$ for $N=2$ and $q<0.5$ for $N=3$) the increase of the dimension $d$ of the target system is detrimental to the Holevo information transmitted by the quantum switch. In contrast, in the moderately depolarized region ($q>0.3$ for $N=2$ and $q>0.5$ for $N=3$) the Holevo information increases both with $q$ and $d$, as expected a maximum (not shown) for completely clean channels.

• •

In the strongly depolarized region, increasing the number of channels to $N=3$ is definitively advantageous for information extraction. For instance, in the case of totally depolarized channels ($q=0$), the Holevo information is approximately doubled with $N=3$ with respect to $N=2$ for all values of the dimension $d$ calculated up to $d=10$

In fact, Table 1 gives the values of the ratio $\chi_{\text{Q}3S}/\chi_{\text{Q}2S}$, finding that the Holevo information is approximately doubled for $N=3$ with respect to $N=2$.

V Conclusions

Communication enhancement is a challenging task in quantum information processing due to imperfection of communication channels subjected to depolarization. Causal order has been proposed as a disruptive procedure to improve communication, compression of quantum information, bringing the quantum possibilities into a new frontier. We have analyzed the quantum control of $N$ operators in the context of the second-quantized Shannon theory and in the specific case of superposition of causal orders, extending the results in the current literature. We obtained a general expression for ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},\ldots,{\mathcal{N}}_{N})$for the quantum $N$-switch for an arbitrary number of channels with any depolarizing strength thus providing an operational formula enabling the exploration of communication channels controlled by causal orders. This formula is useful to explore computationally the cases with an increasing $N$. A detailed analysis to assess the information transmission for the cases of $N=2$ and $N=3$ channels is presented: an increasing number of channels improves the transmission of information. In particular, we remarkably found that the Holevo information is doubled when the number of channels goes from $N=2$ to $N=3$.

We give the matrices corresponding to quantum $N$-switches ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},\ldots,{\mathcal{N}}_{N})$ as a function of the number of channels, depolarization strengths, dimension of the target system. We obtain other general properties for the general case of $N$ channels such as the symmetric properties of matrices $Q^{k,k^{\prime}}_{A_{z}}$ and thus of ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},\ldots,{\mathcal{N}}_{N})$. We also demonstrate that ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},\ldots,{\mathcal{N}}_{N})$ is always a linear combination of $\rho$ and $\mathbb{I}_{t}$, whatever $N$ channels. Expressions for the Holevo limit are equally accessible from our expressions and methodology. Besides, we showed that the depolarizing strengths can be used as control parameters to modify the information transmission on demand. We shall develop elsewhere the analysis of control via selected combinations of the $N!$ available causal order enabled by the present work.

Acknowledgements

F. Delgado acknowledges professor Jesús Ramírez-Joachín for the fruitful discussions and teaching in 1983 about combinatorics required in this work. L.M. Procopio wishes to thank Alastair A. Abbott for commenting this manuscript. F. Delgado and M. Enríquez acknowledge the support from CONACyT and from School of Engineering and Science of Tecnológico de Monterrey in the developing of this research work. L.M. Procopio acknowledges the support from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłukodowska-Curie grant agreement No 800306. This work is also supported by a public grant overseen by the French National Research Agency (ANR) as part of the “Investissements d’Avenir” program (Labex NanoSaclay, reference: ANR-10-LABX-0035) and by the Sitqom ANR project (reference : ANR-SITQOM-15-CE24-0005).

Appendix A Completeness property for $W_{\mathbf{i}}$

We demonstrate here the completeness property

[TABLE]

for the generalized Kraus operators $W_{\mathbf{i}}$ for the full quantum $N$-switch channel by relying on the reordering of the sums obtained by grouping terms with indices $i_{s}$ equal to zero. $W_{\mathbf{i}}:=W_{i_{1}i_{2}\ldots i_{N}}=\sum_{k=1}^{N!}K_{\pi_{k}}\otimes\left|k\right>\left<k\right|$ and $K_{\pi_{k}}:=\pi_{k}(K_{i_{1}}^{1}\cdots K_{i_{N}}^{N})$ where $\pi_{k}$ acts on the subscripts $j$ of the Kraus operators $K_{i_{j}}^{j}$. In the sum $\{i_{j}\}|_{j=1}^{N}$, each index in the set of indices $\{i_{1},i_{2},\ldots,i_{N}\}$ is associated to a channel ${\mathcal{N}}_{j}$ where $j\in\{1,2,\ldots,N\}$ and varies from [math] to $d^{2}$. By introducing the definition of the Kraus operators $K_{i_{j}}^{j}=\frac{\sqrt{1-q_{j}}}{d}U_{i_{j}}^{j}$ into $W_{{\mathbf{i}}}$, the left side from Equation (A.1) can be re-written as

[TABLE]

where $U_{\pi_{k}}:=\pi_{k}(U_{i_{1}}^{1}\cdots U_{i_{N}}^{N})$ and $h_{N}=\prod_{j=1}^{N}1-q_{j}$. As $U_{i}U_{i}^{\dagger}=\mathbb{I}$, $\forall i>0$, the product $U_{\pi_{k}}U_{\pi_{k}}^{\dagger}$ reduces almost to the identity $\mathbb{I}_{t}$ except for the factors $U_{0}U_{0}^{\dagger}=\frac{d^{2}q_{j}}{1-q_{j}}\mathbb{I}_{t}$. Using (4) in the above expression, it can be rearranged into:

[TABLE]

where the sum over $A_{z}$ is the sum of terms $h_{A_{z}}$ over all the elements of ${\bf A}_{z}^{N}$, the set of all subset of $z$ elements in $\{1,2,\ldots,N\}$. This yields the factor $f_{A_{z}}=h_{N}h_{A_{z}}d^{2(z-N)}$ in equation (6).

To prove equation (A.1), we then apply the total probability property $\sum_{z=0}^{N}\sum_{a\in A_{z}}h_{A_{z}}=\frac{1}{h_{N}}$ together with the property $\sum_{b\in B_{z}}d^{2(z-N)}=1$. Thus (A.1) is then proved.

Appendix B Relations to evaluate coefficients $Q^{k,k^{\prime}}_{A_{z}}$

We recall below the relations needed to deduce explicitly matrices $\mathcal{S}_{z}$ and then $\mathcal{S}$ for the quantum $N$-switch from the sums and products of the $Q^{k,k^{\prime}}_{A_{z}}$ factors

[TABLE]

where $X$ is any $d\times d$ matrix and $U_{i}$ an orthonormal basis for the $d\times d$ matrices. Applying equation (B.1) to $X=\mathbb{I}$,

[TABLE]

Applying equation (B.1) to $X=\rho$, such that ${\rm Tr}(\rho)=1$, we get a uniform randomization over the set of unitaries $U_{i\neq 0}$ that completely depolarizes the state $\rho$, thus giving the relation $\sum_{i}U_{i}\rho U_{i}^{\dagger}=d\mathbb{I}$.

Appendix C Matrices $\mathcal{S}_{z}$ for the quantum 3-switch

Matrix ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},{\mathcal{N}}_{3})(\rho\otimes\rho_{c})$ is a $6\times 6$ block-matrix, block-symmetric matrix whose matrix elements are matrices of dimension $d\times d$. Since $Q^{k,k^{\prime}}_{A_{z}}=Q^{k^{\prime},k}_{A_{z}}$, for all $A_{z}$, then ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},{\mathcal{N}}_{3})(\rho\otimes\rho_{c})$ is symmetric with respect to the main diagonal. Now, by expanding equations $\mathcal{S}_{0}$, $\mathcal{S}_{1}$, $\mathcal{S}_{2}$ and $\mathcal{S}_{3}$ in the control qudit basis, $\{\left|1\right>,\left|2\right>,\left|3\right>,\left|4\right>,\left|5\right>,\left|6\right>\}$, we can found the quantum 3-switch matrix ${\mathcal{S}}({\mathcal{N}}_{1},{\mathcal{N}}_{2},{\mathcal{N}}_{3})(\rho\otimes\rho_{c})$ having diagonal elements as

[TABLE]

and off-diagonal elements

[TABLE]

Those matrix elements are the entries of the matrix (34) in the main text.