# Multicomponent Skyrmion lattices and their excitations

## Abstract

We study quantum Hall ferromagnets with a finite density topologically
charged spin textures in the presence of internal degrees of freedom such as
spin, valley, or layer indices, so that the system is parametrised by a -component complex spinor field. In the absence of anisotropies, we find
formation of a hexagonal Skyrmion lattice which completely breaks the
underlying symmetry. The ground state charge density modulation,
which inevitably exists in these lattices, vanishes *exponentially* in . We compute analytically the complete low-lying excitation spectrum,
which separates into gapless acoustic magnetic modes and a
magnetophonon. We discuss the role of effective mass anisotropy for -valley Skyrmions relevant for experiments with AlAs quantum wells. Here, we
find a transition, which breaks a six-fold rotational symmetry of a
triangular lattice, followed by a formation of a square lattice at large
values of anisotropy strength.

###### pacs:

73.43.-f, 71.10.-w, 73.43.Lp, 73.21.-b 81.05.Uw*Introduction.* Skyrmions in quantum Hall ferromagnets (1)
present an early example of topological excitations in condensed matter
physics, a topic of persistently large interest, which by now has grown to
encompass itinerant magnets with spin-orbit coupling (2); (3)
and spinor condensates (4). In the quantum Hall effect (QHE),
these topological textures carry a quantised charge, endowing them with
stability and providing possible probes via Coulomb interactions and charge
transport. The low-energy modes of these defects have possibly been observed
in NMR measurements (5). The textures are also exceptionally
tunable, as their density, and hence the relative strength of interactions,
can be modified by a gate potential or magnetic field. The possibility of
controlling spin by coupling to electron charge makes these systems
interesting candidates for spintronics applications.

In this work, we discuss the physics of Skyrmions with an enlarged internal space. This space minimally consists of two components which are, most simply, represented by physical up and down spins – but other degrees of freedom frequently play a role. Early examples are the layer index in double quantum wells of GaAs, where interlayer phase coherence spontaneously develops (6); (7), or the valley degeneracy in semiconductors (8). Recently, there has been a further proliferation, in the form of spin, layer and valley degeneracies in graphene (9), and cold atoms (4), where internal degrees of freedom and Skyrmions have emerged as interesting topics.

At finite density of defects and in presence of Coulomb interactions, Skyrmions can form a crystal (SC), analogous to Skyrme crystals of nuclear physics (10). The properties of SCs have been extensively discussed in the case, and a large number of results in a wide variety of settings have been obtained using Hartree-Fock (HF) theory for systems (11). Most of these studies are based on the numerical solution of the HF equations of motion, being hard to perform due to presence of several lengthscales.

This motivates us to develop an *analytical* approach to the general
case of Skyrmion crystals with an internal space parametrised locally by a –component complex spinor. The simplifying feature of our model is that
the spinors live in the complex projective space , which supports
exact multi-Skyrmion solutions irrespective of the value of , and, most
crucially, the assumption of large value of spin stiffness.

Our programme consists of two parts. First, we derive a set of largely analytical results on the case of Skyrmion crystals which are described locally by spinors. We find that macroscopic degeneracy, present in the noninteracting case (12), is lifted in favour of a hexagonal SC with the amplitude of charge density modulations decaying exponentially with . The excitation spectrum, obtained by studying small perturbation of the crystal, exhausts the number of gapless modes consistent with the original degeneracy. Here, of these modes correspond to generators of the Lie algebra, in accordance with the Goldstone theorem, supplemented by a magnetophonon excitation, appearing as a result of broken translational invariance. The magnetophonon has a dispersion reminiscent of a 2D Wigner crystal in a magnetic field (13).

Second, as an application of our theory, we consider explicitly a model of multi-component quantum Hall ferromagnet, motivated by experiments on semiconductor quantum wells. The effects of anisotropies, present in these systems, which break symmetry, are studied. With increasing , the number of possible types of anisotropies grows quickly along with the generators of , and we concentrate on the simplest nontrivial, yet instructive, and hitherto unstudied case of Skyrmions. These are of experimental interest for quantum wells grown in a [111] direction, in particular a nascent effort on AlAs, where one finds the appearance of three almost degenerate valleys (8). The latter are subject to a ‘nematic’ anisotropy recently discussed in the context of valley degrees of freedom (14), which we generalise to the case of an arbitrary number of valleys. We show how this anisotropy changes the nature of the crystalline ground state by combining spontaneous symmetry breaking with the externally imposed one, leading to a collection of different crystal structures. We also determine the gaps appearing in the low-energy excitation spectrum.

The remainder of the paper is organised as follows. First, we introduce the model describing spatial variations of internal degrees of freedom, taking into account interactions between electrons. An exact HF ground state in terms of spinors living in an analytic subspace is presented. Starting from this state we derive linearized equations of motion describing excitations of Skyrmion lattices. A generalization of Bloch theorem leads, in the isotropic case, to a reduction of a system of equations into decoupled pairs of linear equations, which analytical solution provides the spectrum together with the eigenfunctions of the excitations. Next, we consider a model of semiconductor, possessing three degenerate valleys, which takes into account the ‘nematic’ anisotropy, and present results for the phase diagram and excitations.

*The model.* We consider a quantum Hall ferromagnet at filling factor with internal degrees of freedom described by a -component
spinor , which corresponds to the Slater determinant
state where the spinor is normalized as at every point , creates an electron in the -th spin
component in the lowest Landau level (LLL) eigenstate with position . The latter are given in the coherent state representation, and is the electron vacuum. For example, in the case of a bilayer
quantum Hall system, components of the spinor correspond to layer and spin
indices, and in graphene to spin and valley. The Hamiltonian of the system
is given by two contributions , where is proportional to the spin stiffness , and describes Coulomb interactions. A generic form of , invariant under rotation and containing only first order
gradients reads (15); (7)

(1) |

where the integral is taken over the two dimensional plane, and the scalar product is defined as . The Hamiltonian has a local “gauge” symmetry, in other words, multiplying a spinor by an arbitrary function, , does not affect the energy. The electron interactions with a potential are described by

(2) |

in the lowest order of the gradient expansion. Here is the topological charge density, which is proportional to the electronic charge density in the quantum Hall effect (1); (15), and is related to a Pontryagin index of the normalized spinor field (12). In the following we consider textures in the form of Skyrmion crystals and fix the topological charge, given by the integral of inside the unit cell (u.c.), to . The two lattice vectors of the unit cell, which define a parallelogram with opening angle can be chosen as and , where (real/imaginary parts of correspond to components of the vectors). We measure the lengths in units of .

The energy is minimized by spinors, with entries analytic functions of , whose number of zeros inside the u.c. fixes the topological charge sector. Without interactions there is a macroscopic degeneracy coming from the arbitrariness in the choice of positions for these zeros. In presence of this degeneracy is lifted, which can lead to formation of a Skyrmion crystal. It is impossible to write an analytic expression for the ground state in the interacting case, thus we propose a variational approach in which the trial state is still expressed in terms of analytic functions . Although the Euler-Lagrange equations, derived from , are highly nonlinear, our procedure is reminiscent of the LLL projection, familiar from the theory of Abrikosov vortex lattices and the QHE physics. Our conjecture is that the zero modes of the Hessian of , spanning the analytic subspace, are separated from the higher (non-analytic) ones by a spectral gap of the order , because the average topological charge density, analogous to a magnetic field strength, is independent of . Further, the energy scale associated with the Coulomb interaction is , where is the Skyrmion density and describes the magnetic length. This scale is smaller than provided that , i.e. if the filling factor remains close to one. This justifies our variational treatment.

*The method.* We choose a basis of analytic functions defined by their
quasi-periodicity where are complex numbers and . Using the properties
of these functions under translations and , we obtain the basis of
linearly-independent theta-functions with a given topological charge inside
the u.c.

(3) |

where runs from to , in agreement with Riemann-Roch theorem (16). From the numerical minimization procedure, which was further supported by a linear stability analysis, we find that the energy minimum is achieved for a spinor formed from basis functions with equal amplitudes (up to a global rotation) , and . This corresponds to a charge distribution with the full symmetries of a hexagonal lattice, although the individual spinor components possess a lower symmetry.

In order to calculate the excitations spectra of the SC we use the basis formed from a set of quasi-periodic analytic functions obtained from by translations by vector with components and , where and and is the number of unit cells on a torus. An arbitrary analytic spinor can be written in terms of as

(4) |

see (17) and a related work (18), where an analogous approach was formulated to a case for the problem of vortex lattices in Bose-Einstein condensates. Here we introduce arbitrary complex square matrices of coefficients with linear dimension We consider time-dependent HF, with dynamics generated by a Lagrangian given as a sum of the energetic term and a Berry’s phase term, which is first order in time derivatives of the spinor components. Variation of the Lagrangian with respect to parameters gives a system of coupled equations describing the motion of a Skyrmion lattice, constrained to the analytic subspace

(5) |

Equation (5) represents one of the main results of our paper.

*LSWT and Bloch theorem.* The system of equations (5) is very
general and is not amenable to analytic treatment, mostly because the
velocity field is not explicit. However, due to
the high degree of symmetry of our reference SC state, it is possible to
obtain full analytical results at arbitrary for small perturbations of
the lattice, i.e. for linearized equations of motion. The matrices of
coefficients now can be written as a sum where are small. The quasi-periodicity of theta-functions and the
periodicity (with elementary periods and ) of
the topological charge density of the SC leads to a generalized Bloch
theorem, and the system of equations fully decouples into pairs of
linear equations. These are readily solved, resulting e.g. in the excitation
spectrum for a case presented in Fig. 2.

*Results.* The magnetophonon mode at small energies has a dispersion , with in the case of Coulomb
interactions, like the dispersion of a 2D Wigner crystal in a magnetic
field. For a short-range interaction while the spectrum is
linear at small energies, , for logarithmic (2D Coulomb)
interactions. In the large- limit, topological charge density variations
of a SC have an exponentially small amplitude, for example for a square
lattice we find

Coulomb interactions suppress charge density fluctuations, and, for a large number of components, it becomes possible to arrange for a nearly-uniform twist of the spinor. Interestingly, similar questions arise in the studies of superconductors (19).

*SU(3) case.* Let us discuss application of our general methods to a
potentially realizable case of Skyrmion crystals in AlAs quantum wells,
where the spin degree of freedom is quenched due to applied magnetic field,
and the valley-degeneracy leads to the three-component spinors with
coefficients , which we describe in zeroth order by a
symmetric theory. Effective mass anisotropy within the valleys with relative
angle generates a symmetry-breaking term in the Hamiltonian
describing a ‘nematic’ anisotropy

(6) |

where is the exchange energy in the isotropic case, and is the magnetic length. With denoting the complete elliptic integral of the first kind, and the coefficients are given by

(7) |

where and,

and stands for the anisotropy parameter. Notice that , describing repulsive interaction between the different components of the spinor. This reduces to the form derived in (14) for the case. In the isotropic limit, , the ground state, for a given filling factor, is a triangular Skyrmion lattice, which transforms into a square one with increasing value of the anisotropy strength, see Fig. 3. This transition is due to the relative scaling of Coulomb and nematic energies with . We do not find the nematic phase transition of (14) even at large values of anisotropy strength. Note that although valley polarization is favoured by the anisotropic term, its Coulomb energy cost turns out to be so large that the different spinor amplitudes in our variational ground-state remain equal, independently of the value of .

Finite gaps appear in the spectrum in presence of anisotropy. The first three modes remain gapless; they correspond to magnetophonon and to the two remaining continuous symmetry generators associated with a two-dimensional Cartan subgroup of . The other six modes are now gapped, so that the collective spectrum forms the pattern shown in Fig.3. The first two gaps scale linearly with the anisotropy strengths, while the behavior is for the other four.

*Energy scales.* In the SU(3) valley-degenerate case we obtain at small
anisotropies , and with . The anisotropy energy is much smaller than the spin stiffness scale where we assumed
for AlAs at and electron density
(20). The characteristic Coulomb energy is and can be made small compared to the stiffness by changing the
filling factor. For relatively small values of we
estimate the value of dimensionless parameter which
corresponds to a square lattice, see Fig.3.

*Conclusions*. We have studied quantum Hall ferromagnets with
finite density of topological defects (Skyrmions). In presence of Coulomb
interactions the ground state of the system is a triangular Skyrmion
lattice, which breaks all the symmetries of the internal as well as
translational symmetry, generating Goldstone modes, of
which correspond to breaking of the , and the remaining one being a
magnetophonon mode. We have explored the phase diagram of a QHE
ferromagnet, which is expected to be relevant to semiconductor
nanostructures with valley-degeneracies (8). The lattice tilts
continuously as a function of anisotropy strength and becomes square through
a phase transition at a critical anisotropy value. In this case, three of
the excitation branches remain gapless, while others acquire gaps exhibiting
different scaling with anisotropy strength.

It would be interesting to extend our theory to the case of large, but finite, spin-stiffness, and to derive an effective sigma-model for the Skyrmion lattice in the projected subspace. The difficulty which one immediately faces is that long-wavelength rotations are not preserved under projection onto the analytic subspace. Another important question is related to entanglement of internal degrees of freedom, as studied for simple Skyrmions structures in (21).

*Acknowledgements*. We are grateful to J. T. Chalker for
insightful comments. D. K. acknowledges hospitality of the Theoretical
physics department at Oxford.

### References

- S. L. Sondhi, A. Karlhede, and S. A. Kivelson, E. H. Rezayi, Phys. Rev. B 47, 16419 (1993).
- X. Z. Yu, Y. Onose, N. Kanazawa, J. H. Park, J. H. Han, Y. Matsui, N. Nagaosa and Y. Tokura, Nature 465, 901 (2010); S. Mühlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, A. Neubauer, R. Georgii, P. Böni, Science 323, 915 (2009).
- A. Bogdanov and A. Hubert, Journ. of Magn. and Magn. Mater. 138, 255 (1994); A. Bogdanov, and U. Rößler, Phys. Rev. Lett. 87, 037203 (2001). J. H. Han, J. Zang, Z. Yang, J.-H. Park, and N. Nagaosa, Phys. Rev. B 82, 094429 (2010).
- R. W. Cherng and E. Demler, Phys. Rev. A 83, 053613 (2011); Phys. Rev. A 83, 053614 (2011).
- R. Tycko, S. E. Barrett, G. Dabbagh, L. N. Pfeiffer, and K. W. West, Science 268, 1460 (1995).
- S. Q. Murphy, J. P. Eisenstein, G. S. Boebinger, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 72, 728 (1994).
- S. Girvin, “The Quantum Hall Effect: Novel Excitations and Broken Symmetries”, arXiv:cond-mat/9907002.
- S. Prabhu-Gaunkar, S. Birner, S. Dasgupta, C. Knaak, and M. Grayson, Phys. Rev. B 84, 125319 (2011); K. Eng, R. N. McFarland, B. E. Kane, Phys. Rev. Lett. 99, 016801 (2007); F. Herzog, M. Bichler, G. Koblmüller, S. Prabhu-Gaunkar, W. Zhou, M. Grayson, Appl. Phys. Lett. 100, 192106 (2012).
- Y. Sakurai and D. Yoshioka, Phys. Rev. B 85, 045108 (2012); M. O. Goerbig, Rev. Mod. Phys. 83, 1193 (2011).
- I. Klebanov, Nuclear Physics B 262, 133, (1985); T. H. R. Skyrme, Proc. Roy. Soc. 260, 127 (1961).
- R. Côté, D. B. Boisvert, J. Bourassa, M. Boissonneault, and H. A. Fertig, Phys. Rev. B 76, 125320 (2007).
- R. Rajaraman, Solitons and Instantons North-Holland, Amsterdam, 1982; A. M. Polyakov, Gauge Fields and Strings, Harwood Academic Publishers, 1987; Z. F. Ezawa, Quantum Hall Effects: Field Theoretical Approach and Related Topics, World Scientific, 2008.
- H. Fukuyama, Solid State Comm., v. 17, 10, 1323 (1975); R. Côté and A. H. MacDonald, Phys. Rev. B 44, 8759 (1991); A. G. Green, I. I. Kogan and A. M. Tsvelik, Phys. Rev. B 54, 16838 (1996).
- D. A. Abanin, S. A. Parameswaran, S. A. Kivelson, S. L. Sondhi, Phys. Rev. B 82, 035428 (2010).
- K. Moon, H. Mori, Kun Yang, S. M. Girvin, A. H. MacDonald, L. Zheng, D. Yoshioka, Shou-Cheng Zhang, Phys. Rev. B 51, 5138 (1995).
- O. Debarre, “Complex tori and Abelian varieties”, American Mathematical Society, Texts and Monographs, (1999).
- F. D. M. Haldane, E. H. Rezayi, Phys. Rev. B. 31, 2529 (1985); G. Eilenberger, Phys. Rev. 164, 628 (1967).
- S. I. Matveenko and G. V. Shlyapnikov, Phys. Rev. A 83, 033604 (2011).
- M. A. Moore, T. J. Newman, A. J. Bray, and S-K. Chin, Phys. Rev. B 58, 936 (1998).
- Y. P. Shkolnikov, S. Misra, N. C. Bishop, E. P. De Poortere, and M. Shayegan, Phys. Rev. Lett. 95, 066809 (2005).
- B. Douçot, M. O. Goerbig, P. Lederer, and R. Moessner, Phys. Rev. B 78, 195327 (2008).