Ellipsoidal Collapse Model

How do the vast, intricate structures of our universe, from galaxies to the cosmic web, come into being? The answer lies in gravity's relentless work on the tiny density fluctuations of the early cosmos. A first attempt to model this process, the spherical collapse model, provides a useful but incomplete picture by assuming these fluctuations are perfectly spherical. This simplification overlooks a crucial aspect of nature: asymmetry. Real cosmic structures are shaped by the complex pushes and pulls of their environment, leading to a collapse that is far from uniform. This article delves into the more realistic and powerful ​​ellipsoidal collapse model​​, which embraces this complexity. The first chapter, "Principles and Mechanisms," will uncover how this anisotropic collapse unfolds sequentially, from three-dimensional perturbations into sheets, filaments, and finally the halos that host galaxies. Following this, the "Applications and Interdisciplinary Connections" chapter will explore how this theory allows us to predict the properties of halos, understand their distribution across the universe, and even explain the shapes of the galaxies we see today.

The story of how the universe built its vast and intricate structures—galaxies, clusters of galaxies, and the great cosmic web that connects them—is a story of gravity. But it's not the simple story you might remember from introductory physics, of an apple falling to the Earth or a planet orbiting the Sun. It's a grander, more complex, and far more beautiful tale of cosmic evolution. To understand it, we must abandon our comfortable notions of perfection and embrace the power of asymmetry.

Physicists love spheres. They are simple, elegant, and mathematically convenient. When we first try to model the formation of a galaxy or a dark matter halo, we naturally start with the simplest possible picture: a perfectly spherical, uniform patch of matter that is slightly denser than the cosmic average. In an expanding universe, this "top-hat" overdensity will fight against the cosmic expansion. If it's dense enough, its own gravity will win. It will slow down, stop expanding, turn around, and collapse into a compact, bound object. This is the essence of the ​​spherical collapse model​​. It’s a wonderfully useful starting point and gives us a crucial number: the ​​critical overdensity​​, $\delta_c \approx 1.686$. This number is a cosmic deadline. In the language of linear theory—a way of tracking the growth of the initial, tiny density fluctuations—any region whose density contrast was destined to exceed $1.686$ has, by today, collapsed.

But the universe, in its magnificent complexity, doesn't seem to share our fondness for perfect spheres. Imagine a truly colossal stellar explosion. If it were perfectly, impossibly symmetric—expanding outwards in a perfect sphere—it would be utterly silent in the realm of gravitational waves. The universe's fabric would remain undisturbed. Why? Because the generation of gravitational waves requires a changing asymmetry, a changing mass quadrupole moment. A perfectly symmetric motion has none. Any slight deviation from a perfect sphere, however, and the explosion will send ripples through spacetime.

This is a profound clue. Nature's most dramatic events are often messy and asymmetric. The collapse of matter to form cosmic structures is no different. The initial seeds of structure were not perfect spherical "top-hats." They were irregular, lumpy patches of slightly higher density, stretched and squeezed by the gravitational pull of their neighbors. To truly understand how a galaxy forms, we must move beyond the sphere and embrace the ellipsoid.

Let's imagine one of these slightly lumpy, overdense regions in the early universe. The best way to approximate its shape isn't a sphere, but a triaxial ellipsoid—something like a slightly squashed and elongated potato. This object has three principal axes of different lengths. Now, gravity gets to work. But it doesn’t pull equally in all directions. The pull is strongest along the shortest axis and weakest along the longest. The result is not a uniform, spherical implosion. It is an ​​anisotropic collapse​​, a sequence of events, a symphony in three movements.

To describe this process, cosmologists use a brilliantly effective tool called the ​​Zel'dovich approximation​​. It describes the motion of matter not as a uniform contraction, but as a displacement field acting on the initial grid of particles. For an ellipsoidal region, this field can be characterized by three crucial numbers, the eigenvalues $\lambda_1, \lambda_2, \lambda_3$ of the initial deformation tensor, which we can order as $\lambda_1 > \lambda_2 > \lambda_3$. You can think of these eigenvalues as representing the "strength of the gravitational squeeze" along each of the three principal axes. The largest eigenvalue, $\lambda_1$, corresponds to the strongest compression and thus the shortest axis.

The collapse then unfolds sequentially:

1. ​​The Pancake:​​ The matter collapses first along the direction of the strongest gravitational pull—the axis corresponding to $\lambda_1$. The initial ellipsoid is squashed into a two-dimensional, sheet-like structure. Cosmologists whimsically call this a "pancake." At the exact moment this first collapse happens, the other two axes are still finite in size. The Zel'dovich approximation tells us that the ratio of their lengths is determined purely by the initial eigenvalues: the ratio of the longest axis to the intermediate one is $(\lambda_1 - \lambda_3) / (\lambda_1 - \lambda_2)$. The initial shape dictates the shape of the resulting pancake.

2. ​​The Filament:​​ The pancake doesn't live forever. Gravity continues its work, now primarily within the plane of the sheet. The next collapse happens along the second-strongest direction of compression, the axis corresponding to $\lambda_2$. This squeezes the pancake into a one-dimensional, thread-like structure—a "filament." The time delay between the pancake and filament formation is a direct consequence of the initial anisotropy. For a small initial difference between the collapse strengths, the delay $\Delta t_{23}$ between the second and third axis collapses is proportional to the initial anisotropy parameter $\epsilon_a$. A more anisotropic initial state leads to a more staggered collapse.

3. ​​The Halo:​​ Finally, the matter along the filament collapses along the last remaining direction, the one with the weakest initial compression (smallest eigenvalue $\lambda_3$). This crunches the filament down into a compact, roughly spherical, gravitationally bound object—what we call a ​​dark matter halo​​, the cradle in which galaxies are born. The entire process, from the initial perturbation to the final halo, is complete only when collapse has occurred along all three axes, a moment determined by the smallest eigenvalue.

This sequential collapse—from 3D to 2D to 1D and back to a compact 3D object—is the fundamental mechanism of the ​​ellipsoidal collapse model​​. It paints a far richer and more realistic picture of structure formation than the simple spherical model.

The spherical model gave us a simple, universal rule: collapse happens when the linearly extrapolated overdensity reaches $\delta_c \approx 1.686$. But if the collapse itself is not spherical, should the rule for when it happens be so simple?

The ellipsoidal model forces us to rethink this. If a proto-halo is already elongated, it has a "head start" on collapsing along its shortest axis. It shouldn't need to be, on average, as dense as a perfect sphere to collapse by the same time. The collapse criterion itself must depend on the object's shape.

We can build a beautifully simple analytical model to capture this. Let's define the shape of our initial ellipsoid by two parameters: ​​ellipticity​​ ($e$), which measures how elongated it is, and ​​prolateness​​ ($p$), which measures whether it's more like a cigar (prolate) or a discus (oblate). These parameters can be defined directly from the initial eigenvalues $\lambda_1, \lambda_2, \lambda_3$. The central assumption of this model is that the "main" collapse—the event that forms the final halo—occurs when the linear theory overdensity along the shortest axis reaches a certain universal threshold. By calibrating this model with the known spherical case (where $e=0$ and $p=0$), we arrive at a stunningly elegant new rule for collapse. The critical overdensity is no longer a constant, but a function of the initial shape, with the key result that for a non-spherical object the required overdensity is lower than for a perfect sphere. This confirms our intuition: elongated or flattened objects need a lower average initial overdensity to collapse than a perfect sphere does. They are "easier" to collapse. This shape-dependent collapse criterion is a major triumph of the ellipsoidal model. It explains why we see a spread in the properties of halos of the same mass—they were born from initial perturbations with different shapes. We can even use this framework to predict what the linearly-evolved overdensity would look like at key stages, like filament formation, just by knowing the initial shape ratios of the collapsing object.

This raises a final, crucial question: where does this initial shape come from? Why would a patch of the early universe be ellipsoidal rather than spherical? The answer is that no patch of the universe is an island. Every overdense region is surrounded by other over- and under-dense regions, all pulling on it gravitationally. This complex gravitational environment creates what we call a ​​tidal field​​.

Imagine our proto-halo. On one side, a giant neighboring overdensity might be pulling matter away from it, stretching it in that direction. On another side, it might be getting squeezed between two other concentrations of matter. This stretching and squeezing by the external environment is what imprints the initial ellipsoidal shape onto the proto-halo. The principal axes of the ellipsoid are simply the principal axes of the local tidal field.

The eigenvalues $\lambda_1, \lambda_2, \lambda_3$ that govern the anisotropic collapse are, in fact, a direct measure of the strength of these tidal forces. A positive eigenvalue corresponds to a stretching direction (a tidal tension), while a negative eigenvalue corresponds to a squeezing direction (a tidal compression). The effect of these tides is to modify the conditions for collapse. A weak external tidal compression, for example, can help an overdensity collapse faster than it would in isolation, effectively lowering the critical overdensity required.

This provides the final, unifying piece of the puzzle. The ellipsoidal collapse model is not just a description of an isolated object; it is a theory of how objects form within their cosmic environment. The sequential, anisotropic collapse into pancakes, filaments, and halos is the fundamental process that weaves the ​​cosmic web​​—the vast, interconnected network of structures we observe in the universe today. The great voids are the regions where tidal forces have stretched matter in all directions. The sheets and walls are the "pancakes" that collapsed along one direction. The filaments that thread between them are the second stage of collapse. And the dense, compact clusters of galaxies are the halos that have completed their collapse in all three directions, sitting at the intersections of this great cosmic network. The ellipsoidal collapse model, in all its elegance, gives us the principles and mechanisms to understand this grand cosmic architecture.

We have spent some time appreciating the subtle mechanics of ellipsoidal collapse, understanding that the universe, in its grand construction project, is not a builder of perfect spheres. The initial, minuscule anisotropies in the primordial cosmic soup—the slight pushes and pulls from neighboring regions—are not mere imperfections to be averaged away. On the contrary, gravity seizes upon these asymmetries and amplifies them, sculpting the magnificent, complex structures we observe today. Having grasped the principles, we now ask the most exciting question a physicist can ask: So what? What does this elegant piece of theory allow us to understand about the real world?

It turns out that the answer is: a great deal. Moving beyond the idealized spherical cow, the ellipsoidal collapse model serves as a master key, unlocking insights into the very anatomy of cosmic structures, their demographics across the universe, their place in the vast cosmic web, and even the luminous galaxies they host.

The most immediate and intuitive consequence of ellipsoidal collapse concerns the final shape of a dark matter halo. Imagine a slightly lumpy, overdense patch of matter in the early universe. The simple spherical model would treat it as a ball, collapsing uniformly until it virializes—settles into a stable, spherical halo. But the ellipsoidal model tells a more interesting story.

This patch is not an island; it feels the gravitational tidal forces from all the matter around it. These forces stretch and squeeze the patch, defining three principal axes. Along one axis, the collapse might be accelerated, while along another, it is slightly hindered. Gravity is a runaway process; these initial differences in the collapse speed become magnified over billions of years. The axis that collapses the fastest and earliest ends up as the shortest axis of the final, virialized object. The slowest collapsing axis becomes the longest. Thus, an initially near-spherical region naturally evolves into a triaxial ellipsoid—a shape something like a flattened potato or a football.

Our model allows us to make this quantitative. By characterizing the initial anisotropy, whether through the eigenvalues of the primordial tidal field tensor or the initial peculiar velocity field, we can directly predict the final axis ratios of the virialized halo. This is a remarkable achievement: it forges a direct, mathematical link between the "genetic code" of a perturbation at the dawn of time and the final "adult" morphology of the dark matter halo it becomes.

But the story doesn't end with a static shape. The internal dynamics—the very "soul" of the halo—are also imprinted with this history of anisotropic collapse. A halo is not a solid object; it's a swarm of dark matter particles held together by gravity, supported against further collapse by their random motions, or velocity dispersion. The ellipsoidal model predicts that this cloud of particles is not "hot" in the same way in all directions. The particles tend to have more kinetic energy—a higher velocity dispersion—along the axes that collapsed more violently and ended up shorter. This theoretical connection between the halo's geometry and its internal kinematics provides a richer, more detailed picture of these fundamental objects, a picture we can test with observations and simulations.

From the anatomy of a single halo, we now lift our gaze to the entire cosmos. One of the most fundamental predictions of any structure formation theory is the ​​halo mass function​​: a cosmic census that answers the question, "For any given mass $M$, how many dark matter halos of that mass should exist in the universe?"

The standard tool for this is the "excursion set formalism," a wonderfully intuitive picture. Imagine the smoothed density fluctuation at a point, $\delta$. As we smooth over smaller and smaller scales (corresponding to smaller masses), this value of $\delta$ performs a random walk. A halo is said to form when this random walk first crosses a critical barrier, $\delta_c$. In the spherical collapse model, this barrier is a constant value—a fixed wall that the random walker must hit.

Here, the ellipsoidal collapse model introduces a profound and beautiful complication. The barrier is no longer fixed; it moves. The logic is this: for a halo to form, it must not only collapse along its shortest axis but must eventually collapse along all three axes. For smaller mass halos, which typically form in denser environments, the external tidal forces are stronger. This shear and squashing from the outside world make it harder for the proto-halo to collapse; it needs an extra gravitational push. Therefore, the effective critical density required for collapse is higher for these smaller objects. In our random walk analogy, the wall moves farther away as the walker (representing smaller mass scales) approaches it. This "moving barrier" is a direct consequence of considering ellipsoidal, rather than spherical, dynamics.

This single idea has enormous consequences. By incorporating a moving barrier whose shape is motivated by the physics of ellipsoidal collapse, theorists developed refined mass functions, the most famous of which is the Sheth-Tormen mass function. This model, which accounts for the influence of tidal shear on the collapse threshold, provides a dramatically better fit to the results of large-scale computer simulations than the simpler spherical collapse model. It has become a cornerstone of modern cosmology, used in nearly every analysis that seeks to connect the underlying theory of dark matter to observations of galaxies and galaxy clusters.

Knowing how many halos exist is one thing; knowing where they are is another. We observe that galaxies and clusters are not sprinkled randomly through space. They trace a magnificent, filamentary structure known as the cosmic web. This means that the locations of halos are "biased" relative to the underlying distribution of matter. Denser regions of the universe are more than proportionally filled with halos.

Ellipsoidal collapse deepens our understanding of this bias. It tells us that halo formation is sensitive not just to the local density, but also to the shape of the large-scale environment. A proto-halo trying to form inside a large, filamentary structure that is stretching it in one direction will have a different fate than one forming in a massive cluster that is squeezing it from all sides. This environmental influence is, once again, the effect of the large-scale tidal field.

Using a powerful technique called the "peak-background split," we can calculate how the number of halos changes in response to an external tidal shear. This gives rise to a new kind of bias, known as ​​tidal bias​​. The ellipsoidal collapse model provides the crucial physical input, telling us precisely how the collapse threshold is modified by the presence of an external shear field. From this, we can derive the tidal bias parameters that quantify this effect. This is of immense practical importance. When we map the universe using galaxy surveys, we are not seeing the smooth matter field, but this biased distribution of halos. Understanding tidal bias is essential for correctly interpreting these maps to learn about the fundamental properties of our universe, such as the nature of dark energy and the total mass of neutrinos.

Thus far, our discussion has centered on the invisible realm of dark matter. But the ultimate goal of cosmology is to explain the universe we can see. The glorious culmination of this line of reasoning is its connection to the formation of galaxies themselves. Galaxies are not just floating in space; they are born in the gravitational potential wells at the very centers of dark matter halos. It stands to reason, then, that the properties of the parent halo should be imprinted upon the galaxy it nurtures.

And so they are. The same process of anisotropic collapse that forges a triaxial dark matter halo also orchestrates the collapse of the baryonic gas that will form an elliptical galaxy at its heart. The tidal forces that shaped the halo also shape the galaxy. The ellipsoidal collapse model, therefore, can be extended to predict the intrinsic shapes and axis ratios of elliptical galaxies, connecting them directly to the properties of the primordial tidal field from which they arose. This provides a physical origin for the observed range of galaxy shapes, transforming what was once a purely descriptive classification (Hubble's "tuning fork") into a predictive science rooted in fundamental physics.

From the shape of a single halo to the grand cosmic census, and from the clustering in the cosmic web to the visible forms of galaxies, the ellipsoidal collapse model proves its worth. It is a testament to a deep principle in physics: often, the most profound understanding comes not from ignoring the complexities and asymmetries of the world, but from embracing them and discovering the elegant laws that govern them.