Why Stellar Sequences Turn Over: Fixed Points, Instability, and Equation-of-State Universality

arXiv:2603.26973 · astro-ph.HE, gr-qc, nucl-th · Submitted 2026-03-27 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Today's paper: "Why Stellar Sequences Turn Over".

Jocelyn: The gist: The maximum mass of stellar sequences arises from the existence of a fixed point in the relativistic regime,

Vera: First, who's behind it and why it matters.

Title and authors: Vera: So we're looking at this paper today, "Why Stellar Sequences Turn Over: Fixed Points, Instability, and Equation-of-State Universality." It’s by Legred and Yunes. Basically, they’re taking the complex equations for how stars are structured and treating them like a dynamical system.

Jocelyn: That sounds intense. So what's the main idea here for someone just tuning in? Are we talking about finding a neat mathematical pattern that explains why the heaviest stable stars have a limit on how big they can get?

Subrahmanyan: Exactly. The core concept is that the maximum mass of stellar sequences comes from finding a fixed point in the relativistic regime. This fixed point is directly linked to where you see that turnover in the mass-radius curve.

Vera: Right, so instead of just solving for one star at a time, they’re looking at an entire family of stars as trajectories in some kind of phase space. This approach lets them find these organizing principles that aren't just dependent on the specific equation of state details.

Jocelyn: So they’re trying to see if this fixed point idea applies across different types of stellar matter, whether we’re talking about neutron stars or something else.

Subrahmanyan: They start by using the enthalpy-based formulation of the Tolman–Oppenheimer–Volkoff equations and rewriting that system as a two-dimensional, non-autonomous flow in a compact set of variables. The state vector they use is energy density and compactness, defined as x(τ) = (˜e(τ), v(τ)) where tau is the logarithm of enthalpy.

Vera: That sounds like a lot to track. So what happens when you analyze this flow? Does it actually lead to something concrete about the star's structure?

Subrahmanyan: They identify a fixed point in that energy density and compactness plane where both are dimensionless versions of the energy density and the enclosed mass. This location is determined by just a few thermodynamic combinations, specifically the ratio of pressure to energy density, or w, and the speed of sound squared.

Jocelyn: That means instead of needing every single microphysical detail for every star, you only need those ratios to pinpoint where the structure settles. That’s a big simplification for observation.

Vera: And they find that high-mass relativistic stars generically evolve near this fixed point over the part where most of their mass is built up. This evolution creates that spiral and turnover behavior we see in the mass-radius curve, and it explains why the termination of these sequences is so insensitive to EoS details once you force them to be physical.

Title and authors: Subrahmanyan: The real parts of the eigenvalues of the Jacobian matrix are complex conjugate pairs for physical equations of state with a speed of sound squared greater than zero. This means that while the frozen-time equilibrium is an unstable focus when you evolve forward in the logarithm of enthalpy, it acts as an attractor along the outward integration direction.

Jocelyn: An unstable focus that still pulls things towards a limit—that’s interesting. But what about when we move to less extreme scenarios, like Newtonian physics?

Vera: That’s where they test if this same structure exists in Newtonian equations for weakly gravitating stars. They find that periodic or spiraling behavior can appear there, but only under very special conditions where the equation of state approaches a constant sound speed.

Subrahmanyan: In practice, this mechanism is not generic in Newtonian physics because it requires tuning the EoS so that the speed of sound squared stays approximately constant as you go through the logarithm of enthalpy. This leads them to identify what they call the "compressible limit," where a common internal structure emerges that doesn't depend much on microphysical details.

Jocelyn: So, if we look at Newtonian physics, it seems like the universe settles into a standard shape regardless of minor changes in the stuff inside, provided we are in that compressible limit.

Vera: But they do derive a separate quasi-universal relation for the maximum mass in terms of central quantities specifically within that limit. It’s a different structure than what you get from general relativistic physics.

Subrahmanyan: The fixed-point geometry gives them compact scaling relations for the maximum mass and radius based on two things: first, a high-density stiffness scale, like c2s,typ or wtyp, which sets up the fixed point; and second, an exit density e0 that marks where the star leaves that near-fixed-point region.

Jocelyn: So if we want to estimate a maximum mass for a star without knowing its exact equation of state all over its volume, we only need those two characteristic scales.

Vera: The leading-order estimate they get for the maximum mass is given by Mmax ≃ c4 G v to the three over two square root of four pi G e0. This scaling relation holds true even when tested against a Gaussian ensemble of equations of state, showing it’s insensitive to order ten percent variation.

Subrahmanyan: Astrophysical constraints actually require a lower bound on this maximum mass, which suggests that the sound speed in the dense matter has to increase quite rapidly between about two and four times nuclear saturation density.

Title and authors: Jocelyn: That means the physics at those extreme densities really matters for setting those limits, not just some small correction factor. So what’s next? How do we see this in action?

Vera: The paper distinguishes between two ways things can terminate: one controlled by the relativistic fixed-point geometry and another controlled by near-critical structure in the weak-field regime.

Subrahmanyan: The ultrarelativistic comparison curve, which they derive from rewriting that fixed-point estimate, is different from the Newtonian compressible limit relation because they come from different limiting mechanisms. This distinction is important for understanding where these universal relations apply.

Jocelyn: So we have two potential limits to consider when thinking about stellar stability: one relativistic and one Newtonian, and they behave differently at their extremes.

Vera: The paper shows that for generic neutron stars, the dynamical systems picture isn't super insightful because the trajectory in that energy density space doesn't always perfectly follow the fixed point. But near the maximum mass, dynamically stable configurations do show a strong reduction to just that two-dimensional fixed point location.

Subrahmanyan: The conclusion is that it’s those compact scaling relations for Mmax and Rmax based on the high-density stiffness scale and exit density e0 that are the main thing here. The dimensioness correction factor alpha encodes the finite spiral evolution before that star leaves the fixed-point regime.

Jocelyn: So, to wrap up, this paper gives us a way to predict maximum mass based on core stiffness and an exit density for stars in both relativistic and Newtonian contexts. It connects microphysics to macroscopic limits.

Vera: Precisely. The main implication is that the maximum mass isn't just some lucky accident of the EoS; it’s controlled by that fixed point geometry and those few scales that determine when a star exits the high-density region.

Subrahmanyan: I think what this really highlights is how core stiffness, specifically related to c2s,typ and wtyp, primarily dictates the maximum mass scaling, schematically scaling like c2s to the three over two.

Jocelyn: So we’re looking at how that core stiffness sets the stage for these universal mass limits across different gravitational regimes. That gives us a concrete tool for testing different nuclear models against observational constraints.

Vera: It does, and it helps us focus our EoS searches on the right regions of high-density physics. We'll see what this means for how we interpret future pulsar data next time.

The paper's summary: Vera: So we’ve looked at how these stars are structured as trajectories in some phase space, and now we’re going to break down what they found about those fixed points and why they matter for stellar sequences.

Jocelyn: Basically, the main thing is that the maximum mass of a star sequence isn't just random; it comes from this specific mathematical point where things stop changing in a certain way.

Subrahmanyan: Exactly, Jocelyn. They show that this fixed point in the relativistic regime is tied directly to where you see that turnover happening on the mass-radius curve. It’s the organizing principle for why things cap out at a certain mass before they fall off again.

Vera: And they use this framework to explain how those sequences behave, suggesting that if you look at a whole family of stars, their evolution follows these predictable paths around this fixed point. It moves the focus away from just solving one star model at a time.

Jocelyn: That’s the big shift for us observing pulsars too; it means we don't have to test every tiny equation of state variation individually to predict where the mass limit lies.

Subrahmanyan: Right, and they find that this mechanism is really robust across different physics, even when you switch from general relativity over neutron stars to simpler Newtonian gravity for less extreme cases.

Vera: It’s interesting how they found that in the Newtonian limit, you can still see some spiral behavior but only under very specific conditions where the equation of state has a constant sound speed. It’s not as straightforward as the relativistic version.

Jocelyn: So, it means for our surveys looking at neutron stars, we have to be careful about which limit we're actually in when we try to estimate those maximum masses based on what we see.

Subrahmanyan: That’s a fair point, Jocelyn. The paper doesn't just give us one answer; it shows that the relativistic fixed point and the Newtonian compressible limit are two different ways termination can happen, so you have to know which mechanism is at play.

Vera: And they give us some very concrete scaling relations for that maximum mass based on core stiffness and an exit density, like e0. It’s not just some abstract theory; it gives us a formula we can actually use to put constraints on the nuclear physics inside those stars.

Jocelyn: So, if I'm listening—whether I'm looking at data or just following the sky—the point is that the maximum mass scales primarily with how stiff the matter gets in the core.

Subrahmanyan: That’s what they emphasize. The leading-order estimate for maximum mass depends mainly on a characteristic speed of sound in that dense region, which scales like c2s to the three over two.

Vera: And they showed this scaling relation holds up even when you test it against a wide range of possible equations of state, which is pretty significant for testing our nuclear models.

Jocelyn: It sounds like they’ve given us a better way to look at the constraints we're putting on these exotic objects, moving from just checking individual stars to looking at the underlying mathematical structure governing all of them.

Subrahmanyan: Precisely. And because this fixed point analysis also helps us diagnose what happens if there’s a strong phase transition in that matter, it opens up a new way for us to test those extreme physics scenarios.

The paper's improvements: Tom: So we've talked about how the core physics dictates these maximum masses, and now we're going to look at how those researchers think they can make this whole dynamical systems approach even better.

Vera: Basically, the paper suggests a few ways to refine the analysis of those fixed points and their stability. They’re not just stating facts; they’re suggesting improvements to the mathematical setup.

Jocelyn: Right, so what are these suggested improvements? Are we talking about better numbers or a different way of looking at the flow?

Subrahmanyan: They're focusing on making that two-dimensional flow even more explicit, using tools like the Jacobian matrix to classify solutions more thoroughly. This helps them pinpoint exactly where things become unstable or stable in that energy density space.

Vera: That makes sense because they found earlier that the fixed point might not always be perfectly described by a simple trajectory, so mapping out those local behaviors is crucial for understanding why some stars stop and others keep going.

Jocelyn: So, instead of just saying "it's unstable," they're showing us *how* it becomes unstable based on the specific parameters of the equation of state. That gives us a better handle on the physics driving that turnover.

Subrahmanyan: Yes, and they’re exploring how a strong first-order phase transition in the matter could actually shift where that fixed point is located, which would substantially change how we predict the maximum mass limit.

Vera: That’s interesting because if a phase transition happens, it introduces a new physics element that could completely alter the predicted outcome of the spiral evolution before it terminates. It gives us a diagnostic tool for what kind of microphysics we should expect to see in those extreme objects.

Jocelyn: So, if we were trying to use this paper to test different nuclear models, these suggestions give us a clearer signal for when a specific type of phase transition might be the reason for an observed mass limit.

Subrahmanyan: Right. And they're also looking at how these results translate into actual astrophysical constraints, applying that scaling relation to current pulsar data to see what it tells us about the core stiffness parameters like c2s and wtyp.

Vera: So, the paper is essentially saying we need better mathematical tools—like a more detailed look at the flow dynamics and phase transition effects—to fully harness this idea of equation-of-state universality.

Jocelyn: It sounds like they’re moving beyond just finding a scaling law and trying to build a more comprehensive picture of the physical mechanisms that cause stellar sequences to terminate.

Subrahmanyan: That’s the direction they're pushing in, Jocelyn. They want to connect these abstract dynamics back into observable limits on the structure of neutron stars themselves.

Conclusion: Tom: So we’re wrapping up this look at "Why Stellar Sequences Turn Over: Fixed Points, Instability, and Equation-of-State Universality." We’ve covered how these fixed points control the maximum mass of stars across different gravity regimes.

Vera: Exactly. The main takeaway is that we can use this dynamical systems view to get compact scaling relations for maximum mass based on core stiffness and an exit density, e0. It’s a pretty clean way to relate microscopic physics to macroscopic limits for neutron stars.

Jocelyn: I think what this really means for us observing pulsars is that when we try to constrain the equation of state, we have a specific mathematical target—that fixed point geometry—to aim for in our models.

Subrahmanyan: That’s right. And the paper shows that because these scaling relations are so robust, they give us a much more reliable way to test different nuclear models against what we actually see in the sky.

Vera: It really boils down to this: whether you're looking at relativistic or Newtonian physics, there's a universal structure governing how those stars turn over at their maximum mass.

Jocelyn: And it changes how we frame the questions we ask about dense matter; instead of just asking about one star’s structure, we start asking about the phase space trajectory of the whole sequence.

Subrahmanyan: It helps us see that core stiffness, which is described by parameters like c2s and wtyp, is the primary driver for setting those mass limits in this framework.

Vera: So to summarize, "Why Stellar Sequences Turn Over: Fixed Points, Instability, and Equation-of-State Universality" gives us a consistent mathematical language to predict maximum stellar masses using only a few core stiffness numbers.

Jocelyn: It’s a useful tool for our surveys because it helps us define the boundaries of what we should expect to see in terms of mass limits.

Subrahmanyan: And the future work they point toward is using this framework to diagnose specific microphysical behaviors, like strong phase transitions, that would alter those predictions.

Illinois Center for Advanced Studies of the Universe · Department of Physics, University of Illinois at Urbana-Champaign

astro-ph.HE, gr-qc, nucl-th

Submitted: 2026-03-27

Updated: 2026-10-09

Comments: 30 Pages, 13 Figures

Journal ref: Phys. Rev. D 114, 064006 Published 2 September, 2026

DOI: 10.1103/y75z-bbyk

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 83/100

The gist: The gist: The maximum mass of stellar sequences arises from the existence of a fixed point in the relativistic regime, which is directly associated with the turnover of the mass–radius curve.

Key concepts

Fixed Point
In dynamical systems terms, a fixed point is a stable state where the system stops changing over time. For these stars, this occurs in the relativistic regime when the energy density and mass compactness reach specific dimensionless values. This fixed point dictates the maximum mass achievable before instability sets in, directly causing the stellar sequence to turn over.
Compressible Limit
This limit describes a universal internal structure that emerges when stellar sequences are analyzed under Newtonian gravity. It occurs because, under these conditions, the behavior becomes largely independent of the specific microphysical details of the equation-of-state. This provides a common framework for understanding how stars behave near their maximum mass.
Relativistic Fixed Point Mechanism
This mechanism uses a two-dimensional flow in phase space to model stellar evolution. The fixed point is where the system exhibits spiraling behavior, which manifests as the turnover in the mass-radius curve. This geometric structure explains why high-mass stars evolve toward this point, linking relativistic dynamics to observable stellar properties.

Terminology

Summary

The gist: The maximum mass of stellar sequences arises from the existence of a fixed point in the relativistic regime, which is directly associated with the turnover of the mass–radius curve.

Theoretical Framework

The authors reformulate stellar structure equations in dynamical systems and show that this approach reveals organizing principles that transcend EoS details and provide a unified explanation for quasi-universal behavior <ref:2603.26973#pg2> The basic idea is to stop thinking of a stellar model as a single numerical solution and instead view an entire sequence of stars as a family of trajectories in a low-dimensional phase space <ref:2603.26973#pg2> They begin by working in the enthalpy-based formulation of the relativistic TOV equations and rewriting the system as a two-dimensional, non-autonomous flow in a compact set of variables <ref:2603.26973#pg2> The state vector is defined as x(τ) = (˜e(τ), v(τ)) where τ = ln h, ˜e is the dimensionless energy-density variable, and v is the compactness of the enclosed mass <ref:2603.26973#pg2>

Relativistic Fixed Point Mechanism

The analysis identifies a fixed point of the frozen system in the (˜e, v) plane where ˜e ∝ r squared e and v = Gm(r)/(rc 2) are dimensionless versions of the energy density and the enclosed mass <ref:2603.26973#pg2> The location of this fixed point is set by a small set of thermodynamic combinations, essentially the ratio of pressure to energy density, w = p/e, and the speed of sound squared, c 2s = dp/de <ref:2603.26973#pg2> High-mass relativistic stars generically evolve near this fixed point over the region where the bulk of the mass is accumulated, producing spiraling and turnover behavior in the M–R curve <ref:2603.26973#pg2> The real parts of the eigenvalues of the Jacobian J⋆ are complex conjugate pairs for physical EoSs with c 2s > 0, meaning the frozen-time equilibrium is an unstable focus when evolved forward in ln h, but acts as an attractor along the outward integration direction <ref:2603.26973#pg2>

Newtonian and Compressible Limits

The authors ask whether the same organizing structure is present in Newtonian equations, showing that periodic or spiraling behavior can appear only under special conditions where the EoS approaches a constant sound speed <ref:2603.26973#pg2> In practice, this mechanism is non-generic in Newtonian physics because it requires tuning the EoS so that c 2s is approximately constant in ln h <ref:2603.26973#pg2> This leads to the identification of the compressible limit, where a common internal structure emerges that is largely independent of microphysical details <ref:2603.26973#pg2> In this limit, they derive a separate quasi-universal relation for the maximum mass in terms of central quantities <ref:2603.26973#pg2>

Maximum Mass Scaling and Constraints

The fixed-point geometry yields compact scaling relations for Mmax and Rmax in terms of (i) a high-density stiffness scale (c 2s,typ, wtyp) that determines the fixed point, and (ii) an exit density e0 that marks where the star leaves the near-fixed-point regime <ref:2603.26973#pg2> The leading-order estimate for Mmax is given by Mmax ≃ c 4 G v(3/2) sqrt(4πGe0) <ref:2603.26973#pg2> This scaling relation is shown to be EoS insensitive to order 10% variation when evaluated against a Gaussian-process ensemble of EoSs <ref:2603.26973#pg2> Astrophysical constraints imply a lower bound on Mmax, suggesting that the dense-matter sound speed must rise rapidly between ∼ 2–4 ρnuc <ref:2603.26973#pg2>

Comparison of Limits and Implications

The analysis distinguishes between two mechanisms for termination: one controlled by the relativistic fixed-point geometry and one controlled by near-critical structure in the weak-field regime <ref:2603.26973#pg2> The ultrarelativistic comparison curve, derived from rewriting the fixed-point estimate, differs from the Newtonian compressible-limit relation because they arise from different limiting mechanisms <ref:2603.26973#pg2> The results show that for generic neutron stars, the dynamical systems picture is not particularly insightful because the orbit of the trajectory in (˜e, v) space is never truly well-described by the fixed-point <ref:2603.26973#pg2> However, near the maximum mass, dynamically stable configurations show strong dimensionality reduction to just the location of the (2D) fixed point <ref:2603.26973#pg2>

Conclusion

The fixed-point geometry yields compact scaling relations for Mmax and Rmax in terms of (i) a high-density stiffness scale (c 2s,typ, wtyp) that determines the fixed point, and (ii) an exit density e0 that marks where the star leaves the near-fixed-point regime <ref:2603.26973#pg2> The dimensioness correction factor α encodes the finite spiral evolution before that exit <ref:2603.26973#pg2> The maximum mass depends primarily on a characteristic sound speed in the core, schematically scaling like ⟨c 2s⟩(3/2) <ref:2603.26973#pg2>

How it works

The authors recast the TOV system as a two-dimensional flow whose explicit EoS dependence enters through only a small set of thermodynamic combinations <ref:2603.26973#pg2> They define the state vector and time variable as x(τ) = (˜e(τ), v(τ)) where τ = ln h <ref:2603.26973#pg2> The TOV equations become dx/dτ = F(x, τ) <ref:2603.26973#pg2>

Fixed Points and Linearization

The analysis identifies a fixed point of the frozen system in the (˜e, v) plane where ˜e ∝ r squared e and v = Gm(r)/(rc 2) are dimensionless versions of the energy density and the enclosed mass <ref:2603.26973#pg2> Near this fixed point, the local behavior is governed by the eigenvalues of the Jacobian J0, which classify solutions as a focus if det J0 > 0 and (tr J0) squared This focus structure is identified as the geometric engine behind the turnover of the M–R curve and the appearance of spirals at high central density <ref:2603.26973#pg2>

Spiral Evolution toward Fixed Point

The linearized system around x⋆ yields eigenvalues λ±(τfr) = σ(τfr) ± iω(τfr), where σ(τfr) is the real part and ω(τfr) is the imaginary part <ref:2603.26973#pg2> The frozen-time equilibrium is an unstable focus when evolved forward in ln h, but it acts as an attractor along the physical (outward) integration direction s:= ln hc - ln h <ref:2603.26973#pg2> The amplitude of the spiral can be characterized by a multiplier µ:= exp(-γT) without solving the ODE in closed form <ref:2603.26973#pg2>

Compressible Limit and Universal Core Behavior

In the Newtonian limit, sequences terminate because equilibrium configurations become strongly centrally condensed, leading to the compressible limit where a common internal structure is largely independent of microphysical details <ref:2603.26973#pg2> In this limit, near-core expansion yields u¯(s) = 6s + 9/5 c 2s,c s squared + O(s 3) <ref:2603.26973#pg2> This structure is crucial for understanding the quasi-universal scaling of Mbarmax in Eq. (78) <ref:2603.26973#pg2>

Post-Newtonian Extension and Astrophysical Constraints

The post-Newtonian extension of the compressible limit provides a simple quasi-universal relation for neutron stars that connects the maximum mass to the central stiffness of the maximum-mass configuration <ref:2603.26973#pg2> The analysis shows that stable stars have (wc, Mbar) values that lie between the incompressible and compressible limits <ref:2603.26973#pg2> Astrophysical constraints effectively impose a lower cut on the minimum value of Mbarmax <ref:2603.

Improvements for AI systems

  1. Precise dynamical systems modeling for stellar evolution: The improved AI system can explicitly recast stellar structure equations into a two-dimensional, non-autonomous flow in enthalpy space to identify fixed points and analyze their local stability using the Jacobian matrix (Eq. 34). This allows the system to predict the spiraling and turnover behavior of the mass–radius curve based on EoS characteristics, as shown in Section III.

  2. Universal scaling relation prediction: The AI can utilize Equation (48) to derive compact scaling relations for maximum mass, specifically identifying that the maximum-mass configuration has a trajectory (˜e(ln h), v(ln h)) that, by the time sound speed drops to a “small” value, at ln h = ln h0, has evolved to very nearly its maximum value of v√˜e, thereby predicting Mmax based on high-density stiffness and an exit density e0.

  3. Differentiation between termination mechanisms: The system can distinguish between relativistic and Newtonian instabilities by analyzing the dynamical system's behavior in different limits; it will predict that the spiral mechanism that organizes relativistic termination almost never controls the end of Newtonian stellar sequences (Section III A).

  4. Diagnostic for strong phase transitions: The AI can serve as a diagnostic tool for EoS physics by determining when a strong first-order phase transition can shift the effective fixed point and substantially alter the termination of the spiral (Section III C), providing a clean signal for such microphysics.

  5. Quasi-universal relation application in astrophysics: The system can apply Equation (86) to current astrophysical constraints, mapping observed pulsar properties onto the theoretical Mmax–wc relation to test EoS models and constrain high-density stiffness parameters like c2s,typ and wtyp (Section V).

Abstract

We reformulate the stellar structure equations in the language of dynamical systems and show that the maximum mass of stellar sequences arises from the existence of a fixed point in the relativistic regime. In an appropriate representation of the Tolman-Oppenheimer-Volkoff equations, this fixed point becomes manifest and is directly associated with the turnover of the mass-radius curve. The existence of a fixed point implies an effective reduction in dimensionality near the onset of instability, which provides a simple explanation for several equation-of-state-insensitive relations and predicts new ones. In the weakly relativistic limit, we identify a complementary universal structure shared by stellar sequences at their maximum mass, which we term the "compressible limit," and derive distinct universal relations governing the maximum mass in the Newtonian and post-Newtonian regimes. Combining these theoretical results with current astrophysical constraints, we show that the J0740+6620 pulsar is unlikely to lie near the Tolman-Oppenheimer-Volkoff maximum mass unless the equation of state exhibits a strong first-order phase transition at densities just above its central density.

Sources

Related papers