{"id":"35d431e0-7cb0-49ff-b01a-34dcf0a6817f","arxiv_id":"2506.04570","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A broad-band imaging method finds 24 candidate Lyman-alpha nebulae around z=1.9-3.0 quasars, with tentative null results on morphology evolution and environmental dependence.","lead":"This pilot survey used existing wide-field ultraviolet and optical images to search for giant glowing gas clouds around 483 distant quasars, and found 24 candidate nebulae. The same broad-band technique could be applied to the Rubin Observatory's LSST survey, covering thousands of times more sky.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 24 candidate nebulae may be dominated by systematic residuals: the color-term relations (Eqs. 1-6) are calibrated on non-quasar galaxies, and no control sample demonstrates that the broad-band selection has a low false-positive rate.","rationale":"The reader's verdict is CONDITIONAL, and my stress-test does not move it. The paper is a pilot survey with appropriately hedged language: it calls the objects 'candidate quasar nebulae,' states that spectroscopic confirmation is missing (Sec. 6), and provides processed images as online material so readers can inspect individual cases. The Monte Carlo covering-fraction analysis and comparison with previous IFU radial profiles are genuine independent checks of the method's recovery properties. My concern targets a different, unquantified property: the false-positive rate of the selection itself. The reader identified the same root assumption (the color-term calibration on non-quasar galaxies in Section 3), and the paper gives no control experiment that would show the 5.9% detection rate exceeds what pure systematics would produce. This matters because the abstract's strongest claim is not just that 24 objects exist, but that broad-band photometry alone can select extended Lyα nebulae; without a false-positive control, the method's selectivity is unproven. I also note a separate numeric inconsistency: the conclusions bullet reports 39 and 17 candidates with effective areas >40 arcsec2 (56 total), whereas Table 2 and Section 5.1 consistently give 24; this typo should be corrected but it does not affect the main concern. Nevertheless, the central issue is a resolvable empirical question rather than an internal inconsistency, and the authors explicitly frame the result as a feasibility demonstration requiring follow-up. Hence CONDITIONAL (UNCHANGED) is the right verdict.","tokens_in":29188,"tokens_out":11459,"duration_ms":135301,"concrete_test":"Run the identical pipeline on a matched control sample: ~500 bright stars and ~500 random sky positions in the same HCD-JF fields, with the same seeing matching, photo-z masking, u/g continuum subtraction (Eqs. 1-6), 2σ segmentation, and Area_eff>40 arcsec2 threshold. If the control 'nebula' fraction is statistically comparable to 24/409 (≈5.9%), the candidate list and the LSST extrapolation are not evidence for selective Lyα-nebula detection; if it is much lower (e.g., <1%), the color-term systematics are likely subdominant and the candidates are strengthened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"For the abstract's central claim to be correct, the residual u/g flux at 50-170 kpc around the quasars must be dominated by Lyα emission rather than continuum-subtraction error or unmasked interlopers. Section 3 constructs the continuum using linear combinations of g,r (or r,i) whose coefficients are empirical fits to 1,541-1,822 spec-z galaxies, explicitly excluding SDSS/eBOSS quasars (Eqs. 1-6). Quasar rest-UV continua are frequently dominated by AGN power laws rather than the galaxy SEDs used in the fit, so the fitted galaxy color terms can be systematically offset for the actual targets. Because the subtraction is applied pixel-by-pixel over the whole cutout, any such offset, together with color-dependent PSF wings and diffraction spikes from the central quasar, produces a spatially extended residual that survives the central 2.5 arcsec mask and mimics a low-surface-brightness nebula. The photo-z masking is also incomplete: Eq. 8 gives an outlier rate of ~13% for i<25 sources at z=1.9-3.0, so unmasked projected neighbors can create the asymmetric, elongated features used to identify the 'giant' candidates. No control experiment (stars, random sky positions, or non-quasar galaxies) is presented to quantify the false-positive rate of the Area_eff>40 arcsec2 selection, so the 24/409 detection fraction (5.9%) has no demonstrated excess over systematics. This is the load-bearing weakness: if the false-positive rate is comparable, the method is not yet shown to be selective for Lyα nebulae and the LSST projection in Section 7 loses its empirical support. The paper's own caveats (Secs. 3 and 6) acknowledge missing spectroscopic confirmation and shallow depth, but they do not substitute for a false-positive control.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a pilot broad-band search for extended Lyman-alpha nebulae around 483 SDSS/BOSS quasars at z=1.9-3.0, using deep u-to-K imaging from the CLAUDS-HSC-SSP-DUNES2 joint dataset. The method constructs a continuum image from two broad bands (Eqs. 1-6), subtracts it from the Ly-alpha band, masks foreground/background objects with photometric redshifts, and selects 24 candidate quasar nebulae with effective Ly-alpha area Area_eff>40 arcsec2. The authors then measure radial profiles, covering fractions, asymmetry, and environmental overdensities, concluding that there is no significant redshift evolution of asymmetry and no environmental dependence, while noting the lack of spectroscopic confirmation and various data limitations.","tokens_in":29485,"tokens_out":3974,"duration_ms":42646,"significance":"If validated, the broad-band Ly-alpha imaging technique would be a powerful way to discover giant Ly-alpha nebulae over very wide fields, including the upcoming Rubin/LSST survey. The paper is careful in tone: candidates are explicitly called candidates, the absence of spectroscopic confirmation is acknowledged, and several limitations are stated in the text. The public release of the candidate catalogue and processed images is a strength. However, the central claim that 24 candidate nebulae exhibit genuine u/g-band excess rests on an empirical continuum extrapolation and has not been tested against control samples or systematics. The paper is therefore a useful pilot and method demonstration, but its main quantitative result is not yet firmly established.","major_comments":[{"comment":"The selection of 24 quasar nebulae via Area_eff>40 arcsec2 is not accompanied by any control experiment using stars, random sky positions, or non-quasar galaxies processed through the same pipeline. The 5.9% detection fraction (24/409) therefore has no demonstrated excess over the expected rate of spurious extended residuals from PSF wings, colour-term mismatch, or unmasked projected neighbours. The central claim of the abstract depends on this comparison, so a control experiment should be added.","section":"Section 4.1, Area_eff selection"},{"comment":"The colour-term relations are fitted to 1,541-1,822 spectroscopic galaxies, with SDSS/eBOSS quasars explicitly excluded. Quasar rest-UV continua are often AGN power laws rather than galaxy SEDs, so these galaxy-calibrated relations may be systematically offset when applied to the quasar targets. Because the continuum subtraction is performed pixel-by-pixel over the full cutout, even a small systematic offset will produce a spatially extended residual that survives the 2.5-arcsec central mask and mimics a low-surface-brightness nebula. The paper should validate the colour terms on quasar spectra, on quasars whose Ly-alpha does not fall in the Ly-alpha band, or on stars, and quantify the resulting spatial residual.","section":"Section 3, Eqs. (1)-(6)"},{"comment":"The photometric-redshift outlier rate is f_outlier=0.13 for i<25 sources at z=1.9-3.0, and the mask only removes objects with |delta z|>0.1. A substantial fraction of foreground/background galaxies therefore remains unmasked and can create asymmetric extended features in the Ly-alpha residual images. This concern is directly relevant to the morphology analysis: Section 5.1 notes that three of the nebulae with alpha<0.5 have bright companions in the Ly-alpha images, suggesting that their high asymmetry is not intrinsic. The paper should quantify the probability that an Area_eff>40 arcsec2 detection is produced by an unmasked interloper, for example by injecting fake galaxies into the images.","section":"Section 3, Eq. (8) and Section 5.1"},{"comment":"The environmental analysis uses photometric redshifts with outlier rates of 0.3-0.4 at |delta z|>0.05, and the control and quasar samples are both subject to the same projection effects. The KS tests therefore have limited power, and the conclusion that there is 'no environmental dependence' is not strongly supported. This is acknowledged in the text, but given that the absence of environmental dependence is one of the three headline results, the limitation should be stated more prominently and the density measurement should include an uncertainty estimate arising from photo-z scatter.","section":"Section 5.2"}],"minor_comments":[{"comment":"The Conclusions state that the authors obtain 39 and 17 quasars with large Ly-alpha nebulae in the two redshift intervals, but Section 4.1 and Table 2 report 8, 10, and 6 in the three selections, totaling 24. The abstract also says 24. This numerical inconsistency should be corrected.","section":"Section 6"},{"comment":"The notation for the background standard deviation is inconsistent: the text uses both bg_std and bg_std_z, and Figure 6 uses bd_std_z. Please unify the notation.","section":"Section 4.1"},{"comment":"The definition of 'Type I ELAN' from Li et al. (2024) is cited as '>100 pkpc with M_UV<-22', but M_UV is not defined at that point. Please define the absolute UV magnitude and specify the band or wavelength.","section":"Section 5.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is appropriate for a journal like MNRAS as a pilot methodology paper. The main issue is the absence of control experiments and validation of the colour-term extrapolation for quasar continua; these are fixable within the scope of the paper. I would be willing to reconsider after the authors add a control analysis and correct the numerical inconsistency in the Conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper before trusting the 24 candidate nebulae. First, the analysis is careful, transparent, and appropriately hedged — the authors call them candidates and openly list the limitations. Second, the selection is not yet demonstrated to be selective: there is no control experiment showing that the broad-band subtraction yields a low false-positive rate. That missing control is the load-bearing soft spot.\n\nWhat is actually new: the extension of Shimakawa 2022 to lower redshift using u/u*-band data, the addition of photo-z masking, and the first comparison of nebula morphology and environment in this context. The data are deep, the sample selection is clearly described, and the processed images are released as online material, which is good practice. The radial-profile comparison with IFU surveys is a reasonable sanity check, and the covering-fraction Monte Carlo is a useful estimate of detection loss.\n\nThe soft spots are real but proportionate. The color-term relations (Eqs. 1–6) are fitted to spec-z galaxies with quasars excluded, so the residual u/g flux around quasars could be affected by AGN power-law continua, PSF wing differences, or diffraction spikes. The photo-z outlier rate of ~13% means unmasked foreground/background neighbors can create asymmetric, elongated features that mimic giant nebulae. The paper does not present a control sample of stars, random positions, or non-quasar galaxies to quantify the contamination rate of the Area_eff > 40 arcsec2 selection. Without that, the 5.9% detection fraction (24/409) has no demonstrated excess over systematics.\n\nThat said, this is a pilot survey, and the authors say so. The abstract uses the word “candidate,” and the conclusions call for follow-up spectroscopy. The LSST projection in Section 7 is speculative but not unreasonable; it just inherits the same false-positive uncertainty. I do not think the central argument collapses — the method is plausible and was previously tested in a different context — but the evidence is not yet strong enough to claim the 24 objects are reliably extended Ly-alpha nebulae.\n\nThis paper deserves peer review. A serious referee should ask for a false-positive control (e.g., running the same selection on stars or random sightlines) or at least a quantitative estimate of the expected contamination rate. The candidate list is worth pursuing spectroscopically, and the methodology discussion is useful for LSST preparations. I would bring it to a reading group focused on high-z gas or survey methods, and I would cite it as a pilot if I were working on broad-band Ly-alpha searches.","headline":"A careful pilot of a known broad-band Ly-alpha selection method, useful for LSST planning, but the candidate list is not yet proven and a false-positive control is missing.","tokens_in":30182,"tokens_out":2357,"would_cite":true,"duration_ms":30968,"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":"Extended Lyman-alpha nebulae around $z=1.9$--$3.0$ quasars can be found in ordinary broad-band images; this pilot finds 24 candidates across 13 square degrees and paves the path to a Rubin-LSST-wide census.","keywords":["Lyman-alpha nebulae","quasar nebulae","broad-band selection","ELAN","circumgalactic medium","Rubin LSST","photometric redshifts","HSC-SSP"],"falsifier":"Run the identical colour-term subtraction and nebula selection on thousands of random sky positions or on non-quasar galaxies matched to the same photometric-redshift distribution and imaging depth; if the rate of 50-170 kpc excess regions approaches the roughly 5% rate found around quasars, the candidates are dominated by colour-term systematics. Spectroscopy settles it directly: integral-field or narrow-band observations of the 24 candidates should show extended Ly-alpha with rest-frame equivalent width above roughly 100-200 A and spatial morphology matching the broad-band excess.","tokens_in":28940,"feed_emoji":"🌌","tokens_out":9278,"duration_ms":79572,"temperature":0.7,"pith_summary":"The paper sets out to show that huge Lyman-alpha nebulae around quasars can be found without expensive narrow-band or integral-field observations, using only the broad-band images that large surveys already produce. Around 483 spectroscopically confirmed quasars at $z=1.9$--$3.0$ in 13 square degrees of deep $u$-to-$K$ imaging, the authors extrapolate a continuum image from two redder bands into the $u$ or $g$ band, subtract it, and look for residual extended emission. They find 24 candidate nebulae with $u$- or $g$-band excess spanning 50--170 kpc, three of which rival known giant nebulae in size and asymmetry. The Rubin Observatory's LSST will cover 18,000 square degrees in the same bands, so the same method could scan hundreds of thousands of quasars and, for the first time, measure how common these giant gas reservoirs are and where they live.","feed_headline":"24 extended quasar nebulae found in broad-band images alone","feed_subtitle":"A 13-square-degree pilot traces Ly-alpha gas around z~2-3 quasars via u/g-band excess, priming the Rubin LSST era","key_machinery":"The broad-band colour selection is the mechanism: three images per quasar — one Ly$\\alpha$-band ($u$, $u^*$, or $g$) and two redder continuum bands ($g,r$ or $r,i$) — are combined through empirical colour-term relations (Eqs. 1--6) to build a pseudo-continuum image of the Ly$\\alpha$ band, which is then subtracted to leave a Ly$\\alpha$ map. The colour terms and their small redshift dependences are fitted to thousands of spec-$z$ galaxies from the same imaging, and an SED-fitting photometric-redshift catalogue is used to mask projected neighbours outside $|\\delta z| > 0.1$, rejecting about 95% of $i<25$ sources as irrelevant foreground or background objects. Residual flux is converted to Ly$\\alpha$ surface brightness using the filter width, the filter transmission at the quasar's Ly$\\alpha$ redshift, and a pixel area, then scaled by $(1+z)^4$ to a common $z=2.3$ for comparison. Detection requires a 2$\\sigma$ excess in binned pixels, and a nebula is defined by an effective area above 40 arcsec$^2$ after masking.","core_discovery":"The central claim is that extended Lyman-$\\alpha$ nebulae around $z=1.9$--$3.0$ quasars are detectable as a residual excess in a single broad-band image after subtracting an empirically extrapolated continuum. For quasars at $z=1.90$--$2.23$ the Ly$\\alpha$ line falls in the $u$ (or $u^*$) band, and for $z=2.34$--$3.00$ in the $g$ band; the continuum at that wavelength is estimated from two redder bands using colour-term relations calibrated on 1,541--1,822 galaxies with spectroscopic redshifts, with a small redshift correction. Photometric-redshift masking removes foreground and background objects that would mimic nebulae, and the central quasar is masked. Applying this to 483 quasars yields 24 candidates with effective Ly$\\alpha$ areas above 40 arcsec$^2$, corresponding to sizes of roughly 50--170 kpc, whose radial surface-brightness profiles agree with the exponential profiles measured by integral-field surveys, with scale lengths of about 15--17 kpc. The authors also report no significant difference in nebular asymmetry between the $z\\sim2.1$ and $z\\sim2.8$ samples, and no statistically significant environmental overdensity around quasars with large nebulae compared to control samples, although the largest nebula sits near a density peak.","pith_inferences":["The null environmental result may be a depth effect: the Monte Carlo simulations show the method recovers only about 10% of the Ly$\\alpha$ flux beyond 40 pkpc, so any environmental trends acting on the outer nebular gas would largely escape this selection.","A direct test of the systematics would be to apply the same colour-term subtraction to spectroscopically confirmed non-quasar galaxies at the same redshifts; the false-positive rate of 'nebulae' around those galaxies would calibrate how much of the 5% detection rate is continuum mismatch rather than real Ly$\\alpha$ emission.","Because the quasar sample is fainter than the integral-field-survey quasars ($-27.7 < M_i < -24.9$ versus $M_i < -27$), the absence of morphological evolution between $z\\sim2.1$ and $z\\sim2.8$ cannot yet be compared with earlier claims of evolving asymmetry; LSST's larger sample would allow a luminosity-matched comparison.","The same broad-band trick could in principle be pushed to the $g$-band at $z=3$--$3.5$ or the $i$-band at $z\\sim4$, but the paper's own colour-term analysis shows IGM absorption increasingly corrupts the continuum extrapolation beyond $z=3$, setting a practical boundary for the method."],"forward_implications":["If the method is valid, the Rubin LSST's 18,000 square degrees in the same $u,g,r,i$ bands would allow the selection to run over roughly 200,000 quasars, expanding the surveyed volume by a factor of hundreds relative to current narrow-band surveys.","A statistically meaningful census of giant (ELAN-like) nebulae becomes possible, including their number density and their association with protoclusters, which presently rests on a handful of objects.","The $u$-band path extends Ly$\\alpha$ searches to $z=1.9$--$2.3$, where the filter is about two times narrower than the $g$-band, giving better line contrast and reaching surface brightnesses as low as roughly $0.5\\times10^{-17}\\,\\mathrm{erg\\,s^{-1}\\,cm^{-2}\\,arcsec^{-2}}$.","The three giant candidates, with projected extents beyond 130 kpc and asymmetric shapes, would join the rare Type I ELAN class if confirmed, providing targets for detailed studies of gas around massive haloes."],"supporting_citations":[{"why":"Supplies the broad-band Ly-alpha selection method this paper extends to u-band selections at lower redshift.","marker":"Shimakawa 2022"},{"why":"Provides the DR16Q quasar catalogue that defines the 483 targets and their redshifts and absolute magnitudes.","marker":"Lyke et al. 2020"},{"why":"Delivers the deep CLAUDS u/u*-band imaging that makes the z=1.9-2.3 search reach surface brightnesses near 0.5e-17 erg/s/cm2/arcsec2.","marker":"Sawicki et al. 2019"},{"why":"Supplies the HSC-SSP g,r,i,z,y imaging and catalogue system the joint dataset is built on.","marker":"Aihara et al. 2019"},{"why":"Provides the SED-fitting photometric redshifts used to mask neighbouring contaminants and to measure environmental overdensities.","marker":"Tanaka et al. 2018"},{"why":"Offers the IFU radial-profile benchmark and the roughly one-percent ELAN abundance estimate the results are compared against.","marker":"Arrigoni Battaia et al. 2019"},{"why":"Supplies the z~2.25 IFU comparison sample and the claim of evolving Ly-alpha morphology that this study tests.","marker":"Cai et al. 2019"},{"why":"Defines the Type I ELAN criterion (>100 pkpc, M_UV < -22) used to classify the three giant candidates.","marker":"Li et al. 2024"}],"fun_headline_variants":["24 quasar nebulae found via broad-band excess in pilot","Pilot survey finds 24 Ly-alpha nebulae around quasars","Broad-band imaging reveals 24 quasar nebulae for Rubin","24 candidate quasar nebulae from single-band images","No density bias seen in 24 quasar nebulae from pilot"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The colour-term relations that predict Ly-alpha-band flux from redder bands were fitted to ordinary galaxies, not quasars; the method assumes those relations hold for quasar host galaxies and their surroundings, so that any leftover u- or g-band light is genuinely extended Ly-alpha emission rather than a mismatch in the continuum extrapolation or light from other ultraviolet lines.","fun_headline_variants_meta":{"raw":{"variants":["24 quasar nebulae found via broad-band excess in pilot","Pilot survey finds 24 Ly-alpha nebulae around quasars","Broad-band imaging reveals 24 quasar nebulae for Rubin","24 candidate quasar nebulae from single-band images","No density bias seen in 24 quasar nebulae from pilot"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000939,"raw_usage":{"total_tokens":4158,"prompt_tokens":1231,"completion_tokens":2927,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":847,"completion_tokens_details":{"reasoning_tokens":2841}},"tokens_in":847,"tokens_out":2927,"duration_ms":21213,"temperature":1.0,"reasoning_tokens":2841,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:39:08.973395+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the identical colour-term subtraction and nebula selection on thousands of random sky positions or on non-quasar galaxies matched to the same photometric-redshift distribution and imaging depth; if the rate of 50-170 kpc excess regions approaches the roughly 5% rate found around quasars, the candidates are dominated by colour-term systematics. Spectroscopy settles it directly: integral-field or narrow-band observations of the 24 candidates should show extended Ly-alpha with rest-frame equivalent width above roughly 100-200 A and spatial morphology matching the broad-band excess.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the broad-band Ly-alpha selection method this paper extends to u-band selections at lower redshift."}],"review_version":1}