{"id":"c27765c9-e325-43b1-866c-b3418aa92b9f","arxiv_id":"1908.01934","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Chiral NLO S=-2 interactions, processed through G-matrix and local-density approximation, produce a Xi-nucleus potential with an attractive surface that binds Xi- in 12C and 14N at energies close to KEK emulsion values, though the (K-,K+) quasi-free peak is predicted too low in energy.","lead":"This paper predicts the force felt by a Xi hyperon inside atomic nuclei, using a chiral effective field theory interaction. The result gives a repulsive core and an attractive surface, which can bind Xi- particles in carbon and nitrogen, and it offers predictions for upcoming kaon experiments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The attractive surface pocket and the claimed Xi- binding energies rest on a low-density interpolation that the paper itself shows is ambiguous; recomputing Table II with the kF-spline is required before 'conformable' is robust.","rationale":"The reader's weakest assumption and this concern coincide: the bound states live in the surface region where the G-matrix potentials are interpolated, and the two natural spline variables give different attractions. The paper is honest about this limitation, but it remains the least secure link in the central claim. The 14N comparison is partly circular because the interaction parameters were updated using that datum [15], but the interpolation issue also affects 12C and the shape of the potential, so it is the more general concern. The 9Be spectrum gives some independent support for absolute cross-section normalization, yet the quasi-free peak sits too low in Xi energy, which the author reads as preference for a more repulsive potential; this does not settle the low-density interpolation ambiguity. Because the paper itself flags the interpolation dependence and the reader already returned CONDITIONAL, no verdict change is needed; the concern is exactly why acceptance should be conditional rather than unconditional. The concrete check is a straightforward recalculation that would either confirm the Table II numbers are robust or quantify the shift.","tokens_in":14902,"tokens_out":3584,"duration_ms":39256,"concrete_test":"Recompute the finite-nucleus potentials and Table II bound states using a cubic-spline interpolation of U_Xi(E,kF) as a function of kF (the dashed curves in Fig. 6), instead of as a function of density, keeping all G-matrix input values at kF = 0.8, 0.9, 1.07, 1.2, 1.35, and 1.5 fm^-1 and imposing U=0 at kF=0. Compare the 0s and 0p energies for 12C and 14N with Table II. If the 12C 0s energy shifts by more than about 1 MeV, or if the 14N 0s energy moves outside the 4.38 +- 0.25 MeV window, the claimed conformability is not robust to the interpolation choice.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central qualitative result -- an attractive Xi-nucleus surface pocket that binds 0s states near 4-5 MeV in 12C and 14N -- is controlled by the Xi potential at local Fermi momenta below about 0.8 fm^-1. In this region the G-matrix calculation is not self-consistent and the potential is obtained by interpolating from kF >= 0.8 fm^-1 down to U=0 at kF=0. The paper states this explicitly in Sec. III and shows in Fig. 6 that a cubic spline in density rho gives a different low-density potential than a cubic spline in kF, with the kF-spline being more attractive at low density and producing deeper 0s and 0p bound states than Table II. Because the quoted agreement with the KEK emulsion values (3.89 +- 0.24 MeV for 12C, 4.38 +- 0.25 MeV for 14N) is a quantitative claim, an interpolation ambiguity of several MeV in the surface region is load-bearing. No empirical constraint is offered to prefer one interpolation over the other, and the 9Be quasi-free peak being systematically lower than data suggests the potential may need to be more repulsive, which adds tension rather than resolving the ambiguity.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript by M. Kohno calculates Ξ single-particle potentials in symmetric nuclear matter by solving the baryon-channel coupled G-matrix equation with the chiral NLO S = −2 baryon-baryon interactions of Ref. [5] using the updated parameters of Ref. [15]. The momentum-dependent nuclear-matter potential is converted to an energy-dependent coordinate-space optical potential for finite nuclei via a local-density approximation followed by Gaussian folding with β = 1 fm, and is then parametrized as a sum of attractive and repulsive Woods-Saxon terms. The resulting Ξ− potentials in 12C and 14N yield 0s bound states at about −4 and −5 MeV, respectively, which the paper argues are conformable with KEK emulsion binding-energy determinations, and also produce 0p shifts of the atomic states. The same parametrized potentials are used in a semiclassical distorted-wave calculation of K+ spectra for (K−,K+) Ξ production on 9Be and 12C; the 9Be spectrum reproduces the absolute magnitude of the data but places the quasi-free peak at a lower K+ momentum (higher Ξ energy) than observed, while the 12C spectrum is presented as a prediction for upcoming KEK data.","tokens_in":15261,"tokens_out":7687,"duration_ms":106876,"significance":"If the results are robust, this would be one of the first nuclear-medium Ξ potentials derived from chiral NLO S = −2 interactions, connecting the microscopic interaction to emulsion bound states and (K−,K+) spectra. The partial-wave decomposition in Fig. 2 and the identification of the T = 1 3S1 baryon-channel coupling as the source of weak attraction are informative, and the paper is transparent about its approximations: the angle-averaged Pauli operator, the cutoff factor, the on-shell treatment, and the local-density convolution. The parametrized potentials in Table I are a practical output for future experimental comparisons. However, the central quantitative claim of agreement with the KEK binding energies currently rests on an unvalidated low-density interpolation, and part of the 14N agreement follows from the input parametrization having been tuned to that datum. These issues limit the conclusiveness of the 'conformable' statement and require quantification before the central claim can be accepted.","major_comments":[{"comment":"The quantitative claim of agreement with the KEK emulsion binding energies rests on the extrapolation of U_Ξ(E,k_F) below k_F ≈ 0.8 fm^-1, where the G-matrix calculation is not self-consistent and the potential is obtained by interpolating computed points to U=0 at k_F=0. The paper itself shows in Fig. 6 that a cubic-spline interpolation in density ρ differs from one in k_F, the latter being more attractive at low densities, and states that with the k_F interpolation the 0s and 0p bound states appear at lower energies than in Table II. Because the bound-state energies (e.g., 14N 0s at −5.40 MeV, close to the experimental 4.38 ± 0.25 MeV) are presented as a main result, the sensitivity of Table II to the interpolation choice must be quantified. Please recompute the bound states with the k_F interpolation and report both sets of energies, or otherwise justify the choice of the density interpolation as the physical one.","section":"Sec. III, Fig. 6 and Table II"},{"comment":"The parameters of the chiral NLO S = −2 interactions were updated in Ref. [15] by considering the experimental evidence for a Ξ− bound state in 14N (Nakazawa et al. [9]). The present calculation uses these parameters and then finds a 14N 0s state at about −5.4 MeV, which the paper describes as 'nearly matching' the same experimental value. This comparison is not independent: the agreement is partially built into the input interaction. The manuscript should explicitly acknowledge this circularity and base its claim of conformability on the 12C states and the (K−,K+) spectra, treating the 14N comparison as a consistency check rather than an independent validation.","section":"Sec. I and Sec. III.B"},{"comment":"The Gaussian resolution function in Eq. (10) as written is not normalized: its integral over E is π/ln2 ≈ 4.53, not 1. The correct normalization factor for the kernel exp(−ln2 (E/∆E)^2) is (1/∆E)√(ln2/π), not 1/(∆E√(ln2/π)). Since the 9Be spectrum is compared with data 'without any multiplicative factor' in Sec. IV.A, this normalization error directly affects the claimed reproduction of the absolute cross section in Fig. 7. Please correct Eq. (10), rerun the spectra, and confirm that the absolute normalization remains consistent with the experimental data.","section":"Sec. IV, Eq. (10)"},{"comment":"The computed quasi-free peak for 9Be sits at a lower K+ momentum (higher Ξ energy) than the data, and the paper concludes that 'the ChEFT potential may need more repulsive character.' This introduces a tension with the attractive surface potential that generates the bound states in Table II. The author should discuss whether the additional repulsion suggested by the 9Be spectrum would materially change the predicted 12C and 14N bound-state energies, and ideally quantify this with a test case, since the same potential is used for both the bound-state and the quasi-free analyses.","section":"Sec. IV.A, Fig. 7"}],"minor_comments":[{"comment":"The word 'Physiks' appears twice and should be 'Physics'.","section":"Abstract and Introduction"},{"comment":"Several entries appear to contain typographical errors: for 12C, 'r1 = {min(2.86−0.0008, 2.45+0.008E)' should likely be 'min(2.86−0.0008E, 2.45+0.008E)'; for 14N, 'a1 = 0.55 + 0.00038' should likely be '0.55 + 0.00038E'; similar missing E factors may affect other lines and should be checked.","section":"Table I"},{"comment":"In the sentence 'It is possible for the NLO ChEFT Ξ potential to generate a 0s Ξ0 bound state in 12C and 14N, but no 0p bound state exists', the phrase 'no 0p bound state exists' refers only to Ξ0, since Table II lists a Ξ− 0p state. The wording should be made unambiguous.","section":"Sec. III.B"},{"comment":"Please state explicitly how the momentum resolution (∆p/p)_K+ = 1% is converted to the energy resolution ∆E = 6.1 MeV used in Eq. (10), so that readers can reproduce the smearing.","section":"Sec. IV.A"},{"comment":"The captions contain '∆E−2 MeV', which should read '∆E = 2 MeV'.","section":"Captions of Figs. 8 and 9"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a standard and largely transparent G-matrix calculation, and the author is clearly capable of addressing the technical points. The main substantive issue is the unquantified low-density interpolation ambiguity, which affects the central bound-state claim; this should be resolved by a sensitivity study. The normalization error in Eq. (10) is a concrete, easily correctable mistake but one that touches on a quantitative comparison. There is also a partial circularity in the 14N comparison that should be acknowledged. None of these issues requires rejection; they are fixable within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take on arXiv:1908.01934. The genuinely new pieces are the finite-nucleus LDA potentials and the 0s/0p bound-state energies for 12C and 14N from the chiral NLO S=-2 interaction, plus the first realistic (K-,K+) spectrum on 9Be. The SNM Xi potential is a re-evaluation of Haidenbauer-Meissner with a different self-consistency choice, so not new by itself.\n\nWhat the paper does well: it is unusually transparent. The G-matrix and SCDW machinery is standard, and the author explicitly lists the approximations (angle-averaged Pauli, cutoff, on-shell, Gaussian folding). More importantly, Sec. III includes a frank remark that below kF ~ 0.8 fm^-1 the G-matrix is not self-consistent and the low-density potential has to be interpolated, and Fig. 6 shows the density-spline and kF-spline diverge there. The 9Be spectrum is compared with data without a normalization fudge, and the author admits the quasi-free peak sits at lower Xi energy than the data, which points toward a more repulsive potential. That is honest.\n\nThe soft spots are the ones the paper half-admits. First, the 14N 'conformable' comparison is mildly circular, since the NLO parameters in Ref. [15] were updated to give the 14N bound state. Second, the bound 0s states are controlled by exactly the low-density region where the interpolation ambiguity lives. If the kF-spline is used, the states get deeper; Table II changes by an unknown amount, likely a few MeV. The paper says this but does not quantify it, and then continues to use 'conformable' in the abstract and conclusions. The 9Be data add tension in the same direction: a more repulsive potential would push the bound states up, not down. So I would not treat the quoted -4.2 and -5.4 MeV as solid predictions.\n\nNone of this is fatal. The qualitative picture—an attractive surface pocket from baryon-channel coupling that can bind Xi-—is probably robust, and the paper gives concrete cross-section predictions for the upcoming KEK 12C analysis. The citation pattern is clean.\n\nI'd send this to a serious referee. The main request would be a revised Table II with the kF-spline interpolation and a more careful discussion of the circularity. The paper deserves publication after that. Bring it to a reading group if anyone works on strangeness; otherwise it's a solid but narrow contribution.","headline":"Useful exploratory Xi-hypernucleus calculation, but the claimed binding-energy agreement is partly circular and partly interpolation-dependent.","tokens_in":15789,"tokens_out":3899,"would_cite":true,"duration_ms":37703,"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 paper claims that the NLO chiral effective field theory interaction in the strangeness S=-2 sector, run through a coupled-channel G-matrix in nuclear matter and a Gaussian-folded local-density approximation for finite nuclei, produces…","keywords":["Xi hyperons","chiral effective field theory","strangeness -2 baryon interactions","Xi-nucleus potential","Brueckner G-matrix","local density approximation","Xi hypernuclear bound states","(K-,K+) production"],"falsifier":"A high-resolution measurement of the 12C(K-,K+)Xi- spectrum near threshold would settle it: the calculation predicts a bump from the 0s bound state around -4 MeV, whereas a null or purely repulsive potential shows no such bump; likewise, if the quasi-free peak in 9Be(K-,K+) is measured at the position the data already indicate rather than at the lower Xi energy this calculation gives, the repulsive-core/attractive-surface picture would need revision.","tokens_in":14667,"feed_emoji":"⚛️","tokens_out":11218,"duration_ms":105408,"temperature":0.7,"pith_summary":"The paper tries to show that the same next-to-leading-order chiral effective field theory (ChEFT) interaction used for baryons with strangeness -2, when evaluated in nuclear matter and then mapped to finite nuclei, yields a Xi-nucleus potential with a repulsive core and an attractive surface at low Xi energy. That surface attraction is enough to support shallow Xi- bound 0s states in 12C and 14N, with energies around -4 and -5 MeV, which match the binding energies reported from emulsion experiments. If the claim is right, the ChEFT description of the S=-2 sector can account for the observed Xi hypernuclear candidates without an artificially deep potential, and the same potential makes a concrete prediction for the near-threshold (K-,K+) spectrum on 12C that a forthcoming higher-resolution experiment can check. The paper also reports that the calculated 9Be(K-,K+) cross-section magnitude is right but the quasi-free peak sits at a lower Xi energy than the data, suggesting the potential may need extra repulsion.","feed_headline":"Surface attraction puts Xi- bound states at measured energies","feed_subtitle":"A repulsive-core, attractive-surface Xi potential yields 0s states in 12C and 14N close to emulsion binding energies.","key_machinery":"The load-bearing machinery is the coupled-channel G-matrix equation of lowest-order Brueckner theory, a resummation that handles short-range correlations, applied to the NLO ChEFT S=-2 interaction: the XiN channel is coupled nonperturbatively to Lambda-Lambda, Lambda-Sigma, and Sigma-Sigma channels, and this coupling is what turns an otherwise repulsive Xi potential into a weakly attractive one at low momentum. To reach finite nuclei, the density-dependent potential is converted to an energy-dependent coordinate-space optical potential by the local-density approximation and smoothed with a Gaussian form factor of range 1 fm to restore finite-range effects. The resulting numerically calculated potential is well reproduced by a sum of attractive and repulsive Woods-Saxon parts, and that parametrized potential is what carries the bound-state and (K-,K+) cross-section calculations.","core_discovery":"Using the NLO ChEFT S=-2 baryon-baryon interaction with the updated parameters and a 550 MeV cutoff, the paper computes Xi single-particle potentials in symmetric nuclear matter by solving the baryon-channel coupled G-matrix equation in lowest-order Brueckner theory. The resulting potential is weakly attractive at low momenta, an effect generated by coupling to Lambda-Lambda, Lambda-Sigma, and Sigma-Sigma channels, most notably in the T=1 3S1 XiN-LambdaSigma-SigmaSigma channel. Transformed to finite nuclei through a local-density approximation with a Gaussian folding of 1 fm range, the potential is repulsive in the central region and attractive in the surface at low Xi energy. It produces Xi- 0s states at about -4.2 MeV in 12C and -5.4 MeV in 14N, values the paper judges conformable with the emulsion binding energies of about 3.9 and 4.4 MeV; the alternative higher-energy assignments correspond to the 0p state, which is only slightly shifted from the atomic level. The same parametrized potential, used in a semiclassical distorted-wave calculation, reproduces the absolute magnitude of the 9Be inclusive spectrum but places the quasi-free peak lower than measured, and for 12C it predicts a near-threshold bump whose location is sensitive to the potential strength.","pith_inferences":["If the repulsive-core/attractive-surface shape is confirmed, fitting Xi-nucleus data with a single Woods-Saxon well of around 14 MeV would misrepresent the interaction: the same shallow bound states can come from a surface pocket rather than a deep central attraction.","The interpolation ambiguity at low density could be tested theoretically by computing the Xi potential directly in finite nuclei, bypassing the homogeneous-matter assumption where G-matrix self-consistency is lost.","Because the attractive T=0 s-wave contribution is absent in neutron-rich matter, this mechanism implies that Xi hyperons may be less abundant in neutron star interiors than models with universally attractive Xi potentials assume, with consequences for the equation of state.","The Gaussian-folded local-density approach used here could be extended to heavier nuclei such as 40Ca once the non-monotonic density distributions are handled, giving predictions for Xi- atomic level shifts and bound states that the paper leaves for future work."],"forward_implications":["The NLO ChEFT Xi-nucleus potential supports a Xi- 0s bound state in both 12C and 14N, so the emulsion binding energies can be explained without invoking an ad hoc deep attractive well.","No 0p bound state is predicted; the 0p state is only slightly shifted downward from the atomic level, so an experiment that can distinguish orbital assignments would separate the two possible readings of the 14N event.","The canonical Woods-Saxon potential with depth 14 MeV gives deeper 0s states than the emulsion candidates, so unless deeper states are found, that standard depth is disfavored by this calculation.","The near-threshold 12C(K-,K+) cross section is sensitive to the strength of the Xi potential, so the forthcoming higher-resolution data can discriminate between the ChEFT potential, a null potential, and a deep attractive well.","In symmetric nuclear matter the weak attraction is generated by baryon-channel coupling, while at higher densities the T=0 s-wave attraction is absent, so the same interaction implies a repulsive Xi potential in neutron-rich high-density matter."],"supporting_citations":[{"why":"supplies the NLO ChEFT baryon-baryon interaction in the S=-2 sector (Lambda-Lambda, Lambda-Sigma, Sigma-Sigma, and Xi-N channels) used as the bare two-body input.","marker":"[5]"},{"why":"provides the updated parameter set of that interaction, constrained by the Xi-14N bound-state candidate, with the 550 MeV cutoff used throughout.","marker":"[15]"},{"why":"reports the emulsion event interpreted as a deeply bound Xi-14N state with binding energy 4.38 or 1.11 MeV, the main experimental comparison for the 14N calculation.","marker":"[9]"},{"why":"reports candidate Xi-12C bound states at 3.89, 2.84, and 0.82 MeV used to compare with the calculated 12C levels.","marker":"[10]"},{"why":"provides the Lambda and Sigma single-particle potentials in nuclear matter used in the baryon propagators of the coupled G-matrix equation.","marker":"[8]"},{"why":"provides the preceding nuclear-matter calculation with ChEFT interactions that supplies the nucleon single-particle potential and saturation inputs.","marker":"[16]"},{"why":"introduces the Gaussian-folded improved local-density approximation used to convert the matter potential into finite-nucleus potentials.","marker":"[23]"},{"why":"provides the semiclassical distorted-wave method used to evaluate the (K-,K+) Xi production spectra.","marker":"[27]"}],"fun_headline_variants":["Xi hyperons get surface attraction, match emulsion binding","Surface attraction yields Xi- bound states in 12C and 14N","Chiral NLO Xi potential: repulsive core, attractive skin","Xi- states in 12C and 14N match emulsion energies","Surface pocket binds Xi- hyperons to measured levels"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The surface attraction that creates the bound states comes from the Xi potential at low nuclear density (Fermi momentum below about 0.8 $fm^{-1}$), where the G-matrix self-consistency is not reliable and the potential is filled in by interpolation; the two natural interpolation choices, as a function of density or of Fermi momentum, give different low-density attractions and therefore different bound-state energies.","fun_headline_variants_meta":{"raw":{"variants":["Xi hyperons get surface attraction, match emulsion binding","Surface attraction yields Xi- bound states in 12C and 14N","Chiral NLO Xi potential: repulsive core, attractive skin","Xi- states in 12C and 14N match emulsion energies","Surface pocket binds Xi- hyperons to measured levels"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000347,"raw_usage":{"total_tokens":1983,"prompt_tokens":1112,"completion_tokens":871,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":728,"completion_tokens_details":{"reasoning_tokens":783}},"tokens_in":728,"tokens_out":871,"duration_ms":7893,"temperature":1.0,"reasoning_tokens":783,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:59:17.143549+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-resolution measurement of the 12C(K-,K+)Xi- spectrum near threshold would settle it: the calculation predicts a bump from the 0s bound state around -4 MeV, whereas a null or purely repulsive potential shows no such bump; likewise, if the quasi-free peak in 9Be(K-,K+) is measured at the position the data already indicate rather than at the lower Xi energy this calculation gives, the repulsive-core/attractive-surface picture would need revision.","supporting_citations":[{"cited_title":"Haidenbauer, U.-G","cited_arxiv_id":null,"evidence_quote":"supplies the NLO ChEFT baryon-baryon interaction in the S=-2 sector (Lambda-Lambda, Lambda-Sigma, Sigma-Sigma, and Xi-N channels) used as the bare two-body input."},{"cited_title":"Haidenbauer and U.-G","cited_arxiv_id":null,"evidence_quote":"provides the updated parameter set of that interaction, constrained by the Xi-14N bound-state candidate, with the 550 MeV cutoff used throughout."},{"cited_title":"Nakazawa et al","cited_arxiv_id":null,"evidence_quote":"reports the emulsion event interpreted as a deeply bound Xi-14N state with binding energy 4.38 or 1.11 MeV, the main experimental comparison for the 14N calculation."},{"cited_title":"Aoki et al","cited_arxiv_id":null,"evidence_quote":"reports candidate Xi-12C bound states at 3.89, 2.84, and 0.82 MeV used to compare with the calculated 12C levels."},{"cited_title":"Kohno, Phys","cited_arxiv_id":null,"evidence_quote":"provides the preceding nuclear-matter calculation with ChEFT interactions that supplies the nucleon single-particle potential and saturation inputs."},{"cited_title":"Jeukenne, A","cited_arxiv_id":null,"evidence_quote":"introduces the Gaussian-folded improved local-density approximation used to convert the matter potential into finite-nucleus potentials."},{"cited_title":"Hashimoto, M","cited_arxiv_id":null,"evidence_quote":"provides the semiclassical distorted-wave method used to evaluate the (K-,K+) Xi production spectra."}],"review_version":1}