Ising Superconductivity in 2D Materials

• Ising superconductivity is a phenomenon where spin–orbit coupling pins electron spins out-of-plane, greatly suppressing pair breaking by in-plane magnetic fields.
• Monolayer TMDs like MoS₂ and NbSe₂ exhibit in-plane critical fields several times above the Pauli limit, highlighting the robust nature of their spin–valley locked superconducting state.
• Engineering approaches extend Ising superconductivity via spin–orbital locking in centrosymmetric crystals and heterostructures, enabling mixed singlet–triplet pairing and potential topological phases.

Ising superconductivity is a superconducting regime in which spin–orbit coupling generates an effective internal Zeeman field that pins electronic spins predominantly perpendicular to a two-dimensional plane, thereby strongly suppressing paramagnetic pair breaking by external in-plane magnetic fields. In its canonical form, realized in monolayer transition-metal dichalcogenides (TMDs), the effect is tied to spin–valley locking at the $K$ and $K'$ valleys; more recent work has generalized the concept to centrosymmetric two-dimensional crystals through spin–orbital locking, to acentric bulk layered polymorphs, and to interface-engineered thin films. Across these settings, the central phenomenology is a conventional or predominantly spin-singlet superconducting condensate formed on spin-split bands with strong out-of-plane polarization, accompanied by anomalously large in-plane upper critical fields, mixed singlet–triplet correlations, and multiple routes to topological superconductivity (Yuan et al., 2016, Wang et al., 2019, Patra et al., 11 Jun 2025, Yi et al., 6 Sep 2025).

1. Symmetry setting and microscopic origin

In monolayer TMDs such as MoS$_2$ and NbSe$_2$, heavy transition-metal atoms provide strong atomic spin–orbit coupling, while the monolayer crystal structure lacks inversion symmetry and, in the standard low-energy description, supports in-plane electric fields. The resulting SOC is commonly summarized schematically by

$$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$

so that the effective SOC field points mainly along $\pm \hat z$, yielding an Ising term of the form

$H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$

Because time-reversal symmetry maps $\mathbf{k}$ to $-\mathbf{k}$, the SOC field reverses sign between opposite valleys, producing spin–valley locking: one spin orientation is favored at $K$, the opposite at $K'$0 (Yuan et al., 2016, Zhou et al., 2015).

The low-energy Hamiltonian near a TMD valley is often written in the form

$K'$1

with $K'$2 for $K'$3, $K'$4, and $K'$5 the Ising SOC strength. In gated MoS$K'$6, the fitting value $K'$7 implies a band-edge spin splitting $K'$8, while the broader TMD literature described here places conduction- or valence-band spin splittings in the $K'$9 range or larger (Yuan et al., 2016).

A common misconception is that this mechanism requires only the absence of inversion symmetry and valley degrees of freedom. That statement is accurate for the original, now often called type-I, form of Ising superconductivity, where the SOC is written as

$_2$0

with $_2$1 labeling $_2$2. However, later work identified a distinct type-II mechanism in multi-orbital two-dimensional crystals with rotational symmetry $_2$3 ($_2$4), where the same algebraic structure arises from spin–orbital rather than spin–valley locking: $_2$5 Here the “valley index” is replaced by an orbital pseudospin $_2$6, the protection is strongest near time-reversal-invariant momenta such as $_2$7, and inversion symmetry may be present or absent. In SnH, the fitted value $_2$8 corresponds to an effective Zeeman field of order $_2$9 T, and high-throughput first-principles screening identified about one hundred candidate type-II materials (Wang et al., 2019).

2. Pairing structure and superconducting correlations

The defining superconducting state in the original TMD setting is not an exotic interaction-driven triplet phase but a predominantly conventional spin-singlet condensate living on bands already split by Ising SOC. In Nambu basis, the mean-field Bogoliubov–de Gennes Hamiltonian is

$_2$0

and the basic pairing matrix is taken as momentum-independent $_2$1-wave singlet pairing,

$_2$2

This is the sense in which Ising superconductivity is usually defined: conventional pairing on top of strongly out-of-plane spin-polarized, spin–valley-locked normal-state bands (Yuan et al., 2016).

Non-centrosymmetric SOC nevertheless mixes singlet and triplet pairing correlations in the anomalous propagator. Writing the Gor’kov pairing matrix in singlet–triplet form,

$_2$3

the Ising case yields a $_2$4-vector along the Ising axis, $_2$5. If the spin quantization axis is rotated into the plane by an angle $_2$6, the pairing matrix becomes

$_2$7

so the triplet sector corresponds to equal-spin Cooper pairs whose spin polarization points in-plane. This point is central: Ising superconductors can be interaction-wise $_2$8-wave and still host equal-spin triplet correlations in their pairing amplitudes (Zhou et al., 2015, Yuan et al., 2016).

That distinction also clarifies what the term does not mean. It does not denote an Ising order parameter in the Landau sense, and it does not imply purely triplet superconductivity. Rather, the “Ising” label refers to the highly anisotropic spin structure of the paired quasiparticles. This interpretation extends to later generalized settings: in type-II Ising superconductivity the locking is orbital rather than valley based, but the superconducting state is still described primarily as time-reversal-invariant pairing of Kramers partners with strongly suppressed in-plane Zeeman splitting (Wang et al., 2019).

3. Critical fields and experimental hallmarks

The canonical quantitative benchmark is the Pauli paramagnetic limit. In the cited weak-coupling treatments it is written either as

$_2$9

or

$$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$0

with field in tesla when $$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$1 is in kelvin. In atomically thin superconductors under in-plane field, orbital depairing is strongly suppressed, so a conventional spin-singlet state would ordinarily be limited by this Pauli scale. Ising SOC alters the situation because an in-plane Zeeman term $$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$2 acts transverse to a much larger internal SOC field $$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$3, so the spins cant only weakly and the effective pair-breaking splitting is drastically reduced (Yuan et al., 2016, Lu et al., 2015).

This mechanism is directly reflected in experiment. In gated MoS$$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$4, samples with $$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$5, $$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$6, and $$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$7 showed in-plane $$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$8 up to about six times the Pauli limit, and theoretical curves using $$\mathbf{S}(\mathbf{k}) \sim \mathbf{k}\times \mathbf{E},$$9 plus a gate-dependent Rashba term reproduced the data. Monolayer NbSe$\pm \hat z$0 exhibited a similar $\pm \hat z$1 enhancement, while out-of-plane critical fields remained much smaller because they are limited by orbital effects and vortex physics. In gated MoS$\pm \hat z$2, angle-dependent $\pm \hat z$3 displayed the cusp characteristic of two-dimensional superconductivity, and the extracted $\pm \hat z$4 in one device established Berezinskii–Kosterlitz–Thouless behavior (Yuan et al., 2016, Lu et al., 2015).

Intrinsic metallic TMDs push this further. In monolayer TaS$\pm \hat z$5, $\pm \hat z$6 implies $\pm \hat z$7, yet $\pm \hat z$8 exceeded $\pm \hat z$9 for $H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$0, and fitting yielded

$H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$1

about $H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$2. Monolayer NbSe$H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$3 displayed a smaller but still very large enhancement, and few-layer TaS$H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$4 and NbSe$H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$5 retained $H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$6 by factors of $H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$7–$H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$8 or more, showing that local Ising protection can survive when spin–orbit coupling dominates interlayer hopping (Barrera et al., 2017).

The same hallmark now appears beyond monolayer TMDs. In acentric bulk 4H-NbSe$H_{\text{Ising}} \propto S(\mathbf{k})\,\sigma_z.$9, two-band analysis gave

$\mathbf{k}$0

with anisotropy $\mathbf{k}$1. Using $\mathbf{k}$2, the Pauli limit is $\mathbf{k}$3, so the in-plane critical field exceeds the Pauli limit by nearly a factor of three, while the out-of-plane critical field stays below it. In the graphene/trilayer Ga/SiC heterostructure, electrical transport found

$\mathbf{k}$4

approximately $\mathbf{k}$5 times the Pauli paramagnetic limit $\mathbf{k}$6, and a Ginzburg–Landau extrapolation gave $\mathbf{k}$7 (Patra et al., 11 Jun 2025, Yi et al., 6 Sep 2025).

4. Material platforms and variants

The original Ising-superconducting platforms were monolayer or effectively monolayer TMDs. Gated MoS$\mathbf{k}$8 realizes a surface-confined superconducting layer with carriers concentrated at the topmost layer; monolayer NbSe$\mathbf{k}$9 and TaS$-\mathbf{k}$0 are intrinsic superconductors in the atomic-layer limit; and the combination of strong SOC, broken inversion symmetry, and suppressed orbital depairing makes these systems the cleanest realizations of type-I Ising superconductivity (Lu et al., 2015, Barrera et al., 2017).

Later work broadened the materials landscape in several orthogonal directions. Type-II Ising superconductivity showed that inversion breaking is not essential if multi-orbital rotational symmetry generates spin–orbital locking near a time-reversal-invariant momentum. Bulk non-centrosymmetric 4H-NbSe$-\mathbf{k}$1 demonstrated that Ising phenomenology can survive in an intrinsically three-dimensional crystal if the stacking itself is acentric. Interface engineering then extended the idea to light-element systems, as in graphene/trilayer Ga/SiC, where orbital hybridization with the substrate creates split Fermi surfaces with Ising-type spin textures at $-\mathbf{k}$2 and $-\mathbf{k}$3 (Wang et al., 2019, Patra et al., 11 Jun 2025, Yi et al., 6 Sep 2025).

Platform	Dominant mechanism	Representative observation
Gated MoS$-\mathbf{k}$4	Type-I spin–valley locking	In-plane $-\mathbf{k}$5 up to about six times $-\mathbf{k}$6 (Yuan et al., 2016)
Monolayer NbSe$-\mathbf{k}$7, TaS$-\mathbf{k}$8	Type-I Ising SOC in intrinsic TMDs	$-\mathbf{k}$9 in monolayer TaS$K$0 (Barrera et al., 2017)
SnH and related 2D crystals	Type-II spin–orbital locking	$K$1; about one hundred candidates identified (Wang et al., 2019)
4H-NbSe$K$2	Bulk acentric valley-selective spin polarization	$K$3 (Patra et al., 11 Jun 2025)
NbSe$K$4F$K$5	Symmetry-preserving fluorination of monolayer NbSe$K$6	$K$7 vs $K$8 in the CDW phase (Wu et al., 2024)
Graphene/trilayer Ga/SiC	Interfacial orbital hybridization-induced Ising texture	$K$9 at $K'$00 (Yi et al., 6 Sep 2025)

Materials engineering has become a central theme. In monolayer NbSe$K'$01, symmetric fluorination on both sides preserves out-of-plane mirror symmetry and retains a large spin splitting at $K'$02: $K'$03 in the CDW phase and $K'$04 in NbSe$K'$05F$K'$06. At the same time, fluorination suppresses the CDW, pushes a van Hove singularity close to the Fermi level, increases the density of states from $K'$07 to $K'$08 states/eV per formula unit, raises the isotropic EPC constant from $K'$09 to $K'$10, and yields an anisotropic Migdal–Eliashberg $K'$11, compared with $K'$12 for the CDW monolayer. The same calculations show that $K'$13 remains well above the Pauli limit and is larger than in pristine monolayer NbSe$K'$14 (Wu et al., 2024).

5. Topological, impurity, and proximity phenomena

A major consequence of the Ising-pinned spin structure is the emergence of topological superconductivity under proximity or field engineering. Because Ising superconductors carry equal-spin triplet correlations with spins pointing in-plane, they can induce effective $K'$15-wave pairing in a spin-polarized one-dimensional system. One proposal places a paramagnetic nanowire such as InSb on an Ising TMD under an in-plane field: the wire is spin-polarized by the field, while the TMD contributes equal-spin triplet correlations, producing a Kitaev-like one-dimensional topological superconductor. Numerical calculations in that setting show a full induced gap in the spectral function of an infinite wire and zero-energy end-localized Majorana modes in a finite wire (Zhou et al., 2015, Yuan et al., 2016).

A related route uses magnetic atoms. In monolayer NbSe$K'$16, sufficiently strong in-plane field can drive the host itself into a nodal topological superconducting phase when the Zeeman energy exceeds the superconducting gap but remains below $K'$17. The bulk spectrum then develops nodal points, and the edge hosts Majorana flat bands connecting their surface projections. Tunneling into such an edge is predicted to yield zero-bias conductance peaks, while the equal-weight electron–hole character of the edge Majorana modes makes the boundary act as a perfect spin filter (Yuan et al., 2016).

Magnetic impurities also provide unusually sensitive local probes. For a single magnetic impurity in an Ising superconductor, Yu–Shiba–Rusinov states display a strongly anisotropic magnetic-field response: zero-bias conductance peaks split readily for out-of-plane fields but remain robust to much larger in-plane fields, mirroring the Ising protection of the condensate itself. For a chain of magnetic impurities with moments parallel to the plane of the Ising superconductor, the low-energy YSR band can host topological superconductivity and Majorana fermions as a direct manifestation of Ising SOC-induced topological effects (Sharma et al., 2016).

More recent impurity work has pushed this logic in two directions. First, a local-probe theory of YSR states in Ising superconductors with in-plane field found distinct bound-state spectra and tunneling responses that differentiate the state from conventional superconductors and encode the underlying mixed pairing structure (Hein et al., 12 May 2026). Second, a theory for monolayer NbSe$K'$18 dosed with magnetic $K'$19 atoms predicted that the critical temperature can be slightly increased by an in-plane magnetic field in Cr-dosed samples: when easy-axis impurity spins are reoriented by the field, spin-conserving scattering is converted into spin-flip scattering, and because band spin splitting makes spin-flip processes require finite momentum transfer, pair breaking can be reduced rather than enhanced (Fang et al., 2023).

6. Conceptual scope, misconceptions, and open directions

The term “Ising superconductor” is most securely attached to superconductors whose normal state has strong out-of-plane spin locking and whose superconducting state is correspondingly protected against in-plane paramagnetic depairing. Within that established definition, several misunderstandings recur. Ising superconductivity is not synonymous with purely triplet pairing; the standard TMD realization is modeled with conventional $K'$20-wave spin-singlet pairing, although the resulting pairing correlations contain both singlet and triplet components. Nor is it restricted to monolayer non-centrosymmetric TMDs: type-II Ising superconductivity shows that inversion symmetry can be retained, while 4H-NbSe$K'$21 demonstrates that bulk acentric crystals can also realize the key phenomenology (Yuan et al., 2016, Wang et al., 2019, Patra et al., 11 Jun 2025).

Several open problems remain central. The microscopic pairing mechanism in gated MoS$K'$22 is still treated phenomenologically as $K'$23-wave electron–phonon pairing in most modeling, but its exact nature remains to be fully clarified. The predicted nodal topological phase in monolayer NbSe$K'$24 under in-plane field still awaits direct confirmation by ARPES, STM/STS, or spin-resolved tunneling. TMD-based heterostructures that induce Ising-type SOC into adjacent conventional superconductors, thereby enhancing $K'$25 or driving nodal topological phases, remain a promising but incompletely realized direction. Disorder, electron–electron interactions, and the full effect of Rashba components under gating or substrate asymmetry also remain incompletely understood (Yuan et al., 2016, Wickramaratne et al., 2023).

Defect and alloy physics have added a further layer of complexity. In monolayer niobium dichalcogenide alloys NbS$K'$26Se$K'$27, first-principles work found that both the density of states at the Fermi level and the proximity to magnetism are reduced relative to NbSe$K'$28, and that Se vacancies—likely magnetic pair-breaking defects—may form in large concentrations in NbSe$K'$29. This was advanced as an alternative explanation for the reported non-monotonic $K'$30 dependence on sulfur content, without invoking multifractality. The result suggests that defect chemistry and spin fluctuations can tune $K'$31 substantially without destroying the underlying Ising spin–valley structure (Wickramaratne et al., 2021).

More recent work has also generalized the label “Ising superconductivity” beyond SOC-driven TMD settings. In $K'$32-wave magnets with zero net magnetization but nonrelativistic, collinear spin splitting in momentum space, superconductivity has been argued to support only an Ising mixed-parity state in which each Cooper pair is a 50:50 singlet–triplet mixture. A follow-up study of a two-dimensional $K'$33-wave magnet found a leading $K'$34 Ising instability and showed that when the triplet amplitude exceeds the singlet one, the system enters a nodal topological superconducting phase with Majorana edge modes, while a Zeeman field perpendicular to the exchange field can induce a $K'$35 topological superconducting phase (Khodas et al., 27 Jan 2026, Kim et al., 3 May 2026).

Taken together, these developments define Ising superconductivity less as a single materials class than as a symmetry-protected superconducting response of spin-split electronic states. Its original realization in monolayer TMDs remains paradigmatic, but the concept now spans spin–valley and spin–orbital locking, bulk acentric crystals, impurity and heterostructure engineering, and even exchange-driven generalizations. The unifying criterion is the same throughout: superconductivity forms on bands whose spins are pinned predominantly out of plane, so that in-plane magnetic fields are unusually ineffective at breaking Cooper pairs.