{"id":"66363aed-69e1-4833-b859-f143559f3651","arxiv_id":"2506.23600","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"A linear-system construction for the symmetric logarithmic derivative is derived from the master equation and applied to temperature and relaxation-rate metrology in a squeezed Gaussian state.","lead":"The authors show how to build the optimal measurement operator for estimating environment properties directly from a quantum master equation, without first solving for the full quantum state. The method is checked on equilibrium thermometry and then applied to a non-equilibrium squeezed state to estimate temperature and relaxation rate.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Appendix A drops the ∂tΛ terms in Eq. (41) without justification beyond 'excluding a time dependence'; yet the c_i(t) plotted in Figs. 1–4 are time-dependent, so the transient construction in Eq. (10) does not follow from the master equation.","rationale":"The reader's weakest assumption identifies the same load-bearing defect: the derivation of Eq. (10) relies on dropping the time derivative of the SLD in Eq. (41), even though the paper's own transient coefficients c_i(t) are time-dependent. I find no additional fatal flaw beyond this, and I confirm that it is decisive. In the steady-state validation of Sec. 3 the assumption ∂tΛ = 0 is satisfied, and the reproduction of the known equilibrium result is genuine evidence that the algebraic machinery works there. The deposited Mathematica scripts also support reproducibility. However, the abstract and Secs. 2 and 4 claim validity for non-equilibrium and transient regimes, and there the missing ∂tΛ terms are not negligible by construction; the paper never quantifies the error or compares against an exact transient solution. Therefore the central novelty is unproven as written, and the reader's rejection is appropriate. The concrete test proposed above would settle the matter: for the Gaussian case an exact transient SLD is computable, and for a simple two-level channel it is fully analytic. If the exact result agrees with Fig. 1, the concern is refuted; if not, the central claim fails.","tokens_in":13805,"tokens_out":5391,"duration_ms":63918,"concrete_test":"For the same Caldeira–Leggett master equation used in Sec. 4, compute the exact SLD for the squeezed Gaussian state by two independent routes: (i) solve the covariance-matrix equations for the Gaussian state to obtain explicit ρ(t,T), then evaluate the SLD from the standard Gaussian-state formula or by numerically solving the Lyapunov equation; (ii) re-derive Eq. (41) while keeping the ∂tΛ terms, which converts Eq. (10) into a differential equation for c_i(t), and integrate it for the parameters of Fig. 1 (T/m = 4, γ/m = 1/8, ω/m = 4/3 and 2/3). If either comparison deviates from the algebraic coefficients plotted in Fig. 1 at early times, the dropped ∂tΛ terms materially change the result and the transient claim fails; if all agree, the omission is harmless.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the SLD for an environmental parameter can be constructed directly from the master equation via the algebraic system Eq. (10). The derivation in Appendix A differentiates the SLD definition Eq. (1) with respect to time. The correct identity contains the terms (1/2)(∂tΛ ρ + ρ ∂tΛ) and an additional contribution Tr[ρ ∂tΛ] inside ∂t⟨Λ⟩. Equation (41) instead contains neither, with the sole justification 'excluding a time dependence of the SLD operator.' This assumption fails in the transient regime that is the paper's stated novelty: the coefficients c_i(t) shown in Figs. 1–4 are time-dependent while the basis operators A_i are time-independent, so the operator Λ(t) must itself be time-dependent. The algebraic system Eq. (10) is therefore at best an uncontrolled approximation in the non-equilibrium setting, and no error estimate or independent benchmark is provided for that regime. The equilibrium validation in Sec. 3 is consistent because ∂tΛ = 0 in the steady state, but it cannot license the transient construction. Thus the claim that optimal observables are obtained exactly in the transient regime without solving the master equation is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a method to construct the symmetric logarithmic derivative (SLD) for an environmental parameter θ directly from a master equation of an open quantum system, without solving for the full density matrix. The technique expands the SLD in a fixed operator basis, uses the equality of mixed partial derivatives (Schwarz's theorem) to relate ∂θ∂tρ and ∂t∂θρ, and reduces the problem to a linear system for the expansion coefficients, Eq. (10). The method is tested in two applications: the equilibrium steady state of Caldeira–Leggett quantum Brownian motion, where the known temperature SLD of Ref. [13] is reproduced, and a non-equilibrium squeezed Gaussian state, for which SLDs for temperature and for the relaxation rate are constructed using time-dependent expectation values from the Ehrenfest equations. The paper concludes that its approach enables optimal-observable construction in transient and non-equilibrium regimes.","tokens_in":14034,"tokens_out":11832,"duration_ms":121017,"significance":"If the proposed construction were valid, it would be a practical tool for quantum metrology in open systems, where SLDs are usually obtained from explicit density matrices. The equilibrium validation in Sec. 3 is a useful consistency check and the published Mathematica scripts (Ref. [24]) support reproducibility. However, the central novelty—the transient and non-equilibrium construction—rests on a derivation that explicitly assumes a time-independent SLD, which is inconsistent with the time-dependent coefficients shown in the paper's own figures. The non-equilibrium results are not benchmarked against any independent computation, and the assertion that the Caldeira–Leggett master equation has no closed-form solution in that regime is incorrect. The main claim is therefore not established.","major_comments":[{"comment":"The derivation of the linear system (10) drops the terms (1/2)(∂tΛ ρ + ρ ∂tΛ) in Eq. (41), with the explicit justification 'excluding a time dependence of the SLD operator'. In the non-equilibrium application, however, the coefficients c_i(t) shown in Figs. 1–4 are time-dependent while the basis operators A_i are time-independent, so Λ(t) must itself be time-dependent. The linear system (10) therefore does not follow from the master equation in the transient regime, and the paper's primary claim of constructing optimal observables out of equilibrium is unsupported. The equilibrium validation in Sec. 3 is consistent because the steady-state SLD is time-independent, but it cannot license the transient construction, for which no error estimate or independent benchmark is provided.","section":"Section 2 and Appendix A, Eq. (41)"},{"comment":"The paper states that 'there does not exist a closed form explicit solution of the density matrix to the CL master equation in this regime', but the Caldeira–Leggett master equation for a harmonic oscillator with Gaussian initial states admits analytic solutions, e.g., the Wigner-function solutions in Refs. [23] and [16], which the paper itself cites. This makes the claimed novelty of the non-equilibrium application questionable and, more importantly, means that an independent benchmark for the transient SLD is available but has not been used. Without such a benchmark, the coefficients in Figs. 1–3 and the QFI in Fig. 4 are uncontrolled.","section":"Section 4, beginning of Sec. 4"},{"comment":"The claim that the approach 'does not rely on the explicit solution of the master equation' is overstated. In Sec. 4, the matrices M and D are evaluated using time-dependent expectation values obtained from the Ehrenfest equations (35)–(37). For the harmonic Caldeira–Leggett model these equations close exactly, but for a general master equation the required moment hierarchy will generally not close, and the paper does not discuss how the necessary expectation values are to be obtained in such cases. The generality claimed in the introduction and Sec. 2 is therefore not supported beyond the specific Gaussian examples treated.","section":"Section 2 (method's scope) and Section 4 (Ehrenfest equations)"}],"minor_comments":[{"comment":"The expression '⟨p^4⟩ = 3⟨p^2⟩' is missing the exponent on the right-hand side and should read '⟨p^4⟩ = 3⟨p^2⟩^2'.","section":"Eq. (19)"},{"comment":"There is a typo: 'valubale' should be 'valuable'.","section":"Introduction, last paragraph"},{"comment":"The phrase 'They key relation' should be 'The key relation'.","section":"Section 3, first paragraph"},{"comment":"The caption refers to 'blue solid' and 'orange solid' lines, but the text preceding the figure describes 'green solid' and 'purple solid' for the same ratios; the colors should be made consistent.","section":"Figure 3 caption"},{"comment":"The theorem is usually attributed as 'Schwarz's theorem' (with one 't'), while the paper repeatedly writes 'Schwartz'; this should be corrected.","section":"Throughout"},{"comment":"The susceptibility χ_T(ˆx^2) is written twice in the text; the second instance should be χ_T(ˆp^2).","section":"Eq. (3)"}],"recommendation":"reject","confidential_remarks":"The equilibrium part of the paper is a clean derivation of a known result and the code deposit is a positive feature. The main problem is that the transient non-equilibrium construction, which is the advertised novelty, is not justified by the derivation: the explicit time-independence assumption in Appendix A is violated by the time-dependent coefficients plotted in the application. This is not a local issue that can be fixed with a minor revision; the algebraic system (10) would need to be replaced by a differential equation for the SLD, which is a substantial reformulation. I also note the incorrect statement about the absence of closed-form solutions for the Caldeira–Leggett master equation, which weakens the novelty claim and leaves the non-equilibrium results without the independent check that is actually available."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe interesting part of this paper is the first step: using Schwarz's theorem to trade ∂θ∂tρ for ∂t∂θρ and then substituting the SLD definition to get a linear system for the SLD coefficients. That is a neat bridge, and I have not seen it in the cited SLD constructions. The equilibrium validation in Sec. 3 is also careful—under γ ≪ 2πT they reproduce the known temperature SLD for quantum Brownian motion, and the deposited Mathematica scripts back the algebra. Full credit for that.\n\nThe problem is the transient claim, which is the paper's stated novelty. In Appendix A, Eq. (41) is written without the (1/2)(∂tΛ ρ + ρ ∂tΛ) terms, with the justification 'excluding a time dependence of the SLD operator.' That is an assumption, not a derivation. It happens to be valid in the steady-state limit, where expectation values are constant and Λ is fixed, but in the very transient regime they advertise, the coefficients c_i(t) plotted in Figs. 1–4 are explicitly time-dependent. So the assumption is violated precisely where the new results live. The internal inconsistency is hard to miss: the text says 'the SLD is independent of time' while the figures show its expansion coefficients varying with time.\n\nThe consequence is that the central result, Eq. (10), is at best an uncontrolled approximation out of equilibrium. No error estimate or independent benchmark is provided for the transient regime, and the equilibrium test cannot license it. The abstract's claim of not relying on an explicit solution is also weaker than it looks: the coefficients require evolving expectation values, which for Gaussian states means solving a closed hierarchy of moments—useful, but not the same as avoiding the master equation entirely.\n\nWho is this for? People working on open-system metrology will find the equilibrium construction and the moment-hierarchy approach worth a look, and the Schwarz trick is a genuinely reusable idea. But as it stands, the non-equilibrium results are not established. I would send it to review—the idea is worth refereeing and the flaw is fixable—but the referee should be asked to either extend the derivation to include the ∂tΛ terms or restrict the claims to steady states and clearly label the transient construction as an approximation with a validity estimate.\n\nRegards,","headline":"The Schwarz-mixed-derivative bridge to the SLD is a genuinely new technique and the equilibrium validation is clean, but the transient construction drops the time derivative of the SLD and is not established.","tokens_in":14585,"tokens_out":5315,"would_cite":false,"duration_ms":56864,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81P50","81S22"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper constructs the symmetric logarithmic derivative directly from the master equation, via the linear system $\\sum_i M_{ji} c_i = D_j$, without first solving for the density matrix.","keywords":["quantum metrology","symmetric logarithmic derivative","quantum Fisher information","master equation","open quantum systems","Caldeira-Leggett model","quantum Brownian motion","non-equilibrium thermometry"],"falsifier":"Numerically solve the Caldeira-Leggett master equation for the squeezed Gaussian initial state used in Section 4, reconstruct $\\hat{\\rho}(t,T)$ on a fine grid, solve Eq. (1) for the exact time-dependent $\\hat{\\Lambda}_T(t)$, and expand it in the same five-operator basis. If the resulting coefficients disagree with the $c_T^{(i)}(t)$ shown in Fig. 1 in the transient regime beyond numerical error, the dropped $\\partial_t \\hat{\\Lambda}$ terms are significant; if they agree, the time-independent-SLD ansatz is validated.","tokens_in":13566,"feed_emoji":"🎯","tokens_out":10730,"duration_ms":96424,"temperature":0.7,"pith_summary":"The paper aims to make optimal-observable construction available for open quantum systems that have not reached equilibrium. Its central claim is that the symmetric logarithmic derivative, the operator whose measurement minimizes estimation error for an environmental parameter $\\theta$, can be computed from the master equation alone: Schwarz's theorem allows the parameter derivative and the time derivative to be exchanged, and an operator-basis expansion turns the SLD's defining equation into a finite linear system $\\sum_i M_{ji} c_i = D_j$. Explicit knowledge of the density matrix $\\hat{\\rho}(t,\\theta)$ is not required. If correct, this opens transient and strongly coupled regimes where the density matrix is only available numerically. The paper validates the construction by reproducing the known equilibrium temperature SLD in the Caldeira-Leggett model and then builds time-dependent SLDs for temperature and for the relaxation rate in a squeezed Gaussian state.","feed_headline":"No density matrix needed: master equation yields optimal quantum sensors","feed_subtitle":"Construction gives the optimal observable straight from the generator, covering transient and non-equilibrium metrology.","key_machinery":"The key machinery is the mixed-partial-derivative identity (Schwarz's theorem), used as a bridge between dynamics and parameter sensitivity, combined with a fixed operator-basis expansion $\\hat{\\Lambda}_\\theta = \\sum_i c_i \\hat{A}_i$. Inserting the expansion into the $\\theta$-derivative of the master equation and tracing against test operators $\\hat{A}_j$ produces the linear system of Eq. (10); the Caldeira-Leggett generator, Weyl-symbol correspondences (an operator-ordering bookkeeping device), and Wick factorization for Gaussian states supply the expectation values that fill $M$ and $D$.","core_discovery":"The discovery is a constructive route to the SLD that bypasses the density matrix. Starting from a linear master equation $\\partial_t \\hat{\\rho} = \\mathcal{L}[\\theta]\\hat{\\rho}$ and the mixed-derivative identity $\\partial_\\theta \\partial_t \\hat{\\rho} = \\partial_t \\partial_\\theta \\hat{\\rho}$, the paper rewrites the implicit Lyapunov equation for $\\hat{\\Lambda}_\\theta$ as a linear system $\\sum_i M_{ji} c_i = D_j$. The matrix $M$ is built from expectation values of operator bilinears generated by $\\mathcal{L}$, and the source vector $D$ from the explicit $\\theta$-dependence of the generator; solving it gives the coefficients of $\\hat{\\Lambda}_\\theta = \\sum_i c_i \\hat{A}_i$. In the Caldeira-Leggett model the method reproduces the known equilibrium temperature SLD, and in a squeezed Gaussian state it yields time-dependent SLDs for both temperature and relaxation rate, with the temperature Fisher information growing toward its late-time maximum while the relaxation-rate Fisher information decays to zero.","pith_inferences":["Extension: because the method needs only expectation values of a chosen operator basis, it can be paired with any numerical solver that propagates those values, which would extend it beyond Gaussian states to lattice or tensor-network simulations.","Extension: in a multi-parameter setting the same matrix $M$ can be reused with different source vectors $D_j$, suggesting a practical route to joint estimation of temperature and relaxation rate, although multi-parameter bounds are not derived in the paper.","Extension: a direct numerical computation of the exact time-dependent SLD for the same squeezed state would quantify how much the time-independence assumption used in Appendix A affects the transient coefficients."],"forward_implications":["The symmetric logarithmic derivative becomes computable for transient, driven, or strongly coupled probes from the master-equation generator plus evolving expectation values, without solving for $\\hat{\\rho}(t,\\theta)$.","In the Caldeira-Leggett steady state the construction reproduces the known temperature SLD and shows that the anticommutator $\\{\\hat{x},\\hat{p}\\}$ must be retained in the basis even where its expectation value vanishes.","For a squeezed Gaussian state, the temperature-SLD coefficients evolve in time and asymptote to their equilibrium values, meaning the optimal temperature estimator is a time-dependent weighted combination of position and momentum variances.","The relaxation-rate SLD is intrinsically non-equilibrium: its quantum Fisher information decays to zero at late times, while the temperature quantum Fisher information grows to its steady-state maximum.","The same linear system works for any environmental parameter appearing in a linear master equation, with different parameters entering through different source vectors $D_j$."],"supporting_citations":[{"why":"Supplies the known equilibrium SLD for a heavy bosonic impurity in terms of $\\hat{x}^2$ and $\\hat{p}^2$ coefficients that the paper reproduces in the late-time Caldeira-Leggett limit.","marker":"[13]"},{"why":"Provides the non-equilibrium dephasing-impurity thermometry example whose SLD requires derivatives of the explicit density-matrix solution, motivating the derivative-free construction.","marker":"[14]"},{"why":"Provides the open-quantum-systems formalism and the Ehrenfest equations used to evolve the expectation values in the Caldeira-Leggett examples.","marker":"[7]"},{"why":"Gives the analytic treatment of the Caldeira-Leggett master equation used to argue twice differentiability and to characterize transient density-matrix behavior.","marker":"[16]"},{"why":"Source of the Caldeira-Leggett master equation used as the working model in both the equilibrium validation and the non-equilibrium applications.","marker":"[19, 20]"},{"why":"Supplies the Weyl-symbol correspondences used to evaluate the operator-bilinear expectation values entering the matrix $M$.","marker":"[21]"},{"why":"Provides Wick's theorem, which factors higher moments of the Gaussian state and closes the expectation values needed for $M$ and $D$.","marker":"[22]"},{"why":"Contains the supplementary coefficient expressions and scripts that make the non-equilibrium SLD results reproducible.","marker":"[24]"}],"fun_headline_variants":["Master equation yields optimal quantum observables without solving it","Skip the density matrix: master equation builds optimal sensors","Optimal metrology from master equation, equilibrium or not","Direct master equation route to optimal quantum sensitivity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The construction assumes the optimal-observable operator itself does not change with time when the parameter derivatives are taken, even though the paper applies it to transient regimes where the operator's expansion coefficients are time dependent; if the SLD genuinely varies during the evolution, the central linear system misses terms.","fun_headline_variants_meta":{"raw":{"variants":["Master equation yields optimal quantum observables without solving it","Skip the density matrix: master equation builds optimal sensors","Optimal metrology from master equation, equilibrium or not","Direct master equation route to optimal quantum sensitivity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000192,"raw_usage":{"total_tokens":1304,"prompt_tokens":861,"completion_tokens":443,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":477,"completion_tokens_details":{"reasoning_tokens":381}},"tokens_in":477,"tokens_out":443,"duration_ms":4696,"temperature":1.0,"reasoning_tokens":381,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:37:08.910546+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically solve the Caldeira-Leggett master equation for the squeezed Gaussian initial state used in Section 4, reconstruct $\\hat{\\rho}(t,T)$ on a fine grid, solve Eq. (1) for the exact time-dependent $\\hat{\\Lambda}_T(t)$, and expand it in the same five-operator basis. If the resulting coefficients disagree with the $c_T^{(i)}(t)$ shown in Fig. 1 in the transient regime beyond numerical error, the dropped $\\partial_t \\hat{\\Lambda}$ terms are significant; if they agree, the time-independent-SLD ansatz is validated.","supporting_citations":[{"cited_title":"Using Polarons for sub-nK Quantum Nondemolition Thermometry in a Bose-Einstein Condensate","cited_arxiv_id":null,"evidence_quote":"Supplies the known equilibrium SLD for a heavy bosonic impurity in terms of $\\hat{x}^2$ and $\\hat{p}^2$ coefficients that the paper reproduces in the late-time Caldeira-Leggett limit."},{"cited_title":"In Situ Ther- mometry of a Cold Fermi Gas via Dephas- ing Impurities","cited_arxiv_id":null,"evidence_quote":"Provides the non-equilibrium dephasing-impurity thermometry example whose SLD requires derivatives of the explicit density-matrix solution, motivating the derivative-free construction."},{"cited_title":"The theory of open quantum sys- tems","cited_arxiv_id":null,"evidence_quote":"Provides the open-quantum-systems formalism and the Ehrenfest equations used to evolve the expectation values in the Caldeira-Leggett examples."},{"cited_title":"Positivity violations of the density operator in the caldeira-leggett mas- ter equation","cited_arxiv_id":null,"evidence_quote":"Gives the analytic treatment of the Caldeira-Leggett master equation used to argue twice differentiability and to characterize transient density-matrix behavior."},{"cited_title":"On the Function in Quantum Mechanics which Corresponds to a Given Function in Classical Mechanics","cited_arxiv_id":null,"evidence_quote":"Supplies the Weyl-symbol correspondences used to evaluate the operator-bilinear expectation values entering the matrix $M$."},{"cited_title":"The evaluation of the collision matrix","cited_arxiv_id":null,"evidence_quote":"Provides Wick's theorem, which factors higher moments of the Gaussian state and closes the expectation values needed for $M$ and $D$."},{"cited_title":"url: https:// doi.org/10.5281/zenodo.15770098","cited_arxiv_id":null,"evidence_quote":"Contains the supplementary coefficient expressions and scripts that make the non-equilibrium SLD results reproducible."}],"review_version":1}