arXiv:2004.13055v2 [quant-ph] 15 May 2020 loeteeyrbs ihrsett atcels,re- loss, particle to any if respect even entangled with maining robust extremely also esbsn n multipartite and bosons dispersion- with less in- excitation nonlocal spinless-fermion a a of describing teraction model lattice (1D) dimensional [22– systems 28]. various in preparation their motivates which 19–21], to [5, themselves protocols lend -information They in applications [18]. particle their about mation omdit ahohrvalcloeain n classical and operations [17]. local (LOCC-inequivalence) via communication other each into formed ui ttsaemxmlyetnldGreenberger- maximally-entangled and [15] are (GHZ) Horne-Zeilinger states photons and prepara- SC [14]. of accurate alike states them, allow quantum various among 13], of tion [12, and, systems [8–11] continu-circuit-QED circuits su- their (SC) properties to physical perconducting owing and diverse scalability particular, utilizing improving In ously by years progress [5–7]. recent tantalizing platforms context, in this made In was 4]. [3, quantum technologies next-generation of development the for requisite en nter“on tts npriua,i sknown is it ones particular, In remaining the states. a state, that “down” “up” their its in in being qubit one exactly qubit ne osdrto a eepotdfrafs,deter- model fast, the a of for property of exploited this preparation be how ministic can here excita- consideration shown bare is under a is It of state tion. its Bloch strength zero-quasimomentum coupling the e-b effective un- value an the critical a has of below – Namely, types state. cou- ground breathing-mode conventional (e-b) and regimes excitation-boson relevant Peierls the includes of of plings which one – in model that this notion the is parture iuao fti oe [29]. model analog this an of , simulator and qubits SC inductively-coupled hsLte salse oncinbtenaone- a between connection a establishes Letter This h otpoietcasso nage many- entangled of classes prominent most The pre- a is 2] [1, engineering quantum-state Sophisticated h tt-rprto rtclpooe ee based here, proposed protocol state-preparation The aeectto rudsaeo pnesfrin–boson – spinless-fermion a of state ground Bare-excitation W W tt steeulsproiino ttswith states of superposition equal the is state tt n t H onepr antb trans- be cannot counterpart GHZ its and state ntttfu neadePyi,TcnclUiest fD of University Technical Physik, f¨ur Angewandte Institut rpryo h atc oe ne osdrto a eexp be can consideration multipartite excitati under of undressed) model tion completely lattice (i.e., the of bare property a of state ground state exact Bloch the strength Namel t coupling mechanisms. excitation-boson coupled nonlocal-coupling (excitation) different fermion two through spinless itinerant single a uisadmcoaeresonators. microwave and qubits niern na ra fsprodcigqbt n reson and qubits superconducting of array an in engineering hsLte naesa neetn rpryo one-dimen a of property interesting an unravels Letter This N -qubit W N W − W tts t on fde- of point Its states. W ate oeteinfor- the lose parties 2 ttsi edl elzbesse finductively-coupl of system realizable readily a in states ttsi naryof array an in states tts[6.An [16]. states W ldmrM Stojanovi´c M. Vladimir ttsare states Dtd a 9 2020) 19, May (Dated: N - fismi xeietlko netra cflux. dc external an values – of knob window experimental state large main ground parametrically its the a of is particularly in state system protocol target the its this of that makes fact number What the the is on all. robust depend at not do qubits they with of here unfavorably size, scale system times the preparation control, typical quantum-state in where situation of the times Unlike coherence qubits. achievable SC currently the than shorter multipartite obtain W to one allows pumping, microwave on eatn ato t oa aitna nldsteexci- terms: the free-boson includes and Hamiltonian kinetic-energy- total tation its nonin- of The part mechanisms. through teracting coupling bosons nonlocal different dispersionless two ex- with spinless-fermion interacting single citation a describes consideration der Here ( aest,while site, same h neatn eb ati ie by given is part (e-b) interacting The incpue h niymti opigo h excita- the equa- site of last at coupling the density antisymmetric tion of the (rhs) captures right-hand-side tion the on term first edneo h ffcieectto-opn amplitude de- excitation-hopping linear sites the effective between for the accounts of term second pendence The [30]. pling) sites neighboring the on where lcmns(eel-yeculn)[31–33]. coupling) (Peierls-type placements eei oetmsaeform momentum-space generic n H oe n t rudstate ground its and Model h opigHamiltonian coupling The 1 = ttswti iefae he reso magnitudes of orders three frames time within states e-b H eodmninl(iprines bosons (dispersionless) zero-dimensional o n ti eosrtdhr o hslast this how here demonstrated is It on. c fti oe stezero-quasimomentum the is model this of 0 rsat -48 amtd,Germany Darmstadt, D-64289 armstadt, N , . . . , H g n † ,blwaciia au fteeffective the of value critical a below y, = = ( + otdfrafs,dtriitcprepara- deterministic fast, a for loited e-b stedmninesebculn tegh The strength. coupling e-b dimensionless the is ( ∗ c − n g inlltiemdlta describes that model lattice sional = t c rae dsry)a xiaina site at excitation an (destroys) creates ) ~ e n † ω ), +1 X n √ b n 1 b N X and c n † t ( n n e c ( X H + n † k,q n steectto opn amplitude. hopping excitation the is b h +1 n n c ihtelclbsndisplacements boson local the with oo ihfrequency with boson a ) n † γ c ntersetv oo dis- boson respective the on 1 + . c e-b n c n n . oe and model )( H + ( dsuperconducting ed ± ( b ,q k, b n † betigmd-yecou- (breathing-mode-type 1 n † –Te1 atc oe un- model lattice 1D The .– − +1 . H 1 c ) . + e-b + ) c + k † b + b n ators q a ercs nthe in recast be can ~ n − c +1 ω k 1 b ( − b − X W − † n b q b n † n † -state + +1 b n † − b b − q n b ω ) n . b b . ) n i tthe at +1 , ) (1) (3) (2) n 2

Its corresponding e-b vertex function depends on both units of Φ0/2π, where Φ0 hc/(2e) is the flux quantum]. the excitation- and boson quasimomenta (k and q, re- In particular, the upper loop≡ is subject to ac-driving with u spectively, here expressed in units of the inverse lattice the flux π cos(ω0t). The other contribution to φn origi- period) and is given by nates from the modes of resonators n and n + 1 and is † † given by φn,res = δθ[(bn+1 + bn+1) (bn + bn)], where ~ 1/2 − γe-b(k, q)=2ig ωb [ sin k + sin q sin(k + q)] . (4) δθ = [2eAeff/(~d0c)] (~ω /C0) , with Aeff being the − × c effective coupling area, C0 the capacitance, and The ground state of H = H0 + He-b undergoes a sharp d0 the effective spacing in the resonator [37]. Therefore, level-crossing transition [34] at a critical value λc e-b ∼ unlike the much more common capacitive coupling [38], 1 [cf. Fig. 2] of the effective coupling strength λe-b the qubit-resonator coupling in the system at hand is in- 2 ~ c ≡ 2g ωb/te. For λe-b < λe-b the ground state is the K =0 l ductive [11]. Similarly, φn includes an ac contribution, eigenvalue of the total quasimomentum Ktot = † † c given by (π/2) cos(ω0t), and a dc part φdc, the main kc ck + qb bq. For λe-b λ , on the other hand, − Pk k Pq q ≥ e-b experimental knob in this system. the ground state is twofold-degenerate and corresponds J The Josephson energy of Bn is given by Hn = to K = Kgs (Kgs = 0). 3 ± 6 Ei cos ϕi , with ϕi being the respective phase A ground state with K = 0 is by no means unusual – − Pi=1 J n n drops on the three Josephson junctions within Bn and in fact, an overwhelming majority of coupled e-b mod- i 1 2 EJ their energies, here chosen such that EJ = EJ EJ els have such ground states. Yet, the model at hand 3 ≡ and EJ = EJb = EJ . Using the flux-quantization has the peculiar property that its K = 0 ground state 6 J c rules [9], the total Josephson energy n Hn can be ex- for λe-b < λe-b is the k = 0 bare-excitation Bloch state P † pressed in terms of the gauge-invariant phase variables Ψk=0 c 0 e 0 b, where 0 e and 0 b are the exci- k=0 ϕn of SC islands of different qubits. The latter enter |tationi≡ and boson| i vacuum⊗| i states.| i In what| i follows, it will this energy through terms of the type cos(ϕn ϕn+1), first be demonstrated explicitly that Ψ =0 is an exact − | k i which in the regime of interest for transmons/gatemons eigenstate of H for an arbitrary value of λe-b. It will sub- s s (EJ EC ) can be recast (up to an additive constant) as sequently be shown numerically (see Fig. 2 below) that 2 ≫+ − − + z z 2 c δϕ0 σn σn+1 + σn σn+1 (σn + σn+1)/2 , where δϕ0 for λe-b < λe-b this state is the ground state of H. s s 1/2 −  ≡ (2EC /EJ ) [38]. Given that Ψk=0 is an eigenstate of H0, to prove that it is an eigenstate| ofi the total Hamiltonian H it suffices to show that it is also an eigenstate of He-b. Indeed, by acting with He-b [cf. Eq. (3)] on this state and making † use of the fact that c c 0 e δ 0 0 e, one obtains k 0 | i ≡ k, | i

1 † † He-b Ψ =0 = γe-b(k =0, q) c 0 e b 0 b . | k i √ X q| i ⊗ −q| i N q (5) Because here γe-b(k = 0, q) = 0 for an arbitrary q [cf. Eq. (4)], each term in the sum on the rhs of Eq. (5) van- ishes, implying that He-b Ψ =0 = 0. Therefore, Ψ =0 | k i | k i is an eigenstate of He-b (for an arbitrary λe-b), the corre- sponding eigenvalue being equal to zero. This concludes the proof that Ψ =0 is an exact eigenstate of H. | k i Qubit-resonator system.– The analog simulator of the model under consideration [see Fig. 1(a)] consists of SC qubits (Qn) with the energy splitting εz, microwave res- onators (Rn) with the photon frequency ωc, and coupler circuits (Bn) [35] which mediate both qubit-qubit and qubit-resonator interactions in this system. The simula- s s FIG. 1: (Color online)(a) Schematic of the qubit-resonator tor can be realized with transmons [12] (EJ /EC 100, s s ∼ system, whose n-th repeating unit (indicated by the dashed where EC and EJ are the single-qubit charging- and Josephson energies) or gatemons [36] (Es /Es 25). Its rectangle) comprises SC qubit Qn, resonator Rn, and coupler J C circuit Bn. The fluxes from the resonator modes n and n + 1 n-th repeating unit is described by the free Hamiltonian∼ thread the upper loops of the coupler circuit Bn, effectively 0 z ~ † Hn = (εz/2)σn + ωc bnbn, where the pseudospin-1/2 op- giving rise to an indirect inductive qubit-resonator coupling. erators σn represent qubit n and the bosonic operators (b) Pictorial illustration of the effective lattice model of the † (bn,bn) photons in the n-th resonator. system, with the excitation hopping amplitude t0(φdc) and The upper and lower loops of Bn are threaded by mag- each lattice site hosting dispersionless bosons with the fre- u l quency δω = ωc − ω0. netic fluxes φn and φn, respectively [both are expressed in 3

Further analysis is carried out in the rotating frame 0 of the drive, where δω ωc ω0 is the effective bo- -2 t ( φ ) son frequency and the Josephson-coupling≡ − term becomes 0 dc time-dependent. This time dependence can, however, u be disregarded due to its rapidly-oscillating character, in line with the rotating-wave approximation (RWA). ) / E -10

J e-b The remaining part of Hn can succinctly be written as c λ λ λ ¯ J J e-b = e-b Hn = n (φdc, φn,res) cos(ϕn ϕn+1), where (

−E − gs J = E (1 + cos φdc) E J1(π/2)φ res , (6) E En Jb − J n, -20 δω / 2π = 300 ΜΗz and Jm(x) are Bessel functions of the first kind whose presence in this expression stems from the use of the Jacobi-Anger expansion [39] in conjunction with the 1 2 3 4 5 λ RWA. In what follows, without significant loss of gen- e-b erality EJb is chosen to be given by 2EJ J0(π/2). The above expression for cos(ϕ ϕ +1) in terms of n − n FIG. 2: Ground-state energy of the system with δω/2π = the operators σn implies that the effective interaction 300 MHz as a function of the effective coupling strength λe-b. c between adjacent qubits in this system is of XY type. For λe-b <λe-b ≈ 0.72 (i.e., φdc < 0.972π) the ground state of the system corresponds to a bare excitation, while for λe-b ≥ Through the flux φn,res, the interaction strength acquires c † λe-b it corresponds to a heavily-dressed (polaronic) one. a dependence on the boson displacements un bn + bn whose form is equivalent to that of the XY -Peierls∝ model [40]. The spinless-fermion – boson coupling that The ground-state energy of the system, expressed in results from this interaction via the JW transformation units of E 10−3 δϕ2 E , was evaluated through is nonlocal in nature, in contrast to other examples of u 0 J Lanczos-type≡ exact diagonalization [29, 43] and illus- such couplings in various solid-state systems [28, 41]. trated (without the constant-energy contribution) in Effective Hamiltonian and its ground states.– To show Fig. 2. For λ λc the system has a polaron-like that the effective system Hamiltonian consists of contri- e-b e-b ground state (strongly≥ boson-dressed excitation), with butions akin to H and H [cf. Eqs. (1) and (2)], one 0 e-b its energy showing a rather weak dependence on λ . switches to the spinless-fermion representation using the e-b On the other hand, for λ < λc the ground state cor- Jordan-Wigner (JW) transformation. The latter reads e-b e-b † responds to Ψk=0 , i.e., a bare excitation with k = 0. σz =2c† c 1, σ+ =2c† eiπ Pl

In the following, a microwave-pumping based protocol tween the initial and target states (cf. Fig. 2). for the preparation of an N-qubit W state is proposed The proposed protocol is deterministic in nature and assuming that the system is initially in the vacuum state generates W -type entanglement of all the qubits in par- 0 0 e 0 ph. The external driving required for this allel. Moreover, in contrast to the typical situation in | i ≡ | i ⊗ | i purpose is assumed to be represented by the operator quantum-state control, where state-preparation times of- ten scale unfavorably with the system size, here τprep does N ~β(t) not depend on the system size at all. Finally, because Ω (t)= σ+e−iqdn + σ−eiqdn , (10) qd √ X n n  they represent ground states of the system, multipartite N n=1 W states prepared by this protocol can be expected to where β(t) describes its time dependence and the factors be extremely robust. e±iqdn account for the possibility that flip operations on Besides allowing W -state preparation, the proposed different qubits are applied with a phase difference. It system features an XY -type qubit-qubit interaction, is important to stress that the general form of driving in which opens the possibility for a universal quantum com- Eq. (10) allows one to prepare – through different choices putation [50, 51]. Because the strength of this interac- of qd and β(t) – various states of the proposed system, in- tion depends dynamically on the boson degrees of free- cluding its strongly boson-dressed ground states realized dom (photons), this system bears a formal similarity to when φdc is above the critical value. certain trapped-ion systems in which the role of bosons

The transition matrix element of the operator Ωqd (t) is played by collective motional modes (phonons) [52]. between the initial state 0 and the target state Ψk=0 Compared to its trapped-ion counterparts, this system | i~ | i≡ WN 0 ph evaluates to β(t) δqd,k=0, which indicates has an added advantage that it merely involves disper- |thati the ⊗ | preparationi of this particular target state re- sionless bosons of one single frequency, which circumvents quires only a global driving field [i.e., qd = 0 in Eq. (10)]. the spectral crowding problem resulting from the quasi- Thus, in contrast to some other schemes for W -state continuous character of phonon spectra in large trapped- preparation [22], the present one does not require a lo- ion chains [53]. cal qubit control [45, 46]. By assuming that β(t) = Robustness to losses and feasibility.– It is pertinent to 2βp cos(ωdt), where ~ωd is the energy difference between briefly address the robustness of the system at hand to the two relevant states, in the RWA these states are Rabi- possible deleterious effects of losses. To this end, it is coupled with the effective Rabi frequency βp [47, 48]. worthwhile to first note that qubit-state flips and dis- Thus, starting from the state 0 , the desired state N- placements of the resonator modes are the two leading qubit W state will be prepared| withini a time interval of sources of decoherence in this system. In addition to the duration τprep = π~/(2βp), which does not depend on N. very long T2 times of transmon (gatemon) qubits, the Taking the pumping amplitude to be βp/(2π~) = damping time of microwave photons in coplanar waveg- 10 MHz, one finds τprep 25ns, which is three orders uide resonators can reach the same order of magnitude ≈ of magnitude shorter than typical coherence times of SC as T2, with the corresponding quality factor being larger 7 qubits (e.g., for transmons T2 20 100 µs [10]). Thus than 10 [54]. Besides, the relevant excitation- and pho- ∼ − the proposed protocol should not be affected by a loss ton energy scales in this system (δω, gδω, t0/~), ex- of coherence in the system. At the same time, the ob- pressed in frequency units, are all of the order of several tained τ is sufficiently long that a leakage outside of 2π 100MHz. Thus, they far exceed the decoherence prep × the computational subspace of a single qubit can be ne- rates whose state-of-the-art values in this type of systems glected. Namely, due to the multilevel character of SC are γ 0.01 MHz [10]. Finally, in this system thermal ∼ qubits, a finite anharmonicity α E12 E01 (where Eij excitations – which at temperatures typical for such SC- is the energy difference between≡ qubit− states j and i) is qubit setups (T 100 mK) have characteristic energies ∼ required. In order to avoid such a leakage, the minimal of a few GHz – can be safely neglected. Therefore, the pulse duration of tp ~/ α is necessary. For transmons loss mechanisms do not pose obstacles to realizing the (α 200MHz), even∼ a| few-nanoseconds-long| pulse is proposed system. frequency∼ − selective enough that such a leakage is negligi- Conclusions.– The present paper proposes a scheme ble [13]. The obtained τprep 25 ns suffices even in the for a fast, deterministic creation of a large-scale W -type case of gatemons, whose typical∼ anharmonicity is by a entanglement [55] in a system of inductively-coupled su- factor of two smaller than that of transmons [49]. perconducting qubits and microwave resonators. The Importantly, the large energy separation ~δω between mechanism behind this scalable entanglement resource the target state and the lowest-lying excited state of the – which allows one to engineer W states with the prepa- system ensures that the proposed W -state preparation ration times independent of the system size – stems will not be hampered by an inadvertent population of un- from the unconventional ground-state properties of a one- desired states. For instance, for δω/2π = 200 (300) MHz dimensional model describing a nonlocal coupling of a this energy separation is equal to 2Eu (3Eu), which rep- spinless fermion to zero-dimensional bosons. The feasi- resents a significant fraction of the energy difference be- bility and robustness of the underlying state-preparation 5 protocol – which only requires a global driving field – is and R. J. Schoelkopf, Charge-insensitive qubit design de- demonstrated with realistic system parameters. rived from the Cooper pair box, Phys. Rev. A 76, 042319 This study can be viewed as being complementary to (2007). that of Ref. 22, where the preparation of W states of pho- [13] For an introduction, see S. M. Girvin, in Lecture Notes on Strong Light-Matter Coupling: from Atoms to Solid-State tons – rather than qubits – was proposed. The common Systems (World Scientific, Singapore, 2013), pp. 155-206. denominator of these two proposals is that they both rely [14] See, e.g., B. Vlastakis, G. Kirchmair, Z. Leghtas, S. E. on superconducting systems and an in-situ tunability of a Nigg, L. Frunzio, S. M. Girvin, M. Mirrahimi, M. H. hopping amplitude, albeit being based on completely dif- Devoret, and R. J. Schoelkopf, Deterministically Encod- ferent physical mechanisms. These schemes are far more ing Using 100-Photon Schr¨odinger 342 scalable than the conventional ones in which resonator- Cat States, Science , 607 (2013). mediated qubit-qubit interactions are utilized to control- [15] D. M. Greenberger, M. A. Horne, and A. Zeilinger, Bell’s Theorem, Quantum Theory, and Conceptions of the Uni- lably entangle multiple qubits; such an approach was re- verse (Kluwer Academic, Dordrecht, 1989), pp. 73-76. cently used to prepare a GHZ state of 10 superconduct- [16] W. D¨ur, G. Vidal, and J. I. Cirac, Three qubits can be ing qubits [7] – the largest entanglement demonstrated so entangled in two inequivalent ways, Phys. Rev. A 62, far in solid-state architectures. Thus, the need to demon- 062314 (2000). strate the envisioned W -state preparation is compelling. [17] M. A. Nielsen, Conditions for a Class of Entanglement 83 Acknowledgments.– The author dedicates this paper Transformations, Phys. Rev. Lett. , 436 (1999). to the memory of P. W. Anderson (1923 - 2020). This [18] M. Koashi, V. Buˇzek, and N. Imoto, Entangled webs: Tight bound for symmetric sharing of entanglement, research was supported by the Deutsche Forschungsge- Phys. Rev. A 62, 050302(R) (1999). meinschaft (DFG) – SFB 1119 – 236615297. [19] J. Joo, Y.-J. Park, S. Oh, and J. Kim, Quantum telepor- tation via a W state, New J. Phys. 5, 136 (2003). [20] H. H¨affner et. al., Scalable multiparticle entanglement of trapped ions, Nature (London) 438, 643 (2005). [21] R. McConnell, H. Zhang, J. Z. Hu, S. Cuk,´ and V. ∗ Electronic address: [email protected] Vuleti´c, Entanglement with negative Wigner function [1] A. M. Zagoskin, Quantum Engineering: Theory and De- of almost 3,000 atoms heralded by one photon, Nature sign of Quantum Coherent Structures (Cambridge Uni- (London) 519, 439 (2015). versity Press, Cambridge, UK, 2011). [22] A. A. Gangat, I. P. McCulloch, and G. J. Milburn, De- [2] Y. Makhlin, G. Sch¨on, and A. Shnirman, Quantum-state terministic Many-Resonator W Entanglement of Nearly engineering with Josephson-junction devices, Rev. Mod. Arbitrary Microwave States via Attractive Bose-Hubbard Phys. 73, 357 (2001). Simulation, Phys. Rev. X 3, 031009 (2013). [3] J. P. Dowling and G. J. Milburn, “Quantum Technology: [23] C. Li and Z. Song, Generation of Bell, W , The Second Quantum Revolution,” Phil. Trans. R. Soc. and Greenberger-Horne-Zeilinger states via exceptional Lond. A 361, 16551674 (2003). points in non-Hermitian quantum spin systems, Phys. [4] W. P. Schleich et. al., Quantum Technology: from re- Rev. A 91, 062104 (2015). search to application, Appl. Phys. B 122, 130 (2016). [24] F. Reiter, D. Reeb, and A. S. Sørensen, Scalable Dissi- [5] F. Haas, J. Volz, R. Gehr, J. Reichel, and J. Est`eve, pative Preparation of Many-Body Entanglement, Phys. Entangled states of more than 40 atoms in an optical Rev. Lett. 117, 040501 (2016). fiber cavity, Science 344, 180 (2014). [25] Y.-H. Kang, Y.-H. Chen, Z.-C. Shi, J. Song, and Y. Xia, [6] N. Friis et. al., Observation of Entangled States of a Fully Fast preparation of W states with superconducting quan- Controlled 20-Qubit System, Phys. Rev. X 8, 021012 tum interference devices by using dressed states, Phys. (2018). Rev. A 94, 052311 (2016). [7] C. Song et. al., 10-Qubit Entanglement and Parallel [26] J. Chen, H. Zhou, C. Duan, and X. Peng, Preparing Logic Operations with a Superconducting Circuit, Phys. Greenberger-Horne-Zeilinger and W states on a long- Rev. Lett. 119, 180511 (2017). range Ising spin model by global controls, Phys. Rev. [8] Z.-L. Xiang, S. Ashhab, J. Q. You, and F. Nori, Hybrid A 95, 032340 (2017). quantum circuits: Superconducting circuits interacting [27] B. Fang, M. Menotti, M. Liscidini, J. E. Sipe, and V. with other quantum systems, Rev. Mod. Phys. 85, 623 O. Lorenz, Three-Photon Discrete-Energy-Entangled W (2013). State in an Optical Fiber, Phys. Rev. Lett. 123, 070508 [9] For an introduction, see U. Vool and M. Devoret, Intro- (2019). duction to quantum electromagnetic circuits, Int. J. Circ. [28] See, e.g., F. M. Souza, P. A. Oliveira, and L. Sanz, Quan- Theor. Appl. 45, 897 (2017). tum entanglement driven by -vibrational mode [10] See, e.g., G. Wendin, Quantum information process- coupling, Phys. Rev. A 100, 042309 (2019), and refer- ing with superconducting circuits: a review, Rep. Prog. ences therein. Phys. 80, 106001 (2017). [29] V. M. Stojanovi´cand I. Salom, Quantum dynamics of [11] For a recent review, see P. Krantz, M. Kjaergaard, F. the small-polaron formation in a superconducting analog Yan, T. P. Orlando, S. Gustavsson, and W. D. Oliver, A simulator, Phys. Rev. B 99, 134308 (2019). Quantum Engineer’s Guide to Superconducting Qubits, [30] C. Slezak, A. Macridin, G. A. Sawatzky, M. Jarrell, and Appl. Phys. Rev. 6, 021318 (2019). T. A. Maier, Spectral properties of Holstein and breath- [12] J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, D. I. ing polarons, Phys. Rev. B 73, 205122 (2006). Schuster, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, 6

[31] V. M. Stojanovi´c, P. A. Bobbert, and M. A. J. Michels, phonon coupling: A platform for observing sharp small- Nonlocal electron-phonon coupling: Consequences for polaron transitions, Phys. Rev. B 89, 144508 (2014). the nature of polaron states, Phys. Rev. B 69, 144302 [45] R. Heule, C. Bruder, D. Burgarth, and V. M. Stojanovi´c, (2004). Controlling qubit arrays with anisotropic XXZ Heisen- [32] V. M. Stojanovi´c and M. Vanevi´c, Quantum- berg interaction by acting on a single qubit, Eur. Phys. entanglement aspects of polaron systems, Phys. Rev. B J. D 63, 41 (2011). 78, 214301 (2008). [46] V. M. Stojanovi´c, Feasibility of single-shot realizations of [33] K. Hannewald, V. M. Stojanovi´c, and P. A. Bobbert, conditional three-qubit gates in exchange-coupled qubit A note on temperature-dependent band narrowing in arrays with local control, Phys. Rev. A 99, 012345 oligo-acene crystals, J. Phys: Condens. Matter 16, 2023 (2019). (2004). [47] See, e.g., M. Russ, D. M. Zajac, A. J. Sigillito, F. Borjans, [34] For a general discussion of such transitions, see V. J. M. Taylor, J. R. Petta, and G. Burkard, High-fidelity M. Stojanovi´c, Entanglement-spectrum characterization quantum gates in Si/SiGe double quantum dots, Phys. of ground-state nonanalyticities in coupled excitation- Rev. B 97, 085421 (2018). phonon models, Phys. Rev. B 101, 134301 (2020). [48] See, e.g., D. M. Zajac, A. J. Sigillito, M. Russ, F. Borjans, [35] For other types of tunable couplers, see M. R. Geller, E. J. M. Taylor, G. Burkard, and J. R. Petta, Resonantly Donate, Y. Chen, C. Neill, P. Roushan, and J. M. Mar- driven CNOT gate for electron spins, Science 359, 439 tinis, Tunable coupler for superconducting Xmon qubits: (2018). Perturbative nonlinear model, Phys. Rev. A 92, 012320 [49] A. Kringhøj, L. Casparis, M. Hell, T. W. Larsen, F. (2015). Kuemmeth, M. Leijnse, K. Flensberg, P. Krogstrup, J. [36] T. W. Larsen, F. Kuemmeth, T. S. Jespersen, Nyg˚ard, K. D. Petersson, and C. M. Marcus, Anhar- P. Krogstrup, J. Nyg˚ard, and C. M. Marcus, monicity of a superconducting qubit with a few-mode Semiconductor-Nanowire-Based Superconducting Qubit, Josephson junction, Phys. Rev. B 97, 060508(R) (2018). Phys. Rev. Lett. 115, 127001 (2015). [50] N. Schuch and J. Siewert, Natural two-qubit gate for [37] T. P. Orlando and K. A. Delin, Introduction to Applied quantum computation using the XY interaction, Phys. (Addison-Wesley, Reading, MA, 1991). Rev. A 67, 032301 (2003). [38] F. Mei, V. M. Stojanovi´c, I. Siddiqi, and L. Tian, Ana- [51] D. M. Abrams, N. Didier, B. R. Johnson, M. P. da log superconducting for Holstein po- Silva, and C. A. Ryan, Implementation of the XY larons, Phys. Rev. B 88, 224502 (2013). interaction family with calibration of a single pulse, [39] M. Abramowitz and I. A. Stegun, Handbook of Mathe- arXiv:1912.04424. matical Functions (Dover Publications, New York, 1972). [52] M. L. Wall, A. Safavi-Naini, and A. M. Rey, Boson- [40] See, e.g., L. G. Caron and S. Moukouri, mediated quantum spin simulators in transverse fields: Renormalization Group Applied to the Ground State of XY model and spin-boson entanglement, Phys. Rev. A the XY Spin-Peierls System, Phys. Rev. Lett. 76, 4050 95, 013602 (2017). (1996), and references therein. [53] K. A. Landsman, Y. Wu, P. H. Leung, D. Zhu, N. M. [41] See, e.g., A. P´alyi, P. R. Struck, M. Rudner, K. Flens- Linke, K. R. Brown, L. Duan, and C. Monroe, Two-qubit berg, and G. Burkard, Spin-Orbit-Induced Strong Cou- entangling gates within arbitrarily long chains of trapped pling of a Single Spin to a Nanomechanical Resonator, ions, Phys. Rev. A 100, 022332 (2019). Phys. Rev. Lett. 108, 206811 (2012), and references [54] Z. Wang, S. Shankar, Z. K. Minev, P. Campagne-Ibarcq, therein. A. Narla, and M. H. Devoret, Cavity Attenuators for [42] P. Coleman, Introduction to Many-Body Physics (Cam- Superconducting Qubits, Phys. Rev. Appl. 11, 014031 bridge University Press, Cambridge, UK, 2015). (2019). [43] J. K. Cullum and R. A. Willoughby, Lanczos Algo- [55] F. Fr¨owis, P. Sekatski, N. Gisin, W. D¨ur, and N. San- rithms for Large Symmetric Eigenvalue Computations gouard, Macroscopic quantum states: Measures, fragility, (Birkh¨auser, Boston, 1985). and implementations, Rev. Mod. Phys. 90, 025004 [44] See, e.g., V. M. Stojanovi´c, M. Vanevi´c, E. Demler, and (2018). L. Tian, Transmon-based simulator of nonlocal electron-