Dephasing-Assisted Macrospin Transport

Transport phenomena are ubiquitous in physics, and it is generally understood that the environmental disorder and noise deteriorates the transfer of excitations. There are, however, cases in which transport can be enhanced by fluctuations. In the present work, we show, by means of micromagnetics simulations, that transport efficiency in a chain of classical macrospins can be greatly increased by an optimal level of dephasing noise. We also demonstrate the same effect in a simplified model, the dissipative Discrete Nonlinear Schrödinger equation, subject to phase noise. Our results point towards the realization of a large class of magnonics and spintronics devices, where disorder and noise can be used to enhance spin-dependent transport efficiency.


Introduction
Dissipation and fluctuations due to interaction between a system and the environment play a crucial role in transport and relaxation processes.The understanding of the resulting non-equilibrium states is one of the main aims of contemporary statistical mechanics [1].In addition, small open systems display many special properties with respect to their macroscopic counterparts due to finite-size fluctuations, reduced dimensionality, boundary effects, disorder, etc. [2,3].Among countless examples of applications, we may mention heat transfer in nanoscale systems and nanophononics [4].In addition, in the last three decades, magnonics and spintronics emerged as new research fields that aim at using the spin degree of freedom in electronics devices [5,6] in the form of spin waves and spin dependent electrical currents.This field has grown enormously in recent years, since it offers promising opportunities for reliable ultrafast nano-size electronics that can be controlled by a variety of means, including electrical currents, magnetic fields, and temperature gradients.
Spin-caloritronics [7][8][9] concerns precisely the coupled spin/heat transport in systems with non-uniform temperature.The wide interest in these types of setups is due to their potential for energy efficient electronics, where heat flows can be used to control information [10,11].However, in magnonic devices, especially arrays of nano disks and spin-transfer nano oscillators [12], the sample-to-sample variability and the noise from the environment are deleterious for the synchronization and transport performance.
On general grounds, there are many situations in which stochastic fluctuations may play a constructive role in enhancing the response of nonlinear systems to an external coherent driving.Known examples are the enhancement of the decay time from a meta-stable state (noise-enhanced stability [13]), the synchronization with a weak periodic input signal (stochastic resonance [14]), the regularizaton of the response at a given optimal noise intensity like coherence resonance [15]) and resonant activation [16].
In the context of quantum transport, the interplay between coherence, structural disorder and noise is a prominent topic of investigation [17][18][19][20].A decade ago, it was shown that transport in quantum dissipative networks can be enhanced by the presence of pure dephasing noise [21,22].This counter-intuitive phenomenon occurs in disordered systems, where different local frequencies suppress coherent transport, and a certain level of dephasing noise can broaden the resonance lines of neighboring sites, allowing for the population transfer.Moreover, the effect appears as a general feature of quantum systems, which can be described using the language of (Markovian or not) quantum master equations.These discoveries fostered an intense research activity, aimed at understanding and exploiting dephasing in a variety of systems, including quantum dots [23], optical fibers [24], photosynthetic reactions [25][26][27] and other biological systems [28].Recently, noise-assisted transport has been studied experimentally in a chain of trapped atomic ions [29].
In this paper, we show that dephasing-assisted spin transport occurs in a chain of classical coupled (macro)spins governed by the Landau-Lifshitz-Gilbert (LLG) equation.The system is described in Section 2.In this context the effect emerges from the competition between phase coherence and phase noise and does not require intrinsically quantum-mechanical effects.The effect will be first demonstrated in Section 3 by micromagnetic simulations of the full LLG equation.To get further insights, we give in Section 4 a reduced effective description through the non-equilibrium discrete nonlinear Schrödinger equation (DNLS) [30][31][32][33] that constitutes an approximation of the LLG equation for small amplitude deviations of the spins from their equilibrium position [12].

The System: Microscopic and Macroscopic Descriptions
The system studied here consists of a chain of N = 10 Permalloy (Py) nano disks coupled through dipolar magnetic interaction, shown in Figure 1.Hereafter, we briefly describe the dynamics of the chain; we refer to References [12,34,35] for a thorough discussion.The dynamics of the magnetization M(r, t) depends on the position r inside each of the n = 1, ..., N disks and is described by the LLG equation of motion: The first term in Equation ( 1), proportional to the gyromagnetic ratio γ, describes the precession of the magnetization vector M(r, t) around the effective field H eff , while the second term accounts for energy dissipation at a rate proportional to the phenomenological Gilbert damping parameter α.The saturation magnetization M s is the norm of the magnetization, conserved during the dynamics, which depends on the material properties and sample geometry.
The effective field is the sum of four contributions.First comes the applied field H ext = H ẑ, which defines the precession axis of the magnetization ẑ.Then, the exchange field H exc (proportional to the exchange stiffness A) is the short-range interaction that accounts for the coherent precession of the magnetization inside each disk.The third term is the dipolar field H dip , which contains contributions from volume and surface charges.The dipolar field acts as a demagnetizing field in each disk and is responsible for the coupling between the disks and the nonlinearity of the dynamics [12].Thermal fluctuations in each disk are modeled by the stochastic field Here, T is the temperature of the thermal bath and where k B is the Boltzmann constant, µ 0 the vacuum magnetic permeability, and V the elementary volume containing the magnetization vector M. In finite-element micromagnetics simulations, the latter corresponds to the volume of the mesh elements.Each component of the thermal field is modeled as a Gaussian random variable with zero average and correlation h th i (r, t), h th j (r, t ) = δ(r − r )δ(t − t )δ ij with i, j = x, y, z.
In this paper, we consider a chain of Permalloy nanodisks with thickness t = 3 nm, radius R = 20 nm, and an interlayer distance d = 3 nm.The applied field H ext = 1 T defines the precession axis of the magnetization along the ẑ direction.The exchange stiffness A = 10 −11 J/m corresponds to that of Permalloy, while the other micromagnetics parameters are M s = 0.94T/µ 0 , α = 8 × 10 −3 , and γ = 1.873 × 10 11 rad s −1 T −1 .Those parameters are taken from Reference [36].The evolution of the disk chain was performed at the microscopic level by solving Equation (1) with the NMAG software [37], using a tetrahedral finite element mesh with maximum size of 3 nm, of the order of the Permalloy exchange length.
In nano disks with radius of a few hundreds of nm, the magnetization precession typically exhibits several precession normal modes [36,38].For significantly smaller nanodisks, like in our case, one can verify that the only mode that is effectively active is the uniform one (i.e., a uniform, phase coherent magnetization precession over the disk volume) [38].In this regime the system can be reduced to an ensemble of coupled macrospins defined as Note that here V n is the volume of each disk and not the volume element that enters the diffusion constant in Equation ( 4).
It is thus convenient to introduce the complex Spin Wave (SW) amplitude [12]: By writing Equation ( 5) as ψ n = p n (t)e iφ n (t) , one can see that the phase φ n (t) describes the precession of M n in the x-y plane, while p n = |ψ n | 2 is the local SW power.From these definitions, it follows that the equilibrium solution M n = M s ẑ corresponds to a vanishing SW field with p n = 0.More in general, the macroscopic spin dynamics of the disk chain can be mapped to the evolution of a system of damped oscillators described by complex variables [35], as discussed in Section 4.
Before presenting the main results, it is useful to understand the nature of the collective excitations of the disk array.A way to evaluate empirically the spin-wave spectrum is to first determine the main frequency peaks ω n appearing in the Fourier transform of the total magnetization ∑ n M n , as done in Reference [35].From micromagnetics simulations with the parameters given above, we found that there are five main components at frequencies ω 1 = 18.0, ω 2 = 19.7,ω 3 = 20.5, ω 4 = 21.8, and ω 5 = 24.0GHz.Since each disk behaves as a magnetic dipole, aligning the disks in a chain gives a structure where the intensity of the dipolar field, which controls the frequencies, is symmetric around the centre of the chain, yielding a double degeneracy of the spectrum.
To determine the spatial structure of the spin-wave modes in an effective way, we apply a uniform time-dependent magnetic field oscillating at the frequencies of the modes and measure the average profile of the spin-wave powers p n in the steady state.The simulations were performed at zero temperature.The resulting profiles are displayed in Figure 2. The profiles are symmetric with respect to the chain center.Moreover the mode with the highest-frequency ω 5 = 24.0GHz differs from the others, being localized on the two central disks.Although the disk chain is not an intrinsically disordered system, the quite well localized nature of its collective modes indicates that effective SW transport is limited to a restricted dynamical regime where the macrospins perform coherent nonlocal oscillations.On the other hand, spatially localized excitations are quickly damped by the Gilbert dissipative term in Equation (1).
Figure 2. Time-average of the SW powers of each collective mode, obtained by exciting the dynamics with a uniform time-dependent magnetic field with the frequencies of the modes until the system reaches a steady state.Simulations are performed at zero temperature.The total SW power of each profile is normalized to one for better comparison.

Micromagnetics Simulations with Dephasing Noise
Here, we are interested on the effect of pure dephasing noise on the disks dynamics.In principle, this could be physically realized by introducing a suitable random modulation of the external applied field, varying on much faster time-scales with respect to the typical precession period of the macrospins.The strength of the randomness should be the main accessible control parameter.To implement dephasing in the model, at each time step τ of the numerical solution of the LLG equation, we transform the magnetization in each disk as M n → R z (Θ n )M n , where R z is a rotation matrix by a random angle Θ n around the z axis (R z (0) being the identity matrix).The angles Θ n are independent, identically distributed Gaussian variables with zero average and a standard deviation equal to θc √ τ, where θ is an non-dimensional parameter which controls the strength of the dephasing noise and c = 1 deg ps −1/2 .By construction, such a process exactly conserves the spin powers.
Note that, at variance with previous studies [35,39], here we do not consider off-equilibrium steady states, but we study transport in a transient regime.In particular, we start from a condition where only the magnetization of the first disk is not aligned with the ẑ axis, and we monitor the evolution of the system as the energy pulse propagates through the chain, until M is aligned with the z axis in all the disks.More precisely, the initial condition is generated by tilting the magnetization of the first disk in the x direction by a given angle, leaving the other magnetizations aligned with ẑ and monitoring the time evolution for about 10 ns with a time step of 1 ps.Unless otherwise stated, all the simulations were performed setting the bath temperature T = 0 in Equation ( 3) and changing the noise strength θ.
To understand better the nature of the dephasing process, we first compare the dynamics of the system in contact with the thermal bath only, as given by Equation (3) (θ = 0, T = 0), with the one in presence of a pure dephasing noise (θ = 0, T = 0).Figure 3 compares the time evolution of the total SW power P = ∑ n p n in the two cases.From Figure 3a, one can see that the strength of the dephasing process hardly changes the relaxation to the equilibrium state P = 0, which is solely determined by intrinsic dissipation.On the other hand, thermal fluctuations affect the dynamics in a completely different fashion by making the precession amplitudes fluctuate steadily as expected.One can see that in the first case P drops to zero and the magnetization aligns with the z axis, since the dephasing conserves the total power.In the second case, the bath temperature excites the dynamics and the system thermalizes with P increasing with the bath temperature T.
Let us now turn to the effect of the dephasing on transport.Figure 4 shows the time evolution of the SW power of the last disk as a function of time for increasing θ.One can see clearly that the pulse increases up to θ ≈ 4 and then decreases.To quantify transport efficiency, we compare the total power that flows through the last disk, and we compute A value E > 1 is indicative of an enhancement or energy transfer induced by the dephasing process with respect to the noise-free case.In Figure 5a, we report the efficiency E, defined in Equation ( 6) as a function of the parameter θ; one can see that the E reaches the maximum between θ = 4 and 6, and then decreases.In other words, there is an optimal range of θ where transport is enhanced with respect to the noiseless case.In previous studies [10,39], it has been shown that the current increase at increasing phase synchronization and transport is due to phase coherence.However, in the present case, transport is mostly incoherent.This is seen in Figure 5b, which shows the Kuramoto synchronization parameter [40,41], defined as K = 1 N | ∑ N n=1 e iφ n (t) | , where • and • denote, respectively, ensemble and time averages and φ n (t) is the time-dependent phase of the n th oscillator.At variance with the efficiency, K does not display any maximum but rather decreases monotonically as dephasing destroys the phase coherence among magnetization vectors of the the disks.We remark that, as already illustrated from Figure 3, the increase of transport efficiency is not trivially due to the fact that more SW power is injected into the system, since dephasing noise conserves the total power.This indicates that a different mechanism needs to be invoked to explain the increase of the efficiency with noise.
In order to get some insight on such mechanisms, we consider the spectral content of the magnetization in the different regimes.A convenient way to analyze non-stationary signals is to use wavelet analysis in the time-frequency domain.This method allows to detect transient frequency components appearing at specific times and lasting for finite lapses of time.In this work, we computed the Gabor transform [42] of the complex magnetization of a given disk n, namely It is interesting to compare the behavior of the wavelet signal at two different positions of the chain, namely the bulk and the output region.Figure 6 shows the average density maps |G n (ω, t)| 2 for different values of the noise strength θ in the interval [1,10] for disks n = 5 and n = 10.In the whole temperature range, the wavelet signal lasts approximately for 5 ns before being damped, which is the same typical damping time observed in Figure 4 for the local power at the chain end.For small θ, we can distinguish a dominant mode excited at ω = ω 5 = 24 GHz localized on the central site n = 5 and a mode at ω 1 = 18 GHz localized on n = 10.This structure reflects the pattern of collective modes shown in Figure 2. At larger θ values, the number of transmitting modes in the frequency range around ω = 20 GHz increases.A further increment of the noise strength broadens and weakens the resonance lines, yielding the efficiency reduction described above.
A more quantitative analysis can be performed by computing the contribution of a certain frequency ω to the wavelet signal in Figure 6.Specifically, we computed the parameter g n ( ω) which is the average of |G n (ω, t)| 2 over the whole observation time and over a small frequency interval centered on ω and length equal to δω ω. Figure 7 shows the behavior of g n ( ω) for the set of characteristic frequencies ω 1 , • • • , ω 5 for sites n = 5 and n = 10 and δω = 0.25 GHz (solid lines).A large nonmonotonic contribution of ω = ω 1 is found on site n = 10 which is absent for site n = 5.For δω = 0.25 GHz (black solid line), the maximum amplitude is found for θ = 2, which is slightly smaller than the noise amplitude that maximizes the efficiency in Figure 5.By enlarging the averaging frequency window to δω = 2 GHz (black dashed line), we observe a better agreement with the efficiency curve.Therefore, we can conclude that the relevant contribution to the increase of transport efficiency of the chain is related to the excitation of a relatively narrow set of frequencies in the interval ω = 18 ± 1 GHz on the output disk.

Comparison with Coupled-Oscillators Model
For small precession angles, the dynamics of the array can be effectively described by a simplified model.Considering again the complex variables ψ n defined in Equation ( 5), and performing a suitable expansion retaining only terms up to the order α and p n ≡ |ψ n | 2 , the LLG equation for the chain can be approximated as an ensemble of coupled nonlinear oscillators which, in absence of any noise source, is of the form: [12,43] The first two terms on the right hand side of Equation ( 8) are, respectively, the nonlinear frequencies ω n (p n ) = γ|H eff • z| and damping rates Γ n (p n ).In our case, nonlinear effects are due to the dipolar (demagnetizing) field in each disk, and they are taken into account by expanding into powers of p n the frequencies and damping rates, respectively, as ω n (p n ) ≈ ω 0 n + νp n and Γ n (p n ) = αω n (p n ) [12,35].Finally, the term C nn is the interlayer (complex) coupling due to magneto-dipolar interaction.In order to keep the model as simple as possible, in the following, we will consider the simple case of a uniform nearest-neighbor interaction which amounts to retain only terms containing ψ n±1 in Equation ( 8) and set C nn ≡ J(1 + α)δ n,n±1 .Altogether, based on Equation (8) and the above simplifying assumption, we consider the following model where we also added a pure-dephasing (multiplicative) noise controlled by the parameter θ that acts independently on each oscillator, as well as where η n (t) is a Gaussian noise with zero mean and unit variance satisfying η n (t)η m (t ) = δ mn δ(t − t ) (noise can be interpreted in the usual Stratonovich sense).In the θ = 0 limit, the model is the (dissipative) Discrete Nonlinear Schrödinger equation (with site-dependent frequencies).The form of dissipation in Equation ( 9) ensures that the system reaches thermal equilibrium when put in contact with a Langevin thermal bath at temperature T [31,44].
To account for the mirror symmetry exhibited by the micromagnetic system, we consider a set of linear frequencies ω 0 n such that ω 0 n = ω 0 N−n .Upon normalizing the frequency ω 0 1 of the first oscillator to 1, we choose non-dimensional parameters J = 0.1, α = 0.008 and ν = 1.These parameters provide a quite reasonable coarse-grained description of the out-of-equilibrium dynamics of macrospin system, as shown in Reference [35].
The DNLS transport efficiency E is shown in Figure 8, as a function of the noise strength θ .Simulations are performed by initialising the DNLS chain in a state where the first site is excited with norm |ψ 1 | 2 = 1 and the rest of the chain is empty.The system is then evolved in the presence of dephasing noise until the total norm of the system goes below the threshold value 10 −5 .It was verified that this threshold provides a reasonable approximation of the integral in Equation (6).A clear maximum is visible as in the micromagnetic case, confirming that, in spite of all the crude approximations, the DNLS model captures the basic features of the original setup.Within this simplified dynamical picture, the onset of optimal transport efficiency in the disk chain can be rationalized as follows.For vanishing θ , the spatial heterogeneity of the local frequencies ω 0 n produces spatially localized normal modes [45], which are almost entirely damped by the internal dissipation of the chain, with no net transport to the output region.On the other hand, in the limit θ → ∞, the norm current J Im(ψ * n ψ n+1 ) [35] is effectively suppressed by the effect of the local dephasing noise irrespective of the disks dynamics.Optimal transport is, therefore, attained for intermediate values of the noise amplitude which provide a partial suppression of the localization mechanism [27,46].

Conclusions
In conclusion, we demonstrated a new phenomenon, dephasing-assisted spin transport, which shows how environmental noise can contribute to increasing the transfer of excitations in out-of-equilibrium spin systems.
We argued that the effect is due to the subtle effect of noise on the dynamics of spin-wave modes, as evidenced by the spectral analysis via the wavelet transform.The optimal transmission occurs at noise level large enough to excite more extended modes without degrading their overall coherence.This phenomenology is analogous to environment-assisted transport enhancement observed in open quantum systems [22,29] and originates from the emergence of an effective macroscopic Spin Wave amplitude which evolves according to a dissipative DNLS equation.On the one hand, this points towards new energy-efficient spintronics and magnonics devices.On the other hand, the generality of the DNLS model strongly suggests that dephasing-assisted transport is a general phenomenon that can be observed in a large class of classical oscillating systems.
Although our study has focused entirely on a transient nonequilibrium dynamics, we expect that dephasing-assisted transport can show up also in stationary conditions, i.e., when the disk chain steadily exchanges energy and magnetization between two external reservoirs [35].In the latter setup, which appears to be more accessible to experimental tests, stationary currents should depend non-monotonously on the dephasing noise strength.
Admittedly, our results are mostly phenomenological so far .A more detailed theoretical study is needed and will be the subject of future work.

Figure 1 .
Figure 1.Sketch of the system, consisting of 10 Py nano-disks coupled through the magneto-dipolar interaction.The magnetization in the first disk (input) is tilted away from equilibrium.The transport efficiency is the integrated Spin Wave (SW) power in the last disk (output).

Figure 3 .
Figure 3.Total SW power P vs time for different values of the dephasing noise amplitude θ (a) and of the bath temperature T (b).One can see that in the first case P drops to zero and the magnetization aligns with the z axis, since the dephasing conserves the total power.In the second case, the bath temperature excites the dynamics and the system thermalizes with P increasing with the bath temperature T.

Figure 4 .
Figure 4. Effect of the dephasing noise on the dynamics of the macrospin chain: time evolution of the local SW power p 10 of the last disk for different values of θ.Transmitted power is maximized for an optimal value around θ ≈ 4. Simulation parameters as given in the text.

Figure 5 .
Figure 5. (a) Efficiency E and (b) Kuramoto parameter K versus θ.E increases of a factor 3 until θ = 6 and then decreases again, showing that transport can be effectively promoted by dephasing.On the other hand, K decreases monotonically with θ.Thus, in the present case, transport is not related to phase synchronization.

Figure 6 .
Figure 6.Wavelet analysis of the complex spin amplitudes ψ n on the central disk (n = 5, left panels) and on the last disk (n = 10, right panels) for different values of θ.Each plot shows the density map of the average square modulus |G n (ω, t)| 2 of the Gabor transform (7) averaged over a sample of 32 independent realizations of the dyanamics.The parameter a has been set equal to 7.5 ns 2 , optimized so as to maximize the resolution in both the time and frequency domains.Notice the difference in the density scales.