{"id":"6dd11cb3-bd5f-439d-a6f3-eb482e9137ee","arxiv_id":"2602.01141","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A phase-space basis ('Wigner function shapelets') is constructed for galaxy images, with explicit formulas for weak-lensing responses and a proposed parity-violation probe.","lead":"Astronomers usually analyze galaxy images either in real space or in Fourier space. This paper builds a new set of analysis tools that describe a galaxy in both spaces at once, using a quantum-mechanics tool called the Wigner function.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"WFS basis is mathematically complete, but it is defined on the unobservable complex amplitude ψ; intensity-only images do not fix the Wigner function, so the claimed estimators lack an observational basis.","rationale":"The paper is a mathematically coherent formalism: the WFS basis is orthogonal and complete, and the expansion coefficients coincide with products of LG-mode coefficients (Eq. 69). I checked the orthogonality via the Moyal identity and found no flaw. The weak-lensing and parity-violating response formulas are plausible but not fully derived in the main text. The load-bearing concern is observability. The Wigner function is defined on the complex field ψ, but detectors record only its modulus. The paper acknowledges this but provides no route from I to W. The ψ vs ψ* ambiguity is a concrete, unavoidable example: the same intensity yields different Wigner functions and hence different WFS coefficients. This means the paper's headline observable W_kℓ is not an estimator of standard astronomical images. This does not invalidate the mathematical contributions, but it does make the claimed applications conditional on phase recovery or interferometric measurements that are not demonstrated. The reader's weakest assumption points to exactly this issue, and my assessment matches. Therefore the reader's CONDITIONAL verdict is appropriate; no change is needed.","tokens_in":39476,"tokens_out":14572,"duration_ms":146431,"concrete_test":"Take a non-real complex field ψ that is a superposition of LG modes, e.g., a single mode with s=1, so ψ*(θ)=Ψ_{j,-s}^{LG}* has the same intensity |ψ|² but opposite orbital angular momentum. (i) Verify |ψ|²=|ψ*|² numerically. (ii) Compute the WFS coefficients c_{j,s,j',s'} from Eq. (69) for both fields, and the corresponding W_kℓ(Q0,Q2) via Eq. (106). (iii) Show that the two sets differ. This demonstrates that the mapping from intensity I to WFS observables is not single-valued, so the formalism cannot be applied to intensity-only images without supplying a phase-recovery prescription.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central mathematical claim—orthogonality/completeness of the WFS basis—is correct and follows from the Moyal/Plancherel identity (Eq. 21) applied to the LG modes. The problem is the gap between the formalism and the data. The Wigner function in Eq. (2) is a bilinear functional of the complex amplitude ψ(θ), but astronomical images supply only I(θ)=|ψ(θ)|², which is the θ-marginal of W (Eq. 15). The paper itself flags this in Sec. II.A ('ψ is often a real positive function as it loses phase information of electromagnetic fields, but we treat ψ as an arbitrary field that can be complex in general') and in Sec. V.D ('reconstruct the image of galaxies within the framework of quantum information theory ... which few have addressed similar algorithms'). No phase-recovery method is provided. This is not a merely practical difficulty: I(θ) does not uniquely determine W(θ,p). For any complex field ψ, its complex conjugate ψ* has the same intensity, yet W_{ψ*}(θ,p)=W_ψ(θ,-p). Under the WFS expansion, the coefficients transform as c'_{j,s,j',s'} = c_{j',s',j,s}^*, which in general yields different W_kℓ(Q0,Q2). Thus the central observable of the paper is not a function of the measured intensity. Without an explicit phase-retrieval step that picks out one of the many fields consistent with I, the weak-lensing and parity-violating estimators of Sec. IV cannot be evaluated from standard CCD images.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces Wigner Function Shapelets (WFS), an orthogonal and complete basis for the Wigner function of a galaxy image, built from the cross-Wigner functions of Laguerre-Gaussian modes. It develops the Hopf-torus projection W_kℓ(Q0,Q2), derives closed forms for the basis, and connects the formalism to BiPoSHs and FPFS. It also derives linear weak-lensing shear and flexion responses, a cosmic-birefringence response, and parity-violating galaxy-shape correlators, and it frames PSF convolution and noise as quantum channels. The stated aim is to make quantum-information phase-space methods available for astronomical image analysis.","tokens_in":39795,"tokens_out":9109,"duration_ms":108360,"significance":"The mathematical core is a legitimate and explicit extension of Simon & Agarwal: orthogonality/completeness follows from the Moyal/Plancherel identity, and the closed forms in Eqs. (95) and (106) are valuable. The Sp(4)/SU(2)/Hopf-torus structure provides an elegant organizing principle for shape modes and lensing selection rules. However, the astrophysical impact is conditional on an unaddressed observational issue: the Wigner function is defined from a complex field, while CCD images deliver only intensity. The paper itself acknowledges this gap in Sec. V.D. If that gap can be closed, the framework could become a genuinely new phase-space morphology tool; as written, its estimators lack a demonstrated data connection.","major_comments":[{"comment":"The Wigner function in Eq. (2) is defined from the complex amplitude ψ(θ), but astronomical images supply I(θ)=|ψ(θ)|², which is only the θ-marginal (Eq. 15). The WFS coefficients in Eq. (69), c_{js,j's'} = ψ_js ψ*_j's', are not functions of I alone: ψ and ψ* give identical intensity but W(θ,p) versus W(θ,−p), producing different W_kℓ and different shear responses. Section V.D concedes that image reconstruction in this framework is 'non-trivial' and that few practical algorithms exist. Without an explicit phase-retrieval or coherence-measurement prescription, the estimators in Secs. IV.A–C (Eqs. 128, 131, 132) have no observational handle. This is load-bearing for the paper's application claims.","section":"Sec. II.A / V.D and Eqs. (2), (15), (69)"},{"comment":"The text asserts as general properties that W is invariant under phase-space reflection and under (θ,p)→(θ,−p), with the justification that W is real. Both claims are false for generic complex ψ: realness does not imply W(θ,p)=W(θ,−p), and a Wigner function is not invariant under spatial reflection unless the underlying field has the corresponding symmetry. The later conditional selection rules in Sec. III.C are derived correctly as symmetry constraints, but they contradict the universal 'invariance' statements in Sec. II.A. The manuscript must correct these property statements, since they bear on the parity-violation discussion.","section":"Sec. II.A, 'Parity' and 'Chirality' bullets"},{"comment":"The cosmic-birefringence response W^I_kℓ − W^S_kℓ = ω(k−ℓ)ρ^S_kℓ is asserted without derivation. The symbol ρ^S_kℓ is undefined (presumably W^S_kℓ), and the physical mechanism is unclear: cosmic birefringence rotates the polarization angle, not the scalar intensity amplitude, so it is not obvious why the scalar Wigner coefficient responds as (k−ℓ)W^S_kℓ. Since this is one of the headline cosmological outputs, it needs a derivation or an explicit reference to a calculation.","section":"Sec. IV.C, Eq. (132)"}],"minor_comments":[{"comment":"Typos: 'inrtoducing' (Sec. III.A), 'Correpspondence' (Contents and Sec. III.F), 'electrocmagnetic' (Sec. II.A), 'Swinger' (Appendix F title should be 'Schwinger'). Section II.A has 'funciton' and other small spelling errors.","section":"General"},{"comment":"Eq. (132) uses ρ^S_kℓ; define it and keep notation consistent with W_kℓ used elsewhere in the same section.","section":"Sec. IV.C"},{"comment":"Figure 1 is not referenced in the main text; the caption should explain the (Q0,Q2) normalization and the color scale. Table I is large; consider moving to an appendix or summarizing the count formula in the text.","section":"Fig. 1 and Table I"},{"comment":"The variable λ is at times called a 'scale' and at others a dimensionless phase-space cell area. The derivation leading to λ≈1.37 is heuristic; clarify that this is an order-of-magnitude estimate, not a precision choice.","section":"Sec. V.B"},{"comment":"The discussion of speckle statistics is difficult to follow and arguably tangential to static galaxy imaging. A brief remark distinguishing coherent speckle from incoherent galaxy-intensity statistics would help.","section":"Sec. II.B.1"}],"recommendation":"major_revision","confidential_remarks":"The mathematical construction is sound and may be publishable as a formal phase-space shapelet framework. The main risk is that the paper is positioned as an astrophysical method while lacking a phase-recovery route; the authors should either provide an explicit observational prescription or clearly re-scope the paper as a mathematical formalism. The incorrect universal parity/chirality claims in Sec. II.A must also be fixed. If the journal accepts purely formal contributions, this is a useful addition after revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is best read as a pure formalism: it takes the known Wigner function of Laguerre–Gaussian modes, extends it to cross-Wigner functions, and shows that these form an orthogonal and complete basis for any square-integrable Wigner function. That claim is solid, following from the Moyal/Plancherel identity and the completeness of LG modes. The author is honest about the lineage, citing Simon & Agarwal and framing the contribution as a transposition to astronomical morphology rather than claiming the diagonal Wigner function as new. The derived shear response, the Hopf-torus observable W_kℓ, and the parity-violating trispectrum estimators are genuinely new in this context, and the quantum-channel treatment of PSF effects is a useful way to think about systematics.\n\nThe soft spots are real and one is load-bearing. The stress-test note is right: the Wigner function needs the complex amplitude ψ, but astronomical images give only I=|ψ|². The paper itself flags this in Sec. II.A and again in Sec. V.D, where it concedes that reconstruction 'few have addressed similar algorithms.' That is not a minor technical aside. For any ψ with the same intensity, ψ* gives a different Wigner function (mirrored in momentum), and the WFS coefficients transform non-trivially. So the central observable W_kℓ is not a function of the measured image. Without an explicit phase-retrieval step—or a restriction to interferometric or polarimetric data where phases are available—the shear and parity estimators in Sec. IV cannot be evaluated on standard CCD images. The paper neither provides that step nor shows that any measurable statistic is insensitive to the ambiguity.\n\nOther concerns are proportionate. Some formulas, like the birefringence response Eq. (132), are asserted without derivation; the flexion response is at least derived in an appendix, so the reader’s complaint is partially addressed. There is no code or simulation, which is fine for a formalism paper but does mean practical confidence is limited. There are also minor typos (e.g., 'Swinger' oscillators, 'Correpspondence') that a referee would catch.\n\nIf the phase-recovery gap is resolved, this could become a genuinely useful framework for weak-lensing and morphology analyses. As it stands, the paper is a mathematically careful construction looking for an observational handle. I would send it to peer review, but with a strong request that the authors either supply a path from intensity data to the Wigner function, or explicitly limit the claims to data where phase information is available. The core math deserves referee time; the astronomical significance does not yet.","headline":"A mathematically sound phase-space shapelet formalism whose central observable is not accessible from intensity-only images—worth refereeing as a formalism paper, but the authors must confront the phase-recovery problem before any astronomical estimator can be taken seriously.","tokens_in":40374,"tokens_out":2582,"would_cite":false,"duration_ms":31547,"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 cross-Wigner functions of Laguerre-Gaussian modes form an orthogonal, complete basis for the Wigner function of any galaxy image, and that the resulting phase-space 'band structure' separates morphology from rotation,","keywords":["Wigner function","shapelets","phase space","Laguerre-Gaussian modes","weak gravitational lensing","cosmic shear","parity violation","Hopf fibration"],"falsifier":"Simulate a galaxy image with a known complex field, apply a small shear γ, compute the Wigner function, expand it in WFS, and check whether the coefficients W_kℓ change exactly as Eq. (128) predicts; then repeat using only intensity data, attempting phase retrieval or forward-model fitting. If the intensity-only route cannot recover W_kℓ, the observable loses its observational handle regardless of the mathematical completeness of the basis.","tokens_in":39265,"feed_emoji":"🔭","tokens_out":5537,"duration_ms":56715,"temperature":0.7,"pith_summary":"This paper extends the familiar shapelet decomposition of galaxy images from separate position and Fourier spaces into a single four-dimensional phase space. It argues that the cross-Wigner functions of Laguerre-Gaussian modes form an orthogonal and complete basis for the Wigner function of a galaxy image, so every square-integrable image can be expanded in them, with coefficients equal to products of the usual Laguerre-Gaussian mode amplitudes. From this basis the paper builds a two-dimensional observable, W_kℓ(Q0,Q2), a kind of 'band structure' labelled by torus winding numbers that separates morphological charge from rotational charge and turns weak-lensing and parity-violating effects into explicit selection rules. If the formalism is correct, it would give galaxy image analysis, cosmic shear measurement, and parity tests on the sky a symmetry-preserving phase-space representation with rigorous completeness.","feed_headline":"Wigner shapelets give galaxy images a complete phase-space basis","feed_subtitle":"The expansion separates shape from rotation and turns weak lensing and parity checks into selection rules.","key_machinery":"The central object is the cross-Wigner function of two Laguerre-Gaussian modes, defined as W_ab(θ,p) = (2πλ)^{-2} ∫ d²ξ ψ_a(θ+ξ/2) ψ_b*(θ−ξ/2) e^{-ip·ξ/λ}; the collection {W^LG_{(j,s),(j',s')}} spans the Hilbert–Schmidt space of phase-space functions. The load-bearing identity is the Moyal–Plancherel relation, which converts mode orthogonality into orthogonality of the cross-Wigner basis. The Hopf spinor parametrisation (u,v) of phase space, with the level-set invariants Q0 (harmonic energy) and Q2 (axial angular momentum), organizes the phase space into Hopf tori; the torus harmonics χ_kℓ(φ_u,φ_v) and the projected observables W_kℓ(Q0,Q2) carry the gauge charges that separate morphology (J=","core_discovery":"The paper establishes that the set of cross-Wigner functions of Laguerre-Gaussian modes, W^LG_{(j,s),(j',s')}, is an orthogonal and complete basis for the Hilbert–Schmidt space of Wigner functions in the four-dimensional phase space of an image. Therefore any quadratically integrable Wigner function can be expanded in this basis, and the expansion coefficients coincide with products of the LG-mode expansion coefficients of the image field. It further introduces Hopf-torus observables W_kℓ(Q0,Q2), which project the Wigner function onto torus harmonics labelled by winding numbers (k,ℓ), with selection rules: statistical isotropy forces k=ℓ, parity forces even k+ℓ, and chirality imposes W_kℓ(Q0","pith_inferences":["Editorial inference: the same cross-Wigner machinery suggests a natural way to combine multi-band or polarised images, since mixed-index cross-Wigner functions could capture colour or polarisation coherence that single-band intensity analysis integrates out; the paper sketches but does not pursue this.","Editorial inference: the W_kℓ band map could serve as a compact morphological summary statistic for large surveys, analogous to power spectra; a testable extension would be to measure W_kℓ from image simulations and check whether it separates morphology classes more cleanly than current shapelet coefficients.","Editorial inference: because real detectors record intensity only, the formalism's observational reach depends on recovering phase information; a concrete practical test would be to reconstruct ψ(θ) from intensity-plus-speckle data and see whether the measured W_kℓ reproduces the predicted shear and flexion responses.","Editorial inference: framing PSF deconvolution as a quantum channel suggests that optimal image reconstruction could be posed as a quantum error-correction problem; the next step would be to construct the Petz–Wiener recovery map from the fidelity kernel and benchmark its shear recovery against existing methods in simulated surveys."],"forward_implications":["Any galaxy image whose Wigner function is square-integrable can be represented exactly in the WFS basis, with coefficients inherited directly from the Laguerre-Gaussian expansion of the image field.","The (Q0,Q2) band map gives a sharp diagnostic: on a statistically isotropic, parity-symmetric ensemble only even-k diagonal modes W_kk survive, so off-diagonal or odd-k power flags systematics, anisotropy, or parity violation.","Cosmic shear and flexion act on the basis through explicit selection rules (Δs=±2 for shear; Δs=±1,±3 for flexion), so lensing responses can be computed locally in phase space rather than through global shape moments.","The parity-odd correlator ⟨W^E W^B*⟩ has zero disconnected contribution for parity-symmetric skies, making the connected four-point function a direct, contamination-free probe of parity-violating physics.","PSF convolution is represented as a quantum channel, and the entanglement-fidelity kernel M identifies which morphology and rotation charges can be recovered, guiding optimal deconvolution in a symmetry-adapted basis."],"fun_headline_variants":["Wigner shapelets give images a complete phase-space basis","Phase-space basis from Wigner shapelets for astronomical images","Wigner shapelets expand images in symplectic phase space","Galaxy images get a full Wigner-function phase-space basis","Wigner shapelets: orthogonal basis for image phase space"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the full complex image field ψ(θ), or equivalently its phase correlations, is accessible from observations, since the Wigner function is defined from ψ(θ) while ordinary detectors record only the intensity |ψ(θ)|².","fun_headline_variants_meta":{"raw":{"variants":["Wigner shapelets give images a complete phase-space basis","Phase-space basis from Wigner shapelets for astronomical images","Wigner shapelets expand images in symplectic phase space","Galaxy images get a full Wigner-function phase-space basis","Wigner shapelets: orthogonal basis for image phase space"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000312,"raw_usage":{"total_tokens":1702,"prompt_tokens":928,"completion_tokens":774,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":672,"completion_tokens_details":{"reasoning_tokens":705}},"tokens_in":672,"tokens_out":774,"duration_ms":8372,"temperature":1.0,"reasoning_tokens":705,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T05:44:47.868399+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate a galaxy image with a known complex field, apply a small shear γ, compute the Wigner function, expand it in WFS, and check whether the coefficients W_kℓ change exactly as Eq. (128) predicts; then repeat using only intensity data, attempting phase retrieval or forward-model fitting. If the intensity-only route cannot recover W_kℓ, the observable loses its observational handle regardless of the mathematical completeness of the basis.","supporting_citations":[],"review_version":1}