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.

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.

Universal chiral Luttinger liquid behavior in a graphene fractional quantum Hall point contact.

One-dimensional conductors are described by Luttinger liquid theory, which predicts a power-law suppression of the single-electron tunneling density of states at low voltages. The scaling exponent is predicted to be quantized when tunneling into a single isolated chiral edge state of the fractional quantum Hall effect. We report conductance measurements across a point contact linking integer and fractional quantum Hall edge states (at fillings 1 and [Formula: see text], respectively). At weak coupling, we observe the predicted universal quadratic scaling with temperature and voltage. At strong coupling, we demonstrate perfect Andreev reflection of fractionalized quasiparticles at the point contact. We use the strong coupling physics to realize a nearly dissipationless direct current voltage step-up transformer, whose gain arises directly from topological fractionalization of electrical charge.

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.

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.

Local spectroscopy of a gate-switchable moiré quantum anomalous Hall insulator.

In recent years, correlated insulating states, unconventional superconductivity, and topologically non-trivial phases have all been observed in several moiré heterostructures. However, understanding of the physical mechanisms behind these phenomena is hampered by the lack of local electronic structure data. Here, we use scanning tunnelling microscopy and spectroscopy to demonstrate how the interplay between correlation, topology, and local atomic structure determines the behaviour of electron-doped twisted monolayer-bilayer graphene. Through gate- and magnetic field-dependent measurements, we observe local spectroscopic signatures indicating a quantum anomalous Hall insulating state with a total Chern number of ±2 at a doping level of three electrons per moiré unit cell. We show that the sign of the Chern number and associated magnetism can be electrostatically switched only over a limited range of twist angle and sample hetero-strain values. This results from a competition between the orbital magnetization of filled bulk bands and chiral edge states, which is sensitive to strain-induced distortions in the moiré superlattice.

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.

Inter-valley coherent order and isospin fluctuation mediated superconductivity in rhombohedral trilayer graphene.

Superconductivity was recently discovered in rhombohedral trilayer graphene (RTG) in the absence of a moiré potential. Superconductivity is observed proximate to a metallic state with reduced isospin symmetry, but it remains unknown whether this is a coincidence or a key ingredient for superconductivity. Using a Hartree-Fock analysis and constraints from experiments, we argue that the symmetry breaking is inter-valley coherent (IVC) in nature. We evaluate IVC fluctuations as a possible pairing glue, and find that they lead to chiral unconventional superconductivity when the fluctuations are strong. We further elucidate how the inter-valley Hund's coupling determines the spin-structure of the IVC ground state and breaks the degeneracy between spin-singlet and triplet superconductivity. Remarkably, if the normal state is spin-unpolarized, we find that a ferromagnetic Hund's coupling favors spin-singlet superconductivity, in agreement with experiments. Instead, if the normal state is spin-polarized, then IVC fluctuations lead to spin-triplet pairing.

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.

Visualizing broken symmetry and topological defects in a quantum Hall ferromagnet.

The interaction between electrons in graphene under high magnetic fields drives the formation of a rich set of quantum Hall ferromagnetic (QHFM) phases with broken spin or valley symmetry. Visualizing atomic-scale electronic wave functions with scanning tunneling spectroscopy (STS), we resolved microscopic signatures of valley ordering in QHFM phases and spectral features of fractional quantum Hall phases of graphene. At charge neutrality, we observed a field-tuned continuous quantum phase transition from a valley-polarized state to an intervalley coherent state, with a Kekulé distortion of its electronic density. Mapping the valley texture extracted from STS measurements of the Kekulé phase, we could visualize valley skyrmion excitations localized near charged defects. Our techniques can be applied to examine valley-ordered phases and their topological excitations in a wide range of materials.

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.

Visualizing delocalized correlated electronic states in twisted double bilayer graphene.

The discovery of interaction-driven insulating and superconducting phases in moiré van der Waals heterostructures has sparked considerable interest in understanding the novel correlated physics of these systems. While a significant number of studies have focused on twisted bilayer graphene, correlated insulating states and a superconductivity-like transition up to 12 K have been reported in recent transport measurements of twisted double bilayer graphene. Here we present a scanning tunneling microscopy and spectroscopy study of gate-tunable twisted double bilayer graphene devices. We observe splitting of the van Hove singularity peak by ~20 meV at half-filling of the conduction flat band, with a corresponding reduction of the local density of states at the Fermi level. By mapping the tunneling differential conductance we show that this correlated system exhibits energetically split states that are spatially delocalized throughout the different regions in the moiré unit cell, inconsistent with order originating solely from onsite Coulomb repulsion within strongly-localized orbitals. We have performed self-consistent Hartree-Fock calculations that suggest exchange-driven spontaneous symmetry breaking in the degenerate conduction flat band is the origin of the observed correlated state. Our results provide new insight into the nature of electron-electron interactions in twisted double bilayer graphene and related moiré systems.

Charged skyrmions and topological origin of superconductivity in magic-angle graphene.

Topological solitons, a class of stable nonlinear excitations, appear in diverse domains as in the Skyrme model of nuclear forces. Here, we argue that similar excitations play an important role in a remarkable material obtained on stacking and twisting two sheets of graphene. Close to a magic twist angle, insulating behavior is observed, which gives way to superconductivity on doping. Here, we propose a unifying description of both observations. A symmetry breaking condensate leads to the ordered insulator, while topological solitons in the condensate-skyrmions-are shown to be charge 2e bosons. Condensation of skyrmions leads to a superconductor, whose physical properties we calculate. More generally, we show how topological textures can mitigate Coulomb repulsion and provide a previously unexplored route to superconductivity. Our mechanism not only clarifies why several other moiré materials do not show superconductivity but also points to unexplored platforms where robust superconductivity is anticipated.

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.

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.

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.

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.

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.

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.