pith. sign in

arxiv: 2605.16465 · v1 · pith:NAPUBYBOnew · submitted 2026-05-15 · ✦ hep-ph · astro-ph.CO· physics.comp-ph

Freeze-in Warm Dark Matter via Dimension-6 Operators in 3-3-1 Models

Fuwei Liu , Fujun Liu This is my paper

Pith reviewed 2026-05-20 17:20 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.COphysics.comp-ph
keywords freeze-in dark matterwarm dark matterdimension-six operators3-3-1 modelsZ13 symmetryPeccei-Quinn symmetrysterile singletrelic abundance
0
0 comments X

The pith

An exact Z13 symmetry suppresses all low-dimensional portals so dark matter interacts with the Standard Model only through a dimension-six operator.

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

The paper shows how to resolve the fine-tuning issue in freeze-in dark matter models by placing a sterile singlet inside a 3-3-1 electroweak extension. An exact Z13 discrete gauge symmetry is imposed to eliminate all operators of dimension lower than six. The effective coupling then arises automatically from the large separation between the electroweak scale and the high scale of Peccei-Quinn symmetry breaking. Boltzmann equation integration confirms that this setup produces a keV-mass warm dark matter particle with the correct relic density of 0.12. The candidate also stays compatible with Lyman-alpha bounds and can help explain galactic structure puzzles.

Core claim

By imposing an exact Z13 discrete gauge symmetry in a 3-3-1 model, all low-dimensional portals are formally suppressed so that the dark sector communicates with the Standard Model exclusively through a dimension-six operator. This structure lets the small coupling required for freeze-in production emerge naturally from the hierarchy between the electroweak scale and the ultra-high Peccei-Quinn symmetry breaking scale. Numerical integration of the Boltzmann equations shows that the sterile singlet can be produced to match the observed relic abundance Ω_S h² = 0.12.

What carries the argument

The exact Z13 discrete gauge symmetry that forbids lower-dimensional interactions, forcing the dark matter to couple to the Standard Model only via a dimension-six operator.

If this is right

  • The sterile singlet is produced through the infrared freeze-in mechanism to achieve the observed relic density.
  • The keV-scale warm dark matter remains consistent with Lyman-alpha forest constraints.
  • It provides a viable solution to galactic-scale issues such as the cusp-core and missing satellites problems.
  • The framework unifies dark matter production with the strong CP problem solution without fine-tuning parameters.

Where Pith is reading between the lines

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

  • The same Z13 protection could be applied in other gauge extensions to generate naturally small couplings for hidden sector particles.
  • Future cosmological observations of small-scale structure might further test the keV warm dark matter prediction.
  • Collider searches in 3-3-1 models could look for signatures of the extended gauge sector that accompanies this dark matter setup.

Load-bearing premise

The Z13 symmetry stays exact at all scales and does not produce any extra lower-dimensional operators or mixing that would ruin the dimension-six suppression.

What would settle it

A calculation or observation demonstrating that the Z13 symmetry is violated at some scale, allowing a lower-dimensional portal that would either overproduce the dark matter or violate current constraints.

read the original abstract

We propose a natural resolution to the fine-tuning problem inherent in the freeze-in dark matter paradigm by embedding a sterile singlet within a 3-3-1 electroweak extension. By imposing an exact $Z_{13}$ discrete gauge symmetry, we formally suppress all low-dimensional portals to ensure that the dark sector communicates with the Standard Model (SM) exclusively through a dimension-six operator. This theoretical structure allows the extraordinarily small coupling required for dark matter production to emerge naturally from the profound hierarchy between the electroweak scale and the ultra-high Peccei-Quinn symmetry breaking scale. Detailed numerical integration of the Boltzmann equations demonstrates that the sterile singlet can be produced via the infrared freeze-in mechanism to match the observed relic abundance of $\Omega_S h^2 = 0.12$. The resulting keV-scale warm dark matter candidate remains consistent with stringent Lyman-alpha forest constraints while offering a viable solution to galactic-scale discrepancies such as the cusp-core and missing satellites problems. Ultimately, this framework provides a self-consistent unification of dark matter genesis and the strong CP solution that is completely independent of ad hoc parameter adjustments.

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

2 major / 2 minor

Summary. The paper proposes embedding a sterile singlet S as warm dark matter in a 3-3-1 electroweak extension, protected by an exact Z13 discrete gauge symmetry. This symmetry is asserted to forbid all portals of dimension 5 or lower, so that the leading interaction is a dimension-6 operator suppressed by the high Peccei-Quinn scale; numerical integration of the Boltzmann equations is stated to reproduce the observed relic density Ω_S h² = 0.12 for a keV-scale candidate consistent with Lyman-alpha constraints.

Significance. If the Z13 charge assignments are shown to eliminate all lower-dimensional operators and the numerical results prove robust under explicit parameter choices, the work supplies a concrete mechanism for natural freeze-in production tied to the electroweak–PQ hierarchy while simultaneously addressing the strong-CP problem inside a 3-3-1 framework.

major comments (2)
  1. [Section describing the Z13 symmetry and field content] The load-bearing naturalness claim requires that the chosen Z13 charges forbid every gauge-invariant operator of dimension ≤5 involving S. The manuscript must supply an exhaustive list of all such operators (including those constructed with the additional scalars that break the 3-3-1 symmetry) together with the explicit Z13 charge assignments for every field; without this enumeration the assertion that the effective coupling emerges automatically from the scale hierarchy remains unverified.
  2. [Section on Boltzmann-equation integration and relic-density calculation] The central numerical result—that freeze-in via the dimension-6 operator yields Ω_S h² = 0.12—depends on the precise form of the operator, the integration limits, and the treatment of the 3-3-1 gauge-boson and scalar contributions. The manuscript should present the Boltzmann equation(s) used, the numerical values of the effective coupling and PQ scale, and any convergence or error estimates.
minor comments (2)
  1. [Introduction and model-definition section] Clarify the precise relation between the PQ breaking scale and the 3-3-1 symmetry-breaking scales; the abstract refers to an “ultra-high” PQ scale but does not quantify how it is decoupled from the TeV-scale 3-3-1 vevs.
  2. [Symmetry section] The phrase “discrete gauge symmetry” should be accompanied by a brief statement on anomaly cancellation or consistency conditions for the Z13 assignment.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address the major comments point by point below and will revise the paper to incorporate the requested clarifications.

read point-by-point responses

  1. Referee: The load-bearing naturalness claim requires that the chosen Z13 charges forbid every gauge-invariant operator of dimension ≤5 involving S. The manuscript must supply an exhaustive list of all such operators (including those constructed with the additional scalars that break the 3-3-1 symmetry) together with the explicit Z13 charge assignments for every field; without this enumeration the assertion that the effective coupling emerges automatically from the scale hierarchy remains unverified.

    Authors: We agree that an explicit enumeration of the forbidden operators would make the naturalness argument more transparent and verifiable. In the revised manuscript we will add a new subsection that lists all gauge-invariant operators of dimension 5 or lower involving the singlet S (including those built with the additional 3-3-1 scalars), together with a table of the Z_{13} charge assignments for every field. This will allow direct verification that the Z_{13} symmetry indeed eliminates all lower-dimensional portals, leaving only the dimension-6 operator suppressed by the Peccei-Quinn scale. revision: yes

  2. Referee: The central numerical result—that freeze-in via the dimension-6 operator yields Ω_S h² = 0.12—depends on the precise form of the operator, the integration limits, and the treatment of the 3-3-1 gauge-boson and scalar contributions. The manuscript should present the Boltzmann equation(s) used, the numerical values of the effective coupling and PQ scale, and any convergence or error estimates.

    Authors: We accept that additional technical details on the numerical calculation are needed for reproducibility. In the revised version we will explicitly display the Boltzmann equation(s) employed for the infrared freeze-in yield, quote the numerical values chosen for the effective dimension-6 coupling and the Peccei-Quinn scale, specify the integration temperature range, and include a brief discussion of numerical convergence and estimated uncertainties in the final relic density. These additions will substantiate the reported result Ω_S h² = 0.12. revision: yes

Circularity Check

1 steps flagged

Relic density matching via Boltzmann integration reduces to fitting the PQ scale or effective coupling to the observed value

specific steps
  1. fitted input called prediction [Abstract]
    "Detailed numerical integration of the Boltzmann equations demonstrates that the sterile singlet can be produced via the infrared freeze-in mechanism to match the observed relic abundance of Ω_S h² = 0.12."

    The integration is presented as demonstrating successful production, but reproducing the specific numerical value Ω_S h² = 0.12 requires choosing the dimension-6 coefficient (set by the PQ scale in the hierarchy) to reproduce the observed density; this converts the claimed natural emergence into a parameter adjustment calibrated to the target rather than a first-principles output independent of the observed abundance.

full rationale

The paper's derivation imposes Z13 to enforce dim-6-only communication, then invokes the EW/PQ hierarchy for natural smallness, followed by numerical Boltzmann integration that reproduces Ω_S h² = 0.12. This matching step requires selecting the scale or coefficient to hit the target abundance, rendering the 'demonstration' a fit rather than an independent prediction. The symmetry choice itself is a standard construction and does not reduce by definition to the target result. No self-citation chains or other enumerated circular patterns appear in the provided text. The central naturalness claim retains independent content from the discrete symmetry and 3-3-1 embedding, so the circularity is partial and localized to the relic-density step.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The central claim rests on the 3-3-1 gauge extension and the exact Z13 symmetry as domain assumptions, plus the sterile singlet as a new postulated particle; the small effective coupling is presented as emerging from scale hierarchy but still requires numerical matching to the observed density.

free parameters (1)
  • effective coupling or PQ breaking scale
    The hierarchy is invoked to generate the small coupling, yet the numerical integration is used to match the specific relic density value of 0.12, implying adjustment of at least one scale parameter.
axioms (1)
  • domain assumption Exact Z13 discrete gauge symmetry suppresses all low-dimensional portals
    Invoked in the abstract to ensure communication occurs exclusively through the dimension-six operator.
invented entities (1)
  • sterile singlet S no independent evidence
    purpose: keV-scale warm dark matter candidate
    New particle introduced in the model with no independent evidence supplied beyond the proposal itself.

pith-pipeline@v0.9.0 · 5728 in / 1643 out tokens · 167915 ms · 2026-05-20T17:20:54.486104+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

27 extracted references · 27 canonical work pages · 1 internal anchor

  1. [1]

    Aghanim et al.,Planck 2018 results

    N. Aghanim et al.,Planck 2018 results. vi. cosmological parameters,Astron. Astrophys641 (2020) A6

    work page 2018

  2. [2]

    Bullock and M

    J.S. Bullock and M. Boylan-Kolchin,Small-scale challenges to theλcdm paradigm,Annual Review of Astronomy and Astrophysics55(2017) 343

    work page 2017

  3. [3]

    M. Viel, J. Lesgourgues, M.G. Haehnelt, S. Matarrese and A. Riotto,Constraining warm dark matter candidates including sterile neutrinos and light gravitinos with wmap and the lyman-α forest,Physical Review D—Particles, Fields, Gravitation, and Cosmology71(2005) 063534

    work page 2005

  4. [4]

    L.J. Hall, K. Jedamzik, J. March-Russell and S.M. West,Freeze-in production of fimp dark matter,Journal of High Energy Physics2010(2010) 1

    work page 2010

  5. [5]

    Bernal, M

    N. Bernal, M. Heikinheimo, T. Tenkanen, K. Tuominen and V. Vaskonen,The dawn of fimp dark matter: a review of models and constraints,International Journal of Modern Physics A 32(2017) 1730023

    work page 2017

  6. [6]

    Frampton,Chiral dilepton model and the flavor question,Physical Review Letters69 (1992) 2889

    P.H. Frampton,Chiral dilepton model and the flavor question,Physical Review Letters69 (1992) 2889

    work page 1992

  7. [7]

    Liu and F

    F. Liu and F. Liu,Gravitational robustness of automatic peccei-quinn symmetry in 3-3-1 models,Physics of the Dark Universe(2026) 102340

    work page 2026

  8. [8]

    Dong, T.D

    P.V. Dong, T.D. Tham and H.T. Hung,3-3-1-1 model for dark matter,Phys. Rev. D87(2013) 115003. – 12 –

    work page 2013

  9. [9]

    Ib´ a˜ nez and G.G

    L.E. Ib´ a˜ nez and G.G. Ross,Discrete gauge symmetries and the origin of baryon and lepton number conservation in supersymmetric versions of the standard model,Nuclear Physics B368 (1992) 3

    work page 1992

  10. [10]

    Raffelt,Astrophysical axion bounds, inAxions: Theory, Cosmology, and Experimental Searches, pp

    G.G. Raffelt,Astrophysical axion bounds, inAxions: Theory, Cosmology, and Experimental Searches, pp. 51–71, Springer (2008)

    work page 2008

  11. [11]

    Bode, J.P

    P. Bode, J.P. Ostriker and N. Turok,Halo formation in warm dark matter models,The Astrophysical Journal556(2001) 93

    work page 2001

  12. [12]

    Bernal, F

    N. Bernal, F. Elahi, C. Maldonado and J. Unwin,Ultraviolet freeze-in and non-standard cosmologies,Journal of Cosmology and Astroparticle Physics2019(2019) 026

    work page 2019

  13. [13]

    Aebischer, W

    J. Aebischer, W. Altmannshofer, E.E. Jenkins and A.V. Manohar,Dark matter effective field theory and an application to vector dark matter,Journal of High Energy Physics2022(2022) 86

    work page 2022

  14. [14]

    Ibanez and G.G

    L.E. Ibanez and G.G. Ross,Discrete gauge symmetries and the origin of baryon and lepton number conservation in supersymmetric versions of the standard model,Nuclear Physics B368 (1992) 3

    work page 1992

  15. [15]

    Higgs-field Portal into Hidden Sectors

    B. Patt and F. Wilczek,Higgs-field portal into hidden sectors,arXiv preprint hep-ph/0605188 (2006)

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  16. [16]

    Lebedev, H.M

    O. Lebedev, H.M. Lee and Y. Mambrini,Vector higgs portal dark matter and the invisible higgs,Physics Letters B707(2012) 570

    work page 2012

  17. [17]

    Frigerio, T

    M. Frigerio, T. Hambye and E. Masso,Sub-gev dark matter as pseudo-nambu-goldstone bosons from the seesaw scale,Physical Review X1(2011) 021026

    work page 2011

  18. [18]

    Saikawa and S

    K. Saikawa and S. Shirai,Primordial gravitational waves, precisely: The role of thermodynamics in the standard model,Journal of Cosmology and Astroparticle Physics2018(2018) 035

    work page 2018

  19. [19]

    Gondolo and G

    P. Gondolo and G. Gelmini,Cosmic abundances of stable particles: Improved analysis,Nuclear Physics B360(1991) 145

    work page 1991

  20. [20]

    Irˇ siˇ c, M

    V. Irˇ siˇ c, M. Viel, M.G. Haehnelt, J.S. Bolton, S. Cristiani, G.D. Becker et al.,New constraints on the free-streaming of warm dark matter from intermediate and small scale lyman-αforest data,Physical Review D96(2017) 023522

    work page 2017

  21. [21]

    Pal and L

    P.B. Pal and L. Wolfenstein,Radiative decays of massive neutrinos,Physical Review D25 (1982) 766

    work page 1982

  22. [22]

    Roach, K.C

    B.M. Roach, K.C. Ng, K. Perez, J.F. Beacom, S. Horiuchi, R. Krivonos et al.,Nustar tests of sterile-neutrino dark matter: new galactic bulge observations and combined impact,Physical Review D101(2020) 103011

    work page 2020

  23. [23]

    Boyarsky, A

    A. Boyarsky, A. Neronov, O. Ruchayskiy and M. Shaposhnikov,Constraints on sterile neutrinos as dark matter candidates from the diffuse x-ray background,Monthly Notices of the Royal Astronomical Society370(2006) 213

    work page 2006

  24. [24]

    Fields, K.A

    B.D. Fields, K.A. Olive, T.-H. Yeh and C. Young,Big-bang nucleosynthesis after planck, Journal of Cosmology and Astroparticle Physics2020(2020) 010

    work page 2020

  25. [25]

    Gildener,Gauge-symmetry hierarchies,Physical Review D14(1976) 1667

    E. Gildener,Gauge-symmetry hierarchies,Physical Review D14(1976) 1667

    work page 1976

  26. [26]

    G. Aad, B. Abbott, D.C. Abbott, A.A. Abud, K. Abeling, D. Abhayasinghe et al.,Search for high-mass dilepton resonances using 139 fb- 1 of pp collision data collected at s= 13 tev with the atlas detector,Physics Letters B796(2019) 68

    work page 2019

  27. [27]

    Alves, L

    A. Alves, L. Duarte, S. Kovalenko, Y. Oviedo-Torres, F. Queiroz and Y. Villamizar, Constraining 3-3-1 models at the lhc and future hadron colliders,Physical Review D106 (2022) 055027. – 13 –

    work page 2022