arXiv:2006.13529v2 [quant-ph] 1 Jul 2020 utrwith ductor es[e Fig. [see ters ae nteplrnmse qain eciigopen- is describing dy- study a equation, of relaxation Our master dynamics dissipative polaron mode the counterpart. edge on Markovian based intriguing no to with namics rise give can o-akva aiypeevn neato fatopo- a logical of interaction parity-preserving non-Markovian a dissipation. non-Markovian and topo- unique matter between the of interplay are state what the logical to of as consequences challenge con- dynamical Markovian the a raises sys- in and the absent into text, strictly scale is time that to critical dynamics and tem system new a the a at generates from back This flows lost but is environment non-Markovian information the the in The effects from memory arise processes: quantum differences of key presence scenarios, Markovian to oper- density equa- system master reduced Lindblad-form [ time-evolved ator a the by for described tion dynam- is open-dissipative that inducing ics coupled bath, is Markovian system a topological to the where is scenario matic eso htafiieqatmmmr losfrsub- for allows memory quantum finite a that times, show long we for features topological destroys typically that [ gases Bose into ultracold mat- immersed Such when are condensed occur atoms also in environments. can acoustic situations non-Markovian nanostructured e.g. with Lind- fails, the ter usually thus formalism and approximation blad Markovian the this In case, cor- states. the of density fine-tuning frequency-dependent by responding environment strategies the design restructure nanotechnological that on rely often matter oia ttso atri h a [ lab topo- the exploit in and matter probe of realize, states to effort logical the in central ee edmntaeta h unu eoyfrom memory quantum the that demonstrate we Here, xlrn oooia rprisoto qiiru is equilibrium of out properties topological Exploring 12 p – wv uecnuo ihsronigphonons surrounding with supercondutor -wave 22 eoyCiia yaia ulu fPoo-rse Majo Phonon-Dressed of Buildup Dynamical Memory-Critical .Yt oi-tt elztoso topological of realizations solid-state Yet, ]. ndnmcleulbimwt hnn,wt eryfl bui full nearly dressed with new phonons, t a with the approach equilibrium to can dynamical close system in the is depth, Hamiltonian memory polaron sub critical of the plateau when a even to inc behavior, monotonically depth relaxing memory from the transits when ics that show We dynamics. open-system u eoy eseiclysuyaoedmninlpolaron one-dimensional a study specifically scale phonon-dressed time We crucial and memory. new a tum to rise gives which dissipation, phonon-renormalized 1 eivsiaetednmclitrlybtentopologica between interplay dynamical the investigate We .I otatt akva decoherence Markovian to contrast In ]. ihlnaeOtkudQatnlkrnk Hardenbergstr Quantenelektronik, und Optik Nichtlineare nttt fLsrSetocp,Sax nvriy Taiyu University, Shanxi Spectroscopy, Laser of Institute 2 tt e aoaoyo unu pisadQatmOtc De Optics Quantum and Optics Quantum of Laboratory Key State polaronic p 1 wv arn,we attmeaueices nsurround in increase temperature fast a when pairing, -wave ehiceUiesta eln ntttf¨ur Theoretisc Institut Universit¨at Berlin, Technische lvrKaestle, Oliver 3 hniUiest,Tiun hni000,China 030006, Shanxi Taiyuan, University, Shanxi olbrtv noainCne fEteeOptics, Extreme of Center Innovation Collaborative aitna parame- Hamiltonian oooia supercon- topological 1 ae time later 23 – 11 – 25 .Aparadig- A ]. .Compared ]. 1, Dtd uy2,2021) 28, July (Dated: [ ∗ 26 igHu, Ying , 27 ]. rudsae w nardMjrn demodes edge Majorana unpaired two state, ground Hamil- total The by reservoir. denoted phonon is structured tonian 3D a to un- Fig. it [see or relaxation plateau, collapse-and-buildup edge a a to dergoes Majorana characteristic relax the the monotonically memory), can quantum (i.e., dynamics the depth of scale memory time the on pending iue1 oaoi iavcan ihphonon-dressed renormalized with a chain, Kitaev exhibits polaronic fermions, A spinless 1. Figure aesprodco [ superconductor wave limit. coupling weak per- and a static as the act goes in phonons steady-state turbation where polaronic frameworks with typical topological beyond equilibrium This dynamical phonons. in superconductor fully logical a nearly to can corresponding correlation up, a build Majorana above the increases value, depth critical memory the when Remarkably, h ae lutae asinpol of profile Gaussian a illustrates panel otetplgclpaebudr seFig. [see boundary phase close topological is from Hamiltonian the polaron to far the properties when even topological equilibrium, of preservation stantial setx) epciey for respectively, text), (see mre h coupling The emerge. ue oedpnec ihaseta width spectral a with mode-dependence tures avHamiltonian taev H hi of chain a ,3 2, B . oceey ecnie h aaimtcKitaev paradigmatic the consider we Concretely, c i ttemperature at ∆ . lxne Carmele Alexander ] − notesse yaisdet t quan- its to due dynamics system the into tnilvleit collapse-and-buildup a into value stantial µ dpo aoaacorrelation. Majorana of ldup tt fmte n non-Markovian a and matter of state l P plgclpaebudr.Aoea Above boundary. phase opological e3,163Bri,Germany Berlin, 10623 36, ae N tt ftplgclsuperconductor topological of state ctplgclsprodco with superconductor topological ic l N =1 n hni000,China 030006, Shanxi an, ess h aoaaeg dynam- edge Majorana the reases, sites c l † c l H ePhysik, he l ecie pnesfermions spinless describes ihanaetniho tunneling nearest-neighbor a with , k g T H k = n hnn nue an induces phonons ing 28 1 0 otesrcue hnnbt fea- bath phonon structured the to seE.( Eq. [see σ P ,wt ueomccoupling superohmic a with ], = 0 = N l vices, =1 H e oaoi state polaronic new − . k 1 gen and (green) 2 [( + − aaFermions rana 1 H ] ntetopological the In )]. Jc g k b l † nmmnu space momentum in ( c ~ l +1 ≡ p ∆ + wv pairing -wave σ σ ) h Ki- The 1). 2 h right The . 0 = b] De- (b)]. c c l ftopo- of l . c , (red). 6 γ l 2 +1 L c l † , (c)]. + ) γ on p R - 2 amplitude J ∈ R, pairing amplitude ∆ ∈ R, and chem- in a polaronic Kitaev chain featuring phonon-dressed ical potential µ. When |µ| < 2J and ∆ 6= 0, the su- p-wave pairing, with phonon-induced quantum fluctua- perconductor is in the topological regime featuring un- tions. Note that the polaron transformation also renders † † paired Majorana edge modes γL = γL and γR = γR an energy renormalization which has been canceled by a at two ends [29], which exhibit a nonlocal correlation counter term in Eq. (1)[38, 43]. This is justified given the θ = −ihγLγRi = ±1 that corresponds to the fermionic coupled system is initially in equilibrium in the present parity distinguishing the two degenerate ground states. study. Before tracing out the phononic degrees of free- The entire chain is coupled to a 3D structured phonon dom, we rewrite Eq. (1) as Hp = Hp,s + Hp,I + Hp,b with † reservoir with a parity-preserving interaction, described Hp,b = dk ωkrkrk for the reservoir. To recover the 3k † N † † by Hb = d ωkrkrk + l=1 gkcl cl(rk + rk) [30, 31]. bare KitaevR Hamiltonian dynamics for the limiting case R  (†)  Here, operator rk annihilatesP (creates) phonons with gk → 0, we introduce a Franck-Condon renormalization momentum k and frequencies ωk = csk, where cs is the of Hp,I satisfying TrB{[Hp,I,ρ(t)]} = 0, with ρ(t) denot- sound velocity of the environment. We choose a generic ing the total density operator. The renormalized system superohmic fermion-phonon coupling gk which features Hamiltonian Hp,s is given by a frequency dependence modeled as a Gaussian function, 2 2 2 N−1 N gk = fph k/(V σ )exp −k /σ with a width σ and a † † † † dimensionlessp amplitude f . Such
fermion-phonon in- Hp,s= −Jcl cl+1 +∆hBicl+1cl +H.c. −µ cl cl, (2) ph X h i X teractions can represent the coupling of electrons with l=1 l=1 acoustic phonons in relevant solid-state experiments and where the pairing renormalization factor hBi is given protocols [32–34], and can arise in ultracold atoms from explicitly below. The system-reservoir interaction in coupling a fermionic lattice to a Bose-Einstein conden- N−1 −2(R−R†) the polaron picture reads Hp,I = ∆ [(e − sate inducing p-wave superconductivity [35–37]. We con- l=1 hBi)c† c† + H.c.]. Crucially, in theP following we treat sider the topological superconductor and bath are ini- l+1 l H perturbatively in second-order Born theory, as we tially in equilibrium at low temperatures, before a fast p,I are interested in phonon equilibration time scales much increase in the bath temperature induces an open-system faster than the system dynamics, so that dynamical de- dynamics. coupling effects cannot occur [38, 40]. We then derive Polaron master equation. The interactions with the the polaron master equation for the reduced system den- structured reservoir result in non-Markovian dynamics, sity matrix ρ (t) of the polaron chain, obtaining (Suppl. whose description - particularly on long time scales - S Mat.) is a challenge due to unfavorable scaling of the Hilbert space dimension. At the heart of our following solu- ∂tρS(t)= −i [Hp,s,ρS(t)] tion lies the polaron representation of the coupled sys- t tem [see Fig. 1]: We derive a master equation in second- 2 − hBi dτ (cosh [φ(τ)] − 1) [Xa,Xa(−τ)ρS(t)] order perturbation theory of the dressed-state system- Z0 n reservoir Hamiltonian tracing out phonons, whilst retain- − sinh [φ(τ)] [Xb,Xb(−τ)ρS (t)] + H.c. . (3) ing the coherent process (i.e., higher order contributions) o from the fermion-phonon interaction through phonon- N−1 † † renormalized Hamiltonian parameters [38–42], thus effi- Here Xa = −J l=1 (cl cl+1 + cl+1cl) and N−1 † † P ciently accounting for the non-Markovian character of the Xb = J l=1 (cl cl+1 − cl+1cl) denote collective dynamics in the long-time limit not captured in typical system operators,P whose dynamics obeys a time-reversed second-order Born approximation of the bare interaction unitary evolution governed by the renormalized Hamil- −iHp,sτ iHp,sτ Hamiltonian Hb [30]. tonian Hp,s, Xa,b(−τ) ≡ e Xa,be including Defining collective bosonic operators R† = the full density matrix (Suppl. Mat.). The φ(τ) † dk (gk/ωk)rk, we apply a polaron transforma- represents the phonon correlation function, φ(τ) = 3 2 R N † † d k |2gk(σ)/ωk| {coth [~ωk/(2kBT )] cos (ωkτ) − tion Up = exp[ l=1 cl cl(R − R)], which results † † † Ri sin(ω τ)}, with k the Boltzmann constant. The in c → e−(R−RP)c that describes phonon dress- k B l l renormalization factor hBi in Eq. (2) is determined by ing of fermions. The transformed total Hamiltonian the initial phonon correlation, hBi = exp[−φ(0)/2], with H ≡ U H U −1 is derived as p p 0 p temperature dependence. N−1 Equation (3) provides the key equation for our study. † −2(R−R†) † † Hp = − Jcl cl+1 + ∆e cl cl + H.c. It features a phonon-renormalized Hamiltonian H de- h +1 i p,s Xl=1 scribing polarons which exhibit temperature-dependent N pairing, and a non-Markovian dissipative term in the † 3 † − µ c c + d k ω r rk. (1) l l Z k k form of a memory kernel which involves both reservoir Xl=1 and system correlators φ(τ) and Xa,b(−τ), respectively, Thus the considered fermion-phonon interaction results with a finite memory depth determined by the lifetime 3

(a) (b) (c) 10 1 1 Re[Φ(τ)] Im[Φ(τ)] Abs[Φ(τ)] 5 0.5 0.5

(t) (t) 1 θ θ ) ∞

0 0 0 (t

σ = 0.4 θ

0.2 σ = 0.2 0.2 σ 0.8 -5 -0.5 -0.5 0 25 50 0 0.5 1 0 1 2 Jτ 1000 Jt 1000 Jt

Figure 2. Non-Markovian dynamics of the polaronic topological superconductor. (a) Phonon correlation function φ(τ) for bandwidth σ = 0.2 (green) and σ = 0.6 (red) of the fermion-phonon coupling gk. (b) Comparisons of Majorana correlation θ(t) calculated using, respectively, time-independent Lindblad-type master equation for dephasing (solid black), Eq. (3) in the Markovian limit (blue), Eq. (3) with full account of memory (orange). The dashed black line shows asymptotic θ(t) in a coherent quench scenario ∆ → ∆hBi. (c) Non-Markovian dynamics of θ(t) for various bandwidths σ of phonon coupling. The corresponding steady-state value θ(t∞) is shown as a function of σ in the inset. In all plots, the amplitude of gk is taken as −1 fph = 0.1, and phonon modes within k ∈ [0.0, 4.0] nm are considered. Due to the numerically very expensive size of the density matrix and its memory kernel, computations are performed for N = 4 sites. of φ(τ). The memory depth critically depends on the system and reservoir unfolds. The orange line in Fig. 2 bandwidth σ of fermion-phonon coupling gk [see Fig. 2 (b) shows the non-Markovian dynamics for σ =0.6 corre- (a)]: A larger σ results in a faster decay of φ(t) and sponding to hBi =0.07: The Majorana edge correlation hence a smaller memory size, but it also leads to a smaller still decays in an oscillatory manner [orange line], but it φ(0) and thus a larger hBi in both Hp,s and the prefac- relaxes to a substantial value as opposed to the Marko- tor of the memory kernel. In the limit when the system vian limit [blue line]. Considering the small system size evolves much slower than the phonons, one can approxi- in our computation, we have verified that such asymp- mate Xa,b(−τ) ≈ Xa,b, andEq.(3) transits into a Marko- totic non-local correlation is genuinely of topological ori- vian type of master equation with a time-dependent de- gin, rather than phonon-mediated long range correlations phasing rate γ(t). (Suppl. Mat.). Below we investigate the Majorana edge correlation The long lived and substantial Majorana correlation in θ(t) = −iTr [ρs(t)γLγR] based on Eq. (3), starting from Fig. 2 (b) is quite remarkable, given that Hp,s is near the a factorizing system-bath state. For concreteness, we as- topological phase boundary due to a significantly sup- sume the system is initially at zero temperature in the pressed renormalized pairing ∆ hBi ≪ ∆. Indeed, with- ground state of a dressed Kitaev Hamiltonian with ∆ = J out the dissipation in Eq. (3), the dynamics formally re- and µ = 0, exhibiting even parity θ(0) = 1. Then, a fast duces to that of a coherent quench in the pairing from increase of temperature to T = 4 K results in ∆ → ∆hBi ∆ to ∆ hBi. There, the Majorana correlation would ap- of the dressed Hamiltonian, thus inducing the dynamics proach an asymptotic value determined by the overlap of the polaron chain for times t> 0. of the edge mode wave functions for the pre- and post- Markovian limit. In the Markovian limit Xa,b(−τ) ≈ quench topological Hamiltonians [18], which is small if Xa,b(0) of Eq. (3), the edge correlation θ(t) rapidly de- the post-quench Hamiltonian is close to the phase bound- cays to a very small but finite value [blue line in Fig. 2 ary [see dashed black line in Fig. 2 (b)]. This differs sig- (b)]. This is expected because dressing fermions with fast nificantly from the non-Markovian behavior in Fig. 2 (b) phonons induces dephasing, as has been familiar from and underscores the essential role of the memory effect. Markovian decoherence based on a time-independent A unique feature of the memory effect is that it, Lindblad operator [see e.g., Ref. [44] and black line in because of the dependence on both φ(τ) and the Fig. 2 (b)]. Still, the non-Markovian character of the reversed dynamics of system correlations Xa,b(−τ), structured bath leads to a small residue correlation. simultaneously introduces decoherence and backflow Memory: loss vs. rephasing of topological properties. of coherence. The dynamical consequence of these However, the picture drastically changes when taking two competing processes can be intuitively under- into account the full memory including the system’s past stood as follows: The dressed Kitaev wire initially in Xa,b(−τ), as a highly non-Markovian interplay between its ground state is perturbed by a temperature in- 4 crease to T , resulting in a renormalization of the po- 1 laron chain towards the phase boundary via φ(0) = k 2 ~ d |2gk(σ)/ωk| coth ( ωk/(2kBT )), which generates 0.5 significantR bulk excitations and populates the Majorana (t) edge mode, changing the parity of Majorana states. θ Combined with phonon-assisted dephasing, this leads to 0 strong decoherence in the polaronic Kitaev wire. On the U = 0.002 U = -0.002 other hand, the reversed dynamics of Xa,b(−τ) acts to U = 0.000 U = -0.004 -0.5 reinstate coherence of the p-wave pairing that is the key 0 0.5 1 ingredient for a topological wire. Such a rephasing ef- 1000 Jt fect is marginal at times smaller than the characteristic time of the memory, so an irreversible loss of parity in- Figure 3. Non-Markovian Majorana dynamics for renormal- ized Hamiltonian Hp,s in the presence of weak p-wave interac- formation dominates the short-time dynamics. Yet φ(τ) N−1 † † tion, Hint = U (c c − 1/2)(c c +1 − 1/2). Other pa- proceeds to decay over time, as its sine and cosine func- Pl=1 l l l+1 l rameters, as well as initial conditions, are the same as Fig. 2. tions tend to cancel each other. Once φ(τ) approaches zero at large times, the memory reaches its full depth and the rephasing of topological properties grows due the system asymptotically approaches a new polaronic to Xa,b(−τ), giving considerable Majorana correlation in Fig. 2 (b). steady-state, in dynamical equilibrium with phonons at Critical memory depth. We find the edge dynamics can T = 4 K, which can remarkably exhibit θ ≈ 1. The crit- exhibit distinct relaxation behavior depending crucially ical value of σ in our case is between 0.21 and 0.20 cor- on the memory depth, tunable through the bandwidth σ responding to hBi = 0.01, but it is model specific. We of fermion-phonon coupling. Fig. 2 (c) presents the non- note that the fermion-phonon coupling bandwidth can Markovian dynamics of θ(t) for various σ. Compared to be controlled, such as in solid state setups by nanotech- σ = 0.6 [orange line of Fig. 2 (b)], an initial decrease nological design, e.g. alloys, impurities and confinement of σ results in a steeper monotonic decay of θ(t) and a potentials [45]. smaller asymptotic value [blue line in Fig. 2 (c)]. How- Concluding discussions. The central results of our ever, when σ decreases further, the monotonic relaxation work shown in Fig. 2 are found to be robust for initially transits into a non-monotonic one: While the short-time non-ideal Kitaev chains and other forms of superohmic decoherence is accelerated, a buildup of edge correlation coupling (Suppl. Mat.). Moreover, we find that the pres- nonetheless occurs at some large times [purple line]. Such ence of a weak attractive p-wave interaction can strongly buildup becomes stronger with decreased σ, approach- suppress the phonon-induced short-time parity loss, thus ing an asymptotic value larger than the σ = 0.6 case augmenting the buildup of Majorana correlation for sub- [turquoise and green lines]. Strikingly, once σ surpasses critical memory depth. In Fig. 3, we calculate the non- a critical value, the asymptotic Majorana correlation ap- Markovian evolution of θ(t) by including a weak p-wave N−1 † † proaches θ(t∞) → 1 [orange and dotted red lines]. Fig. 2 interaction, Hint = U l=1 (cl cl − 1/2)(cl+1cl+1 − 1/2) (c) summarizes the non-monotonic variation of θ(t∞): with interaction strengthP |U| ≪ 1, in Hp,s of Eq. (3) for When σ decreases from a large value, θ(t∞) first de- σ =0.6. Compared to the U = 0 case [blue line], adding a creases to a minimum and then increases toward unity. weak attractive interaction U < 0 damps out oscillations Insights into above intriguing phenomena can be ob- and, in particular, strongly mitigates the transient decay. tained from the fact that reducing σ leads to an increased Depending on |U|, the asymptotic Majorana correlation memory size at the cost of a smaller hBi = exp[−φ(0)/2] can be significantly increased [orange and red lines]. An [see Fig. 2 (a)]; the former enhances the time scale of intuitive understanding can be obtained by noting that the time-reversed evolution of Xa,b(−τ) and hence the the attractive interaction is energetically favorable for co- rephasing of pairing, whereas the latter further sup- herent formation of superconductive pairing, which pro- presses the superconducting gap, rendering Hp,s closer vides a mechanism to counteract aforementioned phonon- to the phase boundary as well as weakening the mem- induced dephasing. This is consistent with the observa- ory strength. When σ is initially decreased from 0.6, tion that for U > 0, θ(t) significantly declines from the the latter effect dominates, aggravating the decay. With U = 0 case at long times [purple line], as repulsive inter- further reduction of σ, however, the former rephasing actions energetically suppress pairing. effect grows, allowing phonons and fermions to synchro- Summarizing, we have demonstrated memory-critical nize and hence inducing backflow of parity information. edge dynamics in a topological superconductor with non- Consequently, a new dressed state manages to emerge, Markovian interaction with phonons. We show this in- with buildup of Majorana correlation starting to domi- triguing phenomenon uniquely arises from the interplay nate over decoherence. Importantly, the existence of a between the phonon-renormalized topological Hamilto- critical σ indicates a critical memory depth, above which nian and the quantum memory effect that simultaneously 5 induces dephasing and information backflow. Our analy- Phys. Rev. Lett. 102, 065703 (2009). sis is based on the Kitaev chain, but we expect the essen- [13] S. Diehl, E. Rico, M. A. Baranov, and P. Zoller, tial physics to occur for a wide class of topological materi- Nature Physics 7, 971 (2011). 85 als coupled to a superohmic reservoir. These discussions [14] D. Rainis and D. Loss, Phys. Rev. B , 174533 (2012); M. J. Schmidt, D. Rainis, and D. Loss, are different from recent work [46–49], e.g., where an im- Phys. Rev. B 86, 085414 (2012); M. Ippoliti, purity or a qubit acts as a non-Markovian quantum probe M. Rizzi, V. Giovannetti, and L. Mazza, of a topological reservoir. Our result is relevant to ongo- Phys. Rev. A 93, 062325 (2016). ing efforts aimed at realizing topological computations in [15] L. Mazza, M. Rizzi, M. D. Lukin, and J. I. Cirac, experimentally realistic conditions where non-Markovian Phys. Rev. B 88, 205142 (2013). effects are inevitable. It further opens an appealing new [16] C.-E. Bardyn, M. A. Baranov, C. V. Kraus, ˙ prospect as to the explorations and control of memory- E. Rico, A. Imamo˘glu, P. Zoller, and S. Diehl, 15 dependent topological phenomena. New Journal of Physics , 085001 (2013). [17] F. L. Pedrocchi and D. P. DiVincenzo, O.K. and A.C. gratefully acknowledge support from Phys. Rev. Lett. 115, 120402 (2015). the Deutsche Forschungsgemeinschaft (DFG) through [18] Y. Hu, Z. Cai, M. A. Baranov, and P. Zoller, SFB 910 project B1 (project number 163436311). Y.H. Phys. Rev. B 92, 165118 (2015). acknowledges National Natural Science Foundation of [19] J. C. Budich, P. Zoller, and S. Diehl, Phys. Rev. A 91, 042117 (2015). China (Grant No. 11874038). [20] V. V. Albert, B. Bradlyn, M. Fraas, and L. Jiang, Phys. Rev. X 6, 041031 (2016). [21] Z. Gong, S. Higashikawa, and M. Ueda, Phys. Rev. Lett. 118, 200401 (2017). [22] S. Lieu, M. McGinley, and N. R. Cooper, ∗ [email protected] Phys. Rev. Lett. 124, 040401 (2020). [1] N. Lindner, G. Refael, and V. Galitski, [23] T. H. Johnson, M. Bruderer, Y. Cai, Nature Physics 7, 490 (2011). S. R. Clark, W. Bao, and D. Jaksch, [2] T. Kitagawa, M. A. Broome, A. Fedrizzi, M. S. Rudner, EPL (Europhysics Letters) 98, 26001 (2012). E. Berg, I. Kassal, A. Aspuru-Guzik, E. Demler, and [24] S. Peotta, D. Rossini, M. Polini, F. Minardi, and A. G. White, Nature Communications 3, 882 (2012). R. Fazio, Phys. Rev. Lett. 110, 015302 (2013). [3] G. Jotzu, M. Messer, R. Desbuquois, M. Le- [25] F. Grusdt, K. Seetharam, Y. Shchadilova, and E. Dem- brat, T. Uehlinger, D. Greif, and T. Esslinger, ler, Phys. Rev. A 97, 033612 (2018). Nature 515, 237 (2014). [26] H.-P. Breuer, E.-M. Laine, J. Piilo, and B. Vacchini, [4] D. Aasen, M. Hell, R. V. Mishmash, A. Higgin- Rev. Mod. Phys. 88, 021002 (2016). botham, J. Danon, M. Leijnse, T. S. Jespersen, J. A. [27] I. de Vega and D. Alonso, Folk, C. M. Marcus, K. Flensberg, and J. Alicea, Rev. Mod. Phys. 89, 015001 (2017). Phys. Rev. X 6, 031016 (2016). [28] A. Y. Kitaev, Physics-Uspekhi 44, 131 (2001). [5] Y. Hu, P. Zoller, and J. C. Bu- [29] The Majorana edge modes are of the form: γL/R = dich, Phys. Rev. Lett. 117, 126803 (2016); Pj fL/R,j aj in the Majorana representation a2j−1 = J. C. Budich, Y. Hu, and P. Zoller, † † c + c and a2 = −i(c − c ) of the fermionic opera- Phys. Rev. Lett. 118, 105302 (2017). j j j j j tors, with fL/R,j being exponentially localized near the [6] N. Fl¨aschner, D. Vogel, M. Tarnowski, B. S. Rem, left (L) and right (R) edges. D.-S. L¨uhmann, M. Heyl, J. C. Budich, L. Mathey, 14 [30] G. D. Mahan, Many-particle physics (Springer Science & and K. Sengstock, Nature Physics , 265 (2018); Business Media, 2013); M. Kira and S. W. Koch, Semi- ¨ M. Tarnowski, F. N. Unal, N. Fl¨aschner, B. S. Rem, conductor quantum optics (Cambridge University Press, A. Eckardt, K. Sengstock, and C. Weitenberg, 2011); V. May and O. K¨uhn, Charge and energy trans- 10 Nature Communications , 1728 (2019). fer dynamics in molecular systems (John Wiley & Sons, 81 [7] M. Heyl, Reports on Progress in Physics , 054001 (2018). 2008); M. A. Stroscio and M. Dutta, Phonons in nanos- [8] W. Sun, C.-R. Yi, B.-Z. Wang, W.-W. Zhang, B. C. tructures (Cambridge University Press, 2001). Sanders, X.-T. Xu, Z.-Y. Wang, J. Schmiedmayer, [31] J. F¨orstner, C. Weber, J. Danckwerts, and Y. Deng, X.-J. Liu, S. Chen, and J.-W. Pan, A. Knorr, Phys. Rev. Lett. 91, 127401 (2003); 121 Phys. Rev. Lett. , 250403 (2018); B. Song, C. He, B. Krummheuer, V. M. Axt, and S. Niu, L. Zhang, Z. Ren, X.-J. Liu, and G.-B. Jo, T. Kuhn, Physical Review B 65, 195313 (2002); 15 Nature Physics , 911 (2019); L. Zhang, L. Zhang, A. Carmele, A. Knorr, and F. Milde, 63 S. Niu, and X.-J. Liu, Science Bulletin , 1385 (2018). New Journal of Physics 15, 105024 (2013); [9] Z. Shi, P. W. Brouwer, K. Flensberg, L. I. Glazman, and N. N´emet, S. Parkins, A. Knorr, and F. von Oppen, (2020), arXiv:2002.01539. A. Carmele, Phys. Rev. A 99, 053809 (2019); [10] J. W. McIver, B. Schulte, F.-U. Stein, T. Mat- D. Reiter, T. Kuhn, and V. M. Axt, suyama, G. Jotzu, G. Meier, and A. Cavalleri, Advances in Physics: X 4, 1655478 (2019). 16 Nature Physics , 38 (2020). [32] Y. Oreg, G. Refael, and F. von Oppen, [11] A. Elben, J. Yu, G. Zhu, M. Hafezi, Phys. Rev. Lett. 105, 177002 (2010); P. W. Brouwer, F. Pollmann, P. Zoller, and B. Vermersch, M. Duckheim, A. Romito, and F. von Oppen, 6 Science advances , 3666 (2020). Phys. Rev. Lett. 107, 196804 (2011); V. Mourik, [12] M. S. Rudner and L. S. Levitov, 6

K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. Bakkers, K. S¨a¨askilahti, Y. A. Kosevich, H. Han, D. Donadio, and and L. P. Kouwenhoven, Science 336, 1003 (2012); S. Volz, Phys. Rev. Lett. 117, 025503 (2016); A. A. Ba- L. P. Rokhinson, X. Liu, and J. K. Furdyna, landin, Journal of Nanoscience and Nanotechnology 5, Nature Physics 8, 795 (2012); A. Das, Y. Ronen, 1015 (2005); A. A. Balandin, E. P. Pokatilov, and D. L. Y. Most, Y. Oreg, M. Heiblum, and H. Shtrik- Nika, Journal of Nanoelectronics and Optoelectronics 2, man, Nature Physics 8, 887 (2012); S. Nadj-Perge, 140 (2007). I. K. Drozdov, J. Li, H. Chen, S. Jeon, J. Seo, [46] F. Grusdt, N. Y. Yao, D. Abanin, M. Fleischhauer, and A. H. MacDonald, B. A. Bernevig, and A. Yazdani, E. Demler, Nature Communications 7, 11994 (2016). Science 346, 602 (2014). [47] A. Camacho-Guardian, N. Goldman, P. Massignan, and [33] H. O. H. Churchill, V. Fatemi, K. Grove-Rasmussen, G. M. Bruun, Phys. Rev. B 99, 081105 (2019). M. T. Deng, P. Caroff, H. Q. Xu, and C. M. Mar- [48] G. L. Giorgi, S. Longhi, A. Cabot, and R. Zambrini, cus, Phys. Rev. B 87, 241401 (2013); A. D. K. Finck, Annalen der Physik 531, 1900307 (2019). D. J. Van Harlingen, P. K. Mohseni, K. Jung, and X. Li, [49] A. Ricottone, M. S. Rudner, and W. A. Coish, (2020), Phys. Rev. Lett. 110, 126406 (2013). arXiv:2003.00350. [34] V. Fatemi, S. Wu, Y. Cao, L. Bretheau, Q. D. Gib- son, K. Watanabe, T. Taniguchi, R. J. Cava, and P. Jarillo-Herrero, Science 362, 926 (2018); D. Laroche, D. Bouman, D. J. van Woerkom, A. Proutski, C. Murthy, D. I. Pikulin, C. Nayak, R. J. J. van Gulik, J. Ny- grd, P. Krogstrup, L. P. Kouwenhoven, and A. Geresdi, Nature Communications 10, 245 (2019); A. Fornieri, A. M. Whiticar, F. Setiawan, E. Portols, A. C. C. Drachmann, A. Keselman, S. Gronin, C. Thomas, T. Wang, R. Kallaher, G. C. Gardner, E. Berg, M. J. Manfra, A. Stern, C. M. Marcus, and F. Nichele, Nature 569, 89 (2019); H. Ren, F. Pientka, S. Hart, A. T. Pierce, M. Kosowsky, L. Lunczer, R. Schlereth, B. Scharf, E. M. Hankiewicz, L. W. Molenkamp, B. I. Halperin, and A. Yacoby, Nature 569, 93 (2019). [35] L. Jiang, T. Kitagawa, J. Alicea, A. R. Akhmerov, D. Pekker, G. Refael, J. I. Cirac, E. Demler, M. D. Lukin, and P. Zoller, Phys. Rev. Lett. 106, 220402 (2011). [36] S. Nascimb`ene, Journal of Physics B: Atomic, Molecular and Optical Physics 46, 134005 (2013). [37] Y. Hu and M. A. Baranov, Phys. Rev. A 92, 053615 (2015). [38] H. P. Breuer and F. Petruccione, The theory of open quantum systems (Oxford University Press, 2002). [39] I. Wilson-Rae and A. Imamo˘glu, Phys. Rev. B 65, 235311 (2002); S. Weiler, A. Ul- haq, S. M. Ulrich, D. Richter, M. Jetter, P. Michler, C. Roy, and S. Hughes, Phys. Rev. B 86, 241304 (2012); R. Manson, K. Roy-Choudhury, and S. Hughes, Phys. Rev. B 93, 155423 (2016); C. K. Lee, J. Moix, and J. Cao, The Journal of Chemical Physics 136, 204120 (2012). [40] D. P. S. McCutcheon and A. Nazir, New Journal of Physics 12, 113042 (2010); H.- T. Chang, P.-P. Zhang, and Y.-C. Cheng, The Journal of Chemical Physics 139, 224112 (2013). [41] A. Carmele and S. Reitzenstein, Nanophotonics 8, 655 (2019); M. Strauß, A. Carmele, J. Schleibner, M. Hohn, C. Schnei- der, S. H¨ofling, J. Wolters, and S. Reitzenstein, Phys. Rev. Lett. 122, 107401 (2019). [42] E. V. Denning, J. Iles-Smith, and J. Mork, Phys. Rev. B 100, 214306 (2019); E. V. Den- ning, J. Iles-Smith, N. Gregersen, and J. Mork, Opt. Mater. Express 10, 222 (2020). [43] A. Caldeira and A. Leggett, Physica A: Statistical Me- chanics and its Applications 121, 587 (1983). [44] A. Carmele, M. Heyl, C. Kraus, and M. Dalmonte, Phys. Rev. B 92, 195107 (2015). [45] W. Kim, J. Mater. Chem. C 3, 10336 (2015); S. Xiong,