  Degenerate Fermi gases of atoms near a Feshbach resonance show universal
thermodynamic properties, which are here calculated with the geometry of
thermodynamics, and the thermodynamic curvature $R$. Unitary thermodynamics is
expressed as the solution to a pair of ordinary differential equations, a
"superfluid" one valid for small entropy per atom $z\equiv S/N k_B$, and a
"normal" one valid for high $z$. These two solutions are joined at a
second-order phase transition at $z=z_c$. Define the internal energy per atom
in units of the Fermi energy as $Y=Y(z)$. For small $z$, $Y(z)=y_0+y_1
z^{\alpha}+y_2 z^{2 \alpha}+\cdots,$ where $\alpha$ is a constant exponent,
$y_0$ and $y_1$ are scaling factors, and the series coefficients $y_i$ ($i\ge
2$) are determined uniquely in terms of $(\alpha, y_0, y_1)$. For large $z$ the
solution follows if we also specify $z_c$, with $Y(z)$ diverging as $z^{5/3}$
for high $z$. The four undetermined parameters $(\alpha,y_0,y_1,z_c)$ were
determined by fitting the theory to experimental data taken by a Duke
University group on $^6$Li in an optical trap with a Gaussian potential. The
very best fit of this theory to the data had $\alpha=2.1$, $z_c=4.7$,
$y_0=0.277$, and $y_1=0.0735$, with $\chi^2=0.95$. The corresponding Bertsch
parameter is $\xi_B=0.462(40)$.
