\section{Introduction}

Two companion technical papers
\cite{Hamilton:2010b,Hamilton:2010c},
hereafter Papers~2 and 3,
present
conformally stationary, axisymmetric, conformally separable
solutions for
the interior structure of a rotating black hole.
Paper~2 deals with uncharged black holes,
while Paper~3 generalizes to charged black holes.
The claim is that these papers provide a solution to a problem that has
remained outstanding since
Kerr's 1963 \cite{Kerr:1963}
discovery of the exterior geometry of a rotating black hole.
The purpose of the present paper is to
give an abbreviated derivation of the solution
for an uncharged black hole,
and to summarize the principal features of the solution.
A Mathematica notebook containing many details of the calculations
is at \cite{Hamilton:notebook}.

The papers consider only classical general relativity,
not alternate theories of gravity,
nor speculative quantum processes
that might for example destabilize the outer horizon.

\penrosekerrinflationfig

The Penrose diagram of the analytically extended Kerr geometry,
Figure~\ref{penrosekerrinflation},
provides a good starting point
for understanding where and how the interior Kerr geometry fails.
A spherical charged (Reissner-Nordstr\"om) black hole
has a similar interior structure,
with essentially the same Penrose diagram,
and much of the literature has focussed on this simpler case.
The Kerr geometry, and more generally the Kerr-Newman geometry,
has two inner horizons that are gateways to regions of unpredictability,
signalled by the presence of timelike singularities.
In 1968, Penrose
\cite{Penrose:1968}
pointed out that an observer passing through the outgoing inner horizon
(the Cauchy horizon)
of a spherical charged black hole
would see the outside Universe infinitely blueshifted,
and he suggested that the infinite blueshift would destabilize the inner horizon.
The infinite blueshift is plain from the Penrose diagram,
Figure~\ref{penrosekerrinflation},
which shows that a person passing through the outgoing inner horizon
sees the entire future of the outside Universe go by in a finite time.
Perturbation theory,
much of it expounded in Chandrasekhar's (1983) monograph
\cite{Chandrasekhar:1983},
confirmed that waves from the outside Universe would
amplify to a diverging energy density on the outgoing inner horizon
of a spherical charged black hole.
The result was widely interpreted as indicating the instability
of the inner horizon.

It was not until 1990 that
the full nonlinear nature of the instability at the inner horizon
was eventually clarified by Poisson \& Israel
\cite{Poisson:1990eh}.
Poisson \& Israel
showed that if ingoing
and outgoing
streams are simultaneously present just above the inner horizon
of a spherical charged black hole,
then cross-flow between the two streams would lead to an
exponential growth of the interior mass.
They called the instability ``mass inflation.''
The inflationary instability in spherical charged black holes
was confirmed analytically and numerically
in several studies,
as reviewed by \cite{Hamilton:2008zz}.

The physical reason for the inflationary instability
can be seen in the Penrose diagram, Figure~\ref{penrosekerrinflation}.
In the Black Hole region between the outer and inner horizons,
the time coordinate (the one that expresses time translation symmetry)
is spacelike,
so that it is possible to go either forward or backward in time.
Inside the inner horizon, the time coordinate reverts to being timelike.
Ingoing
particles want to fall into a region where
the time coordinate is progressing forwards,
while outgoing
particles want to fail into a region where
the time coordinate is progressing backwards.
The Penrose diagram, Figure~\ref{penrosekerrinflation},
shows that indeed there are two distinct ingoing and outgoing inner horizons
which disgorge on to two causally separated pieces of spacetime
where the time coordinate is pointed in opposite directions.
To achieve this causal separation,
the ingoing and outgoing streams must exceed the speed of light
relative to each other.
This is Penrose's infinite blueshift.
In reality, if ingoing and outgoing streams are present,
then their attempt to exceed the speed of light relative to each other
produces a counter-streaming energy and pressure
that, however tiny the initial streams may be,
inevitably grows to the point that it becomes a significant source of gravity.
As expounded by \cite{Hamilton:2008zz},
the counter-streaming pressure produces a gravitational force
that is in opposite directions for ingoing and outgoing streams,
accelerating the streams ever faster through each other,
in turn increasing the counter-streaming pressure.
The inflationary instability
thus grows exponentially.

Almost the only prior work on the inflationary instability inside
rotating black holes is that of Ori
\cite{Ori:1992,Ori:2001pc}
(see \cite{Hamilton:2009hu} for a full list of references).
Ori considered the situation in which
a black hole undergoes collapse,
and inflation is driven by Price tails of ingoing and outgoing
gravitational radiation generated by the collapse.
The situation is more complicated than that considered
in the present series of papers,
and Ori was able to advance only qualitative arguments about the
character of the instability.

\section{Approach}

The strategy adopted in the present papers
is motivated by two key physical insights.

The first insight is that,
as shown in \S{V} of Paper~2 \cite{Hamilton:2010b},
collisionless ingoing and outgoing streams
falling towards the inner horizon of the Kerr-Newman geometry
become highly focussed into twin narrow, intense beams
pointed along the ingoing and outgoing principal null directions.
The focussing is along these two special directions
regardless of the initial distributions of orbital parameters of the streams.
This implies that the energy-momentum tensor of the ingoing and outgoing streams
takes a simple and predictable form near the inner horizon.
(We use the term inner horizon to describe
the narrow region where inflation takes place,
even though the inner horizon is destroyed by inflation,
and therefore does not actually exist.)

As first shown by
\cite{Carter:1968c},
the Kerr-Newman geometry
(and some other electrovac geometries)
is Hamilton-Jacobi separable in a tetrad aligned with the
principal null directions.
The fact that collisionless streams focus near the inner horizon
along precisely these principal null directions
suggests that the spacetime might continue to be separable
in the presence of inflation.

\kerrrcontourfig

The second insight is that
the geometry of a spherical charged black hole undergoing inflation
at (just above) its inner horizon
has a step-function character.
The spacetime is well-approximated by the electrovac
(Reissner-Nordstr\"om) geometry down to just above the inner horizon.
Then, in a tiny interval of radius and proper time near the inner horizon,
the center-of-mass counter-streaming energy and pressure,
the Weyl curvature, and the interior mass all inflate
to exponentially huge values.
Counter-intuitively,
the smaller the incident ingoing and outgoing streams,
the more rapidly quantities exponentiate
\cite{Hamilton:2008zz}.
In the limit of tiny accretion rate,
the geometry tends to a step-function.
This suggests that inflationary spacetimes might be found
by looking for solutions with a step-like character.
The step-like character of the inflationary solutions
can be seen in the sharp turns at the inner horizon
in the contours of constant radius and latitude in Figure~\ref{kerrrcontour}.

As an aside, it is worth commenting on the challenges and pitfalls
of computing inflationary spacetimes numerically
rather than analytically.
One major numerical challenge arises from the fact that
during inflation
physical quantities inflate to exponentially huge values
over tiny intervals of distance and time.
One potential pitfall is that inflation requires
ingoing and outgoing streams that can stream relativistically
through each other.
A code, or indeed analytic model,
that treats the matter as a single fluid with a sound speed
less than the speed of light
artificially suppresses inflation
by disallowing the relativistic counter-streaming that drives it.

\section{Summary derivation}
\label{derivation}

This section summarizes the derivation of the solution
for an uncharged rotating black hole.
Complete details are given in Paper~2.

The solutions presented in this series of papers are
conformally stationary, axisymmetric, and conformally separable.
Conformally stationarity combines the assumption of
conformal time translation invariance (self-similarity)
with an infinitesimal expansion rate.
Let $x^\mu \equiv \{ x , t , y , \phi \}$ be coordinates
in which $t$ is conformal time,
$\phi$
is the azimuthal angle,
and $x$ and $y$ are radial and angular coordinates.
As shown in Appendix~A of Paper~2 \cite{Hamilton:2010b},
the line-element may be taken to be
\begin{equation}
\label{lineelement}
  d s^2
  =
  \rho^2
  \left[
  {d x^2 \over \Delta_x}
  -
  {\Delta_x \over \sigma^4}
  \left( d t - \varpi \, d \phi \right)^2
  +
  {d y^2 \over \Delta_y}
  +
  {\Delta_y \over \sigma^4}
  \left( d \phi - \omega \, d t \right)^2
  \right]
  \ ,
\end{equation}
where
\begin{equation}
  \sigma
  \equiv
  \sqrt{ 1 - \omega \, \varpi }
  \ .
\end{equation}
The determinant of the $2 \times 2$ submatrix
of $t$--$\phi$ metric coefficents
defines the radial and angular horizon functions
$\Delta_x$ and $\Delta_y$:
\begin{equation}
  g_{tt} g_{\phi\phi} - g_{t\phi}^2
  =
  -
  {\rho^4 \over \sigma^4}
  \Delta_x \Delta_y
  \ .
\end{equation}
Horizons occur when one or other of the horizon functions
$\Delta_x$ and $\Delta_y$
vanish.
The focus here is the region near the inner horizon
where the radial horizon function $\Delta_x$ is negative and tending to zero.
The line-element~(\ref{lineelement}) defines a tetrad
that is aligned with the principal null directions.
Since the radial coordinate $x$ is timelike
near the inner horizon,
it is convenient to take $x$ as the time coordinate of the tetrad,
and to choose the sign of $x$ so that it increases inwards,
the direction of advancing time.

The conformal factor $\rho$
is a product of separable (electrovac)
$\rho_{\rm s}$,
time-dependent
$e^{\text{\sl v} t}$,
and
inflationary
$e^{-\xi}$
factors,
\begin{equation}
\label{rho}
  \rho
  =
  \rho_{\rm s}
  e^{\text{\sl v} t - \xi}
  \ .
\end{equation}
Conformal time translation symmetry is expressed by the fact
that the spacetime expands conformally
(that is, without changing shape)
by factor $\rho \rightarrow e^{\text{\sl v} \Delta t} \rho$
when the conformal time increases by $t \rightarrow t + \Delta t$.
Conformal stationarity means taking the limit
of small expansion rate, or small accretion rate,
{\em after\/} calculations are complete,
\begin{equation}
  \text{\sl v} \rightarrow 0
  \ .
\end{equation}
This is not the same as stationarity,
which sets $\text{\sl v}$ to zero at the outset.
A feature of inflation is that the smaller the accretion rate,
the faster inflation exponentiates.
Even in the limit of infinitesimal accretion rate,
inflation drives the
center-of-mass streaming density and pressure,
and the Weyl curvature, to exponentially huge values.
Mathematically,
Einstein's
equations contain terms of order $\sim \text{\sl v} / \Delta_x$
that grow large at the inner horizon $\Delta_x \rightarrow -0$
however small the accretion rate $\text{\sl v}$ may be.

This paper adopts a collisionless fluid as the source of
energy-momentum that ignites and then drives inflation.
Because collisionless particles stream hyper-relativistically
through each other during inflation,
the trajectories of massive freely-falling particles are
well-approximated by those of massless particles.
Conformal separability posits that the equations of motion
of freely-falling massless particles are Hamilton-Jacobi separable,
which implies that a conformal Killing tensor exists.
As shown in Appendix~A of Paper~2,
conformal separability requires that
\begin{equation}
\label{fnxy}
  \begin{array}{ccl}
  \omega
  \ ,
  &
  \Delta_x

  &
  \mbox{~are functions of $x$ only}
  \ ,
  \\
  \varpi
  \ ,
  &
  \Delta_y

  &
  \mbox{~are functions of $y$ only}
  \ .
  \end{array}
\end{equation}
Unlike strict separability (for massive as well massless particles),
conformal separability does not impose
conditions on the conformal factor $\rho$.

Conformally separable inflationary solutions
are obtained by separating the Einstein equations systematically.
Homogeneous solution of the Einstein components
$G_{xy}$, $G_{t\phi}$, and $G_{xx} + G_{tt}$ and $G_{yy} - G_{\phi\phi}$
leads to the usual electrovac solutions
for the electrovac conformal factor $\rho_{\rm s}$
and the vierbein coefficients $\omega$ and $\varpi$:
\begin{equation}
\label{rhoxy}
  \rho_{\rm s}
  =
  \sqrt{
  \rho_x^2
  +
  \rho_y^2
  }
  \ , \quad
  \rho_x
  =
  \sqrt{
  {g_0 - g_1 \omega
  \over
  ( f_0 g_1 + f_1 g_0 )
  ( f_0 + f_1 \omega ) }
  }
  \ , \quad
  \rho_y
  =
  \sqrt{
  {g_1 - g_0 \varpi
  \over
  ( f_0 g_1 + f_1 g_0 )
  ( f_1 + f_0 \varpi )}
  }
  \ ,
\end{equation}
\begin{equation}
\label{domegadvarpi}
  {d \omega \over d x}
  =
  2
  \sqrt{
  \left( f_0 + f_1 \omega \right)
  \left( g_0 - g_1 \omega \right)
  }

  \ , \quad
  {d \varpi \over d y}
  =
  2
  \sqrt{
  \left( f_1 + f_0 \varpi \right)
  \left( g_1 - g_0 \varpi \right)
  }

  \ ,
\end{equation}
where $f_0$, $f_1$, $g_0$, and $g_1$ are constants determined
by boundary conditions.
Equations~(\ref{rhoxy}) and (\ref{domegadvarpi})
continue to hold throughout inflation and collapse.
The inflationary solutions are generic,
applying wherever a separable electrovac spacetime
has an inner horizon,
so the specific choice of constants
$f_0$, $f_1$, $g_0$, and $g_1$
does not affect the argument.

The most important Einstein equations,
since they lead to equations governing the evolution of
the inflationary exponent $\xi$
and the horizon function $\Delta_x$,
are those for the Einstein components
$G_{xx} - G_{tt}$ and $G_{yy} + G_{\phi\phi}$.
The collisionless source of both these components can be treated
as negligible, in the conformally stationary limit of small accretion rate.
The angular components are negligible because inflation amplifies
the radial, not angular components;
and the trace of the collisionless energy-momentum remains small
because it depends on the rest mass of the particles,
which is unchanged by inflation.
Define $U_x$, $U_y$, $X_x$, $X_y$, $Y_x$, and $Y_y$ by
\begin{subequations}
\begin{align}
\label{Uxy}
  U_x
  \equiv
  -
  {\partial \xi \over \partial x} \Delta_x
  \ &, \quad
  U_y
  \equiv
  {\partial \xi \over \partial y} \Delta_y
  \ ,
\\
\label{Xxy}
  X_x
  \equiv
  {\partial U_x \over \partial x}
  +
  2 {U_x^2 - \text{\sl v}^2 \over \Delta_x}
  \ &, \quad
  X_y
  \equiv
  {\partial U_y \over \partial y}
  -
  2 {U_y^2 + \text{\sl v}^2 \varpi^2 \over \Delta_y}
  \ ,
\\
\label{Yxy}
  Y_x
  \equiv
  {d \Delta_x \over d x}
  +
  3 U_x
  -
  \Delta_x
  {d \over d x}
  \ln
  \left[ ( f_0 {+} f_1 \omega ) {d \omega \over d x} \right]
  \ &, \quad
  Y_y
  \equiv
  {d \Delta_y \over d y}
  -
  3 U_y
  -
  \Delta_y
  {d \over d y}
  \ln
  \left[ ( f_1 {+} f_0 \varpi ) {d \varpi \over d y} \right]
  \ .
\end{align}
\end{subequations}
In terms of these quantities, the Einstein components
$G_{xx} - G_{tt}$ and $G_{yy} + G_{\phi\phi}$ are
\begin{subequations}
\label{GUv}
\begin{align}
\label{GxxmttUv}
  \rho^2 \left( G_{xx} - G_{tt} \right)
  &=
  {1 \over \sigma^2}
  \left(
  Y_x {d \ln \omega \over d x}
  -
  Y_y {d \ln \varpi \over d y}
  \right)
  -
  2 X_x
  +
  Y_x
  {d
  \over d x}
  \ln
  \left(
  {f_0 {+} f_1 \omega \over \omega}
  \right)
  +
  X_y
  -
  {\partial Y_y \over \partial y}
  +
  Y_y
  {d
  \over d y}
  \ln
  \left[
  {\varpi ( f_1 {+} f_0 \varpi ) \over d \varpi / d y}
  \right]
\nonumber
  \\
  & \quad
  + \,
  U_x
  {\partial
  \over \partial x}
  \ln
  \left[ \sigma^2 ( f_0 {+} f_1 \omega ) \right]
  -
  U_y
  {\partial \over \partial y}
  \ln \left[
  {( g_1 {-} g_0 \varpi ) \over \sigma^2} {d \varpi \over d y}
  \right]
  \ ,
\\
\label{GyypphiphiUv}
  \rho^2 \left( G_{yy} + G_{\phi\phi} \right)
  &=
  {1 \over \sigma^2}
  \left(
  Y_x {d \ln \omega \over d x}
  -
  Y_y {d \ln \varpi \over d y}
  \right)
  -
  2 X_y
  -
  Y_y
  {d \over d y}
  \ln \left( {f_1 {+} f_0 \varpi \over \varpi} \right)
  +
  X_x
  +
  {\partial Y_x \over \partial x}
  -
  Y_x
  {d \over d x}
  \ln
  \left[
  {\omega ( f_0 {+} f_1 \omega ) \over d \omega / d x}
  \right]
\nonumber
  \\
  & \quad
  + \,
  U_y
  {\partial \over \partial y}
  \ln \left[ \sigma^2 ( f_1 {+} f_0 \varpi ) \right]
  -
  U_x
  {\partial \over \partial x}
  \ln \left[
  {( g_0 {-} g_1 \omega ) \over \sigma^2} {d \omega \over d x}
  \right]
  \ .
\end{align}
\end{subequations}
Homogeneous solutions of these equations can be found by supposing that
$U_x$, $X_x$, and $Y_x$ are all functions of radius $x$,
while
$U_y$, $X_y$, and $Y_y$ are all functions of radius $y$,
and by separating each of the equations as
\begin{equation}
\label{sepG}
  {1 \over \sigma^2}
  \left(
  {f_0 h_0 {+} h_2 \omega {+} f_1 h_1 \omega^2 \over \omega}
  -
  {f_1 h_1 {+} h_2 \varpi {+} f_0 h_0 \varpi^2 \over \varpi}
  \right)
  -
  {f_0 h_0 {+} h_3 \omega \over \omega}
  +
  {f_1 h_1 {+} h_3 \varpi \over \varpi}
  =
  0
  \ ,
\end{equation}
for some constants $h_0$, $h_1$, $h_2$, and $h_3$.
If one attempts to separate equations~(\ref{GUv}) exactly,
then the attempt fails unless $U_x$ and $U_y$ are identically zero,
which is the usual electrovac case.
But if $U_x$ is taken to be small but finite,
then separation succeeds,
and inflation emerges.
If $U_x$ and $U_y$ on the second lines of equations~(\ref{GUv})
are treated as negligibly small,
then separating the first lines of each of equations~(\ref{GUv})
according to the pattern of equation~(\ref{sepG})
leads to the homogeneous solutions
\begin{subequations}
\label{XYxhomog}
\begin{align}
\label{Xxhomog}
  X_x
  =
  0
  \ &, \quad
  X_y
  =
  0
  \ ,
\\
\label{Yxhomog}
  Y_x
  =
  {( f_0 + f_1 \omega ) ( h_0 + h_1 \omega ) \over d \omega / d x}
  \ &, \quad
  Y_y
  =
  {( f_1 + f_0 \varpi ) ( h_1 + h_0 \varpi ) \over d \varpi / d y}
  \ .
\end{align}
\end{subequations}
If $U_x = U_y = 0$,
then solution of equations~(\ref{Yxy}) and (\ref{Yxhomog})
for $Y_x$ and $Y_y$,
subject to appropriate boundary conditions,
yields the radial and angular horizon functions $\Delta_x$ and $\Delta_y$
of the Kerr line-element.
The result is easily generalized to other electrovac spacetimes
by admitting appropriate sources for $Y_x$ and $Y_y$.

The quantity $X_y$ defined by equation~(\ref{Xxy})
determines the evolution of $U_y$,
and the solution $X_y = 0$, equation~(\ref{Xxhomog}),
then implies that inflation leaves $U_y$ unchanged,
and sensibly equal to its electrovac value of zero,
in the conformally stationary limit.
Thus inflation leaves the angular horizon function $\Delta_y$
unchanged from its electrovac value.

On the other hand, inflation drives $U_x$ away from zero
however small it might initially be.
In the vicinity of the inner horizon, where $\Delta_x \rightarrow -0$,
the solutions~(\ref{XYxhomog}) for $X_x$ and $Y_x$
defined by equations~(\ref{Xxy}) and (\ref{Yxy})
imply the evolution equations
\begin{subequations}
\label{UDx}
\begin{align}
\label{Ux}
  {\partial U_x \over \partial x}
  +
  2 {U_x^2 - \text{\sl v}^2 \over \Delta_x}
  &=
  0
  \ ,
\\
\label{Dx}
  {d \Delta_x \over d x}
  +
  3 U_x
  &=
  \Delta_x^\prime
  \ ,
\end{align}
\end{subequations}
where
$\Delta_x^\prime \equiv \left. d \Delta_x / d x \right|_{x_-}$
is the (positive) derivative of the electrovac
horizon function at the inner horizon
$x = x_-$.
Below, equation~(\ref{xinf}),
it will be found that the radius $x$ remains frozen
at its inner horizon value $x_-$ throughout inflation and collapse,
so the right hand side of equation~(\ref{Dx}),
which is the electrovac solution for $Y_x$
evaluated at the inner horizon,
is constant during inflation and collapse.
The evolution equation~(\ref{Ux}) for $U_x$
involves a term inversely proportional to the horizon function $\Delta_x$,
which diverges at the inner horizon $\Delta_x \rightarrow -0$,
driving $U_x$ away from zero
however small $U_x$ might initially be.

Equation~(\ref{Ux}) leads to instability only at the inner horizon,
where $\Delta_x \rightarrow -0$.
At the outer horizon, where $\Delta_x \rightarrow 0$,
solutions of equation~(\ref{Ux}) decay rather than grow.

The separation of the Einstein components~(\ref{GUv})
that leads to the evolution equations~(\ref{UDx})
was premised on the assumption that the terms proportional
to $U_x$ and $U_y$ on the second lines of equations~(\ref{GUv})
could be neglected.
However,
the separation continues to remain valid during inflation
and collapse when $U_x$ grows huge.
The reason for this is that the dominant terms
in the Einstein components~(\ref{GUv})
during inflation and collapse
are of order $U_x^2 / \Delta_x$,
coming from the expression~(\ref{Xxy}) for $X_x$.
Thus,
once $U_x$ ceases to be negligble,
the condition for the validity of the separation becomes
$U_x \ll U_x^2 / | \Delta_x |$,
or equivalently
$| \Delta_x | \ll U_x$.
Consequently the condition for the validity of the separation
of the Einstein components~(\ref{GUv}) is
\begin{equation}
\label{sepcondition}
  \mbox{either}
  \quad
  U_x \ll 1
  \quad
  \mbox{or}
  \quad
  | \Delta_x | \ll U_x
  \ .
\end{equation}
Condition~(\ref{sepcondition})
holds from electrovac through inflation and collapse,
provided that the accretion rate is small,
as conformal stationarity prescribes.
The fact that condition~(\ref{sepcondition}) suffices
is verified in \S{VIII\,I} of Paper~2,
where the Einstein equations are solved to next order in $\Delta_x / U_x$,
and it is shown that the effect on the evolution of the inflationary
exponent $\xi$ and horizon function $\Delta_x$ is negligible.

The evolution equations~(\ref{UDx}) are solved
in the next section, \S\ref{inflationcollapse},
but first it is necessary to attend to the other Einstein equations.

In the conformally separable geometry,
freely-falling
collisionless ingoing and outgoing streams
become highly focussed along
the principle ingoing and outgoing null directions
as they
approach the inner horizon.
Inflation accelerates the streams even faster along the same null directions,
causing the $x$ and $t$ components of the tetrad-frame momenta
of freely-falling collisionless streams to grow exponentially.
Consequently the collisionless energy-momentum is dominated by
its $x$-$t$ components.
The associated components of the Einstein tensor are
\begin{equation}
\label{GxtxUv}
  \rho^2
  \left(
  {G_{xx} + G_{tt} \over 2}
  \,\pm\,
  G_{xt}
  \right)
  =
  ( U_x \mp \text{\sl v} )
  \left[
  {Y_x \pm \text{\sl v} \over - \Delta_x}
  -
  {d \over d x}
  \ln \left( d \omega \over d x \right)
  \right]
  +
  X_x
  \ .
\end{equation}
Since $X_x = 0$,
and $Y_x = \Delta_x^\prime$,
and the term proportional to
${d \ln ( d \omega / d x ) / d x}$
is sub-dominant
(in fact the term disappears when the Einstein components~(\ref{GUv})
and corresponding evolution equations~(\ref{UDx})
are solved to next order in $\Delta_x / U_x$; see Paper~2),
equation~(\ref{GxtxUv}) simplifies to
\begin{equation}
\label{GxtxUvsimp}
  \rho^2
  \left(
  {G_{xx} + G_{tt} \over 2}
  \,\pm\,
  G_{xt}
  \right)
  =
  {1 \over - \Delta_x}
  ( U_x \mp \text{\sl v} )
  ( \Delta_x^\prime \pm \text{\sl v} )
  \ .
\end{equation}
The right hand side of equation~(\ref{GxtxUv}) agrees with $8\pi$
times the energy-momentum tensor of two collisionless streams,
one ingoing ($+$) and one outgoing ($-$),
\begin{equation}
  T_{kl}
  =
  N^+ p^+_k p^+_l
  +
  N^- p^-_k p^-_l
  \ ,
\end{equation}
with densities
\begin{equation}
\label{Npm}
  N^\pm
  =
  {1 \over 16\pi}
  ( U_x \mp \text{\sl v} )
  ( \Delta_x^\prime \pm \text{\sl v} )
  \ ,
\end{equation}
and tetrad-frame momenta
\begin{equation}
\label{ppm0}
  p^\pm_k
  =
  {1 \over \rho}
  \left\{
  -
  {1 \over \sqrt{- \Delta_x}}
  \, , \
  \mp
  {1 \over \sqrt{- \Delta_x}}
  \, , \
  0
  \, , \
  0
  \right\}
  \ .
\end{equation}
The tetrad-frame momenta~(\ref{ppm0}) are null vectors
pointed along the principal ingoing and outgoing null directions.
That equations~(\ref{Npm}) and (\ref{ppm0}) describe correctly
the behaviour of freely-falling streams can be shown by
solving the Hamilton-Jacobi and collisionless Boltzmann equations
(see Paper~2),
and can be confirmed by checking that the densities and momenta
satisfy, to requisite accuracy, covariant number conservation,
$D^k N^\pm p^\pm_k = 0$,
and the geodesic equation
$d p^\pm_k / d \lambda = 0$,
where $\lambda$ is an affine parameter.
It might seem somewhat miraculous that the
$x$-$t$ components of the Einstein equations
are satisfied with a collisionless source,
but it is no coincidence.
Einstein's equations enforce covariant energy-momentum conservation,
$D^k T_{kl} = 0$.
Since the angular components are sub-dominant,
only the 3 distinct $x$-$t$ components of the energy-momentum tensor
are important.
The 3 components are subject to 2 energy-momentum conservation equations,
but in the present instance the 2 conservation equations are redundant,
so there is effectively 1 conservation equation.
However, the energy-momentum conservation equations for
ingoing and outgoing streams are symmetrically related to each other by
$\text{\sl v} \rightarrow - \text{\sl v}$,
so conservation of the sum of their energies,
as enforced by Einstein,
implies conservation of both.
The two conservation equations,
coupled with solution of the Einstein equation for
$G_{xx} - G_{tt}$, equation~(\ref{GxxmttUv}),
leads to a complete, self-consistent set of equations.

The angular motions of the freely-falling streams are small
compared to their radial motions, but not necessarily zero.
Next in order of magnitude,
after the 3 radial ($x$-$t$) components
of the energy-momentum tensor,
are its 4 off-diagonal radial-angular components.
The corresponding components of the Einstein tensor are
\begin{subequations}
\label{GxtyphiUv}
\begin{align}
\label{GxtyUv}
  \rho^2
  \left(
  G_{xy}
  \,\pm\,
  G_{ty}
  \right)
  &=
  -
  {1 \over \sqrt{- \Delta_x \Delta_y}}
  ( U_x \mp \text{\sl v} )
  \left(
  \Delta_y
  {\partial \ln \rho_{\rm s}^2 \over \partial y}
  -
  2 U_y
  \right)
  -
  \sqrt{- \Delta_x \over \Delta_y}
  \left(
  U_y
  {\partial \ln \rho_{\rm s}^2 \over \partial x}
  \, \pm \,
  {\text{\sl v} \varpi \over \sigma^2} {d \varpi \over d y}
  \right)

  \ ,
\\
\label{GxtphiUv}
  \rho^2
  \left(
  G_{x\phi}
  \,\pm\,
  G_{t\phi}
  \right)
  &=
  \pm
  {1 \over \sqrt{- \Delta_x \Delta_y}}
  ( U_x \mp \text{\sl v} )
  \left(
  {\Delta_y \over \sigma^2}
  {d \omega \over d x}
  \mp
  2 \text{\sl v} \varpi
  \right)
  \,\mp\,
  \sqrt{- \Delta_x \over \Delta_y}
  \left(
  {U_y \over \sigma^2} {d \varpi \over d y}
  \mp
  \text{\sl v} \varpi {\partial \ln \rho_{\rm s}^2 \over \partial x}
  \right)
  \ .
\end{align}
\end{subequations}
Since $U_y$
and the terms proportional to $\text{\sl v} \sqrt{- \Delta_x}$ are negligible,
equations~(\ref{GxtyphiUv}) simplify to
\begin{subequations}
\label{GxtyphiUvsimp}
\begin{align}
\label{GxtyUvsimp}
  \rho^2
  \left(
  G_{xy}
  \,\pm\,
  G_{ty}
  \right)
  &=
  -
  {1 \over \sqrt{- \Delta_x \Delta_y}}
  ( U_x \mp \text{\sl v} )
  \Delta_y
  {\partial \ln \rho_{\rm s}^2 \over \partial y}
  \ ,
\\
\label{GxtphiUvsimp}
  \rho^2
  \left(
  G_{x\phi}
  \,\pm\,
  G_{t\phi}
  \right)
  &=
  \pm
  {1 \over \sqrt{- \Delta_x \Delta_y}}
  ( U_x \mp \text{\sl v} )
  \left(
  {\Delta_y \omega^\prime \over \sigma^2}
  \mp
  2 \text{\sl v} \varpi
  \right)
  \ ,
\end{align}
\end{subequations}
where
$\omega^\prime \equiv \left. d \omega / d x \right|_{x_-}$,
which is effectively constant throughout inflation and collapse,
is the derivative of $\omega$ at the inner horizon $x = x_-$,
equation~(\ref{domegadvarpi}).
The right hand sides of equations~(\ref{GxtyphiUvsimp})
agree with $8\pi$ times the energy-momentum tensor of
ingoing and outgoing streams with the same densities $N^\pm$
as before, equation~(\ref{Npm}), but with tetrad-frame momenta
$p^\pm_k$
having finite rather than zero angular components:
\begin{equation}
\label{ppm}
  p^\pm_k
  =
  {1 \over \rho}
  \left\{
  -
  {1 \over \sqrt{- \Delta_x}}
  \, , \
  \mp
  {1 \over \sqrt{- \Delta_x}}
  \, , \
  {1 \over \sqrt{\Delta_y}}
  \left(
  {\Delta_y
  \partial \ln \rho_{\rm s}^2 / \partial y
  \over
  \Delta_x^\prime \pm \text{\sl v}}
  \right)
  \, , \
  \mp
  {1 \over \sqrt{\Delta_y}}
  \left(
  {\Delta_y \omega^\prime / \sigma^2 \mp 2 \text{\sl v} \varpi
  \over
  \Delta_x^\prime \pm \text{\sl v}}
  \right)
  \right\}
  \ .
\end{equation}
The tetrad-frame momenta~(\ref{ppm})
satisfy the Hamilton-Jacobi equations
with constant Hamilton-Jacobi parameters
along the path of the streams.
The angular components of the momenta
are small compared to the radial components as long as
\begin{equation}
\label{Dsmall}
  | \Delta_x | \ll 1
  \ .
\end{equation}
The momentum~(\ref{ppm}) is hyper-relativistic,
and $p^\pm_t = \pm p^\pm_x$ to an excellent approximation
so long as condition~(\ref{Dsmall}) is true.
If the condition~(\ref{Dsmall}) is violated,
then it signifies that angular motions are becoming important,
and the solution is breaking down.

It should be emphasized that,
as long as condition~(\ref{Dsmall}) holds,
the purely radial ($x$-$t$) Einstein equations hold regardless of
angular motions,
and thus the radial solution is unaffected by angular motions.
However,
if one requires that the sub-dominant radial-angular Einstein equations
are also satisfied,
then the angular motion of the collisionless streams must be
as given by equation~(\ref{ppm}).
One might perhaps have expected that conformally separable solutions
would require that the collisionless streams
would move exactly along the principal null directions,
but equation~(\ref{ppm}) shows that this is not true.
The fact that some motion is required in the angular $y$ direction (latitude)
can be attributed to the fact that surfaces of constant conformal factor $\rho$
do not coincide with surfaces of constant radius $x$.

Again, it might seem remarkable that the radial-angular Einstein equations
are satisfied by collisionless ingoing and outgoing streams.
And again, this coincidence results from energy-momentum conservation.
There are 4 radial-angular Einstein components,
subject to 2 energy-momentum conservation equations.
The energy-momentum conservation equations for the ingoing and outgoing
streams are symmetrically related by $\text{\sl v} \rightarrow - \text{\sl v}$,
so conservation of their sum implies conservation of both.

The final, sub-sub-dominant, components of the energy-momentum tensor
are the 3 purely angular ($y$-$\phi$) components.
The component $G_{yy} + G_{\phi\phi}$ component has already been addressed,
equation~(\ref{GyypphiphiUv}).
The remaining 2 components are
\begin{equation}
\label{GyphiyUv}
  \rho^2
  \left(
  {G_{yy} - G_{\phi\phi} \over 2}
  \,\pm\,
  i
  G_{y\phi}
  \right)
  =
  ( U_y \mp i \text{\sl v} \varpi )
  \left[
  {- \,
  Y_y
  \pm
  i \text{\sl v} \varpi \over \Delta_y}
  -
  {d \over d y}
  \ln \left( d \varpi \over d y \right)
  \right]
  +
  X_y
  \,\mp\,
  i \text{\sl v} {d \varpi \over d y}
  \ .
\end{equation}
Since $X_y = 0$,
and $U_y$ and $\text{\sl v}$ are negligibly small,
and there are no denominators of the radial horizon function $\Delta_x$,
equation~(\ref{GyphiyUv}) simplifies to
\begin{equation}
\label{GyphiyUvsimp}
  \rho^2
  \left(
  {G_{yy} - G_{\phi\phi} \over 2}
  \,\pm\,
  i
  G_{y\phi}
  \right)
  =
  0
  \ .
\end{equation}
During inflation,
the collisionless streams have negligble angular components
of energy-momentum
because the densities of the accreting streams are negligible,
in the conformally stationary limit,
and inflation does not amplify angular motions.
During collapse, the conformal factor $\rho$ shrinks,
and angular motions grow.
However, as long as the angular motions are sub-dominant,
which is true as long as condition~(\ref{Dsmall}) is satisfied,
the angular components of the energy-momentum can be neglected consistently:
the Einstein equations for the purely radial components,
and for the radial-angular components,
are unaffected by the angular components~(\ref{GyphiyUvsimp})
(the angular component $G_{yy} + G_{\phi\phi}$,
along with $G_{xx} - G_{tt}$,
equations~(\ref{GUv}),
determined the angular horizon function $\Delta_y$,
which inflation leaves unaltered from its electrovac form).
Equation~(\ref{GyphiyUvsimp}) requires that the $2 \times 2$
angular submatrix of the energy-momentum tensor be isotropic,
proportional to the $2 \times 2$ unit matrix.
As discussed in Paper~2,
it is possible to arrange the angular energy-momentum to be isotropic
by admitting multiple ingoing and outgoing streams,
with mean momenta set by equation~(\ref{ppm})
and isotropic mean squared momenta.
Treating the diagonal components of the angular energy-momentum
requires taking equations~(\ref{GUv}) to next order in $\Delta_x / U_x$,
but this can be done.

Eventually however,
the angular components do become important,
when $| \Delta_x | \sim 1$,
and the solution fails.

\inflationaryhorizonfig

\section{Inflation and collapse}
\label{inflationcollapse}

Denote the initial value of $U_x$, equation~(\ref{Uxy}),
incident on the inner horizon by $\text{\sl u}$,
a small parameter of order $\text{\sl v}$,
\begin{equation}
  U_x = \text{\sl u}
  \quad
  \mbox{initially}
  \ .
\end{equation}
The densities $N^\pm$ of ingoing and outgoing streams
incident on the inner horizon are proportional to $\text{\sl u} \mp \text{\sl v}$,
equation~(\ref{Npm}).
Inflation is driven by counter-streaming
between ingoing and outgoing streams,
so both streams must be present for inflation to occur,
but even the tiniest amount suffices to trigger inflation.
Positivity of both ingoing and outgoing densities requires that
\begin{equation}
  \text{\sl u} > \text{\sl v} > 0
  \ ,
\end{equation}
the condition $\text{\sl v} > 0$ coming from the fact that the black hole
must expand as it accretes.
The case $\text{\sl v} = 0$ is the stationary (or homogeneous)
approximation  of
\cite{Burko:1997xa}.
The densities $N^\pm$, equation~(\ref{Npm}),
are also proportional to $\Delta_x^\prime \mp \text{\sl v}$,
where $\Delta_x^\prime$ is the positive derivative of the electrovac
horizon function at the inner horizon.
Positivity of both ingoing and outgoing densities
requires $\Delta_x^\prime$ to be strictly positive,
which excludes extremal black holes,
whose inner and outer horizons coincide,
and which have $\Delta_x^\prime = 0$ at the horizon.

Solution of the evolution equations~(\ref{UDx})
for $U_x$ and $\Delta_x$ yields
\begin{subequations}
\label{UDxinf}
\begin{align}
  U_x
  &=
  \sqrt{
  \text{\sl v}^2
  +
  ( \text{\sl u}^2 - \text{\sl v}^2 )
  e^{4 \xi}
  }
  \ ,
\\
  \Delta_x
  &=
  -
  \left( {U_x^2 - \text{\sl v}^2 \over \text{\sl u}^2 - \text{\sl v}^2} \right)^{3/4}
  \left[
  {
  ( U_x + \text{\sl v} )
  ( \text{\sl u} - \text{\sl v} )
  \over
  ( U_x - \text{\sl v} )
  ( \text{\sl u} + \text{\sl v} )
  }
  \right]^{\Delta_x^\prime / ( 4 \text{\sl v} )}
  \ ,
\\
\label{xinf}
  x - x_-
  &=
  -
  \int
  {\Delta_x \, d U_x \over 2 ( U_x^2 - \text{\sl v}^2 )}
  \ .
\end{align}
\end{subequations}
The integral on the right hand side of equation~(\ref{xinf})
can be expressed analytically as an incomplete beta function,
but the expression is not useful.
Physically, equation~(\ref{xinf}) says
that the radius $x$ is frozen at its inner horizon value $x_-$
during inflation and collapse,
where $U_x$ is growing, while $\Delta_x$ remains small.

Figure~\ref{inflationaryhorizon}
illustrates the evolution of $U_x$ and the horizon function $\Delta_x$
as a function of the inflationary exponent $\xi$,
for parameters $\text{\sl v} = 0.001$ and $\text{\sl u} = 0.002$.
The value $\text{\sl v} = 0.001$,
which physically represents the velocity
with which a distant observer sees the characteristic radius
of the black horizon expand,
is large compared to a typical astronomical accretion rate,
but a large value is needed to avoid numerical overflow.
Inflation ignites near the inner horizon
as the horizon function $| \Delta_x |$ drops below $\text{\sl u}$.
During inflation,
the horizon function $| \Delta_x |$ decreases exponentially,
while $U_x$ increases slowly.
During inflation, the inflationary exponent $\xi$
in the conformal factor $\rho$, equation~(\ref{rho}),
satisfies
\begin{equation}
\label{expinfderivs}
  \xi
  \ll {d \xi \over d x}
  \ll {\dd^2\mkern-1.5mu \xi \over d x^2}
  \ ,
\end{equation}
which is the behaviour characteristic of a step-function.
The inequalities~(\ref{expinfderivs}) essentially say
that the acceleration $\dd^2\mkern-1.5mu \xi / d x^2$
of the inflationary exponent, which is driven by the
radial energy-momentum of the collisionless streams,
is much larger than the velocity $d \xi / d x$,
which in turn is much larger than the distance moved $\xi$.

Inflation ends when
the absolute value of the horizon function reaches a minimum,
at an exponentially tiny value,
\begin{equation}
  | \Delta_x | \sim e^{- 1 / \text{\sl v}}
  \ ,
\end{equation}
at which point the spacetime collapses.
During collapse,
the horizon function $| \Delta_x |$ increases,
while the conformal factor
$\rho \propto e^{- \xi}$
shrinks exponentially,
no longer satisfying the inequalities~(\ref{expinfderivs}).
The radial coordinate $x$ remains frozen
even while the conformal factor $\rho$ is shrinking.
That the spacetime collapses rather than leads to a null singularity
accords with the conclusion of
\cite{Hamilton:2008zz}
that the outcome of inflation is collapse
when a black holes continues to accrete,
as is ensured in the present case by the assumption of conformal time
translation invariance (self-similarity).

During collapse the horizon function increases back to of order unity,
$| \Delta_x | \sim 1$.
At this point the angular motion of the freely-falling ingoing and outgoing
streams becomes comparable to their radial motion,
and the solution breaks down.
This happens when the conformal factor has collapsed to an exponentially
tiny value,
\begin{equation}
  \rho \sim e^{- 1 / \text{\sl v}}
  \ .
\end{equation}

The right panel of Figure~\ref{inflationaryhorizon}
shows the magnitudes of the radial, radial-angular, and angular components
$T_{xx}$, $T_{xy}$, and $T_{yy}$
of the tetrad-frame energy-momenta of the collisionless streams.
During inflation,
the radial energy-momentum grows fastest,
reaching an exponentially huge value
\begin{equation}
  T_{xx} \sim e^{1 / \text{\sl v}}
  \ .
\end{equation}
During collapse,
the angular energy-momentum grows fastest.

The Weyl curvature tensor has only a spin-$0$ component,
which classifies the spacetime as Petrov Type~D.
The right panel of Figure~\ref{inflationaryhorizon}
shows minus the polar (real) part
of the spin-$0$ component of the Weyl curvature.

It is notable that
the smaller the accretion rate $\text{\sl v}$,
the more rapidly inflation exponentiates,
and the larger the energy-momentum and curvature grow,
in agreement with the conclusions of
\cite{Hamilton:2004aw,Hamilton:2008zz}.

\section{Boundary conditions}

The solutions are determined by boundary conditions of the
collisionless streams incident on the inner horizon.
Since the solution above the inner horizon is well-approximated
by the Kerr (or other electrovac) solution,
the behaviour of gas above the inner horizon does not affect the solution.

The requirement of conformal separability imposes special boundary conditions.
The densities $N^\pm$ of ingoing ($+$) and outgoing ($-$)
streams incident on the inner horizon are, equation~(\ref{Npm}),
since $U_x = \text{\sl u}$ initially,
\begin{equation}
\label{Npminit}
  N^\pm
  =
  {1 \over 16\pi}
  ( \text{\sl u} \mp \text{\sl v} )
  ( \Delta_x^\prime \pm \text{\sl v} )
  \ .
\end{equation}
This is just a constant, independent of angular position on the inner horizon.
Thus conformal separability requires that the incident flow
of ingoing and outgoing streams be ``monopole,''
independent of latitude.
It makes physical sense that conformal separability
would require this high degree of symmetry of the accretion flow.

As emphasized in \S\ref{derivation},
because the radial motions of collisionless streams dominate
their angular motions during inflation and collapse
(up until the angular motions become important, at $| \Delta_x | \sim 1$),
the radial Einstein equations are unaffected
by the angular motion, and the boundary condition~(\ref{Npminit})
is all that is needed to ensure conformal separability
with sufficient accuracy.
However, if it is required that the sub-dominant radial-angular
components of the Einstein equations are also satisfied,
which is a more stringent constraint on conformal separability,
then the tetrad-frame momenta $p^\pm_k$ of the streams
must have finite angular components,
satisfying equation~(\ref{ppm}).
Figure~\ref{rotinforbit} illustrates the required flow pattern
for a black hole of spin parameter $a = 0.96 M_\bullet$,
The energy per unit mass of infalling particles,
which is unspecified by boundary conditions,
is chosen here to be $E/m = \pm 1$.

\rotinforbitfig

As described in Paper~2,
the required angular flow pattern cannot be achieved with collisionless
streams falling from outside the outer horizon.
Streams that fall from outside the outer horizon
must necessarily be ingoing at the outer horizon,
which eliminates half the phase space available to collisionless streams,
making it impossible to satisfy the required conditions
on the angular motion of the streams.
Thus, if the angular conditions are imposed,
then the collisionless ingoing and outgoing streams must be considered
as being delivered ad hoc to just above the inner horizon.
This is a defect, but not fatal.
Firstly, the collisionless assumption is likely to fail above the
inner horizon, where baryons are collisional.
Secondly,
the angular conditions are not essential to
the radial part of the solution,
which  holds provided that the monopole condition~(\ref{Npminit}) is satisfied.
If the angular conditions are relaxed,
then the required boundary conditions~(\ref{Npminit})
can easily be satisfied by collisionless streams
falling from outside the outer horizon.
Although particles are necessarily ingoing at the outer horizon,
they can switch from ingoing to outgoing between the inner and outer horizons
if they have a large enough angular momentum
in the same direction as the black hole.

\section{Conclusion}

This paper has presented conformally stationary, axisymmetric,
conformally separable solutions for the interior structure of an uncharged
rotating black hole
that undergoes inflation just above its inner horizon, then collapses.
It has long been known that linear perturbations diverge at
the inner horizon of the Kerr geometry,
and it has been suspected that perturbations
would develop nonlinearly similarly to the inflationary instability
\cite{Poisson:1990eh}
known to operate in spherical charged black holes
\cite{Hamilton:2008zz}.
The self-consistent nonlinear solutions found here
confirm that,
at least in the conformally separable special case considered here,
the inflationary instability develops
in rotating black holes as anticipated.

A feature of the Kerr geometry (and other separable electrovac geometries)
is that, as
freely-falling ingoing and outgoing particles approach the inner horizon,
they become highly focussed along the ingoing and outgoing null directions,
regardless of their initial angular motion.
However small the accretion rate may be,
eventually the energy and pressure of the twin beams of particles
counter-streaming hyper-relativistically along the ingoing and outgoing
null directions grows large enough to be a source of gravity.
The gravity produced by the counter-streaming acts to accelerate
ingoing and outgoing stream even faster through each other,
leading to an exponential growth in the streaming density and pressure,
and in the curvature.
This is inflation.
The huge gravitational acceleration produced by the counter-streaming
is in the inward direction, to smaller radius,
but each stream thinks that they are moving in the inward direction,
so the streams are accelerated in opposite directions.

Inflation takes place over an extremely short interval of proper time.
Inflation is like a bullet fired in the chamber of a gun:
an explosion accelerates the bullet,
and shortly after the bullet achieves high velocity,
but still the bullet has hardly moved
(see the inequalities~(\ref{expinfderivs})).
Inflation does in due course alter the geometry,
but in a predictable way:
the conformal factor,
having been accelerated to huge velocity inward,
proceeds to shrink rapidly.
The geometry collapses.

During inflation,
the ingoing and outgoing streams were accelerated along the principal
null directions, without amplifying the angular motion.
During collapse, the angular motions grow.
At an exponentially tiny scale,
the angular motions become comparable to the radial motion,
and the solutions considered in this paper break down.
What happens then is undetermined.

The existence of conformally separable solutions for
the inflationary zones of rotating black holes is not surprising.
Ingoing and outgoing streams focus along
the principal null directions as they approach the inner horizon,
and the streaming energy and pressure generated by the radial beams
accelerates the streams along the same null directions.
The acceleration depends on the accretion rate,
being faster for smaller accretion rates.
If the densities of ingoing and outgoing streams incident on
the inner horizon are uniform, independent of latitude,
then inflation accelerates the beams at the same rate at all latitudes.
When the geometry begins to collapse, it does so uniformly,
preserving conformal separability.
What happens in the more general case when the accretion flow
on to the inner horizon varies with angular position
remains to be seen.
But that inflation will occur is physically inevitable.

\begin{acknowledgements}
This work was supported by NSF award
AST-0708607.
\end{acknowledgements}

\section*{References}

\bibliographystyle{unsrt}


\section{Introduction}
Optimal control of molecular alignment is crucial for applications of
liquid crystal films.
Thus far, many sophisticated techniques  have been developed for obtaining
homogeneous alignment in a desired direction: mechanical rubbing treatment
\cite{berreman,fukuda1,fukuda2},
surface-polarization effect \cite{barbero1,barbero2},
surface coating \cite{nesrullajev} and
photoalignment \cite{ichimura,hwang,ho,kawatsuki} are only a few examples
to mention.
When employing these techniques, topological defects are regarded as
unwanted matters and thus they are eliminated artificially.
Under certain conditions, however, topological defects can have beneficial
effects on the alignment control of liquid crystal molecules
\cite{chiccoli1,berggren,chiccoli2,nelson1,vitelli1,santangelo1,fernandez,shin,xing,lopezleon}
as well as of guest particles embedded in a liquid
crystal medium \cite{ruan,musevic,yoon,skacej,guo,fukuda3,ravnik,bates}.
The positive effects result largely from topological constraints for the
molecular orientational field ($=$ director field) that is confined in a
closed curved substrate.
In fact, the director field over closed curved surfaces strongly depend
on the positions and charges of topological defects distributed on the
surfaces, which implies the controllability of the director field through
the defect manipulation.
Such the utilities of topological defects may be enhanced when combined
with the geometric potential effect.
It is known that topological defects in curved liquid crystal films
experience a surface-curvature-induced potential, attractive or repulsive
depending on the local geometry of the surface
\cite{nelson2,seung,park,vitelli2,vitelli3,gomez}
\footnote{Similar curvature effects were discussed in the literature of
quantum transport \cite{shima1,ono} and classical spin models \cite{baek} on
curved surfaces.}.
Therefore, a better understanding of the correlation between curvature
variation of the underlying surface and favourable configurations of
topological defects will shed new light on molecular alignment control of
liquid crystal films.

A Gaussian bump is one of the exemplary curved surfaces for studying the
correlation mentioned above.
It has a continuous rotation symmetry around a vertical axis, showing a
positive (negative) Gaussian curvature close to (far from) the top of the
bump [see Fig. \ref{fig:geometry}(a)].
The coexistence of a positively curved region and a negatively one on the
same surface allows a simultaneous observation of the curvature effects at
the two distinct regions with oppositely-signed Gaussian curvature.
Earlier work suggested that when a nematic liquid crystal film is placed on
a Gaussian bump, a point disclination with $+1$ charge is trapped to the
top of the bump due to the curvature potential
\cite{vitelli2,vitelli3}.
Contrariwise, we have found in the previous work that point disclinations
are attracted to the inflection point of the bump regardless of their
charges \cite{isaku}.
These two results appear controversial; a possible clue to resolve the
controversy is the difference in the assumed bump widths between the two
existing studies.
In the former study, the bump width is sufficiently small, so that the
bump is covered almost entirely by a single disclination as depicted in
Fig. \ref{fig:geometry}(c).
This situation is expressed by saying that the coherence length
\cite{kleman1,wong}, $\xi$, of the nematic order fluctuation exceeds the
scale of the bump width (see \S \ref{sec:method} for the definition of
$\xi$).
In contrast, the latter study was based on the assumption that the bump
width is larger enough than the coherence length and the system can involve
many small disclinations [see Fig. \ref{fig:geometry}(b)].
It is thus reasonable to conjecture that the ratio of the bump width to the
coherence length plays a decisive role in determining the stable position of
disclinations on a Gaussian bump.
Furthermore, the two contrasting results imply the existence of a
threshold ratio across which the stable position of disclinations alters.

In the present Letter, we prove our conjecture by showing a sudden shift
in the stable position of a $+1$ disclination with increasing the bump width.
Once the bump width exceeds a threshold value, the stable position is
suddenly switched from the top to the inflection point of the bump.
The positional shift is described by a power law function of the bump width,
as is analogous to the second-order phase transitions of various physical
systems whose properties are characterized by power-law behaviors.
The present result thus indicates a novel kind of phase transition that
is peculiar to curved nematic films.

\begin{figure}[t]
\begin{center}
\includegraphics[width=340pt]{Fig1.eps}
\end{center}
\caption{\label{fig:geometry} (a) Sketch of a Gaussian bump characterized
by the width $\Delta$ and the hight $v$.
The region inside (outside) the dashed circle with radius $\Delta$ possesses
a positive (negative) Gaussian curvature.
(b, c) Disclination embedded in the Gaussian bump.
The coherence length $\xi$ of the nematic order fluctuation measures the
spatial extent of desclinations (hatched regions) in which the director
field $\bm n({\bm R})$ is arranged in a radial manner.}
\end{figure}

\section{Method}
\label{sec:method}
We consider point disclinations in a nematic liquid crystal membrane having a
Gaussian bump geometry.
The bump is described by a position vector
${\bm R(r,\phi)} = r \cos \phi \bm{i} + r \sin \phi \bm{j}
                     + v \exp (-r^2/2 \Delta^2) \bm{k}$,
where $\bm{i}$, $\bm{j}$, $\bm{k}$ are right-handed orthogonal basis vectors
[see Fig. \ref{fig:geometry}(a)].
The constant $\Delta$, which we call the bump width, determines the
inflection point of the bump at
which $(\partial^2 {\bm R} / {\partial r}^2) \cdot {\bm k} = 0$.
A dimensionless parameter $\alpha \equiv v/\Delta$ quantifies the
convexity of the bump; $\alpha = 0$ for the flat plane, and $\alpha \gg 1$
for a spiky bump.
We fix $\alpha \equiv 0.1$ throughout the paper, as a result of which
spatial variation of the surface curvature is smooth and sufficiently
slight to be accessible in experiments.

Our aim is to evaluate the stable position of a disclination constrained on
the Gaussian bump.
It is solved by calculating the elastic free energy $F_{\rm d}$ of the
director field ${\bm n}({\bm R})$ around the disclination core, since
$F_{\rm d}$ takes a minimum at the stable position to be evaluated.
We assign the core of a $+1$ disclination at
${\bm R}_{\rm c} = {\bm R(r_{\rm c},\phi_{\rm c})}$ with arbitrary
$\phi_{\rm c}$, and assume that the disclination is endowed with the coherence
length $\xi$ of the nematic order fluctuation.
The latter assumption ensures a radial pattern of $\bm n({\bm R})$ within
a circular region $A$ having the radius $\xi$ centered at the core [i.e.,
the hatched regions depicted in Figs. \ref{fig:geometry}(b) and
\ref{fig:geometry}(c)] \footnote{Precisely speaking, the area $A$ slightly
deviates from an exact circle due to the finite curvature of the surface.
In fact, we define the perimeter of $A$ by the locus of points whose
distances from the core are equal to $\xi$ in the sense of geodesic distance.
Therefore, the perimeter is distorted a bit from the exact circle.}.
The radial pattern of $\bm n({\bm R})$ is gradually distorted with increasing
distance from the core; in fact, two-point correlation in the director
orientation disappears if the points are well separated by a distance larger
than $\xi$.
A typical value of $\xi$ is on the order of tens micrometers as can be
estimated from polarization microscope images of liquid crystal films
\cite{nehring1,shiwaku}.

The continuum elastic theory states that $F_{\rm d}$ associated with the
disclination reads as \cite{santangelo2,kleman2}
\begin{eqnarray}
F_{\rm d} &=& h \int_A dA \left\{ \frac{K_1}{2} ({\bm D} \cdot {\bm n})^2
+ \frac{K_3}{2} ({\bm D} \times {\bm n})^2 + \frac{K_{24}}{2}
{\bm D} \cdot [({\bm n} \cdot {\bm D}) {\bm n} -
{\bm n} ({\bm D} \cdot {\bm n})] \right\}
\mbox{,}
\label{eq:energy}
\end{eqnarray}
where $\bm D$ is a vector operator whose components $D_{\mu}$ are the
covariant derivative \cite{shima2} on the curved surface.
The constant $h$ is the thickness of the liquid crystal film; $K_1$,
$K_3$ and $K_{24}$ are the elastic constants associated with splay, bend
and saddle-splay distortions of the field $\bm n$, respectively
\cite{stewart}.
The integration range in Eq. (\ref{eq:energy}) is restricted to the finite
circular region $A$, despite of slow decay in the elastic energy density with
distance from the core; we will revisit this issue in \S \ref{sec:discussions}.
To avoid the divergence in $F_{\rm d}$ at the core, we introduce an inner
cut-off $A_{\rm {core}}$ around the core according to the formula
$A_{\rm {core}} = \pi K_1/2 k_{\rm B} (T_{\rm c} - T)$, in which
$T_c$ is the nematic-isotropic transition temperature and $T$ is the
temperature of the system \cite{cladis}.

In actual calculations, we set $\xi = 20 {\rm {\mu m}}$
and $h = 4 \  {\rm {\mu m}}$ by referring to the experimental
observation \cite{wong,shah}.
The values of $K_1$ and $K_3$ are those of
4-methoxybenzylidene-4$^{\prime}$-butylaniline (MBBA): $K_1 = 6$ pN,
$K_3 = 7.5$ pN at room temperature \cite{stewart}.
Since the corresponding value of $K_{24}$ was not yet measured, we set
$K_{24} = 4.9$ pN on the basis of the relation $K_{24} = (K_1 + K_2)/2$
that involves the twist coefficient $K_2 = 3.8$ pN for MBBA \cite{nehring2}.

\begin{figure*}[t]
\begin{center}
\includegraphics[width=330pt]{Fig2.eps}
\end{center}
\caption{\label{fig:energyplot} (a) Disclination energy $F_{\rm d}$ in
deformed MBBA films as a function of the core position $r_{\rm c}$.
(b) Sudden shift in the stable position $r_{\rm c}^{\rm min}$ as a function
of the bump width $\Delta$.
The threshold $\Delta^{\rm th}$ is estimated to be
$\Delta^{\rm th}/\xi = 0.95$.
Inset: The same data in the log-log plot, showing a power law
$r_{\rm c}^{\rm min} \propto \left| \Delta - \Delta^{\rm th} \right|^{\beta}$
with the exponent $\beta \sim 0.2$.}
\end{figure*}

\section{Results}
Figure \ref{fig:energyplot}(a) shows the disclination energy $F_{\rm d}$
as a function of the core position $r_{\rm c}$.
The normalized bump width, $\Delta/\xi$, ranges from $0.8$ to $2.0$.
For $\Delta/\xi = 2.0$, the curve of $F_{\rm d}$ is convex downward, showing a
minimum at $r_{\rm c} \simeq \Delta$.
The energy scale of the potential depth corresponds to ${k_{\rm B} T}$ at a
temperature of $T = 10^4 {\rm K}$, which is by far larger than the thermal
excitation energy at room temperature.
Hence, this downward peak implies that a disclination is
trapped to the annulus region having the radius of $r = \Delta$
surrounding the bump top [see Fig. \ref{fig:halocon}(a)].
We call this annular distribution of disclinations by the
\lq\lq halo\rq\rq\ phase,
according to the discussion in Ref. \cite{isaku}.

What deserves attention is the $\Delta$-driven transition of the downward
peak profile.
It follows from Fig. \ref{fig:energyplot}(a) that both the position and
magnitude of the peak are insensitive to the decrease in $\Delta/\xi$ until
it reaches a specific value of $\Delta/\xi = 0.95$.
A drastic change takes place at $\Delta/\xi = 0.95$; the downward peak in
$F_{\rm d}$ is replaced by a plateau that extends from $r_{\rm c} = 0$ to
$r_{\rm c} \sim 0.5 \Delta$.
This plateau means that the disclination can move freely within a finite
circular region centered at the bump top.
It is also found that a further decrease in $\Delta/\xi$ gives rise to a new
downward peak in $F_{\rm d}$ at $r_{\rm c} = 0$.
That is, when the bump width $\Delta$ is comparable to (or smaller than) the
coherence length $\xi$, then the disclination is trapped at the bump top
[see Fig. \ref{fig:halocon}(b)]; this result is in complete agreement with
the earlier observation \cite{vitelli2,vitelli3}.
We call this defect configuration by the \lq\lq convergence\rq\rq\ phase.

We have seen that the peak of $F_{\rm d}$ appearing at
$r_{\rm c} = r_{\rm c}^{\rm min}$ is
gradually shifted from the inflection point
$(r_{\rm c}^{\rm min} \sim \Delta)$ to the top
($r_{\rm c}^{\rm min} = 0$) of the bump with increasing $\Delta$.
Figure 2(b) shows the $\Delta$-dependence of the normalized
peak position $r_{\rm c}^{\rm min}/\Delta$.
For a large bump width ($\Delta/\xi \gg 1$),
$F_{\rm d}$ takes the minimum at $r_{\rm c}^{\rm min}/\Delta \sim 1$
implying the halo phase.
The disclination halo shrinks with decreasing $\Delta$,
and finally it converges at the bump top when
$\Delta = \Delta^{\rm th} \equiv 0.95 \xi$.
These results indicate the continuous phase transition
(i.e., of the second order)
of the stable disclination configuration on the Gaussian bump;
the system undergoes the transition from the halo phase
to the convergence phase across the threshold $\Delta^{\rm th}/\xi$.
This is the main finding of the present Letter.

It is interesting to note that the $\Delta$-dependence of
$r_{\rm c}^{\rm min}$ in the halo phase is described by a power law
of the form
\begin{equation}
r_{\rm c}^{\rm min} \propto \left| \Delta - \Delta^{\rm th} \right|^{\beta}.
\label{eq_shima}
\end{equation}
The inset of Fig.~2(b) evidences the power-law behavior of
Eq.~(\ref{eq_shima})
with the exponent $\beta \sim 0.2$,
which holds over the two orders of the lateral axis.
We thus believe that the present result is a precursor of an
untouched configurational phase transition peculiar to curved liquid
crystal systems, though further studies are needed to establish its
universality and physical interpretation.

\section{Discussions}
\label{sec:discussions}
The transition of the stable defect configuration
we have demonstrated is a consequence of the following two
$F_{\rm d}$ variations caused by reducing $\Delta$: a slight
$F_{\rm d}$-increase to the right of the arrow $P$ marked in Fig.~2(a),
and a rapid $F_{\rm d}$-decrease to the left of $P$.
In this section, we give an outline of the physical mechanisms
of the two $F_{\rm d}$ variations, while a detailed formulation
that describes them will be given elsewhere.

The rapid decrease in $F_{\rm d}$ to the left of $P$ originates from
a suppression of the splay deformation of the director field
${\bm n} (\bm{R})$ around the bump top.
Suppose that a $+1$ disclination is placed near the top
of a broad bump as depicted in Fig.~\ref{fig:geometry}(c).
Since $\xi$ and $\alpha = v/\Delta$ are constants, a reduction in $\Delta$
leads to a lowering of the bump hight $v$ remaining the disclination area
$A$ be fixed.
As a result, the disclination is allowed to cover the whole bump.
In the latter situation, adjacent directors (thick line segments in the
plots) distributed close to the outer edge of $A$ tend to be parallel to
each other.
Eventually, the first integrated term in Eq. (\ref{eq:energy}),
$\propto (\bm{D}\cdot\bm{n})^2$,
which represents the contribution from the splay deformation of $\bm{n}$,
is suppressed by reducing $\Delta$.
This is the reason why $F_{\rm d}$ at $r_{\rm c}=0$ drops off when
$\Delta$ decreases.

A parallel argument to the above accounts for the slight increase
in $F_{\rm d}$ to the right of $P$ in Fig.~\ref{fig:energyplot}(a).
Suppose a $+1$ disclination on the skirt (i.e., at $r_{\rm c} \sim \Delta$)
of the bump.
If the disclination area $A$ intersects the circle $r=\Delta$ around the
bump top,
then it covers a positively-curved region ($r<\Delta$)
and a negatively-curved one ($r>\Delta$).
At the negatively curved region, a reduction in $\Delta$ yields an
enhancement of the splay deformation which leads to the slight increase
in $F_{\rm d}$ at $r\sim \Delta$.
A mathematical proof of the splay enhancement has been given in
Ref. \cite{isaku}.

The data shown in Figs. \ref{fig:energyplot}(a) and \ref{fig:energyplot}(b)
are derived from a specific Gaussian bump shape characterized by the parameter
$\alpha = v/\Delta \equiv 0.1$ (see Fig. \ref{fig:geometry}(a)).
We also confirmed that, when $0.05 < \alpha < 0.5$, the system exhibits the
same configurational transition across $\Delta^{\rm th}$.
A slight increase in $\Delta^{\rm th}$ with $\alpha$ was observed within the range
of $\alpha$, reflecting further splay suppressions and enhancements in
disclinations whose cores locate at $r_{\rm c} = 0$ and $r_{\rm c} = \Delta$,
respectively.
In contrast, the value of $\beta$ is insensitive to the change in $\alpha$.

It should be stated that our discussion is based on the model given by
Eq. (\ref{eq:energy}), in which contributions from the outer region than the
circle $A$ are ignored.
In actual nematic layers, a $+1$ disclination has a profile of the elastic
energy density that decays logarithmically with distance from the core.
Hence, contributions far from the core may be relevant to the stable position
of the disclination, thus should be taken in account in order to obtain a
quantitatively accurate value of $\Delta^{\rm th}$.
The geometric potential approach developed in Ref. \cite{vitelli4} will be a
clue to solve this issue.

Before closing the Letter, we suggest the following system could allow
experimental tests of our theoretical predictions; it is a thin liquid crystal
film coating a randomly corrugated substrate.
The film involves many disclinations and bumps with various sizes; therefore,
we will observe simultaneously the two effects, halo and convergence patterns
depicted in Fig. \ref{fig:halocon}, by tuning the surface corrugation.
We also mention that such the corrugated substrate can be synthesized by polymer
adhesion with patterned surfaces as demonstrated in Ref. \cite{chan}.

\begin{figure*}[t]
\begin{center}
\includegraphics[width=350pt]{Fig3.eps}
\end{center}
\caption{\label{fig:halocon} Distribution patterns of disclinations in the
halo phase (a) and convergence phase (b).
In either phase, disclinations are trapped in the dark region.}
\end{figure*}

\section{Conclusion}
We have found a continuous transition of defect configurations in nematic
films confined on a Gaussian bump.
Depending on the ratio $\Delta/\xi$, the system exhibits either of two
defect distribution patterns: the halo phase ($\Delta/\xi>1$) and
convergence phase ($\Delta/\xi \le 1$).
The stable position $r_{\rm c}^{\rm min}$ of defects turned out to obey the
power law
$r_{\rm c}^{\rm min} \propto \left| \Delta - \Delta^{\rm th} \right|^{\beta}$
with the exponent $\beta \sim 0.2$ under the present condition.
Our findings suggest the possibility of defect manipulation via artificial
deformation of liquid crystal films.

\section*{Acknowledgments}
We would like to thank K.~Yakubo, H. Orihara, J. Fukuda and Y. H. Na for
fruitful discussions.
IH is thankful for the financial support from the Japan Society for the
Promotion of Science for Young scientists.
HS acknowledges M. Arroyo for his hospitality during the stay in UPC.
This work is supported by the Kazima foundation and a Grant-in-Aid for
Scientific Research from the Japan Ministry of Education, Science, Sports and
Culture.