⚛️ nuclear theory

Stability of hybrid stars via catastrophe theory

This paper applies catastrophe theory to analyze the mass-radius topology of hybrid stars, demonstrating how a single first-order phase transition between nuclear and quark matter governs the emergence and stability of different compact star branches through bifurcation structures.

Original authors: Eduardo S. Fraga, Sergio E. Jorás

Published 2026-09-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Eduardo S. Fraga, Sergio E. Jorás

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 within the heart of a dead star, where gravity crushes matter to densities that defy everyday experience, the laws of physics undergo a dramatic transformation. These stellar remnants, known as neutron stars, are so dense that a single teaspoon of their material would weigh a billion tons on Earth. For decades, physicists have wondered what happens when the pressure inside such a star becomes so extreme that the protons and neutrons themselves begin to break apart, dissolving into a soup of their constituent parts called quarks. This transition from ordinary nuclear matter to quark matter creates a "hybrid star," a celestial object with a core of exotic quark fluid surrounded by a shell of traditional nuclear material. Understanding whether these hybrid stars can exist stably, or if they would collapse under their own weight, is crucial for mapping the limits of matter in the universe. To solve this, scientists must analyze the delicate balance between the crushing force of gravity and the internal pressure of the star, a relationship usually described by complex equations that are difficult to solve when a sudden phase change occurs inside the star.

In a recent study, researchers Eduardo S. Fraga and Sergio E. Jorás from the Federal University of Rio de Janeiro approached this problem not by calculating every specific detail of the star's interior, but by looking at the broader shape of the solutions themselves. They applied a branch of mathematics called catastrophe theory, which studies how small changes in a system's conditions can lead to sudden, dramatic shifts in its behavior. Imagine a landscape where the height represents the stability of a star; as you change the conditions, the hills and valleys of this landscape shift. Sometimes, a valley where a stable star can sit simply vanishes, or two hills merge and disappear. The authors used this mathematical framework to map out the possible shapes of the relationship between a hybrid star's mass and its radius. They found that the existence of stable hybrid stars depends entirely on the structure of these mathematical singularities, or "fold points," where different possible states of the star meet and merge.

The researchers began by simplifying the problem to its most basic form, imagining a star made of two distinct, unchangeable layers of fluid, much like a ball with a dense core and a lighter shell. Even in this simplified model, they discovered that the transition between the two layers acts as a critical point. As they adjusted the pressure required to trigger the change from nuclear to quark matter, and the difference in density between the two states, the number of possible stable stars changed. They observed that for certain combinations of these conditions, a stable branch of stars would appear or disappear entirely. This happens because the mathematical curve that describes the star's mass and radius develops a sharp turn or a "fold." When the conditions cross a specific threshold, the stable solution and an unstable solution crash into each other and vanish, leaving no room for that type of star to exist. This phenomenon explains why some theoretical models predict the existence of "twin stars"—two different types of stars with the same mass but different sizes—while others do not.

Moving beyond the simplified model, the authors applied these insights to more realistic scenarios where the density of the star changes gradually rather than in sharp steps. They found that the complexity of the system increases, requiring a more intricate mathematical shape to describe the possible outcomes. In this more detailed view, the landscape of possibilities becomes a "swallowtail" shape, a specific geometric form known in mathematics that allows for four different equilibrium states to exist simultaneously in certain regions. Within this complex structure, there are zones where two stable branches of stars can coexist, separated by a region of instability. The researchers showed that the specific properties of the phase transition—how much the density jumps and at what pressure it occurs—determine which of these zones the star occupies. If the transition is too sharp or the density jump is too large, the stable branch for hybrid stars may never form, or it may merge with an unstable branch and disappear, preventing the star from existing in a stable state.

The study concludes that the stability of these exotic stars is not just a matter of the specific materials inside them, but is fundamentally dictated by the topology of the mathematical functions that describe them. By treating the star's stability as a landscape that can fold and merge, the authors provided a new way to understand why some hybrid stars are possible while others are not. They demonstrated that the appearance or disappearance of entire families of stars is a direct consequence of how the underlying mathematical potential changes as the conditions of the phase transition vary. While the paper does not yet provide a definitive list of which specific equations of state will produce these stars, it establishes a clear rule: the existence of stable hybrid stars is governed by the geometric structure of the solutions, specifically the points where stable and unstable states coalesce. This approach offers a powerful tool for future research, allowing scientists to predict the stability of these cosmic objects by analyzing the shape of the mathematical curves rather than getting lost in the details of every possible equation of state.

Technical Summary: Stability of Hybrid Stars via Catastrophe Theory

Problem Statement
The stability of hybrid stars—compact objects containing a core of deconfined quark matter surrounded by a hadronic mantle—depends critically on the Equation of State (EoS) derived from Quantum Chromodynamics (QCD). However, a nonperturbative first-principles solution for QCD at finite baryon density remains elusive due to the Sign Problem in lattice QCD. While previous studies have investigated how phase transition properties (such as transition pressure $p_{tr}$ and energy density discontinuity $\Delta\epsilon$ ) influence the mass-radius relation, the structural topology of these solutions and the mechanisms governing the emergence or disappearance of stable branches have not been fully analyzed through the lens of singularity theory. The authors posit that the stability of hybrid stars may be less sensitive to specific EoS details and more dependent on the structural properties of the underlying potential, specifically its singularity profile in the case of a sharp first-order phase transition.

Methodology
The authors employ Catastrophe Theory to analyze the topologies of the mass-radius relation for hybrid stars. The methodology proceeds through the following steps:

Theoretical Framework: The study utilizes the Tolman-Oppenheimer-Volkov (TOV) equations to describe hydrostatic equilibrium in General Relativity. For hybrid stars, a sharp first-order phase transition is assumed between nuclear and quark matter, characterized by a pressure $p_{tr}$ and an energy gap $\Delta\epsilon$ . This avoids the need for a Gibbs construction.
Analytical Warm-up: To illustrate the basic mechanisms, the authors first analyze "bi-homogeneous hybrid stars" composed of two incompressible fluids. This allows for analytical solutions where the mass $M$ is expressed as a function of dimensionless control parameters: $x \equiv p_{tr}/\epsilon_{out}$ and $y \equiv \Delta\epsilon/\epsilon_{out}$ .
Catastrophe Theory Application: The authors map the stellar equilibrium configurations to the extrema of an effective potential function $V(z)$ $V (z)$ , where $z$ $z$ represents the central energy density (or pressure).
- For standard relativistic stars, the mass-density relation near the maximum mass is approximated by a fold catastrophe (cubic potential), where stable and unstable branches coalesce.
- For hybrid stars with a first-order phase transition, the authors identify a richer structure. By analyzing the mass-central density diagram (referencing prior work [13]), they determine that the system requires a potential with four extrema (two stable, two unstable) to describe the full set of solutions. This corresponds to a swallowtail catastrophe (quintic potential) with three control parameters.
Bifurcation Analysis: The study investigates the bifurcation set—the locus of points in the parameter space where extrema coalesce. The authors track how the number and nature of equilibrium solutions (stable vs. unstable branches) change as the control parameters ( $p_{tr}$ and $\Delta\epsilon$ ) vary.

Key Contributions and Results

Topological Classification: The paper demonstrates that the topology of the mass-radius diagram for hybrid stars is defined by the singularity structure of an underlying catastrophe potential. Different branches of compact stars emerge or disappear based on the bifurcation structure of this potential.
Seidov Criterion and Instability: The authors revisit the Seidov criterion, showing how a first-order phase transition can drive instability. In the incompressible fluid model, they derive the condition under which the derivative of mass with respect to central pressure changes sign, signaling the onset of instability.
The Swallowtail Structure: For realistic hybrid stars, the authors propose that the equilibrium surface is governed by a swallowtail catastrophe.
- Region A: In the central region of the parameter space, the potential possesses four extrema: two local maxima (unstable branches) and two local minima (stable branches). This corresponds to the existence of "twin star" solutions and disconnected branches.
- Coalescence: As control parameters change, pairs of neighboring extrema coalesce.
  - Crossing into regions $C_1$ or $C_2$ results in the disappearance of one pair of extrema (leaving one stable and one unstable branch).
  - Crossing into region $B$ results in the coalescence of all extrema, leaving no hybrid star solutions (only the hadronic branch remains).
Parameter Dependence: The analysis shows that the relative positions of the threshold points (where $dM/d\epsilon_c = 0$ ) depend on the dimensionless parameters $x$ and $y$ . The coalescence of these points occurs only when $p_{tr} \to 0$ or $\Delta\epsilon \to 0$ , corresponding to homogeneous stars.

Significance and Claims
The paper claims to offer a different perspective on hybrid star stability that has been overlooked in previous analyses. By shifting the focus from specific EoS parametrizations to the structural properties of singularities, the authors argue that the emergence of complex mass-radius topologies (such as disconnected "twin star" branches) is a generic feature of the bifurcation structure of the underlying catastrophe potential.

The authors explicitly state that their work provides a qualitative discussion explaining the (dis)appearance of stable and unstable branches as control parameters change. They do not claim to provide quantitative predictions for specific EoS models in this study. Instead, they establish a framework where the stability structure is defined by the normal form of the swallowtail catastrophe. The paper concludes by noting that future work will apply this catastrophe theory framework to realistic equations of state to provide quantitative predictions regarding the different sectors of the bifurcation set and to determine if regions $C_1$ and $C_2$ yield physically distinct stellar characteristics.

More like this