Quantum computation and simulation with vibrational modes of trapped ions

Wentao Chen(陳文濤) Jaren Gan Jing-Ning Zhang(張靜寧)Dzmitry Matuskevich and Kihwan Kim(金奇奐)

1Department of Physics,Tsinghua University,Beijing 100084,China

2Centre for Quantum Technologies,National University of Singapore,3 Science Drive 2,117543 Singapore,Singapore

3Beijing Academy of Quantum Information Sciences,Beijing 100193,China

4Department of Physics,National University of Singapore,2 Science Drive 3,117551 Singapore,Singapore

Keywords: quantum computation,quantum simulation,trapped ions,vibrational modes

1. Introduction

The trapped-ion system is a leading platform for the construction of universal quantum computers,offering long coherence time,high-fidelity quantum operations,and reliable state preparation and measurement.[1–5]Recently, the coherence time of a single ion-qubit has been enhanced from 600 s[6]to 5400 s,[7]a universal set of high-fidelity gates reached to 99.999% and 99.9% for single qubit and two-qubit gates,respectively,[8–10]the number of entangled qubits also has been increased to 24 qubits,[11,12]and the number of qubits in quantum simulation has been reported to up to 53 qubits.[13,14]

While previous research focused mostly on internal statebased quantum information processing, recent experimental advances in the precise control of the motional modes of trapped ions has led to motion-based quantum computation and quantum simulation as well, exploiting an infinitedimensional Hilbert space and thus potentially yielding greater computational power. Specifically,Fock states of the motional modes with～10–100 excitations or more,[15,16]also called phonons,have been prepared,manipulated,and measured with high fidelity in ion-trap experiments. With 3 motional modes and up to 10 phonons each, the dimension of available motional Hilbert space of a single trapped ion is on the order of 103,which is roughly equivalent to that of ten qubits. It is thus natural to anticipate that one would require a smaller number of ions for the motion-based trapped-ion quantum computer to outperform classical computers in some specialized tasks,compared to its internal state-based counterpart.[17,18]

The idea that the states of the harmonic oscillators can be used for quantum simulation and quantum computations was first proposed in optics.[19]The mode of the electromagnetic field can be described as a harmonic oscillator, and it is possible to implement several forms of interaction between the modes with linear optic elements, such as beamsplitters. It was recently suggested[20]that these systems were able to carry out computations for various problems,such as linear component analysis and solving linear equations, and even demonstrating quantum advantage using the boson sampling approach.[21]However, the technical challenges faced in optical systems, such as photon losses, imperfect photon detectors,non-deterministic single photon generation, and weak nonlinear interaction between light modes make it difficult to scale up photonic systems[22,23]to reveal the quantum advantages.[24]Recently,instead of the standard boson sampling with the input of photon-number states, a Gaussian boson sampling with the input of squeezed vacuum states over 50 modes into a 100-mode interferometer has been demonstrated.[25]The Gaussian boson sampling might be less problematic for photon loss, but it has not yet been fully resolved.

In contrast, trapped ions systems are well developed and widely adopted for the preparation and characterization of non-classical states of harmonic oscillators.[1,26–30]Moreover, it is possible to introduce strong nonlinear interactions,essential for universal quantum computation and quantum simulation,[19,31–33]among vibrational modes of trapped ions through the laser-induced internal state-motion coupling[1,34]as well as via anharmonicity of the Coulomb interaction between trapped ions.[35–38]In conclusion,the above mentioned features have made motional states of trapped ions a highly attractive platform for the construction of a bosonic quantum computer and quantum simulations.

In this article,we review some recent progress in quantum computation and quantum simulation with vibrational modes of trapped ions. In Section 2,we provide a general framework for the motion-based trapped-ion quantum computer. In Section 3,we introduce the trapped-ion system and the initialization and detection techniques for the motional states of trapped ions. Then we elaborate on Gaussian and non-Gaussian operations on the motional states in Sections 4 and 5, respectively. In Section 6,we describe some applications of motionbased trapped-ion quantum computation and quantum simulation. Finally, we conclude the article with some discussion and outlook.

2. Universal quantum computation and simulation with phonons

Quantum computation and quantum simulation with phonons can be regarded as a physical realization of continuous variable(CV)quantum information processing,which has been developed in various physical platforms such as optical systems[39–41]or microwave modes of electromagnetic radiation and spin ensembles.[42]Recently,there have been theoretical developments[43]for quantum computation with continuous variables that includes a proposal for universal computing independent of encoding.[44]

Similar to the qubit-based circuit model,CV quantum information processing also consists of three essential steps including initialization(or state preparation),unitary manipulation (or operation) and detection (or measurement) of the final quantum states.[19]It has been clearly shown that in the CV model, it is necessary to include any form of nonlinearity to perform universal quantum computation or non-trivial quantum simulation, which cannot be simply computed by classical computers.[31]The nonlinearity can arise from the choice of initial states, number resolving measurements, or non-Gaussian operations.

Nonlinear states refer to Fock states,or certain superpositions thereof,that cannot be generated by Gaussian operations.Gaussian operations,such as displacement and squeezing operations,are unable to provide the necessary nonlinearity.Single phonon detection and number resolving detection can also introduce the necessary nonlinearity in the protocol.Additionally, non-Gaussian operations can be achieved either by coupling motional modes to the internal spin state of the ion,or by implementing nonlinear interactions among motional modes.We will discuss these requirements in details in the following sections.

3. Encoding, initialization, and detection of phonon number states

3.1. A trapped ion system

In a trapped ion system, a qubit can be realized by using two stable energy levels separated by an optical transition (quadrupole or octupole transition), a hyperfine transition, or a Zeeman transition in an atomic ion. In this review,the atomic ion species171Yb+is used as a concrete example, but all the discussion here can be directly applied to any other ion-qubits. The two levels in the hyperfine structure of2S1/2manifold are usually used to realize a qubit, which are denoted as|↓〉≡|F=0,mF=0〉and|↑〉≡|F=1,mF=0〉with a frequency difference ofωHF=(2π)12.6428 GHz. Typically,the ion is trapped by a radio-frequency Paul trap,which provides a harmonic trapping potential with a trap frequencyωm(={ωX,ωY,ωZ})in the MHz range. The state of the system is represented in Fock state basis as|σ,nm〉, whereσis the state of the qubit andnmis the phonon number in the vibrational mode. Two laser beams from a pico-second pulsed laser detuned from the2S1/2→2P1/2transition are used to generate a stimulated Raman process to drive the carrier and motional sideband transitions of the ion[45]when the beatnote frequency between the Raman laser beams is resonant with the qubit and sideband frequencies,respectively.These transitions are described by the time evolution of the following interacting HamiltoniansHcandHm,respectively:[1]

where?candηm?mare the Rabi frequencies of carrier and blue(red)sideband transitions,σ+=|↑〉〈↓|andσ?=|↓〉〈↑|are the spin ladder operators,anda?m(am)is the creation(annihilation)operator of the motional modem.By tuning the beatnote frequencies to address the 2nd order sidebands, a twophonon creation(annihilation)operation with internal state exchange can be realized,

3.2. Ground-state cooling

Generally,manipulation of quantum states of motion begins with preparing all the relevant modes in the ground state.Initialization to the ground state is usually carried out with a series of laser cooling techniques. These cooling methods provide a practical way of directly reducing entropy of the trapped ion system, by transferring energy from the motion to the scattered photons. Laser cooling was first demonstrated using a velocity-dependent radiative force,[46,47]called Doppler cooling. The final temperature that can be achieved with Doppler cooling is limited by the natural linewidth of the transition involved.Typically,after Doppler cooling,the mean phonon number reaches roughly a few tens in typical trap frequencies of a few MHz. Sisyphus coolings[48,49]and EIT coolings[50–57]can further lower the temperature,going below the Doppler limit. For a single ion,mean phonon numbers between 1 and 0.1 have been reported after Sisyphus cooling and EIT cooling, respectively. The final step to achieving ground state is then realized by resolved-sideband cooling.[58,59]

Fig.1. Raman scheme and level diagram of 171Yb+. (a)Raman transitions are carried out via laser beams(blue arrows)directed to the ion(filled circle)after passing through acousto-optical modulators(AOMs),which are responsible for controlling the frequency difference between pairs of beams. The magnetic field is directed along a pair of counter-propagating beams,and various combinations of polarization for each Raman beam are allowed in this configuration. (b) In addition to an rf source driving two diagonally opposite rods,dc voltages are also applied to three of the linear Paul rf trap rods.This allows further control of the radial mode frequencies,ωX and ωY. LPF:low-pass filter. (c) Raman transitions (blue dashed arrows) are carried out between two levels in the 2S1/2 manifold of 171Yb+ with a frequency splitting of ωHF, denoted as |↓〉≡|F =0,mF =0〉 and |↑〉≡|F =1,mF =0〉.The motional modes X,Y,and Z with corresponding mode frequencies ωX,ωY,and ωZ can also be involved in the Raman process by controlling the frequency difference between the two Raman beams. In this manner,single and two mode quantum operations can be implemented. Adapted from Ref.[36].

3.3. Encoding scheme of phonons and initialization

There are various initial phonon states depending on the types of quantum simulation and quantum computation. For boson sampling and Gaussian boson-sampling experiments,the initial state can be a Fock state,a coherent state,a squeezed vacuum state,a displaced squeezed state,or any combination of them. For the machine learning algorithm,a classical vectorαnofncomplex numbers is mapped to the quantum state∑n αn|n〉,where|n〉denotes the Fock state. For quantum computation with bosonic modes,a Schr¨odinger cat state,[60]a binomial state,[61]and Gottesman, Kitaev and Preskill (GKP)states[62]have been used for logical qubits.

Most of these initial states can be implemented by applying either a coherent laser beam, or Raman laser beams on the ions, which creates a Jaynes–Cummings type of coupling Eq. (2) that implements displacement and squeezing operations.The details of these operations are elaborated in the next section. An excited Fock state(n/=0)and an arbitrary superposition of Fock state with phonons can be efficiently generated by the scheme demonstrated in the Be+ion system.[63]With the assistance of the ion’s internal state, an arbitrary quantum superposition of Fock states can thus be created using interactions of the Jaynes–Cummings type. By applying the displacement operation on the vacuum state, a coherent state can be obtained.[64]The Schr¨odinger cat state can be generated by an internal state dependant displacement operation,[60]typically known as the spin-dependant force.[34,65]Additionally,the squeezing operation generates squeezed vacuum-state and displaced-squeezed state if it is applied to a vacuum, or coherent state,[1,64,66]respectively. However,we note that these initial states can also be generated by specially designed bathengineering.[67]

3.4. Detection of phonon states

In an optical system, photon detectors or photon counters are used to measure photon states. Based on the measured photon results, various tomographies of photon states including Wigner function,Q-function,P-function,and Fockstate tomography are performed. However, for phonons in a trapped ion system, such phonon detectors or counters do not exist. The phonon number states, and hence its population distribution, are detected through their interaction with the internal states of the ions. Then, similar to the optical system, the measured phonon-number states can be used to reconstruct the phonon density matrix,Wigner-functions,andQ-functions. The detection methods of phonon numbers can be categorized as follows:

Phonon state detection by interferenceThe traditional phonon state detection is performed by observing the phononnumber dependent Rabi oscillations, typically by driving the blue-sideband transition. More precisely, the probability of occupying a phonon number state|n〉,Pn, can be determined from a Fourier transform of the time evolution ofP↑(t), the probability of detecting the ion in the|↑〉state[64]

Phonon number resolving detection by non-linearityProjective measurements of the motional state of ions can be realized through an effective cross-Kerr coupling to a second motional mode. Direct single-shot measurements of phonon number states and the parity of the harmonic oscillator (defined as the difference between the even and odd phonon number state occupations)can be realized by utilizing the nonlinear interaction among motional modes. This detection scheme involving parity measurements has been used to reconstruct the Wigner function of a motional state.[36]The ability to measure parity of the oscillator is discussed in more details in section 4 and Ref.[38].

Fig. 2. (a) The phonon subtraction operation that changes phonon state from|n〉to|n?1〉is composed of resonant π-pulses on the carrier transition(black arrow)followed by an adiabatic drive on the blue-sideband(grey arrow). The subtraction|n〉→|n?1〉is performed for every Fock state apart from |n=0〉, which is driven to the state |↑,n=0〉. (b) On application of the detection laser(c),the fluorescence observed is attributed to the population in the state|↓,n=0〉prior to phonon subtraction. (d)No fluorescence is detected from the other Fock states,which have had one quanta removed after the operation,and are shelved to the|↓〉state. The initial state of|↓,n〉thus arrives at |↑,n=0〉 after (n+1) repetitions of (a) and (b) procedures,and is detected via fluorescence. Adapted from Ref.[26].

Phonon number resolving detection by phonon subtraction operationsBy implementing a repeated sequence of phonon subtraction and qubit state detection, the phonon number state can be determined (and projected) by counting the number of repetitions before fluorescence is first observed.Upon detection of fluorescence in the(n+1)-th iteration,the phonon number state will have been the Fock state|n〉. The phonon subtraction|n〉→|n?1〉is performed by the successive application of a carrierπ-pulse driving|↓,n〉 →|↑,n〉,and a uniform blue-sideband transition that drives|↑,n〉 →|↓,n?1〉as shown in Fig.2. After a single instance of phonon subtraction, only the state|↓,0〉is transferred to the bright|↑,0〉state that scatters photons at the qubit-detection step.This scheme has been extended to map phonon number states to auxiliary internal states.[70]We note that if the sequence is applied only once, the population of the ground-state|0〉is measured, which allows reconstruction of theQ-function.[28]The method of projective measurement can further be extended to enable correlation measurements of multiple modes.For theNmotional mode detection,the use ofNqubits[26–28]is required.

4. Gaussian operations

A Gaussian operation transforms Gaussian states into Gaussian states. Typical examples of Gaussian operations are displacement and squeezing operations. There exist many schemes to realize these Gaussian operations for vibrational modes of trapped ions.[1,60,64,71]For example,a displacement operation can be implemented by applying an electric field for a certain duration in a desired direction. A squeezing operation can be realized by modulating the trapping potential at twice the trap frequency, also known as a parametric drive.These Gaussian operations can also be realized by applying two Raman laser beams with suitable polarization and frequency difference between them. The implementation of displacement or squeezing operations can be state-independent if the frequency differences areωmor 2ωm,similar to the electric field displacement or parametric drive.[1,72]However, these state-independent operations are challenging to incorporate together with blue- or red-sideband transitions, when welldefined phase relations between them are required. Therefore, the state-dependant operations are more appropriate for the general manipulation of phonons for quantum computation and quantum simulations.

The state-dependant Gaussian operations require one frequency in a laser denoted Raman 1,ωR1, and two frequencies and phases in the other laser denoted Raman 2,ωR2,1andωR2,2, where it is assumed that Raman 1 and Raman 2 are propagating towards the ion from two different directions(Fig. 1(a)). The Hamiltonian of an ion in a harmonic oscillator with the Raman laser beams can be expressed with the following equation:

whereαm=t(1/2)m?andφB=φ2?φ1=0. This is also known as theσx-dependent displacement operation.[34,65]Theσz-dependent displacement operation can also be realized with additionalπ/2 carrier rotation pulses(alongσyandσ?yaxis)before and after theσx-dependent displacement operation. By preparing the qubit in an eigenstate(σx),the displacement operation can be implemented.

Similarly, by settingδ1=ωm ?δSandδ2=?ωm ?δS,as shown in Fig. 3(b), the squeezing operationof a single mode[30,64,66]can be realized as

TypicallyδSis chosen to be a few times that ofηm?. Similarly,the squeezing operation can be implemented by preparing the qubit in an eigenstate(σz).

Fig. 3. Trapped-ion implementation of the quantum optical operations. The Hilbert space is composed of two phonon-modes, 1 and 2, and the internal electronic state, |↑〉 and |↓〉. The quantum operations are implemented via control of the frequencies ωR2,1, ωR2,2 and phases φ1, φ2.(a)Coherent displacement operation and(b)squeezing operation on the mode 1 as an example. (c)Rotation operation between mode 1 and mode 2. Adapted from Ref.[30].

4.1. Rotations

4.1.1. Rotations of collective modes

The rotation of collective modes can be realized by similar Raman processes.[30]By settingδ1=?ω1?δRandδ2=?ω2?δRas shown in Fig.3(c),where the magnitude ofδRis on the order of a few times ofη1,2?1,2,the rotation operation can be written as

4.1.2. Rotations of local modes

In a Coulomb crystal made up of ions,each one of them experiences a confining harmonic potential arising from the electrodes, as well as a repulsive potential due to the other ions. The motion of ions at each site can be treated separately from the others,and these are known as local modes.[73,74]For a linear ion chain,considering the motions in the transverse direction, the Hamiltonian of such a system can be represented as

whereωmdenotes the transverse trap frequency,νiis a frequency shift determined by the position of thei-th ion,andκij=e2/(2Mωmd3i j) describes the hopping rate of phonons[75,76]between ionsiandjwith a separation distancedi j.eandMare the charge and mass of a single ion, respectively. To study the interaction between the local modes,the time taken for both manipulation and detection operations should be much shorter than the time taken for propagation of mechanical waves between ions. This implies that, when compared to using collective motional modes,ions should be separated by a larger distance,which typically results in a hopping rate of kHz level.

The hopping of phonons from site to site can be viewed as an energy exchange among local modes, and this effect in turn finds applications in quantum simulation and quantum information processes. It is also possible to generate entangled states between ions in this manner. The first experimental demonstration between local modes was carried out in a linear trap,with two ions held in two separate potential wells.[73,74]With precise control of electric potentials, the harmonic frequency of two ions can be made equal, thus eliminating the energy gap between them. The Hamiltonian describing the dipole–dipole interaction between the two ions is

In general, it is challenging to control the operations on local modes in a trapped-ion system. However, similar to a photonic system, a blockade[77,79]can be realized by driving a resonant sideband operation that effectively induces a frequency shift to the motional mode. This induced shift in the levels increases the energy cost of phonons occupying the site,and hence decreases the rate of hopping.Phonon blockade can be applied to assist in state preparation of a local mode,or to artificially freeze hopping between specific neighbouring ions in a controlled manner.

5. Non-Gaussian operations

To achieve universal quantum computation and simulation, it is necessary to implement non-Gaussian operations.Such operations typically arise from nonlinear interactions in the system, either directly between motional modes, or by coupling motion with the ions’ internal spin degree of freedom. With non-Gaussian operations, one would be able to realize arbitrary Hamiltonian, and hence achieve universal quantum computation and simulation. In the following section,we review and discuss various experimental methods and techniques for implementing non-Gaussian operations that are available to the trapped ion system.

5.1. Nonlinear coupling between qubit and harmonic oscillators

5.1.1. Controlled and exponential SWAP operations

An experimental demonstration of a non-Gaussian operation resulting from nonlinear interactions between spin and motion was realized recently,[80]where a controlled beam splitter (CBS) gate was achieved. Depending on the internal spin state of the ion,the beam splitter operation acting on two motional modesaiandajis implemented according to the Hamiltonian where|n〉and|m〉correspond to Fock states of modesaiandaj, respectively. The CBS is similar to the well-studied controlled swap gate (CSWAP), which is used in a number of applications including purity measurement, overlap test,[81]quantum fingerprinting,[82]and quantum machine learning.[20]With an additional ancillary modecprepared in the vacuum state, and judicious choices of the phaseυ, the CSWAP gate can be obtained from a combination of the conditional beam splitter operation[83]

The Hamiltonian(12)is implemented by means of a modulating state-dependent optical dipole force,where the modulation frequencyωLexactly matches the frequency difference between the modesaiandaj. Interfering a pair of laser beams that have mutually orthogonal linear polarizations and a frequency differenceωLat the position of the ion generates an optical lattice with a resultant polarization that oscillates from left to right-handed circular with a frequencyωL. The state|e〉experiences a modulating ac Stark shift, which results in an optical dipole force modulated with frequencyωL, while the state|↓〉does not.[84]The choiceωL=|ωai ?ωa j|realizes Hamiltonian(12).

Fig. 4. Truth table of Fredkin gate. By limiting two modes of the ion to a single phonon each,the controlled beam splitter(CBS)is able to simulate a quantum Fredkin gate. For each of the eight basis states prepared as input,the probabilities of all outcomes after applying the CBS gate are measured via a series of projective measurements. Each data square represents an average of 10000 experiments. Without correcting for state preparation and measurement errors, an average gate success probability of 0.82±0.01 is obtained. Adapted from Ref.[80].

The Fredkin gate[85–90]is a universal 3-bit gate which swaps the values of the last two bits depending on the first.It can be implemented with a single application of the CBS gate by constraining the accessible Fock states of both motional modes to either|0〉or|1〉. Experimentally, the results are shown in Fig. 4. To determine the output of the gate after applying it to an input state, projective measurements are used to obtain the probability of each basis state. These measurements are repeated for each basis state as an input. An average of 10000 projective measurements are carried out for each data square,and without correcting for state preparation and measurement errors, an average gate success probability of 0.82±0.01 is obtained.

In addition to the Fredkin gate, the CBS gate also sees applications in other scenarios.[80]Single-shot parity measurements,[36]and hence measurements of Wigner functions, can be carried out by applying the parity operation in Eq. (14). When coupled with qubit rotation gates, the CBS gate enables SWAP test measurements which efficiently measure the overlap between two quantum states. Last but not least, a combination of qubit rotation and CSWAP gates allows one to implement the exponential SWAP gate to generate arbitrary entanglement between two modes.

5.1.2. Phonon arithmetic

Another type of nonlinear operation can be realized by a simple conventional arithmetic operation, which can be used for phonon number resolving detection. The conventional addition and subtraction of a particle can be written as

Fig.5. (a)For the uniform blue-sideband transition whose frequency is independent of motional quantum number n, ?(t) (Rabi frequency of bluesideband transition) and ?(t) (laser detuning from a motional sideband)are controlled as the red and blue curves. The phase iβ in ?(t) is the counter-diabatic term to suppress the transition during the evolution. Here?0=(2π)38.5 kHz,β =0.075,and ?0=1.6?0. (b)Experimental demonstration of the uniform blue-sideband transitions. The total time to drive the transition is 91μs for any|n〉up to n=5,which is about 7 times the π-pulse for an n=0 blue-sideband transition. Adapted from Ref.[27].

The total duration of a uniform blue-sideband transition for a complete population-transfer is about 7 times longer than theπ-pulse of the blue-sideband transition on|n=0〉,which is comparatively faster and more efficient than the standard adiabatic scheme[94]and rapid adiabatic passages.[96]Figure 5(b)shows the experimental performance of the uniform transfer of population from|↓,n〉to|↑,n+1〉in the range fromn=0 to 5.

The addition(subtraction)scheme deterministically adds(removes)one quanta independent of the initial phonon number state while maintaining quantum coherence. Various initial states including a superposition of two Fock states and coherent states have been used to test the performance of the addition and the subtraction operations.[27]For the example of Fock-state superposition shown in Figs. 6(b)–6(d), coherences represented by the off-diagonal terms of the density matrix clearly remain after the addition and subtraction operations. The density matrix is constructed by using the iterative maximum-likelihood algorithm[97]after displacing the state by 8 different angles with an amplitude ofα ～0.8, followed by determining the phonon number distributions(Eq.(4)).[71]

We note that the arithmetic operations discussed are non-Gaussian operations. However, it is not clear whether the nonlinearity introduced by them is strong enough to implement universal quantum computation and quantum simulation with phonons. These operations can make up a part of the phonon number resolving detection as discussed in Section 3,which enables phononic systems to perform algorithms of the boson-sampling type, without facing limitations of detection efficiency.

Fig. 6. (a) Sequence of an addition operation, composed of a uniform blue-sidebandtransitionfollowedbya π-pulse drivingthecarriertransition.(b)Additionsonasuperposition state(|n=0〉+|n=1〉).The reconstructed density matrices indicate a fidelity of 0.99(0.01)for the initially prepared state,and a corresponding fidelity of 0.96(0.01)and 0.92(0.01)for one and two additions, respectively. (c) Sequence of a subtraction operation,realized by the reverse of addition and fluorescence detection. (d)Subtraction on a superposition state (|n=2〉+|n=3〉). The reconstructed density matrices indicate a fidelity of 0.96 (0.01) for the initially prepared state,and a corresponding fidelity of 0.77(0.02)and 0.83(0.01)after one and two subtractions,respectively. Adapted from Ref.[27].

5.2. Nonlinear coupling between harmonic oscillators

Anharmonicity gives rise to effects which are studied in the field of nonlinear quantum optics,and are described by the well known models of degenerate and non-degenerate parametric amplifiers, as well as the cross Kerr effect. However,unlike the case in optical domain, this nonlinear coupling is strong and deterministic, and manifests itself at the level of a single quanta. In this section we show how all of these interactions can be observed in the system of trapped ions.

5.2.1. Degenerate parametric interaction

To estimate the strength of the nonlinear coupling, consider a simple case of two co-trapped ions having the same massMand chargeein a linear Paul trap,with single-ion secular frequenciesωX,ωY,ωZ. The potential energy of the system is[99,100]

whereε0is the permittivity of free space,X,Y,Zare the centerof-mass coordinates,andx,y,zare half the separation between the ions along the direction of the principle trap axes. WhenωZ ＜(ωX,ωY), the ions crystallize along the axial (z) direction at an equilibrium distancez0from the trap center. According to Eq. (16), the motion of the center-of-mass modes is harmonic, but the out-of-phase modes are coupled by the Coulomb interaction. For small axial displacementu=z?z0and keeping only terms up to the third order that contribute to the coupling between thexandzmodes, the potential energy becomes[100]

The coupling is nonlinear: if the resonance conditionδdegen=0 is met, one axial phonon is converted into a pair of radial phonons and vice versa. It is described by a Hamiltonian identical to the degenerate parametric down conversion Hamiltonian in quantum optics. Similar coupling also exist between modes of motions of three trapped ions.

Whenξd/=0, the “bare” energy eigenstates of the axial and radial modes are mixed such that eigenstates of Hamiltonian (17) have non-zero projections along both modes.This coupling manifests itself in the form of an avoided crossing[36,38]asδdegen→0,shown in Fig.7.

Fig. 7. Observation of avoided level crossing between various motional eigenstates as a function of the degenerate parametric interaction detuning from resonance δ. Shown are the probability to detect the ion in the |↑〉state as a function of detuning δdegen and Raman beams. Detunings of Raman beams are scanned around the (a), (b) first order and (c), (d) second order blue sideband transitions of the axial mode, for an initially prepared(a) |1a,0b〉, (b) |1a,1b〉, (c) |2a,0b〉, and (d) |2a,1b〉 state. Figure adapted from Ref.[38].

5.2.2. Cross Kerr effect

If the resonance conditions are not met and energy exchange between the modes is not possible,this coupling manifests itself as a shift in the frequency of a motional mode,which is proportional to the number of phonons in the other motional mode and can be described by an effective cross-Kerr interaction.[100]A small shift (about 20 Hz/phonon) of this origin was first observed in Ref. [101] using a Ramsey type experiment. By tuning the mode frequencies close to the parametric resonance,an order of magnitude increase in sensitivity to phonon number(～300 Hz/phonon)was observed in Ref.[38].

When the axial and radial modes are far detuned from each other, i.e.,δdegen?ξd, the coupling manifests itself as frequency shifts of the modes. In the lowest nonvanishing order of perturbation theory,shift of the axial modeδE=δEna+1,nb ?δEna,nb=?2(2nb+1)ξ2/δdegenis proportional to the number of phononsnbin the radial mode.[100]HereδEna,nb=∑ka,kb|〈na,nb|H|ka,kb〉|2/(Ena,nb ?Eka,kb),andEja,jb=jaωa+jbωb. This shift deviates from a linear relation due to finite value ofδdegen. However, exact diagonalization of Hamiltonian(17)confirms that the frequency shift is still a monotonous function of the radial phonon numbernb.[38]tional modes can be addressed by the Raman lasers,the bright ion is always positioned at the edge of the crystal.

Fig. 8. (a) Full frequency scan of the axial blue sideband when the system is prepared in Fock states|0a,0b〉(blue)and|0a,10b〉(green). For the state |0a,10b〉, the reduced height of the main peak and residual population detected in the states with nb ＜10 are due to imperfections of the state preparation that,according to the fit,populates the state|0a,10b〉with probability p10 =0.80(2),and leaves significant populations p9 =0.06(2),and p8 =0.06(2)in the states|0a,9b〉and|0a,8b〉. Preparation of state|0a,0b〉populates other Fock states with the probability ＜0.04. (b)Measured position of the axial motional sideband for the Fock states|0b〉to|10b〉prepared in the radial mode. The axial mode is initially prepared in the vacuum state|0a〉.Only part of the data around the largest peak is shown for each state.[38]

As shown in Fig.8,the observed frequency shifts are on the order of 300 Hz/phonon and are in good agreement with theoretical predictions calculated by diagonalization of Hamiltonian(17). The peaks are clearly separated from each other,which makes it possible to efficiently detect the phonon number distribution for the radial mode.

5.2.3. Nondegenerate parametric interaction

Symmetry and energy conservation restricts the set of modes that are coupled. The center-of-mass modes cannot be coupled to any other mode, since they are purely determined by the trap potential and do not depend on the Coulomb interaction between the ions. The energy conservation condition is fulfilled in particular by tuning the mode frequencies such that they satisfy

To observe this effect, three171Yb+ions are trapped in a standard linear rf-Paul trap[35,36]with single ion secular frequencies(ωX,ωY,ωZ)=2π×(1042,979,587)kHz.The axial trapping frequencyωZis fixed and the radial trapping frequenciesωX,ωYare fine tuned by adjusting offset voltages applied to the trap electrodes. For convenience, two of the ions are optically pumped to a metastable long-lived2F7/2state and remain dark throughout the experiment. Three Raman beams in total are used to address all nine motional modes along all three principal axes of the trap for sideband cooling,motional state preparation, and detection. To ensure that all the mo-

Fig.9. Coherent energy exchange between three modes of motion coupled by the non-degenerate parametric interaction. The purple,red,and blue dots represent the probabilities to drive the corresponding red sideband transitions of the axial zig-zag(ωa =1414 kHz), the radial tilt(ωb =878 kHz),and the radial zig-zag(ωc =536 kHz)modes, respectively, as functions of three-mode interaction time τ. The purple, red, and blue lines are the corresponding fits. The oscillations decay with a time constant of 0.11(2) s,which corresponds to more than 300 cycles of energy exchange.[102]

A signature of the cross-mode coupling is a coherent energy exchange between different modes under the resonance condition.[102]We start with one phonon in the axial mode and zero phonon in both radial modes, and then adjust the detuningδnondegen=0(resonance). After a wait durationτ,we tune the detuning back to its initial valueδnondegen=?2π×44 kHz and measure the phonon populations in all the three modes.We observe a coherent energy exchange between the axial and the two radial modes with a frequency 2.801(2)kHz(Fig.9).The amplitude of oscillations reduces to 1/e on a time scale of 0.11(2)s,which corresponds to about 330 oscillation cycles.

6. Examples of applications

6.1. Demonstration of the GKP code

Quantum error correction[103,104]is essential for the realization of scalable universal quantum computation.[105]The essence of various error-correction codes is to encode a logical qubit with a larger Hilbert space so as to protect quantum information by redundancy. Instead of using multi-qubit systems, Gottesman, Kitaev, and Preskill(GKP)[62]proposed an error-correction code that encodes a logical qubit in the Hilbert space of a single harmonic oscillator. In this error-correction code,the logical code words are coherent superpositions of an infinite series of infinitely squeezed states, i.e., eigenstates of the position or momentum operators, and their Wigner functions form an infinite lattice in the phase-space of the oscillator.

Recently, the encoding, quantum control, and readout of a GKP logical qubit have been demonstrated using a trapped40Ca+ion.[106]The logical qubit is encoded in a vibrational mode that is manipulated by laser-induced coupling with the internal states of the trapped ion. To ensure normalization and minimal effects from motional dephasing, the logical computational basis states are constructed from a finite series of squeezed states. Specifically,the approximate state to the logical GKP|0〉Lstate is the authors prepared the three-component|0〉Land the fourcomponent|1〉Lstates and measured the marginal position and momentum probability distributions.The average logical state preparation and readout fidelity was reported to be 87.3%.

The authors also demonstrated the implementation of arbitrary control and measurement of logical operators. First of all, the displacement operators, i.e., the logical Pauli operators,are implemented by applying a resonant voltage signal to one of the trapping electrodes.[71]Arbitrary control is realized by the gate teleportation technique,[108]which is widely used in error-correction codes. Here in the trapped-ion context,the gate teleportation is implemented by performing a single-qubit rotation on the internal states of the trapped ion before the MVM sequence. The measurement of logical operators, including Pauli operators and stabilizers, can be mapped to the readout of the MVM sequence by performing appropriate logical operations. The authors characterized the performance of the logical operations by quantum process tomography. The process fidelities for the logical Pauli operators are 97% and those for the partial rotations are 90%.

6.2. NOON state

A maximally entangled state of bosonic particles can be in a form of a NOON state. The NOON state of identical bosons[109,110]provides the ultimate Heisenberg limit,which,forNbosons,has the form[111]

The preparation of|0〉Lstarts with obtaining|r〉by reservoir engineering,[67]withr ?0.9. The superposition of displacement operation is then realized by modular variable measurement(MVM).[107]Two essential ingredients of the MVM sequence are the laser-induced state-dependent force,and fluorescence detection of the internal states.Repeating the MVM sequence twice and conditioned on detecting no fluorenscence,

where the relative phase?Sbetween the two terms of the state is linearly proportional toNin an interferometer, showing the Heisenberg scaling for parameter estimation. For photonic systems, experiments have demonstrated NOON states with particle numbers up toN=5.[112–119]For distinguishable particles, up to 10 photons and 14 ions have been prepared in the closely-related GHZ states.[11,120]NOON states have also been demonstrated in nuclear spins (NMR),[121]atomic spin waves,[122]and microwave photons in superconducting systems.[123]

A deterministic scheme to generate phononic NOON states with arbitrary number of bosonsN[29]has been developed based on blue-sideband transitions of Eq. (2), where two radial modes of an171Yb+ion trapped in a 3D harmonic oscillator have been used, but the generation scheme can be similarly applied to different atomic ions with more modes. In Fig.10(a), the generation scheme is illustrated for the case ofN= 3 NOON state,|3,0〉+|0,3〉, as an example. With the pulse sequence, the NOON state up toN=9 has been generated, which is mainly limited by experimental imperfections.[29]

The sensitivity of the phase estimation with the NOON states increases as the number of phononsNincreases. The phase betweenXandYmodes can be measured by the interference through the beam splitting operation of Eq.(9). In the experiment, the beam splitting operation is realized by linear combinations of input modes.[29]In the experiment,the parity of phonons in the output modes of the generated state is measured depending on the value of?, which provides the phase sensitivity proportional toNas

Figure 10(b)shows the experimental results of the parity oscillations fromN=1 toN=9 of the generated NOON states.As shown in the fitting parameterk,the enhancement of the phase sensitivity is in agreement withNwithin 2.6%deviation. AsNincreases,the contrastCPdecreases due to experimental imperfections. However, it is shown that even forN=9, the contrast is over 0.5,which indicates the existence of quantum entanglement in the state.

Fig. 10. (a) Generation sequence of the NOON state of N =3. The blue arrows indicate blue-sideband transitions, the red arrows indicate carrier transitions, and the green arrows indicate composite-pulse operations. The numbers on the arrows denote the sequence of the operations. (b) Parity oscillations of the generated NOON states from N =1 to N =9. The blue dots are experimental data, and the red lines are fitting curves with〈Π(?)〉=Acosk?+Bsink?+C,and CP= (c)The experimental results of fidelity as well as CP and PN,0+P0,N. The error bars of CP are derived from fitting error and those of PN,0+P0,N from shot-noise error.(d)The quantum Fisher information of the generated states. Figure adapted from Ref.[29].

6.3. Quantum thermodynamics

6.3.1. Test of quantum Jarzynski equality

There is increasing interest in non-equilibrium dynamics at the microscopic scale, spanning across various aspects of quantum physics,thermodynamics,and information theory as technologies of the experimental control at such scales have rapidly been developing. Most principles in non-equilibrium process are represented in the form of an inequality. In the classical regime, a remarkable equality[126]that relates the free-energy difference to the exponential average of the work done on the system was found:

whereWis the work done on the system,?Fis the free-energy difference,Tis the initial temperature of the system in thermal equilibrium,andkBis the Boltzmann constant.Independent of the protocols in which parameters of a system are varied, the Jarzynski equality Eq.(23)always holds,even when the driving is arbitrarily far from equilibrium. This is related in part to the Crooks fluctuation theorem.[127]Experimental tests of the classical Jarzynski equality have been successfully performed in various systems.[128–137]In the quantum regime, however,it is difficult to determine the work done in the system due to Heisenberg’s uncertainty principle. In an isolated quantum system, the work can be measured by two-point projective measurement over the energy eigenstates,[138,139]and the classical Jarzynski equality can be extended to the quantum regime.[140–142]In an harmonic oscillator of a trapped ion system,the projective measurements discussed in Section 3 have been applied to experimentally test the Jarzynski equality in quantum regime.

The work is performed on the system by shifting the center of the potential by a certain distance in theX-direction,which is described by the following time dependent Hamiltonian:

whereMis the mass of the ion, andωXis the trap frequency along theX-axis. The external forcef(t), which can be applied by an electric field or by properly adjusted laser beams, shifts the trap center by?f(t)/Mω2Xand reduces the minimum energy byf2(t)/2Mω2X. When the forcef(t) is applied adiabatically, the final state distribution is conserved in the new basis and remains in thermal equilibrium. If the force is instantly applied to the same maximum value, the final states are highly excited, which reveals farfrom equilibrium dynamics. The Jarzynski equality in the quantum regime can be verified when the average of the exponentiated work〈exp(?W/kBT)〉, which is measured by∑P(W)exp(?W/kBT),is independent of any protocol of applying the work including both cases above.

The essential part of the experimental test in the quantum regime is to perform a two-point projective measurement over the energy eigenstates,which are Fock-states for the harmonic oscillator. In detail, the experimental procedure is composed of the following four stages: 1. preparation of thermal state;2. projection to an energy eigenstate; 3. application of work on the eigenstate; 4. measurement of the final phonon distribution. The procedure is repeated to obtain statistically meaningful results.

Table 1 summarizes the test results of the Jarzynski equality applied to the quantum regime and shows that the Jarzynski estimation gives good prediction of free energy within experimental uncertainties under three different protocols, from three different initial temperatures. In the experiment, the free energy differences ?Fare?2.63kBTeff(Teff=316 nK),?2.13kBTeff(Teff=390 nK),and?1.73kBTeff(Teff=480 nK)depending on the initial temperatures. The experimental data clearly demonstrates the validity of Jarzynski equality when other estimations deviate from the ideal values in far-from equilibrium regime.[126]It is noted that the quantum Jarzynski equality has also been tested in an open quantum system under a decohering heat bath for a spin system of a trapped ion.[143]

Table 1. Summary for the experimental test of the quantum Jarzynski equality.

6.3.2. Quantum absorption refrigerator

A quantum absorption refrigerator is the simplest example of a quantum heat machine.It consists of 3 heat reservoirs:the hot, cold, and work. They interact in such a way that the heat flow from the work to the hot reservoir cools down the cold counterpart. In the case of trapped ions, collective motional modes of the three-ion crystal serve as the heat reservoirs,and nonlinear coupling among the normal modes of motion is induced by the anharmonicity of Coulomb interactions between the ions.The non-degenerate parametric coupling between the three modes of motion described earlier is sufficient for an experimental realization of a quantum absorption refrigerator utilizing three modes of motion of trapped ytterbium ions, one axial and two radials, as the heat bodies (Fig. 11).The radial“zig-zag”(cold)mode represents the heat body that is in contact with a cold environment to be refrigerated,while the axial“zig-zag”(hot)mode of higher frequency represents the heat body in contact with the ambient temperature. The third radial ”rocking” (work) mode serves as the heat source that drives the refrigeration of the cold mode, replacing the work reservoir of a conventional refrigerator.

The interaction Hamiltonian in the system of three ions,induced by anharmonicity of the Coulomb repulsion between the ions,has the form[99,144]

Fig. 11. Illustration of heat flow in an absorption refrigerator. Due to the nature of the interaction among the three bodies (Eq. 25), when the work body is at a higher temperature than the hot body, the resulting heat flow causes heat from the cold body to flow to the hot body as well. Thus heating the work body, and hence increasing its temperature, causes the cold body to cool down.

The operation of the refrigerator consists of three major steps. First, the Raman beams cool down all modes to the ground state of the trap and starting from there selectively prepare the desired states in the modes corresponding to refrigerator bodies. Then, with the lasers switched off, the resonant trilinear interaction is switched on for some time by tuning the mode frequencies to satisfy the resonance condition such that energy exchange between the modes is allowed. Finally, the mode frequencies are brought back to the initial values and measurement of the resulting state in one of the modes is performed with the help of the Raman beams.

The refrigeration itself occurs during the second step,and to see how it works,[144]consider how the resonant interaction Hamiltonian(25)redistributes energy between the modes.The work mode(w)can only remove one of its excess thermal phonons by creating a hot mode (h) phonon and simultaneously annihilating a phonon in the cold mode (c). Hence the transfer of energy from (w) to (h) is always accompanied by energy transfer from (c) to (h). This can result in the cooling of the cold mode when the temperature of the work mode is higher than that of the hot mode, and energy tends to flow from the former to the latter. This process is balanced by the flow of the energy in opposite direction at some temperature,leading to an equilibrium.

To experimentally demonstrate the equilibrium performance of the refrigerator, all modes are prepared in various thermal states. The system is then allowed to evolve for long interaction timeτ ?, and the mean phonon number of each mode is measured to get an estimate of the asymptotic steady state energy.

This experiment also allows exploring quantum effects in the operation of heat machines,by comparing the cooling performance starting from a squeezed thermal state of the work mode[145]to the case where the mode is prepared in a thermal state with the same mean phonon number. The single shot cooling[146,147]method, where the interaction is switched off at the right moment such that the evolution halts at a transient state with a lower mean phonon numberthan the long-time average,was also demonstrated.

6.4. Vibronic sampling

Because of their complexity and the emergent role of quantum nature, molecules are one of the most demanding quantum systems for quantum computers to simulate. For molecular vibronic spectroscopy,it carries the vibrational transitions between nuclear manifolds belonging to two electronic states.[148,149]Upon an electronic transition,a molecule undergoes structural deformation, vibrational frequency shifts, and rotation of normal modes; within a harmonic approximation to the electronic potential energy surfaces,these are equivalent to the displacement(),squeezing(),and rotation()operations in quantum optics,respectively. These effects were simplified as a Duschinsky linear relation.[150]Doktorovet al.[151]decomposed the relation in terms of the following operators:

The process of molecular vibronic spectroscopy can be understood as a modified boson sampling with Gaussian input states such as thermal and squeezed vacuum states. The Gaussian boson sampling, which is classified as classically hard problem in the computational complexity perspective,[153,154]requires additional quantum optical operations on top of the beam splitting and phase shifting operations for standard boson sampling.

Using a trapped ion, the quantum simulation of molecular vibronic spectroscopy with a particular example of photoelectron spectroscopy of sulfur dioxide (SO2)[155,156]is performed. The two vibrational modes of the molecule are mapped to the two radial phonon modes of an ion. Then,the molecular spectroscopy is simulated through the following procedure: (i)the ion is first initialized in the motional ground state,(ii)the quantum optical operations of Doktorov[151]are sequentially applied,and(iii)finally,the vibronic spectrum is then constructed using the collective projection measurements on the transformed state.

After preparing the motional ground-state by Doppler cooling and resolved sideband cooling methods,[58,59]Doktorov operations are performed as the required displacement,squeezing, and rotation operations with the conversion of molecular parameters to the corresponding device parameters asα,ζ′,U,andζ,respectively.The molecular parameters can be obtained via conventional quantum chemical calculations with available program packages(e.g.,Ref.[157]).

Fig. 12. The simulation result of photoelectron spectra of SO2 and SO?2.The two vibrational frequencies of harmonic potential for SO+2, SO2, and SO?2 are(1112.7,415.0),(1178.4,518.9),and(989.5,451.4)cm?1,respectively. (a)The displacement vector α is(?0.026,1.716);rotation angle θ is 0.189.(b)The displacement vector α is(1.360,?0.264);rotation angle θ is 0.065. The results of theoretical calculations are intentionally broadened by convolution with a Gaussian function having a width of 50 cm?1[156]for the comparison. Here Ni0 denotes the i-phonon excitation on N-th mode from the vibrational ground state |0〉, and accordingly, 000 locates at the off-set energy ω0?0=0. Figure adapted from Ref.[30].

As discussed in the previous sections, the quantum optical operations (displacement, squeezing, and rotation)are implemented by Raman laser beams. The spectrum at zero Kelvin from the output measurements of the trapped-ion simulator,reconstructed from the ground to excited states transition intensities,are aligned according to the transition frequencies.Finally,the collective quantum-projection measurement of the final state|nX,nY〉is performed.[27,29]

7. Conclusion and outlook

In conclusion, motional degrees of freedom of trapped ions are valuable resources that give access to the large Hilbert space encoded in bosonic modes.Ultimately,we do not expect any fundamental problems of increasing the number of vibrational modes, which can be manipulated in quantum regime,to the order of hundreds with a few phonons per mode in the transverse direction. In a few examples discussed above,this already leads to various applications in quantum simulations, quantum thermodynamics, and even quantum computations. In addition, the quantum simulation that involves phonon states, or combination of the phonon and internal states of the ions[158–162]has addressed various physics problems, including even the problems of relativistic quantum mechanics,[163,164]quantum field theory,[165,166]and others.

However, several challenges in control of the motional states, both technical and conceptual, currently hinder experimental progress in this area. First limitation in the ion trap experiments comes from detection of motional states,especially when measurement and feedback protocols are employed. Typically, the state of motion in trapped ion experiment is detected by coupling the motional and internal states of the ion,and subsequent fluorescent detection of the ion internal state. Unfortunately during the latter procedure, scattering of thousands or even millions of photons occurs when the internal state is a “bright state”. Most of the time recoil destroys the motional state of all the modes and the experiment has to be restarted. However, it is sometimes possible to harness the recoil provided by the scattered photons to implement feedback protocols.[167]In addition,even though the motional state is characterised by continuous variables that describe the position and momentum of the ions,the information extracted though the measurement is mapped to internal states of the ions and corresponds to discrete variables.[70]Large and efficient collections of scattered photons can also enable us to perform the repetitive detection of phonon number states before the number states are perturbed.[26,27]

Anomalous heating and dephasing of the modes typically limit coherence time of harmonic oscillator for the|0〉+|1〉superpositions to around 10 ms,with dephasing being the dominant decoherence mechanism in most of the experiments. As the impact of both mechanisms increases with the phonon number, this limits the Hilbert space dimension that one can experimentally control,and the number of useful quantum operations that can be performed. Using rf source with stable amplitude for driving the trap electrodes can help to mitigate dephasing.[168]Cleaning of the trap before assembly[169]and placing the trap into a cryogenic environment[170]were shown to reduce heating of the modes.In addition,the effects of heating are typically significantly lower for non-center-of-mass modes.Furthermore,the fact that fluctuation of the radial trapping potential affects some modes in a similar manner, with induced errors cancel each other when these modes are manipulated together, gives rise to an extended coherence time compared to that of the individual modes.[78,102]

In the future, as the number of simultaneously controllable modes increases, more advanced applications will be possible. These include universal quantum computation with continuous variables,various machine learning algorithms,[20]such as matrix inversion, principal component analysis, and implementations of support vector machines.[171]By exploiting a greater number of modes, quantum simulations stand to benefit as well, which includes, for instance, simulations of molecular vibronic spectra and quantum simulation of polaron together with internal states. Ultimately,increase in the number of modes that can be controlled simultaneously may eventually lead to Hilbert spaces whose dimensions exceed the simulation ability of classical computers.