Arxiv:2006.13529V2 [Quant-Ph] 1 Jul 2020 Utrwith Ductor Es[E Fig
Total Page:16
File Type:pdf, Size:1020Kb
Memory-Critical Dynamical Buildup of Phonon-Dressed Majorana Fermions Oliver Kaestle,1, ∗ Ying Hu,2, 3 Alexander Carmele1 1Technische Universit¨at Berlin, Institut f¨ur Theoretische Physik, Nichtlineare Optik und Quantenelektronik, Hardenbergstrae 36, 10623 Berlin, Germany 2State Key Laboratory of Quantum Optics and Quantum Optics Devices, Institute of Laser Spectroscopy, Shanxi University, Taiyuan, Shanxi 030006, China 3Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China (Dated: July 28, 2021) We investigate the dynamical interplay between topological state of matter and a non-Markovian dissipation, which gives rise to a new and crucial time scale into the system dynamics due to its quan- tum memory. We specifically study a one-dimensional polaronic topological superconductor with phonon-dressed p-wave pairing, when a fast temperature increase in surrounding phonons induces an open-system dynamics. We show that when the memory depth increases, the Majorana edge dynam- ics transits from relaxing monotonically to a plateau of substantial value into a collapse-and-buildup behavior, even when the polaron Hamiltonian is close to the topological phase boundary. Above a critical memory depth, the system can approach a new dressed state of topological superconductor in dynamical equilibrium with phonons, with nearly full buildup of Majorana correlation. Exploring topological properties out of equilibrium is stantial preservation of topological properties far from central in the effort to realize, probe and exploit topo- equilibrium, even when the polaron Hamiltonian is close logical states of matter in the lab [1–11]. A paradig- to the topological phase boundary [see Fig. 2 (b)]. De- matic scenario is where the topological system is coupled pending on the memory depth (i.e., the characteristic to a Markovian bath, inducing open-dissipative dynam- time scale of the quantum memory), the Majorana edge ics that is described by a Lindblad-form master equa- dynamics can monotonically relax to a plateau, or it un- tion for the time-evolved reduced system density oper- dergoes a collapse-and-buildup relaxation [see Fig. 2 (c)]. ator [12–22]. Yet, solid-state realizations of topological Remarkably, when the memory depth increases above a matter often rely on nanotechnological design strategies critical value, the Majorana correlation can nearly fully that restructure the environment by fine-tuning the cor- build up, corresponding to a new polaronic state of topo- responding frequency-dependent density of states. In this logical superconductor in dynamical equilibrium with case, the Markovian approximation and thus the Lind- phonons. This topological polaronic steady-state goes blad formalism usually fails, e.g. in condensed mat- beyond typical frameworks where phonons act as a per- ter with nanostructured acoustic environments. Such turbation in the static and weak coupling limit. non-Markovian situations can also occur when ultracold Concretely, we consider the paradigmatic Kitaev p- atoms are immersed into Bose gases [23–25]. Compared wave superconductor [28], with a superohmic coupling to Markovian scenarios, key differences arise from the to a 3D structured phonon reservoir. The total Hamil- presence of quantum memory effects in non-Markovian tonian is denoted by H0 = Hk + Hb (~ ≡ 1). The Ki- processes: The information is lost from the system to N−1 † taev Hamiltonian Hk = l=1 [(−Jcl cl+1 + ∆clcl+1)+ the environment but flows back at a later time [26, 27]. N † P † H.c.] − µ l=1 cl cl describes spinless fermions cl, cl on This generates a new and critical time scale into the sys- a chain ofPN sites l, with a nearest-neighbor tunneling arXiv:2006.13529v2 [quant-ph] 1 Jul 2020 tem dynamics that is strictly absent in a Markovian con- text, and raises the challenge as to what are the unique dynamical consequences of the interplay between topo- logical state of matter and non-Markovian dissipation. Here, we demonstrate that the quantum memory from a non-Markovian parity-preserving interaction of a topo- logical p-wave supercondutor with surrounding phonons can give rise to intriguing edge mode relaxation dy- Figure 1. A polaronic Kitaev chain, with phonon-dressed namics with no Markovian counterpart. Our study is spinless fermions, exhibits a renormalized p-wave pairing based on the polaron master equation, describing open- hBi∆ at temperature T [see Eq. (1)]. In the topological dissipative dynamics of a polaronic topological supercon- ground state, two unpaired Majorana edge modes γL, γR ductor with phonon-renormalized Hamiltonian parame- emerge. The coupling gk to the structured phonon bath fea- ters [see Fig. 1]. In contrast to Markovian decoherence tures mode-dependence with a spectral width σ. The right that typically destroys topological features for long times, panel illustrates a Gaussian profile of gk in momentum space we show that a finite quantum memory allows for sub- (see text), respectively, for σ = 0.2 (green) and σ = 0.6 (red). 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 † † † iRsin(ω τ)}, 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.