Magnetism induced by nonlocal spin-entangled electrons in a superconducting spin-valve

In the traditional view, the magnetic moment appearing in the superconducting region is induced by equal-spin triplet superconducting correlations in superconductor ($S$) ferromagnet ($F$) heterostructure with noncollinear magnetization. In this paper, we represent that in $NSF_1F_2$ ($N$--normal-metal) spin-valve structure the induced magnetic moment emerging in both the $S$ and $N$ regions can also be generated by Cooper pair splitting: one electron coherently tunnels from the $S$ layer into the $F_1$ layer, and the other one stays in the $S$ layer or tunnels into the $N$ layer. Two electrons are spatially separated from each other but their total spin ground state is entangled in this process. In contrast, the magnetic moment induced by the equal-spin triplet correlations hardly penetrates from the $S$ layer into the $N$ layer. In particular, by tuning the size of the exchange field and the thickness of the $F_1$ layer, one may control the direction of the induced magnetic moment in the $N$ layer. This interesting phenomenon can be attributed to the phase-shift obtained by the spin-entangled electrons. Our theoretical proposal will offer an effective way to control the entanglement of the nonlocal electrons, and also may provide possible explanations for previous and recent experimental observations [Stamopoulos et al 2005 Phys. Rev. B 72 212514; Ovsyannikov et al 2016 J. Exp. Theor. Phys. 122 738; Flokstra et al 2016 Nat. Phys. 12 57].


Introduction
The interplay between superconductivity and ferromagnetism in hybrid structures has currently attracted considerable attention because of unusual physical phenomena [1,2,3,4] and potential practical applications [5,6,7,8,9]. Various device applications of superconducting hybrid structures have been widely studied in recent years [10,11,12]. Much effort has been devoted to obtaining a better understanding of these phenomena appeared in heterostructures involving superconductor (S) and ferromagnet (F ). It is well known that ferromagnetism and conventional superconductivity are two antagonistic orders, as ferromagnetism favors a parallel spin alignment, while Cooper pairs consist of electrons with antiparallel aligned spins. The interaction of these two orders leads to the proximity effect and the inverse proximity effect.
The proximity effect causes the superconducting correlations to penetrate into the ferromagnetic region. When the superconductor is adjacent to a homogeneous ferromagnet, the Cooper pairs entering the ferromagnetic region acquire a finite centerof-mass momentum Q ≃ 2h 0 / v F due to the exchange splitting of the ferromagnet.
Here h 0 and v F are the exchange field strength and the Fermi velocity, respectively. The wave function of the Cooper pair oscillates and decays in the ferromagnetic region as a function of QR, where R is the distance from the S/F interface [6,7,8,9]. This oscillatory nature leads to periodic π-phase shifts across the SF S junction [13]. Another unusual effect highlighted in the SF heterostructures is the fact that inhomogeneous textures of the magnetization in the ferromagnet may lead to the creation of spin triplet superconducting correlations with equal spin projections on the quantization axis [8,9,14]. These triplet correlations are not destroyed by the exchange field and can propagate a long distance in the strongly spin-polarized ferromagnet and half-metal. The experimental evidence of such triplet correlations is revealed by the recent observations of the long-range Josephson current [15,16,17,18,19] and the superconducting transition temperature in the S/F spin-valve heterostructure [20,21,22,23].
In contrast, the inverse proximity effect induces ferromagnetic correlation in the superconducting region near the S/F interface. It was proposed theoretically that in a SF bilayer the magnetic moment could be induced in the superconducting region and its direction is opposite to the magnetization in the bulk of the ferromagnet [24,25]. This induced magnetic moment also displays an oscillatory sign-changing behavior changing with the product h 0 d of the ferromagnetic thickness d [26,27]. If the magnetization distribution is inhomogeneous in the ferromagnetic region, then the equalspin triplet correlations give rise to an induced magnetic moment in the superconducting layer [28,29,30]. Previously, Stamopoulos et al. [31] offered an experimental result in the S/manganite multilayers/F structures. They observed that the inhomogeneous magnetization structure of the manganite multilayers, which was modulated by the external magnetic field, could efficiently switch the direction of the magnetic moment appearing in the superconducting region. Recently, Ovsyannikov et al. [32] reported experimentally that the magnetic moment emerged in the superconducting part of the heterostructure that consists of a cuprate superconductor and a ferromagnetic spin valve when the magnetization vectors of the ferromagnetic films have noncollinear orientation.
On the other hand, it is known that the conventional superconductor in principle is considered as a natural source of spin entanglement [33,34]. When two spatially separated normal-metal electrodes form two separate junctions with a superconductor, with the junction separation of the order of the superconducting coherence length of the material, the Cooper pair in the superconductor can break up into two nonlocal entangled electrons that enter into different electrodes via so-called the crossed Andreev reflection or Cooper pair splitting (CPS) [35]. This effect can achieve the nonlocal electronic correlations [36,37,38,39,40,41]: two entangled electrons are localized in the different normal-metal electrodes, and while being spatially separated from each other. What happens if one of the normal-metal electrodes is replaced by an inhomogeneous ferromagnet? Can the electronic spin in the normal-metal electrode be dominated by the characteristics of the ferromagnet? Based on these questions, Flokstra et al. [42] represented a surprising experimental observation in NSF 1 F 2 multilayered structure. The induced magnetic moment appears inside the normal-metal, but not in the adjacent superconducting layer. In particular, the magnetic moment exhibits a spin-valve effect: a significant change in magnitude depends on the mutual orientation of magnetization in two ferromagnetic layers. Nevertheless, a key issue still needs to be resolved and further studied: which component, the spin singlet correlation or the equal-spin triplet correlation, makes a significant contribution to the magnetic moment induced in the normal-metal?
In this paper, we demonstrate that the CPS effect could induce a magnetic moment inside the superconductor and the normal-metal in the NSF 1 F 2 structure with noncollinear magnetization. The equal-spin triplet pairs within the superconducting region also produce a magnetic moment in the different direction, but this magnetic moment hardly penetrates from superconductor to normal-metal. In this spin-valve structure, the induced magnetic moment in both the superconductor and the normalmetal is different from the homogenous SF structure, in which the magnetic moment occurs in the superconducting region only when the exchange field of the ferromagnet is weak enough. However, in our structure the induced magnetic moment still exists in the superconductor and the normal-metal even if the F 2 layer converts into a strongly spinpolarized ferromagnet. In particular, the direction of the induced magnetic moment could be reversed by changing the exchange field and the thickness of the F 1 layer. This feature can be attributed to the phase-shift obtained by the entangled electrons. However, the change of the F 2 layer triggers an amplitude oscillation of the induced magnetic moment inside the N region but does not reverse its direction. When the exchange field or the thickness of the F 2 layer becomes large, the oscillation will decrease accordingly. Moreover, the induced magnetic moment depends on the temperature and the misorientation angle between the two ferromagnetic layers. This effect can be used for engineering the cryoelectronic devices manipulating the magnetic moment in the normal-metal.

Model and formula
The NSF 1 F 2 junction we consider is shown schematically in figure 1(a). We denote the layer thicknesses by L N , L S , L 1 , and L 2 , respectively. The origin of the y axis, that is perpendicular to the layer interfaces, locates at the outer N surface. The magnetization of the F 1 layer is tilted by misorientation angle θ from the z-axis in the x-z plane, while the magnetization of the F 2 layer is aligned along the z-axis. We also assume that the whole system satisfies translational invariance in the x-z plane. The BCS mean-field effective Hamiltonian [6,43] is where H e = − 2 ∇ 2 /2m − E F , ψ † α ( r) and ψ α ( r) represent creation and annihilation operators with spin α, and the vector σ = (σ x , σ y , σ z ) is composed of Pauli spin matrices. m is the effective mass of the quasiparticles in the system, and E F is the Fermi energy. ∆( r) describes the usual superconducting pair potential. The exchange field exists only in the ferromagnetic region. It can be described by wherex (ẑ) is the unit vector along the x (z) direction. In order to diagonalize the effective Hamiltonian, we make use of the Bogoliubov transformation ψ α ( r)= n [u nα (y)γ n +v * nα (y)γ † n ] and take into account the anticommutation relations of the quasiparticle annihilation operatorγ n and creation operatorγ † n . The resulting Bogolubov-de-Gennes (BdG) equation can be expressed as [43] whereĤ(y) = H e1 − h z (y)σ z − h x (y)σ x and1 is the unity matrix. Besides,û n (y) = [u n↑ (y), u n↓ (y)] T andv n (y) = [v n↑ (y), v n↓ (y)] T are quasiparticle and quasihole wave functions, respectively.
To acquire the characteristics of the physical quantities in the system, we solve the BdG equation (2) by Bogoliubov's self-consistent field method [43,44,45,46]. The NSF 1 F 2 junction is placed in a one-dimensional square potential well with infinitely high walls, then the corresponding quasiparticle amplitudes can be expanded in terms of a set of basis vectors of the stationary states [47], u nα (y) = p u α np ζ p (y) and v nα (y) = p v α np ζ p (y) with ζ p (y) = 2/L sin(pπy/L). Here p is a positive integer and L = L N + L S + L 1 + L 2 . The pair potential in the BdG equation (2) satisfies the self-consistency condition [43] where k B and T are the Boltzmann constant and the temperature, respectively. The superconducting oupling parameter g(y) will be taken to be constant in the superconducting region and zero elsewhere. The primed sum is over the eigenstates within range |E n | ≤ ω D , where ω D is the Debye cutoff energy. The matrix elements in equation (2) are then written as where ε ⊥ in equation (4) is the continuous energy in the transverse direction. The BdG equation (2) is solved by an iterative schedule. One first starts from the stepwise approximation for the pair potential and iterations are performed until the change in value obtained for ∆(y) does not exceed a small threshold value. The local magnetic moment of the NSF 1 F 2 geometry has two components [45] where µ B and f n are the effective Bohr magneton and the Fermi function, respectively. It is convenient to normalize these two components to −µ B . Because the local magnetic moment in the normal-metal is rather small, we define an effective magnetic moment dy to describe the total effect of the magnetic moment induced in the normal-metal region.
The amplitude functions of the spin triplet pair with zero spin projection and the equal-spin triplet pair are defined as follows [45] where the sum of E n is in general performed over all positive energies, and η n (t) = cos(E n t) − i sin(E n t) tanh(E n /2k B T ). Additionally, the amplitude function of the spin singlet pair can be written as In this paper the singlet and triplet amplitude functions are all normalized to the value of the singlet pairing amplitude in a bulk superconducting material. The LDOS is obtained according to [45] where f ′ (ǫ) = ∂f /∂ǫ is the derivative of the Fermi function. The amplitude of the LDOS is normalized by its value at ǫ = 2∆ beyond which the LDOS is almost constant.

Results and Discussions
In this section we present and discuss the results obtained by solving numerically the BdG equation (2). All lengths and the exchange field strengths are measured in units of the inverse Fermi wave vector k F and the Fermi energy E F , respectively. The superconducting coherence length and the Debye cutoff energy are taken as k F ξ S =100 and ω D /E F =0.1, respectively. The N layer thickness is assumed to be k F L N =70 and all interfaces in the system are fully transparent.
In figure 2 we show the detailed dependence of the effective magnetic moments M N x and M N z inside the N region on the exchange field h 1 and the thickness L 1 . It is clear see that M N x does not exist in the parallel configuration (θ=0), but it displays an oscillating sign-changing behavior upon increasing h 1 or L 1 in the perpendicular configuration (θ=0.5π). This indicates that the direction of M N x can be reversed by varying the strength of h 1 or L 1 . By contrast, M N z changes irregularly and its magnitude is usually very small in the above two cases. Since the exchange field h 1 is hardly variable in the experiment, one may hope to observe the oscillations of M N x performing measurements on samples with the different thickness L 1 . It needs to emphasize that in figure 2 (and figure 8 depicted below) the intervals between data points have been taken to be large enough for the sake of beauty. In fact, there are other small oscillation peaks in the interval, which may be ascribed to other mechanisms. For simplicity, we ignore these additional oscillation peaks.
In general, since the entangled electrons penetrate into the ferromagnetic region, their spins will be arranged along magnetization direction of the ferromagnet. When  both ferromagnetic layers are all parallel to z-axis, m z will be induced in both the S and N regions for the weak exchange field. This behavior will be discussed in detail below. If the total exchange field of two ferromagnetic layers is chosen to be strong enough, which causes a depairing effect, it is difficult for the entangled electrons to tunnel into these ferromagnetic regions, then the CPS effect can hardly occur. Correspondingly, m z shows a quite small amplitude. When the magnetization directions of both ferromagnets are perpendicular to each other (θ=π/2), the physical picture will become more complicated. Two effects will cause two magnetization components along different directions. Firstly, the equal-spin triplet pair generates a component m z in the S region but not in the N region. Secondly, the CPS effect produces the other component m x in both the S and N regions. The direction of m x will be reversed with varying characters of F 1 layer because the entangled electrons obtain an additional phase in this situation. In the following, we will demonstrate the above descriptions in detail using the numerical calculation results.
In figure 3 (a) we show the spatial profile of the induced magnetic moment m x for two special exchange fields h 1 /E F =0.035 and 0.1, which corresponds to the points A and B in figure 2 (a), respectively. In the former case (h 1 /E F =0.035), the negative m x decays gradually as a function of distance from the S/F 1 interface. When m x penetrates from the S region into the N region, its sign and magnitude almost remain unchanged  inside the whole N region. This phenomenon can be ascribed to the CPS effect [35]. As illustrated in figure 1 (b), two entangled electrons forming a Cooper pair transport from the S region into the N region and the F 1 region, or one electron stays in the S region and the other one tunnels to the F 1 region. The accumulation of these electrons leads to the magnetic moment in the N and S regions, whose direction is opposite to the exchange field in the F 1 layer. In the latter case (h 1 /E F =0.035), the direction of m x may turn from negative to positive, which is due to the π phase-shift acquired by the entangled electrons. The detailed analysis is described below: the spin singlet state |↑↓ −|↓↑ originating from the S layer can be converted into |↑↓ x e iQ·R −|↓↑ x e −iQ·R in the F 1 layer due to the exchange splitting of the energy bands. This state can be rewritten as a mixture of the spin singlet state and the spin triplet state with zero spin projection on the direction of magnetization: (|↑↓ −|↓↑ ) x cos(QR)+i(|↑↓ +|↓↑ ) x sin(QR). This opposite-spin triplet state (|↑↓ +|↓↑ ) x is equivalent to the equal-spin triplet state −(|↑↑ −|↓↓ ) z when viewed with respect to the z-axis [8,9]. During this process, the magnetization of the F 1 layer is preferable to make the spin of one entangled electron align the same direction, then the spin of the other electron will be antiparallel to the F 1 direction. The entangled electron inside the F 1 (N or S) layer can be described as |↑ x e ik ↑ R (|↓ x e −ik ↓ R ), where k ↑(↓) =k F +(−)Q/2. The contribution of two entangled electrons to the pair amplitude is proportional to e ±i(k ↑ −k ↓ )R =e ±iQR . If the exchange field and the thickness of the F 1 layer are weak and short enough, the phase acquired by the spin singlet pair and the opposite-spin triplet pair satisfies the condition QR < π. The magnetic moment induced in the N and S regions is antiparallel to the magnetization direction of the F 1 layer. By contrast, if the thickness L 1 or the exchange field h 1 becomes large, the above phase is possible to meet the condition QR > π. In this case, the entangled electrons residing on both sides of the S/F 1 interface simultaneously get a π phase-shift, and then the m x direction is reversed. As shown in figures 3 (b) and 3 (d), the pattern configuration of f 0 and f 3 will be reversed as h 1 /E F increases from 0.035 to 0.1. This changing behavior can demonstrate the results mentioned above. Additionally, as described in reference [48] the equalspin triplet pair can be modulated through varying the strength of the exchange field and the thickness of the F 1 layer. It can be seen from figure 4 that the change of f 1 corresponds to the reversal of the m z direction. Two m z components decay in the S region and their amplitudes become quite small in the N region. This illustrates the fact that m z almost cannot be induced in the N region. Taken as a whole, if one changes the strength of the exchange field (or the thickness) of the F 1 layer, the direction of total magnetic moment induced in the S region will be reversed. This theoretical result can be used to explain the experimental findings in reference [31]. In S/manganite multilayers/F hybrid structure, the manganite multilayers including 15 bilayers [La 0.33 Ca 0.67 MnO 3 /La 0.60 Ca 0.40 MnO 3 ] 15 could offer noncollinear magnetization. Because the thickness of the manganite multilayers is about 120 nm and the total effective exchange field is strong enough, its original properties satisfy the condition QR > π, in which case the magnetic moment induced in the S region has a fixed direction. Since the applied external magnetic field exceeds the coercive field, the noncollinear magnetization will be reduced and it matches the condition QR < π, then the induced magnetic moment will change its direction. In the experiment, the equal-spin triplet state f 1 appearing in the system may be manifested through a zero energy conductance peak (ZECP) in the differential conductance spectrum, in principle, which is proportional to the LDOS. In order to further demonstrate the distribution of the equal-spin triplet pair, we show in figure 5 the LDOS in different positions for two special exchange fields h 1 . In the F 1 layer  (k F y=130) there is a sharp ZECP in the LDOS. This ZECP decreases as entering the S layer (k F y=95). Meanwhile, two subgap peaks appear on both sides of the ZECP, which are caused by the inverse proximity effect of the N layer. In the N layer (k F y=50), the subgap peaks are enhanced but the ZECP does not appear. This means that the equalspin triplet pair hardly penetrates into the N layer. Additionally, when the exchange field h 1 takes two different values, the variation characteristics of the LDOS are almost identical except that the peaks at h 1 /E F =0.035 are higher than that at h 1 /E F =0.1.
Next, let us consider the effect of the superconducting thickness on the induced magnetic moment. From figures 6 (a) and 6 (b), we can see that the induced magnetic moment inside the N region decreases with an increase of the thickness L S . It is known that the CPS effect only occurs when the superconducting thickness is less than the superconducting coherence length ξ S . So the induced magnetic moment disappears under the condition that k F L S is larger than k F ξ S =100. By contrast, with the increase of L S the spin singlet state f 3 in the N region will increase (see figures 6 (c) and 6(d)), which indicates that the increase of the superconducting proximity effect enhances the subgap in the N region. This feature can be demonstrated by the enhancement of the subgap peaks in the LDOS, as illustrated in figure 7 (a). The above result shows that the magnetic moment cannot be induced in the N region since two entangled electrons ǫ/∆ forming the Cooper pair tunnel simultaneously into the N region. Moreover, we show in figure 7 (b) the LDOS at different positions when the superconducting thickness is more than the superconducting coherence length. In the S region adjacent to the F 1 layer (k F y=160), there is a distinct ZECP due to the tunneling of the equal-spin triplet pair. Near the center of the S region (k F y=120) the ZECP vanishes but the subgap peaks appear. As the position shifts from the S region into the N region, in which case k F y changes from 120 to 50, the subgap peaks continue to increase until reaching the maximum. In the actual measuring device, similar to the experiment in [42], we propose that Nb and Cu can be selected as possible candidates for superconducting and normal metallic materials, respectively. We now turn to discuss the dependence of the effective magnetic moments M N x and M N z inside the N region on the characteristics of the F 2 layer. As shown in figure 8, for the parallel configuration (θ=0) M N x does not appear and M N z oscillates around the zero value, whose amplitude decays with respect to the external field h 2 (or thickness L 2 ). This shows that under the weak exchange field h 2 one entangled electron tunnels into the F 2 region, and the other electron penetrates into the N region. Since the phase of the entangled electron in the F 2 region shifts with h 2 (or L 2 ), the phase of the other entangled electron inside the N region changes simultaneously. The phase-shift acquired by the entangled electrons can cause the M N z direction to reverse. If the value of h 2 (or L 2 ) becomes large, which leads to an increase of the total magnetization of the F 2 layer, the tunneling of the entangled electron into the F 2 layer will be suppressed, and then M N z decreases accordingly. For the perpendicular configuration(θ=0.5π), M N z presents the same oscillatory behavior, but M N x shows a negative value with superimposed oscillations. These oscillations tend to decrease with increasing h 2 . The reason may be explained by the following picture. For weak h 2 the CPS effect can be divided into two parts: one part occurs between the N and F 1 regions and the other part occurs between the N and F 2 regions. It should be noted that the CPS effect between the N and F 2 regions can induce two components M N z and M N x inside the N region. The phase of the entangled electron within the F 2 layer will shift with increasing h 2 (or L 2 ), which results in the oscillation of M N z and M N x . As the strength of h 2 (or L 2 ) grows to a higher value, the entangled electron tunneling is suppressed, and then the contribution of the entangled electron to M N x and M N z decreases. Therefore, M N z decreases and M N x tends to a constant value. In the actual measurement and application device we recommend that the F 2 layer can be chosen to be a strongly spin-polarized ferromagnet, e.g. Co, whose exchange field is 309 meV [49], and the F 1 layer is taken to be weak one, e.g. Pd 90 Ni 10 , whose exchange field is only 35 meV [49] and the polarized direction can be modulated easily by the applied magnetic field. In such a case, only the CPS effect between the N and F 1 regions still exists, which is beneficial for experimental observation.
The entangled electrons also induce two magnetic components m z and m x in the F 2 region when they tunnel into this region. The m z component is superimposed on the original magnetic moment produced by the exchange field h 2 , and the m x component is influenced by the exchange field h 2 . We show in figure 9 the spatial profiles of m x , f 3 and the imaginary part of f 0 and f 1 in two cases h 2 /E F =0.105 and 0.205, which correspond to two particular points C and D in figure 8 (a). We can see that the pattern configurations of f 3 and f 0 in the F 2 region will be reversed when h 2 /E F increases from 0.105 to 0.205. This behavior gives rise to the change of m x in the F 2 region. Meanwhile, we find that f 1 changes slightly in the ferromagnetic region, which can be manifested through the ZECP in the LDOS. As mentioned previously f 1 may induce the magnetic moment in the S region, but it can not generate an effect in the N region. The dependence of the effective magnetic moment M N x inside the N region on the misorientation angle θ and the temperature T is represented in figures 10. It is possible to find that M N x turns from zero to a negative value as the two exchange fields switch from parallel to orthogonal polarization, which means that M N x is directly related to x-component of magnetization in the F 1 layer. This effect may be used for engineering cryoelectronic devices to manipulate the induced magnetic moment in the N region. Moreover, M N x decreases rapidly with increasing temperature. Particularly, if the temperature rises high enough, M N x completely disappears. This demonstrates that M N x arises from the entangled electron tunnelling into the N region but is not caused by the magnetization leakage from the F 1 region into the N region.
We compare our theoretical result with Flokstra et al 's experiment [42]. Our result shows that the magnetic moment could be induced in both the S and N regions. However, in the experiment [42] the magnetic moment does not appear in the S region but only occurs in the adjacent N region. Here we propose two kinds of conjectures to resolve this contradiction. The first possibility is that the magnetic moment induced in the S region is suppressed by other mechanisms, such as the spinorbit interaction and orbital effects (Meissner currents) [50,51]. From figure 2b in reference [42], one can see that a portion of the magnetic moment still presents in the S region although its amplitude is relatively small. We speculate that the accumulation of the entangled electrons will induce a faint Zeeman field in the superconducting region. The flux caused by the Zeeman field will be expelled by the Meissner effect. So the induced magnetic moment can not be measured in the experiment [42]. Very recently, an anomalous Meissner screening has been observed in normal-metal/superconductor (Cu/Nb) and normal-metal/superconductor/ferromagnet (Cu/Nb/Co) thin films [52]. Comparing with isolated superconductor (50-nm-thick Nb film), the flux expulsion becomes significantly enhanced when adding an adjacent normal-metal (40nm Cu layer). This indicates that the added normal-metal is effective to help the superconductor expel flux. Moreover, a further significant enhancement of the flux expulsion is observed when adding a ferromagnet (2.4nm Co layer) to the other side of the superconductor (Nb). From this experiment, we can infer that the magnetic flux expulsion in the NSF is stronger than that in the SF . So the experimental results in reference [52] may account for the fact that the magnetic moment induced in the S region can be observed in the SF 1 F 2 [32] but not in the NSF 1 F 2 [42]. The second possibility is that the direction of the magnetic moment induced in the S region deviates from that in the N region, hence the magnetic moment inside the S region is not easily measured in the experiment. As mentioned before, the magnetic moment along the x-axis is simultaneously induced in both the N and S regions by the CPS effect, but the magnetic moment along the z-axis is only produced in the S region by the equal-spin triplet pairs. Therefore, the direction of the total magnetic moment in the S region deviates from the x-axis. On the other hand, there are two same effects in our result and Flokstra et al 's experiment [42]: (i) The induced magnetic moment in the N region exhibits a spin-valve effect, which means that the magnetic moment depends on the mutual orientation of the two ferromagnetic layers. (ii) The induced magnetic moment can be controlled by temperature. Below superconducting critical temperature T c , the magnetic moment appears in the N region for the orthogonal arrangement, but above T c the corresponding magnetic moment disappears.
Finally, we give a brief discussion on the disorder effect. In the present work, the clean limit is taken by considering that the characteristic size of the total NSF 1 F 2 structure is smaller than the electron mean free path. The good agreement between our recent work [48] and experiments [15,16] indicate that the clean limit used here is a good approximation, and the disorder due to impurity scattering and interface roughness is of no importance to the changes in the physical characteristics. In our calculation, the non-local electron tunneling plays an important role in the CPS effect. With the increase in barrier strength at the N/S, S/F 1 , and F 1 /F 2 interfaces, the transportation of the entangled electrons will be suppressed. For the same reason, the disorder also gives rise to a reduction in the CPS effect, which is unfavorable to the appearance of the induced magnetic moment.

Conclusion
In this article, we have investigated the magnetic moment induced in the N and S regions in the NSF 1 F 2 junctions with the misorientation magnetization. We find that the induced magnetic moment arises from the spin-entangled electrons accumulating in these regions. In contrast, the equal-spin triplet pair within the S region also can induce a magnetic moment in the different direction, but this magnetic moment hardly penetrates from the S region into the N region. Moreover, the variation of the exchange field and the thickness of the F 1 layer leads to the reversal of the induced magnetic moment inside the N region. This behavior is attributed to the phase-shift acquired by the spin-entangled electrons tunneling into the F 1 layer. On the other hand, the change of the F 2 layer causes the oscillation of the induced magnetic moment, and this oscillating character will be suppressed as the increase of the exchange field in the F 2 layer. In addition, the induced magnetic moment in the N region increases with misorientation angle and reaches a maximum for the orthogonal arrangement, but it decreases with increasing temperature. Thus, all these findings may provide new insight into the physical mechanism for the induced magnetic moment inside the N and S regions.