{"id":"086416b7-dbe1-4f72-95b7-88bd5b8973a7","arxiv_id":"2507.01626","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In void environments, the galaxy luminosity function is fitted with Schechter parameters that vary with the void density contrast and the growth function, improving on the McNaught-Roberts model.","lead":"This paper proposes a new way to describe how many galaxies of each brightness live inside cosmic voids, using the void's density contrast and the growth of cosmic structures over time. The authors test the model on two galaxy surveys and report that it fits the data better than the previous standard formula.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 14's D mixes nonlinear observed density contrasts with a linear EST barrier, so the fitted exponents lose their claimed physical meaning unless Eq. 2 is applied.","rationale":"The reader's weakest-assumption analysis already identified the linear/nonlinear D ambiguity, and my reading of the manuscript confirms that this is the load-bearing point. The paper's central claim is not merely that a three-parameter formula fits void luminosity functions, but that the environment dependence enters through the EST two-barrier ratio D, giving the fit a theory-based meaning. That meaning breaks down if the observed Eulerian density contrasts are used without the Eq. 2 mapping, or if the fixed linear threshold is used without environment variation. The concern is internal to the paper: the text invokes both 'linear density contrast of the environment' in the abstract and 'nonlinear critical density contrast' in Table 2 without reconciling them. Independent support is limited: there is no released code or data, the redshift evolution in Eq. 15 has a fitted η rather than a predicted value, and the environmental analysis is restricted to 2dFGRS with only three void/density bins. These limitations do not by themselves refute the empirical fit, but they reinforce that the physical interpretation is fragile. The appropriate response is to require a clarification and refit under the correct linear mapping, so the reader's conditional verdict stands. I therefore recommend no change to the verdict: it remains CONDITIONAL pending the D-coordinate check and the accompanying justification.","tokens_in":10582,"tokens_out":4746,"duration_ms":52160,"concrete_test":"Recompute the three D inputs from Table 2 using Eq. 2: for each nonlinear δ_NL, solve for the linear δ_L via δ_L = 1.68647 − 1.35/(1+δ_NL)^(2/3) − 1.12431/(1+δ_NL)^(1/2) + 0.78785/(1+δ_NL)^0.58661, then set D = |δ_L|/(1.686+|δ_L|). Refit Eq. 14 to the 2dFGRS data with these D values and compare ζ2, ζ4, ζ5, χ², AIC, and BIC with Table 3 and Fig. 2. If the best-fit parameters or the model-comparison statistics shift by more than the quoted 68.3% uncertainties, the current D inputs are coordinate-mismatched and the central physical claim is unsupported as written.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The environmental model's only environment-dependent input is D = |δv|/(δc+|δv|) in Eq. 14, which is motivated in §2.1 as the two-barrier ratio built from the linear void-formation threshold δv ≈ −2.81 and halo-collapse threshold δc ≈ 1.686. The data used for the fit, however, are labeled in Table 2 and Fig. 1 as nonlinear (Eulerian) density contrasts: δv = −0.90, −0.75, −0.43. The paper provides the Lagrangian-to-Eulerian mapping in Eq. 2 but never applies it to convert these values before forming D. If the nonlinear values are inserted directly, D is not the EST barrier ratio and the fitted exponents ζ2, ζ4, ζ5 in Table 3 do not carry the claimed theoretical meaning; if instead the constant linear threshold δv = −2.81 is used, D becomes environment-independent and the model has no environmental dependence at all. This ambiguity propagates directly into the reported χ², AIC, and BIC improvements and into the 'McNaught-modified' comparison, so it is not a minor bookkeeping issue but the central interpretive step of the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an excursion-set-inspired parameterization of the Schechter luminosity function for galaxies in voids, with the Schechter parameters written as functions of the two-barrier ratio D = |δv|/(δc+|δv|) and, separately, a redshift dependence of the form Φ(M,z) ∝ exp(−η/D(z)^2)/D(z), where D(z) is the linear growth function. The environmental model is fitted to 2dFGRS luminosity functions in three underdense environments and compared against the McNaught-Roberts et al. (2014a) model and a modified version of it, reporting lower χ², AIC, and BIC. The redshift model is fitted to GAMA void luminosity functions in three redshift bins. The central claim is that EST-motivated environmental and redshift dependencies provide a better description of the void galaxy luminosity function than previous empirical parameterizations.","tokens_in":10839,"tokens_out":8393,"duration_ms":94027,"significance":"If the environmental model holds, it offers a compact, EST-inspired description of void galaxy luminosity functions and a modest improvement over an existing empirical model using the same data and the same number of free parameters; the AIC/BIC comparison in Fig. 2 is therefore a fair model-selection exercise. The redshift part is more exploratory: the GAMA data provide only one void density interval, as the paper acknowledges, and the single fitted parameter η makes Fig. 4 a demonstration of a functional form rather than an out-of-sample prediction. The paper ships no code or machine-checked derivations; its value rests on the clarity and correctness of the physical mapping between the fitted density contrasts and the EST barrier ratio, which is currently unresolved.","major_comments":[{"comment":"The parameter D = |δv|/(δc+|δv|) is introduced in §2.1 as an EST two-barrier ratio built from the linear thresholds δv ≈ −2.81 and δc ≈ 1.686, but Table 2 and Fig. 1 label the fitted environments by nonlinear (Eulerian) density contrasts δv = −0.90, −0.75, −0.43. The manuscript never applies the Lagrangian-to-Eulerian mapping of Eq. (2) to these values before forming D. If the nonlinear values are inserted directly, D is not the EST barrier ratio and the fitted exponents ζ2, ζ4, ζ5 in Table 3 lose their stated theoretical meaning; if one instead uses the constant linear threshold δv = −2.81, D becomes environment-independent and Eq. (14) has no environmental dependence at all. Please report the D values actually used, state explicitly whether δv is linear or nonlinear, and if the values in Table 2 are nonlinear, convert each one with Eq. (2) and refit; the χ², AIC, and BIC results in Fig. 2 must be rechecked under the corrected D.","section":"§3.1, Eq. (14), Table 2, Eq. (2)"},{"comment":"The redshift model is presented as a test of Eq. (15), but η is a free parameter fitted to the same GAMA void luminosity functions shown in Fig. 4; the agreement is therefore a fitting result, not an independent validation. The derivation of Eq. (15) from Eq. (5) also relies on unstated assumptions, namely that δc(z) − δv(z) scales as 1/D(z) and that the relevant variance combination is redshift-independent. The paper itself acknowledges a deviation at the faint end. Please provide a quantitative comparison of Eq. (15) against, for example, a no-evolution model or the redshift-dependent parameterization used by Loveday et al. (2012), and state explicitly whether the model is intended as a description or as a predictive test until independent data are used.","section":"§3.2, Eq. (15), Fig. 4"},{"comment":"The halo-in-void mass function in Eq. (5) is introduced as the theoretical basis, but it is not used in the fits: Eq. (9) employs the global Sheth-Tormen mass function n(M), and the environmental dependence in Eq. (14) enters only through the ad hoc replacement of the density contrast by D. The mapping from Eq. (5) to the specific functional forms in Eq. (14) is not derived. Either derive Eq. (14) from the halo-in-void mass function under stated approximations, or clearly present Eq. (14) as an empirical fitting function whose parameters are merely inspired by EST; the current text claims an EST grounding that the equations do not actually provide.","section":"§2.2–§2.3, Eqs. (5), (7), (9), (14)"}],"minor_comments":[{"comment":"The symbol D denotes two different quantities: the barrier ratio D ≡ |δv|/(δc+|δv|) in Eq. (14) and the linear growth function D(z) in Eq. (15). This is confusing in §3.2, where both appear nearby; please rename one of them, for example using B for the barrier ratio.","section":"§3.1–§3.2, Eqs. (14) and (15)"},{"comment":"Eq. (1) defines n(M,z) as the cumulative number density of halos with mass greater than M, but Eqs. (7) and (9) treat n(M) as a differential mass function. Please clarify the notation and use dn/dM where appropriate.","section":"§2.1, Eq. (1)"},{"comment":"The table and figure block around Fig. 2 repeats the AIC/BIC values and contains a stray “Model 2” label; please fix the typesetting so the reader can unambiguously identify which row corresponds to the new model, the McNaught-modified model, and the original McNaught-Roberts model.","section":"§3.1, Fig. 2"},{"comment":"The caption of Fig. 4 states the void selection as −1 ≤ δv ≤ −0.75, but Table 2 lists discrete values −0.90, −0.75, and −0.43; please clarify how the redshift-binned GAMA sample maps onto the environment definitions used in the 2dFGRS analysis.","section":"§3.2, Fig. 4"},{"comment":"The references to McNaught-Roberts et al. (2014a) and (2014b) appear to cite the same MNRAS volume and page; please verify that the two citations are indeed distinct works and cite them correctly.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the central model-comparison idea is reasonable, but the environmental claim rests on the D interpretation, which is currently ambiguous between linear and nonlinear density contrasts. This is fixable by applying Eq. (2) and refitting, so I would not recommend rejection; however, the redshift section needs to be framed as a fit with unverified assumptions rather than as a validated theory."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper gives a new, compact parameterization of the void galaxy luminosity function in terms of the EST two-barrier ratio D, and it genuinely outperforms the McNaught-Roberts model on the same 2dFGRS data with the same number of free parameters. That part holds up. The weak spot is that the D values used in the fit are built from nonlinear density contrasts while the EST barriers are linear, so the theoretical interpretation is shaky.\n\nWhat is actually new: the D-dependence of the Schechter parameters (Eq. 14) and the exp(-eta/D(z)^2)/D(z) redshift form (Eq. 15) are not in the cited McNaught-Roberts model. The authors also show that simply swapping delta_m for D in the old model cuts chi^2 by 35%, which is a nice consistency check. The AIC/BIC comparison is fair: same data, same number of free parameters. They use real survey data (2dFGRS, GAMA) and report error bars on the fits.\n\nThe soft spots, in proportion: first, the D ambiguity. Table 2 lists the environment density contrasts as nonlinear (Eulerian) values, delta_v = -0.90, -0.75, -0.43, while the EST two-barrier ratio is derived from linear thresholds, delta_v ~ -2.81 and delta_c ~ 1.686. The paper gives the Lagrangian-to-Eulerian mapping in Eq. 2 but never applies it. If you plug the nonlinear values directly into D, you are not computing the quantity EST motivates. If you instead use the constant linear threshold, D stops being environment-dependent and the model has no environmental lever at all. This is not just bookkeeping; it determines what the fitted exponents zeta2, zeta4, zeta5 mean physically. The empirical fit may still be better than the alternatives, but the theoretical story needs to be repaired, either by converting to linear contrasts or by reframing the model as a purely empirical fitting form.\n\nSecond, the redshift model is a fit, not a prediction. eta is a free parameter fitted to the same GAMA data that Fig. 4 compares against, so the agreement is a fitting result. The authors are honest about the faint-end deviation, but the word 'theory' is doing too much work here. This is a minor issue if the model is presented as an empirical fitting formula; it is a larger issue if the goal is to use this as a cosmological probe.\n\nThird, no code or data release, which makes it harder to check the fitting procedure. That is a smaller point.\n\nBottom line: this is a subfield paper with limited scope, but the empirical claim is real and the data comparison is fair. The density-coordinate issue is the main thing standing between a conditional acceptance and a clean one. I'd send it to a competent referee focused on large-scale structure and voids.\n\nRecommendation: accept with major revision, or conditional accept, pending the authors either applying the linear-to-nonlinear mapping or explicitly recasting Eq. 14 as an empirical fit with no EST-based physical interpretation.","headline":"A modest but real empirical improvement to void luminosity function modeling, undermined by an unresolved linear/nonlinear density-coordinate ambiguity in the D parameter.","tokens_in":11349,"tokens_out":3458,"would_cite":false,"duration_ms":37916,"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":"The Schechter parameters of void galaxy luminosity functions follow power laws in the excursion-set two-barrier ratio, and this parameterization fits 2dFGRS data better than the previous density-based model.","keywords":["galaxy luminosity function","cosmic voids","excursion set theory","Schechter function","environment dependence","redshift evolution","growth function","2dFGRS"],"falsifier":"Use Eq. (2) to convert the nonlinear contrasts $\\delta_v=-0.90,-0.75,-0.43$ to linear values, recompute $D$, refit Eq. (14) to the same 2dFGRS data, and compare $\\chi^2$, AIC, and BIC; if the advantage over the density-based model disappears, the central environmental claim is not supported.","tokens_in":10367,"feed_emoji":"🔭","tokens_out":8351,"duration_ms":87235,"temperature":0.7,"pith_summary":"The paper tries to establish that the galaxy luminosity function in cosmic voids is controlled by the excursion-set two-barrier ratio $D = |\\delta_v|/(\\delta_c+|\\delta_v|)$, not merely by the local density contrast. Writing the Schechter parameters as power laws in $(1+D)$ and fitting to 2dFGRS data in three underdense environments gives $\\chi^2 = 20.8$, compared with 22.2 for the previous model with $D$ substituted for the density and 34.4 for the original density-based model, and the AIC/BIC ranking follows the same order. The paper further claims that the redshift evolution of the void luminosity function enters through the linear growth function $D(z)$ as $\\Phi(M,z)\\propto \\exp(-\\eta/D(z)^2)/D(z)$, with $\\eta$ fitted to GAMA void data. If these claims hold, the luminosity function becomes a direct observable link between galaxy properties in voids and the physics of halo and void formation, and a possible cosmological probe.","feed_headline":"Void galaxy luminosity fits the two-barrier ratio better","feed_subtitle":"Excursion-set barrier ratio D reshapes the Schechter parameters and beats the standard density model on 2dFGRS data.","key_machinery":"The load-bearing object is the ratio $D\\equiv|\\delta_v|/(\\delta_c+|\\delta_v|)$, which measures how close the void-formation barrier $\\delta_v\\approx -2.81$ is to the halo-collapse barrier $\\delta_c\\approx 1.686$; in excursion set theory this ratio controls which voids survive being swallowed by collapsing halos. The paper makes $D$ the environmental variable in the Schechter function, so the three Schechter parameters become power laws in $(1+D)$. The second mechanism is the linear growth function $D(z)$, imported through the halo mass function, which turns the redshift dependence of the luminosity function into $\\Phi(M,z)\\propto\\exp(-\\eta/D(z)^2)/D(z)$ and pushes all of that dependence into $\\phi_*$.","core_discovery":"On the paper's own terms, the discovery is that the environment dependence of the void galaxy luminosity function is not an arbitrary empirical trend but follows from the distance between the two formation barriers in excursion set theory. The claimed result is $\\phi_* = \\zeta_1(1+D)^{\\zeta_2}$, $M_* = \\zeta_3(1+D)^{\\zeta_4}$, $\\alpha = \\zeta_5(1+D)$, where $D\\equiv|\\delta_v|/(\\delta_c+|\\delta_v|)$ and the $\\zeta_i$ are fitted to 2dFGRS; this form outperforms the earlier density-based Schechter model by $\\chi^2$, AIC, and BIC. For redshift, the claim is that the halo mass function conveys the evolution through the growth factor, giving $\\Phi(M,z)\\propto \\exp(-\\eta/D(z)^2)/D(z)$ and concentrating the redshift dependence in $\\phi_*$, with $\\eta = 0.068^{+0.017}_{-0.018}$ from GAMA void luminosity functions.","pith_inferences":["The same $D$-parameterization could be tested in larger void samples from current spectroscopic surveys to see whether the fitted exponents $\\zeta_2$ and $\\zeta_4$ are universal or survey-dependent; this is an extension, not a paper claim.","Because the paper never applies its own linear-to-nonlinear mapping to the Table 2 density contrasts, the numerical value of $D$ feeding the fit is ambiguous; a recalibrated fit might shift the exponents even if the qualitative ranking of models survives.","If the redshift formula holds, the luminosity function of voids could be used as a growth-rate probe, complementing clustering-based growth measurements; this direction is implicit in the paper rather than demonstrated there."],"forward_implications":["Void luminosity functions in sparse, populous, and underdense regions can be described by one shared parameterization in the two-barrier ratio rather than fitted separately.","The redshift dependence of the void luminosity function is concentrated in $\\phi_*$, so measuring $\\phi_*$ in redshift bins is a direct route to the growth function and hence to the cosmological model.","The reported AIC/BIC improvement implies that in underdense regions the two-barrier ratio is a more informative environmental coordinate than the Eulerian density contrast.","Neglecting satellite galaxies is a safe simplification for voids, so the model connects cleanly to the halo mass function there."],"supporting_citations":[{"why":"Supplies the two-barrier void formation threshold $\\delta_v \\approx -2.81$ and the void survival distribution that motivates the ratio $D$.","marker":"Sheth & Van De Weygaert (2004)"},{"why":"Provides the baseline density-dependent Schechter model and the GAMA void luminosity function data used for the redshift fit.","marker":"McNaught-Roberts et al. (2014a)"},{"why":"Supplies the 2dFGRS luminosity function data in three underdense environments used for the environmental fit.","marker":"Croton et al. (2005)"},{"why":"Supplies the Schechter functional form whose parameters the paper makes environment- and redshift-dependent.","marker":"Schechter (1976)"},{"why":"Supplies the mass-luminosity relation that lets the paper convert the halo mass function into a luminosity function.","marker":"Vale & Ostriker (2004)"},{"why":"Supplies the excursion set framework and the first-crossing distribution behind the halo mass function and the two-barrier argument.","marker":"Bond et al. (1991)"}],"fun_headline_variants":["Excursion set theory rewires void galaxy luminosity","Barrier ratio D beats density model for void galaxies","Void galaxy brightness traces two-barrier distance","New void galaxy luminosity model wins on BIC and AIC","Void galaxy luminosity: excursion set barrier ratio works"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the observed nonlinear void density contrasts in Table 2 can be identified directly with the linear barrier $\\delta_v$ when forming $D=|\\delta_v|/(\\delta_c+|\\delta_v|)$; if the linear-to-nonlinear mapping must be applied first, the fitted D values and the claimed model superiority lose their stated physical meaning.","fun_headline_variants_meta":{"raw":{"variants":["Excursion set theory rewires void galaxy luminosity","Barrier ratio D beats density model for void galaxies","Void galaxy brightness traces two-barrier distance","New void galaxy luminosity model wins on BIC and AIC","Void galaxy luminosity: excursion set barrier ratio works"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000165,"raw_usage":{"total_tokens":1292,"prompt_tokens":1027,"completion_tokens":265,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":643,"completion_tokens_details":{"reasoning_tokens":189}},"tokens_in":643,"tokens_out":265,"duration_ms":3642,"temperature":1.0,"reasoning_tokens":189,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:46:57.319062+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Use Eq. (2) to convert the nonlinear contrasts $\\delta_v=-0.90,-0.75,-0.43$ to linear values, recompute $D$, refit Eq. (14) to the same 2dFGRS data, and compare $\\chi^2$, AIC, and BIC; if the advantage over the density-based model disappears, the central environmental claim is not supported.","supporting_citations":[{"cited_title":"K., & Van De Weygaert, R","cited_arxiv_id":null,"evidence_quote":"Supplies the two-barrier void formation threshold $\\delta_v \\approx -2.81$ and the void survival distribution that motivates the ratio $D$."},{"cited_title":"J., Farrar, G","cited_arxiv_id":null,"evidence_quote":"Supplies the 2dFGRS luminosity function data in three underdense environments used for the environmental fit."},{"cited_title":"1976, Astrophysical Journal, Vol","cited_arxiv_id":null,"evidence_quote":"Supplies the Schechter functional form whose parameters the paper makes environment- and redshift-dependent."},{"cited_title":"2004, Monthly Notices of the Royal Astronomical Society, 353, 189","cited_arxiv_id":null,"evidence_quote":"Supplies the mass-luminosity relation that lets the paper convert the halo mass function into a luminosity function."},{"cited_title":"1991, ApJ, Part 1 (ISSN 0004-637X), vol","cited_arxiv_id":null,"evidence_quote":"Supplies the excursion set framework and the first-crossing distribution behind the halo mass function and the two-barrier argument."}],"review_version":1}