{"id":"17094e04-d168-4500-bda7-3b8a25c3f506","arxiv_id":"2509.10431","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A quartic scalar field collapses with nearly the same critical exponent as a radiation fluid during primordial black hole formation, supporting the near universality of type II critical collapse.","lead":"This paper uses computer simulations to study how black holes form in the early universe when the energy is carried by a self-interacting scalar field rather than a hot radiation fluid. It finds the two cases collapse almost identically near the formation threshold, supporting a common approximation used in cosmology.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'extremely diluted' perfect fluid used to define the Misner–Sharp frame is never quantified or tested; if it develops order-unity fractional density perturbations, it can affect the lapse and contaminate the scalar-field critical exponent.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the diluted perfect fluid used to define the comoving frame could contaminate the scalar-field collapse. The central claim is that a quartic scalar field collapses with a critical exponent statistically consistent with radiation, which requires that the scalar-field runs measure the scalar field alone. The paper never states the initial fluid energy fraction nor demonstrates convergence as it is reduced. The mechanism I identified makes the concern concrete: Eq. (12h) lets the fluid pressure gradient set the lapse, and Eq. (12e) sources relative density perturbations whose amplitude does not scale with the fluid's initial density. Thus an 'extremely diluted' fluid is not automatically dynamically negligible, and the scalar-field exponent could be contaminated. The proposed test, varying the initial fluid fraction and refitting gamma_sf, would settle this directly. Secondary concerns reinforce the need for this test: the quoted errors are fit-only, the radiation exponent differs from Musco et al. by roughly 2.4 sigma of the quoted error, and the quartic coupling is fixed at one value with no variation. These issues do not by themselves invalidate the result, and the paper includes constraint monitoring and resolution checks, so CONDITIONAL remains the appropriate verdict.","tokens_in":14895,"tokens_out":9580,"duration_ms":91469,"concrete_test":"Repeat the scalar-field critical-collapse runs with the initial fluid energy fraction decreased by factors of 10 (for example 10^-4, 10^-6, 10^-8) while keeping all other parameters fixed, and refit gamma_sf in each case. If gamma_sf shifts by more than the quoted 0.0071 uncertainty as the fluid fraction tends to zero, or if the fluid's fractional density perturbation grows to order unity near criticality, the diluted fluid is not a test fluid and the reported exponent is contaminated. As a complementary check, monitor max(|delta-rho_pf|/rho_pf) during the approach to criticality and compare it with the scalar-field density contrast.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II.C.2 states that the scalar-field runs include a homogeneous, isotropic, extremely diluted perfect fluid to define the Misner–Sharp comoving frame, but no value for its initial energy fraction is given and no convergence test is reported. This matters because the fluid is not a passive observer: it enters the total stress-energy and, through Eq. (12h), phi' = -P'_pf/(rho_pf + P_pf), it directly sets the lapse. Even if the initial fluid density is tiny, the evolution equation (12e) sources fractional density perturbations of order unity, because the source term is proportional to (rho_pf + P_pf), making delta-rho_pf/rho_pf independent of the fluid's initial amplitude. Pressure gradients inherited from these perturbations can then make P'_pf/(rho_pf + P_pf) order unity, so the fluid's influence on phi need not vanish as the initial fluid density tends to zero. Since phi appears in every scalar-field evolution equation, a contaminated lapse can shift the measured exponent gamma_sf = 0.3401 +/- 0.0071 and the claimed 2-sigma separation from the radiation value. The paper also omits the fluid density perturbation from the scalar-field initial-data list and from the diagnostics shown in Fig. 4, so the reader cannot assess whether the reported scaling is a pure scalar-field result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses a Misner–Sharp, spherically symmetric, fully nonlinear relativistic code to study critical collapse of a quartic self-interacting scalar field and a radiation fluid in an expanding universe. It reports type II critical behavior in both cases, with critical exponents γ = 0.3474 ± 0.004 (radiation) and γ = 0.3401 ± 0.0071 (scalar field), and interprets the difference as about 2σ, supporting the near universality of the critical exponent in primordial black hole formation. The numerical methodology includes horizon detection, excision, Hamiltonian-constraint monitoring, and resolution tests.","tokens_in":15195,"tokens_out":5134,"duration_ms":47438,"significance":"If the result is robust, this is a useful first direct numerical test of the scalar-field/perfect-fluid correspondence in the critical collapse regime, with implications for PBH mass calculations in scalar-field-dominated early-universe scenarios. The paper's strengths are that it uses a fully nonlinear evolution with constraint monitoring, reports resolution dependence, and compares two matter models with the same code. However, the central quantitative claim depends on an unquantified artificial fluid in the scalar-field runs, and the stated 2σ separation is not supported by the quoted errors.","major_comments":[{"comment":"The diluted perfect fluid is not demonstrated to be dynamically negligible. The fluid is not a passive marker: it contributes to the stress-energy tensor and, through Eq. (12h), it directly fixes the lapse gradient. Eq. (12e) sources δρ_pf/ρ_pf through terms independent of the fluid's initial amplitude because the source is proportional to (ρ_pf+P_pf), so fractional perturbations of order unity can develop even for an 'extremely diluted' initial fluid. The manuscript gives no value for ρ_b0pf, does not include fluid perturbations in the initial-data list, and shows no diagnostic for the fluid density in Fig. 4. Please quantify ρ_b0pf and demonstrate by explicit convergence tests, for example by varying ρ_b0pf over several orders of magnitude, that γ_sf is unaffected; otherwise the measured exponent may be contaminated by the gauge fluid.","section":"II.C.2 and Eqs. (12e), (12h)"},{"comment":"The claimed 'about 2σ' difference between the exponents is not consistent with the quoted errors. The difference is 0.3474 − 0.3401 = 0.0073; the quadrature-summed error is sqrt(0.004^2 + 0.0071^2) ≈ 0.0082, giving about 0.9σ. Please correct this statement and avoid interpreting the result as evidence for a statistically significant difference.","section":"Abstract and Sections V.B, VI"},{"comment":"The quoted uncertainties in γ are fit errors only. The conclusion of 'near universality' rests on a small difference between exponents, so systematic uncertainties from the choice of fitting range, resolution, threshold bracket in Eq. (36), horizon-mass extrapolation, and the excision criterion need to be estimated. Without such estimates, the comparison between γ_sf and γ_rad cannot be assigned a reliable significance.","section":"IV.B and V.A"}],"minor_comments":[{"comment":"There are typos: 'dominateed' in the Introduction, 'collaps' in the Introduction, and 'threashold' in Section V.A.","section":"Throughout"},{"comment":"The initial condition for the scalar field is normalized as ~Π_b0^2/2 = 1 − ρ_b0pf, but ρ_b0pf is never specified; please state its value and justify that it is too small to alter the background dynamics.","section":"II.C.2, Eq. (16)"},{"comment":"The threshold is defined as the midpoint between p_min_bh and p_max_no-bh, but the bracketing procedure and how close to threshold the simulations are run are not described; please specify the bisection or bracketing algorithm.","section":"V.A, Eq. (36)"},{"comment":"The set of resolutions used for the convergence tests is not stated in the captions; please list the dA values and the corresponding line styles.","section":"Figs. 2 and 5"},{"comment":"The notation for the horizon velocity criterion is inconsistent: Section IV.B uses v_c − 1 < 3 × 10^-3, while Section III.C and Appendix A define v_AH; please unify the notation.","section":"III.C and IV.B"},{"comment":"The footnote acknowledges that w = 1/3 is only an effective equation of state for the quartic scalar field; near criticality the field may not oscillate rapidly, so please discuss how this could affect the interpretation of the comparison with radiation.","section":"IV.A, footnote 2"}],"recommendation":"major_revision","confidential_remarks":"The main risk is the artificial fluid in the scalar-field runs; if the authors cannot show that γ_sf is robust to the value of ρ_b0pf, the central claim fails. The '2σ' statement should also be corrected. I would not reject at this stage because the numerical framework is sound and the issue is testable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the central numerical claim—that a quartic self-interacting scalar field shows type II critical collapse with an exponent statistically indistinguishable from radiation—is plausible and mostly supported by the simulations. What is genuinely new is the first extraction of the exponent for this matter model in a cosmological setting. The authors deserve credit for the basics: resolution studies, Hamiltonian constraint monitoring, and a clean scaling law. This is not careless numerics.\n\nThe soft spots are real but narrower than the abstract implies. First, the claimed ~2σ separation between the two exponents is not consistent with the quoted errors: 0.3401±0.0071 and 0.3474±0.004 differ by about 0.9σ. That is a presentational error, easy to fix, but it changes the conclusion from \"slightly different\" to \"consistent within errors.\"\n\nSecond—and this is the one I would push on—the homogeneous, \"extremely diluted\" perfect fluid used to define the comoving frame is never quantified. Eq. (12h) makes the lapse gradient directly proportional to the fluid's pressure gradient divided by its density. That ratio is independent of how small the initial fluid energy density is, once the fluid develops order-unity fractional perturbations, which the evolution equations can produce. So the fluid is not automatically a passive spectator. The paper should show that the measured gamma is insensitive to the initial fluid density, or that the fluid's fractional perturbations stay small. As written, the scalar-field exponent could be contaminated by the gauge-defining fluid. This is the difference between a good conference talk and a result I would build on.\n\nThere are also smaller gaps: uncertainties are fit errors only, the quartic coupling is fixed at λ̂4 = 10, no code or data is released, and the dependence on the initial profile width is not explored. None of these are fatal, but together they mean this is a numerical confirmation, not a definitive verification.\n\nWho this is for: people studying PBH formation from oscillating scalar fields and anyone using the radiation-fluid approximation in critical collapse. It deserves a serious referee. I would send it out with the diluted-fluid question and the sigma discrepancy as required revisions. If the fluid test comes back clean, the core result stands.","headline":"Useful first measurement of the quartic-scalar critical exponent, but the abstract overstates the difference from radiation and the gauge-defining fluid is never tested.","tokens_in":15701,"tokens_out":4390,"would_cite":false,"duration_ms":44489,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A quartic scalar field and a radiation fluid collapse into primordial black holes with nearly identical critical exponents.","keywords":["primordial black holes","critical gravitational collapse","critical exponent","quartic scalar field","type II collapse","Misner-Sharp formalism","radiation fluid","near universality"],"falsifier":"Run the same scalar-field collapse with the homogeneous background fluid density reduced by factors of 10 while holding the scalar perturbation fixed; if the fitted exponent $\\gamma$ moves, or if the fluid's share of the Misner-Sharp mass near horizon formation grows, the quoted scalar-field exponent is contaminated by the frame-defining fluid.","tokens_in":14704,"feed_emoji":"🕳️","tokens_out":8524,"duration_ms":69338,"temperature":0.7,"pith_summary":"The paper asks whether a universe dominated by a self-interacting scalar field with a quartic potential forms primordial black holes in the same critical way as a radiation-filled universe. Using fully nonlinear, spherically symmetric simulations in the Misner-Sharp formalism, it finds that the black hole mass near threshold follows the same type II power law in both cases, with critical exponents agreeing to within about $2\\sigma$. If correct, this extends the common shortcut of modelling oscillating scalar fields as perfect fluids into the strongly nonlinear collapse regime and supports the near universality of the critical exponent in primordial black hole formation. The paper reports a small residual difference between the exponents rather than claiming exact equality.","feed_headline":"Quartic scalar collapse matches radiation's critical exponent","feed_subtitle":"Numerical relativity finds near-universal type II scaling for primordial black hole formation across matter models.","key_machinery":"The central object is the power-law scaling of the black hole mass near threshold, $M_{\\rm BH}\\propto(p-p_{\\rm th})^\\gamma$, in type II critical collapse. The machinery that extracts it is the Misner-Sharp formulation of spherical relativistic collapse, evolved with a fourth-order Runge-Kutta method and finite differences, with a scalar field plus an extremely diluted homogeneous perfect fluid that defines the comoving frame. Black hole formation is detected through the apparent-horizon condition $\\Theta_+=0$; the final mass is read from the Misner-Sharp mass once the apparent-horizon velocity satisfies $v_{\\rm AH}\\approx 1$, and $\\gamma$ is obtained by fitting the mass against $p-p_{\\rm th}$.","core_discovery":"The paper's central claim is that a quartic self-interacting scalar field, $V(\\psi)\\propto\\psi^4$, exhibits type II critical collapse (black hole mass becoming arbitrarily small near threshold and scaling as a power law) with critical exponent $\\gamma = 0.3401 \\pm 0.0071$, statistically close to the radiation-fluid value $\\gamma = 0.3474 \\pm 0.004$ measured with the same code. The threshold amplitudes are also close: $p_{\\rm th}^{\\rm sf}=0.0271205$ versus $p_{\\rm th}^{\\rm rad}=0.0270895$. The authors present this as the first explicit verification that the scalar-field/perfect-fluid correspondence extends into the critical collapse regime, at least for the quartic potential, while emphasising that the two exponents differ by about $2\\sigma$ and that the quadratic case departs from dust-like behaviour.","pith_inferences":["Varying the potential exponent between the quadratic and quartic cases, or changing the quartic coupling, could map how $\\gamma$ drifts with the effective equation of state and reveal whether the $2\\sigma$ gap is a physical trend rather than a numerical artefact.","The frame-defining homogeneous fluid could be eliminated by reformulating the evolution in a different gauge; if the same exponent survives, the near-universality claim would no longer rest on that technical crutch.","The compaction function oscillates near threshold in the scalar-field runs, so peak-theory estimates of primordial black hole abundance that rely on a single maximum of the compaction function may need oscillation-aware prescriptions."],"forward_implications":["For quartic-potential scenarios, primordial black hole abundance calculations that approximate the field as a radiation fluid should capture the near-threshold scaling, with the $\\sim2\\%$ exponent difference entering as a small systematic uncertainty.","The same code reproduces the known radiation critical exponent, so the comparison baseline is not an artefact of the numerical setup.","The quadratic scalar-field case departs from dust-like behaviour, so the fluid analogy is potential-dependent and cannot be assumed for all scalar-field models.","The $2\\sigma$ separation between the two exponents defines a concrete target: higher-resolution runs will either close the gap or establish a genuine matter-model dependence of $\\gamma$."],"supporting_citations":[{"why":"It establishes the universal power-law scaling in gravitational collapse that this paper tests for a quartic scalar field.","marker":"[1]"},{"why":"It supplies the radiation-era critical exponent benchmark ($\\gamma=0.357$) and the earlier demonstration of critical collapse in the radiative era.","marker":"[7]"},{"why":"It defines the type I/type II classification of critical collapse that frames the paper's claims.","marker":"[12]"},{"why":"It provides the numerical code, initial-condition construction, and earlier scalar-field collapse results that this work extends.","marker":"[21]"},{"why":"It formulates the spherically symmetric relativistic collapse equations used throughout the simulations.","marker":"[28]"},{"why":"It provides the method for constructing super-horizon growing-mode initial conditions.","marker":"[34]"},{"why":"It gives the gradient-expansion initial conditions for perfect-fluid dominated universes that the paper adapts to its setup.","marker":"[35]"},{"why":"It supplies the causal classification of apparent horizons used to decide when the black hole mass has stabilised.","marker":"[39]"},{"why":"It provides the Kreiss-Oliger dissipation scheme that controls numerical oscillations in the finite-difference evolution.","marker":"[43]"}],"fun_headline_variants":["Quartic scalar collapse: critical exponent within 2σ of radiation","Scalar field vs radiation: near-universal PBH critical scaling","Numerical relativity: quartic scalar collapse near radiation's critical exponent","Type II collapse nearly universal: quartic scalar matches radiation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The scalar-field simulations carry a tiny homogeneous perfect fluid purely to define the comoving frame, and the paper does not show that this fluid stays dynamically negligible at the critical threshold.","fun_headline_variants_meta":{"raw":{"variants":["Quartic scalar collapse: critical exponent within 2σ of radiation","Scalar field vs radiation: near-universal PBH critical scaling","Numerical relativity: quartic scalar collapse near radiation's critical exponent","Type II collapse nearly universal: quartic scalar matches radiation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000797,"raw_usage":{"total_tokens":3482,"prompt_tokens":892,"completion_tokens":2590,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":508,"completion_tokens_details":{"reasoning_tokens":2516}},"tokens_in":508,"tokens_out":2590,"duration_ms":17497,"temperature":1.0,"reasoning_tokens":2516,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:54:13.651738+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same scalar-field collapse with the homogeneous background fluid density reduced by factors of 10 while holding the scalar perturbation fixed; if the fitted exponent $\\gamma$ moves, or if the fluid's share of the Misner-Sharp mass near horizon formation grows, the quoted scalar-field exponent is contaminated by the frame-defining fluid.","supporting_citations":[{"cited_title":"Long- wavelength nonlinear perturbations of a complex scalar field,","cited_arxiv_id":null,"evidence_quote":"It provides the numerical code, initial-condition construction, and earlier scalar-field collapse results that this work extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It establishes the universal power-law scaling in gravitational collapse that this paper tests for a quartic scalar field."},{"cited_title":"Computa- tions of primordial black-hole formation,","cited_arxiv_id":null,"evidence_quote":"It supplies the radiation-era critical exponent benchmark ($\\gamma=0.357$) and the earlier demonstration of critical collapse in the radiative era."},{"cited_title":"Boson stars driven to the brink of black hole formation,","cited_arxiv_id":null,"evidence_quote":"It defines the type I/type II classification of critical collapse that frames the paper's claims."},{"cited_title":"Inflaton clusters and inflaton stars,","cited_arxiv_id":null,"evidence_quote":"It formulates the spherically symmetric relativistic collapse equations used throughout the simulations."},{"cited_title":"Identifying the most crucial parameters of the initial curvature profile for primordial black hole formation,","cited_arxiv_id":null,"evidence_quote":"It provides the method for constructing super-horizon growing-mode initial conditions."},{"cited_title":"Identifying the most crucial parameters of the ini- tial curvature profile for primordial black hole forma- tion,","cited_arxiv_id":null,"evidence_quote":"It supplies the causal classification of apparent horizons used to decide when the black hole mass has stabilised."},{"cited_title":"Horizon boundary condition for black hole space- times,","cited_arxiv_id":null,"evidence_quote":"It provides the Kreiss-Oliger dissipation scheme that controls numerical oscillations in the finite-difference evolution."}],"review_version":2}