Antiferromagnetic and ferromagnetic spintronics and the role of in-chain and inter-chain interaction on spin transport in the Heisenberg ferromagnet

Spin-transport and current-induced torques in ferromagnet heterostructures given by a ferromagnetic domain wall are investigated. Furthermore, the continuum spin conductivity is studied in a frustrated spin system given by the Heisenberg model with ferromagnetic in-chain interaction \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J_1<0$$\end{document}J1<0 between nearest neighbors and antiferromagnetic next-nearest-neighbor in-chain interaction \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J_2>0$$\end{document}J2>0 with aim to investigate the effect of the phase diagram of the critical ion single anisotropy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_c$$\end{document}Dc as a function of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J_2$$\end{document}J2 on conductivity. We consider the model with the moderate strength of the frustrating parameter such that in-chain spin-spin correlations that are predominantly ferromagnetic. In addition, we consider two inter-chain couplings \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J_{\perp ,y}$$\end{document}J⊥,y and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J_{\perp ,z}$$\end{document}J⊥,z, corresponding to the two axes perpendicular to chain where ferromagnetic as well as antiferromagnetic interactions are taken into account.

From the end of 80 decade up to now, the spintronics has witnessed a variety of spin related phenomena such as spin transfer torque, tunneling magnetoresistance 1 and so on. The spintronics demands the spin transport study, where the spin current plays a central role in order of spintronics phenomena to occur in magnetic materials and also in various other materials including semiconductors and oxides 2,3 . So, the spin current can be used to control the magnetization via spin-transfer-torque and spin-orbit-torque in several magnetic (and also nonmagnetic) materials [4][5][6][7][8][9] .
On the other hand, the frustrated antiferromagnetism has been also a rich topic nowadays since to very intriguingly phenomena like topological phase transitions that may occur in this class of systems. The frustration leads to an importance of quantum effects because of the classical order is suppressed and novel phases may occur and to govern the physics at low-energy. Moreover, these systems can exhibit a nematic ground state induced by a spontaneous symmetry breaking induced by terms such as frustrating interactions in the Heisenberg model 10 . This nematic ordering can occur when spin fluctuations are taken at some axis without any direction being chosen 11 . In a general way, the frustration is present in materials as LiVCuO 4 , LiVCu 2 O 4 , LiZrCuO 4 that are adequately described by ferromagnetic nearest-neighbor interactions J 1 and antiferromagnetic next-nearestneighbor interactions J 2 [12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29][30] . All of these materials are of spin-1/2. For spin-1, there are a small number of materials as example the material NiCl 2 4SC(NH 2 ) 2 which is a quasi-one-dimensional antiferromagnet with easy-plane anisotropy dominating the exchange interaction [31][32][33] .
The first generation of spintronic devices are based on spin transport, that utilises the magneto-transport being invented in 2001 34 . The injection efficiency depends on the spin polarisation of the ferromagnet and the spin scattering at the ferromagnetic/non-magnetic interface, where it is also important to eliminate any other effects, namely a stray field from a ferromagnet, that distorts the estimation of the injection efficiency.
The aim of this paper is to analyze the effect of scattering among electrons with a ferromagnetic wall domain on spin transport by electrons and to analyse the spin transport by ions on the lattice in a two-dimensional frustrated spin lattice model with the aim to analyse the effect of quantum phase transition (QPT) on longitudinal

Spin polarized current
The spin transport by electrons in a system can be expressed by a spin current in the form where e is the electron charge. In the same way the charge current J c is given by J c = J ↑ + J ↓ . Both currents obey to a diffusion equation given by where µ ↑,↓ is the magnetic momentum of each electron and Ŵ = √ Dτ is the diffusion coefficient, being D the diffusion constant and τ the spin flip time 35,36 . Here, we consider the model described by the Hamiltonian where J denotes the exchange integral and V(r) is the potential (nonmagnetic) of the lattice. The last term in the Equation above V (r) = Js · S(r ′ ) provides the interaction between electron with the ferromagnetic domain wall. We consider a homogeneous magnetic domain wall with a collinear magnetization. The interaction among the electrons spins with the spins of the ferromagnetic wall is represented in Fig. 1, where the potential of interaction among the electron spins with the spins of the wall domain is given by V (r) = Js · S(r) . The the aim here is to verify the influence of this on spin wave function of the electrons. The transmission of electron through a domain wall was discussed in Refs. 37,38 . The purpose here is to use the Born's expansion for f ↑,↓ (k, k ′ ) to calculate the effect of interaction electron spin-domain wall on spin wave function of the electrons of the current J c . The wave function of the electron at large distance from the wall domain is given by where ψ in (r) = ψ(−∞) is the wave function of the electron after the scattering with the wall. Consequently, far from domain wall, we have ψ out (r) given by where |ψ ↑,↓ (k, r)� = S|ψ ↑,↓ (k, r)� and S is the scattering matrix.
where µ B is the Bohr's magneton and T is given by the Lipmann-Schwinger's equation ∇ 2 ψ(r) + V(r)ψ(r) + Js · S(r)ψ(r). where the A(k) coefficient is given by The integral above was solved approximately as where δ is the width of the wall. The potential of interaction among electron spin with the spins of the domain wall V(x) has the form V (x) = J sd S cos(4 arctan(e −δx )) . The shape of the potential is displayed in Fig. 2. We consider the expansion of the Eq. (9) up to first order. An analysis considering terms of superior order will generate a large quantity of terms in the Eq. (14) and should not generate any change in the scattering. We obtain a very complicated expression for the wave function of the electron after the scattering with the ferromagnetic wall domain however, in a combination of two polarization states. The presence of the coefficient f (k, k ′ ) in the second term making the control of the state of polarization of each electron after the scattering with the domain wall a very difficult problem.
The electron Hamiltonian interacting with the ferromagnetic domain wall can be written as 39,40 where H 0 is the Hamiltonian of the free electron, H sw is the electron-domain-wall Hamiltonian and H w is the Hamiltonian of the wall domain.
Making the transformation of the spin operators where H ′ contains terms of four or more operators a i and A i . The contribution of the interaction between each electron with the domain wall on the electric current operator J sw is given by We use the Matsubara's Green function method at finite temperature 39,40 to determine the contribution of the interaction of the electrons with the wall domain for the regular part of the electric conductivity or continuum conductivity, σ reg (ω) , given by where f k = �n i↑ � = �n i↓ � = (e βω k + 1) −1 is the fermion occupation number and N k = (e βW k − 1) −1 is the boson occupation number associated with the spin waves of the wall domain and β = 1/T . We have that in the low energy limit where v is the Fermi's velocity. In the Fig. 3, we present the behavior of the contribution of the interaction among electrons with the domain wall, σ reg sw (ω) . How the electric resistance is the inverse of the electric conductivity, the inverse of σ reg sw (ω) provides the information about the electric resistance generated by the ferromagnetic domain wall. Our results show a peak of resonance in the contribution of the spin electron-ferromagnetic wall domain at range ω ≃ 2.5J , which indicates a peak in the electric conductivity at this range of ω.
cos k x , + cos k y + cos k z www.nature.com/scientificreports/

Meissner effect for the spin supercurrent
It is a fact well known that an example in the nature of local spontaneous breaking of gauge symmetry is the superconductivity. How the charge conductivity is a response to a time-dependent electric field given by Ohm's law, J (x, t) = σ E(x, t) , in a similar way, we have that the spin current flows in response to a magnetic-field gradient following the Fick's laws for ∇B(x, t) as �J (x, t)� = σ ∇B(x, t) . Hence, so as for electric superconductors, we should have here σ → ∞ for the case of a spin superconductor, where in a finite system with N sites we must have a finite number of spins. As in general, the spin conductivity cannot be infinite consequently, the gradient of the external magnetic field ∇B(x, t) must be zero inside of the a spin superconductor as in the electron superconductivity. So, if B is zero in the beginning, it must be zero inside the superconductor even if we apply a gradient of an external magnetic field outside the superconductor. This means that the applied external magnetic field must not depend on x inside the superconductor. Consequently, the spins of a spin superconductor should generate a current that screens the external gradient of the magnetic field. The action of the system becomes invariant under the gauge transformation We introduce the Goldstone boson field φ(x) that has the property φ(x) = φ(x) + �(x) and where q is the charge. Thus, the magnons in the Heisenberg model must be described by a charged scalar field so as the cooper pairs in the electric superconductor where the Lagrangian describing the interaction between this scalar field with the gradient of the electromagnetic field being given by being L m a not well known functional. The Proca's equation is given by with Thus, the scalar field ψ produces a mass for the photons following the Higg's mechanism where the scalar field that describes the spin waves plays the role of a Higg's field developing a antiferromagnetic vacuum expectation value. So, the photon must acquire a mass inside the spin superconductor like in the electric superconductor with the wave equation becoming the Klein-Gordon equation for the quantity B defined by B = ∇B(x, t) given by the massive Klein-Gordon equation www.nature.com/scientificreports/ Thus, if we apply a gradient of an external magnetic field, the solution of equation above will become ∇ 2 B = 0 for x < 0 and (∇ 2 − � 2 )B = 0 for x > 0 . The solution for x > 0 gives a penetration length given by l = 1/� , being the inverse of the mass that the photon acquires inside of the spin superconductor.

Spin transport in quantum frustrated spin-1/2 ferromagnet
We discuss the role of the interchain coupling on spin transport in coupled frustrated spin-1/2 magnets with a ferromagnet NN in-chain coupling J 1 < 0 and an AFM NNN in-chain coupling J 2 > 0 , where the chains are aligned along the x axis, and they are coupled along the y and z axes by J ⊥,y and J ⊥,z , respectively. Thus the model is given by where i, j , x, y, z labels NN bonds along the corresponding axis, and [i, j], x labels NNN bonds along the chain. Furthermore, we consider J 1 < 0 and J 2 ≤ 0 , whereas no sign restrictions are valid for J ⊥,y and J ⊥,z .
In the linear response theory, the spin conductivity is the response to an actual frequency-dependent gradient of magnetic field is given by [41][42][43] where we have introduced the spin-flip part of the exchange interaction along the x direction. The Eq. (35) exhibits the desired structure �J � = σ ∇B z , where the formula for the spin conductivity is defined as the linear spin-current response to a uniform, k = 0 , frequency-dependent gradient of the magnetic field 44 . In general, the k = 0 conductivity at T = 0 may be written as Re σ βγ (ω) = D S δ(ω) + σ reg (ω) , where D S is the Drude's weight. Therefore, beyond the k = 0 the relation between the "twist conductivity" and the response to an inhomogeneous magnetic field is not clear.
The Green's function at zero temperature is defined as 42 where T is the time ordering operator. The current-response function G(k, ω) at finite temperature is given by G(k = 0, ω → 0) is the susceptibility or retarded Green's function 42,43 .
The operator for spin current from site j to site j + x is defined by [44][45][46][47][48] where j + x is the nearest-neighbor site of the site j in the positive x direction. Furthermore, the spin-current operator can be expressed as We find the spin current operator in terms of boson operators α and β given by where the higher-order terms in Eq. (39) involves terms of four or more boson operators and have been discarded. We find the Green's function given by where g kk ′ (ω) is obtained by applying the Wick's theorem following of the Fourier transform after performing a tedious calculation. The Eq. (40) corresponds to the lowest approximation (noninteracting magnetic excitations) which replaces the magnon propagators by the free propagators G 0 and hence, it is valid only in the mean-field approach. Thus, the Green's function for ψ is �α(t)α † (0)� → G 0 (t) and �β † (t)β(0)� →G 0 (t) , where G 0 is the bare propagator. www.nature.com/scientificreports/ and Furthermore, we employ the formula 42 We obtain the retarded Green's function and σ reg (ω) , using the SU(3) Schwinger boson theory and applying the Gree-Kubo formula we obtain the continuum conductivity given by We find the spin conductivity using the SU(3) Schwinger boson theory given as a second-rank tensor, being different from results in literature obtained using the Dyson-Maleev representation which is given by a scalar or zero-order tensor 42,43 . This difference implies in a different response of the spin current to the gradient of the external magnetic field ∇B. We can improve SU(3) Schwinger boson mean-field formalism including the fluctuations around the meanfield result 49 . We can consider the phase fluctuations φ ij around the mean-field results for A as A ij =Ãe iφ ij taking into account in the action for the Hamiltonian Eq. (34) which is invariant under the gauge transformation: In Fig. 4, we obtain σ reg (ω) at T = 0 using the SU(3) Schwinger boson formalism. We perform the calculations for the following values of parameters: J 1 = −1.0 , J y = 0.1 , J z = 0.0 and J 2 = 0.2 , where the system is near to line of phase transition in the graphic D C vs. J 2 : 0.251 ≤ J 2C ≤ 0.252 , where D C goes continuously to zero D C → 0 33 . There is a discontinuous drop of D C at J 2C which is an effect of the dimension of the system. Being the behavior for three-dimensional and two-dimensional model are similar as was found in Ref. 32 . Where below D C the system is ordered and the magnetization is non-zero and the magnetization goes to zero at D C . We obtain a large peak for the spin conductivity at ω ≈ 0.03 and a finite spin conductivity at ω → 0 and none divergence at DC limit, which is due to the behavior of the expression for σ reg (ω) . Furthermore, this peak is due to behavior of dispersion relation ω k in the Brillouin zone, where the main effect of the frustration occurs at the long-wavelength limit where the gap of the dispersion relation, in the disordered phase vanishes at critical parameter D C → 0 , at a wave vector signaling the magnetic order that appear below D C . Experimental results can be compared with our results when available. As far as I know, there is no result for the spin transport for the model studied here. As In Fig. 5, we obtain σ reg (ω) as a function of J 2 and for the values: J 1 = −1.0 and J y = 0.1 . We perform the calculations for the ω near to the peak of the spin conductivity obtained in Fig. 4, ω 0 = 0.03 . Due to symmetry of equation for conductivity, we should obtain the same behavior for the conductivity as a function of J y . In Fig. 6, we obtain σ reg (ω) as a function of J 2 for J 1 = −1.0 and J y = 0.1 . We perform the calculations also near to www.nature.com/scientificreports/ the peak of the spin conductivity, ω ≈ ω 0 = 0.03 . In this case, we obtain a monotonically increasing behavior for the conductivity as expected.

Conclusions
In brief, we analyze the effect of scattering among electrons with a ferromagnetic wall domain on spin transport by electrons and analyze the spin transport by quasi particles in the frustrated Heisenberg ferromagnet with inchain and inter-chain interactions which is very important in the field of antiferromagnetic and ferromagnetic spintronics. Additionaly, the Meissner effect for the spin superconductivity obtained for many other frustrated spin systems is also proposed. We analyze the case of a lattice model described by the isotropic Heisenberg model with a ferromagnetic in-chain interaction J 1 < 0 between nearest neighbors and an antiferromagnetic next-nearest-neighbor in-chain coupling J 2 > 0 . We obtain a large variation of the spin conductivity with the frustration parameters: J 2 , J ⊥,y and J ⊥,z . There is also a Drude weight for the spin conductivity which exists in many other spin models. However, the purpose here is to analyze the effect of phase transition on continuum conductivity ω = 0 where the Drude weight term does not generate any influence. The peak of the spin conductivity can be determined by measuring of magnetization current 44,50 . In Ref. 51 , some experimental techniques were proposed and seem to be feasible. Another experimental technique that can be used is the nuclear magnetic www.nature.com/scientificreports/ relaxation (NMR). The spin transport in the compound AgVP 2 S 6 was experimentally investigated using this technique (NMR) in Ref. 52 , where the experiment was performed at high temperatures where the behavior for the spin transport is diffusive.

Appendix: SU(3) Schwinger bosons
We consider the model given by The SU(3) Schwinger boson theory is a theory proposed in Refs. 53,54 and is well adequate to study the anisotropic XXZ model with single-ion anisotropy D in the range D > D c , for N → ∞ limit where three boson operators t x , t y and t z are defined as where |v� is the vacuum state. In terms of these t γ operators ( γ = x, y, z ), we write the spin-S operators as We introduces other boson operators u † and d † given by

Thus
In addition, we impose the constraint condition u † u + d † d + t † z t z = 1 . We employ the further condition that �t z � = t 54 and introduce the Lagrange multiplier µ j (T) that is a chemical potential of bosons depending on each site j of the lattice and temperature. From the mean-field approximation, we let µ j (T) = µ(T) . The other parameters in theory: t 2 , µ , p 1 , p 2 , p y and p z are obtained numerically by a set of integral equations as following After performing a Fourier transformation followed by the Bogoliubov transformation with we find (45) H = J 1 �i,j�,x S i · S j + J 2 [i,j],x S i · S j + J ⊥,y �i,j�,y S i · S j + J ⊥,z �i,j�,y S i · S j + D i (S i ) 2 .