BCS-BEC crossover in trapped one-dimensional Fermi-Hubbard chains: entanglement and correlation signatures from DMRG and effective-pairing theory

This paper characterizes the BCS-BEC crossover in harmonically confined one-dimensional Fermi-Hubbard chains by combining DMRG simulations with effective-pairing theory to reveal how spatial confinement reshapes correlation patterns, leading to coexisting insulating and superfluid regions distinguishable via conditioned correlation functions.

The Gist

Imagine a long, narrow hallway where tiny, invisible dancers (electrons) are trying to move around. In a perfect, endless hallway, these dancers follow strict rules: sometimes they dance alone, and sometimes they pair up to waltz together. Physicists call this the "BCS-BEC crossover." It's a spectrum where the dancers go from being loosely connected partners (BCS) to being tightly glued together as a single unit (BEC).

But in the real world, hallways aren't endless; they have walls. In this paper, the researchers study what happens when these dancers are trapped in a curved hallway (a harmonic trap) that gets narrower in the middle and wider at the ends. This confinement changes everything.

Here is the story of their findings, explained simply:

1. The Setup: A Crowded, Curved Dance Floor

The researchers used a super-powerful computer simulation (called DMRG) to watch these electrons. They also built simpler "toy models" (effective theories) to understand the physics without getting lost in the math.

• The Trap: Imagine the hallway is shaped like a bowl. The dancers naturally want to sit in the deepest part (the center).
• The Interaction: The dancers can either ignore each other, push each other away, or be strongly attracted to each other. The researchers cranked up the "attraction" to see how the pairs formed.

2. The Two Extreme Dances

The paper explores two main ways the electrons behave:

• The "Loose Waltz" (BCS Regime): When the attraction is weak, the electrons form pairs, but they are like long-distance partners holding hands across the room. They are spread out and move somewhat independently.
• The "Glued Twins" (BEC Regime): When the attraction is very strong, the electrons snap together so tightly they act like a single, heavy object. They are glued to the same spot.

3. The Surprise: The "Insulating Core" and "Superfluid Wings"

In a normal, endless hallway, the whole floor would behave the same way. But because of the curved trap, the paper discovered a strange, split personality in the system:

• The Center (The Insulator): As the hallway gets crowded, the dancers in the very center get so packed together that they stop moving entirely. They freeze into a solid block. The researchers call this an insulating region. It's like a traffic jam where no one can move.
• The Edges (The Superfluid): Here is the magic. Even though the center is frozen, the dancers at the edges of the hallway keep dancing freely. They form a "superfluid" (a frictionless flow).
• The Result: You get a sandwich: a frozen, stuck core surrounded by a flowing, dancing shell. The paper calls this a composite INS+SF phase.

4. How They Spotted the Difference

How do you tell if the dancers are doing a "Loose Waltz" or acting as "Glued Twins"? The researchers invented a new way to look at the data:

• The "RMS Distance" (The Size of the Pair): They measured how far apart the two dancers in a pair usually are.

  • In the BCS mode, the pair is huge (like holding hands across the room).
  • In the BEC mode, the pair is tiny (glued at the same spot).
  • By watching how this distance shrinks as they turned up the attraction, they could clearly see the transition from one dance style to the other.

• The "Entanglement" (The Connection): They also looked at how "connected" the left half of the hallway is to the right half.

  • When the center freezes (becomes an insulator), the connection between the left and right sides suddenly snaps. It's like cutting a bridge; the two sides can no longer "talk" to each other. This sudden snap tells them exactly when the insulating core forms.

5. Why the Center Freezes

Why does the middle get stuck?

• The "Effective" Trap: When the electrons are glued together (BEC), they act like heavy bosons. The researchers found that the trap effectively feels stronger to these glued pairs. It's as if the bowl gets deeper and steeper for the pairs than for single dancers.
• The Repulsion: Even though the pairs are attracted to each other, the "glued" nature of the BEC pairs makes them repel their neighbors slightly. This pushes them away from the center, creating a weird oscillation where the density goes up and down near the edges of the frozen core.

Summary of the Discovery

The paper shows that when you trap these quantum dancers in a curved space:

1. Strong attraction makes them glue together (BEC).
2. Crowding makes the center freeze into a solid block (Insulator).
3. The edges stay fluid and dancing (Superfluid).
4. The transition between "loose pairs" and "glued pairs" isn't just a smooth slide; it leaves a clear fingerprint in how the pairs are sized and how connected the system is.

The researchers successfully mapped out exactly where these different behaviors happen, creating a "map" (phase diagram) that tells you: "If you have this much crowding and this much attraction, you will get a frozen center with dancing wings." They proved that their simple "toy models" matched perfectly with their complex computer simulations, giving them a unified picture of how quantum matter behaves when squeezed into a trap.

Technical Summary

Technical Summary: BCS-BEC Crossover in Trapped One-Dimensional Fermi–Hubbard Chains

Problem Statement
The Bardeen–Cooper–Schrieffer (BCS) to Bose–Einstein condensation (BEC) crossover is a well-established phenomenon in homogeneous Fermi–Hubbard systems, characterized by the evolution from weakly bound Cooper pairs to tightly bound local pairs as interaction strength increases. However, experimental realizations using ultracold atoms in optical lattices often involve spatial confinement (e.g., harmonic traps), which breaks translational symmetry. This inhomogeneity complicates the identification of the BCS-BEC crossover, as standard probes based on pair correlation lengths, momentum distributions, or the Leggett criterion become ambiguous. In trapped one-dimensional (1D) systems, the loss of integrability prevents exact analytical solutions, and existing numerical or mean-field approaches often struggle to provide physically intuitive criteria for classifying pairing regimes in strongly inhomogeneous many-body systems.

Methodology
The authors investigate the BCS-BEC crossover in a 1D harmonically confined Fermi–Hubbard chain at zero temperature using a dual approach:

1. Numerical Simulations: They employ the Density Matrix Renormalization Group (DMRG) method in its matrix product state (MPS) formulation. Simulations are performed for spin-balanced systems with open boundary conditions, utilizing the ITensor library. The accuracy is controlled by the bond dimension, with initial values set to $\chi_{max}=1500$ and progressively increased.
2. Effective Models: To provide physical intuition and analytical benchmarks, the authors derive effective models for different interaction regimes:
   • Strong Coupling ($|U| \gg t$): Using second-order perturbation theory (Löwdin perturbation theory), they map the system to a confined tight-binding model of hard-core bosons (representing tightly bound on-site pairs). This model includes renormalized hopping ($\tau \approx 2t^2/|U|$), an enhanced effective confinement strength ($\tilde{k} = k|U|/2t$), and nearest-neighbor repulsion.
   • Weak to Intermediate Coupling ($|U| \lesssim 2t$): A mean-field (MF) approximation is used, decoupling the interaction term to solve for Bogoliubov quasiparticles self-consistently.
   • Non-interacting Limit ($U=0$): Analytical solutions for a confined tight-binding model are used to establish reference points for the spectrum and density profiles.

Key Contributions and Results

• Phase Diagram and Transitions: The study constructs a phase diagram as a function of interaction strength $|U|$ and average density $n$. It identifies three primary regimes:

  • BEC Regime: Dominated by tightly bound pairs, occurring at low densities or very strong interactions.
  • BCS Regime: Characterized by extended Cooper pairs, occurring at intermediate densities and weaker interactions.
  • Composite INS+SF Phase: At high densities, the trap center becomes insulating (INS) due to the formation of localized states, while superfluid (SF) states persist at the edges. This "superfluid wings" phenomenon is robust across interaction strengths.

• New Diagnostic Tools:

  • Conditioned Correlation Function: The authors introduce a four-particle correlation function, $P_j(r) = \langle c^\dagger_{j+r,\uparrow} c_{j,\uparrow} n_{j,\downarrow} \rangle$, which probes the local binding between opposite spins. They demonstrate that the root-mean-square distance ($r_{rms}$) of this distribution exhibits two distinct power-law scalings: a slow decay ($\nu_1 \sim 0.1$) in the BCS regime and a rapid decay ($\nu_2 > 1$) in the BEC regime. The crossover is identified by the transition between these scaling behaviors.
  • Entanglement-Based Probes: The half-chain entanglement entropy ($S_{L/2}$) is shown to serve a dual purpose. A sharp drop in $S_{L/2}$ signals the superfluid-to-insulator transition (formation of a central insulating barrier). Furthermore, the magnitude of $S_{L/2}$ relative to a reference value ($2 \ln 2$) distinguishes between BCS-like (weakly correlated, $S_{L/2} \ge 2 \ln 2$) and BEC-like (strongly bound pairs, $S_{L/2} < 2 \ln 2$) regimes.

• Physical Mechanisms:

  • Edge Superfluidity: The effective models reveal that in the strong-coupling limit, effective repulsion between BEC-like quasiparticles (arising from the breaking of translational symmetry by the trap) pushes pairs away from the trap center. This mechanism explains the persistence of superfluid states at the edges even when the center is insulating.
  • Confinement Enhancement: The study confirms that attractive interactions effectively enhance the trapping potential ($\tilde{k} > k$), leading to higher central densities and shifting the superfluid-insulator transition to lower densities as $|U|$ increases.
  • Unusual BCS Pairs: In the intermediate regime, mean-field analysis reveals that BCS-like pairs can avoid the trap center, localizing at the interface between empty sites and the central region, consistent with the density oscillations observed in DMRG results.

Significance
The paper claims to provide a unified picture of the BCS-BEC crossover in harmonically confined geometries by combining robust numerical results with effective theoretical descriptions. Its primary significance lies in:

1. Resolving Identification Challenges: It offers new, physically transparent criteria (based on power-law scaling of correlation functions and entanglement entropy) to distinguish between BCS and BEC regimes in inhomogeneous systems where standard probes fail.
2. Understanding Edge States: It provides a microscopic explanation for the robustness of superfluid states at the edges of trapped chains, attributing them to effective repulsion between pairs and the specific spatial structure of pairing in confined potentials.
3. Consistency: The work demonstrates qualitative and quantitative consistency between DMRG simulations and effective models (including mean-field and hard-core boson descriptions), validating the use of these simplified models to capture essential physics in the strong-coupling limit.
4. Experimental Relevance: The findings are presented as qualitatively consistent with existing experimental observations of local occupation probabilities in quasi-1D harmonic traps and recent 2D experiments, particularly regarding the suppression of unpaired atoms and the emergence of charge-density-wave-like modulations.

The authors conclude that the features identified originate primarily from the breaking of translational symmetry rather than the specific form of the harmonic confinement, suggesting that their qualitative findings may persist in more general trapping potentials.