{"id":"64b695d7-5166-4dcb-a888-e129b701bafb","arxiv_id":"1908.08671","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"In the NJL model, spherical MIT boundary conditions produce stronger finite-size suppression of the quark condensate than antiperiodic boundary conditions.","lead":"The authors study how quark mass and chiral symmetry breaking change when the NJL model is placed in a sphere with MIT bag boundary conditions. They report that this spherical boundary condition suppresses the quark mass much more strongly than an antiperiodic box of equal volume.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (13) omits the 2|κ| magnetic degeneracy of spherical MIT bag modes; the mode density is therefore not the one used in the infinite-volume comparison unless the numerical code silently included it.","rationale":"The reader's weakest assumption is exactly the mode-sum replacement in Eq. (13), and that is also the most load-bearing point in my reading. The spherical MIT bag spectrum has a magnetic degeneracy 2|κ| per κ-shell; the paper does not include it explicitly, and its only justification for the 1/(2V) prefactor concerns the two signs of κ, not the m degeneracy. Because this replacement controls the density of modes that drives chiral symmetry breaking, a literal reading of Eq. (13) would invalidate the quantitative comparison with the antiperiodic box and could change the qualitative conclusion. At the same time, the paper's Fig. 2 strongly suggests the numerical calculation may have used the proper degeneracy, since the displayed integrated mode density is said to fluctuate around the infinite-volume limit. That ambiguity makes the appropriate verdict CONDITIONAL, not REJECT: the authors should state the mode-sum normalization explicitly, report mode counts, or release code. Since the reader already reached CONDITIONAL for the same reason, my pass does not change the verdict. I found no independent fatal flaw; the inhomogeneous-condensate neglect is a secondary uncontrolled assumption but not needed for the concern I raise.","tokens_in":6847,"tokens_out":10380,"duration_ms":110520,"concrete_test":"Recompute the zero-temperature gap equation using the explicit replacement (1/V)∑_{κ=±1}^{±∞}(2|κ|)∑_n e^{-τ p_{nκ}^2}, where p_{nκ} are the roots of Eq. (12), and compare M(R) at R=1,3,5,7,10,14 fm with Fig. 1. In parallel, compute the integrated mode count N(P) with and without the 2|κ| factor for R=0.985 fm and 1.97 fm and compare both to V∫_0^P d^3p/(2π)^3 = V P^3/(6π^2). If the no-degeneracy N(P) differs from the continuum curve by orders of magnitude at P≈1–2.5 GeV, Eq. (13) as written is not a valid replacement; if the degeneracy-included N(P) reproduces Fig. 2 and shifts M(R) by less than ~10%, the issue is a notation gap rather than a numerical error.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claim depends on the replacement (13). In the spherical MIT bag, each κ labels a Dirac angular-momentum shell with j=|κ|-1/2 and magnetic degeneracy 2|κ|; the mode momentum p_{nκ} is independent of m, so a mode sum must contain the factor 2|κ|. Eq. (13) as written sums over κ>0 and κ<0 with equal weight, and the accompanying comment only explains the factor 2 in the denominator as 'nondegeneracy of κ and −κ states'. If the sums over p_k are implemented literally, the integrated density of states scales like P^2R^2 instead of P^3R^3 and cannot approach the infinite-volume limit; the comparison with the antiperiodic box in Fig. 1 would be invalid. The paper's Fig. 2 claims the integrated mode density fluctuates around the infinite-volume limit, which would be true only if the 2|κ| degeneracy was included in the computation. The text never states this, reports no mode counts, and provides no code, so the reader cannot tell whether Eq. (13) is just badly notated or the numbers are wrong. The constant-condensate assumption after Eq. (9) is also uncontrolled, but the mode-sum normalization is the load-bearing step.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the two-flavor NJL model at finite temperature and zero density inside a sphere with the MIT boundary condition. The authors replace the infinite-volume momentum integral by sums over the discrete spherical MIT bag modes, solve the mean-field gap equation with proper-time regularization, and compare the resulting constituent quark mass with the antiperiodic-boundary box result. They report that the MIT boundary condition produces much stronger finite-size effects: the constituent quark mass returns to the infinite-volume value only for a sphere radius of about 14 fm, whereas an antiperiodic box is effectively infinite for L > 3 fm. The paper also presents finite-temperature masses, chiral susceptibility, and integrated mode densities to support this claim.","tokens_in":57,"tokens_out":7600,"duration_ms":198295,"significance":"If the quantitative claim is correct, the paper is valuable because it shows that the choice of spatial boundary condition can substantially change effective-model estimates of finite-size effects in heavy-ion physics, and it provides a direct mode-sum treatment of a spherical boundary rather than an asymptotic multiple-reflection expansion. The use of proper-time regularization and the standard Matsubara treatment are sound in the infinite-volume limit. However, the central numerical comparison depends on the normalization of the mode-sum replacement in Eq. (13), which is asserted rather than derived and appears to omit the 2|kappa| angular degeneracy of MIT bag modes; until this issue is resolved, the 14 fm scale and the comparison with the antiperiodic box are not established.","major_comments":[{"comment":"The replacement in Eq. (13) is the load-bearing step of the paper, but its normalization is not derived and, as written, it omits the angular degeneracy of spherical MIT bag modes. For each kappa, the Dirac angular-momentum quantum number is j = |kappa|-1/2 and there are 2|kappa| magnetic substates with the same radial momentum p_{n,kappa}; a mode sum must carry this multiplicity. Equation (13) sums over kappa>0 and kappa<0 with weight one per kappa and divides by 2V, while the accompanying sentence only explains the factor 2 by 'nondegeneracy of kappa and -kappa states', which neither supplies the missing 2|kappa| nor clarifies why the two sign families are combined with a factor 1/2. If Eq. (13) is implemented literally, the integrated density of states will scale like P^2 R^2 rather than P^3 R^3 for large P, so it cannot fluctuate around the infinite-volume limit as claimed in Fig. 2, and the comparison in Fig. 1 would be invalid. The text reports no mode counts and provides no code, so the reader cannot tell whether Eq. (13) is merely badly notated or whether the numerical results are wrong. The authors must state the multiplicity explicitly, correct the replacement if needed, and repeat the numerical analysis.","section":"§2, Eq. (13)"},{"comment":"After Eq. (9), the authors state that the condensate is inhomogeneous in a finite system in general but that they 'neglect the inhomogeneous effects' and treat <psi psi> as constant. In a small confining sphere this is an uncontrolled mean-field ansatz, and the MIT bag boundary condition in particular is known to produce strong spatial variation of scalar densities near the surface. Since the paper's central message concerns the quantitative size of finite-size effects, this approximation needs at least a quantitative estimate of its error, for example by comparing with a spatially dependent mean-field profile or with a heat-kernel/MRE estimate. Without such a check, the magnitude of the reported 14 fm scale remains uncertain even after the mode-sum normalization is corrected.","section":"§2, after Eq. (9)"}],"minor_comments":[{"comment":"The two angular-momentum labels in Eq. (12) are typeset with the same symbol, making the eigenvalue equation ambiguous; they should be distinguished with different notation, such as l_kappa and \\bar{l}_kappa.","section":"§2, Eq. (12)"},{"comment":"The statement that the constituent quark mass 'gets very close' to the infinite-volume value at R about 14 fm is not quantified; the authors should specify a tolerance (for example 1% or 5%) and plot M(R)/M(infinity) so the threshold is well defined.","section":"§3, Fig. 1"},{"comment":"The numerical implementation is not described: there is no algorithm for finding the roots of Eq. (12), no statement of how many modes were included, and no convergence checks. A short table of the first roots for one or two radii would greatly improve reproducibility.","section":"§3"},{"comment":"There are several typographical errors, including 'relavent' in Section I, 'intergral' after Eq. (9), and 'volumn' in the caption of Fig. 2; these should be corrected in a revised version.","section":"throughout"},{"comment":"No sensitivity analysis with respect to the model parameters m, G, and tau_UV is presented; since the finite-size crossover scale may depend on the proper-time cutoff, a brief parameter scan would strengthen the quantitative claim.","section":"§3"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things. First, this is the first paper I know of that feeds the exact spherical MIT bag spectrum into the NJL gap equation, rather than using the multiple reflection expansion. That is a real step forward for effective finite-size QCD studies. Second, the headline result — that the MIT sphere gives much stronger finite size effects than the antiperiodic box — rests on a mode-counting step that the paper never actually pins down.\n\nThe method itself is sensible. The proper-time regularization avoids the awkwardness of a sharp momentum cutoff in a cavity, the MIT boundary condition is the natural confining choice for fermions on a sphere, and comparing with the antiperiodic box at equal volume is a legitimate question. The figures are clear and the narrative is easy to follow. As a transparent model calculation, the paper earns credit.\n\nThe soft spot is Eq. (13). The replacement of the momentum integral by 1/(2V) times the sum over positive and negative kappa is stated with a one-line comment about the factor 2 coming from \"nondegeneracy of κ and −κ states.\" But in the spherical MIT spectrum, every κ mode carries a magnetic degeneracy 2|κ|, since j = |κ| − 1/2. The written sum has no such factor. If the numerical implementation literally followed Eq. (13), the integrated mode density would scale like P^2 R^2 instead of P^3 R^3 and could not approach the infinite-volume limit. That would invalidate Fig. 2 and Fig. 1. The authors never state whether the 2|κ| degeneracy is included in the sums, report no mode counts, and give no code. They also do not show any check that the discrete mode density reproduces the infinite-volume continuum density, aside from the qualitative \"fluctuates around\" claim in Fig. 2.\n\nIs the paper wrong? I cannot tell. The most plausible reading is that the written formula is sloppy notation and the code silently included the degeneracy — but the text gives the reader no way to confirm this. The constant-condensate assumption after Eq. (9) is also a real approximation in a finite confining sphere, though it is standard in the NJL finite-size literature. Borrowing the parameters from earlier box calculations is fine, but the authors do not discuss whether the same coupling and cutoff remain appropriate in a sphere.\n\nMy read: the qualitative message — confining boundary conditions can enhance finite-size effects relative to periodic/antiperiodic boxes — may survive a corrected calculation. But the specific 14 fm versus 3 fm comparison is not established by this manuscript as written.\n\nRecommendation: send it to peer review. A referee can ask for the derivation of Eq. (13), a clear statement of how degeneracies enter, a mode-count check against the infinite-volume density, and ideally a reproducibility note. This is not a desk reject; it is a potentially useful model study that needs one solid revision.","headline":"First MIT-bag mode-sum NJL calculation in a sphere, but the load-bearing mode counting in Eq. (13) is asserted rather than derived, so the quantitative comparison to the antiperiodic box is unverified.","tokens_in":7623,"tokens_out":2728,"would_cite":false,"duration_ms":30739,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Lg","12.38.Mh","64.60.an"],"model":"deepseek-v4-flash","headline":"Spherical confinement, not just volume, controls how far finite-size effects reach in quark matter.","keywords":["Nambu-Jona-Lasinio model","chiral phase transition","finite size effects","MIT boundary condition","spherical cavity","constituent quark mass","heavy-ion collisions","proper time regularization"],"falsifier":"Recalculate the gap equation with the full mode multiplicity of the spherical MIT spectrum, counting each $\\kappa$ level $2|\\kappa|$ times, and compare M(R) with the paper's threshold; if the corrected density of states makes M saturate well before R = 14 fm, the paper's quantitative comparison with the antiperiodic box fails, while if the threshold persists, the claim is supported.","tokens_in":6611,"feed_emoji":"⚛️","tokens_out":8688,"duration_ms":83374,"temperature":0.7,"pith_summary":"This paper computes the two-flavor Nambu-Jona-Lasinio model inside a hard sphere with the MIT bag boundary condition and compares it with the usual antiperiodic box. It tries to establish that the spherical confining boundary produces much stronger finite-size effects: the constituent quark mass only nears its infinite-volume value when the sphere radius reaches about 14 fm, while an antiperiodic box reaches that limit for L > 3 fm. The authors care because heavy-ion collision fireballs are closer to spheres than to boxes, so realistic finite-size corrections to chiral symmetry restoration may be larger than previous box-based estimates.","feed_headline":"Quark spheres show finite-size effects to 14 fm","feed_subtitle":"MIT boundary condition in a spherical NJL model keeps quark masses small far beyond the 3 fm antiperiodic-box threshold","key_machinery":"The central object is the mode-sum replacement of Eq. (13), in which the continuum momentum integral in the gap equation is replaced by a sum over the discrete momentum modes allowed by the spherical MIT boundary condition, weighted by 1/(2V) over both signs of the Dirac quantum number $\\kappa$. Those allowed momenta are the solutions of the spherical-cavity eigen-equation for the free Dirac equation, involving spherical Bessel functions and the sign of $\\kappa$. The sum is inserted into the proper-time-regularized mean-field gap equation, so finite size enters through missing and shifted modes rather than through a momentum cutoff. This discrete spectrum is what makes M(R) recover so slowly with radius.","core_discovery":"On its own terms, the paper's central discovery is that the spatial boundary condition, not just the volume, controls how strongly finite size acts on dynamical chiral symmetry breaking. With the MIT boundary condition on a sphere, the zero-temperature constituent quark mass M(R) remains well below the infinite-volume value until R ≈ 14 fm; with antiperiodic boundary conditions in a box, the mass already approaches the infinite-volume value for L > 3 fm. At finite temperature the same pattern appears: at equal volume, the spherical-MIT constituent mass is smaller, and the chiral susceptibility peak is smoothed. Since quark-gluon plasma droplets in heavy-ion collisions are estimated to be 2–10 fm in size, the paper concludes that those droplets experience considerable finite-size effects, stronger than earlier box-based estimates.","pith_inferences":["Editorial inference: if the correct mode multiplicity ($2|\\kappa|$ per level) is included, the 14 fm threshold could shift substantially; this is the paper's most testable quantitative assumption.","Editorial inference: the assumption of a uniform condensate inside a small confining sphere is uncontrolled; a spatially varying condensate with surface enhancement or depletion could change the effective finite-size scale.","Editorial inference: if stronger spherical finite-size effects survive a corrected mode count, they would also affect other observables, such as pion properties and the location of the chiral critical endpoint, not just the constituent mass.","Editorial inference: combining the spherical MIT mode sum with a gap equation that allows inhomogeneous condensates would give a sharper test of whether the 2–10 fm fireball indeed sits in the finite-size-dominated regime."],"forward_implications":["At fireball sizes of 2–10 fm, the chiral condensate is substantially suppressed compared with the infinite-volume limit, so finite-size effects should be included in heavy-ion phenomenology.","The antiperiodic-box approximation underestimates the finite-size correction; other boundary-condition choices deserve systematic comparison in effective QCD models.","The same brute-force mode-sum method can be applied to other confining shapes to separate shape effects from volume effects.","Chiral-susceptibility peaks are smoothed by the spherical boundary, meaning finite-size signatures near the phase transition are visible at larger volumes than box calculations suggest."],"supporting_citations":[{"why":"Supplies the 2–10 fm size estimate for quark-gluon plasma droplets that motivates the finite-size study.","marker":"[1]"},{"why":"Provides the two-flavor NJL Lagrangian, the mean-field gap equation, and the proper-time regularization used in the calculation.","marker":"[35]"},{"why":"Gives the box-mode-sum method and regularization scheme that the paper extends to the spherical MIT case.","marker":"[36]"},{"why":"Introduces the MIT bag boundary condition for fermions in a spherical cavity.","marker":"[39]"},{"why":"Completes the MIT bag model boundary-condition setup for the sphere.","marker":"[40]"},{"why":"Supplies the spherical-cavity eigen-equation that determines the allowed discrete quark momenta under the MIT boundary condition.","marker":"[42]"}],"fun_headline_variants":["MIT sphere keeps quark mass low up to 14 fm","Boundary condition, not volume, delays quark mass growth","Sphere with MIT condition shrinks chiral symmetry breaking","NJL sphere: MIT boundary extends finite-size effects to 14 fm"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the mode sum in Eq. (13) correctly counts the spherical quark states—each $\\kappa$ level is summed once, without separately including its $2|\\kappa|$-fold angular degeneracy—and that the chiral condensate stays uniform inside the sphere; if either assumption fails, the 14 fm threshold changes.","fun_headline_variants_meta":{"raw":{"variants":["MIT sphere keeps quark mass low up to 14 fm","Boundary condition, not volume, delays quark mass growth","Sphere with MIT condition shrinks chiral symmetry breaking","NJL sphere: MIT boundary extends finite-size effects to 14 fm"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00131,"raw_usage":{"total_tokens":5243,"prompt_tokens":749,"completion_tokens":4494,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":365,"completion_tokens_details":{"reasoning_tokens":4426}},"tokens_in":365,"tokens_out":4494,"duration_ms":30468,"temperature":1.0,"reasoning_tokens":4426,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:37:01.104685+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recalculate the gap equation with the full mode multiplicity of the spherical MIT spectrum, counting each $\\kappa$ level $2|\\kappa|$ times, and compare M(R) with the paper's threshold; if the corrected density of states makes M saturate well before R = 14 fm, the paper's quantitative comparison with the antiperiodic box fails, while if the threshold persists, the claim is supported.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the 2–10 fm size estimate for quark-gluon plasma droplets that motivates the finite-size study."},{"cited_title":"Wang, Y .-H","cited_arxiv_id":null,"evidence_quote":"Gives the box-mode-sum method and regularization scheme that the paper extends to the spherical MIT case."},{"cited_title":"Chodos, R","cited_arxiv_id":null,"evidence_quote":"Introduces the MIT bag boundary condition for fermions in a spherical cavity."},{"cited_title":"Greiner, S","cited_arxiv_id":null,"evidence_quote":"Supplies the spherical-cavity eigen-equation that determines the allowed discrete quark momenta under the MIT boundary condition."}],"review_version":1}