Source author record

Dieter Suess

Dieter Suess appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2026arXiv

Magnetization reversal mechanism of double-helix nanowires probed by dark-field magneto-optical Kerr effect

Double-helix (DH) nanowires provide a platform to study the influence of geometric chirality on spin chirality. Their three-dimensional (3D) helical architecture and tunable inter-strand coupling enable control of spin chirality, including the stabilization of topological 3D magnetic states such as helical domains and domain walls, topological stray fields, and extended helical vortex/skyrmion tubes. So far, the study of these and other 3D nanostructures is usually confined to a limited number of magnetic microscopy experiments in large facilities. Here, we investigate the reversal mechanism of a single DH nanowire using Dark-Field magneto-optical Kerr effect (DF-MOKE) magnetometry under external 3D magnetic fields. By analyzing the angular dependence of the DF-MOKE signal, we fit the reversal process using established models for domain-wall nucleation and propagation, finding a characteristic behavior similar to that reported for cylindrical nanowires. Micromagnetic simulations indicate that the reversal process goes through nucleation of the helical vortex tube in a curling manner while ptychographic X-ray magnetic circular dichroism data reveal that this helical vortex tube state forms through a mixed nucleation-propagation process. These observations provide a consistent microscopic picture of reversal mediated by a helical vortex tube extending along the nanowire. Our work provides a comprehensive characterization of magnetization reversal in DH nanowires and demonstrates that DF-MOKE magnetometry is effective for probing reversal mechanisms in single 3D nanostructures. This lab-based approach expands the range of accessible experiments beyond large-scale facilities, enabling extensive exploration of the rich spin states supported by 3D nano-geometries.

preprint2022arXiv

Generation and annihilation of skyrmions and antiskyrmions in magnetic heterostructures

We demonstrate the controlled generation and annihilation of (anti)skyrmions with tunable chirality in magnetic heterostructures by means of micromagnetic simulations. By making use of magnetic (anti)vortices in patterned ferromagnetic layer, we stabilize full lattices of (anti)skyrmions in an underlying skyrmionic thin film in a reproducible manner. The stability of the (anti)skyrmion depends on the polarization of the (anti)vortex, whereas their chirality is given by those of the (anti)vortices. Furthermore, we demonstrate that the core coupling between the (anti)vortices and (anti)skyrmions allows to annihilate the spin-objects in a controlled fashion by applying short pulses of in-plane external magnetic fields, representing a new key paradigm in skyrmionic devices.

preprint2022arXiv

In-situ alignment of anisotropic hard magnets of 3D printed magnets

Within this work, we demonstrate in-situ easy-axis alignment of single-crystal magnetic particles inside a polymer matrix using fused filament fabrication. Two different magnetic materials are investigated: (i) Strontium hexaferrite inside a PA6 matrix, fill grade: 49 vol% and (ii) Samarium iron nitride inside a PA12 matrix, fill grade: 44 vol%. In the presence of the external alignment field, the strontium hexaferrite particles inside the PA6 matrix can be well aligned with a ratio of remanent magnetization to saturation magnetization of 0.7. No significant alignment for samarium iron nitride could be achieved. The results show the feasibility to fabricate magnets with arbitrary and locally defined easy axis using fused filament fabrication since the permanent magnets used for the alignment (or alternatively an electromagnet) can be mounted on a rotatable platform.

preprint2021arXiv

Full analytical solution for the magnetic field of uniformly magnetized cylinder tiles

We present an analytical solution for the magnetic field of a homogeneously magnetized cylinder tile and by extension solutions for full cylinders, rings, cylinder sectors and ring segments. The derivation is done by direct integration in the magnetic surface charge picture. Results are closed-form expressions and elliptic integrals. All special cases are treated individually, which enables the field computation for all possible position arguments. An implementation is provided in Python together with a performance analysis. The implementation is tested against numerical solutions and applied to compute the magnetic field in a discrete Halbach cylinder.

preprint2021arXiv

Magnetostatics and micromagnetics with physics informed neural networks

Partial differential equations and variational problems can be solved with physics informed neural networks (PINNs). The unknown field is approximated with neural networks. Minimizing the residuals of the static Maxwell equation at collocation points or the magnetostatic energy, the weights of the neural network are adjusted so that the neural network solution approximates the magnetic vector potential. This way, the magnetic flux density for a given magnetization distribution can be estimated. With the magnetization as an additional unknown, inverse magnetostatic problems can be solved. Augmenting the magnetostatic energy with additional energy terms, micromagnetic problems can be solved. We demonstrate the use of physics informed neural networks for solving magnetostatic problems, computing the magnetization for inverse problems, and calculating the demagnetization curves for two-dimensional geometries.

preprint2021arXiv

Proposal for a micromagnetic standard problem: domain wall pinning at phase boundaries

We propose a novel micromagnetic standard problem calculating the coercive field for unpinning a domain wall at the interface of a multiphase magnet. This problem is sensitive to discontinuities in material parameters for the exchange interaction, the uniaxial anisotropy, and the spontaneous magnetization. We derive an explicit treatment of jump conditions at material interfaces for the exchange interaction in the finite-difference discretization. The micromagnetic simulation results are compared with analytical solutions and show good agreement. The proposed standard problem is well-suited to test the implementation of both finite-difference and finite-element simulation codes.

preprint2020arXiv

Coexistence of distinct skyrmion phases observed in hybrid ferromagnetic/ferrimagnetic multilayers

Materials hosting magnetic skyrmions at room temperature could enable new computing architectures as well as compact and energetically efficient magnetic storage such as racetrack memories. In a racetrack device, information is coded by the presence/absence of magnetic skyrmions forming a chain that is moved through the device. The skyrmion Hall effect that would eventually lead to an annihilation of the skyrmions at the edges of the racetrack can be suppressed for example by anti-ferromagnetically-coupled skyrmions. However, avoiding modifications of the inter-skyrmion distances in the racetrack remains challenging. As a solution to this issue, a chain of bits could also be encoded by two different solitons such as a skyrmion and a chiral bobber. The major limitation of this approach is that it has solely been realized in B20-type single crystalline material systems that support skyrmions only at low temperatures, thus hindering the efficacy for future applications. Here we demonstrate that a hybrid ferro/ferri/ferromagnetic multilayer system can host two distinct skyrmion phases at room temperature. By matching quantitative magnetic force microscopy data with micromagnetic simulations, we reveal that the two phases represent tubular skyrmions and partial skyrmions (similar to skyrmion bobbers). Furthermore, the tubular skyrmion can be converted into a partial skyrmion. Such multilayer systems may thus serve as a platform for designing skyrmion memory applications using distinct types of skyrmions and potentially for storing information using the vertical dimension in a thin film device.

preprint2020arXiv

Hysteresis-free magnetization reversal of exchange-coupled bilayers with finite magnetic anisotropy

Exchange-coupled structures consisting of ferromagnetic and ferrimagnetic layers become technologically more and more important. We show experimentally the occurrence of completely reversible, hysteresis-free minor loops of [Co(0.2 nm)/Ni(0.4 nm)/Pt(0.6 nm)]$_N$ multilayers exchange-coupled to a 20 nm thick ferrimagnetic Tb$_{28}$Co$_{14}$Fe$_{58}$ layer, acting as hard magnetic pinning layer. Furthermore, we present detailed theoretical investigations by means of micromagnetic simulations and most important a purely analytical derivation for the condition of the occurrence of full reversibility in magnetization reversal. Hysteresis-free loops always occur if a domain wall is formed during the reversal of the ferromagnetic layer and generates an intrinsic hard-axis bias field that overcomes the magnetic anisotropy field of the ferromagnetic layer. The derived condition further reveals that the magnetic anisotropy and the bulk exchange of both layers, as well as the exchange coupling strength and the thickness of the ferromagnetic layer play an important role for its reversibility.

preprint2020arXiv

Learning time-stepping by nonlinear dimensionality reduction to predict magnetization dynamics

We establish a time-stepping learning algorithm and apply it to predict the solution of the partial differential equation of motion in micromagnetism as a dynamical system depending on the external field as parameter. The data-driven approach is based on nonlinear model order reduction by use of kernel methods for unsupervised learning, yielding a predictor for the magnetization dynamics without any need for field evaluations after a data generation and training phase as precomputation. Magnetization states from simulated micromagnetic dynamics associated with different external fields are used as training data to learn a low-dimensional representation in so-called feature space and a map that predicts the time-evolution in reduced space. Remarkably, only two degrees of freedom in feature space were enough to describe the nonlinear dynamics of a thin-film element. The approach has no restrictions on the spatial discretization and might be useful for fast determination of the response to an external field.

preprint2020arXiv

Microscopic origin of magnetization reversal in exchange-coupled ferro-/ferrimagnetic bilayers

In this study, the magnetic reversal process of exchange-coupled bilayer systems, consisting of a ferrimagnetic TbFeCo alloy layer and a ferromagnetic [Co/Ni/Pt]N multilayer, was investigated. In particular, minor loop studies, probing solely the reversal characteristics of the softer ferromagnetic layer, reveal two distinct reversal mechanisms, which depend strongly on the thickness of the ferromagnetic layer. For thick layers, irreversible switching of the macroscopic minor loop is observed. The underlying microscopic origin of this reversal process was studied in detail by high-resolution magnetic force microscopy, showing that the reversal is triggered by in-plane domain walls propagating through the ferromagnetic layer. In contrast, thin ferromagnetic layers show a hysteresis-free reversal, which is nucleation-dominated due to grain-to-grain variations in magnetic anisotropy of the Co/Ni/Pt multilayer and an inhomogeneous exchange coupling with the magnetically hard TbFeCo layer, as confirmed by micromagnetic simulations.

preprint2020arXiv

Spin-canting effects in GMR sensors with wide dynamic field range

Magnetoresistive (xMR) sensors find extensive application in science and industry, replacing Hall sensors in various low field environments. While there have been some efforts in increasing the dynamic field range of xMR sensors, Hall sensors remain to dominate high field applications due to their wide linear range. Using a perpendicular magnetized reference system and an in-plane free layer allows us to overcome this disadvantage of xMR sensors, and, furthermore, investigate spin-canting effects in interlayer exchange coupled perpendicular synthetic antiferromagnets (p-SAF). We created p-SAFs with exchange coupling fields of up to 10 kOe, based on magnetic Co/Pt multilayer systems. The p-SAFs are either designed as "single" p-SAFs, where two Co/Pt multilayers are interlayer exchange coupled via a 4 Å thick Ru spacer, or as "double" p-SAFs, where an additional Co layer is interlayer exchange coupled to the top multilayer. These p-SAFs are used for giant magnetoresistance (GMR) sensors with wide dynamic field range. By using a p-SAF as the reference system and employing an in-plane magnetic layer as the GMR's free layer, the linear range can be effectively increased limited only by the p-SAF's switching fields. Additionally, the magnetic anisotropy of the in-plane free layer is fully controlled, which allows saturation fields by design. Different configurations were investigated, ranging from free layer magnetic saturation at lower to far higher fields than the p-SAF's switching fields. We can show through micromagnetic simulations that certain GMR transfer curves are dominated by spin-canting effects in the interlayer exchange coupled reference system. Finally, our simulation results lay out the correlation of the p-SAF's design parameters and its magnetization reversal behavior.

preprint2019arXiv

3D printing of polymer-bonded anisotropic magnets in an external magnetic field and by a modified production process

The possibility of producing polymer-bonded magnets with the aid of additive processes, such as 3D printing, opens up a multitude of new areas of application. Almost any structures and prototypes can be produced cost-effectively in small quantities. Extending the 3D printing process allows the manufacturing of anisotropic magnetic structures by aligning the magnetic easy axis of ferromagnetic particles inside a paste-like compound material along an external magnetic field. This is achieved by two different approaches: First, the magnetic field for aligning the particles is provided by a permanent magnet. Secondly, the 3D printing process itselfs generates an anisotropic behavior of the structures. An inexpensive and customizable end-user fused filament fabrication 3D printer is used to print the magnetic samples. The magnetical properties of different magnetic anisotropic Sr ferrite and SmFeN materials are investigated and discussed.

preprint2019arXiv

Hybrid FFT algorithm for fast demagnetization field calculations on non-equidistant magnetic layers

In micromagnetic simulations, the demagnetization field is by far the computationally most expensive field component and often a limiting factor in large multilayer systems. We present an exact method to calculate the demagnetization field of magnetic layers with arbitrary thicknesses. In this approach we combine the widely used fast-Fourier-transform based circular convolution method with an explicit convolution using a generalized form of the Newell formulas. We implement the method both for central processors and graphics processors and find that significant speedups for irregular multilayer geometries can be achieved. Using this method we optimize the geometry of a magnetic random-access memory cell by varying a single specific layer thickness and simulate a hysteresis curve to determine the resulting switching field.

preprint2019arXiv

Learning magnetization dynamics

Deep neural networks are used to model the magnetization dynamics in magnetic thin film elements. The magnetic states of a thin film element can be represented in a low dimensional space. With convolutional autoencoders a compression ratio of 1024:1 was achieved. Time integration can be performed in the latent space with a second network which was trained by solutions of the Landau-Lifshitz-Gilbert equation. Thus the magnetic response to an external field can be computed quickly.

preprint2019arXiv

Polymer-bonded anisotropic SrFe$_\text{12}$O$_\text{19}$ filaments for fused filament fabrication

In this publication we describe the extrusion process and the properties of polymer-bonded anisotropic SrFe$_\text{12}$O$_\text{19}$ filaments for fused filament fabrication (FFF). Highly filled polyamide 12 filaments with a filling fraction from 40 vol.% to 55 vol.% are mixed and extruded into filaments with a diameter of 1.75 mm. Such filaments are processable with a conventional FFF 3D printer. No modifications of the 3D printer are necessary. Detailed mechanical and magnetic investigations of printed samples are performed and discussed. In the presence of an external alignment field, the Sr ferrite particles inside the PA12 matrix can be aligned along an external magnetic field. The remanence can be increased by 40 % by printing anisotropic structures. For the 55 vol.% filled filament, a remanence of 212.8 mT and a coercivity of 307.4 mT are measured. The capabilities of printing magnetic anisotropic structures in a complex external field are presented with a Halbach-array arrangement. By the aim of an inverse field model, based on a finite element method, the orientation of the particles and the quality of the print can be estimated by a nondestructive method.

preprint2019arXiv

Systematic parameterization of heat-assisted magnetic recording switching probabilities and the consequences for the resulting SNR

The signal-to-noise ratio (SNR) of a bit series written with heat-assisted magnetic recording (HAMR) on granular media depends on a large number of different parameters. The choice of material properties is essential for the obtained switching probabilities of single grains and therefore for the written bits' quality in terms of SNR. Studies where the effects of different material compositions on transition jitter and the switching probability are evaluated were done, but it is not obvious, how significant those improvements will finally change the received SNR. To investigate that influence, we developed an analytical model of the switching probability phase diagram, which contains independent parameters for, inter alia, transition width, switching probability and curvature. Different values lead to corresponding bit patterns on granular media, where a reader model detects the resulting signal, which is finally converted to a parameter dependent SNR value. For grain diameters between 4 and 8nm, we show an increase of ~10dB for bit lengths between 4 and 12nm, an increase of ~9dB for maximum switching probabilities between 0.64 and 1.00, a decrease of ~5dB for down-track-jitter parameters between 0 and 4nm and an increase of ~1dB for reduced bit curvature. Those results are furthermore compared to the theoretical formulas for the SNR. We obtain a good agreement, even though we show slight deviations.

preprint2019arXiv

The superior role of the Gilbert damping on the signal-to-noise ratio in heat-assisted magnetic recording

In magnetic recording the signal-to-noise ratio (SNR) is a good indicator for the quality of written bits. However, a priori it is not clear which parameters have the strongest influence on the SNR. In this work, we investigate the role of the Gilbert damping on the SNR. Grains consisting of FePt like hard magnetic material with two different grain sizes $d_1=5\,$nm and $d_2=7\,$nm are considered and simulations of heat-assisted magnetic recording (HAMR) are performed with the atomistic simulation program VAMPIRE. The simulations display that the SNR saturates for damping constants larger or equal than 0.1. Additionally, we can show that the Gilbert damping together with the bit length have a major effect on the SNR whereas other write head and material parameters only have a minor relevance on the SNR.

preprint2019arXiv

Write head design for effective curvature reduction in heat-assisted magnetic recording by topology optimization

The reduction of the transition curvature of written bits in heat-assisted magnetic recording (HAMR) is expected to play an important role for the future areal density increase of hard disk drives. Recently a write head design with flipped write and return poles was proposed. In this design a large spatial field gradient of the write head was the key to significantly reduce the transition curvature. In this work we optimized the write pole of a heat-assisted magnetic recording head in order to produce large field gradients as well as large fields in the region of the heat pulse. This is done by topology optimization. The simulations are performed with dolfin-adjoint. For the maximum field gradients of $8.1\,$mT/nm, $8.6\,$mT/nm and $11.8\,$mT/nm, locally resolved footprints of an FePt like hard magnetic recording medium are computed with a coarse-grained Landau-Lifshitz-Bloch (LLB) model and the resulting transition curvature is analysed. Additional simulations with a bilayer structure with $50\%$ hard and $50\%$ soft magnetic material are computed. The results show that for both recording media, the optimized head design does not lead to any significant improvement of the written track. Thus, we analyse the transition curvature for the optimized write heads theoretically with an effective recording time window (ERTW) model. Moreover, we check how higher field gradients influence the curvature reduction. The results show that a simple optimization of the conventional head design design is not sufficient for effective curvature reduction. Instead, new head concepts will be needed to reduce transition curvature.

preprint2016arXiv

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.

preprint2016arXiv

Basic noise mechanisms of heat-assisted-magnetic recording

Heat-assisted magnetic recording (HAMR) is expected to be a key technology to significantly increase the areal storage density of magnetic recording devices. At high temperatures thermally induced noise becomes a major problem, which must be overcome in order to reliably write magnetic bits with narrow transitions. We propose an elementary model based on the effective recording time window (ERTW) to compute the switching probability of bits during HARM of bit-patterned media. With few assumptions this analytical model allows to gain deeper insights into the basic noise mechanisms like AC and DC noise. Finally, we discuss strategies to reduce noise and to increase the areal storage density of both bit-patterned as well as granular media.

preprint2016arXiv

Grain structure dependence of coercivity in thin films

We investigated coercive fields of 200nm x 1200nm x 5nm rectangular nanocrystalline thin films as a function of grain size D using finite elements simulations. To this end, we created granular finite element models with grain sizes ranging from 5nm to 60nm, and performed micromagnetic hysteresis calculations along the y-axis (easy direction) as well as along the x-axis (hard direction). We then used an extended Random Anisotropy model to interpret the results and to illustrate the interplay of random anisotropy and shape induced anisotropy, which is coherent on a much larger scale, in thin films.

preprint2016arXiv

Impact of lattice dynamics on the phase stability of metamagnetic FeRh: Bulk and thin films

We present phonon dispersions, element-resolved vibrational density of states (VDOS) and corresponding thermodynamic properties obtained by a combination of density functional theory (DFT) and nuclear resonant inelastic X-ray scattering (NRIXS) across the metamagnetic transition of B2 FeRh in the bulk material and thin epitaxial films. We see distinct differences in the VDOS of the antiferromagnetic (AF) and ferromagnetic (FM) phase which provide a microscopic proof of strong spin-phonon coupling in FeRh. The FM VDOS exhibits a particular sensitivity to the slight tetragonal distortions present in epitaxial films, which is not encountered in the AF phase. This results in a notable change in lattice entropy, which is important for the comparison between thin film and bulk results. Our calculations confirm the recently reported lattice instability in the AF phase. The imaginary frequencies at the $X$-point depend critically on the Fe magnetic moment and atomic volume. Analyzing these non vibrational modes leads to the discovery of a stable monoclinic ground state structure which is robustly predicted from DFT but not verified in our thin film experiments. Specific heat, entropy and free energy calculated within the quasiharmonic approximation suggest that the new phase is possibly suppressed because of its relatively smaller lattice entropy. In the bulk phase, lattice degrees of freedom contribute with the same sign and in similar magnitude to the isostructural AF-FM phase transition as the electronic and magnetic subsystems and therefore needs to be included in thermodynamic modeling.

preprint2016arXiv

Influence of grain size and exchange interaction on the LLB modeling procedure

Reliably predicting bit-error rates in realistic heat-assisted magnetic recording simulations is a challenging task. Integrating the Landau-Lifshitz-Bloch (LLB) equation can reduce the computational effort to determine the magnetization dynamics in the vicinity of the Curie temperature. If one aims that these dynamics coincide with trajectories calculated from the atomistic Landau-Lifshitz-Gilbert equation, one has to carefully model required temperature dependent material functions such as the zero-field equilibrium magnetization as well as the parallel and normal susceptibilities. We present an extensive study on how these functions depend on grain size and exchange interactions. We show that, if the size or the exchange constant of a reference grain is modified, the material functions can be scaled, according to the changed Curie temperature, yielding negligible errors. This is shown to be valid for volume changes of up to $\pm 40$ % and variations of the exchange constant of up to $\pm10$ %. Besides the temperature dependent material curves, computed switching probabilities also agree well with probabilities separately determined for each system. Our study suggest that there is no need to recalculate the required LLB input functions for each particle. Within the presented limits it is sufficient to scale them to the Curie temperature of the altered system.

preprint2016arXiv

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 is demonstrated.

preprint2015arXiv

A three-dimensional spin-diffusion model for micromagnetics

We implement a finite-element scheme that solves the Landau-Lifshitz-Gilbert equation coupled to a diffusion equation accounting for spin-polarized currents. The latter solves for the spin accumulation not only in magnetic materials but also in nonmagnetic conductors. The presented method incorporates the model by Slonczewski for the description of spin torque in magnetic multilayers as well as the model of Zhang and Li for the description of current driven domain-wall motion. Furthermore it is able to do both resolve the time evolution of the spin accumulation or treat it in an adiabatic fashion by the choice of sufficiently large time steps.

preprint2015arXiv

Areal density optimizations for heat-assisted-magnetic recording of high density bit-patterned media

Heat-assisted-magnetic recording (HAMR) is hoped to be the future recording technique for high density storage devices. Nevertheless, there exist several realizations strategies. With a coarse-grained Landau-Lifshitz-Bloch (LLB) model we investigate in detail benefits and disadvantages of continuous and pulsed laser spot recording of shingled and conventional bit-patterned media. Additionally we compare single phase grains and bits having a bilayer structure with graded Curie temperature, consisting of a hard magnetic layer with high $T_{\mathrm{C}}$ and a soft magnetic one with low $T_{\mathrm{C}}$, respectively. To describe the whole write process as realistic as possible a distribution of the grain sizes and Curie temperatures, a displacement jitter of the head and the bit positions are considered. For all these cases we calculate bit error rates of various grain patterns, temperatures and write head positions to optimize the achievable areal storage density. Within our analysis shingled HAMR with a continuous laser pulse moving over the medium reaches the best results, and thus having the highest potential to become the next generation storage device.

preprint2015arXiv

Heat assisted magnetic recording of bit patterned media beyond 10 Tb/in$^2$

The limits of the areal storage density as can be achieved with heat assisted magnetic recording (HAMR) are still an open issue. We want to address this central question and present the design of a possible bit patterned medium with an areal storage density above 10 Tb/in$^2$. The model uses hard magnetic recording grains with 5 nm diameter and 10 nm height. It assumes a realistic distribution of the Curie temperature of the underlying material as well as a realistic distribution of the grain size and the grain position. In order to compute the areal density we analyze the detailed switching behavior of a recording bit under different external conditions, which allows to compute the bit error rate of a recording process (shingled and conventional) for different grain spacings and write head positions. Hence, we are able to optimize the areal density of the presented medium.

preprint2015arXiv

Macroscopic Simulation of Isotropic Permanent Magnets

Accurate simulations of isotropic permanent magnets require to take the magnetization process into account and consider the anisotropic, nonlinear, and hysteretic material behaviour near the saturation configuration. An efficient method for the solution of the magnetostatic Maxwell equations including the description of isotropic permanent magnets is presented. The algorithm can easily be implemented on top of existing finite element methods and does not require a full characterization of the hysteresis of the magnetic material. Strayfield measurements of an isotropic permanent magnet and simulation results are in good agreement and highlight the importance of a proper description of the isotropic material.

preprint2015arXiv

Superior bit error rate and jitter due to improved switching field distribution in exchange spring magnetic recording

We report two effects that lead to a reduction of the switching field distribution in exchange spring media. The first effect relies on a subtle mechanism of the interplay between exchange coupling between soft and hard layers and anisotropy that allows significant reduction of the switching field distribution in exchange spring media. This effect reduces the switching field distribution by about 30% compared to single-phase media. A second effect is that due to the improved thermal stability of exchange spring media over single-phase media, thermal fluctuation leads to reduced fundamental transition jitter. The influence of this overall improved switching field distribution on the transition jitter in granular recording and the bit error rate in bit-patterned magnetic recording is discussed. The transition jitter in granular recording for a distribution of K1 values of 3% in the hard layer, taking into account thermal fluctuations during recording, is estimated to be a = 0.78 nm, which is similar to the best reported calculated jitter in optimized heat-assisted recording media.

preprint2014arXiv

A Device model for xMR Sensors based on the Stoner-Wohlfarth model

The Stoner-Wohlfarth model provides an efficient analytical model to describe the behavior of magnetic layers within xMR sensors. Combined with a proper description of magneto-resistivity an efficient device model can be derived, which is necessary for an optimal electric circuit design. Parameters of the model are determined by global optimization of an application specific cost function which contains measured resistances for different applied fields. Several application cases are examined and used for validation of the device model. Furthermore the applicability of the SW model is verified by comparison with micromagnetic energy minimization results.

preprint2014arXiv

A full-fledged micromagnetic code in less than 70 lines of NumPy

We present a complete micromagnetic finite-difference code in less than 70 lines of Python. The code makes largely use of the NumPy library and computes the exchange field by finite differences and the demagnetization field with a fast convolution algorithm. Since the magnetization in finite-difference micromagnetics is represented by a multi-dimensional array and the NumPy library features a rich interface for this data structure, the presented code is a good starting point for the development of novel algorithms.

preprint2014arXiv

Calculating thermal stability and attempt frequency of advanced recording structures without free parameters

Ensuring a permanent increase of magnetic storage densities is one of the main challenges in magnetic recording. Conventional approaches based on single phase grains are not suitable to achieve this goal, because their grain volume is limited due to the superparamagnetic limit. Grains with graded anisotropy are the most promising candidates to overcome this limit, providing magnetic memory bits with small volumes, low coercivity and high thermal stability at the same time. Combining micromagnetic simulations with forward flux sampling (FFS), a computational method for rare events that has been recently applied to the magnetic nanostructures, we have determined thermal escape rates and attempt frequencies of a graded media grain and two single phase grains of the same geometry. We find that graded anisotropy can increase the thermal stability of a grain by 12 orders of magnitudes from tens of milliseconds to centuries without changing the coercive field.

preprint2014arXiv

Efficient Energy-minimization in Finite-Difference Micromagnetics: Speeding up Hysteresis Computations

We implement an efficient energy-minimization algorithm for finite-difference micromagnetics that proofs especially useful for the computation of hysteresis loops. Compared to results obtained by time integration of the Landau-Lifshitz-Gilbert equation, a speedup of up to two orders of magnitude is gained. The method is implemented in a finite-difference code running on CPUs as well as GPUs. This setup enables us to compute accurate hysteresis loops of large systems with a reasonable computational effort. As a benchmark we solve the μMag Standard Problem #1 with a high spatial resolution and compare the results to the solution of the Landau-Lifshitz-Gilbert equation in terms of accuracy and computing time.

preprint2014arXiv

Fundamental limits in heat assisted magnetic recording and methods to overcome it with exchange spring structures

The switching probability of magnetic elements for heat assisted recording is investigated. It is found that FePt elements with a diameter of 5 nm and a height of 10nm show, at a field of 0.5 T, thermally written in errors of 12 percent, which is significant too large for bit patterned magnetic recording. Thermally written in errors can be decreased if larger head fields are applied. However, larger fields lead to an increase the fundamental thermal jitter. This leads to a dilemma between thermally written in errors and fundamental thermal jitter. This dilemma can be partly relaxed by increasing the thickness of the FePt film up to 30nm. For realistic head fields, it is found that the fundamental thermal jitter is in the same order of magnitude of the fundamental thermal jitter in conventional recording, which is about 0.5 to 0.8 nm. Composite structures consisting of high Curie top layer and FePt as hard magnetic storage layer can reduce the thermally written in errors to be smaller than 10-4 if the damping constant is increased in the soft layer. Large damping may be realized by doping with rare earth elements. Similar to single FePt grains also in composite structure an increase of switching probability is sacrifices by an increase of thermal jitter. Structures utilizing first order phase transitions breaking the thermal jitter and writeability dilemma are discussed.

preprint2014arXiv

Landau-Lifshitz-Bloch equation for exchange coupled grains

Heat assisted recording is a promising technique to further increase the storage density in hard disks. Multilayer recording grains with graded Curie temperature is discussed to further assist the write process. Describing the correct magnetization dynamics of these grains, from room temperature to far above the Curie point, during a write process is required for the calculation of bit error rates. We present a coarse grained approach based on the Landau-Lifshitz-Bloch (LLB) equation to model exchange coupled grains with low computational effort. The required temperature dependent material properties such as the zero-field equilibrium magnetization as well as the parallel and normal susceptibilities are obtained by atomistic Landau-Lifshitz-Gilbert (LLG) simulations. Each grain is described with one magnetization vector. In order to mimic the atomistic exchange interaction between the grains a special treatment of the exchange field in the coarse grained approach is presented.

preprint2014arXiv

Multiscale modeling in micromagnetics: existence of solutions and numerical integration

Various applications ranging from spintronic devices, giant magnetoresistance sensors, and magnetic storage devices, include magnetic parts on very different length scales. Since the consideration of the Landau-Lifshitz-Gilbert equation (LLG) constrains the maximum element size to the exchange length within the media, it is numerically not attractive to simulate macroscopic parts with this approach. On the other hand, the magnetostatic Maxwell equations do not constrain the element size, but cannot describe the short-range exchange interaction accurately. A combination of both methods allows to describe magnetic domains within the micromagnetic regime by use of LLG and also considers the macroscopic parts by a non-linear material law using the Maxwell equations. In our work, we prove that under certain assumptions on the non-linear material law, this multiscale version of LLG admits weak solutions. Our proof is constructive in the sense that we provide a linear-implicit numerical integrator for the multiscale model such that the numerically computable finite element solutions admit weak $H^1$-convergence (at least for a subsequence) towards a weak solution.

preprint2014arXiv

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 non-magnetic 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.

preprint2013arXiv

Breaking the thermally induced write error in heat assisted recording by using low and high Tc materials

Heat assisted recording is believed as a key technology in order to further increase the areal density of magnetic recording. In the work of Richter et al. [Richter et al. J. Appl. Phys. 111, 033909 (2012)] it is stated that storage densities will be limited to 15 to 20 Tbit/in2 due to thermally induced write errors. In this letter we propose a composite structure consisting of two materials with different Curie temperatures. A hard magnetic layer is on top of a high Tc soft magnetic layer. In this composite material the thermal write error is negligible up to areal densities of 50 to 100 Tbit/in2. It is shown that the effective thermal field gradient, which is reduced in this composite structure, is not relevant for a possible increase of the transition jitter. The transition jitter is dominated by a small distribution of the Curie temperature. The smallest jitter is obtained for the composite structure.

preprint2013arXiv

magnum.fe: A micromagnetic finite-element simulation code based on FEniCS

We have developed a finite-element micromagnetic simulation code based on the FEniCS package called magnum.fe. Here we describe the numerical methods that are applied as well as their implementation with FEniCS. We apply a transformation method for the solution of the demagnetization-field problem. A semi-implicit weak formulation is used for the integration of the Landau-Lifshitz-Gilbert equation. Numerical experiments show the validity of simulation results. magnum.fe is open source and well documented. The broad feature range of the FEniCS package makes magnum.fe a good choice for the implementation of novel micromagnetic finite-element algorithms.

preprint2012arXiv

FFT-based Kronecker product approximation to micromagnetic long-range interactions

We derive a Kronecker product approximation for the micromagnetic long range interactions in a collocation framework by means of separable sinc quadrature. Evaluation of this operator for structured tensors (Canonical format, Tucker format, Tensor Trains) scales below linear in the volume size. Based on efficient usage of FFT for structured tensors, we are able to accelerate computations to quasi linear complexity in the number of collocation points used in one dimension. Quadratic convergence of the underlying collocation scheme as well as exponential convergence in the separation rank of the approximations is proved. Numerical experiments on accuracy and complexity confirm the theoretical results.

preprint2011arXiv

Fundamental switching field distribution of a single domain particle derived from the Néel-Brown model

We present an analytical derivation of the switching field distribution for a single domain particle from the Néel-Brown model in the presence of a linearly swept magnetic field and influenced by thermal fluctuations. We show that the switching field distribution corresponds to a probability density function and can be obtained by solving a master equation for the not-switching probability together with the transition rate for the magnetization according to the Arrhenius-Néel Law. By calculating the first and second moments of the probability density function we succeed in modeling rate-dependent coercivity and the standard deviation of the coercive field. Complementary to the analytical approach, we also present a Monte Carlo simulation for the switching of a macrospin, which allows us to account for the field dependence of the attempt frequency. The results show excellent agreement with results from a Langevin dynamics simulation and therefore point out the importance to include the relevant dependencies in the attempt frequency. However, we conclude that the Néel-Brown model fails to predict switching fields correctly for common field rates and material parameters used in magnetic recording from the loss of normalization of the probability density function. Investigating the transition regime between thermally assisted and dynamic switching will be of future interest regarding the development of new magnetic recording technologies.

preprint2008arXiv

Micromagnetics of single and double point contact spin torque oscillators

In this paper we numerically conduct micromagnetic modelling to optimize computational boundaries of magnetic thin-film elements applicable to single and double point contact spin torque nano-oscillators. Different boundary conditions have been introduced to compensate spin waves reflections at boundaries that are based on extended layers, absorbing boundaries, and focal point methods and are compared with a technique based on scattering theory. A surface roughness boundary model is presented which is modelled according to the Rayleigh criterion to minimize specular reflections at computational boundaries. It is shown that the surface roughness model disperses the reflected spin waves and improves the signal to background noise ratio. The model is tested in comparison to conventional approaches such as extended layer systems, variable damping constant and focal point methods for double point contacts. The surface roughness model gives solutions that are stable in time, in qualitative agreement with experiments and capable to reproduce phenomena such as phase locking in double point contacts.