{"id":"309f80bd-931a-46cf-ac32-084f4038b8fa","arxiv_id":"2608.04539","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"For 1/|ω| spectra, any passive Lindblad network with N modes has the same optimal approximation error as N independent damped modes, given by the degree-2N Zolotarev bound.","lead":"This paper proves that adding coherent coupling between auxiliary modes does not improve how accurately a finite set of damped modes can reproduce a 1/f -type noise spectrum. The result gives exact limits on how many modes are needed to approximate such spectra to a given tolerance, which is useful for quantum simulation and noise modeling.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the central equality is well supported; the only load-bearing external input is the Zolotarev parity/type lemma cited to Ref. [27], which the explicit formulas and numerical cross-checks corroborate.","rationale":"The paper's central claim is an exact equality between the optimal errors of the general coupled class and the uncoupled diagonal subclass. The lower bound E_coup >= e_2N follows cleanly from the rational-degree bound (Lemma S1) and does not depend on any external structural result. The upper bound E_diag <= e_2N does depend on the cited Zolotarev parity/type lemma, which is the reader's identified weakest assumption and the only genuinely load-bearing external input. I examined the proof chain around this lemma: the explicit Jacobi product in Lemma S3, the pole-zero interlacing, and the residue-sign calculation in Eq. (S32) are internally consistent, and the numerical verifications in Sec. S7 provide strong independent support that the analytic minimax error coincides with the constructed rational approximant. The parity/type theorem is classical and the paper's reliance on a published reference for it is standard practice. I therefore do not see a concrete route by which the central claim fails, and the correct verdict remains the reader's ACCEPT.","tokens_in":17284,"tokens_out":35299,"duration_ms":408043,"concrete_test":"Independently verify Lemma S3 without invoking Eq. (27) of Ref. [27]: (i) re-derive from Refs. [25,26] that, for even degree 2N, the minimax error over RR_{2N,2N}([-1,-k] U [k,1]) is attained by an odd rational function of type (2N-1,2N); and (ii) as a numerical spot-check, run a rational minimax solver on the full non-odd type-(4,4) class for k=0.1 (R=10) and compare the optimum with e_4(0.1) from Eq. (17). If the non-odd optimum is strictly below e_4(0.1), or if the analytical re-derivation fails, the upper-bound construction in Lemma S3 would not support the theorem.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The upper-bound step E_diag <= e_2N depends entirely on Lemma S3 (Supplemental Sec. S3, Eq. S19): the minimax sign approximant over the full type-(2N,2N) class can be chosen odd of type (2N-1,2N), so that after extracting the factor xi the remaining rational function lies in RR_{N-1,N}([1,R^2]) and has the positive-residue partial-fraction form needed in Eq. (13). This lemma is quoted from Ref. [27] and is not re-derived in the paper. If the odd extremal did not exist, a non-odd type-(2N,2N) approximant could in principle achieve an error below e_2N while the best physical type-(2N-1,2N) spectrum is strictly larger, breaking the chain of inequalities in Eq. (15). I found no internal inconsistency, and the concern is mitigated: the explicit Jacobi-product formula is stated, the paper verifies the implied equalities to 70 significant digits for 35 parameter pairs, and direct optimization in nine cases agrees with the analytic construction. The parity/type statement is also a classical Zolotarev result. On balance, this is a citation-level soft spot rather than a demonstrated flaw, so it does not change the reader's accept verdict.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies whether coherent intermode coupling can reduce the number of auxiliary modes needed to approximate a 1/|ω| spectrum over a finite two-sided frequency band in passive, number-conserving Gaussian Lindblad networks. The central result, Eq. (12), states that for any mode budget N and dynamic range R>1, the optimal maximum relative error for the general coupled class equals that of the uncoupled diagonal subclass, and both equal the degree-2N Zolotarev error e_{2N}(R^{-1}) for sign approximation on two intervals. The authors prove the lower bound via a rational-degree constraint imposed by the mode count, construct an explicit upper-bound realization using N independent damped auxiliary modes at zero detuning, and then derive closed-form inversions for the minimum mode count and maximum dynamic range in terms of elliptic integrals. The paper includes a supplemental material with detailed proofs of the rational-degree bound, the sign/inverse-square-root isometry, the positive-residue partial-fraction construction, and extensive numerical cross-checks.","tokens_in":17475,"tokens_out":14981,"duration_ms":164200,"significance":"If the central theorem holds, it settles a structural question in pseudomode and non-Markovian bath modeling: despite having O(N^2) real parameters, coupled passive Lindblad networks offer no advantage over independent damped modes for this canonical 1/|ω| benchmark. The paper's strengths include a clean reduction of the physical problem to a classical Zolotarev problem, an explicit construction with positive residues that is directly realizable as a CPTP Lindblad network, and exact elliptic-function formulas for mode count and dynamic range. The numerical verification is unusually thorough: the authors report 70-digit agreement between the Jacobi product formula and the modular equation in 35 cases, exact agreement of the closed-form mode-count formula with sequential search in all 84 tested cases, and agreement with direct nonconvex optimization to better than 10^{-4}. The derivation is parameter-free in the sense that no parameters are fitted to the target spectrum; all quantities are determined by the band edges and tolerance.","major_comments":[],"minor_comments":[{"comment":"The upper-bound chain in Eq. (15) relies entirely on Lemma S3, which asserts that the degree-2N Zolotarev problem admits an odd extremal of type (2N-1,2N) and that after factoring out the argument the remaining rational function has the required type. This lemma is quoted from the arXiv v1 version of Ref. [27] rather than proved in the paper. Because this is the only load-bearing step not derived self-containedly, I recommend that the authors either supply a proof in the Supplemental Material or state the precise theorem with equation numbers from the published version of Ref. [27]. The 70-digit numerical checks and agreement with direct optimization make a genuine error unlikely, but the citation should be made verifiable.","section":"Supplemental Sec. S3, Lemma S3"},{"comment":"In the direct-optimization cross-check, the text says that 'a linear program computes the minimax error level on the current exchange set' but does not give the variables or constraints of that linear program. Please add a one-sentence description of the LP (e.g., minimize the maximal deviation at the exchange points subject to the rational ansatz) so the procedure is reproducible.","section":"Supplemental Sec. S7"},{"comment":"The notation δα(R) should be typeset as δ_α(R) to avoid confusion with a Kronecker delta, and the subscript α should be defined explicitly as the power-law exponent.","section":"Supplemental Eq. (S59)"},{"comment":"The sentence 'Returning to the frequency variable ω then defines the corresponding spectrum' could be clearer: the spectrum is defined by J_opt(ω)=r_opt(ω^2), which is then identified with a sum of Lorentzians. Please rephrase to avoid the appearance of a circular definition.","section":"Main text, Eq. (13) and surrounding text"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is within the scope of the journal and the central claim appears sound. The only substantive concern is the reliance on an external Zolotarev lemma without a self-contained proof; this is local and fixable. I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi — quick read on arXiv:2608.04539. The result is real and worth taking seriously: for passive, number-conserving Gaussian Lindblad networks, coherent coupling between auxiliary modes doesn't change the best achievable error when approximating 1/|ω| over a finite band at fixed mode count. The paper pins that optimal error to the degree-2N Zolotarev error and shows the optimum is reached by N independent damped modes at zero detuning. That's a clean answer to a question the pseudomode literature mostly sidesteps, since prior work focused on how many modes suffice for a given time window, not on optimality at fixed budget.\n\nThe proof is structurally sound. The rational-degree bound from the mode count is elementary but correct, including the nonnormal/Jordan-block case. The lower bound via enlarging to the Zolotarev class is tight because the positivity and residue structure of the extremal odd rational function lets them realize it as a sum of Lorentzians. The inclusion argument closes the chain of inequalities cleanly. The explicit Jacobi product formula, the 70-digit checks in 35 cases, and the direct optimization agreement in nine cases all support the construction.\n\nThe one genuine soft spot is that the upper-bound step leans on Lemma S3, quoted from Ref. [27], that the degree-2N Zolotarev optimum can be chosen odd of type (2N-1,2N) with positive partial-fraction residues. That lemma isn't proved in the paper. It is a classical Zolotarev result, and the paper's own formulas and numerics corroborate it, so this is a citation-level gap rather than a demonstrated flaw. Still, for a theorem this clean, I'd like the supplement to either re-derive the needed parity/type statement or state it explicitly with a fuller reference, so the paper is self-contained on its load-bearing fact.\n\nThe scope is honestly stated: passive, number-conserving, vacuum state, single Hermitian bath operator, no white-noise feedthrough. Those restrictions matter; the result shouldn't be oversold beyond that class. Within that class, the paper delivers exact mode-count and dynamic-range formulas practitioners can use directly.\n\nThis paper deserves a serious referee. I'd send it out, asking the referee to check Lemma S3's provenance and the supplemental proof of positive residues. If that holds, accept.","headline":"Exact no-advantage theorem for 1/|ω| spectra: the proof holds up, with only a citation-level soft spot in the Zolotarev lemma.","tokens_in":18072,"tokens_out":2691,"would_cite":true,"duration_ms":30661,"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":"For passive Gaussian Lindblad networks, coherent intermode coupling cannot beat N independent damped modes at zero detuning when fitting a 1/|ω| spectrum: the optimal error is the degree-2N Zolotarev error.","keywords":["open quantum systems","non-Markovian baths","pseudomodes","Lindblad networks","1/f noise","Zolotarev approximation","rational approximation","noise spectroscopy"],"falsifier":"Take a small instance such as R=4, N=2 and globally search the full general coupled class of stable passive Lindblad networks for a maximum relative error over $[-4,-1]\\cup[1,4]$ strictly below $e_4(1/4)$; finding such a network would disprove the equality. A complementary check is to verify the imported positive-residue property for the known Zolotarev approximants: if some degree-2N optimal sign approximant necessarily has a negative residue, the physical realization step would fail.","tokens_in":17028,"feed_emoji":"⚛️","tokens_out":9172,"duration_ms":102212,"temperature":0.7,"pith_summary":"The paper asks whether, at a fixed number of auxiliary modes N, adding coherent intermode coupling, collective dissipation, or nonnormal structure can reduce the best possible maximum relative error when fitting a 1/|ω| spectrum over a finite two-sided frequency band. It proves the answer is no: the optimal error of the general coupled class equals the optimal error of the much simpler uncoupled diagonal subclass, and both equal the degree-2N Zolotarev error for rational approximation of the sign function. The optimum is achieved by N independent damped auxiliary modes at zero detuning, i.e., by a sum of Lorentzians. As a result, the O(N²) parameters available to a coupled network do not reduce the number of modes needed for a prescribed tolerance; they only move poles and residues within a fixed rational-degree constraint. This exact relation yields closed-form answers for the minimum mode count for a prescribed dynamic range and tolerance, and for the maximum dynamic range attainable at a fixed mode budget.","feed_headline":"Coupling cannot reduce auxiliary-mode count for 1/|ω| noise","feed_subtitle":"Even with all the extra coupling parameters, the best fit is N simple damped modes at zero detuning.","key_machinery":"The load-bearing object is the rational-degree bound: an $n$-mode passive Gaussian Lindblad network yields a spectrum $J(\\omega)$ whose numerator and denominator degrees are at most $2n-2$ and $2n$ after cancellation, so $\\omega J(\\omega)$ has rational type at most $(2n-1,2n)$, regardless of whether the drift matrix $M=\\Gamma+iH$ is normal or diagonalizable. This converts spectral fitting into the fourth Zolotarev problem of approximating $\\operatorname{sgn}(\\omega)$ on two disjoint intervals. The construction uses Zolotarev's odd optimal rational approximant of even degree $2N$; after factoring out $\\omega$, the remaining rational function of $\\omega^2$ lies exactly in the class $\\mathcal{RR}_{N-1,N}([1,R^2])$ for approximating $x^{-1/2}$. The crucial structural fact imported from Zolotarev theory is that this approximant can be chosen with all partial-fraction residues positive, which is what makes the realization by independent damped modes possible. The equality then follows by sandwiching: Zolotarev error is a lower bound for the coupled class, the diagonal subclass attains it, and the coupled class cannot do worse than the diagonal subclass.","core_discovery":"The central claim is the equality $E^{\\mathrm{coup}}_N(R)=E^{\\mathrm{diag}}_N(R)=e_{2N}(R^{-1})$, where $E^{\\mathrm{coup}}_N(R)$ and $E^{\\mathrm{diag}}_N(R)$ are the optimal maximum relative errors over $\\Omega_R=[-R,-1]\\cup[1,R]$ for the general coupled class and its uncoupled diagonal subclass, and $e_{2N}(R^{-1})$ is the optimal error of the degree-$2N$ fourth Zolotarev problem for approximating $\\operatorname{sgn}(\\omega)$ on $[-1,-1/R]\\cup[1/R,1]$. The proof has three legs: a rational-degree bound showing that any stable passive $n$-mode network produces $\\omega J(\\omega)$ of type at most $(2n-1,2n)$; an isometry showing that sign approximation on the two-sided band is minimax equivalent to relative-error approximation of $x^{-1/2}$ on $[1,R^2]$; and a realization argument showing that the optimal rational function has an all-positive partial-fraction expansion, so it is exactly a sum of Lorentzians generated by $N$ independent damped modes at zero detuning. Because the optimum lies inside the uncoupled diagonal subclass, the full coupled class cannot beat it.","pith_inferences":["Editorial extension: the no-advantage result relies on the combination of a single Hermitian bath operator and the exact $1/|\\omega|$ target; for multi-channel baths, asymmetric detailed-balance spectra, or a white-noise feedthrough term, the rational-degree bound changes and coherent coupling may plausibly help. This is not claimed in the paper.","Editorial extension: the exact mode-count formula provides a practical lower-bound test for any pseudomode or environment model: if a measured $1/|\\omega|$ noise band is claimed to be reproduced with fewer modes than the paper's $N_{\\min}$, the model either lies outside this passive class or the claim should be re-examined.","Editorial extension: one could apply the same rational-degree-plus-realization strategy to targets like $1/|\\omega|^{\\alpha}$ with $\\alpha\\neq 1$; the sign-function isometry would be replaced by a relative-error problem on one interval, and it is an open question whether positive-residue realizations persist for all such exponents."],"forward_implications":["No passive coupled Lindblad network with N modes can attain a maximum relative error below the Zolotarev error $e_{2N}(R^{-1})$; the rational-degree bound is the only constraint that matters for this benchmark.","The minimum number of auxiliary modes needed for a prescribed tolerance and dynamic range is given by an exact inversion of the error relation, with a large-range, small-tolerance asymptotic of the form $\\ln(4R)\\ln(4/\\varepsilon)/\\pi^2$.","The maximum positive-frequency dynamic range attainable with a fixed mode budget and tolerance follows by the same inversion; for 1%, 5%, and 10% tolerances each added mode asymptotically buys about 0.72, 0.98, and 1.16 decades of span.","The optimum is physically simple: N independent damped modes at zero detuning suffice, so one need not engineer coherent intermode coupling or collective dissipation to reach the fundamental limit.","Any candidate network's optimality can be certified by checking equioscillation of $\\omega J(\\omega)$ at the $2N+1$ alternation points, as done numerically in the paper."],"supporting_citations":[{"why":"Supplies the explicit odd optimal Zolotarev approximant of type (2N-1,2N), with the product form used to prove positive partial-fraction residues.","marker":"[27]"},{"why":"Gives the classical minimax rational approximation theory, the fourth Zolotarev problem, and the alternation criterion that certify the lower bound.","marker":"[25]"},{"why":"Provides the standard formulation and error parametrization of the fourth Zolotarev problem used in Eq. (10).","marker":"[26]"},{"why":"Quantum regression theorem used to identify the bath correlation function with $g^\\dagger e^{-Mt}g$ and hence the spectrum $J(\\omega)$.","marker":"[29]"},{"why":"Supports the statement that for linear system-bath coupling the bath statistics seen by the system are determined by the correlation function of the bath operator.","marker":"[6]"}],"fun_headline_variants":["Coupling can't cut mode count for 1/|ω| noise","No advantage to coupling in 1/|ω| Lindblad networks","Best 1/|ω| fit: uncoupled damped modes, not coupled","Zolotarev bound: coupled auxiliary modes don't help"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof's load-bearing assumption is an imported result from classical approximation theory: the optimal rational sign approximant can be chosen so that, when expanded in simple fractions, every coefficient is positive and every pole lies on the negative real axis. This is what turns a mathematical optimum into a physical array of damped modes; the paper cites this structural fact but does not prove it here.","fun_headline_variants_meta":{"raw":{"variants":["Coupling can't cut mode count for 1/|ω| noise","No advantage to coupling in 1/|ω| Lindblad networks","Best 1/|ω| fit: uncoupled damped modes, not coupled","Zolotarev bound: coupled auxiliary modes don't help"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000632,"raw_usage":{"total_tokens":2984,"prompt_tokens":1076,"completion_tokens":1908,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":692,"completion_tokens_details":{"reasoning_tokens":1826}},"tokens_in":692,"tokens_out":1908,"duration_ms":13513,"temperature":1.0,"reasoning_tokens":1826,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:25:23.634138+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a small instance such as R=4, N=2 and globally search the full general coupled class of stable passive Lindblad networks for a maximum relative error over $[-4,-1]\\cup[1,4]$ strictly below $e_4(1/4)$; finding such a network would disprove the equality. A complementary check is to verify the imported positive-residue property for the known Zolotarev approximants: if some degree-2N optimal sign approximant necessarily has a negative residue, the physical realization step would fail.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the explicit odd optimal Zolotarev approximant of type (2N-1,2N), with the product form used to prove positive partial-fraction residues."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the classical minimax rational approximation theory, the fourth Zolotarev problem, and the alternation criterion that certify the lower bound."},{"cited_title":"Istace and J.-P","cited_arxiv_id":null,"evidence_quote":"Provides the standard formulation and error parametrization of the fourth Zolotarev problem used in Eq. (10)."},{"cited_title":"Tamascelli, A","cited_arxiv_id":null,"evidence_quote":"Supports the statement that for linear system-bath coupling the bath statistics seen by the system are determined by the correlation function of the bath operator."}],"review_version":1}