Magnetic Vortex Core String Gyrotropic Oscillations in Thick Cylindrical Dots

: The nonuniform magnetic vortex gyrotropic oscillations along the cylindrical dot thickness were calculated. A generalized Thiele equation was used for describing the vortex core motion including magnetostatic and exchange forces. The magnetostatic interaction was accounted for in a local form. This allowed reducing the Thiele equation of motion to the Schrödinger differential equation and analytically determining the spin eigenmode spatial proﬁles and eigenfrequencies using the Liouville–Green method for the high-frequency modes. The mapping of the Schrödinger equation to the Mathieu equation was used for the low-frequency gyrotropic mode. The lowest-frequency gyrotropic mode transformed to the dot faces localized mode, increasing the dot thickness. The vortex gyrotropic modes are described for a wide range of the dot thicknesses according to the concept of the turning points in the magnetostatic potential. This approach allows treating the vortex localized modes (turning points) and nonlocalized modes within a uniﬁed picture.


Introduction
There is a longstanding interest in the topological spin textures in small ferromagnetic particles and thin films, such as magnetic vortices [1,2], skyrmions [3][4][5], and hopfions [6,7]. The magnetic vortex is one of the simplest topologically nontrivial and stable spin textures in condensed matter physics. Since the vortex state of magnetization was discovered as the ground state of patterned soft magnetic dots [1], the dynamics of magnetic vortices attract considerable attention.
The vortex gyrotropic frequency and trajectory in circular magnetic dots were predicted by Guslienko et al. [8]. Being displaced from its equilibrium position in the dot center, the vortex core reveals sub-GHz frequency oscillations with a narrow linewidth. Then, existence of the vortex low-frequency gyrotropic mode dominated by the dipolar interaction was observed experimentally by different experimental techniques [9][10][11]. The experimental study of magnetic vortex dynamics was started in 2003 using time-resolved Kerr microscopy [9] and X-ray imaging [10]. Then, Novosad et al. [11] demonstrated that, upon being shifted from its equilibrium position by a microwave field, the vortex core experiences oscillations with the eigenfrequency proportional to the dot thickness/radius aspect ratio β = L/R. Later, it was shown that similar microwave-frequency oscillations with extremely narrow linewidth could be excited by DC spin polarized current [12,13]. Lastly, a large microwave generation from the current-driven magnetic vortex oscillators was observed in magnetic tunnel junctions [14]. These results allow considering the magnetic vortex gyrotropic motion as an important excitation mode of the spin torque nano-oscillators (STNO) in layered nanopillars and nanocontacts [15].
The magnetic vortex excitations in patterned films have been studied extensively for the last decades. However, in all experiments conducted in 2003-2013 the observed vortex gyrotropic modes had a uniform profile along the dot thickness accounting for the dot thickness being relatively small (typically, 20-30 nm). The nonuniform along-thickness vortex core mode was first simulated micromagnetically by Boust and Vukadinovic for small-radius circular dots [16]. The mode frequency decreases with dot thickness L increasing approximately as ∼ 1/L 2 . A model of the vortex core position inhomogeneous along the dot thickness was developed [17], which principally explained this behavior as a consequence of the exchange interaction in the out-of-plane z-direction. Then, the existence of higher-order vortex gyrotropic modes was demonstrated by micromagnetic simulations for thick (60-80 nm) submicron rectangular prisms [18]. These modes were described schematically as oscillations of the vortex core line with several nodes along the dot thickness (or oscillation modes of the asymmetric Bloch wall). Then, it was proven experimentally and confirmed by simulations and calculations that, in thick circular dots, these kinds of the vortex dynamical deformations, i.e., flexure oscillations of the vortex core string with n nodes along the dot thickness, can exist [19,20].
In this paper, a route from 2D to 3D magnetic soliton dynamics is considered conducting generalization of the 2D magnetic vortex strings to 3D topological spin textures. A third spatial dimension is added to 2D magnetic vortices, exploring the vortex gyrotropic modes in thick cylindrical dots and nanowires and considering the 3D vortex strings (tubes) oscillations. Specifically, the thickness-dependent magnetic vortex gyrotropic dynamics is calculated in the thick circular dots. A new analytic approach to the problem is developed, allowing the vortex eigenfrequencies and gyrotropic mode profiles to be calculated. The proposed analytical theory is based on the Thiele equation of motion of the soliton center and turning points of the vortex core motion along the dot thickness. The theory describes the whole vortex gyrotropic spectra in a circular dot, accounting for the exchange and magnetostatic energies.

Stability of the Magnetic Vortex State in Cylindrical Geometry
Let us consider a cylindrical dot (wire) of radius R and thickness (wire length) L made of soft magnetic material like Permalloy (NiFe alloy). Let us also assume that the dot (wire) radius R > 2L e , where L e = √ 2A/M s is the material exchange length, to allow the vortex state energy minimum (A is the exchange stiffness, and M s is the saturation magnetization). Then, following Usov et al. [21], the total dot magnetic energy, consisting of the exchange energy and magnetostatic energy of the face charges, can be minimized with respect to the vortex core radius, c = R c /R. The equilibrium vortex core radius at L > L e is approximately equal to R c = uL e (L/L e ) 1/3 , where u =0.685. We are interested to find a border of the vortex state stability increasing the dot thickness, at L > 2L e , when the vortex state gradually transforms to a quasi-uniform out-of-plane magnetized state (quasi-uniform state magnetized along the nanowire length if L L e , β = L/R 1). The simplest estimation of the critical thickness L c for such a transition is R c = R, i.e., L c (R)/L e = (R/uL e ) 3 . However, more careful analysis based on the vortex equilibrium magnetic energy and the vortex "half-skyrmion" ansatz [21] allows writing down a more exact value of the critical thickness (wire length) as L c (R)/L e = (R/uL e ) 3 / ln 2. For the typical material parameters L e = 14 nm and the nanodot (nanowire) radius R = 60 nm, the critical length is L c ≈ 2.36 µm. It is assumed here that the dot (wire) static magnetization is uniform along the thickness (length) coordinate. If any nonuniformity along the thickness is accounted for, then the critical value of L c can be smaller. An example of such nonuniformity is an "end-vortex" state, when the magnetization configuration is vortex-like near the nanowire ends (dot faces) and almost uniform (along the wire length) in the wire main part. Nevertheless, for the given dot radius R, there is a wide interval L e < L < L c (R) of the dot thickness (wire length) where the vortex state is stable. The minimal value of L c (R) for the given dot radius R is several hundred nanometers. The typical dot radius is 100-300 nm. Therefore, we can safely consider that the dots with the aspect ratio β = L/R = 0.2 to 2.0 are in the single vortex state and consider excitations over this ground (metastable) magnetization background. Such large values of the dot thickness L = 40 to 400 nm correspond to the extension of the vortex stability region above the value of L/L e = 3.5 calculated in [22].

The Generalized Thiele Equation of Motion of the Magnetic Soliton Center
Circular dots (wires) made of soft magnetic material have a vortex state of magnetization as the ground state for a certain dot radius R and thickness L. To calculate the vortex state excitation frequencies, the Landau-Lifshitz-Gilbert (LLG) equation of motion is used, which describes the evolution of magnetization M(r, t) in a ferromagnet in the presence of an effective magnetic field H and damping. It is written in the following form: where γ is the gyromagnetic ratio, α is the dimensionless damping constant, and M s is the saturation magnetization. The effective field H = −δE/δM is a combination of the external magnetic field, exchange field, the demagnetizing field, etc., where E is the total magnetic energy density. It was shown by Thiele [23] that, for a description of the dynamics of magnetic domains (solitons), the LLG equation can be rewritten in another form that allows simplifying the calculations. Thiele's approach is applicable to describe the dynamics of stable magnetization configurations (magnetic solitons) that can be characterized by a position of their center X(t) (a collective coordinate) that can vary with time. Let us use a cylindrical coordinate system with the cylinder axis directed along the z-axis. To consider the thickness-dependent vortex excitations, it is assumed that the magnetization as a function of the coordinates r = (ρ, z), ρ = (ρ, ϕ) and time can be written in the form M(r, t) = M(ρ, X(z, t)). The two-dimensional vector X(z, t) = (X(z, t), Y(z, t)) represents oscillations of the vortex core string in the dot xy-plane for a given value of z. Then, neglecting the damping, we can rewrite the LLG equation as the Thiele equation, where W = dVE is the total magnetic energy of the dot, and α, β = x, y. Introducing the unit vector of magnetization m = M/M s , the components of the gyrotensor density per unit dot thicknessĝ =Ĝ/L related to 2D topological charge of the spin texture can be defined as follows: We can apply Thiele's approach for the magnetic vortex motion in a circular dot (cylindrical wire). The vortex core position is defined in dimensionless units as s = X/R (s = s x + is y ).
The magnetization m has the following components: m x + im y = 2W/ 1 + WW and m z = p 1 − WW / 1 + WW , where W = W ζ, ζ is a complex function of the complex coordinate ζ = x + iy and its complex conjugate ζ, and p = m z (0) is the vortex core polarization [2]. Inside the vortex core the vortex configuration is described as a 2D topological soliton, i.e., W ζ, ζ = minimal value of ( ) for the given dot radius R is several hundred nanometers. The typical dot radius is 100-300 nm. Therefore, we can safely consider that the dots with the aspect ratio = = 0.2 to 2.0 ⁄ are in the single vortex state and consider excitations over this ground (metastable) magnetization background. Such large values of the dot thickness = 40 to 400 nm correspond to the extension of the vortex stability region above the value of = 3.5 ⁄ calculated in [22].

The Generalized Thiele Equation of Motion of the Magnetic Soliton Center
Circular dots (wires) made of soft magnetic material have a vortex state of magnetization as the ground state for a certain dot radius R and thickness L. To calculate the vortex state excitation frequencies, the Landau-Lifshitz-Gilbert (LLG) equation of motion is used, which describes the evolution of magnetization ( , ) in a ferromagnet in the presence of an effective magnetic field and damping. It is written in the following form: where γ is the gyromagnetic ratio, is the dimensionless damping constant, and is the saturation magnetization. The effective field = − ⁄ is a combination of the external magnetic field, exchange field, the demagnetizing field, etc., where is the total magnetic energy density. It was shown by Thiele [23] that, for a description of the dynamics of magnetic domains (solitons), the LLG equation can be rewritten in another form that allows simplifying the calculations. Thiele's approach is applicable to describe the dynamics of stable magnetization configurations (magnetic solitons) that can be characterized by a position of their center ( ) (a collective coordinate) that can vary with time. Let us use a cylindrical coordinate system with the cylinder axis directed along the z-axis. To consider the thickness-dependent vortex excitations, it is assumed that the magnetization as a function of the coordinates = ( , ), = ( , ) and time can be written in the form where = ∫ is the total magnetic energy of the dot, and , = , . Introducing the unit vector of magnetization = ⁄ , the components of the gyrotensor density per unit dot thickness ̑= ⁄ related to 2D topological charge of the spin texture can be defined as follows: We can apply Thiele's approach for the magnetic vortex motion in a circular dot (cylindrical wire). The vortex core position is defined in dimensionless units as = ⁄ ( = + ). The magnetization has the following components: is a complex function of the complex coordinate = + and its complex conjugate ̅ , and = (0) is the vortex core polarization [2]. Inside the vortex core the vortex configuration is described as a 2D topological soliton, i.e., ( , ) = ( ), where an analytic function ( ) satisfies the condition | ( )| ≤ 1. The vortex core position s in the xy-plane corresponds to a zero of the function ( ). Outside the vortex core region, where | ( )| > 1, the magnetization distribution is described by the function , = ( ) | ( )| ⁄ . The func-(ζ), where an analytic function minimal value of ( ) for the given dot radius R is several hundred nanometers. T typical dot radius is 100-300 nm. Therefore, we can safely consider that the dots with aspect ratio = = 0.2 to 2.0 ⁄ are in the single vortex state and consider excitati over this ground (metastable) magnetization background. Such large values of the thickness = 40 to 400 nm correspond to the extension of the vortex stability reg above the value of = 3.5 ⁄ calculated in [22].

The Generalized Thiele Equation of Motion of the Magnetic Soliton Center
Circular dots (wires) made of soft magnetic material have a vortex state of magn ization as the ground state for a certain dot radius R and thickness L. To calculate vortex state excitation frequencies, the Landau-Lifshitz-Gilbert (LLG) equation of m tion is used, which describes the evolution of magnetization ( , ) in a ferromagne the presence of an effective magnetic field and damping. It is written in the follow form: where γ is the gyromagnetic ratio, is the dimensionless damping constant, and the saturation magnetization. The effective field = − ⁄ is a combination of external magnetic field, exchange field, the demagnetizing field, etc., where is the to magnetic energy density. It was shown by Thiele [23] that, for a description of the namics of magnetic domains (solitons), the LLG equation can be rewritten in anot form that allows simplifying the calculations. Thiele's approach is applicable to descr the dynamics of stable magnetization configurations (magnetic solitons) that can characterized by a position of their center ( ) (a collective coordinate) that can v with time. Let us use a cylindrical coordinate system with the cylinder axis directed alo the z-axis. To consider the thickness-dependent vortex excitations, it is assumed that magnetization as a function of the coordinates = ( , ), = ( , ) and time can written in the form ( , ) = , ( , ) . The two-dimensional vector ( , ( , ), ( , ) represents oscillations of the vortex core string in the dot xy-plane fo given value of z. Then, neglecting the damping, we can rewrite the LLG equation as Thiele equation, where = ∫ is the total magnetic energy of the dot, and , = , . Introduc the unit vector of magnetization = ⁄ , the components of the gyrotensor den per unit dot thickness ̑= ⁄ related to 2D topological charge of the spin texture be defined as follows: We can apply Thiele's approach for the magnetic vortex motion in a circular (cylindrical wire). The vortex core position is defined in dimensionless units as = ( = + ). The magnetization has the following components: is a complex fu tion of the complex coordinate = + and its complex conjugate ̅ , and = is the vortex core polarization [2]. Inside the vortex core the vortex configuration is scribed as a 2D topological soliton, i.e., ( , ) = ( ), where an analytic function satisfies the condition | ( )| ≤ 1. The vortex core position s in the xy-plane correspo to a zero of the function ( ).
Outside the vortex core region, where | ( )| > 1, magnetization distribution is described by the function , = ( ) | ( )| ⁄ . The fu (ζ) satisfies the condition | e of ( ) for the given dot radius R is several hundred nanometers. The dius is 100-300 nm. Therefore, we can safely consider that the dots with the = = 0.2 to 2.0 ⁄ are in the single vortex state and consider excitations und (metastable) magnetization background. Such large values of the dot 40 to 400 nm correspond to the extension of the vortex stability region lue of = 3.5 ⁄ calculated in [22].

ralized Thiele Equation of Motion of the Magnetic Soliton Center
dots (wires) made of soft magnetic material have a vortex state of magnete ground state for a certain dot radius R and thickness L. To calculate the excitation frequencies, the Landau-Lifshitz-Gilbert (LLG) equation of mowhich describes the evolution of magnetization ( , ) in a ferromagnet in of an effective magnetic field and damping. It is written in the following e gyromagnetic ratio, is the dimensionless damping constant, and is n magnetization. The effective field = − ⁄ is a combination of the netic field, exchange field, the demagnetizing field, etc., where is the total ergy density. It was shown by Thiele [23] that, for a description of the dyagnetic domains (solitons), the LLG equation can be rewritten in another ows simplifying the calculations. Thiele's approach is applicable to describe s of stable magnetization configurations (magnetic solitons) that can be by a position of their center ( ) (a collective coordinate) that can vary t us use a cylindrical coordinate system with the cylinder axis directed along consider the thickness-dependent vortex excitations, it is assumed that the n as a function of the coordinates = ( , ), = ( , ) and time can be the form ( , ) = , ( , ) . The two-dimensional vector ( , ) = ) represents oscillations of the vortex core string in the dot xy-plane for a of z. Then, neglecting the damping, we can rewrite the LLG equation as the on, ∫ is the total magnetic energy of the dot, and , = , . Introducing or of magnetization = ⁄ , the components of the gyrotensor density thickness ̑= ⁄ related to 2D topological charge of the spin texture can follows: apply Thiele's approach for the magnetic vortex motion in a circular dot ire). The vortex core position is defined in dimensionless units as = ⁄ ). The magnetization has the following components: is a complex funcomplex coordinate = + and its complex conjugate ̅ , and = (0) core polarization [2]. Inside the vortex core the vortex configuration is de-2D topological soliton, i.e., ( , ) = ( ), where an analytic function ( ) condition | ( )| ≤ 1. The vortex core position s in the xy-plane corresponds the function ( ).
Outside the vortex core region, where | ( )| > 1, the n distribution is described by the function , = ( ) | ( )| ⁄ . The func- The vortex core position s in the xy-plane corresponds to a zero of the function 3 e of ( ) for the given dot radius R is several hundred nanometers. The dius is 100-300 nm. Therefore, we can safely consider that the dots with the = = 0.2 to 2.0 ⁄ are in the single vortex state and consider excitations und (metastable) magnetization background. Such large values of the dot 40 to 400 nm correspond to the extension of the vortex stability region ue of = 3.5 ⁄ calculated in [22].

alized Thiele Equation of Motion of the Magnetic Soliton Center
dots (wires) made of soft magnetic material have a vortex state of magnetground state for a certain dot radius R and thickness L. To calculate the xcitation frequencies, the Landau-Lifshitz-Gilbert (LLG) equation of mohich describes the evolution of magnetization ( , ) in a ferromagnet in of an effective magnetic field and damping. It is written in the following e gyromagnetic ratio, is the dimensionless damping constant, and is magnetization. The effective field = − ⁄ is a combination of the etic field, exchange field, the demagnetizing field, etc., where is the total rgy density. It was shown by Thiele [23] that, for a description of the dygnetic domains (solitons), the LLG equation can be rewritten in another ws simplifying the calculations. Thiele's approach is applicable to describe of stable magnetization configurations (magnetic solitons) that can be by a position of their center ( ) (a collective coordinate) that can vary t us use a cylindrical coordinate system with the cylinder axis directed along consider the thickness-dependent vortex excitations, it is assumed that the as a function of the coordinates = ( , ), = ( , ) and time can be he form ( , ) = , ( , ) . The two-dimensional vector ( , ) = ) represents oscillations of the vortex core string in the dot xy-plane for a f z. Then, neglecting the damping, we can rewrite the LLG equation as the n, is the total magnetic energy of the dot, and , = , . Introducing r of magnetization = ⁄ , the components of the gyrotensor density hickness ̑= ⁄ related to 2D topological charge of the spin texture can follows: apply Thiele's approach for the magnetic vortex motion in a circular dot ire). The vortex core position is defined in dimensionless units as = ⁄ ). The magnetization has the following components: is a complex funcmplex coordinate = + and its complex conjugate ̅ , and = (0) core polarization [2]. Inside the vortex core the vortex configuration is de-2D topological soliton, i.e., ( , ) = ( ), where an analytic function ( ) ondition | ( )| ≤ 1. The vortex core position s in the xy-plane corresponds the function ( ).
Outside the vortex core region, where | ( )| > 1, the distribution is described by the function , = ( ) | ( )| ⁄ . The func-(ζ). Outside the vortex core region, where | Magnetism 2022, 2, FOR PEER REVIEW 3 minimal value of ( ) for the given dot radius R is several hundred nanometers. The typical dot radius is 100-300 nm. Therefore, we can safely consider that the dots with the aspect ratio = = 0.2 to 2.0 ⁄ are in the single vortex state and consider excitations over this ground (metastable) magnetization background. Such large values of the dot thickness = 40 to 400 nm correspond to the extension of the vortex stability region above the value of = 3.5 ⁄ calculated in [22].

The Generalized Thiele Equation of Motion of the Magnetic Soliton Center
Circular dots (wires) made of soft magnetic material have a vortex state of magnetization as the ground state for a certain dot radius R and thickness L. To calculate the vortex state excitation frequencies, the Landau-Lifshitz-Gilbert (LLG) equation of motion is used, which describes the evolution of magnetization ( , ) in a ferromagnet in the presence of an effective magnetic field and damping. It is written in the following form: where γ is the gyromagnetic ratio, is the dimensionless damping constant, and is the saturation magnetization. The effective field = − ⁄ is a combination of the external magnetic field, exchange field, the demagnetizing field, etc., where is the total magnetic energy density. It was shown by Thiele [23] that, for a description of the dynamics of magnetic domains (solitons), the LLG equation can be rewritten in another form that allows simplifying the calculations. Thiele's approach is applicable to describe the dynamics of stable magnetization configurations (magnetic solitons) that can be characterized by a position of their center ( ) (a collective coordinate) that can vary with time. Let us use a cylindrical coordinate system with the cylinder axis directed along the z-axis. To consider the thickness-dependent vortex excitations, it is assumed that the magnetization as a function of the coordinates = ( , ), = ( , ) and time can be written in the form ( , ) = , ( , ) . The two-dimensional vector ( , ) = ( , ), ( , ) represents oscillations of the vortex core string in the dot xy-plane for a given value of z. Then, neglecting the damping, we can rewrite the LLG equation as the Thiele equation, where = ∫ is the total magnetic energy of the dot, and , = , . Introducing the unit vector of magnetization = ⁄ , the components of the gyrotensor density per unit dot thickness ̑= ⁄ related to 2D topological charge of the spin texture can be defined as follows: We can apply Thiele's approach for the magnetic vortex motion in a circular dot (cylindrical wire). The vortex core position is defined in dimensionless units as = ⁄ ( = + ). The magnetization has the following components: is a complex function of the complex coordinate = + and its complex conjugate ̅ , and = (0) is the vortex core polarization [2]. Inside the vortex core the vortex configuration is described as a 2D topological soliton, i.e., ( , ) = ( ), where an analytic function ( ) satisfies the condition | ( )| ≤ 1. The vortex core position s in the xy-plane corresponds to a zero of the function ( ).
Outside the vortex core region, where | ( )| > 1, the magnetization distribution is described by the function , = ( ) | ( )| ⁄ . The func-(ζ)| > 1, the magnetization distribution is described by the function W ζ, ζ = ism 2022, 2, FOR PEER REVIEW 3 minimal value of ( ) for the given dot radius R is several hundred nanometers. The typical dot radius is 100-300 nm. Therefore, we can safely consider that the dots with the aspect ratio = = 0.2 to 2.0 ⁄ are in the single vortex state and consider excitations over this ground (metastable) magnetization background. Such large values of the dot thickness = 40 to 400 nm correspond to the extension of the vortex stability region above the value of = 3.5 ⁄ calculated in [22].

The Generalized Thiele Equation of Motion of the Magnetic Soliton Center
Circular dots (wires) made of soft magnetic material have a vortex state of magnetization as the ground state for a certain dot radius R and thickness L. To calculate the vortex state excitation frequencies, the Landau-Lifshitz-Gilbert (LLG) equation of motion is used, which describes the evolution of magnetization ( , ) in a ferromagnet in the presence of an effective magnetic field and damping. It is written in the following form: where γ is the gyromagnetic ratio, is the dimensionless damping constant, and is the saturation magnetization. The effective field = − ⁄ is a combination of the external magnetic field, exchange field, the demagnetizing field, etc., where is the total magnetic energy density. It was shown by Thiele [23] that, for a description of the dynamics of magnetic domains (solitons), the LLG equation can be rewritten in another form that allows simplifying the calculations. Thiele's approach is applicable to describe the dynamics of stable magnetization configurations (magnetic solitons) that can be characterized by a position of their center ( ) (a collective coordinate) that can vary with time. Let us use a cylindrical coordinate system with the cylinder axis directed along the z-axis. To consider the thickness-dependent vortex excitations, it is assumed that the magnetization as where = ∫ is the total magnetic energy of the dot, and , = , . Introducing the unit vector of magnetization = ⁄ , the components of the gyrotensor density per unit dot thickness ̑= ⁄ related to 2D topological charge of the spin texture can be defined as follows: We can apply Thiele's approach for the magnetic vortex motion in a circular dot (cylindrical wire). The vortex core position is defined in dimensionless units as = ⁄ ( = + ). The magnetization has the following components: is a complex function of the complex coordinate = + and its complex conjugate ̅ , and = (0) is the vortex core polarization [2]. Inside the vortex core the vortex configuration is described as a 2D topological soliton, i.e., ( , ) = ( ), where an analytic function ( ) satisfies the condition | ( )| ≤ 1. The vortex core position s in the xy-plane corresponds to a zero of the function ( ).
Outside the vortex core region, where | ( )| > 1, the magnetization distribution is described by the function , = ( ) | ( )| ⁄ . The func- minimal value of ( ) for the given dot radius R is several hundred nanometers. The typical dot radius is 100-300 nm. Therefore, we can safely consider that the dots with the aspect ratio = = 0.2 to 2.0 ⁄ are in the single vortex state and consider excitations over this ground (metastable) magnetization background. Such large values of the dot thickness = 40 to 400 nm correspond to the extension of the vortex stability region above the value of = 3.5 ⁄ calculated in [22].

The Generalized Thiele Equation of Motion of the Magnetic Soliton Center
Circular dots (wires) made of soft magnetic material have a vortex state of magnetization as the ground state for a certain dot radius R and thickness L. To calculate the vortex state excitation frequencies, the Landau-Lifshitz-Gilbert (LLG) equation of motion is used, which describes the evolution of magnetization ( , ) in a ferromagnet in the presence of an effective magnetic field and damping. It is written in the following form: where γ is the gyromagnetic ratio, is the dimensionless damping constant, and is the saturation magnetization. The effective field = − ⁄ is a combination of the external magnetic field, exchange field, the demagnetizing field, etc., where is the total magnetic energy density. It was shown by Thiele [23] that, for a description of the dynamics of magnetic domains (solitons), the LLG equation can be rewritten in another form that allows simplifying the calculations. Thiele's approach is applicable to describe the dynamics of stable magnetization configurations (magnetic solitons) that can be characterized by a position of their center ( ) (a collective coordinate) that can vary with time. Let us use a cylindrical coordinate system with the cylinder axis directed along the z-axis. To consider the thickness-dependent vortex excitations, it is assumed that the magnetization as where = ∫ is the total magnetic energy of the dot, and , = , . Introducing the unit vector of magnetization = ⁄ , the components of the gyrotensor density per unit dot thickness ̑= ⁄ related to 2D topological charge of the spin texture can be defined as follows: We can apply Thiele's approach for the magnetic vortex motion in a circular dot (cylindrical wire). The vortex core position is defined in dimensionless units as = ⁄ ( = + ). The magnetization has the following components: is a complex function of the complex coordinate = + and its complex conjugate ̅ , and = (0) is the vortex core polarization [2]. Inside the vortex core the vortex configuration is described as a 2D topological soliton, i.e., ( , ) = ( ), where an analytic function ( ) satisfies the condition | ( )| ≤ 1. The vortex core position s in the xy-plane corresponds to a zero of the function ( ).
Outside the vortex core region, where | ( )| > 1, the magnetization distribution is described by the function , = ( ) | ( )| ⁄ . The func-(ζ)|. The function w ζ, ζ corresponds to the description of magnetic vortex as a nonlocalized topological soliton. The magnetostatic energy plays an essential role for the magnetic vortices in restricted geometries. Therefore, for describing the vortex dynamics, a two-vortex model (no side surface charges induced in the course of motion) is used with the analytical function Magnetism 2022, 2, FOR PEER REVIEW minimal value of ( ) for the given dot radius R is several hundred nanom typical dot radius is 100-300 nm. Therefore, we can safely consider that the do aspect ratio = = 0.2 to 2.0 ⁄ are in the single vortex state and consider over this ground (metastable) magnetization background. Such large values thickness = 40 to 400 nm correspond to the extension of the vortex stab above the value of = 3.5 ⁄ calculated in [22].

The Generalized Thiele Equation of Motion of the Magnetic Soliton Center
Circular dots (wires) made of soft magnetic material have a vortex state ization as the ground state for a certain dot radius R and thickness L. To ca vortex state excitation frequencies, the Landau-Lifshitz-Gilbert (LLG) equa tion is used, which describes the evolution of magnetization ( , ) in a ferr the presence of an effective magnetic field and damping. It is written in th form: where γ is the gyromagnetic ratio, is the dimensionless damping constant the saturation magnetization. The effective field = − ⁄ is a combina external magnetic field, exchange field, the demagnetizing field, etc., where magnetic energy density. It was shown by Thiele [23] that, for a description namics of magnetic domains (solitons), the LLG equation can be rewritten form that allows simplifying the calculations. Thiele's approach is applicable the dynamics of stable magnetization configurations (magnetic solitons) t characterized by a position of their center ( ) We can apply Thiele's approach for the magnetic vortex motion in a c (cylindrical wire). The vortex core position is defined in dimensionless units ( = + ). The magnetization has the following components: is a com tion of the complex coordinate = + and its complex conjugate ̅ , and is the vortex core polarization [2]. Inside the vortex core the vortex configur scribed as a 2D topological soliton, i.e., ( , ) = ( ), where an analytic fu satisfies the condition | ( )| ≤ 1. The vortex core position s in the xy-plane c to a zero of the function ( ).
Outside the vortex core region, where | ( magnetization distribution is described by the function , = ( ) | ( )| ⁄ (ζ) being written as 3 alue of ( ) for the given dot radius R is several hundred nanometers. The radius is 100-300 nm. Therefore, we can safely consider that the dots with the o = = 0.2 to 2.0 ⁄ are in the single vortex state and consider excitations round (metastable) magnetization background. Such large values of the dot = 40 to 400 nm correspond to the extension of the vortex stability region alue of = 3.5 ⁄ calculated in [22].
the gyromagnetic ratio, is the dimensionless damping constant, and is ion magnetization. The effective field = − ⁄ is a combination of the agnetic field, exchange field, the demagnetizing field, etc., where is the total nergy density. It was shown by Thiele [23] that, for a description of the dymagnetic domains (solitons), the LLG equation can be rewritten in another llows simplifying the calculations. Thiele's approach is applicable to describe ics of stable magnetization configurations (magnetic solitons) that can be ed by a position of their center ( ) (a collective coordinate) that can vary Let us use a cylindrical coordinate system with the cylinder axis directed along To consider the thickness-dependent vortex excitations, it is assumed that the tion as a function of the coordinates = ( , ), = ( , ) and time can be the form ( , ) = , ( , ) . The two-dimensional vector ( , ) = , ) represents oscillations of the vortex core string in the dot xy-plane for a e of z. Then, neglecting the damping, we can rewrite the LLG equation as the ation, ∫ is the total magnetic energy of the dot, and , = , . Introducing ctor of magnetization = ⁄ , the components of the gyrotensor density t thickness ̑= ⁄ related to 2D topological charge of the spin texture can as follows: n apply Thiele's approach for the magnetic vortex motion in a circular dot l wire). The vortex core position is defined in dimensionless units as = ⁄ ). The magnetization has the following components: is a complex funccomplex coordinate = + and its complex conjugate ̅ , and = (0) x core polarization [2]. Inside the vortex core the vortex configuration is dea 2D topological soliton, i.e., ( , ) = ( ), where an analytic function ( ) e condition | ( )| ≤ 1. The vortex core position s in the xy-plane corresponds of the function ( ).
Outside the vortex core region, where | ( )| > 1, the tion distribution is described by the function , = ( ) | ( )| ⁄ . The func- [2,8], where C = ±1 is the vortex chirality, and c = R c /R is the reduced vortex core radius. It is assumed that the dependence of magneti-zation on the thickness coordinate z is accounted for via the function s(z). The vortex core with the radius R c = uL e (L/L e ) 1/3 should be accounted for explicitly in the case of dots (nanowires) with radii R of several tens of nanometers.
The righthand-side term in Equation (2) is the variational derivative of the total magnetic energy W(X) of the dot with a displaced magnetic vortex centered in a point s = s x , s y taken with respect to the vortex core coordinate s = X/R. The total magnetic energy W consists of the magnetostatic, exchange, and Zeeman energies and can be calculated as a function of s. The magnetostatic energy consists of the energy of volume magnetic charges ρ(ρ, X) = −M s divm(ρ, X), energy of the side surface charges, and energy of the dot face (wire end) charges. The side surface charges are absent for the selected vortex model. The face charges related to the magnetization component m z contribute only to the boundary conditions at the dot faces z = ±L/2. Therefore, the magnetostatic energy is related with the volume magnetic charges and can be written in the following form: The exchange energy is W ex (X) = A dV ∑ α (∇m α (ρ, X)) 2 with A being the exchange stiffness parameter. Integration over the in-plane dot coordinates and neglecting the vortex core contribution (c 1) yields, in the quadratic approximation on the small parameter |s| 1, the nonlocal magnetostatic energy, where the magnetostatic kernel is κ(z, z ) = 8πM 2 s β dke −βk|z−z | I 2 (k), I(k) = 1 c dρρJ 1 (kρ), β = L/R, and V = πR 2 L is the dot volume. The thickness coordinate z is in units of L.
The exchange energy is where µ = 2M 2 s (L e /R) 2 is the in-plane and η = M 2 s (L e /L) 2 ln 1 c + 5 4 + c 2 is the out-ofplane exchange stiffness coefficient calculated within the model, and c = R c /R is the reduced vortex core radius.
The generalized Thiele equation of motion (Equation (2)) is where g = g xyẑ is the gyrovector per unit dot thickness, g = |g| = 2πM s /γ, W = V dzw/π, and W = W m + W ex . Assuming circular motion of the vortex core . s = ω × s with the frequency ω = |ω| (or making a Fourier transform s(z, t) = s(z, ω) exp(iωt) with respect to time t for the variable s = s x + is y ), the Thiele equation can be written as a two-component differential-integral equation containing nonlocal terms, for the vector s(z) = s x (z), s y (z) components, where s(z) = s(z, ω), and λ = 2ωM s /γ. Equation (8a) should be amended by the boundary conditions at the dot faces z = ±1/2. This equation can be written in the compact formΓs(z) = λs(z), where the linear integraldifferential operatorΓ is defined by the righthand side of Equation (8a). Equation (8a) allows finding eigenmodes of the vortex core string oscillations s n (z) and corresponding eigenfrequencies ω n , which are expressed via the eigenvalues λ n . It is convenient to write Equation (8a) in the dimensionless form, dividing both parts of the equation by 8πM 2 s . Then, the eigenvalue λ corresponds to the dimensionless frequency λ = ω/ω M , where ω M = γ4πM s . The value of ω M /2π is 30 GHz for a typical composition of Ni x Fe 1−x alloy with x ≈ 0.80 (M s = 800-860 kA/m).
It is convenient to introduce a complex circular variable s(z) = s x (z) + is y (z). The equation of motion for s(z) has the form The solution of Equation (8b) yields the complex eigenfunctions s n (z) and eigenfrequencies ω n of the vortex string oscillations. The absolute value u n (z) = |s n (z)| and the phase of the n-th eigenmode as s n (z) = u n (z)exp[iΦ n (z)] are introduced. It is determined from the symmetry of Equation (8b) that the eigenfunctions satisfy the condition s n (−z) = s n (z), which is equivalent to the conditions u n (−z) = u n (z) and Φ n (−z) = −Φ n (z). Then, the time dependence of the vortex displacement vector components s n (z, t) = s nx (z, t) + is ny (z, t) is given by the equations s nx (z, t) = u n (z)cos[ω n t + Φ n (z)], s ny (z, t) = u n (z)sin[ω n t + Φ n (z)]. Equation (8b) has pure real (even) solutions or pure imaginary (odd) solutions s n (z), if the phases Φ n (z) = 0 or Φ n (z) = ±π/2, respectively.
The aim is to find analytical solutions of Equation (8a), the eigenfrequencies λ n and eigenfunctions s n (z) of the operatorΓ numbered by the integer index n = 0, 1, 2, . . . The problem is that there is no theory of the linear integral-differential equations. Therefore, we reduce Equation (8a) to a local form, i.e., a differential equation for the function s n (z). The kernel κ(z, z ) has a sharp maximum at z ≈ z if the dot aspect ratio β is large enough. This allows substituting s(z ) in the magnetostatic term containing the kernel κ(z, z ) to s(z) and writing the Thiele equation in the form of the differential Schrödinger equation with an effective potential energy, where I(z) = 1/2 −1/2 dz κ(z, z ). The function I(z) can be calculated explicitly and has the form I(z, β) = 2 ∞ 0 dkk −1 I 2 (k)[1 − exp(−βk/2)cosh(βkz)]. I(z) has a smooth maximum at z = 0, in the dot middle plane, and minima at the dot faces, z = ±1/2. The potential I(z) ≈ const at a small dot aspect ratio β ∼ 0.1. However, upon increasing β, the function I(z, β) reveals more and more pronounced minima at z = ±1/2. Such minima result in the low-frequency (n = 0) eigenmode localization near the dot faces z = ±1/2 increasing the dot thickness or the dot aspect ratio β at the fixed dot radius R. The accuracy of the local approximation, Equation (9), increases with β increasing.
We can rewrite Equation (9) in the standard oscillator-like form for both vortex core displacement s components, s x , s y .
where a small thickness-independent coefficient µ is included in the eigenvalue λ.

The Lowest Frequency Gyrotropic Mode
Starting from Equation (10), it is possible to satisfy, for a given λ, the condition f (z, β) = 0, which corresponds to the classical turning points ±z * located somewhere near the dot faces. Other turning points z = ±1/2 appear naturally due to the restricted dot geometry. More detailed analysis shows that the turning points exist only for the lowest eigenvalue λ 0 . In other words, the lowest frequency eigenmode is localized between the turning points z * and 1 2 (or −z * and − 1 2 ) near the dot top/bottom surfaces. The mode frequency ω 0 (β) is within the interval I(1/2, β) < ω 0 (β)/ω M < I(0, β) shown in Figure 1. This interval is very narrow for small β ∼ 0.1 and rapidly increases with β increasing. Another estimation of the lowest mode gyrofrequency valid at small β << 1 can be obtained assuming the mode uniformity along the dot thickness, i.e., putting s(z) = const in Equation (10). This frequency was calculated in [8] as ω G (β)/ω M = Starting from Equation (10), it is possible to satisfy, for a given λ, the condition ( , ) = 0, which corresponds to the classical turning points ± * located somewhere near the dot faces. Other turning points = ± 1 2 ⁄ appear naturally due to the restricted dot geometry. More detailed analysis shows that the turning points exist only for the lowest eigenvalue . In other words, the lowest frequency eigenmode is localized between the turning points * and ½ (or − * and −½) near the dot top/bottom surfaces. The mode frequency ( ) is within the interval ℐ(1 2, ⁄ ) < ( ) < ⁄ ℐ(0, ) shown in Figure 1. This interval is very narrow for small ~0.1 and rapidly increases with increasing. Another estimation of the lowest mode gyrofrequency valid at small << 1 can be obtained assuming the mode uniformity along the dot thickness, i.e., putting s ( ) = in Equation (10). This frequency was calculated in [8] as Figure 1. The lowest vortex gyrofrequency (solid line 3) calculated using Equation (12). The dashed lines 1 and 4 describe the limiting frequencies corresponding to the frequency window ℐ(1 2, ⁄ ) < ( ) < ⁄ ℐ(0, ) of the mode localization between the turning points. The frequency (solid line 2) is the gyrotropic frequency calculated in [8] assuming the gyromode uniformity along the dot thickness extrapolated for large values of β. The reduced dot radius = = 10 ⁄ .
We approximate the magnetostatic potential ℐ( ) in Equation (10) as ℐ( ) = ℐ(0) ( ). This allows rewriting Equation (10) in the form of the Mathieu equation [24,25], We approximate the magnetostatic potential I(z) in Equation (10) as I(z) = I(0)cos(pz). This allows rewriting Equation (10) in the form of the Mathieu equation [24,25], where ξ = pz/2, a = 4λ/η p 2 , and q = 2I (0)/η p 2 . The fitting parameter p is defined as p = 2arccos(I (1/2)/I(0)) and is a function only of β. The solutions of the Mathieu equation are known. We are interested to use only the lowest eigenvalue a 0 of Equation (11), which is the increasing function of the parameter q. q is an increasing function on β (the dot radius R is fixed). For the values of q > 1, the eigenvalue a 0 of the Mathieu equation is defined approximately by the function a 0 (q) = 2q − 2 √ q + 1/4 − 1/32 √ q [24,25], where q(β, r) = 16πβ 2 r 2 I(0, β)/p 2 (β) ln(1/c) + 5 4 + c 2 , r = R/L e is the reduced dot radius, and c = 0.685β 1/3 /r 2/3 . For instance, if the value of r = 10 is used, then q > 1.82 at β > 0.2, and q(β = 1) = 56.64. The eigenvalue a 0 corresponds to the lowest gyrotropic mode eigenfrequency (see Figure 1), and the corresponding eigenfunction is the zero-order Mathieu even function ce 0 (ξ, q). For the large values of the parameter q >> 1, such an eigenfunction (not normalized) can be represented in the following form [24,25]: The even function s 0 (ξ) for the typical dot sizes is shown in Figure 2, where the normalized eigenfunctions N −1/2 0 s 0 (z), N 0 = dzs 2 0 (z) are plotted for different values of β. The degree of the lowest gyrotropic mode localization at the dot faces increases with increasing dot thickness. The lowest mode eigenfrequency given by Equation (12a) is asymptotically correct at the large values of the parameter q, i.e., for thick enough dots (q > 12.8 at β > 0.5, r = 10). The frequency ω 0 (β, r) is lower than the gyrotropic frequency ω G (β) calculated within uniform along thickness mode approximation in [8] at intermediate values of β and practically coincides with ω G (β) at large values of β. The approximation of ω 0 (β, r) by the Mathieu method (Equation (12a)) is not applicable at small values of q (small β < 0.2). It leads to underestimation of the lowest mode frequency in comparison with the homogeneous solution ω G (β).

High-Frequency Vortex Gyrotropic Modes
Let us consider the equation of motion (Equation (10)) for higher frequencies, when there are no turning points for the function ( ). If, additionally, we assume that th function is smooth enough, then the solution of Equation (10) can be written in the Liou ville-Green approximation [26] as

High-Frequency Vortex Gyrotropic Modes
Let us consider the equation of motion (Equation (10)) for higher frequencies, when there are no turning points for the function f (z). If, additionally, we assume that the function is smooth enough, then the solution of Equation (10) can be written in the Liouville-Green approximation [26] as where C + , C − are some constants, which can be defined by the boundary conditions at the dot faces, z = ±1/2. The function u(z) = ( f (z)) −1/4 has sense of the amplitude, and Φ(z) = z 0 dz f 1 2 (z) + Φ(0) has sense of the phase of the vortex string displacement s(z) = s x (z), s y (z) used in the complex form s = s x + is y . It is assumed that there is no surface magnetic anisotropy and no dipolar pining at the dot faces. Therefore, the natural exchange boundary conditions, (∂s/∂z) = 0 at z = ±1/2, are used. Equation (13) combined with these boundary conditions allows finding the vortex oscillation eigenfrequencies and explicit form of the eigenfunctions. It can be shown that solving the linear system of the equations for the coefficients C + , C − and the eigenfrequencies λ n can be found from the equation and the corresponding eigenfunctions are where ψ(λ) = Φ(1/2, λ) − Φ(−1/2, λ), f (z, λ), and σ = sign(cos(ψ)) are taken at λ = λ n . The function G(λ) = 2g(λ)/ g 2 (λ) − 1 rapidly decreases with the frequency λ = ω/ω M increasing. It is different from zero solely due to the magnetostatic potential I(z) nonuniformity along the dot thickness. For thin dots (small values of β 1), ∂ f (1/2, λ)/∂z ≈ 0 and G(λ) ≈ 0, tan(ψ(λ)) ≈ 0, ψ(λ) ≈ nπ for all eigenmode numbers n = 1, 2, . . . For the large values of β = 0.5 − 2, the condition ψ(λ) ≈ nπ is satisfied with a high degree of accuracy for all eigenmodes except the n = 1 mode. The phase difference ψ(β)/π for n = 1 mode decreases from 0.998 to 0.849 increasing β from 0.2 to 1.0. Apparently, the condition ψ(λ) ≈ nπ corresponds to the simple exchange eigenfunctions cos(nπz), sin(nπz) and corresponding exchange-dominated eigenfrequencies can be written in the form of Equation (4) of the paper by Ding et al. [19]. λ n = ω n /ω M .
The eigenfunctions s n (z) are even at σ = +1 or odd at σ = −1 with respect to the dot middle plane z = 0 according to the general approach. The non-normalized eigenfunctions can be written in the form s n (z) = f −1/4 (z)exp(iΦ(0))cos z 0 dz f 1 2 (z) for the even modes and s n (z) = f −1/4 (z)exp[i(Φ(0) + π/2)]sin z 0 dz f 1 2 (z) for the odd modes taken at λ = λ n . The normalized eigenfrequencies λ n = ω n /ω M are shown in Figure 3 and the normalized eigenfunctions are plotted in Figure 4 for typical dot parameters. The following normalization was used: N −1/2 n s n (z), N n = dz|s n (z)| 2 . The eigenfrequencies ω n (β) with n = 2, 3, . . . rapidly decrease with the dot thickness L increasing at small β as approximately ∼ 1/L 2 and approach some constant values at large β. The thickness dependence of the n = 1 mode ω 1 (β) is different. It reveals a smooth minimum at intermediate values of β. Another important peculiarity of the n = 1 vortex gyrotropic mode is that the mode frequency λ 1 = ω 1 (β)/ω M becomes close to the magnetostatic integral I(0, β), increasing the aspect ratio β up to approximately 1.0; therefore, the function f (z, λ) becomes very small, and the Liouville-Green approximation is no longer applicable. This corresponds to the regime of the frequency crossing between the lowest n = 0 (ω 0 (β)) and n = 1 (ω 1 (β)) vortex gyromodes and will be considered elsewhere.

The Role of the Vortex Mass in the Gyrotropic Dynamics
The relationship of the eigenvalue parameter λ with the vortex gyrotropic frequency ω is more complicated if other high-order time derivatives of the vortex core displacement are accounted for in the Thiele equation of motion (Equation (7)). Formally, this can be achieved assuming a velocity-dependent soliton ansatz ( , ) = ( , ( , ), ( , )⁄ ). It was shown in [27] that the mass term is, in general, nonlocal. If the vortex mass is accounted in a local approximation, then the equation of motion of the vortex center position is where = ⁄ is the vortex mass per unit dot thickness, and is the dimension- . The normalized eigenfunctions of the high-frequency vortex gyrotropic modes with the numbers n = 1, 2, and 3 calculated using the Liouville-Green approximation, Equation (15). The dot aspect ratio β = 1, and the reduced dot radius r = R/L e = 10.

The Role of the Vortex Mass in the Gyrotropic Dynamics
The relationship of the eigenvalue parameter λ with the vortex gyrotropic frequency ω is more complicated if other high-order time derivatives of the vortex core displacement are accounted for in the Thiele equation of motion (Equation (7)). Formally, this can be achieved assuming a velocity-dependent soliton ansatz M(r, t) = M(ρ, X(z, t), dX(z, t)/dt). It was shown in [27] that the mass term is, in general, nonlocal. If the vortex mass is accounted in a local approximation, then the equation of motion of the vortex center position is where M υ = m υ /γ 2 is the vortex mass per unit dot thickness, and m υ is the dimensionless vortex mass.
In the case where Equation (16) is valid, the eigenvalue λ is represented by the function λ(ω) = (ω/ω M ) + 2m υ (ω/ω M ) 2 . If the spectrum λ n of the operator Γ Equation (8a) is found, then the eigenfrequency spectrum can be represented as ω n /ω M = √ 1 + 8m υ λ n − 1 /4m υ . The reduced eigenfrequencies ω n /ω M < λ n for any finite value of the vortex mass m υ , i.e., accounting for the vortex mass, results in a reduction in all vortex gyrotropic frequencies ω n . The vortex mass M υ = m υ /γ 2 can be calculated via the vortex string coupling with the azimuthal spin wave modes having the azimuthal indices m = ±1 [27] according to the concept of the topological gauge field related to the vortex motion [28]. The interaction of the moving vortex core with the azimuthal spin waves was confirmed experimentally [29]. It was recently found that some extra azimuthal spin-wave modes having a curled structure at the dot top and bottom faces appear in the spin excitation spectrum when increasing the dot thickness [30]. Therefore, the calculation of the vortex dynamical mass in thick dots is a complicated problem. The frequency correction due to the vortex mass is small for thin dots [28]. However, the importance of the mass term drastically increases with the dot thickness L increasing in the range of 50-100 nm [27]. The Thiele equation of motion of the soliton center position X and the concept of the soliton mass were recently used to interpret the low-frequency skyrmion excitation modes in CoPt circular dots [31], where a giant value of the mass order of 10 −21 kg was introduced.

Discussion and Conclusions
The magnetic vortex gyrotropic modes in thick cylindrical soft magnetic dots were described for a wide range of the dot thickness according to the concept of the turning points in the magnetostatic potential. This approach allows treating the vortex localized modes (turning points) and nonlocalized modes within a unified picture.
The lowest-frequency vortex gyrotropic mode (ω 0 ) was calculated for the dot aspect ratios β = 0.2 to 2.0 within the approach. The calculated frequency ω 0 (β, r) was lower than the gyrotropic frequency ω G (β) calculated within uniform along thickness mode approximation in [8] at intermediate values of β (<1) and practically coincided with ω G (β) at large values of β (1 < β < 2) (see Figure 1). The good accuracy of the expression for ω G (β) at large β is a surprise because the lowest-frequency eigenmode depicted in Figure 1 is strongly inhomogeneous, localized at the dot faces at such values of β. The frequency ω 0 (β, r) monotonically increases with β increasing. This is in some contradiction with the measurements and simulations conducted in [19,27], where a smooth maximum of the dependence ω 0 (β) was obtained at β ≈ 0.7. We can speculate that the maximum arises due to repulsing of the frequency branches ω n (β) of the lowest (n = 0) and the next (n = 1) vortex gyrotropic modes when the frequencies ω 0 (β) and ω 1 (β) are close to their crossing point [32]. The low-frequency eigenmode localization at the dot faces z = ±1/2 at large values of β (1 < β < 2) is due to strong nonuniformity of the magnetostatic potential I(z, β) near the dot top and bottom faces.
The calculated eigenfrequencies of the nonlocalized vortex string modes (n = 1, 2 . . . ) follow the exchange-dominated decrease ∼ 1/L 2 only for relatively small dot thicknesses and approach some constant values at large values of the dot aspect ratio β (see Figure 3). This effect is beyond the measurements/simulations conducted in [19,21] and should be verified. The eigenmodes s n (z) are odd for n = 1, 3, . . . and even for n = 2, 4, . . . functions of the thickness coordinate z with respect to the dot middle plane z = 0 (see Figure 4). The n = 1 mode profile s 1 (z) deviates essentially from the simple function sin(πz) used in [19], and the deviation increases with the dot aspect ratio β increasing. The profiles of the high-order vortex gyrotropic modes s n (z), n = 2, 3, . . . are close to the standing plane waves supported by the dominating exchange interaction [19].
The vortex gyrotropic modes calculated here for thick circular cylinders are not specific to this dot shape. Similar magnetic vortex excitations should exist also for other dot shapes such as a dome shape or square/rectangular shape. The point is that the dot thickness should be large enough to allow 3D magnetization texture excitations. The inhomogeneous gyrotropic oscillations of the vortex core string can be considered as a step toward understanding the magnetic topological soliton dynamics increasing the system dimensionality from 2D to 3D.