Direct observation of a magnetic-field-induced Wigner crystal.

Wigner predicted that when the Coulomb interactions between electrons become much stronger than their kinetic energy, electrons crystallize into a closely packed lattice1. A variety of two-dimensional systems have shown evidence for Wigner crystals2-11 (WCs). However, a spontaneously formed classical or quantum WC has never been directly visualized. Neither the identification of the WC symmetry nor direct investigation of its melting has been accomplished. Here we use high-resolution scanning tunnelling microscopy measurements to directly image a magnetic-field-induced electron WC in Bernal-stacked bilayer graphene and examine its structural properties as a function of electron density, magnetic field and temperature. At high fields and the lowest temperature, we observe a triangular lattice electron WC in the lowest Landau level. The WC possesses the expected lattice constant and is robust between filling factor ν ≈ 0.13 and ν ≈ 0.38 except near fillings where it competes with fractional quantum Hall states. Increasing the density or temperature results in the melting of the WC into a liquid phase that is isotropic but has a modulated structure characterized by the Bragg wavevector of the WC. At low magnetic fields, the WC unexpectedly transitions into an anisotropic stripe phase, which has been commonly anticipated to form in higher Landau levels. Analysis of individual lattice sites shows signatures that may be related to the quantum zero-point motion of electrons in the WC lattice.

Imaging tunable Luttinger liquid systems in van der Waals heterostructures.

One-dimensional (1D) interacting electrons are often described as a Luttinger liquid1-4 having properties that are intrinsically different from those of Fermi liquids in higher dimensions5,6. In materials systems, 1D electrons exhibit exotic quantum phenomena that can be tuned by both intra- and inter-1D-chain electronic interactions, but their experimental characterization can be challenging. Here we demonstrate that layer-stacking domain walls (DWs) in van der Waals heterostructures form a broadly tunable Luttinger liquid system, including both isolated and coupled arrays. We have imaged the evolution of DW Luttinger liquids under different interaction regimes tuned by electron density using scanning tunnelling microscopy. Single DWs at low carrier density are highly susceptible to Wigner crystallization consistent with a spin-incoherent Luttinger liquid, whereas at intermediate densities dimerized Wigner crystals form because of an enhanced magneto-elastic coupling. Periodic arrays of DWs exhibit an interplay between intra- and inter-chain interactions that gives rise to new quantum phases. At low electron densities, inter-chain interactions are dominant and induce a 2D electron crystal composed of phased-locked 1D Wigner crystal in a staggered configuration. Increased electron density causes intra-chain fluctuation potentials to dominate, leading to an electronic smectic liquid crystal phase in which electrons are ordered with algebraical correlation decay along the chain direction but disordered between chains. Our work shows that layer-stacking DWs in 2D heterostructures provides opportunities to explore Luttinger liquid physics.

Quantum textures of the many-body wavefunctions in magic-angle graphene.

Interactions among electrons create novel many-body quantum phases of matter with wavefunctions that reflect electronic correlation effects, broken symmetries and collective excitations. Many quantum phases have been discovered in magic-angle twisted bilayer graphene (MATBG), including correlated insulating1, unconventional superconducting2-5 and magnetic topological6-9 phases. The lack of microscopic information10,11 of possible broken symmetries has hampered our understanding of these phases12-17. Here we use high-resolution scanning tunnelling microscopy to study the wavefunctions of the correlated phases in MATBG. The squares of the wavefunctions of gapped phases, including those of the correlated insulating, pseudogap and superconducting phases, show distinct broken-symmetry patterns with a √3 × âˆš3 super-periodicity on the graphene atomic lattice that has a complex spatial dependence on the moiré scale. We introduce a symmetry-based analysis using a set of complex-valued local order parameters, which show intricate textures that distinguish the various correlated phases. We compare the observed quantum textures of the correlated insulators at fillings of ±2 electrons per moiré unit cell to those expected for proposed theoretical ground states. In typical MATBG devices, these textures closely match those of the proposed incommensurate Kekulé spiral order15, whereas in ultralow-strain samples, our data have local symmetries like those of a time-reversal symmetric intervalley coherent phase12. Moreover, the superconducting state of MATBG shows strong signatures of intervalley coherence, only distinguishable from those of the insulator with our phase-sensitive measurements.

Continuous symmetry breaking in a two-dimensional Rydberg array.

Spontaneous symmetry breaking underlies much of our classification of phases of matter and their associated transitions1-3. The nature of the underlying symmetry being broken determines many of the qualitative properties of the phase; this is illustrated by the case of discrete versus continuous symmetry breaking. Indeed, in contrast to the discrete case, the breaking of a continuous symmetry leads to the emergence of gapless Goldstone modes controlling, for instance, the thermodynamic stability of the ordered phase4,5. Here, we realize a two-dimensional dipolar XY model that shows a continuous spin-rotational symmetry using a programmable Rydberg quantum simulator. We demonstrate the adiabatic preparation of correlated low-temperature states of both the XY ferromagnet and the XY antiferromagnet. In the ferromagnetic case, we characterize the presence of a long-range XY order, a feature prohibited in the absence of long-range dipolar interaction. Our exploration of the many-body physics of XY interactions complements recent works using the Rydberg-blockade mechanism to realize Ising-type interactions showing discrete spin rotation symmetry6-9.

Electrical Control of Spin and Valley in Spin-Orbit Coupled Graphene Multilayers.

Electrical control of magnetism has been a major technological pursuit of the spintronics community, owing to its far-reaching implications for data storage and transmission. Here, we propose and analyze a new mechanism for electrical switching of isospin, using chiral-stacked graphene multilayers, such as Bernal bilayer graphene or rhombohedral trilayer graphene, encapsulated by transition metal dichalcogenide (TMD) substrates. Leveraging the proximity-induced spin-orbit coupling from the TMD, we demonstrate electrical switching of correlation-induced spin and/or valley polarization, by reversing a perpendicular displacement field or the chemical potential. We substantiate our proposal with both analytical arguments and self-consistent Hartree-Fock numerics. Finally, we illustrate how the relative alignment of the TMDs, together with the top and bottom gate voltages, can be used to selectively switch distinct isospin flavors, putting forward correlated Van der Waals heterostructures as a promising platform for spintronics and valleytronics.

Energy Gap of the Even-Denominator Fractional Quantum Hall State in Bilayer Graphene.

Bernal bilayer graphene hosts even-denominator fractional quantum Hall states thought to be described by a Pfaffian wave function with non-Abelian quasiparticle excitations. Here, we report the quantitative determination of fractional quantum Hall energy gaps in bilayer graphene using both thermally activated transport and by direct measurement of the chemical potential. We find a transport activation gap of 5.1 K at B=12 T for a half filled N=1 Landau level, consistent with density matrix renormalization group calculations for the Pfaffian state. However, the measured thermodynamic gap of 11.6 K is smaller than theoretical expectations for the clean limit by approximately a factor of 2. We analyze the chemical potential data near fractional filling within a simplified model of a Wigner crystal of fractional quasiparticles with long-wavelength disorder, explaining this discrepancy. Our results quantitatively establish bilayer graphene as a robust platform for probing the non-Abelian anyons expected to arise as the elementary excitations of the even-denominator state.

Absolutely Stable Time Crystals at Finite Temperature.

We show that locally interacting, periodically driven (Floquet) Hamiltonian dynamics coupled to a Langevin bath support finite-temperature discrete time crystals (DTCs) with an infinite autocorrelation time. By contrast to both prethermal and many-body localized DTCs, the time crystalline order we uncover is stable to arbitrary perturbations, including those that break the time translation symmetry of the underlying drive. Our approach utilizes a general mapping from probabilistic cellular automata to open classical Floquet systems undergoing continuous-time Langevin dynamics. Applying this mapping to a variant of the Toom cellular automaton, which we dub the "π-Toom time crystal," leads to a 2D Floquet Hamiltonian with a finite-temperature DTC phase transition. We provide numerical evidence for the existence of this transition, and analyze the statistics of the finite temperature fluctuations. Finally, we discuss how general results from the field of probabilistic cellular automata imply the existence of discrete time crystals (with an infinite autocorrelation time) in all dimensions, d≥1.

Detecting Symmetry Breaking in Magic Angle Graphene Using Scanning Tunneling Microscopy.

A growing body of experimental work suggests that magic angle twisted bilayer graphene exhibits a "cascade" of spontaneous symmetry-breaking transitions, sparking interest in the potential relationship between symmetry breaking and superconductivity. However, it has proven difficult to find experimental probes which can unambiguously identify the nature of the symmetry breaking. Here, we show how atomically resolved scanning tunneling microscopy can be used as a fingerprint of symmetry-breaking order. By analyzing the pattern of sublattice polarization and "Kekulé" distortions in small magnetic fields, order parameters for each of the most competitive symmetry-breaking states can be identified. In particular, we show that the "Kramers intervalley coherent state," which theoretical work predicts to be the ground state at even integer fillings, shows a Kekulé distortion which emerges only in a magnetic field.

Preparation of Low Entropy Correlated Many-Body States via Conformal Cooling Quenches.

We propose and analyze a method for preparing low entropy many-body states in isolated quantum optical systems of atoms, ions, and molecules. Our approach is based upon shifting entropy between different regions of a system by spatially modulating the magnitude of the effective Hamiltonian. We conduct two case studies, on a topological spin chain and the spinful fermionic Hubbard model, focusing on the key question: can a "conformal cooling quench" remove sufficient entropy within experimentally accessible timescales? Finite-temperature, time-dependent matrix product state calculations reveal that even moderately sized bath regions can remove enough energy and entropy density to expose coherent low-temperature physics. The protocol is particularly natural in systems with long-range interactions, such as lattice-trapped polar molecules and Rydberg-excited atoms, where the magnitude of the Hamiltonian scales directly with the interparticle spacing. To this end, we propose simple, near-term implementations of conformal cooling quenches in systems of atoms or molecules, where signatures of low-temperature phases may be observed.

Strain-Induced Quantum Phase Transitions in Magic-Angle Graphene.

We investigate the effect of uniaxial heterostrain on the interacting phase diagram of magic-angle twisted bilayer graphene. Using both self-consistent Hartree-Fock and density-matrix renormalization group calculations, we find that small strain values (ε∼0.1%-0.2%) drive a zero-temperature phase transition between the symmetry-broken "Kramers intervalley-coherent" insulator and a nematic semimetal. The critical strain lies within the range of experimentally observed strain values, and we therefore predict that strain is at least partly responsible for the sample-dependent experimental observations.

Phases of the (2+1) Dimensional SO(5) Nonlinear Sigma Model with Topological Term.

We use the half-filled zeroth Landau level in graphene as a regularization scheme to study the physics of the SO(5) nonlinear sigma model subject to a Wess-Zumino-Witten topological term in 2+1 dimensions. As shown by Ippoliti et al. [Phys. Rev. B 98, 235108 (2019)PRBMDO2469-995010.1103/PhysRevB.98.235108], this approach allows for negative sign free auxiliary field quantum Monte Carlo simulations. The model has a single free parameter U_{0} that monitors the stiffness. Within the parameter range accessible to negative sign free simulations, we observe an ordered phase in the large U_{0} or stiff limit. Remarkably, upon reducing U_{0} the magnetization drops substantially, and the correlation length exceeds our biggest system sizes, accommodating 100 flux quanta. The implications of our results for deconfined quantum phase transitions between valence bond solids and antiferromagnets are discussed.

Universal Tripartite Entanglement in One-Dimensional Many-Body Systems.

Motivated by conjectures in holography relating the entanglement of purification and reflected entropy to the entanglement wedge cross section, we introduce two related non-negative measures of tripartite entanglement g and h. We prove structure theorems which show that states with nonzero g or h have nontrivial tripartite entanglement. We then establish that in one dimension these tripartite entanglement measures are universal quantities that depend only on the emergent low-energy theory. For a gapped system, we argue that either g≠0 and h=0 or g=h=0, depending on whether the ground state has long-range order. For a critical system, we develop a numerical algorithm for computing g and h from a lattice model. We compute g and h for various CFTs and show that h depends only on the central charge whereas g depends on the whole operator content.

Experimental Determination of the Energy per Particle in Partially Filled Landau Levels.

We describe an experimental technique to measure the chemical potential µ in atomically thin layered materials with high sensitivity and in the static limit. We apply the technique to a high quality graphene monolayer to map out the evolution of µ with carrier density throughout the N=0 and N=1 Landau levels at high magnetic field. By integrating µ over filling factor ν, we obtain the ground state energy per particle, which can be directly compared to numerical calculations. In the N=0 Landau level, our data show exceptional agreement with numerical calculations over the whole Landau level without adjustable parameters as long as the screening of the Coulomb interaction by the filled Landau levels is accounted for. In the N=1 Landau level, a comparison between experimental and numerical data suggests the importance of valley anisotropic interactions and reveals a possible presence of valley-textured electron solids near odd filling.

Isometric Tensor Network States in Two Dimensions.

Tensor-network states (TNS) are a promising but numerically challenging tool for simulating two-dimensional (2D) quantum many-body problems. We introduce an isometric restriction of the TNS ansatz that allows for highly efficient contraction of the network. We consider two concrete applications using this ansatz. First, we show that a matrix-product state representation of a 2D quantum state can be iteratively transformed into an isometric 2D TNS. Second, we introduce a 2D version of the time-evolving block decimation algorithm for approximating of the ground state of a Hamiltonian as an isometric TNS-which we demonstrate for the 2D transverse field Ising model.

Mechanism for Anomalous Hall Ferromagnetism in Twisted Bilayer Graphene.

Motivated by the recent observation of an anomalous Hall effect in twisted bilayer graphene, we use a lowest Landau level model to understand the origin of the underlying symmetry-broken correlated state. This effective model is rooted in the occurrence of Chern bands which arise due to the coupling between the graphene device and its encapsulating substrate. Our model exhibits a phase transition from a spin-valley polarized insulator to a partial or fully valley unpolarized metal as the bandwidth is increased relative to the interaction strength, consistent with experimental observations. In sharp contrast to standard quantum Hall ferromagnetism, the Chern number structure of the flat bands precludes an instability to an intervalley coherent phase, but allows for an excitonic vortex lattice at large interaction anisotropy.

Lattice Homotopy Constraints on Phases of Quantum Magnets.

The Lieb-Schultz-Mattis (LSM) theorem and its extensions forbid trivial phases from arising in certain quantum magnets. Constraining infrared behavior with the ultraviolet data encoded in the microscopic lattice of spins, these theorems tie the absence of spontaneous symmetry breaking to the emergence of exotic phases like quantum spin liquids. In this work, we take a new topological perspective on these theorems, by arguing they originate from an obstruction to "trivializing" the lattice under smooth, symmetric deformations, which we call the "lattice homotopy problem." We conjecture that all LSM-like theorems for quantum magnets (many previously unknown) can be understood from lattice homotopy, which automatically incorporates the full spatial symmetry group of the lattice, including all its point-group symmetries. One consequence is that any spin-symmetric magnet with a half-integer moment on a site with even-order rotational symmetry must be a spin liquid. To substantiate the claim, we prove the conjecture in two dimensions for some physically relevant settings.

Filling-Enforced Gaplessness in Band Structures of the 230 Space Groups.

Nonsymmorphic symmetries like screws and glides produce electron band touchings, obstructing the formation of a band insulator and leading, instead, to metals or nodal semimetals even when the number of electrons in the unit cell is an even integer. Here, we calculate the electron fillings compatible with being a band insulator for all 230 space groups, for noninteracting electrons with time-reversal symmetry. Our bounds are tight-that is, we can rigorously eliminate band insulators at any forbidden filling and produce explicit models for all allowed fillings-and stronger than those recently established for interacting systems. These results provide simple criteria that should help guide the search for topological semimetals and, also, have implications for both the nature and stability of the resulting nodal Fermi surfaces.

Space Group Symmetry Fractionalization in a Chiral Kagome Heisenberg Antiferromagnet.

The anyonic excitations of a spin liquid can feature fractional quantum numbers under space group symmetries. Detecting these fractional quantum numbers, which are analogs of the fractional charge of Laughlin quasiparticles, may prove easier than the direct observation of anyonic braiding and statistics. Motivated by the recent numerical discovery of spin-liquid phases in the kagome Heisenberg antiferromagnet, we theoretically predict the pattern of space group symmetry fractionalization in the kagome lattice SO(3)-symmetric chiral spin liquid. We provide a method to detect these fractional quantum numbers in finite-size numerics which is simple to implement in the density matrix renormalization group. Applying these developments to the chiral spin liquid phase of a kagome Heisenberg model, we find perfect agreement between our theoretical prediction and numerical observations.

Constraints on topological order in mott insulators.

We point out certain symmetry induced constraints on topological order in Mott insulators (quantum magnets with an odd number of spin 1/2 moments per unit cell). We show, for example, that the double-semion topological order is incompatible with time reversal and translation symmetry in Mott insulators. This sharpens the Hastings-Oshikawa-Lieb-Schultz-Mattis theorem for 2D quantum magnets, which guarantees that a fully symmetric gapped Mott insulator must be topologically ordered, but is silent about which topological order is permitted. Our result applies to the kagome lattice quantum antiferromagnet, where recent numerical calculations of the entanglement entropy indicate a ground state compatible with either toric code or double-semion topological order. Our result rules out the latter possibility.

Topological characterization of fractional quantum Hall ground states from microscopic Hamiltonians.

We show how to numerically calculate several quantities that characterize topological order starting from a microscopic fractional quantum Hall Hamiltonian. To find the set of degenerate ground states, we employ the infinite density matrix renormalization group method based on the matrix-product state representation of fractional quantum Hall states on an infinite cylinder. To study localized quasiparticles of a chosen topological charge, we use pairs of degenerate ground states as boundary conditions for the infinite density matrix renormalization group. We then show that the wave function obtained on the infinite cylinder geometry can be adapted to a torus of arbitrary modular parameter, which allows us to explicitly calculate the non-Abelian Berry connection associated with the modular T transformation. As a result, the quantum dimensions, topological spins, quasiparticle charges, chiral central charge, and Hall viscosity of the phase can be obtained using data contained entirely in the entanglement spectrum of an infinite cylinder.