{"id":"f59c620d-8f73-4a5a-8628-174522d1ebb5","arxiv_id":"1908.02142","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The optimal birefringence distribution for star test polarimetry, minimizing PSF size while encoding polarization, is nearly the same as that of a three-point stress-engineered optic.","lead":"This paper uses a variational calculation to find the birefringence pattern in a camera pupil that encodes polarization while keeping the point-spread function as small as possible. It finds that the optimal pattern is nearly identical to a window squeezed at three points, the stress-engineered optic already used in star test polarimetry.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed optimum is only proven within the separable q3=0 ansatz of Section 7; no argument excludes a better q3≠0 or non-separable distribution, so 'optimal' in the title is stronger than demonstrated.","rationale":"The reader's weakest_assumption identifies exactly the load-bearing concern: the variational optimum is sought inside a separable q3 = 0 family, and no proof is given that the global optimum lies in that family. The paper's derivation of the Euler-Lagrange equations (16) is general, and the ansatz solution is a legitimate stationary point, but the title and abstract claim global optimality. Because the problem is nonconvex, a stationary point of the restricted family need not be the global minimizer; the paper does not supply a second-variation analysis, a uniqueness theorem, or a numerical search over the full field of birefringence distributions. The natural q3 ≠ 0 extension mentioned above satisfies two of the three constraints automatically and changes the energy functional only through positive Θ-gradient terms and the cos²Θ factor on the azimuthal term; this possibility should be checked rather than assumed away. The SEO near-equivalence conclusion may still be correct, but it would be on firmer ground if the optimality claim were either proven within the full space or explicitly qualified as optimal within the considered class. Since the reader already recommended a conditional acceptance on these grounds, the verdict should remain CONDITIONAL.","tokens_in":14115,"tokens_out":45746,"duration_ms":466019,"concrete_test":"Perform a full 2D numerical optimization of the variational problem: discretize q = (q1, q2, q3) on a polar grid with q0 = (1 - |q|²)^(1/2), minimize κ∆r² subject to ⟨β⟩ = 0 and the boundary condition Eq. (17), without imposing separability or q3 = 0. Compare the resulting minimum ∆r to the paper's 0.220λ/NA. If the full optimum is lower, the ansatz is not global; if it matches, the claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The decisive gap is in Section 7: the 'optimal' solution is obtained by postulating q(u,φ) = sinδ(u)[cos(mφ), sin(mφ), 0], i.e., q3 = 0 and a perfectly m-fold-symmetric, separable form. Substituting this ansatz into the general Euler-Lagrange equations (16) produces a valid stationary point of the full problem, but the paper gives no argument that the global minimizer of κ∆r² under the constraints must have this form. The objective is not convex (unit-sphere constraint q0² + |q|² = 1 and nonlinear constraints ⟨β⟩ = 0), so stationarity plus a scan over integer m does not imply global optimality. A concrete unexplored family is q = sinδ(u)(cosΘ(u) cosφ, cosΘ(u) sinφ, sinΘ(u)), which automatically satisfies ⟨β1⟩ = ⟨β2⟩ = 0 and reduces the azimuthal gradient energy by a factor cos²Θ, at the cost of extra Θ gradient terms; even if the optimal Θ turns out to be zero, this must be demonstrated. The paper's abstract and title claim 'the optimal birefringence distribution,' but the proof establishes at most optimality within the proposed symmetric q3 = 0 family. If a q3 ≠ 0 or non-separable distribution achieves a smaller ∆r, the near-equivalence conclusion for the SEO (5–10% in Table 1) could change.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a variational calculation to find the birefringence distribution of a pupil-plane mask that minimizes the RMS width increase of the point-spread function (PSF) in star test polarimetry, subject to constraints that ensure equal power in the two circular polarization components and hence a polarization-independent total power. The authors derive the Euler-Lagrange equations (Eqs. 16) and boundary conditions (Eq. 17), first analyze a simplified problem ignoring boundary conditions (Section 6), then propose a separable ansatz q(u,φ)=sinδ(u)[cos(mφ),sin(mφ),0] in Section 7. They solve the resulting nonlinear ODE numerically, find the optimal |m|=1 solution, and give a quadratic approximation δ(v)≈1.856v−0.922v². This solution is compared with the stress-engineered optic (SEO) with δ(v)=1.166v, showing 5–10% differences in FWHM and Strehl ratio. The paper concludes that the SEO is nearly optimal for star test polarimetry and supports this with a Fisher information analysis of measurement accuracy.","tokens_in":14396,"tokens_out":3760,"duration_ms":38496,"significance":"If the optimality claim were fully established, this paper would provide a fundamental benchmark for star test polarimetry and a rigorous explanation of why the simple stress-engineered optic is a near-ideal mask. The work has clear strengths: the variational derivation is careful and internally consistent, the Fisher information analysis is well executed, and the final comparison is quantitative and reproducible from the given equations. The numerical solution is transparent and the approximate fit is useful for practitioners. However, the central claim of global optimality is stronger than what is demonstrated, since the optimization is restricted to a specific symmetric, q3=0 family without proof that the global minimizer lies in this family. This gap affects the title and abstract claims, and also the near-equivalence conclusion for the SEO, which depends on the comparison being made against the true optimum.","major_comments":[{"comment":"The global optimality claim is not proven. The paper proposes a separable solution constrained to the q1q2 plane, q(u,φ)=sinδ(u)[cos(mφ),sin(mφ),0], and solves the resulting ODE (Eq. 25) for the half-retardance. This yields a stationary point of the full variational problem only within this restricted family. The objective κ∆r² and the constraints ⟨β⟩=0 are nonconvex (the unit-sphere condition q0²+|q|²=1 is nonlinear), so stationarity within an ansatz does not imply global optimality. A concrete unexplored family is q=sinδ(u)(cosΘ(u)cosφ, cosΘ(u)sinφ, sinΘ(u)), which automatically satisfies ⟨β1⟩=⟨β2⟩=0 and reduces the azimuthal gradient energy by a factor cos²Θ at the cost of additional Θ-gradient terms. Even if the optimal Θ turns out to be zero, that must be demonstrated. As written, the paper establishes at most optimality within the proposed symmetric q3=0 family, and the abstract and title claim 'the optimal birefringence distribution' is therefore too strong. I recommend either providing a symmetry or convexity argument that the global minimizer has this form, or substantially tempering the claims to 'optimal within the proposed family.' Since the SEO near-equivalence conclusion (Table 1) rests on this comparison, the gap is load-bearing.","section":"Section 7, Eq. (21) and the surrounding text"},{"comment":"The statement 'there is no loss of generality in this choice, because other solutions can be found through cascading with uniform birefringent plates' is unsupported. Cascading the BM with uniform wave plates can indeed generate other solutions, but it is not shown that every candidate distribution with q3≠0 (or non-separable azimuthal dependence) can be mapped to a q3=0 solution with equal or smaller ∆r while still satisfying the constraints ⟨β⟩=0. The constraints and the RMS width increment are not obviously invariant under such cascading. This matters because the ansatz q3=0 is a key restriction in the derivation of the claimed optimum. The paper should either prove the reduction or explicitly state that q3≠0 distributions are outside the scope of the optimization.","section":"Section 6, text following Eq. (21)"}],"minor_comments":[{"comment":"The phrase 'of a a glass window' contains a duplicated article 'a'.","section":"Section 10, concluding remarks"},{"comment":"The word 'birefingence' is a typo; it should be 'birefringence'.","section":"Section 10, concluding remarks"},{"comment":"When stating that 'the first one, c = 1.166/NA' is the root that minimizes the width increase, it would be clearer to specify that this is the first nonzero root of Eq. (26), and to note that the root search is over cNA>0.","section":"Section 8, Eq. (26)"},{"comment":"The quadratic approximation 1.856v−0.922v² is presented with R²=0.9996, but the paper does not state how the fitting was performed or over what range. Since the numerical solution is already available, the fit is only a presentation aid and could be described as such.","section":"Section 7, Eq. (25) and the fit"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope and the core derivations are sound, but the overclaiming of global optimality is a substantive issue that must be addressed before publication. If the authors can justify the restriction to the separable q3=0 family (for example, via a symmetry argument or by extending the numerical search to a broader class of distributions), the paper would be a solid contribution. If not, the title and abstract should be revised to reflect optimality within the proposed family. I would not recommend rejection, because the engineering conclusion about the SEO's near-optimality is plausible and the variational machinery is valuable; however, the current claim is too strong for the evidence presented."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper does a real variational calculation, gets a defensible practical conclusion, and the flaw is an over-strong title/claim rather than a wrong result.\n\nWhat's new and good: previous star test polarimetry work took the SEO as given; this paper asks what the best birefringence distribution is if you minimize the RMS width increase of the PSF under equal-power constraints. The Euler-Lagrange setup is careful, the Fisher-information follow-up is a nice check, and the final comparison table is honest: the SEO is 5-10% off the derived optimum in FWHM and Strehl. The near-equivalence is a first-principles explanation of something that was empirical. That's a genuine contribution, and the paper is well written.\n\nThe soft spot: global optimality isn't established. Section 7 postulates q(u,φ)=sinδ(u)[cos(mφ), sin(mφ),0], i.e., q3=0 and a separable azimuthal form, then solves the ODE for δ within that family. That's a restricted ansatz. The objective is not convex, and the constraint set is nonlinear, so stationarity within a family doesn't give the global minimum. The stress-test note's alternative q = sinδ(u)(cosΘ(u)cosφ, cosΘ(u)sinφ, sinΘ(u)) is a concrete unexplored direction; even if the optimum turns out to be Θ=0, it needs a sentence or an appendix. So 'optimal' in the title and abstract is stronger than what's proven. 'Optimal within the separable m=1 class' would be accurate; or add a proof/justification.\n\nAlso minor: no code or numerical data, but equations are complete enough that replication is straightforward. Self-citation to [12] for the width formula is justified since it's a published derivation.\n\nMy recommendation: this deserves peer review. The core calculation is sound, the practical conclusion (SEO near-optimal) is well supported for the considered family, and the over-claim is fixable with wording or a modest extra argument. The paper will be useful to anyone working in single-shot imaging polarimetry or fluorescence microscopy with birefringent pupils. I'd encourage review with a request to address the ansatz gap.","headline":"A clean variational derivation showing the stress-engineered glass optic is near-optimal for star test polarimetry; just don't read 'optimal' as a proven global result.","tokens_in":14947,"tokens_out":1752,"would_cite":true,"duration_ms":17329,"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 derives the optimal birefringence distribution for star test polarimetry and shows a stressed-glass optic nearly achieves it.","keywords":["star test polarimetry","birefringence mask","stress-engineered optic","point-spread function","variational optimization","Stokes parameters","Fisher information","single-shot polarimetry"],"falsifier":"Solve the full variational equations, Eqs. (16), without the separable ansatz, allowing arbitrary azimuthal dependence and $q_3\\neq0$, and look for a stationary solution with $\\langle\\boldsymbol{\\beta}\\rangle_A=0$ and $\\Delta r < 0.220\\lambda/\\mathrm{NA}$; a numerical search that finds any admissible mask with a smaller RMS-width increase would disprove the claimed optimum.","tokens_in":13851,"feed_emoji":"🔭","tokens_out":7993,"duration_ms":76210,"temperature":0.7,"pith_summary":"Star test polarimetry encodes the polarization of an incoming field into the shape of the point-spread function by placing a spatially varying birefringent mask in the pupil. This paper treats the choice of that mask as a variational problem: minimize the increase in RMS PSF width caused by the birefringence, subject to the constraint that the two circular polarization components carry equal power. The optimal solution has half-retardance closely fitted by $1.856v - 0.922v^2$ across the normalized pupil radius, reaching a width increase of $0.220\\lambda/\\mathrm{NA}$. That optimum is within 5 to 10 percent in FWHM and Strehl ratio of the linear profile $\\delta = 1.166v$ produced by a stress-engineered optic, so the widely used SEO is nearly optimal. The result matters because it shows that a simple stressed-glass window, not a nanofabricated device, is close to the best possible pupil mask for single-shot imaging polarimetry.","feed_headline":"A stressed-glass pupil mask is nearly the optimal polarimetry optic","feed_subtitle":"A variational calculation places the best birefringence profile within 5–10% of the stressed glass already in use.","key_machinery":"The central object is the Poincaré-sphere/quaternion representation of spatially varying birefringence, $\\boldsymbol{\\hat{q}}(u)=(q_0,q_1,q_2,q_3)$ with $|\\boldsymbol{\\hat{q}}|=1$, in which the half-retardance and axis orientation are encoded in the three-vector $\\mathbf{q}$. The variational functional $\\Delta r^2=\\kappa^{-1}\\int A^2\\|\\nabla\\boldsymbol{\\hat{q}}\\|^2d^2u$ ties the PSF broadening directly to the gradient of this unit vector, and the Lagrange-multiplier constraint enforces equal power in the two polarization components. Substitution of the azimuthal ansatz turns the problem into the ordinary differential equation $\\bar\\delta''(v)+\\bar\\delta'(v)/v+(\\mathrm{NA}^2\\Lambda_3 - 1/(2v^2))\\sin[2\\bar\\delta(v)]=0$, whose solution is the fitted quadratic profile. A secondary mechanism is the Fisher information matrix, which converts the PSF shapes into predicted uncertainties in the retrieved Stokes parameters.","core_discovery":"Working in the circular polarization basis and representing the Jones matrix by the unit quaternion $\\boldsymbol{\\hat{q}}=(q_0,q_1,q_2,q_3)$, the paper derives Euler-Lagrange equations for the birefringence distribution that minimizes $\\Delta r^2 = \\kappa^{-1}\\int A^2 \\|\\nabla \\boldsymbol{\\hat{q}}\\|^2\\,d^2u$ under the equal-power constraint $\\langle\\boldsymbol{\\beta}\\rangle_A=0$. For a circular hard aperture, the separable ansatz $\\mathbf{q}(u,\\phi)=\\sin\\delta(u)[\\cos\\phi,\\sin\\phi,0]$ reduces these equations to a single nonlinear boundary-value problem for the half-retardance $\\delta(v)$, Eq. (25), with boundary conditions $\\delta(0)=0$ and $\\delta'(1)=0$. The numerical solution is essentially quadratic, $\\bar\\delta(v)\\approx 1.856v-0.922v^2$, with an RMS-width increase of $0.220\\lambda/\\mathrm{NA}$. Compared with the stress-engineered optic's linear profile $\\bar\\delta(v)=1.166v$, the optimal mask improves FWHM and Strehl by only 5 to 10 percent, and the Fisher-information analysis gives essentially identical expected Stokes-parameter errors, so the conclusion is that the SEO is a near-optimal implementation.","pith_inferences":["If the separable $q_1q_2$ ansatz is relaxed to allow $q_3\\ne0$ or non-separable azimuthal dependence, the true optimum could beat $0.220\\lambda/\\mathrm{NA}$, but the 5 to 10 percent gap to the SEO suggests any further gain would be small.","The fitted quadratic profile $1.856v-0.922v^2$ could be tested directly by fabricating a mask with that retardance profile and comparing its PSF FWHM to an SEO's; the paper's numbers predict almost no visible difference.","The same variational machinery could be reapplied to apodized pupils, high-NA vectorial focusing, or other PSF-size metrics; each would shift the optimal $\\delta(v)$ slightly, likely preserving the SEO's near-optimality."],"forward_implications":["A stress-engineered optic with threefold edge stress is, to within 5 to 10 percent in FWHM and Strehl, an optimal birefringence mask for star test polarimetry, so no nanofabricated or programmable mask is needed for near-optimal performance.","The optimal solution and the SEO both satisfy the equal-power constraint, so the Stokes parameters can be read from a single polarization component's PSF shape without calibrating total power.","Encoding polarization information costs roughly a 60 percent increase in FWHM and a 50 percent drop in Strehl relative to the diffraction-limited PSF; that is the price of single-shot polarimetry.","When both circular components are imaged, the four PSF contributions are nearly orthonormal under the weight $1/I_0^{(1)}(x)$, which simplifies maximum-likelihood retrieval of the Stokes parameters.","Azimuthal orders $|m|>1$ produce larger PSFs in this framework, so the unit-vortex geometric phase written by the mask is the compact choice for polarization encoding."],"supporting_citations":[{"why":"supplies the Poincaré-sphere/quaternion representation of the birefringent Jones matrix and the RMS-width formula on which the variational calculation is built.","marker":"[12]"},{"why":"introduces star test polarimetry with stress-engineered optics and supplies the experimental setup the optimal solution is compared with.","marker":"[6]"},{"why":"gives the linear half-retardance form $\\delta(u)=cu$ for the SEO and the imaging demonstration that defines the performance baseline.","marker":"[7]"},{"why":"provides the stereographic-projection derivation of the separable solution family $\\mathbf{q}(u,\\phi)=2bu[\\cos\\phi,\\sin\\phi,0]/(1+b^2u^2)$ used as the starting point for the constrained optimum.","marker":"[13]"},{"why":"supplies the maximum-likelihood/Fisher-information formalism used to estimate Stokes-parameter uncertainties for the optimal and SEO masks.","marker":"[15]"}],"fun_headline_variants":["Stressed glass gets within 10% of optimal polarimetry","Variational calc finds stressed glass near-optimal for star test","Glass stress pattern nearly optimal for imaging polarimetry","Optimal birefringence profile is nearly that of stressed glass"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation only searches among masks whose birefringence axis stays in the $q_1q_2$ plane and winds once around the pupil; if the true optimum needs a different axis orientation or a higher-order azimuthal variation, the claimed optimality could fail.","fun_headline_variants_meta":{"raw":{"variants":["Stressed glass gets within 10% of optimal polarimetry","Variational calc finds stressed glass near-optimal for star test","Glass stress pattern nearly optimal for imaging polarimetry","Optimal birefringence profile is nearly that of stressed glass"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00052,"raw_usage":{"total_tokens":2513,"prompt_tokens":937,"completion_tokens":1576,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":553,"completion_tokens_details":{"reasoning_tokens":1506}},"tokens_in":553,"tokens_out":1576,"duration_ms":12388,"temperature":1.0,"reasoning_tokens":1506,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:23:23.878024+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the full variational equations, Eqs. (16), without the separable ansatz, allowing arbitrary azimuthal dependence and $q_3\\neq0$, and look for a stationary solution with $\\langle\\boldsymbol{\\beta}\\rangle_A=0$ and $\\Delta r < 0.220\\lambda/\\mathrm{NA}$; a numerical search that finds any admissible mask with a smaller RMS-width increase would disprove the claimed optimum.","supporting_citations":[{"cited_title":"Poincaré sphere representation for spatially varying birefringence,","cited_arxiv_id":null,"evidence_quote":"supplies the Poincaré-sphere/quaternion representation of the birefringent Jones matrix and the RMS-width formula on which the variational calculation is built."},{"cited_title":"Star test polarimetry using stress-engineered optical elements,","cited_arxiv_id":null,"evidence_quote":"introduces star test polarimetry with stress-engineered optics and supplies the experimental setup the optimal solution is compared with."},{"cited_title":"Imaging the polarization of a light field,","cited_arxiv_id":null,"evidence_quote":"gives the linear half-retardance form $\\delta(u)=cu$ for the SEO and the imaging demonstration that defines the performance baseline."},{"cited_title":"Description and applications of space-variant polarization states and elements,","cited_arxiv_id":null,"evidence_quote":"provides the stereographic-projection derivation of the separable solution family $\\mathbf{q}(u,\\phi)=2bu[\\cos\\phi,\\sin\\phi,0]/(1+b^2u^2)$ used as the starting point for the constrained optimum."},{"cited_title":"Tutorial: Maximum likelihood estimation in the context of an optical measurement","cited_arxiv_id":"1806.04503","evidence_quote":"supplies the maximum-likelihood/Fisher-information formalism used to estimate Stokes-parameter uncertainties for the optimal and SEO masks."}],"review_version":1}