{"id":"fd1a2d74-dcd6-4f5c-aa9e-e7fda7d9df9a","arxiv_id":"2601.03326","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Higher-order moment tensors and Gaussian-times-polynomial expansions give rotation-invariant shape descriptors that can be made arbitrarily detailed, extending PCA.","lead":"This paper proposes describing shapes with higher-order moment tensors and polynomial-times-Gaussian expansions, then using graph-based tensor contractions as rotation-invariant features. It extends PCA's ellipsoid approximation to arbitrary detail, aiming at molecular shape descriptors and rotation-invariant object recognition without expensive rotation search.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Higher-order graph invariants are only necessary conditions (Sec. III.B); without completeness or experiments, the claim of describing 'shape modulo rotation' is unsupported.","rationale":"The reader's weakest assumption correctly identifies the load-bearing issue: the invariants for orders ≥3 are only necessary conditions, and the paper provides no evidence they are sufficiently discriminative. My stress-test agrees with that assessment. The paper has real strengths: the polynomial-times-Gaussian representation using Hermite functions is mathematically sound, and the order-2 covariance invariants are well grounded. The paper also honestly flags the completeness problem in Section III.B and in the conclusion lists it as an open question. However, the title and abstract claim 'shape modulo rotation' description, and the applications (molecular descriptors, recognition, similarity metrics) rely on the invariants acting as a faithful signature. Without either a completeness proof or a practical evaluation on a recognition task, this central claim is not supported. The concern is not an internal inconsistency or an attack on the author; it is an unvalidated assumption that is explicitly load-bearing. A concrete computational test on the simplest case (d=2, r=3) can settle whether the proposed graph invariants even separate orbits in that setting, which would either refute the general claim or narrow the open problem. Since the reader's verdict is already CONDITIONAL and my analysis does not move it, I recommend UNCHANGED.","tokens_in":5809,"tokens_out":11144,"duration_ms":109520,"concrete_test":"For d=2, r=3 (symmetric order-3 tensors), compute all graph invariants generated by the contractions in Sec. III.B up to, say, 4 vertices. Independently compute the complete rotation-invariant set from the Fourier decomposition p(cos θ, sin θ) = Σ_{k=-3,-1,1,3} A_k e^{ikθ}, with invariants |A_3|, |A_1|, and arg(A_3) − 3 arg(A_1). Then solve for two distinct tensors that have identical graph-invariant values but different complete invariant values. If such a pair exists, the proposed descriptor is not injective modulo rotation even in the simplest case; if it does not, the concern shifts to higher dimensions/orders, but the paper's stated incompleteness must be addressed for the general claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the proposed invariants describe shapes modulo rotation. For orders ≥3, however, the paper itself states in Section III.B: 'while agreement of such invariants is necessary conditions, to be certain that [p]∼[q] we would also need sufficient condition: a complete set of invariants ... for higher orders it seems a difficult open problem.' The proposed similarity metric is a distance between finite vectors of graph-based invariants, and any finite set of necessary-only invariants can map two non-congruent shapes to the same feature vector, producing false matches. No experiment or argument shows that the invariants are sufficiently discriminative in practice, despite the applications (molecular shape, 2D/3D recognition) depending on it. Concretely, in d=2 for order-3 tensors, the tensor space has 4 dimensions and SO(2) orbits have dimension 1, so 3 independent invariants are needed. The graph invariants shown in Fig. 3/4 for r=3 are quadratic (Frobenius and a self-loop contraction), and a degree-3 invariant is impossible because the number of tensor indices is odd; the invariant ring is generated in higher degree, and it is not shown that the finite set actually separates orbits. The paper's honest admission plus the absence of discriminative validation makes 'shape modulo rotation' an overstatement relative to what is established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes to extend PCA's covariance-matrix shape descriptor to higher-order central moments and to polynomial-times-Gaussian expansions, and to use graph-based tensor contractions as rotation-invariant features. These invariants are gathered into a feature vector intended to describe shapes modulo rotation and to define an inexpensive shape-similarity metric. The paper illustrates the polynomial representation on MNIST digits but provides no empirical evaluation, and it explicitly acknowledges that for order ≥3 a complete set of rotation invariants is an open problem.","tokens_in":6054,"tokens_out":8543,"duration_ms":74776,"significance":"The underlying mathematics is largely sound: traces of powers determine eigenvalues for symmetric matrices, and graph contractions indeed yield rotation invariants. The idea of extending PCA beyond order-2 moments is natural and could be of practical use in chemoinformatics and 2D/3D recognition. However, the central claim—that the proposed features describe shape modulo rotation—is not established for order ≥3, because the graph invariants are only necessary conditions, not sufficient ones. The paper is honest about this limitation, but the abstract and conclusion overstate the result. The MNIST figure demonstrates representation fidelity, not invariance discriminative power.","major_comments":[{"comment":"The paper admits that for order ≥3 a complete set of invariants is a difficult open problem and that agreement of graph invariants is only a necessary condition. Thus the proposed finite invariant vectors cannot be guaranteed to separate rotation orbits, and two non-congruent shapes may share the same descriptor. In d=2, the space of order-3 tensors has dimension 4 while SO(2) has dimension 1, so at least 3 independent invariants are needed; the paper provides no such set nor any argument that its chosen invariants are sufficiently discriminative. This gap directly undermines the 'shape modulo rotation' claim in the title and abstract.","section":"Section III.B"},{"comment":"The proposed shape similarity metric is not actually specified. The text says only 'some distance between two such vectors' without defining a concrete distance, the set of features to use, or how to handle the redundancy and weighting of invariants. Without a specific metric, claims that the approach enables 'inexpensive comparison' and avoids costly rotation optimization cannot be tested or reproduced.","section":"Section III.C"},{"comment":"No experiments are reported. Despite the stated applications (molecular shape, 2D/3D object recognition), there is no empirical evidence that the higher-order invariants separate non-congruent shapes or that the feature vector is robust to noise and deformations. Given that the theoretical completeness is open, a synthetic benchmark measuring orbit separation or false match rates for d=2 and d=3 tensors is essential to support the practical significance of the method.","section":"General"}],"minor_comments":[{"comment":"The Hermite polynomial list contains a typo: the degree-5 polynomial should be 32x^5 - 160x^3 + 120x, not '32x^2 - 160x^3 + 120x'. Also, the notation for Hermite polynomials is inconsistent (h_i vs f_j).","section":"Section II.E"},{"comment":"The notation [p] is used for the covariance matrix but also for tensors in later sections; this should be clarified to avoid confusion.","section":"Section II.B"},{"comment":"The formula for p_ab is incomplete, ending with an ellipsis. The abstract would benefit from a cleaner presentation of the central equations.","section":"Abstract"},{"comment":"Typo: 'mane of them' should be 'many of them'.","section":"Section III.C"},{"comment":"Figures 3 and 4 are referenced but not described in sufficient detail in the text; the graph construction should be defined with explicit equations so that the invariants are unambiguous.","section":"Section III.B"}],"recommendation":"major_revision","confidential_remarks":"The paper is a proposal with a mathematically correct core, but the gap between its claims and what is demonstrated is substantial. If the authors can supply a completeness result for small d (e.g., d=2 or d=3, where classical invariant theory provides complete sets) or provide strong empirical evidence of discriminative power, the paper could become acceptable. Otherwise, the claims should be substantially weakened to reflect the necessary-only nature of the invariants. The paper's reliance on the author's own preprints for key constructions is also a concern for novelty assessment."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Jarek's paper proposes extending PCA's covariance ellipsoid to higher-order moment tensors and to polynomial-times-Gaussian expansions, with rotation invariants from graph contractions. The genuinely useful part is the framing: a single recipe that covers central moments, decodable polynomial approximations, and a similarity metric based on invariant vectors, plus the honest discussion of what is known. The MNIST reconstruction figures give a nice intuition for how expressive the polynomial representation is.\n\nThe mathematics that is stated is mostly correct. Traces of powers of the covariance matrix do determine the matrix up to rotation, and the diagrammatic contractions are rotation invariant. The paper also clearly flags that a complete set of invariants for order ≥3 is open and that the independence problem is hard. That is more candid than most.\n\nThe soft spots are the ones the paper itself names. The central claim of describing shape modulo rotation is not established for order ≥3. A finite vector of necessary conditions can easily collide for non-congruent shapes; no experiment or argument shows the invariants are discriminative in practice. The stress-test point about degree-3 invariants in d=2 is a concrete instance: the graph invariants shown are quadratic, and it is not shown that they separate SO(2) orbits. Also, the diagrammatic machinery is from the author's earlier [7], and classical moment invariants (Hu, Zernike) are not cited. So the novelty is mostly the application and the polynomial-Gaussian decoding idea, not the invariant construction.\n\nThat said, the paper is explicitly an initial proposal and says so. The open-problem section is honest, and the proposed framework is coherent enough that someone could test it. It does not overclaim in the conclusion; it lists exactly what remains.\n\nWho is this for? Someone working on rotation-invariant shape descriptors who wants a compact summary of this particular approach and its open problems. It is not a validation paper. If it came to a journal, I would send it to a careful referee, mainly because the proposal is concrete and the limitations are stated well enough that a referee could ask for targeted experiments. But it should not be published as-is; it needs either a proof of separation for a useful invariant set or empirical demonstration on a recognition/similarity task.\n\nI would cite it if I were writing about higher-order moment invariants, because it cleanly lays out the recipe and the open question, and I would bring it up in a reading group as an example of a paper that is honest about its own limits.","headline":"Honest proposal, but the headline claim outruns what is actually shown: the invariants are necessary, not sufficient, and no experiment demonstrates they separate shapes.","tokens_in":6556,"tokens_out":1698,"would_cite":true,"duration_ms":16759,"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":"Rotation-invariant shape features go beyond PCA's ellipsoid.","keywords":["rotation invariants","shape descriptors","central moments","tensors","polynomial expansion","Hermite polynomials","shape similarity","PCA"],"falsifier":"Construct two non-congruent 3D point sets with the same covariance ellipsoid but different third-order structure, then compute all graph-based invariants up to a chosen order; if the two sets produce identical invariant vectors, the descriptor fails to separate shape modulo rotation.","tokens_in":5624,"feed_emoji":"🔄","tokens_out":4418,"duration_ms":38911,"temperature":0.7,"pith_summary":"This paper proposes extending PCA's covariance-ellipsoid shape descriptor to higher-order central moment tensors (order 3 and up) and to polynomial-times-Gaussian expansions. The aim is a rotation-invariant shape description that can be made as detailed as needed and that can be decoded back into a shape. If the approach works, two shapes that differ only by rotation can be recognized and compared by computing a fixed vector of invariants, without costly optimization over rotations. This matters for molecular shape description in chemoinformatics, 2D image and 3D scan recognition, and cheap shape-similarity metrics where rotation should be ignored.","feed_headline":"Rotation-invariant shape features go beyond PCA's ellipsoid","feed_subtitle":"Higher-order moment tensors and polynomial-times-Gaussian expansions make shape description modulo rotation as detailed as needed.","key_machinery":"The load-bearing object is the symmetric tensor of central moments of order r, defined by averaging products of r centered coordinates, together with the graph-based invariants obtained by contracting such tensors over paired indices; each contracted edge remains invariant under a common orthogonal rotation because O^T O = I. For decodable and complete representations, the paper uses an orthonormal basis of Hermite polynomials times a Gaussian, whose coefficients are expectations of basis functions over the shape.","core_discovery":"The central claim is that extending PCA to order-r central moments gives rotation-invariant features that describe a shape modulo rotation with as much detail as needed, and that polynomial-times-Gaussian expansions make the description decodable and complete in the limit as the polynomial degree grows. The paper shows how to build these tensors from densities, volumes, or weighted point sets, and how to generate graph-based invariants—such as sums of squared entries of the order-3 tensor—that are unchanged by rotation because each contracted edge pairs matching indices. It notes that these invariants are necessary conditions, and that a complete set for order 3 and higher is an open problem","pith_inferences":["Because completeness is open for order 3 and higher, a practical risk is that two non-congruent shapes share the same invariant set; this could be probed empirically by testing discrimination on datasets of rotated versus non-congruent 2D and 3D shapes.","A natural extension is to train machine-learning models directly on these invariant vectors, using them as fixed-dimensional features alongside existing descriptors.","In d=2 and d=3, where the rotation group has only 1 or 3 dimensions, even redundant invariant sets may be discriminative, so a direct comparison against alignment-based methods could test their practical sufficiency."],"forward_implications":["Shape similarity modulo rotation can be computed as a distance between invariant feature vectors, with no iterative rotation search.","The polynomial-times-Gaussian representation gives a continuous, reconstruction-capable shape descriptor whose accuracy increases with the maximum polynomial degree.","The same invariants apply to densities, volumes, weighted point sets (e.g. atoms), and grayscale images.","Variability of a shape, such as molecular dynamics, can be encoded as a distribution of invariant vectors rather than a single vector."],"fun_headline_variants":["Higher-order moments give detailed rotation-invariant shape descriptors","Beyond PCA: tensor moments handle complex shapes modulo rotation","Rotation-invariant shape detail from higher-order central moments","Polynomial-Gaussian expansions refine rotation-invariant shape features","Shape description modulo rotation gets sharper with tensor moments"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The argument depends on a finite set of graph-based invariants for order-3 and higher tensors being discriminative enough—only necessary conditions are known, and the paper itself says a complete set is a difficult open problem; if two non-congruent shapes share all tested invariants, the descriptor no longer describes shape modulo rotation.","fun_headline_variants_meta":{"raw":{"variants":["Higher-order moments give detailed rotation-invariant shape descriptors","Beyond PCA: tensor moments handle complex shapes modulo rotation","Rotation-invariant shape detail from higher-order central moments","Polynomial-Gaussian expansions refine rotation-invariant shape features","Shape description modulo rotation gets sharper with tensor moments"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000398,"raw_usage":{"total_tokens":1894,"prompt_tokens":692,"completion_tokens":1202,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":436,"completion_tokens_details":{"reasoning_tokens":1126}},"tokens_in":436,"tokens_out":1202,"duration_ms":9445,"temperature":1.0,"reasoning_tokens":1126,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T12:21:19.976227+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Construct two non-congruent 3D point sets with the same covariance ellipsoid but different third-order structure, then compute all graph-based invariants up to a chosen order; if the two sets produce identical invariant vectors, the descriptor fails to separate shape modulo rotation.","supporting_citations":[],"review_version":1}