# Confinement of Skyrmion states in noncentrosymmetric magnets

###### Abstract

Skyrmionic states in noncentrosymmetric magnets with Lifshitz invariants are investigated within the phenomenological Dzyaloshinskii model. The interaction between the chiral Skyrmions, being repulsive in a broad temperature range, changes into attraction at high temperatures. This leads to a remarkable confinement effect: near the ordering temperature Skyrmions exist only as bound states, and Skyrmion lattices are formed by an unusual instability-type nucleation transition. Numerical investigations on two-dimensional models demonstrate the confinement and the occurrence of different Skyrmion lattice precursor states near the ordering transition that can become thermodynamically stable by anisotropy or longitudinal softness in cubic helimagnets. We introduce a new fundamental parameter, confinement temperature separating the peculiar region with ”confined” chiral modulations from the main part of the phase diagram with regular helical and Skyrmion states.

###### pacs:

75.10.-b 75.30.Kz 61.30.Mp 03.75.LmIn certain noncentrosymmetric magnetic systems, the asymmetric Dzyaloshinskii-Moriya (DM) exchange results in chiral couplings that can stabilize long-period, non-collinear modulations of the magnetization with a fixed sense of rotation Dz64 (); PRB02 (). These chiral couplings are phenomenologically described by Lifshitz invariants that destroy the homogeneity of ordered phases Dz64 (). In highly symmetric systems with Lifshitz invariants multiple modulations occur as textures with localized twists in two or more spatial directions Hornreich82 (); Wright89 (). In noncentrosymmetric magnetic systems with Lifshitz invariants these double-twist motifs exist as smooth (baby)-Skyrmions JETP89 (); JMMM94 (); Nature06 (), static solitonic textures localized in two spatial directions, which can be extended into the third direction as Skyrmion strings or Hopfions. These magnetic Skyrmions are stabilized solely by the chiral DM couplings Nature06 (); JMMM94 (), which prevent a spontaneous collapse into topological singularities. Skyrmionic matter is created by the condensation of these solitonic units, similar to vortex matter in type-II superconductors JETP89 (). Just such chiral Skyrmions have been recently observed in thin layers of noncentrosymmetric ferromagnet (Fe,Co)Si Yu10 (). Skyrmionic states stabilized by Lifshitz-type invariants may exist and form extended mesophases in various condensed matter systems, as chiral liquid crystals, ferroelectrics, multiferroics, and in confined achiral systems (e.g., thin magnetic layers) PRL01 (); Bode07 (); Wright89 (); Tkalec09 ().

Skyrmionic textures usually form through nucleation, following a classification introduced by DeGennes DeGennes75 () for (continuous) transitions into incommensurate modulated phases. As demonstrated in Refs. JETP89 (); JMMM94 (), for the low-temperature micromagnetic model of chiral magnets, at the transition from the field-driven Skyrmion lattices into the polarized homogeneous state the lattice period diverges and the Skyrmions are set free as localized excitations. As the Skyrmions retain their size and axisymmetric shape, there is a full spectrum of lattice modes up to the transition, in contrary to an instability type transition where the amplitude of one fundamental mode acts as small parameter of the transition DeGennes75 ().

Here, we show for the standard model of chiral isotropic ferromagnets Dz64 (); Bak80 (); Nature06 () that Skyrmions are confined very close to the ordering temperature. In that part of the phase diagram, the creation of Skyrmions as stable units and their condensation into extended textures occurs simultaneously through a rare case of an instability-type nucleation transition Felix86 (), but the confined Skyrmions as discernible units may arrange in different mesophases. This is a consequence of the coupling between the magnitude and the angular part of the order parameter. Thus, near the ordering transitions, the local magnetization is not only multiply twisted but also longitudinally modulated. From numerical investigations on 2D models of isotropic chiral ferromagnets, a staggered half-Skyrmion square lattice at zero and low fields and a hexagonal Skyrmion lattice at larger fields are found in overlapping regions of the phase diagram near the transition temperature. Furthermore, the thermodynamic stability of Skyrmionic states can be favoured with respect to one-dimensional modulations by supplementing the model with cubic exchange anisotropy or in a modified model for metallic chiral magnets Nature06 ().

Within the standard phenomenological (Dzyaloshinskii) theory Dz64 () the magnetic energy density of a noncentrosymmetric ferromagnet can be written in the dimensionless form Bak80 (); Nature06 ()

(1) |

Here, reduced values of the spatial variable , the magnetization , and the applied magnetic field (, is a unity vector along ), are expressed via the parameters of the energy density Bak80 () (where , and ): , , , , . DM energy consists of Lifshitz invariants Dz64 ()

(2) |

For noncentrosymmetric uniaxial ferromagnets DM functionals are listed in JETP89 (). In particular, for isotropic and cubic helimagnets Bak80 () (for details see appendix). Functional (1) includes only the basic (isotropic) interactions essential to stabilize chiral modulations. This depends on three internal variables (components of the magnetization vector ) and two control parameters, the reduced magnetic field amplutude and the ”effective” temperature . In chiral magnetism models (1) play a fundamental role and are similar to the Frank energy in liquid crystals and Landau-Ginzburg functional in superconductivity. We consider 2D chiral modulations homogeneous along the applied field and modulated in the plane perpendicular .

Isolated and embedded Skyrmions. The equations minimizing functional (1) include solutions for axisymmetric localized states (isolated Skyrmions), , and , where we use spherical coordinates for the magnetization and cylindrical coordinates for the spatial variables . The solutions for all noncentrosymmetric classes have been derived in JETP89 (). For cubic helimagnets and uniaxial systems with 22 symmetry, (Fig. 1 a). Profiles and of isolated Skyrmions are derived from numerical solutions of the Euler equations, as given in Ref. Nature06, and common for all noncentrosymmetric classes with Lifshitz invariants, with the boundary conditions , , , ( is the magnetization in the saturated state). For extended textures of two-dimensional models, the functional has been investigated by numerical energy minimization using finite-difference discretization on rectangular grids with adjustable grid spacings Nature06 ().

Isolated Skyrmions (Fig. 1 (a)) exist only below a critical line and condense into a hexagonal lattice (Fig. 2) below a field , which marks the transition between the Skyrmion lattice (SL) and the homogeneous paramagnetic state (see Fig. 2 (a), (b)). Near the ordering temperature a square half-Skyrmion lattice (Fig. 1 (c)) has lower energy than the hexagonal lattice. Half-SLs consist of cells with up and down magnetization in the center and in-plane magnetization along the cell boundaries. Such cells have a topological charge . In the hexagonal SLs the magnetization at the boundaries (center) is (anti)parallel to the applied field. The cells bear unit topological charges.

Confinement. By solving the linearized Euler equations for the asymptotics of isolated Skyrmions () one finds three distinct regions in the magnetic phase diagram with different character of Skyrmion-Skyrmion interactions (Fig. 2 (a)): repulsive interactions in a broad temperature range (area (I)) are changed to attractive interaction at higher temperatures (area (II)). Finally in area (III) near the ordering temperature strictly confined Skyrmions exist. These regions are separated by the line

(3) |

with turning points (), (), and (-0.5, 0) (dashed line in Fig. 2 (a)). In the major part of the phase diagram, the evolution of SLs under a magnetic field closely agrees with the behavior studied earlier for the low-temperature limit JETP89 (); JMMM94 () and the transition mechanism at the high-field limit is of the nucleation type with isolated Skyrmion excitations appearing below the instability line . The temperature of the turning point (), a confinement temperature determines the frontier between low- () and high-temperature () chiral modulations. Usually and high-temperature modulations are restricted by a narrow vicinity of the Neel temperature. According to calculations in Nature06 () for MnSi high-temperature Skyrmions exist in temperature interval = 2 K below the Curie temperature .

Evolution of confined Skyrmions. In the vicinity of the ordering temperature the Skyrmion lattices exist within the confinement pocket (III). In this region Skyrmion states drastically differ from those in the main part of the phase diagram. Due to the ”softness” of the magnetization modulus the field-driven transformation of the Skyrmion lattices evolves by distortions of the modulus profiles both in the hexagonal and square Skyrmion lattices while the equilibrium periods of the lattices do not change strongly with increasing applied field (Fig. 3). Despite the strong transformation of their internal structures the Skyrmion lattices preserve axisymmetric distribution of the magnetization near the centers of the Skyrmion lattice cells (Fig. 3). This remarkable property reflects the basic physical mechanism underlying the formation of Skyrmion lattices. The local energetic advantage of Skyrmion lattices over helicoids is due to a larger energy reduction in the “double-twisted” Skyrmion cell core compared to “single-twisted” helical states Nature06 (); Wright89 (). This explains the unusual axial symmetry of the cell cores and their stability. An increasing magnetic field gradually suppresses the antiparallel magnetization in the cell core reducing the energetic advantage of the “double-twist” and increases the overall energy of the condensed Skyrmion lattice. At critical field a first-order transition occurs into the polarized paramagnetic state. In the interval the hexagonal Skyrmion lattice exists as a metastable state. At the lability field the magnetization modulus in the cell center becomes zero (see magnetization profile for in Fig. 3 (c)). At this field, the ”double-twist” region is suppressed, and the lattice loses its stability.

At zero field with increasing temperature the magnetization modulus in hexagonal and square lattices gradually decreases to zero at the ordering temperature. This is the instability-type nucleation transition into the paramagnetic phase : the order parameter becomes zero at the transition point (as in the instability mode), however, the lattices retain their symmetry and the arrangement of axisymmetric Skyrmions up to the critical point. In finite fields in the region (II) the attractive Skyrmion-Skyrmion interaction means that multi-Skyrmions can always form bound states. Near such clusters of Skyrmions are more stable than the isolated Skyrmions. This is seen in Fig. 2 (a),(b). For temperatures above point , , the metastable SL as the densest infinite cluster is more stable than isolated Skyrmions. Thus, in region (II) of the phase diagram, there is clustering and, for higher temperatures, confinement of isolated Skyrmion excitations.

The predicted Skyrmionic textures and the pre-cursor phenomena associated with the confinement of Skyrmions near the ordering transition are observable in magnetic materials with appropriate symmetry. In a large group of uniaxial noncentrosymmetric magnets the DM energy is described by gradients only along directions perpendicular to the axis (e.g. multiferroic BiFeO, space group and antiferromagnetic BaCuGeO, space group PRB02 (); Sosnowska82 ()). In such magnets one-dimensional modulations with the propagation vector in the basal plane (helicoids) exist in broad ranges of the magnetic fields Dz64 (); Mukamel85 (). In the confinement pocket of the phase diagram (III) the helicoids also exist only as bound states Mukamel85 (). They have the lowest energy at lower fields, and the hexagonal SL corresponds to the global energy minimum at higher fields (Fig. 2 (c)).

Cubic helimagnets. In other groups of chiral magnets the DM energy includes contributions with gradients along all three spatial directions. This stabilizes chiral modulations with propagation along the direction of an applied field as cone phases Bak80 (). For the isotropic model (1) the cone phase solution with the fixed magnetization modulus and rotation of around the applied magnetic field: , is the global energy minimum in the whole region where modulated states exist. For cubic helimagnets, the energy density (1) has to be supplemented by anisotropic contributions, , where and are reduced values of exchange and cubic anisotropies Bak80 (). These anisotropic interactions impair the ideal harmonic twisting of the cone phase and lead to the thermodynamic stability of Skyrmion states, as shown in the equilibrium phase diagram Fig. 4 (a). The difference between the energy of the hexagonal Skyrmion lattice and of the cone phase calculated for the isotropic model, , has minima along a curve which reaches the critical point as ). Weak exchange anisotropy of a cubic helimagnet, therefore, creates a pocket around , where the hexagonal Skyrmion lattice becomes the global energy minimum in a field (Fig. 4 (a)). This case is realized in cubic helimagnets with negative exchange anisotropy () as in MnSi Bak80 (). This anisotropy effect provides a basic mechanism, by which a Skyrmionic texture is stabilized in applied fields, as observed experimentally as so-called ”A-phase” in MnSi Gregory92 (); Lebech95 (); Muhlbauer09 (). The exchange anisotropy also leads to the thermodynamic stability of half-SLs (Fig. 4 (a)). The stabilization of such textures may be responsible for anomalous precursor effects in cubic helimagnets in zero field Pfleiderer04 (); Nature06 (); Pappas09 ().

Chiral modulations in non-Heisenberg models. A generalization of isotropic chiral magnets proposed in Nature06 () replaces the usual Heisenberg-like exchange model by a non-linear sigma-model coupled to a modulus field with different stiffnesses. This yields a generalized gradient energy for a chiral isotropic system with a vector order parameter, which is equivalent to the phenomenological theory in the director formalism Wright89 (); Nature06 () = . Parameter equals unity for a “Heisenberg” model, in chiral nematics Wright89 (). Confined chiral modulations are very sensitive to values of . The magnetic phase diagram calculated for includes pockets with square half-SLs, hexagonal SLs with the magnetization in the center of the cells parallel to the applied field, and helicoids with propagation transverse to the applied field (Fig. 4 b).

The confinement effects on chiral Skyrmions strongly changes the picture of the formation and evolution of chiral modulated textures and shed new light on the problem of precursor states observed as blue phases in chiral nematics Wright89 () and in chiral magnets Gregory92 (); Lebech95 (); Pfleiderer04 (); Nature06 (); Pappas09 (); Muhlbauer09 (). The results show that confinement / deconfinement transitions of localized string-like solitons can be realized in these condensed matter systems. They provide counterparts of formation mechanisms for extended microscopic matter from topological solitons, as devised originally in the Skyrme model Skyrme61 ().

Acknowledgment. We thank K. v. Bergmann, G. Bihlmayer, S. Blügel, H. Eschrig, S.V. Grigoriev, J.-H. Han, S. Heinze, R. Mössner, C. Pappas, W. Selke, H. Wilhelm, and M. Zhitomirsky for discussions. Support by DFG project RO 2238/9-1 is gratefully acknowledged.

## Appendix A APPENDIX:Dzyaloshinskii theory of chiral helimagnetism

The appendix includes technical comments on the mathematical structure of the interaction functional (1), on the mathematical relation between this model and those for different classes of noncentrosymmetric magnetic systems and chiral liquid crystals. The equations for isolated and bound Skyrmions are provided along with a short information on the numerical methods to solve these problems.

### a.1 1. Basic model.

In noncentrosymmetric ferromagnets chiral asymmetry of exchange coupling related to the quantum-mechanical Dzyaloshinskii-Moriya interactions Dz57 () induces long-period modulations of the magnetization vector with a fixed rotation sense Dz64 (). The equilibrium chiral modulations are derived by minimization of the magnetic energy which can be written in the following general form Dz64 ()

(4) |

where is the exchange stiffness constant, is an applied magnetic field, collects short-range magnetic interactions independent on spatial derivatives of the magnetization, and Dzyaloshinsky-Moriya energy is composed of antisymmetric terms linear with respect to first spatial derivatives of the magnetization vector

Functional forms of energy contributions are determined by crystallographic symmetry of a noncentrosymmetric magnetic crystal Dz64 (); Bak80 (); JETP89 (). Particularly, for important cases of noncentrosymmetric ferromagnets belonging to cubic (23, 432) and uniaxial crystallographic classes (mm, 2m) the DM energy contributions have the following form Dz64 (); Bak80 (); JETP89 ()

(5) |

where and is a Dzyaloshinskii constant. For other noncentrosymmetric classes this energy includes two or more Dzyaloshinskii constants JETP89 ().

Near the ordering temperatures the magnetization amplitude varies under the influence of the applied field and temperature. Commonly this process is described by including into the magnetic energy an additional term Bak80 ()

(6) |

This model is fundamental to the systematic phenomenological description of magnetic states in noncentrosymmetric magnetic systems Dz64 (). Interaction functional (4) plays in chiral magnetism a similar role as the Frank energy in liquid crystals DeGennes () or Ginzburg-Landau functional in physics of superconductivity Brandt95 (). Its form is dictated by the natural extension of Landau’s approach to ordering transtions for modulated systems which break the third Lifshitz condition, as pioneered by Dzyaloshinskii in particular in magnetism Toledano87 (). As a result, this type of theory constitutes the canonical form of statistical field theories for condensed matter systems, where particular couplings and symmetry enable Lifshitz invariants to induce modulated phases. Extensions including higher-order secondary effects in magnetism (like anisotropies or dipolar interactions) are widely used to describe evolution of modulated states in many chiral magnets Izyumov84 (); PRB02 (); JMMM94 (); Nature06 (), including metamagnetic transitions in cubic helimagnets induced by high pressure Bak80 (); MnSi (); Plumer81 ().

By rescaling the spatial variable in (1),(6) , the magnetic field , the magnetization

(7) |

energy (6) can be written in the following reduced form (Eq. (1))

(8) |

where and is a unity vector along the applied magnetic field. Functional (8) includes three internal variables (components of the magnetization vector ) and two control parameters, the reduced magnetic field amplutude and the ”effective” temperature (7). By direct minimization of Eq. (8) we have derived one-dimensional (kinks, helicoids, conical helices or spirals) and two-dimensional solutions (isolated and bound Skyrmions) for arbitrary values of the control parameters. These results are collected in the phase diagram of solutions (Fig. 2).

Energy (8) includes only basic interactions essensial to stabilize Skyrmionic and helical phases. Solutions for chiral modulated phases and their most general features attributed to all chiral ferromagnets are determined by interactions functional (8). Generically, there are only small energy differences between these different modulated states. On the other hand, weaker energy contributions (as magnetic anisotropy, stray-fields, magneto-elastic couplings) impose distortions on solutions of model (8) which reflect crystallographic symmetry and values of magnetic interactions in individual chiral magnets. These weaker interactions determine the stability limits of the different modulated states (Fig. 4). The fact that thermodynamical stability of individual phases and conditions of phase transfomations between them are determined by magnetocrystalline anisotropy and other relativistic or weaker interactions means that (i) the basic theory only determines a set of different and unusual modulated phases, while (ii) the transitions between these modulated states, and their evolution in magnetization processes depends on symmetry and details of magnetic secondary effects in chiral magnets, in particular the strengths of relativistic magnetic interactions. Thus functional (8) is the generic model for a manifold of interaction functionals describing different groups of noncentrosymmetric magnetic crystals, because it allows to identify the basic modulated structures that may be found in them.

### a.2 2. Equations for isolated and bound Skyrmions

The Euler equations for functional (8) have solutions for axisymmetric isolated structures of type , , (we introduce here spherical coordinates for the magnetization and cylindrical coordinates for the spatial variables ) Nature06 (). Solutions are known for all uniaxial and cubic noncentrosymmetric ferromagnets JETP89 (). Particularly, for cubic helimagnets, and for uniaxial ferromagnets with and symmetry, correspondingly. After substitution of solutions into (8) and integration with respect to the total energy of a Skyrmion (per unit length along ) is with energy density

(9) |

where a common convention is applied. The Euler equations for the functional (9)

(10) |

with boundary conditions , , describe the structure of isolated Skyrmions (the magnetization of the homogeneous phase is derived from equation ). Typical solutions , of Eqs. (10) are plotted in Fig. 1.

For present work, we have evaluated numerically the solutions for the energy functionals of type Eq. (1) within a standard approach using finite differences for gradient terms and adjustable grids to accommmodate modulated states with periodic boundary conditions. Search for modulated states and energy minimization was done by Monte Carlo simulated annealing. For 2-dimensional models, this approach can give converged results, as checked by comparison with analytical results. Including secondary effects, the approach is able and has been used by us to calculate thermodynamic stability of phases, and to study the complete temperature-field phase diagrams in terms of equilbrium (mean-field) states of the energy functionals Eq. (1). As we are investigating here long-period modulated phases, the phase diagrams are expected to give qualitatively the correct results for these phenomenological models.

## References

- (1) I.E. Dzyaloshinskii, Sov. Phys. JETP 19, 960.
- (2) A. N. Bogdanov et al., Phys. Rev. B 66, 214410 (2002);
- (3) R.M. Hornreich et al., Phys. Rev. Lett. 48, 1404 (1982).
- (4) D. C Wright, N. D. Mermin, Rev. Mod. Phys. 61, 385 (1989).
- (5) A.N. Bogdanov and D.A. Yablonsky, Sov. Phys. JETP 68, 101 (1989).
- (6) A. Bogdanov, A. Hubert, J. Magn. Magn. Mater. 138, 255 (1994); ibid. 195, 182 (1999); A. Bogdanov, JETP Lett. 62, 247 (1995); A.N. Bogdanov et al., Physica B 359, 1162 (2005).
- (7) U.K. Rößler et al., Nature 442, 797 (2006).
- (8) X. Z. Yu et al., Nature 465, 901 (2010).
- (9) A. N. Bogdanov, U. K. Rößler, Phys. Rev. Lett. 87, 037203 (2001).
- (10) M. Bode et al., Nature 447 190 (2007).
- (11) U. Tkalec et al., Phys. Rev. Lett. 103, 127801 (2009).
- (12) P.G. de Gennes, in Fluctuations, Instabilities, and Phase transitions, ed. T. Riste, NATO ASI Ser. B, vol. 2 (Plenum, New York, 1975).
- (13) P. Bak and M. H. Jensen, J. Phys.C: Solid State Phys. 13, L881 (1980); O. Nakanishi et al., Solid State Comm. 35, 995 (1980).
- (14) J. W. Felix et al., Phys. Rev. Lett. 57, 2180 (1986).
- (15) I. Sosnowska et al., J. Phys. C 15, 4835 (1982); A. Zheludev et al. Phys. Rev. Lett. 78, 4857 (1997).
- (16) B. Schaub, D. Mukamel, Phys. Rev. B 32, 6385 (1985).
- (17) C. I. Gregory et al., J. Magn. Magn. Mater. 104-107, 689 (1992); D. Lamago et al., Physica B 385-386, 385 (2006).
- (18) B. Lebech et al., J. Magn. Magn. Mater. 140, 119 (1995).
- (19) S. Muhlbauer et al., Science, 323, 915 (2009).
- (20) C. Pfleiderer et al., Nature 427, 227 (2004).
- (21) C. Pappas et al., Phys. Rev. Lett. 102, 197202 (2009).
- (22) T. H. R. Skyrme, Proc. Roy. Soc. Lon. 260, 127 (1961).
- (23) I. E. Dzyaloshinskii, Sov. Phys. JETP 5, 1259 (1957), T. Moriya, Phys. Rev. 120, 91 (1960).
- (24) P. G. De Gennes and J. Prost, The Physics of Liquid Crystals (Oxford University Press, Oxford, 1993), 2nd ed., D. C Wright, N. D. Mermin, Rev. Mod. Phys. 61, 385 (1989).
- (25) E. H. Brandt, Rep. Prog. Phys. 58, 1465 (1995), Phys. Rev. B, 68, 054506 (2003).
- (26) See, e.g., J.-C. Tolédano, P. Tolédano, The Landau Theory of Phase Transitions, World Scientific, Singapore 1987, chap. V and chap. VI,6.
- (27) Yu. A. Izyumov, Sov. Phys. Usp. 27, 845 (1984).
- (28) G. G. Lonzarich, L. Taillefer, J. Phys. C. 18, 4339 (1985), M. Yamada et al. J. Alloys Comp. 364, 37, (2004).
- (29) M. L. Plumer, M. B. Walker, J. Phys. C. 14, 4689 (1981), M. L. Plumer, J. Phys.:Cond Matt. 2, 7503 (1990).