Identifying superconducting states of matter without prior assumptions is a central challenge in strongly correlated electron systems. We introduce a canonical framework for diagnosing the formation of Cooper pair condensates based on the Penrose–Onsager criterion, in which superconducting order is encoded in the spectral properties of the two-particle reduced density matrix (2RDM). Within this formulation, the symmetry and structure of the condensate are obtained by projecting the 2RDM onto irreducible representations of the underlying symmetry group, enabling an unbiased identification of both conventional and exotic superconducting states. We demonstrate the power and versatility of the approach through applications to the two-dimensional Hubbard model, using both auxiliary-field quantum Monte Carlo (AFQMC) and the density matrix renormalization group (DMRG). For attractive interactions without a magnetic field, we reveal a clear finite-size scaling of the condensate fraction on square lattices of size up to 20 × 20 . The framework further provides direct access to the internal structure and extent of Cooper pairs, which we track across the BCS–BEC crossover. Moreover, it enables a clean diagnosis of the finite-momentum Fulde–Ferrell–Larkin–Ovchinnikov (FFLO) phase in a magnetic field. Finally, we apply the approach to a supersolid phase in the repulsive Hubbard model with an additional next-nearest-neighbor hopping t ′ , where a charge-density wave coexists with superconductivity. We confirm the fragmented nature of the condensate and uncover substantial pairing correlations in the triplet channel with p -wave spatial symmetry in addition to the dominant singlet d -wave pairing. Our results establish the 2RDM-based Penrose–Onsager framework as a broadly applicable and unbiased tool for characterizing superconducting order in correlated quantum matter.
Identifying superconducting states of matter without prior assumptions is a central challenge in strongly correlated electron systems. We introduce a canonical framework for diagnosing the formation of Cooper pair condensates based on the Penrose–Onsager criterion, in which superconducting order is encoded in the spectral properties of the two-particle reduced density matrix (2RDM). Within this formulation, the symmetry and structure of the condensate are obtained by projecting the 2RDM onto irreducible representations of the underlying symmetry group, enabling an unbiased identification of both conventional and exotic superconducting states. We demonstrate the power and versatility of the approach through applications to the two-dimensional Hubbard model, using both auxiliary-field quantum Monte Carlo (AFQMC) and the density matrix renormalization group (DMRG). For attractive interactions without a magnetic field, we reveal a clear finite-size scaling of the condensate fraction on square lattices of size up to 20 × 20 . The framework further provides direct access to the internal structure and extent of Cooper pairs, which we track across the BCS–BEC crossover. Moreover, it enables a clean diagnosis of the finite-momentum Fulde–Ferrell–Larkin–Ovchinnikov (FFLO) phase in a magnetic field. Finally, we apply the approach to a supersolid phase in the repulsive Hubbard model with an additional next-nearest-neighbor hopping t ′ , where a charge-density wave coexists with superconductivity. We confirm the fragmented nature of the condensate and uncover substantial pairing correlations in the triplet channel with p -wave spatial symmetry in addition to the dominant singlet d -wave pairing. Our results establish the 2RDM-based Penrose–Onsager framework as a broadly applicable and unbiased tool for characterizing superconducting order in correlated quantum matter.
We report a systematic study of the superconducting ground state in the two-dimensional attractive Hubbard model, employing the {\it numerically exact} auxiliary-field quantum Monte Carlo method. We focus primarily on the pairing order and condensate fraction across the BCS-BEC crossover. Our numerical calculations cover a broad and physically relevant parameter regime, with the interaction strength U/t=2∼12 (where t is the nearest-neighbor hopping amplitude) and fermion density n=0.25∼1.0. We systematically characterize the dependence of the pairing order and condensate fraction on both the interaction strength and fermion density, and perform finite-size extrapolations to obtain their thermodynamic-limit values. As complementary quantities, we also provide high-precision numerical results for the double occupancy and total energy. Furthermore, we investigate the effects of next-nearest-neighbor hopping on the superconducting properties of the model. Our results offer a comprehensive set of unbiased numerical data for the attractive Hubbard model, establishing quantitative benchmarks for understanding superconductivity across the BCS-BEC crossover and for validating future theoretical and numerical approaches.
Condensation of bosons in Bose-Einstein condensates or Cooper pairs in superconductors refers to a macroscopic occupation of a few single- or two-particle states. A condensate is called "fragmented" if not a single, but multiple states are macroscopically occupied. While fragmentation is known to occur in particular Bose-Einstein condensates, we propose that fragmentation naturally takes place in striped superconductors. To this end, we investigate the nature of the superconducting ground state realized in the two-dimensional t-t^{'}-J model. In the presence of charge density modulations, the condensate is shown to be fragmented and composed of partial condensates located on the stripes. The fragments of the condensates hybridize to form an extended macroscopic wave function across the system. The results are obtained from evaluating the singlet-pairing two-particle density matrix of the ground state on finite cylinders computed via the density matrix renormalization group method. Our results shed light on the intricate relation between stripe order and superconductivity in systems of strongly correlated electrons.
d-wave FFLO state and charge-2e supersolidity in the t-t′-J model under Zeeman fields
Unconventional superconductivity under strong Zeeman fields--particularly beyond the Pauli paramagnetic limit--remains a central challenge in condensed matter physics. The exotic Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, in particular, remains in need of definitive study within fundamental electronic models. Here we employ state-of-the-art finite-temperature and ground-state tensor network approaches to systematically explore the superconducting (SC) phase diagram of the t-t′-J model subjected to Zeeman fields. We find that zero-momentum d-wave superconductivity persists until the spin gap closes, coexisting with charge density waves. A novel d-wave FFLO phase emerges under a higher Zeeman field even above the Pauli limit, concomitant with a field-enhanced spin density waves. We identify these phases, characterized by the simultaneous presence of pairing condensate and density wave orders, as charge-2e supersolids. Analysis of Matsubara Green's function reveals that the FFLO pairing momentum is locked to the underlying Fermi surface. Our results provide microscopic insights into field-induced unconventional pairing mechanisms and reveal the long-sought FFLO state in a fundamental correlated electron model, offering a promising route for its realization in ultracold atom optical lattice.
We investigate odd-in-time—or —pairing of fermions in equilibrium systems within the particle-number-conserving framework of Penrose, Onsager, and Yang, where superfluid order is defined by macroscopic eigenvalues of reduced density matrices. We show that odd-frequency pair correlations are synonymous with even fermion-exchange symmetry in a time-dependent correlation function that generalises the two-body reduced density matrix. Macroscopic even-under-fermion-exchange pairing is found to emerge from conventional Penrose-Onsager-Yang condensation in two-body or higher-order reduced density matrices through the symmetry-mixing properties of the Hamiltonian. We identify and characterize a matrix responsible for producing macroscopic even fermion-exchange correlations that coexist with a conventional Cooper-pair condensate, while a matrix is shown to be responsible for creating macroscopic even fermion-exchange correlations from hidden orders such as a multiparticle condensate. The transformer scenario is illustrated using the spin-balanced s -wave superfluid with Zeeman splitting as an example. The generator scenario is demonstrated by the composite-boson condensate arising for itinerant electrons coupled to magnetic excitations. Structural analysis of the transformer and generator matrices is shown to provide general conditions for odd-frequency pairing order to arise in a given system. Our formalism facilitates a fully general derivation of the Meissner effect for odd-frequency superconductors that holds also beyond the regime of validity for mean-field theory. Published by the American Physical Society 2024
The following thesis is divided in two main parts. Chapters 2, 3 and 4 are devoted to the study of the so called pair-density-wave (PDW) superconducting state and some of its connections to electronic liquid crystal (ELC) phases, its topological aspects in a one dimensional model and its appearance in a quasi-one dimensional system. On the other hand, chapter 5 is focused on the investigation of the classical statistical mechanics properties of dimers, in particular, the dimer model on the Aztec diamond graph and its relation with the octahedron equation. In chapter 2 we present a theory of superconducting states where the Cooper pairs have a nonzero center-of-mass momentum, inhomogeneous superconducting states known as a pair-density-waves (PDWs) states. We show that in a system of spin-1/2 fermions in two dimensions in an electronic nematic spin-triplet phase where rotational symmetry is broken in both real and spin space PDW phases arise naturally in a theory that can be analysed using controlled approximations. We show that several superfluid phases that may arise in this phase can be treated within a controlled BCS mean field theory, with the strength of the spin-triplet nematic order parameter playing the role of the small parameter of this theory. We find that in a spin-triplet nematic phase, in addition to a triplet p-wave and spin-singlet d-wave (or s depending on the nematic phase) uniform superconducting states, it is also possible to have a d-wave (or s) PDW superconductor. The PDW phases found here can be either unidirectional, bidirectional, or tridirectional depending on the spin-triplet nematic phase and which superconducting channel is dominant. In addition, a triple-helix state is found in a particular channel. We show that these PDW phases are present in the weak-coupling limit, in contrast to the usual Fulde-Ferrell-Larkin-Ovchinnikov phases, which require strong coupling physics in addition to a large magnetic field (and often both). In chapter 3 we show that the pair-density-wave (PDW) superconducting state emergent in extended Heisenberg-Hubbard models in two-leg ladders is topological in the presence of an Ising spin symmetry and supports a Majorana zero mode (MZM) at an open boundary and at a junction with a uniform d-wave one-dimensional superconductor. Similarly to a conventional finite-momentum paired state, the order parameter of the PDW state is a charge-2e field with finite momentum. However, the order parameter here is a quartic electron operator and conventional mean-field theory cannot be applied to study this state. We use bosonization to show that the 1D PDW state has a MZM at a boundary. This superconducting state is an exotic topological phase supporting Majorana fermions with finite-momentum pairing fields and charge-4e superconductivity. In chapter 4 we provide a quasi-one-dimensional model which can support a PDW state. The model consists of an array of strongly-interacting one-dimensional systems, where the one-dimensional systems are coupled to each other by local interactions.Within the interchain mean-field theory (MFT), we find several SC states from the model, including a conventional uniform SC state, PDW SC state, and a coexisting phase of the uniform SC and PDW states. In this quasi-1D regime we can treat the strong correlation physics essentially exactly using bosonization methods and the crossover to the 2D system by means of interchain MFT. The resulting critical temperatures of the SC phases generically exhibit a power-law scaling with the coupling constants of the array, instead of the essential singularity found in weak-coupling BCS-type theories. Electronic excitations with an open Fermi surface, which emerge from the electronic Luttinger liquid systems below their crossover temperature to the Fermi liquid, are then coupled to the SC order parameters via the proximity effect. From the Fermi surface thus coupled to the SC order parameters, we calculate the quasiparticle spectrum in detail. We show that the quasiparticle spectrum can be fully gapped or nodal in the uniform SC phase and also in the coexisting phase of the uniform SC and PDW parameters. In the pure PDW state, the excitation spectrum has a reconstructed Fermi surface in the form of Fermi pockets of Bogoliubov quasiparticles. In chapter 5 we study the octahedron relation (also known as the A∞ T-system), obeyed in particular by the partition function for dimer coverings of the Aztec Diamond graph. For a suitable class of doubly periodic initial conditions, we find exact solutions with a particularly simple factorized form. For these, we show that the density function that measures the average dimer occupation of a face of the Aztec graph, obeys a system of linear recursion relations with periodic coefficients. This allows us to explore the thermodynamic limit of the corresponding dimer models and to derive exact “arctic” curves separating the various phases of the system.
In a Fermi system near or at half-filling, a specific superconducting pairing channel, if not explicitly included in the Hamiltonian, can be boosted by suppressing a competing pairing channel; this is exemplified by the enhancement of extended s-wave correlations upon suppressing s-wave Cooper pairing. This phenomenon, originally found by the use of generalized uncertainty relations is referred to as order by projection. The case of zero on-site Coulomb interaction in the thermodynamic limit, confirms this mechanism through the analytical solution. In this study, we go further and systematically investigate this mechanism for a strongly correlated fermionic Hubbard model, now with finite on-site interaction, on a square lattice with an extended set of hopping parameters. We explore the behaviors of different pairing channels when one of them is suppressed, utilizing density matrix renormalization group calculations. Our findings provide numerical evidence supporting the existence of order by projection in the strongly correlated system we studied. We also investigate the effect of the strength of Hubbard U, next-nearest-neighbor ${t}^{\ensuremath{'}}$, hole doping, as well as finite-size scaling approaching the thermodynamic limit.