{"id":"e8931fcd-510c-4920-a8c4-95141b549395","arxiv_id":"2601.16251","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"First dark-matter exclusion limits from a PbWO4 cryogenic detector (13 g, 32.4 g·d) are reported, validating the RES-NOVA concept.","lead":"A 13-gram PbWO4 crystal made from archaeological lead was operated underground as a cryogenic detector, yielding the first dark-matter exclusion limits ever obtained with this material. The measurement is a proof-of-principle for the RES-NOVA detector concept, though the limits are far weaker than current best experiments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"DM exclusion limits depend on an efficiency curve measured with a ~500 keV gamma-template pulse injection; its validity for low-energy nuclear recoils and for all event positions is unvalidated.","rationale":"The reader's weakest-assumption analysis identifies the efficiency curve in Fig. 6 as the critical unvalidated element, and I agree. This is the single most load-bearing concern because the exclusion limits are obtained by folding the expected recoil spectrum with both the measured efficiency and energy resolution. The efficiency is measured using simulated pulses generated from a high-energy (~500 keV) template; this does not establish that the pulse shape—and therefore the efficiency—is correct for low-energy nuclear recoils. The manuscript's own text acknowledges that events near threshold have signal-to-noise ratio close to unity and that the pulse-shape consistency criterion affects the spectrum down to ~2.5 keV, so the efficiency is sensitive to the very shape assumption being used. No neutron calibration or low-energy shape validation is presented. I considered other potential concerns: the SD limit relies on the collaboration's own prior form-factor calculation (Ref. [13]), but self-citation alone is not evidence of error; the lack of data/code sharing is a reproducibility issue rather than a correctness issue; and the Yellin optimum-interval method is robust to the unmodelled background in the ROI, which the paper explicitly notes. Therefore the efficiency validation is the condition on which the reported exclusion limits should be accepted. The reader's CONDITIONAL verdict is appropriate: if the authors produce a low-energy or nuclear-recoil efficiency validation and release the data, the result would be acceptable. My analysis does not change that verdict, so I recommend UNCHANGED.","tokens_in":12467,"tokens_out":4175,"duration_ms":49842,"concrete_test":"Re-measure ε(E) in Fig. 6 using a template built from low-energy events (e.g., the 46 keV 210Pb line or a 55Fe source) instead of the 500 keV average, and re-apply the same injected-pulse procedure. Then compare the two efficiency curves in the 2.5–10 keV band. If they differ by more than 20% at any ROI energy, recompute the Fig. 8 exclusion limits with both efficiencies and report the resulting shift in excluded cross-section. Alternatively, perform a neutron-source calibration (e.g., 252Cf) to directly measure nuclear-recoil efficiency and compare with Fig. 6.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section 4, the detection efficiency ε(E) of Fig. 6 is measured by injecting simulated mono-energetic pulses constructed from a signal template obtained by averaging ~1000 particle-induced pulses around 500 keV, then running them through the full analysis chain. This template is assumed to also describe low-energy nuclear recoils from DM and to be position-independent. However, cryogenic calorimeters can exhibit energy- and particle-dependent pulse shapes (e.g., thermistor non-linearity, position-dependent phonon collection), and the acceptance cut requires the OF amplitude and the MLE template fit to agree within 10%. If real low-energy nuclear recoil pulses deviate from the 500 keV gamma template, both the optimal-filter amplitude and the cut efficiency will be biased. No neutron-source calibration or low-energy electron/gamma shape measurement is reported; the only low-energy validation is the visual comparison of two events in Fig. 4. Since the 2.5–10 keV ROI contains many noise-dominated events and the limit is directly scaled by ε(E), an unvalidated efficiency is the weakest link in the central exclusion-limit claim. The manuscript itself flags shape uncertainty (Section 4: 'signals whose pulse shape is inconsistent with the detector response model'), but does not quantify its effect on the reported limits.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports the operation of a 13 g PbWO4 crystal grown from archaeological Pb, read out by a Ge NTD thermistor at ~12 mK in the LNGS Hall-C Ieti dilution cryostat. A triggerless analysis chain based on optimum filtering and MLE pulse-shape consistency is used to process 32.4 g·d of exposure. The detection efficiency is measured by injecting simulated mono-energetic pulses, built from a ~500 keV signal template, into the raw data stream. Using Yellin's optimum-interval method, the paper derives 90% C.L. exclusion limits on spin-independent and spin-dependent (on 207Pb and 17O neutrons) DM interactions. It also reports in-situ vibrational characterization with cryogenic geophones and presents the detector as a proof of principle for the RES-NOVA concept.","tokens_in":12781,"tokens_out":8624,"duration_ms":105030,"significance":"If the quantitative limits are accepted, this is the first DM exclusion limit obtained with PbWO4 as the target material, and the paper demonstrates several useful R&D advances: triggerless DAQ, a fully characterized cryogenic analysis chain, two-point gamma calibration, in-situ cryogenic vibration monitoring, and the use of archaeological Pb in an active detector element. The analysis is standard and mostly conservative: expected rates are computed with the external wimprates package, and Yellin's optimum-interval method avoids a background model. The main caveat is that the central limits inherit an unquantified systematic from the detection-efficiency measurement, which is based on a high-energy gamma/electron template and has no nuclear-recoil or position-dependent calibration. This makes the quantitative limits provisional until that uncertainty is bounded.","major_comments":[{"comment":"The detection efficiency ε(E) is measured by injecting simulated pulses built from a template averaged over ~1000 pulses at ~500 keV, then applying the full chain including the 10% OF-vs-MLE consistency cut. No nuclear-recoil calibration (e.g., neutron source) or low-energy electron/gamma shape measurement is reported, and no position dependence is measured for the 0.7×0.7×4 cm³ crystal with a single NTD thermistor. Real low-energy nuclear recoils may have a different pulse shape (thermal non-linearity, position-dependent phonon collection), which would bias both the OF amplitude and the acceptance cut. Section 4 itself states that some low-energy events have 'pulse shape inconsistent with the detector response model,' but the effect on ε(E) is not quantified. Because the limits in Fig. 8 are scaled by ε(E), this is a load-bearing systematic. Please provide a conservative efficiency enve","section":"§4, Fig. 6"},{"comment":"The absolute energy scale is anchored to the 2615 keV 208Tl line and cross-checked with the 210Pb onset at 46 keV; there is no calibration point below 46 keV. The 2.5–10 keV ROI therefore relies on an assumed linearity of the OF amplitude. The simulated line at 2.4 keV (Fig. 7) is reconstructed with the same template and thus does not validate the absolute scale for real low-energy events. A small gain nonlinearity would shift the ROI boundaries and distort the recoil-spectrum shape used in the Yellin limit. Please state the linearity assumption explicitly and quantify its impact (e.g., with an energy-scale uncertainty or a low-energy line measurement).","section":"§4–§5, energy calibration"},{"comment":"The band in Fig. 8 is described as the variation from changing the pulse-shape consistency tolerance between 5% and 20%, with a nominal 10%. This is a selection-sensitivity check, not a full systematic uncertainty estimate. The band does not cover template-mismatch, position-dependent efficiency, or energy-scale nonlinearity. Please clarify in the figure caption and text that the shown band represents only the tolerance variation, and, if possible, add the dominant efficiency/energy-scale uncertainties separately. Otherwise, the figure may be read as overstating the robustness of the limits.","section":"§5, Fig. 8"}],"minor_comments":[{"comment":"The caption uses 'amplitude spectral density (ASP)' while the text uses 'amplitude spectral density (ASD)'; please make the terminology consistent.","section":"Fig. 2 caption"},{"comment":"'as shown in see Fig. 4' should read 'as shown in Fig. 4'.","section":"§4, first paragraph"},{"comment":"The DOI '10.1103/wcvd-rk1f' appears to be a placeholder or is not in standard PRD format; please verify.","section":"Ref. [13]"},{"comment":"The abstract says 'spin-dependent interactions on neutrons' without specifying isotopes; the main text and figure caption are clearer. The statement that 183W is neglected 'due to its lower natural isotopic abundance' should also note that the omission is conservative (it can only reduce the expected rate), since spin-structure factors also matter.","section":"§5, SD limits"},{"comment":"The grey 'bulk contaminations' histogram is compared with the measured spectrum, but the normalization of this Monte Carlo component is not specified. Please state the scaling (e.g., exposure and measured activities) so the comparison is quantitative.","section":"Fig. 5 right"},{"comment":"A 30% tolerance is listed among the tested values but is not included in the Fig. 8 band (5–20%). Please explain why the 30% value is excluded from the band.","section":"§4, tolerance"}],"recommendation":"major_revision","confidential_remarks":"The paper is a legitimate R&D contribution and the 'first limit with PbWO4' claim is defensible, but the central quantitative limits depend on an efficiency curve whose applicability to low-energy nuclear recoils and to all event positions is unvalidated. If the authors can quantify this systematic (or at least provide a conservative efficiency envelope), I would support publication. The 5–20% tolerance band should not be presented as a complete systematic uncertainty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the first actual dark-matter exclusion limit with PbWO4 as the target. The earlier RES-NOVA paper (Ref [13]) was only a sensitivity projection; here they ran a 13 g crystal in a dilution fridge at LNGS for 32.4 g·d and produced real numbers. The limits are far from XENONnT/LZ/CRESST, but they are new and they validate that archaeological lead can be grown into a workable low-background crystal for cryogenic detection.\n\nWhat's done well: the analysis is careful and transparent. Trigger-less acquisition, optimal filtering, MLE fits, and Yellin's optimum interval method—no background model needed. Efficiency is estimated by injecting simulated pulses into the raw stream, and the calibration is anchored on a 2615 keV gamma line plus the 46 keV 210Pb onset. They also show in-situ geophone vibration measurements at 7 mK, which is a nice instrumentation contribution for cryogenic rare-event detectors. The paper is honest about its limitations: small exposure, high background from the cryostat, no TES readout.\n\nThe main soft spot is the efficiency curve. It's measured using a signal template built from ~1000 pulses around 500 keV (electron/gamma events), and the same efficiency is applied to low-energy nuclear recoils in the 2.5–10 keV ROI. No neutron-source calibration or low-energy line is presented to validate that template for nuclear recoils. If real recoil pulses have a different shape—position-dependent phonon collection or thermistor non-linearity—both the optimal-filter amplitude and the cut efficiency could be biased. The authors test the tolerance in the OF-vs-MLE consistency cut (5–20%) and show a band, but they don't quantify the effect of using a high-energy, electron/gamma template. This is the weakest link in the limit-setting chain.\n\nThat said, the efficiency concern is unlikely to change the big picture. The limits are many orders of magnitude away from existing bounds, and the paper's main claim is a proof of principle, not a competitive exclusion. Even a factor-of-two systematic would leave the conclusion intact.\n\nMinor notes: the SD limit relies on 207Pb form factors from the collaboration's own sensitivity paper (Ref [13]). That's a physics input, not circular, but independent verification would be nice. No data or code files are provided; the exposed code (wimprates) is public.\n\nBottom line: this is a solid instrumentation paper with a genuine first result. It deserves a serious referee. The referee should push for a discussion—and ideally a measurement—of the efficiency systematic for nuclear recoils, but I would not reject it over that.","headline":"First PbWO4 DM limits from a 13 g archaeological-Pb crystal: a careful proof-of-principle whose limits hinge on an unvalidated efficiency extrapolation.","tokens_in":13490,"tokens_out":3895,"would_cite":true,"duration_ms":48524,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.35.+d","07.20.Mc","29.40.Vj"],"model":"deepseek-v4-flash","headline":"First dark-matter limits from a lead-tungstate crystal detector","keywords":["dark matter","direct detection","cryogenic detector","PbWO4","lead tungstate","archaeological lead","RES-NOVA","bolometer"],"falsifier":"Expose the same crystal to a neutron source producing nuclear recoils of known energy in the 2–10 keV range and run the full analysis pipeline. If the measured nuclear-recoil efficiency at, say, 3 keV differs from the injected-pulse efficiency by more than the systematic band quoted in the paper, the exclusion limits would need to be revised.","tokens_in":12394,"feed_emoji":"❄️","tokens_out":7596,"duration_ms":74592,"temperature":0.7,"pith_summary":"The paper establishes that lead tungstate (PbWO4) — a heavy crystal containing both lead and oxygen nuclei — can be operated as a cryogenic calorimeter for direct dark-matter searches. Using a 13-gram crystal grown from archaeological lead and cooled to 12 mK in an underground laboratory, the authors achieved a low-energy threshold near 2.5 keV and produced the first dark-matter exclusion limits ever derived with a PbWO4 target. The limits, covering both spin-independent and spin-dependent interactions, are not competitive with leading experiments, but the demonstration of stable, low-vibration operation and a complete triggerless analysis chain is the paper's intended contribution. The work is presented as a proof of principle for the RES-NOVA detector concept, which aims to detect supernova neutrinos and dark matter with larger lead-tungstate detectors.","feed_headline":"First dark-matter limits from a lead-tungstate detector","feed_subtitle":"A 13-gram archaeological-lead crystal at 12 mK yields the first PbWO4 dark-matter limits and validates RES-NOVA.","key_machinery":"The PbWO4 crystal is the central object: its high-mass lead nuclei give strong coherent scattering for heavy dark matter, while its oxygen nuclei extend sensitivity to lower masses and, through 17O, to spin-dependent interactions on neutrons. The thermal pulses are read out with a germanium thermistor, and the analysis uses an optimum filter built from the measured noise spectrum and a ~500 keV signal template, together with a maximum-likelihood amplitude fit. A pulse-shape consistency cut rejects pathological events, and the surviving efficiency is measured by processing injected mono-energetic simulated pulses through the entire chain. The exclusion limits are derived by the optimum-interv","core_discovery":"The central claim is that a 13 g PbWO4 crystal, grown from archaeological lead to suppress intrinsic radioactivity, can serve as a working cryogenic target for rare-event searches. Operated at approximately 12 mK and read out by a germanium thermistor, the detector reached a baseline energy resolution of σ ≈ 234 eV and an analysis threshold near 2.5 keV. By injecting simulated pulses into the raw data stream, the authors measured the detection efficiency as a function of energy, and with an exposure of 32.4 g·day they set 90% confidence upper limits on dark-matter scattering using an optimum-interval method that requires no background model. These are the first such limits for a PbWO4 target","pith_inferences":["The efficiency curve is built from a signal template averaged from ~500 keV particle pulses; if low-energy nuclear recoils produce a different pulse shape, the efficiency and therefore the limits could be biased, an effect not addressed in the paper.","The energy scale is anchored to a 2615 keV gamma line and a 46 keV lead feature; sub-keV linearity of the thermistor is assumed, so the threshold region may carry unquantified energy-scale uncertainty.","Because the background in the region of interest is dominated by the cryogenic infrastructure, a dedicated low-background cryostat would likely improve the sensitivity by orders of magnitude without changing the crystal.","The simulated-pulse injection method could be extended to map position-dependent response in larger crystals, which would be needed to validate the final RES-NOVA detector geometry."],"forward_implications":["PbWO4 becomes a demonstrated target material for direct dark-matter detection, with the heavy lead and light oxygen nuclei enabling sensitivity across roughly four orders of magnitude in dark-matter mass.","The low-vibration dry cryostat and triggerless analysis chain are directly applicable to coherent elastic neutrino-nucleus scattering detectors, including real-time supernova monitoring.","The archaeological-lead radiopurity measurements support using such crystals in future ultra-low-background experiments.","Scaling the prototype to a 200 kg target mass with more sensitive thermal sensors is projected to probe previously unexplored dark-matter parameter space.","The efficiency-measurement technique described here can be reused for any new detector geometry or readout technology."],"fun_headline_variants":["13g cryogenic PbWO4 yields first dark matter limits","First dark matter limits validate RES-NOVA prototype","PbWO4 target sets first dark matter exclusion limits","13g detector at 12 mK yields dark matter limits","Archaeological lead crystal sets first dark matter limits"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The dark-matter limits assume that the detection efficiency measured using simulated pulses derived from a ~500 keV electron/gamma template applies unchanged to the low-energy nuclear recoils that dark matter would produce, at every position in the crystal.","fun_headline_variants_meta":{"raw":{"variants":["13g cryogenic PbWO4 yields first dark matter limits","First dark matter limits validate RES-NOVA prototype","PbWO4 target sets first dark matter exclusion limits","13g detector at 12 mK yields dark matter limits","Archaeological lead crystal sets first dark matter limits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001375,"raw_usage":{"total_tokens":5582,"prompt_tokens":714,"completion_tokens":4868,"prompt_tokens_details":{"cached_tokens":640},"prompt_cache_hit_tokens":640,"prompt_cache_miss_tokens":74,"completion_tokens_details":{"reasoning_tokens":4790}},"tokens_in":74,"tokens_out":4868,"duration_ms":322186,"temperature":1.0,"reasoning_tokens":4790,"cache_read_input_tokens":640,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T06:15:58.479280+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Expose the same crystal to a neutron source producing nuclear recoils of known energy in the 2–10 keV range and run the full analysis pipeline. If the measured nuclear-recoil efficiency at, say, 3 keV differs from the injected-pulse efficiency by more than the systematic band quoted in the paper, the exclusion limits would need to be revised.","supporting_citations":[],"review_version":1}