{"id":"0c7d61c5-c253-426e-96bb-6bdc8d04fc5c","arxiv_id":"2510.06929","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Three first-law definitions give substantially different energy, heat, and work for two coupled thermal oscillator baths, with the minimal-dissipation approach showing coupling-enhanced work peaks.","lead":"This paper compares three ways to assign internal energy, heat, and work to two coupled sets of quantum oscillators that both act as thermal reservoirs. It finds the definitions disagree strongly even for weak coupling and large detuning, and that one approach produces work spikes that grow with coupling.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Minimal-dissipation work peaks coincide with near-singular inversion of the reduced propagator; they grow as Γ²/Δ in the ultrastrong limit but vanish at Δ=0, suggesting an artifact rather than a physical signature.","rationale":"The reader’s weakest_assumption identifies exactly this issue: the minimal-dissipation quantities depend on inverting the reduced propagator, and the reported peaks coincide with near-singularities of that inversion. I agree with the reader that this is the main caveat. The concern is real but not fatal. The paper’s central comparative claim—that the three definition sets differ substantially even in the dispersive regime—does not rely on the ultrastrong peak phenomenon; it is supported by the exact homogeneous formulas (Eqs. 54–58) and the numerical dispersive-regime results (Fig. 2). Moreover, the paper consistently labels the peaks as features of the minimal-dissipation definition, not as established physical work. Still, the abstract’s emphasis on work peaks that “increase with coupling strength” would be sharpened by an explicit acknowledgment that in the exact Δ=0 limit the peaks disappear, so the “unbounded” behavior is a non-commuting-limit artifact of the inversion. This is a clarifying caveat rather than a reason to reject. The reader’s ACCEPT verdict therefore stands unchanged.","tokens_in":19905,"tokens_out":10237,"duration_ms":88633,"concrete_test":"Evaluate the exact homogeneous formulas (56)–(58) at t = π/Γ for Δ = 10^{-k} Γ, k = 1,…,6, with fixed Γ and thermal occupations, and compare with the value at Δ = 0. If the peak height grows as Γ²/Δ but fails to approach the finite Δ=0 value, the advertised peaks are inversion artifacts. Separately, repeat the distributed-frequency numerical calculation with a regularized inverse (Φ_t + ε I)^{-1} for small ε; if the secondary peaks vanish or saturate on a scale set by ε, they are numerical/inversion artifacts rather than robust physical features.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline phenomenon of Sec. 5.1.2—secondary internal-energy peaks that are “virtually unbounded the deeper the ultrastrong regime”—is generated by B(t) (Eq. 58) and enters ΔU_md through Eq. 56 and δW_md through Eq. 69. For any finite Δ≠0, the reduced propagator Φ_t remains invertible, so these are not actual finite-time singularities. However, at the nominal peak time t≈π/Γ, the exact B(t) ≈ −Γ²/[2Δ²(1+π²/4)] for Δ≪Γ, so the contribution to ΔU_md scales as (ν1−ν2)B(t) ≈ −Γ²/Δ. Thus the peak height diverges as Δ→0. At exactly Δ=0, the prefactor (ν1−ν2)=0 in Eq. 56 kills the divergent term, and the energy variation is finite. The limit Δ→0 of the minimal-dissipation energy is therefore discontinuous. This non-uniform limit is a byproduct of inverting Φ_t, which becomes singular at isolated times in the resonant limit, not a feature of the physical model. The paper does not acknowledge this discontinuity or separate the inversion artifact from the genuine content of the minimal-dissipation definition. Since the abstract advertises these work peaks, this is the most load-bearing weakness.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript compares three sets of definitions of internal energy, heat, and work for a closed bipartite system consisting of two large but finite sets of linearly coupled harmonic oscillators, each initially in a thermal state. The authors derive exact analytical expressions for a homogeneous-frequency model and use exact diagonalization for distributed frequencies. They report that the two common asymmetric definitions (\"interaction\" and \"bare\") disagree with each other and with the minimal-dissipation approach even in the dispersive regime, that minimal-dissipation internal energy develops secondary peaks in the ultrastrong regime which are traced to work contributions, and that none of the definitions satisfies an energy balance between the two subsystems.","tokens_in":20158,"tokens_out":14605,"duration_ms":122690,"significance":"If the claims hold, the paper makes a useful and nontrivial point: in a fully autonomous, symmetric two-reservoir setup, dispersive conditions do not by themselves select a unique thermodynamic weak-coupling limit, and the choice of energy assignment affects even qualitative features. The analytical solution in Appendix A is explicit and checkable; the numerical results state their parameters; and the comparison is not circular, since the minimal-dissipation quantities are computed from the exact reduced propagator rather than assumed. The main advertised new feature, however, the growth of work peaks with coupling, is entangled with a near-singular inversion of the reduced one-body propagator, and the paper does not currently separate that artifact from a physical signature.","major_comments":[{"comment":"Equation (56) as printed has a sign error. For x=1 it gives +(ν1−ν2)B[n1+G]=+ΔB[...], while the Appendix (A.58)–(A.59) and the expansion in Eq. (65) are consistent with ΔU_md1 = +ν1G − ΔB[n1+G] and ΔU_md2 = −ν2G + ΔB[n2−G]. The correct general form is ΔU_md_x = ν_xG_x − (ν_x−ν_\\bar{x})B[n_x+G_x]. Since B(t)≤0, the printed sign reverses the stated dependence on the sign of Δ and would predict peaks where Fig. 3 shows dips. This is load-bearing: all subsequent qualitative statements about the sign and structure of minimal-dissipation energies inherit this error. It must be corrected.","section":"Sec. 4, Eq. (56)"},{"comment":"The secondary-peak phenomenon is not separated from an inversion artifact of the minimal-dissipation construction. The peak term arises from B(t), whose denominator is Δ²+Γ²cos²(Ωt/2). For finite Δ≠0, at t=π/Γ the exact contribution to ΔU_md1 is ≈Γ²/(2Δ) after the sign correction, which grows without bound as Δ→0; at exactly Δ=0 the prefactor (ν1−ν2)=0 removes the term, so the limit is discontinuous. This is precisely the regime where the factors α1(t), α2(t) entering Φ_t vanish and Φ_t is not invertible. Because the abstract advertises peaks that are \"virtually unbounded the deeper the ultrastrong regime,\" the paper should quantify this non-uniform limit, acknowledge that the peaks are properties of the chosen minimal-dissipation decomposition rather than of the underlying dynamics, and discuss whether they survive any physical regularization of the inversion.","section":"Sec. 5.1.2, Eqs. (58), (65), (69)"}],"minor_comments":[{"comment":"The expansion B(t) ≈ −(1/2)tan²(Ωt) appears to miss a half-angle: the exact B(t) in the ultrastrong limit behaves like −(1/2)tan²(Ωt/2), not tan²(Ωt). Please correct the argument.","section":"Sec. 4.2, Eq. (62)"},{"comment":"The caption says \"Energy (left), heat (middle) and work (left)\"; the last should presumably be \"work (right)\".","section":"Fig. 10 caption"},{"comment":"The notation −Δ/(2|α(t)|²) could be misread as −Δ|α(t)|²/2. Parentheses would improve clarity.","section":"Appendix A, Eq. (A.44)"}],"recommendation":"major_revision","confidential_remarks":"The sign error in Eq. (56) is likely a typo given the appendix and figures, but it sits in the central equation and must be fixed. The more substantive concern is the non-uniform Δ→0 limit of the minimal-dissipation peaks; I do not think it requires rejection, because the dispersive-regime incompatibility and the balance analysis are independent of that feature. A revision that fixes the sign, qualifies the ultrastrong-peak claim, and explicitly discusses the invertibility of Φ_t would make the paper acceptable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, this paper does something useful: it takes a fully solvable bipartite bosonic model with no natural system/bath split and compares three first-law schemes — interaction, bare, and minimal dissipation — on both sides. The analytical solution in Appendix A is explicit, and the numerics are exact diagonalization with clearly stated parameters. The central qualitative result, that the three schemes disagree substantially even in a far-detuned, weak-coupling-looking regime, and that the bare and interaction approaches do not converge to a common thermodynamic description, is well supported. That alone is worth a careful read.\n\nWhere I get less comfortable is the advertised secondary work peaks in the minimal dissipation approach. Those peaks come from the effective Hamiltonian built through the reduced propagator via L_t = dPhi_t/dt times Phi_t^{-1}, and in the homogeneous case Phi_t becomes singular at isolated times in the resonant limit. The peak height scales as Gamma^2/Delta and diverges as Delta goes to 0, but at exactly Delta = 0 the prefactor (nu1 - nu2) kills the divergent term, so the Delta-to-0 limit is discontinuous. The paper does not discuss this non-uniformity or separate an inversion artifact from physical content. The abstract advertises this as a feature of the approach, and Sec. 5.1.2 calls the peaks 'virtually unbounded.' That is the load-bearing weakness. It may be that the peaks are a genuine feature of the MD definition, but the current presentation does not establish that the divergence is not an artifact of the inversion. A referee should push on this.\n\nThe rest is solid: the interaction vs bare comparison is exact and clean, the dispersive-limit result that MD heat reduces to weak-coupling heat is a nice observation, and the collective-regime detuning reversal is a good example of how effective parameters matter. Citation practice is normal; the MD framework is the authors' own, but they compute it exactly here rather than assuming it.\n\nBottom line: send it to peer review. It deserves referee time. The referees should ask for a discussion of the invertibility of Phi_t and the Delta-to-0 discontinuity, and possibly a softened abstract. For my own work, the bare-vs-interaction disagreement in the dispersive regime is citable; the MD work peaks I would not cite without qualification.","headline":"A careful, analytically explicit comparison of three first-law definitions in a symmetric two-bath model; the dispersive-regime disagreement is solid, but the headline work peaks in the minimal-dissipation approach look like inversion artifacts and need scrutiny.","tokens_in":20668,"tokens_out":2992,"would_cite":true,"duration_ms":24833,"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":"The paper argues that in a two-reservoir model where each side is a thermal harmonic-oscillator bath, three established definitions of internal energy, heat, and work disagree even in the weakly coupled dispersive regime, and that the minim","keywords":["quantum thermodynamics","first law","strong coupling","internal energy ambiguity","minimal dissipation","harmonic oscillators","ultrastrong coupling","energy balance"],"falsifier":"Compute the minimal-dissipation heat and work for the same homogeneous model with a small regularizing term that keeps Phi_t invertible (e.g., adding a tiny dephasing or broadening the single-frequency limit), and check whether the tan² peaks at t=(2n+1)π/Γ become finite; or measure the renormalized local energy of one oscillator ensemble at those times in an ultrastrong-coupling experiment and look for the predicted divergent work spikes.","tokens_in":19745,"feed_emoji":"⚛️","tokens_out":4856,"duration_ms":39469,"temperature":0.7,"pith_summary":"This paper asks how to assign internal energy, heat, and work when neither side of a bipartition is clearly \"the system\" or \"the bath\". Using two large sets of coupled thermal harmonic oscillators, it compares two asymmetric definitions from the literature with a third, symmetric \"minimal dissipation\" definition. It finds that all three disagree substantially even when the coupling is weak and the detuning large, so dispersive conditions do not restore a unique first law. It also finds that the minimal-dissipation internal energy develops secondary peaks, traced to work-like contributions, whose magnitude grows without bound as the ultrastrong-coupling regime deepens. The conclusion matters because unambiguous thermodynamic bookkeeping at strong coupling is a prerequisite for quantum heat engines and thermal devices.","feed_headline":"Ultrastrong coupling yields unbounded work peaks in two-bath setup","feed_subtitle":"Even in the dispersive regime, three standard definitions of heat, work, and energy cannot be reconciled.","key_machinery":"The argument runs through the time-local reduced generator L_t = dot(Phi_t) Phi_t^(-1), where Phi_t is the exact evolution matrix for the first moments of one oscillator ensemble. Extracting its anti-Hermitian part gives the renormalized Hamiltonian K_t=(L_t^dagger - L_t)/(2i), which defines the minimal-dissipation internal energy, heat, and work. The principle of minimal dissipation selects this K_t as the unique Hamiltonian part. In the homogeneous model, Phi_t is built from phase factors alpha_1(t) and alpha_2(t); their zeros at t=(2n+1)π/Γ cause the tan² divergences in work.","core_discovery":"The central claim is that, in a fully autonomous model of two thermal oscillator ensembles interacting linearly, no agreement exists between the interaction, bare, and minimal-dissipation definitions of first-law quantities. Concretely, the variation of interaction energy is of the same order as the bare energy changes in the dispersive regime, so the bare and interaction definitions cannot be compatible; the minimal-dissipation definition also differs from both, and its heat reduces to the standard weak-coupling heat in that regime. In the ultrastrong regime the minimal-dissipation internal energy exhibits secondary peaks that diverge at times t=(2n+1)π/Γ, and the analytic expansion shows t","pith_inferences":["An extension the authors do not pursue: if the inversion singularity is an artifact, the unbounded peaks may be specific to the minimal-dissipation construction rather than physically realizable work, though the qualitative trend of growing peaks could survive in a regularized version with finite Phi_t.","A direct experimental probe could measure the renormalized level shifts spectroscopically in the ultrastrong regime and compare them with K_t; observation of the predicted tan² peaks would turn a formal divergence into a testable signature.","The result suggests that strong-coupling quantum thermal machines need an operational convention for \"work\" before efficiencies can be quoted, not just a Hamiltonian of the bare system.","The homogeneous special case is fully soluble, so extending the same comparison to non-Gaussian couplings or few-mode finite systems would show whether the convention dependence is generic."],"forward_implications":["Dispersive or weak-coupling conditions do not by themselves make interaction energy negligible; thermodynamic predictions depend on which first-law convention is chosen.","The minimal-dissipation heat is the only one that reduces to the standard weak-coupling heat in the dispersive limit, making it a candidate for a consistent extension.","The unbounded secondary peaks in ultrastrong coupling are work-like, so any attempt to extract work from these systems must contend with the chosen convention.","No set of definitions conserves energy between the two reservoirs; reported heat and work flows carry convention-dependent offsets of order 20–50% of the energy exchange.","The effective detuning set by collective eigenvalues, not the bare frequency difference, controls the direction and magnitude of energy flows."],"fun_headline_variants":["Quantum heat definitions clash even in weak-coupling regime","Two-bath model shows irreconcilable heat and work definitions","Work peaks diverge in ultrastrongly coupled thermal baths","Even dispersive coupling reveals irreconcilable quantum thermodynamic definitions","Local energy assignment impossible in coupled thermal reservoirs"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The minimal-dissipation results assume the reduced evolution map Phi_t is invertible at every time; at t=(2n+1)π/Γ its factors vanish, and the unbounded work peaks come from that inversion, so if the non-invertibility is an artifact the main quantitative claim weakens.","fun_headline_variants_meta":{"raw":{"variants":["Quantum heat definitions clash even in weak-coupling regime","Two-bath model shows irreconcilable heat and work definitions","Work peaks diverge in ultrastrongly coupled thermal baths","Even dispersive coupling reveals irreconcilable quantum thermodynamic definitions","Local energy assignment impossible in coupled thermal reservoirs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000521,"raw_usage":{"total_tokens":2331,"prompt_tokens":692,"completion_tokens":1639,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":436,"completion_tokens_details":{"reasoning_tokens":1558}},"tokens_in":436,"tokens_out":1639,"duration_ms":11056,"temperature":1.0,"reasoning_tokens":1558,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T11:03:50.090580+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the minimal-dissipation heat and work for the same homogeneous model with a small regularizing term that keeps Phi_t invertible (e.g., adding a tiny dephasing or broadening the single-frequency limit), and check whether the tan² peaks at t=(2n+1)π/Γ become finite; or measure the renormalized local energy of one oscillator ensemble at those times in an ultrastrong-coupling experiment and look for the predicted divergent work spikes.","supporting_citations":[],"review_version":1}