Topological networks for quantum communication between distant qubits

Efficient communication between qubits relies on robust networks, which allow for fast and coherent transfer of quantum information. It seems natural to harvest the remarkable properties of systems characterized by topological invariants to perform this task. Here, we show that a linear network of coupled bosonic degrees of freedom, characterized by topological bands, can be employed for the efficient exchange of quantum information over large distances. Important features of our setup are that it is robust against quenched disorder, all relevant operations can be performed by global variations of parameters, and the time required for communication between distant qubits approaches linear scaling with their distance. We demonstrate that our concept can be extended to an ensemble of qubits embedded in a two-dimensional network to allow for communication between all of them. Qubits are extremely fragile and tend to decohere rapidly; nonetheless they can travel through long chains of resonators unperturbed if the latter are coupled in a very special—topological—way. Researchers from the University of Stuttgart, Germany, contrived a scheme to transfer quantum information reliably between the endpoints of a chain of resonators. Analytical and numerical methods reveal that the transfer becomes resilient to disorder that respects a particular symmetry if the coupled resonators are in a so called topological phase. Such phases are characterized by a winding number (the “topological invariant”) which explains their exceptional stability and, as a consequence, the high fidelity of the transfer. In the future, this method may be employed to connect the separate quantum registers of a quantum processor.


INTRODUCTION
Systems characterized by topological invariants are well known to exhibit unique properties with potential applications in quantum information processing and engineering. 1 Ever since the first experimental observation of the integer quantum Hall effect, 2-4 many other condensed matter systems have been identified and experimentally characterized, such as fractional quantum Hall fluids 5-7 and topological insulators and superconductors. [8][9][10][11][12][13][14][15][16] The latter belong to a particularly well-understood family of topological systems described by non-interacting fermions, where topological invariants can be defined on classes of random matrices. [17][18][19][20] This concept can be straightforwardly generalized to bosonic setups, as well as classical systems, [21][22][23] where the topological features still give rise to intriguing properties, such as localized and chiral edge modes. Here, we are interested in such systems: We demonstrate that their topological properties can be harvested for robust and efficient transfer of quantum information over large distances.
Several different platforms for the realization of topological systems of artificial matter with bosonic degrees of freedom are currently explored: The construction of topological band structures and the observation of edge states has been achieved with photonic circuits in the optical 21,24,25 and the radio-frequency 26 regime, as well as with classical coupled harmonic oscillators, 22,23,27 and with cold atomic gases. [28][29][30][31][32][33][34] These experimental advances have been prepared and are supported by many theoretical proposals, e.g. refs. [35][36][37][38][39][40][41][42][43][44][45]. Several of the above platforms are suitable to carry a single, quantized excitation with low losses and dissipation along protected edge channels, which opens the opportunity to harvest topological phenomena for guiding and transmitting quantum information reliably. First approaches in this direction have been proposed 46,47 and primarily focus on the transmission of excitations along protected edge modes on the boundary of a two-dimensional, topologically non-trivial medium.
Here, we show that a linear network of coupled bosonic degrees of freedom, characterized by topological bands, can be employed for the highly efficient exchange of quantum information over large distances. We demonstrate the superiority of this setup over its topologically trivial counterparts and exemplify its application with the implementation of a robust quantum phase gate. Our proposal is based on a (quasi) one-dimensional setup, characterized by a Z topological index, 48 and derived from paradigmatic systems such as Kitaev's Majorana chain 49 and the Su-Schrieffer-Heeger (SSH) model. 50 It features symmetry protected, localized edge modes, the extend and overlap of which can be tuned via coupling parameters to facilitate controllable communication between them. Important features of our setup are that relevant operations can be performed by global variations of parameters, its robustness with respect to the pulse shapes used for the transfer protocol, and that the time for the transfer scales favorably with the separation of the qubits. This high gain in performance is bought by more complex preparation schemes as the coupling parameters have to respect symmetries protecting the topological invariants. 21,22,26,27 Finally, we demonstrate that our concept can be extended to an ensemble of qubits embedded in a two-dimensional network of local bosonic degrees of freedom to allow for communication between all of them.
We consider macroscopically separated qubits that are coupled by a linear quantum network, see Fig. 1a for an example. The quantum network itself is constructed from bosonic degrees of freedom with only local couplings between them, and generically described by the Hamiltonian Here, b y i (b i ) are bosonic creation (annihilation) operators accounting for the mode at site i with H ij the coupling amplitudes. The network is designed such (see below) that at the end p of each branch, a localized bosonic edge modeb p emerges with a controllable coupling between this mode and a local qubit. The conceptually simplest setup to envisage is an optical network coupled to a single atom with the level structure shown in Fig. 1b. There the coupling Hamiltonian for the target qubit T takes the form (within the rotating wave approximation) The coupling g T (t) is controlled by external laser fields and allows for the application of π-pulses between the qubit state 1 j i T and the edge modeb T , i.e., the emission of a photon into the edge modeb T from state 1 j i T is accompanied by a transition into the auxiliary state a j i T ; in the following, we denote such a π-pulse at edge p by the unitary operation Π p . Note that the Hamiltonian H C ðtÞ for the control qubit C is similar, with the role of 1 j i C and a j i C exchanged.
Several fundamental quantum information processing tasks between the qubits reduce to the transfer of edge excitations within the linear network; we denote the corresponding unitary operation that describes the transfer of excitations between edges p and q as T p↔q . As an example, the protocol for a controlledphase (CP) gate between a control qubit at position C and a target qubit at position T reads Another example, the transport of a control qubit to a target position, is simply described by the protocol Here, A p denotes the exchange of the two states 1 j i p and a j i p . Note that this operation even performs the full exchange of the two qubits due to the linearity of the network; a detailed discussion of these operations can be found at the end of the manuscript.
Motivated by these observations, we are in the following interested in the efficient transfer (T C↔T ) of edge excitations within the linear network. The basic idea is most conveniently illustrated for two qubits coupled by a one-dimensional network as illustrated in Fig. 1a: The structure of the network gives rise to topological bands with a gapped dispersion relation and entails the existence of degenerate and localized edge modes within the bulk gap. As the existence of edge modes is topologically protected, it is robust against disorder. In a finite system, the degeneracy of the topological edge states is only lifted exponentially in the edge separation. However, globally tuning the quantum network closer to the topological phase transition into the trivial phase increases overlap and finite size splitting of the edge states, and eventually allows for a π-pulse T C↔T between the two edges. This simple idea is the core of our protocol; compared to topologically trivial systems, it features several superior properties regarding the speed of state transfer and its robustness against disorder, as well as the absence of individual addressing of each part of the network.
Finally, we would like to stress that our scheme is generic and one can envisage various experimental platforms for its implementation. In addition to the discussed optical network, 21 alternative setups are coupled optical cavities and circuit quantum electrodynamics systems, 51,52 as well as trapped polar molecules or Rydberg atoms with a coupling mediated by dipolar exchange interactions, 45 while the local qubits can be artificial atoms, 53 nitrogen-vacancy centers in diamond, 54 or trapped ions. 55

Topological network
We start with a description of the requirements on the quantum networkĤ n to exhibit topologically protected edge modes in a one-dimensional chain with two edges, as illustrated in Fig. 1a; the generalization to 2D networks is discussed at the end of the manuscript.
The most prominent paradigmatic model in one-dimension is the Majorana chain, which exemplifies the concept of symmetry protected topological phases, originally formulated for spinless fermions with a mean-field p-wave pairing term. 49 This model is closely related to the Su-Schrieffer-Heeger (SSH) model 50 in the single-particle picture. It turns out that the necessary steps to translate these models into our bosonic quantum network  Fig. 1 Topological quantum networks. a One-dimensional chain of bosonic modes b i with globally tunable onsite couplings wðtÞ (orange); derived from the SSH chain. It features a topological band structure with localized edge modes (in the topological phase), which can be coupled to local qubits T and C. b The local qubits are realized as three-level systems with logical states 0 j i p , 1 j i p and the auxiliary state a j i p . The state a j i C 1 j i T À Á can decay into 1 j i C a j i T À Á via an off-resonant transition and thereby emits an excitation into the localized edge modeb C b T . Tunneling excitations between the two edges is facilitated by tuning the chain close to the phase transition via wðtÞ and adiabatically decoupling bulk from edge modes. c Possible 2D generalization of the network on a dimerized honeycomb lattice. Scattering-free transport is guaranteed by topological protection which requires a sublattice symmetry. The latter is realized by directly coupling only "even" (filled circles) with "odd" modes (filled squares). State transfer between qubits (empty circles) of different (the same) type is possible (impossible), illustrated by the bold yellow (gray) path. Stray couplings in the bulk (shown for the upper path) are not detrimental to the transfer fidelity. d Instead of locating the edge modes (filled colored circles and squares) with their qubits at the boundary, emanating SSH chains can be used to separate the qubits from each other and the 2D bulk. There is no need to trace out a specific path as in c, but a weak addressability of the individual chains is sufficient whereas the couplings of the bulk can be tuned globally. Details are given at the end of the manuscript with the two fermion operators c i and c i on each site. Note that the indices i label the "upper" fermionic modes whereas bar-ed indices i denote the "lower" ones, see Fig. 2a; we will use uppercase indices I if we refer to both indifferently. Here, w i and t i are the hopping amplitudes. For a uniform system with w i w and t i ¼ t, one obtains a gapless point for w ¼ t separating a topological phase for w<t from the trivial phase for w>t. The former features topologically protected edge modes, which are fermionic in nature. The second-quantized Hamiltonian can be encoded by a matrix H SSH viaĤ SSH = Φ † H SSH Φ with the pseudospinor Φ = c 1 ; c 1 ; ; c L ; c L À Á T . For real hopping amplitudes w i and t i , the Hamiltonian exhibits time reversal symmetry T ¼ K, where K denotes complex conjugation. Furthermore, time reversal T together with the sublattice symmetry S ¼ U SSH C , represented by the unitary yields the particle-hole (PH) symmetry C ¼ KU SSH C with C 2 = +1. Hence, the SSH chain is in symmetry class BDI of the Altland-Zirnbauer classification. [17][18][19][20] In one dimension, this allows for the definition of a Z topological invariant, [11][12][13][14] which is responsible for the emergence of the disorder-resilient edge modes bound to the open ends of the chain in the topological phase.
The implementation of an analog system with bosonic degrees of freedom is straightforward: We replace the fermionic operators by bosonic ones, i.e., c positive energy shift δ > 0 is required to enforce positivity on the matrix H bSSH and accounts for the energy ω I of each bosonic mode b I . The bosonic HamiltonianĤ bSSH features the same single-particle band structure as the original fermionic chain, and exhibits the same topological properties and topological quantum numbers. Therefore, it gives rise to the same edge modes. Note that these are statements about single-particle physics where statistics is not relevant. To satisfy the PH symmetry, one must respect the sublattice symmetry (6) that protects the topological invariant. This is equivalent to the constraints for all sites i and i, i.e., all bosonic modes must have the same energy. Note that there are no constraints on the couplings w i and t i .
To complete the picture, we point out that the bosonic realization of the Majorana chain H bMC is unitarily equivalent to that of the SSH chain H bSSH ; this is shown in the Supplementary Information Section 2. However, despite their unitary equivalence, the bosonic networks impose different symmetry constraints on the coupling Hamiltonians. Therefore, depending on the experimental constraints on protecting the symmetry, it may be advantageous to implement one or the other of the unitary equivalent models.
Protocol for state transfer Next, we discuss the protocol for state transfer. One of the key features of the protocol is that we require only global (translational invariant), time-dependent tuning of the hopping amplitudes w i ; in particular, single-site addressability and control is not required. The goal of the protocol is to coherently transfer a single, quantized excitation from one of the localized edge modes to the other by means of an adiabatic variation of the couplings w in H bSSH .
The crucial point we exploit for state transfer is that in finite systems and in the topological phase (for 0<w<t), there is a finite overlap between the edge modes due to their exponential extension into the bulk. While deep in the topological phase this overlap is exponentially suppressed, it can be strongly enhanced by tuning w closer to t from below, allowing for tunneling Network inspired by the Majorana chain (unitarily equivalent to the SSH setup) and described byĤ bMC t ð Þ; the homogeneous eigenfrequency differences δωðtÞ ω À ðtÞ À ω þ ðtÞ are tunable. For details we refer the reader to Supplementary Information Section 1. c Model with two artificial edge modes separated by a simple tunneling barrier of tunable eigenfrequencies ω barrier (t); described by the HamiltonianĤ bB ðtÞ. d The simplest model, based on free propagation of excitations and described byĤ bP ðtÞ; all couplings tðtÞ are tuned simultaneously Topological networks for quantum communication N Lang and HP Büchler between the macroscopically separated edge modes. In order to prevent scattering into bulk modes, edge physics and bulk physics have to be adiabatically decoupled. This can be achieved by tuning w smoothly (and slowly, see below) towards the topological transition and return to the "sweet spot" w ¼ 0 afterwards to relocalize (and thereby decouple) the edge modes. To this end, we introduce a time-dependent hopping rate wðtÞ ¼ w max F ðtÞ giving rise to the time-dependent network Hamiltonian H bSSH ðtÞ with perfectly localized edge modes at t = 0 and t = τ. For simplicity, we choose for the adiabatic process the smooth pulse shape Here, τ denotes the characteristic time scale for the pulse. The exact pulse shape does have influence on the performance of the protocol, and setup-specific optimizations may yield quantitatively better results, see below.
We analyze the transfer efficiency and its dependence on the parameters τ and w max by numerically evaluating the full unitary time evolution. We start with an excitation in the left edge modẽ b C , and are interested in the transfer to the right edge modeb T , i.e., the state 0; 0; 1 j i. The transfer is characterized by the overlap Here, O ! 0 denotes the transfer fidelity, while φ is the relative phase accumulated during the adiabatic process. U τ w max ð Þ is the unitary time evolution operator at time t = τ, which depends parametrically on w max . To quantify the degree of adiabaticity, we introduce another characteristic parameter which describes the total edge mode population, Deviations of E from unity indicate undesired losses into bulk modes.
The qualitative results for the transfer are shown in Fig. 3a for an optimized set of parameters τ and w max , i.e., slow transfer with τ ) t À1 (ħ = 1). As expected, we find perfect transfer and decoupling of edge and bulk modes. The overall performance is quantified by O, φ and E, and depends on how close the protocol parameter w max is to the critical value t, the size of the system L, and the global time scale τ; see Fig. 4a, b for a chain of length L = 5. The edge weight E (gray background tiles) equals unity almost everywhere, except for very fast protocols and tiny bulkedge gaps. We observe quite generally that for a smooth pulse shape like F , the adiabatic bulk-edge decoupling is rather generically established in the topological setup. The size of colored squares denotes the transfer O, while the color accounts for the value of the phase φ accumulated during the transfer (measured in the rotating frame of the localized edge modes). We find several disjoint branches with O % 1 corresponding to an increasing number of round trips of the excitation; see Fig. 4b. The outermost branch allows for the fastest and most robust transfer, and is therefore the desired parameter regime to perform quantum operations. However, the most striking property of this setup is a fixed phase φ accumulated during a transfer, i.e., φ = ±π/ 2. The sign depends on the number of round trips and on the parity of the chain length L, see Supplementary Information Sections 4 and 5 for an explanation. This remarkable feature is a peculiarity of the PH symmetric topological setup and in general violated for other setups (see comparison below). A motivation for the relation of PH symmetry and fixed phase is presented in Supplementary Information Section 4.
As a concluding remark, note that there may be residual couplings w min ( w max that cannot be switched off for t < 0 and t > τ. Weak residual couplings w min (compared to t) can be tolerable on the relevant timescales as they are exponentially suppressed with the qubit distance L in the topological setup, while controlled coupling is always possible for w max ! t.

Scaling and adiabaticity
An important aspect for quantum information processing over large distances is the scalability of the protocol with separation L between the qubits. We identify the two relevant time scales of the transfer protocol: The inverse edge mode splitting ΔE À1 edge , which determines the time for a state transfer between the two edge states, and second, the inverse of the bulk-edge separation ΔE À1 bulk , which gives a lower bound on the protocol time scale due to the required adiabatic bulk-edge decoupling.
We start by considering the scaling of these energies when the topological phase transition w ¼ t is approached from the topological phase w<t. In the limit L → ∞, the eigenvalues of H bSSH derive from the transcendental equation (see Supplementary Information Section 5 for the derivation)   Fig. 5a, b. Notably, both energies scale as~1/L with the result that their ratio R = ΔE bulk /ΔE edge saturates for large L. For a fixed ratio R of the two energy scales, it is therefore required to approach the critical point as t À w $ 1=L. E. g., for R = 10 one finds t À w % 3:3=L, see Fig. 5b.
This result demonstrates that if one requires an adiabatic protocol wðtÞ with a fixed minimum ratio R min ¼ ΔE min bulk =ΔE max edge at the minimal distance Δw min ¼ t À w max from the critical coupling, then the time τ for the protocol scales as τ $ L; see Fig. 5c, d for simulations. The latter corresponds to the optimal scaling achievable since the Lieb-Robinson bound predicts a finite propagation speed for information. 56 However, we still need to adiabatically decouple bulk from edge modes, since losses to the bulk cannot be refocused in edge modes via a global tuning of parameters. A common (and conservative) estimate for adiabaticity then reads τ ≳ ΔE min bulk À Á À2 , which leads to the non-optimal scaling condition τ~L 2 . We demonstrate in the following, that a much better scaling is achievable by a rigorous estimation of the adiabaticity condition.
Then, the non-adiabatic losses to the bulk can be rigorously upper bounded 57 by where with I 1 À P, C 1,2 numerical constants, and ε L ¼ Δw 0 min = L À Δw 0 min À Á ; see the Supplementary Information Section 6 for details. Note that Ið1=2Þ ¼ 0 and ε L $ 1=L so that C L ½P diverges for L → ∞ in general. In order to bound the bulk losses, the scaling of τ has to match the scaling of C L ½P.
For PðsÞ ¼ FðsÞ ¼ sin 2 ðπsÞ we find C L ½F $ L 1þ 1 2 so that a scaling of τ $ L 1þ 1 2 is necessary for bulk-edge decoupling in the limit of long chains, see Fig. 5e. This is better than the quadratic scaling expected from the minimal gap ΔE min bulk . Unfortunately, the optimal scaling τ $ L allowed by the Lieb-Robinson bound cannot be reached by the unoptimized pulse F . However, in Supplementary Information Section 6 we prove that there is a sequence of polynomial pulses P n such that C L ½P n $ L 1þ 1 n for n ≥ 2 even integers, i.e., the scaling can be drastically improved by pulse optimization so that linear scaling can be approached to an arbitrary degree. Numerical simulations of the bulk losses for various scalings of τ reveal that they indeed follow the prescribed scaling of the rigorous upper bounds, see Fig. 5f. As a final remark, we stress that the coefficients in C L ½P n become larger with n, i.e., there is a pay-off between scaling and offset. Thus one may even benefit from pulses with poor scaling if only chains of fixed length are considered.
Benchmarking against topologically trivial setups To unveil the characteristic features of the topological setup, we contrast it with two similar but topologically trivial networks, see Fig. 2c, d. The simplest approach to envisage is based on initially decoupled modes at fixed frequency ω I ω, a homogeneous tuning of all couplings w i ¼ t i tðtÞ ¼ t max F ðtÞ, and employs the free bulk propagation of the initially localized edge modes. As shown in Fig. 3d, this approach fails to relocalize the excitation at the opposite edge due to the propagation via bulk modes; such a protocol would require either fine tuning of the pulse shape via optimal control and/or local addressability of all couplings within the network. It is therefore not competitive against the topological setup.
A more sophisticated approach mimics the presence of localized edge modes by a large tunnel barrier: The two modes at the edge have fixed frequency ω 1 ¼ ω L ¼ ω edge ¼ const, and are separated from each other by a "potential barrier" of modes with tunable frequencies ω barrier (t) and fixed couplings t ¼ w. In analogy to the topological setup, this network exhibits exponentially localized edge modes. Transfer is again achieved by lowering the excitation gap to the bulk modes to allow for tunneling Fig. 4 State transfer-quantitative results. In a, c we plot the figures of merit for transfers driven byĤ bSSH ðtÞ andĤ bB ðtÞ in dependence of the protocol timescale τ and the distance from criticality, namely Δw min ¼ t À w max and Δω min barrier ¼ ω min barrier À ω edge . The diameter of the gray background squares encodes the edge weight E; a thinning of the gray background tiling therefore indicates a loss of adiabaticity. In the shown parameter regimes, however, the edge weight is almost everywhere close to unity as there is barely any loss to bulk excitations (except for regions of fast protocols close to criticality). The diameter (color) of the colored squares encodes the transfer O (phase φ) after the protocol reached its final state (φ is measured in the rotating frame of the localized edge modes). a shows results for a topological SSH setup of size L = 5. c shows the corresponding data for a trivial tunneling barrier setup of size L = 5. In b, d we plot O, φ, and E along the dashed slices in a, c, respectively. Note that the phase is fixed for the topological setup: φ = ±π/2 Topological networks for quantum communication N Lang and HP Büchler between the edges. The protocol of this scheme reads ω barrier ðtÞ ¼ ω max barrier þ ω min barrier À ω max barrier À Á Á FðtÞ (16) where ω max barrier ) ω edge will be kept fixed and ω min barrier >ω edge is a tunable protocol parameter. The bosonic network Hamiltonian of this scheme is denoted byĤ bB ðtÞ.
As shown in Fig. 3c, the tunneling approach still allows for near perfect transfer for optimal parameters and long times. However, the quantitative comparison in Fig. 4c shows that the trivial tunneling approach requires longer time scales of the protocol and is more sensitive to bulk losses. Even then the adiabatic decoupling is much harder to achieve withĤ bB ðtÞ than witĥ H bSSH ðtÞ, as the plots of E along the dashed cuts in Fig. 4b, d reveal.
However, the most striking difference is the phase accumulated during the protocol: for the trivial setup, it is highly sensitive to both parameters. This is expected for a generic adiabatic protocol and is in stark contrast to the topological setup. The reason for this qualitative difference is rooted in the PH symmetry of the SSH setup which gives rise to the symmetric band structure depicted in Fig. 3a, as opposed to the asymmetric band structure of the barrier setup in Fig. 3c. As a consequence, we find that even for the ideal, topologically trivial setup, adiabatic protocols are Fixing the ratio R = ΔE bulk / ΔE edge (here R = 10) determines Δw 0 via the intersection marked with a circle (here Δw 0 % 3:3). c Simulations of transfer fidelity O (solid red) and edge weight E (dashed black) for system size L = 5, protocol timescale τ = 50 and pulse wðtÞ ¼ t À Δw min ð ÞÁF ð tÞ as function of Δw min . Optimal transfer for fixed L and τ is found numerically for Δw min % 0:26 with bulk loss 1 À E % 1 Á 10 À4 . d The same as in c for doubled size L = 10 and timescale τ = 100. Now, optimal transfer is achieved for Δw min % 0:26=2 ¼ 0:13 with bulk loss 1 À E % 2 Á 10 À5 . e Rigorous upper bounds τ À1 C L ½P for PðsÞ ¼ FðsÞ ¼ sin 2 ðπsÞ and τ = τ 0 ⋅ L 1+α with α ¼ 0; 1 2 ; 1 and τ 0 = 100, Δw 0 min ¼ 3:3. A scaling τ $ L 1þ 1 2 yields constant bulk losses for L → ∞. f Simulations of the bulk losses 0 1 À E 1 for the parameters in e without tuning for optimal transfer. We find that the scaling follows the corresponding upper bounds. Note that the loss was chosen large (~50% for L = 10) for illustrative purposes and can be controlled via τ 0 (here τ 0 = 1, 0.3, 0.1 for α ¼ 0; 1 2 ; 1) unsuitable as the sensitivity of the phase increases for longer wires: a transfer preserving quantum coherence requires fine tuning of the shape of the transfer pulse. This effect becomes even more drastic in the presence of disorder.

Effects of disorder and symmetry protection
The unique features of the topological setup become even more apparent in the presence of disorder and/or imperfections in the preparation. Here, we focus on quenched disorder on the time scales for a transfer. The disorder is described as Gaussian noise with dimensionless standard deviation p acting on the parameters in the Hamiltonian, i.e., for the onsite hopping we have w i h i h i¼ w and w 2 i À w 2 ¼ p 2 w 2 with h i h i the disorder average. In the following, two classes of disorder will be of interest: PH symmetric disorder affects only mode couplings, but assumes perfect mode frequencies; recall Eq. (8). In contrast, PH breaking disorder affects both mode couplings and frequencies.
For the topologically trivial setup, both types of disorder give rise to Anderson localization of the "artificial" edge modes, see Fig. 6a-c. The transfer protocol, however, relies on the delocalization of the edge modes, which is prohibited by Anderson localization. As a consequence, disorder leads to a significant reduction of transfer fidelity and increased bulk losses, Fig. 6d. Furthermore, the phase φ accumulated during the transfer strongly fluctuates for each disorder realization; more details on this aspect are given in Supplementary Information Section 4.
In contrast, for the topological SSH setup, the PH symmetric disorder respects the protecting symmetry. Then, Anderson localization of the edge modes is forbidden by a topological obstruction 11,13,14 and the required overlap between the two edge modes can be established, see Fig. 6b. As a consequence, the transfer can still be performed perfectly with a fixed phase φ = ±π/ 2. However, this requires that for each disorder realization one is allowed to adapt the transfer time τ of the protocol. In an experimental setup this corresponds, for example, to imperfections in sample preparation, which can be overcome by calibrating the setup and the transfer protocol beforehand. In turn, the PH breaking disorder leads also in the topological setup to localization of the edge modes and a reduction of transfer fidelity, see Fig. 6d.  Application: CP gate As an application, we demonstrate how the proposed state transfer, protected by PH symmetry, can be employed for a CP gate between two remote qubits that are coupled to the local edge modes of the topological SSH network. The complete protocol for a CP gate between the target qubit T and the control qubit C is based on a well-known scheme that makes use of auxiliary levels a j i T /C . 58 We focus on the setup shown in Fig. 1b for the qubits and the auxiliary states with the coupling Hamiltonian between the qubits and the edge modes of the SSH chain given by Eq. (2). The full protocol for the CP gate follows the procedure in Eq. (3) and is described in the following, see Fig. 7.
First, the full sequence U CP leaves the states 0 j i C 0 j i T and 1 j i C 0 j i T invariant because there are no excitations in the network. On the other hand, for 0 j i C 1 j i T and 1 j i C 1 j i T the first π-pulse Π T on the target qubit maps the state 1 j i T to a bosonic excitation in the right edge modeb T with the phase −π/2. It is important that this operation is performed slowly compared to the energy gap to bulk excitations in the SSH chain: Then, energy conservation allows one to only address the coupling to the edge states and suppress admixture of bulk excitations. The subsequent transfer of the excitation to the left edge T C↔T implies an additional phase π/ 2 + Lπ. The full Rabi cycle Π 2 C provides a phase π if and only if the control qubit is in state 1 j i C . The subsequent transfer back T C↔T and the π-pulse Π T provide additional phases π/2 + Lπ and −π/2. Therefore, the full protocol implements the mapping 1 j i C 1 j i T → − 1 j i C 1 j i T while all orthogonal states remain invariant. Given quantum coherence during the protocol, this realizes a controlled phase gate with phase π. A full numerical time evolution for the state 1 j i C 1 j i T is shown in Fig. 7c, d, confirming the above argumentation.
Finally, we point out that the linearity of the network implies that the transfer takes place for each excitation of edge modes independently. I.e., if n C ; n T j idenotes the state with n C excitations in the left edge modeb C and n T excitations in the right edge modeb T , the transfer operation T C↔T implements the mapping n C ; n T j i! ð ± iÞ nCþnT n T ; n C j i (for even/odd L). This observation immediately implies that the unitary operation U SWAP in Eq. (4) swaps the qubits.
Extension to 2D networks of coupled qubits An important aspect of the SSH chain is its symmetry class BDI (due to the real hopping amplitudes w i and t i ), so that the setup is characterized by a Z topological invariant. This allows us to extend the analysis to two-dimensional setups by placing several SSH chains parallel to each other and adding couplings between them, see Fig. 1c for a possible realization. As long as these couplings are real and respect the sublattice symmetry, the setup is still topologically protected and each chain endpoint carries an edge mode. Possible relaxations of the symmetry constraints are presented in Supplementary Information Section 3.
Due to the sublattice symmetry, there are two types of modes, "even" and "odd" ones, depending on how they transform under S ¼ U SSH C [recall Eq. (6)]. Edge mode pairs of different types can communicate efficiently by tuning the couplings along a connecting path that resembles the one-dimensional SSH setup; this is illustrated in Fig. 1c. Because of the bulk gap, this procedure is very robust and couplings that deviate from the desired path ("stray couplings") have no detrimental effect on the state transfer, as long as they respect the symmetries and do not couple to other edge modes. This can be guaranteed by a modification of the setup such that qubits and edge modes are relocated at the end of one-dimensional chains that emanate from the 2D network, see Fig. 1d. Remarkably, the setup allows for an enhancement of the edge mode overlap by tuning the couplings of the 2D bulk globally instead of tracing out a particular path that connects the qubits. Therefore, the minimal experimental requirement is the individual addressability of each branch that connects an edge mode to the 2D bulk. Since the edge modes can now be separated, this constraint is very weak, and in general already satisfied by the requirement of local gate operations on the qubit.
Note that coupling edge modes of the same type is obstructed by the sublattice symmetry, as indicated by the lower path in Fig. 1c and proven in Supplementary Information Section 3. However, the implementation of an exchange of qubits via U SWAP facilitates the application of the controlled phase gate U CP between any pair in the network: If the two qubits couple to edge modes of different types, one can directly perform the CP gate between them. Conversely, if the qubits couple to edge modes of the same type, one first performs an exchange U SWAP with an arbitrary qubit of the opposite type, applies the CP gate, and maps the qubit back by another exchange.

DISCUSSION
We have demonstrated that a topological network consisting of linearly coupled bosonic degrees of freedom, capable of carrying single, quantized excitations, allows for efficient quantum communication between distant qubits. Weak addressability of each branch within the network is sufficient and no local addressability of individual sites is required. Remarkably, the time scale for the operations scales almost linearly with the distance between the qubits. Furthermore, topological protection guarantees robustness  Fig. 1a, b. b The complete pulse sequence consisting of two π pulses Π T to map the target qubit to and from edge modeb T , a 2π pulse Π 2 C to perform the actual CP gate, and two-edge mode tunneling pulses T C↔T . c Numerical singleparticle evolution for the two-qubit basis state 1 j i C 1 j i T for a chain of length L = 10. The density plot encodes the squared single-particle amplitude, where the upper and lower edges correspond to the logical states 1 j i C and 1 j i T with their adjacent boundary modes 1 j ib C ¼ 1; 0; 0 j i¼b y C 0 j i and 1 j ib T ¼ 0; 0; 1 j i¼b y T 0 j i. d Square of the absolute value of the overlaps with 1 j i T (solid black), 1 j ib T (solid red), 1 j ib C (solid blue), and 1 j i C (dashed black). The relative phases w.r.t. the target qubit are shown as insets. Note that the scales for the Rabi pulses and the topological state transfer differ by a factor of τ Topological networks for quantum communication N Lang and HP Büchler against quenched disorder in the setup by evading Anderson localization of the edge modes. In summary, we have shown that the unique properties of (quasi) one-dimensional topological systems can be harvested for efficient quantum communication between qubits. These benefits come with the price of higher complexity in realization and preparation as the coupling parameters have to respect the symmetries protecting the topological invariants.

METHODS
For more details on the derivations, see the Supplementary Information.

Data availability
The numerical code and the datasets used for this study are available from the corresponding author on reasonable request.