{"id":"41b70c93-b045-4c60-8f47-b05f7c1bd4cf","arxiv_id":"2501.12409","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A sensitivity study finds that a proposed 170 kg archaeological-lead PbWO4 cryogenic detector could reach spin-independent dark matter cross-sections near 2 x 10^-46 cm^2 at 20 GeV/c^2 and probe spin-dependent interactions via 207Pb.","lead":"This design study asks whether a proposed cryogenic detector made from ancient Roman lead, built to catch supernova neutrinos, could also spot dark matter particles recoiling off lead and oxygen nuclei. It estimates that after one year the detector could reach dark matter cross-sections near 2 x 10^-46 cm^2 for a 20 GeV/c^2 particle, if its energy threshold and background rejection work as assumed.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The abstract's 2 GeV/c^2 reach of 1e-43 cm^2 is in tension with the stated 1 keV threshold: 2 GeV DM produces oxygen recoils mostly below 1 keV, so the claimed sensitivity needs an unstated lower threshold or efficiency.","rationale":"The reader's weakest assumption was the 1 keV threshold as an experimental performance risk. My concern is a sharper, more specific version: the stated threshold and the quoted 2 GeV/c^2 reach are internally in tension, because oxygen recoils from 2 GeV/c^2 DM are mostly below 1 keV and Pb/W recoils are kinematically negligible. This directly affects the abstract's flagship low-mass number, making it load-bearing. I do not think it changes the verdict from CONDITIONAL, because the issue is fixable: the paper needs to state and apply the actual efficiency function and correct the low-mass projection if it was computed with a lower threshold. I considered other weaknesses—the absent numerical S_n(q) for 207Pb and the assumed 100% e/gamma rejection—but they are secondary or already framed as assumptions; the 2 GeV inconsistency sits inside the strongest claim. The 20 GeV/c^2 reach and the demonstrator concept may survive, so REJECT would be too harsh; a conditional acceptance with a corrected 2 GeV projection is appropriate.","tokens_in":13484,"tokens_out":12790,"duration_ms":125724,"concrete_test":"Recompute the m_chi = 2 GeV/c^2 SI projection of Fig. 6 using the same wimprates adaptation, but with an explicit detector response consistent with Sec. VII: Gaussian energy smearing of sigma = 200 eV, trigger efficiency eps(E) = 0.5 * erfc((1 keV - E)/(sqrt(2)*0.2 keV)) for E > 0.2 keV and 0 below, and the NR-only background spectrum of Fig. 4 (right). Derive the 90% CL upper limit with PyHF as in Sec. VI. If the resulting cross-section exceeds 1e-43 cm^2 by, say, a factor of 3 or more, the abstract's 2 GeV/c^2 claim is not supported. Also report the expected number of accepted O-recoil events at 1e-43 cm^2 and, if available, the efficiency function actually used to produce Fig. 6.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. VII states the sensitivity calculation assumes 'a 50% trigger efficiency at a nuclear recoil energy of 1 keV (with a hard cut at 200 eV)'. The abstract's headline SI projection includes 1e-43 cm^2 at 2 GeV/c^2. For a 2 GeV/c^2 WIMP, the only target in PbWO4 with observable recoil energy is O: with the paper's v_esc=544 km/s, E_max(O) = 2 mu^2 v^2 / m_O = 1.30 keV, while Pb and W recoils cap at about 0.13 keV. The standard halo spectrum for O is thus concentrated below the 1 keV threshold. Folding in a 50% efficiency at 1 keV and a hard zero below 200 eV leaves only a small high-velocity tail; at sigma = 1e-43 cm^2 the accepted signal is roughly an order of magnitude short of the ~2.3 events needed for a 90% upper limit. The 2 GeV/c^2 curve in Fig. 6 therefore cannot follow from the stated parameters. This is not only a concern that the 1 keV threshold may not be achieved: if the threshold is achieved exactly as described, the quoted low-mass reach is still unsupported unless the efficiency function is much more permissive below 1 keV than the text implies. The paper does not report the efficiency function actually used in the wimprates-based calculation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the dark matter search potential of the RES-NOVA technology demonstrator, a 170 kg PbWO4 cryogenic detector made from archaeological lead and operated at LNGS. The authors develop a Geant4-based background Monte Carlo with a list of internal, shield, surface, and external sources, and use the PyHF statistical framework together with an adapted version of the wimprates package to project 90% confidence-level sensitivities for spin-independent and spin-dependent WIMP interactions after one year of exposure. The headline results are a spin-independent reach of 1e-43 cm^2 at 2 GeV/c^2 and 2e-46 cm^2 at 20 GeV/c^2 in the most optimistic background scenario, with a claimed mass coverage from sub-GeV/c^2 to TeV/c^2.","tokens_in":13791,"tokens_out":7422,"duration_ms":74154,"significance":"If the projections were fully supported, the paper would be a useful contribution to the direct-detection literature: it is the first proposal to use archaeological Pb as an active DM target, it exploits 207Pb for spin-dependent scattering, and it is built on public tools (PyHF, wimprates) with explicit halo parameters and a documented background simulation. The strengths are the transparent use of standard rate formulas, the explicit list of background sources in Table I, and the code-based approach to the sensitivity calculation. However, the central low-mass projection is kinematically questionable, and the missing detail about the efficiency function prevents independent verification, so the significance of the headline claim is currently limited.","major_comments":[{"comment":"The quoted 90% C.L. sensitivity of 1e-43 cm^2 at m_chi = 2 GeV/c^2 is not supported by the stated detector response. With the halo parameters in Sec. VI (v_esc = 544 km/s), a 2 GeV WIMP on oxygen has maximum recoil energy E_max = 2 mu^2 v_esc^2 / m_O ~ 1.3 keV, and the recoil spectrum is concentrated at considerably lower energies. The text in Sec. VII states a 50% trigger efficiency at 1 keV with a hard cut at 200 eV; the accepted signal above 1 keV then comes only from the small high-velocity tail of the halo and is insufficient to yield a 90% upper limit at the quoted cross-section. The efficiency function actually used in the wimprates-based calculation is not reported, so this apparent inconsistency cannot be resolved from the manuscript as written.","section":"Section VII and Fig. 6, left panel"},{"comment":"The sensitivity analysis is not reproducible. The paper does not specify the energy binning, the efficiency function epsilon(E_R) used as input to wimprates, or the background counts per bin after applying the anti-coincidence logic and, in the optimistic case, the assumed 100% rejection of electron and gamma events. The statement that the input spectrum was obtained by 'simply scaling the results presented in Fig. 4' does not define the scaling procedure or the treatment of the 200 eV hard cut. Because the central projections depend directly on these choices, the authors should document them or release the analysis configuration.","section":"Section VI (Eq. 8) and Section VII"},{"comment":"The spin-dependent projection rests on the nuclear spin structure of 207Pb, but the manuscript does not report the adopted spin matrix elements S_n(q) or the momentum-transfer form factor used in the calculation. Reference [30] is cited, yet the actual numerical input and any quenching procedure are not given. Without this information the SD exclusion curve cannot be reproduced or meaningfully compared with the existing limits plotted in Fig. 6.","section":"Section III (Eq. 6) and Fig. 6, right panel"}],"minor_comments":[{"comment":"The claim of sensitivity 'from sub-GeV/c^2 to TeV/c^2' is not supported by Fig. 6, whose horizontal axis appears to start at 1 GeV/c^2; the mass range actually shown should be stated, or the sub-GeV portion should be displayed if it is covered.","section":"Abstract and Section VIII"},{"comment":"The label 'realistic prediction' appears in the figure but is not defined in the text or caption; the reader cannot tell which curve or band it refers to.","section":"Fig. 6, left panel"},{"comment":"The conversion '10^-3 c/keV/ton/s equivalent to 0.086 c/keV/kg/d' is presented as a footnote; it would be clearer to state this in the main text where the background goal is introduced.","section":"Section V, footnote 46"},{"comment":"The sensitivity definition uses the nonstandard notation 'I P' for probabilities; conventional probability notation would make the equation easier to read.","section":"Section VI, Eq. (7)"},{"comment":"The text notes that activity limits at 90% C.L. were used as central values in the Monte Carlo, but it does not state whether these were treated as fixed values or as Gaussian means with the quoted uncertainties; this should be specified.","section":"Table I and Section V"}],"recommendation":"major_revision","confidential_remarks":"The threshold problem in the low-mass reach is the principal concern. I would ask the authors to provide the efficiency curve and the exact wimprates configuration used for the 2 GeV/c^2 point, or to recompute the low-mass sensitivity under a clearly documented efficiency model. If the 2 GeV/c^2 reach cannot be recovered, the abstract and summary should be revised to remove that specific claim, since it is the paper's headline result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plain take: this is a useful sensitivity projection for using archaeological-lead PbWO4 bolometers as a DM direct-detection target. The genuinely new piece is the 207Pb spin-dependent channel with finite momentum transfer, and the qR argument for why that matters is clearly made. The background Monte Carlo is careful, the use of public tools (wimprates, pyhf) is transparent, and the conservative/optimistic bracket is honest. The neutrino-fog comparison is a nice addition. Credit is earned for those.\n\nSoft spots come in three sizes. First, the abstract promises \"detailed calculations of nuclear matrix elements,\" but the text gives a shell-model argument and a qR plot—not a numerical S_n(q). For a paper whose SD novelty rests on those elements, that is a missing piece. Second, calling RES-NOVA \"the first experiment employing Pb as an active material\" is too strong for a projection; there are other Pb-target concepts floating around, even if not precisely for DM. Tone that down.\n\nThird, and most important: the 2 GeV/c^2 reach quoted in the abstract. With the stated 1 keV threshold and \"50% trigger efficiency at 1 keV (hard cut at 200 eV)\", a 2 GeV WIMP's oxygen recoil spectrum ends near 1.3 keV and most of it sits below 1 keV. The paper never reports the actual efficiency curve used in the wimprates calculation, so the reader cannot check whether 1e-43 cm^2 at 2 GeV follows from their own inputs. If the efficiency is as restrictive as the text implies, that number is not supported. This is fixable—spell out the efficiency shape or soften the low-mass claim—but as written it is a load-bearing gap in the headline result.\n\nRecommendation: send it to peer review, expecting major revision. The detector concept and background work deserve a serious referee, and the SD channel is a legitimate complement to existing targets. But the low-mass projection needs to be re-derived with a documented efficiency function, and the SD matrix elements need to be shown or the abstract claim withdrawn. I would not cite the 2 GeV number as it stands.","headline":"A worthwhile detector-concept paper with one genuinely new SD idea (207Pb at finite momentum transfer), but the abstract oversells the low-mass reach and the SD branch lacks the promised matrix-element detail.","tokens_in":14603,"tokens_out":6151,"would_cite":true,"duration_ms":59115,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.35.+d","29.40.-n"],"model":"deepseek-v4-flash","headline":"A 170 kg PbWO4 array made from archaeological lead could, in one year, probe dark-matter nucleon cross-sections down to $1\\times10^{-43}\\,\\mathrm{cm}^2$ at 2 GeV and $2\\times10^{-46}\\,\\mathrm{cm}^2$ at 20 GeV.","keywords":["dark matter direct detection","archaeological lead","PbWO4 cryogenic detectors","WIMP","spin-dependent scattering","lead-207","CEvNS","low-background bolometer"],"falsifier":"Assemble a kg-scale archaeological-lead PbWO$_4$ crystal with the transition-edge-sensor readout and measure the trigger efficiency on nuclear recoils down to 1 keV; if the 50% point lands above 1 keV, or if the measured background in the 1\\,keV to 10\\,keV nuclear-recoil band exceeds $10^{-3}\\,\\mathrm{counts}/\\mathrm{keV}/\\mathrm{ton}/\\mathrm{s}$, the quoted one-year projections do not hold.","tokens_in":13237,"feed_emoji":"🔭","tokens_out":9174,"duration_ms":91642,"temperature":0.7,"pith_summary":"This paper argues that the same cryogenic lead-tungstate (PbWO$_4$) detectors built to catch supernova neutrinos could double as a direct dark-matter search. The key is archaeological lead: it is extremely radiopure, heavy, and rich in $^{207}$Pb, which carries unpaired nuclear spin. The authors build a Monte Carlo background model for a 170 kg demonstrator and compute sensitivity for both spin-independent and spin-dependent WIMP scattering. They find that with full electron/gamma rejection and a 1 keV nuclear-recoil threshold, one year of data could probe dark-matter-nucleon cross-sections down to $1\\times10^{-43}\\,\\mathrm{cm}^2$ at $2\\,\\mathrm{GeV}/c^2$ and $2\\times10^{-46}\\,\\mathrm{cm}^2$ at $20\\,\\mathrm{GeV}/c^2$, covering about four orders of magnitude in dark-matter mass.","feed_headline":"Ancient lead detector targets dark matter at 1e-43 cm2","feed_subtitle":"A 170 kg PbWO4 array built from archaeological lead would rival xenon detectors in one year of data.","key_machinery":"The central object is the archaeological-lead PbWO$_4$ bolometer module, in which each crystal is read out by a transition-edge sensor (a superconducting thermometer measuring the tiny temperature rise from a nuclear recoil) together with a germanium light detector that records the scintillation light. This dual heat-and-light readout allows event-by-event rejection of electron and gamma backgrounds, leaving only nuclear recoils and alphas in the region of interest. The argument also rests on two cross-section scalings: coherent elastic neutrino-nucleus scattering grows roughly with the square of the neutron number $N$, while spin-independent dark-matter scattering grows with the square of the atomic number $A$, so lead's large nucleus amplifies the signal; for spin-dependent scattering, $^{207}$Pb's simple $2p_{1/2}$ neutron-hole structure relative to doubly magic $^{208}$Pb keeps the spin matrix element reliable beyond the zero-momentum-transfer approximation.","core_discovery":"The paper's central claim is that a low-background array of PbWO$_4$ bolometers made from archaeological lead, originally proposed as a supernova-neutrino observatory, is also a viable direct dark-matter detector. Using a GEANT4-based background model that includes bulk contaminations, surface contaminants, shields, and environmental neutrons, muons, and gammas, the authors project the sensitivity of a 170 kg demonstrator after one year of operation. In the optimistic scenario with 100% rejection of electron and gamma events and a 1 keV nuclear-recoil threshold, the detector would be able to exclude spin-independent dark-matter-nucleon cross-sections down to $1\\times10^{-43}\\,\\mathrm{cm}^2$ at $2\\,\\mathrm{GeV}/c^2$ and to $2\\times10^{-46}\\,\\mathrm{cm}^2$ at $20\\,\\mathrm{GeV}/c^2$. For spin-dependent interactions, the natural abundance of $^{207}$Pb provides a neutron-sensitive target with a calculable momentum-dependent spin structure, making this the first projection with lead as the active dark-matter target material.","pith_inferences":["The paper does not single it out, but the oxygen in PbWO$_4$ is the main carrier of the sub-GeV reach: for light dark matter, lead's large mass suppresses the recoil energy, so the low-mass sensitivity quoted for the crystal depends on oxygen scattering as much as on lead's radiopurity.","The optimistic scenario assumes 100% rejection of electron and gamma events; in a real detector the sensitivity would lie between the paper's two curves, and even small leakage would shift the quoted limits upward. A graded rejection-efficiency scan would make the projection easier to test.","Because a dark-matter signal would persist while a supernova burst lasts only seconds, the same one-year dataset could be searched for the roughly 5% annual modulation the paper mentions; that analysis is not included in the projected limits.","If the planned 1.8 t upgrade operates at the same background per mass, the paper's scaling argument implies it would push sensitivity into the neutrino fog, where extra exposure alone stops helping and spectral or modulation information becomes necessary."],"forward_implications":["A 170 kg-year exposure with full electron/gamma rejection would set 90% C.L. limits of $1\\times10^{-43}\\,\\mathrm{cm}^2$ at $2\\,\\mathrm{GeV}/c^2$ and $2\\times10^{-46}\\,\\mathrm{cm}^2$ at $20\\,\\mathrm{GeV}/c^2$ for spin-independent dark-matter scattering.","The $^{207}$Pb component opens a spin-dependent channel on neutrons with momentum-dependent nuclear-structure calculations, complementing targets that must rely on the zero-momentum approximation.","After one year, the optimistic sensitivity is not background-limited, so increasing the exposure by an order of magnitude (for example to 2 ton-year) would improve the $30\\,\\mathrm{GeV}/c^2$ limit by more than an order of magnitude.","The solar and atmospheric neutrino fog sets the ultimate floor for PbWO$_4$; no amount of exposure can beat that floor, though the demonstrator remains below it.","Without particle identification, the projected limits weaken substantially, so the simultaneous heat-and-light readout is what makes the quoted dark-matter reach possible."],"supporting_citations":[{"why":"Defines the RES-NOVA concept of an archaeological-lead neutrino observatory and sets the target detector parameters that this paper adopts for the dark-matter search.","marker":"[4]"},{"why":"Supplies the measured radiopurity of archaeological Roman lead that drives the background-model assumptions for the detector and shields.","marker":"[6]"},{"why":"Establishes the low internal contamination of a kg-scale PbWO4 crystal produced from archaeological lead, giving the bulk radioactivity values used in the simulation.","marker":"[24]"},{"why":"Reports the 15 g PbWO4 proof-of-principle with transition-edge-sensor readout that motivates the projection of a 1 keV threshold on kg-scale absorbers.","marker":"[33]"},{"why":"Provides the gram-scale cryogenic calorimeter scaling law used to argue that the 1 keV threshold is achievable for the larger demonstrator crystals.","marker":"[35]"},{"why":"Supplies the recoil-rate calculation code that the authors adapted to include Pb, W, and O isotopes for the dark-matter sensitivity projections.","marker":"[50]"},{"why":"Defines the halo parameters and reporting conventions adopted for the WIMP limits, including local density and velocity distributions.","marker":"[51]"},{"why":"Defines the neutrino-fog calculation that sets the ultimate sensitivity floor shown in the projections.","marker":"[62]"}],"fun_headline_variants":["Archaeological lead bolometers target dark matter","Ancient lead array seeks dark matter at 1e-43 cm2","First dark matter projection with lead target","RES-NOVA: Lead detector probes DM to 2e-46 cm2","Sub-GeV dark matter search using archaeological Pb"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The projected sensitivity depends on the as-yet-undemonstrated ability of a kilogram-scale PbWO$_4$ bolometer to detect 1 keV nuclear recoils with 50% efficiency and a hard cutoff at 200 eV; without that threshold, the sub-GeV to few-GeV reach, including the $2\\,\\mathrm{GeV}/c^2$ limit, degrades sharply.","fun_headline_variants_meta":{"raw":{"variants":["Archaeological lead bolometers target dark matter","Ancient lead array seeks dark matter at 1e-43 cm2","First dark matter projection with lead target","RES-NOVA: Lead detector probes DM to 2e-46 cm2","Sub-GeV dark matter search using archaeological Pb"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000225,"raw_usage":{"total_tokens":1501,"prompt_tokens":1022,"completion_tokens":479,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":399}},"tokens_in":638,"tokens_out":479,"duration_ms":5333,"temperature":1.0,"reasoning_tokens":399,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:29:29.693348+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Assemble a kg-scale archaeological-lead PbWO$_4$ crystal with the transition-edge-sensor readout and measure the trigger efficiency on nuclear recoils down to 1 keV; if the 50% point lands above 1 keV, or if the measured background in the 1\\,keV to 10\\,keV nuclear-recoil band exceeds $10^{-3}\\,\\mathrm{counts}/\\mathrm{keV}/\\mathrm{ton}/\\mathrm{s}$, the quoted one-year projections do not hold.","supporting_citations":[{"cited_title":"RES-NOVA: A new neutrino observatory based on archaeological lead","cited_arxiv_id":"2004.06936","evidence_quote":"Defines the RES-NOVA concept of an archaeological-lead neutrino observatory and sets the target detector parameters that this paper adopts for the dark-matter search."},{"cited_title":"Radiopurity of an archeological Roman Lead cryogenic detector","cited_arxiv_id":"1904.04040","evidence_quote":"Supplies the measured radiopurity of archaeological Roman lead that drives the background-model assumptions for the detector and shields."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the low internal contamination of a kg-scale PbWO4 crystal produced from archaeological lead, giving the bulk radioactivity values used in the simulation."},{"cited_title":"Operation of an archaeological lead PbWO$_4$ crystal to search for neutrinos from astrophysical sources with a Transition Edge Sensor","cited_arxiv_id":"2111.07638","evidence_quote":"Reports the 15 g PbWO4 proof-of-principle with transition-edge-sensor readout that motivates the projection of a 1 keV threshold on kg-scale absorbers."},{"cited_title":"Aalbers, B","cited_arxiv_id":null,"evidence_quote":"Supplies the recoil-rate calculation code that the authors adapted to include Pb, W, and O isotopes for the dark-matter sensitivity projections."}],"review_version":1}