## Scattering theory of transport through disordered magnets

Martin F. Jakobsen, Alireza Qaiumzadeh, Arne Brataas Center for Quantum Spintronics, Department of Physics, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway

We present a scattering theory of transport through noncollinear disordered magnetic insulators. For concreteness, we study and compare the random field model (RFM) and the random anisotropy model (RAM). The RFM and RAM are used to model random spin disorder systems and amorphous materials, respectively. We utilize the Landauer-B¨ttiker formalism to compute the transmission u probability and spin conductance of one-dimensional disordered spin chains. The RFM and the RAM both exhibit Anderson localization, which means that the transmission probability and spin conductance decay exponentially with the system length. We define two localization lengths based on the transmission probability and the spin conductance, respectively. Next, we numerically determine the relationship between the localization lengths and the strength of the disorder. In the limit of weak disorder, we find that the localization lengths obey power laws and determine the critical exponents. Our results are expressed via the universal exchange length and are therefore expected to be general.

## I. INTRODUCTION

conductive to diffusive [29, 43], which has important consequences for magnonics applications in low dimensions.

In magnonics [1-8], the primary focus has recently been on the propagation of spin waves through various types of magnetic insulators. A particular emphasis has been on ordered systems, such as (anti)ferromagnets, and ferrimagnets. An advantage is that the spin current may suffer less Joule heating compared to electric currents, making insulator-magnonics applications potentially much more energy efficient [3, 9]. Numerous successful experiments have generated and manipulated spin currents using the spin Hall effect and the inverse spin Hall effect [10]. A common experimental setup consists of sandwiching a magnetic insulator between two conductors and using the spin Hall effect to generate a spin current in the left conductor that propagates through the magnetic insulator and into the right conductor. The spin current in the right conductor is converted into a charge current via the inverse spin-Hall effect. This provides a useful method to measure the spin current and infer the spin-transport properties of the magnetic insulator [11-25].

A class of materials that has recently attracted attention in the spintronics community is disordered magnetic insulators [26-31]. Notably, a recent experiment claimed that a spin current flowing through a sample of amorphous yttrium-iron-garnet could travel tens of micrometers [32]. This distance is comparable to the spin current propagation length in a crystalline (anti)ferromagnet [25, 33]. More generally, it is crucial to study disordered magnetic materials because almost all materials contain some degree of disorder, which will affect the functional properties of magnonic devices. When the disorder is sufficiently strong, the eigenstates become trapped in a finite spatial region, completely suppressing the transport properties. This phenomenon is known as Anderson localization, and the first discussion of this phenomenon in magnetic systems began in the 1960s [34-42]. Furthermore, it has been shown that even with a small onset of disorder, the transport properties change from

The common sources of quenched disorder in magnetic insulators are i ) randomness due to anisotropies, local fields, and amorphous structure and ii ) frustration due to competing long-range exchange interactions. In this paper we focus on magnetic insulators with quenched disorder due to i ). Two models with these properties are the random field model (RFM) and the random anisotropy model (RAM), where the disorder is caused by the competition between the exchange interaction and the coupling to local random fields and anisotropies respectively. The RFM and RAM is used to model quenched spin disorder and amorphous magnets respectively [44-49]. Experimental realizations of such systems are plentiful [5054].

Furthermore, there are two types of RFM/RAM spin models. The first is the Ising model, where the spins are scalars S i = ± 1 and are randomly pointing either parallel or antiparallel to each other in the ground state [55-62]. The second type is the Heisenberg model, where the spins are vectors S i that in the ground state are pointing noncollinearly in random directions [63-68].

Because the ground state in the RFM/RAM Ising model is relatively simple, it can often be studied efficiently with analytical methods. For example, one can either solve the equations of motion by a transfer matrix approach paralleling Anderson's celebrated work on disordered fermionic systems [69, 70] or one can use fieldtheory methods, particularly the replica trick, replica symmetry breaking, and mean-field theory [54, 71, 72]. Although the RFM/RAM Ising models are analytically accessible, they are only simplified idealizations of a real disordered magnet where the spins are noncollinear. In this work, we wish to focus on systems with noncollinear spins that are harder to describe analytically but exhibit more realistic spin-wave dynamics.

Disordered magnetic insulators with a noncollinear ground state are a notoriously difficult system to describe. Due to their complexity, it is often useful to study

the classical spin waves of the system. Our work is related to a recent study [29, 73] in which the micromagnetic Landau-Lifshitz-Gilbert (LLG) equation was solved using a quasimonochromatic Gaussian wave packet as the initial condition. They found that the width of the wave packet increases in time until it saturates around the localization length of the system; a hallmark of Anderson localization. In systems that exhibit Anderson localization, the localization length decreases as the system becomes more disordered. However, the exact relationship between the localization length and the strength of disorder is far from being well established in noncollinear disordered magnetic insulators. In this work, we attempt to shed some light on these issues.

The localization effect in spin models depends on the dimensionality of the system, similar to disordered fermionic systems [74]. For fermionic systems in one dimension, there is Anderson localization; in two dimensions, the effect remains present but much weaker, while in three dimensions, there is the possibility of both a localized and a delocalized phase. The same observations have been established for disordered magnets [50, 75]. We focus on one-dimensional spin chains. With more computational time, this method can also be applied to two- and three-dimensional systems.

The numerical method that we develop is based on the Landauer-B¨ttiker formalism u [76, 77], which has proven to be extremely useful in studying the transport properties of electronic systems. To the best of our knowledge, such a method has not previously been applied to the RFM/RAM Heisenberg model. In this paper, we investigate the effect of Anderson localization on the spinwave transport properties of a disordered magnetic insulator. To this end, we first determine the relationship between the system size and the transmission probability for different strengths of disorder and then calculate the spin conductance. With this knowledge, we can investigate how the localization length of the system scales with the strength of the disorder. In particular, we calculate and compare the critical exponents of the RFM and the RAM. These quantities provide us with direct insights into how the transport properties of the spin waves are affected by the localization phenomenon that is present due to quenched disorder.

We hope that this theoretical investigation may inspire an experimental investigation into the transport properties of disordered magnetic nanowires [78-82]. In particular, it would be interesting to compare the experimental relationship between the localization length and strength of disorder to the critical exponents that we determine in this work.

The paper is organized as follows. In Sec. II, we introduce the RFM and the RAM Hamiltonians and discuss its ground state. In section III we find the linearized equations of motion, and derive expressions for the spin current and the spin conductance. Section IV. contains our numerical calculations of the scattering properties of the system. In Sec. IV we summarize our results.

## II. THEORETICAL MODEL

In this section we carefully introduce the model we are interested in studying. We start by presenting the Hamiltonians for the RFM and the RAM, and introduce the geometry. We conclude this section by presenting a method to calculate the classical metastable states.

## A. Hamiltonian

To investigate the transport properties of onedimensional disordered noncollinear spin chains, we use the Hamiltonian

$$H ^ { \kappa } = - J \sum _ { i } \mathbf S _ { i } \cdot \mathbf S _ { i + 1 } - K \sum _ { i } \left ( \mathbf S _ { i } \cdot \mathbf n _ { i } \right ) ^ { \kappa + 1 }, \quad ( 1 )$$

where κ = 0 and κ = 1 represent the RFM and the RAM, respectively.

The dimensionless spins S i are attached to a onedimensional lattice with lattice spacing d . The exchange interaction with J &gt; 0 attempts to align the spins. The terms proportional to K encapsulate the quenched disorder of the system, and we choose K &gt; 0 without loss of generality. Each spin S i is coupled to a local random vector n i . The competition between the exchange and the random interactions in Eq. (1) results in a noncollinear disordered ground state. We use the parameter K/J to characterize the strength of disorder.

## B. Geometry

We consider a one-dimensional chain with N lattice sites. The chain is split into three regions that we call i ) the left lead, ii ) the random region, and iii ) the right lead; see Fig. 1.

FIG. 1: A disordered magnet (blue) is sandwiched between two ferromagnets (red). In regions i ) and iii ), n i = ˆ, while in region z ii ) n i , is uniformly distributed on the unit sphere. Consequently, the spins in region ii ) point in random directions, while the spins deep inside regions i ) and iii ) point in the z -direction. The spins close to the two interfaces rotate similar to the spins in a domain wall. The length of the domain-wall-like region is illustrated and given by the exchange length l ex = √ J/Kd .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00817v3-figure-95ee6720-ed63-4805-b405-cdb5b9911719.jpg)

In regions i ) and iii ), we let the number of spins be equal

to N L and N R , respectively. In addition, we let n i point in the ˆ-direction. z In region ii ), we let the number of spins be equal to N rand and n i to point in some random direction uniformly distributed on the unit sphere. Note that far away from the random region (deep inside of the leads), the spins point in the ˆ-direction, z while in the random region, the spins are oriented randomly. In the regions close to the interface, the spins are rotating in a domain-wall-like fashion. The length of this domain-wall region is given by the exchange length l ex = √ J/Kd .

The scattering problem that we are interested in studying can now be realized by exciting coherent spin waves in the left lead propagating towards the random region. As the spin wave approaches the random region, it will be scattered either back into the left lead (reflection) or into the right lead (transmission). We assume semi-infinite leads such that N L and N R - ∞ → .

## C. The ground state

Determining the true ground state of a disordered magnet (collinear or noncollinear) is a very challenging problem. The primary reason is that the randomness results in free energy with many nearly degenerate minima, separated by high energy barriers. The problem of determining the exact ground state of disordered systems is its own research field, and we do not wish to address that problem here [83-94]. However, due to the high energy barriers, the probability of tunneling between different metastable states is small. Hence, in an experiment, the disordered magnet becomes trapped in a state that may differ from the exact ground state when the system is cooled down, depending on the history. Thus, in this paper, we study the transport properties of disordered magnets around classical metastable states.

We can find a classical metastable state of the system by treating the spins as classical vectors obeying the LLG equation of the form

$$\frac { d \mathbf S _ { i } } { d t } = - \gamma \mathbf S _ { i } \times \mathcal { H } _ { i } ^ { \kappa } - \lambda \mathbf S _ { i } \times ( \mathbf S _ { i } \times \mathcal { H } _ { i } ^ { \kappa } ) \,. \quad ( 2 )$$

Here, the first term with γ &gt; 0 describes the spin S i precessing around its instantaneous effective field H κ i = -δH /δ κ S i , while the second term describes the damping towards the direction of the instantaneous effective field. The metastable state is then obtained by specifying some arbitrary initial configuration and allowing the spins to evolve according to this equation for sufficiently long times t - ∞ → .

## III. SCATTERING THEORY

In this section we outline the theoretical approach that we will use to determine the transport properties of the RFM and the RAM. We start by determining the linearized equations of motion, and formulate the scatter- ing problem. Finally, we derive the expressions for the spin current and spin conductance in the linear response regime.

## A. Hamiltonian in terms of spin-wave operators

To study the transport properties of the system, we can perform a Holstein-Primakoff expansion around one of the metastable states. Let us at each site i define a local coordinate system { ˆ ( ) e x i , e ˆ ( ) y i , e ˆ ( ) z i } such that ˆ ( ) is parallel to the spin at site e z i i in the ground state. The spin operator in a low-lying excited state can then be written as

$$\mathbf S _ { i } = \hat { e } _ { z } ( i ) S _ { i } ^ { z } + \hat { e } _ { x } ( i ) S _ { i } ^ { x } + \hat { e } _ { y } ( i ) S _ { i } ^ { y }. \quad \ \ ( 3 )$$

We perform a Holstein-Primakoff transformation of the form

$$S _ { i } ^ { x } \approx \sqrt { \frac { S } { 2 } } \underset { - } { \left ( a _ { i } + a _ { i } ^ { \dagger } \right ) } \,,$$

$$S _ { i } ^ { y } \approx - i \sqrt { \frac { S } { 2 } } \left ( a _ { i } - a _ { i } ^ { \dagger } \right ),$$

$$S _ { i } ^ { z } = S - a _ { i } ^ { \dagger } a _ { i }.$$

In Eq. (4), we have only included the lowest-order terms because we are not interested in studying the interactions between the spin waves. If we substitute Eqs. (3) and (4) into Eq. (1) and introduce the notations ˆ e ± ( ) i = ˆ ( ) e x i ± ie ˆ ( ) and y i n ± i = n x i ± in y i , we obtain a Hamiltonian of the form

$$H ^ { \kappa } = \sum _ { i j } A _ { i j } ^ { \kappa } a _ { i } ^ { \dagger } a _ { j } + B _ { i j } ^ { \kappa } a _ { i } a _ { j } + h. c, \quad \ \ ( 5 )$$

where

$$& \text{where} \\ \text{LG} \\ & \quad A _ { i j } ^ { \kappa } = \delta _ { i, j } \left \{ J S \hat { e } _ { z } ( i ) \hat { e } _ { z } ( i + 1 ) + \frac { 1 } { 2 } K n _ { i } ^ { \kappa } \\ & \quad + \kappa [ K S ( n _ { i } ^ { \kappa } ) ^ { 2 } - \frac { 1 } { 2 } K S n _ { i } ^ { - } n _ { i } ^ { + } - \frac { 1 } { 2 } K n _ { i } ^ { \kappa } ] \right \} \\ ( 2 ) & \quad - \frac { J S } { 2 } \delta _ { i, j + 1 } \hat { e } _ { - } ( i ) \hat { e } _ { + } ( j ), \\ \kappa _ { i } = \quad B _ { i j } ^ { \kappa } = - \kappa \frac { K S } { 2 } ( n _ { i } ^ { - } ) ^ { 2 } \delta _ { i, j } - \frac { J S } { 2 } \hat { e } _ { - } ( i ) \hat { e } _ { - } ( i + 1 ) \delta _ { i, j + 1 }. \\ \text{type} \\ & \quad \text{In the following, we will study the spin waves associated}$$

In the following, we will study the spin waves associated with the Hamiltonian of Eq. (5).

## B. Equations of motion

The equations of motion for the spin-wave operators can now be calculated from the Heisenberg equation

$$\frac { d } { d t } a _ { i } ^ { \pm } = \frac { i } { \hbar } \left [ H ^ { \kappa }, a _ { i } ^ { \pm } \right ].$$

For clarity, we reinstate the spin operators { S , S x i y i } using Eq. (4) and cast the equation of motion in the form

$$\hbar { \frac { d S ^ { g } _ { j } } { d t } } & = J S \left \{ \hat { e } _ { z } ( j ) [ \hat { e } _ { z } ( j + 1 ) + \hat { e } _ { z } ( j - 1 ) ] S ^ { y } _ { j } \\ & - \hat { e } _ { y } ( j - 1 ) \hat { e } _ { y } ( j ) S ^ { y } _ { j - 1 } - \hat { e } _ { x } ( j - 1 ) \hat { e } _ { y } ( j ) S ^ { x } _ { j - 1 } - \hat { e } _ { y } ( j ) \hat { e } _ { y } ( j + 1 ) S ^ { y } _ { j + 1 } - \hat { e } _ { y } ( j ) \hat { e } _ { x } ( j + 1 ) S ^ { x } _ { j + 1 } \right \} \\ & + K n ^ { z } _ { j } S ^ { y } _ { j } + \kappa \left \{ 2 K S ( n ^ { z } _ { j } ) ^ { 2 } S ^ { y } _ { j } - 2 K S ( n ^ { y } _ { j } ) ^ { 2 } S ^ { y } _ { j } - 2 K S n ^ { x } _ { j } n ^ { y } _ { j } S ^ { x } _ { j } - K n ^ { z } _ { j } S ^ { y } _ { j } \right \},$$

$$\hbar { \frac { d S ^ { \mu } _ { j } } { d t } } & = J S \left \{ - \hat { e } _ { z } ( j ) [ \hat { e } _ { z } ( j + 1 ) + \hat { e } _ { z } ( j - 1 ) ] S ^ { z } _ { j } \\ & + \hat { e } _ { x } ( j - 1 ) \hat { e } _ { x } ( j ) S ^ { z } _ { j - 1 } + \hat { e } _ { y } ( j - 1 ) \hat { e } _ { x } ( j ) S ^ { y } _ { j - 1 } + \hat { e } _ { x } ( j ) \hat { e } _ { x } ( j + 1 ) S ^ { z } _ { j + 1 } + \hat { e } _ { x } ( j ) \hat { e } _ { y } ( j + 1 ) S ^ { y } _ { j + 1 } \right \} \\ & - K n ^ { z } _ { j } S ^ { z } _ { j } - \kappa \left \{ 2 K S ( n ^ { z } _ { j } ) ^ { 2 } S ^ { z } _ { j } - 2 K S ( n ^ { z } _ { j } ) ^ { 2 } S ^ { y } _ { j } - 2 K S n ^ { z } _ { j } n ^ { y } _ { j } S ^ { y } _ { j } + K n ^ { z } _ { j } S ^ { y } _ { j } \right \}.$$

Eqs. (8) and (9) are identical to the linearized classical Landau-Lifshitz equations expressed in the local coordinate system { ˆ ( e x j ) , e ˆ ( y j ) , e ˆ ( z j ) } . Since we are only interested in studying how the intrinsic disorder affects the transport properties of the system, we have not included a Gilbert damping term.

## C. Scattering problem and solution ansatz

Eqs. (8) and (9) can be solved numerically in the classical regime, where we treat the spin operators as classical vectors. The spin-wave solutions are the normal modes of the system and precess with the same frequency ω . Therefore, we can factorize out the time dependence of the spin operators as e -iωt .

Deep inside the leads, the spins at neighboring sites are pointing in the z -direction; see Fig. 1. This considerably simplifies the equations of motions in the leads

$$\mu \quad & \quad \mu \\ & - i \hbar { \omega } S _ { j } ^ { x } = J S \left ( 2 S _ { j } ^ { y } - S _ { j - 1 } ^ { y } - S _ { j + 1 } ^ { y } \right ) \quad \text{                                                                                                                                                                                                         in                                                                                                                                                                                                       
                                                                                                                                                                                                        <t h                                                                                                                                                                                                        </t h                                                                                                                                                                                                     
  & + K \left ( 1 + \kappa \left ( 2 S - 1 \right ) \right ) S _ { j } ^ { y }, \quad \text{                                                                                                                                                                                                   <t h  
                                                                                                                                                                                                    
   & - i \hbar { \omega } S _ { j } ^ { y } = J S \left ( - 2 S _ { j } ^ { x } + S _ { j - 1 } ^ { x } + S _ { j + 1 } ^ { x } \right ) \quad \text{                                                                                                                                                                                              
         & - K \left ( 1 + \kappa \left ( 2 S - 1 \right ) \right ) S _ { j } ^ { x }.$$

The system behaves as a ferromagnet with an external field or intrinsic anisotropy in the z -direction. The solutions are therefore circularly polarized plane waves traveling with a fixed wavenumber q and frequency ω . The dispersion relation can be determined by substituting the ansatze S j x = e iqjd and S j y = -ie iqjd into Eq. (10). The result is

$$\epsilon = \hbar { \omega } = 2 J S \left ( 1 - \cos \ q d \right ) + K \left ( 1 + \kappa \left ( 2 S - 1 \right ) \right ). \ \left ( 1 1 \right )$$

Let us now formulate the scattering problem. Deep inside the regions i ) and iii ) in Fig. 1, we know that the solution must have the form

$$S _ { j } ^ { x } & = e ^ { i q j d } + r _ { x } e ^ { - i q j d }, \\ S _ { j } ^ { y } & = - i \left ( e ^ { i q j d } + r _ { y } e ^ { - i q j d } \right ),$$

and

$$S _ { j } ^ { x } & = t _ { x } e ^ { i q j d }, \\ S _ { j } ^ { y } & = - i t _ { y } e ^ { - i q j d },$$

respectively. Inside region ii ), we know that the spin components must satisfy Eqs. (8) and (9). Using the ansatze as boundary conditions, we have found a finite set of algebraic equations that we can solve numerically to determine the reflection and transmission amplitudes { r x , r y , t x , t y } as functions of glyph[epsilon1] .

## D. Spin current and conductance

Once we know the reflection and transmission amplitude, we can calculate the spin conductance of the disordered magnet utilizing the Landauer-B¨ttiker formalism u in the linear response regime. In this section, we derive the expression for spin conductance.

In the leads, the Hamiltonian in Eq. (5) simplifies to

$$H ^ { \kappa } & = \sum _ { i } \left ( 2 J S + K \left ( 1 + \kappa \left ( 2 S - 1 \right ) \right ) \right ) a _ { i } ^ { \dagger } a _ { i } \\ & \quad - J S \left ( a _ { j } a _ { j + 1 } ^ { \dagger } + a _ { j } ^ { \dagger } a _ { j + 1 } \right ).$$

From the equation of motion,

$$\begin{smallmatrix}. \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text
 \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{
 \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text
\text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ }
\text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ }\text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \end{smallmatrix}$$

where N i = a a † i i is the number operator, and we can extract the spin current from site j to j +1 as

$$I _ { j, j + 1 } = i J S \left ( a _ { j + 1 } ^ { \dagger } a _ { j } - a _ { j } ^ { \dagger } a _ { j + 1 } \right ). \quad \ \ ( 1 6 )$$

Now consider the situation in Fig. 2, where two reservoirs in thermodynamic equilibrium are attached to two leads with a scattering region between them. If the spin accumulation in the left reservoir µ L is greater than the spin accumulation in the right reservoir µ R , the spin current in Eq. (16) will flow from the left to the right reservoir. We define the operators α L R , ( q ) and β L R , ( q ) injecting and removing magnons with wavenumbers q into the leads, respectively. The relationship between these operators is given by the scattering matrix

$$\begin{pmatrix} \beta _ { L } ( q ) \\ \beta _ { R } ( q ) \end{pmatrix} = \begin{pmatrix} r & t ^ { \prime } \\ t & r ^ { \prime } \end{pmatrix} \begin{pmatrix} \alpha _ { L } ( q ) \\ \alpha _ { R } ( q ) \end{pmatrix}, \quad \quad ( 1 7 ) \quad r _ { y } \quad \vdots \\ \underset { \text{insi} } { \text{insi} }$$

where r ( r ′ ) and t ( t ′ ) are the reflection and transmission amplitudes, respectively, for a spin wave originating from the left (right) lead.

FIG. 2: A disordered magnet (blue) is sandwiched between two ferromagnetic leads (red). The leads are connected to two spin reservoirs (green) with spin accumulations µ L and µ R . The reservoirs are in thermodynamic equilibrium such that the magnon population is characterized by the Bose-Einstein distribution. A spin current is induced when there is a nonzero spin bias δµ = µ L -µ R .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00817v3-figure-f2c2e26c-1128-42ff-a80d-333195643f96.jpg)

In the left lead, we can express a j as [95]

$$a _ { j } = \int _ { 0 } ^ { \pi / d } \frac { d q } { 2 \pi / d } \left [ e ^ { i q j d } \alpha _ { \mathrm L } ( q ) + e ^ { i q j d } \beta _ { \mathrm L } ( q ) \right ]. \quad ( 1 8 ) \quad \text{This} \quad \text{the $f$} \\ \text{div}$$

If we substitute Eq. (18) and its complex conjugate into Eq. (16) and utilize that the leads are in thermal equilibrium with the reservoirs such that 〈 α † L R , ( q 1 ) α L R , ( q 2 ) 〉 = 2 π d δ q ( 1 -q 2 ) f L R , ( q 1 ), we find that

$$\langle I _ { j, j + 1 } \rangle = \frac { 1 } { 2 \pi } \int _ { \epsilon _ { \min } } ^ { \epsilon _ { \max } } d \epsilon \, T ( \epsilon ) \left ( f _ { \text{L} } ( \epsilon ) - f _ { \text{R} } ( \epsilon ) \right ). \quad ( 1 9 )$$

In this expression, f L R , ( glyph[epsilon1] ) represents the Bose-Einstein distributions in the left and right reservoirs, respectively, and T glyph[epsilon1] ( ) = | t | 2 . The integration limits are obtained from Eq. (11).

Assume that the spin accumulation in the left lead is µ L = µ + δµ and that the spin accumulation in the right lead is µ R = µ , where δµ/µ glyph[lessmuch] 1. We find that in the linear response, the spin conductance is given by

$$G = \frac { 1 } { 2 \pi } \int _ { \tilde { \epsilon } _ { \min } } ^ { \tilde { \epsilon } _ { \max } } d \tilde { \epsilon } \, T ( \tilde { \epsilon } ) \left ( - \frac { d f } { d \tilde { \epsilon } } \right ). \quad \ \ ( 2 0 ) \quad \text{ and}$$

This result can also be derived using Green's functions [96]. In Eq. (20), we are integrating over the dimensionless energies ˜ = glyph[epsilon1] glyph[epsilon1]/J . Energies outside of the integration interval result in the spin waves in Eqs. (12) and (13) becoming evanescent waves that do not contribute to the spin conductance.

## IV. RESULTS AND DISCUSSION

For each realization of the system, we find that r x = r y ≡ r and that t x = t y ≡ t , reflecting the fact that inside the leads, the spin waves are circularly polarized. Furthermore, we define R = | r | 2 and T = | t | 2 as the reflection and transmission probabilities, respectively, and find that R + T = 1. Since R and T depend on the realization of the system, we must perform an ensemble average 〈 . . . 〉 to obtain physically meaningful quantities. In our calculations, we used 10 3 different realizations for the random vectors n i . In Fig. 3, we plotted 〈 T 〉 and 〈 ln T 〉 as a function of ˜ for different values of glyph[epsilon1] K/J and a fixed system length L = N rand d . In the remainder of this paper we set d = 1 for convenience.

## A. Transmission probability

As the system becomes more disordered, the transmission probability decreases for both the RFM and the RAM. However, Fig. 3 demonstrates that the quantitative behavior of the localization is significantly different in the two models. In both models, as K/J increases, the maxima 〈 T 〉 max and 〈 ln T 〉 max shift towards higher ˜, but glyph[epsilon1] in the RAM, this shift is greater than that in the RFM. In addition, the peak in the transmission probability is wider in the RAM compared to the RFM for small K/J . Thus, a broader range of spin waves can pass through the RAM compared to the RFM in the limit of weak disorder.

FIG. 3: The behaviors of 〈 T 〉 and 〈 ln T 〉 as a function of ˜ glyph[epsilon1] and K/J .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00817v3-figure-f21013af-9d9c-4889-8115-1ba76df2762e.jpg)

We can understand the difference in width from the Hamiltonian in Eq. (1). In the RAM, the term causing disorder is ( S i · n i ) 2 ; thus, the spin S i wants to point either parallel or antiparallel to n i . The spin is also coupled to its neighbors through the exchange interaction. Therefore, in the RAM, whether the spin S i chooses to point parallel or antiparallel to n i depends on the neighboring spins. Meanwhile, in the RFM, the term causing disorder is S i · n i , and the spin wants to only point parallel to n i . The ability to select whether to point parallel or antiparallel to n i leads to the spin chains in the RAM being less disordered than the spin chains in the RFM, which in turn leads to a broader peak in the transmission probability.

## B. Self-averaging

In disordered systems, certain quantities are not selfaveraging in the thermodynamic limit. This is well known in disordered fermionic systems and is expected to be a general feature of a broad spectrum of disordered systems [97]. A test to determine whether a quantity O is self-averaging is to check whether the relative variance RV O = Var( O / O ) 〈 〉 2 vanishes (or is sufficiently small) in the limit L - ∞ → . For the fermionic 1D Anderson model with on-site disorder, one finds that the transmission probability and hence the conductance are not self-averaging [98]. In two and three dimensions, one finds that the logarithms ln( T , ) ln( G ) are self-averaging such that RV ln G ∼ L -D ( D = 2 3) [97, 99]. , In one dimension at finite temperature, one finds that ln G is only marginally self-averaging because RV ln G decays logarithmically with L [99-101].

FIG. 4: The length dependence of the relative variances RV , RV G ln G , RV , and RV T ln T for the RFM and the RAM. The strength of disorder is K/J = 0 4. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00817v3-figure-663aeefd-1f9f-49aa-8dc4-eb46594c2046.jpg)

As expected, we find similar results in this work. Fig. 4 shows that the relative variances RV T and RV G increase with the length of the system. In addition, the relative variances RV ln T and RV ln G decrease with the length of the system. Hence, as the length of the system increases, the fluctuations in T and G become much greater than the corresponding expectation values, meaning that they are not representative variables in the thermodynamic limit. Therefore, we use { ln( T , ) ln( G ) } , rather than { T, G } , to calculate the localization lengths of the system.

## C. The localization length

In this work, it is natural to define two types of localization lengths. The first is based on the maximum of the transmission amplitude 〈 ln T 〉 max in Fig. 3. The second is based on the conductance 〈 ln G 〉 . We refer to these localization lengths as ˜ L ln T and ˜ L ln G , respectively.

In Fig. 5, we plotted 〈 ln T 〉 max and 〈 ln G 〉 as a function of the system length L for a fixed K/J and temperature ˜ T = kT/J . We have performed a curve fit with the functions

$$\langle \ln T \rangle _ { \max } = \frac { L } { \tilde { L } _ { \ln T } } + A \, \text{ \quad \ \ } ( 2 1 )$$

and

$$\langle \ln G \rangle = \frac { L } { \tilde { L } _ { \ln G } } + B \, \text{ \quad \ \ } \text{(2)}$$

such that the localization length can be extracted as the gradient of the straight lines in Fig. 5. In this particular case, we found a coefficient of determination R 2 with the value R 2 = 0 95 indicating a good fit. . To determine

FIG. 5: The length dependence of 〈 ln G 〉 and 〈 ln T 〉 max for the RFM can be approximated with a linear fit. The strength of disorder is K/J = 0 4, and the temperature is . ˜ = 0 05. T .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00817v3-figure-c49ccacd-7dfe-4b31-b822-be950dbacee7.jpg)

the localization lengths as a function of K/J , we performed straight line curve fits for graphs such as those

found in Fig. 5 but with different K/J and ˜ . T In all cases, we found that the coefficient of determination was in the range (0 9 . , 1) and that the average coefficient was 〈 R 2 〉 = 0 95, indicating reasonably good straight-line fits. . By then calculating the gradient of these straight lines, we can estimate the localization lengths as a function of K/J .

## 1. Localization length from transmission

In Fig. 6, we plot the localization length ˜ L ln T and the 95% -confidence interval for the RAM and RFM, respectively. In both cases, we have performed a curve fit with the function

$$\tilde { L } _ { \ln T } = \eta \left ( \frac { K } { J } \right ) ^ { \nu } + \xi. \quad \quad ( 2 3 ) \quad \stackrel { ( 0. } { T }. \\ \text{ish}$$

The parameters (with confidence intervals) are displayed in Table I. Similar to fermionic systems [102], we find that the localization length decays monotonically as a power law as we increase the strength of disorder. Our result can be made more universal by introducing the exchange length such that

$$\tilde { L } _ { \ln T } = \eta \left ( l _ { \text{ex} } \right ) ^ { - 2 \nu } + \xi. \text{ \quad \ \ \ } \left ( 2 4 \right ) \text{ \ \ } \text{$sys$}$$

Note that for weak disorder, the localization length is greater in the RAM than in the RFM. This is a consequence of the fact that the spin chains are less disordered in the RAM compared to the RFM, as we discussed at the end of section IV A.

FIG. 6: The behavior of ˜ L ln T as a function of K/J for the RAM and the RFM, respectively. The line represents the numerical fit in Eq. (23), the dashed lines represent the 95% confidence interval, and the points with error bars represent the localization length calculated from Eqs. (8) and (9) with standard error.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00817v3-figure-230e0d68-6dd8-43dd-9512-ea3f6d7d7de0.jpg)

TABLE I: The numerical values of the parameters in Eq. (23) for the RFM and the RAM. The brackets ( . . . ) give the 95% confidence interval.

| RFM                         | RAM                             |
|-----------------------------|---------------------------------|
| η 1 . 3 (0 . 7 , 2 . 0)     | 0 . 2 (0 . 1 , 0 . 3)           |
| ν - 1 . 2 ( - 1 . 4 , - 1 . | 0) - 2 . 2 ( - 2 . 4 , - 2 . 0) |
| ξ 1 . 1 (0 . 3 , 1 . 9)     | 1 . 6 (1 . 0 , 2 . 1)           |

## 2. Localization length from conductance

In Fig. (7), we plot ˜ L ln G as a function of K/J for different temperatures ˜ . T There is an interval K/J ≈ (0 5 . , 2) where the localization length increases for small ˜ . T Furthermore, for sufficiently large ˜ , this interval vanT ishes such that the localization length decays monotonically for all K/J . This nontrivial behavior arises because there is a competition between the temperature dependence of the broadening function -df/dglyph[epsilon1] ˜and the disorder dependence of the transmission probability T glyph[epsilon1] (˜) in Eq. (20). As the temperature increases, the broadening function excites an increasing number of magnons, which in turn leads to a greater conductance. Meanwhile, as the system becomes more disordered, the transmission probability T glyph[epsilon1] ( ˜ ) decreases, resulting in a smaller conductance. On the interval K/J ≈ (0 5 . , 2), the increase in conductance due to temperature is greater than the decrease in conductance due to disorder, which results in an increase in localization length. Furthermore, in this interval, the localization length is comparable to the lattice spacing d , which means that there may be complicated microscopic details of the model that may further enhance this effect.

Due to the complicated temperature and disorder dependence, it is numerically challenging to determine a closed formula such as the one in Eq. (23) for the localization length ˜ L ln G . However, in the weak-disorder limit K/J -→ 0, it is reasonable to assume that the localization length decays as a power law of the form

$$L _ { \ln G } \sim \left ( \frac { K } { J } \right ) ^ { \gamma } = ( l _ { \text{ex} } ) ^ { - 2 \gamma }$$

where γ is the critical exponent. Fig. 8 shows the result of such a curve fit for the RFM and the RAM for different temperatures. The corresponding critical exponents γ RFM and γ RAM are given in Table II.

In our simulations, we kept the temperature below the Curie temperature ˜ = 1, where the temperature fluctuaT tions of the spins are negligible. For temperatures ˜ T ≈ 1, there will be additional temperature-induced disorder. This issue has previously been investigated [73, 103, 104] by including a temperature-dependent stochastic field in the effective field H κ i in Eq. (2), and it was found that temperature fluctuations shorten the localization length and enhance the Anderson localization.

FIG. 7: The temperature dependence of ˜ L ln G for strongly disordered magnetic insulators.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00817v3-figure-aac999b5-79c7-4d02-a20b-91a5032cc565.jpg)

TABLE II: The numerical values of the critical exponent for the RFM and the RAM for different temperatures. The brackets ( . . . ) provide the 95% confidence interval.

|   ˜ T | γ RFM γ RAM                                               |
|-------|-----------------------------------------------------------|
|  0.05 | - 3 . 8 ( - 4 . 3 , - 3 . 4) - 1 . 9 ( - 2 . 3 , - 1 . 6) |
|  0.1  | - 3 . 7 ( - 4 . 2 , - 3 . 2) - 1 . 7 ( - 2 . 0 , - 1 . 4) |
|  0.2  | - 3 . 6 ( - 4 . 1 , - 3 . 1) - 1 . 3 ( - 1 . 6 , - 1 . 0) |
|  0.5  | - 3 . 3 ( - 3 . 8 , - 2 . 8) - 0 . 7 ( - 1 . 0 , - 0 . 4) |

## V. SUMMARY AND CONCLUSIONS

In this paper, we have applied the Landauer-B¨ttiker u formalism to noncollinear disordered magnetic insulators. We have considered both amorphous magnets and magnets with spin disorder modeled by the RAM and the RFM respectively. We calculated the self-averaging quantities 〈 ln T 〉 and 〈 ln G 〉 as a function of system length L for a broad range of disorder strengths K/J . Consistent with the literature, we found evidence for Anderson localization such that 〈 ln T 〉 and 〈 ln G 〉 were linear functions of the system length L . This allowed us to define two localization lengths ˜ L ln G and ˜ L ln T based on the conductance and the maximum transmission probability, respectively. In the limit of weak disorder, the localization lengths obeyed power laws, and we calculated the

FIG. 8: The temperature dependence of ˜ L ln G in the limit of weak disorder K/J -→ 0 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00817v3-figure-185a5a50-cad1-4b3d-8da7-4196d2dcdc59.jpg)

relevant critical exponents. We expect our results to be general because they are expressed through the universal exchange length l ex .

We found that the Anderson localization is more prominent in the RFM than in the RAM. The reason for this result is that the competition between the exchange interaction and the disorder term leads to more disordered spin chains in the RFM than in the RAM. The spin chains in the RAM are less disordered because the disorder arises from a random anisotropy, where the spin can point either parallel or antiparallel to the anisotropy with the same energy cost. Whether the spin points parallel or antiparallel to the anisotropy is determined by the neighboring spins through the exchange interaction, and consequently, the configuration with the least disorder will be chosen by the system.

The results obtained here are valid in the limit of quenched disorder, i.e., ˜ T glyph[lessmuch] 1, where the random field and anisotropy are temperature independent. If the temperature is close to the Curie temperature of the system, one must include temperature fluctuations in the Landau-Lifshitz equations. Such effects have been considered in other works, and it has been shown that temperature fluctuations decrease the localization length.

To experimentally verify the critical exponents obtained in this paper, we propose a setup in which a disordered magnetic nanowire is sandwiched between two normal metals. Similar setups for ordered magnets have

been considered in other works [11-25]. By applying a charge current in the left metal, the spin Hall effect generates a spin current through the disordered nanowire and into the right metal. This will give rise to a spin wave propagating through the hybrid structure and into the right metal, where the spin current is converted into a charge current via the inverse spin Hall effect.

Alternatively, we can instead sandwich a disordered magnet between two ferromagnetic leads. We can excite spin waves in the left ferromagnet by applying a microwave with the ferromagnet resonance frequency. This spin wave will then propagate through the disordered insulator and into the right ferromagnet, where the resulting spin current can be measured. By measuring the spin

- [1] V. V. Kruglyak, S. O. Demokritov, and D. Grundler, J. Phys. D: Appl. Phys. 43 , 264001 (2010).
- [2] A. V. Chumak, V. I. Vasyuchka, A. A. Serga, and B. Hillebrands, Nat. Phys. 11 , 453 EP (2015).
- [3] A. A. Serga, A. V. Chumak, and B. Hillebrands, J. Phys. D: Appl. Phys. 43 , 264002 (2010).
- [4] V. Kruglyak and R. Hicken, J. Magn. Magn. Mater. 306 , 191 (2006).
- [5] A. Qaiumzadeh, L. A. Kristiansen, and A. Brataas, Phys. Rev. B 97 , 020402(R) (2018).
- [6] E. G. Tveten, A. Qaiumzadeh, and A. Brataas, Phys. Rev. Lett. 112 , 147204 (2014).
- [7] S. Neusser and D. Grundler, Adv. Mater. 21 , 2927 (2009).
- [8] Y. Tserkovnyak, A. Brataas, and G. E. W. Bauer, Phys. Rev. Lett. 88 , 117601 (2002).
- [9] G. E. W. Bauer, E. Saitoh, and B. J. van Wees, Nat. Mater. 11 , 391 EP (2012).
- [10] J. Sinova, S. O. Valenzuela, J. Wunderlich, C. H. Back, and T. Jungwirth, Rev. Mod. Phys. 87 , 1213 (2015).
- [11] R. Karplus and J. M. Luttinger, Phys. Rev. 95 , 1154 (1954).
- [12] M. Dyakonov and V. Perel, Phys. Lett. A 35 , 459 (1971).
- [13] S. Zhang, Phys. Rev. Lett. 85 , 393 (2000).
- [14] S. Murakami, N. Nagaosa, and S.-C. Zhang, Science 301 , 1348 (2003).
- [15] J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, and A. H. MacDonald, Phys. Rev. Lett. 92 , 126603 (2004).
- [16] T. Jungwirth, J. Wunderlich, and K. Olejn´ ık, Nat. Mater. 11 , 382 EP (2012).
- [17] Y. Kajiwara, K. Harii, S. Takahashi, J. Ohe, K. Uchida, M. Mizuguchi, H. Umezawa, H. Kawai, K. Ando, K. Takanashi, S. Maekawa, and E. Saitoh, Nature 464 , 262 EP (2010).
- [18] E. Saitoh, M. Ueda, H. Miyajima, and G. Tatara, Appl. Phys. Lett. 88 , 182509 (2006).
- [19] S. O. Valenzuela and M. Tinkham, Nature 442 , 176 (2006).
- [20] T. Kimura, Y. Otani, T. Sato, S. Takahashi, and S. Maekawa, Phys. Rev. Lett. 98 , 156601 (2007).
- [21] J.-i. Ohe, A. Takeuchi, and G. Tatara, Phys. Rev. Lett. 99 , 266603 (2007).

current propagating through the disordered region, one should be able to characterize the localization length in terms of the critical exponents of the system.

## ACKNOWLEDGMENTS

The authors would like to express their thanks to R. A. Duine for fruitful discussions of the subject.

This research was supported by the European Research Council via Advanced Grant No. 669442 'Insulatronics' and the Research Council of Norway through its Centres of Excellence funding scheme, Project No. 262633, 'QuSpin'.

- [22] T. Seki, Y. Hasegawa, S. Mitani, S. Takahashi, H. Imamura, S. Maekawa, J. Nitta, and K. Takanashi, Nat. Mater. 7 , 125 EP (2008).
- [23] K. Ando, Y. Kajiwara, S. Takahashi, S. Maekawa, K. Takemoto, M. Takatsu, and E. Saitoh, Phys. Rev. B 78 , 014413 (2008).
- [24] K. Ando, M. Morikawa, T. Trypiniotis, Y. Fujikawa, C. H. W. Barnes, and E. Saitoh, Appl. Phys. Lett. 96 , 082502 (2010).
- [25] L. J. Cornelissen, J. Liu, R. A. Duine, J. B. Youssef, and B. J. van Wees, Nat. Phys. 11 , 1022 EP (2015).
- [26] H. Ochoa, R. Zarzuela, and Y. Tserkovnyak, Phys. Rev. B 98 , 054424 (2018).
- [27] Y. Tserkovnyak and H. Ochoa, Phys. Rev. B 96 , 100402(R) (2017).
- [28] N. Arakawa and J.-i. Ohe, Phys. Rev. B 98 , 014421 (2018).
- [29] M. Evers, C. A. M¨ller, u and U. Nowak, Phys. Rev. B 97 , 184423 (2018).
- [30] H. V. Gomonay, R. V. Kunitsyn, and V. M. Loktev, Phys. Rev. B 85 , 134446 (2012).
- [31] P. Buczek, S. Thomas, A. Marmodoro, N. Buczek, X. Zubizarreta, M. Hoffmann, T. Balashov, W. Wulfhekel, K. Zakeri, and A. Ernst, J. Phys: Condens. Matter 30 , 423001 (2018).
- [32] D. Wesenberg, T. Liu, D. Balzar, M. Wu, and B. L. Zink, Nat. Phys. 13 , 987 (2017).
- [33] R. Lebrun, A. Ross, S. A. Bender, A. Qaiumzadeh, L. Baldrati, J. Cramer, A. Brataas, R. A. Duine, and M. Kl¨ aui, Nature 561 , 222 (2018).
- [34] B. I. Halperin and W. M. Saslow, Phys. Rev. B 16 , 2154 (1977).
- [35] C. L. Henley, H. Sompolinsky, and B. I. Halperin, Phys. Rev. B 25 , 5849 (1982).
- [36] P. W. Anderson, B. I. Halperin, and C. M. Varma, Philos. Mag. 25 , 1 (1972).
- [37] P. Anderson, Mater. Res. Bull. 5 , 549 (1970).
- [38] I. Y. Korenblit and E. F. Shender, Sov. Phys. Usp. 21 , 832 (1978).
- [39] R. N. Costa Filho and M. G. Cottam, J. Appl. Phys. 81 , 5749 (1997).
- [40] R. Bhargava and D. Kumar, Phys. Rev. B 19 , 2764 (1979).
- [41] R. Weber, Zeitschrift f¨r Physik A Hadrons and nuclei u

- 223 , 299 (1969).
- [42] S. E. Barnes, J. Phys. F: Met. Phys. 11 , L249 (1981).
- [43] N. Arakawa and J.-i. Ohe, Phys. Rev. B 97 , 020407(R) (2018).
- [44] C. De Dominicis and I. Giardina, Random Fields and Spin Glasses: A Field Theory Approach (Cambridge University Press, 2006).
- [45] G. Herzer, in Properties and Applications of Nanocrystalline Alloys from Amorphous Precursors , edited by B. Idzikowski, P. ˇ vec, S and M. Miglierini (Springer Netherlands, Dordrecht, 2005) pp. 15-34.
- [46] R. Alben, J. J. Becker, and M. C. Chi, J. Appl. Phys. 49 , 1653 (1978).
- [47] R. Harris, M. Plischke, and M. J. Zuckermann, Phys. Rev. Lett. 31 , 160 (1973).
- [48] J. Herbrych and J. Kokalj, Phys. Rev. B 95 , 125129 (2017).
- [49] A. A. Fedorenko, Phys. Rev. E 86 , 021131 (2012).
- [50] S. F. Edwards and P. W. Anderson, J. Phys. F: Met. Phys. 5 , 965 (1975).
- [51] D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35 , 1792 (1975).
- [52] D. Mattis, Phys. Lett. A 56 , 421 (1976).
- [53] A. P. Young, Spin Glasses and Random Fields (WORLD SCIENTIFIC, 1997).
- [54] G. Parisi, Phys. Rev. Lett. 50 , 1946 (1983).
- [55] J. Smith, A. Lee, P. Richerme, B. Neyenhuis, P. W. Hess, P. Hauke, M. Heyl, D. A. Huse, and C. Monroe, Nat. Phys. 12 , 907 EP (2016).
- [56] D. Belanger and A. Young, J. Magn. Magn. Mater. 100 , 272 (1991).
- [57] N. G. Fytas, V. Mart´ ın-Mayor, M. Picco, and N. Sourlas, Phys. Rev. Lett. 116 , 227201 (2016).
- [58] G. Tarjus and M. Tissier, Phys. Rev. Lett. 93 , 267008 (2004).
- [59] J. Z. Imbrie, Phys. Rev. Lett. 53 , 1747 (1984).
- [60] A. Malakis and N. G. Fytas, Phys. Rev. E 73 , 016109 (2006).
- [61] S. Fishman and A. Aharony, J. Phys. C: Sol. Stat. Phys. 12 , L729 (1979).
- [62] J. L. Cardy, Phys. Rev. B 29 , 505 (1984).
- [63] D. J. Luitz, N. Laflorencie, and F. Alet, Phys. Rev. B 91 , 081103(R) (2015).
- [64] A. Pal and D. A. Huse, Phys. Rev. B 82 , 174411 (2010).
- [65] A. D. Luca and A. Scardicchio, Eur. Phys. Lett. 101 , 37003 (2013).
- [66] A. Nanduri, H. Kim, and D. A. Huse, Phys. Rev. B 90 , 064201 (2014).
- [67] M. ˇ nidariˇ, T. c. v. Prosen, Z c and P. Prelovˇ sek, Phys. Rev. B 77 , 064426 (2008).
- [68] I. Khait, S. Gazit, N. Y. Yao, and A. Auerbach, Phys. Rev. B 93 , 224205 (2016).
- [69] P. W. Anderson, Phys. Rev. 109 , 1492 (1958).
- [70] D. Vollhardt and P. W¨lfle, Phys. Rev. Lett. o 48 , 699 (1982).
- [71] A. J. Bray and M. A. Moore, Phys. Rev. Lett. 41 , 1068 (1978).
- [72] G. Parisi, Phys. Rev. Lett. 43 , 1754 (1979).
- [73] M. Evers, C. A. M¨ller, u and U. Nowak, Phys. Rev. B 92 , 014411 (2015).
- [74] E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, Phys. Rev. Lett. 42 , 673 (1979).
- [75] V. Cannella and J. A. Mydosh, Phys. Rev. B 6 , 4220 (1972).
- [76] R. Landauer, IBM. J. Res. Dev. 1 , 223 (1957).
- [77] C. W. J. Beenakker, Rev. Mod. Phys. 69 , 731 (1997).
- [78] K. R. Sapkota, R. Dulal, B. R. Dahal, I. L. Pegg, and J. Philip, J. Vac. Sci. Tech. B 33 , 051807 (2015).
- [79] R. R. Galazka, S. Nagata, and P. H. Keesom, Phys. Rev. B 22 , 3344 (1980).
- [80] S. Venugopalan, A. Petrou, R. R. Galazka, A. K. Ramdas, and S. Rodriguez, Phys. Rev. B 25 , 2681 (1982).
- [81] T. Giebultowicz, B. Lebech, B. Buras, W. Minor, H. Kepa, and R. R. Galazka, J. Appl. Phys. 55 , 2305 (1984).
- [82] T. Kaneyoshi, Phys. B: Condens. Mater. 462 , 34 (2015).
- [83] S. F. Edwards and F. Tanaka, J. Phys. F: Met. Phys. 10 , 2471 (1980).
- [84] B. Drossel and M. A. Moore, Phys. Rev. B 70 , 064412 (2004).
- [85] Z. Ma, J. Wang, Z.-Y. Dong, J. Zhang, S. Li, S.-H. Zheng, Y. Yu, W. Wang, L. Che, K. Ran, S. Bao, Z. Cai, P. Cerm´ ak, ˇ A. Schneidewind, S. Yano, J. S. Gardner, X. Lu, S.-L. Yu, J.-M. Liu, S. Li, J.-X. Li, and J. Wen, Phys. Rev. Lett. 120 , 087201 (2018).
- [86] T. Castellani and A. Cavagna, J. Stat. Mech. 2005 , P05012 (2005).
- [87] M. Mezard, G. Parisi, and M. Virasoro, Spin Glass Theory and Beyond (WORLD SCIENTIFIC, 1986).
- [88] A. J. Bray and M. A. Moore, J. Phys. C: Sol. Stat. Phys. 14 , 1313 (1981).
- [89] T. Li, Phys. Rev. B 24 , 6579 (1981).
- [90] C. M. Newman and D. L. Stein, Phys. Rev. E 60 , 5244 (1999).
- [91] Z. Burda, A. Krzywicki, O. C. Martin, and Z. Tabor, Phys. Rev. E 73 , 036110 (2006).
- [92] Z. Burda, A. Krzywicki, and O. C. Martin, Phys. Rev. E 76 , 051107 (2007).
- [93] B. Waclaw and Z. Burda, Phys. Rev. E 77 , 041114 (2008).
- [94] A. Heuer, J. Phys: Condens. Matter 20 , 373101 (2008).
- [95] An analogous expression can be used for the right lead.
- [96] B. Wang, J. Wang, J. Wang, and D. Y. Xing, Phys. Rev. B 69 , 174403 (2004).
- [97] A. Aharony and A. B. Harris, Phys. Rev. Lett. 77 , 3700 (1996).
- [98] E. Abrahams, 50 Years of Anderson Localization (WORLD SCIENTIFIC, 2010).
- [99] R. A. Serota, R. K. Kalia, and P. A. Lee, Phys. Rev. B 33 , 8441 (1986).
- [100] M. A. Aamir, P. Karnatak, A. Jayaraman, T. P. Sai, T. V. Ramakrishnan, R. Sensarma, and A. Ghosh, Phys. Rev. Lett. 121 , 136806 (2018).
- [101] P. A. Lee, Phys. Rev. Lett. 53 , 2042 (1984).
- [102] J. Billy, V. Josse, Z. Zuo, A. Bernard, B. Hambrecht, P. Lugan, D. Cl´ment, L. Sanchez-Palencia, P. Bouyer, e and A. Aspect, Nature 453 , 891 EP (2008).
- [103] U. Ritzmann, D. Hinzke, and U. Nowak, Phys. Rev. B 89 , 024409 (2014).
- [104] A. A. Starikov, Y. Liu, Z. Yuan, and P. J. Kelly, Phys. Rev. B 97 , 214415 (2018).
Martin Fonnum Jakobsen Alireza Qaiumzadeh Arne Brataas
Scattering theory of transport through disordered magnets
We present a scattering theory of transport through noncollinear disordered
magnetic insulators. For concreteness, we study and compare the random field
model (RFM) and the random anisotropy model (RAM). The RFM and RAM are used to
model random spin disorder systems and amorphous materials, respectively. We
utilize the Landauer-Buttiker formalism to compute the transmission probability
and spin conductance of one-dimensional disordered spin chains. The RFM and the
RAM both exhibit Anderson localization, which means that the transmission
probability and spin conductance decay exponentially with the system length. We
define two localization lengths based on the transmission probability and the
spin conductance, respectively. Next, we numerically determine the relationship
between the localization lengths and the strength of the disorder. In the limit
of weak disorder, we find that the localization lengths obey power laws and
determine the critical exponents. Our results are expressed via the universal
exchange length and are therefore expected to be general.