## Signature of point nodal superconductivity in the Dirac semimetal PdTe

C. S. Yadav , 1 ∗ Sudeep Kumar Ghosh , 2 † Pankaj Kumar , A. Thamizhavel , and P.L. Paulose 1 3 3 1 School of Physical Sciences, Indian Institute of Technology Mandi, Kamand, Mandi-175075 (H.P.) India 2 Department of Physics, Indian Institute of Technology Kanpur, Kanpur- 208016 (U.P.), India and 3 DCMP &amp; MS, Tata Institute of Fundamental Research, Mumbai -400005 (Maharashtra) India

(Dated: August 20, 2024)

Recent Angle-Resolved Photo-emission Spectroscopy (ARPES) experiments [Phys. Rev. Lett. 130, 046402 (2023)] on PdTe, a 3D-Dirac semimetal and a superconductor with the transition temperature T c ∼ 4.3 K, have revealed compelling evidence of the presence of bulk nodes in the superconducting order parameter. To investigate the validity of this proposition, here we present a detailed investigation of the magnetic field dependence of the specific heat of PdTe down to temperatures ∼ 58 mK. We observed that the low temperature specific heat of PdTe with an externally applied magnetic field exhibits a power-law field dependence, a characteristic of unconventional superconductivity. Furthermore, the zero-field low-temperature electronic specific heat follows a cubic temperature dependence, which is a signature of the presence of bulk point nodes in PdTe. These intriguing observations suggest that PdTe is a rare and fascinating topological material that exhibits both Dirac semimetallic properties and superconductivity with point nodal gap symmetry.

## I. INTRODUCTION

Topological semimetals represent a fascinating class of materials characterized by topologically nontrivial bulk band crossings and robust surface states that exhibit unique electronic properties, making them highly desirable for practical applications [1]. The realm of topological semimetals becomes even more captivating when combined with superconductivity, as it opens up the possibility of realizing intrinsic topological superconductivity. Recent discoveries of materials exhibiting both topological semimetallic behaviour and superconductivity have further bolstered this research direction, providing new avenues for exploring and understanding the intriguing interplay between topology and superconductivity. Examples include, superconducting 3D Dirac semimetals (Ta,Nb)OsSi [2], 3D Weyl semimetals LaPt(Si,Ge) [3], Z 2 TMs (Cs,Rb,K)V Sb 3 5 [4, 5] and ZrOsSi [6], Weyl nodal-line SMs La(Ni,Pt)Si [7] and Kramers nodal-line SMs (Ti, Nb, Hf, Ta)RuSi [3].

state with surface Bogoliubov-Dirac cones and a secondorder time-reversal invariant topological superconductor with helical Majorana hinge modes [13]. Recent angleresolved photo-emission spectroscopy (ARPES) measurements [14] provide compelling evidence that the chalcogenide superconductor PdTe exemplifies such intriguing material characteristics as well.

Among the transition metal chalcogenides, FeSe shows unconventional superconductivity below ∼ 8 2 . K. The Scanning Tunneling Microscopy (STM) experiments on FeSe crystalline films suggested the existence of nodal superconductivity in stoichiometric FeSe [8]. Furthermore, the identification of C 4 rotation-symmetry protected Dirac semimetallic behavior, featuring distinctive surface Dirac cones, in the Fe(Se,Te) superconductors [9, 10] has significantly advanced research in this class of materials. Several theoretical investigations [11-13] on these compounds have emphasized the potential for realizing various topological phases in Dirac semimetal superconductors. These phases include topological crystalline superconductivity [11], nodal higher-order topological superconductivity [12] and, novel gapped superconducting

∗ shekhar@iitmandi.ac.in

† skghosh@iitk.ac.in

PdTe is a hexagonal superconductor with T c ∼ 4 3 K . [15-17] that has a layered structure similar to that of the tetragonal Fe chalcogenide superconductors. However, the 4 d -orbital character of Pd and the face-shared octahedral environment suppress electronic correlations and promote three-dimensionality in PdTe [18, 19]. This enhanced three-dimensionality, in turn, gives rise to weaker Fermi surface nesting and diminished long-range correlations [18]. Based on first principles calculations, it was proposed that there is a diminished charge transfer energy between Pdd and Tep states in PdTe when compared to Fe chalcogenides. The Pdd states in PdTe are almost fully occupied, leaving no states available to form a large Hund's moment and there is no long-ranged magnetic correlations present in PdTe [18].

In recent ARPES measurements on PdTe [14], a noteworthy revelation emerged that a bulk Dirac point is situated just below the Fermi level accompanied by Fermi arc surface states on the (010) surface. These investigations also uncovered compelling evidence of a rare phenomenon: the coexistence of nodeless superconductivity on the surface and nodal superconductivity in the bulk of PdTe. Inspired by the distinctive characteristics exhibited by this material, here we embarked on comprehensive field- and temperature-dependent investigations of the heat capacity of PdTe, down to ultra-low temperatures ( ∼ 58mK). We found that the zero-field electronic specific heat data at these low temperatures in the superconducting state aligns remarkably well with the Bardeen-Cooper-Schrieffer (BCS)-type model featuring a p-wave gap symmetry with point nodes in the bulk. Furthermore, our findings gain additional support from

/s40/s101/s41

/s32

/s32

/s32

Figure 1. Characterization of PdTe samples : (a) An image of the PdTe ingot showing the shiny nature of crystallites. (b) XRD pattern of powdered PdTe sample showing a clean hexagonal phase (space group: P6 3 /mmc ). (c) The temperature dependence of the electrical resistivity. The inset shows a zoomed region around the superconducting transition (1.8 ≤ T ≤ 5 K). (d) DC magnetic susceptibility as a function of temperature at H = 20 Oe. The inset shows H c 1 ( T ) obtained from the M -H isotherms, and the fit with H c 1 ( T ) = H c 1 (0)(1 -( T/T C ) 2 ). (e) AC magnetic susceptibility measured at different external hydrostatic pressures p = 0, 4, 6, and 9 kbar. (f) Total specific heat C(T) fitted with the Debye model (Eq. (1)).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00216v2-figure-762d0d01-8647-4693-86b9-700a55643ea8.jpg)

the magnetic field dependence of the electronic heat capacity, which distinctly points toward the unconventional nature of superconductivity in PdTe.

## II. RESULTS AND DISCUSSION

PdTe samples were prepared by solid-state reaction method using high purity ( ≥ 99 9%) Pd and Te pow-. ders inside an evacuated quartz tubes at 1080 C for 24 0 hours. The ingot were found to consist of large crystallites of 1-2 mm size as shown in Fig. 1(a). The Rietveld refinement of the X-ray diffraction (XRD) data shown in Fig. 1b revealed the formation of a clean NiAs-type hexagonal phase in PdTe with the space group symmetry P6 3 /mmc (No. 194) and with the lattice parameters, a = 4.15 ˚ and c = 5.67 ˚ . A A

is observed at T ∼ 4 3 K (inset of Fig. 1(c)), . with full diamagnetic signal (4 πχ = -1) as shown in Fig. 1(d). We note that ρ T ( ) follows normal Fermi liquid behaviour ( ρ T ( ) ∝ T ) up to temperatures 2 ∼ 42 K, which is more than the 1/5 th of the Debye-temperature (Θ D ∼ 200 K) of PdTe. The inset of Fig. 1(d) shows temperature dependence of the lower critical field ( H c 1 ( T )) obtained from the M -H isotherms. Fitting of H c 1 ( T ) curve with the empirical relation H c 1 ( T ) = H c 1 (0)(1 -( T/T C ) 2 ) gives H c 1 (0) ≈ 165 ± 10 Oe, which is slightly lower than the reported value ( ≈ 200 Oe) [17].

The electrical characterizations (under an ac current 10 mA at frequency 79 Hz) and the heat capacity measurements were performed using a Physical Properties measurement System (PPMS), and the magnetic measurements were performed using a Magnetic property measurement system (SQUID-MPMS). The temperature dependence of the electrical resistivity ρ T ( ) of PdTe is shown in Fig. 1(c). A sharp superconducting transition

The AC susceptibility data measured at different externally applied hydrostatic pressures ( p ) = 0, 4, 6 and 8 kbar are shown in Fig. 1(e). We note that the T c of PdTe gets suppressed with increasing pressure (∆T /∆ c p ∼ -0 01 K/kbar). . The low value of the coefficient of the suppression of T c of PdTe with p may imply that its density of states (DOS) is nearly flat at the Fermi energy. This is consistent with the band structure calculations showing that the Fermi energy is located very close to the local maxima of the DOS [18]. The decrease in T with pressure for PdTe contrasts with unusually large c p -dependence on T c shown by the Fe chalcogenide superconductors [20]. Note that the high-pressure AC susceptibility measurements were conducted using a pressure

cell, making it challenging to calculate the sample's exact demagnetization factor. While we couldn't estimate the precise superconducting volume fraction under pressure, the data clearly indicates that it decreases as pressure increases.

The variation of the total heat capacity (C) in zero magnetic field with temperature for PdTe is shown in shown in the Fig. 1(f). We note that in the normal state at high temperatures the specific saturates to C(T = 250 K) ∼ 50 J/mole-K which is close to classical DulongPetit limit of C = 3 nR = 6 R = 49 9 J/mole-K due to . the lattice vibrations, as expected. The low-T normal state heat capacity data is fitted with a combination of the electronic and lattice contributions (Debye model of specific heat) using the expression [21]

$$C = \gamma T + 9 n R \left ( \frac { T } { \Theta _ { D } } \right ) ^ { 3 } \int _ { 0 } ^ { y } \frac { x ^ { 4 } e ^ { x } } { ( e ^ { x } - 1 ) ^ { 2 } } d x \quad ( 1 ) \quad \overline { \Xi } \,.$$

where y =Θ D / Tand Θ D is the Debye Temperature. The best possible fit gives the Sommerfeld coefficient γ ∼ 6 32 . mJ/mole-K 2 and Θ D ∼ 200 K. The values of γ in normal state, obtained from the low Temperature linear fit of C/T vs T 2 data is also ∼ 6 30 mJ/mole-K . The DOS . 2 at Fermi energy N(E F ) obtained using the relation γ = π 2 k 2 B N(E F )/3 is N(E F ) ∼ 2 67 states/eV-f.u. . This DOS N(E F ) contains the enhancement factor due to electronphonon coupling constant λ el -ph , and can be related to the bare-band DOS at the Fermi level N (E 0 F ) by the expression N(E F ) = N (E 0 F )(1 + λ el -ph ). The calculated value of the λ el -ph using the MacMillan's theory, corresponding to Θ D = 200 K is λ el -ph ∼ 0 7, . which is similar to that of a weak coupling BCS superconductor [22]. However, considering the modified Macmillan formula with the logarithmic averaged phonon frequency it was reported that λ el -ph ≈ 1.4 in the Ref. [15].

Further, we estimate the Fermi velocity v F in PdTe from the N(E F ) using a single band free electron gas model as v F = π 2 ℏ 3 m ∗ 2 V f . u . N E ( F ), where V f.u. =V cell /2 [21]. Assuming m ∗ = m , we have e v F = 4 8 . × 10 8 cm/sec. The mean free path l = v F τ , ( τ = mean relaxation time) obtained using the formula l = 3 π 2 ( ℏ e 2 ρ 0 )( ℏ m v ∗ F ) 2 is estimated to be l ∼ 43.8 nm [21]. This is an order of magnitude smaller than l = 305.4 nm (estimated using l = (2m e v F )/(ne 2 ρ 0 ), where n is carrier concentration obtained from Hall measurement and ρ 0 is the resistivity at T = 2 K [15]) corresponding to the v F ∼ 7 × 10 7 cm/sec estimated by considering the bare band density of states obtained from the band structure calculation reported in the Ref. [15].

Now, we focus on the low temperature ( T ∼ 0 058 . -6 K) behavior of the specific heat C(T) in externally applied magnetic fields. The variation of the specific heat in this regime as a function of temperature for different magnetic fields up to ∼ 1200 Oe are shown in the Fig. 2(a). The gradual suppression of the peak height and of T c in the presence of magnetic field is clearly evident from this figure. The variation of the upper critical field H c 2 (T) with temperature is shown in the inset of Fig. 2(a). We have fitted this H c 2 (T) vs T data using the generalized Ginzburg-Landau (GL) equation: H c 2 (T) = H c 2 (0)[(1-t 2 )/(1+t )] [23-25] where t = 2 T/T c is the reduced temperature, the zero-temperature

/s32

/s50

Figure 2. Field dependence of specific heat : (a) Variation of the total specific heat (C) with temperature at various magnetic fields: H = 0 - 1200 Oe, The inset shows the temperature of the upper critical field H c 2 . The dashed red line in the inset shows a fitting of the H c 2 data using the generalized Ginzburg-Landau equation: H c 2 (T) = H c 2 (0)[(1t 2 )/(1+t )] [23-25]. 2 (b) Magnetic field dependence of the relative electronic specific heat ∆ γ given by Eq. (2). The fittings shown are with the form ∆ γ ∝ ( H -H 0 ) z in the mixed state of PdTe for H &gt; H 0 with H 0 ≈ H c 1 (0). The best is fit found with the exponent z = z 0 = 0 644 . ± 0 059 which is clearly . different from the linear ( z = 1) and square-root ( z = 0 5) . dependences shown.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00216v2-figure-2a477cec-8fd7-44c9-9b32-cb151f63fa70.jpg)

/s32

upper critical field is H c 2 (0), and T c is 4.47 K. The value of H c 2 (0) as obtained from the fit is H c 2 (0) = 1473.16 Oe. Further, we have extracted the zerotemperature orbital limited upper critical field H orb c 2 (0) from the Werthamer-Helfand-Hohenberg formula given as µ H 0 c 2 (0) = -0 693 . T c ( µ dH 0 c 2 /dT ) T = T c [26]. The calculated value of µ H 0 c 2 (0) is 1332 Oe, which is lower than 1473.155 Oe, obtained from the GL fit. The coherence length ξ calculated using ξ = ( ϕ 0 / (2 π µ 0 H c 2 )) 1 / 2 ∼ 47 2 . nm which is close to the mean free path l ∼ 43 8 . nm of PdTe indicating that PdTe is close to the dirty limit [27, 28].

We have extracted the relative electronic specific heat defined as:

$$\Delta \gamma ( H ) = \lim _ { T \rightarrow 0 } \left [ \frac { C ( H, T ) } { T } - \frac { C ( H = 0, T ) } { T } \right ] \quad ( 2 ) \quad \text{iter} \quad \stackrel { \tilde { k } } { \text{isort} } \quad.$$

from the low temperature specific heat data (shown in Fig. 2(a)) at different magnetic fields. The field dependence of ∆ γ H ( ) for type-II superconductors in the mixed state is expected to be ∝ H for a conventional s -wave BCS gap and a power law ∝ H z with 0 5 . &lt; z &lt; 1 for an unconventional gap symmetry in general [29-32] with z = 0 5 usually implying the presence of line nodes in the . order parameter for a two dimensional material [29, 30]. The variation of ∆ γ H ( ) with the applied magnetic fields in the mixed state of the PdTe superconductor is shown in Fig. 2(b). We have fitted the data for ∆ γ H ( ) in the mixed state for 200 ≲ H ≲ 1200 Oe by a power law ∼ ( H -H 0 ) z for H &gt; H 0 where H 0 ≈ H c 1 (0). We find a good fit with the exponent z = 0 644 . ± 0 059 as . shown in Fig. 2(b). Such a power law field dependence of ∆ γ H ( ) clearly points to the unconventional nature of superconductivity in PdTe consistent with previously reported results [17]. However, to clearly determine the nature and position of the nodes in the superconducting gap, orientation-dependent measurements of the low temperature specific heat using good quality single crystals are needed [33, 34].

The low temperature total specific heat data C total ( T ) in the zero field of PdTe in the normal state ( T c ≤ T ≤ 10K) is fitted with the form C n ( T ) = γ N T + C ph , where the first term is the electronic contribution to the specific heat in the normal state and the second term is the phononic contribution to the specific heat C ph = β T 3 + δ T . The fitting shown in the inset of Fig. 3(a) gives 5 γ N = 6.303 mJ/mole-K , 2 β = 0.748 mJ/mole-K 4 and δ = 0.0019 mJ/mole-K 6 for PdTe. Then we extract the low temperature electronic specific heat (C el ) data of PdTe by defining C el ( T ) = C total ( T ) -C ph ( T ). The variation of the C el (T) at zero-field ( H = 0) at low temperatures is shown in the Fig. 3. The low temperature behaviour of C el can be used to investigate symmetry of the superconducting gap realized in a particular material. To this end, we have defined γ e ( T ) = C el /T and fitted it with a generalized single-band BCS-type form applicable for different types of superconducting gap symmetries given by [35]:

$$\stackrel { ( \mathbf v ) } { z e r o - } } _ { \stackrel { [ o r b } { c 2 } ( 0 ) } } \gamma _ { e } = \frac { N ( E _ { F } ) } { 2 k _ { B } T ^ { 3 } } \left \langle \int _ { 0 } ^ { \infty } d \epsilon \, \text{sech} ^ { 2 } \left ( \frac { E _ { k } } { k _ { B } T } \right ) \left [ E _ { k } ^ { 2 } - \frac { T } { 2 } \frac { d \Delta _ { k } ^ { 2 } ( T ) } { d T } \right ] \right \rangle _ { F S } \, ( 3 )$$

where E k = √ ϵ 2 +∆ ( 2 k T ) and ⟨⟩ FS represents an average over the Fermi surface (assumed to be spherical). The gap function (∆ ( k T )) defined as ∆ ( k T ) = ∆ ( ′ T g ) ( k ) has a universal temperature dependence ∆ ( ′ T ) = ∆ tanh 0 [ 1 82 . { 1 018 (T . c /T -1) 0 51 . }] and all the information about its symmetry is contained in the function g ( k ) = g θ, ϕ ( ) with θ and ϕ being azimuthal and polar angles respectively [36]. Then for a particular gap symmetry we would have only one fitting parameter ∆ 0 k B T c in Eq. (3). In an s-wave superconductor, the isotropic gap at the Fermi surface leads to the exponential decay of the electronic specific heat on decreasing temperature. However, this behaviour of the electronic specific heat gets affected by the presence of nodes in the superconducting gap at the Fermi surface and for T &lt; T c a power-law decay of electronic specific heat is observed that goes as ∝ T 2 for line nodes and ∝ T 3 in the presence of point nodes [37-39].

We have fitted the experimental data of C el for PdTe by the form given in Eq. (3) assuming a spherical Fermi surface, for three different types of models corresponding to the superconducting gap symmetries: (i) s -wave with g ( k ) = 1 having a full gap, (ii) p -wave with g ( k ) = sin( θ ) having only point nodes, and (iii) d -wave with g ( k ) = sin(2 θ ) having a line node at the equator ( θ = π/ 2) that dominates the low temperature behaviour [36]. We first make sure that the jump at T c in the electronic specific heat is accurately reproduced from all the fitted models. The fittings of C /T el vs. T, which enables us to look into the low temperature region carefully, are shown in Fig. 3(a). We note from this Fig. 3(a) that the specific heat behaves as a power-law expected for a nodal order parameter instead of being exponentially suppressed expected for a fully gapped order parameter as a function of temperature in the low temperature region 0 ≤ T/ T c ≲ 1 / 3. As a result, the s -wave does not give good fit to the data. The d -wave model also gives a poor fitting while a reasonably good fit to the experimental electronic specific heat data throughout the temperature range was obtained for the p -wave superconducting gap function. To exemplify the low temperature behavior of the specific heat we have also shown a zoomed version of Fig. 3(a) in Fig. 3(b). The values of the fitting parameter ∆ 0 k B T c for the s -, p -and d -wave models shown in Fig. 3 are ∼ 1 83, . ∼ 2 29 . and ∼ 2 56 . respectively. Hence, we conclude that the electronic specific heat of PdTe is well described by a weak-coupling p -wave generalized BCS-type order parameter with ∆ 0 k B T c ∼ 2 29 . which is close to the conventional weak-coupling BCS value ∼ 1 76. . This conclusion is further validated by explicitly fitting the low temperature behaviour of the specific heat with temperature power laws expected for

Figure 3. Zero-field electronic specific heat : (a) The electronic specific heat (C el ) data as a function of temperature without any external field present. The experimental data is fitted with three different models of the order parameter: a fully gapped s -wave BCS (pink), a point nodal p -wave (red) and a line-node dominated d -wave (blue). The inset shows behavior of the total specific heat ( C total ), the total specific heat in the normal state ( C n ) extrapolated to zero temperature and the electronic specific heat ( C el ) obtained after subtracting the phonon contribution from C total . (b) The zoom of the marked low temperature region of the specific heat figure in (a). We note that clearly the p -wave model fits the specific heat data much better throughout the temperature window than the other two models. (c) The low temperature specific heat shows a clear ∼ T 3 temperature dependence characteristic of the presence of bulk point nodes as in the p -wave model. Here, the values of the constant parameters used in the fitting are: a = 15 49 and . b = 6 51. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00216v2-figure-5f3eac2a-6d1d-4e7d-9677-36f406405e80.jpg)

the nodal order parameters as shown in the Fig. 3(c). We note that there is a clear ∼ T 3 temperature dependence in the region 0 ≤ T/ T c ≤ 1 / 3 characteristic of the presence of bulk point nodes in the superconducting order parameter of PdTe. This C el ∝ T 3 dependence at low temperatures for PdTe is quite remarkable and conforms with the evidence of the bulk point nodes in the superconducting order parameter reported based on ARPES measurements [14].

PdTe is a multiband superconductor with the point group D 6 h which has both non-degenerate and degen- erate irreducible representations. As a result, several symmetry allowed singleand multi-component superconducting order parameters are allowed in PdTe and in general it is difficult to clearly identify which superconducting instability channel is realized in PdTe. However, mean-field calculations using a Slater-Koster type tightbinding model for PdTe, considering only the p x and p y orbitals of Te and the effect of spin-orbit coupling such that the Type-I Dirac crossing observed in ARPES along the Γ -A direction is qualitatively captured, also suggest that a superconducting gap function with a bulk point node is indeed possible in PdTe even in the most symmetric superconducting instability channel [14] .

The jump in the electronic specific heat (∆ C ) at T c in zero field ( H = 0) for PdTe shown in the Fig. 3(a) is ∆ C/ γ ( N T c ) ≈ 1 49 which is also similar to the weak . coupling BCS value ∼ 1 43 [22]. . Our finding is consistent with previous studies ∆ C/ γ ( N T c ) ≈ 1 33 [17] . and ≈ 1 67 [15] but a little less compared to the value . ≈ 2 1 . reported in the Ref. [19]. The value of ∆ C/ γ ( N T c ) can, however, depend on the sample quality and the superconducting volume fraction in the samples. Moreover, some superconductors such as many of the cuprate and Fe-chalcogenide superconductors do not exhibit jump in the specific heat as well [40-42]. We note that the electronic specific heat data of PdTe was reported to be fitted with an unconventional two-gap s-wave type model in the Ref.[19] where the corresponding samples had atleast ∼ 16% non-superconducting volume fraction. We also note that another recent work on PdTe reported multigap behaviour from thermal conductivity data [43]. However, the thermal conductivity data may not be able to properly determine the nodal superconducting behaviour due to the corresponding small size and irregular shapes of PdTe single crystals used [43].

## III. SUMMARY AND CONCLUSION

Nodal superconducting materials [36, 44-46] represent a rare and intriguing class of unconventional superconductors, captivating the interest of researchers seeking to unravel their nature and pairing mechanisms. Based on recent angle-resolved photo-emission spectroscopy (ARPES) measurements [14] it was proposed that PdTe, a 3D Dirac semimetal in its normal state, also exhibits point nodal superconductivity in the bulk. To verify the unconventional nature of superconducting order parameter in PdTe, we have conducted extensive fieldand temperature-dependent specific heat measurements in the mixed state of PdTe.

In our field dependent specific heat data, we have observed a power-law field dependence in the specific heat that deviates from the usual linear behavior expected for an isotropic fully gapped order parameter [32]. Our zero-field specific heat data offers additional insights, remarkably revealing a clear ∼ T 3 temperature dependence in the low-temperature regime ( T &lt; T c / 3) that

strongly suggests the presence of bulk point nodes in the superconducting order parameter of PdTe. Notably, a weak-coupling BCS-type superconducting order parameter with point nodes yields an excellent fit to the zero-field specific heat data throughout the temperature range. Thus, our bulk specific heat measurements in the superconducting state of the 3D Dirac semimetal PdTe point towards its rare identification as a point nodal superconductor as proposed in the Ref. [14].

PdTe stands out as a unique 3D Dirac semimetal material [14] because it: i) features Dirac points near the Fermi level, ii) has an accessible transition temperature ( T c ), iii) is cleavable to measure superconducting Fermi arc states, and iv) has unconventional point nodal superconducting ground state. Further confirmation of the point nodal superconducting order parameter in PdTe will require the use of complementary experimental techniques, such as solid-state nuclear magnetic

- [1] B. Q. Lv, T. Qian, and H. Ding, Experimental perspective on three-dimensional topological semimetals, Rev. Mod. Phys. 93 , 025002 (2021).
- [2] S. K. Ghosh, P. K. Biswas, C. Xu, B. Li, J. Z. Zhao, A. D. Hillier, and X. Xu, Time-reversal symmetry breaking superconductivity in three-dimensional Dirac semimetallic silicides, Phys. Rev. Res. 4 , L012031 (2022).
- [3] T. Shang, J. Zhao, L.-H. Hu, J. Ma, D. J. Gawryluk, X. Zhu, H. Zhang, Z. Zhen, B. Yu, Y. Xu, Q. Zhan, E. Pomjakushina, M. Shi, and T. Shiroka, Unconventional superconductivity in topological kramers nodalline semimetals, Sci. Adv. 8 , eabq6589 (2022).
- [4] K. Jiang, T. Wu, J.-X. Yin, Z. Wang, M. Z. Hasan, S. D. Wilson, X. Chen, and J. Hu, Kagome superconductors AV3Sb5 (A = K, Rb, Cs), Natl. Sci. Rev. 10 , nwac199 (2022).
- [5] Z. Shan, P. K. Biswas, S. K. Ghosh, T. Tula, A. D. Hillier, D. Adroja, S. Cottrell, G.-H. Cao, Y. Liu, X. Xu, Y. Song, H. Yuan, and M. Smidman, Muon spin relaxation study of the layered kagome superconductor CsV Sb , Phys. Rev. Res. 3 5 4 , 033145 (2022).
- [6] S. K. Ghosh, B. Li, C. Xu, A. D. Hillier, P. K. Biswas, X. Xu, and T. Shiroka, ZrOsSi: a Z 2 topological metal with a superconducting ground state, Front. Phys. 11 , 1256166 (2023).
- [7] T. Shang, S. K. Ghosh, M. Smidman, D. J. Gawryluk, C. Baines, A. Wang, W. Xie, Y. Chen, M. O. Ajeesh, M. Nicklas, et al. , Spin-triplet superconductivity in weyl nodal-line semimetals, npj Quantum Mater. 7 , 35 (2022).
- [8] C.-L. Song, Y.-L. Wang, P. Cheng, Y.-P. Jiang, W. Li, T. Zhang, Z. Li, K. He, L. Wang, J.-F. Jia, H.-H. Hung, C. Wu, X. Ma, X. Chen, and Q.-K. Xue, Direct Observation of Nodes and Twofold Symmetry in FeSe Superconductor, Sci. 332 , 1410 (2011).
- [9] P. Zhang, K. Yaji, T. Hashimoto, Y. Ota, T. Kondo, K. Okazaki, Z. Wang, J. Wen, G. D. Gu, H. Ding, and S. Shin, Observation of topological superconductivity on the surface of an iron-based superconductor, Sci. 360 ,

resonance (NMR) to measure spin-susceptibility [47] and tunnel-diode oscillator (TDO) to measure the London penetration depth [48] in combination with orientationdependent measurements with good quality single crystals [33, 34].

## IV. ACKNOWLEDGMENT

CSY acknowledges SERB-DST (India) for the CRG grant (CRG/2021/002743). Advanced Material Research Center (AMRC), IIT Mandi is acknowledged for the experimental facility. SKG gratefully acknowledges financial support from SERB, Government of India via the Startup Research Grant: SRG/2023/000934, as well as from IIT Kanpur through the Initiation Grant (IITK/PHY/2022116). The authors also acknowledge Dibyendu Samanta for discussions.

## 182 (2018).

- [10] P. Zhang, Z. Wang, X. Wu, K. Yaji, Y. Ishida, Y. Kohama, G. Dai, Y. Sun, C. Bareille, K. Kuroda, et al. , Multiple topological states in iron-based superconductors, Nat. Phys. 15 , 41 (2019).
- [11] T. Kawakami and M. Sato, Topological crystalline superconductivity in Dirac semimetal phase of iron-based superconductors, Phys. Rev. B 100 , 094520 (2019).
- [12] Z. Wu and Y. Wang, Nodal higher-order topological superconductivity from a C -symmetric Dirac semimetal, 4 Phys. Rev. B 106 , 214510 (2022).
- [13] M. Kheirkhah, Z.-Y. Zhuang, J. Maciejko, and Z. Yan, Surface Bogoliubov-Dirac cones and helical Majorana hinge modes in superconducting Dirac semimetals, Phys. Rev. B 105 , 014509 (2022).
- [14] X. P. Yang, Y. Zhong, S. Mardanya, T. A. Cochran, R. Chapai, A. Mine, J. Zhang, J. S´nchez-Barriga, Z.-J. a Cheng, O. J. Clark, et al. , Coexistence of bulk-nodal and surface-nodeless Cooper pairings in a superconducting Dirac semimetal, Phys. Rev. Lett. 130 , 046402 (2023).
- [15] A. B. Karki, D. A. Browne, S. Stadler, J. Li, and R. Jin, PdTe: A strongly coupled superconductor, J. Condens. Matter Phys. 24 , 055701 (2012).
- [16] A. B. Karki, V. O. Garlea, R. Custelcean, S. Stadler, E. Plummer, and R. Jin, Interplay between superconductivity and magnetism in Fe 1 -x Pd Te, Proc. Natl. Acad. x Sci. U.S.A. 110 , 9283 (2013).
- [17] B. Tiwari, R. Goyal, R. Jha, A. Dixit, and V. P. S. Awana, PdTe: a 4.5 K type-II BCS superconductor, Supercond. Sci. Technol. 28 , 055008 (2015).
- [18] C. E. Ekuma, C.-H. Lin, J. Moreno, W. Ku, and M. Jarrell, First-principles Wannier function analysis of the electronic structure of PdTe: weaker magnetism and superconductivity, J. Condens. Matter Phys. 25 , 405601 (2013).
- [19] R. Chapai, P. S. Reddy, L. Xing, D. E. Graf, A. B. Karki, T.-R. Chang, and R. Jin, Evidence for unconventional superconductivity and nontrivial topology in pdte, Sci.

Rep. 13 , 6824 (2023).

- [20] D. C. Johnston, The puzzle of high temperature superconductivity in layered iron pnictides and chalcogenides, Adv. Phys. 59 , 803 (2010).
- [21] C. Kittel, Introduction to Solid State Physics , 8th ed. (Wiley, 2004).
- [22] W. L. McMillan, Transition Temperature of StrongCoupled Superconductors, Phys. Rev. 167 , 331 (1968).
- [23] K. Maki, Effect of pauli paramagnetism on magnetic properties of high-field superconductors, Phys. Rev. 148 , 362 (1966).
- [24] N. R. Werthamer, E. Helfand, and P. C. Hohenberg, Temperature and purity dependence of the superconducting critical field, H c 2 . iii. electron spin and spin-orbit effects, Phys. Rev. 147 , 295 (1966).
- [25] L. Fang, Y. Wang, P. Y. Zou, L. Tang, Z. Xu, H. Chen, C. Dong, L. Shan, and H. H. Wen, Fabrication and superconductivity of na Tas x 2 crystals, Phys. Rev. B 72 , 014534 (2005).
- [26] N. Momono and M. Ido, Evidence for nodes in the superconducting gap of La 2 -x Sr x CuO - T 4 2 dependence of electronic specific heat and impurity effects, Phys. C: Supercond. 264 , 311 (1996).
- [27] C. S. Yadav and P. L. Paulose, Upper critical field, lower critical field and critical current density of FeTe 0 60 . Se 0 40 . single crystals, New J. Phys. 11 , 103046 (2009).
- [28] A. Pourret, L. Malone, A. B. Antunes, C. S. Yadav, P. L. Paulose, B. Fauqu´ e, and K. Behnia, Strong correlation and low carrier density in Fe 1+ y Te 0 6 . Se 0 4 . as seen from its thermoelectric response, Phys. Rev. B 83 , 020504 (2011).
- [29] G. E. Volovik, Superconductivity with lines of GAP nodes: density of states in the vortex, Soviet Journal of Experimental and Theoretical Physics Letters 58 , 469 (1993).
- [30] N. Nakai, P. Miranovi´ c, M. Ichioka, and K. Machida, Field dependence of the zero-energy density of states around vortices in an anisotropic-gap superconductor, Phys. Rev. B 70 , 100503 (2004).
- [31] Y. Bang, Volovik effect in the ± s -wave state for the ironbased superconductors, Phys. Rev. Lett. 104 , 217001 (2010).
- [32] M. Zehetmayer, A review of two-band superconductivity: materials and effects on the thermodynamic and reversible mixed-state properties, Superconductor Science and Technology 26 , 043001 (2013).
- [33] P. Miranovi´ c, N. Nakai, M. Ichioka, and K. Machida, Orientational field dependence of low-lying excitations in the mixed state of unconventional superconductors, Phys. Rev. B 68 , 052501 (2003).
- [34] I. Vekhter, P. J. Hirschfeld, J. P. Carbotte, and E. J. Nicol, Anisotropic thermodynamics of d -wave superconductors in the vortex state, Phys. Rev. B 59 , R9023 (1999).
- [35] J. Xing, H. Lin, Y. Li, S. Li, X. Zhu, H. Yang, and H.H. Wen, Nodal superconducting gap in tetragonal FeS, Phys. Rev. B 93 , 104520 (2016).
- [36] P. K. Biswas, S. K. Ghosh, J. Zhao, D. A. Mayoh, N. Zhigadlo, X. Xu, C. Baines, A. Hillier, G. Balakrishnan, and
- M. Lees, Chiral singlet superconductivity in the weakly correlated metal LaPt 3 P, Nat. Commun. 12 , 2504 (2021).
- [37] G. Goll, Unconventional superconductors: experimental investigation of the order-parameter symmetry , Vol. 214 (Springer Science &amp; Business Media, 2006).
- [38] D. Pines and P. Monthoux, d x 2 -y 2 pairing and spin fluctuations in the cuprate superconductors: A progress report, J. Phys. Chem. Solids. 56 , 1651 (1995).
- [39] S. K. Ghosh, M. Smidman, T. Shang, J. F. Annett, A. D. Hillier, J. Quintanilla, and H. Yuan, Recent progress on superconductors with time-reversal symmetry breaking, Journal of Physics: Condensed Matter 33 , 033001 (2020).
- [40] L. Ding, C. He, J. K. Dong, T. Wu, R. H. Liu, X. H. Chen, and S. Y. Li, Specific heat of the iron-based highT c superconductor SmO 1 -x F FeAs, Phys. Rev. B x 77 , 180510 (2008).
- [41] M. Tropeano, A. Martinelli, A. Palenzona, E. Bellingeri, E. Galleani d'Agliano, T. D. Nguyen, M. Affronte, and M. Putti, Thermal properties of SmFeAsO 1 -x F x as a probe of the interplay between electrons and phonons, Phys. Rev. B 78 , 094518 (2008).
- [42] W. Schnelle, N. Knauf, J. Bock, E. Preisler, and J. H¨deu pohl, The specific heat anomaly of Bi 2 Sr CaCu O 2 2 8+ δ at Tc in dependence on the oxygen content, Physica C: Superconductivity 209 , 456 (1993).
- [43] C. Zhao, X. Liu, J. Wang, C. Xu, B. Wang, W. Xia, Z. Yu, X. Jin, X. Zhang, J. Wang, D. Dai, C. Tu, J. Nie, H. Wang, Y. Jiao, D. Duong, S. Huang, R. Jin, Z. Xu, Y. Guo, X. Xu, and S. Li, Multigap nodeless superconductivity in the topological semimetal pdte, Phys. Rev. B 109 , 174502 (2024).
- [44] W. J. Gannon, W. P. Halperin, C. Rastovski, K. J. Schlesinger, J. Hlevyack, M. R. Eskildsen, A. B. Vorontsov, J. Gavilano, U. Gasser, and G. Nagy, Nodal gap structure and order parameter symmetry of the unconventional superconductor UPt , 3 New J. Phys. 17 , 023041 (2015).
- [45] A. Bhattacharyya, D. T. Adroja, K. Panda, S. Saha, T. Das, A. J. S. Machado, O. V. Cigarroa, T. W. Grant, Z. Fisk, A. D. Hillier, and P. Manfrinetti, Evidence of a nodal line in the superconducting gap symmetry of noncentrosymmetric ThCoC , Phys. Rev. Lett. 2 122 , 147001 (2019).
- [46] T. Shang, M. Smidman, A. Wang, L.-J. Chang, C. Baines, M. K. Lee, Z. Y. Nie, G. M. Pang, W. Xie, W. B. Jiang, M. Shi, M. Medarde, T. Shiroka, and H. Q. Yuan, Simultaneous nodal superconductivity and timereversal symmetry breaking in the noncentrosymmetric superconductor CaPtAs, Phys. Rev. Lett. 124 , 207001 (2020).
- [47] D. F. Smith and C. P. Slichter, The study of mechanisms of superconductivity by nmr relaxation, in Novel NMR and EPR techniques , edited by J. Dolinˇ sek, M. Vilfan, and S. Zumer (Springer Berlin Heidelberg, Berlin, Heiˇ delberg, 2006) pp. 243-295.
- [48] R. Giannetta, A. Carrington, and R. Prozorov, London penetration depth measurements using tunnel diode resonators, J. Low Temp. Phys. 208 , 119 (2022).
C. S. Yadav Sudeep Kumar Ghosh Pankaj Kumar A. Thamizhavel P. L. Paulose
Signature of point nodal superconductivity in the Dirac semimetal PdTe
Recent Angle-Resolved Photo-emission Spectroscopy (ARPES) experiments [Phys.
Rev. Lett. 130, 046402 (2023)] on PdTe, a 3D-Dirac semimetal and a
superconductor with the transition temperature Tc ~ 4.3 K, have revealed
compelling evidence of the presence of bulk nodes in the superconducting order
parameter. To investigate the validity of this proposition, here we present a
detailed investigation of the magnetic field dependence of the specific heat of
PdTe down to temperatures ~ 58 mK. We observed that the low temperature
specific heat of PdTe with an externally applied magnetic field exhibits a
power-law field dependence, a characteristic of unconventional
superconductivity. Furthermore, the zero-field low-temperature electronic
specific heat follows a cubic temperature dependence, which is a signature of
the presence of bulk point nodes in PdTe. These intriguing observations suggest
that PdTe is a rare and fascinating topological material that exhibits both
Dirac semimetallic properties and superconductivity with point nodal gap
symmetry.