{"id":"221c0336-4978-420f-8731-f20aea911a9c","arxiv_id":"2507.16994","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":15,"one_line_summary":"Neural network quantum states, extended to include Lambda hyperons, reproduce hypernuclear separation energies to within roughly 9% and predict the observed proton-radius shrinkage in 7ΛLi.","lead":"This paper calculates the binding energies, sizes, and internal densities of light hypernuclei, nuclei containing a strange Lambda particle, using a neural-network-based quantum simulation. It is a step toward reliable ab initio predictions for medium-mass hypernuclei and for the role of strange matter in neutron stars.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The few-percent predictive claim is not yet secured: the fixed-range LO pionless Hamiltonian carries no truncation-error estimate, and the largest p-shell BΛ deviations (12ΛC, 12ΛB, 13ΛC ≈ 12%) are attributed to missing p-wave terms without being quantified by a subleading-order calculation.","rationale":"I agree with the reader's weakest-assumption analysis. The paper has real independent support: the Pfaffian-Jastrow NQS ansatz is benchmarked against SVM for hypernuclei up to A = 5, the finite-range regularization substantially improves on the naive LO pionless expectation, and the 7ΛLi core shrinkage is a qualitative prediction that does not rely on the fitted BΛ values. However, none of this closes the truncation-error gap. The paper's own diagnostics identify missing p-wave contributions as the likely source of the largest deviations in both the nucleonic and hypernuclear sectors, yet the central few-percent claim is made without estimating those contributions. The concern is not that the Hamiltonian is outside current consensus, but that the fixed-range construction deliberately gives up the standard EFT check of cutoff variation, so the theory error is simply absent. A p-wave-inclusive rerun of the non-fitted p-shell hypernuclei would settle whether the missing terms are the dominant omission; until then, the conditional verdict is appropriate and no change is needed.","tokens_in":20850,"tokens_out":6311,"duration_ms":74889,"concrete_test":"Repeat the BΛ calculation for 12ΛC, 13ΛC, and 7ΛLi with the same VMC-NQS pipeline but with leading p-wave NN and NΛ interactions added (e.g., following Ref. [80] for NN and a consistent NLO NΛ contact term), refitting the LECs to the same scattering inputs and Table II/III calibration data. If the computed BΛ values shift by more than about 1 MeV—equivalently, if the 12ΛC discrepancy decreases from roughly 12% to below 3%—the missing p-wave terms are the dominant omission and the current good agreement is not robust evidence for the LO Hamiltonian. If the shift is below the Monte Carlo errors, the truncation concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing premise is that the fixed-range, leading-order pionless EFT Hamiltonian of Eqs. (2)–(6), with Gaussian regulators fixed to scattering effective ranges (Table I) and LECs/regulators fitted to binding energies and separation energies (Tables II–III), is accurate enough to turn VMC-NQS energies into few-percent BΛ predictions. The paper states in Section IV (after Fig. 5) that the truncation error is not included because a cutoff-variation study is infeasible with the improved interaction. This matters because the same Hamiltonian already underbinds p-shell nuclei by about 9% (11B, 11C in Fig. 5), and the largest hypernuclear deviations—12ΛB, 12ΛC, and 13ΛC—are about 12% and exceed the combined statistical and model uncertainties shown in Fig. 6. The authors attribute these to missing p-wave terms but do not test that hypothesis within the EFT. Without a subleading-order calculation or another estimate of the omitted p-wave contributions, the claimed few-percent agreement for non-fitted hypernuclei is not distinguishable from an accidental outcome of the fitted central interaction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the variational Monte Carlo method with neural network quantum states (VMC-NQS), based on a Pfaffian-Jastrow ansatz with backflow, to single-Λ hypernuclei. The Hamiltonian is an improved leading-order pionless EFT with fixed finite-range regulators; the NN and NΛ two-body LECs and regulators are fitted to scattering data, and the NNN and NNΛ three-body LECs and regulators are fitted to binding and separation energies using Gaussian-process emulators of SVM and NQS energies followed by MCMC sampling. The NQS solver is benchmarked against SVM for A ≤ 5 hypernuclei (Table IV) with agreement at the few-keV level. The authors compute Λ separation energies for selected s- and p-shell hypernuclei up to 16ΛO, report an average discrepancy below 9% for systems not included in the NNΛ fit, and find a proton-radius shrinkage of δr ≈ −0.13 in 7ΛLi relative to 6Li, consistent with experiment.","tokens_in":21367,"tokens_out":7208,"duration_ms":77333,"significance":"The methodological advance is genuine and valuable: this is the first application of NQS to hypernuclear systems, it provides pure estimators for densities and radii, and it is benchmarked against SVM for light hypernuclei with excellent agreement. The Gaussian-process-plus-MCMC framework for propagating LEC uncertainties is a useful template. If the Hamiltonian truncation error were quantified, the approach would offer a scalable path to medium-mass hypernuclei. However, the paper's central few-percent predictive claim is conditional on the fixed-range leading-order Hamiltonian, and the largest p-shell deviations are about 12%, with no subleading-order estimate to show whether those deviations are a controlled EFT correction or a limitation of the model.","major_comments":[{"comment":"The central few-percent predictive claim is not secured because the truncation error of the fixed-range improved-LO Hamiltonian is not estimated. The text explicitly states after Fig. 5 that a cutoff-variation study is not feasible and that truncation error is not included. The same Hamiltonian underbinds the p-shell parents 11B and 11C by about 9%, and the largest BΛ deviations in Fig. 6 (12ΛB, 12ΛC, 13ΛC) are about 12%, larger than the combined statistical and NNΛ model uncertainties. The paper attributes these deviations to missing p-wave terms, but it does not test this hypothesis with a subleading-order calculation, an estimate of omitted p-wave contributions, or any range-sensitivity study. As a result, the agreement for non-fitted p-shell hypernuclei is not yet distinguishable from an accidental outcome of the fitted central interaction.","section":"Section IV (text after Fig. 5 and Fig. 6)"},{"comment":"The statement that agreement \"remains excellent—within theoretical and experimental uncertainties—even for the 4ΛHeS=1 and 4ΛHeS=0 states\" is not supported by the numbers in Fig. 6. For 4ΛHeS=0, the NQS value is BΛ = 2.129 ± 0.042 MeV versus the experimental 2.347 ± 0.036 MeV, a discrepancy of 0.218 MeV, which is several times the combined uncertainty and the error bars do not overlap. This system is not included in the NNΛ fit, so it is a direct test of the model, and the discrepancy should be reported quantitatively and included in the accuracy assessment.","section":"Section IV, Fig. 6"},{"comment":"Per-system accuracy for non-fitted p-shell hypernuclei is weaker than the paper's summary language. The abstract and conclusions say \"few percent\" and \"average discrepancy below 9%\", but individual non-fitted systems deviate by 12–13% (12ΛB ≈ −12.3%, 12ΛC ≈ −11.9%, 13ΛC ≈ −12.8%, 4ΛHe1 ≈ +13.1%), so the average is not representative of the worst cases. A table listing each system with its fitted/non-fitted status, BΛ, experimental value, and signed percentage deviation would make the predictive claim transparent and is needed before the summary statements can be assessed.","section":"Section IV and Section V (Conclusions)"}],"minor_comments":[{"comment":"There is a typo: \"VNC-NQS\" should read \"VMC-NQS\".","section":"Section V, Conclusions"},{"comment":"In the definition of ψ(xi), the arguments of uψ and vψ are written as (xi, xj), but xj is not defined for a single-particle orbital; they should depend only on xi.","section":"Eq. (18)"},{"comment":"The exponent for η(xi, xj) has an unbalanced parenthesis: \"a tanh(uη(xi, xj)/a) + iπ vη(xi, xj))\" contains one closing parenthesis too many.","section":"Eq. (15)"},{"comment":"The caption contains a duplicated article: \"Joint distributions of the the NNN potential\". Similar small grammar issues appear elsewhere, e.g., \"The overlaying best-fit values and corresponding energies, are indicated\" in Fig. 2.","section":"Fig. 1 caption"},{"comment":"The Gaussian-process emulator is used to replace expensive few-body and NQS calculations in the fit, but no cross-validation or estimate of GP interpolation error is reported. Even if the GP uncertainty is small relative to experimental uncertainties, stating this explicitly would strengthen the reliability of the fitted LECs and posterior distributions.","section":"Section II A"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a genuine technical advance. The Pfaffian-Jastrow VMC-NQS method is extended to single-Lambda hypernuclei, with the Lambda treated as distinguishable, and the authors carefully benchmark against SVM for A <= 5, getting few-keV agreement. Pushing to 16_Lambda O and computing radii and densities, including the 7_Lambda Li proton shrinkage, is real new content. The GP-based LEC fitting is sensible, and the paper is honest about several limitations.\n\nWhat it does well: the benchmark table is solid, the radii results are an independent prediction that goes beyond energy reproduction, and the explicit statement that the truncation error is not included (Section IV, after Fig. 5) shows the authors are not trying to hide the main caveat.\n\nThe soft spots are real, but they are the kind one expects in a first methods paper. First, there is no truncation-error estimate: the regulators are fixed by the fit, so cutoff variation is infeasible, and the few-percent claim carries no EFT error bar. Second, the abstract's \"remarkably good agreement\" leans on systems that are in the fit. The non-fitted B_Lambda values average below 9% deviation, but the p-shell cases 12_Lambda B, 12_Lambda C, and 13_Lambda C sit around 12%, larger than the combined statistical and model uncertainties. The authors attribute this to missing p-wave terms, which is plausible but not tested by a subleading-order calculation. Third, the p.Lambda scattering input is from 1968 and comes without uncertainties, and the 2.5% inflation of experimental errors on binding energies is ad hoc. Fourth, no code or data are released, which would help reproducibility but is not disqualifying.\n\nThe stress-test note is essentially correct: the few-percent predictive claim is not fully secured. But I would not call this a load-bearing flaw for the paper's main purpose. The method is credible, the benchmark proves it, and the physics conclusions are qualitative and reasonable. The next step should be quantified p-wave contributions or an explicit statement that the truncation error dominates the quoted uncertainties.\n\nThis is a paper for people working on ab initio hypernuclear structure, NQS methods, and the hyperon puzzle. It deserves a serious referee. My recommendation: send to review, require a clearer separation of fitted versus predicted results, and ask for a subleading-order estimate or an explicit bound on the truncation error.","headline":"First NQS calculations with a Lambda: a genuine method extension, well benchmarked, but the few-percent predictive claim lacks an EFT truncation error and partially rests on fitted systems.","tokens_in":21870,"tokens_out":1584,"would_cite":true,"duration_ms":20753,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V35","81V70","82B80"],"pacs":["21.80.+a","21.10.Dr","02.70.Uu"],"model":"deepseek-v4-flash","headline":"This paper extends neural-network quantum state Monte Carlo to include Λ hyperons, computing hypernuclei up to 16ΛO and matching experimental separation energies within a few percent, including the proton-radius shrinkage of 7ΛLi.","keywords":["hypernuclei","neural network quantum states","variational Monte Carlo","pionless effective field theory","Lambda separation energy","proton radius shrinkage","Gaussian process","Pfaffian-Jastrow"],"falsifier":"A concrete test: refit the NN$\\Lambda$ three-body couplings with a second, larger regulator range (for example one 20% larger) while still reproducing the same $np$ and $p\\Lambda$ scattering data and fitted $B_\\Lambda$ values, then recompute the separation energies of unconstrained hypernuclei; if ${}^{9}_{\\Lambda}\\mathrm{Be}$, ${}^{11}_{\\Lambda}\\mathrm{B}$, ${}^{12}_{\\Lambda}\\mathrm{C}$, ${}^{13}_{\\Lambda}\\mathrm{C}$, and ${}^{15}_{\\Lambda}\\mathrm{N}$ shift by more than the few-percent agreement claimed, the apparent accuracy is tied to the particular range rather than to the EFT expansion. An independent cross-check would be to compute ${}^{40}_{\\Lambda}\\mathrm{Ca}$ with this Hamiltonian and compare its $B_\\Lambda$ with a many-body calculation based on a chiral hypernuclear interaction, exposing whether the few-percent agreement persists at medium mass.","tokens_in":20615,"feed_emoji":"⚛️","tokens_out":20686,"duration_ms":177027,"temperature":0.7,"pith_summary":"The paper sets out to show that an ab initio, many-body description of light and medium-mass hypernuclei is achievable with a machine-learning solver. It extends variational Monte Carlo with neural network quantum states (VMC-NQS), previously applied to nucleons, to treat a $\\Lambda$ hyperon as a distinguishable particle alongside protons and neutrons, and couples this ansatz with an improved leading-order pionless effective field theory Hamiltonian. The low-energy constants of that Hamiltonian are fitted to two-body scattering data and to the binding and separation energies of a few light systems, with uncertainties propagated through Gaussian Process surrogates and Markov Chain Monte Carlo. Using this setup the authors compute ground-state energies, $\\Lambda$ separation energies, single-particle densities, and radii for $s$- and $p$-shell single-$\\Lambda$ hypernuclei up to ${}^{16}_{\\Lambda}\\mathrm{O}$. For hypernuclei not included in the fit, the predicted separation energies agree with experiment to within an average of 9%, and the calculation reproduces the experimentally observed shrinkage of the proton radius in ${}^{7}_{\\Lambda}\\mathrm{Li}$ relative to ${}^{6}\\mathrm{Li}$.","feed_headline":"Match hypernuclear separation energies within a few percent","feed_subtitle":"Neural-network quantum states, extended to Λ hyperons, match measured BΛ and confirm 7ΛLi core shrinkage.","key_machinery":"The machinery that carries the argument is the hypernuclear Pfaffian-Jastrow neural-network quantum state: a permutation-invariant Jastrow factor, now including nucleon–$\\Lambda$ pairs, multiplied by a Pfaffian that antisymmetrizes only the nucleons and by a separate single-particle orbital for the $\\Lambda$ hyperon, all built from message-passing backflow-transformed coordinates. This ansatz is paired with an improved leading-order pionless-EFT Hamiltonian whose two-body regulators and low-energy constants are fixed to reproduce $np$ and $p\\Lambda$ scattering lengths and effective ranges, and whose NNN and NN$\\Lambda$ three-body couplings are calibrated through Gaussian Process surrogates trained on high-accuracy few-body calculations, with posterior uncertainties sampled by Markov Chain Monte Carlo. Because the $\\Lambda$ is never antisymmetrized into the nuclear wave function, the variational state can treat it as a distinguishable particle while backflow and Jastrow correlations still keep the nucleon–$\\Lambda$ interactions fully coupled to the nuclear dynamics.","core_discovery":"The central claim is that a single, deliberately simple Hamiltonian — contact interactions at leading order in pionless effective field theory, with finite ranges and couplings fixed by two-body scattering data and a handful of few-body energies — is sufficient, when solved with a flexible neural-network variational wave function, to reproduce the $\\Lambda$-separation energies of hypernuclei up to ${}^{16}_{\\Lambda}\\mathrm{O}$ within a few percent of experiment and to capture the structural role of the hyperon. The key new element is the wave function: the Pfaffian-Jastrow ansatz is extended so that the $\\Lambda$, being distinguishable from nucleons, enters through its own single-particle orbital rather than through antisymmetrization, and this ansatz reproduces exact few-body Stochastic Variational Method benchmarks in light hypernuclei. The authors further identify a two-regime pattern in hyperon-induced size changes: in the lightest systems the weakly bound $\\Lambda$ forms a halo that pulls nucleons outward, whereas from ${}^{7}_{\\Lambda}\\mathrm{Li}$ onward the $\\Lambda$ occupies the central $1s$ orbital and acts as a glue that contracts the core — the proton radius of ${}^{7}_{\\Lambda}\\mathrm{Li}$ is computed to shrink by about 13% relative to ${}^{6}\\mathrm{Li}$ — an effect that fades to about a percent by ${}^{16}_{\\Lambda}\\mathrm{O}$.","pith_inferences":["The quoted uncertainties include Monte Carlo statistics and the spread of the NN$\\Lambda$ parameter posterior, but not the EFT truncation error, because the fixed regulator range makes a cutoff-variation study impossible; the few-percent agreement therefore includes an unquantified model-error component.","The Gaussian Process plus MCMC calibration pipeline would transfer almost unchanged to chiral hypernuclear interactions, making the workflow itself a reusable tool for sensitivity analysis rather than a single-Hamiltonian study.","The sharp halo-to-glue transition between $A=5$ and $A=7$ implies a specific, testable signature: hypernuclear and parent-nucleus charge radii should cross somewhere in that mass range, an observable that would be accessible to hypernuclear laser spectroscopy if such measurements become feasible.","A direct check of the paper's attribution of the largest deviations would be to confirm that the same $p$-wave correction that fixes the parent nuclei (${}^{11}\\mathrm{B}$, ${}^{11}\\mathrm{C}$) also removes the excess underbinding in ${}^{12}_{\\Lambda}\\mathrm{B}$, ${}^{12}_{\\Lambda}\\mathrm{C}$, and ${}^{13}_{\\Lambda}\\mathrm{C}$, without requiring new hypernuclear physics."],"forward_implications":["Hypernuclei not included in the fit, up to ${}^{16}_{\\Lambda}\\mathrm{O}$, have predicted $\\Lambda$ separation energies within an average of 9% of experiment, so the Hamiltonian ansatz is predictive rather than merely interpolative.","The same solver and interaction can be applied to medium-mass hypernuclei such as ${}^{40}_{\\Lambda}\\mathrm{Ca}$ and ${}^{48}_{\\Lambda}\\mathrm{Ca}$, giving access to the isospin dependence of the N$\\Lambda$ interaction; the paper states these systems have been measured and are within reach of the method's scaling.","In ${}^{16}_{\\Lambda}\\mathrm{O}$ the $\\Lambda$ localizes in the central $1s$ orbital, meaning the NN$\\Lambda$ interaction is probed near saturation density, the regime relevant to neutron-star cores and the hyperon puzzle.","The two-regime pattern in hyperon-induced radius changes — expansion in the lightest hypernuclei, contraction from ${}^{7}_{\\Lambda}\\mathrm{Li}$ onward, fading by ${}^{16}_{\\Lambda}\\mathrm{O}$ — is a systematic, mass-dependent prediction for how $\\Lambda$ particles modify nuclear sizes.","Residual underbinding in $p$-shell systems is attributed to missing $p$-wave contributions, so adding leading $p$-wave contact terms to both NN and N$\\Lambda$ interactions is the paper's stated route to removing the largest deviations."],"supporting_citations":[{"why":"Supplies the improved pionless-EFT interaction for nucleons and the fitting strategy of fixing regulators to experimental effective ranges.","marker":"[46]"},{"why":"Establishes the improved-action finite-range leading-order Hamiltonian for nuclei, the paradigm this work extends to the strange sector.","marker":"[55]"},{"why":"Provides the improved-action Hamiltonian benchmarks against the Stochastic Variational Method, including the agreement level the hypernuclear extension is expected to match.","marker":"[58]"},{"why":"Provides the experimental Λ separation energies used both to fit the NNΛ three-body couplings and to benchmark the predictions.","marker":"[62]"},{"why":"Introduces the Pfaffian-Jastrow neural-network quantum state used as the starting point of the hypernuclear ansatz.","marker":"[44]"},{"why":"Generalizes the Pfaffian-Jastrow ansatz to nuclear systems, the wave function extended here to include a Λ hyperon.","marker":"[45]"},{"why":"Supplies the Stochastic Variational Method used to generate the accurate few-body training data for the Gaussian Process fits and the benchmark for the neural-network ansatz.","marker":"[64]"},{"why":"Provides the experimental neutron–proton scattering lengths and effective ranges used to fix the NN low-energy constants.","marker":"[48]"},{"why":"Provides the experimental Λp scattering lengths and effective ranges used to fix the NΛ low-energy constants.","marker":"[49]"},{"why":"Reports the experimental observation of the proton-radius shrinkage in 7ΛLi that the calculation reproduces.","marker":"[14]"}],"fun_headline_variants":["Neural-network quantum states match hypernuclear binding energies","Quantum neural states reproduce Λ-hypernuclei separation energies","AI wave functions confirm hypernuclear core shrinkage","Neural states nail Λ-hypernuclei binding within a few percent"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a fixed-range, leading-order pionless-EFT Hamiltonian, with regulators and couplings calibrated to two-body scattering and a handful of few-body energies, is accurate enough to predict unconstrained hypernuclear separation energies to within a few percent — a premise the paper cannot test by varying the cutoff, because the regulator range is frozen by the fit, so the EFT truncation error is never explicitly included in the quoted uncertainties.","fun_headline_variants_meta":{"raw":{"variants":["Neural-network quantum states match hypernuclear binding energies","Quantum neural states reproduce Λ-hypernuclei separation energies","AI wave functions confirm hypernuclear core shrinkage","Neural states nail Λ-hypernuclei binding within a few percent"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000717,"raw_usage":{"total_tokens":3257,"prompt_tokens":1015,"completion_tokens":2242,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":2172}},"tokens_in":631,"tokens_out":2242,"duration_ms":18199,"temperature":1.0,"reasoning_tokens":2172,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:59:31.162754+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test: refit the NN$\\Lambda$ three-body couplings with a second, larger regulator range (for example one 20% larger) while still reproducing the same $np$ and $p\\Lambda$ scattering data and fitted $B_\\Lambda$ values, then recompute the separation energies of unconstrained hypernuclei; if ${}^{9}_{\\Lambda}\\mathrm{Be}$, ${}^{11}_{\\Lambda}\\mathrm{B}$, ${}^{12}_{\\Lambda}\\mathrm{C}$, ${}^{13}_{\\Lambda}\\mathrm{C}$, and ${}^{15}_{\\Lambda}\\mathrm{N}$ shift by more than the few-percent agreement claimed, the apparent accuracy is tied to the particular range rather than to the EFT expansion. An independent cross-check would be to compute ${}^{40}_{\\Lambda}\\mathrm{Ca}$ with this Hamiltonian and compare its $B_\\Lambda$ with a many-body calculation based on a chiral hypernuclear interaction, exposing whether the few-percent agreement persists at medium mass.","supporting_citations":[{"cited_title":"Chart of Hypernu- clides — Hypernuclear Structure and Decay Data,","cited_arxiv_id":null,"evidence_quote":"Provides the experimental Λ separation energies used both to fit the NNΛ three-body couplings and to benchmark the predictions."},{"cited_title":"Suzuki and K","cited_arxiv_id":null,"evidence_quote":"Supplies the Stochastic Variational Method used to generate the accurate few-body training data for the Gaussian Process fits and the benchmark for the neural-network ansatz."}],"review_version":1}