Front. Phys. Frontiers in Physics Front. Phys. 2296-424X Frontiers Media S.A. 651163 10.3389/fphy.2021.651163 Physics Mini Review Time-Reversal Symmetry Breaking in Re-Based Superconductors: Recent Developments Shang and Shiroka TRS Breaking in Re-Based Superconductor Shang Tian 1 2 * Shiroka Toni 3 4 * Key Laboratory of Polar Materials and Devices (MOE), School of Physics and Electronic Science, East China Normal University, Shanghai, China Laboratory for Multiscale Materials Experiments, Paul Scherrer Institut, Villigen, Switzerland Laboratory for Muon-Spin Spectroscopy, Paul Scherrer Institut, Villigen, Switzerland Laboratorium für Festkörperphysik, ETH Zürich, Zürich, Switzerland

Edited by: Yuji Muro, Toyama Prefectural University, Japan

Reviewed by: Amitava Bhattacharyya, Ramakrishna Mission Vivekananda Educational and Research Institute, India

Jess H. Brewer, University of British Columbia, Canada

 *Correspondence: Tian Shang, tshang@phy.ecnu.edu.cn; Toni Shiroka, tshiroka@phys.ethz.ch

This article was submitted to Condensed Matter Physics, a section of the journal Frontiers in Physics

 24 05 2021 2021 9 651163 08 01 2021 29 04 2021 Copyright © 2021 Shang and Shiroka. 2021 Shang and Shiroka

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

In the recent search for unconventional- and topological superconductivity, noncentrosymmetric superconductors (NCSCs) rank among the most promising candidate materials. Surprisingly, some of them—especially those containing rhenium—seem to exhibit also time-reversal symmetry (TRS) breaking in their superconducting state, while TRS is preserved in many other isostructural NCSCs. To date, a satisfactory explanation for such discrepant behavior, albeit crucial for understanding the unconventional superconductivity of these materials, is still missing. Here we review the most recent developments regarding the Re-based class, where the muon-spin relaxation (μSR) technique plays a key role due to its high sensitivity to the weak internal fields associated with the TRS breaking phenomenon. We discuss different cases of Re-containing superconductors, comprising both centrosymmetric- and noncentrosymmetric crystal structures, ranging from pure rhenium, to ReT (T = 3d-5d early transition metals), to the dilute-Re case of ReBe22. μSR results suggest that the rhenium presence and its amount are two key factors for the appearance and the extent of TRS breaking in Re-based superconductors. Besides summarizing the existing findings, we also put forward future research ideas regarding the exciting field of materials showing TRS breaking.

 time-reversal symmetry breaking noncentrosymmetric superconductors unconventional superconductivity muon-spin spectroscopy rhenium compounds Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung10.13039/501100001711 1 Introduction

The combination of intriguing fundamental physics with far-reaching potential applications has made unconventional superconductors one of the most studied classes of materials. Standing out among them are the noncentrosymmetric superconductors (NCSCs) [1], whose crystal structures lack the inversion symmetry. As a consequence, in NCSCs, the strict symmetry-imposed requirements are relaxed, allowing mixtures of spin-singlet and spin-triplet Copper pairing channels, thus setting the scene for a variety of exotic properties, as e.g., upper critical fields beyond the Pauli limit, nodes in the superconducting gaps, etc. (see Refs. [1, 2, 3] for an overview). The degree of mixing in such combined pairings is related to the strength of the antisymmetric spin-orbit coupling (ASOC) and to other microscopic parameters, still under investigation. Currently, NCSCs rank among the foremost categories of superconducting materials in which to look for topological superconductivity (SC) or to realize the Majorana fermions, pairs of the latter potentially acting as noise-resilient qubits in quantum computing [4–11].

In general, the various types of NCSCs can be classified into two classes. One consists of strongly correlated materials, as e.g., CePt3Si [12], or Ce(Rh,Ir)Si3 [13], which belong to the heavy-fermion compounds. Owing to the strong correlation and the interplay between d- and f-electrons, these materials often exhibit rich magnetic and superconducting properties. Since their superconductivity is most likely mediated by spin fluctuations, this implies an unconventional (i.e., non phonon-related) pairing mechanism. Conversely, the other class consists mainly of weakly correlated materials, which are free of “magnetic” f-electrons, as e.g., LaNiC2, La7Ir3, CaPtAs, or ReT (T = 3d-5d early transition metals) [14–20]. Obviously, their superconductivity is not mediated by the electrons’ spin fluctuations. Hence, they lead themselves as prototype parent systems where one can study the intrinsic pairing mechanisms in NCSCs.

Recently, superconductivity with broken time-reversal symmetry (TRS) has become a hot topic in NCSCs. The main reason for this is the discovery of TRS breaking in some weakly-correlated NCSCs using muon-spin relaxation (μSR). Surprisingly, the superconducting properties of the latter largely resemble those of conventional superconductors, i.e., their properties are clearly distinct from those of the above mentioned strongly-correlated NCSCs. To date, only a handful of NCSC families have been shown to exhibit TRS breaking in the superconducting state, including LaNiC2 [14], La7(Rh,Ir)3 [15, 21], Zr3Ir [22], CaPtAs [16], and ReT [14–20]. Except for the recently studied CaPtAs, where coexisting TRS breaking and superconducting gap nodes were observed below T c , in most of the above cases the superconducting properties evidence a conventional s-wave pairing, characterized by a fully opened superconducting gap. This leads to an interesting fundamental question: does the observed TRS breaking in NCSCs originate from an unconventional superconducting mechanism (i.e., from a pairing other than that mediated by phonons), or it can occur also in presence of conventional pairing (via some not yet understood mechanism) [2, 3]? Why, among the many different NCSCs families, only a few exhibit a broken TRS in the superconducting state, also remains an intriguing open question.

In general, the causes behind the TRS breaking at the onset of superconductivity are mostly unknown. In particular, the α-Mn-type noncentrosymmetric ReT (T = Ti, Nb, Zr, and Hf) superconductors have been widely studied and demonstrated to show a superconducting state with broken TRS [17–20]. Yet, TRS seems to be preserved in the isostructural (but Re-free) Mg10Ir19B16 and Nb0.5Os0.5 [23, 24]. Further, depending on the synthesis protocol, Re3W is either a centro- (hcp-Mg-type) or a noncentrosymmetric (α-Mn-type) superconductor, yet neither is found to break TRS [25]. In case of binary Re-Mo alloys, depending on the Re/Mo ratio, the compounds can exhibit up to four different crystal structures, including both centrosymmetric and noncentrosymmetric cases. Most importantly, all these alloys become superconductors at low temperatures [26]. A comparative μSR study of Re-Mo alloys, covering all the different crystal structures, reveals that the spontaneous magnetic fields occurring below T c (an indication of TRS breaking) were only observed in elementary rhenium and in Re0.88Mo0.12. By contrast, TRS was preserved in the Re-Mo alloys with a lower Re-content (below ∼ 88 % ), independent of their centro- or noncentrosymmetric crystal structures [27]. Since both pure rhenium and Re0.88Mo0.12 have a simple (hcp-Mg-type) centrosymmetric structure, this strongly suggests that a noncentrosymmetric structure and the accompanying ASOC effects are not essential in realizing the TRS breaking in ReT superconductors. The μSR results regarding the Re-Mo family, as well as other Re-free α-Mn-type superconductors, clearly imply that not only the Re presence, but also its amount are crucial for the appearance and the extent of TRS breaking in the ReT superconductors. How these results can be understood within a more general framework requires further experimental and theoretical investigation.

This short review article focuses mostly on the experimental study of Re-based binary superconductors. In Section 2, we discuss the basic principles of our probe of choice, the μSR, here used to detect the TRS breaking and to characterize the superconducting properties. Section 3 describes the possible crystal structures and superconducting transition temperatures of ReT binary alloys. In Section 4, we focus on the upper critical fields and the order parameter in ReT superconductors. Section 5 discusses the TRS breaking in ReT superconductors and its possible origins. Finally, in the last section, we outline some possible future research directions.

 2 Muon-Spin Relaxation and Rotation

Initially considered as an “exotic” technique, over the years muon-spin rotation, relaxation, and resonance (known as μSR), has become one of the most powerful methods to study the magnetic and superconducting properties of matter. This follows from a series of fortunate circumstances, related to the muon properties as a fundamental particle. Most notably, these include the 100% initial muon-spin polarization, following the two-body decay from pions, and the subsequent preservation of such information through the weak decay into positrons. In the search for unconventional superconductivity, as well as for TRS breaking effects, the very high sensitivity of the μSR technique to tiny magnetic fields is especially important [28]. Below we briefly outline the basics of the μSR technique and direct the reader to other references for more detailed information [29–31].

2.1 Principles of the μSR Technique

Central to the μSR method is the availability of polarized positive muon ( μ + ) beams, obtained by collecting the muons produced in the two-body decay of positive pions, π + → μ + + ν μ (with ν μ the muon neutrino), decaying at rest in the laboratory frame. Since pions have no intrinsic angular momentum and neutrinos have a fixed helicity (relative orientation of spin and linear momentum), the resulting muon beam is 100% spin polarized, with the muon spins directed antiparallel to the linear momentum (see Figure 1A). Having an energy of ∼ 4 . 12   MeV , muons can penetrate a sample between 0.1 and 1 mm, depending on the sample density. Once implantated, the monoenergetic muons decelerate within 100 ps through ionization processes (which do not perturb the muon spin) and finally come to rest at an interstitial site, practically without loss of their initial spin polarization. From this moment on, if subject to magnetic interactions, the muon-spin polarization P ( t ) evolves with time (the muon spin precesses around the local magnetic field), thus providing important information on the sample’s magnetism. The detection of the P ( t ) evolution is made possible by the parity-violating weak-decay interaction μ + → e + + ν e + ν ¯ μ ( e + , ν e , and ν ¯ μ are the positron, electron neutrino, and muon antineutrino, respectively), which implies a preferential emission of positrons along the muon-spin direction at the time of decay (see Figure 1A, which depicts also the anisotropic positron-emission pattern). Thus, by detecting the spatial distribution of positrons as a function of time, one can determine the time evolution of the muon-spin polarization P ( t ) .

Principle of the time-differential μSR experiment. (A) An incoming polarized muon (with spin S μ antiparallel to momentum p μ ) is first detected by a (thin) muon detector, which starts the electronic clock. In the sample, the muon spin precesses in the internal/external field until the muon decays into a positron. This is emitted preferentially along the muon-spin direction and hits one of the positron detectors [here forward (F) or backward (B)], whose signal stops the clock. The gray curve depicts the anisotropic positron-emission pattern at the moment of muon implantation. This pattern rotates rigidly with the muon spin (initially pointing toward the B detector) at an angular frequency ω μ = γ μ B ext . (B) Detected positron counts in the F and B detectors as a function of time after ca. 10 7 events. Inset: The asymmetry signal, obtained as the normalized difference between the F and B counts, eliminates the inessential exponential decay and highlights the signal decay (here assumed to be Gaussian) reflecting the nature of the sample.

A schematic diagram of a time-differential μSR experiment is shown in Figure 1A. The incoming muon triggers a clock that defines the starting time t 0 . Once implanted, the muon spin precesses about the local magnetic field B ( r ) with a Larmor frequency ω μ = γ μ B ( r ) , where γ μ / 2 π = 135.53 MHz/T is the muon gyromagnetic ratio. The clock stops when, after a mean lifetime of 2.197 μs, the muon decays into a positron e + , registered as an event by one of the positron detectors. The measured time intervals for ca. 10–50 millions of such events are stored in a histogram, given by (see Figure 1B): N ( t ) = N 0 exp ( − t / τ μ ) [ 1 + A 0 P ( t ) ] + C . Here, the exponential factor accounts for the radioactive muon decay, N 0 is the initial count rate at time t 0 , while C is a time-independent background (due to uncorrelated start and stop counts). As shown in the inset of Figure 1B, unlike the inessential exponential decay, the physical information in a μSR experiment is contained in the A ( t ) = A 0 P ( t ) term (often known as the μSR spectrum). Here, A 0 is the so-called initial asymmetry (typically 0.3, depending on the detector’s solid angle and efficiency), while P ( t ) is the muon-spin depolarization function, here given by the projection of P ( t ) on the unit vector describing the detector. Since P ( t ) represents the autocorrelation function of the muon spin S , i.e., P ( t ) = 〈 S ( t ) S ( 0 ) 〉 / S ( 0 ) 2 , it depends on the average value, the distribution, and the time evolution of the internal magnetic fields, thus reflecting the physics of the magnetic interactions in the sample under study. To access the μSR signal we need to remove the extrinsic decay factor by combining the positron counts from pairs of opposite-lying detectors, for instance, N F and N B (for forward and backward), and obtain the asymmetry A ( t ) = [ N F ( t ) − α N B ( t ) ] / [ N F ( t ) + α N B ( t ) ] . Clearly, A ( t ) behaves as a normalized “contrast”, proportional to A 0 . The parameter α is introduced to take into account the different efficiencies of the positron detectors and has to be determined by calibration.

 2.2 Transverse-Field μSR

Depending on the reciprocal orientation of the external magnetic field B with respect to the initial muon-spin direction S ( 0 ) , in a μSR experiment, two different configurations are possible. i) In transverse-field (TF) μSR the externally applied field B is perpendicular to S ( 0 ) and the muon spin precesses around B (see Figure 1A). ii) In a longitudinal field (LF) configuration the applied field is parallel to S ( 0 ) , generally implying only a relaxing μSR signal.

Although, in principle, the TF scheme shown in Figure 1A works fine, strong transverse fields perpendicular to the muon momentum ( p μ ) would deviate the muon beam too much from its original path. The resulting Lorentz force can be zeroed by applying the field along the muon momentum. At the same time, to maintain the transverse geometry, the initial muon spin is rotated by 90 ∘ (in the x or y direction) by using a so-called Wien filter, consisting of crossed electric and magnetic fields. Such a configuration is also known as transverse muon-spin mode, while Figure 1A plots the longitudinal muon-spin mode (i.e., p μ ∥ S μ ).

Since muons are uniformly implanted in the sample, they can detect the coexistence of different domains, characterized by distinct P i ( t ) functions, whose amplitudes A i represent a measure of the respective volume fractions. In case of superconductors, one can thus extract the SC volume fraction. More importantly, in a TF-μSR experiment one can directly probe the SC flux-line lattice (FLL). In this case, at the onset of superconductivity, the muon-spin precession in a TF field loses coherence, reflecting the magnetic field modulation (i.e., broadening) due to the FLL. The shape of the field distribution arising from the FLL can be analyzed and eventually used to extract the magnetic penetration depth λ and the coherence length ξ [32]. In many type-II superconductors, the simple relation, σ sc 2 / γ μ 2 = 0.00371 Φ 0 2 / λ eff 4 , connects the muon-spin depolarization rate in the superconducting phase, σ sc , with the effective magnetic penetration depth, λ eff (here Φ 0 is the magnetic flux quantum) [33, 34]. In case of superconductors with relatively low upper critical fields, the effects of the overlapping vortex cores with increasing field ought to be considered when extracting the magnetic penetration depth λ eff from σ sc . Since λ ( T ) is sensitive to the low-energy excitations, its evolution with temperature is intimately related to the structure of the superconducting gap. Hence, μSR allows us to directly study the symmetry and value of the superconducting gap.

More in detail, in a TF-μSR experiment, the time evolution of the asymmetry can be modeled by: A TF ( t ) = ∑ i = 1 n A i cos ( γ μ B i t + ϕ ) e − σ i 2 t 2 / 2 + A bg cos ( γ μ B bg t + ϕ ) . Here A i , A bg and B i , B bg are the asymmetries and local fields sensed by the implanted muons in the sample and the sample holder, γ μ is the muon gyromagnetic ratio, ϕ is a shared initial phase, and σ i is the Gaussian relaxation rate of the ith component. The number of required components is material dependent, typically in the 1 ≤ n ≤ 5 range. In general, for superconductors with a large Ginzburg-Landau parameter κ ( ≫ 1), the magnetic penetration depth is much larger than the coherence length. Hence, the field profiles of each fluxon overlap strongly, implying a narrow field distribution. Consequently, a single-oscillating component is sufficient to describe A ( t ) . In case of a small κ ( ≥ 1 / 2 ), the magnetic penetration depth is comparable with the coherence length. Here, the small penetration depth implies fast-decaying fluxon field profiles and a broad field distribution, in turn requiring multiple oscillations to describe A ( t ) . The choice of n can be determined from the fast-Fourier-transform (FFT) spectra of the TF-μSR, which is normally used to evaluate the goodness of the fits. In case of multi-component oscillations, the first term in Eq. 2 describes the field distribution as the sum of n Gaussian relaxations [35]: p ( B ) = γ μ ∑ i = 1 n A i σ i   exp [ − γ μ 2 ( B − B i ) 2 2 σ i 2 ] . The first- and second moments of the field distribution in the sample can be calculated by: 〈 B 〉 = ∑ i = 1 n A i B i A tot   and   〈 B 2 〉 = σ eff 2 γ μ 2 = ∑ i = 1 n A i A tot [ σ i 2 γ μ 2 + ( B i − 〈 B 〉 ) 2 ] , where A tot = ∑ i = 1 n A i . The total Gaussian relaxation rate σ eff in Eq. 4 includes contributions from both a temperature-independent relaxation, due to nuclear moments ( σ n ), and a temperature-dependent relaxation, related to the FLL in the superconducting state ( σ sc ). The σ sc values are then extracted by subtracting the nuclear contribution following σ sc = σ e f f 2 − σ n 2 .

To get further insights into the superconducting gap value and its symmetry, the temperature-dependent superfluid density ρ sc ( T ) [proportional to λ eff − 2 ( T ) ] is often analyzed by using a general model: ρ sc ( T ) = λ 0 2 λ eff 2 ( T ) = 1 + 2   〈 ∫ Δ k ∞ E E 2 − Δ k 2   ∂ f ∂ E 〉 FS Here, λ 0 is the effective magnetic penetration depth in the 0-K limit, f = ( 1 + e E / k B T ) − 1 is the Fermi function and 〈 〉 FS represents an average over the Fermi surface [36]. Δ k ( T ) = Δ ( T ) A k is an angle-dependent gap function, where Δ is the maximum gap value and A k is the angular dependence of the gap, equal to 1, cos 2 ϕ , and sin θ for an s-, d-, and p-wave model, respectively, with ϕ and θ being the azimuthal angles. The temperature dependence of the gap is assumed to follow the relation Δ ( T ) = Δ 0 tanh { 1.82 [ 1.018 ( T c / T − 1 ) ] 0.51 } [36, 37], where Δ 0 is the 0-K gap value.

 2.3 Zero-Field μSR

A particular case of LF, is that of zero-field (ZF) μSR, characterized by the absence of an external magnetic field. In this configuration the frequency of the μSR signal is exclusively proportional to the internal magnetic field, making it possible to determined the size of the ordered moments and, hence, the magnetic order parameter. Unlike various techniques, which require an external field to polarize the probe, μSR is unique in its capability of studying materials unperturbed by externally applied fields and in accessing their spontaneous magnetic fields. Due to the large muon magnetic moment ( μ μ = 8.89   μ N ), ZF-μSR can sense even very small internal fields ( ∼ 10 − 2   mT ), and, thus, can probe local magnetic fields of either nuclear or electronic nature. In addition, since the muon is an elementary spin-1/2 particle, it acts as a purely magnetic probe, i.e., free of quadrupole interactions. All these features make ZF-μSR an ideal technique for detecting TRS breaking in the superconducting state. The latter corresponds to the appearance (at the onset of SC) of spontaneous magnetic moments, whose magnitude can be very small, often lacking a proper magnetic order. As we show further on, in case of TRS breaking, we expect the appearance of an additional enhancement of μSR relaxation below T c , reflecting the occurrence of such weak spontaneous fields. During the ZF-μSR measurements, to exclude the possibility of stray magnetic fields (typically larger than the weak internal fields), the magnets are quenched before starting the measurements, and an active field-nulling facility is used to compensate for stray fields down to 1 μT.

If the amplitudes of the local fields reflect a Gaussian distribution with zero average (a rather common circumstance), the μSR signal consists of overlapping oscillations with different frequencies. While at short times the spin dephasing is limited, at long times it becomes relevant and gives rise to a so-called Kubo-Toyabe (KT) relaxation function [31, 38]. Two different models are frequently used to analyze the ZF-μSR data: A ZF = A s [ 1 3 + 2 3 ( 1 − σ ZF 2 t 2 − Λ ZF t )   e ( − σ ZF 2 t 2 2 − Λ ZF t ) ] + A bg , or A ZF = A s [ 1 3 + 2 3 ( 1 − σ ZF 2 t 2 )   e − σ ZF 2 t 2 2 ] e − Λ ZF t + A bg . Equation 6 is also known as a combined Gaussian- and Lorentzian Kubo-Toyabe function, with the additional exponential relaxation describing the electronic contributions present in many real materials. In polycrystalline samples, the 1/3-non-relaxing and the 2/3-relaxing components of the asymmetry correspond to the powder average of the internal fields with respect to the initial muon-spin direction (statistically, with a 1/3 probability, the directions of the muon spin and of the local field coincide). Clearly, in the case of single crystals, such prefactors might be different. The σ ZF and Λ ZF represent the zero-field Gaussian and Lorentzian relaxation rates, respectively. Typically, Λ ZF shows an almost temperature-independent behavior. Hence, an increase of σ ZF across T c can be attributed to the spontaneous magnetic fields which break the TRS, as e.g., in ReT [17, 18, 20]. In case of diluted nuclear moments, σ ZF is practically zero, hence, the TRS breaking is reflected in an increase of Λ ZF below T c , as e.g., in Zr3Ir and CaPtAs [16, 22].

 3 Re-Based Superconductors

In this section, we review the different phases of the binary ReT alloys. These are obtained when rhenium reacts with various early transition metals (see Figure 2A) and show rich crystal structures. Representative examples are shown in Figures 2C–F, including the hexagonal hcp-Mg- ( P 6 3 / m m c , No. 194), cubic α-Mn- ( I 4 ¯ 3 m , No. 217), tetragonal β-CrFe- ( P 4 2 / m n m , No. 136), and cubic bcc-W-type ( I m 3 ¯ m , No. 229). Among these the cubic α-Mn-type structure is noncentrosymmetric, while the rest are centrosymmetric [39]. Besides the above cases, a few other crystal structures have also been reported, including the cubic CsCl- (Pm-3m, No. 221), cubic Cr3Si- (Pm-3n, No. 223), and trigonal Mn21Zn25-type (R-3c, No. 167) [39]. As for the pure elements listed in Figure 2A, both Re and Os have an hcp-Mg-type structure, and show superconductivity below 2.7 and 0.7 K, respectively [20, 40]; while V, Nb, Mo, Ta, and W all adopt a bcc-W-type structure, and become superconductors at ∼5.4, 9.0, 1.0, 4.5, and 0.015 K, respectively [40]. Unlike the above cases, Ti, Zr, and Hf can form either high-temperature bcc-W-type or low-temperature hcp-Mg-type structures, with T c ∼ 0.4, 0.6, and 0.13 K, respectively [40].

Crystal structures of rhenium transition-metal (ReT) superconductors. (A) List of 3d, 4d, and 5d early transition metals, which can react with rhenium to form different crystal structures. (B) Binary phase diagram for the typical case of Re-Mo alloys (data adopted from Ref. [39]). (C–F) Unit cells of four most representative crystal structures of ReT binary compounds. Among these the cubic α-Mn type ( I 4 ¯ 3 m , No. 217) in (D) is noncentrosymmetric, while the hexagonal hcp-Mg ( P 6 3 / m m c , No. 194), tetragonal β-CrFe ( P 4 2 / m n m , No. 136), and cubic bcc-W ( I m 3 ¯ m , No. 229) are centrosymmetric. The atomic coordinates for each structure can be found in Refs. [26, 93].

For T = Ti 3d metal, the known binary compounds are Re24Ti5, Re6Ti, and ReTi [41, 42]. The former two adopt a noncentroysmmetric α-Mn-type structure and become superconductors below T c = 6 K [19, 40, 43], while the latter one crystallizes in a cubic CsCl-type structure. To date, no detailed physical properties have been reported for ReTi. For T = V, superconductivity has been reported in hcp-Mg-type Re0.9V0.1 ( T c = 9.4 K ), β-CrFe-type Re0.76V0.24 ( T c = 4.5 K ), and bcc-W-type Re0.6V0.4 ( T c = 2.2 K ) [40, 44]. Also for them, to date a microscopic study of their SC is still missing. The cubic Cr3Si-type Re0.71V0.29 has also been synthesized, but its physical properties were never characterized [45].

For T = Zr 4d metal, the α-Mn-type Re24Zr5 ( T c = 5 K ) and Re6Zr ( T c = 6.7 K ) have been investigated via both nuclear quadrupole resonance and μSR techniques [17, 46]. Except for the α-Mn-type Re-Zr alloys, the MgZn2-type Re2Zr (similar to hcp-Mg-type) and Mn21Zn25-type Re25Zr21 have been synthesized [39]. Re2Zr exhibits a T c value of ∼6–7 K [40, 47], while Re25Zr21 has not been studied. For T = Nb, depending on Re/Nb concentration, four different solid phases including hcp-Mg-, α-Mn-, β-CrFe-, and bcc-W-type have been reported. On the Re-rich side, the hcp-Mg-type Re-Nb alloys are limited to less than 3% Nb concentration [39], whereas many α-Mn-type Re-Nb binary alloys have been grown and widely studied by various techniques [20, 48–51], with the highest T c reaching 8.8 K in Re24Nb5 (denoted as Re0.82Nb0.18 in the original paper [20]). At intermediate Re/Nb values, for example, in β-CrFe-type Re0.55Nb0.45, T c s in the range of 2–4 K [40] have been reported, but no microscopic studies yet. As for the Nb-rich side (Nb concentration larger than 60%), here the Re-Nb alloys exhibit the same structure as that of pure Nb, but much lower T c values than Nb [39, 40]. For T = Mo, the binary Re-Mo phase diagram (see Figure 2B) covers also four different solid phases [39]. The binary Re-Mo alloys have been characterized by different techniques and all of them become superconductors at low temperatures [26, 27]. The T c varies nonmonotonically upon changing the Mo concentration, giving rise to three distinct superconducting regions. On the Re-rich side, the first SC region shows the highest T c ∼ 9.4   K in the hcp-Mg-type Re0.77Mo0.23. The same material but with an α-Mn-type structure can also be grown, with a T c value about 1 K lower than the hcp-Mg-type. In the second superconducting region, where the alloys adopt a β-CrFe-type structure, the superconducting transition temperature T c ∼ 6.3   K is almost independent of Mo content. Finally, on the Mo-rich side, all Re-Mo alloys display a cubic bcc-W-type structure and form a third superconducting region with the highest T c reaching 12.4 K in Re0.4Mo0.6.

For T = Hf 5d metal, the Re-Hf alloys show a similar phase diagram to Re-Zr. With only ∼3% Hf substitution, T c increases from <3 to 7.3 K in the hcp-Mg-type Re-Hf alloys [40]. Both the α-Mn-type Re6Hf and the MgZn2-type Re2Hf become superconductors below T c ∼ 6 K [18, 40, 47, 52, 53], whereas the physical properties of Mn21Zn25-type Re25Hf21 remain largely unknown. On the Hf-rich side, the bcc-W-type alloys exhibit relatively low T c s , e.g., T c = 1.7   K for Hf0.875Re0.125 [40]. For T = Ta, although the four different structures shown in Figure 2C–F can be synthesized, only the α-Mn-type Re-Ta alloys have been well studied. For example, Re3Ta and Re5.5Ta show T c values of 4.7 and 8 K, respectively [54, 55]. On the Ta-rich side, the bcc-W-type Re-Ta alloys become superconducting at T c < 3.5   K , lower than the T c of pure Ta [56]. We note that in case of the β-CrFe-type Re-Ta alloys, no superconducting transition has been observed down to 1.8 K in either Re0.5Ta0.5 or Re0.6Ta0.4. For T = W, the Re-W alloys show a very similar phase diagram to Re-Mo in Figure 2B. As the W concentration increases, the highest T c values reach ∼8, 9, 6, and 5 K in the hcp-Mg-, α-Mn-, β-CrFe-, and the bcc-W-type alloys, respectively [39, 40]. Among them, only the hcp-Mg- and the α-Mn-type Re3W have been investigated [25, 57]. Finally, in case of T = Os, the Re-Os alloys show a rather monotonous phase diagram, since only hcp-Mg-type compounds with T c values below 2 K can be synthesized [39, 40].

 4 Upper Critical Field and Nodeless Superconductivity

As mentioned in the introduction, due to the mixture of singlet- and triplet paring, some NCSCs may exhibit relatively high upper critical fields, often very close to or even exceeding the Pauli limit, as e.g., CePt3Si [12], Ce(Rh,Ir)Si3 [58, 59], and recently (Ta,Nb)Rh2B2 [60]. Therefore, the upper critical field can provide valuable clues about the nature of superconductivity. To investigate the temperature evolution of the upper critical field H c 2 ( T ) , in general, the temperature- (or field-) dependent electrical resistivity ρ, magnetic susceptibility χ, and specific heat C / T at various magnetic fields (or at various temperatures) are measured [19, 20, 27]. As an example, Figure 3A shows the H c 2 ( T ) for Re24Nb5 (α-Mn-type) and Re0.4Mo0.6 (bbc-W-type) versus the normalized temperature T / T c ( 0 ) . To obtain the upper critical field in the zero-temperature limit, H c 2 ( 0 ) , the Werthamer-Helfand-Hohenberg (WHH) or the Ginzburg-Landau (GL) models are usually applied when analyzing the H c 2 ( T ) data of ReT superconductors. Both models can adequately describe single-gap superconductors. Here, in case of Re24Nb5 and Re0.4Mo0.6, the WHH model (solid line in Figure 3A) reproduces the data very well and gives μ 0 H c 2 ( 0 ) =   15 . 6  T, and  3 . 08  T , respectively. Figure 3B summarizes the μ 0 H c 2 ( 0 ) values of the ReT and α-Mn-type NbOs2 superconductors. As discussed in Section 3, most of the previous studies focused exclusively on α-Mn-type ReT superconductors, the physical properties of the other ReT superconductors being practically neglected and requiring further studies. Unlike other ReT, all Re-Mo alloys belonging to four different structures have been studied via macro- and microscopic techniques [26, 27]. The μ 0 H c 2 ( 0 ) of centrosymmetric Re-Mo alloys, including hcp-Mg-, β-CrFe-, and bcc-W-type, are far away from the Pauli limit μ 0 H P = 1.86   T c (indicated by a dashed line in Figure 3B). Conversely, the α-Mn-type ReT and NbOs2 both exhibit large upper critical fields, very close to or even exceeding the Pauli limit, despite their different T c values. For example, μ 0 H c 2 ( 0 ) = 15.6  and  16 . 5   T for Re24Nb5 and Re5.5Ta, while their μ 0 H P ( 0 ) are 16.4 and 14.9 T, respectively. The hcp-Mg-type Re3W also exhibits a relatively high H c 2 , as determined from electrical resistivity data. However, its H c 2 value might be overestimated since, e.g., at 9 T, no zero resistivity could be observed down to 2 K. Therefore, other bulk techniques, including magnetization- or heat capacity measurements are required to determine the intrinsic H c 2 . In general, it would be interesting to know the H c 2 values of other centrosymmetric ReT superconductors. Overall, the upper critical fields in Figure 3B indicate the possibility of singlet-triplet mixing in the noncentrosymmetric α-Mn-type superconductors.

Upper critical field and superconducting energy gap. (A) The upper critical field H c 2 , as determined from electrical resistivity-, heat capacity-, and magnetic susceptibility measurements, as a function of the reduced superconducting transition temperature T/Tc (0) for Re24Nb5 and Re0.4Mo0.6. Solid-lines represent fits to the Werthamer-Helfand-Hohenberg (WHH) model. (B) Zero-temperature H c 2 versus the superconducting transition temperature T c for α-Mn type NbOs2 and all ReT superconductors. The shaded region in (B) marks the noncentrosymmetric α-Mn type superconductor, while the dashed line indicates the Paul limit (i.e., μ 0 H P = 1.86 T c ). (C) Superfluid density vs reduced temperature T/Tc (0) for Re24Nb5 and Re0.4Mo0.6. Lines are fits to a fully-gapped s-wave model. The insert shows the TF-μSR spectra for Re0.4Mo0.6 measured in a field of 60 mT in the normal- (16 K) and the superconducting state (1.5 K). Solid lines are fits to Eq. 2. (D) Zero-temperature superconducting energy gap Δ 0 (in k B T c units) as a function of T c for ReT and α-Mn-type NbOs2 and TaOs superconductors. Here, the dashed line represents the BCS superconducting gap in the weak-coupling limit (i.e., 1.76 k B T c ). Data were taken from Refs. [17, 18, 19, 20, 25, 27, 43, 52, 53, 55, 57, 61, 62, 63].

Transverse-field μSR represents one of the most powerful techniques to investigate the superconductivity at a microscopic level. To illustrate this, in the inset of Figure 3C we show two typical TF-μSR spectra for bcc-W-type Re0.4Mo0.6 in the normal and the superconducting states. Below T c , the fast decay induced by FLL (encoded into σ sc ) is clearly visible, while the slow decay in the normal state is attributed to the randomly oriented nuclear magnetic moments. By comparing the two spectra, one can also determine the superconducting volume fraction of a superconductor. As an example, the main panel of Figure 3C shows the normalized superfluid density calculated from σ sc ( T ) , which is proportional to [ λ ( T ) / λ ( 0 ) ] − 2 (see details in Section 2.2), as a function of the reduced temperature T / T c ( 0 ) for Re24Nb5 and Re0.4Mo0.6 [20, 27]. The low-T superfluid density is practically independent of temperature, clearly suggesting a lack of low-energy excitations and a fully-gapped superconductivity. Contrarily, such excitations exist in case of nodes in the superconducting gap, implying a temperature-dependent superfluid density below ∼ T c / 3 . As shown by solid lines in Figure 3C, the ρ sc ( T ) of ReT superconductors is described very well by a fully-gapped s-wave model (see Eq. 5). The other α-Mn-type ReT, TaOs, and NbOs2 exhibit similar temperature-invariant superfluid densities below T c /3 [17, 18, 19, 25, 43, 55, 61, 62]. Although ReT alloys adopt different crystal structures (i.e., centrosymmetric or noncentrosymmetric, see Figure 2C–F) and have different T c values, they regularly exhibit low-T superfluid densities which are independent of temperature [27]. Except for T = Mo (and for some α-Mn structures), a systematic microscopic study of superconductivity in other ReT superconductors is still missing. Clearly, it would be interesting to know if their SC behavior is similar to that of Re-Mo alloys. The nodeless SC scenario in ReT alloys is also supported by other techniques, as the electronic specific heat, the magnetic penetration depth measured via the tunnel-diode-oscillator-based technique, or the point-contact Andreev reflection [19, 20, 27, 52, 53, 57, 63, 64, 65]. In addition, some studies have found evidence of two-gap SC in Re0.82Nb0.18 and Re6Zr [48, 65].

Figure 3D summarizes the zero-temperature superconducting energy gap value for ReT and α-Mn-type NbOs2 and TaOs superconductors as a function of their critical temperatures. Most of them exhibit a Δ0/kB T c ratio larger than 1.76 (see dashed line in Figure 3D), the value expected for a weakly coupled BCS superconductor, which indicates a moderately strong coupling in these superconductors. In addition, the specific-heat discontinuity at T c (i.e., Δ C / γ T c ) is larger than the conventional BCS value of 1.43, again indicating an enhanced electron-phonon coupling [19, 20, 27, 52, 53, 57, 63]. As mentioned above, it is worth noting that the superconducting parameters of all the other ReT materials (except for α-Mn-type and T = Mo) are missing, prompting further research efforts in this direction.

As discussed in the introduction, the lack of inversion symmetry in the NCSCs often induces an ASOC. This splits the Fermi surface by lifting the degeneracy of the conduction electrons, thus allowing admixtures of spin-singlet and spin-triplet pairing. In general, the strength of ASOC determines the degree of such an admixture and thus the superconducting properties of NCSCs [1, 2]. A fully-gapped superconductor (i.e., dominated by spin-singlet pairing) can be tuned into a nodal superconductor (dominated by spin-triplet pairing) by increasing the strength of ASOC. Such mechanism has been successfully demonstrated, e.g., in weakly-correlated Li2Pt3B ( E SOC / k B T c ∼ 831 ) [66, 67], CaPtAs ( E SOC / k B T c ∼ 800 ) [16, 68], and in strongly-correlated CePt3Si ( E SOC / k B T c ∼ 3095 ) superconductors [12, 69], all exhibiting a relatively large band splitting E SOC compared to their superconducting energy scale k B T c . In the α-Mn-type ReT alloys, the density of states (DOS) near the Fermi level is dominated by the 5d orbitals of rhenium atoms, while contributions from the d orbitals of T atoms are negligible [70–72]. Therefore, a possible enhancement of SOC due to 3d-(e.g., Ti, V) up to 5d-electrons (e.g., Hf, Ta, W, Os) will, in principle, neither increase the band splitting E SOC nor affect the pairing admixture and thus the superconducting properties of α-Mn-type ReT superconductors. According to band-structure calculations, in Re6Zr, the SOC-induced band splitting is about 30 meV [72], implying a very small ratio E SOC / k B T c ∼ 25 , comparable to that of fully-gapped Li2Pd3B, Mo3P, and Zr3Ir superconductors [22, 67, 73]. Therefore, despite the relatively large SOC of rhenium atoms, its effects are too weak to significantly influence the bands near the Fermi level. This might explain why all the α-Mn-type ReT superconductors exhibit nodeless superconductivity, more consistent with a spin-singlet dominated pairing [17, 18, 19, 20, 25, 54]. However, we recall that often, due to the similar magnitude and same-sign of the order parameter on the spin-split Fermi surfaces, a possible mixed-pairing superconductor may be challenging to detect or to distinguish from a single-gap s-wave superconductor [74]. The almost spherical symmetry of the Fermi surface in these materials may also explain their BCS-like superconducting states [71]. As for the other centrosymmetric ReT alloys, in most of them the Re and T atoms occupy the same atomic positions in the unit cell. In this case, as the T-content increases, the contribution of T d orbitals to the DOS is progressively enhanced, at the expense of the Re 5 d orbitals. Therefore, the chemical substitution of Re by another 3d, 4d, or 5d T metal (see Figure 2), should significantly tune the SOC and, hence, the band splitting, an interesting hypothesis waiting for theoretical confirmation. However, even for T = Hf, Ta, W, and Os, the maximum E SOC should still be comparable to that of α-Mn-type ReT alloys. Finally, irrespective of the strength of SOC, due to their centrosymmetric crystal structures, these compounds may exhibit either singlet- or triplet-pairing, but not an admixture of both. According to the TF-μSR results (see Figure 3C), despite a change in SOC and of the different crystal structures (see Figure 2C–F), all ReT superconductors exhibit fully-gapped superconducting states. This finding strongly suggests that, in the ReT superconductors, spin-singlet pairing is dominant.

 5 Time-Reversal Symmetry Breaking

Owing to its very high sensitivity (see details in Section 2.3), ZF-μSR has been successfully used to search for spontaneous magnetic fields, reflecting the breaking of TRS in the superconducting states of different types of superconductors, as e.g., Sr2RuO4, UPt3, PrOs4Sb12, LaNiGa2, LaNiC2, La7(Rh,Ir)3, and α-Mn-ReT [14, 15, 17, 18, 19, 20, 21, 75, 76, 77, 78, 79, 80]. The latter three are typical examples of weakly-correlated NSCSs, to be contrasted with strongly-correlated NCSCs, where either the TRS is broken by a coexisting long-range magnetic order, or the tiny TRS-breaking signal is very difficult to detect due to the presence of strong magnetic fluctuations [28]. In the former case, the broken TRS is unrelated to the superconductivity, while in the later case, a genuine TRS breaking effect is masked by the much faster muon-spin relaxation caused by magnetic fluctuations. Therefore, in general, a TRS breaking effect is more easily (and reliably) detected in weakly-correlated- or non-magnetic superconductors using μSR techniques. Normally, in the absence of external fields, the onset of superconductivity does not imply changes in the ZF-μSR relaxation rate. However, in presence of a broken TRS, the onset of a tiny spontaneous polarization or of currents gives rise to associated (weak) magnetic fields, readily detected by ZF-μSR as an increase in the relaxation rate. Given the tiny size of such effects, the ZF-μSR measurements are usually performed in both the normal- and the superconducting state with a relatively high statistics, at least twice that of the TF-μSR spectra. As an example, Figure 4A plots the ZF-μSR spectra of α-Mn-type Re25Nb5, with the other ReT superconductors showing a similar behavior. The ZF-μSR spectra collected below- and above T c (at 1.5 and 12 K) exhibit small yet measurable differences. The lack of any oscillations in the spectra, implies the non-magnetic nature of ReT superconductors. Further, longitudinal-field μSR measurements under a relatively small applied field (typically a few tens of mT) in the superconducting state are usually performed to check if the applied field can fully decouple the muon spins from the weak spontaneous magnetic fields, and thus exclude extrinsic effects. In non-magnetic materials in the absence of external magnetic fields, the muon-spin relaxation is mostly determined by the interaction between the muon spins and the randomly oriented nuclear magnetic moments. Therefore, the spontaneous magnetic fields due to the TRS breaking will be reflected in an additional increase of muon-spin relaxation. The ZF-μSR asymmetry can be described by means of a Gaussian- or a Lorentzian Kubo-Toyabe relaxation, or a combination thereof (see Eqs. 6, 7). Figure 4B summarizes the Gaussian relaxation rate σ ZF versus the reduced temperature T / T c for the α-Mn-type Re24Nb5, Re6Zr, and Re0.77Mo0.23, and the hcp-Mg-type elementary Re. Above T c , all the samples show a temperature-independent σ ZF . Except for Re0.77Mo0.23, a small yet clear increase of σ ZF ( T ) below T c indicates the onset of spontaneous magnetic fields, which represent the signature of TRS breaking in the superconducting state [17, 20, 27]. The other α-Mn-type superconductors, e.g., Re6Ti, and Re6Hf [18, 43], show similar σ ZF ( T ) to Re24Nb5 and Re6Zr, and thus the breaking of TRS in the superconducting state. At the same time, in the isostructural Re3Ta, Re5.5Ta, and Re3W cases, there is no clear increase in σ ZF ( T ) upon crossing T c , implying a preserved TRS [25, 54, 55].

ZF-μSR and evidence for TRS breaking. (A) ZF-μSR spectra for Re24Nb5 collected in the superconducting and normal states. Top: additional μSR data collected at 1.5 K in a 15-mT longitudinal field, are also shown. The solid lines are fits using Eq. 6. (B) Gaussian relaxation rate Δ σ ZF vs T / T c for Re24Nb5, Re6Zr, Re, and Re0.77Mo0.23 — here Δ σ ZF ( T ) = σ ZF ( T ) − σ ZF ( T > T c ) . While for the first three there is a clear increase of Δ σ ZF across T c (hence a breaking of TRS), no changes occur in the last case (TRS is preserved). Data were taken from Refs. [17, 19, 20, 27].

Recently, the breaking of TRS and the presence of nodes in the SC gap, attributed to an admixture of singlet- and triplet paring, has been reported in the noncentrosymmetric CaPtAs superconductor [16]. In general, however, the breaking of TRS in the superconducting state and a lack of space-inversion symmetry in the crystal structure are independent events, not necessarily occurring together. For instance, the unconventional spin-triplet pairing is expected to break TRS below T c , as has been shown, e.g., in Sr2RuO4, UPt3, and UTe2 triplet superconductors [75–77, 79, 81–85]. An s + i d spin-singlet state was proposed to account for the TRS breaking in some iron-based high- T c superconductors [86], where a nodal gap is also expected. The frequent occurrence of TRS breaking in the fully-gapped (i.e., dominated by spin-singlet pairing) ReT superconductors (see Section 4) is, therefore, rather puzzling. A similarly surprising result is the report that elementary rhenium also exhibits signatures of TRS breaking in its superconducting state (see Figure 4B), with Δ σ ZF ( T ) being comparable to that of Re6Zr [20, 27]. Since elementary rhenium adopts a centrosymmetric hcp-Mg crystal structure (see Figure 2C), this indicates that a lack of inversion symmetry and the accompanying ASOC effects are not crucial factors for the occurrence of TRS breaking in ReT superconductors. Further on, a comparison of ZF-μSR measurements on Re-Mo alloys with different Re/Mo contents, covering almost all the crystal structures reported in Figure 2, shows that only Re and Re0.88Mo0.12 exhibit a broken TRS in the superconducting state, while those with a higher Mo-content (∼23–60%), including both the centrosymmetric- and noncentrosymmetric Re0.77Mo0.23, preserve the TRS. Considering the preserved TRS in Mg10Ir19B16, and Nb0.5Os0.5 [23, 24], all of which share the same α-Mn-type structure, this implies that TRS breaking in ReT superconductors is clearly not related to the noncentrosymmetric crystal structure or to a possible mixed pairing but, most likely, is due to the presence of rhenium and to its amount. Such conclusion is further reinforced by the preserved TRS in many Re-based superconductors, whose Re-content is below a certain threshold. Such cases include, e.g., Re3W, Re3Ta, Re-Mo (with Mo-content higher than 12%) [25, 27, 54], the recently reported Re-B superconductors [87], and the diluted ReBe22 superconductor [88]. Moreover, by comparing the ZF-μSR relaxation across various ReT superconductors, a clear positive correlation between Δ σ ZF (i.e., spontaneous fields) and the size of the nuclear magnetic moments μ n was identified [20]. For instance, among the ReT superconductors, Re24Nb5 shows the largest spontaneous fields below T c (see Figure 4B), a fact compatible with the large nuclear magnetic moment of niobium, practically twice that of rhenium (6.17 vs. 3.2 µ N). However, the correlation between μ n and Δ σ ZF alone cannot explain TRS breaking, considering that elementary Nb itself, despite having the highest μ n , does not break TRS. Clearly, the origin of such correlation is not yet understood and it requires further experimental and theoretical studies.

If SOC can be ignored, an alternative mechanism, which can account for the TRS breaking in ReT superconductors in presence of a fully-opened superconducting gap, is the internally-antisymmetric nonunitary triplet (INT) pairing. The INT pairing was originally proposed to explain the TRS breaking and nodeless SC in centrosymmetric LaNiGa2 [3, 80, 89] and noncentrosymmetric LaNiC2 [14, 90], both exhibiting a relatively weak SOC. In case of INT pairing, the superconducting pairing function is antisymmetric with respect to the orbital degree of freedom, while remaining symmetric in the spin- and crystal-momentum channels [14, 80, 89, 90]. Since in ReT superconductors, too, the SOC interaction is relatively weak (∼30 meV, see Section 4) [72] and since neither TRS breaking nor the nodeless SC are related to the symmetry of ReT crystal structures, the effect of SOC to the observed TRS breaking is insignificant. This could, therefore, explain why a lack of inversion symmetry (essential to SOC) is not a precondition for TRS breaking in ReT superconductors. Moreover, the occurrence of an INT state relies on the availability of a local-pairing mechanism driven by Hund’s rules, e.g., by Ni 3d-electrons in LaNiC2 and LaNiGa2 [3, 14, 80, 89, 90]. Such local-pairing mechanism may also occur in ReT superconductors, since rhenium too can be magnetic [91, 92]. This consideration is also in good agreement with the observation that TRS breaking depends on Re content, but not on a noncentrosymmetric crystal structure.

 6 Conclusion

In this short review we focused on recent experimental studies of ReT superconductors, where time-reversal symmetry breaking effects are often present and whose superconductivity can, therefore, be considered as unconventional. Due to its high sensitivity to the weak internal fields associated with TRS breaking, μSR represents one of the key techniques in the search for TRS-breaking effects in the superconducting state. Nonetheless, in certain cases, the amplitude of the spontaneous magnetic fields (the fingerprint of TRS breaking) may still be below the resolution of the μSR technique ( ∼ 10 − 2   mT ). Hence, the future use of other techniques, e.g., based on the optical Kerr effect [11], another very sensitive probe of spontaneous fields in unconventional superconductors, remains crucial. Due to their rich crystal structures, covering both centro- and noncentrosymmetric cases, and the pervasive presence of superconductivity at low temperatures, the nonmagnetic Re-based materials are the ideal choice for investigating the origin of TRS breaking. Here, we reviewed different cases of Re-containing superconductors, ranging from elementary rhenium, to ReT (T = 3d-5d early transition metals), to the dilute-Re case of ReBe22, all of which were investigated through both macroscopic and microscopic techniques. By a comparative study of ReT with different T metals mostly using the μSR technique, we could demonstrate the secondary role played by SOC and why the spin-singlet pairing is dominant in ReT superconductors. This, however, brings up the question of reconciling the occurrence of TRS breaking with a fully-gapped SC state (spin-singlet pairing). A possible solution to this apparent contradiction is offered by the so-called INT model, which requires an antisymmetric pairing function involving the orbital degree of freedom, making it insensitive to the presence (or lack) of inversion symmetry and SOC. Overall, the reported results suggest that the rhenium presence and its amount are two key factors for the appearance and the extent of TRS breaking in the Re-based superconductors. These key observations, albeit important, demand new experimental and theoretical investigations to further generalize them.

To date, as nearly all current studies have focused exclusively on α-Mn-type ReT superconductors (except for the Re-Mo series considered here), the superconducting properties of most other ReT alloys remain basically unexplored. Hence, the synthesis and characterization of non-α-Mn-type ReT alloys, including the study of their electrical, magnetic, and thermodynamic properties, is of clear interest. Similarly, systematic μSR measurements, crucial for detecting the presence of TRS breaking in Re-based superconductors, are in high demand. For instance, although both Re-Zr and Re-Nb alloys exhibit rich crystal structures and superconducting phase diagrams, only their α-Mn-type phase has been explored. In addition, most of the original measurements were performed only on polycrystalline samples. Hence, the synthesis of single crystals will be essential in the comprehensive search for possible superconducting nodes and, thus, for mixed singlet-triplet pairing. Finally, it would be of interest to extend the μSR studies on elementary rhenium from the bulk-to its thin-film form, where inversion symmetry is artificially broken. By checking if the TRS breaking is maintained or not, will help us to further clarify the rhenium conundrum.

 Author Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

 Funding

This work was supported by the start funding from East China Normal University (ECNU), the Swiss National Science Foundation (Grant No. 200021-169455) and the Sino-Swiss Science and Technology Cooperation (Grant No. IZLCZ2-170075).

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

We thank M. Shi for the fruitful discussion. We thank M. Medarde for the assistance during the electrical resistivity and magnetization measurements, and D. J. Gawryluk and E. Pomjakushina for synthesizing the materials. We acknowledge the allocation of beam time at the Swiss muon source (SμS) (Dolly, GPS, and LTF spectrometers).

 References Bauer E Sigrist M ., editors. Non-Centrosymmetric Superconductors, 847. Berlin: Springer-Verlag (2012). Smidman M Salamon MB Yuan HQ Agterberg DF . Superconductivity and Spin–Orbit Coupling in Non-centrosymmetric Materials: A Review. Rep Prog Phys (2017) 80:036501. 10.1088/1361-6633/80/3/036501 Ghosh SK Smidman M Shang T Annett JF Hillier AD Quintanilla J Recent Progress on Superconductors with Time-Reversal Symmetry Breaking. J Phys Condens Matter (2021) 33:033001. 10.1088/1361-648x/abaa06 Kim H Wang K Nakajima Y Hu R Ziemak S Syers P Beyond Triplet: Unconventional Superconductivity in a Spin-3/2 Topological Semimetal. Sci Adv (2018) 4. 10.1126/sciadv.aar796910.1126/sciadv.aao4513 eaao4513 Sun Z Enayat M Maldonado A Lithgow C Yelland E Peets DC Dirac Surface States and Nature of Superconductivity in Noncentrosymmetric BiPd. Nat Commun (2015) 6:6633. 10.1038/ncomms7633 Ali MN Gibson QD Klimczuk T Cava RJ . Noncentrosymmetric Superconductor with a Bulk Three-Dimensional Dirac Cone Gapped by Strong Spin-Orbit Coupling. Phys Rev B (2014) 89:020505. 10.1103/PhysRevB.89.020505 Sato M Fujimoto S . Topological Phases of Noncentrosymmetric Superconductors: Edge States, Majorana Fermions, and Non-abelian Statistics. Phys Rev B (2009) 79:094504. 10.1103/PhysRevB.79.094504 Tanaka Y Mizuno Y Yokoyama T Yada K Sato M . Anomalous Andreev Bound State in Noncentrosymmetric Superconductors. Phys Rev Lett (2010) 105:097002. 10.1103/PhysRevLett.105.097002 Sato M Ando Y . Topological Superconductors: A Review. Rep Prog Phys (2017) 80:076501:076501. 10.1088/1361-6633/aa6ac7 Qi X-L Zhang S-C . Topological Insulators and Superconductors. Rev Mod Phys (2011) 83:1057–110. 10.1103/RevModPhys.83.1057 Kallin C Berlinsky J . Chiral Superconductors. Rep Prog Phys (2016) 79:054502:054502. 10.1088/0034-4885/79/5/054502 Bauer E Hilscher G Michor H Paul C Scheidt EW Gribanov A Heavy Fermion Superconductivity and Magnetic Order in Noncentrosymmetric CePt3Si. Phys Rev Lett (2004) 92:027003. 10.1103/PhysRevLett.92.027003 Muro Y Eom D Takeda N Ishikawa M . Contrasting Kondo-Lattice Behavior in CeTSi3 and CeTGe3 (T = Rh and Ir). J Phys Soc Jpn (1998) 67:3601–4. 10.1143/JPSJ.67.3601 Hillier AD Quintanilla J Cywinski R . Evidence for Time-Reversal Symmetry Breaking in the Noncentrosymmetric Superconductor LaNiC2 . Phys Rev Lett (2009) 102:117007. 10.1103/PhysRevLett.102.117007 Barker JAT Singh D Thamizhavel A Hillier AD Lees MR Balakrishnan G Unconventional Superconductivity in La7Ir3 Revealed by Muon Spin Relaxation: Introducing a New Family of Noncentrosymmetric Superconductor that Breaks Time-Reversal Symmetry. Phys Rev Lett (2015) 115:267001. 10.1103/PhysRevLett.115.267001 Shang T Smidman M Wang A Chang L-J Baines C Lee MK Simultaneous Nodal Superconductivity and Time-Reversal Symmetry Breaking in the Noncentrosymmetric Superconductor CaPtAs. Phys Rev Lett (2020) 124:207001. 10.1103/PhysRevLett.124.207001 Singh RP Hillier AD Mazidian B Quintanilla J Annett JF Paul DM Detection of Time-Reversal Symmetry Breaking in the Noncentrosymmetric Superconductor Re6Zr Using Muon-Spin Spectroscopy. Phys Rev Lett (2014) 112:107002. 10.1103/PhysRevLett.112.107002 Singh D Barker JAT Thamizhavel A Paul DM Hillier AD Singh RP . Time-reversal Symmetry Breaking in the Noncentrosymmetric Superconductor Re6Hf : Further Evidence for Unconventional Behavior in the α-Mn Family of Materials. Phys Rev B (2017a) 96:180501. 10.1103/PhysRevB.96.180501 Shang T Pang GM Baines C Jiang WB Xie W Wang A Nodeless Superconductivity and Time-Reversal Symmetry Breaking in the Noncentrosymmetric Superconductor Re24Ti5 . Phys Rev B (2018) 97:020502. 10.1103/PhysRevB.97.020502 Shang T Smidman M Ghosh SK Baines C Chang LJ Gawryluk DJ Time-reversal Symmetry Breaking in Re-based Superconductors. Phys Rev Lett (2018) 121:257002. 10.1103/PhysRevLett.121.257002 Singh D Scheurer MS Hillier AD Adroja DT Singh RP . Time-reversal-symmetry Breaking and Unconventional Pairing in the Noncentrosymmetric Superconductor La7Rh3 . Phys Rev B (2020) 102:134511. 10.1103/PhysRevB.102.134511 Shang T Ghosh SK Zhao JZ Chang L-J Baines C Lee MK Time-reversal Symmetry Breaking in the Noncentrosymmetric Zr3Ir superconductorIr Superconductor. Phys Rev B (2020) 102:020503. 10.1103/PhysRevB.102.020503 Aczel AA Williams TJ Goko T Carlo JP Yu W Uemura YJ Muon Spin Rotation/relaxation Measurements of the Noncentrosymmetric Superconductor Mg10Ir19B16 . Phys Rev B (2010) 82:024520. 10.1103/PhysRevB.82.024520 Singh D Barker JAT Thamizhavel A Hillier AD Paul DM Singh RP . Superconducting Properties and μSR Study of the Noncentrosymmetric Superconductor Nb0.5Os0.5 . J Phys Condens Matter (2018) 30:075601:075601. 10.1088/1361-648X/aaa376 Biswas PK Hillier AD Lees MR Paul DM . Comparative Study of the Centrosymmetric and Noncentrosymmetric Superconducting Phases of Re3W Using Muon Spin Spectroscopy and Heat Capacity Measurements. Phys Rev B (2012) 85:134505. 10.1103/PhysRevB.85.134505 Shang T Gawryluk DJ Verezhak JAT Pomjakushina E Shi M Medarde M Structure and Superconductivity in the Binary Re1−xMox Alloys. Phys Rev Mater (2019) 3:024801. 10.1103/PhysRevMaterials.3.024801 Shang T Baines C Chang L-J Gawryluk DJ Pomjakushina E Shi M Re1−xMox as an Ideal Test Case of Time-Reversal Symmetry Breaking in Unconventional Superconductors. npj Quan Mater. (2020) 5:76. 10.1038/s41535-020-00279-1 Amato A . Heavy-fermion Systems Studied by μSR Technique. Rev Mod Phys (1997) 69:1119–80. 10.1103/RevModPhys.69.1119 Blundell SJ . Spin-polarized Muons in Condensed Matter Physics. Contemp Phys (1999) 40:175–92. 10.1080/001075199181521 Brewer JH . Muon Spin Rotation/relaxation/resonance. In: Trigg GL ., editor. Digital Encyclopedia of Applied Physics. Weinheim: Wiley VCH eap258. (2003). 10.1002/3527600434.eap258 Yaouanc A de Réotier PD . Muon Spin Rotation, Relaxation, and Resonance: Applications to Condensed Matter. Oxford: Oxford University Press (2011). Sonier JE Brewer JH Kiefl RF . μSR Studies of the Vortex State in Type-II Superconductors. Rev Mod Phys (2000) 72:769–811. 10.1103/RevModPhys.72.769 Barford W Gunn JMF . The Theory of the Measurement of the London Penetration Depth in Uniaxial Type II Superconductors by Muon Spin Rotation. Physica C: Superconductivity (1988) 156:515–22. 10.1016/0921-4534(88)90014-7 Brandt EH . Properties of the Ideal Ginzburg-Landau Vortex Lattice. Phys Rev B (2003) 68:054506. 10.1103/PhysRevB.68.054506 Maisuradze A Khasanov R Shengelaya A Keller H . Comparison of Different Methods for Analyzing μSR Line Shapes in the Vortex State of Type-II Superconductors. J Phys Condens Matter 21, 075701 (2009) 075701. 10.1088/0953-8984/21/7/075701 Tinkham M . Introduction to Superconductivity. 2nd ed. Mineola, NY: Dover Publications) (1996). Carrington A Manzano F . Magnetic Penetration Depth of MgB2 . Physica C: Superconductivity (2003) 385:205–14. 10.1016/S0921-4534(02)02319-5 Kubo R Toyabe T . A Stochastic Model for Low-Field Resonance and Relaxation. In: Blinc R ., editor. Magnetic Resonance and Relaxation. Proceedings of the XIVth Colloque Ampère. Amsterdam: North-Holland) (1967). p. 810–23. Massalski TB Okamoto H Kacprzak L Subramanian PR . Binary Alloy Phase Diagrams . 2 edn. Materials Park, OH: ASM International (1996). Roberts BW . Survey of Superconductive Materials and Critical Evaluation of Selected Properties. J Phys Chem Reference Data (1976) 5:581–822. 10.1063/1.555540 Murray JL . The Re−Ti (Rhenium-Titanium) System. J Phase Equilib (1982) 2:462–6. 10.1007/BF02876164 Philip TV Beck PA . CsCl-type Ordered Structures in Binary Alloys of Transition Elements. JOM J Min Met Mat S (1957) 9:1269–71. 10.1007/BF03398305 Singh D K. P. S Barker JAT Paul DM Hillier AD Singh RP . Time-reversal Symmetry Breaking in the Noncentrosymmetric Superconductor Re6Ti. Phys Rev B (2018b) 97:100505. 10.1103/PhysRevB.97.100505 Jorda JL Muller J . The Vanadium-Rhenium System: Phase Diagram and Superconductivity. J Less Common Met (1986) 119:337–45. 10.1016/0022-5088(86)90694-6 Eremenko VN Velikanova T . Intrusion Phases Based on Metallides in Ternary Systems of Transition Metals with Carbon. Sov Prog Chem (1990) 56:21–8. Matano K Yatagai R Maeda S Zheng G-q. Full-gap Superconductivity in Noncentrosymmetric Re6Zr, Re27Zr5, and Re24Zr5 . Phys Rev B (2016) 94:214513. 10.1103/PhysRevB.94.214513 Giorgi AL Szklarz EG . Superconductivity and Lattice Parameters of the Dirhenides and Ditechnides of Thorium, Hafnium and Zirconium. J Less Common Met (1970) 22:246–8. 10.1016/0022-5088(70)90027-5 Cirillo C Fittipaldi R Smidman M Carapella G Attanasio C Vecchione A Evidence of Double-Gap Superconductivity in Noncentrosymmetric Nb0.18Re0.82 Single Crystals. Phys Rev B (2015) 91:134508. 10.1103/PhysRevB.91.134508 Chen J Jiao L Zhang JL Chen Y Yang L Nicklas M BCS-like Superconductivity in the Noncentrosymmetric Compounds NbxRe1−x(0.13≤x≤0.38). Phys Rev B (2013) 88:144510. 10.1103/PhysRevB.88.144510 Karki AB Xiong YM Haldolaarachchige N Stadler S Vekhter I Adams PW Physical Properties of the Noncentrosymmetric Superconductor Nb0.18Re0.82 . Phys Rev B (2011) 83:144525. 10.1103/PhysRevB.83.144525 Lue CS Su TH Liu HF Young B-L . Evidence Fors-Wave Superconductivity in Noncentrosymmetric Re24Nb5 from 93Nb NMR Measurements. Phys Rev B (2011) 84:052509. 10.1103/PhysRevB.84.052509 Chen B Guo Y Wang H Su Q Mao Q Du J Superconductivity in the Noncentrosymmetric Compound Re6Hf. Phys Rev B (2016) 94:024518. 10.1103/PhysRevB.94.024518 Singh D Hillier AD Thamizhavel A Singh RP . Superconducting Properties of the Noncentrosymmetric Superconductor Re6Hf. Phys Rev B (2017) 96:064521. 10.1103/PhysRevB.96.06452110.1103/physrevb.96.180501 Barker JAT Breen BD Hanson R Hillier AD Lees MR Balakrishnan G Superconducting and Normal-State Properties of the Noncentrosymmetric Superconductor Re3Ta. Phys Rev B (2018) 98:104506. 10.1103/PhysRevB.98.104506 Arushi SD Singh D Biswas PK Hillier AD Singh RP . Unconventional Superconducting Properties of Noncentrosymmetric Re5.5Ta. Phys Rev B (2020) 101:144508. 10.1103/PhysRevB.101.144508 Mamiya T Nomura K Masuda Y . Superconductivity of Tantalum-Rhenium Alloys. J Phys Soc Jpn (1970) 28:380–9. 10.1143/JPSJ.28.380 Biswas PK Lees MR Hillier AD Smith RI Marshall WG Paul DM . Structure and Superconductivity of Two Different Phases of Re3W. Phys Rev B (2011) 84:184529. 10.1103/PhysRevB.84.184529 Kimura N Ito K Aoki H Uji S Terashima T . Extremely High Upper Critical Magnetic Field of the Noncentrosymmetric Heavy Fermion Superconductor CeRhSi3 . Phys Rev Lett (2007) 98:197001. 10.1103/PhysRevLett.98.197001 Sugitani I Okuda Y Shishido H Yamada T Thamizhavel A Yamamoto E Pressure-induced Heavy-Fermion Superconductivity in Antiferromagnet CeIrSi3 without Inversion Symmetry. J Phys Soc Jpn (2006) 75:043703. 10.1143/JPSJ.76.05100910.1143/jpsj.75.043703 Carnicom EM Xie W Klimczuk T Lin J Górnicka K Sobczak Z TaRh2B2 and NbRh2B2: Superconductors with a Chiral Noncentrosymmetric Crystal Structure. Sci Adv (2018) 4:eaar7969:eaar7969. 10.1126/sciadv.aar7969 Singh D K. P. S Marik S Hillier AD Singh RP . Superconducting and Normal State Properties of the Noncentrosymmetric Superconductor NbOs2 . Investigated by Muon Spin Relaxation and Rotation. Phys Rev B (2019) 99:014516. 10.1103/PhysRevB.99.014516 Singh D Sajilesh KP Marik S Hillier AD Singh RP . Superconducting Properties of the Noncentrosymmetric Superconductor TaOs. Supercond Sci Technol (2017) 30:125003. 10.1088/1361-6668/aa8f8e Mayoh DA Barker JAT Singh RP Balakrishnan G Paul DM Lees MR . Superconducting and Normal-State Properties of the Noncentrosymmetric Superconductor Re6Zr. Phys Rev B (2017) 96:064521. 10.1103/PhysRevB.96.064521 Pang GM Nie ZY Wang A Singh D Xie W Jiang WB Fully Gapped Superconductivity in Single Crystals of Noncentrosymmetric Re6Zr with Broken Time-Reversal Symmetry. Phys Rev B (2018) 97:224506. 10.1103/PhysRevB.97.224506 Parab P Singh D Haram S Singh RP Bose S . Point Contact Andreev Reflection Studies of a Non-centro Symmetric Superconductor Re6Zr. Sci Rep (2019) 9:2498. 10.1038/s41598-019-39160-y Yuan HQ Agterberg DF Hayashi N Badica P Vandervelde D Togano K S-wave Spin-Triplet Order in Superconductors without Inversion Symmetry:Li2Pd3B and Li2Pt3B. Phys Rev Lett (2006) 97:017006. 10.1103/PhysRevLett.97.017006 Nishiyama M Inada Y Zheng G-q. Spin Triplet Superconducting State Due to Broken Inversion Symmetry in Li2Pt3B. Phys Rev Lett (2007) 98:047002. 10.1103/PhysRevLett.98.047002 Xie W Zhang P Shen B Jiang W Pang G Shang T CaPtAs: a New Noncentrosymmetric Superconductor. Sci China Phys Mech Astron (2020) 63:237412. 10.1007/s11433-019-1488-5 Samokhin KV Zijlstra ES Bose SK . CePt3Si: An Unconventional Superconductor without Inversion Center. Phys Rev B (2004) 69:094514. 10.1103/PhysRevB.69.094514 Suetin DV Ivanovskii AL . Comparative Study of Electronic Structure of Cubic and Hexagonal Phases of Re3W as Non-centrosymmetric and Centrosymmetric Low-T c Superconductors. Intermetallics (2013) 34:101–5. 10.1016/j.intermet.2012.11.015 Winiarski MJ . Electronic Structure of Non-centrosymmetric Superconductors Re24(Nb;Ti)5 by Ab Initio Calculations. J Alloys Comp (2014) 616:1–4. 10.1016/j.jallcom.2014.07.081 Mojammel AK Karki AB Samanta T Browne D Stadler S Vekhter I Complex Superconductivity in the Noncentrosymmetric Compound Re6Zr. Phys Rev B (2016) 94:144515. 10.1103/PhysRevB.94.144515 Shang T Philippe J Verezhak JAT Guguchia Z Zhao JZ Chang L-J Nodeless Superconductivity and Preserved Time-Reversal Symmetry in the Noncentrosymmetric Mo3P Superconductor. Phys Rev B (2019) 99:184513. 10.1103/PhysRevB.99.184513 Yip S . Noncentrosymmetric Superconductors. Annu Rev Condens Matter Phys (2014) 5:15–33. 10.1146/annurev-conmatphys-031113-133912 Luke GM Keren A Le LP Wu WD Uemura YJ Bonn DA Muon Spin Relaxation in UPt3 . Phys Rev Lett (1993) 71:1466–9. 10.1103/PhysRevLett.71.1466 Luke GM Fudamoto Y Kojima KM Larkin MI Merrin J Nachumi B Time-reversal Symmetry-Breaking Superconductivity in Sr2RuO4 . Nature (1998) 394:558–61. 10.1038/29038 Xia J Maeno Y Beyersdorf PT Fejer MM Kapitulnik A . High Resolution Polar Kerr Effect Measurements of Sr2RuO4: Evidence for Broken Time-Reversal Symmetry in the Superconducting State. Phys Rev Lett (2006) 97:167002. 10.1103/PhysRevLett.97.167002 Aoki Y Tsuchiya A Kanayama T Saha SR Sugawara H Sato H Time-Reversal Symmetry-Breaking Superconductivity in Heavy-fermion PrOs4Sb12 Detected by Muon-Spin Relaxation. Phys Rev Lett (2003) 91:067003. 10.1103/PhysRevLett.91.067003 Schemm ER Gannon WJ Wishne CM Halperin WP Kapitulnik A . Observation of Broken Time-Reversal Symmetry in the Heavy-Fermion Superconductor UPt3 . Science (2014) 345:190–3. 10.1126/science.1248552 Hillier AD Quintanilla J Mazidian B Annett JF Cywinski R . Nonunitary Triplet Pairing in the Centrosymmetric Superconductor LaNiGa2 . Phys Rev Lett (2012) 109:097001. 10.1103/PhysRevLett.109.097001 Ran S Eckberg C Ding Q-P Furukawa Y Metz T Saha SR Nearly Ferromagnetic Spin-Triplet Superconductivity. Science (2019) 365:684–7. 10.1126/science.aav8645 Ishida K Mukuda H Kitaoka Y Asayama K Mao ZQ Mori Y Spin-triplet Superconductivity in Sr2RuO4 Identified by 17O Knight Shift. Nature (1998) 396:658–60. 10.1038/25315 Tou H Kitaoka Y Ishida K Asayama K Kimura N O¯nuki Y Nonunitary Spin-Triplet Superconductivity in UPt3: Evidence from 195Pt Knight Shift Study. Phys Rev Lett (1998) 80:3129–32. 10.1103/PhysRevLett.80.3129 Mackenzie AP Maeno Y . The Superconductivity of Sr2RuO4 and the Physics of Spin-Triplet Pairing. Rev Mod Phys (2003) 75:657–712. 10.1103/RevModPhys.75.657 Joynt R Taillefer L . The Superconducting Phases of UPt3 . Rev Mod Phys (2002) 74:235–94. 10.1103/RevModPhys.74.235 Lee W-C Zhang S-C Wu C . Pairing State with a Time-Reversal Symmetry Breaking in FeAs-Based Superconductors. Phys Rev Lett (2009) 102:217002. 10.1103/PhysRevLett.102.217002 Sharma S Motla K Beare J Nugent M Pula M Munsie T Fully Gapped Superconductivity in Centrosymmetric and Noncentrosymmetric Re-B Compounds Probed with μSR. Phys Rev B (2021) 103:104507. 10.1103/PhysRevB.103.104507 Shang T Amon A Kasinathan D Xie W Bobnar M Chen Y Enhanced T c and Multiband Superconductivity in the Fully-Gapped ReBe22 Superconductor. New J Phys (2019) 21:073034:073034. 10.1088/1367-2630/ab307b Weng ZF Zhang JL Smidman M Shang T Quintanilla J Annett JF Two-Gap Superconductivity in LaNiGa2 with Nonunitary Triplet Pairing and Even Parity Gap Symmetry. Phys Rev Lett (2016) 117:027001. 10.1103/PhysRevLett.117.027001 Quintanilla J Hillier AD Annett JF Cywinski R . Relativistic Analysis of the Pairing Symmetry of the Noncentrosymmetric Superconductor LaNiC2 . Phys Rev B (2010) 82:174511. 10.1103/PhysRevB.82.174511 Yang S Wang C Sahin H Chen H Li Y Li S-S Tuning the Optical, Magnetic, and Electrical Properties of ReSe2 by Nanoscale Strain Engineering. Nano Lett (2015) 15:1660–6. 10.1021/nl504276u Kochat V Apte A Hachtel JA Kumazoe H Krishnamoorthy A Susarla S Re Doping in 2D Transition Metal Dichalcogenides as a New Route to Tailor Structural Phases and Induced Magnetism. Adv Mater (2017) 29:1703754. 10.1002/adma.201703754 Cenzual K Parthé E Waterstrat RM . Zr21Re25, a New Rhombohedral Structure Type Containing 12Å-Thick Infinite MgZn2 (Laves)-type Columns. Acta Crystallogr C (1986) 42:261–6. 10.1107/S0108270186096555