## Abstract

Although statistical mechanics describes thermal equilibrium states, these states may or may not emerge dynamically for a subsystem of an isolated quantum many-body system. For instance, quantum systems that are near-integrable usually fail to thermalize in an experimentally realistic time scale, and instead relax to quasi-stationary prethermal states that can be described by statistical mechanics, when approximately conserved quantities are included in a generalized Gibbs ensemble (GGE). We experimentally study the relaxation dynamics of a chain of up to 22 spins evolving under a long-range transverse-field Ising Hamiltonian following a sudden quench. For sufficiently long-range interactions, the system relaxes to a new type of prethermal state that retains a strong memory of the initial conditions. However, the prethermal state in this case cannot be described by a standard GGE; it rather arises from an emergent double-well potential felt by the spin excitations. This result shows that prethermalization occurs in a broader context than previously thought, and reveals new challenges for a generic understanding of the thermalization of quantum systems, particularly in the presence of long-range interactions.

## INTRODUCTION

In the classical world, thermalization is expected in all but special cases where conserved quantities or hidden symmetries prevent the ergodic exploration of phase space. Because the classical world is ultimately composed of quantum systems, we therefore expect that a closed quantum system will also reach thermal equilibrium (*1*–*9*). Although quantum dynamics are unitary, measurements made within a subsystem trace over the rest of the system and appear thermal because the rest of the system acts as a thermal bath (*10*–*13*).

However, this is not always the case. For integrable models, an extensive number of conserved quantities prevent the efficient exploration of phase space (*14*), and the system relaxes to a steady state predicted by a generalized Gibbs ensemble (GGE) (*4*, *15*–*18*) specified by the initial values of the integrals of motion. For near-integrable systems, such as weakly interacting ultracold gases, thermalization can still occur but only over extremely long time scales beyond current experimental reach (*16*, *17*). However, it is possible to observe quasi-stationary states, often called prethermal (*19*), which emerge within an experimentally accessible time scale.

Previous observations of prethermal states have focused on those described by a GGE associated with the integrable part of the model (*16*, *17*). Here, we observe a different type of prethermalization (*20*), where a system of interacting spins rapidly evolves to a quasi-stationary state that cannot be predicted by a standard GGE. This type of prethermal state arises, even in the thermodynamic limit, when a system has long-range interactions and open boundaries such that the translational invariance is broken. As a result, spin excitations feel an emergent double-well potential whose depth grows with interaction range (Fig. 1). Memory of the initial state is preserved by this emergent potential but is eventually lost because of weak quantum tunneling between the two wells and the interactions between spin excitations.

## RESULTS

Effective spin-1/2 particles are encoded in the ^{2}S_{1/2} |*F* = 0, *m*_{F} = 0〉 and |*F* = 1, *m*_{F} = 0〉 hyperfine “clock” states of a ^{171}Yb^{+} ion, denoted | ↓ 〉_{z} and | ↑ 〉_{z} (*21*). We confine a chain of ions in a linear radio frequency Paul trap and apply optical dipole forces to generate the effective spin-spin coupling (*22*, *23*) of an Ising Hamiltonian(1)where are the Pauli matrices acting on the *i*th spin, *J*_{ij} is the coupling between spins *i* and *j*, and *B*/(2π) = 10 kHz is a uniform effective transverse field. We use units in which Planck’s constant equals 1. The spin-spin interaction is long-range and can be described by a power-law decay, , where *J*_{max} is the maximum coupling strength, which ranges from 0.45 to 0.98 kHz. We tune the power-law exponent α between 0.55 and 1.33 by changing the axial confinement of the ions. With long-range interactions, *H* is, in general, nonintegrable [in contrast to the nearest-neighbor case where the one-dimensional (1D) model is integrable (*24*)], and thermalization is anticipated in the long-time limit according to the eigenstate thermalization hypothesis (*25*–*28*). However, we can map each spin excitation along the *z* direction into a bosonic particle to turn Eq. 1 into a bosonic model with two parts: An integrable part made from noninteracting bosons that will be used to construct a GGE, and an integrability-breaking part consisting of interactions among the bosons, which is responsible for the eventual thermalization (see the Supplementary Materials for details). When the initial state has a low spin/bosonic excitation density, we expect the bosonic model to be near-integrable because the interactions are weak.

We initialize the chain of seven spins by optically pumping all spins to the |↓〉_{z} state, and then use a tightly focused individual-ion addressing laser to excite a single spin on one end of the chain to the | ↑ 〉_{z} state as seen in Fig. 1A (*29*). The spins then evolve under Eq. 1, and we measure the time evolution of the spin projection in the *z*-basis. For the shortest-range interactions we realize (α = 1.33), the system rapidly evolves to a prethermal state predicted by the GGE associated with the integrals of motion corresponding to the eigenmode occupation numbers of the noninteracting bosons, which does not preserve memory of the initial spin excitation location (Fig. 1B).

However, in the long-range interacting case (α = 0.55), we see that the position of the spin excitation reaches an equilibrium value that retains a memory of the initial state (Fig. 1D) out to the longest experimentally achievable time of 25/*J*_{max}. This prethermal state is in obvious disagreement with both a thermal state and the GGE prediction, which both maintain the right-left symmetry of the system.

The dynamics of the spin-wave boson model for short-range interactions are similar to those of a free particle in a square-well potential. However, the presence of long-range interactions distorts the square-well into a double-well potential (Fig. 1A). Here, we emphasize that the double-well potential emerges, despite the fact that our model (Eq. 1) is transitionally invariant without boundaries. The spin-wave boson model has an extensive number of near-degenerate eigenstates that are symmetric and antisymmetric superpositions of spin excitations in the left and right potential wells. For seven spins, we calculate the energy difference, Δ*E*_{mn}, between all pairs of eigenstates and plot them with respect to α in Fig. 1C, where the amplitudes, |〈*m*|ψ_{0}〉〈ψ_{0}|*n*〉|^{2}, are products of the overlaps of the eigenstates |*m*〉 and |*n*〉 with the initial state |ψ_{0}〉. With α = 0.55, the two lowest energy states are almost degenerate, with an energy difference approximately 1000 times smaller than *J*_{max}. This is due to the tunneling rate between the ground states of each well, which is exponentially small in the barrier height. As a result, the spin excitation will remain in its initial well until it tunnels across the potential barrier at much longer times. Note that even for α = 1.33 (Fig. 1B), a small potential still exists, but tunneling already happens during the experiment.

To better characterize the “location” of the spin excitation, we construct the observable(2)where *N* is the number of ions. The expectation value of *C* varies between −1 and 1 for a spin excitation on the left and right ends, respectively. Because of the spatial inversion symmetry of the spin-wave boson model and Eq. 1, both the GGE predicted and thermal values of 〈*C*〉 are zero. In Fig. 2, we plot 〈*C*〉 along with its cumulative time average for α = 1.33 (short range) and α = 0.55 (long range). This averaging smooths out fast temporal fluctuations in our small systems. To further accentuate the asymmetry of the prethermalization, we prepare initial states with a single-spin excitation on the left or right ends of the chain.

With short-range interactions, the system quickly relaxes to the prethermal value predicted by the GGE, where 〈*C*〉 is near zero irrespective of the initial state. But with long-range interactions, the prethermal state that emerges retains a clear memory of the initial conditions and is different than both the thermal and GGE predictions.

In addition to showing the relaxation behavior of the global observable 〈*C*〉, we also observe the dynamics of local observables with single-site resolution. In Fig. 2, we plot the cumulative time average of the deviation of the single-spin magnetizations, , from the equilibrium value predicted by the GGE. Here, we postselect the data for the correct number of spin excitations to eliminate the effects of imperfect state preparation and small deviations from our model Hamiltonian due to unwanted excitations of the phonon modes (see Supplementary Materials). For the short-range interactions, we see that the cumulative time average for each spin quickly converges to the GGE predicted value. However, with the long-range interactions, we see that the individual spins reach a steady state that does not match the GGE because the emergent double well prevents the efficient transfer of the excitation across the chain.

We show the prethermal state’s robustness to weak interactions, similar to how many-body localization is robust to interactions (*7*, *30*), by preparing initial states with two spin excitations. In this case, the multiple spin flips increase the size of the integrability-breaking part of the Hamiltonian, which represents weak interactions between the spin-wave bosons. However, the prethermal state still persists. For better contrast between the prethermal and GGE predicted values of 〈*C*〉, we flip the second and fourth spins such that |ψ_{0}〉 = | ↓↑↓↑↓↓↓ 〉. We emphasize that the observed prethermalization and departure from the GGE prediction are not sensitive to the specific choice of initial state in the thermodynamic limit, as long as the initial state is not symmetric. But for a small-sized system, these initial states offer us the maximum signal. As before, we also prepare the mirror image of the initial state by exciting the fourth and sixth ions ( = | ↓↓↓↑↓↑↓ 〉). We observe relaxation to the value predicted by the GGE for short-range interactions, but with long-range interactions, we see a prethermal state that strongly deviates from the GGE (bottom panel of Fig. 2).

In Fig. 3, we plot the experimental evolution of the prethermal state for both the double and single spin flip initial states along with a long-time numerical simulation under Eq. 1, accounting for known experimental noise. The plots demonstrate excellent agreement between numerical simulations and experimental data and confirm that the prethermal states persist well beyond the current experimental time limit. Because of the nonconservation of the number of spin-wave bosons for any finite *B* and the interactions between them, the system will eventually relax to the thermal equilibrium in the thermodynamic limit; however, relaxation to the GGE may or may not be seen depending on the range of interactions.

To demonstrate that the prethermal state we observe is not sensitive to system size, we repeat the experiment in a chain of 22 ions, the largest ion chain used for quantum simulation in the literature to date. The power-law interaction range is α ≈ 0.9. The time evolution and cumulative time average of 〈*C*〉 are depicted in Fig. 4. Experimentally, as well as in the analytic result in the thermodynamic limit (see the Supplementary Materials), we see that the system relaxes to a similar quasi-equilibrium state as before that clearly has memory of the initial state. This is because the number of near-degenerate pairs of (localized) eigenstates is an extensive quantity, and the observed prethermalization will not disappear as the system size grows.

## DISCUSSION

We point out that the observed prethermalization and deviation from the GGE should disappear if the system is subject to periodic boundary conditions, regardless of system size. Here, we emphasize that the long-range interactions make the boundary conditions relevant for bulk properties, and thus, changing the boundary conditions can affect an extensive number of eigenstates (*31*). The effects of long-range interactions on many-body dynamics are far richer than the effect discussed in the current experiment. For sufficiently long-range interactions, the notion of locality breaks down and quasi-particles in the system can travel at divergent velocities for thermodynamic systems, potentially leading to markedly different thermalization/prethermalization time scales in certain systems (*32*–*34*). We believe that the current experiment, as well as the platform it is built upon, will pave the way to a more complete understanding of the fundamental role long-range interactions play in the quench dynamics and emergent statistical physics of quantum many-body systems.

## MATERIALS AND METHODS

### Effective Hamiltonian generation

We generated spin-spin interactions by applying spin-dependent optical dipole forces to ions confined in a three-layer linear Paul trap with a 4.8-MHz radial frequency. Two off-resonant laser beams with a wave vector difference δ*k* along a principal axis of transverse motion globally address the ions and drive stimulated Raman transitions. The two beams contain a pair of beat note frequencies symmetrically detuned from the resonant transition at ν_{0} = 12.642819 GHz by a frequency μ, comparable to the transverse motional mode frequencies. In the Lamb-Dicke regime, this results in the Ising-type Hamiltonian in Eq. 1 (*22*, *23*, *35*) with(3)where Ω is the global Rabi frequency, *ħ* is the reduced Planck’s constant, *V*_{i,m} is the normal-mode matrix (*36*), and ω_{m} are the transverse mode frequencies. The coupling profile may be approximated as a power-law decay *J*_{ij} ≈ *J*_{max}/|*i* − *j*|^{α}, where, in principle, α can be tuned between 0 and 3 by varying the laser detuning μ or the bandwidth of ω_{m} (*23*, *35*). For the seven-ion data in this work, α was tuned to 0.55 (long-range interactions) and 1.33 (short-range interactions) by changing the bandwidth of ω_{m}. By asymmetrically adjusting the laser beat note detuning μ about the carrier by a value of 2*B*, we apply a uniform effective transverse magnetic field of (*37*).

### Single spin flip initialization and site-resolved detection

We initialized individual spin excitations using a tightly focused laser beam to imprint a fourth-order ac Stark shift (*29*) in conjunction with a Ramsey or Rabi sequence. When the ion spacing is larger than the beam waist of the individual-ion addressing laser as is the case for the seven-ion short-range data, we used a Ramsey method. This consists of first optically pumping the spins to | ↓ 〉_{z}. Then, we globally performed a π/2 rotation so that all of the spins are in | ↓ 〉_{x}. Using the individual-ion addressing beam, we applied a Stark shift to the spins to be flipped, and then we allowed the chain to evolve until these spins are π out of phase compared to the spins without an applied Stark shift. Afterward, a global π/2 rotation brings the spins back into the *z*-basis. With this method, individual spin flips can be prepared with a fidelity of ~0.97, whereas *N* spin flips can be achieved with a fidelity of ~(0.97)^{N}.

We used the Rabi method for the long-range interacting data because site-resolved Stark shifts can no longer be applied as the ion separation is smaller than the beam waist. Here, we applied a large Stark shift to all of the spins, except the ones to be flipped, and a global π pulse at the hyperfine splitting between the two effective spin levels. Thus, only the ions without an applied Stark shift were flipped. This approach has a single and *N* spin flip fidelity of ~0.85 and ~(0.85)^{N}, respectively.

After quenching to and allowing time evolution under our spin Hamiltonian, we measured the spin projections of each ion along the *z* direction of the Bloch sphere. We exposed the ions to a laser beam that addresses the cycling transition ^{2}S_{1/2} |*F* = 1〉 to ^{2}P_{1/2}|*F* = 0〉 for 3 ms. Ions fluoresce only if they are in the state | ↑ 〉_{z}. This fluorescence was collected through an NA (numerical aperture) = 0.23 objective and imaged using an intensified charge-coupled device (CCD) camera with single-site resolution.

To discriminate between “bright” and “dark” states (| ↑ 〉_{z} and | ↓ 〉_{z}, respectively), we begin by calibrating the camera with 1000 cycles each of all-bright and all-dark states. For the bright states, the projection of the 2D CCD image onto a 1D row gives a profile composed of Gaussian distributions at each ion location. We performed fits to locate the center and fluorescence width of each ion.

We achieved single-shot discrimination of individual ion states in the experimental data by fitting the captured 1D profile to a series of Gaussian distributions with calibrated widths and positions but freely varying amplitudes. These extracted values for each ion were then compared with a threshold found via Monte Carlo simulation to determine whether the measured state was bright or dark. Our discrimination protocol also gives an estimate of the detection error (for example, misdiagnosing a bright ion as dark), which is typically of order ~5%. Corrected state probabilities (along with their respective errors) were found following the method outlined by Shen and Duan (*38*), which also takes into account errors due to quantum projection noise.

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/3/8/e1700672/DC1

Experimental noise sources and their influence on the thermalization dynamics

Measuring the spin-spin coupling matrix

Justification for postselection

The spin-boson mapping and the GGE

Single-particle properties of *H*_{0}

Discussion

fig. S1. We directly measure the spin-spin coupling matrix with seven ions for both short-range (left matrix) and long-range (right matrix) interactions and see if it is symmetric.

fig. S2. Numerical calculation to illustrate the origin of the double-well potential.

Reference (*39*)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is **not** for commercial advantage and provided the original work is properly cited.

## REFERENCES AND NOTES

**Acknowledgments:**

**Funding**

**:**This work is supported by the Army Research Office (ARO) Atomic and Molecular Physics Program, the Air Force Office of Scientific Research (AFOSR) Multidisciplinary University Research Initiative (MURI) on Quantum Measurement and Verification, the Intelligence Advanced Research Projects Activity LogiQ program, the NSF Physics Frontier Center at Joint Quantum Institute, NSF Quantum Information Science, AFOSR, Army Research Laboratory Center for Distributed Quantum Information, and ARO MURI.

**Author contributions:**B.N., Z.-X.G., P.R., and C.M. conceived and designed the experiment; J.Z., J.S., A.C.L., and P.W.H. performed the experiment; B.N., Z.-X.G., J.S., and A.C.L. performed numerical simulations; and Z.-X.G. and A.V.G. performed theoretical analysis. All authors contributed to the writing of the manuscript.

**Competing interests:**C.M. is a founding scientist of ionQ Inc. All other authors declare that they have no competing interests.

**Data and materials availability:**All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Original data are available upon reasonable request.

- Copyright © 2017 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).