\section{Introduction}
Let $\Sigma \subset \mathbb{R}^3$ be a smooth, oriented minimal surface. In this paper, we introduce a new meromorphic quadratic differential $P$ on $\Sigma$, which we call the \emph{entropy differential}. We use $P$ to characterize several classical surfaces -- including Enneper's surface, the catenoid and the helicoid. In particular, subsets of Enneper's surface are the only minimal surfaces on which $P$ vanishes -- a fact which we use to prove a novel curvature estimate for embedded minimal surfaces with small entropy differential and an associated compactness result.

The differential $P$ arises naturally in a number of different contexts.  Of particular interest is the realization of $T=\re P$, which we call the \emph{entropy form}, as a conservation law for the diffeomorphism invariant functional
\begin{equation*}
 \mathcal{E}[g]=\int_{\Sigma} K_g \log K_g \mu_g.
\end{equation*}
 This functional, which is a type of entropy for the curvature has been previously considered by R. Hamilton in the context of the Ricci flow on surfaces ~\cite{Hamilton2Sphere}.  In particular, we show that if $g$ is a minimal surface metric (i.e. the metric induced by a smooth minimal immersion) for which $K_g\neq 0$, then the metric $\hat{g}=(-K_g)^{3/4} g$ is a critical point of $\mathcal{E}$ with respect to compactly supported conformal deformations.  The crucial fact used here is the observation -- due to Ricci~\cite{Ricci} -- that such minimal surface metrics satisfy the so-called \emph{Ricci condition}:
\begin{equation*}
 \Delta_g \log |K_g|=4K_g.
\end{equation*}
The differential $P$ also arises as certain geometric Schwarzian derivative of the Gauss map -- a point of view we discuss in Section \ref{WeierstrassSec} and which has antecedents in work of~\cite{DurenEtAl} and~\cite{KusnerConversation}.  Moreover, $P$ seems to arise naturally when studying minimal surfaces from the point of view of integrable systems -- we will investigate this connection more thoroughly in a forthcoming paper~\cite{BernsteinMettler}.

A key observation of the present paper is that, modulo rigid motions, a minimal surface $\Sigma$ is determined, up to a three-parameter family, by its Hopf differential $Q$ and its entropy differential $P$.  This allows one to characterize several natural minimal surfaces in terms of simple relationships between the Hopf differential and entropy differential:
\setcounter{section}{4}
\begin{thm}\label{IntroMainCharacterizationThm}
 Let $\Sigma$ be a smooth oriented non-flat minimal surface in $\mathbb R^3$ with entropy differential $P$. We have:
\begin{enumerate}
 \item \label{IntroEnneperCase} If $P\equiv 0$, then up to a rigid motion and homothety, $\Sigma$ is contained in Enneper's surface;
 \item \label{IntroCatCase} If $\lambda \neq 0$ and  $P\equiv \lambda Q$, then,  up to a rigid motion and homothety, $\Sigma$ is contained in a surface $C\in\mathcal{C}$.  If $\Sigma$ is properly embedded, then it is the catenoid;
\item \label{IntroHelCase} If $\lambda \neq 0$ and $P \equiv  i \lambda Q$, then, up to a rigid motion and homothety, $\Sigma$ is contained  in a surface $H\in \mathcal{H}$.  If $\Sigma$ is properly embedded, then it is the helicoid.
\end{enumerate}
\end{thm}
The families $\mathcal{C}$ and $\mathcal{H}$ are, respectively, the deformed catenoids and deformed helicoids.  These are one parameter families of surfaces containing, respectively, the catenoid and the helicoid -- their geometry is discussed thoroughly in Section \ref{CharSec}.

A consequence of Item \eqref{IntroEnneperCase} of Theorem \ref{IntroMainCharacterizationThm} is a novel curvature estimate for embedded minimal surfaces.  Namely, we introduce a certain natural scale invariant norm on the entropy form and use standard blow-up arguments to derive a curvature bound and a corresponding compactness theorem.  Specifically, for a smooth minimal surface $\Sigma$ with entropy form $T$, we define the weighted $L^{1/2}$ norm of $T$ by:
\begin{equation*}
||T||_{\hat{L}^{1/2}(\Sigma)}^{1/2}:=\int_{\Sigma} |T|^{1/2}_g \sqrt{-K_g} \mu_g.
\end{equation*}
This norm is scale invariant, which allows us to use blow-up arguments to obtain the following curvature estimate:
\setcounter{section}{5}
\begin{thm}\label{IntroCurvEstThm}
 There exist constants  $\epsilon>0$ and $C>0$ so that: if $\Sigma$ is a properly embedded minimal surface in $B_{2R}$ and
$$||T||_{\hat{L}^{1/2}(\Sigma)}<\epsilon,$$
 then
\begin{equation*}
 R^2 \sup_{B_{R}\cap \Sigma} |A|^2 \leq C^2.
\end{equation*}
\end{thm}
In Section \ref{GapSec}, we address the question of the best possible $\epsilon$ and some partial results are obtained.
Such a curvature estimate gives, using standard techniques, a compactness result:
\setcounter{thm}{7}
\begin{thm}
 Let $\Sigma_i$ be a sequence of properly embedded minimal surfaces in an open region $\Omega\subset \mathbb R^3$ with entropy forms $T^{\Sigma_i}$ satisfying
$$||T^{\Sigma_i}||_{\hat{L}^{1/2}(\Sigma_i)}\leq \bar{C}<\infty.$$
 Then there is a subsequence of the $\Sigma_i$ and  a finite (possibly empty) set of points $p_1, \ldots, p_N\in \Omega$ so that:
\setcounter{thm}{1}
\begin{enumerate}
  \item \label{IntroCpctItemOne} On each compact set $K\subset \subset \Omega\backslash \set{p_1, \ldots, p_N},$
     $$\sup_{K\cap \Sigma_i} |A|\leq C(K)<\infty;$$
\item\label{IntroCpctItemTwo}  $\epsilon N^2<\bar{C}$ where $\epsilon>0$ is given by Theorem \ref{IntroCurvEstThm};
\item \label{IntroCpctItemThree} The $\Sigma_i$ converge in $\Omega\backslash \set{p_1, \ldots, p_N}$ to a smooth minimal lamination $\mathcal{L}$ of $\Omega\backslash \set{p_1, \ldots, p_N}$.  Moreover, the closure of $\overline{\mathcal{L}}$ of $\mathcal{L}$ in $\Omega$ is a smooth lamination of $\Omega$.
\end{enumerate}
\end{thm}
We conclude the paper with Appendix \ref{RicciFlowSec}, wherein the entropy form is used to make a connection between minimal surfaces in $\mathbb R^3$ and gradient Ricci soliton metrics on surfaces.

\subsection*{Acknowledgements}
The authors would like to thank Rob Kusner for several stimulating discussions regarding the topics of this paper. The second author also would like to thank Daniel Fox for interesting discussions.

\setcounter{section}{1}

\section{A Geometric Entropy Functional for Surfaces}

\subsection{Definitions}
We assume $\mathbb{R}^3$ to be equipped with the standard Euclidean metric $g_E$ and orientation. Let $M$ be an open orientable smooth $2$-manifold. For a smooth immersion $\mathbf{x} : M \to \mathbb{R}^3$ let $\Sigma=\mathbf{x}(M)$ and $g=\mathbf{x}^*g_E$ be the first fundamental form. We write $\gnabla$ for the Levi-Civita connection, $K_g$ for the Gauss curvature, and $\mu_g$ for the area form of $g$. The integrable almost complex structure on $M$ induced by $g$ and the orientation will be denoted by $J$. We will use the same letter $J$ to denote the induced map on the space $\Gamma(S^2(T^*M))$ of smooth sections of $S^2(T^*M) \to M$, i.e. for $A \in \Gamma(S^2(T^*M))$, we define
$
(JA)(X,Y)=A(JX,Y)
$
where $X,Y \in \Gamma(TM)$ are smooth vector fields on $M$. Here, as usual, $S^2(T^*M)$ denotes the second symmetric power of the cotangent bundle of $M$. In particular, the map $A\mapsto A+i JA$ embeds the space of symmetric trace-free $2$-forms on $M$ into the space of quadratic differentials on $M$. Furthermore, it is easy to check that $A+i JA$ is holomorphic if and only if $A$ is divergence-free.

Let $\mathbf{n}$ denote the orientation compatible Gauss map of $\mathbf{x}$ taking values in $\partial B_1$, the unit-sphere in $\mathbb{R}^3$ centered at $0$. The second fundamental form of $\mathbf{x}$ will be denoted by $A$ and its trace with respect to $g$, the \textit{mean curvature}, by $H$. A point $p\in M$ at which the eigenvalues of $A_p$ agree is called \emph{umbilic} and we define $\hat{M}\subset M$ to be the open submanifold of non-umbilic points.

Recall that the Bonnet theorem states that a pair $(g,A)$ on a simply connected $2$-manifold $N$ arises as the first and second fundamental form via an immersion $\mathbf{x} : N \to \mathbb{R}^3$
if and only if the pair satisfies
the Gauss equations
\begin{equation*}
 |A|_g^2+2K_g =(\tr_g\! A)^2
\end{equation*}
and the Codazzi equations
\begin{equation*}
 \gnabla_X A (Y,Z)=\gnabla_Y A(X,Z)
\end{equation*}
where $X,Y,Z\in \Gamma(TN)$. Furthermore, $\mathbf{x}$ is unique up to composition with a rigid motion of $\mathbb{R}^3$. For this reason we refer to the triple $(M,g,A)$ as \emph{geometric data} of $\mathbf{x}$.

\subsection{The Ricci condition} We suppose from now on that $\mathbf{x} : M \to \mathbb{R}^3$ is minimal.  That is $H\equiv 0$. The Gauss equations imply that $K_g\leq 0$ and that $K_g$ is negative on $\hat{M}$. It follows with the Codazzi equations that the second fundamental form is divergence free with respect to $g$. This yields \textit{Simons' identity}
\begin{equation*}
 \Delta_g A =-|A|_g^2 A
\end{equation*}
where here $\Delta_g$ is the rough Laplacian. From
$$
4 |A|_g^2 |\gnabla A|^2_g=|\gnabla |A|_g^2|^2,
$$
Simons' identity, and the Gauss equations we obtain that on $\hat{M}$ the following \textit{Ricci condition}
\begin{equation} \label{RicIdent}
 \Delta_g \log |{K}_g| = 4 K_g
\end{equation}
holds. Abbreviate $u_g=-\frac{1}{4}\log|K_g|$, then the Ricci condition becomes
$$
\Delta_gu_g=e^{-4u_g}.
$$

Conversely, Ricci~\cite{Ricci} showed that if $g$ is a Riemannian metric of strictly negative Gauss curvature $K_g$ on a simply connected $2$-manifold $N$ satisfying the Ricci condition, then there exists a minimal immersion $\mathbf{x} : N \to \mathbb{R}^3$, unique up to rigid motion, with $\mathbf{x}^*g_E=g$. A proof of this fact using modern language may be found in~\cite{MR655419}.

\subsection{The entropy functional}
We will study a certain natural functional $\mathcal{E}$ defined  on the space $\mathcal{M}^+(M)$ of smooth Riemannian metrics on $M$ which have positive Gauss curvature.
Define for $g\in \mathcal{M}^+(M)$
\begin{equation*}
\mathcal{E}[g]= \int_M K_g \log K_g \mu_g
\end{equation*}
This functional has been applied to the study of Ricci flow on surfaces by Hamilton~\cite{Hamilton2Sphere} and Chow~\cite{Chow1991} -- in particular Hamilton observed that it is monotonically increasing along the Ricci flow on spheres with positive Gauss curvature.  See Appendix~\ref{RicciFlowSec} of the present paper for additional connections to Ricci solitons.

We compute the Euler-Lagrange equations associated to $\mathcal{E}$.
To do so let $f$ be a smooth symmetric $2$-form and write $F=\tr_g f$ for its trace with respect to $g$. Let $g_t=g+tf$, then (cf.~\cite[pg. 99]{ChowBook})
\begin{align*}
 \left.\frac{\partial}{\partial t}\right|_{t=0} K_{g_t} &= - \frac{1}{2} \Delta_g F + \frac{1}{2} \Div_g(\Div_g\! f)- \frac{1}{2} K_g F\\
&= - \frac{1}{4} \Delta_g F + \frac{1}{2} \Div_g\left(\Div_g\! \mathring{f}\right)- \frac{1}{2} K_g F
\end{align*}
where $\mathring{f}$ is the trace-free part of $f$.
Hence,
\begin{equation*}
\left.\frac{\partial}{\partial t}\right|_{t=0} K_{g_t} \log K_{g_t}  =  \left(-\frac{1}{4} \Delta_g F + \frac{1}{2} \Div_g\left(\Div_g\! \mathring{f}\right)- \frac{1}{2} K_g F\right) (\log K_{g_t} -1)
\end{equation*}
and so
\begin{equation*}
 \delta_f\mathcal{E}[g]= \frac{1}{2}\int_{M}F K_g \log K_g+\left(\Div_g\left(\Div_g\! \mathring{f}\right)-\frac{1}{2} \Delta_g F - K_g F\right) (\log K_{g_t} -1)\mu_g.
\end{equation*}
If $f$ is compactly supported
then integrating by parts yields
\begin{equation*}
 \delta_f\mathcal{E} [g]= -\frac{1}{4}\int_{M} F \left(\Delta_g \log K_g +2  K_g   \right)  -2\langle f, \tensor[^g]{\mathring{\nabla}}{}^2 \log K_g\rangle_g \mu_g,
\end{equation*}
where $\tensor[^g]{\mathring{\nabla}}{}^2$ denotes the trace-free Hessian and $\langle a,b\rangle_g$ the natural bilinear pairing on elements $a,b\in \Gamma(S^2(T^*M))$ obtained via $g$.

We will say that $g\in \mathcal{M}^+(M)$ is an $\mathcal{E}$-critical metric if $\mathcal{E}$ is stationary at $g$ with respect to compactly supported conformal deformations.
Hence,  $g$ is an $\mathcal{E}$-critical metric if and only if the Gauss curvature $K_g$ of $g$ satisfies
\begin{equation*}
 \Delta_g \log K_g = -2 K_g.
\end{equation*}

As $\mathcal{E}$ is computed purely in terms of geometric quantities it is manifestly diffeomorphism invariant, that is
if $\phi: M \to M$ is a diffeomorphism we have
\begin{equation*}
 \mathcal{E}[\phi^* g] =\mathcal{E}[g].
\end{equation*}
By Noether's principle this invariance leads to a conservation law for $\mathcal{E}$-critical metrics.  Indeed, let $X$ be a compactly supported vector field on $M$ and $\phi_t$ the flow of $X$. We have that
\begin{equation*}
 \phi_t^* g=g+ t L_X g +o(t).
\end{equation*}
Recall, the Lie derivative, $L_X g$ can be computed as
\begin{equation*}
 (L_X g)(Y,Z)=g( \gnabla_Y X,Z) + g(\gnabla_Z X, Y).
\end{equation*}
where $Y,Z \in \Gamma(TM)$. By the diffeomorphism invariance we have at an $\mathcal{E}$-critical metric that
\begin{align*}
 0 &= \delta_{L_X g} \mathcal{E} [g] \\ &= \frac{1}{2} \int_{M} \langle L_X g, \tensor[^g]{\mathring{\nabla}}{}^2 \log K_g\rangle_g \mu_g\\
   &=\int_{M} \Div_g\! \left( (\tensor[^g]{\mathring{\nabla}}{}^2 \log K_g )(X, \cdot)\right) - \langle X, \Div_g\! \tensor[^g]{\mathring{\nabla}}{}^2 \log K_g \rangle_g \mu_g\\
   &=-\int_{M} \langle X, \Div_g\! \tensor[^g]{\mathring{\nabla}}{}^2 \log K_g \rangle_g \mu_g
\end{align*}
where we used that $X$ has compact support and the divergence theorem.
As $X$ is arbitrary,
\begin{equation}\label{divfree}
 \Div_g\! \tensor[^g]{\mathring{\nabla}}{}^2 \log K_g =0.
\end{equation}
In other words, the quantity
\begin{equation*}
V_g=\tensor[^g]{\mathring{\nabla}}{}^2 \log K_g
\end{equation*}
 is a trace-free divergence free symmetric $2$-form, i.e.~a conservation law for the $\mathcal{E}$ functional.

\subsection{The entropy form}
Let $g$ be a smooth Riemannian metric and $\omega$ a smooth real-valued function on $M$. We note the following standard formula for the trace-free Hessian and the Laplacian operating on $u\in C^{\infty}(M)$
\begin{equation}\label{conformaltrafo}
\aligned
\Delta_{e^{2\omega}g}u&=e^{-2\omega}\Delta_gu,\\
^{e^{2\omega}g}\mathring{\nabla}^2&u=\tensor[^g]{\mathring{\nabla}}{}^2 u-\left(\d u \otimes \d \omega+\d \omega \otimes \d u-g(\gnabla u,\gnabla \omega)g\right).
\endaligned
\end{equation}
Also, the Gauss-curvature transforms under conformal rescaling as
$$
K_{e^{2\omega}g}=e^{-2\omega}\left(K_g-\Delta_g \omega \right).
$$
We let $\mathcal{R}^{\pm}_{C}$ denote the space of smooth positively (or negatively) curved Riemannian metrics on $M$ satisfying the \textit{generalized Ricci condition}
$$
\Delta_g\log|K_g|=CK_g
$$
for some real constant $C$. In particular, the $\mathcal{E}$-critical metrics are the elements of $\mathcal{R}^+_{-2}$. For $g \in \mathcal{R}^{\pm}_{C}$ and $\alpha \in \mathbb{R}$ let $g_{\alpha}=|K_g|^{2\alpha}g$, then $g_{\alpha}$ has Gauss-curvature
$$
K_{g_{\alpha}}=(1-C\alpha)|K_g|^{-2\alpha}K_g
$$
which, for $\alpha \neq \frac{1}{C}$, satisfies
\begin{equation*}
 \Delta_{g_{\alpha}} \log |K_{g_\alpha}| = \left(\frac{2\alpha -1}{\alpha-\frac{1}{C}}\right)K_{g_\alpha}=C_{\alpha}K_{g_{\alpha}}.
\end{equation*}
It follows that for $\alpha>\frac{1}{C}$ the map $\varphi_{\alpha}$ sending $g$ to $g_{\alpha}$ satisfies
$$
\varphi_{\alpha} : \mathcal{R}^{\pm}_{C} \to \mathcal{R}^{\mp}_{C_{\alpha}}
$$
whereas for $\alpha < \frac{1}{C}$
$$
\varphi_{\alpha} : \mathcal{R}^{\pm}_{C} \to \mathcal{R}^{\pm}_{C_{\alpha}}.
$$
Note that the choice $\alpha=\frac{1}{C}$ maps the elements of $\mathcal{R}^{\pm}_{C}$ to flat metrics and the choice $\alpha=\frac{1}{2}$ (assuming $C\neq 2$) maps the elements of $\mathcal{R}^{\pm}_{C}$ to metrics of non-zero constant Gauss curvature.

Suppose $g \in \mathcal{R}^{-}_{4}$, then $\hat{g}=g_{3/8}$ is an $\mathcal{E}$-critical metric with Gauss curvature
$$
K_{\hat{g}}=\frac{1}{2}|K_g|^{1/4}.
$$
It follows with \eqref{divfree} that the symmetric trace-free \textit{entropy form}
$$
\aligned
T&=V_{\hat{g}}={}^{\hat{g}}\mathring{\nabla}^2\log K_{\hat{g}}={}^{\hat{g}}\mathring{\nabla}^2\left(\frac{1}{4}\log |K_g|\right)\\&=-{}^{\hat{g}}{\mathring{\nabla}}^2u_g=-\tensor[^g]{\mathring{\nabla}}{}^2u_g-3\left((\d u_g)^2-\frac{1}{2}g(\gnabla u_g,\gnabla u_g)g\right)
\endaligned
$$
is divergence-free with respect to $\hat{g}$.

If $M$ is oriented, then we call the associated quadratic differential $P=T+i J T$ the \textit{entropy differential}. Since the condition on a symmetric $2$-form on $M$ to be trace-free and divergence free is conformally invariant, we obtain:
\begin{thm}
Let $(M,g)$ be a smooth oriented Riemannian $2$-mani\-fold with $K_g<0$ and $g$ satisfying the Ricci condition. Then the entropy differential
$P=T+i J T$ is holomorphic.
\end{thm}

\begin{rem}
Note that a metric of negative Gauss-curvature on a surface arising via a constant mean $2$ immersion into hyperbolic $3$-space $\mathbb{H}^3$ also satisfies the Ricci condition (recall that with our convention the `mean' curvature is the sum of the principal curvatures). Besides satisfying the Ricci condition, these so-called \textit{Bryant surfaces} share many properties with minimal surfaces in Euclidean $3$-space, the most important being that they possess a Weierstrass representation~\cite{bryant}.
\end{rem}

\subsection{The inverse problem}
Suppose we are given a Riemann surface $(M,J)$ and a holomorphic quadratic differential $P$ on $M$.  We ask whether we can locally find a $J$-compatible metric $g$ of negative Gauss curvature on $M$ which satisfies the Ricci condition and so that the entropy differential of $g$ is $P$.

Let $z: V \to \mathbb{C}$ be local holomorphic coordinates on $(M,J)$. It is easy to check that if the real-valued function $u$ solves Liouville's equation
\begin{equation}\label{liouville}
4\partial^2_{z\bar z}u=e^{-2u},
\end{equation}
then the metric
$
g=e^{2u}\d z \circ \d \bar{z}
$
satisfies the Ricci condition and $u_g=-\frac{1}{4}\log |K_g|=u$. Now a straightforward computation yields
$$
\aligned
\operatorname{Re}(P)=T&=-{}^{g_0}\mathring{\nabla}^2u-\d u^2 +\frac{1}{2}g_0({}^{g_0}\nabla u,{}^{g_0}\nabla u)g_0\\
&=-2\operatorname{Re}\left(\left(\partial^2_{zz}u+(\partial_z u)^2\right)\d z^2\right)
\endaligned
$$
where $g_0=\d z \circ \d \bar z$. Writing $P=\frac{\rho}{2}\,\d z^2$ for some holomorphic function $\rho$ on $V$, we are thus interested in the solutions $u$ of the system
$$
4\partial^2_{z\bar z}u=e^{-2u}, \quad \partial^2_{zz} u+(\partial_z u)^2=-\frac{\rho}{4}.
$$
\begin{lem}\label{cplxode}
Let $V\subset \mathbb{C}$ be a simply-connected domain and $\rho$ a holomorphic function on $V$.  We let $z$ be the usual complex coordinate on $\mathbb{C}$.  Then there exist holomorphic functions $w_1, w_2$ on $V$ solving the equation
\begin{equation} \label{linqeqn}
 \partial_{zz}^2 w+\frac{\rho}{4} w=0
\end{equation}
and with Wronskian satisfying
\begin{equation} \label{WronskianCond}
W(w_1,w_2)=w_1 \partial_z w_2-w_2 \partial_z w_1=\frac{1}{2}.
\end{equation}
If $\hat{w}_1, \hat{w}_2$ is another pair of holomorphic solutions to \eqref{linqeqn} satisfying \eqref{WronskianCond}, then there is a unique matrix $B \in \mathrm{SL}(2,\mathbb{C})$
so that $\hat{\wW}=B\wW$ where
\begin{equation*}
 \wW= \begin{pmatrix} w_1 \\ w_2 \end{pmatrix} \mbox{ and }\hat{\wW}=\begin{pmatrix} \hat{w}_1\\ \hat{w}_2 \end{pmatrix}.
\end{equation*}
\end{lem}
\begin{proof}
See for instance~\cite{Hille}.
\end{proof}

We now have the following:
\begin{prop} \label{SurfaceFromAandTProp}
 Let $V\subset \mathbb{C}$ be a simply-connected domain and suppose that $\rho$ is a holomorphic function on $V$. Then every real-valued function $u\in C^{\infty}(V)$ that satisfies the system
\begin{equation}\label{frobsys}
4 \partial^2_{z\bar{z}}u=e^{-2u}, \quad  \partial_{zz}^2 u+(\partial_z u)^2=-\frac{\rho}{4},
\end{equation}
is of the form
$$
 u_{\wW}=\log |\wW|^2=\log \left(|w_1|^2 +|w_2|^2\right),
$$
where
\begin{equation*}
 \wW= \begin{pmatrix} w_1 \\ w_2 \end{pmatrix}
\end{equation*}
and $w_1, w_2$ satisfy \eqref{linqeqn} and \eqref{WronskianCond}.  Hence, for each $\rho$ there is a three-dimensional space $\mathrm{SL}(2,\mathbb{C})/\mathrm{SU}(2)$ of solutions $u$.
\end{prop}
\begin{proof}
On $X=V\times \mathbb{R}\times \mathbb{C}$ with coordinates $(z,u,q)$ consider the rank $2$ subbundle $E\subset TX$ defined by the common kernel of the $1$-forms
$$
\varphi_1=\d u -q \d z -\bar q\d \bar z, \quad \varphi_2=\d q +\left(\frac{\rho}{4}+q^2\right)\d z -\frac{1}{4} e^{-2u}\d \bar z.
$$
Now
$$
\aligned
\d\varphi_1&=\d z\wedge \varphi_2+ \d\bar{z}\wedge \overline{\varphi_2},\\
\d\varphi_2&=-\frac{1}{2} e^{-2u} \d\bar{z}\wedge \varphi_1-2q \d z \wedge \varphi_2 ,\\
\endaligned
$$
hence $E$ is Frobenius integrable. Furthermore, the $1$-graph
$$
z \mapsto (z,u(z),\partial_z u(z))
$$ of a solution $u$ to \eqref{frobsys} is an integral manifold of $E$. Consequently, a solution $u$ to \eqref{frobsys} is uniquely determined by specifying $u$ and $\partial_z u$ at some point $z_0 \in V$. Simple computations show that for any $(z_0,u_0,q_0) \in X$ there exist holomorphic map $\wW: V \to \mathbb{C}^2$ satisfying (\ref{linqeqn}, \ref{WronskianCond}) so that $u_{\wW}=\log |\wW|^2$ solves \eqref{frobsys} and satisfies
$$
u(z_0)=u_0, \quad \partial_z u(z_0)=q_0.
$$
Clearly, if $\hat{\wW}=U\wW$ for $U \in \mathrm{SU}(2)$, then $u_{\hat{\wW}}=u_{\wW}$.
\end{proof}
\begin{cor}\label{MinSurfaceFromHillsEqnCor}
Let $V\subset \mathbb{C}$ be a simply-connected domain and suppose that $\rho$ is a holomorphic function on $V$. Let
\begin{equation*}
g_{\wW}=|\wW|^4 \d z \circ \d \bar z \quad \text{and} \quad A=\operatorname{Re}(\d z^2)
\end{equation*}
where
\begin{equation*}
 \wW= \begin{pmatrix} w_1 \\ w_2 \end{pmatrix}
\end{equation*}
and $w_1, w_2$ satisfy \eqref{linqeqn} and \eqref{WronskianCond}. Then there is a minimal immersion $\xX_\wW:V\to \Sigma\subset \mathbb R^3$ with geometric data $(V, g_{\wW}, A)$ and entropy differential $P=\frac{\rho}{2} \d z^2$.
\end{cor}
\begin{proof}
 This is an immediate consequence of Proposition \ref{SurfaceFromAandTProp} and the fundamental theorem of submanifold geometry.
\end{proof}

\section{Weierstrass Representation}\label{WeierstrassSec}
In this section we express the entropy differential $P$ in terms of the Weierstrass data of the minimal surface $\Sigma$ -- this allows us to compute $P$ more readily and to easily analyze its singular and asymptotic behavior.

\subsection{The Weierstrass Representation}
Recall, to an oriented minimal surface $\Sigma$ in $\mathbb R^3$ with parametrization $\xX_\Sigma:M \to \Sigma$ one can associate \emph{Weierstrass data} which encodes the surface and parametrization $\xX_\Sigma$ in complex analytic data. More precisely, the Weierstrass data associated to $\xX_\Sigma$ is the quadruple $(M, J, G, \eta)$ where $(M,J)$ is a Riemann surface, $G$ is a meromorphic function on $(M,J)$ and $\eta$ a holomorphic one form on $(M,J)$.  The data is determined as follows:
\begin{enumerate}
 \item $J$ is the almost-complex structure induced by $\xX_\Sigma$;
\item $G=S\circ \nN$ where $\nN$ is the Gauss map and
$$S:\partial B_1 \backslash (0,0,-1)\to \mathbb{C}$$
 is stereographic projection;
\item $\xX_\Sigma^* \d x_3= \re \eta$.
\end{enumerate}
The Weierstrass data allows one to reconstruct $\xX_\Sigma$ by the means of the \emph{Weierstrass representation}:
\begin{equation} \label{WeierstrassRep}
 \xX_\Sigma(p)-\xX_\Sigma(p_0)= \re \int^p_{p_0} \left( \frac{1}{2} (G^{-1}-G), \frac{i}{2 }( G^{-1}+G) , 1\right) \eta.
\end{equation}
Conversely, given any quadruple $(M, J, G, \eta)$ we may use \eqref{WeierstrassRep} to construct a parametrization $\xX_\Sigma$ of a branched minimal surface $\Sigma$ provided:
\begin{enumerate}
 \item \label{WRCond1}Both $G \eta$ and $G^{-1} \eta$ are holomorphic;
 \item \label{WRCond2}For any 1-cycle $\gamma$  in $M$:
  \begin{equation*}
   \int_\gamma \left( \frac{1}{2} (G^{-1}-G), \frac{i}{2 }( G^{-1}+G) , 1\right) \eta\in i \mathbb R^3.
  \end{equation*}
\end{enumerate}
Condition \eqref{WRCond2} is known as the period condition.
\begin{rem}
 The parametrizing map $\xX_\Sigma$ is an immersion if and only if $G\eta, G^{-1} \eta$, and $\eta$ do not all simultaneously vanish at any point of $M$.
\end{rem}

It is convenient to choose a local complex coordinate patch $(V, z)$ on $M$ and to write $\eta=h\d z$ and $G=G(z)$.  We write $f'$ for $\partial_z f$ for any function $f\in C^{1}(V, \mathbb{C})$.
Standard computations (see for instance~\cite{KarcherNotes}) give:
the metric as
\begin{equation*}
 g=\xX_\Sigma^* g_E= \frac{1}{4} (|G|+|G|^{-1})^2 \eta \circ \overline{\eta}=\frac{|h|^2}{4} (|G|+|G|^{-1})^2 \d z \circ \d\overline{z},
\end{equation*}
the Hopf differential as
\begin{equation*}
Q = -\frac{1}{G}\d G\circ \eta = -\frac{h G'}{G} \d z^2,
\end{equation*}
and the Gauss curvature
\begin{equation*}
  K_{g} =- \frac{16  |G G'|^2}{ |h| (1+|G|^{2})^4}.
\end{equation*}
Hence,
\begin{equation*}
 u_g =-\log 2 -\frac{1}{4} \log |h^{-1} G G'|^2 + \log (1+|G|^2).
\end{equation*}

\subsection{Computing $P$ in terms of Weierstrass data}
We now compute the entropy differential $P$ in terms of the Weierstrass data.
\begin{prop} \label{WeierstrassCompProp}
 Let $\Sigma$ be an oriented minimal surface in $\mathbb R^3$ with Weierstrass data $(M, J, G, \eta)$. If $(U, z)$ is a coordinate chart of $M$ on which  $K_g<0$ and we write $\eta=h\d z$, $G=G(z)$, then $P=\frac{\rho}{2}\d z^2$ with
\begin{align*}
 \rho&= \left(\frac{G'''}{G'}  +\frac{G''}{2G}-\frac{3(G')^2}{4G^2} -\frac{7(G'')^2}{4(G')^2}+\frac{G'' h'}{2G' h} -\frac{G' h'}{2Gh}  -\frac{h''}{h}+\frac{5(h')^2}{4h^2} \right).
\end{align*}
If $Q=\d z^2$, then
\begin{align*}
 P &= \left(\left(\frac{G''}{G'}\right)' -\frac{1}{2} \left( \frac{G''}{G'}\right)^2\right) \d z^2\\
   &= \lbrace G,z\rbrace \d z^2,
\end{align*}
where $\lbrace G, z \rbrace$ is the Schwarzian derivative of $G$.
\end{prop}
\begin{rem}
 The Schwarzian derivative of $G$ has also been studied from a different perspective by Duren, Chuaqui and Osgood~\cite{DurenEtAl} (see also~\cite{MR1656822} for a coordinate free definition of the Schwarzian derivative).
\end{rem}
\begin{proof}
If  $K_g<0$ on $V$, then $ \frac{h G'}{G}  \neq 0$ on $V$.  Hence, if $V$ is simply connected there is global square root of $-\frac{hG'}{G}$. Indeed,  there is a function $w$ on $V$ so that
\begin{equation*}
 \d w= \sqrt{-\frac{h G'}{G}} \d z
\end{equation*}
and so
\begin{equation*}
 Q=\d w^2.
\end{equation*}
The exact one-form $\d w$ is well-defined up to multiplication by $\pm 1$.
In particular, we have that the entropy differential is given by
\begin{equation*}
 P= -2 \left( \partial^2_w u_g +(\partial_w u_g)^2\right) \d w^2.
\end{equation*}

In order to express $P$ in terms of the Weierstrass data we note that:
\begin{equation*}
 \partial_w = \sqrt{ -\frac{G}{h G'}} \partial_z
\end{equation*}
and so
\begin{equation*}
 \partial_w^2 = -\frac{G}{h G'} \frac{\partial^2}{\partial z^2}-\frac{G}{2 h G' } \left(\frac{G'}{G}-\frac{G'' }{G' } -\frac{h'}{h}\right) \frac{\partial}{\partial z}.
\end{equation*}
Hence,
\begin{align*}
 \sqrt{- \frac{h G'}{G}}\partial_w u_g& =-\frac{1}{4}\frac{ (h^{-1} G G')'}{h^{-1} G G'} +\frac{ G' \bar{G}}{1+|G|^2}\\
                              &= \frac{1}{4}\left(\frac{h'}{h} -\frac{G'}{G} -\frac{G''}{G'} \right) +\frac{ G' \bar{G}}{1+|G|^2}
\end{align*}
and
\begin{align*}
 -\frac{h G'}{G}\partial^2_w u_g=& \frac{1}{4} \left( \frac{h''}{h}-\left(\frac{h'}{h}\right)^2-\frac{G''}{G}+\left(\frac{G'}{G}\right)^2-\frac{G'''}{G'}+\left(\frac{G''}{G'}\right)^2\right)+\\&+\frac{G'' \bar{G}}{1+|G|^2}-\left(\frac{G' G}{1+|G|^2}\right)^2+\frac{1}{8} \left( \frac{G''}{G'}\right)^2-\frac{1}{8}\left( \frac{h'}{h}-\frac{G'}{G}\right)^2-\\& -\frac{1}{2}\left( \frac{G''}{G'}+\frac{h'}{h}-\frac{G'}{G} \right) \frac{ G' \bar{G}}{1+|G|^2}.
\end{align*}
We note that both these expressions are independent of replacing $w$ by $-w$ and so hold even if $V$ is not simply-connected.
Combining the above we determine that
$
P=\frac{\rho}{2}\,\d z^2
$
with
\begin{align*}
 \rho&= \left(\frac{G'''}{G'}  +\frac{G''}{2G}-\frac{3(G')^2}{4G^2} -\frac{7(G'')^2}{4(G')^2}+\frac{G'' h'}{2G' h} -\frac{G' h'}{2Gh}  -\frac{h''}{h}+\frac{5(h')^2}{4h^2} \right).
\end{align*}
as claimed.
If $Q=\d z^2$, then
\begin{equation*}
 h=-\frac{G}{G'}
\end{equation*}
and so
\begin{equation*}
 \frac{h'}{h}=\frac{G'}{G}-\frac{ G''}{G'}
\end{equation*}
 and
\begin{equation*}
 \frac{h''}{h}=-\frac{ G''}{G}+2\frac{ (G'')^2}{(G')^2}.
\end{equation*}
Plugging these into the formula for $P$ gives
\begin{align*}
  P &= \left( \left(\frac{G''}{G'}\right)'-\frac{1}{2}\left(\frac{G''}{G'}\right)^2\right) \d z^2.
\end{align*}
\end{proof}

As an application of the previous computation, we determine the behavior  of the entropy differential at umbilic points of $\Sigma$:
\begin{cor}\label{DoublePoleCor}
 If $\Sigma$ is a minimal surface in $\mathbb R^3$ and $p\in \Sigma$ an isolated umbilic point, then $P$, the entropy differential of $\Sigma$, has a double pole at $p$.  Indeed, there is a complex coordinate $z$ around $p$ satisfying $z(p)=0$ and so
\begin{equation*}
 P= -\left( \frac{3n^2+4n}{8} \right) \frac{\d z^2}{z^2} +O(1),
\end{equation*}
where $n$ is the order of vanishing of the Hopf differential $Q$ at $p$.
\end{cor}
\begin{proof}
By rotating $\Sigma$ in $\mathbb R^3$, we may assume that $\nN(p)=\eE_1$ where $(\eE_1,\eE_2,\eE_3)$ denotes the standard basis of $\mathbb{R}^3$. Hence, there is $p$-neighbor\-hood $V$ with a $p$-centered complex coordinate $z$, together with Weierstrass data $(V, J, G, \eta)$ parametrizing $\Sigma$ near $p$ which satisfies $\eta=\d z$ and $ G(z)=1+o(1)$.
In fact, there is an $a,b\in \mathbb{C}$ with $a\neq 0$ so that
\begin{equation*}
 G(z)=1+a z^{n+1} +b z^{n+2} +O(z^{n+2}),
\end{equation*}
because the umbilic point is isolated.
Indeed,
\begin{equation*}
 Q=-\frac{h G'}{G} \d z^2= -a (n+1) z^{n} \d z^2 +O(z^n)
\end{equation*}
and $n$ is the order of vanishing of $Q$ at $p$.

We let $V^*=V\backslash\set{p}$ and apply Proposition \ref{WeierstrassCompProp} to compute that
\begin{equation*}
 P= \frac{1}{2}\left( -\left(\frac{3}{4} n^2+n\right)z^{-2}-\frac{3}{2} \frac{n(n+2)}{n+1} \frac{b}{a} z^{-1} \right) \d z^2 +O(1).
\end{equation*}
However, by changing coordinates to $z\to z+c z^2$ for an appropriate choice of $c$ we obtain $P$ in the desired form.
\end{proof}
We may also use Proposition \eqref{WeierstrassCompProp} to compute the entropy differential at branch points.
\begin{cor} \label{BranchPtCor}
 Suppose that $(M, J)$ is a Riemann surface and $\xX:M\to \Sigma\subset \mathbb R^3$ is a non-flat branched minimal immersion. Let $p\in M$ be a branch point of $M$ of order $n$ and index $k$.
\begin{enumerate}
 \item If $n-k+1\neq 0$, then the entropy differential, $P$, has a double pole at $p$ and there is a complex coordinate patch $(V,z)$ about $p$ with $z(p)=0$ so that
\begin{equation*}
 P=\left(\frac{(n+k+1)^2-4 k^2}{8 } \right)\frac{\d z^2}{z^2}+O(1);
\end{equation*}
\item If $n-k+1=0$, then $P$ has at most a simple pole at $p$.
\end{enumerate}
\end{cor}
\begin{proof}
 We may pick a complex coordinate patch $(V,z)$ about $p$ so that $z(p)=0$ and on $V^*=V\backslash \set{p}$ the parameterization $\xX$ is a smooth immersion. Let $(z(V^*), J, G, \eta)$ be the Weierstrass data of this immersion where here $J$ is the usual complex structure.  As $\xX$ has an order $n\geq 1$ branch point with index $k$ at $z(p)=0$, up to an ambient rotation of $\mathbb R^3$ and a re-parameterization the data has the form
\begin{equation*}
 G(z)= z^k
\end{equation*}
for $k\geq 1$ and
\begin{equation*}
 \eta=\left(a z^{n+k}+bz^{n+k+1}\right) \d z+O(z^{n+k+2})
\end{equation*}
where $a\neq 0$.
Computing gives
\begin{align*}
 P&=\left(\frac{(n+k+1)^2-4 k^2}{8 z^2}+\frac{\frac{b}{a}(n+k-1)}{2z}\right) \d z^2+O(1)\\
 &=\left(\frac{(n-k+1)(n+3k+1)}{8 z^2}+\frac{\frac{b}{a}(n+k-1)}{2z}\right) \d z^2+O(1).\\
\end{align*}
The corollary follows by noting that if $n-k+1=0$, then $P$ has at most a simple pole at $p$ as claimed. If $n-k+1\neq 0$, then $P$ has a double pole and may be put in the claimed form by replacing $z$ by $z+cz^2$ for an appropriate choice of $c$.
\end{proof}
\begin{rem}
 We do not distinguish between true and false branch points.  However, any false branch point of a smooth minimal surface at a point with non-vanishing curvature must have order of vanishing $n$ and index $n+1$.
\end{rem}

\subsection{Hill's equation and the (spinor) Weierstrass representation}
We conclude by relating the solutions $w_1,w_2$ from Proposition \ref{SurfaceFromAandTProp} to the Weierstrass data. We observe a connection with the spinorial Weierstrass representation of~\cite{Kusner} but do not explore this in depth.
\begin{prop}\label{WeierstrassHillEqnProp}
Fix a simply-connected domain $V\subset \mathbb{C}$.
Suppose $(V, J_{std}, G, h\,\d z)$ is the Weierstrass data of a minimal immersion with Hopf differential $Q=\d z^2$ and entropy differential $P=\frac{\rho}{2}dz^2$, then
\begin{align*}
 w_1(z)&=\frac{\sqrt{2}}{2}  \sqrt{-G^{-1}(z) h(z)}\\
 w_2(z)&=\frac{\sqrt{2}}{2}  \sqrt{-G(z) h(z)}
\end{align*}
are single-valued and satisfy \eqref{linqeqn}.  Furthermore, $w_1$ and $w_2$ satisfy \eqref{WronskianCond} provided the branches of the square-root are chosen so $\frac{w_2}{w_1}=G$.
\end{prop}
\begin{proof}
As $G$ and $h\, \d z$ is the Weierstrass data of a minimal immersion,  $G h$ or $G^{-1} h$ do not have a pole on $V$.  Moreover, if either function vanished at a point $z_0$, then $h(z_0)=0$. As $Q=\d z^2$,  $-\frac{G'}{G} h=1$. Because $G$ is meromorphic,  $h$ has at most a simple zero at $z_0$ and so $G$ has either a simple pole or a simple zero at $z_0$. Hence, at $z_0$ either $G h\neq 0$  and $G^{-1} h$ has a double zero or $G^{-1} h\neq 0$  and $G h$ has a double zero.  Taken together this implies that $w_1$ and $w_2$ are single-valued.

A straightforward computation gives that $w_1, w_2$ satisfy the Wronskian condition~\eqref{WronskianCond}.  Differentiating \eqref{WronskianCond} once, gives that
$$\frac{w_1''}{w_1}=\frac{w_2''}{w_2}=-\frac{\hat{\rho}}{4}$$
for a meromorphic function $\hat{\rho}$.
It is a classical fact -- see for instance~\cite{Hille} -- that if $\hat{w}_1, \hat{w}_2$ solve $w''+\frac{\hat{\rho}}{4} w=0$, then $\hat{G}=\frac{\hat{w}_1}{\hat{w}_2}$ satisfies $\{\hat{G}, z\}=\frac{\hat{\rho}}{2}$. As $\frac{w_2}{w_1}= G$, this implies that $\hat{\rho}=2 \{ G, z \}=\rho$ and so $w_1, w_2$ satisfy \eqref{linqeqn}.
\end{proof}
\begin{cor} \label{WeierstrassFromHillsEqnCor}
Let $V\subset \mathbb{C}$ be a fixed simply-connected domain.
 If $\rho$ is a holomorphic function on $V$ and
$\wW=(w_1, w_2)^\top$ satisfies \eqref{linqeqn} and \eqref{WronskianCond}, then the minimal immersion $\xX_\wW$ of Corollary \ref{MinSurfaceFromHillsEqnCor} may be be chosen to have Weierstrass data  $\left(V, J_{std},G, \eta\right)$, where
\begin{equation*}
G= \frac{w_2}{w_1}\quad \text{and}\quad \eta=-2w_1 w_2 \d z.
\end{equation*}
\end{cor}
\begin{rem}
 If we let $s_i=w_i \sqrt{\d z} $ be holomorphic spinors, then the $s_i$ are (up to choices of normalization) the spinor Weierstrass data of~\cite{Kusner}.
\end{rem}
\begin{proof}
Set ${G}=\frac{{w}_2}{{w}_1}$ and  $\eta=-2{w}_1{w}_2 \d z$ and let $\xX_{{\wW}}$ be the  minimal immersion corresponding to this data.
As $w_1=\frac{\sqrt{2}}{2}  \sqrt{-G^{-1}(z) h(z)}$ and $w_2=\frac{\sqrt{2}}{2}  \sqrt{-G(z) h(z)}$, Proposition \ref{WeierstrassHillEqnProp} implies that the entropy differential of $\xX_\wW$ is $\frac{\rho}{2} dz^2$.
A direct computation and \eqref{WronskianCond} imply that the Hopf differential of $\xX_{\wW}$ is $dz^2$.
Finally,
$$
\xX_{\wW}^*g_E=\frac{1}{4} |h|^2\left( |G|+|G^{-1}|\right)^2 dz \otimes d\bar{z}=|\wW|^4 \d z \otimes \d\bar{z} =g_{\wW}.
$$
Hence, $\xX_{\wW}$ satisfies the conclusions of Corollary \ref{MinSurfaceFromHillsEqnCor} which verifies the claim.
\end{proof}

\section{Characterization of Minimal Surfaces in Terms of the Entropy Differential}\label{CharSec}
In this section we characterize a number of classical minimal surfaces in terms of the entropy form. In particular, we show that the entropy form vanishes if and only if the surface is contained in Enneper's surface.  The catenoid and helicoid are also characterized in terms of a simple relationship between the entropy form and the second fundamental form.

\subsection{Deformed Catenoids and Helicoids}
In order to get a complete characterization we must introduce two one-parameter families of surfaces, $\mathcal{C}$ and $\mathcal{H}$, which we call, respectively, deformed catenoids and deformed helicoids.
Specifically, $\mathcal{C}$ is the family of surfaces $C_t$ with Weierstrass data
\begin{equation*}
\left(\mathbb{C}, J, \frac{t-e^z}{1-t e^z}, \frac{1}{1-t^2} (1-t e^{-z})(1-t e^{z})\d z\right).
\end{equation*}
Similarly, $\mathcal{H}$ is the family of surfaces $H_t$ with Weierstrass data
\begin{equation*}
\left(\mathbb{C}, J, \frac{t- e^{z}}{1- t e^{z}}, \frac{-i}{1-t^2} (1-t e^{-z})(1-t e^{z})\d z\right).
\end{equation*}
In both cases, $z$ is the usual coordinate on $\mathbb{C}$, $J$ the usual complex structure and $t\in (-1,1)$.  In particular,  $C_0$ is the vertical catenoid and $H_0$ is the vertical helicoid.  Computing as in the preceding section we obtain that for surfaces in $\mathcal{C}$
\begin{equation*}
 P=-\frac{1}{2} \d z^2=\frac{1}{2} Q,
\end{equation*}
and for surfaces in $\mathcal{H}$
\begin{equation*}
 P=-\frac{1}{2} \d z^2=\frac{i}{2} Q.
\end{equation*}
We remark that $\mathcal{C}$ and $\mathcal{H}$ are obtained from $C_0$ and from $H_0$ by applying the one parameter family of M\"{o}bius transforms
$$
B_t:z\mapsto \frac{t+z}{1-t z}
$$
to the Gauss maps of $C_0$ and $H_0$.

Writing $z=x+iy$ and integrating \eqref{WeierstrassRep} gives the parameterizations of $C_t \in \mathcal{C}$:
\begin{align*}
 \FF_t^{\mathcal{C}}(x,y)&=  \FF_0^{\mathcal{C}}(x,y)+\frac{2 t}{1-t^2} \left( 0,- y+t \cosh x \sin y   , t x-\sinh x \cos y\right) ;\\
\FF_0^{\mathcal{C}}(x,y)&=\left(\cosh x \cos y , \cosh x \sin y, x\right).
\end{align*}
Here $\FF_0^{\mathcal{C}}$ is a parameterization of (an infinite cover of) the catenoid.  By inspection,  $C_t$ has $\Pi_2=\set{x_2=0}$ and $\Pi_3=\set{x_3=0}$ as planes of reflectional symmetry. Moreover,
$$\FF_t^{\mathcal{C}}(x,y+2\pi)=\FF_t^{\mathcal{C}}(x,y)-\frac{4\pi t}{1-t^2} \eE_2$$
 and so $C_t$ is singly-periodic.
When $t \neq 0$, it is straightforward to see that $C_t$ is not embedded. Suppose $E_\theta$ is the rotation of the upper half of $C_0$ by $\theta$ around the $x_2$-axis. One verifies that $C_t$ is close to the union of translates of $E_\theta$ and of $E_{\pi-\theta}$ where here $\theta=\tan^{-1}\left( \frac{2t}{1-t^2}\right)$.

Similarly, elements of $\mathcal{H}$ are parametrized by
\begin{align*}
 \FF_t^{\mathcal{H}}(x,y)&=  \FF_0^{\mathcal{H}}(x,y)+\frac{2t}{1-t^2} \left(0,  x+t\sinh x \cos y, t y- \cosh x \sin y\right);\\
\FF_0^{\mathcal{H}}(x,y)&= (\sinh x \sin y,- \sinh x \cos y, y).
\end{align*}
For $t=0$ this is a parametrization of the helicoid. Note that the image of $\set{x=0}$ is a the $x_3$-axis while the the image of $\set{y=n \pi}$ for $n$ an integer are the set of parallel lines $\set{x_1=0, x_3=\frac{1+t^2}{1-t^2} n \pi}$ contained in the $\set{x_1=0}$ plane.  Moreover,
$$\FF_t^{\mathcal{H}}(x,y+2\pi)= \FF_t^{\mathcal{H}}(x,y)+2\pi \frac{1-t^2}{1+t^2} \eE_3$$
 so $H_t$ is singly-periodic.
 For $t\neq 0$, $H_t$ is not embedded. However, if we denote by $H^\pm_t$ the two components of $H_t \backslash \set{x_1=x_2=0}$, then each $H^\pm_t$ is embedded.  In fact, each is a multi-valued graphs over the plane $\Pi_\theta$ which contains the $x_2$-axis and makes an angle $\theta=\tan^{-1}\left( \frac{2t}{1-t^2}\right)$ with the plane $\Pi_3=\set{x_3=0}$.  In particular, rotating $H^\pm_t$ by $\theta$ around the $x_2$-axis gives a surface that looks (roughly) like a sheered copy of $H_0^\pm$.

\subsection{Characterization of minimal surfaces in terms of $P$ and $Q$}
We now characterize surfaces in terms of simple relationships between $P$ and $Q$. In light of Proposition \ref{SurfaceFromAandTProp}, we expect there to be a three-parameter family of surfaces for any fixed of $P$ and $Q$.  However, in simple settings two of these parameters correspond to re-parameterizations.
\begin{thm}\label{MainCharacterizationThm}
 Let $\Sigma$ be a smooth oriented non-flat minimal surface in $\mathbb R^3$ with Hopf differential $Q$ and entropy differential $P$. We have:
\begin{enumerate}
 \item \label{EnneperCase} If $P\equiv 0$, then up to a rigid motion and homothety, $\Sigma$ is contained in Enneper's surface;
 \item \label{CatCase} If $\lambda \neq 0$ and  $P\equiv \lambda Q$, then,  up to a rigid motion and homothety, $\Sigma$ is contained in a surface $C\in\mathcal{C}$.  If $\Sigma$ is properly embedded, then it is the catenoid;
\item \label{HelCase} If $\lambda \neq 0$ and $P \equiv  i Q$, then, up to a rigid motion and homothety, $\Sigma$ is contained  in a surface $H\in \mathcal{H}$.  If $\Sigma$ is properly embedded, then it is the helicoid.
\end{enumerate}
\end{thm}
\begin{rem}
If $\Sigma$ is an oriented minimal surface with Hopf differential $Q$ and entropy differential $P$, then for any $\lambda>0$ the rescaling scaling of $\lambda \Sigma$ has Hopf differential $\lambda Q$ and entropy differential $P$.  Reversing the orientation of $\Sigma$ changes $Q$ to $-Q$ but leaves $P$ unchanged.
\end{rem}
\begin{proof}
After possibly rescaling $\Sigma$ and reversing the orientation, we may assume that $P=-\frac{1}{2} \alpha^2 \d z^2$ where $\alpha=0$ in Case \eqref{EnneperCase}, $\alpha^2=1$ in Case \eqref{CatCase} and $\alpha^2=i$ in Case \eqref{HelCase}.
As $\Sigma$ is smooth and non-flat, the second fundamental form has no singularities and $P$ can only have isolated singularities. Hence, by Corollary \ref{DoublePoleCor}, in all cases $P$ has no singularities and $Q$ has no zeros on $\Sigma$.
Hence, for any point $p\in \Sigma$ there is a simply connected neighborhood $V$ of $p$ and complex coordinate $z:V \to \mathbb{C}$ so that the Hopf differential satisfies $Q=-\d z^2$. That is, $P=\frac{\alpha^2}{2} dz^2$.
By Corollary \ref{WeierstrassFromHillsEqnCor}, in order to recover the surface it is enough to understand the holomorphic solutions on  $z(V)$ to the Hill's equation:
\begin{equation}\label{HillsEqnConst}
 \frac{\partial^2}{\partial z^2} w-\frac{\alpha^2}{4}w =0.
\end{equation}
Clearly, this equation makes sense on all of $\mathbb{C}$ (with $z$ as the usual coordinate) and analytic continuation implies that all solutions are obtained by restricting global solutions to $z(V)$.  Let $\wW(z)=(w_1(z), w_2(z))^\top$ be a pair of solution to the Hill's equation with Wronskian $W(w_1,w_2)=\frac{1}{2}$.

We note there are two natural actions on the space of solutions.  The first is the natural action of $\mathrm{SL}(2,\mathbb{C})$ of Proposition \ref{SurfaceFromAandTProp} which is transitive.  The second is an action of $\mathbb{C}$ that arises from the translation invariance of \eqref{HillsEqnConst}.  Specifically, let $\mathbb{C}$ act on $\wW$ by $\tau\mapsto \wW(z+\tau)$. The translation invariance of \eqref{HillsEqnConst} and of the Wronskian condition implies that this is a well defined action.  By Proposition \ref{SurfaceFromAandTProp}, the action of $\mathrm{SU}(2)\subset \mathrm{SL}(2,\mathbb{C})$ does not change the geometry of the surface. Likewise, the action of $\mathbb{C}$ amounts to a change of coordinates and also does not change the geometry.  Our goal is to determine all geometrically distinct solutions.

First, note that the Gram-Schmidt procedure implies that any matrix $B\in \mathrm{SL}(2,\mathbb{C})$ may be factored as
\begin{equation*}
 B=U L
\end{equation*}
where $U\in \mathrm{SU}(2)$ and $L\in \mathrm{SL}(2,\mathbb{C})$ is lower triangular with positive entries on the diagonal and $\det L=1$.  This is sometimes called the $QR$ (or in this case $QL$) factorization.
We write any such $L$ as
\begin{equation*}
L=\begin{bmatrix} \mu & 0 \\ \nu & \mu^{-1} \end{bmatrix}
\end{equation*}
where $\mu>0$ and $\nu \in \mathbb{C}$.
We now treat the case $\alpha=0$ and $\alpha \neq 0$ separately.
{\bf Case \eqref{EnneperCase}:}
By inspection a pair of solutions to \eqref{HillsEqnConst} with $\alpha=0$ and satisfying the Wronskian condition are
\begin{equation*}
 w_1(z)=1 \mbox{ and } w_2(z)=\frac{1}{2} z.
\end{equation*}
Hence, by the $QL$ factorization, the functions
\begin{equation*}
 w_1(z)=\mu \mbox{ and } w_2(z)=\nu+\frac{1}{2}\mu^{-1} z,
\end{equation*}
with $\mu>0$ and $\nu\in \mathbb{C}$ give all geometrically distinct solutions to \eqref{HillsEqnConst}.
Applying the translation action with  $\tau=-2\mu \nu $
gives all geometrically distinct solutions in the form
\begin{equation*}
 w_1(z)=\mu \mbox{ and } w_2(z)= \frac{1}{2}\mu^{-1} z.
\end{equation*}
By Corollary \ref{WeierstrassFromHillsEqnCor} the Gauss map the associated minimal surfaces maybe chosen so
\begin{equation*}
G(z)=\frac{w_2}{w_1}=\frac{ z}{2\mu^2}.
\end{equation*}
Moreover, as $Q=-\d z^2$ the height differential is $\eta=z \d z$. This is precisely the Weierstrass data of a rescaling of Enneper's surface proving the claim in this case.

{\bf Case \eqref{CatCase} and \eqref{HelCase}:}
As $\alpha\neq 0$, a pair of solutions to \eqref{HillsEqnConst} that satisfy the Wronskian condition are
\begin{equation*}
 w_1(z)=\frac{1}{\sqrt{2}\alpha } e^{-\frac{\alpha}{2} z} \mbox{ and } w_2(z)=\frac{1}{\sqrt{2}\alpha} e^{\frac{\alpha}{2} z}.
\end{equation*}
Hence, by the $QL$ factorization, we may write all geometrically distinct solutions  to \eqref{HillsEqnConst} in the form
\begin{equation*}
 w_1(z)= \frac{\mu}{\sqrt{2}\alpha } e^{-\frac{\alpha}{2} z} \mbox{ and } w_2(z)= \frac{1}{\sqrt{2}\alpha }\left(\nu e^{-\frac{\alpha}{2} z}  +\mu^{-1} e^{\frac{\alpha}{2} z}\right)
\end{equation*}
with $\mu>0$ and $\nu\in \mathbb{C}$.
The translation action
allows us to express all geometrically distinct solutions as
\begin{equation*}
 w_1(z)=i\frac{e^{-i \theta/2}}{\sqrt{2}\alpha}  e^{-\frac{\alpha}{2} z} \mbox{ and } w_2(z)=i\frac{e^{i \theta/2}}{\sqrt{2}\alpha} \left( \gamma e^{-\frac{\alpha}{2} z}  - e^{\frac{\alpha}{2} z}\right) .
\end{equation*}
where $\gamma\geq 0$ and $\theta\in [0,2\pi)$.  Indeed, either $\nu=0$ and we take $\gamma=\theta=0$ or $\nu\neq 0$ and we write $\nu=\gamma \mu e^{i\theta}$. In both cases, we act by $\tau=\frac{1}{\alpha} i(\theta-{\pi}) +\frac{2}{\alpha} \ln \mu$.
Let $\phi\in (-\pi/4,\pi/4)$ satisfy
\begin{equation*}
 \tan 2 \phi = \gamma.
\end{equation*}
The matrix
\begin{equation*}
 \begin{pmatrix} \cos \phi & -\sin \phi \\ \sin \phi & \cos \phi \end{pmatrix} \begin{pmatrix} -i e^{i\theta/2} & 0 \\ 0 & -ie^{-i\theta/2} \end{pmatrix}
\end{equation*}
is the product of two elements of $\mathrm{SU}(2)$ and so is in $ \mathrm{SU}(2)$.  Acting by this matrix gives that all geometrically distinct solutions can be put in the form
\begin{equation*}
 w_1(z)=\frac{\cos \phi }{\sqrt{2}\alpha \cos 2 \phi}  e^{-\frac{\alpha}{2} z}-\frac{\sin \phi }{\sqrt{2}\alpha }  e^{\frac{\alpha}{2} z}
\end{equation*}and
\begin{equation*} w_2(z)=\frac{\sin \phi}{\sqrt{2}\alpha \cos 2 \phi} e^{-\frac{\alpha}{2} z}  - \frac{\cos \phi }{\sqrt{2}\alpha }e^{\frac{\alpha}{2} z},
\end{equation*}
where $\phi\in (-\pi/4,\pi/4)$. By applying the translation action with $\tau= -\frac{1}{\alpha} \ln \cos 2\phi$, all  geometrically distinct solutions can be put in the simplified form
\begin{equation*}
 w_1(z)=\frac{\cos \phi e^{-\frac{\alpha}{2} z}-\sin \phi   e^{\frac{\alpha}{2} z} }{\sqrt{2 \cos 2 \phi}\alpha}  \mbox{ and } w_2(z)=\frac{\sin \phi e^{-\frac{\alpha}{2} z}  - \cos \phi e^{\frac{\alpha}{2} z} }{\sqrt{2 \cos 2 \phi}\alpha}.
\end{equation*}
By Corollary \ref{WeierstrassFromHillsEqnCor} the Gauss map the associated minimal surfaces may be chosen so
\begin{equation*}
G(z)=\frac{w_2}{w_1}= \frac{ \tan\phi - e^{\alpha z}}{1- \tan \phi e^{\alpha z}}.
\end{equation*}
Set $t= \tan \phi$.
If $\alpha^2=1$, then we may take $\alpha=1$ and as $Q=-\d z^2$ we see that $\eta = \frac{1}{1-t^2} (1-t e^{-z})(1-t e^{z})\d z$ which together with $G(z)$ is precisely the data of a deformed catenoid.
If $\alpha^2=i$, then we write $\zeta=\alpha z$.  In this case $Q=-\frac{1}{\alpha^2} \d\zeta^2=i \d\zeta^2$ and so
$\eta=- \frac{i}{1-t^2} (1-t e^{-\zeta})(1-t e^{\zeta})\d\zeta$ which together with $G(\zeta)$ is precisely the data of a deformed helicoid.
\end{proof}

\section{Curvature Estimate for Embedded Minimal Surfaces in Terms of $T$}
An interesting problem is to make the characterizations of Theorem \ref{MainCharacterizationThm} effective.  For instance, to show that a minimal surface with ``small'' entropy form must be close to a rescaling of a piece of Enneper's surface.  A major challenge is to determine an appropriate notion of smallness for the entropy form -- something made more difficult by the need to account for the possible singularities of $T$.  We propose a certain weighted $L^p$ norm of the entropy form as a natural way to measure this smallness and as an application give a novel curvature estimate for embedded minimal surfaces.

Before introducing the norm we note the following consequence of Corollary \ref{DoublePoleCor}.
\begin{lem} \label{weightedEntLem}
 Let $\Sigma$ be a smooth minimal surface with metric $g$ and entropy form $T$. We define,  $\hat{T}$,  the \emph{weighted entropy form} of $\Sigma$ by  $\hat{T}\equiv 0$ if $\Sigma$ is flat and by
\begin{equation*}
 \hat{T}=K_g T
\end{equation*}
otherwise.
In either case, the function $|\hat{T}|_g$ is continuous on $\Sigma$.
\end{lem}
\begin{proof}
 If $\Sigma$ is flat then there is nothing to prove as $\hat{T}$ is identically zero.  Otherwise, by Corollary \ref{DoublePoleCor}, $T$ is smooth away from the isolated poles where $K_g$ has a zero. As $\Sigma$ is smooth and $K_g\leq 0$, if $K_g(p)=0$ at a point $p$, then $\nabla_g K_g (p)=0$.  In particular $K_g=O(r^2)$ where $r$ is the distance to $p$. On the other hand, by Corollary \ref{DoublePoleCor}, $P$ has a double pole at $p$ and indeed $|T|_g=C r^{-2} +o(r^{-2})$ for some constant $C$.  Hence, $|\hat{T}|_g$ extends continuously across $p$.  Since $p$ was an arbitrary singularity of $T$, this proves the lemma.
\end{proof}
We now define the following weighted $L^{1/2}$ of the entropy differential $T$ of a smooth minimal surface $\Sigma$ as
\begin{equation*}
 ||T||^{1/2}_{\hat{L}^{1/2}(\Sigma)}= \int_{\Sigma} |\hat{T}|_g^{1/2} \mu_g= \int_{\Sigma} |T|_g^{1/2} \sqrt{-K_g} \mu_g.
\end{equation*}
If $\Sigma_0\subset \Sigma$, then we obviously have a domain monotonicity property
\begin{equation*}
  ||T||_{\hat{L}^{1/2}(\Sigma_0)}\leq  ||T||_{\hat{L}^{1/2}(\Sigma)}.
\end{equation*}
By Lemma \ref{weightedEntLem}, if $\Sigma$ is a smooth minimal surface and $\Sigma_0$ is pre-compact in $\Sigma$, then
\begin{equation*}
   ||T||_{\hat{L}^{1/2}(\Sigma_0)}<\infty
\end{equation*}
Finally, If $\hat{T}^\Sigma$ is the weighted entropy form of $\Sigma$, then  $\hat{T}^{\lambda \Sigma}=\lambda^{-2} \hat{T}^\Sigma$ is the weighted entropy form of $\lambda \Sigma$.  Hence, the weighted $L^{1/2}$ norm of the entropy form is scale invariant. That is
\begin{equation*}
||{T}^{\lambda \Sigma}||_{\hat{L}^{1/2}(\lambda \Sigma)} =||{T}^\Sigma||_{\hat{L}^{1/2}(\Sigma)}
\end{equation*}
for all $\lambda>0$.

\subsection{The Curvature Estimate}
We now use the scale invariance of $||T^\Sigma||_{\hat{L}^{1/2}}$ and Theorem \ref{MainCharacterizationThm} to prove an $\epsilon$-regularity result:
\begin{thm}\label{CurvEstThm}
 There are constants $\epsilon>0$ and a $C>0$ so that: if $\Sigma$ is a properly embedded minimal surface in $B_{2R}$ and
$$||T||_{\hat{L}^{1/2}(\Sigma)}<\epsilon,$$
 then
\begin{equation*}
 R^2 \sup_{B_{R}\cap \Sigma} |A|^2 \leq C^2.
\end{equation*}
\end{thm}
\begin{rem}
 The embeddedness condition is essential as can by seen by considering an appropriate rescaling of Enneper's surface.
\end{rem}

We begin with a Lemma which is crucial to the blow-up argument.
\begin{lem}\label{BlowUpLem}
  Fix $C>0$, $p \in \mathbb{R}^3$ and suppose $\Sigma$ is a properly embedded smooth surface in $B_{R}(p)\subset\mathbb R^3$ satisfying
\begin{equation*}
\sup_{{B}_{R}(p)\cap \Sigma} |A|^2 \geq 16 C^2 R^{-2}.
\end{equation*}
  Then there is a point $q\in \Sigma$ and scale $s>0$ so that $B_{Cs}(q)\subset {B}_{2R}(p)$ and
\begin{equation*}
 \sup_{B_{Cs}(q)\cap \Sigma} |A|^2\leq 4 s^{-2}=4|A|^2(q).
\end{equation*}
\end{lem}
\begin{proof}
 With $r(x)=|x-p|$ define the function
\begin{equation*}
 F(x)=\left(r(x)-\frac{3}{2}R\right)^2 |A|^2.
\end{equation*}
This is a Lipschitz function on $B_{\frac{3}{2}R}(p)\cap \Sigma$ that vanishes on $\partial B_{\frac{3}{2}R}(p)\cap \Sigma$.  As $F$ is continuous, non-negative and vanishes on $\partial B_{\frac{3}{2}R}(p)\cap \Sigma$, $F$ achieves its positive maximum at a point $q\in B_{\frac{3}{2}R}(p)\cap \Sigma$.

The lower bound
$$\sup_{{B}_{R}(p)\cap \Sigma} |A|^2 \geq 16 C^2 R^{-2} $$ implies that
$ F(q)\geq 4C^2.$
Set $s = |A|^{-1}(q)$ and
$ \sigma=\frac{3}{2}R- r(q) $
and note that $2 C s\leq \sigma$. Furthermore, if $x\in B_{\sigma/2}(q)$, then
$$r(x)\leq \frac{3}{2} R-\frac{\sigma}{2}<\frac{R}{2}$$ and so $\sigma^2 \leq 4(r(x)-\frac{3}{2}R)^2$ and $B_{\sigma/2}(q)\subset B_{\frac{3}{2} R}(p)$. Combining these facts,
\begin{equation*}
 \sup_{B_{Cs}(q)\cap \Sigma} \frac{\sigma^2}{4} |A|^2\leq  \sup_{B_{\sigma/2}(q)\cap \Sigma} \frac{\sigma^2}{4} |A|^2\leq\sup_{B_{\sigma/2}(q)\cap \Sigma} F\leq F(q)=\sigma^2 |A|^2(q).
\end{equation*}
Which verifies the claim.
\end{proof}

We also note the following well-known fact:
\begin{prop} \label{GeomArzAscoProp}
 Suppose that $R_i \nearrow \infty$ and that $\Sigma_i$ are properly embedded minimal surfaces in $B_{R_i}$ so that
\begin{enumerate}
 \item $0\in \Sigma_i$ and $|A^{\Sigma_i}|(0)=1$;
 \item For each $K\subset \subset \mathbb R^3$ if $K\subset B_{R_i}$, then
$$ \sup_{K\cap \Sigma_i} |A^{\Sigma_i}|\leq C(K)<\infty;$$
\end{enumerate}
then up to passing to a subsequence, the $\Sigma_i$ converge smoothly and with multiplicity one to a properly embedded minimal surface $\Sigma$ in $\mathbb R^3$ so that $0\in \Sigma$ satisfies $|A^\Sigma|(0)=1$.
\end{prop}
\begin{proof}
 Up to passing to a subsequence, the $\Sigma_i$ converge to a smooth  minimal lamination $\mathcal{L}$ of $\mathbb R^3$. As $0\in \Sigma_i$ for each $i$,  there is a leaf $L$ of the lamination containing $0$, moreover $|A^L|(0)=1$ and so $L$ is not flat. In particular, by the strong half-space theorem~\cite{Hoffman1990a} $L$ is the only leaf of $\mathcal{L}$.  Finally, if the convergence is with multiplicity greater than one, then $L$ would be stable and hence flat by~\cite{Fischer-Colbrie1980}.
\end{proof}

\begin{proof}[Proof of Theorem \ref{CurvEstThm}]
By rescaling we may take $R=1$.  Assume the theorem were false, then there was a sequence of minimal surfaces $\Sigma_i$ properly embedded in $B_2$ so  $||{T}^{\Sigma_i}||_{\hat{L}^{1/2}(\Sigma_i)}\to 0$ and $\sup_{B_{1}\cap \Sigma_i} |A^\Sigma_i|^2 \to \infty$.
By Lemma \ref{BlowUpLem}, there exist a sequence of $C_i\to \infty$, points $q_i \in \Sigma_i$ and scales $s_i\to 0$ so
$B_{C_i s_i}(q_i)\subset B_{2}$ and
 \begin{equation*}
 \sup_{B_{C_is_i}(q_i)\cap \Sigma_i} |A^{\Sigma_i}|^2\leq 4 s^{-2}_i=4|A^{\Sigma_i}|^2(q_i).
\end{equation*}
We set $\hat{\Sigma}_i=s^{-1}_i(\Sigma_i\cap B_{C s_i}(q_i)- q_i )$. The scaling properties of  $||{T}||_{\hat{L}^{1/2}(\Sigma)}$ and  domain monotonicity together imply that
$$||{T}^{\hat{\Sigma}_i}||_{\hat{L}^{1/2}(\hat{\Sigma}_i)}\leq ||{T}^{\Sigma_i}||_{\hat{L}^{1/2}(\Sigma_i)}\to 0.
$$
Moreover, each $ \hat{\Sigma}_i$ is properly embedded in $B_{C_i}$, contains $0$ and satisfies
 \begin{equation*}
 \sup_{B_{C_i}(0)\cap \hat{\Sigma}_i} |A^{\hat{\Sigma}_i}|^2\leq 4 =4|A^{\hat{\Sigma}_i}|^2(0).
\end{equation*}
Hence, by Proposition \ref{GeomArzAscoProp}, up to passing to a subsequence, the $\hat{\Sigma}_i$ converge to a smooth properly embedded minimal surface $\hat{\Sigma}$ in $\mathbb R^3$. The convergence is with multiplicity one and $0\in \hat{\Sigma}$ satisfies
$|A^{\hat{\Sigma}}|^2(0)=1$.  By smoothness of $\hat{\Sigma}$ and the monotonicity formula, there is a $\rho>0$ so that in $B_{\rho}(0)\cap \hat{\Sigma}$ one has $|A^{\hat{\Sigma}}|> \frac{1}{2}$ and so $\pi \rho^2<Area(\hat{\Sigma}\cap B_{\rho})\leq 2\pi \rho^2$. As the
$\hat{\Sigma}_i$ converge smoothly and with multiplicity one to $\hat{\Sigma}$, there is an $i_0$ large so that  $i>i_0$ implies
$|A^{\hat{\Sigma_i}}|> \frac{1}{4}$ and $\frac{\pi}{2} <Area(\hat{\Sigma}\cap B_{\rho})\leq 3\pi \rho^2$. As $A^{\Sigma_i}\neq 0$  , $\hat{T}^{\hat{\Sigma}_i}$ is smooth in $B_{\rho}\cap \hat{\Sigma}_i $ for $i>i_0$ and converges smoothly to $\hat{T}^{\hat{\Sigma}}$ in $B_{\rho}$.
However,
 $||{T}^{\hat{\Sigma}_i}||_{\hat{L}^{1/2}(B_\rho\cap \hat{\Sigma}_i)} \to 0$, hence $\hat{T}^{\hat{\Sigma}}\equiv 0$ on $B_{\rho}\cap \hat{\Sigma}$.  Together with $A^{\hat{\Sigma}}\neq 0$ on $B_\rho\cap \hat{\Sigma}$ this implies $T^{\hat{\Sigma}}\equiv 0$ on $B_{\rho}\cap \hat{\Sigma}$ and so $B_{\rho}\cap \hat{\Sigma}$ is contained in a finite rescaling Enneper's surface by Theorem \ref{MainCharacterizationThm}. It then follows from the strong unique continuation property of smooth minimal surfaces  that $\hat{\Sigma}$ is a finite rescaling of Enneper's surface, contradicting that $\hat{\Sigma}$ is properly embedded and proving the theorem.
\end{proof}

\subsection{Gap properties of the entropy form}\label{GapSec}
In light of Theorem \ref{CurvEstThm} an interesting question is to determine the optimal constant $\epsilon$ in Theorem \ref{CurvEstThm}.  This is equivalent to determining a lower bound for $||T||_{\hat{L}^{1/2}(\Sigma)}$ when $\Sigma$ is a non-flat properly embedded minimal surface in $\mathbb R^3$.  We present some partial results in this direction as well as pose a question about the expected behavior.

We begin with a simple computation:
\begin{lem}
 If $C$ is the catenoid then
\begin{equation*}
 ||{T}||_{\hat{L}^{1/2}(C)}=2\sqrt{2} \pi^4.
\end{equation*}
\end{lem}
\begin{proof}
 Weierstrass data for the catenoid is $(\mathbb{C}/\langle 2\pi i\rangle, J, - e^z, \d z)$ where $J$ is the usual complex structure. Using this data and writing $z=x+iy$, we have
$$
P=-\frac{1}{2} \d z^2, \quad g=\cosh^2 x \,\d z \circ \d\bar{z}, \quad K_g=-
\frac{1}{\cosh^4 x}.
$$
As a consequence, $|T|_g=\frac{\sqrt{2}}{2\cosh^2 x}$ and $|\hat{T}|_g=\frac{\sqrt{2}}{2\cosh^6 x}$, hence
\begin{equation*}
 ||{T}||_{\hat{L}^{1/2}(C)}^{1/2}= 2^{\frac{3}{4}}\pi \int_{-\infty}^\infty \frac{1}{\cosh x} dx=2^{\frac{3}{4}}\pi^2.
\end{equation*}
\end{proof}
A consequence of Theorem \ref{CurvEstThm} and~\cite{MeeksIIIa} is that the optimal surface must have finite total curvature:
\begin{prop}
 If $\Sigma$ is a properly embedded minimal surface in $\mathbb R^3$ and
$$ ||{T}||_{\hat{L}^{1/2}(\Sigma)}<\infty,$$
then
$$\int_{\Sigma} |A|^2 \mu_g<\infty.$$
\end{prop}
\begin{proof}
 Since $||{T}||_{\hat{L}^{1/2}(\Sigma)}<\infty$ there is a value $R>0$ so that
\begin{equation*}
 \int_{\Sigma\backslash B_R} |\hat{T}|_g^{1/2} \mu_g <\epsilon^{1/2}
\end{equation*}
where $\epsilon$ is given by Theorem \ref{CurvEstThm}. This implies that there is a constant $C>0$ given by Theorem \ref{CurvEstThm} so that for $p\in \Sigma \backslash \overline{B}_{2R}$ we have $B_{\frac{1}{2}|p|}{p}\subset \mathbb R^3\backslash B_{R}$ and so
\begin{equation*}
 |A|^2(p)\leq \frac{16 C^2}{|p|^2}.
\end{equation*}
That is, $\Sigma$ has quadratic extrinsic decay of curvature.  Hence, by Theorem 1.3 of~\cite{MeeksIIIa}, $\Sigma$ has finite total curvature.
\end{proof}
We pose the following question:
\begin{ques}
 Let $\Sigma$ be a non-flat properly embedded minimal surface in $\mathbb R^3$ and let $C$ be the catenoid.  Is it true that
\begin{equation*}
 ||T^{\Sigma}||_{\hat{L}^{1/2}(\Sigma)}\geq ||T^{C}||_{\hat{L}^{1/2}(C)}=2\sqrt{2}\pi^4
\end{equation*}
with equality only if $\Sigma$ is a catenoid?
\end{ques}

\subsection{Compactness properties for uniform bounds on $T$}
We conclude with a compactness result for sequences of properly embedded minimal surfaces $\Sigma_i$ which admit a uniform bound on the entropy differential.  This is a standard consequence of Theorem \ref{CurvEstThm} and the removable singularities result of~\cite{MeeksIIIa}.
\begin{thm}
 Let $\Sigma_i$ be a sequence of properly embedded minimal surfaces in an open region $\Omega\subset \mathbb R^3$ with entropy forms $T^{\Sigma_i}$ satisfying
$$||T^{\Sigma_i}||_{\hat{L}^{1/2}(\Sigma_i)}\leq \bar{C}<\infty.$$
 Then there is a subsequence of the $\Sigma_i$ and  a finite (possibly empty) set of points $p_1, \ldots, p_N\in \Omega$ so that:
\setcounter{thm}{1}
\begin{enumerate}
  \item \label{CpctItemOne} On each compact set $K\subset \subset \Omega\backslash \set{p_1, \ldots, p_N},$
     $$\sup_{K\cap \Sigma_i} |A|\leq C(K)<\infty;$$
\item\label{CpctItemTwo}  $\epsilon N^2<\bar{C}$ where $\epsilon>0$ is given by Theorem \ref{IntroCurvEstThm};
\item \label{CpctItemThree} The $\Sigma_i$ converge in $\Omega\backslash \set{p_1, \ldots, p_N}$ to a smooth minimal lamination $\mathcal{L}$ of $\Omega\backslash \set{p_1, \ldots, p_N}$.  Moreover, the closure of $\overline{\mathcal{L}}$ of $\mathcal{L}$ in $\Omega$ is a smooth lamination of $\Omega$.
\end{enumerate}
\end{thm}
\begin{proof}
 We define a sequence of Radon measures, $\mu_i$, on $\Omega$ by setting
\begin{equation*}
 \mu_i(U)=\int_{\Sigma_i\cap U} |\hat{T}^{\Sigma_i}|^{1/2} \mu_g
\end{equation*}
so
\begin{equation*}
 \mu_i(\Omega)=||T^{\Sigma_i}||_{\hat{L}^{1/2}(\Sigma_i)}^{1/2}\leq \bar{C}^{1/2}<\infty.
\end{equation*}
By the standard compactness theorem for Radon measures, up to passing to a subsequence, the $\mu_i$ weak* converge to a Radon measure $\mu$. It follows with Theorem \ref{CurvEstThm} that if for $p\in \Omega$ there is an $r>0$ so that $B_{2r}(p)\subset \Omega$ and $\mu(B_{2r}(p))< \frac{1}{2} \epsilon^{1/2}$, then there is a constant $C>0$ so that
\begin{equation*}
 \sup_{B_{r}(p)\cap \Sigma_i} |A|^2\leq \frac{C^2}{r^2}<\infty.
\end{equation*}
By standard covering arguments and the pigeonhole principle one concludes that there are at most $N$ points $p_1, \ldots, p_N\in \Omega$ with $N\epsilon^{1/2}<\bar{C}^{1/2}$ so that no such $r$ exists.  It follows that for any compact set $K\subset \Omega \backslash \set{p_1, \ldots, p_N}$ we have the curvature estimate:
 $$\sup_{K\cap \Sigma_i} |A|\leq C(K)<\infty.$$
This verifies Items \eqref{CpctItemOne} and \eqref{CpctItemTwo}.

To prove Item \eqref{CpctItemThree}, we note that the uniform curvature estimates of Item \eqref{CpctItemTwo} and standard compactness results -- see Appendix B of \cite{Colding2004} -- imply that, up to passing to a further subsequence, the $\Sigma_i$ converge in $\Omega \backslash \set{p_1, \ldots, p_N}$ to a smooth minimal lamination, $\mathcal{L}$, of $\Omega\backslash \set{p_1, \ldots, p_N}$.
We claim that near each $p_i$ the lamination has quadratic curvature decay.
To prove this we apply the Lebesgue decomposition theorem to $\mu$ and to $L_\Omega$, Lebesgue measure on $\mathbb R^3$ restricted to $\Omega$. This implies
\begin{equation*}
 \mu=\mu_{reg}+\mu_{sing}.
\end{equation*}
where $\mu_{reg}$ is absolutely continuous with respect to ${L}_\Omega$ while  $\mu_{sing}\perp {L}_\Omega $. In fact, the support of $\mu_{sing}$ is $\set{p_1, \ldots, p_n}$, because $\mathcal{L}$ is a lamination of $\Omega\backslash \set{p_1, \ldots, p_N}$.  Hence,
$$
\lim_{\rho\to 0} \mu(B_{\rho}(p_i)\backslash \set{p_i}) = \lim_{\rho\to 0} \mu_{reg}(B_{\rho}(p_i)\backslash \set{p_i}) = 0
$$
and so there is a $\delta>0$ so that $\mu(B_{2\delta}(p_i)\backslash \set{p_i})<\epsilon^{1/2}$.
Hence, for $p\in B_{\delta}(p_i)\backslash \set{p_i}$ Theorem \ref{CurvEstThm} yields
\begin{equation*}
 |A|^2(p)\leq \frac{4 C^2}{|p-p_i|^2}.
\end{equation*}
Theorem 1.2 of~\cite{MeeksIIIa} then implies that each $p_i$ is a removable singularity of $\mathcal{L}$ which concludes the proof of Item \eqref{CpctItemThree}.
\end{proof}


\section{Introduction}
Searching for a quantum theory of gravity has been a fruitful field for theoretical physics in the last century. In this sense, Ho\v{r}ava proposed in Ref. \cite{Horava:2008jf} a new quantizable theory of gravity, using the ideas from solid state physics. This theory was originally called Ho\v{r}ava-Lifshitz quantum gravity, that is a power-counting renormalizable theory with consistent ultraviolet (UV)  behavior and, on the other hand, the  theory has one fixed point in the infrared limit (IR) namely General Relativity (GR) \cite{Horava:2008jf, Horava:2008ih, Horava:2009if, Horava:2009uw}. Thus, this theory leads to a modification of the Einstein's general relativity at  high energies producing interesting features in cosmology.   The first ideas about this subject have been presented in Refs. \cite{Calcagni:2009ar, Kiritsis:2009sh, Takahashi:2009wc} where the cosmological consequences of Ho\v{r}ava-Lifshitz (HL)  cosmology were studied in a vast detail,  see also Refs. \cite{Nojiri:2010wj, Clifton:2011jh, Brandenberger:2009yt, Mukohyama:2009gg, Cai:2009dx, Saridakis:2009bv, Mukohyama:2009zs, Mukohyama:2009mz, Cai:2009in, Wang:2009azb, Leon:2009rc, Minamitsuji:2009ii, Carloni:2009jc, Gao:2009wn, Kobayashi:2010eh, Maeda:2010ke, Saridakis:2011pk, Bertolami:2011ka}. The new theory has virtues and defects, and its basis  ideas are not free of controversy \cite{Henneaux:2009zb}. In this article we would like to address  an interesting point concerning HL cosmology that is still an open problem. In particular, we are seeing a behavior of HL cosmology for different kinds of matter content of the universe such as dust, cosmological constant and phantom matter. The present paper is organized as follows. In Sec. I we discuss the HL cosmology for different kinds of the matter and also we present general settings and constrains in the theory. In Sec. II we analyze the role of stiff-like matter taken as the dominant at early times. Finally, in Sec. III we give some final remarks.

\section{Ho\v{r}ava-Lifshitz Cosmology}
In this section we would like to review the main results of Hl cosmology and its implications for different kind of mater. Also, we are doing the comparison between the behaviors of HL and GR cosmology.
Let us first to introduce a cosmological model under
consideration based on Ref. \cite{Mukohyama:2010xz}.
It is well-known that in Ho\v{r}ava-Lifshitz gravity we do not have the
Hamiltonian constraint and therefore, from the formal cosmological point of
view, also we do not have one Friedmann equation as in standard cosmology. In this theory the starting point is the dynamical equation for flat Friedmann-Robertson-Walker (FRW) spacetime
\begin{equation}
\eta \left( 2\dot{H}+3H^{2}\right) =-p,  \label{eq.1}
\end{equation}%
where we introduced a parameter
\begin{equation}
\eta =\frac{1}{2}\left( 3\lambda -1\right) ,  \label{eq.2}
\end{equation}%
and $\lambda $ represents a dimensionless constant. The value of $\lambda$ is fixed by the diffeomorsphism invariance of four-dimensional general relativity. At this point, we would like to emphasize that in HL cosmology any value of the constant $\lambda$ is consistent with a foliation that preserves the diffeomorphism invariance. For example,  it was shown in Ref. \cite{Gumrukcuoglu:2011xg} that for $1/3<\lambda
<1$ (or equivalently $0<\eta <1$) the scalar graviton is a ghost and  is
ghosts-free if $\lambda <1/3$ (that is $\eta <0$) and  $\lambda
>1$  (or $\eta >1$).
The field equation for matter field is described satisfies a
non-conservation equation, and at high energy limit, or early times limit, we
can write it in the following form
\begin{equation}
\dot{\rho}+3H(\rho +P)=-Q,  \label{eq.3}
\end{equation}
where $Q$ represents the rate of energy non-conservation, and  the low energy limit  can be recovered only if $Q\rightarrow
0 $. We can conclude that the energy non-conservation
 is an effect of high-energy physics. Now,  equations (\ref{eq.1})
and (\ref{eq.3}) contain enough information  to describe a
cosmological evolution of the spacetime, together with the equation
of state (EoS) $p=\omega \rho $. Indeed, using this dynamical equation plus the
non-conservation equation we can obtain one first integral in the form
\begin{equation}
3H^{2}=\frac{1}{\eta }\rho +\frac{C\left( t\right) }{a^{3}},  \label{eq.4}
\end{equation}%
and this result  can be interpreted as a Friedmann equation in HL gravity where
the integration "constant" is given by
\begin{equation}
C\left( t\right) =C_{0}+\frac{1}{\eta }\int_{t_{0}}^{t}d\tau a^{3}\left(
\tau \right) Q\left( \tau \right) .  \label{eq.c}
\end{equation}%
In order to obtain a cosmological description from this kind of theories, we
would like to discuss first some well known results from GR in the context of HL
cosmology. The subscript $0$ indicates quantities today ($t_{0}$).

 \textbf{ $\bullet$ \,Dust matter:\ $p=0\longleftrightarrow \omega =0$
} \vspace{0.5cm}

 In this case from (\ref{eq.1})
we find exactly the same solution for the Hubble parameter as in GR, i.e.,
\begin{equation}
H\left( t\right) =H_{0}\left[ 1+\frac{3}{2}H_{0}\left( t-t_{0}\right) \right]
^{-1},  \label{eq.5}
\end{equation}%
and for the cosmic scale factor reads
\begin{equation}
a\left( t\right) =a_{0}\left[ 1+\frac{3}{2}H_{0}\left( t-t_{0}\right) \right]
^{2/3}.  \label{eq.6}
\end{equation}%
It is well known in GR, when we have $\omega =0$,  an evolution is described
by an energy density proportional to $a^{-3}$ ( this result is obtained from the conservation energy density equation $\dot{\rho}+3H(\rho +P)=0$)
and the formal solution of the cosmological evolution in GR is
identical to Eqs. (\ref{eq.5}) and (\ref{eq.6}), obtained for HL cosmology.  Although the formal solutions for the Hubble parameter and the scale factor from HL and GR are identical in this case, the dynamical behavior of the energy density is absolutely different, because, in
HL cosmology the energy density is not fixed by the non-conservation
equation  (\ref{eq.3}). Finally for the dust matter case, we would like to mention that (\ref%
{eq.5}) and (\ref{eq.6}) are also solutions of (\ref{eq.1}) if we set $\rho =0$
and $\eta $ is finite. Thus,  these equations can also be seen as a
self-decelerated evolution.

In a nutshell, both $\omega =0$ and $\rho =0$  lead
 to the same solutions described by Eqs. (\ref{eq.5}) and (\ref{eq.6}). Therefore, in  case of dust matter we  found that the effect of high energy
physics is reflected in the fact that during the early stages, the rate of energy non-conservation is not constant  $Q(t)\neq 0$,
and we have an interchange of energy the between the mater  of the
universe and some unknown source.

 \textbf{ $\bullet$ \,Cosmological constant: \textit{$\omega =-1$. }
} \vspace{0.5cm}

In this case, the conservation
law

\begin{equation}
\dot{\rho}=-Q, \label{constant1}
\end{equation}
therefore the behavior of cosmological constant at early universe would be different that the respective one in GR. Then, the interchange of energy between the matter ( described in this case for the cosmological constant kind) and the unknown source gives a non constant behavior for the energy density for our cosmological constant type. We note that  the GR  limit is  recovered if we take  $Q\left( a\rightarrow \infty \right)
\rightarrow 0$, i.e., $\rho =const.$. Then the "role" of the
cosmological constant in HL gravity is different with the respective one in GR.

 \textbf{ $\bullet$ \,Phantom Evolution.
} \vspace{0.5cm}

 In this case, from (\ref{eq.1}) we
can write the following formal solution for the Hubble parameter
\begin{equation}
H^{2}\left( a\right) =\frac{1}{a^{3}}\left( a_{0}^{3}H_{0}^{2}-\frac{\omega
}{\eta }\int_{a_{0}}^{a}daa^{2}\rho \left( a\right) \right) ,  \label{eq.16s}
\end{equation}%
and if we take one phantom Ansatz for the energy density given by $\rho \left( a\right) =\rho _{0}\left( a/a_{0}\right)
^{\beta }$ with  $\beta >1$, we get  the explicit form of the formal solution
\begin{equation}
H^{2}\left( a\right) =H_{0}^{2}\left( \frac{a_{0}}{a}\right) ^{3}\left[ 1-%
\frac{\omega \rho _{0}/3H_{0}^{2}}{\eta \left( 1+\beta /3\right) }\left[
\left( \frac{a}{a_{0}}\right) ^{\beta +3}-1\right] \right] ,  \label{eq.18s}
\end{equation}%
and
the future behavior of the Hubble parameter reads as follow
\begin{equation}
H^{2}\left( a\rightarrow \infty \right) =-\frac{\omega \rho _{0}/3}{\eta
\left( 1+\beta /3\right) }\left( \frac{a_{0}}{a}\right) ^{\beta }.
\end{equation}
One realistic model implies the following restriction for the barotropic index, $\omega <0$. If we use the following setting
\begin{equation}
\left\vert \omega \right\vert \rho _{0}/3H_{0}^{2}=\left( 1+\beta /3\right)
\eta ,  \label{eq.19s}
\end{equation}%
we obtain one phantom solution for the scale factor given by
\begin{equation}
a\left( t\right) =a_{0}\left( \frac{2}{\beta H_{0}}\right) ^{2/\beta }\left(
t_{s}-t\right) ^{-2/\beta }\text{ \ \ }and\text{ \ \ }t_{s}=t_{0}+\frac{2}{%
\beta H_{0}}.  \label{eq.20s}
\end{equation}%
It is worthwhile noticing that the solution (\ref{eq.20s}) has a new region for
the phantom evolution for the barotropic index $\omega <0$ instead of the
usual, more restrictive one $\omega <-1$ that appears in GR.

We  also notice that if $\eta
\rightarrow \infty$, then $\beta \rightarrow -3$ and we can write
\begin{equation}
H^{2}\left( a\right) _{\eta \rightarrow \infty }\rightarrow H_{0}^{2}\left(
\frac{a_{0}}{a}\right) ^{3},
\end{equation}
so that the phantom approach has no sense. According to
\begin{equation}
\eta =\left\vert \omega \right\vert \left( \rho _{0}/3H_{0}^{2}\right)
\left( 1+\beta /3\right) ^{-1}>1, \label{lol}
\end{equation}
or
\begin{equation}
\left\vert \omega \right\vert \left( \rho _{0}/3H_{0}^{2}\right) >1+\beta
/3,
\end{equation}
we conclude that the observational data do not allow fulfill the last inequality. On the other hand, if $%
0<\eta <1$ (ghost scalar graviton) we have $\left\vert \omega \right\vert
\left( \rho _{0}/3H_{0}^{2}\right) <1+\beta /3$ and in this case we can have
a phantom phase provided  we can accept an existence of a ghost scalar graviton.

\textbf{ $\bullet$ \,General Settings.
} \vspace{0.5cm}
We are now interested in  the limit  where the matter sector is decoupled from the gravity sector, where dark matter as an integration constant dominated the evolution of the universe, that its $\eta \rightarrow \infty \left( \lambda
\rightarrow \infty \right) $ \cite{Gumrukcuoglu:2011xg}, in this limit  we can
distinguish several situations. For instance, if $p$ is finite and $\eta
\rightarrow \infty $, from (\ref{eq.1}) we obtain exactly the same solution  given by (\ref{eq.5})
and (\ref{eq.6}) and the Hubble parameter becomes in this limit
\begin{equation}
3H^{2}=\frac{C_{0}}{a^{3}}, \label{heta}
\end{equation}
thus, a dust-like evolution but $p$  and the barotropic index $\omega$ do not vanish. Also, we would like to notice that if $t\rightarrow \infty $ and $\eta
$ is finite then
\begin{equation}
3H^{2}\left( t\rightarrow \infty \right) \rightarrow \frac{1}{\eta }\rho
\left( t\rightarrow \infty \right) +\frac{1}{a^{3}\left( t\rightarrow \infty
\right) }\left[ C_{0}+\frac{1}{\eta }\int_{t_{0}}^{\infty }d\tau a^{3}\left(
\tau \right) Q\left( \tau \right) \right] =\frac{1}{\eta }\rho \left(
t\rightarrow \infty \right) , \label{hgrand}
\end{equation}
if at least $Q\left( t\right) $ decreases as $a^{-4}\left( t\right) $ when $%
a\left( t\rightarrow \infty \right) \rightarrow \infty $.  This condition on $Q$ is consistent with the recovering of local invariance for the matter sector, where it is demanded that  $Q\left( a\rightarrow
\infty \right) \rightarrow 0$ and $C\left( t\rightarrow \infty \right)
\rightarrow C_{0}$ (see \cite{Gumrukcuoglu:2011xg}). So, we have the following relation between observational parameters
\begin{equation}
3H_{0}^{2}=\frac{1}{\eta }\rho _{0}+\frac{C_{0}}{a_{0}^{3}}.
\end{equation}

On the other hand, from  (\ref{eq.1}) and (\ref{eq.4}) we can write
\begin{equation}
\dot{H}=-\frac{1}{2}\frac{C_{0}}{a^{3}}-\frac{1}{2\eta }\left[ \left(
1+\omega \right) \rho +\frac{1}{a^{3}}\int_{t_{0}}^{t}d\tau a^{3}\left( \tau
\right) Q\left( \tau \right) \right] , \label{hdot}
\end{equation}
so that if we take the limit $\eta \rightarrow \infty $, and keep finite second term in (\ref{hdot}), we obtain

\begin{equation}
H\left( a\right) =H_{0}%
\sqrt{1+\frac{C_{0}a_{0}^{-3}}{3H_{0}^{2}}\left[ \left( \frac{a_{0}}{a}%
\right) ^{3}-1\right] }. \label{hdot2}
\end{equation}
Furthermore, if $C_{0}a_{0}^{-3}/3H_{0}^{2}<1$, we find

\begin{equation}
H\left( a\rightarrow \infty \right) \rightarrow H_{0}\sqrt{1-\frac{%
C_{0}a_{0}^{-3}}{3H_{0}^{2}}}<H_{0}. \label{hlimit}
\end{equation}
To conclude, in  case with  $\eta \rightarrow \infty $ we find a de Sitter phase at late times. We also notice that this behavior can also be  obtained  if
we choose $\omega =-1$ but keep $\eta $ finite, that is,
\begin{equation}
\dot{H}=-\frac{1}{2}\frac{C\left( t\right) }{a^{3}}\rightarrow \dot{H}\left(
t\rightarrow \infty \right) \rightarrow 0,
\end{equation}
given that in this limit $C\left( t\rightarrow \infty \right) \rightarrow
C_{0}$. In another words, taking either  $\eta \rightarrow \infty $, for all $t$, and or $t\rightarrow \infty
$ and $\eta $ finite, we arrive to a de Sitter phase at late times.

 To complete this point we would like to  comment about the acceleration of universe. By using the expression for $\dot{H}$, as well as (\ref{eq.4}), we can
write
\begin{equation}
\dot{H}+H^{2}=-\frac{1}{6\eta }\left[ \left( 1+3\omega \right) \rho +\eta
\frac{C\left( t\right) }{a^{3}}\right] ,
\end{equation}
and it is straightforward to check that $\eta =1$ plus $C\left( t\right) =0$ implies the
standard expression for the acceleration in GR. Also, if $\dot{H}%
+H^{2}>0$ we must have $\omega <-1/3$ , to avoid violating  the weak
energy condition (WEC) $\rho >0$. Then, we can have quintessence, cosmological
constant or phantom schemes. Now, if $\dot{H}+H^{2}<0$ then $\omega >-1/3$  the WEC is fulfilled, then we have an evolution driven by dark matter.
Finally when $\dot{H}+H^{2}=0$ we obtain the solution
\begin{equation}
H\left( t\right) =H_{0}\left[ 1+H_{0}\left( t-t_{0}\right) \right] ^{-1},
\end{equation}
i.e., $\omega =-1/3$ (string gas) as in GR. As a curiosity we notice that
\begin{equation}
\left( \dot{H}+H^{2}\right) _{\omega =-1/3}=\frac{1}{3}\left( \dot{H}\right)
_{\omega =-1}.
\end{equation}
Finally, by setting  $\dot{H}+H^{2}=0$, we can write the following
expression for the density of energy
\begin{equation}
\rho \left( t\right) =\rho _{0}\left( a/a_{0}\right) ^{-3}+\frac{1}{\eta
\left\vert 1+3\omega \right\vert }\frac{1}{a^{3}}\int_{t_{0}}^{t}d\tau
a^{3}\left( \tau \right) Q\left( \tau \right) ,
\end{equation}
where $C_{0}=\left\vert 1+3\omega \right\vert \rho _{0}a_{0}^{-3}$ and $%
\omega <-1/3$, and we find that $\rho \left( t\rightarrow \infty \right)
\rightarrow \rho _{0}\left( a/a_{0}\right) ^{-3}$, i. e., a dust-like
behavior at late times is driven by a $\omega $ that is dark energy-like.

\section{Stiff-like matter as the dominant term at early times}
When the higher curvature terms are included in the cosmological evolution the dynamic equation for the Hubble parameter is given by \cite{Mukohyama:2010xz},
\begin{equation}
\eta \left( 2\dot{H}+3H^{2}\right) =-\omega \rho +\frac{\alpha k^{3}}{a^{6}} +\frac{\alpha' k^{2}}{a^{4}}-\frac{k}{a^{2}}+ \Lambda,  \label{eq234}
\end{equation}
where $\alpha$ and $\alpha'$ are constants and $k$ is the spatial
curvature. If we consider the early time limit, the dominant contribution is the term of the form  $\sim a^{-6}$, which could be interpreted as stiff matter like, and also would be explain by the supposition that stiff matter could be one important matter content at the very early universe  \cite{Banks:2004eb, Banks:2008ep}.
We write (\ref{eq234}) in the form
\begin{equation}
\eta \left( 2\dot{H}+3H^{2}\right) =-\omega \rho +\frac{\alpha k^{3}}{a^{6}}%
\left[ 1+\frac{a^{2}}{\alpha }\left( \alpha `k-a^{2}\right) \right] ,
\label{eq234b}
\end{equation}%
and the incorporation of $\Lambda $ will be discussed later.
Therefore, we have
 a new equation to handle that is
\begin{equation}
3\eta H^{2}=\rho -\frac{\alpha k^{3}}{a^{6}}.  \label{eq.32ss}
\end{equation}
Similarly as before, the parameter $\eta $ can  take any value. The dynamical equation then can be written as
\begin{equation}
\eta \left( 2\dot{H}+3H^{2}\right) =-\omega \rho +\frac{\alpha k^{3}}{a^{6}},
\label{eq.33ss}
\end{equation}%
The aproximation at early times is justified if the scale factor is kept
around $a^{2}\sim \left\vert \alpha `\right\vert $, if $k=1$ and for $k=-1$
we have a more restrictive condition: $a^{2}\left( \alpha `+a^{2}\right)
<<\alpha $ . From Eqs. (\ref{eq.32ss}) and (\ref{eq.33ss}) we obtain the
conservation law (the term proportional to $C\left( t\right) /a^{3}$,
i.e., proportional to $Q$, is not present here given that only we are
holding the dominant contribution $\sim a^{-6}$)
\begin{equation}
\dot{\rho}+3H\left( 1+\omega \right) \rho =0,  \label{eq.34ss}
\end{equation}%
as well as the equation
\begin{equation}
\dot{H}+3H^{2}=\frac{1}{2\eta }\left( 1-\omega \right) \rho,
\label{eq.35ss}
\end{equation}%
and from these last equations we note that if we do $\omega =1$, we obtain
the usual solution as in GR, i.e. $\rho \sim a^{-6}$ and. We find the
following solution for the Hubble parameter and the scale factor,
respectively,
\begin{equation}
H\left( t\right) =H_{0}\left[ 1+3H_{0}\left( t-t_{0}\right) \right] ^{-1},
\label{eq.36ss}
\end{equation}%
and
\begin{equation}
a\left( t\right) =a_{0}\left[ 1+3H_{0}\left( t-t_{0}\right) \right] ^{1/3}.
\label{eq.37ss}
\end{equation}

In order to see the implications of this result, we are comparing it to standard
GR, and the equations $3H^{2}=\rho $ and $\dot{\rho}+6H\rho =0$ ($\omega =1$)
that lead (\ref{eq.36ss}) and (\ref{eq.37ss}). Thus replacing (\ref{eq.36ss})
and (\ref{eq.37ss}) in (\ref{eq.32ss}), we obtain for $\omega =1$
\begin{equation}
3\eta H^{2}=\left( \rho _{0}-\frac{\alpha k^{3}}{a_{0}^{6}}\right) \left(
\frac{a_{0}}{a}\right) ^{6},   \label{eq.38ss}
\end{equation}%
and for $k=1$ we must to have $\rho _{0}a_{0}^{6}>\alpha $ and $\eta >0$ (do
not forget that if $0<\eta <1$ we have a ghost scalar graviton and If $\eta
>1$ we are ghosts-free). If $\rho _{0}a_{0}^{6}<\alpha $, then $\eta <0$ and
we are ghosts-free. For $k=-1$ we must have $\eta >0$.
We want now to discuss an early universe that contains a mixture of stiff
matter-like and cosmological constant,
\begin{equation}
3\eta H^{2}=\rho -\frac{\alpha k^{3}}{a^{6}}+\Lambda ,    \label{s0}
\end{equation}%
and
\begin{equation}
\eta \left( 2\dot{H}+3H^{2}\right) =-\omega \rho +\frac{\alpha k^{3}}{a^{6}}%
+\Lambda .  \label{s1}
\end{equation}%
Then we obtain the $k$-independent equation
\begin{equation}
\dot{H}+3H^{2}=\frac{1}{2\eta }\left[ \left( 1-\omega \right) \rho +2\Lambda %
\right] ,    \label{s2}
\end{equation}%
and for $\omega =1$ and $\eta >0$, the formal solution is
\begin{equation}
H\left( t\right) =\sqrt{\Lambda /3\eta }\left( \frac{1+\Delta _{0}\exp \left[
-2\sqrt{3\Lambda /\eta }\left( t-t_{0}\right) \right] }{1-\Delta _{0}\exp %
\left[ -2\sqrt{3\Lambda /\eta }\left( t-t_{0}\right) \right] }\right) ,
  \label{s3}
\end{equation}%
where we have denoted
\begin{equation}
\Delta _{0}=\left( H_{0}-\sqrt{\Lambda /3\eta }\right) \left( H_{0}+\sqrt{%
\Lambda /3\eta }\right) ^{-1},    \label{s4}
\end{equation}%
and we observe that $\Delta _{0}=0\rightarrow H=\sqrt{\Lambda /3\eta }$,
i.e., an usual early de Sitter phase (old inflation-like). If $\Delta
_{0}\neq 0$, the solution (40) is a reasonable one solution for $t<<t_{0}$
(very early times), i. e., $H\left( t<<t_{0}\right) \rightarrow $constant
and we note also that there is no singularity in $H\left( t\right) $ for $%
\Delta _{0}>0$. For completeness, with $\omega =1$ and $\eta <0$, the Eq.
(39) has the solution
\begin{equation}
H\left( t\right) =\sqrt{\Lambda /3\left\vert \eta \right\vert }\tan \left[
\arctan \left( \frac{H_{0}}{\sqrt{\Lambda /3\left\vert \eta \right\vert }}%
\right) -6\sqrt{\Lambda /3\left\vert \eta \right\vert }\left( t-t_{0}\right) %
\right] ,    \label{s5}
\end{equation}%
and this solution will be reasonable only if $t<<t_{0}$ (very early times
and, for instance, $0<\arctan \left( H_{0}/\sqrt{\Lambda /3\left\vert \eta
\right\vert }\right) +6\sqrt{\Lambda /3\left\vert \eta \right\vert }%
t_{0}<\pi /2\Longrightarrow H\left( t\right) >0$), i. e., $H\left(
t<<t_{0}\right) \rightarrow $constant. Thus, even the presence of $\omega =1$%
, is a cosmological constant the dominant component at early times. In other
words, even accepting the presence of stiff matter at early times, in HL
cosmology this is a little relevant fact (nevertheless, a fluid for which $%
p=\rho $ could play a relevant role at early times if we are thinking in to
build an holographic approach to cosmology, see [29] and [30]).
Finally, we compare the deceleration parameter as given in GR and HL
cosmology. This parameter is defined by $q=-\left( 1+\dot{H}/H^{2}\right) $
so that, in GR reads

\begin{equation}
q=\frac{1}{2}\left( 1+3\omega \right) -\frac{1}{2}\left( 1+\omega \right)
\frac{\Lambda }{H^{2}},  \label{sl1}
\end{equation}%
and in HL cosmology

\begin{equation}
q=\frac{1}{2}\left( 1+3\omega \right) -\left( 1-\omega \right) \frac{\alpha
k^{3}}{\eta H^{2}a^{6}}-\frac{1}{2}\left( 1+\omega \right) \frac{\Lambda }{%
\eta H^{2}},  \label{sl2}
\end{equation}%
and for $\omega =1$ we have $q=2-\Lambda /H^{2}$ and $q=2-\Lambda /\eta
H^{2} $, respectively, and we note that when $\omega =1$, the curvature term
$\sim \alpha k^{3}$ disappears.

\section{Final Remarks}

We have studied some aspects of the HL cosmology where the emphasis has been
put on discuss some cosmological solutions which are present too in GR, in
particular, an evolution driven by dust is the same in both theories (flat
case). A possible phantom stage has been discussed where we have found a
condition less restrictive over the $\omega $-parameter ($\omega <0$, not $%
\omega <-1$). At late times, the energy density exhibits a like dust
behavior nevertheless the $\omega $-parameter satisfies the inequality $%
\omega <-1/3$. If we do infinite the parameter which preserves the
diffeomorphism invariance in the present theory, a late de Sitter phase can
be obtained without consider the usual scheme $\omega =-1$. The combined
effect of a curvature-dependent term, which is stiff matter reminiscent and
dominant at early times, beside the usual stiff matter component ($\omega =1$%
, possible important role at early evolution) and a cosmological constant
has been discussed and solutions has been found which exhibit an early de
Sitter phase and this fact shows that in HL cosmology the role of stiff
matter is of little relevance.

\begin{acknowledgments}
This work has been supported by COMISI\'ON NACIONAL DE CIENCIAS Y TECNOLOG\'IA
through FONDECYT Grants 1110076 (JS and SL), 1090613 and 1110230 (JS). This work was also
partially supported by PUCV-VRIEA grant No. 037.492/2013  (SL) and PUCV grant No. 123.713/2012.
(JS).
\end{acknowledgments}