{"id":"40169fde-dee6-412a-9ed0-f11a26e368fc","arxiv_id":"2506.00510","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Quark Wigner maps in boost-invariant 3D impact-parameter space are centrally peaked, symmetric, and show orbital-like nodal rings.","lead":"Scientists modeled a quark as a 'dressed' quark with a cloud of gluons, then computed where the quark is likely to be found in a three-dimensional, boost-invariant position space. The maps are centered on the quark's average position, symmetric, and show ring-shaped empty regions the authors compare to atomic orbitals.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The orbital-like zero-density rings are the central claim, but they depend on an unverified Fourier measure and an untested k_perp cutoff; if either is corrected, the nodal pattern may shift or disappear.","rationale":"The reader's conditional verdict is appropriate. My main concern overlaps with the reader's weakest assumption about the k_perp cutoff, but I would put at least as much weight on the unverified Jacobian/Fourier-measure factors connecting Eq. (12) to Eqs. (17)-(24). If the measure is wrong, every plotted distribution is affected and the zero locations in particular are not trustworthy. I do not see grounds for rejection: the model construction is a straightforward extension of the authors' earlier published GTMD work [32], and the missing items are checkable rather than irreparable. The absence of parameter values and convergence tests, however, genuinely prevents the central 'atomic-orbital' conclusion from being assessed as stated. Therefore the verdict should remain CONDITIONAL, with the condition being a corrected, fully specified numerical computation plus a k_max sensitivity study.","tokens_in":8688,"tokens_out":7866,"duration_ms":77398,"concrete_test":"Recompute the unpolarized rho_UU map (Eq. 17) with an independently re-derived Fourier measure: track D_perp = Delta_perp/(1-xi^2) through Eq. (10) and verify whether the (1-xi^2)^(-3/2) factor in Eq. (17) is correct or should instead be (1-xi^2)^(-2). Separately, repeat the k_perp integration for k_max = 0.5, 1, 2, and 5 GeV (or with a smooth UV regulator) at x=0.3 with explicitly stated m and g. If the zero-density rings persist at the same (b_perp, sigma) locations under both changes, the orbital claim survives; if they shift or vanish, Section 6 does not support the claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central physical claim is that the Wigner distributions in (b_perp, sigma) space contain symmetrically placed zero-density regions resembling atomic orbitals. The only evidence is Figs. 1-3, generated from Eqs. (17)-(24). Two inputs control the nodal pattern and neither is adequately supported. First, the Fourier measure: Eq. (10) uses d^2D_perp/(2pi)^2 with D_perp = Delta_perp/(1-xi^2), so substituting D_perp gives d^2Delta_perp/[(2pi)^2(1-xi^2)^2]. Eq. (12), however, writes d^2Delta_perp/(1-xi^2), and Eq. (17) introduces an unexplained factor (1-xi^2)^(-3/2). No derivation of this factor is provided, yet it changes the phase weight that determines where the oscillatory integrals have zeros. Second, the ultraviolet behavior: Section 6 states only x=0.3 and integrates k_perp over (0,0.5) GeV. The GTMD denominators from [32] behave like [m^2 - (m^2+q^2)/x - q^2/(1-x)]^2, so the k_perp integral is not manifestly convergent; a hard cutoff at 0.5 GeV can impose artificial oscillations. The quark mass, coupling, and grid parameters are not stated, so the maps are not reproducible. Because the zero-density 'orbital' pattern is the load-bearing conclusion and its stability under the correct measure and a converged k_perp integral is not demonstrated, the central claim is currently unverified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes quark Wigner distributions in a boost-invariant three-dimensional position space (transverse impact parameter b_perp and longitudinal coordinate sigma) within the dressed quark model. The distributions are obtained by Fourier transforming GTMDs with respect to the skewness xi and the transverse momentum transfer Delta_perp, using GTMD expressions from the authors' earlier paper [32]. Numerical plots are presented for various target and quark polarizations, showing central concentration, reflection symmetry, and zero-density regions that the authors compare to atomic orbitals. The claimed novelty is the extension of Wigner distribution studies to a frame-independent 3D coordinate space.","tokens_in":9055,"tokens_out":7411,"duration_ms":62768,"significance":"If the derivation and numerics were sound, the paper would provide a useful model calculation of quark phase-space structure in a boost-invariant 3D setting, extending prior 2D impact-parameter studies. The apparent 'atomic orbital' zero-density pattern, if genuine, would be an interesting prediction of spatial quantization. However, the central claim is not currently established because of an inconsistent Fourier measure and an untested transverse-momentum cutoff. The paper also benefits from a clear connection to the existing GTMD literature, and the qualitative symmetry and central-peak features are robust expectations for such model distributions.","major_comments":[{"comment":"The Fourier measure is inconsistent between Eq. (10), Eq. (12), and Eqs. (17)-(24). Since D_perp = Delta_perp/(1-xi^2), the measure in Eq. (10) is d^2 Delta_perp / [(2 pi)^2 (1-xi^2)^2], but Eq. (12) writes d^2 Delta_perp / (1-xi^2), and Eqs. (17)-(24) contain a factor (1-xi^2)^(-3/2) without derivation. This factor multiplies the phase of the oscillatory integrals and therefore controls the positions of the reported zero-density regions, so the central claim is not supported by a self-consistent derivation; please derive the Jacobian and any spinor-normalization factors explicitly.","section":"Section 5, Eqs. (10)-(24)"},{"comment":"The transverse momentum integral is truncated to |k_perp| in (0, 0.5) GeV with no sensitivity check. The GTMD denominators from Ref. [32] are not manifestly convergent in k_perp, so this hard cutoff can produce artificial oscillations and zero-density nodes. To establish that the orbital-like pattern is a property of the model rather than of the cutoff, show convergence as the cutoff increases and provide a physical justification for the chosen scale.","section":"Section 6"},{"comment":"The numerical results are not reproducible because the quark mass m and coupling g are never stated. Since the nodal structure of the Wigner maps may depend on these parameters, provide the values used in Figs. 1-3 and study the parameter dependence of the claimed zero-density regions.","section":"Section 6"},{"comment":"The GTMD inputs are taken from Ref. [32] without being reproduced or summarized in this manuscript. Because the Wigner distributions are obtained as Fourier transforms of those GTMDs, the paper is not self-contained and the reader cannot verify the transformation without consulting the prior work; please include the relevant GTMD formulas (or an appendix) and specify the numerical integration grid and domains.","section":"Sections 4-6"}],"minor_comments":[{"comment":"The paper refers to 'quark probability density' and 'zero-density regions,' but the Wigner distribution is a quasi-probability that can be negative; the atomic orbital analogy should be qualified accordingly.","section":"Section 6"},{"comment":"Several formulas have ambiguous typesetting of fractions, e.g., '-i/m2(1-xi^2)^(3/2)'; please ensure all expressions are typeset unambiguously.","section":"Equations (17)-(24)"},{"comment":"The integration range is stated as '(0 -> 0.5) GeV'; please specify that this is |k_perp| and give the full integration domains for xi and Delta_perp.","section":"Section 6"},{"comment":"The symbol rho is used for both the Wigner distribution and the correlator; please use distinct notation or define the conventions clearly.","section":"Section 4"},{"comment":"The definition of sigma = (1/2) b^- P^+ appears without derivation; it would help to explain why this variable is boost-invariant.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a short numerical application of the authors' own prior GTMD calculation; the primary new element is the 3D Wigner map in (b_perp, sigma). The central claim is interesting but is currently undermined by the inconsistent Fourier measure and the untested cutoff. If the authors can fix the derivation and provide sensitivity and reproducibility details, the paper could become suitable; the novelty is modest but acceptable for a model calculation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent but incremental follow-up to the authors' own GTMD calculations. The new piece is the joint (σ, b⊥) Wigner maps in a boost-invariant 3D space, which I haven't seen elsewhere. The qualitative features – concentration near the center and inversion symmetry – are solid and likely robust. The paper's headline claim, however, is the appearance of zero-density regions that resemble atomic orbitals. That claim rests on two wobbly supports: an unexplained factor in the Fourier measure and an arbitrary transverse-momentum cutoff. Neither is adequately supported, so the orbital analogy is currently unverified.\n\nWhat the paper does well: it systematically covers all polarization configurations for both target and struck quark, and it clearly explains the kinematic setup. The plots are readable and the central-peak behavior is consistent with physics expectations. The approach is a natural extension of [31–33], and the authors are honest that the GTMD inputs come from their earlier work.\n\nThe soft spots are real. First, the measure. Eq. (10) defines the Fourier transform over D⊥ = Δ⊥/(1−ξ²). Substituting D⊥ in terms of Δ⊥ should give a Jacobian (1−ξ²)^{-2}. Eq. (12) writes the Δ⊥ integral with a single (1−ξ²)^{-1} and no (2π)^{-2}, and Eqs. (17)–(24) use (1−ξ²)^{-3/2}. No derivation is given for these powers. The difference matters because the zeros of the oscillatory integrals are controlled by the phase weight; a wrong power can shift or wipe out the zeros. This is the load-bearing issue for the orbital claim.\n\nSecond, the numerics. The quark transverse momentum is integrated over (0, 0.5) GeV with no physical justification and no sensitivity test. The GTMD denominators from [32] don't obviously make the k⊥ integral converge fast, so the hard cutoff can impose artificial oscillations. The quark mass and coupling are not stated, so the plots are not reproducible without guessing. The orbital analogy is also purely qualitative – there's no comparison to actual orbital nodal patterns or any quantitative metric.\n\nThat said, the paper isn't a waste of time. The central-peak and symmetry observations would survive even if the measure factor is corrected. But the atomic-orbital conclusion, the part that makes the paper interesting, depends on exactly those unverified ingredients.\n\nWho's this for? People working on light-front models and Wigner distributions of hadrons. It won't change experimental or phenomenological practice. Still, it's a reasonable model study that could be worth publishing after revision. I'd send it to a referee, not desk-reject. The referee should ask for a clean derivation of the measure factors, a cutoff sensitivity study, and explicit parameter values.","headline":"The 3D Wigner maps are a nice new visualization, but the atomic-orbital claim rests on an unexplained Fourier measure factor and an untested momentum cutoff, so it's not yet trustworthy.","tokens_in":9571,"tokens_out":8016,"would_cite":false,"duration_ms":66107,"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":"A quark's 3D position map shows atomic-orbital-like nodes.","keywords":["Wigner distribution","dressed quark model","GTMDs","light-front coordinates","boost-invariant impact parameter","skewness","spatial quantization","quark phase-space density"],"falsifier":"Recompute the Wigner distributions with the $k_\\perp$ integral extended well beyond 0.5 GeV, or with a smooth cutoff, while scanning the quark mass and coupling; if the zero-density rings move, disappear, or lose their symmetry, the claimed spatial quantization is a numerical artifact rather than a feature of the dressed quark model.","tokens_in":8494,"feed_emoji":"⚛️","tokens_out":8687,"duration_ms":75634,"temperature":0.7,"pith_summary":"The paper tries to establish that a quark in a dressed quark model has a structured spatial probability distribution when viewed in a frame-independent three-dimensional position space built from the transverse impact parameter $\\mathbf{b}_\\perp$ and the boost-invariant longitudinal impact parameter $\\sigma$. The computed Wigner distributions peak at the center, fall off in both directions, and contain symmetrically placed regions where the quark probability density is exactly zero, a pattern the authors compare to atomic orbitals. The result matters because it suggests that a quark's phase-space image inside a target is not a featureless smear but carries a spatial quantization structure, and it offers a concrete model-based target for future hadron-structure calculations. The qualitative shape is stable across polarization choices, although polarization sharpens the peaks and troughs.","feed_headline":"A quark's 3D position map shows atomic-orbital-like nodes","feed_subtitle":"Model maps quark position in frame-independent 3D space; nodes resemble atomic orbitals.","key_machinery":"The central object is the three-dimensional Wigner distribution built from two Fourier transforms: one taking the skewness $\\xi$ to the boost-invariant longitudinal impact parameter $\\sigma=\\tfrac{1}{2}b^-P^+$, and one taking the transverse momentum transfer $\\Delta_\\perp/(1-\\xi^2)$ to the transverse impact parameter $\\mathbf{b}_\\perp$. The integrand is the quark\\textendash quark correlator, expressed through generalized transverse momentum\\textendash dependent parton distributions (GTMDs) $F_{1,i}$, $G_{1,i}$, $H_{1,j}$ at nonzero skewness; in the dressed quark model these GTMDs descend from the two-particle light-front wavefunction of the target state. The named objects are GTMDs, the generalization of parton distributions that keep both transverse momentum and transverse position information. This double-Fourier machinery converts the model's momentum-space wavefunctions into a spatial map that can be read as a probability density for finding the quark at a given $\\sigma$ and $\\mathbf{b}_\\perp$.","core_discovery":"Within the dressed quark model, the authors report that the quark Wigner distributions $\\rho_{XY}(x,\\sigma,\\mathbf{b}_\\perp)$ in boost-invariant three-dimensional space, after integrating over the quark transverse momentum $k_\\perp$, are concentrated near $\\mathbf{b}_\\perp=0$ and $\\sigma=0$, symmetric under reflection in both coordinates, and interspersed with regions of zero density that are symmetrically distributed around the origin. They interpret these nodal regions as resembling atomic orbitals, where occupation probability is concentrated in lobes separated by nodes, and take the pattern as evidence of spatial quantization around the target center. Their plots are made at fixed $x=0.3$ with the $k_\\perp$ integral truncated to the interval $(0, 0.5)$ GeV, and the same qualitative features appear for unpolarized, longitudinally polarized, and transversely polarized targets and struck quarks.","pith_inferences":["The zero-density rings are demonstrated only under a transverse momentum cutoff at 0.5 GeV with no sensitivity study; a natural next step is to check whether the nodal pattern survives as the cutoff is raised or replaced by a smooth regulator.","If the nodes are physical, they should also appear as zeros in the DVCS amplitude in the same $(\\sigma,\\mathbf{b}_\\perp)$ conjugate space, giving an indirect experimental avenue to test the model.","The analogy to atomic orbitals suggests that each node may trace back to a zero of the two-particle light-front wavefunction; locating those zeros explicitly would turn the resemblance into a quantitative prediction."],"forward_implications":["If the reported pattern is real, the quark at fixed $x=0.3$ is most likely found at the center of the target, with equal probability on either side in both the longitudinal and transverse directions.","The symmetric zero-density regions imply that the spatial density has genuine nodes, not merely a smooth exponential falloff, which is a sharp qualitative signature of the dressed quark model.","Because the overall shape barely changes when the target or quark polarization is changed, the peak-and-node structure is a property of the unpolarized Wigner distribution, not an accident of one spin configuration.","The same construction can be applied to other targets, such as hadrons and mesons, and to gluons, using the corresponding light-front wavefunctions, as the authors propose at the end of the paper."],"supporting_citations":[{"why":"Supplies the analytical nonzero-skewness GTMD expressions that are the input to all Wigner distribution formulas.","marker":"[32]"},{"why":"Establishes the dressed quark model framework and earlier zero-skewness Wigner distributions for this target.","marker":"[11]"},{"why":"Introduces the boost-invariant longitudinal position space and Wigner distributions at nonzero skewness used in Eq. (11).","marker":"[31]"},{"why":"Defines the three-dimensional invariant coordinate space $(\\mathbf{b}_\\perp,\\sigma)$ from DVCS, which this paper adopts.","marker":"[26]"},{"why":"Introduces the $\\sigma$ coordinate conjugate to skewness in DVCS, motivating the longitudinal position space.","marker":"[25]"},{"why":"Provides the GTMD parametrization of the quark\\textendash quark correlator used in Eqs. (7)\\textendash(9).","marker":"[30]"},{"why":"Connects Wigner distributions to GTMDs through the Fourier transform used in Eq. (10).","marker":"[10]"}],"fun_headline_variants":["Frame-independent quark Wigner map shows orbital nodes","Quark Wigner map reveals atomic-orbital-like nodes","3D quark position map mirrors atomic orbital lobes","Quark Wigner density in 3D resembles orbital nodes","Frame-independent quark map exhibits orbital-like nodes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The pattern of zero-density regions rests on a truncated transverse momentum integral ($k_\\perp$ from 0 to 0.5 GeV) together with unreported quark mass and coupling parameters; if the truncation or parameter choice distorts the integrand, the orbital-like nodes could be artifacts rather than quark structure.","fun_headline_variants_meta":{"raw":{"variants":["Frame-independent quark Wigner map shows orbital nodes","Quark Wigner map reveals atomic-orbital-like nodes","3D quark position map mirrors atomic orbital lobes","Quark Wigner density in 3D resembles orbital nodes","Frame-independent quark map exhibits orbital-like nodes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000186,"raw_usage":{"total_tokens":1255,"prompt_tokens":802,"completion_tokens":453,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":418,"completion_tokens_details":{"reasoning_tokens":377}},"tokens_in":418,"tokens_out":453,"duration_ms":4557,"temperature":1.0,"reasoning_tokens":377,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:03:31.561950+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the Wigner distributions with the $k_\\perp$ integral extended well beyond 0.5 GeV, or with a smooth cutoff, while scanning the quark mass and coupling; if the zero-density rings move, disappear, or lose their symmetry, the claimed spatial quantization is a numerical artifact rather than a feature of the dressed quark model.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the analytical nonzero-skewness GTMD expressions that are the input to all Wigner distribution formulas."},{"cited_title":"Mukherjee, S","cited_arxiv_id":null,"evidence_quote":"Establishes the dressed quark model framework and earlier zero-skewness Wigner distributions for this target."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the boost-invariant longitudinal position space and Wigner distributions at nonzero skewness used in Eq. (11)."},{"cited_title":"Brodsky, D","cited_arxiv_id":null,"evidence_quote":"Defines the three-dimensional invariant coordinate space $(\\mathbf{b}_\\perp,\\sigma)$ from DVCS, which this paper adopts."},{"cited_title":"Brodsky, D","cited_arxiv_id":null,"evidence_quote":"Introduces the $\\sigma$ coordinate conjugate to skewness in DVCS, motivating the longitudinal position space."},{"cited_title":"Meißner, A","cited_arxiv_id":null,"evidence_quote":"Provides the GTMD parametrization of the quark\\textendash quark correlator used in Eqs. (7)\\textendash(9)."},{"cited_title":"Lorce, B","cited_arxiv_id":null,"evidence_quote":"Connects Wigner distributions to GTMDs through the Fourier transform used in Eq. (10)."}],"review_version":1}