Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Dipole-Obstructed Cooper Pairing: Theory and Application to j=3/2𝑗32j=3/2italic_j = 3 / 2 Superconductors

Penghao Zhu zhu.3711@osu.edu Department of Physics, The Ohio State University, Columbus, OH 43210, USA    Rui-Xing Zhang ruixing@utk.edu Department of Physics and Astronomy, The University of Tennessee, Knoxville, TN 37996, USA Department of Materials Science and Engineering, The University of Tennessee, Knoxville, TN 37996, USA
Abstract

Like electrons, Cooper pairs can carry a monopole charge if the pairing electrons come from two or more Fermi surfaces with different Chern numbers. In such an instance, a superconductor is necessarily nodal due to an inherent topological pairing obstruction. In this work, we show that a similar obstruction is also possible when there is only one Fermi surface involved in the pairing process. By developing a Chern-vorticity theorem, we have identified a class of Fermi surfaces with a quantized dipolar Berry flux pattern, where all intra-Fermi-surface Cooper pairings are “dipole-obstructed” and nodal. As a real-world application, we find that the dipole obstruction plays a crucial role in stabilizing the superconducting nodal structure for j=3/2𝑗32j=3/2italic_j = 3 / 2 half-Heusler compounds.

Introduction.— Superconductors (SCs) are generally classified by the symmetry pattern of their Cooper pairs. By definition, electron pairing in a conventional s𝑠sitalic_s-wave SC is spatially isotropic, generating a uniform energy gap of the normal-state Fermi surface (FS) [1]. Meanwhile, an unconventional SC such as the cuprates can feature a gapless spectrum when its anisotropic pairing order Δⁱ(đ€)Î”đ€\Delta({\bf k})roman_Δ ( bold_k ) has symmetry-enforced zeros in the momentum space [2, 3, 4]. In experiments, the SC gap structure can be feasibly probed by angle-resolved photoemission spectroscopy [5], scanning tunneling spectroscopy [6], penetration depth measurements [7], etc. Such gap information often offers valuable insights into unraveling the nature of Cooper pairs for a new SC candidate.

The pairing symmetry, however, does not fully account for the gap structure of SCs. For example, a SC is found to be necessarily nodal, once the electrons forming a Cooper pair come from two FSs with opposite Berry monopole charges [8, 9, 10, 11]. Such an inability to develop a full energy gap is intrinsic to the topological texture of the FS, which holds even when Δⁱ(đ€)Î”đ€\Delta({\bf k})roman_Δ ( bold_k ) is constant. This striking phenomenon has sparked a growing research interest in uncovering similar mechanisms of topologically obstructed nodal pairing orders for FSs carrying an Euler index [12, 13, 14] or a â„€2subscriptâ„€2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT index [15]. Notably, most existing theories have assumed the obstructed pairings to involve multiple FSs, while generalizations to intra-FS electron pairings are less explored.

In this work, we have identified a new class of FSs where all intra-FS Cooper pairings exhibit topology-enforced zeros. The FS of our interest encloses a quantized dipole of the Berry curvature [16], i.e., there exists a Berry flux of ±2âąÏ€plus-or-minus2𝜋\pm 2\pi± 2 italic_π through either half of the FS, while the net Berry flux vanishes. As a proof of concept, we consider a minimal model with such a Berry-dipole FS and find all intra-FS Cooper pairings are dipole-obstructed from fully gapping out the FS. Specifically, the resulting SC state always exhibits zero-energy Weyl nodes and/or nodal loops in the Bogoliubov-de Gennes (BdG) spectrum. Finally, we revisit the half-Heusler SCs such as YPtBi [17, 4] and LuPdBi [18, 19, 20], which are believed to feature a mixed-parity singlet-septet pairing [21]. We find that the dipole obstruction naturally exists in the septet pairing channel and further clarify its contribution to the BdG line nodes observed in experiments.

Chern-Vorticity Theorem and Pairing Zeros.— We start by presenting a Chern-vorticity theorem that will guide us to the target FSs. As shown in Fig. 1(a), we consider two 2D closed or effectively closed manifolds ℳ1subscriptℳ1\mathcal{M}_{1}caligraphic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and ℳ2subscriptℳ2\mathcal{M}_{2}caligraphic_M start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT in đ€đ€{\bf k}bold_k-space, as well as two Bloch states |ψ1⁹(đ€1)⟩ketsubscript𝜓1subscriptđ€1|\psi_{1}(\mathbf{k}_{1})\rangle| italic_ψ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ⟩ and |ψ2⁹(đ€2)⟩ketsubscript𝜓2subscriptđ€2|\psi_{2}(\mathbf{k}_{2})\rangle| italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ⟩ with đ€i∈ℳisubscriptđ€đ‘–subscriptℳ𝑖{\bf k}_{i}\in{\cal M}_{i}bold_k start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∈ caligraphic_M start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT. Here, we define an open manifold to be effectively closed if the Bloch states over each of its boundaries are identical. Generally, ℳ1subscriptℳ1\mathcal{M}_{1}caligraphic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and ℳ2subscriptℳ2\mathcal{M}_{2}caligraphic_M start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are related by a đ€đ€\mathbf{k}bold_k-space transformation {g|𝐭}conditional-set𝑔𝐭\{g|{\bf t}\}{ italic_g | bold_t } with đ€2=gâąđ€1+𝐭.subscriptđ€2𝑔subscriptđ€1𝐭\mathbf{k}_{2}=g\mathbf{k}_{1}+\mathbf{t}.bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = italic_g bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + bold_t . Here g𝑔gitalic_g denotes a point-group operation (e.g., rotation and mirror) that satisfies g−1=gT∈ℝsuperscript𝑔1superscript𝑔𝑇ℝg^{-1}=g^{T}\in\mathbb{R}italic_g start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = italic_g start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ∈ blackboard_R and 𝐭𝐭\mathbf{t}bold_t is a translation in đ€đ€\mathbf{k}bold_k-space. We now consider the matrix element of a general two-particle operator O^^𝑂\hat{O}over^ start_ARG italic_O end_ARG bridging electrons on ℳ1subscriptℳ1\mathcal{M}_{1}caligraphic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and ℳ2subscriptℳ2\mathcal{M}_{2}caligraphic_M start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT:

đ’Ș⁹(đ€1)=⟚ψ1⁹(đ€1)|O^|ψ2⁹(đ€2)⟩=|đ’Ș⁹(đ€1)|⁹eiâąÏ†âą(đ€1).đ’Șsubscriptđ€1quantum-operator-productsubscript𝜓1subscriptđ€1^𝑂subscript𝜓2subscriptđ€2đ’Șsubscriptđ€1superscript𝑒𝑖𝜑subscriptđ€1{\cal O}({\bf k}_{1})=\langle\psi_{1}(\mathbf{k}_{1})|\hat{O}|\psi_{2}(\mathbf% {k}_{2})\rangle=|{\cal O}({\bf k}_{1})|e^{i\varphi({\bf k}_{1})}.caligraphic_O ( bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = ⟹ italic_ψ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | over^ start_ARG italic_O end_ARG | italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ⟩ = | caligraphic_O ( bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | italic_e start_POSTSUPERSCRIPT italic_i italic_φ ( bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_POSTSUPERSCRIPT . (1)

When Ï†đœ‘\varphiitalic_φ displays a vortex pattern around đ€1=đ€visubscriptđ€1subscriptđ€subscript𝑣𝑖{\bf k}_{1}={\bf k}_{v_{i}}bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = bold_k start_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT, the value of đ’Ș⁹(đ€vi)đ’Șsubscriptđ€subscript𝑣𝑖{\cal O}({\bf k}_{v_{i}})caligraphic_O ( bold_k start_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) necessarily vanishes and is thus topologically obstructed. Such a vortex visubscript𝑣𝑖v_{i}italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is captured by a vorticity index around đ€visubscriptđ€subscript𝑣𝑖{\bf k}_{v_{i}}bold_k start_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT, defined as Îœi=1/(2âąÏ€)âąâˆźđ‘‘đ€1â‹…âˆ‚đ€1Ï†âą(đ€1)∈℀subscript𝜈𝑖12𝜋contour-integral⋅differential-dsubscriptđ€1subscriptsubscriptđ€1𝜑subscriptđ€1â„€\nu_{i}=1/(2\pi)\oint d\mathbf{k}_{1}\cdot\partial_{\mathbf{k}_{1}}\varphi({% \bf k}_{1})\in\mathbb{Z}italic_Îœ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = 1 / ( 2 italic_π ) ∟ italic_d bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ⋅ ∂ start_POSTSUBSCRIPT bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_φ ( bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ∈ blackboard_Z. The Chern-vorticity theorem proven in the Supplemental Material (SM) [22] dictates that the net vorticity on ℳ1subscriptℳ1{\cal M}_{1}caligraphic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, Îœ=∑iÎœiâˆˆâ„€đœˆsubscript𝑖subscriptđœˆđ‘–â„€\nu=\sum_{i}\nu_{i}\in\mathbb{Z}italic_Îœ = ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_Îœ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∈ blackboard_Z, is determined by

Îœđœˆ\displaystyle\nuitalic_Îœ =\displaystyle== 𝒞1−(det⁥g)ⁱ𝒞2+ℐφ,subscript𝒞1det𝑔subscript𝒞2subscriptℐ𝜑\displaystyle{\cal C}_{1}-({\operatorname{det}g}){\cal C}_{2}+{\cal I_{\varphi% }},caligraphic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - ( roman_det italic_g ) caligraphic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + caligraphic_I start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ,
ℐφsubscriptℐ𝜑\displaystyle{\cal I}_{\varphi}caligraphic_I start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT =\displaystyle== ∑∂ℳ1∼∂ℳ1dâąđ€12âąÏ€â‹…âˆ‚đ€1Ï†âą(đ€1),subscriptsubscriptℳ1subscriptcontour-integralsubscriptℳ1⋅𝑑subscriptđ€12𝜋subscriptsubscriptđ€1𝜑subscriptđ€1\displaystyle\sum_{\partial{\cal M}_{1}}\oint_{\partial{\cal M}_{1}}\frac{d% \mathbf{k}_{1}}{2\pi}\cdot\partial_{\mathbf{k}_{1}}\varphi(\mathbf{k}_{1}),∑ start_POSTSUBSCRIPT ∂ caligraphic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ∟ start_POSTSUBSCRIPT ∂ caligraphic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT divide start_ARG italic_d bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_π end_ARG ⋅ ∂ start_POSTSUBSCRIPT bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_φ ( bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) , (2)

where đ’žÎ±subscriptđ’žđ›Œ{\cal C}_{\alpha}caligraphic_C start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT is the single-particle Chern number of |uα⁹(đ€Î±)âŸ©â‰Ąe−iâąđ€Î±â‹…đ«Î±âą|ψ2⁹(đ€Î±)⟩ketsubscriptđ‘ąđ›Œsubscriptđ€đ›Œsuperscript𝑒⋅𝑖subscriptđ€đ›Œsubscriptđ«đ›Œketsubscript𝜓2subscriptđ€đ›Œ|u_{\alpha}(\mathbf{k}_{\alpha})\rangle\equiv e^{-i{\bf k}_{\alpha}\cdot{\bf r% }_{\alpha}}|\psi_{2}(\mathbf{k}_{\alpha})\rangle| italic_u start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT ) ⟩ ≡ italic_e start_POSTSUPERSCRIPT - italic_i bold_k start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT ⋅ bold_r start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT end_POSTSUPERSCRIPT | italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT ) ⟩ on ℳαsubscriptâ„łđ›Œ{\cal M}_{\alpha}caligraphic_M start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT. The loop integral ℐφsubscriptℐ𝜑{\cal I}_{\varphi}caligraphic_I start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT is summed over all possible boundaries of ℳ1subscriptℳ1{\cal M}_{1}caligraphic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, which, by default, vanishes for a closed ℳ1subscriptℳ1{\cal M}_{1}caligraphic_M start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT.

Focusing on SCs, we choose |ψ1⁹(đ€1)⟩=|Ο1(e)⁹(đ€1)⟩ketsubscript𝜓1subscriptđ€1ketsubscriptsuperscript𝜉𝑒1subscriptđ€1|\psi_{1}(\mathbf{k}_{1})\rangle=|\xi^{(e)}_{1}(\mathbf{k}_{1})\rangle| italic_ψ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ⟩ = | italic_Ο start_POSTSUPERSCRIPT ( italic_e ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ⟩ and |ψ2⁹(đ€2)⟩=|Ο2(e)⁣⋆ⁱ(đ€2)⟩ketsubscript𝜓2subscriptđ€2ketsubscriptsuperscript𝜉𝑒⋆2subscriptđ€2|\psi_{2}(\mathbf{k}_{2})\rangle=|\xi^{(e)\star}_{2}(\mathbf{k}_{2})\rangle| italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ⟩ = | italic_Ο start_POSTSUPERSCRIPT ( italic_e ) ⋆ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ⟩ in order to properly project the pairing operator O^^𝑂\hat{O}over^ start_ARG italic_O end_ARG, since O^^𝑂\hat{O}over^ start_ARG italic_O end_ARG physically describes hoppings between electrons and holes. Here, |Οα(e)⁹(đ€Î±)⟩ketsubscriptsuperscriptđœ‰đ‘’đ›Œsubscriptđ€đ›Œ|\xi^{(e)}_{\alpha}(\mathbf{k}_{\alpha})\rangle| italic_Ο start_POSTSUPERSCRIPT ( italic_e ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT ) ⟩ is the Bloch state of electrons on the Fermi surface ℳαsubscriptâ„łđ›Œ{\cal M}_{\alpha}caligraphic_M start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT. It is crucial to note that 𝒞2subscript𝒞2{\cal C}_{2}caligraphic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, the Chern number for |ψ2⁹(đ€2)⟩ketsubscript𝜓2subscriptđ€2|\psi_{2}(\mathbf{k}_{2})\rangle| italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ⟩, is exactly opposite to 𝒞2(e)superscriptsubscript𝒞2𝑒{\cal C}_{2}^{(e)}caligraphic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_e ) end_POSTSUPERSCRIPT, the electronic Chern number for |Ο2(e)⁹(đ€2)⟩ketsubscriptsuperscript𝜉𝑒2subscriptđ€2|\xi^{(e)}_{2}(\mathbf{k}_{2})\rangle| italic_Ο start_POSTSUPERSCRIPT ( italic_e ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ⟩ on ℳ2subscriptℳ2{\cal M}_{2}caligraphic_M start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT [22]. As a concrete example, we immediately arrive at Îœ=2𝜈2\nu=2italic_Îœ = 2 when (i) ℳ1,2subscriptℳ12\mathcal{M}_{1,2}caligraphic_M start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT are both closed FSs with 𝒞1=𝒞2=+1subscript𝒞1subscript𝒞21{\cal C}_{1}={\cal C}_{2}=+1caligraphic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = caligraphic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = + 1 and (ii) g𝑔gitalic_g is the spatial inversion with det⁥g=−1det𝑔1{\operatorname{det}g}=-1roman_det italic_g = - 1, which exactly corresponds to the obstructed nodal pairing for monopole SCs in Ref. [9]. Therefore, the Chern-vorticity theorem manifests as a natural generalization of Ref. [9], but its application transcends the realm of SC systems with inter-FS pairings [23, 24, 25].

Refer to caption
Figure 1: (a) The Chern-vorticity theorem informs the phase vortex of O^^𝑂\hat{O}over^ start_ARG italic_O end_ARG when projecting onto two closed or effectively closed manifolds ℳ1,2subscriptℳ12{\cal M}_{1,2}caligraphic_M start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT. (b) Berry-dipole FSs for h0subscriptℎ0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, where the arrows denote the Berry curvature vectors on the FS. A change of ÎŒđœ‡\muitalic_ÎŒ or ÎŁÎŁ\Sigmaroman_ÎŁ can modify the FS topology, while the quantization of Berry dipole always remains robust.

Berry-Dipole Fermi Surface.— The above Chern-vorticity theorem suggests that for a single Fermi surface ℳ≡ℳ+âˆȘℳ−ℳsubscriptℳsubscriptℳ{\cal M}\equiv{\cal M}_{+}\cup{\cal M}_{-}caligraphic_M ≡ caligraphic_M start_POSTSUBSCRIPT + end_POSTSUBSCRIPT âˆȘ caligraphic_M start_POSTSUBSCRIPT - end_POSTSUBSCRIPT, the intra-FS pairing between two patches ℳ±subscriptℳplus-or-minus{\cal M}_{\pm}caligraphic_M start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT can be topologically obstructed if each patch is effectively closed and carries a nonzero quantized Berry flux. Further requiring the net Berry flux to be zero, we are thus looking for a FS with a quantized dipolar texture of Berry curvature, i.e., a Berry-dipole FS [26, 16, 27, 28, 29].

As discussed in Refs. [26, 16] and reviewed in the SM [22], quantization of a Berry-dipole can be achieved by enforcing a mirror symmetry along the dipole axis. As a concrete example, we consider a minimal model with Berry-dipole physics proposed in Ref. [16],

h0=2⁹kz⁹(kxâąÏƒx+kyâąÏƒy)+(kx2+ky2−kz2+ÎŁ)âąÏƒz.subscriptℎ02subscript𝑘𝑧subscriptđ‘˜đ‘„subscriptđœŽđ‘„subscript𝑘𝑩subscript𝜎𝑩superscriptsubscriptđ‘˜đ‘„2superscriptsubscript𝑘𝑩2superscriptsubscript𝑘𝑧2ÎŁsubscript𝜎𝑧h_{0}=2k_{z}(k_{x}\sigma_{x}+k_{y}\sigma_{y})+(k_{x}^{2}+k_{y}^{2}-k_{z}^{2}+% \Sigma)\sigma_{z}.italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) + ( italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_ÎŁ ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT . (3)

Besides mirror symmetry along z𝑧zitalic_z-direction Mz=iâąÏƒzsubscript𝑀𝑧𝑖subscript𝜎𝑧M_{z}=i\sigma_{z}italic_M start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = italic_i italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT, h0subscriptℎ0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT respects a continuous rotation symmetry around z𝑧zitalic_z with CΞ=exp⁹[−i⁹Jz⁹Ξ]subscriptđ¶đœƒexpdelimited-[]𝑖subscriptđœđ‘§đœƒC_{\theta}=\text{exp}[-iJ_{z}\theta]italic_C start_POSTSUBSCRIPT italic_Ξ end_POSTSUBSCRIPT = exp [ - italic_i italic_J start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_Ξ ], where Jz=diag⁹(3/2,1/2)subscriptđœđ‘§diag3212J_{z}=\text{diag}(3/2,1/2)italic_J start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = diag ( 3 / 2 , 1 / 2 ) is a diagonal matrix. Regularizing h0subscriptℎ0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT will reduce the symmetry group to C4⁹hsubscriptđ¶4ℎC_{4h}italic_C start_POSTSUBSCRIPT 4 italic_h end_POSTSUBSCRIPT, which will be exploited for later pairing analysis. When ÎŁ>0ÎŁ0\Sigma>0roman_ÎŁ > 0, h0subscriptℎ0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT describes a 3D minimal Weyl semimetal (WSM) with two Weyl nodes at kz=±Σsubscript𝑘𝑧plus-or-minusÎŁk_{z}=\pm\sqrt{\Sigma}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = ± square-root start_ARG roman_ÎŁ end_ARG, respectively. At ÎŁ=0ÎŁ0\Sigma=0roman_ÎŁ = 0, the Weyl nodes merge at ΓΓ\Gammaroman_Γ to form a quadratic band touching. The annihilation of Weyl nodes, however, does not lead to a mass generation for the low-energy electrons. Instead, we note that h0subscriptℎ0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT with ÎŁ<0ÎŁ0\Sigma<0roman_ÎŁ < 0 features a 1D doubly degenerate nodal loop in the kz=0subscript𝑘𝑧0k_{z}=0italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0 plane, i.e., a nodal-loop semimetal (NLSM). This unexpected robustness of Weyl nodes arises from the Berry-dipole delicate topological charges [26, 16].

Let us generally denote the FS patch with kz>0subscript𝑘𝑧0k_{z}>0italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT > 0 (kz<0subscript𝑘𝑧0k_{z}<0italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT < 0) as ℳ+subscriptℳ\mathcal{M}_{+}caligraphic_M start_POSTSUBSCRIPT + end_POSTSUBSCRIPT (ℳ−subscriptℳ\mathcal{M}_{-}caligraphic_M start_POSTSUBSCRIPT - end_POSTSUBSCRIPT). For |ÎŒ|>|ÎŁ|𝜇Σ|\mu|>|\Sigma|| italic_ÎŒ | > | roman_ÎŁ |, we always find a single spherical-like FS ℳℳ{\cal M}caligraphic_M regardless of the value of ÎŁÎŁ\Sigmaroman_ÎŁ, where ℳ+subscriptℳ\mathcal{M}_{+}caligraphic_M start_POSTSUBSCRIPT + end_POSTSUBSCRIPT and ℳ−subscriptℳ\mathcal{M}_{-}caligraphic_M start_POSTSUBSCRIPT - end_POSTSUBSCRIPT are hence the north and south hemispheres, respectively. Undergoing the Lifshitz transition with |ÎŒ|<|ÎŁ|𝜇Σ|\mu|<|\Sigma|| italic_ÎŒ | < | roman_ÎŁ |, the WSM with ÎŁ>0ÎŁ0\Sigma>0roman_ÎŁ > 0 comprises a pair of closed spherical FSs (i.e., ℳ±subscriptℳplus-or-minus{\cal M}_{\pm}caligraphic_M start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT), while a single torus-like FS is found for the NLSM (ÎŁ<0ÎŁ0\Sigma<0roman_ÎŁ < 0) with ℳ±subscriptℳplus-or-minus\mathcal{M}_{\pm}caligraphic_M start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT now being both open and of an annulus shape. Evidently, the boundaries of ℳ+subscriptℳ\mathcal{M}_{+}caligraphic_M start_POSTSUBSCRIPT + end_POSTSUBSCRIPT and ℳ−subscriptℳ\mathcal{M}_{-}caligraphic_M start_POSTSUBSCRIPT - end_POSTSUBSCRIPT, if present, are always sitting in the kz=0subscript𝑘𝑧0k_{z}=0italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0 plane regardless of ÎŁÎŁ\Sigmaroman_ÎŁ. Consequently, the eigenstates of h0subscriptℎ0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT along each boundary can always be uniform, owing to the Mzsubscript𝑀𝑧M_{z}italic_M start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT symmetry. Hence, we conclude that ℳ+subscriptℳ\mathcal{M}_{+}caligraphic_M start_POSTSUBSCRIPT + end_POSTSUBSCRIPT and ℳ−subscriptℳ\mathcal{M}_{-}caligraphic_M start_POSTSUBSCRIPT - end_POSTSUBSCRIPT must be either closed or effectively closed for all choices of ÎŒđœ‡\muitalic_ÎŒ and ÎŁÎŁ\Sigmaroman_ÎŁ. This guarantees the FS(s) of h0subscriptℎ0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT to always feature a quantized Berry dipole, as explicitly confirmed in Fig. 1(b).

Dipole-Obstructed Pairing.— We are now ready to explore Cooper pairing physics on the above Berry-dipole FS by updating h0subscriptℎ0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT to a BdG Hamiltonian,

H⁹(đ€)=(h0⁹(đ€)âˆ’ÎŒÎ”âą(đ€)Δ†ⁱ(đ€)Ό−h0T⁹(âˆ’đ€)).đ»đ€matrixsubscriptℎ0đ€đœ‡Î”đ€superscriptÎ”â€ đ€đœ‡subscriptsuperscriptℎ𝑇0đ€H(\mathbf{k})=\begin{pmatrix}h_{0}(\mathbf{k})-\mu&\Delta(\mathbf{k})\\ \Delta^{{\dagger}}(\mathbf{k})&\mu-h^{T}_{0}(-\mathbf{k})\end{pmatrix}.italic_H ( bold_k ) = ( start_ARG start_ROW start_CELL italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( bold_k ) - italic_ÎŒ end_CELL start_CELL roman_Δ ( bold_k ) end_CELL end_ROW start_ROW start_CELL roman_Δ start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( bold_k ) end_CELL start_CELL italic_ÎŒ - italic_h start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( - bold_k ) end_CELL end_ROW end_ARG ) . (4)

The pairing matrix generally takes the form Δⁱ(đ€)=d0⁹(đ€)âąÏƒ0+𝐝ⁱ(đ€)â‹…đˆÎ”đ€subscript𝑑0đ€subscript𝜎0â‹…đđ€đˆ\Delta(\mathbf{k})=d_{0}(\mathbf{k})\sigma_{0}+\mathbf{d}(\mathbf{k})\cdot\bm{\sigma}roman_Δ ( bold_k ) = italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( bold_k ) italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + bold_d ( bold_k ) ⋅ bold_italic_σ and the Fermi statistics requires Δⁱ(đ€)=−ΔT⁹(âˆ’đ€)Î”đ€superscriptÎ”đ‘‡đ€\Delta({\bf k})=-\Delta^{T}(-{\bf k})roman_Δ ( bold_k ) = - roman_Δ start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( - bold_k ). For our purpose, we consider expanding di⁹(đ€)subscriptđ‘‘đ‘–đ€d_{i}({\bf k})italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( bold_k ) up to đ’Ș⁹(đ€2)đ’Șsuperscriptđ€2{\cal O}({\bf k}^{2})caligraphic_O ( bold_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) and classify all pairing channels based on the irreducible representations (irreps) of C4⁹hsubscriptđ¶4ℎC_{4h}italic_C start_POSTSUBSCRIPT 4 italic_h end_POSTSUBSCRIPT group in Table. 1 [22]. For |ÎŒ|>|ÎŁ|𝜇Σ|\mu|>|\Sigma|| italic_ÎŒ | > | roman_ÎŁ |, Δⁱ(đ€)Î”đ€\Delta({\bf k})roman_Δ ( bold_k ) describes the intra-FS pairings (as there is only one FS). Following Eq. 1, the projected Cooper pairing onto ℳ±subscriptℳplus-or-minus{\cal M}_{\pm}caligraphic_M start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT is given by Δeff=⟚Ο+(e)⁹(đ€)|Δⁱ(đ€)|Ο−(e)⁣⋆ⁱ(âˆ’đ€)⟩subscriptΔeffquantum-operator-productsubscriptsuperscriptđœ‰đ‘’đ€Î”đ€subscriptsuperscriptđœ‰đ‘’â‹†đ€\Delta_{\text{eff}}=\langle\xi^{(e)}_{+}(\mathbf{k})|\Delta(\mathbf{k})|\xi^{(% e)\star}_{-}(-\mathbf{k})\rangleroman_Δ start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT = ⟹ italic_Ο start_POSTSUPERSCRIPT ( italic_e ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( bold_k ) | roman_Δ ( bold_k ) | italic_Ο start_POSTSUPERSCRIPT ( italic_e ) ⋆ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ( - bold_k ) ⟩, where |Ο±(e)⟩ketsubscriptsuperscript𝜉𝑒plus-or-minus|\xi^{(e)}_{\pm}\rangle| italic_Ο start_POSTSUPERSCRIPT ( italic_e ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ⟩ are normal states over ℳ±subscriptℳplus-or-minus\mathcal{M}_{\pm}caligraphic_M start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT.

We first note that a general pairing with ℐφ=0subscriptℐ𝜑0{\cal I}_{\varphi}=0caligraphic_I start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT = 0 must be topologically obstructed since Îœ=C1+C2=2⁹sgn⁹(ÎŒ)𝜈subscriptđ¶1subscriptđ¶22sgn𝜇\nu=C_{1}+C_{2}=2\text{sgn}(\mu)italic_Îœ = italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 2 sgn ( italic_ÎŒ ), which directly applies to all σx,ysubscriptđœŽđ‘„đ‘Š\sigma_{x,y}italic_σ start_POSTSUBSCRIPT italic_x , italic_y end_POSTSUBSCRIPT pairings. To see this, note that |Ο±(e)⁹(kz=0)⟩ketsuperscriptsubscript𝜉plus-or-minus𝑒subscript𝑘𝑧0|\xi_{\pm}^{(e)}(k_{z}=0)\rangle| italic_Ο start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_e ) end_POSTSUPERSCRIPT ( italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0 ) ⟩ must be eigenstates of σzsubscript𝜎𝑧\sigma_{z}italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT, since h0|kz=0∌σzsimilar-toevaluated-atsubscriptℎ0subscript𝑘𝑧0subscript𝜎𝑧h_{0}|_{k_{z}=0}\sim\sigma_{z}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0 end_POSTSUBSCRIPT ∌ italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT. As a result, any Δⁱ(đ€)∌σx,ysimilar-toÎ”đ€subscriptđœŽđ‘„đ‘Š\Delta({\bf k})\sim\sigma_{x,y}roman_Δ ( bold_k ) ∌ italic_σ start_POSTSUBSCRIPT italic_x , italic_y end_POSTSUBSCRIPT will vanish on the equator upon projection, further leading to ℐφ=0subscriptℐ𝜑0{\cal I}_{\varphi}=0caligraphic_I start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT = 0. Meanwhile, the σ±subscript𝜎plus-or-minus\sigma_{\pm}italic_σ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT component of Δⁱ(đ€)Î”đ€\Delta({\bf k})roman_Δ ( bold_k ) can contribute to ℐφsubscriptℐ𝜑{\cal I}_{\varphi}caligraphic_I start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT by the winding phase of d+subscript𝑑d_{+}italic_d start_POSTSUBSCRIPT + end_POSTSUBSCRIPT (d−subscript𝑑d_{-}italic_d start_POSTSUBSCRIPT - end_POSTSUBSCRIPT) around the equator for ÎŒ>0𝜇0\mu>0italic_ÎŒ > 0 (ÎŒ<0𝜇0\mu<0italic_ÎŒ < 0), so long as d±⁹(kz=0)≠0subscript𝑑plus-or-minussubscript𝑘𝑧00d_{\pm}(k_{z}=0)\neq 0italic_d start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0 ) ≠ 0. Here, we have defined σ±=(σ0±σz)/2subscript𝜎plus-or-minusplus-or-minussubscript𝜎0subscript𝜎𝑧2\sigma_{\pm}=(\sigma_{0}\pm\sigma_{z})/2italic_σ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT = ( italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ± italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ) / 2 and d±⁹(đ€)=d0⁹(đ€)±dz⁹(đ€)subscript𝑑plus-or-minusđ€plus-or-minussubscript𝑑0đ€subscriptđ‘‘đ‘§đ€d_{\pm}({\bf k})=d_{0}({\bf k})\pm d_{z}({\bf k})italic_d start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( bold_k ) = italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( bold_k ) ± italic_d start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( bold_k ). Based on Table. 1, it is straightforward to see that Ausubscript𝐮𝑱A_{u}italic_A start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT and Busubscriptđ”đ‘ąB_{u}italic_B start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT pairings are the only two pairing channels that can contribute to ℐφsubscriptℐ𝜑{\cal I}_{\varphi}caligraphic_I start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT. Since d±⁹(đ€)subscript𝑑plus-or-minusđ€d_{\pm}({\bf k})italic_d start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( bold_k ) must be an odd function of k±=kx±i⁹kysubscript𝑘plus-or-minusplus-or-minussubscriptđ‘˜đ‘„đ‘–subscript𝑘𝑩k_{\pm}=k_{x}\pm ik_{y}italic_k start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT = italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ± italic_i italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT for both Ausubscript𝐮𝑱A_{u}italic_A start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT and Busubscriptđ”đ‘ąB_{u}italic_B start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT irreps, the corresponding phase winding of d±subscript𝑑plus-or-minusd_{\pm}italic_d start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT around the equator must be (4⁹m+1)âąÏ€4𝑚1𝜋(4m+1)\pi( 4 italic_m + 1 ) italic_π with mâˆˆâ„€đ‘šâ„€m\in\mathbb{Z}italic_m ∈ blackboard_Z [22]. This necessarily leads to a non-zero odd-integer-valued Îœđœˆ\nuitalic_Îœ for all relevant Ausubscript𝐮𝑱A_{u}italic_A start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT and Busubscriptđ”đ‘ąB_{u}italic_B start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT pairings. Since all other pairings feature ℐφ=0subscriptℐ𝜑0{\cal I}_{\varphi}=0caligraphic_I start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT = 0 and Îœ=2𝜈2\nu=2italic_Îœ = 2, we conclude that all intra-FS pairing channels for the Berry-dipole FS are obstructed and nodal.

Agsubscript𝐮𝑔A_{g}italic_A start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT Bgsubscriptđ”đ‘”B_{g}italic_B start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT Eg1superscriptsubscript𝐾𝑔1{}^{1}E_{g}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT italic_E start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT Eg2superscriptsubscript𝐾𝑔2{}^{2}E_{g}start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_E start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT Ausubscript𝐮𝑱A_{u}italic_A start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT Busubscriptđ”đ‘ąB_{u}italic_B start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT Eu1superscriptsubscript𝐾𝑱1{}^{1}E_{u}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT italic_E start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT Eu2superscriptsubscript𝐾𝑱2{}^{2}E_{u}start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_E start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT
(kx2−ky2)âąÏƒysuperscriptsubscriptđ‘˜đ‘„2superscriptsubscript𝑘𝑩2subscript𝜎𝑩(k_{x}^{2}-k_{y}^{2})\sigma_{y}( italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT Δ0âąÏƒysubscriptΔ0subscript𝜎𝑩\Delta_{0}\sigma_{y}roman_Δ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT kz⁹kâˆ’âąÏƒysubscript𝑘𝑧subscript𝑘subscript𝜎𝑩k_{z}k_{-}\sigma_{y}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT kz⁹k+âąÏƒysubscript𝑘𝑧subscript𝑘subscript𝜎𝑩k_{z}k_{+}\sigma_{y}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT kÂ±âąÏƒÂ±subscript𝑘plus-or-minussubscript𝜎plus-or-minusk_{\pm}\sigma_{\pm}italic_k start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT kâˆ“âąÏƒÂ±subscript𝑘minus-or-plussubscript𝜎plus-or-minusk_{\mp}\sigma_{\pm}italic_k start_POSTSUBSCRIPT ∓ end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT kzâąÏƒâˆ’subscript𝑘𝑧subscript𝜎k_{z}\sigma_{-}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT - end_POSTSUBSCRIPT kzâąÏƒ+subscript𝑘𝑧subscript𝜎k_{z}\sigma_{+}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT
kx⁹kyâąÏƒysubscriptđ‘˜đ‘„subscript𝑘𝑩subscript𝜎𝑩k_{x}k_{y}\sigma_{y}italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT kzâąÏƒxsubscript𝑘𝑧subscriptđœŽđ‘„k_{z}\sigma_{x}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT kâˆ’âąÏƒxsubscript𝑘subscriptđœŽđ‘„k_{-}\sigma_{x}italic_k start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT k+âąÏƒxsubscript𝑘subscriptđœŽđ‘„k_{+}\sigma_{x}italic_k start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT
Table 1: Pairing classification for Berry-dipole FS following the C4⁹hsubscriptđ¶4ℎC_{4h}italic_C start_POSTSUBSCRIPT 4 italic_h end_POSTSUBSCRIPT irreps. For Bgsubscriptđ”đ‘”B_{g}italic_B start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT-pairing, dy=Δ0subscript𝑑𝑩subscriptΔ0d_{y}=\Delta_{0}italic_d start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = roman_Δ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is a constant.

As a concrete example, let us focus on Δ=dx⁹(đ€)âąÏƒxΔsubscriptđ‘‘đ‘„đ€subscriptđœŽđ‘„\Delta=d_{x}({\bf k})\sigma_{x}roman_Δ = italic_d start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ( bold_k ) italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT to explicitly illustrate the dipole-induced obstruction and relegate the discussions on other pairing channels in the SM [22]. Without loss of generality, we set ÎŁ=0ÎŁ0\Sigma=0roman_ÎŁ = 0 and the FS at Ό≠0𝜇0\mu\neq 0italic_ÎŒ ≠ 0 is a sphere of radius |ÎŒ|𝜇\sqrt{|\mu|}square-root start_ARG | italic_ÎŒ | end_ARG, which can be parameterized by a polar angle Ξ∈[0,π]𝜃0𝜋\theta\in[0,\pi]italic_Ξ ∈ [ 0 , italic_π ] and an azimuthal angle ϕ∈[0,2âąÏ€)italic-ϕ02𝜋\phi\in[0,2\pi)italic_ϕ ∈ [ 0 , 2 italic_π ). For ÎŒ<0𝜇0\mu<0italic_ÎŒ < 0, we find the electron wavefunction to be |Ο(e,I)⁹(đ€)⟩=(−cos⁥Ξ⁹e−iâąÏ•,sin⁥Ξ)Tketsuperscriptđœ‰đ‘’đŒđ€superscript𝜃superscript𝑒𝑖italic-Ï•đœƒđ‘‡|\xi^{(e,I)}({\bf k})\rangle=(-\cos\theta e^{-i\phi},\sin\theta)^{T}| italic_Ο start_POSTSUPERSCRIPT ( italic_e , italic_I ) end_POSTSUPERSCRIPT ( bold_k ) ⟩ = ( - roman_cos italic_Ξ italic_e start_POSTSUPERSCRIPT - italic_i italic_ϕ end_POSTSUPERSCRIPT , roman_sin italic_Ξ ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT, which is a constant spinor (0,1)Tsuperscript01𝑇(0,1)^{T}( 0 , 1 ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT at the equator. Under this gauge choice, ℳ±subscriptℳplus-or-minus{\cal M}_{\pm}caligraphic_M start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT are both effectively closed, while |Ο(e,I)⁹(đ€)⟩ketsuperscriptđœ‰đ‘’đŒđ€|\xi^{(e,I)}({\bf k})\rangle| italic_Ο start_POSTSUPERSCRIPT ( italic_e , italic_I ) end_POSTSUPERSCRIPT ( bold_k ) ⟩ becomes singular at the north pole (Ξ=0𝜃0\theta=0italic_Ξ = 0). This implies an obstruction to define a globally smooth gauge, thanks to the non-zero Chern number on each hemisphere.

To ensure the wavefunction is merely locally singular, we consider |Ο(e,I)⁹(đ€)⟩ketsuperscriptđœ‰đ‘’đŒđ€|\xi^{(e,I)}({\bf k})\rangle| italic_Ο start_POSTSUPERSCRIPT ( italic_e , italic_I ) end_POSTSUPERSCRIPT ( bold_k ) ⟩ for Ξ∈[π4,3âąÏ€4]𝜃𝜋43𝜋4\theta\in[\frac{\pi}{4},\frac{3\pi}{4}]italic_Ξ ∈ [ divide start_ARG italic_π end_ARG start_ARG 4 end_ARG , divide start_ARG 3 italic_π end_ARG start_ARG 4 end_ARG ], and choose a different gauge choice for Ξ∈[0,π4]âˆȘ[3âąÏ€4,π]𝜃0𝜋43𝜋4𝜋\theta\in[0,\frac{\pi}{4}]\cup[\frac{3\pi}{4},\pi]italic_Ξ ∈ [ 0 , divide start_ARG italic_π end_ARG start_ARG 4 end_ARG ] âˆȘ [ divide start_ARG 3 italic_π end_ARG start_ARG 4 end_ARG , italic_π ] with |Ο(e,I⁹I)⁹(đ€)⟩=(−cos⁥Ξ,sin⁥Ξ⁹eiâąÏ•)Tketsuperscriptđœ‰đ‘’đŒđŒđ€superscript𝜃𝜃superscript𝑒𝑖italic-Ï•đ‘‡|\xi^{(e,II)}({\bf k})\rangle=(-\cos\theta,\sin\theta e^{i\phi})^{T}| italic_Ο start_POSTSUPERSCRIPT ( italic_e , italic_I italic_I ) end_POSTSUPERSCRIPT ( bold_k ) ⟩ = ( - roman_cos italic_Ξ , roman_sin italic_Ξ italic_e start_POSTSUPERSCRIPT italic_i italic_ϕ end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT. Straightforward calculations lead to:

Δeff(I)=dx⁹eiâąÏ•âąsin⁥2⁹Ξ,Δeff(I⁹I)=dx⁹e−iâąÏ•âąsin⁥2⁹Ξ.formulae-sequencesuperscriptsubscriptΔeffđŒsubscriptđ‘‘đ‘„superscript𝑒𝑖italic-ϕ2𝜃superscriptsubscriptΔeffđŒđŒsubscriptđ‘‘đ‘„superscript𝑒𝑖italic-ϕ2𝜃\Delta_{\text{eff}}^{(I)}=d_{x}e^{i\phi}\sin 2\theta,\ \ \Delta_{\text{eff}}^{% (II)}=d_{x}e^{-i\phi}\sin 2\theta.roman_Δ start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_I ) end_POSTSUPERSCRIPT = italic_d start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ϕ end_POSTSUPERSCRIPT roman_sin 2 italic_Ξ , roman_Δ start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_I italic_I ) end_POSTSUPERSCRIPT = italic_d start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_ϕ end_POSTSUPERSCRIPT roman_sin 2 italic_Ξ . (5)

Apparently, ΔeffsubscriptΔeff\Delta_{\text{eff}}roman_Δ start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT features two zeros on ℳ+subscriptℳ{\cal M}_{+}caligraphic_M start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, one at the north pole and another at the equator. Based on the form of Δeff(I⁹I)superscriptsubscriptΔeffđŒđŒ\Delta_{\text{eff}}^{(II)}roman_Δ start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_I italic_I ) end_POSTSUPERSCRIPT, it is easy to see that the vorticity of the north-pole zero is Îœ0=−1subscript𝜈01\nu_{0}=-1italic_Îœ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = - 1. Meanwhile, the vorticity for the equator zero is also −11-1- 1, which can be achieved by performing a clockwise loop integral of the phase of Δeff(I)superscriptsubscriptΔeffđŒ\Delta_{\text{eff}}^{(I)}roman_Δ start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_I ) end_POSTSUPERSCRIPT at Ξ=π2âˆ’Ï”đœƒđœ‹2italic-Ï”\theta=\frac{\pi}{2}-\epsilonitalic_Ξ = divide start_ARG italic_π end_ARG start_ARG 2 end_ARG - italic_ϔ [22]. Together, we find the net vorticity Îœ=−2𝜈2\nu=-2italic_Îœ = - 2 for ℳ+subscriptℳ{\cal M}_{+}caligraphic_M start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, which is consistent with the prediction of the Chern-vorticity theorem.

We now make a few remarks. First of all, the vorticity-induced pairing zeros always exist, regardless of the detailed form of dx⁹(đ€)subscriptđ‘‘đ‘„đ€d_{x}({\bf k})italic_d start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ( bold_k ). The zeros of dx⁹(đ€)subscriptđ‘‘đ‘„đ€d_{x}({\bf k})italic_d start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ( bold_k ) itself may lead to additional zeros of ΔeffsubscriptΔeff\Delta_{\text{eff}}roman_Δ start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT, beyond the topologically obstructed ones. Second, we highlight that the counting of Îœđœˆ\nuitalic_Îœ for intra-FS pairing does depend on a special gauge choice with which ℳ±subscriptℳplus-or-minus{\cal M}_{\pm}caligraphic_M start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT are effectively closed. Crucially, Mzsubscript𝑀𝑧M_{z}italic_M start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT symmetry guarantees the existence of such a gauge in our case [22], ensuring the Îœđœˆ\nuitalic_Îœ counting is always possible. Apparently, the nodal structure of ΔeffsubscriptΔeff\Delta_{\text{eff}}roman_Δ start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT must be gauge-invariant. Therefore, applying the Chern-vorticity theorem under a proper gauge thoroughly informs the obstructed pairing zeros that are gauge-independent.

Refer to caption
Figure 2: Superconducting gap on the FS. (a)-(c) shows the projected pairing gap for NLSM (ÎŁ=−0.1ÎŁ0.1\Sigma=-0.1roman_ÎŁ = - 0.1) with ÎŒ=−0.2,−0.1,−0.05𝜇0.20.10.05\mu=-0.2,-0.1,-0.05italic_ÎŒ = - 0.2 , - 0.1 , - 0.05, respectively. (d)-(f) shows the projected pairing gap for WSM (ÎŁ=0.1ÎŁ0.1\Sigma=0.1roman_ÎŁ = 0.1) with ÎŒ=−0.2,−0.1,−0.05𝜇0.20.10.05\mu=-0.2,-0.1,-0.05italic_ÎŒ = - 0.2 , - 0.1 , - 0.05, respectively. Note that a Lifshitz transition occurs at |ÎŒ|=|ÎŁ|𝜇Σ|\mu|=|\Sigma|| italic_ÎŒ | = | roman_ÎŁ | where the number of FSs changes.

Finally, let us take the Eu2superscriptsubscript𝐾𝑱2{}^{2}E_{u}start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_E start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT pairing dx⁹(đ€)=k+subscriptđ‘‘đ‘„đ€subscript𝑘d_{x}({\bf k})=k_{+}italic_d start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ( bold_k ) = italic_k start_POSTSUBSCRIPT + end_POSTSUBSCRIPT as an example to numerically explore the nodal nature of the obstructed SC phases. In Fig. 2, we plot the pairing gap function on the corresponding FS and track its evolution as a function of ÎŁÎŁ\Sigmaroman_ÎŁ and ÎŒđœ‡\muitalic_ÎŒ, where the BdG Weyl nodes and nodal loops are highlighted in black. When |ÎŒ|>|ÎŁ|𝜇Σ|\mu|>|\Sigma|| italic_ÎŒ | > | roman_ÎŁ |, we find a single sphere-like FS for both NLSM (ÎŁ<0ÎŁ0\Sigma<0roman_ÎŁ < 0) and WSM (ÎŁ>0ÎŁ0\Sigma>0roman_ÎŁ > 0), which exhibits both a BdG nodal loop at the equator and a pair of BdG Weyl nodes at the poles. This directly follows our analysis of vorticity counting. Further reducing |ÎŒ|𝜇|\mu|| italic_ÎŒ | triggers a Lifshitz transition for NLSM, which makes the point nodes merge. Surprisingly, the merging of these oppositely charged Weyl nodes leads to a second nodal loop at kz=0subscript𝑘𝑧0k_{z}=0italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0. This is in contrast with the expected pair annihilation process in conventional Weyl SCs, as shown in the SM [22]. Meanwhile, a similar Lifshitz transition for WSM shrinks the BdG nodal loop to a point, which further splits into a pair of Weyl nodes. This robust conversion between Weyl nodes and nodal loops for the BdG bands directly arises from the effective closedness of ℳ±subscriptℳplus-or-minus{\cal M}_{\pm}caligraphic_M start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT [22], which is reminiscent of similar behaviors in its normal state h0subscriptℎ0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [26].

Application to j=3/2𝑗32j=3/2italic_j = 3 / 2 SCs. — As a real-world example, we now demonstrate that the BdG line nodes of half-Heusler-based j=3/2𝑗32j=3/2italic_j = 3 / 2 SCs such as YPtBi and LuPdBi are deeply rooted in the mechanism of dipole obstruction.

We start with a generalized Luttinger-Kohn Hamiltonian as a normal state of interest [30, 31],

hL(α)⁹(đ€)superscriptsubscriptâ„Žđżđ›Œđ€\displaystyle h_{L}^{(\alpha)}({\bf k})italic_h start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_α ) end_POSTSUPERSCRIPT ( bold_k ) =\displaystyle== −3⁹α⁹(kx2−ky2)⁹γ1−2⁹3⁹α⁹kx⁹ky⁹γ23đ›Œsuperscriptsubscriptđ‘˜đ‘„2superscriptsubscript𝑘𝑩2subscriptđ›Ÿ123đ›Œsubscriptđ‘˜đ‘„subscript𝑘𝑩subscriptđ›Ÿ2\displaystyle-\sqrt{3}\alpha(k_{x}^{2}-k_{y}^{2})\gamma_{1}-2\sqrt{3}\alpha k_% {x}k_{y}\gamma_{2}- square-root start_ARG 3 end_ARG italic_α ( italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_Îł start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - 2 square-root start_ARG 3 end_ARG italic_α italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT italic_Îł start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (6)
−\displaystyle-- 2⁹3⁹kz⁹kx⁹γ3−2⁹3⁹kz⁹ky⁹γ4+M⁹(đ€)⁹γ5,23subscript𝑘𝑧subscriptđ‘˜đ‘„subscriptđ›Ÿ323subscript𝑘𝑧subscript𝑘𝑩subscriptđ›Ÿ4đ‘€đ€subscriptđ›Ÿ5\displaystyle 2\sqrt{3}k_{z}k_{x}\gamma_{3}-2\sqrt{3}k_{z}k_{y}\gamma_{4}+M({% \bf k})\gamma_{5},2 square-root start_ARG 3 end_ARG italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_Îł start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - 2 square-root start_ARG 3 end_ARG italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT italic_Îł start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT + italic_M ( bold_k ) italic_Îł start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT ,

where Îł1=τxâąÏƒ0subscriptđ›Ÿ1subscriptđœđ‘„subscript𝜎0\gamma_{1}=\tau_{x}\sigma_{0}italic_Îł start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_τ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, Îł2=τyâąÏƒ0subscriptđ›Ÿ2subscript𝜏𝑩subscript𝜎0\gamma_{2}=\tau_{y}\sigma_{0}italic_Îł start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = italic_τ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, Îł3=τzâąÏƒxsubscriptđ›Ÿ3subscript𝜏𝑧subscriptđœŽđ‘„\gamma_{3}=\tau_{z}\sigma_{x}italic_Îł start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT = italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT, Îł4=τzâąÏƒysubscriptđ›Ÿ4subscript𝜏𝑧subscript𝜎𝑩\gamma_{4}=\tau_{z}\sigma_{y}italic_Îł start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT = italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT, Îł5=τzâąÏƒzsubscriptđ›Ÿ5subscript𝜏𝑧subscript𝜎𝑧\gamma_{5}=\tau_{z}\sigma_{z}italic_Îł start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT = italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT, and M⁹(đ€)=kx2+ky2−2⁹kz2đ‘€đ€superscriptsubscriptđ‘˜đ‘„2superscriptsubscript𝑘𝑩22superscriptsubscript𝑘𝑧2M({\bf k})=k_{x}^{2}+k_{y}^{2}-2k_{z}^{2}italic_M ( bold_k ) = italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 2 italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. When α=1đ›Œ1\alpha=1italic_α = 1, hL(1)subscriptsuperscriptℎ1𝐿h^{(1)}_{L}italic_h start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT recovers the standard isotropic Luttinger semimetal, which nicely captures the low-energy normal state of general j=3/2𝑗32j=3/2italic_j = 3 / 2 half-Heusler SCs. While several candidate pairing channels have been theoretically proposed [32, 33, 34, 35, 36], a recent experiment reports a doping-dependent nodeless-to-nodal transition in LuPdBi [20]. This observation strongly supports the mixed-parity singlet-septet pairing Δ~⁹(đ€)=Δ~s+Δ~p~Î”đ€subscript~Δ𝑠subscript~Δ𝑝\tilde{\Delta}({\bf k})=\tilde{\Delta}_{s}+\tilde{\Delta}_{p}over~ start_ARG roman_Δ end_ARG ( bold_k ) = over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT + over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT [21], with Δ~s=iⁱΔsâąÏ„xâąÏƒysubscript~Δ𝑠𝑖subscriptΔ𝑠subscriptđœđ‘„subscript𝜎𝑩\tilde{\Delta}_{s}=i\Delta_{s}\tau_{x}\sigma_{y}over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = italic_i roman_Δ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT and

Δ~p=Δpsubscript~Δ𝑝subscriptΔ𝑝\displaystyle\tilde{\Delta}_{p}=\Delta_{p}over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = roman_Δ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT (34k−τzσ++34k+τzσ−+32kzτ0σx\displaystyle\bigg{(}\frac{3}{4}k_{-}\tau_{z}\sigma_{+}+\frac{3}{4}k_{+}\tau_{% z}\sigma_{-}+\frac{\sqrt{3}}{2}k_{z}\tau_{0}\sigma_{x}( divide start_ARG 3 end_ARG start_ARG 4 end_ARG italic_k start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT + divide start_ARG 3 end_ARG start_ARG 4 end_ARG italic_k start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT - end_POSTSUBSCRIPT + divide start_ARG square-root start_ARG 3 end_ARG end_ARG start_ARG 2 end_ARG italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT (7)
+34k+τxσ+−34k−τxσ−),\displaystyle+\frac{\sqrt{3}}{4}k_{+}\tau_{x}\sigma_{+}-\frac{\sqrt{3}}{4}k_{-% }\tau_{x}\sigma_{-}\bigg{)},+ divide start_ARG square-root start_ARG 3 end_ARG end_ARG start_ARG 4 end_ARG italic_k start_POSTSUBSCRIPT + end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT - divide start_ARG square-root start_ARG 3 end_ARG end_ARG start_ARG 4 end_ARG italic_k start_POSTSUBSCRIPT - end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) ,

which belongs to the A1subscript𝐮1A_{1}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT irrep of Tdsubscript𝑇𝑑T_{d}italic_T start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT group. Specifically, with just Δ~psubscript~Δ𝑝\tilde{\Delta}_{p}over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT, the BdG spectrum of hL(1)subscriptsuperscriptℎ1𝐿h^{(1)}_{L}italic_h start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT displays point nodes at kx,y,zsubscriptđ‘˜đ‘„đ‘Šđ‘§k_{x,y,z}italic_k start_POSTSUBSCRIPT italic_x , italic_y , italic_z end_POSTSUBSCRIPT axes [c.f. Fig. 3(b)]. Turning on Δ~ssubscript~Δ𝑠\tilde{\Delta}_{s}over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT further inflates the point nodes into nodal loops, eventually leading to a full gap when Δ~ssubscript~Δ𝑠\tilde{\Delta}_{s}over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT dominates. However, why the p𝑝pitalic_p-wave pairing can lead to point nodes can be quite puzzling. For example, let us set Δs=0subscriptΔ𝑠0\Delta_{s}=0roman_Δ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0 and focus on the BdG physics along the kzsubscript𝑘𝑧k_{z}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT axis. While Δ~p∌kzâąÏ„0âąÏƒx≠0similar-tosubscript~Δ𝑝subscript𝑘𝑧subscript𝜏0subscriptđœŽđ‘„0\tilde{\Delta}_{p}\sim k_{z}\tau_{0}\sigma_{x}\neq 0over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ∌ italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ≠ 0 for kz≠0subscript𝑘𝑧0k_{z}\neq 0italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ≠ 0, the existence of point nodes implies that the projection of Δ~p⁹(0,0,kz)subscript~Δ𝑝00subscript𝑘𝑧\tilde{\Delta}_{p}(0,0,k_{z})over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ( 0 , 0 , italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ) onto the FS must vanish.

To trace the origin of these point nodes, we note that when α=0đ›Œ0\alpha=0italic_α = 0, hL(0)⁹(đ€)∌τz⊗h0⁹(đ€)similar-tosuperscriptsubscriptℎ𝐿0đ€tensor-productsubscript𝜏𝑧subscriptℎ0đ€h_{L}^{(0)}({\bf k})\sim\tau_{z}\otimes h_{0}({\bf k})italic_h start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( bold_k ) ∌ italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ⊗ italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( bold_k ) exactly describes a Dirac-dipole semimetal [37] that comprises two decoupled copies of the Berry-dipole Hamiltonian at ÎŁ=0ÎŁ0\Sigma=0roman_ÎŁ = 0, up to some parameter rescaling. The block-diagonal nature of hL(0)⁹(đ€)superscriptsubscriptℎ𝐿0đ€h_{L}^{(0)}({\bf k})italic_h start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( bold_k ) suggests an emergent spin conservation symmetry U⁹(1)s𝑈subscript1𝑠U(1)_{s}italic_U ( 1 ) start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT generated by Sz=τzâąÏƒ0subscript𝑆𝑧subscript𝜏𝑧subscript𝜎0S_{z}=\tau_{z}\sigma_{0}italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, with which the two Berry-dipole blocks are carrying opposite spins. Note that an α≠0đ›Œ0\alpha\neq 0italic_α ≠ 0 generally reduces U⁹(1)s𝑈subscript1𝑠U(1)_{s}italic_U ( 1 ) start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT to the time-reversal symmetry 𝒯=iâąÏ„xâąÏƒyⁱ𝒩𝒯𝑖subscriptđœđ‘„subscript𝜎𝑩𝒩\mathcal{T}=i\tau_{x}\sigma_{y}{\cal K}caligraphic_T = italic_i italic_τ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT caligraphic_K, with 𝒩𝒩{\cal K}caligraphic_K the complex conjugation. Besides, hL(0)⁹(đ€)superscriptsubscriptℎ𝐿0đ€h_{L}^{(0)}({\bf k})italic_h start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( bold_k ) fully inherits the C4⁹hsubscriptđ¶4ℎC_{4h}italic_C start_POSTSUBSCRIPT 4 italic_h end_POSTSUBSCRIPT symmetry of h0subscriptℎ0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT with an updated rotation generator Jz=diag⁹(32,12,−12,−32)subscriptđœđ‘§diag32121232J_{z}=\text{diag}(\frac{3}{2},\frac{1}{2},-\frac{1}{2},-\frac{3}{2})italic_J start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = diag ( divide start_ARG 3 end_ARG start_ARG 2 end_ARG , divide start_ARG 1 end_ARG start_ARG 2 end_ARG , - divide start_ARG 1 end_ARG start_ARG 2 end_ARG , - divide start_ARG 3 end_ARG start_ARG 2 end_ARG ) and Mz=iâąÏ„0âąÏƒzsubscript𝑀𝑧𝑖subscript𝜏0subscript𝜎𝑧M_{z}=i\tau_{0}\sigma_{z}italic_M start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = italic_i italic_τ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT. For a finite ÎŒđœ‡\muitalic_ÎŒ, the Dirac-dipole semimetal features a pair of decoupled FSs, one for each spin sector. Hence, all intra-FS pairings are intra-spin pairings, and vice versa.

Notably, the first three terms of Δ~p⁹(đ€)subscript~Î”đ‘đ€\tilde{\Delta}_{p}(\mathbf{k})over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ( bold_k ) describe the intra-spin pairing process. Within each spin sector, these pairings correspond to the Busubscriptđ”đ‘ąB_{u}italic_B start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT pairings in Table 1, all of which are hence dipole-obstructed for hL(0)⁹(đ€)superscriptsubscriptℎ𝐿0đ€h_{L}^{(0)}({\bf k})italic_h start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( bold_k ). Applying the Chern-vorticity theorem in Eq. 2, we find this set of intra-spin pairing terms to carry a nontrivial pairing vorticity of Îœ=−1𝜈1\nu=-1italic_Îœ = - 1 [22], which strictly enforces pairing zeros (or equivalently point nodes) on both the north and south poles of the FS. Meanwhile, the remaining pairing terms of Δ~p⁹(đ€)subscript~Î”đ‘đ€\tilde{\Delta}_{p}(\mathbf{k})over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ( bold_k ) are inter-spin and exhibit no obvious vortex pattern. Nonetheless, both inter-spin pairings vanish on the kzsubscript𝑘𝑧k_{z}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT axis and their existence is thus invisible to the dipole-obstructed pairing zeros. In Fig. 3(a), we numerically confirm the predicted pairing zeros of Δ~p⁹(đ€)subscript~Î”đ‘đ€\tilde{\Delta}_{p}(\mathbf{k})over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ( bold_k ) by mapping out the pairing gap on the FS of the Dirac-dipole semimetal.

Refer to caption
Figure 3: Superconducting gap over the FS of (a) a Dirac-dipole semimetal; (b) a Luttinger semimetal with (Δs,Δp)=(0,0.05)subscriptΔ𝑠subscriptΔ𝑝00.05(\Delta_{s},\Delta_{p})=(0,0.05)( roman_Δ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT , roman_Δ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) = ( 0 , 0.05 ); and (c) a Luttinger semimetal with (Δs,Δp)=(0.01,0.05)subscriptΔ𝑠subscriptΔ𝑝0.010.05(\Delta_{s},\Delta_{p})=(0.01,0.05)( roman_Δ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT , roman_Δ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) = ( 0.01 , 0.05 ), at ÎŒ=−0.2𝜇0.2\mu=-0.2italic_ÎŒ = - 0.2. (d) Spectrum along a path on FS indicated by the green arrow in (c). The purple dots highlight the band degeneracies induced by the dipole obstruction.

We now revisit the isotropic limit (α=1đ›Œ1\alpha=1italic_α = 1) that the j=3/2𝑗32j=3/2italic_j = 3 / 2 SCs follow. Since all Î±đ›Œ\alphaitalic_α-relevant terms in hL(α)superscriptsubscriptâ„Žđżđ›Œh_{L}^{(\alpha)}italic_h start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_α ) end_POSTSUPERSCRIPT vanish on the kzsubscript𝑘𝑧k_{z}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT axis, turning on Î±đ›Œ\alphaitalic_α will not disrupt the obstructed pairing zeros in the Dirac-dipole limit. Moreover, hL(1)⁹(đ€)superscriptsubscriptℎ𝐿1đ€h_{L}^{(1)}({\bf k})italic_h start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ( bold_k ) respects a larger symmetry group (e.g., Tdsubscript𝑇𝑑T_{d}italic_T start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT) than that of hL(0)⁹(đ€)superscriptsubscriptℎ𝐿0đ€h_{L}^{(0)}({\bf k})italic_h start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( bold_k ), which features a three-fold rotation C3,111subscriptđ¶3111C_{3,111}italic_C start_POSTSUBSCRIPT 3 , 111 end_POSTSUBSCRIPT that permutes kx,kysubscriptđ‘˜đ‘„subscript𝑘𝑩k_{x},k_{y}italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT, and kzsubscript𝑘𝑧k_{z}italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT. Hence, BdG point nodes should also emerge at both kxsubscriptđ‘˜đ‘„k_{x}italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT and kysubscript𝑘𝑩k_{y}italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT axes, which is explicitly confirmed in Fig. 3(b). In other words, C3,111subscriptđ¶3111C_{3,111}italic_C start_POSTSUBSCRIPT 3 , 111 end_POSTSUBSCRIPT enables us to designate any of the kx,y,zsubscriptđ‘˜đ‘„đ‘Šđ‘§k_{x,y,z}italic_k start_POSTSUBSCRIPT italic_x , italic_y , italic_z end_POSTSUBSCRIPT-axes as the C4subscriptđ¶4C_{4}italic_C start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT-rotation axis for the previous Dirac-dipole discussion. In this case, all of the six p𝑝pitalic_p-wave-induced point nodes for j=3/2𝑗32j=3/2italic_j = 3 / 2 SCs have a dipole-obstructed origin.

Finally, let us turn on a subdominant Δ~ssubscript~Δ𝑠\tilde{\Delta}_{s}over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT to achieve the BdG line nodes, as shown in Fig. 3(c). Specifically, Δ~ssubscript~Δ𝑠\tilde{\Delta}_{s}over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT lifts the spin-degeneracy of FSs by breaking the inversion symmetry, which leads to accidental inter-FS band crossings that manifest as line nodes. While the dipole-obstructed pairing nodes [purple dots in Fig. 3(d)] now show up at a finite energy E=±Δ~s𝐾plus-or-minussubscript~Δ𝑠E=\pm\tilde{\Delta}_{s}italic_E = ± over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT, their existence ensures the BdG line nodes to necessarily show up even for an arbitrarily small Δ~ssubscript~Δ𝑠\tilde{\Delta}_{s}over~ start_ARG roman_Δ end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT. Note that we have omitted a small asymmetric spin-orbit coupling (ASOC) term in Eq. 6, which is intrinsic to zinc-blende materials. While the ASOC is key to mixing s𝑠sitalic_s and p𝑝pitalic_p-wave pairings, its effect on the BdG band spectrum is similar to that of the s𝑠sitalic_s-wave pairing. A detailed discussion on the ASOC term can be found in the SM [22].

Discussions.— To summarize, we have established a general theoretical framework for comprehending topologically obstructed nodal pairing over a single FS. This leads us to uncover a class of Berry-dipole FSs where the SC state induced by any intra-FS pairing will be dipole-obstructed and nodal. The hidden Berry-dipole physics in the Luttinger-Kohn model further motivates us to explore the j=3/2𝑗32j=3/2italic_j = 3 / 2 pairings in half-Heusler SCs. Focusing on the singlet-septet pairing channel, we find that the dipole obstruction contributes significantly to the BdG line nodes observed in experiments. Notably, Berry and Dirac dipoles, as well as Luttinger semimetals, have been recently established as critical points for delicate topological phases [16, 37, 31]. Therefore, our dipole-obstructed pairing offers the first example of how delicate topological bands can decisively impact the correlated electronic orders in real-world quantum materials.

We emphasize that our Chern-vorticity theorem goes beyond SCs and applies to general electronic orders such as charge-density waves, excitons, magnetism, etc. It would be interesting to explore possible dipole-obstructed phenomena in other non-superconducting Luttinger semimetals such as HgTe [38] and Pr2Ir2O7 [39]. Besides, we find many inter-spin pairings for the j=3/2𝑗32j=3/2italic_j = 3 / 2 system to feature emergent nodal structures in the BdG spectrum, despite carrying no vortex structure. As elaborated in the SM [22] for the Dirac-dipole model, these “unobstructed” pairing zeros arise from the geometric textures of the Bloch states on the FSs, which are beyond the scope of the Chern-vorticity theorem. How to interpret these geometry-relevant nodal pairings is an absolutely intriguing question for future research.

Acknowledgements.
We thank J. Yu, Y. Wang, L.-H. Hu, A. Alexandradinata, and Y.-M. Lu for helpful discussions. P.Z was primarily supported by the Center for Emergent Materials, an NSF MRSEC, under award number DMR-2011876. R.-X.Z. is supported by a start-up fund of the University of Tennessee.

References

  • Bardeen et al. [1957] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Microscopic theory of superconductivity, Phys. Rev. 106, 162 (1957).
  • Sigrist and Ueda [1991] M. Sigrist and K. Ueda, Phenomenological theory of unconventional superconductivity, Rev. Mod. Phys. 63, 239 (1991).
  • Tsuei and Kirtley [2000] C. C. Tsuei and J. R. Kirtley, Pairing symmetry in cuprate superconductors, Rev. Mod. Phys. 72, 969 (2000).
  • Kim et al. [2018] H. Kim, K. Wang, Y. Nakajima, R. Hu, S. Ziemak, P. Syers, L. Wang, H. Hodovanets, J. D. Denlinger, P. M. Brydon, et al., Beyond triplet: Unconventional superconductivity in a spin-3/2 topological semimetal, Science advances 4, eaao4513 (2018).
  • Damascelli et al. [2003] A. Damascelli, Z. Hussain, and Z.-X. Shen, Angle-resolved photoemission studies of the cuprate superconductors, Rev. Mod. Phys. 75, 473 (2003).
  • Fischer et al. [2007] O. Fischer, M. Kugler, I. Maggio-Aprile, C. Berthod, and C. Renner, Scanning tunneling spectroscopy of high-temperature superconductors, Rev. Mod. Phys. 79, 353 (2007).
  • Prozorov and Giannetta [2006] R. Prozorov and R. W. Giannetta, Magnetic penetration depth in unconventional superconductors, Superconductor Science and Technology 19, R41 (2006).
  • Murakami and Nagaosa [2003] S. Murakami and N. Nagaosa, Berry phase in magnetic superconductors, Phys. Rev. Lett. 90, 057002 (2003).
  • Li and Haldane [2018] Y. Li and F. D. M. Haldane, Topological nodal cooper pairing in doped weyl metals, Phys. Rev. Lett. 120, 067003 (2018).
  • Wang and Ye [2016] Y. Wang and P. Ye, Topological density-wave states in a particle-hole symmetric weyl metal, Phys. Rev. B 94, 075115 (2016).
  • Muñoz et al. [2020] E. Muñoz, R. Soto-Garrido, and V. Juričić, Monopole versus spherical harmonic superconductors: Topological repulsion, coexistence, and stability, Phys. Rev. B 102, 195121 (2020).
  • Yu et al. [2022] J. Yu, Y.-A. Chen, and S. Das Sarma, Euler-obstructed cooper pairing: Nodal superconductivity and hinge majorana zero modes, Phys. Rev. B 105, 104515 (2022).
  • Yu et al. [2023] J. Yu, M. Xie, F. Wu, and S. Das Sarma, Euler-obstructed nematic nodal superconductivity in twisted bilayer graphene, Phys. Rev. B 107, L201106 (2023).
  • Wang et al. [2024] Y.-J. Wang, G.-D. Zhou, S.-Y. Peng, B. Lian, and Z.-D. Song, Molecular pairing in twisted bilayer graphene superconductivity, arXiv preprint arXiv:2402.00869  (2024).
  • Sun and Li [2020] C. Sun and Y. Li, ℀2subscriptâ„€2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT topologically obstructed superconducting order, arXiv preprint arXiv:2009.07263  (2020).
  • Nelson et al. [2022] A. Nelson, T. Neupert, A. Alexandradinata, and T. c. v. BzduĆĄek, Delicate topology protected by rotation symmetry: Crystalline hopf insulators and beyond, Phys. Rev. B 106, 075124 (2022).
  • Butch et al. [2011] N. P. Butch, P. Syers, K. Kirshenbaum, A. P. Hope, and J. Paglione, Superconductivity in the topological semimetal yptbi, Phys. Rev. B 84, 220504 (2011).
  • Nakajima et al. [2015] Y. Nakajima, R. Hu, K. Kirshenbaum, A. Hughes, P. Syers, X. Wang, K. Wang, R. Wang, S. R. Saha, D. Pratt, et al., Topological r pdbi half-heusler semimetals: A new family of noncentrosymmetric magnetic superconductors, Science advances 1, e1500242 (2015).
  • Pavlosiuk et al. [2015] O. Pavlosiuk, D. Kaczorowski, and P. Wiƛniewski, Shubnikov-de haas oscillations, weak antilocalization effect and large linear magnetoresistance in the putative topological superconductor lupdbi, Scientific reports 5, 9158 (2015).
  • Ishihara et al. [2021] K. Ishihara, T. Takenaka, Y. Miao, Y. Mizukami, K. Hashimoto, M. Yamashita, M. Konczykowski, R. Masuki, M. Hirayama, T. Nomoto, R. Arita, O. Pavlosiuk, P. Wiƛniewski, D. Kaczorowski, and T. Shibauchi, Tuning the parity mixing of singlet-septet pairing in a half-heusler superconductor, Phys. Rev. X 11, 041048 (2021).
  • Brydon et al. [2016] P. M. R. Brydon, L. Wang, M. Weinert, and D. F. Agterberg, Pairing of j=3/2𝑗32j=3/2italic_j = 3 / 2 fermions in half-heusler superconductors, Phys. Rev. Lett. 116, 177001 (2016).
  • [22] See the Supplemental Material for (i) a detailed proof of the Chern-vorticity theorem; (ii) detailed analysis for general pairings in Berry-dipole semimetals and inter-spin pairings in Dirac-dipole semimetals; and (iii) case studies on WSM, Berry dipole SC, and perturbations from ASOC.
  • Bobrow et al. [2020] E. Bobrow, C. Sun, and Y. Li, Monopole charge density wave states in weyl semimetals, Phys. Rev. Res. 2, 012078 (2020).
  • Bultinck et al. [2020] N. Bultinck, S. Chatterjee, and M. P. Zaletel, Mechanism for anomalous hall ferromagnetism in twisted bilayer graphene, Phys. Rev. Lett. 124, 166601 (2020).
  • Zhu and Alexandradinata [2023] P. Zhu and A. Alexandradinata, Anomalous shift and optical vorticity in the steady photovoltaic current, arXiv preprint arXiv:2308.08596  (2023).
  • Sun et al. [2018] X.-Q. Sun, S.-C. Zhang, and T. c. v. BzduĆĄek, Conversion rules for weyl points and nodal lines in topological media, Phys. Rev. Lett. 121, 106402 (2018).
  • Graf and PiĂ©chon [2023] A. Graf and F. PiĂ©chon, Massless multifold hopf semimetals, Phys. Rev. B 108, 115105 (2023).
  • Zhuang et al. [2024] Z.-Y. Zhuang, C. Zhang, X.-J. Wang, and Z. Yan, Berry-dipole semimetals, arXiv preprint arXiv:2404.10049  (2024).
  • Tyner and Sur [2024] A. C. Tyner and S. Sur, Dipolar weyl semimetals, Phys. Rev. B 109, L081101 (2024).
  • Hu and Zhang [2023] L.-H. Hu and R.-X. Zhang, Topological superconducting vortex from trivial electronic bands, Nature Communications 14, 640 (2023).
  • Zhu and Zhang [2023] P. Zhu and R.-X. Zhang, Delicate topology of luttinger semimetal, arXiv preprint arXiv:2308.05793  (2023).
  • Yang et al. [2017] W. Yang, T. Xiang, and C. Wu, Majorana surface modes of nodal topological pairings in spin-3232\frac{3}{2}divide start_ARG 3 end_ARG start_ARG 2 end_ARG semimetals, Phys. Rev. B 96, 144514 (2017).
  • Roy et al. [2019] B. Roy, S. A. A. Ghorashi, M. S. Foster, and A. H. Nevidomskyy, Topological superconductivity of spin-3/2323/23 / 2 carriers in a three-dimensional doped luttinger semimetal, Phys. Rev. B 99, 054505 (2019).
  • Boettcher and Herbut [2018] I. Boettcher and I. F. Herbut, Unconventional superconductivity in luttinger semimetals: Theory of complex tensor order and the emergence of the uniaxial nematic state, Phys. Rev. Lett. 120, 057002 (2018).
  • Yu and Liu [2018] J. Yu and C.-X. Liu, Singlet-quintet mixing in spin-orbit coupled superconductors with j=32𝑗32j=\frac{3}{2}italic_j = divide start_ARG 3 end_ARG start_ARG 2 end_ARG fermions, Phys. Rev. B 98, 104514 (2018).
  • Venderbos et al. [2018] J. W. F. Venderbos, L. Savary, J. Ruhman, P. A. Lee, and L. Fu, Pairing states of spin-3232\frac{3}{2}divide start_ARG 3 end_ARG start_ARG 2 end_ARG fermions: Symmetry-enforced topological gap functions, Phys. Rev. X 8, 011029 (2018).
  • Zhu et al. [2023] P. Zhu, A. Alexandradinata, and T. L. Hughes, ℀2subscriptâ„€2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT spin hopf insulator: Helical hinge states and returning thouless pump, Phys. Rev. B 107, 115159 (2023).
  • Piotrzkowski et al. [1965] R. Piotrzkowski, S. Porowski, Z. Dziuba, J. Ginter, W. Giriat, and L. Sosnowski, Band structure of hgte, physica status solidi (b) 8, K135 (1965).
  • Kondo et al. [2015] T. Kondo, M. Nakayama, R. Chen, J. Ishikawa, E.-G. Moon, T. Yamamoto, Y. Ota, W. Malaeb, H. Kanai, Y. Nakashima, et al., Quadratic fermi node in a 3d strongly correlated semimetal, Nature communications 6, 10042 (2015).