{"id":"eb8425ab-0398-404f-a4b4-06a7692f41b1","arxiv_id":"2608.09461","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Frequency-resolved thermometry of baths in a non-equilibrium steady state shows that bath-state dispersion tracks the turnover of quantum heat current, with effective temperatures reverting to initial bath values in the strong-coupling limit.","lead":"This paper introduces a frequency-selective thermometer, a weakly coupled two-level probe, to measure an effective temperature spectrum of each heat bath in a non-equilibrium steady state. Using numerically exact HEOM simulations of two quantum transport models, the authors find that the spectral dispersion tracks the heat-current turnover and that effective temperatures return to the initial bath values at strong coupling.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The turnover explanation rests on the strong-coupling flattening of T_eff toward T_in, but the HEOM truncation (N unspecified, J_k=1 or 2) is not shown to be converged in that regime; the flattening may be a truncation artifact.","rationale":"I read the paper as proposing a new operational thermometer and using it to give a physical picture of the turnover. The protocol is clearly described and the two models are appropriate. The weakest point is not the concept but the evidence base: all conclusions are drawn from HEOM data, and the error control for the numerics is asserted, not demonstrated. In the strong-coupling regime the hierarchy couples to many ADOs; with only one or two Padé terms per bath and no stated N, it is plausible that the hierarchy is truncated too shallowly. Because the claim is precisely that T_eff approaches T_in at strong coupling, a spurious flattening would produce the claimed explanation even if the true bath state remains out of equilibrium. The tiny probe coupling also raises conditioning questions. These issues are checkable by rerunning the numerics; they do not invalidate the concept. The reader's verdict of CONDITIONAL is appropriate, so I do not change it.","tokens_in":11889,"tokens_out":16000,"duration_ms":183833,"concrete_test":"Rerun Model I at λ=1 (and at the turnover maximum λ≈0.03) for representative probe frequencies ω=0.5 and 1.0 with hierarchy tiers N=4,5,6 and Padé orders J_k=2,3,4; require T_eff/T_in and Q_ss to change by <1% as N and J_k increase. Also compute Q_ss and probe populations by explicit long-time integration to t ≫ 1/η and compare with the nullspace solution; if T_eff/T_in at λ=1 moves by more than a few percent with N/J_k, the observed strong-coupling flattening is a truncation artifact and the central explanation fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim — that the turnover is explained by the bath state returning to its initial temperature in the strong-coupling limit — depends on the numerically computed T_eff(ω,λ) in Figs. 2 and 3. The HEOM implementation described in Sec. II A uses a finite hierarchy tier N and Padé order J_k (J_k=1 for Model I, J_k=2 for Model II). The paper states only that 'the convergence of this truncation is verified by systematically increasing N until the dynamics stabilizes,' but no convergence plots, N values, or error bars are presented. In the strong-coupling part of Fig. 2 (λ≈1), system-bath correlations are strongest, so the required N and J_k are largest. If N is too small, the hierarchy truncation artificially suppresses these correlations, which would mimic an 'effective decoupling' and make T_eff collapse to T_in. The same flattening is the sole evidence for the proposed mechanism. An additional numerical hazard is the probe coupling η=10^-8, which creates a ~10^-8 relaxation rate for the probe; the steady-state nullspace computation may be ill-conditioned unless special care is taken. Without convergence and error data, the central explanation is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a frequency-resolved thermometry protocol for characterizing non-equilibrium baths in quantum heat transport. A tunable two-level probe is weakly coupled to each bath, and the steady-state probe population ratio is converted into a frequency-dependent effective temperature T_eff(ω) via a detailed-balance relation. The protocol is applied with the hierarchical equations of motion (HEOM) to two models: a spin-boson system coupled to two baths and a two-qubit system with independent baths, both with Drude-Lorentz spectral densities. The authors compute the steady-state heat current using the Kato-Tanimura expression and observe the familiar turnover as a function of system-bath coupling λ. Their central physical claim is that the turnover is explained by the bath state itself: at strong coupling, T_eff(ω,λ) becomes nearly frequency-independent and approaches the initial bath temperature, indicating that the baths effectively decouple from the system and hence the current is suppressed. The protocol is presented as platform-agnostic, with HEOM used as the demonstration tool.","tokens_in":12146,"tokens_out":4342,"duration_ms":49160,"significance":"If the numerical results are reliable, the paper offers a genuinely operational diagnostic for non-equilibrium bath states and a new, physically intuitive explanation for the heat-current turnover: rather than being a purely system-side effect, the turnover is accompanied by a return of the bath's local thermal state toward its initial equilibrium. This could be a useful complement to existing explanations based on the quantum Zeno effect or system-bath hybridization. The strengths are the use of a standard, well-documented HEOM framework, a thermodynamically consistent heat-current definition, and a clear two-model demonstration. The main weaknesses are that the entire evidence is numerical, with no convergence data, error bars, or code/data provided, and that a key methodological assumption—single-probe simulation being equivalent to two-probe simulation—is asserted without proof or numerical verification. Because the central mechanism is inferred from small deviations in T_eff at strong coupling, the absence of convergence and uncertainty quantification is load-bearing.","major_comments":[{"comment":"The central evidence for the proposed mechanism is the flattening of T_eff(ω,λ) toward the initial bath temperature at strong coupling. This evidence is entirely numerical, yet the manuscript does not report the hierarchy truncation tier N, the convergence of the results with N, or the convergence with respect to the Padé order J_k (J_k=1 for Model I and J_k=2 for Model II). The statement in Sec. II A that 'the convergence of this truncation is verified by systematically increasing N until the dynamics stabilizes' is not backed by any data. At large λ, where system-bath correlations are strongest, an insufficient hierarchy depth could artificially suppress correlations and produce exactly the observed collapse of T_eff to T_in. Please provide convergence plots (or tables) for representative λ values across the full range shown in Fig. 2, including the largest λ, with N and J_k varied, and report estimated error bars on T_eff and Q_ss. Without this, the strong-coupling flattening cannot be distinguished from a truncation artifact.","section":"Sec. II A and Figs. 2–3"},{"comment":"The manuscript states that 'we simulate the system with a single probe attached at a time; this approach yields identical results to a full two-probe simulation.' This is a substantive methodological claim: the protocol requires that probes on different baths do not affect each other or the system-bath steady state, and that the presence of one probe does not alter the bath state measured by another. No proof or numerical comparison is given. Please demonstrate the equivalence explicitly, for both models, by comparing a full two-probe simulation with sequential single-probe simulations over the relevant λ and ω ranges, and by checking invariance of the extracted T_eff as the probe coupling η is varied (including values above and below η=10^-8). If the equivalence only holds asymptotically in η, state the asymptotic regime and quantify the residual error.","section":"Sec. II B, near Eq. (14)"},{"comment":"The protocol's validity rests on the probe being minimally invasive, and the paper states that this 'is verified by confirming that bare system observables and inter-reservoir heat currents remain invariant within numerical tolerance as probe-bath coupling (η) is varied.' However, no such verification is shown. Please provide numerical evidence: for a representative set of ω and λ, plot system observables and Q_ss versus η, and show that the variations are below the numerical tolerance and below the observed T_eff changes. This is particularly important because η=10^-8 creates a very slow probe relaxation channel, and it is not obvious that the steady-state null-space computation is well conditioned at that scale.","section":"Sec. II C, 'Minimal Invasiveness'"},{"comment":"The paper interprets the convergence of T_eff(ω,λ) to T_in at strong coupling as evidence of 'effective decoupling' and uses this to explain the turnover. As presented, this is a correlation between two computed quantities (T_eff and Q_ss), not a demonstrated causal mechanism. Moreover, the observed deviations in Fig. 3 are very small (of order 0.1–1%), and without error bars it is unclear whether the trend is significant. Please provide a quantitative connection: for example, show that the frequency dispersion of T_eff, defined by a suitable measure such as max_ω |T_eff(ω)-T_in| or the variance across ω, peaks at the same λ as the heat current and decays on the same scale as the current suppression. Alternatively, state explicitly what test would distinguish the proposed bath-return mechanism from the quantum Zeno or hybridization explanations mentioned in Sec. IV. As written, the explanatory claim is stronger than the evidence supports.","section":"Sec. III and Sec. IV"}],"minor_comments":[{"comment":"The title contains a typo: 'T urnover' should be 'Turnover'. Also, in Sec. I the phrase 'are are presented' should be 'are presented'.","section":"Title and Abstract"},{"comment":"The axis labels in Fig. 2 appear as 'W eak', 'Intermediate', 'Strong'; these spacing artifacts should be corrected for readability.","section":"Fig. 2 labels"},{"comment":"The definition of B_k,k' as (i/ħ)^2 [[V_k, V_k'], V_k] is asymmetric; please clarify whether the intended object is [[V_k,V_k'],V_k'] or a symmetrized version, since the last terms involve commutators with V_k'.","section":"Eq. (18)"},{"comment":"The Padé decomposition formulas are standard but the notation for the matrix Λ and the eigenvalues λ_i conflicts with the system-bath coupling strength λ_k used in the main text; this is a minor but potentially confusing notational clash.","section":"Appendix A"},{"comment":"The heat current expression in Eq. (18) is stated to follow from Kato and Tanimura, but the intermediate steps are not shown; a brief derivation or reference to the specific equation in Ref. 41 would help readers verify the terminator contributions.","section":"Sec. II D"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely topic and the proposed thermometry protocol is interesting, but the numerical evidence is not yet at the standard required for the central claim. I would encourage the editor to request the missing convergence and invasiveness data, as well as a quantitative link between T_eff dispersion and heat current, before considering publication. The 'single probe equals two probes' assertion is particularly important to verify, as it is central to the protocol and currently unsupported. If the authors can provide the requested numerical checks, the paper could become a valuable contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this paper gives you a way to measure a frequency-resolved effective temperature of a non-equilibrium bath using a weakly coupled two-level probe, and it reports a new observation—in the strong-coupling limit the temperature spectrum flattens back to the initial bath temperature, right where the heat current turns over. If that observation holds, it is a genuinely useful diagnostic and a plausible explanation for the turnover.\n\nWhat is actually new: the probe idea itself comes from earlier work (refs 21,22), but building the full T_eff(ω) spectrum from HEOM steady states and watching it collapse with coupling is new. The paper uses standard, thermodynamically consistent heat-current definitions, and it demonstrates the effect in two different models, which helps generality. The protocol is platform-agnostic, which is a practical plus.\n\nThe soft spot is the one the reader flagged: the entire evidence base is numerical, and the numerics are under-documented. No convergence data, no hierarchy depth N, no error bars, no code. That matters here because the interesting regime is strong coupling, where HEOM truncation is most likely to distort the result. The flattening of T_eff to T_in could in principle be an artifact of an under-converged hierarchy, and the paper gives you no way to rule that out. The probe coupling at η=10^-8 also deserves a stability check, since the steady-state nullspace can be ill-conditioned at that scale. The single-probe equivalence is asserted without demonstration; it is probably true, but it is easy to check and should be shown.\n\nI don't think these problems sink the paper. The central protocol is reasonable and the observation is interesting enough to be worth a serious referee. But the turnover explanation is an interpretation, not a derivation: you see T_eff flatten and the current drop, and the paper calls that 'effective decoupling.' That is a fair hypothesis, but it needs the convergence data to back it up, and ideally a check that the flattening is not an artifact of the probe or the truncation.\n\nFor a reader in quantum thermodynamics, this is a useful contribution and a good starting point for further work. I would cite it once the numerics are verified. Recommendation: send to peer review, and ask the authors for convergence plots, hierarchy parameters, and the code or data. If those check out, the paper is fine.\n\nThat's my read.","headline":"A useful new probe for non-equilibrium baths, with a striking observation that deserves careful numerical verification before the turnover explanation is taken as established.","tokens_in":12693,"tokens_out":2537,"would_cite":true,"duration_ms":24872,"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":"A frequency-resolved thermometer protocol shows that the turnover of the steady-state heat current is caused by the bath spectrum flattening to its initial temperature at strong coupling.","keywords":["quantum heat transport","turnover effect","effective temperature","frequency-resolved thermometry","hierarchical equations of motion","non-equilibrium steady state","spin-boson model","strong coupling"],"falsifier":"Recompute the strong-coupling steady states with progressively larger HEOM hierarchy tiers and Padé orders and check whether $T_{\\rm eff}(\\omega)$ still flattens to the initial bath temperature; if the flattening recedes or shifts as the truncation is tightened, the decoupling signature is numerical rather than physical. In parallel, a tunable-coupling experiment could measure the probe's population ratio together with the heat current and look for the simultaneous return of the equilibrium ratio and suppression of current.","tokens_in":11656,"feed_emoji":"🌡️","tokens_out":7678,"duration_ms":68258,"temperature":0.7,"pith_summary":"This paper introduces a way to 'see inside' a heat bath that drives a quantum system, by attaching a tunable two-level thermometer and reading its steady-state population ratio at each probe frequency. That ratio defines a frequency-resolved effective temperature $T_{\\rm eff}(\\omega)$: flat in equilibrium, dispersive when the bath is out of equilibrium. Applied to a spin coupled to two reservoirs and to two coupled spins each in its own reservoir, the protocol attributes the heat-current turnover to the bath's own thermal state. At weak coupling the effective temperature peaks at the system's resonant frequencies; at intermediate coupling the peaks broaden and the spectrum spreads widest; at strong coupling every probe frequency reports the original bath temperature, meaning system and bath have effectively decoupled. The paper concludes this decoupling, read off the bath spectrum itself, is what suppresses the steady-state heat current.","feed_headline":"Quantum heat-current turnover traced to bath decoupling","feed_subtitle":"A frequency-resolved probe sees the bath spectrum flatten to its original temperature exactly where the current fades.","key_machinery":"The central object is the frequency-resolved effective temperature spectrum $T_{\\rm eff}(\\omega)$, extracted from a weakly coupled two-level probe through the detailed-balance relation $P_e/P_g=\\exp(-\\hbar\\omega/k_B T_{\\rm eff}(\\omega))$. The protocol is platform-agnostic: the steady-state populations come from a numerically exact simulation, here the hierarchical equations of motion (HEOM) with Padé decomposition of the Drude-Lorentz bath correlation functions. The probe coupling $\\eta$ is kept much smaller than the system-bath coupling $\\lambda$ so the measurement is minimally invasive, and the heat current is evaluated from first-tier auxiliary density operators. The spectral dispersion — the variation of $T_{\\rm eff}$ across probe frequencies — is the diagnostic that carries the argument: flatness signals equilibrium, dispersion signals non-equilibrium, and the collapse to a flat spectrum at the initial bath temperature signals the decoupling that produces the turnover.","core_discovery":"The central claim is that the turnover effect in non-equilibrium quantum heat transport can be explained by how the bath's local thermal state evolves with system-bath coupling strength, witnessed by a frequency-selective thermometer. A weakly coupled two-level probe with transition frequency $\\omega$ defines $T_{\\rm eff}(\\omega)=\\hbar\\omega/[k_B\\ln(P_g/P_e)]$ from its steady-state populations, and scanning $\\omega$ gives a spectrum. For both the non-equilibrium spin-boson model and the two-qubit model, the spectrum shows three regimes: at weak coupling, deviations from the initial bath temperature concentrate near the bare system transitions; at intermediate coupling, the resonances broaden and shift while the spectral dispersion is maximal; at strong coupling, $T_{\\rm eff}(\\omega)$ collapses to the initial bath temperature at every frequency. The authors interpret that collapse as effective decoupling between the system and its baths, which is the physical origin of the declining branch of the turnover curve.","pith_inferences":["A natural extension is to probe the bath with a multi-level or harmonic thermometer, which would yield a fuller effective spectral function and could reveal whether the strong-coupling flattening is exact or only asymptotic in $\\lambda$.","The framework suggests that the critical coupling of the turnover is set by the coupling where the bath-state deviation is maximal, so the peak location should be predictable from the growth and then shrinkage of the spectral dispersion alone.","The same population-ratio thermometer could be applied to fermionic or phononic junctions, where an equivalent energy-resolved probe would test whether the decoupling picture generalizes beyond bosonic baths."],"forward_implications":["The frequency dispersion of $T_{\\rm eff}$ can serve as a quantitative, probe-based measure of how far a bath is from equilibrium in any non-equilibrium steady state.","In the intermediate-coupling regime, where dispersion is largest, the bath state differs markedly from its initial thermal state, so methods that freeze the bath state will misestimate transport there.","The strong-coupling collapse of $T_{\\rm eff}$ to the initial bath temperature indicates that the turnover's declining branch is a decoupling effect that a weak-coupling or frozen-bath theory cannot capture.","Because the protocol only requires the steady-state population of a weakly coupled probe, it can be attached to any exact solver, not just HEOM.","The three identified regimes — resonant, non-resonant, and decoupled — give a concrete map of the coupling ranges where common approximations such as Redfield or the non-interacting-blip approximation can be trusted."],"supporting_citations":[{"why":"introduces the hierarchical equations of motion that supply the numerically exact steady states used throughout the paper.","marker":"[26]"},{"why":"provides the HEOM implementation with the Drude-Lorentz spectral density decomposition used for the simulations.","marker":"[29]"},{"why":"supplies the frequency-resolved effective-temperature idea via a weakly coupled probe.","marker":"[21]"},{"why":"gives the thermodynamically consistent heat-current expression evaluated from first-tier auxiliary density operators.","marker":"[40]"},{"why":"reports the two-qubit turnover effect this paper's protocol is designed to explain.","marker":"[16]"},{"why":"gives the analytic representation of Drude-Lorentz bath correlation functions used in the decomposition.","marker":"[39]"},{"why":"introduces the Padé spectrum decomposition that accelerates HEOM convergence.","marker":"[44]"}],"fun_headline_variants":["Quantum heat turnover pinned to bath thermal collapse","Frequency-resolved probe spots bath decoupling in heat current","Heat current decline traced to bath spectrum flattening","Bath thermometer reveals why quantum heat current turns over"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument rests on the hierarchical-equations-of-motion calculation being numerically converged at every coupling strength studied; if the finite hierarchy depth and Padé order are too shallow at large $\\lambda$, the reported flattening of $T_{\\rm eff}$ toward the initial bath temperature could be a truncation artifact.","fun_headline_variants_meta":{"raw":{"variants":["Quantum heat turnover pinned to bath thermal collapse","Frequency-resolved probe spots bath decoupling in heat current","Heat current decline traced to bath spectrum flattening","Bath thermometer reveals why quantum heat current turns over"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000107,"raw_usage":{"total_tokens":1018,"prompt_tokens":892,"completion_tokens":126,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":508,"completion_tokens_details":{"reasoning_tokens":65}},"tokens_in":508,"tokens_out":126,"duration_ms":2352,"temperature":1.0,"reasoning_tokens":65,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:05:26.214846+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the strong-coupling steady states with progressively larger HEOM hierarchy tiers and Padé orders and check whether $T_{\\rm eff}(\\omega)$ still flattens to the initial bath temperature; if the flattening recedes or shifts as the truncation is tightened, the decoupling signature is numerical rather than physical. In parallel, a tunable-coupling experiment could measure the probe's population ratio together with the heat current and look for the simultaneous return of the equilibrium ratio and suppression of current.","supporting_citations":[{"cited_title":"Alicki \\ and\\ author D","cited_arxiv_id":null,"evidence_quote":"supplies the frequency-resolved effective-temperature idea via a weakly coupled probe."},{"cited_title":"Kato \\ and\\ author Y","cited_arxiv_id":null,"evidence_quote":"gives the thermodynamically consistent heat-current expression evaluated from first-tier auxiliary density operators."},{"cited_title":"Kato \\ and\\ author Y","cited_arxiv_id":null,"evidence_quote":"reports the two-qubit turnover effect this paper's protocol is designed to explain."},{"cited_title":"Ritschel \\ and\\ author A","cited_arxiv_id":null,"evidence_quote":"gives the analytic representation of Drude-Lorentz bath correlation functions used in the decomposition."}],"review_version":1}