Spontaneous Edge Accumulation of Spin Currents in Finite-Size Two-Dimensional Diffusive Spin-Orbit Coupled SFS Heterostructures

We theoretically study spin and charge currents through finite-size two-dimensional $s$-wave superconductor/uniform ferromagnet/s-wave superconductor (SFS) junctions with intrinsic spin-orbit interactions (ISOIs) using a quasiclassical approach. Considering experimentally realistic parameters, we demonstrate that the combination of spontaneously broken time-reversal symmetry and lack of inversion symmetry can result in spontaneously accumulated spin currents at the edges of finite-size two-dimensional magnetic SF hybrids. Due to the spontaneous edge spin accumulation, the corners of the F wire host the maximum spin current density. We further reveal that this type edge phenomena are robust and independent of either the actual type of ISOIs or exchange field orientation. Moreover, we study spin current-phase relations in these diffusive spin-orbit coupled SFS junctions. Our results unveil net spin currents, not accompanied by charge supercurrent, that spontaneously accumulate at the sample edges through a modulating superconducting phase difference. Finally, we discuss possible experimental implementations to observe these edge phenomena.


I. INTRODUCTION
Spintronics devices operate by spin transport mechanisms 1-6 rather than by utilizing charged carriers, as is done typically in conventional electronics devices. The use of spin currents can result in higher speeds and reduced dissipation 2,6 while exhibiting weak sensitivity to nonmagnetic impurities and temperature. [4][5][6] For functional spin-based devices, it is necessary to manipulate and generate spin-currents in a practical and efficient manner. For this reason, many investigations have focused on harnessing the spin-orbit interactions 12,13 (SOIs) present in many materials, including semiconductors. [7][8][9][10][11][14][15][16][17] The SOI is a quantum relativistic phenomenon that can be divided into two categories: i) intrinsic (originating from the electronic band structure of the material) and ii) extrinsic (originating from spin-dependent scattering of impurities). [4][5][6] The intrinsic spin-orbit interactions (ISOIs) such as Rashba 12 and Dresselhaus 13 , are experimentally controllable via tuning a gate voltage 16,[18][19][20][21][22][23][24] . This particular attribute has proliferated efforts striving for high-performance spin-based devices, including transistors, and new routes in information storage and transport. 3,15,16,20,[25][26][27] Similarly, ferromagnet (F ) and superconductor (S) heterostructures have received renewed interest lately due to the possibility of generating spin polarized triplet supercurrents [28][29][30][31][32][33][34][35][36][37] that can be used for practical purposes 32 . By considering a ferromagnet with an ISOI, the spin orbit interaction can couple with the magnetic exchange field, resulting in modified superconducting proximity effects and additional venues for new spin phenomena. Indeed, the ISOI can induce long-range proximity effects in uniformly magnetized S/F structures due to the momentum-dependence of the effective exchange field 73 . It is therefore of fundamental importance not only to find a simple, experimentally accessible structure that can support spin currents in F/S systems, but it is also crucial to determine the spatial behavior of the spin currents near the boundaries of the superconducting hybrids.
Many past works are based on the application of external electric or magnetic fields. One of the earliest such cases involved the combination of SOIs and an external electric field, giving rise to an accumulation of spin currents at the edges 41,42 (the so called spin-Hall effect 39,40 ). The spin currents generally tend to peak near the sample boundaries and vanish at the electrode/sample interfaces. [44][45][46][47] These theoretical predictions were later experimentally observed in semiconductor samples 43 . The spin-Hall phenomenon was also extensively studied in superconducting heterostructures where various types of spin-orbit coupling (SOC) play key roles. 23,[50][51][52][53][54][55][56][57][58][59][60][61][62][63][64][65][66][67] For example, the out-of-plane component of the spin density was theoretically investigated 50 in a spin-orbit coupled S/N/S junction [with normal metal (N ) interlayer] subject to an inhomogeneous external magnetic field. It was found that the spin density varies along the transverse direction, leading to a longitudinal phase difference between the S electrodes. The influence of extrinsic SOIs on the critical supercurrent in diffusive magnetic hybrid structures was also studied. 71,72 In nonmagnetic S/N /S Josephson junctions with SOC subject to a magnetic field, 0-π transitions may be induced. 29 In an earlier work 74 , singlet-triplet pair conversion was numerically investigated using a lattice model in a ballistic half metal ferromagnetic Josephson junction with an interfacial Rashba SOC. Several optimal configurations have also been theoretically studied for generating and detecting the predicted longrange triplet correlations in experiments. 73 In this paper, we study the local spin currents in uniformly magnetized S/F /S Josephson junctions with spatially uniform intrinsic SOIs, avoiding any external electric or magnetic fields. We employ a two-dimensional quasiclassical Keldysh-Usadel approach that incorporates a generic spin-dependent vector potential to study the behavior of the spin current components. We consider two types of ISOCs: i) Rashba and ii) Dresselhaus SOC, and the magnetization of the F wire can take arbitrary orientations. We find that the coupling of the F wire's exchange field and ISOIs leads to edge spin currents with three nonzero components. The spin current density components peak near the edges of the F strip and sharply decline when moving away from the edges. Therefore, the maximal spin current accumulation takes place near the F wire's corners. This phenomenon can be also observed in ISO coupled S/N /S junctions with a single spin active N /S interface as demonstrated in Ref. 75. Our results show that the spin current can be switched on or off at the S/F contacts, depending on the magnetization direction. The spatially averaged spin current components reveal a 2π-periodicity and even-functionality in ϕ, the phase difference between the S terminals. This is in contrast to the charge supercurrent which is a 2π-periodic odd function of ϕ (and consistent with previous studies 68 ). Note that for such junctions, the argument in the current phase relation for some situations can become modified by a simple ϕ 0 shift. 66,69,70 The simple hybrid structure proposed here relies solely on the intrinsic properties of the system, in contrast to other structures that rely inextricably on external fields to observe the edge spin currents [39][40][41][42][43][44][45][46]49,50,67 . Our device consists of a finitesize intrinsic SO coupled F wire (with uniform magnetization) sandwiched between two S banks. The spin currents then spontaneously accumulate at the sample edges, without the application of an external electric or magnetic field to the system. We demonstrate that the device allows for the realization of spin currents in the absence of charge supercurrent by simply modulating ϕ. The edge spin accumulation is a signature of the spin Hall effect 67 , and hence can be experimentally measured by optical experiments for instance 43 (see the discussions in Sec. III). Also, we discuss the symmetries present among the spin current components when varying the magnetization orientation with Rashba or Dresselhaus SOC present. Moreover, we compare our results with the charge and spin currents found in a nonmagnetic diffusive S/N /S Josephson junction with Rashba and/or Dresselhaus SOC. We find that the spin currents vanish in the S/N /S devices, consistent with previous works 67,75 , and that the charge current displays a spatially uniform profile without any transverse component, indicating conservation of charge current.
The paper is organized as follows. We outline the theoretical framework used to study hybrid structures with ISOCs in Sec. II. In Sec. III, the results of diffusive S/F /S Josephson junctions are presented for the case of Rashba ISOC. We next utilize the symmetries in the spin currents to give a simple prescription for finding the corresponding results for the Dresselhaus spin-orbit interaction. We finally present concluding remarks in Sec. IV.

II. THEORETICAL FORMALISM
The intrinsic SOI is a consequence of the moving carriers' spin interaction with an atomic potential V (r). Therefore, the total Hamiltonian of a moving electron in such an atomic potential can be expressed as, 20,76 where m 0 = 0.51Mev, is the free electron mass and and c is the velocity of light in vacuum. We define P to represent the electron's momentum vector, and ∂ ≡ (∂ x , ∂ y , ∂ z ). The vector of Pauli matrices, τ , is given in Appendix A. It has been shown that the linearized SOC term can be simply accounted for as an effective background field that follows SU(2) gauge symmetries. 47,73,77,78 Hence, it is sufficient to replace partial derivatives, appearing in the quasiclassical formalism, by their covariants. 47,73 Another advantage of the SU(2) approach is the convenient definition of physical quantities such as spin currents. 79 We start with the Usadel equations that enable us to study the charge and spin transport through diffusive F/S systems with the ferromagnetic regions having arbitrary magnetization patterns h(r) = h x (r), h y (r), h z (r) : 28,73,80 where ρ and σ denote vectors comprised of 4 × 4 and 2 × 2 Pauli matrices (see Appendix A), and D represents the diffusive constant of the ferromagnetic medium. We have denoted the quasiparticles' energy by ε which is measured from the Fermi surface ε F . Throughout this work, we focus on the low proximity limit of the diffusive regime 28 . In this limit, the normal and anomalous components of the Green's function can be approximated by, F no (r) 1 and F (r) 1, respectively. Thus, the advanced component of total Green's function,Ĝ(r), takes the following form: where each entry stands for a 2×2 matrix block. Considering the Taylor expansion, the advanced component can be given by:Ĝ Here we restrict our calculations to the equilibrium situations where the Retarded and Keldysh blocks of total Green's function are obtained by:Ĝ A (r) = − ρ 3Ĝ R (r)ρ 3 † , and Here, k B and T denote the Boltzmann constant and system temperature, respectively. The Usadel equation, Eq. (2), leads to sixteen coupled complex partial differential equations in the low proximity limit that become highly complicated with the presence of intrinsic SOI terms. Unfortunately, the resultant system of coupled differential equations can only be simplified and decoupled under very limiting conditions, 28,29 leading to analytical results. However, for the systems considered in this paper, numerical methods are the most appropriate, and often the only possible routes to investigate the relevant transport properties. 73 The differential equations must be supplemented by the appropriate boundary conditions to properly capture the transport characteristics of S/F /S hybrid structures. We thus employ the Kupriyanov-Lukichev boundary conditions at the S/F interfaces 81 and control the intensity of induced proximity correlations using the barrier resistance parameter, ζ: wheren is a unit vector, directed perpendicular to a given interface. The solutions to Eqs. (2) for a bulk, even-frequency s-wave superconductor reads, , (6) in which, is defined in terms of the superconducting gap ∆. Here the macroscopic phase of the bulk superconductor is denoted by θ, so that the difference between the macroscopic phases of the left and right S electrodes are given by θ l − θ r = ϕ. For more compact expressions in our subsequent calculations, we define the following piecewise functions: where Θ(x) stands for the usual step function. It is clear that the general boundary conditions given by Eq. (5) do not permit current flow through the hard wall boundaries of the finite-size two-dimensional S/F /S Josephson junction, shown in Fig. 1.
To study the influence of differing types of ISOI on the system transport characteristics, we adopt a spin-dependent tensor vector potential A(r) = A x (r), A y (r), A z (r) , as follows: 47,73,75,77,78 FIG. 1. (Color online) Schematic of a finite-size two-dimensional magnetic S/F /S Josephson junction. The superconducting electrodes and rectangular ferromagnetic nano-wire are labelled S and F , respectively. We assume that the quasiparticle current experiences an intrinsic spin-orbit interaction (ISOI) solely inside the F region.
The thickness and width of the ferromagnetic strip are labeled dF and WF , respectively. The junction is located in the xy plane and the S/F interfaces are along the y axis. The F region has a uniform exchange field denoted by h and can take arbitrary orientations h x , h y , h z .
Using the above vector potential, we define the covariant derivatives by;∂ Accordingly, the brackets seen in the Usadel equation, Eq. (2), and the boundary conditions, Eq. (5), (as well as the charge and spin currents that shall be discussed below, Eqs. (10) and (11)) take the following form: The spin and charge currents are key quantities that lend insight into the fundamental system transport aspects that provide valuable and crucial information for nanoscale elements in superconducting spintronics devices, as described in the introduction. Under equilibrium conditions, the vector charge ( J c ) and spin ( J sγ ) current densities can be expressed by the Keldysh block as follows: 47,78 where J c 0 = N 0 eD/4, J s 0 = J c 0 /2e, and N 0 is the number of states at the Fermi surface. The vector current densities determine the local direction and amplitude of the currents as a function of coordinates inside the F strip. In other words, J(r), provides a spatial map to the currents inside the system. We designate γ = x, y, z for the three components of spin current, J sγ . The matrices we use throughout our derivations are given in Appendix A. To obtain the total Josephson charge current flowing through the magnetic strip, an additional integration over the y direction should be performed on Eq. (10) (see Fig. 1). The spin-dependent fields yield lengthly and cumbersome expressions, the details of which are not presented here for clarity. Having now outlined the theoretical approach utilized in this paper, we can now present our findings in the next section.

III. RESULTS AND DISCUSSIONS
In our computations below, we consider a uniform and coordinate-independent vector potential, A(r), i.e. ∂ · A(r) = 0, so that the spin vector potential is constant in the entire F region. A specific choice for the constant spin vector potential that results in Rashba (α) 12 and Dresselhaus (β) 13 types of SOC is, By substituting the above set of parameters into Eqs. (7), we arrive at, The Rashba SOI 20 can be described through spatial inversion asymmetries while the Dresselhaus SOI 13 is described by bulk inversion asymmetries in the crystal structure. 20,21 Crystallographic inversion asymmetries 87 or lack of structural inversion symmetries 16,17,84,87 in heterostructures may cause the ISOIs (Color online) Spatial profile of the spin current in a uniformly magnetized S/F /S Josephson (see Fig. 1) junction without ISOI. The magnetic exchange field is oriented along z, h = (0, 0, h z ), and therefore, solely the z component of spin current J sz x (x, y) is nonvanishing. The junction length and width are set to dF = 2.0ξS and WF = 2.0ξS, respectively. The top panel exhibits the spin current variations along the x-position (the junction length) at four differing locations along the junction width, y = 0.25ξS, 0.5ξS, 0.75ξS, 1.0ξS. The bottom panel shows J sz x (x, y) as a function of y-position along the junction width, at x = 0.25ξS, 0.5ξS, 0.75ξS, 1.0ξS. considered here. For example, strain can induce such inversion asymmetries 43,87,92,93 and thus, ISOIs, or the adjoining of two differing materials may generate the requisite interfacial SOIs 16,17,73,84,87 . Nonetheless, there is no straightforward method to measure SOIs in a hybrid structure. One possible approach would be first principle calculations 85 in conjunction with spin transfer torque experiments 73,86,87 . The intrinsic SOIs are often given by the first-order quasiparticle momentum, which is locked to their spins. This linearized approach is a simplification to the more generic picture dealing with higher orders of momentum, 21,82,[92][93][94][95] which can be observed in e.g., engineered materials. 92, 93 We here assume that ISOCs can be described by linear terms in the carriers' momentum. 13,20 Candidate materials to support spontaneous broken time-reversal and broken inversion symmetries include electron liquids with ISOIs, which naturally tend to have a Stoner-type magnetism at low densities, and a magnetically doped topological insulator surface (or by directly coating a topological insulator surface with magnetic insulators). [88][89][90] Other promising candidates involve the ferromagnetic semiconductors (Ga,Mn)As, where both the electronic structure and inherent magnetism make these materials well suited for experimental studies. 7,8,91 Our quasiclassical approach allows us to study systems involving nontrivial magnetizations and spin vector potentials with arbitrary spatial patterns. We thus consider a finite-sized, uniformly magnetized F wire whose exchange field can take arbitrary orientations. In order to determine systematically the behaviors of the spin and charge currents, we consider three orthogonal magnetization directions, namely along the x, y, and z axes. In addition, we incorporate pure Rashba (α = 0, β = 0) and Dresselhaus (β = 0, α = 0) SOCs that allow isolation of their effects relative to the physical quantities under study. When finding solutions to the Usadel equation, Eq. (2), and the corresponding current densities [Eqs. (10) and (11)], we have added a small imaginary part, δ ≈ 0.01∆ 0 , to the quasiparticles' energy, ε → ε + iδ, to enhance stability of the numerical solutions. The imaginary part can be physically viewed as accounting for inelastic scatterings. 68 Due to the presence of the finite parameter δ, we take the modulus of the currents in Eqs. (10) and (11). We normalize the quasiparticles' energy, ε, and exchange field h by the gap, ∆ 0 , at T = 0. Also, all lengths are measured in units of the superconducting coherence length ξ S . In our computations, we adopt natural units, so that k B = = 1.
To begin, we consider for comparison purposes, an S/F /S Josephson junction in the absence of SOCs. 28,29 The schematic of the S/F /S structure is depicted in Fig. 1. The parameters ζ = 4, | h| = 10∆ 0 and d F = 2.0ξ S , ensure the validity of low proximity limit considered throughout the paper. To have absolute comparisons, we set h = (0, 0, h z ) and compute the charge and spin currents using Eqs. (10) and (11), respectively. Figure 2 exhibits the spatial map of the spin current for W F = 2.0ξ S (see Fig. 1). Since the magnetization orientation is fixed along the z direction, J sz x (x, y = y 0 ) is the only nonvanishing component of spin current for a given fixed location y 0 . The top panel of Fig. 2 illustrates the spatial variations of J sz x (x, y = y 0 ) along the junction length in the x direction x (x, y), J sy x (x, y), and J sz x (x, y) in a uniformly magnetized Rashba S/F /S junction. The exchange field of the ferromagnetic strip points along the z direction: h = (0, 0, h z ) (see Fig. 1). The panels in the top row show the spin current components, J sγ x (x, y), as a function of x at four differing locations along the junction width: y=0.25ξS, 0.5ξS, 0.75ξS, and 1.0ξS. The bottom row exhibits J sγ x (x, y) versus y at x=0.25ξS, 0.5ξS, 0.75ξS, and 1.0ξS.
at differing positions along the junction width: y 0 = 0.25ξ S , 0.5ξ S , 0.75ξ S , and 1.0ξ S . The macroscopic phase difference between the S electrodes is set at a representative value, i.e., ϕ = π/2. The bottom panel in Fig. 2 shows J sz x (x = x 0 , y) as a function of y, at x 0 = 0.25ξ S , 0.5ξ S , 0.75ξ S , 1.0ξ S . The results demonstrate that the spin current is y independent in such hybrid junctions, namely J sz x (x = x 0 , y) = const. (we also have found J sz y (x, y) = 0). In other words, it is appropriate to view this type of system as an effectively onedimensional junction. The variation of J sz x (x, y) along the x direction is a consequence of spin torque transfer, and hence the spin current is not a conserved quantity. 44,45,47,59 The spin current is maximal at the S/F interfaces and vanishes at the middle of junction, x = 1.0ξ S = d F /2. This is contrast to the charge supercurrent in the F region, which is conserved, and thus has a constant value within the entire F strip (not shown).
To identify some of the salient features in Fig. 2, we consider now a simplified one-dimensional S/F /S system, which permits analytical expressions for the spin current density. To this end, we linearize the Usadel equation, and incorporate the Kupriyanov boundary conditions, where the superconducting electrodes have strong scattering impurities. We also still assume that the magnetization is oriented along z: h = (0, 0, h z ). Correspondingly, we define the dimensionless quantity, λ ± = 2i(ε±h z )/ε T , in which ε T is the Thouless energy, and the dimensionless x coordinate,x = x/d F ∈ [0, 1]. After some straightforward calculations, we obtain the follow-ing expression for the charge current [Eq. (10)]: The charge current in this case is seen to exhibit the usual sin ϕ odd-functionality in the superconducting phase difference. Likewise, by substituting the solutions into Eq. (11), we arrive at the following expressions for the spin current components: Equations (15a)-(15c) clearly demonstrate that the only nonvanishing component of spin current is J sz x , which is consistent with the exchange field aligned along z. 75 From Eq. (15c), it is also evident that J sz x is an odd function of the coordinatex relative to the middle of the junction (and thus vanishes there), and an even function of the phase difference, ϕ. These features are entirely consistent with the numerical results seen in Fig. 2.
We now incorporate Rashba and Dresselhaus SOCs, while keeping the magnetization orientation intact along the z direction. The ISOIs are confined within the F region and are not present within the S electrodes. Through exhaustive numerical investigations, we have found several symmetries among the components of spin current (discussed below) at three particular directions of the exchange field. Due to the symmetries available among the spin current components, we focus here on Rashba SOC. We emphasize that similar conclusions can be drawn for Dresselhaus SOC through the symmetries described below. Figure 3 exhibits the spatial profiles for the spin current density components, J sx x (x, y), J sy x (x, y), and J sz x (x, y). A square ferromagnetic strip is considered, with d F = W F = 2.0ξ S , and the superconducting phase difference is equal to ϕ = π/2. The Rashba SOC coefficient is set to a representative value α = 2.0ξ S , without loss of generality. 73 The top set of panels show J sγ x (x, y = y 0 ) [γ = x, y, z] as a function of the x coordinate at y 0 = 0.25ξ S , 0.5ξ S , 0.75ξ S , and 1.0ξ S . Whereas the bottom panels represent the same quantities, but now as a function of y at x 0 = 0.25ξ S , 0.5ξ S , 0.75ξ S , and 1.0ξ S . As seen in Fig. 1, the junction length and width are parallel to the x and y axes, respectively. The components J sx x (x, y = y 0 ) and J sy x (x, y = y 0 ), shown in the top row of Fig. 3, demonstrate that these spin current densities vanish at the S/F contacts. This finding is consistent with previous works involving nonsuperconducting heterojunctions [43][44][45]47,67 . The z component, J sz x (x, y = y 0 ), however exhibits opposite behavior, and is nonzero at the S/F contacts due to the exchange field, which is oriented along the z axis. Similarly, as seen in Fig. 2, J sz x (x, y = y 0 ) is finite at the S/F interfaces near the S reservoirs. One of the most important features of the results is seen in the top panels of Fig. 3, where two peaks in J sγ x (x, y) emerge near the S/F contacts. We restrict the spatial profiles to 0 < x < d F /2 and 0 < y < W F /2, since the results are symmetrical with respect to x = 1.0ξ S = d F /2 and y = 1.0ξ S = W F /2, so that the maxima of J sγ x (x, y = y 0 ) occurs near the edges of the F wire [at x = 0, and x = d F ]. Turning to the bottom row of panels in Fig. 3, we see that J sγ x (x = x 0 , y) are nonzero at the vacuum boundaries, y = 0, and y = W F . Here also the largest values in the spin current density components take place near the transverse edges of the F wire (y = 0, and y = W F ). The magnitude of the spin current densities at x = 0.25ξ S are generally larger than the other x positions, in agreement with the results of J sγ x (x, y = y 0 ) shown in the top row of panels.
We now consider the effects of changing the magnetization alignment in the ferromagnet. Thus, Fig. 4 represents the same Rashba spin-orbit coupled S/F /S junction as in Fig. 3, except the magnetization of the F wire is now ori-ented along the y axis. This specific direction of h leads to J sx x (x, y = y 0 ) = J sz x (x, y = y 0 ) = 0 at the S/F interfaces and the spin current densities peak near the edges of F wire. The spin current density J sy x (x, y = y 0 ) however is nonzero at the S/F contacts similarly to J sz x (x, y = y 0 ) in the configuration where the magnetization points along the z direction (Fig. 3). As mentioned earlier, this nonvanishing behavior is directly related to the exchange field direction which lies now parallel to the y axis. Examining J sγ x (x = x 0 , y) in the bottom row panels in Fig. 4, the maximal values of J sγ x (x = x 0 , y) take place near the vacuum boundaries, i.e. y = 0, and y = W F . Our investigations demonstrate similar qualitative trends for the components of J sγ x (x = x 0 , y) when the magnetization resides along the x axis. Note that the transverse components of the spin currents are nonzero inside the ISO coupled ferromagnetic wire, i.e., J sγ y (x, y) = 0, and vanish at the vacuum boundaries (y = 0, W F ). We mainly focus here on the J sγ x (x, y), since the longitudinal components contain the relevant information needed to describe and understand the accumulation of spin current densities at the edges of the structures.
Considering now the previous characterization of the spin current components in systems with ISOCs, we schematically summarize the spatial maps in Fig. 5 for J sγ (x, y). The largest amplitudes of J sγ (x, y) reside near the edges of the F strip, i.e. x = 0, d F and y = 0, W F . We have qualitatively marked these regions by light yellow "ribbons". Therefore, the overlap of maximal amplitudes take place near the corner regions of the F strip. We have marked these areas by dashed curves and with a deeper yellow color. The spatial profiles found here are qualitatively similar to the existence of edge spin currents found in nonsuperconducting heterojunctions with ISOIs, 39-46 except with one crucial difference: here spin accumulation at the edges arises in the absence of an external field. As mentioned in other works 67 , the spin accumulation is a signature of the spin Hall effect. Therefore, the predicted spin accumulation in this paper might be measurable through optical experiments 43 , such as through Kerr rotation microscopy 43 , where spatial profiles of the spin polarizations near the edges can be imaged. An alternate experimental proposal involves multiterminal devices 44,45 . When transverse leads are attached to the lateral edges of a two-dimensional S/F /S junction (borders at y = 0, y = W N in Fig. 1(a)), the spin accumulations at the F wire's edges inject spin currents into the leads. 44,45 The transversely injected spin currents into the lateral leads in turn may induce a voltage drop between the additional leads. 44,45 The Josephson effect is a significant example of a macroscopic quantum phenomenon and is of fundamental importance in determining the properties of dissipationless coherent transport. Thus, the behavior of the spin currents upon varying the macroscopic phase difference is crucial to experiments and applications utilizing spin-Hall effects and spin transport. In Fig. 6, we therefore study the spin current components as a function of the macroscopic phase difference, ϕ, between the S banks. We consider the parameter set used in Fig. 4, including h = (0, h y , 0). In the top set of panels, the spatial variations of J sγ x (x, y) are plotted at ϕ = 0, 0.2π, 0.4π, and x (x, y), J sy x (x, y), and J sz x (x, y) in an S/F /S system. The Rashba ferromagnetic wire's width and length are equal to WF = dF = 2.0ξS. The exchange field of the ferromagnetic strip is fixed along the y direction, h = (0, h y , 0). Top row: spatial behavior of J sγ x (x, y) along the junction length, x, at y = 0.25ξS, 0.5ξS, 0.75ξS, 1.0ξS. Bottom row: spatial variations of J sγ x (x, y) along the junction width in the y direction at x = 0.25ξS, 0.5ξS, 0.75ξS, 1.0ξS.
1.0π. We have also chosen a representative position along the junction width, corresponding to y = 1.0ξ S = W F /2, which simplifies the analysis while maintaining the generality of the discussion. Although J sγ x (x, y) has a minimum at y = 1.0ξ S , it exhibits the same trends as a function of ϕ compared to the other positions inside the F wire. By increasing the superconducting phase difference from 0 to π, the amplitudes of the spin current components decrease overall. In the bottom set of panels of Fig. 4, we illustrate J sγ x (x, y) x as a function of ϕ, at y = 0.25ξ S , 0.5ξ S , 0.75ξ S , and 1.0ξ S . Here, we denote the spatial average over the x coordinate from 0 to d F by ... x . In order to better visualize the averaged profiles, we use logarithmic scales in the vertical axes of the bot- The light yellow ribbons display edge regions with maximum spin current densities. The induced spin currents can be considered as a response of the intrinsic spin-orbit coupled system to the presence of an exchange field (the combination of spontaneously broken time-reversal symmetries and the lack of inversion symmetries). As shown in Fig. 1, the exchange field of the ferromagnetic wire is uniform and can take arbitrary directions. The regions that carry maximal accumulation of spin currents are qualitatively shown by the semicircular regions. tom row of panels. As seen, the three components of spin current J sγ x (x, y) x are even-functions of ϕ, with a period of 2π, namely J sγ (2nπ + ϕ) = J sγ (−ϕ), n ∈ Z. This is contrary to the charge supercurrent which is an odd-function of ϕ, i.e., J c (2nπ + ϕ) = −J c (−ϕ) regardless of a finite phase-shift ϕ 0 66,69,70 . These findings are entirely consistent with previous studies of S/F /S Josephson junctions with inhomogeneous magnetization patterns 68 . We here remark that an additional phase-shift ϕ 0 may appear in such junctions due to the coupling of exchange field and ISOIs. 66,69,70 Nonetheless, the explicit current-phase relations simply undergo a shift in ϕ 0 . 66,69 According to the current-phase relations, the charge supercurrent vanishes at certain ϕ that is quite different than the behavior of the spin current components which clearly show nonzero values at the same ϕ. Therefore, these differences in charge and spin currents allows for an examination of edge spin currents without any net charge current in an ISO coupled F wire sandwiched between two S banks.
We are now in a position to discuss symmetries that may arise among the spin current density components for differing magnetization orientations in systems with either Rashba or Dresselhaus SOCs. Our investigations have found that the out-of-plane spin current, J sz (x, y), remains unchanged upon exchanging the Rashba and Dresselhaus SOCs, regardless of the magnetization orientation. This follows from the form of the spin vector potential discussed at the beginning of this section. However, this picture changes for the in-plane J s{x,y} (x, y) components. The x and y components of the spin current become interchanged when transforming from one type of spin-orbit interaction to another. Precisely speaking, by going from Rashba to Dresselhaus spin-orbit coupling, x (x, y) as a function of x along the junction length at ϕ = 0, 0.2π, 0.4π, and 1.0π. The location along the junction width is fixed at the middle of the junction, y = WF /2 = 1.0ξS. The bottom row of panels represents the spatially averaged spin current components over the junction length (denoted by J sγ x (x, y) x) vs ϕ. The average is performed along four positions: y=0.25ξS, 0.5ξS, 0.75ξS, and 1.0ξS. The ferromagnetic wire is a square strip with dF = WF = 2.0ξS, and exchange field h = (0, h y , 0). one simply needs to exchange indices x and y in the components of both the exchange field and the spin current. Otherwise, everything stays the same. By making use of the simple transformation rules described, one can easily deduce the results of Dresselhaus spin-orbit coupled systems from the plots presented for Rashba spin-orbit coupled S/F /S systems shown in Figs. 3, 4, and 6.
To conclude this section, we briefly discuss the importance of having a magnetic element in the Josephson junction for the effect of spin current edge accumulation to take place spontaneously. We thus take the limiting case of h = 0 in our previous calculations above involving S/F /S junctions. Using otherwise the same geometrical and material parameters, this case was found to produce no spin current, J sγ (x, y) = 0, in the presence of Rashba (α = 0, β = 0) and/or Dresselhaus (β = 0, α = 0) SOIs. These findings are consistent with previous works, 67 where several simplifying approximations were employed for Rashba-based S/N /S systems. Examining also the charge supercurrent, J c (x, y), for both the Rashba and Dresselhaus interactions, we observed a uniform spatial map for the charge current density for all ϕ, with J c x = const., and J c y = 0. In other words, the spin-dependent fields cannot induce transverse charge supercurrents in a diffusive S/N /S junction. This is in stark contrast to its ballistic S/N /S counterpart, where a transverse charge supercurrent (that is, equivalent to a supercurrent flowing along the y direction in our configuration depicted in Fig. 1) was theoretically predicted due to the presence of intrinsic SOIs 83 .

IV. CONCLUSIONS
We have theoretically studied the behavior of spin and charge currents in a finite-size two-dimensional S/F /S Josephson junction with intrinsic spin-orbit couplings. We utilized a two-dimensional Keldysh-Usadel quasiclassical approach that incorporates a generic spin-dependent vector potential. Our results demonstrate that the combination of a uniform magnetization and ISOIs drives the spin currents which spontaneously accumulate at the F wire's edges. The corners of the F wire were shown to host the maximum density of spin currents. (As demonstrated in Ref. 75, similar edge phenomena can be found in finite-size two-dimensional intrinsically spin orbit coupled S/N /S junctions with a single spin active interface. Additionally, it was shown that maximum singlettriplet conversions take place at the corners of N wire nearest the spin active interfaces 75 .) Our investigations show that the spontaneous edge accumulation of the spin currents are robust and can exist at all magnetization orientations, independent of the actual type of ISOIs. Our investigations have also found several symmetries among the spin current components upon varying magnetization orientations coupled to a Rashba or Dresselhaus SOI. By varying the superconducting phase difference, ϕ, between the S banks, we determined the spin and charge currents as a function of phase difference. We have found that net spin currents therefore emerge and accumulate spontaneously at the edges, in the absence of charge flow, when properly modulating ϕ in finite-size two-dimensional intrinsically spin-orbit coupled S/F /S hybrid structures. This work can be viewed as complementary to previous studies involving edge spin currents in non-superconducting spinorbit coupled structures where externally imposed fields were required [42][43][44][45]47,50,67 . We have shown that remarkably, edge spin currents can be spontaneously driven by the coupling of intrinsic properties of a system, i.e. spontaneously broken time-reversal and the lack of inversion symmetries in the absence of any externally imposed field.