{"id":"209dfc7e-90f2-4d12-983b-795129db9c9f","arxiv_id":"2608.00550","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"For a photon redshifted in curved spacetime and subjected to amplitude-damping or Ohmic-like dephasing, the Nagaoka bound is the tightest multiparameter estimation limit, with best precision in the strong-coupling and sub-Ohmic regimes.","lead":"This paper calculates the best possible precision for measuring a photon's phase, amplitude, and gravitational redshift strength together, when the photon also loses information to noise. It finds the tightest error bound is the Nagaoka bound, and that precision is highest in the strong-coupling and sub-Ohmic noise regimes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (27) defines P_t as the amplitude rather than the squared amplitude, so all amplitude-damping error bounds in Sec. IV.A are based on an incorrect P_t that diverges in the Markovian limit.","rationale":"The most load-bearing concern is a concrete mathematical error: Eq. (27) missing the square on the amplitude factor, which invalidates the quantitative expressions for all amplitude-damping bounds if taken literally. The reader's flagged beam-splitter composition/phase issue is less severe: the loss and dephasing channels commute on this qubit subspace, and the beam-splitter phase φ(χ) cancels when the orthogonal mode is traced out, so neither affects the reduced density matrices (25) and (38). The P_t error is directly testable and fixable; it does not necessarily overturn the qualitative claim (bounds still increase with λ in both versions), so a conditional acceptance with a required correction is appropriate. The reader's verdict of CONDITIONAL therefore remains, but for a different reason than the one stated.","tokens_in":16691,"tokens_out":23518,"duration_ms":259067,"concrete_test":"Recompute P_t by integrating γ_t of Eq. (26): P_t = exp(-∫_0^t γ_{t'} dt'). At λ=2, γ0=1, t=1, the exact value is P_t≈0.54, while Eq. (27) gives ≈0.27. Recalculate C_S from Eq. (30) at α=θ=π/3 with both values; the discrepancy is roughly a factor of 2. Then check whether the curves in Fig. 2(a) match the squared or unsquared version. If the squared version is required to reproduce the figure, Eq. (27) is a typo and must be corrected.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. IV.A, Eq. (27) gives P_t = e^{-λt}[cosh(dt/2)+λ/d sinh(dt/2)]. For a two-level system in a Lorentzian reservoir, the excited-state amplitude is c(t)=e^{-λt/2}[cosh(dt/2)+λ/d sinh(dt/2)], so the population survival probability is |c(t)|², i.e., the square of the bracket. The paper uses P_t as the population in Eq. (25) and in all derived bounds (e.g., C_S = (1+csc²α)sec²θ/P_t). In the Markovian limit (λ≫γ0), the printed Eq. (27) gives P_t≈e^{-λt}e^{(λ-γ0)t/2}=e^{-(λ+γ0)t/2}→0 for fixed t, making every error bound diverge with λ, which is unphysical. The correct expression is P_t≈e^{-γ0 t}, yielding finite bounds. If this typo propagates into Eqs. (29)–(32) and Figs. 2–3, the central quantitative claim of enhanced precision in the strong-coupling regime is not supported as written.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies multiparameter quantum estimation for a photon wave packet subject to gravitational redshift, modeled as a beam-splitter mode mixing, followed by either amplitude-damping or Ohmic-like dephasing noise. It computes the SLD-CRB, RLD-CRB, Holevo CRB, and Nagaoka bound for two-parameter estimation (initial weight and phase), and adds the Nagaoka-Hayashi bound for three-parameter estimation (adding the redshift strength). Analytic expressions are given for the two-parameter qubit bounds using known formulas, while the three-parameter bounds are computed numerically via SDP. The central claims are that the Nagaoka bound is the tightest among the considered bounds and that estimation precision is enhanced in the strong-coupling (amplitude damping) and sub-Ohmic (dephasing) regimes.","tokens_in":17102,"tokens_out":12640,"duration_ms":127173,"significance":"If the quantitative results hold, the paper extends multiparameter estimation bounds to a curved-spacetime photonic setting and demonstrates regimes where the SLD-CRB is non-tight. The explicit two-parameter analytic formulas are a useful contribution, and the paper correctly identifies that the NB-tightest result is consistent with the known bound hierarchy. However, the main quantitative claim for the amplitude-damping channel is compromised by an amplitude/probability error in P_t, and the three-parameter numerical results are not reproducible from the text. The NB-tightest observation is largely a consistency check of Eq. (13) rather than an independent finding.","major_comments":[{"comment":"P_t is defined as an amplitude, not a probability, yet it is used as the excited-state population in Eq. (25) and in all amplitude-damping bounds (Eqs. (30)–(31), Figs. 2–3). For a qubit in a Lorentzian reservoir the excited-state amplitude is c(t)=e^{-λt/2}[cosh(dt/2)+λ/d sinh(dt/2)], so the survival probability is |c(t)|²=e^{-λt}[cosh(dt/2)+λ/d sinh(dt/2)]². Eq. (27) prints c(t). Consequently, in the Markovian limit λ≫γ0 the printed P_t≈e^{-(λ+γ0)t/2} vanishes, making C_S∝1/P_t diverge as e^{(λ+γ0)t/2}; the correct P_t≈e^{-γ0 t} is finite. In the strong-coupling regime (d imaginary) the printed P_t can even become negative, so ρ_red^A(t) in Eq. (25) is not positive semidefinite. This propagates into the central claim that strong coupling improves estimation precision. Please replace P_t by |c(t)|², recompute Eqs. (29)–(32) and the figures, and state whether the qualitative conclusions","section":"§IV.A, Eq. (27)"},{"comment":"The three-parameter results are stated to be computed 'using an SDP', but no SDP formulation, solver, tolerances, or data are given. Moreover, Eq. (11) is not a complete definition of the NHB: as printed it only imposes Hermiticity of L and L≥XX^T; it omits the unbiasedness constraints on X and the admissible-measurement conditions, so the displayed optimization is not the NHB. Without the full SDP and its implementation, the numerical comparisons underlying the three-parameter claims (strong vs weak coupling, sub-Ohmic vs Ohmic) cannot be checked. Please supply the complete SDP, code, and output data, or at least a fully specified mathematical definition and solver details.","section":"§IV.A/§IV.B, Figs. 3 and 5"},{"comment":"The frequency-ratio formula is stated as χ²=Ω_B/Ω_A=sqrt(f(r_B)/f(r_A)). For r_B>r_A, sqrt(f(r_B)/f(r_A))>1, so this predicts Ω_B>Ω_A, i.e., a blueshift, contradicting the sentence immediately above that 'the frequency Ω_B observed by Bob is lower than Ω_A' and the assertion that χ>1 corresponds to redshift. The standard relation is χ²=Ω_A/Ω_B=sqrt(f(r_B)/f(r_A)) (equivalently Ω_B/Ω_A=sqrt(f(r_A)/f(r_B))). Since χ enters Eq. (20) and hence the redshift parameter θ used throughout, the physical calibration of θ is affected. Correct the relation and re-examine the numerical value (χ−1)=3.5×10^{-10}.","section":"§III, Eq. (17)"}],"minor_comments":[{"comment":"The text says 'In Sec. V, we investigate...' but the multiparameter analysis is in Sec. IV; Sec. V is the conclusions. Please correct the cross-reference.","section":"Introduction"},{"comment":"The normalization denominator '4√(2πσ²)' appears to be a typo; a normalized Gaussian amplitude should use (2πσ²)^{1/4} or the equivalent normalization factor.","section":"Eq. (20)"},{"comment":"The reference list contains duplicates: [19] repeats [13], [24] repeats [15], and [53] repeats [56]. Please merge to avoid duplicate entries.","section":"References"},{"comment":"The statement that NB 'unequivocally' establishes the tightest bound is partly a consistency check with Eq. (13), since C_N≥C_H≥max(C_S,C_R) holds by construction. Please phrase this as consistency with the established hierarchy rather than an independent numerical discovery.","section":"§II and Figs. 2–5"},{"comment":"In Eq. (35) γ_t is a rate, while in q_t=e^{-γ_t/2} of Eq. (38) it must be the integrated decoherence factor. The notation is confusing; please distinguish the instantaneous decay rate from the integrated decoherence exponent.","section":"Eqs. (35) and (38)"}],"recommendation":"major_revision","confidential_remarks":"The central quantitative claim is currently not supported because of the P_t error in Eq. (27) and the undocumented three-parameter SDP. If the corrected P_t preserves the qualitative behavior, a revised version could be publishable, but the current text is not reproducible. The novelty over Refs. [16–18] is mostly incremental: the estimation-theory tools (Suzuki explicit formulas, SDP) are standard, and the NB-tightest result follows from the known hierarchy. The gravitational-redshift application is the main motivation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent application of standard multiparameter estimation bounds (SLD, RLD, HCRB, NB/NHB) to a photon qubit under gravitational redshift plus amplitude damping or Ohmic dephasing. The two-parameter analytic expressions are a genuine addition relative to the cited single-parameter studies, but there is a likely population/amplitude error in the amplitude-damping section that needs fixing before the quantitative results can be trusted.\n\nThe useful part: Eqs. (30), (31), (41), and (43) give closed-form SLD-CRB, RLD-CRB, HCRB, and NB for two-parameter estimation (weight and phase), which are absent from earlier redshift work. The checks that the SLD operators do not commute and the Uhlmann curvature is nonzero are consistent with the known non-tightness. The Ohmic dephasing section looks cleaner; the sub-Ohmic versus Ohmic/super-Ohmic trend is plausible and the formulas are internally consistent.\n\nThe problem: in Sec. IV.A, Eq. (27) defines P_t as e^{-λt}[cosh(dt/2)+λ/d sinh(dt/2)], which is the excited-state amplitude, not the population. For the density matrix in Eq. (25) to be correct, P_t must be the population, i.e., the square of that bracket. As printed, P_t ≈ e^{-(λ+γ0)t/2} in the Markovian limit λ≫γ0, which vanishes with λ for fixed t; the correct Markovian population is e^{-γ0 t}. Since all amplitude-damping bounds scale like 1/P_t, the curves in Figs. 2–3 and the claim about strong-coupling enhancement are quantitatively suspect. The qualitative ranking may survive, but the numbers are not supported as written.\n\nOther soft spots: the three-parameter numerical results (Figs. 3 and 5) depend on an unstated SDP implementation; no code or data is provided, so those numbers are not independently checkable. The conclusion that NB is the tightest bound is a consistency check of the known hierarchy (13), not a new result, and the abstract should frame it that way. The physical model also inherits the zero-phase beam-splitter approximation (φ(χ)=0) from earlier work; if that phase matters, the whole calculation shifts. That is a limitation worth stating explicitly.\n\nWho this is for: people interested in open-system multiparameter metrology applied to gravity. It is a solid case study, not a breakthrough. A referee would give useful feedback, and the amplitude-damping error should be caught in review. I would send it to review, but only with the expectation of a correction to Eq. (27) and a rerun of the affected results.","headline":"Workmanlike application of known multiparameter bounds to a gravitational-redshift qubit model; two-parameter analytic formulas are new, but Eq. (27) appears to use the amplitude instead of the population, which undermines Sec. IV.A as written.","tokens_in":17561,"tokens_out":5454,"would_cite":false,"duration_ms":62550,"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":"For a photon distorted by gravitational redshift and then damped or dephased, the tightest attainable two-parameter error bound is the Nagaoka bound, and precision improves in strong-coupling and sub-Ohmic regimes.","keywords":["multiparameter quantum estimation","gravitational redshift","Nagaoka bound","Holevo Cramér–Rao bound","amplitude-damping channel","Ohmic-like dephasing","semidefinite program","quantum Fisher information"],"falsifier":"Recompute or experimentally realize the two-parameter bounds with the relative phase kept at its physical value instead of zero: if the Nagaoka bound drops below the analytic values in Eq. (31) or (43), or if the bound ordering of Eq. (13) is violated, the central tightness claim fails. A tabletop beam-splitter plus engineered noise can serve as the test.","tokens_in":16655,"feed_emoji":"🛰️","tokens_out":9775,"duration_ms":105490,"temperature":0.7,"pith_summary":"This paper asks what the fundamental precision limit is when one simultaneously estimates several properties—the weight, phase, and gravitational redshift strength—of a photon that has climbed out of a gravitational potential and then been disturbed by noise. It argues that for two-parameter estimation the standard quantum Cramér–Rao bound is not attainable, and that the Nagaoka bound is the tightest achievable error bound for single-copy measurements. The same conclusion is carried over to three-parameter estimation, where the Nagaoka–Hayashi bound plays that role. The paper also reports that, in both noise models studied, precision is better in the memory-retaining regimes: strong coupling for amplitude damping and sub-Ohmic spectra for dephasing. If true, this matters for quantum metrology in curved spacetime, because it identifies which error bound a realistic receiver should be compared against.","feed_headline":"Nagaoka bound sets tightest precision limit for redshifted photons","feed_subtitle":"For a redshifted photon, single-copy precision is capped by the Nagaoka bound, not the quantum Fisher bound.","key_machinery":"The load-bearing object is the Nagaoka bound (NB), a lower bound on the mean-square error matrix optimized over single-copy measurements; for two-parameter qubit systems it is known to be tight, unlike the quantum Cramér–Rao bounds. The argument feeds the redshifted, noise-affected density matrices (25) and (38) through the bound hierarchy (13) and uses the semidefinite-program formulations of the Holevo, Nagaoka, and Nagaoka–Hayashi bounds for numerical evaluation. The beam-splitter overlap cosθ of the redshifted wave packet is the parameter that turns spacetime curvature into a quantum channel.","core_discovery":"The paper considers a photonic qubit initially prepared in cos(α/2)|0⟩ + e^{iφ} sin(α/2)|1⟩, sent through Schwarzschild spacetime from Alice to Bob. Gravitational redshift is represented as a beam-splitter mixing of the signal mode with an orthogonal mode, with overlap cosθ determined by the wave-packet deformation and with relative phase set to zero. After tracing out the orthogonal mode, the photon evolves under either amplitude damping or Ohmic-like dephasing, yielding the two-qubit density matrices (25) and (38). For the two parameters α and φ, the authors find the SLD operators do not commute and the mean Uhlmann curvature is nonzero, so the quantum Cramér–Rao bound is not tight even as","pith_inferences":["Because the Nagaoka bound is attainable with single-copy measurements, a laboratory experiment that simulates gravitational redshift with a beam splitter and adds engineered amplitude damping or dephasing can test the predicted hierarchy and regime ordering without a space link; the quantum-optics part is the load-bearing physics, not the spacetime propagation.","The zero-phase choice for the beam splitter is a simplification. If the physical relative phase is nonzero, the compatibility of α and φ may change, so the tightness ranking should be re-checked before relying on these bounds for satellite-based quantum communication.","The improved precision in memory-retaining regimes suggests the same tools could be applied to other gravity-induced noise models, for example gravity acting as a universal dephasing channel for qubits, but those extensions are not established by this paper.","The numerical equality of the RLD-CRB and HCRB in the three-parameter case is reported as a numerical observation; a general proof of that equivalence would be a useful follow-up."],"forward_implications":["The SLD-CRB cannot be used as the precision benchmark for simultaneous weight-and-phase estimation of a redshifted single photon; the Nagaoka bound is the correct single-copy benchmark and is larger.","Strong-coupling (non-Markovian) amplitude damping improves the attainable precision relative to weak coupling, so memory effects can be treated as a resource for multiparameter estimation.","Sub-Ohmic dephasing environments consistently give tighter attainable bounds than Ohmic or super-Ohmic environments, for both two- and three-parameter estimation.","In three-parameter estimation, the RLD-CRB and HCRB coincide, meaning both are asymptotically attainable, while the Nagaoka–Hayashi bound remains the tightest among the single-copy bounds considered.","The known hierarchy of multiparameter bounds—most informative bound, Nagaoka, Holevo, then SLD/RLD—holds quantitatively in this gravitational-redshift setting."],"supporting_citations":[{"why":"Supplies the two reduced density matrices (25) and (38) that combine gravitational redshift with each noise channel; all estimation bounds are computed on these states.","marker":"[18]"},{"why":"Introduces the gravitational-redshift beam-splitter model for phase estimation that this paper extends to multiparameter estimation.","marker":"[16]"},{"why":"Provides the beam-splitter treatment of gravitational redshift used to set up the photon state before the noise channels.","marker":"[17]"},{"why":"Used to justify setting the beam-splitter relative phase to zero, a premise on which the computed density matrices rest.","marker":"[25]"},{"why":"Defines the Nagaoka bound and establishes its tightness for two-parameter qubit estimation, the paper's central claim.","marker":"[45]"},{"why":"Supplies the explicit two-parameter qubit Holevo bound formula used to write analytic HCRB and NB expressions in Eqs. (31) and (43).","marker":"[67]"},{"why":"Defines the Nagaoka–Hayashi bound via a semidefinite program for separable measurements, used in the three-parameter numerical comparison.","marker":"[44]"},{"why":"Gives the semidefinite-program approach used to evaluate the Holevo Cramér–Rao bound numerically.","marker":"[42]"},{"why":"Provides the general multiparameter bound hierarchy and the compatibility conditions showing when the SLD-CRB is not tight.","marker":"[28]"}],"fun_headline_variants":["Nagaoka bound provides tightest multiparameter precision for redshifted photons","Cramer-Rao bound not tight for multiparameter redshifted photon estimation","Multiparameter precision for redshifted photons: Nagaoka bound is tightest","Gravitational redshift: Nagaoka bound caps multiparameter precision"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The whole calculation assumes gravitational redshift can be modeled as a beam-splitter mode-mixing with zero relative phase, and that the redshift reduction followed by the noise channel is the correct order; if either fails, the computed bounds are for a different physical situation.","fun_headline_variants_meta":{"raw":{"variants":["Nagaoka bound provides tightest multiparameter precision for redshifted photons","Cramer-Rao bound not tight for multiparameter redshifted photon estimation","Multiparameter precision for redshifted photons: Nagaoka bound is tightest","Gravitational redshift: Nagaoka bound caps multiparameter precision"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001142,"raw_usage":{"total_tokens":4613,"prompt_tokens":820,"completion_tokens":3793,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":3709}},"tokens_in":564,"tokens_out":3793,"duration_ms":29834,"temperature":1.0,"reasoning_tokens":3709,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T00:40:35.108577+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute or experimentally realize the two-parameter bounds with the relative phase kept at its physical value instead of zero: if the Nagaoka bound drops below the analytic values in Eq. (31) or (43), or if the bound ordering of Eq. (13) is violated, the central tightness claim fails. A tabletop beam-splitter plus engineered noise can serve as the test.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two reduced density matrices (25) and (38) that combine gravitational redshift with each noise channel; all estimation bounds are computed on these states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the beam-splitter treatment of gravitational redshift used to set up the photon state before the noise channels."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Used to justify setting the beam-splitter relative phase to zero, a premise on which the computed density matrices rest."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Nagaoka bound and establishes its tightness for two-parameter qubit estimation, the paper's central claim."},{"cited_title":"Bradshaw, P","cited_arxiv_id":null,"evidence_quote":"Gives the semidefinite-program approach used to evaluate the Holevo Cramér–Rao bound numerically."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the general multiparameter bound hierarchy and the compatibility conditions showing when the SLD-CRB is not tight."}],"review_version":1}