I'm a physical chemist and I am involved in “colloidal dice”. These are small, cube-like particles with a really nice, regular shape. These particles are not really cubic, but more rounded, much like a dice. I've got a neat way to quantify their size and "roundness" and I'm interested in their volume and surface area. I've already found an expression for the former, but I'm still looking for the latter.


The (three-dimensional) particles have a diameter $r$ and can be described as being a special case of superellipsoids. My superellipsiod is a centrally symmetric particle defined by


\begin{equation}\label{eq:superellipsoid}
\left| \frac{x}{r} \right|^{m} +\left| \frac{y}{r} \right|^{m} + \left| \frac{z}{r} \right|^{m} \leq 1
\end{equation}


where $x$, $y$ and $z$ are Cartesian coordinates, $r$ is the radius and $m$ is the deformation parameter. Here, $2 \leq m \leq \infty$, where $m=2$ represents a sphere and $m=\infty$ a sharp cube. Most of my dice have a roundness $m=3.5$ and a radius $r\sim80$ nm.



> 
> The question: What is the surface are of a centrally symmetric superellipsoid, given the radius $r$ and deformation parameter $m$. A complete proof would be appreciated.
> 
> 
> 


Additional information: The volume of such a superellipsoid is given by


\begin{equation}\label{eq:volumecuboid}
V(r,m) = 8 r^{3} \frac{\displaystyle\left[ \Gamma \left(1+\frac{1}{m}\right) \right]^{3}}{\displaystyle\Gamma \left(1+\frac{3}{m}\right)}
\end{equation}


and the $\Gamma$-function has the property that


\begin{equation}
\Gamma(n) = (n-1)!
\end{equation}


where $n$ must be a positive integer.


The two important properties of the gamma function which will be used here are $x \Gamma(x)=\Gamma(x+1)$ and $\Gamma(\frac{1}{2})=\sqrt{\pi}$. Using these properties, it can be shown that for a sphere ($m=2$) with radius $r$:


\begin{equation}
V\_{\textrm{sphere}} (r,2) = \frac{4}{3}\pi r^{3}
\end{equation}


and for a cube ($m = \infty$) with edges $2r$


\begin{equation}
V\_{\textrm{cube}} (r,\infty) = (2r)^{3} .
\end{equation}


But now for the surface area. Thank you for reading.



Surface area of superellipsoid (dice)
You can't expect a closed formula for this surface area. The perimeter of an ellipse, much less the perimeter of a superellipse or the surface area of an ellipsoid or a superellipsoid, is already an integral that doesn't have a formula in the usual sense of an elementary formula. Instead, people did what they always do when an integral is meaningful but has no formula: They named it. That one is called an elliptic integral, but it's just for the case $m=2$. For the $m=2$ case, there is a way out in 3 dimensions: After you define elliptic functions to be the answer in 2 dimensions (for not just the full perimeter but the perimeter of arcs), you can express the 3-dimensional surface area in terms of them without naming any more functions. However, I don't expect that to happen for general $m$, which in any case only reduces the dimension of the problem by 1 instead of solving it.


It's true that your superellipsoids are in the special case that the three semiradii are equal, but I don't think that that rescues you for general values of $m$. Obviously it does rescue you when $m=2$.


What you really want to do for your problem is integrate numerically. You might want a convenient initial change of variables such as
$$(a,b,c) = (|x/r|^m,|y/r|^m,|z/r|^m),$$
so that you then get an integral over the triangle $a+b+c=1$. After that there are a variety of numerical integration methods that converge quickly. For instance, you can do the double integral by doing Gaussian integration twice. This particular strategy is not robust as $m \to \infty$, but you didn't say that you needed that limit. An alternative which I guess is robust in that limit, in the case that the semiradii are all equal, is to integrate over one 1/24 region of the spherical angle with the change of coordinates:
$$(x,y,z) \propto (1,u,v)$$
with $0 \le u \le v \le 1$. That is also an integral over a triangle.


