## Abstract

We formulate elasticity theory with microrotations using the framework of gauge theories, which has been developed and successfully applied in various areas of gravitation and cosmology. Following this approach, we demonstrate the existence of particle-like solutions. Mathematically, this is due to the fact that our equations of motion are of sine-Gordon type and thus have soliton-type solutions. Similar to Skyrmions and Kinks in classical field theory, we can show explicitly that these solutions have a topological origin.

## 1. Introduction

### (a) Elasticity and gauge theories

In classical elasticity, one assumes that the elements of a material continuum are structureless points which are characterized by their positions only, whereas their internal properties (such as orientation, deformation, etc.) are disregarded. This picture is generalized to the concept of continua with microstructure when a material body or fluid is formed of a large number of deformable particles whose (micro)elastic properties contribute to the macroscopic dynamics of such a material medium. The corresponding mechanical model was developed very early by the Cosserat brothers (Cosserat & Cosserat 1909), who replaced a structureless material point by a point that carries a triad of axes (an oriented frame, in the modern geometrical language). This idea has led to a rich variety of models that are known by various different names such as oriented or multipolar medium, asymmetric elasticity, micropolar elasticity (Grioli 1960, 1976, 2003; Truesdell & Toupin 1960; Ericksen 1962*a*,*b*, 1967; Toupin 1962; Eringen & Suhubi 1964*a*,*b*; Green 1964; Mindlin 1964; Toupin 1964; Schaefer 1967*a*,*b*, 1970; Nowacki 1985; Capriz 1989; Eringen 1999; Dyszlewicz 2004). This field is also related to other areas of continuum mechanics, like the theory of ferromagnetic materials, cracked media, liquid crystals (De Gennes & Prost 1993), superfluids (Shankar 1977; Mermin 1979; Vollhardt & Wölfle 1990; Volovik 2003) and granular media (see the application of the Cosserat model to the description of the grain-size effects in metallic polycrystal (Zeghadi *et al.* 2005)).

In general case, the 3-frame attached to material elements is deformable, thus introducing nine new degrees of freedom in addition to the classical elasticity (three rotations, one volume expansion and five shear deformations); this is usually called a micromorphic continuum (Eringen 1976, 1999). A special case is the microstretch continuum for which the 3-frame can rotate and change volume (4 degrees of freedom: three microrotations and one microvolume expansion). Finally, in the micropolar continuum, the structured material point is rigid, thus possessing three microrotation degrees of freedom. The last is the subject of the Cosserat elasticity that considers media which can experience displacements and microrotations. In turn, this theory has two important limiting cases, one with no microrotations, which is the classical elasticity, and the second case where one assumes that the medium only experiences microrotations and no displacements. Models of the latter type have in fact a long history which can be traced back to MacCullagh (1848) (see Whittaker (1951) and Unzicker (2002)). Recently, models have been investigated in the fully nonlinear setting (Böhmer 2010; Böhmer *et al.* 2011), and plane wave-type solutions were explicitly constructed. The existence of such solutions is a highly non-trivial fact. In a linearized setting, similar solutions were investigated in Natroshvili *et al.* (2003, 2006, 2007).

Recent years have seen a revival of elasticity, in particular, in the theory of dislocations which has been analysed from a gauge-theoretical point of view. The gauge approach was developed and successfully applied in high energy physics and in the theories of gravity (Gronwald & Hehl 1993; Hehl *et al.* 1995, 2007). The gauge approach in elasticity also has a long and rich history (see, for example, Kadic & Edelen 1983; Edelen & Lagoudas 1988; Katanaev & Volovich 1992; Malyshev 2000, 2007; Katanaev 2003, 2005; Lazar 2009, 2011; Lazar & Anastassiadis 2009; Lazar & Hehl 2010). Of particular interest is the fact that the dislocations behave very much like particles. From a more formal point of view, one of the open questions in this field of research is whether or not it is possible to find soliton-type solutions which could justify the particle interpretation from a mathematical point of view. The main result of this paper is that we are able to find soliton-type solutions in a particular model of elasticity motivated by gauge theories of gravity.

### (b) Geometry and elasticity

Our notation and conventions follow the book by Hehl & Obukhov (2003) where the modern differential geometry of manifolds is presented. Here, we briefly summarize the notions and define objects needed for our current study.

We consider an elastic medium which occupies the whole of , and we can identify material points with points in space. We use the Latin alphabet to label the coordinate indices, *i*,*j*,*k*,…=1,2,3, say; we write *x*^{i} for the spatial (holonomic) coordinates of , a three-dimensional simply connected manifold which we identify with the material points of the continuum. This manifold occupies the entire three-dimensional Euclidean space . Greek indices refer to the (co)frame indices (also tangent space or anholonomic indices), we use *α*,*β*,…=1,2,3 to label the frame covectors (1 forms) which we can refer to as the distortion in elasticity. We follow the Einstein summation convention whereby we sum over twice repeated indices. We denote the frame by , such that . The *e*_{α} form the basis vectors of the tangent manifold , and ⌋ denotes the interior product.

The Euclidean metric *g* is defined on as a scalar product . As usual, the values of this scalar product on the basis vectors introduce the metric tensor *g*_{ij}=*g*(∂_{i},∂_{j}), the components of which for the Cartesian coordinates read *g*_{ij}=*δ*_{ij}=diag(+1,+1,+1). The components of the metric *g*_{αβ}=*g*(*e*_{α},*e*_{β}) in an arbitrary basis are obviously related to the frame components via . One can think of *ϑ*^{α} as a set of three mutually orthogonal and normal covectors which form a basis at any point . Each such basis covector has components (*ϑ*^{α})_{i} and thus we can also view the object as an orthogonal matrix. We can think of the metric of the deformed medium as the Cauchy–Green tensor in elasticity. Technically, the metric is used to raise and lower the indices, for example we obtain a covector-valued 1-form from a vector-valued coframe as *ϑ*_{α}=*g*_{αβ}*ϑ*^{β}, etc. The components of the coframe comprise the 3×3 matrix which is inverse to matrix composed of the frame components . From the definitions, we directly obtain the following identity .

One can draw a close analogy between gravity and elasticity; see the instructive discussion of Gronwald & Hehl (1993). The force and the hyperforce (or the couple-force), as well as the stress (force per area) and the hyperstress (couple per area), caused by the force and hyperforce, respectively, are the basic concepts in the continuum mechanics of structured media. The direct gravitational counterparts of the stress and hyperstress are the energy–momentum current and the hypermomentum (that includes spin, dilation and shear currents). The space–time itself is geometrically modelled as a continuum deformed under the action of the sources with non-trivial ‘stress’ (energy–momentum) and ‘hyperstress’ (e.g. spin current). Using this correspondence, it is possible to put both gravity theory and elasticity theory into the same differential-geometric framework in which the metric, coframe and connection are the fundamental field variables. We refer the reader for more formalism and technical detail to the comprehensive review (Hehl *et al.* 1995).

When considering deformations of an elastic material, one places certain compatibility conditions on the induced strain. These so-called Saint-Venant compatibility conditions are a form of integrability conditions so that the strain can be expressed and the symmetric derivative of the displacement. In a more geometrical language, by deforming the medium, we do not want to induce any curvature into the material. This means that one assumes the Riemann curvature tensor of the deformed medium to vanish identically. It is well known that these two views are equivalent (see Sokolnikoff 1956; Amenzade 1979). In the continuum theory of media with defects (Kröner 1981), the elastic properties of a crystal are described by means of the non-Riemannian geometries in which the notion of curvature is supplemented by a new geometrical quantity called torsion that was first introduced by Cartan (1922). Therefore, we will analyse elastic materials using the language of the Riemann–Cartan geometries with torsion and we will stick to the well-established requirement of having vanishing total curvature. The latter condition is sometimes called a second basic law in the differential geometry of media with defects (Kröner 1981), and it has the physical meaning that the crystal structure is uniquely defined everywhere, see also Kondo (1952, 1964) and Bilby *et al.* (1955). In the following subsection, we will collect the most important facts of such geometries.

### (c) Geometries with vanishing Riemann-Cartan curvature tensor

Let us describe each material point of our elastic material by a coframe *ϑ*^{α} (a vector-valued 1-form) and a connection *Γ*_{α}^{β} (a matrix-valued 1-form). Recall that the latter arises from the computation of the covariant derivative ∇_{v} along a vector field *v* of the frame vectors, ∇_{v}*e*_{β}=*Γ*_{β}^{α}(*v*) *e*_{α}, see Hehl & Obukhov (2003), §C.1. The connection is assumed to be metric-compatible which means the vanishing non-metricity *Dg*_{αβ}=0, hence the connection 1-form is skew-symmetric, *Γ*_{αβ}=−*Γ*_{βα} (where *Γ*_{αβ}=*Γ*_{α}^{γ}*g*_{βγ}). The coframe specifies the orientation of the orthonormal basis vectors at this point, while the connection determines how an arbitrary vector is parallelly transported near this point. The vanishing of the Riemann–Cartan curvature tensor (a matrix valued 2-form) takes the form
1.1
where *d* denotes the exterior derivative and ∧ is the exterior (alternating) product. In geometries, where the Riemann–Cartan curvature identically vanishes the connection is often referred to as a Weitzenböck connection. The condition *R*_{α}^{β}≡0 is known in the literature by various names like teleparallel, fernparallel of distant parallel geometry. Geometrically, this condition means that the notion of parallelism is no longer a local statement, but can be defined globally, this means on the entire manifold. A straight line connecting two material points defines an absolute notion of parallelism. Physically, the condition (1.1) has the meaning that the structure of an elastic medium is uniquely defined everywhere in , (Kondo 1952, 1964; Bilby *et al.* 1955; Kröner 1981). The only non-vanishing geometrical object in such geometries is Cartan's torsion (a vector valued 2-form) defined by
1.2

Let *ε*=*ϑ*^{1}∧*ϑ*^{2}∧*ϑ*^{3} denote the volume 3-form, then the dual forms of the products of the coframe are defined by
1.3
1.4
and
1.5
The operator *, introduced by the formulas above, is called the Hodge operator. The dual forms satisfy the following useful identities
1.6
1.7
and
1.8
The object *ε*_{αβγ} is the totally antisymmetric Levi-Civita symbol with *ε*_{123}=1.

### (d) Elastic invariants

Since *T*^{α} is the only non-trivial geometric object, it is therefore the only non-trivial object characterizing the deformations of the medium. It is natural to decompose it into its irreducible pieces which will serve as our building blocks of elastic invariants. The torsion is represented as a rank 3 tensor which is skew-symmetric in one pair of indices, *T*_{μν}^{α}=−*T*_{νμ}^{α}. It has nine independent components. Since we work in , it is also natural to consider the Hodge dual of torsion which can be regarded as a 3×3 matrix with no *a priori* symmetries. At any point we can view a 3×3 matrix as an element of a nine-dimensional vector space *V* . Performing rotations which leave the metric invariant will change the components of this matrix, thus the group SO(3) acts on *V* and we are interested in identifying invariant subspaces. Those give rise to the so-called irreducible pieces of either this matrix, or its Hodge dual, the torsion tensor. Therefore, let us define .

The three irreducible pieces of torsion are defined by
1.9
1.10
and
1.11
As we already explained, the metric is used for the raising and lowering of the indices, i.e. *ε*^{α}=*g*^{αβ}*ε*_{β} and *ϑ*_{β}=*g*_{βγ}*ϑ*^{γ}.

The irreducible pieces (1.9)–(1.11) can be conveniently understood in tensor language as the skew-symmetric part, the trace part and the traceless symmetric part of the matrix, see the discussion below and, in particular, the explicit formulas (2.8)–(2.10).

It is important to note that there is a gauge freedom in choosing the variables {*ϑ*^{α},*Γ*_{α}^{β}}. All coframes are related to each other by means of local rotations. Therefore, the same elastic medium can be described equivalently by another pair {*ϑ*^{′α},*Γ*_{α}^{′β}} which is related to the former by means of the gauge transformation
1.12
and
1.13
where *L*_{β}^{α} is an arbitrary 3×3 rotation matrix, this means *L*∈SO(3). The covariant condition *R*_{α}^{β}≡0 is preserved by this gauge transformations, whereas torsion *T*^{α} transforms covariantly, *T*^{′α}=*L*_{σ}^{α}*T*^{σ}.

Making use of the local transformations (1.12)–(1.13), we can eliminate either of the variables in the pair {*ϑ*^{α},*Γ*_{α}^{β}}. The ability to do this is closely to our assumption of a globally flat manifold (1.1). In particular, in Böhmer *et al.* (2011), the teleparallel connection is gauged away *Γ*_{α}^{β}=0, and the coframe *ϑ*^{α} is left as the only variable. Alternatively in Böhmer (2010), the coframe is gauged to the standard constant values and the only variable is the connection, in which case once can write
1.14
Here *Λ*∈SO(3) parametrizes a teleparallel connection. We will refer to this choice of variables as the connection gauge since the only dynamical variable is the connection.

## 2. Basic equations of the model

### (a) The potential energy

The natural potential energy is based on the sum of the three quadratic microrotational elastic invariants constructed from the non-vanishing torsion
2.1
where *c*_{i}, *i*=1,2,3, are non-negative constants which are physically interpreted as elastic moduli. The non-negativity of the material constants *c*_{i} is a consequence of the condition of the positive semi-definiteness of the microrotational elastic energy, *V* ≥0, in complete agreement with Lazar & Anastassiadis (2009).

In general, *V* =*V* (*T*^{α}) can depend arbitrarily on the torsion, but the choice of the quadratic function (2.1) yields the linear constitutive law (Lazar & Anastassiadis 2009; Lazar & Hehl 2010; Lazar 2011). When our medium occupies the whole of , the three quadratic terms in (2.1) are not independent because of the teleparallel condition (1.1). Recall that there is an identity (see Hehl *et al.* 1995 eqn. (5.9.18)) that is valid in *n*-dimensional space:
2.2
The tilde denotes the Riemannian geometric objects, i.e. those constructed from the Christoffel symbols of the corresponding Riemannian metric of the manifold.

Taking into account that the Riemann–Cartan curvature vanishes owing to the teleparallel constraint (1.1) and that by assumption the elastic medium is embedded in a flat Euclidean space with , for *n*=3 the identity (2.2) allows us to express the square of the first (tensor) irreducible part in terms of the trace and axial trace squares:
2.3

In tensor language, using (1.9), (1.10), and (1.6), (1.7), we have 2.4 and 2.5

It is straightforward to establish relation between the irreducible parts of the torsion 2-form *T*^{α} and its dual 1-form . We have , denote the matrix by **T**, and we find explicitly the decomposition of the second rank tensor
2.6
and
2.7
into the trace, antisymmetric and traceless symmetric parts
2.8
2.9
and
2.10

### (b) The choice of the gauge

In the recent papers (Böhmer 2010; Chervova & Vassiliev 2010; Böhmer *et al.* 2011), the propagation of waves was studied in the structured medium under consideration. The results obtained describe the elastic waves in the two complementary pictures which arise owing to the two different choices of the variables allowed by the gauge freedom (1.12)–(1.13). In Böhmer *et al.* (2011), the coframe gauge is used when the connection is zero (completely gauged out). The components of the coframe are the only dynamical variables. Geometrically, they represent an arbitrary orthogonal 3×3 matrix. Since the connection is trivial in this gauge, the torsion 2-form reduces to *T*^{α}=*dϑ*^{α}. Using the components of the coframe and frame, , , we have for the holonomic torsion
2.11
The dual 1-form then reads
2.12
Here, is the Levi-Civita tensor with respect to the holonomic coordinate basis.

A complementary approach is developed in Böhmer (2010), where the only dynamical variable is the connection, whereas the coframe is fixed to its trivial constant value . The torsion (1.2) reduces in this gauge to *T*^{α}=*Γ*_{β}^{α}∧*ϑ*^{β}, and since in the Weitzenböck space, the connection is given by (1.14), we have
2.13
Denoting , we thus have
2.14
This is completely equivalent to (2.11), with the difference that the orthogonal matrix , representing the coframe, is replaced with the orthogonal matrix , representing the connection.

In Böhmer (2010), there are two technical deviations when compared with the model developed in Böhmer *et al.* (2011). The first deviation is the different choice of the variables mentioned above. Noticing that is skew-symmetric in *i*,*j* (being an element of the Lie algebra of the orthogonal group), in Böhmer (2010) one chooses as a basic variable a new 3×3 matrix
2.15
which has no *a priori* symmetries. This ☆ denotes the dualization operator which relates skew symmetric matrices to vectors, similar to the Hodge * operator. Writing out the indices explicitly, this matrix can be written in the following way
2.16
This quantity is known in polar elasticity theory under various different names: the Nye tensor (Kröner 1967), the wryness tensor (Ericksen & Truesdell 1957), the second Cosserat deformation tensor (Eringen 1999), the third-order right micropolar curvature tensor (Neff 2004) or torsion–curvature tensor (note that this is somewhat a misnomer since it neither directly relates to torsion nor curvature in the sense of differential geometry). The inverse of (2.16) yields the components of the torsion tensor (2.14)
2.17
Substituting this into (2.14), we can find a relation between the second rank tensor **T** (which we can view as a matrix in ) to the matrix **A**, so that
2.18

The second deviation of Böhmer (2010) from Böhmer *et al.* (2011) is a different structure of the Lagrangian. Namely, the Lagrangian of a rotationally elastic medium is assumed in Böhmer (2010) to be a function of orthogonal matrix only.

### (c) Kinetic energy and total Lagrangian

The Lagrangian 4-form *V* =*L* *dt*∧*ε* is constructed by taking *L*=*L*(*u*,∂_{0}*u*,∂_{i}*u*) as a quadratic function of the irreducible parts of *A*_{lk}. In addition, the kinetic term is chosen to be
2.19
This form of the kinetic energy is well motivated, since when linearized, it yields angular kinetic energy or rotation energy. Strictly speaking, in general, the kinetic energy should be constructed by taking the microinertia into account (Maugin 1998, 2011; Eringen 1999).

Introducing, analogously to (2.16), the velocity of the deformations of the material continuum
2.20
and denoting *A*^{lt}=*g*^{lk}*A*_{kt}, we can recast (2.19) as
2.21

Taking into account the identity (2.3), the potential energy contains just two quadratic invariants. As a result, the general Lagrangian reads
2.22
The coupling constants and *λ*_{2}=*c*_{1}+*c*_{2} are constructed from the elastic moduli (2.1).

One can use different parametrizations of the orthogonal matrices. For example, in Böhmer *et al.* (2011), the spinor parametrization was used. Here, we find it more convenient to describe an arbitrary orthogonal matrix with the help of the three real functions *β*^{i} which is based on the identification of SO(3) with the three-sphere. From these independent constituents, an orthogonal matrix is constructed as follows:
2.23
Here, *β*^{i}=*δ*^{ij}*β*_{j}, *α*^{2}+*β*^{2}=1 (with *β*^{2}=*δ*_{ij}*β*^{i}*β*^{j}). The inverse matrix reads
2.24
Substituting (2.23) and (2.24) into (2.16) and (2.20), we find
2.25
and
2.26
For small *β*, the above formulas can be linearized, so that *A*_{lk}≈−2∂_{k}*β*_{l}, *A*_{lt}≈−2∂_{t}*β*_{l}, and the linearized model is described by the Lagrangian
2.27
Here as usual *divβ*=∂_{i}*β*^{i}, and (curl*β*)_{1}=∂_{2}*β*_{3}−∂_{3}*β*_{2}, etc.

### (d) The nonlinear equations of motion

The complete nonlinear equations are more non-trivial. Denote the derivatives 2.28 Then the field equations read 2.29 where we introduced 2.30 2.31 and 2.31

The equations (2.29) are valid for *any* Lagrangian *L*=*L*^{kin}(*A*_{lt})−*L*^{pot}(*A*_{lk}) with an arbitrary dependence on the variables *A*_{lt},*A*_{lk}. However, for the specific choice of the quadratic Lagrangian (2.21) and (2.22), we have explicitly
2.33

The tensor (2.30) is invertible and the inverse reads
2.34
One can easily check that . Multiplying (2.29) with this inverse, we obtain another convenient form of the fields equations:
2.35
where we have introduced
2.36
and
2.37
Note that both objects are antisymmetric in *i*,*j*.

The resulting system of nonlinear differential equations (obtained after substituting (2.33), (2.30)–(2.32) and (2.25), (2.26) into (2.29)) is quite non-trivial. It is possible, however, to find simple solutions under the additional assumptions.

## 3. Spherically-symmetric solutions: the soliton

Let us now look for the configurations with a three-dimensional spherical symmetry. The corresponding ansatz reads
3.1
where the scalar function *w*=*w*(*t*,*r*) depends on time and on the radial variable .

Then (2.25) and (2.26) yield 3.2 and 3.3 Hereafter, the dot and the prime denote the derivatives with respect to time and radius, respectively.

Accordingly, the trace and the skew-symmetric parts of (3.2) are 3.4 and 3.5 As a result, we find 3.6 and 3.7 It is straightforward to see that 3.8 and 3.9

Using these results, we find that the field equations (2.35) reduce to
3.10
where
3.11
Since the two terms in the round brackets are the Laplacian in spherical coordinates, this equation of motion can also be written in the neat form . Let us introduce a new function *φ*:=*wr*, and next let us rescale , , then (3.10) becomes
3.12
The resulting equation is closely related to the spherical ‘sine-Gordon’ equation (see Derrick 1964; Olsen & Samuelsen 1980) and references therein. It is worthwhile to mention that the model of Kontorova-Frenkel (Frenkel & Kontorowa 1938; Kontorova & Frenkel 1938; Braun & Kivshar 1998) was the first theory discussing media with dislocations where the ‘sine-Gordon’ equation and its soliton solutions naturally appeared. The main difference between previously studied spherical sine-Gordon equations and our equation (3.12) is that the nonlinearity carries an extra factor of 1/*r*. This is similar to angular momentum when studying the spherical Schrödinger equation or Newton's equation. This additional factor has some interesting implications. At large distances from the centre, this term becomes negligible and asymptotically we recover the wave equations. This analysis renders a strong support for the existence of soliton solutions in our model with a localized configuration near the centre. Below we confirm this by numerical integration.

In the static case, the equations of motion reduce to the single ODE
3.13
For *λ*_{2}=2*λ*_{1}, this further reduces to , whereas when *λ*_{2}=*λ*_{1} one is left with . The qualitative analysis of the equation (3.13) reveals the existence of static solutions that vanish at *r*=0 and approach asymptotically *π*/4 at infinity for .

Numerical integration is straightforward. The form of the solution depends on the coupling constants and on the initial value of *w*′(0). However, the qualitative behaviour remains the same. As a specific example, figure 1 presents the soliton for *λ*_{2}=*λ*_{1} when *w*′(0)=1.

In the static case, we can introduce the new function *f*, defined by
3.14
which transforms (3.10) into an autonomous second-order differential equation
3.15
where . This equation can now be analysed using standard techniques from ordinary differential equations or dynamical systems, and we find that this static system has three equilibrium points
3.16
with eigenvalues (0,−1) in all three cases which in turn yields interesting (in)stability properties.

## 4. Discussion

The rotational elasticity model has many features similar to the model suggested by Skyrme (1961). Although the Lagrangians are different, the dynamics looks qualitatively the same. In particular, the remarkable feature of the rotational elasticity is the existence of solitons that are close relatives to the Skyrmions. The topological nature of such configurations (which we demonstrate below) guarantees the stability of solutions obtained in the sense that the variations of the coframe and connection cannot change the value of the conserved charge. Moreover, the moduli space of these solutions is obviously compact.

Recently, it was noticed (Randono & Hughes 2011) that in the framework of the gauge approach in gravity and elasticity (that underlies the model under consideration, see §1*c*) one can define an identically conserved current 3-form that gives rise to the topological charge that naturally classifies the field configurations. Such a 3-current ‘lives’ on the four-dimensional space–time that is constructed as a foliation where the *t*=*const* slices coincide with the manifold . Specializing to the case of the Weitzenböck geometry with the flat curvature (1.1), this topological current reduces to
4.1
Here, is the contortion defined as the difference between the Riemann–Cartan and the Riemannian connections. This 3-form is identically conserved on the 4-space–time, d*J*^{top}≡0 in view of (1.1). As a result, we can construct the topological charge
4.2
The integral is taken over the 3-manifold . The topological nature of the charge (4.2) and current (4.1) is manifested in the fact that the corresponding conservation law does not depend on any dynamics of the fundamental fields *ϑ*^{α} and *Γ*_{α}^{β} (Randono & Hughes 2011).

The integral (4.2) is taken over the whole 3-space and the result is a constant for configurations that we studied in the previous section when *w* vanishes at the origin and approaches constant value at infinity. Specializing to the connection gauge (1.14), we have (the coframe is ‘gauged away’: ) and hence
4.3
Fixing the asymptotic boundary conditions so that approaches unity near spatial infinity allows for the standard compactification (Randono & Hughes 2011) of the space . We thus verify that the charge (4.2), (4.3) is a degree (‘winding number’) of the maps which are classified by the integer numbers that are the elements of the third homotopy group . By direct computation, we can verify that for the soliton solution described above.

Construction of a multi-soliton generalization is an interesting problem. One can approach this along the lines proposed in Randono & Hughes (2011). As a first step, we can replace the ansatz (3.1) by taking instead of (2.23) the orthogonal matrix that is a product of *N* factors
4.4
The dot denotes the usual matrix product. Such an ansatz would obviously generate a topological solution with the higher charge . The construction of the multi-soliton configurations is a more non-trivial problem which one could try to solve by using the instanton technique.

The realization of the current theoretical model in the condensed matter systems is an interesting physical problem. The geometric constraint (1.1) of the vanishing curvature rules out the disclinations. It was shown in Randono & Hughes (2011) that dislocations are also not relevant to the point-like soliton solutions. The solitons of this type are mentioned by Unzicker (2002), who calls them a ‘Shankar monopole’ (Shankar 1977) that describes a point-like defect in a *A*-phase superfluid Helium-3, (Mermin 1979; Vollhardt & Wölfle 1990; Volovik 2003). As we mentioned already, our microrotational elasticity model belongs to the class of the so-called micropolar elasticity theories. The relevant physical continua describe, in particular, liquid crystals with rigid molecules, superfluids, rigid suspensions, blood fluid with rigid cells, magnetic fluids, dust fluids, etc. Eringen (1999) (see also references therein) extensively studied micromorphic, microstretch and micropolar elasticity models with other defects beyond the dislocations and disclinations. We refer again to the paper of Randono & Hughes (2011), where the point-like topological solution (called there a torsional monopole) is depicted and a possible physical manifestations of such structures are qualitatively related to certain non-trivial electronic behaviour of condensed matter systems with strong spin-orbit coupling such as (1+3)-dimensional topological insulators. A careful quantitative analysis of such possibilities will be left for a future study.

## Acknowledgements

We thank Dmitri Vassiliev and Friedrich Hehl for discussions and advice. We also thank the referees for their valuable reports.

- Received December 8, 2011.
- Accepted January 5, 2012.

- This journal is © 2012 The Royal Society