[Paper Review] Cosmic crystallography: the euclidean isometries

Armando Bernui, A. F. F. Teixeira|ArXiv.org|Mar 16, 2000

Scientific Research and Discoveries4 references3 citations

TL;DR

This paper derives exact analytical probability density functions for pair separation distances in Euclidean isometries—translations, screw motions, rotations, reflections, and glide reflections—used in cosmic crystallography to probe the topology of the universe. By replacing noisy computer-simulated histograms with continuous, normalized probability densities, the study enables more precise comparison with observational data, particularly identifying sharp spikes and discontinuities linked to specific isometries.

ABSTRACT

Exact expressions for probability densities of conjugate pair separation in euclidean isometries are obtained, for the cosmic crystallography.These are the theoretical counterparts of the mean histograms arising from computer simulation of the isometries. For completeness, also the isometries with fixed points are examined, as well as the orientation reversing isometries.

Motivation & Objective

• To provide exact analytical probability density functions for pair separation distances in Euclidean isometries, serving as theoretical counterparts to simulated histograms in cosmic crystallography.
• To eliminate statistical noise in simulated histograms by replacing them with continuous, normalized probability densities derived from exact geometric and topological principles.
• To examine all Euclidean isometries—including those with fixed points and orientation-reversing types—for completeness, despite current focus on orientation-preserving, fixed-point-free isometries.
• To enable more accurate detection of cosmic topology by identifying unique spectral features (e.g., spikes, discontinuities) in the separation density functions.

Proposed method

• Derives normalized pair separation probability densities $ \mathcal{P}_g^\mathcal{B}(l) $ for each Euclidean isometry $ g $, based on the intersection of a solid ball $ \mathcal{B} $ and its image $ \mathcal{B}_g $ under $ g $.
• Uses geometric constraints: the distance $ m $ between centers $ C $ and $ C_g $ must be less than $ 2a $ for intersection, with $ a $ the ball radius.
• Applies rotational symmetry to reduce the 3D problem to 2D cross-sections, computing densities via areas of circular or annular regions in the intersection.
• For screw motions, derives $ \mathcal{P}_g^\mathcal{B}(l) $ as a function of parameters $ a, b, t, \omega $, using auxiliary variables like $ r $ and $ \rho $ to describe intersection geometry.
• For reflections and glide reflections, derives densities that diverge at minimum displacement $ l = t $, yet remain integrable with total area 1.
• Validates results by comparing analytical densities with simulated histograms, confirming discontinuities and spikes observed in prior simulations.

Experimental results

Research questions

• RQ1 How can exact analytical probability densities for pair separations be derived for all Euclidean isometries in cosmic crystallography?
• RQ2 What specific features (e.g., spikes, discontinuities) in the separation density functions correspond to different isometries, and how do they differ from simulated histograms?
• RQ3 How do fixed-point isometries (rotations, reflections) and orientation-reversing isometries (reflections, glide reflections) affect the pair separation distribution compared to free isometries?
• RQ4 Can the noise inherent in simulated histograms be eliminated by replacing them with exact continuous probability densities?

Key findings

• The probability density for pure translations is a Dirac delta function $ \mathcal{P}_t(l) = \delta(l - t) $, reflecting exact displacement by $ t $.
• For screw motions, the density $ \mathcal{P}_g^\mathcal{B}(l) $ is piecewise-defined and depends on parameters $ a, b, t, \omega $, with distinct behaviors based on intersection geometry.
• For pure reflections, the density $ \mathcal{P}_g^\mathcal{B}(l) $ is proportional to $ a^2 - (b + l/2)^2 $, peaking at $ l = 0 $ and diverging near $ l = 0 $ when $ b \to a $.
• For glide reflections, the density diverges at $ l = t $, the minimum displacement, due to a high concentration of points near the reflection plane, yet remains normalized to 1.
• Discontinuities in the density functions—such as at $ l \sim 0.7 $ or $ l = 1 $—correspond to sharp features observed in prior simulations and are robust against statistical noise.
• The transition from glide reflection to pure reflection is continuous: as $ t \to 0 $, the divergence at $ l = t $ shrinks and vanishes, yielding the reflection density.