{"id":"ca441aec-1d80-4f91-bf8a-020c2cc42e39","arxiv_id":"2501.02826","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A relativistic Brueckner-Hartree-Fock calculation with leading-order covariant chiral hyperon-nucleon and nucleon-nucleon forces reproduces the empirical Lambda single-particle potential in nuclear matter.","lead":"This paper computes how a lambda hyperon feels inside nuclear matter using a relativistic many-body method with a leading-order chiral interaction. It finds that the computed binding potential matches the empirical value without requiring higher-order corrections, a first for this framework.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Neglecting the Σ in-medium self-energy in the coupled-channel G-matrix is load-bearing; a moderate Σ potential could move U_Λ out of the empirical band.","rationale":"The paper is a careful first application of LO covariant chiral YN and NN forces in RBHF. The agreement with the empirical U_Λ is a genuine prediction, since the LECs are fitted only to free-space YN scattering data (36 points, χ²≈15.8) and not to hypernuclear binding energies. The cutoff variation at saturation density yields U_Λ(p=0) from -25.1 to -29.3 MeV, with the 700 MeV result inside the empirical -27 to -30 MeV band. The use of a microscopic interaction with no in-medium parameters is a strength, though the absence of code/data and LEC uncertainties weakens reproducibility. The most load-bearing assumption is the neglect of the Σ in-medium self-energy in the coupled-channel Bethe–Goldstone equation (5). Even though the Σ density is zero, the ΣN intermediate channel is essential; the ΛN–ΣN transition contributes substantially to the Λ potential. Under the continuous choice, intermediate-state energies should include self-consistent single-particle potentials for all active baryons; omitting U_Σ for the Σ is an internal inconsistency. A repulsive Σ potential, as suggested by phenomenological nuclear-matter studies, would raise the ΣN threshold and reduce the coupled-channel attraction, potentially shifting U_Λ by several MeV. Since the LO(550) result is already outside the empirical band, a few MeV shift in LO(700) could invalidate the agreement. The paper does not quantify this effect. The proposed sensitivity test—scanning constant U_Σ values—would settle whether this approximation is benign. I considered other concerns (LEC uncertainties, cutoff dependence at higher densities), but these are secondary because the central claim is made at saturation density, where the cutoff band overlaps the empirical range. My verdict remains conditional, pending the Σ self-energy check and ideally artifact release.","tokens_in":16796,"tokens_out":7328,"duration_ms":75705,"concrete_test":"Re-run the RBHF self-consistency for ΛF=700 MeV and k_F=1.35 fm^-1 including a Σ single-particle potential in the intermediate-state energy denominator of Eq. (5). As a fast diagnostic, repeat the calculation with U_Σ fixed to constant values of -20, 0, +20, and +40 MeV, and report U_Λ(p=0) for each. If U_Λ changes by more than ~3 MeV compared to the paper's -29.3 MeV, the neglect of Σ medium effects is load-bearing and the empirical agreement is not robust; if the change is below ~1 MeV, the assumption is safe. A full self-consistent treatment (iterating U_Σ from the ΣN G-matrix) would be the definitive check.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that LO covariant chiral YN+NN forces reproduce the empirical U_Λ(p=0) at saturation density—depends on the treatment of ΣN intermediate states in Eq. (5). The text states (Sec. 2.1): 'the Σ hyperon density is assumed to be zero; thus, the in-medium effects on the Σ baryon are not considered.' However, the Σ appears as an intermediate channel in the Bethe–Goldstone equation, and under the adopted continuous choice, the energy denominator E_B5B6 should include the self-consistent Σ single-particle potential U_Σ. Omitting U_Σ (using free Σ energies) violates the continuous choice and shifts the ΛN–ΣN coupling strength. Because this coupling contributes significantly to the attraction/repulsion balance in U_Λ (the Jul05 model shows a strong ΛN–ΣN tensor component), even a moderate Σ potential—phenomenologically repulsive in nuclear matter—could alter U_Λ by several MeV. At ΛF=550 MeV, the reported -25.1 MeV is already 1.9 MeV outside the empirical -27 to -30 MeV band; at ΛF=700 MeV, -29.3 MeV sits near the edge. A repulsive U_Σ of +20 MeV could plausibly move the LO(700) value outside the band, undermining the claimed agreement. The paper provides no estimate of this systematic uncertainty. This is an internal approximation, not an external inconsistency, but its magnitude is unquantified and it directly affects the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes the in-medium Lambda-nucleon interaction in the relativistic Brueckner-Hartree-Fock (RBHF) framework, using for the first time the leading-order covariant chiral hyperon-nucleon (YN) and nucleon-nucleon (NN) forces. The authors re-derive the LO covariant chiral YN interaction with physical baryon masses, fit its 12 low-energy constants (LECs) to the 36 YN scattering data, and then solve the coupled-channel Bethe-Goldstone equation for the LambdaN-SigmaN system to obtain the Lambda single-particle potential U_Lambda(p=0) in symmetric nuclear matter at saturation density. They report U_Lambda(0) = -25.1 MeV for cutoff Lambda_F=550 MeV and -29.3 MeV for Lambda_F=700 MeV, comparing these with the empirical value of about -27 to -30 MeV. The paper argues that, in contrast to non-relativistic chiral calculations that require higher-order two-body forces, the leading-order covariant chiral forces already describe both the scattering data and the empirical Lambda potential, and that the resulting density dependence is consistent with a repulsive high-density Lambda potential.","tokens_in":17099,"tokens_out":2860,"duration_ms":31678,"significance":"If the central claim holds, the paper provides an important step: a microscopic, chiral-symmetry-based input for relativistic in-medium hyperon-nucleon interactions, with a genuine external benchmark (the empirical U_Lambda) that is not used in the LEC fit. The cross-section agreement in Fig. 1 and the comparison with several non-relativistic and phenomenological potentials are useful and mostly clearly presented. The paper is also honest about the cutoff variation and about the assumption that the Sigma density vanishes. However, the central quantitative claim is currently supported by only two cutoff points, with no propagated LEC uncertainties and no assessment of the omitted Sigma in-medium self-energy; these issues need to be addressed before the claim is fully established.","major_comments":[{"comment":"The treatment of the Sigma hyperon as a free intermediate state is a load-bearing approximation. The text states that the Sigma density is assumed to be zero and therefore in-medium effects on the Sigma are not considered. However, in the coupled-channel Bethe-Goldstone equation the intermediate SigmaN state enters through the energy denominator E_{B5B6}, and under the adopted continuous choice this denominator should contain the self-consistent Sigma single-particle potential U_Sigma. Omitting U_Sigma changes the LambdaN-SigmaN coupling strength, which contributes significantly to the attraction/repulsion balance in U_Lambda (the Jul05 model, for example, shows a strong LambdaN-SigmaN tensor component). A moderate repulsive U_Sigma, as phenomenologically expected in nuclear matter, could plausibly shift the reported U_Lambda by several MeV. At Lambda_F=550 MeV the result -25.1 MeV is already outside the empirical -27 to -30 MeV band, and at Lambda_F=700 MeV the value -29.3 MeV sits near the edge, so this systematic uncertainty directly affects the central claim. The authors should at least estimate the sensitivity by repeating the calculation with a simple phenomenological U_Sigma or by providing a quantitative argument for its smallness.","section":"Sec. 2.1, Eq. (5)"},{"comment":"The LEC fit uncertainties are not propagated to U_Lambda. Table 1 lists the 12 LECs and chi^2 values for two cutoffs, but gives no uncertainties or covariance matrix, and the quoted U_Lambda values in Table 2 are presented as single numbers with only the cutoff variation as an uncertainty estimate. Since the central claim is agreement with the empirical range, the authors should either propagate the fit uncertainties through the G-matrix calculation or, failing that, state clearly that no statistical uncertainty is included and provide an estimate based on, e.g., resampling or alternative fit strategies. Without this, the apparent agreement at Lambda_F=700 MeV cannot be distinguished from a fluctuation of the LEC fit.","section":"Sec. 3, Table 1 and Table 2"},{"comment":"The cutoff dependence itself is large enough to matter for the claim: U_Lambda(0) changes from -25.1 MeV at Lambda_F=550 MeV to -29.3 MeV at Lambda_F=700 MeV, and the partial-wave decomposition shows that the 3S1+3D1 contribution changes from -23.5 MeV to -12.2 MeV while the 1S0 contribution changes by about 5 MeV. The authors note that the cutoff dependence is a lower bound on theoretical uncertainties, but they do not discuss whether a larger cutoff range or a different regulator form would move the result outside the empirical band. A more systematic regulator study, or at least a discussion of the expected higher-order corrections, would strengthen the central claim.","section":"Sec. 3, Fig. 2 and Table 2"}],"minor_comments":[{"comment":"The rows labeled NonRel.-NLO19(500) and NonRel.-NLO19(650) list total potentials as 39.3 and 29.2 MeV without a minus sign, which is inconsistent with the text and with the empirical comparison; these should presumably be negative values and should be corrected.","section":"Table 2"},{"comment":"The figure caption says the red line represents the Lambda density used in Fig. 2, but the plotted curves are not labeled in the figure itself; adding a legend or explicit labels would improve readability.","section":"Fig. 5"},{"comment":"Several reference titles contain typographical errors: Ref. [76] has 'Hatree-Fock' instead of 'Hartree-Fock' and 'Bruckner' in the title, and Ref. [32] uses 'Bruckner' instead of 'Brueckner'; these should be corrected.","section":"References"},{"comment":"In Eq. (7), the expression for U_Lambda(p_Lambda) as a matrix element of U is clear, but the notation U_Lambda^S and U_Lambda^0 is introduced without explicitly defining them as the scalar and timelike-vector components used in Eq. (2); a sentence connecting the two would improve clarity.","section":"Sec. 2.1, Eq. (7)"},{"comment":"The conclusion states that the Lambda single-particle potential 'can serve as a crucial input' for hypernuclear structure and neutron star studies; it would be helpful to mention that the present calculation uses a small Lambda fraction and that the density dependence shown in Fig. 2 is for vanishing Lambda density, so extrapolation to neutron-star conditions requires additional assumptions.","section":"Sec. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of Physics Letters B and presents a genuinely new application of covariant chiral YN forces in RBHF. The main concern is not circularity—the empirical U_Lambda is an external benchmark—but rather the unquantified systematic effect of the omitted Sigma self-energy and the lack of LEC uncertainty propagation. These are fixable in a revision, so I recommend major revision rather than rejection. I also note that the manuscript's title and abstract emphasize agreement with the empirical value, but the LO(550) point lies outside the quoted band; the revision should present the uncertainty budget more carefully."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is the first RBHF study that uses LO covariant chiral YN and NN forces together, and it re-derives the YN interaction with physical baryon masses rather than an average mass. Second, the headline result—that LO alone gives U_Lambda consistent with the empirical -27 to -30 MeV—is not as clean as the abstract implies. The cutoff band spans -25.1 (Lambda_F=550) to -29.3 (Lambda_F=700) MeV, so only the upper end actually sits inside the empirical window.\n\nWhat the paper does well: the LECs are fitted to 36 YN scattering data points, not to hypernuclear binding energies or U_Lambda, so the U_Lambda comparison is a genuine external benchmark. The cross-section description is reasonable and comparable to phenomenological potentials. The partial-wave decomposition and comparison with Jul94, NSC97f, and non-relativistic chiral results are useful. The paper is honest and clearly written.\n\nThe soft spots are real but not disqualifying. The biggest one is the treatment of the Sigma. The coupled-channel G-matrix includes SigmaN intermediate states, but the Sigma self-energy is set to zero because the Sigma density is zero. With the continuous choice, the intermediate-state energies should contain the self-consistent U_Sigma. The stress-test concern is fair: the LambdaN-SigmaN coupling contributes to the attraction/repulsion balance, and a repulsive U_Sigma could shift U_Lambda by several MeV. The authors give no estimate of this systematic uncertainty. Given that the Lambda_F=550 result is already 1.9 MeV outside the empirical band, this could matter for the central claim.\n\nSecond, only cutoff variation is quoted as uncertainty. The LEC fit errors are not propagated. A chi-square of ~16 for 36 data points is okay, but the covariance information would give a more honest error bar. Third, \"Data will be made available on request\" is not adequate for a calculation-heavy paper like this.\n\nThe unexplained cutoff dependence at higher densities is also a loose end, though not central to the saturation-density claim.\n\nBottom line: this is a credible first step in an established program. The empirical comparison is a real prediction, and the paper deserves a serious referee. But the claim that LO reproduces the empirical U_Lambda needs to be re-framed: the cutoff band marginally overlaps the empirical window, and the Sigma approximation is a load-bearing simplification that should be quantified. I would engage with it, and I would send it to review, but ask for an explicit discussion of the Sigma self-energy and a release of the calculation details.","headline":"A first, useful RBHF calculation with LO covariant chiral YN+NN forces, but the U_Lambda agreement with the empirical band is less clean than claimed and the neglected Sigma self-energy is a real unquantified systematic.","tokens_in":17657,"tokens_out":3202,"would_cite":true,"duration_ms":29360,"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 a relativistic many-body calculation using only leading-order covariant chiral hyperon-nucleon and nucleon-nucleon forces reproduces the empirical Lambda single-particle potential of about -27 to -30 MeV at nuclear…","keywords":["Lambda hypernuclei","hyperon-nucleon interaction","relativistic Brueckner-Hartree-Fock","covariant chiral effective field theory","Lambda single-particle potential","Bethe-Goldstone equation","neutron star equation of state","nuclear saturation density"],"falsifier":"Repeat the relativistic Brueckner-Hartree-Fock calculation with a nonzero Sigma density and a self-consistent Sigma self-energy; if the Lambda single-particle potential at zero momentum leaves the empirical band of -27 to -30 MeV at saturation, the central agreement depends on that assumption. A tighter empirical determination of the Lambda potential depth, for example from a systematic analysis of Lambda hypernuclear binding energies, would also settle it.","tokens_in":16566,"feed_emoji":"⚛️","tokens_out":10636,"duration_ms":94159,"temperature":0.7,"pith_summary":"The paper aims to show that the depth of the Lambda single-particle potential in nuclear matter can be obtained from a relativistic many-body calculation using only leading-order covariant chiral hyperon-nucleon and nucleon-nucleon forces. This matters because that potential governs hypernuclear binding, Lambda flow in heavy-ion collisions, and the density at which hyperons appear in neutron stars. In nonrelativistic Brueckner-Hartree-Fock calculations, reproducing the empirical depth usually requires higher-order two-body chiral forces; the authors claim the relativistic treatment removes that need. Their result at saturation density is about -25 MeV for cutoff 550 MeV and -29 MeV for cutoff 700 MeV, overlapping the empirical band of -27 to -30 MeV.","feed_headline":"Lambda potential hits empirical depth with leading-order chiral forces","feed_subtitle":"Relativistic many-body calculation reproduces the -27 to -30 MeV Lambda depth without higher-order corrections.","key_machinery":"The carrying mechanism is the relativistic Brueckner-Hartree-Fock self-consistency loop: the in-medium Dirac equation defines baryon spinors through scalar and vector self-energies; the Bethe-Goldstone equation with an angle-averaged Pauli operator and a continuous intermediate spectrum builds the in-medium LambdaN G-matrix with the LambdaN-SigmaN coupled channel; and the G-matrix defines the Lambda single-particle potential whose scalar and vector parts feed back into the Dirac equation. The bare interaction is the leading-order covariant chiral hyperon-nucleon potential, contact terms plus one-pseudoscalar-meson exchange, regularized with an exponential cutoff from 550 to 700 MeV and re-derived with physical Lambda, Sigma, and nucleon masses.","core_discovery":"The central discovery claimed is that a re-derived leading-order covariant chiral hyperon-nucleon interaction with physical baryon masses, solved through the relativistic Brueckner-Hartree-Fock G-matrix with coupled LambdaN-SigmaN channels, yields a Lambda single-particle potential at nuclear saturation density consistent with the empirical value of about -27 to -30 MeV. The same interaction also reproduces the low-energy Lambda-proton and Sigma-proton scattering cross sections. Because the nucleon side uses the matching leading-order covariant chiral nucleon-nucleon force, the calculation is internally consistent, and the agreement is reached without the higher-order two-body chiral forces that nonrelativistic calculations typically require. The authors interpret this as evidence that the relativistic framework already captures much of the in-medium repulsion.","pith_inferences":["The paper does not pursue this, but letting the Sigma density be nonzero in the same coupled-channel loop would show how much of the leading-order agreement relies on a cancellation between LambdaN and SigmaN contributions.","An implication not drawn in the paper is that a direct comparison of the in-medium G-matrix from the relativistic and nonrelativistic frameworks using the same bare chiral interaction would isolate whether the relativistic scalar-vector cancellation is what removes the need for higher-order two-body forces.","The momentum dependence of the Lambda potential is a sharper discriminator than its depth: the 550 MeV cutoff result changes sign near a Lambda momentum of about 1.4 inverse femtometers, so Lambda directed-flow data at moderate momenta could distinguish the cutoff choices."],"forward_implications":["The empirical Lambda depth is reproduced with leading-order two-body chiral forces, so the relativistic framework weakens the usual motivation for adding next-to-leading-order two-body hyperon-nucleon terms in in-medium studies.","The predicted density dependence, with repulsion setting in at moderate densities and strengthening at high density, pushes Lambda onset to higher density and therefore favors massive neutron star formation.","The extracted scalar and vector Lambda self-energies can serve as microscopic inputs for covariant density functionals and for hypernuclear structure calculations.","The cutoff spread from 550 to 700 MeV sets a lower bound on theoretical uncertainty, with the largest spread at high density and high Lambda momentum."],"supporting_citations":[{"why":"Supplies the leading-order covariant chiral nucleon-nucleon interaction and the nucleon scalar and vector potentials used as the background in the relativistic Brueckner-Hartree-Fock loop.","marker":"[76]"},{"why":"Introduced the leading-order covariant chiral hyperon-nucleon contact and one-meson-exchange potentials that this work re-derives with physical baryon masses.","marker":"[80]"},{"why":"Provided the average-mass version of the covariant chiral interactions whose cross sections are compared with the physical-mass results.","marker":"[82]"},{"why":"Supplies the nonrelativistic next-to-leading-order chiral benchmark whose in-medium Lambda potential is compared in Table 2 and Figure 2.","marker":"[41]"},{"why":"Supplies the nonrelativistic next-to-leading-order chiral interaction used as a comparison for the saturation-depth result.","marker":"[47]"},{"why":"Supplies the Nijmegen NSC97f meson-exchange potential whose Lambda single-particle potential is a comparison point.","marker":"[51]"},{"why":"Supplies the Jul05 potential whose density dependence, driven by the LambdaN-SigmaN tensor force, is compared with the covariant result.","marker":"[53]"},{"why":"Provides the meson-exchange hyperon-nucleon G-matrix formalism and the Jul94 potential used as the main relativistic Brueckner-Hartree-Fock comparison.","marker":"[32]"},{"why":"Established the metastable treatment of a small Lambda fraction and the Lambda-density dependence that this work follows.","marker":"[77]"},{"why":"Provides the experimental hyperon-nucleon scattering data used to fix the twelve low-energy constants.","marker":"[86–91]"}],"fun_headline_variants":["RBHF with leading-order chiral LambdaN nails empirical depth","Leading-order chiral force matches Lambda depth in RBHF","Relativistic Lambda potential hits -30 MeV without extra terms","One chiral force fits Lambda scattering and nuclear depth"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the Sigma hyperon feels no in-medium self-energy because its density is set to zero, even though the Sigma appears as an intermediate state in the coupled LambdaN-SigmaN channel.","fun_headline_variants_meta":{"raw":{"variants":["RBHF with leading-order chiral LambdaN nails empirical depth","Leading-order chiral force matches Lambda depth in RBHF","Relativistic Lambda potential hits -30 MeV without extra terms","One chiral force fits Lambda scattering and nuclear depth"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000314,"raw_usage":{"total_tokens":1732,"prompt_tokens":844,"completion_tokens":888,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":460,"completion_tokens_details":{"reasoning_tokens":823}},"tokens_in":460,"tokens_out":888,"duration_ms":7264,"temperature":1.0,"reasoning_tokens":823,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:03:52.000879+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the relativistic Brueckner-Hartree-Fock calculation with a nonzero Sigma density and a self-consistent Sigma self-energy; if the Lambda single-particle potential at zero momentum leaves the empirical band of -27 to -30 MeV at saturation, the central agreement depends on that assumption. A tighter empirical determination of the Lambda potential depth, for example from a systematic analysis of Lambda hypernuclear binding energies, would also settle it.","supporting_citations":[{"cited_title":"Neutron-proton effective mass splitting in neutron-rich matter","cited_arxiv_id":"2304.13333","evidence_quote":"Supplies the leading-order covariant chiral nucleon-nucleon interaction and the nucleon scalar and vector potentials used as the background in the relativistic Brueckner-Hartree-Fock loop."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Nijmegen NSC97f meson-exchange potential whose Lambda single-particle potential is a comparison point."},{"cited_title":"Holzenkamp, K","cited_arxiv_id":null,"evidence_quote":"Supplies the Jul05 potential whose density dependence, driven by the LambdaN-SigmaN tensor force, is compared with the covariant result."},{"cited_title":"Reuber, K","cited_arxiv_id":null,"evidence_quote":"Provides the meson-exchange hyperon-nucleon G-matrix formalism and the Jul94 potential used as the main relativistic Brueckner-Hartree-Fock comparison."}],"review_version":1}