\section{Introduction}

There exists a close analogy of geometry of spacetime in General
Relativity (GR) with a structure of defects in a crystal
\cite{1,2,3}. The crystal's defects present a special version of the
curved Riemannian spacetime - the Riemann-Cartan spacetime which
has an additional geometric property called "torsion" \cite{4}. In
the absence of matter, the world crystal is a model for Einstein's
gravitation with a new type of gauge symmetry in which a zero
torsion is a special gauge, while a zero connection (and the
absence of the Cartan curvature) is another equivalent gauge
leading to the Einstein's theory of "teleparallelism" \cite{5}. The
relation of torsion to electromagnetism \cite{6} makes photon
massive. If torsion cannot be seen experimentally, then photon is
massless. Limits on the range of magnetic fields emerging from
planets and stars give the upper bounds on the photon mass which
is extremely small: $m_{\gamma} < 10^{-27}$ eV (e.g. $m_{\gamma}<
10^{-60}$ g). Laboratory experiments yield much weaker bound:
$m_{\gamma} < 3\times 10^{-16}$ eV (e.g. $m_{\gamma} < 10^{-49}$
g). The recent critical discussion can be found in Ref.~\cite{7}.
Later in Refs.~\cite{8,9,10} the relationship between torsion and
the spin density of the gravitational field was considered.

\section{The structure of defects in a crystal}

\subsection{The structure of defects in the 3-dimensional crystal}

Crystalline defects may be generated via so called "Volterra
process", when layers or sections of matter are cut from a crystal
with a subsequent smooth rejoining of the cutting surfaces.

A crystal can have two different types of topological defects
\cite{1,2,3}, which are line-like defects in the 3-dimensional space
and form world surfaces in the 4-dimensional crystal (they may be
objects of string theory).

A first type of topological defects are translational defects
called {\it dislocations}: a single-atom layer is removed from the
crystal and the remaining atoms relax to equilibrium under the
elastic forces (see Fig.~1).

\begin{figure} [htb]
\centering
\includegraphics[height=50mm,keepaspectratio=true,angle=0]{bdlnfig1_potrait.eps}
\caption{ \it Formation of a dislocation line in a crystal.}
\end{figure}

A second type of defects is of the rotation type and called {\it
disclinations}. They arise by removing an entire wedge from the
crystal and re-gluing the free surfaces (see Fig.~2).

\begin{figure} [htb]
\centering
\includegraphics[height=45mm,keepaspectratio=true,angle=0]{bdlnfig2.eps}
\caption{ \it Formation of disclinations in a crystal.}
\end{figure}

Considering the displacement field $u_i(\vec{x})$ in the
3-dimensional crystal, we can calculate the dislocation density
given by the following tensor:
\begin{equation} \alpha_{ij} =
\epsilon_{ikl}\bigtriangledown_k\bigtriangledown_lu_j(\vec x).
                              \label{1}
\end{equation}
The local rotation field
\begin{equation}
\omega_i = \frac 12 \epsilon_{ijk}[\bigtriangledown_ju_k(\vec
x) - \bigtriangledown_ku_j(\vec x)]      \label{2} \end{equation}
determines the disclination density:
\begin{equation}
     \Theta_{ij} = \epsilon_{ikl}\bigtriangledown_k
     \bigtriangledown_l\omega_j(\vec x).
                              \label{3} \end{equation}
The defect densities satisfy the conservation laws:
\begin{equation} \bigtriangledown _i\Theta_{ij} = 0, \quad
\bigtriangledown_i\alpha_{ij} = - \epsilon_{jkl} \Theta_{kl}.
                              \label{4} \end{equation}
The plastic displacement field $u_i^p(\vec{x})$ doesn't obey the
integrability conditions of Schwarz, i.e. near the plastic
deformations the second derivatives don't commute:
\begin{equation}
   (\bigtriangledown_i\bigtriangledown_j -
    \bigtriangledown_j\bigtriangledown_i)u_k^p(\vec x) \neq 0.
                                      \label{5} \end{equation}
Eqs.~(\ref{4}) are Bianchi identities if densities $\Theta_{ij}$
and $\alpha_{ij}$ are expressed in terms of the plastic gauge
fields $\beta_{ij}^p$ and $\phi_{ij}^p$ \cite{3}:
\begin{equation}   \beta_{ij}^p \equiv \bigtriangledown_iu_j^p(\vec x)
        \label{6} \end{equation}
- for dislocation, and
\begin{equation}   \phi_{ij}^p \equiv \bigtriangledown_i\omega_j^p(\vec x)
        \label{7} \end{equation}
- for disclination. Then the defect densities
\begin{equation}
    \Theta_{ij} = \epsilon_{ikl}\bigtriangledown_k\phi_{lj}^p
    \quad {\mbox{and}}\quad \alpha_{ij} =
    \epsilon_{ikl}\bigtriangledown_k\beta_{lj}^p + \delta_{ij}
    \phi_{kk}^p - \phi_{ji}^p
                                     \label{8} \end{equation}
are invariant under the following gauge transformations:
\begin{equation} \beta_{ij}^p \to \beta_{ij}^p + \bigtriangledown_iu_j^p(\vec
x) - \epsilon_{ijk}\omega_k^p(\vec x)\quad {\mbox{and}}\quad
 \phi_{ij}^p \to  \phi_{ij}^p +  \bigtriangledown_i\omega_j^p.
                                     \label{9} \end{equation}
Here we have chosen to use the same symbols for the two gauge
functions $u_k^p(\vec{x})$ and $\omega_k^p(\vec{x})$ as we used
for the fields describing the state of the (Kleinert) crystal,
because these quantities have the same character. But it is of
course a different meaning to a gauge function that can be chosen
arbitrarily and the true field describing the crystal. Therefore,
\begin{equation}
     h_{ij} \equiv \beta_{ij}^p + \epsilon_{ijk}\omega_k^p
              \quad {\mbox{and}}\quad
      A_{ijk} \equiv \phi_{il}^p\epsilon_{ljk}
                                   \label{10} \end{equation}
are {\it translational} and {\it rotational} defect gauge fields
in the crystal.

In the Riemann-Cartan geometry of the 3-dimensional crystal the
{\it connection} is:
\begin{equation}   \Gamma_{ijk} = \bigtriangledown_i\bigtriangledown_ju_k^p,
                         \label{11} \end{equation}
the {\it torsion} is:
\begin{equation}   S_{ijk} = \frac 12(\Gamma_{ijk} - \Gamma_{jik}),
                                 \label{12} \end{equation}
and the {\it curvature} is:
\begin{equation}  R_{ijkl} = (\bigtriangledown_i\bigtriangledown_j -
    \bigtriangledown_j\bigtriangledown_i)\bigtriangledown_ku_l^p
     =  \bigtriangledown_i\Gamma_{jkl} - \bigtriangledown_j\Gamma_{ikl}.
                                    \label{13} \end{equation}
In the 3-dimensional crystal the disclination density represents
the Einstein tensor $G_{ij}$ associated with the curvature
(\ref{13}):
\begin{equation} G_{ij} = \Theta_{ij}(\vec x),
                     \label{14} \end{equation}
and the dislocation density represents the torsion:
\begin{equation}     \alpha_{ij} = \epsilon_{ikl}S_{lkj}(\vec x).
                                    \label{15} \end{equation}

\subsection{Topological defects in the 4-dimensional crystal}

The geometry of the 4-dimensional Riemann-Cartan spacetime is
described by the direct generalizations of the translational and
rotational defect gauge fields $h_{ij}$ and $A_{ijk}$ on the
vierbein field $\theta_{\mu}^I$ and the spin connection
$A_{\mu}^{IJ}$ with $I,J=(0,1,2,3)$ and $I=(0,i)$ where $i=1,2,3$
is the spatial index.

In four dimensions the indices $I$ and $J$ which function as the
``flat'' index on the vierbein may be chosen to be identified with
the directions in the Kleinert-lattice (= world-chrystal) meaning
the index $k$ on the $u_k^p(\vec{x})$, so that we may call it now
- in four dimensions - $I$ instead of $k$, and write:
\begin{equation} \theta^I_{\mu} = \partial_{\mu}u^p_I(x).  \label{15a}  \end{equation}
Similarly and analogously to the $A_{ijk}$ of equation (\ref{10})
we can in four dimensions write:
\begin{equation} A^{IJ}_{\mu} = \partial_{\mu}\omega_{IJ}. \label{15b}\end{equation}
Here we have analogously to (\ref{2}) defined
\begin{equation} \omega_{I J} =\nabla_I u_J^p(x)- \nabla_J u_I^p(x). \label{15c}
\end{equation}
In the traditional literature on gravity with spinning particles,
a special role is played by single-valued vierbein fields
$h^{\alpha}_{\mu}(x)$.

Considering the coordinate transformation:
\begin{equation}   x^{\alpha} = x^{\alpha}(x^{\mu}),
                           \label{16} \end{equation}
we go from rectilinear coordinates $x^{\alpha}$ with
$\alpha=0,1,2,3$ to arbitrary curvilinear ones $x^{\mu}$ with
$\mu=0,1,2,3$, therefore, from the flat Minkowski metric
$\eta_{\alpha\beta}=(1,-1,-1,-1)$ to the induced metric:
\begin{equation} g_{\mu\nu} = \eta_{\alpha\beta}
h_{\mu}^{\alpha}h_{\nu}^{\beta},
             \label{17} \end{equation}
where
\begin{equation}
 h_{\mu}^{\alpha}\equiv \frac{\partial x^{\alpha}}{\partial x^{\mu}}
                        \label{18} \end{equation}
are the vierbein fields (tetrads). They define local coordinate
differentials $dx^{\alpha}$ by a transformation:
\begin{equation}
 dx^{\alpha} = h^{\alpha}_{\mu}(x)dx^{\mu}.
                                    \label{18a} \end{equation}
Mathematically, the cutting and joining of the matter in the world
crystal may be described by {\it nonholonomic mappings} \cite{2} of
the next-neighbor atomic distance vectors. This mapping is not
integrable to a global coordinate transformation (\ref{16}). It is
described by a local transformation:
\begin{equation}     dx^I=\theta^I_{\mu}(x)dx^{\mu},
                        \label{19} \end{equation}
whose coefficients  $\theta^I_{\mu}(x)$ have a nonvanishing curl:
\begin{equation}  \partial_{\mu}\theta^I_{\nu} -
\partial_{\nu}\theta^I_{\mu}\neq 0,
                        \label{20} \end{equation}
what means that any coordinate transformation $x^I = x^I(x^{\mu})$
has to disobey the integrability conditions of Schwarz, and gives
the noncommutativity:
\begin{equation}   (\partial_{\mu}\partial_{\nu} -
\partial_{\nu}\partial_{\mu})x^I(x_{\mu}))\neq 0.  \label{21} \end{equation}
Therefore, now the functions $x^I(x^{\mu})$ are {\it multivalued}
\cite{2}. This is in contrast to the curl of the usual vierbein
fields $h^{\alpha}_{\mu}(x)$ which is equal to zero and determines
purely Riemannian spacetime, i.e. spacetime without torsion.

A set of new nonholonomic coordinates $dx^I$ are related to
$dx^{\alpha}$ by a multivalued Lorentz transformation:
\begin{equation}
 dx^I = \Lambda^I_{\alpha}(x)dx^{\alpha},
                              \label{21a} \end{equation}
and they are related to the physical $dx^{\mu}$ by the multivalued
tetrad fields:
\begin{equation} dx^I = \theta^I_{\mu}(x)dx^{\mu}\equiv
\Lambda^I_{\alpha}(x)h^{\alpha}_{\mu}(x)dx^{\mu}.
                            \label{21b} \end{equation}
This procedure transforms the physical laws from the flat space to
spaces with curvature and torsion.

The condition (\ref{21}) is not enough to describe all topological
defects in a crystal. Also the coefficients $\theta^I_{\mu}$ in
Eq.~(\ref{19}) themselves must violate the Schwarz conditions
leading to the noncommutative derivatives:
\begin{equation}
  (\partial_{\mu}\partial_{\nu} -
\partial_{\nu}\partial_{\mu})\theta^I_{\lambda}\neq 0.
 \label{22} \end{equation}
 %
They are called {\it maltivalued tetrads} \cite{2}. The
multivaluedness distinguishes them from the well-known
single-valued tetrads (or vierbeins) $e_{\mu}^I\equiv h_{\mu}^I$
considered in the standard literature on gravity.\\
The field strength of $A_{\mu}^{IJ}$ is:
\begin{equation}
     F_{\mu\nu}^{IJ} = \partial_{\mu}A_{\nu}^{IJ} -
         \partial_{\nu}A_{\mu}^{IJ} - [A_{\mu}, A_{\nu}]^{IJ},
                               \label{23} \end{equation}
and determines the Riemann-Cartan curvature \cite{4}:
\begin{equation}
       R_{\mu\nu\lambda\kappa} =
 \theta_{\lambda}^I \theta_{\kappa}^J
 F_{\mu\nu}^{IJ}.
                  \label{24} \end{equation}
The field strength of $\theta_{\mu}^I$ is the torsion \cite{4}:
\begin{equation}
     S_{\mu\nu}^I = D_{\mu}^{IJ}\theta_{\nu}^J - D_{\nu}^{IJ}\theta_{\mu}^J ,
                                  \label{25} \end{equation}
where
\begin{equation}   D_{\mu}^{IJ} = \delta^{IJ}\partial_{\mu} - A_{\mu}^{IJ}
                             \label{26} \end{equation}
is a covariant derivative.

\section{Plebanski formulation of gravity}

The translational and rotational crystalline defect gauge fields
-- the tetrads $\theta_{\mu}^I$ (dislocations) and the spin
connection $A_{\mu}^{IJ}$ with $I,J=(0,1,2,3)$ (disclinations) --
were used by J.~Plebanski \cite{1a} for the construction of the
gravitational action. The main idea of Plebanski's formulation of
the 4-dimensional theory of gravity was to present the
gravitational action without metrics as a product of two 2-forms
\cite{1a,2a,3a,4a,5a,6a,7a}. These 2-forms were constructed using
the connection $A^{IJ}$ and tetrads $\theta^I$ as independent
dynamical variables.  The tetrads $\theta^I$ were used instead of
the metric $g_{\mu\nu}$.

Both $A^{IJ}$ and $\theta^I$ are 1-forms:
\begin{equation}   A^{IJ} = A_{\mu}^{IJ}dx^{\mu}  \quad {\mbox{and}}\quad
       \theta^I = \theta_{\mu}^Idx^{\mu}.
                           \label{1r} \end{equation}
In Eq.~(\ref{1}) the indices $I,J = 0,1,2,3$ refer to the flat
space-time with Minkowski metric $\eta^{IJ} = {\rm diag}(1, -1,
-1, -1)$, which is tangential to the curved space with the metric
$g_{\mu\nu}$. The world interval is represented as
\begin{equation} ds^2 =  \eta_{IJ}\theta^I \otimes \theta^J,  \label{3r} \end{equation}
i.e.
\begin{equation} g_{\mu\nu} = \eta_{IJ} \theta^I_{\mu}\otimes \theta^J_{\nu}.
 \label{4r} \end{equation}
In the Plebanski's BF-theory \cite{1a}, the gravitational action
with nonzero cosmological constant $\Lambda$ is presented by the
integral:
\begin{equation}  I_{GR} = \frac{1}{\kappa^2}\int
\epsilon^{IJKL}\left(B^{IJ}\wedge
    F^{KL} + \frac{\Lambda}{4}B^{IJ}\wedge B^{KL}\right),
                                  \label{5r} \end{equation}
where $\kappa^2=8\pi G$ ($G$ is the gravitational constant).

Below we use units $\kappa=1$.

The topological sector of gravity is described by the following
integral \cite{8a}:
\begin{equation}  I_{TG} = 2 \int\left (B^{IJ}\wedge
    F^{IJ} + \frac{\Lambda}{4}B^{IJ}\wedge B^{IJ}\right) =
    \int \epsilon^{IJKL}\left(B^{IJ}\wedge F^{*KL}
    + \frac{\Lambda}{4}B^{IJ}\wedge B^{*KL}\right),
                                  \label{6r} \end{equation}
where $F^{*IJ} = \frac 12 \epsilon^{IJKL} F^{KL}$ is the dual
tensor. Here we have the topological sector only in the flat
space, and the dual tensor $F^{*IJ}$ corresponds only to the flat
indexes $I,J$.

The antisymmetric tensor $F^{IJ}$ can be split into a self-dual
component $F^+$ and an anti-self-dual component $F^-$, according
to the relation:
\begin{equation}
           F^{\pm} = \frac 12 (F \pm F^*).
                                 \label{6rr} \end{equation}

Here $B^{IJ}$ and $F^{IJ}$ are the following 2-forms:
\begin{equation}
      B^{IJ} = \theta^I\wedge \theta^J = \frac 12
      \theta_{\mu}^I\theta_{\nu}^Jdx^{\mu}\wedge dx^{\nu} ,
                    \label{7r} \end{equation}
\begin{equation}
      F^{IJ} = \frac 12 F_{\mu\nu}^{IJ}dx^{\mu}\wedge dx^{\nu}.
                    \label{8r} \end{equation}
The tensor $F_{\mu\nu}^{IJ}$ is the field strength (\ref{23}) of
the spin connection $A_{\mu}^{IJ}$.

As it was shown in Refs.~\cite{1a,2a,3a,4a,5a,6a,7a}, actions
$I_{GR}$ and $I_{TG}$, given by Eqs.~(\ref{5r}) and (\ref{6r}),
respectively, can be presented in terms of the self-dual
"left-handed" and anti-self-dual "right-handed" gravity
\cite{7a,9a}:
\begin{equation}  I_{GR} =  \int [\Sigma^i\wedge F^i - {\bar{\Sigma}}^i\wedge
{\bar F}^i + \Lambda(\Sigma^i\wedge \Sigma^i -
{\bar{\Sigma}}^i\wedge {\bar{\Sigma^i}})],
                      \label{9r} \end{equation}
and
\begin{equation}  I_{TG} =  \int [\Sigma^i\wedge F^i + {\bar{\Sigma}}^i\wedge
{\bar F}^i + \Lambda(\Sigma^i\wedge \Sigma^i +
{\bar{\Sigma}}^i\wedge {\bar{\Sigma^i}})],
                      \label{10r} \end{equation}
where $F\equiv F^+,\quad \bar F\equiv F^-,\quad \Sigma\equiv
\Sigma^+, \quad \bar{\Sigma}\equiv
 \Sigma^-,$
and the left-handed ("+") and right-handed ("--") $F^{\pm
i}_{\mu\nu}$ and $\Sigma^{\pm i}$ are given by:
\begin{equation} F^{\pm i}_{\mu\nu} =
\partial_{\mu}A^{\pm i}_{\nu} -
\partial_{\nu}A^{\pm i}_{\mu} + \epsilon^{ijk}A^{\pm j}_{\mu}A^{\pm k}_{\nu},
\label{11r} \end{equation}
\begin{equation}
    - \Sigma^{\pm i} = - i \theta^0\wedge \theta^i \pm \frac 12
     \epsilon^{ijk}\theta^j\wedge \theta^k.
                          \label{12r} \end{equation}
Choosing the correct gauge in his gravitational actions, Plebansky
introduced  the Lagrange multipliers $\psi_{ij}$, which are
considered in theory as an auxiliary fields, symmetric and
traceless.

Finally, the resulting actions of the Plebanski gravity are
\cite{1a,2a,3a,4a,5a,6a}:
\begin{equation}  I_{GR} =  \int [\Sigma^i\wedge F^i - {\bar{\Sigma}}^i\wedge
{\bar F}^i + (\Psi^{-1})_{ij}(\Sigma^i\wedge \Sigma^j -
{\bar{\Sigma}}^i\wedge {\bar{\Sigma^j}})],
                      \label{13r} \end{equation}
and
\begin{equation}  I_{TG} = \int [\Sigma^i\wedge F^i + {\bar{\Sigma}}^i\wedge
{\bar F}^i + (\Psi^{-1})_{ij}(\Sigma^i\wedge \Sigma^j +
{\bar{\Sigma}}^i\wedge {\bar{\Sigma^j}})],
                      \label{14r} \end{equation}
where
\begin{equation}   (\Psi^{-1})_{ij} = \Lambda\delta_{ij} + \psi_{ij}.
                     \label{15r} \end{equation}
The explanation of the introduction of the auxiliary field
$\psi_{ij}$ is connected with the condition for metricity. The
'area metric' (see \cite{9a}) is:
\begin{equation}    m = \frac 12 B^{IJ}\otimes B^{IJ},      \label{49} \end{equation}
which can be written in terms of $\Sigma^{\pm i}$:
\begin{equation}
   m^{\pm} = \Sigma^{\pm i}\otimes \Sigma^{\pm i} - \frac 14 V^{\pm}
                       \label{50} \end{equation}
with
\begin{equation}  V^{\pm} = \pm 4\Sigma^{\pm i}\wedge \Sigma^{\pm i},
 \label{51} \end{equation}
where $ V^{\pm}$ both are equal to the 4-volume form $V= \theta^0
\wedge \theta^1 \wedge \theta^2 \wedge \theta^3$. Here $m^+$ and
$V^+$ define a "left" geometry and $m^-$ and $V^-$ define a
"right" geometry of spacetime. These geometries are defined, but
they can be in disagreement. The stationarity with respect to
$\psi_{ij}$ requires the equality $m^+=m^-$ of the left and right
area metrics $m^{\pm}$, which have no totally antisymmetric part,
and leaves the totally antisymmetric parts of $\Sigma^{\pm
i}\otimes \Sigma^{\pm i}$, which form $V$, unconstrained. We see
that $m^+ - m^-$ has 20 linearly independent components. The
equality $m^+=m^-$ provides the correct number of constraints,
reducing the 36 degrees of freedom of ($\Sigma^+,\Sigma^-$) to the
16 degrees of freedom of tetrads $\theta_{\mu}^I$.

The main assumption of Plebanski was that our world, in which we
live, is a self-dual left-handed gravitational world described by
the action:
\begin{equation}
      I_{(selfdual\,\, GR)}(\Sigma,A) =  \int [\Sigma^i\wedge F^i +
 (\Psi^{-1})_{ij}\Sigma^i\wedge \Sigma^j] ,
                      \label{52} \end{equation}
and the anti-self-dual right-handed gravitational world is absent
in Nature ($\bar F=0$ and $\bar \Sigma = 0$), and we have the
equality \cite{7a}:
\begin{equation}     I_{GR} = I_{TG},
 \label{18r} \end{equation}
i.e. the gravitational sector of gravity coincides with
topological phase of gravity.

Postulating the existence of the Mirror (or Hidden) World
\cite{11,12,13,14,15,16}, we assumed in Ref.~\cite{7a}, that the
anti-self-dual right-handed gravity is the "mirror gravity" given
by the equation:
\begin{equation}
 I_{(antiselfdual\,\, GR)}(\bar \Sigma, \bar A) = \int [{\bar{\Sigma}}^i\wedge {\bar F}^i
 + (\Psi^{-1})_{ij}{\bar{\Sigma}}^i\wedge {\bar{\Sigma^j}}].
                      \label{53} \end{equation}
This "mirror gravity" describes the gravity in the Mirror World.

\section{Torsion}

The gravitational theory with torsion can be presented by the
following integral:
\begin{equation} I_S = 2 \int (2 S^I\wedge S^I + \frac{\Lambda}{4}B^{IJ}\wedge
B^{IJ}),  \label{54} \end{equation}
where the 2-form
\begin{equation}  S^I =\frac 12 S_{\mu\nu}^Idx^{\mu}\wedge dx^{\nu}
            \label{19r} \end{equation}
contains the torsion $S_{\mu\nu}^I$, which is the field strength
(\ref{25}) of the tetrads $\theta_{\mu}^I$.

Using the partial integration and putting
\begin{equation} \int \partial_{\mu}(T_{\nu\kappa\lambda})dx^{\mu}\wedge
dx^{\nu}\wedge dx^{\kappa}\wedge dx^{\lambda}=0,
 \label{16x} \end{equation}
it is not difficult to show that
\begin{equation} \int B^{IJ}\wedge
    F^{IJ} = 2\int S^I\wedge S^I.  \label{17x} \end{equation}
According to Eqs.~(\ref{6r}) and (\ref{54}), we have:
\begin{equation}
 I_{TG} = 2 \int \left(B^{IJ}\wedge
    F^{IJ} + \frac{\Lambda}{4}B^{IJ}\wedge B^{IJ}\right) =
2 \int \left(2 S^I\wedge S^I + \frac{\Lambda}{4}B^{IJ}\wedge
B^{IJ}\right) = I_S. \label{18x} \end{equation}
This means that in the self-dual Plebanski formulation of
gravitational theory (i.e. in our world) all sectors of gravity
coincide \cite{7a}:
\begin{equation}
  I_{GR} = I_{TG} = I_S.
 \label{55} \end{equation}
Here we see a close analogy of the  geometry of spacetime in GR
with a structure of defects in a crystal \cite{1,2,3}. Torsion is
presented by the translational defects --dislocations, and
curvature is given by the rotational defects -- disclinations. But
these defects are not independent of each other: a dislocation is
equivalent to a disclination-antidisclination pair, and a
disclination gives a string of dislocations.

The crystalline defects present a special version of the curved
spacetime - the Riemann-Cartan spacetime which has an additional
geometric property - "torsion" \cite{4}, described by dislocations.
The world crystal is a model for the Einstein's gravitation which
has a new type of gauge symmetry with a zero torsion as a special
gauge, while a zero connection (with zero Cartan curvature) is
another equivalent gauge with nonzero torsion which corresponds to
the Einstein's theory of "teleparallelism" \cite{5}. But any
intermediate choice of the gauge for the field $A^{IJ}_{\mu}$ is
also allowed \cite{3}. In the present investigation we have shown
that in the Plebanski formulation the phase of gravity with
torsion is equivalent to the ordinary or topological gravity, and
we can exclude a torsion as a separate dynamical variable. This
explains why the Einstein's theory of gravity described only by
curvature corresponding to the disclination defects of a crystal
can be rewritten as a teleparallel theory of gravity described
only by torsion which corresponds to the dislocation defects.

\section{Summary and Conclusions}

In this paper we have shown the close analogy of geometry of
space-time in GR with a structure of defects in a crystal, given
by H.~Kleinert in his model \cite{1,2,3}. This model was related
with the Plebanski formulation of gravity. In the Section 1 we
have discussed a role of torsion. We have reviewed the Kleinert's
model of a crystal with defects in the Section 2 considering the
translational defects -- dislocations, and the rotational defects
-- disclinations -- in the 3- and 4-dimensional crystals. In the
Section 3 we have explained the main idea of Plebanski \cite{1a} to
construct the 4-dimensional gravitational action without metric.
The integrand of this action contains a product of two 2-forms,
which are constructed from the tetrads $\theta^I$ and the
connection $A^{IJ}$ considered as independent dynamical variables.
Here we have presented the relation between the Kleinert's model
of a crystal with defects and Plebanski's formulation of gravity.
It was shown that the tetrads $\theta^I_{\mu}$ are dislocation
gauge fields, and the connections $A^{IJ}_{\mu}$ are the
disclination gauge fields, the curvature is the field strength of
connections $A^{IJ}_{\mu}$, and the torsion is the field strength
of tetrads $\theta^I_{\mu}$. Both $A^{IJ}$ and $\theta^I$ are
1-forms. The tetrads $\theta^I_{\mu}$ were used instead of the
metric $g_{\mu\nu}$.

Section 4 is devoted to the torsion. The crystalline defects
represent a special version of the curved space-time -- the
Riemann--Cartan space-time with torsion \cite{4}. The world crystal
gives a model with a new type of gauge symmetry for GR \cite{3}. The
Einstein's gravitation, which has a zero torsion, corresponds to a
special gauge, while the theory with a zero connection (i.e. zero
Cartan curvature) is another equivalent gauge with nonzero torsion
which corresponds to the Einstein's theory of "teleparallelism"
\cite{5}. In this paper we have showed that in the Plebanski
formulation of gravity the torsional gravitational phase (phase
described only by nonzero torsion) is equivalent to the ordinary
or topological gravity, and we can exclude torsion as a separate
dynamical variable of gravity. In the Plebanski formulation of
4-dimensional gravity the Einstein's gravity with a space-time
corresponding to the structure of a crystal with disclinations,
i.e. rotational defects, coincides with the Einstein's
"teleparallel" gravitational theory \cite{5} corresponding to a
crystal having only the translational defects -- dislocations.

\section{Acknowledgments}

H.B.N. thanks for being alluded hospitality and office to finish
this work. L.V.L. thanks the Niels Bohr Institute for the
wonderful hospitality and financial support. The authors deeply
thank M.~Chaichian and A.~Tureanu from Helsinki University for
useful discussions and help. CRD acknowledges a scholarship from
the Funda\c c\~ ao para a Ci\^ encia e a Tecnologia (FCT, Portugal)
(ref.~SFRH/BPD/41091/2007) and this work was
partially supported by the FCT projects CERN/FP/123580/2011,
PTDC/FIS/098188/2008 and CFTP-FCT Unit 777 which are partially funded
through POCTI (FEDER).


\section{Classical M\"obius transformations and\\ Poincar\'e distance on the quaternionic unit ball} \label{classicM}

A classical topic in quaternionic analysis is the study of M\"obius transformations. It is well known that the set of \emph{linear fractional transformations} of the extended quaternionic space ${\mathbb{H}} \cup \{\infty\} \cong \mathbb{HP}^1$
\begin{equation}\label{IndianaM}
\mathbb{G}=\left\{g(q)=(aq+b)\cdot (cq+d)^{-1} \ \left|\ \begin{bmatrix}
a & b \\
c & d
\end{bmatrix} \in GL(2,{\mathbb{H}}) \right\}\right.
\end{equation}
is a group with respect to the composition operation. We recall that $GL(2,{\mathbb{H}})$ denotes the group of $2 \times 2$ invertible quaternionic matrices, and that $SL(2,{\mathbb{H}})$ denotes the subgroup of those such matrices which have unit Dieudonn\'e determinant (for details, see \cite{poincare} and references therein). It is known in literature that $\mathbb{G}$ is isomorphic to $PSL(2,{\mathbb{H}})= SL(2,{\mathbb{H}}) / \{\pm {\mathrm{Id}}\}$ and that all of its elements are conformal maps. Among the works that treat this matter, even in the more general context of Clifford Algebras, let us mention \cite{ahlforsmoebius, maass, vahlen}.
The alternative representation
\begin{equation}\label{InvM}
\mathbb{G}=\left\{F_A(q)=(qc+d)^{-1} \cdot (qa+b) \ \left|\ A=\begin{bmatrix}
a & c \\
b & d
\end{bmatrix} \in GL(2,{\mathbb{H}}) \right\}\right.
\end{equation}
is also possible, and the anti-homomorphism $GL(2,{\mathbb{H}}) \to \mathbb{G}$ mapping $A$ to $F_A$ induces an anti-isomorphism between $PSL(2,{\mathbb{H}})$ and $\mathbb{G}$.
The group $\mathbb{G}$ is generated by the following four types of transformations:
\begin{itemize}
\item[(i)] $L_1(q)=q+b,\,\,\, b \in \mathbb{H};$
\item[(ii)] $L_2(q)=q \cdot a, \,\,\, a \in \mathbb{H}, |a|=1;$
\item[(iii)] $L_3(q)= r \cdot q = q \cdot r, \,\,\, r \in \mathbb{R}^+ \setminus \{ 0 \};$
\item[(iv)] $L_4(q)=q^{-1}.$
\end{itemize}
Moreover, if $\mathcal{S}_i$ is the family of all real $i$-dimensional spheres, if $\mathcal{P}_i$ is the family of all real $i$-dimensional affine subspaces of $\mathbb{H}$ and if $\mathcal{F}_i=\mathcal{S}_i \cup \mathcal{P}_i$ then $\mathbb{G}$ maps elements of $\mathcal{F}_i$ onto elements of $\mathcal{F}_i,$ for $i=3,2,1$. At this regard, see \cite{wilker}; detailed proofs of all these facts can be found in \cite{poincare}.

The subgroup $\mathbb{M}\leq\mathbb{G}$ of \emph{(classical) M\"obius transformations} mapping the quaternionic open unit ball
$$\mathbb{B}=\{q \in {\mathbb{H}}\ |\  |q|<1\}$$
onto itself, has also been studied in detail. Let us denote $H =
\begin{bmatrix}
1 & 0 \\
0 & -1
\end{bmatrix}$ and
$$Sp(1,1) = \{C \in GL(2, {\mathbb{H}}) \ |\  \overline C ^t H C = H\}\subset SL(2, {\mathbb{H}}).$$

\begin{theorem}
An element $g \in \mathbb{G}$ is a classical M\"obius transformation of ${\mathbb{B}}$ if and only if $g(q)=(qc+d)^{-1} \cdot (qa+b)$ with
$\begin{bmatrix}
a & c\\
b & d
\end{bmatrix}\in Sp(1,1)$. This is equivalent to
\begin{equation*}
g(q) = v^{-1} (1-q\bar q_0)^{-1}(q-q_0)u
\end{equation*}
for some $u,v \in \partial {\mathbb{B}}, q_0 \in {\mathbb{B}}$.
\end{theorem}

For a proof, see \cite{poincare}. Hence, $\mathbb{M}$ is anti-isomorphic to $Sp(1,1)/\{\pm {\mathrm{Id}}\}$.
Since $\mathbb{G}$ leaves invariant the family $\mathcal{F}_1$ of circles and affine lines of $\mathbb{H}$, and since the elements of $\mathbb{G}$ are conformal, the group $\mathbb{M}$ of classical M\"obius transformations of ${\mathbb{B}}$ preserves the following class of curves.
\begin{definition}\label{non-euclid}
If $q_1 \neq q_2 \in {\mathbb{B}}$ are $\mathbb{R}-$linearly dependent, then the diameter of ${\mathbb{B}}$ through $q_1,q_2$ is called the \emph{non-Euclidean line} through $q_1$ and $q_2$. Otherwise, the \emph{non-Euclidean line} through $q_1$ and $q_2$ is defined as the unique circle through $q_1,q_2$ that intersects $\partial \mathbb{B}= {\mathbb{S}}^3$ orthogonally.
\end{definition}

\begin{theorem}
The formula
\begin{equation}\label{poincaredist}
\delta_{\mathbb{B}}(q_1,q_2)
=\frac{1}{2}\log\left(\frac{1+|1-q_1\overline{q_2}|^{-1}|q_1-q_2|}{1-|1-q_1\overline{q_2}|^{-1}|q_1-q_2|}\right)
\end{equation}
(for $q_1,q_2 \in{\mathbb{B}}$) defines a distance that has the non-Euclidean lines as its geodesics. The elements of $\mathbb{M}$ and the map $q \mapsto \overline{q}$ are all isometries for $\delta_{\mathbb{B}}$.
\end{theorem}

We refer the reader to \cite{ahlfors1988}; a detailed presentation can be found in \cite{poincare}.

So far, we mentioned properties of the classical M\"obius transformations that are completely analogous to the complex case. However, the analogy fails on one crucial point. The group $\mathbb{M}$ is not included in the best known analog of the class of holomorphic functions: the set of Fueter regular functions, i.e., the kernel of $\frac{\partial}{\partial x_0}+i\frac{\partial}{\partial x_1}+ j\frac{\partial}{\partial x_2}+k\frac{\partial}{\partial x_3}$ (see \cite{sudbery}). For instance, none of the rotations $q \mapsto aq$ with $a \in {\mathbb{H}}, a \neq 0$ is Fueter regular, nor are any of the transformations listed as (i),(ii),(iii),(iv) in our previous discussion. The variant of the Fueter class considered in \cite{perotti2009}, defined as the kernel of $\frac{\partial}{\partial x_0}+i\frac{\partial}{\partial x_1}+ j\frac{\partial}{\partial x_2}-k\frac{\partial}{\partial x_3}$, includes part of the group, for instance the rotations $q \mapsto qb$ for $b \in {\mathbb{H}}, b \neq 0$, but not all of it (for instance $q \mapsto kq$ is not in the kernel, nor is $q \mapsto q^{-1}$).

A more recent theory of quaternionic functions, introduced in \cite{cras,advances}, has proven to be more comprehensive. The theory is based on the next definition.

\begin{definition}
Let $\Omega$ be a domain in ${\mathbb{H}}$ and let $f : \Omega \to {\mathbb{H}}$ be a function. For all $I \in {\mathbb{S}} = \{q \in {\mathbb{H}} \ |\  q^2 = -1\}$, let us denote $L_I = {\mathbb{R}} + I {\mathbb{R}}$, $\Omega_I = \Omega \cap L_I$ and $f_I = f_{|_{\Omega_I}}$.
The function $f$ is called (Cullen or) \emph{slice regular} if, for all $I \in {\mathbb{S}}$, the function $\bar \partial_I f : \Omega_I \to {\mathbb{H}}$ defined by
$$
\bar \partial_I f (x+Iy) = \frac{1}{2} \left( \frac{\partial}{\partial x}+I\frac{\partial}{\partial y} \right) f_I (x+Iy)
$$
vanishes identically.
\end{definition}

By direct computation, the class of slice regular functions includes all of the generators we listed as (i),(ii),(iii),(iv). It does not contain the whole group $\mathbb{G}$ (nor its subgroup $\mathbb{M}$), because composition does not, in general, preserve slice regularity. However, \cite{moebius} introduced the new classes of \emph{(slice) regular fractional transformations} and \emph{(slice) regular M\"obius transformations} of ${\mathbb{B}}$, which are nicely related to the classical linear fractional transformations. They are presented in detail in sections \ref{regularfractional} and \ref{regularmoebius}.

One of the purposes of the present paper is, in fact, to compare the slice regular fractional transformations with the classical ones. Furthermore, we take a first glance at the role played by slice regular M\"obius transformations in the geometry of ${\mathbb{B}}$. In section \ref{differential}, we undertake a first study of their differential properties: we prove that they are not, in general, conformal, and that they do not preserve the Poincar\'e distance $\delta_{\mathbb{B}}$. In section \ref{schwarzpick}, we overview the main results of \cite{schwarzpick}. Among them is a quaternionic analog of the Schwarz-Pick lemma, which discloses the possibility of using slice regular functions in the study of the intrinsic geometry of ${\mathbb{B}}$.

\section{Regular fractional transformations}\label{regularfractional}

This section surveys the algebraic structure of slice regular functions, and its application to the construction of regular fractional transformations. From now on, we will omit the term `slice' and refer to these functions as regular, {\it tout court}.
Since we will be interested only in regular functions on Euclidean balls
$$B(0,R) = \{q \in {\mathbb{H}} \ |\  |q| <R\},$$
or on the whole space ${\mathbb{H}} = B(0,+\infty)$, we will follow the presentation of \cite{zeros,poli}. However, we point out that many of the results we are about to mention have been generalized to a larger class of domains in \cite{advancesrevised}.

\begin{theorem}
Fix $R$ with $0<R\leq + \infty$ and let $\mathcal{D}_R$ be the set of regular functions $f:B(0,R)\to {\mathbb{H}}$. Then $\mathcal{D}_R$  coincides with the set of quaternionic power series $f(q) =\sum_{n \in {\mathbb{N}}} q^n a_n$ (with $a_n \in {\mathbb{H}}$) converging in $B(0,R)$. Moreover, defining the \emph{regular multiplication} $*$ by the formula
\begin{equation}
\left(\sum_{n \in {\mathbb{N}}} q^n a_n\right) *  \left(\sum_{n \in {\mathbb{N}}} q^n b_n\right)= \sum_{n \in {\mathbb{N}}} q^n \sum_{k=0}^n a_k b_{n-k},\end{equation}
we conclude that $\mathcal{D}_R$ is an associative real algebra with respect to $+,*$.
\end{theorem}

The ring $\mathcal{D}_R$ admits a classical ring of quotients
$$\mathcal{L}_R = \{f^{-*}*g \ |\  f,g \in \mathcal{D}_R, f \not \equiv 0\}.$$
In order to introduce it, we begin with the following definition.

\begin{definition}\label{conjugate}
Let $f(q) = \sum_{n \in {\mathbb{N}}} q^n a_n$ be a regular function on an open ball $B = B(0,R)$. The \textnormal{regular conjugate} of $f$, $f^c : B \to {\mathbb{H}}$, is defined as $f^c(q) = \sum_{n \in {\mathbb{N}}} q^n \bar a_n$ and the \textnormal{symmetrization} of $f$, as $f^s = f * f^c = f^c*f$.
\end{definition}

Notice that $f^s(q) = \sum_{n \in {\mathbb{N}}} q^n r_n$ with $r_n = \sum_{k = 0}^n a_k \bar a_{n-k} \in {\mathbb{R}}$. Moreover, the zero-sets of $f^c$ and $f^s$ have been fully characterized.

\begin{theorem}\label{conjugatezeros}
Let $f$ be a regular function on $B = B(0,R)$. For all $x,y \in {\mathbb{R}}$ with $x+y{\mathbb{S}} \subseteq B$, the regular conjugate $f^c$ has as many zeros as $f$ in $x+y{\mathbb{S}}$. Moreover, the zero set of the symmetrization $f^s$ is the union of the $x+y{\mathbb{S}}$ on which $f$ has a zero.
\end{theorem}

We are now ready for the definition of regular quotient. We denote by
$${\mathcal{Z}}_h = \{q \in B \ |\  h(q) = 0\}$$
the zero-set of a function $h$.

\begin{definition}\label{quotient}
Let $f,g : B = B(0,R) \to {\mathbb{H}}$ be regular functions. The \emph{left regular quotient} of $f$ and $g$ is the function $f^{-*} * g$ defined in $B \setminus {\mathcal{Z}}_{f^s}$ by
\begin{equation}
f^{-*} * g (q) = {f^s(q)}^{-1} f^c * g(q).
\end{equation}
Moreover, the \emph{regular reciprocal} of $f$ is the function $f^{-*} = f^{-*} * 1$.
\end{definition}

Left regular quotients proved to be regular in their domains of definition.
If  we set $(f^{-*}*g)*(h^{-*}*k) = (f^{s}h^{s})^{-1} f^c*g*h^c*k$ then $(\mathcal{L}_R,+,*)$ is a division algebra over $\mathbb{R}$ and it is the classical ring of quotients of $(\mathcal{D}_R,+,*)$ (see \cite{rowen}). In particular, $\mathcal{L}_R$ coincides with the set of \emph{right regular quotients}
$$g*h^{-*} (q) = {h^s(q)}^{-1} g * h^c(q).$$
The definition of regular conjugation and symmetrization is extended to $\mathcal{L}_R$ setting $(f^{-*}*g)^c = g^c*(f^c)^{-*}$ and $(f^{-*}*g)^s(q) = {f^s(q)}^{-1}g^s(q)$.
Furthermore, the following relation between the left regular quotient $f^{-*} * g(q)$ and the quotient $f(q)^{-1} g (q)$ holds.

\begin{theorem}\label{quotients}
Let $f,g$ be regular functions on $B=B(0,R)$. Then
\begin{equation}
f*g(q) = f(q) g(f(q)^{-1}qf(q)),
\end{equation}
and setting $T_f(q) = f^c(q)^{-1} q f^c(q)$ for all $q \in B \setminus {\mathcal{Z}}_{f^s}$,
\begin{equation}
f^{-*}*g(q) = f(T_f(q))^{-1} g(T_f(q)),
\end{equation}
for all $q \in B \setminus {\mathcal{Z}}_{f^s}$. For all $x,y \in {\mathbb{R}}$ with $x+y{\mathbb{S}}\subset B \setminus {\mathcal{Z}}_{f^s}$, the function $T_f$ maps $x+y{\mathbb{S}}$ to itself (in particular $T_f(x) = x$ for all $x \in {\mathbb{R}}$). Furthermore, $T_f$ is a diffeomorphism from $B \setminus {\mathcal{Z}}_{f^s}$ onto itself, with inverse $T_{f^c}$.
\end{theorem}

We point out that, so far, no simple result relating $g*h^{-*}(q)$ to $g(q) h(q)^{-1}$ is known.

This machinery allowed the introduction in \cite{moebius} of regular analogs of linear fractional transformations, and of M\"obius transformations of ${\mathbb{B}}$. To each
$A=\begin{bmatrix}
a & c \\
b & d
\end{bmatrix} \in GL(2, {\mathbb{H}})$ we can associate the \emph{regular fractional transformation}
$${\mathcal{F}}_A (q) = (qc+d)^{-*}*(qa+b).$$
By the formula $(qc+d)^{-*}*(qa+b)$ we denote the aforementioned left regular quotient $f^{-*}*g$ of $f(q) = qc+d$ and $g(q) = qa +b$. We denote the $2 \times 2$ identity matrix as ${\mathrm{Id}}$.

\begin{theorem}
Choose $R>0$ and consider the ring of quotients of regular quaternionic functions in $B(0,R)$, denoted by $\mathcal{L}_R$.
Setting
\begin{equation}
f.A = (f c + d)^{-*}*(f a+ b)
\end{equation}
for all $f \in \mathcal{L}_R$ and for all $A =
\begin{bmatrix}
a & c \\
b & d
\end{bmatrix} \in GL(2,{\mathbb{H}})$
defines a right action of $GL(2,{\mathbb{H}})$ on $\mathcal{L}_R$. A left action of $GL(2,{\mathbb{H}})$ on $\mathcal{L}_R$ is defined setting
\begin{equation}
^{t}A.f=(a*f+b)*(c*f+d)^{-*}.
\end{equation}
The stabilizer of any element of $\mathcal{L}_R$ with respect to either action includes the normal subgroup $N=\left \{t \cdot {\mathrm{Id}}\ |\  t \in {\mathbb{R}} \setminus \{0\} \right \}\unlhd GL(2,{\mathbb{H}})$. Both actions are faithful when reduced to $PSL(2,{\mathbb{H}}) = GL(2,{\mathbb{H}})/N$.
\end{theorem}

The statements concerning the right action are proven in \cite{moebius}, the others can be similarly derived. The two actions coincide in one special case.

\begin{proposition}\label{hermitian}
For all Hermitian matrices $A = \begin{bmatrix}
a & \bar b \\
b & d
\end{bmatrix}$ with $a,d \in {\mathbb{R}}, b \in {\mathbb{H}}$,
$$f.A = (f \bar b + d)^{-*}*(f a+ b) =(a*f+b)*(\bar b*f+d)^{-*}= A^{t}.f$$
\end{proposition}

\begin{proof}
We observe that
$$(f \bar b + d)^{-*}*(f a+ b) =(a*f+b)*(\bar b*f+d)^{-*}$$
if, and only if,
$$(f a+ b)*(\bar b*f+d) =(f \bar b + d)*(a*f+b),$$
which is equivalent to
$$a f* \bar b*f + |b|^2 f+ ad f + db = a f* \bar b*f +ad f+ |b|^2 f +db.$$
\end{proof}

In general, a more subtle relation holds between the two actions.

\begin{remark}\label{leftright}
For all $A \in GL(2,{\mathbb{H}})$ and for all $f\in \mathcal{L}_R$, by direct computation
\begin{equation*}
\left(f.A\right)^c=\bar{A}^t.f^c
\end{equation*}
As a a consequence, if $A\in GL(2, {\mathbb{H}})$ is Hermitian then $(f.A)^c= f^c.\bar{A}$.
\end{remark}

Interestingly, neither action is free, not even when reduced to $PSL(2,{\mathbb{H}})$. Indeed, the stabilizer of the identity function with respect to either action of $GL(2,{\mathbb{H}})$ equals
$$\{c\cdot {\mathrm{Id}}\ |\  c \in {\mathbb{H}} \setminus \{0\} \},$$
a subgroup of $GL(2,{\mathbb{H}})$ that strictly includes $N$ and is not normal. As a consequence, the set of regular fractional transformations
$${\mathfrak{G}} = \{{\mathcal{F}}_A \ |\  A \in GL(2,{\mathbb{R}})\} = \{{\mathcal{F}}_A \ |\  A \in SL(2,{\mathbb{R}})\},$$
which is the orbit of the identity function $id = {\mathcal{F}}_{{\mathrm{Id}}}$ under the right action of $GL(2,{\mathbb{H}})$ on $\mathcal{L}_\infty$, does not inherit a group structure from $GL(2,{\mathbb{H}})$.

\begin{lemma}
The set ${\mathfrak{G}}$ of regular fractional transformations is also the orbit of $id$ with respect to the left action of $GL(2,{\mathbb{H}})$ on $\mathcal{L}_\infty$.
\end{lemma}

\begin{proof}
Let $A=\begin{bmatrix}
a & c \\
b & d
\end{bmatrix} \in GL(2,{\mathbb{H}})$, and let us prove that ${\mathcal{F}}_A = id.A$ can also be expressed as $C.id$ for some $C \in GL(2,{\mathbb{H}})$. If $c=0$ then
$${\mathcal{F}}_A(q) = d^{-1}*(qa+b) = (d^{-1}a)*q+ d^{-1}b,$$
else
$${\mathcal{F}}_A(q) = {\mathcal{F}}_{c^{-1}A}(q) = (q-p)^{-*}*(q\alpha +\beta) = [(q-p)^s]^{-1} (q-\bar p)*(q\alpha + \beta)$$
for some $p, \alpha, \beta \in {\mathbb{H}}$. If  $p= x+Iy$ then there exists $\tilde p \in x+y{\mathbb{S}}$ and $\gamma, \delta \in {\mathbb{H}}$ such that $ (q-\bar p)*(q\alpha + \beta) = (q \gamma + \delta)*(q-\tilde p)$; additionally, $(q-p)^s = (q-\tilde p)^s$ (see \cite{zeros} for details). Hence,
$${\mathcal{F}}_A(q) = [(q-\tilde p)^s]^{-1} (q \gamma + \delta)*(q-\tilde p) = (q \gamma + \delta) * (q-\overline{\tilde p})^{-*} =  (\gamma* q + \delta) * (q-\overline{\tilde p})^{-*},$$
which is of the desired form.
Similar manipulations prove that for all $C \in GL(2,{\mathbb{H}})$, the function $C.id$ equals ${\mathcal{F}}_A = id.A$ for some $A \in GL(2,{\mathbb{H}})$.
\end{proof}

We now state an immediate consequence of the previous lemma and of remark \ref{leftright}.

\begin{remark}
The set ${\mathfrak{G}}$ of regular fractional transformations is preserved by regular conjugation.
\end{remark}

\section{Regular M\"obius transformations of ${\mathbb{B}}$}\label{regularmoebius}

The regular fractional transformations that map the open quaternionic unit ball ${\mathbb{B}}$ onto itself, called \emph{regular M\"obius transformations  of ${\mathbb{B}}$}, are characterized by two results of \cite{moebius}.
\begin{theorem}
For all $A \in SL(2, {\mathbb{H}})$, the regular fractional transformation ${\mathcal{F}}_A$ maps ${\mathbb{B}}$ onto itself if and only if $A \in Sp(1,1)$, if and only if there exist (unique) $u \in \partial {\mathbb{B}}, a \in {\mathbb{B}}$ such that
\begin{equation}
{\mathcal{F}}_A(q) = (1-q \bar a)^{-*}*(q-a)u.
\end{equation}
In particular, the set $\mathfrak{M}= \{f \in {\mathfrak{G}} \ |\  f({\mathbb{B}}) = {\mathbb{B}}\}$ of the regular M\"obius transformations of ${\mathbb{B}}$ is the orbit of the identity function under the right action of $Sp(1,1)$.
\end{theorem}

\begin{theorem}
The class of regular bijective functions $f : {\mathbb{B}} \to {\mathbb{B}}$ coincides with the class $\mathfrak{M}$ of regular M\"obius transformations of ${\mathbb{B}}$.
\end{theorem}

As a consequence, the right action of $Sp(1,1)$ preserves the class of regular bijective functions from ${\mathbb{B}}$ onto itself.
We now study, more in general, the effect of the actions of $Sp(1,1)$ on the class
$$\frak{Reg}(\B,\B) = \{f:{\mathbb{B}} \to {\mathbb{B}}\ |\ f \mathrm{\ is\ regular}\}.$$

\begin{proposition}
If $f \in \frak{Reg}(\B,\B)$ then for all $a \in {\mathbb{B}}$
\begin{equation}
(1-f \bar a)^{-*}*(f-a) =(f-a)*(1-\bar a*f)^{-*}.
\end{equation}
Moreover, the left and right actions of $Sp(1,1)$ preserve $\frak{Reg}(\B,\B)$.
\end{proposition}

\begin{proof}
The fact that $(1-f\bar a)^{-*}*(f-a)=(f-a)*(1-\bar a*f)^{-*}$ is a consequence of proposition \ref{hermitian}.

Let us turn to the second statement, proving that for all $a \in {\mathbb{B}}, u,v \in \partial {\mathbb{B}}$, the function
$$\tilde f =v^{-1}*(1-f\bar a)^{-*}*(f-a) u = (v-f\bar av)^{-*}*(f-a) u$$
is in $\frak{Reg}(\B,\B)$. The fact that for all $a \in {\mathbb{B}}, u,v \in \partial {\mathbb{B}}$ the function $u*(f-a)*(1-\bar a*f)^{-*} *v^{-*}$ belongs to $\frak{Reg}(\B,\B)$ will then follow from the equality just proven.

The function $\tilde f$ is regular in ${\mathbb{B}}$ since $h=v- f \bar av$ has no zero in ${\mathbb{B}}$ (as a consequence of the fact that $|a|<1$, $|f|<1$ and $|v|=1$). Furthermore,
$$\tilde f = (v- (f \circ T_h) \bar av)^{-1}(f\circ T_h-a)u=$$
$$v^{-1}(1- (f \circ T_h) \bar a)^{-1}(f\circ T_h-a)u =v^{-1}(M_a \circ f \circ T_h)u$$
where $T_h$ and $M_a(q) = (1-q \bar a)^{-1}(q-a)$ map ${\mathbb{B}}$ to itself bijectively and $u,v \in \partial{\mathbb{B}}$.
Hence, $\tilde f=v^{-*}*(1-f\bar a)^{-*}*(f-a) *u \in \frak{Reg}(\B,\B)$, as desired.
\end{proof}

As a byproduct, we obtain that the orbit of the identity function under the left action of $Sp(1,1)$ equals $\mathfrak{M}$.

\begin{proposition}
If $f\in \frak{Reg}(\B,\B)$ then its regular conjugate $f^c$ belongs to $\frak{Reg}(\B,\B)$ as well. Furthermore, $f^c$ is bijective (hence an element of $\mathfrak{M}$) if and only if $f$ is.
\end{proposition}

\begin{proof}
Suppose $f^c(p) = a \in {\mathbb{H}} \setminus {\mathbb{B}}$ for some $p = x+Iy \in {\mathbb{B}}$. Then $p$ is a zero of the regular function $f^c-a$. By theorem \ref{conjugatezeros}, there exists $\tilde p \in x+y{\mathbb{S}} \subset {\mathbb{B}}$ such that $(f^c-a)^c = f-\bar a$ vanishes at $\tilde p$. Hence, $f({\mathbb{B}})$ includes $\bar a \in {\mathbb{H}} \setminus{\mathbb{B}}$, a contradiction with the hypothesis $f({\mathbb{B}}) \subseteq {\mathbb{B}}$.

As for the second statement, $f^c$ is bijective if and only if, for all $a \in {\mathbb{B}}$, there exists a unique $p \in {\mathbb{B}}$ such that $f^c(p) =a$. Reasoning as above, we conclude that this happens if and only if for all $a \in {\mathbb{B}}$, there exists a unique $\tilde p \in {\mathbb{B}}$ such that $f(\tilde p) = \bar a$. This is equivalent to the bijectivity of $f$.
\end{proof}

\section{Differential and metric properties of regular M\"obius transformations}\label{differential}

The present section is concerned with two natural questions:
\begin{enumerate}
\item[(a)] whether the regular M\"obius transformations are conformal (as the classical M\"obius transformations);
\item[(b)] whether they preserve the quaternionic Poincar\'e distance defined on ${\mathbb{B}}$ by formula (\ref{poincaredist}).
\end{enumerate}
For a complete description of the Poincar\'e metric, see \cite{poincare}. In order to answer question (a), we will compute for a regular M\"obius transformation the series development introduced by the following result of \cite{expansion}. Let us set
$$U(x_0+y_0{\mathbb{S}},r) = \{q \in {\mathbb{H}} \ |\  |(q-x_0)^2+y_0^2| < r^2\}$$
for all $x_0,y_0 \in {\mathbb{R}}$, $r>0$.

\begin{theorem}
Let $f$ be a regular function on $\Omega=B(0,R)$, and let  $U(x_0+y_0{\mathbb{S}},r) \subseteq \Omega$. Then for each $q_0 \in x_0+y_0 {\mathbb{S}}$ there exists $\{A_n\}_{n \in {\mathbb{N}}} \subset {\mathbb{H}}$ such that
\begin{equation}\label{expansion}
f(q) = \sum_{n \in {\mathbb{N}}} [(q-x_0)^2+y_0^2]^n [A_{2n}+(q-q_0) A_{2n+1}]
\end{equation}
for all $q \in U(x_0+y_0{\mathbb{S}},r)$. As a consequence,
\begin{equation}
\lim_{t\to 0} \frac{f(q_0+tv)-f(q_0)}{t} = v A_1 + (q_0v - v\bar q_0)A_2.
\end{equation}
for all $v \in T_{q_0}\Omega\cong {\mathbb{H}}$.
\end{theorem}

If $q_0 \in L_I$ and if we split the tangent space $T_{q_0}\Omega\cong {\mathbb{H}}$ as ${\mathbb{H}} = L_I \oplus L_I^\perp$, then the differential of $f$ at $q_0$ acts by right multiplication by $A_1$ on $L_I^\perp$ and by right multiplication by $A_1+2Im(q_0)A_2$ on $L_I$.
Furthermore, if for all $q_0 \in \Omega$ the \emph{remainder} $R_{q_0}f$ is defined as
$$R_{q_0}f(q) = (q-q_0)^{-*}*(f(q)-f(q_0))$$
then the coefficients of \eqref{expansion} are computed as $A_{2n} = (R_{\bar q_0}R_{q_0})^nf(q_0)$ and $A_{2n+1}= R_{q_0}(R_{\bar q_0}R_{q_0})^nf(\bar q_0)$.

Let us recall the definition of the \emph{Cullen derivative} $\partial_cf$, given in \cite{advances} as
\begin{equation}\label{cullen}
\partial_cf(x+Iy)=\frac{1}{2}\left(\frac{\partial}{\partial x}-I\frac{\partial}{\partial y}\right)f(x+Iy)
\end{equation}
for $I \in {\mathbb{S}},\ x,y \in {\mathbb{R}}$, as well as the definition of the \emph{spherical derivative}
\begin{equation}\label{spherical}
\partial_sf(q) = (2Im(q))^{-1}(f(q)-f(\bar q))
\end{equation}
given in \cite{perotti}. We can make the following observation.

\begin{remark}
If $f$ is a regular function on $B(0,R)$ and if \eqref{expansion} holds then $\partial_cf(q_0) = R_{q_0}f(q_0) = A_1+2Im(q_0)A_2$ and $\partial_sf(q_0)= R_{q_0}f(\bar q_0) = A_1$.
\end{remark}

In the case of the regular M\"obius transformation
$${\mathcal{M}}_{q_0}(q) = (1-q\bar q_0)^{-*}*(q-q_0) = (q-q_0)*(1-q\bar q_0)^{-*},$$
clearly $R_{q_0}{\mathcal{M}}_{q_0}(q) = (1-q\bar q_0)^{-*}$, so that we easily compute the coefficients $A_n$.

\begin{proposition}
Let $q_0 = x_0+y_0I \in {\mathbb{B}}$. Then the expansion \eqref{expansion} of ${\mathcal{M}}_{q_0}$ at $q_0$ has coefficients
\begin{eqnarray}
A_{2n-1} &=& \frac{\bar q_0^{2n-2}}{(1-|q_0|^2)^{n-1}(1-\bar q_0^2)^{n}}\\
A_{2n} &=& \frac{\bar q_0^{2n-1}}{(1-|q_0|^2)^n(1-\bar q_0^2)^n}
\end{eqnarray}
for all $n\geq 1$. As a consequence, for all $v \in {\mathbb{H}}$,
\begin{equation}
\frac{\partial{\mathcal{M}}_{q_0}} {\partial v}(q_0)=v (1-\bar q_0^2)^{-1}+(q_0 v-v \bar q_0) \frac{\bar q_0}{(1-|q_0|^2)(1-\bar q_0^2)}.
\end{equation}
\end{proposition}

\begin{proof}
We have already observed that $R_{q_0}{\mathcal{M}}_{q_0}(q) = (1-q\bar q_0)^{-*}$, so that $A_1 = R_{q_0}{\mathcal{M}}_{q_0}(\bar q_0) = (1-\bar q_0^2)^{-1}$. Moreover,
$$R_{\bar q_0}R_{q_0}{\mathcal{M}}_{q_0}(q) = (q-\bar q_0)^{-*}*\left[R_{q_0}{\mathcal{M}}_{q_0}(q)- A_{1}\right]=$$
$$=(q-\bar q_0)^{-*}*\left[(1-q\bar q_0)^{-*} - A_{1}\right] =$$
$$= (1-q\bar q_0)^{-*}* (q-\bar q_0)^{-*}*\left[(1-\bar q_0^2)- (1-q\bar q_0)\right]A_1=$$
$$= (1-q\bar q_0)^{-*} \bar q_0 A_1.$$
The thesis follows by induction, proving that for all $n \geq 1$
$$(R_{\bar q_0}R_{q_0})^{n}{\mathcal{M}}_{q_0}(q) =(1-q\bar q_0)^{-*} \bar q_0 A_{2n-1},$$
$$R_{q_0}(R_{\bar q_0}R_{q_0})^{n}{\mathcal{M}}_{q_0}(q) = (1-q\bar q_0)^{-*} \bar q_0 A_{2n}$$
by means of similar computations.
\end{proof}

We are now in a position to answer question (a).

\begin{remark}
For each $q_0 = x_0 + I y_0\in {\mathbb{B}} \setminus {\mathbb{R}}$, the differential of ${\mathcal{M}}_{q_0}$ at $q_0$ acts by right multiplication by $\partial_c{\mathcal{M}}_{q_0}(q_0)=(1-|q_0|^2)^{-1}$ on $L_I$ and by right multiplication by $\partial_s{\mathcal{M}}_{q_0}(q_0)=(1-\bar q_0^2)^{-1}$  on $L_I^\perp$. Since these quaternions have different moduli, ${\mathcal{M}}_{q_0}$ is not conformal.
\end{remark}

We now turn our attention to question (b): whether or not regular M\"obius transformations preserve the quaternionic Poincar\'e metric on ${\mathbb{B}}$ described in section \ref{classicM} and in \cite{poincare}. We recall that this metric was constructed to be preserved by the classical (non regular) M\"obius transformations of ${\mathbb{B}}$. Thanks to theorem \ref{quotients}, we observe what follows.

\begin{remark}
If ${\mathcal{M}}_{q_0}(q) = (1-q\bar q_0)^{-*}*(q-q_0)$ and $M_{q_0}(q) = (1-q\bar q_0)^{-1}(q-q_0)$ then
\begin{equation}
{\mathcal{M}}_{q_0}(q)= M_{q_0}(T(q)).
\end{equation}
where $T(q) = (1-q q_0)^{-1} q (1-q q_0)$ is a diffeomorphism of ${\mathbb{B}}$ with inverse $T^{-1}(q) = (1-q \bar q_0)^{-1} q (1- q \bar q_0)$.
\end{remark}

Thus, a generic regular M\"obius transformation of ${\mathbb{B}}$
$$q \mapsto {\mathcal{M}}_{q_0}(q) u = M_{q_0}(T(q)) u$$
(with $u \in \partial {\mathbb{B}}$) is an isometry if and only if $T$ is. An example shows that this is not the case.

\begin{example}
Let $q_0=\frac{I_0}{2}$ for some $I_0 \in {\mathbb{S}}$. Then $T(q) = (1-q q_0)^{-1} q (1-q q_0)$ is not an isometry for the Poincar\'e metric defined by formula \eqref{poincaredist}. Indeed, if $J_0 \in {\mathbb{S}}, J_0 \perp I_0$ then setting $q_1=\frac{J_0}{2}$ we have
$$\delta_{\mathbb{B}}(q_0,q_1) > \delta_{\mathbb{B}}(T(q_0),T(q_1))$$
since
$$\frac{|q_1-q_0|^2}{|1-q_1\bar q_0|^2} = \frac{|2J_0-2I_0|^2}{|4+J_0I_0|^2} =\frac{8}{17}$$
while, computing $T(q_0) = q_0 = \frac{I_0}{2}$ and $T(q_1) = \frac{8I_0+15J_0}{34}$, we conclude that
$$\frac{|T(q_1)-q_0|^2}{|1-T(q_1)\bar q_0|^2} = 4\frac{|-3I_0+5J_0|^2}{|20-5K_0|^2}= \frac{8}{25}.$$
Thus, the regular M\"obius transformations do not have a definite behavior with respect to $\delta_{\mathbb{B}}$: we have proven that $T$ (hence ${\mathcal{M}}_{q_0}$) is not an isometry, nor a dilation; the same computation proves that  $T^{-1}(q) = (1-q \bar q_0)^{-1} q (1- q \bar q_0)$ (hence ${\mathcal{M}}_{\bar q_0}$) is not a contraction.
\end{example}

The previous discussion proves that the study of regular M\"obius transformations cannot be framed into the classical study of ${\mathbb{B}}$, and that it requires further research. On the other hand, the theory of regular functions provides working tools that were not available for the classical M\"obius transformations. These tools lead us in \cite{schwarzpick} to an analog of the Schwarz-Pick lemma, which we will present in the next section.

\section{The Schwarz-Pick lemma for regular functions}\label{schwarzpick}

In the complex case, holomorphic functions play a crucial role in the study of the intrinsic geometry of the unit disc $\Delta = \{z \in {\mathbb{C}}\ |\  |z|<1\}$ thanks to the Schwarz-Pick lemma \cite{pick2,pick1}.

\begin{theorem}
Let $f:\Delta \to \Delta$ be a holomorphic function and let $z_0 \in \Delta$. Then
\begin{equation}
\left|\frac{f(z)-f(z_0)}{1-\overline{f(z_0)}f(z)}\right| \leq \left|\frac{z-z_0}{1-\bar z_0 z}\right|,
\end{equation}
for all $z \in \Delta$ and
\begin{equation}
\frac{|f'(z_0)|}{1-|f(z_0)|^2} \leq \frac{1}{1-|z_0|^2}.
\end{equation}
All inequalities are strict for $z \in \Delta \setminus \{z_0\}$, unless $f$ is a M\"obius transformation.
\end{theorem}

This is exactly the type of tool that is not available, in the quaternionic case, for the classical M\"obius transformations. To the contrary, an analog of the Schwarz-Pick lemma is proven in \cite{schwarzpick} for quaternionic regular functions.
To present it, we begin with a result concerning the special case of a function $f \in \frak{Reg}(\B,\B)$ having a zero.

\begin{theorem} \label{teo1}
If $f : {\mathbb{B}} \to {\mathbb{B}}$ is regular and if $f(q_0)=0$ for some $q_0 \in {\mathbb{B}}$, then
\begin{equation}
|{\mathcal{M}}_{q_0}^{-*}*f(q)| \leq 1
\end{equation}
for all $q \in {\mathbb{B}}$. The inequality is strict, unless ${\mathcal{M}}_{q_0}^{-*}*f(q) \equiv u$ for some $u \in \partial {\mathbb{B}}$.
\end{theorem}

A useful property of the moduli of regular products is also proven in \cite{schwarzpick}:

\begin{lemma}\label{modulusproduct}
Let $f,g,h: B(0,R) \to {\mathbb{H}}$ be regular functions. If $|f|\leq |g|$ then $|h*f| \leq |h*g|$. Moreover, if $|f|< |g|$ then $|h*f| < |h*g|$ in $B \setminus {\mathcal{Z}}_h$.
\end{lemma}

The property above allows us to derive from theorem \ref{teo1} the perfect analog of the Schwarz-Pick lemma in the special case $f(q_0)=0$. We recall that $\partial_c f$ denotes the Cullen derivative of $f$, defined by formula \eqref{cullen}, while $\partial_s f$ denotes the spherical derivative, defined by formula \eqref{spherical}.

\begin{corollary}\label{schwarzp}
If $f : {\mathbb{B}} \to {\mathbb{B}}$ is regular and if $f(q_0)=0$ for some $q_0 \in {\mathbb{B}}$ then
\begin{equation}
|f(q)| \leq |{\mathcal{M}}_{q_0}(q)|
\end{equation}
for all $q \in {\mathbb{B}}$. The inequality is strict at all $q \in {\mathbb{B}} \setminus \{q_0\}$, unless there exists $u \in \partial {\mathbb{B}}$ such that $f(q) = {\mathcal{M}}_{q_0}(q) \cdot u$ at all $q \in {\mathbb{B}}$.
Moreover, $|R_{q_0}f(q)| \leq |(1-q\bar q_0)^{-*}|$ in ${\mathbb{B}}$ and in particular
\begin{eqnarray}
|\partial_c f (q_0)| \le \frac{1}{1-|q_0|^2}\\
|\partial_s f(q_0)|
\leq \frac{1}{|1-\overline{q_0}^2|}.
\end{eqnarray}
These inequalities are strict, unless $f(q) = {\mathcal{M}}_{q_0}(q) \cdot u$ for some $u \in \partial {\mathbb{B}}$.
\end{corollary}

Finally, the desired result is obtained in full generality.

\begin{theorem}[Schwarz-Pick lemma] \label{MainSchwarzPick}
Let $f : {\mathbb{B}} \to {\mathbb{B}}$ be a regular function and let $q_0 \in {\mathbb{B}}$. Then
\begin{equation}\label{MainSchwarzPickEq}
|(f(q)-f(q_0))*(1-\overline{f(q_0)}*f(q))^{-*}| \leq |(q-q_0)*(1-\bar q_0*q)^{-*}|,
\end{equation}
\begin{equation}
|R_{q_0}f(q)*(1-\overline{f(q_0)}*f(q))^{-*}| \leq |(1-\bar q_0*q)^{-*}|
\end{equation}
in ${\mathbb{B}}$. In particular,
\begin{equation} \label{culderiv}
|\partial_c f *(1-\overline{f(q_0)}*f(q))^{-*}|_{|_{q_0}} \leq \frac{1}{1-|q_0|^2}.
\end{equation}
Apart from \eqref{MainSchwarzPickEq} at $q_0$ (which reduces to $0\leq0$), all inequalities are strict unless $f$ is a regular M\"obius transformation.
\end{theorem}

This promising result makes it reasonable to expect that regular functions play an important role in the intrinsic geometry of the quaternionic unit ball. Therefore, it encourages to continue the study of regular M\"obius transformations and of their differential or metric properties.