For the past decade there has been a considerable debate about the existence of chaos in the mixmaster cosmological model. The debate has been hampered by the coordinate, or observer, dependence of standard chaotic indicators such as Lyapunov exponents. Here we use coordinate-independent, fractal methods to show the mixmaster universe is indeed chaotic.[S0031-9007 (97)02322-3] PACS numbers: 98.80.Hw, 05.45. + bThe origin of the Universe and the fate of collapsing stars are two of the great mysteries in nature. In general relativity without exotic matter, the singularity theorems of Hawking and Penrose [1] argue that the gravitational collapse of very massive stars ends singular and that the Universe was born singular. These singular settings force gravity to face quantum mechanics. As well as exposing the fundamental laws of physics, the singular cores of black holes and the origin of the cosmos draw deep connections to the laws of thermodynamics and, as we will discuss here, to chaos.Earlier, Khalatnikov and Lifshitz [2] argued that singular solutions were the exception rather than the rule, putting them at odds with the singularity theorems. Their argument was that deformations in spacetime would be amplified during collapse and, consequently, would fight the formation of a singularity. This implied that the known symmetric singular solutions were atypical. The conflict was resolved when they realized the singularity in a collapsed star could be chaotic [3]. They conjectured that a generic singularity drives spacetime to churn and oscillate chaotically. Independently, Misner [4] suggested a chaotic approach to an early universe singularity. In his mixmaster universe, the different directions in three-space alternate in cycles of anisotropic collapse and expansion. A popular account of these developments may be found in Thorne's recent book [5].While the emergence of chaos helped our understanding of singularities in general relativity, building a resilient theory of relativistic chaos has become a task of its own. A debate has raged over whether or not the mixmaster universe is chaotic. Studies of the mixmaster dynamics using both approximate maps [6,7] and numerical integrations [8,9] have each yielded conflicting results as to the existence of positive Lyapunov exponents-a standard chaotic indicator. Finally, it was realized [9,10] that Lyapunov exponents are coordinate dependent and the conflicting results were a consequence of the different coordinate systems. In short, Lyapunov exponents are not reliable indicators of chaos in general relativity. Using a different approach, it was shown that the mixmaster equations fail the Painlevé test [11]. This suggests that the mixmaster may be chaotic, but the Painlevé test is also inconclusive.A detailed review of the mixmaster debate can be found in Ref. [12].In this Letter we show that the mixmaster universe is indeed chaotic by using coordinate independent, fractal methods. A fractal set of self-similar universes is uncovered by numerically solving Einstein's equ...