Dipolar skyrmions and antiskyrmions of arbitrary topological charge at room temperature.

Magnetic skyrmions are localized, stable topological magnetic textures that can move and interact with each other like ordinary particles when an external stimulus is applied. The efficient control of the motion of spin textures using spin-polarized currents opened an opportunity for skyrmionic devices such as racetrack memory and neuromorphic or reservoir computing. The coexistence of skyrmions with high topological charge in the same system promises further possibilities for efficient technological applications. In this work, we directly observe dipolar skyrmions and antiskyrmions with arbitrary topological charge in Co/Ni multilayers at room temperature. We explore the dipolar-stabilized spin objects with topological charges of up to 10 and characterize their nucleation process, their energy dependence on the topological charge and the effect of the material parameters on their stability. Furthermore, our micromagnetic simulations demonstrate spin-transfer-induced motion of these spin objects, which is important for their potential device application.

magnum.np: a PyTorch based GPU enhanced finite difference micromagnetic simulation framework for high level development and inverse design.

magnum.np is a micromagnetic finite-difference library completely based on the tensor library PyTorch. The use of such a high level library leads to a highly maintainable and extensible code base which is the ideal candidate for the investigation of novel algorithms and modeling approaches. On the other hand magnum.np benefits from the device abstraction and optimizations of PyTorch enabling the efficient execution of micromagnetic simulations on a number of computational platforms including graphics processing units and potentially Tensor processing unit systems. We demonstrate a competitive performance to state-of-the-art micromagnetic codes such as mumax3 and show how our code enables the rapid implementation of new functionality. Furthermore, handling inverse problems becomes possible by using PyTorch's autograd feature.

Topology optimization of a heat-assisted magnetic recording write head to reduce transition curvature using a binary optimization algorithm utilizing the adjoint method.

In this work, the possibility to reduce transition curvature in heat-assisted magnetic recording, using a conventional write head design, by shaping the recording field to counteract the circular profile of the heat pulse is investigated. Topology optimization of the head tip is performed in order to create the desired cross-track field profile for increasing distances from the write head tip. For the topology optimization, the adjoint method is utilized to calculate the necessary gradients and a binary optimization scheme is proposed. The optimizations are performed considering linearized material parameters reducing the computational complexity and the results are compared to optimizations incorporating the full non-linear material behavior. The optimized field profiles are evaluated for their influence on the read-back process. To do so, switching probability phase diagrams are calculated and the curvature parameter, the signal to noise ratio and the channel bit density are extracted. The presented results show that while transition curvature can be reduced by shaping the cross-track profile of the write field, this alone does not consequently lead to an improvement of the read back process. Therefore, completely new head designs, considering additional parameters have to be investigated.

Domain Wall Automotion in Three-Dimensional Magnetic Helical Interconnectors.

The fundamental limits currently faced by traditional computing devices necessitate the exploration of ways to store, compute, and transmit information going beyond the current CMOS-based technologies. 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 that are 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 demonstrates a possible mechanism for efficient transfer of magnetic information in three dimensions.

A topology optimization algorithm for magnetic structures based on a hybrid FEM-BEM method utilizing the adjoint approach.

A method to optimize the topology of hard as well as soft magnetic structures is implemented using the density approach for topology optimization. The stray field computation is performed by a hybrid finite element-boundary element method. Utilizing the adjoint approach the gradients necessary to perform the optimization can be calculated very efficiently. We derive the gradients using a "first optimize then discretize" scheme. Within this scheme, the stray field operator is self-adjoint allowing to solve the adjoint equation by the same means as the stray field calculation. The capabilities of the method are showcased by optimizing the topology of hard as well as soft magnetic thin film structures and the results are verified by comparison with an analytical solution.

Complex free-space magnetic field textures induced by three-dimensional magnetic nanostructures.

The design of complex, competing effects in magnetic systems-be it via the introduction of nonlinear interactions1-4, or the patterning of three-dimensional geometries5,6-is an emerging route to achieve new functionalities. In particular, through the design of three-dimensional geometries and curvature, intrastructure properties such as anisotropy and chirality, both geometry-induced and intrinsic, can be directly controlled, leading to a host of new physics and functionalities, such as three-dimensional chiral spin states7, ultrafast chiral domain wall dynamics8-10 and spin textures with new spin topologies7,11. Here, we advance beyond the control of intrastructure properties in three dimensions and tailor the magnetostatic coupling of neighbouring magnetic structures, an interstructure property that allows us to generate complex textures in the magnetic stray field. For this, we harness direct write nanofabrication techniques, creating intertwined nanomagnetic cobalt double helices, where curvature, torsion, chirality and magnetic coupling are jointly exploited. By reconstructing the three-dimensional vectorial magnetic state of the double helices with soft-X-ray magnetic laminography12,13, we identify the presence of a regular array of highly coupled locked domain wall pairs in neighbouring helices. Micromagnetic simulations reveal that the magnetization configuration leads to the formation of an array of complex textures in the magnetic induction, consisting of vortices in the magnetization and antivortices in free space, which together form an effective B field cross-tie wall14. The design and creation of complex three-dimensional magnetic field nanotextures opens new possibilities for smart materials15, unconventional computing2,16, particle trapping17,18 and magnetic imaging19.

Dipolar-stabilized first and second-order antiskyrmions in ferrimagnetic multilayers.

Skyrmions and antiskyrmions are topologically protected spin structures with opposite vorticities. Particularly in coexisting phases, these two types of magnetic quasi-particles may show fascinating physics and potential for spintronic devices. While skyrmions are observed in a wide range of materials, until now antiskyrmions were exclusive to materials with D2d symmetry. In this work, we show first and second-order antiskyrmions stabilized by magnetic dipole-dipole interaction in Fe/Gd-based multilayers. We modify the magnetic properties of the multilayers by Ir insertion layers. Using Lorentz transmission electron microscopy imaging, we observe coexisting antiskyrmions, Bloch skyrmions, and type-2 bubbles and determine the range of material properties and magnetic fields where the different spin objects form and dissipate. We perform micromagnetic simulations to obtain more insight into the studied system and conclude that the reduction of saturation magnetization and uniaxial magnetic anisotropy leads to the existence of this zoo of different spin objects and that they are primarily stabilized by dipolar interaction.

Non-Planar Geometrical Effects on the Magnetoelectrical Signal in a Three-Dimensional Nanomagnetic Circuit.

Expanding nanomagnetism and spintronics into three dimensions (3D) offers great opportunities for both fundamental and technological studies. However, probing the influence of complex 3D geometries on magnetoelectrical phenomena poses important experimental and theoretical challenges. In this work, we investigate the magnetoelectrical signals of a ferromagnetic 3D nanodevice integrated into a microelectronic circuit using direct-write nanofabrication. Due to the 3D vectorial nature of both electrical current and magnetization, a complex superposition of several magnetoelectrical effects takes place. By performing electrical measurements under the application of 3D magnetic fields, in combination with macrospin simulations and finite element modeling, we disentangle the superimposed effects, finding how a 3D geometry leads to unusual angular dependences of well-known magnetotransport effects such as the anomalous Hall effect. Crucially, our analysis also reveals a strong role of the noncollinear demagnetizing fields intrinsic to 3D nanostructures, which results in an angular dependent magnon magnetoresistance contributing strongly to the total magnetoelectrical signal. These findings are key to the understanding of 3D spintronic systems and underpin further fundamental and device-based studies.

Strayfield calculation for micromagnetic simulations using true periodic boundary conditions.

We present methods for calculating the strayfield in finite element and finite difference micromagnetic simulations using true periodic boundary conditions. In contrast to pseudo periodic boundary conditions, which are widely used in micromagnetic codes, the presented methods eliminate the shape anisotropy originating from the outer boundary. This is a crucial feature when studying the influence of the microstructure on the performance of composite materials, which is demonstrated by hysteresis calculations of soft magnetic structures that are operated in a closed magnetic loop configuration. The applied differential formulation is perfectly suited for the application of true periodic boundary conditions. The finite difference equations can be solved by a highly efficient Fast Fourier Transform method.

Intrinsic DMI-free skyrmion formation and robust dynamic behaviors in magnetic hemispherical shells.

We performed finite-element micromagnetic simulations to examine the formation of skyrmions without intrinsic Dzyaloshinskii-Moriya interaction (DMI) in magnetic hemispherical shells. We found that curvature-induced DM-like interaction allows for further stabilization of skyrmions without the DMI in curved-geometry hemispherical shells for a specific range of uniaxial perpendicular magnetic anisotropy (PMA) constant Ku. The larger the curvature of the shell, the higher the Ku value required for the formation of the skyrmions. With well-stabilized skyrmions, we also found in-plane gyration modes and azimuthal spin-wave modes as well as an out-of-plane breathing mode, similarly to previously found modes for planar geometries. Furthermore, additional higher-frequency hybrid modes were observed due to coupling between the gyration and azimuthal modes. This work provides further physical insight into the static and dynamic properties of intrinsic DMI-free skyrmions formed in curved-geometry systems.

Magnetic Position System Design Method Applied to Three-Axis Joystick Motion Tracking.

This manuscript discusses the difficulties with magnetic position and orientation (MPO) system design and proposes a general method for finding optimal layouts. The formalism introduces a system quality measure through state separation and reduces the question "How to design an MPO system?" to a global optimization problem. The latter is then solved by combining differential evolution algorithms with magnet shape variation based on analytical computations of the field. The proposed formalism is then applied to study possible realizations of continuous three-axis joystick motion tracking, realized with just a single magnet and a single 3D magnetic field sensor. The computations show that this is possible when a specific design condition is fulfilled and that large state separations as high as 1mT/∘ can be achieved under realistic conditions. Finally, a comparison to state-of-the-art design methods is drawn, computation accuracy is reviewed critically, and an experimental validation is presented.

Control of the noncollinear interlayer exchange coupling.

Interlayer exchange coupling in transition metal multilayers has been intensively studied for more than three decades and is incorporated into almost all spintronic devices. With the current spacer layers, only collinear magnetic alignment can be reliably achieved; however, controlling the coupling angle has the potential to markedly expand the use of interlayer exchange coupling. Here, we show that the coupling angle between the magnetic moments of two ferromagnetic layers can be precisely controlled by inserting a specially designed magnetic metallic spacer layer between them. The coupling angle is controlled solely by the composition of the spacer layer. Moreover, the biquadratic coupling strength, responsible for noncollinear alignment, is larger than that of current materials. These properties allow for the fabrication and study of not yet realized magnetic structures that have the potential to improve existing spintronic devices.

Micromagnetic Simulations of Submicron Vortex Structures for the Detection of Superparamagnetic Labels.

We present a numerical investigation on the detection of superparamagnetic labels using a giant magnetoresistance (GMR) vortex structure. For this purpose, the Landau-Lifshitz-Gilbert equation was solved numerically applying an external z-field for the activation of the superparamagnetic label. Initially, the free layer's magnetization change due to the stray field of the label is simulated. The electric response of the GMR sensor is calculated by applying a self-consistent spin-diffusion model to the precomputed magnetization configurations. It is shown that the soft-magnetic free layer reacts on the stray field of the label by shifting the magnetic vortex orthogonally to the shift direction of the label. As a consequence, the electric potential of the GMR sensor changes significantly for label shifts parallel or antiparallel to the pinning of the fixed layer. Depending on the label size and its distance to the sensor, the GMR sensor responds, changing the electric potential from 26.6 mV to 28.3 mV.

Additive Manufactured and Topology Optimized Passive Shimming Elements for Permanent Magnetic Systems.

A method to create a highly homogeneous magnetic field by applying topology optimized, additively manufactured passive 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 a simple example of two axial symmetric cylindrical magnets, which generates a high inhomogeneous magnetic field. In this case, the maximum magnetic field density is 25 mT and the the homogeneity can be increased by a factor of 35 or down to 6‰. Simulation and measurement results point out a good conformity. Additional topology optimizations of more than one shimming element layer show the opportunity to make the manufactured magnetic system even more homogeneous.

3D Printing of Polymer-Bonded Rare-Earth Magnets With a Variable Magnetic Compound Fraction for a Predefined Stray Field.

Additive manufacturing of polymer-bonded magnets is a recently developed technique, for single-unit production, and for structures that have been impossible to manufacture previously. Also, new possibilities to create a specific stray field around the magnet are triggered. The current work presents a method to 3D print polymer-bonded magnets with a variable magnetic compound fraction distribution. This means the saturation magnetization can be adjusted during the printing process to obtain a required external field of the manufactured magnets. A low-cost, end-user 3D printer with a mixing extruder is used to mix permanent magnetic filaments with pure polyamide (PA12) filaments. The magnetic filaments are compounded, extruded, and characterized for the printing process. To deduce the quality of the manufactured magnets with a variable magnetic compound fraction, an inverse stray field framework is developed. The effectiveness of the printing process and the simulation method is shown. It can also be used to manufacture magnets that produce a predefined stray field in a given region. This opens new possibilities for magnetic sensor applications. This setup and simulation framework allows the design and manufacturing of polymer-bonded permanent magnets, which are impossible to create with conventional methods.

Solving Large-Scale Inverse Magnetostatic Problems using the Adjoint Method.

An efficient algorithm for the reconstruction of the magnetization state within magnetic components is presented. The occurring inverse magnetostatic problem is solved by means of an adjoint approach, based on the Fredkin-Koehler method for the solution of the forward problem. Due to the use of hybrid FEM-BEM coupling combined with matrix compression techniques the resulting algorithm is well suited for large-scale problems. Furthermore the reconstruction of the magnetization state within a permanent magnet as well as an optimal design application are demonstrated.

A self-consistent spin-diffusion model for micromagnetics.

We propose a three-dimensional micromagnetic model that dynamically solves the Landau-Lifshitz-Gilbert equation coupled to the full spin-diffusion equation. In contrast to previous methods, we solve for the magnetization dynamics and the electric potential in a self-consistent fashion. This treatment allows for an accurate description of magnetization dependent resistance changes. Moreover, the presented algorithm describes both spin accumulation due to smooth magnetization transitions and due to material interfaces as in multilayer structures. The model and its finite-element implementation are validated by current driven motion of a magnetic vortex structure. In a second experiment, the resistivity of a magnetic multilayer structure in dependence of the tilting angle of the magnetization in the different layers is investigated. Both examples show good agreement with reference simulations and experiments respectively.

A three-dimensional spin-diffusion model for micromagnetics.

We solve a time-dependent three-dimensional spin-diffusion model coupled to the Landau-Lifshitz-Gilbert equation numerically. The presented model is validated by comparison to two established spin-torque models: The model of Slonzewski that describes spin-torque in multi-layer structures in the presence of a fixed layer and the model of Zhang and Li that describes current driven domain-wall motion. It is shown that both models are incorporated by the spin-diffusion description, i.e., the nonlocal effects of the Slonzewski model are captured as well as the spin-accumulation due to magnetization gradients as described by the model of Zhang and Li. Moreover, the presented method is able to resolve the time dependency of the spin-accumulation.

Spin-polarized transport in ferromagnetic multilayers: An unconditionally convergent FEM integrator.

We propose and analyze a decoupled time-marching scheme for the coupling of the Landau-Lifshitz-Gilbert equation with a quasilinear diffusion equation for the spin accumulation. This model describes the interplay of magnetization and electron spin accumulation in magnetic and nonmagnetic multilayer structures. Despite the strong nonlinearity of the overall PDE system, the proposed integrator requires only the solution of two linear systems per time-step. Unconditional convergence of the integrator towards weak solutions is proved.