⚛️ nuclear theory

Signed Sound-Speed Deformations of BSk24-Anchored Barotropes: Neutron-Star Response

This paper quantifies how localized, thermodynamically consistent deformations of the squared sound speed in a BSk24-anchored neutron-star barotropic model induce mass-dependent, sign-asymmetric changes in stellar structure and tidal deformability, demonstrating that such phenomenological interventions propagate through reconstruction to alter central states and create distinct core or shell features without implying specific microscopic phase transitions.

Original authors: Ioannis Papathanasiou

Published 2026-08-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Ioannis Papathanasiou

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Deep in the heart of a neutron star, matter is squeezed so tightly that a single teaspoon of it would weigh as much as a mountain. These stellar remnants are the densest objects in the universe, and their internal structure is governed by a simple but profound rule: how much pressure is needed to hold up a given amount of energy. Physicists call this relationship the equation of state. It acts like a blueprint for the star, determining how big it is, how heavy it can get before collapsing, and how it ripples when shaken by a cosmic collision. One key measure of this stiffness is the speed at which sound would travel through the star's core. In the extreme environment of a neutron star, this sound speed is not just a number; it is a direct indicator of how the pressure changes as the density increases. If scientists can understand how a small, localized change in this sound speed ripples through the entire star, they can better interpret the signals from real neutron stars observed by telescopes and gravitational wave detectors.

A new study by Ioannis Papathanasiou at the Aristotle University of Thessaloniki investigates exactly this question. The researcher asked what would happen if the speed of sound inside a neutron star were suddenly sped up or slowed down in a specific region, while keeping everything else about the star exactly the same. To do this, he used a highly detailed, mathematically precise model of neutron star matter known as BSk24. This model serves as a reliable baseline, representing the best current understanding of how matter behaves under such extreme conditions without any strange, unknown features. The study did not try to discover new particles or prove that a specific phase transition exists. Instead, it treated the sound speed as a dial that could be turned up or down in a controlled, localized bump, and then watched how the entire star responded.

The process began by taking the standard model and adding a smooth, bell-shaped bump to the sound speed at a specific energy density. This bump could be positive, making the sound travel faster, or negative, making it slower. The researcher then checked to ensure that this change did not break the fundamental laws of physics, such as allowing sound to travel faster than light or causing the pressure to drop below zero. Once a valid change was confirmed, the computer recalculated the entire internal structure of the star. Because the pressure at any point depends on the history of the sound speed from the center outward, a local change in sound speed creates a permanent shift in the pressure throughout the star. This shift forces the star to rearrange itself to maintain its weight, moving its center of mass and changing its overall size.

The results showed that even a small, localized change in the sound speed has a massive, non-local effect on the star. When the sound speed was increased, the star became stiffer and expanded, becoming larger in radius and more resistant to tidal forces. When the sound speed was decreased, the star softened and shrank. However, the response was not perfectly symmetrical. A positive change that stiffened the star produced a different magnitude of effect than a negative change of the same size that softened it. For a typical neutron star with a mass of 1.4 times that of our Sun, increasing the sound speed by a specific amount increased the star's tidal deformability—a measure of how easily it can be stretched—by about 35 percent. Decreasing the sound speed by the same amount reduced this deformability by about 46 percent. For heavier stars of 2.0 solar masses, these differences became even more dramatic, with the stiffened star showing an 87 percent increase in deformability and the softened star showing a 63 percent decrease.

The study also revealed where these changes actually sit inside the star. For the stiffened stars, the region where the sound speed was altered remained connected to the very center of the star, effectively becoming part of the core. In contrast, for the softened stars, the altered region was pushed outward, forming a shell that did not touch the center, especially in the most massive stars. This means that the same physical change in the material properties can look completely different depending on the mass of the star and whether the change made the matter stiffer or softer. The researchers found that these effects are driven by the way the pressure integrates the sound speed changes from the center to the surface, creating a "memory" of the change that persists even after the local bump in sound speed has faded away.

It is important to note that these findings are specific to the mathematical model used and the particular shape of the sound-speed change applied. The study does not claim that real neutron stars actually contain these specific bumps in their sound speed, nor does it prove that a specific type of particle or phase transition is responsible for them. Instead, it provides a clear map of how a neutron star would react if such a feature existed. By holding the baseline model and the shape of the change constant, the study isolates the pure response of the star to the sign and size of the deformation. This helps astronomers understand that if they observe a neutron star behaving in a certain way, they cannot simply look at a single feature in the sound speed to explain it; they must consider how that feature interacts with the star's total mass and its global structure. The work confirms that the interior of a neutron star is a tightly coupled system where a local change in one property inevitably reshapes the entire object.

Technical Summary: Signed Sound-Speed Deformations of BSk24-Anchored Barotropes

Problem Statement
The internal structure of neutron stars is governed by the equation of state (EoS), specifically the relationship between pressure $P$ and total energy density $\epsilon$ . The local stiffness of this relationship is defined by the squared sound speed, $c_s^2 = dP/d\epsilon$ . While localized structures in $c_s^2$ (such as bumps, steps, or kinks) have been studied in various contexts to infer microscopic physics or phase transitions, it remains unclear how a single, signed feature propagates through the stellar structure when all other design choices are held constant. Specifically, the study addresses the need to quantify the response to positive and negative localized sound-speed changes within a thermodynamically consistent framework, ensuring that modifications to $dP/d\epsilon$ are integrated into a valid barotrope over its complete declared domain.

Methodology
The study employs the unified BSk24 equation of state as a fixed analytic baseline, utilizing its published Appendix-C representation to ensure explicit derivatives and domain limits. The methodology involves the following steps:

Deformation Geometry: A localized deformation is introduced to the squared sound speed $c_s^2$ using a Gaussian profile multiplied by a quintic smootherstep activation ramp. This ensures the deformation is exactly zero below a match energy density ( $\epsilon_{match} = 80 \text{ MeV fm}^{-3}$ ) and possesses continuous first and second derivatives. The only variable parameter is the dimensionless signed amplitude $A$ .
Thermodynamic Reconstruction: For each proposed $A$ , the pressure $P_A(\epsilon)$ is reconstructed by integrating the modified sound speed. The reconstruction is anchored at $\epsilon_{match}$ to preserve the baseline state below this point.
Admission Criteria: Before stellar calculations, raw proposals are strictly vetted. Cases violating $P > 0$ , $0 < c_s^2 \leq 1$ (causality), or $\epsilon > 0$ anywhere in the retained domain are rejected. No repairs or extrapolations are performed on failed cases.
Stellar Modeling: Accepted barotropes are evolved using the Tolman–Oppenheimer–Volkoff (TOV) equations and tidal deformation equations. Observables (Radius $R$ , Tidal Deformability $\Lambda$ , Love number $k_2$ , central pressure $P_c$ , and central energy density $\epsilon_c$ ) are compared at fixed gravitational masses ($0.2$ to $2.0 M_\odot$ ).
Radial Support Analysis: The study tracks the radial location of the deformation feature within the star, classifying it as "core-connected" (starting at $r=0$ ) or an "off-center shell" based on whether the central energy density $\epsilon_c$ intersects the deformation interval.

Key Contributions

Controlled Comparative Framework: The paper isolates the effect of a signed sound-speed feature by holding the analytic baseline, deformation geometry, and comparison protocol fixed, varying only the amplitude $A$ .
Strict Thermodynamic Consistency: It implements a "fail-closed" admission protocol where any violation of causality or positivity in the raw proposal leads to immediate rejection, ensuring all reported results are physically valid effective barotropes.
Exact Zero-Amplitude Control: The study utilizes an exact $A=0$ control that is array-identical to the analytic baseline, providing a precise reference for quantifying deviations.
Radial Topology Diagnosis: It introduces a diagnostic to map how a fixed density-space feature shifts relative to the stellar center as mass and stiffness change, distinguishing between core-connected and shell-like occurrences.

Results

Non-Local Pressure Memory: A localized change in $c_s^2$ results in a persistent pressure offset ( $\Delta P$ ) that integrates over the domain. This "pressure memory" causes a permanent shift in the equation of state even after the local sound-speed feature decays.
Sign-Asymmetric Response: The stellar response is not mirror-symmetric for equal magnitudes of positive and negative $A$ $A$ .
- Stiffening ( $A > 0$ ): Increases pressure support, allowing a fixed mass to be supported at lower central densities and pressures. This results in larger radii and significantly increased tidal deformability ( $\Lambda$ ).
- Softening ( $A < 0$ ): Reduces pressure support, requiring higher central densities and pressures to support the same mass. This leads to smaller radii and decreased $\Lambda$ .
Quantitative Impact: At $1.4 M_\odot$ , a deformation of $A = -0.8$ reduces $\Lambda$ by $46.02\%$ , while $A = +0.8$ increases it by $35.33\%$ . At $2.0 M_\odot$ , the asymmetry is more pronounced: $A = -0.8$ yields a $-63.52\%$ change, while $A = +0.8$ yields a $+87.17\%$ change.
Radial Topology Shifts:
- For stiffening ( $A = +2$ ), the feature remains "core-connected" across the mass range, as the lowered central density keeps the feature near the center.
- For softening ( $A = -1.6$ ), the required increase in central density pushes the feature outward. At high masses ( $2.0 M_\odot$ ), the feature exists as an off-center shell ( $r/R \approx 0.84\text{--}0.94$ ), distinct from the core.
Admissibility Limits: The study identifies strict bounds for $A$ ( $A_{min} \approx -1.97$ , $A_{max} \approx 4.37$ ) where causality or positivity is violated. Selected reporting endpoints ( $A = \pm 0.8$ , $A = -1.6$ , $A = +2$ ) lie safely within these bounds.

Significance and Claims
The paper explicitly frames its findings as characterizing phenomenological BSk24-anchored effective barotropes, not as predictions of the BSk24 functional itself, nor as evidence for a specific microscopic composition or phase transition.

Baseline Dependence: The results demonstrate that even a localized change in $dP/d\epsilon$ has nonlocal barotropic and stellar consequences due to pressure integration and central-state relocation. These findings are specific to the selected BSk24 baseline and the chosen geometry.
Scope of Results: The study characterizes the response of these specific effective barotropes but explicitly does not establish a universal response law independent of the baseline, nor does it infer microscopic constituents or identify phase transitions. The results are strictly limited to the defined finite-domain surface convention and the assumptions of general relativity, spherical symmetry, and zero temperature.
Limitations: The author emphasizes that the numerical values are specific to the selected BSk24 baseline and the ramped-Gaussian geometry. The study does not infer microscopic constituents, identify phase transitions, or establish a universal response law independent of the baseline. The results are strictly limited to the defined finite-domain surface convention and the assumptions of general relativity, spherical symmetry, and zero temperature.

Technical Summary: Signed Sound-Speed Deformations of BSk24-Anchored Barotropes

More like this