{"id":"c4ee07e6-7c34-4b44-968d-8abf94976266","arxiv_id":"2607.15160","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Lossy Fock boson sampling probabilities are finite Fourier series in each phase, so their exact gradients can be recovered from shifted photon-count measurements; Gaussian boson sampling under general loss admits no such finite rule.","lead":"This paper derives exact parameter-shift rules — a way to compute gradients from a handful of shifted measurements — for lossy Fock boson sampling circuits, and shows the standard recipe cannot work for Gaussian boson sampling under general loss. Practical use: gradient-based tuning of photonic quantum chips without the noise sensitivity of finite differences.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Exactness of the Fock PSR rests on Eq. (31)'s theta-independent transmission sqrt(eta); the paper's physical justification is an order-of-magnitude argument, not a device-level guarantee.","rationale":"The central positive claim—exact order-n PSR for Fock boson sampling under arbitrary loss—is mathematically sound conditional on the phase-shifter model of Eq. (31). The permanent/hafnian derivation is internally consistent, and the frequency-support argument for [-n,n] can be verified by block decomposing the permanent. The most load-bearing external assumption is that the phase-shifter amplitude transmission is θ-independent and that each θ appears exactly once. The paper's justification is explicitly heuristic (Sec. II.B), not a measured device property. If real phase shifters exhibit amplitude-phase coupling, the exactness of the PSR on hardware fails, and the experimental validation in Sec. IV would not detect this because the comparison is against the same lossy model. The GBS no-go is also somewhat overclaimed, but it is a secondary negative result; the positive Fock claim and the hardware demonstration are more central. The reader identified the same weakest assumption, and my concrete test targets exactly that gap. Since the paper states the assumption explicitly and the math follows from it, the appropriate verdict remains conditional rather than rejection.","tokens_in":20562,"tokens_out":11090,"duration_ms":118350,"concrete_test":"Characterize the actual phase shifters on the Belenos QPU (or an identical test structure) by sending coherent light through a single phase shifter and measuring transmitted power |Φ(θ)|^2 as θ is stepped in fine increments over [0,2π). If |Φ(θ)|^2 deviates from a constant by more than the shot-noise-limited uncertainty of the gradient measurements, or if its discrete Fourier spectrum has components at frequencies |m|>0 above noise, Eq. (31) is violated and the order-n PSR is not exact on that device. A complementary numerical check—replacing Eq. (31) by Φ(θ)=√η(θ)e^{iθ} with η(θ)=η0(1+ε cos2θ) and recomputing the four-mode example of Fig. 3—would quantify the resulting gradient error, but the hardware measurement is the decisive test.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The Fock PSR result, Eq. (66), follows from Eq. (64), which assumes every tunable phase enters T as a e^{iθ}+b with θ-independent a,b. This in turn relies on Eq. (31): Φ_i^(j)=√η_{i,j} e^{iθ_{i,j}}, i.e., a phase shifter whose energy transmission is independent of the applied phase, and with no crosstalk so each θ appears once. The paper's support for this (Sec. II.B) is that the optical path modification is at most of order the wavelength, so transmission is 'not significantly impacted'; this is an order-of-magnitude plausibility statement, not a characterization of the actual Belenos/Perceval phase shifters. In realistic thermo-optic or electro-optic phase shifters, η can vary with the phase setpoint (e.g., free-carrier absorption, mode-overlap changes, heater-induced index/loss coupling). If η=η(θ), the entries of T become sums of terms a(θ)e^{iθ}+b(θ); the permanent in Eq. (66) is no longer guaranteed to be a finite Fourier series of degree n, and the n-th order PSR acquires a model error set by the Fourier bandwidth of η(θ). This does not contradict the mathematics under the stated assumption, but it means the headline claim—exact gradients under arbitrary loss on a real device—is not yet established for the hardware class used.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives parameter-shift rules (PSRs) for gradient estimation in photonic boson sampling circuits with loss. For Fock-state inputs, it claims that each transition probability is a finite Fourier series in any tunable phase with frequencies bounded by the number of input photons n, so an n-th order PSR gives the exact derivative even under arbitrary loss. For Gaussian boson sampling (GBS), it claims that no finite-order PSR exists under general loss; a finite-order rule of order 2|j| exists only when the loss commutes through the interferometer (T = UΔ). The paper also presents numerical simulations and a hardware demonstration on Quandela's Belenos QPU comparing PSR with finite differences.","tokens_in":20877,"tokens_out":8685,"duration_ms":99870,"significance":"If the Fock-state result is correct, it is a useful and nontrivial generalization of PSRs: it removes the need for post-selection or unitary-only models when computing gradients on lossy photonic hardware. The derivation is largely self-contained, and the permanent/hafnian formulas are cross-checked against known results. The GBS no-go result, if made rigorous, would be an important limitation for variational methods with squeezed states. The hardware demonstration is a valuable practical check, although it does not by itself validate the exactness of the Fock PSR under the most general loss model.","major_comments":[{"comment":"The claim that no finite-order PSR exists for GBS under general loss is not rigorously established. The two supporting arguments — that det(Λ'_uu)^(-1/2) is not polynomial in the entries of T, and that Eq. (97) contains an infinite sum over photon numbers — are suggestive but not proofs. A non-polynomial dependence on T does not logically exclude cancellation that leaves a finite Fourier series in θ, and Eq. (97) is only an example (loss after the unitary), not a general impossibility. Since the abstract advertises this as a main result, please either supply a concrete counterexample (e.g., a simple two-mode circuit with generic loss where the probability has infinite Fourier support) or soften the claim to 'the standard PSR construction does not apply in general.'","section":"Sec. III.B, Eqs. (94)–(97)"},{"comment":"The exactness of the Fock PSR rests on the assumption Φ_i^(j) = √η_{i,j} e^{iθ_{i,j}}, i.e., that each phase shifter's amplitude transmission is independent of the applied phase and that each θ appears exactly once. If η depends on θ (thermo-optic or electro-optic phase shifters can exhibit loss-phase coupling), then T entries no longer have the form a e^{iθ} + b, and the frequency bound in Eq. (66) can fail. The paper's physical justification is an order-of-magnitude statement about optical path length, not a device-level guarantee. Please either provide phase-shifter characterization data from Belenos/Perceval showing θ-independent transmission over the relevant range, or explicitly scope the headline claim to the model in Eq. (31).","section":"Sec. II.B, Eq. (31)"},{"comment":"The derivation of the frequency support [-n, n] is only sketched as 'three observations.' For a reader, the bound is not immediate because the permanent in Eq. (61) has size |i|+|j|, not |i|. I recommend explicitly stating the block-counting argument: in B^{i⊕j}, with m=|j|, every permutation must contain exactly m T-factors and m T†-factors and (|i|-m) E-factors; since T-factors and E-factors contribute at most +1 frequency and T†-factors contribute at most 0, the maximum positive frequency is m + (|i|-m) = |i|. This would make the main result easier to verify.","section":"Sec. III.A, Eq. (66)"}],"minor_comments":[{"comment":"The prefactor is written as '1/⃗j!' but the state is |⃗k⟩; it should be '1/⃗k!'.","section":"Eq. (94)"},{"comment":"The vertical-axis labels are garbled ('| Pr( )|'); please fix the mathematical notation in the captions.","section":"Fig. 3 and Fig. 4 captions"},{"comment":"The phrase 'without loss of generality' before the assumption that phase shifters do not affect transmission is too strong; it is an assumption of the model, not a consequence of generality.","section":"Sec. II.B"},{"comment":"The ordering of indices in the linear system is unusual but acceptable; please add a sentence clarifying that the DFT choice μ_m = 2πm/(2n+1) yields an invertible system for all n, not just for the displayed ordering.","section":"Sec. II.C after Eq. (38)"},{"comment":"The hardware experiment uses a two-mode HOM setup with uniform loss; this is a valid proof-of-principle but does not exercise the 'arbitrary loss' scenario of the main theorem. A sentence acknowledging this limitation would be appropriate.","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The central Fock PSR result appears mathematically sound under the stated model, and the paper is a good candidate for publication after revision. The GBS no-go claim, however, is currently argued rather than proven, and the abstract overstates it. The phase-shifter loss-independence assumption is also load-bearing for the 'exact on real hardware' language. Both issues are fixable within the manuscript's scope, so I do not recommend rejection, but they need to be addressed before the paper can be accepted. Note also that Ref. [42] shares an author with the present work; the manuscript does cite it properly, and the derivations here are self-contained, so I do not see a novelty-disclosure problem."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Core result is real: for Fock boson sampling under arbitrary loss, each transition probability is a finite Fourier series in any given phase, with order n=|i| (Eq. 66), so an n-th order parameter-shift rule is exact. I checked the duplicated-permanent decomposition (Eqs. 59-61) and the frequency-support argument; it holds, and the derivation is self-contained, agreeing with Ref. [19] without the Choi-Jamiolkowski route. The GBS boundary result — no finite-order rule under general loss, with a 2|j|-th order rule for T=UΔ — is new relative to the cited literature, including Ref. [42] which shares an author. Good work.\n\nSoft spots, in order of importance. First, the Fock PSR is exact only under Eq. (31): phase-shifter transmission η independent of θ, no crosstalk. The Sec. II.B justification is an order-of-magnitude argument about optical path length, not a device-level guarantee. On thermo-optic or electro-optic shifters, η can vary with setpoint; then T entries become a(θ)e^{iθ}+b(θ), and the frequency bound fails. This does not contradict the mathematics under the stated assumption, but the abstract's \"exact gradients under arbitrary loss on a real device\" is stronger than what is established. Second, the GBS no-go is argued, not proven: the non-polynomial dependence of det^{-1/2} and the hafnian, plus the infinite-support example, constitute strong evidence, but not an impossibility proof for all lossy GBS circuits. The abstract states it categorically; the body is more careful. Third, the novelty boundary vs. Refs. [42,43] is never explicitly drawn; \"going further and beyond\" should be itemized. The experiments are a useful sanity check — PSR beats FD in simulation and on a small two-mode Belenos run — but single-run, no code/data release, and cannot validate the main theorem.\n\nWho this is for: researchers doing gradient-based tuning and variational training of lossy photonic devices with Fock inputs. The GBS boundary is of independent interest to theory groups. The paper deserves a serious referee; the core is sound and the new result is worth publishing. The revision is about precision of claims, not existence of the result. Send to peer review.","headline":"The Fock PSR result is real and the GBS boundary is new, but the paper overclaims on hardware exactness and the GBS no-go; worth refereeing with revision requests on scope and claims.","tokens_in":21504,"tokens_out":2959,"would_cite":true,"duration_ms":31875,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In Fock boson sampling under loss, transition probabilities are finite Fourier series in each phase, enabling exact gradients via n-th order parameter-shift rules; general GBS admits no such rule.","keywords":["parameter-shift rules","Fock boson sampling","Gaussian boson sampling","photon loss","gradient estimation","permanent","hafnian","variational quantum algorithms"],"falsifier":"Measure the complex transmission matrix T(θ) of a real interferometer as a function of one phase (coherent-state input and heterodyne detection). If any entry deviates from a e^{iθ}+b — for example, amplitude varies with θ — the n-th order shift rule will show a systematic bias that scales with that coupling. Alternatively, for the GBS claim, measure a squeezed-state output probability versus θ under general loss and fit a finite Fourier series; if a finite order fits exactly, the claimed impossibility is wrong.","tokens_in":20382,"feed_emoji":"⚛️","tokens_out":7282,"duration_ms":72519,"temperature":0.7,"pith_summary":"Fock boson sampling on a lossy interferometer retains a hidden structure: every transition probability is a finite Fourier series in any one tunable phase, with frequencies bounded by the number of input photons. This paper shows that structure makes the exact gradient computable by an n-th order parameter-shift rule — evaluating the circuit at shifted phases and linearly combining the results — even when the device is lossy and noisy. For Gaussian boson sampling the same strategy generally fails, because squeezed states have unbounded photon-number support and loss mixes in infinitely many frequency terms; an exact finite-order rule survives only when loss can be pushed before the interferometer. The authors verify the rule on simulated and real hardware, where it is markedly more stable than finite differences. This matters for variational photonic algorithms, which need reliable gradients to optimize circuits on the device itself.","feed_headline":"Finite shift rule yields exact gradients in lossy Fock boson sampling","feed_subtitle":"Transition probabilities are finite Fourier series in each phase, so gradients stay exact under loss and phase noise.","key_machinery":"The load-bearing identity is the permanent representation of the lossy Fock transition probability, (1/(i! j!)) Perm[B^{i⊕j}], together with the fact that each entry of the transmission matrix is a single-harmonic function of any given phase, T_mn = a e^{iθ}+b. From this, the entries of the permanent matrix B take the form a' e^{iθ}+b' e^{−iθ}+c', and because the permanent is a polynomial of degree equal to the total photon number n, the whole probability becomes a finite Fourier series with frequencies confined to [−n,n]. The parameter-shift rule exploits this structure by differentiating the series and solving a small linear system (Eqs. 38–43) to get exact derivative coefficients from shi","core_discovery":"For a lossy interferometer with transmission matrix T, when each tunable phase θ enters entries as T_mn = a e^{iθ}+b with a,b independent of θ, the Fock-state transition probability from input i to output j equals (1/(i! j!)) Perm[B^{i⊕j}], where B is constructed from T and E = I − T†T. Because a permanent of size n is a polynomial in its entries, the probability is a finite Fourier series Σ_{m=−n}^{n} k_m e^{imθ}, and its derivative can be reconstructed exactly from shifted circuit evaluations via an n-th order parameter-shift rule. The analogous Gaussian boson sampling probability is a hafnian expression divided by the square root of a determinant (Eqs. 90–94); under general loss these fac","pith_inferences":["A testable extension suggested by the argument: the same finite-Fourier logic should apply to any photonic observable that is a degree-n polynomial in single-harmonic matrix entries, such as photon-number moments, multi-mode correlation functions, or averages of observables under Fock inputs.","The GBS no-go points to a practical truncation strategy: an approximate finite-order rule could be built by truncating the infinite Fock support of squeezed states, with an error that decays as the tail of the squeezed-state distribution; the paper does not explore this.","Because exactness rests on amplitude–phase decoupling, a hardware-level characterization of T(θ) along the lines of coherent-state tomography would provide a direct engineering criterion for when the method is applicable; the paper itself offers only an order-of-magnitude argument.","The permanent formula is derived without the Choi–Jamiołkowski trick, so the same derivation route might generalize to partial distinguishability or mixed input states, though the paper does not address those cases."],"forward_implications":["On any device where the phase-shifter model holds, the gradient of a Fock-state transition probability is computable exactly from 2n+1 (or n with symmetric shifts) circuit evaluations, independent of the number of modes or the loss level.","Variational photonic algorithms can replace finite differences with this shift rule, removing the step-size trade-off and gaining robustness to phase noise, as demonstrated numerically and on Quandela's QPU.","For Gaussian boson sampling, the result is a no-go: under general loss no exact finite-order rule exists; users must either engineer loss before the interferometer (where a rule of order 2|j| applies) or accept approximate alternative gradient estimators.","Threshold-detector probabilities inherit the same finite Fourier structure in the Fock case, since they are finite sums of Fock probabilities, so the rule extends to click patterns.","The order bound n = number of input photons gives a direct cost estimate: gradient evaluation scales linearly with input photon number, not with Hilbert-space dimension."],"fun_headline_variants":["Shift rules give exact gradients in lossy Fock boson sampling","Shift rules make Fock boson gradients exact under loss","Parameter-shift rules for exact gradients in Fock boson sampling","Exact gradients via parameter-shift rules in Fock boson sampling","Finite-order shift rules yield exact boson sampling gradients"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The rule is exact only if each tunable phase appears exactly once in the circuit and its amplitude transmission is independent of the phase value, so that T entries are a e^{iθ}+b; if real phase-shifters couple amplitude to phase, the finite Fourier bound and the n-th order rule break down.","fun_headline_variants_meta":{"raw":{"variants":["Shift rules give exact gradients in lossy Fock boson sampling","Shift rules make Fock boson gradients exact under loss","Parameter-shift rules for exact gradients in Fock boson sampling","Exact gradients via parameter-shift rules in Fock boson sampling","Finite-order shift rules yield exact boson sampling gradients"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001735,"raw_usage":{"total_tokens":6660,"prompt_tokens":675,"completion_tokens":5985,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":419,"completion_tokens_details":{"reasoning_tokens":5897}},"tokens_in":419,"tokens_out":5985,"duration_ms":43571,"temperature":1.0,"reasoning_tokens":5897,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T00:01:26.662134+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the complex transmission matrix T(θ) of a real interferometer as a function of one phase (coherent-state input and heterodyne detection). If any entry deviates from a e^{iθ}+b — for example, amplitude varies with θ — the n-th order shift rule will show a systematic bias that scales with that coupling. Alternatively, for the GBS claim, measure a squeezed-state output probability versus θ under general loss and fit a finite Fourier series; if a finite order fits exactly, the claimed impossibility is wrong.","supporting_citations":[],"review_version":1}