{"id":"87fcef00-5cee-4b66-a119-cf20f377eba1","arxiv_id":"1908.01522","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"IRTS data reveal a theta^-1 power-law decline in near-infrared background fluctuations from 2 to 20 degrees, implying a 1-degree bump when joined with sub-degree CIBER measurements.","lead":"This paper reports the first measurement of near-infrared extragalactic background light fluctuations at angular scales from 2 to 20 degrees, using 1995 IRTS satellite data. The fluctuations decline as a smooth power law, and connecting with smaller-scale data suggests a broad bump around 1 degree that known sources cannot explain.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Zodiacal-light residual after the single-slope Kelsall fit is not directly ruled out at 2-20 deg and could fake the reported power law.","rationale":"The reader's weakest_assumption correctly identifies the Kelsall-model ZL subtraction as the most dangerous step. I agree that the single fitted slope cannot validate the model's spatial structure at the scales being measured. Although the paper also infers the 1-degree bump by joining non-overlapping IRTS and CIBER datasets, the text uses hedged language ('connect well', 'appear to have'), and the genuinely new measurement is the 2-20 deg spectrum; if the ZL residual is significant, even that measurement fails. The paper has real strengths: the absolute NIREBL brightness from the y-intercept is consistent with Matsumoto et al. (2015), the error budget is detailed, and the analysis uses publicly available tools (PolSpice, POKER). However, none of these independently validates the spatial smoothness of the ZL residual at 2-20 deg. The proposed independent-ZL-model test is a standard robustness check that would settle the concern. Therefore the reader's CONDITIONAL verdict remains appropriate; no change is needed.","tokens_in":17125,"tokens_out":7377,"duration_ms":77786,"concrete_test":"Recompute the 1.6 and 2.2 um fluctuation spectra with the Kelsall et al. (1998) ZL model replaced by an independent zodiacal-light model (e.g., Wright 1998 or a DIRBE-based empirical ZL map), keeping the same per-band slope-fitting, clipping, mask, and PolSpice pipeline. If the 2-20 deg power amplitudes change by more than the quoted total errors, the reported NIREBL fluctuations are not robust to ZL modelling; if they are unchanged, the ZL-residual concern is empirically retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing assumption is that the Kelsall et al. (1998) zodiacal-light model, after a single per-band multiplicative correction fitted to the same IRTS data (slopes 0.811 and 1.146, Section 4.2.3), is free of angular structure at 2-20 deg at the required precision. ZL is more than 90% of the raw sky brightness (Section 4.1), so a residual of order 0.5-1% of ZL has amplitude comparable to the reported NIREBL fluctuations. The single slope can fix overall normalization and band mismatch, but it cannot remove a spatially varying error in the model's ecliptic-latitude gradient, in its seasonal/time-dependent terms, or in the 1.6 um color extrapolated from 1.25 um assuming a fixed spectral shape. The cited large-scale smoothness (Abraham et al. 1997; Pyo et al. 2012) does not quantify the residual after this specific fit, and the paper does not measure ZL fluctuation power at 2-20 deg (Section 8). Because the slope is fitted to the same data that produce the final map, the subtraction absorbs only the component correlated with the model; an uncorrelated ZL residual remains. Such a residual could directly create a theta^-1-like power spectrum and would shift the claimed connection to the 1-degree bump.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports measurements of the spatial fluctuations of the near-infrared extragalactic background light (NIREBL) at 1.6 and 2.2 μm using IRTS/NIRS data, covering angular scales from 2° to 20°. After subtracting foregrounds (zodiacal light, diffuse Galactic light, and integrated star light), the authors estimate fluctuation power spectra with the PolSpice algorithm, subtract noise via Monte Carlo simulations, and derive an error budget including statistical and systematic terms. They find a monotonic decline of F(√(l(l+1)Cl/2π)) roughly as θ^{-1} over 2°–20°, and, by connecting to CIBER sub-degree measurements, they infer a broad bump centered at about 1°. They examine candidate sources—normal galaxies, intra-halo light, DGL, and high-redshift objects—and conclude that none can explain the excess around 1°.","tokens_in":17389,"tokens_out":6167,"duration_ms":59385,"significance":"If the measurement is robust, this is the first determination of NIREBL fluctuations at degree scales, bridging the gap between sub-degree space-based measurements and the large-scale regime. The paper has clear strengths: it uses a partial-sky power-spectrum estimator (PolSpice), performs Monte Carlo noise subtraction, presents a detailed error budget, and checks consistency with DIRBE and with the previously measured absolute NIREBL brightness. The connection to CIBER data, if taken as an inference rather than a direct measurement, offers a plausible synthesis of the current data. However, the central claims rest on the adequacy of the zodiacal-light subtraction and on the interpretation of the 1° bump as an inferred feature; both points require further support before the conclusions can be accepted as quantitative results.","major_comments":[{"comment":"The zodiacal-light subtraction is the load-bearing step in the measurement, and the paper does not demonstrate that the residual ZL after the per-band slope correction is free of angular fluctuations at 2°–20°. ZL is more than 90% of the raw sky brightness (§4.1); a residual of order 0.5–1% of ZL has amplitude comparable to the reported NIREBL fluctuations. The correction factors 0.811 and 1.146 are fitted to the same IRTS data in §4.2.3, so they can only absorb the component of the observed sky that is proportional to the Kelsall model. A spatially varying mismatch in the ecliptic-latitude dependence, in the seasonal terms, or in the 1.25→1.6 μm color extrapolation would remain in the residual. Section 8 states that \"since the ZL is based only on the model, we do not evaluate the ZL fluctuation in this work,\" and the error budget in Table 1 contains no ZL subtraction uncertainty. I request a quantitative test, for example recomputing C_l after fitting the ZL normalization in independent angular patches, after masking the highest-ZL-gradient regions, or after substituting an independent zodiacal-light model, and reporting the spread of the resulting power spectra. Without such a test, the reported θ^{-1} power law could be produced by a residual ZL structure.","section":"§4.2.3, §8"},{"comment":"The claim of a \"wide bump with a center at around 1°\" is an inference, not a direct measurement from the IRTS data. The IRTS spectra cover only 2°–20°; the bump is obtained by connecting them to the CIBER sub-degree measurements across a gap with no data. The abstract and summary state this as an observational result, which overstates the evidence. In addition, the 2.2 μm CIBER spectrum shown in Fig. 11 is constructed by scaling the CIBER 1.6 μm spectrum using the IRTS 1.6/2.2 color ratio from Fig. 12, so part of the apparent agreement at 2.2 μm is built into the comparison. The text should explicitly label the bump as an inference from comparison with other experiments, and the location of the peak should be given with the caveat that it is not constrained by the IRTS data themselves.","section":"Abstract, §7, §9"},{"comment":"The statement that the fluctuation spectrum follows \"nearly a single power-law spectrum (i.e. F(√(l(l+1)Cl/2π)) ~ θ^{-1})\" is made without a quantitative fit. No best-fit power-law index, its uncertainty, or a goodness-of-fit statistic is reported for either band. Because the power-law description is part of the central claim, the authors should fit the data over the 2°–20° range and report the index and its errors, including the effect of the systematic error budget.","section":"§7, Figs. 10–11"}],"minor_comments":[{"comment":"In the sentence beginning \"Second, their exists systematic uncertainty,\" \"their\" should be \"there.\"","section":"§4.2.3"},{"comment":"The word \"analized\" should be \"analyzed.\"","section":"§8"},{"comment":"The phrase \"Laplace's spherical harmonics\" is awkward; consider simply \"spherical harmonics\" or \"Laplace's spherical harmonic functions.\"","section":"§5, Eq. (3)"},{"comment":"The caption states that the ZL fluctuation is not shown because the ZL is model-based and expected to have very small fluctuation; this is precisely the assumption that needs to be tested, so the caption should refer the reader to the analysis requested in the major comment on ZL subtraction.","section":"Fig. 13 caption"},{"comment":"In \"The final data at 24 discrete bands covers 1% of whole sky,\" \"covers\" should be \"cover\" and \"whole sky\" should be \"the whole sky.\"","section":"§3"},{"comment":"The Knox formula is described as \"empirically determined\"; it is an analytic approximation for the sample variance of C_l. Please describe it more precisely.","section":"§6"},{"comment":"The absolute NIREBL brightness values 56.032 and 28.228 nW m^{-2} sr^{-1} are quoted with more significant figures than the systematic uncertainties justify; rounding to at most one decimal place would be more appropriate.","section":"§7"}],"recommendation":"major_revision","confidential_remarks":"The main risk to the paper is the zodiacal-light subtraction: the ZL model is corrected with a single per-band factor fitted to the same data, and the paper explicitly declines to evaluate ZL fluctuations in Section 8. This is a load-bearing gap. The claim of a 1-degree bump is also overstated in the abstract, since it is an inference across a data gap. Both issues are addressable with additional analysis and rewording, so I recommend major revision rather than rejection. The editor may also wish to note that the absence of any ZL-related contribution in Table 1 should be explicitly justified or corrected."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is genuinely the first measurement of NIR EBL fluctuations at 2-20 degrees at 1.6 and 2.2 um, using IRTS archival data. That alone makes it worth a serious look. The power-spectrum pipeline is careful — PolSpice for partial sky, Monte Carlo noise subtraction, a real error budget, and direct checking of DGL with AKARI NEP data. The result at 2-20 deg is a clean power-law with F ~ theta^-1. That part holds up.\n\nThe soft spot is the rest of the storyline. The 'bump at 1 degree' is not measured by IRTS; the data begin at 2 degrees. It is inferred by connecting to CIBER sub-degree points. The authors are fairly honest about this — they say 'infer' and 'appear to have' — but the abstract and summary state it as a conclusion. A reader needs to see the gap.\n\nThe bigger concern is zodiacal light. ZL is over 90% of raw sky brightness. The Kelsall model is scaled by a single slope per band fitted to the same IRTS data (0.811 and 1.146). That absorbs normalization but not spatially varying model error. The paper explicitly says they do not evaluate ZL fluctuation at 2-20 deg. A residual of ~0.5-1% of ZL would be comparable to the reported signal and could create a theta^-1-like spectrum. This is not a fatal flaw — previous work supports large-scale ZL smoothness and the authors cite it — but the burden is on them to show the residual is small. They don't.\n\nAlso, the 1.6/2.2 color ratio used to scale CIBER to 2.2 um comes from the same maps, so the comparison at 2.2 um is not independent. Minor, but worth noting.\n\nNet: the 2-20 degree measurement is a real contribution, and the ZL issue is a known hard problem in this field, not a careless error. The paper deserves peer review and likely conditional acceptance after the authors either quantify the ZL residual or soften the claim about the 1-degree bump.","headline":"First measurement of near-IR EBL fluctuations at 2-20 degrees, but the 1-degree bump is an inference across a gap and the zodiacal-light subtraction has an unquantified residual.","tokens_in":17999,"tokens_out":2005,"would_cite":true,"duration_ms":19407,"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 near-infrared extragalactic background shows a broad, unexplained fluctuation bump centered near 1 degree, discovered by the first measurement at 2–20 degree scales.","keywords":["near-infrared extragalactic background light","NIREBL","fluctuation power spectrum","IRTS","zodiacal light","cosmic infrared background","degree-scale anisotropies","intra-halo light"],"falsifier":"Re-observe the same IRTS sky fields at a different time of year (or use a different zodiacal-light model) and recompute the 1.6 and 2.2 $\\mu$m power spectra; if the amplitude or shape of the ~1-degree bump changes with season or model choice, the excess is zodiacal rather than extragalactic. Alternatively, a firm IRTS–Planck cross-correlation detection at degree scales would tie the bump to dusty galaxies, whereas a null cross-correlation with tight limits would force an unknown component.","tokens_in":16909,"feed_emoji":"🌌","tokens_out":6255,"duration_ms":58091,"temperature":0.7,"pith_summary":"This paper reports the first measurement of spatial fluctuations in the near-infrared extragalactic background light (NIREBL) at angular scales from 2 to 20 degrees, at 1.6 and 2.2 microns, using data from the IRTS satellite. The authors subtract foreground emission from zodiacal light, diffuse Galactic light, and integrated starlight, then compute the angular power spectrum of the residual sky. They find that the fluctuation amplitude declines as roughly $\\theta^{-1}$ across these large scales, and that this decline connects smoothly to previous sub-degree measurements, implying a broad bump centered near 1 degree. The paper argues that known sources—normal galaxies, high-redshift objects, intra-halo light, and the far-infrared cosmic background—cannot explain this excess. If correct, the result points to an unidentified component of the near-infrared sky that fluctuates on degree scales.","feed_headline":"Unexplained 1-degree bump in near-infrared sky glow","feed_subtitle":"IRTS data connect 2–20 degree scales to sub-degree surveys; known galaxies and dust can't explain it.","key_machinery":"The analysis rests on the angular power spectrum of a foreground-subtracted NIREBL map: fluctuations are expressed as $F(\\sqrt{l(l+1)C_l/2\\pi})$ vs angular scale $\\theta \\approx 180^\\circ/l$, computed with the pseudo-power-spectrum method on a HEALPix map of the IRTS sky. Foreground removal is the load-bearing step: the zodiacal light, which is over 90% of the raw sky brightness, is modeled with a standard zodiacal model and rescaled by per-band correction factors (0.811 and 1.146 for 1.6 and 2.2 $\\mu$m) derived from correlation fits; diffuse Galactic light is estimated from the 100 $\\mu$m dust map with an empirical scale factor; and integrated starlight is built from 2MASS and TRILEGAL Galactic model stars. The power spectrum of the residual sky is then compared with sub-degree measurements from CIBER, AKARI, and Spitzer to reveal the degree-scale bump.","core_discovery":"The central discovery is that the near-infrared extragalactic background fluctuation spectrum, measured for the first time at 2–20 degrees, joins with existing sub-degree measurements to form a broad bump centered around 1 degree. At scales larger than 2 degrees the fluctuation amplitude follows $F(\\sqrt{l(l+1)C_l/2\\pi}) \\sim \\theta^{-1}$, i.e., a nearly single power law, which indicates a random, structureless distribution of sources. Because the authors rule out foregrounds and known extragalactic contributors, the bump at approximately 1 degree is presented as an unexplained excess that needs new physics or a new population of sources.","pith_inferences":["If the zodiacal model is slightly wrong on degree scales—for example, if the interplanetary dust cloud has patchiness at the few-percent level—the apparent NIREBL bump could be a zodiacal artifact; this can be tested by comparing the fluctuation amplitude across seasons, since the zodiacal signal changes with Earth's orbit while a true extragalactic signal does not.","A possible physical origin not considered in detail here is a population of faint, unresolved sources at intermediate redshifts whose clustering peaks near degree scale; a testable prediction would be that the bump's amplitude should scale with wavelength according to the source spectral energy distribution, and that cross-correlating with future wide-area submillimeter maps would recover the sign","The same approach applied to the DIRBE time-ordered data, which cover the full sky at 1.25, 2.2, 3.5 and 4.9 $\\mu$m, could confirm or refute the bump with an independent instrument and much larger sky coverage."],"forward_implications":["The 1-degree bump, if real, implies an unknown source of near-infrared background fluctuations that previous sub-degree surveys could not see because of their limited sky coverage.","The $\\theta^{-1}$ power law at 2–20 degrees implies the large-scale NIREBL fluctuations are dominated by a random, unclustered spatial distribution rather than by known large-scale structure.","The smooth connection between IRTS and CIBER spectra means the 1-degree bump is a common feature across 1.6 to 3.6 $\\mu$m.","The IRTS–Planck cross-correlation, although an upper limit, suggests that some degree-scale fluctuations could be produced by dusty, star-forming galaxies at $z < 0.8$; better data could confirm this association.","Future deep observations such as MIRIS near the North Ecliptic Pole can directly probe the 1-degree to several-degree range and test the bump."],"supporting_citations":[{"why":"Provides the zodiacal light model that is the dominant foreground and the starting point for the per-band rescaling.","marker":"Kelsall et al. 1998"},{"why":"Supplies the CIBER sub-degree fluctuation measurements that the IRTS data connect to, forming the 1-degree bump.","marker":"Zemcov et al. 2014"},{"why":"Gives the IRTS absolute brightness measurements used to validate the NIREBL level and the correlation-based ZL correction method.","marker":"Matsumoto et al. 2015"},{"why":"Provides Spitzer fluctuation measurements at larger sub-degree scales that are compared with IRTS and scaled to 2.2 microns.","marker":"Kashlinsky et al. 2012"},{"why":"Supplies wide-field AKARI fluctuation data and the DGL measurement toward the NEP used for comparison and DGL assessment.","marker":"Seo et al. 2015"},{"why":"Provides the 2MASS star catalog used to estimate integrated starlight from Galactic stars.","marker":"Cohen et al. 2003"},{"why":"Supplies the TRILEGAL Galactic model used to estimate the contribution of stars fainter than the 2MASS limit.","marker":"Girardi et al. 2005"},{"why":"Describes the pseudo-power-spectrum estimator used to recover the NIREBL power spectrum from partial-sky observations.","marker":"Chon et al. 2004"}],"fun_headline_variants":["IRTS reveals unexplained bump in near-infrared background","Near-infrared background: 1-degree bump challenges known sources","Large-scale infrared fluctuations show excess at 1 degree","Cosmic near-IR background: known sources can't explain bump","New data connect scales, uncover mystery in near-IR glow"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The zodiacal light model, after multiplying by a single per-band correction factor fitted to the IRTS data, leaves no significant angular fluctuations at 2–20 degrees; if real interplanetary dust emission has degree-scale structure, those fluctuations directly contaminate the measured NIREBL power spectrum.","fun_headline_variants_meta":{"raw":{"variants":["IRTS reveals unexplained bump in near-infrared background","Near-infrared background: 1-degree bump challenges known sources","Large-scale infrared fluctuations show excess at 1 degree","Cosmic near-IR background: known sources can't explain bump","New data connect scales, uncover mystery in near-IR glow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00043,"raw_usage":{"total_tokens":2197,"prompt_tokens":944,"completion_tokens":1253,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":560,"completion_tokens_details":{"reasoning_tokens":1174}},"tokens_in":560,"tokens_out":1253,"duration_ms":13565,"temperature":1.0,"reasoning_tokens":1174,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:10:14.174569+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-observe the same IRTS sky fields at a different time of year (or use a different zodiacal-light model) and recompute the 1.6 and 2.2 $\\mu$m power spectra; if the amplitude or shape of the ~1-degree bump changes with season or model choice, the excess is zodiacal rather than extragalactic. Alternatively, a firm IRTS–Planck cross-correlation detection at degree scales would tie the bump to dusty galaxies, whereas a null cross-correlation with tight limits would force an unknown component.","supporting_citations":[{"cited_title":"L., & Franz, B","cited_arxiv_id":null,"evidence_quote":"Provides the zodiacal light model that is the dominant foreground and the starting point for the per-band rescaling."},{"cited_title":"2014, Sci, 346, 732","cited_arxiv_id":null,"evidence_quote":"Supplies the CIBER sub-degree fluctuation measurements that the IRTS data connect to, forming the 1-degree bump."},{"cited_title":"G., Pyo, J","cited_arxiv_id":null,"evidence_quote":"Gives the IRTS absolute brightness measurements used to validate the NIREBL level and the correlation-based ZL correction method."},{"cited_title":"G., & Ashby, M","cited_arxiv_id":null,"evidence_quote":"Provides Spitzer fluctuation measurements at larger sub-degree scales that are compared with IRTS and scaled to 2.2 microns."},{"cited_title":"J., Lee, H","cited_arxiv_id":null,"evidence_quote":"Supplies wide-field AKARI fluctuation data and the DGL measurement toward the NEP used for comparison and DGL assessment."},{"cited_title":"A., & Megeath, S","cited_arxiv_id":null,"evidence_quote":"Provides the 2MASS star catalog used to estimate integrated starlight from Galactic stars."}],"review_version":1}