pith. sign in

arxiv: 1907.01330 · v1 · pith:IZ3PDV2Inew · submitted 2019-07-02 · 🌌 astro-ph.HE · nucl-th

Gravitational waves from phase transition of NS to QS

Prasad R , Ritam Mallick This is my paper

Pith reviewed 2026-05-25 10:55 UTC · model grok-4.3

classification 🌌 astro-ph.HE nucl-th
keywords gravitational wavesphase transitionneutron starhybrid starquark mattercombustionshock wave2D simulation
0
0 comments X

The pith

The phase transition from a neutron star to a hybrid star generates gravitational waves with amplitude 10^{-22} lasting tens of microseconds.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper runs a two-dimensional simulation of combustion inside a neutron star that converts nuclear matter to two-flavor quark matter. A central density fluctuation creates a shock that starts the combustion front, which travels outward but halts at 6 km where nuclear matter remains more stable than quark matter; the shock wave itself continues to the surface. The resulting gravitational-wave strain reaches 10^{-22}, persists for tens of microseconds, oscillates mildly, and shows power-spectrum peaks at high frequencies. A reader would care because this pattern supplies a distinct observational marker that could confirm the phase change occurred and reveal whether the star ends as a hybrid star.

Core claim

We perform a 2-d simulation of combustion of neutron star (NS) to hybrid star (HS). We assume that a sudden density fluctuation at the center of the NS initiates a shock discontinuity near the center of the star. This shock discontinuity deconfines NM to 2-f QM, initiating combustion of the star. This combustion front propagates from the center to the surface converting NM to 2-f QM. This combustion stops at a radius of 6 km inside the star, as at this density the NM is much stable than QM. Beyond 6 km although the combustion stops but the shock wave propagates to the surface. We study the gravitational wave signal for such a PT of NS to HS. We find that such PT has unique GW strain of anam

What carries the argument

The outward-moving combustion front initiated by a central shock discontinuity that converts neutron matter to two-flavor quark matter and halts at 6 km while the continuing shock generates the waves.

If this is right

  • The gravitational-wave strain reaches an amplitude of 10^{-22}.
  • The signals persist for a few tens of microseconds and display small oscillations.
  • The power spectrum shows distinct peaks at fairly high frequencies.
  • The conversion leaves a unique signature that can identify the phase transition and the final state of the star.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • High-frequency gravitational-wave observatories could search for these brief bursts to test whether any neutron stars undergo this conversion.
  • Varying the stopping radius or initial fluctuation strength in follow-up simulations would map how sensitive the waveform is to those choices.
  • The same combustion setup might be compared with other quark-matter formation channels to isolate what makes this gravitational-wave pattern distinctive.

Load-bearing premise

A sudden central density fluctuation creates a shock that deconfines nuclear matter into two-flavor quark matter and drives combustion that stops at 6 km radius.

What would settle it

Absence of short-duration gravitational-wave bursts with amplitude near 10^{-22} and high-frequency spectral peaks from any observed neutron-star event would contradict the predicted signal.

Figures

Figures reproduced from arXiv: 1907.01330 by Prasad R, Ritam Mallick.

Figure 1
Figure 1. Figure 1: FIG. 1. (Color online) The chemical potential is shown as [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The density variation along the star radial distance [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (Color online) The evolution of density as a function [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (Color online) The density fluctuation results also i [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. (Color online) The matter velocity at initial time is [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The velocity of the combustion front and the shock [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. PT in neutron star (PLZ-M1) for different time instanc [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. The change in the Quadruple moment of the star as [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. The amplitude spectrcal density is shown as a [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Amplitude spectral density for GW produced from [PITH_FULL_IMAGE:figures/full_fig_p010_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. We show the strain for the fours stage of combusting N [PITH_FULL_IMAGE:figures/full_fig_p011_13.png] view at source ↗
read the original abstract

In this article, we perform a 2-d simulation of combustion of neutron star (NS) to hybrid star (HS). We assume that a sudden density fluctuation at the center of the NS initiates a shock discontinuity near the center of the star. This shock discontinuity deconfines NM to 2-f QM, initiating combustion of the star. This combustion front propagates from the center to the surface converting NM to 2-f QM. This combustion stops at a radius of $6 km$ inside the star, as at this density the NM is much stable than QM. Beyond $6 km$ although the combustion stops but the shock wave propagates to the surface. We study the gravitational wave signal for such a PT of NS to HS. We find that such PT has unique GW strain of amplitude $10^{-22}$. These signals last for few tens of $\mu s$ and shows small oscillating behaviour. The power spectrum consists of peaks and at fairly high frequency range. The conversion to NS to HS has a unique signature which would help in defining the PT and the fate of the NS.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript presents a 2D hydrodynamic simulation of the phase transition from a neutron star (NS) to a hybrid star (HS) via combustion. A central density fluctuation initiates a shock discontinuity that deconfines nuclear matter (NM) to two-flavor quark matter (QM); the combustion front propagates outward but is stated to terminate at a radius of 6 km where NM is more stable than QM, while the shock continues to the surface. The resulting gravitational-wave (GW) signal is extracted and reported to have a unique strain amplitude of 10^{-22}, duration of a few tens of microseconds with small oscillations, and a power spectrum with peaks at high frequencies. The authors conclude that this provides a distinctive signature for identifying NS-to-HS conversions.

Significance. If the modeling assumptions prove robust and the simulation is shown to be converged, the work would supply a concrete, falsifiable GW template for a specific class of compact-object phase transitions, which could be tested against high-frequency GW data from future detectors. The approach of coupling combustion hydrodynamics directly to GW extraction in the NS context is relevant to multi-messenger astrophysics, though the current lack of sensitivity studies limits the strength of the uniqueness claim.

major comments (3)
  1. [Abstract] Abstract: the combustion termination radius of 6 km is asserted because 'at this density the NM is much stable than QM,' yet no explicit stability criterion, bag constant, chemical-potential comparison, or EOS table is supplied to locate this radius for the chosen initial NS profile; this choice directly fixes the converted mass and therefore sets the reported GW amplitude, duration, and spectrum.
  2. [Abstract] Abstract: the 2D simulation reports no convergence tests, resolution studies, or sensitivity analysis with respect to the initiation mechanism (sudden central density fluctuation) or the 6 km cutoff; without these, the headline values (strain 10^{-22}, duration tens of μs, high-frequency peaks) cannot be established as robust outputs of the phase-transition dynamics rather than artifacts of the chosen parameters.
  3. [Abstract] Abstract: the claim that the GW signal is 'unique' and 'would help in defining the PT and the fate of the NS' rests on a single run with an ad-hoc termination radius; a quantitative demonstration that the waveform features survive variations in the cutoff radius or EOS parameters is required to support the uniqueness assertion.
minor comments (2)
  1. [Abstract] Abstract: '2-f QM' is used without prior definition; 'shows small oscillating behaviour' and 'consists of peaks and at fairly high frequency range' are grammatically unclear and should be rephrased for precision.
  2. [Title] Title vs. abstract: the title refers to 'NS to QS' while the text discusses conversion to a hybrid star (HS); the distinction should be clarified or the title adjusted.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive comments. We agree that the manuscript requires additional details on the EOS and stability criterion used for the termination radius, as well as convergence and sensitivity studies, to support the reported GW signals and uniqueness claim. We will revise the manuscript to address these points. Point-by-point responses to the major comments are provided below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the combustion termination radius of 6 km is asserted because 'at this density the NM is much stable than QM,' yet no explicit stability criterion, bag constant, chemical-potential comparison, or EOS table is supplied to locate this radius for the chosen initial NS profile; this choice directly fixes the converted mass and therefore sets the reported GW amplitude, duration, and spectrum.

    Authors: We agree that the current manuscript does not supply the explicit stability criterion or EOS parameters. In the revised version we will add the bag constant value, the specific EOS tables for NM and QM, and the chemical-potential comparison that identifies the 6 km radius at which NM becomes more stable than QM for the chosen initial profile. This will be placed in the methods section so that the converted mass and resulting GW quantities can be reproduced. revision: yes

  2. Referee: [Abstract] Abstract: the 2D simulation reports no convergence tests, resolution studies, or sensitivity analysis with respect to the initiation mechanism (sudden central density fluctuation) or the 6 km cutoff; without these, the headline values (strain 10^{-22}, duration tens of μs, high-frequency peaks) cannot be established as robust outputs of the phase-transition dynamics rather than artifacts of the chosen parameters.

    Authors: The original submission indeed contains no convergence or resolution studies. We will incorporate these in the revision by presenting results from at least two additional grid resolutions and by testing a modest range of central density fluctuation amplitudes. The GW amplitude, duration, and high-frequency spectral peaks will be shown to remain consistent within the reported precision. revision: yes

  3. Referee: [Abstract] Abstract: the claim that the GW signal is 'unique' and 'would help in defining the PT and the fate of the NS' rests on a single run with an ad-hoc termination radius; a quantitative demonstration that the waveform features survive variations in the cutoff radius or EOS parameters is required to support the uniqueness assertion.

    Authors: We acknowledge that the uniqueness statement is based on a single simulation. While a full parameter survey lies outside the scope of this initial study, the revision will qualify the claim by noting that the short duration and high-frequency content arise directly from the rapid combustion front. We will add a brief qualitative discussion of how modest changes in cutoff radius affect the signal and will label the uniqueness as suggestive rather than exhaustive. revision: partial

Circularity Check

0 steps flagged

No circularity: GW signals are direct numerical outputs from stated model assumptions

full rationale

The paper performs a 2D hydrodynamical simulation of combustion initiated by a central density fluctuation, with the combustion front explicitly stated to halt at 6 km because NM is more stable than QM at that density. The reported GW strain amplitude of 10^{-22}, duration of tens of μs, and high-frequency spectral peaks are computed outputs of this evolution. No equations, fitted parameters, or self-citations are shown to reduce the GW quantities to the inputs by construction; the stopping radius is presented as an input assumption rather than a derived result that is then renamed as a prediction. The derivation chain is therefore self-contained as a forward simulation under the given model.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The model depends on domain assumptions about combustion initiation and termination that are introduced to enable the simulation rather than derived from first principles or external benchmarks.

free parameters (1)
  • combustion termination radius = 6 km
    Set to 6 km on the basis of the density at which NM is stated to be more stable than QM.
axioms (2)
  • domain assumption A sudden central density fluctuation initiates a shock discontinuity that deconfines nuclear matter to two-flavor quark matter.
    This is the stated trigger for the combustion front in the abstract.
  • domain assumption The combustion front stops at 6 km because nuclear matter is much more stable than quark matter at that density.
    This fixes the radial extent of the phase transition.

pith-pipeline@v0.9.0 · 5715 in / 1722 out tokens · 76223 ms · 2026-05-25T10:55:49.735029+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages

  1. [1]

    The star merger has revealed itself in the form of GW, neutrino production and also in other various spectra of electromagnetic range

    of a neutron star merger has shed new light not only in the field of GW astronomy but also to multimessenger as- tronomy. The star merger has revealed itself in the form of GW, neutrino production and also in other various spectra of electromagnetic range. With the rejuvenation of GW astronomy and the GW detection, several new gravitational wave detectors ...

    work page

  2. [2]

    This smoothing of variation suppresses the numerical noise in gravitational wave strain

    filter to perform differentiation, smoothing the data by choosing appropriate window size such that it dimin- ishes variations in data in 10 − 7 s, and differentiation is carried out for variations in quadrupole moment which are above 10 − 6s. This smoothing of variation suppresses the numerical noise in gravitational wave strain. We have shown curves for th...

    work page 2015

  3. [3]

    B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Phys. Rev. Lett. 116, 061102 (2016)

    work page 2016

  4. [4]

    B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Phys. Rev. Lett. 119, 161101 (2017)

    work page 2017

  5. [5]

    Annala, T

    E. Annala, T. Gorda, A. Kurkela, and A. Vuorinen, Phys. Rev. Lett. 120, 172703 (2018)

    work page 2018

  6. [6]

    Margalit and B

    B. Margalit and B. D. Metzger, Astrophys. J. Lett. 850, L19 (2017)

    work page 2017

  7. [7]

    Radice, Astrophys

    D. Radice, Astrophys. J. Lett. 838, L2 (2017)

    work page 2017

  8. [8]

    M. Ruiz, S. L. Shapiro, and A. Tsokaros, Phys. Rev. D 97, 021501 (2018)

    work page 2018

  9. [9]

    Shibata, S

    M. Shibata, S. Fujibayashi, K. Hotokezaka, K. Kiuchi, K. Kyutoku, Y. Sekiguchi, and M. Tanaka, Phys. Rev. D 96, 123012 (2017)

    work page 2017

  10. [10]

    E. R. Most, L. R. Weih, L. Rezzolla, and J. Schaffner- Bielich, Phys. Rev. Lett. 120, 261103 (2018)

    work page 2018

  11. [11]

    Ferrari and R

    A. Ferrari and R. Ruffini, Astrophys. J. 158, L71 (1969)

    work page 1969

  12. [12]

    Zimmerman and E

    M. Zimmerman and E. Szedenits, Phys. Rev. D 20, 351 (1979)

    work page 1979

  13. [13]

    Chandrasekhar, Phys

    S. Chandrasekhar, Phys. Rev. Lett. 24, 611 (1970)

    work page 1970

  14. [14]

    Friedman and B

    J. Friedman and B. Schutz, Astrophys. J. 222, 281 (1978)

    work page 1978

  15. [15]

    Thorne and A

    K. Thorne and A. Campolattaro, Astrophys. J. 149, 591 (1967)

    work page 1967

  16. [16]

    Prasad & R

    R. Prasad & R. Mallick, Astrophys. J. 859, 57 (2018)

    work page 2018

  17. [17]

    Itoh, Prog

    N. Itoh, Prog. Theor. Phys. 44, 291 (1970)

    work page 1970

  18. [18]

    A. R. Bodmer, Phys. Rev. D 4, 1601 (1971)

    work page 1971

  19. [19]

    Witten, E., PRD 30, 272 (1984)

    work page 1984

  20. [20]

    Parenti & I

    Bombaci, I. Parenti & I. Vidana, AstroPhys. J. 614, 314 (2004)

    work page 2004

  21. [21]

    Drago, A

    A. Drago, A. Lavagno & G. Pagliara, Eur. Phys. J. A 19, 197 (2004)

    work page 2004

  22. [22]

    B. W. Mintz, E. S. Fraga, G. Pagliara & J. Schaffner- Bielich, Phys. Rev. D 81, 123012 (2010)

    work page 2010

  23. [23]

    R., Oertel, M., & Mar- gueron, J., Phys.Rev

    Gulminelli, F., Raduta, Ad. R., Oertel, M., & Mar- gueron, J., Phys.Rev. C 87, 055809 (2013)

    work page 2013

  24. [24]

    Bombaci & B

    I. Bombaci & B. Datta, AstroPhys. J. 530, L69 (2000)

    work page 2000

  25. [25]

    Berezhiani, I

    Z. Berezhiani, I. Bombaci, A. Drago, F. Frontera & A. Lavagno, AstroPhys. J. 586, 1250 (2003)

    work page 2003

  26. [26]

    Mallick & P

    R. Mallick & P. K. Sahu, Nucl. Phys. A 921, 96 (2014)

    work page 2014

  27. [27]

    L. M. Lin, K. S. Cheng, M. C. Chu & W. -M. Suen, Astrophys. J. 639, 382 (2006)

    work page 2006

  28. [28]

    E. B. Abdikamalov, H. Dimmelmeier, L. Rezzolla & J. C. Miller, Mon. Not. R. Astron. Soc. 392, 52 (2009)

    work page 2009

  29. [29]

    Reinhard, Z

    P.-G. Reinhard, Z. Phys. A 329, 257 (1988)

    work page 1988

  30. [30]

    Chodos, R

    A. Chodos, R. L. Jaffe, K. Johnson, C. B. Thorn, and V. F. Weisskopf, Phys. Rev. D 9, 3471 (1974)

    work page 1974

  31. [31]

    Bhattacharyya, S

    A. Bhattacharyya, S. K. Ghosh, P. S. Joarder, R. Mallick & S. Raha, Phys. Rev. C 74, 065804 (2006)

    work page 2006

  32. [32]

    Komatsu, Y

    H. Komatsu, Y. Eriguchi, and I. Hachisu, Mon. Not. R. Astron. Soc. 237, 355 (1989)

    work page 1989

  33. [33]

    Stergioulas and J

    N. Stergioulas and J. L. Friedman, Astrophys. J. 444, 30 6 (1995)

    work page 1995

  34. [34]

    G. B. Cook, S. L. Shapiro, and S. A. Teukolsky, Astro- phys. J. 398, 203 (1992)

    work page 1992

  35. [35]

    OConnor and C

    E. OConnor and C. D. Ott, Class. Quantum Grav., 27, 114103 (2009)

    work page 2009

  36. [36]

    Font J A, Miller M, Suen W-M and Tobias M, Phys. Rev. D. 61, 044011 (2000)

    work page 2000

  37. [37]

    Hyman J M, The method of lines solution of par- tial differential equations Technical Report COO-3077- 139 (ERDA Mathematics and Computing Laboratory, Courant Institute of Mathematical Sciences, New York University (1976))

    work page 1976

  38. [38]

    J., 462, 839 (1996)

    Romero J V, Ibez J M, Mart J M and Miralles J A, Astrophys. J., 462, 839 (1996)

    work page 1996

  39. [39]

    Colella P, & Woodward P R, J. Comput. Phys. 54, 174 (1984)

    work page 1984

  40. [40]

    A., Mueller E., Astron

    Dimmelmeier H., Font J. A., Mueller E., Astron. & As- trophys. 393, 523 (2002)

    work page 2002

  41. [41]

    & Mueller, E., Astron

    Zwerger, T. & Mueller, E., Astron. Astrophys., 320, 209 (1997)

    work page 1997

  42. [42]

    H., Teukolsky, S

    Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in C: The art of scientific programming (Cambridge, U. K.: Cambridge 13 University Press)

    work page 1992

  43. [43]

    Savitzky and M

    A. Savitzky and M. J. E. Golay, Analytical Chemistry, 36, 8 (1964)

    work page 1964

  44. [44]

    Aasiel et al., AstroPhys

    J. Aasiel et al., AstroPhys. J. 785, 119 (2014)

    work page 2014

  45. [45]

    P. D. Lansky, Pub. Astron. Soc. Aus. (PASA) 32, e034 (2015)

    work page 2015