A tractable framework for phase transitions in phase-fluctuating disordered 2D superconductors: applications to bilayer MoS$_2$ and disordered InO$_x$ thin films

This paper develops a self-consistent microscopic thermodynamic framework that unifies fermionic quasiparticles, Coulomb-regularized phase fluctuations, and BKT vortices to explain how disorder and carrier density drive the separation between the superconducting gap and transition temperature in 2D systems, successfully reproducing experimental results for bilayer MoS$_2$ and disordered InO$_x$ films.

The Gist

Imagine you are trying to build a perfect, synchronized dance troupe. Every dancer (an electron) must move in perfect unison, holding hands to form a single, giant, flowing entity. In the world of physics, this synchronized state is called superconductivity, where electricity flows with zero resistance.

For decades, scientists believed that if you just had enough dancers and they were paired up correctly, the dance would happen automatically, regardless of how messy the dance floor was. This was the "Mean-Field Theory" view: if the pairs exist, the dance exists.

However, in the 2D world (like a single sheet of paper or a thin film), things get chaotic. This new paper by Yang and Chen introduces a new way to understand why these 2D dances sometimes fail, even when the pairs are there. They built a "tractable framework"—essentially a new rulebook—that accounts for three specific types of chaos that happen in 2D superconductors.

Here is the breakdown using simple analogies:

1. The Three Types of Chaos

In a 2D superconductor, the dancers aren't just standing still; they are constantly wobbling. The paper treats three types of wobbling as equally important:

• The "Smooth Wobble" (NG Modes): Imagine the dancers trying to sway gently in a long line. In a 3D room, long-range forces (like Coulomb interactions) act like a heavy ceiling that stops this swaying from getting too wild. But in 2D, the ceiling is gone. The paper shows that while long-range forces still help calm these waves down, they don't stop them completely.
• The "Vortex Storm" (BKT Fluctuations): This is the big one. Imagine a tornado forming in the middle of the dance floor. In 2D, these tiny tornadoes (called vortices) can pop up easily. If too many tornadoes form, they rip the synchronized dance apart, even if the dancers are still holding hands. This is the Berezinskii–Kosterlitz–Thouless (BKT) effect.
• The "Messy Floor" (Disorder): Real dance floors aren't perfect; they have bumps, holes, and sticky spots (impurities). In 3D, the dancers can just step over these. In 2D, a messy floor makes the dancers stumble, which weakens their ability to stay in sync.

2. The Big Discovery: Two Different Temperatures

The most exciting finding of this paper is that there are two different temperatures where things break down, not just one.

• Temperature A ($T^*$): The "Pairing" Temperature. This is when the dancers actually find a partner and hold hands. They are paired up, but they are still wobbling and spinning out of sync. They are a "mush" of pairs, but not a synchronized dance yet.
• Temperature B ($T_c$): The "Superconducting" Temperature. This is the lower temperature where the dancers finally stop wobbling, the tornadoes die down, and the whole troupe moves as one perfect unit.

The Analogy: Think of it like a crowd at a concert.

• At $T^*$, everyone has found a friend (pairing), but the crowd is just a chaotic jumble of people talking to their friends.
• At $T_c$, the music starts, and suddenly, everyone starts clapping in perfect unison (superconductivity).
• The "gap" between $T^*$ and $T_c$ is a Pseudogap. It's a weird state where the pairs exist, but the magic of zero-resistance electricity hasn't started yet because the "tornadoes" (vortices) are still ruining the rhythm.

3. Why Disorder and Density Matter

Old theories said that if you added more dancers (higher density) or cleaned up the floor (less disorder), the dance would be perfect. This paper says: Not so fast.

• Low Density: If you have fewer dancers, it's harder to keep the line straight. The "tornadoes" get stronger, pushing $T_c$ way down, even if the dancers are still holding hands.
• Messy Floor (Disorder): If the floor is bumpy, the dancers stumble. This makes the "tornadoes" worse and the "smooth wobbles" stronger. This actually shrinks the gap between the dancers (the superconducting gap) and lowers the temperature where the dance works.

4. Testing the Theory

The authors didn't just write math; they tested their new rulebook on two real-world materials:

1. Bilayer MoS2 (Molybdenum Disulfide): A thin, 2D material that can be tuned like a radio dial (by adding electricity) to change how many dancers are on the floor.
2. Disordered InOx (Indium Oxide): A messy, glass-like film where the floor is very bumpy.

The Result: Their new framework predicted exactly what experiments saw. It explained why the superconducting temperature drops as the material gets messier or emptier, and it accurately calculated the "Pseudogap" region where pairs exist but don't conduct electricity yet.

The Takeaway

This paper is like upgrading the map for a 2D superconductor. Previously, scientists thought the map was a straight line: "If you have pairs, you have superconductivity."

Yang and Chen showed the map is actually a two-step journey. First, you find your partner (Pairing). Then, you have to survive the storm of wobbling and tornadoes to finally start dancing in sync (Superconductivity).

This new framework is a powerful tool because it allows scientists to predict exactly how much "mess" a 2D superconductor can handle before the dance falls apart, which is crucial for building future quantum computers and ultra-efficient electronics.

Technical Summary

Here is a detailed technical summary of the paper "A tractable framework for phase transitions in phase-fluctuating disordered 2D superconductors: applications to bilayer MoS2 and disordered InO$_x$ thin films" by F. Yang and L. Q. Chen.

1. Problem Statement

The paper addresses a critical gap in the theoretical understanding of superconductivity in two-dimensional (2D) disordered systems.

• Limitations of Mean-Field Theory: Conventional BCS mean-field theory, protected by the Anderson theorem, predicts that non-magnetic disorder does not affect the superconducting gap ($\Delta$) or the transition temperature ($T_c$). It also assumes pairing and phase coherence occur simultaneously.
• The 2D Reality: In 2D systems, long-range order is forbidden by the Hohenberg–Mermin–Wagner–Coleman (HMWC) theorem due to gapless Nambu–Goldstone (NG) phase fluctuations. Furthermore, topological Berezinskii–Kosterlitz–Thouless (BKT) vortex fluctuations drive a transition where global phase coherence is lost at a temperature $T_c$ significantly lower than the pairing (gap-closing) temperature $T^*$.
• The Missing Link: Existing theories often treat fermionic quasiparticles, bosonic phase fluctuations (NG modes), topological defects (BKT vortices), and disorder separately or in limiting cases. There is no unified, self-consistent microscopic framework that treats all these degrees of freedom on equal footing to explain the density- and disorder-dependent suppression of superconductivity observed in experiments (e.g., in MoS$_2$ and InO$_x$).

2. Methodology

The authors develop a fully microscopic, self-consistent thermodynamic framework derived from a path-integral approach starting with a microscopic $s$-wave superconducting Hamiltonian.

• Unified Degrees of Freedom: The framework simultaneously incorporates:
  1. Fermionic Bogoliubov quasiparticles.
  2. Bosonic NG phase fluctuations (long-wavelength, smooth deformations).
  3. Topological BKT vortex–antivortex fluctuations.
  4. Long-range Coulomb interactions.
  5. Disorder (via elastic scattering).
• Self-Consistent Equations:
  • Gap Equation: The superconducting gap equation is modified to include the Doppler shift caused by the superfluid momentum ($\mathbf{p}_s$). The NG phase fluctuations enter as a renormalization of the gap via the mean-squared superfluid momentum $\langle p_{s,\parallel}^2 \rangle$, which includes both thermal ($S_{th}$) and zero-point ($S_{zo}$) contributions.
  • Superfluid Stiffness: The bare superfluid density ($n_s$) is calculated self-consistently, incorporating disorder effects via a Tinkham-like interpolation factor $(1 + \xi/l)^{-1}$, where $\xi$ is the coherence length and $l$ is the mean free path.
  • BKT Renormalization: The bare $n_s$ serves as the input for standard BKT Renormalization Group (RG) flow equations to determine the fully renormalized superfluid density ($\bar{n}_s$) and the critical temperature $T_c$.
• Coulomb Regularization: The framework explicitly accounts for the 2D Anderson-Higgs mechanism, where long-range Coulomb interactions convert the gapless NG dispersion from linear ($\omega \propto q$) to plasmonic ($\omega \propto \sqrt{q}$). This eliminates infrared divergences and suppresses thermal NG fluctuations at $T > 0$.
• Disorder Treatment: While the gap equation remains robust against non-magnetic disorder (consistent with Anderson's theorem), the superfluid stiffness is strongly suppressed by impurity scattering, which amplifies both quantum and thermal phase fluctuations.

3. Key Contributions

• Beyond Mean-Field: The theory demonstrates that phase fluctuations, driven by reduced carrier density and increased disorder, generate a density- and disorder-dependent zero-point gap $\Delta(0)$ and transition temperatures, violating the Anderson theorem's prediction of disorder independence for the gap.
• Separation of Scales: It naturally explains the emergence of a pseudogap regime ($T_c < T < T^*$) where Cooper pairs exist but lack global phase coherence. The width of this regime ($T^* - T_c$) is shown to be controlled by the competition between NG quantum fluctuations (affecting $\Delta$) and BKT fluctuations (affecting $T_c$).
• Quantitative Predictive Power: The framework requires only a minimal set of physical parameters (effective mass, pairing strength, and scattering time) to quantitatively reproduce complex experimental trends without ad-hoc fitting of transition temperatures.
• Resolution of HMWC Paradox: It clarifies how 2D superconductivity survives the HMWC theorem: Coulomb interactions suppress thermal NG fluctuations, while the system remains stable against quantum fluctuations only until disorder/density reduction drives the superfluid stiffness to zero.

4. Results and Applications

The framework was applied to two distinct experimental systems, showing excellent quantitative agreement:

A. Gate-Tunable Bilayer MoS$_2$

• System: A multi-valley system where superconductivity arises primarily from the $Q$-valley.
• Findings:
  • The theory reproduces the carrier-density dependence of $T_c$, $T^*$, and $\Delta(0)$.
  • As carrier density decreases, the superfluid stiffness drops, enhancing NG quantum fluctuations. This leads to a significant reduction in the zero-point gap $\Delta(0)$ and a widening of the pseudogap window ($T^* - T_c$).
  • The model correctly predicts the critical carrier density below which superconductivity vanishes due to the depopulation of the $Q$-valley.
  • Calculated critical current $J_c(T)$ matches experimental data without additional fitting parameters.

B. Disordered InO$_x$ Thin Films

• System: Amorphous films where disorder is tuned via oxygen stoichiometry, driving a superconductor-to-insulator transition (SIT).
• Findings:
  • Weak Disorder: As disorder increases, superfluid stiffness is suppressed, leading to a reduction in $T^*$ (due to NG quantum fluctuations) and a rapid separation between $T^*$ and $T_c$ (due to BKT fluctuations).
  • Strong Disorder: In the regime approaching the SIT ($R > 15$ k$\Omega$/sq), the theory predicts a collapse of $T_c$ and $\Theta$ while the local pairing gap $\Delta$ persists, consistent with the formation of a Cooper-pair glass.
  • The model quantitatively reproduces the experimental hierarchy reversal where the pairing scale exceeds the phase coherence scale ($\Theta < \Delta$) in the disordered limit.

5. Significance

• Unified Theory: This work provides the first tractable, self-consistent microscopic framework that unifies fermionic pairing, bosonic phase fluctuations, topological defects, and disorder.
• Experimental Relevance: It resolves long-standing discrepancies between mean-field predictions and experimental observations in 2D materials, specifically explaining why $T_c$ is suppressed by disorder and why a pseudogap phase exists.
• Design Tool: The framework offers a practical tool for predicting the phase diagrams of emerging 2D superconductors (e.g., twisted graphene, transition metal dichalcogenides) and understanding the mechanisms driving the superconductor-insulator transition.
• Limitations and Future Work: The authors note that the current diffusive framework does not capture the physics of the fully localized Cooper-pair glass phase (Anderson localization), which requires a future extension to include spatial inhomogeneity and localization effects explicitly.

In summary, the paper establishes that phase fluctuations are the dominant mechanism governing the stability and critical behavior of 2D disordered superconductors, and it provides a rigorous, quantitative method to model these effects across the entire phase diagram.