The behavior of strongly interacting electrons in bands with Berry curvature is a problem of wide interest. In this paper, we study this problem by numerically studying a fluxed Hubbard-type model on square lattice. Using this model, we demonstrate a metallic ferromagnet in electron bands equipped with Berry curvature can develop non-coplanar spin order in which spin polarization axes at different position span finite solid angles. We find spin chirality can emerge in this setting by doping or adding gauge flux on top of a collinear ferromagnet. This result supports the prediction of spin chirality occurring through an emergent spin orbital interaction. Meanwhile, our result shows that, on top a ferromagnetic background, the spin chirality emerges at a finite threshold value of orbital magnetization, resembling the predicted behavior in theory.

Bernal bilayer graphene has recently been discovered to exhibit a wide range of unique ordered phases resulting from interaction-driven effects and encompassing spin and valley magnetism, correlated insulators, correlated metals, and superconductivity. This Letter reports on a novel family of correlated phases characterized by spin and valley ordering, distinct from those reported previously. These phases emerge in electron-doped bilayer graphene where the energy bands are exceptionally flat, manifested through an intriguing nonlinear current-bias behavior that occurs at the onset of the phases and is accompanied by an insulating temperature dependence. These characteristics align with the presence of charge- or spin-density-wave states that open a gap on a portion of the Fermi surface or fully gapped Wigner crystals, resulting in an exceptionally intricate phase diagram.

Spin chirality in electron gases is mandated by symmetry to be nonzero whenever time reversal symmetry (TRS) is broken. In metals, both normal and superconducting, spin chirality is generated by orbital currents present in the ground state as a result of TRS breaking. Nonzero chirality is found to emerge in the absence of any spin-dependent interactions, even when the ground state remains spin-unpolarized. We present a diagrammatic derivation of the chirality effect and illustrate its generality for problems such as valley magnetism and topological superconductivity. Spin chirality results in a chiral three-spin RKKY interaction between localized spins coupled to carriers by s-d interaction, an effect that can be detected by local probes such as spin-sensitive STM. In systems where detecting TRS breaking by macroscopic measurements is challenging, such as topological superconductors, local detection of spin chirality can serve as a reliable diagnostic of broken time reversal symmetry.

This paper aims to clarify the nature of a surprising ordered phase recently reported in biased Bernal bilayer graphene that occurs at the phase boundary between the isospin-polarized and unpolarized phases. Strong nonlinearity of transport at abnormally small currents, with dI/dV vs. I sharply rising and then falling back, is typical for a charge/spin-density-wave state (CDW or SDW) sliding transport. Here, however, it is observed at an isospin-order phase boundary, prompting a question about the CDW/SDW mechanism and its relation to the quantum critical point. We argue that the observed phase diagram cannot be understood within a standard weak-coupling picture. Rather, it points to a mechanism that relies on an effective interaction enhancement at a quantum critical point. We develop a detailed strong-coupling framework accounting for the soft collective modes that explain these observations.

Controlling the bandstructure of Dirac materials is of wide interest in current research but has remained an outstanding challenge for systems such as monolayer graphene. In contrast, Bernal bilayer graphene (BLG) offers a highly flexible platform for tuning the bandstructure, featuring two distinct regimes. In one regime, which is well established and widely used, a tunable bandgap is induced by a large enough transverse displacement field. Another is a gapless metallic band occurring near charge neutrality and at not too strong fields, featuring rich 'fine structure' consisting of four linearly-dispersing Dirac cones with opposite chiralities in each valley and van Hove singularities. Even though BLG was extensively studied experimentally in the last two decades, the evidence of this exotic bandstructure is still elusive, likely due to insufficient energy resolution. Here, rather than probing the bandstructure by direct spectroscopy, we use Landau levels as markers of the energy dispersion and carefully analyze the Landau level spectrum in a regime where the cyclotron orbits of electrons or holes in momentum space are small enough to resolve the distinct mini Dirac cones. We identify the presence of four distinct Dirac cones and map out complex topological transitions induced by electric displacement field. These findings introduce a valuable addition to the toolkit for graphene electronics.

Stoner magnetism in bands endowed with Berry curvature is shown to be influenced by the coupling between spin-chirality density s (partial derivative(x)s x partial derivative (y) s ) and Berry's orbital magnetization. The key effect is that carriers moving in the presence of a spin texture see it as a source of a geometric magnetic field coupled to the carrier's orbital motion through a spin-dependent Aharonov-Bohm effect. This interaction was recently predicted to enable chiral magnons propagating along system boundaries. Here we show that it also favors chiral spin textures such as skyrmions-the topologically protected objects with particle-like properties, stabilized in the ground state. Unlike previously studied systems, here skyrmion textures can arise in the absence of microscopic spin-dependent interactions such as spin-orbit coupling or Zeeman coupling. The threshold for Stoner instability is found to soften, rendering chiral spin-ordered phases accessible under realistic conditions. We present a detailed analysis of the chiral effect for Bernal bilayer graphene and discuss the unique properties of skyrmion textures in graphene multilayers.

Numerous correlated many-body phases, both conventional and exotic, have been reported in magic-angle twisted bilayer graphene (MATBG)1-24. However, the dynamics associated with these correlated states, crucial for understanding the underlying physics, remain unexplored. Here we combine exciton sensing and optical pump-probe spectroscopy to investigate the dynamics of isospin orders in MATBG with WSe2 substrate across the entire flat band, achieving sub-picosecond resolution. We observe remarkably slow isospin dynamics in a broad filling range around ν = 2 and between ν = -3 and -2, with lifetimes of up to 300 ps that decouple from the much faster cooling of electronic temperature (about 10 ps). This non-thermal behaviour demonstrates the presence of abnormally long-lived modes in the isospin degrees of freedom. This observation, not anticipated by theory, implies the existence of long-range propagating collective modes, strong isospin fluctuations and memory effects and is probably associated with an intervalley coherent or incommensurate Kekulé spiral ground state. We further demonstrate non-equilibrium control of the isospin orders previously found around integer fillings. Specifically, through ultrafast manipulation, it can be transiently shifted away from integer fillings. Our study demonstrates a unique probe of collective excitations in MATBG and paves the way for actively controlling non-equilibrium phenomena in moiré systems.

Quantum oscillations originating from the quantization of electron cyclotron orbits provide sensitive diagnostics of electron bands and interactions. We report on nanoscale imaging of the thermodynamic magnetization oscillations caused by the de Haas–van Alphen effect in moiré graphene. Scanning by means of superconducting quantum interference device (SQUID)–on-tip in Bernal bilayer graphene crystal axis-aligned to hexagonal boron nitride reveals large magnetization oscillations with amplitudes reaching 500 Bohr magneton per electron in weak magnetic fields, unexpectedly low frequencies, and high sensitivity to superlattice filling fraction. The oscillations allow us to reconstruct the complex band structure, revealing narrow moiré bands with multiple overlapping Fermi surfaces separated by unusually small momentum gaps. We identified sets of oscillations that violate the textbook Onsager Fermi surface sum rule, signaling formation of broad-band particle-hole superposition states induced by coherent magnetic breakdown.

Spin excitations that soften near the onset of magnetic order have long been known to act as `paramagnon' pairing glue that can drive spin-triplet superconductivity. Recent findings of superconductivity in graphene bilayers and trilayers, occurring in the proximity of different itinerant ordered phases polarized in isospin (spin and valley), have motivated us to conduct a comprehensive investigation of an isospin extension of the paramagnon pairing mechanism in the vicinity of spin/isospin orders. In each case, we identify a soft mode, associated with the order parameter fluctuations, that mediates pairing interaction. We develop an approach that relates the soft mode described through summation of the contributions most strongly divergent at the onset of spin/valley isospin orders. This interaction is not always attractive, but if it is, it gives rise to an enhancement of superconducting T_c in an appropriate pairing channel. In the cases when the pairing interaction is attractive, it leads to the formation of a superconducting state which can be either spin-triplet and valley-singlet or vice versa, depending on the specific isospin order type. These findings demonstrate that the occurrence of superconductivity in the vicinity of an itinerant magnetic phase is a generic phenomenon, closely mirroring experimental observations.

Electron bands in the untwisted bilayer graphene flatten out in a transverse electric field, offering a promising platform for correlated electron physics. We predict that the spin/valley isospin magnetism, resembling that seen in moiré bands, coexists with momentum-polarized phases occurring via a ``flocking transition'' in momentum space. This transition results in the electron momentum distribution being spontaneously displaced relative to the $K$ and $K'$ valley centers. These phases feature unusual observables such as persistent currents in the ground state. Momentum-polarized carriers ``sample'' the Berry curvature of the conduction band, resulting in a unique behavior of the anomalous Hall conductivity and other effects that do not occur in previously studied systems.

Different superconducting pairing mechanisms are markedly distinct in the underlying Cooper pair kinematics. Quantum-critical soft modes drive pairing interactions in which the pair scattering processes are highly collinear and can be classified into two categories: forward scattering and backscattering. Conversely, in conventional phonon mechanisms, Cooper pair scattering is of a generic noncollinear character. In this study, we present a method to discern the kinematic type by observing the evolution of superconductivity while adjusting the Fermi surface geometry. To demonstrate our approach, we utilize the recently reported phase diagrams of untwisted graphene multilayers. Our analysis connects the emergence of superconductivity at "ghost crossings" of Fermi surfaces in distinct valleys to the pair kinematics of a backscattering type. Together with the observed nonmonotonic behavior of superconductivity near its onset (sharp rise followed by a drop), it lends strong support to a particular quantum-critical superconductivity scenario in which pairing is driven by intervalley coherence fluctuations. These findings offer direct insights into the genesis of pairing in these systems, providing compelling evidence for the electron-electron interactions driving superconductivity. More broadly, our work highlights the potential of tuning bands via ghost crossings as a promising means of boosting superconductivity.

The quest for unconventional superconductivity governed by Coulomb repulsion between electrons rather than phonon attraction received new momentum with the advent of moir\'e graphene. Initially, delineating the phonon and Coulomb-repulsion-based pairing mechanisms has proven to be a challenging task, however the situation has changed after recent discovery of superconductivity in non-twisted graphene bilayers and trilayers. Superconductivity occurring at the phase boundaries of spin and valley polarized orders calls for non-phonon scenarios, yet the specific pairing mechanisms remain to be understood. Here we analyze a striking example -- superconductivity in graphene bilayers occurring at the onset of valley-polarized order induced by a magnetic field. We describe an attraction-from-repulsion mechanism for pairing mediated by a quantum-critical mode, which fully explains the observed phenomenology. While it is usually notoriously difficult to infer the pairing mechanism from the observed superconducting phases, this case presents a rare exception, allowing for a fairly unambiguous identification of the origin of the pairing glue. A combination of factors such as the location of superconducting phase at the onset of isospin-polarized phase, a threshold in a magnetic field, above which superconductivity occurs, and its resilience at high magnetic fields paints a clear picture of a triplet superconductivity driven by quantum-critical fluctuations.

We argue that spin- and valley-polarized metallic phases recently observed in graphene bilayers and trilayers support chiral edge modes that allow spin waves to propagate ballistically along system boundaries without backscattering. The chiral edge behavior originates from the interplay between the momentum-space Berry curvature in Dirac bands and the geometric phase of a spin texture in position space. The edge modes are weakly confined to the edge, featuring dispersion that is robust and insensitive to the detailed profile of magnetization at the edge. This unique character of edge modes reduces their overlap with edge disorder and enhances the mode lifetime. The mode propagation direction reverses upon reversing valley polarization, an effect that provides a clear testable signature of geometric interactions in isospin-polarized Dirac bands.

In thermodynamic equilibrium, current in metallic systems is carried by electronic states near the Fermi energy, whereas the filled bands underneath contribute little to conduction. Here, we describe a very different regime in which carrier distribution in graphene and its superlattices is shifted so far from equilibrium that the filled bands start playing an essential role, leading to a critical-current behavior. The criticalities develop upon the velocity of electron flow reaching the Fermi velocity. Key signatures of the out-of-equilibrium state are current-voltage characteristics that resemble those of superconductors, sharp peaks in differential resistance, sign reversal of the Hall effect, and a marked anomaly caused by the Schwinger-like production of hot electron-hole plasma. The observed behavior is expected to be common to all graphene-based superlattices.

We present a theory of superconducting pairing originating from soft critical fluctuations near isospin-polarized states in rhombohedral trilayer graphene. Using a symmetry-based approach, we determine possible isospin order types and derive the effective electron-electron interactions mediated by isospin fluctuations. Superconductitivty arising due to these interactions has symmetry and order parameter structure that depend in a unique way on the "mother" isospin order. This model naturally leads to a superconducting phase adjacent to isospin-ordering phase transition, which mimics the behavior observed in experiment. The symmetry of the paired state predicted for the isospin order type inferred in experiments matches the observations. These findings support a scenario of superconductivity originating from electron-electron interactions.

Dragging of light by moving media was predicted by Fresnel1 and verified by Fizeau's celebrated experiments2 with flowing water. This momentous discovery is among the experimental cornerstones of Einstein's special relativity theory and is well understood3,4 in the context of relativistic kinematics. By contrast, experiments on dragging photons by an electron flow in solids are riddled with inconsistencies and have so far eluded agreement with the theory5-7. Here we report on the electron flow dragging surface plasmon polaritons8,9 (SPPs): hybrid quasiparticles of infrared photons and electrons in graphene. The drag is visualized directly through infrared nano-imaging of propagating plasmonic waves in the presence of a high-density current. The polaritons in graphene shorten their wavelength when propagating against the drifting carriers. Unlike the Fizeau effect for light, the SPP drag by electrical currents defies explanation by simple kinematics and is linked to the nonlinear electrodynamics of Dirac electrons in graphene. The observed plasmonic Fizeau drag enables breaking of time-reversal symmetry and reciprocity10 at infrared frequencies without resorting to magnetic fields11,12 or chiral optical pumping13,14. The Fizeau drag also provides a tool with which to study interactions and nonequilibrium effects in electron liquids.

Achieving Bloch oscillations of free carriers under a direct current, a long-sought-after collective many-body behavior, has been challenging due to stringent constraints on the band properties. We argue that the flat bands in moiré graphene fulfill the basic requirements for observing Bloch oscillations, offering an appealing alternative to the stacked quantum wells used in previous work aiming to access this regime. Bloch-oscillating moiré superlattices emit a comblike spectrum of incommensurate frequencies, a property of interest for converting direct currents into high-frequency currents and developing broadband amplifiers in terahertz domain. The oscillations can be synchronized through coupling to an oscillator mode in a photonic or plasmonic resonator. Phase-coherent collective oscillations in the resonant regime provide a realization of current-pumped terahertz lasing.

We discuss superconducting pairing in a narrow conduction band sandwiched between unoccupied and occupied bands, an arrangement that enables an unconventional pairing mechanism governed by Coulomb repulsion. Pairing interaction originates from repulsion-assisted scattering between far-out pair states in the higher-energy bands and those at the Fermi level. Optimizing the bandstructure design and carrier density in order to bring plasma frequency below the bandgap renders the repulsion unscreened for the processes with a large frequency transfer. This allows the pairing to fully benefit from the pristine Coulomb repulsion strength. The repulsion-induced attraction is particularly strong in two dimensions and is assisted by a low density of carriers and the resulting low plasma frequency values. We assess the possible connection of this mechanism to superconductivity in magic-angle twisted bilayer graphene where the bandstructure features wide dispersive upper and lower minibands. We use a simple model to illustrate the importance of the far-out pairs in these bands and predict testable signatures of this superconductivity mechanism.