Showing 1–25 of 25 results for author: Uzdin, V M

Optimal magnetization switching via spin-orbit torque on the surface of a topological insulator

Authors: Grzegorz J. Kwiatkowski, Ivan P. Miranda, Cecilia M. Holmqvist, Carlo M. Canali, Igor S. Lobanov, Valery M. Uzdin, Andrei Manolescu, Pavel F. Bessarab, Sigurður I. Erlingsson

Abstract: An optimal protocol for the current-induced switching of a perpendicularly magnetized nanoelement placed on the surface of a topological insulator is presented. The time dependence of both in-plane components of the surface current that induces the magnetization reversal via Dirac spin-orbit torque with minimal Joule heating is derived analytically as a function of the required switching time and… ▽ More An optimal protocol for the current-induced switching of a perpendicularly magnetized nanoelement placed on the surface of a topological insulator is presented. The time dependence of both in-plane components of the surface current that induces the magnetization reversal via Dirac spin-orbit torque with minimal Joule heating is derived analytically as a function of the required switching time and material properties. It is demonstrated that a particularly energy-efficient switching is realized for vanishing dampinglike torque. The optimal reversal time providing a tradeoff between the switching speed and energy efficiency is derived. The obtained switching protocol is contrasted with the one realized in heavy-metal systems. Topological insulators provide a highly tunable platform for the realization of energy-efficient magnetization switching. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Heliknoton in a film of cubic chiral magnet

Authors: Vladyslav M. Kuchkin, Nikolai S. Kiselev, Filipp N. Rybakov, Igor S. Lobanov, Stefan Blügel, Valery M. Uzdin

Abstract: Cubic chiral magnets exhibit a remarkable diversity of two-dimensional topological magnetic textures, including skyrmions. However, the experimental confirmation of topological states localized in all three spatial dimensions remains challenging. In this paper, we investigate a three-dimensional topological state called a heliknoton, which is a hopfion embedded into a helix or conic background. We… ▽ More Cubic chiral magnets exhibit a remarkable diversity of two-dimensional topological magnetic textures, including skyrmions. However, the experimental confirmation of topological states localized in all three spatial dimensions remains challenging. In this paper, we investigate a three-dimensional topological state called a heliknoton, which is a hopfion embedded into a helix or conic background. We explore the range of parameters at which the heliknoton can be stabilized under realistic conditions using micromagnetic modeling, harmonic transition state theory, and stochastic spin dynamics simulations. We present theoretical Lorentz TEM images of the heliknoton, which can be used for experimental comparison. Additionally, we discuss the stability of the heliknoton at finite temperatures and the mechanism of its collapse. Our study offers a pathway for future experimental investigations of three-dimensional topological solitons in magnetic crystals. △ Less

Submitted 20 April, 2023; originally announced April 2023.

Comments: 16 pages, 5 figures

Optimal protocol for spin-orbit torque switching of a perpendicular nanomagnet

Authors: Sergei M. Vlasov, Grzegorz J. Kwiatkowski, Igor S. Lobanov, Valery M. Uzdin, Pavel F. Bessarab

Abstract: It is demonstrated by means of the optimal control theory that the energy cost of the spin-orbit torque induced reversal of a nanomagnet with perpendicular anisotropy can be strongly reduced by proper shaping of both in-plane components of the current pulse. The time-dependence of the optimal switching pulse that minimizes the energy cost associated with Joule heating is derived analytically in te… ▽ More It is demonstrated by means of the optimal control theory that the energy cost of the spin-orbit torque induced reversal of a nanomagnet with perpendicular anisotropy can be strongly reduced by proper shaping of both in-plane components of the current pulse. The time-dependence of the optimal switching pulse that minimizes the energy cost associated with Joule heating is derived analytically in terms of the required reversal time and material properties. The optimal reversal time providing a tradeoff between the switching speed and energy efficiency is obtained. A sweet-spot balance between the field-like and damping-like components of the spin-orbit torque is discovered; it permits for a particularly efficient switching by a down-chirped rotating current pulse whose duration does not need to be adjusted precisely. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 6 pages, 3 figures

Lifetime of skyrmions in discrete systems with infinitesimal lattice constant

Authors: M. N. Potkina, I. S. Lobanov, H. Jónsson, V. M. Uzdin

Abstract: Topological protection of chiral magnetic structures is investigated by taking a two-dimensional magnetic skyrmion as an example. The skyrmion lifetime is calculated based on harmonic transition state theory for a discrete lattice model using various values of the ratio of the lattice constant and the skyrmion size. Parameters of the system corresponding to exchange, anisotropy and Dzyaloshinsky-M… ▽ More Topological protection of chiral magnetic structures is investigated by taking a two-dimensional magnetic skyrmion as an example. The skyrmion lifetime is calculated based on harmonic transition state theory for a discrete lattice model using various values of the ratio of the lattice constant and the skyrmion size. Parameters of the system corresponding to exchange, anisotropy and Dzyaloshinsky-Moriya interaction are chosen in such a way as to keep the energy and size of the skyrmion unchanged for small values of the lattice constant, using scaling relations derived from continuous micromagnetic description. The number of magnetic moments included in the calculations reaches more than a million. The results indicate that in the limit of infinitesimal lattice constant, the energy barrier for skyrmion collapse approaches the Belavin-Polyakov lower bound of the energy of a topological soliton in the $σ$-model, the entropy contribution to the pre-exponential factor in the Arrhenius rate expression for collapse approaches a constant and the skyrmion lifetime can, for large enough number of spins, correspond to thermally stable skyrmion at room temperature even without magnetic dipole-dipole interaction. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Hysteresis and Freedericksz thresholds for twisted states in chiral nematic liquid crystals: Minimum-energy path approach

Authors: Semen S. Tenishchev, Ivan M. Tambovtcev, Alexei D. Kiselev, Valery M. Uzdin

Abstract: We study minimum-energy pathways (MEPs) between the branches of metastable helical structures in chiral nematic liquid crystals (CNLCs)subjected to the electric field applied across the cell. By performing stability analysis we have found that, for the branches with non-vanishing half-turn number, the threshold (critical) voltage of the Freedericksz transition is an increasing function of the free… ▽ More We study minimum-energy pathways (MEPs) between the branches of metastable helical structures in chiral nematic liquid crystals (CNLCs)subjected to the electric field applied across the cell. By performing stability analysis we have found that, for the branches with non-vanishing half-turn number, the threshold (critical) voltage of the Freedericksz transition is an increasing function of the free twisting wave number. The curves for the threshold voltage depend on the elastic anisotropy and determine the zero-field critical free twisting number where the director out-of-plane fluctuations destabilize the CNLC helix. For each MEP passing through a first order saddle point we have computed the energy barrier as the energy difference between the saddle-point and the initial structures at different values of the applied field. In our calculations, where the initial approximation for a MEP at the next step was determined by the MEP obtained at the previous step, the electric field dependence of the energy barrier is found to exhibit the hysteresis. This is the hysteresis of electrically driven transition of the saddle-point configuration between the planar and the tilted structures involving out-of-plane director deformations. It turned out that, by contrast to the second-order Freedericksz transition, this transition is first order and we have studied how it depends on the zenithal anchoring energy strength. △ Less

Submitted 4 November, 2020; originally announced November 2020.

Comments: revtex 4-2, 10 pages, 7 figures, a continuation of arXiv:1911.03987 [cond-mat.soft]

The lifetime of micron scale topological chiral magnetic states with atomic resolution

Authors: I. S. Lobanov, V. M. Uzdin

Abstract: A new method for the numerical computation of the lifetimes of magnetic states within harmonic transition state theory (HTST) has been developed. In the simplest case, the system is described by a Heisenberg-like Hamiltonian with short-range interaction. Calculations are performed in Cartesian coordinates. Constraints on the values of magnetic moments are taken into account using Lagrange multipli… ▽ More A new method for the numerical computation of the lifetimes of magnetic states within harmonic transition state theory (HTST) has been developed. In the simplest case, the system is described by a Heisenberg-like Hamiltonian with short-range interaction. Calculations are performed in Cartesian coordinates. Constraints on the values of magnetic moments are taken into account using Lagrange multipliers. The pre-exponential factor in the Arrhenius law in HTST is written in terms of the determinants of the Hessian of energy at the minima and saddle points on the multidimensional energy surface. An algorithm for calculating these determinants without searching for eigenvalues of the Hessian but using recursive relations is proposed. The method allows calculating determinants for systems containing millions of magnetic moments. This makes it possible to calculate the pre-exponential factor and estimate the lifetimes of micron-scale topological structures with atomic resolution, which until now has been impossible using standard approaches. The accuracy of the method is demonstrated by calculating 2D and 3D skyrmionic structures. △ Less

Submitted 3 September, 2021; v1 submitted 15 August, 2020; originally announced August 2020.

Journal ref: Computer Physics Communications, Volume 269, December 2021, page 108136

Stability of in-plane and out-of-plane chiral skyrmions in epitaxial MnSi(111)/Si(111) thin films: surface twists versus easy-plane anisotropy

Authors: Andrey O. Leonov, Ivan M. Tambovtcev, Igor S. Lobanov, Valery M. Uzdin

Abstract: The revisited theoretical phase diagrams for thin films of cubic helimagnets with the easy-plane anisotropy are shown to have different topology as previously reported [Phys. Rev. B 85, 094429 (2012)]. For both in-plane and out-of-plane directions of an applied magnetic field, the phase diagrams exhibit extensive areas of stable skyrmions, which overlap for a wide range of anisotropy parameters. A… ▽ More The revisited theoretical phase diagrams for thin films of cubic helimagnets with the easy-plane anisotropy are shown to have different topology as previously reported [Phys. Rev. B 85, 094429 (2012)]. For both in-plane and out-of-plane directions of an applied magnetic field, the phase diagrams exhibit extensive areas of stable skyrmions, which overlap for a wide range of anisotropy parameters. Although the existence of the out-of-plane skyrmions was contradicted within the previous theoretical models, we prove that additional surface twists lead to their stability, while the moderate easy-plane anisotropy increases the stability range of in-plane skyrmions. Moreover, the interplay between the anisotropy and the surface twists gives rise to a stable spiral state canted with respect to the surfaces. Being absent in bulk helimagnets, this oblique spiral occupies vast areas at the phase diagrams in thin-film nanosystems and serves as a connecting-link between cones and helicoids. Our theory gives clear directions for renewed experimental studies of in-plane and out-of-plane skyrmions in epitaxial MnSi(111)/Si(111) thin films. △ Less

Submitted 15 July, 2020; originally announced July 2020.

Journal ref: Phys. Rev. B 102, 174415 (2020)

Skyrmions in antiferromagnets: thermal stability and the effect of external field and impurities

Authors: Maria N. Potkina, Igor S. Lobanov, Hannes Jónsson, Valery M. Uzdin

Abstract: Calculations of skyrmions in antiferromagnets (AFMs) are presented, and their properties compared with skyrmions in corresponding ferromagnets (FMs). The rates of skyrmion collapse and escape through the boundary of a track, as well as the binding to and collapse at a non-magnetic impurity, are calculated as a function of applied magnetic field. The activation energy for skyrmion annihilation is t… ▽ More Calculations of skyrmions in antiferromagnets (AFMs) are presented, and their properties compared with skyrmions in corresponding ferromagnets (FMs). The rates of skyrmion collapse and escape through the boundary of a track, as well as the binding to and collapse at a non-magnetic impurity, are calculated as a function of applied magnetic field. The activation energy for skyrmion annihilation is the same in AFMs and corresponding FMs in the absence of an applied magnetic field. The pre-exponential factor in the Arrhenius rate law is, however, different because skyrmion dynamics is different in the two systems. An applied magnetic field has opposite effects on skyrmions in the two types of materials. In AFMs the rate of collapse of skyrmions as well as the rate of escape through the edge of a magnetic strip decreases slightly with increasing field, while these rates increase strongly for a skyrmion in the corresponding FMs when the field is directed antiparallel to the magnetization in the center of the skyrmion. A non-magnetic impurity is less likely to trap a skyrmion in AFMs especially in the presence of a magnetic field. This, together with the established fact that a spin polarized current moves skyrmions in AFMs in the direction of the current, while in FMs skyrmions move at an angle to the current, demonstrates that skyrmions in AFMs have several advantageous properties over skyrmions in FMs for memory and spintronic devices. △ Less

Submitted 19 May, 2020; v1 submitted 31 March, 2020; originally announced April 2020.

Efficient Optimization Method for Finding Minimum Energy Paths of Magnetic Transitions

Authors: Aleksei V. Ivanov, Damjan Dagbartsson, Julien Tranchida, Valery M. Uzdin, Hannes Jónsson

Abstract: Efficient algorithms for the calculation of minimum energy paths of magnetic transitions are implemented within the geodesic nudged elastic band (GNEB) approach. While an objective function is not available for GNEB and a traditional line search can, therefore, not be performed, the use of limited memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) and conjugate gradient algorithms in conjunction with… ▽ More Efficient algorithms for the calculation of minimum energy paths of magnetic transitions are implemented within the geodesic nudged elastic band (GNEB) approach. While an objective function is not available for GNEB and a traditional line search can, therefore, not be performed, the use of limited memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) and conjugate gradient algorithms in conjunction with orthogonal spin optimization (OSO) approach is shown to greatly outperform the previously used velocity projection and dissipative Landau-Lifschitz dynamics optimization methods. The implementation makes use of energy weighted springs for the distribution of the discretization points along the path and this is found to improve performance significantly. The various methods are applied to several test problems using a Heisenberg-type Hamiltonian, extended in some cases to include Dzyaloshinskii-Moriya and exchange interactions beyond nearest neighbours. Minimum energy paths are found for magnetization reversals in a nano-island, collapse of skyrmions in two-dimensional layers and annihilation of a chiral bobber near the surface of a three-dimensional magnet. The LBFGS-OSO method is found to outperform the dynamics based approaches by up to a factor of 8 in some cases. △ Less

Submitted 26 November, 2020; v1 submitted 28 January, 2020; originally announced January 2020.

Journal ref: J. Phys.: Condens. Matter 32 345901 (2020)

Multiple minimum energy paths and scenarios of unwinding transitions in chiral nematic liquid crystals

Authors: Semen S. Tenishchev, Alexei D. Kiselev, Aleksei V. Ivanov, Valery M. Uzdin

Abstract: We apply the minimum energy paths (MEPs) approach to study the helix unwinding transition in chiral nematic liquid crystals. A mechanism of the transition is determined by a MEP passing through a first order saddle point on the free energy surface. The energy difference between the saddle point and the initial state gives the energy barrier of the transition. Two starting approximations for the pa… ▽ More We apply the minimum energy paths (MEPs) approach to study the helix unwinding transition in chiral nematic liquid crystals. A mechanism of the transition is determined by a MEP passing through a first order saddle point on the free energy surface. The energy difference between the saddle point and the initial state gives the energy barrier of the transition. Two starting approximations for the paths are used to find the MEPs representing different transition scenarios: (a) the director slippage approximation with in-plane helical structures; and (b) the anchoring breaking approximation that involves the structures with profound out-of-plane director deviations. It is shown that, at sufficiently low voltages, the unwinding transition is solely governed by the director slippage mechanism with the planar saddle point structures. When the applied voltage exceeds its critical value below the threshold of the Freedericksz transition, the additional scenario through the anchoring breaking transitions is found to come into play. For these transitions, the saddle point structure is characterized by out-of-plane deformations localized near the bounding surface. The energy barriers for different paths of transitions are computed as a function of the voltage and the anchoring energy strengths. △ Less

Submitted 10 November, 2019; originally announced November 2019.

Comments: revtex4-1, 13 pages, 12 figures

Journal ref: Phys. Rev. E 100, 062704 (2019)

Real-space observation of skyrmion clusters with mutually orthogonal skyrmion tubes

Authors: Hayley. R. O. Sohn, Sergei M. Vlasov, Valeriy M. Uzdin, Andrey O. Leonov, Ivan I. Smalyukh

Abstract: We report the discovery and direct visualization of skyrmion clusters with mutually-orthogonal orientations of constituent isolated skyrmions in chiral liquid crystals and ferromagnets. We show that the nascent conical state underlies attracting inter-skyrmion potential, whereas an encompassing homogeneous state leads to the repulsive skyrmion-skyrmion interaction. The crossover between different… ▽ More We report the discovery and direct visualization of skyrmion clusters with mutually-orthogonal orientations of constituent isolated skyrmions in chiral liquid crystals and ferromagnets. We show that the nascent conical state underlies attracting inter-skyrmion potential, whereas an encompassing homogeneous state leads to the repulsive skyrmion-skyrmion interaction. The crossover between different regimes of skyrmion interaction could be identified upon changing layer thickness and/or the surface anchoring. We develop a phenomenological theory describing two types of skyrmions and the underlying mechanisms of their interaction. We show that isolated horizontal skyrmions with the same polarity may approach a vertical isolated skyrmion from both sides and thus constitute two energetically-different configurations which are also observed experimentally. In an extreme regime of mutual attraction, the skyrmions wind around each other forming compact superstructures with undulations. We also indicate that our numerical simulations on skyrmion clusters are valid in a parameter range corresponding to the A-phase region of cubic helimagnets. △ Less

Submitted 17 June, 2019; originally announced June 2019.

Journal ref: Phys. Rev. B 100, 104401 (2019)

Stability of Long-lived Antiskyrmions in Mn-Pt-Sn Material

Authors: M. N. Potkina, I. S. Lobanov, O. A. Tretiakov, H. Jónsson, V. M. Uzdin

Abstract: The lifetime of antiskyrmions at room temperature in a Mn-Pt-Sn tetragonal Heusler material has been calculated using an atomic scale representation including nearly a million spins. The evaluation of the pre-exponential factor in the Arrhenius rate expression for this large system is made possible by an implementation of harmonic transition state theory that avoids evaluation of the eigenvalues o… ▽ More The lifetime of antiskyrmions at room temperature in a Mn-Pt-Sn tetragonal Heusler material has been calculated using an atomic scale representation including nearly a million spins. The evaluation of the pre-exponential factor in the Arrhenius rate expression for this large system is made possible by an implementation of harmonic transition state theory that avoids evaluation of the eigenvalues of the Hessian matrix.The parameter values in the extended Heisenberg Hamiltonian, including anisotropic Dzyaloshinskii-Moriya interaction, are chosen to reproduce experimental observations [A. K. Nayak $\it{et\,al.}$, Nature $\textbf{548}$, 561 (2017)], in particular the 150 nm diameter. The calculated results are consistent with the long lifetime observed in the laboratory and this exceptional stability of the antiskyrmions is found to result from large activation energy for collapse due to strong exchange coupling while the pre-exponential factor in the Arrhenius expression for the lifetime is found to have a typical magnitude of 10$^{-12}$ s, despite the large number of spins. The long lifetime is, therefore, found to result from energetic effects rather than entropic effects in this system. △ Less

Submitted 4 October, 2020; v1 submitted 14 June, 2019; originally announced June 2019.

Journal ref: Phys. Rev. B 102, 134430 (2020)

Magnetic Skyrmion Annihilation by Quantum Mechanical Tunneling

Authors: Sergei M. Vlasov, Pavel F. Bessarab, Igor S. Lobanov, Mariia N. Potkina, Valery M. Uzdin, Hannes Jónsson

Abstract: Magnetic skyrmions are nano-scale magnetic states that could be used in various spintronics devices. A central issue is the mechanism and rate of various possible annihilation processes and the lifetime of metastable skyrmions. While most studies have focused on classical over-the-barrier mechanism for annihilation, it is also possible that quantum mechanical tunneling through the energy barrier t… ▽ More Magnetic skyrmions are nano-scale magnetic states that could be used in various spintronics devices. A central issue is the mechanism and rate of various possible annihilation processes and the lifetime of metastable skyrmions. While most studies have focused on classical over-the-barrier mechanism for annihilation, it is also possible that quantum mechanical tunneling through the energy barrier takes place. Calculations of the lifetime of magnetic skyrmions in a two-dimensional lattice are presented and the rate of tunneling compared with the classical annihilation rate. A remarkably strong variation in the crossover temperature and the lifetime of the skyrmion is found as a function of the values of parameters in the extended Heisenberg Hamiltonian, i.e. the out-of-plane anisotropy, Dzyaloshinskii-Moriya interaction (DMI) and applied magnetic field. Materials parameters and conditions are identified where the onset of tunneling could be observed on a laboratory time scale. In particular, it is predicted that skyrmion tunneling could be observed in the PdFe/Ir(111) system when an external magnetic field on the order of 6 T is applied. △ Less

Submitted 28 April, 2020; v1 submitted 1 June, 2019; originally announced June 2019.

Fast and Robust Algorithm for the Energy Minimization of Spin Systems Applied in an Analysis of High Temperature Spin Configurations in Terms of Skyrmion Density

Authors: Aleksei V. Ivanov, Valery M. Uzdin, Hannes Jónsson

Abstract: An algorithm for the minimization of the energy of magnetic systems is presented and applied to the analysis of thermal configurations of a ferromagnet to identify inherent structures, i.e. the nearest local energy minima, as a function of temperature. Over a rather narrow temperature interval, skyrmions appear and reach a high temperature limit for the skyrmion density. In addition, the performan… ▽ More An algorithm for the minimization of the energy of magnetic systems is presented and applied to the analysis of thermal configurations of a ferromagnet to identify inherent structures, i.e. the nearest local energy minima, as a function of temperature. Over a rather narrow temperature interval, skyrmions appear and reach a high temperature limit for the skyrmion density. In addition, the performance of the algorithm is further demonstrated in a self-consistent field calculation of a skyrmion in an itinerant magnet. The algorithm is based on a geometric approach in which the curvature of the spherical domain is taken into account and as a result the length of the magnetic moments is preserved in every iteration. In the limit of infinitesimal rotations, the minimization path coincides with that obtained using damped spin dynamics while the use of limited-memory quasi-newton minimization algorithms, such as the limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) algorithm, significantly accelerates the convergence. △ Less

Submitted 26 November, 2020; v1 submitted 4 April, 2019; originally announced April 2019.

Journal ref: Comput. Phys. Commun. 260, 107749 (2021)

Duplication, collapse and escape of magnetic skyrmions revealed using a systematic saddle point search method

Authors: Gideon P. Müller, Pavel F. Bessarab, Sergei M. Vlasov, Fabian R. Lux, Nikolai S. Kiselev, Stefan Blügel, Valery M. Uzdin, Hannes Jónsson

Abstract: Various transitions that a magnetic skyrmion can undergo are found in calculations using a method for climbing up the energy surface and converging onto first order saddle points. In addition to collapse and escape through a boundary, the method identifies a transition where the skyrmion divides and forms two skyrmions. The activation energy for this duplication process can be similar to that of c… ▽ More Various transitions that a magnetic skyrmion can undergo are found in calculations using a method for climbing up the energy surface and converging onto first order saddle points. In addition to collapse and escape through a boundary, the method identifies a transition where the skyrmion divides and forms two skyrmions. The activation energy for this duplication process can be similar to that of collapse and escape. A tilting of the external magnetic field for a certain time interval is found to induce the duplication process in a dynamical simulation. Such a process could turn out to be an important avenue for the creation of skyrmions in future magnetic devices. △ Less

Submitted 3 August, 2018; v1 submitted 26 July, 2018; originally announced July 2018.

Journal ref: Phys. Rev. Lett. 121, 197202 (2018)

Topological Hall effect for electron scattering on skyrmions in external magnetic field

Authors: K. S. Denisov, I. V. Rozhansky, M. N. Potkina, I. S. Lobanov, E. Lahderanta, V. M. Uzdin

Abstract: We consider topological Hall effect (THE) in thin ferromagnetic films due to electron scattering on magnetic skyrmions in the presence of the relatively strong external magnetic field. We account for the effect of the magnetic field on a skyrmion structure and describe the hallmarks of THE differentiating it from ordinary and anomalous Hall effects. We have found that, although in typical ferromag… ▽ More We consider topological Hall effect (THE) in thin ferromagnetic films due to electron scattering on magnetic skyrmions in the presence of the relatively strong external magnetic field. We account for the effect of the magnetic field on a skyrmion structure and describe the hallmarks of THE differentiating it from ordinary and anomalous Hall effects. We have found that, although in typical ferromagnets the variation of magnetic field in the range 1-5 T substantially affects the skyrmion size, THE changes rather weakly remaining quite robust. Therefore, the magnitude of THE is primarily determined by the skyrmions sheet density $n_{sk}$, being comparable to the magnitude of the ordinary Hall effect (OHE) at $n_{sk}=10^{11}$ cm$^{-2}$. The sign of THE is opposite to that of OHE for skyrmions with positive vorticity, while for antiskyrmions the signs are the same. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. B 98, 214407 (2018)

Energy surface and lifetime of magnetic skyrmions

Authors: Valery M. Uzdin, Maria N. Potkina, Igor S. Lobanov, Pavel F. Bessarab, Hannes Jónsson

Abstract: The stability of skyrmions in various environments is estimated by analyzing the multidimensional surface describing the energy of the system as a function of the directions of the magnetic moments in the system. The energy is given by a Heisenberg-like Hamiltonian that includes Dzyaloshinskii-Moriya interaction, anisotropy and external magnetic field. Local minima on this surface correspond to th… ▽ More The stability of skyrmions in various environments is estimated by analyzing the multidimensional surface describing the energy of the system as a function of the directions of the magnetic moments in the system. The energy is given by a Heisenberg-like Hamiltonian that includes Dzyaloshinskii-Moriya interaction, anisotropy and external magnetic field. Local minima on this surface correspond to the ferromagnetic and skyrmion states. Minimum energy paths (MEP) between the minima are calculated using the geodesic nudged elastic band method. The maximum energy along an MEP corresponds to a first order saddle point on the energy surface and gives an estimate of the activation energy for the magnetic transition, such as creation and annihilation of a skyrmion. The pre-exponential factor in the Arrhenius law for the rate, the so-called attempt frequency, is estimated within harmonic transition state theory where the eigenvalues of the Hessian at the saddle point and the local minima are used to characterize the shape of the energy surface. For some degrees of freedom, so-called 'zero modes', the energy of the system remains invariant. They need to be treated separately and give rise to temperature dependence of the attempt frequency. As an example application of this general theory, the lifetime of a skyrmion in a track of finite width for a PdFe overlayer on a Ir(111) substrate is calculated as a function of track width and external magnetic field. Also, the effect of non-magnetic impurities is studied. Various MEPs for annihilation inside a track, via the boundary of a track and at an impurity are presented. The attempt frequency as well as the activation energy has been calculated for each mechanism to estimate the transition rate as a function of temperature. △ Less

Submitted 1 July, 2017; originally announced July 2017.

Lifetime of racetrack skyrmions

Authors: P. F. Bessarab, G. P. Müller, I. S. Lobanov, F. N. Rybakov, N. S. Kiselev, H. Jónsson, V. M. Uzdin, S. Blügel, L. Bergqvist, A. Delin

Abstract: The skyrmion racetrack is a promising concept for future information technology. There, binary bits are carried by nanoscale spin swirls -- skyrmions -- driven along magnetic strips. Stability of the skyrmions is a critical issue for the realization of this technology. Here we demonstrate that the racetrack skyrmion lifetime can be calculated from first principles as a function of temperature, mag… ▽ More The skyrmion racetrack is a promising concept for future information technology. There, binary bits are carried by nanoscale spin swirls -- skyrmions -- driven along magnetic strips. Stability of the skyrmions is a critical issue for the realization of this technology. Here we demonstrate that the racetrack skyrmion lifetime can be calculated from first principles as a function of temperature, magnetic field and track width. Our method combines harmonic transition state theory extended to include Goldstone modes, with an atomistic spin Hamiltonian parametrized from density functional theory calculations. We demonstrate that two annihilation mechanisms contribute to the skyrmion stability: At low external magnetic field, escape through the track boundary prevails, but a crossover field exists, above which the collapse in the interior becomes dominant. Considering a Pd/Fe bilayer on an Ir(111) substrate as a well-established model system, the calculated lifetime is found to be consistent with reported experimental measurements. Our results open the door for predictive simulations, free from empirical parameters, to aid the design of skyrmion-based information technology. △ Less

Submitted 2 October, 2017; v1 submitted 22 June, 2017; originally announced June 2017.

The effect of confinement and defects on the thermal stability of skyrmions

Authors: Valery M. Uzdin, Maria N. Potkina, Igor S. Lobanov, Pavel F. Bessarab, Hannes Jónsson

Abstract: The stability of magnetic skyrmions against thermal fluctuations and external perturbations is investigated within the framework of harmonic transition state theory for magnetic degrees of freedom. The influence of confined geometry and atomic scale non-magnetic defects on the skyrmion lifetime is estimated. It is shown that a skyrmion on a track has lower activation energy for annihilation and hi… ▽ More The stability of magnetic skyrmions against thermal fluctuations and external perturbations is investigated within the framework of harmonic transition state theory for magnetic degrees of freedom. The influence of confined geometry and atomic scale non-magnetic defects on the skyrmion lifetime is estimated. It is shown that a skyrmion on a track has lower activation energy for annihilation and higher energy for nucleation if the size of the skyrmion is comparable with the width of the track. Two mechanisms of skyrmion annihilation are considered: inside the track and escape through the boundary. For both mechanisms, the dependence of activation energy on the track width is calculated. Non-magnetic defects are found to localize skyrmions in their neighborhood and strongly decrease the activation energy for creation and annihilation. This is in agreement with experimental measurements that have found nucleation of skyrmions in presence of spin-polarized current preferably occurring near structural defects. △ Less

Submitted 8 May, 2017; originally announced May 2017.

Energy surface and minimum energy paths for Fréedericksz transitions in bistable cholesteric liquid crystals

Authors: A. V. Ivanov, P. F. Bessarab, E. V. Aksenova, V. P. Romanov, V. M. Uzdin

Abstract: The multidimensional energy surface of a cholesteric liquid crystal in a planar cell is investigated as a function of spherical coordinates determining the director orientation. Minima on the energy surface correspond to the stable states with particular director distribution. External electric and magnetic fields deform the energy surface and positions of minima. It can lead to the transitions be… ▽ More The multidimensional energy surface of a cholesteric liquid crystal in a planar cell is investigated as a function of spherical coordinates determining the director orientation. Minima on the energy surface correspond to the stable states with particular director distribution. External electric and magnetic fields deform the energy surface and positions of minima. It can lead to the transitions between states, known as the Fréedericksz effect. Transitions can be continuous or discontinuous depending on parameters of the liquid crystal which determine an energy surface. In a case of discontinuous transition when a barrier between stable states is comparable with the thermal energy, the activation transitions may occur and it leads to the modification of characteristics of the Fréedericksz effect with temperature without explicit temperature dependencies of liquid crystal parameters. Minimum energy path between stable states on the energy surface for the Fréedericksz transition is found using geodesic nudged elastic band method. Knowledge of this path, which has maximal statistical weight among all other paths, gives the information about a barrier between stable states and configuration of director orientation during the transition. It also allows one to estimate the stability of states with respect to the thermal fluctuations and their lifetime when the system is close to the Fréedericksz transition. △ Less

Submitted 7 December, 2015; originally announced December 2015.

Journal ref: Phys. Rev. E 93, 042708 (2016)

Qualitative Insight and Quantitative Analysis of the Effect of Temperature on the Coercivity of a Magnetic System

Authors: M. Moskalenko, P. F. Bessarab, V. M. Uzdin, Hannes Jónsson

Abstract: The temperature dependence of the response of a magnetic system to an applied field can be understood qualitatively by considering variations in the energy surface characterizing the system and estimated quantitatively with rate theory. In the system analysed here, Fe/Sm-Co spring magnet, the width of the hysteresis loop is reduced to a half when temperature is raised from 25~K to 300~K. This narr… ▽ More The temperature dependence of the response of a magnetic system to an applied field can be understood qualitatively by considering variations in the energy surface characterizing the system and estimated quantitatively with rate theory. In the system analysed here, Fe/Sm-Co spring magnet, the width of the hysteresis loop is reduced to a half when temperature is raised from 25~K to 300~K. This narrowing can be explained and reproduced quantitatively without invoking temperature dependence of model parameters as has typically been done in previous data analysis. The applied magnetic field lowers the energy barrier for reorientation of the magnetization but thermal activation brings the system over the barrier. A 2-dimensional representation of the energy surface is developed and used to gain insight into the transition mechanism and to demonstrate how the applied field alters the transition path. Our results show the importance of explicitly including the effect of thermal activation when interpreting experiments involving the manipulation of magnetic systems at finite temperature. △ Less

Submitted 7 December, 2015; originally announced December 2015.

Method for finding mechanism and activation energy of magnetic transitions, applied to skyrmion and antivortex annihilation

Authors: Pavel F. Bessarab, Valery M. Uzdin, Hannes Jónsson

Abstract: A method for finding minimum energy paths of transitions in magnetic systems is presented and used to determine the mechanism and estimate the activation energy of skyrmion and antivortex annihiliation in nano-systems. The path is optimized with respect to orientation of the magnetic vectors while their magnitudes are fixed or obtained from separate calculations. The curvature of the configuration… ▽ More A method for finding minimum energy paths of transitions in magnetic systems is presented and used to determine the mechanism and estimate the activation energy of skyrmion and antivortex annihiliation in nano-systems. The path is optimized with respect to orientation of the magnetic vectors while their magnitudes are fixed or obtained from separate calculations. The curvature of the configuration space is taken into account by: (1) using geodesics to evaluate distances and displacements of the system during the optimization, and (2) projecting the path tangent and the magnetic force on the tangent space of the manifold defined by all possible orientations of the magnetic vectors. The method, named geodesic nudged elastic band (GNEB), and its implementation are illustrated with calculations of complex transitions involving annihilation and creation of skyrmion and antivortex states. The lifetime of the latter was determined within harmonic transition state theory using a noncollinear extension of the Alexander-Anderson model. △ Less

Submitted 17 February, 2015; originally announced February 2015.

Calculations of magnetic states and minimum energy paths of transitions using a noncollinear extension of the Alexander-Anderson model and a magnetic force theorem

Authors: Pavel F. Bessarab, Valery M. Uzdin, Hannes Jónsson

Abstract: Calculations of stable and metastable magnetic states as well as minimum energy paths for transitions between states are carried out using a noncollinear extension of the multiple-impurity Alexander-Anderson model and a magnetic force theorem which is derived and used to evaluate the total energy gradient with respect to orientation of magnetic moments -- an important tool for efficient navigation… ▽ More Calculations of stable and metastable magnetic states as well as minimum energy paths for transitions between states are carried out using a noncollinear extension of the multiple-impurity Alexander-Anderson model and a magnetic force theorem which is derived and used to evaluate the total energy gradient with respect to orientation of magnetic moments -- an important tool for efficient navigation on the energy surface. By using this force theorem, the search for stable and metastable magnetic states as well as minimum energy paths revealing the mechanism and activation energy of transitions can be carried out efficiently. For Fe monolayer on W(110) surface, the model gives magnetic moment as well as exchange coupling between nearest and next-nearest neighbors that are in good agreement with previous density functional theory calculations. When applied to nanoscale Fe islands on this surface, the magnetic moment is predicted to be 10\% larger for atoms at the island rim, explaining in part an experimentally observed trend in the energy barrier for magnetization reversal in small islands. Surprisingly, the magnetic moment of the atoms does not change much along the minimum energy path for the transitions, which for islands containing more than 15 atom rows along either $[001]$ or $[1\bar{1}0]$ directions involves the formation of a thin, temporary domain wall. A noncollinear magnetic state is identified in a $7\times 7$ atomic row Fe island where the magnetic moments are arranged in an antivortex configuration with the central ones pointing out of the $(110)$ plane. This illustrates how the model can describe complicated exchange interactions even though it contains only a few parameters. The minimum energy path between this antivortex state and the collinear ground state is also calculated and the thermal stability of the antivortex state estimated. △ Less

Submitted 1 July, 2014; v1 submitted 1 April, 2014; originally announced April 2014.

The Effect of Hydrogen Adsorption on the Magnetic Properties of a Surface Nanocluster of Iron

Authors: Pavel F. Bessarab, Valery M. Uzdin, Hannes Jónsson

Abstract: The effect of hydrogen adsorption on the magnetic properties of an Fe$_3$ cluster immersed in a Cu(111) surface has been calculated using densifty functional theory and the results used to parametrize an Alexander-Anderson model which takes into account the interaction of d-electrons with itinerant electrons. A number of adatom configurations containing one to seven H-atoms were analyzed. The sequ… ▽ More The effect of hydrogen adsorption on the magnetic properties of an Fe$_3$ cluster immersed in a Cu(111) surface has been calculated using densifty functional theory and the results used to parametrize an Alexander-Anderson model which takes into account the interaction of d-electrons with itinerant electrons. A number of adatom configurations containing one to seven H-atoms were analyzed. The sequential addition of hydrogen atoms is found to monotonically reduce the total magnetic moment of the cluster with the effect being strongest when the H-atoms sit at low coordinated sites. Decomposition of the charge density indicates a transfer of 0.4 electrons to each of the H-atoms from both the Fe-atoms and from the copper substrate, irrespective of adsorption site and coverage. The magnetic moment of only the nearest neighbor Fe-atoms is reduced and mainly due to increased population of minority spin d-states. This can be modeled by increased indirect coupling of d-states via the conduction s-band in the Alexander-Anderson model. △ Less

Submitted 25 June, 2013; originally announced June 2013.

Journal ref: Phys. Rev. B 88, 214407 (2013)

A drop of hyperfine field at Sn in Fe/Cr/Sn/Cr multilayers

Authors: A. K. Arzhnikov, L. V. Dobysheva, D. V. Fedorov, V. M. Uzdin

Abstract: The magnetism of the Fe/Cr/Sn/Cr multilayers was studied by first-principles functional theory. The calculations by LAPW and SKKR methods showed that two solutions exist in the $Fe_9/Cr_{14}/Sn/Cr_2$ system. One of them is originated from the antiferromagnetic order in the bulk Cr, and the other is connected with incommensurate spin density wave (ISDW) in Cr. A characteristic feature of this sys… ▽ More The magnetism of the Fe/Cr/Sn/Cr multilayers was studied by first-principles functional theory. The calculations by LAPW and SKKR methods showed that two solutions exist in the $Fe_9/Cr_{14}/Sn/Cr_2$ system. One of them is originated from the antiferromagnetic order in the bulk Cr, and the other is connected with incommensurate spin density wave (ISDW) in Cr. A characteristic feature of this system is realization of ISDW in the Cr film thinner than the half length period of a common ISDW in Cr. In the Cr layers with width lower than a quarter of the wave, the ISDW cannot be the cause of the sharp drop of the HFF at Sn. The calculations of the system $Fe_3/Cr_8/Sn/Cr_8$ with and without allowance for the roughness at the Fe/Cr interface showed that the roughness leads to a significant decrease in the HFF at Sn. △ Less

Submitted 24 November, 2004; originally announced November 2004.

Comments: 5 pages, 4 figures, 2 tables