Periodic Boundary Conditions for Finite-Differentiation-Method Fast-Fourier-Transform Micromagnetics

We describe an accurate periodic boundary condition （PBC） called the symmetric PBC in the calculation of the magnetostatie interaction field in the finite-differentiation-method fast-Fourier-transform （FDM-FFT） micromagneties. The micromagnetic cells in the regular mesh used by the FDM-FFT method are finite-sized elements, but not geometrical points. Therefore, the key PBC operations for FDM-FFT methods are splitting and relocating the micromagnetic cell surfaces to stay symmetrically inside the box of half-total sizes with respect to the origin. The properties of the demagnetizing matrix of the split micromagnetic cells are discussed, and the sum rules of demagnetizing matrix are fulfilled by the symmetric PBC.

Microwave permeability of single cobalt nanotube studied by micromagnetics simulations and generalized Snoek's law

Microwave permeability spectra of single Co nanotube under equilibrium state have been studied by micromagnetics simulation.More than four obvious resonance peaks have been found(11.72,24.20,33.18 and 39.55 GHz).Such large resonance frequency cannot be found in other traditional magnetic materials.The configurations of magnetic moments along the nanotube have been simulated.The results show that the top end of nanotube has a"flow-out"pattern of magnetic moments configuration.The bottom end has a"flow-in"pattern of magnetic moments configuration.The magnetic moments within the main body of nanotube are aligned perfectly along the length of nanotube.The magnitude of natural resonance peak is strongly related to the volume fraction of a zone,which has the same orientation of magnetic moments.Large microwave permeability values have been found for single nanotube.The generalized Snoek’s law has been used to validate the micromagnetics simulations in this paper.

HIGH ACCURACY NUMERICAL METHOD OF THIN-FILM PROBLEMS IN MICROMAGNETICS

In this paper, a new high accuracy numerical method for the thin-film problems of micron and submicron size ferromagnetic elements is proposed. For the computation of stray field, we use the finite element method(FEM) by introducing a semi-discrete artificial boundary condition [1, 2]. In our numerical experiments about the domain patterns and their movement, we can see that the results are accordant to that of experiments and other numerical methods. Our method are very convenient to deal with arbitrary shape of thin films such as a polygon with high accuracy.

Mapping the antiparallel aligned domain rotation by microwave excitation

The evolution process of magnetic domains in response to external fields is crucial for the modern understanding and application of spintronics.In this study,we investigated the domain rotation in stripe domain films of varying thicknesses by examining their response to microwave excitation in four different orientations.The resonance spectra indicate that the rotation field of stripe domain film under an applied magnetic field approaches the field where the resonance mode of sample changes.The saturation field of the stripe domain film corresponds to the field where the resonance mode disappears when measured in the stripe direction parallel to the microwave magnetic field.The results are reproducible and consistent with micromagnetic simulations,providing additional approaches and techniques for comprehending the microscopic mechanisms of magnetic domains and characterizing their rotation.

Shape-influenced non-reciprocal transport of magnetic skyrmions in nanoscale channel

Skyrmions, with their vortex-like structures and inherent topological protection, play a pivotal role in developing innovative low-power memory and logic devices. The efficient generation and control of skyrmions in geometrically confined systems are crucial for the development of skyrmion-based spintronic devices. In this study, we focus on investigating the non-reciprocal transport behavior of skyrmions and their interactions with boundaries of various shapes. The shape of the notch structure in the nanotrack significantly affects the dynamic behavior of magnetic skyrmions. Through micromagnetic simulation, the non-reciprocal transport properties of skyrmions in nanowires with different notch structures are investigated in this work.

Creation and annihilation of artificial magnetic skyrmions with the electric field

Recent theory and experiments show that artificial magnetic skyrmions can be stabilized at room temperature without the need for the external magnetic field,casting strong potentials for the device applications.In this work,we study the electric field manipulation of artificial magnetic skyrmions imprinted by Co disks on CoPt multilayers utilizing the micromagnetic simulations.We find that the reversible annihilation and creation of skyrmions can be realized with the electric field via the strain mediated magnetoelastic coupling.In addition,we also demonstrate controllable manipulation of individual skyrmion,which opens a new platform for constructing magnetic field-free and low-energy dissipation skyrmion based media.

Consistency between domain wall oscillation modes and spin wave modes in nanostrips

Investigations on domain wall(DW) and spin wave(SW) modes in a series of nanostrips with different widths and thicknesses have been carried out using micromagnetic simulation. The simulation results show that the frequencies of SW modes and the corresponding DW modes are consistent with each other if they have the same node number along the width direction. This consistency is more pronounced in wide and thin nanostrips, favoring the DW motion driven by SWs.Further analysis of the moving behavior of a DW driven by SWs is also carried out. The average DW speed can reach a larger value of ~ 140 m/s under two different SW sources. We argue that this study is very meaningful for the potential application of DW motion driven by SWs.

MicroMagnetic.jl:A Julia package for micromagnetic and atomistic simulations with GPU support

MicroMagnetic.jl is an open-source Julia package for micromagnetic and atomistic simulations.Using the features of the Julia programming language,MicroMagnetic.jl supports CPU and various GPU platforms,including NVIDIA,AMD,Intel,and Apple GPUs.Moreover,MicroMagnetic.jl supports Monte Carlo simulations for atomistic models and implements the nudged-elastic-band method for energy barrier computations.With built-in support for double and single precision modes and a design allowing easy extensibility to add new features,MicroMagnetic.jl provides a versatile toolset for researchers in micromagnetics and atomistic simulations.

Investigation on mechanism of magnetization reversal for nanocrystalline Pr-Fe-B permanent magnets by micromagnetic finite element methods

Magnetization configurations were calculated under various magnetic fields for nanocrystalline Pr-Fe-B permanent magnets by micromagnetic finite element method.According to the configurations during demagnetization process, the mechanism of magnetization reversal was analyzed.For the Pr2Fe14B with 10 nm grains or its composite with 10vol.% α-Fe, the coercivity was determined by nucleation of reversed domain that took place at grain boundaries.However, for Pr2Fe14B with 30 nm grains, coercivity was controlled by pinning of the nucle-ated domain.For Pr2Fe14B/α-Fe with 30vol.% α-Fe, the demagnetization behavior was characterized by continuous reversal of α-Fe moment.

Remanence Characteristic of Nanostructure of Hard／Soft Magnetic Multilayered Systems

The relation between microscopic properties (e.g., layer thickness, easy axis orientation) and the macroscopic magnetic properties such as remanent magnetization of the ferromagnetic multilayer system is investigated based on a simple micromagnet approach. We concentrate on a multilayer design with periodic boundary condition, where alternating soft/hard layers build a nanostructured multilayer. For any easy axis direction in the soft and hard layers a simple explicit expression of remanence of the system has been derived analytically. We find that the remanence clearly depends on the thickness of the soft magnetic layer and is nearly independent of the thickness of hard magnetic layer. On the other hand, the remanence increases upon reducing the angle enclosed by the saturation magnetization and the easy axis directions of soft magnetic layer. However, it is unsensitive to the easy axis direction of hard magnetic layer, but there exists a maximum remanence for a certain easy axis direction of hard magnetic layer.

Three-and two-dimensional calculations for the interface anisotropy dependence of magnetic properties of exchange-spring Nd2Fe(14)B/α-Fe multilayers with out-of-plane easy axes

Hysteresis loops,energy products and magnetic moment distributions of perpendicularly oriented Nd2Fe(14)B/α-Fe exchange-spring multilayers are studied systematically based on both three-dimensional(3D)and one-dimensional(1D)micromagnetic methods,focused on the influence of the interface anisotropy.The calculated results are carefully compared with each other.The interface anisotropy effect is very palpable on the nucleation,pinning and coercive fields when the soft layer is very thin.However,as the soft layer thickness increases,the pinning and coercive fields are almost unchanged with the increment of interface anisotropy though the nucleation field still monotonically rises.Negative interface anisotropy decreases the maximum energy products and increases slightly the angles between the magnetization and applied field.The magnetic moment distributions in the thickness direction at various applied fields demonstrate a progress of three-step magnetic reversal,i.e.,nucleation,evolution and irreversible motion of the domain wall.The above results calculated by two models are in good agreement with each other.Moreover,the in-plane magnetic moment orientations based on two models are different.The 3D calculation shows a progress of generation and disappearance of vortex state,however,the magnetization orientations within the film plane calculated by the 1D model are coherent.Simulation results suggest that negative interface anisotropy is necessarily avoided experimentally.

Micromagnetic Simulation of Magnetization Reversal in SmCo_5/Sm_2Co_(17) Magnets

A three-dimensional finite element micromagnetic algorithm was developed to study the magnetization reversal of the SmCo 5/Sm 2Co 17 based magnets. The influences of the microstructure and magnetic parameters on the coercivity were studied based on the model consisting of 64 irregular cells according to the experimental microstructure. Numerical results show that the coercivity increases with increasing the 2∶17-type cell size. Large cell boundary thickness leads to small coercivity. The drop of anisotropy constant of 1∶5 phase leads to the coercivity reducing, while the effect of exchange constant of 1∶5 phase on coercivity is contrary to that of exchange constant. The calculated field dependence of coercivity can be predicted by an inhomogeneous domain-wall pinning model. The microstructure parameter was analyzed by comparing the calculated coercivity.

Effect of exchange coupling on magnetic property in Sm–Co/α-Fe layered system

The hysteresis loops as well as the spin distributions of Sm-Co/a-Fe bilayers have been investigated by both three- dimensional （3D） and one-dimensional （1D） micromagnetic calculations, focusing on the effect of the interface exchange coupling under various soft layer thicknesses ts. The exchange coupling coefficient Alas between the hard and soft ,layers varies from 1.8 x10-6 erg/cm to 0.45 x 10-6 erg/cm, while the soft layer thickness increases from 2 nm to 10 nm. As the exchange coupling decreases, the squareness of the loop gradually deteriorates, both pinning and coercive fields rise up monotonically, and the nucleation field goes down. On the other hand, an increment of the soft layer thickness leads to a significant drop of the nucleation field, the deterioration of the hysteresis loop squareness, and an increase of the remanence. The simulated loops based on the 3D and 1D methods are consistent with each other and in good agreement with the measured loops for Sm-Co/a-Fe multilayers.

Computational Design of Rare-Earth Reduced Permanent Magnets

Multiscale simulation is a key research tool in the quest for new permanent magnets.Starting with first principles methods,a sequence of simulation methods can be applied to calculate the maximum possible coercive field and expected energy density product of a magnet made from a novel magnetic material composition.Iron(Fe)-rich magnetic phases suitable for permanent magnets can be found by means of adaptive genetic algorithms.The intrinsic properties computed by ab initio simulations are used as input for micromagnetic simulations of the hysteresis properties of permanent magnets with a realistic structure.Using machine learning techniques,the magnet’s structure can be optimized so that the upper limits for coercivity and energy density product for a given phase can be estimated.Structure property relations of synthetic permanent magnets were computed for several candidate hard magnetic phases.The following pairs(coercive field(T),energy density product(kJ·m^-3))were obtained for iron-tin-antimony(Fe3Sn0.75Sb0.25):(0.49,290),L10-ordered iron-nickel(L10 FeNi):(1,400),cobalt-iron-tantalum(CoFe6Ta):(0.87,425),and manganese-aluminum(MnAl):(0.53,80).

SIMULATION OF THE MAGNETIZATION REVERSAL PROCESS OF RECTANGLE-SHAPED NiFe FILM ELEMENTS UNDER AN ORTHOGONAL MAGNETIC FIELD

The magnetization reversal process of nano-size rectangle-shaped NiFe film elements with different aspect ratios have been investigated under the orthogonally applied magnetic fields by micromagnetic simulation. Different magnetization reversal modes can appear depending on whether the bias field is applied or not. When there is no bias field, double “C” state is the initial reversal state. However, when there is a bias field, “S” state is the starting mode. The larger the aspect ratio is, the larger the switching field is. But, when the aspect ratio is larger than 3, the increase of the switching field ceases. These results can provide useful information to the application of the patterned NiFe film with rectangular elements.

Magnetization Reversal for Ni Nanowires Studied by Micromagnetic Simulations

The magnetization reversal mechanisms for Ni nanowires with different diameters were investigated by micromagnetic simulations. The results show that the reversal mechanisms are significantly dependeht on the diameter of wire. For very thin wires, the reversal occurs by pseudo-coherent rotation. With increasing diameter, magnetization reversal takes place via different nucleation （the transverse domain wall and the vortex domain wall） and subsequent propagation. The reason of transition from the transverse domain wall to the vortex domain wall is given by analytical studies. With further increase of the diameter, the reversal nuclear domain wall becomes tundishoshaped form. As the diameter increases, the width of wall becomes larger.

Micromagnetic simulations of reversal magnetization in cerium-containing magnets

Single-grain models with different cerium contents or structural parameters have been introduced to investigate the reversal magnetization behaviors in cerium-containing magnets. All the micromagnetic simulations are carried out via the object oriented micromagnetic framework(OOMMF). As for single(Nd,Ce)_2 Fe_(14)B type grain, the coercivity decreases monotonously with the increase of the cerium content. Four types of grain structure have been compared: single(Nd,Ce)_2 Fe_(14)B type, core((Nd,Ce)_2 Fe_(14)B)-shell(Nd_2 Fe_(14)B) type with 2 nm thick shell, core(Ce_2 Fe_(14)B)-shell(Nd_2 Fe_(14)B) type, and core(Nd_2 Fe_(14)B)-shell(Ce_2 Fe_(14)B) type. It is found that core((Nd,Ce)_2 Fe_(14)B)-shell(Nd_2 Fe_(14)B)type grain with 2 nm thick shell always presents the largest coercivity under the same total cerium content. Furthermore,the relationship between the coercivity and the shell thickness t in core((Nd,Ce)_2 Fe_(14)B)-shell(Nd_2 Fe_(14)B) type grain has been studied. When the total cerium content is kept at 20.51 at.%, the analyzed results show that as t varies from 1 nm to 7 nm, the coercivity gradually ascends at the beginning, then quickly descends after reaching the maximum value when t = 5 nm. From the perspective of the positions of nucleation points, the reasons why t affects the coercivity are discussed in detail.

Micromagnetic simulation of Sm–Co/α-Fe/Sm–Co trilayers with various angles between easy axes and the film plane

Hysteresis loops and energy products have been calculated systematically by a three-dimensional （3D） software OOMMF for Sm-Co/α-Fe/Sm-Co trilayers with various thicknesses and β, where β is the angle between the easy axis and the field applied perpendicular to the film plane. It is found that trilayers with a perpendicular anisotropy possess considerably larger coercivities and smaller remanences and energy products compared with those with an in-plane anisotropy. Increase of β leads to a fast decrease of the maximum energy product as well as the drop of both remanence and coercivity. Such a drop is much faster than that in the single-phased hard material, which can explain the significant discrepancy between the experiment and the theoretical energy products. Some modeling techniques have been utilized with spin check procedures performed, which yield results in good agreement with the one-dimensional （1D） analytical and experimental data, justifying our calculations. Further, the calculated nucleation fields according to the 3D calculations are larger than those based on the 1D model, whereas the corresponding coercivity is smaller, leading to more square hysteresis loops and better agreement between experimental data and the theory.

Magnetic vortex gyration mediated by point-contact position

Micromagnetic simulation is employed to study the gyration motion of magnetic vortices in distinct permalloy nanodisks driven by a spin-polarized current. The critical current density for magnetic vortex gyration, eigenfrequency, trajectory, velocity and the time for a magnetic vortex to obtain the steady gyration are analyzed. Simulation results reveal that the magnetic vortices in larger and thinner nanodisks can achieve a lower-frequency gyration at a lower current density in a shorter time. However, the magnetic vortices in thicker nanodisks need a higher current density and longer time to attain steady gyration but with a higher eigenfrequency. We also find that the point-contact position exerts different influences on these parameters in different nanodisks, which contributes to the control of the magnetic vortex gyration. The conclusions of this paper can serve as a theoretical basis for designing nano-oscillators and microwave frequency modulators.

Micromagnetic study of magnetization reversal in inhomogeneous permanent magnets

Macroscopic magnetic properties of magnets strongly depend on the magnetization process and the microstructure of the magnets.Complex materials such as hard-soft exchange-coupled magnets or just real technical materials with impurities and inhomogeneities exhibit complex magnetization behavior.Here we investigate the effects of size,volume fraction,and surroundings of inhomogeneities on the magnetic properties of an inhomogeneous magnetic material via micromagnetic simulations.The underlying magnetization reversal and coercivity mechanisms are revealed.Three different demagnetization characteristics corresponding to the exchange coupling phase,semi-coupled phase,and decoupled phase are found,depending on the size of inhomogeneities.In addition,the increase in the size of inhomogeneities leads to a transition of the coercivity mechanism from nucleation to pinning.This work could be useful for optimizing the magnetic properties of both exchange-coupled nanomagnets and inhomogeneous single-phase magnets.