Low-energy theory of a two-dimensional itinerant chiral magnet

Low-energy theory of a two-dimensional itinerant chiral magnet

Armin Rahmani, Rodrigo A. Muniz, and Ivar Martin Theoretical Division, T4 and CNLS, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
International Institute of Physics - UFRN, Natal, RN 59078-400, Brazil
Department of Physics and Institute for Optical Sciences, University of Toronto, Toronto, ON, M5S 1A7, Canada
Materials Science Division, Argonne National Laboratory, Argonne, IL 60439, USA
Abstract

Effective actions, such as nonlinear sigma models, are important tools in studying low-energy fluctuations of magnetically ordered states. Here we derive an effective action for the smooth order-parameter distortions of noncoplanar magnetic textures [with an SO(3) order parameter, as opposed to the SOSO(2) of collinear magnets] coupled to itinerant electrons. Noncoplanar magnetic textures commonly arise in the presence of electron-mediated long-range spin exchange interactions, giving rise to the emergence of rich electronic phenomena such as anomalous quantum Hall effect. We parameterize the smooth order-parameter distortions of such a magnetic texture by certain non-Abelian fields, and derive the action in terms of these fields by integrating out the itinerant electrons. As a concrete example, we perform explicit calculations for a triangular-lattice model with tetrahedral magnetic ordering. The action encodes the charge and spin quantum-Hall responses as well as the energetics of twisting the magnetic texture.

I introduction

Beyond the familiar collinear magnetically ordered states (such as Néel antiferromagnets), a myriad of complex magnetic structures can emerge in frustrated magnetic systems even in the classical large- limit: The next level of complexity arises when the order parameter is coplanar but not collinear (such as the -degree ordering of a Heisenberg antiferromagent on the triangular lattice). More exotic noncoplanar magnetic textures are in some cases realized by classical magnetic moments with short-range interactions. Momoi et al. (1997); Messio et al. (2011); Lapa and Henley (2013) Electron-mediated interactions in itinerant systems (which have a long-range character), on the other hand, commonly give rise to such exotic orders. Indeed, many models of large- local moments, residing on geometrically frustrated lattices and coupled to itinerant electrons, exhibit energetically stable phases characterized by noncoplanar magnetic textures. Martin and Batista (2008); Akagi and Motome (2010); Kumar and van den Brink (2010); Kato et al. (2010); Chern (2010); Li (2012); Yu and Li (2012); Venderbos et al. (2012); Chern et al. (2011) Due to the presence of a nonvanishing scalar spin chirality, such magnetically ordered itinerant systems may exhibit rich electronic phenomena including spontaneous quantized integer quantum Hall effect Ohgushi et al. (2000); Shindou and Nagaosa (2001); Martin and Batista (2008) and fractionalization on topologically stable defects. Rahmani et al. (2013)

In addition to specifying the stable magnetic structures, a complete characterization of these exotic phases requires the identification of low-energy magnetic excitations. Such understanding can be obtained from effective actions for long-wave-length fluctuations around magnetically ordered states [the classic examples being the nonlinear sigma models with (without) topological terms for Néel antiferromagnet in one (two) dimensions Haldane (1983); Affleck (1985); Haldane (1988); Ioffe and Larkin (1988); Fradkin and Stone (1988); Dombre and Read (1988); Wen and Zee (1988)]. Such effective action is not currently available for noncoplanar spin textures stabilized by electron-mediated interactions. In this paper, we take a step toward characterizing the low-energy magnetic excitations, by deriving an effective action for smooth distortions of the order-parameter. (In addition to such smooth distortions, noncoplanar magnetic textures may also be distorted by fast fluctuating modes, which are not addressed in the present paper.) The action (i) determines the energetics of smooth distortions, i.e., it allows us to compute the excess energy of any smoothly distorted texture and hence can be used to evaluate, e.g. the interaction potential between two vortices, and, (ii) encodes the electronic responses to time- and position-dependent perturbations, including the Hall response.

Unlike Néel states, where the order parameter can be represented by a unit vector , here, the order parameter is specified by a full three-dimensional rotation matrix, which can be parametrized by a unit vector and a scalar (respectively representing an axis and an angle of rotation of a reference noncoplanar configuration). Kawamura and Miyashita (1984); Dombre and Read (1989) We show in this paper that it is convenient to parameterize the smooth distortions of the order parameter in terms of non-Abelian SU-gauge-like fields with , (we limit ourselves here to two-dimensional systems) where , , are the Pauli matrices and . This characterization of spin fluctuations not only encodes the relevant energetics, but also makes the electronic responses of the system transparent. To derive an effective action for the fields above, we need to integrate out the fermionic degrees of freedom.

To leading order in , the structure of the long-wavelength effective action obtained from the fermionic integration is as follows:

(1)

where the actual electromagnetic vector potential is treated on the same footing as the the fields parametrizing the distortions for . As we show here, the energetics of the distortions are encoded in the mass-term coefficients , while the response-term coefficients determine the electronic responses of the system.

Focusing on an explicit large- Kondo-lattice model on the triangular lattice, Martin and Batista (2008) which forms a common nocoplanar texture, known as all-out or tetrahedral, Momoi et al. (1997), we explicitly derive the coefficients above, and discuss their implications. Our derivation parallels previous work on the effective action of classical fields coupled to Dirac fermions, Volovik and Yakovenko (1989); Yakovenko (1990); Hlousek et al. (1990); Abanov and Wiegmann (2000); Ryu et al. (2009) but is done directly on a multisite-unit-cell lattice with nonlinear dispersion. In case of a one-dimensional large- Kondo lattice model with Néel order, a nonlinear sigma model, which also incorporates fast fluctuations, has been derived by integrating out Dirac fermions. Tsvelik (1994)

The outline of the paper is as follows. In Sec. II, we briefly introduce the model of Ref. Martin and Batista, 2008 on the triangular lattice. In Sec. III, we first introduce the non-Abelian fields, which encode the smooth order-parameter distortions of the magnetic medium, and then derive the effective action through explicitly integrating out the fermions. We comment on the physical interpretation of different terms in the action, and, finally, close the paper in Sec. IV with a discussion.

Ii model and integer quantum Hall response

Consider a Kondo-lattice model in the limit of large (), where the local moments can be treated as classical. Even though there may be a direct interaction between the spins (such as nearest-neighbor Heisenberg), the most novel aspects of the physics of these itinerant systems stem from the electron-mediated interactions, originating from the Kondo coupling of the local moments to itinerant electrons. The Hamiltonian is

(2)

where is the fermion annihilation operator on site with spin , is a vector of Pauli matrices, is the electronic hopping between sites and , and is a classical magnetic moment on site (we consider magnetic moments of unit length with their large amplitude absorbed in the coupling ). The magnetic moments explore different classical configurations with energetics determined by the quantum fermionic Hamiltonian (the configuration of the magnetic moments can be thought of as external parameters in such Hamiltonians). As a function of the electronic filling and , the energetically stable configurations of the local moments determine the magnetic phase diagram of such systems. Although many interesting phases, e.g., stripes, can emerge in such systems, Akagi and Motome (2010) throughout this paper, we focus on noncopanar textures.

The Hamiltonian above may describe a system with two species of itinerant electrons with local spin moments as in Kondo lattices, or may alternatively arise in a mean-field decoupling of he Hubbard model, where represents the spin-density-wave (SDW) order parameter . In the latter case as well, the system may exhibit quantum-Hall responses, and similar energetics for smooth distortions of the SDW order parameter.

As a concrete example, consider the above Hamiltonian on the triangular lattice with nearest-neighbor electronic hopping. It has been shown that in some regions of the magnetic phase diagram, the moments form an all-out tetrahedral noncoplanar texture, which has a magnetic unit cell consisting of four sites as shown in Fig. 1Martin and Batista (2008); Akagi and Motome (2010); Kumar and van den Brink (2010); Kato et al. (2010) We represent the local moments by their components in a fixed Cartesian coordinate system. Note that because there is no spin-orbit coupling in the models we study, this frame is independent of the real-space coordinate system.

Figure 1: All-out tetrahedral order. For site in sublattice , . The four moments point from the center to the four vertices of a regular tetrahedron. The components are written explicitly in a fixed reference frame. The vectors are the lattice vectors.

Once this texture is stabilized by the electron-mediated interactions (no additional interactions are necessary in this model at and filling fraction), the itinerant electrons experience a nontrivial Berry phase, which results in a gapped integer-quantum-Hall electronic state at precisely the same and filling fractions.

As any noncoplanar structure (in this case the tetrahedral structure formed by four nearby moments on the triangular lattice) can rotate around any axis by any angle, while preserving the order (and the energy), the order-parameter space is SO corresponding to the rotation of a solid object in three dimensions. One can represent the order parameter by a unit vector (axis of rotation) and a scalar (angle of rotation) of a particular reference state (e.g., the texture shown in Fig. 1 for the triangular lattice model). A global rotation does not change the energy so the effective action should vanish for uniform and . For smooth distortions, a low-energy effective action can be obtained by a gradient expansion in and . We emphasize again the distinction of the present problem with the widely studied collinear magnets: in collinear states, the direction of the magnetic moments themselves serves as an order parameter. Since a SO rotation around the collinear moment does not change the texture, the order-parameter is the quotient SOSO(2) (a two-sphere) instead of full SO of our case (see Refs. Dombre and Read, 1989; Kawamura and Miyashita, 1984)

Before proceeding, let us comment on the limitations of our approach. Our action describes the energetics of smooth distortions of the order parameter, i.e., smooth twisting of the magnetic texture around an slowly changing axis. It describes how such distortions couple to one another as well as to the electromagnetic vector potential in the energy functional of the system, and yields various electronic responses of the system (e.g., charge and spin Hall) to such twists. However, our action does not account for all possible fluctuations of the magnetic texture. In case of the triangular-lattice model of Fig. 1, the unit cell has four sublattices, and each local moment can point in any direction (characterized by two angles). So the total number of modes for a unit cell is equal to eight. The SO order parameter is only characterized by three numbers so in addition to the space- and time-dependent orientation of the tetrahedral order, which changes smoothly to remain in the low-energy sector, there are five gapped modes that derive from the intra-unit cell magnetic distortions. Inclusion of these gapped modes can be done in analogy to the inclusion of the ferromagnetic fluctuations in the Néel order case Fradkin (1991); Wen (2004). A similar issue emerges also in the case of 120-degree order on the triangular lattice, Dombre and Read (1989) where overall magnetization fluctuations are included in addition to the SO(3) order parameter. A full theory of low-energy magnetic fluctuations in our case(tetrahedral magnetic structure in the triangular lattice) must also include the five gapped modes in addition to the order-parameter distortions. In analogy with Ref. Dombre and Read, 1989, we expect that such fast modes only couple to the smooth order-parameter distortions in the Berry-phase terms and not in the energy functional. As such, they can affect the dynamics but not the energetics of smoothly twisting the magnetic texture or the electronic responses. Similarly, deriving a full nonlinear sigma model for noncoplanar textures remains and open problem.

Iii effective action for distortions

In this section, we derive an effective action for smooth distortions of the magnetic medium around the all-out state of Fig. 1. First, we argue that certain non-Abelian gauge-like fields provide a convenient set of parameters for writing out the effective action. We then derive the action in terms of these fields by explicitly integrating out the fermions as discussed in the proceeding subsections.

iii.1 Parameterizing the distortions

In the absence of an electromagnetic external gauge field, the Hamiltonian (2) leads to the following action:

(3)

where is a two-component (for spin up and down) Grassman variable and is a smooth SO rotation matrix related to and through

(4)

where is the cross-product matrix of the unit vector , which yields the cross product when acting on any vector , and is the projection matrix of the direction. Both and are smooth functions of and . The lattice vectors represent the the bonds connecting the nearest-neighbor sites (without loss of generality, let us consider the triangular lattice model where ) and the moments correspond to the noncopanar texture stabilized in the magnetic phase under study. Note that although the rotation matrix is defined on each site (as opposed to unit cell), the constraint of smoothness of the rotation matrix restricts the possible spin orientations to smooth variations of the order parameter. In case of the triangular-lattice example, the magnetic moments take on the values shown in Fig. 1 on the four sublattices in the tetrahedral phase. The three lattice vectors are also shown in Fig. 1 for the triangular-lattice model. On each lattice site, the rotation matrix , which characterizes the distortions of a reference ordered state, is described by three independent real parameters (two for and one for ).

We now make a change of Grassman variable using an SU transformation : . The change of variable is chosen so as to absorb the rotation of the magnetic moments in the definition of new fermions Rahmani et al. (2013) In other words, we have

(5)

which leads to a term in the Lagrangian of the form . The effect of the distortions now appears in the hopping terms in the action. We can then write the action as

(6)

where is the action of a uniform noncoplanar texture in the absence of the distortions characterized by ( has the same form as Eq. (3) with the substitution and ), and comes from the inserting into the hopping terms. Notice that just like the rotation matrix , the matrix is also characterized by three real parameters.

Let us now write out corresponding to nonuniform distortions. From the first term in Eq. (3), we get the contribution to the Lagrangian . Similarly, the second term in Eq. (3) gives contributions of the form . For smooth (over length scales of the order of the lattice spacing), we can expand the above term in gradients. Up to second order, we have

(7)

where summation over repeated indices is implied.

Since , we can then write upon differentiation:

(8)

which, after substituting in Eq.(7), leads to

(9)

with the right-hand side computed at . Notice that due to the unitarity of , the second derivatives are related to the first ones and the second-order expansion (9) can be written in terms of quantities (). We thus find that the matrices provide a convenient parameterization of the smooth distortions of a uniform noncoplanar magnetic texture (such as the tetrahedral texture of Fig. 1 on the triangular lattice). With a little algebra, we can show that is related to and and their derivatives through

(10)

Therefore deriving an effective action for immediately yields an action for the original geometric fields and through substituting Eq. (10).

Let us define the coefficients of the above expansion in Pauli matrices as

(11)

where summation over is implicit. In this notation, the electromagnetic gauge fields can be added on the same footing (even though we do not have full SU symmetry) as the gauge-like fields for , which were constructed in order to characterize the distortions of the noncoplanar magnetic medium. It is important to distinguish the origin of the fields with and : while the former is a real gauge field with a well-defined gauge transformation, the latter merely resembles a gauge field. The fields for do not have their own dynamics (Maxwell term) and a gauge structure (in fact if we did have a gauge structure, these fields would be pure gauge and could be gauged away); they are just some parameters characterizing the smooth distortion of the magnetic medium around a reference noncoplanar texture. For , on the other hand, we have a Maxwell term determining the dynamics, and the standard gauge symmetry.

Putting Eqs. (11) and  (9) together, we can finally write [see Eq. (6)] as

(12)

where

(13)

Note that represents the identity matrix. Also notice that the explicit dependence of and in the above expressions all on is suppressed for brevity.

In this paper, we discuss how an effective action for can be derived by

  1. Fourier transforming the action.

  2. Integrating out the fermionic degrees of freedom.

  3. Expanding the resulting action in (justified for smooth distortions due to the presence of derivatives in the definition of the fields).

  4. Expanding in momenta as we are interested in the long-distance low-energy behavior.

We will carry out the above procedure explicitly for the triangular-lattice model. However, as can be seen from the generic structure, the form (1) and the physical implication of the effective action are not specific to this model. The value of certain nonuniversal coefficients, which we compute explicitly, however, are special properties of the triangular-lattice model.

iii.2 Full action in momentum space

In order to integrate out the fermions, it is helpful to write the action in momentum space. If the unit cell of the noncoplanar structure of interest has sites, the Hamiltonian in momentum space can be simply written as a matrix, where the factor of accounts for the two spin species. The diagonal elements of the Hamiltonian are given by , where is the local moment on sublattice in a reference texture. The off-diagonal hopping terms are given by , where is the lattice vector connecting sublattice to sublattice . In case of the triangular lattice model of Fig. 1, we have the following Hamiltonian

(14)

where the matrix is given by

(15)

in the basis given by the eight-component vector

(16)

A similar structure arises in the general case for a -component and an matrix .

To write the action, we replace the above creation operators (and similarly for the annihilation operators) by Grassman variables and endow each Grassman variable with time-dependence, i.e., . Upon Fourier transformation, we can now write as

(17)

where the inverse Green’s function and the ( in case of the triangular lattice model of Fig. 1) matrix represents the momentum-space Hamiltonian matrix [see Eq. (14)]. We can now similarly write [see Eq. (12)] in momentum space. Setting factors of volume to unity, we have , with

(18)

where is the identity matrix. In the above expression, summation over repeated indices is implied ( and are summed over and and over ). Notice that unless one of the or indices is equal to zero, the last term vanishes for as the two different Pauli matrices anticommute.

iii.3 Integrating out fermions

It is well-known that for Grassman variables and a matrix , we have . Treating the momentum and frequency dependence of the fields as matrix indices, performing the fermionic path integral over , and using the matrix identity leads to the following exact effective action

(19)

where indicates a trace over the matrix indices as well as and (we shall represent a trace over only the matrix indices by “tr”). Notice that the Green’s function is diagonal in the and indices. As mentioned before, in the limit of smooth distortions, the fields are small due to the derivatives in their definitions. Therefore, we can compute the effective action through an expansion in powers of , and truncating the expansion at a given order. Here, we go up to second order.

To obtain such expansion in , we observe that the matrix [see Eq. (18)] is comprised of first and second order terms in :

(20)

Using the above and the following expansion of the logarithm , we obtain

(21)
(22)

where the first (second) term on the right-hand side of Eq. (22) comes from the first (second) order term in the expansion of the logarithm. The equations above provide the starting point of the calculation of the effective action.

Before focusing on the triangular lattice case, where we can disentangle the spin and subsystem indices due to a symmetry, we comment on the general structure of the calculation through an example. Let us first simplify the notation by introducing , where and , and the following shorthand notation:

(23)

We now define operators , such that from Eq. (18), we have

(24)

Notice that is independent of .

As a concrete example of the computations involved, let us consider the term , which appears in [see Eq. (22)]. Each or is, in the most general case, labeled by two indices . As is diagonal in these indices, we need to compute an integral over and of the following trace of a matrix: . The frequency in such integration runs from to and the sum over discrete momenta reduces to an integral over the Brillouin zone in the thermodynamic limit. Inserting Eq. (24) into the above expression gives .

To write an action for , it is natural to make a change of variables and . The trace multiplying can then be expanded in (as we are interested in a long-wavelength action). The main calculation then involves an integration over to find the coefficients. Such calculation can be performed by diagonalizing the Hamiltonian so that the Green’s function is written as , where and are the eigenvalues and eigenvectors of the Hamiltonian. The integral over can then be exactly performed by contour integration, leaving us with just a finite momentum integral over the Brillouin zone. It then turns out that the general structure of the effective action, to the leading order we study, is given by Eq. (1). In the following section, we explicitly compute the coefficients in Eq. (1) for the triangular lattice model of Fig. 1 with tetrahedral ordering, and discuss their implications.

iii.4 Explicit calculation of the coefficients on the triangular lattice

The noncoplanar texture of Fig. 1 has a symmetry (a combination of lattice translation and spin rotation), Martin and Batista (2008) which allows us to write the Hamiltonian in a block-diagonal form (with two identical locks) using a unitary transformation [see Eq. (14) and Appendix. A]:

(25)

This transformation allows us to write the bare action as

(26)

where

(27)

Interestingly, one can check that, for all Pauli matrices , the same transformation acting on and [see Eq. (18)] gives

(28)
(29)

where the diagonal matrices are defined as

We can then cast Eq. (18) to the following form:

(30)

with :

(31)

Notice that for is independent of .

iii.4.1 First order in

As a warmup, let us start by the first-order term, which comes from Eq. (21) with the substitution [see Eq. (31)] and [see Eq. (27)] due to the transformation (25):

(32)

In deriving the above expression we have made use of the following trace identity:

(33)

The above expression for for obviously vanishes due to the zero traces of the Pauli matrices. However, for , as discussed in Appendix. B.1, we get

(34)

This expression simply corresponds to a chemical potential proportional to the scalar electromagnetic potential , which is not relevant to the distortion of the magnetic structure and therefore not included in Eq. (1).

iii.4.2 Second order in

As for the second-order term , we have two contributions according to Eq. (22): one from at second order in the expansion of the logarithm, , and one from at first order, . Using our general method, we can write these terms as follows:

(35)
(36)

where , and we have once again made use of identity (33). For the simplicity of notation, it is implied that the derivatives only act on a single proceeding term, i.e., is shorthand for . By a simple change of variables, and using , we can write as

(37)

As we are interested in the long-distance behavior, we can expand the trace in the expression above in . To zeroth order, we obtain a contribution to , while the first order term contributes to [see Eq. (1)]. The term , on the other hand has no dependence on in the trace, and therefore only contributes to .

Mass terms: Let us first focus on the mass-term coefficient . The discussion above implies that

(38)

where and are respectively the contributions of and , explicitly given by

(39)
(40)

Using and integration by parts, we can the write

(41)

Comparing with Eq. (39) indicates that for (), the two terms cancel out and the coefficient of in the effective action vanishes as expected (for instance nonzero would correspond to superconducting response). The coefficient of the mass term is then given by

(42)

We show these coefficients in Fig. 2 as a function of (for details of the calculation, see Appendix. B.2). Other coefficients not shown in Fig. 2 vanish.

The nonvanishing terms above describe the energetics of twisting the magnetic texture around the axis (as mentioned before due to the absence of spin-orbit coupling the coordinate system for spin components is independent of the coordinate system of the real-space lattice). Consider a region in real space where every local moment is rotated around the axis. If this rotation is uniform then all the spins are rotated together and there is no energetic cost with respect to the ground state. However, there is a stiffness against nonuniform rotations. Consider an example texture, where as we move in the direction, the rotation angle of the spins around the fixed axis linearly increases. This leads to a constant nonvanishing proportional to the rate of change of the rotation angle in the direction. The excess energy density (with respect to the ground state energy) required for such twist then goes as .

For a more general twist, if we decompose the rotations into rotations around three orthogonal axes, there is no cross term and the energies add up (due to the in ). If the rotation angle changes with both and , however, there are cross terms , which is expected as the underlying lattice is not symmetric with respect to the and directions. The above stiffness coefficient can be used to study the energetics of distorted configurations of interest (e.g., a state with several vortices) and compute the force between different vortices, or the force between a vortex and a boundary. The terms proportional to describe the work we need to perform to rotate the whole spin texture around the axis with a constant angular velocity (spin stiffness against time-dependent rotations).

Figure 2: The nonvanishing coefficients of the mass term as a function of . These coefficients saturate to constants for large .

Response terms: Thus far we have computed the nonquantized coefficients of the mass terms in Eq. (42), which encode the energetic cost (to second order in perturbation theory) of a twisted magnetic texture. We now turn to the

(43)

coefficients [see Eq. (1)], which originate from expanding Eq. (37) to the first order in . As we will see, these terms encode different quantized and nonquantized quantum-Hall response functions in charge and spin sectors. We start by expanding the trace in Eq. (37) in (and ). There are three terms that depend on : , , and . However, using the cyclic property of the trace and transformations and , we find that the contributions from expanding the last two terms cancel out. We then consider the expansion of , which gives a term

We can then write

(44)

The above integral can be computed as explained in Appendix. B.3. It turns out, as expected from the integer-quantum-Hall response of this system, that

(45)

Moreover, we have the following quantized coefficients in the spin sector:

(46)

which encode the transverse spin current response to a twist of the magnetic texture around the axis. Interestingly, in the spin sector , we obtain nonquantized (dependent on ) coefficients when the index in is nonzero (corresponding to position derivative of the field ). We have the following relationships between these nonquantized coefficients:

(47)

Given the relationships above, we just need to specify for as a function of . The results are summarized in Fig. 3. Notice that if a coefficient (or its real or imaginary part) is symmetric in indices and (as opposed to antisymmetric), they cancel out in the effective action. In listing the above nonvanishing coefficients, we have taken this cancellation into account.

Figure 3: The nonvanishing coefficients as a function of , which determine all other nonquantized nonzero .

The above coefficients lead to several interesting electronic spin responses. From the antisymmetry in and , we can explicitly write

(48)

Which leads to the following expectation values for spin currents:

(49)
(50)

The above expressions imply that making a time-dependent twist of the magnetic texture (around the axis) in the () direction leads to an expectation value for the spin- current in the () direction. Moreover if the velocity of spin rotation changes as we move in some real-space direction, we also get a contribution to the expectation value of the spin current.

Iv Conclusions

In this paper, using an effective-action approach, we studied smooth distortions of noncoplanar magnetic textures with interactions mediated by itinerant electrons. We argued that parameterizing such distortions in terms of certain non-Abelian fields allows us to simultaneously study the energetics of the distortions as well as the responses of electrons to such twists. The electromagnetic filed can also be readily included in our formalism. Our work constitutes a first step toward a complete theory of low-energy magnetic excitations of noncoplanar textures coupled to itinerant electrons.

The action we obtain has some mass terms, which fully characterize the energetics of twisting the magnetic texture. Moreover, it yields the expectation values of charge and spin currents for insulator electronic states coupled to the magnetic texture. Topologically stable vortices are expected to play an important role in chiral magnets, and our action can serve as a basis for studying the energetics of a collection of vortices. In particular, by providing an energy functional, our action makes it possible to apply Langevin-type simulations to study the relaxation of magnetic textures.

An extension of our formalism to include the required number of fast fluctuating modes (five modes in case of the triangular lattice with tetrahedral ordering), is expected to give rise to additive mass-like terms in the action characterizing the energetics of these modes (analogous to the ferromagnetic fluctuations in an antiferromagnet). In the Wess-Zumino terms, which contribute to the spin dynamics, however, the fast modes play a more important role.

Acknowledgements.
We thank G.-W. Chern, O. Motrunich, D. Mozyrsky, C. Nayak, Y. Nishida, and S. Ryu for enlightening discussions. This work was supported by the U. S. Department of Energy under the LANL/LDRD program (A. R, I. M, and, R. M). Work performed at Argonne National Laboratory (by I. M.) was supported by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357. R. M. also thanks CNPq (Brazil) for financial support.

Appendix A Unitary Transformation for Block Diagonalizing the Hamiltonian

The explicit form of the unitary transformation (25) is given by

(51)

It can be shown that for the Hamiltonian of Eq. (14), we have , where

(52)

with

(53)
(54)
(55)
(56)

The properties (28) and (29) can also be obtained by direct matrix multiplication.

Appendix B Calculation of the Coefficients

b.1 The first-order term

To compute Eq. (32) for , we use the method of Sec. III.3 to write , where and <