Thermal Phase Structure of the Attractive Fermi Hubbard Model with Imaginary Chemical Potential

This paper investigates the BCS-BEC crossover in the large-$N$ attractive Fermi-Hubbard model on a one-dimensional lattice using mean-field theory with an imaginary chemical potential, revealing that the crossover is governed by three key parameters and identifying a specific thermal window at unitarity where the superconducting gap vanishes while the fermion number exhibits extrema, signaling a transition between particle-like and hole-like dominant regimes.

The Gist

Imagine a crowded dance floor where thousands of tiny dancers (electrons) are trying to find partners. In the world of quantum physics, these dancers can pair up in two very different ways, depending on how "cold" the room is and how much they like each other.

• The BCS Style (The Cool, Distant Waltz): When the dancers are cold and don't interact much, they form pairs that are far apart. They move in a synchronized, elegant waltz. This is like superconductivity, where electricity flows with zero resistance.
• The BEC Style (The Hot, Tight Embrace): When the attraction is strong, the dancers grab each other tightly and form a tight-knit group, almost like a single giant blob. This is Bose-Einstein Condensation, a state of matter where particles act as one.

Usually, scientists study how to switch from the "Waltz" to the "Tight Embrace" by changing the temperature or the strength of the attraction. But this paper introduces a weird, invisible twist to the dance floor: an Imaginary Chemical Potential.

The Invisible Twist: The "Ghost" Music

In this study, the author adds a special, invisible rule to the dance floor. Imagine the music isn't just playing; it's also spinning the dancers in a circle as they move. This isn't a physical force you can touch; it's more like a ghostly phase shift. In physics, we call this an "imaginary chemical potential."

Think of it like a twisted rubber band connecting the start and end of the dance floor. As a dancer completes a lap, they don't just return to the start; they return slightly "out of sync" with where they began, as if they stepped into a different dimension for a split second.

The Discovery: The "Thermal Window"

The author discovered that when you tune this invisible twist to two very specific angles (mathematically, 120 degrees and 240 degrees, or $2\pi/3$ and $4\pi/3$), something magical happens.

At these specific angles, the dance floor enters a "Thermal Window."

1. The Gap Vanishes: Usually, to keep the dancers paired up, you need a certain amount of energy (a "gap") to keep them from falling apart. But at these special angles, that gap disappears completely. The dancers are in a state of pure, chaotic potential.
2. The Balance Point: At the same time, the number of dancers moving forward (particles) versus backward (holes) hits a perfect peak or valley. It's as if the dance floor is holding its breath, perfectly balanced between two states.

The Three Control Knobs

The paper explains that this whole phenomenon is controlled by three main knobs:

1. The Twist ($\theta$): How much we spin the dancers (the imaginary chemical potential).
2. The Temperature ($T$): How hot or cold the dance floor is.
3. The Attraction ($\delta u$): How much the dancers want to hold hands.

The author found that when you turn the "Twist" knob to exactly 120 degrees or 240 degrees, the system becomes incredibly sensitive. It's like a seesaw that is perfectly balanced; a tiny nudge in the temperature or attraction will instantly flip the dancers from the "Waltz" style to the "Tight Embrace" style.

Why is this important?

You might wonder, "Who cares about imaginary numbers in a dance?"

Well, these "imaginary" twists are actually real tools used in modern physics to simulate complex systems, like:

• Topological materials: Materials that conduct electricity on their surface but are insulators inside.
• Ultracold atoms: Scientists use lasers to create "synthetic" magnetic fields that act just like this imaginary twist.
• High-temperature superconductors: Understanding how these pairs form could help us build better electronics.

The Big Picture

The author is essentially saying: "We found a secret door in the quantum world."

By tuning the invisible "twist" in the system to these specific angles ($120^\circ$ and $240^\circ$), we can create a special zone where the rules of pairing change. In this zone, the system is undecided, hovering right on the edge between being a superconductor and a condensate.

It's like finding a specific frequency on a radio where the static clears up, and you can suddenly hear the music of the universe a little more clearly. This paper gives us a new map to navigate the strange, twisted landscape of quantum matter, showing us exactly where to look to understand how particles pair up in the most extreme conditions.

Technical Summary

1. Problem Statement

The paper investigates the BCS–BEC crossover (the transition from weakly coupled Bardeen-Cooper-Schrieffer superconductivity to strongly coupled Bose-Einstein condensation) within the context of the large-$N$ attractive Fermi-Hubbard model on a one-dimensional lattice.

The central challenge addressed is understanding how this crossover behaves in the presence of an imaginary chemical potential ($\mu = i\theta$). While imaginary chemical potentials are not directly observable in standard lattice theories, they arise naturally in the canonical ensemble when constraining the fermion number. Physically, this corresponds to a background $U(1)$ gauge field creating a holonomy (an Aharonov-Bohm phase) along the thermal circle. The author seeks to determine how this "temporal gauge field" influences the pairing gap, the fermion number, and the phase boundaries between BCS and BEC regimes.

2. Methodology

The study employs the following theoretical framework and approximations:

• Model: The spin-1/2 attractive Fermi-Hubbard model with interaction strength $U > 0$.
• Large-$N$ Expansion: The fermions are extended to $N$ flavors ($\alpha = 1, \dots, N$) to allow for a controlled saddle-point approximation.
• Path Integral Formalism: The canonical partition function is expressed using a projector on the fermion number, introducing the imaginary chemical potential $\mu = i\theta$.
• Hubbard-Stratonovich Transformation: The quartic interaction is decoupled in the pairing channel by introducing a complex bosonic field $\Delta$ (the superconducting gap).
• Mean-Field Approximation: The fermionic fields are integrated out to derive an effective action $S_{eff}[\Delta, \theta]$. The static saddle-point configuration is assumed ($\Delta_i(\tau) = \Delta$).
• Gap and Number Equations: The authors derive the gap equation (minimizing $S_{eff}$ with respect to $\Delta^*$) and the number equation (minimizing with respect to $\theta$) in the thermodynamic limit.
• Renormalization: To handle logarithmic divergences at the Brillouin zone edges, a renormalized gap equation is formulated by subtracting the critical coupling $1/U_c$.

3. Key Contributions and Theoretical Framework

A. The Thermal Kernel and Control Parameters

The analysis identifies three governing parameters for the crossover:

1. Coupling Parameter ($\delta_u$): Defined as $\delta_u = \frac{1}{U} - \frac{1}{U_c}$.
   • $\delta_u > 0$: BCS regime.
   • $\delta_u < 0$: BEC regime.
   • $\delta_u = 0$: Unitarity point.
2. Imaginary Chemical Potential ($\theta$): Represented dimensionlessly as $\phi = \beta\theta$.
3. Thermal Kernel ($g$): A function $g(x, \phi) = \frac{\sinh x}{\cosh x + \cos \phi}$ (where $x = \beta E_k$) that dictates the contribution of thermal fluctuations to the gap equation.

B. Discovery of the "Thermal Window"

The most significant theoretical contribution is the identification of a "thermal window" at specific angles:
$$\phi = \beta\theta = \frac{2\pi}{3} \quad \text{and} \quad \frac{4\pi}{3}$$
At these angles, $\cos \phi = -1/2$. The authors demonstrate that at the unitarity point ($\delta_u = 0$), the superconducting gap $\Delta$ vanishes, and the fermion number $N_f$ reaches local extrema.

C. Mechanism of the Thermal Window

The condition $\cos \phi = -1/2$ causes the thermal kernel $g(x, \phi) - 1$ to change sign at a specific energy scale $x^* = \ln 2$.

• For low energies ($x < \ln 2$), the kernel suppresses pairing ($g < 1$).
• For high energies ($x > \ln 2$), the kernel enhances pairing ($g > 1$).
  At the critical coupling, these competing contributions balance exactly, leading to a vanishing gap. This represents a redistribution of spectral weight controlled by the imaginary chemical potential.

4. Key Results

A. Phase Structure and Gap Behavior

• Vanishing Gap: At unitarity ($U=U_c$) and $\phi = 2\pi/3, 4\pi/3$, the gap equation yields $\Delta = 0$.
• Phase Diagrams: Numerical plots (Figures 1–3) show that the boundaries of the thermal window are universal; they remain fixed at $2\pi/3$ and $4\pi/3$ regardless of the coupling strength $\delta_u$.
• Regime Selection: Inside this thermal window, the system is highly sensitive to coupling changes, allowing a selection between BCS and BEC dominance based on the sign of $\delta_u$.

B. Fermion Number ($N_f$) Extrema

The fermion number $N_f$, which quantifies the balance between particle-like and hole-like excitations, exhibits critical behavior at the special angles:

• At $\phi = 2\pi/3$: $N_f$ reaches a local maximum. This corresponds to a dominance of particle-like excitations.
• At $\phi = 4\pi/3$: $N_f$ reaches a local minimum. This corresponds to a dominance of hole-like excitations.
• Curvature: The second derivative of $N_f$ with respect to $\phi$ changes sign at these points, confirming them as extrema rather than inflection points.

C. Analytical Limits for $N_f(T)$

The paper derives exact analytical expressions for the fermion number as a function of the temperature-to-hopping ratio ($T/t$):

1. Low Temperature Limit ($T \ll t$):
   $$N_f\left(\frac{2\pi}{3}, T\right) \approx \frac{2}{3} \frac{T}{t}$$
   The fermion number grows linearly with temperature.

2. High Temperature Limit ($T \gg t$):
   $$N_f\left(\frac{2\pi}{3}, T\right) \approx \sqrt{3} \left[ 1 - 2\left(\frac{t}{T}\right)^2 + \dots \right]$$
   As $T \to \infty$, $N_f$ approaches a constant value of $\sqrt{3} \approx 1.732$, independent of material parameters. This indicates that the imaginary chemical potential induces a fixed particle-hole asymmetry even at infinite temperature.

3. Crossover Point:
   The transition between the linear low-$T$ regime and the constant high-$T$ regime occurs numerically at $T/t \approx 2$.

5. Significance and Implications

• Universal Angles: The discovery that $\phi = 2\pi/3$ and $4\pi/3$ are universal boundaries for the thermal window suggests a deep algebraic structure (related to the primitive third roots of unity) organizing the BCS, Unitarity, and BEC regimes.
• Connection to Topology: The condition $\cos \phi = -1/2$ appears in other contexts, such as topological band theories with staggered magnetic fluxes and Floquet-engineered gauge fields. This work bridges the gap between lattice many-body physics and these topological phenomena.
• Experimental Relevance: The results provide a theoretical framework for ultracold atomic gas experiments using synthetic gauge fields or Floquet engineering, where such imaginary chemical potentials (or equivalent phases) can be realized.
• New Window into Pairing: The study reveals that thermal fluctuations can actively compete with interaction scales to control the stability of the normal state, offering a new perspective on how pairing correlations emerge in lattice systems.

Conclusion

The paper establishes that the BCS-BEC crossover in the attractive Fermi-Hubbard model is governed by a delicate interplay between coupling strength, temperature, and an imaginary chemical potential. The identification of the "thermal window" at $\beta\theta = 2\pi/3, 4\pi/3$ provides a precise theoretical prediction for where the superconducting gap vanishes and fermion number asymmetry is maximized, offering a new tool for understanding and manipulating quantum many-body phases.