Source author record

Thomas Schrefl

Thomas Schrefl 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

29works
10topics
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

29 published item(s)

preprint2022arXiv

Description of collective magnetization processes with machine learning models

This work introduces a latent space method to calculate the demagnetization reversal process of multigrain permanent magnets. The algorithm consists of two deep learning models based on neural networks. The embedded Stoner-Wohlfarth method is used as a reduced order model for computing samples to train and test the machine learning approach. The work is a proof of concept of the method since the used microstructures are simple and the only varying parameters are the magnetic anisotropy axes of the hard magnetic grains. A convolutional autoencoder is used for nonlinear dimensionality reduction, in order to reduce the required amount of training samples for predicting the hysteresis loops. We enriched the loss function with physical information about the underlying problem which increases the accuracy of the machine learning approach. A deep learning regressor is operating in the latent space of the autoencoder. The predictor takes a series of previously encoded magnetization states and outputs the next magnetization state along the hysteresis curve. To predict the complete demagnetization loop, we apply a recursive learning scheme.

preprint2022arXiv

Exploring the hysteresis properties of nanocrystalline permanent magnets using deep learning

We demonstrate the use of model order reduction and neural networks for estimating the hysteresis properties of nanocrystalline permanent magnets from microstructure. With a data-driven approach, we learn the demagnetization curve from data-sets created by grain growth and micromagnetic simulations. We show that the granular structure of a magnet can be encoded within a low-dimensional latent space. Latent codes are constructed using a variational autoencoder. The mapping of structure code to hysteresis properties is a multi-target regression problem. We apply deep neural network and use parameter sharing, in order to predict anchor points along the demagnetization curves from the magnet's structure code. The method is applied to study the magnetic properties of nanocrystalline permanent magnets. We show how new grain structures can be generated by interpolation between two points in the latent space and how the magnetic properties of the resulting magnets can be predicted.

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

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

Microstructure Role in Permanent Magnet Eddy Current Losses

The impact of granular microstructure in permanent magnets on eddy current losses are investigated. A numerical homogenization procedure for electrical conductivity is defined. Then, an approximated simple analytical model for the homogenized conductivity able to capture the main features of the geometrical and material dependences is derived. Finally eddy current losses analytical calculations are given, and the two asymptotic expressions for losses in the stationary conduction limit and advanced skin effect limit are derived and discussed.

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.

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

Grain-size dependent demagnetizing factors in permanent magnets

The coercive field of permanent magnets decreases with increasing grain size. The grain size dependence of coercivity is explained by a size dependent demagnetizing factor. In Dy free Nd$_2$Fe$_{14}$B magnets the size dependent demagnetizing factor ranges from 0.2 for a grain size of 55 nm to 1.22 for a grain size of 8300 nm. The comparison of experimental data with micromagnetic simulations suggests that the grain size dependence of the coercive field in hard magnets is due to the non-uniform magnetostatic field in polyhedral grains.

preprint2016arXiv

Influence of defect thickness on the angular dependence of coercivity in rare-earth permanent magnets

The coercive field and angular dependence of the coercive field of single-grain Nd$_{2}$Fe$_{14}$B permanent magnets are computed using finite element micromagnetics. It is shown that the thickness of surface defects plays a critical role in determining the reversal process. For small defect thicknesses reversal is heavily driven by nucleation, whereas with increasing defect thickness domain wall de-pinning becomes more important. This change results in an observable shift between two well-known behavioral models. A similar trend is observed in experimental measurements of bulk samples, where a Nd-Cu infiltration process has been used to enhance coercivity by modifying the grain boundaries. When account is taken of the imperfect grain alignment of real magnets, the single-grain computed results appears to closely match experimental behaviour.

preprint2016arXiv

Micromagnetic simulation of exchange coupled ferri-/ferromagnetic composite in bit patterned media

Ferri-/ferromagnetic exchange coupled composites are promising candidates for bit patterned media because of the ability to control the magnetic properties of the ferrimagnet by its composition. A micromagnetic model for the bilayer system is presented where we also incorporate the microstructural features of both layers. Micromagnetic finite element simulations are performed to investigate the magnetization reversal behaviour of such media. By adding the exchange coupled ferrimagnet to the ferromagnet, the switching field could be reduced by up to $40\,\%$ and also the switching field distribution is narrowed. To reach these significant improvements, an interface exchange coupling strength of $2\,\mathrm{mJ/m^2}$ is required.

preprint2016arXiv

Switching field distribution of exchange coupled ferri-/ferromagnetic composite bit patterned media

We investigate the switching field distribution and the resulting bit error rate of exchange coupled ferri-/ferromagnetic bilayer island arrays by micromagnetic simulations. Using islands with varying microstructure and anisotropic properties, the intrinsic switching field distribution is computed. The dipolar contribution to the switching field distribution is obtained separately by using a model of a triangular patterned island array resembling $1.4\,\mathrm{Tb/in}^2$ bit patterned media. Both contributions are computed for different thickness of the soft exchange coupled ferrimagnet and also for ferromagnetic single phase FePt islands. A bit patterned media with a bilayer structure of FeGd($5\,\mathrm{nm}$)/FePt($5\,\mathrm{nm}$) shows a bit error rate of $10^{-4}$ with a write field of $1.16\,\mathrm{T}$.

preprint2016arXiv

Thermal activation in permanent magnets

The coercive field of permanent magnets decays with temperature. At non-zero temperature the system can overcome a finite energy barrier through thermal fluctuations. Using finite element micromagnetic simulations, we quantify this effect, which reduces coercivity in addition to the decrease of the coercive field associated with the temperature dependence of the anisotropy field, and validate the method through comparison with existing experimental data.

preprint2014arXiv

Micromagnetics of shape anisotropy based permanent magnets

In the search for rare-earth free permanent magnets, various ideas related to shape anisotropy are being pursued. In this work we assess the limits of shape contributions to the reversal stability using micromagnetic simulations. In a first series of tests we altered the aspect ratio of single phase prolate spheroids from 1 to 16. Starting with a sphere of radius $4.3$ times the exchange length $ L_{\mathrm{ex}}$ we kept the total magnetic volume constant as the aspect ratio was modified. For a ferromagnet with zero magnetocrystalline anisotropy the maximum coercive field reached up to $0.5$ times the magnetization $M_{\mathrm{s}}$. Therefore, in materials with moderate uniaxial magnetocrystalline anisotropy, the addition of shape anisotropy could even double the coercive field. Interestingly due to non-uniform magnetization reversal there is no significant increase of the coercive field for an aspect ratio greater than 5. A similar limit of the maximum aspect ratio was observed in cylinders. The coercive field depends on the wire diameter. By decreasing the wire diameter from $8.7 L_{\mathrm{ex}}$ to $2.2 L_{\mathrm{ex}}$ the coercive field increased by 40%. In the cylinders nucleation of a reversed domain starts at the corners at the end. Smoothing the edges can improve the coercive field by about 10%. In further simulations we compacted soft magnetic cylinders into a bulk-like arrangement. Misalignment and magnetostatic interactions cause a spread of $0.1 M_{\mathrm{s}}$ in the switching fields of the rods. Comparing the volume averaged hysteresis loops computed for isolated rods and the hysteresis loop computed for interacting rods, we conclude that magnetostatic interactions reduce the coercive field by up to 20%.

preprint2014arXiv

Non-uniform FFT for the finite element computation of the micromagnetic scalar potential

We present a quasi-linearly scaling, first order polynomial finite element method for the solution of the magnetostatic open boundary problem by splitting the magnetic scalar potential. The potential is determined by solving a Dirichlet problem and evaluation of the single layer potential by a fast approximation technique based on Fourier approximation of the kernel function. The latter approximation leads to a generalization of the well-known convolution theorem used in finite difference methods. We address it by a non-uniform FFT approach. Overall, our method scales O(M + N + N log N) for N nodes and M surface triangles. We confirm our approach by several numerical tests.

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

Guided self-assembly of magnetic beads for biomedical applications

Micromagnetic beads are widely used in biomedical applications for cell separation, drug delivery, and hypothermia cancer treatment. Here we propose to use self-organized magnetic bead structures which accumulate on fixed magnetic seeding points to isolate circulating tumor cells. The analysis of circulating tumor cells is an emerging tool for cancer biology research and clinical cancer management including the detection, diagnosis and monitoring of cancer. Microfluidic chips for isolating circulating tumor cells use either affinity, size or density capturing methods. We combine multiphysics simulation techniques to understand the microscopic behavior of magnetic beads interacting with Nickel accumulation points used in lab-on-chip technologies. Our proposed chip technology offers the possibility to combine affinity and size capturing with special antibody-coated bead arrangements using a magnetic gradient field created by Neodymium Iron Boron permanent magnets. The multiscale simulation environment combines magnetic field computation, fluid dynamics and discrete particle dynamics.

preprint2013arXiv

LaBonte's method revisited: An effective steepest descent method for micromagnetic energy minimization

We present a steepest descent energy minimization scheme for micromagnetics. The method searches on a curve that lies on the sphere which keeps the magnitude of the magnetization vector constant. The step size is selected according to a modified Barzilai-Borwein method. Standard linear tetrahedral finite elements are used for space discretization. For the computation of static hysteresis loops the steepest descent minimizer is faster than a Landau-Lifshitz micromagnetic solver by more than a factor of two. The speed up on a graphic processor is 4.8 as compared to the fastest single-core CPU implementation.

preprint2013arXiv

Simulation of magnetic active polymers for versatile microfluidic devices

We propose to use a compound of magnetic nanoparticles (20-100 nm) embedded in a flexible polymer (Polydimethylsiloxane PDMS) to filter circulating tumor cells (CTCs). The analysis of CTCs is an emerging tool for cancer biology research and clinical cancer management including the detection, diagnosis and monitoring of cancer. The combination of experiments and simulations lead to a versatile microfluidic lab-on-chip device. Simulations are essential to understand the influence of the embedded nanoparticles in the elastic PDMS when applying a magnetic gradient field. It combines finite element calculations of the polymer, magnetic simulations of the embedded nanoparticles and the fluid dynamic calculations of blood plasma and blood cells. With the use of magnetic active polymers a wide range of tunable microfluidic structures can be created. The method can help to increase the yield of needed isolated CTCs.

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.

preprint2012arXiv

Finite Element Modelling of Impedance Spectroscopy Data in Composite Electroceramics

A time domain finite element numerical study of impedance spectroscopy in composite electroceramics is presented. The simulations take into account the complexity of the realistic three dimensional granular structure including grains and grain boundaries. Two types of finite elements are used to represent the bulk (tetrahedrons) and the grain boundary (prisms) phases, which allow using extremely high bulk to grain boundary ratios. Simulations at various frequencies reveal a strong influence of the grain shape on the impedance data of the sample, both at low and at high frequencies.

preprint2012arXiv

Numerical Methods for the Stray-Field Calculation: A Comparison of recently developed Algorithms

Different numerical approaches for the stray-field calculation in the context of micromagnetic simulations are investigated. We compare finite difference based fast Fourier transform methods, tensor grid methods and the finite-element method with shell transformation in terms of computational complexity, storage requirements and accuracy tested on several benchmark problems. These methods can be subdivided into integral methods (fast Fourier transform methods, tensor-grid method) which solve the stray field directly and in differential equation methods (finite-element method), which compute the stray field as the solution of a partial differential equation. It turns out that for cuboid structures the integral methods, which work on cuboid grids (fast Fourier transform methods and tensor grid methods) outperform the finite-element method in terms of the ratio of computational effort to accuracy. Among these three methods the tensor grid method is the fastest. However, the use of the tensor grid method in the context of full micromagnetic codes is not well investigated yet. The finite-element method performs best for computations on curved structures.

preprint2011arXiv

A simple model for calculating magnetic nanowire domain wall fringing fields

We present a new approach to calculating magnetic fringing fields from head-to-head type domain walls in planar magnetic nanowires. In contrast to calculations based on micromagnetically simulated structures the descriptions of the fields are for the most part analytic and thus significantly less time and resource intensive. The models presented begin with an intuitive picture of domain walls, which is built upon in a phenomenological manner. Comparisons with fields calculated using micromagnetic methods show good quantitative agreement.

preprint2011arXiv

A tunable cancer cell filter using magnetic beads: cellular and fluid dynamic simulations

In the field of biomedicine magnetic beads are used for drug delivery and to treat hyperthermia. Here we propose to use self-organized bead structures to isolate circulating tumor cells using lab-on-chip technologies. Typically blood flows past microposts functionalized with antibodies for circulating tumor cells. Creating these microposts with interacting magnetic beads makes it possible to tune the geometry in size, position and shape. We develop a simulation tool that combines micromagnetics, discrete particle dynamics and fluid dynamics, in order to design micropost arrays made of interacting beads. For the simulation of blood flow we use the Lattice-Boltzmann method with immersed elastic blood cell models. Parallelization distributes large fluid and particle dynamic simulations over available resources to reduce overall calculation time.

preprint2011arXiv

Fast stray field computation on tensor grids

A direct integration algorithm is described to compute the magnetostatic field and energy for given magnetization distributions on not necessarily uniform tensor grids. We use an analytically-based tensor approximation approach for function-related tensors, which reduces calculations to multilinear algebra operations. The algorithm scales with N^(4/3) for N computational cells used and with N^(2/3) (sublinear) when magnetization is given in canonical tensor format. In the final section we confirm our theoretical results concerning computing times and accuracy by means of numerical examples.

preprint2011arXiv

Manipulating ultracold atoms with a reconfigurable nanomagnetic system of domain walls

The divide between the realms of atomic-scale quantum particles and lithographically-defined nanostructures is rapidly being bridged. Hybrid quantum systems comprising ultracold gas-phase atoms and substrate-bound devices already offer exciting prospects for quantum sensors, quantum information and quantum control. Ideally, such devices should be scalable, versatile and support quantum interactions with long coherence times. Fulfilling these criteria is extremely challenging as it demands a stable and tractable interface between two disparate regimes. Here we demonstrate an architecture for atomic control based on domain walls (DWs) in planar magnetic nanowires that provides a tunable atomic interaction, manifested experimentally as the reflection of ultracold atoms from a nanowire array. We exploit the magnetic reconfigurability of the nanowires to quickly and remotely tune the interaction with high reliability. This proof-of-principle study shows the practicability of more elaborate atom chips based on magnetic nanowires being used to perform atom optics on the nanometre scale.

preprint2011arXiv

Nanomagnetic engineering of the properties of domain wall atom traps

We have used the results of micromagnetic simulations to investigate the effects of nanowire geometry and domain wall magnetization structure on the characteristic parameters of magnetic atom traps formed by domain walls in planar ferromagnetic nanowires. It is found that when traps are formed in the near-field of a domain wall both nanowire geometry and wall structure have a substantial effect on trap frequency and adiabaticity. We also show that in certain regimes a trap's depth depends only on the amplitude of an externally applied rotating magnetic field, thus allowing it to be tuned independently of the trap's other critical parameters.

preprint2011arXiv

Self-organizing magnetic beads for biomedical applications

In the field of biomedicine magnetic beads are used for drug delivery and to treat hyperthermia. Here we propose to use self-organized bead structures to isolate circulating tumor cells using lab-on-chip technologies. Typically blood flows past microposts functionalized with antibodies for circulating tumor cells. Creating these microposts with interacting magnetic beads makes it possible to tune the geometry in size, position and shape. We developed a simulation tool that combines micromagnetics and discrete particle dynamics, in order to design micropost arrays made of interacting beads. The simulation takes into account the viscous drag of the blood flow, magnetostatic interactions between the magnetic beads and gradient forces from external aligned magnets. We developed a particle-particle particle-mesh method for effective computation of the magnetic force and torque acting on the particles.

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.