# Exact non-circular symmetric -skyrmions in helical magnets without inversion symmetry

Topological objects of intricate structures have been found in a wide range of systems, like cosmology to liquid crystal, DNA chains, superfluid (Ref. 3), quantum Hall magnets, Bose-Einstein condensates etc. The topological spin texture in helical magnets, observed in recently, makes a new entry into this fascinating phenomenon. The unusual magnetic behaviour of helimagnet MnSi, noticed in recent years, prompted the suspicion that the magnetic states arising in such crystals are of topological nature. Experiments based on the topological Hall effect confirmed such topologically nontrivial states as the skyrmions, located on a plane perpendicular to the applied magnetic field. However, the available models close to MnSi, investigating the formation of skyrmion states, are based mostly on approximate or numerical methods. We present here a theoretical model for chiral magnets with competing exchange and Dzyaloshinskii-Moriya type interaction, which leads to an exact skyrmionic solution with integer topological charge . Such topologically stable spin states with analytic solution on a two-dimensional plane show helical structures of partial order, without inversion and circular symmetry. These exact N-skyrmions, though represent higher excited states, correspond to the lowest energy stable configuration in each topological sector and are likely to appear in MnSi under suitable experimental conditions. The present exact topological solitons, with explicit non-circular symmetry could be applicable also to other fields, where skyrmions are observed, especially in natural systems with less symmetries.

Topological properties can be revealed through a mapping from a continuum space to a differentiable manifold. Therefore for describing topological objects one has to shift from the lattice to a continuum picture. In a magnetic model with the spin configuration varying slowly over the lattice spacing, one can approximate using a long wavelength description, the dimensional lattice to a continuum space . At the same time the associated spin would go to a vector field of unit length , at the classical limit, after a proper renormalisation, which might induce topological invariants of different nature, at different dimensions under suitable boundary conditions .

We intend to use this construction on a two-dimensional -plane perpendicular to the applied magnetic field, to simulate the result on MnSi found through the topological quantum Hall experiment. The physically motivated boundary condition demands that at large distances: , the spin field should go to its vacuum solution, orienting itself to a fixed vector , along the applied field. This condition in turn introduces a nontrivial topology by identifying the infinities of the coordinateÂ space to a single point, which compactifies the vector space to a sphere and defines thus a mapping: , linked to the nontrivial homotopy group (see Fig. 1). Each topological sector is labeled by an integer -valued topological charge , which defined as the degree of this mapping, shows how many times such a skyrmionic field configuration winds up or covers the target sphere, when the coordinate space is swept once. The corresponding charge density can be expressed through the spin vector field as

(1) |

Figure 1 The induction of nontrivial topology. The two-dimensional coordinate space , due to the boundary condition compactifies to a 2-sphere , with its north-pole identified to all points at the space-infinities. The unit vector field , therefore describes a sphere to sphere mapping with the integer-valued topological charge , defined as the degree or the winding number of this mapping.

Another important idea, that we use for constructing our magnetic model is to maintain a lower bound for energy of the system through its topological charge . Note that, such an energy bound would guarantee the stability of the finite energy solution in each topological sector, since the topological charge is conserved independent of the dynamics of the system and can not be changed by unwinding the spin configuration by smooth transformations. Moreover, when the energy reaches its lowest value saturating this bound: , called Bogomolny limit, some extraordinary thing happens in almost all known models : it yields an exact solution to the field model, in some cases for all topological sectors.

The unusual properties of the helical magnets like MnSi and especially their partial helical order created by the skyrmion spin states are believed to be due to the competing forces between the ferromagnetic exchange and an effective Dzyaloshinskii-Moriya (DM) type interaction. The ferromagnetic interaction tries to align the spins parallel to each other, while the DM-term with the broken inversion symmetry tends to orient them perpendicular to each other, settling finally to a partial helical order of topological origin. Moreover, the skyrmionic spin texture in MnSi, as observed in the recent topological Hall experiments, appears in the plane perpendicular to the applied magnetic field. Therefore to model such a system we construct our Hamiltonian from a standard Heisenberg ferromagnet in two-dimensions: where the spin , located at site interact with its nearest-neighbour at , in the 2-dimensional lattice. We include next the DM term: which in spite of its traditional form, bears some crucial differences in the structure of its Dzyaloshinskii (D)-vector , which would result to our field model allowing exact skyrmion solution. Firstly, we take the D-vector aligned along the applied field, which leaves nontrivial only its third-component , breaking the rotational symmetry in the internal spin-space. Our next deciding assumption is the nonlinear nature of this vector, given through a spin-dependent structure: . Note that the second factor obeys the usual antisymmetric exchange of the D-vector, while the first factor, which we take in the explicit form manifestly breaks the space-inversion symmetry of the vector: inducing the same property to the Hamiltonian of the system as usual through the DM interaction.

Our total Hamiltonian, which has the required broken inversion-symmetry and the competing spin interactions of opposite trends, should now be tested for its topological properties, for which we have to go for the continuum limit . For constructing the corresponding field model from the 2-dimensional lattice model, we have to take carefully the vanishing limit of the lattice constants in both the lattice directions , which would result, for example, etc. This procedure reduces our lattice spin-model to a field Hamiltonian with

(2) |

It is not difficult to see that, at this continuum limit the chirality of three non-coplanar spins gets linked to the topological charge density , expressed through the -field as (1).

It is important to note, that by tuning the coupling constants , our field Hamiltonian , as can be shown by using a simple school-level geometric inequality, would acquire the crucial lower bound : , through the topological charge , the significance of which for the topological stability of the solutions is emphasized above. Focusing now on the conjecture, that the Bogomolny limit, at which the energy bound is saturated, should yield the exact finite-energy soliton solution, we look for the situation when this limit could be reached for nontrivial contributions from both parts of the Hamiltonian. We find that this is possible for a field configuration , where the space-inversion as well as the circular symmetry is lost. More precisely, the spin texture should be such that, the directions and , defining the helical order, would cross each other at an angle , with three field vectors: , and forming an isosceles triangle. Under this special geometric condition, the energy bound is exactly saturated attaining its lowest value in each sector with . Analysing this intriguing geometry we can extract the exact magnetic field solution as

(3) |

in polar coordinates with a constant scaling parameter . This exact solution, which can be checked by direct insertion, obeys the required boundary condition and yields the topological charge . The distribution of the charge density , which coincides also with the energy density at the solution manifold is shown in Fig. 2, together with the magnetization and the spin component , calculated for the exact result (3), specific to the first excited state with . The figure shows that at the skyrmion axis the charge density becomes maximum, vanishing gradually with the distance, while the magnetisation and the spin component exhibit a toroidal structure without circular and space-inversion symmetry. Their maximum/minimum is attained at a distance , varying with , which describes an intriguing form of the torus, with its radius spreading out with increasing polar angle. As revealed in quantum Hall experiments, the peculiarly large Hall conductance is created by an effective magnetic field proportional to the topological charge density. Therefore the detailed and exact form of the charge density together with the magnetization and the spin field component, found here based on the exact skyrmion solution, should serve as a guiding resource for the future precision experiment.

a

b c

Figure 2 Exact skyrmion state with . a, The charge density which coincides also with the energy density for the exact skyrmion. b, The magnetization and c, Spin component , showing toroidal structure with explicitly broken inversion and circular symmetry.

The graphical result based on the numerical calculation on a model for proposed in Ref 13, shows that the skyrmion lattice derived as a spontaneous ground state solution, preserves the circular as well as the space-inversion symmetry, when projected on two-dimensions. However, our exact result (3) based on our magnetic model with the broken space-inversion and the spin rotational symmetry, induced by the DM interaction, gives a different picture. It shows, as depicted in Fig. 3, that the partial helical order of the spin texture on the -plane manifestly breaks both the circular and the inversion symmetry. The future experiments capable of identifying the finer details of the topological spin texture should be able to settle this issue. At the theoretical level, it is significant to note also that, the exact topological soliton without circular symmetry, presented here, is a nontrivial generalisation over the well known solution of the nonlinear -model, exhibiting this symmetry. In fact all other known exact solutions of the field models, obtained at the Bogomolny limit, are basically spherical symmetric, except only the present one.

Another valuable information, e.g. the detail spin structure in higher topological sectors with arbitrary integer charge , can be found easily from our exact solution (3). Although the solution with higher topological charge was investigated by using a combination of analytic and numerical methods in connection with the Hall skyrmions and recently for the bilayer graphine, no such study, as far as we know, is undertaken for the magnetic models. Therefore the present result, designed specifically for the magnets like MnSi. and more significantly having purely analytic result with manifestly broken inversion symmetry, is novel in many respect. Such an exact skyrmion with higher topological charge , is shown in Fig. 3b. These stable states in higher sectors with more intricate texture, as we see in the figure, would be challenging to detect in topological Hall experiments.

The topological soliton with non-circular symmetry, presented here, is exceptional as an exact field theoretic solution and would be interesting for its detection and applications in other systems, lacking such a symmetry.

a b

Figure 3 Skyrmion spin texture with partial helical order,
described by the exact solution of the spin vector field a , with
the topological charge and
b, . Both the figures demonstrate manifestly broken space-inversion and the
circular symmetry on the -plane, with a distinct difference in their
helical pattern, describing different degrees of mapping in these
two cases.

## References

- [1] Bouick, M. J., Chandar, L., Schiff, E. A. & Srivastava, A. M. Cosmological Kibble mechanism in the laboratory: string formation in liquid crystal. Science 263 943-945 (1994)
- [2] Travers, A. DNA topology: dynamic DNA looping.Current Biology 16 R838-R840 (2006)
- [3] Bäuerie, C. et al. Laboratory simulation of cosmic string formation in the early universe using superfluid He. Nature 382 332-334 (1996)
- [4] Sondhi, S. L., Karlhede, A. & Kivelson, S. A. Skyrmion and the crossover from integer to fractional quantum Hall effect at small Zeeman energies. Phys. Rev. B 47, 16419-19426 (1993).
- [5] Lillihöök, D., Lejnell, K., Karlhede, A. & Sondhi, S. L. Quantum Hall skyrmions with higher topological charge. cond-mat.9704121
- [6] Al’Khawaja, U. & Stoof H. T. C. Skyrmions in a ferromagnetic Bose-Einstein condensate. Nature 411, 918-920 (2001)
- [7] Lee, M., Kang, W., Onose, Y., Tokura, Y.& Ong, N. P. Unusual Hall effect anomaly in MnSi under pressure. Phys. Rev. Lett. 102, 186601 1-4 (2009).
- [8] Neubauer, A. et. al. Topological Hall effect in the A phase of MnSi. Phys. Rev. Lett. 102, 186602 1-4 (2009).
- [9] Pfleiderer, C., Julian, S. R. & Lonzarich, G. G. Non-Fermi-liquid nature of the normal state of itinerant-electron ferromagnets. Nature 414, 427 - 430 (2001).
- [10] Pfleiderer, C. Non-Fermi-liquid puzzle of MnSi. Physica B 328, 100 - 104 (2003).
- [11] Dolron-Leyraud, N. et. al. Fermi-liquid breakdown in the parametric phase of a pure metal. Nature 425, 595 - 599 (2003).
- [12] Pfleiderer, C. et. al. Partial order in the non-Fermi-liquid phase of MnSi. Nature 427, 227 - 231 (2004).
- [13] Röler, U. K., Bogdanov, A. N. & Pfleiderer, C. Spontaneous skyrmion ground states in magnetic metals. Nature 442, 797 - 801 (2006).
- [14] Binz, B. & Vishwanath, A. Chirality induced anomalous - Hall effect in helical spin crystals. Physica B 403, 1336 - 1340 (2008).
- [15] Mühlbauer, S. et. al. Skyrmion lattice in a chiral magnet. Science 323, 915 - 919 (2009).
- [16] Skyrme, T. H. A nonlinear field theory, Proc. R. Soc. Lond. Ser. A 260, 127 - 138 (1961).
- [17] Belavin, A. A. & Polyakov, A. M. Metastable states of two-dimensional Heisenberg ferromagnet. JETP Lett. 22, 245 - 246 (1975).
- [18] Honerkamp, J., Patani, A., Schlindwein, M. & Shafi, Q. On kink solution in two space dimensions. Nuovo Cimento Lett. 15, 97-100 (1976).
- [19] Kundu, A. Instanton solutions in the anisotropic - model. Phys. Lett. B 110, 61-63 (1982).
- [20] Fogedby, H.C. Classical and quantum aspects of the continuous Heisenberg chain at T=0. Z. Phys. B 41, 115-127 (1981)
- [21] Duff, M. J. & Isham, C. J. Soliton and vortex type solutions in nonlinear chiral theories. Nucl. Phys. B 108, 130 - 140 (1976).
- [22] Faddeev, L. & Niemi, A. J. Stable knot-like structures in classical field theory Nature 387, 58-61 (1997)
- [23] Felsager, B & Leinaas, J. M. A conformally invariant -model in 4-dimensions. Phys. Lett. B 94, 192-194 (1980)
- [24] Bogomolny, E. B. The stability of classical solutions. Sov. J. Nucl. Phys. 24, 449-454 (1976)
- [25] Actor, A. Rev. Mod. Phys. 51, 451-525 (1979) Classical solutions of Yang-Mills theories.
- [26] Abanin, D. A., Parameswaran, S. A. & Sondhi, S. L. Charge 2e skyrmions in Bilayer graphine. arXiv:0904.0040 [cond-mat.mes-hall] (1 Apr. 2009)