Previous Article in Journal
Stability of Noncentrosymmetric Square Skyrmion Crystals with Easy-Axis and Easy-Plane Magnetic Anisotropy
Previous Article in Special Issue
Can We Still Find an Ideal Memristor?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Magnetic Hopfions: A Review

by
Konstantin Guslienko
1,2,3
1
Departamento de Polimeros y Materiales Avanzados: Fisica, Quimica y Tecnologia, Universidad del Pais Vasco, UPV/EHU, 20018 San Sebastian, Spain
2
EHU Quantum Center, University of the Basque Country, UPV/EHU, 48940 Leioa, Spain
3
IKERBASQUE, the Basque Foundation for Science, 48009 Bilbao, Spain
Magnetism 2024, 4(4), 383-399; https://doi.org/10.3390/magnetism4040025
Submission received: 25 September 2024 / Revised: 17 November 2024 / Accepted: 18 November 2024 / Published: 20 November 2024
(This article belongs to the Special Issue Mathematical Modelling and Physical Applications of Magnetic Systems)

Abstract

:
Recent advances in the research area of 3D magnetic topological solitons (hopfions) in restricted geometries are reviewed. The description of the magnetic solitons is based on a macroscopic micromagnetic approach and the Landau–Lifshitz equation of the magnetization motion. The concepts of the gauge emergent vector potential and emergent magnetic field are widely used to calculate the 3D topological charge (the Hopf index) of magnetic textures. The relation of the magnetic hopfions with classical field theory is demonstrated, and a special role of the curvilinear toroidal coordinates in the description of the hopfions is underlined. The hopfion stability and dynamics in ferromagnetic films and dots are considered. A critical discussion of calculations of the magnetization emergent magnetic field and the Hopf index of the toroidal magnetic hopfions in restricted geometries is presented.

1. Introduction

Non-linear configurations of the continuous classical field (solitons) play an important role in the different fields of modern physics such as classical field theory, optics, condensed matter physics, etc. The solitons are localized, particle-like objects, and, in many cases, the solitons can be classified by integer numbers related to their topology (topological charges) [1].
Let us assume that a physical system can be described by a classical 3D vector field n r of unit length n 2 = 1 (the order-parameter space is a unit sphere S 2 in 3D space) depending on three spatial coordinates, represented by the vector r . For instance, in the area of ferromagnetism, such a field is the unit magnetization field [2]. The vector field in optics is the Stokes vector field [3,4], whereas in liquid crystals, it is a director field [5,6], etc. There are topologically non-trivial configurations of the vector field n r , which describes an order parameter in the media. Such field configurations can be classified by using a mapping r n r from the coordinate space (r) to the order parameter space n r [7].
The macroscopic equation to describe a magnetization equilibrium configuration and magnetization dynamics in ordered magnetic media (the Landau–Lifshitz equation) is essentially nonlinear and allows for several soliton-type solutions including topological magnetic solitons [2]. Topologically non-trivial magnetization configurations in ferromagnets and ferri- and antiferromagnets, such as domain walls, vortices, skyrmions, and hopfions, etc., are currently a focus of researchers working in the area of solid-state magnetism. However, the existence of the different stable 3D magnetization configurations M(r) and the role of 3D (Hopf index) and 2D (skyrmion number) topological charges in the magnetization dynamics are still subject of intensive research. Nowadays, complicated 3D magnetization configurations in ferromagnetic media can be observed experimentally using electron holography or X-ray magnetic imaging [8,9].
The magnetization field M(r) is the vector order parameter of a ferromagnet. The magnetization configuration M(r) in 3D space is represented by the unit field vector m(r) = M(r)⁄|M(r)| depending on three spatial coordinates r   x 1 , x 2 , x 3 . There are particular cases when the magnetization configurations depend only on one spatial coordinate (domain walls) or two spatial coordinates (magnetic vortices and skyrmions in thin magnetic films or in flat magnetic dots). For stable magnetization configurations m r , it is possible to calculate topological charges, which describe the degrees of mappings (homotopy invariants) of a 1D- ( R 1 ) or 2D coordinate space ( R 2 ) to the unit sphere m2 = 1 in the magnetization space S2(m), i.e., R 1 S2(m), R 2 S2(m). The mapping of the 3D coordinate space R 3 (r) to the magnetization unit sphere m2 = 1 (S2) is more complicated and is related to the Hopf index (3D topological charge). The corresponding magnetization configurations with a non-zero Hopf charge are called magnetic hopfions. The 3D magnetization textures are topologically equivalent if they have the same degree of mapping (the Hopf index). The magnetic energy of the given field configuration E m r can be calculated as a functional of the unit magnetization field m(r) and spatial derivatives of the vector m(r). Such energy includes the exchange, magnetostatic, magnetic anisotropy, etc., energy contributions. The minimization of the energy functional E m r yields some equilibrium (stable or metastable) magnetization field configurations m r . Then, the magnetization dynamics of these equilibrium configurations can be considered on the base of the Landau–Lifshits equation of the magnetization motion, and any parameters, including the topological charges, can be calculated.
In this review, I consider 3D magnetic topological solitons—hopfions and Bloch points. In particular, toroidal hopfions, which were introduced in field theory, are analyzed in detail. Such magnetic solitons reveal intriguing and novel physical properties due to their topologically non-trivial 3D magnetization configurations m r . The article is organized as follows. The basic properties of the hopfions in classical field theory are considered in Section 2. The concept of the emergent electromagnetic field and definition of 3D topological charges of the magnetic textures are presented in Section 3. The static and slow dynamics behaviors of the magnetic hopfions are considered in Section 4. The fast linear and nonlinear hopfion dynamics are considered in Section 5. The review is concluded by a summary in Section 6.

2. Hopfions in the Classical Field Theory

Topological three-dimensional solitons with a non-zero Hopf charge, named hopfions, were first introduced in classical field theory in connection with the non-linear σ -model [7,10,11,12]. The explicit form of the mapping of a 3D space ( R 3 ) to the surface of the unit sphere n 2 = 1 of the classical field n r , R 3 S 2 , was introduced by Hopf in 1931 [13]. Later, it was found that the Hopf charge (a degree of the mapping R 3 S 2 n ) can be expressed as some integral of the function composed by a continuous classical field n r (for instance, the magnetization field m r ) and spatial derivatives of the field n r [14].
Faddeev [15] suggested a new Lagrangian to describe the field n r (sometimes referred to as the Faddeev’s or Faddeev–Skyrme´s Lagrangian), which has stable soliton solutions for the three-component classical field n r in 3D coordinate space. The stable solutions conserve a 3D topological charge (the Hopf charge). The Faddeev–Skyrme´s Lagrangian is a linear combination of two invariants of the rotation group O(3) represented by the spatial derivatives of the vector field n r . Then, de Vega [16] explicitly found the field n r components (calling the field configurations “closed vortices”) for the simplest non-trivial unit Hopf charge in the toroidal coordinates. Such hopfions are now called toroidal hopfions or torus-like vortex rings. Nicole [17] suggested an analytical form of the Hopf mapping of 3D-coordinate space R 3 to the surface of the unit sphere n 2 = 1 , R 3 S 2 n , which allowed for the explicit introduction of analytic equations for the soliton field configuration with the unit Hopf index. Later, motivated by the excellent paper by Faddeev and Niemi [7], a series of papers on the toroidal hopfions in classical field theory were published [10,11,12,18]. Gladikowski et al. [10] introduced the field configurations n r of the toroidal hopfions with an arbitrary Hopf index and presented explicit expressions for the emergent magnetic field and vector potential components in the toroidal coordinates. The toroidal coordinates play a special role in the theory of the hopfions because the analytical equations are essentially simplified in these curvilinear coordinates. The toroidal hopfions attracted considerable attention from researchers because such topological solitons are special stable solutions of the Faddeev–Skyrme´s model [15]. A recent review can be found in Ref. [1].
For calculations of the soliton topological charges, it is important to distinguish localized and non-localized solitons. The field n r approaches a constant vector value n 0 at infinity r for the localized solitons. The field is inhomogeneous at infinity r for non-localized solitons. The condition n r n 0 means that the Hopf index, which characterizes the different homotopy classes for the mapping R 3 S 2 n , is an integer number π 3 S 2 = Z in infinite media [7]. Therefore, the toroidal hopfions, which are described by an integer Hopf index Q H = 0 ,   ± 1 , ± 2   , represent a class of the hopfions in ordered media. The Hopf index of the toroidal hopfions can be represented as a product of two winding numbers [10], the planar winding (the azimuthal vorticity) and the twisting of the field configuration n r along the hopfion tube (the poloidal vorticity), respectively. Due to the property that the vector field is asymptotically trivial, n r n 0 at r , it is possible to compactify the three-dimensional coordinate space to the four-dimensional unit sphere surface, R 3 S 3 . The toroidal hopfions as localized solitons of a classical three-dimensional vector field resemble particle-like objects.
To illustrate the toroidal hopfion field n r , we introduce the toroidal coordinates r η , β , φ . The connection between the cylindrical ρ , φ , z and toroidal η , β , φ coordinates is ρ = a τ 1 s i n h η , z = a τ 1 s i n β , φ = φ , and τ = c o s h η c o s β , where the toroidal parameter η varies from 0 to , the poloidal angle β varies from π to π , the azimuthal angle φ varies from 0 to 2 π , and a is a scale parameter (the hopfion radius in the plane z = 0 ). The toroidal hopfion field components are [19]
n z η = p 1 c o s h 2 m η t a n h 2 n η 1 + c o s h 2 m η t a n h 2 n η , n x r + i n y r = 1 n z 2 η e x p i n φ + m β ,
where p = n z η = 0 is the hopfion polarity, p = ± 1 , and the integer numbers m, n are the hopfion poloidal and azimuthal vorticities, respectively.
The profiles of the out-of-plane field component, n z ρ , of the toroidal hopfions are shown in Figure 1 for the different hopfion vorticities m , n . The profiles can be plotted by substituting the poloidal coordinate η ρ , z = a t a n h 2 a ρ / ρ 2 + z 2 + a 2 to Equation (1). The component n z 0 = n z = 1 and n z a = 1 at the hopfion radius ρ = a in the basal plane z = 0 if the hopfion polarity p = + 1 .
Nowadays, different kinds of hopfions are investigated in the field of condensed matter physics (magnetic media [20], liquid crystals and colloids [5], ferroelectrics [21]) and in electromagnetism and gravitation [22] photonics [3], optics [4], etc. There is a growing interest in 3D inhomogeneous magnetization textures classified by a linking number of the preimages of two different points on an S2(m) unit sphere in the 3D coordinate space ( R 3 ), i.e., by the non-zero Hopf index. For a given magnetization texture m r , it is possible to plot preimages, the 3D curves in real space, r m 1 and r m 2 of two points m 1 and m 2 on the magnetization unit sphere m 2 = 1 and visually check the number of their crossings (a linking number of the preimages). The linking number introduced solely for the toroidal hopfions is an integer number by definition and is equal to the Hopf index of a magnetization texture m r . However, aside from the toroidal hopfions, there are other kinds of hopfions, and the Hopf index is, in general, not integer and cannot be calculated as the linking number of preimages. Although, 3D localized topological magnetic solitons were introduced by Dzyaloshinskii et al. [23] a long time ago, and the reincarnation of interest to such magnetization field 3D textures started after the publication of the papers by Sutcliffe [24,25] relatively recently.
The simplest magnetic toroidal hopfions with the Hopf index Q H = 1   were considered in infinite ferromagnetic films [26,27] and cylindrical dots [25,28,29,30] using the vector field toroidal hopfion ansatz [10,11]. It was shown numerically that the toroidal hopfions with Q H = 1 can be the ground state of circular ferromagnetic nanodots [29] assuming a strong surface uniaxial magnetic anisotropy. The first experimental observation of the 3D magnetization configurations interpreted as magnetic hopfions was carried out in the Ir/Co/Pt multilayer films [31]. It is shown in Ref. [32] that aside from the toroidal magnetic hopfions, there is, in soft magnetic materials, another class of the magnetization textures with a non-zero Hopf index—the Bloch points. The skyrmion-like winding magnetic hopfions are considered in Section 5.
For the simplest magnetic hopfion with the vorticities m = 1 ,   n = 1 , Q H = 1   , the hopfion magnetization field components (1) in the cylindrical coordinates are [11]
m z = 1 8 ρ 2 a 2 a 2 + ρ 2 + z 2 2 ,   m x + i m y = 4 a ρ e x p i φ + φ 0 2 z a + i ρ 2 + z 2 a 2 a 2 + ρ 2 + z 2 2 .
It is evident from Equation (2) that the toroidal hopfion magnetization field m r approaches a constant vector at infinity r , m = 0 , 0 , 1 , as it should be for the localized solitons. The constant angle φ 0 describes a freedom of choosing the directions of the O x and O y coordinate axes in the basal x O y plane. According to Ref. [27], we can distinguish Bloch hopfions ( φ 0 = 0 , π ) and Neel hopfions φ 0 = ± π / 2 analogously to the Bloch and Neel skyrmions [33]. The magnetization rotates in the φ ^ z ^ plane ( m ρ = 0 ) for the Bloch hopfions or in the ρ ^ z ^ plane ( m φ = 0 ) for the Neel hopfions when one goes from the coordinate system origin r = 0 to the hopfion radius ρ = a along the radial direction ρ ^ in the x O y plane.
The in-plane magnetization configurations ( z = 0 ) of the toroidal Bloch and Neel hopfions plotted according to Equation (2) are shown in Figure 2. The spatial distribution of the magnetization components m x , m z of the toroidal Bloch hopfion φ 0 = 0   with the vorticities m , n = 1 , 1 plotted in the x z plane is shown in Figure 3.

3. Emergent Electromagnetic Field and the Hopf Index

An inhomogeneous and moving magnetization texture m r , t results in the appearance of effective magnetic and electric fields, which act on other subsystems. These so-called emergent magnetic and electric fields are related to the spatial and time derivatives of the magnetization field m r , t . Initially, this effect was calculated by Korenman et al. [34] and Volovik [35] as a result of the exchange interaction of the local spins with the spins of itinerant electrons in ferromagnetic metals. Then, it was shown [36] that this is a more general, pure geometrical effect related to the choice of the local moving coordinate frame with the Oz axis parallel to the local instant magnetization m r , t .
The emergent electromagnetic field tensor can be written as
F μ ν m = m · μ m × ν m ,
where μ = / x μ denotes spatial and time derivatives. The indices μ ,   ν = 0 ,   1 ,   2   , 3 , where x 0 = c t and r = x 1 , x 2 , x 3 , correspond to the components of a 3D radius–vector r in an orthogonal coordinate system. The field tensor (represented in the units of Φ 0 / 4 π , where Φ 0 = h / e is the magnetic flux quantum) is related to the emergent field four-component vector potential A = A 0 , A 1 , A 2 , A 3 as F μ ν m = μ A ν m ν A μ m . The emergent magnetic field (sometimes called gyrocoupling density) B = × A , A = A 1 , A 2 , A 3 can be defined as in classical electrodynamics B i m = ε i j k F j k m / 2 [36], where the indices i, j, k mark the spatial coordinates x i . It is important that the emergent field B is a divergence-free vector field, · B = 0 . The flux of the emergent magnetic field B = × A through a closed surface defines a 2D topological charge (skyrmion number).
The emergent magnetic field introduced above differs from the real magnetic field defined in standard electrodynamics and is “fictious” in some sense. Nevertheless, the emergent field results in some experimentally measurable effects. Prominent examples of such effects are the topological Hall effect (influence of the emergent magnetic field on trajectories of the conductivity electrons) and skyrmion Hall effect (appearance of the magnetic soliton velocity component transverse to the driving force direction due to the gyroforce). A theoretical approach to exploiting the topological Hall resistance to electrically detect the magnetic hopfion 3D magnetization textures was recently suggested by Göbel et. al. [30]. The flux of the emergent magnetic field trough the plane z = const in 2D ferromagnets determines the gyrovector, an important parameter to describe the motion of 2D topological solitons, magnetic vortices, and skyrmions. The consequences of the non-zero gyrovector in 2D nanostructures have been investigated many times experimentally. A prominent example is the vortex/skyrmion gyrotropic excitation mode immediately related to the gyrovector (see a recent review [33]).
The dot product A · B of the emergent field vector potential A and the emergent magnetic field B defines a 3D topological charge or the Hopf index [37]
Q H = 1 4 π 2 d 3 r A m · B m
of a 3D magnetization texture m r . There is a gauge freedom in choosing the vector potential A . To have a physical sense, the Hopf index should be gauge invariant. The necessary condition for infinite samples is either the nullification of the emergent magnetic field at the sample borders S, B S = 0 , or the field B being tangential to the sample surface, B · σ = 0 , where σ is the external normal to the surface. The condition B S = 0 is satisfied for the localized solitons (toroidal hopfions, for instance) in infinite samples. The condition B · σ = 0 , in general, is not satisfied at the sample surface. It was shown in Ref. [19] for the toroidal hopfions that only a particular choice of A = 2 γ + A e , A e = 1 c o s Θ Φ , related to the hopfion helicity γ leads to the integer and invariant values of the Hopf indices Q H in infinite samples, allowing for considering them as the degrees of mapping (the linking numbers or integer numbers of crossings of the magnetization configuration m r preimages). Here, Θ r = Θ ρ , z , Φ r = n φ + γ ρ , z , are the spherical angles Θ , Φ of the magnetization, and the hopfion helicity is represented by the poloidal angle β as γ ρ , z = m β ρ , z , β ρ , z = a t a n 2 a z / ρ 2 + z 2 a 2 [19]. The gauge invariance should be considered with respect to this particular choice of the vector potential A = 2 γ + A e . In the case of finite samples and magnetization textures different from the toroidal hopfions, the situation is more complicated. It is difficult to choose a proper form of the emergent field vector potential and prove the gauge invariance of the Hopf index. The important question about an integer Hopf index and its gauge invariance should be carefully investigated for each magnetization texture in a restricted geometry.
Recently, Zheng et al. [38] reported on a direct experimental observation of the magnetic toroidal hopfions forming coupled states with skyrmion tube strings in the submicron FeGe plates. They also provided a theoretical interpretation of the observed hopfions (hopfion rings). Zheng et al. [38] used the hopfion ring magnetization to calculate the Hopf invariant (Hopf index) by applying the concept of the emergent magnetic field. The hopfion topological charge ( Q H ) was calculated for confined samples. However, the large values of Q H = 5 ,   6 ,   10 ,   12 are very strange from the point of view of the theory of 3D magnetization configurations (including the magnetic hopfions [7,19,27]). The calculation method resulted in such large integer values of Q H , and gauge invariance of the calculated values of the Hopf index Q H should be carefully analyzed.
Obviously, the magnetization configurations considered by Zheng et al. are not the so-called toroidal hopfions (torus-like vortex rings) introduced in field theory for infinite media; see Refs. [19,27] and references therein. It was proved that the Hopf index is an integer only for the toroidal hopfions in infinite space (see Section 2). In [38], Zheng et al. considered neither toroidal magnetic hopfions nor infinite samples. Therefore, there is no ground for the speculations about the integer Hopf index and its equivalence to the linking number of their 3D magnetization textures. The Hopf index can be calculated by Equation (4) and is, in general, non-integer for arbitrary 3D magnetization configurations in a finite domain like the ones considered in the paper by Zheng et al. [38].
The strong statement about the integer Hopf index of an arbitrary magnetization configuration in a finite sample was made by Zheng et al. without proof or a proper reference. In modern hopfion theory, such a statement is only valid for the toroidal hopfions in infinite media. The simulations of the equilibrium magnetization configurations m r and the corresponding Hopf indices were conducted by Zheng et al. for finite samples: 0.5 μm diameter circular dot with a thickness of 180 nm; the square plate, 180 × 1000 × 1000 nm or “bulk” sample 700 × 350 × 350 nm. None of the samples used for the simulations in Ref. [38] is an infinite sample.
The Hopf index can be calculated only if an explicit expression for the vector potential A of the emergent magnetic field is known. The authors of Ref. [38] suggested a heuristic equation to calculate the vector potential A. However, they did not discuss the important problem of the vector potential gauge and the gauge invariance of the Hopf index. It was shown in Ref. [27] that to make the Hopf index (4) gauge invariant, the emergent field B(r) components should decay fast enough when approaching zero and increasing the absolute value of the radius vector r up to infinity r . This is possible to implement only in an infinite sample. In other words, the Hopf index is invariant and does not depend on the vector potential gauge in infinite samples. Even if one introduces the artificial boundary condition for the magnetization m in a finite sample at the sample borders ( m = m 0 = c o n s t ) as Zheng et al. [38] suggested, such boundary conditions do not guarantee that the emergent magnetic field B goes to zero at the borders because this field is defined via the magnetization space derivatives, not via the magnetization m itself. There is a specific case of the axially symmetric magnetization configuration m r = m ρ , z in a cylindrical sample [19,27] ( r = ρ , φ , z are the cylindrical coordinates), for which B = 0 at the sample surface if m = m 0 = z ^ ( z ^ is a unit vector in the out-of-plane direction) at the borders. None of the magnetization configurations simulated by Zheng et al. is axially symmetric. Moreover, in many cases, the square plates or “bulk” rectangular samples were used for the simulations. The boundary condition m = z ^ used in Ref. [38] means a strong surface uniaxial magnetic anisotropy. It has no physical sense for such magnetic material as FeGe with a cubical crystal structure [39]. Therefore, we can exclude this particular case of axial symmetry from the consideration. It is important that if B 0 at the sample borders, then the Hopf index is not gauge invariant [27] and cannot be used to characterize a magnetization texture m r . Even if one can suggest an artificial expression for the vector potential like the equation used by Zheng et al. to satisfy the equation B = × A , this does not mean that the introduced vector potential is correct.
An expression, alternative to Equation (4) for the Hopf index, is suggested in Refs. [40,41]. This expression is based on the double volume integral from the emergent magnetic field B(r) components and can be derived as a result of the application of the Helmholtz theorem [42] (the fundamental theorem of vector calculus) to the field B(r). An important consequence of the Helmholtz theorem used in Refs. [40,41] is the integral representation of the vector potential A via the emergent field B, namely
A r = 1 4 π d 3 r R × B r R 3 ,       R = r r .
This expression for A r (5) is different from the expression for the emergent field vector potential given by Zheng et al. [38]. Therefore, the heuristic vector potential introduced by Zheng et al. contradicts the Helmholtz theorem [42] and the integral expressions for the Hopf index used in Refs. [40,41]. The similar wrong heuristic expressions for the vector potential (represented as one-fold integral from the emergent magnetic field B) were used in Refs. [20,43] without any justification.
An expression for the vector potential was used by Liu et al. [28] for calculations of the Hopf index of the circular magnetic dots in the momentum representation. Liu et al. calculated the Hopf index to be 0.96 for their dot parameters. However, Liu et al. [26,28] used the Coulomb gauge for the emergent field vector potential, · A = 0 , which is not compatible with Equation (5) for the potential A r given above. The gauge · A = 0 is unphysical for the particular case of the toroidal hopfions. This can be checked if one uses the toroidal vector potential components calculated explicitly in Ref. [19].
The value of the Hopf index represented by Equation (4) depends essentially on the choice of the vector potential A of the emergent magnetic field. Recently, it was shown [19] that a definite gauge of the vector potential A should be chosen for the toroidal hopfions in infinite samples to guarantee the integer values of the Hopf index Q H . The heuristic vector potential given by Zheng et al. [38] can guarantee neither the integer Hopf indices nor their stability with respect to the vector potential gauge change in the finite samples. It seems that the integer values of Q H simulated in the paper by Zheng et al. [38] are some arbitrary and ungrounded numbers resulting in a confusion of the readers who are not specialists in the theory of 3D magnetic textures (hopfions, in particular). Therefore, the methods and calculations used in Ref. [38] should be reconsidered to bring them in accordance with the hopfion theory.

4. Static Properties of Magnetic Hopfions

Recently, the stability of the toroidal hopfions was studied numerically in chiral ferromagnets FeGe [28,29] and in confined ferroelectric nanoparticles (nanospheres of PbZrTi oxide) [21]. Based on the previous theoretical predictions [28,29], the first experimental observation of magnetic toroidal hopfions was carried out in the Ir/Co/Pt multilayer systems [31] with ultrathin Co layers. A restricted sample cylindrical geometry was used for the chiral ferromagnets FeGe and Ir/Co/Pt with a strong Dzyaloshinskyi-Moriya interaction and out-of-plane uniaxial magnetic anisotropy on the sample faces. The calculated Hopf index Q H for chiral ferromagnets was equal to 1 [29,31] or close to 1 [28]. It was simulated in Ref. [28] that the toroidal hopfions are metastable states of a chiral nanodisk existing at a large-enough disk thickness and radius. Their magnetic energy E m r is higher than the energy of the monopole–antimonopole pair 3D magnetization configuration.
The hopfions were also considered in ferroelectrics. The Hopf index of the polarization field P was defined using a non-standard equation in Ref. [21]. If one assumes that the polarization field P(r) for ferroelectrics is similar to the magnetization field M(r) for ferromagnets, then the definition of the Hopf index should be the emergent magnetic field B(r), not the field M(r) (or P(r)). Therefore, the parameter H defined in Ref. [21] as a volume integral from the dot product P · A is not a Hopf index, although the spatial configuration of the field P(r) is similar to the one for the toroidal hopfion of the magnetization field M(r). It is not a surprise that the Hopf index H is not an integer for the considered spherical samples. The calculated values of the Hopf index are not listed in Ref. [21].
The magnetization configurations of the toroidal magnetic hopfions in infinite media with the integer Hopf index Q H = m n and the arbitrary poloidal and azimuthal hopfion vorticities m, n were explicitly calculated in Ref. [19] using the toroidal and cylindrical coordinates. The calculation method was based on the Hopf mapping’s explicit definition and the concept of the emergent magnetic field defined via spatial derivatives of the magnetization field m r . It was shown that the Hopf index density can be represented as a Jacobian of the transformation from the toroidal to the cylindrical coordinates [19]. The calculated components of the emergent magnetic field and emergent field vector potential can be used, in particular, for explanations of the topological and skyrmion Hall effects of the toroidal magnetic hopfions.
The separate problem is the stability of the different 3D magnetization textures m r . It was proven within field theory [44,45] that for any physical systems with a squared gradient field term in the Lagrangian, there are no stable, stationary, localized solutions in the 3D case for any form of the potential. This statement is known as the Hobart–Derrick theorem. However, stable localized solutions (localized solitons) may exist if any energy contributions linear with respect to spatial derivatives or with higher-order spatial derivatives of the magnetization field are present in the Lagrangian [46,47]. In the theory of magnetism, there are some specific energy terms with the first-order derivatives, so-called Lifshitz invariants, accounting for the Dzyaloshinskii–Moriya interactions (DMI) in ferro-magnetic materials with broken inversion symmetry. Another opportunity to obtain a stable 3D magnetization field configurations is accounting for the higher-order spatial derivatives in the Lagrangian. The term quartic in spatial derivatives was introduced by Skyrme [48] and Faddeev [15] a long time ago within classical field theory. It was shown that the Faddeev–Skyrme Lagrangian
L = 1 2 d 3 r μ n 2 + g F μ ν n 2
has stable 3D localized soliton field solutions in the form of toroidal hopfions [7,10,11,12]. Here, g is the coupling constant.
The toroidal hopfions belong to a class of the localized topological solitons and are characterized by integer values of a 3D topological charge (Hopf index) [13,14]. It was shown recently [20] that the classical Heisenberg model with competing long-range exchange interactions results in quadratic terms in the second spatial derivatives of the magnetization field. Although such Heisenberg model is beyond the standard theory of micromagnetism, it may lead to the stabilization of the toroidal magnetic hopfions. The question is this: is it possible to stabilize the toroidal magnetic hopfions in a ferromagnet with the standard micromagnetic exchange energy (avoiding exotic exchange interactions) due to non-zero DMI terms and/or magnetostatic energy? Such micromagnetic energy contributions are beyond the field theory frameworks and, therefore, the applicability of the Hobart–Derrick theory to the evaluation of the stability of the magnetization field configurations should be reconsidered. The simple scaling analysis accounting for the DMI energy terms was conducted in Ref. [49]. However, this scaling analysis ignored the finite sample sizes and the magnetostatic interaction (which is unavoidable in real ferromagnetic samples). The complicated non-local magnetostatic energy is usually not accounted for in the theory of magnetic skyrmions and hopfions or accounted for in a simplified form. The skyrmions are usually considered in the bulk ferromagnetic crystals without inversion symmetry or in ultrathin magnetic films. In both cases, the magnetostatic energy is reduced to a local form and accounted for as an extra contribution to the magnetic anisotropy energy. Accounting for the magnetostatic energy in relatively thick ferromagnetic dots [50] allows us to stabilize quasi-2D magnetic skyrmions without the presence of any DMI terms if a small out-of-plane magnetic uniaxial anisotropy is included in the energy functional. The magnetostatic energy was not included in the energy functional in Refs. [25,29,30] describing the toroidal magnetic hopfions in thick cylindrical dots. The hopfion’s constant phase angle φ 0 introduced in Section 2 does not contribute to the exchange energy. However, the magnetostatic and DMI energies depend on the particular value of φ 0 and are different for the Bloch and Neel toroidal hopfions. We note that the magnetostatic energy can also lead to the stabilization of other kinds of 3D magnetization textures: the Bloch point hopfions with non-zero Hopf charge or half-hedgehog (3D quasi-Neel skyrmion) magnetization textures even in soft magnetic materials with no DMI [32,51,52].
It is demonstrated in Ref. [53] that the Bloch toroidal hopfion magnetization texture is a metastable state of a thick cylindrical dot or a long cylindrical wire of a finite radius R. The existence of this metastable state is a result of a competition of the exchange and magnetostatic energies. The Dzyaloshinskyi–Moriya exchange energy and magnetic uniaxial anisotropy are of the second importance for the toroidal hopfion stabilization. The important role of the magnetostatic interaction in the toroidal hopfion stabilization in a cylindrical dot was confirmed by simulations [43].

5. Dynamics of Magnetic Hopfions

A numerical study of the spin excitation spectra of the toroidal magnetic hopfions was performed in recent papers [54,55]. Very recently, theoretical papers on the magnetic hopfion dynamics [26,43,56,57] were published.
The spin-polarized electrical current-induced 3D dynamics of the toroidal magnetic hopfions were studied both analytically (the collective coordinates approach) and numerically in Ref. [26] in an infinite frustrated ferromagnet (different signs of the exchange integrals for the nearest neighbor spins). The hopfions exhibit complicated dynamics including a longitudinal motion along the current direction, a transverse motion with respect to the current direction, a rotational motion, and dilation [26]. The skyrmion Hall effect (due to non-zero hopfion gyrovector) is clearly seen in the solutions, although the gyrovector was not explicitly defined. The calculated hopfion velocity is quite small, <10 m/s, for the typical values of the current intensity. The results by Liu et al. [26] are in some contradiction with the hopfion dynamics calculated by Wang et al. in Ref. [27], where the hopfion’s internal structure dynamics were ignored. Wang et al. did not observe any skyrmion Hall effect in the course of the toroidal hopfion motion in magnetic stripes. The hopfions moved under the drive of the spin-transfer torque or spin Hall torque along the nanostripe with velocities up to 20 m/s. The origin of this discrepancy is still unclear. The Hopf index was Q H = 1 in both cases.
The hopfion dynamics depend essentially on the phase angle φ 0 defining the toroidal hopfion static magnetization (2). The spin wave excitation spectra of the Bloch and Neel toroidal hopfions in FeGe nanodisks were investigated numerically in Ref. [43]. It was shown that the Bloch hopfions reveal very dense eigenfrequencies in the low-frequency region (<5 GHz). Only a few eigenfrequencies were detected for the Neel hopfions. The eigenmode spatial distributions were found for some selected spin excitation modes. The modes were mainly localized at the disk edge for the Bloch hopfions or near the dot center for Neel hopfions [43]. The eigenmode classification according to their symmetry or number of the nodes in different spatial directions were not presented. It seems that the Bloch hopfion frequency spectra are numerical artifacts and further careful analysis is necessary. The Hopf index of the Neel hopfions calculated using the concept of the emergent electromagnetic is above 0.9. It is unclear why the Neel hopfions transform from a circularly symmetric shape (azimuthal symmetry) to a square shape, increasing the axial magnetic field [43]. There is no reason to break the azimuthal symmetry of the static magnetization configurations.
The spin excitation modes of the Bloch toroidal hopfions in nanodisks (the radius R is 100 nm and thickness L is 70 nm) were simulated in Ref. [54] in a wide frequency range of 0–20 GHz. The spin eigenmodes were classified according to their symmetry, number of the nodes in the azimuthal and out-of-plane directions, and degree of the mode localization. The dot edge, dot middle, and dot center localized modes were identified. No simple rule establishing a connection of the mode localization with its frequency was found. The edge and middle localized modes had large azimuthal indices. Such complicated inhomogeneous mode patterns mean that the ferromagnetic resonance intensity (average mode volume magnetization) of the modes is small and it will be difficult to observe these modes by the standard ferromagnetic resonance spectroscopy technique. It was found numerically [55] that the toroidal hopfions have distinctly less resonance peaks in comparison with skyrmion tubes. It was also found in Ref. [55] that the hopfion breathing excitation modes (oscillations of the hopfion radius) and rotating spin modes (magnetization rotates in the hopfion basal plane x O y ) hybridize, applying the oscillation magnetic field along the hopfion axis 0z.
As it is well known, the skyrmion Hall effect (existence of a gyroforce, which is perpendicular to the soliton velocity) for magnetic topological solitons is attributed to a non-zero gyrovector. The gyrovector, in the case of the magnetic vortices and skyrmions, determines their low-frequency dynamics in both 3D and 2D cases [58,59]. Therefore, the problem of the calculation of the magnitude and direction of the gyrovector of the magnetic hopfions is of practical importance. The components of the emergent magnetic field in the cylindrical coordinate system, defining the gyrovector G of the toroidal magnetic hopfion, are presented in Ref. [27]. It is proven unambiguously [27] that the gyrovector out-of-plane component G z of an axially symmetric toroidal magnetic hopfion with the Hopf index Q H = 1 is equal to zero. Concerning the zero values of the two in-plane gyrovector components, only a plausible assumption is considered that they have to vanish due to the toroidal hopfion and system symmetry. Nevertheless, this assumption was accepted in a number of subsequent articles, where the equality to zero of all components of the toroidal hopfion gyrovector G was mentioned as a proven fact [56,57,60]. Meanwhile, this issue needs to be clarified, since the possible non-zero values of the gyrovector components affect essentially both the three-dimensional current-induced hopfion translational motion [26] and the spin waves excited over the hopfion background [56] by the external magnetic field. Liu et al. [57] excluded the toroidal hopfion gyrovector from consideration and represented the hopfion dynamics solely in terms of the hopfion emergent magnetic field toroidal T = 1 / 2 d V r × B   and octupole moments. Such an approach resulted in a specific toroidal hopfion dynamic presented in Ref. [57]. The static toroidal hopfion ansatz for the particular case of the Hopf index Q H = 1 in Ref. [26] was initially written incorrectly, and then it was corrected in a paper (Ref. [60]) by the same group.
The calculation approach suggested in Ref. [61] is based on the concept of the emergent magnetic field B r introduced in Section 3 and a calculation of the emergent field components and their volume averages for the toroidal hopfion magnetization texture in the appropriate curvilinear coordinate systems. The simplest non-trivial toroidal hopfion with the Hopf index Q H = 1 in the cylindrical magnetic dot was considered, and the dependencies of the Hopfion gyrovector components on the dot sizes were calculated. It was demonstrated by analytical calculations that the magnetic hopfion gyrovector G is not equal to zero and does not vanish even in the limit of an infinite sample. Namely, two components of the gyrovector in the curvilinear cylindrical coordinates, G z and G φ , are finite. The out of-plane z-component of the hopfion gyrovector ( G z ) goes to zero, increasing the dot radius; however, the in-plane gyrovector φ-component ( G φ ) remains finite [61]. The calculated components of the hopfion emergent magnetic field and gyrovector can be used for calculations of the topological and skyrmion Hall effect of the toroidal magnetic hopfions, respectively. It was recently shown [62] that the toroidal hopfions reveal a Hall motion (skyrmion Hall effect) under the current pulses, while the skyrmionium moves only along the current direction. This is because the gyrovector (2D topological charge) of skyrmionium vanishes, whereas it is not equal to zero for the toroidal magnetic hopfions in accordance with Ref. [61]. The fractional Hopf index Q H = 0.8 was calculated in Ref. [62] due to the restricted geometry of FeGe samples. The current pulses were employed to drive the dynamics of the fractional hopfions. An asymmetric Hall motion of the hopfions with respect to the current direction was detected [62].
Space–time magnetic hopfions are suggested in Ref. [63]. They are treated as 2D magnetic textures (skyrmions) excited by an oscillating magnetic or electric field. The authors considered the coupled dynamics of the skyrmion radius and helicity using the model of the exchange interaction quadratic and quartic on the magnetization spatial derivatives. The authors of Ref. [63] believe that the introduced hopfion topological invariant, the space–time Hopf charge ( Q H = ± 1 ), can be tuned by the applied electric field. The emergent magnetic field B i m = ε i j k m · j m × k m  (i, j, k = x, y, t) and gauge vector potential Awere introduced using the standard equations. However, these equations are valid only if the indices i, j, k used in their definition [63] mark the spatial coordinates, for instance, the Cartesian coordinates (x, y, z). However, one of the indices i, j, k corresponds to time (t) according to Knapman et al. [63]. If the index i = t, then the emergent field component B t has no physical sense because B= B x , B y , B z is a vector in three-dimensional coordinate space (x, y, z). If the indices j = t or k = t, then the vector B is not the emergent magnetic field, and the equation B = × A is not valid. In the case of j = t or k = t, the field B is an emergent electric field E proportional to the time derivative of the unit magnetization vector m. The relation of the emergent electric field defined as E i = F i 0 ( i = x , y , z )   or E = A 0 A / c t and the four-vector potential A = A 0 , A is essentially different from the equation B = × A for the emergent magnetic field. Only the space Hopf index defined by Equation (4) via the dot product A · B has physical sense. The time derivative of the magnetization m and emergent electric field E ~ m / t cannot be used for the calculation of the Hopf index defined by Equation (4). Therefore, the equations defining the space–time Hopf index used by Knapman et al. [63] are wrong and should be corrected.
The hopfions considered by Knapman et al. are a particular case of 3D winding hopfions introduced by Kobayashi et al. [64] in the R 2 × S 1 coordinate space, where the third coordinate S 1 (z) is directed along the sample thickness, L. Knapman et al. [64] substituted the thickness coordinate (z) to time ( t ). There is the standard angular parameterization of the magnetization m r components via the spherical angles Θ r , Φ r . The hopfion magnetization spherical angles should be written in an axially symmetric form, Θ r = Θ ρ , z , Φ r = n φ + γ ρ , z , as it was performed for the azimuthally symmetric toroidal hopfions in Refs. [19,53]. If following Ref. [64] one introduces the hopfion helicity as γ ρ , z = P g z , where the phase angle g z satisfies the boundary condition g L g 0 = 2 π at the sample upper/down faces and P is an integer, then the calculated Hopf index of the winding hopfion is an integer. Kobayashi et al. assumed that the angle Θ ρ , z does not depend on the thickness coordinate z . Therefore, the dependence Θ ρ describes a 2D skyrmion profile in the x O y plane.
Following Kobayashi et al. [64], we introduce the winding hopfion magnetization components as
m z ρ = c o s Θ ρ ,   m x r + i m y r = s i n Θ ρ e i n φ P g z + φ 0 ,
where Θ ρ is the skyrmion magnetization polar angle satisfying the boundary conditions Θ 0 = π , Θ = 0 , n is the azimuthal vorticity, and φ 0 is a constant phase angle.
To plot the winding hopfion magnetization, we take the skyrmion profile in the Belavin–Polyakov form [33] c o s Θ ρ = ρ 2 1 / ρ 2 + 1 ( ρ is taken in the units of the skyrmion radius ρ c , m z ρ c , = 0 ) and the function g z = 2 π z / L in the simplest form satisfying the conditions g L g 0 = 2 π . The magnetization components of the winding hopfions are plotted in Figure 4. We note that the skyrmions at the film surfaces z = 0 ,   L are of the Neel type at φ 0 = 0 (Figure 4A) or Bloch type at φ 0 = π / 2 (Figure 4B).
The calculation by Equation (4) yields Q H = Q P , where Q = n is the skyrmion number [36] (2D topological charge). However, in the finite samples, Q is not an integer and the condition g L g 0 = 2 π , in general, is not satisfied. Therefore, the Hopf index of the winding hopfions can be an arbitrary number. The winding hopfions are very different from the toroidal magnetic hopfions, which have almost integer Hopf indices even in finite samples [28,29].
Very recently, the paper by Saji et al. [65] was published, where the authors rewrote the toroidal hopfion magnetization in the spherical coordinates regardless of whether the hopfion has a cylindrical (azimuthal) symmetry. The validity of the ansatz derived in Ref. [65] for the Hopf index Q H > 1 should be proven. It is also unclear why the emergent vector potential calculated by Saji et al. [65] differs from the vector potential presented in Ref. [59] calculated on the base of field theory [10,11]. In Refs. [26,60,65], the toroidal hopfion magnetization, emergent magnetic field, and vector potential are presented via an arbitrary smooth function f r satisfying the boundary conditions f 0 = 0 , f = π to ensure an integer Hopf index Q H in an infinite sample. Therefore, these important local parameters m r , B r , A r are not completely defined and cannot be used for calculations of any effects related to the toroidal magnetic hopfions. Using the function f r is not necessary because the exact explicit expressions for m r , B r , A r of the toroidal magnetic hopfions in infinite media were found in Ref. [19] for any values of the 3D topological charge Q H .
The authors of Ref. [65] demonstrated that the spin waves excited over the hopfion background experience the emergent electromagnetic field generated by the toroidal hopfion magnetization texture m r . In particular, it was shown that the spin waves propagating along the toroidal hopfion symmetry axis (0z) are deflected by the hopfion magnetic texture. This effect using an analogy with the skyrmion topological Hall effect was named as “magnonic Hall effect”.

6. Summary

I briefly review above the current status of research on 3D magnetic topological solitons (hopfions) in restricted and infinite geometries. The hopfions are initially introduced, similarly to the skyrmions, in classical field theory. Then, these complicated 3D textures are calculated and measured in ferromagnetic media. It is shown that the concepts of the gauge emergent vector potential and emergent magnetic field are absolutely necessary to describe the magnetic hopfions. The toroidal hopfions with the integer Hopf indices found in the field theory occupy a special place among the magnetic hopfions. However, one should be careful using the equations developed for infinite media to calculate the hopfions in finite samples. Except for the toroidal hopfions, other hopfions with a non-zero and non-integer Hopf index are also possible (the Bloch points and winding hopfions, for instance). I note that the Hopf index (3D topological charge) can be used as an important parameter (invariant) to characterize a given 3D magnetization texture. The magnetization textures with the same values of the Hopf index are topologically equivalent. The Hopf index has less physical sense in comparison to a 2D topological charge (skyrmion number) because the latter defines the soliton gyrovector resulting, in particular, in the vortex/skyrmion gyrotropic excitation mode observed experimentally.
The main articles on the toroidal hopfion stability and dynamics in ferromagnetic films and dots are reviewed. Both topics are considered mainly theoretically (numerically) because it is difficult to observe stable hopfions and their dynamics experimentally. No applications of the magnetic hopfions have been realized to the moment. I believe that it can be achieved in the near future.
A critical discussion on the calculations and simulations of the magnetization, emergent magnetic field, and Hopf index of the toroidal magnetic hopfions and other kinds of magnetic hopfions in restricted geometries is presented. The critical point is using proper equations for the hopfion emergent magnetic field and vector potential and proof of gauge invariance of the calculated Hopf index.

Funding

The author acknowledges support by IKERBASQUE (the Basque Foundation for Science). The research was funded in part by the Spanish Ministry of Science, Innovation and Universities grant PID2022-137567NB-C21/AEI/10.13039/501100011033, by the Basque Country government under the scheme “Ayuda a Grupos Consolidados” (Ref. IT1670-22), and by the National Science Centre of Poland, project no UMO-2023/49/B/ST3/02920.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting reported results can be found in the preprints by K.Y. Guslienko, “Emergent magnetic field and vector potential of the toroidal magnetic hopfions”, 6 February 2024, DOI: 10.21203/rs.3.rs-3930623/v1, and “3D Magnetization Textures: Toroidal Magnetic Hopfion Stability in Cylindrical Samples”, 12 December 2023, DOI: 10.20944/preprints202312.0816.v1.

Acknowledgments

The author thanks O. Chubykalo-Fesenko, P. Fischer, A. Hoffmann, G.N. Kakazei, N. Kiselev, R. Knapman, M. Krawczyk, A. Lyberatos, R. Rama-Eiroa, F. Rybakov, I. Smalyukh, E.V. Tartakovskaya, G. Volovik, and J. Zang for stimulating discussions on the magnetic hopfions.

Conflicts of Interest

The author declares no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Shnir, Y.M. Chapter 8: Hopfions. In Topological and Non-Topological Solitons in Scalar Field Theories; Cambridge University Press: Cambridge, UK, 2018; pp. 221–246. [Google Scholar]
  2. Kosevich, A.M.; Ivanov, B.A.; Kovalev, A.S. Magnetic solitons. Phys. Rep. 1990, 194, 117–238. [Google Scholar] [CrossRef]
  3. Shen, Y.; Yu, B.; Wu, H.; Li, C.; Zhu, Z.; Zayats, A.V. Topological transformation and free-space transport of photonic hopfions. Adv. Photonics 2023, 5, 015001. [Google Scholar] [CrossRef]
  4. Sugic, D.; Droop, R.; Otte, E.; Ehrmanntraut, D.; Nori, F.; Ruostekoski, J.; Denz, C.; Dennis, M.R. Particle-like topologies in light. Nat. Commun. 2021, 12, 6785. [Google Scholar] [CrossRef] [PubMed]
  5. Smalyukh, I.I. Review: Knots and other new topological defects in liquid crystals and colloids. Rep. Prog. Phys. 2020, 83, 106601. [Google Scholar] [CrossRef]
  6. Tai, J.S.B.; Ackerman, P.J.; Smalyukh, I.I. Topological transformations of Hopf solitons in chiral ferromagnets and liquid crystals. Proc. Natl. Acad. Sci. USA 2018, 15, 921–926. [Google Scholar] [CrossRef]
  7. Faddeev, L.; Niemi, A. Stable knot-like structures in classical field theory. Nature 1997, 387, 58–61. [Google Scholar] [CrossRef]
  8. Donnelly, C.; Scagnoli, V. Imaging three-dimensional magnetic systems with X-rays. J. Phys. Cond. Matter 2020, 32, 213001. [Google Scholar] [CrossRef]
  9. Donnelly, C.; Metlov, K.L.; Scagnoli, V.; Guizar-Sicairos, M.; Holler, M.; Bingham, N.S.; Raabe, J.; Heyderman, L.J.; Cooper, N.R.; Gliga, S. Experimental observation of vortex rings in a bulk magnet. Nat. Phys. 2021, 17, 316–321. [Google Scholar] [CrossRef]
  10. Gladikowski, J.; Hellmund, M. Static solitons with non-zero Hopf number. Phys. Rev. D 1997, 56, 5194–5199. [Google Scholar] [CrossRef]
  11. Hietarinta, J.; Salo, P. Faddeev-Hopf knots: Dynamics of linked un-knots. Phys. Lett. B 1999, 451, 60–67. [Google Scholar] [CrossRef]
  12. Hietarinta, J.; Salo, P. Ground state in the Faddeev-Skyrme model. Phys. Rev. D 2000, 62, 081701. [Google Scholar] [CrossRef]
  13. Hopf, H. Über die Abbildungen der dreidimensionalen Sphäre auf die Kugelfläche. Math. Analen 1931, 104, 637–665. [Google Scholar] [CrossRef]
  14. Whitehead, J.H.C. An expression of Hopf ’s invariant as an integral. Proc. Natl. Acad. Sci. USA 1947, 33, 117–123. [Google Scholar] [CrossRef] [PubMed]
  15. Faddeev, L.D. Some comments on the many-dimensional solitons. Lett. Math. Phys. 1976, 1, 289–293. [Google Scholar] [CrossRef]
  16. de Vega, H.J. Closed vortices and the Hopf index in classical field theory. Phys. Rev. D 1978, 18, 2945. [Google Scholar] [CrossRef]
  17. Nicole, D.A. Solitons with non-vanishing Hopf index. J. Phys. G Nucl. Phys. 1978, 4, 1363–1369. [Google Scholar] [CrossRef]
  18. Sutcliffe, P. Knots in the Skyrme−Faddeev model. Proc. R. Soc. A 2007, 463, 3001–3020. [Google Scholar] [CrossRef]
  19. Guslienko, K.Y. Emergent magnetic field and vector potential of the toroidal magnetic hopfions. Chaos Solitons Fractals 2023, 174, 113840. [Google Scholar] [CrossRef]
  20. Rybakov, F.N.; Kiselev, N.S.; Borisov, A.B.; Döring, L.; Melcher, C.; Blügel, S. Magnetic hopfions in solids. APL Mater. 2022, 10, 111113. [Google Scholar] [CrossRef]
  21. Luk’yanchuk, I.; Tikhonov, Y.; Razumnaya, A.; Vinokur, V.M. Hopfions emerge in ferroelectrics. Nat. Commun. 2020, 11, 2433. [Google Scholar] [CrossRef]
  22. Thompson, A.; Wickes, A.; Swearngin, J.; Bouwmeester, D. Classification of electromagnetic and gravitational hopfions by algebraic type. J. Phys. A Math. Theor. 2015, 48, 205202. [Google Scholar] [CrossRef]
  23. Dzyaloshinskii, I.E.; Ivanov, B.A. Localized topological solitons in a ferromagnet. JETP Lett. 1979, 29, 540–542. [Google Scholar]
  24. Sutcliffe, P. Skyrmion knots in frustrated magnets. Phys. Rev. Lett. 2017, 118, 247203. [Google Scholar] [CrossRef] [PubMed]
  25. Sutcliffe, P. Hopfions in chiral magnets. J. Phys. A Math. Theor. 2018, 51, 375401. [Google Scholar] [CrossRef]
  26. Liu, Y.; Hou, W.; Han, X.; Zang, J. Three-dimensional dynamics of a magnetic hopfion driven by spin transfer torque. Phys. Rev. Lett. 2020, 124, 127204. [Google Scholar] [CrossRef]
  27. Wang, X.S.; Qaiumzadeh, A.; Brataas, A. Current-driven dynamics of magnetic hopfions. Phys. Rev. Lett. 2019, 123, 147203. [Google Scholar] [CrossRef]
  28. Liu, Y.; Lake, R.K.; Zang, J. Binding a hopfion in a chiral magnet nanodisk. Phys. Rev. B 2018, 98, 174437. [Google Scholar] [CrossRef]
  29. Tai, J.-S.B.; Smalyukh, I.I. Static Hopf solitons and knotted emergent fields in solid-state non-centrosymmetric magnetic nanostructures. Phys. Rev. Lett. 2018, 121, 187201. [Google Scholar] [CrossRef]
  30. Göbel, B.; Akosa, C.A.; Tatara, G.; Mertig, I. Topological Hall signatures of magnetic hopfions. Phys. Rev. Res. 2020, 2, 013315. [Google Scholar] [CrossRef]
  31. Kent, N.; Reynolds, N.; Raftrey, D.; Campbell, I.T.G.; Virasawmy, S.; Dhuey, S.; Chopdekar, R.V.; Hierro-Rodriguez, A.; Sorrentino, A.; Pereiro, E.; et al. Creation and observation of Hopfions in magnetic multilayer systems. Nat. Commun. 2021, 12, 1562. [Google Scholar] [CrossRef]
  32. Tejo, F.; Hernández Heredero, R.; Chubykalo-Fesenko, O.; Guslienko, K.Y. The Bloch point 3D topological charge induced by the magnetostatic interaction. Sci. Rep. 2021, 11, 21714. [Google Scholar] [CrossRef] [PubMed]
  33. Guslienko, K.Y. Magnetic vortices and skyrmions. J. Magn. 2019, 24, 549–567. [Google Scholar] [CrossRef]
  34. Korenman, V.; Murray, J.L.; Prange, R.E. Local-band theory of itinerant ferromagnetism. Phys. Rev. B 1977, 16, 4032. [Google Scholar] [CrossRef]
  35. Volovik, G.E. Linear momentum in ferromagnets. J. Phys. C 1987, 20, L83–L87. [Google Scholar] [CrossRef]
  36. Guslienko, K.Y. Gauge and emergent electromagnetic fields for moving magnetic topological solitons. Europhys. Lett. 2016, 113, 67002. [Google Scholar] [CrossRef]
  37. Faddeev, L.D. 40 Years in Mathematical Physics; World Scientific Publishing: Singapore, 1995; pp. 374–377. [Google Scholar]
  38. Zheng, F.; Kiselev, N.S.; Rybakov, F.N.; Yang, L.; Shi, W.; Blügel, S.; Dunin-Borkowski, R.E. Hopfion rings in a cubic chiral magnet. Nature 2023, 623, 718–723. [Google Scholar] [CrossRef]
  39. Yu, X.Z.; Kanazawa, N.; Onose, Y.; Kimoto, K.; Zhang, W.Z.; Ishiwata, S.; Matsui, Y.; Tokura, Y. Near room-temperature formation of a skyrmion crystal in thin-flms of the helimagnet FeGe. Nat. Mater. 2011, 10, 106–109. [Google Scholar] [CrossRef]
  40. Kundu, A.; Rybakov, Y.P. Closed-vortex-type solitons with Hopf index. J. Phys. A Math. Gen. 1982, 15, 269–275. [Google Scholar] [CrossRef]
  41. Komineas, S.; Papanicolaou, N. Topology and dynamics in ferromagnetic media. Phys. D Nonlinear Phenom. 1996, 99, 81–107. [Google Scholar] [CrossRef]
  42. Arfken, G.B.; Weber, H.J. Chapter 1. Vector analysis. In Mathematical Methods for Physicists, 6th ed.; Academic Press: San Diego, CA, USA, 2005; pp. 1–95. [Google Scholar]
  43. Khodzhaev, Z.; Turgut, E. Hopfion dynamics in chiral magnets. J. Phys. Condens. Matter 2022, 34, 225805. [Google Scholar] [CrossRef]
  44. Hobart, R.H. On the instability of a class of unitary field models. Proc. Phys. Soc. 1963, 82, 201–203. [Google Scholar] [CrossRef]
  45. Derrick, G.H. Comments on nonlinear wave equations as models for elementary particles. J. Math. Phys. 1964, 5, 1252–1254. [Google Scholar] [CrossRef]
  46. Dzyaloshinskii, I.E. Theory of helicoidal structures in antiferromagnets. Sov. Phys. JETP 1964, 19, 960–971. [Google Scholar]
  47. Bogdanov, A.N.; Yablonskii, D.A. Thermodynamically stable “vortices” in magnetically ordered crystals. The mixed state of magnets. Sov. Phys. JETP 1989, 68, 101–103. [Google Scholar]
  48. Skyrme, T.H.R. A non-linear field theory. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1961, 260, 127–138. [Google Scholar] [CrossRef]
  49. Bogdanov, A. New localized solutions of the nonlinear field equations. JETP Lett. 1995, 62, 247–251. [Google Scholar]
  50. Guslienko, K.Y. Skyrmion state stability in magnetic nanodots with perpendicular anisotropy. IEEE Magn. Lett. 2015, 6, 4000104. [Google Scholar] [CrossRef]
  51. Berganza, E.; Jaafar, M.; Fernandez-Roldan, J.A.; Goiriena-Goikoetxea, M.; Pablo-Navarro, J.; García-Arribas, A.; Guslienko, K.; Magén, C.; De Teresa, J.M.; Chubykalo-Fesenko, O.; et al. Half-hedgehog spin textures in sub-100 nm soft magnetic nanodots. Nanoscale 2020, 12, 18646–18653. [Google Scholar] [CrossRef]
  52. Berganza, E.; Fernandez-Roldan, J.A.; Jaafar, M.; Asenjo, A.; Guslienko, K.; Chubykalo-Fesenko, O. 3D quasi-skyrmions in thick cylindrical and dome-shape soft nanodots. Sci. Rep. 2022, 12, 3426. [Google Scholar] [CrossRef]
  53. Guslienko, K. 3D Magnetization Textures: Toroidal Magnetic Hopfion Stability in Cylindrical Samples. Nanomaterials 2024, 14, 125. [Google Scholar] [CrossRef]
  54. Sobucki, K.; Krawczyk, M.; Tartakivska, O.; Graczyk, P. Magnon spectrum of Bloch hopfion beyond ferromagnetic resonance. APL Mater. 2022, 10, 091103. [Google Scholar] [CrossRef]
  55. Bo, L.; Ji, L.; Hu, C.; Zhao, R.; Li, Y.; Zhang, J.; Zhang, X. Spin excitation spectrum of a magnetic hopfion. Appl. Phys. Lett. 2021, 119, 212408. [Google Scholar] [CrossRef]
  56. Raftrey, D.; Fischer, P. Field-Driven Dynamics of Magnetic Hopfions. Phys. Rev. Lett. 2021, 127, 257201. [Google Scholar] [CrossRef] [PubMed]
  57. Liu, Y.; Watanabe, H.; Nagaosa, N. Emergent magnetomultipoles and nonlinear responses of a magnetic hopfion. Phys. Rev. Lett. 2022, 129, 267201. [Google Scholar] [CrossRef]
  58. Guslienko, K.Y. Magnetic vortex core string gyrotropic oscillations in thick cylindrical dots. Magnetism 2022, 2, 239–250. [Google Scholar] [CrossRef]
  59. Guslienko, K.Y.; Gareeva, Z.V. Magnetic skyrmion low frequency dynamics in thin circular dots. J. Magn. Magn. Mater. 2017, 442, 176–182. [Google Scholar] [CrossRef]
  60. Pershoguba, S.S.; Andreoli, D.; Zang, J. Electronic scattering off a magnetic hopfion. Phys. Rev. B 2021, 104, 075102. [Google Scholar] [CrossRef]
  61. Popadiuk, D.; Tartakovskaya, E.; Krawczyk, M.; Guslienko, K. Emergent magnetic field and non-zero gyrovector of the toroidal magnetic Hopfion. Phys. Status Solidi (RRL)–Rapid Res. Lett. 2023, 17, 2300131. [Google Scholar] [CrossRef]
  62. Yu, X.; Liu, Y.; Iakoubovskii, K.V.; Nakajima, K.; Kanazawa, N.; Nagaosa, N.; Tokura, Y. Realization and Current-Driven Dynamics of Fractional Hopfions and Their Ensembles in a Helimagnet FeGe. Adv. Mater. 2023, 35, 2210646. [Google Scholar] [CrossRef]
  63. Knapman, R.; Tausendpfund, T.; Díaz, S.A.; Everschor-Sitte, K. Spacetime magnetic hopfions from internal excitations and braiding of skyrmions. Commun. Phys. 2024, 7, 151. [Google Scholar] [CrossRef]
  64. Kobayashi, M.; Nitta, M. Winding Hopfions on R2 × S1. Nucl. Phys. B 2016, 876, 605–618. [Google Scholar] [CrossRef]
  65. Saji, C.; Troncoso, R.E.; Carvalho-Santos, V.L.; Altbir, D.; Nunez, A.S. Hopfion-Driven Magnonic Hall Effect and Magnonic Focusing. Phys. Rev. Lett. 2023, 131, 166702. [Google Scholar] [CrossRef]
Figure 1. The profiles of the out-of-plane field component, n z ρ , of the toroidal hopfions in the plane z = 0. The green solid line corresponds to m , n = 1 , 1 , the red dashed line— m , n = 1 , 2 , and the blue dotted line— m , n = 2 , 1 . The hopfion polarity p = + 1 .
Figure 1. The profiles of the out-of-plane field component, n z ρ , of the toroidal hopfions in the plane z = 0. The green solid line corresponds to m , n = 1 , 1 , the red dashed line— m , n = 1 , 2 , and the blue dotted line— m , n = 2 , 1 . The hopfion polarity p = + 1 .
Magnetism 04 00025 g001
Figure 2. The in-plane magnetization components m x , m y of the toroidal hopfions with the vorticities m , n = 1 , 1 in the z = 0 plane calculated by Equation (2): (A) the Bloch hopfion, φ 0 = 0 ; (B) the Neel hopfion, φ 0 = π / 2 .
Figure 2. The in-plane magnetization components m x , m y of the toroidal hopfions with the vorticities m , n = 1 , 1 in the z = 0 plane calculated by Equation (2): (A) the Bloch hopfion, φ 0 = 0 ; (B) the Neel hopfion, φ 0 = π / 2 .
Magnetism 04 00025 g002aMagnetism 04 00025 g002b
Figure 3. The magnetization components m x , m z of the toroidal Bloch hopfion φ 0 = 0   with the vorticities m , n = 1 , 1 in the x z plane calculated by Equation (2).
Figure 3. The magnetization components m x , m z of the toroidal Bloch hopfion φ 0 = 0   with the vorticities m , n = 1 , 1 in the x z plane calculated by Equation (2).
Magnetism 04 00025 g003
Figure 4. The magnetization components m x , m z of the winding hopfions in the x z plane calculated by Equation (7): (A) φ 0 = 0 ; (B) φ 0 = π / 2   . The azimuthal vorticity n = 1 . The parameter P = 1 . The vertical coordinate z 0 , 1 is plotted in units of the film thickness L . The horizontal coordinate x varies from 3 to 3 in units of the skyrmion radius ρ c .
Figure 4. The magnetization components m x , m z of the winding hopfions in the x z plane calculated by Equation (7): (A) φ 0 = 0 ; (B) φ 0 = π / 2   . The azimuthal vorticity n = 1 . The parameter P = 1 . The vertical coordinate z 0 , 1 is plotted in units of the film thickness L . The horizontal coordinate x varies from 3 to 3 in units of the skyrmion radius ρ c .
Magnetism 04 00025 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guslienko, K. Magnetic Hopfions: A Review. Magnetism 2024, 4, 383-399. https://doi.org/10.3390/magnetism4040025

AMA Style

Guslienko K. Magnetic Hopfions: A Review. Magnetism. 2024; 4(4):383-399. https://doi.org/10.3390/magnetism4040025

Chicago/Turabian Style

Guslienko, Konstantin. 2024. "Magnetic Hopfions: A Review" Magnetism 4, no. 4: 383-399. https://doi.org/10.3390/magnetism4040025

APA Style

Guslienko, K. (2024). Magnetic Hopfions: A Review. Magnetism, 4(4), 383-399. https://doi.org/10.3390/magnetism4040025

Article Metrics

Back to TopTop