NSSM-I-SC Junction: Theory and Conductance

• The NSSM-I-SC junction is a hybrid structure combining a nodal surface semimetal, a thin insulating barrier, and an s-wave superconductor that enables unusual quasiparticle tunneling phenomena.
• Theoretical analysis using the Bogoliubov–de Gennes framework yields explicit reflection amplitudes and π-period conductance oscillations arising from the momentum-dependent barrier parameter.
• Distinct transport signatures, tunable via gating and high-frequency irradiation, differentiate NSSM junctions from graphene-based systems and open paths for optoelectronic quantum devices.

A Nodal Surface Semimetal–Insulator–Superconductor (NSSM-I-SC) junction is a hybrid quantum structure comprising a nodal surface semimetal, an intervening thin insulating barrier, and an adjoining $s$-wave superconductor. The distinctive band topology of the NSSM—marked by two-dimensional nodal surfaces in momentum space—imparts unique electronic and thermal transport features at the junction interface. These features manifest in the behavior of quasiparticle tunneling, Andreev and normal reflections, and the resulting electrical and thermal conductance, which crucially distinguish NSSM-based junctions from their graphene or silicene counterparts. Recent theoretical analyses articulate the full Bogoliubov–de Gennes (BdG) framework for this system, produce explicit reflection amplitudes and conductance formulas, and describe clear experimental signatures for barrier-tunable quantum transport (Pandit et al., 6 Apr 2025, Pandit et al., 25 Jan 2026).

1. Hamiltonian Framework and Junction Composition

The junction consists of three regions along the $x$-axis:

• Region I ($x< -d$): A topological nodal-surface semimetal with chemical potential $\mu$. The low-energy Hamiltonian, centered around $k_0 = (0,0,\pi)$, is given by

$$H_{\rm NSSM}(\mathbf q) = q_z(q_x \sigma_x + q_y \sigma_y) + q_z \sigma_z$$

where $\sigma_i$ are Pauli matrices in orbital space, and $\mathbf q = \mathbf k - \mathbf k_0$.

• Region II ($-d < x < 0$): A thin insulating barrier of height $V_0$, width $d$.
• Region III ($x > 0$): An $s$-wave superconductor with pair potential $\Delta\,\Theta(x)$ and chemical potential shifted by $-U_0$ due to gating.

The overall system is described in the Nambu spinor basis $\Psi(x)$, and the BdG Hamiltonian within each region takes the block form

$$\mathcal{H}_{\rm BdG}(x) = \begin{pmatrix} H_{\rm NSSM}(-i\nabla) + U(x) - \mu & \Delta(x) \ \Delta^*(x) & -H_{\rm NSSM}(-i\nabla) - U(x) + \mu \end{pmatrix}$$

with $U(x) = 0, V_0, -U_0$ for $x < -d$, $-d < x < 0$, $x > 0$, respectively (Pandit et al., 25 Jan 2026).

2. Wavefunction Structure and Boundary Matching

Quasiparticle states are constructed in each region:

• NSSM ($x < -d$): Incident, reflected, and Andreev-reflected components are

$$\Psi_{N}(x) = \psi^{e+} e^{i q_x^{e} x} + r_{N}\,\psi^{e-} e^{-i q_x^{e} x} + r_{A}\,\psi^{h-} e^{-i q_x^{h} x}$$

Here, $q_x^e$ and $q_x^h$ are obtained from the NSSM dispersion for electrons and holes, with spinorial structure determined by

$$\chi_{11} = \sqrt{\frac{E + \mu - q_z}{E + \mu + q_z}}\,e^{i\theta_e}, \quad \chi_{22} = \sqrt{\frac{E - \mu + q_z}{E - \mu - q_z}}\,e^{i\theta_A}$$

The angles satisfy a “Snell’s law”: $q_\rho^e \sin\theta_e = q_\rho^h \sin\theta_A$.

• Insulator ($-d < x < 0$): The general state is a superposition of electron- and hole-like waves, with barrier-modified momenta. In the thin-barrier or “delta function” limit ($V_0\to\infty$, $d\to0$ with $V_0 d$ finite), the relevant dimensionless barrier strength is $\beta = V_0 d / q_z$.
• Superconductor ($x > 0$): Bogoliubov quasiparticles with coherence factors

$$u^2 = \frac{1}{2}\left(1 + \frac{\sqrt{E^2 - \Delta^2}}{E}\right),\quad v^2 = \frac{1}{2}\left(1 - \frac{\sqrt{E^2 - \Delta^2}}{E}\right)$$

and phase factors analogous to the normal region.

Boundary conditions require:

• Continuity of the Nambu spinor at $x=-d$ and $x=0$,
• Conservation of quasiparticle flux.

Solving the resulting linear system yields explicit amplitudes $r_N$, $r_A$, $t_e$, $t_h$ for normal and Andreev reflection and transmission (Pandit et al., 25 Jan 2026).

3. Reflection Amplitudes and Conductance Oscillations

The reflection coefficients, for arbitrary incidence, take the form

$$r_N = \frac{A_N + B_N\,e^{2i\beta}}{C + D\,e^{2i\beta}}, \qquad r_A = \frac{A_A}{C + D\,e^{2i\beta}}$$

where $A_N, B_N, A_A, C, D$ are determined by system and excitation parameters.

In the thin-barrier regime, at normal incidence,

$$r_A(\theta_e=0) = \frac{\Delta\,\exp(-i\varphi)}{E + iZ\sqrt{\Delta^2 - E^2}}, \quad Z = \tan\beta$$

Analogous to the Blonder–Tinkham–Klapwijk (BTK) result, but with $Z$ set by the combination $V_0 d / q_z$—reflecting the anisotropic velocity due to the NSSM’s dispersion rather than a simple Fermi velocity as in graphene.

For subgap energy ($E < \Delta$), Andreev and normal reflection coefficients satisfy the unitarity relation

$$|r_N|^2 + \frac{\cos\theta_A}{\cos\theta_e} |r_A|^2 = 1$$

The zero-temperature differential electrical conductance is given by

$$G(eV) = G_0 \int_{-\pi/2}^{+\pi/2} \Bigl[1 - |r_N(E=eV,\theta_e)|^2 + |r_A(E=eV,\theta_e)|^2 \frac{\cos\theta_A}{\cos\theta_e}\Bigr] \cos\theta_e\,d\theta_e$$

In both the high- and low-doping limits, $G(eV)$ and the thermal conductance $\kappa(T)$ oscillate as functions of $\beta$ with period $\pi$:

$G(0) \propto \cos^2\beta$

so the conductance is maximal at $\beta = n\pi$ and vanishes at $\beta = (n+1/2)\pi$. This is in clear contrast to graphene or silicene NIS junctions, where the oscillation period is $\pi/2$ and the maxima/minima positions are shifted (Pandit et al., 25 Jan 2026).

4. Andreev Physics and Floquet Extensions

The NSSM dispersion $E \sim q_z\sqrt{1 + q_\rho^2}$ enables both retro- and specular Andreev reflection depending on the doping regime: retro-AR dominates when $\mu \gg \Delta$, while specular AR appears for $\mu \ll \Delta$, similar to the case in graphene. This duality directly impacts the angular dependence and energy scaling of the Andreev reflection probability $R_a = |r_h|^2$, with steep changes at critical angles or energies determined by mode-matching conditions.

Under high-frequency irradiation (Floquet regime), the stroboscopic Hamiltonian acquires polarization-dependent renormalizations:

• Circular polarization:

$$q_z \mapsto q_z(1 + q_1 q_z), \quad q_1 = \frac{(e E_0)^2}{(\hbar\omega)^3}$$

• Linear polarization (along $y$):

$$q_z \mapsto q_z(1 - q_2 q_z^2), \quad q_2 = \frac{(e E_0)^2}{2(\hbar\omega)^4}$$

These modifications systematically alter the reflection and transmission probabilities, leading to phenomena such as total suppression of subgap AR, enhanced $R_n+R_a \to 1$ for $E > \Delta$ (blocking superconducting quasiparticle transmission), and angular asymmetry in reflectances (Pandit et al., 6 Apr 2025).

5. Comparison to Graphene and Silicene NIS Junctions

The periodicity and phase of conductance oscillations in NSSM-I-SC junctions differ fundamentally from those in conventional normal metal–insulator–superconductor (NIS) junctions on graphene or silicene:

• Periodicity: NSSM-I-SC: $\pi$; graphene/silicene NIS: $\pi/2$.
• Conductance maxima: In NSSM, maxima occur at zero barrier and at integer multiples of $\pi$ in $\beta$; in graphene, the maximum is shifted, and zero barrier yields a minimum.
• Origin: In NSSM, the barrier parameter $Z = \tan\beta$ is set by the nodal-surface momentum $q_z$, introducing a momentum-dependent velocity scale. In graphene, $Z$ is determined by $V_0 d / (\hbar v_F)$, with $v_F$ a constant.

Table: Comparison of Barrier Oscillations

Material	Periodicity	Position of Maximum $G$	Barrier Parameter Origin
NSSM-I-SC	$\pi$	$\beta = n\pi$	$Z = \tan(V_0 d/q_z)$
Graphene NIS	$\pi/2$	$Z = (2n+1)\pi/2$	$Z = V_0 d/(\hbar v_F)$
Silicene NIS	$\pi/2$	Electric-field tunable	As for graphene plus gate tunability

This suggests the doubled period and the distinct phase origin in NSSM junctions constitute direct hallmarks of the underlying band topology (Pandit et al., 25 Jan 2026).

6. Experimental Signatures and Material Considerations

Potential NSSM materials include ZrSiS-type compounds, $\mathrm{Ca}_3\mathrm{P}_2$, and nonsymmorphic oxides (e.g., SrIrO$_3$). Mesoscopic junctions with thin tunnel barriers (e.g., Al$_2$O$_3$ of $d\sim1$ nm and $V_0 \sim$ eV) yield experimentally relevant $\beta$ values. Gating allows independent tuning of $\mu$ (in the NSSM) and $U_0$ (beneath the superconductor).

Key experimental predictions are:

• Electrical conductance ($dI/dV$): Clear $\pi$-periodic oscillations in barrier strength, with maxima at zero barrier and at integer multiples of $\pi$ in $\beta$.
• Thermal conductance ($\kappa(T)$): Similar $\pi$-periodic oscillations, robust across temperature and doping regimes.
• Distinguishing NSSM from graphene: Both the oscillation period and its dependence on $q_z$ (via $\beta$) serve as signatures. Tilting the junction with respect to the nodal plane, modifying $q_z$, directly shifts the period, evidencing the momentum-dependent (anisotropic) dispersion.
• Optoelectronic control: High-frequency irradiation can selectively modulate AR probabilities, enabling photonic gating of Andreev physics (Pandit et al., 6 Apr 2025, Pandit et al., 25 Jan 2026).

A plausible implication is that the NSSM-I-SC platform offers unique routes to control quantum transport with both electrostatic and optical means, informed by the topological features of the underlying semimetal.

7. Relevance and Future Directions

The NSSM-I-SC junction framework provides a comprehensive platform for exploring topologically-enhanced quantum transport, tunable Andreev phenomena, and distinctive barrier-controlled interference effects. The predicted $\pi$-period oscillations in conductance and thermal response, their dependence on doping, barrier geometry, and irradiation, and the sharp contrast with graphene/silicene signatures, establish a roadmap for both fundamental investigations and device-oriented studies in nodal-surface systems (Pandit et al., 6 Apr 2025, Pandit et al., 25 Jan 2026). Continued research may further elucidate the interplay of symmetry, topology, and electron correlations in such hybrid structures, and enable precision tests via advanced spectroscopy and transport experiments in candidate NSSMs.