UMD-PP-021-05
Engineering an Effective Three-spin Hamiltonian in Trapped-ion Systems for Applications in Quantum Simulation
1 2 1 3, 4 5 6, B´arbaraAndrade, Zohreh Davoudi, Tobias Graß, Mohammad Hafezi, Guido Pagano, and Alireza Seif ∗ 1ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Castelldefels (Barcelona) 08860, Spain. 2Maryland Center for Fundamental Physics and Department of Physics, University of Maryland, College Park, MD 20742, USA. 3Joint Quantum Institute and Department of Physics, University of Maryland, College Park, MD 20742, USA. 4Department of Electrical and Computer Engineering and Institute for Research in Electronics and Applied Physics, University of Maryland, College Park, MD 20742, USA. 5Department of Physics and Astronomy, Rice University, Houston, TX 77005, USA. 6Pritzker School of Molecular Engineering University of Chicago, Chicago, Illinois 60637, USA. (Dated: August 3, 2021) Trapped-ion quantum simulators, in analog and digital modes, are considered a primary candidate to achieve quantum advantage in quantum simulation and quantum computation. The underlying controlled ion-laser interactions induce all-to-all two-spin interactions via the collective modes of mo- tion through Cirac-Zoller or Mølmer-Sørensen schemes, leading to effective two-spin Hamiltonians, as well as two-qubit entangling gates. In this work, the Mølmer-Sørensen scheme is extended to in- duce three-spin interactions via tailored first- and second-order spin-motion couplings. The scheme enables engineering single-, two-, and three-spin interactions, and can be tuned via an enhanced protocol to simulate purely three-spin dynamics. Analytical results for the effective evolution are presented, along with detailed numerical simulations of the full dynamics to support the accuracy and feasibility of the proposed scheme for near-term applications. With a focus on quantum sim- ulation, the advantage of a direct analog implementation of three-spin dynamics is demonstrated via the example of matter-gauge interactions in the U(1) lattice gauge theory within the quantum link model. The mapping of degrees of freedom and strategies for scaling the three-spin scheme to larger systems, are detailed, along with a discussion of the expected outcome of the simulation of the quantum link model given realistic fidelities in the upcoming experiments. The applications of the three-spin scheme go beyond the lattice gauge theory example studied here and include studies of static and dynamical phase diagrams of strongly interacting condensed-matter systems modeled by two- and three-spin Hamiltonians.
I. INTRODUCTION freedom of the physical system, can be taken into account in enhancing and extending the simulator’s toolkit. With Quantum simulation is expected to be a promising av- a focus on the analog mode of operation of the quan- enue to revealing the rich dynamics of quantum many- tum device, such an objective is applied in this work to body systems, from strongly correlated electron systems trapped-ion quantum simulators [10–12]. These systems in material [1], to dense nuclear matter in the interior have demonstrated their value in recent years in success- of neutron stars [2,3], to coherent neutrino propagation ful simulation of quantum many-body systems [13–23], from core-collapse supernovae and early universe [4,5], some of which have started to challenge classical numer- ical methods. These simulations were mostly based on arXiv:2108.01022v1 [quant-ph] 2 Aug 2021 to strongly interacting gauge field theories of the Stan- dard Model of particle physics [6–9]. Quantum simu- an effective two-spin Ising Hamiltonian with single-spin lators can perform in analog, digital, or hybrid modes, field terms. Motivated by the need for an enhanced dy- with the underlying hardware architecture ranging from namics as exhibited in certain exotic spin systems as well condensed-matter-based systems to optical, atomic, and as in lattice gauge theories, we generalize this effective molecular systems. As the vibrant field of quantum sim- Hamiltonian to a Hamiltonian with simultaneous single-, ulation moves toward making large-scale noise-resilient two-, and three-spin dynamics, and examine its accuracy quantum simulators and quantum computers a reality, it through a numerical simulation with realistic experimen- is crucial to pursue a co-design process in which the speci- tal parameters. fications of the theoretical problem subject to the simula- While low-energy descriptions of interactions in na- tion, such as characteristics of interactions and degrees of ture primarily rely on two-body couplings, three- and higher-body couplings occur in many physical systems and impact the dynamics in a nontrivial manner. In high- ∗ The authors’ list is alphabetically ordered. energy physics, quantum field theories exhibit interacting 2 terms beyond two-field couplings, including gauge-matter scheme [102]. Effective interactions proportional to + + + + and gauge-field self interactions in gauge theories of the σi σj− + h.c. and σi σj σk + h.c. are generated, along z Standard Model [24–27]. In nuclear physics, three and with a single-spin Hamiltonian proportional to σi . Off- higher-nucleon interactions that are effectively generated resonant contributions constitute interactions that entan- from quantum chromodynamics (QCD) have proven to gle the spin and phonon degrees of freedom and are en- be important in providing an accurate description of nu- sured to be kept small with a careful choice of laser and clei and dense matter [28–32]. While present-day quan- trap parameters as will be discussed. Our strategy for tum simulators are not yet suited for simulating quan- tuning the relative importance of the two- and three-spin tum fields and their complex interactions in general (see dynamics is different from Ref. [100] and does not require e.g., Refs. [7–9, 23, 33–76] for some recent progress), the addressing the axial modes of motion. In particular, by possibility of inducing simultaneous interactions among introducing an almost “mirror” copy of the drives that more than two degrees of freedom of a quantum simulator induce effective single-, two-, and three-spin interactions, can open up many interesting possibilities. A more im- the single- and two-spin Hamiltonians can be significantly mediate application is in the realm of condensed matter suppressed. To validate the effective dynamics, we go physics where it is argued that spin systems with two- beyond a qualitative discussion of requisite experimen- and three-spin Hamiltonians exhibit novel phase dia- tal parameters and present a thorough numerical simu- grams and quantum phase transitions, signified by unique lation of a three-spin coupling scheme in the trapped- entanglement characteristics that are absent in systems ion system by examining the exact dynamics (within a with only two-spin couplings [77–90]. Such dynamics are highly accurate rotating-wave approximation), including interesting in the context of spin glass physics [86–88] the effects of spin-phonon contamination. High fidelities and quantum statistics [91]. Furthermore, quantum link for the three-spin dynamics, as well as the effectiveness models that exhibit local symmetry constraints, and are of the “mirror” drive in suppressing lower-body interac- argued to approximate infinite-dimensional gauge field tions, are established for realistic experimental values for theories in certain limits [92, 93], are spin systems (after laser intensities and frequencies, and for the axial trap mapping fermions to hard-core bosons) with multi-spin frequency that controls the relative strength of contribu- interactions. tions from multiple normal modes to the dynamics. Given such rich physics exhibited in systems that are To demonstrate an application of our scheme in simpli- governed by multi-spin interactions, and given the in- fying near-term quantum simulations of physical models terest in demonstrating the near-term benefits of ana- of interest, the example of the lattice Schwinger model log (and hybrid) quantum simulators, multitude of pro- within its simplest quantum-link-model representation is posals and strategies are put forward to directly engi- examined and feasible strategies are proposed for scaling neer correlated three-qubit interactions. These ideas in- the simulations given undesired contamination from pro- clude the use of Floquet engineering in superconduct- liferation of the normal modes of motions. Furthermore, ing qubits [94], optimization of periodic drives in weakly the question of the severity of local gauge-symmetry vi- driven quantum systems such as superconducting cir- olations in light of imperfect fidelity of realistic exper- cuits and molecular nanomagnets [95], perturbative gen- iments is analyzed through a crude model of interac- eration of multi-spin interactions with tunable couplings tions. The underlying scheme for Hamiltonian engineer- in a triangular configuration of an optical lattice of two ing of this work, through analog or gate-based implemen- atomic species [96], natural generation of three-body cou- tations, can be extended analogously to engineer three- pling in polar molecules driven by microwave fields [97], spin Hamiltonians with higher-dimensional spin degrees adiabatic passages with tunneling interactions among of freedom, as well as four- and higher-body spin inter- atoms in a one-dimensional optical lattice [98], coupling actions, with direct applications in lattice gauge-theory Rydberg-pair interactions and collective motional modes simulations. In the context of bosonic quantum field the- in trapped Rydberg ion systems [99], and resonantly driv- ories such as gauge theories, another interesting extension ing three-spin interactions through an adiabatic elimina- of the scheme is the engineering of effective spin-spin- tion of the off-resonant transitions involving spin-phonon phonon dynamics which will be explored in future work. couplings in trapped-ion systems [100]. The present paper is organized as follows. SectionII More specifically, the possibility of engineering a three- contains a detailed description of the scheme and its spin Hamiltonian in a trapped-ion quantum simulator numerical simulation. In particular, Sec.IIA offers a was proposed by Bermudez et al. in Ref. [100], where the semi-qualitative description of the underlying mechanism use of first and second order spin-phonon couplings en- for the generation of an effective three-spin Hamiltonian abled the generalization of the two-spin “phase gate”[101] without full technical details. A detailed derivation of z z z z z Hamiltonian σi σj to a three-spin Hamiltonian σi σj σk, the effective Hamiltonian for single- and multi-mode sce- where i, j, k are∼ the ion indices. An analogous strategy is narios and using single- and multi-drive schemes to con- applied in this work to induce effective correlated spin- trol the three-spin dynamics is provided in Sec.IIB. Sec- flipping transitions among (quasi)spins of three ions us- tionIIC presents the numerical simulations of the exact ing nearly resonant single- and double-sideband transi- unitary time evolution of the three-ion system, includ- tions, extending upon the well-known Mølmer-Sørensen ing phonon contributions to the dynamics, and explores 3 the range of laser intensities and frequencies, as well ions used in present-day quantum simulators. In a linear as trap characteristics, that result in maximum fideli- trap with a common confining potential and with long- ties. To present an application of the engineered three- range Coulomb interactions among the ions, the motion spin Hamiltonian in the context of quantum simulation, of the ions can be described in terms of a set collec- the example of the quantum-link-model representation of tive normal modes with quantized excitations, i.e., the U(1) lattice gauge theory is studied in Sec.III, and in- phonons. As a result, the displacement of each ion from cludes a discussion of an exact mapping of the degrees its equilibrium position, ∆ri, can be expressed in terms of freedom and the Hamiltonian and the scaling specifi- of phonon creation (am† ) and annihilation (am) operators cations and limitations in Sec.IIIA, an analysis of the associated with the motion along the three orthogonal expected dynamics in the quantum link model assum- principal axes of the trap, x, y and z, with commutation ing imperfect gauge-violating interactions in Sec.IIIB, relations [am, am0 ] = [am† , am† 0 ] = 0 and [am, am† 0 ] = δm,m0 . and a brief note on how the fully digital implementation Here, x and y denote the most-confined directions in the of the model may compare with its fully analog simula- trap, which will have the same normal-mode spectra for tion within the scheme of this work in Sec.IIIC. Conclu- symmetric traps commonly used. These will be denoted sions and a discussion of future directions are presented as transverse directions. The least-confined direction is in Sec.IV. denoted as z and is named the axial direction.2 In gen- eral, the indices m run from 1 to 3N, but for the coupling scheme described below, only the modes along one spa- II. ENGINEERING AN EFFECTIVE tial direction, selected by a light field, will be relevant. THREE-SPIN HAMILTONIAN WITH TRAPPED The index m will then run from 1 to N to denote the IONS relevant set of modes. One or multiple pairs of counter-propagating laser The effective Hamiltonian that couples three spins simul- beams, indexed by L, can be introduced with wave- taneously can be engineered using tailored ion-laser inter- vector difference ∆kL, frequency difference (beatnote) actions. SectionIIA presents basic theoretical aspects of ∆ωL, and phase difference ∆φL, to stimulate the two- the ion-laser dynamics in a linear Paul trap [103], along photon Raman transition among the two states of the with a qualitative description of an enhanced Mølmer- qubit with a Rabi frequency ΩL. The beams can ad- Sørensen scheme that leads to an effective three-spin dress the ions globally, or for more flexibility in devising Hamiltonian. SectionIIB presents the exact relations for interactions individually, in which case both amplitude the effective Hamiltonians, including the desired three- and frequency control of the beams are required (and an spin Hamiltonian and accompanying single- and two-spin index i must be attached to laser’s parameters). Hamiltonians, along with their detailed derivation and a The Hamiltonian describing the ion-laser interactions description of the errors. SectionIIC presents a numer- can be written as [12] ical study of the full evolution compared with the one anticipated from the effective picture in an experimen- N nL i∆ωLt+i∆ϕL+i∆kL ∆ri tally realistic three-ion system, and identifies regimes in Hint. = ΩLe− · the tunable parameters of the scheme that are best suited i=1 L=1 X X x y z for suppressing the single- and two-spin terms, giving rise (α0Ii + α1σi + α2σ + α3σi ) + h.c. (2) × i to predominantly three-spin dynamics. Here, index L runs over nL pairs of Raman beams. α0, α1, α2, and α3 are constants related to spin- A. Ion-laser dynamics and a qualitative description dependent forces on the two states of the qubit [12] of the scheme and are controlled by the intensity, geometry, and po- larization of the laser beams. These are set to α0 = α = α = 0 and α = 1 , associated with a com- Consider a system of N ions confined in a radio- 2 3 1 2 mon choice. The operator ∆kL ∆ri can be written frequency Paul trap. The Hamiltonian of the system in 3N · as m=1 ηm,j(am + am† ), where ηm,i are Lamb-Dicke pa- absence of any interactions with external lasers can be 2 1 ∆k written as rameters defined as η = | | b . Here, b P m,i 2Mionωm m,i m,i N 3N are the (normalized) normal-modeq eigenvector compo- ω0 z 1 H = σ + ω (a† a + ). (1) nents between ion i and mode m, and Mion denotes free 2 i m m m 2 i=1 m=1 the mass of the ion. For optimal control, the ion-laser X X system must operate in the Lamb-Dicke regime where Here, σ is a Pauli operator acting on the space of the ions’ (quasi)spin, namely the qubit. ω0 is the qubit fre- quency which is of the order of a few GHz for typical 2 x, y and z in san-serif font correspond to spatial physical coor- dinates of the trap, hence labeling the Cartesian components of the wave-vector of the laser beams. x, y, and z in italic font 1 Planck’s constant ~ is set to unity throughout. correspond to the Bloch-sphere axes in the qubit Hilbert space. 4
H x x Carrier (a) / i j (b)
Beatnote frequency
+ + + H i j (c) / k (d) Carrier
Beatnote frequency
FIG. 1. (a,b) Traditional Mølmer-Sørensen scheme based on a pair of bichromatic laser beatnotes off-resonantly driving first- order spin-phonon couplings with symmetric detuning (±δ), giving rise to an effective spin-spin interaction. The two-ion case is shown for simplicity. (c,d) Generalized Mølmer-Sørensen scheme to generate an effective three-spin coupling. A second-order blue sideband is driven with twice the detuning (2δ) as the first-order red (−δ) sideband. As shown in (c), this process creates two virtual phonons with a second-order process and annihilates the same number of phonons through two first-order processes. Note that only two out of several possibilities are depicted. In all subfigures, Ωr and Ωb are the Rabi frequencies of the red and blue beatnotes, respectively. ω0 is the qubit frequency, and ω [≡ ωcom] is the transverse center-of-mass frequency.
2 1/2 3 η (a + a† ) 1. To exclusively couple the Ra- arise interactions proportional to m,i h m m i man beams to one set of normal-mode excitations, ∆kL + + can be set along only a single principal axis. σi am, σi−am† , σi am† , σi−am, (4) 1 x y In the interaction picture with respect to the free where σ± = (σ iσ ). These describe spin excitation i 2 i ± i Hamiltonian Hfree, the interacting Hamiltonian becomes of ion i accompanied by either absorption or emission of a quantum of motion in the mth mode, the so-called blue and red sideband transitions, respectively. On the N nL P3N −iωmt † iωmt ΩL i η ,i(a e +a e ) other hand, when µ = ωm but µ ωm ηm,iΩL, the H0 = e m=1 m m m 6 ± | − | int. 2 × motional modes are only virtually excited, and the first- i=1 X LX=1 order contributions in the Lamb-Dicke parameter give i(∆ωL ω0)t+i∆φL + e− − σi + h.c., (3) rise to an effective Hamiltonian proportional to
σxσx. (5) where it is assumed that ∆ω ω ω . The prime i j | L − 0| 0 on H is to denote that this Hamiltonian is in the inter- As will be shown in Sec.IIB, in the limit where the de- action picture. Tuning the Raman-beams beatnote to ω0 tuning from the center-of-mass mode, δ, satisfies: δ | | ≡ (ω0 ωm) leads to carrier (sideband) transitions, which µ ωcom ∆m for a typical mode splitting ∆m, the two- ± 0 | − | 2 can be derived by expanding Eq. (3) to (η )( (η)), ηcomΩr Ωb O O spin coupling can be approximated as δ , neglect- where η is the short-hand notation for the Lamb-Dicke 2 th ∆k parameter defined above. Similarly, an n -order side ing counter-rotating terms. Here, η = | | com 2MionωcomN band transition can be achieved by expanding Eq. (3) to q (ηn) and setting ∆ω = ω nω . O L 0 ± m Consider one set of Raman beams. Setting the detun- 3 ing µ ∆ω ω of the Raman beams to ω , there Upon tunig the Raman-beam phase differences properly. ≡ L − 0 ± m 5 is the Lamb-Dicke parameter associated to the center-of- we introduce two pairs of Raman beams, as described in mass normal mode, and Ωr(Ωb) are the Rabi frequencies the previous section. The Rabi and beatnote frequencies associated with the red and blue sideband beatnotes, re- associated with pair I are denoted as Ω and µ ∆ω r r ≡ 0 − spectively. This is the familiar Mølmer-Sørensen scheme, ωI , respectively, and those associated with pair II are that is schematically described in Fig.1(a-b). In par- denoted as Ωb and µb ω0 ∆ωII . Additionally, the ticular, it is shown how phonon emission and absorption following relations ≡ − cooperate to create a coherent coupling between the |↓↓i and two-ion states. µr = ωcom δ, (8) |↑↑i − − The key feature of this work, which resembles that µb = 2µr = 2ωcom + 2δ, (9) first proposed in Ref. [100], is the generalization of the − Mølmer-Sørensen scheme to create an effective three-spin will be applied to ensure resonant three-spin transitions, Hamiltonian. This scheme relies on second-order contri- as will be seen shortly. This choice of beatnote frequen- butions in the Lamb-Dicke parameter, interactions that cies justifies the use of ‘r’ and ‘b’ subscripts for each pair, are proportional to corresponding to (single) red and (double) blue sideband + + + + transitions of the qubit, respectively. With this setup, σi aman, σi aman†, σi am† an, σi aman, and ignoring contributions from counter-rotating waves σi−am† an†, σi−am† an, σi−aman†, σi−am† an†. (6) oscillating with a frequency of ωm or higher, the Hamil- tonian of the system at (η) in the interaction picture Here, a spin excitation can be associated with the cre- becomes O ation or annihilation of two phonons at the same time. As shown in Fig.1(c-d), this scheme amounts to driv- N N RWA i i(ωm+µr )t + ing the ion chain using two beatnotes with asymmetric H = ηm,iΩre a† σ + h.c., int. (η) 2 m i detunings (for example, 2δ for the blue and δ for the O i=1 m=1 − X X red sideband) so that a single second-order spin-phonon (10) excitation process can be combined with two first-order and that at (η2) reads absorption processes, giving rise to a coherent resonant O coupling between and states. This leads to |↓↓↓i |↑↑↑i N N N an effective three-spin Hamiltonian proportional to RWA 1 i(ω +ω µb)t H = η η Ω e− m n− int. (η2) −4 m,i n,i b + + + O i=1 m=1 n=1 σi σj σk , σj−σj−σk−, (7) X X X + amanσ + h.c., (11) 4 2 i ηcomΩr Ωb with a coupling that can be approximated by δ2 , in the limit of small detuning δ from the center-of-mass where RWA in the superscripts refers to the rotating- mode. This process generates resonant single- and two- wave approximation applied. Only the transverse normal spin interactions as well, and if the goal is to achieve a modes along the x direction of the trap are addressed with pure three-spin Hamiltonian, these contributions must be both pairs of the beams, hence the sums over modes m suppressed, particularly given that their strength is larger and n run from 1 to N. than that of the three-spin interactions. For this reason, Performing a Magnus expansion [104] of the time- a second drive is added to cancel out the undesired single- evolution operator, keeping only contributions that 1 are resonant and further ignoring contributions from and two-spin contributions that are proportional to δ and will enhance the three-spin contribution that is propor- counter-rotating terms (see discussions below), give rise 1 to an effective (time-independent) spin Hamiltonian. Ex- tional to δ2 . This additional drive is composed of two sets of Raman beams with asymmetric detunings of op- plicitly, posite sign compared with the first drive i.e., 2δ for i R t H0 (τ) +H0 (τ) dτ the blue and δ for the red sideband. This cancellation∼ − U(t, 0) = e− 0 int. O(η) int. O(η2) ∼ T is only exact for the center-of-mass mode contribution. i H(σ) +H(σσ)+ H(σσσ) t e− eff eff eff , (12) The exact relations in the next section will provide fur- ≈ ther essential detail of the single- and two-drive schemes where denotes time ordering. The single-spin Hamil- in the single- and multi-mode scenarios, and the numeri- tonianT has the form cal study of Sec.IIC will investigate the accuracy of the 2 effective three-spin dynamics. (σ) 1 Ω 1 H = η2 r (n + ) eff 4 m,i ω + µ m 2 i m m r X X 1 Ω2 B. Derivation of the effective Hamiltonian b η2 η2 (n + n + 1) + − 4 ω + ω + 2µ m,i n,i m n m,n m n r To obtain the exact relations for the effective dynamics X z 0 0 † 0 †0 of the ion-laser system within the scheme just described, ηm,iηn,jηm ,iηn ,jδωm0 +ωn0 ,ωm+ωn am an aman σi . one can start from the interaction-picture Hamiltonian in m0,n0 X Eq. (3). To engineer an effective three-spin Hamiltonian, (13) 6
The two-spin Hamiltonian reads certain lattice gauge theories. The single- and two-spin dynamics arising from Eqs. (13) and (14) can in gen- 2 (σσ) 1 Ωr 1 eral be significant compared with the three-spin dynam- Heff = ηm,iηm,j + −4 ωm + µr 2 ics from Eq. (15). Nonetheless, there is a simplifying i i=j m m,n X X6 X X limit in which a scheme can be devised to cancel the un- 2 Ωb + desired single- and two-spin contributions entirely. Be- ηm,iηn,iηm,jηn,j (nm + nn + 1) σi σj−, ωm + ωn + 2µr fore introducing this scheme, let us consider the following (14) simplification. If the spacing among the modes is much larger than δ, the only mode contributing to the effective and the desired three-spin Hamiltonian takes the form dynamics is that close to the lasers beatnote frequency. This can be achieved for the transverse modes if trap’s (σσσ) 2 + + + Heff = ηm,iηn,jηm,kηn,k ΩrΩb σi σj σk axial confinement is increased. In the scheme adopted m,n i,j,kX X above, this is the center-of-mass mode, see Eqs. (8) and 3µ + ω + 2ω (9). Ignoring contributions from all other modes, the ef- r m n . (15) 24(µr + ωm)(µr + ωn)(2µr + ωm + ωn) fective single-, two-, and three-spin Hamiltonians become 2 (σ) 1 ηcom 2 1 In these relations, nm(n) a† am(n) is the occupation- H = Ω (n + ) ≡ m(n) eff −4 δ r com 2 number operator for the corresponding phonon mode, i X and the condition µb = 2µr is imposed. In an 1 − η2 Ω2 (n2 + n + 1) σz, (16) individual-addressing scheme, the Rabi frequencies can − 8 com b com com i be set independently at the location of each ion, hence the following replacements must be taken into account: 2 2 2 2 (σσ) 1 η Ω Ω in Eq. (13), Ω Ωr(b),iΩr(b),j in com 2 r(b) r(b),i r(b) Heff = Ωr + → 2 → 4 δ Eq. (14), and Ω Ωb Ωr,iΩr,jΩb,k in Eq. (15). i j=i r X X6 Corrections to the→ above picture can be attributed to 1 2 2 1 + at least three sources. First, off-resonant contributions, η Ω (ncom + ) σ σ−, (17) 2 com b 2 i j i.e., those introducing oscillatory time dependence in the exponent of U(t, 0) (instead of linear time dependence 4 that is a feature of resonant contributions), introduce (σσσ) η H = com Ω2Ω σ+σ+σ+ + h.c., (18) corrections to the effective Hamiltonians. They include eff 16 δ2 r b i j k first-order terms in the Lamb-Dicke parameter and can i,j,kX iδt have a time-dependence as slow as e . At long times and respectively. ncom is the occupation-number operator for 4 as long as ηΩr(b)/δ remains small, these contributions the center-of-mass mode, i.e., ncom acom† acom. can be ignored. Such a slow time dependence occurs at The above single-mode approximation≡ makes it clear higher orders in the Lamb-Dicke parameter as well, but that the single- and two-spin couplings are odd under δ they will average to zero compared with the resonant δ while the three-spin coupling is even. Therefore, the→ contributions at long times, i.e., when t 1/δ. Second, net− single- and two-spin contributions can be canceled, there are higher-order resonant and off-resonant terms by introducing two extra Raman-beam pairs (I0 and II0) that are suppressed by positive powers of Lamb-Dicke pa- each with Rabi and beatnote frequencies Ωr0 = Ωr, µr0 = rameter compared with the contributions retained. The ωcom + δ, and Ωb0 = Ωb, µb0 = 2ωcom 2δ, respectively. third source of corrections to the effective Hamiltonians −However, such detunings give rise to resonant− sideband are from counter-rotating terms. While being resonant, transitions when applied along with the first set of lasers. these generally scale with inverse powers of mode fre- To avoid such undesired contributions and still achieve quency or mode-frequency separation (and combinations the desired cancellation, one can set the detuning of the of), which are suppressed compared with those included second set of lasers from the center-of-mass mode, δ0, to in Eqs. (13)-(15), which scale as inverse powers of δ, pro- be incommensurate with respect to that of the first set, δ. vided that δ is the smallest frequency scale in the set- The cancellation of single- and two-spin couplings can be ting. Given these sources of corrections to the effective still achieved by scaling the Rabi frequencies accordingly. picture above, it is important to benchmark the scheme Explicitly, assuming δ0 = qδ for non-rational constant q in a complete numerical simulation, and obtain the range Ω0 = √qΩ , µ0 = ω + qδ, of experimental parameters with which the effective dy- r r r − com namics above can be achieved. Such an investigation is Ωb0 = √qΩb, µb0 = 2ωcom 2qδ, (19) conducted and described in detail in the next section. − While the simultaneous one-, two-, and three-spin dy- namics may be relevant in describing a generic quantum 4 1 Note that for the center of mass mode, bcom,i = √ {1, 1, ··· , 1}, many-body system, the overarching goal of our study is N to achieve a pure three-spin Hamiltonian for applications hence ηcom,i can be replaced with ηcom for each i, where ηcom is in quantum computing, and in quantum simulation of defined after Eq. (5) 7
Once again, given the approximation made, it is impor- the two pairs needs to be introduced through an appro- tant to investigate how close the simplified dynamics is to priate scale factor q = 1. With this, the Hamiltonian the full dynamics, as will be studied in the next section. including both drives6 is given by The original scheme introduced earlier will be denoted N in the following as a “single-drive” scheme while the en- com iΩr iδt iqδt + H = η e− + √qe a† σ hanced scenario is denoted as a “two-drive” scheme. 2drive 2 com com i i=1 X Ωb 2 2iδt 2iqδt 2 + ηcom e + √qe− acomσ + h.c. (21) C. Numerical simulations − 4 i For convenience, we choose Ωb = 2Ωr/ηcom, such that The effective Hamiltonian obtained from the Magnus the red- and blue-sideband processes appear with equal expansion relies on the assumption that only the reso- strengths. Additionally, the scale factor is set to q = 1.3, nant combination of second- and third-order processes which as will be evident shortly, appears to be a good contribute to the system dynamics. In the present sec- choice. Assuming a transverse trap frequency ωx = 5 tion, we validate and corroborate the analytical results 2 × 2π MHz, and a recoil energy of ωrec ∆k /(2Mion) = of Sec.IIB by performing a numerical simulation of the 26 2π kHz, the Lamb-Dicke parameter≡ | for| N = 3 ions system dynamics. × 5 evaluates to ηcom = ωrec/(Nωx) 0.0416. For dif- The numerical treatment requires discretization of ≈ ferent choices of detuning δ and Rabi frequency Ωr, the time, with time steps smaller than the time scale of the dynamics are simulatedp under the two-drive single-mode fastest processes in the Hamiltonian. Therefore, to fa- Hamiltonian in Eq. (21), with the results shown in Figs.2 cilitate the simulation, the rotating-wave approximation and3. in Eqs. (10)-(11) is assumed, which is justified by the Specifically, Fig.2 shows in black the evolution of the fact that δ µ ωcom ωm. The dynamics gener- z | | ≡ | − | spin expectation value σi as a function of time if the ated by these Hamiltonians can differ from the desired system is initially preparedh i in a state 0 corre- effective spin dynamics due to off-resonant terms which sponding to all spins pointing down and|↓↓↓i the phonon ⊗ | i vac- are neglected in the effective picture. These terms man- uum. In all cases, this evolution exhibits an oscillatory ifest themselves by changes in the phonon number dur- behavior which is a manifestation of the three-spin inter- ing the time evolution. Moreover, as explained above, actions. In addition, blue and red lines in Fig.3 show the effective dynamics becomes a pure three-spin inter- the evolution if the system is initially prepared in the action when a second drive cancels the undesired res- 0 state. In this case, the dynamics is strongly onant single- and two-spin contributions in Eqs. (16)- suppressed,|↑↓↓i ⊗ | i demonstrating that the cancellation of two- (17). However, this cancellation scheme is complete only body terms indeed works very well. with respect to the contributions from the center-of-mass To compare more quantitatively the three-spin dy- mode. Therefore, the presence of more than one mode namics in Fig.2 with the prediction from Eq. (18), yields another source of error if the purpose is the design one can consider the period of the oscillations in the of a pure three-spin Hamiltonian. full and the effective dynamics. In the effective Hamil- tonian, the three-spin coupling from a single drive is 4 2 1drive ηcomΩr Ωb J 2 , and therefore, from the combination 1. Single-mode approximation 3 16δ of the two≡ drives, one expects three-spin Rabi oscillations with a period In order to discern between these two sources of er- rors, one can first study the effect of off-resonant terms π T3 = 2drive , (22) in the single-mode approximation, where the contribu- 6J3 tions from all modes except for the center-of-mass mode 2drive 1drive 1 are discarded. Up to second order in η, the Hamiltonian where J J (1+ q− ) and the factor 6 in the 3 ≡ 3 of one Raman pair reads denominator accounts for all permutations of the three spins. For the parametersp specified in Fig.2(a-d), the N expected oscillations periods are: (a) T3 22.8 ms, (b) com iΩr iδt + ≈ H = ηcome− a† σ T 19.7 ms, (c) T 11.4 ms, (d) T 19.2 ms. These 1drive 2 com i 3 ≈ 3 ≈ 3 ≈ i=1 values agree well with the observed period, showing for X Ω all cases the presence of three-spin interactions with ap- b η2 e2iδt a2 σ+ + h.c. (20) − 4 com com i proximately the strength predicted by Eq. (18). Another A second Raman pair with opposite detuning would can- cel the undesired single- and two-spin terms from the first pair, but it would also produce a new resonant second- 5 Values that are accessible in current experimental trapped-ion order term from the combination of the two pairs. There- systems operating using Ytterbium ions in radio-frequency Paul fore, as introduced in Eq. (19), an asymmetry between traps. 8
1 1 0.8 0.8 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -0.8 -1 -1 0 5 10 15 20 25 30 35 40 0 5 10 15 20 25 30 35 40
1 1 0.8 0.8 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -0.8 -1 -1 0 5 10 15 20 25 30 35 40 0 5 10 15 20 25 30 35 40
z com FIG. 2. Evolution of hσi i under the single-mode two-drive Hamiltonian H2drive in Eq. (21), where i = 1, 2, 3 is the ion index, for different choices of detuning δ and Rabi frequency Ωr. The evolution of the spins from an initial state |↓↓↓i ⊗ |0i probes the system’s three-spin interactions, and it is shown by the black curve (for ion 1), as well as fully hidden (identical) curves for ions 2 and 3. feature of the evolution plots is that, qualitatively, the os- spin dynamics occurs for the largest detuning employed, cillations are quite different in the four scenarios shown in which is consistent with the highest fidelity of three-spin Fig.2. Clean sinusoidal oscillations, with an amplitude oscillations in Fig.2(d). ranging from -1 to 1, are seen in panel (d), which among For a direct comparison between single- and two-drive all shown cases has the largest detuning (δ = 4 2π kHz). × schemes, in Fig.4 we have simulated the dynamics of In this case, the suppression of two-body terms, as shown the Hamiltonian with only one Raman pair, Eq. (20). in Fig.3(d), also works best. In panels (a) and (b) of As expected, when the system is initialized in a state Fig.2, similarly fast three-spin dynamics as in panel (d) with all-parallel spins, the single-drive scheme in Fig.4(a) are achieved, but at smaller values of the detuning. It is leads to similar three-spin oscillations as the two-drive seen that for δ = 2 2π kHz in panel (b) the results are × scheme in Fig.2(d), yet with approximately twice the still in good accordance with sine-like oscillations, but period. However, when the initial state contains anti- the fidelity of the evolution becomes significantly poorer parallel spins, as considered in Fig.4(b), the single-drive at δ = 1 2π kHz in panel (a). Similarly, the fidelity × Hamiltonian exhibits very fast two-spin dynamics which decreases if the oscillation are sped up by increasing Ωr. is absent in the two-drive scheme in Fig.3(d). In this respect, one can compare panels (b) and (c), that This interpretation of the spin dynamics is corrobo- exhibit identical detuning, but different Rabi frequencies. rated by the number of phonons generated during the Clearly, the price for the speed-up in (c) is a significant evolution, as a quantitative measure for the deviation loss of fidelity, also resulting in an increase of spurious from the ideal evolution in which the systems remains in two-body interactions, as shown in Fig.3(c). Specifi- the phonon vacuum throughout. The average number of cally, in Fig.3 we demonstrate the two-spin dynamics by phonons over a time duration of 40 ms is given in Ta- initializing the system in a mixed-spin state 0 . blesI andII for the scenarios shown in Figs.2,3 and4. As is seen in Fig.3(d) the maximum suppression|↓↑↓i of ⊗ two- | i Smaller detunings and/or larger Rabi frequencies lead to 9
1 1 0.8 0.8 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -0.8 -1 -1 0 5 10 15 20 25 30 35 40 0 5 10 15 20 25 30 35 40
1 1 0.8 0.8 0.6 0.6 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -0.8 -1 -1 0 5 10 15 20 25 30 35 40 0 5 10 15 20 25 30 35 40
z com FIG. 3. Evolution of hσi i under the single-mode two-drive Hamiltonian H2drive in Eq. (21), where i = 1, 2, 3 is the ion index, for different choices of detuning δ and Rabi frequency Ωr. The evolution of the spins from an initial state |↓↑↓i ⊗ |0i probes the system’s two-body dynamics, and it is shown by the yellow, blue and red curves for ions 1, 2 and 3, respectively. Here, the yellow curve is fully hidden by the red one. larger deviations from the phonon vacuum. The number As a result of these considerations, one might conclude of phonons oscillates around this time-averaged value but that a larger detuning is the better choice. However, the its peak remains below one at all times. anticipation from the more realistic multi-mode scheme is that a small detuning of the center-of-mass mode, to-
1 2π {δ, Ωr}[kHz] {1, 6}{2, 10}{2, 12}{4, 16}
1 2π {δ, Ωr}[kHz] {4, 16} |↓↓↓i ⊗ |0i 0.43 0.19 0.44 0.08
|↓↓↓i ⊗ |0i 0.05 |↓↑↓i ⊗ |0i 0.14 0.08 0.14 0.05
|↓↑↓i ⊗ |0i 0.05 TABLE I. The average values of the center-of-mass phonon occupation over a time duration of 40 ms for the three-spin evolution revealed using an all-spin-parallel initial state and TABLE II. The average values of the center-of-mass phonon the two-spin evolution revealed using a mixed-spin initial occupation over a time duration of 40 ms for the three- and state. The evolution occurs under the single-mode two-drive two-body evolutions revealed by different initial states under com com Hamiltonian H2drive in Eq. (21) for the detunings and Rabi the single-mode single-drive Hamiltonian H1drive in Eq. (20) frequencies used in Figs.2 and3. The Hilbert space is re- for the parameters considered in Fig.4. The Hilbert space is stricted by truncating the phonon occupation of the center- restricted by truncating the phonon occupation of the center- of-mass mode at 6. of-mass mode at 6. 10
Dicke parameter reads
N N 1 i 0 multi iωmt iµr t iµr t + H2drive = ηm,ie Ωre + Ωr0 e am† σi 0.8 2 i=1 m=1 0.6 X X N 0.4 0 1 iµbt iµ t ηm,iηn,i Ωbe + Ω0 e b 0.2 − 4 b n=1 0 X -0.2 i(ωm+ωn)t + e− amanσi + h.c., (23) -0.4 × -0.6 -0.8 Here, Ωr0 ,Ωb0 , µr0 , and µb0 concern the second Raman- -1 0 5 10 15 20 25 30 35 40 beam pair. They are defined in Eq. (19) in terms of the parameter q. Simulations of the evolution under the multi-mode Hamiltonian were performed under the same conditions 1 as for the single-mode case, i.e., with q = 1.3, ω = 0.8 rec 0.6 26 2π kHz, and N = 3 ions. In addition to the trans- verse× trap frequency, ω = 5 2π MHz, the axial trap 0.4 x × 0.2 frequency ωz also needs to be fixed, as it determines 0 the bandwidth of the phonon spectrum, i.e., transverse -0.2 modes separations. The results for different ωz values are -0.4 shown in Fig.5 for an all-spin-parallel initial spin state, -0.6 and in Fig.6 for different mixed-spin initial states. For -0.8 the detuning and the Rabi frequencies, those parameters -1 which worked well in the single-mode approximation have 0 5 10 15 20 25 30 35 40 been chosen, i.e., δ = 2 2π kHz with Ωr = 10 2π kHz and δ = 4 2π kHz with× Ω = 16 2π kHz. × × r × From an initial state 0, 0, 0 (where the latter |↓↓↓i ⊗ | i z ket corresponds to the vacuum of three phonon modes), FIG. 4. Evolution of hσi i under a single-mode (center- of-mass) Hamiltonian with a single drive corresponding to one observes three-spin oscillations with good fidelity, com H1drive in Eq. (20), where i = 1, 2, 3 is the ion index, for de- as shown in Fig.5(a)-(c). This is not surprising, since tuning and Rabi frequency shown. (a) The evolution from an two-spin interactions are automatically suppressed by the initial state |↓↓↓i ⊗ |0i probes three-spin interactions. They choice of initial state, and off-resonant contributions from give rise to three-spin oscillations, represented by the black the other modes are negligible due to their large detun- curve for ion 1. Ions 2 and 3 behave identically. (b) The two- ing. Notably, the evolution of the three ions is almost spin dynamics is probed by considering the evolution from an exactly the same in these cases, that is, one can hardly initial state |↓↑↓i ⊗ |0i. Yellow, blue, and red curves corre- discern the black curve in all panels of Fig.5 for ion 1, spond to ions 1, 2, 3, respectively. Note that the yellow curve from the green and the purple curves for ion 2 and ion 3. is fully hidden by the red one. These plots should be com- pared with Figs.2(d) and3(d), that is the evolution under An interesting observation is that the oscillation am- two-drive scenario using the same laser parameters. plitude when evolving under the multi-mode scenario is larger than that for the single-mode approximation with similar parameters (compare e.g., Fig.2(b) and Fig. gether with a large phonon bandwidth, will be a feasi- 5(b)). This might seem counterintuitive given the ex- ble strategy to suppress the undesired effect of the other act cancellation of single- and two-spin dynamics in the modes. In particular, within the two-drive scheme these single-mode approximation. An explanation of this be- modes spoil the cancellation of the two-spin dynamics. havior could be a destructive interference of the two fac- The optimal set of parameters for an experimental im- tors which contribute to the suppression of the oscilla- plementation should, therefore, be concluded from the tions’ amplitude hence decrease in fidelity: i) the effective analysis of the multi-mode scenario, as will be presented single-spin terms in Eq. 13, and ii) noise, i.e., any terms next. which are not accounted for by the effective Hamiltonian, that involve various spin-phonon couplings. In the single- mode calculation, only the latter contribution appears as the former is fully suppressed in the two-drive scheme. 2. Multi-mode simulation In the multi-mode case, both factors contribute, and it is plausible that they can partially cancel each other out. Next, let us turn our attention to the multi-mode Let us now probe the two-body dynamics by initializ- Hamiltonian, which up to second order in the Lamb- ing the system in a mixed-spin state. Since the second 11