{"id":"72c5efc4-9093-4173-a365-bad18983ef44","arxiv_id":"2412.19898","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":17,"one_line_summary":"A new halo occupation model that uses galaxy luminosity as well as clustering jointly fits [O II] emitters at z=1.193 and 1.471 and links them to Milky-Way-scale halos.","lead":"This paper fits a new model linking [O II]-emitting galaxies to the dark matter halos they live in, using both how clustered the galaxies are and how bright they are. The model matches Subaru telescope data at two redshifts and suggests these galaxies may grow into Milky-Way-sized halos by today.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The LHMR is assumed to be a deterministic double power law fitted to the same LF it is asked to predict; if the true L(Mh) relation has different shape or scatter, the inferred M_LF^min and satellite fractions can be biased.","rationale":"The reader's weakest assumption is exactly the point that carries the strongest weight: the deterministic double-power-law LHMR, with no scatter, is fitted to the same LFs it is used to reproduce, and no independent or out-of-sample check validates the shape. I read the paper in good faith: the fits are acceptable, the no-Gaussian ablation is a real internal consistency test, and the paper does not claim machine-level verification or ship code. However, the central quantitative claims - M_LF^min ~ 11.6, f_sat ~ 0.2-0.3, and the EPS-based Milky-Way-descendant statement - all propagate through Eq. 11 into the LF and occupation predictions. A flexible-shape or scatter test is the single check that would settle whether this assumption is responsible for the inferred parameters. The reader's conditional verdict remains appropriate: the paper is promising but needs this validation before its strong claims are accepted. I therefore recommend no change to the verdict, while emphasizing that this specific test should be run.","tokens_in":44725,"tokens_out":4979,"duration_ms":57819,"concrete_test":"Re-run the HOD MCMC for both NB816 and NB921 with the deterministic double power law of Eq. 11 replaced by a flexible non-parametric spline in log10 Mh, or augmented by an explicit intrinsic lognormal scatter in log10 L at fixed Mh, using the same ACF and LF likelihoods and the same flat priors. If the posterior of log10(M_LF^min) shifts by more than ~0.2 dex, or f_sat by more than ~0.05, relative to Table 3, the double-power-law functional form is not validated by the data and the claimed constraints are not robust. Report the resulting ACF and LF chi2 values to see whether the simpler deterministic form is actually preferred.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The inferred ELG-halo connection rests on the luminosity-to-halo-mass relation f_LHMR(Mh) in Eq. 11, a deterministic double power law with no explicit scatter. This function enters every central occupation through Eq. 10 and every LF prediction through Eqs. 10, 14, and 21. In Eq. 10 only the error-function term depends on the threshold Lth; the Gaussian term cancels identically in the LF difference in Eq. 21, so the LF constraints act almost entirely on the width/slope of the erfc transition and on the satellite term, with the LHMR shape assumed rather than independently tested. The LF is used both as a fitting constraint and as the reproduced output, and the bright-end LF uncertainties are large, so the double-power-law shape at the massive end is weakly constrained; Section 5.2.1 concedes that this end is statistically difficult to constrain. If the true relation deviates from the assumed form or has intrinsic scatter that is degenerate with sigma_logL, the posterior central occupations, M_LF^min, satellite fractions, and the Milky-Way-descendant conclusion in Section 5.1.1 could all be biased.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new HOD framework for [O II] emission-line galaxies in which the central occupation is tied to a double-power-law luminosity-to-halo-mass relation (LHMR, Eq. 11). This allows the model to predict luminosity functions as well as angular correlation functions. The authors fit the model jointly to ACFs and LFs of HSC NB816 and NB921 [O II] emitters at z=1.193 and z=1.471, compare with the Geach et al. (2012) HOD model, and report good fits (best-fit reduced chi2 around 0.8-0.9 in Table 3, with a text value of 0.79/0.86 in Section 4.3.2). They infer characteristic halo masses log10(M_LF^min/h^-1 Msun) ~ 11.5-11.7, satellite fractions ~0.2-0.3, and use extended Press-Schechter evolution to argue that the z~1.4 sample may evolve into Milky-Way-like halos. The paper also compares derived LHMR, baryon conversion efficiency, and SHMR with literature and with IllustrisTNG satellite fractions.","tokens_in":45081,"tokens_out":6930,"duration_ms":74621,"significance":"If the central claim is sustained, the framework is a useful step toward ELG HOD models that connect line luminosity to halo mass, with clear applications to mock catalogues for PFS, DESI, Euclid, and Rubin. The observational analysis is careful in several respects: jackknife covariance matrices with Hartlap correction, integral-constraint correction, explicit contamination modelling via f_ELG(M_h), and consistency checks against spectroscopic contamination estimates and IllustrisTNG. The public availability of the underlying halomod and emcee infrastructure also aids reproducibility. However, the load-bearing claim that the LF constraint makes the HOD constraints 'more robust' is not demonstrated by an ablation, and the luminosity-function fit is not an independent validation because the same data enter the likelihood that generates the model LF. These gaps prevent the paper from fully establishing its headline claim, although they are addressable with additional analysis.","major_comments":[{"comment":"The Gaussian term in the central occupation (Eq. 10) is independent of the threshold L_th, so it cancels exactly in the LF difference N_tot,ELG(M_h|>L1) - N_tot,ELG(M_h|>L2). Consequently the LF likelihood in Eq. (38) provides no constraint on F_Gauss, L_Gauss, or sigma_logM_h; those parameters are constrained only by the ACF (and the abundance through f_ELG). The text should state this limitation explicitly and discuss what the LF actually constrains: the error-function transition, the massive-end decay phi_d(M_h), and the satellite term. This is not fatal because the Gaussian component turns out to be small, but it is essential for the paper's claim that luminosity information makes the central occupation more robust.","section":"§3.3.3, Eq. (21)"},{"comment":"The central claim that joint ACF+LF fitting provides 'more robust constraints' than ACF+number-density fitting is not demonstrated by any ablation. The comparison with the Geach et al. (2012) model changes both the HOD parameterisation and the data entering the likelihood, so it cannot isolate the effect of the LF. I request an ablation in which the proposed HOD model is fitted (a) to the ACF plus the observed number density with a Gaussian penalty analogous to Eq. (33), and (b) to the ACF plus the LF, with a comparison of posterior widths and pulls for M_LF^min, f_sat, and b_g. In addition, the statement that the best-fitting LFs are 'fully consistent with the observed LFs' is not an independent check: the same observed LF data enter Eq. (38) and are generated by Eq. (21) from the fitted LHMR. A posterior-predictive or cross-validation test, such as predicting the NB921 LF from the NB816 fit or vice versa, would materially support the claim that the LF constraint adds information.","section":"§4.3.2 and §4.4"},{"comment":"The Milky-Way-descendant conclusion follows from M_LF^min (Eq. 39) evolved with EPS. M_LF^min is derived from the central occupation built on the deterministic double power-law LHMR of Eq. (11), with no explicit scatter and a functional form that is not validated independently. The paper itself notes in Section 5.2.1 that the bright/massive end of the LHMR is statistically difficult to constrain. A different LHMR shape, a redshift-dependent LHMR, or intrinsic scatter degenerate with sigma_logL could shift the inferred M_LF^min and therefore the z=0 descendant masses. Please test the sensitivity of M_LF^min and the Milky-Way connection to (i) an alternative LHMR parametrisation and (ii) scatter about the LHMR, or compare the predicted central occupation against hydrodynamical ELG catalogues (e.g., IllustrisTNG or MillenniumTNG) to bound the systematic error.","section":"§5.1.1, Fig. 11"}],"minor_comments":[{"comment":"The text states minimum reduced chi2 = 0.79 (NB816) and 0.86 (NB921), whereas Table 3 reports best-fit chi2/dof = 6.93/15 = 0.46 and 12.50/14 = 0.89 for the model with Gaussian central occupation. Please reconcile these values.","section":"§4.3.2 vs Table 3"},{"comment":"The text says the two power-law slopes satisfy 0<β<γ, but the no-Gaussian best fit for NB816 in Table 3 has β=1.03 and γ=1.03. Please clarify whether this inequality is a prior assumption actually imposed in the MCMC or an expectation stated after the fit.","section":"§3.1, Eq. (11)"},{"comment":"The variation of the fixed ELG fraction endpoints (f_min_ELG and f_max_ELG) is described as having little impact, but no results are shown. A short table or appendix giving the tested range and the resulting changes in M_LF^min, f_sat, and b_g would strengthen this claim.","section":"§4.3.1"},{"comment":"The y-axis range of the LHMR plot spans 38-52 in log10(L_ELG/erg/s), which is much wider than the data and literature points; narrowing the axis would improve readability.","section":"Figure 13"},{"comment":"The abstract contains the phrase 'mock catalogues of for future surveys'; please correct this typo.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this paper if you work on ELG clustering or mock catalogues. The genuinely new thing is putting a double-power-law luminosity-to-halo-mass relation directly into the central halo occupation so the model predicts luminosity functions and can be fit jointly to ACFs and LFs. That is a real step beyond Geach et al. (2012) and beyond Schechter-based CLF models.\n\nWhat the paper does well: the fits are statistically acceptable (reduced chi2 around 0.8 for both redshift slices), and the treatment of covariance, integral constraint, and contamination is careful. The no-Gaussian robustness test is good practice. The inferred M_LF_min ~ 11.6 and satellite fractions 0.2-0.3 are within the range of other ELG studies, and the LHMR comparisons to literature look sensible.\n\nThe main weakness is that the LF is used both as a fitting constraint and as a reported output, and the Gaussian central term cancels exactly in the LF difference (Eq. 21). So the LF constraints only shape the erfc transition and the satellite term; they do not constrain the Gaussian component. The double-power-law form of f_LHMR is assumed, not independently validated, and the bright end is weakly constrained—the paper admits this in Section 5.2.1. More importantly, the claim that LF constraints make the HOD more robust is never tested. They compare to the Geach model, but they do not run their own model with only ACF+number density to show what the LF actually adds. That is a fixable omission. No code is shipped, though public packages and public data make reproduction feasible.\n\nThese are moderate concerns, not fatal. The paper is honest about its limitations, and the resulting occupation picture—smooth rise around 10^11.5-12 h^-1 M_sun, small Gaussian component—is a reasonable outcome. The Milky Way descendant claim is speculative but clearly framed as such.\n\nWho gets value: people building ELG mocks for DESI, PFS, Euclid, and anyone interpreting [O II] clustering at z~1.2-1.5. It deserves a serious referee. I would send it to review with a request for an ablation (fit with and without LF) and a sharper statement about which parameters the LF actually constrains.","headline":"A useful HOD extension for [O II] emitters that joins ACFs and luminosity functions, but the added value of the LF constraint is asserted more than demonstrated.","tokens_in":45614,"tokens_out":2539,"would_cite":true,"duration_ms":24986,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper establishes that a halo occupation distribution model embedding a double-power-law luminosity-to-halo-mass relation can jointly reproduce the angular clustering and luminosity functions of [O II] emission-line galaxies, giving…","keywords":["halo occupation distribution","emission-line galaxies","[O II] emitters","luminosity function","angular correlation function","galaxy–halo connection","dark matter halos"],"falsifier":"One could measure the [O II] luminosities of central galaxies in halos of known mass (e.g., via galaxy groups or weak lensing) at z ~ 1.2–1.5 and check whether the relation follows a double power law; a plateau, strong redshift dependence, or significant scatter would falsify the model. Alternatively, a deeper survey constraining the bright end of the [O II] luminosity function at these redshifts would test the extrapolation: a bright-end excess or deficit relative to the double-power-law prediction would rule out the assumed LHMR shape.","tokens_in":2082,"feed_emoji":"🌌","tokens_out":1910,"duration_ms":102564,"temperature":0.7,"pith_summary":"This paper argues that the connection between [O II] emission-line galaxies (ELGs) and their host dark matter halos can be captured more completely by folding galaxy luminosity directly into the halo occupation distribution (HOD). The authors introduce a double-power-law luminosity-to-halo-mass relation (LHMR) that assigns each halo a predicted line luminosity, which lets the model predict luminosity functions as well as angular correlation functions. Applied to [O II] emitters at z = 1.193 and z = 1.471 from narrow-band imaging, the joint fit reproduces both measured angular clustering (reduced chi-squared around 0.8) and the observed luminosity functions. If correct, this yields a smoother, more gradual occupation of dark halos than earlier ELG models, characteristic halo masses around $10^{11}$.6 $h^{-1}$ solar masses, and satellite fractions of about 0.2–0.3. The authors further use the inferred halo masses to argue that the z ~ 1.4 emitters are plausible progenitors of present-day Milky-Way-like halos.","feed_headline":"Double power law ties [O II] emitters to 10^11.6 solar-mass halos","feed_subtitle":"A double-power-law luminosity-to-halo relation lets one model reproduce both clustering and galaxy luminosity data.","key_machinery":"The load-bearing object is the luminosity-to-halo-mass relation (LHMR), a double power law f_LHMR(M_h) = (L_0/2)[(M_h/M_trans)^β + (M_h/M_trans)^γ] (Eq. 11), which converts a halo mass into the predicted [O II] line luminosity of its central galaxy. This relation is embedded in the central occupation N_c(M_h | > L_th) (Eq. 10) through a Gaussian term plus an error-function term, and its inverse sets the satellite mass scale (Eq. 14). Inserting the occupation into Eq. 21 yields the model luminosity function, so the likelihood (Eq. 36) includes both the angular correlation function and the LF. The double-power-law shape, with positive slopes 0 < β < γ, is motivated by the stellar-to-halo-mass relation and observed ELG luminosity–halo-mass correlations; it is what allows the model to predict the LF without prescribing a Schechter or other analytic form.","core_discovery":"The central claim is that a deterministic double-power-law relation between [O II] line luminosity and halo mass, f_LHMR(M_h) = (L_0/2)[(M_h/M_trans)^β + (M_h/M_trans)^γ], when inserted into the HOD central occupation function, gives the model enough information to predict the differential luminosity function without assuming a functional form for it. Fitting jointly to the angular correlation function and the luminosity function of [O II] emitters at two redshifts (z = 1.193 and z = 1.471) yields reduced chi-squared values of 0.79 and 0.86, with the best-fitting luminosity functions fully consistent with observations. The inferred median halo masses of central ELGs are log10(M_LF^min / $h^{-1}$ M_sun) = 11.60 (+0.19, -0.20) and 11.66 (+0.18, -0.19) for the two samples, and satellite fractions are 0.31 ± 0.08 and 0.20 ± 0.06. The paper shows that the Gaussian component of the central occupation is negligible, so a simplified model without it fits equally well, and it interprets the smooth occupation as a more realistic representation of ELG environments compared with a sharp mass threshold.","pith_inferences":["The framework is naturally extendable to other line tracers (Hα, Lyα) or broad-band-selected populations: one only needs to replace the LHMR with the appropriate luminosity–halo-mass calibration, and the LF constraint becomes available.","A direct test of the deterministic LHMR is to compare the model's predicted zero scatter in central luminosity with measurements from group or weak-lensing samples; adding a lognormal scatter is a straightforward generalization.","The joint ACF+LF likelihood may break degeneracies that clustering alone leaves unresolved, especially between the satellite occupation amplitude and the central–satellite boundary, which could sharpen forecasts for ELG redshift-space distortion analyses.","The model ignores assembly bias by construction; hydrodynamical or high-resolution N-body simulations that introduce formation-time dependence could shift the inferred M_LF^min and satellite fractions, so quantifying that sensitivity would be a useful next step."],"forward_implications":["If the model is right, ELG clustering analyses gain a new constraint: the luminosity function carries differential number-density information that the integrated number density alone does not, making halo-mass inferences less degenerate.","The inferred occupation is smooth rather than sharply thresholded, implying ELGs populate halos across a broad mass range rather than switching on abruptly at a single mass.","The characteristic halo masses of about 10^11.6 h^-1 M_sun and satellite fractions of 0.2–0.3 at z ~ 1.2–1.5 provide calibration targets for ELG mock catalogues for upcoming redshift surveys.","Extrapolating the inferred halo masses forward with the extended Press–Schechter formalism suggests the z ~ 1.5 [O II] emitters evolve into Milky-Way-sized halos by z = 0.","The measured satellite fraction at z ~ 1.47 is more than 1σ below the IllustrisTNG simulation prediction, hinting that the simulation may overproduce star-forming satellites at that epoch."],"supporting_citations":[{"why":"Provides the baseline ELG HOD model with Gaussian and error-function central occupations that this paper extends and compares against.","marker":"Geach et al. (2012)"},{"why":"Supplies the [O II] ELG catalogues and the observed luminosity functions used as constraints in the joint fit.","marker":"Hayashi et al. (2020)"},{"why":"Previous HOD analysis of the same [O II] samples with the Geach model, giving the comparison values for bias, satellite fraction, and fake fraction.","marker":"Okumura et al. (2021)"},{"why":"Observational evidence for a halo mass–line luminosity correlation that motivates the double-power-law LHMR shape.","marker":"Khostovan et al. (2019)"},{"why":"The standard HOD central occupation with an error-function transition, which the error-function term in Eq. 10 adopts.","marker":"Zheng et al. (2005)"},{"why":"Simulation-based ELG HOD predictions from IllustrisTNG used for the satellite-fraction comparison.","marker":"Osato & Okumura (2023)"},{"why":"Provides the empirical relation between M_cut and M_sat that reduces the satellite occupation parameters.","marker":"Conroy et al. (2006)"}],"fun_headline_variants":["Double power law links [O II] emitters to 10^11.6 solar-mass halos","New HOD model uses luminosity to pin down [O II] galaxy halos","Luminosity-aware HOD reveals [O II] emitters seed Milky Way-like halos","HOD with double-power-law luminosity fits [O II] clustering and LF","Subaru HSC finds [O II] emitters occupy ~10^11.6 solar-mass halos"],"cache_read_input_tokens":47616,"weakest_assumption_plain":"The load-bearing premise is that the true relation between [O II] line luminosity and halo mass is exactly the deterministic double power law of Eq. 11, with no scatter and no assembly bias, and that its parameters can be learned from the same luminosity function data used in the likelihood.","fun_headline_variants_meta":{"raw":{"variants":["Double power law links [O II] emitters to 10^11.6 solar-mass halos","New HOD model uses luminosity to pin down [O II] galaxy halos","Luminosity-aware HOD reveals [O II] emitters seed Milky Way-like halos","HOD with double-power-law luminosity fits [O II] clustering and LF","Subaru HSC finds [O II] emitters occupy ~10^11.6 solar-mass halos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000737,"raw_usage":{"total_tokens":3410,"prompt_tokens":1176,"completion_tokens":2234,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":792,"completion_tokens_details":{"reasoning_tokens":2115}},"tokens_in":792,"tokens_out":2234,"duration_ms":14059,"temperature":1.0,"reasoning_tokens":2115,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T23:49:09.766753+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One could measure the [O II] luminosities of central galaxies in halos of known mass (e.g., via galaxy groups or weak lensing) at z ~ 1.2–1.5 and check whether the relation follows a double power law; a plateau, strong redshift dependence, or significant scatter would falsify the model. Alternatively, a deeper survey constraining the bright end of the [O II] luminosity function at these redshifts would test the extrapolation: a bright-end excess or deficit relative to the double-power-law prediction would rule out the assumed LHMR shape.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the [O II] ELG catalogues and the observed luminosity functions used as constraints in the joint fit."}],"review_version":1}