---
title: '2H-NbSe2-xSx: Mixed-Chalcogen Superconductivity'
url: https://www.emergentmind.com/topics/2h-nbse2-xsx
type: topic
---

# 2H-NbSe2-xSx: Mixed-Chalcogen Superconductivity

2H-NbSe\(_{2-x}\)S\(_x\) is a mixed-chalcogen layered dichalcogenide in the 2H polytype, spanning pure 2H-NbSe\(_2\) at \(x=0\), intermediate Se/S-substituted compositions, and the mixed-chalcogen analogue 2H-NbSeS at \(x\approx 1\). Across this series, random Se/S substitution preserves the global hexagonal \(P6_3/mmc\) structure while modifying the normal-state band structure, suppressing long-range charge-density-wave coherence, and reshaping the superconducting response to both nonmagnetic and magnetic disorder. Recent work has emphasized two closely related themes: first, that 2H-NbSeS realizes a fully gapped dirty-limit two-band superconductor with two nodeless gaps of different magnitudes; second, that in Fe-doped 2H-NbSe\(_{2-x}\)S\(_x\), gapless superconductivity can emerge at remarkably low magnetic impurity concentrations once \(x>0\) [2508.00435] [2606.19107].

## 1. Structural framework and substitution series

2H-NbSe\(_{2-x}\)S\(_x\) crystallizes in the 2H polytype with hexagonal space group \(P6_3/mmc\). Each layer contains one Nb plane sandwiched between two Se/S planes, and the layers are weakly coupled by van der Waals forces [2508.00435]. This layered architecture is central to both the electronic anisotropy and the long-range spatial extent of impurity-induced states discussed later.

For the substitution series characterized by XRD and EDS, the reported lattice parameters evolve systematically with sulfur content:

| \(x\) | \(a=b\) (nm) | \(c\) (nm) |
|---|---:|---:|
| 0.0 | 0.34451 | 1.2542 |
| 0.2 | 0.34327 | 1.2506 |
| 0.4 | 0.34197 | 1.2474 |
| 0.6 | 0.34099 | 1.2423 |
| 0.8 | 0.33984 | 1.2387 |

The same measurements report \(T_c=7.2\) K for \(x=0.0\), \(6.2\) K for \(x=0.2\), \(5.0\) K for \(x=0.4\), \(4.6\) K for \(x=0.6\), and \(3.3\) K for \(x=0.8\) [2508.00435]. In a separate study of single-crystalline 2H-NbSeS, the material becomes superconducting below \(6.0\) K with moderate magnetic anisotropy [2606.19107]. A plausible implication is that the superconducting transition temperature in this family is sensitive not only to nominal sulfur content but also to the disorder landscape and the experimental context used to define the compound.

Random Se/S occupation does not change the nominal electron count or the global \(P6_3/mmc\) symmetry in 2H-NbSeS, but it does introduce strong local potential fluctuations [2606.19107]. This distinction between preserved average crystallography and altered local electronic environment is a recurring feature of the series.

## 2. Electronic structure, charge-density-wave suppression, and nesting

The electronic structure evolves substantially under Se\(\rightarrow\)S substitution. In pure 2H-NbSe\(_2\), density functional theory on a 5-layer slab with PBE + Grimme-D3 finds that the central Se-derived pocket at \(\Gamma\) is nearly spherical and that the Nb-derived tubular pockets exhibit strong \(k_x\)–\(k_y\) warping along \(k_z\), so that the \(\Gamma\)–K–M and A–H–L dispersions differ [2508.00435]. For \(x=0.2\), the Nb sub-bands become more two-dimensional: \(k_z\)-warping is strongly suppressed, the dispersions at \(\Gamma\)–K–M and A–H–L become nearly identical, and the Se/S pocket shrinks slightly. As \(x\) increases up to \(0.8\), the normal-state DOS at \(E_F\) is reduced and van Hove singularities are smeared by substitutional disorder; a tight-binding model captures this trend [2508.00435].

A complementary description is provided for 2H-NbSeS. In pristine 2H-NbSe\(_2\), an incommensurate CDW sets in below \(T_{CDW}\approx 33\) K, reconstructing large portions of the Fermi surface and opening a partial gap \(\Delta_{CDW}\approx 3\) meV on nested bands around the K–M region. In 2H-NbSeS, random Se/S substitution destroys the phase coherence necessary for long-range CDW order without changing the nominal electron count or the global symmetry [2606.19107]. Experimentally, the absence of any anomaly in \(\rho(T)\) and the lack of superlattice Bragg reflections in XRD confirm CDW suppression; a Bloch–Grüneisen fit yields \(\rho_0=11.68\,\mu\Omega\,\mathrm{cm}\), \(\Theta_D=195\) K, and \(\mathrm{RRR}\approx 9.5\) [2606.19107].

Within that same study, the three-sheet Nb-4d/Se-4p Fermi surface predicted by DFT for 2H-NbSe\(_2\)—two cylindrical Nb-derived sheets at \(\Gamma\)-A and K-H plus a small Se-\(p_z\) pocket at \(\Gamma\)—remains essentially unmodified except for the removal of CDW-induced backfolded replicas [2606.19107]. Taken together, the two studies indicate that Se/S substitution both suppresses CDW reconstruction and reshapes the detailed dispersion. This suggests that the relevant low-energy physics cannot be reduced either to a rigid-band picture or to disorder-averaged CDW phenomenology alone.

Scanning tunneling microscopy further resolves how the dominant in-gap scattering channels change with substitution. In pure NbSe\(_2\), in-gap conductance maps around isolated Fe impurities are dominated by the atomic lattice vector \(q_{at}\) and the CDW vector \(q_{CDW}=0.34\,\mathrm{\AA}^{-1}\). For \(x=0.2\), new wavevectors appear at biases of \(0.7\)–\(1.0\) mV: \(q_1\) and \(q_2\) connect M–L points of the Nb Fermi surface, \(q_3\) lies along \(\Gamma\)–M nesting of the hexagonal Nb pocket and is identified with the Johannes–Mazin nesting vector, and \(q_4\) corresponds to near-circular scattering from the central Se/S pocket [2508.00435]. DFT correlates \(q_3\) with flattened hexagonal segments, \(q_4\) with the central pocket, and \(q_1/q_2\) with M–L curvature. Extended Data 8 shows an exponential CDW peak decay, leading to the conclusion that S substitution decouples layers, suppresses long-range CDW coherence, and replaces pure CDW scattering with band-structure-specific nesting vectors that dominate pair-breaking quasiparticle interference [2508.00435].

## 3. Multiband superconductivity in 2H-NbSeS

The superconducting state of 2H-NbSe\(_2\) has been the subject of sustained debate. Angle-resolved photoemission spectroscopy data have been interpreted as evidence for multiband superconductivity, whereas scanning tunneling microscope experiments have been interpreted in terms of strongly anisotropic single-band superconductivity, with the CDW mimicking multigap character by reconstructing the Fermi surface and modifying the superconducting gap distribution [2606.19107]. In 2H-NbSeS, where long-range CDW order is suppressed but the layered \(P6_3/mmc\) structure is preserved, the same issue can be examined with reduced CDW-related ambiguity.

Several probes converge on a fully gapped superconducting state with two nodeless gaps of different magnitudes in 2H-NbSeS [2606.19107]. The upper critical field \(H_{c2}(T)\), determined by the 90%-onset criterion in \(\rho(H,T)\) for both \(H\parallel ab\) and \(H\perp ab\), exhibits pronounced upward curvature that cannot be described within a single-band framework. A conventional single-band WHH fit fails at low temperature. By contrast, the dirty-limit two-band Usadel model of Gurevich captures the data through the implicit equation
$$
a_0\bigl[\ln t + U(h)\bigr]\bigl[\ln t + U(\eta h)\bigr]
+a_1\bigl[\ln t + U(h)\bigr]
+a_2\bigl[\ln t + U(\eta h)\bigr] = 0,
$$
with
$$
t=\frac{T}{T_c},\qquad
h=\frac{H D_1}{2\Phi_0 T},\qquad
\eta=\frac{D_2}{D_1},\qquad
U(x)=\psi\bigl(x+\tfrac12\bigr)-\psi(\tfrac12),
$$
and
$$
a_0=\lambda_{11}\lambda_{22}-\lambda_{12}\lambda_{21},\qquad
a_1=\lambda_{11}+\lambda_{22}-(\lambda_{12}+\lambda_{21}),\qquad
a_2=\lambda_{12}+\lambda_{21}.
$$

Fitting yields, for \(H\parallel ab\), \(\lambda_{11}=1.53\), \(\lambda_{22}=0.87\), \(\lambda_{12}=\lambda_{21}=0.29\), and \(\eta=6.17\); for \(H\perp ab\), \(\lambda_{11}=1.22\), \(\lambda_{22}=1.09\), \(\lambda_{12}=\lambda_{21}=0.14\), and \(\eta=7.85\) [2606.19107]. The zero-temperature critical fields extrapolate to \(H_{c2}^{\parallel ab}(0)=14.55\) T and \(H_{c2}^{\perp ab}(0)=4.68\) T. Within anisotropic GL theory, these correspond to \(\xi_{ab}\approx 8.4\) nm, \(\xi_c\approx 2.7\) nm, and \(\gamma=\xi_{ab}/\xi_c\approx 3.1\) [2606.19107]. The in-plane upper critical field exceeds the weak-coupling Pauli limit, which the same study states as \(1.84\,T_c\approx 11\) T.

Independent thermodynamic and electrodynamic probes support the same two-gap interpretation. From lower-critical-field data, the superfluid density is modeled as
$$
H_{c1}(T)=H_{c1}(0)\,\rho_s(T),\qquad
\rho_s(T)=x_1\,\rho_1(T)+x_2\,\rho_2(T),
$$
with
$$
\rho_i(T)=1+2\int_{\Delta_i(T)}^\infty \frac{\partial f}{\partial E}\,
\frac{E\,dE}{\sqrt{E^2-\Delta_i^2(T)}},
$$
and
$$
\Delta_i(T)=\Delta_i(0)\tanh\!\Bigl[1.82\,(1.018/t-1)^{0.51}\Bigr].
$$
For \(H\parallel ab\), the fit gives \(2\Delta_1/k_B T_c=3.98\), \(2\Delta_2/k_B T_c=1.60\), \(x_1=0.80\), \(x_2=0.20\), and \(H_{c1}(0)=169.7\) Oe; for \(H\perp ab\), \(2\Delta_1/k_B T_c=4.00\), \(2\Delta_2/k_B T_c=1.58\), \(x_1=0.84\), \(x_2=0.16\), and \(H_{c1}(0)=75.1\) Oe [2606.19107]. A single-gap model cannot reproduce the pronounced low-temperature upward curvature of \(\rho_s(T)\), whereas the two-gap form yields excellent agreement.

Tunnel-diode-oscillator penetration-depth measurements likewise require a two-gap description over the full temperature range. At low \(T\), a single-gap expression,
$$
\Delta\lambda(T)=\lambda(0)\sqrt{\frac{\pi\Delta_0}{2k_B T}\;
\exp\!\Bigl(-\frac{\Delta_0}{k_B T}\Bigr)},
$$
gives \(2\Delta_0/k_B T_c=4.08\), but the full-range fit yields \(2\Delta_1/k_B T_c=4.24\), \(2\Delta_2/k_B T_c=1.46\), \(x_1=0.85\), and \(x_2=0.15\) [2606.19107]. Electronic specific heat at zero field is also described by a two-gap fit,
$$
C_e/\gamma_n T_c = t\,\frac{d}{dt}\Bigl[\;-\frac{6}{\pi^2}\,\frac{\Delta(0)}{k_B T_c}\!\int (f\ln f + (1-f)\ln(1-f))\,dy\Bigr],
$$
with extracted values \(2\Delta_1/k_B T_c=3.68\), \(2\Delta_2/k_B T_c=1.28\), \(x_1=0.87\), and \(x_2=0.13\) [2606.19107]. The normal-state fit yields \(\gamma_n=12.40\,\mathrm{mJ/mol\,K^2}\) and \(\Theta_D=236\) K, and only the \(s+s\) fit captures both \(\Delta C/\gamma_n T_c\approx 1.4\) and the low-temperature behavior.

These results are interpreted as evidence that multiband pairing is intrinsic to the Nb-4d Fermi surface and not a CDW artifact in 2H-NbSeS [2606.19107].

## 4. Extremely dilute magnetic disorder and gapless superconductivity

A distinct aspect of 2H-NbSe\(_{2-x}\)S\(_x\) emerges when the system is doped with extremely dilute Fe impurities. In the STM study, samples with sulfur substitution \(x=0\ldots 0.8\) were doped with Fe at \(n\approx 0.020\) at.%—approximately one Fe per 3,000 unit cells—and measured at \(T\lesssim 0.05\,T_c\), down to below 100 mK [2508.00435].

In clean 2H-NbSe\(_2\) (\(x=0\)) with Fe, the response is that expected for localized magnetic impurity states in a gapped superconductor: well-defined Yu–Shiba–Rusinov peaks appear at \(\pm 0.2\,\Delta_0\), with \(\Delta_0\approx 5\) meV, while the DOS at zero bias vanishes far from impurities [2508.00435]. For any \(x>0\), however, the phenomenology changes qualitatively. At Fe concentrations as low as \(n\approx 0.014\)–\(0.032\) at.%, zero-bias conductance maps show a finite, spatially homogeneous DOS at \(E_F\) even in regions many coherence lengths, \(\xi\approx 7\) nm, from any Fe. Tunneling spectra far from impurities evolve from a full gap with \(G(0)\approx 0\) to \(G(0)\approx 0.05\)–\(0.15\), normalized to the above-gap conductance, signaling gapless superconductivity [2508.00435].

The reported threshold Fe concentration for the gapless state is \(n\approx 0.01\)–\(0.02\) at.%, which is far below the few-percent level normally required in conventional \(s\)-wave superconductors [2508.00435]. The same study explicitly contrasts this with the few-percent threshold for gapless behavior in bulk conventional \(s\)-wave systems discussed by Shiba and by Hauser et al. This constitutes the central anomaly of the series: apparently negligible magnetic-impurity densities can produce a bulk-like finite DOS at the Fermi level once Se/S substitution is present.

The spatial homogeneity of the zero-bias DOS is particularly significant. The measurements show that the finite \(E_F\) DOS is not confined to local impurity neighborhoods but persists many coherence lengths from identifiable Fe sites [2508.00435]. This suggests a collective impurity-band or overlap mechanism rather than a simple superposition of isolated in-gap bound states.

## 5. Microscopic modeling of disorder, impurity states, and in-gap scattering

The theoretical description of the gapless state uses a real-space BCS Hamiltonian on a triangular Nb lattice with hopping up to fifth nearest neighbors,
\[
\{t_1,\ldots,t_5\}=\{-23,\,-128.75,\,-2.2,\,7.5,\,-3\}\ \mathrm{meV},
\]
and chemical potential \(\mu=-275\) meV [2508.00435]. Magnetic impurities are modeled as local exchange potentials \(V_m\) at random sites \(r_m\), while nonmagnetic Se–S disorder is represented by local repulsive potentials \(V_p\) at random sites \(r_p\). In the Nambu basis \(\Psi(r)=(c_\uparrow,c_\downarrow^\dagger)^T\), the BdG Hamiltonian is
$$
H_{\mathrm{BdG}}=
\begin{bmatrix}
H_0-\mu & \Delta(r)\\
\Delta^*(r) & -H_0+\mu
\end{bmatrix}.
$$
The gap is determined self-consistently through
$$
\Delta(r)=V_{SC}\sum_{E_n>0}u_n(r)v_n^*(r)\tanh(E_n/2k_B T),
$$
with \(V_{SC}=-90\) meV chosen so that the homogeneous \(\Delta_0\approx 5\) meV. The local DOS is
$$
N(r,\omega)= -\frac{1}{\pi}\,\mathrm{Im}\sum_{n,\sigma}
\left[
\frac{|u^n_\sigma(r)|^2}{\omega-E_n+i\eta}
+
\frac{|v^n_\sigma(r)|^2}{\omega+E_n+i\eta}
\right],
$$
using \(\eta\approx 0.02\,\Delta_0\) as numerical broadening [2508.00435].

The key BdG result is that even at \(n\) up to \(0.06\)–\(0.6\) at.% and \(x\) up to \(0.8\), dilute magnetic and nonmagnetic disorder cooperate to fill in the gap, producing finite \(\langle N(\omega=0)\rangle\) and six-fold-smeared YSR states [2508.00435]. The experimentally observed finite zero-energy DOS appears already at \(n\approx 0.014\)–\(0.02\) at.%, while theoretically comparable values require \(n\approx 0.06\)–\(0.2\) at.%; the study nonetheless emphasizes that both concentrations remain extremely small [2508.00435].

The proposed origin of this high sensitivity is that YSR bound states in layered TMDs have long real-space tails, with \(r_0\approx\) a few nm. These tails overlap at very low impurity densities and, together with S-induced band-structure changes, generate a homogeneous in-gap background [2508.00435]. This mechanism links the disorder response to low dimensionality and to the material-specific Fermi surface rather than to impurity concentration alone.

The same framework is supported by the in-gap QPI data. In pure NbSe\(_2\), impurity scattering is dominated by \(q_{at}\) and \(q_{CDW}\), whereas for \(x=0.2\) the dominant in-gap wavevectors are \(q_1\)–\(q_4\), which trace nesting and curvature of the modified Nb and Se/S pockets [2508.00435]. The conclusion is that pair breaking for \(x>0\) is dictated by substitution-modified band structure rather than by the dominant CDW interactions of pure 2H-NbSe\(_2\).

## 6. Comparison with pure 2H-NbSe\(_2\) and broader implications

Relative to representative literature values for pure 2H-NbSe\(_2\), the 2H-NbSeS study reports \(T_c=6.0\) K rather than \(7.2\) K, \(H_{c1}^{\parallel ab}(0)\approx 170\) Oe rather than \(\sim 120\) Oe, \(H_{c1}^{\perp ab}(0)\approx 75\) Oe rather than \(\sim 60\) Oe, \(H_{c2}^{\parallel ab}(0)\approx 14.6\) T rather than \(\sim 12.2\) T, and \(H_{c2}^{\perp ab}(0)\approx 4.7\) T rather than \(\sim 3.8\) T [2606.19107]. The coherence lengths remain similar, \(\xi_{ab}(0)\sim 8.4\) nm and \(\xi_c(0)\sim 2.7\) nm versus literature values \(\sim 8.0\) nm and \(\sim 2.5\) nm, while the anisotropy remains close to \(\gamma\sim 3.1\)–\(3.2\). The large diffusivity ratio, however, changes substantially: the literature value for pure 2H-NbSe\(_2\) is approximately \(1\)–\(3\), whereas 2H-NbSeS gives \(6\)–\(8\), which is described as strongly dirty [2606.19107].

The same comparison emphasizes several trends upon Se\(\rightarrow\)S substitution: suppression of CDW entirely, removal of Fermi-surface reconstruction and CDW-induced gap anisotropy, modest reduction of \(T_c\), increased residual scattering rate, enhancement of \(H_{c2}^{\parallel ab}(0)\) above the weak-coupling Pauli limit, and robust two-gap ratios and weights [2606.19107]. In that view, multiband pairing is an intrinsic property of the Nb-4d-derived electronic structure.

The Fe-doped study extends the significance of the substitution series beyond multiband phenomenology. It demonstrates that even extremely dilute magnetic moments can render an \(s\)-wave state gapless in the presence of nonmagnetic substitutional disorder if the normal-state band structure favors extended impurity states and nesting [2508.00435]. This is presented as evidence for the inadequacy of Anderson’s theorem once realistic band structures and low-dimensionality effects are included. The same work further suggests that the robustness of superconducting gaps in other layered TMDs, and by extension 2D superconductors, must be reassessed by incorporating material-specific DFT dispersions and the interplay of CDW, electron–phonon coupling, and magnetic correlations [2508.00435].

A common misconception is that random chalcogen substitution merely adds nonmagnetic disorder while leaving superconductivity qualitatively conventional. The available results point instead to a more specific conclusion. In 2H-NbSe\(_{2-x}\)S\(_x\), Se/S substitution suppresses long-range CDW order, alters dimensionality and nesting, produces strong band-dependent scattering, and changes how dilute magnetic moments couple to the condensate [2508.00435] [2606.19107]. This suggests that the series is best understood not as a simple disorder-broadened version of 2H-NbSe\(_2\), but as a platform in which band-structure engineering and disorder cooperate to expose both intrinsic multiband superconductivity and an unusual route to gapless superconductivity.

Source: https://www.emergentmind.com/topics/2h-nbse2-xsx