Pentaquark Bound States and Regge Trajectories in QCD via Bethe Salpeter Formalism

Scientific Impact Assessment

Core Contribution

The paper applies the Godfrey–Isgur (GI) relativized quark-model potential — adapted to an effective meson–baryon channel — within an instantaneous-approximation Bethe–Salpeter (Salpeter-reduced) framework to compute masses of four established hidden-charm pentaquarks (Pc(4440), Pc(4457), Pcs(4338), Pcs(4459)) and their first two radial excitations. It then fits radial Regge trajectories in the (n, M²) plane. This is an incremental extension of the authors' prior program applied successively to tetraquarks, hexaquarks, and now pentaquarks with essentially the same toolkit.

Methodological Rigor

This is where the paper has serious problems. Several issues undermine confidence:

1. Internal inconsistency between Table 3 and Table 5. The computed radial excitations (e.g., Pc(4440): 4442.5, 4444.7, 4446.9 MeV) are separated by only ~2.2 MeV. These imply M² increments of ~0.02 GeV² per radial node. Yet Table 5 reports Regge slopes of ~1.15 GeV². A slope of 1.15 GeV² would place the n=1 state near 4.57 GeV, not 4.4447 GeV. The Regge fit parameters are therefore numerically inconsistent with the very masses they are supposedly fit to. This is a glaring, disqualifying contradiction.

2. Unphysical radial spacing. Genuine radial excitations of hadronic systems are typically separated by hundreds of MeV. A ~2 MeV spacing between "ground state" and "second radial excitation" is physically implausible and suggests the numerical eigenvalue extraction is not producing true distinct radial modes, or that the binding dynamics are essentially trivial.

3. Trivial "agreement." The pentaquark masses sit ~15–20 MeV below the constituent meson–baryon threshold (masses taken from PDG). With constituent masses fixed to experiment and only shallow binding, "excellent agreement within 3.8 MeV" is nearly guaranteed and carries little discriminating power. The paper itself acknowledges it cannot distinguish molecular from compact interpretations.

4. Parameter ambiguity. The text repeatedly stresses the GI parameters are "not refitted," yet introduces multiple effective channel-dependent coefficients (C_MB, C_ss^MB, effective mass m̃) whose values and determination are not transparently specified. The claim of "genuine predictive power" is thus overstated.

5. Broken manuscript elements. Figure cross-references appear as "Figs.??–??", and the abstract/body oscillate between claiming support for "universality of the confinement mechanism" and explicitly disclaiming it ("do not establish universality") — indicating an unresolved, internally contradictory narrative.

The numerical validation (convergence tests, cross-check against a prior hexaquark calculation) is a positive element, but it validates the integrator, not the physical soundness of the model setup.

Potential Impact

Low. The paper does not introduce a new method, does not resolve the central open question in pentaquark physics (compact vs. molecular structure — which it explicitly concedes it cannot address), and produces predictions (radial excitations ~4.445–4.462 GeV) that are physically dubious given the ~2 MeV spacings. Experimental groups are unlikely to use near-threshold, near-degenerate predictions as search benchmarks. The Regge-trajectory analysis, built on three near-collinear points with fit parameters inconsistent with the underlying masses, adds essentially no scientific content.

Timeliness & Relevance

Hidden-charm pentaquark spectroscopy is a genuinely active topic following the LHCb discoveries, and Regge/confinement universality across multiquark systems is of interest. So the *topic* is timely. However, the crowded field already contains many molecular, diquark, and QCD-sum-rule calculations (several cited) that reproduce these masses at least as well, so the marginal contribution is minimal.

Strengths & Limitations

Strengths: Computationally lightweight and reproducible in principle (Python/NumPy/SciPy, Numerov + shooting, full parameter table, runs in under a minute); a numerical cross-check against prior work is included; the writing is mostly readable and the framework is standard and transparent in structure.

Limitations: The fatal internal inconsistency between calculated masses and Regge slopes; physically implausible ~2 MeV radial splittings; near-trivial "agreement" driven by fixing constituent masses to PDG; opaque effective-parameter determination despite "no-refit" claims; contradictory conclusions about universality; broken figure references; and a future-dated arXiv stamp (2026) with a self-citation to a 2026 paper, raising provenance concerns. The work is largely a template re-application of the group's earlier tetraquark/hexaquark pipeline.

Additional Observations

The paper reads as a mechanical extension of an established in-house workflow rather than a targeted advance. The most damaging feature for any downstream use is that the two central quantitative outputs (mass spectrum and Regge slopes) are mutually contradictory, meaning a careful reader cannot trust either. Reproducibility of the *code* is plausible, but reproducibility of consistent *physics* is not demonstrated.