Out-of-Equilibrium Spin Transport in Mesoscopic Superconductors

# Out-of-Equilibrium Spin Transport in Mesoscopic Superconductors

C.H.L. Quay and M. Aprili
###### Abstract

The excitations in conventional superconductors, Bogoliubov quasiparticles, are spin-1/2 fermions but their charge is energy-dependent and, in fact, zero at the gap edge. Therefore, in superconductors (unlike normal metals) the spin and charge degrees of freedom may be separated. In this article, we review spin injection into conventional superconductors and focus on recent experiments on mesoscopic superconductors. We show how quasiparticle spin transport and out-of-equilibrium spin dependent superconductivity can be triggered using the Zeeman splitting of the quasiparticle density of states in thin-film superconductors with small spin-mixing scattering. Finally, we address the spin dynamics and the feedback of quasiparticle spin imbalances on the strength of the superconducting energy gap.

rsta.royalsocietypublishing.orgResearch

Article submitted to journal

Subject Areas:

low-temperature physics, materials science, mesoscopics, nanotechnology, quantum physics, spintronics

Keywords:

superconductivity, out-of-equilibrium, spintronics, mesoscopic physics

Author for correspondence:

Out-of-Equilibrium Spin Transport in Mesoscopic Superconductors

C.H.L. Quay and M. Aprili

Laboratoire de Physique des Solides, CNRS, Univ. Paris-Sud, University Paris-Saclay, 91405 Orsay Cedex, France

\@xsect

Out-of-equilibrium superconductivity has received a lot of attention since the 1970s [1, 2], with the Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity providing a framework for describing the excitations of a ‘conventional’ superconducting condensate of singlet Cooper pairs [3]. The advent of micro- then nanotechnology and their continuous improvement greatly facilitated the fabrication of devices and circuits with dimensions on the order of the excitations’ relaxation lengths of 1–100[4]. Out-of-equilibrium states have thus become achievable with small external perturbations, e.g. external voltage biases, and also more easily detected. In parallel, spin-dependent superconducting proximity effects in hybrid devices at equilibrium were discovered; this was enabled by progress in engineering superconductor/(ferro)magnet interfaces. Oscillations of the superconducting order parameter, ferromagnetic -junctions [5] as well as long-range odd-frequency triplet correlations [6] have been observed. These results belied received wisdom about the incompatibility between superconductivity and magnetism, and called for a better understanding of spin transport and polarization in conventional superconductors.

Work on out-of-equilibrium superconductivity has concentrated mainly on the energy and charge excitation modes associated, respectively, with heat and charge transport [7] and their applications in cryogenic detection of radiation and thermometry [8]. These excitations correspond to an excess of quasiparticles with respect to thermal ones (energy mode) [9] and an imbalance between the hole-like and electron-like excitations (charge mode) [10, 11]. Both energy and charge modes are spin degenerate, though quasiparticles are spin-1/2 fermions and obey Fermi-Dirac statistics at equilibrium — their distribution function is then given by , with the energy, the temperature and Boltzmann’s constant. (The Fermi energy is taken to be equal to zero.)

The energy mode can be excited by charge-neutral perturbations such as electromagnetic radiation whose frequency is larger than the superconducting energy gap ; the absorption of such radiation breaks pairs and creates quasiparticules [12]. The charge mode, on the other hand, can be excited by injecting charged carriers (i.e. electrons or holes) from a normal electrode into a superconductor, where they become quasiparticles [13]. As quasiparticles are not instantaneously converted into Cooper pairs, their chemical potential is shifted up or down with respect to that of Cooper pairs. This has been measured as a voltage drop between the superconductor and a normal electrode in contact with the superconductor via a tunnel barrier [14, 15]. If electrons or holes are injected at energy , both charge and energy modes are excited. The relaxation time for the energy mode is the inelastic scattering time [8] while charge imbalance relaxes over the pair-breaking time (or equivalently quasiparticle recombination time) [16].

In a visionary theoretical paper published in 1976, A. Aronov [17] introduced the concept of spin injection into superconductors. The main idea of his paper was to use spin injection to produce an internal magnetic field in order to perform an NMR experiment in the superconducting state. (The Meissner effect prevents external magnetic fields from penetrating into the volume of superconductors.)

In a pioneering experiment, Johnson et al. [18] showed in 1985 that spin injection from a ferromagnetic electrode into a normal metal is possible by applying a voltage bias across the interface between the two. The out-of-equilibrium magnetization created in the normal metal is detected electrically, by measuring the voltage between it and a second ferromagnetic electrode [19]. This ‘nonlocal’ signal is directly proportional to the shift in the chemical potential, , of spin up (down) electrons due to spin accumulation [20], [21]. (spin up and down chemical potentials shift by the same amount, but in opposite directions.) The spin relaxation length measured in high purity light metals (which have low spin-orbit coupling) can reach  [18] while the spin relaxation time is 50ns. Subsequently Johnson et al. extended spin injection into a superconductor close to its critical temperature  [22]. This not only raised the question of spin transport in superconductors, the spins being carried by the out-of-equilibrium quasiparticles, but also motivated a search for superconducting states in which singlet Cooper pairs coexist with spin-polarized quasiparticles [23].

Among excitation techniques besides those mentioned above, we note in particular charge and spin currents generated by the magnetic losses of the precessing magnetisation of a ferromagnet also called ‘spin pumping’ [24, 25, 26, 27, 28, 30, 29]. In the case of a ferromagnetic insulator in contact with a superconductor, a pure spin current with no charge should be injected [31]. In addition, as this ‘spin pumping’ technique injects spins at low energies, in contrast to spin-polarised current injection, it should in principle also result in little or no excitation of the charge mode. Nevertheless, while spin pumping into metals has been demonstrated experimentally, this technique has not been applied to superconductors.

\@xsect

Spin injection from a ferromagnetic electrode as in Johnson’s work [18] is always associated with charge transport and charge imbalance as it arises from a charge current of spin-polarised carriers. However, unlike normal metals, in which charge imbalance cannot exist, in superconductors both charge and spin imbalances are possible in the quasiparticle population. Fundamentally, this is because the superconducting condensate acts as a ‘charge reservoir’, such that quasiparticle charge need not be conserved [32]. (The total charge of both quasiparticles and Cooper pairs is conserved. Charge imbalance in the quasiparticle population thus means that the quasiparticle chemical potential is shifted with respect to that of the pairs [10, 13, 14, 32]). Note, in addition, that the charge of Bogouliobov spin-1/2 quasiparticles is energy dependent: ; therefore, while quasiparticles at the gap edge are neutral, high energy quasiparticles are (hole-like) or (electron-like) excitations [33]. This results in an effective spin-charge separation in the quasiparticle population as pointed out by Kivelson et al. [34]. As spin and charge can relax via different mechanisms and thus over different timescales, this leads to a peculiar situation in diffusive out-of-equilibrium superconductors in which charge and spin excitation modes are spatially separated.

A superconductor in good metallic (i.e. high transparency) contact with a ferromagnet experiences a strong suppression of superconductivity in a region up to a few superconducting coherence lengths from the interface, even at equilibrium [35, 36]. (The superconducting coherence length is in the clean limit and in the diffusive limit, with the diffusion constant, the Fermi velocity and Planck’s constant.) This phenomenon, also known as the inverse proximity effect, is due to the strong pair-breaking effect of the exchange field on singlet Cooper pairs [37] that leak from the superconductor into the ferromagnet, which feeds back through the boundary conditions on superconductivity at the interface [36].

Placing a tunnel barrier between the superconductor and the ferromagnet decouples the spin and Cooper pair reservoirs and suppresses undesirable proximity and inverse proximity effects [38], including the appearance of Andreev bound states in the superconducting gap [39]. This also affords a technical advantage: as a tunnel barrier results in a finite interface resistance, the interface can now be voltage biased and thus the quasiparticle injection energy chosen. Injection at energies lower than the superconducting energy gap is not possible due to the lack of quasiparticle states; the bias voltage must be thus larger than . Finally, tunneling through an insulator preserves spin orientation; this turns out to be essential for efficient spin injection [40]. In the following, ‘spin injection’ should be taken to mean spin injection through a tunnel barrier. As the spin and charge degrees of freedom are not directly coupled, we focus mainly on spin imbalance. Charge imbalance phenomena are treated in Ref.s [41].

Assuming negligible spin accumulation in the ferromagnetic electrode as well as a uniform temperature for the whole system and no charge imbalance (in the superconductor), the spin, ‘electrical charge’ and ‘quasiparticle charge’ currents across a ferromagnet-insulator-superconductor junction obtained from the Fermi Golden rule read [42, 43]:

 Ie(V)=2πeℏN0|T|2∫+∞−∞NBCS(E){NF[f(E−eV)−f(E+eV)]−δNF[f(E−μs)−f(E+μs)]}dE,Is(V)=πN0|T|2∫+∞−∞NBCS(E){δNF[f(E−eV)−f(E+eV)]−NF[f(E−μs)−f(E+μs)]}dE,IQP(V)=πℏN0|T|2∫+∞−∞q2(E)NBCS(E){NF[f(E−eV)−f(E+eV)]+δNF[f(E+μs)−f(E−μs)]}dE. (2.1)

Here is the voltage (applied or measured) across the junction; the electron charge, the density of states of the superconductor at the Fermi level when it is in the normal state; the BCS quasiparticle density of states; a chemical potential shift (of opposite sign for opposite spins ); the quasiparticle charge; and the junction transmission assumed to be identical for each spin. and account for the different density of states, and at the Fermi level for spin up and down in the ferromagnet . (The polarisation of the ferromagnet is encoded in ; thus, in the case where the ferromagnet is replaced by a normal metal electrode .) The assumption of a uniform system temperature implies fast thermalization [44].

The difference between the electrical charge current and the quasiparticle charge current is carried by the superconducting condensate [45]. These equations apply equally to both injector (‘local’) and detector (‘nonlocal’) junctions: At the injector the currents are determined by the applied voltage , while at the latter a voltage is measured in presence of spin imbalance. can be found by inverting the Equations 2.1 and imposing the condition , i.e. zero charge current through the detector as is usual for a voltage probe.

Spin accumulation appears in the superconductor — in a dynamic equilibrium between injection and relaxation — and can be described by: , with the spin imbalance and the spin imbalance relaxation time. In principle, this equation is coupled to Equations 2.1 as is directly linked to , and the whole system of equations must be solved self-consistently. However if the superconductor is only slightly out-of-equilibrium (i.e. , ) can be neglected in Equations 2.1 and in the stationary limit (i.e. ) we obtain , being the polarisation of the ferromagnetic leads (assumed to be the same for the injector and the detector). Unlike in normal metals, spin accumulation in superconductors is a non-linear function of the injection voltage , reflecting the voltage dependence of the junction conductance set by the BCS quasiparticle density of states.

The spin imbalance contribution to the nonlocal voltage can be isolated by measuring the latter when the magnetisations of injector and detector are aligned () and anti-aligned (). Inverting Equations 2.1 one gets, in the limit , :

 VAPNL(VL)−VPNL(VL)=P2ρsgNSΩIe(VL) (2.2)

where is the local charge current, the spin resistivity, the normalised conductance of the detector junction and the injection volume [41]. We can define a spin resistance . The result of Eq. 2.2 for spin accumulation arising from a spin injection current agrees with earlier calculations by Maekawa et al. [43], and by Zhao and Hershfeld. [42]. This signal is directly proportional to . Thus, just as the electrical resistance probes the momentum relaxation time, so also the spin resistance probes .

The question arises as to the microscopic mechanisms which determine the spin imbalance relaxation time in a superconductor. In metals, spin relaxation is generally due to spin-flip scattering off magnetic impurities and/or spin-orbit scattering [46, 47, 48]. (Note that the spin-orbit interaction is time-reversal invariant and must therefore be combined with spatial or momentum-space inhomogeneities in order to flip spins.) In superconductors, interactions between quasiparticles and Cooper pairs also may affect spin relaxation. Calculations in the limit in which the out-of-equilibrium distribution function is practically unchanged from the equilibrium one (i.e. much smaller than all other energy scales) have shown a monotonic increase of the spin-flip time for spin-orbit scattering when temperature decreases [49, 50]. shows instead a minimum below when spin-flip by magnetic impurities is the dominant relaxation mechanism [42]. (These calculations are similar to those for in nuclear magnetic resonance [51].) Calculations within the quasi-classical approximation pointed to an enhanced spin-flip scattering at the gap edge leading to a strong suppression of the spin imbalance at low temperature [52].

Early experiments on superconducting Al [53] and Nb [54] junctions of macroscopic dimensions (area ) showed that spin accumulation produces an excess of quasiparticles. This was interpreted as a longer recombination time due to a very long quasiparticle spin lifetime (i.e. , the inelastic scattering time). In other words, spin relaxation was proposed as the bottleneck for quasiparticle relaxation.

\@xsect

The long spin relaxation time measured in light metals at low temperatures (compared to that in heavy metals), and also indications of an increase of the spin lifetime in the superconducting state, suggest that the superconducting spin relaxation length could be on the order of m, a length scale accessible to standard electron lithography. The spatial decay of spin imbalances can thus be probed directly with detector electrodes placed at different distances from the injection point. Spin-polarised electrons diffuse from the injection point (or area) and relax over a length . When both injector and detector are ferromagnets, the device is known as a ‘lateral spin valve’. In superconductors, as mentioned above, charge and spin imbalances can diffuse (from the injection point) and relax over different length scales.

Spin imbalance in lateral spin valves depends on the distance (of the detector) from the injection point. The spin resistance decays exponentially with and is given by:

 Rs(x)=P2ρsgNSAxse−x/xs (3.1)

where is the cross-sectional area of the metallic or superconducting wire perpendicular to the main direction of diffusion (The system is assumed to be quasi-one-dimensional.)

Generalising the Johnson technique to lateral spins valves Jedema et al. [55, 56] showed that the spin imbalance relaxes exponentially as expected. They obtained characteristic decay lengths of 650 nm at 4.2 K and 350 nm at room temperature. This corresponds to spin imbalance relaxation times of 100 ps and 45 ps respectively in good agreement with previous experiments. A measurement of the spin resistance in a lateral spin valve and its -dependence is shown in Figure 1, with the Al in the normal state.

Following Johnson [57], Jedema et al. also used the Hanle effet (coherent spin precession around a magnetic field perpendicular to the spin polarisation) [56] to directly measure the spin imbalance relaxation time.

 Rs=P2ρsgNSAxs∫+∞0P(t)cos(ωLt)e−x/xsdt (3.2)

where , with the Bohr magneton and the applied magnetic field. is the probability that, once injected at the injector F1 (cf. Figure 1), an electron arrives at the detector F2 after a diffusion time . For a diffusive wire .

Unfortunately the Hanle effect cannot be measured in superconducting thin films as superfluid screening (the Meissner effect) prevents perpendicular magnetic fields from uniformly penetrating the film. In addition, (nonuniform) perpendicular magnetic fields strongly modify the density of states in the film complicating data interpretation. (Field penetration would be uniform in superconducting nanowires whose diameters are much smaller than the penetration depth, regardless of field orientation; however, we are not not aware of spin injection experiments in these systems.)

Measurements on superconducting lateral spin valves also suffer from stray fields due to the magnetisation of the ferromagnetic electrodes; this can strongly affect superconductivity in a spatially nonuniform manner. As spin resistance is defined as the difference in the nonlocal signal between the parallel and anti-parallel detector-injector magnetisation configurations, it is of the essence to be able to distinguish between out-of-equilibrium spin-dependent effects and equilibrium magnetic landscape effects.

Designed to minimize stray field effects while switching between anti-parallel and parallel configurations, the experiment performed by Poli et al. [58] yielded the nonlocal spin resistance in the superconducting state of a lateral spin-valve similar to that used by Jedema et al. [56] and Valenzuela et al. [59]. The spin relaxation length was observed to decrease in the superconducting state with respect to the normal state, consistent with spin relaxation dominated by spin-flip scattering from magnetic impurities (cf. Ref. [52]). In contrast, a recent experiment by Wakamura et al. [60] showed that the spin relaxation time increases by a factor of at least 4 with respect to the normal state in Nb, a superconductor with strong spin-orbit coupling. In this experiment, a (superconducting) wire between the injector and detector of a lateral spin valve serves as a spin current sink; this ‘spin absorption’ technique [61, 62] allows an estimate of the spin imbalance relaxation time in materials for which is not lithographically accessible.)

\@xsect

An external magnetic field can be used not only to control the magnetisation orientation of ferromagnetic electrodes but also to change the quasiparticle excitation spectrum in the superconductor. In thin films whose thickness is smaller than  [63], orbital screening effects due to an in-plane magnetic field can be neglected. This condition on , known as the Pauli limit, comes from imposing that the diamagnetic energy be smaller that the Zeeman energy. It results in smaller than the condition the magnetic field penetration depth, which ensures a homogeneous magnetic field in the superconductor. The (dominant) Zeeman effect then shifts quasiparticle energies up or down depending on their spin orientation with respect to the applied field , , with the Zeeman energy [64]. The resulting spin-splitting of the BCS density of states was first observed by Tedrow and Meservey in tunnelling spectroscopy measurements [65].

Building on the concept introduced by Huertas-Hernando et al. [66] of the ‘absolute spin valve’, Giazotto et al. [67] pointed out that the Zeeman spitting on the BCS density of states can turn a superconductor into an excellent spin source. Current polarisations of up to 100% have been predicted in the tunneling limit (low transparency) for a bias voltage .

Two experiments have followed this idea [68, 69]. Introducing a Zeeman-split BCS DOS in Equations 2.1; assuming , , and ; and expanding in , the spin resistance is:

 Rs(VL)=P2gnsΩ[Pρs2+dNBCS(VL)/dVLNBCS(VL)μBHρs1+...] (4.1)

This expression is particularly suggestive: the spin imbalance can be seen here to depend clearly on the polarisation of the injector electrode (first term) as well as on the Zeeman splitting of the DOS in the superconductor (second term). Here we have introduced two different spin resistivities and associated with two different relaxation times and in the first and second terms, respectively. Physically this corresponds to different relaxation mechanisms for spin imbalances in the presence or absence of a Zeeman-induced splitting of the BCS DOS. The shape of the signal measured in the experiment (Figure 2.) and its amplitude as a function of the applied magnetic field in the limit of small field agrees well with Eq. 4.1 provided the off-set produced by the first term is negligible. obtained from the data is 10ns much longer than the spin relaxation time measured in the normal state. Furthermore  [70] consistent with spin-charge separation as suggested theory [42].

Note that, as spin polarisation is provided mainly by the Zeeman splitting of the superconducting density of states, the injector electrode no longer needs to be a ferromagnet. This can be seen in Figure 3, where we present data from a device, similar to that shown in Figure 1, but in which the injector is normal. The thickness of the superconducting Al wire is reduced to 8.5 nm to increase the critical field from 650 mT to 2100 mT. While the spin accumulation first increases with increasing field (Figure 3c), the charge imbalance monotonically vanishes with field as expected for increasing orbital depairing (Figure 3 d). Note also that the signals from spin and charge imbalance have different symmetries, even and odd respectively, in agreement with Equation 2.1.

Spin currents carried by quasiparticles can be converted into transverse charge currents by the inverse spin Hall effect [71], resulting in charge accumulation at the transverse edges of the superconductor [72]. The spin Hall effect originates from side jump and skew scattering induced by spin-orbit coupling [74, 62, 73]. Quasiparticle-mediated spin Hall and inverse spin Hall effects are expected to be greatly enhanced in the superconducting state compared to the normal state [72], as the nonlocal spin resistance is directly proportional to the quasiparticle density that decreases exponentially with temperature below , as . Very recently an increase of three orders of magnitude in the nonlocal inverse spin Hall resistance of a NbN wire well below (compared to the normal state) was reported [75]. The inverse Hall signal was observed only when injector and detector electrodes were closer together than the charge relaxation length.

\@xsect

The experiments by Quay et al.[69] and Hübler et al. [68] show unequivocally that a magnetization appears in a mesoscopic superconductor in the Pauli limit when the superconductor is driven out-of-equilibrium by a tunneling current from both a ferromagnetic and a normal electrode (see also Ref.s [76] and [77]).

These experiments are consistent with — and were interpreted in the framework of — the spin imbalance model presented above, in which spin imbalance appears as a chemical potential shift which is opposite for opposite spins, thus bracketing out the question of the spin-dependent energy distribution function.

In a simple picture, non-equilibrium excitations in Zeeman-split superconductors can be described by four parameters together quantifying the charge and spin degrees of freedom of the system: , , and . This presupposes that the electron-electron interaction time the electron-phonon time, and that the distribution functions can thus be approximated by Fermi-Dirac distributions, and generalises the even and odd modes of the non-Zeeman-split case by taking spin into account. The conceptually simplest spin imbalances are then due to (or associated with) either and as in Ref.s [68, 69, 45], or and both different from the equilibrium (bath) temperature (Figure 4). In both cases, a finite voltage will be measured between the superconductor and a ferromagnetic detector with which it is in contact via a tunnel barrier. The latter scenario, which corresponds to a spin-dependent thermoelectric effect, is a result of the broken electron-hole symmetry of the DOS of each spin in the superconductor due to the Zeeman field.

A complete description of the non-equilibrium excitations requires the knowledge of the spinful distribution functions and , where the need not be Fermi-Dirac not share the same chemical potential.

Indeed, by solving Boltzmann diffusion equations, and taking into account local modification of superconductivity due to the out-of-equilibrium distribution functions, Krishtop et al. [78] found that nonlocal signals with similar bias dependence to those measured in Ref.s [68, 69] are also expected in the second case above [78], as well as the more general case where is not a Fermi-Dirac distribution. Bobkova et al. [79] and Silaev et al. [80] find similar behaviour for the nonlocal signal using the quasiclassical Keldysh-Usadel formalism. (See also Ref.s [81, 82].)

Silaev et al. [80] further proposed a description of the out-of-equilibrium state in the following basis: charge, spin, energy and spin energy, with spin and energy modes contributing to the observed signal. The former relaxes via elastic spin-flip processes; the latter relaxes via inelastic processes (electron-electron or electron-phonon scattering) and is unaffected by elastic scattering.

All of this recent theoretical work point to the electron-phonon relaxation time as an upper bound for the spin imbalance relaxation time. This is a natural explanation for the long times measured by Quay et al.[69] and Hübler et al. [68]: however, the precise ‘modes’ excited and out-of-equilibrium distribution functions produced remain to be experimentally investigated.

The spin imbalance relaxation length is observed to increase with the applied magnetic field, particularly when the orbital depairing is strong, and is explained by Bobkova et al. [79] as follows: Quasiparticle scattering occurs primarily at the gap edge with rates which go as for electron-electron scattering and for electron-phonon scattering. A reduced gap results in a diminished inelastic scattering rate and thus an increased lifetime. The upper limit for the spin relaxation time becomes the normal state energy relaxation time at that temperature. However this characteristic time scale corresponds to the lifetime of the out-of-equilibrium magnetization in the superconductor rather than the spin-flip time of each quasiparticle. This is similar to in nuclear magnetic resonance terminology [51].

\@xsect

Spin resonance experiments allow a direct measurement of the spin coherence time,  [83, 84]. In a typical electron spin resonance (ESR) experiment, electrons are immersed in an external homogenous static magnetic field, . Microwave radiation creates a perturbative transverse magnetic field (perpendicular to the static field) of frequency . The power absorbed by the spins from the microwave field is determined, usually by measuring the fraction of the incident microwaves that is not absorbed, i.e. either transmitted or reflected [85]. When is tuned to its resonance value, —with the gyromagnetic ratio —the electron spins precess around and is maximal. if proportional to the full-width at half-maximum of the resonance [83].

In a very recent experiment Quay et al. [86] measured the spin resonance of out-of-equilibrium quasiparticles using two ‘on-chip’, local microwave powermeters: a supercurrent measurement and the conductance across a normal metal-insulator-superconductor junction. At resonance, microwave radiation is absorbed by the precessing out-of-equilibrium quasiparticle magnetization.

Some data from Ref. [86] are shown in Figure 5. Figures 5a and Figure 5b show spin resonances obtained for two devices with superconducting Al bars of different thicknesses ( and ). The resonant field as a function of the microwave frequency from these and other measurements is reported in Figure 5c. A linear fit to the data gives a g-factor of 1.950.2, consistent with previous measurements of electrons in Al in the normal state (Figure 5c). The spin relaxation time is also shown in Figure 5c. = 9520ps (7015ps) for Al which was 8.5nm (6nm) thick. Both the order of magnitude of , as well as the fact that it is inversely proportional to the film thickness, are consistent with spin coherence limited primarily by the Elliott-Yafet spin-orbit scattering time  [46, 47]. Further, these figures agree with an independent, equilibrium measurement of (Figure 5d).

We note that, unlike normal metals where the spin imbalance relaxation time, in superconductors is much larger than , and that is of the order of the spin imbalance relaxation time measured in the normal state [56, 69], consistent with plural sources and relaxation mechanisms of spin imbalance in the superconducting state [87].

\@xsect

The influence of out-of-equilibrium quasiparticle distributions, including spin imbalance, is in principle accounted for by the self-consistent BCS gap equation [32]:

 (7.1)

where is the Debye energy and the BCS interaction strength.

The solution in the presence of a Zeeman/Pauli (as opposed to orbital) magnetic field was given by Sarma in 1963 [23]. In the case where a spin imbalance can be described by , while the quasiparticle temperature is the bath temperature, acts as an effective magnetic field and simply adds to the applied field  [88, 89, 90]. The effects of the two thus add or subtract depending on their relative sign. The case of cancellation has been described as a ‘recovery’ of superconductivity through spin injection [89]. We note, however, that in the case described here, where spin imbalance is provided by Zeeman-splitting due to the same magnetic field, the two effects have the same sign and there is no recovery.

More generally, as already pointed out above, the non-equilibrium spin-dependent distribution functions due to spin injection have to be calculated by solving kinetic equations, which take spin diffusion into account. This is particularly important for mesoscopic devices, as quasiparticles thermalise internally, which is to say recover a Fermi-Dirac distribution, over distances larger than the electron-electron scattering length (typically of the the order m at very low temperature) [8].

In order to observe the predicted modification of the gap due to spin imbalance, a number of experiments have been performed on high-Tc superconductors. Indeed, it has been predicted that spin injection results in a greater suppression of d-wave superconductivity than of conventional s-wave superconductivity, mainly due to injection of low-energy excitations at the order-parameter nodes [91]. Most experiments in this direction have focused on reductions of the critical current with spin injection [92]. V.A. Vas’ko et al. [93] showed complete suppression of superconductivity in LaSrMnO/LaCuO/DyBaCuO heterostructures. Similar results have been obtained in Au/YBaCuO/LaAlO/NdSrMnO heterostructures [94, 95]. There is nevertheless some debate around these results, and in particular whether they are due to or simply to spin-independent out-of-equilibrium quasiparticles [96]. To our knowledge, there has not yet been a spin-resolved measurement of the distribution function correlated with spectroscopic evidence of gap reduction by spin injection.

Coming back to the case of fast thermalization (i.e. quasiparticles thermalise quickly amongst themselves and with the environment compared to the elastic spin-flip time, and the spin imbalance is given by ), at zero temperature the gap is independent of for while at finite temperature the system goes through a first order transition to the normal state at  [88, 97]. In addition, Sarma pointed out that for Zeeman-split superconductors at zero temperature the self-consistent gap equation admits a solution corresponding to conventional Cooper pairs occupying the energy region while a shell around the Fermi surface is occupied only by spin down electrons [23]. This solution appears as a mixture of Pauli paramagnetism and a BCS superconductor with a field-dependent gap.

The Sarma state and the related ‘breached pair’ [98] state have received a great deal of attention in the context of BCS superfuidity in cold atoms where the imbalance between spin up and down densities can be finely controlled [99, 100, 101]. Evidence of a core-shell configuration where the BCS superfluid occupies the core while the shell is a gas of spin-polarised atoms has been shown by microscopy [102]. There has, however, yet to be theoretical or experimental work on inhomogeneous superconducting states in which spins and Cooper pairs are separated in real or momentum space. We note that, in addition to cold atoms, BCS-like theories with spin/population imbalance have also been applied to dense quark matter [103, 104]; they are thus of quite general interest.

\@xsect

Out-of-equilibrium superconductivity has been an active area of research since the 1970s; however, comparatively little attention has been paid to spinful excitations. The past few years have seen a revival of interest in the latter, in the wake of experiments demonstrating out-of-equilibrium magnetisation, spin-charge separation and quasiparticle resonance in superconducting aluminium in the Pauli limit. With close analogues in cold atom systems and dense quark matter, spin imbalanced superconductors are of interest to the broader physics community beyond condensed matter.

Nevertheless, comprehensive theoretical study of out-of-equilibrium superconductivity where spin degeneracy is lifted by Zeeman splitting, which includes the self-consistent calculation of the spatial dependence of the superconducting gap, is still a work in progress. Apart from traditional quasiclassical approaches, attempts to apply hydrodynamical methods to the problem are also underway [105]. Among open questions are the mutual interactions between the condensate and spinful excitations and the existence of novel ground states, including those due to the influence of spin imbalances on the internal spin structure of the superconducting order parameter [106, 107].

Authors’ Contributions. Both authors contributed equally to this manuscript; all data shown have been published elsewhere (see citations).

Competing Interests. The authors declare that they have no competing interests.

Funding. The authors’ and collaborators’ work mentioned in this review was funded by European Research Council Starting Independent Researcher (NANO-GRAPHENE 256965) and Synergy Grants, a C’NANO grant (DYNAH) from the Ile-de-France region, ANR Blanc grants (DYCOSMA and MASH) from the French Agence Nationale de Recherche, and the Netherlands Organization for Scientific Research (NWO/OCW).

Acknowledgements. The authors’ work mentioned in this review benefitted from collaboration with Cristina Bena, Denis Chevallier, Yann Chiffaudel, Clément Dutreix, Marine Guigou, Christoph Strunk, Mircea Trif and Maximilian Weideneder. We also acknowledge helpful discussions with Detlef Beckmann, Jaroslav Fabian, Julien Gabelli, Manuel Houzet, Hervé Hurdequint, Tatiana Krishtop, Brigitte Leridon, Jérôme Lesueur, Julia S. Meyer, Yuli Nazarov, Bertrand Reulet and Stanislas Rohart.

## References

• 1 Chang JJ & Scalapino DJ 1978 Nonequilibrium superconductivity. Journal of Low Temperature Physics 31, 1–32.
• 2 Gray KE, ed. 1981 Nonequilibrium Superconductivity, Phonons, and Kapitza Boundaries. Springer.
• 3 Kopnin N. 2009 Theory of Nonequilibrium Superconductivity, volume 110 of International Series of Monographs on Physics. Oxford University Press.
• 4 Kaplan SB et al. 1976 Quasiparticle and phonon lifetimesin superconductors. Physical Review B 14, 4854-4873.
• 5 Aprili M, Kontos T, Della Rocca ML, Lesueur J, Guichard W, Gandit P, Bauer A, Strunk C. 2006 Superconductivity and macroscopic quantum effects in superconducting/ferromagnetic hybrid nanostructures. Compte-rendus de l’Academie des Sciences 7, 116-125.
• 6 Linder J, Robinson JWA. 2015 Superconducting spintronics. Nature Physics 11, 307–315.
• 7 Schmid A, Schon G. 1975 Linearized kinetic equations and relaxation processes of a superconductor near Tc. Journal of Low Temperature Physics 20, 207–227.
• 8 Giazotto F, Heikkila TT, Luukanen A, Savin AM, Pekola JP. 2006 Opportunities for mesoscopics in thermometry and refrigeration: Physics and applications. Reviews of Modern Physics 78, 217–274.
• 9 Parker WH. 1975 Modified heating theory of nonequilibrium superconductors. Physical Review B 12, 3667-3672.
• 10 Tinkham M, Clarke J. 1972, Theory of Pair-Quasiparticle Potential Difference in Nonequilibrium Superconductors, Physical Review Letters 28, 1366-1369.
• 11 Owen CS, Scalapino DJ. 1972 Superconducting states under the influence of External Dynamic Pair Breaking. Physical Review Letters 28, 1559–1561.
• 12 de Visser PJ, Baselmans JJA, Diener P, Yates SJC, Endo A, Klapwijk TM 2011 Number Fluctuations of Sparse Quasiparticles in a Superconductor Phys. Rev. Lett. 106, 167004.
• 13 Tinkham M. 1972, Tunneling Generation, Relaxation and Tunneling Detection of Hole-Electron Imbalance in Superconductors, Physical Review B 6, 1747–1756.
• 14 Clarke, J. 1972, Experimental Observation of Pair-Quasiparticle Potential Difference in Nonequilbrium Superconductors, Physical Review Letters 28, 1363–1366.
• 15 Hübler F, Lemyre JC, Beckmann D, v Löhneysen H. 2010 Charge imbalance in superconductors in the low-temperature limit. Physical Review B 81, 184524.
• 16 Lemberger TR. 1984 One-to-one correspondence of charge-imbalance relaxing mechanisms with pair-breaking mechanisms in superconductors. Physical Review B 29, 4946–4950.
• 17 Aronov AG. 1976 Spin injection and polarization of excitations and nuclei in superconductors. Journal of Experimental and Theoretical Physics 44, 193.
• 18 Johnson M, Silsbee RH. 1985 Interfacial charge-spin coupling: Injection and detection of spin magnetization in metals. Physical Review Letters 55, 1790–1793.
• 19 Johnson M, Silsbee RH. 1988 Spin-injection experiment. Physical Review B 37, 5326–5335.
• 20 van Son PC, van Kempen H, Wyder P. 1987 Boundary Resistance of the Ferromagnetic-Nonferromagnetic Metal Interface. Physical Review Letters 58, 2271–2273.
• 21 Valet T, Fert A. 1993 Theory of the perpendicular magnetoresistance in magnetic multilayers. Physical Review B 48, 7099–7113.
• 22 Johnson M. 1994 Spin coupled resistance observed in ferromagnet/superconductor/ferromagnet trilayers. Applied Physics Letters 65, 1460–1462.
• 23 Sarma G. 1963 On the influence of a uniform exchange field acting on the spins of the conduction electrons in a superconductor. Journal of Physics and Chemistry of Solids 24, 1029–1032.
• 24 Tserkovnyak Y, Brataas A, Bauer GEW. 2002 Enhanced Gilbert Damping in Thin Ferromagnetic Films. Physical Review Letters 88, 117601.
• 25 Trif M, Tserkovnyak Y. 2013 Dynamic Magnetoelectric Effect in Ferromagnet/Superconductor Tunnel Junctions. Physical Review Letters 111, 087602.
• 26 Richard C, Houzet M, Meyer JS. 2012 Andreev Current Induced by Ferromagnetic Resonance. Physical Review Letters 22, 057002.
• 27 Moriyama T, Cao R, Fan X, Xuan G, Nikolic BK, Tserkovnyak Y, Kolodzey J, Xiao JQ. 2008 Tunnel Barrier Enhanced Voltage Signal Generated by Magnetization Precession of a Single Ferromagnetic Layer. Physical Review Letters 100, 067602.
• 28 Xiao J, Bauer GEW, Brataas A. 2008 Charge pumping in magnetic tunnel junctions: Scattering theory. Physical Review B 77, 180407.
• 29 Tserkovnyak Y, Moriyama T, Xiao JQ. 2008 Tunnel-barrier-enhanced dc voltage signals induced by magnetization dynamics in magnetic tunnel junctions. Physical Review B 78, 020401.
• 30 Tserkovnyak Y, Brataas A, Bauer GEW, Halperin BI. 2005 Nonlocal magnetization dynamics in ferromagnetic heterostructures. Reviews of Modern Physics 77, 1375–1421.
• 31 Kajiwara Y, Harii K, Takahashi S, Ohe J, Uchida K, Mizuguchi M, Umezawa H, Kawai H, Ando K, Takanashi K, Maekawa S, Saitoh E. 2010 Transmission of electrical signals by spin-wave interconversion in a magnetic insulator. Nature 464, 262–266.
• 32 Tinkham M. 1996 Introduction to Superconductivity. Mineola: Dover, 2 edition.
• 33 Ronen Y, Cohen Y, Kang JH, Haim A, Rieder MT, Heiblum M, Mahalu D, Shtrikman H. 2015 What is the charge of a quasiparticle in a superconductor? arXiv:1505.03147 [cond-mat]
• 34 Kivelson SA, Rokhsar DS. 1990 Bogoliubov quasiparticles, spinons, and spin-charge decoupling in superconductors. Physical Review B 41, 11693–11696.
• 35 Strunk C, Suergers C., Paschen U. and H. v. Loehneysen 1994 Superconductivity in Layered Nb/Gd films Physical Review B 49, 4053.
• 36 Buzdin A. 2005 Proximity effects in superconductor-ferromagnet heterostructures Review of Modern Physics 77, 935.
• 37 Kontos T, Aprili M, Lesueur J, Grison X, Dumoulin L. 2004 Superconducting Proximity Effect at the Paramagnetic-Ferromagnetic Transition. Physical Review Letters 93, 137001.
• 38 Latempa R., Aprili M. and Petkovic I. 2009 Quasiparticle trapping in three terminal ferromagnetic tunneling devices J. Appl. Phys. 106, 103925.
• 39 Blonder GE, Tinkham M, Klapwijk TM. 1982 Transition from metallic to tunneling regimes in superconducting microconstrictions: Excess current, charge imbalance, and supercurrent conversion. Physical Review B 25, 4515–4532.
• 40 Tedrow PM, Meservey R. 1971 Spin-Dependent Tunneling into Ferromagnetic Nickel. Physical Review Letters 26, 192–195.
• 41 Clarke, J. 1981 Charge Imbalance in Nonequilibrium Superconductivity, Phonons, and Kapitza Boundaries, Gray KE, ed. Springer.
• 42 Zhao HL, Hershfield S. 1995 Tunneling, relaxation of spin-polarized quasiparticles, and spin-charge separation in superconductors. Physical Review B 52, 3632–3638.
• 43 Maekawa S, Takahashi S, Imamura H. 2002 Spin injection into superconductors. Journal of Physics D: Applied Physics 35, 2452.
• 44 Tserkovnyak Y, Brataas A. 2002 Current and spin torque in double tunnel barrier ferromagnet-superconductor-ferromagnet systems. Physical Review B 65, 094517.
• 45 Chevallier D, Trif M, Dutreix C, Guigou M, Quay CHL, Aprili M, Bena C. 2014 Superconductor spintronics: Modeling spin and charge accumulation in out-of-equilibrium NS junctions subjected to Zeeman magnetic fields.
• 46 Elliott RJ. 1954 Theory of the Effect of Spin-Orbit Coupling on Magnetic Resonance in Some Semiconductors. Physical Review 96, 266–279.
• 47 Yafet Y. 1963 g Factors and Spin-Lattice Relaxation of Conduction Electrons. In Solid State Physics (ed.s F Seitz and D Turnbull), volume 14, pp. 1–98. Academic Press.
• 48 Fabian J, Das Sarma S. 1999 Phonon-Induced Spin Relaxation of Conduction Electrons in Aluminum. Physical Review Letters 83, 1211–1214.
• 49 Yafet Y, Vier DC, Schultz S. 1984 Conduction electron spin resonance and relaxation in the superconducting state. Journal of Applied Physics 55, 2022–2024.
• 50 Maki K. 1973 Theory of Electron-Spin Resonance in Gapless Superconductors. Physical Review B 8, 191–199.
• 51 Hebel LC, Slichter CP. 1959 Nuclear Spin Relaxation in Normal and Superconducting Aluminum. Physical Review 113, 1504–1519.
• 52 Morten JP, Brataas A, Belzig W. 2004 Spin transport in diffusive superconductors. Physical Review B 70, 212508.
• 53 Yang H, Yang SH, Takahashi S, Maekawa S, Parkin SSP. 2010 Extremely long quasiparticle spin lifetimes in superconducting aluminium using MgO tunnel spin injectors. Nature Materials 9, 586–593.
• 54 Leridon B, Lesueur J, Aprili M. 2005 Spin-bottleneck due to spin-charge separation in a superconductor. Physical Review B 72, 180505.
• 55 Jedema FJ, Filip AT, van Wees BJ. 2001 Electrical spin injection and accumulation at room temperature in an all-metal mesoscopic spin valve. Nature 410, 345–348.
• 56 Jedema FJ, Heersche HB, Filip AT, Baselmans JJA, van Wees BJ. 2002 Electrical detection of spin precession in a metallic mesoscopic spin valve. Nature 416, 713–716.
• 57 Johnson M, Silsbee RH. 1988 Coupling of electronic charge and spin at a ferromagnetic-paramagnetic metal interface. Physical Review B 37, 5312.
• 58 Poli N, Morten JP, Urech M, Brataas A, Haviland DB, Korenivski V. 2008 Spin Injection and Relaxation in a Mesoscopic Superconductor. Physical Review Letters 100, 136601.
• 59 Valenzuela SO, Tinkham M. 2004 Spin-polarized tunneling in room-temperature mesoscopic spin valves. Applied Physics Letters 85, 5914–5916.
• 60 Wakamura T, Hasegawa N, Ohnishi K, Niimi Y, Otani Y. 2014 Spin Injection into a Superconductor with Strong Spin-Orbit Coupling. Physical Review Letters 112, 036602.
• 61 Morota M, Niimi Y, Ohnishi K, Wei DH, Tanaka T, Kontani H, Kimura T, Otani Y. 2011 Indication of intrinsic spin Hall effect in $4d$ and $5d$ transition metals. Physical Review B 83, 174405.
• 62 Hoffmann A. 2013 Spin Hall Effects in Metals. Magnetics, IEEE Transactions on 49, 5172–5193.
• 63 Fulde P. 1973 High field superconductivity in thin films. Advances in Physics 22, 667–719.
• 64 Maki K, Tsuneto T. 1964 Pauli Paramagnetism and Superconducting State. Progress of Theoretical Physics 31, 945–956.
• 65 Meservey R, Tedrow PM, Fulde P. 1970 Magnetic Field Splitting of the Quasiparticle States in Superconducting Aluminum Films. Physical Review Letters 25, 1270–1272.
• 66 Huertas-Hernando D, Nazarov YV, Belzig W. 2002 Absolute Spin-Valve Effect with Superconducting Proximity Structures. Physical Review Letters 88, 047003.
• 67 Giazotto F, Taddei F. 2008 Superconductors as spin sources for spintronics. Physical Review B 77, 132501.
• 68 Hübler F, Wolf MJ, Beckmann D, v Löhneysen H. 2012 Long-range spin-polarized quasiparticle transport in mesoscopic al superconductors with a Zeeman splitting. Physical Review Letters 109, 207001.
• 69 Quay CHL, Chevallier D, Bena C, Aprili M. 2013 Spin imbalance and spin-charge separation in a mesoscopic superconductor. Nature Physics 9, 84–88.
• 70 Kleine A, Baumgartner A, Trbovic J, Golubev DS, Zaikin AD, Schönenberger C. 2010 Magnetic field and contact resistance dependence of non-local charge imbalance. Nanotechnology 21, 274002.
• 71 Dyakonov MI, Perel VI. 1971 Current-induced spin orientation of electrons in semiconductors. Physics Letters A 35, 459–460.
• 72 Takahashi S, Maekawa S. 2011 Spin Hall Effect in Superconductors. Japanese Journal of Applied Physics 51, 010110.
• 73 Sinova J, Valenzuela SO, Wunderlich J, Back C, Jungwirth T. 2015 Spin Hall effects. Reviews of Modern Physics 87, 1213–1260.
• 74 Žutić I, Fabian J, Das Sarma S. 2004 Spintronics: Fundamentals and applications. Reviews of Modern Physics 76, 323–410.
• 75 Wakamura T, Akaike H, Omori Y, Niimi Y, Takahashi S, Fujimaki A, Maekawa S, Otani Y. 2015 Quasiparticle-mediated spin Hall effect in a superconductor. Nature Materials 14, 675–678.
• 76 Wolf MJ, Hübler F, Kolenda S, v Löhneysen H, Beckmann D. 2013 Spin injection from a normal metal into a mesoscopic superconductor. Physical Review B 87.
• 77 Quay CHL, Dutreix C, Chevallier D, Bena C, Aprili M. 2016 Frequency-Domain Measurement of the Spin Imbalance Lifetime in Superconductors. Physical Review B, in press as a Rapid Communication.
• 78 Krishtop T, Houzet M, Meyer JS. 2015 Nonequilibrium spin transport in Zeeman-split superconductors. Physical Review B 91, 121407.
• 79 Bobkova IV, Bobkov AM. 2015, Long-Range Spin Imbalance in Mesoscopic Superconductors under Zeeman splitting, JETP Letters 101, 118-124.
• 80 Silaev M, Virtanen P, Bergeret FS, Heikkiä TT. 2015 Long-Range Spin Accumulation from Heat Injection in Mesoscopic Superconductors with Zeeman Splitting. Physical Review Letters 114, 167002.
• 81 Machon P, Eschrig M, Belzig W. 2013 Nonlocal Thermoelectric Effects and Nonlocal Onsager relations in a Three-Terminal Proximity-Coupled Superconductor-Ferromagnet Device. Physical Review Letters 110, 047002.
• 82 Ozaeta A, Virtanen P, Bergeret FS, Heikkilä TT. 2014 Predicted very large thermoelectric effect in ferromagnet-superconductor junctions in the presence of a spin-splitting magnetic field. Physical Review Letters 112, 057001.
• 83 Pottier N. 2007 Physique statistique hors d’équilibre : Processus irréversibles linéaires. Les Ulis: EDP Sciences.
• 84 Vier DC, Schultz S. 1983 Observation of conduction electron spin resonance in both the normal and superconducting states of niobium. Physics Letters A 98, 283–286.
• 85 Wiemann Y, Simmendinger J, Clauss C, Bogani L, Bothner D, Koelle D, Kleiner R, Dressel M, Scheffler M. 2015 Observing electron spin resonance between 0.1 and 67 GHz at temperatures between 50 mK and 300 K using broadband metallic coplanar waveguides. Applied Physics Letters 106, 193505.
• 86 Quay CHL, Weideneder M, Chiffaudel Y, Strunk C, Aprili M. 2015 Quasiparticle spin resonance and coherence in superconducting aluminium. Nature Communications 6, 8660.
• 87 Yamashita T, Takahashi S, Imamura H, Maekawa S. 2002 Spin transport and relaxation in superconductors. Physical Review B 65, 172509.
• 88 Takahashi S, Imamura H, Maekawa S. 1999 Spin Imbalance and Magnetoresistance in Ferromagnet/Superconductor/Ferromagnet Double Tunnel Junctions. Physical Review Letters 82, 3911–3914.
• 89 Bobkova IV, Bobkov AM. 2011 Recovering the superconducting state via spin accumulation above the pair-breaking magnetic field of superconductor/ferromagnet multilayers. Physical Review B 84, 140508.
• 90 Jin B. 2008 Response of spin-accumulated ferromagnet/superconductor/ferromagnet double tunnel junction to applied magnetic field. Journal of Magnetism and Magnetic Materials 320, 603–607.
• 91 Bhattacharjee S, Sardar M. 2000 Spin polarized carrier injection into high- superconductors: A test for the superconductivity mechanism. Physical Review B 62, R6139–R6142.
• 92 Koller D, Osofsky MS, Chrisey DB, Horwitz JS, Jr RJS, Stroud RM, Eddy CR, Kim J, Auyeung RCY, Byers JM, Woodfield BF, Daly GM, Clinton TW, Johnson M. 1998 Suppression of superconductivity by injection of spin-polarized current. Journal of Applied Physics 83, 6774–6776.
• 93 Vas’ko VA, Larkin VA, Kraus PA, Nikolaev KR, Grupp DE, Nordman CA, Goldman AM. 1997 Critical Current Suppression in a Superconductor by Injection of Spin-Polarized Carriers from a Ferromagnet. Physical Review Letters 78, 1134–1137.
• 94 Dong ZW, Ramesh R, Venkatesan T, Johnson M, Chen ZY, Pai SP, Talyansky V, Sharma RP, Shreekala R, Lobb CJ, Greene RL. 1997 Spin-polarized quasiparticle injection devices using Au/YBaCuO/LaAlO/NdSrMnO heterostructures. Applied Physics Letters 71, 1718–1720.
• 95 Yeh NC, Vasquez RP, Fu CC, Samoilov AV, Li Y, Vakili K. 1999 Nonequilibrium superconductivity under spin-polarized quasiparticle currents in perovskite ferromagnet-insulator-superconductor heterostructures. Physical Review B 60, 10522–10526.
• 96 Gim Y, Kleinsasser AW, Barner JB. 2001 Current injection into high temperature superconductors: Does spin matter? Journal of Applied Physics 90, 4063–4077.
• 97 Zheng Z, Xing DY, Sun G, Dong J. 2000 Andreev reflection effect on spin-polarized transport in ferromagnet/superconductor/ferromagnet double tunnel junctions. Physical Review B 62, 14326–14330.
• 98 Liu WV, Wilczek F. 2003 Interior Gap Superfluidity. Physical Review Letters 90, 047002.
• 99 Bedaque PF, Caldas H, Rupak G. 2003 Phase Separation in Asymmetrical Fermion Superfluids. Physical Review Letters 91, 247002.
• 100 Sheehy DE, Radzihovsky L. 2006 BEC-BCS Crossover in “Magnetized” Feshbach-Resonantly Paired Superfluids. Physical Review Letters 96, 060401.
• 101 Gubbels KB, Romans MWJ, Stoof HTC. 2006 Sarma Phase in Trapped Unbalanced Fermi Gases. Physical Review Letters 97, 210402.
• 102 Shin Y, Zwierlein MW, Schunck CH, Schirotzek A, Ketterle W. 2006 Observation of Phase Separation in a Strongly Interacting Imbalanced Fermi Gas. Physical Review Letters 97, 030401.
• 103 Casalbuoni R, Nardulli G. 2004 Inhomogeneous superconductivity in condensed matter and QCD. Reviews of Modern Physics 76, 263–320.
• 104 Shovkovy I, Huang M. 2003 Gapless two-flavor color superconductor. Physics Letters B 564, 205–211.
• 105 Shpielberg O, Akkermans E. 2015 Le Chatelier principle for out-of-equilibrium and boundary driven systems: application to dynamical phase transitions. arXiv:1510.05254 [cond-mat.stat-mech]
• 106 Zapata I, Sols F, Demler E. 2012 Triplet Pair Correlations in $s$-Wave Superfluids as a Signature of the Fulde-Ferrell-Larkin-Ovchinnikov State. Physical Review Letters 109, 155304.
• 107 Ho TL, Zhai H. 2007 Homogeneous Fermion Superfluid with Unequal Spin Populations. Journal of Low Temperature Physics 148, 33–41.
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters