arxiv: 2604.06038 · v1 · submitted 2026-04-07 · 🌌 astro-ph.CO · astro-ph.GA

Recognition: 2 theorem links

· Lean Theorem

Lyman-α Forest Signatures of Mixed Fuzzy and Cold Dark Matter

Yourong Frank Wang

Authors on Pith no claims yet

Pith reviewed 2026-05-10 18:14 UTC · model grok-4.3

classification 🌌 astro-ph.CO astro-ph.GA

keywords Lyman-alpha forestfuzzy dark mattermixed dark mattervelocity power spectrumflux power spectrumSchrödinger-Poisson simulationwave mechanicsdark matter dynamics

0 comments

The pith

Mixed fuzzy and cold dark matter produce 10 percent differences in Lyman-alpha flux power spectra through velocity suppression, even when matter power spectra are nearly identical.

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

The paper compares two simulations of mixed fuzzy and cold dark matter with identical initial conditions: one using full wave-mechanical Schrödinger-Poisson evolution and the other using standard N-body particles. Despite near-degeneracy in the nonlinear matter power spectrum, the Lyman-alpha flux power spectra differ at the 10 percent level on intermediate scales. This discrepancy stems from a strong suppression of small-scale velocity power in the wave evolution, which affects the gas distribution in ways the matter power spectrum alone does not capture. A sympathetic reader would care because the result indicates that Lyman-alpha forest observations can reveal dynamical details of dark matter models that power-spectrum probes miss, offering a path to break degeneracies between different dark matter compositions.

Core claim

Evolving the dark matter distribution from redshift 120 to 2 for an axion mass of 0.01 and 10 percent fuzzy dark matter fraction, the Schrödinger-Poisson simulation produces strong suppression of small-scale velocity power compared to the particle-only N-body run. As a result, the corresponding Lyman-alpha flux power spectra computed via the Fluctuating Gunn-Peterson Approximation differ at the 10 percent level on intermediate scales even though the matter power spectra are nearly degenerate. This demonstrates that the flux statistics depend sensitively on the dynamical evolution of the velocity field, and that wave-mechanical effects in fuzzy dark matter leave distinct kinematic imprints in

What carries the argument

Comparison of Schrödinger-Poisson wave evolution versus N-body particle evolution with identical initial conditions, isolating the effect of wave dynamics on velocity power suppression.

If this is right

Lyman-alpha flux power spectra differ at the 10 percent level on intermediate scales despite near-degeneracy in the nonlinear matter power spectrum.
Flux statistics depend sensitively on the dynamical evolution of the velocity field rather than the matter power spectrum alone.
Wave-mechanical effects in fuzzy dark matter produce distinct kinematic imprints in Lyman-alpha observables beyond those from initial-condition suppression.
Mixed dark matter models can break degeneracies present in standard structure-based probes.

Where Pith is reading between the lines

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

Full hydrodynamical simulations are needed to check whether these velocity-induced flux differences survive realistic gas physics.
Existing or future Lyman-alpha forest datasets could be used to place new constraints on the fuzzy dark matter fraction.
Similar velocity-dependent distinctions might appear in other observables sensitive to small-scale motions, such as redshift-space distortions in galaxy surveys.

Load-bearing premise

The Fluctuating Gunn-Peterson Approximation applied to the simulation outputs accurately captures the Lyman-alpha flux statistics without full hydrodynamical gas treatment.

What would settle it

A full hydrodynamical simulation of the same mixed fuzzy and cold dark matter initial conditions that shows whether the 10 percent flux power spectrum difference between wave and particle evolution persists or disappears.

Figures

Figures reproduced from arXiv: 2604.06038 by Yourong Frank Wang.

**Figure 1.** Figure 1: Initial density transfer functions for the FDM and CDM components for 𝑚22 = 0.01 and 𝑓A = 0.1. The curves are normalized against the CDM transfer function 𝑇CDM (𝑘). The baryonic transfer function is not explicitly used, but the BAO signatures are visible in the total matter and particles curves. The Fundamental and Nyquist wavenumbers of simulation box are marked on the 𝑥-axis. The transfer functions for… view at source ↗

**Figure 2.** Figure 2: A slice view of the initial overdensities generated by MUSIC, showing the axion’s smoothing effects. 𝑁 = 10243 . The three panels share the same color scale. scales, setting the physical baryon density to 𝜌b (x) ≡ 𝑓b 𝜌m(x) (1 + 𝑧) 3 , 𝑓b = Ωb/Ωm. (5) We apply a phenomenological baryonic pressure smoothing to the density and velocity fields using a Gaussian filter with width 𝜆F ≈ 0.2 Mpc ( 1 1 + 𝑧 )1/2 , (6… view at source ↗

**Figure 3.** Figure 3: Dimensionless matter power spectra at 𝑧 = 120, 30, 4, and 2 (top panel), and the ratios of mixed dark matter (MDM) simulations to the ΛCDM baseline at 𝑧 = 120 vs. 𝑧 = 2 (bottom panel). Solid lines correspond to the full Schrödinger–Poisson evolution (MDM Full), and dotted lines show the particle-only approximation (MDM N-body). At 𝑧 = 2, both MDM models have less suppression at small scales, as the CDM com… view at source ↗

**Figure 4.** Figure 4: Projected density fields at 𝑧 = 2 within the 30 Mpc/ℎ simulation box. The left panels show a zoomed-in subregion (width ≈ 5 cMpc around a halo) highlighting the distinct morphologies of the particle (CDM) and FDM components. The right panel displays the total projected density across the full simulation volume, with the location of the zoomed subregion indicated on it. While the CDM follows the standard co… view at source ↗

**Figure 7.** Figure 7: Line-of-sight velocity power spectra at 𝑧 = 120 and 𝑧 = 2, shown as ratios to the ΛCDM baseline. At 𝑧 = 120, both MDM treatments show comparable mild suppression inherited from the initial conditions. By 𝑧 = 2, the N-body approximation retains ∼ 70 per cent of the ΛCDM velocity power on small scales, while the full Schrödinger–Poisson treatment drives it to near zero —reflecting the hard kinematic cutoff i… view at source ↗

**Figure 5.** Figure 5: Ratio of the Lyman-𝛼 1D flux power spectrum in the two MDM simulations to the corresponding ΛCDM result at 𝑧 = 4 and 𝑧 = 2. Mean fluxes are calibrated to the same fiducial value at each redshift. Shaded bands show the standard deviation across the sample slabs. The common lower scale gives 𝑘 in the unit ℎ Mpc−1 , while the upper scales show the corresponding speed ranges in s/km; the conversion is redshift… view at source ↗

**Figure 6.** Figure 6: The Lyman-𝛼 transmission power spectrum compared with the corresponding LCDM result as the IGM temperature at mean density, T0 changes between 6 000, 10 000, and 15 000 K. For each value of T0, the mean flux is independently calibrated so that the LCDM box matches the fiducial mean opacity. 10 0 10 1 k [h/Mpc] 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 Pv||(k, z)/Pv||, LCDM(k, z) kJ(z = 120) kJ(z = 2) z = 2 z = … view at source ↗

read the original abstract

We investigate Lyman-alpha forest flux statistics in mixed fuzzy dark matter (FDM) and cold dark matter (CDM) cosmologies using the Fluctuating Gunn-Peterson Approximation (FGPA) applied to hybrid Schr\"odinger-Poisson and N-body simulations. We evolve the dark matter distribution from z = 120 to z = 2 for an axion mass ( m_22 = 0.01) and FDM fraction (f_A = 0.1), and compare two realizations with identical initial conditions: one evolved with a particle-only approximation and one with full wave--mechanical dynamics. We find that, despite near-degeneracy in the nonlinear matter power spectrum, the corresponding Ly {\alpha} flux power spectra differ at the 10 percent level on intermediate scales. This discrepancy arises from a strong suppression of small-scale velocity power in the Schr\"odinger--Poisson evolution, which is not captured by N-body treatments with matched initial transfer functions. As a result, the flux statistics cannot be fully characterized by the matter power spectrum alone, but depend sensitively on the dynamical evolution of the velocity field. These results demonstrate that wave-mechanical effects in FDM leave distinct kinematic imprints in Ly {\alpha} observables beyond those associated with initial-condition suppression. While our analysis is based on an idealized FGPA framework, it isolates a mechanism by which mixed dark matter models can break degeneracies present in standard structure-based probes, motivating further investigation with full hydrodynamical simulations.

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.

Desk Editor's Note private letter to a colleague

The paper shows that wave dynamics in mixed FDM-CDM can shift Lyα flux power by 10% through velocity suppression even when nonlinear matter power spectra match, but this rests on FGPA without hydro checks.

read the letter

The main takeaway is that velocity fields from Schrödinger-Poisson evolution leave a measurable imprint on Lyman-alpha flux statistics that the matter power spectrum alone does not capture. With identical initial conditions, the hybrid run and the N-body run produce near-degenerate nonlinear matter spectra at m_22=0.01 and f_A=0.1, yet the flux power differs at the 10 percent level on intermediate scales because small-scale velocity power is suppressed in the wave treatment. That separation is the concrete new piece here, and the direct comparison of the two realizations makes the point cleanly without relying on fitted parameters or circular definitions. The work is useful for anyone thinking about how to break degeneracies in future Lyα constraints on fuzzy dark matter fractions. It isolates a kinematic channel that standard transfer-function arguments miss. The limitation is straightforward: the flux statistics come from the Fluctuating Gunn-Peterson Approximation applied straight to the dark-matter fields. That approximation omits thermal broadening, Jeans smoothing, and the fact that baryons do not trace dark-matter velocities perfectly on the relevant scales. Those effects are known to shape the Lyα power spectrum and could easily reduce or erase the reported 10 percent difference. No convergence tests, error bars, or full hydrodynamical runs are shown, so the result stays conditional on the validity of the idealized mapping. The paper is for people already working on small-scale structure probes and mixed dark-matter models. A reader who wants to see whether the velocity imprint survives realistic gas physics will get value from it. It is not yet ready to change how surveys interpret data, but the mechanism is worth checking. I would send it to peer review so referees can press on the approximation and ask for hydro validation or at least a clearer statement of the uncertainty range.

Referee Report

2 major / 1 minor

Summary. The manuscript investigates Lyman-α forest flux statistics in mixed fuzzy dark matter (FDM) and cold dark matter (CDM) models. It applies the Fluctuating Gunn-Peterson Approximation (FGPA) to hybrid Schrödinger-Poisson and N-body simulations evolved from z=120 to z=2 with parameters m_22=0.01 and f_A=0.1, using two realizations with identical initial conditions. Despite near-degenerate nonlinear matter power spectra, the Lyα flux power spectra differ at the 10% level on intermediate scales; this is attributed to strong suppression of small-scale velocity power in the wave-mechanical evolution, implying that flux statistics depend sensitively on velocity-field dynamics beyond the matter power spectrum alone.

Significance. If the reported difference survives more complete modeling, the work would show that wave-mechanical effects in FDM produce distinct kinematic imprints in Lyα observables not captured by initial-condition suppression or matter power spectrum comparisons. A clear strength is the controlled comparison of Schrödinger-Poisson versus N-body evolution with matched initial conditions, which isolates the dynamical velocity contribution without confounding transfer-function differences.

major comments (2)

The central claim (abstract) that 'the flux statistics cannot be fully characterized by the matter power spectrum alone' rests on a 10% flux-power discrepancy arising from velocity suppression. However, this discrepancy is reported at a single unspecified scale range without error bars, convergence tests, or quantification of the exact k-interval and ratio, undermining assessment of whether the effect is robust or numerical.
The results rely on direct application of the idealized FGPA to the DM density and velocity fields (abstract and simulation description). This approximation omits thermal broadening (~10 km/s at z=2), Jeans smoothing, and baryon-DM velocity decoupling, all of which are known to dominate Lyα power-spectrum shape on the intermediate scales of interest and could reduce or erase the reported 10% difference.

minor comments (1)

The abstract states 'near-degeneracy' in the nonlinear matter power spectrum but does not quantify the level of agreement (e.g., percent difference as a function of k); adding this detail would strengthen the contrast with the flux-power result.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the careful reading of our manuscript and for highlighting both the strengths of our controlled comparison and the areas needing clarification. We address each major comment below and indicate the revisions we will make.

read point-by-point responses

Referee: The central claim (abstract) that 'the flux statistics cannot be fully characterized by the matter power spectrum alone' rests on a 10% flux-power discrepancy arising from velocity suppression. However, this discrepancy is reported at a single unspecified scale range without error bars, convergence tests, or quantification of the exact k-interval and ratio, undermining assessment of whether the effect is robust or numerical.

Authors: We agree that the manuscript would benefit from greater quantitative precision. In the revised version we will explicitly report the wavenumber interval over which the ~10% difference appears, provide the scale-dependent ratio of the flux power spectra, and include error bars derived from the two realizations. We will also add a brief discussion of numerical convergence based on the simulation resolution and box size used. These additions will allow readers to assess the robustness of the reported difference directly from the data. revision: yes
Referee: The results rely on direct application of the idealized FGPA to the DM density and velocity fields (abstract and simulation description). This approximation omits thermal broadening (~10 km/s at z=2), Jeans smoothing, and baryon-DM velocity decoupling, all of which are known to dominate Lyα power-spectrum shape on the intermediate scales of interest and could reduce or erase the reported 10% difference.

Authors: This is a substantive limitation of the current analysis. The FGPA was chosen precisely because it isolates the kinematic effect of the dark-matter velocity field while keeping all other aspects identical between the Schrödinger-Poisson and N-body runs. We will expand the caveats section to discuss the expected magnitude of thermal broadening, Jeans smoothing, and relative baryon-DM velocities on the relevant scales, citing the relevant literature. We will also state more clearly that the 10% difference we report is a lower bound on the possible signature and that only full hydrodynamical simulations can determine whether the kinematic imprint survives these additional physics. Such simulations lie beyond the scope of the present work. revision: partial

standing simulated objections not resolved

Whether the reported 10% flux-power difference survives the inclusion of thermal broadening, Jeans smoothing, and baryon-DM velocity decoupling cannot be answered without performing full hydrodynamical simulations, which are outside the scope of this study.

Circularity Check

0 steps flagged

No significant circularity: direct simulation comparison with independent evolution methods

full rationale

The paper's central result follows from evolving identical initial conditions under two distinct dynamical prescriptions (Schrödinger-Poisson wave mechanics versus N-body particle evolution), then applying the standard FGPA mapping to the resulting density and velocity fields. The reported 10% difference in flux power is an output of those separate integrations, not a redefinition or fit of the input power spectra. No step reduces by construction to the matter power spectrum alone, and FGPA is invoked as an external approximation rather than derived from the present data. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The claim depends on the validity of the FGPA mapping from dark matter density and velocity to flux statistics, plus the assumption that the chosen simulation parameters and hybrid evolution capture the relevant physics at z=2.

free parameters (2)

m_22 = 0.01
Axion mass parameter set to 0.01 for the simulation runs
f_A = 0.1
Fraction of fuzzy dark matter set to 0.1

axioms (1)

domain assumption FGPA provides a sufficiently accurate mapping from dark matter distribution to Lyman-alpha flux statistics
Invoked to compute flux power spectra from the evolved density and velocity fields

pith-pipeline@v0.9.0 · 5571 in / 1402 out tokens · 62299 ms · 2026-05-10T18:14:14.126238+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

despite near-degeneracy in the nonlinear matter power spectrum, the corresponding Ly α flux power spectra differ at the 10 percent level... arises from a strong suppression of small-scale velocity power in the Schrödinger–Poisson evolution
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

flux statistics cannot be fully characterized by the matter power spectrum alone, but depend sensitively on the dynamical evolution of the velocity field

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

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

[1]

Cosmological analysis of the DESI DR1 Lyman alpha 1D power spectrum

Almgren A. S., Bell . B., Lijewski M. J., Lukić Z., Andel E. V ., 2013, ApJ, 765, 39 Almgren A., et al., 2019, AMReX-Codes/amrex: AMReX 19.09, doi:10.5281/zenodo.2555438 Bird S., 2017, FSFE: Fake Spectra Flux Extractor, Astrophysics Source Code Library, record ascl:1710.012 (ascl:1710.012) Chaves-Montero J., et al., 2026, arXiv e-prints, p. arXiv:2601.214...

work page internal anchor Pith review Pith/arXiv arXiv doi:10.5281/zenodo.2555438 2013