Isolated quantum many-body systems with integrable dynamics generically do not thermalize when taken far from equilibrium. As one perturbs such systems away from the integrable point, thermalization sets in, but the nature of the crossover from integrable to thermalizing behavior is an unresolved and actively discussed question. We explore this question by studying the dynamics of the momentum distribution function in a dipolar quantum Newton’s cradle consisting of highly magnetic dysprosium atoms. This is accomplished by creating the first one-dimensional Bose gas with strong magnetic dipole-dipole interactions. These interactions provide tunability of both the strength of the integrability-breaking perturbation and the nature of the near-integrable dynamics. We provide the first experimental evidence that thermalization close to a strongly interacting integrable point occurs in two steps: prethermalization followed by near-exponential thermalization. Exact numerical calculations on a two-rung lattice model yield a similar two-timescale process, suggesting that this is generic in strongly interacting near-integrable models. Moreover, the measured thermalization rate is consistent with a parameter-free theoretical estimate, based on identifying the types of collisions that dominate thermalization. By providing tunability between regimes of integrable and nonintegrable dynamics, our work sheds light on the mechanisms by which isolated quantum many-body systems thermalize and on the temporal structure of the onset of thermalization.

Biopolymer translocation is a key step in viral infection processes. The transfer of information-encoding genomes allows viruses to reprogram the cell fate of their hosts. Constituting 96% of all known bacterial viruses [A. Fokine and M. G. Rossmann, Molecular architecture of tailed double-stranded DNA phages, Bacteriophage 4, e28281 (2014)], the tailed bacteriophages deliver their DNA into host cells via an “ejection” process, leaving their protein shells outside of the bacteria; a similar scenario occurs for mammalian viruses like herpes, where the DNA genome is ejected into the nucleus of host cells, while the viral capsid remains bound outside to a nuclear-pore complex. In light of previous experimental measurements of in vivo bacteriophage $\lambda$ ejection, we analyze here the physical processes that give rise to the observed dynamics. We propose that, after an initial phase driven by self-repulsion of DNA in the capsid, the ejection is driven by anomalous diffusion of phage DNA in the crowded bacterial cytoplasm. We expect that this two-step mechanism is general for phages that operate by pressure-driven ejection, and we discuss predictions of our theory to be tested in future experiments.

The electronic and motional degrees of freedom of trapped ions can be controlled and coherently coupled on the level of individual quanta. Assembling complex quantum systems ion by ion while keeping this unique level of control remains a challenging task. For many applications, linear chains of ions in conventional traps are ideally suited to address this problem. However, driven motion due to the magnetic or radio-frequency electric trapping fields sometimes limits the performance in one dimension and severely affects the extension to higher-dimensional systems. Here, we report on the trapping of multiple barium ions in a single-beam optical dipole trap without radio-frequency or additional magnetic fields. We study the persistence of order in ensembles of up to six ions within the optical trap, measure their temperature, and conclude that the ions form a linear chain, commonly called a one-dimensional Coulomb crystal. As a proof-of-concept demonstration, we access the collective motion and perform spectrometry of the normal modes in the optical trap. Our system provides a platform that is free of driven motion and combines advantages of optical trapping, such as state-dependent confinement and nanoscale potentials, with the desirable properties of crystals of trapped ions, such as long-range interactions featuring collective motion. Starting with small numbers of ions, it has been proposed that these properties would allow the experimental study of many-body physics and the onset of structural quantum phase transitions between one- and two-dimensional crystals.

The quantum Rabi model, involving a two-level system and a bosonic field mode, is arguably the simplest and most fundamental model describing quantum light-matter interactions. Historically, due to the restricted parameter regimes of natural light-matter processes, the richness of this model has been elusive in the lab. Here, we experimentally realize a quantum simulation of the quantum Rabi model in a single trapped ion, where the coupling strength between the simulated light mode and atom can be tuned at will. The versatility of the demonstrated quantum simulator enables us to experimentally explore the quantum Rabi model in detail, including a wide range of otherwise unaccessible phenomena, as those happening in the ultrastrong and deep strong-coupling regimes. In this sense, we are able to adiabatically generate the ground state of the quantum Rabi model in the deep strong-coupling regime, where we are able to detect the nontrivial entanglement between the bosonic field mode and the two-level system. Moreover, we observe the breakdown of the rotating-wave approximation when the coupling strength is increased, and the generation of phonon wave packets that bounce back and forth when the coupling reaches the deep strong-coupling regime. Finally, we also measure the energy spectrum of the quantum Rabi model in the ultrastrong-coupling regime.

The eigenstate thermalization hypothesis (ETH) posits that the reduced density matrix for a subsystem corresponding to an excited eigenstate is “thermal.” Here we expound on this hypothesis by asking: For which class of operators, local or nonlocal, is ETH satisfied? We show that this question is directly related to a seemingly unrelated question: Is the Hamiltonian of a system encoded within a single eigenstate? We formulate a strong form of ETH where, in the thermodynamic limit, the reduced density matrix of a subsystem corresponding to a pure, finite energy density eigenstate asymptotically becomes equal to the thermal reduced density matrix, as long as the subsystem size is much less than the total system size, irrespective of how large the subsystem is compared to any intrinsic length scale of the system. This allows one to access the properties of the underlying Hamiltonian at arbitrary energy densities (or temperatures) using just a single eigenstate. We provide support for our conjecture by performing an exact diagonalization study of a nonintegrable 1D quantum lattice model with only energy conservation. In addition, we examine the case in which the subsystem size is a finite fraction of the total system size, and we find that, even in this case, many operators continue to match their canonical expectation values, at least approximately. In particular, the von Neumann entanglement entropy equals the thermal entropy as long as the subsystem is less than half the total system. Our results are consistent with the possibility that a single eigenstate correctly predicts the expectation values of all operators with support on less than half the total system, as long as one uses a microcanonical ensemble with vanishing energy width for comparison. We also study, both analytically and numerically, a particle-number conserving model at infinite temperature that substantiates our conjectures.

Collective migration dominates many phenomena, from cell movement in living systems to abiotic self-propelling particles. Focusing on the early stages of tumor evolution, we enunciate the principles involved in cell dynamics and highlight their implications in understanding similar behavior in seemingly unrelated soft glassy materials and possibly chemokine-induced migration of $\mathrm{CD}{8}^{+}\mathrm{T}$ cells. We performed simulations of tumor invasion using a minimal three-dimensional model, accounting for cell elasticity and adhesive cell-cell interactions, as well as cell birth and death, to establish that cell-growth-rate-dependent tumor expansion results in the emergence of distinct topological niches. Cells at the periphery move with higher velocity perpendicular to the tumor boundary, while the motion of interior cells is slower and isotropic. The mean-square displacement $\mathrm{\Delta }(t)$ of cells exhibits glassy behavior at times comparable to the cell cycle time, while exhibiting superdiffusive behavior, $\mathrm{\Delta }(t)\approx t^{\alpha }$ ($\alpha >1$), at longer times. We derive the value of $\alpha \approx 1.33$ using a field theoretic approach based on stochastic quantization. In the process, we establish the universality of superdiffusion in a class of seemingly unrelated nonequilibrium systems. Superdiffusion at long times arises only if there is an imbalance between cell birth and death rates. Our findings for the collective migration, which also suggest that tumor evolution occurs in a polarized manner, are in quantitative agreement with in vitro experiments. Although set in the context of tumor invasion, the findings should also hold in describing the collective motion in growing cells and in active systems, where creation and annihilation of particles play a role.

Capacitive mixing (CapMix) and capacitive deionization (CDI) are currently developed as alternatives to membrane-based processes to harvest blue energy—from salinity gradients between river and sea water—and to desalinate water—using charge-discharge cycles of capacitors. Nanoporous electrodes increase the contact area with the electrolyte and hence, in principle, also the performance of the process. However, models to design and optimize devices should be used with caution when the size of the pores becomes comparable to that of ions and water molecules. Here, we address this issue by simulating realistic capacitors based on aqueous electrolytes and nanoporous carbide-derived carbon (CDC) electrodes, accounting for both their complex structure and their polarization by the electrolyte under applied voltage. We compute the capacitance for two salt concentrations and validate our simulations by comparison with cyclic voltammetry experiments. We discuss the predictions of Debye-Hückel and Poisson-Boltzmann theories, as well as modified Donnan models, and we show that the latter can be parametrized using the molecular simulation results at high concentration. This then allows us to extrapolate the capacitance and salt adsorption capacity at lower concentrations, which cannot be simulated, finding a reasonable agreement with the experimental capacitance. We analyze the solvation of ions and their confinement within the electrodes—microscopic properties that are much more difficult to obtain experimentally than the electrochemical response but very important to understand the mechanisms at play. We finally discuss the implications of our findings for CapMix and CDI, both from the modeling point of view and from the use of CDCs in these contexts.

Many complex systems are modular. Such systems can be represented as “component systems,” i.e., sets of elementary components, such as LEGO bricks in LEGO sets. The bricks found in a LEGO set reflect a target architecture, which can be built following a set-specific list of instructions. In other component systems, instead, the underlying functional design and constraints are not obvious a priori, and their detection is often a challenge of both scientific and practical importance, requiring a clear understanding of component statistics. Importantly, some quantitative invariants appear to be common to many component systems, most notably a common broad distribution of component abundances, which often resembles the well-known Zipf’s law. Such “laws” affect in a general and nontrivial way the component statistics, potentially hindering the identification of system-specific functional constraints or generative processes. Here, we specifically focus on the statistics of shared components, i.e., the distribution of the number of components shared by different system realizations, such as the common bricks found in different LEGO sets. To account for the effects of component heterogeneity, we consider a simple null model, which builds system realizations by random draws from a universe of possible components. Under general assumptions on abundance heterogeneity, we provide analytical estimates of component occurrence, which quantify exhaustively the statistics of shared components. Surprisingly, this simple null model can positively explain important features of empirical component-occurrence distributions obtained from large-scale data on bacterial genomes, LEGO sets, and book chapters. Specific architectural features and functional constraints can be detected from occurrence patterns as deviations from these null predictions, as we show for the illustrative case of the “core” genome in bacteria.

We address the question of whether the super-Heisenberg scaling for quantum estimation is indeed realizable. We unify the results of two approaches. In the first one, the original system is compared with its copy rotated by the parameter-dependent dynamics. If the parameter is coupled to the one-body part of the Hamiltonian, the precision of its estimation is known to scale at most as $N^{-1}$ (Heisenberg scaling) in terms of the number of elementary subsystems used $N$. The second approach compares the overlap between the ground states of the parameter-dependent Hamiltonian in critical systems, often leading to an apparent super-Heisenberg scaling. However, we point out that if one takes into account the scaling of time needed to perform the necessary operations, i.e., ensuring adiabaticity of the evolution, the Heisenberg limit given by the rotation scenario is recovered. We illustrate the general theory on a ferromagnetic Heisenberg spin chain example and show that it exhibits such super-Heisenberg scaling of ground-state fidelity around the critical value of the parameter (magnetic field) governing the one-body part of the Hamiltonian. Even an elementary estimator represented by a single-site magnetization already outperforms the Heisenberg behavior providing the $N^{-1.5}$ scaling. In this case, Fisher information sets the ultimate scaling as $N^{-1.75}$, which can be saturated by measuring magnetization on all sites simultaneously. We discuss universal scaling predictions of the estimation precision offered by such observables, both at zero and finite temperatures, and support them with numerical simulations in the model. We provide an experimental proposal of realization of the considered model via mapping the system to ultracold bosons in a periodically shaken optical lattice. We explicitly derive that the Heisenberg limit is recovered when the time needed for preparation of quantum states involved is taken into account.

We study the dynamics of a supersonically expanding, ring-shaped Bose-Einstein condensate both experimentally and theoretically. The expansion redshifts long-wavelength excitations, as in an expanding universe. After expansion, energy in the radial mode leads to the production of bulk topological excitations—solitons and vortices—driving the production of a large number of azimuthal phonons and, at late times, causing stochastic persistent currents. These complex nonlinear dynamics, fueled by the energy stored coherently in one mode, are reminiscent of a type of “preheating” that may have taken place at the end of inflation.

The collective plasmonic modes of a metal comprise a simple pattern of oscillating charge density that yields enhanced light-matter interaction. Here we unveil that beneath this familiar facade plasmons possess a hidden internal structure that fundamentally alters its dynamics. In particular, we find that metals with nonzero Hall conductivity host plasmons with an intricate current density configuration that sharply departs from that of ordinary zero Hall conductivity metals. This nontrivial internal structure dramatically enriches the dynamics of plasmon propagation, enabling plasmon wave packets to acquire geometric phases as they scatter. At boundaries, these phases accumulate allowing plasmon waves that reflect off to experience a nonreciprocal parallel shift. This plasmon Hall shift, tunable by Hall conductivity as well as plasmon wavelength, displaces the incident and reflected plasmon trajectories and can be readily probed by near-field photonics techniques. Anomalous plasmon geometric phases dramatically enrich the nanophotonics toolbox, and yield radical new means for directing plasmonic beams.

Roughly every 2–10 min, a pair of stellar-mass black holes merge somewhere in the Universe. A small fraction of these mergers are detected as individually resolvable gravitational-wave events by advanced detectors such as LIGO and Virgo. The rest contribute to a stochastic background. We derive the statistically optimal search strategy (producing minimum credible intervals) for a background of unresolved binaries. Our method applies Bayesian parameter estimation to all available data. Using Monte Carlo simulations, we demonstrate that the search is both “safe” and effective: it is not fooled by instrumental artifacts such as glitches and it recovers simulated stochastic signals without bias. Given realistic assumptions, we estimate that the search can detect the binary black hole background with about 1 day of design sensitivity data versus $\approx 40$ months using the traditional cross-correlation search. This framework independently constrains the merger rate and black hole mass distribution, breaking a degeneracy present in the cross-correlation approach. The search provides a unified framework for population studies of compact binaries, which is cast in terms of hyperparameter estimation. We discuss a number of extensions and generalizations, including application to other sources (such as binary neutron stars and continuous-wave sources), simultaneous estimation of a continuous Gaussian background, and applications to pulsar timing.

Nonlocality and contextuality are at the root of conceptual puzzles in quantum mechanics, and they are key resources for quantum advantage in information-processing tasks. Bell nonlocality is best understood as the incompatibility between quantum correlations and the classical theory of causality, applied to relativistic causal structure. Contextuality, on the other hand, is on a more controversial foundation. In this work, I provide a common conceptual ground between nonlocality and contextuality as violations of classical causality. First, I show that Bell inequalities can be derived solely from the assumptions of no signaling and no fine-tuning of the causal model. This removes two extra assumptions from a recent result from Wood and Spekkens and, remarkably, does not require any assumption related to independence of measurement settings—unlike all other derivations of Bell inequalities. I then introduce a formalism to represent contextuality scenarios within causal models and show that all classical causal models for violations of a Kochen-Specker inequality require fine-tuning. Thus, the quantum violation of classical causality goes beyond the case of spacelike-separated systems and already manifests in scenarios involving single systems.

Unitary transformations are routinely modeled and implemented in the field of quantum optics. In contrast, nonunitary transformations, which can involve loss and gain, require a different approach. In this work, we present a universal method to deal with nonunitary networks. An input to the method is an arbitrary linear transformation matrix of optical modes that does not need to adhere to bosonic commutation relations. The method constructs a transformation that includes the network of interest and accounts for full quantum optical effects related to loss and gain. Furthermore, through a decomposition in terms of simple building blocks, it provides a step-by-step implementation recipe, in a manner similar to the decomposition by Reck et al. [Experimental Realization of Any Discrete Unitary Operator, Phys. Rev. Lett. 73, 58 (1994)] but applicable to nonunitary transformations. Applications of the method include the implementation of positive-operator-valued measures and the design of probabilistic optical quantum information protocols.

We report an experimental demonstration of the nonclassical stationary light pulse (SLP) in a cold atomic ensemble. A single collective atomic excitation is created and heralded by detecting a Stokes photon in the spontaneous Raman scattering process. The heralded single atomic excitation is converted into a single stationary optical excitation or the single-photon SLP, whose effective group velocity is zero, effectively forming a trapped single-photon pulse within the cold atomic ensemble. The single-photon SLP is then released from the atomic ensemble as an anti-Stokes photon after a specified trapping time. The second-order correlation measurement between the Stokes and anti-Stokes photons reveals the nonclassical nature of the single-photon SLP. Our work paves the way toward quantum nonlinear optics without a cavity.

As a fundamental phenomenon of fluid mechanics, recent studies suggested laminar-turbulent transition belonging to the universality class of directed percolation. Here, the onset of a laminar separation bubble on an airfoil is analyzed in terms of the directed percolation model using particle image velocimetry data. Our findings indicate a clear significance of percolation models in a general flow situation beyond fundamental ones. We show that our results are robust against fluctuations of the parameter, namely, the threshold of turbulence intensity, that maps velocimetry data into binary cells (turbulent or laminar). In particular, this percolation approach enables the precise determination of the transition point of the laminar separation bubble, an important problem in aerodynamics.

Random quantum circuits yield minimally structured models for chaotic quantum dynamics, which are able to capture, for example, universal properties of entanglement growth. We provide exact results and coarse-grained models for the spreading of operators by quantum circuits made of Haar-random unitaries. We study both $1+1\mathrm{D}$ and higher dimensions and argue that the coarse-grained pictures carry over to operator spreading in generic many-body systems. In $1+1\mathrm{D}$, we demonstrate that the out-of-time-order correlator (OTOC) satisfies a biased diffusion equation, which gives exact results for the spatial profile of the OTOC and determines the butterfly speed $v_B$. We find that in $1+1\mathrm{D}$, the “front” of the OTOC broadens diffusively, with a width scaling in time as $t^{1/2}$. We address fluctuations in the OTOC between different realizations of the random circuit, arguing that they are negligible in comparison to the broadening of the front within a realization. Turning to higher dimensions, we show that the averaged OTOC can be understood exactly via a remarkable correspondence with a purely classical droplet growth problem. This implies that the width of the front of the averaged OTOC scales as $t^{1/3}$ in $2+1\mathrm{D}$ and as $t^{0.240}$ in $3+1\mathrm{D}$ (exponents of the Kardar-Parisi-Zhang universality class). We support our analytic argument with simulations in $2+1\mathrm{D}$. We point out that, in two or higher spatial dimensions, the shape of the spreading operator at late times is affected by underlying lattice symmetries and, in general, is not spherical. However, when full spatial rotational symmetry is present in $2+1\mathrm{D}$, our mapping implies an exact asymptotic form for the OTOC, in terms of the Tracy-Widom distribution. For an alternative perspective on the OTOC in $1+1\mathrm{D}$, we map it to the partition function of an Ising-like statistical mechanics model. As a result of special structure arising from unitarity, this partition function reduces to a random walk calculation which can be performed exactly. We also use this mapping to give exact results for entanglement growth in $1+1\mathrm{D}$ circuits.

Thermalization and scrambling are the subject of much recent study from the perspective of many-body quantum systems with locally bounded Hilbert spaces (“spin chains”), quantum field theory, and holography. We tackle this problem in 1D spin chains evolving under random local unitary circuits and prove a number of exact results on the behavior of out-of-time-ordered commutators (OTOCs) and entanglement growth in this setting. These results follow from the observation that the spreading of operators in random circuits is described by a “hydrodynamical” equation of motion, despite the fact that random unitary circuits do not have locally conserved quantities (e.g., no conserved energy). In this hydrodynamic picture, quantum information travels in a front with a “butterfly velocity” $v_{\mathrm{B}}$ that is smaller than the light-cone velocity of the system, while the front itself broadens diffusively in time. The OTOC increases sharply after the arrival of the light cone, but we do not observe a prolonged exponential regime of the form $\sim e^{{\lambda }_{\mathrm{L}}(t-x/v)}$ for a fixed Lyapunov exponent ${\lambda }_{\mathrm{L}}$. We find that the diffusive broadening of the front has important consequences for entanglement growth, leading to an entanglement velocity that can be significantly smaller than the butterfly velocity. We conjecture that the hydrodynamical description applies to more generic Floquet ergodic systems, and we support this idea by verifying numerically that the diffusive broadening of the operator wavefront also holds in a more traditional nonrandom Floquet spin chain. We also compare our results to Clifford circuits, which have less rich hydrodynamics and consequently trivial OTOC behavior, but which can nevertheless exhibit linear entanglement growth and thermalization.

We generate and characterize entangled states of a register of 20 individually controlled qubits, where each qubit is encoded into the electronic state of a trapped atomic ion. Entanglement is generated amongst the qubits during the out-of-equilibrium dynamics of an Ising-type Hamiltonian, engineered via laser fields. Since the qubit-qubit interactions decay with distance, entanglement is generated at early times predominantly between neighboring groups of qubits. We characterize entanglement between these groups by designing and applying witnesses for genuine multipartite entanglement. Our results show that, during the dynamical evolution, all neighboring qubit pairs, triplets, most quadruplets, and some quintuplets simultaneously develop genuine multipartite entanglement. Witnessing genuine multipartite entanglement in larger groups of qubits in our system remains an open challenge.

Information-theoretic approaches provide a promising avenue for extending the laws of thermodynamics to the nanoscale. Here, we provide a general fundamental lower limit, valid for systems with an arbitrary Hamiltonian and in contact with any thermodynamic bath, on the work cost for the implementation of any logical process. This limit is given by a new information measure—the coherent relative entropy—which accounts for the Gibbs weight of each microstate. The coherent relative entropy enjoys a collection of natural properties justifying its interpretation as a measure of information and can be understood as a generalization of a quantum relative entropy difference. As an application, we show that the standard first and second laws of thermodynamics emerge from our microscopic picture in the macroscopic limit. Finally, our results have an impact on understanding the role of the observer in thermodynamics: Our approach may be applied at any level of knowledge—for instance, at the microscopic, mesoscopic, or macroscopic scales—thus providing a formulation of thermodynamics that is inherently relative to the observer. We obtain a precise criterion for when the laws of thermodynamics can be applied, thus making a step forward in determining the exact extent of the universality of thermodynamics and enabling a systematic treatment of Maxwell-demon-like situations.