Absence of localization in interacting spin chains with a discrete symmetry.

Novel paradigms of strong ergodicity breaking have recently attracted significant attention in condensed matter physics. Understanding the exact conditions required for their emergence or breakdown not only sheds more light on thermalization and its absence in closed quantum many-body systems, but it also has potential benefits for applications in quantum information technology. A case of particular interest is many-body localization whose conditions are not yet fully settled. Here, we prove that spin chains symmetric under a combination of mirror and spin-flip symmetries and with a non-degenerate spectrum show finite spin transport at zero total magnetization and infinite temperature. We demonstrate this numerically using two prominent examples: the Stark many-body localization system (Stark-MBL) and the symmetrized many-body localization system (symmetrized-MBL). We provide evidence of delocalization at all energy densities and show that delocalization persists when the symmetry is broken. We use our results to construct two localized systems which, when coupled, delocalize each other. Our work demonstrates the dramatic effect symmetries can have on disordered systems, proves that the existence of exact resonances is not a sufficient condition for delocalization, and opens the door to generalization to higher spatial dimensions and different conservation laws.

Thermalization dynamics of a gauge theory on a quantum simulator.

Gauge theories form the foundation of modern physics, with applications ranging from elementary particle physics and early-universe cosmology to condensed matter systems. We perform quantum simulations of the unitary dynamics of a U(1) symmetric gauge field theory and demonstrate emergent irreversible behavior. The highly constrained gauge theory dynamics are encoded in a one-dimensional Bose-Hubbard simulator, which couples fermionic matter fields through dynamical gauge fields. We investigated global quantum quenches and the equilibration to a steady state well approximated by a thermal ensemble. Our work may enable the investigation of elusive phenomena, such as Schwinger pair production and string breaking, and paves the way for simulating more complex, higher-dimensional gauge theories on quantum synthetic matter devices.

Quantum aging and dynamical universality in the long-range O(N→∞) model.

Quantum quenches to or near criticality give rise to the phenomenon of aging, manifested by glassylike dynamics at short times and far from equilibrium. The recent surge of interest in the dynamics of quantum many-body systems has rejuvenated interest in this phenomenon. Motivated by the ubiquitous long-range interactions in emerging experimental platforms, it is vital to study quantum aging in such settings. In this paper, we investigate the dynamical universality and aging in the d-dimensional O(N) model with the long-range coupling 1/x^{d+σ} and in the mean-field limit N→∞ that allows an exact treatment. An immediate consequence of long-range coupling is the emergence of nonlinear light cones. We focus on the correlation and response functions, and identify a rich scaling behavior depending on how the corresponding space-time positions are located relative to each other, via a local light cone, and to the time of the quench via a global quench light cone. We determine the initial-slip exponent that governs the short-time dependence of two-point functions. We highlight the qualitative features of aging due to the long-range coupling, in particular in the region outside the light cones. As an important consequence of long-range coupling, the correlation function decays as 1/x^{d+σ} outside the quench light cone while increasing polynomially with the total time after quench. This is while, for short-time differences, the two-time response function "equilibrates" at all distances even outside this light cone. Our analytic findings are in excellent agreement with exact numerics, and provide a useful benchmark for modern experimental platforms with long-range interactions.

Observation of gauge invariance in a 71-site Bose-Hubbard quantum simulator.

The modern description of elementary particles, as formulated in the standard model of particle physics, is built on gauge theories1. Gauge theories implement fundamental laws of physics by local symmetry constraints. For example, in quantum electrodynamics Gauss's law introduces an intrinsic local relation between charged matter and electromagnetic fields, which protects many salient physical properties, including massless photons and a long-ranged Coulomb law. Solving gauge theories using classical computers is an extremely arduous task2, which has stimulated an effort to simulate gauge-theory dynamics in microscopically engineered quantum devices3-6. Previous achievements implemented density-dependent Peierls phases without defining a local symmetry7,8, realized mappings onto effective models to integrate out either matter or electric fields9-12, or were limited to very small systems13-16. However, the essential gauge symmetry has not been observed experimentally. Here we report the quantum simulation of an extended U(1) lattice gauge theory, and experimentally quantify the gauge invariance in a many-body system comprising matter and gauge fields. These fields are realized in defect-free arrays of bosonic atoms in an optical superlattice of 71 sites. We demonstrate full tunability of the model parameters and benchmark the matter-gauge interactions by sweeping across a quantum phase transition. Using high-fidelity manipulation techniques, we measure the degree to which Gauss's law is violated by extracting probabilities of locally gauge-invariant states from correlated atom occupations. Our work provides a way to explore gauge symmetry in the interplay of fundamental particles using controllable large-scale quantum simulators.

Reliability of Lattice Gauge Theories.

Currently, there are intense experimental efforts to realize lattice gauge theories in quantum simulators. Except for specific models, however, practical quantum simulators can never be fine-tuned to perfect local gauge invariance. There is thus a strong need for a rigorous understanding of gauge-invariance violation and how to reliably protect against it. As we show through analytic and numerical evidence, in the presence of a gauge invariance-breaking term the gauge violation accumulates only perturbatively at short times before proliferating only at very long times. This proliferation can be suppressed up to infinite times by energetically penalizing processes that drive the dynamics away from the initial gauge-invariant sector. Our results provide a theoretical basis that highlights a surprising robustness of gauge-theory quantum simulators.

Dynamical Quantum Phase Transitions: A Geometric Picture.

The Loschmidt echo is a purely quantum-mechanical quantity whose determination for large quantum many-body systems requires an exceptionally precise knowledge of all eigenstates and eigenenergies. One might therefore be tempted to dismiss the applicability of any approximations to the underlying time evolution as hopeless. However, using the fully connected transverse-field Ising model as an example, we show that this indeed is not the case and that a simple semiclassical approximation to systems well described by mean-field theory is, in fact, in good quantitative agreement with the exact quantum-mechanical calculation. Beyond the potential to capture the entire dynamical phase diagram of these models, the method presented here also allows for an intuitive geometric interpretation of the fidelity return rate at any temperature, thereby connecting the order parameter dynamics and the Loschmidt echo in a common framework. Videos of the postquench dynamics provided in Supplemental Material visualize this new point of view.

Universal equilibrium scaling functions at short times after a quench.

By analyzing spin-spin correlation functions at relatively short distances, we show that equilibrium near-critical properties can be extracted at short times after quenches into the vicinity of a quantum critical point. The time scales after which equilibrium properties can be extracted are sufficiently short so that the proposed scheme should be viable for quantum simulators of spin models based on ultracold atoms or trapped ions. Our results, analytic as well as numeric, are for one-dimensional spin models, either integrable or nonintegrable, but we expect our conclusions to be valid in higher dimensions as well.

Probing the anomalous dynamical phase in long-range quantum spin chains through Fisher-zero lines.

Using the framework of infinite matrix product states, the existence of an anomalous dynamical phase for the transverse-field Ising chain with sufficiently long-range interactions was first reported in J. C. Halimeh and V. Zauner-Stauber [Phys. Rev. B 96, 134427 (2017)10.1103/PhysRevB.96.134427], where it was shown that anomalous cusps arise in the Loschmidt-echo return rate for sufficiently small quenches within the ferromagnetic phase. In this work we further probe the nature of the anomalous phase through calculating the corresponding Fisher-zero lines in the complex time plane. We find that these Fisher-zero lines exhibit a qualitative difference in their behavior, where, unlike in the case of the regular phase, some of them terminate before intersecting the imaginary axis, indicating the existence of smooth peaks in the return rate preceding the cusps. Additionally, we discuss in detail the infinite matrix product state time-evolution method used to calculate Fisher zeros and the Loschmidt-echo return rate using the matrix product state transfer matrix. Our work sheds further light on the nature of the anomalous phase in the long-range transverse-field Ising chain, while the numerical treatment presented can be applied to more general quantum spin chains.

Photorealistic rendering of unidirectional free-space invisibility cloaks.

Carpet or ground-plane invisibility cloaks hide an object in reflection and inhibit transmission experiments by construction. This concept has significantly reduced the otherwise demanding material requirements and has hence enabled various experimental demonstrations. In contrast, free-space invisibility cloaks should work in both reflection and transmission. The fabrication of omnidirectional three-dimensional free-space cloaks still poses significant challenges. Recently, the idea of the carpet cloak has been carried over to experiments on unidirectional free-space invisibility cloaks that only work perfectly for one particular viewing direction and, depending on the design, also for one linear polarization of light only. Here, by using photorealistic ray tracing, we visualize the performance of four types of such unidirectional cloaks in three dimensions for different viewing directions and different polarizations of light, revealing virtues and limitations of these approaches in an intuitive manner.

Photorealistic ray tracing of free-space invisibility cloaks made of uniaxial dielectrics.

The design rules of transformation optics generally lead to spatially inhomogeneous and anisotropic impedance-matched magneto-dielectric material distributions for, e.g., free-space invisibility cloaks. Recently, simplified anisotropic non-magnetic free-space cloaks made of a locally uniaxial dielectric material (calcite) have been realized experimentally. In a two-dimensional setting and for in-plane polarized light propagating in this plane, the cloaking performance can still be perfect for light rays. However, for general views in three dimensions, various imperfections are expected. In this paper, we study two different purely dielectric uniaxial cylindrical free-space cloaks. For one, the optic axis is along the radial direction, for the other one it is along the azimuthal direction. The azimuthal uniaxial cloak has not been suggested previously to the best of our knowledge. We visualize the cloaking performance of both by calculating photorealistic images rendered by ray tracing. Following and complementing our previous ray-tracing work, we use an equation of motion directly derived from Fermat's principle. The rendered images generally exhibit significant imperfections. This includes the obvious fact that cloaking does not work at all for horizontal or for ordinary linear polarization of light. Moreover, more subtle effects occur such as viewing-angle-dependent aberrations. However, we still find amazingly good cloaking performance for the purely dielectric azimuthal uniaxial cloak.

Time-of-flight imaging of invisibility cloaks.

As invisibility cloaking has recently become experimental reality, it is interesting to explore ways to reveal remaining imperfections. In essence, the idea of most invisibility cloaks is to recover the optical path lengths without an object (to be made invisible) by a suitable arrangement around that object. Optical path length is proportional to the time of flight of a light ray or to the optical phase accumulated by a light wave. Thus, time-of-flight images provide a direct and intuitive tool for probing imperfections. Indeed, recent phase-sensitive experiments on the carpet cloak have already made early steps in this direction. In the macroscopic world, time-of-flight images could be measured directly by light detection and ranging (LIDAR). Here, we show calculated time-of-flight images of the conformal Gaussian carpet cloak, the conformal grating cloak, the cylindrical free-space cloak, and of the invisible sphere. All results are obtained by using a ray-velocity equation of motion derived from Fermat's principle.

Newtonian photorealistic ray tracing of grating cloaks and correlation-function-based cloaking-quality assessment.

Grating cloaks are a variation of dielectric carpet (or ground-plane) cloaks. The latter were introduced by Li and Pendry. In contrast to the numerical work involved in the quasi-conformal carpet cloak, the refractive-index profile of a conformal grating cloak follows a closed and exact analytical form. We have previously mentioned that finite-size conformal grating cloaks may exhibit better cloaking than usual finite-size carpet cloaks. In this paper, we directly visualize their performance using photorealistic ray-tracing simulations. We employ a Newtonian approach that is advantageous compared to conventional ray tracing based on Snell's law. Furthermore, we quantify the achieved cloaking quality by computing the cross-correlations of rendered images. The cross-correlations for the grating cloak are much closer to 100% (i.e., ideal) than those for the Gaussian carpet cloak.

Conformal carpet and grating cloaks.

We introduce a class of conformal versions of the previously introduced quasi-conformal carpet cloak, and show how to construct such conformal cloaks for different cloak shapes. Our method provides exact refractive-index profiles in closed mathematical form for the usual carpet cloak as well as for other shapes. By analyzing their asymptotic behavior, we find that the performance of finite-size cloaks becomes much better for metal shapes with zero average value, e.g., for gratings.

Optical microscopy of 3D carpet cloaks:ray-tracing calculations.

In a recent publication (T. Ergin et al., Science 328, 337 (2010)), three-dimensional broadband dielectric carpet cloaks have been fabricated and experimentally characterized by optical bright-field and dark-field microscopy using unpolarized light from an incandescent lamp. A direct comparison with theory has not been provided so far. In the present work, we treat the carpet cloak as well as the entire optical microscope within the ray-optics approximation and the cloak within the effective-medium approximation. We find good qualitative agreement between experimental results and our calculations.

Photorealistic images of carpet cloaks.

Using home-built dedicated ray-tracing software, we simulate photorealistic images of sceneries in three dimensions including dielectric carpet cloaks--i.e., continuously varying refractive-index distributions that allow for invisibility cloaking of a bump in a metallic carpet. Results for the ideal and for a simplified cloak are shown. The presented material gives a visual and intuitive impression of the performance of different arrangements and might be ideally suited for communicating the concepts of transformation optics to the general public.