New Types of Three-Dimensional Vortices in the Heisenberg Model

The Heisenberg model for an isotropic magnet is considered in this work. A substitution that reduces the corresponding equations to equations with a simpler geometric interpretation is applied. One of the solutions of the latter describes a new magnetic structure comprised of two straight intersecting vortex filaments, which change the topological charge after intersection.


INTRODUCTION
Among the various structures in magnetic media investigated in the last decade, there has been a particular interest in vortex structures. The important role of vortices in the description of magnetic topological phase transitions has been noted by many researchers (see, e.g., [1][2][3][4][5]). Such structures are not only of academic interest. Vortices in chiral magnets (chiral skyrmions), predicted two decades ago [6,7], can find important applications in running memory [8][9][10]. To date, two-dimensional vortices in ferro-and antiferromagnets have been studied in sufficient detail [6,7,[11][12][13][14][15]. Note that, despite the years-long study of vortices in "large" (bulk) samples, three-dimensional vortex structures have been insufficiently studied theoretically and experimentally. However, a new experimental technique: off-axis electron holography, has recently appeared [16], which makes it possible to study three-dimensional localized structures. It was revealed that the two-dimensional chiral skyrmions predicted in [5,6] have a three-dimensional structure [17]. Of the three-dimensional structures, only magnonic drops [10], magnetic "hedgehogs", and spiral structures [18] have been theoretically described to date. Therefore, the study of three-dimensional structures by analytical methods is of great importance.
Such methods are generally applicable to models with a high degree of symmetry. Although the applicability of analytical methods is limited, their scientific value is not in doubt. They allow one to thoroughly investigate the structure of the nuclei of nonlinear structures and qualitatively take into account the influence of other interactions not included in the initial model. In addition, exact solutions are "seed" functions in computer simulations of energy density minimization, including magnetic fields, anisotropy energy, etc. Vortex structures in a ferromagnet are formed mainly by exchange interaction (Heisenberg model), and the objective of this work is to find threedimensional vortices in this model.
We discuss new types of three-dimensional vortex structures in the Heisenberg model in this paper. The article is planned as follows. In the first section, in the equations of an isotropic magnet, we use a substitution that reduces the model to the pendulum equation and four equations with a simple geometric interpretation. The latter are reduced to two equations for a complex function We find the general solution of these equations, which depends on an arbitrary function. The simplest solution of these equations is discussed in Section 2. It is shown that it describes a new magnetic structure, which is comprised of two intersecting rectilinear vortex filaments, which change the topological charge after the intersection. At the end of the section, the experimental implementation of the structures found is discussed.

SIMPLIFYING THE EQUATIONS
OF THE HEISENBERG MODEL The Heisenberg model is based on using the Hamiltonian H of the exchange interaction, where the summation is performed over all pairs {i, j} of different sites of the crystal accomodating ions with spins S i and S j and the constants J ij characterize the exchange interaction between these ions. Taking into account the isotropic exchange interaction between ( ) nearest neighbors in the continuum limit, the density of the Hamiltonian of a ferromagnet has a simple form: (1) For the ferromagnetic Heisenberg model, J > 0, and, for the antiferromagnetic model, J < 0 [19]. Note that decomposition (1) takes into account only quadratic terms in the spin gradients. In the general case, to describe even the exchange interaction in some systems, it is necessary to take into account in (1) terms of a higher order in including the biquadratic terms for the spin In addition, we further investigate the case of zero magnetic field, when the ground state of the Heisenberg model is a uniformly magnetized ferromagnet without any structures.
Then, the three-dimensional Landau-Lifshitz equations can be written in the form (2) where n is the unit vector of magnetization S = M 0 n and M 0 is the spontaneous magnetization. Hereinafter, is the three-dimensional Laplace operator. The vector n is parameterized by the fields and In these variables, Eq. (2) turns into a system of nonlinear differential equations: Equation (2) is invariant with respect to the group of spin and spatial rotations. This symmetry makes it possible to find a wide class of exact solutions. The analytical solution of Eq. (2) is possible only in certain classes of solutions. To single out one of them, it is necessary to generalize the procedure proposed in [18] and assume that the field locally depends on the auxiliary field: Then, by direct calculations, it is easy to verify that the equations As we will see below, such substitutions lead to a wide class of exact solutions of the nonintegrable model (2).
Note that Eqs. (5) have a simple geometric interpretation. They describe two surfaces: and intersecting at each point, the normals to which, and , respectively, have unit length and are orthogonal and divergence-free: Let us turn to solving Eqs. (4) and (5). We introduce a complex field and write (5) as a system of two equations for the field S: This system has a remarkable property of invariance with respect to arbitrary changes of the field S, which we will use below. It is easy to verify that, if the field is a solution of this system, then an arbitrary function is also its solution. To solve system (6) and (7), we use a procedure proposed in [18]. We introduce a new complex field as (8) Here, the subscript denotes differentiation with respect to the corresponding variable: etc. The substitution of (8) into (6) determines The consistency condition for (8) and (9) gives a nonlinear equation for the field T: From the theory of first-order partial differential equations, it follows that the field is determined by the implicit equation (10) where and are integrals of the characteristic system of equations for the coordinates x z S x y z T x y z S x y z

t T x y z t t x t T x y z T x y z t t
The integrals have the form (11) From this wide class of solutions (10) and (11), we choose one with the function in the form Then, the field T satisfies the equation and, therefore, It turns out that the field (12) satisfies not only Eq. (6), but also Eq. (7) For , expression (14) is simplified to (15) Further on, we analyze formulas (13) and (15) of the new structure in an isotropic magnet.

ANALYSIS OF THE VORTEX STRUCTURE
Note first that substitution (4) and (5) and system (6) and (7) are a generalization of the instanton theory [20] in the two-dimensional case for an isotropic magnet. A wide class of solutions found in [20] is described by the simple formula with an analytic function F. In the two-dimensional case, Eqs. (6) and (7) ln , Fz z x y S = S(z), and the relationship with the field w, according to (15), is established by the simple formula Despite the apparent simplicity, the expression for the field (13) has a rich three-dimensional vortex structure. Previously investigated [6,7,11,12] twodimensional vortices comprised a straight filament with dependences and in a polar coordinate system Formula (13) describes two intersecting vortex filaments. One of them is a rectilinear vortex filament located on the axis and, in the polar coordinate system ( ) in the plane , has the form (16) with a topolgical charge equal to Here, is an arbitrary closed contour in the plane enclosing the vortex center. Integration is performed counterclockwise. Formula (16) has the simplest form near the vortex center ( ): and describes a vortex for and an antivortex for , with topological charges , respectively. From the boundary conditions and a numerical experiment, it follows that (13) describes two domain walls with nonexponential behavior at infinity. The centers of these domain walls (the locus of the points with the maximum value of the derivative) coincides, respectively, with the axes and Such vortex formation by the intersection of domain walls was studied in detail for the sine-Gordon model in [21,22].
The vortex structure is more evident in the distribution of the vector field This field in the planes and is shown in Figs. 1 and 2, respectively. It can be seen from Fig. 1 that, at a constant , the vector Z rotates from to counterclockwise if y < 0 and clockwise if y > 0; therefore, after crossing the domain boundaries, they change sign. Expression (13) is not defined at −π 2 π 2 = 0, x therefore, in the plane , it describes the second straight vortex filament. In the coordinate system , the vortex structure described by the expression is a vortex (if ) with a topological charge Q and an antivortex (if ) with the same charge, which is also formed by the intersection of domain boundaries.
The vortex structure is shown in Figs. 3 and 4, in which we clearly see jumps of the field by (Fig. 3, ) and (Fig. 4, ).
The structure of the field at in the standard polar system, depends on three spatial variables.
Although, at , the field tends to the ground state ( ), at , the azimuthal angle depends on the polar angle . t Therefore, the field n (2) has the following asymptotics: it is not singular at the origin and (17) The energy density of the Heisenberg model, under constraints (4) and (5) has a simple form and the total energy is proportional to the size of the system, L: This energy is significantly lower than the energy of a two-dimensional vortex [11,12] in an easy-plane ferromagnet, which is proportional to ( is the order of the lattice constant) per atomic layer.
The vortex structures predicted can nucleate in a cylinder with surface anisotropy, by thermal fluctuations, or an alternating magnetic field. In this case, the magnetization on the lateral surfaces takes the form (17). ACKNOWLEDGMENTS We are grateful to F.N. Rybakov for discussion and useful comments.

FUNDING
This work was carried out as part of state assignment of the Ministry of Education and Science of the Russian Federation (topic "Quant", no. AAAA-A18-118020190095-4).

OPEN ACCESS
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.