Abstract
Sr2RuO4 is the best candidate for spin-triplet superconductivity, an unusual and elusive superconducting state of fundamental importance. In the last three decades, Sr2RuO4 has been very carefully studied and despite its apparent simplicity when compared with strongly correlated high-T c cuprates, for which the pairing symmetry is understood, there is no scenario that can explain all the major experimental observations, a conundrum that has generated tremendous interest. Here, we present a density-functional-based analysis of magnetic interactions in Sr2RuO4 and discuss the role of magnetic anisotropy in its unconventional superconductivity. Our goal is twofold. First, we access the possibility of the superconducting order parameter rotation in an external magnetic field of 200âOe, and conclude that the spinâorbit interaction in this material is several orders of magnitude too strong to be consistent with this hypothesis. Thus, the observed invariance of the Knight shift across T c has no plausible explanation, and casts doubt on using the Knight shift as an ultimate litmus paper for the pairing symmetry. Second, we propose a quantitative double-exchange-like model for combining itinerant fermions with an anisotropic Heisenberg magnetic Hamiltonian. This model is complementary to the Hubbard-model-based calculations published so far, and forms an alternative framework for exploring superconducting symmetry in Sr2RuO4. As an example, we use this model to analyze the degeneracy between various p-triplet states in the simplest mean-field approximation, and show that it splits into a single and two doublets with the ground state defined by the competition between the âIsingâ and âcompassâ anisotropic terms.
Similar content being viewed by others
Introduction
Superconductivity in Sr2RuO4, even though it occurs at a rather low temperature, has been attracting attention comparable to that attached to high-temperature superconductors.1 For many years the dominant opinion was that it represents a unique example of a chiral triplet pairing state.2,3,4,5 Interestingly, the original premise that led to this hypothesis was the presumed proximity of Sr2RuO4 to ferromagnetism, and thus it was touted as a three-dimensional analogue of 3He.6, 7 It was soon discovered, first theoretically,8 and then experimentally,9 that the leading instability occurs in an antiferromagnetic, not ferromagnetic channel, and thus a spin-fluctuation exchange (FLEX) in the Berk-Schrieffer spirit would normally lead to a d-wave, not p-wave superconductivity.
The issue seems to have been decided conclusively when the Knight shift on Ru was shown to be temperature-independent across T c,2 and later also on O,3 and the neutron-measured spin-susceptibility was found to be roughly constant across the transition as well.4 The chiral p-wave state with an order parameter \({\bf{d}} = {{const}}\left( {x + iy} \right){\bf{\hat z}},\) where the Cooper pair spins can freely rotate in-plane, is the only state that could have this property. Moreover, since in this state spins are confined in the xy plane, the Knight shift in a magnetic field parallel to \({\bf{\hat z}}\) is supposed to drop below T c in pretty much the same manner as in singlet superconductors. Nonetheless, when eventually this experiment was performed,10 it appeared that K z is also independent of temperature. The authors of ref. 10 attempted to reconcile the accepted pairing symmetry with their experiment, by assuming that the experimental magnetic field of 200âOe is affecting drastically the pairing state and converting it to \({\bf{d}} = f(x,y){\bf{\hat y}}\) (or the corresponding x â y partner state). One goal of our paper is to estimate whether this hypothesis is tenable with realistic material parameters.
It is worth noting that the invariance of the in-plane susceptibility is the only experiment consistent exclusively with a chiral p-state (CpS). Some probes indicate chirality (μSR detects spontaneous currents below T c 5), while others indicate breaking of time-reversal symmetry,11 but the triplet parity is not, in principle, necessary to explain these experiments. For instance, the singlet chiral state Îâ=âconst(xzâ+âiyz), or even Îâ=âconst(x 2â+ây 2â+âiαxy), which is not chiral (although this second state would require two-phase transitions with decreasing temperature, which has never been detected), but breaks the time reversal symmetry, are other admissible candidates. Josephson junction experiments12 suggested that the order parameter changes sign under the (x, y) â (âx, ây) transformation, which is consistent with a CpS, but also with other order parameters.13
Spinâorbit (SO) coupling plays an important role not only in selecting between different triplet states (chiral vs. planar), but also in the structure of the chiral state itself. For instance, in Cu x Bi2Se3, instead of the expected chiral state, a nematic spin-triplet state was observed.14, 15 Indeed, in Cu x Bi2Se3 the large SO coupling necessarily implies that a \({\bf{d}} = const(x + iy){\bf{\hat z}}\) induces also an in-plane d-vector component \(const\left({{\hat{\bf x}} + i{\bf{\hat y}}} \right)z\).16 This in-plane component leads to a non-unitary pairing state, which is not energetically favored in weak coupling,17 and, instead, the lower-symmetry nematic state with \({\bf{d}} = {c_1}x{\bf{\hat z}} + {c_2}z{\hat{\bf x}}\) is realized. In principle, similar physics must occur in Sr2RuO4, but there the corresponding induced in-plane d-vector component should be much smaller, thus allowing for a chiral p-wave state to exist. However, this is a quantitative, not qualitative difference, and needs a better understanding of the role of SO coupling.
Finally, recent years have brought about an array of experiments that are actually inconsistent with the CpS. One prediction of a CpS is the existence of edge states at boundaries and at domain walls.18,19,20 However, no evidence for these edge states has been found.19, 21 There are a variety of predictions about the response of CpS to in-plane magnetic fields that have not been observed experimentally. In particular, it is known that a finite in-plane magnetic field should lead to two superconducting transitions as temperature is reduced22, 23 and that the slope of the upper critical field with temperature at T c should depend on the in-plane field direction (this is only true for pairing states that can break time-reversal symmetry).23, 24 In addition, several different probes indicate behavior resembling substantial Pauli paramagnetic effects (see ref. 25 for discussion and original references). The latest cloud on the CpS sky appeared because of the uniaxial strain experiments. For the CpS (or, in fact, any other two-component state) the critical temperature, T c, under an orthorhombic stress must change linearly with the strain (the \(x{\bf{\hat z}}\) and \(y{\bf{\hat z}}\) state will not degenerate any more, and the splitting is linear in strain). In the experiment26 T c varies at least quadratically (more likely, quartically), whereas the linear term is absent within the experimental accuracy, and only one, very well-expressed specific heat jump, ÎC, has been observed, with no trace of a second transition even while the critical temperature changes a lot (Note: Li, Gibbs, Mackenzie, Hicks, and Nicklas (2017). Heat capacity measurements of Sr2RuO4 under uniaxial stress. Unpublished (reported at the APS March Meeting, New Orlean)). Moreover, it was established that both T c and the ÎC variations trace the changes in the density of states, and peak when the Fermi level passes the van Hove singularities at the X or Y points. This observation is particularly important, because, by symmetry, the superconducting gap in a triplet channel in a tetragonal superconductor is identically zero at X and Y (it need not be zero at a finite k z, but in a highly 2D material like Sr2RuO4 it will be still very small by virtue of continuity). Correspondingly, one expects these van Hove singularities to have little effect on superconductivity. A slightly more subtle, but even more convincing argument against triplet pairing in ref. 26 is related to the reduced critical field anisotropy. Finally, a recent detailed study of thermal conductivity has concluded that a d-wave state is by far better consistent with the thermal transport than the CpS.27
In fact, only one fact unambiguously points toward the CpS: the invariance of the spin susceptibility in the in-plane magnetic fieldâbut, as discussed above, the analogous experiment for the out-of-plane field also shows such an invariance. Thus, our acceptance of the nuclear magnetic resonance data as an ultimate proof of the CpS hinges upon the possibility of a magnetic field Bâââ200âOe (0.02âT, or 13âmK in temperature units) to overcome the energy difference between the helical \(\left( {{\bf{d}} \bot {\bf{\hat z}}} \right)\) and chiral \(\left( {\left. {\bf{d}} \right\|{\bf{\hat z}}} \right)\) states. One can show (the derivation is presented below) that this implies that the two states, whose energy difference comes from the SO interaction, are nearly degenerate with the accuracy δâââ10â7âKâââ10â10â λ, where λâââ100âmeV is the SO constant. Moreover, it is often claimed that the solution of other paradoxes outlined above may be obtained (although nobody has convincingly succeeded in that) in a formalism where the relativistic effects would be fully accounted for, since the separation between singlet and triplet channels is only possible in terms of the full angular moment, rather than just electron spins.
Results and discussion
In order to illustrate how SO coupling affects the core assumption of the field-induced d-vector rotation, let us show a simple back-of-the-envelope calculation: suppose that the one-electron Hamiltonian has a relativistic term of the order of \(\kappa M_z^2\). The physical meaning of this term is that in the normal state when n electron spins are confined in the xy plane (as opposed to be parallel to z), this affects the exchange part of the effective crystal potential, and, correspondingly, one-electron energies. The change is proportional to n, and so is the number of affected one-electron states, leading to an energy loss of the order of κn 2, where κ is the magnetic anisotropy scale that is determined by the SO coupling. One way in which this energy contribution manifests itself is the conventional magnetic anisotropy in a spin-ordered state in which case nâââM/μ B . However, the same âfeedbackâ effect must be present in a triplet superconducting state. The number of electrons bound in Cooper pairs and thus forced to be either parallel or perpendicular to z can be estimated as nâ~âÎN, where Î is some average superconducting gap, and N is the density of states, which has been experimentally measured to be about 8âstates/spin/Ru/eV.1 Assuming Î ~7.5âK, we estimate n ~0.005âe/Ru. If the magnetic anisotropy scale κ is of the order of 10âK (we will show later that this is the case), then the total energy loss incurred by rotating the spins of the Cooper pairs is ÎE scâââ2âÃâ10â4âK (this is smaller than various model estimates of the change in T c, as reviewed in ref. 28; we use the above estimate because we wanted to have a conservative lower bound on ÎE sc and a model-independent estimate of the energy, and not simply a critical temperature difference, since the latter may, in principle, dramatically differ from the former). This seems like a small number, but we shall compare it with the energy gained by allowing screening of an external field of 200âOe by Cooper pairs, which is \(\Delta {E_{{\rm{mag}}}} \approx \mu _B^2{B^2}N \approx {10^{ - 7}}\,{\rm{K}}\). This is four orders of magnitude smaller than the estimated loss of superconducting energy. In other words, to allow for the presumed d-vector rotation, various relativistic effects must fortuitously cancel each other with a 10â3 accuracy. Note that in ref. 28, instead, ÎT c was compared with the Zeeman splitting, μ B B, but this comparison is hardly relevant at all for the problem at hands; the correct way is to compare the energy gain with the energy loss.
This simple estimate emphasizes the importance of getting a handle of the type and scale of relativistic effects in Sr2RuO4. So far all efforts in this direction have been performed either within simplified models or by educated guesses from the experiment.29,30,31,32,33 The goal of this paper is to address the issue from a first-principle perspective. It is known that this approach correctly describes (only slightly underestimating) the SO interactions34 (our SO splitting is exactly the same as calculated in that reference, 90âmeV), and, by comparing the Fourier transform of the calculated exchange interaction with the experimentally measured q-dependent spin susceptibility, we observe that the latter is also well reproduced. The only serious problem with this approach is that it overestimates the tendency to magnetic ordering for a given set of magnetic interactions because of the mean field nature of the density functional theory (DFT). Thus, we start with a realistic paramagnetic state of Sr2RuO4, using the alloy analogy model in the first-principles DFT framework and calculate the isotropic exchange interactions (see âMethodsâ). The Fourier transform of these interactions gives us the shape of the full spin susceptibility in the momentum space; as expected, this is peaked at the nesting vector \({{\bf{q}}_3} = \left( {1,1,0} \right)\frac{{2\pi }}{{3a}}\), in agreement with the experiment. Next, we calculate the mean-field energy of several ordered magnetic states, all characterized by the same wave vector qâ =âq 3, and degenerate without SO interaction. This shall allow us to calculate nearest neighbor (NN) relativistic Ising terms (see below). Finally, we calculate magnetic anisotropy for the \({\bf{q}} = \left( {1,0,0} \right)\frac{\pi }{a}\) states, which breaks the tetragonal symmetry, and from there we extract the NN compass exchange (see âMethodsâ). The energy scale of magnetic anisotropy appears rather large, which not only renders the hypothesis of a d-vector rotation unlikely, but also supports the idea that anisotropic interactions must be properly accounted for before drawing conclusions from the experiment. The set of interactions that we derived should serve as a launching pad and testbed for model calculation of the superconducting properties. We maintain that a model where all interelectron interactions are absorbed into spinâspin interactions (with Hundâs coupling between the spins and non-interacting electrons) is complementary to the widely used Hubbard-model and at least as realistic.
Experimentally, Sr2RuO4 shows no sign of magnetic ordering down to the low temperatures. However, neutron diffraction studies have revealed35,36,37 spin-fluctuations in the paramagnetic state with a characteristic nearly commensurate wave vector \({\bf{q}} = \left( {0.3,0.3,0} \right)\frac{{2\pi }}{a}\), close to \({{\bf{q}}_3} = \left( {1,1,0} \right)\frac{{2\pi }}{{3a}}\), which persist even at the room temperature.38 The DFT, being a static mean field theory (by some criteria, the best such theory possible), overestimates the tendency to magnetism. In its generalized gradient approximation flavor DFT stabilizes even ferromagnetic order, albeit with small moments.39 Unsurprisingly, spin-density waves (SDWs) with qâ=âq 3 are even lower in energy. This deficiency of the DFT can, however, be put to a good use by mapping the DFT (i.e., mean field) energetics onto a spin-Hamiltonian, as it is often done, for instance, for Fe-based superconductors.40 Since the isotropic and anisotropic magnetic interactions entails completely different energy scales, and require different level of accuracy, we have chosen two different techniques to calculate them; as discussed below, the isotropic calculations were performed perturbatively, allowing us to fully account for the long-range, nesting-driven interaction, while the NN exchange interactions were calculated by brute force comparing highly accurate energy values in different magnetic configurations.
First, we have calculated the Heisenberg part of the Hamiltonian, defined as:
where M i is the Ru moment on the site i, and the summation is performed over all bonds up to a given coordination sphere. The parameters are calculated in the disordered local moments (DLM) approximation,41 which is used to model the paramagnetic state of Sr2RuO4 (see âMethodsâ for more details and employed approximations).
The results presented in the Fig. 1 are derived for the Ru local moment being fixed to 1âμ B in the DLM state. The obtained values of the exchange constants, however, are fairly independent of the values of the local moment fixed in the DLM state; the minimum of the Fourier transform is always at qâ=â(α, α, 0) with αâ=â0.3âââ0.31. Note that the interplane exchanges nearly vanish, indicating an almost perfect 2D character of the magnetism in Sr2RuO4. For instance, the NNs between-the-planes J 001âââ0.5âââ1âK/\(\mu _B^2\) (ferromagnetic), or about â0.01J 200.
The leading term is the in-plane third NN antiferromagnetic interaction J 200, which is quite counterintuitive from the point of view of the Hubbard-model and superexchange mechanism that is often employed as a starting point. This is a consequence of the Ru electrons itinerancy, since Sr2RuO4 is a metal. The lattice Fourier transform, J(q), of the calculated interactions is shown in the Fig. 1b. J(q) has a meaning of a measure of the energy (J(q)·M 2) of the spin-density fluctuations with a wave-vector q and a given amplitude M [the quantitity that is directly related to the static zero-temperature spin susceptibility is 1/J(q)]. The deep minima of J(q) at \({\bf{q}} = \left( {0.31,0.31,0} \right)\frac{{2\pi }}{a}\) suggest that the spin-fluctuations with the wave vector q will be dominant in the paramagnetic state of Sr2RuO4. The position of these minima is indeed in perfect agreement with the sharp maxima of the integrated magnetic scattering intensity, experimentally observed in neutron diffraction.38 Thus, both our calculation and the experiment suggest the dominance of the spin-fluctuations with the wave vector q 3 in the excitation spectra of Sr2RuO4.
In order to extract the relevant anisotropic exchange interaction parameters, we used direct calculations of the total energy in different magnetic configurations compatible with the ordering vector q 3. Note that anisotropic magnetic interactions appear exclusively due to the SO coupling (see âMethodsâ for the description of codes and approximations used in these calculations). Allowed anisotropic terms for the NN terms are absorbed in the following Hamiltonian (simplified compared to a more complete expression discussed in the âMethodsâ section):
where the first term is given by Eq. (1), the second is Ising exchange (sometimes called the Kitaev interaction), and the last two represent the compass term. Summation in the last two terms is over all horizontal and all vertical bonds, respectively, while in the Ising term it is over all inequivalent bonds. Note that Dzyaloshinskii-Moriya terms42, 43 are not allowed by symmetry.
The six most energetically favorable states are depicted in Fig. 2. The first three states can be described as harmonic SDWs:
where \({{\bf{A}}_a} = \left[ { - \frac{1}{2},\frac{{\sqrt 3 }}{2},0} \right]\), \({{\bf{A}}_b} = \left[ {\frac{i}{{2\sqrt 2 }}, - \frac{i}{{2\sqrt 2 }},\frac{1}{2}} \right]\), \({{\bf{A}}_c} = \left[ {\frac{i}{{2\sqrt 2 }},\frac{i}{{2\sqrt 2 }},\frac{1}{2}} \right]\), with m hardly varying between the three states and equal to 0.76âμ B . The fourth to sixth states are collinear where the amplitude of the moments varies along each of the crystallographic directions 100, 010, and 110 as mâ², âmâ²/2, âmâ²/2 (more precisely, 1.07, â0.56, â0.56âμ B ). Note that mâ² is very close to \(\sqrt 2 m\) in the harmonic SDWs, and the average ãM 2ã is the same in all these states (within a 1.3% error). In this collinear state the direction of the magnetization can be selected in three inequivalent ways, namely, along 110, 1\({{\bar 1}}\)0, or 001. Upon inclusion of the SO term, the 001 collinear up-up-down structure is the ground state (Table 1).
Next, we fit the energy differences in Table 1 to the Hamiltonian (Eq. (2)), extracting J zz and J xy (the fitting procedure included more parameters than in Eq. (2), and is discussed in the âMethodsâ section). All isotropic (Heisenberg) parts of the exchange interactions are included in the H H . The compass parameter J xy does not affect the states with q â (1, 1, 0), and was extracted from a separate set of calculations with \({{\bf{q}}_2} = \left( {0,1,0} \right)\frac{\pi }{a}\) and \({{\bf{M}}_{ijk}} = m{\bf{A}}\,{\rm{exp}}\left( { - i{{\bf{R}}_{ijk}} \cdot {{\bf{q}}_2}} \right)\), where A â¥â=â(1, 0, 0) and A ||â=â(0, 1, 0), and m was fixed to be equal to its value in the spiral states, 0.76âμ B . These wave vectors define so-called single stripe antiferromagnetic order, well known in Fe-based superconductors.
Thus, obtained parameters are J zzâ=ââ1.2â±â0.6âmeV/\(\mu _B^2\) and J xyâ=â1.0âmeV/\(\mu _B^2\) (J zz m 2â=ââ0.70â±â0.35âmeV, J xyâ=â0.57âmeV, for mâ=â0.76âμ B ). The details of the fitting are described in the âMethodsâ section. Note that J xy does not have an error bar not because it was accurately determined, but because we did not have enough calculations to estimate the error. First, one observes that the scale of the anisotropy induced by SO is of the order of 10âK. As discussed in the introduction, this renders the explanation of the invariance of the Knight shift below T c in term of the order parameter rotation10 untenable and shakes the main argument in favor of the chiral triplet superconductivity in Sr2RuO4. Second, our fitting provides a powerful tool for modeling normal and especially superconducting properties of Sr2RuO4 from an entirely different perspective. Compared to the generally accepted models based on the Hubbard-Hund Hamiltonians, our new approach is based entirely on first-principles calculations, and emphasizes the role of magnetic interactions. The corresponding DFT-inspired model Hamiltonian reads:
where the first term is the non-interacting energy, with the band (spin) indices α (s), and the second is the Hundâs rule (Stoner, in the DFT parlance) coupling. All electronâelectron interactions carried by spin fluctuations are absorbed in the local Hundâs interaction I and the intersite magnetic interactions H rH , while interactions due to charge fluctuations are not included in Eq. (5), but can be added separately, if needed (or just collected in one Coulomb pseudopotential μ *, as in the Eliashberg theory). Equation (5) can be understood as a generalized double-exchange Hamiltonian.44 Indeed, this model, inspired by DFT calculations, entails electrons moving in the same effective potential as used in other techniques, and described by the same tight-binding parameters. However, as it is usual in DFT, all electronâelectron interactions are implicitly integrated out. Instead, we introduce quasi-local magnetic moments that interact with the electrons via the local Hundâs rule coupling (parameterized as the Stoner parameter in DFT), while the moments interact among themselves according to the sum of the long-range Heisenberg and the short-range anisotropic Hamiltonian (Eq. (2)). The former part incorporates implicitly all Fermi surface effects, including nesting at \({\bf{q}} = \left\{ {0.3,0.3,0} \right\}\frac{{2\pi }}{a}\), while the latter selects between different triplet states. It is important not to attempt to integrate out the free carriers c k αs in Eq. (5) in order to extract additional interaction between the local moments M; that would have been incorrect, because all such interactions had been computed previously and embedded in H rH . On the contrary, the intended solution of these equations is integrating out the Mâs in order to obtain the effective pairing interaction, as illustrated below.
It might be instructive to demonstrate how Eqs. (4) and (5) can be reduced to a Hamiltonian including only the itinerant electrons (as convenient for analyzing superconductivity). We can safely assume that all Js are much smaller than I, introduce the itinerant spin polarization \({{\bf{s}}_{i\alpha }} = \mathop {\sum}\nolimits_{ss'} {c_{i\alpha s}^\dag {\sigma _{ss'}}{c_{i\alpha s'}}} \), and single out the terms relevant to the pairwise interaction between s iα and s jβ :
In the lowest order in J, the mean field solution requires that M i and s iα be parallel, \({E_{ij,\alpha \beta }} = - 2{\rm{IMs}} - {J_{ij}}{M^2}{{\bf{\hat s}}_{i\alpha }} \cdot {{\bf{\hat s}}_{i\beta }}\), and the effective pairwise interaction can be written as \( - {J_{ij}}{M^2}{{\bf{\hat s}}_{i\alpha }} \cdot {{\bf{\hat s}}_{i\beta }}\) (note that essentially the same Hamiltonian, only written in the orbital basis rather than the band basis, which can also be done in this case, was applied to Fe-based superconductors in several papers, for instance, in ref. 45; after summation of the total energy over the band indices α, β these approaches become equivalent). In principle, one can easily derive the next-order correction to the interaction, which is \( + \left( {J_{ij}^2{M^3}{\rm{/}}Is} \right){\left( {{{{\bf{\hat s}}}_{i\alpha }} \cdot {{{\bf{\hat s}}}_{i\beta }}} \right)^2}.\)
As an example of how this Hamiltonian can be used to address superconductivity, we solve in the simplest mean field approximation the problem of the relative energetics of the five unitary p-triplet states. In particular, beginning with \(H = - {J_{ij}}{M^2}{{\bf{\hat s}}_{i\alpha }} \cdot {{\bf{\hat s}}_{i\beta }}\) and restricting the electronic spins to a single band for simplicity (generalizing onto three bands with realistic dispersions is straighforward), we find that the Ising and compass exchange modify the pairing interaction δV in different pairing channels differently, as shown in Table 2.
Thus, in this approximation the five states split into two planar doublets (of course, this degeneracy is not driven by symmetry, and will be lifted in more sophisticated calculations, but likely the splitting will be small) and a CpS singlet, which is located between the doublets if J zz â>ââ|J xy | and below both of them otherwise (note that we found J zz to be negative). In other words, we have shown that selection between chiral and planar superconductivity is driven by the competition between the Ising and compass anisotropic exchange. Of course, this is just an illustration of principle; in principle, this approach should be applied to the true three-band electronic structure and extended to singlet as well as triplet states, but this is beyond the scope of this paper.
We reiterate that we do not insist that this approach is superior to the Hubbard-Hamiltonian, but it is different and complementary, having the potential to uncover new physics. Similar to the former, it can be used in the contexts of, e.g., random-phase approximation, FLEX, or functional renormalization group calculations.
A final note relates to the recent experiments on strained Sr2RuO4. This is a large topic mostly outside of the scope of this paper. However, we would like to make one comment in this regard. The fact that T c rapidly grows with the strain and peaks at the strain corresponding to the Lifshits transition (where the γ band touches the X-point) can be explained by either a DOS effect (van Hove singularity) or by a change in pairing interaction. The former explanation, as mentioned before, is realistic for singlet, but not triplet pairing symmetries. The latter would be viable if the changes in DOS were sufficient to shift the balance between the AF and FM tendencies toward the latter. To verify that, we have repeated the calculations of the Heisenberg parameters in the strained case. However, we found that the main effect of the strain is not related to the van Hove singularity, and that the average exchange coupling does not become more ferromagnetic. Instead, the strain introduces a splitting between J 1a and J 1b , while the average value barely changes, as shown in Fig. 3. These results therefore indicate that the peak in T c is directly related to the peak in DOS, and not via enhanced pairing interaction. This conclusion is supported by recently reported thermodynamic results (see Note), which strongly suggest that not only T c, but also ÎC/T c is peaked at the van Hove singularity.
To summarize, we have presented first-principles calculations of the leading isotropic and anisotropic magnetic interactions in Sr2RuO4. Our results indicate that rotating a p-wave superconducting order parameter during measurements of the Knight shift is impossible by several orders of magnitude, and thus the invariance of the Knight shift across the transition remains an unresolved puzzle. We further proposed a model framework, based on a double-exchange type Hamiltonian, and incorporating the calculated magnetic interactions in their entirety, and present an example of using this framework for addressing superconducting pairing symmetry.
Methods
First-principles calculations
For relativistic total energy calculations we have employed the projector augmented wave method46 as implemented in the Vienna Ab initio Simulation Package,47 including SO coupling.48 We have used the DFT within the PerdewâBurkeâErnzerhof parametrization for the exchange and correlation potential,49 and the experimental lattice structure is employed in all calculations. The energy cutoff was set to 400âeV with convergence criteria of 10â6âeV. We used up to 1386 irreducible k-points, reduced to 900 for the four formula units cell. For Ru, a pseudopotential with p-states included as valence states was selected.
For the calculation of the isotropic exchange constants we used the KorringaâKohnâRostokker method within the atomic sphere approximation50 and the Green function-based magnetic-force theorem.51 The implementation of this technique has been described elsewhere.52 Physically, this technique can be considered to be a magnetic analogue of the disordered alloys theory based on coherent potential approximation52 and is known as the DLM approximation.41, 53 Upon fixing the Ru magnetic moments in the DLM state we achieved self-consistency using 115 irreducible k-points in the Brillouin zone, and then used an extended set of k-points (1529) to compute the isotropic exchange constants in the framework of the magnetic force theorem.
Fitting procedure
The full equation used to describe the calculated energies, including all bilinear terms up to the second neighbors, reads:
Here, for completeness, we have included the single-site anisotropy term K; since it always enters in the same combination with J zz, they cannot be decoupled within this set of calculations. While this term is, in principle, allowed because of itinerancy, we note that the calculated magnetization is close to the Sâ=â1/2 and therefore we expect \(K \ll {J^{zz}}\). This approximation was used in the main text. We have also included, besides the NN anisotropic interaction \(J_1^{zz}\) and \(J_1^{xy},\) the corresponding second NN interactions \(J_2^{zz}\) and \(J_2^{xy}\). The latter distinguishes between the collinear state polarized along the (110) tetragonal direction and the one polarized along (100), and the transverse and rolling spirals. We found it to be relatively small, 0.17â±â0.05âmeV/\(\mu _B^2\). The second NN Ising interaction \(J_2^{zz}\) simply adds to J zz, and therefore was absorbed into the latter in the fitting procedure. The difference in energies between the planar spiral and the (100) collinear structure, 0.24âmeV, is likely related to the fact that the isotropic exchange constants enter these two state differently. Our non-relativistic calculations find them degenerate within the computational accuracy, apparently, fortuitously. Since SOC also affects the isotropic constants, it is no surprise that relativistic effects break this accidental degeneracy.
One can calculate KâââJ zz and \(J_2^{xy}\) either from the set of spiral calculations, or from collinear calculations; the results differ by ±30%. It is unlikely that this is due to computational inaccuracy, but rather to other interactions not accounted for, such as third neighbors (which is the leading isotropic exchange) or anisotropic biquadratic coupling.
The full summary of the magnetic patterns and their energies used for the fitting, as well as the expressions for the total energies in terms of the parameters in Eq. (8), are presented in Table 3.
Mean-field comparison of pairing energies
To find the interactions in Table 2, we begin with the following Hamiltonian H int that includes charge and spin fluctuations. As an example of how this approach can be used we ask a relatively simple question of how the magnetic anisotropy we have found affects spin-triplet pairing states. To this end, we generalize ref. 54 and consider only a single band with the following Hamiltonian with charge, Ï(q), and spin, S i (q), interactions:
Focussing on superconductivity with zero momentum Cooper pairs, H int can be rewritten as:
where \({s_k} = \mathop {\sum}\nolimits_{s,s'} {{{\left( {i{\sigma _y}} \right)}_{s,s'}}{c_{ - k,s}}{c_{k,s'}}} \) and \({t_{i,k}} = \mathop {\sum}\nolimits_{s,s'} {{{\left( {i{\sigma _i}{\sigma _y}} \right)}_{s,s'}}{c_{ - k,s}}{c_{k,s'}}} \) are the possible singlet and triplet Cooper pair operators, and the effective interactions for the different pairing channels are found to be
This result reduces to that found when spin interactions are isotropic54, 55 or have uniaxial symmetry.56 In our case, the specific form of the spin anisotropy is
Expressing H int with the above spin anisotropy in terms of irreducible representations of tetragonal symmetry for the Cooper pairs leads to Table 2.
Data availability
The authors declare that all source data supporting the findings of this study are available within the paper.
References
Mackenzie, A. P. & Maeno, Y. The superconductivity of Sr2RuO4 and the physics of spin-triplet pairing. Rev. Mod. Phys. 75, 657â712 (2003).
Ishida, K. et al. Anisotropic pairing in superconducting Sr2RuO4: Ru NMR and NQR studies. Phys. Rev. B 56, R505âR508 (1997).
Ishida, K. et al. Spin-triplet superconductivity in Sr2RuO4 identified by 17O Knight shift. Nature 396, 658â660 (1998).
Duffy, J. A. et al. Polarized-neutron scattering study of the cooper-pair moment in Sr2RuO4. Phys. Rev. Lett. 85, 5412â5415 (2000).
Luke, G. M. et al. Time-reversal symmetry-breaking superconductivity in Sr2RuO4. Nature 394, 558â561 (1998).
Rice, T. M. & Sigrist, M. Sr2RuO4: an electronic analogue of 3He? J. Phys. Condens. Matter 7, L643âL648 (1995).
Mazin, I. I. & Singh, D. J. Ferromagnetic spin fluctuation induced superconductivity in Sr2RuO4. Phys. Rev. Lett. 79, 733â736 (1997).
Mazin, I. I. & Singh, D. J. Competitions in layered ruthenates: ferromagnetism versus antiferromagnetism and triplet versus singlet pairing. Phys. Rev. Lett. 82, 4324â4327 (1999).
Braden, M. et al. Inelastic neutron scattering study of magnetic excitations in Sr2RuO4. Phys. Rev. B 66, 064522 (2002).
Murakawa, H., Ishida, K., Kitagawa, K., Mao, Z. Q. & Maeno, Y. Measurement of the 101Ru-Knight shift of superconducting Sr2RuO4 in a parallel magnetic field. Phys. Rev. Lett. 93, 167004 (2004).
Xia, J., Maeno, Y., Beyersdorf, P. T., Fejer, M. M. & Kapitulnik, A. High resolution polar Kerr effect measurements of Sr2RuO4: evidence for broken time-reversal symmetry in the superconducting state. Phys. Rev. Lett. 97, 167002 (2006).
Nelson, K. D., Mao, Z. Q., Maeno, Y. & Liu, Y. Odd-parity superconductivity in Sr2RuO4. Science 306, 1151â1154 (2004).
ŽutiÄ, I. & Mazin, I. Phase-sensitive tests of the pairing state symmetry in Sr2RuO4. Phys. Rev. Lett. 95, 217004 (2005).
Matano, K., Kriener, M., Segawa, K., Ando, Y. & Zheng, G.-Q. Spin-rotation symmetry breaking in the superconducting state of Cu x Bi2Se3. Nat. Phys. 12, 852â854 (2016).
Yonezawa, S. et al. Thermodynamic evidence for nematic superconductivity in Cu x Bi2Se3. Nat. Phys. 13, 123â126 (2017).
Yip, S.-K. Models of superconducting Cu:Bi2Se3: single- versus two-band description. Phys. Rev. B 87, 104505 (2013).
Fu, L. Odd-parity topological superconductor with nematic order: application to Cu x Bi2Se3. Phys. Rev. B 90, 100509 (2014). (R).
Matsumoto, M. & Sigrist, M. Quasiparticle states near the surface and the domain wall in a p x ± ip y -wave superconductor. J. Phys. Soc. Jpn. 68, 994â1007 (1999).
Kallin, C. & Berlinsky, A. J. Is Sr2RuO4 a chiral p-wave superconductor? J. Phys. Condens. Matter 21, 164210 (2009).
Scaffidi, T. & Simon, S. H. Large Chern number and edge currents in Sr2RuO4. Phys. Rev. Lett. 115, 087003 (2015).
Kirtley, J. R. et al. Upper limit on spontaneous supercurrents in Sr2RuO4. Phys. Rev. B 76, 014526 (2007).
Agterberg, D. F. Vortex lattice structures of Sr2RuO4. Phys. Rev. Lett. 80, 5184â5187 (1998).
Mineev, V. P. Superconducting phase transition of Sr2RuO4 in a magnetic field. Phys. Rev. B 89, 134519 (2014).
Gorâkov, L. P. Anisotropy of the upper critical field in exotic superconductors. JETP Lett. 40, 1155â1158 (1984).
Amano, Y., Ishihara, M., Ichioka, M., Nakai, N. & Machida, K. Pauli paramagnetic effects on mixed-state properties in a strongly anisotropic superconductor: application to Sr2RuO4. Phys. Rev. B 91, 144513 (2015).
Hicks, C. W. et al. Strong increase of T c of Sr2RuO4 under both tensile and compressive strain. Science 344, 283â285 (2014).
Hassinger, E. et al. Vertical line nodes in the superconducting gap structure of Sr2RuO4. Phys. Rev. X 7, 011032 (2017).
Maeno, Y., Kittaka, S., Nomura, T., Yonezawa, S. & Ishida, K. Evaluation of spin-triplet superconductivity in Sr2RuO4. J. Phys. Soc. Jpn. 81, 011009 (2012).
Ng, K. K. & Sigrist, M. The role of spin-orbit coupling for the superconducting state in Sr2RuO4. Europhys. Lett. 49, 473â479 (2000).
Eremin, I., Manske, D. & Bennemann, K. H. Electronic theory for the normal-state spin dynamics in Sr2RuO4: anisotropy due to spin-orbit coupling. Phys. Rev. B 65, 220502 (2002). (R).
Annett, J. F., Györffy, B. L., Litak, G. & WysokiÅski, K. I. Magnetic field induced rotation of the d-vector in the spin-triplet superconductor Sr2RuO4. Phys. Rev. B 78, 054511 (2008).
Cobo, S., Ahn, F., Eremin, I. & Akbari, A. Anisotropic spin fluctuations in Sr2RuO4: role of spin-orbit coupling and induced strain. Phys. Rev. B 94, 224507 (2016).
Scaffidi, T., Romers, J. S. & Simon, S. H. Pairing symmetry and dominant band in Sr2RuO4. Phys. Rev. B 89, 220510 (2014).
Veenstra, C. N. et al. Spin-orbital entanglement and the breakdown of singlets and triplets in Sr2RuO4 revealed by spin-and angle-resolved photoemission spectroscopy. Phys. Rev. Lett. 112, 127002 (2014).
Sidis, Y. et al. Evidence for incommensurate spin fluctuations in Sr2RuO4. Phys. Rev. Lett. 83, 3320â3323 (1999).
Servant, F. et al. Magnetic excitations in the normal and superconducting states of Sr2RuO4. Phys. Rev. B 65, 184511 (2002).
Braden, M. et al. Anisotropy of the incommensurate fluctuations in Sr2RuO4: a study with polarized neutrons. Phys. Rev. Lett. 92, 097402 (2004).
Iida, K. et al. Inelastic neutron scattering study of the magnetic fluctuations in Sr2RuO4. Phys. Rev. B 84, 060402 (2011). (R).
de Boer, P. K. & de Groot, R. A. Electronic structure of magnetic Sr2RuO4. Phys. Rev. B 59, 9894â9897 (1999).
Glasbrenner, J. K. et al. Effect of magnetic frustration on nematicity and superconductivity in iron chalcogenides. Nat. Phys. 11, 953â958 (2015).
Gyorffy, B. L., Pindor, A. J., Staunton, J., Stocks, G. M. & Winter, H. A first-principles theory of ferromagnetic phase transitions in metals. J. Phys. F Met. Phys. 15, 1337â1386 (1985).
Dzyaloshinskii, I. E. A thermodynamic theory of weak ferromagnetism of antiferromagnetics. J. Phys. Chem. Solids 4, 241â255 (1958).
Moriya, T. Anisotropic superexchange interaction and weak ferromagnetism. Phys. Rev. 120, 91â98 (1960).
Khomskii, D. I. Basic Aspects of the Quantum Theory of Solids (Cambridge University Press, 2010).
Seo, K., Bernevig, B. A. & Hu, J. P. Pairing symmetry in a two-orbital exchange coupling model of oxypnictides. Phys. Rev. Lett. 101, 206404 (2008).
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953â17979 (1994).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169â11186 (1996).
Steiner, S., Khmelevskyi, S., Marsmann, M. & Kresse, G. Calculation of the magnetic anisotropy with projected-augmented-wave methodology and the case study of disordered Fe1âx Co x alloys. Phys. Rev. B 93, 224425 (2016).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865â3868 (1996).
Ruban, A. V. & Skriver, H. L. Calculated surface segregation in transition metal alloys. Comp. Mater. Sci. 15, 143â199 (1999).
Liechtenstein, A. I., Katsnelson, M. I., Antropov, V. P. & Gubanov, V. A. Local spin density functional approach to the theory of exchange interactions in ferromagnetic metals and alloys. J. Magn. Magn. Mater. 67, 65â74 (1987).
Ruban, A. V., Shallcross, S., Simak, S. I. & Skriver, H. L. Atomic and magnetic configurational energetics by the generalized perturbation method. Phys. Rev. B 70, 125115 (2004).
Cyrot, M. Phase transition in Hubbard model. Phys. Rev. Lett. 25, 871â874 (1970).
Miyake, K., Schmitt-Rink, S., & Varma, C. M. Spin-fluctuation-mediated even-parity pairing in heavy-fermion superconductors. Phys. Rev. B 34, 6554 (1986).
Cho, W., Thomale, R., Raghu, S. & Kivelson, S. A. Band structure effects on the superconductivity in Hubbard models. Phys. Rev. B 88, 064505 (2013).
Kuwabara, T., & Ogata, M. Spin-Triplet Superconductivity due to Antiferromagnetic Spin-Fluctuation in Sr2RuO4. Phys. Rev. Lett. 85, 4586â4589 (2000).
Acknowledgements
B.K., S.K., and C.F. were supported by the joint FWF and Indian Department of Science and Technology (DST) project INDOX (I1490-N19), and by the FWF-SFB ViCoM (Grant No. F41). I.I.M. is supported by ONR through the NRL Basic Research Program. D.F.A. was supported by the National Science Foundation Grant DMREF-1335215.
Author information
Authors and Affiliations
Contributions
S.K. and I.M. conceived the research; B.K. has carried out most of the numerical calculations with contributions by S.K. and I.M.; the superconductivity-related discussion was authored by D.A., who also performed the sample mean-field calculations, and by I.M. All authors participated in the discussions and contributed to writing the paper; C.F. supervised the Vienna part of the project.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Additional information
Publisherâs note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kim, B., Khmelevskyi, S., Mazin, I.I. et al. Anisotropy of magnetic interactions and symmetry of the order parameter in unconventional superconductor Sr2RuO4 . npj Quant Mater 2, 37 (2017). https://doi.org/10.1038/s41535-017-0041-8
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41535-017-0041-8
This article is cited by
-
Unconventional superconductivity in Cr-based compound Pr3Cr10âxN11
npj Quantum Materials (2024)
-
Cooper Pairing in A Doped 2D Antiferromagnet with Spin-Orbit Coupling
Scientific Reports (2018)