$r$-Process as Production of Lanthanides from Compact Object Mergers

The universe is a vast chemical factory, constantly forging the heavy elements that make up our world, from the gold in jewelry to the iron in our blood. While lighter elements like carbon and oxygen are born in the fiery hearts of stars, the heaviest elements in the periodic table require a more violent and exotic origin. These heavy atoms are created through a rapid sequence of neutron captures, a process known as the rapid neutron-capture process, or r-process. For decades, astronomers have debated where exactly this cosmic alchemy takes place. The leading suspects are the catastrophic collisions of compact stellar remnants: binary neutron star mergers, where two ultra-dense stars crash together, and black hole-neutron star mergers, where a black hole consumes a neutron star. Understanding which of these events dominates the production of heavy elements, and how they have enriched our galaxy over billions of years, is essential to explaining why the stars and planets around us contain the specific mix of materials they do.

A team of researchers has now constructed a comprehensive digital simulation to trace the history of these heavy elements from the moment stars are born to the present day. By combining three powerful computational tools, they modeled the life cycles of billions of binary star systems, calculated the nuclear reactions that occur when these systems collide, and tracked how the resulting debris spreads through the Milky Way. Their goal was to determine the relative contributions of binary neutron star mergers versus black hole-neutron star mergers to the galaxy's inventory of heavy elements, particularly the lanthanides, a group of rare earth metals that includes europium. The study did not rely on a single set of assumptions; instead, it tested five different theoretical models for the internal structure of neutron stars and three different spin rates for the black holes involved, ensuring the results were robust against the uncertainties of current physics.

The researchers began by simulating the evolution of one million binary star systems, tracking how they interact, lose energy, and eventually merge. They found that while most supernova explosions happen quickly, within the first 50 million years, the collisions of compact objects are more spread out over time. About half of the binary neutron star mergers occur within the first 100 million years, but a significant portion happens much later. When these collisions occur, they eject material into space at tremendous speeds. The team discovered that the specific internal structure of the neutron stars, defined by their equation of state, did not significantly change the overall pattern of heavy elements produced. However, the type of merger mattered greatly depending on the mass of the elements. For elements lighter than a specific threshold, the material ejected from the swirling disk of gas surrounding the new object was the primary source. For the heaviest elements, the material thrown out directly during the violent collision itself was the dominant contributor.

A key finding of the study concerns the role of the black hole's spin. In a black hole-neutron star merger, a rapidly spinning black hole can pull the neutron star apart before it is swallowed, creating a much larger amount of ejected material. The simulations showed that higher spin rates lead to more ejected mass, which in turn produces more heavy elements. Despite this, the researchers found that the spin of the black hole did not significantly alter the specific mix or "fingerprint" of the elements created; it simply changed the total quantity. When comparing the two types of mergers, the study revealed a distinct division of labor. Binary neutron star mergers were found to be the superior producers of the heaviest elements, particularly those with atomic masses greater than 130. In contrast, black hole-neutron star mergers were less efficient at creating these heavyweights, though they still contributed to the lighter end of the spectrum.

The researchers then placed these findings into a model of the Milky Way's chemical evolution to see how these events shaped the galaxy over time. They treated the galaxy as a collection of independent gas clouds, each capable of being enriched by a single merger event. Because binary neutron star mergers happen frequently enough and produce a consistent amount of heavy elements, they act as a steady, reliable source that gradually builds up the abundance of elements like europium in the later stages of the galaxy's life. Black hole-neutron star mergers, being rarer and more variable in their output, create a patchy distribution of heavy elements in the early universe, leading to significant differences in chemical composition from one star to another. The simulations suggest that while black hole-neutron star mergers play a role, binary neutron star mergers are the dominant force responsible for the heavy r-process elements found in the Milky Way today.

However, the study also highlighted a gap between their simulations and the actual observations of our solar system. While binary neutron star mergers successfully explained the abundance of the heaviest elements, they fell short in producing enough of the lighter heavy elements, specifically those around the second peak of the periodic table. This discrepancy suggests that other cosmic events, such as certain types of supernova explosions or perhaps the outflows from rapidly spinning neutron stars known as magnetars, must also contribute significantly to the galaxy's chemical makeup. The work confirms that the collisions of compact objects are indeed the primary engines for creating the universe's heaviest ingredients, but it also reminds us that the full story of cosmic chemistry likely involves a chorus of different astrophysical events working in concert.