Generalized quasiclassical theory of the long-range proximity effect and spontaneous currents in superconducting heterostructures with strong ferromagnets

We present the generalized quasiclassical theory of the long-range superconducting proximity effect in heterostructures with strong ferromagnets, where the exchange splitting is of the order of Fermi energy. In the ferromagnet the propagation of spin-triplet Cooper pairs residing on the spin-split Fermi surfaces is shown to be governed by the spin-dependent Abelian gauge field which results either from the spin-orbital coupling or from the magnetic texture. The additional gauge field enters into the quasiclassical equations in superposition with the usual electromagnetic vector potential and results in the generation of spontaneous superconducting currents and phase shifts in various geometries which provide the sources of long-range spin-triplet correlations. We derive the Usadel equations and boundary conditions for the strong ferromagnet and consider several generic examples of the Josephson systems supporting spontaneous currents.


I. INTRODUCTION
Effective gauge theories has been introduced in many condensed matter systems including spin-triplet superfluid 3 He 1 , cold atom systems 2,3 and magnetic materials 4,5 . In spatially inhomogeneous magnetic textures the additional spin-dependent U (1) gauge field of topological origin affects the motion of conduction electrons in the same way as the external electromagnetic field 4, [6][7][8][9] . That results in the topological Hall effect and emergent electrodynamics 5,10 observed recently in the chiral magnets with skyrmion lattices [11][12][13][14] .
Geometric flux associated with the spin-dependent gauge field was predicted to generate spontaneous spin and charge currents in mesoscopic rings with spatiallyinhomogeneous texture of the Zeeman field [15][16][17][18] . Up to now these effects have not been observed. Experimental detection of persistent currents in normal metals is in general rather challenging 19,20 since their magnitude is determined by the single-level contribution which is rather small and highly sensitive to the details of disorder potential 21,22 .
The situation is completely different in the superconducting state where the locally broken U (1) gauge symmetry leads to the Meissner effect, i.e. the generation of persistent condensate currents in response to the external magnetic field. However effects associated with the geometric spin-dependent flux [15][16][17][18] have not been identified in usual superconducting systems since the condensate of spin-singlet Cooper pairs is not sensitive to the Zeeman field rotations.
In the present paper we show that the superconducting condensate in fact can be coupled to the spin-dependent gauge fields emerging in superconductor/ferromagnet (SC/FM) hybrids. In such systems the interplay of superconducting and magnetic orderings results in the generation of the spin-triplet correlations through the proximity Gauge theory of the long-range proximity effect and spontaneous currents in superconducting heterostructures with strong ferromagnets.
We present the generalized quasiclassical theory of the long-range superconducting proximity effect in heterostructures with strong ferromagnets, where the exchange splitting is of the order of Fermi energy. In the ferromagnet the propagation of spin-triplet Cooper pairs residing on the spin-split Fermi surfaces is shown to be governed by the spin-dependent Abelian gauge field which results either from the spin-orbital coupling or from the magnetic texture. The additional gauge field enters into the quasiclassical equations in superposition with the usual electromagnetic vector potential and results in the generation of spontaneous superconducting currents and phase shifts in various geometries which provide the sources of long-range spin-triplet correlations. We derive the Usadel equations and boundary conditions for the strong ferromagnet and consider several generic examples of the Josephson systems supporting spontaneous currents.
Effective gauge theories has been introduced in many condensed matter systems including spin-triplet superfluid 3 He ? , cold atom systems ? ? and magnetic materials ? ? . In spatially inhomogeneous magnetic textures the additional spin-dependent U (1) gauge field of topological origin affects the motion of conductions electrons in the same way as the external electromagnetic field ? ? ? ? ? . That results in the topological Hall effect and emergent electrodynamics ? ? observed recently in the chiral magnets with skyrmion lattices ? ? ? ? .
Geometric flux associated with the spin-dependent gauge field was predicted to generate spontaneous spin and charge currents in mesoscopic rings with spatiallyinhomogeneous texture of the Zeeman field ? ? ? ? . Yet up to now these effects have not been observed. Experimental detection of persistent currents in normal metals is in general rather challenging ? ? since their magnitude is determined by the single-level contribution which is rather small and highly sensitive to the details of disorder potential ? ? .
The situation is completely different in the superconducting state where the locally broken U (1) gauge symmetry leads to the Meissner effect, i.e. the generation of persistent condensate currents in response to the external magnetic field. However in usual superconductors the condensate of spin-singlet Cooper pairs is not sensitive to the Zeeman field rotations. Therefore effects associated with the geometric spin-dependent flux ? have not been identified in superconducting systems.
In the present paper we show that the superconducting condensate in fact can be coupled to the emergent spin-dependent gauge fields in hybrid ferromagnet/superconductor (FS) structures consisting of the ordinary spin-singlet superconductors and strong ferromagnets which have substantial exchange energy splitting of spin subbands. In such systems the interplay of superconducting and magnetic orderings results in the generation of the spin-triplet correlations through the proxim-   To understand the behaviour of ESC we develop the gauge theory formalism to treat proximity effect in SC/FM systems with spin-textured strong ferromagnets. So far, proximity and transport calculations in SC/FM hybrids have mostly concentrated on either fully polarized systems, so called half metals [24][25][26][27][28] , or in the opposite limit of weakly polarized systems 23,29 , where the difference between spin subbands is completely neglected. However, most FMs have an intermediate exchange splitting of the energy bands of the order of but less than the Fermi energy. By now the quasiclassical theory for this regime has been formulated for the case of homogeneous magnetization of the strong ferromagnet 30 . To describe the general situation we go beyond those limitations and consider SC/FM structures with arbitrary large and spatially inhomogeneous exchange field.
Our approach relies on the adiabatic approximation for spin transport 31 which has been extensively used for studying transport phenomena in spin-textured magnets 4,5 . We derive the quasiclassical equations describing ESC interacting with the spin-dependent U (1) gauge field which can be induced either by the magnetic texture or spin-orbit coupling (SOC). The phases picked up by the S z = ±1 Cooper pairs in response to this gauge field generate spontaneous superconducting currents through strong FMs.

A. The Model
We consider the Gor'kov equations in the presence of spin-dependent gradient terms, the exchange field h and the external vector potential A: HereΠ =p − eAτ 3 ,Ǧ =Ǧ(r, r ) is the matrix Green function (GF) in spin-Nambu space, {, } is the anticommutator added to have the hermitian Hamiltonian in the system with space-dependent field M kj , µ is the chemical potential, ω is the Matsubara frequency, m is the effective mass which is equal to m F in the ferromagnet and to m S in the superconductor, e is the electron charge, p = (p x ,p y ,p z ) is the momentum differential operator, σ k andτ k are the spin and Nambu Pauli matrices. The self-energy termΣ includes the effects related to disorder scattering as well as the non-diagonal superconducting potential. The spin-dependent term M kj can be associated either with the SOC or with the pure gauge SU (2) field M kj = −iTr σ kÛ † ∇ jÛ /2m where the transfor-mationÛ (r) = e iσθ(r)/2 rotates spin axes to the local frame where h z. It is parametrized by the spin vector θ = θn defined by the spatial texture of the ex-change field distribution h(r) =R(θ(r))h, whereR is the spatially-dependent rotation matrix and we choose h = hz. Therefore, Eq. (1) is written in the local reference frame, where the quantization axis is aligned with the local exchange field. We assume that the exchange field rotates slowly, on the large scales as compared to the atomic distances. For this reason we neglect secondorder spatial derivatives of the exchange field. In the framework of this approach the inhomogeneity of the exchange field enters the equations as the pure gauge SU (2) field.
In the general case of large exchange splitting |h| ∼ µ the spin-dependent Gor'kov Eq. (1) is rather complicated and most importantly one cannot apply here the quasiclassical theory. The quasiclassical approximation is violated by the mixed-spin correlations (MSC) residing in spin-split subbands [ Fig. 1b] which are characterized by the spatial length scale of the order of the Fermi wavelength λ F = 2m/µ. Therefore MSC yield vanishingly small contribution to the momentum-averaged observables at the distances much larger than the atomic length from the FM/SC interface. Such correlations can be incorporated to the effective boundary conditions as the source terms for the ESC [ Fig. 1a]. The ESC survive in the ferromagnet at much larger distances and can be treated within quasiclassics considered separately for each of the spin-up and spin-down Fermi surfaces.

B. Adiabatic approximation
To develop the quasiclassical approximation we divide the GF into the parts corresponding ESC of the spinup/down states (see Fig.1a) and the one corresponding to the MSC (see Fig.1b) Then from the Gor'kov equation (1) one can see that the amplitude of MSC is in general proportional tǒ Therefore MSC are small as compared to that of the ESC if the adiabatic criterion is satisfied |M ij /λ F h| 1. Within the adiabatic approximation neglecting the MSC in Eq.(1) and substituting the expansion of momentum operatorp 2 = p 2 − 2ip∇, we obtain the quasiclassical equation for ESC part: Here v ± = 2(µ ± h)/m F are the spin-dependent Fermi velocities determined on each of the spin-split Fermi surfaces labelled by the subscript σ = ±. We introduce the projection operators to spin-up and spin-down stateš γ + =τ ↑σ↑ +τ ↓σ↓ andγ − =τ ↑σ↓ +τ ↓σ↑ , respectively, The paired states on each of the Fermi surfaces are given by the corresponding parts of the equal-spin correlator:Ǧ ± =γ ±ǦES . This decomposition allows to introduce quasiclassical propagators separately for spinup and spin-down blockŝ where ξ pσ = p 2 /2m F + σh − µ and the notation means that the integration takes into account the poles of GF neat the corresponding Fermi surface. From Eq.(5) we obtain generalized Eilenberger equations fro the spin-less quasiclassical propagators where the covariant operator iŝ One can see that the Eilenberger-type equations for the spin-up/down correlations contain an additional U(1) gauge field Z which is added to the usual electromagnetic vector potential A with the opposite effective charges for spin-up and spin-down Cooper pairs. The U(1) field is obtained by projecting the initial SU(2) field σ k M ki to the basis of spin-triplet pairing states: Z = −i(Û −1 ∇Û ) 11 . This reduction means that we neglect spin-flip transitions between the spin-up and spin-down Cooper pairs induced by the SU(2) potential. On a qualitative level it is equivalent to the adiabatic approximation in the single-particle problems that allows to describe the quantum system evolution in terms of the Berry gauge fields 31 . Finally, the quasiclassical expression for the charge current is given by where ν σ are the spin-resolved DOS and .. denotes the averaging over the spin-split Fermi surface.

C. Usadel equation for ESC
Let us consider the system with large impurity scattering rate as compared to the superconducting energies determined by the bulk energy gap ∆. In this experimentally relevant diffusive limit it is possible to derive the generalized Usadel theory with the help of the normalization conditionĝ 2 σ = 1 which holds due to the commutator structure of the quasiclassical equations (8).
The impurity self-energy in the Born approximation is given byΣ σ = ĝ σ /2iτ σ . In the dirty limit we have The solution of Eq. (11) can be found asĝ σ = ĝ σ + g a σ pσ pσ , where the anisotropic part of the solutionĝ a σ is small with respect to ĝ σ . Making use of the relation { ĝ σ ,ĝ a σ } = 0, which follows from the normalization condition, one obtainŝ Substituting to Eq. (7) and omitting the angle brackets we get the diffusion equation where D σ are the spin-dependent diffusion coefficients, in the isotropic case given by D σ = τ σ v 2 σ /3. The current is obtained substituting expansion (12) to the Eq.(10) Equations (13,14) together with the boundary conditions derived in the next section provide the framework to study proximity effect-related phenomena in strong ferromagnets with large amount of disorder. This regime describes a different physical situation as compared to the previous works, where other approximations have been used. That concerns papers studying Rashba superconductor with large SOC 32 and that dealing with weak exchange fields and SOC within SU(2)-covariant formulation [33][34][35][36] . The applicability condition of our theory requires that the exchange splitting has to be much larger than all other energy scales except the Fermi energy. In particular, in the bulk superconductor we have omitted the spin-flipping terms produced by various sources like SOC, magnetic texture or magnetic impurity scattering. However the spin-flipping terms are still important in strong ferromagnets/half metals since they provide the source of ESC. Such correlations appear due to the conversion of spin-singlet Cooper pairs leaking from the superconductor electrode into the longrange spin-triplet ones. Below we study this conversion and boundary conditions for ESC using a generic model of SC/FM interface.

III. BOUNDARY CONDITIONS
The quasiclassical Eq. (7) deals with the long-range ESC transformed adiabatically in space under the action of spin gauge fields. In FM/SC systems with usual spinsinglet superconductors such correlations can appear only in result of the non-adiabatic spin flip which converts MSC into the ESC 26,28,37 . This process occurs within the thin layer near FM/SC interface and can be described by the effective boundary conditions. To derive them we consider the simple microscopic model of the boundary with the non-magnetic potential barrier. The FM/SC interface is located at x = 0, the normal n = n x x directed from SC to FM and n x = ±1. The momentum components p in yz plane parallel to the boundary are conserved. The effective masses m S and m F in the superconducting and ferromagnet regions are assumed to be different and the FM/SC boundary is characterized by the interfacial potential barrier of the strength V δ(x) added to the HamiltonianĜ −1 0 in Eq.(1). Let us outline the general strategy to deriving the boundary conditions for quasiclassical ESC propagators. We need to solve the exact Gor'kov equations (1) near the boundary with accuracy up to the first order in SU (2) terms, which provide the conversion of MSC to ESC. Therefore we will use an expansion by the small parameter |M ij p j /h| 1. In result we will find the slow component of the anomalous functionF (x, x ) = (Ǧ) 21 , where the index corresponds to the Nambu space. This component does not contain fast oscillations as a function of the center of mass coordinate X = (x + x )/2. Let us denote such components asF(x, x ) . Below we will show that these correlations have the form The spin-upF + and spin-downF − pairing amplitudes are given bỹ where p σn = 2m F (µ + σh) − p 2 || , and S σ⊥ = (M xi + iσM yi )n i /h is the combination of SU(2) field components that generate ESC in the ferromagnet near the superconducting interface. Here the coefficient K σ incorporates the dependence of the pairing amplitude on the interface barrier strength, order parameter and effective masses.
By writing Eq.(16) we assume the averaging over the directions of the in-plane momentum. This is enough in the dirty limit although in the clean case additional important effects resulting from the in-plane gradients of the exchange field h can be obtained beyond this approximation. Eq.(16) is valid for ω > 0 and for ω < 0 the amplitude can be obtained using symmetry relations, as discussed below.
We need to find the GF near the FM/SC interface determined by the 1D Gorkov' equation (1) along x coordinate. Since each GF is the 2x2 matrix in spin space they can be represented in the form of two spinorŝ where the spinor elements areû k = (u k1 , u k2 ) T andv k = (v k1 , v k2 ) T . Let us consider the componentsû 1 ,v 1 in detail. The other pairû 2 ,v 2 is given byû 2 = σ xû1 and v 2 = σ xv1 and changing the sign of the fields h and M yi . For the Nambu spinorψ = (û 1 ,v 1 ) T from (1) we obtain the equation in the ferromagnet and in the superconductor where µ S,F = µ − p 2 /2m S,F . These equations look similar to the Bogolubov-de Gennes equation but taken at the imaginary frequency. The boundary condition are obtained by integrating Eq.(1) near the singularity points x = x and x = 0. In this way we obtain the boundary conditions at x = x : ant at x = 0 Here we neglect the impurity scattering and the inverse proximity effect in the superconducting region. The impurity self energy can be included into the consideration in case of the tunnelling limit when the surface barrier is strong enough to suppress the inverse proximity effect in the superconductor. This consideration demonstrates that the impurity scattering does not change boundary conditions for quasiclassical propagators. In the opposite case of weak barriers one should strictly speaking take into account how the impurity self energy in the superconductor is modified by the inverse proximity effect, which in principle can affect the boundary conditions.
Let's assume that n x < 0. Then the solution for electron wave in bulk normal ferromagnet can be written aŝ v = 0 and where the labels > (<) denote right-and left-going waves, p σn = 2m F (µ F + σh + iω) and we neglect SO corrections to the momenta. Expressions (26) are valid for the half-metal when h > µ F as well. In this case one should take p −n = i 2m F (h − µ F ). To the first order by the small parameter |M ij p j /h| 1 we get expressions for the amplitudes in Eq.(26) where S σ = (M xy +iσM yy )/h and S σ⊥ is defined above. The zero -order amplitudes are given by We neglect the second term in Eq.(26)Â −> since its amplitude is much smaller than a (0) +> due to the prefactors p σn S σ⊥ and p y S σ . Therefore up to the first order in these small parameters we should take into account the reflected holes generated byÂ −> without spin flip which have the same wave vector as the incident wave and therefore does not contribute to the slow-varying correlation (15).
The solution (26) can be considered as the incident electronic wave at the FM/SC interface. The reflected wave consists of electronicû r and holev r components having the form Here the structure ofD σ , is similar to that ofÂ σ< The reflected hole-like wave is given by (34) with the amplitudeŝ We are interested in the waveB − because the corresponding contribution to the reflected hole amplitudê v r given by the second term in Eq.(34) does not contain fast oscillations as the function of X = (x + x )/2. Hence it provides the source of the long-range superconducting correlations. On a qualitative level the waveB − determines the spin-flip Andreev reflection leading to the generation of spin-triplet Cooper pairs.
Hence we obtain the component of the ESC in the form . The other component is obtained from the other pair of spinors (û 2 ,v 2 ) and has the formF − (x, x ) = b + e i(p−nx−p * −n x ) . Averaging over in-plane momentum directions we get bσ = S σ⊥ K σ in accordance with the Eq.(16) where we put n x = −1.
General expressions for the amplitudes K σ are rather involved (see Appendix A). However in the tunnelling limit they read where F bcs = ∆/ ω 2 + |∆| 2 and v Sn = 2µ S /m S . For the half-metal Eq.(39) is valid with the imaginary mo- The other Nambu component of the anomalous functionF (x, x ) = (Ǧ) 12 can be obtained from the general particle-hole symmetrŷ Now having in hand the expression for the slowlyvarying amplitude (16) we can derive the boundary conditions for the components of quasiclassical propagator. Following Ref. 38 we write them in the form whereτ ± = (τ x ± iτ y )/2. The quasiclassical propagators can be obtained by taking the Fourier transform of the slow-varying exact GF components (16) and then using the definition (6). In this way we obtain the propagators as functions of the momentum directionp σ , which determines the quasiclassical trajectories in each of the spin subbands. We will use the notationsf σ,in(out) for the "incoming" ("outgoing") trajectories withp σ n < (>)0 . Then Eq. (16) yields the quasiclassical ESC propagators at the interface x = 0: where v σn = p σn /m F is the spin-dependent Fermivelocity. At the same timef σ,out (x = 0, ω > 0) = 0.
The advantage of the boundary conditions (42) is that they give an explicit value of the Green's function at the ferromagnetic side of the interface. The price, which we have paid for it, is that, strictly speaking, they are only valid for an isolated interface, because the asymptotic conditions at the infinity were essentially used in the derivation. However, they can be safely applied to the dirty systems with more than one interface if the distances between the interfaces are large as compared to the mean free path. Below we are only interested in the dirty case.
The boundary conditions to the Usadel equations can be obtained from Eq. (42) in a straightforward way. As usual, one can show 40 that in the isotropization region near the interface the matrix current is where ... −(+) means the averaging over the part of the ferromagnet FS corresponding to p σ n < (>)0 and real values of p σn and p sn , l σ = v σ τ σ is the mean free path and∂ n = n∂ R . The boundary condition to the Usadel equation is obtained using (43), taken at x = 0 witĥ g σ,in(out) from Eq. (42) and the symmetry relations discussed above: where κ σ = −i K σ v 2 σn /v σ − . From the boundary condition (44) one can see that the generation of ESC is determined by the non-adiabatic spin-flipping terms near the boundary. In case if these terms are of the SOC origin, e.g. having the Rashba form the magnitude of ESC correlations is given by S σ⊥ ∼ α/h, where α is the SOC constant. Otherwise if it comes from the magnetic texture with the characteristic scale ξ θ the estimation is S σ⊥ ∼ 1/(m F ξ θ h). However as shown below, that smallness affects only the overall amplitude of the critical current in the generic SC/FM/SC Josephson systems , but not the spontaneous phase shift of the current-phase relation. Indeed the emergent gauge field Z which drives the spontaneous supercurrents through strong ferromagnets does not contain any small parameter. Therefore the anomalous current at zero phase difference across the junction can be of the order of the critical one.

IV. SPONTANEOUS JOSEPHSON CURRENT THROUGH STRONG FERROMAGNETS
Having in hand the machinery of the generalized quasiclassical theory described above we can calculate Josephson current-phase relations for different systems with spin-dependent fields. Example of such a system with magnetic helix texture is shown in Fig. 1c .
Here we work in the dirty limit using the linearized (with respect to the anomalous Greens function) version of the spin-less Usadel equations (13) and boundary conditions (44). This simplification is adequate if the proximity effect at the SC/FM interface is weak, for example when the interface is low-transparent. The absolute value of the order parameter is assumed to be the same in the superconducting leads, while there is the phase difference χ between them. The electric current in the ferromagnetic interlayer can be calculated according to Eq. (14), which in the linearised form is reduced to The anomalous GF f σ andf σ should be calculated from the linearized version of the Usadel equation (13): The solution of these equations takes the form: The coefficients C σ,± andC σ,± are to be found from the boundary conditions (44) taken at the S/F interfaces x = ∓d/2. The resulting expressions take the form: andC where λ σ = 2|ω|/D σ . Substituting the anomalous GF from (47) with the coefficients C σ,± andC σ,± from (48) and (49) into Eq. (45) we obtain the general currentphase relation (CPR): where χ is the Josephson phase difference R σ = 1/(e 2 ν σ D σ ) is the spin-resolved resistivity and λ σ = 2|ω|/D σ . The spin gauge field Z x = 0 and finite spin splitting D + = D − in the CPR (50) lead to the spontaneous current at zero phase difference known as the anomalous Josephson effect. The ground state phase difference χ 0 can be found from the zero-current condition I(χ = χ 0 ) = 0 The spontaneous phase shift of Josephson current has been obtained in several FM/SC systems 25,27,30,35,[41][42][43][44][45][46][47][48][49][50][51][52][53][54] .
Here we demonstrate that this effect is essential only for the case of strong ferromagnets. When the ferromagnet is weak and treated within the usual quasiclassical approximation, the difference between I − and I + is neglected and the anomalous Josephson effect disappears.
The spin gauge field Z is generated by the spin helix shown schematically in Fig.  1. Recently the proximity effect in helical magnets has been observed experimentally 55,56 . In this case the magnetization texture is described by h = h(cos α, sin α cos θ, sin α sin θ), where we assume that the angle α is spatially independent and θ = θ(x).
The general theory developed above describes the proximity effect in homogeneous ferromagnet with a linear in momentum SOC 33,34 . For example let us consider the SC/FM/SC junction through the quasi-2D ferromagnet in xz plane, interfaces in yz planes and the exchange field in the plane of the ferromagnet h z . In case of the Rashba SOC in the ferromagnetic region this system is characterized by the spin-dependent fields M zx = −M xz = −α/2 which leads to while the Dresselhaus SOC yields M zz = −M xx = β/2 and therefore In each of these cases the Josephson CPR can be found substituting the fields into the general Eqs. (50,51). Since the ground state phase shift χ 0 is determined by the component Z x parallel to the Josephson current, for h z we have χ 0 = 0, π only for the Rashba SOC, but not for the Dresselhaus one. This is natural, because in general case of a magneto-electric effect the spontaneous current and the magnetization are perpendicular to each other for the Rashba SOC 41,57,58 , but they are parallel for the Dresselhaus SOC. Therefore, in order to get the anomalous Josephson current for the Dresselhaus SOC, h is to have a component parallel to the current. Comparing Eqs. (51) and (53,54) one can see that in the considered geometry the Dresselhaus SOC produce the long-range ESC even in the case of the homogeneous magnet (while in general the both Rashba and Dresselhause SOC can produce ESC 34 ). However their amplitudes are determined by the SOC constants which in general are rather small in metals. Although the anomalous phase shift is also determined by SOC, it can become rather large for sufficiently long junction, i.e. when |Z x d| ≥ 1. Therefore even the weak SOC leads to the significant phase shifts of the CPR although the overall critical current amplitude is rather small . In reality however the long-range ESC can be generated by the magnetic inhomogeneity near the interface 26 . Let us consider the following model : h = h(sin θ, 0, cos θ), U = e −i(θ/2)σy , (55) where θ = θ(x) changes linearly in the region ξ θ d near the interfaces and θ = 0 in the bulk FM. Neglect the effect of SOC we get S σ⊥ = −iσn x ∂ x θ/(2m F h) in the CPR Eq. (50), where the spin gauge field Z x is determined by the SOC in the bulk FM. For small exchange fields h µ the anomalous phase shift χ 0 is determined by the prefactor (I + − I − )/(I + + I − ) ∼ h/µ. If we assume additionally (in general it can be not the case) that Z x d 1, then χ 0 ∼ m F (hd/µ)α, in agreement with Ref. 41.

V. CONCLUSION
To conclude, we have developed the generalized quasiclassical formalism to calculate the behaviour of longrange ESC in the ferromagnet. These correlations can be generated at the ferromagnetic/superconductor interface in the presence of either magnetization inhomogeneity or SOC. The general conditions for ESC generation are derived in terms of the SU(2) gauge fields. In the ferromagnetic material the behaviour of ESC is shown to be governed by the adiabatic spin gauge field which generates spontaneous superconducting currents through strong FMs with magnetic texture or SOC. These results demonstrate that spontaneous superconducting currents exist as a robust and experimentally observable phenomenon in many superconducting/ferromagnetic systems studied in connection to the superconducting spintronics 59,60 .

VI. ACKNOWLEDGEMENTS
We thank A.S. Mel'nikov for stimulating discussions. This work was supported by the Academy of Finland. I.V.B. and A.M. B. acknowledge the support by the Program of Russian Academy of Sciences Electron spin resonance, spin-dependent electron effects and spin technologies and Russian-Greek Project N 2017-14-588-0007-011 Experimental and theoretical studies of physical properties of low-dimensional quantum nanoelectronic systems.

Appendix A: Spin-flip Andreev reflection coefficient
Here we find the reflection coefficients of the electronic and hole-like waves in Eqs. (33,34). For this purpose we use Eq.(26) as the incident wave coming from the ferromagnet to the FM/SC interface. To apply the boundary conditions (24,25) we write solution in the superconductor, as the superposition of two terms decaying at x → ∞ : where q 1 = 2m S (µ S + iΩ) and q 2 = − 2m S (µ S − iΩ) and Ω = ω 2 + |∆| 2 . Below we will neglect the imaginary part of q 1,2 and use q 1 ≈ /q 2 ≈ p Sn where p Sn = √ 2m S µ S .
First, let us find reflection coefficients without spin-flip, in zero order by small parameter M + x,y p x,y /h. For this purpose we obtain the following system of equations : where The solutions are a (0) The spin-flip reflection amplitude b − can be found taking into account first-order corrections in M + x,y p x,y /h when matching the electron and hole waves at FM/SC boundary. In this way we get the linear system where the coefficientsα = α 0 −M kxσk andα = α * 0 +M kxσk take into account the correction to the boundary condition from the effective spin-orbital term. The spinorsÂ The solution of this system readsB wherê Z = Z + + 4iΩV M kxσk (A20) where v Σ = v +n + v −n and v d = v +n − v −n . In this way we obtain the spin-flip Andreev reflection amplitude in the form b − = K x+ S +⊥ + K y+ S + with (ii) Let us now consider the other pair of spinors in Eq.(17)û 2 ,v 2 which determine the correlation functionF − . They can be obtained fromû 1 ,v 1 by transforming the HamiltonianȞ → σ xȞ σ x which flips the spin index σ. The reflected hole-like states are given by the same Eq.(37), but this time we are interested in the waveB + (x ). This wave has the form (37) with the amplitude b + = K x− S −⊥ + K y− S − , where where the amplitudes d (iii) Limiting cases: Consider large barrier V v S , v σ , v σ , h/p σ so that α ≈ iV and Z 0 ≈ Y ≈ ΩV 2 . Then from Eqs. we get d Substituting to the Eq.(A23) and taking into account (A9) we get These expressions are valid for the half-metal when h > µ F as well. In this case one should take p −n = i 2m F (h − µ F ).