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.

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.

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.

Spin Torque Efficiency and Analytic Error Rate Estimates of Skyrmion Racetrack Memory.

In this paper, the thermal stability of skyrmion bubbles and the critical currents to move them over pinning sites were investigated. For the used pinning geometries and the used parameters, the unexpected behavior is reported that the energy barrier to overcome the pinning site is larger than the energy barrier of the annihilation of a skyrmion. The annihilation takes place at boundaries by current driven motion, as well as due to the excitation over energy barriers, in the absence of currents, without forming Bloch points. It is reported that the pinning sites, which are required to allow thermally stable bits, significantly increase the critical current densities to move the bits in skyrmion-like structures to about jcrit = 0.62 TA/m². The simulation shows that the applied spin transfer model predicts experimentally obtained critical currents to move stable skyrmions at room temperature well, which is in contrast to simulations based on spin orbit torque that predict significantly too low critical currents. By calculating the thermal stability, as well as the critical current, we can derive the spin torque efficiency η = ΔE/Ic = 0.19 kBT300/µA, which is in a similar range to the simulated spin torque efficiency of MRAM structures. Finally, it is shown that the stochastic depinning process of any racetrack-like device requires an extremely narrow depinning time distribution smaller than ~6% of the current pulse length to reach bit error rates smaller than 10-9.

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.

Combining micromagnetism and magnetostatic Maxwell equations for multiscale magnetic simulations.

Magnetostatic Maxwell equations and the Landau-Lifshitz-Gilbert (LLG) equation are combined to a multiscale method, which allows to extend the problem size of traditional micromagnetic simulations. By means of magnetostatic Maxwell equations macroscopic regions can be handled in an averaged and stationary sense, whereas the LLG allows to accurately describe domain formation as well as magnetization dynamics in some microscopic subregions. The two regions are coupled by means of their strayfield and the combined system is solved by an optimized time integration scheme.

Removal of earth's magnetic field effect on magnetoelastic resonance sensors by an antisymmetric bias field.

Magnetoelastic sensors are used in a wide field of wireless sensing applications. The sensing element is a low-cost magnetostrictive ribbon whose resonant frequency depends on the measured quantity. The accuracy of magnetoelastic sensors is limited by the fact that the resonant frequency is also affected by the earth's magnetic field. In this paper we present a technique to minimize this effect by applying an antisymmetric magnetic bias field to the ribbon. The ribbon's response to external perturbation fields was measured and compared to a conventional sensor design. Our results show that the influence of the earth's magnetic field could be reduced by 77%.

A nano-positioning system for macromolecular structural analysis.

Very often, the positions of flexible domains within macromolecules as well as within macromolecular complexes cannot be determined by standard structural biology methods. To overcome this problem, we developed a method that uses probabilistic data analysis to combine single-molecule measurements with X-ray crystallography data. The method determines not only the most likely position of a fluorescent dye molecule attached to the domain but also the complete three-dimensional probability distribution depicting the experimental uncertainty. With this approach, single-pair fluorescence resonance energy transfer measurements can now be used as a quantitative tool for investigating the position and dynamics of flexible domains within macromolecular complexes. We applied this method to find the position of the 5' end of the nascent RNA exiting transcription elongation complexes of yeast (Saccharomyces cerevisiae) RNA polymerase II and studied the influence of transcription factor IIB on the position of the RNA.

Single-molecule tracking of mRNA exiting from RNA polymerase II.

Single-pair fluorescence resonance energy transfer was used to track RNA exiting from RNA polymerase II (Pol II) in elongation complexes. Measuring the distance between the RNA 5' end and three known locations within the elongation complex allows us determine its position by means of triangulation. RNA leaves the polymerase active center cleft via the previously proposed exit tunnel and then disengages from the enzyme surface. When the RNA reaches lengths of 26 and 29 nt, its 5' end associates with Pol II at the base of the dock domain. Because the initiation factor TFIIB binds to the dock domain and exit tunnel, exiting RNA may prevent TFIIB reassociation during elongation. RNA further extends toward the linker connecting to the polymerase C-terminal repeat domain (CTD), which binds the 5'-capping enzyme and other RNA processing factors.