Competition and coexistence of superconductivity and nematic order in a two-dimensional electron gas with quadrupolar interactions

This paper investigates a two-dimensional electron gas model with competing pairing and quadrupolar interactions, revealing that while nematic order strongly competes with $d$-wave superconductivity to induce a first-order phase transition, it can coexist with $s$-wave pairing to form an anisotropic superconducting state, with finite-temperature effects further enabling the simultaneous emergence of $s$-wave, $d$-wave, and nematic orders.

The Gist

Imagine a crowded dance floor filled with electrons. Usually, these electrons move around chaotically, bumping into each other but mostly ignoring one another. However, under certain conditions, they can start dancing in perfect sync, forming a superconductor (where electricity flows with zero resistance) or organizing themselves into a nematic phase (where they all face the same direction, like a school of fish turning together, but without forming a rigid grid).

This paper is a theoretical study that asks a big question: What happens when you try to force these two different dance styles to happen at the same time?

Here is a simple breakdown of their findings using everyday analogies.

1. The Two Types of Dancers

The researchers looked at a 2D "dance floor" (a thin layer of material) where electrons have two main ways to interact:

• The Pairing Dance (Superconductivity): Electrons pair up and move together. They can do this in two styles:
  • S-wave: A simple, round, uniform hug. Everyone pairs up the same way, no matter which direction they are facing.
  • D-wave: A more complex, directional dance. The "hug" changes depending on the angle. It has a specific shape (like a four-leaf clover).
• The Nematic Dance (Nematic Order): Here, the electrons don't pair up, but they all decide to "lean" in a specific direction. Imagine a crowd of people who suddenly all tilt their heads to the right. The crowd is still fluid, but it has lost its rotational symmetry (it's no longer the same in all directions).

2. The Great Conflict: When Styles Clash

The paper explores what happens when the electrons are forced to choose between these dances.

Scenario A: The D-Wave vs. The Nematic (The Rivalry)

• The Analogy: Imagine a dance-off between a group doing the "Four-Leaf Clover" dance (D-wave) and a group trying to "tilt their heads" (Nematic).
• The Problem: Both of these styles rely on the same specific directionality. They are fighting for the same "real estate" on the dance floor.
• The Result: They cannot coexist peacefully. It's like trying to wear a left shoe and a right shoe on the same foot. The paper finds that when the "tilting" force gets too strong, it violently kicks out the "Four-Leaf Clover" dance. The system snaps from one state to the other in a sudden, first-order transition (like a light switch flipping). One wins, the other loses.

Scenario B: The S-Wave vs. The Nematic (The Unlikely Roommates)

• The Analogy: Now, imagine the "Round Hug" dance (S-wave) trying to coexist with the "Head-Tilting" crowd.
• The Difference: The "Round Hug" doesn't care about direction; it's uniform. The "Head-Tilting" crowd cares about direction. Because they are so different, they don't fight as hard.
• The Result: They can actually live together! The electrons can form their uniform pairs while the whole crowd leans in one direction. The result is a hybrid state: a superconductor that flows perfectly, but on a distorted, anisotropic (stretched) dance floor. It's like a perfectly synchronized marching band that is marching on a slanted floor.

3. The Temperature Factor

The researchers also looked at what happens when the room gets hotter (adding thermal energy).

• The "Melting" Effect: Heat usually breaks these delicate dances. However, the paper found something surprising. If you have a mix of interactions, heating the system doesn't just kill the superconductivity immediately.
• The Triple Threat: In some cases, the heat and the interactions create a "sweet spot" where all three things happen at once: S-wave pairing, D-wave pairing, and Nematic tilting. It's a chaotic but stable party where three different dance styles are happening simultaneously before the music (the order) eventually stops as it gets too hot.

4. Why Does This Matter?

You might wonder, "Why do we care about electrons tilting their heads?"

• Real-World Connection: This isn't just theory. Materials like iron-based superconductors and cuprates (the stuff used in high-temperature superconductors) show signs of this exact behavior.
• The Takeaway: Understanding how these "dances" compete or cooperate helps scientists design better materials. If we can figure out how to make the "hybrid roommates" (S-wave + Nematic) stable, we might be able to create superconductors that work at higher temperatures or carry more current, which is the holy grail of energy technology.

Summary

Think of this paper as a study of social dynamics in a crowd of electrons:

1. If two groups want to do the same directional dance, they fight, and one kicks the other out.
2. If one group does a directionless dance and the other does a directional dance, they can share the floor.
3. By tweaking the "music" (temperature and interaction strength), you can create complex, multi-layered states where different orders coexist.

The authors built a mathematical "dance floor" to simulate these interactions, proving that the symmetry (the shape of the dance) is the most important factor in deciding whether the electrons will fight or make peace.

Technical Summary

1. Problem Statement

The paper addresses the complex interplay between superconductivity (SC) and electronic nematic order in strongly correlated electron systems. Nematic order is characterized by the spontaneous breaking of rotational symmetry while preserving translational symmetry, a phenomenon observed in materials like iron-based superconductors, cuprates, and quantum Hall systems.

The central theoretical challenge is understanding how these two orders compete or coexist, particularly when they share similar symmetry properties. While previous studies have explored these phases separately or in specific lattice contexts, this work aims to construct a minimal, unified framework for a two-dimensional (2D) electron gas with continuous rotational symmetry. It specifically investigates the competition between:

• s-wave and d-wave superconducting channels.
• Nematic order driven by quadrupolar forward-scattering interactions.

The authors seek to determine how interaction strengths and temperature dictate the phase diagram, specifically looking for regimes of direct competition (mutual exclusion) versus coexistence.

2. Methodology

The authors employ a mean-field (MF) theoretical approach to analyze a 2D electron gas model.

• Hamiltonian: The system is described by a Hamiltonian containing kinetic energy, a BCS-type pairing interaction (involving both s-wave and d-wave channels), and a particle-hole quadrupolar forward-scattering interaction.
  • The pairing interaction $V^{SC}_{kk'}$ includes coupling constants $V_s$ (s-wave) and $V_d$ (d-wave).
  • The nematic interaction $V^N_{kk'}$ is governed by a quadrupolar coupling constant $f_2$, which drives a Pomeranchuk instability.
• Symmetry Considerations: The model assumes continuous rotational symmetry (perturbative lattice effects). The geometric factor $f_{kk'} = \cos(2[\theta_k - \theta_{k'}])$ ensures the nematic interaction shares the same symmetry as the d-wave superconducting order parameter.
• Mean-Field Decoupling: The interaction terms are decoupled into order parameters:
  • $\Delta_s$: s-wave superconducting gap.
  • $\Delta^c_d, \Delta^s_d$: Two components of the d-wave superconducting gap.
  • $\Delta^c_n, \Delta^s_n$: Two components of the nematic order parameter.
• Free Energy Calculation: The authors derive the mean-field free energy density $F$ by diagonalizing the Hamiltonian using the Bogoliubov–Valatin–de Gennes transformation. The quasiparticle dispersion $E_k(\theta)$ incorporates both the superconducting gap and the anisotropic chemical potential induced by nematicity.
• Numerical Solution: The phase diagrams are obtained by minimizing the free energy density with respect to the five order parameters. This leads to a set of five coupled self-consistent equations solved numerically for various interaction strengths ($V_s, V_d, f_2$) and temperatures ($T$).

3. Key Contributions

• Unified Framework: The paper provides a minimal model that treats s-wave, d-wave, and nematic orders on equal footing within a continuous rotational symmetry framework.
• Symmetry-Driven Topology: It explicitly demonstrates how the relative symmetry of the competing order parameters dictates the topology of the phase diagram. Specifically, it contrasts the behavior of nematic order against d-wave (same symmetry) versus s-wave (different symmetry) superconductivity.
• Finite-Temperature Effects: The study reveals that quadrupolar interactions can induce additional superconducting components at finite temperatures, leading to complex coexistence regimes not present at $T=0$.

4. Key Results

A. Zero-Temperature ($T=0$) Phase Diagrams

1. s-wave vs. d-wave (No Nematicity): When $f_2 = 0$, s-wave and d-wave superconductivity are mutually exclusive. They are separated by a first-order phase transition, with the s-wave phase generally favored energetically near the boundary.
2. d-wave SC vs. Nematicity:
   • Since both d-wave SC and nematic order break rotational symmetry in the same way (sharing the same symmetry channel), they strongly compete.
   • The transition between the d-wave superconducting phase and the nematic phase is a direct first-order transition.
   • Strong attractive quadrupolar interactions ($f_2 < 0$) suppress d-wave superconductivity discontinuously.
3. s-wave SC vs. Nematicity:
   • Because s-wave SC preserves rotational symmetry while nematic order breaks it, they belong to different symmetry classes.
   • Coexistence Phase: For sufficiently strong s-wave coupling ($|V_s|$), a coexistence phase emerges. In this phase, a uniform s-wave superconducting gap exists on an anisotropic (nematic) Fermi surface.
   • The transition from the disordered state to the coexistence phase can be continuous (second-order) or first-order depending on the path in parameter space.

B. Finite-Temperature Effects

• Induced Components: At finite temperatures, quadrupolar interactions can stabilize a metastable d-wave component even when the ground state is s-wave.
• Triple Coexistence: In regimes with strong quadrupolar interactions and specific coupling strengths, the system enters a phase where s-wave, d-wave, and nematic orders coexist simultaneously.
• Critical Points:
  • Tetra-critical point: Observed in the s-wave dominant regime (Fig. 4), where three second-order lines meet a first-order line.
  • Tri-critical point: Observed in the d-wave dominant regime (Fig. 5), where two second-order lines meet a first-order line.
• Thermal Suppression: As temperature increases, superconductivity is suppressed first, leaving a purely nematic phase, which eventually undergoes a continuous transition to a disordered state.

5. Significance and Implications

• Mechanism for Intertwined Orders: The results provide a theoretical mechanism for how nematicity and superconductivity can coexist in materials like iron-based superconductors and cuprates. The finding that s-wave and nematic orders can coexist on an anisotropic Fermi surface offers a potential explanation for experimental observations of nematicity coexisting with superconductivity.
• Role of Symmetry: The work highlights that the nature of the phase transition (first-order vs. coexistence) is fundamentally determined by whether the competing orders share the same symmetry representation.
• Beyond Mean-Field: While the study uses a mean-field approximation, the authors acknowledge the Mermin-Wagner theorem (which forbids long-range order in 2D at finite $T$ due to Goldstone modes). However, they argue that in realistic systems, weak coupling to the underlying lattice explicitly breaks continuous symmetry, opening a gap in the Goldstone mode spectrum and stabilizing the mean-field phase structure qualitatively.
• Future Directions: The paper sets the stage for future investigations into fluctuation effects beyond mean-field, the role of specific lattice structures, and the application of this framework to more realistic band structures in correlated materials.

In summary, this paper establishes that quadrupolar interactions act as a crucial tuning parameter that can drive first-order transitions between competing orders or stabilize complex coexistence phases, depending heavily on the symmetry matching between the superconducting channel and the nematic order.