Upper critical in-plane magnetic field in quasi-2D layered superconductors

This paper presents an analytical framework to model the upper critical in-plane magnetic field in quasi-2D layered superconductors, which, when applied to recent bilayer graphene experiments, reveals a discrepancy in spin-orbit parameters that is likely resolved by an enhanced Landé g-factor.

The Gist

The Big Picture: Superconductors vs. Magnetic Fields

Imagine a superconductor as a super-highway where electrons (the cars) can drive without any friction or traffic jams. This is amazing because it allows electricity to flow with zero energy loss.

However, there is a villain in this story: Magnetic Fields. Usually, if you bring a strong magnet near a superconductor, it acts like a "force field" that pushes the cars off the highway, destroying the frictionless flow. There is a specific limit, called the Pauli Limit, which is the maximum magnetic strength a superconductor can usually handle before it quits.

The Mystery: Recently, scientists discovered some special "quasi-2D" materials (like stacked layers of graphene) that are incredibly tough. They can withstand magnetic fields far stronger than the Pauli Limit should allow. In fact, in some cases, the magnetic field seems to help the superconductivity rather than kill it. This is like a car that drives faster when a hurricane hits it.

What This Paper Does: The "Rulebook" for Tough Superconductors

The authors of this paper (Ma, Chichinadze, and Lewandowski) wanted to understand why these materials are so tough. They created a new mathematical rulebook (a framework) to predict how these materials behave under magnetic stress.

Think of their rulebook as a video game physics engine.

• The Players: Electrons trying to pair up (like dance partners) to form the superconducting state.
• The Obstacles: Magnetic fields trying to rip the partners apart.
• The Secret Weapons: Two types of "spin-orbit coupling" (SOC).
  • Ising SOC: Imagine this as a magnetic seatbelt. It locks the electrons' "spin" (their internal compass) in a specific direction, making it very hard for the external magnetic field to flip them and break the pair.
  • Rashba SOC: Imagine this as a wobbly dance floor. It makes the electrons' spins wobble and change direction depending on how fast they are moving. This usually makes it easier for the magnetic field to break the pairs.

The Main Discovery: The "Magic Seatbelt"

The researchers applied their rulebook to real-world experiments involving Bernal Bilayer Graphene (a specific type of stacked carbon) wrapped in a material called WSe2.

Here is what they found:

1. The Seatbelt Works: The data showed that the "Ising seatbelt" is the hero here. It is very strong and explains why the material can handle such huge magnetic fields. The "wobbly dance floor" (Rashba) was surprisingly weak and didn't play a big role.
2. The Glitch in the Matrix (The G-Factor Problem): When they tried to fit their math to the real data, something weird happened. To make the numbers match, they had to assume the electrons were acting as if they had a super-powerful magnetic personality (a "g-factor" much larger than normal).
   • Analogy: Imagine you are trying to predict how heavy a suitcase is based on how hard it is to lift. You do the math, and the result says the suitcase weighs 500 lbs. But you know the suitcase is made of light plastic. The only way your math works is if you assume the suitcase has been secretly filled with lead.
   • The Paper's Conclusion: The electrons in these graphene layers seem to be "heavier" magnetically than they should be. The authors suggest this might be because the electrons are interacting with each other in a way that amplifies their magnetic sensitivity, or perhaps the way we measure the "start temperature" of the superconductivity is slightly off.

Why This Matters

This paper is like a decoder ring for the future of quantum technology.

• For Scientists: It gives them a reliable tool to figure out what kind of "dance" (pairing mechanism) the electrons are doing. Are they dancing in a standard way (singlet) or a weird, twisted way (triplet)? The shape of the curve in their graphs tells the story.
• For the Future: If we can understand how to make materials that resist magnetic fields this well, we can build:
  • Better MRI machines that are stronger and cheaper.
  • Quantum computers that don't crash when exposed to magnetic noise.
  • Lossless power grids that can operate in environments where magnets are everywhere.

Summary in One Sentence

The authors built a new mathematical model to explain why certain ultra-thin graphene materials are super-resistant to magnetic fields, discovering that while a specific "locking" mechanism (Ising SOC) is the main hero, the electrons are behaving with a mysterious, amplified magnetic strength that we still need to fully understand.

Technical Summary

1. Problem Statement

Recent discoveries of superconductivity in 2D van der Waals heterostructures, particularly Bernal bilayer graphene (BBG) proximitized with transition metal dichalcogenides (TMDs) like WSe$_2$, have revealed a significant violation of the Pauli paramagnetic limit. These systems exhibit high resilience to in-plane magnetic fields and, in some cases, magnetic field-induced superconductivity.

However, a comprehensive theoretical framework to interpret the upper critical field ($H_{c2}$) data in these systems is lacking. The complexity arises from the interplay of:

• Spin-Orbit Coupling (SOC): Specifically Ising and Rashba types induced by the substrate.
• Depairing Mechanisms: Zeeman splitting (spin) and orbital coupling.
• Pairing Symmetry: The uncertainty regarding whether the superconducting order parameter is spin-singlet or spin-triplet.
• Experimental Discrepancies: Previous fits to experimental data often yielded inconsistent or unphysical parameters (e.g., Rashba SOC magnitudes or Landé $g$-factors).

The authors aim to resolve these challenges by developing a general, analytically tractable framework to calculate $H_{c2}(T)$ and extract microscopic parameters from experimental data.

2. Methodology

The authors developed an effective low-energy microscopic model tailored for quasi-2D layered superconductors with broken inversion symmetry.

• Hamiltonian Construction:

  • They start from a faithful 4-band (8-band with SOC) microscopic model for BBG.
  • Using a systematic expansion (treating interlayer coupling $\gamma_1$ and velocity $v$ as large parameters, and SOC as perturbations), they project the Hamiltonian onto a low-energy 2-band subspace near the trigonal warping pockets.
  • The effective Hamiltonian includes:
    • Quasiparticle dispersion $\epsilon_\xi(k)$.
    • Ising SOC ($g_I \propto \lambda_I \hat{z}$): Out-of-plane spin splitting.
    • Rashba SOC ($g_R \propto \lambda_R (k_y, -k_x, 0)$): In-plane spin splitting.
    • Zeeman term ($b \propto g\mu_B B \hat{x}$): In-plane magnetic field.
    • Orbital coupling ($z_0 \propto g_{orb}\mu_B B$): Inter-layer coupling induced by the field.

• Gap Equation Formulation:

  • They derive a linearized gap equation for both spin-singlet and spin-triplet pairing channels.
  • The gap function is parameterized as $\hat{\Delta} = (d_0 s_0 + \mathbf{d} \cdot \mathbf{s}) i s_y$.
  • The equation is solved using Matsubara summation, resulting in a unified implicit equation involving the digamma function $\Psi(x)$ and Fermi surface averages.
  • A key parameter $\chi$ is introduced to distinguish singlet and triplet behaviors based on the misalignment of electron spins in the Cooper pair.

• Analytical and Numerical Limits:

  • The authors analyze limits for $T \to 0$ and $T \to T_c$.
  • They derive scaling laws for the Pauli Limit Violation Ratio (PVR), defined as $PVR = B_{c2}/B_P$ (where $B_P$ is the standard Pauli limit).
  • They perform numerical fits to experimental data from four recent BBG/WSe$_2$ experiments [1–4].

3. Key Contributions

1. Unified Theoretical Framework: The paper provides a general analytical framework that separates the effects of Ising SOC, Rashba SOC, orbital depairing, and Zeeman splitting on $H_{c2}(T)$ for both singlet and triplet order parameters.
2. Distinct Behavior of Triplet Pairing: The authors demonstrate that spin-triplet pairing exhibits qualitatively different $H_{c2}(T)$ curves depending on the orientation of the $\mathbf{d}$-vector relative to the magnetic field and SOC.
   • Specifically, for $\mathbf{d}$-vectors perpendicular to the field and SOC (e.g., $d_y, d_z$ with Ising SOC), the system can exhibit non-monotonic behavior, where superconductivity is induced or enhanced by the magnetic field at intermediate fields.
3. Microscopic Derivation for BBG: They explicitly derive how the microscopic parameters of BBG (trigonal warping, displacement field $U$) renormalize the effective SOC and $g$-factor in the low-energy limit.
   • Crucially, they show that the effective Ising SOC remains largely unrenormalized, while the effective Rashba SOC is suppressed by the displacement field.
4. Resolution of Experimental Discrepancies: By applying their model to experimental data, they identify a systematic discrepancy in previous interpretations regarding the Landé $g$-factor.

4. Key Results

• Ising SOC Dominance: The analysis confirms that Ising SOC is the dominant mechanism enhancing the upper critical field in BBG/WSe$_2$, consistent with experimental interpretations.
• Negligible Effective Rashba SOC: The fitted values for the effective Rashba SOC ($\lambda_R$) are found to be negligible compared to the bare layer-induced values. This aligns with their theoretical derivation showing that the projection to the low-energy pockets suppresses the effective Rashba term.
• Enhanced Landé $g$-Factor:
  • The most significant finding is that fitting the experimental $H_{c2}(T)$ data requires a Landé $g$-factor significantly larger than 2 (ranging from $g/g_0 \approx 2.5$ to $3.5$).
  • Standard non-interacting renormalization due to the displacement field can only account for an $\sim 8\%$ increase.
  • The authors propose that this enhancement arises from electron-electron interaction effects enhancing the Zeeman energy scale, or potentially a mismatch between the BCS pair-formation temperature ($T_c^{BCS}$) used in the model and the Berezinskii-Kosterlitz-Thouless (BKT) transition temperature observed experimentally.
• Singlet vs. Triplet Signatures:
  • Singlet: Ising SOC increases the PVR at all temperatures.
  • Triplet: The response is highly anisotropic. For certain $\mathbf{d}$-vector orientations, the system shows a "re-entrant" superconducting phase where $H_{c2}$ increases with temperature or field in specific regimes, a feature absent in singlet pairing.

5. Significance

• Diagnostic Tool: The framework serves as a powerful tool to distinguish between singlet and triplet superconducting order parameters in 2D materials by analyzing the shape of the $H_{c2}(T)$ curve.
• Understanding Field-Induced SC: The theoretical prediction of non-monotonic $H_{c2}$ curves for triplet pairing provides a mechanism to understand "field-induced" or "field-reinforced" superconductivity observed in recent experiments.
• Implications for Quantum Materials: The finding of an interaction-enhanced $g$-factor in these 2D systems suggests that many-body correlations play a critical role in their magnetic response, which is vital for designing future spintronic and topological quantum devices based on graphene and TMD heterostructures.
• General Applicability: While focused on BBG, the methodology is applicable to other 2D superconductors, including rhombohedral multilayer graphene and TMDs, provided they can be mapped to the effective Hamiltonian derived.

In summary, this paper bridges the gap between complex microscopic models and experimental $H_{c2}$ data, revealing that the robustness of superconductivity in these 2D systems is driven by strong Ising SOC and likely enhanced by many-body interactions that significantly modify the effective $g$-factor.