Unusual ergodic and chaotic properties of trapped hard rods

We investigate ergodicity, chaos, and thermalization for a one-dimensional classical gas of hard rods confined to an external quadratic or quartic trap, which breaks microscopic integrability. To quantify the strength of chaos in this system, we compute its maximal Lyapunov exponent numerically. The approach to thermal equilibrium is studied by considering the time evolution of particle position and velocity distributions and comparing the late-time profiles with the Gibbs state. Remarkably, we find that quadratically trapped hard rods are highly nonergodic and do not resemble a Gibbs state even at extremely long times, despite compelling evidence of chaos for four or more rods. On the other hand, our numerical results reveal that hard rods in a quartic trap exhibit both chaos and thermalization, and equilibrate to a Gibbs state as expected for a nonintegrable many-body system.

Out-of-time-ordered correlator in the one-dimensional Kuramoto-Sivashinsky and Kardar-Parisi-Zhang equations

The out-of-time-ordered correlator (OTOC) has emerged as an interesting object in both classical and quantum systems for probing the spatial spread and temporal growth of initially local perturbations in spatially extended chaotic systems. Here, we study the (classical) OTOC and its "light cone" in the nonlinear Kuramoto-Sivashinsky (KS) equation, using extensive numerical simulations. We also show that the linearized KS equation exhibits a qualitatively similar OTOC and light cone, which can be understood via a saddle-point analysis of the linearly unstable modes. Given the deep connection between the KS (deterministic) and the Kardar-Parisi-Zhang (KPZ, which is stochastic) equations, we also explore the OTOC in the KPZ equation. While our numerical results in the KS case are expected to hold in the continuum limit, for the KPZ case it is valid in a discretized version of the KPZ equation. More broadly, our work unravels the intrinsic interplay between noise/instability, nonlinearity, and dissipation in partial differential equations (deterministic or stochastic) through the lens of OTOC.

Many-Body Resonances in the Avalanche Instability of Many-Body Localization

Many-body localized (MBL) systems fail to reach thermal equilibrium under their own dynamics, even though they are interacting, nonintegrable, and in an extensively excited state. One instability toward thermalization of MBL systems is the so-called "avalanche," where a locally thermalizing rare region is able to spread thermalization through the full system. The spreading of the avalanche may be modeled and numerically studied in finite one-dimensional MBL systems by weakly coupling an infinite-temperature bath to one end of the system. We find that the avalanche spreads primarily via strong many-body resonances between rare near-resonant eigenstates of the closed system. Thus we find and explore a detailed connection between many-body resonances and avalanches in MBL systems.

Many-Body Quantum Chaos and Emergence of Ginibre Ensemble

We show that non-Hermitian Ginibre random matrix behaviors emerge in spatially extended many-body quantum chaotic systems in the space direction, just as Hermitian random matrix behaviors emerge in chaotic systems in the time direction. Starting with translational invariant models, which can be associated with dual transfer matrices with complex-valued spectra, we show that the linear ramp of the spectral form factor necessitates that the dual spectra have nontrivial correlations, which in fact fall under the universality class of the Ginibre ensemble, demonstrated by computing the level spacing distribution and the dissipative spectral form factor. As a result of this connection, the exact spectral form factor for the Ginibre ensemble can be used to universally describe the spectral form factor for translational invariant many-body quantum chaotic systems in the scaling limit where t and L are large, while the ratio between L and L_{Th}, the many-body Thouless length is fixed. With appropriate variations of Ginibre models, we analytically demonstrate that our claim generalizes to models without translational invariance as well. The emergence of the Ginibre ensemble is a genuine consequence of the strongly interacting and spatially extended nature of the quantum chaotic systems we consider, unlike the traditional emergence of Hermitian random matrix ensembles.

Self-dual quasiperiodic percolation

How does the percolation transition behave in the absence of quenched randomness? To address this question, we study two nonrandom self-dual quasiperiodic models of square-lattice bond percolation. In both models, the critical point has emergent discrete scale invariance, but none of the additional emergent conformal symmetry of critical random percolation. From the discrete sequences of critical clusters, we find fractal dimensions of D_{f}=1.911943(1) and D_{f}=1.707234(40) for the two models, significantly different from D_{f}=91/48=1.89583... of random percolation. The critical exponents ν, determined through a numerical study of cluster sizes and wrapping probabilities on a torus, are also well below the ν=4/3 of random percolation. While these new models do not appear to belong to a universality class, they demonstrate how the removal of randomness can fundamentally change the critical behavior.

Probing site-resolved correlations in a spin system of ultracold molecules

Synthetic quantum systems with interacting constituents play an important role in quantum information processing and in explaining fundamental phenomena in many-body physics. Following impressive advances in cooling and trapping techniques, ensembles of ultracold polar molecules have emerged as a promising platform that combines several advantageous properties^1-11. These include a large set of internal states with long coherence times^12-17 and long-range, anisotropic interactions. These features could enable the exploration of intriguing phases of correlated quantum matter, such as topological superfluids¹⁸, quantum spin liquids¹⁹, fractional Chern insulators²⁰ and quantum magnets^21,22. Probing correlations in these phases is crucial to understanding their properties, necessitating the development of new experimental techniques. Here we use quantum gas microscopy²³ to measure the site-resolved dynamics of quantum correlations of polar ²³Na⁸⁷Rb molecules confined in a two-dimensional optical lattice. By using two rotational states of the molecules, we realize a spin-1/2 system with dipolar interactions between particles, producing a quantum spin-exchange model^21,22,24,25. We study the evolution of correlations during the thermalization process of an out-of-equilibrium spin system for both spatially isotropic and anisotropic interactions. Furthermore, we examine the correlation dynamics of a spin-anisotropic Heisenberg model engineered from the native spin-exchange model by using periodic microwave pulses^26-28. These experiments push the frontier of probing and controlling interacting systems of ultracold molecules, with prospects for exploring new regimes of quantum matter and characterizing entangled states that are useful for quantum computation^29,30 and metrology³¹.

Many-body quantum chaos and space-time translational invariance

We study the consequences of having translational invariance in space and time in many-body quantum chaotic systems. We consider ensembles of random quantum circuits as minimal models of translational invariant many-body quantum chaotic systems. We evaluate the spectral form factor as a sum over many-body Feynman diagrams in the limit of large local Hilbert space dimension q. At sufficiently large t, diagrams corresponding to rigid translations dominate, reproducing the random matrix theory (RMT) behaviour. At finite t, we show that translational invariance introduces additional mechanisms via two novel Feynman diagrams which delay the emergence of RMT. Our analytics suggests the existence of exact scaling forms which describe the approach to RMT behavior in the scaling limit where both t and L are large while the ratio between L and L_Th(t), the many-body Thouless length, is fixed. We numerically demonstrate, with simulations of two distinct circuit models, that the resulting scaling functions are universal in the scaling limit.

Emergence of Fermi's Golden Rule

Fermi's golden rule applies in the limit where an initial quantum state is weakly coupled to a continuum of other final states overlapping its energy. Here we investigate what happens away from this limit, where the set of final states is discrete, with a nonzero mean level spacing; this question arises in a number of recently investigated many-body systems. For different symmetry classes, we analytically and/or numerically calculate the universal crossovers in the average decay of the initial state as the level spacing is varied, with the golden rule emerging in the limit of a continuum. Among the corrections to the exponential decay of the initial state given by Fermi's golden rule is the appearance of the spectral form factor in the longtime regime for small but nonzero level spacing.

Field Theory of Charge Sharpening in Symmetric Monitored Quantum Circuits

Monitored quantum circuits (MRCs) exhibit a measurement-induced phase transition between area-law and volume-law entanglement scaling. MRCs with a conserved charge additionally exhibit two distinct volume-law entangled phases that cannot be characterized by equilibrium notions of symmetry-breaking or topological order, but rather by the nonequilibrium dynamics and steady-state distribution of charge fluctuations. These include a charge-fuzzy phase in which charge information is rapidly scrambled leading to slowly decaying spatial fluctuations of charge in the steady state, and a charge-sharp phase in which measurements collapse quantum fluctuations of charge without destroying the volume-law entanglement of neutral degrees of freedom. By taking a continuous-time, weak-measurement limit, we construct a controlled replica field theory description of these phases and their intervening charge-sharpening transition in one spatial dimension. We find that the charge fuzzy phase is a critical phase with continuously evolving critical exponents that terminates in a modified Kosterlitz-Thouless transition to the short-range correlated charge-sharp phase. We numerically corroborate these scaling predictions also hold for discrete-time projective-measurement circuit models using large-scale matrix-product state simulations, and discuss generalizations to higher dimensions.

Operator Scaling Dimensions and Multifractality at Measurement-Induced Transitions

Repeated local measurements of quantum many-body systems can induce a phase transition in their entanglement structure. These measurement-induced phase transitions (MIPTs) have been studied for various types of dynamics, yet most cases yield quantitatively similar critical exponents, making it unclear how many distinct universality classes are present. Here, we probe the properties of the conformal field theories governing these MIPTs using a numerical transfer-matrix method, which allows us to extract the effective central charge, as well as the first few low-lying scaling dimensions of operators at these critical points for (1+1)-dimensional systems. Our results provide convincing evidence that the generic and Clifford MIPTs for qubits lie in different universality classes and that both are distinct from the percolation transition for qudits in the limit of large on-site Hilbert space dimension. For the generic case, we find strong evidence of multifractal scaling of correlation functions at the critical point, reflected in a continuous spectrum of scaling dimensions.

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

We present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝∑_{i,j=1(i≠j)}^{N}|x_{i}-x_{j}|^{-k}∀k>1 where x_{i} is the position of the ith particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the "light cone." We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1<k<3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

Quantum Ergodicity in the Many-Body Localization Problem

We generalize Page's result on the entanglement entropy of random pure states to the many-body eigenstates of realistic disordered many-body systems subject to long-range interactions. This extension leads to two principal conclusions: first, for increasing disorder the "shells" of constant energy supporting a system's eigenstates fill only a fraction of its full Fock space and are subject to intrinsic correlations absent in synthetic high-dimensional random lattice systems. Second, in all regimes preceding the many-body localization transition individual eigenstates are thermally distributed over these shells. These results, corroborated by comparison to exact diagonalization for an SYK model, are at variance with the concept of "nonergodic extended states" in many-body systems discussed in the recent literature.

Transport, correlations, and chaos in a classical disordered anharmonic chain

We explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size N, disorder strength Δ, and temperature T. The conductivity κ_{N}, obtained for finite length (N), saturates to a value κ_{∞}>0 in the large N limit, for all values of disorder strength Δ and temperature T>0. We show evidence that for any Δ>0 the conductivity goes to zero faster than any power of T in the (T/Δ)→0 limit, and find that the form κ_{∞}∼e^{-B|ln(CΔ/T)|^{3}} fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be κ_{N}<κ_{∞} due to boundary resistance when the bulk conductivity is high (the weak disorder case), or κ_{N}>κ_{∞} due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

Scalable Probes of Measurement-Induced Criticality

We uncover a local order parameter for measurement-induced phase transitions: the average entropy of a single reference qubit initially entangled with the system. Using this order parameter, we identify scalable probes of measurement-induced criticality that are immediately applicable to advanced quantum computing platforms. We test our proposal on a 1+1 dimensional stabilizer circuit model that can be classically simulated in polynomial time. We introduce the concept of a "decoding light cone" to establish the local and efficiently measurable nature of this probe. We also estimate bulk and surface critical exponents for the transition. Developing scalable probes of measurement-induced criticality in more general models may be a useful application of noisy intermediate scale quantum devices, as well as point to more efficient realizations of fault-tolerant quantum computation.

Localization as an Entanglement Phase Transition in Boundary-Driven Anderson Models

The Anderson localization transition is one of the most well studied examples of a zero temperature quantum phase transition. On the other hand, many open questions remain about the phenomenology of disordered systems driven far out of equilibrium. Here we study the localization transition in the prototypical three-dimensional, noninteracting Anderson model when the system is driven at its boundaries to induce a current carrying nonequilibrium steady state. Recently we showed that the diffusive phase of this model exhibits extensive mutual information of its nonequilibrium steady-state density matrix. We show that this extensive scaling persists in the entanglement and at the localization critical point, before crossing over to a short-range (area-law) scaling in the localized phase. We introduce an entanglement witness for fermionic states that we name the mutual coherence, which, for fermionic Gaussian states, is also a lower bound on the mutual information. Through a combination of analytical arguments and numerics, we determine the finite-size scaling of the mutual coherence across the transition. These results further develop the notion of entanglement phase transitions in open systems, with direct implications for driven many-body localized systems, as well as experimental studies of driven-disordered systems.

Bad metallic transport in a cold atom Fermi-Hubbard system

Strong interactions in many-body quantum systems complicate the interpretation of charge transport in such materials. To shed light on this problem, we study transport in a clean quantum system: ultracold lithium-6 in a two-dimensional optical lattice, a testing ground for strong interaction physics in the Fermi-Hubbard model. We determine the diffusion constant by measuring the relaxation of an imposed density modulation and modeling its decay hydrodynamically. The diffusion constant is converted to a resistivity by using the Nernst-Einstein relation. That resistivity exhibits a linear temperature dependence and shows no evidence of saturation, two characteristic signatures of a bad metal. The techniques we developed in this study may be applied to measurements of other transport quantities, including the optical conductivity and thermopower.

Light-Cone Spreading of Perturbations and the Butterfly Effect in a Classical Spin Chain

We find that the effects of a localized perturbation in a chaotic classical many-body system-the classical Heisenberg chain at infinite temperature-spread ballistically with a finite speed even when the local spin dynamics is diffusive. We study two complementary aspects of this butterfly effect: the rapid growth of the perturbation, and its simultaneous ballistic (light-cone) spread, as characterized by the Lyapunov exponents and the butterfly speed, respectively. We connect this to recent studies of the out-of-time-ordered commutators (OTOC), which have been proposed as an indicator of chaos in a quantum system. We provide a straightforward identification of the OTOC with a natural correlator in our system and demonstrate that many of its interesting qualitative features are present in the classical system. Finally, by analyzing the scaling forms, we relate the growth, spread, and propagation of the perturbation with the growth of one-dimensional interfaces described by the Kardar-Parisi-Zhang equation.

Weyl Semimetal to Metal Phase Transitions Driven by Quasiperiodic Potentials

We explore the stability of three-dimensional Weyl and Dirac semimetals subject to quasiperiodic potentials. We present numerical evidence that the semimetal is stable for weak quasiperiodic potentials, despite being unstable for weak random potentials. As the quasiperiodic potential strength increases, the semimetal transitions to a metal, then to an "inverted" semimetal, and then finally to a metal again. The semimetal and metal are distinguished by the density of states at the Weyl point, as well as by level statistics, transport, and the momentum-space structure of eigenstates near the Weyl point. The critical properties of the transitions in quasiperiodic systems differ from those in random systems: we do not find a clear critical scaling regime in energy; instead, at the quasiperiodic transitions, the density of states appears to jump abruptly (and discontinuously to within our resolution).

Obtaining highly excited eigenstates of the localized XX chain via DMRG-X

We benchmark a variant of the recently introduced density matrix renormalization group (DMRG)-X algorithm against exact results for the localized random field XX chain. We find that the eigenstates obtained via DMRG-X exhibit a highly accurate l-bit description for system sizes much bigger than the direct, many-body, exact diagonalization in the spin variables is able to access. We take advantage of the underlying free fermion description of the XX model to accurately test the strengths and limitations of this algorithm for large system sizes. We discuss the theoretical constraints on the performance of the algorithm from the entanglement properties of the eigenstates, and its actual performance at different values of disorder. A small but significant improvement to the algorithm is also presented, which helps significantly with convergence. We find that, at high entanglement, DMRG-X shows a bias towards eigenstates with low entanglement, but can be improved with increased bond dimension. This result suggests that one must be careful when applying the algorithm for interacting many-body localized spin models near a transition.This article is part of the themed issue 'Breakdown of ergodicity in quantum systems: from solids to synthetic matter'.

Thermal inclusions: how one spin can destroy a many-body localized phase

Many-body localized (MBL) systems lie outside the framework of statistical mechanics, as they fail to equilibrate under their own quantum dynamics. Even basic features of MBL systems, such as their stability to thermal inclusions and the nature of the dynamical transition to thermalizing behaviour, remain poorly understood. We study a simple central spin model to address these questions: a two-level system interacting with strength J with N≫1 localized bits subject to random fields. On increasing J, the system transitions from an MBL to a delocalized phase on the vanishing scale Jc(N)∼1/N, up to logarithmic corrections. In the transition region, the single-site eigenstate entanglement entropies exhibit bimodal distributions, so that localized bits are either 'on' (strongly entangled) or 'off' (weakly entangled) in eigenstates. The clusters of 'on' bits vary significantly between eigenstates of the same sample, which provides evidence for a heterogeneous discontinuous transition out of the localized phase in single-site observables. We obtain these results by perturbative mapping to bond percolation on the hypercube at small J and by numerical exact diagonalization of the full many-body system. Our results support the arguments that the MBL phase is unstable in systems with short-range interactions and quenched randomness in dimensions d that are high but finite.This article is part of the themed issue 'Breakdown of ergodicity in quantum systems: from solids to synthetic matter'.

Spin-imbalance in a 2D Fermi-Hubbard system

The interplay of strong interactions and magnetic fields gives rise to unusual forms of superconductivity and magnetism in quantum many-body systems. Here, we present an experimental study of the two-dimensional Fermi-Hubbard model—a paradigm for strongly correlated fermions on a lattice—in the presence of a Zeeman field and varying doping. Using site-resolved measurements, we revealed anisotropic antiferromagnetic correlations, a precursor to long-range canted order. We observed nonmonotonic behavior of the local polarization with doping for strong interactions, which we attribute to the evolution from an antiferromagnetic insulator to a metallic phase. Our results pave the way to experimentally mapping the low-temperature phase diagram of the Fermi-Hubbard model as a function of both doping and spin polarization, for which many open questions remain.

Two Universality Classes for the Many-Body Localization Transition

We provide a systematic comparison of the many-body localization (MBL) transition in spin chains with nonrandom quasiperiodic versus random fields. We find evidence suggesting that these belong to two separate universality classes: the first dominated by "intrinsic" intrasample randomness, and the second dominated by external intersample quenched randomness. We show that the effects of intersample quenched randomness are strongly growing, but not yet dominant, at the system sizes probed by exact-diagonalization studies on random models. Thus, the observed finite-size critical scaling collapses in such studies appear to be in a preasymptotic regime near the nonrandom universality class, but showing signs of the initial crossover towards the external-randomness-dominated universality class. Our results provide an explanation for why exact-diagonalization studies on random models see an apparent scaling near the transition while also obtaining finite-size scaling exponents that strongly violate Harris-Chayes bounds that apply to disorder-driven transitions. We also show that the MBL phase is more stable for the quasiperiodic model as compared to the random one, and the transition in the quasiperiodic model suffers less from certain finite-size effects.

Thermal Equilibrium of a Macroscopic Quantum System in a Pure State

We consider the notion of thermal equilibrium for an individual closed macroscopic quantum system in a pure state, i.e., described by a wave function. The macroscopic properties in thermal equilibrium of such a system, determined by its wave function, must be the same as those obtained from thermodynamics, e.g., spatial uniformity of temperature and chemical potential. When this is true we say that the system is in macroscopic thermal equilibrium (MATE). Such a system may, however, not be in microscopic thermal equilibrium (MITE). The latter requires that the reduced density matrices of small subsystems be close to those obtained from the microcanonical, equivalently the canonical, ensemble for the whole system. The distinction between MITE and MATE is particularly relevant for systems with many-body localization for which the energy eigenfuctions fail to be in MITE while necessarily most of them, but not all, are in MATE. We note, however, that for generic macroscopic systems, including those with MBL, most wave functions in an energy shell are in both MATE and MITE. For a classical macroscopic system, MATE holds for most phase points on the energy surface, but MITE fails to hold for any phase point.

Slowest local operators in quantum spin chains

We numerically construct slowly relaxing local operators in a nonintegrable spin-1/2 chain. Restricting the support of the operator to M consecutive spins along the chain, we exhaustively search for the operator that minimizes the Frobenius norm of the commutator with the Hamiltonian. We first show that the Frobenius norm bounds the time scale of relaxation of the operator at high temperatures. We find operators with significantly slower relaxation than the slowest simple "hydrodynamic" mode due to energy diffusion. Then we examine some properties of the nontrivial slow operators. Using both exhaustive search and tensor network techniques, we find similar slowly relaxing operators for a Floquet spin chain; this system is hydrodynamically "trivial," with no conservation laws restricting their dynamics. We argue that such slow relaxation may be a generic feature following from locality and unitarity.

Thermalization of entanglement

We explore the dynamics of the entanglement entropy near equilibrium in highly entangled pure states of two quantum-chaotic spin chains undergoing unitary time evolution. We examine the relaxation to equilibrium from initial states with either less or more entanglement entropy than the equilibrium value, as well as the dynamics of the spontaneous fluctuations of the entanglement that occur in equilibrium. For the spin chain with a time-independent Hamiltonian and thus an extensive conserved energy, we find slow relaxation of the entanglement entropy near equilibration. Such slow relaxation is absent in a Floquet spin chain with a Hamiltonian that is periodic in time and thus has no local conservation law. Therefore, we argue that slow diffusive energy transport is responsible for the slow relaxation of the entanglement entropy in the Hamiltonian system.

Testing whether all eigenstates obey the eigenstate thermalization hypothesis

We ask whether the eigenstate thermalization hypothesis (ETH) is valid in a strong sense: in the limit of an infinite system, every eigenstate is thermal. We examine expectation values of few-body operators in highly excited many-body eigenstates and search for "outliers," the eigenstates that deviate the most from ETH. We use exact diagonalization of two one-dimensional nonintegrable models: a quantum Ising chain with transverse and longitudinal fields, and hard-core bosons at half-filling with nearest- and next-nearest-neighbor hopping and interaction. We show that even the most extreme outliers appear to obey ETH as the system size increases and thus provide numerical evidences that support ETH in this strong sense. Finally, periodically driving the Ising Hamiltonian, we show that the eigenstates of the corresponding Floquet operator obey ETH even more closely. We attribute this better thermalization to removing the constraint of conservation of the total energy.

Three species of Schrödinger cat states in an infinite-range spin model

We explore a transverse-field Ising model that exhibits both spontaneous symmetry breaking and eigenstate thermalization. Within its ferromagnetic phase, the exact eigenstates of the Hamiltonian of any large but finite-sized system are all Schrödinger cat states: superpositions of states with "up" and "down" spontaneous magnetization. This model exhibits two dynamical phase transitions within its ferromagnetic phase: In the lowest-temperature phase the magnetization can macroscopically oscillate between up and down. The relaxation of the magnetization is always overdamped in the remainder of the ferromagnetic phase, which is divided into phases where the system thermally activates itself over the barrier between the up and down states, and where it quantum tunnels.

Unexpected z-direction Ising antiferromagnetic order in a frustrated spin-1/2 J1-J2 XY model on the honeycomb lattice

Using the density matrix renormalization group on wide cylinders, we study the phase diagram of the spin-1/2 XY model on the honeycomb lattice, with first-neighbor (J1=1) and frustrating second-neighbor (J2>0) interactions. For the intermediate frustration regime 0.22≲J2≲0.36, we find a surprising antiferromagnetic Ising phase, with ordered moments pointing along the z axis, despite the absence of any S(z)S(z) interactions in the Hamiltonian. Surrounding this phase as a function of J2 are antiferromagnetic phases with the moments pointing in the x-y plane for small J2 and a close competition between an x-y plane magnetic collinear phase and a dimer phase for large values of J2. We do not find any spin-liquid phases in this model.

Ballistic spreading of entanglement in a diffusive nonintegrable system

We study the time evolution of the entanglement entropy of a one-dimensional nonintegrable spin chain, starting from random nonentangled initial pure states. We use exact diagonalization of a nonintegrable quantum Ising chain with transverse and longitudinal fields to obtain the exact quantum dynamics. We show that the entanglement entropy increases linearly with time before finite-size saturation begins, demonstrating a ballistic spreading of the entanglement, while the energy transport in the same system is diffusive. Thus, we explicitly demonstrate that the spreading of entanglement is much faster than the energy diffusion in this nonintegrable system.

Weak plaquette valence bond order in the S = 1/2 honeycomb J1 - J2 Heisenberg model

Using the density matrix renormalization group, we investigate the S = 1/2 Heisenberg model on the honeycomb lattice with first (J(1)) and second (J(2)) neighbor interactions. We are able to study long open cylinders with widths up to 12 lattice spacings. For J(2)/J(1) near 0.3, we find an apparently paramagnetic phase, bordered by an antiferromagnetic phase for J(2) ≲ 0.26 and by a valence bond crystal for J(2) ≳ 0.36. The longest correlation length that we find in this intermediate phase is for plaquette valence bond order. This correlation length grows strongly with cylinder circumference, indicating either quantum criticality or weak plaquette valence bond order.

Nonequilibrium dynamic critical scaling of the quantum Ising chain

We solve for the time-dependent finite-size scaling functions of the one-dimensional transverse-field Ising chain during a linear-in-time ramp of the field through the quantum critical point. We then simulate Mott-insulating bosons in a tilted potential, an experimentally studied system in the same equilibrium universality class, and demonstrate that universality holds for the dynamics as well. We find qualitatively athermal features of the scaling functions, such as negative spin correlations, and we show that they should be robustly observable within present cold atom experiments.

Zero- and low-temperature behavior of the two-dimensional ±J Ising spin glass

Scaling arguments and precise simulations are used to study the square lattice ±J Ising spin glass, a prototypical model for glassy systems. Droplet theory explains, and our numerical results show, entropically stabilized long-range spin-glass order at zero temperature, which resembles the energetic stabilization of long-range order in higher-dimensional models at finite temperature. At low temperature, a temperature-dependent crossover length scale is used to predict the power-law dependence on temperature of the heat capacity and clarify the importance of disorder distributions.

Spin-liquid ground state of the S = 1/2 kagome Heisenberg antiferromagnet

We use the density matrix renormalization group to perform accurate calculations of the ground state of the nearest-neighbor quantum spin S = 1/2 Heisenberg antiferromagnet on the kagome lattice. We study this model on numerous long cylinders with circumferences up to 12 lattice spacings. Through a combination of very-low-energy and small finite-size effects, our results provide strong evidence that, for the infinite two-dimensional system, the ground state of this model is a fully gapped spin liquid.

Trimers, molecules, and polarons in mass-imbalanced atomic Fermi gases

We consider the ground state of a single "spin-down" impurity atom interacting attractively with a "spin-up" atomic Fermi gas. By constructing variational wave functions for polarons, molecules, and trimers, we perform a detailed study of the transitions between these dressed bound states as a function of mass ratio r=m↑/m↓ and interaction strength. Crucially, we find that the presence of a Fermi sea enhances the stability of the p-wave trimer, which can be viewed as a Fulde-Ferrell-Larkin-Ovchinnikov molecule that has bound an additional majority atom. For sufficiently large r, we find that the transitions lie outside the region of phase separation of the imbalanced Fermi gas and should thus be observable in experiment, unlike the well-studied equal-mass case.

Collective states of interacting Fibonacci anyons

We show that chains of interacting Fibonacci anyons can support a wide variety of collective ground states ranging from extended critical, gapless phases to gapped phases with ground-state degeneracy and quasiparticle excitations. In particular, we generalize the Majumdar-Ghosh Hamiltonian to anyonic degrees of freedom by extending recently studied pairwise anyonic interactions to three-anyon exchanges. The energetic competition between two- and three-anyon interactions leads to a rich phase diagram that harbors multiple critical and gapped phases. For the critical phases and their higher symmetry end points we numerically establish descriptions in terms of two-dimensional conformal field theories. A topological symmetry protects the critical phases and determines the nature of gapped phases.

Quasi-one-dimensional polarized fermi superfluids

We calculate the zero-temperature (T=0) phase diagram of a polarized two-component Fermi gas in an array of weakly coupled parallel one-dimensional (1D) "tubes" produced by a two-dimensional optical lattice. Increasing the lattice strength drives a crossover from three-dimensional (3D) to 1D behavior, stabilizing the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) modulated superfluid phase. We argue that the most promising regime for observing the FFLO phase is in the quasi-1D regime, where the atomic motion is largely 1D but there is weak tunneling in the other directions that stabilizes long-range order. In the FFLO phase, we describe a phase transition where the quasiparticle spectrum changes from gapless near the 3D regime to gapped in quasi-1D.

Thermodynamics of an Incommensurate Quantum Crystal

We present a theory of the thermodynamics of an incommensurate quantum solid. The ground state of the solid is assumed to be an incommensurate crystal, with quantum zero-point vacancies and interstitials and thus a non-integer number of atoms per unit cell. We show that at low temperature T, the variation of the net vacancy concentration should be as T4 and that the first correction to the specific heat due to this varies as T7; these are quite consistent with experiments on solid helium-4. We also make some observations about the recent experimental reports of "supersolidity" in solid helium-4 that motivate a renewed interest in quantum crystals.

Ordering mechanisms in two-dimensional sphere-forming block copolymers

We study the coarsening dynamics of two-dimensional hexagonal patterns formed by single microdomain layers of block copolymers, using numerical simulations. Our study is focused on the temporal evolution of the orientational correlation length, the interactions between topological defects, and the mechanisms of coarsening. We find no free disclinations in the system; rather, they are located on large-angle grain boundaries, commonly where such boundaries bifurcate. The correlation lengths determined from the scattering function, from the density of dislocations, and from the density of disclinations exhibit similar behavior and grow with time according to a power law. The orientational correlation length also grows following a power law, but with a higher exponent than the other correlation lengths. The orientational correlation length grows via annihilation of dislocations, through preferential annihilation of small-angle grain boundaries due to poor screening of the strain field around dislocations located on small-angle grain boundaries. Consequently, the patterns are characterized by large-angle grain boundaries. The most commonly observed mechanism of coarsening is the collapse of smaller grains residing on the boundary of two larger grains delimited by large-angle grain boundaries. Simulations agree remarkably well with experimental results recently obtained.

Optimizing the ensemble for equilibration in broad-histogram Monte Carlo simulations

We present an adaptive algorithm which optimizes the statistical-mechanical ensemble in a generalized broad-histogram Monte Carlo simulation to maximize the system's rate of round trips in total energy. The scaling of the mean round-trip time from the ground state to the maximum entropy state for this local-update method is found to be O ( [N ln N](2) ) for both the ferromagnetic and the fully frustrated two-dimensional Ising model with N spins. Our algorithm thereby substantially outperforms flat-histogram methods such as the Wang-Landau algorithm.

Counterflow measurements in strongly correlated GaAs hole bilayers: evidence for electron-hole pairing

We study interacting GaAs bilayer hole systems, with very small interlayer tunneling, in a counterflow geometry where equal currents are passed in opposite directions in the two, independently contacted layers. At low temperatures, both the longitudinal and Hall counterflow resistances tend to vanish in the quantum Hall state at total bilayer filling nu=1, demonstrating the pairing of oppositely charged carriers in opposite layers. The counterflow Hall resistance decreases much more strongly than the longitudinal resistances as the temperature is reduced.

Coulomb and liquid dimer models in three dimensions

We study classical hard-core dimer models on three-dimensional lattices using analytical approaches and Monte Carlo simulations. On the bipartite cubic lattice, a local gauge field generalization of the height representation used on the square lattice predicts that the dimers are in a critical Coulomb phase with algebraic, dipolar correlations, in excellent agreement with our large-scale Monte Carlo simulations. The nonbipartite fcc and Fisher lattices lack such a representation, and we find that these models have both confined and exponentially deconfined but no critical phases. We conjecture that extended critical phases are realized only on bipartite lattices, even in higher dimensions.

Gaussian superconducting fluctuations, thermal transport, and the nernst effect

We calculate the contribution of superconducting fluctuations to thermal transport in the normal state, at low magnetic fields. We do so in the Gaussian approximation to their critical dynamics which is also the Aslamazov-Larkin approximation in the microscopics. Our results for the thermal conductivity tensor and the transverse thermoelectric response are new. The latter compare favorably with the data of Ong and collaborators on the Nernst effect in the cuprates.

Permutation-symmetric multicritical points in random antiferromagnetic spin chains

We present a general theory of a class of multicritical points in the phase diagrams of random antiferromagnetic spin chains. We show that low-energy properties of these points are almost completely determined by a permutation symmetry of the effective theory not shared by the microscopic Hamiltonian. One case provides an analytic theory of the quantum critical point in the random spin-3/2 chain, studied in a recent work by Refael, Kehrein, and Fisher.

Multiple field-induced phase transitions in the geometrically frustrated dipolar magnet: Gd(2) Ti(2)O(7)

Field-driven phase transitions generally arise from competition between Zeeman energy and exchange or crystal-field anisotropy. Here we present the phase diagram of a frustrated pyrochlore magnet Gd(2)Ti(2)O(7), where crystal-field splitting is small compared to the dipolar energy. We find good agreement between zero-temperature critical fields and those obtained from a mean-field model. Here, dipolar interactions couple real space and spin space, so the transitions in Gd(2)Ti(2)O(7) arise from field-induced "cooperative anisotropy," reflecting the broken spatial symmetries of the pyrochlore lattice.

Dynamics of pattern coarsening in a two-dimensional smectic system

We have followed the coarsening dynamics of a single layer of cylindrical block copolymer microdomains in a thin film. This system has the symmetry of a two-dimensional smectic. The orientational correlation length of the microdomains was measured by scanning electron microscopy and found to grow with the average spacing between +/-1/2 disclinations, following a power law xi2(t) approximately t(1/4). By tracking disclinations during annealing with time-lapse atomic force microscopy, we observe dominant mechanisms of disclination annihilation involving tripoles and quadrupoles (three and four disclinations, respectively). We describe how annihilation events involving multiple disclinations result in similarly reduced kinetic exponents as observed here. These results map onto a wide variety of physical systems that exhibit similarly striped patterns.

Mechanisms of ordering in striped patterns

We have studied the ordering dynamics of the striped patterns of a single layer of cylindrical block copolymer microdomains in a thin film. By tracking disclinations during annealing with time-lapse atomic force microscopy, we observe a dominant mechanism of disclination annihilation involving three or four disclinations (quadrupoles). Pairwise disclination annihilation events are suppressed as a result of the topological constraints in this system. The kinetic scaling laws with exponents observed here are consistent with topologically allowed annihilation events involving multiple disclinations. The results provide insight into two-dimensional pattern formation and may lead to the successful application of block copolymer lithography.