System-size dependence of bottomonium suppression from Pb-Pb to light-ion collisions

In the heart of the universe's most extreme environments, matter behaves in ways that defy our everyday experience. When protons and neutrons are smashed together at nearly the speed of light, they melt into a primordial soup known as the quark-gluon plasma. This state of matter existed only microseconds after the Big Bang, a time when the fundamental building blocks of the universe were free to roam rather than being locked inside atomic nuclei. Physicists recreate this fleeting state in massive particle accelerators, firing heavy ions into one another to generate temperatures trillions of degrees hotter than the center of the sun. One of the most reliable ways to detect this invisible fire is by watching how it affects heavy particles called bottomonium. These particles are like tiny, dense anchors made of a heavy quark and its antimatter partner. In normal conditions, they hold together tightly, but if they pass through the quark-gluon plasma, the intense heat and pressure can pull them apart. By measuring how many of these particles survive a collision, scientists can map the temperature and density of the plasma they traversed.

For decades, researchers have studied this phenomenon in collisions between lead atoms, which are large and create a substantial, long-lasting fireball. However, a new question has emerged: does this same physics apply when the collision involves much smaller atoms, like oxygen or neon? These lighter collisions create a smaller, shorter-lived version of the quark-gluon plasma, offering a controlled way to test if the rules governing this exotic matter change with the size of the system. A team of physicists at the Technical University of Munich has now tackled this question by running sophisticated computer simulations that track the fate of bottomonium particles as they move through the plasma created in lead, oxygen, and neon collisions. Their work suggests that the microscopic rules describing how the plasma breaks apart these heavy particles remain consistent, regardless of whether the collision involves a massive lead nucleus or a much smaller neon one.

The researchers used a framework that treats the bottomonium particle not as a static object, but as a quantum system interacting with a chaotic environment. In this view, the particle is constantly jostled by the surrounding plasma, a process that can be described by two key numbers that characterize how the medium resists and disturbs the particle's motion. These numbers were previously determined by matching simulations to data from lead-lead collisions. The team's goal was to see if they could use those exact same numbers to predict what would happen in the smaller oxygen and neon collisions, without adjusting the settings to fit the new data. They employed a powerful software tool called QTraj to solve complex equations that describe the real-time evolution of the bottomonium as it travels through the plasma. This simulation accounts for the fact that the plasma is not uniform; it expands and cools rapidly, and the bottomonium particle experiences different temperatures and densities as it moves through it.

To make these predictions, the team first modeled the creation and expansion of the plasma for each type of collision. For the oxygen and neon collisions, which occurred at the Large Hadron Collider at an energy of 5.36 TeV, they calculated how the initial energy density translates into a temperature and how long the fireball lasts. They found that while the neon collisions produce a slightly hotter and denser initial state than oxygen, the overall lifetime of the plasma is still very short. They then ran their simulations to see how many of the different types of bottomonium states would survive. There are several versions of these particles, ranging from the most tightly bound ground state to more loosely bound excited states. The theory predicts a specific hierarchy: the loosely bound states should dissolve easily in the heat, while the tightly bound ones should survive in greater numbers. This pattern of sequential melting is a hallmark of the quark-gluon plasma.

The results of the simulation showed a striking agreement with recent experimental measurements. When the team compared their predictions for oxygen collisions to data collected by the CMS and LHCb experiments, the simulated survival rates matched the observed pattern of suppression. The model correctly predicted that the excited states would be significantly more suppressed than the ground state, reproducing the sequential hierarchy seen in the data. This was achieved without changing the fundamental transport coefficients that describe the interaction between the particle and the plasma. The fact that the same microscopic rules that worked for the massive lead collisions also worked for the much smaller oxygen collisions suggests a deep consistency in the nature of the quark-gluon plasma. It implies that the plasma behaves as a universal medium, where the way it disrupts heavy particles depends on its temperature and density, not on the size of the system that created it.

The researchers also tested how sensitive their results were to the exact conditions of the simulation, particularly the temperature at which they stopped the calculation. They found that lowering this final temperature led to stronger suppression, as the particles spent more time in the hot medium. However, even with these variations, the relative difference between the oxygen and neon results remained small. The simulations predicted that the suppression in neon collisions would be only slightly different from that in oxygen, roughly 11 to 13 percent smaller for one specific ratio of survival rates. This prediction stands in contrast to the central values of the latest experimental data, which hint at a larger difference between the two systems. The authors note that the experimental uncertainty on the neon measurement is currently quite large, so the data does not yet definitively rule out their simulation. If future, more precise measurements confirm a larger gap between oxygen and neon, it would suggest that the current models of how the plasma evolves are missing something, perhaps related to the initial fluctuations of the collision or other nuclear effects not yet included in the calculation.

Ultimately, this study reinforces the idea that the quark-gluon plasma is a well-defined state of matter with consistent properties across different collision sizes. By extending the analysis from lead to lighter ions, the researchers demonstrated that the microscopic description of bottomonium suppression is robust. The same temperature-dependent parameters that explain the behavior of heavy particles in the largest collisions also explain their behavior in the smallest ones. This consistency supports the interpretation that a quark-gluon plasma is indeed formed in these light-ion collisions, and that its ability to dissolve heavy particles is governed by universal laws. While the current data leaves room for refinement, the work provides a strong foundation for understanding the fundamental dynamics of the early universe, showing that the rules of this extreme matter hold true even when the stage is scaled down.