A proposal for reconciling diverse experiments on the superconducting state in Sr2RuO4

A variety of precise experiments have been carried out to establish the character of the superconducting state in Sr2RuO4. Many of these appear to imply contradictory conclusions concerning the symmetries of this state. Here, we propose that these results can be reconciled if we assume that there is a near-degeneracy between a d_{x^2-y^2} (B_{1g} in group theory nomenclature) and a g_{xy(x^2-y^2)} (A_{2g}) superconducting state. From a weak-coupling perspective, such an accidental degeneracy can occur at a point at which a balance between the on-site and nearest-neighbor repulsions triggers a d-wave to g-wave transition.


I. INTRODUCTION
Sr 2 RuO 4 should be the poster child for the theory of "unconventional" superconductors: it is a clean, stoichiometric crystalline material, it settles into a Fermi liquid "normal state" in a range of temperatures, T , that extends roughly a factor of 20 above the superconducting T c , the normal state fermiology is well established, the superconducting transition is sharp and of mean field type, and there is no doubt that the superconducting state is unconventional in the sense of a sign-changing order parameter with a mean equal to or near 0 along the Fermi surface [1][2][3][4]. Of all currently known unconventional superconductors, Sr 2 RuO 4 is the one in which a BCS-like weak coupling mean-field-theory treatment is most easily justified, meaning that the microscopic "mechanism" of superconductivity should be well defined. Instead, a large number of careful experiments have lead to a set of observations which cannot be accounted for even at the scenario level by any current theory -35 years after the first phenomenological theories [5,6] we are still debating the symmetry of the superconducting order parameter [4].
Indeed, some sets of observations almost appear to be mutually contradictory: The evidence [7][8][9] of time-reversal symmetry breaking (TRSB) in the superconducting (SC) state is difficult to reconcile with the existence of gapless nodal lines in the quasiparticle spectrum [10][11][12][13][14][15]. The evidence of a two component superconducting order parameter [16,17] is difficult to reconcile with the lack of singular T c dependence on symmetry breaking strain [18] or in-plane magnetic field [19]. The evidence of triplet (odd parity [20]) pairing is difficult to reconcile with the evidence of singlet pairing from NMR [21,22] and from the fact that the in-plane critical field appears to be Pauli-limited [23].
It is always possible that some set of experimental observations are misleading -perhaps they reflect an extrinsic effect or have an alternative interpretation that has not been considered. Here, we will take various observations at face value, and attempt to construct a phenomenological scenario that can reconcile as many of them as possible. To do so, we have introduced a single "fine -tuned" assumption -that Sr 2 RuO 4 happens to be close to a multi-critical point at which two symmetry-distinct pairing symmetries are degeneratein particular, a d x 2 −y 2 (B 1g in group theory nomenclature) and a g xy(x 2 −y 2 ) (A 2g ) state. As we will discuss later in the paper, even with this assumption, there are some observations that either require a more complete theory to interpret (i.e. the interpretation does not follow simply from general symmetry considerations) or an alternative interpretation of the experiment.
The article is organized as follows: we first outline the primary experimental facts we have in mind. Those can be considered the axioms of our phenomenological theory. In Section III we set forward the Landau-Ginzburg theory describing a near-degeneracy between a d-wave and a g-wave SC state, and discuss some features of the resulting phase diagram. In Section IV we briefly address more features of the phases and phase transitions at a level that would be expected from a BCS-type treatment of the problem; while these considerations could straightforwardly be extended to include the effects of a realistic treatment of the multiband character of the SRO fermiology, here we offer a stripped-down discussion that focusses on qualitative aspects of the problem. In particular, we explore the specific heat signatures of the various transitions expected in the theory. In Section VII, we discuss some of the experimental observations that are not immediately accounted for by our scenario. This particularly concerns reconciling our proposal with the observed Kerr signal, which requires including considerations of explicit breaking of spatial symmetries and effects of disorder.
In both Sections II and VII, we allude to proposed future work that could more seriously take account of specific microscopic aspects of SRO and thus could lead to further avenues for comparing theory with experiment. In Section VIII, we conclude by discussing some possible experiments that could critically test the applicability of our proposed scenario.

II. EXPERIMENTAL "FACTS" IN SRO
There have been a remarkably large number of extremely careful experiments carried out to explore the character of the superconducting (SC) state in Sr 2 RuO 4 (SRO). We will focus on a subset of these -that we will accept as experimental facts with clear implications.
• 1) There is a single superconducting (SC) transition with a well-defined T c , and the transition is largely mean-field like. This inference follows from the fact that the transition as seen in transport and in magnetization is sharp, and the specific heat shows a sharp mean-field-like anomaly at T c . There is no indication of hysteresis (as could happen if the SC transition were first order), significant regimes of fluctuation SC above T c , or any glassy feature of the SC transition [4].
• 2) The superconducting phase is "unconventional," that is the SC order parameter changes sign in such a way that when averaged over the Fermi surface it is approximately zero. This inference follows [3] from the extreme sensitivity of T c to nonmagnetic disorder, i.e. T c → 0 when the elastic mean-free path, , approaches the clean limit SC coherence length, ξ 0 . (The existence of nodal quasiparticles, discussed below, constitute additional evidence of unconventional pairing.) • 3) To a high level of precision, there are gapless nodal quasiparticles along line nodes that are more or less perpendicular to the Ru-O plane, i.e. they extend along the tetragonal c-axis. This follows from the low temperature behavior of the specific heat and from the low temperature thermal conductivity [13]. Recent STM studies further corroborate this conclusion [15]. While the location (in k-space) of the nodes is not entirely clear, both the thermal conductivity experiments and the STM results are strongly suggestive that the line nodes lie along wave-vectors that are at 45 o to the Ru-O bond direction, that is along the tetragonal (1,1,0) and (1,-1,0) directions. If there is any gap at the putative nodal points, it has been estimated [12,13] that it must be less than 10µeV, or less than 3% of the total gap.
• 4) The SC T c ( ) is a smooth and non-singular function of symmetry breaking (shear) strain, , near = 0. In the case in which a stress is applied along (1,1,0), T c is a weak linear function of the stress; since by symmetry T c must be an even function of a purely shear strain, the linear stress dependence presumably reflects the fact that, in addition to a xy (B 2g ) component, the uniaxial stress produces a significant symmetric component of symmetric strain (with A 1g symmetry) 0 . There is a nearly symmetric variation of T c in response to uniaxial stress along (1,0,0). This suggests that, in this case, the induced symmetric strain is small, and that what is being observed is primarily the response to a B 1g shear strain, x 2 −y 2 . Local measurements [24] (using scanning magnetometry) of T c vs stress show a quadratic minimum at each location, where, however, T c (0) and¯ vary from location to location in the crystal; the variation of¯ is, presumably, associated with a random quenched distribution of local strains. Indeed, an interpretation in terms of random strain is strongly indicated based on the fact that some regions of the crystal are seen to have local T c 's that are of order 10% higher than the zero-strain bulk T c ; given the sensitivity of the SC state to impurity scattering [3], if what was being seen were regions with higher and lower than average impurity concentration, there would be a sharply defined maximum T c (associated with all sufficiently clean portions of the crystal) and then a distribution of more disordered regions with reduced local T c 's. [51] • 5) Time reversal symmetry is spontaneously broken below a temperature T trsb that is (in unstrained crystals) approximately equal to the SC T c . This inference follows directly from the observation of an anomalous Kerr signal [8], and from µSr experiments [7]; it is also consistent with experiments on the geometry dependence of the Josephson relation [9], although beyond supporting the time-reversal symmetry breaking, this last set of experiments is still not fully digested.
• 5b) On symmetry grounds, equilibrium currents are possible in a TRSB SC state when not forbidden on the basis of purely spatial symmetries. This is particularly true for the case of chiral superconductors that break both time-reversal and inversion symmetry. Theoretical estimates of how large these expected currents should be vary greatly depending on details of the assumed microscopic structure of the state, but initially suggested that they should be detectable by high sensitivity magnetic scanning microscopy. However, careful searches [25] for the magnetic fields produced by such currents have failed to detect any such currents, despite sensitivity that is at least two orders of magnitude greater than what would be required to detect the currents implied by the most straightforward theoretical analysis (based on an assumed chiral p x + ip y state).
• 6) There is a sharp drop in the Knight-shift upon cooling through T c , which continues to the lowest temperatures and fields that have been probed [21,22]. While it is still necessary to pursue this experiment to lower fields (relative to H c2 ) and lower T (relative to T c ), already the existing observations are strongly suggestive of a spinsinglet order parameter. The NMR spectrum shows no evidence of a phase transition (within the SC state) as a function of strain; given this, results at high strain (where both H c2 and T c are greatly enhanced) can be included in this analysis, making the evidence for singlet pairing considerably stronger.
• 7) There is a discontinuity at T c in the elastic modulus c 66 associated with the shear B 2g ( xy ) strain, as inferred from ultrasound experiments [16,17]. While the elastic moduli associated with A 1g compressional strains exhibit discontinuities at all second-order phase transitions, such a discontinuity in a shear elastic constant is most straightforwardly associated with a two-component SC order parameter [26]. This follows from the requirement that a gauge-invariant bilinear can be constructed from the order parameter that transforms as xy under the point group operations -such a bilinear cannot be constructed for any one-component SC order parameter, but can be constructed for a number of two-component order parameters. It is also important to note that no discontinuity has been observed in the elastic modulus associated with B 1g ( x 2 −y 2 ) strain.
However, the interpretation of this split transition is made somewhat subtle by the observation that in high resolution specific heat measurements, while the critical anomaly at T c ( ) is still clear and relatively sharp and of the expected magnitude, no thermodynamic signature of any feature associated with a transition at T trsb ( ) has been detected [28].
Despite a long history of interpreting results in SRO in terms of a p-wave order parameter (or more technically, an odd parity order parameter) and of the existence of a variety of other experimental results that provide circumstantial support to this idea, we will interpret the new NMR results as ruling this out [29]. Focusing on possible singlet (even parity) SC states, consistent with the tetragonal symmetry of the crystal there is one possible symmetry-related two-component state, which corresponds to {d xz , d yz } (E g ) pairing, and four possible single We discard the symmetry-related two-component option (E g ) [30][31][32] from consideration for two reasons: Firstly, the thermal conductivity and STM indicate vertical line nodes while these states have symmetry protected horizontal line nodes. This is not entirely definitive, since such a state could have accidental vertical line nodes as well, but this seems a rather unnatural occurrence. Secondly, regarding the non-onsite pairing E g proposal [30], such a state would have a vanishing pair-wave function for two electrons in the same Ru-O plane, which seems implausible on theoretical grounds, given the extreme quasi-2D nature of the Fermi surfaces. Regarding the local pairing proposal which involves the multi-orbital nature of the ruthenates [31,32], it would produce a jump in the elastic modulus for both the x 2 −y 2 and xy type strain, whereas only the latter shows a discontinuity in experiments. Moreover, there is currently no clearly established signatures of the predicted BdG Fermi surface that arises in such a state.
Arguments discriminating among the remaining options are slightly more subtle. The d x 2 −y 2 and g xy(x 2 −y 2 ) have symmetry protected vertical line nodes along the (1,1,0) direction, while for any of the other cases, such nodes would necessarily be accidental nodes. It is not yet clear how precisely the observed line nodes are aligned with this symmetry axis -the more precise the experimental constraints on this, the more solid is the identification of these as the relevant states.
Having focused on two distinct 1D representations, the natural next step would be to determine which one is favored. However, several of the observations summarized above require a two-component order parameter -specifically 5, 7 and 8. We are thus driven to assume that there is an accidental near-degeneracy between d x 2 −y 2 and g xy(x 2 −y 2 ) pairing, i.e. we assume that SRO happens to be fine-tuned close to the boundary between a regime in which one of these is dominant, and the other in which the other dominates.

III. LANDAU GINZBURG TREATMENT
Thus, in the following, we will assume that at zero applied strain, SRO happens to be tuned to a multi-critical point at which these two SC states have (nearly) the same T c .
Depending on other details, the SC phase below T c can have various patterns of broken symmetry. This is most easily explored by considering the Landau-Ginzburg effective free energy density expressed in terms of the complex (charge 2e) scalar fields, D and G, that represent the local value of the pair-fields of the stated symmetries: where V 0 includes the terms in the effective potential (through quartic order in the fields) in the absence of disorder, external strain (i.e. it assumes tetragonal symmetry), and under the assumption of an exact Z 2 symmetry between the D and G orders, while V 1 includes to quadratic (leading) order in the fields terms related to breaking of this accidental Z 2 symmetry, as well as the effects of disorder and externally applied strain (breaking of tetragonal symmetry), K is a quadratic form in the gauge invariant gradients of the order parameter, and . . . refers to higher order terms in powers of the fields and their derivatives. Specifically, and Note that V 0 has an enlarged SO(4) symmetry for the special case in which γ 1 = γ 2 = 0.
Here, the principal T dependence is incorporated in a T dependence of α 0 , such that it changes from positive to negative as T drops from above to below T c . The multicritical point that serves as the focus point for our analysis arises when h 1 = 0; otherwise, the mean value of h 1 encodes the preference for one or the other form of pairing, with d-wave pairing favored for h 1 > 0. Below T c , the nature of the SC state is determined by the signs and relative magnitudes of the quartic terms. To simplify the discussion, we will consider the case in which γ 0 is positive and larger in magnitude than the remaining terms. In the uniform SC state when h 1 = 0, the preferred form of ordering is determined by the relative sign and magnitude of γ 1 and γ 2 . For γ 1 < 0 and γ 2 > γ 1 the SC state is either pure d or pure g, and thus preserves all spatial symmetries as well as time-reversal symmetry. The case that will be of primary interest here is that in which γ 1 > 0 and γ 2 > 0, so the SC state breaks both time-reversal symmetry as well as some lattice symmetries; we refer to this as a TRSB SC, and it corresponds to d x 2 −y 2 ± ig xy(x 2 −y 2 ) pairing. Finally, if γ 1 > γ 2 and γ 2 < 0 the SC state breaks various lattice symmetries but preserves time-reversal symmetry; this is a nematic superconductor, which has two possible symmetry related states d x 2 −y 2 ± g xy(x 2 −y 2 ) , depending on the relative sign of the two components.
Turning to the terms that appear in V 1 , h 0 has the interpretation as a local T c , h 1 as the local deviation from multicritcality (i.e. a local preference for d or g), and h 2 has the symmetry of a random local strain, xy . Importantly, such a local strain favors one orientation of the nematic order over the other. In the absence of an externally applied strain, the configuration (or spatial) average of all the random fields can be assumed to vanish, h j ( r) = 0, but in the presence of applied strain, where a sum over a is implicit, Λ a j and Λ ab j are appropriate coupling constants, the indices a and b (which are implicitly summed over) label the components of the strain tensor -which we decompose according to symmetry as a = s for A 1g (although in a tetragonal system, there are properly two distinct components of A 1g strain), a = (x 2 − y 2 ) for B 1g , and a = xy for B 2g , and . . . indicates higher order terms in powers of the strain. All components of Λ a j = 0 by symmetry other than Λ s 0 , Λ s 1 , and Λ xy 2 . Similarly, all terms in Λ ab j = 0 other than the diagonal terms, Λ a,a 0 and Λ a,a 1 (for any a), and the off-diagonal term, Λ s,xy 2 = Λ xy,s 2 . In other words, there generically is a shift of the mean T c and of the splitting between the d and g wave order that is linear in the isotropic strain and quadratic in each component of the shear strain. However, the degeneracy between d x 2 −y 2 + g xy(x 2 −y 2 ) and d is lifted only in the presence of non-zero B 2g strain. The detailed shape and even topology of the resulting phase diagrams in the T -strain plane depend on the relative magnitudes of various parameters in the Landau-Ginzburg free energy; representative phase diagrams as a function of strain (assuming degeneracy between d and g at zero strain) are shown in Fig.   1 When we come to discuss disorder, we will define the root mean squared disorder strengths, σ j , as (where ν is the volume of a unit cell). Because the most relevant disorder is likely that associated with inhomogeneous strain, it may be important to consider the range of the disorder correlations, as well.

IV. QUASIPARTICLE PROPERTIES; BCS CONSIDERATIONS
While we are free to normalize the order parameter fields that appear in the Landau-Ginzburg effective field theory in any way we like, the same freedom does not pertain when In keeping with the present scenario, we assume that the D and G wave orders are essentially degenerate at zero strain. All panels were simulated with LG parameters γ 0 = 2, γ 1 = 1, γ 2 = 0.5. we relate these fields to the gap parameter that governs the quasi-particle spectrum. Interpreting the various possible states through the lense of the BdG equations of BCS mean-field theory, we relate the gap parameter to the order parameter fields according to where θ is the overall SC phase, φ is the relative phase between the d and g-wave components, and where we normalize the real dimensionless gap functions F a ( k) so that Max[F a ( k)] = 1; thus D and G are, respectively, the maximum gap in the d-wave and g-wave state, respectively. F a can have an arbitrarily complicated k dependence (depending on microscopic details), so long as they transform in the desired fashion under the crystalline symmetries.
In keeping with the spirit of the Landau-Ginzburg approach, at least in a range of T near to T c , we can imagine that all the singular T and strain dependences in the problem are inherited from those of the order parameter fields, D and G. From here, as usual, it follows (still assuming the validity of BCS mean-field theory) that the quasiparticle spectrum is where ε( k) is the normal state (Fermi liquid) dispersion.
At the level of BCS mean-field theory, once the values of D(T ) and G(T ) are determined, a host of dynamic and thermodynamic quantities can be computed in the usual way. One particularly important quantity is the jump in the specific heat at the various transitions.
Consider, for example, the case in which h 1 > 0 and γ 1 , γ 2 > 0, so that upon lowering the temperature there is first a transition at T c to a pure d-wave SC state, and then at a lower temperature, T trsb < T c , a transition to a d x 2 −y 2 + ig xy(x 2 −y 2 ) state. If we assume the simplest possible mean-field dependence of the gaps the specific heat jump at T c can be expressed as and the specific heat jump at T trsb as where f and [[f ]] signify, respectively, the average of f over the Fermi surface and over the Brillouin zone, E d ( k) is the dispersion in the pure d-wave SC state, and F is a dimensionless function defined implicitly in Eq. 11.
In Fig. 2 we have computed F as a function of T trsb /T c for a band-structure corresponding to nearest-neighbor hopping on a square lattice with a mean electron density n = 1.38, per site and assumed F a have the simplest forms consistent with symmetry, , and a typical mean-field value of D 0 /T c = 2.5. For comparison, we have performed the same calculation foo the case in which the second order parameter is a simple s-wave, F g → F s = 1.
It is striking that the specific heat jump is much smaller for the transition to the d+ig state, and especially so the smaller T trsb /T c . There are several contributing factors to this. Firstly, because the d-wave gap grows so singularly with decreasing temperature below T c , the total remaining entropy in the system is rapidly reduced relative to that in the normal statethis effect would apply for any "second" transition in the SC state. However, the effect is particularly dramatic due to the fact that the g-wave gap also has nodes in the one region of the Fermi surface that remains ungapped (nodal) in the d-wave SC state. Thus, the density of thermally excited quasiparticles is only slightly (further) reduced upon entering the d+ig -TRSB phase.

V. A PLAUSIBLE SCENARIO
The set of experimental observations that we have enumerated above have, for some time, seemed difficult to reconcile with each other. With the assumption that SRO is tuned close to the conjectured multi-critical point, much of this seeming difficulty is readily eliminated. For h 1 = 0, there is a single SC transition, and the SC state is unconventional, in keeping with points 1 and 2 . The SC order is even parity, spin-singlet pairing, consistent with point 6. Because both the d and the g wave components of the order parameter are odd under any spatial transformation that exchanges k x and k y , consistent with point 3, the existence of nodal lines in the quasi-particle spectrum along the (1,1,0) direction is protected by symmetry, even in the TRSB SC phase; this is distinct from the case for any other order parameter symmetries that have been considered to date. Importantly, at the assumed multicritical point, the d and g-wave components of the order parameters behave like a two-component order parameter with precisely the requisite relative symmetries to account for point 7. Specifically, the bilinear formed of B 1g and A 2g order parameters is of the B 2g representation, thus producing a discontinuity in c 66 at T c . This scenario also provides a natural explanation for why no discontinuity is seen in (c 11 − c 12 )/2-the B 1g elastic modulus-since such a bilinear cannot be formed. This is in contrast with other two-component superconducting order parameters, such as d xz,yz or p x,y states, where a discontinuity in (c 11 − c 12 )/2 is expected in general, and the lack of experimental evidence for such a discontinuity would have to be explained away in some other way.
Since T c appears to be most strongly sensitive to B 1g strain, it is important to stress that no singular non-analytic changes in either T c or T trsb as a function of The TRSB in the SC state discussed as point 5 is consistent with the present scenario so long as we assume all γ j > 0. Precisely at the point at which the two orders are degenerate, the onset of TRSB would occur simultaneously with T c , consistent with present observations. However, as is exhibited in the various schematic phase diagrams in Fig. 1, any strain is expected to drive the system away from this fine-tuned condition, resulting in a phase diagram in which T c > T trsb , consistent with point 8. The fact that one set of the g-wave nodes coincides with the d-wave nodes means that the specific heat jump at the transition from a pure d-wave SC to a d x 2 −y 2 + ig xy(x 2 −y 2 ) state is significantly reduced, which may partially account for the fact that no signature of T trsb has yet been observed in specific heat measurements. (Note, in a regime in which the higher T transition is to a pure g-wave SC, one would expect a stronger thermodynamic signature of the lower temperature transition to the d x 2 −y 2 + ig xy(x 2 −y 2 ) state since one set of g-wave nodes would be gapped below T trsb .
There is good reason to expect small or even vanishing equilibrium currents in this state.
Since the SC order of a d x 2 −y 2 + ig xy(x 2 −y 2 ) state preserves parity, it is non-chiral, i.e..the associated Bogoliubov bands have zero Chern number. Thus there is no topological reason to expect edge currents. Moreover, the fact that reflection through a mirror plane perpendicular to the (110) direction is respected even in the d x 2 −y 2 +ig xy(x 2 −y 2 ) state is sufficient to establish that edge currents must vanish at any surface that respects this symmetry.
However, despite these successes, one is hard-pressed in the context of the present scenario to simultaneously account for a portion of the observations in point 4 (no non-analytic xy strain dependence of T c in the limit of vanishing stain). To account for this in a consistent fashion we invoke the effect of weak disorder that is equivalent to random local strains. Note that while it is an established fact that the relevant materials are, to an extraordinary level of approximation, nearly perfectly crystalline, with normal state mean-free paths in excess of 1 µm (point 2), long-wave-length disorder (i.e. strain disorder) does not make much contribution to determining the mean free path, and is not expected to lead to a substantial reduction of T c , even in an unconventional SC. Indeed, as discussed in point 4, there is some direct experimental evidence of the existence of the presence of random strain [24].
There are several different scenarios that can (at least roughly) reconcile the remaining observations by incorporate the effects of weak, random strain -that is to say non-zero σ j as defined in Eq. 5. We plan to address this issue in a forthcoming work. In general, for the same reasons discussed in Ref. 33

VI. ASPECTS OF A MICROSCOPIC PICTURE
Because the normal state of SRO is a good Fermi liquid, and because the SC phase transition appears to be mean-field like, it makes sense to think about the mechanism of SC from a weak-coupling BCS perspective. Various studies using an asymptotically exact weakcoupling approach to the problem taking into account only on-site interactions have found that the SRO band-structure is conducive to SC states with p-wave symmetry and d x 2 −y 2wave symmetry [34][35][36]. Similar conclusions have been drawn on the basis of a random phase approximation ansatz [37] and multi-orbital functional renormalization group (FRG) studies which better account for the strong antiferromagnetic fluctuations seen in neutron scattering experiments on SRO [38,39]. In all cases, the pair-wave-function approximately vanishes on site, thus avoiding the effects of a repulsive on-site Hubbard U . Given the uncertainty concerning the correct relative magnitudes of the various interactions, and the intrinsic issues concerning the applicability of such approaches to real materials (where interactions are never weak), even thought the p-wave instability is often found to be dominant, it would be difficult to argue strongly on purely theoretical grounds which of these is preferred.
However, as far as we know, no previous study of SRO has found indication of a significant tendency to pairing in a g-wave channel.
However, in a related study, the effect of a repulsive nearest-neighbor repulsion, V , on the character of the weak coupling SC for the Hubbard model on a square lattice was analyzed [40]. In the range of doping in which the dominant pairing instability is in the d (x 2 −y 2 ) channel at V = 0, there is a transition to a g-wave state for V greater than a critical value, V c ∼ U 2 ρ(E F ). This is natural in the sense that the pair wave-function in the g-wave state not only vanishes on-site, but also on any pair of nearest-neighbor sites.
Thus, for V = V c , precisely the sort of accidental degeneracy we have invoked arises. (We are currently in the process of extending this analysis to the case of the SRO bandstructure [41].) Alternatively, even though probably more contrived than the latter microscopic reasoning, a proximity of d (x 2 −y 2 ) -wave and g-wave pairing propensity had been found in an effective twodimensional electronic model inspired by the ruthenates which was assumed to be subject to three-dimensional phonon coupling [42].
Another important insight from the weak coupling approach concerns the form of the gap function. All such calculations [34,35,38,[43][44][45] indicate a tendency toward pairwave functions with substantially more complex structure than the minimal functions we considered above. In fact, deep gap minima not required by symmetry arise ubiquitously in such studies. Such "accidental near-nodes" were invoked in earlier attempts to reconcile p x + ip y and d+is SC states with experiments that suggest the existence of line nodes. The possibility of d x 2 −y 2 + is pairing near a point of accidental degeneracy of the two ordering tendencies has been considered previously in the context of iron pnictides [46], and more recently in the context of SRO [47]. While in our present scenario, there are symmetry protected true line-nodes even in the d x 2 −y 2 + ig xy(x 2 −y 2 ) SC phase that would be absent in its d x 2 −y 2 + is-wave cousin, it is possible that further structure of the gaps could be invoked to produce near-nodes at the requisite locations. ( [48], of odd-parity pairing from phase sensitive Josephson tunneling [20], of horizontal line nodes from field dependent specific heat measurements [49], and of half quantum vortices from cantilever magnetometry [50]. All these observatioons pose some challenge to be reconciled with our ansatz, and thus deserve further analysis and revisiting from both an experimental and a theoretical viewpoint.

VIII. PROPOSALS FOR FUTURE EXPERIMENTS
Because the proposed degeneracy between d and g is accidental, it is not plausible that it would turn out to be exact under conditions of zero strain and zero disorder. Thus, we would expect there to be a small but non-zero splitting between T c and T trsb -however, if this is small, it might be difficult to detect under circumstances in which the root mean square variance of h 0 or h 1 (σ 0 or σ 1 ) produce variations of T c and T trsb that are larger than the intrinsic splitting. A possibly more directly testable aspect of this proposal, however, is that the application of any strain -regardless of its symmetry -would be expected to enhance this splitting. Specifically, in the absence of disorder, and under the assumption that the degeneracy is effective exact at zero strain, we would expect a strain-induced splitting In particular, the prediction that a symmetry preserving strain, s , will lead to a splitting of the transition directly reflects the assumption that it is an accidental degeneracy -rather than a symmetry protected one -that is responsible for the two-component character of the order.
Similarly, disorder (especially short-range correlated disorder) is expected to have a substantial impact on T c for any unconventional SC. In the present case, because they are not symmetry related, it should have a different effect on the d and g wave components, thus leading to a split transition. Since the g-wave order has more (symmetry related) nodes on the Fermi surface, one would probably expect disorder to suppress the g-wave component more strongly, thus leading to a case in which there is pure d-wave pairing below T c and a lower temperature transition to a d x 2 −y 2 + ig xy(x 2 −y 2 ) SC at T trsb < T c .
[51] Note that for large strain, x 2 −y 2 , there is a large enhancement of T c (from 1.5K to 3.5K) near a critical strain, which is thought to be associated with a Lifshitz transition in the underlying band structure. The behavior at these larger values of the strain will not be part of the present discussion.
[52] The lack of an observable signature of T trsb in the specific heat is a quantitative issue that will need further investigation. However, as we have shown, the fact that the nodal structure of the d x 2 −y 2 + ig xy(x 2 −y 2 ) SC is the same as for the pure d-wave SC leads to a significant reduction in the expected size of the specific heat jump at T trsb . Moreover, effects of random strain can be expected to broaden and further reduce the magnitude of any specific heat signature of T trsb .
[53] The scanning SQUID measurements of Ref. 24 establish that even locally, the dependence of T c on x 2 −y 2 strain is quadratic, such experiments have not been carried out for xy strain.