Recurrences in an isolated quantum many-body system

See allHide authors and affiliations

Science  20 Apr 2018:
Vol. 360, Issue 6386, pp. 307-310
DOI: 10.1126/science.aan7938

Recurring coherence

A finite isolated system should return almost to its initial state if it evolves for long enough. For a large system, “long enough” is often unfeasibly long. Rauer et al. found just the right conditions to observe the recurrence of the initial state in a system of two one-dimensional superfluids with thousands of atoms in each. The superfluids were initially coupled—locking their quantum mechanical phases together—and then allowed to evolve independently. After the uncoupling, the researchers observed their phases regaining coherence two more times.

Science, this issue p. 307


The complexity of interacting quantum many-body systems leads to exceedingly long recurrence times of the initial quantum state for all but the smallest systems. For large systems, one cannot probe the full quantum state in all its details. Thus, experimentally, recurrences can only be determined on the level of the accessible observables. Realizing a commensurate spectrum of collective excitations in one-dimensional superfluids, we demonstrate recurrences of coherence and long-range order in an interacting quantum many-body system containing thousands of particles. Our findings will enable the study of the coherent dynamics of large quantum systems even after they have reached a transient thermal-like state.

The expectation that a nonequilibrium system evolves toward thermal equilibrium is deeply rooted in our daily experience and is one of the foundations of statistical mechanics (1). However, as formulated by Poincaré and Zermelo, a finite isolated physical system will recur arbitrarily close to its initial state after a sufficiently long but finite time (2, 3). The reconciliation of these seemingly contradicting statements is at the heart of the emergence of irreversible processes from reversible microscopic mechanics (4).

The above discussion can be transferred to the quantum domain. In analogy to Boltzmann’s conjecture, von Neumann formulated a quasi-ergodic theorem for the evolution of the wave function (5), and the equilibration of isolated quantum systems grew into an active field of research (6). However, a general recurrence theorem can be proven in quantum mechanics as well (7, 8) and shows explicitly that the wave function returns arbitrarily close to its initial state. Experimentally, recurrent behavior has been observed only in small quantum systems. Examples are the revival dynamics in the Jaynes-Cummings model (9) or few interacting atoms in isolated sites of an optical lattice (10, 11). These studies accomplished coherent dynamics over long times and thereby revealed properties of the underlying Hamiltonian spectrum. Enabling recurrent behavior in more complex quantum many-body systems opens a window into their long-term coherent evolution by providing a strong measurable signal at times far beyond the initial relaxation. Further, it is of fundamental interest as it provides insight into the emergence of statistical ensembles from coherent unitary evolution.

For larger systems, however, it becomes exponentially difficult to observe the eigenstates directly, and the complexity of the spectrum of eigenstates leads to exceedingly long recurrence times, in general prohibiting their observation. Nevertheless, for a large class of systems, the essential dynamical features can be described by effective field theories with a much simpler structure (12, 13). This reflects the distinction between relevant and irrelevant operators in the microscopic model and sharply reduces the complexity of the problem from a large number of constituents to a much smaller number of populated modes. Within this level of description, quantum recurrences manifest through the return of observables dominated by these modes. Designing the whole system such that the collective excitations follow a simple commensurate spectrum makes the observation of recurrences feasible, even for many-body systems containing thousands of interacting particles.

Ultracold gases (14) are an ideal starting point for studying these fundamental phenomena as they can be well isolated from the environment, and excellent tools are available to prepare, manipulate, and probe them.

As a model system, we study coherence in one-dimensional (1D) superfluids. In a first approximation, the corresponding many-body physics can be mapped to an effective low-energy description: a bosonic Luttinger liquid (1518) where the collective excitations are free phonons. These phonons are directly related to the phase fluctuations observed in the interference of two 1D superfluids (19).

The dephasing of these collective phononic excitations leads to a loss of coherence and long-range order, proceeding in a light-cone–like fashion (2023). Whereas the fully dephased state is described by a generalized Gibbs ensemble (24), the long time behavior crucially depends on the spectrum of the collective phononic modes. In a harmonic longitudinal confinement with trap frequency Embedded Image, these phonon frequencies are noncommensurate: Embedded Image (25), with j being the mode index. If the atoms are confined to a box-shaped trap, the phonon frequencies become commensurate with Embedded Image, and recurrences should be observable at short times (26).

We implement a boxlike confinement in an atom chip double-well setup (27) by adding hard walls to the very weak longitudinal harmonic confinement with the help of a blue-detuned optical dipole potential (28). To observe recurrent dynamics, we prepare a thermal equilibrium state in a coupled double-well potential (29). Typical samples have a linear density of about 70 atoms per micrometer and, depending on the box size, between 2300 and 4800 atoms in each well, resulting in an interaction energy per atom close to 1 kHz. The tunneling coupling Embedded Image is tuned to a regime where the phases Embedded Imageof the two superfluids lock, creating a system with a strongly correlated relative phase field Embedded Image. To initiate the nonequilibrium dynamics, the coupling is rapidly ramped to zero, leaving the two gases to evolve independently (Fig. 1).

Fig. 1 Schematics of the experiment and the measurement process.

(A) Two coupled 1D superfluids Embedded Image in a double-well potential are taken out of equilibrium by a sudden quench of the tunneling coupling Embedded Image to zero. Along the longitudinal direction Embedded Image, the system is confined to a box-shaped potential of variable length Embedded Image. (B) The evolution of the decoupled system is measured through matter-wave interference in time-of-flight. A typical interference picture showing the atomic density Embedded Image is given. A sinusoidal function is fitted to the local interference fringes for each position z. On the left, we show an example of such a fit (red) for one slice (blue), indicated by the dashed line in the image; the Embedded Image coordinate being the spatial direction of the double-well separation. These fits give access to the spatially resolved relative phase Embedded Image between the two superfluids (bottom).

We observe the subsequent dynamics by matter-wave interferometry (27), which gives direct access to the spatially resolved relative phase Embedded Image between the two superfluids. To study the coherence and long-range order, we evaluate the two-point correlation function (23)Embedded Image(1)with the expectation value taken over many experimental realizations. Further, we average over all points Embedded Image within the central part of the clouds that fulfill Embedded Image.

A typical temporal evolution of the phase correlations in a box trap of length Embedded Image is shown in Fig. 2A. Before the quench at Embedded Image, the relative phase between the superfluids is locked and correlations are close to unity over the whole sample. Immediately after decoupling, this long-range order decays, reflecting the dephasing of the collective excitations. At the first minimum in Fig. 2B, the initial state is completely dephased, and the system is indistinguishable from a thermal state (see inset). For a system with incommensurately spaced modes, this dephased state persists for a long time, showing the emergence of statistical properties from the unitary quantum evolution (22). In our system, however, mode frequencies are designed to be commensurate, and two partial recurrences of phase coherence are clearly visible in the subsequent evolution.

Fig. 2 Dynamics after decoupling.

(A) Temporal evolution of the phase correlation function Embedded Image after decoupling in a 49-μm-long box potential. (B) A cut at Embedded Image, with the error bars giving the 68% confidence interval obtained from a bootstrap. The first and second recurrences are clearly visible and occur at the expected times Embedded Image and Embedded Image, as indicated by the vertical blue and red bars, respectively. The distance Embedded Image is chosen such that it is considerably smaller than the size of the system but long enough that the correlations after the initial dephasing are low; this is to facilitate a high recurrence visibility. The inset shows the phase correlation function at the first (blue) and second (red) recurrence and compares them to the correlations in between the recurrences (green) and the correlations in the initial state (orange). The corresponding points in (B) are colored accordingly. A thermal fit to the correlations in between the recurrences is given by the solid black line. (C) Evolution of the coherence factor Embedded Image in the center of the trap at Embedded Image, showing no sign of a recurrence (see main text for details). Owing to the near homogeneity of the system, the behavior is similar over the whole sample. Therefore, in contrast to the phase correlation function, the averaged interference picture shows no recurrent behavior, as shown in the inset (integrated along the z axis).

To understand the recurrence time, we need to look at the dispersion relation of excitations. For a perfect hard-wall box confinement, Embedded Image, where Embedded Image is the speed of sound and j is the mode index (28). These modes are not only commensurate but also equally spaced, facilitating a recurrence at the earliest possible time Embedded Image, with Embedded Image being the energy spacing between the modes. At this time, the lowest-lying mode has finished a full rotation whereas the higher-energy modes have all performed an integer number of turns, bringing all excitations back to their initial configuration. Half-way to this full recurrence, the system rephases to the mirrored initial state. As we initially start from a nearly flat relative phase profile and our observable Embedded Image (Eq. 1) is insensitive to the transformation Embedded Image, this point is equivalent to the full recurrence. Therefore, the expected recurrence time for the correlations is given by Embedded Image (blue and red bars in Fig. 2B), which agrees well with the observed peaks in coherence.

Although two-point correlations return close to the initial state, the relative phase field does not recur. The time evolution of the coherence factor Embedded Image shows no recurrence (Fig. 2C). It relaxes during the initial dephasing dynamics and stays close to zero from then on. This is equivalent to observing that the ensemble-averaged interference picture shows no revival of high-contrast fringes (Fig. 2C, inset). At the time of the recurrences, the interference fringes return close to their straight initial state, displaying long-range coherence throughout the system; however, for each distinct realization, they are shifted by a random global phase. The reason behind this global phase is a small random atom number imbalance between the two wells. This imbalance originates from the thermal fluctuations of the initial state, imperfections in the experiment, and quantum noise relevant for lower temperatures. It leads to an inevitable population of the Embedded Image mode (30), resulting in a global phase accumulation that is different for each realization, and therefore to vanishing interference contrast in the ensemble average. In contrast, the phase correlations (Eq. 1) are insensitive to a global offset of the phase Embedded Image, and the recurrences of excitations on top of the field can be observed. This illustrates that recurrent behavior can be hidden below a global phase diffusion, which necessitates the measurement of, at least, two-point correlations to reveal the underlying coherent dynamics.

To confirm the scaling of the recurrence time with the size of the system, we vary the length of the box potential by changing the position of the dipole trap walls. As the system size is increased, the recurrence is shifted to later times (Fig. 3A). For this comparison, the time axis was rescaled by the theoretical prediction for the speed of sound (28) to make measurements with slightly different atomic densities comparable. Extracting the exact times of the first and second recurrence by fitting the peaks in the data (28) reveals the linear scaling with Embedded Image (Fig. 3B).

Fig. 3 Comparison of the recurrences for different box lengths.

(A) Phase correlations for three different box lengths Embedded Image and Embedded Image (top to bottom). The time axis is rescaled with the theoretical prediction for the speed of sound c. (B) Recurrence time over the box length extracted from the phase correlations at Embedded Image (28). For each box length, the time of the first (blue) and second (red) recurrence is plotted and compared to the ideal linear scaling (dashed lines). The error bars give the 68% confidence interval obtained from a bootstrap, whereas the shaded bars indicate the predictions of the Luttinger liquid model for our particular trap (28). The vertical extension of the bars corresponds to the uncertainty of the decoupling time, whereas the horizontal extension is chosen arbitrarily.

Although the first two recurrences are clearly visible, for all these measurements, the height of the recurrences is rapidly damped (see Fig. 2, A and B), and observing a third recurrence becomes infeasible in most cases. To probe the recurrence decay in more detail, we studied the evolution for different initial temperatures in the Embedded Image trap. Increasing the initial temperature, the height of the observed recurrence with respect to the correlations of the initial state decreases rapidly (Fig. 4). The temperatures for this analysis are extracted from the full distribution functions of interference contrasts (31) for the completely dephased state in between the recurrences (28).

Fig. 4 Temperature dependence of the recurrence height.

The data points represent measurements of the height Embedded Image of the first (blue) and second (red) recurrence in a box of length Embedded Image for different effective temperatures Embedded Image. The height is extracted from fitting the peaks of the correlation function at Embedded Image and the error bars give the 68% confidence interval obtained from a bootstrap. A fit of the full distribution function of measured interference contrasts at Embedded Image gives the effective temperature. The solid lines are results of GPE simulations analyzed in the same way as the experiment. The shaded area indicates the uncertainty caused by the limited experimental statistics (1σ deviation). The dashed lines are the predictions of the Luttinger liquid model for the experimental parameters (28). The experimental shot-to-shot fluctuations of the atom number as well as the inhomogeneous confinement are incorporated in both theoretical calculations and constitute the only source of damping in the Luttinger liquid model.

To understand the origin of this damping with temperature, we considered different theoretical descriptions of our system. We first investigated the low-energy effective description by solving the Luttinger liquid Hamiltonian (dashed lines in Fig. 4). This model describes the free propagation of phononic excitations on top of a stationary background density. For a homogeneous background, it would give perfect recurrences of phase coherence. For the comparison to the experimental data, we consider an inhomogeneous background density that reflects our boxlike potential. In addition, we took into account the typical spread in total particle number and particle number imbalance between the wells, measured independently for the respective samples (28). The shot-to-shot fluctuations of the speed of sound induced by this spread together with the inhomogeneous background density constitute the only source of damping within the model. Although the recurrence times are well described, as seen in Fig. 3B, the damping in this description is too weak to explain the experimental observations even for the lowest temperatures.

As a second model (solid lines in Fig. 4), we numerically simulated the dynamics using the Gross-Pitaevskii equation (GPE) for finite-temperature initial states (28). This description goes beyond the free phononic excitations of the Luttinger liquid and takes interactions between these quasiparticles into account. It agrees well with our experimental findings (Fig. 4). Together with the failure of the Luttinger liquid model, this agreement shows that physics beyond the low-energy description becomes relevant and indicates that phonon-phonon interactions are responsible for the observed damping. As the temperature is increased, these processes become more important as higher-energy modes are populated.

This study illustrates that the observation of recurrences enables insight into the nonequilibrium dynamics of quantum many-body systems (32). Their shape and magnitude are sensitive probes into coherent dynamics far beyond the time scales of the initial relaxation. This feature makes them a versatile tool for the verification of coherence in quantum simulators and to test the regimes of validity for approximate models. Moreover, the long coherent evolution probed through recurrences magnifies the effects of small perturbations whose influence on the short-term dynamics is negligible. In our model case, the decay of recurrences clearly shows that physics beyond the Luttinger liquid model sets in rather fast. A combined study looking at recurrences and mode occupations will shed light onto the fundamental quantum processes in the relaxation of 1D systems.

Supplementary Materials

Materials and Methods

Figs. S1 to S7

References (3350)

Data File S1

References and Notes

  1. See supplementary materials.
Acknowledgments: We acknowledge discussions with I. Mazets, J. Berges, and T. Gasenzer. Funding: This work was supported by the European Union through the EU-FET Proactive grant AQuS, Project no. 640800, and the European Research Council advanced grant QuantumRelax, by the J. Templeton foundation Project no. 60478 and by the Austrian Science Fund (FWF) through the doctoral program CoQuS (W1210) (B.R., T.S., F.C.). This work is supported by the SFB 1225 ‘ISOQUANT’ financed by the German Research Foundation (DFG) and the FWF. Author contributions: B.R., T.S., F.C., and M.T. conducted the experiments and the data analysis. S.E. and B.R. performed the theoretical calculations. J.S. designed the experiment and supervised the project. All authors contributed to interpreting the data and writing the manuscript. Competing interests: None declared. Data and materials availability: The data presented in this paper are available in the supplementary materials.

Stay Connected to Science

Navigate This Article