# Observation of Hopping and Blockade of Bosons in a Trapped Ion Spin Chain

## Abstract

The local phonon modes in a Coulomb crystal of trapped ions can represent a Hubbard system of coupled bosons. We selectively prepare single excitations at each site and observe free hopping of a boson between sites, mediated by the long-range Coulomb interaction between ions. We then implement phonon blockades on targeted sites by driving a Jaynes-Cummings interaction on individually addressed ions to couple their internal spin to the local phonon mode. The resulting dressed states have energy splittings that can be tuned to suppress phonon hopping into the site. This new experimental approach opens up the possibility of realizing large-scale Hubbard systems from the bottom up with tunable interactions at the single-site level.

Trapped atomic ions are an excellent medium for quantum computation and quantum simulation, acting as a many-body system of spins with programmable and reconfigurable Ising couplings Islam2013 (); Monz2016 (); Debnath2016 (). In this system, the long-range spin-spin interaction is mediated by the collective motion of an ion chain and emerges over time scales longer than the propagation time of mechanical waves or phonons through the crystal Zhu2006a (); Zhu2006 (). On the other hand, at shorter timescales, such a chain represents a bosonic system of phonon modes that describe the local motion of individual ions. Here each local mode is defined by the harmonic confinement of a particular ion with all other ions pinned. In this picture, phonons hop between the local modes due to the long-range Coulomb interaction between ions Porras2004 (); Ivanov2009 (); Mering2009 (). This intrinsic hopping in trapped ion crystals makes it a viable candidate for simulating many-body systems of bosons Porras2004 (); Ivanov2009 (), boson interference Toyoda2015 () and applications such as boson sampling Shen2014 ().

Such a system of local oscillators can be approximated to the lowest order of the transverse ion displacement by the phonon Hamiltonian (),

(1) |

Here the local mode frequency of each ion is expressed as a sum of the common mode transverse trap frequency and a position-dependent frequency shift experienced by the th ion Porras2004 (); Ivanov2009 (). The local mode bosonic creation and annihilation operators are and , respectively. The long-range hopping term is determined by the distance between ions and , where and are the charge and mass of a single ion.

By applying external controls to the system, on-site interactions between phonons lead to the simulation of Hubbard models of bosons. For instance, applied position-dependent Stark shifts can result in effective phonon-phonon interactions Porras2004 (). Combined with phonon hopping between sites, such a system follows the Bose-Hubbard model. In the approach considered here, the internal spin is coupled to the external phonon mode by driving the spin resonance on a motion-induced sideband transition Leibfried2003 (). This gives rise to nonlinear on-site interactions between spin-phonon excitations (polaritons). Such a system simulates the Jaynes-Cummings-Hubbard model, which describes an array of coupled cavities Ivanov2009 (); Mering2009 (); Greentree2006 (); Hartmann2006 (); Toyoda2013 ().

In order to study the dynamics of such bosonic systems, the local phonon modes must be manipulated and detected faster than the hopping rate. Addressing these modes requires fields that target individual ions in space and each local mode in frequency. Previous experiments have observed hopping only between two sites by either using the collective motion of ions in separate but nearby trapping zones Brown2011 (); Harlander2011 () or by varying the spacing between ions in the same trap to spatially resolve each site Haze2012 (). In contrast, the direct addressing of each local mode in a single ion crystal circumvents such overheads and provides a complete toolbox for implementing larger bosonic systems by simply increasing the number of ions.

In this Letter, we report the observation of free phonon hopping in an ion chain, and study its suppression by applying targeted phonon blockades on individual sites. We access all local motional modes along the transverse direction of a static linear chain of ions. Phonons are prepared and measured by driving sideband transitions on each mode faster than the rate of hopping in the chain, with an overall fidelity of . This is achieved by setting the transverse confinement to be much larger than the axial confinement of the chain, which results in suitably low hopping rates due to both a relatively large inter-ion distance (m) and high mass of the ion.

A phonon blockade is implemented by resonantly driving red-sideband transitions on each site, which couple the internal spin of individually addressed ions to their local phonon mode via a Jaynes-Cummings interaction. In the rotating frame of the free spin and the transverse motional common-mode Hamiltonian (), this interaction is represented by the blockade Hamiltonian as

(2) |

Here, the spin- ‘ground’ and ‘excited’ states of the th ion are represented by and , respectively, with energy splitting , and spin raising and lowering operators and . A local motional red sideband is driven at a Rabi frequency and detuned from resonance by .

Phonon blockades are applied on individual sites that have ions prepared in the ground state of spin and motion , where the second index denotes the local mode phonon number. Upon applying the Jaynes-Cummings interaction at resonance (), a maximal energy splitting of occurs between and the next excited polaritonic states . This energy cost suppresses phonons from entering the targeted sites and thereby creates a blockade (see Fig. 1a). This scheme is analogous to implementing photon blockades using single-atom cavity QED systems Hamsen2017 ().

The experiment consists of a linear chain of three ions, each with an internal spin defined by a pair of hyperfine ‘clock’ states as and of the electronic ground level with a hyperfine energy splitting of GHz Olmschenk2007 (). Here, and denote the quantum numbers associated with the total atomic angular momentum and its projection along the quantization axis defined by an applied magnetic field of G. The external motion of the trapped ions is defined by a linear rf-Paul trap with transverse (X,Y) and axial (Z) harmonic confinement at frequencies MHz such that the ion chain is aligned along Z with a distance of m between adjacent ions. During an experiment, we excite local phonons in the transverse modes along X, which can then hop between the ion sites. The inherent hopping rates are approximately kHz and , respectively. The combined effect of the transverse (X) harmonic confinement and repulsion between ions (determined by ) define the position-dependent local mode frequency shifts . Fig. 1a represents the local modes with frequencies in a frame rotating at the common mode frequency .

Coherent control of the spin and motion of each ion is implemented with stimulated Raman transitions using a nm mode-locked laser Hayes2010 (), where pairs of Raman beams couple the spin of an ion to its transverse motion Debnath2016 (). A global beam illuminates the entire chain, and a counterpropagating array of individual addressing beams is focused to a waist of m at each ion. The beat note between the Raman beams can then be tuned to to implement a “carrier” transition for coherent spin flips, or tuned to to drive a blue- or red-sideband transition involving local phonon modes. The individual addressing beams are modulated independently using a multi-channel acousto-optic modulator ^{1}^{2}

A typical experimental sequence, as shown in Fig. 1b, starts with the preparation of each ion in state by Doppler cooling and subsequent Raman sideband cooling of each of the transverse modes. A single phonon excitation is introduced at a single site by resonantly driving a blue-sideband and carrier pulse to prepare the state . In order to minimize the effect of hopping during this process, the sideband and carrier pulses are kept short (s and s, respectively). Phonon blockades are applied to particular ions, initially prepared in the state, by resonantly driving the red-sidebands of their respective local modes. Finally, the single phonon occupancy denoted by states and is measured at each site using a red-sideband pulse on each ion, which coherently projects it to spin states and , respectively. The spin-dependent fluorescence can then be detected using a multi-channel photomultiplier tube, thereby measuring a binary phonon occupancy of 0 or 1 for each site Olmschenk2007 (); Debnath2016 ().

Parameter | Fitted value | Measured value |
---|---|---|

— | ||

— | ||

Figure 2 shows the hopping dynamics. During free hopping, a single excitation is observed to hop predominantly to the neighboring site. The extent of hopping is indicated by the amplitude of the oscillations in phonon occupancy. This is determined by the strength of hopping relative to the energy splitting between local modes . We observe different hopping rates between ions 1 and 2 compared to that between 2 and 3, which indicates an asymmetry in the local mode energy differences, . This is likely due to a stable non-linearity in the transverse confinement of the ion trap. We also note that the sign of the local mode energy difference is critical in governing next-nearest neighbor hopping in systems with three or more modes. This is due to a Raman-type hopping process where appropriate energy splittings between the local modes can facilitate hopping between ion 1 and 3 via ion 2 (see supplementary material).

Phonon hopping is also observed in the presence of a blockade applied on neighboring sites (Fig.2d-g). Here, we resonantly drive on the red sideband, creating a ladder of Jaynes-Cummings eigenstates , where is a spin-phonon dressed state with polaritonic excitation number (Fig.1a). Since the blockade ion is initially in eigenstate , hopping into this site is suppressed when the energy splitting of the first excited states is much larger than the hopping rate. This implies that when , hopping is suppressed from the th to th site, where the tunable blockade strength is set by the red-sideband Rabi frequency at the blockaded site . For ions and , a higher suppression is observed compared to ions and , where the large phonon mode splitting results in some residual hopping despite the applied blockade (see supplementary material).

In Fig.2, we further observe nonzero phonon occupancies at sites prepared in the motional ground state owing to imperfect initial sideband cooling. Based on the non-oscillatory near-zero phonon occupancies of the hopping data (Fig.2 a, c, and g for ions , , and , respectively), we estimate average local mode phonon numbers of . This also leads to an imperfect preparation and measurement of state , which additionally suffers from residual phonon hopping over the finite duration of sideband -pulses.

We fit the hopping data to a theoretical model using the Hamiltonian presented in Eq.1 and Eq.2, with , and as free parameters. The steady-state phonon occupancy, in the absence of hopping (see Fig.3), is used to characterize the systematic errors in the preparation and measurement of the states and . This is subsequently incorporated into all theoretical curves in Fig.2 Shen2012 () (see supplementary material). The effect of imperfect sideband cooling is directly included by starting from a thermal phonon distribution with mean phonon numbers . Both the spin and the motional degrees of freedom are considered in the theoretical simulation of the full hopping dynamics. A single set of parameters (see Table 1) is used to fit all data in Fig.2, which shows a very good agreement between experimental data and the theoretical description of the system. The hopping and blockade strengths given by the parameter values of and are consistent with those measured directly from the inter-ion distance and red-sideband spectroscopy, respectively (see supplementary Fig.S1b).

A straightforward way to improve fidelities of local phonon operations is setting the sideband Rabi frequency much larger than the hopping strength. In this experiment, since the Raman transition drives both transverse modes, the maximum Rabi frequency is limited by the spacing between the modes kHz. In future experiments, this can be resolved by rotating the two principle axes of transverse motion (X and Y) with respect to the Raman beams such that only one mode is excited. For multi-phonon hopping experiments, it is also necessary to implement fast measurement of phonon distributions (instead of single phonon occupancy) using cascaded sideband pulses, as has been implemented on single ions An2014 (); Um2016 ().

Tunable local phonon blockades can be a useful tool in studying energy transport in ion chains Ramm2014 (). It can also, in principle, be used in mechanically isolating pairs of ions in a crystal for the implementation of fast entangling gates mediated by their local phonon modes Zhu2006 (); Wong-Campos2017 ().

This work was supported by the ARO with funds from the IARPA LogiQ program, the AFOSR MURI program on Quantum Measurement and Verification, the ARO MURI program on Modular Quantum Circuits, and the NSF Physics Frontier Center at JQI.

## Appendix A Supplementary Materials

### a.1 1. Spectroscopy of Local Phonon Modes.

Figure S1 shows the spectroscopy of the motional sidebands of local modes of each ion in a chain of three. The blue sideband spectrum (Fig.S1a) shows the two transverse motional modes of each ion separated by kHz. We choose to excite the mode, which has the higher mode frequency. The red sideband spectrum is obtained by preparing each ion in the state with a single phonon excited in the mode followed by a red sideband pulse and fluorescence detection (Fig.S1b). The red sideband Rabi frequencies of the local modes are extracted from the spectroscopy and are shown in Table 1.

The transverse common mode trap frequencies (, ) are set by the amplitude of the ion trap rf signal. This can drift over time causing the local mode frequencies to vary. Therefore, in order to mitigate errors in the spectral addressing of each local mode, the transverse confinement is actively stabilized Johnson2016 ().

### a.2 2. State Preparation and Measurement (SPAM) Error Correction.

In the experiment, the state preparation and measurement are not perfect, due to single ion fluorescence detection error, imperfect cooling (), and mapping of phonon occupancy. For detection, the phonon number states and are mapped onto the spin states and respectively.

We consider two separate single-ion SPAM errors, for and for . In order to determine these, we prepare the ion in state and and measure the fraction of bright events. Our readout errors are and . The experimental readout and are then related to theoretical predictions by Shen2012 ()

(S1) |

We now include the effect of imperfect cooling. After sideband cooling, the phonon distribution is assumed to be thermal with an average phonon number of . Then, we can consider the initial state to be an incoherent mixture of different phonon numbers, with probabilities . Hence, the effect of imperfect cooling can be directly incorporated in the theoretical model.

Finally, we extend the error model to include another layer of SPAM error in the detection of phonon occupancy due to imperfect mapping to spin states with errors and . Assuming the ground state preparation is limited only by fluorescence detection error and the effect of the thermal phonon distribution (), i.e., , we have

(S2) |

where the theoretical values and depend on the mean phonon number . From Fig. 3 of the main text, the experimentally-measured probabilities are and for the state and , respectively (averaged over time and ions). Fitting the theoretical values to those, we have two variables, and , and two equations, from which we can extract the values as and . The SPAM error matrix is

(S3) |

which is then used to scale all theoretical curves.

### a.3 3. Long-range Hopping and Raman-like Hopping.

The phonon hopping amplitude decays as , where is the distance between ions and , so direct long-range hopping . Hence, we expect the effect of long-range hopping in the three ion chain to be unimportant. A natural approach is to use the Rabi flopping of the free hopping data in Fig.2a of the main text to extract the local mode energy difference and . Similarly, and are expected to be extracted in the same way from Fig.2c. However, this picture of two-site (two-level) approximation is incomplete, primarily due to a Raman-like process of hopping, i.e., phonons hopping from ion 1 to ion 2 and then from ion 2 to ion 3. We have confirmed our observation by using the two-ion approximation to extract out the experimental parameters; the fitted parameters deviate substantially from those in Table I of the main text and the former produce poor fits to other panels of data.

In addition, the relative sign between and does not matter in the two-ion picture. However, considering the dynamics of the full system of three modes, the relative sign is important. One can again understand why this is the case via an analogy to a Raman process. Phonons can hop from ion 1 to ion 3 via ion 2. If and had opposite signs, the magnitude of would be much smaller (analogous to the two-photon detuning in Raman transitions). This would lead to a Raman-like hopping from ion 1 to 3, which would induce much stronger oscillations than are observed in the experiment.

Therefore, we need to consider phonon hopping among three ions to correctly model their full dynamics. In the theoretical fitting, we use seven free parameters to fit all experimental data in Fig.2 of the main text, taking into account the SPAM error and the thermal phonon distribution after sideband cooling.

### a.4 4. Additional Explanation on the Phonon Blockade Data.

This section provides some additional explanation on a few non-intuitive behaviors of the blockade data. To understand it better, we include a simplified three-level picture. Consider Fig.2d of the main text. Three collective states are predominantly involved: , , and , where the subscript labels the ion number. In this subspace, the Hamiltonian can be approximated as

(S4) |

where kHz, kHz, and kHz (the common energy is subtracted out). Instead of suppressing hopping, the local mode frequency difference facilitates hopping due to the fact that it partially offsets the blockade strength . Figures S2a and S2b compare the scenarios in which kHz and kHz, respectively. It is evident that hopping is further suppressed for the case kHz. Since , this explains why the blockade data in Fig.2d seems to have stronger oscillation than that in Fig.2g. In addition, from Fig. S2, we also see that much of the state population is stored in , which explains the observation that, for Fig. 2d, ion 1 seems to have a stronger oscillation than ion 2; the apparent non-conservation of the total phonon numbers is due to our detection method of mapping phonons to spins via a single red sideband pulse.

In addition, we can understand how the red sideband Rabi frequencies come into play for the blockade data. For the free hopping data, red sideband Rabi frequency does not play a role other than some small hopping during the short pulse period. So the free hopping data (Fig.2a-c) is not sensitive to . However, for the blockade data (Fig.2d-g), due to a much longer hopping period and the offset between local mode frequency difference and the blockade strength, the observations depend sensitively on . Therefore, we also fit the red sideband blockade strengths as free parameters in the theoretical model.

### Footnotes

- Model H-601 Series 32-Channel UV Acousto-Optic Modulator, PN: 66948-226460-G01, Harris Corporation
- Model WX1284C-1 1.25 GS/s Four Channel Arbitrary Waveform Generator, PN: 126182, Tabor Electronics Ltd.

### References

- R. Islam, C. Senko, W. C. Campbell, S. Korenblit, J. Smith, a. Lee, E. E. Edwards, C.-C. J. Wang, J. K. Freericks, and C. Monroe, “Emergence and frustration of magnetism with variable-range interactions in a quantum simulator.,” Science, vol. 340, no. 6132, pp. 583–587, 2013.
- T. Monz, D. Nigg, E. A. Martinez, M. F. Brandl, P. Schindler, R. Rines, S. X. Wang, I. L. Chuang, and R. Blatt, “Realization of a scalable Shor algorithm,” Science (80-. )., vol. 351, pp. 1068–1070, mar 2016.
- S. Debnath, N. M. Linke, C. Figgatt, K. A. Landsman, K. Wright, and C. Monroe, “Demonstration of a small programmable quantum computer with atomic qubits,” Nature, vol. 536, pp. 63–66, aug 2016.
- S.-L. Zhu, C. Monroe, and L.-M. Duan, “Trapped Ion Quantum Computation with Transverse Phonon Modes,” Phys. Rev. Lett., vol. 97, p. 050505, aug 2006.
- S. L. Zhu, C. Monroe, and L. M. Duan, “Arbitrary-speed quantum gates within large ion crystals through minimum control of laser beams,” Europhys. Lett., vol. 73, no. 4, pp. 485–491, 2006.
- D. Porras and J. I. Cirac, “Bose-Einstein Condensation and Strong-Correlation Behavior of Phonons in Ion Traps,” Phys. Rev. Lett., vol. 93, p. 263602, dec 2004.
- P. A. Ivanov, S. S. Ivanov, N. V. Vitanov, A. Mering, M. Fleischhauer, and K. Singer, “Simulation of a quantum phase transition of polaritons with trapped ions,” Phys. Rev. A, vol. 80, p. 60301, dec 2009.
- A. Mering, M. Fleischhauer, P. A. Ivanov, and K. Singer, “Analytic approximations to the phase diagram of the Jaynes-Cummings-Hubbard model,” Phys. Rev. A, vol. 80, p. 53821, nov 2009.
- K. Toyoda, R. Hiji, A. Noguchi, and S. Urabe, “Hong-Ou-Mandel interference of two phonons in trapped ions,” Nature, vol. 527, no. 7576, pp. 74–77, 2015.
- C. Shen, Z. Zhang, and L.-M. Duan, “Scalable Implementation of Boson Sampling with Trapped Ions,” Phys. Rev. Lett., vol. 112, p. 050504, feb 2014.
- D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, “Quantum dynamics of single trapped ions,” Rev. Mod. Phys., vol. 75, pp. 281–324, mar 2003.
- A. D. Greentree, C. Tahan, J. H. Cole, and L. C. L. Hollenberg, “Quantum phase transitions of light,” Nat. Phys., vol. 2, pp. 856–861, dec 2006.
- M. J. Hartmann, F. G. S. L. Brandão, and M. B. Plenio, “Strongly interacting polaritons in coupled arrays of cavities,” Nat. Phys., vol. 2, pp. 849–855, dec 2006.
- K. Toyoda, Y. Matsuno, A. Noguchi, S. Haze, and S. Urabe, “Experimental Realization of a Quantum Phase Transition of Polaritonic Excitations,” Phys. Rev. Lett., vol. 111, p. 160501, oct 2013.
- K. R. Brown, C. Ospelkaus, Y. Colombe, a. C. Wilson, D. Leibfried, and D. J. Wineland, “Coupled quantized mechanical oscillators,” Nature, vol. 471, pp. 196–199, mar 2011.
- M. Harlander, R. Lechner, M. Brownnutt, R. Blatt, and W. Hänsel, “Trapped-ion antennae for the transmission of quantum information.,” Nature, vol. 471, pp. 200–203, 2011.
- S. Haze, Y. Tateishi, A. Noguchi, K. Toyoda, and S. Urabe, “Observation of phonon hopping in radial vibrational modes of trapped ions,” Phys. Rev. A, vol. 85, p. 031401, mar 2012.
- C. Hamsen, K. N. Tolazzi, T. Wilk, and G. Rempe, “Two-Photon Blockade in an Atom-Driven Cavity QED System,” Phys. Rev. Lett., vol. 118, no. 13, pp. 1–6, 2017.
- S. Olmschenk, K. C. Younge, D. L. Moehring, D. N. Matsukevich, P. Maunz, and C. Monroe, “Manipulation and detection of a trapped Yb+ hyperfine qubit,” Phys. Rev. A, vol. 76, p. 052314, nov 2007.
- D. Hayes, D. N. Matsukevich, P. Maunz, D. Hucul, Q. Quraishi, S. Olmschenk, W. Campbell, J. Mizrahi, C. Senko, and C. Monroe, “Entanglement of atomic qubits using an optical frequency comb,” Phys. Rev. Lett., vol. 104, no. 14, pp. 1–4, 2010.
- Model H-601 Series 32-Channel UV Acousto-Optic Modulator, PN: 66948-226460-G01, Harris Corporation.
- Model WX1284C-1 1.25 GS/s Four Channel Arbitrary Waveform Generator, PN: 126182, Tabor Electronics Ltd.
- A. C. Lee, J. Smith, P. Richerme, B. Neyenhuis, P. W. Hess, J. Zhang, and C. Monroe, “Engineering large Stark shifts for control of individual clock state qubits,” Phys. Rev. A, vol. 94, p. 042308, oct 2016.
- C. Shen and L.-M. Duan, “Correcting detection errors in quantum state engineering through data processing,” New J. Phys., vol. 14, p. 053053, may 2012.
- S. An, J.-N. Zhang, M. Um, D. Lv, Y. Lu, J. Zhang, Z.-Q. Yin, H. T. Quan, and K. Kim, “Experimental test of the quantum Jarzynski equality with a trapped-ion system,” Nat. Phys., vol. 11, no. 2, pp. 193–199, 2014.
- M. Um, J. Zhang, D. Lv, Y. Lu, S. An, J.-N. Zhang, H. Nha, M. S. Kim, and K. Kim, “Phonon arithmetic in a trapped ion system,” Nat. Commun., vol. 7, p. 11410, apr 2016.
- M. Ramm, T. Pruttivarasin, and H. Häffner, “Energy transport in trapped ion chains,” New J. Phys., vol. 16, p. 063062, jun 2014.
- J. D. Wong-Campos, S. A. Moses, K. G. Johnson, and C. Monroe, “Demonstration of two-atom entanglement with ultrafast optical pulses,” Arxiv 1709.05179, sep 2017.
- K. G. Johnson, J. D. Wong-Campos, A. Restelli, K. A. Landsman, B. Neyenhuis, J. Mizrahi, and C. Monroe, “Active stabilization of ion trap radiofrequency potentials,” Rev. Sci. Instrum., vol. 87, no. 5, p. 053110, 2016.