Search

Magnonics is a field of science that addresses the physical properties of spin waves and utilizes them for data processing. Scalability down to atomic dimensions, operations in the GHz-to-THz frequency range, utilization of nonlinear and nonreciprocal phenomena, and compatibility with CMOS are just a few of many advantages offered by magnons. Although magnonics is still primarily positioned in the academic domain, the scientific and technological challenges of the field are being extensively investigated, and many proof-of-concept prototypes have already been realized in laboratories. This roadmap is a product of the collective work of many authors that covers versatile spin-wave computing approaches, conceptual building blocks, and underlying physical phenomena. In particular, the roadmap discusses the computation operations with Boolean digital data, unconventional approaches like neuromorphic computing, and the progress towards magnon-based quantum computing. The article is organized as a collection of sub-sections grouped into seven large thematic sections. Each sub-section is prepared by one or a group of authors and concludes with a brief description of the current challenges and the outlook of the further development of the research directions., Comment: 74 pages, 57 figures, 500 references

The fundamental limits currently faced by traditional computing devices necessitate the exploration of new ways to store, compute and transmit information. Here, we propose a three-dimensional (3D) magnetic interconnector that exploits geometry-driven automotion of domain walls (DWs), for the transfer of magnetic information between functional magnetic planes. By combining state-of-the-art 3D nanoprinting and standard physical vapor deposition, we prototype 3D helical DW conduits. We observe the automotion of DWs by imaging their magnetic state under different field sequences using X-ray microscopy, observing a robust unidirectional motion of DWs from the bottom to the top of the spirals. From experiments and micromagnetic simulations, we determine that the large thickness gradients present in the structure are the main mechanism for 3D DW automotion. We obtain direct evidence of how this tailorable magnetic energy gradient is imprinted in the devices, and how it competes with pinning effects due to local changes in the energy landscape. Our work also predicts how this effect could lead to high DW velocities, reaching the Walker limit during automotion. This work provides new possibilities for efficient transfer of magnetic information in three dimensions.

The Landau-Lifshitz equation governing magnetization dynamics is written in terms of the amplitudes of normal modes associated with the micromagnetic system's appropriate ground state. This results in a system of nonlinear ordinary differential equations (ODEs), the right-hand side of which can be expressed as the sum of a linear term and nonlinear terms with increasing order of nonlinearity (quadratic, cubic, etc.). The application of the method to nanostructured magnetic systems demonstrates that the accurate description of magnetization dynamics requires a limited number of normal modes, which results in a considerable improvement in computational speed. The proposed method can be used to obtain a reduced-order dynamical description of magnetic nanostructures which allows to adjust the accuracy between low-dimensional models, such as macrospin, and micromagnetic models with full spatial discretization. This new paradigm for micromagnetic simulations is tested for three problems relevant to the areas of spintronics and magnonics: directional spin-wave coupling in magnonic waveguides, high power ferromagnetic resonance in a magnetic nanodot, and injection-locking in spin-torque nano-oscillators. The case studies considered demonstrate the validity of the proposed approach to systematically obtain an intermediate order dynamical model based on normal modes for the analysis of magnetic nanosystems. The time-consuming calculation of the normal modes has to be done only one time for the system. These modes can be used to optimize and predict the system response for all possible time-varying external excitations (magnetic fields, spin currents). This is of utmost importance for applications where fast and accurate system simulations are required, such as in electronic circuits including magnetic devices.

We investigate the effect of curvature on the energy and stability of domain wall configurations in curved cylindrical nanotubes and nanowires. We use micromagnetic simulations to calculate the phase diagram for the transverse wall (TW) and vortex wall (VW) states in tubes, finding the ground state configuration and the metastability region where both types of walls can exist. The introduction of curvature shifts the range for which the TW is the ground state domain wall to higher diameters, and increases the range of metastability. We interpret this behavior to be primarily due to the curvature-induced effective Dzyaloshinskii-Moriya term in the exchange energy. Furthermore, we demonstrate qualitatively the same behavior in solid cylindrical nanowires. Comparing both tubes and wires, we observe how while in tubes curvature tends to suppress the transformation from the TW to VW, in wires it promotes the transformation of the VW containing the Bloch point into the TW. These findings have important implications in the fundamental understanding of domain walls in 3D geometries, and the design of future domain wall devices.

We perform a direct comparison between Kramers' method in many dimensions -- i.e., Langer's theory -- adapted to magnetic spin systems, and a path sampling method in the form of forward flux sampling, as a means to compute collapse rates of metastable magnetic skyrmions. We show that a good agreement is obtained between the two methods. We report variations of the attempt frequency associated with skyrmion collapse by three to four orders of magnitude when varying the applied magnetic field by 5$\%$ of the exchange strength, which confirms the existence of a strong entropic contribution to the lifetime of skyrmions. This demonstrates that in complex systems, the knowledge of the rate prefactor, in addition to the internal energy barrier, is essential in order to properly estimate a lifetime., Comment: 5 pages, 5 figures (main text), 8 pages including supplemental material

A method to create a highly homogeneous magnetic field by applying topology optimized, additively manufactured shimming elements is investigated. The topology optimization algorithm can calculate a suitable permanent and nonlinear soft magnetic design that fulfills the desired field properties. The permanent magnetic particles are bonded in a polyamide matrix, and they are manufactured with a low-cost, end-user 3D printer. Stray field measurements and an inverse stray field simulation framework can determine printing and magnetization errors. The customized shimming elements are manufactured by a selective melting process which produces completely dense soft magnetic metal parts. The methodology is demonstrated on an example of two axial symmetric cylindrical magnets. In this case, the homogeneity can be increased by a factor of 35. Simulation and measurement results point out a good conformity.

We compute annihilation rates of metastable magnetic skyrmions using a form of Langer's theory in the intermediate-to-high damping (IHD) regime. For a N\'eel skyrmion, a Bloch skyrmion, and an antiskyrmion, we look at two possible paths to annihilation: collapse and escape through a boundary. We also study the effects of a curved vs. a flat boundary, a second skyrmion and a non-magnetic defect. We find that the skyrmion's internal modes play a dominant role in the thermally activated transitions compared to the spin-wave excitations and that the relative contribution of internal modes depends on the nature of the transition process. Our calculations for a small skyrmion stabilized at zero-field show that collapse on a defect is the most probable path. In the absence of a defect, the annihilation is largely dominated by escape mechanisms, even though in this case the activation energy is higher than that of collapse processes. Escape through a flat boundary is found more probable than through a curved boundary. The potential source of stability of metastable skyrmions is therefore found not to lie in high activation energies, nor in the dynamics at the transition state, but comes from entropic narrowing in the saddle point region which leads to lowered attempt frequencies. This narrowing effect is found to be primarily associated with the skyrmion's internal modes., Comment: 14 pages, 9 figures

Permanent magnet materials require a high uniaxial magneto-crystalline anisotropy. Exchange coupling between small crystallites with easy-plane anisotropy induces an effective uniaxial anisotropy if arranged accordingly. Nanostructuring of materials with easy-plane anisotropy is an alternative way to create hard-magnetic materials. The coercivity increases with decreasing feature size. The resulting coercive field is about 12 percent of the anisotropy field for a crystal size of 3.4 times the Bloch parameter.

Conjugate gradient methods for energy minimization in micromagnetics are compared. When the step length in the line search is controlled, conjugate gradient techniques are a fast and reliable way to compute the hysteresis properties of permanent magnets. The method is applied to investigate demagnetizing effects in NdFe12 based permanent magnets. The reduction of the coercive field by demagnetizing effects is 1.4 T at 450 K.

We investigate the correlation between roughness, remanence and coercivity in Co/Ni films grown on Cu seed layers of varying thickness. Increasing the Cu seed layer thickness of Ta/Cu/8x[Co/Ni] thin films increases the roughness of the films. In-plane magnetization loops show that both the remanance and coercivity increase with increasing seed layer roughness. Polar Kerr microscopy and magnetic force microscopy reveal that the domain density also increases with roughness. Finite element micromagnetic simulations performed on structures with periodically modulated surfaces provide further insight. They confirm the connection between domain density and roughness, and identify the microsocpic structure of the domain walls as the source of the increased remanence in rough films. The simulations predict that the character of the domain walls changes from Bloch-like in smooth films to N\'eel-like for rougher films.

We derive a hierarchy of stochastic evolution equations for pure states (quantum trajectories) to efficiently solve open quantum system dynamics with non-Markovian structured environments. From this hierarchy of pure states (HOPS) the exact reduced density operator is obtained as an ensemble average. We demonstrate the power of HOPS by applying it to the Spin-Boson model, the calculation of absorption spectra of molecular aggregates and energy transfer in a photosynthetic pigment-protein complex.

We report a finite temperature micromagnetic method (FTM) that allows for the calculation of the coercive field of arbitrary shaped magnetic nanostructures at time scales of nanoseconds to years. Instead of directly solving the Landau-Lifshitz-Gilbert equation, the coercive field is obtained without any free parameter by solving a non linear equation, which arises from the transition state theory. The method is applicable to magnetic structures where coercivity is determined by one thermally activated reversal or nucleation process. The method shows excellent agreement with experimentally obtained coercive fields of magnetic nanostructures and provides a deeper understanding of the mechanism of coercivity., Comment: submitted to Phys. Rev. B

A numerical implementation of the transition state theory (TST) is presented which can be used to calculate the attempt frequency $f_{0}$ of arbitrary shaped magnetic nanostructures. The micromagnetic equations are discretized using the finite element method. The climbing image nudged elastic band method is used to calculate the saddle point configuration, which is required for the calculation of $f_{0}$. Excellent agreement of the implemented numerical model and analytical solutions is obtained for single domain particles. The developed method is applied to compare $f_{0}$ for single phase and graded media grains of advanced recording media. $f_{0}$ is predicted to be comparable if the maximum anisotropy is the same in these two media types.

The potential of exchange spring bilayers and graded media is reviewed. An analytical model for the optimization of graded media gives an optimal value of the magnetic polarization of Js = 0.8 T. The optimum design allows for thermally stable grains with grain diameters in the order of 3.3 nm, which supports ultra high density up to 5 to 10 Tbit per inch2. The switching field distribution is significantly reduced in bilayer media and graded media compared to single phase media. For the graded media the switching field distribution is reduced by about a factor of two. For bilayer media the minimum switching field distribution is obtained for soft layer anisotropies about one fifth of the hard layer anisotropy. The influence of precessional switching on the reversal time and the reversal field is investigated in detail for magnetic bilayers. Exchange spring bilayers can be reversed with field pulses of 20 ps., Comment: submitted to JMMM, 'Current Perspectives; Perpendicular recording'

We perform a qualitative analysis of phase locking in a double point-contact spinvalve system by solving the Landau-Lifshitz-Gilbert-Slonzewski equation using a hybrid-finite-element method. We show that the phase-locking behaviour depends on the applied field angle. Starting from a low field angle, the locking-current difference between the current through contact A and B increases with increasing angle up to a maximum of 14 mA at 30 degree and it decreases thereafter until it reaches a minimum of 1 mA at 75 degree. The tunability of the phase-lock frequency with current decreases linearly with increasing out of plane angle from 45 to 21 MHz/mA., Comment: 11 pages, 3 figures, submitted to APL