Unlike conventional charge-2e superconductors, a charge-4e superconductor exhibits long-range coherence of electron quartets rather than Cooper pairs. Clear zero-temperature realizations of charge-4e superconductivity remain rare. Here, we investigate the zero-temperature phase diagram of the attractive SU(4) Hubbard model with numerically exact, large-scale quantum Monte Carlo (QMC) simulations overcoming major technical hurdles. We identify both charge-2e and charge-4e superconducting phases. Upon increasing interaction, charge-2e correlations are suppressed and eventually vanish, while the charge-4e correlations remain robust and converge with system size, signaling the onset of a quartet-condensed phase. Interestingly, across the charge-2e--charge-4e transition, single electrons remain gapped, while charge-2e correlations exhibit a scaling behavior inconsistent with a conventional Landau description. These features are naturally captured by a fractionalized framework in which the physical charge-2e order parameter is a composite field coupled to an emergent non-Abelian gauge structure. We formulate an Sp(4) gauge-Higgs theory that realizes deconfined quantum pseudocriticality between the Higgs (charge-2e) phase and the confined (charge-4e) phase. The Sp(4) gauge-Higgs theory yields pseudocriticality through a fixed-point collision, and its one-loop collision-point exponents quantitatively track the QMC results. Our results establish charge-4e superconductivity as a bona fide zero-temperature phase, provide a simple model for future studies in a numerically exact framework, and reveal an unconventional route to superconducting criticality.
Unlike conventional charge-2e superconductors, a charge-4e superconductor exhibits long-range coherence of electron quartets rather than Cooper pairs. Clear zero-temperature realizations of charge-4e superconductivity remain rare. Here, we investigate the zero-temperature phase diagram of the attractive SU(4) Hubbard model with numerically exact, large-scale quantum Monte Carlo (QMC) simulations overcoming major technical hurdles. We identify both charge-2e and charge-4e superconducting phases. Upon increasing interaction, charge-2e correlations are suppressed and eventually vanish, while the charge-4e correlations remain robust and converge with system size, signaling the onset of a quartet-condensed phase. Interestingly, across the charge-2e--charge-4e transition, single electrons remain gapped, while charge-2e correlations exhibit a scaling behavior inconsistent with a conventional Landau description. These features are naturally captured by a fractionalized framework in which the physical charge-2e order parameter is a composite field coupled to an emergent non-Abelian gauge structure. We formulate an Sp(4) gauge-Higgs theory that realizes deconfined quantum pseudocriticality between the Higgs (charge-2e) phase and the confined (charge-4e) phase. The Sp(4) gauge-Higgs theory yields pseudocriticality through a fixed-point collision, and its one-loop collision-point exponents quantitatively track the QMC results. Our results establish charge-4e superconductivity as a bona fide zero-temperature phase, provide a simple model for future studies in a numerically exact framework, and reveal an unconventional route to superconducting criticality.
The condensation of electron quartets, known as charge-4e superconductivity (SC), represents a novel quantum state of matter beyond the standard paradigm of Cooper pairing. However, concrete microscopic models realizing this phase in two dimensions remain a central challenge. Here, we introduce a non-engineered and sign-problem-free model, unambiguously demonstrating the emergence of a robust and high-temperature charge-4e SC phase using unbiased quantum Monte Carlo simulations. At zero temperature, the phase diagram reveals that charge-4e SC is the primary ground state in the strong-coupling regime. At finite temperature in the absence of charge-2e SC, we identify charge-4e SC through a Berezinskii-Kosterlitz-Thouless transition, marked by a universal jump in the superfluid stiffness consistent with a condensate of charge 4e. Remarkably, the transition temperature Tc increases nearly linearly with interaction strength, providing a robust mechanism for high-Tc quartet superconductivity. Furthermore, spectral analysis reveals a prominent pseudogap above Tc arising from strong phase fluctuations. Our results establish a canonical and numerically exact model system for charge-4e superconductivity, offering crucial guidance for its realization in experimental platforms such as moiré materials and ultracold atomic systems.
Many features of charge-4e superconductors remain unknown because even the ``mean-field Hamiltonian'' describing them is an interacting model. Here we introduce an interacting model to describe a charge-4e superconductor (SC) deep in the superconducting phase and explore its properties using quantum Monte Carlo (QMC) simulations. The QMC is sign-problem-free but only when a Majorana representation is employed. As a function of the chemical potential we observe two sharply-distinct behaviors: a ``strong'' quarteting phase in which charge-4e quartets are tightly bound (like molecules) so that charge-2e pairing does not occur even in the temperature $T\ensuremath{\rightarrow}0$ limit, and a ``weak'' quarteting phase in which a further transition to a charge-2e superconducting phase occurs at a lower critical temperature. Analogous issues arise in a putative Z4 spin liquid with a pseudo-Fermi surface and other interacting models with composite order parameters. Under certain circumstances, we also identified a stable T=0 charge-4e SC phase with gapless nodal quasiparticles. We further discuss possible relevance of our results to various experimental observations in 81-doped LBCO.
Superconductivity is usually understood as a phase in which charge-2e Cooper pairs are condensed. Charge-4e superconductivity has largely been discussed as a vestigial order at finite temperature emerging from charge-2e states. Primary charge-4e superconducting phases at zero temperature remain scarce in both experiments and microscopic models. Here we argue that a doped featureless Mott insulator with SU(4) symmetry provides a natural platform for primary charge-4e superconductivity, based on perturbative renormalization group arguments and group theoretic considerations. As a concrete realization, we construct a bilayer Hubbard model with tunable onsite SU(4) and Sp(4) symmetries that exhibits a featureless Mott insulating phase at half filling. Its low energy physics is captured by a generalized ESD model, featuring an effective Hamiltonian that is purely kinetic within the constrained Hilbert space. Using density matrix renormalization group (DMRG) simulations, we find a primary charge-4e superconducting phase in the SU(4) ESD model and a conventional primary charge-2e phase in the Sp(4) case. We further characterize the corresponding normal states and discuss the resulting finite temperature phase diagram.
A direct charge-4e superconductor exhibits coherent four-electron order without condensation in any charge-2e pairing channel. We establish three complementary results. First, a particle-number-conserving Hamiltonian with electron hopping and two-body interactions moves bound quartets efficiently: for electron- and pair-hopping amplitudes t and K and local gap U0, virtual pair transfer produces quartet hopping at order K2/U0, whereas electron hopping alone contributes only at order t4/U03. Second, the microscopic theory identifies an exactly solvable parent Hamiltonian at the phase-separation boundary. Its ground state combines quartet off-diagonal long-range order and finite phase stiffness without pair condensation, while rigorous bounds establish quartet order beyond the parent in a family with explicit quartet transfer. Third, a single convergent strong-coupling expansion simultaneously controls quartet kinetics and thermodynamic stability and rules out condensation in every pair channel of the microscopic model. On bipartite hypercubic lattices at half quartet filling, it yields a parametrically controlled homogeneous phase for d≥2, with quartet order, a linear collective mode, and gapped charge-2e excitations. Virtual pair transfer thus enables a direct charge-4e superconducting ground state without pair condensation.
We consider a SU(2) lattice gauge theory on the square lattice, with a single fundamental complex fermion and a single fundamental complex boson on each lattice site. Projective symmetries of the gauge-charged fermions are chosen so that they match with those of the spinons of the π -flux spin liquid. Global symmetries of all gauge-invariant observables are chosen to match with those of the particle-hole symmetric electronic Hubbard model at half-filling. Consequently, both the fundamental fermion and fundamental boson move in an average background π -flux, their gauge-invariant composite is the physical electron, and eliminating gauge fields in a strong gauge-coupling expansion yields an effective extended Hubbard model for the electrons. The SU(2) gauge theory displays several confining/Higgs phases: a nodal d -wave superconductor, and states with Néel, valence-bond solid, charge, or staggered current orders. There are also a number of quantum phase transitions between these phases that are very likely described by ( 2 + 1 ) -dimensional deconfined conformal gauge theories, and we present large flavor expansions for such theories. These include the phenomenologically attractive case of a transition between a conventional insulator with a charge gap and Néel order, and a conventional d -wave superconductor with gapless Bogoliubov quasiparticles at four nodal points in the Brillouin zone. We also apply our approach to the honeycomb lattice, where we find a bicritical point at the junction of Néel, valence bond solid (Kekulé), and Dirac semimetal phases.
The low-temperature properties of the two-dimensional attractive Hubbard model are strongly influenced by the fermion density. Away from half-filling, there is a finite-temperature transition to a phase with s-wave pairing order. However, Tc is suppressed to zero at half-filling, where long-range charge-density-wave order also appears, degenerate with superconductivity. This paper presents determinant quantum Monte Carlo simulations of the attractive Hubbard model in the presence of a confining potential Vtrap which makes the fermion density $\ensuremath{\rho}$ inhomogeneous across the lattice. Pair correlations are shown to be large at low temperatures in regions of the trapped system with incommensurate filling, and to exhibit a minimum as the local density $\ensuremath{\rho}(\mathbit{i})$ passes through one fermion per site. In this ring of $\ensuremath{\rho}(\mathbit{i})=1$, charge order is enhanced. A comparison is made between treating Vtrap within the local-density approximation (LDA) and in an ab initio manner. It is argued that certain sharp features of the LDA result at integer filling do not survive the proximity of doped sites. The critical temperature of confined systems of fixed characteristic density is estimated.