Phase transition structure of scalarized neutron stars: the effect of rotation and linear coupling

This paper extends the Landau theory analysis of scalarization phase transitions in neutron stars to include general linear and quadratic couplings and stellar rotation, demonstrating that the former reveals overlooked solution branches while the latter primarily shifts transition masses to higher values without altering the qualitative phase structure.

The Gist

Imagine the universe as a giant, cosmic kitchen. In this kitchen, the chefs are neutron stars—the incredibly dense, heavy leftovers of exploded stars. For decades, physicists thought these stars followed a strict, boring recipe called General Relativity (Einstein's theory of gravity). They believed that no matter how you cooked them, they would always taste the same.

But recently, scientists discovered a secret ingredient: a mysterious "flavor" called a scalar field. When added to certain neutron stars, this ingredient causes a dramatic transformation called Spontaneous Scalarization. It's like a plain vanilla cake suddenly turning into a rich, chocolate lava cake just because the oven got hot enough.

This paper is about exploring exactly how and when this transformation happens, and what happens if you add two new variables to the recipe: rotation (spinning the star) and a linear twist (a new way the flavor mixes).

Here is the breakdown of their findings using simple analogies:

1. The Phase Transition: A Light Switch vs. A Dimmer

In the past, scientists thought this transformation was like a dimmer switch. You slowly turn up the "gravity heat" (mass), and the star slowly, gradually starts to develop this new flavor. It's a smooth, continuous change (a "second-order" transition).

However, this paper confirms a newer, more exciting idea: sometimes, it's not a dimmer switch; it's a light switch.

• The "First-Order" Jump: The star sits there looking normal, and then SNAP! It instantly flips into a highly flavored state.
• The "Metastable" Zone: Even crazier, there is a middle ground where the star is undecided. It could be normal or flavored, and both states are stable. It's like a ball sitting on a hilltop; a tiny nudge could send it rolling down into the "flavored" valley or back to the "normal" valley. This creates a situation where two different types of stars could exist with the exact same amount of matter.

2. The New Ingredient: The "Linear" Twist

The original recipes for these stars only used a "quadratic" ingredient (where the flavor grows with the square of the amount). This paper asks: "What if we add a linear ingredient?"

• The Analogy: Imagine you are mixing paint. The old recipe said, "Add red paint, and the color gets deeper." The new recipe says, "Add red paint, but also add a little bit of blue paint automatically."
• The Result: This "linear twist" (represented by the symbol $\alpha$ in the paper) breaks the symmetry. In the old world, the star could be "flavored" in two opposite directions (positive or negative). With the new twist, one direction becomes much stronger, and the other becomes weak or disappears entirely.
• The Surprise: The authors used a "Landau Model" (a fancy physics map) to predict that this twist could change the number of possible star types from five down to just one. It's like discovering that a menu with five different dishes suddenly only has one option left if you change the spice level. This helps them find "hidden" stars that computer simulations usually miss because they are looking in the wrong places.

3. The Spin: Spinning the Cake

Neutron stars in real life spin incredibly fast (some hundreds of times a second). The paper asked: "Does spinning the star change the recipe?"

• The Analogy: Imagine a spinning ice skater. When they spin, they feel heavier and their shape changes.
• The Finding: Spinning the star does make the "flavor jump" happen at a higher mass. It's like you need a bigger, heavier cake before the light switch flips.
• The Catch: While spinning helps, it doesn't help enough. The "flavor jump" still mostly happens for very small, light stars (less than the mass of our Sun). Since most real neutron stars we see are heavy (1.4 times the Sun's mass), these spinning stars don't quite solve the problem of making the phenomenon relevant to the heavy stars we actually observe in the sky.

4. Why This Matters

The main takeaway is that the universe is more complex than we thought.

• Hidden Diversity: By understanding the "phase transition" (the switch from normal to flavored), the scientists found that there are many more types of neutron stars hiding in the math than we previously realized.
• A New Tool: They showed that using "phase transition" thinking is a better way to find these stars than just brute-forcing the math with computers. It's like having a treasure map instead of just digging randomly in the sand.
• Observation: While the spinning stars are a bit more massive, they are still likely too light to be the main stars we see today. However, the discovery of some very light neutron stars recently suggests that we might actually be able to spot these "flavor-switching" events in the future.

Summary

This paper is a guide to the "flavor" of neutron stars. It tells us that:

1. The switch from normal to "flavored" stars is often a sudden, dramatic jump, not a slow change.
2. Adding a new type of interaction (the linear term) changes the menu, sometimes removing options and sometimes revealing hidden ones.
3. Spinning the stars makes them heavier, but not heavy enough to match most of the stars we see in the sky yet.

Ultimately, the authors are saying: "Don't just look for the stars you expect; look for the ones that are hiding in the gaps of the phase transition, because that's where the new physics is."

Technical Summary

1. Problem Statement

The paper addresses the phenomenon of spontaneous scalarization in scalar-tensor theories (STT) of gravity, specifically focusing on neutron stars. While scalarization was historically viewed as a continuous (second-order) phase transition from General Relativity (GR) to a scalarized state, recent work (e.g., Unlütürk et al.) has established that first-order (discontinuous) transitions are actually more common in the parameter space of the Damour-Esposito-Farèse (DEF) model.

However, previous studies were limited by two main constraints:

1. Coupling Function: They primarily considered quadratic scalar-matter coupling terms ($\alpha=0$ in the conformal factor), ignoring linear terms. In generic scalar-tensor theories, linear coupling terms are expected to exist.
2. Symmetry: They assumed spherically symmetric (static) stars, neglecting the ubiquitous rotation of astrophysical neutron stars.

The authors aim to investigate how the inclusion of a linear coupling term ($\alpha \neq 0$) and stellar rotation ($\Omega \neq 0$) alters the phase transition structure, the number of equilibrium solutions, and the astrophysical relevance of these phenomena.

2. Methodology

Theoretical Framework

• Action: The study utilizes the Einstein-frame action for scalar-tensor theories with a massive scalar field $\phi$ and a potential $V(\phi) = \frac{1}{2}m_\phi^2 \phi^2$.
• Conformal Factor: The coupling between the scalar field and matter is defined via the conformal factor $A(\phi) = e^{\alpha\phi + \frac{1}{2}\beta\phi^2}$.
  • $\beta$: Controls the quadratic coupling (responsible for tachyonic instability and scalarization).
  • $\alpha$: Controls the linear coupling (breaks the $\phi \to -\phi$ symmetry).
• Matter: Neutron stars are modeled as perfect fluids using the MPA1 equation of state (EOS).
• Geometry: Both static and rapidly rotating (uniform rotation) configurations are solved using quasi-isotropic coordinates.

Analytical Approach: Landau Theory

The authors employ Landau theory of phase transitions to construct a phenomenological ansatz for the ADM mass ($M_{ADM}$) as a function of baryon mass ($M_b$) and an order parameter $Q$ (representing scalar field strength):
$$M_{ADM}(M_b, Q) = M_0(M_b) + a(M_b)Q^2 + \frac{1}{2}b(M_b)Q^4 + \frac{1}{3}c(M_b)Q^6$$

• $\alpha = 0$: The theory is symmetric ($\phi \to -\phi$), so only even powers of $Q$ exist.
• $\alpha \neq 0$: The symmetry is broken, introducing odd powers of $Q$ (specifically a linear term). This fundamentally alters the energy landscape, potentially reducing the number of equilibrium branches.

Numerical Implementation

• Static Stars: Solved using a 1D shooting method code.
• Rotating Stars: Solved using a modified RNS code (based on the Komatsu-Eriguchi-Hachisu scheme) adapted for scalar-tensor theories.
• Validation: The authors cross-verified results between the shooting method and a relaxation-based code to ensure accuracy, particularly for high coupling values ($|\beta|$) and scalar masses.

3. Key Contributions

1. Generalization of Coupling: The paper is the first to systematically explore the phase transition structure of scalarization with linear coupling terms ($\alpha \neq 0$) in the conformal factor.
2. Systematic Solution Discovery: The authors demonstrate that the Landau phenomenological model acts as a powerful predictive tool. It allows researchers to systematically identify multiple branches of scalarized solutions (including those with positive and negative scalar field signs) that are frequently missed by standard numerical searches which tend to follow a single branch.
3. Rotation Effects: The study quantifies how uniform rotation shifts the critical masses for phase transitions, addressing whether rotation can make first-order scalarization observable in typical neutron star mass ranges.

4. Key Results

A. Effect of Linear Coupling ($\alpha \neq 0$)

• Symmetry Breaking: The introduction of $\alpha \neq 0$ breaks the $\phi \to -\phi$ symmetry. Consequently, the "unscalarized" GR solution ($Q=0$) no longer exists as an exact solution; it becomes a slightly scalarized state.
• Reduction of Solution Branches:
  • Deformed Second-Order: For small $|\beta|$, the system typically exhibits 3 equilibrium solutions (1 unstable, 2 stable) at low $\alpha$. As $|\alpha|$ increases, the number of solutions drops to 1 (a single stable, highly scalarized branch).
  • Deformed First-Order: For large negative $\beta$ (where first-order transitions dominate in the $\alpha=0$ case), the system can exhibit up to 5 equilibrium solutions (3 stable, 2 unstable) at low $\alpha$. As $|\alpha|$ increases, this reduces to 3, and eventually to 1.
• Phenomenological Prediction: The Landau model successfully predicted the existence of these multiple branches. Without this theoretical guidance, numerical searches would likely miss the "hidden" branches (e.g., the positively scalarized branch or the metastable states).

B. Effect of Rotation ($\Omega \neq 0$)

• Mass Shift: Rotation increases the stellar mass for a given central energy density. This shifts the bifurcation points (where scalarization begins) and the local minima of the mass curves to higher baryon masses.
• Quantitative Impact: For the parameters studied, rotation increased the critical mass by up to ~5% (for angular momentum $J \approx 0.2$).
• Qualitative Impact: Rotation did not qualitatively change the phase transition structure (i.e., it did not turn a first-order transition into a second-order one, nor did it create new types of transitions).
• Astrophysical Relevance: While rotation increases the mass of scalarized stars, the shift is insufficient to bring the first-order transition mass (typically $< 1 M_\odot$) into the range of the most massive observed neutron stars ($\sim 2 M_\odot$). However, it makes the phenomenon slightly more relevant to lower-mass neutron stars (e.g., the recently discovered $0.77 M_\odot$ star).

C. Stability Analysis

• Stability Criteria: Stability was determined by the sign of $dM_b/d\epsilon_c$ (where $\epsilon_c$ is central energy density).
  • $dM_b/d\epsilon_c > 0$: Hydrostatically stable.
  • $dM_b/d\epsilon_c < 0$: Unstable.
• Ground State: The global energy minimum (ground state) was consistently found to be the branch with the most negative scalar field amplitude.
• Metastability: In the first-order transition regime, multiple locally stable configurations exist at the same baryon mass. Transitions between these states could be triggered by astrophysical perturbations, potentially leading to observable gravitational wave signatures.

5. Significance and Conclusion

• Theoretical Insight: The paper confirms that spontaneous scalarization is a rich phase transition phenomenon. The inclusion of linear coupling terms ($\alpha$) drastically simplifies the solution space at high values, effectively suppressing the complex multi-branch structure seen in pure DEF models ($\alpha=0$).
• Methodological Advancement: The study highlights the critical importance of combining phenomenological Landau analysis with numerical relativity. Relying solely on numerical shooting methods risks missing entire branches of solutions, whereas the Landau approach provides a roadmap for a complete exploration of the solution space.
• Observational Implications:
  • The existence of metastable states and discontinuous transitions offers new channels for observation, such as gravitational wave emission during transitions between scalarized and unscalarized (or differently scalarized) states.
  • While rotation enhances the mass of scalarized stars, it does not fully resolve the tension between the low-mass theoretical predictions of first-order scalarization and the high masses of most observed neutron stars. However, the discovery of low-mass neutron stars ($< 1 M_\odot$) suggests these effects remain astrophysically relevant.
• Future Directions: The authors suggest extending this analysis to black holes, different equations of state (including phase transitions in the core), and differential rotation to further explore the astrophysical signatures of scalar-tensor gravity.