{"id":"f36c1272-3e91-42e5-9db2-bed8a06f7235","arxiv_id":"2606.07257","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"Gerchberg-Saxton phase retrieval is shown to be exactly unit gradient descent on amplitude least-squares loss, with global and local probabilistic interpretations.","lead":"This paper shows that the classic Gerchberg-Saxton magnitude replacement step in phase retrieval is mathematically identical to one unit step of gradient descent on an amplitude least-squares loss. The finding lets older iterative methods plug directly into modern automatic-differentiation pipelines used in computational imaging.","discovery_kind":"unification","skeptic_critique":{"model":"grok-4.3","headline":"No significant objection identified","rationale":"The reader’s concern was limited to the abstract’s lack of explicit conditions. The full text supplies the derivation; the identity holds exactly under the stated loss and step size, so the abstract-level caveat does not affect the correctness of the central claim.","tokens_in":1639,"tokens_out":279,"duration_ms":12761,"concrete_test":"Recompute the gradient of L = ½∑(|z_i| − a_i)² in Wirtinger calculus and verify that one step with η = 1 produces exactly the phase-preserving magnitude projection; the algebra is a few lines and settles the claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the GS magnitude replacement is exactly the unit-step GD update on the amplitude least-squares loss L = ½‖|z| − a‖². Direct computation of the (Wirtinger) gradient yields ∇L = (|z| − a)(z/|z|), so the update z − ∇L collapses algebraically to a·(z/|z|). This identity is unconditional once the loss and step size are fixed; the paper’s title correctly flags the “conditional” qualifier. No hidden assumptions about the propagation operator or noise model are required for the replacement step itself.","agreement_with_reader":"disagree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript establishes that the magnitude replacement step in the Gerchberg-Saxton (GS) algorithm is mathematically identical to a single unit-step gradient descent update on the amplitude least-squares loss L = ½‖|z| − a‖². Direct Wirtinger differentiation yields ∇L = (|z| − a)(z/|z|), so the update z − ∇L algebraically simplifies to a·(z/|z|). The paper additionally supplies two probabilistic readings: the amplitude loss as negative log-likelihood under Gaussian amplitude noise, and each projection as a local Bayesian update with the propagated field as prior. These observations are presented as enabling integration of classical phase retrieval into differentiable physics pipelines.","tokens_in":1731,"tokens_out":358,"duration_ms":12161,"significance":"The result supplies a parameter-free, direct algebraic bridge between a classical iterative method and modern gradient-based optimization, allowing GS steps to be dropped unchanged into autodiff frameworks. The probabilistic interpretations supply qualitative guidance for step-size relaxation without introducing new parameters or data fits. Because the identity follows immediately once the loss and step size are fixed, the contribution is self-contained and falsifiable by direct substitution.","major_comments":[],"minor_comments":[{"comment":"The abstract states the equivalence but does not explicitly name the Wirtinger gradient or the precise loss function; adding one sentence with the expression for ∇L would improve immediate readability for readers outside phase retrieval.","section":null},{"comment":"Notation for the complex field z and measured amplitude a is introduced without a dedicated symbols table; a short table in §2 would eliminate any ambiguity when the same symbols appear in the probabilistic sections.","section":null}],"recommendation":"accept","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their positive review and recommendation to accept the manuscript.","responses":[],"tokens_in":1174,"tokens_out":34,"duration_ms":5671,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The central result is the algebraic identity between the GS magnitude replacement and a unit step of gradient descent on L = 1/2 || |z| - a ||^2. Direct computation of the gradient gives exactly the replacement operation, so the equivalence is exact under those conditions.\n\nThis is new in the sense that it makes the connection explicit and adds the probabilistic readings. The global view as negative log-likelihood under Gaussian amplitude noise is standard, but tying the local step to a Bayesian update is a nice way to think about relaxation.\n\nThe paper does this cleanly. The math is simple and checks out, as the stress test shows. No circularity or fitting involved.\n\nSoft spots: the equivalence is conditional on the specific loss and unit step size, but that's flagged in the title and not a flaw. The local probabilistic interpretation is more heuristic than rigorous, but presented as such. Without the full text I can't check if there are examples or edge cases explored, but the core claim stands.\n\nThis is for researchers in phase retrieval and differentiable imaging who want to combine the two toolkits. A reader working on hybrid algorithms would get immediate value from the identity.\n\nIt deserves peer review because the observation is precise and the implications are practical for the field.","headline":"The GS magnitude replacement is exactly unit GD on amplitude loss, a direct algebraic identity that bridges classical and differentiable phase retrieval.","tokens_in":2211,"tokens_out":328,"would_cite":true,"duration_ms":13746,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"The Gerchberg-Saxton magnitude replacement step is exactly a unit gradient descent update on the amplitude least-squares loss.","keywords":["phase retrieval","Gerchberg-Saxton","gradient descent","amplitude least-squares","differentiable imaging","computational optics","Bayesian update"],"falsifier":"Direct numerical evaluation of the gradient of the amplitude least-squares loss followed by a unit step, compared against the field produced by magnitude replacement under identical conditions, would show whether the two operations coincide.","tokens_in":2509,"feed_emoji":"","tokens_out":639,"duration_ms":12851,"temperature":0.7,"pith_summary":"The paper shows that the magnitude replacement operation in the Gerchberg-Saxton algorithm matches a single gradient descent step of size one when the objective is the amplitude least-squares loss. This identity is derived by direct comparison of the two update rules under the stated conditions. A reader cares because it removes the apparent separation between classical projection-based phase retrieval and modern gradient-based methods in computational imaging. The same identity also yields two probabilistic readings: the loss as negative log-likelihood under Gaussian amplitude noise, and each step as a Bayesian update with the propagated field as prior.","feed_headline":"Gerchberg-Saxton step equals unit gradient descent on amplitude loss","feed_subtitle":"The identity lets classical projection methods plug directly into gradient-based optimization pipelines for imaging.","key_machinery":"The conditional mathematical identity between the Gerchberg-Saxton magnitude replacement operation and a unit-step gradient descent update on the amplitude least-squares loss.","core_discovery":"The GS magnitude replacement step is exactly a unit gradient descent step on an amplitude least-squares loss. This equivalence enables seamless integration of classical phase retrieval with differentiable physics pipelines. Two complementary probabilistic interpretations follow: globally, the amplitude loss is the negative log-likelihood under Gaussian amplitude noise; locally, each projection step arises as a Bayesian update with the propagated field as prior. The local view supplies qualitative guidance for relaxation in iterative phase retrieval.","pith_inferences":["The equivalence suggests that any phase retrieval variant using a different loss could be rewritten as a modified gradient step, opening a route to new algorithms.","Implementation in machine-learning frameworks becomes direct: the magnitude replacement can be replaced by an autograd call with fixed step size.","The Bayesian-update view may extend to other noise models, yielding analogous projection steps for Poisson or other intensity statistics."],"forward_implications":["Classical phase retrieval methods become native components of automatic-differentiation pipelines without custom implementation.","The amplitude least-squares loss functions as the canonical objective that unifies projection and gradient approaches.","Probabilistic readings allow relaxation parameters in iterative algorithms to be interpreted as prior strengths.","Hybrid algorithms can alternate or combine projection steps with gradient steps while preserving the same loss surface."],"fun_headline_variants":["GS step equals unit gradient descent on amplitude loss","Gerchberg-Saxton matches amplitude loss gradient step","Phase retrieval projections are gradient descent steps","GS magnitude replacement is unit GD on amplitude loss"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The identity holds only when the objective is exactly the amplitude least-squares loss and the descent step size is precisely one.","fun_headline_variants_meta":{"raw":{"variants":["GS step equals unit gradient descent on amplitude loss","Gerchberg-Saxton matches amplitude loss gradient step","Phase retrieval projections are gradient descent steps","GS magnitude replacement is unit GD on amplitude loss"]},"model":"grok-4.3","cost_usd":0.002902,"raw_usage":{"total_tokens":1556,"prompt_tokens":575,"num_sources_used":0,"completion_tokens":56,"cost_in_usd_ticks":29024500,"prompt_tokens_details":{"text_tokens":575,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":925,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":575,"tokens_out":56,"duration_ms":8170,"temperature":1.0,"reasoning_tokens":925,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-27T20:59:27.720725+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Direct numerical evaluation of the gradient of the amplitude least-squares loss followed by a unit step, compared against the field produced by magnitude replacement under identical conditions, would show whether the two operations coincide.","supporting_citations":[],"review_version":1}