{"id":"8d5c6402-68ac-423c-bb56-0e9165a9cad0","arxiv_id":"1908.08657","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Analytic pion TMDs for f1 and Boer-Mulders are derived in a spin-improved light-front holographic model, and f1 is evolved from 0.316 GeV to 1 GeV.","lead":"This paper computes the two leading-twist transverse-momentum-dependent parton distributions of the pion in a light-front holographic model with spin-improved wave functions. These analytic curves are candidate model inputs for future pion Drell-Yan and semi-inclusive scattering analyses.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Direct substitution of the stated LFWFs into Eq. (12) does not reproduce Eq. (13); the central unpolarized TMD is missing an x^3(1-x)^3 factor and is unverified as printed.","rationale":"The reader's conditional verdict already notes an unexplained mismatch between the wave-function mapping and the final unpolarized TMD. This stress-test makes the mismatch concrete: direct evaluation of the printed overlap, Eq. (12), from the printed LFWFs, Eqs. (4) and (7), produces an expression that differs from Eq. (13) by a factor 1/[x^3(1-x)^3]. Since Eq. (13) feeds every quantitative statement of the paper, including the comparison plots, the positivity check, and the TMD evolution in Sec. IV, the central claim is not internally consistent as printed. This is independent of the reader's stated weakest assumption about the model-scale parameters: even if kappa = 523 MeV, m = 330 MeV, and M_pi = 139 MeV are accepted, the derivation of the central formula does not close. The issue is fixable by correcting either the LFWF mapping or the analytic result and then rechecking the figures, so the appropriate disposition remains conditional rather than outright rejection. I mark agreement as partial because the reader's formal weakest-assumption field points to parameter inheritance, while the same concern appears in the reader's rationale; the load-bearing defect is the algebraic mismatch.","tokens_in":11983,"tokens_out":14378,"duration_ms":140849,"concrete_test":"Use a computer algebra system (SymPy or Mathematica) to substitute Eq. (4) and Eq. (7) into Eq. (12) and simplify. If the result is not Eq. (13) but rather N0^2/(pi*kappa^2) * [k_perp^2 + (m + M_pi x(1-x))^2] * exp[-(k_perp^2 + m^2)/(kappa^2 x(1-x))], then the printed formula and its figures are unsupported by the stated wave functions. As a cross-check, integrate both candidate expressions over k_perp and x with the N0 used in Sec. III: with a fixed N0 they cannot both satisfy the normalization condition of Eq. (14) unless an unstated x-dependent factor is introduced, and the text states none.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central unpolarized TMD formula, Eq. (13), is the basis for Figs. 1-2, the positivity check, and the evolution in Sec. IV, but it does not follow from the wave functions stated in the paper. Inserting Eq. (4) and Eq. (7) into the overlap Eq. (12) gives, after collecting constants, f1 = N0^2/(pi*kappa^2) * [k_perp^2 + (m + M_pi x(1-x))^2] * exp[-(k_perp^2 + m^2)/(kappa^2 x(1-x))]. The printed Eq. (13) instead carries the prefactor N0^2/(pi*kappa^2 x^3(1-x)^3). No step in Sec. III introduces such an x-dependent denominator; its only possible origin is the 1/[x(1-x)] denominator in the spin wave function (Eq. 6), which is absent from the mapped LFWFs in Eq. (7). Thus either the LFWF mapping or the analytic TMD, and therefore all numerical results and comparisons, is wrong as printed. This is a derivation-level check, not a disagreement about model choices, so it cannot be repaired by adjusting kappa, m, M_pi, or the comparison set.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes the leading-twist unpolarized TMD f1π(x,k⊥²) and the Boer-Mulders function h1π⊥(x,k⊥²) for the pion using light-front holographic wave functions with a spin-improved spin wave function. The authors quote analytic expressions in Eqs. (13) and (18), compare their results with light-front constituent and soft-wall AdS/QCD models, check the positivity bound, and perform a leading-order TMD evolution of f1 from a model scale Q0=0.316 GeV to a higher scale. The parameters κ, m, and Mπ are taken from earlier fits to pion properties. The central new objects are the closed-form TMDs and the evolved f1.","tokens_in":12325,"tokens_out":8241,"duration_ms":76996,"significance":"If the formulas were correct, the paper would provide simple analytic pion TMDs with parameters fixed by the pion decay constant, charge radius, and form factor, extending the light-front holographic program to transverse-momentum-dependent structure. The positivity check and the comparison with existing model calculations are conceptually useful. However, the main unpolarized result does not follow from the stated wave functions as written, and the Boer-Mulders derivation omits essential specifications; the quantitative content of the paper is therefore not currently reliable.","major_comments":[{"comment":"The central formula Eq. (13) does not follow from the preceding equations. Direct substitution of Eq. (4) and Eq. (7) into Eq. (12) gives |Ψπ^(0)|² = (8π²N0²/κ²)(m - Mπ x(1-x))² exp[-(k⊥²+m²)/(κ²x(1-x))] and k⊥²|Ψπ^(1)|² = (8π²N0²/κ²) k⊥² exp[-(k⊥²+m²)/(κ²x(1-x))], so that f1π(x,k⊥²) = N0²/(πκ²)[k⊥² + (m - Mπ x(1-x))²] exp[-(k⊥²+m²)/(κ²x(1-x))]. The denominator 1/[x³(1-x)³] in Eq. (13) is not produced by this overlap; it appears to be inherited from the spin wave function in Eq. (6) rather than from the mapped LFWFs in Eq. (7). Since Eq. (13) underlies Figs. 1, 2, 5, and 6, the normalization condition Eq. (14), and the qualitative conclusions, all numerical results and model comparisons must be recomputed after correcting this prefactor.","section":"Section III.A, Eqs. (4), (7), (12), (13)"},{"comment":"The derivation of the Boer-Mulders function is not reproducible as written. Equation (17) contains Ψ^(0)(x,k⊥) and Ψ^(1)(x',k'⊥) with k'⊥ = k⊥ - q⊥, but x' is never defined and no relation between x and x' is stated; the q⊥ integral is not evaluated in the text; and the step from Eq. (17) to Eq. (18) is not shown. In addition, the statement g² = 4πα_s(μ0) is not accompanied by any numerical value for α_s(μ0) or g², although this parameter controls the overall magnitude in Figs. 3-5 and in the positivity check in Fig. 5. Please specify x', carry out or reference the q⊥ integration, and quote the coupling value used.","section":"Section III.B, Eqs. (17)-(18)"},{"comment":"The evolution formula Eq. (20) is incomplete as stated. The TMD is written without a rapidity scale (ζ), the Sudakov exponent ~S([b*; μ_b, μ) is not defined explicitly, and the μ0 dependence appears only through the g_K ln(μ/μ0) term. Without specifying the rapidity evolution and the exact form of ~S, Eq. (20) cannot be evaluated and Fig. 6 is not reproducible. Please supply the missing definitions or cite the exact convention used in Refs. [71-74].","section":"Section IV, Eq. (20) and Fig. 6"}],"minor_comments":[{"comment":"The normalization constant N0 is never specified; it should be stated explicitly, or its determination from Eq. (14) should be shown, so that the plotted curves can be reproduced.","section":"Section III.A, Eq. (13)"},{"comment":"The caption of Fig. 6 states that the evolution is performed up to Q = 1 GeV, while the labels embedded in the figure show Q0 = 0.316 and Q = 5 GeV, with 'g1 = 0.13' instead of 'g2 = 0.13'; this inconsistency should be resolved.","section":"Section IV, Fig. 6"},{"comment":"The text says the unpolarized TMD is plotted as a function of k⊥², but the horizontal axis of Fig. 6 is labeled k⊥ [GeV]; the vertical axis is labeled k⊥ f(...) while Eq. (13) gives f(...). Please make the axes and the labels consistent with the quantity shown.","section":"Section IV, Fig. 6"},{"comment":"The spin wave function components in Eq. (6) contain explicit 1/[x(1-x)] factors, while the mapped LFWFs in Eq. (7) do not; this discrepancy is likely the source of the erroneous prefactor in Eq. (13) and should be clarified or removed.","section":"Section II, Eq. (6) and Eq. (7)"}],"recommendation":"major_revision","confidential_remarks":"The paper's main formulas cannot be relied upon without a full re-derivation and re-plotting: Eq. (13) does not follow from the stated LFWFs, and Eq. (18) lacks essential specifications. The issues appear to be technical slips rather than intentional misrepresentation, but they affect every numerical result, so a major revision is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. The paper is a genuine new application of the Ahmady–Mondal–Sandapen spin-improved light-front holographic wave functions to pion TMDs; nobody had computed f1 and h1^\\perp with these wave functions before. But as printed, the central unpolarized TMD does not follow from its own inputs. Direct substitution of Eq. (7) into Eq. (12) gives f1 = N0^2/(pi kappa^2) [k^2 + (m + M_pi x(1-x))^2] exp[-(k^2+m^2)/(kappa^2 x(1-x))], with no x^3(1-x)^3 denominator. Eq. (13) carries exactly that denominator, and no stated step in Section III produces it. The stress-test note is right.\n\nThe positive side: the model setup is coherent, the physical motivation is reasonable, and the comparisons with the light-front constituent model and the soft-wall AdS/QCD model provide useful context. A corrected version would probably still show the same qualitative physics—for fixed x, the k-perp dependence is unchanged; what changes is the x shape, so Figs. 1–2 and the normalization condition need to be redone.\n\nThe soft spots go beyond the main equation. The Boer-Mulders derivation is under-specified: Eq. (17) writes x' without ever saying x' = x, the q-perp integral is not shown, and g^2 is defined through alpha_s but no value is given. The evolution section is a sketch: the model scale 0.316 GeV is asserted rather than derived, the Sudakov form factor is written only formally, and the transformation from b space back to k space is omitted. There is also a caption/body mismatch: the text says evolution to Q = 1 GeV, the figure caption says Q = 5 GeV, and one panel label says g1 while the text says g2.\n\nI would not desk-reject this. The errors are algebraic and presentational, not a broken physical idea, and a referee could require the complete derivation, parameter values, and corrected figures. But in the current form the numerical results are unverified and I wouldn't cite it. For a reading group, it's useful mainly as a reminder to check formulas by substitution.","headline":"The paper's model application is new and worth a look, but the central unpolarized TMD as printed does not follow from its own wave functions, so all numerical results inherit an unverified x-dependent prefactor.","tokens_in":12906,"tokens_out":8308,"would_cite":false,"duration_ms":80956,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The light-front holographic model, with a spin-improved wave function, predicts analytic pion TMDs and a measurable Sudakov broadening of $f_{1\\pi}$ between $Q_0=0.316$ GeV and $Q=1$ GeV.","keywords":["pion","transverse momentum dependent parton distributions","TMDs","Boer-Mulders function","light-front holographic model","light-front wave functions","Sudakov evolution","leading twist"],"falsifier":"Measure the unpolarized pion TMD through pion-induced Drell-Yan or SIDIS at $Q\\approx 1$ GeV and compare the evolved $f_{1\\pi}(x=0.3,k_\\perp)$ from Eq. (13): the model predicts the $k_\\perp$ peak moves from $0.2$ GeV at $Q_0=0.316$ GeV to roughly $0.4$ GeV ($g_2=0.09$) or $0.5$ GeV ($g_2=0.13$) and the distribution broadens. If the observed peak location or the double-peak shape at small $k_\\perp$ is absent, the spin-improved holographic input or the $g_2$ dependence of the evolution is ruled out.","tokens_in":11754,"feed_emoji":"⚛️","tokens_out":7079,"duration_ms":61059,"temperature":0.7,"pith_summary":"This paper claims that the light-front holographic model, upgraded with a spin-dependent wave function, gives analytic closed forms for the two leading-twist transverse-momentum-dependent parton distributions (TMDs) of the pion: the unpolarized distribution $f_{1\\pi}(x,k_\\perp^2)$ and the Boer-Mulders function $h_{1\\pi}^\\perp(x,k_\\perp^2)$. These are the functions that map where quarks sit inside a pion in both longitudinal momentum fraction and transverse momentum. A sympathetic reader should care because pion TMDs are needed to interpret semi-inclusive deep inelastic scattering and Drell-Yan data, and model predictions with a simple analytic form can be evolved to experimental scales and compared. The paper also performs a leading-order Sudakov evolution of $f_{1\\pi}$ from model scale $Q_0=0.316$ GeV to $Q=1$ GeV and shows the distribution broadens in transverse momentum, with the peak shifting from $k_\\perp\\approx 0.2$ GeV to roughly $0.4$-$0.5$ GeV depending on the nonperturbative parameter $g_2$.","feed_headline":"Pion quark distributions derived in closed form","feed_subtitle":"Two leading-twist pion momentum distributions are computed analytically and evolved to scales experiments can reach.","key_machinery":"The load-bearing object is the spin-improved light-front holographic wave function of the pion, built as the product of the holographic momentum-space wave function $\\psi_\\pi(x,k_\\perp)=\\frac{4\\pi N_0}{\\kappa}\\sqrt{x(1-x)}\\exp[-(k_\\perp^2+m^2)/(2\\kappa^2 x(1-x))]$ and a spin wave function $\\varphi_\\pi(\\lambda_1,\\lambda_2)$ with parameters $A=B=1$. From this wave function the paper obtains two pion light-front wave-function amplitudes, $\\Psi_\\pi^{(0)}$ for $L_z=0$ and $\\Psi_\\pi^{(1)}$ for $|L_z|=1$, whose overlap yields the unpolarized TMD in Eq. (13); the Boer-Mulders function in Eq. (18) follows from expanding the Wilson line to one gluon exchange between the struck quark and spectator. The evolution machinery is the $b_\\perp$-space TMD evolution formula with the Sudakov form factor, using $b_*(b_\\perp)$ and a nonperturbative function $g_K(b_\\perp)=-g_2 b_\\perp^2/2$ with $g_2=0.09$ or $0.13$.","core_discovery":"Using the spin-improved light-front holographic wave function, which includes the $L_z=0$ and $|L_z|=1$ orbital angular momentum components of the valence $|q\\bar q\\rangle$ Fock state, the paper derives the unpolarized pion TMD as $f_{1\\pi}(x,k_\\perp^2)$ in Eq. (13) and the Boer-Mulders function as $h_{1\\pi}^\\perp(x,k_\\perp^2)$ in Eq. (18). The $f_{1\\pi}$ expression is an overlap of the two light-front wave-function amplitudes and reduces to a Gaussian in $k_\\perp$ modulated by $[k_\\perp^2+(m+x(1-x)M_\\pi)^2]$; the Boer-Mulders function comes from a one-gluon-exchange final-state interaction through the Wilson line and satisfies the model-independent positivity bound $f_{1\\pi}\\ge (k_\\perp/M_\\pi)|h_{1\\pi}^\\perp|$ at all $x$ and $k_\\perp$. Both distributions are symmetric under $x\\leftrightarrow 1-x$, yielding a double-peak structure at small $k_\\perp$ that merges into a single Gaussian at larger $k_\\perp$, and the evolution of $f_{1\\pi}$ from $Q_0=0.316$ GeV to $Q=1$ GeV in $b_\\perp$-space exhibits Sudakov broadening whose size depends on $g_2$.","pith_inferences":["The same spin-improved holographic wave function could be applied to other pseudoscalars such as the kaon by changing quark masses; if the $A=B=1$ spin structure is universal, one would predict analogous double-peak TMDs at the kaon's lower model scale.","The double-peak signature at small $k_\\perp$ is essentially a valence-structure effect inherited from the $x\\leftrightarrow 1-x$ symmetry; one might test whether any observed pion TMD asymmetry in $z$ or $Q^2$ tracks this symmetry before invoking higher Fock states.","The evolution from such a low model scale ($0.316$ GeV) assumes the TMD is purely nonperturbative at that scale; a lattice QCD extraction of the pion TMD at moderate virtualities could settle whether the $g_2=0.09$-$0.13$ range is realistic or merely tuned."],"forward_implications":["The analytic $f_{1\\pi}$ and Boer-Mulders expressions can be used directly as model-scale inputs for phenomenological studies of pion SIDIS and Drell-Yan azimuthal asymmetries.","The model satisfies the positivity bound $h_{1\\pi}^\\perp$ at all $x$ and $k_\\perp$, so it provides a controlled template where the unpolarized quark probability always dominates the transversely polarized one.","The leading-order Sudakov evolution predicts a measurable broadening: at $x=0.3$ the peak of $k_\\perp f_{1\\pi}$ moves from $0.2$ GeV to $0.4$-$0.5$ GeV when $Q$ rises from $0.316$ GeV to $1$ GeV.","Because $f_{1\\pi}$ and $h_{1\\pi}^\\perp$ are symmetric under $x\\leftrightarrow 1-x$, the model predicts a double peak in $x$ at small $k_\\perp$ and a single Gaussian-shaped peak at larger $k_\\perp$; future data at low transverse momentum can test this shape.","At large $k_\\perp$ the model's $f_{1\\pi}$ converges with light-front constituent and soft-wall AdS/QCD results, suggesting the high-transverse-momentum tail is insensitive to the detailed wave function."],"supporting_citations":[{"why":"Supplies the spin-improved holographic wave function and the parameters $\\kappa=523$ MeV, $m=330$ MeV, $A=B=1$ used in Eqs. (13) and (18).","marker":"[66]"},{"why":"Earlier version of the spin-improved wave function with $A=1$, $B=1$; motivates the spin structure and parameter choice.","marker":"[65]"},{"why":"Light-front constituent model results for $f_{1\\pi}$ and $h_{1\\pi}^\\perp$ used for comparison, and the one-gluon-exchange treatment of the Boer-Mulders function.","marker":"[47]"},{"why":"Soft-wall AdS/QCD pion TMD results used as the second comparison for $f_{1\\pi}$.","marker":"[55]"},{"why":"Provides the nonperturbative $g_2$ values $0.09$ and $0.13$ and the $b$-space TMD evolution framework used in Section V.","marker":"[54]"},{"why":"Defines the general pion light-front state with spin wave function, used to set up the $L_z=0$ and $|L_z|=1$ amplitudes.","marker":"[67]"},{"why":"Supplies the $b_\\perp$-space TMD evolution formula with Sudakov form factor used to evolve $f_{1\\pi}$ to $Q=1$ GeV.","marker":"[71]"}],"fun_headline_variants":["Pion TMDs get closed-form holographic derivation","Holographic model yields pion quark momentum distributions","Pion TMDs: spin-improved holographic wave function delivers exact forms","Closed-form pion TMDs from light-front holography","Pion Boer-Mulders function derived analytically"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole prediction inherits its absolute size and shape from an earlier fit of the model parameters $\\kappa=523$ MeV, $m=330$ MeV, $M_\\pi=139$ MeV and $A=B=1$ to pion decay constant, charge radius, and form factor; no pion TMD data anchor these curves, so if that fit is not valid for TMDs at leading twist the predicted normalization and $x,k_\\perp$ dependence shift accordingly.","fun_headline_variants_meta":{"raw":{"variants":["Pion TMDs get closed-form holographic derivation","Holographic model yields pion quark momentum distributions","Pion TMDs: spin-improved holographic wave function delivers exact forms","Closed-form pion TMDs from light-front holography","Pion Boer-Mulders function derived analytically"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000562,"raw_usage":{"total_tokens":2683,"prompt_tokens":978,"completion_tokens":1705,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":594,"completion_tokens_details":{"reasoning_tokens":1621}},"tokens_in":594,"tokens_out":1705,"duration_ms":11228,"temperature":1.0,"reasoning_tokens":1621,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:35:10.133615+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the unpolarized pion TMD through pion-induced Drell-Yan or SIDIS at $Q\\approx 1$ GeV and compare the evolved $f_{1\\pi}(x=0.3,k_\\perp)$ from Eq. (13): the model predicts the $k_\\perp$ peak moves from $0.2$ GeV at $Q_0=0.316$ GeV to roughly $0.4$ GeV ($g_2=0.09$) or $0.5$ GeV ($g_2=0.13$) and the distribution broadens. If the observed peak location or the double-peak shape at small $k_\\perp$ is absent, the spin-improved holographic input or the $g_2$ dependence of the evolution is ruled out.","supporting_citations":[{"cited_title":"Ahmady, C","cited_arxiv_id":null,"evidence_quote":"Supplies the spin-improved holographic wave function and the parameters $\\kappa=523$ MeV, $m=330$ MeV, $A=B=1$ used in Eqs. (13) and (18)."},{"cited_title":"Pasquini and P","cited_arxiv_id":null,"evidence_quote":"Light-front constituent model results for $f_{1\\pi}$ and $h_{1\\pi}^\\perp$ used for comparison, and the one-gluon-exchange treatment of the Boer-Mulders function."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Soft-wall AdS/QCD pion TMD results used as the second comparison for $f_{1\\pi}$."},{"cited_title":"Bacchetta, S","cited_arxiv_id":null,"evidence_quote":"Provides the nonperturbative $g_2$ values $0.09$ and $0.13$ and the $b$-space TMD evolution framework used in Section V."},{"cited_title":"Ahmady, F","cited_arxiv_id":null,"evidence_quote":"Defines the general pion light-front state with spin wave function, used to set up the $L_z=0$ and $|L_z|=1$ amplitudes."},{"cited_title":"Meissner, A","cited_arxiv_id":null,"evidence_quote":"Supplies the $b_\\perp$-space TMD evolution formula with Sudakov form factor used to evolve $f_{1\\pi}$ to $Q=1$ GeV."}],"review_version":1}