Spectroscopic signatures of time-reversal symmetry breaking superconductivity

The collective mode spectrum of a symmetry-breaking state, such as a superconductor, provides crucial insight into the nature of the order parameter. In this work, we study two collective modes which are unique to unconventional superconductors that spontaneously break time reversal symmetry. We show that these modes are coherent and underdamped for a wide variety of time-reversal symmetry breaking superconducting states. By further demonstrating that these modes can be detected using a number of existing experimental techniques, we propose that our work can be leveraged as a form of “collective mode spectroscopy” that drastically expands the number of experimental probes capable of detecting time-reversal symmetry breaking in unconventional superconductors. Variations in unconventional superconductors are increasing in number and diversity in recent years, and so it is increasingly important to develop analytical methods to distinguish and categorise their different features. Here, the authors study two types of collective modes of the order parameter unique to time-reversal symmetry breaking superconductors and propose using them to identify these exotic states.

T here is a rich and constantly expanding taxonomy of unconventional superconducting phases, and increasingly sophisticated probes are needed to distinguish one such phase from another. A defining feature of any superconducting state is its collective mode spectrum, which encodes the dynamics of the order parameter. The collective modes of conventional (phonon-mediated) superconductors are well-established, consisting of two modes which correspond to fluctuations in the amplitude and phase of the order parameter. The first, called the Higgs mode 1,2 , is massive and resides at the edge of the quasiparticle continuum 3,4 , while the second, called the Anderson-Bogoliubov-Goldstone (ABG) mode 5,6 , is massless in a neutral system, in accordance with Goldstone's theorem, but is lifted to the plasma frequency in the presence of long-ranged Coulomb interactions 7 , and is thus indistinguishable from the usual plasmon in real materials.
However, systems with more complex order parameters can exhibit a rich collective mode spectrum featuring other modes, such as additional Higgs modes in anisotropically gapped (e.g., d-wave) superconductors 8 , or Leggett modes in multi-band superconductors 9 . The presence of these additional modes in the spectrum can then be taken as a fingerprint of the underlying order parameter symmetry. That is, the dynamics of the order parameter can be studied to gain insight into its equilibrium structure. Such a scheme has only recently been proposed in the context of anisotropically gapped superconductors, where the spectrum of non-equilibrium Higgs modes can be used to deduce the orbital symmetry of the order parameter 10,11 .
In this article, we generalize this notion of "collective mode spectroscopy" to a particularly exotic class of unconventional superconductors, namely those which spontaneously break timereversal symmetry in addition to global U(1) symmetry at the superconducting transition. These time-reversal symmetry breaking (TRSB) superconducting states are the subject of considerable current interest, and are believed to be realized in a number of bulk materials, including Sr 2 RuO 4 12,13 , UPt 3 14,15 , URu 2 Si 2 16 , UTe 2 17,18 , PrOs 4 Sb 12 19,20 , and K-doped BaFe 2 As 2 21 , as well as engineered structures such as Bi/Ni bilayers 22 and SnTe nanowires 23 . Moreover, recent theoretical proposals have suggested that such states may also be realized in moiré heterostructures 24,25 , and could be generically engineered in twisted bilayers of anisotropically gapped superconductors 26,27 .
In what follows, we identify two collective modes, the "generalized clapping modes," which are unique to TRSB superconductors and subsequently derive their spectrum from a generally applicable microscopic weak-coupling theory. Our first key finding is that these modes are coherent collective excitations for a wide variety of TRSB states, even those which have point or line nodes in the superconducting gap. Owing to the universality of these modes in TRSB superconductors, they may serve as clear spectroscopic signatures which identify a TRSB superconducting state. To this end, we discuss a variety of existing experimental probes that can couple to these modes, and hence can be used as means to detect TRSB superconductivity in quantum materials. We pay special attention to the superconducting state of Sr 2 RuO 4 , where the generalized clapping mode spectrum could distinguish the two current leading candidate order parameters.

Results
Generalized clapping modes. TRSB superconducting states are characterized by a doubly-degenerate complex multi-component order parameter of the form Δ = Δ 1 ± iΔ 2 (see supplementary Note 1), and can be divided into two classes: (1) systems where both components Δ 1 and Δ 2 belong to the same multidimensional irreducible representation (irrep) of the crystalline point group, in which case |Δ 1 | = |Δ 2 | is required by symmetry, and (2) "mixed symmetry" systems where the two components belong to different irrep's, which can arise due to either an accidental degeneracy between two pairing channels or two successive superconducting transitions (see Supplementary Note 1), in which case the ratio |Δ 1 |/|Δ 2 | is unconstrained.
Given this internal orbital structure of the order parameter, we would expect that in addition to the usual ABG and Higgs modes, there should be a massive mode corresponding to fluctuations of the relative phase between the two order parameter components around its equilibrium value of ±π/2. Further, we expect that there should be a second amplitude mode which corresponds to fluctuations of the relative amplitude |Δ 1 |/|Δ 2 | so there are a total of four real modes.
Historically, similar modes corresponding to "internal vibrations of the structure of the order parameter" were first recognized in 3 He-A 28 , one of which was named the "clapping mode." This was subsequently extended to the two-dimensional chiral p-wave (p + ip) superconductor 29 , where there are two such clapping modes which are degenerate and reside at a frequency Ω ¼ ffiffi ffi 2 p Δ 0 , with Δ 0 the magnitude of the order parameter.
As argued above, we can anticipate analogous modes, which we call the generalized clapping modes, for any multi-component TRSB state on symmetry grounds alone 30 . Although it is intuitively obvious that the generalized modes should exist in principle, there is no reason to expect a priori that they are not overdamped by quasiparticle excitations. For a fully gapped superconductor, this requires that the frequencies of both modes lie below the quasiparticle continuum or, in the case of a nodal quasiparticle gap function, requires that the spectral function of each mode retains sharp features despite the presence of nodal quasiparticles. To establish that the generalized clapping modes are coherent, and thus experimentally detectable, for a generic TRSB state requires a derivation of their spectrum starting from a microscopic theory, which we furnish below. By studying the generalized clapping mode spectrum for a wide variety of TRSB order parameters we are able to establish the general features of these modes, which ultimately enable their use in the collective mode spectroscopy of real materials.
Weak-coupling theory. We begin with a single band of fermions subject to the attractive interaction V kk 0 ¼ À∑ '¼1;2 g ' χ ' k χ ' k 0 where g ℓ > 0 are coupling constants and χ ' k are form factors which encode the orbital symmetry of the interaction. We take these to be real and normalized according to the inner product We assume pairing in the S z tot ¼ 0 sector, but within this sector our results are applicable to both singlet and m = 0 triplet pairing. We treat this system within the imaginarytime path-integral formalism by introducing a Hubbard-Stratonovich decoupling field Δ ℓ in each pairing channel and integrating out the fermions. One then arrives at the effective action for the order parameters Δ ℓ of The inverse fermion propagator is G À1 Àq χ ' k τ À where τ i are the Pauli matrices in Nambu space, τ ± ¼ 1 2 τ x ± iτ y , and ξ k = k 2 /2m − μ is the single-particle energy measured with respect to the Fermi level. We have also combined fermionic/bosonic Matsubara frequencies and momenta into the four-vectors k = (iω n , k) and q = (iΩ m , q), where q corresponds to the center-of-mass momentum of the fermion pair and k corresponds to the relative momentum.
We assume g 1 , g 2 are such that we find a saddle-point configuration in which both Δ (1) and Δ (2) are condensed with a relative phase of π/2, breaking time-reversal symmetry as discussed above. It will be convenient to change basis from (Δ (1) , Δ (2) and η 1,2 quantify the relative magnitude of each order parameter component. We choose to normalize them such that η 2 1 þ η 2 2 ¼ 1 and can express both in terms of a "mixing angle" as η 1 ¼ cos η and η 2 ¼ sin η.
Expanding around the saddle point with Δ þ q ¼ Δ 0 ¼ const. and Δ À q ¼ 0, the mean-field equations are where E 2 k ¼ ξ 2 k þ jΔ k j 2 and the angle-dependent quasiparticle gap function is Δ k ¼ Δ 0 χ þ k . Given a particular set of pairing symmetries and coupling constants, these equations can be solved to determine the equilibrium values of η and Δ 0 which characterize the condensate.
Effective action for fluctuations. Now, we move on to consider the fluctuations around this saddle point, which we parameterize as where θ is the ABG phase mode, h is the Higgs mode, and the a and b modes are fluctuations in the relative amplitude and phase of the two order parameter components, i.e., the generalized clapping modes (see Fig. 1a). This parameterization suggests that we can equivalently think of the generalized clapping modes as being fluctuations in the degenerate time-reversed Δ − pairing channel, which are formally similar to Bardasis-Schreiffer modes 31,32 , as illustrated in Fig. 1b. To organize our calculations in a manifestly gauge-invariant way, we minimally couple the system to an external (classical) gauge field and perform a unitary transformation G À1 ! UG À1 U y with U ¼ e Àiθτ z such that the ABG mode and gauge field only appear together as the gauge-invariant vector field V 0 = A 0 + ∂ τ θ and V = A − ∂θ (where we have set the electron charge equal to one). To quadratic order in these fields and setting q → 0, the action is In the above, Π 00 is the electronic compressibility, n ij s is the superfluid density, D h;a;b are the propagators for the Higgs, relative amplitude, and relative phase modes, andΠ ha and Π 0b are linear couplings between the Higgs/relative amplitude and ABG/ relative phase modes which are non-vanishing in the q → 0 limit. This action is derived and the correlation functions which appear in it are evaluated in Supplementary Note 2.
To gain intuition, we begin by considering the familiar case of p + ip pairing in two dimensions. When the p x and p y components of the order parameter occur with equal amplitudes (η = π/4), as is dictated by symmetry in most cases of physical interest, we find that after analytically continuing to real-time the propagators for both the a and b modes have a pole at Ω ¼ ffiffi ffi 2 p Δ 0 , i.e., the two modes are degenerate. So, we see that the usual clapping modes previously studied in p + ip superconductors and 3 He-A are indeed a special case of the generalized clapping modes a and b studied in this work. Moreover, at this point the gap is isotropic and the couplingsΠ ha and Π 0b vanish so that the (generalized) clapping modes decouple from both quasiparticle excitations and other collective modes and thus are infinitely long-lived at zero temperature.
In real materials, anisotropic crystal fields can lead to deviations from the equal-amplitude η = π/4 state 33 , making it interesting to consider the collective mode spectrum of the p + ip state for general mixing angles. We plot the spectral functions of the a and b modes in Fig. 2a, b as a function of the mixing angle η and frequency Ω. We see that the two modes split as the mixing angle deviates from π/4, so that for generic mixing angles there are two generalized clapping modes corresponding to the relative amplitude and phase fluctuations.
Two-dimensional systems. We will now investigate the generalized clapping mode spectrum for several even parity two-component TRSB order parameters which are potentially relevant to experimental systems: the d x 2 Ày 2 þ id xy state, originally studied in the context of the cuprate high-temperature superconductors 34,35 , and now the subject of renewed interest due to its potential relevance to a number of moiré systems 24,25,36 , other heterostructures 37 , and its proposed realization in twisted bilayers of cuprates 26 and other unconventional superconductors 27 ; the s + id xy state, which has long been of interest in relation to the iron pnictide high-temperature superconductors [38][39][40] ; and the d x 2 Ày 2 þ ig xyðx 2 Ày 2 Þ state which has recently been proposed as the order parameter of Sr 2 RuO 4 41,42 . The basic properties of each order parameter studied in this work are listed in Table 1.
Because these are all mixed symmetry states where the mixing angle is unconstrained by point group symmetries, it is important to survey the generalized clapping mode spectra over the full range of η. We first turn our attention to the relative phase mode, the spectra of which we plot for each of the above order parameters in Fig. 3a-f. Crucially, we observe that the relative phase mode is a well-defined, coherent mode for all values of the mixing angle, with a sharp spectral function that can be seen in the line cuts in Fig. 3d-f. This even remains true for nodal order parameters, e.g., the d + ig state.
Next, we turn our attention to the relative amplitude mode. As we show in Supplementary Note 2, this mode is coupled to the Higgs mode even at zero momentum. The relative amplitude and Higgs modes then hybridize to form two orthogonal amplitude modes, which we call A + and A − . We calculate the propagators D A± (Ω) for these modes in Supplementary Note 2, which we use to plot the spectral functions of each mode for various pairing symmetries, as shown in Fig. 4a-f.
For all of the pairing symmetries studied, the A − mode resides at or slightly below the gap edge, much like the conventional Higgs mode in a single-component superconductor. More interesting is the A + mode, which lies well below the quasiparticle continuum for a wide range of mixing angles. This low-frequency amplitude mode represents a second novel collective excitation characteristic of TRSB superconducting states.
We also note that at η = π/4, the d þ id 0 state is fully gapped and chiral, and both generalized clapping modes are degenerate with a frequency of Ω ¼ ffiffi ffi 2 p Δ 0 , much like the chiral p-wave state in Fig. 2 43 . As we will show below, this degeneracy is a general feature of chiral order parameters.
Three-dimensional systems. We now consider several TRSB superconducting states in three-dimensional systems: the 3d p x + ip y state (i.e., a p x + ip y order parameter defined over a spherical Fermi surface), better known as the Anderson-Brinkman-Morel (ABM) state of 3 He-A 44 , which is also a minimal model for a Weyl superconductor, and thus may bear some qualitative similarities to the superconducting state of UTe 2 18 and other candidate Weyl systems; the three-dimensional d x 2 Ày 2 þ id xy "double Weyl" superconducting state possibly realized in the "hidden order" phase of URu 2 Si 2 45,46 or SrPtAs 47,48 ; and finally the d xz + id yz state which is a candidate order parameter for Sr 2 RuO 4 as well as the hidden order phase of URu 2 Si 2 49,50 .
All of these states are most naturally considered as equaladmixtures, i.e., with η = π/4, in which case the generalized clapping modes decouple from the Higgs mode. As seen in Fig. 5a-d, the generalized clapping modes are degenerate with one another in all cases. Unlike the two-dimensional case, however, the mode frequency is not ffiffi ffi 2 p Δ 0 . By inspection of the analytic form of the generalized clapping mode propagators (see Supplementary Note 2), it is evident that the mode frequency is only ffiffi ffi 2 p Δ 0 if the order parameter is both chiral and isotropically gapped, whereas the ABM and 3d d þ id 0 (d + id) states exhibit point (line) nodes. We also note that this demonstrates that the generalized clapping modes remain coherent even for a system with line nodes in the superconducting gap, i.e., these modes' survival seems largely insensitive to the nodal structure of the order parameter. Combined with the prior results for the chiral p + ip and d þ id 0 states in two-dimensions, we see that chiral order parameters are generically characterized by degenerate generalized clapping modes. In this way, the generalized clapping mode spectrum can be used to differentiate chiral states from mixed symmetry TRSB states, as further discussed below.
Altogether, the results we have presented so far provide strong evidence that the generalized clapping modes are always Fig. 1 Generalized clapping modes. a Schematic of the collective mode spectrum for a time-reversal symmetry breaking superconductor, featuring the massless Anderson-Bogoliubov-Goldstone mode (acoustic plasmon), massive relative phase mode (optical plasmon), and two massive amplitude modes. The A + mode generically lies below the quasiparticle continuum. b Illustration of the collective modes: once the system condenses into e.g., the Δ 1 + iΔ 2 ground state, the corresponding order parameter Δ + resides in a Mexican hat potential V½Δ þ with a finite expectation value. The Higgs (ABG) mode is then an amplitude (phase) fluctuation in this condensed pairing channel. In contrast, the time-reversed Δ − channel remains undcondensed and experiences a parabolic potential V½Δ À centered around the minimum Δ − = 0. The generalized clapping modes are then the in-phase and out-of-phase fluctuations in this time-reversed uncondensed pairing channel. well-defined sub-gap excitations in TRSB superconducting states. A rigorous proof of this conjecture would be an interesting direction for future research, but is beyond the scope of the present work. At a practical level, this universality is essential to these modes' application in collective mode spectroscopy.
Having established that the relative phase and amplitude modes are typically well-defined collective excitations in TRSB superconductors, we now propose several means for their experimental detection.
Detection of relative phase modes. From the action (6), we see that the relative mode contributes to the charge density as ρ~i∂ t b due to the non-vanishing linear coupling Π 0b between b and the scalar potential. Since the relative phase mode resides at finite frequency, it is not affected by Coulomb interactions in the usual way and is not lifted up to the plasma frequency (similar to ref. 51 ).
One may then integrate b out of this action, renormalizing the electronic compressibility Π 00 (Ω) → Π 00 (Ω) + δΠ 00 (Ω) where is the relative phase mode's contribution to Π 00 (Ω). We plot this function for several pairing symmetries in Fig. 6a-d, where we see clear features at the relative phase mode frequency. Taking into account Coulomb interactions and integrating out the ABG phase mode, the fully renormalized longitudinal dielectric function is ϵ L ðΩ; qÞ ¼ 1 þ n s mΠ 00 ðΩÞ n s q 2 =m À Ω 2Π 00 ðΩÞ whereΠ 00 ðΩÞ ¼ Π 00 ðΩÞ þ δΠ 00 ðΩÞ is the renormalized compressibility and n s is the superfluid density. For frequencies much smaller than the plasma frequency, the dielectric function is ϵ L ðΩ; qÞ % 1 þΠ 00 ðΩÞ=q 2 , which includes a feature at the relative SrPtAs, URu 2 Si 2 For each order parameter we consider in this work, we list the corresponding irreducible representation (irrep) of the tetragonal point group, a representative basis function, the type of time-reversal symmetry breaking state, i.e., whether it belongs to a multi-dimensional irrep or is a "mixed symmetry" state (see Supplementary Note 1), the spatial dimension of the Fermi surface it exists on, and candidate materials where such a state is believed to be realized. phase mode frequency that can be detected by probes sensitive to the density-density response function, as discussed below. For most candidate TRSB superconductors, the superconducting gap, and hence the relative phase mode frequency, is in the terahertz regime. The ac electronic compressibility can be directly measured at THz frequencies using existing experimental techniques such as momentum-resolved electron energy loss spectroscopy (M-EELS) [52][53][54][55] , which enables direct experimental detection of the relative phase mode.
In such an experiment, the relative phase mode can be distinguished from trivial non-electronic modes by its dispersion and the fact that the peak corresponding to this mode should vanish above T c . It can be shown that the relative phase mode disperses as where Ω b is the frequency of the mode and α is a constant dependent on the mixing angle and orbital symmetries of the order parameter. In light of these considerations, the observation of a sub-gap peak in the charge response, measured via M-EELS, which disperses with an electronic-scale wavelength would constitute smokinggun evidence for the relative phase mode in a TRSB superconductor.
Alternatively, these modes may also be detected by sub-gap peaks in the microwave power absorption, as shown in refs. 56,57 , and via magneto-optical measurements 58 . In fact, a collective mode of unknown origin has been observed in microwave power absorption measurements of the heavy-fermion superconductor UBe 13 59 . To date, there is no evidence that the superconducting state of this system is TRSB, but a split transition has been reported in specific heat measurements of its Th-doped relative U 1−x Th x Be 13 60,61 , which is suggestive of a multi-component order. In light of this, it could be interesting to revisit the order parameter symmetry of these compounds and whether time- reversal symmetry is broken. If so, ref. 59 might represent the first measurement of a generalized clapping mode.
In addition to the equilibrium probes discussed above, the relative phase mode should also be detectable out of equilibrium using time-resolved THz spectroscopy 62 . In fact, a recent experiment has claimed to observe the Leggett phase mode in the multi-band superconductor MgB 2 using this technique 63 . However, another THz study of MgB 2 attributed the same experimental signatures to other modes, complicating the identification of the Leggett mode 64 . In any case, further work is necessary to establish the precise experimental signatures of the relative phase mode in such an experiment.
If one were to observe a peak near the relative phase mode frequency in any of the experiments listed above, its origin can be further distinguished with the application of an external magnetic field. For most tetragonal two-component superconductors, there is a symmetry-allowed linear coupling to an external magnetic In terms of the relative phase φ   6 ac compressibility. Imaginary part of the relative phase mode's contribution to the electronic compressibility, Im δΠ 00 , given by Eq. (7), plotted as a function of frequency Ω (scaled to the typical gap 2Δ 0 ) for a p + ip, b d þ id 0 , c s + id, and d d + ig orders, in units of e 2 ν. All exhibit clear features at the clapping mode frequency. For the p + ip and d þ id 0 cases, the clapping mode peak is small (and is indicated by the arrows in the figure) near η = π/4 due to the suppression of the coupling to the scalar potential at η = π/4. For the same two cases, the large shoulder at higher frequencies arises due to the relative phase mode's decay into, and coupling to, quasiparticles. We also note that within our model (where the basis functions are taken to be spherical harmonics), the results for the p + ip and d þ id 0 states are identical, which is an artifact of our highly symmetric model. In the s + id and d + ig cases, which are not subject to the same selection rules as the chiral states, the principal peaks are due to the relative phase mode.
between Δ 1 and Δ 2 , this can be written as δf ¼ gH sin φ % ÀgHδφ 2 where we have assumed that the fluctuations δφ are small around the equilibrium value of φ 0 = π/2. Then, on very general grounds 30 , this coupling shifts the relative phase mode frequency from Ω b,0 at zero field to Ω 2 b $ Ω 2 b;0 þ gH. That is, the frequency of an experimentally observed peak should vary linearly with small applied magnetic fields if it truly corresponds to the relative phase mode, and in this way, one may unambiguously assign an observed peak to the relative phase mode.
Detection of amplitude modes. The A ± modes couple to electromagnetic fields in the same manner as the Higgs mode in a conventional single-component superconductor, through the non-linear coupling δS~A 2 A ± . One method to detect amplitude modes using this coupling is via the observation of resonantly enhanced third-harmonic generation in a non-linear THz spectroscopy experiment, which has been successfully performed for both conventional 66 and high-temperature superconductors 67 . Amplitude modes have also been detected in THz pump, optical probe experiments for conventional superconductors 68 and cuprates 69 . However, despite the rapid development of THz spectroscopy of a probe of Higgs modes in superconductors, the impact of charge density fluctuations 70 , disorder [71][72][73] , and other modes 74 on the THz response continue to be actively discussed, and the identification of the Higgs mode in unconventional superconductors remains controversial.
The A + mode in TRSB superconductors is attractive on account of its low frequency, which should make it less heavily damped than conventional Higgs modes and easier to disentangle from the charge density fluctuations which onset at the gap edge 75 . Nonetheless, in light of the aforementioned controversy surrounding the THz detection of conventional Higgs modes, further work is necessary to establish the detailed THz response of the amplitude modes discussed in our work to facilitate comparison with potential future experiments.
Alternatively, the amplitude modes can be detected using microwave spectroscopy, or within linear response as a resonance in the optical conductivity in the presence of a background dc supercurrent 76 , as has been demonstrated in NbN 77 .
Like the relative phase mode discussed above, the observation of the low-frequency A + amplitude mode is a direct signature of the TRSB superconducting state. Thus, the detection of this mode represents a different avenue in the emerging field of "Higgs spectroscopy" 10,11,67 where non-equilibrium amplitude oscillations are used to gain insight into the symmetry of the condensate. However, in contrast to these pioneering works, generalized clapping modes can be detected in equilibrium and using a wide variety of experimental probes, as discussed above.

Discussion
To summarize, we have shown that fluctuations in the relative amplitude and phase of two-component TRSB order parameters are well-defined collective modes for a wide class of model order parameters with both spin-singlet and spin-triplet pairing (and possibly mixed-parity systems as well). Moreover, even for nodal gap functions, we have found that these modes are not overdamped by low-energy quasiparticles in the T → 0 limit. The frequency of each mode depends strongly on both the orbital symmetries of the two order parameter components and the relative amplitude of the two components in equilibrium.
Further, we have proposed a number of means to experimentally detect generalized clapping modes: the relative phase mode can be directly detected via measurement of the ac electronic compressibility (using e.g., M-EELS), the amplitude modes can be detected using ultrafast and non-linear THz spectroscopy as well as optical conductivity measurements (in the presence of a dc supercurrent), and both modes can be detected in microwave power absorption measurements. The observation of these modes in a given material would constitute robust evidence of a twocomponent TRSB order parameter.
Our work enables a variety of existing experimental techniques including ultrafast and non-linear optics, electron scattering, and microwave spectroscopy to be used as direct probes of TRSB in unconventional superconductors. These measurements can be thought of as a form of "collective mode spectroscopy," where one obtains information about the structure of the order parameter from its collective mode spectrum. Given the relative scarcity of probes directly sensitive to TRSB in superconductors (previously limited to only Kerr rotation and muon spin relaxation), this represents a substantial expansion of the experimental tools available to characterize these exotic TRSB superconducting states.
Detection of the generalized clapping modes also offers the unique ability to estimate the relative magnitudes of each order parameter component through the frequency at which the mode resides (as seen in all of our results, the phase mode frequency is maximum for an equal admixture of order parameter components, and decreases as one component or the other becomes dominant). Moreover, observing these modes using ultrafast THz techniques would allow one to assess TRSB superconductivity in driven non-equilibrium superconductors. The exploration and manipulation of generalized clapping modes also represent a distinct direction in the wider, rapidly developing field of THz control of unconventional superconductors and their collective excitations 74,78-83 .
As mentioned above, the generalized clapping mode spectrum can be used to unambiguously distinguish chiral TRSB states from non-chiral states, as a generic TRSB superconductor has two non-degenerate generalized clapping modes, whereas it is only in the special case of a chiral state that the two modes are degenerate.
This capability is particularly well-suited to clarify the structure of the order parameter in Sr 2 RuO 4 -a problem of tremendous current interest. At the time of writing, the two leading candidate order parameters for this system are the mixed symmetry d + ig state (see Figs. 3c, f and 4e, f) and the chiral d + id state (see Fig. 5d), where the generalized clapping mode spectrum is plotted for the cylindrical Fermi surface relevant to Sr 2 RuO 4 ). Quasiparticle interference 84 and thermal transport 85 measurements have demonstrated that the order parameter exhibits vertical line nodes along the zone diagonals (i.e., along the [110] direction), consistent with the d + ig scenario, and counter to the horizontal line nodes expected for a d + id state. In contrast, recent muon spin relaxation measurements failed to detect a splitting of the critical temperature under hydrostatic pressure 86 , which is only consistent with the chiral d + id order. Given these seemingly conflicting results, a measurement of the generalized clapping mode spectra, via e.g., microwave absorption measurements, could provide crucial insight into the chirality, or lack thereof, of the order parameter, and help the community converge on one candidate order parameter over the other.
Beyond its utility as a form of spectroscopy, the detection of generalized clapping modes is also interesting from a fundamental physics perspective, as the analog of the clapping mode in 3 He-A has yet to be realized in any electronic system. Previously, the existence and detection of clapping modes were considered primarily in p + ip superconductors 29,87,88 , which have proven elusive to realize experimentally 89 . This work opens up a number of other materials platforms 12-27 as candidate systems to finally realize these exotic collective modes in a solid-state system.
Finally, we speculate that our work may also be extended to the enigmatic pseudogap phase of the cuprate high-temperature superconductors where a non-zero Kerr rotation has been reported 90 , suggesting the existence of a TRSB phase above the superconducting transition. It has recently been suggested 36,[91][92][93][94] that the relative phase between two order parameter components can acquire a phase stiffness before either order parameter becomes phase coherent and condenses. In this scenario, the relative phase mode studied in this work could persist even above the superconducting transition, representing a distinctive collective excitation in a TRSB metallic phase.

Code availability
All relevant code used in this study is available from the corresponding authors upon reasonable request.