Observe novel tricritical phenomena in self-organized Fermi gas induced by higher order Fermi surface nesting

This paper investigates 1D fermionic superradiance in optical lattices, revealing novel tricritical phenomena and multistability arising from higher-order Fermi surface nesting that distinguish 1D from 2D systems and offer new insights into the relationship between quantum and classical phase transitions.

The Gist

Imagine a crowded dance floor where the dancers are tiny, invisible particles called fermions (like electrons). Usually, these particles are very shy; they refuse to stand next to each other or do the same move at the same time. This is a rule of nature called the "Pauli Exclusion Principle."

Now, imagine this dance floor is inside a giant, mirrored box (an optical cavity). We shine a bright laser light on them, and the mirrors bounce the light back and forth, creating a standing wave of light.

Here is the story of what happens when these shy dancers interact with the light, told in simple terms:

1. The "Super-Flash" (Superradiance)

Normally, if you shine a light on a crowd, they might just scatter the light randomly. But in this experiment, something magical happens. If the light is strong enough, the dancers suddenly stop acting individually. They all decide to move in perfect unison, creating a giant, synchronized wave.

When they do this, they dump all their energy into the light, creating a massive, blinding flash of photons (light particles) inside the box. This is called Superradiance. It's like a choir where everyone suddenly decides to sing the exact same note at the exact same volume, creating a sound so powerful it shakes the walls.

2. The "Perfect Fit" (Fermi Surface Nesting)

Why does this happen? The paper explains it using a concept called Fermi Surface Nesting.

Imagine the dancers are arranged in a circle. If the light wave has a specific rhythm, it might be able to push a dancer from one side of the circle to the exact opposite side with zero effort. It's like a "perfect fit" or a "magic shortcut."

• In 2D (a flat dance floor): The dancers have too much room to wiggle. Even if the light tries to push them, they can just dodge it slightly. The "perfect fit" is rare, so the synchronized flash is hard to trigger in a special way.
• In 1D (a narrow hallway): The dancers are stuck in a line. They cannot dodge. If the light pushes them, they must jump to the perfect spot on the other side. This "perfect fit" happens everywhere, making the system extremely sensitive.

3. The "Tipping Point" (Tricritical Phenomena)

The paper discovers a very special "tipping point" called a Tricritical Point.

Think of a seesaw.

• Normal Critical Point: Usually, if you add a little weight (light), the seesaw tips slowly and smoothly. The dancers gradually start dancing together. This is a "second-order" transition.
• Tricritical Point: But at this special spot, the seesaw behaves differently. If you add just a tiny bit more weight, the seesaw doesn't just tip; it snaps violently to the other side. The dancers don't gradually sync up; they suddenly jump into a new, chaotic state.

The paper shows that in the narrow hallway (1D), this "snapping" behavior exists because of the "perfect fit" (nesting) and the fact that the dancers are packed so tight they create a mathematical "infinity" (infrared divergence) that doesn't happen in the flat room (2D).

4. The "Stuck Door" (Multistability and Hysteresis)

The researchers also found that the system can get "stuck."

Imagine a door that can be either open or closed.

• If you push the door open gently, it stays open.
• If you push it closed, it stays closed.
• But in the middle, there is a "gray zone" where the door could be either open or closed, depending on how you pushed it last.

This is called Multistability. If you slowly increase the light and then slowly decrease it, the system doesn't retrace its steps. It follows a different path, creating a loop called Hysteresis. It's like a thermostat that turns the heat on at 70°F but doesn't turn it off until 60°F. The system has a "memory" of where it came from.

5. The "Goldilocks Temperature"

Finally, the team looked at what happens when the room gets warm (finite temperature).

• You might think that heat would ruin the synchronized dance, making it impossible.
• Surprisingly, they found that there is a "Goldilocks Temperature." It's not too cold (where quantum rules are strict) and not too hot (where chaos rules). There is a specific, warm-ish temperature where the synchronized flash is actually easiest to trigger.

The Big Takeaway

This paper is like a map for a new kind of quantum playground. It tells us:

1. Shape matters: How you arrange the particles (1D vs. 2D) changes the rules of the game completely.
2. Snap vs. Slide: Sometimes things change smoothly; sometimes they snap violently, and we now know exactly where that "snap" happens.
3. Memory: These quantum systems can have "memories," staying in different states depending on their history.
4. Heat helps: Sometimes, a little bit of warmth is actually the key to unlocking these amazing quantum effects.

This research helps scientists understand how to build better quantum computers and simulate complex materials, using light and cold atoms as the building blocks.

Technical Summary

1. Problem Statement

The paper investigates the superradiant phase transition in a system of spinless fermions confined in a 1D optical lattice and coupled to a cavity mode. While superradiance in Bose-Einstein condensates (BECs) and the Dicke model is well-studied, the behavior of fermionic superradiance remains less explored, particularly regarding:

• The emergence of tricritical points (where second-order and first-order phase transitions meet).
• The role of higher-order Fermi surface nesting (FSN) in strong pumping regimes.
• The existence of multistability and hysteresis in dissipative cavities.
• The evolution of critical phenomena at finite temperatures, specifically the scaling behavior of the critical boundary near absolute zero.

A key open question is how dimensionality (1D vs. 2D) and temperature affect the stability of these phases and the nature of the phase transitions.

2. Methodology

The authors employ a combination of analytical perturbation theory, mean-field approximations, and numerical simulations:

• Model Hamiltonian: They consider a macroscopic number of fermions in a cavity pumped by a strong standing-wave laser. Using adiabatic elimination of the excited atomic states, they derive an effective Hamiltonian describing the coupling between the cavity field and the fermions via photon scattering.
• Landau Theory Expansion: The free energy ($F$) is expanded in terms of the order parameter (cavity field amplitude $\psi$) up to the fourth order:
  $$F = \tilde{\omega}|\psi|^2 + \chi(\psi + \psi^*)^2 + \eta(\psi + \psi^*)^4 + \dots$$
  where $\chi$ (2nd order) and $\eta$ (4th order) are calculated perturbatively.
• Dimensionality Comparison: The study explicitly compares 1D and 2D systems. The authors analyze the infrared divergence properties of the perturbation integrals in both dimensions.
• Finite Temperature Analysis: Analytical expressions for $\chi$ and $\eta$ are derived for finite temperatures ($T > 0$), accounting for the temperature dependence of the chemical potential $\mu$.
• Dissipative Dynamics: For dissipative cavities (finite decay rate $\kappa$), the authors solve the master equation using a numerical iterative method to find stable states and simulate quench dynamics to observe hysteresis.
• Scaling Analysis: The critical scaling laws near the tricritical point are derived and verified numerically.

3. Key Contributions

A. Origin of Tricriticality in 1D Systems

The paper identifies that the tricritical point in 1D fermionic superradiance arises from the interplay between:

1. First-order FSN: Leading to the standard critical boundary ($\chi = 0$).
2. Second-order FSN (Higher-order): Leading to the fourth-order coefficient $\eta$.

• Crucial Finding: In 1D, the second-order FSN contribution causes an infrared divergence in the coefficient $\eta$ at specific filling factors ($k_F \sim 1/2$ and $k_F \sim 1$). This allows $\eta$ to change sign (become negative), creating a tricritical point where the transition changes from second-order to first-order.
• Contrast with 2D: In 2D systems, the additional integration dimension renders the fourth-order coefficient convergent and strictly positive ($\eta > 0$). Consequently, tricritical points are absent in 2D under these conditions.

B. Multistability and Hysteresis

• The authors map out the stable phase diagram for dissipative cavities. They identify a region (beyond the tricritical point) where three stable states coexist: one normal phase (NP) and two superradiant phases (SRP).
• Numerical simulations of quench dynamics (forward and backward pulse sequences) demonstrate hysteresis, confirming that the final state depends on the initial conditions, a hallmark of multistability.

C. Finite-Temperature Phase Diagrams and Scaling

• Tricritical Curves: As temperature increases, the zero-temperature tricritical point extends into a tricritical curve.
• Two Types of Tricritical Points: The study distinguishes between:
  • Quantum-type: Occurs at low temperatures where the Fermi surface is preserved; dominated by 2nd-order FSN effects.
  • Classical-type: Occurs at higher temperatures where the Fermi distribution broadens; the 2nd-order FSN effect is submerged by thermal fluctuations.
• Non-trivial Scaling Law: The authors prove that the critical pumping strength $B_c$ is independent of temperature to the first order ($dB_c/dT|_{T \to 0} = 0$). This leads to a scaling law:
  $$\Delta B_c \sim \Delta T^\nu \quad \text{with} \quad \nu > 1$$
  This indicates a very flat critical boundary near $T=0$, differing from standard critical behaviors.

D. Optimal Temperature for Observation

Contrary to the intuition that quantum effects are best observed at $T=0$, the paper finds that there exists an optimal finite temperature where the minimum coupling strength required to observe superradiance is lowest. This facilitates experimental observation of these critical phenomena.

4. Key Results

• Phase Diagrams:
  • 1D: Shows a rich phase diagram with a tricritical point at $k_F \approx 0.9425$ (at $T=0$). The transition is second-order for $k_F < k_F^t$ and first-order for $k_F > k_F^t$.
  • 2D: Shows only second-order transitions; no tricritical point exists.
• Scaling Exponents:
  • Normal critical point: $\nu = 1/2$.
  • Tricritical point: $\nu = 1/4$ (for the order parameter vs. pumping strength).
• Multistability: Verification of a region supporting 1, 2, or 3 stable states depending on the pumping strength and filling factor.
• Temperature Dependence: The critical boundary remains flat ($\nu > 1$) near $T=0$, and the tricritical point evolves into a curve separating quantum and classical regimes.

5. Significance

• Theoretical Insight: The work provides a rigorous theoretical framework for understanding how dimensionality and higher-order scattering processes (FSN) dictate the nature of phase transitions in fermionic quantum gases. It resolves why tricriticality is a 1D phenomenon in this context.
• Experimental Guidance: By identifying an "optimal temperature" and predicting specific filling factors ($k_F$) where multistability and tricriticality occur, the paper offers concrete targets for experimentalists working with ultracold atoms in optical cavities.
• New Scaling Laws: The discovery of the scaling rate $\nu > 1$ for the critical boundary near zero temperature adds a new dimension to the understanding of quantum criticality in open quantum systems.
• Quantum Simulation: The findings broaden the scope of quantum simulation using cold atoms, demonstrating their capability to simulate complex many-body phenomena like quantum frustration, spin liquids, and non-equilibrium dynamics mediated by long-range cavity interactions.

In summary, this paper establishes that higher-order Fermi surface nesting in 1D systems is the key mechanism driving tricritical phenomena and multistability in fermionic superradiance, a feature absent in 2D, and provides a comprehensive map of these behaviors across temperature and interaction strength.