{"id":"d6f5ebe8-9cec-4f8d-8c76-da12a8cd4818","arxiv_id":"1908.09739","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"A first-principles coupled-cluster calculation predicts the weak form factor and neutron skin of 40Ar, validated by electron scattering data.","lead":"This paper computes how neutrinos scatter off argon-40 nuclei using a first-principles nuclear theory called coupled-cluster, and compares the results to electron scattering data. It predicts the neutron radius and neutron skin of argon-40, quantities relevant to neutrino detectors like DUNE and COHERENT.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Charge-form-factor validation is proton-dominated; the 40Ar neutron radius and skin rest on an untested correlation and could be biased despite the good F_ch agreement.","rationale":"The reader's verdict is CONDITIONAL, and my analysis agrees that the paper is solid but needs additional error quantification. I identify the most load-bearing concern as the proton-only validation of the neutron-dominated observables, which is the same underlying issue the reader flags in the weakest_assumption. However, I judge the neglect of two-body weak currents to be less critical for the quoted Rn and Rskin because those are density observables, not current-dependent; the missing two-body current affects the CEνNS cross-section prediction and the extraction of Rn from future data, but not the ab initio density result. Similarly, propagation of the experimental Rp uncertainty is a relatively minor error-bar issue that would widen the quoted range but not overturn the central claim. The more fundamental risk is that the charge form factor comparison cannot test the neutron density, and the Rn-Rp correlation method, while clever and partly supported by DFT, has not been directly validated for 40Ar. An independent many-body calculation with the same Hamiltonian would settle whether the DCE-EOM-CC open-shell treatment is reliable for the neutron radius. If the independent method agrees within the quoted uncertainty, the concern is resolved; if it does not, the paper's central Rn and Rskin predictions would need revision. This is consistent with retaining the CONDITIONAL verdict, as the requested conditions address exactly this missing validation.","tokens_in":11989,"tokens_out":14846,"duration_ms":163628,"concrete_test":"Compute the 40Ar ground state with the same NNLOsat Hamiltonian using an independent open-shell ab initio method (e.g., IM-SRG or self-consistent Green's function) and compare the point-neutron radius Rn and the weak form factor FW(q) to the DCE-EOM-CC results. If Rn differs by more than ~0.05 fm, or FW differs by more than 2% at q=100 MeV/c, the quoted ranges understate the many-body uncertainty and the proton-only validation is insufficient.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central Rn and Rskin predictions are presented as first-principles results validated by the 40Ar charge form factor, but the measured F_ch provides almost no constraint on the neutron distribution. Because G_En(0)=0, the neutron contribution to F_ch enters only through the small neutron charge form factor, so the Ottermann data mainly validates the point-proton density. The quoted Rn range (3.36-3.45 fm) and Rskin range (0.035-0.09 fm) are obtained by intersecting the Hamiltonian-spread correlation band with the experimental Rp, relying on a tight linear Rn-Rp correlation. If the isovector part of the chiral Hamiltonian (e.g., the three-nucleon force) shifts Rn without shifting Rp, the Rp anchor does not capture that uncertainty. The DFT points in Fig. 3 show vertical scatter relative to the ab initio band, indicating the correlation is not perfectly tight. The D/T-1 difference is used as the many-body truncation uncertainty, but the convergence of the open-shell DCE-EOM-CC neutron density is not separately checked against any neutron-sensitive observable. Thus the neutron-dominated weak form factor and the quoted radii could carry a systematic error invisible to the charge-form-factor validation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents coupled-cluster calculations of the coherent elastic neutrino-nucleus scattering (CEvNS) observables for 40Ar, using several chiral effective field theory Hamiltonians. The authors compute the charge form factor and compare it to electron scattering data from Ottermann et al., validate the many-body convergence between D and T-1 levels, and then predict the weak form factor, the CEvNS cross section, and the point-proton and point-neutron radii. By exploiting a strong correlation between Rp and Rn, and intersecting the computed band with the experimental Rp, they quote a neutron radius range 3.36-3.45 fm and a neutron skin thickness range 0.035-0.09 fm. The results are found to be consistent with density functional theory predictions.","tokens_in":12169,"tokens_out":2475,"duration_ms":28446,"significance":"If the predictions are reliable, this is a useful first-principles result for neutrino experiments using liquid argon (COHERENT, DUNE) and for constraining neutron distributions. The paper's strengths include the use of multiple chiral Hamiltonians calibrated to independent data, an explicit check of D vs T-1 convergence for the charge form factor, and a transparent presentation of the correlation between Rp and Rn. The consistency with density functional theory results adds confidence. However, the central validation argument rests on charge form factor data, which is proton-dominated, while the key predictions concern neutron-dominated observables; this transferability is the main correctness risk. The manuscript is honest about its reliance on one-body electroweak currents but does not quantify the associated uncertainty.","major_comments":[{"comment":"The statement 'This comparison validates the theory' is too strong for the neutron observables that follow. The electron scattering data constrain the charge form factor, which is dominated by the proton distribution because the neutron charge form factor is suppressed. The subsequent predictions for Fn, Rn, and Rskin rely on an untested transferability of the proton-sensitive validation to neutron-dominated quantities. Please either provide a quantitative argument for why the agreement in Fch constrains Fn (e.g., using the measured 48Ca neutron skin or parity-violating electron scattering constraints) or explicitly state that the neutron predictions currently rest on the assumed reliability of the chiral Hamiltonians rather than on direct validation.","section":"Results, Fig. 1"},{"comment":"The neglect of two-body electroweak currents is asserted with a reference to their expected smallness, but no estimate is given for 40Ar or for the momentum range relevant to CEvNS. At q=50-100 MeV, where the quoted spread among Hamiltonians is 2-6%, two-body currents could contribute at a comparable level. Please quantify the expected size of two-body current contributions to the weak form factor of a medium-mass nucleus, or provide a conservative uncertainty band that includes their possible effect.","section":"Interactions, 'For electroweak operators...'"},{"comment":"The quoted constraint 3.36 <= Rn <= 3.45 fm is obtained by intersecting the correlation band with the experimental Rp, but the uncertainty on the experimental Rp from Ref. [46] is not propagated into the final range, and the construction of the 'symmetric spread' band is not fully specified. The figure also shows that the DFT points scatter vertically relative to the ab initio band, which indicates that the linear Rn-Rp correlation is not exact. Please propagate the experimental Rp uncertainty and define the band construction quantitatively, or weaken the quoted ranges accordingly.","section":"Results, Fig. 3"},{"comment":"The convergence between D and T-1 levels is demonstrated for the charge form factor, but the neutron density and weak form factor are the key outputs. No convergence check is shown for the neutron-dominated observables against any neutron-sensitive data or against higher-order coupled-cluster truncations. Please state whether the D/T-1 difference for Fn and Rn was computed (and if so, report it), or explicitly list this as a limitation of the systematic uncertainty estimate.","section":"Results, Fig. 2(a) and Fig. 3"}],"minor_comments":[{"comment":"The text 'via q2 = sqrt(2E_nu M T / (E_nu - T)) approx sqrt(2 M T)' appears to contain a typo: the expression with the square root should be for q, not q^2. Please correct this to avoid dimensional inconsistency.","section":"Results, CEνNS cross section section"},{"comment":"The phrase 'The neutron-skin thickness of 40Ar40' contains a duplicated isotope label; it should read '40Ar'.","section":"Abstract"},{"comment":"The reference to 'N. Schunk, private communication' likely should be 'N. Schunck', consistent with Ref. [5].","section":"References, Ref. [47]"},{"comment":"The text states that the ab initio weak form factor agrees with density functional theory [5], but the corresponding DFT curve is not shown in Fig. 2. Please indicate where this agreement is visible or provide the comparison.","section":"Results, Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely topic and the coupled-cluster calculations appear carefully executed. My main concern is the gap between the proton-sensitive validation and the neutron-sensitive claims; this is fixable with additional analysis or more cautious wording, but it is load-bearing for the central predictions. The manuscript would also benefit from a quantitative statement on two-body currents. I recommend major revision rather than rejection because the underlying method and results are likely sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper does what it says: it computes the CEνNS weak form factor, neutron radius, and neutron skin of 40Ar from coupled-cluster theory with chiral EFT Hamiltonians, and it validates the proton side against electron scattering. The D/T-1 convergence checks and the spread over six Hamiltonians give a credible systematic range. The result that Rskin is 0.035–0.09 fm and Rn is 3.36–3.45 fm is new for 40Ar and consistent with DFT, which is reassuring.\n\nThe main soft spot, as the stress-test note says, is that the charge form factor validation is proton-dominated. G_En is tiny, so the Ottermann data pins down the point-proton density, not the neutron density. The neutron radius therefore leans on the Rp–Rn correlation across Hamiltonians plus the experimental Rp anchor. That is a legitimate way to narrow the range, but it assumes the isovector part of the Hamiltonian is captured by that correlation. If a three-nucleon force shifted Rn without moving Rp, the quoted range would miss it. The paper should say this more plainly and, ideally, show a neutron-sensitive comparison if any exists.\n\nSecond, the two-body current contribution is dismissed with one sentence and a citation. At the CEνNS-relevant low q it may indeed be negligible, but since this is a first-principles claim, a rough estimate or a citation that quantifies it for a similar nucleus would close the gap. Third, the experimental uncertainty on Rp from Angeli and Marinova is not propagated into the Rn range; the spread is Hamiltonian-driven only. This is minor because the Hamiltonians dominate the uncertainty, but it is a missing error bar.\n\nThe circularity burden is low: these Hamiltonians are fit to independent data, not to 40Ar, so the prediction is not circular. The paper is honest about what it does and does not validate.\n\nBottom line: this is a serious, useful calculation for the COHERENT, DUNE, and dark-matter communities. It deserves a careful referee. I would accept it with minor revisions asking for the two-body current estimate and the Rp uncertainty propagation, plus a sentence that the neutron radius rests on the correlation rather than direct neutron-sensitive data.\n\nWould cite: yes. Reading group: maybe—good for a nuclear-structure/neutrino interface discussion.","headline":"Solid first-principles 40Ar weak form factor with honest caveats; the neutron radius range rests partly on an untested Rp-Rn correlation, but this is a credible, useful contribution.","tokens_in":12766,"tokens_out":1919,"would_cite":true,"duration_ms":20399,"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":"Coupled-cluster calculations from chiral Hamiltonians predict 40Ar's weak form factor and a neutron skin of 0.035–0.09 fm, with proton-scattering data as validation.","keywords":["coherent elastic neutrino-nucleus scattering","coupled-cluster theory","40Ar","weak form factor","neutron skin","neutron radius","chiral effective field theory","charge form factor"],"falsifier":"Measure the 40Ar neutron radius directly, for example with parity-violating electron scattering at momentum transfers around 0.5–1.0 fm$^{-1}$, and compare with the predicted range $R_n = 3.36$–$3.45$ fm and $R_{\\mathrm{skin}} = 0.035$–$0.09$ fm; a result outside this band while the charge form factor remains reproduced would refute the assumption.","tokens_in":11753,"feed_emoji":"⚛️","tokens_out":11746,"duration_ms":110633,"temperature":0.7,"pith_summary":"The paper aims to establish that the neutron-sensitive weak form factor of 40Ar can be computed reliably from first principles, which matters because liquid-argon detectors are major targets for coherent elastic neutrino-nucleus scattering and dark matter searches. The authors use coupled-cluster theory with chiral effective field theory Hamiltonians, validating their wave functions by reproducing the measured electron-scattering charge form factor, then predicting the weak form factor, neutron radius, and neutron skin. They find a neutron radius of 3.36–3.45 fm and a neutron skin of 0.035–0.09 fm, consistent with density functional theory. They also find that the coherent scattering cross section varies only mildly across Hamiltonians, so precise measurements would be needed to constrain nuclear structure; a liquid-argon detector's neutrino signal is dominated by neutron counting.","feed_headline":"First-principles calculation sets 40Ar's neutron skin at 0.035–0.09 fm","feed_subtitle":"Because the weak form factor is nearly all neutrons, this anchors liquid-argon neutrino detectors.","key_machinery":"The load-bearing object is the weak form factor $$F_W($q^{2}$)=\\frac{N F_n($q^{2}$)-(1-4\\$sin^{2}$\\theta_W)Z F_p($q^{2}$)}{Q_W},$$ which, because $1-4\\sin^2\\theta_W\\simeq0.0457$, is almost exactly the neutron form factor $F_n(q^2)$ at the low momentum transfers relevant to coherent scattering. The wave functions come from coupled-cluster theory: a similarity transformation of the Hamiltonian, truncated at singles, doubles, and linearized triples, with a double-charge-exchange equation-of-motion operator that turns two protons of 40Ca into two neutrons and maps the closed-shell reference state onto open-shell 40Ar. The second mechanism is the strong $R_p$–$R_n$ correlation observed in the coupled-cluster results: intersecting the computed band with the known experimental proton radius narrows the neutron radius to 3.36–3.45 fm, while the correlated radii keep the skin thickness small.","core_discovery":"The paper's central claim is that coupled-cluster theory with Hamiltonians from chiral effective field theory can produce predictive weak form factors for 40Ar from first principles, without fitting anything to 40Ar neutron data. The validation is the charge form factor: the same wave functions reproduce electron-scattering measurements up to about $q=2$ fm$^{-1}$, and the electroweak operator is taken as the one-body current alone, with two-body currents argued to be negligible. Since the weak charge of the proton is suppressed, the predicted $F_W$ is essentially the neutron form factor, leading to $R_n=3.36$–$3.45$ fm and $R_{\\mathrm{skin}}=0.035$–$0.09$ fm. The paper further shows that the coherent neutrino-scattering cross section on 40Ar is only mildly sensitive to the choice of Hamiltonian, varying by 2–6% over the relevant momentum range.","pith_inferences":["If the proton-validated wave functions transfer to the neutron channel, the same pipeline should give trustworthy weak form factors for other open-shell nuclei that have electron-scattering data, which could be checked against future CEνNS spectra.","The tight $R_p$–$R_n$ covariance means one independent neutron-radius measurement on 40Ar would discriminate among the chiral Hamiltonians; the paper's own argument already shows most of the spread is a common shift in both radii.","Since the low-$q$ cross section moves by only 2–6% across Hamiltonians, CEνNS at current precision is better suited to testing Standard Model parameters, while nuclear-structure extraction would require pushing measurements to higher momentum transfers.","For argon-based dark matter detectors, the constrained weak form factor directly shrinks the uncertainty on the neutrino-floor background they must subtract."],"forward_implications":["The weak form factor $F_W(q^2)$ is essentially the neutron form factor, so the CEνNS rate on 40Ar scales with $N^2$ and is a neutron-distribution observable.","The calculation constrains the 40Ar neutron radius to $R_n = 3.36$–$3.45$ fm and the neutron skin to $R_{\\mathrm{skin}} = 0.035$–$0.09$ fm, consistent with density functional theory.","The predicted CEνNS cross section varies only about 2% at $q = 50$ MeV and 6% at $q = 100$ MeV across the Hamiltonians, so existing-level experiments probably cannot distinguish the interactions.","The first minimum of the weak form factor sits about $0.035$ fm$^{-1}$ below the charge form factor's minimum, showing that the neutron distribution extends beyond the proton distribution.","Precision CEνNS measurements on argon could, in turn, discriminate among chiral effective field theory Hamiltonians."],"supporting_citations":[{"why":"Supplies the double-charge-exchange equation-of-motion method that turns the 40Ca closed-shell reference into the 40Ar ground state.","marker":"[31]"},{"why":"Supplies the electron-scattering charge form factor data that the calculation must reproduce, and which the paper uses to validate the theory.","marker":"[33]"},{"why":"One of the chiral Hamiltonians used for the predictions; its results set part of the spread in form factors and radii.","marker":"[37]"},{"why":"The delta-full chiral potential used as a second, independent Hamiltonian for the form-factor and radius predictions.","marker":"[38]"},{"why":"Defines the family of softened chiral potentials whose spread is shown as a band in the weak form factor.","marker":"[43]"},{"why":"Provides the experimental 40Ar charge (proton) radius used to intersect the $R_p$–$R_n$ correlation and narrow $R_n$.","marker":"[46]"},{"why":"Establishes the $R_p$–$R_n$ correlation method in 48Ca that the paper transfers to 40Ar.","marker":"[50]"},{"why":"Density functional theory neutron-density results that the predicted radius and skin must be and are consistent with.","marker":"[5]"},{"why":"Cited grounds for neglecting two-body electroweak currents, supporting the one-body operators used in the form-factor predictions.","marker":"[44,45]"}],"fun_headline_variants":["First-principles 40Ar neutron skin: 0.035–0.09 fm","Coupled-cluster predicts 40Ar neutron skin for CEvNS","Chiral EFT pins 40Ar weak form factor for neutrino detectors","40Ar neutron skin from first principles, validated by electron scattering"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that matching the measured proton charge form factor, using only one-body electroweak currents, is enough to guarantee that the predicted neutron-dominated weak form factor—for which no direct 40Ar neutron-sensitive data exist—is correct.","fun_headline_variants_meta":{"raw":{"variants":["First-principles 40Ar neutron skin: 0.035–0.09 fm","Coupled-cluster predicts 40Ar neutron skin for CEvNS","Chiral EFT pins 40Ar weak form factor for neutrino detectors","40Ar neutron skin from first principles, validated by electron scattering"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000269,"raw_usage":{"total_tokens":1566,"prompt_tokens":833,"completion_tokens":733,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":449,"completion_tokens_details":{"reasoning_tokens":652}},"tokens_in":449,"tokens_out":733,"duration_ms":7222,"temperature":1.0,"reasoning_tokens":652,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:02:44.754011+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the 40Ar neutron radius directly, for example with parity-violating electron scattering at momentum transfers around 0.5–1.0 fm$^{-1}$, and compare with the predicted range $R_n = 3.36$–$3.45$ fm and $R_{\\mathrm{skin}} = 0.035$–$0.09$ fm; a result outside this band while the charge form factor remains reproduced would refute the assumption.","supporting_citations":[{"cited_title":"How robust is the n = 34 subshell closure? first spectroscopy of 52Ar,","cited_arxiv_id":null,"evidence_quote":"Supplies the double-charge-exchange equation-of-motion method that turns the 40Ca closed-shell reference into the 40Ar ground state."},{"cited_title":"Elastic electron scattering from 40ar,","cited_arxiv_id":null,"evidence_quote":"Supplies the electron-scattering charge form factor data that the calculation must reproduce, and which the paper uses to validate the theory."},{"cited_title":"Sta- tistical uncertainties of a chiral interaction at next-to- next-to leading order,","cited_arxiv_id":null,"evidence_quote":"One of the chiral Hamiltonians used for the predictions; its results set part of the spread in form factors and radii."},{"cited_title":"Jiang et al","cited_arxiv_id":null,"evidence_quote":"The delta-full chiral potential used as a second, independent Hamiltonian for the form-factor and radius predictions."},{"cited_title":"Improved nuclear matter calculations from chiral low-momentum interactions,","cited_arxiv_id":null,"evidence_quote":"Defines the family of softened chiral potentials whose spread is shown as a band in the weak form factor."},{"cited_title":"Table of experimental nu- clear ground state charge radii: An update,","cited_arxiv_id":null,"evidence_quote":"Provides the experimental 40Ar charge (proton) radius used to intersect the $R_p$–$R_n$ correlation and narrow $R_n$."},{"cited_title":"Neutron and weak-charge distributions of the 48Ca nucleus,","cited_arxiv_id":null,"evidence_quote":"Establishes the $R_p$–$R_n$ correlation method in 48Ca that the paper transfers to 40Ar."},{"cited_title":"Neutrino-nucleus coherent scattering as a probe of neu- tron density distributions,","cited_arxiv_id":null,"evidence_quote":"Density functional theory neutron-density results that the predicted radius and skin must be and are consistent with."}],"review_version":1}