{"id":"642189f4-28d5-4af4-bb0a-fbe841250ffc","arxiv_id":"2505.16669","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For a three-oscillator chain with two unequal-temperature baths, the local GKSL master equation is numerically closer to the exact dynamics for weak inter-oscillator coupling, while the global equation wins above a temperature-dependent critical coupling.","lead":"The authors compare two standard ways of deriving Markovian master equations, global and local, against an exactly solvable model of three coupled harmonic oscillators attached to two thermal baths at different temperatures. They report a temperature-dependent critical inter-oscillator coupling strength that decides which approximation is closer to the exact dynamics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed gc crossover rests on a single unvalidated finite-bath proxy: M=50 at t=50 with no convergence or spectral-density specification.","rationale":"The analytical framework is internally coherent: the Gaussian formalism, the global steady state, and the perturbative local steady state are derived without obvious algebraic errors, and the physical scenario is plausible. However, the paper's headline result is a numerical crossover extracted from a single finite-bath unitary simulation. The reader's weakest assumption identifies exactly the load-bearing issue: M=50 and t=50 are used as a proxy for the infinite-bath Markovian limit, yet no convergence tests or error bars are provided. My read agrees with this concern and adds precision: the missing specification of the bath spectral density and mode discretization makes it impossible to verify that the exact simulation implements the same J(ω) used in the master equations, and the mode spacing may exceed the global frequency splittings relevant near the claimed crossover. These issues do not invalidate the derivations, but they do prevent the central claim from being accepted as established. Since the reader already returned CONDITIONAL, no change to that verdict is needed.","tokens_in":20588,"tokens_out":7116,"duration_ms":65562,"concrete_test":"Specify the spectral density (e.g., flat J(ω)=γ in [0,ωc]) and a deterministic quadrature for the bath mode frequencies and couplings, then recompute gc as the g at which Δ(t)=F_loc−F_glb changes sign, for M=50, M=100 and M=200, using both the single-time value at t=50 and a time average over a quasi-stationary window (e.g., t∈[40,60] identified by a plateau criterion). If gc shifts by more than 20% with M or with the time-averaging choice, the claimed temperature-dependent crossover is not established by the current numerics.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, Eqs. (109)-(110) and Figure 4, is determined entirely by comparing F(ρ∞_loc, e^{tL_exc}[ρ∞_loc]) and F(ρ∞_glb, e^{tL_exc}[ρ∞_glb]) at one finite time t=50 in a unitary simulation with M=M'=50 modes per bath. The paper asserts without demonstration that this simulation effectively mimics infinite reservoirs and reaches a quasi-stationary plateau before recurrence; no convergence in M, no plateau-window analysis, and no error bars are given. This is load-bearing because any shift in the crossing of the two fidelity curves with M or t directly changes the reported gc and the claimed temperature dependence. The problem is compounded by an unspecified spectral density and discretization: the bath coupling functions h_{α,k} and Jα(ω) are never defined, so the exact simulation is not tied to the rates used in the local and global GKSL generators. In particular, with a hard cutoff ωc=3 and M=50 the mode spacing is about 0.06, larger than the global frequency splitting √2g≈0.014 near g=0.01, so the finite bath may not resolve the transitions on which the global master equation's superiority for larger g depends. A numerical artifact of the discretization could therefore produce the apparent crossover.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a three-site harmonic oscillator chain whose end sites are weakly coupled to two independent thermal baths at different temperatures. Within a Gaussian-state formalism, the authors derive the exact unitary dynamics of system plus finite baths, the global and local GKSL master equations, and their steady states (analytic for the global, second-order perturbative in the inter-oscillator coupling g for the local). They compare the fidelity between each approximate steady state and its evolution under the exact dynamics and claim the existence of a temperature-dependent critical coupling gc such that the local approach is more accurate for g<gc and the global one for g>gc. The numerical evidence for gc comes from finite-bath simulations with M=50 modes per bath at time t=50.","tokens_in":20823,"tokens_out":3552,"duration_ms":30084,"significance":"The question of local versus global master equations is of current interest, and the solvable Gaussian model offers a clean testbed. The authors provide explicit analytical formulas for the global steady state, a perturbative local steady state, and a Gaussian framework that is a valuable contribution. The claimed crossover, if backed by solid numerics, would be a useful practical criterion. However, the numerical validation is presently incomplete, and the paper would be considerably strengthened by convergence checks and a fully specified bath model.","major_comments":[{"comment":"The central claim about gc rests solely on finite-bath simulations with M=M'=50 modes per bath evaluated at t=50 for each g. The authors assert that when the baths contain a sufficiently large number of modes the dynamics reaches an almost stationary behavior before periodicity, but they do not demonstrate convergence in M, do not show a quasi-stationary plateau for every g used in Fig. 4, and report no error bars. If the crossing of the two fidelity curves shifts with M or with the evaluation time t, the reported gc values and their temperature dependence may be numerical artifacts. The authors should provide convergence checks (e.g., M=50 vs 100 vs 200) and verify that the ordering in Eqs. (109)-(110) is stable over a time window, not just at t=50.","section":"Section VII, Fig. 4, Eqs. (109)-(110)"},{"comment":"The bath spectral density J_alpha(omega) and the coupling functions h_{alpha,k} (or mode frequencies omega_{alpha,k}) used in the exact simulations are never specified. Equations (20) and (24) define the interaction via h_{ell,k}, h_{r,k}, but no explicit choice is given. The local and global master equations depend on J_alpha(omega_0) and J_alpha(epsilon_i), so without specifying the spectral density and its discretization, the exact simulation is not reproducibly tied to the GKSL rates. Moreover, with a hard cutoff omega_c=3 and M=50 modes per bath, the mode spacing is about 0.06, which is larger than the global frequency splitting sqrt(2) g near g=0.01; the finite bath may fail to resolve transitions on which the global master equation relies. Please specify the spectral density, the discretization prescription, and test sensitivity to the number of modes.","section":"Sections III, IV, VII"},{"comment":"The quasi-stationary behavior is shown for only three values of g (0.01, 0.1, 0.6), and no quantitative criterion is given for identifying the plateau. For g=0.01 the two fidelity curves are very close, as the authors note, so the sign of the inequality in Eq. (109) may be sensitive to the chosen evaluation time. The authors should establish the time window where the ordering of the two fidelities is stable and show that the same ordering persists over that window, not merely at t=50.","section":"Section VII, Fig. 3"}],"minor_comments":[{"comment":"The text after Fig. 4 says the green curve is for T_l=15, while the Fig. 4 caption says the green line is for T_l=20. This inconsistency should be corrected.","section":"Section VII, Fig. 4 caption and text"},{"comment":"The notation F^4(rho_loc^infty, rho_glb^infty) is confusing: it appears to denote the fourth power of the fidelity, but the exponent placement could be misinterpreted. Please clarify the notation throughout.","section":"Eq. (107)"},{"comment":"The arXiv title contains a typo: \"dyn amics\" should be \"dynamics\". There are also minor grammatical issues in the abstract, e.g., \"the behavior of fidelity between them versus inter-oscillator coupling depends on the two bath temperatures\" could be rephrased for clarity.","section":"Title and Abstract"},{"comment":"It would be helpful to state explicitly how the reduced covariance matrix of the three-oscillator chain is extracted from the total covariance matrix C_tot(t), since this is the quantity used in subsequent fidelity calculations.","section":"Section IV, Eq. (33)"},{"comment":"The paper would benefit from a more explicit discussion of how its results extend earlier exactly solvable comparisons of local and global master equations, e.g., Ref. [30] by Gonz\\'alez et al. and related works on bosonic chains.","section":"Introduction, references"}],"recommendation":"major_revision","confidential_remarks":"The paper has merit and the analytical parts appear sound. The main weakness is that the central numerical claim, the existence and temperature dependence of gc, is supported only by an unvalidated finite-bath proxy with M=50 modes per bath at a single time t=50, and the bath spectral density is left unspecified. This is fixable within the manuscript's scope by adding convergence analysis, defining the spectral density and discretization, and re-evaluating gc under those checks. I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper has a genuinely useful analytical core and a plausible but under-supported numerical claim. The authors derive local and global GKSL generators for a three-oscillator chain with asymmetric thermal baths, give an analytic global steady state, and produce a second-order perturbative local steady state. The equal-temperature fidelity expression and the point that the local steady state does not reduce to the global one as g→0 unless the baths are symmetric are clean, correct-looking results. Those parts alone make the paper worth having.\n\nThe new claim is the temperature-dependent critical coupling gc, below which the local master equation's steady state tracks the exact dynamics better than the global one, and above which the global wins. The comparison criterion—evolve each steady state under the exact unitary dynamics and compare fidelities—is sensible and not circular. The problem is that the numerical evidence for gc is under-specified where it matters.\n\nFirst, the bath spectral density J_alpha(omega) is never defined. The master equation rates depend on J at omega0 and at the global eigenfrequencies, but the finite-bath simulation uses arbitrary coupling constants h_{alpha,k}. Unless those h's come from the same J, the exact simulation and the GKSL generators are not describing the same physical bath. That is a reproducibility gap.\n\nSecond, there is no convergence analysis for the finite-bath simulation. M=50 and t=50 are asserted to mimic infinite reservoirs, and Figure 3 shows plateaus for three g values, but we don't see that the crossing in Figure 4 is stable as M grows. The stress-test worry is on target: with omega_c=3 and M=50 the bath mode spacing is ~0.06, larger than the global splitting sqrt(2)g≈0.014 at g=0.01. If the discrete bath can't resolve the transitions the global generator is built on, the small-g comparison may be biased toward the local approach, making the apparent crossover a numerical artifact. Not a fundamental flaw—the analytical framework is fine—but the right question to put to the authors.\n\nI largely agree with the reader's conditional verdict. The derivations are internally consistent, the literature is engaged honestly, and the claim is worth checking rigorously. Before publication the authors need to specify J(omega) and h_{alpha,k}, show convergence in M and t, and check the gc crossing under different discretizations.\n\nSend it to a serious referee. The analytical part deserves publication; the headline crossover needs the numerics pinned down first.","headline":"Solid analytic core, plausible crossover claim, but the numerical evidence is under-specified and the missing spectral density undermines reproducibility.","tokens_in":21334,"tokens_out":5569,"would_cite":true,"duration_ms":42382,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.30.-d","03.65.Yz","03.65.-w"],"model":"deepseek-v4-flash","headline":"Local master equations beat global ones for weak chain coupling, but global wins above a temperature-dependent critical strength.","keywords":["open quantum systems","GKSL master equation","local approach","global approach","harmonic chain","Gaussian states","fidelity","critical coupling"],"falsifier":"Run the same fidelity comparison with $M=100$ and $M=200$ modes per bath at the same parameters ($\\lambda=0.1$, $\\omega_0=1$, $\\omega_c=3$, $T_\\ell=10$, $T_r=8$). If the curves of $F(\\rho_{\\rm loc}^{\\infty}, e^{t\\mathcal{L}_{\\rm exc}}[\\rho_{\\rm loc}^{\\infty}])$ and the global analogue cross at a value of $g$ that shifts by more than the spacing between the $g$ values plotted, or if the apparent quasi-stationarity at $t=50$ breaks down at larger $M$, the reported $g_c$ values in Fig. 4 would be an artifact of the finite bath and finite-time cutoff rather than a physical crossover.","tokens_in":20332,"feed_emoji":"🔗","tokens_out":1743,"duration_ms":16305,"temperature":0.7,"pith_summary":"The paper studies a three-oscillator harmonic chain whose end sites each couple to an independent thermal bath at different temperatures. Because the total Hamiltonian is quadratic, the exact system-plus-bath dynamics can be solved numerically as a Gaussian map, and the authors compare that exact evolution with two standard approximations: the 'global' GKSL master equation, built from the full chain eigenmodes, and the 'local' one, built from the uncoupled oscillators. The central claim is that there exists a temperature-dependent critical inter-oscillator coupling strength $g_c$ such that for $g < g_c$ the local steady state deviates less from the exact dynamics, while for $g > g_c$ the global steady state is more accurate. If correct, this gives a practical rule for choosing between the two approaches in open harmonic networks: use the global equation once internal coupling is strong, and the simpler local equation for weak coupling, with the crossover shifting upward as the bath temperatures become more unequal.","feed_headline":"Temperatures pick the right master equation for an open chain","feed_subtitle":"A critical coupling strength $g_c$ decides when local beats global GKSL dynamics, and the threshold grows with bath-temperature imbalance.","key_machinery":"The analysis is carried out in the Gaussian-state formalism: quadratic Hamiltonians generate symplectic evolution, so all three approaches (exact, local, global) are reduced to linear equations for the covariance matrix $C(t)$ and displacement vector $d(t)$ of the form $\\dot{C} = M C + C M^\\dagger + N$. The central object is the covariance-matrix steady-state equation $M C_\\infty + C_\\infty M^\\dagger = -N$, which is solved exactly for the global approach in the normal-mode basis and perturbatively for the local approach by treating the inter-oscillator coupling $g$ as a small parameter. The comparison quantity is the Gaussian-state fidelity, evaluated between each approximate steady state and its own time-evolved state under the exact dynamics.","core_discovery":"The paper claims to reveal a temperature-dependent critical inter-oscillator coupling strength $g_c$ that sorts the regimes of validity of the local and global GKSL master equations for a dissipative harmonic chain. They compute the steady states of the local and global reduced dynamics (analytically for the global one, perturbatively to second order in $g$ for the local one), and compare how each steady state evolves under the exact unitary dynamics generated by the coupled chain plus two finite thermal baths. The finding is that when $g < g_c$, the fidelity $F(\\rho_{\\rm loc}^{\\infty}, e^{t\\mathcal{L}_{\\rm exc}}[\\rho_{\\rm loc}^{\\infty}])$ exceeds the analogous global fidelity, so the local steady state stays closer to the exact evolution; when $g > g_c$, the inequality reverses and the global approach is the better approximation. The critical value $g_c$ grows as the temperature difference between the two baths increases, and vanishes when the bath temperatures are equal, at which point the global approach is preferred for all $g$ considered.","pith_inferences":["The paper's crossover criterion is stated for a single three-oscillator chain; a natural extension is to test whether an analogous temperature-dependent $g_c$ appears for longer chains or for chains with different inter-site couplings, where the local approach is often the only tractable option.","The criterion compares only the steady states under exact evolution; one could extend this to full-time fidelity curves, which may reveal that the preferred approximation depends on the time window of interest, not just on $g$ and temperatures.","The claim that $g_c$ vanishes for equal bath temperatures suggests a symmetry-breaking mechanism: the temperature imbalance is what makes the local approximation competitive, and a testable prediction is that increasing the imbalance should monotonically raise the crossover coupling.","Because the exact dynamics is run with finite $M=50$ bath modes, a stronger test of the claimed crossover would be to check whether the curves of Fig. 3 remain stable as $M$ increases, before applying the result to infinite-reservoir settings."],"forward_implications":["If the claim is correct, practitioners can select between local and global GKSL master equations for an open harmonic chain based solely on whether the inter-oscillator coupling lies below or above a temperature-dependent threshold $g_c$.","The result implies that the computationally cheaper local approach is not merely an approximation to the global approach; for weak coupling it can be a genuinely more accurate description of the exact reduced dynamics.","The finding suggests that in transport calculations for bosonic chains, where bath temperatures differ, a crossover in accuracy will occur as internal coupling is increased, so conclusions drawn from a single master-equation choice may flip sign across that crossover.","The analytical forms of the steady-state covariance matrices (exact in the global case, second-order perturbative in the local case) can be reused as benchmarks for other approximations in open bosonic systems."],"supporting_citations":[{"why":"Compares local and global GKSL equations against an exactly solvable model, supplying the benchmark methodology that this paper extends to a three-oscillator chain with two baths.","marker":"[30]"},{"why":"Examines local versus global master equations for quantum thermal machines, providing the framework and motivation for judging which approximation better captures reduced dynamics.","marker":"[31]"},{"why":"Tests local versus global master equations with common and separate baths, a direct counterpart for the two-bath configuration studied here.","marker":"[33]"},{"why":"Structures the Gaussian-state formalism used to solve the exact dynamics of the quadratic system-environment model.","marker":"[39]"},{"why":"Supplies the theory of Gaussian states and maps that lets the paper reduce all three approaches to covariance-matrix equations.","marker":"[40]"},{"why":"Provides the background on Gaussian channels used to compute exact and approximate reduced dynamics.","marker":"[41]"},{"why":"Gives the solution formula for covariance-matrix dynamics of Gaussian semigroups, used to obtain both local and global time evolutions.","marker":"[43]"},{"why":"Provides the formula for Gaussian-state fidelity used in all numerical comparisons between steady states and exact-evolved states.","marker":"[45]"}],"fun_headline_variants":["Critical coupling picks local vs global master equation","Bath temperatures set the line where local beats global","Temperature gap shifts the local-global master equation crossover","Open chain: temperature decides local or global dynamics","When does local beat global? Ask the bath temperatures"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The finite-bath simulation with $M=50$ modes per bath and evaluation at time $t=50$ is treated as a faithful stand-in for the infinite-reservoir quasi-stationary dynamics, without checking convergence in $M$ or that $t=50$ lies inside the quasi-stationary window for every coupling $g$.","fun_headline_variants_meta":{"raw":{"variants":["Critical coupling picks local vs global master equation","Bath temperatures set the line where local beats global","Temperature gap shifts the local-global master equation crossover","Open chain: temperature decides local or global dynamics","When does local beat global? Ask the bath temperatures"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000224,"raw_usage":{"total_tokens":1514,"prompt_tokens":1050,"completion_tokens":464,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":666,"completion_tokens_details":{"reasoning_tokens":391}},"tokens_in":666,"tokens_out":464,"duration_ms":3908,"temperature":1.0,"reasoning_tokens":391,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:57:03.077611+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same fidelity comparison with $M=100$ and $M=200$ modes per bath at the same parameters ($\\lambda=0.1$, $\\omega_0=1$, $\\omega_c=3$, $T_\\ell=10$, $T_r=8$). If the curves of $F(\\rho_{\\rm loc}^{\\infty}, e^{t\\mathcal{L}_{\\rm exc}}[\\rho_{\\rm loc}^{\\infty}])$ and the global analogue cross at a value of $g$ that shifts by more than the spacing between the $g$ values plotted, or if the apparent quasi-stationarity at $t=50$ breaks down at larger $M$, the reported $g_c$ values in Fig. 4 would be an artifact of the finite bath and finite-time cutoff rather than a physical crossover.","supporting_citations":[{"cited_title":"Levy and R","cited_arxiv_id":null,"evidence_quote":"Compares local and global GKSL equations against an exactly solvable model, supplying the benchmark methodology that this paper extends to a three-oscillator chain with two baths."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Examines local versus global master equations for quantum thermal machines, providing the framework and motivation for judging which approximation better captures reduced dynamics."},{"cited_title":"Benatti, R","cited_arxiv_id":null,"evidence_quote":"Structures the Gaussian-state formalism used to solve the exact dynamics of the quadratic system-environment model."},{"cited_title":"Scali, J","cited_arxiv_id":null,"evidence_quote":"Supplies the theory of Gaussian states and maps that lets the paper reduce all three approaches to covariance-matrix equations."},{"cited_title":"Konopik and E","cited_arxiv_id":null,"evidence_quote":"Provides the background on Gaussian channels used to compute exact and approximate reduced dynamics."},{"cited_title":"Serafini, Quantum continuous variables: a primer of theoretical methods, first issued in paperback ed","cited_arxiv_id":null,"evidence_quote":"Provides the formula for Gaussian-state fidelity used in all numerical comparisons between steady states and exact-evolved states."}],"review_version":1}