Revisiting vestigial order in nematic superconductors: gauge-field mechanisms and model constraints

Through large-scale Monte Carlo simulations of a three-dimensional Ginzburg-Landau model, this study confirms that standard models of nematic superconductors do not support vestigial nematic phases, but demonstrates that such order can emerge under restrictive conditions via gauge-field-mediated intercomponent coupling or strong correlations.

The Gist

The Big Idea: The "Ghost" of Superconductivity

Imagine a dance floor where couples (electrons) are dancing in perfect sync. This is a superconductor. In a normal superconductor, the couples are just holding hands and spinning. But in a special type called a nematic superconductor, the couples also have a specific "facing direction." They might all be facing North, or East, or South. This directionality is called nematic order.

Usually, when the music stops (the material gets too hot), the couples break up, the spinning stops, and the directionality disappears all at once. The dance floor becomes a chaotic mess of people walking randomly.

The Mystery:
Scientists wondered: Is it possible for the couples to break up and stop spinning first, but for the "facing direction" to stay for a little while longer?

• Phase 1: Everyone is dancing in a specific direction (Superconducting + Nematic).
• Phase 2 (The "Vestigial" Phase): The couples break up and stop spinning, but they are still all facing the same direction. It's like a crowd of people who have stopped dancing but are still standing in a perfect line, all looking North.
• Phase 3: Eventually, even the direction is lost, and everyone walks randomly.

This "Phase 2" is called a vestigial nematic phase. It's a "ghost" of the superconducting order that lingers after the main event is over.

The Problem: The Math Said "No"

Scientists had been arguing about whether this "ghost phase" could actually exist in real materials (like the famous $Bi_2Se_3$ crystals).

• Team A said: "Yes! The math shows the direction should stick around even after the superconductivity breaks."
• Team B (How and Yip) said: "No. If you look closely at the standard math models, the direction disappears at the exact same time the superconductivity does. There is no 'in-between' phase."

The Investigation: A Giant Digital Simulation

The authors of this paper (Maccari, Babaev, and Carlström) decided to settle the debate by building a massive digital simulation. Instead of doing simple math on paper, they used a supercomputer to simulate millions of electrons interacting in a 3D grid, step-by-step, as they heated up the material.

Think of it like a video game where you control the temperature and watch to see if the "directional line" survives the breakup.

The Findings: Two Different Stories

Their results were a mix of "You were right" and "But wait, there's a catch."

1. The Standard Model (No Magic Fields)

When they simulated the material using the standard, most common rules (ignoring complex magnetic interactions), Team B was right.

• The Result: As soon as the superconductivity broke, the directional order vanished instantly.
• The Analogy: Imagine a line of soldiers. If you tell them to drop their weapons (stop superconducting), they immediately scatter and lose their formation. There is no moment where they are unarmed but still standing in a perfect line.
• Conclusion: In the simplest models, the "ghost phase" does not exist.

2. The Twist: Introducing the "Gauge Field"

However, the authors didn't stop there. They asked: "What if we add a strong magnetic influence (a gauge field) or strong interactions between the electrons?"

• The Result: Suddenly, the "ghost phase" did appear, but only under very specific, strict conditions.
• The Analogy: Imagine the soldiers are in a very strong wind (the gauge field). If the wind is just a breeze, they scatter when they drop their weapons. But if the wind is a hurricane, it might push them into a line and hold them there even after they drop their weapons.
• The Catch: The wind (the magnetic coupling) has to be extremely strong. In the simulation, the material had to be pushed to a point where the magnetic forces were huge (represented by a value $e > 3.5$).

Why Does This Matter?

This paper tells us two important things:

1. Don't expect to find this easily: If you look at standard nematic superconductors (like doped $Bi_2Se_3$) without doing anything special, you probably won't see this "vestigial" phase. The math and the simulation agree that it likely doesn't happen naturally in those simple models.
2. It's possible, but you have to force it: To see this strange "ghost" state where direction remains after superconductivity dies, you might need to apply a very strong external magnetic field. This field would melt the superconducting "dance" but leave the "directional line" intact.

The Bottom Line

The authors successfully simulated the behavior of these exotic materials. They confirmed that in the "real world" of standard models, the vestigial phase is likely a myth. However, they also proved that it can exist if you create a very specific, high-pressure environment (strong magnetic coupling).

It's like proving that while a house of cards will always fall when you pull the bottom card, if you glue the cards together with super-strong glue (the gauge field), you might be able to pull the bottom card and still have a standing structure for a moment. But that glue is hard to find in nature.

Technical Summary

1. Problem Statement

The paper addresses the existence and stability of vestigial nematic order in three-dimensional (3D) nematic superconductors.

• Context: In certain superconductors (e.g., $Bi_2Se_3$-based materials), the superconducting order parameter breaks both $U(1)$ gauge symmetry and a discrete rotational symmetry ($Z_3$).
• The Vestigial Phase Hypothesis: It has been proposed that above the superconducting critical temperature ($T_c$), the $U(1)$ symmetry might be restored while the $Z_3$ nematic symmetry remains broken. This would result in a "vestigial" phase characterized by a composite order parameter (formed by four fermionic fields, e.g., charge-4e or nematic order) without long-range superconductivity.
• The Conflict: Recent analytical work by How and Yip (2023) concluded that standard Ginzburg-Landau (GL) models for these systems do not support such vestigial phases; instead, the $U(1)$ and $Z_3$ transitions occur simultaneously.
• The Gap: Analytical approaches often rely on mean-field approximations or Gaussian fluctuations, which may fail to capture the complex interplay of topological defects (vortices, domain walls) and strong fluctuations in 3D. The authors aim to resolve this discrepancy using large-scale numerical simulations.

2. Methodology

The authors employ large-scale Monte Carlo (MC) simulations on a discretized 3D lattice to study a two-component Ginzburg-Landau model.

• Model Hamiltonian:
  The system is described by a two-component order parameter $\vec{\Delta} = (\Delta_1, \Delta_2)$. The free energy includes:

  • Gradient terms with a gauge field $\vec{A}$ (coupling constant $e$).
  • A potential $V_{nem}$ containing quadratic, quartic, and higher-order terms.
  • Crucially, the potential includes a $Z_3$-symmetry breaking term (coefficient $\lambda$) and a term controlling the relative phase/amplitude (coefficient $K = 2\beta_2$).
  • The model allows for the study of both the neutral limit ($e=0$, extreme Type-II) and the charged limit ($e \neq 0$).

• Simulation Techniques:

  • Algorithm: Metropolis-Hastings sweeps combined with parallel tempering (replica exchange) to overcome energy barriers and ensure ergodicity.
  • System Size: Simulations performed on lattices up to $L=32$ (with extrapolation to the thermodynamic limit).
  • Observables:
    • Superconducting Transition ($U(1)$): Measured via the helicity modulus sum ($\Upsilon_+$) for $e=0$ and dual stiffness ($\rho$) for $e \neq 0$ (to account for gauge invariance).
    • Nematic Transition ($Z_3$): Measured via the nematic order parameter ($O_{Nem}$) and the Binder cumulant ($U_{nem}$).
    • Thermodynamics: Specific heat ($C_v$) and energy distributions to detect phase transition order (first vs. second).
  • Finite-Size Scaling (FSS): Critical temperatures ($T_c$) are extracted by analyzing the crossing points of observables for different system sizes $L$, using scaling ansatz appropriate for 3D XY (for $U(1)$) and 3-state Potts (for $Z_3$) universality classes.

3. Key Contributions

1. Validation of Analytical Results: The authors provide numerical confirmation that standard, widely used models of nematic superconductors (without strong gauge coupling) do not exhibit a vestigial nematic phase. The $U(1)$ and $Z_3$ transitions coincide.
2. Identification of a Gauge-Field Mechanism: They demonstrate that a vestigial phase can emerge, but only under a specific mechanism: intercomponent coupling mediated by the electromagnetic gauge field.
3. Dimensionality Constraints: The paper clarifies why this mechanism is difficult to realize in 3D compared to 2D, highlighting the role of vortex line tension and the energy cost of topological defects.
4. Proposal for Experimental Realization: The authors propose that realizing this phase in materials like $Bi_2Se_3$ may require external suppression of superconductivity (e.g., via magnetic fields) to melt a vortex lattice, thereby decoupling the transitions.

4. Key Results

A. Zero Gauge Coupling Limit ($e = 0$)

• Scenario: Corresponds to a neutral system or extreme Type-II superconductor.
• Finding: For all investigated values of the nematic coupling $\lambda$, the system undergoes a single phase transition.
• Evidence:
  • The specific heat ($C_v$) shows a single peak.
  • The crossing points of the nematic Binder cumulant ($U_{nem}$) and the helicity modulus ($\Upsilon_+$) converge to the same critical temperature ($\beta_c$) in the thermodynamic limit.
  • Conclusion: No vestigial phase exists in the standard neutral model. This aligns with the analytical findings of How and Yip.

B. Finite Gauge Coupling ($e \neq 0$)

• Scenario: The electromagnetic field couples to the supercurrents.
• Finding: A vestigial nematic phase emerges, but only for very strong gauge coupling ($e > 3.5$ in their dimensionless units).
• Evidence:
  • Splitting of Transitions: For $e > 3.5$, the specific heat exhibits two distinct peaks, indicating two separate phase transitions.
  • Order of Transitions: The $Z_3$ nematic transition occurs at a higher temperature ($T_{Z3} > T_{U(1)}$) than the superconducting transition.
  • Intermediate Phase: Between these temperatures, the system is in a state where $Z_3$ symmetry is broken (nematic order exists) but $U(1)$ symmetry is restored (no superconductivity). This is the vestigial nematic phase.
• Mechanism: The gauge field reduces the energy cost of skyrmionic vortices (which carry integer flux). Their proliferation destroys superconductivity ($U(1)$) while leaving the nematic order ($Z_3$) intact, provided the coupling is strong enough to decouple the defect proliferation.

C. Nature of the Transition

• As the system approaches the bicritical point where the two transitions split, the transition acquires a strong first-order character.
• This is evidenced by a bimodal energy distribution and significant finite-size effects near the splitting point.

5. Significance and Implications

• Resolution of Controversy: The paper resolves the conflict between earlier analytical predictions and experimental claims of vestigial order. It confirms that the standard GL models do not naturally produce vestigial phases in 3D.
• Mechanism for Composite Order: It identifies a specific, albeit restrictive, mechanism (strong gauge coupling) that can stabilize composite orders. This suggests that vestigial order is not a generic feature of all nematic superconductors but requires specific conditions.
• Experimental Guidance:
  • The requirement for "very strong" gauge coupling in the model suggests that in real materials (like $Bi_2Se_3$), the effect might be too weak to observe spontaneously.
  • Proposed Solution: The authors suggest that applying an external magnetic field could induce the necessary conditions. A magnetic field creates a vortex lattice; melting this lattice (via thermal fluctuations) could restore gauge symmetry while preserving nematic order, effectively "catalyzing" the vestigial phase.
• Theoretical Insight: The work highlights the critical importance of topological defects (skyrmions, fractional vortices) and their interactions in determining the phase diagram of multicomponent superconductors, which mean-field theories often miss.

In summary, the paper establishes that while vestigial nematic order is theoretically possible in 3D nematic superconductors, it is not a robust feature of standard models. It requires specific, strong coupling mechanisms (likely induced externally) to decouple the nematic and superconducting transitions.