Clustering effects on the Dark Siren determination of $H_0$: A simulation study

The universe is expanding, stretching the space between galaxies like a rising loaf of bread. For nearly a century, astronomers have tried to measure the exact speed of this expansion, a value known as the Hubble constant. This number is fundamental to our understanding of how the cosmos began and how it will end. However, a stubborn disagreement has emerged in modern cosmology. When scientists look back at the early universe using the afterglow of the Big Bang, they calculate one speed. When they look at the nearby universe using exploding stars and other cosmic beacons, they calculate a different, faster speed. This conflict, known as the Hubble tension, suggests that either our measurements are flawed or our understanding of the universe's physics is incomplete. To solve this, researchers need a new, independent way to measure cosmic expansion that does not rely on the same assumptions as previous methods.

Gravitational waves offer a promising path forward. These are ripples in the fabric of space-time caused by violent cosmic events, such as the collision of two black holes. When these waves reach Earth, detectors can measure the distance to the source with high precision, essentially acting as a cosmic ruler. The problem is that these detectors cannot pinpoint exactly where the waves came from on the sky, and for black hole collisions, there is no accompanying flash of light to help identify the host galaxy. Without knowing the distance and the speed at which the host galaxy is moving away from us, astronomers cannot calculate the expansion rate. This is where the "dark siren" method comes in. Instead of finding the single galaxy that hosted the event, scientists look at all the galaxies within the large, fuzzy area where the waves originated. By statistically analyzing the distances of these many potential host galaxies, they can infer the expansion rate of the universe.

A new study by Marios Kalomenopoulos and his colleagues investigates how well this statistical method works when faced with the messy reality of the universe. In the real world, galaxies are not scattered randomly; they clump together in groups and clusters, separated by vast empty voids. Furthermore, our telescopes cannot see every galaxy; they are limited by brightness, meaning many faint, distant galaxies are missing from our catalogs. The researchers wanted to know if these two factors—clustering and missing data—would ruin the dark siren method or if it could still work. To find out, they did not wait for new telescope data. Instead, they built a sophisticated computer simulation of the universe, complete with realistic clusters of galaxies and the physical laws that govern them. They then simulated thousands of gravitational wave events within this virtual universe to see how accurately the method could recover the true expansion rate.

The team discovered that the clumpy nature of the universe is actually a help, not a hindrance. When galaxies are clustered together, the statistical method converges on the correct answer much faster. In their simulations, the researchers found that they could recover the true expansion rate with just forty gravitational wave events, provided those events occurred in a universe where galaxies were grouped naturally. In contrast, if they assumed galaxies were spread out perfectly evenly, as some older models did, even one hundred events were not enough to pin down the answer. The clustering provides a structural context that allows the math to work more efficiently, effectively compensating for the lack of a single, identified host galaxy.

Perhaps the most surprising finding concerns the completeness of galaxy catalogs. The researchers tested what happens when the galaxy lists are incomplete, mimicking real-world scenarios where telescopes miss faint galaxies. They found that the method is remarkably robust. Even when the catalog was missing seventy-five percent of the galaxies, leaving only twenty-five percent of the true population, the results remained competitive with a fully complete list. This suggests that astronomers do not need to wait for perfect, all-encompassing maps of the entire universe to make progress. The presence of clustering means that even a partial view of the cosmic web contains enough information to solve the puzzle. The study indicates that improving the sensitivity of gravitational wave detectors to find more events is a more effective strategy than trying to build more complete galaxy catalogs.

The researchers also examined other potential pitfalls, such as the specific weights given to galaxies based on their mass or brightness, and the size of the error margins in the gravitational wave measurements. They found that adding complex physical weights to the galaxies did not significantly improve the results beyond what the clustering already provided. They also confirmed that the method is not biased by the range of values the scientists initially assume the expansion rate might have. The most critical factor for precision turned out to be the distance of the events; collisions happening closer to Earth provided much sharper results than those happening at the edge of the observable universe.

Ultimately, this work provides a clear roadmap for future cosmology. It demonstrates that the "dark siren" technique is a viable and powerful tool for resolving the Hubble tension, even with imperfect data. The study reassures the scientific community that the natural clustering of galaxies, which was once seen as a complication, is actually a feature that strengthens the measurement. As gravitational wave detectors like the Laser Interferometer Gravitational-Wave Observatory and its international partners come online and detect thousands of black hole collisions in the coming years, this method will likely become a cornerstone of our understanding of the universe's expansion. The path to solving the Hubble tension does not require a perfect map of every galaxy, but rather a steady stream of gravitational wave signals and a statistical approach that embraces the clumpy, structured reality of the cosmos.