⚛️ general relativity

Low-Mass Neutron Stars and Effective Phase Transitions from a Hybrid Van der Waals-Polytropic Equation of State

This paper presents an analytically tractable hybrid Van der Waals-polytropic equation of state that successfully reproduces phase-transition signatures and predicts low-mass neutron star configurations within the range of 0.99 to 2.05 solar masses.

Original authors: P. H. F. Arruda, S. D. Campos

Published 2026-06-26

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: P. H. F. Arruda, S. D. Campos

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

Imagine a neutron star not as a solid ball of rock, but as a giant, cosmic pressure cooker. Inside, matter is squeezed so tightly that a single teaspoon would weigh a billion tons. For decades, scientists have tried to write a "recipe" (called an Equation of State) to describe how this matter behaves under such extreme pressure.

This paper by Arruda and Campos is like a chef trying to create a simplified, easy-to-read recipe book for these cosmic pressure cookers. Instead of using complex, messy ingredients that require a supercomputer to mix, they built a model using two simple, smooth formulas that snap together.

Here is the breakdown of their work using everyday analogies:

1. The Two-Part Recipe (The Core and the Crust)

The authors split the neutron star into two distinct zones, like a layered cake:

The Crust (The Outer Layer): This is the cooler, less dense outer shell. They describe this using a "polytropic" rule, which is a fancy way of saying the pressure rises in a predictable, standard way as you squeeze it.
The Core (The Center): This is the super-hot, ultra-dense heart of the star. Here, they use a modified version of the Van der Waals equation. You might know this from chemistry class as the formula that explains how real gases (like the air in a balloon) behave when they get crowded. It accounts for the fact that gas molecules take up space and push against each other.

By stitching these two recipes together, they created a "hybrid" model that is mathematically simple enough to solve with pen and paper (mostly) but complex enough to show interesting behavior.

2. The "Phase Transition" Mystery

In physics, a "phase transition" is like water turning into ice. It's a sudden, dramatic change in how matter behaves. Scientists suspect that deep inside neutron stars, matter might undergo a similar shift—perhaps turning from normal nuclear matter into something exotic like "quark soup."

The authors wanted to see if their simple two-part recipe could mimic this dramatic shift without needing a million complex variables.

The Analogy: Imagine driving a car. Usually, pressing the gas pedal makes you go faster smoothly. But what if, at a certain speed, the car suddenly feels like it's driving through mud, then suddenly feels like it's on ice? That "bump" or change in how the car responds is like a phase transition.
The Finding: They found that by tweaking just two numbers in their core recipe (which control how much the matter "pushes back" due to heat and how much it "pushes back" due to particle interactions), they could create these "bumps." Sometimes, the matter gets "stiff" (hard to squeeze), and sometimes it gets "soft" (easy to squeeze), creating a signature that looks very much like a phase transition.

3. The "Low-Mass" Surprise

One of the biggest questions in astronomy is: What is the smallest a neutron star can be?

The Expectation: Most models suggest neutron stars are heavy, usually at least 1.4 times the mass of our Sun.
The Discovery: The authors' model showed that if you adjust the "ingredients" just right, you can build stable neutron stars that are surprisingly light—ranging from about 0.99 to 2.05 times the mass of the Sun.
The Metaphor: Think of it like building a house. Most people think a house needs a massive foundation to stand up. This paper says, "Actually, if you arrange the bricks in this specific, curved way, you can build a sturdy house that is much smaller than we thought possible." This helps explain recent observations of very light neutron stars that were previously a mystery.

4. The "Blurred" Edge

In their model, the line between the crust and the core isn't a sharp, hard wall like the edge of a table. Instead, they describe it as a blurred transition zone.

The Analogy: Imagine dipping a sponge into water. The sponge doesn't go from "dry" to "wet" instantly; there's a gradient where it gets progressively wetter.
The Result: The authors found that in this "blurred" zone, the density changes smoothly. They even calculated a "conductivity" factor (a measure of how easily the transition happens) that shows the core and crust talk to each other gently rather than abruptly.

5. The Chemical Potential "Fingerprint"

To prove their model was working, they looked at something called "chemical potential." In simple terms, this is a measure of how much "energy cost" it takes to add more matter to the star.

The Expectation: If there were a true, sharp phase transition (like water freezing), this cost would stay flat for a while (a plateau).
The Reality: They didn't see a flat plateau. Instead, they saw a smooth curve with a distinct "kink" or change in slope deep inside the star.
The Conclusion: This tells us that while the star isn't undergoing a sudden, sharp explosion of change, it is experiencing a subtle, smooth shift in its internal structure deep in the core, far away from the surface.

Summary

This paper is a "proof of concept." The authors didn't try to recreate the exact, messy physics of the real universe. Instead, they built a simplified, mathematical sandbox. They showed that even with a very simple, two-part recipe, you can naturally generate:

Low-mass neutron stars (solving a recent observational puzzle).
Signatures of phase transitions (showing how matter changes behavior deep inside).
Smooth transitions between layers (rather than sharp, unrealistic walls).

Their work suggests that we don't need incredibly complex models to understand these phenomena; sometimes, a few well-chosen numbers in a simple equation can reveal the hidden secrets of the densest objects in the universe.

Technical Summary: Low-Mass Neutron Stars and Effective Phase Transitions from a Hybrid Van der Waals–Polytropic Equation of State

Problem Statement
The internal structure and Equation of State (EoS) of neutron stars (NS) remain central challenges in high-energy and nuclear astrophysics. While observational data has confirmed the existence of massive pulsars ( $\sim 2 M_\odot$ ), the lower mass limit of neutron stars remains an open question, with recent candidates suggesting masses as low as $\sim 0.77\text{--}0.98 M_\odot$ . A primary difficulty in modeling these objects is that realistic microphysical EoS models often involve complex many-body techniques and large parameter sets, making them analytically intractable. This complexity hinders the systematic identification of specific components driving softening, stiffening, or phase-transition-like behavior. Furthermore, while phase transitions in dense matter are theoretically expected, their signatures in macroscopic observables are often obscured in fully numerical models. The authors aim to address whether a simplified, analytically tractable EoS can reproduce key phenomenology associated with phase transitions and support low-mass neutron star configurations without relying on detailed microphysics.

Methodology
The authors construct a piecewise, hybrid Equation of State that couples a modified Van der Waals (vdW)-type core to a polytropic crust. The model is designed to be analytically tractable while incorporating nonlinear density dependence.

Equation of State Construction:
- Core ( $\varepsilon > \varepsilon_0$ ): Modeled by a modified vdW expression:
  $p(\varepsilon) = \frac{\alpha_1}{(\varepsilon_0/\varepsilon) - 1} \left(\frac{\varepsilon}{\varepsilon_0}\right)^{\tau_1} + \beta_1 \left(\frac{\varepsilon}{\varepsilon_0}\right)^{\sigma_1}$
  Here, $\alpha_1$ controls repulsive thermal contributions, $\beta_1$ represents interaction-driven high-density terms, and exponents $\tau_1$ and $\sigma_1$ govern the density dependence of these terms. The temperature is assumed to have a mild density dependence $T(\varepsilon) \propto \varepsilon^{\tau_1}$ .
- Crust ( $\varepsilon \le \varepsilon_0$ ): Modeled by a two-term polytropic form:
  $p(\varepsilon) = \alpha_2 \left(\frac{\varepsilon}{\varepsilon_0}\right)^{\tau_2} + \beta_2 \left(\frac{\varepsilon}{\varepsilon_0}\right)^{\sigma_2}$
- Interface: The transition between the core and crust is treated as a "blurred" region around a radius $r_1$ , where pressure and energy density are matched approximately rather than enforcing a sharp discontinuity.
Theoretical Framework:
- Thermodynamic Constraints: The model imposes thermal stability ( $dp/d\varepsilon \ge 0$ ) and causality ( $v_s^2 \le c^2$ ).
- Critical Densities: The authors derive conditions for critical energy densities ( $\varepsilon_c$ ) where the curvature of the $p(\varepsilon)$ relation changes. These are determined by the parameters $(\tau_1, \sigma_1)$ .
- Non-Relativistic Limit: The authors derive generalized Lane–Emden equations for the core and crust. Near the interface, they define an effective coupling function $A(\eta)$ (acting as a density-dependent conductivity) to describe the transition layer as a diffusive boundary rather than a discontinuity.
- Relativistic Treatment: The Tolman–Oppenheimer–Volkoff (TOV) equations are integrated numerically across various parameter regimes defined by the critical density analysis.
Diagnostic Tools:
- The study analyzes mass–radius ( $M-R$ ) sequences.
- The baryonic chemical potential $\mu(\varepsilon)$ is computed to identify signatures of phase transitions (e.g., plateaus or changes in slope).

Key Contributions and Results

Parameter Space Classification: The authors map the $(\tau_1, \sigma_1)$ parameter space to identify distinct regimes based on the existence and number of positive critical densities ( $\varepsilon_c$ ).
- Regimes with no positive roots correspond to monotonic, stiff EoS.
- Regimes with one root show mild curvature changes (softening).
- Regimes with two positive roots exhibit "spinodal-like" behavior with two inflection points, mimicking the thermodynamic signatures of a first-order phase transition (metastability and mixed phases) without requiring a true microphysical phase transition.
Low-Mass Neutron Star Configurations: Numerical integration of the TOV equations yields stable stellar configurations with gravitational masses in the range of $0.99\text{--}2.05 M_\odot$ . Specifically, the model successfully reproduces configurations with masses around $1.0\text{--}1.1 M_\odot$ and radii $\approx 11$ km, consistent with recent low-mass NS candidates.
Nature of the "Phase Transition":
- The analysis of the chemical potential $\mu(\varepsilon)$ reveals that while the model exhibits strong nonlinearities and curvature changes, it does not produce the flat plateaus characteristic of a textbook first-order phase transition.
- Instead, the "transition-like" behavior manifests as pronounced changes in the slope of $\mu(\varepsilon)$ and $dp/d\varepsilon$ at densities significantly higher than the nominal core-crust matching point ( $\varepsilon_0$ ).
- The "physical" transition region is driven by the competition between the vdW-like core terms and is located deep within the stellar core, not at the artificial piecewise junction.
Interface Dynamics: The non-relativistic analysis demonstrates that the core-crust interface behaves as a thin diffusive layer governed by the effective coefficient $A(\eta)$ . In regimes with two critical densities, $A(\eta)$ develops a minimum near the interface, indicating a localized reduction in effective conductivity where the competition between pressure terms is strongest.

Significance and Claims
The paper claims that its primary significance lies in providing a controlled, fully analytic testbed for exploring the emergence of phase-transition signatures and low-mass configurations. The authors explicitly state that the model is not intended to reproduce detailed microphysics or specific nuclear interactions. Instead, it serves to demonstrate that:

A small number of parameters in a simplified, piecewise EoS can control the emergence of complex thermodynamic behaviors (such as effective phase transitions) and low-mass NS configurations.
The location of thermodynamic signatures (diagnosed via chemical potential) does not necessarily coincide with the structural junction of a piecewise EoS; rather, it reflects internal nonlinearities at high densities.
Low-mass neutron stars ( $M \lesssim 1.1 M_\odot$ ) are compatible with simplified phenomenological models, suggesting their existence cannot be ruled out solely by gross structural arguments, though their actual formation depends on supernova dynamics and microphysics.

The work concludes that while the model is phenomenological, it effectively traces the impact of nonlinearities in the pressure-energy density relation from the EoS to global stellar structure, offering a transparent framework for investigating the conditions under which phase-transition-like features arise in analytic models.