{"id":"b21d45a8-9814-404a-aac8-0403394924a4","arxiv_id":"2506.02844","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The rho meson decay constants f_parallel = 215 MeV, f_perp = 173 MeV, f_A = 122 MeV, and f_S = 61 MeV, together with twist-2, twist-3, and twist-4 distribution amplitudes, are extracted consistently from vector, tensor, axial-vector, and scalar currents in a self-consistent light-front quark model.","lead":"Using a quark model of hadrons built on the light front, this paper extracts the rho meson's two decay constants and its internal momentum distributions, including subleading corrections, from four different quark currents, and shows that all the extractions agree once the model's mass prescription is applied.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing assumption is the uniform M→M0 replacement inside the integrals, including α=2m/M0 in Eqs.","rationale":"Read in good faith, the paper's objective is to show that a self-consistent LFQM with on-shell constituents can extract rho decay constants and higher-twist DAs independently of the computational choices. The strongest evidence is algebraic: Table I exhibits identical operator forms across current components after the BT substitution, f∥ reproduces the PDG value, and the chiral-limit shapes match the known asymptotic forms. These give real support. However, the entire extraction program hinges on the rule that the physical mass M in Lorentz factors is replaced by M0 inside the integrals. The paper directly says this rule is needed for consistency and that the fixed physical α fails. That makes the central claim conditional on the correctness of the Type II link. The reader identified this same weakest assumption. My concern is not that the model disagrees with QCD sum rules (that would be a model-dependence issue) but that the key step is stipulated rather than derived. The recommended test is analytic: derive the M0 replacement from the covariant BS regularization, not from the requirement that results agree. Until that is done, the paper should be read as a consistent model calculation with the M0 prescription as an explicit assumption, which is exactly the CONDITIONAL verdict. I therefore keep the reader's verdict unchanged.","tokens_in":23250,"tokens_out":7649,"duration_ms":84769,"concrete_test":"Independently re-derive Eqs. (30) and (34) from the covariant BS model defining the Type II link (Ref. [30], Eq. (49); Refs. [54,55]), keeping the physical mass M in P_A^μ and P_S and applying the same Pauli–Villars regularization. If the constituent on-shell limit reproduces the uniform M→M0 replacement, including α=2m/M0 inside the integrand, for all current components and in a moving frame with P⊥≠0, the central claim is supported; if the limit leaves α at the physical mass or leaves frame-dependent residues, the M0 prescription is an unproven ansatz and the verdict should stay CONDITIONAL.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. III.B defines f_A and f_S through the nonlocal matrix elements (21), whose Lorentz factors P_A^μ and P_S contain the physical mass M. Equations (29)–(34) then replace M by M0 inside the integrand to obtain f_S and f_A, and Eq. (35) is 'verified' by moving the same replacement into the mixing coefficient α=2m/M0 (Eqs. 36–38). The paper states that using the fixed physical value α=2m/M leads to inconsistencies across current components, so the M0 choice is what makes the relations hold. This is a consistency prescription, not a derivation: the only antecedent is the Type II link of Ref. [30] (Eq. 49). If that link is not a faithful zero-mode resummation, then the claimed equality f_A = f∥ − 2m/M0 f⊥ and f_S = f⊥ − 2m/M0 f∥ redefines the QCD EOM relation (35) rather than testing it, and the 'self-consistency' is partially circular. The framework's valence-only Fock space and Gaussian radial wave function also force the twist-4 DAs to coincide with twist-2 (Sec. IV), which limits the advertised 'beyond leading twist' content, but the M0 prescription is the more load-bearing condition for the central extraction claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the authors' 'self-consistent light-front quark model' based on the Bakamjian-Thomas (BT) construction to compute the rho meson decay constants and distribution amplitudes (DAs) up to twist 4. The central technical claim is that the decay constants f_parallel and f_perp, together with the twist-3 and twist-4 chiral-even and chiral-odd DAs, can be extracted from local and nonlocal matrix elements of the vector, tensor, axial-vector, and scalar currents in a way that is independent of current component, polarization, and reference frame, provided the physical meson mass M is replaced by the internal invariant mass M0 inside all integrands, including the mixing coefficient alpha = 2m/M0. The paper presents explicit integral formulas and numerical results: f_parallel = 215 MeV (in agreement with experiment), f_perp = 173 MeV, f_A = 122 MeV, f_S = 61 MeV, and a full set of DAs with their xi-moments and Gegenbauer moments. The authors also compare their DAs with QCD sum rules, DSE, and lattice results.","tokens_in":23442,"tokens_out":7844,"duration_ms":71675,"significance":"If the framework is accepted, the paper would provide a complete, internally consistent set of rho meson DAs up to twist 4 from a simple valence-quark model, with the notable feature that the previously problematic axial-vector and scalar current extractions are regularized by the BT mass substitution. The numerical prediction for f_parallel matches experiment, and the model yields the SU(6) relation f_pi + f_parallel = 2 f_perp. The paper is transparent about its assumptions and provides explicit integral formulas that allow independent numerical reproduction, although no code is supplied. The main caveat is that the advertised 'beyond leading twist' content is partially diminished by the model's prediction that the twist-4 DAs coincide with twist-2 DAs in the equal-mass case, and the verification of the QCD equation-of-motion relations is contingent on the same M-to-M0 prescription used to define the extracted constants.","major_comments":[{"comment":"The demonstration that the model satisfies the QCD equation-of-motion relations (35) is circular. The tilde quantities \\tilde f_\\perp and \\tilde f_\\parallel are defined with the factor 2m/M0 evaluated inside the integrals, and f_A and f_S are computed after the same replacement in Eqs. (30) and (34). Consequently, the equalities f_A = f_\\parallel - \\tilde f_\\perp and f_S = f_\\perp - \\tilde f_\\parallel hold by construction; Table II does not constitute an independent verification of Eq. (35). The text's statement that \"to ensure that Eq. (35) is satisfied self-consistently within our model, we should take M as M0\" makes the desired outcome an input. The authors should either derive the M0 substitution from the BT construction or explicitly describe Eq. (35) as a consistency constraint used to fix the mixing coefficient, not as a benchmark being tested.","section":"Sec. III.B, Eqs. (36)-(38), Table II"},{"comment":"The uniform replacement of the physical mass M by the invariant mass M0 inside the integrands is an ad hoc prescription, inherited from Ref. [30] but not derived in this paper. Because the claimed current-component and polarization independence of f_A and f_S depends entirely on this replacement (the paper notes that using a fixed alpha = 2m/M would lead to inconsistencies across current components), the central extraction claim rests on an input assumption. Please provide a detailed derivation of this substitution from the covariant BS model or from the BT mass operator, or at least a discussion of its range of validity. Without this, the 'self-consistency' is a property of the prescription rather than a result of the model.","section":"Sec. III.B, Eqs. (29)-(34)"},{"comment":"The title and abstract advertise 'beyond leading twist' DAs, but in the equal-mass case the twist-4 DAs \\phi^\\parallel_{4;V} and \\phi^\\perp_{4;T} coincide identically with their twist-2 counterparts because the corresponding operators in Table I are the same. The text calls this 'fortuitous,' yet it means that the only genuinely higher-twist predictions are the twist-3 DAs. This limitation should be stated prominently in the abstract and conclusions, or the model should be extended (e.g., with unequal masses or higher Fock states) to generate distinct twist-4 DAs, before the claim of a comprehensive twist-4 analysis can be sustained.","section":"Sec. IV, Table III and Eqs. (45)-(53)"}],"minor_comments":[{"comment":"The word 'chrial-even' should be 'chiral-even'; in Table IV, 'Gegenbaur' should be 'Gegenbauer'.","section":"Sec. IV.B"},{"comment":"The symbol M is used both for the meson mass and for the combination M2 - M1 in the spin-orbit wave functions; please rename the latter (e.g., \\Delta M) to avoid confusion.","section":"Eq. (10)"},{"comment":"The notation for the decay constants is inconsistent: the text uses both f_A^\\rho and f_\\rho^A (similarly for f_S); please adopt a single convention.","section":"Throughout"},{"comment":"The legends in Figure 1 ('Ball99', 'Gao14', 'Fors10') do not match the references cited in the text (QCD SRs [14], DSE [67], HERA-fit [68]); please align the labels with the bibliography.","section":"Fig. 1"},{"comment":"The statement about the heavy-quark limit that 'all DAs converge to a common form proportional to \\delta(x - 1/2)' is not fully supported by Eq. (56) alone, since the DAs also involve helicity-dependent spin factors from Table I that may introduce further x dependence; please clarify the argument.","section":"Sec. IV.C"},{"comment":"The 'Expt.' row lists values for f^\\parallel but not for f^\\perp; a reference for an experimental or lattice determination of f^\\perp would be useful for comparison.","section":"Table II"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid continuation of the authors' established LFQM program, and the numerics are internally consistent. The main concern is the circularity of the QCD equation-of-motion 'verification': if the authors can reframe this as a consistency constraint and be more explicit about the ad hoc nature of the M-to-M0 prescription, the paper would be acceptable. The equality of twist-2 and twist-4 DAs should also be highlighted as a limitation. The paper fits the scope of the journal, but the load-bearing status of the M0 prescription requires attention before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThis is a careful, internally coherent extension of the authors' self-consistent light-front quark model, and it delivers more than the abstract alone suggests: the extraction of f_A and f_S from nonlocal axial-vector and scalar currents, the chiral-odd twist-3 DAs, and the explicit resolution of the f_parallel/f_perp mixing through the Bakamjian-Thomas construction. The numerics close (f_A = 215 - 93 = 122 MeV, f_S = 173 - 112 = 61 MeV), f_parallel matches experiment, and the chiral-limit shapes are the expected asymptotic forms. The paper is transparent about its machinery: the M -> M0 replacement is flagged as a BT requirement, and the authors admit that frame independence holds by construction.\n\nThe soft spots are real, though not fatal. The advertised \"beyond leading twist\" content is thinner than the title implies: the twist-4 DAs coincide with twist-2 because of the valence-only Gaussian wave function, so there is no genuinely new twist-4 information here. The verification of the QCD equation-of-motion relations (Eq. 35) is the weakest link: alpha is set to 2m/M0 inside the integrals, and the paper notes that using the physical alpha = 2m/M would break current-component independence. That means the relation is being imposed by the model's prescription rather than independently tested. I'd like to see the authors check whether Eq. (35) with alpha = 2m/M0 is an algebraic identity forced by the spinor structure; if it is, the \"verification\" is a consistency check, not a prediction. The disagreement with QCD sum rules for twist-3 DAs is model dependence that deserves a more prominent caveat, and the absence of uncertainty estimates is a practical nuisance for anyone who wants to use the numbers.\n\nBottom line: this deserves a serious referee, but the right verdict is conditional, not clean accept. I'd ask for a revision that states the twist-2/twist-4 degeneracy up front, tests whether the EOM relation is an identity in this model, and adds at least a rough uncertainty estimate. If the authors clarify the status of the M0 substitution—or even if they don't—this is a useful reference for vector-meson phenomenology. Send it to review.","headline":"Careful, internally coherent LFQM extension delivering f_A, f_S, and chiral-odd DAs for the rho, but twist-4 collapses to twist-2 and the EOM check leans on the M→M0 prescription.","tokens_in":24144,"tokens_out":3236,"would_cite":true,"duration_ms":29003,"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":"This paper claims that replacing the physical rho mass by the internal quark-antiquark mass in light-front integrals yields decay constants and distribution amplitudes through twist 4 that are independent of current component…","keywords":["rho meson","decay constants","distribution amplitudes","light-front quark model","Bakamjian-Thomas construction","higher twist","chiral-even and chiral-odd","light-front zero modes"],"falsifier":"A calculation using the fixed physical value $\\alpha = 2m/M$ outside the integrals that still yields current-component-independent $f_A$ and $f_S$ would undercut the paper's central justification; alternatively, a lattice QCD determination of $f_A$ and $f_S$ at a scale near 1 GeV that disagrees with 122 MeV and 61 MeV, or with the equation-of-motion mixing relations evaluated at $M_0$, would falsify the extraction prescription.","tokens_in":22911,"feed_emoji":"⚛️","tokens_out":14284,"duration_ms":115668,"temperature":0.7,"pith_summary":"This paper aims to establish that the rho meson's two decay constants and its eight distribution amplitudes (DAs, the functions describing how the quark and antiquark share the meson's longitudinal momentum) through twist 4 can be extracted from light-front matrix elements in a way that is independent of the chosen current component, polarization, and reference frame. The key claim is that this holds when the Bakamjian-Thomas construction is applied uniformly: the physical rho mass $M$ is replaced by the quark-antiquark invariant mass $M_0$ everywhere inside the integrals, including in the mixing coefficient $\\alpha = 2m/M_0$. With that substitution, the axial-vector and scalar current extractions give $f_A = 122$ MeV and $f_S = 61$ MeV, and the QCD equation-of-motion relations $f_A = f_\\parallel - (2m/M_0) f_\\perp$ and $f_S = f_\\perp - (2m/M_0) f_\\parallel$ are satisfied. If the paper is right, a single light-front quark model supplies a complete, internally consistent set of rho DAs through twist 4, useful as nonperturbative input for exclusive processes and for comparisons with lattice QCD.","feed_headline":"Rho meson decay constants and DAs turn out frame-independent","feed_subtitle":"Replacing rho mass by quark-antiquark invariant mass inside integrals gives f∥=215 MeV and the full twist-4 DA set","key_machinery":"The central object is the Bakamjian-Thomas (BT) construction, implemented as the uniform substitution $M \\to M_0$, from the physical rho mass to the internal invariant mass of the quark-antiquark pair, applied both to the Lorentz prefactors and inside the integrands of the matrix elements, including the mixing coefficient $\\alpha = 2m/M_0$. This substitution absorbs the light-front zero-mode contamination that otherwise appears in the 'bad' current components and is what makes the extracted decay constants and DAs independent of current components, polarizations, and frames. The nonlocal axial-vector and scalar matrix elements are handled by integrating over $z^-$; the derivative of the delta function in Eq. (24) converts the $x$-moment into the integrated forms for $f_A$ and $f_S$. The Gaussian radial wave function with a single variational parameter $\\beta$ supplies the nonperturbative profile, and in the equal-mass case it makes the twist-2 and twist-4 DAs coincide by construction.","core_discovery":"The central claim is that in the standard light-front quark model built on the Bakamjian-Thomas construction, the rho-meson decay constants $f_\\parallel$ and $f_\\perp$ and the chiral-even DAs $(\\phi^\\parallel_{2;\\mathrm{V}}, \\phi^\\parallel_{3;\\mathrm{V}}, \\psi^\\perp_{3;\\mathrm{A}}, \\phi^\\parallel_{4;\\mathrm{V}})$ and chiral-odd DAs $(\\phi^\\perp_{2;\\mathrm{T}}, \\phi^\\perp_{3;\\mathrm{T}}, \\psi^\\parallel_{3;\\mathrm{S}}, \\phi^\\perp_{4;\\mathrm{T}})$ are all obtainable from the local and nonlocal matrix elements $\\langle 0|\\bar{q}(z)\\,\\Gamma\\, q(-z)|\\rho(P,h)\\rangle$ with $\\Gamma = (\\gamma^\\mu, \\sigma^{\\mu\\nu}, \\gamma^\\mu\\gamma_5, \\mathbf{1})$, and the results do not depend on current components, polarizations, or frames. The novelty lies in the axial-vector and scalar channels, where $f_\\parallel$ and $f_\\perp$ mix; the paper shows the mixing is resolved only when the physical mass $M$ is replaced by the internal invariant mass $M_0$ inside the integrals, including in $\\alpha = 2m/M_0$. The resulting values $f_\\parallel = 215$ MeV, $f_\\perp = 173$ MeV, $f_A = 122$ MeV, and $f_S = 61$ MeV satisfy the QCD equation-of-motion mixing relations written above. In the chiral limit all twist-2 and twist-4 DAs (and the two twist-3 partners $\\psi^\\perp_{3;\\mathrm{A}}$ and $\\psi^\\parallel_{3;\\mathrm{S}}$) reduce to $6x(1-x)$; $\\phi^\\parallel_{3;\\mathrm{V}}$ and $\\phi^\\perp_{3;\\mathrm{T}}$ approach $\\frac{3}{4}(1+\\xi^2)$ and $3\\xi^2$, consistent with QCD sum-rule expectations.","pith_inferences":["This suggests the same $M \\to M_0$ substitution should yield frame-independent results for unequal-mass vector mesons ($K^*$, $\\phi$, $D^*$, $B^*$), for which the paper itself notes that the twist-2 and twist-4 DAs no longer coincide; those systems would be a sharper test of the prescription.","A lattice QCD computation of $f_A$ and $f_S$ at a scale near 1 GeV would independently test the 122 MeV and 61 MeV predictions and the $M_0$-evaluated mixing relations.","One could compare the BT-substitution results with an explicit Bethe-Salpeter light-front calculation that includes zero-mode diagrams; agreement would indicate the substitution is equivalent to resumming those contributions."],"forward_implications":["The decay constants $f_\\parallel = 215$ MeV, $f_\\perp = 173$ MeV, $f_A = 122$ MeV, and $f_S = 61$ MeV are concrete predictions; $f_\\parallel$ already matches the experimental value $216(5)$ MeV from $e^+e^-$ annihilation.","The full twist-4 set of chiral-even and chiral-odd DAs gives the nonperturbative input for collinear-factorization calculations of exclusive rho production, radiative, and semileptonic processes.","Satisfying the equation-of-motion mixing relations with $\\alpha = 2m/M_0$ inside the integrals confirms that the model reproduces the QCD constraint without residual zero-mode contamination.","In the chiral limit the DAs reduce to the expected asymptotic shapes $6x(1-x)$, $\\frac{3}{4}(1+\\xi^2)$, and $3\\xi^2$, providing a consistency check against QCD sum rules."],"supporting_citations":[{"why":"First proposed the Type II link, i.e., replacing M by M0 in vector-meson decay constant analyses; this paper extends that prescription to axial-vector and scalar currents.","marker":"[30]"},{"why":"Previously demonstrated current-component, polarization, and frame independence for f_parallel and f_perp in the same model; the basis for the present extension.","marker":"[39]"},{"why":"Supplied the z^- integration identity and the twist-4 decay-constant analysis for pseudoscalar mesons that is adapted here for the rho.","marker":"[60]"},{"why":"Defined the twist-3 DA formalism and the QCD equation-of-motion mixing relations the paper verifies.","marker":"[13]"},{"why":"Provided the twist-4 parametrizations of the nonlocal vector and tensor matrix elements used to assign the DAs to twists 2 through 4.","marker":"[16]"},{"why":"The Bakamjian-Thomas construction itself, which ensures the Poincare algebra for the mass-operator formulation.","marker":"[58]"},{"why":"World-average experimental value f_parallel = 216(5) MeV used to benchmark the predicted 215 MeV.","marker":"[65]"}],"fun_headline_variants":["Rho decay constants: frame-independent in light-front model","Self-consistent rho DAs: all twists, no frame dependence","Mixing resolved: rho f_∥ and f_⊥ from nonlocal currents","Bakamjian-Thomas model yields full twist-4 rho DAs","Rho meson: decay constants and DAs without frame bias"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the physical rho mass $M$ can be replaced by the internal quark-antiquark mass $M_0$ uniformly inside the integrals, including inside the mixing coefficient $\\alpha = 2m/M_0$; the paper justifies this by internal consistency rather than deriving it from QCD.","fun_headline_variants_meta":{"raw":{"variants":["Rho decay constants: frame-independent in light-front model","Self-consistent rho DAs: all twists, no frame dependence","Mixing resolved: rho f_∥ and f_⊥ from nonlocal currents","Bakamjian-Thomas model yields full twist-4 rho DAs","Rho meson: decay constants and DAs without frame bias"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000248,"raw_usage":{"total_tokens":1782,"prompt_tokens":1417,"completion_tokens":365,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":1033,"completion_tokens_details":{"reasoning_tokens":268}},"tokens_in":1033,"tokens_out":365,"duration_ms":4036,"temperature":1.0,"reasoning_tokens":268,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:17:02.001836+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A calculation using the fixed physical value $\\alpha = 2m/M$ outside the integrals that still yields current-component-independent $f_A$ and $f_S$ would undercut the paper's central justification; alternatively, a lattice QCD determination of $f_A$ and $f_S$ at a scale near 1 GeV that disagrees with 122 MeV and 61 MeV, or with the equation-of-motion mixing relations evaluated at $M_0$, would falsify the extraction prescription.","supporting_citations":[{"cited_title":"Hwang, Analyses of decay constants and light-cone distribution amplitudes for s-wave heavy meson, Phys","cited_arxiv_id":null,"evidence_quote":"First proposed the Type II link, i.e., replacing M by M0 in vector-meson decay constant analyses; this paper extends that prescription to axial-vector and scalar currents."},{"cited_title":"Dhiman, H","cited_arxiv_id":null,"evidence_quote":"Previously demonstrated current-component, polarization, and frame independence for f_parallel and f_perp in the same model; the basis for the present extension."},{"cited_title":"Bakamjian and L","cited_arxiv_id":null,"evidence_quote":"Supplied the z^- integration identity and the twist-4 decay-constant analysis for pseudoscalar mesons that is adapted here for the rho."},{"cited_title":"Beneke, G","cited_arxiv_id":null,"evidence_quote":"Defined the twist-3 DA formalism and the QCD equation-of-motion mixing relations the paper verifies."},{"cited_title":"Ball and V","cited_arxiv_id":null,"evidence_quote":"Provided the twist-4 parametrizations of the nonlocal vector and tensor matrix elements used to assign the DAs to twists 2 through 4."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The Bakamjian-Thomas construction itself, which ensures the Poincare algebra for the mass-operator formulation."},{"cited_title":"Cheng, C.-K","cited_arxiv_id":null,"evidence_quote":"World-average experimental value f_parallel = 216(5) MeV used to benchmark the predicted 215 MeV."}],"review_version":1}