{"id":"31c1c24d-daeb-42c5-977b-fbdfb83515ea","arxiv_id":"2501.05575","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The distribution of Lyman-alpha optical depths in 71 XQ-100 quasar spectra is consistent with no helium-reionization temperature fluctuations, giving sigma(ln T) < 0.29 at 2 sigma on 100 Mpc scales at z=3.76.","lead":"This paper uses 71 quasar spectra to search for large-scale temperature variations in the gas between galaxies at redshifts 3.7 to 4.2, the epoch when helium reionization should heat the universe unevenly. It finds no such variations and rules out temperature contrasts larger than about 30 percent at 100 megaparsec scales, the tightest limit so far.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline limit assumes full 100 Mpc coherence of temperature fluctuations; if the real coherence length is shorter, the same optical-depth data allow larger local temperature contrasts, so the abstract's '<30%' interpretation is not yet robust.","rationale":"The reader's conditional verdict identifies the same load-bearing assumption. The paper is internally consistent: the numerical pipeline is careful (PCA continuum errors, DLA masking, noise injection, mean-flux calibration), and the 0.29 limit is a valid statement about the coherent-injection model. The soft spot is external validity: the statistic used (tau_eff scatter) is sensitive to the amplitude of sightline-mean temperature offsets, and the model converts that observable into a physical temperature contrast by assuming each sightline is entirely hot or cold. The theory prediction being tested is the contrast between recently-ionized and neutral regions, which need not be coherent over a full 100 Mpc bin. McQuinn et al. 2009 and later simulations place the coherence scale at tens of Mpc, and the paper's own 50 Mpc limits are weaker (Table 3: <0.40 vs <0.29 at z=3.76). A simple decomposition of a 100 Mpc sightline into independent chunks predicts the allowed per-patch sigma scales as sqrt(N), so local contrasts of 30-50% would escape the headline limit for N approximately 3-4. The proposed N-chunk rerun settles this directly with the same code and data; if the scaling is much weaker than sqrt(N), the coherence concern would be largely retired. The paper is transparent about the 'rather simplistic' model (Sec. 3.5), and no ad hominem is warranted. The appropriate verdict remains conditional: accept the measurement as a constraint on fully coherent 100 Mpc-scale temperature fluctuations, but do not yet read it as a general bound on helium-reionization temperature contrasts.","tokens_in":19364,"tokens_out":21286,"duration_ms":212835,"concrete_test":"Rerun the z=3.76, L=100 Mpc likelihood of Sec. 3.5.1 with finite-coherence injection: divide each simulated 100 Mpc sightline into N equal subsegments (N=2, 5, 10), draw independent Delta ln T ~ N(0, sigma^2) per subsegment, convert locally with Delta ln tau = -0.352 Delta ln T, compute the full-segment tau_eff, and repeat the full forward-modeling and likelihood (including mean-flux calibration) to find the 2 sigma upper limit on sigma as a function of N. If the limit rises from 0.29 at N=1 toward roughly sqrt(N)*0.29 at N>=5, the headline bound constrains only the 100 Mpc-smoothed temperature and the abstract's '<30% local contrast' conclusion is model-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central constraint is derived in Sec. 3.5.1 by drawing a single Delta ln T per 100 Mpc sightline and applying it as a flat rescaling of the transmitted flux via Delta ln tau = -0.352 Delta ln T (Bolton et al. 2005). This is a fully coherent temperature fluctuation. The statistic actually measured is the sightline-to-sightline scatter of tau_eff, and a coherent offset is the most efficient way to generate that scatter. If real helium-reionization temperature structure has a shorter correlation length, a 100 Mpc sightline averages over many independent patches, so a given local temperature contrast produces a smaller scatter in tau_eff; the same observations would then permit a larger physical Delta T/T. The paper's own Table 3 shows the sensitivity: at z=3.76 the 2 sigma limit loosens from <0.29 (100 Mpc assumed coherence) to <0.40 (50 Mpc assumed coherence), and for shorter coherence the degradation continues roughly as sqrt(L/l_c). Section 5 concedes the temperature model is 'rather simplistic' and that density-correlated gamma changes would 'slightly reduce' the effect, but no finite-coherence or radiative-transfer run is used to calibrate the mapping. The abstract's conclusion that the IGM temperature contrast is 'less than ~30%' follows only under the unverified assumption that the fluctuations are coherent on 100 Mpc scales.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper introduces a new method to constrain large-scale IGM temperature fluctuations during helium reionization by measuring the distribution of effective Lyman-α optical depths toward 71 quasars from the XQ-100 survey at z ≈ 3.76–4.19. The observed distributions are compared to Nyx hydrodynamical simulations with no temperature fluctuations, and the simulations are then post-processed to include coherent, lognormal temperature offsets per sightline with amplitude σ(ln T). The central result is an upper limit σ(ln T) < 0.29 (0.40 at 3σ) at z = 3.76 for 100 comoving Mpc averaging, with weaker constraints at higher redshifts and at 50 Mpc. The authors interpret this as implying that helium reionization had not imprinted temperature contrasts larger than about 30% at z = 3.76, or that it had not yet significantly started.","tokens_in":19635,"tokens_out":13063,"duration_ms":111006,"significance":"The measurement is carefully executed: the PCA continuum reconstruction has quantified accuracy, the sample is restricted to a flux-calibration-safe wavelength range, DLAs and bad pixels are masked, bootstrap resampling is used for uncertainties, and the simulations are forward-modeled with instrumental noise and continuum errors. The no-fluctuation model matches the observed cumulative distributions within 2σ, and the new statistic (the distribution of effective optical depths) is shown to have constraining power comparable to forecasts from power-spectrum analyses. If the model-dependent interpretation is validated, this would be the tightest constraint to date on large-scale temperature fluctuations from helium reionization. However, the physical interpretation relies on the assumption that temperature fluctuations are coherent over the sightline length and that the adopted τeff ∝ T^0.352 mapping captures the relevant physics; the paper itself labels the model 'rather simplistic,' which limits the astrophysical conclusion.","major_comments":[{"comment":"The headline limit σ(ln T) < 0.29 is derived for a model in which each 100 Mpc sightline receives a single coherent lognormal temperature offset. If real helium reionization temperature fluctuations have a shorter coherence length, a given local temperature contrast produces less sightline-to-sightline scatter in τeff, so the same data allow larger physical temperature contrasts. The paper's own Table 3 shows the 2σ limit degrades from <0.29 at 100 Mpc to <0.40 at 50 Mpc, and the degradation continues for shorter coherence lengths approximately as sqrt(L/L_c). Section 5 concedes the fluctuations are 'unlikely to be so solidly coherent' but does not quantify the effect. The abstract's conclusion that the IGM temperature contrast is 'less than ~30%' is therefore not robust to finite coherence. Please either restrict the physical interpretation to the coherent model parameter, or calibrate the mapping using a more realistic reionization temperature field (e.g., a radiative-transfer simulation) and report the resulting limits on the local temperature contrast.","section":"§3.5.1, Table 3, abstract"},{"comment":"The description of the injection step is ambiguous: the text states that a single ΔlnT is drawn and converted to 'excess optical depth, Δlnτ = −0.352ΔlnT', and is then introduced 'via a flat rescaling of the transmitted flux of the entire sightline.' A constant multiplicative rescaling of the flux changes τeff by an additive constant, not by a multiplicative factor; the relation Δlnτeff = −0.352ΔlnT requires multiplying the optical depths by a constant (or applying a flux factor that depends on the sightline's τeff). As written, the implementation is not uniquely defined, so the injected amplitude—and hence the limits in Table 3—is not reproducible. Please provide the exact transformation and verify that the resulting effective optical depth scales as τeff ∝ T^0.352 as intended from Bolton et al. (2005) Eq. (4).","section":"§3.5.1 (Adding temperature fluctuations)"},{"comment":"The paper states that incorporating γ fluctuations correlated with T fluctuations would 'slightly reduce their impact,' but no calculation or simulation is shown to support this. Since the aim is to constrain physical temperature contrasts, the possible cancellation between temperature and temperature-density-relation changes is a source of systematic uncertainty that should be quantified (e.g., by applying a simple correlated T–γ model to the simulated sightlines) or explicitly folded into the reported limits. Without this, the robustness of the upper limits to realistic helium reionization heating is not established.","section":"§5 (Discussion)"}],"minor_comments":[{"comment":"There are several typographical errors, including 'sigthline' (§3.5.1), 'uncertainity' (§3.4), 'continum' (§3.1), and 'wavelenghts' (§3.2). These should be corrected before publication.","section":"Throughout"},{"comment":"The caption reads 'blue mean and1/2σ contours'; the meaning of '1/2σ' is unclear. It likely should be '1σ and 2σ' or '1–2σ' shaded contours.","section":"Figure 3 caption"},{"comment":"The kernel density estimation bandwidth is not specified; the resulting p-values may depend on the chosen bandwidth. Please state the bandwidth or the rule used to set it.","section":"§3.6"},{"comment":"The statement 'Our measurements are the only such constraints on temperature fluctuations from helium reionization thus far' is stronger than the cited literature supports; consider softening to 'the first direct constraints using this statistic' or adding a comparison to existing indirect constraints.","section":"§5"},{"comment":"Using the z = 4 snapshot for all redshift bins is an approximation; the effect of structure growth between z = 4.19 and z = 3.76 on the τeff distribution width is not quantified. Please add a sentence noting the expected magnitude of this systematic.","section":"§3.5"}],"recommendation":"major_revision","confidential_remarks":"The paper's own Sec. 5 acknowledges the 'rather simplistic' temperature model, yet the abstract presents the constraint as a physical limit on the IGM temperature contrast. This mismatch between the abstract and the model limitations is the main reason I am recommending major revision rather than minor. The measurement itself appears careful and the method is novel; if the authors address the coherence and implementation issues, the paper will be a solid contribution to the field."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nPunchline: This is a genuinely new null result—first direct constraints on large-scale IGM temperature fluctuations during helium reionization from the effective optical depth distribution—but the headline '<30%' interpretation is only as strong as the coherence assumption built into the injection model. The authors know this and say so. The measurement itself is careful and honest.\n\nWhat's new and good: They apply the tau_eff distribution statistic, previously used for hydrogen reionization at z~5-6, to He reionization at z~3.7-4.2 with 71 XQ-100 sightlines. The pipeline is solid: PCA continuum reconstruction with empirical uncertainties, DLA masking, and a real effort to handle a known X-Shooter flux-calibration problem by discarding the UV arm and re-defining bins. Forward modeling includes instrumental noise and continuum errors, and they calibrate the mean flux per bin. The no-temperature-fluctuation simulation matches the observed CDFs within 2σ at all redshifts. The resulting upper limit at z=3.76 on 100 Mpc scales, sigma(ln T)<0.29 at 2σ, is a real number under their model.\n\nSoft spots: The stress-test concern is correct. The model injects one Delta ln T per 100 Mpc sightline, applied as a flat rescaling of the entire transmitted flux. That is the most efficient way to produce sightline-to-sightline scatter in tau_eff. Real helium reionization presumably has shorter coherence lengths and density-correlated temperature changes, both of which would dilute the observable signal. Table 3 shows the limit loosens from <0.29 to <0.40 at 50 Mpc; if the true coherence is even shorter, the same data allow even larger physical temperature contrasts. The abstract's 'less than ~30%' thus follows only under the assumed 100 Mpc coherence. This is a model-dependence concern, not a fatal flaw—the paper concedes as much in Sec. 3.5 and 5. I'd like to see a more realistic helium reionization simulation (or at least a coherence-length scan) before accepting the physical interpretation at face value.\n\nOther minor issues: they use a single Nyx snapshot at z=4 rescaled to the redshifts of the bins, and there's no systematic error budget on the UV background or the Bolton et al. scaling. The p-value approach is reasonable for a null measurement, though the 'tentative detection' at z>4 is appropriately not claimed.\n\nWho it's for: IGM and reionization folk. It's a useful data paper and a method paper, and it deserves a serious referee.\n\nRecommendation: send to peer review. The referees should push on the coherence-length question and ask for the constraints to be re-expressed in a more model-robust form or demonstrated to be insensitive to the assumed coherence. I would not desk-reject this.","headline":"A careful new null measurement of IGM temperature fluctuations during helium reionization, but the headline upper limit depends on an unverified 100 Mpc coherence assumption that the authors themselves flag.","tokens_in":20182,"tokens_out":4038,"would_cite":true,"duration_ms":35726,"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":"Quasar spectra cap IGM temperature swings at 100 Mpc scales to below 30 percent","keywords":["intergalactic medium","helium reionization","Lyman-alpha forest","effective optical depth","temperature fluctuations","XQ-100 survey","quasar spectra"],"falsifier":"A measurement that would settle the claim is a comparison at z=3.76 on 100 Mpc scales using several hundred independent quasar sightlines with known continuum uncertainties: if the observed tau_eff distribution remains consistent with the no-fluctuation simulation, the sigma(ln T) < 0.29 limit is confirmed; if the distribution shows excess width beyond the density-only prediction, the limit is falsified. A second decisive test is a direct simulation of helium reionization with full radiative transfer that predicts a coherent temperature contrast above 30 percent at z=3.76; if such a model reproduces the observed narrow optical depth distribution, then the paper's temperature-opacity mapping is wrong rather than the temperature contrast being small.","tokens_in":88,"feed_emoji":"🔭","tokens_out":1681,"duration_ms":61188,"temperature":0.7,"pith_summary":"The last great heating of intergalactic gas, helium reionization, should leave hot and cold patches of gas side by side, with temperature contrasts of tens of percent. This paper argues that those contrasts would widen the spread of Lyman-alpha forest optical depths seen toward background quasars. Using 71 high-quality X-Shooter spectra from the XQ-100 survey, they measure that spread in four redshift bins from z=3.76 to z=4.19 and compare it to a hydrodynamical simulation that contains no temperature fluctuations. The observed distribution matches the no-fluctuation simulation, so no helium-reionization heating signature is detected, and by injecting artificial fluctuations until the model disagrees they set upper limits, the tightest so far, on the allowed temperature contrast.","feed_headline":"Quasar spectra cap IGM temperature swings below 30 percent","feed_subtitle":"The spread of Lyman-alpha absorption across 71 sightlines rules out large helium-reionization heating at z=3.76.","key_machinery":"The central object is the distribution of effective Lyman-$\\alpha$ optical depths, tau_eff, measured in 100 Mpc (and 50 Mpc) bins along 71 quasar sightlines. The mechanism is the temperature-opacity relation: hotter gas recombines more slowly, so a higher IGM temperature lowers the Lyman-$\\alpha$ effective optical depth, and coherent large-scale temperature fluctuations rescale whole sightlines, broadening the observed tau_eff distribution beyond what the density field alone produces. The argument is carried by a lognormal temperature-fluctuation model in which each simulated sightline is rescaled by $\\Delta$ ln tau = -0.352 $\\Delta$ ln T, using the Bolton et al. (2005) calibration tau_eff proportional to $T^{0}$.352, followed by forward modeling that adds instrumental noise, continuum uncertainty, and a mean-flux calibration.","core_discovery":"At redshift z=3.76, on 100 comoving Mpc scales, the rms large-scale temperature fluctuation of the intergalactic medium is constrained to $\\sigma$(ln T) < 0.29 at 2 $\\sigma$ (and < 0.40 at 3 $\\sigma$), corresponding to a temperature contrast between ionized and neutral regions of roughly $\\Delta$ T / T ~ 34% (49%). The observed effective optical depth distribution is consistent within 2 $\\sigma$ with a cosmological hydrodynamical simulation that includes no temperature fluctuations, and the limits weaken at higher redshifts (z=3.90: <0.32; z=4.04: <0.74; z=4.19: <0.64 at 2 $\\sigma$) because fewer sightlines are available. The paper concludes that either helium reionization had not yet imprinted temperature contrasts larger than roughly 30 percent at z=3.76, or the process had not significantly started at that redshift.","pith_inferences":["The quoted upper limit on sigma(ln T) is only as strong as the coherence assumption; real helium reionization produces density-dependent, spatially structured temperature fields, so the limit should be read as constraining the coherent part of the temperature contrast rather than the total physical contrast.","Because the same statistic is sensitive to any large-scale opacity fluctuation, the null result can be read as a consistency check of the standard cosmological model and of the UV background model used in the simulations.","A testable extension would repeat the measurement on the much larger, lower signal-to-noise samples from DESI, WEAVE-QSO, and 4MOST, where the reduced continuum precision and increased noise are the main challenges to overcome.","The scaling estimate sigma(ln T) ~ 1/sqrt(N_qso) suggests that pushing the same method to higher redshift with more sightlines could distinguish between a late-start and a low-contrast helium reionization history."],"forward_implications":["If the central claim holds, helium reionization's large-scale temperature contrast at z=3.76 is at most about 30 percent, ruling out the most extreme heating scenarios at that epoch.","The effective optical depth distribution is a competitive and much simpler probe than the Lyman-alpha forest power spectrum for constraining reionization-induced temperature fluctuations.","Roughly 400 to 700 quasar sightlines of similar quality would be needed to reach 2 sigma sensitivity to sigma(ln T) = 0.1, which would probe the temperature fluctuations predicted during most of helium reionization.","The same method would be sensitive to any physical process that adds large-scale fluctuations to quantities modifying Lyman-alpha opacity, such as UV background fluctuations or the matter distribution.","The tentative preference for non-zero temperature fluctuations at z=4.04 and z=4.19, if confirmed with larger samples, could indicate the onset of helium reionization at z > 4."],"supporting_citations":[{"why":"Supplies the calibration tau_eff proportional to T^0.352 used to convert temperature fluctuations into optical depth fluctuations.","marker":"Bolton et al. 2005"},{"why":"Provides the XQ-100 survey, the 71 X-Shooter quasar spectra that constitute the observed optical depth distributions.","marker":"López et al. 2016"},{"why":"Supplies the Nyx cosmological hydrodynamical code used to run the baseline simulation of density-driven Lyman-alpha opacity fluctuations.","marker":"Almgren et al. 2013"},{"why":"Defines the simulation setup, skewer extraction, and Lyman-alpha opacity computation that produces the no-temperature-fluctuation model.","marker":"Lukić et al. 2015"},{"why":"Provides the PCA continuum reconstruction method used to predict each quasar's intrinsic Lyman-alpha forest continuum, with empirically calibrated uncertainties.","marker":"Bosman et al. 2021"},{"why":"Provides the literature measurements of mean effective optical depth against which the XQ-100 measurements are validated.","marker":"Becker et al. 2013"},{"why":"Supplies the radiative transfer prediction of peak temperature fluctuations Delta T/T ~ 0.2 on 50-150 Mpc scales, the benchmark that the new constraints approach or rule out.","marker":"McQuinn et al. 2009"}],"fun_headline_variants":["IGM temperature swings held below 30% by quasar spectra","Lyman-alpha forest narrows IGM temperature fluctuations to under 30%","Helium reionization leaves no large temperature imprints in IGM","On 100 Mpc scales, IGM temperature contrast capped at 0.29","No sign of helium reionization heating in IGM at z=3.76"],"cache_read_input_tokens":22272,"weakest_assumption_plain":"The argument assumes that a large-scale temperature fluctuation acts as a coherent, sightline-wide rescaling of the effective optical depth with tau_eff proportional to $T^{0}$.352, whereas real helium reionization produces temperature fluctuations that are density-dependent, spatially structured, and also alter the temperature-density relation.","fun_headline_variants_meta":{"raw":{"variants":["IGM temperature swings held below 30% by quasar spectra","Lyman-alpha forest narrows IGM temperature fluctuations to under 30%","Helium reionization leaves no large temperature imprints in IGM","On 100 Mpc scales, IGM temperature contrast capped at 0.29","No sign of helium reionization heating in IGM at z=3.76"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000784,"raw_usage":{"total_tokens":3552,"prompt_tokens":1128,"completion_tokens":2424,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":744,"completion_tokens_details":{"reasoning_tokens":2336}},"tokens_in":744,"tokens_out":2424,"duration_ms":16322,"temperature":1.0,"reasoning_tokens":2336,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:13:06.320993+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement that would settle the claim is a comparison at z=3.76 on 100 Mpc scales using several hundred independent quasar sightlines with known continuum uncertainties: if the observed tau_eff distribution remains consistent with the no-fluctuation simulation, the sigma(ln T) < 0.29 limit is confirmed; if the distribution shows excess width beyond the density-only prediction, the limit is falsified. A second decisive test is a direct simulation of helium reionization with full radiative transfer that predicts a coherent temperature contrast above 30 percent at z=3.76; if such a model reproduces the observed narrow optical depth distribution, then the paper's temperature-opacity mapping is wrong rather than the temperature contrast being small.","supporting_citations":[],"review_version":1}