{"id":"ac00a01a-2afa-49fa-9da6-f9734d58fa1e","arxiv_id":"1909.02568","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"XENON1T S2 data constrain the solar scalar-electron coupling below 2e-15 at 90% CL, about three times weaker than the best stellar cooling bound.","lead":"A new analysis of XENON1T and LUX data looks for a light scalar particle emitted from the Sun's core and absorbed by xenon atoms. It sets new direct bounds on the particle's coupling to electrons, strongest in the lower-threshold S2 ionization channel, though still a few times weaker than indirect stellar cooling limits.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The S2-only headline bound rests on the unvalidated resonant production formula Eq. (12), not on the verified non-resonant relations; this is the least-supported element of the central claim.","rationale":"The paper's S1 bound is built on Eq. (4) applied to documented photon absorption data and public spectra; there the reader's Eq. (4) concern is mitigated because the relevant electrons are non-relativistic and the authors explicitly confirmed Eqs. (5) and (9) by direct computation. The S2 bound, however, is the headline number and is sensitive almost entirely to the near-resonance flux at omega ~ 0.2 keV described by Eq. (12). That equation is imported from Refs. [26,32,33] without an independent derivation or a stated value of Gamma_L in this work, and the paper's confirmation note covers only the non-resonant relations. This is a concrete validation gap rather than an internal inconsistency. A moderate error in the resonant flux would shift g by a factor of order 10^{1/4}, and a large error could move the S2 bound closer to or beyond the stellar cooling reference. The reader's CONDITIONAL verdict already captures the need for additional documentation of approximations, so my read does not change the verdict; I only sharpen the specific formula that should be checked.","tokens_in":9416,"tokens_out":28281,"duration_ms":314368,"concrete_test":"Compute the solar scalar flux independently for omega in [0.1, 0.3] keV using the RPA longitudinal dielectric function for a non-relativistic electron-ion plasma with the solar model of Ref. [31], and compare the height and integrated area of the resonance peak with the prediction of Eq. (12) for g = 1e-14. In parallel, rerun the XENON1T S2-only likelihood with the resonant term in Eq. (12) set to zero; if the resulting 90% CL bound is weaker than about 6e-15, the headline result is controlled by the unvalidated resonance term.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline result, g_phi_ee < 2e-15 (Eq. 17), is driven by the resonant enhancement of the solar scalar flux at omega ~ omega_p ~ 0.2 keV, visible as the peak in Fig. 2. The production rate used for this resonance is Eq. (12), a Lorentzian in (1 - omega_p^2/omega^2)^2 whose prefactor is the same scalar-photon relation as Eq. (4) but whose width Gamma_L is the damping rate of the longitudinal plasmon. The paper states that Eqs. (5) and (9) were confirmed by direct computation in the non-relativistic limit, but no such confirmation or derivation is reported for Eq. (12), and the numerical implementation of Gamma_L from the solar opacity data is not specified. If the resonant flux is over- or under-estimated by a factor of 10, the S2 event rate, which scales as g^4, shifts the 90% CL bound by a factor of about 1.8; a larger error in the Lorentzian prefactor or width would move the headline number more substantially. The central claim therefore rests on the least-documented element of the calculation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper computes the flux of sub-keV scalar particles produced in the solar core through electron-ion Bremsstrahlung, Compton-like scattering, atomic transitions, and resonant longitudinal-plasmon excitation, and estimates their absorption rate in liquid xenon through a scaled photoelectric cross-section. Using LUX and XENON1T data, the authors quote 90% CL upper limits g_φee < 7e-15 (XENON1T S1-only) and g_φee < 2e-15 (XENON1T S2-only), with the S2 limit valid for m_φ ≲ 0.1 keV and driven by the resonant peak near ω_p ~ 0.2 keV. These are translated into Higgs-mixing bounds for a relaxion-like scalar, reaching sinθ < 7e-10 for the standard electron-Higgs coupling.","tokens_in":9559,"tokens_out":14211,"duration_ms":155904,"significance":"If the central result holds, the S2-only bound is the strongest direct detector limit on a solar scalar coupled to electrons and probes natural Higgs-portal parameter space. The paper has genuine strengths: the non-resonant production formulas are cross-checked by a direct non-relativistic computation, the flux calculation uses tabulated solar opacity via detailed balance without free normalizations, and the S1 analysis includes several explicit systematic checks with reported <5% shifts. The claims are carefully delimited by mass ranges, and the scaling of the sensitivity with exposure and with g^4 is made explicit. However, the headline S2 limit rests on the least-documented part of the calculation, namely the resonant production formula Eq. (12), and on a partial background model, so the quantitative claim needs substantiation before it can be accepted.","major_comments":[{"comment":"The S2-only headline limit Eq. (17) is produced by the resonant peak at ω ~ 0.2 keV shown in Fig. 2. The resonant production rate is given by Eq. (12), but this formula is presented without derivation and, more importantly, without specifying how the longitudinal damping rate Γ_L is computed from the solar opacity data or from any other input. The direct non-relativistic cross-check claimed in the text covers only Eqs. (5) and (9), not Eq. (12). Since the event rate scales as g^4, a factor-of-10 error in the resonant flux changes the quoted bound by about 1.8, so the current text leaves the headline number unreproducible. The authors should derive Eq. (12), give the explicit expression and numerical implementation of Γ_L (including the relevant plasma absorption channels and their relation to k(ω)), and validate the resonant flux against Ref. [32] or a dedicated independent calculation.","section":"Section II, Eq. (12)"},{"comment":"The absorption cross-section used for both S1 and S2 limits is obtained by applying the free-electron matrix-element relation Eq. (4) to the measured photoelectric cross-section of xenon. The authors state that Eq. (4) was verified in the non-relativistic limit for the bremsstrahlung and Compton production amplitudes, but Eq. (13) involves bound atomic electrons and final-state interactions in the photoelectric process, which are not validated by that check. Since all quoted limits scale directly with this cross-section, the relation needs either a dedicated derivation for atomic absorption (e.g., with xenon wavefunctions and a treatment of the bound-state matrix elements) or an explicit estimate of the error incurred by using the free-electron ratio. Without this, the absolute normalization of every quoted limit is uncertain.","section":"Section III, Eq. (13)"},{"comment":"The 90% limit Eq. (17) is obtained with a 'partial background model' from XENON1T, but the manuscript does not specify which background components are included, how the model is normalized, or which systematic uncertainties are propagated. In the resonant region the signal is a modest contribution on top of a sizable background, so a few-percent background normalization or shape error can shift the limit by more than the precision implied by Eq. (17). The authors should either use the full background model of Ref. [42] or repeat the fit with nuisance parameters and report the resulting systematic error on g_φee.","section":"Section III, S2-only analysis and Fig. 2"},{"comment":"The mass-dependent limits shown in Fig. 3 are not documented. Section II states that the flux is computed in the massless limit, while the results claim bounds for m_φ up to 1 keV and present upper limits for masses above 0.01 keV. The factors v_φ in Eqs. (4)-(13) and the kinematic lower bound ω > m_φ depend on the mass. The authors should state explicitly how the scalar mass is implemented in the flux and cross-section used to produce Fig. 3, including the treatment of the phase-space threshold.","section":"Section II and Fig. 3"}],"minor_comments":[{"comment":"The function F(a) is defined in Eq. (7), but the parameter a is not defined in the displayed equation; from Eq. (6) it is a = ω/T. Please define this explicitly.","section":"Section II, Eq. (7)"},{"comment":"The text describes Γ_bb as production from transitions of 'bounded electrons'; this should read 'bound electrons'.","section":"Section II"},{"comment":"The number of bins and the exact S2 energy window used for the S2-only profile likelihood are not stated. Please specify the S2 range and how the 186 eV threshold is mapped into S2, so that the analysis is reproducible.","section":"Section III, S2-only analysis"},{"comment":"The abstract quotes sinθ < 2e-9 (7e-10) without labeling whether these correspond to the S1 or S2 limits; the body text Eq. (19) gives 7e-10 and 1e-12 for κ_e = 1 and 600. Please make the correspondence between the abstract values and Eq. (19) explicit.","section":"Abstract and Eq. (19)"},{"comment":"The axis labels and caption of Fig. 3 are garbled in the version provided; the final submission should have readable labels for the coupling and scalar mass axes.","section":"Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely and interesting question, and the authors have made a serious attempt to derive both the flux and the detection rate from first principles. The main problem is that the headline S2 limit depends on Eq. (12) and on the background treatment, neither of which is documented sufficiently for the result to be reproduced or certified. These are likely fixable within the scope of the manuscript, so I recommend major revision rather than rejection. I would also ask the editor to ensure the requested derivation and numerical specifications are actually added rather than deferred to a companion paper or code release."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nRead this one if you care about light scalar searches in direct detection. It takes the established solar-axion machinery, applies it to a scalar coupled to electrons, and produces the first direct detector limits on g_phi ee from LXe data: 7e-15 from XENON1T S1-only, 2e-15 from the S2-only analysis (valid below ~0.1 keV). Those numbers are a factor few weaker than the current stellar cooling bound with in-medium mixing, which the authors state openly. That honesty is typical of the whole paper.\n\nWhat is genuinely new: the solar scalar flux computation including bremsstrahlung, Compton-like scattering, atomic transitions, and resonant plasma production, and the use of the XENON1T low-energy S2 data to capture the resonant peak near 0.2 keV. The non-resonant rates are cross-checked against a direct non-relativistic computation, and the absorption cross-section is tied to measured photoelectric cross-sections. The analysis uses public data and published detector response models, and the statistical procedure is standard profile likelihood.\n\nSoft spots, in proportion: the S1-only XENON1T limit rests on a single bin with an assumed flat background (66 expected, 70 observed) — that is rough, and the paper admits it. More importantly, the resonant production formula, Eq. (12), is the least-documented piece of the calculation. It is not unvalidated—it is imported from Hardy and Lasenby—but the implementation of the longitudinal plasmon damping rate from solar opacity data is not spelled out, and a factor of several in that rate would shift the headline S2 bound by a fourth-root factor. The solar model and opacity systematic uncertainties are not propagated into the limits, which I would have liked to see quantified. None of these are fatal; the central result as an order-of-magnitude constraint is sound.\n\nThis paper deserves peer review as is. It is a clean, useful addition to the direct-detection literature, and the S2 analysis is worth having on record even if the bound is not the world's best. I would bring it to a reading group and would cite it in work on Higgs-portal scalars. Recommendation: accept with minor revision requests asking for the Gamma_L implementation details and a sensitivity estimate for solar-model uncertainties.","headline":"Solid, honest limits paper: first LXe bounds on solar scalars coupled to electrons, with a defensible S2 result from resonant production; a few unpropagated systematics but the central numbers hold.","tokens_in":10203,"tokens_out":2392,"would_cite":true,"duration_ms":24610,"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 Sun's core would shine scalars that liquid xenon can absorb; XENON1T and LUX data bound the scalar–electron coupling at $g_{\\phi ee} < 2\\times10^{-15}$.","keywords":["light scalar","relaxion","solar production","electron coupling","liquid xenon detectors","XENON1T","LUX","resonant plasma production"],"falsifier":"Compute the exact electron-ion Bremsstrahlung rate for scalar emission in a screened solar plasma at energies $\\omega\\approx0.1$–$0.3$ keV and compare it with the photon rate scaled by $g_{\\phi ee}^2 v_\\phi^2/(4\\pi\\alpha_{\\rm em})$; a deviation by an order of magnitude would move the S2 coupling bound by about $10^{1/4}\\approx1.8$ and could erase the claimed exclusion.","tokens_in":9151,"feed_emoji":"☀️","tokens_out":7039,"duration_ms":63935,"temperature":0.7,"pith_summary":"The paper argues that a sub-keV scalar particle that couples to electrons would be produced in the solar core—mainly by electron-ion Bremsstrahlung—and would be absorbed by liquid xenon atoms through an analogue of the photoelectric effect, leaving a detectable scintillation (S1) or ionization (S2) signal. Using public XENON1T and LUX data, the authors derive 90% confidence upper limits on the scalar–electron coupling $g_{\\phi ee}$: $7\\times10^{-15}$ from S1-only analysis and $2\\times10^{-15}$ from the S2-only analysis. The S2 result is the strongest direct-detector bound on a solar electron-coupled scalar and translates into a Higgs-mixing bound $\\sin\\theta<7\\times10^{-10}$ for a relaxion-like scalar. It remains about a factor of three weaker than the strongest indirect bound from stellar cooling, but it probes parameter space that astrophysical limits constrain only through model-dependent plasma effects.","feed_headline":"Solar scalar–electron coupling capped at 2×10⁻¹⁵","feed_subtitle":"XENON1T's low-threshold ionization search beats earlier direct probes and approaches stellar-cooling limits.","key_machinery":"The load-bearing device is the scalar-photon matrix-element relation $|\\mathcal{M}(e\\to e+\\phi)|^2/|\\mathcal{M}(e\\to e+\\gamma)|^2 \\simeq g_{\\phi ee}^2 v_\\phi^2/(4\\pi\\alpha_{\\rm em})$, where $v_\\phi=k/\\omega$ is the scalar velocity and $\\alpha_{\\rm em}$ is the fine-structure constant. This single ratio lets the authors convert well-measured photon processes into scalar processes: electron-ion Bremsstrahlung and Compton production in the Sun, photoelectric absorption in xenon, and, through detailed balance with tabulated solar opacities, recombination and bound-bound production from heavy elements. The relation is stated as confirmed only in the non-relativistic limit, yet it carries the entire flux and cross-section calculation, so the predicted event rate scales as $g_{\\phi ee}^4$ and every bound inherits its uncertainty.","core_discovery":"The central claim is that existing liquid xenon dark matter detectors can serve as scalar telescopes: a light scalar $\\phi$ produced in the Sun with energies of order the solar core temperature (keV) reaches Earth and is absorbed by xenon electrons, depositing its full energy as an electronic recoil. The absorption rate is obtained by rescaling the photoelectric cross-section by $g_{\\phi ee}^2 v_\\phi^2/(4\\pi\\alpha_{\\rm em})$, the same ratio that governs solar production. Applying this to XENON1T and LUX, the paper finds $g_{\\phi ee}<2\\times10^{-15}$ at 90% CL (S2-only, $m_\\phi \\lesssim 0.1$ keV), with the sensitivity boosted by resonant production of scalars at the solar plasma frequency $\\omega\\approx\\omega_p\\approx 0.1$–$0.3$ keV. In relaxion/Higgs-portal language this is $\\sin\\theta<7\\times10^{-10}$ for a standard electron-Higgs coupling, and it is a direct noble-liquid bound that comes within an order of magnitude of stellar-evolution limits.","pith_inferences":["If the matrix-element ratio holds at sub-keV energies in a plasma, the same rescaling applies to any electron target, so semiconductor or gaseous detectors with thresholds below 100 eV could test the resonant flux peak and potentially surpass stellar-cooling limits.","The strong sensitivity of the S2 bound to the resonant peak implies that the exact solar plasma frequency profile, not just the total flux, controls the limit; a refined solar model with position-dependent Debye screening would be a direct test of the calculation.","The same Sun-as-source logic could be applied to scalars with muon or nucleon couplings, though the absorption signal in xenon would have to be replaced by a different detection channel.","If future low-threshold experiments see a signal consistent with the resonant peak, comparing spectral shape across two targets would distinguish a solar scalar from background or from other new particles such as axions."],"forward_implications":["The XENON1T S2-only analysis excludes $g_{\\phi ee}>2\\times10^{-15}$ for scalars below about 0.1 keV, a stronger direct limit than any previous laboratory search for solar scalars coupled to electrons.","For a relaxion or Higgs-portal singlet, the bound becomes $\\sin\\theta<7\\times10^{-10}$, which lies inside the naturalness window $\\sin\\theta\\lesssim7\\times10^{-9}(m_\\phi/100\\,{\\rm eV})$, so the experiment probes theoretically preferred parameters.","Because the event rate scales as $g_{\\phi ee}^4$, a 20 ton-year XENONnT exposure improves the coupling limit only to about $3\\times10^{-15}$, a modest gain rather than a dramatic jump.","The resonant solar production peak near $\\omega\\approx\\omega_p\\approx0.1$–$0.3$ keV is what makes the S2 analysis powerful; detectors with even lower ionization thresholds would capture more of this peak.","The remaining gap to the stellar-cooling bound ($g_{\\phi ee}<7\\times10^{-16}$) is a factor of three, so a modest improvement in exposure or background could make direct detection the leading constraint."],"supporting_citations":[{"why":"Supplies the central matrix-element ratio $|\\mathcal{M}(e\\to e+\\phi)|^2/|\\mathcal{M}(e\\to e+\\gamma)|^2 \\simeq g_{\\phi ee}^2 v_\\phi^2/(4\\pi\\alpha_{\\rm em})$, from which both solar production and xenon absorption rates are rescaled.","marker":"[16]"},{"why":"Provides the LUX 95-day, 118-kg scintillation data set and its background model for the S1-only limit.","marker":"[19]"},{"why":"Provides tabulated solar opacity data used (via detailed balance) to compute scalar production from bound-bound and recombination transitions in heavy elements.","marker":"[27]"},{"why":"Provides the solar model (density, temperature, composition profiles) over which the differential production rate is integrated to obtain the flux at Earth.","marker":"[31]"},{"why":"Computes the in-medium mixing of the scalar with the longitudinal photon, giving the resonant production peak and the improved stellar-cooling bound used for comparison.","marker":"[32]"},{"why":"Supplies the XENON1T 1-ton-year WIMP-search data set and background estimate for the S1-only analysis.","marker":"[36]"},{"why":"Supplies the XENON1T 22-tonne-day S2-only low-threshold data set and the assumption that sub-186 eV events are undetectable.","marker":"[42]"},{"why":"Gives the stellar-cooling bound on $g_{\\phi ee}$ against which the new direct limits are compared.","marker":"[44]"}],"fun_headline_variants":["XENON1T and LUX set direct bound on solar scalar coupling","Solar scalars: XENON1T strengthens coupling limit to 2e-15","Noble liquid detectors home in on solar relaxions","Direct solar scalar search nears stellar-evolution limits","XENON1T S2 analysis yields strongest direct solar scalar bounds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole calculation rests on a single proportionality: a scalar is emitted or absorbed exactly like a photon, with its rate reduced by the factor $g_{\\phi ee}^2 v_\\phi^2/(4\\pi\\alpha_{\\rm em})$; the paper verifies this relation only in the non-relativistic limit, so if the solar plasma or sub-keV energies alter it, the limits shift by an unquantified factor.","fun_headline_variants_meta":{"raw":{"variants":["XENON1T and LUX set direct bound on solar scalar coupling","Solar scalars: XENON1T strengthens coupling limit to 2e-15","Noble liquid detectors home in on solar relaxions","Direct solar scalar search nears stellar-evolution limits","XENON1T S2 analysis yields strongest direct solar scalar bounds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000345,"raw_usage":{"total_tokens":1914,"prompt_tokens":986,"completion_tokens":928,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":602,"completion_tokens_details":{"reasoning_tokens":849}},"tokens_in":602,"tokens_out":928,"duration_ms":9863,"temperature":1.0,"reasoning_tokens":849,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:47:46.359356+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact electron-ion Bremsstrahlung rate for scalar emission in a screened solar plasma at energies $\\omega\\approx0.1$–$0.3$ keV and compare it with the photon rate scaled by $g_{\\phi ee}^2 v_\\phi^2/(4\\pi\\alpha_{\\rm em})$; a deviation by an order of magnitude would move the S2 coupling bound by about $10^{1/4}\\approx1.8$ and could erase the claimed exclusion.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the central matrix-element ratio $|\\mathcal{M}(e\\to e+\\phi)|^2/|\\mathcal{M}(e\\to e+\\gamma)|^2 \\simeq g_{\\phi ee}^2 v_\\phi^2/(4\\pi\\alpha_{\\rm em})$, from which both solar production and xenon absorption rates are rescaled."},{"cited_title":"Up-dated opacities from the Opacity Project","cited_arxiv_id":"astro-ph/0410744","evidence_quote":"Provides tabulated solar opacity data used (via detailed balance) to compute scalar production from bound-bound and recombination transitions in heavy elements."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the stellar-cooling bound on $g_{\\phi ee}$ against which the new direct limits are compared."}],"review_version":1}