Multiband Superconductivity in the Exactly Solvable Hatsugai-Kohmoto Model

This paper investigates multiband superconductivity within the exactly solvable Hatsugai-Kohmoto model by classifying symmetry-allowed gap structures in a two-orbital system and computing critical temperatures and order parameters to establish a systematic framework for analyzing the interplay of strong correlations, orbital structure, and pairing symmetry.

The Gist

Imagine a bustling city where electrons are the citizens. In most materials, these citizens move around freely, like people walking through open streets. But in certain special materials called superconductors, these electrons decide to pair up and dance in perfect unison, allowing electricity to flow without any resistance.

This paper explores a specific, highly complex type of superconductor where the "citizens" (electrons) have multiple identities or "jobs" (called orbitals) they can hold, and they are also very "social" (strongly correlated), meaning their behavior is heavily influenced by their neighbors.

Here is a breakdown of the paper's story using simple analogies:

1. The Setting: The "Hatsugai-Kohmoto" City

The authors use a mathematical model called the Hatsugai-Kohmoto (HK) model. Think of this model as a simplified, perfectly organized city map.

• The Special Rule: In this city, every citizen interacts with everyone else instantly, regardless of distance. It's like if you could hear a whisper from the other side of the world instantly.
• Why use it? Because of this weird rule, the city is "exactly solvable." This means the authors can calculate exactly how the citizens behave without needing to make messy approximations. It's a perfect laboratory to test ideas about how strong social pressure (correlations) affects dancing (superconductivity).

2. The Twist: Adding "Orbitals" (Multiple Jobs)

Previous studies looked at electrons with just one "job" (one orbital). This paper upgrades the model to a two-orbital system.

• The Analogy: Imagine the citizens now have two hats they can wear: a "Red Hat" and a "Blue Hat." They can switch between them or wear them in combinations.
• The Challenge: Now, when two electrons decide to dance (pair up), they have to coordinate not just their steps (spin) but also which hats they are wearing (orbitals). This creates a much richer, more complex landscape of possible dances.

3. The Goal: Classifying the Dances (Symmetry)

The first major part of the paper is like a dance instructor cataloging every possible way these two-hatted citizens can pair up while obeying the city's laws (symmetry rules).

• The Laws: The city has a specific shape (a square grid with specific symmetries). The laws say that if you rotate the city or flip it, the dance must look consistent.
• The Result: The authors created a massive "menu" of allowed dances. They found that electrons can pair up in many new ways:
  • Spin Singlet/Triplet: How their internal spins align (like holding hands vs. high-fiving).
  • Orbital Singlet/Triplet: How their "hats" align (both red, both blue, or mixed).
  • They listed specific patterns (like $A_{1g}$, $E_u$, etc.) that act as the "choreography" for these dances.

4. The Experiment: Turning Up the Heat and Pressure

In the second half, the authors simulate what happens when they change the conditions:

• The Interaction Strength ($U$): This is like turning up the volume on the citizens' gossip. When the gossip is low, they dance easily. When it gets very loud (strong correlation), they might stop moving entirely (a "Mott transition," where they get stuck in place).
• The Pairing Strength ($g$): This is how much the citizens want to dance.

What they found:

• The "Mott" Wall: There is a critical point where the gossip becomes so loud that the citizens freeze. The authors found that superconductivity behaves very differently before and after this freezing point.
• Sudden Jumps vs. Smooth Slides:
  • In some dance styles, as the temperature rises, the dancing slows down smoothly until it stops (a normal transition).
  • In other styles, especially when the gossip is very loud (in the "Mott regime"), the system acts strangely. It might be dancing, then suddenly stop, then start dancing again at a different temperature. It's like a light switch that flickers before turning off, rather than a dimmer switch. This is called a first-order phase transition.
• The Sweet Spot: The "best" dancing (highest critical temperature) doesn't happen when the citizens are totally free or totally frozen. It happens at a medium level of gossip. If the interaction is too weak or too strong, the superconductivity dies out.

5. The Takeaway

This paper doesn't invent a new superconductor for your phone or a new medical device. Instead, it provides a theoretical map.

It tells us that when you have electrons with multiple identities (orbitals) that are strongly influenced by each other, the rules for how they pair up become incredibly complex. The authors have written down the "rulebook" for these complex dances and shown that the transition from "dancing" to "frozen" can be abrupt and surprising, depending on how strong the interactions are.

In short: They built a perfect, solvable toy city to understand how complex electron dances work when the electrons are very social and have multiple identities, revealing that the path to superconductivity can be bumpy and full of sudden jumps, not just a smooth slide.

Technical Summary

Problem Statement
Multiband superconductivity presents a complex landscape of pairing states, particularly when strong electron correlations are involved. While the Hatsugai-Kohmoto (HK) model is an exactly solvable framework for studying correlated electrons with momentum-local interactions, its application has largely been restricted to single-band scenarios or specific contexts like the Mott transition. A gap exists in the systematic understanding of how strong correlations, orbital degrees of freedom, and pairing symmetries interplay in a multiband setting. Specifically, there is a need to classify the symmetry-allowed superconducting gap structures in an orbital HK model and to determine how the Mott transition and interaction strength influence the critical temperature and the nature of the superconducting phase transition.

Methodology
The authors extend the exactly solvable HK model to a two-orbital system (the "orbital HK model") defined on a square lattice with $p$-orbitals, possessing $D_{4h}$ point-group symmetry. The model incorporates nearest- and second-nearest-neighbor hopping, resulting in a two-band Bloch Hamiltonian. To introduce superconductivity, a pairing term is added and treated within a mean-field approximation. This approach preserves the momentum-local structure of the HK interaction, allowing the Hamiltonian to be decoupled into independent momentum sectors that can be diagonalized exactly.

The study proceeds in two main stages:

1. Symmetry Classification: The authors perform a comprehensive classification of superconducting gap structures by considering spin, orbital, and momentum degrees of freedom. They derive the constraints imposed by Fermi antisymmetry ($SPO = -1$, where $S$, $P$, and $O$ denote spin, parity, and orbital exchange behaviors) and classify the resulting gap functions according to the irreducible representations of the $D_{4h}$ point group. This involves distinguishing between spin-singlet/triplet and orbital-singlet/triplet sectors.
2. Numerical Analysis: For selected representative pairing channels (specifically $A_{1g}$, $E_u$, $B_{2g}$, and $B_{1g}$), the authors compute the mean-field free energy and the critical temperature ($T_c$) as functions of interaction strength ($U$) and pairing strength ($g$). The calculations are performed at half-filling, with the chemical potential fixed by the normal-state problem. The free energy landscapes are analyzed to identify global minima, local extrema, and the nature of phase transitions (continuous vs. first-order).

Key Contributions

• Systematic Symmetry Classification: The paper provides a complete classification of symmetry-allowed superconducting basis functions for a two-orbital HK model under $D_{4h}$ symmetry. This generalizes previous single-band results to include orbital-singlet and orbital-triplet states, yielding basis functions for all irreducible representations (including $A_{1g}, A_{2g}, B_{1g}, B_{2g}, E_g$ for even parity and $A_{1u}, A_{2u}, B_{1u}, B_{2u}, E_u$ for odd parity).
• Exact Diagonalization in Multiband Context: The authors demonstrate that the mean-field treatment of the orbital HK model retains exact solvability in each momentum sector, enabling the calculation of thermodynamic quantities without relying on approximate many-body techniques like dynamical mean-field theory (DMFT) for the correlation part.
• Free Energy Topology Analysis: The study reveals distinct free-energy topologies depending on the interaction regime. Below the Mott transition, the system exhibits conventional BCS-like continuous transitions. Within the Mott regime, the authors identify first-order phase transitions characterized by coexisting minima and metastable states.
• Channel-Dependent Behavior: The analysis highlights that the critical temperature and the stability of the superconducting state are highly dependent on the specific pairing channel. For instance, the $E_u$ channel displays anomalous free-energy structures (persistent local maxima down to $T=0$) in the metallic regime, which are suppressed in the Mott regime.

Results

• Mott Transition: The model exhibits a Mott transition at a critical interaction strength $U_c = W$ (where $W$ is the non-interacting bandwidth), consistent with the single-band HK model.
• Phase Transition Nature: In the pairing channels $A_{1g}$, $B_{2g}$, and $B_{1g}$, the transition is continuous below $U_c$ but becomes first-order within the Mott regime, where the global minimum at $\Delta \neq 0$ is separated from the normal state by a local maximum.
• Critical Temperature: $T_c$ is maximized at finite interaction strengths, well below the Mott transition point. Increasing the pairing strength $g$ generally enhances $T_c$. However, the behavior varies by channel: the $A_{1g}$ and $B_{2g}$ channels show broad regions of finite $T_c$, whereas the $E_u$ and $B_{1g}$ channels exhibit regions where the global minimum remains at $\Delta = 0$ until a threshold in $g$ is reached.
• Anomalous Behavior: The $E_u$ channel (spin-singlet, orbital-triplet) shows a unique free-energy landscape where a local maximum persists down to zero temperature in the metallic regime, a feature that disappears as the interaction strength increases into the Mott regime.

Significance
The paper establishes a systematic framework for analyzing superconductivity in correlated multiband systems using the exactly solvable orbital HK model. By generalizing symmetry-based approaches to include orbital degrees of freedom, the work provides a minimal yet rigorous setting to explore the interplay between strong correlations and pairing symmetry. The authors claim that while their specific classification is tailored to the $p$-orbital square lattice model, the HK construction is generalizable to other band structures and lattice geometries. Consequently, this work serves as a concrete example for classifying and analyzing superconducting order in orbital HK models, offering a foundation for understanding correlated multiband superconductivity in more complex materials.