Visualizing the Odd-Parity Superconducting Order Parameter and Its Quasiparticle Surface Band in UTe2

TL;DR

This paper uses advanced imaging techniques to visualize the unique superconducting properties of UTe2, revealing a non-chiral topological superconducting state.

Contribution

The study provides direct visualization of the quasiparticle surface band and confirms the odd-parity, spin-triplet superconducting order in UTe2.

Findings

An intense zero-energy Andreev conductance peak is observed at the UTe2 (0–11) crystal termination.

Quasiparticle interference patterns reveal a sextet of interference wavevectors, indicating a specific dispersion of Bogoliubov quasiparticles.

The findings are consistent with a spin-triplet, time-reversal conserving, odd-parity superconducting order in UTe2.

Abstract

A distinctive identifier of nodal intrinsic topological superconductivity (ITS) would the appearance of an Andreev bound state on crystal surfaces parallel to the nodal axis, in the form of a topological quasiparticle surface band (QSB) appearing only for \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T < T_{C}$$\end{document}T<TC. Moreover, the theory shows that specific QSB characteristics observable in tunneling to an s-wave superconductor can distinguish between chiral and non-chiral ITS order parameter \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Delta_{{\boldsymbol{k}}}$$\end{document}Δk. To search for such phenomena in UTe2, s-wave superconductive scan-tip scanning tunneling microscopy (STM) imaging was employed. It reveals an intense zero-energy Andreev conductance maximum at the UTe2 (0–11) crystal termination. The development of the zero-energy Andreev conductance peak into two finite-energy particle-hole symmetric conductance maxima as the tunnel barrier is reduced and then signifies that UTe2 superconductivity is non-chiral. Quasiparticle interference imaging (QPI) for an ITS material should be dominated by the QSB for energies within the superconductive energy gap \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\left| E \right| \le {\Delta }$$\end{document}E≤Δ, so that bulk \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Delta_{{\boldsymbol{k}}}$$\end{document}Δk characteristics of the ITS can only be detected excursively. Again using a superconducting scan-tip, the in-gap quasiparticle interference patterns of the QSB of UTe2 were visualized. Specifically, a band of Bogoliubov quasiparticles appears as a characteristic sextet \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\boldsymbol{q}}_{i} :i = 1 - 6{ }$$\end{document}qi:i=1-6 of interference wavevectors, showing that QSB dispersions \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\boldsymbol{k}$$\end{document}k(E) occur only for energies \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\left| E \right| \le \Delta_{\max }$$\end{document}E≤Δmax and only within the range of Fermi momenta projected onto the (0–11) crystal surface. In combination, these phenomena are consistent with a bulk \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Delta_{{\boldsymbol{k}}}$$\end{document}Δk exhibiting spin-triplet, time-reversal conserving, odd-parity, a-axis nodal, B3u symmetry in UTe2.

Figures

Schematic of SIP (*s*-wave to *p*-wave) tunneling through a superconductive topological surface state. This is the basic technique implemented throughout this paper

Full text

Intrinsic Topological Superconductivity

For spin-triplet superconductors [1–4], the order parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}} = \left( {\begin{array}{*{20}c} {{\Delta }_{{{\boldsymbol{k}} \uparrow \uparrow }} } & {{\Delta }_{{{\boldsymbol{k}} \uparrow \downarrow }} } \\ {{\Delta }_{{{\boldsymbol{k}} \downarrow \uparrow }} } & {{\Delta }_{{{\boldsymbol{k}} \downarrow \downarrow }} } \\ \end{array} } \right)$$\end{document} with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{ - {\boldsymbol{k}}}}^{{\mathrm{T}}} = - {\Delta }_{{\boldsymbol{k}}}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{\boldsymbol{k}} = {\Delta }_{\boldsymbol{k}}^{{\mathrm{T}}}$$\end{document} is also represented in the d-vector notation as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}} \equiv {\Delta }_{0} \left( {{\boldsymbol{d}} \cdot {\boldsymbol{\sigma}}} \right){{i}}\sigma_{2}$$\end{document} where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma_{i}$$\end{document} are the Pauli matrices. In principle, such systems are ITS whose signature is the existence of an odd-parity bulk superconducting energy gap, along with the presence of symmetry-protected gapless topological quasiparticle surface bands of Bogoliubov quasiparticles within that gap. When such superconductors are topological [5], it is not because of electronic band-structure topology but because \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} itself exhibits topologically non-trivial properties [6]. The search for technologically viable ITS is now a forefront of quantum matter research [7]. Here, we summarize recent scanning tunneling microscopy studies, specifically using superconducting scan-tips in the Josephson and Andreev modes as explained below, of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} and associated quasiparticle surface bands in the candidate ITS material, UTe_2_.

Prevenient Charge-Density-Wave State in UTe2

Well above the superconducting critical temperate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{C}$$\end{document} , three charge-density-wave (CDW) states with distinct wavevectors are observed [8] at the equivalent (0–11) cleave surface of UTe_2_ where our studies are carried out. These states have not been detected in bulk [9, 10]. Upon entering the superconductive state, three pair-density-wave (PDW) states with distinct wavevectors are observed, through their periodically modulating superconducting energy gap, at the identical wavevectors as the prevenient CDWs [11]. These phenomena are all consistent with induction of these UTe_2_ PDW states due to the interactions between the prevenient CDW states and the superconductivity. At present, these UTe_2_ PDW states have not played a role in determination of the symmetry of the bulk superconductor order parameter.

Superconducting Order Parameter of UTe2

The recently discovered superconductor UTe_2_ is the leading candidate to be a 3D nodal spin-triplet superconductor [12, 13] and thus an ITS. The crystal symmetry point-group is D2h so that there are four possible odd-parity order-parameter symmetries designated Au, B1u, B2u and B3u. All of these preserve time-reversal symmetry: Au is fully gapped, whereas B1u, B2u and B3u have nodes in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta_{{\boldsymbol{k}}}$$\end{document} , whose axial alignment is along lattice vectors \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{c}},{\boldsymbol{b}}\;{\mathrm{or}}\;{\boldsymbol{a}}$$\end{document} , respectively. If they are accidentally degenerate, linear combinations of these order parameters are also possible, which break point-group and time-reversal symmetries, resulting in a chiral QSB along with persistent surface supercurrents orthogonal to the nodal axis [6, 7]. For UTe_2_, there are two chiral states of particular interest with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta_{{\boldsymbol{k}}}$$\end{document} nodes aligned with the crystal c-axis, and two with nodes aligned with the a-axis. Identifying which (if any) of these \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta_{{\boldsymbol{k}}}$$\end{document} exists in UTe_2_ is key to demonstrating and utilizing the novel physics of this material.

However, this has proven a challenging objective [14]. For example, a magnetic susceptibility upon entering the superconducting phase that is equivalent to Pauli paramagnetism is deduced from minimal suppressions of the Knight shift [15] and used to adduce spin-triplet pairing. Some NMR studies measuring the change of the spin susceptibility across \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{{\mathrm{c}}}$$\end{document} report a decrease in the Knight shift in all directions and hypothesize the Au state [15], whereas other NMR studies detect a reduction in the Knight shift along the b and c axes only, thence hypothesizing B3u state [16]. Magnetic field orientation of the thermal conductivity indicates point nodes parallel to the crystal a-axis [17], whereas other field-oriented thermal conductivity measurements [18] report isotropic results and hypothesize an Au symmetry. Field-oriented specific heat measurements reveal peaks around the crystal a-axis implying point nodes oriented along this direction and hypothesize an order parameter with chiral Au + iB3u or helical B3u symmetries [19]. Some electronic specific heat studies report two specific heat peaks and hypothesize a chiral Au + iB1u or B2u + iB3u order parameter [20], whereas other specific heat studies detect only a single specific heat peak and thus hypothesize a single-component order parameter [21]. London penetration depth measurements of superfluid density report anisotropic saturation consistent with nodes along the a-axis suggesting B3u symmetry pairing for a cylindrical Fermi surface [22], while other penetration depth measurements exhibiting an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n \le 2$$\end{document} power law dependence of the penetration depth on temperature motivate a hypothesis of B3u + iAu pairing symmetry [23]. Scanning tunneling microscopy experiments in the (0–11) plane parallel to a-axis show energy-reversed particle-hole symmetry breaking of some electronic-structure elements at opposite UTe_2_ step edges [24] with the consequent hypothesis of a chiral surface state B1u + iB2u whose nodes are aligned to the a-axis. Polar Kerr effect measurements report a field-induced Kerr rotation indicating the presence of time-reversal symmetry breaking and hypothesize chiral B2u + iB3u or Au + iB1u pairing [20] with nodes aligned to the c-axis, whereas other polar Kerr effect measurements report no detectable spontaneous Kerr rotation [25]. However, until recently, no tunneling spectroscopic measurements of Δ_k_, which could differentiate directly between these scenarios, had been reported.

Scanned Andreev Tunneling Microscopy

Ideally, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta_{{\boldsymbol{k}}}$$\end{document} of UTe_2_ might be established by using Bogoliubov quasiparticle interference imaging, a recognized technique for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta_{{\boldsymbol{k}}}$$\end{document} determination in complex superconductors [26–34]. However, odd-parity superconductors should support a topological [28, 35] QSB on crystal termination surfaces only for energies within the superconductive energy gap \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| E \right| \le {\Delta }$$\end{document} . Classically, odd-parity superconductors exhibit zero-energy surface Andreev bound states [36–40], which are generated by the universal \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi$$\end{document} -phase-shift during Andreev reflections from the odd-parity pair potential \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} . More intriguingly, ITS [5, 41] exists most simply in the case of odd-parity spin-triplet superconductors. Thus, a definitive characteristic [28] of an ITS would be a topological quasiparticle surface band with momentum–energy relationship k(E) existing only for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T < T_{c}$$\end{document} and for energies \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| E \right| \le {\Delta }$$\end{document} within the maximum superconducting energy gap [28, 42–52].

QPI is a tunneling spectroscopic technique for establishing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} in unconventional superconductors. However, normal-tip QPI had also proven ineffective for determining \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} of UTe_2_ because even at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T/T_{c} \underset{\raise0.3em\hbox{\smash{\scriptscriptstyle\thicksim}}}{ < } 1/7$$\end{document} a typical quasiparticle density-of-states spectrum \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( {E \le {\Delta }_{0} } \right)$$\end{document} is essentially metallic with only tenuous hints of opening the bulk \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} [8, 11]. The classic QPI signature [26] of scattering interference between k-space locations defined by a bulk superconductive \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} had proven impossible to detect, apparently because the extraordinarily high \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( {E \le {\Delta }_{0} } \right)$$\end{document} of the QSB overwhelms any tunneling conductance signal from the 3D quasiparticles. Hence, the possibility of using a superconductive tip [53–58] to carry out tunneling spectroscopy in the Andreev mode was explored. There are two channels for conduction from the fully gapped s-wave superconductive tip to a nodal spin-triplet superconductor: (a) single-electron tunneling for which the minimum voltage bias is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| V \right| > \Delta_{{{\mathrm{tip}}}} /e$$\end{document} due to the energy cost \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{{\mathrm{tip}}}}$$\end{document} of creating an unpaired electron in the superconducting tip; (b) Andreev reflection of pairs of sub-gap quasiparticles allowing the transfer of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2e$$\end{document} across the junction, thus generating strong conductance at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| V \right| < \Delta_{{{\mathrm{tip}}}} /e$$\end{document} . Hence, in principle, there are strong advantages to using scanned Andreev tunneling spectroscopy for ITS studies, especially that QSB quasiparticles at the interface between sample and tip predominate the Andreev process (as shown schematically in Fig. 1) and that the order-parameter symmetry difference between sample and tip does not preclude the resulting zero-bias Andreev conductance.Fig. 1. Schematic of SIP (s-wave to p-wave) tunneling through a superconductive topological surface state. This is the basic technique implemented throughout this paper

Generally, in superconductive-tip scanned Josephson tunneling microscopy, the electron-pair density in a superconductor*, * \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho_{P} \left( {\boldsymbol{r}} \right)$$\end{document} , is visualized by measuring Josephson critical-current \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$I_{J} \left( {\boldsymbol{r}} \right){ }$$\end{document} from a superconducting STM tip [59], since \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho_{P} \left( {\boldsymbol{r}} \right) \propto I_{J}^{2} \left( {\boldsymbol{r}} \right)R_{N}^{2} \left( {\boldsymbol{r}} \right)$$\end{document} where RN is the normal-state junction resistance [60, 61]. However, thermal fluctuation energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{B} T$$\end{document} typically greatly exceeds the Josephson energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{J}$$\end{document} so that the tip-sample Josephson junction exhibits a phase-diffusive [62–64] steady-state electron-pair current \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$I_{P} \left( V \right) = \frac{1}{2}I_{J}^{2} ZV/\left( {V^{2} + V_{c}^{2} } \right)$$\end{document} at voltage \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V{ }$$\end{document} , where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V_{c} = 2eZk_{B} T/\hbar$$\end{document} and Z is the high-frequency junction impedance. In this case, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$dI_{P} /dV \equiv g\left( V \right) = \frac{1}{2}I_{J}^{2} Z(V_{c}^{2} - V^{2} )/\left( {V^{2} + V_{c}^{2} } \right)^{2}$$\end{document} so that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g\left( 0 \right) \propto I_{J}^{2}$$\end{document} . The key consequence is that spatially resolved measurements of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g\left( {{\boldsymbol{r}}, 0} \right){ }$$\end{document} using superconductive-tip STM at sub-kelvin temperatures now provide a practical technique to visualize electron-pair density \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho_{P} \left( {\boldsymbol{r}} \right) \propto g\left( {{\boldsymbol{r}},0} \right)R_{N}^{2} \left( {\boldsymbol{r}} \right)$$\end{document} at the atomic scale. Technically closely related is scanned Andreev tunneling microscopy (SATM)[65], which, in theory, is highly advantageous for studying ITS. SATM measures the differential Andreev conductance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right) = dI/dV\left( {{\boldsymbol{r}},V} \right)$$\end{document} and, in the case of ITS, exhibits unique phenomena due to the fact that tunneling occurs from an s-wave scan-tip to a p-wave ITS through its QSB (Fig. 1).

Modeling Andreev Tunneling Spectroscopy for ITS

Novel models are required to understand SATM from an even-parity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} (e.g., s-wave) superconducting scan-tip to an odd-parity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} (e.g., p-wave) ITS sample. Most simply, a nodal spin-triplet p-wave superconductor on a spherical Fermi surface within a cubic 3D Brillouin zone (BZ) exhibits two nodal points at ± kn. Its Hamiltonian is:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H = \mathop \sum \limits_{{k_{x} }} \mathop \sum \limits_{{{\boldsymbol{k}}_{ \bot } }} \psi^{ + } \left( {k_{x} ,{\boldsymbol{k}}_{ \bot } } \right) h\left( {k_{x} , {\boldsymbol{k}}_{ \bot } } \right)\psi \left( {k_{x} ,{\boldsymbol{k}}_{ \bot } } \right).$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi^{T} \left( {\boldsymbol{k}} \right) = \left( {c_{{{\boldsymbol{k}} \uparrow }} ,c_{{{\boldsymbol{k}} \downarrow }} ,c_{{ - {\boldsymbol{k}} \uparrow }}^{ + } ,c_{{ - {\boldsymbol{k}} \downarrow }}^{ + } } \right)$$\end{document} is the Nambu fermion operator, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h\left( {k_{x} ,{\boldsymbol{k}}_{ \bot } } \right)$$\end{document} is a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4 \times 4$$\end{document} matrix, containing the information on both band structure and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} [65]. Considering only a particular 2D slice of the 3D Brillouin zone with a fixed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{x}$$\end{document} , its Hamiltonian \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h\left( {k_{x} , {\boldsymbol{k}}_{ \bot } } \right)$$\end{document} is that of a 2D superconductor within a 2D Brillouin zone spanned by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{k}}_{ \bot }$$\end{document} . The 2D states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| {k_{x} } \right| < \left| {k_{n} } \right|$$\end{document} are topological and those \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|k_{x} | > |k_{n} |$$\end{document} are non-topological. The essential signature of such physics is the presence of a QSB also termed an Andreev bound state [6], on the edges of each 2D slice for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|k_{x} \left| { < |k_{n} } \right|$$\end{document} . The 2D Brillouin zone of any crystal surface parallel to the nodal axis of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} has a line of zero-energy QSB states, the so-called Bogoliubov-Fermi Arc, that should in theory connect the two points representing the projections of the nodal wavevectors \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm {\boldsymbol{k}}_{n}$$\end{document} onto this 2D zone. Calculation of the density of QSB quasiparticle states versus energy, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( E \right)$$\end{document} from the QSB dispersion \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{k}}\left( E \right)$$\end{document} , yields a continuum in the range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- {\Delta }_{0} \le E \le {\Delta }_{0}$$\end{document} , with a sharp central peak at E = 0 due to the Bogoliubov–Fermi arc. In this picture, the presence or absence of a gapless QSB on a given surface of a 3D crystal, a zero-energy peak in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( E \right)$$\end{document} from the QSB Fermi-arcs, and the response of the QSB to breaking specific symmetries can reveal the symmetry of the 3D \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} .

Consider an s-wave superconducting tip (Nb) and a nodal p-wave superconductor (UTe_2_), which sustain a QSB within the interface and are connected by tunneling (SIP model). The Hamiltonian of the SIP model has three elements: H = HNb +  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{{{\mathrm{UTe}}_{2} }}$$\end{document}  + HT. Here HNb is the Hamiltonian for an ordinary s-wave superconductor given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{{{\mathrm{Nb}}}} \left( {\boldsymbol{k}} \right) = \left( {\begin{array}{*{20}c} {\epsilon_{{{\mathrm{Nb}}}} \left( {\boldsymbol{k}} \right)\sigma_{0} } & {{\Delta }_{{{\mathrm{Nb}}}} \left( {i\sigma_{2} } \right)} \\ {{\Delta }_{{{\mathrm{Nb}}}}^{*} \left( { - i\sigma_{2} } \right)} & { - \epsilon_{{{\mathrm{Nb}}}} \left( { - {\boldsymbol{k}}} \right)\sigma_{0} } \\ \end{array} } \right)$$\end{document} . Here \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon_{Nb} \left( {\boldsymbol{k}} \right)$$\end{document} is the band-structure model for Nb and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{{\mathrm{Nb}}}}$$\end{document} is the Nb superconducting order parameter, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma_{0,1,2,3}$$\end{document} are the four components of Pauli matrices. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{{{\mathrm{UTe}}_{2} }}$$\end{document} is the Hamiltonian of the putative p-wave superconductor with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {\begin{array}{*{20}l} {\epsilon_{{{\mathrm{UTe}}_{2} }} \left( {\boldsymbol{k}} \right)\sigma_{0} } \hfill & {{\Delta }_{{{\mathrm{UTe}}_{2} }} \left( {\boldsymbol{k}} \right)} \hfill \\ {{\Delta }_{{{\mathrm{UTe}}_{2} }}^{ + } \left( {\boldsymbol{k}} \right)} \hfill & { - \epsilon_{{{\mathrm{UTe}}_{2} }} \left( { - {\boldsymbol{k}}} \right)\sigma_{0} } \hfill \\ \end{array} } \right)$$\end{document} . Here \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon_{{{\mathrm{UTe}}_{2} }} \left( {\boldsymbol{k}} \right)$$\end{document} is the band-structure model containing the relevant Fermi surface, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{{\mathrm{UTe}}_{2}}}\left({\boldsymbol{k}} \right)$$\end{document} is a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2 \times 2$$\end{document} spin-triplet pairing matrix given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{{\mathrm{UTe}}_{2} }} \left( {\boldsymbol{k}} \right) \equiv {\Delta }_{{{\mathrm{UTe}}_{2} }} i\left( {{\boldsymbol{d}} \cdot {\boldsymbol{\sigma}}} \right)\sigma_{2}$$\end{document} . HT is the tunneling Hamiltonian between the two superconductors \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{{\mathrm{T}}} = - \left| M \right|\mathop \sum \limits_{{{\boldsymbol{k}}\parallel }} [\psi_{{{\mathrm{Nb}},{\boldsymbol{k}}\parallel }}^{*} \sigma_{3} \otimes \sigma_{0} \psi_{{{\mathrm{UTe}}_{2},{\boldsymbol{k}}\parallel }} \left( {\boldsymbol{k}} \right) + h.c.]$$\end{document} ; \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{k}}_{\parallel }$$\end{document} is the momentum in the plane parallel to the interface, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi$$\end{document} is the four-component fermion field localizing on the adjacent planes of Nb and UTe_2_, and | \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M$$\end{document} | is the tunneling matrix element. To simplify calculation, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon_{{{\mathrm{Nb}}}} \left( {\boldsymbol{k}} \right)$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon_{{{\mathrm{UTe}}_{2} }} \left( {\boldsymbol{k}} \right)$$\end{document} are approximated as single bands via a nearest-neighbor tight-binding dispersion.

For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{{{\mathrm{UTe}}_{2} }}$$\end{document} , two scenarios were then considered: (1) chiral pairing state Au + iB3u with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{d}}\left( {\boldsymbol{k}} \right) = \left( {0,k_{y} + ik_{z} ,ik_{y} + k_{z} } \right)$$\end{document} and (2) non-chiral pairing state B3u with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{d}}\left( {\boldsymbol{k}} \right) = \left( {0,k_{z} ,k_{y} } \right)$$\end{document} . In both examples, the two nodes of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} lie along the a-axis as in Fig. 1, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{{\mathrm{UTe}}_{2} { }}} = \frac{1}{5}{\Delta }_{{{\mathrm{Nb}}}}$$\end{document} approximate the ratio of maximum energy gaps of Nb and UTe_2_. First, for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right| = 0$$\end{document} the spectrum of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{{{\mathrm{UTe}}_{2} }}$$\end{document} was solved exactly. The quasiparticle eigenstates \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E\left( {k_{x} = 0,k_{y} } \right)$$\end{document} versus \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{y}$$\end{document} have been predicted for the chiral, time reversal symmetry breaking, p-wave order parameter with Au + iB3u symmetry. Here, a chiral QSB spans the full energy range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- {\Delta }_{{{\mathrm{UTe}}_{2} }} \le E \le {\Delta }_{{{\mathrm{UTe}}_{2} }}$$\end{document} , crossing the Fermi level (E = 0) and generating a finite density of quasiparticle states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( {\left| E \right| < {\Delta }_{{{\mathrm{UTe}}_{2} }} } \right).$$\end{document} The quasiparticle spectrum versus \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{y}$$\end{document} at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{x} = 0$$\end{document} were predicted for non-chiral, time-reversal symmetry conserving p-wave order parameter with B3u symmetry. Here, two non-chiral QSBs also span the full energy range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- {\Delta }_{{{\mathrm{UTe}}_{2} }} \le E \le {\Delta }_{{{\mathrm{UTe}}_{2} }}$$\end{document} , and feature E = 0 states, thus generating a finite \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( {\left| E \right| < {\Delta }_{{{\mathrm{UTe}}_{2} }} } \right)$$\end{document} . Although these QSBs have dispersion in both the positive and negative \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{y}$$\end{document} directions and can backscatter, their gaplessness is protected by time-reversal symmetry with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T^{2} = - I$$\end{document} .

To distinguish a chiral from non-chiral \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} by using SATM within the SIP model requires quantitative calculation of the Andreev conductance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( V \right) = dI/dV|_{{{\mathrm{SIP}}}}$$\end{document} between Nb and UTe_2_ using the QSB to demonstrate that a sharp \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( V \right)$$\end{document} peak should occur surrounding zero-bias [65]. Because Andreev reflection of QSB quasiparticles allows highly efficient transfer of charge \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2e$$\end{document} across the junction, its sharpness is robust, meaning that Andreev transport between s-wave/p-wave electrodes through a QSB makes scanned Andreev tunneling spectroscopy an ideal new approach for studying superconductive topological quasiparticle surface bands of ITS. In the limit where the tunneling matrix element to the s-wave electrode \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right| \to 0$$\end{document} , these phenomena are indistinguishable but, as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right|$$\end{document} increases, the wavefunctions of the Nb overlap those of UTe_2_ allowing detection of the QSB quasiparticles at the s-wave electrode. The quasiparticle bands within the SIP interface between Nb and UTe_2_ for the chiral order parameter A_u_ + iB3u symmetry as a function of increasing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right|$$\end{document} were predicted. With increasing | \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M|\sim 1/R$$\end{document} where R is the SIP tunnel junction resistance, the proximity effect of the s-wave electrode generates two chiral QSBs for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| E \right| < {\Delta }_{{{\mathrm{UTe}}_{2} { }}}$$\end{document} , both of which cross E = 0. Hence, for the chiral \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} , the zero-energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( E \right)$$\end{document} will be virtually unperturbed by increasing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right|$$\end{document} . Likewise, the QSB within the SIP interface as a function of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right|$$\end{document} for the non-chiral order parameter with B3u symmetry was also predicted. When \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right| \to 0$$\end{document} , the non-chiral QSB crosses E = 0. But, with increasing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right|\sim 1/R$$\end{document} , time-reversal symmetry breaking due to interactions with the s-wave electrode split the QSB of quasiparticle into two, neither of which cross \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E = 0$$\end{document} . This reveals that the zero-energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( E \right)$$\end{document} peak must split as the zero-energy quasiparticles of the QSB disappear, generating two particle-hole symmetric \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( E \right)$$\end{document} maxima at finite energy. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( 0 \right)$$\end{document} is quantitatively predicted to split into two particle-hole symmetric \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( E \right)$$\end{document} maxima as a function of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right|$$\end{document} for a chiral \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} but not for a non-chiral \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} . Thus, in theory, Andreev tunneling between an s-wave electrode and a p-wave superconductor through the latter’s QSB allows a non-chiral pairing state to be distinguished from a chiral pairing state [65].

Modeling the QPI signature of the QSB was the next challenge. Here it is the normal state electronic structure of UTe_2_ forms the basis upon which \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta_{{\boldsymbol{k}}}$$\end{document} phenomenology emerges at lower temperatures. Atomic-resolution differential tunneling conductance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g\left( {{\boldsymbol{r}},V} \right) \equiv dI/dV\left( {{\boldsymbol{r}},V} \right)$$\end{document} imaging visualizes the density-of-states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( {{\boldsymbol{r}},E} \right)$$\end{document} , and its Fourier transform \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g\left( {{\boldsymbol{q}},E} \right) \propto N\left( {{\boldsymbol{q}},E} \right)$$\end{document} can be used to establish electronic-structure characteristics. Hence, a conventional model of the bulk first BZ of UTe_2_ sustains a two-band Fermi surface (FS) as now widely hypothesized [66, 67]. Quantitative predictions for the normal state QPI in UTe_2_ then require a Hamiltonian \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{{{\mathrm{UTe}}_{2} }} = \left( {\begin{array}{*{20}c} {H_{U - U} } & {H_{U - Te} } \\ {H_{U - Te}^{ + } } & {H_{Te - Te} } \\ \end{array} } \right)$$\end{document} such that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{U - U}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{Te - Te}$$\end{document} describe, respectively, the two uranium and tellurium orbitals and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{U - Te}$$\end{document} their hybridization. From this, one anticipates strong scattering interference with a sextet of wavevectors \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{p}}_{i} :i = 1 - 6$$\end{document} viewed from the (001) plane, where a is the x-axis unit-cell distance and b is the y-axis unit-cell distance.Wavevector \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{p}}_{1}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{p}}_{2}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{p}}_{3}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{p}}_{4}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{p}}_{5}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{p}}_{6}$$\end{document} Coordinate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {\frac{2\pi }{a},\frac{2\pi }{b}} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.29,0} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.43, 1} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.29,2} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0, 2} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( { - 0.14, 1} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.57,0} \right)$$\end{document}

However, the natural cleave surface of UTe_2_ crystal is not (001) but rather (0–11), here shown schematically in Fig. 2a, and it is this surface that the scan-tip approaches perpendicularly. To clarify the normal-state band structure and quasiparticle interference viewed from (0–11) plane, the k-space joint density of states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J\left( {{\boldsymbol{q}},E} \right)$$\end{document} was calculated at the (011) plane using the UTe_2_ FS that takes into account the uranium f orbital spectral weight. The sextet of scattering wavevectors \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{p}}_{i} :i = 1 - 6$$\end{document} derived heuristically above is then revealed as primary peaks in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J\left( {{\boldsymbol{q}},E} \right)$$\end{document} . Here, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J\left( {{\boldsymbol{q}},E} \right)$$\end{document} for the same band-structure model has been calculated but viewed along the normal to the (0–11) plane [68], where the y-coordinates of the (0–11) sextet become \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{1,y} = {\boldsymbol{p}}_{1,y} sin\theta$$\end{document} where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta =$$\end{document} 24° and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c^{*}$$\end{document} is the (0–11) surface y:z-axis lattice periodicity, as indicated by the colored arrows in Ref. [68].Fig. 2a Schematic of (0–11) cleave surface of UTe_2_ shown in relative orientation to the STM tip tunneling direction. Measured high-resolution T(r) at low junction resistance (Is = 3 nA, Vs = 5 mV), clarifying two types of Te atom in light and dark blue. The U atoms are indicated in red. b Typical topographic image T(r) of UTe_2_ (0–11) surface measured with a superconducting tip at T = 280 mK (Is = 0.5 nA, Vs = 30 mV). Inset: measured T(q), the Fourier transform of T(r) in b, with the surface reciprocal-lattice points labeled as dashed red circlesWavevector \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{1}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{2}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{3}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{4}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{5}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{6}$$\end{document} Coordinate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {\frac{2\pi }{a},\frac{2\pi }{{c^{*} }}} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.29,0} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.43, 0.5} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.29,1} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0, 1} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( { - 0.14, 0.5} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {0.57,0} \right)$$\end{document}

In UTe_2_, the Au state should be completely gapped on both Fermi surfaces whereas B2u and B3u states could exhibit point nodes along the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{z} -$$\end{document} axis, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{y} -$$\end{document} axis and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{x} -$$\end{document} axis, respectively. These bulk Bogoliubov eigenstates are described by the dispersion

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{{\boldsymbol{k}}} = \sqrt {\xi_{{\boldsymbol{k}}}^{2} + {\Delta }^{2} \left( {\left| {{\boldsymbol{d}}\left( {\boldsymbol{k}} \right)} \right|^{2} \pm \left| {{\boldsymbol{d}}\left( {\boldsymbol{k}} \right) \times {\boldsymbol{d}}^{*} \left( {\boldsymbol{k}} \right)} \right|} \right)}$$\end{document}

so that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{k}}$$\end{document} -space locations of energy-gap zeros are defined in general by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\left. {{\boldsymbol{d}}\left( {\boldsymbol{k}} \right)} \right|^{2} \pm \left| {{\boldsymbol{d}}\left( {\boldsymbol{k}} \right) \times {\boldsymbol{d}}^{*} \left( {\boldsymbol{k}} \right)} \right| = 0$$\end{document} . Thus, although Au supports no energy-gap nodes by definition and B1u exhibits no energy-gap nodes in this model, there are numerous nodes in highly distinct k-space nodal locations for B2u and B3u. The bulk FSs have energy-gap nodal locations for B2u and B3u from Eq. (4). QPI predictions for the QSB in UTe_2_ used the Hamiltonian

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H\left( \boldsymbol{k} \right) = \left( {\begin{array}{*{20}l} {H_{{{\mathrm{UTe}}_{2} }} \left( \boldsymbol{k} \right) \otimes I_{2} } \hfill {{\Delta }_{\boldsymbol{k}} \otimes I_{4} } \\ {{\Delta }^{ + }_{\boldsymbol{k}} \otimes I_{4} } \hfill { - H_{{{\mathrm{UTe}}_{2} }}^{*} \left( { - \boldsymbol{k}} \right) \otimes I_{2} } \hfill \\ \end{array} } \right)$$\end{document}

where the order parameter is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{ \boldsymbol{k}} = {\Delta }_{0} \left( {{\boldsymbol{d}} \cdot {\boldsymbol{\sigma}}} \right)i\sigma_{2}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$I_{2} , I_{4}$$\end{document} are the unit matrices. The focus primarily was on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_{2u}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_{3u}$$\end{document} :

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{d}}_{{B_{2u} }} = \left( {C_{1} sin\left( {k_{z} c} \right), C_{0} sin\left( {k_{x} a} \right)sin\left( {k_{y} b} \right)sin\left( {k_{z} c} \right),C_{3} sin\left( {k_{x} a} \right)} \right)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{d}}_{{B_{3u} }} = \left( {C_{0} sin\left( {k_{x} a} \right)sin\left( {k_{y} b} \right)sin\left( {k_{z} c} \right),C_{2} sin\left( {k_{z} c} \right),C_{3} sin\left( {k_{y} b} \right)} \right)$$\end{document}

where a, b, c are lattice constants, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_{0} = 0,C_{1} = 300\;{\upmu} {\mathrm{eV}},C_{2} = 300\;{\upmu} {\mathrm{eV}},{\mathrm{and}}\;C_{3} = 300\;{\upmu} {\mathrm{eV}}$$\end{document} . The unperturbed bulk Green’s function is then: G0(k, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E$$\end{document} ) = [( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E$$\end{document}  + iη)I − H(k)]^−1^ (η =  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$100\;\upmu {\mathrm{eV}}$$\end{document} ) with the corresponding unperturbed spectral function: A0(k, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E$$\end{document} ) = −1/π Im G0(k, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E$$\end{document} ). The surface Green’s function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G_{\mathrm{s}} \left( {{\boldsymbol{k}},E} \right)$$\end{document} characterizes a semi-infinite system with broken translation symmetry and therefore cannot be calculated directly. A novel technique was used to model the surface using a strong planar impurity [69–71]. In the limit of an infinite impurity potential, the impurity plane splits the system into two semi-infinite spaces. So the only wavevectors in the (0–11) plane remain good quantum numbers. The effect of the planar impurity can then be exactly calculated using the T-matrix formalism, which gives one access to the surface Green’s function of the semi-infinite system. For Bogoliubov QPI predictions at the (0–11) surface of UTe_2_, a localized impurity potential \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{V} = V\tau_{z} \otimes I_{8}$$\end{document} where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V = 0.2\ {\mathrm{eV}}$$\end{document} was used to determine the surface Green’s function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{\mathrm{s}}\left( {{\boldsymbol{q}},{\boldsymbol{k}},E} \right)$$\end{document} using the T-matrix \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T\left( E \right) = \left( {I - \hat{V}\int {\frac{{d^{2} {\boldsymbol{k}}}}{{S_{BZ} }}} G_{s} \left( {{\boldsymbol{k}},E} \right)} \right)^{ - 1} \hat{V}$$\end{document} . Then, the QPI patterns for the UTe_2_ QSB are predicted directly using

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( {{\boldsymbol{q}},E} \right) = \frac{i}{2\pi }\int {\frac{{d^{2} {\boldsymbol{k}}}}{{S_{\mathrm{BZ}} }}Tr\left[ {g_{\mathrm{s}} \left( {{\boldsymbol{q}},{\boldsymbol{k}},E} \right)} \right]}$$\end{document}

where

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{\mathrm{s}} \left( {{\boldsymbol{q}},{\boldsymbol{k}},E} \right) = G_{\mathrm{s}} \left( {{\boldsymbol{q}},E} \right)T\left( E \right)G_{\mathrm{s}} \left( {{\boldsymbol{q}} - {\boldsymbol{k}},E} \right) - G_{\mathrm{s}}^{*} \left( {{\boldsymbol{q}} - {\boldsymbol{k}},E} \right)T^{*} \left( E \right)G_{\mathrm{s}}^{*} \left( {{\boldsymbol{q}},E} \right).$$\end{document}

SATM Experiments on UTe2

To explore UTe_2_ for such conjectured ITS phenomenology, single-crystal samples are introduced to a superconductive-tip scanning tunneling microscope [53–58], cleaved at 4.2 K in cryogenic ultrahigh vacuum, inserted to the scan head, and cooled to T = 280 mK. A typical topographic image \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T\left( {\boldsymbol{r}} \right)$$\end{document} of the (0–11) cleave surface as measured by a superconductive Nb tip is shown in Fig. 2b with atomic periodicities defined by vectors a*, b*, where a* = a = 4.16 Å is the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{x}$$\end{document} -axis unit-cell vector and b* = 7.62 Å is a vector in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{y}:\hat{z}$$\end{document} plane. As the temperature is reduced, a sharp zero-energy peak appears within the overall energy gap in the spectrum (Fig. 3a). This robust zero-bias \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$dI/dV|_{{{\mathrm{SIP}}}}$$\end{document} peak is observed universally, as exemplified, for example, by Figs. 3b, c. One sees that these phenomena are not due to Josephson tunneling because the zero-bias conductance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( 0 \right)$$\end{document} of Nb/UTe_2_ is many orders of magnitude larger than it could possibly be due to Josephson currents through the same junction [65], and because \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( 0 \right)$$\end{document} grows linearly with falling R before diminishing steeply as R is further reduced while \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g\left( 0 \right)$$\end{document} due to Josephson currents grow continuously as 1/R^2^. Moreover, the SIP model predicts quantitatively that such an intense \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( 0 \right)$$\end{document} peak should occur if \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} of UTe_2_ supports a QSB within the interface (Fig. 1) and because Andreev transport due to these QSB quasiparticles allows a strong zero-bias conductance to the Nb electrode.Fig. 3a Typical SIP Andreev conductance spectrum \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( V \right) \equiv dI/dV|_{{{\mathrm{SIP}}}}$$\end{document} measured with Nb scan-tip on UTe_2_ (0–11) surface for junction resistance R = 6 MΩ and T = 280 mK. A high intensity zero-bias \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$dI/dV|_{{{\mathrm{SIP}}}}$$\end{document} peak is observed. b Typical topographic image T(r) of (0–11) surface (Is = 0.2 nA, Vs = 5 mV). c Evolution of measured \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right)$$\end{document} across the (0–11) surface of UTe_2_ indicated by the arrow in b for junction resistance R = 6 MΩ and T = 280 mK. The zero-bias \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$dI/dV|_{{{\mathrm{SIP}}}}$$\end{document} peaks are universal and robust, indicating that the zero-energy ABS is omnipresent

To determine spectroscopically whether the UTe_2_ order parameter is chiral, the evolution of Andreev conductance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( V \right)$$\end{document} at T = 280 mK was measured as a function of decreasing junction resistance R or equivalently increasing tunneling matrix element \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right|$$\end{document} . Figure 4a shows the strong energy splitting \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta E$$\end{document} observable in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( V \right)$$\end{document} , that first appears and then evolves with increasing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1/R$$\end{document} . Figure 4b shows the measured \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right)$$\end{document} splitting across the (0–11) surface of UTe_2_ along the arrow indicated in Fig. 3b, demonstrating that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right)$$\end{document} split-peaks are pervasive. Decisively, we plot in Fig. 4c the measured \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta E$$\end{document} between peaks in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right)$$\end{document} at T = 280 mK versus \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1/R$$\end{document} . On the basis of predictions for energy splitting \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta E$$\end{document} within the SIP model [65] for chiral \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} and non-chiral \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} , the chiral \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} appears ruled out.Fig. 4a Measured evolution of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( V \right) \equiv dI/dV|_{{{\mathrm{SIP}}}}$$\end{document} at T = 280 mK in UTe_2_ as a function of decreasing junction resistance R and thus increasing tunneling matrix element \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right| \sim 1/R$$\end{document} . When the junction resistance falls below R ~ 5 MΩ, the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( V \right)$$\end{document} spectra start to split. b Evolution of measured \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right)$$\end{document} splitting across the (0-11) surface of UTe_2_ at junction resistance R = 3 MΩ and T = 280 mK, demonstrating that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right)$$\end{document} split-peaks are pervasive at low junction resistance R and high tunneling matrix \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| M \right|$$\end{document} . c Measured energy splitting of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( 0 \right)$$\end{document} at T = 280 mK in UTe_2_ versus \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1/R$$\end{document} . These data may be compared with predictions of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( V \right)$$\end{document} splitting for Au + iB3u and B3u order parameters [65] of UTe_2_

For QSB QPI studies, Fig. 5a shows a typical 66-nm-square field-of-view (FOV) topography of the (0–11) cleave surface, which can be studied both in the normal and superconducting states. Figure 5b shows typical \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$dI/dV$$\end{document} spectra measured with a superconductive tip in both the normal state at 4.2 K and the superconducting state at 280 mK, far below TC. In the latter case, two intense joint-coherence peaks are located at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E = {\Delta }_{{{\mathrm{Nb}}}} + {\Delta }_{{{\mathrm{UTe}}_{{2}} }}$$\end{document} . More importantly, a high density of QSB quasiparticles allows efficient creation and annihilation of Cooper pairs in both superconductors, thus generating intense Andreev differential conductance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right) \equiv dI/dV|_{A} \left( {{\boldsymbol{r}},V} \right)$$\end{document} for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| V \right| < \Delta_{_{{\mathrm{UTe}}_{{2}}}}/e \sim 300\;\upmu {\mathrm{V}}$$\end{document} as indicated by yellow shading. Compared to conventional NIS tunneling using a normal metallic tip, this Andreev conductance provides a significant improvement in the energy resolution ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta E\sim 10$$\end{document} µeV) of QSB scattering interference measurements. Comparing measured \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g\left( {{\boldsymbol{r}},V} \right)$$\end{document} : \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g\left( {{\boldsymbol{q}},V} \right)$$\end{document} recorded in the normal state at 4.2 K (Fig. 5c) with measured \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right):a\left( {{\boldsymbol{q}},V} \right)$$\end{document} in the superconducting state at 280 mK (Fig. 5d), both with identical FOV and junction characteristics, allows determination of which phenomena at the (0–11) surface emerge only due to superconductivity. Several peaks of the sextet are present in the normal state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g\left( {{\boldsymbol{q}},V} \right)$$\end{document} in Fig. 5c as they originate from scattering of the normal state band structure [68]. The complete predicted QPI sextet \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{i} :i = 1 - 6$$\end{document} are only detected in the superconducting state and appear to rely on scattering between QSB states. The sextet wavevectors are highlighted by colored arrows in Fig. 5d. The experimental maxima in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{q}},V} \right)$$\end{document} and the theoretically predicted \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{i}$$\end{document} from Ref. [68] are in excellent quantitative agreement with a maximum 3% difference between all their wavevectors. This demonstrated, for the first time, that the FS, which dominates the bulk electronic structure of UTe_2_, is also what controls QSB k-space geometry at its cleave surface. Furthermore, Fig. 5e reveals how the amplitudes of the superconducting state QPI are enhanced compared to the normal state measurements. The predominant effects of bulk superconductivity are the strongly enhanced arc-like scattering intensity connecting  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}} = 0$$\end{document}  and  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{5}$$\end{document}  and the unique appearance of wavevector  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{1}$$\end{document} .Fig. 5a Typical topographic image T(r) of the (0–11) cleave surface of UTe_2_ where QPI patterns are imaged. b Measured differential conductance in the UTe_2_ normal state g(V) at T = 4.2 K (red curve); and Andreev differential conductance in the superconducting state a(V) at T = 280 mK (blue curve). Intense Andreev conductance is observed at V = 0. c Measured g(r, 0) and g(q, 0) at T = 4.2 K in the UTe_2_ normal state in the identical FOV as a. The setpoint is Vs = 3 mV and I = 200 pA. d Measured a(r, 0) and a(q, 0) at T = 280 mK in the UTe_2_ superconducting state in the identical FOV as a and c. Here a sextet of scattering interference wavevectors qi, i = 1–6 are identified. This experimental detection of the sextet has been repeated multiple times [68]. The setpoint is Vs = 3 mV and I = 200 pA. e Relative amplitudes of the sextet wavevectors in the normal and superconducting states. Comparison of g(q, 0) linecuts at T = 4.2 K and a(q, 0) linecuts measured T = 280 mK. The linecuts are taken horizontally in the q space indicated by red arrow in d. The linecuts have been normalized by their background intensities at 280 mK and 4.2 K. The intensities of q5 and q6 are significantly enhanced in the superconducting state. Most importantly, q1 only appears in the superconducting state

To visualize the QSB dispersion k(E) of UTe_2_, we next use superconductive-tip \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right):{ }a\left( {{\boldsymbol{q}},V} \right)$$\end{document} measurements to image energy resolved QPI at the (0–11) cleave surface. Figure 6a presents the measured \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{r}},V} \right)$$\end{document} at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| V \right| = 0\ {\upmu \mathrm{V}}, 50\ {\upmu \mathrm{V}}, 100\ {\upmu \mathrm{V}},150\ {\upmu \mathrm{V}},200\ {\upmu \mathrm{V}},250\ {\upmu \mathrm{V}}$$\end{document} recorded at T = 280 mK in the identical FOV as Fig. 5a. These data are highly typical of such experiments in UTe_2_ [68]. Figure 6b contains the consequent scattering interference patterns \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( {{\boldsymbol{q}},V} \right)$$\end{document} at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| V \right| = 0\ {\upmu \mathrm{V}}, 50\ {\upmu \mathrm{V}}, 100\ {\upmu \mathrm{V}},150\ {\upmu \mathrm{V}},200\ {\upmu \mathrm{V}},250\ {\upmu \mathrm{V}}$$\end{document} as derived by Fourier analysis of Fig. 6a. Here the energy evolution of scattering interference of the QSB states is obvious. For comparison with theory, detailed predicted characteristics of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( {{\boldsymbol{q}},E} \right)$$\end{document} for a B2u-QSB and B_3u*-*QSB at the (0–11) SBZ were determined in Ref. [68]; here again energies range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| E \right| = 0\;{\upmu} {\mathrm{eV}}, 50\;{\upmu} {\mathrm{eV}}, 100\;{\upmu} {\mathrm{eV}},150\;{\upmu} {\mathrm{eV}},200\;{\upmu} {\mathrm{eV}},250\;{\upmu} {\mathrm{eV}}$$\end{document} . Each QPI wavevector is determined by maxima in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\left( {{\boldsymbol{q}},E} \right)$$\end{document} QPI pattern (Fig. 6b); these phenomena are highly repeatable in multiple independent experiments. The strongly enhanced QPI features occurring along the arc connecting \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}} = 0$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{5}$$\end{document} (Fig. 6b) are characteristic of the theory for a B3u*-*QSB [68]. Most critically, the intense QPI appearing at wavevector \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{1}$$\end{document} (yellow circle in Fig. 6b) is a characteristic of the B3u_ superconducting state, deriving from its geometrically unique nodal structure [68]. The appearance of scattering interference of QSB quasiparticles at q1 in the superconducting state (Figs. 5d and 6b) is precisely as would be anticipated in theory [51, 52] due to projection of B3u energy-gap nodes on the bulk FS [68] onto the (0–11) crystal surface 2D Brillouin zone.Fig. 6a Measured a(r, V) at the (0–11) cleave plane of UTe_2_ at bias voltages \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| V \right| = 0\;{\upmu} {\mathrm{V}}, 50\;{\upmu} {\mathrm{V}}, 100\;{\upmu} {\mathrm{V}},150\;{\upmu} {\mathrm{V}},200\;{\upmu} {\mathrm{V}},250\;{\upmu} {\mathrm{V}}$$\end{document} . The setpoint is Vs = 3 mV and I = 200 pA. b Measured a(q, V) at the (0–11) cleave plane of UTe_2_ at bias voltages \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| V \right| = 0\;{\upmu} {\mathrm{V}}, 50\;{\upmu} {\mathrm{V}}, 100\;{\upmu} {\mathrm{V}},150\;{\upmu} {\mathrm{V}},200\;{\upmu} {\mathrm{V}},250\;{\upmu} {\mathrm{V}}$$\end{document} . The setpoint is Vs = 3 mV and I = 200 pA. Each QPI wavevector in this FOV, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{1}$$\end{document} (yellow), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{2}$$\end{document} (brown) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{5}$$\end{document} (white), is identified as the maxima position (colored circles) in the QPI data. Particularly \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{1}$$\end{document} is a characteristic only of the B3u superconducting state [68], and it only exists inside the energy gap

Conclusions

Overall, the chiral order parameters Au + iB1u and B3u + iB2u proposed for UTe_2_ appear inappropriate because of the observed Andreev conductance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( 0 \right)$$\end{document} splitting when reducing the Nb/UTe_2_ separation [65]. Within the four possible odd-parity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta }_{{\boldsymbol{k}}}$$\end{document} symmetries Au, B1u, B2u and B3u, the isotropic Au order parameter also appears insupportable because its QSB is a Majorana-cone of Bogoliubons with zero density-of-states at zero energy, meaning that Andreev conductance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( 0 \right)$$\end{document} would be highly suppressed. Andreev conductance between Nb (s-wave) and UTe_2_ (putative p-wave) superconductors allows visualization of a powerful zero-energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a\left( V \right) = dI/dV|_{{{\mathrm{SIP}}}}$$\end{document} peak at the UTe_2_ (0–11) surface. And, with enhanced tunneling to an s-wave electrode (Nb) this zero-energy Andreev spectrum splits strongly into two finite-energy conductance maxima [65]. Moreover, visualizing dispersive QSB scattering interference reveals unique in-gap QPI patterns exhibiting a characteristic sextet of wavevectors \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{{\boldsymbol{i}}} :i = 1 - 6$$\end{document} due to projection of the bulk superconductive band structure onto the (0–11) surface [68]. Although \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{2}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{6}$$\end{document} are weakly observable in the normal state, features at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{5}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{6}$$\end{document} become strongly enhanced for superconducting state QPI at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| E \right| < {\Delta }$$\end{document} and QPI appears at wavevector \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{q}}_{1}$$\end{document} uniquely in the superconducting state. This complete phenomenology, by correspondence with theory [65, 68], is most consistent with a 3D, odd-parity, spin-triplet, time-reversal-symmetry conserving, a-axis nodal superconducting order parameter with B3u symmetry in UTe_2_.