{"id":"4fee4956-109d-4e5f-a587-f0057931fad6","arxiv_id":"1909.00013","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Empirical helicity amplitude fits for nine nucleon and Delta resonances are modified near the pseudothreshold using a polynomial in the photon momentum matched at QP^2 = 0.1, 0.3, and 0.5 GeV^2.","lead":"Standard fits for how electrons excite protons into heavier states ignore special mathematical rules that apply when the photon momentum is zero. This paper adds those rules to the fits and shows where they change the predicted curves at low momentum transfer, especially for the Delta(1232) resonance.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Untested truncation of Eq. (4.2): the Delta(1232) conclusion depends on the 4-term polynomial ansatz whose validity over the whole low-Q^2 interval is not established.","rationale":"The reader's weakest assumption identifies the same point, and I agree it is load-bearing. The paper's own hedging for most states (Sec. VI) and the unquantified visual comparison for Delta(1232) support a conditional verdict. The ansatz in Eq. (4.2) is a plausible smooth interpolation, but the conclusion that the pseudothreshold constraints specifically force the observed low-Q^2 behavior is only as strong as the basis choice. The stated validity range of the input parametrizations (0.5–5 GeV^2) makes matching at Q_P^2 = 0.1 and 0.3 especially fragile. The proposed test — changing the interpolation basis while holding data fixed — directly targets whether the \"conclusive\" Delta(1232) statement survives. If it survives, the paper's contribution is solid; if not, the conclusion should be softened to a model-dependent demonstration, which is still worth reporting. No ad hominem concerns; the issue is an under-tested assumption.","tokens_in":25426,"tokens_out":9119,"duration_ms":75314,"concrete_test":"Refit the Delta(1232) extension with (i) a 5-term polynomial adding alpha4 qtilde^8 and (ii) a [2/2] Pade approximant in qtilde^2, using the same matching conditions and the same data as in Sec. V.C for Q_P^2 = 0.1, 0.3, 0.5 GeV^2. If Q_P^2 = 0.3 GeV^2 is no longer uniquely preferred by the low-Q^2 data (e.g., the chi-square difference between Q_P^2 values changes sign or becomes negligible), the claim that pseudothreshold constraints cannot be ignored below 0.3 GeV^2 is an artifact of the truncation in Eq. (4.2).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim — that pseudothreshold constraints \"cannot be ignored below 0.3 GeV^2\" for Delta(1232) — rests on the assumed form Eq. (4.2): A = qtilde^n (alpha0 + alpha1 qtilde^2 + alpha2 qtilde^4 + alpha3 qtilde^6). The exclusion of odd powers is motivated only by the behavior as qtilde→0 (footnote 2), and the truncation at alpha3 is not tested. No proof or numerical convergence check shows that this 4-term polynomial in qtilde^2 can represent the true amplitude over the entire interval from the pseudothreshold up to Q_P^2 = 0.5 GeV^2. Because the coefficients are fixed by matching A and its derivatives at Q_P^2 to the JLab parametrizations (which Ref. [33] states are valid for Q^2 = 0.5–5 GeV^2, while Q_P^2 = 0.1 and 0.3 are outside that range), the extension is a Hermite interpolation with an arbitrary basis. The visual preference for Q_P^2 = 0.3 over 0.1 and 0.5 in Fig. 3 is therefore specific to this polynomial ansatz; a different functional basis could change which Q_P^2 best fits the same data, and with it the conclusion about where pseudothreshold constraints matter.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper addresses the low-Q^2 behavior of empirical parametrizations of the helicity amplitudes for γ*N → N* transitions. It argues that pseudothreshold constraints, which dictate specific powers of the photon three-momentum |q| and correlations among helicity amplitudes, are commonly ignored in data parametrizations and can lead to erroneous low-Q^2 descriptions. The paper proposes a method to modify an existing analytic parametrization below a chosen transition point Q_P^2 by matching the amplitude and its first (and sometimes second/third) derivatives at Q_P^2 to a polynomial expansion in |q|^2 that respects the pseudothreshold behavior. The method is applied to nine resonances using the Jefferson Lab parametrizations of Ref. [33], with Q_P^2 = 0.1, 0.3, and 0.5 GeV^2, and the resulting extensions are compared visually with data. The main quantitative claim is that for the Δ(1232) the pseudothreshold constraints cannot be ignored below 0.3 GeV^2, while for several other resonances the data are insufficient to discriminate among the choices of Q_P^2. The paper includes explicit formulas for the matching coefficients, derivative conversion tables, and full coefficient tables in appendices.","tokens_in":25735,"tokens_out":6038,"duration_ms":49578,"significance":"If the proposed method is sound, it provides a practical and general prescription for incorporating current-algebra constraints into empirical parametrizations of helicity amplitudes, which would be useful for the interpretation of low-Q^2 electroproduction data and for model comparisons. The paper is candid about the inconclusive cases, and the algebraic derivation of the matching conditions is internally consistent; the explicit coefficient tables and derivative relations (Appendices A–C) are a useful resource. The manuscript’s strengths include the generality of the formalism, the explicit treatment of all J^P = 1/2^±, 3/2^± cases, and the transparent discussion of limitations for resonances with scarce data. However, the headline Δ(1232) conclusion rests on an untested truncation of the low-Q^2 ansatz, on matching points that may lie outside the stated validity range of the input parametrization, and on a visual (not quantitative) comparison with data, in a region where the author’s own prior analysis supplies several of the data points. The result is therefore promising but not yet established to the level claimed.","major_comments":[{"comment":"The central claim in Sec. VI that for the Δ(1232) the pseudothreshold constraints ‘cannot be ignored below 0.3 GeV²’ rests on the ansatz A = ũ^n(α_0 + α_1 ũ² + α_2 ũ⁴ + α_3 ũ⁶) being valid over the whole interval from the pseudothreshold to Q_P². Footnote 2 only shows that odd powers vanish as |q|→0; it does not establish that the bracket is a quartic polynomial in ũ² throughout the interval, and no convergence test or comparison with alternative bases (including odd powers or higher-order terms) is given. Because the coefficients are fixed by matching at Q_P², the curves in Fig. 3 are not fits to the low-Q² data, so the visual preference for Q_P² = 0.3 is specific to this basis and cannot, as it stands, support the ‘conclusive’ wording in Sec. VI.","section":"Eq. (4.2) and Sec. VI"},{"comment":"The method assumes (Sec. IV) that the original parametrization ‘describe[s] well the data above Q_P²’, but the Jefferson Lab parametrizations are stated in Sec. I to be valid for Q² = 0.5–5 GeV², so Q_P² = 0.1 and 0.3 GeV² lie outside the region where the input is established; Sec. VI later says the same parametrizations ‘cover the region Q² = 0–5 GeV²’, which is internally inconsistent. The derivatives at Q_P² used to fix the extension coefficients are thus not reliable for the two smaller Q_P² values, and the conclusions drawn from comparing Q_P² = 0.1, 0.3, and 0.5 are correspondingly weakened.","section":"Sec. IV and Sec. V vs. Sec. I/VI"},{"comment":"In Sec. V.C the low-Q² Δ(1232) data (Q² < 0.15 GeV²) are not the original measurements from MAMI and MIT-Bates but are replaced by results from Refs. [41,42] (the author’s own work), converted to helicity amplitudes using MAID 2007. Since the central conclusion concerns which Q_P² best reproduces the data in this region, the use of author-derived pseudo-data introduces a circularity that the text acknowledges but does not quantify; at minimum a comparison with the original MAMI and MIT-Bates data should be shown to demonstrate that the Q_P² preference is not an artifact of the replacement.","section":"Sec. V.C (Delta(1232))"},{"comment":"The comparison in Sec. V and Figs. 1–4 is purely visual; no χ² or other quantitative measure is given. Statements such as ‘the parametrization characterized by Q_P² = 0.3 GeV² is the one that better describes the data’ (Sec. V.C) and the corresponding conclusion in Sec. VI require a quantitative comparison that accounts for the data uncertainties, especially because the data are sparse and the curves differ mainly below Q² = 0.3 GeV².","section":"Figs. 1–4 and Sec. V.C"}],"minor_comments":[{"comment":"In the 1/2⁻ row, the expression for S_{1/2} lists c_3 ũ⁵ twice; it should presumably read c_3 ũ⁷ to follow the stated pattern, and the coefficient c_0 should be bold in the table as indicated by the text.","section":"Table II"},{"comment":"The phrase ‘The states 3/2⁻ are the exception to this role’ should read ‘exception to this rule’.","section":"Sec. IV.D"},{"comment":"The caption refers to ‘tick lines’; this should be ‘thick lines’.","section":"Fig. 2 caption"},{"comment":"The abstract states ‘the invariant four-momentum square became q²’; the verb should be ‘becomes’ or ‘is’.","section":"Abstract"},{"comment":"Reference [33] is a URL; it should be replaced by a published source or by a more complete description of the parametrizations, since several conclusions depend on its validity range and functional forms.","section":"Reference [33]"}],"recommendation":"major_revision","confidential_remarks":"The paper is candid and the formalism is potentially useful, but the headline Δ(1232) claim is more model-dependent than the word ‘conclusively’ suggests. The untested truncation in Eq. (4.2), the use of Q_P² values below the stated validity range of the input parametrization, and the reliance on the author’s own data conversion in the key case all weaken the central claim. A revision that adds a basis-independence check, a quantitative comparison with data, and a test against the original low-Q² data—or that tempers the conclusions accordingly—would bring the paper to a publishable level. I would ask the editor to ensure the revised version addresses these points before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the gluing algorithm: take any smooth empirical parametrization of N* helicity amplitudes valid at moderate Q^2, and replace it below a chosen matching point QP^2 by a short power series in |q|^2 that builds in the pseudothreshold factors and amplitude correlations, matching value and derivatives at QP^2. That is a practical tool, and the paper gives enough detail to re-implement it: derivative conversion formulas in the appendices, and coefficient tables for nine states at three QP^2 values. Credit where it is due: the algebra is transparent, the appendices are useful, and the author is candid that most of the nine states are inconclusive.\n\nThe soft spots are real but not fatal. First, Eq. (4.2) truncates the |q|^2 series at |q|^6 with no sensitivity test. The even-power-only structure follows from regularity in Q^2, but the number of terms is set by how many derivatives you want to match, not by any physical scale. A different polynomial order or a different variable would give different interpolants, and the visual ranking of QP^2 = 0.3 for the Delta could shift. So \"conclusively\" is too strong. Second, the Delta comparison replaces older MAMI/MIT-Bates points below 0.15 GeV^2 with JLab/Hall A data, partly on the strength of the author's own reanalysis. That may be right, but it deserves a sensitivity check rather than a footnote. Third, the JLab input parametrizations are declared valid for Q^2 = 0.5–5 GeV^2, while two of the three matching points (0.1 and 0.3 GeV^2) are below that range, so the derivatives at those points are extrapolations of the fit, not data. Finally, agreement with data is judged by eye, with no chi-square or error propagation.\n\nNone of this kills the method. As an interpolation tool that enforces kinematic constraints, it works, and the coefficient tables are useful. The paper should be read as a methods proposal with illustrative applications, not as a determination that pseudothreshold constraints cut in at 0.3 GeV^2 for the Delta. The author should supply the input derivatives, add a higher-order sensitivity test, and quantify the Delta comparison.\n\nThis deserves peer review. Send it to a serious referee and ask for those revisions; with them, it would be a solid reference for anyone doing low-Q^2 N* fits.","headline":"A useful methods paper that enforces pseudothreshold constraints on empirical N* helicity amplitudes via a matching algorithm, but the Delta(1232) headline claim is stronger than the evidence supports.","tokens_in":26220,"tokens_out":3468,"would_cite":false,"duration_ms":32534,"reading_group":"maybe","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 empirical fits to nucleon-resonance helicity amplitudes that ignore the pseudothreshold kinematics are inconsistent at low $Q^2$, and presents an analytic matching method that enforces the constraints while…","keywords":["helicity amplitudes","pseudothreshold","nucleon resonances","electromagnetic transition form factors","analytic continuation","Siegert theorem","Delta(1232)","low momentum transfer"],"falsifier":"Precise $\\Delta(1232)$ helicity amplitudes at $Q^2=0.05$, $0.1$, and $0.2$ GeV$^2$ from electron-scattering data of the kind used in the paper would settle the issue: the $Q_P^2=0.3$ GeV$^2$ constrained extension and the original parametrization differ visibly in that interval, and data that follow the unconstrained curve would refute the claim that ignoring the pseudothreshold constraints is erroneous below $0.3$ GeV$^2$.","tokens_in":25207,"feed_emoji":"⚛️","tokens_out":10160,"duration_ms":79743,"temperature":0.7,"pith_summary":"At the pseudothreshold, where the photon three-momentum vanishes, the $\\gamma^\\ast N \\to N^\\ast$ helicity amplitudes are forced by current conservation and kinematics to have specific powers of $|\\mathbf q|$ and to satisfy correlations between transverse and longitudinal amplitudes. The paper argues that common empirical parametrizations ignore these constraints and therefore give wrong low-$Q^2$ behavior, especially for resonances whose mass is close to the nucleon mass. It proposes a method that replaces a given parametrization below a matching point $Q_P^2$ with an expansion in powers of $|\\mathbf q|$ that respects the pseudothreshold constraints and joins smoothly, with amplitude plus first, second, and sometimes third derivatives continuous. Applying this to nine resonances, the paper finds that the constraints have a definite impact for the $\\Delta(1232)$ below $0.3$ GeV$^2$ and for the $N(1520)$ below about $0.1$ GeV$^2$, while other cases remain undecided by current data. If correct, the method supplies a consistency requirement that any future low-$Q^2$ fit or theory comparison should meet.","feed_headline":"Delta data prove pseudothreshold limits cannot be ignored","feed_subtitle":"Fits that ignore the photon-momentum power laws at the pseudothreshold misdescribe the Delta below 0.3 GeV^2.","key_machinery":"The central object is the pseudothreshold expansion $A=\\tilde q^n(\\alpha_0+\\alpha_1\\tilde q^2+\\alpha_2\\tilde q^4+\\alpha_3\\tilde q^6)$, where $\\tilde q=|\\mathbf q|/M_R$ is the photon three-momentum in the resonance rest frame normalized by the resonance mass and $n=0,1,2$ encodes the required $|\\mathbf q|$ power behavior near the pseudothreshold. The coefficients are fixed by matching the amplitude, its first derivative, and its second (and sometimes third) derivative at the matching point $Q_P^2$, while the leading coefficient is set by pseudothreshold correlations such as the long-wavelength (Siegert) relation between $A_{1/2}$ and $S_{1/2}$. This machinery turns the kinematic constraints into a smooth boundary condition that any data parametrization can be forced to satisfy without refitting the higher-$Q^2$ data.","core_discovery":"The central claim is that the pseudothreshold point $Q^2=-(M_R-M_N)^2$ is not a remote technicality: for resonances close to the nucleon it sits near $Q^2=0$, and the fixed $|\\mathbf q|$ power laws, such as $A_{1/2}\\propto |\\mathbf q|$ or constant and $S_{1/2}\\propto |\\mathbf q|^2$ or $|\\mathbf q|$, together with correlations like $S_{1/2}\\propto A_{1/2}|\\mathbf q|$, control the amplitudes throughout the low-$Q^2$ region. The paper demonstrates this by constructing, for each of nine resonances, an analytic extension $A=\\tilde q^n(\\alpha_0+\\alpha_1\\tilde q^2+\\alpha_2\\tilde q^4+\\alpha_3\\tilde q^6)$ on the interval from the pseudothreshold up to a matching point $Q_P^2$, with coefficients fixed by continuity of the amplitude and its derivatives at $Q_P^2$ and by the pseudothreshold correlations. Scanning $Q_P^2=0.1$, $0.3$, and $0.5$ GeV$^2$, the paper shows that the $\\Delta(1232)$ data select the $Q_P^2=0.3$ GeV$^2$ extension and that ignoring the constraints below $0.3$ GeV$^2$ gives erroneous amplitudes, while the $N(1520)$ case selects $Q_P^2=0.1$ GeV$^2$ through the $A_{3/2}(0)$ constraint. For the remaining resonances, the available data cannot distinguish the different extensions.","pith_inferences":["The method turns the pseudothreshold constraint into a model-independent boundary condition, so the same matching could be applied to any theoretical transition current to test whether theory curves respect the same low-$Q^2$ behavior.","The $\\Delta(1232)$ case suggests a useful rule of thumb: the constraints matter up to roughly a few times $(M_R-M_N)^2$, so resonances with smaller mass splittings should show the largest pseudothreshold effects near the photon point.","The absence of odd powers of $|\\mathbf q|$ in the ansatz is a testable assumption; extending the expansion with an additional odd term and refitting would quantify the truncation error and reveal whether the data tolerate such a term.","Accurate longitudinal $S_{1/2}$ data below $Q^2=0.3$ GeV$^2$ would be the sharpest discriminator, because the paper shows that $S_{1/2}$ extensions differ most strongly across the chosen matching-point values."],"forward_implications":["Any empirical fit aiming to describe nucleon-resonance amplitudes below about $0.3$ GeV$^2$ should build in the $|\\mathbf q|$ power laws and amplitude correlations from the pseudothreshold, or it will misrepresent the region between $Q^2=0$ and the first data.","The constrained $Q_P^2=0.3$ GeV$^2$ extension for the $\\Delta(1232)$ provides a testable prediction for the amplitudes in the low-$Q^2$ gap.","For the $N(1520)$, only extensions matched at $Q_P^2=0.1$ GeV$^2$ reproduce the measured $A_{3/2}(0)$, so fits matched at higher momentum transfer should be avoided for that resonance.","For the other resonances, the method identifies which amplitudes are most sensitive to the matching point and therefore which new measurements below $Q^2=0.3$ GeV$^2$ would be most informative.","The matching procedure applies to any parametrization whose amplitudes and first derivatives are continuous, so other empirical fits can be made pseudothreshold-consistent without changing their high-$Q^2$ behavior."],"supporting_citations":[{"why":"These supply the gauge-invariant current decomposition and the pseudothreshold power-law constraints that the whole method enforces.","marker":"[23–25]"},{"why":"This provides the empirical parametrizations that the paper modifies and extends toward the pseudothreshold.","marker":"[33]"},{"why":"This supplies the compiled resonance electrocoupling data used to check that the extensions stay consistent with measurements.","marker":"[34]"},{"why":"This is used for the magnetic form factor when converting the low-$Q^2$ $\\Delta(1232)$ data into helicity amplitudes.","marker":"[4]"},{"why":"This derives the pseudothreshold relation and the form-factor representation used for the negative-parity $1/2^-$ states.","marker":"[27]"},{"why":"These derive the electric–Coulomb form-factor relations and the $3/2^\\pm$ amplitude correlations used for the $\\Delta(1232)$ and the $N(1520)$.","marker":"[24, 25]"},{"why":"This provides the recent low-$Q^2$ $\\Delta(1232)$ data that replace older points below $Q^2=0.15$ GeV$^2$ in the comparison.","marker":"[49]"}],"fun_headline_variants":["Pseudothreshold constraints now built into low-Q2 N* fits","N* parametrizations honor pseudothreshold power laws","Delta fits corrected by pseudothreshold constraints","Low-Q2 nucleon resonance fits obey pseudothreshold"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that the simple four-term power series in the photon momentum describes the amplitudes exactly across the whole interval from the pseudothreshold to $Q_P^2$, not just right at the pseudothreshold; if extra terms are needed there, the matched curves and the conclusions about which matching point fits the data would change.","fun_headline_variants_meta":{"raw":{"variants":["Pseudothreshold constraints now built into low-Q2 N* fits","N* parametrizations honor pseudothreshold power laws","Delta fits corrected by pseudothreshold constraints","Low-Q2 nucleon resonance fits obey pseudothreshold"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000326,"raw_usage":{"total_tokens":1961,"prompt_tokens":1216,"completion_tokens":745,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":832,"completion_tokens_details":{"reasoning_tokens":679}},"tokens_in":832,"tokens_out":745,"duration_ms":6932,"temperature":1.0,"reasoning_tokens":679,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:05:15.527086+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Precise $\\Delta(1232)$ helicity amplitudes at $Q^2=0.05$, $0.1$, and $0.2$ GeV$^2$ from electron-scattering data of the kind used in the paper would settle the issue: the $Q_P^2=0.3$ GeV$^2$ constrained extension and the original parametrization differ visibly in that interval, and data that follow the unconstrained curve would refute the claim that ignoring the pseudothreshold constraints is erroneous below $0.3$ GeV$^2$.","supporting_citations":[{"cited_title":"Nucleon resonance contributions to unpolarised inclusive electron scattering","cited_arxiv_id":"1904.08016","evidence_quote":"This provides the empirical parametrizations that the paper modifies and extends toward the pseudothreshold."},{"cited_title":"Drechsel and L","cited_arxiv_id":null,"evidence_quote":"This supplies the compiled resonance electrocoupling data used to check that the extensions stay consistent with measurements."},{"cited_title":"The coefficient α3 can then be de- termined by Eq","cited_arxiv_id":null,"evidence_quote":"This is used for the magnetic form factor when converting the low-$Q^2$ $\\Delta(1232)$ data into helicity amplitudes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This derives the pseudothreshold relation and the form-factor representation used for the negative-parity $1/2^-$ states."}],"review_version":1}