Reservoir Computing is a type of recursive neural network commonly used for recognizing and predicting spatio-temporal events relying on a complex hierarchy of nested feedback loops to generate a memory functionality. The Reservoir Computing paradigm does not require any knowledge of the reservoir topology or node weights for training purposes and can therefore utilize naturally existing networks formed by a wide variety of physical processes. Most efforts to implement reservoir computing prior to this have focused on utilizing memristor techniques to implement recursive neural networks. This paper examines the potential of magnetic skyrmion fabrics and the complex current patterns which form in them as an attractive physical instantiation for Reservoir Computing. We argue that their nonlinear dynamical interplay resulting from anisotropic magnetoresistance and spin-torque effects allows for an effective and energy efficient nonlinear processing of spatial temporal events with the aim of event recognition and prediction.

## I. INTRODUCTION

A great deal has been written about the end of CMOS scaling, continuation of Moore’s Law and the need for alternative models of computing and related technologies. One of the most authoritative discussions on Moore’s Law can be found in the “final” International Technology Roadmap for Semiconductors (ITRS)^{1} published in 2015 and which had been continuously published since 1991. It predicted that CMOS transistors would quit shrinking in 2021 with the 5 nm node and that a great many technical challenges would need to be met for the 5 nm node to be economically viable. Those technical challenges were primarily related to controlling the economic costs associated with lithography, packaging, testing and the process technology itself. The technical challenges of gate leakage, interconnect power losses and material integration were all considered daunting and drove the economic issues to the point where further scaling was possible but not economically viable. The inescapable conclusion of the final ITRS document is that the scientific research community in collaboration with industry must investigate alternative models of computing including those which appear to be quite radical.

The universe of alternative models of computing is enormous and growing. A nonexclusive list includes^{2} membrane computing, DNA computing, immune computing, quantum computing, neuromorphic computing, in-materio computing, swarm computing, analog computing, chaos/edge-of-chaos computing, computational aspects of dynamics of complex systems, self-organizing systems (e.g. multiagent systems, cellular automata, artificial life), and many others. In evaluating this list of alternative computational paradigms, one is reminded of Kroemer’s Law^{3} which states that “the principal applications of any new and innovative technology always have been and will continue to be created by that new technology”. One should therefore not only judge new technologies by how they fit in with present applications but specifically by their potential to create innovative applications themselves. Optimally, however, they will satisfy both criteria.

Predicted^{4} and subsequently discovered^{5} over the past two decades, skyrmions are widely regarded as promising candidates for spintronic applications due to their room-temperature stability^{6–11} and mobility at ultra-low current densities.^{12,13} Based on these alone, skrymions have been proposed for enhancing a spectrum of existing technologies such as racetrack memories,^{14–17} transistors^{18} and logic gates.^{19,20} However, the integration of their intrinsic two-dimensional nature also enables their use for radically new technologies which have not been addressed until now.^{21–25}

In this paper we focus on Reservoir Computing (RC) models^{26–31} implemented with self-organizing neural networks in complex magnetic textures.^{22} The nodes are represented by magnetic skyrmions and the random connectivity by low magnetoresistive pathways in the material. Specifically, we will consider the effects of anisotropic magnetoresistance (AMR) on the conductivity pathways in systems with broken bulk and surface inversion symmetry. For purposes of placement in the taxonomy above, RC models are one category of neuromorphic computing and utilizing complex magnetic systems is an example of in-materio computing.^{32}

The paper is organised as follows. Sec. II contains a selective review of RC focusing on the requirements placed on the physical implementation system. In Sec. III we first briefly review the skyrmion and micromagnetic literature relevant to RC implementation and then present new simulation results to address the suitability of magnetic skyrmion networks for RC. In particular, we address the complex current patterns which form in skyrmion fabrics resulting from the interplay of anisotropic magnetoresistance and spin-torque effects and comment on optimally tuning the applied voltages to maximize the fabric’s susceptibility without actually displacing its magnetic structures. In Sec. IV we conclude how well the capabilities of magnetic substrates meet the needs of RC and delineate areas of future research.

## II. RESERVOIR COMPUTING

A Recursive Neural Network (RNN) is a network of nonlinear processing units (similar to neurons in the brain) with weighted connections between them (synapses) characterized by a flow of information that feeds back in loops. These loops imbue the system with a memory functionality derived from the “echo” of previous input sequences persisting over time. Previously, RNN’s have been trained to perform specific tasks by algorithmically tuning the weights of each connection in the network. However, the presence of feedback loops makes RNN extremely difficult to train because of bifurcation points and the tendency to exhibit chaotic solutions. As such, the same property responsible for their power is also responsible for their limited applications to date.

Reservoir Computing (RC) models address this problem by treating the recurrent part (the reservoir) differently than both the read-in and read-outs from it.^{27} This completely eliminates the need to train the large complex reservoir and only train the output weights. Because the latter are taken to not have any loops, they can be trained by straightforward linear regression techniques, i.e. least square algorithms. The usefulness of this framework derives from the reservoir’s capacity to project different spatial-temporal events into a sparsely populated high dimensional space where they become easier to recognize and categorize. The ultimate separation of such events is a complex function of the network topology, synaptic delays and the variation in the node response functions including but not limited to their nonlinear properties. This has been shown with mathematical rigor,^{29–31} allowing RC to be successfully demonstrated for spatial temporal event recognition and prediction.

Unlike designed computation, where each device has a specific role, computation in RC networks does not rely on specific devices in specific roles, but is encoded in the collective nonlinear dynamics excited by an applied input signal. This most attractive aspect of RC results from the fact it is unnecessary to have any knowledge related to the structure of the reservoir itself. It is not necessary to know the reservoir’s structure, individual node connections, weights, nor any of their nonlinear characteristics. For this reason RC methods may use “found” networks that may be self-assembled autonomously by a wide variety of physical processes.

Reservoir Computing is typically further divided into two categories, Echo State Networks (ESNs) and Liquid State Networks (LSNs).^{29} LSNs are characterized by nodes whose state is considered a continuous-time binary value (on/off) whose spiking frequency is determined by the activity of neighboring nodes. Even though they are not actively considered in this work, we note in passing that they are considered more reminiscent of bio-inspired neuromorphic models since they directly emulate the spiking potential generation and propagation observed in biological systems. On the other hand, ESNs are characterized by nodes whose state is defined by continuous values to be updated in discrete time steps depending on the state of nearby nodes. The ESN and LSN approaches alike both provide a similar function – they project the input spatio-temporal events into a high-dimensional space defined by the *state* of the entire reservoir.

### A. Echo state networks

A generic ESN is shown schematically in Fig. 1. It consists of an input layer, the reservoir network, and the output layer. The individual nodes in the reservoir network are indicated by light blue dots and the nodes of the input (output) layer are represented by light green (red) dots. The arrows indicate connection paths. The input and output nodes represent a feedforward network (black arrows) whereas the reservoir nodes are in general bidirectional forming both single node and multimode loops (pink arrows). The combination of feedforward and feedback connections results in recursive operation.

Following standard nomenclature,^{33} the input signal ** u**(

*n*) is a discrete function of time,

**(**

*x**n*) represents the state of the reservoir, e.g. the electric potential at every node, and

**(**

*y**n*) represents the output signal. The corresponding weighting functions are denoted as $W^in$, $W^$ and $W^out$ for the input, the reservoir, and the output. In general, the weights of the reservoir and the input, i.e. $Win^$ and $W^$ are time-independent. As already stated, the reservoir weighting function $W^$ can be unknown. This allows random networks of self assembled reservoirs produced by a wide variety of physical processes to be used (e.g. networks formed by the skyrmion fabrics discussed in this paper).

The evolution of the reservoir state at step *n* in reaction to its input ** u**(

*n*) and the reservoir state at time step

*n*− 1 is defined as:

where $x\u0303$ is a provisional state variable, sig(⋅) is generically a sigmoidal function biasing the input signal so that the reservoir is excited but not saturated, and λ is a *leakage* parameter characterizing the lossiness of the network’s memory. Note, that for systems without leakage (λ = 0), one has $x(n)=x\u0303(n)$. In a natural reservoir system, the specific form of the sigmoidal function is determined by the system’s physical properties.

The ultimate goal of the ESN is to classify similar input signals into identical outputs. This is achieved by first training the network on a sample set of pre-classified inputs *u*^{train}(*n*). Denoting by *y*^{train}(*n*) and *y*^{target}(*n*) the system’s response to the training inputs and the desired target output respectively, we can calculate the error *E* between them, averaged over all *N*_{y} output nodes and *T* time steps:

The training operation is defined by finding the one-dimensional scalar array of output weights $W^out$ that minimizes the error function *E*. While the specific algorithms implementing the minimization task are beyond the scope of this paper, we would like to emphasize again that Eq. (4) only requires the output weights to be modified. This is inherently more efficient and robust than the methods required for training full RNNs. Furthermore, since training of the output weights does not modify the reservoir in any way, different features of the reservoir can be searched for simultaneously by setting multiple output arrays in parallel. This makes RC well suited for sensor fusion type applications.^{34,35} It should be noted that since the training involves solving a straightforward linear regression equation, the overhead associated with training each feature can be very small relative to the amount of processing that can be done with the system.

To summarize, the reservoir requires a few qualitative key characteristics to properly function in an Echo State computing system:

It must have a short term memory, i.e. be recursively connected and/or use nodes with internal memory. As discussed, this guarantees the reservoir’s sensitivity to the input’s temporal correlations.

The dimensionality of the reservoir’s state space must be much larger than the input array. This corresponds basically to the number of nodes in the reservoir. The larger the reservoir, the greater the separation and probability of the linear classifier being able to successfully recognise specific events.

The response of the reservoir must be a nonlinear function of its inputs and previous states. The stronger the nonlinearity effects the larger the seperation of events in the reservoir’s phase space facilitating classification during training.

The reservoir’s echo state time, defined as the timescale beyond which the reservoir dynamics effectively lose all initialization information, must be much larger than the largest relevant temporal correlations in the input signals. The reservoir’s echo state time can be tuned by varying the leakage parameter λ.

These reservoir properties (and the parameters controlling them) determine the performance of the recursive network for spatial temporal event recognition and prediction. For computational materials like the magnetic textures considered in this paper, these (or an equivalent set) must be ultimately deduced by experimental measurements and tuned for optimal performance by leveraging physical insight over the materials being used. In the following section, we will proceed by reviewing the literature of magnetic skyrmions relevant to a potential RC implementation.

## III. MAGNETIC SKYRMIONS AND “SKYRMION FABRICS”

Magnetic skyrmions are nontrivial topological magnetic textures that were predicted more than two decades ago.^{4} Experimentally, they were first discovered in the form of a skyrmion lattice in 2009^{5} and later also as isolated magnetic textures.^{36} Their presence has been observed in many device-relevant materials and their properties have been extensively summarized in several reviews.^{37–40} Skyrmions are regarded as promising candidates for spintronic applications due to their mobility when driven by ultra-low currents^{12,13} and their room temperature stability.^{6–11} Particularly, the skyrmion racetrack memory has been a significant driver for intensively studying individual skyrmions.^{14–17}

Not much attention has, however, been given to applications involving intermediate skyrmion phases known as “skyrmion fabrics”.^{22} These are phases that interpolate between single skyrmions, skyrmion crystals and magnetic domain walls^{41} (examples shown in the top panels of Fig. 2). In the past, skyrmion fabrics have been studied only to observe how the different phases contribute to transitions between them.^{42}

We claim that skyrmion fabrics can provide a good basis for RC reservoirs in light of their random phase structure. The input signals can be realized via voltage patterns applied directly to the magnetic texture through various nanocontacts. Magnetoresistive effects^{43–45} such as the anisotropic magnetoresistance (AMR) will then guarantee that a certain magnetic texture will result in a unique corresponding current pattern throughout the reservoir (shown in middle panels of Fig. 2). As the input signal is varied, the magnetic texture will deform due to the interplay of AMR-mediated spin-torques and local pinning. These deformations will take place on a timescale that is slower than that taken by the current density pattern to adjust to the magnetic texture. The delay between the two will guarantee both the nonlinear resistance response of the magnetic material as well as the echoing of the spatial-temporal correlations in the dynamical response of the magnetic texture. Note that, for this to work properly, the skyrmion fabric must relax back as close as possible to its initial state upon switching off the voltage inputs. This in turn will guarantee that the reservoir operates reproducibly over separate runs. As such, optimal operation of the reservoir will require to tune voltage intensities and input signal frequencies with the objective to maximize the fabric’s susceptibility without actually displacing its constituting skyrmions and domain walls. Overall, the random skyrmion structure and corresponding current pattern will effectively model the reservoir node and weight structure discussed in the previous section.

The magnetization profiles and the current paths shown in Fig. 2 have been obtained by micromagnetic simulations using Micromagnum^{46} and selfwritten software extensions where the magnetization dynamics and the current paths have been computed selfconsistently.^{22} The magnetization dynamics are given by the Landau-Lifshitz-Gilbert (LLG) equation for the magnetization direction ** m** =

**/**

*M**M*

_{s}with spin-transfer-torque effects:

^{47,48}

Here, *M*_{s} is the saturation magnetization, the factor *ξ* = *Pμ*_{B}/(*eM*_{s}) contains the polarization *P* the electron charge *e* and the Bohr magneton *μ*_{B}. The effective magnetic field is given by $Beff=\u2212Ms\u22121(\delta F[m]/\delta m)$, where the micromagnetic free energy comprising exchange, anisotropy and dipolar interactions is:

and $FDMIB=\u222bDBm\u22c5(\u2207\xd7m)\u2009dV$ describes Bloch and $FDMIN=\u222bDNm\u22c5[(z\xd7\u2207)\xd7m]dV$ Néel DMI.^{49–51} Note that the current density *j* is a function of the applied voltage and the local magnetization. Since current density relaxation happens on a much faster time scale than the magnetization dynamics, it can be calculated self-consistently based on the AMR effect^{52} through ** j**[

*U*,

**] = −**

*m***[**

*σ***] ⋅**

*m***[**

*E**U*]. Here the electric field induced by the applied voltage is calculated by solving the Poisson equation

**= −∇Φ with boundary conditions Φ|**

*E*_{c1}= −Φ|

_{c2}=

*U*at the two contacts, and the conductivity tensor $\sigma [m]=1\rho \u22a51+(1\rho \u2225\u22121\rho \u22a5)m\u2297m$ varies with the local magnetization. We denote by

*ρ*

_{⊥}(

*ρ*

_{∥}) the current resistivities for flows perpendicular (parallel) to the magnetization direction. Based on these we define

*σ*

_{0}= (1/

*ρ*

_{∥}+ 2/

*ρ*

_{⊥})/3, and the AMR ratio $a=2(\rho \u2225\u2212\rho \u22a5)\rho \u2225+\rho \u22a5$.

^{22,52}

Previous work has focused on details of current paths through magnetic nano-ribbons in the presence of single magnetic Bloch and Néel skrymions.^{22} As a result of the AMR effect, Néel skyrmions show a tendency to deflect current flow lines tangentially around their centers while Bloch skyrmions favor current flows through their centers. Furthermore, single skyrmions were shown to exhibit non-linear current-voltage characteristics due to the interplay of magnetoresistive effects and pinning.

This work expands on previous results by exploring the effect of Néel and Bloch skyrmion fabrics on their resulting AMR-mediated current flow. In Fig. 2, we show examples of Bloch (left column) and Néel (right column) skyrmion fabrics. Their complex magnetic texture (top row) reflects on the total current flow between two contacts where a voltage difference is applied (middle row). To isolate the pure AMR effect of the skyrmion fabric we subtract the current response of a trivial out-of-plane ferromagnetic state and plot the resulting differential current flow (bottom row). We find that this artificially constructed AMR-dominated regime exhibits current backflows reminiscent of the recursive connectivity required in a RC network.

Due to the tunable properties of magnetic skyrmions, the skyrmion fabric can be tweaked by application of static magnetic fields. These alter the density of skyrmions, thus tuning the effective node density throughout the sample. Figure 3 demonstrates this by showing how a Bloch magnetic texture and its resulting differential current flows both respond to variations of out-of-plane applied magnetic fields.

## IV. CONCLUSIONS

We have argued that skyrmion fabrics embedded in broken inversion symmetry magnetic substrates are potentially attractive options for implementing Echo State (ES) recognition and prediction systems. The principle result from this paper shows how skyrmion fabrics induce a strongly perturbed current flow through the magnetic texture as compared to the expected current flow through the homogenously magnetized state. We further show that such a perturbed current flow creates conduction paths similar to those required for Echo State Networks. Leveraging results from previous work on nonlinear I-V characteristics of individual skyrmions,^{22} we argue that skyrmion fabrics appear to meet the basic requirements of an Echo State Reservoir computing system. It should be pointed out that the functionality of the proposed system is not directly related to the magnitude of the nonlinear effect. Nonetheless, the total intensity of the AMR effect scales with the net magnetic distortion present in the magnetic sample. As such, complex random fabrics will allow for a more robust operating signal than the single isolated skyrmions previously studied.^{22} We would like to emphasize that the information processing capabilities inherent to such magnetic materials derives from the dynamic interplay of current flow with the complex magnetic texture.^{53}

The present calculation does not address the influence of thermal noise or potential grain structures in the magnetic material. Whereas our study is based solely on the anisotropic magneto-resistive (AMR) effect, other similar magnetization-modulated resistance effects could be included to tune and enhance our results. Furthermore, a properly functioning reservoir must guarantee a basis state across operational cycles. This implies that the skyrmion fabric must restore to its initial state whenever the driving voltages are switched off. Due to the dynamical character of magnetic skyrmions and domain walls, this may not be guaranteed if applied voltage intensities go beyond the threshold voltage for motion. Future work will focus on exploring how material inhomogeneities and impurities can be leveraged to tune for this threshold voltage deterministically.

The size scales of skyrmion fabrics are orders of magnitude smaller than other proposed implementations of reservoirs for ES networks, like memristor or optical networks. Furthermore, the characteristics of the skyrmion fabric (skyrmion density/size, domain wall width, etc ...) can be altered by tuning the material properties and/or applying magnetic fields. These effectively tune both the dimensionality and net nonlinearity of the reservoir’s dynamics, thus determining its performance.

## ACKNOWLEDGMENTS

We acknowledge discussions with Kai Litzius, Diana Prychynenko and Jairo Sinova. G. B. appreciates support and useful discussions with Narayan Srinivasa of Intel. We acknowledge the funding from the German Research Foundation (DFG) under the Project No. EV 196/2-1, the Alexander von Humboldt Foundation and the Transregional Collaborative Research Center (SFB/TRR) 173 Spin+X.

## REFERENCES

*23*(