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.
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.
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.
Fermionic condensation typically occurs via pairing. In recent decades, however, a fundamental question has emerged: whether alternative forms of order exist, such as condensates of fermion quadruplets. These states--including ``charge-4e" superconductors and ``charge-0" counterflow condensates--lie beyond the standard Bardeen-Cooper-Schrieffer framework, and require strong fluctuations and correlation effects that invalidate the BCS mean-field description. This makes the problem notoriously difficult to study numerically at a microscopic level, as it involves both strong interactions and the fermionic sign problem. Here, we present a microscopic fermionic model featuring correlated hopping that significantly mitigates the sign problem, enabling rigorous Monte-Carlo-based analysis. Using large-scale simulations, we demonstrate the existence of a fermion-quadrupling condensate with a transition temperature comparable to the hopping energy scale. These results provide direct numerical evidence for quartic fermionic order in a microscopic system and suggest that these exotic states are also experimentally accessible in ultracold atomic gases.
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.
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.
This paper addresses the transition from the normal to the superfluid state in strongly correlated two dimensional fermionic superconductors and Fermi gases. We arrive at the Berezinskii-Kosterlitz-Thouless (BKT) temperature TBKT as a function of \emph{attractive} pairing strength by associating it with the onset of quasi-condensation in the normal phase. Our approach builds on a criterion for determining the BKT transition temperature for atomic gases which is based on a well established Quantum Monte Carlo analysis of the phase space density. This latter quantity, when derived from BCS-BEC crossover theory for fermions, leads to non-monotonic behavior for TBKT as a function of the attractive interaction or inverse scattering length. In Fermi gases, this implies a robust superconducting dome followed by a long tail from the flat BEC asymptote, rather similar to what is observed experimentally. For lattice systems we find that TBKT has an absolute maximum of the order of 0.1EF. We discuss how our results compare with those derived from the Nelson Kosterlitz criterion based on the mean field superfluid density and the approach to the transition from below. While there is agreement in the strict mean-field BCS regime at weak coupling, we find that at moderate pairing strength bosonic excitations cause a substantial increase in TBKT followed by an often dramatic decrease before the system enters the BEC regime.
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.