Pith. sign in

REVIEW 3 major objections 5 minor 18 references

Ultra-dense quark pellets formed in the early universe could constitute all of dark matter without requiring new particles.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 07:10 UTC pith:C7AXKBUB

load-bearing objection Interesting scenario on paper, but the central TOV equation has a factor-of-3 error and the formation mechanism is explicitly deferred; the dark-matter claim is not supported as written. the 3 major comments →

arxiv 2607.10672 v2 pith:C7AXKBUB submitted 2026-07-12 hep-ph astro-ph.COgr-qc

Do primordial quark pellets solve the dark matter puzzle?

Domenec Espriu This is my paper
classification hep-ph astro-ph.COgr-qc
keywords primordial quark pelletsdark matterquark starsTolman-Oppenheimer-Volkoff equationPeccei-Quinn domain wallsbaryon overdensitiesstrange mattermicrolensing constraints
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper argues that primordial quark pellets — ultra-dense quark-matter mini-stars that form when the universe is about 1 GeV old — are a viable cold dark matter candidate. Solving the Tolman-Oppenheimer-Volkoff equation with a hot quark equation of state yields stable solutions capped at roughly 0.01 solar masses and 100-meter radii. Because each pellet is so heavy, only one maximal pellet per billion horizon volumes (or lighter pellets in correspondingly more horizons) would reproduce the observed dark matter density without disturbing Big Bang nucleosynthesis. The central open question is whether collapsing Peccei-Quinn domain walls or superhorizon bubbles can concentrate baryons enough to seed such objects; the paper explicitly defers a quantitative treatment of that sweep-and-collect mechanism to a companion study.

Core claim

The paper's central claim is that stable hot quark stars can exist in the early universe at T≈1 GeV, with a maximum mass of about 0.013 solar masses and a radius near 88 meters, set by balancing internal quark degeneracy pressure against the ambient radiation pressure of the surrounding plasma. The TOV integration gives a one-parameter family of solutions; the largest stable objects approach the horizon mass at formation. The required baryon overdensities, if produced by collapsing superhorizon bubbles or Peccei-Quinn domain walls, would give a steep mass function dN/dM ∝ M^-2 and formation probabilities of 10^-9 to 10^-4 per horizon volume that account for the full dark matter abundance. On

What carries the argument

The Tolman-Oppenheimer-Volkoff equation integrated with a hot, degenerate u,d,s quark equation of state, with a surface boundary condition at which the internal pressure matches the external radiation pressure at T=1 GeV. Defining x as the quark number density normalized to its surface value, the paper integrates from a central density x_c down to x=1, producing the mass-radius relation. This is paired with a cosmological sweep-and-collect mechanism: collapsing Peccei-Quinn domain walls or superhorizon bubbles concentrate baryons into rare overdensities, yielding a mass function dN/dM_PQP ∝ M_PQP^-2.

Load-bearing premise

The load-bearing premise is that collapsing Peccei-Quinn domain walls or superhorizon bubbles can concentrate baryons by a factor of up to about 10^4 relative to the horizon average at T=1 GeV; the paper explicitly states in Section III.D that the quantitative sweep-and-collect treatment is beyond its scope and deferred to a companion study.

What would settle it

Compute the baryon mass a collapsing Peccei-Quinn domain wall or superhorizon bubble can sweep up at T=1 GeV; if the maximum collected mass falls below about 10^-7 solar masses, the formation probability cannot reach the required 10^-4 per horizon volume and PQPs cannot be the dark matter.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If PQPs are the dark matter, no exotic cold-dark-matter particle is needed; the only non-Standard-Model input is a local baryon-concentration mechanism such as Peccei-Quinn domain walls.
  • Because the excess baryons are sequestered inside PQPs, the radiation-dominated expansion and Big Bang nucleosynthesis proceed as in the standard cosmology, with no entropy dilution required.
  • PQPs with masses below roughly 10^-7 solar masses evade current optical microlensing limits, opening a window in which they could constitute most of the dark matter.
  • If the formation epoch is tied to a first-order phase transition near 1 TeV, the model predicts a stochastic gravitational-wave background peaking in the mHz–Hz band, observable by future space-based interferometers.
  • Depending on the true QCD ground state, the pellets end as primordial mini neutron stars or stable strange-matter nuggets, with identical astronomical signatures.
  • The predicted mass function is steep (dN/dM ∝ M^-2), so the population is numerically dominated by the lightest pellets, making sub-Earth-mass microlensing surveys the most direct observational test.

Where Pith is reading between the lines

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

  • Extension: The rate at which collapsing domain walls sweep up baryons is the single most decisive input; a dedicated hydrodynamical calculation of sweep efficiency would either confirm or kill the scenario before any astronomical observation is needed.
  • Extension: If the Peccei-Quinn transition is the source of the overdensities, then the same sector producing axions would also seed PQPs, linking two seemingly independent dark matter candidates in a single cosmological story.
  • Extension: The TOV boundary condition depends on the formation temperature; the paper's T=400 MeV results show the maximum mass rises as temperature drops, so measuring the upper cutoff of the PQP mass function could reconstruct the formation epoch.
  • Extension: A non-detection of the predicted gravitational-wave background would not by itself falsify PQPs, since the signal is tied to a first-order phase transition rather than to the pellets themselves.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper proposes that primordial quark pellets (PQPs) — compact, hot quark-matter objects formed around T ~ 1 GeV — could account for all dark matter. The author integrates the TOV equation with an analytic hot-quark EOS and an environmental pressure boundary condition, obtaining a maximum mass of 0.013 M☉ and radius ~88 m. Using the observed DM density, the paper infers a required formation probability of 10^-9 to 10^-4 per horizon volume, and argues that a steep mass function with most PQPs below 10^-7 M☉ would evade microlensing limits. The seeding mechanism is proposed to be contracting Peccei-Quinn domain walls or superhorizon bubbles, but the quantitative treatment is explicitly deferred.

Significance. If correct, the scenario would be an interesting and testable dark-matter candidate using known QCD physics plus general relativity, with potentially observable signatures in microlensing and LISA-band gravitational waves. The paper is transparent in laying out its equations and boundary conditions, and the abundance calculation is presented as a consistency check rather than a parameter fit. However, the quantitative core is compromised by concrete errors in the TOV equation and the equation of state, and the formation mechanism, which is load-bearing for the claim that PQPs 'naturally arise', is not actually derived. As it stands, the paper does not establish the viability of the PQP scenario.

major comments (3)
  1. [II, Eq. (9)] The TOV equation as written has a factor-of-three error. For the stated EOS, P = ε/3, so the standard relativistic hydrostatic equilibrium equation dP/dr = -G(ε+P)(m+4πr³P)/(r(r-2Gm)) becomes dP/dr = -G(4ε/3)(m+4πr³P)/(r(r-2Gm)). Equation (9) instead uses 4Gε, corresponding to an effective gravitational source that is three times too strong. Because the quoted maximum mass, 0.013 M☉, the radius 88 m, and the stability boundary dM/dx_c = 0 all come from integrating Eq. (9), these central numerical results are not reliable and must be recomputed.
  2. [II, Eqs. (1)-(2)] The EOS is not internally consistent for the stated three-flavor quark gas. If n is interpreted as baryon density, as implied by the boundary value n_R ≈ 233.8 fm^-3 being described as ~1500 times nuclear saturation density, the n^{4/3} coefficient is consistent with three flavors with color degeneracy, but the T² coefficient 3^{1/3}π^{4/3} is too small by approximately a factor of three. If n is interpreted as quark number density, the n^{4/3} coefficient is incorrect. Since both P and dP/dx enter the TOV integration, the corrected EOS will change the mass-radius curve and the abundance window in Sec. III.C.
  3. [III.D] The baryon-concentration mechanism is the load-bearing ingredient for the entire scenario, but it is not quantified. The paper asserts that collapsing domain walls or superhorizon bubbles sweep up baryons by many orders of magnitude relative to the horizon average, and then asserts scaling relations dN/dR ∝ R^{-4} and dN/dM ∝ M^{-2} without a calculation of sweep efficiency, baryon capture, or survival probability. The text itself states that a full quantitative treatment is 'beyond the scope of this paper'. Without an actual derivation of the formation probability, the central claim that PQPs 'naturally arise' and could constitute all dark matter is not supported; Eq. (18) only specifies what formation probability would be needed.
minor comments (5)
  1. [Abstract / I] The abstract and introduction state T ~ 1 TeV in one place and T ~ 1 GeV elsewhere; the intended formation temperature in the body is ~1 GeV.
  2. [II, Eq. (15)] Equation (15) writes the thermal pressure term as 3.26 x^{2/3} T², but in Eq. (7) the same term was evaluated at T = 1 GeV. The cooling expression should use a dimensionless ratio such as (T/T_ref)², otherwise the equation has inconsistent units.
  3. [III.C, Eq. (17)] Equation (17) uses a comma as a decimal separator ('3,6×10^{-11}'); this should be a period for consistency.
  4. [II] The text following Eq. (8) says 'Schwarzchild' in the phrase 'The Schwarzchild radius of this object is 38 meters'; this should be 'Schwarzschild'.
  5. [Abstract / V] The abstract calls the proposal 'conservative, Standard-Model-based', but the seeding mechanism relies on Peccei-Quinn domain walls or other new physics beyond the Standard Model; the conclusions are more accurate in this respect.

Circularity Check

0 steps flagged

No significant circularity: the TOV-derived mass-radius relation and the abundance consistency check are independent of the fitted dark-matter input; the formation mechanism is explicitly deferred rather than used to derive the result.

full rationale

The paper's central derivation is the TOV integration in Sec. II. The EOS (Eqs. 1-2, 6-7) is an analytic free-quark-gas model with thermal corrections; the only boundary input is the pressure-matching condition (Eq. 3), which fixes n_R, and the temperature T=1 GeV is chosen from microphysical constraints (charm threshold, QCD transition). No parameter in this calculation is fitted to the dark-matter density; the M_max=0.013 M_sun, R=88 m result is a genuine TOV output. In Sec. III.C the observed ratio epsilon_DM/epsilon_r(T0) is used only to compute the required number of PQPs per horizon (Eqs. 16-18). The paper explicitly presents this as a 'required formation abundance' and says 'could suffice,' not as a prediction derived from a formation model. The formation mechanism (Sec. III.D) is openly deferred: 'a full quantitative treatment ... is beyond the scope of this paper and will be addressed in a companion study.' Thus the paper does not claim to derive the abundance from first principles; it derives a consistency requirement. Citations are to external work (Witten, Krnjaic-Rocha-Xiao, Vilenkin-Shellard, etc.), with no load-bearing self-citation or uniqueness theorem imported from the present author's prior work. A possible factor-of-3 issue in Eq. (9) would be a numerical/correctness concern, not a circularity, and does not affect the circularity score.

Axiom & Free-Parameter Ledger

2 free parameters · 5 axioms · 1 invented entities

The central calculation relies on a chosen EOS and formation temperature, while the viability of the scenario depends on an unquantified baryon-concentration mechanism. The visible data do not fix the EOS coefficients or the minimum PQP mass.

free parameters (2)
  • formation temperature T_form = 1 GeV (also 0.4 GeV for comparison)
    Chosen by hand to be below the charm mass and above the QCD transition; maximum mass and radius scale with this choice.
  • minimum PQP mass M_min
    The mass function dN/dM ∝ M^-2 requires a lower cutoff to compute the total dark-matter density; the paper states this cutoff is not established.
axioms (5)
  • domain assumption Quark matter at T~1 GeV is described by a highly degenerate, ultrarelativistic ideal gas of u,d,s quarks with the EOS of Eqs. (1)-(2).
    The entire stellar-structure calculation reduces to this EOS; the T^2 coefficient is inconsistent with a standard 3-flavor Fermi gas.
  • domain assumption The universe at T=1 GeV is radiation-dominated with g*=61.75 and the standard baryon-to-photon ratio.
    Used for the boundary pressure and horizon mass; Sec. III.D's average baryon content is inconsistent with this input.
  • ad hoc to paper PQ domain walls or superhorizon bubbles can concentrate baryons sufficiently to seed PQP formation.
    Section III.D is explicitly qualitative; no hydrodynamics or sweep efficiency is computed.
  • domain assumption After cooling, PQPs remain gravitationally bound as mini neutron stars or strange-matter nuggets.
    Sec. III.B; if hadronization ejects baryons or the strange-matter ground state is wrong, the dark-matter abundance changes.
  • standard math TOV and Friedmann equations are the correct gravitational framework.
    Used in Eqs. (9)-(14).
invented entities (1)
  • Primordial quark pellet (PQP) independent evidence
    purpose: Dark matter candidate composed of hot or cold quark matter in compact objects.
    Not a new fundamental field, but a new macroscopic configuration with falsifiable handles: predicted mass range, microlensing limits, and possible LISA gravitational-wave signature.

pith-pipeline@v1.3.0-alltime-deepseek · 8236 in / 24594 out tokens · 239103 ms · 2026-08-02T07:10:51.871129+00:00 · methodology

0 comments
read the original abstract

We show that primordial quark pellets (PQP), ultra dense quark matter mini-stars formed at around 1 TeV, naturally arise in a radiation dominated universe if rare baryon overdensities are produced by collapsing Peccei-Quinn domain walls or similar superhorizon structures. Solving the Tolman-Oppenheimer-Volkoff equation with a hot quark equation of state, we find stable solutions with a maximum mass of $10^{-2}$ solar masses and radii of approximately 100 meters, although the formation mechanism favours much smaller objects. Once formed, PQPs cool and evolve intro mini neutron stars or stable strange matter nuggets, depending on the QCD ground state. Formation probabilities of $10^{-9}$ to $10^{-4}$ per horizon volume could suffice to reproduce the present dark matter density without altering Big Bang Nucleosynthesis or requiring entropy dilution. PQPs completely evade microlensing constraints if their mass is below $10^{-7}$ solar masses, a range that could easily accommodate most of the dark matter. PQPs thus potentially constitute a conservative, Standard-Model based, and observationally testable solution to the dark matter puzzle.

Figures

Figures reproduced from arXiv: 2607.10672 by Domenec Espriu.

Figure 1
Figure 1. Figure 1: FIG. 1: Mass versus central density resulting from the integration of the relativistic TOV equation [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Mass versus central density resulting from the integration of the relativistic TOV equation [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗

Review history (2 revisions) →

discussion (0)

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

Reference graph

Works this paper leans on

18 extracted references · 3 linked inside Pith

  1. [1]

    Cosmic separation of phases

    H.Terazawa, talk presented at the Symposium on Present Status and Future Prospect of Elementary-Particle Physics and Nuclear Physics, INS, Univ. of Tokyo, March 22–24, 1979; E. Witten, “Cosmic separation of phases",Physical Review D,30, 272 (1984); E. Farhi and R.L. Jaffe, “Strange matter,"Physical Review D,30, 2379 (1984); Y. Bai, A.J. Long and S.Lu, “Da...

    1979

  2. [2]

    A history of dark matter

    G. Bertone and D. Hooper, “A history of dark matter”,Review of Modern Physics,90, 45002 (2018)

    2018

  3. [3]

    Primordial neutron stars

    G. Krnjaic, D. Rocha and H. Xiao, “Primordial neutron stars”, https://arxiv.org/abs/2604.08651

    Pith/arXiv arXiv

  4. [4]

    Big-Bang Nucleosynthesis after Planck,

    B. D. Fieldset al., “Big-Bang Nucleosynthesis after Planck,”JCAP03, 010 (2020)

    2020

  5. [5]

    CP Conservation in the Presence of Pseudoparticles

    R. D. Peccei and H.R.Quinn, “CP Conservation in the Presence of Pseudoparticles",Phys. Rev. Lett.,38, 1440 (1977)

    1977

  6. [6]

    Primordial Black Holes from Cosmic Domain Walls

    L. Liu, Z.-K. Guo and R.-G. Cai, “Primordial Black Holes from Cosmic Domain Walls", arXiv:1901.07072 (2019)

    Pith/arXiv arXiv 1901

  7. [7]

    More axions from diluted domain walls

    K. Harigaya, K. Kawasaki, M. Mukaida and T.T. Yanagida, “More axions from diluted domain walls",Phys. Rev. D,94, 063506 (2016); M. Gorghetto, E. Hardy and G. Villadoro, “Early vs late string networks from a minimal QCD Axion",JHEP,06, 151 (2018)

    2016

  8. [8]

    Cosmicstringsandothertopologicaldefects

    A.VilenkinandE.P.S.Shellard, “Cosmicstringsandothertopologicaldefects”, CambridgeUni- 14 versity Press (1994); I. Baldes, Y. Gouttenoire and F. Sala, “Bubbletrons", arXiv:2306.15555 (2023)

    Pith/arXiv arXiv 1994

  9. [9]

    Formation of galaxies and clusters of galaxies by self-similar gravitational condensation

    W. H. Press and P. Schechter, “Formation of galaxies and clusters of galaxies by self-similar gravitational condensation",Astrophys. J.,187, 425 (1974)

    1974

  10. [10]

    Microlensing constraints on primordial black holes with Subaru/Hyper Suprime-Cam

    H. Niikuraet al.(HSC Collaboration), “Microlensing constraints on primordial black holes with Subaru/Hyper Suprime-Cam”,Nature Astron.3, 524–534 (2019)

    2019

  11. [11]

    Constraints on Earth-mass primordial black holes from OGLE 5-year microlensing events,

    H. Niikuraet al.(OGLE Collaboration), “Constraints on Earth-mass primordial black holes from OGLE 5-year microlensing events,”Phys. Rev. D99, 083503 (2019)

    2019

  12. [12]

    No massive black holes in the Milky Way halo

    P. Mroz et al.,“No massive black holes in the Milky Way halo”,Nature,632, 749 (2024); “Limits on Planetary-mass Primordial Black Holes from the OGLE High-cadence Survey of the Magellanic Clouds”,The Astrophysical Journal Letters,976, L19 (2024)

    2024

  13. [13]

    Baryonic Dark Matter,

    B. Carr, “Baryonic Dark Matter,”Annual Review Astronomy and Astrophysics,32, 531 (1994)

    1994

  14. [14]

    Limits on the MACHO content of the Galactic Halo from the EROS-2 Survey of the Magellanic Clouds

    P. Tisserand et al. (EROS-2 Collaboration), “Limits on the MACHO content of the Galactic Halo from the EROS-2 Survey of the Magellanic Clouds”Astronomy & Astrophysics,469, 387 (2007)

    2007

  15. [15]

    Can Finite Source Sizes Be Resolved in Gravitational Microlensing Events?

    H.J. Witt and S. Mao, “Can Finite Source Sizes Be Resolved in Gravitational Microlensing Events?”The Astrophysical Journal,430, 505 (1994)

    1994

  16. [16]

    Finite Source Effects in Microlensing Events

    A. Gould and C. Gaucherel, “Finite Source Effects in Microlensing Events”Astrophysical Jour- nal,477, 580 (1997)

    1997

  17. [17]

    Science with the space-based interferometer eLISA. II: Gravitational waves from cosmological phase transitions,

    C. Capriniet al., “Science with the space-based interferometer eLISA. II: Gravitational waves from cosmological phase transitions,”JCAP04, 001 (2016)

    2016

  18. [18]

    Detecting gravitational waves from cosmological phase transitions with LISA: an update,

    C. Capriniet al., “Detecting gravitational waves from cosmological phase transitions with LISA: an update,”JCAP03, 024 (2020). 15

    2020