\section{Introduction}
\label{sec:intro}
As is well known \cite{jacksonEM,smythe}, a solenoid of infinite
length will generate a perfectly uniform axial magnetic field
everywhere inside of a cylindrical volume.  As is also well known
\cite{jacksonEM,smythe}, an infinitely-long ``$\cos\theta$ coil''
(i.e., a cylindrical coil with a continuous surface current $\vec{K}
= K \cos\theta \hat{z}$, with the angle $\theta$ defined relative
to the $\hat{y}$-axis as shown in
Fig.\ \ref{fig:cos_theta_coil_schematic}), will generate a perfectly
uniform transverse magnetic field (along the $\hat{x}$-direction in
Fig.\ \ref{fig:cos_theta_coil_schematic}) within a cylindrical volume.
And indeed, experiments requiring a highly uniform magnetic field have
generally employed discretized and finite-length solenoidal or
$\cos\theta$ coils for the generation of highly uniform axial or
transverse magnetic fields.  For example, experimental searches for
non-zero permanent electric dipole moments, such as of the neutron and
the $^{199}$Hg and $^{225}$Ra atoms, require highly uniform magnetic
fields and past and future experiments
\cite{baker06,balashov07,altarev09,griffith09,holt10,perezgalvan11}
have typically used highly-optimized $\cos\theta$ or solenoidal coils.
In these experiments, the volume required by the $\cos\theta$ and
solenoidal cylindrical coil forms is, in general, significantly larger
than the sensitive experimental volume over which the requirements on
the field uniformity and/or field gradient must be satisfied.

\begin{figure}
\begin{center}
\includegraphics[scale=0.45]{cos_theta_coil_schematic_2.pdf}
\caption{Schematic diagram of a $\cos\theta$ coil with a continuous
  surface current distribution $\vec{K} = K \cos\theta\hat{z}$.
  Note that a coil with such a surface current distribution is also
  commonly referred to as a ``sine-phi'' coil \cite{jin98,bidinosti05},
  where the angle $\phi$ is the usual cylindrical coordinate.}
\label{fig:cos_theta_coil_schematic}
\end{center}
\end{figure}

Perhaps less well known or appreciated is the fact that a
\textit{finite-sized} spherical coil with a continuous surface current
$\vec{K} = K\sin\theta \hat{\phi}$, shown schematically in
Fig.\ \ref{fig:spherical_coil_schematic}, where $(\theta,\phi)$ are
the usual spherical coordinates, will generate a perfectly uniform
magnetic field along the sphere's $\hat{z}$-axis within the spherical
volume \cite{jacksonEM,smythe,haus_melcher,everett66,jin98}.  Although
the technical challenges associated with the fabrication of an
experimentally-realizable discretized spherical coil (e.g., precise
wire placement, etc.) may be more significant than those associated
with the fabrication of solenoidal and discretized $\cos\theta$ coils,
we note that the potential appeal of a spherical coil is that even for
the (non-realistic) case of continuous surface currents, a spherical
coil of finite size will generate a perfectly uniform field, whereas
the $\cos\theta$ and solenoidal coils must be of infinite length.
That is, a spherical coil does not suffer from any such
``finite-length'' effects'' that are by definition present in
finite-sized $\cos\theta$ and solenoidal coils.  Therefore, we
conjecture that the ratio of the fiducial volume (i.e., the volume
over which the requirements on the field uniformity or field gradient
must be satisfied) to the total volume occupied by a spherical coil
form will be larger than this ratio for $\cos\theta$ and solenoidal
coil forms.  And, indeed, a new experimental search for the neutron
electric dipole moment will employ a spherical coil \cite{masuda12}.
However, to our knowledge, detailed calculations of the off-axis field
uniformity properties of a discretized spherical coil have not been
presented in the literature.

\begin{figure}
\begin{center}
\includegraphics[scale=0.45]{spherical_coil_schematic_2.pdf}
\caption{Schematic diagram of a spherical coil with a continuous
surface current distribution $\vec{K}=K \sin\theta\hat{\phi}$.}
\label{fig:spherical_coil_schematic}
\end{center}
\end{figure}

In the remainder of this paper we investigate this conjecture.  We
begin, in Section \ref{sec:models}, by describing our models for our
calculations of the magnetic fields for these three coil types,
including our methods for the calculation of the off-axis fields for
the solenoidal and spherical coils.  We then present the numerical
results of our calculations in Section \ref{sec:results}, where we
ultimately compare the ratios of the fiducial volume to the total
volume occupied by each of the coil types for various requirements on
the field uniformity or field gradient.  Finally, we conclude with a
summary of our findings in Section \ref{sec:summary}.  We emphasize
that throughout this paper we are considering standard (i.e.,
non-optimized) versions of the $\cos\theta$ and solenoidal coils
(i.e., truncated versions of the infinitely-long ideal coils).

\section{Models for Discretized and Finite-Sized Coils}
\label{sec:models}

In this section we describe our numerical models for the calculation
of the magnetic fields (on- and off-axis) from discretized and
finite-sized $\cos\theta$, solenoidal, and spherical coils.  Note that
our discretization of the ideal surface currents is such that we use
a single point (i.e., zero radius) line wire to approximate regions
of equal integrated surface current.  We describe our implementation
of such for each of the coil types below.

\subsection{Model for $\cos\theta$ Coil}
\label{sec:models_cos_theta}

\begin{figure}
\begin{center}
\includegraphics[scale=0.50]{cos_theta_coil_discretized.pdf}
\caption{Schematic diagram of our model for a standard discretized and
  finite-sized $\cos\theta$ coil.  The thick lines indicate the
  rectangular current loops.  Note that this choice of a return path
  for the currents is not unique; another well-known approach is to
  employ a ``saddle-shaped'' winding \cite{bidinosti05,perezgalvan11},
  in which the return wires are wound around the perimeters of the
  cylinder ends.}
\label{fig:cos_theta_coil_discretized}
\end{center}
\end{figure}

A schematic diagram of our discretized and finite-length $\cos\theta$
coil is shown in Fig.\ \ref{fig:cos_theta_coil_discretized}.  As shown
in Ref.\ \cite{bidinosti05}, the discretization of an infinitely-long
continuous $\vec{K}\ = K \cos\theta\hat{z}$ surface current
distribution into line currents approximating the continuous surface
current distribution (obtained by apportioning the total integrated
current from $\theta = 0$ to $\pi/2$ into equal discrete parts) is
such that the line currents are spaced at equal intervals along the
indicated $x$-axis.  Thus, as indicated in
Fig.\ \ref{fig:cos_theta_coil_discretized}, our model consists of a
cylindrical shell of radius $a$ and length $L$, upon which are wound
$N$ rectangular current loops, spaced at equal intervals $d = 2a/N$
along the $x$-axis at positions $x = \pm \frac{1}{2}d$, $\pm
\frac{3}{2}d$, \ldots, $\pm \frac{N-1}{2}d$.  The line currents then
flow in the $+z$ ($-z$) direction for $\cos\theta>0$ $(<0)$ along the
length of the cylinder.  Our model then assumes a simple return path
for the currents (i.e., straight line currents along the $\pm y$
direction), such that our closed current loops are rectangular.
Thus, all of the
rectangular current loops are of the same length $L$ along the axis of
the cylindrical form (i.e., the $z$-axis), but are of different
lengths along the $y$-axis, where these lengths are constrained by the
cylindrical form to be $2\sqrt{a^2 - x_i^2}$, where $x_i$ denotes the
$i$th loop's position along the $x$-axis.

The calculations we present later have assumed $N$ even; one could, of
course, consider $N$ odd, with an additional wire located at
$\theta=0$, but this would not modify the symmetry of the problem.
Also, note that although our model does not include any connecting
wires between adjacent current loops (which are, of course, required
if a $\cos\theta$ is to be wound with a single wire), an actual
$\cos\theta$ coil can be wound in such a way so that the contributions
to the field from the connecting wires carrying the current from, say,
loop $i$ to loop $i+1$ to loop $i+2$, etc.\ are largely cancelled by
the return wire to the current source.

It is then straightforward to calculate the magnetic field everywhere
in space in rectangular coordinates via the Biot-Savart law.

\subsection{Model for Solenoidal Coil}
\label{sec:models_solenoidal}

Because the surface current distribution on an ideal solenoid is
uniform in the cylindrical coordinates $\phi$ and $z$ (with $z$
oriented the solenoid axis), such that $\vec{K} = K \hat{\phi}$,
our model for a discretized and
finite-length solenoidal coil consists of a cylindrical coil form of
radius $a$ and length $L$ upon which are wound $N$ circular current
loops of radius $a$ oriented in the indicated $xy$-plane and spaced at
equal intervals $d$ along the length of the cylinder.  Note that each
of our discretized current loops is ``isolated'' (i.e., stand-alone
circular loops, with no connecting wires), such that we do not attempt
to account for the effects of a helical winding inherent to the
winding of a solenoid with one continuous wire.

We then employ cylindrical coordinates to calculate the magnetic field
of our solenoidal coil at some observation point $\vec{x} =
(\rho,\phi,z)$.  By symmetry, the vector potential
$\vec{A}_i(\vec{x})$ of the $i$th circular current loop centered at
$(0,0,z_i)$ includes only a $\phi$-component, and is of the form
\cite{jacksonEM,smythe}
\begin{equation}
A_{\phi,i}(\rho,z) = \frac{\mu_0 I}{\pi \kappa_i} \sqrt{\frac{a}{\rho}}
\left[\left(1 - \frac{1}{2}\kappa_i^2\right)K(\kappa_i) - E(\kappa_i)\right],
\end{equation}
where $K(\kappa_i)$ and $E(\kappa_i)$ denote, respectively, the
complete elliptic integrals of the first and second kind, $\kappa_i^2
\equiv 4a\rho/[(a+\rho)^2 + (\Delta z_i)^2]$, and $\Delta z_i \equiv z
- z_i$.

The resulting $\rho$- and $z$-components of the $\vec{B}_i(\vec{x})$
field at $\vec{x}$ due to the $i$th circular current loop are then
calculated as $\vec{B} = \vec{\nabla} \times \vec{A}$ in cylindrical
coordinates and are of the form
\begin{eqnarray}
B_{\rho,i}(\vec{x}) &=& \frac{\mu_0 I (\Delta z_i)}
{2\pi\rho\sqrt{(a+\rho)^2 + (\Delta z_i)^2}}
\Bigg[ -K(\kappa_i)
+ \frac{a^2 + \rho^2 + (\Delta z_i)^2}
{(a-\rho)^2 + (\Delta z_i)^2}E(\kappa_i)\Bigg], \nonumber \\
\\
B_{z,i}(\vec{x}) &=& \frac{\mu_0 I}{2\pi\sqrt{(a+\rho)^2 + (\Delta z_i)^2}}
\Bigg[K(\kappa_i)
+ \frac{a^2 - \rho^2 - (\Delta z_i)^2}
{(a-\rho)^2 + (\Delta z_i)^2}E(\kappa_i)\Bigg], \nonumber
\end{eqnarray}
with the complete elliptic integrals then expressed as power series
\cite{arfken}
in the parameter $\kappa_i$,
\begin{eqnarray}
K(\kappa_i) &=& \frac{\pi}{2}\left\{1 + \sum_{n=1}^\infty
\left[\frac{(2n-1)!!}{(2n)!!}\right]^2 \kappa_i^{2n}\right\}, \nonumber \\
\\
E(\kappa_i) &=& \frac{\pi}{2}\left\{1 - \sum_{n=1}^\infty
\left[\frac{(2n-1)!!}{(2n)!!}\right]^2 \frac{\kappa_i^{2n}}{2n-1}\right\}.
\nonumber
\end{eqnarray}
Note that the values of $\kappa_i$ are restricted to the range $0 \leq
\kappa_i \leq 1$.  We sum the power series until the contribution of
the $n$th term is $< 10^{-9}$.  The resulting $B_{\rho,i}(\vec{x})$
and $B_{z,i}(\vec{x})$ field components for each circular current loop
are then converted to rectangular field components.  The total field
$\vec{B}(\vec{x})$ at observation point $\vec{x}$ from all $N$
circular current loops comprising the discretized solenoidal coil is
then simply the superposition of the rectangular field components from
each of the current loops.

\subsection{Model for Spherical Coil}
\label{sec:models_spherical}

\begin{figure}
\begin{center}
\includegraphics[scale=0.50]{spherical_coil_discretized.pdf}
\caption{Schematic diagram of our model for a standard discretized
spherical coil.  The thick lines indicate the circular current loops.}
\label{fig:spherical_coil_discretized}
\end{center}
\end{figure}

Finally, a schematic diagram of our discretized spherical coil is
shown in Fig.\ \ref{fig:spherical_coil_discretized}.  As indicated
there, our model consists of a spherical form of radius $R$ upon which
are wound $N$ circular current loops oriented in the $xy$-plane.  The
discretization of the spherical coil's $\vec{K} = K
\sin\theta\hat{\phi}$ continuous surface current distribution into $N$
discrete line wires proceeds in an identical manner to how the
discretization of a $\cos\theta$ coil current is performed in
Ref.\ \cite{bidinosti05}, and is as follows.  Suppose we desire to
divide the integral of the surface current, $\propto
\int\limits_0^{\pi} \sin\theta~R~d\theta \propto 2R$, into $N$ equal
parts, with each part bounded by the angles $[\alpha_{i-1},\alpha_i]$
for $i=1, 2, \ldots, N$, where $\alpha_0 = 0$ and $\alpha_N = \pi$.
We then approximate each of these parts carrying equal currents
$\propto 2R/N$ with a wire.  Then, let $\theta_i$ denote the angular
position of the $i$th wire,
where $\theta_i$ locates the midpoint
of the integrated surface current over the angular interval
$[\alpha_{i-1},\alpha_i]$.  Thus, per these definitions, $\theta_i$
must satisfy the condition
\begin{equation}
\int_{\alpha_{i-1}}^{\theta_i} \sin\theta~R~d\theta =
\int_{\theta_i}^{\alpha_i} \sin\theta~R~d\theta = \frac{1}{2}\frac{2R}{N}.
\end{equation}
From this, it follows that $\cos\alpha_{i-1} - \cos\alpha_i = 2/N$,
and so starting from $\alpha_0 = 0$, we find that
$\cos\alpha_i = 1 - 2i/N$.  From the above integrals, it also follows
that $\cos\theta_i = \cos\alpha_{i-1} - 1/N$, from which it then
follows that the $z$-coordinate for the $i$th discretized wire is
\begin{equation}
z_i = R\cos\theta_i = R\left(1 - \frac{2i-1}{N}\right).
\end{equation}
This shows that a discretized approximation to a spherical coil
consists of circular current loops spaced at equal intervals of $d =
2R/N$ along the $z$-axis at positions $z_i = \pm \frac{1}{2}d$, $\pm
\frac{3}{2}d$, \ldots, $\pm \frac{N-1}{2}d$.  From this, it then
follows that the radius $a_i$ of the $i$th current loop is constrained
to be $a_i = \sqrt{R^2 - z_i^2}$, and the radii of the two smallest
loops, at positions $z_1 = -z_N = \frac{N-1}{2}d$, are $a_1 = a_N =
\sqrt{R^2 - z_{1,N}^2}$.  Again, we note that the above is for $N$
even; however, as before, $N$ odd, with a wire at $\theta=\pi/2$,
would not modify the symmetry of the problem.

To calculate the magnetic field of our discretized spherical coil
everywhere inside of the coil, we now employ $(r,\theta,\phi)$
spherical coordinates, where the origin is located at the center of
our spherical coil.  As shown in Refs.\ \cite{smythe,ferraro}, the
vector potential $A_{\phi,i}(\vec{x})$ of the $i$th circular current
loop (located at $\theta = \theta_i$) at some observation point
$\vec{x} = (r<R, \theta, \phi)$, as derived via the stream function
approach, is of the form
\begin{equation}
A_{\phi,i}(\vec{x}) = \frac{\mu_0 I}{2} \sum_{n=1}^\infty
\frac{\sin\theta_i}{n(n+1)}\left(\frac{r}{R} \right)^n P^1_n(\cos\theta_i)
P^1_n(\cos\theta).
\end{equation}
From this, it then follows that the $r_i$- and
$\theta_i$-components of the magnetic field at $\vec{x}$, as
calculated from $\vec{B} = \vec{\nabla} \times \vec{A}$ in
spherical coordinates, are then for $r < R$ of the explicit form
\begin{equation}
\left.\begin{array}{c}
B_{r,i}(\vec{x}) \\
B_{\theta,i}(\vec{x})
\end{array}\right\} =
\frac{\mu_0 I \sin\theta_i}{2R} \sum_{n=1}^\infty
\left(\frac{r}{R} \right)^{n-1} P^1_n(\cos\theta_i)
\left\{\begin{array}{l}
P_{n}(\cos\theta) \\
\displaystyle{-\frac{1}{n}}P^1_n(\cos\theta),
\end{array}\right.
\end{equation}
where $P_n$ denotes the ordinary Legendre polynomial and $P_n^1$
the associated Legendre polynomial.
We then summed the field components from all of the current loops
comprising the spherical coil, and then converted the
resulting $B_r(\vec{x})$ and $B_\theta(\vec{x})$ field
components in spherical coordinates to rectangular field components,
ultimately for comparison with the other two coil types.

\section{Comparisons of Field Uniformities and Field Gradients}
\label{sec:results}

\begin{figure}
\begin{center}
\includegraphics[scale=0.62]{cos_theta_fractional_fiducial.pdf}
\caption{Fractional fiducial volume for $\beta < 10^{-3}$ (top plot)
  and $<10^{-4}$ (bottom plot) field uniformities for a standard
  discretized and finite-sized $\cos\theta$ coil, as a function of the
  number of rectangular loops and the coil's aspect ratio, or
  length-to-radius ratio $L/a$.}
\label{fig:cos_theta_fractional_fiducial}
\end{center}
\end{figure}

We now compare the field uniformities and gradients from standard
discretized and finite-sized $\cos\theta$, solenoidal, and spherical
coils for two example scenarios.  First, to illustrate the fiducial
volume properties of these coils, we show in
Figs.\ \ref{fig:cos_theta_fractional_fiducial}--\ref{fig:fractional_fiducial_spherical}
calculations of each coil's fractional fiducial volume (i.e., the
fraction of each coil's volume which satisfies the fiducial volume
requirement) for $10^{-3}$ and $10^{-4}$ field uniformities,
where we define the field uniformity $\beta$ at some point $\vec{x}$
relative to the center of the coil, $\vec{x}=0$, to be
\begin{equation}
\beta \equiv \left| \frac{|\vec{B}(\vec{x})| - |\vec{B}(\vec{x}=0)|}
{|\vec{B}(\vec{x}=0)|} \right|.
\end{equation}
The calculations for the $\cos\theta$ coil are shown in
Fig.\ \ref{fig:cos_theta_fractional_fiducial} for $\beta < 10^{-3}$
and $< 10^{-4}$ in two-dimensional parameter space of the number of
rectangular current loops, $N$, and the coil's length-to-radius ratio,
$L/a$.  Those for the solenoidal coil are shown in
Fig.\ \ref{fig:solenoidal_fractional_fiducial} for $\beta < 10^{-3}$
and $< 10^{-4}$ in two-dimensional parameter space in terms of the
number of circular loops, and the solenoid's length-to-radius ratio,
$L/a$.  Finally, those for the spherical coil are shown in
Fig.\ \ref{fig:fractional_fiducial_spherical} for $\beta <
10^{-3}$ and $< 10^{-4}$ as a function of the number of turns, which
is the only free parameter for a standard spherical coil of fixed
radius.

We note that for the $\cos\theta$ and solenoidal coils, further
increases in the number of turns $N$ beyond the upper range of $N=100$
shown in
Figs.\ \ref{fig:cos_theta_fractional_fiducial}--\ref{fig:solenoidal_fractional_fiducial}
does not lead to appreciable increases in their fractional fiducial
volumes.  Thus, our results clearly indicate that the fractional
fiducial volume of a spherical coil can be significantly larger (by up
to an order of magnitude for spherical coils consisting of $> 60$
circular loops) than those of $\cos\theta$ and solenoidal coils.  As a
visual illustration of the uniformity properties of a spherical coil,
Fig.\ \ref{fig:fractional_fiducial_spherical_growth} shows the
fiducial volumes of 1.0-m radius spherical coils consisting of $N=10$,
30, and 50 circular loops which satisfy a $\beta <10^{-3}$ and $< 10^{-4}$
field uniformity requirement.

\begin{figure}
\begin{center}
\includegraphics[scale=0.62]{solenoidal_fractional_fiducial.pdf}
\caption{Fractional fiducial volume for $\beta < 10^{-3}$ (top plot)
  and $\beta < 10^{-4}$ (bottom plot) field uniformities for a
  standard discretized and finite-sized solenoidal coil, as a function
  of the number of circular loops comprising the solenoid and the
  coil's aspect ratio, or length-to-radius ratio $L/a$.}
\label{fig:solenoidal_fractional_fiducial}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[scale=0.50]{fractional_fiducial_spherical.pdf}
\caption{Fractional fiducial volume for $\beta < 10^{-3}$ (open
  circles) and $\beta < 10^{-4}$ (filled circles) field uniformities
  for a standard discretized spherical coil, as a function of the
  number of circular loops comprising the spherical coil.}
\label{fig:fractional_fiducial_spherical}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[scale=0.68]{fractional_fiducial_spherical_growth.jpg}
\caption{Visualization of the fractional fiducial volumes of 1.0-m
  radius spherical coils consisting of $N=10$, 30, and 50 circular
  loops which satisfy $\beta < 10^{-3}$ (dark regions) and
  $<10^{-4}$ (gray regions) field uniformities.  Note that the circular
  current loops are in the $xy$ plane.}
\label{fig:fractional_fiducial_spherical_growth}
\end{center}
\end{figure}

\begin{table}[t!]
\begin{center}
\caption{Geometry parameters for spherical coils (number of turns $N$,
radius of sphere $R$, and the resulting radii of the two smallest
current loops $a_{\mathrm{min}} \equiv a_1 = a_N$; see text in Section
\ref{sec:models_spherical}) which would satisfy a requirement on the
fractional field gradient of $\gamma < 10^{-5}$ everywhere inside of a
$10 \times 10 \times 10$ cm$^3$ cube.  These are compared with
geometry parameters (number of turns $N$, radius of cylindrical coil
form $a$, and length of cylindrical coil form $L$) for benchmark
$\cos\theta$ and solenoidal coils.  The total volume $V$ occupied by
each coil form is also listed.}
\begin{tabular}{lcccc} \hline\hline
Coil& $N$& $R$ [cm]& $a_{\mathrm{min}}$ [cm]& $V$ [m$^3$] \\ \hline
Spherical& 20&  87& 27.2& 2.76 \\
Spherical& 40&  45& 10.0& 0.38 \\
Spherical& 60&  31& 5.64& 0.12 \\
Spherical& 80&  24& 3.78& 0.06 \\
Spherical& 100& 20& 2.82& 0.03 \\ \hline
Coil& $N$& $a$ [cm]& $L$ [cm]& $V$ [m$^3$] \\ \hline
$\cos\theta$& 100& 19& 229& 0.26 \\
Solenoidal& 100& 13& 249& 0.13 \\ \hline
\end{tabular}
\label{tab:examples_geometry_parameters_cube}
\end{center}
\end{table}

As a second example, we consider an experimental requirement on the
fractional field gradient, $\gamma$, which we define to be
\begin{equation}
\gamma \equiv \frac{1}{|\vec{B}(\vec{x}=0)|}
\left|\frac{\partial B_i}{\partial x_i}\right|,
\end{equation}
where $B_i$ denotes the field component along the field's primary
direction (i.e., the $x$-component for the $\cos\theta$ coil and the
$z$-component for the solenoidal and spherical coils).  To illustrate,
suppose a hypothetical requirement is that $\gamma < 10^{-5}$
everywhere inside of a $10 \times 10 \times 10$ cm$^3$ rectangular
volume located at the center of each coil.  In Table
\ref{tab:examples_geometry_parameters_cube} we show examples of
spherical coil geometry parameters which would satisfy this
requirement, and compare these with benchmark examples of geometry
parameters for approximate best-case $\cos\theta$ and solenoidal coils
(i.e., large length-to-radius ratios and number of turns, as suggested
by the fractional fiducial volume calculations presented in
Figs.\ \ref{fig:cos_theta_fractional_fiducial}--\ref{fig:solenoidal_fractional_fiducial}).
Our results indicate that a particular requirement on the fractional
field gradient could be achieved with a spherical coil occupying a
smaller volume than the volumes occupied by $\cos\theta$ and
solenoidal coils.

However, we do note that a potential appeal of the $\cos\theta$-coil
(if wound with a ``saddle winding'' \cite{bidinosti05,perezgalvan11})
and solenoidal-coil geometries is that their winding patterns permit
significant access to their interior regions along their respective
$z$-axes (i.e., via a circular aperture of radius equal to the radius
$a$ of their cylindrical form), whereas access to the interior of a
spherical coil along its $z$-axis is constrained by the radii of the
two smallest current loops, $a_1 = a_N = \sqrt{R^2 - z_{1,N}^2}$,
located at positions $z_1 = -z_N = \frac{N-1}{2}d$, as discussed
previously in Section \ref{sec:models_spherical}.  If a $\cos\theta$
coil is wound with rectangular loops (as in our model calculations),
access to the interior along its $z$-axis is restricted by the
spacing $d=2a/N$ between adjacent rectangular loops (see
Fig.\ \ref{fig:cos_theta_coil_discretized}).

For access to their interior regions along the transverse direction,
the dimensions of the $\cos\theta$ coil's two smallest rectangular loops
are $\frac{2 a}{N} \sqrt{2 N - 1}$ ($y$-direction) $\times$ $L$
($z$-direction; see Fig.\ \ref{fig:cos_theta_coil_discretized}).  The
solenoidal coil can, in principle, be accessed in the transverse
direction via the gaps of size $\frac{L}{N-1}$ in between adjacent
circular windings, and the spherical coil can also be accessed in the
transverse direction in between the circular current loops, which are
spaced at equal intervals of $d=2R/N$ along its $z$-axis (see
Fig.\ \ref{fig:spherical_coil_discretized}).

\section{Summary}
\label{sec:summary}

In this paper we have presented a model for numerical calculations of
the magnetic field of a standard discretized spherical coil.  Because
a discretized spherical coil does not suffer from any ``finite-length
effects'' that are inherent to standard discretized and finite-sized
$\cos\theta$ and solenoidal coils, our calculations indicate that the
fractional fiducial volume of a discretized spherical coil is
potentially significantly larger than those of standard (i.e.,
non-optimized) discretized and finite-length $\cos\theta$ and
solenoidal coils. \\

\noindent\textbf{Acknowledgments} \\ \\
\indent This work was supported in part by the U.\ S.\ Department of
Energy Office of Nuclear Physics under Award Number DE-FG02-08ER41557.
We thank Prof.\ J.\ W.\ Martin for encouraging us to write this paper,
and Prof. C.\ B.\ Crawford for several valuable discussions.  We also
thank an anonymous referee for many excellent suggestions which, we
believe, significantly improved this paper.


\section{Introduction}
SBs are cavities across the interstellar medium (ISM) created by cataclysmic processes such as strong stellar winds and supernova explosions associated with massive stars in an OB association. The swept up material by stellar winds and supernova remnants produce an outer dense shell with an interior filled by low-density hot gas. The gas at the center of the SB is bright in X-rays, the outer shell formed by dust and cool gas is observed in the infrared while the interface between the cool shell and the hot interior is seen in the ultraviolet (UV). O~{\small VI} absorption lines at 1032 \AA\ and 1037.6 \AA\ are diagnostic of the energetic processes of the interface environments of the SBs which are produced by shock heating and are collisionally ionized. The LMC hosts more than 20 SBs which have been observed in X-rays and UV due to their proximity ($\sim$50 kpc). {\em FUSE} has observed O~{\small VI} doublet along many sightlines using early type stars as background objects. We have presented a study of O~{\small VI} absorption at 1032 \AA\ in the SBs of the LMC in order to assess the properties of intermediate environments in them.

\footnotesize{
\begin{table}
\caption{ Details of O~{\small VI} observations in the SBs of the LMC.}
\label{tab1}
\begin{tabular}{c c c c c c c c}
\hline
SBs      & Size of SBs              & {\em FUSE} Targets   & RA (J2000)   & Dec (J2000)    & FWHM       & Integration & N(O~{\small VI})/ \\
        &                         & 	     	   & hr min sec   & deg min sec    & (m\AA)     & limit(km/s) &  (10$^{14}$ atoms/cm$^{2}$) \\
\hline
N11     & \(10'\)$\times$\(7'\)    & PGMW-3070     & 04 56 43.25  &   -66 25 02.0  & 132$\pm$23 &  180, 345   & 1.27$^{+0.43}_{-0.43}$ \\
        &                          & LH103102      & 04 56 45.40  &   -66 24 45.9  & 132$\pm$23 &  180, 330   & 1.46$^{+0.10}_{-0.13}$ \\
        &                          & LH91486       & 04 56 55.58  &   -66 28 58.0  & 266$\pm$29 &  175, 385   & 2.96$^{+0.68}_{-0.68}$ \\
        &                          & PGMW-3223     & 04 57 00.80  &   -66 24 25.3  & 129$\pm$15 &  175, 315   & 1.27$^{+0.18}_{-0.18}$ \\
N51     & \(12'\)$\times$\(10'\)   & Sk-67D106     & 05 26 15.20  &   -67 29 58.3  & 180$\pm$7  &  175, 345   & 2.01$^{+0.50}_{-0.53}$ \\
        &                          & Sk-67D107     & 05 26 20.67  &   -67 29 55.4  & 254$\pm$10 &  160, 360   & 2.85$^{+0.27}_{-0.27}$\\
        &                          & Sk-67D111     & 05 26 47.95  &   -67 29 29.9  & 214$\pm$19 &  175, 365   & 2.20$^{+0.29}_{-0.24}$\\
N57     & \(12'\)$\times$\(7'\)    & Sk-67D166     & 05 31 44.31  &   -67 38 00.6  & 206$\pm$9  &  165, 390   & 2.09$^{+0.25}_{-0.20}$ \\
N70     & \(8'\)$\times$\(7'\)     & SK-67D250     & 05 43 15.48  &   -67 51 09.6  & 316$\pm$33 &  165, 375   & 3.71$^{+0.23}_{-0.23}$ \\
        &                          & D301-NW8      & 05 43 15.96  &   -67 49 51.0  & 228$\pm$30 &  175, 365   & 2.60$^{+0.09}_{-0.13}$ \\
        &                          & D301-1005     & 05 43 08.30  &   -67 50 52.4  & 284$\pm$57 &  165, 385   & 3.37$^{+0.17}_{-0.22}$ \\
N144    & \(13'\)$\times$\(12'\)   & HD36521       & 05 26 30.32  &   -68 50 25.4  & 126$\pm$9  &  175, 340   & 1.18$^{+0.28}_{-0.28}$ \\
        &                          & Sk-68D80      & 05 26 30.43  &   -68 50 26.6  & 303$\pm$16 &  145, 335   & 3.61$^{+0.30}_{-0.26}$\\
N204    & \(14'\)$\times$\(13'\)   & Sk-70D91      & 05 27 33.74  &   -70 36 48.3  & 256$\pm$7  &  160, 365   & 2.69$^{+0.24}_{-0.17}$ \\
N206    & \(9'\)$\times$\(15'\)    & BI184         & 05 30 30.60  &   -71 02 31.3  & 118$\pm$10 &  165, 330   & 1.09$^{+0.22}_{-0.22}$ \\
        &                          &  Sk-71D45     & 05 31 15.55  &   -71 04 08.9  & 194$\pm$9  &  160, 345   & 1.80$^{+0.24}_{-0.21}$ \\
N154    & \(12'\)$\times$\(8'\)    & SK-69D191     & 05 34 19.39  &   -69 45 10.0  & 185$\pm$25 &  165, 340   & 1.65$^{+0.26}_{-0.23}$\\
N158    & \(8'\)$\times$\(7'\)     & HDE269927     & 05 38 58.25  &   -69 29 19.1  & 245$\pm$8  &  160, 320   & 2.62$^{+0.16}_{-0.21}$ \\
30DOR C & \(7'\)$\times$\(6'\)     & MK42          & 05 38 42.10  &   -69 05 54.7  & 228$\pm$24 &  160, 330   & 2.60$^{+0.46}_{-0.45}$ \\
        &                          & SK-69D243     & 05 38 42.57  &   -69 06 03.2  & 307$\pm$15 &  150, 345   & 3.63$^{+0.40}_{-0.45}$ \\
        &                          & 30DOR-S-R136  & 05 38 51.70  &   -69 06 00.0  & 185$\pm$25 &  165, 320   & 1.77$^{+0.20}_{-0.24}$  \\
        &                          & SK-69D246     & 05 38 53.50  &   -69 02 00.7  & 211$\pm$7  &  155, 325   & 2.36$^{+0.18}_{-0.14}$ \\

\hline
\end{tabular}
\end{table}
}

\section{Measurement of O~{\small VI} Column Density}
We selected high resolution {\em FUSE} spectra from the archival data which were observed in the SB lines of sight in the LMC and had well pronounced O~{\small VI} profiles. We found 22 good O~{\small VI} observations covering 10 SBs in the LMC. The stellar continuum of O~{\small VI} profiles are of low order Legendre polynomials ($\le$ 5). The measurement of the optical depth of O~{\small VI} and the equivalent width for the LMC component are done by the apparent optical depth method following the procedures of \cite[Savage \& Sembach (1991)]{Savage91}, \cite[Sembach \& Savage (1992)]{Sembach92} and  \cite[Howk et al. (2002)]{Howk02}. The results are listed in Table \ref{tab1}. The 1$\sigma$ error in the measurements were estimated from the fitting procedures using the uncertainties in the {\em FUSE} data. The Milky Way ($v\lesssim$ 150 km/s) and LMC ($v\gtrsim$ 150 km/s) absorption components are well resolved by {\em FUSE} for most of the sightlines. The average value of the equivalent width of O~{\small VI} profiles of the SBs in the LMC comes out to be 214$\pm$19 m\AA.

\section{Results and Discussions}
N(O~{\small VI}) of SBs varies from (1.09$\pm$0.22)$\times$10$^{14}$ to (3.71$\pm$0.23)$\times$10$^{14}$ atoms/cm$^{2}$. The mean N(O~{\small VI}) for the SBs is found to be (3.07$\pm$0.62)$\times$10$^{14}$ atoms/cm$^{2}$ while the mean N(O~{\small VI}) for the non-SBs lines of sight is (2.10$\pm$0.50)$\times$10$^{14}$ atoms/cm$^{2}$ (\cite[Pathak et al. 2011]{Pathak11}). Thus, the SBs in the LMC show an excess O~{\small VI} abundance of about 46\% in comparison to non-SB regions. Studies for SB N70 (\cite[Danforth et al. 2006]{Danforth06}) found similar results with 60\% more O~{\small VI} than the non-SB targets. Considering N(O~{\small VI}) inside the individual SBs (e.g., N11 and 30 Dor C), a variation of about a factor of 2 to 2.5 in N(O~{\small VI}) was obtained. The thermal conduction between the interior hot, X-ray producing gas and the cool, photoionized shell of SBs is found to be the most favourable mechanism for the production of O~{\small VI} in the SBs (\cite[Danforth et al. 2006]{Danforth06}) as the thermal conduction models (\cite[Weaver et al. 1977]{Weaver77} and \cite[Borkowski et al. 1990]{Borkowski90}) have found N(O~{\small VI}) $\approx$ few $\times$ 10$^{13}$ atoms/cm$^{2}$, which is a close approximation of the observed N(O~{\small VI}). O~{\small VI} in the SBs does not show any correlation with the ISM morphologies such as H$\alpha$ and X-ray surface brightnesses. The temperature estimated from the line width of O~{\small VI} comes to be $\sim$10$^{6}$ K which represents slightly higher FWHM than expected for O~{\small VI} absorption.


\section{Introduction}

\subsection{} Let $F$ be an imaginary quadratic field of discriminant
$D<0$, let $\ensuremath{\mathcal{O}} = \ensuremath{\mathcal{O}}_{D}$ be its ring of integers, and let $\Gamma$ be
the group $\GL_{N} (\ensuremath{\mathcal{O}})$.  The homology and cohomology of $\Gamma$
when $N=2$ --- or rather
its
close cousin the Bianchi group
$\PSL_{2} (\ensuremath{\mathcal{O}})$ --- have been well studied in the literature.  For an
(incomplete) selection of results we refer to
\cite{rahm1,rahm2,berkove,schwermer-vogtmann,vogtmann,sengun,cremona}.
Today we have a good understanding of a wide range of examples, and
one can even compute them for very large discriminants
(cf.~\cite{yasaki}).  For $N>2$, on the other hand, the group $\Gamma$
has not received the same attention.  The first example known to us is
the work of Staffeldt \cite{k3gauss}.  He treated the case $N=3$,
$D=-4$ with the goal of understanding the $3$-torsion in $K_{3}
(\ensuremath{\mathbb{Z}}[\sqrt{-1}])$.  The second example is \cite{coul-herm}, which
investigates the case of the groups $\GL(L)$ for $L$ an $\ensuremath{\mathcal{O}}$-lattice
that is not necessarily a free $\ensuremath{\mathcal{O}}$-module.  This allows the authors
to compute the Hermite constants of those rings in case $D\geq -10$
and $\rank(L) \leq 3$. Our methods apply as well to the non-free case, and
the corresponding cohomology computations would be useful when
investigating automorphic forms over number fields that are not
principal ideal domains (cf.~\cite[Appendix]{SteinBook} for more about
the connection between cohomology of arithmetic groups and automorphic
forms).

In this paper we rectify this situation somewhat by beginning the
first systematic computations for higher rank linear groups over
$\ensuremath{\mathcal{O}}$.  In particular investigate the cohomology of $\Gamma$ for
$N=3,4$ and for a selection of discriminants: $D\geq -24$ when $N=3$,
and $D=-3,-4$ when $N=4$.  We explicitly compute the polyhedral
reduction domains arising from Voronoi's theory of perfect forms, as
generalized by Ash \cite[Ch.~II]{AMRT_SmoothCompactification} and
Koecher \cite{koecher}.  This allows us to compute the integral
cohomology of $\Gamma$ up to $p$-power torsion for small primes $p$.
In a sequel \cite{aim-ktheory} to this paper, we will apply some of
these results to computations with the $K$-groups $K_{4} (\ensuremath{\mathcal{O}}_{D})$,
when $D=-3,-4$.

\subsection{} Here is a guide to the paper (which closely follows the
structure of the first five sections of \cite{PerfFormModGrp}).  In
Section~\ref{sec:polyhedralcone} we recall the explicit reduction
theory we need to build our chain complexes to compute cohomology.  In
Section \ref{sec:cohomologyandhomology} we define the complexes, and
explain the relation between what we compute and the cohomology of
$\Gamma$.  In Section \ref{sec:mass} we describe a ``mass formula''
for the cells in our tessellations that provides a strong
computational check on the correctness of our constructions.  In
Section~\ref{sec:explicit} we give an explicit representative for the
nontrivial class in the top cohomological degree; this construction is
motivated by a similar construction in \cite{PerfFormModGrp,
PEV-HG-CS-announce}.  Finally, in Section~\ref{sec:tables} we give the
results of our computations.

\subsection{Acknowledgments} We thank A.~Ash, P.~Elbaz-Vincent, and
C.~Soul\'e for helpful discussions.  This research, and the research
in the companion papers \cite{aim-voronoi, aim-ktheory} was conducted
as part of a ``SQuaRE'' (Structured Quartet Research Ensemble) at the
American Institute of Mathematics in Palo Alto, California in February
2012.  It is a pleasure to thank AIM and its staff for their support,
without which our collaboration would not have been possible.  Also
our initial computations of cohomology and $K$-groups occurred on the
large shared computers \texttt{parsley} and \texttt{rosemary} at the
University of Georgia mathematics department, and we thank them for
their support in allowing these machines to be used for collaborative
mathematical research projects.

\section{The polyhedral cone}\label{sec:polyhedralcone}
Fix an imaginary quadratic field $F$ of discriminant $D<0$ with ring of
integers $\ensuremath{\mathcal{O}} = \ensuremath{\mathcal{O}}_D$
and define the element
\[\omega :=
\omega_D :=
\begin{cases}
  \sqrt{D/4} & \text{if $D \equiv 0 \bmod{4}$,}\\
  (1 + \sqrt{D})/2 & \text{if $D \equiv 1 \bmod{4}$,}
\end{cases}
\]
so that
$F = \ensuremath{\mathbb{Q}}(\omega)$ and $\ensuremath{\mathcal{O}} = \ensuremath{\mathbb{Z}}[\omega]$.  Throughout we fix a complex
embedding $F \hookrightarrow \ensuremath{\mathbb{C}}$ and consistently identify $F$ with its image in
$\ensuremath{\mathbb{C}}$.  We also extend this identification to vectors and matrices with
coefficients in $F$.

\subsection{Hermitian forms} Let
$\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ denote the $N^2$-dimensional real vector space of $N
\times N$ Hermitian matrices with complex coefficients.  Using the
chosen complex embedding of $F$ we can view $\ensuremath{\mathcal{H}}^N(F)$, the Hermitian
matrices with coefficients in $F$, as a subset of $\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$.
Moreover this embedding allows us to view $\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ as a
$\ensuremath{\mathbb{Q}}$-vector space such that the rational points of $\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ are
exactly $\ensuremath{\mathcal{H}}^N(F)$.

Define a map
$q \colon \ensuremath{\mathcal{O}}^N \to \ensuremath{\mathcal{H}}^N(F)$
by the outer product $q(x) = x {x}^{*}$,
where  ${*}$ denotes conjugate transpose (with conjugation being
the nontrivial complex conjugation automorphism of $F$).
Each $A \in \ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ defines a Hermitian
form $A[x]$ on $\ensuremath{\mathbb{C}}^N$ by the rule
\begin{equation}
A[x] := x^* A x, \quad \text{for $x \in \ensuremath{\mathbb{C}}^N$}.
\end{equation}
Define the
non-degenerate bilinear pairing
\[\ip{\cdot, \cdot} \colon \ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}}) \times \ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}}) \to \ensuremath{\mathbb{C}}\]
by $\ip{A,B} := \Tr(AB)$.  For $x \in \ensuremath{\mathcal{O}}^N$ (identified with its image
in $\ensuremath{\mathbb{C}}^N$) one can easily verify that
\begin{equation}\label{eq:traceform}
A[x] = \Tr(Aq(x)) = \ip{A, q(x)}.
\end{equation}

Let $C_N \subset \ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ denote the cone of positive definite
Hermitian matrices.

\begin{definition}
  For $A \in C_N$, we define the \emph{minimum of $A$} as
\[m(A) := m_D(A) := \inf_{x \in \ensuremath{\mathcal{O}}^N \setminus \set{0}} A[x].\]
Note that $m (A)>0$ since $A$ is positive definite.  A vector $v \in
\ensuremath{\mathcal{O}}^N$ is called a \emph{minimal vector of $A$} if $A[v] = m(A)$.
We denote the
set of minimal vectors of $A$ by $M(A)$.
\end{definition}

It should be emphasized that these notions depend on the fixed choice
of the imaginary quadratic field $F$.  Since $q(\ensuremath{\mathcal{O}}^N)$ is discrete in
$\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$
and the level sets $q(x) = C$ are compact,
the minimum for each $A$ is attained by only finitely many
minimal vectors.

From \eqref{eq:traceform}, we see that each vector $v \in \ensuremath{\mathcal{O}}^N$
gives rise to a linear functional on $\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ defined by $q(v)$.
\begin{definition}\label{def:perfect}
We say
a Hermitian form $A \in C_N$ is a \emph{perfect Hermitian
form over $F$} if
\[\Span_\ensuremath{\mathbb{R}}\set{q(v) \mid v \in M(A)} = \ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}}).\]
\end{definition}

From Definition \ref{def:perfect} it is clear that a form is perfect
if and only if it is uniquely determined by its minimum and its
minimal vectors.  Equivalently, a form $A$ is perfect when $M(A)$
determines $A$ up to a positive real scalar.  It is convenient to
normalize a perfect form by requiring that $m (A)=1$, and we will do
so throughout this paper.  A priori there is no reason to expect that
every perfect form (up to rescaling) can be realized as
an $F$-rational point in $\ensuremath{\mathcal{H}}^{n} (\ensuremath{\mathbb{C}})$; indeed,
when one generalizes these concepts to general number fields this is
too much to expect (cf.~\cite{CubicField}).  However
for imaginary quadratic fields this rationality property does hold:

\begin{theorem}\rftheoremcite{Theorem~3.2}{okuda-yano}
Suppose $F$ is a CM field.  Then if $A \in C_N$ is a perfect Hermitian
form over $F$ such that $m(A) = 1$, we have $A \in \ensuremath{\mathcal{H}}^N(F)$.
\end{theorem}

We now construct a partial compactification of the cone $C_{N}$:
\begin{definition}
A matrix $A \in \ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ is said to have an \emph{$F$-rational
kernel} when the kernel of $A$ is spanned by vectors in $F^N \subset
\ensuremath{\mathbb{C}}^N$.  Let $C^*_N \subset \ensuremath{\mathcal{H}}^N (\ensuremath{\mathbb{C}})$ denote the subset of nonzero
positive semi-definite Hermitian forms with $F$-rational kernel.
\end{definition}

Let $\ensuremath{\mathbf{G}}$ be the reductive group over $\ensuremath{\mathbb{Q}}$ given by the restriction
of scalars $\Res_{F/\ensuremath{\mathbb{Q}}}(\GL_N)$.  Thus $\ensuremath{\mathbf{G}}(\ensuremath{\mathbb{Q}}) = \GL_N(F)$ and
$\ensuremath{\mathbf{G}}(\ensuremath{\mathbb{Z}}) = \GL_N(\ensuremath{\mathcal{O}})$.  The group $\ensuremath{\mathbf{G}}(\ensuremath{\mathbb{R}}) = \GL_N(\ensuremath{\mathbb{C}})$ acts on
$C^*_N$ on the left by
\[g \cdot A = gAg^*,\]
where $g \in \GL_N(\ensuremath{\mathbb{C}})$ and $A \in C^*_N$; one can easily verify that
this action preserves $C_N$.  Let $H = R_d\ensuremath{\mathbf{G}}(\ensuremath{\mathbb{R}})^0$ be the identity
component of the group of real points of the split radical of $\ensuremath{\mathbf{G}}$.
Then $H \simeq \ensuremath{\mathbb{R}}_+$, and as a subgroup of $\ensuremath{\mathbf{G}} (\ensuremath{\mathbb{R}})$ acts on
$C^*_N$ by positive real homotheties.  Voronoi's work \cite{VoronoiI},
generalized by Ash \cite[Ch.~II]{AMRT_SmoothCompactification} and
Koecher \cite{koecher} shows that there are only finitely many perfect
Hermitian forms over $F$ modulo the action of $\GL_N(\ensuremath{\mathcal{O}})$ and $H$.

Let $X^*_N$ denote the quotient $X^*_N = H\backslash C_N^*$, and let
$\pi\colon C^*_N \to X^*_N$ denote the projection.  Then $X_N =
\pi(C_N)$ can be identified with the global Riemannian symmetric space
for the reductive group $H\backslash\GL_N(\ensuremath{\mathbb{C}})$.
\begin{remark}\label{rem:orientable}
  The action of $\ensuremath{\mathbf{G}} (\ensuremath{\mathbb{R}}) = \GL_{N} (\ensuremath{\mathbb{C}})$ on $\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ gives a
representation $\rho\colon \GL_N(\ensuremath{\mathbb{C}}) \to \GL_{N^2}(\ensuremath{\mathbb{R}})$.  Any
element $g$ acts on the orientation of $\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ via the sign of
$\det(\rho(g))$.  Since $\GL_N(\ensuremath{\mathbb{C}})$ is connected, $\GL_N(\ensuremath{\mathbb{C}} )$ and
its subgroup $\GL_{N} (F)$ act on $\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ via orientation
preserving automorphisms.
\end{remark}

\subsection{Two cell complexes}\label{subsec:cell-complexes}
Let $M$ be a finite subset of $\ensuremath{\mathcal{O}}^N \setminus \set{0}$.  The
\emph{perfect cone} of $M$ is the set of nonzero matrices of the form
$\sum_{v \in M} \lambda_v q(v)$, where $\lambda_v\in \ensuremath{\mathbb{R}}_{\geq 0}$; by
abuse of language we also call its image by $\pi$ in $X^*_N$ a perfect
cone.  For a perfect form $A$, let $\sigma(A) \subset X^*_N$ be the
perfect cone of $M(A)$.  One can show
\cite{AMRT_SmoothCompactification, koecher} that the cells $\sigma(A)$
and their intersections, as $A$ runs over equivalence classes of
perfect forms, define a $\GL_{N} (\ensuremath{\mathcal{O}})$-invariant cell decomposition
of $X^*_N$.  In particular, for a perfect form $A \in C_N$ and an
element $\gamma \in \GL_N(\ensuremath{\mathcal{O}})$, we have
\[\gamma \cdot \sigma(A) = \text{perfect cone of $\set{\gamma v
  \mid v \in M(A)}$} = \sigma((\gamma ^*)^{-1}A \gamma ).\]
Endow $X^*_N$ with the CW-topology \cite[Appendix]{hatcher}.

If $\tau$ is a closed cell in $X^*_N$ and $A$ is a perfect form with
$\tau \subset \sigma(A)$, we let $M(\tau)$ denote the set of vectors
$v \in M(A)$ such that $q(v) \in \tau$.  The set $M (\tau)$ is
independent of the (possible) choice of $A$.  Then $\tau$ is the image
in $X^{*}_{N}$ of the cone $C_{\tau}$ generated by $\{q (v)\mid v \in
M(\tau) \}$.  For any two closed cells $\tau$ and $\tau'$ in $X^*_N$,
we have $M(\tau) \cap M(\tau') = M(\tau \cap \tau')$.

Let $\tilde{\Sigma} \subset C^*_N$ be the (infinite) union of all
cones $C_\sigma$ such that $\sigma \in X_N^*$ has nontrivial
intersection with $X_{N}$.  One can verify that the stabilizer of
$C_\sigma$ in $\GL_{N} (\ensuremath{\mathcal{O}})$ is equal to the stabilizer of $\sigma$.
By abuse of notation, we write $M(C_\sigma)$ as $M(\sigma)$.

The collection of cones $\tilde{\Sigma }$ was used by Ash \cite{ash77,
ash-small} to construct the \emph{well-rounded retract}, a
contractible $N^2 - N$ dimensional cell complex $W$ on which
$\GL_N(\ensuremath{\mathcal{O}})$ acts cellularly with finite stabilizers of cells.  More
precisely, a nonzero finite set $M \subset
 \ensuremath{\mathcal{O}}^{N}$ is called \emph{well-rounded} if the $\ensuremath{\mathbb{C}}$-span
of $M$ is $\ensuremath{\mathbb{C}}^N$.  For a well-rounded subset $M$, let $\sigma(M)$
denote the set of forms $A \in C_N$ with $M(A) = M$ and $m(A) = 1$.
It is easy to prove that if $\sigma(M)$ is non-empty then it is convex
and thus topologically a cell.  The well-rounded retract is then
defined to be
\[W = \bigcup_{\text{$M$ well-rounded}} \sigma(M).\]

The space $W$ is dual to the decomposition of $X_N$ in a sense made
precise in Theorem~\ref{theorem:duality} below.  For instance when $N = 2$
and $F = \ensuremath{\mathbb{Q}}$, $X$ can be identified with the upper
half-plane $\ensuremath{\mathfrak{H}} = \{x+iy\mid y>0 \}$.  The Voronoi tessellation is the tiling of $\ensuremath{\mathfrak{H}}$ by
with the $\SL_{2} (\ensuremath{\mathbb{Z}})$-translates of the ideal geodesic triangle
with vertices $\{0,1,\infty \}$.  The well-rounded retract is the dual
infinite trivalent tree (see Figure~\ref{fig:tess-spine}).

\begin{figure}
\includegraphics[scale=0.6]{tess-spine.eps}
  \caption{Voronoi tessellation of $\ensuremath{\mathfrak{H}}$ shown in green, with its dual the
    well-rounded retract shown in black.} \label{fig:tess-spine}
\end{figure}

\begin{lemma}\label{lem:wellrounded-positivedefinite}
If $C_\sigma \in \tilde{\Sigma}$, then $M(\sigma)$ is well-rounded.
\end{lemma}
\begin{proof}
This is proved by McConnell in \cite[Theorem~2.11]{mcconnell-projgeom}
when $F = \ensuremath{\mathbb{Q}}$, and we mimic his proof to yield the analogous result
for imaginary quadratic fields.  Let $M = M(C_\sigma)$ denote the
spanning vectors of $C_\sigma$.  Suppose $M$ is not well-rounded.
Then $M$ does not span $\ensuremath{\mathbb{C}}^N$, and so there exists a non-zero vector
$w \in \ensuremath{\mathbb{C}}^N$ such that $v w^* = 0$ for all $v \in M$.  Then each
$\theta \in C_\sigma$ can be written as a non-zero Hermitian form
$\theta = \sum_{v \in M}a_vq(v)$, where $a_v \geq 0$.  It follows that
\[\theta[w] = \sum_{v \in M}a_v \ip{q(v), q(w)} = \sum_{v \in M}a_v
|vw^*|^2 = 0,
 \]
where $|\cdot|$ is the usual norm on $\ensuremath{\mathbb{C}}^N$.  In particular, $\theta$
is not positive definite, contradicting the
assumption that $C_\sigma \in \tilde{\Sigma}$.
\end{proof}

\begin{theorem}\label{theorem:duality}
Let $M\subset \ensuremath{\mathcal{O}}^{N}$ be well-rounded and let $\sigma(M) \in W$.
Then there is a unique cone $C_\sigma \in \tilde{\Sigma}$ such that
$M(C_\sigma) = M$.  The map $\sigma(M) \mapsto C_\sigma$ is a
canonical bijection $W \to \tilde{\Sigma}$ that is
inclusion-reversing on the face relations.
\end{theorem}

\begin{proof}
It suffices to prove that a well-rounded subset $M \subset\ensuremath{\mathcal{O}}^N$ has
$\sigma(M) \neq \emptyset$ if and only if there is a cone $C_\sigma
\in \tilde{\Sigma}$ such that $M(\sigma) = M$.

Suppose $C_\sigma \in \tilde{\Sigma}$.  By
Lemma~\ref{lem:wellrounded-positivedefinite}, $M = M(C_\sigma)$ is
well-rounded.  Thus it remains to show that $\sigma(M)$ is non-empty.
We do so by constructing a form $B \in \sigma(M)$ with $M(B) =
M(C_\sigma)$.  Since $C_\sigma \in \tilde{\Sigma}$, there is a perfect
form $A$ such that $C_\sigma$ is a face of the cone $S_A =
\pi^{-1}(\sigma(A))$.  Furthermore, $C_\sigma$ can be described as the
intersection of $S_A$ with some supporting hyperplane $\set{\theta
\mid \ip{H,\theta} = 0}$ for some $H \in \ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$.  It follows that
\[
\ip{H,q(v)} = 0  \quad \text{if $v \in M$} \quad \text{and} \quad
\ip{H,q(v)} > 0  \quad \text{if $v \in M(A) \setminus M$.}
\]
Let $B = A + \rho H$.  Since $\ip{H, q(v)} = B[v]$, a standard
argument for Hermitian forms shows that for sufficiently small
positive $\rho$, $B$ is positive definite, $B[v] = m(A)$ for $v \in M$
and $B[v] > m(A)$ for $v \in M(A) \setminus M$.  Thus $M(B) = M$ and
so $B \in \sigma (M)$.  In particular, $\sigma(M)$ is non-empty.

Conversely, suppose $M$ is a well-rounded subset of $\ensuremath{\mathcal{O}}^N$ with
$\sigma(M)$ non-empty.  Choose $B \in \sigma(M)$ so that $M(B) = M$.
If $B$ is perfect, then $M = M(B)$ and we are done.  Otherwise, we can
use the generalization of an algorithm of Voronoi \cite{aim-voronoi}
to find a perfect form $A$ such that $M(B) \subset M(A)$. Let $H = B -
A$.  Then
\[
\ip{H,q(v)} = 0  \quad \text{if $v \in M$} \quad \text{and} \quad
\ip{H,q(v)} > 0  \quad \text{if $v \in M(A) \setminus M$.}
\]
Thus the hyperplane $\set{\theta \mid \ip{H,\theta} = 0}$ is a
supporting hyperplane for the subset of $S_A = \pi^{-1}(\sigma(A))$
spanned by $\set{q(v) \mid v \in M}$.  Therefore $M$ defines a face of
$S_A$ in $\tilde{\Sigma}$ as desired.
\end{proof}

\begin{remark}\label{rem:W-Sigma}
Let $\Sigma_n^*$ denote a set of representatives, modulo the action of
$\GL_N(\ensuremath{\mathcal{O}})$, of $n$-dimensional cells of $X^*_N$ that meet $X_N$.
Let $\Sigma^* = \cup_n \Sigma_n^*$.

The well-rounded retract $W$ is a proper, contractible
$\GL_N(\ensuremath{\mathcal{O}})$-complex.  Modulo $\GL_N(\ensuremath{\mathcal{O}})$, the cells in $W$ are in
bijection with cells in $\Sigma^*$, and the isomorphism classes of the
stabilizers are preserved under this bijection.  To see this, let
$\sigma(M)$ be a cell in $W$.  Then the forms $B \in \sigma(M)$ have
$M(B) = M$ and $m(B) = 1$.  Under the identification in
Theorem~\ref{theorem:duality}, $\sigma(M)$ corresponds to the cone
$C_\sigma'$ with spanning vectors $\set{q(v) \mid v \in M}$.  Let
$\sigma '$ be the corresponding cell.  There is a cell $\sigma \in
\Sigma^*$ that is $\Gamma$-equivalent to $\sigma'$.  If $\gamma \in
\GL_N(\ensuremath{\mathcal{O}})$, then $\gamma \cdot \sigma(M) =
\sigma(\set{(\gamma^*)^{-1}v \mid v \in M}) $.  Thus if $\gamma \in
\Stab(\sigma(M))$, then $(\gamma^*)^{-1} \in \Stab(C_\sigma')$.  It
is therefore clear that $\Stab(C_\sigma') = \Stab(\sigma') \simeq
\Stab(\sigma)$.

The map $g \mapsto (g^*)^{-1}$ is an isomorphism of groups, so the
duality between $W$ mod $\Gamma$ and $\Sigma^*$ in
Theorem~\ref{theorem:duality} allows us to work with $\Sigma^*$ and its
stabilizer subgroups as if we were working with a proper, contractible
$\GL_N(\ensuremath{\mathcal{O}})$-complex.  This fact will be used later to compute the
mass formula (\S\ref{sec:mass}).
\end{remark}

\section{The cohomology and homology}\label{sec:cohomologyandhomology}
The material in this section follows \cite[\S3]{PerfFormModGrp},
\cite[\S2]{soule-3torsion}, and \cite{agm4}, all of which rely on
\cite{brown}.  Recall that $\Gamma = \GL_N(\ensuremath{\mathcal{O}})$.  In this section, we
introduce a complex $\Vor_{N,D} = (V_*(\Gamma), d_*)$ of
$\ensuremath{\mathbb{Z}}[\Gamma]$-modules whose homology is isomorphic to the group
cohomology of $\Gamma$ modulo small primes.  More precisely, for any
positive integer $n$ let $\ensuremath{\mathcal{S}}_{n}$ be the Serre class of finite
abelian groups with orders only divisible by primes less than or equal
to $n$ \cite{serre}.  Then the main result of this section (Theorem
\ref{theorem:groupcoh}) is that, if $n=n (N,D)$ is larger than all the
primes dividing the orders of finite subgroups of $\Gamma$, then
modulo $\ensuremath{\mathcal{S}}_{n}$ the homology of $\Vor_{N,D}$ is isomorphic to the
group cohomology of $\Gamma$.

As above let $\Sigma_n^* = \Sigma^*_n(\Gamma)$ denote a finite set of
representatives, modulo the action of $\Gamma$, of $n$-dimensional
cells of $X^*_N$ which meet $X_N$.  A cell $\sigma$ is called
orientable if every element in $\Stab(\sigma)$ preserves the
orientation of $\sigma$.  By $\Sigma_n = \Sigma_n(\Gamma)$ we denote
the set of orientable cells in $\Sigma_n^*(\Gamma)$ Let $\Sigma^* =
\cup_n \Sigma^*_n$ and let $\Sigma = \cup_n \Sigma_n$.

\subsection{Steinberg homology}
The arithmetic group $\Gamma$---like
any arithmetic group---is a virtual duality group \cite[Theorem
11.4.4]{borel-serre}.  This means that there is a $\ensuremath{\mathbb{Z}}
[\Gamma]$-module $I$ that plays the role of an ``orientation module''
for an analogue of Poincar\'e duality: for all coefficient modules $M$
there is an isomorphism between the cohomology of $\Gamma '$ with
coefficients in $M$ and the homology of $\Gamma '$ with coefficients
in $M\otimes I$, for any torsion-free finite index subgroup $\Gamma
'\subset \Gamma $.  For more details see \cite[\S11]{borel-serre}.
The module $I$ is called the \emph{dualizing module} for $\Gamma$.
Borel and Serre prove that $I$ is isomorphic to
$H^{\nu} (\Gamma , \ensuremath{\mathbb{Z}} [\Gamma])$, where $\nu = \vcd \Gamma$ is the
\emph{virtual cohomological dimension} of $\Gamma$
(cf.~\S\ref{sec:mass}); note that for $\Gamma = \GL_{N}(\ensuremath{\mathcal{O}})$, we have
$\nu = N^{2}-N$.

One can show that $I$ is isomorphic to $\St_{\Gamma}\otimes \Omega $,
where $\St_{\Gamma}$ is the \emph{Steinberg module} (the top
non-vanishing reduced homology group of the Tits building attached to
$\GL_{N}/F$), and $\Omega$ is the \emph{orientation module} (which
records whether or not an element changes the orientation of the
symmetric space).  Since the action of $\GL_N(\ensuremath{\mathcal{O}})$ preserves the
orientation of $\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$, (cf.~Remark~\ref{rem:orientable}), the
module $\Omega$ is trivial, and thus $H^{\nu} (\Gamma , \ensuremath{\mathbb{Z}} [\Gamma])
= \St_\Gamma$

For any $\Gamma$-module $M$, the \emph{Steinberg homology}
\cite[p.~279]{brown}, denoted $H^{\St}_*(\Gamma, M)$, is defined as
$H^{\St}_*(\Gamma, M) = H_*(\Gamma, \St_\Gamma \otimes M)$.  Let
$\hat{H}^*$ denote the \emph{Farrell cohomology} $\hat{H}^*(\Gamma,M)$
of $\Gamma$ with coefficients in $M$ (cf.~\cite[X.3]{brown}).  There
is a long exact sequence
\begin{equation}\label{eq:long-exact}
\cdots \to H^{\St}_{\nu - i} \to H^i \to \hat{H}^i \to  H^{\St}_{\nu
  - (i + 1)} \to H^{i + 1}\to \hat{H}^{i + 1}\to \cdots
\end{equation}
that implies that we can understand the group cohomology by
understanding the Steinberg homology and the Farrell cohomology.  In
particular, if the Farrell cohomology vanishes, (which is the case
when the torsion primes in $\Gamma$ are invertible in
$M$, see \cite[IX.9~et~seq.]{brown}), then the
Steinberg homology is exactly the group cohomology $H^{\St}_{\nu -
k}(\Gamma, M) \simeq H^k(\Gamma, M)$.

We now specialize to the case $M = \ensuremath{\mathbb{Z}}$ with trivial $\Gamma$-action,
and in the remainder of this section omit the coefficients from
homology and cohomology groups.  Note that the well-rounded retract
$W$ is an $(N^2 - N)$-dimensional proper and contractible
$\Gamma$-complex.

\begin{proposition}\label{prop:steinberg-group}
Let $b$ be an upper bound on the torsion primes for $\Gamma =
\GL_N(\ensuremath{\mathcal{O}})$.  Then modulo the Serre class $\ensuremath{\mathcal{S}}_b$, we
have
\[
H^{\St}_{\nu - k}(\Gamma) \simeq H^k(\Gamma).
\]
\end{proposition}

\subsection{The Voronoi complex}\label{ss:vor}
Let $V_n(\Gamma)$ denote the free abelian group generated by
$\Sigma_n(\Gamma)$.  Let $d_n \colon V_n(\Gamma) \to V_{n -
  1}(\Gamma)$ be the map defined in \cite[\S3.1]{PerfFormModGrp}, and
  denote the complex $(V_{*} (\Gamma), d_{*})$ by $\Vor_{N,D}$.

Let $\partial X^*_N$ denote the cells in $X^*_N$ that do not meet $X$.  Then
$\partial X^*_N$ is a $\Gamma$-invariant subcomplex of $X^*_N$.
Let $H_*^{\Gamma }(X^*_N, \partial X_N^*)$ denote the relative
equivariant homology of the
pair $(X^*_N, \partial X^*_N)$ with integral coefficients
(cf.~\cite[VII.7]{brown}).

\begin{proposition}\label{prop:Voronoi-equivariant}
  Let $b$ be an upper bound on the torsion primes for $\Gamma =
  \GL_N(\ensuremath{\mathcal{O}})$.  Modulo the Serre class $\ensuremath{\mathcal{S}}_b$,
\[
  H_n(\Vor_{N,D}) \simeq H^\Gamma_n(X_N^*, \partial X_N^*).
\]
\end{proposition}
\begin{proof}
This result follows from \cite[Proposition~2]{soule-3torsion}.  The
argument is explained in detail for $F = \ensuremath{\mathbb{Q}}$ in
\cite[\S3.2]{PerfFormModGrp} and can be extended to imaginary quadratic $F$.
For the convenience of the reader we recall the argument.

There is a spectral sequence $E^r_{pq}$
converging to the equivariant homology groups
\[
H^\Gamma_{p + q}(X_N^*,
\partial X_N^*)
\]
of the homology pair $(X^*_N, \partial X^*_N)$ such
that
\begin{equation}\label{eq:brown}
E^1_{pq} = \bigoplus_{\sigma \in \Sigma^*_p}H_q(\Stab(\sigma),
\ensuremath{\mathbb{Z}}_\sigma) \Rightarrow H^\Gamma_{p + q}(X_N^*,
\partial X_N^*),
\end{equation}
where $\ensuremath{\mathbb{Z}}_\sigma$ is the orientation module for $\sigma$.

When $\sigma$ is not orientable the homology
$H_0(\Stab(\sigma), \ensuremath{\mathbb{Z}}_\sigma)$
is killed by $2$.
Otherwise, $H_0(\Stab(\sigma), \ensuremath{\mathbb{Z}}_\sigma) \simeq \ensuremath{\mathbb{Z}}_\sigma$.
Therefore modulo $\ensuremath{\mathcal{S}}_2$ we have
\[E^1_{n,0} \simeq \bigoplus_{\sigma \in \Sigma_n} \ensuremath{\mathbb{Z}}_\sigma.\]
 Furthermore, when $q > 0$, $H_q(\Stab(\sigma), \ensuremath{\mathbb{Z}}_\sigma)$ lies in
the Serre class $\ensuremath{\mathcal{S}}_b$, where $b$ is the upper bound on the torsion
primes for $\Gamma$. In other words, modulo $\ensuremath{\mathcal{S}}_b$, the $E^1$ page
of the spectral sequence \eqref{eq:brown} is concentrated in the
bottom row, and there is an identification of the bottom row groups with the
groups $V_n(\Gamma)$ defined above.  Finally,
\cite[Proposition~2]{soule-3torsion} shows that $d^1_n = d_n$.
\end{proof}
\begin{remark}
If we do not work modulo $\ensuremath{\mathcal{S}}_b$, then there are other
entries in the spectral sequence to consider before we get
$H^\Gamma_n(X^*_{N}, \partial X^*_{N})$.  In particular, the small torsion
in $H_*(\Vor_{N,D})$ may not agree in general with the small torsion
in $H^\Gamma_{p + q}(X_N^*, \partial X_N^*)$.
\end{remark}

\subsection{Equivariant relative homology to Steinberg homology}
For $V \subset F^N$ a proper subspace, let $C(V)$ be the set of
matrices $A \in C_N^*$ such that the kernel of $A$ is $V
\otimes_{\ensuremath{\mathbb{Q}} } \ensuremath{\mathbb{R}}$, and let $X(V) = \pi(C(V))$.  The closure
$\overline{C(V)}$ of $C(V)$ in the usual topology induced from
$\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ consists of matrices whose kernel contains $V$.  Let
$\overline{X(V)} = \pi(\overline{C(V)})$.

The following Lemma shows that $\overline{X(V)}$ is contractible.
Soul\'e proves this for $F = \ensuremath{\mathbb{Q}}$ in \cite[Lemma~1]{soule-3torsion}
and \cite[Lemma~2]{soule-addendum}, and the same proof with simple
modifications applies to our setting:

\begin{lemma}\label{lem:contractible}
  For any proper subspace $V \subset F^N$, the CW-complex
  $\overline{X(V)}$ is contractible.
\end{lemma}
\begin{proof}
Let $A$ be a perfect Hermitian form over $F$.  Then $\sigma(A) \cap
\overline{X(V)}$ is the perfect cone of the intersection $M(A) \cap
V^\perp$, where $V^\perp$ is the orthogonal complement to $V$ in
$F^N$.  It follows that $\overline{X(V)}$ is a sub-CW-complex of
$X^*_N$.

Now we show that $\overline{X(V)}$ is contractible.  Since
$\overline{X(V)}$ has the same homotopy type as $X (V)$, it is enough to prove that
$X (V)$ is contractible.  For this it suffices to prove that the
CW-topology on ${X(V)}$ coincides with its usual topology, since
$X(V)$ in this topology is clearly contractible (because it is convex).
To prove this, we argue that the covering of $X (V)$ by the closed
sets of the form $\sigma (A)\cap X (V)$, where $A$ is perfect, is
locally finite.

Given any positive definite Hermitian form $A$, let $A^{\perp}$ be its
restriction to the real space $V_{\ensuremath{\mathbb{R}}}^{\perp} = V^{\perp}\otimes_{\ensuremath{\mathbb{Q}} }
\ensuremath{\mathbb{R}}$.  Then $X (V)$ is isomorphic to the symmetric space for the group
$\Aut (V_{\ensuremath{\mathbb{R}}}^{\perp})$ via the map $A\mapsto A^{\perp}$.  If $\gamma \in
\Gamma$ satisfies $\gamma \cdot X (V)\cap X (V)\not = \emptyset$, then
in fact $\gamma $ stabilizes $V$.  Let $P\subset \Gamma$ be the
stabilizer of $V$ and let $\alpha \colon P\rightarrow \Gamma ' = \Aut
(V^{\perp}_{\ensuremath{\mathbb{R}}}\cap L)$ be the projection map.  Then the
set of cells $\sigma (A)\cap X (V)$, as $A$ ranges over the perfect
Hermitian forms, is finite modulo $\Gamma '$: if $[\sigma (A)\cap X
(V)] \cap [(\gamma \cdot \sigma (A)) \cap X (V)] \not =\emptyset$, then
$\gamma \in P$, and thus $(\gamma \cdot \sigma (A) \cap X (V)) =
\alpha (\gamma) (\sigma (A)\cap X (V))$.

To conclude the argument one uses Siegel sets; we refer to \cite{borel}
for their definition and to \cite[Ch.~II]{AMRT_SmoothCompactification}
for their properties in our setting.  Given any point $x\in X (V)$,
one can find an open set $U\ni x$ and a Siegel set $\mathfrak{S} $ such that
$\mathfrak{S}\supset U$.  Furthermore, any cell of the form $\sigma (A)\cap X
(V)$ is itself contained in another Siegel set $\mathfrak{S}'$.  Thus if $U$
meets $\gamma (\sigma (A)\cap X (V))$ for $\gamma \in \Gamma '$, we
must have that $\gamma \cdot \mathfrak{S} '$ meets $\mathfrak{S} $.  But this is only
possible for finitely many $\gamma$ by the standard properties of Siegel
sets.  Thus the covering of $X (V)$ by the closed sets of the form
$\sigma (A)\cap X (V)$ ($A$ perfect) is locally finite, which
completes the proof.
\end{proof}

\begin{proposition}\label{prop:relative-steinberg}
  For every $n \geq 0$, there are canonical isomorphisms of
  $\Gamma$-modules
\[H_n(X^*_N, \partial X^*_N) \simeq
\begin{cases}
  \St_\Gamma & \text{if $n = N - 1$,}\\
0 & \text{otherwise.}
\end{cases}
\]
\end{proposition}
\begin{proof}
Since $X^*_N$ is contractible, the long exact sequence for the pair
$(X^*_N, \partial X^*_N)$
\[\cdots \to H_n(\partial X^*_N) \to H_n(X^*_N) \to H_n(X^*_N, \partial
X^*_N) \to H_{n - 1}(\partial X^*_N) \to H_{n-1}(X^*_N) \to \cdots\]
implies $H_n(X^*_N, \partial X^*_N) \simeq \tilde{H}_{n-1}(\partial
X^*_N)$ for $n \geq 1$, where $\tilde{H}$
denotes reduced homology.

From properties of Hermitian forms, it is clear that $\overline{X(V)}
\cap \overline{X(W)} = \overline{X(V \cap W)}$
for each pair of proper, non-zero subspaces $V, W  \subset F^N$.  Thus the nerve
of the covering $\mathcal{U}$ of $\partial X^*_N$ by $\overline{X(V)}$
as $V$ ranges
over non-zero, proper subspaces of $F^N$ is the spherical Tits
building $T_{F,N}$.  Since each $\overline{X(V)}$ in $\mathcal{U}$ is
contractible by Lemma~\ref{lem:contractible}, the cover $\mathcal{U}$
is a good cover (i.e.,~non-empty finite intersections are diffeomorphic
to $\ensuremath{\mathbb{R}}^{d}$ for some $d$).  It follows that the relative homology groups are
isomorphic, i.e. $\tilde{H}_n(\partial X^*_N) \simeq \tilde{H}_n(T_{F,N})$.
By the Solomon-Tits theorem \cite{solomon}, the latter is isomorphic to $\St_\Gamma$
if $n = N - 2$ and is trivial otherwise.
\end{proof}

\begin{proposition}\label{prop:equivariant-steinberg} For all $n$, we
have
\[  H^\Gamma_{n}(X^*_N, \partial X^*_N) = H^{\St}_{n - (N - 1)}(\Gamma).\]
\end{proposition}
\begin{proof}
There is a spectral sequence \cite[equation (2)]{soule-3torsion}
computing the relative equivariant homology
\[E^2_{pq} = H_p(\Gamma, H_q(X_N^*, \partial X_N^*)) \Rightarrow
H^\Gamma_{p + q}(X^*_N, \partial X^*_N).\]
Proposition~\ref{prop:relative-steinberg} implies the $E^2$ page of
the spectral sequence is concentrated in the $q = N - 1$ column.  Then
\[H^\Gamma_{p + (N - 1)}(X^*_N, \partial X^*_N) \simeq
 H_p(\Gamma, H_q(X^*_N, \partial X^*_n)) = H_p(\Gamma, \St_\Gamma) =
 H^{\St}_p(\Gamma),\]
and the result follows.
\end{proof}

\begin{theorem} \label{theorem:groupcoh}
Let $b$ be an upper bound on the torsion primes for $\GL_N(\ensuremath{\mathcal{O}})$.
Modulo the Serre class $\ensuremath{\mathcal{S}}_b$,
\[H_n(\Vor_{N,D}) \simeq H^{N^2 - 1 - n}(\GL_N(\ensuremath{\mathcal{O}})).\]
\end{theorem}
\begin{proof}
Let $\Gamma = \GL_N(\ensuremath{\mathcal{O}})$.
Modulo $\ensuremath{\mathcal{S}}_b$, Propositions~\ref{prop:Voronoi-equivariant},
    \ref{prop:equivariant-steinberg}, and \ref{prop:steinberg-group} imply
\[H_n(\Vor_{N,D}) \simeq H_n^\Gamma(X^*_N, \partial X^*_N)
\simeq H^{\St}_{n - (N -  1)}(\Gamma) \simeq H^{N^2 - 1 - n}(\Gamma).\]
\end{proof}

\subsection{Torsion elements in \texorpdfstring{$\Gamma$}{Gamma}}\label{subsec:torsion}
To finish this section we discuss the possible torsion that can arise
in the stabilizer subgroups of cells in our complexes.  This allows us
to make the bound in Theorem~\ref{theorem:groupcoh} effective.

\begin{lemma}\label{lem:cyclotomic}
  Let $p$ be an odd prime, and let $F/\ensuremath{\mathbb{Q}}$ be a quadratic field.  Let
  $\Phi_p = x^{p-1} + x^{p-2} + \dots + x + 1$ be the $p^\text{th}$
  cyclotomic polynomial.  Then $\Phi_p$ factors over $F$ as a product
  of irreducible polynomials of degree $(p-1)/2$ if $F =
  \ensuremath{\mathbb{Q}}(\sqrt{p^*})$, where $p^* = (-1)^{(p-1)/2}p$, and is irreducible
  otherwise.
\end{lemma}
\begin{proof}
Let $\zeta_p$ denote a primitive $p^\text{th}$ root of unity.
Consider the diagram of Galois extensions below:
\[
\begin{tikzpicture}
  \node (Fz) {$F(\zeta_p)$};
  \node (K) [below of = Fz, right of = Fz]{$\ensuremath{\mathbb{Q}}(\zeta_p)$};
  \node (F) [below of = Fz, left of = Fz]{$F$};
  \node (KnF)[below of=K, left of=K] {$F \cap \ensuremath{\mathbb{Q}}(\zeta_p)$};
  \node (Q) [below of=KnF]{$\ensuremath{\mathbb{Q}}$};
  \draw[-] (Fz) to (K);
  \draw[-] (Fz) to (F);
  \draw[-] (KnF) to (Q);
  \draw[-] (F) to (KnF);
  \draw[-] (K) to (KnF);
\end{tikzpicture}
\]

If $F = \ensuremath{\mathbb{Q}}(\sqrt{p^*})$, then $F$ is the unique quadratic subfield of
$\ensuremath{\mathbb{Q}}(\zeta_p)$.  It follows that $\Phi_p$ factors over $F$ as a
product of irreducible polynomials of degree $(p-1)/2$.  If $F \neq
\ensuremath{\mathbb{Q}}(\sqrt{p^*})$, then $F \cap \ensuremath{\mathbb{Q}}(\zeta_p) = \ensuremath{\mathbb{Q}}$.  It follows that
$[F(\zeta_p):\ensuremath{\mathbb{Q}}(\zeta_p)] = [F: \ensuremath{\mathbb{Q}}] =  2$ and so $[F(\zeta_p): F] =
p-1$.  Thus $\Phi_p$ is irreducible over $F$ if $F \neq
\ensuremath{\mathbb{Q}}(\sqrt{p^*})$.
\end{proof}

\begin{lemma}\label{lem:orderp}
  Let $p$ be an odd prime, and let $F/\ensuremath{\mathbb{Q}}$ be a imaginary quadratic
  field.  If $g \in \GL_N(F)$ has order $p$, then
\[p \leq
\begin{cases}
  N + 1 & \text{if $p \equiv 1 \bmod{4}$,}\\
  2N + 1 & \text{otherwise.}
\end{cases}
\]
\end{lemma}
\begin{proof}
If $p \equiv 1 \bmod{4}$, then $p^* > 0$.  In particular, $F \neq
\ensuremath{\mathbb{Q}}(\sqrt{p^*})$ and so by Lemma~\ref{lem:cyclotomic}, $\Phi_p$ is
irreducible over $F$.  Then the minimal polynomial of $g$ is
$\Phi_p$.  By the Cayley-Hamilton Theorem, $\Phi_p$ divides the
characteristic polynomial of $g$.  Therefore $p-1 \leq N$.
Similarly, if $p \not \equiv 1 \bmod{4}$, then  $(p-1)/2 \leq N$.
\end{proof}

Lemmas~\ref{lem:cyclotomic} and \ref{lem:orderp} immediately imply the following:

\begin{proposition}\label{prop:stabprimes}
If $g \in \GL_3(\ensuremath{\mathcal{O}}_{F})$ has prime order $q$, then $q \in
\set{2,3,7}$ for $F = \ensuremath{\mathbb{Q}}(\sqrt{-7})$ and $q\in \set{2,3}$ otherwise.
If $g \in \GL_4(\ensuremath{\mathcal{O}}_{F})$ has prime order $q$, then $q \in \set{2,3,5,7}$ for
$F = \ensuremath{\mathbb{Q}}(\sqrt{-7})$ and $q\in \set{2,3,5}$ otherwise.
\end{proposition}

In Tables~\ref{tab:diff3-3}--\ref{tab:diff4-4} we give the
factorizations of the orders of the stabilizers of the cells in
$\Sigma^{*}$.

\section{A mass formula for the Voronoi complex}\label{sec:mass} The
computation of the cell complex is a relatively difficult task and
there are many ways
that small mistakes in the computation could cause the final answer to be incorrect.
Hence it is very important to have checks that allow us to give strong
evidence for the correctness of our computations.  One is that the
complexes we construct actually are chain complexes, namely that their
differentials square to zero.  Another is the \emph{mass formula},
stated in Theorem~\ref{theorem:mass}.  According to this formula, the
alternating sum over the cells of $\Sigma^*$ of the inverse orders of
the stabilizer subgroups must vanish.  A good reference for this
section is \cite[Ch.~IX, \S\S 6--7]{brown}, and we follow it closely.
The main theorem underlying this computation is due to Harder
\cite{harder}.

\subsection{Euler characteristics}\label{subsec:euler}
We begin by recalling some definitions.  The \emph{cohomological
dimension} $\cd\Gamma$ of a group $\Gamma$ is the largest $n\in \ensuremath{\mathbb{Z}}
\cup \{\infty \}$ such that there exists a $\ensuremath{\mathbb{Z}} \Gamma$-module $M$
with $H^{n} (\Gamma ; M)\not = 0$.  The \emph{virtual cohomological
dimension} $\vcd \Gamma$ of $\Gamma$ is defined to be the
cohomological dimension of any torsion-free finite index subgroup of
$\Gamma$ (one can show that this is well-defined).  We recall that
$\Gamma$ is said to be of \emph{finite homological type} if (i) $\vcd
\Gamma <\infty$ and (ii) for every $\ensuremath{\mathbb{Z}} \Gamma$-module $M$ that is
finitely generated as an abelian group, the homology group $H_{i}
(\Gamma ; M)$ is finitely generated for all $i$.

\begin{definition}
Let $\Gamma$ be a torsion-free group of finite homological type, and
let $H_{*} (\Gamma)$ denote the homology of $\Gamma$ with (trivial)
$\ensuremath{\mathbb{Z}}$-coefficients.  The \emph{Euler characteristic} $\chi(\Gamma)$ is
\[\chi(\Gamma) = \sum_{i}(-1)^i\rank_\ensuremath{\mathbb{Z}}(H_i(\Gamma)).\]
\end{definition}

\begin{proposition} \label{prop:euler-torsion}
\rftheoremcite{Theorem 6.3}{brown} If $\Gamma$ is torsion-free and of finite homological type and
$\Gamma' \subseteq \Gamma$ is a subgroup of finite index, then
\[\chi(\Gamma') = [\Gamma \colon \Gamma'] \cdot \chi(\Gamma).\]
\end{proposition}

One can use Proposition~\ref{prop:euler-torsion} to extend the notion of Euler
characteristic to groups with torsion.  Namely, if $\Gamma$ is an
arbitrary group of finite homological type with a torsion-free subgroup
$\Gamma'$ of finite index, one sets
\begin{equation}\label{eq:tordef}
\chi(\Gamma) = \frac{\chi(\Gamma')}{[\Gamma \colon \Gamma']}.
\end{equation}
By Proposition \ref{prop:euler-torsion}, this is independent of the
choice of $\Gamma '$.

\begin{proposition}\label{prop:euler-finite}
If $\Gamma$ is a finite group, then $\chi(\Gamma) = 1/|\Gamma|$.
\end{proposition}
\begin{proof}
  Take $\Gamma'$ to be the trivial subgroup in \eqref{eq:tordef}.
\end{proof}

\begin{theorem} \label{theorem:euler-characteristic}
\theoremcite{harder}  Let $F$ be a number field with ring of integers $\ensuremath{\mathcal{O}}$, and let
$\zeta_{F} (s)$ be the Dedekind zeta function of $F$.  Then
\[\chi(\SL_N(\ensuremath{\mathcal{O}})) = \prod_{k = 2}^N \zeta_F(1 - k).\]
\end{theorem}
In particular, since
for imaginary quadratic $F$ the Dedekind zeta function
$\zeta_{F} (m)$
vanishes for
all $m \in \mathbb{Z} < 0$,
we have that $\chi(\SL_N(\ensuremath{\mathcal{O}})) = 0$
when $N \geq 2$.

\begin{corollary}\label{cor:euler-characteristic-gl}
Let $F$ be an imaginary quadratic field with ring of integers $\ensuremath{\mathcal{O}}$.
Then $\chi(\GL_N(\ensuremath{\mathcal{O}})) = 0$ when $N \geq 2$.
\end{corollary}
\begin{proof}
This follows immediately from Theorem~\ref{theorem:euler-characteristic}
and the definition of the Euler characteristic since $\SL_N(\ensuremath{\mathcal{O}})$ is of
finite index in $ \GL_N(\ensuremath{\mathcal{O}})$.
\end{proof}

Now we turn to a different concept, the \emph{equivariant Euler
characteristic} $\chi_{\Gamma} (X)$ of $\Gamma$. Here $X$ is any cell
complex with
$\Gamma$ action such that (i) $X$ has finitely many cells mod
$\Gamma$, and (ii) for each $\sigma\in X$, the stabilizer subgroup
$\Stab_{\Gamma } (\sigma)$ is finite.  One defines
\[
\chi_{\Gamma} (X) = \sum_{\sigma \in S} (-1)^{\Dim \sigma}
\chi(\Stab_{\Gamma }(\sigma)),
\]
where $S$ is a set of representatives of cells of $X$ mod $\Gamma$.

\subsection{The mass formula}\label{subsec:mass}
The well-rounded retract $W$ defined in \S\ref{subsec:cell-complexes}
is a proper, contractible $\Gamma$-complex and so its equivariant
Euler characteristic is defined.  We compute its equivariant Euler
characteristic, phrased in terms of cells in $\Sigma^*$ using
Remark~\ref{rem:W-Sigma}, to get a mass formula.

\begin{theorem}[Mass Formula] \label{theorem:mass}
We have
\[\sum_{\sigma \in \Sigma^*_k} (-1)^k\frac{1}{|\Stab_{\Gamma }(\sigma)|} = 0.\]
\end{theorem}
\begin{proof}
Let $\Gamma = \GL_{N} (\ensuremath{\mathcal{O}})$.  By Proposition \ref{prop:euler-finite}
we have $\chi (\Stab_{\Gamma } (\sigma)) = 1/|\Stab_\Gamma (\sigma)|$.
Thus the result follows from Corollary
\ref{cor:euler-characteristic-gl} if we can show $\chi_{\Gamma} (W) =
\chi (\Gamma)$, where $W$ is the $\Gamma$-complex formed from the
well-rounded retract, using Theorem~\ref{theorem:duality} and
Remark~\ref{rem:W-Sigma} to identify cells of $W$ mod $\Gamma$ with
cells in $\Sigma^*$.  But according to \cite[Proposition 7.3
(e')]{brown}, this equality is true if $\Gamma$ has a torsion-free
subgroup of finite index, which is a standard fact for $\Gamma$.
\end{proof}

Using the stabilizer information in
Tables~\ref{tab:diff3-3}--\ref{tab:diff4-4} one can easily verify that
$\chi(\GL_n(\ensuremath{\mathcal{O}})) = 0$ for each of our examples.  For instance, if we
add together the terms $1/|\Stab_{\Gamma }(\sigma)|$ for cells $\sigma$ of the
same dimension to a single term for $\GL_4(\ensuremath{\mathcal{O}}_{-4})$, we find
(ordering the terms by increasing dimension in $\Sigma^{*}$)
\begin{multline*}
-\frac{11}{3072} +\frac{127}{960} -\frac{4187}{2304} +\frac{28375}{2304} -\frac{868465}{18432} +\frac{126127}{1152} -\frac{81945}{512}\\ +\frac{340955}{2304} -\frac{48655}{576} +\frac{16075}{576} -\frac{21337}{4608} +\frac{101}{384} -\frac{17}{92160}  = 0.
\end{multline*}
The other groups have been checked similarly.

\section{Explicit homology classes}\label{sec:explicit}

By Theorem~\ref{theorem:groupcoh} we have $H_{N^2 - 1}(\Vor_{N,D} \otimes
\ensuremath{\mathbb{Q}}) \simeq H^0(\GL_N(\ensuremath{\mathcal{O}}), \ensuremath{\mathbb{Q}})$, which in turn is isomorphic to
$\ensuremath{\mathbb{Q}}$.  This suggests that there should be a canonical generator for
this homology group, a fact already explored in \cite[Section
5]{PerfFormModGrp}.  An obvious choice is the analogue of the chain
presented there, namely
\[
\xi := \xi_{N,D} := \sum_\sigma \frac{1}{|\Stab(\sigma)|} [\sigma],
\]
where $\sigma$ runs through the cells in $\Sigma_{N^2 -
1}(\GL_N(\ensuremath{\mathcal{O}}))$.  In this section we verify that this is true for our
examples.  We should point out that all the cells in $\Sigma_{N^2 -
1}(\GL_N(\ensuremath{\mathcal{O}}))$ are orientable. The reason is that the group
$\GL_n(\ensuremath{\mathbb{C}})$ is connected and so the determinant of its action on
$\ensuremath{\mathcal{H}}^N(\ensuremath{\mathbb{C}})$ is positive.  Since the stabilizers are included in
$\GL_N(\ensuremath{\mathcal{O}})$ and the faces in $\Sigma_{N^2 - 1}(\GL_N(\ensuremath{\mathcal{O}}))$ are
full-dimensional the orientation has to be preserved.

\begin{theorem}\label{theorem:explicit}
When $N = 3$ and $D\geq -24$ or $N = 4$ and $D=-3,-4$, the chain $\xi$
is a cycle and thus generates $H_{N^2 - 1}(\Vor_{N,D} \otimes \ensuremath{\mathbb{Q}})$.
\end{theorem}
\begin{proof}

The proof is an explicit computation with differential matrix $A$
representing the map $V_{N^{2}-1} (\Gamma)\rightarrow V_{N^{2}-2}
(\Gamma)$ (cf.~\S \ref{ss:vor}).  Note that the signs of the entries
of $A$ depend on a choice of orientation for each of the cells in
$\Sigma_{N^2 - 1}$ and $\Sigma_{N^2 - 2}$.  In each non-zero row of
$A$, there are exactly two non-zero entries.  Each non-zero entry
$A_{i,j}$ has absolute value $|\Stab(\sigma_j)|/|\Stab(\tau_i)|$,
where $\sigma_j \in \Sigma_{N^2 - 1}(\GL_N(\ensuremath{\mathcal{O}}))$ and $\tau_i \in
\Sigma_{N^2 - 2}(\GL_N(\ensuremath{\mathcal{O}}))$.  One then checks that there is a choice
of orientations such that the non-zero entries in a given row have
opposite signs.
\end{proof}

For example, consider $N = 4$ and $D = -4$.  The differential matrix is
\[d_{15} = \begin{bmatrix}
0 & 0\\ 1920 & -256
\end{bmatrix},\] with kernel generated by $(2, 15) = 92160
({1}/{46080}, {1}/{6144})$.  The orders of the two stabilizer groups
for the cells in $\Sigma_{15}(\GL_N(\ensuremath{\mathcal{O}}_{-4}))$ are $46080$ and
$6144$, respectively, and thus $\xi_{4,-4}$ is a cycle.

\begin{remark}
It seems likely that Theorem~\ref{theorem:explicit} holds for all $\GL_{N}
(\ensuremath{\mathcal{O}}_{D})$, although we do not presently have a proof.
\end{remark}

\section{Results and tables}\label{sec:tables}

We conclude by presenting the results of our computations.  First we
summarize the most interesting numerical results from the tables,
namely those involving the homology computations for $\GL_*(\ensuremath{\mathcal{O}}_D)$.

\begin{theorem}
Let $0>D\geq -24$ be a fundamental discriminant.
\begin{enumerate}
\item The Voronoi homology $H_n(\Vor_{3,D})$ can have non-trivial rank
only for $n=2,4,5$ and $8$.
\item In each case for $n=8$ we have $H_8(\Vor_{3,D})\cong \ensuremath{\mathbb{Z}}$.
\item In homology degree $2$, $4$ and $5$, we have the following table
  of ranks (empty slots being zero)
\[
\begin{array}{|c|rrrrrrrrrr|}
\hline
D \phantom{\Big|} & -3 &-4 &-7 &-8 &-11 &-15 &-19 &-20 &-23 &-24  \\
\hline
\#\{\text{perf.~forms}\} \phantom{\Big|}  &2&1&2&2&12&90&157&212&870&596  \\
\hline
\rank\ H_2\phantom{\big|} & &&&&&&&1&4&1\\
\rank\ H_4\phantom{\big|} &1 &1&1&1&1&3&2&4&5&5\\
\rank\ H_5\phantom{\big|} & 1&1&2&2&2&5&3&6&1&7\\
\hline
\end{array}
  \]

\medskip
\item The $p$-Sylow subgroups for the primes occurring in (the orders of) the respective
  homology groups are given in the following table, where we use
  the shorthand $n( G)$ to denote that the $p$-Sylow subgroup in
  $H_n(\Vor_{3,D})$ is equal to $G$.
 \def \mathbb Z {\mathbb Z}
\[
\begin{array}{|c|cccccccc|}
\hline
D \phantom{\Big|} & -3&-4&-7&-8&-11&-15&-19&-20\\
\hline
\Syl_2&  & 5(\mathbb Z_2)& & 4(\mathbb Z_2), 5(\mathbb Z_2)& 5(\mathbb Z_4) & 4(\mathbb Z_4) & 5(\mathbb Z_2\times \mathbb Z_8)&3(\mathbb Z_2),4(\mathbb Z_2^4), 7(\mathbb Z_2^2\times \mathbb Z_4)  \\
\Syl_3& 7(\mathbb Z_9) &&7(\mathbb Z_3) && 7(\mathbb Z_3^2)&7(\mathbb Z_3^5)&7(\mathbb Z_3^4)&5(\mathbb Z_3), 7(\mathbb Z_3^5)\\
\Syl_5&  && && &7(\mathbb Z_5)&&\\
\Syl_7&  && 5(\mathbb Z_7)&& &&& 7(\mathbb Z_7)\\
\hline
\end{array}
\]

\[
\begin{array}{|c|cc|}
\hline
D \phantom{\Big|} &-23&-24 \\
\hline
\Syl_2
& 3(\mathbb Z_2),4(\mathbb Z_2^{13}),5(\mathbb Z_2^5),7(\mathbb Z_2^3\times\mathbb Z_4^2)& 3(\mathbb Z_2^7),4(\mathbb Z_2^5),5(\mathbb Z_2^8\times \mathbb Z_4^2),6(\mathbb Z_2),7(\mathbb Z_2^4\times\mathbb Z_4^2) \\
\Syl_3
& 4(\mathbb Z_3^3),7(\mathbb Z_3^7)& 7(\mathbb Z_3^3) \\
\hline
\end{array}
\]

\medskip
\medskip

\begin{comment}

\item The $p$-torsion ranks for $p=2$ and $3$ for the respective
  homology groups are comprised in the following table, where we use
  the shorthand $n( r)$ if the $p$-rank in
  $H_n(\Vor_{3,D})$ equals $r$.
\[
\begin{array}{|c|cccccc|}
\hline
D \phantom{\Big|} & -3&-4&-7&-8&-11&-15\\
\hline
2\text{-}\rank >0&  &5(1)& &4(1),5(1)&5(1)&4(4),7(3)\\
3\text{-}\rank >0& 7(1) &&7(1) && 7(2)&7(5)\\
\hline
\end{array}
\]

\[
\begin{array}{|c|cccc|}
\hline
D \phantom{\Big|} &-19&-20&-23&-24 \\
\hline
2\text{-}\rank >0&5(2),7(3)&3(1),4(4), 7(3) & 3(1),4(13),5(5),7(5)& 3(7),4(5),5(10),6(1),7(6) \\
3\text{-}\rank >0&7(4)&5(1), 7(5)& 3(3),7(7)& 7(5) \\
\hline
\end{array}
\]
\item Finally, the $p^k$-subgroups for $p^{k}>3$ and $k>0$ are as follows.
\smallskip
\begin{itemize}
\item  $D=-3$, $H_7 > \ensuremath{\mathbb{Z}}/9\ensuremath{\mathbb{Z}} $.
\item  $D=-7$, $H_5 > \ensuremath{\mathbb{Z}}/7\ensuremath{\mathbb{Z}}$.
\item  $D=-8$, $H_5 > \ensuremath{\mathbb{Z}}/4\ensuremath{\mathbb{Z}}$.
\item  $D=-11$, $H_5 > \ensuremath{\mathbb{Z}}/4\ensuremath{\mathbb{Z}}$.
\item  $D=-15$, $H_7 > \ensuremath{\mathbb{Z}}/5\ensuremath{\mathbb{Z}}$.
\item  $D=-19$, $H_5 > \ensuremath{\mathbb{Z}}/8\ensuremath{\mathbb{Z}}$.
\item  $D=-20$, $H_7 > \ensuremath{\mathbb{Z}}/7\ensuremath{\mathbb{Z}}$.
\item  $D=-23$, $H_7 > (\ensuremath{\mathbb{Z}}/4\ensuremath{\mathbb{Z}})^2$.
\item  $D=-24$, $H_5\,,H_7  > (\ensuremath{\mathbb{Z}}/4\ensuremath{\mathbb{Z}})^2$.
\end{itemize}
 \end{comment}
\end{enumerate}
\end{theorem}

\begin{theorem}
Let $V_{4}$ be either the Voronoi complex for $\GL_{4} (\ensuremath{\mathcal{O}}_{-3})$ or
$\GL_{4} (\ensuremath{\mathcal{O}}_{-4})$.  Then
 $H_n(V_{4})$ has non-trivial rank only for $n=6,\,7,\,9,\,10,\,12$, and $15$.
 More precisely,
 \begin{enumerate}
\item For $n=3k$ we have $\rank H_{3k}(V_{4})= \begin{cases}1 &\text{if }k>1,\\
0& \text{if }k=1.\end{cases}$
\item For $n=7$ we have $\rank H_{n}(V_{4})=2\,$.
\item For $n=10$ we have $\rank H_{n}(V_{4})=1\,$.
\end{enumerate}
Moreover, the only torsion primes appearing
are 2, 3 and  5.
\end{theorem}

Invoking Theorem~\ref{theorem:groupcoh},  we also get information on
the cohomology:

\begin{theorem}
Let $D = -3$ or $-4$.  Then
 $H^n(\GL_4(\ensuremath{\mathcal{O}}_{D}))$ has non-trivial rank only for
$n = 0,\,3,\,5,\,6,\,8,$ and  $9$.
 More precisely,
 \begin{enumerate}
\item For $n = 3k$ we have
$\rank H^n(\GL_4(\ensuremath{\mathcal{O}}_D))= 1$ if $0 \leq k \leq 3$.
\item For $n=5$ we have $\rank H^n(\GL_4(\ensuremath{\mathcal{O}}_D))=1$.
\item For $n=8$ we have $\rank H^n(\GL_4(\ensuremath{\mathcal{O}}_D))=2$.
\end{enumerate}
\end{theorem}

The remainder of the text is devoted to presenting details about the
Voronoi complexes and full information about their homology.  The
notation is as follows:

\begin{itemize}
\item The first three columns concern the cell decomposition of $X_{N}^{*}$ mod $\Gamma$:
\begin{itemize}
\item $n$ is dimension of the cells in the (partially) compactified
symmetric space $X_{N}^{*}$.
\item $|\Sigma^{*}_{n}|$ is the number of $\Gamma$-orbits in the
cells that meet $X_{N}\subset X_{N}^{*}$.
\item $|\Stab |$ gives the sizes of the stabilizer
subgroups in factored form.  The notation $A (k)$ means that, of the
$|\Sigma^{*}_{n}|$ cells of dimension $n$, $k$ of them have a
stabilizer subgroup of order $A$.
\end{itemize}
\item The next four columns concern the differentials $d_{n}$ of the
Voronoi complex $\Vor_{N,D}$:
\begin{itemize}
\item $|\Sigma_{n}|$ is the number of orientable $\Gamma$-orbits in $\Sigma^{*}_{n}$.
\item $\Omega$ is the number of nonzero entries in the differential $d_{n}\colon V_{n} (\Gamma )\rightarrow V_{n-1} (\Gamma)$.
\item $\rank$ is the rank of $d_{n}$.
\item elem.~div.~gives the elementary divisors of $d_{n}$.  As in the
stabilizer column, the notation $d (k)$ means that the elementary
divisor $d$ occurs with multiplicity $k$.  If the rank of $d_{n}$
is zero then this column is left empty.
\end{itemize}
\item Finally, the last column gives the homology of the Voronoi
complex.  One can easily check that $H_{n}\simeq \ensuremath{\mathbb{Z}}^{r}\oplus
\bigoplus (\ensuremath{\mathbb{Z}}/{d}\ensuremath{\mathbb{Z}} )^{k}$, where $r=|\Sigma_{n}|-\rank
(d_{n})-\rank (d_{n+1})$ and the sum is taken over the elementary
divisors $d (k)$ from row $n+1$.  To save space, we abbreviate $\ensuremath{\mathbb{Z}}
/d\ensuremath{\mathbb{Z}}$ by $\ensuremath{\mathbb{Z}}_{d}$.  By Theorem \ref{theorem:groupcoh} we have
$H_n(\Vor_{N,D}) \simeq H^{N^2 - 1 - n}(\GL_N(\ensuremath{\mathcal{O}}))$ modulo the
torsion primes in $\GL_N(\ensuremath{\mathcal{O}})$.  These primes are visible in the third
column of each table.
\end{itemize}

\newcolumntype{x}[1]{>{\RaggedRight}p{#1}}
\input{gl.table.3.-3}
\input{gl.table.3.-4}
\input{gl.table.3.-7}
\input{gl.table.3.-8}
\input{gl.table.3.-11}
\input{gl.table.3.-15}
\input{gl.table.3.-19}
\input{gl.table.3.-20}
\input{gl.table.3.-23}
\input{gl.table.3.-24}
\input{gl.table.4.-3}
\input{gl.table.4.-4}
\clearpage

\bibliographystyle{amsplain_initials_eprint}