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

Unconventional Discontinuous Transitions in Isospin Systems

We show that two-dimensional fermions with dispersion k^{2} or k^{4} undergo a first-order Stoner transition to a fully spin-polarized state despite the fact that the spin susceptibility diverges at the critical point. We extend our analysis to systems with dispersion k^{2α} and spin and valley isospin and show that there is a cascade of instabilities into fractional-metal states with some electron bands fully depleted; narrow intermediate ranges of partially depleted bands exist for α<1 or α>2. The susceptibility becomes large near each transition. We discuss applications to biased bi- and trilayer graphene and moiré systems.

Nematic Metal in a Multivalley Electron Gas: Variational Monte Carlo Analysis and Application to AlAs

The two-dimensional electron gas is of fundamental importance in quantum many-body physics. We study a minimal extension of this model with C_{4} (as opposed to full rotational) symmetry and an electronic dispersion with two valleys with anisotropic effective masses. Electrons in our model interact via Coulomb repulsion, screened by distant metallic gates. Using variational Monte Carlo simulations, we find a broad intermediate range of densities with a metallic valley-polarized, spin-unpolarized ground state. Our results are of direct relevance to the recently discovered "nematic" state in AlAs quantum wells. For the effective mass anisotropy relevant to this system, m_{x}/m_{y}≈5.2, we obtain a transition from an anisotropic metal to a valley-polarized metal at r_{s}≈12 (where r_{s} is the dimensionless Wigner-Seitz radius). At still lower densities, we find a (possibly metastable) valley and spin-polarized state with a reduced electronic anisotropy.

Emergent Correlated Phases in Rhombohedral Trilayer Graphene Induced by Proximity Spin-Orbit and Exchange Coupling

The impact of proximity-induced spin-orbit and exchange coupling on the correlated phase diagram of rhombohedral trilayer graphene (RTG) is investigated theoretically. By employing ab initio-fitted effective models of RTG encapsulated by transition metal dichalcogenides (spin-orbit proximity effect) and ferromagnetic Cr_{2}Ge_{2}Te_{6} (exchange proximity effect), we incorporate the Coulomb interactions within the random-phase approximation to explore potential correlated phases at different displacement fields and doping. We find a rich spectrum of spin-valley resolved Stoner and intervalley coherence instabilities induced by the spin-orbit proximity effects, such as the emergence of a spin-valley-coherent phase due to the presence of valley-Zeeman coupling. Similarly, proximity exchange removes the phase degeneracies by biasing the spin direction, enabling a magnetocorrelation effect-strong sensitivity of the correlated phases to the relative magnetization orientations (parallel or antiparallel) of the encapsulating ferromagnetic layers.

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.

Stoner Ferromagnetism in a Momentum-Confined Interacting 2D Electron Gas

In this work we investigate the ground state of a momentum-confined interacting 2D electron gas, a momentum-space analog of an infinite quantum well. The study is performed by combining analytical results with a numerical exact diagonalization procedure. We find a ferromagnetic ground state near a particular electron density and for a range of effective electron (or hole) masses. We argue that this observation may be relevant to the generalized Stoner ferromagnetism recently observed in multilayer graphene systems. The collective magnon excitations exhibit a linear dispersion, which originates from a diverging spin stiffness.

Gate-Defined Topological Josephson Junctions in Bernal Bilayer Graphene

Recent experiments on Bernal bilayer graphene (BLG) deposited on monolayer WSe_{2} revealed robust, ultraclean superconductivity coexisting with sizable induced spin-orbit coupling. Here, we propose BLG/WSe_{2} as a platform to engineer gate-defined planar topological Josephson junctions, where the normal and superconducting regions descend from a common material. More precisely, we show that if superconductivity in BLG/WSe_{2} is gapped and emerges from a parent state with intervalley coherence, then Majorana zero-energy modes can form in the barrier region upon applying weak in-plane magnetic fields. Our results spotlight a potential pathway for "internally engineered" topological superconductivity that minimizes detrimental disorder and orbital-magnetic-field effects.

Signatures of Cooper pair dynamics and quantum-critical superconductivity in tunable carrier bands

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.

Shear Strain-Induced Two-Dimensional Slip Avalanches in Rhombohedral MoS2

Slip avalanches are ubiquitous phenomena occurring in three-dimensional materials under shear strain, and their study contributes immensely to our understanding of plastic deformation, fragmentation, and earthquakes. So far, little is known about the role of shear strain in two-dimensional (2D) materials. Here we show some evidence of 2D slip avalanches in exfoliated rhombohedral MoS2, triggered by shear strain near the threshold level. Utilizing interfacial polarization in 3R-MoS2, we directly probe the stacking order in multilayer flakes and discover a wide variety of polarization domains with sizes following a power-law distribution. These findings suggest that slip avalanches can occur during the exfoliation of 2D materials, and the stacking orders can be changed via shear strain. Our observation has far-reaching implications for the development of new materials and technologies, where precise control over the atomic structure of these materials is essential for optimizing their properties as well as for our understanding of fundamental physical phenomena.

Ising superconductivity induced from spin-selective valley symmetry breaking in twisted trilayer graphene

We show that the e-e interaction induces a strong breakdown of valley symmetry for each spin channel in twisted trilayer graphene, leading to a ground state where the two spin projections have opposite sign of the valley symmetry breaking order parameter. This leads to a spin-valley locking in which the electrons of a Cooper pair are forced to live on different Fermi lines attached to opposite valleys. Furthermore, we find an effective intrinsic spin-orbit coupling explaining the protection of the superconductivity against in-plane magnetic fields. The effect of spin-selective valley symmetry breaking is validated as it reproduces the experimental observation of the reset of the Hall density at 2-hole doping. It also implies a breakdown of the symmetry of the bands from C₆ to C₃, with an enhancement of the anisotropy of the Fermi lines which is at the origin of a Kohn-Luttinger (pairing) instability. The isotropy of the bands is gradually recovered, however, when the Fermi level approaches the bottom of the second valence band, explaining why the superconductivity fades away in the doping range beyond 3 holes per moiré unit cell in twisted trilayer graphene.

Stabilizing Fluctuating Spin-Triplet Superconductivity in Graphene via Induced Spin-Orbit Coupling

A recent experiment showed that a proximity-induced Ising spin-orbit coupling enhances the spin-triplet superconductivity in Bernal bilayer graphene. Here, we show that, due to the nearly perfect spin rotation symmetry of graphene, the fluctuations of the spin orientation of the triplet order parameter suppress the superconducting transition to nearly zero temperature. Our analysis shows that both an Ising spin-orbit coupling and an in-plane magnetic field can eliminate these low-lying fluctuations and can greatly enhance the transition temperature, consistent with the recent experiment. Our model also suggests the possible existence of a phase at small anisotropy and magnetic field which exhibits quasilong-range ordered spin-singlet charge 4e superconductivity, even while the triplet 2e superconducting order only exhibits short-ranged correlations. Finally, we discuss relevant experimental signatures.

Functional Renormalization Group Study of Superconductivity in Rhombohedral Trilayer Graphene

We employ a functional renormalization group approach to ascertain the pairing mechanism and symmetry of the superconducting phase observed in rhombohedral trilayer graphene. Superconductivity in this system occurs in a regime of carrier density and displacement field with a weakly distorted annular Fermi sea. We find that repulsive Coulomb interactions can induce electron pairing on the Fermi surface by taking advantage of momentum-space structure associated with the finite width of the Fermi sea annulus. The degeneracy between spin-singlet and spin-triplet pairing is lifted by valley-exchange interactions that strengthen under the RG flow and develop nontrivial momentum-space structure. We find that the leading pairing instability is d-wave-like and spin singlet, and that the theoretical phase diagram versus carrier density and displacement field agrees qualitatively with experiment.

Optical Conductivity as a Probe of the Interaction-Driven Metal in Rhombohedral Trilayer Graphene

Study of the strongly correlated states in van der Waals heterostructures is one of the central topics in modern condensed matter physics. Among these, the rhombohedral trilayer graphene (RTG) occupies a prominent place since it hosts a variety of interaction-driven phases, with the metallic ones yielding exotic superconducting orders upon doping. Motivated by these experimental findings, we show within the framework of the low-energy Dirac theory that the optical conductivity can distinguish different candidates for a paramagnetic metallic ground state in this system. In particular, this observable shows a single peak in the fully gapped valence-bond state. On the other hand, the bond-current state features two pronounced peaks in the optical conductivity as the probing frequency increases. Finally, the rotational symmetry breaking charge-density wave exhibits a minimal conductivity with the value independent of the amplitude of the order parameter, which corresponds precisely to the splitting of the two cubic nodal points at the two valleys into two triplets of the band touching points featuring linearly dispersing quasiparticles. These features represent the smoking gun signatures of different candidate order parameters for the paramagnetic metallic ground state, which should motivate further experimental studies of the RTG.