{"id":"a36cbef0-15ba-4b92-b6e2-70c90c508389","arxiv_id":"2506.11352","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Phase retrieval in paraxial optics is reduced, in the ray optics limit, to an optimal transport problem whose solution equals the leading-order Wigner distribution, yielding vortex-free beam shaping.","lead":"This thesis solves a laser beam shaping problem by computing the phase of a light field from intensity measurements using optimal transport, the mathematics of moving one distribution into another at minimum cost. It also shows that the optimal transport plan is the ray optics limit of the Wigner distribution, giving a new way to think about phase retrieval.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"OT-plan/Wigner identity and the vortex-free guarantee are established only for convex, ray-optics-valid phases; the abstract states them unconditionally, and the cubic experiment shows the OT plan is a convex surrogate, not the retrieved zero-order Wigner distribution.","rationale":"The reader's weakest assumption already identifies convexity and beta >> 1 as the key validity condition. I agree that this is the central load-bearing assumption. I add a sharper formulation: the equality W_f^(0) = Gamma is not merely an approximation that degrades for non-convex phases; it is a selection rule. Optimal transport with quadratic cost always returns the unique convex transport plan, while the zero-order Wigner distribution of a non-convex phase is a different, possibly multi-branched measure with the same marginals. Thus the deeper claim is a convexity prior on the unknown field, not a general property of all fields. The paper does show awareness of this in Section 3.2.3 and in the cubic experiment, so the internal argument is not fraudulent, but the abstract and the vortex-free marketing overstate the scope. The numerical experiments included in the thesis support the conditional conclusion for smooth convex ray-optics targets; the memory-reduction result is a genuine technical contribution; and the one-shot/two-shot beam estimation results are plausible. The correct response is to keep the conditional verdict and require that the claims be explicitly restricted to the convex ray-optics regime, with code and benchmark data provided for verification.","tokens_in":62486,"tokens_out":8974,"duration_ms":117285,"concrete_test":"Reproduce the Section 3.4.3 cubic-phase experiment at alpha=2, beta=-2 on the natural lattice with N=128, using the paper's OT+GS pipeline of Section 4.6. Count the number of phase singularities with nonzero winding in the final OT+GS phase, both before and after Gerchberg-Saxton refinement, and record the intensity loss (B.25). If any singularity appears inside the support of g^2, or if the OT+GS loss is not better than the random-GS baseline of Figure 4.11, then the unconditional vortex-free and state-of-the-art claims fail for non-convex phases. Additionally, compute the L1 distance between the OT plan Gamma and the true zero-order Wigner distribution g^2 delta(phi'(u)-mu); a nonzero distance confirms that Eq. 3.58 is an assumption, not a retrieval, outside the convex regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central theoretical identity, W_f^(0) = Gamma (Eq. 3.58), is derived by a stationary-phase expansion that assumes the unknown phase is strictly convex or concave and that beta = 2*pi*R_in*R_out/(f*lambda) is large (Sections 3.2.2 and 3.2.3). Under these assumptions, the OPTIMAL TRANSPORT plan between g^2 and G^2 coincides with the zero-order Wigner distribution of the unknown field. But the claim is presented in the abstract as a property of 'the unknown field' without the convexity caveat. For any non-convex phase, e.g. the cubic phase in Eq. 3.72, the zero-order Wigner measure g^2 delta(phi'(u)-mu) is not the Brenier OT plan: the OT plan is the unique monotone/convex rearrangement of the marginals, while the true ray-optics phase-space measure may be multiply-branched. The paper's own Figures 3.8-3.10 demonstrate that for beta != 0 the OT-retrieved Wigner distribution differs from the ground truth and the intensity loss increases. Therefore the strong claims 'completely bypasses the creation of phase vortices' and 'state-of-the-art ... in terms of both accuracy and efficiency' are only supported for smooth, convex-phase, ray-optics-valid targets. Section 3.2.3 candidly notes that the method is best for smooth light fields and not for diffraction-limited spot arrays, but the abstract and the vortex-free phrasing do not carry this scope restriction. This is a scoping gap rather than an internal inconsistency, but it is load-bearing because the headline contribution is precisely that the OT plan recovers the Wigner distribution of an unknown field.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a ray-optics reduction of phase retrieval. It derives the Monge-Ampere equation from stationary phase, maps that equation to optimal transport, identifies the zero-order term of the Wigner distribution with the OT transport plan, generalizes the argument to fractional Fourier transforms, and uses entropy-regularized OT with a Gaussian-kernel convolution to reduce the memory cost of hologram computation from O(n^4) to O(n^2). The same machinery is applied to beam estimation, and a final chapter sketches a quantum-learning reinterpretation of the phase-generation problem. The text is written as a Stanford BS thesis and includes substantial review material on wave optics, Wigner distributions, and iterative phase-retrieval algorithms.","tokens_in":62825,"tokens_out":6545,"duration_ms":79816,"significance":"If restricted to its demonstrated regime, the paper is a valuable contribution: the stationary-phase derivation of Monge-Ampere, the explicit error parameter beta = 2*pi*R_in*R_out/(f*lambda), and the identification of the OT plan with the zero-order Wigner term are original and pedagogically useful bridges between holography and optimal transport. The O(n^2)-memory Sinkhorn implementation with Gaussian convolutions is a concrete algorithmic improvement over the earlier O(n^4) formulation, and the open-source implementation makes the claims checkable. However, the significance is bounded by the conditions of the theory (strictly convex or concave phase, beta >> 1) and by the fact that the headline performance claims rest on a small number of favorable numerical examples rather than systematic benchmarks. The central theoretical chain is coherent, but the scope of several claims in the abstract and conclusions exceeds what the manuscript proves.","major_comments":[{"comment":"The identity W_f^(0) = Gamma (Eq. 3.58) and the associated vortex-free claim are established only in the ray-optics regime: the stationary-phase reduction assumes a strictly convex or concave phase and beta = 2*pi*R_in*R_out/(f*lambda) >> 1 (Secs. 3.2.2-3.2.3). The paper's own cubic-phase experiment (Eq. 3.72; Figs. 3.8-3.10) shows visibly degraded retrieval and increased intensity loss for a non-convex phase. The abstract, however, presents these results for 'the unknown complex-valued field' without these restrictions. This is a scoping gap rather than an internal inconsistency, but it is load-bearing because the abstract's unconditional phrasing is the headline of the paper; the scope conditions should be stated there and in the conclusions.","section":"Abstract, Secs. 3.2.3 and 3.4.3"},{"comment":"The claim that OT-initialized GS/MRAF yields state-of-the-art accuracy and efficiency and 'completely bypasses' vortex formation is supported by a small set of numerical examples: the Gaussian-to-ring experiment of Sec. 4.2.2 and the 128x128 comparison from [13] reproduced in Fig. 4.13. In Fig. 4.13(j) the MRAF mixing parameter m is hand-tuned to 0.48 to obtain epsilon = 5.95e-16 at 85.15% efficiency, while OT+GS gives 2.58% at 99.91%; this is a favorable demonstration, not a general benchmark. To support the stated generality, the paper should either explicitly restrict the performance claims to smooth, convex, ray-optics-valid targets or provide systematic sweeps over input-output pairs and the two free parameters (epsilon, m).","section":"Sec. 4.6 and Fig. 4.13"},{"comment":"The claim in Sec. 3.4.4 and the abstract that the method 'can be used to retrieve the phase-space transformation of any unknown quadratic phase system' is not demonstrated. The section derives only the coefficient A/(2B) for a single collimated input and then states that one or two additional diversity images should suffice, with the 'exact implementations of this idea' left to future work. No full ABCD-matrix retrieval algorithm, numerical experiment, or error analysis is given. This statement should be reframed as a conjecture or the missing construction should be supplied.","section":"Sec. 3.4.4"},{"comment":"The quantum-learning equivalence in Problem 5 is proved only in the infinite-N, zero-error, one-dimensional setting, and the text explicitly concedes that the finite-N, error-tolerant version 'we have not proven it yet.' The abstract's mention of a quantum learning framework and the chapter's framing as a 'deep theoretical connection' therefore outrun the proved content. This secondary claim should be explicitly labeled as a conjecture, or the finite-N proof and the tomographic algorithm with guarantees should be added.","section":"Sec. 5, Problem 5"}],"minor_comments":[{"comment":"The stationary-phase formula in Eq. (3.14) omits the standard nondegeneracy condition det(nabla^2 psi) != 0; strict convexity alone does not exclude a vanishing Hessian at an isolated point (e.g., x^4), so the theorem's assumptions should be stated precisely.","section":"Theorem 3.2.1"},{"comment":"There is a missing parenthetical grouping in the term 'cot alpha pi (x^2 + y^2)'; it should read pi cot(alpha) (x^2 + y^2) or an equivalent unambiguous expression.","section":"Eq. (3.41)"},{"comment":"The text refers to 'Parallel's identity' in Eq. (A.5); this should be 'Parseval's identity.'","section":"Appendix A.5"},{"comment":"The footnote stating that the convergence conditions were 'summarized from the paper [15] by GPT-o3' is not an appropriate citation style for a scientific article and should be rewritten as a normal scholarly summary with proper attribution.","section":"Sec. 4.2.1"},{"comment":"The bullet list in Sec. 4.6 says the method has a single hyperparameter epsilon, but when OT is combined with MRAF, the mixing parameter m is also user-selected, as shown in Fig. 4.13(j); the hyperparameter discussion should acknowledge this coupling.","section":"Sec. 4.6"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a BS thesis with substantial review material and with clear overlap with the author's own prior work [13]. For a journal submission, the novelty relative to [13] and to the Peyre/Cuturi convolution formulation should be stated more explicitly, and the thesis-style narrative (including acknowledgments and GPT-generated summary) should be converted to standard article form. I see no evidence of misconduct, but the fit between the current document and a regular journal article is imperfect."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this is a solid methods thesis, not a breakthrough. The core idea—use optimal transport to compute a convex, unwrapped phase that seeds Gerchberg–Saxton or MRAF—is genuinely useful and supported by the numerical experiments shown. The derivation of the Monge–Ampere equation from stationary phase is legitimate, and identifying the zero-order Wigner term with the OT plan is a nice observation, but it only holds for strictly convex phases in the ray-optics limit (beta >> 1). The paper itself admits this in Section 3.2.3 and demonstrates the breakdown on a cubic phase in Figures 3.8–3.10. The abstract, however, states the identity unconditionally and promises 'completely bypasses' vortices and 'state-of-the-art' accuracy without those caveats. That is a scoping gap, not an internal contradiction, but it is load-bearing because it affects how the headline contribution should be read.\n\nWhat is genuinely new and good: the fractional Fourier transform generalization is a real extension of prior OT phase generation, and the O(n^2) memory reduction via the convolution trick is a practical improvement that matters for high-resolution holograms. The Sinkhorn–Knopp convolution formulation is clearly explained and the complexity analysis is honest. The numerical comparison of OT, GS, and MRAF is also credible for smooth targets, and the paper demonstrates an order-of-magnitude loss improvement when OT seeds GS.\n\nSoft spots, in proportion: the quantum learning chapter is explicitly unfinished—the author says the finite-N proof is not done and calls the vortex-free claim for that approach 'hopeful.' The ABCD retrieval section is also a sketch, not a proven method, with the exact implementation left to future work. The abstract's claim about retrieving any quadratic phase system goes beyond what the body supports. There are no error bars or a linked code repository with a commit hash, which makes the 'state-of-the-art' claim hard to verify independently.\n\nWho this is for: people working on beam shaping for neutral atom arrays or ultracold atom experiments will get real value from the OT initialization and the memory-efficient implementation. The theory is also a clean exposition of ray-optics limits of Wigner distributions.\n\nMy recommendation: this deserves a serious referee. It should not be desk-rejected. But it needs major revision to scope the abstract correctly, move the unfinished quantum and ABCD claims to clearly marked speculation, and provide reproducible code and benchmarks. I would not block publication on the theory—the central derivation holds under its stated assumptions—but the claims need to match the evidence.","headline":"A credible OT-based phase initialization with real computational improvements, but the abstract overstates the scope of the Wigner–OT identity and several headline claims are explicitly left unfinished in the body.","tokens_in":63376,"tokens_out":1713,"would_cite":true,"duration_ms":27125,"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 optimal-transport plan between two measured intensities is the ray-optics limit of the light field's Wigner distribution, and seeding phase-retrieval algorithms with it eliminates vortex stagnation.","keywords":["phase retrieval","optimal transport","Wigner distribution","holography","beam shaping","Gerchberg-Saxton","Sinkhorn algorithm","fractional Fourier transform"],"falsifier":"Choose a target phase $\\phi(u,v)=\\alpha(u^2+v^2)+\\beta u^3$ with $\\beta$ comparable to $\\alpha$, propagate exactly to get $G^2$, run OT initialization followed by Gerchberg-Saxton refinement, and count vortices and measure intensity loss; the paper's Figures 3.8-3.10 already show the loss growing with $\\beta$, so a regime where the vortex-free improvement disappears would falsify the claim outside the convex/ray-optics domain.","tokens_in":62244,"feed_emoji":"🌀","tokens_out":12698,"duration_ms":131112,"temperature":0.7,"pith_summary":"The paper claims that the phase-retrieval problem behind laser beam shaping can be approximated, in the ray-optics limit, by an optimal transport (OT) problem, and that the OT plan is exactly the zeroth-order term of the Wigner distribution of the unknown field: $W_f^{(0)} = \\Gamma$. This matters because OT is convex and vortex-free, whereas standard iterative solvers such as Gerchberg-Saxton stagnate by forming phase vortices. The paper reports that using the entropically regularized OT solution as an initialization for Gerchberg-Saxton or MRAF gives an order-of-magnitude improvement in intensity loss with no vortex formation and near-total laser efficiency. The same machinery is extended to fractional Fourier transforms, giving a two-measurement beam estimation protocol and a way to retrieve the phase-space effect of an unknown quadratic-phase optical system. The solver's memory cost is reduced from $O(n^4)$ to $O(n^2)$.","feed_headline":"Optimal transport kills phase vortices in hologram solvers","feed_subtitle":"Seeding phase-retrieval with an optimal-transport plan gives vortex-free, order-of-magnitude more accurate holograms.","key_machinery":"The load-bearing object is the optimal transport plan $\\Gamma$ between the two measured intensities, obtained by entropy-regularized optimal transport. In the ray-optics limit the transport plan coincides with the zero-order Wigner term $\\Gamma(u,\\mu) = g(u)^2\\delta(\\mu-\\nabla\\phi(u))$, a bundle of rays; the gradient of the phase $\\nabla\\phi$ is the transport map, and integrating it gives the convex, unwrapped phase that seeds the iterative refinement. The practical engine is the Sinkhorn-Knopp algorithm in the form of repeated Gaussian convolutions, which enforces the marginal constraints with a single hyperparameter $\\epsilon$ and avoids storing the $n^4$ transport plan.","core_discovery":"The central claim is that the optimal transport plan $\\Gamma$ between the squared input intensity $g^2$ and the squared target intensity $G^2$, under quadratic cost, is exactly the zeroth-order stationary-phase term of the Wigner distribution of the unknown field, $W_f^{(0)} = \\Gamma$, with $\\Gamma(u,\\mu)=g(u)^2\\delta(\\mu-\\phi'(u))$ in one dimension. Solving the entropy-regularized OT problem therefore produces a convex, unwrapped phase with no $2\\pi$ branch cuts, and using that phase to initialize Gerchberg-Saxton or MRAF eliminates the vortex formation that otherwise stalls convergence. The paper demonstrates numerically that the reconstructed Wigner distribution from the OT-seeded phase matches the ground truth for flat and quadratic phases, and shows how the same reduction applies when the two planes are related by a fractional Fourier transform instead of a plain Fourier transform. On the algorithmic side, writing the Sinkhorn-Knopp updates as Gaussian convolutions reduces the memory cost from $O(n^4)$ to $O(n^2)$ and the runtime from $O(n^4/\\epsilon^2)$ to $O(n^2\\log n/\\epsilon^2)$. The same ray-optics reduction yields a two-shot beam-estimation protocol from two fractional Fourier intensity measurements, and the Hermite-Gaussian expansion of the problem is reinterpreted as a quantum state tomography task.","pith_inferences":["The identity $W_f^{(0)}=\\Gamma$ suggests that the difference $W_\\psi-\\Gamma$ between the reconstructed Wigner distribution and the transport plan directly measures diffraction and interference corrections, so the norm of that difference could serve as a validity diagnostic for when the ray-optics initialization is trustworthy.","The convexity limitation suggests a windowed or locally convex extension: decompose a non-convex target into convex patches, solve OT on each patch, and stitch the phases; the paper's cubic-phase experiments indicate exactly where such a scheme would be stress-tested.","The Hermite-Gaussian/quantum learning formulation implies a concrete algorithmic relaxation not proven in the paper: promote the pure-state vector to a rank-$N$ density matrix in the same constraints, which should make the phase-generation problem feasible even when no pure state exists, at the cost of producing a mixed beam."],"forward_implications":["OT-initialized Gerchberg-Saxton produces vortex-free phases on smooth targets and reaches an intensity loss roughly an order of magnitude below random or inverse-Fourier initialization, with no efficiency sacrifice.","Composing OT with MRAF gives a controllable accuracy-efficiency trade-off: near-zero error in a chosen signal region at the cost of a tunable fraction of laser power.","The convolution form of Sinkhorn-Knopp reduces memory from $O(n^4)$ to $O(n^2)$ and runtime from $O(n^4/\\epsilon^2)$ to $O(n^2\\log n/\\epsilon^2)$, making $1024\\times1024$ and larger hologram computations feasible.","Two diversity images at different fractional Fourier angles are sufficient for a closed OT-based estimate of an unknown beam's amplitude and phase, with error set by the angle separation.","The same OT/Wigner machinery can retrieve the phase-space (ABCD) transformation of an unknown quadratic-phase optical system from projective intensity measurements."],"supporting_citations":[{"why":"The earlier coauthored paper whose OT phase-generation and two-shot beam-estimation pipeline this work extends, and whose $O(n^4)$ memory bottleneck is reduced here.","marker":"[13]"},{"why":"Supplies the geometric/ray-optics beam-shaping precedent and the dimensionless validity parameter controlling when the stationary-phase limit is accurate.","marker":"[24]"},{"why":"Provides the fractional Fourier transform, Wigner distribution, and linear-canonical transform formalism the entire reformulation is built on.","marker":"[25]"},{"why":"Gives the theorem that an optimal transport map with quadratic cost is the gradient of a scalar satisfying the Monge-Ampere equation.","marker":"[39]"},{"why":"Introduces entropy-regularized optimal transport, whose dual formulation leads to the Sinkhorn-Knopp solver used in the algorithm.","marker":"[47]"},{"why":"Proves the matrix scaling theorem that underwrites the Sinkhorn-Knopp iterations enforcing the intensity marginals.","marker":"[48]"},{"why":"Defines the Gerchberg-Saxton baseline whose vortex stagnation the OT initialization is designed to cure.","marker":"[14]"},{"why":"Analyzes convergence and failure to converge of Gerchberg-Saxton, motivating the need for a convex initialization.","marker":"[15]"},{"why":"Defines the MRAF algorithm with signal and noise regions, the accuracy-efficiency trade-off that OT seeding improves.","marker":"[1]"}],"fun_headline_variants":["Optimal transport erases vortices in hologram phase retrieval","No more vortices: optimal transport speeds up holography","Ray-optics trick kills vortex stall in phase retrieval","Optimal transport seeds vortex-free phase retrieval","Holography without vortex traps via optimal transport"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the unknown phase is a single smooth, bowl-shaped surface (strictly convex or concave) and that the beam features are large enough relative to focal length and wavelength that diffraction can be neglected; the paper's own cubic-phase simulation shows degraded Wigner estimation and increased intensity loss when a non-convex term is present.","fun_headline_variants_meta":{"raw":{"variants":["Optimal transport erases vortices in hologram phase retrieval","No more vortices: optimal transport speeds up holography","Ray-optics trick kills vortex stall in phase retrieval","Optimal transport seeds vortex-free phase retrieval","Holography without vortex traps via optimal transport"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000291,"raw_usage":{"total_tokens":1744,"prompt_tokens":1034,"completion_tokens":710,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":650,"completion_tokens_details":{"reasoning_tokens":631}},"tokens_in":650,"tokens_out":710,"duration_ms":7182,"temperature":1.0,"reasoning_tokens":631,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:10:16.479337+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Choose a target phase $\\phi(u,v)=\\alpha(u^2+v^2)+\\beta u^3$ with $\\beta$ comparable to $\\alpha$, propagate exactly to get $G^2$, run OT initialization followed by Gerchberg-Saxton refinement, and count vortices and measure intensity loss; the paper's Figures 3.8-3.10 already show the loss growing with $\\beta$, so a regime where the vortex-free improvement disappears would falsify the claim outside the convex/ray-optics domain.","supporting_citations":[{"cited_title":"Van de Graaff, Jan Rudolph, and Jason M","cited_arxiv_id":null,"evidence_quote":"The earlier coauthored paper whose OT phase-generation and two-shot beam-estimation pipeline this work extends, and whose $O(n^4)$ memory bottleneck is reduced here."},{"cited_title":"Dickey.Laser Beam Shaping: Theory and Techniques","cited_arxiv_id":null,"evidence_quote":"Supplies the geometric/ray-optics beam-shaping precedent and the dimensionless validity parameter controlling when the stationary-phase limit is accurate."},{"cited_title":"Murat Ozaktas, M","cited_arxiv_id":null,"evidence_quote":"Provides the fractional Fourier transform, Wigner distribution, and linear-canonical transform formalism the entire reformulation is built on."},{"cited_title":"The Monge–Ampère equation and its link to optimal transportation","cited_arxiv_id":null,"evidence_quote":"Gives the theorem that an optimal transport map with quadratic cost is the gradient of a scalar satisfying the Monge-Ampere equation."},{"cited_title":"Sinkhorndistances: Lightspeedcomputationofoptimaltransportationdistances, 2013","cited_arxiv_id":null,"evidence_quote":"Introduces entropy-regularized optimal transport, whose dual formulation leads to the Sinkhorn-Knopp solver used in the algorithm."},{"cited_title":"A Relationship Between Arbitrary Positive Matrices and Doubly Stochastic Matrices","cited_arxiv_id":null,"evidence_quote":"Proves the matrix scaling theorem that underwrites the Sinkhorn-Knopp iterations enforcing the intensity marginals."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Gerchberg-Saxton baseline whose vortex stagnation the OT initialization is designed to cure."},{"cited_title":"Alternating projections with applications to gerchberg-saxton error reduction, 2021","cited_arxiv_id":null,"evidence_quote":"Analyzes convergence and failure to converge of Gerchberg-Saxton, motivating the need for a convex initialization."},{"cited_title":"A high-accuracy algorithm for designing arbitrary holographic atom traps.Opt","cited_arxiv_id":null,"evidence_quote":"Defines the MRAF algorithm with signal and noise regions, the accuracy-efficiency trade-off that OT seeding improves."}],"review_version":1}