Topological networks for quantum communication between distant qubits

TL;DR

This paper proposes a topological bosonic network for robust, efficient quantum communication between distant qubits, featuring global control, disorder resilience, and near-linear scaling with distance.

Contribution

It introduces a novel topological network architecture enabling scalable, robust quantum communication with minimal local control and high resilience to disorder.

Findings

Robust quantum information transfer over large distances.

Global parameter variations suffice for all operations.

Communication time scales nearly linearly with distance.

Abstract

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.

Equations

$${\hat{H}}_n=\sum _{i,j}{b}_i^{†}{H}_{ij}{b}_j.$$

$${\hat{H}}_{\mathrm{T}}(t)=g_{\mathrm{T}}(t)\left[{\tilde{b}}_{\mathrm{T}}^{†}|a⟩⟨1{|}_{\mathrm{T}}+{\tilde{b}}_{\mathrm{T}}|1⟩⟨a{|}_{\mathrm{T}}\right].$$

$${\mathcal{U}}_{\mathrm{CP}}={\Pi }_{\mathrm{T}}\circ T_{\mathrm{C}\leftrightarrow \mathrm{T}}\circ {\Pi }_{\mathrm{C}}^2\circ T_{\mathrm{C}\leftrightarrow \mathrm{T}}\circ {\Pi }_{\mathrm{T}}.$$

$${\mathcal{U}}_{\mathrm{SWAP}}={\mathcal{A}}_{\mathrm{C}}\circ {\Pi }_{\mathrm{T}}\circ {\Pi }_{\mathrm{C}}\circ T_{\mathrm{C}\leftrightarrow \mathrm{T}}\circ {\Pi }_{\mathrm{C}}\circ {\Pi }_{\mathrm{T}}\circ {\mathcal{A}}_{\mathrm{C}}.$$

$${\hat{H}}_{\mathrm{SSH}}=\sum _{i=1}^L{w}_i{c_i}^{†}{c}_{\overline{i}}+\sum _{i=1}^{L-1}{t}_i{c_{\overline{i}}}^{†}{c}_{i+1}+h.c.$$

$$U_C^{\mathrm{SSH}}=\text{mathds}{1}_{L\times L}\otimes \left(\begin{array}{cc}1& 0\\ 0& -1\end{array}\right),$$

$${\hat{H}}_{\mathrm{bSSH}}\equiv {\mathbf{\xi }}^{†}\left[H_{\mathrm{SSH}}+\delta \text{mathds}1\right]\mathbf{\xi }\equiv {\mathbf{\xi }}^{†}{H}_{\mathrm{bSSH}}\mathbf{\xi }$$

$${\omega }_i=\delta ={\omega }_{\overline{i}}$$

$$\mathcal{F}(t)={\sin }^2\left(\pi t/\tau \right)\text{for}0\le t\le \tau .$$

$$|{\Psi }_0⟩={|1⟩}_1\otimes {|0,\dots ,0⟩}_{\mathrm{bulk}}\otimes {|0⟩}_{\overline{L}}\equiv |1,0,0⟩,$$

$$⟨0,0,1|U_{\tau }({\bar{w}}_{\max })|1,0,0⟩\equiv \sqrt{\mathcal{O}}{e}^{i\phi }.$$

$$\mathcal{E}=\mathcal{O}+|⟨1,0,0|U_{\tau }({\bar{w}}_{\max })|1,0,0⟩{|}^2.$$

$$\frac{\Delta {\bar{w}}^{\prime }-\sqrt{\Delta {\bar{w}}^{\prime 2}-{\lambda }^{\prime 2}}}{\Delta {\bar{w}}^{\prime }+\sqrt{\Delta {\bar{w}}^{\prime 2}-{\lambda }^{\prime 2}}}=e^{-2\sqrt{\Delta {\bar{w}}^{\prime 2}-{\lambda }^{\prime 2}}}$$

$$1-\mathcal{E}\le {\left(\frac{C_L[\mathcal{P}]}{\tau }\right)}^2,$$

$$C_L[\mathcal{P}]={\int }_0^1\mathrm{d}s\frac{C_1|{\mathcal{I}}^{\prime \prime }|}{{\left({\epsilon }_L+\mathcal{I}\right)}^2}+{\int }_0^1\mathrm{d}s\frac{C_2|{\mathcal{I}}^{\prime }{|}^2}{{\left({\epsilon }_L+\mathcal{I}\right)}^3}$$

$${\omega }_{\mathrm{barrier}}(t)={\omega }_{\mathrm{barrier}}^{\max }+({\omega }_{\mathrm{barrier}}^{\min }-{\omega }_{\mathrm{barrier}}^{\max })\cdot \mathcal{F}(t)$$

Full text

I Introduction

Systems characterized by topological invariants are well known to exhibit unique properties with potential applications in quantum information processing and engineering Nayak et al. (2008). Ever since the first experimental observation of the integer quantum hall effect Ando et al. (1975); v. Klitzing et al. (1980); Laughlin (1981), many other condensed matter systems have been identified and experimentally characterized, such as fractional quantum Hall fluids Tsui et al. (1982); Stormer et al. (1983); Willett et al. (1987) and topological insulators and superconductors Konig et al. (2007); Hsieh et al. (2008); Xu et al. (2014); Schnyder et al. (2008); Kitaev (2009); Schnyder et al. (2009); Ryu et al. (2010); Hasan and Kane (2010); Qi and Zhang (2011). 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 Wigner (1951, 1958); Dyson (1962); Altland and Zirnbauer (1997). This concept can be straightforwardly generalized to bosonic setups as well as classical systems Hafezi et al. (2013); Susstrunk and Huber (2015); Nash et al. (2015), 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 Hafezi et al. (2013); Rechtsman et al. (2013); Ling et al. (2015) and the radio-frequency Ningyuan et al. (2015) regime, as well as with classical coupled harmonic oscillators Susstrunk and Huber (2015); Nash et al. (2015); Süsstrunk et al. (2017), and with cold atomic gases Atala et al. (2013); Jotzu et al. (2014); Aidelsburger et al. (2014); Mancini et al. (2015); Stuhl et al. (2015); Duca et al. (2015); Lohse et al. (2015). These experimental advances have been prepared and are supported by many theoretical proposals, e.g. Haldane and Raghu (2008); Koch et al. (2010); Hafezi et al. (2011); Berg et al. (2011); Yannopapas (2012); Kane and Lubensky (2013); Lu et al. (2014); Kariyado and Hatsugai (2015); Wang et al. (2015a, b); Peter et al. (2015). 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 Yao et al. (2013); Dlaska et al. (2017) 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 $\mathbb{Z}$ topological index Zak (1989), and derived from paradigmatic systems such as Kitaev’s Majorana chain Kitaev (2001) and the Su-Schrieffer-Heeger (SSH) model Su et al. (1979). 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 favourably 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 Hafezi et al. (2013); Ningyuan et al. (2015); Susstrunk and Huber (2015); Süsstrunk et al. (2017). 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. 1 (a) 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

[TABLE]

Here, $b^{{\dagger}}_{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 mode $\tilde{b}_{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. 1 (b). There the coupling Hamiltonian for the target qubit $\rm T$ takes the form (within the rotating wave approximation)

[TABLE]

The coupling $g_{\rm\scriptscriptstyle T}(t)$ is controlled by external laser fields and allows for the application of $\pi$-pulses between the qubit state $\left|1\right>_{\rm\scriptscriptstyle T}$ and the edge mode $\tilde{b}_{\rm\scriptscriptstyle T}$, i.e., the emission of a photon into the edge mode $\tilde{b}_{\rm\scriptscriptstyle T}$ from state $\left|1\right>_{\rm\scriptscriptstyle T}$ is accompanied by a transition into the auxiliary state $\left|a\right>_{\rm\scriptscriptstyle T}$; in the following, we denote such a $\pi$-pulse at edge $p$ by the unitary operation $\Pi_{p}$. Note that the Hamiltonian $\hat{H}_{\rm\scriptscriptstyle C}(t)$ for the control qubit $\rm C$ is similar, with the role of $\left|1\right>_{\rm\scriptscriptstyle C}$ and $\left|a\right>_{\rm\scriptscriptstyle 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\leftrightarrow q}$. As an example, the protocol for a controlled phase (CP) gate between a control qubit at position ${\rm C}$ and a target qubit at position ${\rm T}$ reads

[TABLE]

Another example, the transport of a control qubit to a target position, is simply described by the protocol

[TABLE]

Here, $\mathcal{A}_{p}$ denotes the exchange of the two states $|1\rangle_{p}$ and $|a\rangle_{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_{\rm\scriptscriptstyle C\leftrightarrow\rm\scriptscriptstyle 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. 1 (a): 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 $\pi$-pulse $T_{\rm\scriptscriptstyle C\leftrightarrow\rm\scriptscriptstyle 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 Hafezi et al. (2013), alternative setups are coupled optical cavities and circuit QED systems Houck et al. (2012); Carusotto and Ciuti (2013) as well as trapped polar molecules or Rydberg atoms with a coupling mediated by dipolar exchange interactions Peter et al. (2015), while the local qubits can be artificial atoms Lukin et al. (2001), NV centers in diamond Doherty et al. (2013), or trapped ions Duan and Monroe (2010).

II Results

II.1 Topological network

We start with a description of the requirements on the quantum network $\hat{H}_{n}$ to exhibit topologically protected edge modes in a one-dimensional chain with two edges, as illustrated in Fig. 1 (a); 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 Kitaev (2001). This model is closely related to the Su-Schrieffer-Heeger (SSH) model Su et al. (1979) in the single-particle picture. It turns out that the necessary steps to translate these models into our bosonic quantum network language are more conveniently performed for the SSH model; the discussion of the Majorana chain and its relation to the SSH chain is presented in the Supplementary Information Sections I and II.

The SSH model on a chain with $L$ sites and open boundaries is described by the Hamiltonian

[TABLE]

with the two fermion operators $c_{i}$ and $c_{\overline{i}}$ on each site. Note that the indices $i$ label the “upper” fermionic modes whereas bar-ed indices ${\overline{i}}$ denote the “lower” ones, see Fig. 2 (a); we will use upper-case 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}\equiv\bar{w}$ and $t_{i}=\bar{t}$, one obtains a gapless point for $\bar{w}=\bar{t}$ separating a topological phase for $\bar{w}<\bar{t}$ from the trivial phase for $\bar{w}>\bar{t}$. The former features topologically protected edge modes which are fermionic in nature. The second-quantized Hamiltonian can be encoded by a matrix $H_{\rm\scriptscriptstyle SSH}$ via $\hat{H}_{\rm\scriptscriptstyle SSH}=\bm{\Phi}^{\dagger}H_{\rm\scriptscriptstyle SSH}\bm{\Phi}$ with the pseudo spinor $\bm{\Phi}={\left(c_{1},c_{\overline{1}},\dots,c_{L},c_{\overline{L}}\right)}^{T}$. For real hopping amplitudes $w_{i}$ and $t_{i}$, the Hamiltonian exhibits time reversal symmetry $T=\mathcal{K}$ where $\mathcal{K}$ denotes complex conjugation. Furthermore, time reversal $T$ together with the sublattice symmetry $S=U_{C}^{\rm\scriptscriptstyle SSH}$, represented by the unitary

[TABLE]

yields the particle-hole (PH) symmetry $C=\mathcal{K}U_{C}^{\rm\scriptscriptstyle SSH}$ with $C^{2}=+\mathds{1}$. Hence, the SSH chain is in symmetry class BDI of the Altland-Zirnbauer classification Wigner (1951, 1958); Dyson (1962); Altland and Zirnbauer (1997). In one dimension, this allows for the definition of a $\mathbb{Z}$ topological invariant Schnyder et al. (2008); Kitaev (2009); Schnyder et al. (2009); Ryu et al. (2010) 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 analogue system with bosonic degrees of freedom is straightforward: We replace the fermionic operators by bosonic ones, i.e., $c_{I}^{({\dagger})}\rightarrow b_{I}^{({\dagger})}$. The bosonic Hamiltonian takes the form

[TABLE]

with $\bm{\xi}={\left(b_{1},b_{\overline{1}},\ldots,b_{L},b_{\overline{L}}\right)}^{T}$ and $H_{\rm\scriptscriptstyle bSSH}=H_{\rm\scriptscriptstyle SSH}+\delta\mathds{1}$. The constant positive energy shift $\delta>0$ is required to enforce positivity on the matrix $H_{\rm\scriptscriptstyle bSSH}$ and accounts for the energy $\omega_{I}$ of each bosonic mode $b_{I}$. The bosonic Hamiltonian $\hat{H}_{\rm\scriptscriptstyle 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

[TABLE]

for all sites $i$ and ${\overline{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_{\rm\scriptscriptstyle bMC}$ is unitarily equivalent to that of the SSH chain $H_{\rm\scriptscriptstyle bSSH}$; this is shown in the Supplementary Information Section II. 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.

II.2 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 $\bar{w}$ in $H_{\rm\scriptscriptstyle bSSH}$.

The crucial point we exploit for state transfer is that in finite systems and in the topological phase (for $0<\bar{w}<\bar{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 be strongly enhanced by tuning $\bar{w}$ closer to $\bar{t}$ from below, allowing for tunnelling between the macroscopically separated edge modes. In order to prevent scattering into bulk modes, edge- and bulk physics have to be adiabatically decoupled. This can be achieved by tuning $\bar{w}$ smoothly (and slowly, see below) towards the topological transition and return to the “sweet spot” $\bar{w}=0$ afterwards to relocalize (and thereby decouple) the edge modes. To this end, we introduce a time dependent hopping rate $\bar{w}(t)=\bar{w}_{\rm\scriptscriptstyle max}\,\mathcal{F}(t)$ giving rise to the time-dependent network Hamiltonian $\hat{H}_{\rm\scriptscriptstyle bSSH}(t)$ with perfectly localized edge modes at $t=0$ and $t=\tau$. For simplicity, we choose for the adiabatic process the smooth pulse shape

[TABLE]

Here, $\tau$ 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 analyse the transfer efficiency and its dependence on the parameters $\tau$ and $\bar{w}_{\rm\scriptscriptstyle max}$ by numerically evaluating the full unitary time evolution. We start with an excitation in the left edge mode $\tilde{b}_{\rm\scriptscriptstyle C}$,

[TABLE]

and are interested in the transfer to the right edge mode $\tilde{b}_{\rm\scriptscriptstyle T}$, i.e., the state $\left|0,\bm{0},1\right>$. The transfer is characterized by the overlap

[TABLE]

Here, $\mathcal{O}\geq 0$ denotes the transfer fidelity, while $\varphi$ is the relative phase accumulated during the adiabatic process. $U_{\tau}(\bar{w}_{\rm\scriptscriptstyle max})$ is the unitary time evolution operator at time $t=\tau$ which depends parametrically on $\bar{w}_{\rm\scriptscriptstyle max}$. To quantify the degree of adiabaticity, we introduce another characteristic parameter which describes the total edge mode population,

[TABLE]

Deviations of $\mathcal{E}$ from unity indicate undesired losses into bulk modes.

The qualitative results for the transfer are shown in Fig. 3 (a) for an optimized set of parameters $\tau$ and $\bar{w}_{\rm\scriptscriptstyle max}$, i.e., slow transfer with $\tau\gg\bar{t}^{-1}$ ($\hbar=1$). As expected, we find perfect transfer and decoupling of edge and bulk modes. The overall performance is quantified by $\mathcal{O}$, $\varphi$ and $\mathcal{E}$, and depends on how close the protocol parameter $\bar{w}_{\rm\scriptscriptstyle max}$ is to the critical value $\bar{t}$, the size of the system $L$, and the global time scale $\tau$; see Fig. 4 (a) and (b) for a chain of length $L=5$. The edge weight $\mathcal{E}$ (grey background tiles) equals unity almost everywhere, except for very fast protocols and tiny bulk-edge gaps. We observe quite generally that for a smooth pulse shape like $\mathcal{F}$, the adiabatic bulk-edge decoupling is rather generically established in the topological setup. The size of coloured squares denotes the transfer $\mathcal{O}$, while the colour accounts for the value of the phase $\varphi$ accumulated during the transfer (measured in the rotating frame of the localized edge modes). We find several disjoint branches with $\mathcal{O}\approx 1$ corresponding to an increasing number of round trips of the excitation; see Fig. 4 (b). 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 $\varphi$ accumulated during a transfer, i.e., $\varphi=\pm\pi/2$. The sign depends on the number of round trips and on the parity of the chain length $L$, see Supplementary Information Sections IV and V 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 IV.

As a concluding remark, note that there may be residual couplings $\bar{w}_{\rm\scriptscriptstyle min}\ll\bar{w}_{\rm\scriptscriptstyle max}$ that cannot be switched off for $t<0$ and $t>\tau$. Weak residual couplings $\bar{w}_{\rm\scriptscriptstyle min}$ (compared to $\bar{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 $\bar{w}_{\rm\scriptscriptstyle max}\to\bar{t}$.

II.3 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 $\Delta E_{\rm\scriptscriptstyle edge}^{-1}$ which determines the time for a state transfer between the two edge states, and second, the inverse of the bulk-edge separation $\Delta E_{\rm\scriptscriptstyle bulk}^{-1}$ 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 $\bar{w}=\bar{t}$ is approached from the topological phase $\bar{w}<\bar{t}$. In the limit $L\to\infty$, the eigenvalues of $H_{\rm\scriptscriptstyle bSSH}$ derive from the transcendental equation (see Supplementary Information Section V for the derivation)

[TABLE]

with $\Delta\bar{w}^{\prime}/L=\bar{t}-\bar{w}$ measuring the distance to the topological phase transition. The lowest two solutions $\lambda_{0}^{\prime}$ and $\lambda_{1}^{\prime}$ of equation (13) determine the relevant energies $\Delta E_{\rm\scriptscriptstyle edge}=2\lambda_{0}^{\prime}/L$ and $\Delta E_{\rm\scriptscriptstyle bulk}=|\lambda_{1}^{\prime}-\lambda_{0}^{\prime}|/L$ , see Fig. 5 (a) and (b). Notably, both energies scale as $\sim 1/L$ with the result that their ratio $R=\Delta E_{\rm\scriptscriptstyle bulk}/\Delta E_{\rm\scriptscriptstyle 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 $\bar{t}-\bar{w}\sim 1/L$. E.g., for $R=10$ one finds $\bar{t}-\bar{w}\approx 3.3/L$, see Fig. 5 (b).

This result demonstrates that if one requires an adiabatic protocol $\bar{w}(t)$ with a fixed minimum ratio $R^{\rm\scriptscriptstyle min}=\Delta E_{\rm\scriptscriptstyle bulk}^{\rm\scriptscriptstyle min}/\Delta E_{\rm\scriptscriptstyle edge}^{\rm\scriptscriptstyle max}$ at the minimal distance $\Delta\bar{w}_{\rm\scriptscriptstyle min}=\bar{t}-\bar{w}_{\rm\scriptscriptstyle max}$ from the critical coupling, then the time $\tau$ for the protocol scales as $\tau\sim L$; see Fig. 5 (c) and (d) for simulations. The latter corresponds to the optimal scaling achievable since the Lieb-Robinson bound predicts a finite propagation speed for information Lieb and Robinson (1972).

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 $\tau\gtrsim{\left(\Delta E_{\rm\scriptscriptstyle bulk}^{\rm\scriptscriptstyle min}\right)}^{-2}$, which leads to the non-optimal scaling condition $\tau\sim L^{2}$. We demonstrate in the following, that a much better scaling is achievable by a rigorous estimation of the adiabaticity condition. To this end, we parametrize the time with $s=t/\tau$, $s\in[0,1]$ and make the ansatz $\bar{w}(t)=\bar{w}_{\rm\scriptscriptstyle max}\cdot\mathcal{P}(s)$, where the generic pulse $\mathcal{P}:[0,1]\to[0,1]$ and its first derivative vanish for $s=0,1$, and $\mathcal{P}(1/2)=1$.

Then, the non-adiabatic losses to the bulk can be rigorously upper-bounded Jansen et al. (2007) by

[TABLE]

where

[TABLE]

with $\mathcal{I}\equiv 1-\mathcal{P}$, $C_{1,2}$ numerical constants, and $\varepsilon_{L}=\Delta\bar{w}_{\rm\scriptscriptstyle min}^{\prime}/(L-\Delta\bar{w}_{\rm\scriptscriptstyle min}^{\prime})$; see the Supplementary Information Section VI for details. Note that $\mathcal{I}(1/2)=0$ and $\varepsilon_{L}\sim 1/L$ so that $C_{L}[\mathcal{P}]$ diverges for $L\to\infty$ in general. In order to bound the bulk losses, the scaling of $\tau$ has to match the scaling of $C_{L}[\mathcal{P}]$.

For $\mathcal{P}(s)=\mathcal{F}(s)=\sin^{2}(\pi s)$ we find $C_{L}[\mathcal{F}]\sim L^{1+\frac{1}{2}}$ so that a scaling of $\tau\sim L^{1+\frac{1}{2}}$ is necessary for bulk-edge decoupling in the limit of long chains, see Fig. 5 (e). This is better than the quadratic scaling expected from the minimal gap $\Delta E_{\rm\scriptscriptstyle bulk}^{\rm\scriptscriptstyle min}$. Unfortunately, the optimal scaling $\tau\sim L$ allowed by the Lieb-Robinson bound cannot be reached by the unoptimized pulse $\mathcal{F}$. However, in Supplementary Information Section VI we prove that there is a sequence of polynomial pulses $\mathcal{P}_{n}$ such that $C_{L}[\mathcal{P}_{n}]\sim L^{1+\frac{1}{n}}$ for $n\geq 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 $\tau$ reveal that they indeed follow the prescribed scaling of the rigorous upper bounds, see Fig. 5 (f). As a final remark, we stress that the coefficients in $C_{L}[\mathcal{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.

II.4 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. 2 (c) and (d). The simplest approach to envisage is based on initially decoupled modes at fixed frequency $\omega_{I}\equiv\bar{\omega}$, a homogeneous tuning of all couplings $w_{i}=t_{i}\equiv\bar{t}(t)=\bar{t}_{\rm\scriptscriptstyle max}\,\mathcal{F}(t)$, and employs the free bulk propagation of the initially localized edge modes. As shown in Fig. 3 (d), 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 $\omega_{1}=\omega_{\overline{L}}=\omega_{\rm\scriptscriptstyle edge}=\textrm{const}$, and are separated from each other by a “potential barrier” of modes with tunable frequencies $\omega_{\rm\scriptscriptstyle barrier}(t)$ and fixed couplings $\bar{t}=\bar{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 tunnelling between the edges. The protocol of this scheme reads

[TABLE]

where $\omega_{\rm\scriptscriptstyle barrier}^{\rm\scriptscriptstyle max}\gg\omega_{\rm\scriptscriptstyle edge}$ will be kept fixed and $\omega_{\rm\scriptscriptstyle barrier}^{\rm\scriptscriptstyle min}>\omega_{\rm\scriptscriptstyle edge}$ is a tunable protocol parameter. The bosonic network Hamiltonian of this scheme is denoted by $\hat{H}_{\rm\scriptscriptstyle bB}(t)$.

As shown in Fig. 3 (c), the tunnelling approach still allows for near perfect transfer for optimal parameters and long times. However, the qualitative comparison in Fig. 4 (c) shows that the trivial tunnelling 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 $\hat{H}_{\rm\scriptscriptstyle bB}(t)$ than with $\hat{H}_{\rm\scriptscriptstyle bSSH}(t)$, as the plots of $\mathcal{E}$ along the dashed cuts in Fig. 4 (b) and (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. 3 (a), as opposed to the asymmetric band structure of the barrier setup in Fig. 3 (c). As a consequence, we find that even for the ideal, topologically trivial setup, adiabatic protocols are 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.

II.5 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 $\langle\langle w_{i}\rangle\rangle=\bar{w}$ and $\langle\langle w_{i}^{2}-\bar{w}^{2}\rangle\rangle=p^{2}\bar{w}^{2}$ with $\langle\langle\bullet\rangle\rangle$ 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 equation (8). In contrast, PH breaking disorder affects both mode couplings and frequencies.

For the topological trivial setup, both types of disorder give rise to Anderson localization of the “artificial” edge states, see Fig. 6 (a)–(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. 6 (d). Furthermore, the phase $\varphi$ accumulated during the transfer strongly fluctuates for each disorder realization; more details on this aspect are given in Supplementary Information Section IV.

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 Schnyder et al. (2008, 2009); Ryu et al. (2010) and the required overlap between the two edge modes can be established, see Fig. 6 (b). As a consequence, the transfer can still be performed perfectly with a fixed phase $\varphi=\pm\pi/2$. However, this requires that for each disorder realization one is allowed to adapt the transfer time $\tau$ 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. 6 (d).

II.6 Application: Controlled-phase gate

As an application, we demonstrate how the proposed state transfer, protected by PH symmetry, can be employed for a controlled-phase (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 $\rm T$ and the control qubit $\rm C$ is based on a well-known scheme that makes use of auxiliary levels $\left|a\right>_{\rm\scriptscriptstyle T/\rm\scriptscriptstyle C}$ Haack et al. (2010). We focus on the setup shown in Fig. 1 (b) for the qubits and the auxiliary states with the coupling Hamiltonian between the qubits and the edge modes of the SSH chain given by equation (2). The full protocol for the CP gate follows the procedure in equation (3) and is described in the following, see Fig. 7.

First, the full sequence $\mathcal{U}_{\rm\scriptscriptstyle{CP}}$ leaves the states $\left|0\right>_{\rm\scriptscriptstyle C}\left|0\right>_{\rm\scriptscriptstyle T}$ and $\left|1\right>_{\rm\scriptscriptstyle C}\left|0\right>_{\rm\scriptscriptstyle T}$ invariant because there are no excitations in the network. On the other hand, for $\left|0\right>_{\rm\scriptscriptstyle C}\left|1\right>_{\rm\scriptscriptstyle T}$ and $\left|1\right>_{\rm\scriptscriptstyle C}\left|1\right>_{\rm\scriptscriptstyle T}$ the first $\pi$-pulse $\Pi_{\rm\scriptscriptstyle T}$ on the target qubit maps the state $\left|1\right>_{\rm\scriptscriptstyle T}$ to a bosonic excitation in the right edge mode $\tilde{b}_{\rm\scriptscriptstyle T}$ with the phase $-\pi/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_{\rm\scriptscriptstyle C\leftrightarrow\rm\scriptscriptstyle T}$ implies an additional phase $\pi/2+L\pi$. The full Rabi cycle $\Pi_{\rm\scriptscriptstyle C}^{2}$ provides a phase $\pi$ if and only if the control qubit is in state $\left|1\right>_{\rm\scriptscriptstyle C}$. The subsequent transfer back $T_{\rm\scriptscriptstyle C\leftrightarrow\rm\scriptscriptstyle T}$ and the $\pi$-pulse $\Pi_{\rm\scriptscriptstyle T}$ provide additional phases $\pi/2+L\pi$ and $-\pi/2$. Therefore, the full protocol implements the mapping $\left|1\right>_{\rm\scriptscriptstyle C}\left|1\right>_{\rm\scriptscriptstyle T}\,\rightarrow\,-\left|1\right>_{\rm\scriptscriptstyle C}\left|1\right>_{\rm\scriptscriptstyle T}$ while all orthogonal states remain invariant. Given quantum coherence during the protocol, this realizes a controlled phase gate with phase $\pi$. A full numerical time-evolution for the state $\left|1\right>_{\rm\scriptscriptstyle C}\left|1\right>_{\rm\scriptscriptstyle T}$ is shown in Fig. 7 (c) and (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 $\left|n_{\rm\scriptscriptstyle C},n_{\rm\scriptscriptstyle T}\right>$ denotes the state with $n_{\rm\scriptscriptstyle C}$ excitations in the left edge mode $\tilde{b}_{\rm\scriptscriptstyle C}$ and $n_{\rm\scriptscriptstyle T}$ excitations in the right edge mode $\tilde{b}_{\rm\scriptscriptstyle T}$, the transfer operation $T_{\rm\scriptscriptstyle C\leftrightarrow\rm\scriptscriptstyle T}$ implements the mapping $\left|n_{\rm\scriptscriptstyle C},n_{\rm\scriptscriptstyle T}\right>\,\rightarrow\,{(\pm i)}^{n_{\rm\scriptscriptstyle C}+n_{\rm\scriptscriptstyle T}}\,\left|n_{\rm\scriptscriptstyle T},n_{\rm\scriptscriptstyle C}\right>$ (for even/odd $L$). This observation immediately implies that the unitary operation $\mathcal{U}_{\rm\scriptscriptstyle SWAP}$ in equation (4) swaps the qubits.

II.7 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 $\mathbb{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. 1 (c) 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 III.

Due to the sublattice symmetry, there are two types of modes, “even” and “odd” ones, depending on how they transform under $S=U_{C}^{\rm\scriptscriptstyle SSH}$ [recall equation (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. 1 (c). 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. 1 (d). 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. 1 (c) and proven in Supplementary Information Section III. However, the implementation of an exchange of qubits via $\mathcal{U}_{\rm\scriptscriptstyle SWAP}$ facilitates the application of the controlled phase gate $\mathcal{U}_{\rm\scriptscriptstyle{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 $\mathcal{U}_{\rm\scriptscriptstyle SWAP}$ with an arbitrary qubit of the opposite type, applies the CP gate, and maps the qubit back by another exchange.

III 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 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.

IV Acknowledgements

This research has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 681208). This research was supported in part by the National Science Foundation under Grant No. NSF PHY-1125915.