{'fraction_non_alphanumeric': 0.03678374392308953, 'fraction_numerical': 0.015076200360517834, 'mean_word_length': 4.629662658753344, 'pattern_counts': {'":': 0, '<': 16, '<?xml version=': 0, '>': 13, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 16, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Thermal, electric and spin transport in superconductor/ferromagnetic-insulator structures

A ferromagnetic insulator (FI) attached to a conventional superconductor (S) changes drastically the properties of the latter. Specifically, the exchange field at the FI/S interface leads to a splitting of the superconducting density of states. If S is a superconducting film, thinner than the superconducting coherence length, the modification of the density of states occurs over the whole sample. The co-existence of the exchange splitting and superconducting correlations in S/FI structures leads to striking transport phenomena that are of interest for applications in thermoelectricity, superconducting spintronics and radiation sensors. Here we review the most recent progress in understanding the transport properties of FI/S structures by presenting a complete theoretical framework based on the quasiclassical kinetic equations. We discuss the coupling between the electronic degrees of freedom, charge, spin and energy, under non-equilibrium conditions and its manifestation in thermoelectricity and spin-dependent transport.


Introduction
The study of the interaction between ferromagnetism and conventional superconductivity in hybrid systems has attracted a great attention during the past decades [1,2,3]. In principle, these two states of matter are antagonistic: whereas the ferromagnetic correlations try to align the spins of electrons, a conventional singlet Bardeen-Cooper-Schrieffer (BCS) [4] superconductor exhibits a condensate with electron pairs (Cooper pairs) with opposite spins. This simple argument explains why usual ferromagnetic metals do not show transition to a conventional superconducting state at any temperature.
Superconductivity and magnetism can coexist though in hybrid systems consisting of superconductor (S) and ferromagnet (F) layers, where the two states couple through mutual proximity effects. The interplay between them leads to striking novel phenomena not present in either system alone. For example the spatial oscillation of the superconducting correlations induced in F via the proximity effect leads to the π-junctions in S/F/S Josephson systems [5,6]. The exchange field in F also leads to a triplet component of the superconducting condensate [1,2,7].
The effects discussed above are caused by the leakage of the superconducting condensate into the ferromagnet. There is also a reciprocal effect, caused by the interaction of the condensate with the exchange field that can leak back into the S region over distances of the order of the superconducting coherence length ξ s . These triplet correlations induced in the superconductor lead to a finite magnetization and a drastic change of the local density of states [8,9,10,11]. If the superconductor is a thin film with thickness smaller than ξ s , this magnetic proximity effect may cause an almost homogeneous spin splitting of the BCS density of states (DOS) [12]. Such spin splitting plays a central role in this review.
In an all-metallic S/F system however it is difficult to achieve such a homogeneous spin splitting. The reason for this is the leakage of Cooper pairs from S to F. Usual metallic ferromagnets such as Fe, Co or Ni, have intrinsic exchange fields much larger than the characteristic superconducting energy (∼ ∆, the BCS pair potential). For thin S-films such strong exchange fields in F lead eventually to a full suppression of the superconducting state. Indeed, the first observation of the 0−π-transition in S/F/S Josephson could be achieved thanks to the use of ferromagnetic alloys with small exchange fields [6,13]. However, a sharp splitting of the DOS in a superconductor in all-metallic systems has not been observed.
A quite different situation emerges when instead of a metallic ferromagnet one uses a ferromagnetic insulator (FI) adjacent to a thin S layer. The band gap of the FI prevents the penetration of the superconducting condensate into the FI. Electrons are reflected at the S/FI interface. Because the FI is magnetic, this reflection is spin-dependent and leads to the creation of an effective exchange field in the superconductor [14]. Superconductivity with a spin-splitting field has been observed in several spectroscopy experiments on Al films attached to europium chalcogenides, such as EuO, EuS and EuSe [15,16,12].
Interestingly, the spin splitting in the DOS observed in those experiments resembles the spin splitting created by a strong in-plane field applied to a thin superconducting film. This is caused by the fact that when a magnetic field is applied parallel to a thin S film, the orbital diamagnetic effect is weak [17] and the field mainly induces a paramagnetic response [18,19,20,21]. This effect was explored experimentally already in the 1970s [22,23]. The spin-split DOS has been utilized to determine the spin polarization of ferromagnets in S/F structures [24,25,26,27,12].
Even though spectral properties of S/FI have been intensively studied for a long time, only recently it has been realized that spin-split superconductors may find applications in spintronics, thermoelectricity, and sensors [3,28,29,30,31,32,33,34]. To understand the physics underlying these applications one needs to consider nonequilibrium phenomena. Nonequilibrium effects can extend over distances larger than the coherence length and they are related to the deviation of the electron distribution function from its equilibrium form. Such a deviation leads to an imbalance of the electron degrees of freedom: charge, energy and spin. From recent studies it has turned out that non-equilibrium effects in spin-split superconductors is a very rich field of research [35,36,37,38,31,39].
In this work we review the non-equilibrium properties of spin-split superconductors, focusing on thermal, electric and spin transport. We have recently published a Colloquium on this topic [3] in which we gave a qualitative overview of the field and the experimental activity. Our main purpose here is different: We mainly focus on the theoretical aspects of the transport properties of S/FI structures and elaborate deeper on the thermoelectricity in superconducting systems. We provide all the necessary theoretical tools for those readers interested in exploring non-equilibrium effects in novel setups combining superconductors and magnetic materials. Those readers seeking more phenomenological approach are referred to the Colloquium mentioned above.
The review is organized as follows: In Sec. 2 we summarize the main properties of S/FI structures and discuss the coexistence between superconductivity and ferromagnetism in these systems. The main focus is on the spin-splitting induced in the S layer and its manifestation on the spectrum. We also analyze the magnetic susceptibility of superconductors and briefly discuss other effects related to spin-splitting fields, including the Fulde-Ferrel-Larkin-Ovchnikov (FFLO) state, triplet superconducting correlations and the cryptoferromagnetic state in the FI layer.
In Sec. 3 we present the quasiclassical Green's function formalism, with a special focus on the description of non-equilibrium properties of diffusive S/FI systems, the role of inelastic scattering and transport through hybrid interfaces. In this section we introduce the concept of nonequilibrium modes which is essential for understanding the different ways of exciting a spin-split superconductor out of equilibrium. Especially the charge, spin and heat transport properties are determined by these modes and their coupling.
On the basis of the quasiclassical formalism we analyze in Sec. 4 the transport properties of different FI/S structures. We focus in particular on spin injection, diffusion and relaxation, in different setups.
In Sec. 4.3.2 we study the non-equilibrium quasiparticle dynamics of a spinsplit superconductor in an alternating radio frequency field, both in the linear and the non-linear regime. For this we make use of time-dependent quasiclassical equations.
Spin splitting the BCS spectrum leads to a strongly electron-hole asymmetric spin-resolved density of states. The symmetry is recovered by averaging over spin, but the electron-hole asymmetry shows up in the case of spin-polarized tunneling into spin-split superconductors. Section 5 discusses the transport properties of such spin-polarized contacts to spin-split superconductors, with a special emphasis on the thermoelectric response of such systems. We also review thermoelectricity in superconducting systems in general to set this finding into historical context. Finally, Sec. 6 summarizes the main phenomena discussed in the review.

Ferromagnetic insulator-superconductor structures
At interfaces between a superconductor (S) and a magnetic material (F) both states of matter may coexist, however with certain modifications. Such modifications extend over the characteristic correlation lengths in each material. In S this length is the coherence length, which in diffusive systems and at low temperatures is of the order of D/∆ 0 , where D is the diffusion coefficient and ∆ 0 the superconducting gap at zero temperature. If the ferromagnet is a metal, the superconducting condensate (Cooper pairs) can penetrate into it. The superconducting correlation length in usual ferromagnets as Ni, Co or Fe, is approximately given by D/h, where h is the intrinsic exchange field of the ferromagnet. Usually this length is very short (a few nanometers or even less) and therefore one speaks about a short-range proximity effect. Interestingly, the interaction of singlet electron pairs stemming from the S layer with the local exchange field leads to pairs in a triplet state that can get back into superconductor [2]. Such leakage of pairs and the suppression of the density of Cooper pairs in singlet state leads to a suppression of superconductivity [1].
If instead of a metal, a ferromagnetic insulator (FI) is placed adjacent to a superconductor, no leakage of pairs is possible. Conduction electrons are reflected at the FI/S interface and during this process their spin interacts via exchange interaction with the localized magnetic moment of the FI material. If the latter is not too strong and the superconducting layer is thin enough, both magnetic and superconducting orders can coexist. In this case the superconductor behaves as a ferromagnetic superconductor with a density of states showing a spin splitting, similar to the one that an external in-plane magnetic field would induce in a thin S layer [12].
In this section we describe the main features of superconductors with a spinsplit density of states due to the proximity of a FI. Such artificial ferromagnetic superconductors are the main focus of this review. They have been fabricated and measured in several experiments along the past decades [15,16,40,41].
2.1. Induced spin splitting in a superconductor-ferromagnetic insulator structure As mentioned above, a way of creating a superconductor with a spin-split density of states is via the exchange field induced by an adjacent ferromagnetic insulator. This may lead to a spin splitting of the order of the superconducting gap, even without an externally applied magnetic field. The first evidence of an exchange field induced in a superconductor via the magnetic proximity effect was observed in a thin Al film in contact with an EuO film [15] and with EuS [16,42,40]. The equivalent internal field of the Zeeman splitting in Al-EuS systems has been reported to be as large as 5 T [43]. Another S-FI combination explored so far is GdN-NbN [44,45]. This system shows, however, a weaker and less clear spin splitting than the Al-europium chalcogenide combination, due to the influence of the sizable spin-orbit interaction in Nb compared to the Al-based devices.
Because the electrons of the superconductor cannot penetrate into the insulating FI, the superconducting properties are only modified by the induced spin-splitting field at the S/FI interface, and not by the leakage of Cooper pairs into the FI. Moreover, FIs can also be used as spin-filter barriers [43], in some cases with a very high spin-filtering efficiency, and therefore they will play a crucial role for several of the applications discussed in subsequent sections which require strong spin filtering.
The amplitude of the spin splitting in a FI/S structure depends on both the intrinsic properties of the superconductor, such as the amount of magnetic impurities and the strength of spin-orbit coupling, and on the quality of the S/FI interface. The latter is crucial for obtaining a large splitting, as shown for example in Ref. [16].
The spin splitting caused by a ferromagnetic insulator on an adjacent superconductor can be understood with the help of the following model [48,49]. The effective Hamiltonian describing a FI/S bilayer consists of the ferromagnetic coupling within the FI layer, conductor via exchange interaction, Here, the symbolΣ means that we consider only the magnetic moments S r localized at the interface. J ex is an effective parameter which for example describes the s-d or s-f exchange interaction. We consider a ferromagnetic insulator with a Curie temperature much larger than the superconducting transition temperature. Thus we can assume that the magnetization of the FI is only determined by the Hamiltonian in Eq. (1) and not affected by the superconducting state [50,51]. With this assumption Hamiltonian (2) describes conduction electrons interacting with an effective exchange field proportional to the local spin average S r at the FI/S interface. This average can be computed by solving independently the magnetic Hamiltonian (1) [41].
In the superconducting state the localized exchange field leads to a modified density of states characterized by the splitting of the coherent BCS peaks. This modification of the spectral properties of the superconductor is non-local and survives over distances away from the FI/S interface of the order of the coherence length ξ s [14,8]. If the thickness d of the S film is much smaller than ξ s , the spin splitting can be assumed to be homogeneous across the film. Thus the density of states is given by Eq. (4) with an effective exchange field h ef f ≈ J ex S r a/d [52,48,14], where a is the characteristic distance between the localized spins.
In this review we mainly focus on thin S films and hence we adopt this approximation for the description of the uniform magnetic proximity effect when the film is adjacent to a FI layer. For an inhomogeneous magnetic configuration of the FI layer the spatial scale that determines the effective splitting is the superconducting coherence length. Thus the splitting becomes observable if the ferromagnet consists of magnetic domains with sizes larger than the superconducting coherence length. This is for example the case of EuS films with magnetic domains of micrometer size [53] that explain the spin splitting observed in experiments [16] on Al-EuS structures, where the Al layers are a few nanometers thick and with a coherence length of the order of 100 nm. In contrast, a spatially fast changing magnetization averages to zero and results in a vanishing effective splitting [41,54].
It also follows from this description that the strength of the spin splitting depends crucially on the quality of the FI/S interface. The growth of any nonmagnetic oxide between the FI and S layers drastically suppresses the effective exchange interaction J ex and hence reduces the spin splitting [16]. In addition to the quality of the interfaces, the spin splitting also depends on the intrinsic properties of the superconducting film. For example, magnetic disorder may lead to a strong suppression of superconductivity in the S film [55]. But even in the absence of magnetic disorder and magnetic impurities, spin-orbit coupling may lead to a modifications in the DOS. In particular, large spin-orbit scattering rates lead to a rounding of the BCS peaks and to a less sharp spin splitting. This explains why splitting has been observed in materials with relative small atomic number Z, such as Al, Be, V, but not in heavier materials such as Pb.
For a detailed description of the intrinsic properties of different superconductors in this respect, we refer the reader to the excellent review in Ref. [12].
Spin splitting in a superconductor can also be achieved by applying an external magnetic field [22,24,26,27,16,40]. In this case, however, the magnetic field also couples to the orbital motion of the electrons and creates circulating currents, the Meissner effect, [56] that try to expel the field from the bulk (diamagnetic response). By increasing the amplitude of the applied field superconductivity is gradually reduced. We denote below this mechanism of suppression of superconductivity by the orbital depairing effect. This mechanism dominates in bulk samples or in thin films with the field applied perpendicular to the plane of the film. If the amplitude of the applied field reaches a critical value H c , the created currents increase the free energy such that the system undergoes a transition into the normal state [17]. However, in superconducting films with thickness smaller than the London penetration length, the critical value H c for a magnetic field applied in-plane largely exceeds the critical value H c⊥ of a perpendicular field [17,57]. In such a case the spin paramagnetic effect (due to Zeeman effect) dominates with respect to the orbital one. The magnetic field penetrates uniformly the film, screening currents are relatively small and therefore H c is limited by the spin paramagnetic effect that tries to align the spin of the original singlet Cooper pairs, as demonstrated in Refs. [19,18]. At T = 0, and in the absence of spin-orbit coupling and magnetic impurities, the critical field due to the paramagnetic effect is given by where ∆ 0 is the superconducting gap at zero field and temperature and µ B is the Bohr magneton. For conventional BCS superconductors with critical temperature ranging from 1 to 10 K, H cP can be of the order of several Tesla and hence superconductivity exists within a large range of field strengths. In short, the paramagnetic effect in thin superconducting films leads to a Zeeman splitting of the density of states (DOS) [22,24,26,27,16,40], similar to the spin-splitting induced in S-FI structures.

Spin-split density of states in a superconductor
The splitting discussed in the previous section can be detected by measuring the tunneling conductance of the superconductor with a normal or a superconducting probe. For example, if one uses a normal metal (N) tunnel coupled to a superconducting film S, the current through the normal-insulatorsuperconductor (NIS) junction is determined by the tunneling expression [58] where N and N s are the reduced density of states of the normal metal and the superconductor, respectively, n F the Fermi distribution function and G T is the conductance of the insulating tunneling barrier I. We first assume G T to have no spin polarization. Electronic transport occurs mainly at energies close to the Fermi level where the density of states of the normal probe N (ε) can be accurately approximated by a function independent of the energy. In such a case, and at low temperatures the measured differential conductance G = dI/dV is proportional to the density of the states of the superconductor. If the normal metal probe is non-magnetic, i.e., electrons are not spin-polarized, the measured G shows the spin-split peaks and is symmetric with respect to the voltage as sketched in Fig. 1(a). If the tunneling barrier has a spin-dependent transmission, i.e., a spin filter, or a ferromagnetic electrode with a non-vanishing spin polarization, the observed peaks are asymmetric [26] [see Fig. 1 (b)]. The asymmetry is proportional to the spin polarization of the conduction electrons of the electrode or the spinfilter efficiency of the barrier. Indeed spin-split superconductor-ferromagnet bilayers have been used for determining the spin polarization of magnetic metals [26,27,24].
The normalized total density of states of the spin-split superconductor, sketched in Fig. 1 (a), can be written as the sum of the DOS N ↑/↓ of each spin species, 1 where ±h eff is the effective spin-splitting field. Equation (4) for the DOS is a simplified description of spin-split superconductors, which is accurate, as we see below, only in certain limiting cases. The expression does not take into account the effect of magnetic impurities or spin-orbit coupling (SOC). For example, magnetic impurities suppress superconductivity and eventually lead to a gap-less situation [55]. On the other hand the SOC counteracts the magnetic impurities but it may also lead to a broadening of the coherence peaks in the spectrum [59,60,61,23]. Moreover, if the spin-splitting field h eff in Eq. (4) is due to an external magnetic field, also orbital depairing contributes to the pair breaking. We present a more quantitative analysis of these effects in the next section.

Paramagnetic depairing mechanisms and the FFLO state
When writing the DOS of a superconductor in an exchange field, Eq. (4), we assume that the ground state of the S layer corresponds to a spatially homogeneous order parameter ∆. The value of the order parameter is determined from the self-consistency equation, written in terms of the Matsubara frequencies ω n = (2n + 1)πT , 2 and a free-energy functional discussed in Sec. 3.5. Here g 01 is a singlet anomalous part of the Green's function described in more detail in Sec. 3. In the following we describe the behavior of ∆ as a function of the exchange field h. For this, we denote ∆ 0 = 1.76T c0 the order parameter at T = h = 0 without depairing effects. T c0 is the corresponding critical temperature. In the absence of spin relaxation the results coincide with those of Ref. [62]: For large values of h, the paramagnetic depairing drives the superconductor to the normal state. There is a critical value of temperature, T 0 ≈ 0.556T c0 , above which the transition between the normal and superconducting states is of the second order. For T < T 0 , the transition is of the first order, and the self-consistency equation has three solutions, of which one is unstable. The stable solutions correspond to superconducting (∆ = 0) and normal (∆ = 0) states. Below the critical field h c , known as the Chandrasekhar-Clogston limit [18,19], the superconducting state is preferred, and above it the normal state. At T = 0 this critical field is h c = ∆ 0 / √ 2. Finally, for h > ∆ 0 , the superconducting solution does not exist at all.
A non-zero spin relaxation rate leads to a quantitative modification of this behavior [63]. In Fig. 2 we show the computed ∆ as a function of temperature and exchange field for a normal-state spin relaxation rate 1/T c0 τ sn = 0.96. Panels (a) and (b) in Fig. 2 show the effects of spin-flip (β = 1) and spinorbit (β = −1) scattering, respectively. The phase-transition line, ∆(T, h) = 0, is shown in Fig. 2 by the solid curve in the (T, h) plane. The red part on this curve corresponds to the second-order transition where ∆(h) continuously goes to zero with increasing h. This behavior changes to the abrupt first-order transition at T < T 0 . The first-order transition line is shown by circles in the (T, h) plane and it is different from the ∆(T, h) = 0 curve shown by the green line at T < T 0 . Both relaxation mechanisms reduce the range of temperatures for which the first-order transition takes place. In other words, the threshold The transition line ∆(h, T ) = 0 is shown by the solid curve in the (h, T ) plane. Its thick red part at T > T 0 shows the second-order transition line. Its thin green part at T < T 0 shows the points on T, h plane where the unstable branch starts (not shown). The first-order transition is shown by the circles in the (h, T ) plane.
temperature T 0 between first-and second-order phase transitions is reduced as compared to the case without relaxation.
In other respects the modification of the superconducting state strongly depends on the spin relaxation mechanism. Spin-flip scattering breaks the timereversal symmetry and therefore leads to a strong suppression of the superconducting gap and T c . On the other hand, spin-orbit interaction is time-reversal invariant and keeps the T c intact.
A striking effect of the spin-orbit scattering is that it increases the critical field of the Chandrasekhar-Clogston limit [63]. This tendency is shown in Fig. 3, where we compare the dependencies of the superconducting gap on h and relaxation rates 1/τ sn for (a) spin-flip and (b) spin-orbit relaxation, respectively. These two cases are characterized by the opposite behaviors of the critical field h c as a function of 1/τ sn -it is suppressed by spin-flip scattering and enhanced by the spin-orbit one. This effect can be understood from the comparison to the changes in the density of states caused by the two spin-relaxation mechanisms (see Fig. 4): spin-flip scattering primarily lifts the gap in the density of states, whereas spin-orbit scattering acts to nullify the spin splitting without affecting the gap. In both cases there is a threshold value of the scattering rate τ −1 0 when the superconducting phase transition changes from the second order at to the first order at τ −1 sn < τ −1 0 . The assumption of a homogeneous order parameter does, however, not always correspond to the lowest-energy state of the system. In Refs. [20,21] it was predicted that the exchange or Zeeman field can cause a transition to an inhomogeneous superconducting state, with a spatial periodical modulation of the order parameter at the scale of the coherence length. Such an FFLO state exists for temperatures below T 0 and only for exchange fields that satisfy 0.71∆ 0 < h < 0.755∆ 0 when T → 0.
The FFLO state has not been observed in conventional superconductors.

13
One of the main reasons for this is that it is very sensitive to disorder [64] and spin-orbit coupling and hence it is expected to occur only in extremely clean samples. In the present review we mainly focus on diffusive systems and therefore we do not pay further attention to the FFLO state. Moreover, the effective exchange fields we mostly consider are far below the range for the FFLO state to occur.
The FFLO state has been also studied in the context of heavy-fermion superconductors [65] and ultra-cold Fermi gases [66,67,68] especially in the context of the BCS-BEC crossover [69,70].

Static spin susceptibility of conventional superconductors
The spin splitting has a direct consequence on the paramagnetic spin susceptibility of the superconductor. In general, the spin susceptibility χ can be found by calculating the spin polarization S = −χh generated by a small spinsplitting field h. In the superconducting state the susceptibility normalized with respect to that of the normal state χ n can be expressed in terms of the difference between spin-up and spin-down DOS N − = N ↑ − N ↓ and the equilibrium distribution function f eq = tanh(ε/2T ) [71,2], Using Eq. (4) for the spin-split DOS in the absence of magnetic depairing processes one can see that ∞ −∞ N − dε = −4h. Therefore for T → 0 and in the absence of spin-dependent scattering χ(T = 0) = 0. This is a consequence of the lack of polarizability of the condensate consisting of spin-singlet Cooper pairs.
However, if spin-dependent scattering is present, in the form of spin-orbit and/or magnetic impurities causing spin flips, the superconducting condensate may exhibit a nonzero paramagnetic susceptibility for T → 0 [71,72].
A generalization of this result can be obtained by using the microscopic equations introduced in Sec. 3 below. Here we are interested in the final results for the susceptibility shown in Fig. 5. In particular, we show the low-temperature dependence of the spin susceptibility as a function of the normal-state spin scattering rate 1/τ sn = 1/τ so +1/τ sf and the parameter β = (τ so −τ sf )/(τ so +τ sf ) determining the relative strength of spin-flip and spin-orbit scattering. From these plots one can see that both scattering mechanisms result in a nonzero susceptibility at low temperatures. Notice that only the spin-flip scattering generates a strong suppression of the order parameter (see Fig. 2). Therefore, at T T c0 the growth of χ as a function of the rate 1/τ sn towards the normal-state value is much faster when the spin-flip scattering dominates the spin-relaxation in the normal state, i.e. when β > 0. The non-zero susceptibility explains the observed Knight shifts in ordinary superconductors [73].
The static spin susceptibility (6) characterizes the paramagnetic response of the superconductor to the external magnetic field. This expression cannot be used for calculating the non-equilibrium spin accumulation induced by the   spin-dependent chemical potential shift ±δµ when the distribution functions in different spin subbands are given by f ↑ (E) = f 0 (E+δµ) and f ↓ (E) = f 0 (E−δµ). The theory of non-equilibrium spin states in superconductors is explained in detail in Secs. 3 and 4. In particular, there we explicitly show [see the discussion after Eqs. (25,26)] that the non-equilibrium spin accumulation generated by δµ is always exponentially suppressed at low temperatures provided that there is a non-zero gap in the quasiparticle spectrum. This is the case even in the presence of the strong spin-orbit interaction [cf. Fig. 4(b)] although the static spin susceptibility (6) is non-zero at T = 0.

Other effects related to the Zeeman and exchange fields
In the next sections we discuss effects related to non-equilibrium properties of spin-split superconductors in FI-S structures. It is however important to mention other effects that we are not considering but that may appear in the presence of Zeeman and exchange fields.

Fermi liquid effects
The analysis in the previous subsections assumes noninteracting electrons, which is a rather good approximation in metals. Corrections due to interactions can be incorporated by renormalizing certain parameters of the theory with the help of the powerful Landau's Fermi-liquid theory [74]. In the context of the quasiclassical theory such effects have been studied in Ref. [75]. In particular, and importantly for our work, the spin-splitting field amplitude h ≡ |h| entering Eq. (8) is renormalized. According to the Fermi liquid theory it is given by [76,40]: where H ext is the external magnetic field, H F I is the field induced by a FI, and G 0 is the effective antisymmetric Fermi-liquid parameter of Landau's theory. The value of G 0 in the superconducting state depends on temperature, is negligibly small at low temperatures (T T c ), and nonzero for temperatures close to the superconducting critical temperature T c and above it. In what follows we neglect the Fermi-liquid corrections by assuming that G 0 1. This is justified at the low temperatures where majority of our work concentrates. However, when describing equilibrium properties of spin-split superconductors at temperatures close to T c , the Fermi-liquid corrections to the field h in Eq. (8) should also be included for some superconductors. For example, for Al films G 0 ∼ 0.16 − 0.3 [77,76] close to T c and above, whereas in V it is negligibly small at all temperatures [78].

Long-range triplet correlations
In this review we mainly focus on the transport properties of the superconductor in S/F structures. For completeness, however, we briefly discuss in this section the appearance of triplet correlations in ferromagnetic metals induced by the superconducting proximity effect [79,2]. Such correlations may induce non-dissipative spin-polarized currents that could be particularly useful to lower the energy consumption in spintronic devices [80,81,7].
Superconducting correlations can be induced in non-superconducting (N) metals by means of the proximity effect [82]. This can be seen as the "leakage" of Cooper pairs from the superconductor into the normal region. The microscopic origin of this "leakage" is the Andreev reflection [83], where an incident electron from N is retro-reflected as a hole, resulting in the formation of a Cooper pair in S. An important consequence of the Andreev reflection is that the correlation between the electron and the hole leads to a formation of a non-vanishing pair amplitude in N, within a characteristic distance ξ N away from the S-N interface. This distance depends on the intrinsic properties of N, such as the degree of disorder, spin-dependent fields, etc. For example, in a diffusive normal metal without magnetic impurities this length is of the order of the thermal length D/T , where D is the diffusion coefficient. In a superconductor-ferromagnet (S-F) junction the situation is very different. In usual ferromagnets the intrinsic exchange field h is much larger than the superconducting characteristic energies. The proximity-induced superconducting correlations decay and oscillate on the length scale ξ F = Re D/h. The oscillation of the condensate leads to the well-established 0-π transition in S-F-S Josephson junctions predicted in Ref. [5] and observed much later in experiments [6,84] on SFS structures with ferromagnetic alloys with a weak exchange field. In conventional ferromagnets, such as Fe, Ni or Co, the effective exchange field is large (∼ 0.1 − 1 eV) and hence ξ F is very short [85]. This explains why the Josephson effect in S-F-S junctions with transition metals has been only observed for thickness of the F layer smaller than 10 nm [86,87].
The situation is very different if the exchange field is spatially non-homogeneous, as for example in a magnetic domain wall. In such a case, as shown in Ref. [79], triplet components with equal spin components can generated. These are insensitive to the local exchange field and hence can penetrate the F metal over the thermal length. Such long-range triplet component of the superconducting condensate explains the long-range Josephson currents observed in S-F-S junctions made for example of half-metal [88,89], ferromagnetic multilayers [90], and ferromagnets with intrinsic inhomogeneous magnetization [91].
We emphasize that the physics of the triplet component addressed in this section describes the equilibrium state of ferromagnets in contact with superconducting electrodes. The main focus of the present review is however the non-equilibrium properties of superconductors in spin-splitting fields, either generated by an adjacent FI or by an external field. As discussed in Sec. 2.3, in order to preserve the superconducting state, those fields have to be smaller than the superconducting gap, and hence all components of the condensate, singlet and triplets, vary over similar length scales. This in particular means that no distinction between short-and long-range has to be made.

Cryptoferromagnetic state
The way we model the exchange field generated in a thin superconducting film by an adjacent ferromagnetic-insulating layer, see Sec. 2.1, is based on the assumption that the magnetic ordering of the FI is the same in both the normal and superconducting states of the S layer. This requires that the effective exchange constants J and J ex in Eqs. (1) and (2), respectively, satisfy J J ex . In this case the ferromagnetic order is robust against the exchange interaction with the conduction electrons and the superconducting condensate. Roughly speaking, this corresponds to the case when the Curie temperature T Curie of the FI layer is larger than the effective exchange field, whose maximum value in the superconducting state is smaller than the superconducting order parameter ∆. For example, in EuS/Al junctions T Curie ∼ 16.6 K and ∆ 0 ∼ 0.25 meV, thus k B T Curie /∆ 0 ∼ 5.3. In other FIs such as EuO or NbN, the Curie temperature is even larger and hence this ratio increases. Thus the approximation made in Sec. 2.1 is well justified.
In the case of a weak magnetic stiffness of the FI compared to the characteristic superconducting energy ∆, one should take into account the effect of the superconducting condensate on the magnetic ordering mediated by the Hamiltonian in Eq. (2). The competition between superconducting and magnetic ordering in such a case has been first considered in Ref. [92]. They considered a ferromagnetic state mediated via the Ruderman-Kittel-Kasuya-Yosida (RKKY) indirect exchange interaction. Above the superconducting critical temperature T c the ground state of the system is a homogeneous magnetic state. When the temperature is lowered below T c the ground state of the system corresponds to a spatially inhomogeneous magnetic structure characterized by a modulation with a wave vector of the order of (λ 2 F ξ s ) 1/3 , where λ F is the Fermi wavelength and ξ s is the superconducting coherence length. This state was called by Anderson and Suhl the cryptoferromagnetic state and may describe the situation of certain ternary rare earth compounds [93]. Such a magnetic modulation was also predicted for thin metallic [50,51] and insulating [49] ferromagnetic films on top of bulk superconductors.
Although in this review we assume that the condition k B T Curie /∆ 0 1 is satisfied, and hence no cryptoferromagnetic state is induced, it is important to keep in mind that such a state may exist in the case of FIs with weak magnetic stiffness.

Quasiclassical theory of diffusive FI/S structures
In this section, we outline a microscopic theory useful for describing nonequilibrium effects in spin-split superconductors. We concentrate on diffusive superconductors where scattering by non-magnetic impurities results to a mean free path that is small compared to other lengths involved in the problem, with the exception of the Fermi wavelength λ F . This condition is equivalent to the assumption that the inverse of the elastic relaxation time 1/τ is much larger than all characteristic energies apart from the Fermi energy.
Within a quasiclassical theory framework, which operates on length scales large compared to the Fermi wavelength, the superconducting properties can be described by a diffusion-like equation [94], now called the Usadel equation, which describes the behavior of the Keldysh Green's function (GF), Here D is the diffusion coefficient,Σ =Σ(r, ε) is the general self-energy term anď g(r, ε) is the (momentum isotropic part of the) quasiclassical Green's function, obtained by integrating the microscopic Green's functions over the quasiparticle energy. In addition to (8), the GF satisfies the quasiclassical normalization The derivation of the Usadel equation (8) is discussed in several previous works [95,94,96,2], and we do not present it here. For the readers not familiar with the quasiclassical method it is sufficient to understand that the knowledge of the structure of the GF sheds light on the spectral and transport properties of superconductors. In particular, the self-energiesΣ generate the collision integrals for the different scattering processes of the kinetic transport theory, as described in more detail in Sec. 3.1. In Eq. (8), the matrices τ j , j = 1, 2, 3 are Pauli matrices in the Nambu (or electron-hole) space. 3 The matrices σ j (j = 0, 1, 2, 3), are the Pauli matrices in spin space, with σ 0 denoting the unit matrix. The vector h denotes the spinsplitting field generated either by an external field or the magnetic proximity effect in a FI/S junction, and∆ = ∆τ ↑ + ∆ * τ ↓ is the singlet superconducting order parameter.
The Green functions and self-energies are functions of the spatial coordinate and the energy. They are also matrices in the Keldysh⊗Nambu⊗spin space. The 4 × 4 Nambu⊗spin matrix structure corresponds to the Nambu bi-spinor, which we in this review choose as The structure in the Keldysh space can be represented via three independent components. Here we use the representatioň where the retarded (R), advanced (A) and Keldysh (K) components are 4×4 matrices in the Nambu ⊗ spin space. The set of equations is completed by expressions for the self-energyΣ and a self-consistency relation for ∆. The former is discussed in subsequent sections, and the latter for conventional superconductors can be written as Here, λ is the effective coupling constant and Ω D is the Debye cutoff frequency. Quite generally, the R/A components of the Green function are related to spectral properties of the superconductors whereas K also contains information about the quasiparticle kinetics. Indeed, it can be related to the electron distribution functionf ,ĝ In the most general case, the 4×4 matrix distribution functionf can be written asf where f T /L,j are real-valued. Here, we generalize the notation introduced in Ref. [97] to the spin-dependent case. The different components off carry specific information of the nonequilibrium state of the electron system. In equilibrium at temperature T the distribution function isf =1 tanh(E/2T ). We discuss the physical interpretation of each component in the following section. For a homogeneous bulk superconductor, with no spatial dependence in the quantities, we can write the general form of R and A functions aŝ The off-diagonal components, i.e., those proportional to τ 1 , are characteristic of the superconducting state and describe the anomalous GFs. The component of the anomalous GFs proportional to σ 0 describes the usual singlet correlations in the BCS theory. On the other hand, the terms proportional to the Pauli matrices σ j , j = 1, 2, 3, describe the three triplet components of the condensate that appear in the presence of spin-dependent fields [79,2]. The density of states, normalized by its normal state value, can also be written in terms of them: Magnetic behaviour of dirty multiband superconductors near the upper   which for a bulk spin-split superconductor reduces to Eq. (4). Accordingly, the difference of the DOS of the spin split bands is determined by The latter quantity is shown in Fig. 4. The real and imaginary parts of the different components of the retarded Green's function are plotted in Fig. 6 in the case of a homogeneous magnetization pointing in the z direction. In that case, the GF components proportional to σ 1,2 vanish.
In what follows, we discuss the nonequilibrium quasiparticle physics in more detail.

Nonequilibrium quasiparticles
We call the different components f T /Lj (ε), j = 0, 1, 2, 3 of the distribution function the nonequilibrium modes. They provide a description of the different ways in which the electron distribution in the quasiclassical limit can deviate from equilibrium.
The excitation spectrum of a spin-split superconductor has distinct electronand hole-like branches, each with spin either up or down. The exchange field h splits the energies of the spin branches, as illustrated in Fig. 7(a). In a collinear situation with a single spin quantization axis, the electron population can be described with the four distribution function components f T , f L , f T 3 , and f L3 .
The mode f T ≡ f T 0 [ Fig. 7(b)] is well-known in superconductor physics: it is the charge imbalance (or "branch imbalance") mode [98,99,100], which corresponds to an imbalance in the number of quasiparticles above and below the Fermi surface. This mode can be generated by charge injection into the superconductor for example from a normal-state electrode, and measured by observing the local potential µ of the quasiparticles. The mode f L is associated with the energy content of the superconductor, and reflects local changes in the effective temperature. It can generally be induced by any heating mechanism.
The spin imbalance mode f T 3 [ Fig. 7(d)] is related to a spin accumulation µ z , which can be nonzero also in the absence of spin splitting of the spectrum. It can be induced for example by a spin-polarized injection from a ferromagnetic electrode, even in the normal state [101]. Finally, the spin-energy mode f L3 encodes a nonequilibrium state with antisymmetric differences in the electronhole and spin-up/down distributions. We are not aware of existing experiments probing f L3 nonequilibrium.
The spin splitting of the superconductor spectrum also modifies how the modes contribute to observables, which generally depend on the amount of quasiparticles on the different spectral branches. The charge imbalance potential µ then acquires contributions from both electron-hole antisymmetric modes f T and f L3 . The other two modes, f L and f T 3 , contribute to the spin accumulation µ z .
To excite the modes f T , f T 3 and f L3 , one only needs to transfer quasiparticles between the different spectral branches in an elastic process, i.e., between equal-energy states (horizontally in Fig. 7). They can also relax back to equilibrium due to elastic scattering processes (horizontal arrows in Fig. 7). These relaxation mechanisms depend on properties of the material, and also on the superconducting spectrum, and are discussed in subsequent sections.
In contrast, relaxation of the f L mode generally requires inelastic processes, which often are slow compared to elastic processes at low temperatures. As a consequence, perturbations in f L can survive up to longer times and propagate longer distances than those in the other modes. Because the f L mode also contributes to the spin accumulation, this is crucial in understanding long-range spin signals observed in spin-split superconductors, for example in Ref. [31], which is discussed in Sec. 4.

Kinetic equations
The different parts of the matrix distribution functionf are determined by a kinetic equation that describes the balance of different transport and relaxation processes. In a diffusive conductor in the steady-state limit, the kinetic equation can be obtained by inserting Eq. (13) in the Usadel equation (8) and evaluating its different matrix components. This results into a set of diffusion-type equations of the generic tensor form Here, I ab coll = Trτ b σ a [Σ,ǧ] K /8 is the collision integral describing scattering processes with self-energyΣ. The term H ab = Trτ b σ a [−ih·στ 3 ,ĝ K ]/8 describes the Hanle spin precession due to the exchange field, and R ab = Trτ b σ a [∆,ĝ K ]/8 the quasiparticle branch conversion processes involving the superconducting condensate. The detailed forms of the different components of the kinetic equations in spin-split superconductors are discussed in Sec. 4.1. Note that the definition of the collision integrals here depends on the form of the diffusion equation adopted; another definition is discussed in Sec. 4.4.
The above equations are formulated for a stationary situation where ∆ is time-independent, and here and below, we quite generally only consider this case. The equations are written in a gauge where the zero of the electric potential is taken to be the potential of the superconducting condensate. We also consider only length scales where the metal can be assumed locally charge neutral. The formulation still allows describing normal-state electrodes at nonzero electrochemical potential, by including distribution functions of the form [96] f eq,L(T ) (ε) = as boundary conditions. Superconducting electrodes however have V = 0 as the quasiparticle and condensate potentials coincide at equilibrium.

Observables
Local densities of charge, and excess spin and energy can be obtained from the Keldysh Green function. These equations in the quasiclassical framework read for the charge imbalance and spin accumulation and for the local energy and spin-energy accumulations: Here, polarization direction of the spin is encoded in the direction components a = 1, 2, 3. The energy content is written relative to a system at equilibrium at the chemical potential µ S = 0 of the superconducting condensate. The charge imbalance is related to the local charge density by ρ = −ν F e 2 φ− eν F µ, where φ is the electrostatic potential. In the locally charge-neutral limit considered here, ρ ≈ 0 so that −eφ ≈ µ and µ corresponds directly to the local electric potential [102,103].
By using Eqs. (13)(14)(15)(16)(17) one can express the charge and spin accumulations in terms of the distribution functions and the density of states, Similarly for the local energy and spin-energy content: The above equations apply to the collinear situation with a single spin quantization axis.
The observable current densities of charge, energy, spin, and spin-energy can be obtained from the spectral current (19). Explicitly, they read, (i) the charge current: (ii) the spin current density polarized in a-direction (iii) the energy current density and (iv) the spin energy current density Here, σ N = e 2 Dν F is the normal-state conductivity and ν F is the normal-state density of states at the Fermi level.

Elastic relaxation mechanisms
In addition to the impurity scattering that results in the momentum relaxation, there are also other elastic scattering processes that, while preserving the energy of the excitations, change their other quantum numbers such as the spin. For example, in metals and semiconductors there are two main types of mechanisms that relax the spin: scattering from magnetic impurities and spin-orbit coupling. As discussed in Sec. 2, these processes can modify substantially the spectral properties of a superconductor. In addition, the superconducting properties are modified by the orbital effect due to external magnetic fields, which for the quasiparticles can result to an electron-hole branch-mixing process. Within Born approximation, the self-energies for these elastic processes obtain the formsΣ Here, τ so and τ sf are the two scattering times for impurity spin-orbit and spinflip scattering. We parametrize these scattering rates in terms of the total spin relaxation time, τ −1 sn = τ −1 so + τ −1 sf and the relative strength of the two scattering mechanisms β = (τ so −τ sf )/(τ so +τ sf ). Values with β ≈ 1 correspond to dominant spin-flip scattering, whereas β ≈ −1 to dominant spin-orbit scattering.
The orbital self-energy above can be used to describe the orbital depairing effect of (Meissner) screening currents induced by an in-plane magnetic field in a thin film, leading to the suppression of superconductivity [104]. It does not enter explicitly the kinetic equation for the distribution function, but modifies indirectly its energy dependence by affecting the spectral properties of the superconductor, contributing to the relaxation of charge imbalance [97,105]. In the case of a thin magnetic film with thickness d and an in-plane applied magnetic field, the orbital depairing relaxation time is given by [106] 1/τ orb = De 2 B 2 d 2 /6.
In the case of spin-orbit coupling we can distinguish two types of relaxation mechanisms according to the origin, intrinsic or extrinsic. Intrinsic spin-orbit coupling occurs in systems without inversion symmetry, either due to geometry constraints or the crystal potential associated with the electronic band structure. Momentum-dependent spin precession together with random momentum relaxation due to impurities leads to the Dyakonov-Perel relaxation [107]. For example, for a diffusive 2D system with a Rashba spin-orbit coupling α R , the Dyakonov-Perel relaxation time is τ = 1/D(2mα R ) 2 for an in-plane spin-polarization, and τ ⊥ = τ /2 for spins perpendicular to the 2D system.
In contrast, the extrinsic spin-orbit effect originates from a random impurity potential and hence leads to an isotropic spin relaxation known as the Elliott-Yafet relaxation mechanism. It is this component that is described by the self-energy Eq. (33). A detailed discussion of spin relaxation mechanisms originating from the spin-orbit coupling can be found in the review [108] and references therein. Because we here focus on centrosymmetric materials, we only consider the extrinsic relaxation mechanism. Moreover, we do not consider effects related to the spin-charge coupling, such as the spin Hall effect [109] and the spin galvanic effects [110]. Such effects enter the kinetic equations in a higher order of the gradient expansion and can be neglected in the leading order. The spin-charge coupling in diffusive superconductors leads to non-dissipative  and spin-flip (β = 1) relaxation mechanisms. We assume that the normal-state spin relaxation rate is τ −1 sn = 0.2T c0 .
magnetoelectric effects and the appearance of an anomalous phase in Josephson junctions, [111,112,113]. The parameters τ sf , β of the spin-dependent scattering are material specific. In Refs. [114,115], using spin injection experiments, the values τ sn ≈ 100 ps and β ≈ 0.5 were found for Al for which then τ sf < τ so . According to Refs. [116,73], the spin-orbit scattering rate is the momentum relaxation rate times (Zα) 4 , where Z is the atomic number and α is the fine-structure constant. Therefore, the spin-orbit scattering rate grows rapidly as the atomic number grows. Indeed, in Nb τ sn is only 0.2 ps and spin-orbit scattering dominates [117].

Nonequilibrium spin relaxation
In this review we only consider centro-symmetric metals, and therefore we only consider spin-relaxation caused by impurities. In this case, the collision integrals entering the kinetic equation due to spin-dependent scattering are given by where the components a, b refer to the different nonequilibrium modes and spin projections. The self-energies are defined in Eqs. (33) and (34). In the superconducting state, the form of the spin-flip and spin-orbit contributions to the collision integrals differ due to their different properties under time reversal. We now discuss the relaxation of the nonequilibrium modes related with the spin and the spin energy imbalances, i.e., the components f Specifically, the spin relaxation, shown schematically by the arrows in Fig. 7 In the normal state these processes are energy independent (within the quasiclassical approximation), and have the rate τ −1 sn . In the superconducting state, the spin-relaxation is described by Eq. (37), with the energy-dependent rate Here the collision integral is normalized with the spinaveraged density of states N + to counterbalance effects due to the mere changes of the quasiparticle spectrum in the superconducting state. The typical behavior of τ −1 s (ε) is shown in Fig. 8 for dominating spin-flip and spin-orbit mechanisms for (a) zero and (b) nonzero spin splitting.
In the presence of strong spin-orbit scattering (β < 0), superconductivity results to reduction of the spin relaxation rate τ −1 s . Still, there is no significant change in the spin relaxation length λ so : [118] Qualitatively λ so ∼ v g τ s , and the reduction in the quasiparticle momentum scattering cross section is mostly canceled by the decrease in the group velocity v g ∼ v F 1 − |∆| 2 /ε 2 . In contrast, the increase in τ −1 s due to spin-flip scattering does result to a reduced spin relaxation length: [118] this mechanism is not related to the momentum scattering as the interaction with magnetic impurities does not depend on the propagation direction and the quasiparticle spin does not depend on energy.
In the presence of a Zeeman or exchange field, h = 0, this situation changes drastically. In this case the nonequilibrium spin accumulation couples to the energy mode, f L . This mode is robust with respect to elastic spin relaxation processes, resulting to a qualitative change in the observed spin relaxation. This effect, as well as the detailed energy dependence of the spin relaxation lengths, are discussed in detail in Sec. 4.2.
For the relaxation of the spin-energy mode f L3 , the situation is more complicated. As sketched in Fig. 7(e), the relaxation of this mode involves transitions between electron and hole branches of the quasiparticle spectrum. As in the case of charge imbalance relaxation, such inter-branch transitions involve the formation of Cooper pairs via Andreev reflection. This results in two contributions to the relaxation of the spin-energy mode: One stems from the (j = 3, k = 3) component of the collision integral (36) that describes transitions between spin-up and spin-down subbands S s The second contribution originates from the particle-hole relaxation described by the non-diagonal self-energy in the Keldysh-Usadel equation (8), The relaxation of the spin-energy mode in the superconductor is given by the sum R T + S L3 . Note that the particle-hole conversion processes responsible for the spin energy relaxation also lead to the generation of charge imbalance producing a coupling between spin and charge degrees of freedom. These mechanisms are discussed in more detail in Sec. 4.2. At large energies ε ∆ the charge and spin energy nonequilibrium modes are decoupled. In this limit the charge relaxation length tends to infinity because the coherence factor Reg 01 → 0. At the same time, the spin-energy relaxation length approaches its normal-state value equal to that of the spin relaxation length Λ sn = √ Dτ sn .

Inelastic relaxation mechanisms
The inelastic processes that are typically relevant in metallic systems, the particle-phonon and particle-particle collisions, are described by an inelastic self-energyΣ in =Σ eph +Σ ee . These processes do not conserve the energies of colliding quasiparticles, but conserve the total spin. A model self-energy for inelastic relaxation due to electron-phonon scattering is given by [ are parts of Fermi surface averages of the free phonon propagators, T ph is the phonon temperature and we denote X RA = X R − X A for X = D,ĝ. Inelastic particle-particle collision self-energiesΣ ee within the quasiclassical theory have been discussed in Refs. [119,120].

Electron-phonon coupling
The relaxation due to the electron-phonon coupling is an important limiting mechanism for some of the nonequilibrium effects discussed in subsequent sections. Generally, it becomes weaker towards lower temperatures as the phonon density of states decreases at low energies. As a consequence, the temperatures T qp and T ph can decouple and the two subsystems can be in different temperatures, and eventually e-ph interaction is dominated by other relaxation processes. As the interaction is also sensitive to the electronic spectrum, it is modified by superconductivity. [119,103,121] Below we discuss what this results to in the spin-split case.
The collision integral for electron-phonon processes in a spin-split superconductor was discussed in Ref. [122] and can be obtained from the self-energies described in Eq. (41) [123]. The overall collision integral entering the kinetic equations [Keldysh part of Eq. (8)] is obtained from a commutator between the self-energy and Green's function,Ǐ eph = [Σ eph ,ǧ] K . This is still a Nambu-spin matrix, its different matrix components describe the relaxation of the different nonequilibrium modes. Electron-phonon interaction causes inelastic scattering of quasiparticles, which means that energy flow between quasiparticles and phonons becomes possible. That process is described by the collision integral I L ≡ Tr[Î eph ]/8. Also the other matrix components of the collision integral are in general non-vanishing; for example the component Tr[τ 3Îeph ] affecting f T is detailed in [124]. However, since electron-phonon interaction conserves charge and spin, the energy integrals of Tr[τ 3Îeph ] and Tr[σ 3Îeph ] vanish (see Appendix A). The remaining component Tr[τ 3 σ 3Îeph ] describes the relaxation of the difference in the thermal energies of electrons with opposite spins due to electron-phonon interaction [125]. In the following, we only concentrate on energy relaxation. Substituting the self-energy (41) in the equation for the collision integral and assuming thatf = f L1 we can write Typically this collision integral is used to calculate the total heat currentQ eph = 2ν F Ω dεεI L (ε) between the electron and phonon systems in a volume Ω. Assuming that the electron system resides in temperature and thus the heat current is [126,127], Here Σ = 3072ζ(5)ν F k 5 B g eph describes the materials dependent magnitude of the electron-phonon coupling (see tabulated values in [128]) typically used for heat flow, and ζ(x) is the Riemann zeta function. The kernel L E,E depends on the spin-splitting field, as it is where N σ (ε) = Re(g 03 + σg 33 )/2 is the density of states for electrons with spin σ and F σ = Re(g 01 + σg 31 )/2 is the anomalous function for spin σ (see Eq. (15) and below it). The latter formula applies in the absence of spin-flip or spin-orbit scattering.
Strictly speaking, Eq. (43) requires two conditions to be met: (i) A welldefined local quasiparticle temperature that may deviate from the phonon temperature. In other words, the quasiparticles should first equilibrate between themselves. (ii) In the presence of a spin-splitting field, also spin imbalance affects the quasiparticle-phonon heat flow [122,123]. Neglecting it therefore requires that the spin relaxation rate exceeds the injection rate causing the out-of-equilibrium situation. 4 In the general case Eq. (43) needs to be evaluated numerically. However, we may consider two limiting cases. In a normal metal ∆ = 0 implies L ε,ε = 1. In that case the heat current between the two systems becomes [129] Another tractable limit is that of linear response, where The two terms in Eq. (46) describe scattering and recombination processes, respectively. Equation (46) is compared to the exact result in Fig. 9. The recombination process, whose heat conductance is almost independent of the spin splitting (as long as h ∆), dominates at high temperatures, whereas the two become of the same order of magnitude for T ≈ 0.1∆. Note that the exponential suppression of the electron-phonon heat conductance ∼ e −∆/T at low temperatures does not directly make it irrelevant in a superconductor, because any effect related to quasiparticles contains such exponential terms due to the gap in the superconducting density of states. However, as lowering the temperature reduces the phase space available for acoustic phonons, also the prefactor of the exponential becomes small low temperatures. Thus, at very low temperatures quasiparticles decouple from phonons, and other heat conduction mechanisms become relevant.

Particle-particle collisions
The particle-particle collision integrals generally reflect symmetries and conservation laws of the electron system. Below, we only point out conservation laws that a consistent formulation of the collision integrals must imply in a system involving collinear magnetizations. The conservation of charge, spin, and overall energy by spin-rotation symmetric scattering within the particle system can be expressed as: The spin energy, involving the component ετ 3 σ 3 of the collision integral, is generally not conserved in electron-electron collisions, because electrons with opposite spins can exchange energy [125]. The specific forms for I ee for superfluids have been discussed in [119,103,120], however, typically in the absence of spin splitting of the spectrum. The above conservation laws are relevant in terms of studying what happens in the limit of strong electron-electron scattering. In that case, the usually assumed model is that the quasiparticle distribution functions tend to equilibrium forms, but retaining some effective temperature and a spin-dependent chemical potential. This quasiequilibrium limit (without spin accumulation) is often used in analyzing heat transport in superconductors driven out of equilibrium (for example, see [130]).

Description of hybrid interfaces
The spectral current density (19) appearing in the kinetic equation is defined in the bulk of a superconductor (S) or a normal metal (N). In addition, we need a description of interfaces between different materials in order to model transport in hybrid structures. At such interfaces, potentials and materials parameters vary over atomic distances, and therefore they cannot directly be described by quasiclassical equations that are valid only when changes occur over distances much larger than λ F . Rather, such interfaces are included in the theoretical description via the derivation of suitable boundary conditions for the Green's functions. This was first done in Ref. [131], and later extended in Refs. [132,133,134]. Magnetic interfaces were discussed in Ref. [14] for the S/FI case and extended to metallic interfaces and general cases in Refs. [135,136,137,29,138].
To get a qualitative understanding of the phenomena related with spinpolarized tunneling into superconductors, we consider below low-transmissive interfaces between superconductors and normal or magnetic leads, using the description in Ref. [137]. Fixing the z direction in the polarization direction of the interface, the transmission through it can be described by the tunneling matrixΓ = tτ 3 + uσ 3 . It is defined through the normalized transparencies satisfying t 2 + u 2 = 1 and determining the interface polarization via 2ut = P . The possible non-collinearity between the interface and the bulk magnetizations can be described by introducing rotation matrices in spin space. We assume that this interface connects a normal-metal electrode with bulk Green's functionǧ N with a spin-split superconductor with GFǧ. The boundary condition gives an expression for the matrix current density through the interface in terms ofΓ, g N andǧ:ǧ where R characterizes the spin-averaged barrier resistance per unit area. Its precise microscopic definition is given in Refs. [137,138]. In experiments it is usually used as a fitting parameter. Equation (48) is valid for an arbitrary spin polarization P ∈ [−1, 1] of the interface. From this condition and Eqs. (29)(30)(31)(32) we obtain the spectral current densities through spin-polarized barriers as The normal electrode Green's functionǧ N is defined byĝ is the distribution function of the normal metal. The boundary conditions for both spectral GF (R) and distribution functions (K) are then given by the continuity of the matrix current density, n k j a kb = I ab , where the tunneling current I ab and bulk current j a kb are provided by Eqs. (49) and (19), respectively, and n k are the components of the unit vector normal to the interface. For small enough junctions, the approximate conservation of matrix current at short distances implied by Eq. (8) [134] allows neglecting the precise 3D structure of the junction. Hence, in quasi-1D problems, the boundary condition can be integrated over the cross-sectional areas of the wire A and tunnel interface A T and written in the form Aj a 1D,b = A T I ab . Low effective transmissivity of interfaces, which corresponds to a large interface resistance (tunneling limit), also allows for further simplified approximations to be made in the description of nonequilibrium effects. In such a case, changes in the spectrum due to the coupling can be neglected, and the probability for Andreev reflection [83] vanishes. In the lowest order in transmissivity, we can also assume the electrodes to be in local thermal equilibrium. If, in addition, P = 0, the charge current between two electrodes is given by the usual tunneling expression Eq. (3) with G T = A T /R . In Sec. 4 we describe the more general limit of an arbitrary P , and the different components of the spectral current.
It is worth emphasizing that the boundary condition, Eq. (48), when combined with the quasiclassical equations, allows for the description of effects which are beyond the quasiclassical limit and that involve strong ferromagnets with large spin polarization [139,140].

Free energy
The quasiclassical approach can be used to describe the collapse of superconductivity caused by an applied exchange field, including the part of the parameter space where a first-order phase transition occurs. The main physical features of this are discussed in Sec. 2.3. In particular, a necessary condition for the stability of a superconducting state is that the free energy density of the superconducting state compared to the normal state is negative, Ω < 0. Below, we express this condition in the quasiclassical framework, neglecting the FFLO state.
The free energy density of a uniform superconductor relative to the normal state can be written in terms of the quasiclassical Green functions, [95,120,141,142,143]: where the auxiliary Green functionĝ λ [Σ] satisfies Here,Λ = (ω m + ih · σ)τ 3 , and Tr = πT N (0) ωn tr contains a Matsubara sum and Nambu and spin traces. Variation of Ω vs.ĝ produces the relation δΦ[ĝ] = − 1 2 Tr[Σ * δĝ] between the self-consistent self-energiesΣ * and the subtracted quasiclassical Φ-functional, with value Φ[ĝ * ,N ] = 0 in the normal state. Variations and derivatives ofĝ need to retain the conditionĝ 2 = 1, and can generally be expressed in the form δĝ = δXĝ −ĝδX for some δX. This observation, and making use of Eq. (51), allows the evaluation of the λ integral: This enables writing the self-consistent free energy as Note that in a spatially non-uniform situation, gradient terms would complicate the λ-integration. This still allows for including the magnetic orbital effect, given the approximation in Eq. (8) where Eq. (51) retains its spatially uniform form also in the presence of the vector potential [144]. The Φ-functional producing the elastic scattering self-energies discussed in preceding sections reads For τ orb , τ so , τ sf ∆ −1 , direct evaluation yields where w n = ω n +ih. At low temperatures , producing the Chandrasekhar-Clogston result discussed in Sec. 2.3.

Linear response and generalized Onsager relations
A way of exciting the nonequilibrium modes in superconductors is by injecting electrons from either a voltage or temperature biased metallic lead. In this case, and in the language of Eq. (14), the distribution function in the normal metal is given byf ]/2 is the voltage-biased distribution function in the normal-metal electrode and f eq (ε) = tanh(ε/2T N ) is the equilibrium distribution corresponding to the normal-metal temperature T N .
In a more general situation, all components of the distribution function, Eq. (14), have to be taken into account. The different components of the spectral current through a tunneling contact between a normal electrode and a superconductor can be derived from the boundary condition in Eq. (49). In the case of collinear magnetizations the relevant components are the spectral charge I 30 , energy I 00 , spin I 03 and spin energy I 33 currents. They satisfy where the parameter κ = 1/(R e 2 ) describes the interface transparency and , with k = T, L, T 3, L3, between the superconductor and normalmetal electrodes. The response matrix is here described by the spin polarization P and the energy-symmetric and energy-antisymmetric parts of the density of states, N ± defined in Eqs. (16)(17).
Expression (57) is valid for arbitrarily large deviations from the equilibrium state and describes several intriguing effects which come into play in the nonequilibrium situation in junctions combining spin-split superconductors (N − = 0) and a spin-polarized injector, modeled by a nonzero spin polarization P = 0. If both N − = 0 and P = 0 simultaneously, all the nonequilibrium modes are coupled. In turn, this coupling results in a rich variety of cross-couplings between different potentials and currents.
To illustrate the basic phenomena which can be expected in such a tunneling contact, let us consider the linear response limit. The different components of the distribution function can be described by spin-dependent temperatures 5 and voltages introduced according to Ref. [145]: where V and ∆T are the voltage and the temperature bias and V s and ∆T s the spin-dependent biases.
Integrating Eq. (57) over energy allows us to calculate the total charge I = e where the conductance, heat conductance and thermoelectric coefficient are given by respectively. 7 In the absence of spin-orbit or spin-flip scattering, these integrals can be evaluated analytically in the limit T ∆ − h. The results are in Eqs. (147).
Note that the four 2×2 quadrants of the response matrix have different timereversal symmetries: the quadrants on the diagonal are symmetric, whereas the off-diagonal quadrants coupling spin and charge are antisymmetric. This result follows from the fact that G and G th are symmetric whereas α and P are antisymmetric in time reversal.
From Eq. (59) we can, for example, study the spin and charge currents created at the SF interface in response to a temperature bias between the normal lead and the superconductor. These are the thermospin and thermoelectric effects. Interestingly, Eq. (59) demonstrates also that the charge current can be induced by the spin-dependent temperature bias ∆T s even without spin-filtering (P = 0) but in the presence of a spin-split DOS, N − = 0. Qualitatively the corresponding non-equilibrium mode f L3 [ Fig. 7(e)] can be interpreted as the spin-dependent particle-hole imbalance which produces electric signal due to N − = 0.
Let us focus on the case where the nonequilibrium state is generated by applying a voltage and a temperature bias to the injecting normal electrode in the absence of spin-dependent potentials, ∆T s = 0 and V s = 0. If the interface lacks spin polarization, P = 0, then from Eq. (59) it follows that all currents are decoupled from each other, recovering standard expressions for charge and heat currents in terms of the local electrical and thermal conductances I = GV andQ = −G th ∆T /T . In addition, the remaining block of the Onsager matrix leads to nonzero spin current I s = −α∆T /T and spin heat currenṫ Q s = αV . From this linear response analysis we can conclude that in a setup consisting of a normal and a spin-split superconducting electrode, a spin current can be generated even in the absence of ferromagnetic electrodes. It is worth mentioning that this thermospin effect leads to a pure spin current I s that does not carry electric charge and therefore does not produce Joule heating and Ohmic losses. Clearly this situation is very different from the one occurring in normal metal systems, where spin currents can be generated only by injecting a charge current from a ferromagnet -a process corresponding to the relation I s = P GV obtained from Eq. (59).
If P = 0 all four types of currents are injected into the superconducting region just by applying the temperature bias ∆T . In particular a very large thermoelectric effect can be generated, described by I = −P α∆T /T . This effect is discussed in detail in Sec. 5.

Nonequilibrium effects on the superconducting order parameter
Nonequilibrium also affects the superconducting state by modifying its order parameter. This can be seen by inserting the parametrization (13,14) of the nonequilibrium Green's function in Eq. (12): In equilibrium, only the first term contributes with The imaginary terms affect the phase of ∆ and ensure charge current conservation in situations where quasiparticle current is converted to supercurrent (described by the R ab term in Eq. (18)). In the absence of spin splitting, the nonequilibrium modifications in the size of the gap, described by the first term, are the most relevant. Often such modifications are studied within the Rothwarf-Taylor phenomenological model [146] that neglects the energy dependence of the distribution functions, and rather concentrates on the overall number of quasiparticles. This model was derived under some simplifiying assumptions from the full energy dependent kinetic equations in Ref. [147]. It has been used to study the gap suppression due to nonequilibrium injection, see for example Ref. [148]. According to that work, for small changes around equilibrium, such models agree with the predictions of the phenomenological models of Refs. [149,150]. For large changes, the full kinetic equations with energy dependent distribution functions should be used. An extreme example is given in Ref. [151] describing a voltage-driven transition of a superconducting wire to the normal state.
In the spin-dependent case we need to consider both the effects of spin splitting and the spin-dependent distribution functions describing spin accumulation. Disregarding the terms affecting mostly only the phase of the order parameter, we can also write Eq. (63) as where F σ = (g R 01 + σg R 31 )/2 and f σ = (f L + σf T 3 )/2; the different signs of σ = ± representing the two spin directions. Ref. [152] went a step further, assuming that the spin accumulation in a superconductor can be described in terms of a simple spin-dependent chemical potential shift µ s , neglecting all other nonequilibrium effects. In that case f σ = f eq (ε − σµ s ). As a result, by a simple shift of the energy (ε → ε + σµ s ) in Eq. (64), µ s shows up as a shift of the energy of the anomalous function F σ and a small shift of the cutoff Ω D . When Ω D µ s , the latter effect can be disregarded, and the net effect of the spin accumulation is the same as that of the spin splitting field, explained in Sec. 2, eventually leading to the suppression of superconductivity. Moreover, in the presence of both spin splitting and spin accumulation, this model shows how in the special case h = µ s spin accumulation can actually lead to the recovery of superconductivity suppressed by spin splitting [153].
This physics can be probed in a FISIF system where a superconducting island or layer is placed between two ferromagnetic electrodes. In this case, the current induced suppression of the superconducting gap should be larger when the magnetizations of the two ferromagnets are antiparallel than when they are parallel [152]. Only in the previous case the spin accumulation builds up. This effect was measured in Ref. [154]. They indeed found a stronger suppression of superconductivity in the antiparallel configuration. However, they found that the required spin relaxation time for fitting the results to the above theory is much longer than that expected from the normal-state measurements. This discrepancy may result from the somewhat simplified model for the spin accumulation described above.
The suppression of superconducting properties due to spin injection has been measured in the case of high-temperature superconductors [155,156,157,158]. However, these features are typically attributed to the non-conventional character of superconductivity, and are therefore outside the scope of this review.
Besides the superconducting gap, in principle also the induced spin-splitting field can obtain nonequilibrium corrections in the presence of spin injection. The theory for such effects was outlined already in Ref. [75], but to our knowledge such effects have not been thoroughly examined in spin-split superconductors.
3.8. Overall strategy to explore the transport and spectral properties of hybrid superconducting systems Spectral and transport properties of a diffusive metal or superconductor of mesoscopic size attached to electrodes can be fully described by using the theoretical framework presented in this section. Specifically, one should solve the boundary problem defined by the Keldysh-Usadel equation (8), the normalization condition forǧ, Eq. (9), and the boundary conditions at the interfaces with the electrodes, Eq. (48). Once the Green's functions are determined one can compute the currents and potentials from Eqs. (21)(22)(23)(24)(25)(26)(27)(28)(29)(30)(31)(32).
If the system under consideration is in the normal state, the equations for the retarded and advanced GFs are decoupled from each other and from the Keldysh one. The solution of the spectral equations within the quasiclassical approach is trivial and leads toĝ R(A) = ±τ 3 . Thus the problem reduces to solving the equation for the distribution functions.
If the system is a superconductor attached to leads, the equations for the retarded and advanced GFs are coupled to the Keldysh component through the self-consistency equation for the superconducting order parameter, Eq. (12). This complication can be overcome if the superconductor and the electrodes are coupled via tunneling junctions. In such a case the spectral properties of the superconductor remain unchanged within the leading order in the interface transmission. This means that retarded and advanced GFs coincide with those of the homogeneous superconductor. In Sec. 4, we concentrate on the case where the normal-state tunnel junction conductance satisfies G T R E 1, where R E is the normal-state resistance of a wire with length equal to the energy relaxation length (due to electron-phonon or electron-electron scattering) E at the energy around the amplitude of the superconducting gap. On the other hand, in Sec. 5, we mostly consider the cooling of the superconducting island with a starting temperature much below T c . In these cases we can disregard the nonequilibrium effects on the energy gap, and can rather use the self-consistent gap that is calculated from the equilibrium version of self-consistency equation, i.e., Eq. (63) In most of the discussed examples of this review we follow this approach and use the self-consistent ∆ calculated in the decoupled spin-split superconductor. From this we obtain the spectral functions of the spin-split superconductor that yield the coefficients of the kinetic equation for the nonequilibrium distribution functions. The latter is obtained from the Keldysh component of Eq. (8). Only in Sec. 4.5, where we study the effect of an ac field on the spectrum of a spin-split superconductor, we compute the self-consistent gap using the nonequilibrium distribution function (see Fig. 20).
Inspection of the different components in Keldysh space of Eq. (8) provides a first insight about the characteristic lengths involved in the different situations. In a superconductor changes of the spectral properties due, for example, to the magnetic proximity effect with a ferromagnetic insulator, occur over the superconducting coherence length ξ s . On the other hand, this ξ s depends on the temperature, and on the concentration of magnetic and spin-orbit impurities via the self-consistent ∆.
Different length scales govern the decay of nonequilibrium components of the distribution function generated, for example, by injecting a current through a contact with an electrode. In this case the characteristic lengths depend on the nature of the excited mode (charge, energy or spin) and the type of relaxation process in the system. The calculation of these characteristic lengths is one of the main goals of Sec. 4. As we show there, the spin splitting plays a crucial role in determining these length scales, since it couples the different nonequilibrium modes, and hence changes the range over which non-local spin and charge signals can be detected.

Transport properties of ferromagnetic insulator-superconductor heterostructures
Electrical injection of spins into a superconductor was first studied in Ref. [159]. The injection of an electric current from a F electrode into a superconductor creates not only spin, but also charge imbalance. Moreover, such imbalances relax away from the injection point over distances different from those in the normal state.
In the absence of a spin-splitting field in the superconductor the spin imbalance is decoupled from the other modes. However, in spin-split superconductors charge, spin and energy modes couple with each other. The goal of this section is to describe this coupling, the characteristic relaxation lengths of the different nonequilibrium modes and effects that occur in superconducting structures as a consequence of such coupling.
We start by reviewing experiments on charge and spin injection in superconductors and by summarizing the main theoretical works on spin injection in superconductors with no spin splitting. and discuss several possible experimental situations both with collinear and non-collinear spin configurations.

Detection of spin and charge imbalance: non-local transport measurements
Nonequilibrium states in superconducting systems can be experimentally accessed using the non-local transport measurement setup invented in Ref. [98]. The geometry used in that experiment was a vertical structure such as the one sketched in Fig. 10(a). Injection of a current from a normal metal (injector) into the superconductor (blue stripe) generates the charge imbalance mode shown schematically in Fig. 7(b). Qualitatively this mode can be understood as the difference between the chemical potentials of the quasiparticles and the superconducting condensate. This relative shift of chemical potentials induces an electric signal (voltage or current), depending on the measurement setting at the detector electrode [green stripe in Fig. 10(a)]. First experiments demonstrating the charge imbalance induced in superconductors by current injection were performed at temperatures close to the critical temperature T c by using vertical samples with large cross sections, such as the one depicted in Fig. 10(a). More recent experiments have allowed for non-local transport measurements at lower temperatures using lateral structures such as the one sketched in Fig. 10(b) [160,161,162,163,38,115]. Such structures allow for an accurate measurement of the charge imbalance and nonequilibrium energy mode and their spatial dependence by placing contacts (detectors) at different distances from the injector. For example, in an aluminum wire at 100 mK the charge and energy modes decay over 5 and 10 µm, respectively, as reported in Ref. [164].
Substituting the normal injector by a ferromagnet, the injected current becomes spin polarized and the charge injection is accompanied by a spin injection into the superconductor. The nonequilibrium charge and spin relax over different lengths. When the superconductor shown by the blue wire in Fig. 10b is in the normal state, charge accumulation is negligible at the detector and only a spin imbalance contributes to the nonlocal signal. The nonequilibrium spin density induced at the interface is polarized in the direction of the injector's magnetization and it diffuses into the normal wire over the spin-relaxation length, which can be several hundreds of nanometers. The spin accumulation can be detected by measuring the non-local voltage between a ferromagnetic detector and the N wire.
First measurement of the nonequilibrium spin accumulation was reported in Ref. [101] on a single-crystal aluminum bar at temperatures below 77 K. This pioneering experiment did not only demonstrate the spin accumulation associated with the injection from a ferromagnetic electrode, but also the coherent spin precession, the Hanle effect, that occurs when an external magnetic field non-collinear with the injector magnetization is applied. Later these effects have been reported at room temperature in Al-based nanostructures in Refs. [165,166]. Since then, electrical injection of spins have been used in several experiments on metallic spintronics devices, for example for the direct detection of the spin Hall effect in metallic structures [167], electronic spin transport in graphene [168], and modulation of the spin density in metal/ferromagnetic insulator bilayers [169,170].
In the superconducting state the situation drastically changes, as reported in several experiments on lateral superconducting structures with ferromagnetic electrodes. In contrast to the charge imbalance, which has been well understood since the 1980s, spin injection and accumulation in superconductors is a more recent research line. In 1990 Kivelson [171] suggested that charge and spin should exhibit different relaxation times in superconductors, leading to the possibility of separating charge and spin transport. Also Ref. [172] studied the spin injection from a ferromagnet to a bulk superconductor and found theoretically that, while the charge imbalance survives only within the field penetration length from the surface, the spin imbalance may also exist in the bulk. However, early experiments measuring spin diffusion length in Nb did not show any evidence of such a long-range spin signal [173,174]. Thus the possibility of having strongly different relaxation scales for spin accumulation and charge imbalance in superconducting state was not confirmed at that time.
More recently, clearer insight into the spin and charge modes has been obtained in experiments using lateral nanostructures with ferromagnetic injectors and detectors [160,175,176,31,177,163,42,115,38,154]. First, the decrease of spin relaxation length in the superconducting state as compared to the normal one has been observed [115] by measuring the non-local spin-dependent resistances. Later, by applying the external spin-splitting field to the superconducting wire it has been possible to prove experimentally the charge-spin separation [32]. A detailed overview of the recent experiments on charge, energy and spin imbalance in superconductors can be found in the recent topical reviews in Refs. [32,178].
In order to describe theoretically the different modes excited in experiments on lateral S/F structures, it is convenient to use the generalized quasiclassical model introduced in Sec. 3.1. First theoretical works on spin injection into mesoscopic superconductors [118,179] addressed the question of how spin re-laxation changes in the superconducting state and how these changes depend on the spin relaxation mechanism. These works pointed out that in the presence of spin-polarized injected currents additional spin-resolved components of the distribution function may appear. The main conclusion of Refs. [118,179] was that the spin-relaxation length changes significantly in the superconducting state and depends on the energy of the injected quasiparticles. In particular, it has been shown that for electrons with energy close to the superconducting gap the spin relaxation length can decrease in agreement with the experiment [115]. However, two important features observed in subsequent studies of spin transport in superconductors with spin-splitting field could not be understood in terms of that theory. One intriguing observation was the drastic increase of the spin accumulation length in the superconducting state as compared to the normal one [38,31,163]. In addition, it has been observed that spin accumulation in such a setup can be created by the current injection from a non-magnetic electrode [163].
As we discuss in this section, to explain these two observations we need to use the kinetic theory which takes into account the modification of the quasiparticle spectrum due to the spin splitting. If it is caused by the Zeeman effect from the external magnetic field, it is also important to take into account the orbital depairing effect which leads to the suppression of superconductivity. Recent works showed that such a modification of the spectrum of the superconductor leads to an intriguing coupling between the nonequilibrium modes in a superconductor [180,35,123,37,36,181]. The effect of this coupling between the modes on the transport properties of a multi-terminal superconducting structure can be theoretically explored with the quasiclassical formalism developed in Sec. 3. Here we review these theory works and provide a quantitative explanation of the long-range spin accumulation detected in multi-terminal superconducting devices.

Nonequilibrium properties of a superconductor with spin splitting
Spin accumulation and transport in superconductors has been considered in a number of theory papers [182,118,179,183,180,35,123,37,36,181]. Here we follow the approach developed in [35] that enables the description of all non-equilibrium modes, the effect of the spin-splitting field, and the various relaxation mechanisms. As discussed in Sec. 3.1, nonequilibrium conditions lead to the excitation of charge, energy and spin modes in superconducting systems. These modes are related to the different components of the distribution function introduced in Eq. (14). A prototypical experimental geometry to reveal the different nonequilibrium modes in superconductors is the lateral structure shown in Fig. 11. It consists of a mesoscopic superconductor attached to four electrodes. One of them serves as an injector for the current, whereas another one is used as a detector. The other two electrodes serve as separate charge current sinks, allowing for such a non-local measurement.
We assume first that both electrodes, injector and detector, are ferromag- netic and that their magnetizations are collinear as in Fig.11a P inj B P det .
The tunnelling current at the detector can be expressed through the charge imbalance µ and spin accumulation µ z (25,26) where R det = R /A is the detector interface resistance in the normal state, and A is the cross-sectional area of the detector. Here we consider collinear magnetizations along the z-axis and therefore only the spin components f T 3 and f L3 enter these expressions. A non-collinear case is discussed in Sec. 4.3.
If the current through the detector is measured at zero bias V det = 0, then the nonlocal differential conductance is defined as This quantity provides information about the modes that, being excited at the injector via the injected current, reach the detector by diffusing in the superconducting wire. The nonlocal conductance g nl depends on the distance L det between the injector and the detector. This dependence reveals the characteristic length scale over which the modes relax and in practice it is measured by using several detectors at different distances to the injector. According to Eqs. (25,26) the non-local current is determined, in the most general case, by the contributions of all four non-equilibrium modes shown schematically in Fig. 7. Here the quantities µ and µ z characterize arbitrary nonequilibrium states of the superconducting wire, including those that cannot be described by the effective potentials V , ∆T , V s and ∆T s introduced in Sec. 3.6. However, if we restrict the analysis to linear response then Eq. (65) reduces to the expression for the current obtained from Eq. (59) with P = P det , When the distance between the injector and detector electrodes is much larger than the charge imbalance relaxation length, the contribution of µ to the detected current can be neglected [36,181]. In the absence of spin splitting h = 0 (N − = 0), the spin-dependent part of I det is determined only by the f T 3 mode related to the spin-dependent chemical potential shift [182,118,179] (see Eq. (58)). In this case in the linear response regime and at low temperatures |eV s |, T ∆, the spin accumulation generated by a spin-dependent chemical potential shift δµ = eV s /2 is exponentially small. This can be seen from the first term on the r.h.s of Eq. (26), by noticing that if the energy gap ∆ is non-zero then N + (|ε| < ∆) = 0. In this case the nonequilibrium spin state is described by f T 3 = (∂f eq /∂ε)eV s /2 and hence the spin accumulation is given by µ z = −(eV s /2) ∞ ∆ dεN + (∂f eq /∂ε) so that µ z ∝ e −∆/T at T ∆ . The non-trivial behaviour of spin accumulation µ z is determined by the second term in the r.h.s of Eq. (26). It appears when the spin splitting becomes non-zero, i.e., h = 0, N − = 0. This term allows for the spin-charge separation which explains experiments in [38,31,163]. On the qualitative level the physics related to this term can be described as follows. When N − = 0 the spin accumulation has a contribution from energy nonequilibrium mode f L shown schematically in Fig. 7. This mode, once excited, can only relax via inelastic processes, especially mediated by the electron-phonon interaction. At low temperatures in metals the electron-phonon relaxation can be much slower than the spin decay in the normal state. This mechanism explains the long-range spin accumulation observed in experiments. Since the nonequilibrium mode f L can be generated in particular by the temperature rise in the quasiparticle system, the long-range non-local spin signals can be explained in terms of the thermospin effect. Its essence is the generation of spin accumulation by heating up quasiparticles in the superconductor. In practice this effective heating can be achieved with the help of the voltage-biased tunnel junction which can inject the nonequilibrium quasiparticles with energies larger than the gap into the superconducting wire [180,35,123,37,36,181]. Such a mechanism does not require injector electrode to be spin-polarized. Thus the spin accumulation can be generated by an injected current even from a non-ferromagnetic electrode. This can be seen from the second term in the r.h.s. of Eq. (26), which shows that µ z can be non-zero even if the injection occurs from a non-magnetic lead, provided that N − = 0, i.e., if the superconductor shows a spin-split spectrum [163].
In addition, with the help of the qualitative picture described above one can understand the antisymmetric shape of the non-local spin signal in g nl with respect to V inj observed in the experiments [38,31,163]. The origin of such a g nl (V inj ) dependence is again the thermo-spin effect in the superconductor with N − = 0. The spin accumulation generated in this way is an even function of the bias voltage at the injector electrode µ z (V inj ) = µ z (−V inj ) because the effective heating of quasiparticles is not sensitive to the sign of V inj . Hence the non-local spin signal (66), that is a derivative of µ z , is an odd function g nl (V inj ) = −g nl (−V inj ). This is also what has been observed in the experiments [162].

Relaxation of the nonequilibrium modes
To provide a quantitative description of the effects discussed in the previous section, we now calculate the potentials µ, µ z and the non-local conductance by using the kinetic equations for superconductors with spin-split subbands introduced in Sec. 3.1.1.
Our starting point is the general Usadel equation (8). By assuming that the transparencies of the interfaces between the superconductor and the electrodes are small we can neglect the changes of the spectral properties of the superconductor due to the proximity effect. This means that the retarded and advanced GFs correspond to the homogeneous superconducting state (cf. Sec. 2), and we only have to focus on the calculation of the components of the distribution function, Eq. (14). In this discussion we disregard inelastic relaxation, assuming the corresponding scattering length to far exceed the spin relaxation length. We thus setΣ in = 0.
In the collinear situation considered here only the four components of the distribution function entering Eqs. (25)(26) are finite. In a homogeneous superconductor, for example in the absence of the supercurrent they are pairwise coupled. The expressions for different currents can be obtained by combining Eqs. (13,14,18,19). It is instructive to represent these expressions in a matrix form  Here, following the notation in Eq. (18), j c · e k = Dj 0 k3 , j se · e k = Dj 3 k3 , j e · e k = Dj 0 k0 and j s · e k = Dj 3 k0 where e k is the unit vector in direction k = x, y, z. The components f L and f T 3 determine the spectral energy j e and spin j s currents, whereas the components f T and f L3 determine the spectral charge j c and spin-heat j se currents . The spectral coefficients appearing here are defined in terms of the spectral GFs discussed in Eq. (15), and the diffusion coefficient D According to Eq. (8) the modes f T and f L3 satisfy the diffusion equations where R T = 2∆Reg 01 , R L3 = 2∆Reg 31 describe the coupling of the quasiparticles to the superconducting condensate, and the collision integral for spin relaxation S L3 is given by Eq. (38). They are hence of the form of spectral current conservation equations, except for the source/sink terms provided by the collision integrals. The kinetic equations (72,73) are supplemented by the boundary conditions at the injector electrode given by the second and the third rows on the matrix of Eq. (57).
On the other hand, according to Eq.(67) the second pair of modes, f L , f T 3 , determines the energy and spin currents j e and j s , respectively. They satisfy the diffusion equations Here the spin-relaxation term S T 3 is given by Eq. (37). Typical energy dependencies of the kinetic coefficients in Eqs. (72,73,74) in the presence of spin splitting and spin-flip relaxation (β = 1) are shown in Fig. 12.
Charge and spin-energy modes. The two systems of diffusion equations (72-73,74-75) describe the coupled transport of spin, charge and heat in superconductors. We start by solving the system (72,73) describing the coupled charge and spin energy modes. The solution for f T and f L3 is a superposition of two exponentially decaying functions where the coefficients A i have to be determined from the boundary conditions. The inverse lengths k 1,2 are energy dependent. This dependence is shown in Figs. 13(a-d). Notice that both the charge and the spin-heat imbalance relaxation are nonvanishing for all energies, below and above the gap, due to the magnetic pair breaking effects [97,105]. In the absence of spin splitting, h = 0, R L3 = D L3 = 0, and the charge and spin-energy modes are decoupled [cf. Eqs. (72,73)]. Then k 1 and k 2 can be ascribed to the spin energy and charge imbalance relaxation, respectively. The non-zero Zeeman splitting h = 0 leads to the coupling of f T and f L3 , at low energies. At high energies ε ∆ though, when the superconducting correlations become negligible, the decoupled behavior is restored. This can be seen in the asymptotic behavior for ε → ∞ shown in Figs. 13(b,d). In this limit k 2 → 0 corresponding to a vanishing charge relaxation, whereas k 1 → λ −1 sn which is the spin-energy relaxation length in the normal state.

46
ZT (T ) in the presence of spin-orbital relaxation.     Spin and energy modes. We now analyze the system of equations (74)(75) that describes the coupling between the spin and energy modes, i.e., the components f L and f T 3 of the distribution function. The solution of the system can be written as the sum of two qualitatively different terms with coefficients B j determined by the boundary conditions. The first term in (77) describes a decay of the (spectral) spin imbalance with a characteristic length scale k s = S T 3 D L /(D 2 L − D 2 T 3 ). In the absence of a spin-splitting field and depairing mechanisms the expression for k se yields the energy dependent spin-relaxation length in superconductors obtained in Ref. [179]: This length is strongly renormalized by superconductivity if the spin-flip scattering rate is nonzero, β = −1. A more accurate calculation of k se (ε), including the modification of spectral functions in the superconductor due to the spin relaxation, is shown in Fig. 13(c). We see that k s is constant for β = −1 and has a pronounced peak near the self-consistent gap edge for β = 1. Note that the spectral gap is reduced by spin-flip scattering to the values smaller than ∆ 0 . Comparing Figs. 13(c,e) one can see that the relaxation lengths of charge imbalance and spin accumulation are very different even in the absence of a spin-splitting field. For example at subgap energies k s < k 2 . In contrast, for quasiparticles with energies ε ∆ the charge relaxation vanishes whereas the spin relaxation remains non-zero. Hence particle injection at voltages V inj well above the gap leads at large distances from the injector mainly to a charge imbalance.
It is worth noticing that in the presence of a non-zero spin-splitting field the behavior of k s (ε) is very similar for spin-orbit and spin-flip relaxation mechanisms as shown in Fig. 13(d).
The second term in the r.h.s. of Eq. (77). It describes an effective increase of the quasiparticle temperature associated with the f L nonequilibrium mode. This mode can only decay via inelastic scattering which is disregarded in the above analysis. In the following we assume the presence of an electrode at a distance L much larger than the spin relaxation length, and assume that the energy mode relaxes there.

Kinetic equations with supercurrent
The expressions for the current in Eq. (67) and hence the kinetic equations can be modified by driving a supercurrent through the superconducting wire. Its effect can be included by assuming a constant phase gradient ∇φ of the order parameter [184], i.e., ∆(r) = |∆|e iφ(r) . As a result, two new spectral coefficients appear: the spectral supercurrent j E and the spectral spin supercurrent j Es , defined as j E = DTr[τ 3 (ĝ R ∇ĝ R −ĝ A ∇ĝ A )]/(8∇φ) and j Es = DTr[τ 3 σ 3 (ĝ R ∇ĝ R −ĝ A ∇ĝ A )]/(8∇φ), respectively. In their presence all distribution function components are coupled so that the spectral currents appearing in Eqs. (72,73,74,75) Hence in this case the presence of spin-splitting and supercurrent leads to a coupling between all nonequilibrium modes. This is in contrast to the case without supercurrent when the matrix in Eq. (67) has a block-diagonal form so that the modes are coupled only pairwise.
In the absence of spin splitting, this leads to a coupling between charge and energy modes, and hence for example to the possibility of creating charge imbalance from temperature gradients [185,186,187]. For a spin-split superconductor, the supercurrent couples charge and spin modes [184] so that one can be converted into another one. In practice this conversion can be measured by analysing the different symmetry components of the non-local conductance with respect to the injection voltage and the spin polarization of the detector.
Strictly speaking, the assumption of a constant phase gradient together with the excitation of the nonequilibrium modes leads to a non-conserved total (energy integrated) charge current in the setup. Therefore, the non-linear problem of analysing the dynamics of modes for supercurrent needs to be coupled with the self-consistency equation (63), where the presence of the modes f T and f L3 affects the phase of the self-consistent order parameter and therefore the phase gradient.

Non-local conductance measurements
In a typical experiment quasiparticles are injected from a normal or a ferromagnetic metal electrode by applying a voltage V inj between the electrode and the superconductor. In this way the nonequilibrium states described by Eqs. (76,77) can be excited. We determine the coefficients A 1,2 and B 1,2 using the general boundary conditions (57) for the currents, where the voltagebiased normal electrodes are described by the distribution functions f  [38,31,162].
We plot the non-local conductance in three different cases (the three curves in Fig. 14). The blue curve is for a vanishing spin-splitting field and injector polarization. In that case the signal is solely due to charge imbalance. The black curve is with a non-zero spin-splitting field, but vanishing injector polarization, and in the red curve both the spin-splitting field and the injector polarization are non-zero. As discussed below, the former shows the thermally created spin 50 ZT (T ) in the presence of spin-orbital relaxation.
is the normal-state resistance of the wire with length ξ and cross section As.
accumulation, whereas the latter shows a combination of the thermally and electrically created spin accumulations.
In an experiment where the spin-splitting field is controlled by an external magnetic field, we also need to take into account the orbital depairing effect [Eq. (35)] of the magnetic field in addition to the Zeeman effect. Whereas the latter leads to the coupling of the energy and spin modes, the former affects especially the relaxation of charge imbalance. For the results in Fig. 14, we describe the relative strength of the orbital depairing and the spin-splitting field by setting (τ orb T c0 ) −1 = α orb (µ B B/T c0 ) 2 and choosing for the dimensionless quantity α orb a value close to the experiments in [31].
First, in the absence of a Zeeman field, N − = 0, and according to Eqs. (26,66) only the modes f T and f T 3 contribute to g nl . In such a case the contribution stemming from the spin accumulation is nonzero only if P inj = 0, which is the condition to obtain a nonvanishing f T 3 . However, this function decays over the spin diffusion length and therefore is negligibly small at the distances L det > λ sn = √ Dτ sn from the injector. Thus, the detected signal in this case is mostly determined by the charge imbalance µ. This explains the approximate symmetry with respect to the injecting voltage: g nl (V inj ) = g nl (−V inj ) (blue curves in Fig. 14). The charge imbalance contribution to g nl grows monotonically when |eV inj | > ∆ g . This behavior is determined by the increase of the charge relaxation scale at high energies, k −1 2 → ∞ shown in Fig. 13(c). In the presence of an applied magnetic field the charge relaxation is strongly enhanced due to the orbital depairing. As a result, the charge imbalance background signal is strongly suppressed by an increased h. On the other hand the spin imbalance contribution stemming from the f L mode is large. This term describes heat injection in the presence of a finite spin-splitting field h and has a long-range behavior. This contribution leads to the large peaks in g nl (V inj ) shown in Fig. 14. The non-linear heating produced by quasiparticles injected at voltages exceeding the energy gap explains the large electric signal observed in the experiments [38,31,162]. Notice that the peaks do not have exactly the same form so that g nl (V inj ) = −g nl (−V inj ). The deviation from the perfect antisymmetric form is due to the small but finite injector polarization P inj , which modifies the boundary conditions for the spin current.
One important feature shown in Fig. 14 is that the spin polarization peaks exist always in the presence of spin splitting, h = 0, even if the injector electrode is not ferromagnetic, P inj = 0. This explains the non-local conductance measurements of Ref. [163] by using a normal metal as injector.
The results summarized in Fig. 14 can have direct applications in the field of spintronics. On the one hand spin accumulation can be created without using ferromagnets as injectors. On the other hand, due to the coupling between the spin and energy mode, the spin signal can be controlled over a very long range. The only requirement for these two features to occur is to use a superconductor with a spin-split DOS due to a Zeeman or an exchange field.

Spin Hanle effect in normal metals and superconductors
In the previous section we study the spin injection and spin accumulation in a superconductor assuming that the magnetization of the injector and the applied field are collinear. We now lift this assumption and consider non-collinear field and magnetizations. This situation has been widely studied in the normal state. The non-collinearity between the external field and the spin of the injected electrons leads to a precession of the latter. This is the spin Hanle effect [101,165,166,169,188]. This precession can be measured via the non-local conductance in a multi-terminal sample as a function of the applied field.
In order to understand the spin Hanle effect in superconductors it is instructive to first discuss the spin precession in a normal metal within the formalism described above. We then demonstrate that in the superconducting state the spin precession can be either enhanced or suppressed, depending on the spin relaxation mechanism.

Normal-metal non-local spin valve
In the normal state the Usadel equation can be drastically simplified because the retarded and advanced Green's functions do not depend on energy and have the simple formĝ R(A) = ±τ 3 . Thus the reduced density of states equals to unity according to Eqs. (16)- (17). It follows from Eq. (13) that the Keldysh component is then proportional to the distribution matrix:ĝ K = 2τ 3f .
It is then useful to integrate Eq. (8) over energies and take the trace after multiplication with the Pauli matrices. This results in a diffusion equation for the nonequilibrium spin density S given by This is the spin diffusion equation [189,190,191] widely used in spintronics. Generalization of this equation in the presence of spin-orbit coupling has been done in several works [192,193,194,195]. Its right hand side describes the effective spin relaxation time defined as τ −1 sn = τ −1 so + τ −1 sf . The second term on the left hand side describes the torque induced by the external field via the spin splitting field. This torque leads to the spin Hanle effect, studied in the language of quasiclassics in Ref. [196].
Let us consider a lateral spin valve as the one sketched in Fig. 11(b). We assume that a field is applied in the x-direction whereas the injector and detectors are polarized in the z-direction. If the width and thickness of the N wire are smaller than the spin relaxation length the problem reduces to a quasi 1D geometry along the length of the wire. In such a case Eq. (79) leads to two linear coupled equations for the components S y,z , We have defined the spin diffusion length λ 2 sn = Dτ sn , the magnetic length The coefficients A and B have to be determined from the boundary condition, Eq. (48), which in this particular case has the simple form Here V I is the voltage across the injector/N interface and R I the resistance per unit area of the barrier. The resulting spin accumulation in the N wire is where V I and R I are the voltage drop and the resistance of the injector, respectively. A ferromagnetic voltage detector at a distance x = L from the injector detects the spin potential µ z = S z (L)/ν F . The resulting µ z (B) dependence, described by Eq. (80), coincides with the Hanle-shaped curves measured experimentally [165,166].

Hanle effect in spin-split superconductors
If the wire in Fig. 11 is in the superconducting state the non-local signal can change drastically as discussed in Ref. [180]. In contrast to the normal-metal case, the retarded and advanced GFs have non-trivial energy dependencies and therefore one cannot integrate straightforwardly the diffusion equation over the energy.
As in previous sections, we assume that the retarded and advanced GFs are position independent and the density of states is the one shown in Fig. 4. In contrast to the normal case the non-local Hanle signal in a superconducting wire depends on the dominating spin relaxation mechanism. We consider again the two mechanisms: magnetic impurities, Eq. (34), and extrinsic spin-orbit coupling, Eq. (33), and introduce again the parameter β describing the relative strength between them β = (τ so − τ sf )/(τ so + τ sf ).
In principle the Hanle signal in the superconductor would include contributions from the long-range energy mode discussed in Sec. 4.2 for the case of collinear field and magnetizations. To separate that contribution from the bare Hanle effect, we concentrate on the linear response regime of low injector and detector voltages, where the effect of the energy mode can be disregarded (cf. Fig. 14). In such a case the non-local conductance is exponentially small (∼ exp(−∆/T )) and therefore we focus on the non-local resistance V det /I inj , where V det is the voltage at the detector in the absence of a current and I inj is the current at the injector. To uncover the Hanle effect, we furthermore study the difference V S between the voltages for the parallel and anti-parallel orientations of the detector and injector polarizations. The non-local resistance of interest is thus R S = V S /I inj .
If we assume that the external field is applied in the z-direction and disregard the orbital effects, the retarded quasiclassical Green's function is described by only four nonzero components in Eq. (15). The energy dependence of these components is shown in Fig. 6.
Whereas the spectral terms remain constant in the S region, the spin-polarized current I inj that flows through the interface with the injector causes a spin accumulation characterized by the vector µ s . The current measured at the detector is then given by [180] where R det = R /A is the normal-state barrier resistance between the superconducting wire and the detector with cross sectional area A and spin polarization P det , and V det is the voltage measured at the detector with respect to the superconductor. Y = ∞ 0 dεN + ∂ ε f eq describes the effect of the density of states in Eq. (60). Its low-temperature analytical estimate for weak spin relaxation is given in Eq. (147a) below.
To find the voltage induced in the electrically open detector circuit we set I det = 0 in Eq. (81). The spin-dependent part of the non-local resistance R S = [V det (P det ) − V det (−P det )]/I inj can be found using the linear response relation for the injector current I inj = V inj Y /R inj with the normal-state resistance R inj of the injector: where V inj is the voltage applied at the injector in order to inject the current. We then need to determine the spin accumulation µ s (x) induced along the S wire which can be written in terms of the Keldysh component of the GF as with m(ε, x) = Tr(τ 3 Sg K )/8. For simplicity we assume that P inj ⊥ h so that the f L mode is not induced by the thermoelectric coupling at linear response, see Sec. 5.4.3. Therefore the spin accumulation does not contain the long-range contribution discussed in Sec. 4.2. The remaining part of m in Eq. (83) is given 55 by: where the vector f T = (f T 1 , f T 2 , f T 3 ) contains the spin components of the general distribution function (14). The first term in this expression is the usual quasiparticle contribution which also appears in the normal state. It is determined by the nonequilibrium distribution function at energies larger than the superconducting gap, that is where the total density of states is non-zero. In contrast, the second term in Eq. (84) only appears in the superconducting case and is finite for subgap energies due to the prefactor Img 33 (cf. Fig. 6). This term can be related to the difference N − between the DOS in spin subbands using the Kramers-Kronig relation for the retarded GF Img π . Let us assume that the spin-splitting field in the superconductor points in the z-direction, h = he z . Then the transverse components f T 1 , f T 2 satisfy the kinetic equations obtained from Eq. (8), In contrast to the normal case, Eq. (79), the diffusion coefficient is now a tensor with components that depend on energy, The terms on the right hand side of Eq. (85) are X 1 = (S T 1 − H 1 ), X 2 = (S T 2 + H 2 ), with S T 2 = 2τ −1 sn (Img 33 Reg 03 − βImg 01 Reg 31 ).
The H 1,2 terms are the "Hanle" terms that describe the coherent spin rotation and relaxation due to the action of the external field. In the normal case only H 2 and the first term of S T 1 are non-zero and we recover Eq. (79). In order to solve Eqs. (85a-85b) we need the boundary conditions at the interface with the injector. They can be obtained from Eq. (48). By keeping terms to leading order in the interface parameter κ I = 1/(R I σ N ) we obtain [180] These are evaluated at the position of the interface. Equations (85a,85b) can be re-written for the spectral density of spin polarization, Eq. (84), in a more familiar form, similar to the Bloch-Torrey transport equation [197] for the magnetic moment: Here γ = −2 is the electron gyromagnetic ratio and j sk is the spin current density in the k-th spatial direction. The diffusion coefficients are defined as Equations (92)(93) are written in the steady-state limit, which is enough for the description of the spin Hanle effect. In the dynamical case, the magnetization depends explicitly on time m = m(ε, t), and the l.h.s. of Eq. (92) contains additional terms. When the spin dynamics is slow enough, with a characteristic frequency much smaller than the superconducting energy gap, ω ∆, the Bloch-Torrey equation in the superconducting state can be written as 8 In Eqs. (92,96), τ S and h s are the transverse spin relaxation time and a correction to the effective Zeeman field appearing in the presence of spin-relaxation processes. The latter determines the (Larmor) precession frequency of the spins. Specifically, Both the transverse relaxation rate τ −1 S and the effective field shift h s = |h s | are proportional to the normal-state spin relaxation rate τ −1 sn . In the superconducting state τ −1 S and h s depend differently on energy for different spin-relaxation mechanisms. The typical dependencies are shown in Fig. 15 for β = ±1. One can see that τ −1 S has a step-wise behavior as a function of energy when the spinorbit relaxation dominates [ Fig. 15(a)]. In contrast, τ −1 S (ε) shows peaks near the spin-split gap edges for spin-flip relaxation [ Fig. 15(b)]. The effective field shift h s has different signs for β = ∓1 as shown in Fig. 15(c,d). Interestingly, the transverse spin relaxation τ −1 S vanishes at the subgap energies. Equations (85a-85b) imply that the transverse components of f T have the general form 8 For the general case of arbitrary large frequencies in the absence of gradients, see Eq. (108).   where A is an integration constant determined by the boundary conditions, and with Rek T > 0. In contrast to the normal metal case considered in Sec. 4.3.1, there is no straightforward relation between k T and the parameters τ S and h s . This is because the diffusion coefficients in the superconducting state entering Eq. (93) are energy dependent and hence the spectral spin current has a complicated expression that couples the different components of the spin polarization. The real part Rek T determines the spin relaxation length and the imaginary part Imk T gives the quasiparticle spin precession. In the superconducting state both these scales are energy dependent. In Fig. 16 we show how the precession and relaxation scales depend on energy and on the parameter β.
It follows that at low energies Re k T is larger than in the normal-metal case which corresponds to the limit ε ∆. At ε = 0 the characteristic relaxation length k −1 T (0) is determined by the coherence length ξ s = D/∆ 0 rather than by the normal-state spin relaxation length λ sn . In contrast, the imaginary part of k T vanishes at the energies below the spectral gap where the DOS shown by the dashed lined in Fig. 16 is absent, i.e., N + = 0. In the case of dominating spinorbit scattering (β < 0), the precession Imk T has a shallow peak at ε ≈ ∆ + h [see Fig. 16(a)]. If the spin-flip mechanism dominates (β > 0), this peak of Im k T is suppressed, Fig. 16(b).
Putting all this together allows us to evaluate the non-local resistance R S exhibiting the Hanle curves in the superconducting state as shown in Fig. 17.    the detector polarization P det = P det e y so that P det P inj ⊥ h. One can see that both, the precession and decay of the nonlocal signal, disappear at T → 0, whereas the shape of the curves at intermediate temperatures depends on the type of spin relaxation.
If instead of the above configuration we assume that the three vectors (P det , P inj , h) are perpendicular to each other (e.g., P det = P det e x , P inj = P inj e y , h = he z ), the subgap current is absent in the detector circuit and the corresponding spin signal R sx has a strong dependence on h even at the temperatures well below T c . This is shown in Figs. 17(c,d).
We are not aware of the measurement of the spin Hanle signals in the superconducting state. We emphasize that the above results are obtained in the linear response regime of small voltages. They disregard the thermal effects coming into play at higher voltages, especially at those of the order of the energy gap ∆±h. In the linear response regime the transport quantities in superconductors are exponentially suppressed due to the gap in the density of states, and hence the above theory is rather made for the thermally activated quasiparticles. Because of this, the best way to uncover the superconducting effects on the Hanle response would be to study it rather close but below the critical temperature Figure 18: Experimental setup from Ref. [198]. An applied RF field excites quasiparticle spin resonance in a thin Al film (S), which is indirectly probed by observing modification of the I(V ) curve of probe N1. Figure adapted from Ref. [198] (Creative Commons Attribution 4.0 International License).
To summarize: The kinetic equations for the nonequilibrium modes constitute a rather generic theoretical tool that we have used in this section to provide a quantitative description of some transport properties, such as spin injection and precession, occurring in superconducting structures in the presence of a spin-splitting field. The predictive capability of the derived kinetic equations has been proven in this section by contrasting the results with existing experiments especially in the case of collinear magnetizations. This makes these equations an ideal tool to study further effects that involve the coupling between different nonequilibrium modes in superconductors, as for example, the possible thermal spin Hanle effect and other non-linear effects taking place in systems with non-collinear magnetizations.
In section 5 we focus on thermoelectric effects that are also direct manifestations of the spin-splitting fields in the superconductor. However, in contrast to the section here where the focus is on the diffusion of the nonequilibrium modes within the superconducting wire, section 5 deals mainly with interface effects, and disregards the position dependence of the nonequilibrium modes within the superconductors. The next section discusses applications of the quasiclassical method to time-dependent nonequilibrium problems and for example generalizes the Bloch-Torrey equation (92) to large frequencies.
sectionNon-equilibrium quasiparticle dynamics The conduction electrons of a superconductor in an oscillating electromagnetic field can absorb energy via excitation of the existing quasiparticles, and by creation of new quasiparticles from breaking of Cooper pairs. In a stationary state, the excitation is balanced by corresponding relaxation processes. These are sensitive to the spectrum of states available in a spin-split superconductor.
Electrons couple to the electromagnetic field via the Zeeman and the orbital terms. The Zeeman coupling is the source for the conduction electron spin resonance (CESR), which was theoretically considered in early works [199,200,201,202,203]. The corresponding dynamical susceptibility for the spin-split case was discussed in Refs. [204,205,206]. The linewidth of the resonant absorption peak is determined by the spin relaxation time in the normal state, τ sn , and is not af-fected by the slow relaxation of the thermal long-range spin imbalance discussed in Sec. 4.2. The magnetic screening complicates the experimental observation of the resonance. Results however exist in the mixed state of type II superconductors as reported in Refs. [207,208,209], and more recently, also in spin-split Al thin films in Ref. [198] [see Fig. 18]. For the Al films, linewidths of 100 ps were observed, consistent with the expected magnitude of τ sn in aluminum.
The CESR physics is also related to pumping effects, where an externally driven spin precession drives a nonequilibrium state or currents in the superconductor, or currents across a junction. Effects of this type have been considered in different ferromagnet/superconductor structures in several works, in adiabatic [210,211] and nonadiabatic cases [212,213,214,215].
As mentioned above, a rf field also couples to the orbital degree of freedom of the electrons. This coupling drives ac currents that also excite quasiparticles. In the absence of spin splitting, the linear and nonlinear response of superconductors has been extensively studied [17]. In the linear response regime, the complex impedance for spin-split superconducting films was considered theoretically and experimentally in Ref. [216] adding Zeeman energy shifts to BCS theory results [217,218]. In the nonlinear regime, strong driving can modify the observed quasiparticle spectrum [219]. Moreover, at temperatures close to T c , exciting quasiparticles can lead to an increase of the superconducting gap ∆ [220,221,222]. In spin-split superconductors, this effect is modified by the spin-split density of states and spin-flip scattering [123].
In this section we discuss the dynamic response of spin-split superconductors in terms of the time-dependent quasiclassical equations, both in the linear and nonlinear regimes. In addition to these effects, oscillating fields can excite magnetic impurities and nuclear moments. Discussion of such extrinsic resonance effects can be found in the articles [223,224,225,226], and are beyond the scope of this review. We also do not discuss magnetization dynamics and how superconductivity affects it; reviews on this active topic can be found in other works [227,228].
The kinetic equation for the distribution functionf follows from Eq. (101) and the above assumptions, whereĵ R/A ≡ Dĝ R/A •∇ĝ R/A and the collision integral is defined aŝ andẐ Note that the definition of collision integrals in Eq. (103) 92)). Such equations can be derived from the Usadel equation in linear response to h 1 , and from other approaches [204,205,206,230]. Here we briefly discuss the physics within the Usadel framework.
It is useful to separate the corrections to the Keldysh Green's function in Eq. (101) due to the distribution function from those due to the modification of the spectral functions: δĝ K = δĝ K reg + δĝ K an , defining δĝ K reg = δĝ R • f eq − f eq • δĝ A . We split δX K similarly. The equation [X • ,ǧ] K = 0 then reduces in linear response toX The solution to Eq. (106) can be obtained with the Ansatz δĝ K To determine the solution components, we can take the trace 1 8 Tr[(. . .)σ] of the (ε, ε − ω) frequency component of Eq. (106). This results in where m = 1 8 Tr[τ 3 σδĝ K an (ε, ε − ω)] describes the nonequilibrium magnetization (spin accumulation), m = 1 8i Tr[τ 1 σδĝ K an (ε, ε − ω)] describes a correction to spin scattering by superconductivity, and N + = [g 03 (ε) + g 03 (ε − ω) * ]/2, F + = [g 01 (ε) − g 01 (ε − ω) * ]/(2i) are finite-frequency generalizations of the spinaveraged density of states and anomalous functions, respectively. The term ] describes the exciting field. There is redundancy in the representation, which enables writing m in terms of m, resulting in the final Bloch equation. Analogously to the spin Hanle effect discussed in Sec. 4.3.2, the spin relaxation time τ S is renormalized with respect to the one in the normal state, τ sn , and there is also a correction h s to the Zeeman field: where N z = [g 33 (ε)−g 33 (ε−ω) * ]/(2i) and F z = [g 31 (ε)+g 31 (ε−ω) * ]/2. For slow driving, ω → 0, the spin relaxation time τ S and the Zeeman field correction h s coincide with those visible in the Hanle effect (cf. Eq. (97) and Fig. 15), which also involves precession of the transverse spin component. The result (108) describes resonant excitation of the transversal modes, f T j , which correspond to a nonequilibrium contribution to the spin accumulation (22), where δµ reg s arises from the modification of the spectral functions, δg R/A . The result contains the conduction electron spin resonance peak, δµ s ∝ h 1 (ω)/[4(h+ h s ) 2 − (ω + i/τ S ) 2 ] at frequency ω 2|h 0 |. As discussed in Sec. 4.2, the f T j modes can relax due to elastic spin-flip scattering, which determines the peak absorption linewidth 1 τ S ∝ τ −1 sn . The result in Ref. [203] can be obtained from Eq. (108) in the quasiequilibrium approximation δf (∂ E f eq )µ s · σ and in the limit h 0 ∆. Equation (108) also coincides with the result in Ref. [204]. If we focus on the orbital effects, the complex impedance of the superconductor, i.e., the linear response to an oscillating electric field described by A(ω), can also be obtained within the Usadel framework. In the above chosen (London) gauge, the perturbation to time-averagedǧ is of the order D|A(ω)| 2 /ω [220] and can be neglected in linear response at nonzero frequency. The charge current response is then given by where and g ↑/↓,1 = (g 01 ±g 31 )/2, N ↑/↓ = Re[g 03 ±g 33 ]/2. The prefactors in Eqs. (112)- (114) are given by the normal-state conductivity, σ N = 2e 2 ν F D, and the kernels R σ have the BCS form [217,218], and are decoupled for each spin species. In other words, the electric field couples to the orbital motion of the electrons and, in our approach, it conserves spin. Therefore, the dissipative response is zero at frequencies ω < 2∆ in the absence of spin flipping. Phenomenological modification of the above by allowing direct spin flips in the coupling matrix element was considered in Ref. [216], obtaining a reduction in the dissipative pair-breaking threshold frequency from 2∆ to 2∆ − 2h (c.f. Fig. 19). Note however that the derivation of the Usadel diffusion equation discussed above works in leading order in τ el and neglects spin-orbit effects when dealing with the vector potential.

Nonlinear spin imbalance of a spin-split superconductor in a rf field
We now go beyond the linear response regime and focus on the stationary quasiparticle distribution which enters the observables. This is determined by disregarding the gradients, in which case the collision term in Eq. (103) is equal to the contribution from the time-dependent vector potential. This results in This rate equation, similar to that used in Ref. [122], describes the balance between excitation of quasiparticles induced by ac fields (Î ac ), the spin relaxation by the elastic spin-flip scattering (Î sn ), and the inelastic relaxation (Î in ). Each term in Eq. (117) can be described within the quasiclassical approach. The orbital ac term of the collision integral can be obtained by assuming a time-dependent spatially uniform vector potential A(t) = A 0 cos(ωt). Here, it is useful to simplify the problem by expanding in DA 2 0 /ω 1, and reduce it to a form that only involves the dc component of the distribution function. In this

case [220]
for the I ↑/↓ = Tr[(1 ± σ z )Î/2] components. We assume that all spin dependent fields are collinear and hence we can define f σ=↑/↓ = f L ± f T 3 , cf. Eq. (14). Notice that the charged modes f T , f L3 are not excited by a uniform A(t).
The absorption kernel R is the one appearing in the real part of the conductivity, Eq. (115). Qualitatively, this collision integral results in the depletion (accumulation) of quasiparticles in an energy band of width ω above (below) the gap edges |E| > |∆ ± h| (see Fig. 19). At higher energies the electron distribution corresponds to an increased temperature.
The collision integrals for the relaxation processes can be written as, cf. Eq. (37), where we describe inelastic scattering in a relaxation time approximation. The spin-conserving microwave absorption term does not directly generate spin imbalance in the superconductor. However, as illustrated in Fig. 19, spin relaxation converts the accumulation of quasiparticles above the gap to an imbalance in the number of spin-up and spin-down quasiparticles. From the kinetic equation (117) we obtain in the limit 1/τ sf 1/τ in , The first term contributes zero total spin imbalance (26), but the second gives a nonzero contribution of the order of τ in /τ sn . A related effect in quasiparticle injection was discussed in Ref. [122]. The generated spin imbalance can be in principle measured by a spinpolarized probe junction, as discussed in Sec. 4.1. Note, however, that in this type of experiments other processes may need to be also considered. On the one hand, photoelectric signals can also occur in the absence of spin splitting as shown in Refs. [231,232]. These processes, however, generally scale with ∝ τ el rather than ∝ τ in . On the other hand an ac bias over the S/F junction may also result in rectification, as discussed in the next section.
It is interesting to notice that the nonequilibrium quasiparticle accumulation generated by the microwave absorption affects the magnitude of the the superconducting gap ∆, leading to deviations from the results in Ref. [220]. Results for the self-consistent ∆(T ) from numerical calculations [123] are shown in Fig. 20 for h = 0 and h > 0. In the absence of spin-splitting, h = 0, a gap enhancement [220] occurs. For h > ω, an additional instability develops at ∆(T ) = h, corresponding to a coexistence of two solutions where the spinaveraged DOS is either gapless (|∆| < h) or gapped (|∆| > h). The instability requires the presence of strong spin-flip scattering: for τ sn τ in the exchange field does not cause significant qualitative changes, because in that case the quasiparticle accumulations above the gap edges that generate the enhancement, are approximatively described as two independent copies of the effect at h = 0.

S/F tunnel junction dynamics
For adiabatic excitation, ω ∆ dynamics of S/F tunnel junctions can be described by the dc relations discussed in Sec. 3 and 5. That this limit can be reached at time scales shorter than the spin-relaxation and inelastic times, was used in a recent experiment [233] to probe the relaxation dynamics in spin-split thin-film superconductors.
At higher frequencies, photoassisted tunneling breaks the adiabatic description. This can be taken into account via a standard tunneling Hamiltonian approach [234,235,236]. The result for the tunneling current is where T is the tunneling matrix element, G >(<) = 1 2 G K ± 1 2 (G R −G A ) are (nonquasiclassical) Green functions for the superconducting (k) and non-superconducting (q) sides, and the phase φ(t) = e t dt V (t ) is related to the time-dependent voltage V (t) across the junction. The above result applies for spin-conserving tunneling with collinear magnetizations. If the terminals are in an internal equilibrium (no spin or charge imbalance), changing integration variables yields the time-dependent generalization of Eq. (3), where I dc is given by Eq. (3), and the memory kernel is [237] and J 1 is a Bessel function. The resulting current-voltage relation is asymmetric, I dc (−V ) = − I dc (V ) for h = 0. This implies that such junctions rectify ac signals [233]. For an ac signal V (t) = V dc + V ac cos(ωt), from the above we have the average dc current [238] At V dc = 0 and small V ac at uniform temperature T h, ω, the rectified current is I dc ∝ P F V 2 ac h/(R T ∆ 2 ), [123] proportional to the exchange field in S. Measurements attempting to use F probes for the nonequilibrium ac effects discussed in the previous sections need to take into account such rectification, or try to suppress V ac e.g. via a suitable microwave circuit design [239].
While dynamical effects in superconductors have been studied extensively experimentally, we are aware of relatively few studies on spin-split superconductors. Tunnel junction rectification effects were observed for example in Ref. [233]. However, to our knowledge, higher-frequency experiments probing gap enhancement or photoassisted tunneling have not been reported so far.

Thermoelectric effects in superconducting structures
A temperature difference across an electric contact typically leads to heat currents aiming to relax this difference. In some cases it may also lead to observable charge currents. Reciprocally, an electrical voltage may drive a heat current not only in nonlinear response (due to Joule heating) but also for small voltages. This connection of charge and heat currents is called the thermoelectric effect. The traditional view of thermoelectric effects in superconductors is that if they exist, they must be very weak. In bulk superconductors, this is partially because any thermoelectrically generated quasiparticle current is screened by a supercurrent [240]. Ginzburg suggested [241] to measure this supercurrent by using an additional constraint to the phase of the superconducting order parameter in a multiply connected structure. However, even this thermally created phase gradient tends to be weak, owing to the near-complete electronhole symmetry in superconductors [242].
Thermoelectric effects typically require strongly energy dependent density of states of the charge carriers. Such energy dependence is present in the BCS density of states of the superconductors. However, typically this density of states is quite symmetric with respect to the Fermi level, and therefore any contribution of the positive energy excitations ("electrons") on thermoelectric effects is cancelled by the negative energy excitations ("holes"). Breaking this electronhole symmetry would hence allow for the appearance of strong thermoelectric effects. As we discuss in this section, this is what happens in spin-split superconductors, as an exchange field breaks the symmetry in each spin sector, but so that the overall spin-summed energy spectrum remains electron-hole symmetric. Further breaking the spin symmetry in transport through a spin filter can then provide large thermoelectric effects as the two spins are weighed differently [243,29,244]. This prediction has been recently confirmed experimentally in [177,245].
In this section we give an overview of the different types of thermoelectric effects discussed for superconductors, and then concentrate on the ones obtained in superconductors with a spin-splitting field. We discuss both the linear response regime and beyond it, and consider also the limiting features such as the electron-phonon coupling. We show that under suitable conditions, in particular for close-to-optimal spin filters, the efficiency of thermoelectric conversion can become very large and exceed that obtained for best thermoelectric devices operated at or above room temperature. Besides the regular quasiparticle current, in devices coupling two superconductors, one with a spin-splitting field and one without, the thermoelectric effect can also be converted to a phase gradient [30]. This large thermophase effect is discussed in more detail in Ref. [3]. As these effects require low operating temperatures, they obviously cannot be directly used to improve the efficiency of various everyday devices. However, the strong thermoelectric effect may become relevant in other types of applications, such as sensors, where the measured (wide-band) signal consists of heating one part of the system [33]. The thermoelectric conversion can then be used to convert the resulting temperature difference to a charge current.

Thermoelectric effects and heat engines
Biasing a contact with a small voltage V and a small temperature difference ∆T leads to linear-response charge and heat currents I andQ of the form The 2 × 2 conductance matrix in Eq. (127) is a part of the generalized 4 × 4 Onsager matrix (59) connecting different interface currents and potentials as discussed in Sec. 3.6. This matrix is symmetric in the presence of time-reversal symmetry, and also in all particular cases considered in this section. Here we consider a conventional situation when the spin-dependent potentials V s and T s are negligibly small. This is relevant when the spin relaxation rate in the electrodes around the contact exceeds the tunneling rate across the contact. This approximation allows the reduction of the general boundary condition (59) to the simpler one (127) describing thermoelectric response. Whenα is not too small, electrical energy may be converted to heat or cooling (Peltier effect), or reciprocally a temperature difference may be converted to electrical power (Seebeck effect). The efficiency of this conversion can be described by constructing a model for a generic electrical heat engine (see Fig. 21). There, a load with resistance R L is driven by the power drawn from the thermoelectric element across which there is a temperature difference ∆T . The power dissipated on the load is P = IV = R L I 2 , whereas the voltage across the thermoelectric element is −IR L . Plugging this into Eq. (127) yields The extracted power hence is On the other hand, the thermoelectric element extracts heat from the temperature difference, i.e., trying to balance it, with the poweṙ The efficiency of thermoelectric conversion is hence where N =α 2 /(G th GT ) ≤ 1. Alternatively, we can write N in terms of the usual thermoelectric figure of merit, where S =α/(GT ) is the thermopower (Seebeck coefficient) andG th = G th − α 2 /(GT ) is the thermal conductance at a vanishing current. Now we should choose R L to optimize the device. For example, the maximum efficiency is obtained with This result is consistent with that obtained in Ref. [246] in the linear response limit T cold /T hot ≈ 1. On the other hand, optimizing the device to yield a maximum power output requires GR L = 1, corresponding to the limit [247,248] Here η CA ≡ 1 − T cold /T hot ≈ ∆T /(2T ) is the maximum efficiency obtained when ZT → ∞. Both of these efficiencies are maximized when the thermoelectric figure of merit becomes large. Above room temperature, the highest figures of merit are obtained in certain strongly doped semiconductor structures [249,250]. The record values are of the order of ZT 1 . . . 2. The particular value and the optimal temperature where it is obtained results from a competition of two generic temperature dependencies: that of phonon heat conductance, and that of the (typically activated) process yielding the thermopower. Both of these decrease towards low temperatures, but the previous decreases as a power law, whereas the latter decreases exponentially. In an electron system the heat conductance is typically close to the Wiedemann-Franz limit G el th ≈ LGT , where L = π 2 k 2 B /(3e 2 ). The total heat conductance is obtained from the sum of the electronic and phononic contributions, G th = G el th + G ph th . Since ZT ∝ GT /G th = 1/[L + G ph th /(GT )], the only way to improve this is to minimize the phononic contribution. Generally the latter decreases towards low temperatures faster than linearly. However, the thermopower is typically an exponential function of temperature, S ∝ exp(−∆/k B T ) as the best thermoelectrics contain a gapped dispersion with the gap ∆. Therefore, the optimal ZT takes place at a temperature which is some fraction of ∆. As shown in Sec. 5.4 (see in particular Fig. 28), the same is true for a normal metal island coupled to spin-splitted superconductors. In that case, the magnitude of ∆ is just orders of magnitude lower than in semiconductors structures.

Thermoelectric effects in superconductors
Research on thermoelectric effects in superconductors dates back to the 1920's, when Meissner [240] concluded them to be absent because any thermoelectric current is cancelled by a counterflowing supercurrent. Ginzburg [241] showed that this is no longer the case in multiply connected bimetallic superconducting structures. This is because the presence of supercurrent is linked to the gradient of the phase of the order parameter, and in multiply connected structures uniqueness of the order parameter imposes a constraint on the relation between the flux through the ring and the circulating supercurrent. This situation can be visualized as in Fig. 22. For definiteness, let us consider a bimetallic superconducting loop, where the thermoelectric effects mostly take place at the two contacts. The current through contact i = 1, 2 consists of a sum of the supercurrent and the thermoelectric current, i.e., In a closed circuit this current must equal for both junctions, and thereby it yields a relation between the phases ϕ 1,2 . A second relation fixing the phases is obtained in the presence of a flux Φ through the loop with inductance L = L T + L B . Without loss of generality the total loop inductance can be described with a single quantity. Namely, the phases must be fixed so that they minimize the energy where n ∈ Z, Φ is the magnetic flux through the junction and Φ 0 = h/(2e) is the flux quantum. In the following we denote Φ/Φ 0 = f + 2π(m − n) with f ∈ [0, 2π[ and m ∈ Z. For f 1 and I th,i I c,i , the induced phases are small, and we may linearize the current-phase relations. As a result, we obtain a circulating current where L tot = L+ (I −1 c1 +I −1 c2 )/2e is the total inductance of the superconducting loop.
Let us assume heating the lower superconductor of the bimetallic loop so that the temperature increases from T to T + ∆T . Such a heating induces thermoelectric currents I th,1 =α 1 ∆T /T, I th,2 = −α 2 ∆T /T.
This produces a circulating thermoelectric current It can be measured for example by placing a SQUID on top of the bimetallic loop, and measuring the induced flux Φ ind = M I circ , where M is the mutual inductance between the two systems. The size of the coefficientα was calculated in Ref. [242] to bẽ where the latter form is due to the reduction of the quasiparticle density in the superconducting state, andα N is the size of the thermoelectric coefficient in the normal state. The precise value ofα N depends on the exact electronic spectrum of the metals in question and its calculation needs extending the theory beyond the quasiclassical approximation employed in this review, as within that approximationα N = 0. For a quadratic dispersionα N = , where E F is the Fermi energy. At temperatures T ∆/k B ,α is thus a product of two small coefficients,α N ∝ k B T /E F , and G(∆/T ). Such a smallα is not easy to measure quantitatively. Because of many spurious effects in such measurements, it is not simple to make the experiments agree quantitatively with this theory.
Nevertheless, at least close to the critical temperature the thermoelectric flux should be observable, and the first experiments to measure it were done in the 1970s and early 1980s. One set of experiments [251,252] followed the idea of Ref. [253] and measured the thermoelectric voltage across a superconductornormal-metal-superconductor contact as the thermoelectric current exceeded its Josephson critical current. The results of these experiments are in line with the expected magnitude of the thermoelectric signal. On the other hand, Ref. [254] measured the thermoelectrically generated flux in a bimetallic loop formed from superconducting Pb and Ti, close to the critical temperature of the latter. These experiments are also discussed in Ref. [255]. Surprisingly, the experiments demonstrated fluxes five orders of magnitude larger than predicted by theory. This discrepancy annoyed Ginzburg so much that he devoted an entire chapter on the topic in his Nobel colloquium [256].
Recently, Ref. [257] claim to have solved this discrepancy with new experiments performed on much smaller superconducting loops than what was possible in the early 1980s. According to them, the discrepancy originated from the temperature dependence of the inductances L tot and M as well as the flux f (via the effective area of the loop), and hence they change as one end of the bimetallic loop is being heated. This produces additional contributions to the thermoelectric flux, The additional contributions are in practice much larger than the pure thermoelectric effect. According to Ref. [257], whereas the effects from f can be accounted for by measuring the period of oscillations as the external flux is altered, the effects due to the trapped flux, accounted for by m, is much harder to deduce, and was likely the reason for the discrepancy, as in the early measurements a geomagnetic field amounted to m ∼ 10 6 . Using much smaller loops, Ref. [257] were able to control the number of trapped flux quanta, and hence get rid of the spurious effects. The remaining thermoelectric flux that they observed is more or less in accord with the value obtained from the above theory. However, the temperature dependence of the measured flux depends, besides G(x) above, also on the temperature dependent inductance, and on the temperature dependence of the main heat contact between the heated electrons and the phonon bath, so the measurement could not directly deduce G(x). The conclusion from these theory and experimental works is that in conventional superconductors thermoelectric effects can be nonzero, but they are extremely weak, and therefore difficult to access. However, other types of superconducting heterostructures besides the bimetallic loop do contain relatively strong thermoelectric effects. In particular, a supercurrent flowing along a temperature gradient leads to an appearance of charge imbalance [187,258,185], which can be measured [186] for example via a non-local geometry similar to those discussed in Sec. 4. A similar type of an effect was found in Andreev interferometers [259,260,261] composed of a hybrid multiterminal geometry of a normal metal in contact with a superconducting loop. These effects were considered theoretically in Refs. [262,263,264,265]. Similar effects are also predicted in ballistic systems as discussed in Refs. [266,267]. Nevertheless, these effects either require complicated multiterminal geometries, and otherwise the thermoelectric effects are typically very weak.
In what follows we show how this situation can be completely reversed in the presence of the spin-splitting field, provided we add a second ingredient into the theory: spin filtering. In this case such superconductor/ferromagnet hybrids can become almost ideal thermoelectric devices. In particular, Ref. [29,244] showed how three-terminal proximity-coupled superconductor-ferromagnet devices can show non-local thermoelectric effects: in this case the density of states (DOS) in a normal metal coupled both to a superconductor and a ferromagnet becomes spin dependent, and the spin-resolved DOS is also electron-hole asymmetric, resulting in the strong thermoelectric effect. On the other hand, Ref. [268,269] showed how a metallic bilayer consisting of two superconductors or a superconductor and a normal metal, separated by a spin-active interface (i.e., an interface whose transmission properties are characterized by a spin-dependent scattering matrix), can exhibit large thermoelectric response. Around the same time, Ref. [243] showed how a superconductor with a spin-splitting field, tunnel coupled to a normal metal via a spin-polarized interface, exhibits a thermoelectric effect where the figure of merit can become very large. To our knowledge, only this mechanism has been so far accessed experimentally [177,245]. We explain this mechanism in more detail below.
Another way to affect the thermoelectric response via magnetism was discussed in Ref. [232] and later in Ref. [270], who argued how magnetic impurities inside a superconductor enhance the thermoelectric coefficient by a large factor k F compared toα in Eq. (140). Here is the elastic mean free path and k F is the Fermi wavenumber. This effect results from the electron-hole asymmetric Andreev states forming in the vicinity of the magnetic impurities.
The above discussion concerns metallic structures. In semiconductor quantum dots the thermoelectric effects can be large when the spectrum of the quantum dot is electron-hole asymmetric [271], even without the presence of superconductivity or magnetism. Theoretical studies showed that the combination of the latter affects the symmetries and may increase the thermoelectric response in quantum dots coupled to superconducting and magnetic electrodes [272,273]. 5.3. Thermoelectric effects at a spin-polarized interface to a spin-split superconductor Tunneling into superconductors has been long used to probe the superconducting density of states [58], for thermometry [128], and for measuring the nonequilibrium distribution functions [274]. We show below that a spinpolarized tunnel contact to a superconductor with a spin-splitting field also ex-hibits a giant thermoelectric effect. Let us denote the spin-dependent normalstate conductance of the tunnel contact by G ↑/↓ for spin ↑ / ↓. In this case the standard tunneling theory yields the spin-resolved charge and heat currents across the tunnel contact from reservoir R to the spin-split superconductor The formula for the current is thus the spin-resolved version of Eq. (3). Here functions of the reservoirs biased at potentials µ L/R and temperatures T L/R . We assume the reduced density of states N R (ε) (i.e., total density of states at energy ε divided by the one at the Fermi energy) of reservoir R spin independent for simplicity; the possible splitting of the Fermi sea in electrode R is included in G σ . The reduced density of states in the superconductor for spin σ is N σ (ε). The heat currentQ i σ is calculated separately for i = R, S, using the potential µ R/S , because the two heat currents differ by the Joule power I(µ R − µ S )/e. Let us denote the spin-dependent reduced density of states via N + = N ↑ + N ↓ and N − = N ↑ − N ↓ and define the total tunneling conductance G T = G ↑ + G ↓ and spin polarization P = (G ↑ − G ↓ )/G T . In that case the spin-averaged tunnel currents In the following, we first analyze these currents in detail, and then discuss a method of building a near-optimal heat engine based on such junctions. In particular, we use the densities of states described in Sec. 3. We first disregard the spin relaxation effects on the density of states, because this assumption allows for some analytically treatable limits and because it is a fair approximation for example for Al. This assumption is lifted in Fig. 31.

Heat current beyond linear response
Let us consider first the case when the reservoir R is a normal metal and therefore N R (ε) = 1. The heat currentQ describing the cooling of the normal metal is positive at voltages eV = (µ R − µ S ) ≈ ∆ even without spin splitting, i.e., the work done by the voltage source can be used to cool the normal metal [275,276,277]. However, that heat current is even in the voltage, and therefore it does not result from the usual Peltier effect (Eq. (127) forQ) where the cooling power is linear in voltage, and therefore can be reversed by reversing the sign of the voltage. As we show below, in the presence of spin polarization eV =" P and with a non-zero spin-splitting field h in the superconductor, the cooling power obtains also components that are odd in voltage, in analogy with the usual Peltier coolers. This leads to an improved coefficient of performance of electron refrigeration [39,278]. The cooling power from reservoir R as a function of voltage is shown in Fig. 23 with and without an exchange field, assuming the ideal case of unit spin polarization P = 1. The figure also shows how decreasing the normalmetal temperature lowers the cooling power. In the absence of extra heating or energy relaxation processes, the normal-metal temperature at the given voltage would then be fixed to the value nullifying the cooling power.
In an electron refrigerator, the cooled element is an island coupled to two electrodes. Let us consider a normal-metal island playing the role of the reservoir R, and coupled to two spin-split superconductors via ferromagnetic insulators with polarizations P L and P R , respectively. The optimal situation is realized when P L = −P R = ±1 or when the exchange fields in the two superconductors are reversed. The total cooling power from the island isQ N = Q R (V /2, P L ) +Q R (−V /2, P R ). It works against other relaxation mechanisms, so that the stationary temperature T of reservoir R is determined from heat balance [128],Q Here we assume that the dominant heat relaxation mechanism on the reservoir R with volume Ω N is due to electron-phonon interaction with strength Σ N , described by Eq. (45). In addition, we assume that the spin accumulation on the island, produced by the nonequilibrium driving, is negligibly small due to eV =" spin relaxation. The assumptions relevant for this limit are discussed in Sec. 5.5 below. The resulting island temperatures are shown in Fig. 24 for different parameters of the S-FI-N-FI-S junction. We have chosen the parameters so that for h = P = 0 they correspond to a typical T (V ) curve found experimentally (see for example [276]). In the absence of spin filtering, P = 0, temperature obtains a minimum at eV ≈ 2(∆ − h), and T (V ) is a symmetric function of voltage. Generally, spin splitting in this case makes the cooling worse, so that the minimum reached temperature is higher than in its absence.
The behavior of T (V ) changes in the presence of spin filtering, P = 0. In particular, the curve becomes non-symmetric, and there is cooling even in the linear response regime, i.e., low voltages. This Peltier effect is discussed more below. In addition, the minimum reached electron temperature is generally lowered by an increasing P , and for a large P ≈ 1, the lowest temperature may be obtained at a non-zero exchange field. However, this effect seems rather weak for the considered parameters.
One possibly relevant aspect of such a magnetic cooler is the fact that the optimum temperature is obtained at a lower absolute value of the voltage [39]. This translates into a somewhat lowered Joule power injected to the device. As this heat is dumped into the superconductor, at the lowest temperatures the heating of the superconductor becomes the dominant limiting obstacle for cooling instead of the electron-phonon coupling [279,280]. The consequences of this are analyzed in [278].
Besides the normal metal, the magnetic element can be used to refrigerate eV =" the superconductor [278,281]. We illustrate this by plotting the cooling power of the superconductor in the case of an N-FI-S junction in Fig. 25. Without spin splitting and a spin-polarized interface, the cooling power is always negative, i.e., corresponding to heating. However, since the Peltier cooling with a non-zero P is linear in voltage, it has to result to a non-vanishing cooling power also from the superconductor. The maximum cooling power is less than in the case of cooling the normal metal (Fig. 23). However, as indicated in Fig. 26, the resulting temperature drop is nevertheless appreciable because the energy gap also weakens the electron-phonon coupling (see Eq. (46)).

Linear response and a heat engine
In the presence of non-vanishing spin polarization P and spin-splitting field h, the heat current has a linear component in voltage V . This Peltier effect is visible in the h = 0 curves in Figs. 23 and 25. In what follows, let us concentrate on low voltages and small temperature differences eV, k B ∆T ∆, so that the response of the junction can be described by the linear-response scheme as in Eq. (127). At low temperatures k B T ∆ − h the thermoelectric coefficients can be evaluated analytically and they are [243] In an open-circuit configuration (I = 0) one rather measures the thermopower At low temperatures the thermopower can thus become very large, and within the above scheme it would even diverge for k B T → 0. However, in practice this divergence would be cut off via either circuit effects (at T → 0 the conductance would also tend to zero) or for example the spin relaxation neglected above. Nevertheless, with proper circuit design one should be able to measure a thermopower much exceeding k B /e in this setup. The prediction for the strong thermoelectric effect was confirmed by a recent experiment by [177] via the measurement of the thermoelectric current I = −α∆T /T . This thermoelectric effect also provides a partial explanation [35] for the measurements of the long-range non-local spin signal [31,38]. In what Figure 27: Two realizations of a superconductor/ferromagnet heat engine, where a heating power P heat leading to a temperature increase in the absorber can be converted to electrical power P work and dissipated in the "device" connected to the heat engine. In both cases the setup requires two junctions with antiparallel magnetization directions. The ferromagnetic insulator FI has a magnetization direction along with the two spin-polarized contacts. To disregard spin accumulation, the island has to be large compared to the spin relaxation length. a) Heat engine with a normal-metal heat absorber. b) Heat engine with a spin-split superconducting island. The ferromagnets can also be replaced by a normal metal if the interfaces to the superconductor contain a ferromagnetic insulator.
follows, we consider the possibility of realizing a true heat engine with a large figure of merit. We discuss this in the case where either a normal metal or the spin-split superconductor realizes an island and therefore works as the heat absorber, see Fig. 27.
Let us first assume that it is only the electrons of the ferromagnetic island that are heated, instead of both the electrons and phonons. 9 In this case the figure of merit of the heat engine is limited by the spurious heat conduction due to the electron-phonon coupling, with heat conductance G e−ph = 5ΣΩT 4 [see Eq. (45)], and from the heat conductance to the extra superconducting electrodes S, with an exponentially suppressed heat conductance, provided the heating current is small. In this case the heat engine has an efficiency described by a figure of merit that for k B T h and P = 1 tends to is a dimensionless quantity characterizing the strength of electron-phonon coupling. Here we take into account the fact that the thermoelectric effect takes place across both contacts, and hence all quantities except the spurious heat conduction channel are doubled. If we do not describe half-metals, i.e., 1 − P 2 , we can simplify this more by dropping the last term in the denominator. We plot this as a function of temperature in Fig. 28 for some representative parameters. For example, for Ω = 0.005 µm 3 , Σ = 10 9 W µm −3 K −5 and 1/G T = 3 kΩ, g = 100. The figure of merit exhibits a peak at around k B T ≈ 0.1∆, and its peak value depends strongly on the relative strength of the spurious heat conduction channels. As seen in the figure, reaching high values for the figure of merit is quite challenging in this way, as besides a high polarization close to one, it requires quite low normal-state resistance of the junctions, a combination not very easy to reach with ferromagnetic insulators.
Note that often such spurious heat conduction mechanism that limit the highest available ZT are disregarded from the theoretical analysis of the figure of merit, for example in the case of quantum dots [272].
Another possible setting for the heat engine is the one where the island is the spin-split superconductor, and it is connected to normal-metal or ferromagnetic electrodes via spin filters as in Fig. 27(b). In the linear response regime this is otherwise similar to the previous case, except that the electron-phonon coupling inside the island has an exponentially suppressed heat conductance, Eq. (46). In this case, as long as the pure tunneling limit remains valid, and the island volume is not overly large, the low-temperature behavior of ZT is dictated mostly by the deviation of the polarization from unity as shown in Fig. 29. However, because of the exponentially decaying heat conductance, also in this case other spurious processes besides electron-phonon coupling start to limit ZT , making it vanish at T → 0. In the case of Fig. 29, the simplest model for them is due to the small but nonvanishing density of states assumed for the superconductor [243], and described by the Dynes parameter Γ [282]. 10 When the heat input into the heat engine is due to electromagnetic radiation coupled to the device, the conversion of this input power to direct charge k B T =" without calculating ∆ self-consistently. The figure of merit at low temperatures reaches very close to P 2 /(1 − P 2 ) unless P is very close to unity, but the exact temperature scale where this happens depends on the value of polarization. At the lowest temperatures ZT is limited by another spurious heat conduction process, due to nonzero density of states inside the gap, described here by the Dynes Γ parameter energy allows realizing a new type of a superconducting thermoelectric detector of radiation [33,283]. In addition, the thermoelectric effects could be used in non-invasive low-temperature thermometry [34] so that the measured thermovoltage would indicate the temperature profile. At room temperature a scanning thermometer based on thermoelectric effects was realized by [284]. Utilizing spin-polarized tunneling from spin-split superconductors would allow extending this technique to low temperatures.

Effect of spin mixing
Rather than the Dynes parameter, a more relevant limitation of the figure of merit in spin-split superconductors is due to the spin mixing caused by spinorbit and spin-flip scattering that are disregarded in the above results. This becomes an issue especially in heavy-metal based superconductors where the spin-orbit scattering is strong. As an example, Fig. 30 shows how the figure of merit decays as the spin-orbit relaxation becomes stronger. The high efficiencies are obtained only when the spin relaxation rates are small compared to the gap energy ∆. Note that for some materials the stronger spin-orbit relaxation due to the increase in the atomic number Z (say, compared to Al) is partially compensated by the increased critical temperature and therefore the energy gap.
As shown in in Fig. 31, the suppression of ZT is monotonous as a function of an increasing spin-orbit scattering. However, in the presence of a non-vanishing Dynes parameter, the spin-flip scattering may also lead to an increased ZT at low temperatures. For completeness, we also show the effect of spin mixing on the thermopower, i.e., the Seebeck coefficient. The results of the effect of spin-flip scattering are in line with those studied in [285].

Nonlinear response heat engine
The above calculations on the heat engine efficiency have been done in the strict linear response limit ∆T T . In this limit both the output power and the Carnot efficiency, proportional to ∆T /T , are vanishingly small. However, the device can have a rather large output power and efficiency also at nonlinear response. To illustrate this, we plot the maximum extractable output poweṙ W max = max V (−IV ) from a N-FI-S junction in Fig. 32 at large differences T S − T N between the spin-split superconductor and a normal-metal reservoir. The output power is somewhat larger for T S > T N than the opposite situation with the same temperature difference. The plotted quantity corresponds to the maximum power that is possible to extract from the junction after optimizing the load resistance -note, however, that this optimized resistance value depends on the exact temperature difference and other parameters. The corresponding efficiency η =Ẇ max /Q is plotted in Fig. 33. There we take into account only the heat current flowing through the junction itself. The efficiency is reduced if other heat relaxation mechanisms become relevant.

Effect of non-collinear magnetizations
Above we assume that the polarization P = Pû P of the tunneling contacts is collinear with the exchange field h = hû h in the spin-split superconductor. This is the limit that guarantees the largest thermoelectric current, but it is also interesting to consider how non-collinearity would show up in the thermoelectric response, in particular in the coefficientα of Eq. (147c). This coefficient is by definition a scalar quantity. Therefore, it can only depend on a scalar combination of the vectors P and h, or rather,û P andû h . This combination is the inner product,û P ·û h . Therefore, the generalization of Eq. (147c) to the non-collinear case isα nc =û P ·û hαc , whereα c is the thermoelectric coefficient in the collinear case, for example given by Eq. (147c).

Thermally induced spin currents
Besides the large thermoelectric effect discussed above, the contact between spin-split superconductors with other conducting materials can exhibit a large (longitudinal) spin Seebeck effect, where a temperature difference drives spin currents to/from the spin-split superconductor [243]. The spin current I S = I ↑ − I ↓ can be obtained from Eq. (142). In linear response, Eq. (59), it is described by the same coefficientα/P = α as the thermoelectric coefficient, Eq. (147c), but without the spin polarization P of the interface. When either of the two materials realizes an island, the spin current can be converted into a spin accumulation µ z that is determined from the balance between thermally induced spin currents and spin relaxation within the island. For weak spin relaxation, the resulting spin Seebeck coefficient µ z /∆T can be large [i.e., not proportional to the number of thermally excited quasiparticles ∼ e −(∆−h)/(k B T ) as the terms in Eqs. (147a)-(147c)], similar to the thermopower in Eq. (149). The spin currents induced in the case of two spin-split superconductors, and the additional effects of Josephson coupling, magnetization texture and spin-orbit effects are discussed in [286,287].
The above discussion of the thermoelectric cooling and heat engines disregards the spin accumulation caused by the nonequilibrium biases. This as-sumption is valid if the overall spin relaxation within the island, described by a spin-flip conductance G sf = e 2 ν F Ω/τ sn is much larger than the conductances of the tunnel junctions. Here ν F is the normal-state density of states in the Fermi level, Ω is the volume of the island, and τ −1 sn is the spin relaxation rate within the island. In the case of a spin-split superconducting island, both the conductance and the spin relaxation rate contain the same exponential term ∼ e −(∆−h)/(k B T ) , and for an estimate of the overall magnitude of the effect, it is hence enough to compare G sf in the normal state to the normal-state tunnel conductance G T . In other words, disregarding the spin accumulation is justified provided e 2 ν F Ω/(τ sn G T ) 1. On the other hand, most of the discussion disregards the effect of spin relaxation on the density of states of the spin-split superconductor. As shown in Fig. 31, spin relaxation starts to affect the results when /τ sn k B T c . These two constraints can be simultaneously satisfied depending on the precise value of τ sn , provided that G T (k B T c /δ I )e 2 /h, where δ I = 1/ν F Ω is the average energy level spacing on the island. This constraint is always satisfied in tunnel barriers between metallic systems.
To our knowledge, the spin Seebeck effect has not been directly measured in spin-split superconductors driven by a temperature difference. However, outside linear response, this is the effect that causes the long-range spin signal in spin injection experiments as discussed in Sec. 4.2. Besides affecting the nonequilibrium dynamics, the presence of spin accumulation in a superconductor with an exchange field couples to the self-consistency equation for the superconducting order parameter ∆ as in Sec. 3.7. Ref. [288] utilized this fact and predicted that the thermally induced spin accumulation in a spin-split superconductor can result to changes in the critical temperature: besides only suppressing it at low temperatures, it can enhance the critical temperature at some intermediate temperatures, or even lead to a situation where superconductivity shows up only between two critical temperatures.
This spin Seebeck effect should be contrasted to the analogous phenomenon discussed in non-superconducting materials [289]. There, a major contribution to the spin Seebeck signal comes from the thermally induced spin pumping from ferromagnetic insulators [290].

Summary
In this work we review the electronic transport properties of superconducting hybrid structures with a spin-split density of states. We mainly focus on ferromagnetic insulator-superconductor systems (FI-S) in which such a splitting field can be achieved without the need of an applied magnetic field. Best FI-S material combinations studied so far are the europium chalcogenides (EuO, EuS and EuSe) together with Aluminum films. Interesting of those materials is that thin films of EuO or EuS can also be used as highly efficient spin filters and therefore they are the best candidates for the realization of the heat engines proposed in Sec. 5.
One of the goals of this review is to provide a complete theoretical framework to study non-equilibrium effects in superconducting structures with spin-split density of states. In this respect, Sec. 3 gives a detailed introduction to the quasiclassical formalism for diffusive hybrid systems, and a description of the nonequilibrium modes. The advantages of using this formalism are multiple. On the one hand it provides a simple way to identify the quasiparticle nonequilibrium modes. As shown in Sec. 3.1, the combination between superconductivity and magnetism produces, besides the widely studied charge and energy modes, two additional ones: spin and spin energy. The quasiclassical equations show explicitly how all these modes couple to each other in different situations.
On the other hand, quasiclassical kinetic equations are a powerful tool for the study of hybrid multi-terminal systems, with different materials. Effects occurring at hybrid interfaces between different materials can be described by the help of effective boundary conditions. In Sec. 3.4 we discuss this issue and present boundary conditions suitable to describe a S-FI interface. The inclusion of the interfacial spin-orbit effect and the associated coupling between charge and spin would be another interesting further development of the formalism. The extension of the formalism to include additional types of interactions and couplings is in principle rather straightforward. For example, one can include magneto-electric effects associated with the charge-spin coupling in systems with linear-in-momentum spin-orbit coupling, by introducing an effective SU (2) gauge potential in the quasiclassical equations [291,292,112], or, in the case of an extrinsic spin-orbit coupling, by taking into account higher-order terms in the impurity potential [113]. Inclusion of such effects in a non-equilibrium situation and with time-dependent fields, is an interesting further development of the field [293].
We use the power of the quasiclassical approach in Sec. 4, where we study in detail different transport properties of FI-S structures and spin-split superconductors and contrast the results with experimental findings. Both charge and spin dependent properties are determined by the behaviour of the nonequilibrium modes of the supercoductor. Specifically we have computed the non-local conductance in a two terminal geometry and the spin Hanle effect in superconductors.
One further perspective of the present work is the extension of the Keldysh quasiclassical formalism used in this review to include magneto-electric effects associated with the spin-orbit coupling (SOC) [291,292,112,113,111,293], and the coupling magnetization dynamics to the electronic degrees of freedom via the reciprocal effects of spin transfer torque and spin pumping [227].
In Sec. 4.3.2 we discuss the dynamics of spin-split superconductors in rf fields, and review how nonequilibrium effects induced by them manifest in the quasiclassical theory. Historically, magnetic resonance effects in superconductors are well studied, but fewer experiments have probed spin-split thin films.
Another very interesting aspect in studying superconductors with a spinsplit density of states are their thermoelectric properties. In Sec. 5 we focus on thermoelectricity in superconducting structures, both in linear and non-linear regimes and discuss the realization of a heat engine based on FI/S structures. Because of their extremely large thermoelectric figure of merit, using S/FI structures has been proposed for accurate radiation sensing and non-invasive scan-ning thermometry. These recently proposed applications add up to the the rich physics offered by spin-split superconductors, which have been long used as tools to characterize equilibrium properties of magnets, especially their spin polarization. Most likely there are also many other applications to be uncovered, opened by the possibility for realizing a controlled combination of magnetism and superconductivity. The theoretical framework presented in this review should help in this task.
The major challenge in fabricating real devices is to find FI-S combinations with a large superconducting critical temperature and simultaneously a large spin splitting and efficient spin-filtering. Superconductors with high T c usually have larger spin-orbit coupling which leads to a broadening of the spin-split peaks in the density of states, as discussed in Sec. 2. Further material research is needed in this respect in order to find optimal material combinations.
Finally, the phenomena discussed in this review have a direct impact on several fields of condensed matter and quantum technologies, as for example the design of hybrid quantum materials and their use in functional devices. Transport properties of such systems can be controlled by the temperature or by external fields, and to interpret measurements it is essential to understand their nonequilibrium properties, which in turn depend on the quality of interface growth and material combination. This review provides most of the theoretical tools for such analysis.

S
superconductor FI ferromagnetic insulator G T tunneling conductance ν F normal-state density of states (per spin) at the Fermi level D diffusion constant σ N = 2e 2 ν F D normal-state conductivity in a diffusive wirě g 8×8 Green's function matrix in Keldysh-Spin-Nambû g 4×4 Green's function matrix in Spin-Nambu N σ (ε), σ =↑, ↓ normalized density of states per spin N + N ↑ + N ↓ = Re g 03 N − N ↑ − N ↓ = Re g 33 n F (ε) Fermi distribution function f eq = 1 − 2n F quasiclassical equilibrium distribution function f T (ε) non-equilibrium transversal distribution function f L (ε) longitudinal distribution f T j (ε), j = 1, 2, 3 spin transversal function σ j , j = 1, 2, 3 Pauli matrices in spin space τ j , j = 1, 2, 3 Pauli matrices in Nambu space g ij σ i τ j component of the retarded Green function τ elastic relaxation time τ sf spin-flip scattering time τ so relaxation time at impurities with spin-orbit coupling τ sn = [τ −1 sf + τ −1 so ] −1 normal state spin relaxation time R barrier resistance per unit area j c,e,s,se charge,energy,spin,spin-energy density currents I charge currenṫ Q energy current I s spin currenṫ Q s spin-energy current P polarization or spin-filter efficiency of FI interface ∆ superconducting order parameter ∆ 0 value of ∆ at T = h = 0 and without spin relaxation BCS critical temperature for h = 0 and without spin relaxation D Diffusion coefficient G electrical conductance G th thermal conductance ξ s superconducting coherence length