{"id":"2851f979-d091-494c-96a1-3ee287cd2c05","arxiv_id":"2501.06768","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"For coherent light, a saturating photodetector's average current is I=I_max(1-exp(-N/Ñ_sat)); inverting this relation recovers the optical phase in the nonlinear response regime.","lead":"This paper calculates how photodetector saturation changes optical phase measurements in a standard homodyne setup, and proposes a correction that inverts the detector's nonlinear response. The work is a short theoretical note relevant to very high light intensities where detectors saturate.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (8) omits the Poissonian variance of mu_j(n): the law of total variance gives Var(k_j) = sigma_j^2 + Var[mu_j(n)], so the central precision formula is mathematically unsupported and predicts zero noise for ideal detectors.","rationale":"The reader's overall rejection is correct, but I locate the load-bearing problem more precisely in the variance calculation than in the imported detector model. Equation (8) is one of the paper's two headline results, and its derivation contains a textbook error: for a mixture distribution, the total variance includes the variance of the conditional mean. The paper writes delta^2 k_j = sigma_j^2, which would be correct only if mu_j(n) were constant. For a coherent state, mu_j(n) fluctuates with n, and this fluctuation is exactly the shot noise responsible for the standard quantum limit in homodyne detection. Ignoring it makes Eq. (8) underestimate the uncertainty and, in the zero-detector-noise limit, predict perfect phase sensitivity, which is impossible for a coherent state. This is an internal inconsistency, not merely a disagreement with a particular calibration model. The reader's stated weakest assumption -- precise knowledge of Eq. (3) -- is a legitimate external concern, but it is secondary: even with perfect knowledge of the saturation curve, Eq. (8) is still wrong. A quick analytical check at sigma = 0 settles the matter, so the rejection should stand.","tokens_in":5727,"tokens_out":5310,"duration_ms":58221,"concrete_test":"Set sigma = 0 and work in the linear regime N_j << N_sat. Apply the law of total variance to Eq. (4): Var(k_j) = (k_max/N_sat)^2 N_j. Substitute this into the error-propagation formula (7) with F(I) = I/r and I_j = r N_j. The corrected phase error is finite and scales as ~1/sqrt(N), while Eq. (8) gives exactly zero at sigma = 0. Recompute the right-hand side of Eq. (7) both ways at the paper's Table I parameters (for example, N/N_sat = 1, |beta|^2 = 10^15); if the corrected value differs from Eq. (8) by more than a few percent, the precision claim is falsified.","verdict_should_be":"REJECT","load_bearing_attack":"The decisive flaw is in the derivation of Eq. (8). The paper claims that from the mixture distribution ~P_j(k) in Eq. (4) one obtains delta^2 k_j = sigma_j^2. This is false: the law of total variance gives Var(k_j) = E_n[Var(k_j|n)] + Var_n(E[k_j|n]) = sigma_j^2 + Var_n[mu_j(n)]. Because P_j(n) is Poisson with mean N_j, the second term is nonzero; in the linear regime N_j << N_sat it equals (k_max/N_sat)^2 N_j. Dropping this term removes exactly the photon-number (shot-noise) contribution that produces the standard quantum limit. The problem is not a small missing correction: for an ideal detector with sigma = 0, Eq. (8) predicts zero phase uncertainty for a coherent state, whereas balanced homodyne detection has finite phase uncertainty scaling as 1/sqrt(N). Thus Eq. (8) cannot serve as the paper's precision result, and the subsequent claim delta_chi ~ 1/sqrt(MN) is not established. The sensitivity of Eq. (8) to the imported saturation model, as noted by the reader, is real but secondary; even perfect knowledge of Eq. (3) would not repair the variance error.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies optical phase estimation via homodyne detection when photodetectors exhibit saturation. It models the mean photoelectron count with μ_j(n)=k_max(1-e^{-n/N_sat}), derives the average current I_j=I_max(1-e^{-N_j/\\tilde N_sat}) in Eq. (5), and introduces the inverse function F so that F(I_2)-F(I_1)=N_2-N_1, enabling phase extraction in the nonlinear regime via Eq. (6). The authors then use error propagation to derive a phase-uncertainty formula Eq. (8) and claim it scales as 1/√(MN) with photon number and measurement repetitions. They illustrate the bias of the standard linear protocol when the detector saturates and argue that the inverse-function correction improves the phase estimate.","tokens_in":6000,"tokens_out":15242,"duration_ms":133658,"significance":"The mean-current calculation leading to Eq. (5) is straightforward and correct, and the idea of undoing the saturation nonlinearity by applying the inverse of the assumed response to the measured currents is sensible. However, the paper's precision analysis is invalid: the variance of k_j is not σ_j², as claimed, because the photon-number fluctuations contribute through the law of total variance. This missing term is the usual shot noise, and its omission produces an unphysical formula that would give zero phase uncertainty for an ideal detector. The claimed 1/√(MN) scaling is therefore not established. The conceptual contribution of the bias correction is limited to the exact validity of the saturation model, and the quantitative results need substantial revision.","major_comments":[{"comment":"The variance calculation leading to Eq. (8) is incorrect. From the mixture distribution Eq. (4), the law of total variance gives Var(k_j) = E_n[Var(k_j|n)] + Var_n(E[k_j|n]) = σ_j² + Var_n[μ_j(n)]. In the linear regime N_j ≪ N_sat, this second term equals (k_max/N_sat)^2 N_j, which is precisely the photon-number shot noise responsible for the standard quantum limit. The text's claim that δ²k_j = σ_j² omits this term; consequently Eq. (8) predicts zero phase uncertainty for a noiseless detector (σ=0) and does not reproduce the SQL. Eq. (8) and the subsequent scaling claim δχ ∼ 1/√(MN) are therefore mathematically unsupported.","section":"Eq. (8) and the variance derivation following Eq. (7)"},{"comment":"The inversion formula Eq. (6) is definitional: F is constructed as the inverse of the assumed mean-current relation Eq. (5), so F(I_2)-F(I_1)=N_2-N_1 holds by construction. The method's practical utility depends on the saturation model Eq. (3) being exact and on knowing N_sat and k_max precisely. The manuscript does not analyze how parameter uncertainty or model misspecification propagates into the phase estimate, nor does it compare the model with measured detector response data. Without such an analysis, the claimed 'improved estimation' in the saturation regime is not quantitatively supported beyond the idealized model.","section":"Eq. (6) and the saturation model of Eq. (3)"},{"comment":"Even accepting the derivation of Eq. (8), the claimed SQL scaling is not obtained: in the linear response regime N_j ≪ N_sat, Eq. (8) simplifies to δχ ∝ σ/(|αβ|), which depends on the local-oscillator amplitude |β| and can be made arbitrarily small by increasing the LO power while keeping σ fixed. This is not the standard quantum limit, which for balanced homodyne detection should be independent of the LO amplitude once shot noise is included. The omitted term Var(μ_j(n)) is required to cancel this |β| dependence and recover the genuine 1/√N scaling.","section":"Final paragraph and the scaling claim after Eq. (8)"}],"minor_comments":[{"comment":"The equality ⟨μ_j(n)⟩ = μ_j(⟨n⟩) is stated without qualification; it is only approximate, because ⟨μ_j(n)⟩ = k_max(1 - e^{-N_j/\\tilde N_sat}) while μ_j(⟨n⟩) = k_max(1 - e^{-N_j/N_sat}), and the two expressions coincide only in the limit N_sat → ∞. Please state this explicitly.","section":"After Eq. (5)"},{"comment":"The notation for the saturation parameter is confusing: the paper uses N_sat, \\tilde N_sat, and the textual abbreviation ~N_sat at different points. Since \\tilde N_sat = (1-e^{-1/N_sat})^{-1} differs from N_sat, the text should clearly indicate whenever the approximation \\tilde N_sat ≈ N_sat is being used.","section":"Equation (5) and the following text"},{"comment":"Reference [4] contains a typo: 'adn' should be 'and'.","section":"References"},{"comment":"In Fig. 2, the estimator for χ~ is not defined; the authors should specify how the estimated phase is obtained from the measured currents in both the linear and nonlinear protocols.","section":"Figure 2"},{"comment":"The phrase 'This relation is not obvious at all in priori' should read 'a priori'.","section":"After Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The principal quantitative result is invalid as stated, but the underlying idea of inverting the saturation model to correct the bias is coherent and could be repaired with a correct variance calculation and a sensitivity analysis. The current manuscript is not acceptable, and the authors should be required to re-derive the precision formula and re-evaluate the scaling claims before a revised version is considered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the paper's average-current formula for saturating detectors, Eq. (5), is correct: Poisson-averaging the saturation model from Harris et al. gives I = Imax(1 - exp(-N/Nsat)), and inverting that to recover the phase bias is valid if the model and its parameters are exact. Second, the precision result, Eq. (8), is wrong. The paper claims the photoelectron variance is just sigma_j^2, but the law of total variance gives Var(k_j) = sigma_j^2 + Var_n[mu_j(n)], and the second term is the shot-noise contribution from Poisson photon statistics. Dropping it makes Eq. (8) predict zero phase uncertainty for an ideal detector, which contradicts balanced homodyne detection's 1/sqrt(N) limit. This isn't a minor correction; the whole point of Eq. (8) is to quantify the precision after saturation. So the paper's central new claim fails.\n\nWhat is actually good: the derivation of Eq. (5) is clean and the inversion idea in Eq. (6) is a legitimate way to correct bias when the response is known. The numerical example in Table I illustrates the saturation effect on the average current clearly. The paper is readable and the self-reported lack of comparison with real detector response curves is acknowledged implicitly by the model assumption.\n\nSoft spots in proportion: the variance error is the load-bearing flaw, and it is fatal for Eq. (8). The secondary issue is that the whole inversion strategy assumes exact knowledge of N_sat and I_max; the paper doesn't quantify sensitivity or compare with measured curves. That would matter even if the variance were fixed. The false variance equality is not a matter of interpretation; it's an elementary probability error.\n\nWho is this for? Someone working on high-intensity homodyne phase estimation who wants to know whether saturation breaks the standard protocol. The bias-correction part could be useful if the noise were handled correctly, but as written the paper doesn't deliver a consistent precision statement. I wouldn't cite it. For peer review: I would desk reject it in this form; the error is decisive and the remaining content is not substantial enough to warrant referee time. If the authors correct Eq. (8) by including the shot-noise term, the bias-correction part could be a short note.","headline":"Clean bias-correction idea, but the precision formula drops the shot-noise variance and is unsupportable.","tokens_in":6564,"tokens_out":3929,"would_cite":false,"duration_ms":37574,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Optical phase can still be estimated from saturated photodetectors by inverting each detector's saturation curve before taking the current difference.","keywords":["optical phase estimation","homodyne measurement","photodetector saturation","nonlinear detector response","coherent states","standard quantum limit","inverse response function","postselection"],"falsifier":"Take a calibrated photodetector, send coherent pulses of known mean photon number $N$ spanning from well below to well above $N_{\\rm sat}$, and record the average current; the paper's Eq. (5) predicts $I=I_{\\max}(1-e^{-N/\\tilde N_{\\rm sat}})$, so any systematic deviation from that curve at fixed $N$ would falsify the saturation model that Eq. (6) relies on.","tokens_in":5505,"feed_emoji":"🔬","tokens_out":13544,"duration_ms":104179,"temperature":0.7,"pith_summary":"The paper extends optical phase estimation by homodyne measurement from the usual linear detector response to the nonlinear, saturated regime. For coherent light it derives the average detector current $I_j=I_{\\max}(1-e^{-N_j/N_{\\rm sat}})$ and shows that applying the inverse function $F$ to the two detector currents gives $F(I_2)-F(I_1)=2|\\alpha\\beta|\\sin(\\chi-\\phi)$, so the unknown phase $\\chi$ can still be extracted when the detectors are saturated. The accompanying precision formula scales as $1/\\sqrt{MN}$, i.e., the standard quantum limit. This matters because proposed high-intensity laser interferometers for relativistic-gravity effects push photon numbers into the saturation territory of photodetectors, where the standard linear protocol would give a wrong phase.","feed_headline":"Inverting a saturated detector's response restores phase readout","feed_subtitle":"At saturating light intensities, inverting the detector response restores homodyne phase readout at the standard quantum limit.","key_machinery":"The load-bearing object is the saturating exponential response $\\mu_j(n)=k_{\\max}(1-e^{-n/N_{\\rm sat}})$ for the average number of photoelectrons produced by $n$ incident photons. Its key property is that averaging over the Poisson photon-number distribution of a coherent field preserves the exponential form: $\\langle\\mu_j(n)\\rangle=\\mu_j(\\langle n\\rangle)$, yielding the current $I_j=I_{\\max}(1-e^{-N_j/\\tilde N_{\\rm sat}})$. The inverse function $F(I_j)$ then converts measured currents back into mean photon numbers, reducing the saturated homodyne problem to the same difference identity $N_2-N_1=2|\\alpha\\beta|\\sin(\\chi-\\phi)$ that is used in the linear regime.","core_discovery":"The central discovery is that the saturation nonlinearity can be removed by a deterministic inversion. Starting from the model $\\mu_j(n)=k_{\\max}^{(j)}(1-e^{-n/N_{\\rm sat}^{(j)}})$ for the average photoelectron number and averaging it over the Poisson photon statistics of a coherent state, the paper obtains the same exponential form for the current, $I_j=I_{\\max}^{(j)}(1-e^{-N_j/\\tilde N_{\\rm sat}^{(j)}})$, with $\\tilde N_{\\rm sat}^{(j)}\\simeq N_{\\rm sat}^{(j)}$. Defining $F$ as the inverse of this current--photon-number relation, the phase is read from $F(I_2)-F(I_1)=2|\\alpha\\beta|\\sin(\\chi-\\phi)$. The precision of this readout, including shot-to-shot detector fluctuations and $M$ repeated measurements, is $\\delta\\chi\\sim 1/\\sqrt{MN}$---the same standard quantum limit as the linear-regime measurement---so the generalization does not cost precision, provided the saturation curve is known and the detector is not oversaturated.","pith_inferences":["Beyond the paper, the same inversion strategy should apply to any monotone detector nonlinearity, not only the exponential saturation model; the essential requirement is that the current remains an invertible function of the mean photon number.","The paper does not quantify how errors in $N_{\\rm sat}$ or $k_{\\max}$ propagate into the phase estimate; a natural extension is to compute the sensitivity $\\partial\\chi/\\partial N_{\\rm sat}$ and design a calibration procedure.","A testable prediction is that, for a fixed phase $\\chi$, the phase extracted with $F$ should stay flat as the input intensity is swept through saturation, whereas the linear protocol shows a growing bias; this could be checked with an independently calibrated detector.","For non-coherent states the simplification $\\langle\\mu_j(n)\\rangle=\\mu_j(\\langle n\\rangle)$ generally fails, so extending the method to squeezed or thermal light would require a different averaging; the present result is specific to coherent states."],"forward_implications":["For mean photon numbers up to a few times $N_{\\rm sat}$, the inverse-response protocol $F(I_2)-F(I_1)$ recovers the phase $\\chi$ that the standard linear protocol would systematically misestimate.","The precision of the saturated-regime estimate remains $\\sim 1/\\sqrt{MN}$, the standard quantum limit, rather than degrading with saturation, as long as the response curve and its parameters are known.","Once both detectors are fully saturated, the currents carry no phase information and no inversion can recover $\\chi$; the paper identifies this oversaturation regime as an in-principle limitation.","To operate high-intensity lasers without entering oversaturation, the paper points to postselection as a way to keep photodetectors below the saturation threshold.","The derivation relies on coherent-state Poisson statistics, so the simple inversion works for the coherent signal states used throughout the paper."],"supporting_citations":[{"why":"Supplies the saturation-response model $\\mu_j(n)=k_{\\max}(1-e^{-n/N_{\\rm sat}})$ and the experimental context for detector saturation in precision measurement.","marker":"[8]"},{"why":"Supplies the same saturation model and demonstrates saturation-limited metrology, providing the empirical basis for Eq. (3).","marker":"[9]"},{"why":"Motivates the high-photon-number regime with a nonlinear laser interferometer proposal whose photon numbers enter the saturation region.","marker":"[5]"},{"why":"Another relativistic-gravity detection proposal requiring the high-intensity regime where saturation is expected.","marker":"[6]"},{"why":"Extends wormhole parameter estimation beyond the Heisenberg limit, setting the context in which very large photon numbers are required.","marker":"[7]"},{"why":"The authors' recent postselection analysis used to argue that postselection can keep detectors below the saturation threshold even at high intensity.","marker":"[13]"},{"why":"Defines the standard quantum limit $1/\\sqrt{N}$ against which the precision scaling of Eq. (8) is identified.","marker":"[4]"}],"fun_headline_variants":["Saturating detectors still hit standard quantum limit with inversion","Invert detector saturation to recover homodyne phase precision","Nonlinear detector response: invert to keep quantum-limited phase","Phase estimation survives detector saturation via inversion","Detector saturation: inversion restores quantum-limited readout"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands or falls on the assumption that a photodetector's average photoelectron count follows exactly the saturating exponential $\\mu_j(n)=k_{\\max}(1-e^{-n/N_{\\rm sat}})$ and that $N_{\\rm sat}$ and $k_{\\max}$ are known precisely when the measured current is inverted.","fun_headline_variants_meta":{"raw":{"variants":["Saturating detectors still hit standard quantum limit with inversion","Invert detector saturation to recover homodyne phase precision","Nonlinear detector response: invert to keep quantum-limited phase","Phase estimation survives detector saturation via inversion","Detector saturation: inversion restores quantum-limited readout"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000278,"raw_usage":{"total_tokens":1590,"prompt_tokens":818,"completion_tokens":772,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":434,"completion_tokens_details":{"reasoning_tokens":695}},"tokens_in":434,"tokens_out":772,"duration_ms":6042,"temperature":1.0,"reasoning_tokens":695,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:50:30.592435+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a calibrated photodetector, send coherent pulses of known mean photon number $N$ spanning from well below to well above $N_{\\rm sat}$, and record the average current; the paper's Eq. (5) predicts $I=I_{\\max}(1-e^{-N/\\tilde N_{\\rm sat}})$, so any systematic deviation from that curve at fixed $N$ would falsify the saturation model that Eq. (6) relies on.","supporting_citations":[{"cited_title":"Giovannetti, S","cited_arxiv_id":null,"evidence_quote":"Supplies the saturation-response model $\\mu_j(n)=k_{\\max}(1-e^{-n/N_{\\rm sat}})$ and the experimental context for detector saturation in precision measurement."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the same saturation model and demonstrates saturation-limited metrology, providing the empirical basis for Eq. (3)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Motivates the high-photon-number regime with a nonlinear laser interferometer proposal whose photon numbers enter the saturation region."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Another relativistic-gravity detection proposal requiring the high-intensity regime where saturation is expected."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The authors' recent postselection analysis used to argue that postselection can keep detectors below the saturation threshold even at high intensity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the standard quantum limit $1/\\sqrt{N}$ against which the precision scaling of Eq. (8) is identified."}],"review_version":1}