Extracting the Central Sound Speed of Neutron Stars through Intrinsic Stellar Structure Scaling with Reduced Equation-of-State Model Dependence

Deep within the heart of a neutron star, matter is crushed to densities that cannot be replicated anywhere else in the universe. These stellar remnants are the collapsed cores of massive stars, packing more mass than our Sun into a sphere only about the size of a city. Because they are so dense, the atoms that make up ordinary matter are squeezed so tightly that protons and electrons merge into neutrons, creating a fluid of pure nuclear matter. Physicists have long wanted to know how this fluid behaves under such extreme pressure. Specifically, they want to know how fast sound waves would travel through the center of such a star. This speed, known as the sound speed, acts as a measure of how stiff or rigid the matter is. If the matter is very stiff, sound travels quickly; if it is soft, sound travels slowly. Understanding this stiffness is crucial because it determines how large a neutron star can be and how much mass it can hold before collapsing into a black hole.

However, figuring out this stiffness from the outside is incredibly difficult. Astronomers can measure the mass and size of neutron stars, and they can observe how these stars ripple when they collide, but they cannot see inside. To interpret these observations, scientists must rely on theoretical models that describe how nuclear matter behaves. The problem is that many different models can fit the same observational data while predicting completely different conditions inside the star. It is like trying to guess the ingredients of a cake by only tasting the frosting; different recipes could produce the same sweet topping while having very different batters inside. This ambiguity has left a gap in our understanding: we can describe the surface of a neutron star well, but the core remains a mystery that depends heavily on which theoretical recipe we choose to use.

A team of researchers has now developed a new way to look past these theoretical differences to find what the core actually tells us. Instead of trying to pick the single best model for nuclear matter, they asked a different question: if we look at the fundamental laws that govern how a star holds itself together, is there a piece of information about the core that remains the same regardless of which model we use? To answer this, they created a vast collection of different theoretical models. They took a standard framework for describing nuclear forces and systematically varied one specific interaction term, creating six distinct families of models. For each family, they generated one thousand different variations, ensuring that every single one satisfied all known constraints from both laboratory experiments on Earth and observations of neutron stars in space.

The results of this exercise were striking. Even though all these thousands of models agreed with every known observation, they predicted wildly different conditions in the center of a neutron star. Some models suggested the core was extremely stiff, while others suggested it was much softer. Some predicted that the speed of sound would rise to a peak in the middle of the star before falling again, while others showed a steady increase. This confirmed that current observations alone are not enough to tell us which model is correct; the data simply cannot distinguish between these very different internal realities. The researchers realized that if they continued to rely on picking one specific model, their conclusions about the core would always be biased by the assumptions built into that model.

To break this deadlock, the team turned to the mathematical equations that describe the structure of a star itself. These equations, known as the Tolman-Oppenheimer-Volkoff equations, link the pressure and density at the center of a star to its total mass and radius. The researchers analyzed these equations to find a specific relationship that connects the central sound speed to two measurable quantities: the ratio of pressure to energy density at the center, and how the star's total mass changes as the central density changes. They derived a new mathematical tool, an estimator, that uses these two quantities to predict the central sound speed without needing to assume a specific model for the nuclear matter.

When they applied this tool to their vast collection of models, a clear pattern emerged. Despite the fact that the models had different internal rules and predicted different behaviors at high densities, they all lined up along a single, tight relationship. The estimator they created was able to predict the central sound speed for every single model with high accuracy. This meant that the relationship between the star's global structure and its core properties was universal, cutting through the noise of the different theoretical assumptions. By finding the point where all these different model families agreed, the researchers could extract a value for the central sound speed that was far less dependent on any single theory.

The final result points to a specific value for the speed of sound squared in the center of a typical neutron star with a mass of 1.4 times that of our Sun. The researchers found this value to be approximately 0.523, with a small range of uncertainty. This number is significantly higher than the value expected for a fluid that behaves like light, suggesting that the matter in the core is quite stiff. They also determined the central pressure-to-density ratio to be about 0.173. Using these numbers, they calculated that the radius of such a star should be approximately 11.73 kilometers. When they repeated the analysis while excluding the most extreme model family, the results shifted only slightly, to a sound speed of 0.496 and a radius of 11.92 kilometers. This consistency suggests that the finding is robust.

The study also revealed something surprising about the internal structure of these stars. In many of the models, the speed of sound does not simply increase steadily toward the center. Instead, it rises to a peak and then drops slightly before the star's surface. This peak is caused by a competition between different forces inside the nuclear matter. The researchers showed that this behavior is a natural consequence of the way nuclear forces interact at high densities, without needing to invoke exotic new types of particles. The fact that this peak appears in models that satisfy all current observations suggests that such a feature might actually exist in real neutron stars.

By focusing on the shared behavior of the star's structure rather than the specific details of the nuclear models, the researchers have provided a new way to probe the deep interior of these cosmic objects. Their work demonstrates that even when we cannot uniquely determine the laws of physics at extreme densities, we can still extract reliable information about the state of matter by looking for patterns that persist across all possible theories. This approach offers a clearer path forward, allowing astronomers to infer the properties of the neutron star core with greater confidence, independent of the specific theoretical framework used to describe it. The findings suggest that the cores of these dense stars are stiffer than previously thought, with a sound speed that is roughly half the speed of light, providing a concrete benchmark for future observations and theoretical developments.