{"id":"f9fa2405-ce22-4f30-ab61-5f1ba2c70f56","arxiv_id":"2411.18232","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Random polarization fields have a jacobian whose sign domains form patches that map to facets on the Poincare sphere, with folds, cusps, heavy-tailed statistics, and percolation behavior.","lead":"Random light patterns, like laser speckle, have a polarization state at every point; this paper treats the mapping from the plane to the sphere of polarization states as a sheet of paper that folds, creases, and cusps. It shows the folds form patches with predictable statistics and connects them to skyrmionic topological quasiparticles.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (14) gives P(|ρ|>aK2)=(1+a)^{-2}, so the paper's stated area fractions 75% (a=1) and 36% (a=4) are incorrect; these numbers match the share of ⟨|ρ|⟩ contributed by those regions instead.","rationale":"The reader's weakest assumption concerns the unverified origami-manifold conditions in Section 5. I agree that Section 5 is speculative, but the paper itself frames it as a brief outline requiring more work, so it is not the most load-bearing part of the central claim. The most load-bearing problem lies in the quantitative core: Eq. (14) is not derived, and even granting it, the paper misreads its own tail probabilities. The complementary CDF is (1+a)^{-2}, so the area fractions quoted in Section 4 are wrong by large factors; the printed numbers instead equal the shares of ⟨|ρ|⟩ contributed by those regions. This is an internal inconsistency that can be checked immediately, and it undermines the precision of the statistical conclusions about superoscillations. The geometric interpretation of patches as sign-domains of ρ, with fold lines and cusps from Whitney theory, is sound, and the local cusp model is helpful. A smaller technical error worth noting but not selecting as the headline: Section 3.3 calls cusp points 'rank zero', whereas the derivative of the model map (10) has rank one at the cusp; this is a minor correction. Given these issues, the CONDITIONAL verdict remains appropriate, with the requested revisions including a corrected statement of the tail probabilities and a derivation or reference for Eq. (14).","tokens_in":16306,"tokens_out":12515,"duration_ms":113968,"concrete_test":"Recompute the complementary CDF from Eq. (14): P(|ρ|>aK2)=(1+a)^{-2}, giving 25% for a=1 and 4% for a=4, and compare with the manuscript's 75% and 36%. Also compute the contribution share E[|ρ| 1_{|ρ|>aK2}]/E|ρ| = 2∫_a^∞ u(1+u)^{-3}du, which reproduces 75% and 36%; this check will show whether the text should say 'area fraction' or 'share of ⟨|ρ|⟩', and any downstream superoscillation claims must be revised accordingly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4's only analytic statistical result, Eq. (14), is applied incorrectly in the sentence: 'the fraction of the area where |ρ| > K2 is 75% of the plane, and where |ρ| > 4K2 is 36% of the plane.' From Eq. (14), P(|ρ|>aK2) = 2∫_{aK2}^∞ K2^{-1}(1+|ρ|/K2)^{-3} dρ = (1+a)^{-2}. For a=1 this is 25%, and for a=4 this is 4%, not 75% and 36%. The printed values equal the fractions of the mean ⟨|ρ|⟩ contributed by those regions: 2∫_a^∞ u(1+u)^{-3}du gives 0.75 and 0.36. Thus the manuscript conflates 'fraction of area' with 'share of the mean topological charge density'. This is a concrete internal inconsistency, independent of the missing derivation of Eq. (14). The downstream superoscillation narrative, which relies on these numbers, is therefore quantitatively mis-stated; either the values or the interpretation must be corrected. The core geometric picture of patches, folds and cusps is not affected, but the central statistical claim is not yet reliably stated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a random transverse optical polarization field as a continuous Stokes map S: R^2 → S^2, and analyses the Jacobian determinant ρ of this map. It claims that the zero lines ρ=0 are fold preimages ('crease lines'), that rank-zero points on these lines are cusp preimages ('pleat points'), and that connected components of constant sign of ρ ('patches') map to 'facets' on the Poincaré sphere. The paper further proposes a statistical distribution for ρ, Eq. (14), reports numerically estimated densities for crease lines and pleat points, discusses percolation of the patches, and sketches an interpretation in terms of origami manifolds with Delzant polytopes. The central geometric picture is clearly presented and supported by a local model, but the quantitative statistical application contains a concrete error, and the origami-manifold conclusions rely on unverified conditions that the paper itself acknowledges.","tokens_in":16571,"tokens_out":7215,"duration_ms":69392,"significance":"If the statistical statements are corrected, the paper gives a useful and original singularity-theoretic description of random polarization textures: it connects Whitney fold/cusp theory to the Skyrme-density viewpoint and makes a falsifiable prediction for the heavy-tailed distribution of the Jacobian. The local model of Eq. (10) is a clean, checkable demonstration of the cusp normal form in a polarization context, and the numerical density estimates, although approximate, provide concrete numbers that can be compared with future experiments or simulations. The paper is honest about the exploratory status of the origami-manifold section, which is a strength, but that section's conclusions are not yet established. Overall the manuscript's core geometric message is sound and deserving of publication after the statistical and framing issues are fixed.","major_comments":[{"comment":"The statement that 'the fraction of the area where |ρ| > K2 is 75% of the plane, and where |ρ| > 4K2 is 36% of the plane' is inconsistent with Eq. (14). From Eq. (14), P(|ρ| > aK2) = 2∫_{aK2}^{∞} K2^{-1}(1+|ρ|/K2)^{-3} dρ = (1+a)^{-2}, which gives 25% for a=1 and 4% for a=4. The printed values 75% and 36% equal instead the fractional contribution of these regions to ⟨|ρ|⟩, namely 2∫_{a}^{∞} u(1+u)^{-3} du. Please correct both the numbers and the interpretation; the superoscillation narrative should be rephrased accordingly. Note that after the correction the claim is actually more striking: a small area fraction (25%) carries most (75%) of the mean topological charge density.","section":"§4, sentence after Eq. (14)"},{"comment":"Equation (14) is introduced as 'straightforward to show' by 'direct application of the methods in [8,33,34]', but no derivation is given. Since this equation is the only analytic statistical result in the paper and is the basis for the quantitative claims in this section, please provide a derivation in an appendix or a precise pointer to where the result appears. This is particularly important because the same section contains the numerical misapplication described above, and the reader needs to be able to check the normalization and the algebraic tail.","section":"§4, Eq. (14)"},{"comment":"The origami-manifold interpretation requires that the null foliation of the folded symplectic form integrate to oriented circle fibres over a compact base. For the random fields studied here, the crease lines percolate across the plane (as shown in Figure 6c and discussed in §4), so the finite-loop/compact-base condition is not satisfied in general. The manuscript partly acknowledges this ('subject to some technicalities', 'More work is needed'), but the subsequent statements about Delzant polytopes and the origami template are phrased as consequences. Please explicitly label the origami-manifold and Delzant-polytope description as a conjecture or proposal for future work, and make the conditional statement 'if the origami conditions are verified' clear in the text.","section":"§5, definition of origami manifold and template"}],"minor_comments":[{"comment":"Typo: '36% if the plane' should read '36% of the plane'.","section":"§4, sentence after Eq. (14)"},{"comment":"Typo: 'correspnds' should be 'corresponds' in the sentence about large values of ρ.","section":"§3.1"},{"comment":"The caption states that the histograms are 'fitted to a straight line of −187/91 given by percolation theory', but no fitting procedure, fit range, or uncertainty is reported; please add these details so the 'plausible agreement' can be assessed.","section":"Figure 6 caption"},{"comment":"Reference [45] lists 'Macdeff D and Salamon D'; the correct spelling is 'McDuff'.","section":"References"},{"comment":"The term 'polarisation skyrmionic anyons' is introduced without a definition or a concrete criterion distinguishing them from ordinary patches. Please define the term or soften the claim, since as written it is suggestive but not well specified.","section":"Abstract and Discussion"}],"recommendation":"major_revision","confidential_remarks":"The paper is an exploratory but potentially valuable contribution. The main reason for major revision is the concrete numerical error in §4, which affects the quantitative superoscillation claims, and the unverified conditions in the origami-manifold section. The core geometric picture of folds, cusps, and patches is sound and should survive revision. The heavy reliance on the authors' own earlier works for the statistical methods is understandable, but an editor may wish to ask the authors to clarify the novelty of Eq. (14) relative to those references."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe geometric scaffold here is worth your time: viewing a random Stokes map as a folded sheet with patches, creases, folds and cusps is a genuine new way to organize random polarization fields. The local cusp model in Section 3.4 is clear and the percolation analogy to nodal domains is suggestive. That part of the paper is solid and I'd happily discuss it.\n\nBut the statistical centerpiece has a concrete error that the authors need to face. Eq. (14) gives P(|ρ|>aK2)=(1+a)^{-2}. That means the fraction of area with |ρ|>K2 is 25%, and with |ρ|>4K2 is 4%. The paper prints 75% and 36%. Those printed numbers are actually the share of the mean |ρ| contributed by those regions, not area fractions. So the \"superoscillation\" narrative built on those numbers is mis-stated. This isn't a subtle interpretation difference; it's an internal inconsistency between Eq. (14) and the sentence that follows it.\n\nThat alone would be fixable, but the problems don't stop there. Eq. (14) itself is asserted as 'straightforward' with no derivation. For a paper whose main quantitative claim lives in that equation, that's a gap. The numerical estimates of crease and pleat-point densities carry large admitted error bars and rely on indirect counting that the authors themselves say may undercount. And Section 5's origami-manifold interpretation is explicitly an outline; the technical conditions from Cannas da Silva et al. are not verified. 'More work is needed' is not a theorem.\n\nNone of this destroys the conceptual framework. The patch/facet decomposition and the singularity-theory reading seem sound, and the paper is honest about its exploratory status. But the abstract's central statistical claim is not reliably stated, and the derived numbers are wrong. As it stands, I wouldn't cite the quantitative results; I would cite the geometric picture if it survives revision.\n\nFor peer review: yes, send it. A serious referee can separate the salvageable core from the statistical errors. The authors should be asked to derive or fix Eq. (14), correct the 75%/36% mistake, and either verify or explicitly drop the origami-manifold conditions.","headline":"The patch/facet picture of random Stokes maps is a fresh idea, but the paper's headline statistics are self-contradictory and need a major correction before the results can be trusted.","tokens_in":17085,"tokens_out":2807,"would_cite":false,"duration_ms":23495,"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 random polarisation field's Stokes map to the Poincaré sphere is generically folded into facets whose boundaries are folds and cusps, the paper argues.","keywords":["Poincaré sphere","random polarisation field","Stokes map","jacobian","fold and cusp singularities","skyrmionic textures","origami manifold","percolation"],"falsifier":"One decisive check is to simulate many isotropic Gaussian random transverse fields, compute $\\rho = S\\cdot\\partial_x S\\times\\partial_y S$ on a fine grid, and histogram it: the paper predicts $P(\\rho)=K_2^{-1}(1+|\\rho|/K_2)^{-3}$, so a tail that decays faster than $|\\rho|^{-3}$, or a convergent second moment of $\\rho$, would falsify the statistical claim. A geometric falsifier is a generic realisation whose crease set contains a transverse self-intersection or an endpoint that is neither a fold nor a cusp, which would contradict the claimed Whitney structure.","tokens_in":16060,"feed_emoji":"🌀","tokens_out":9403,"duration_ms":80643,"temperature":0.7,"pith_summary":"The paper argues that a generic random transverse polarisation field, such as polarisation speckle, is secretly a folded map rather than an unstructured one. The normalised Stokes map S from the plane to the Poincaré sphere is a generic smooth map between surfaces, so its only singularities are fold lines—the crease lines where the jacobian $\\rho$ vanishes—and cusp points where fold lines meet; the connected regions where $\\rho$ keeps its sign, called patches, map to oriented facets on the sphere with positive or negative solid angle. The paper derives the exact probability density $P(\\rho)=K_2^{-1}(1+|\\rho|/K_2)^{-3}$ for isotropic Gaussian random fields, whose heavy tails come from low-intensity regions where polarisation changes extremely fast, a polarisation analogue of superoscillations. It then proposes that the plane with the folded area form $\\rho(x,y)\\,dx\\wedge dy$ is an origami manifold whose template encodes how patches nest and glue, and presents numerical evidence that patch networks follow two-dimensional percolation scaling. A sympathetic reader would care because this gives random light a skyrmion-like, quasiparticle structure without any designed beam, and connects polarisation singularities to singularity theory, percolation, and symplectic topology.","feed_headline":"Random polarisation fields fold like paper onto the Poincaré sphere","feed_subtitle":"The Stokes map of speckle light splits the plane into sign patches whose spherical images are facets bounded by folds and cusps.","key_machinery":"The load-bearing object is the jacobian $\\rho$ of the Stokes map, the Skyrme density of the polarisation texture. It is the coefficient of the folded area form $\\rho(x,y)\\,dx\\wedge dy$; its zero set is the crease lines, its sign labels the patches, its integral over a patch is the facet's solid angle, and its stationary structure locates the pleat points that map to cusps. Whitney's fold-and-cusp theorem supplies the generic singularity classification that organises this geometry, and the Gaussian random-wave ensemble supplies the statistical distribution of $\\rho$.","core_discovery":"The central claim is that every generic random transverse polarisation field carries a well-defined patch/facet decomposition of the Skyrme type. In an isotropic Gaussian random superposition of transverse plane waves, the jacobian of the Stokes map, $\\rho = S\\cdot\\partial_x S\\times\\partial_y S = \\nabla\\phi\\times\\nabla Z\\cdot e_z$, is a signed scalar whose zero set forms crease lines; $\\operatorname{sign}(\\rho)$ labels the patches of the plane, and each patch maps to a facet on the Poincaré sphere subtending a solid angle given by $\\int \\rho$ over the patch. The statistical distribution of $\\rho$ is sign-symmetric with $P(\\rho)=K_2^{-1}(1+|\\rho|/K_2)^{-3}$, so the mean absolute density is the second spectral moment $K_2$ while all higher moments diverge, reflecting superoscillatory low-intensity regions. The paper identifies the preimages of folds as crease lines and the preimages of cusps as pleat points, gives a local normal form for a crease with a pleat point, and estimates the densities of crease lines and pleat points numerically. Finally, it proposes that $(\\mathbb{R}^2, \\rho(x,y)\\,dx\\wedge dy)$ can be read as an origami manifold, with facets corresponding to Delzant polytopes in an origami template, while explicitly noting that more work is needed to make that identification fully rigorous.","pith_inferences":["I would expect the same fold/cusp analysis to organise momentum-space Stokes maps in topological photonics, where the Chern number is the degree; local folds and facet rearrangements cannot change that integer, so the patch picture gives a geometric way to see why the invariant is stable.","The divergent higher moments of $\\rho$ imply that the solid angle of a randomly selected finite facet has no well-defined variance; experiments comparing facet areas will therefore show slow, sample-dependent convergence unless conditioned on intensity.","If the origami-template conditions are eventually verified, the Delzant-polytope data would provide a combinatorial label for random polarisation textures; a testable consequence is that facet images should obey the template's vertex-gluing rules, which could be checked from measured Stokes fields.","The crease lines, being higher-order singularities independent of special polarisation states, are plausible organisers for C-point creation and annihilation events, playing the role that $\\Omega=0$ lines play for vortex-loop topology in scalar waves."],"forward_implications":["Every random speckle pattern has a patch–facet decomposition, so partial skyrmionic coverings of the Poincaré sphere are generic, not design-dependent.","The exact law $P(\\rho)=K_2^{-1}(1+|\\rho|/K_2)^{-3}$ predicts that the average absolute skyrmionic charge density in a random field is $K_2$ and that rare low-intensity regions dominate the covered solid angle.","Measurable densities follow from the analysis: crease-line density is close to the L-line density, and pleat-point density is roughly 1.79 times the density of fixed-polarisation points.","Patch areas show scaling consistent with two-dimensional percolation, making polarisation speckle a candidate laboratory system for percolation critical phenomena.","Because the description depends only on the jacobian, it is invariant under rotations of the Poincaré sphere and applies equally to maps from tori and to higher-dimensional polarisation fields."],"supporting_citations":[{"why":"Supplies Whitney's theorem that generic maps between surfaces have only fold and cusp singularities, the backbone of the patch/facet picture.","marker":"[21]"},{"why":"Establishes the isotropic random vector-wave ensemble used to model random polarisation fields.","marker":"[5]"},{"why":"Provides the polarisation-singularity statistics and C-point/L-line densities that set up the jacobian analysis.","marker":"[8]"},{"why":"Gives the random-wave second-moment and jacobian results needed to identify $K_2$ and derive $P(\\rho)$.","marker":"[33]"},{"why":"Supplies nodal-density and domain statistics methods used for crease-line density and patch percolation.","marker":"[34]"},{"why":"Proposes two-dimensional percolation universality for sign-symmetric random domains, which the patch-area scaling is compared against.","marker":"[37]"},{"why":"Defines symplectic origami manifolds and templates, the abstract structure proposed for the Stokes map.","marker":"[39]"},{"why":"Provides the Delzant-polytope description of moment images invoked for the facet/template correspondence.","marker":"[44]"}],"fun_headline_variants":["Random polarisation fields fold into skyrmionic facets on the Poincaré sphere","Speckle polarisation: patches become facets with folds and cusps","Random polarisation fields produce origami-like facets on the Poincaré sphere","Stochastic polarisation fields: skyrmion patches and origami folds","Random polarisation fields: skyrmion anyons in origami facets"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a random polarisation speckle field is a generic smooth map between two surfaces, so its only singularities are folds and cusps; the origami-template extension additionally assumes the crease lines foliate into oriented circles over a compact base, a condition the paper leaves unverified.","fun_headline_variants_meta":{"raw":{"variants":["Random polarisation fields fold into skyrmionic facets on the Poincaré sphere","Speckle polarisation: patches become facets with folds and cusps","Random polarisation fields produce origami-like facets on the Poincaré sphere","Stochastic polarisation fields: skyrmion patches and origami folds","Random polarisation fields: skyrmion anyons in origami facets"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000978,"raw_usage":{"total_tokens":4198,"prompt_tokens":1035,"completion_tokens":3163,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":651,"completion_tokens_details":{"reasoning_tokens":3063}},"tokens_in":651,"tokens_out":3163,"duration_ms":21076,"temperature":1.0,"reasoning_tokens":3063,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:23:04.516228+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One decisive check is to simulate many isotropic Gaussian random transverse fields, compute $\\rho = S\\cdot\\partial_x S\\times\\partial_y S$ on a fine grid, and histogram it: the paper predicts $P(\\rho)=K_2^{-1}(1+|\\rho|/K_2)^{-3}$, so a tail that decays faster than $|\\rho|^{-3}$, or a convergent second moment of $\\rho$, would falsify the statistical claim. A geometric falsifier is a generic realisation whose crease set contains a transverse self-intersection or an endpoint that is neither a fold nor a cusp, which would contradict the claimed Whitney structure.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies Whitney's theorem that generic maps between surfaces have only fold and cusp singularities, the backbone of the patch/facet picture."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the isotropic random vector-wave ensemble used to model random polarisation fields."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the polarisation-singularity statistics and C-point/L-line densities that set up the jacobian analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the random-wave second-moment and jacobian results needed to identify $K_2$ and derive $P(\\rho)$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies nodal-density and domain statistics methods used for crease-line density and patch percolation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Proposes two-dimensional percolation universality for sign-symmetric random domains, which the patch-area scaling is compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines symplectic origami manifolds and templates, the abstract structure proposed for the Stokes map."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Delzant-polytope description of moment images invoked for the facet/template correspondence."}],"review_version":1}