This paper has been withdrawn by the authors for future revision.

The cosmic crystallography method of Lehoucq et al. [1] produces sharp peaks in the distribution of distances between the images of cosmic sources. But the method cannot be applied to universes with compact spatial sections of negative curvature. We apply to these a second order crystallographic effect, as evidenced by statistical parameters.

This article investigates the computation of the eigenmodes of the Laplacian operator in multi-connected three-dimensional spherical spaces. General mathematical results and analytical solutions for lens and prism spaces are presented. Three complementary numerical methods are developed and compared with our analytic results and previous investigations. The cosmological applications of these results are discussed, focusing on the cosmic microwave background (CMB) anisotropies. In particular, whereas in the Euclidean case too small universes are excluded by present CMB data, in the spherical case there will always exist candidate topologies even if the total energy density parameter of the universe is very close to unity.

This article gives the construction and complete classification of all three-dimensional spherical manifolds, and orders them by decreasing volume, in the context of multiconnected universe models with positive spatial curvature. It discusses which spherical topologies are likely to be detectable by crystallographic methods using three-dimensional catalogs of cosmic objects. The expected form of the pair separation histogram is predicted (including the location and height of the spikes) and is compared to computer simulations, showing that this method is stable with respect to observational uncertainties and is well suited for detecting spherical topologies.

We present a modified version of the cosmic crystallography method, especially useful for testing closed models of negative spatial curvature. The images of clusters of galaxies in simulated catalogs are ``pulled back'' to the fundamental domain before the set of distances is calculated.

We try to apply the cosmic crystallography method of Lehoucq, Lachieze-Rey, and Luminet to a universe model with closed spatial section of negative curvature. But the sharp peaks predicted for Einstein-de Sitter closed models do not appear in our hyperbolic example. So we turn to a variant of that method, by subtracting from the distribution of distances between images in the closed model, the similar distribution in Friedmann's open model. The result is a plot with much oscillation in small scales, modulated by a long wavelength quasi-sinusoidal pattern.

We briefly summarize the idea of cosmological models with compact, flat spatial sections. It has been suggested that, because of the COBE satellite's maps of the microwave background, such models cannot be small in the sense of Ellis, and hence are no longer interesting. Here we use Lehoucq et al.'s method of cosmic crystallography to show that these models are physically meaningful even if the size of the spatial sections is of the same order of magnitude as the radius of the observational horizon.