Author: Denis Avetisyan
New research leverages Bayesian inference and observations of small, dense stars to explore the possibility of a phase transition to quark matter within neutron stars and the resulting ‘twin star’ phenomenon.

This study investigates the constraints on dense-matter equations of state derived from analyzing compact star populations and their potential for exhibiting a twin star configuration.
The extreme densities within neutron stars pose a fundamental challenge to our understanding of quantum chromodynamics and the equation of state of matter. This is addressed in ‘Bayesian Inference of Dense-Matter Equations of State from Small-Radius Compact Stars with Twin-Star Scenarios’, which investigates whether recent observations of small-radius compact stars can constrain the properties of dense matter and reveal evidence of a phase transition to quark matter. By employing Bayesian inference, this work finds a preferred phase transition density of approximately 2.7-2.8n_0, alongside a significant energy-density jump and relatively high post-transition sound speed, potentially generating a disconnected ‘twin star’ branch in the mass-radius relation. Could future multimessenger observations of these uniquely structured neutron stars definitively confirm the existence of a first-order phase transition and illuminate the nature of matter at the highest densities?
The Weight of Collapsed Stars: A Density Beyond Comprehension
Neutron stars represent an extreme state of matter, collapsing under gravity to achieve densities exceeding that of atomic nuclei – a teaspoonful would weigh billions of tons on Earth. These stellar remnants, formed from the cores of massive stars after supernova explosions, pack the mass of the Sun into a sphere roughly the size of a city. Such immense density pushes the limits of known physics; the strong nuclear force, which normally holds atomic nuclei together, is subjected to conditions never replicated in terrestrial laboratories. Consequently, the behavior of matter within a neutron star remains largely theoretical, demanding innovative approaches to understand the interplay between gravity, nuclear forces, and potentially exotic states of matter like quark-gluon plasmas or hyperons. Studying these objects offers a unique window into fundamental physics, challenging and refining existing models of matter and the universe.
Understanding the interiors of neutron stars hinges on determining a precise Equation of State (EOS), a fundamental relationship that dictates how pressure, density, and temperature interact within these incredibly dense objects. This EOS isn’t simply a mathematical curiosity; it’s the key to predicting a neutron star’s size, its gravitational pull, and even how it responds to collisions. At such extreme densities – exceeding that of atomic nuclei – matter behaves in ways not fully understood by current physics, potentially exhibiting exotic phases like quark-gluon plasma or hyperonic matter. Constructing an accurate EOS requires sophisticated theoretical models, often incorporating insights from quantum chromodynamics and nuclear physics, and then rigorously testing these predictions against observational data, such as the measurement of neutron star masses and radii. The challenge lies in the fact that directly probing the interior of a neutron star is impossible, making the EOS a crucial, yet indirect, window into the behavior of matter under the most extreme conditions in the universe.
Despite decades of research, constructing an accurate Equation of State (EOS) for matter at the extreme densities found within neutron stars remains a formidable challenge. Current theoretical models, attempting to describe the relationship between pressure, density, and temperature, consistently struggle to simultaneously match observed neutron star masses and radii. Discrepancies arise because the physics governing matter at these densities – potentially involving exotic states like quark matter or hyperons – is poorly understood and difficult to model. A larger-than-expected radius for a given mass would suggest a “softer” EOS, implying easier compressibility, while a smaller radius indicates a “stiffer” equation. The ongoing refinement of both theoretical predictions and observational constraints – particularly through gravitational wave astronomy and precise radius measurements – is crucial to narrowing the possibilities and ultimately revealing the true nature of matter in these cosmic laboratories.

Building Models of the Void: Reaching for First Principles
Many-body calculations are foundational for deriving an Equation of State (EOS) from first principles due to their ability to account for the complex interactions between nucleons within a nucleus or nuclear matter. Techniques like Brueckner-Hartree-Fock (BHF) solve the many-body Schrödinger equation approximately by considering two-body interactions and incorporating three-body forces for increased accuracy. Relativistic extensions of BHF, such as the Relativistic Brueckner-Hartree-Fock (RBHF) approach, become essential at higher densities and energies where relativistic effects are significant. These calculations determine the energy per nucleon as a function of density, providing the thermodynamic data needed to construct the EOS and model the behavior of matter in extreme conditions, such as those found in neutron stars or heavy-ion collisions. The computational complexity of these methods necessitates approximations and careful validation against experimental data.
Energy-density-functional (EDF) theories represent an alternative to many-body calculations for determining the equation of state (EOS). These approaches, such as Skyrme-Hartree-Fock, Gogny-Hartree-Fock, and Relativistic Mean-Field, calculate the energy of nuclear matter as a functional of the density, thereby circumventing the explicit treatment of complex many-body interactions. Skyrme and Gogny EDFs are non-relativistic, employing effective nucleon-nucleon interactions, while Relativistic Mean-Field theory incorporates special relativity and utilizes effective meson exchange interactions. Each method differs in its treatment of density dependence and incorporates various parameters fitted to experimental data, such as saturation properties of nuclear matter and finite nuclei characteristics. The computational efficiency of EDF methods allows for calculations across a wider range of conditions compared to more computationally intensive many-body approaches.
Meta-modeling of the Equation of State (EOS) utilizes simplified parameterizations to efficiently explore a wide range of potential EOS behaviors. A common technique involves employing a Constant Speed of Sound (c_s) parameterization, where c_s is treated as a constant or a function of density. This approach reduces the computational cost associated with full many-body or energy-density-functional calculations, allowing for rapid evaluation of the EOS across a broad parameter space relevant to astrophysical simulations and experiments. By varying parameters within the meta-model, researchers can assess the impact of different EOS assumptions on observables without undertaking computationally expensive ab initio calculations for each scenario.

The Weight of Observation: Constraining the Unknown
Accurate determinations of neutron star masses, exemplified by PSR J0740+6620 (2.1 ± 0.03 M_{\odot}) and PSR J0030+0451 (1.46 ± 0.03 M_{\odot}), directly constrain the equation of state (EOS) by establishing minimum mass limits for given radii. These measurements, derived from Shapiro delay and gravitational redshift observations, effectively rule out EOS models that predict lower maximum masses or larger radii for a given mass. Specifically, the observed masses necessitate a stiff enough EOS to prevent gravitational collapse, and the precision of these measurements allows for increasingly stringent tests of theoretical EOS predictions across a range of densities.
Observations of millisecond pulsars and low-mass X-ray binaries exhibiting rapid rotation and relatively low masses – exemplified by XTE J1814-338 (0.70 solar masses) and PSR J0614-3329 (1.44 solar masses) – necessitate EOS constructions that support stable neutron stars at masses previously considered unattainable. Conventional EOS models, particularly those relying on purely nucleonic compositions or relatively “soft” symmetry energies, often predict maximum neutron star masses too low to accommodate these observed objects. The existence of these smaller, rapidly rotating neutron stars implies a sufficiently high pressure support at subnuclear densities, requiring either exotic matter compositions – such as hyperons, quarks, or condensates – or a stiffer EOS than previously assumed. These observations therefore place stringent lower limits on the pressure-density relationship at high densities, guiding the refinement of theoretical models and the exploration of alternative compositions.
Bayesian inference statistically combines data from neutron star observations with predictions from theoretical equations of state (EOS) to constrain the possible EOS parameters. Analyses utilizing data from pulsars HESS J1731-347 and XTE J1814-338, when subjected to Bayesian inference, indicate a transition density – representing a change in the internal composition of neutron stars – of approximately 2.75−0.33+0.47 times the saturation density of nuclear matter. This value is derived from posterior distributions generated by fitting theoretical models to observed mass and radius constraints, and reflects the uncertainty inherent in both the observational data and the theoretical models employed.

Beyond the Nucleon: The Ghosts of Matter’s Past
Within the extreme densities of neutron stars, matter isn’t limited to the familiar protons and neutrons; it can undergo a phase transition to exotic states, potentially including quark-gluon plasma. This isn’t simply a change in temperature or pressure, but a fundamental restructuring of matter itself, where quarks and gluons-typically confined within hadrons-become deconfined and move freely. This transition, driven by immense gravitational pressure, is predicted by theoretical models and could dramatically alter the star’s properties. The resulting quark-gluon plasma represents one of the earliest states of matter in the universe, briefly existing after the Big Bang, and studying its potential presence within neutron stars offers a unique window into the conditions of that primordial era. Furthermore, the formation of such exotic matter impacts the star’s equation of state, influencing its mass, radius, and stability, and creating possibilities for observable signatures distinct from those of traditional neutron stars.
The universe may harbor “Twin Stars”-compact objects resembling neutron stars, yet fundamentally distinct in their internal composition. Current understanding posits neutron stars as primarily composed of neutrons, but extreme densities within these objects could trigger phase transitions, potentially birthing stars made of alternative matter like quark-gluon plasma. A Twin Star, therefore, represents a celestial body with a similar mass and radius to a standard neutron star, but a drastically different equation of state governing its core. Distinguishing these elusive objects requires precise measurements of their properties, particularly mass and radius, and scrutinizing subtle deviations from the expected behavior of ordinary neutron stars. The existence of Twin Stars would not only confirm the exotic nature of ultra-dense matter but also necessitate a reevaluation of current models describing the endpoint of stellar evolution.
Investigating the potential for exotic matter within neutron stars necessitates sophisticated modeling, primarily through the application of the Tolman-Oppenheimer-Volkoff (TOV) equation – a cornerstone of relativistic stellar structure. Recent analysis, incorporating constraints derived from observed tidal deformability, reveals compelling evidence for a dramatic shift in energy density at the star’s core, quantifying a jump of approximately 720.73 MeV. Crucially, the calculated squared sound speed of 0.85 following this transition indicates a remarkably ‘stiff’ equation of state. This stiffness suggests that the post-transition matter – potentially quark-gluon plasma or other exotic phases – resists compression more strongly than standard neutron star material, influencing the star’s overall structure and stability and potentially leading to the formation of objects like twin stars with distinct compositions.
The pursuit of a definitive equation of state for dense matter, as explored in this study of compact stars, feels remarkably akin to peering into the abyss. This research, with its focus on twin star scenarios and phase transitions, attempts to constrain the possibilities, yet the universe often delights in exceeding even the most sophisticated models. As Stephen Hawking once observed, “It is not enough to be right. One must also be able to explain why one is right.” The elegance of a theory, particularly concerning something as fundamentally unknowable as the interior of a neutron star, all looks pretty on paper until you look through a telescope – or, in this case, analyze gravitational waves and radius measurements. The very notion of seeking a ‘true’ equation of state is, perhaps, a testament to humanity’s enduring, and possibly delusional, belief in ultimate understanding.
The Horizon Beckons
The search for the equation of state of dense matter, as exemplified by this work, reveals less about the stars themselves and more about the limits of prediction. A strong first-order phase transition to quark matter, even if suggested by observations of small-radius compact stars, remains a theoretical construct. Any inferred parameter, any carefully calculated equation, is ultimately subject to the crushing weight of gravity – a probability, not a certainty. The ‘twin star’ scenario offers a compelling, testable hypothesis, yet the universe does not offer guarantees; it offers only data, which are inevitably incomplete.
Future work will undoubtedly refine Bayesian techniques and expand observational datasets. However, the fundamental problem persists: the event horizon of a neutron star is also the horizon of knowledge. Each additional decimal point of precision in a mass-radius measurement merely delays the inevitable confrontation with the unknown. It is a temporary reprieve, not a solution. The true equation of state, if it exists as a fixed entity, may be intrinsically unknowable.
The next step isn’t simply more data, or better algorithms. It’s a quiet acceptance that any model, however elegant, is provisional. Black holes don’t argue; they consume. And in the end, all theories, like all matter, are subject to that same fate.
Original article: https://arxiv.org/pdf/2604.21379.pdf
Contact the author: https://www.linkedin.com/in/avetisyan/
See also:
- Adam Levine Looks So Different After Shaving His Beard Off
- After AI Controversy, Major Crunchyroll Anime Unveils Exciting Update
- Gold Rate Forecast
- Tekken 8 Fans Furious as Tifa Tipped for Street Fighter 6 Instead
- Dialoop coming to Switch on June 17
- 10 Movies That Were Banned in Different Countries For Random Reasons
- From season 4 release schedule: When is episode 2 out on MGM+?
- Steam Makes Sci-Fi Game 100% Free for 72 Hours
- Japan’s No. 1 Spring 2026 Anime Is the True Successor to an All-Time Great
- Welcome to Demon School! Iruma-kun season 4 release schedule: When are new episodes on Crunchyroll?
2026-04-25 08:58