  We study minimisers of the $p$-conformal energy functionals, \[
\mathsf{E}_p(f):=\int_\ID \IK^p(z,f)\,dz,\quad f|_\IS=f_0|_\IS, \] defined for
self mappings $f:\ID\to\ID$ with finite distortion and prescribed boundary
values $f_0$. Here \[ \IK(z,f) = \frac{\|Df(z)\|^2}{J(z,f)} =
\frac{1+|\mu_f(z)|^2}{1-|\mu_f(z)|^2}\] is the pointwise distortion functional
and $\mu_f(z)$ is the Beltrami coefficient of $f$. We show that for
quasisymmetric boundary data the limiting regimes $p\to\infty$ recover the
classical Teichm\"uller theory of extremal quasiconformal mappings (in part a
result of Ahlfors), and for $p\to1$ recovers the harmonic mapping theory.
  Critical points of $\mathsf{E}_p$ always satisfy the inner-variational
distributional equation \[ 2p\int_\ID
\IK^p\;\frac{\overline{\mu_f}}{1+|\mu_f|^2}\varphi_\zbar \; dz=\int_\ID \IK^p
\; \varphi_z\; dz,\quad\forall\varphi\in C_0^\infty(\ID ). \] We establish the
existence of minimisers in the {\em a priori} regularity class
$W^{1,\frac{2p}{p+1}}(\ID)$ and show these minimisers have a pseudo-inverse - a
continuous $W^{1,2}(\ID)$ surjection of $\ID$ with $(h\circ f)(z)=z$ almost
everywhere. We then give a sufficient condition to ensure $C^{\infty}(\ID)$
smoothness of solutions to the distributional equation. For instance
$\IK(z,f)\in L^r_{loc}(\ID)$ for any $r>p+1$ is enough to imply the solutions
to the distributional equation are local diffeomorphisms. Further $\IK(w,h)\in
L^1(\ID)$ will imply $h$ is a homeomorphism, and together these results yield a
diffeomorphic minimiser. We show such higher regularity assumptions to be
necessary for critical points of the inner variational equation.
