{"id":"11772773-0324-458d-bfac-8a440a5d644e","arxiv_id":"hep-lat/0312030","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A next-to-next-to-leading-order chiral perturbation theory formula predicts the volume dependence of the nucleon mass with no new free parameters, and it matches Nf=2 lattice QCD data at three box sizes.","lead":"Physicists derived how the mass of the proton and neutron should drop when the simulated box of space is shrunk, using chiral perturbation theory with no new adjustable constants. Comparing that curve to lattice QCD simulations from four collaborations, they find the observed volume dependence matches, which matters for extracting the real-world nucleon mass from finite computer simulations.","discovery_kind":"unclear","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The finite-size curves in Figs. 3-5 are anchored to the largest-volume point and are never recomputed with the alternative LECs of Fit 3/Fit 4, so the 'parameter-free' agreement rests on LECs fit at m_pi up to 800 MeV without a direct stability check.","rationale":"The derivation is transparent and the Appendix cross-check against Luescher's formula is an honest independent test, including the reported factor-of-two discrepancy in the leading exponential term. That discrepancy is numerically irrelevant for the data comparison because at m_pi L ~ 4-8 the exponential term is negligible and the agreement relies on the |n| > 1 terms, as the paper itself notes in Section 6. The central claim is supported by Figs. 3-5, where the central values of the smaller-volume points align well with the anchored curves. The remaining soft spot is the bridge from the Fit 1 LECs to the finite-size curves: the paper shows LEC stability via Fit 3/Fit 4 but does not propagate it to the curves, and the visual agreement is not quantified with an error band. This is a real but not fatal gap, and the reader's CONDITIONAL verdict is appropriate. I would keep the verdict unchanged and note that the paper's Fit 3/Fit 4 provide more mitigation than the reader's weakest_assumption credits, since the LECs are already shown to be stable under a stricter mass cut; what remains undone is the final propagation step to the finite-size curves.","tokens_in":18734,"tokens_out":17021,"duration_ms":319228,"concrete_test":"Recompute the three curves in Figs. 3-5 using the LECs of Fit 3 (m0 = 0.88 GeV, c1 = -0.93 GeV^-1) and Fit 4 (m0 = 0.87 GeV, c1 = -1.11 GeV^-1) instead of Fit 1, and also vary m0 and c1 within their Fit 1 errors (m0 = 0.89 +/- 0.06 GeV, c1 = -0.93 +/- 0.05 GeV^-1). Overlay the resulting bands on the same data points. If the bands remain within the quoted data errors at the smaller-volume points (31, 36, 12, 13, 30, 35), the concern does not land; if they shift outside those errors, the parameter-free agreement is not robust to the LEC choice.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The finite-size formula Eq. (24) depends on m0 and c1 through both Delta_a (Eq. 15) and Delta_b (Eq. 23), and these are fixed by Fit 1 to ten large-volume points with m_pi up to about 800 MeV (Table 4). The paper's own Fit 3 (four smallest-mass points, m0 = 0.88 GeV, c1 = -0.93 GeV^-1) and Fit 4 (m0 = 0.87 GeV, c1 = -1.11 GeV^-1) show the LECs are reasonably stable, which is partial mitigation. But the curves in Figs. 3-5 are computed only with Fit 1 and are never recomputed with Fit 3/Fit 4 or with m0 and c1 varied within their quoted errors (m0 = 0.89(6) GeV, c1 = -0.93(5) GeV^-1). Because each curve is anchored to the largest-volume Monte Carlo point, only the shape of the L-dependence is tested, and the shape's sensitivity to m0 and c1 is not quantified. If the LECs shift within their errors, the predicted curves could move by amounts comparable to the data errors (point 31 carries +/-82 MeV), which would weaken the 'remarkable agreement' of abstract:s3. The paper acknowledges the convergence question in Section 6 but does not close this specific loop by propagating the LEC uncertainty into the finite-size curves.","agreement_with_reader":"agree"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The paper computes the volume dependence of the nucleon mass in relativistic SU(2)_f baryon chiral perturbation theory through O(p^4), obtaining the finite-size formulae of Eqs. (15), (23) and (24). The low-energy constants m0, c1 and er1 are fixed by fitting the infinite-volume O(p^4) formula (Eq. 16) to ten large-volume Nf=2 lattice points (Fit 1, Table 4), with gA, f_pi, c2 and c3 taken from phenomenology. The resulting parameter-free (given the LECs) prediction for mN(L) - mN(infinity) is compared with UKQCD/QCDSF data and with published CP-PACS and JLQCD values at three pion masses (Figs. 3-5), and good agreement is reported. The Appendix maps the |n|=1 truncation of the finite-size formula onto Lüscher's dispersive form and reports a factor-of-two discrepancy in the leading exponential term, attributed to nucleon-exchange diagrams.","tokens_in":19020,"tokens_out":8089,"duration_ms":137818,"significance":"The calculation is a genuine first-principles derivation: the finite-volume formula introduces no new parameters beyond those already fixed from the infinite-volume chiral fit, and the L-dependence is therefore independent structure rather than a fitted shape. The Appendix cross-check against Lüscher's formula is a real independent consistency test, and the paper reports the factor-of-two discrepancy honestly rather than suppressing it. The comparison spans four collaborations and several lattice actions, which is a meaningful consistency check, and the prediction is falsifiable by future simulations at smaller pion mass. These are substantive strengths.","major_comments":[{"comment":"The finite-size curves are normalized so that mN(L) on the largest lattice equals the Monte Carlo value, so only the shape of the L-dependence is tested, and that shape is never recomputed with the alternative LECs of Fit 3 or Fit 4, nor with c3 = -4.7 GeV^-1, nor with m0 and c1 varied within their quoted errors (m0 = 0.89(6) GeV, c1 = -0.93(5) GeV^-1). Because the headline claim of [abstract:s3] is 'remarkable agreement', the reader needs a quantified uncertainty band: please recompute the curves in Figs. 3-5 with the Fit 3 and Fit 4 parameters and with the LECs varied within 1 sigma, and show the resulting spread against the data errors (point 31 alone carries +/-82 MeV).","section":"Section 5, Figs. 3-5","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The paper honestly reports that its |n|=1 truncation reproduces Lüscher's formula only up to a factor of two in the leading exponential term, and attributes the difference to 'nucleon exchange diagrams' in Lüscher's treatment. That attribution is asserted rather than demonstrated. Since the main comparison in Figs. 3-5 uses the full sum over n, the discrepancy does not directly affect those figures, but a reader cannot judge from the text whether the |n|=1 truncation is trustworthy or whether the mapping onto Lüscher's dispersive form is incomplete. Please either demonstrate quantitatively that the missing piece is the nucleon-exchange contribution, or state explicitly why the full sum is unaffected and how the factor-of-two should be interpreted.","section":"Appendix, Eqs. (29) and (50)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The chi^2 values are quoted without the number of degrees of freedom. Fit 1 has chi^2 = 12.18 for ten points and three parameters (dof = 7, chi^2/dof approximately 1.74), which the text describes as describing the data 'very well'; Fit 3 has chi^2 = 0.29 for four points and two parameters (dof = 2, chi^2/dof approximately 0.15), which is small enough to suggest underestimated errors or overfitting. Since the LECs from these fits are the bridge to the finite-size prediction, please quote dof and p-values and state whether the errors in Tables 1 and 2 include all systematics relevant to the fit.","section":"Table 4","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"Points 35 (L = 1.73 fm, mN = 1.510(26) GeV) and 40 (L = 2.16 fm, mN = 1.509(18) GeV) agree within errors despite a 25% change in L; please quote the predicted shift numerically so the reader can judge how much of the agreement is nontrivial.","section":"Section 5, Fig. 5","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The statement that the agreement 'would deteriorate' if only pions with |n| = 1 were kept is a quantitative claim that is not shown; add the |n| = 1 curve to Figs. 3-5 or quote the numbers.","section":"Section 5","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"m_pi is taken from the largest lattice and its own finite-size dependence is neglected; state the expected size of this effect relative to the nucleon-mass effect being predicted.","section":"Section 4 and Table 3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The phrase 'parameter-free prediction of the finite size effects' should be qualified as in Section 5, where the absolute normalization is anchored to the largest-volume Monte Carlo point; the shape is predicted given the LECs, not the absolute mass.","section":"Abstract and Section 5","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Point 13 is excluded from the LEC fit but used in Fig. 4; state this explicitly so the reader does not think the fit and the comparison use the same set.","section":"Section 5","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Point 12 has L = 1.04 fm and m_pi L = 4.4, the smallest box in the comparison; note that this is at the lower edge of the p-expansion regime and discuss whether the agreement there is meaningful.","section":"Tables 1-2 and Fig. 4","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The scale is set with r0 = 0.5 fm, and the paper itself flags the reliability of this procedure (Ref. [15]); consider quoting how much the curves in Figs. 3-5 would move under a plausible variation of r0.","section":"Section 4","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Eq. (3) writes 'ma' while Eq. (6) writes 'Delta a', and 'er1(lambda)' and 'e1' appear in different places; please harmonize the notation.","section":"Notation","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"major_revision"},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"The thing worth knowing: this paper derives the finite-volume correction to the nucleon mass in relativistic SU(2) baryon ChPT through O(p^4) (Eqs. 15, 23, 24), and the derivation is correct as far as I can check. The sum-integral replacement is standard Hasenfratz-Leutwyler, the Bessel reductions are explicit, and the Appendix cross-check against Lüscher's dispersive form is a real independent test rather than decoration. That check is also the most interesting part: they find their leading exponential term is twice Lüscher's, and they say so plainly, attributing the difference to nucleon exchange diagrams in his treatment. I would want to see that resolved by someone, because they note the exponential and the integral nearly cancel in Lüscher's form, so a factor of two is not cosmetic. But they flag it themselves, which is the right instinct. The infinite-volume formula is quoted from Ref. [1]; the finite-volume extension is theirs and is new. The data compilation across UKQCD/QCDSF, CP-PACS and JLQCD is genuinely useful, and the observation that different actions land on one curve is a nice sanity check on lattice artefacts. On the soft spots: the curves in Figs. 3-5 are normalised so that mN(L) matches the Monte Carlo at the largest volume, so what is tested is the shape of the L-dependence, not the absolute prediction. That is a legitimate test but weaker than 'parameter-free prediction' implies, and the agreement is visual, with no error band from the LECs. The stress-test note is right that the curves are never recomputed with Fit 3/Fit 4 or with m0, c1 varied inside their errors. Fit 3 and Fit 4 do show the LECs are reasonably stable, which blunts the concern, but propagating the uncertainty into the curves is a half-day's work and would settle it. The ten-point selection is post hoc, and the r0 = 0.5 fm scale setting is a known weak point that the authors acknowledge. None of this breaks the central claim. Who should read it: anyone doing chiral extrapolations of baryon masses who needs to decide whether small-volume points can be corrected rather than discarded. It deserves referees, and I would cite the formulae.","headline":"A clean, useful finite-volume calculation whose 'parameter-free' agreement is really a test of the shape of L-dependence, anchored at the largest box and never given an error band.","tokens_in":19767,"tokens_out":861,"duration_ms":13723},"pith_extraction":{"msc":[],"pacs":["11.15.Ha","12.38.Gc"],"model":"selfhosted/deepseek-v4.1-flash","headline":"Chiral perturbation theory predicts how the nucleon mass changes when the lattice box is made smaller, with no free parameters left in the finite-volume formula.","keywords":["lattice QCD","chiral perturbation theory","nucleon mass","finite size effects","N_f=2 dynamical fermions","baryon chiral Lagrangian","Lüscher finite-volume formula","chiral extrapolation"],"falsifier":"A lattice study holding the pion mass fixed while varying the box size over a range where Eq. 24 makes a definite prediction, and finding the measured nucleon mass outside the predicted band, would break the claim; equally, refitting the infinite-volume parameters with a substantially stricter pion-mass cut and finding the resulting finite-volume curves move appreciably would show the prediction is not parameter-free in practice.","tokens_in":20066,"feed_emoji":"⚛️","tokens_out":3075,"duration_ms":51257,"temperature":0.7,"pith_summary":"Lattice simulations of QCD place the nucleon inside a finite box, and the measured mass drifts with box size because pions wrap around the periodic boundary. This paper computes that drift analytically in relativistic SU(2) baryon chiral perturbation theory through order $p^4$, obtaining a formula (Eq. 24) whose only inputs are low-energy constants already fixed from the pion-mass dependence of large-volume nucleon masses plus a few phenomenological couplings. The result is a parameter-free prediction: once the large-volume curve is anchored, the volume dependence at smaller $L$ is determined. The authors then confront this prediction with N_f=2 dynamical data from UKQCD, QCDSF, CP-PACS and JLQCD, using matched sets of simulations that differ only in lattice size, and find the predicted curvature tracks the data at pion masses of roughly 545, 717 and 732 MeV. If this holds up, finite-volume effects in nucleon-mass calculations become a computable correction rather than an extrapolation uncertainty.","feed_headline":"Chiral theory predicts how nucleon mass shrinks in small boxes","feed_subtitle":"With couplings fixed from large-volume data, the O(p^4) formula needs no new parameters to match simulations.","key_machinery":"The engine is the standard sum-minus-integral identity for a massive propagator in a periodic box (Eq. 7, following Hasenfratz and Leutwyler), which converts $L^{-3}sum_{p} - (2pi)^{-3}int d^3p$ of $(p^2+M^2)^{-r}$ into a sum over integer winding vectors $n$ of modified Bessel functions $K_{r-3/2}(L|n|M)$. Applied to the pion-nucleon loop after Feynman parametrisation and a Wick rotation, this turns the finite-volume shift into a rapidly converging Bessel series whose terms are labelled by how many times the pion circumnavigates the box. The same identity, applied to the $O(p^4)$ counterterm graph, gives the $K_1$ and $K_2$ terms weighted by $2c_1-c_3$ and $c_2$. The appendix re-expresses the $|n|=1$ pieces in the dispersive form used by Lüscher, which relates them to the pion-proton forward scattering amplitude, and finds the leading exponential term carries a coefficient twice Lüscher's.","core_discovery":"The paper's claim is that the volume dependence of the nucleon mass is under control in chiral effective field theory. Two graphs contribute at the stated order: a one-loop pion-nucleon graph giving the $O(p^3)$ piece $m_N(L)-m_N(0) = 3g_A^2 m_0 m_2^2/(16pi^2 f_pi^2) int_0^infty dx sum'_n K_0(L|n| sqrt(m_0^2 x^2 + m_pi^2(1-x)))$ (Eq. 15), and a counterterm graph giving an $O(p^4)$ piece built from $c_1, c_2, c_3$ and Bessel functions $K_1, K_2$ of $L|n|m_pi$ (Eq. 23). Neither piece introduces a new low-energy constant. Fixing $m_0$, $c_1$ and $e_1^r$ from a fit to ten large-volume, fine-lattice points with $m_pi L > 5$ and $m_pi < 800$ MeV (Fit 1: $m_0 = 0.89(6)$ GeV, $c_1 = -0.93(5)$ GeV$^{-1}$, $e_1^r = 2.8(4)$ GeV$^{-3}$), the authors overlay the resulting curves on three same-coupling, different-volume data sets and report agreement. They also state that truncating the winding sum to $|n|=1$ would spoil that agreement, so pions crossing the box more than once matter numerically.","pith_inferences":["If the chiral finite-volume series is trustworthy at these pion masses, the same machinery should predict volume shifts for other baryon masses and for matrix elements such as $g_A$; agreement or disagreement there would be an independent stress test of the constants extracted here.","The factor-of-two mismatch with Lüscher's leading term is a sharp, checkable discrepancy rather than a rounding issue, since the authors note the exponential and the integral nearly cancel in Lüscher's bookkeeping; a dedicated lattice study at very large $m_pi L$ where only the leading term survives could discriminate between the two bookkeepings.","The paper's reliance on $r_0 = 0.5$ fm to set the scale is flagged by the authors themselves as a weak point; recasting the analysis in dimensionless ratios, for instance $m_N/f_pi$ versus $m_pi/f_pi$, would remove that dependence and could be done with the same data tables."],"forward_implications":["Finite-volume corrections to the nucleon mass at a given pion mass and box size become a calculable number rather than an extrapolation ansatz, provided the low-energy constants are known.","The same constants that describe the pion-mass dependence of $m_N$ also describe its volume dependence, so lattice groups can cross-check their chiral fits against their finite-size systematics.","Because the winding sum must include $|n| > 1$ to match the data, analyses that keep only pions crossing the box once will underestimate the effect.","The comparison with Lüscher's formula isolates a factor-of-two difference in the leading term, attributed by the authors to nucleon-exchange diagrams that their chiral loop does not include."],"supporting_citations":[{"why":"Supplies the infinite-volume $O(p^4)$ nucleon-mass formula (Eq. 16) and the fitting procedure whose parameters are then fed into the finite-volume prediction.","marker":"[1]"},{"why":"Provides the relativistic baryon chiral Lagrangian and the infrared regularisation scheme used to define the loop integrals.","marker":"[9]"},{"why":"Provides the sum-minus-integral identity (Eq. 7) that converts momentum sums in a periodic box into winding-number Bessel sums.","marker":"[4]"},{"why":"Source of the Bessel and integral representations (formulas 3.471.9 and 6.596.3) used to evaluate the loop and counterterm integrals in closed form.","marker":"[11]"},{"why":"Lüscher's finite-volume formula, the benchmark against which the paper's $|n|=1$ terms are recast and compared in the appendix.","marker":"[12]"},{"why":"Source of the UKQCD nucleon masses and of the $r_0/a$ determination used to convert lattice numbers to physical units.","marker":"[17]"},{"why":"CP-PACS N_f=2 data, contributing several of the large-volume points used in the chiral fit and in the volume comparisons.","marker":"[18]"},{"why":"JLQCD N_f=2 data, supplying the same-coupling, different-volume sets at $m_pi approx 545$ and $732$ MeV used in the finite-size comparison.","marker":"[19]"},{"why":"Provides the pion-nucleon Lagrangian and the pseudovector derivative coupling that fixes the one-loop vertex in graph (a).","marker":"[10]"}],"fun_headline_variants":["The finite-volume nucleon mass needs no new low-energy constants","The O(p^4) chiral formula matches nucleon masses across box sizes","Winding pions beyond |n|=1 are needed to match the lattice masses","Fitted large-volume couplings predict the nucleon mass in smaller boxes","Chiral perturbation theory describes the nucleon mass in small boxes"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The low-energy constants fitted to infinite-volume nucleon masses at pion masses up to about 800 MeV are assumed to be the same constants that control the finite-volume formula, so any higher-order contamination of that fit passes straight into the supposedly parameter-free prediction.","fun_headline_variants_meta":{"raw":{"variants":["The finite-volume nucleon mass needs no new low-energy constants","The O(p^4) chiral formula matches nucleon masses across box sizes","Winding pions beyond |n|=1 are needed to match the lattice masses","Fitted large-volume couplings predict the nucleon mass in smaller boxes","Chiral perturbation theory describes the nucleon mass in small boxes"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":2201,"prompt_tokens":1012,"completion_tokens":1189,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":1108}},"tokens_in":1012,"tokens_out":1189,"duration_ms":20118,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":1108,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-11-13T08:57:08.128488+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"pipeline_version":"pith-pipeline@v1.5.0-archive-selfhost","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":"A lattice study holding the pion mass fixed while varying the box size over a range where Eq. 24 makes a definite prediction, and finding the measured nucleon mass outside the predicted band, would break the claim; equally, refitting the infinite-volume parameters with a substantially stricter pion-mass cut and finding the resulting finite-volume curves move appreciably would show the prediction is not parameter-free in practice.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the infinite-volume $O(p^4)$ nucleon-mass formula (Eq. 16) and the fitting procedure whose parameters are then fed into the finite-volume prediction."},{"cited_title":"Becher and H","cited_arxiv_id":null,"evidence_quote":"Provides the relativistic baryon chiral Lagrangian and the infrared regularisation scheme used to define the loop integrals."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the sum-minus-integral identity (Eq. 7) that converts momentum sums in a periodic box into winding-number Bessel sums."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the Bessel and integral representations (formulas 3.471.9 and 6.596.3) used to evaluate the loop and counterterm integrals in closed form."},{"cited_title":"L¨ uscher, in:Progress in gauge field theory , eds","cited_arxiv_id":null,"evidence_quote":"Lüscher's finite-volume formula, the benchmark against which the paper's $|n|=1$ terms are recast and compared in the appendix."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the UKQCD nucleon masses and of the $r_0/a$ determination used to convert lattice numbers to physical units."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"CP-PACS N_f=2 data, contributing several of the large-volume points used in the chiral fit and in the volume comparisons."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"JLQCD N_f=2 data, supplying the same-coupling, different-volume sets at $m_pi approx 545$ and $732$ MeV used in the finite-size comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the pion-nucleon Lagrangian and the pseudovector derivative coupling that fixes the one-loop vertex in graph (a)."}],"lane_accuracy":{"status":"not yet measured","lane":{"model":"selfhosted/deepseek-v4.1-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}