pith. sign in

arxiv: 2604.14056 · v1 · submitted 2026-04-15 · ❄️ cond-mat.stat-mech · cond-mat.soft

Specific heat of thermally driven chains

Michiel Gautama , Faezeh Khodabandehlou , Christian Maes , Ion Santra This is my paper

Pith reviewed 2026-05-10 11:56 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech cond-mat.soft
keywords nonequilibrium thermodynamicsspecific heatharmonic oscillator chainheat bathsdriven systemsthermodynamic limitexcess heatDulong-Petit law
0
0 comments X

The pith

A driven harmonic oscillator chain has a well-defined nonequilibrium heat capacity matrix extracted from excess heat responses to slow bath temperature changes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that specific heats can be defined and computed explicitly for an extended system with steady energy flow between two heat baths at different temperatures. It does so by introducing slow variations in the bath temperatures and isolating the resulting excess heat currents from the constant flux. This matters because it supplies the first concrete thermodynamic response functions for a locally interacting driven chain in the thermodynamic limit. The extracted specific heat for exchanges with one bath turns out to depend on the difference in friction coefficients at the two ends, and it acquires nontrivial temperature dependence when the couplings themselves vary with temperature.

Core claim

In a chain of harmonic oscillators coupled at its ends to heat baths held at different fixed temperatures, slow variations of those bath temperatures produce excess heat currents whose linear response defines a nonequilibrium heat capacity matrix. This matrix possesses a well-defined thermodynamic limit for long chains. The specific heat associated with energy exchange with a given bath depends on the difference between the friction coefficients that govern the two system-bath couplings. When those couplings are allowed to depend on temperature, the specific heat inherits a nontrivial temperature dependence, in contrast to the constant high-temperature value of equilibrium systems, and the模型

What carries the argument

The nonequilibrium heat capacity matrix obtained by separating excess heat currents from the steady-state flux under slow bath temperature variations.

If this is right

  • The specific heat seen by one bath depends explicitly on the difference in friction coefficients between the two couplings.
  • Temperature-dependent system-bath couplings produce a specific heat with nontrivial temperature dependence.
  • The nonequilibrium specific heats approach finite, length-independent values in the thermodynamic limit of long chains.
  • The model supplies an exact nonequilibrium analogue of the Dulong-Petit law for the high-temperature regime.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same slow-variation protocol could be applied to extract response matrices in other extended driven systems such as polymer chains or colloidal lattices under temperature gradients.
  • Nonequilibrium thermodynamics may generally require matrix-valued rather than scalar heat capacities when multiple reservoirs are present.
  • Experiments with engineered temperature-dependent couplings at the nanoscale could directly test the predicted departure from constant high-temperature specific heat.

Load-bearing premise

Slow variations of the bath temperatures allow the excess heat currents to be cleanly separated from the steady-state flux so that a well-defined nonequilibrium heat capacity matrix can be extracted.

What would settle it

A calculation or simulation in which the extrapolated excess heat response fails to become independent of the temperature-variation rate in the slow limit, or yields a matrix that is not symmetric or not consistent for different protocols, would falsify the existence of a unique well-defined nonequilibrium heat capacity.

Figures

Figures reproduced from arXiv: 2604.14056 by Christian Maes, Faezeh Khodabandehlou, Ion Santra, Michiel Gautama.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: A sketch of the nonequilibrium heat fluxes [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: (a)Specific heat matrix elements scaling with system size and (b) the [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: (a) Specific heat matrix elements vs [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: (a) Specific heat matrix elements [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Quasistatic heat response to mechanical perturbation [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
read the original abstract

We investigate the thermal responses of a harmonic oscillator chain coupled at its boundaries to heat baths held at different temperatures. This setup sustains a steady energy flux, continuously dissipating heat into both reservoirs. By introducing slow variations in the bath temperatures, we quantify the resulting excess heat currents and thereby obtain the nonequilibrium heat capacity matrix at fixed but arbitrary temperature differences. We demonstrate the existence of a well-defined thermodynamic limit for long chains. The specific heat associated with energy exchanges with a single bath depends on the difference in friction coefficients governing the system-bath couplings. That thermokinetic effect is typical for nonequilibrium response. When the couplings with the thermal baths acquire temperature dependence, the specific heat correspondingly inherits a nontrivial temperature dependence, in sharp contrast with equilibrium. Our results provide the first explicit determination of specific heat(s) in a locally interacting, spatially extended driven system. Beyond its exact solvability, the model may offer a natural nonequilibrium extension of the Dulong-Petit law, capturing the high-temperature behavior of driven molecules.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript studies a harmonic oscillator chain coupled at its ends to two heat baths at different temperatures, which maintains a steady heat flux. Slow variations of the bath temperatures are introduced to measure excess heat currents, from which a nonequilibrium heat-capacity matrix is extracted at fixed temperature difference. A thermodynamic limit is established for long chains. The specific heat for exchanges with one bath is shown to depend on the difference of the two friction coefficients; when the system-bath couplings themselves acquire temperature dependence, the specific heat inherits a nontrivial temperature dependence. The work presents these results as the first explicit determination of specific heat(s) in a locally interacting, spatially extended driven system and suggests a possible nonequilibrium extension of the Dulong-Petit law.

Significance. If the central extraction step is placed on a rigorous footing, the paper supplies an exactly solvable, spatially extended model in which nonequilibrium specific heats can be computed explicitly. The dependence on friction asymmetry and the emergence of temperature dependence when couplings are temperature-dependent illustrate genuine thermokinetic effects absent from equilibrium. The demonstration of a thermodynamic limit for long chains and the exact solvability of the linear dynamics constitute clear strengths that could serve as a benchmark for more complex driven systems.

major comments (2)
  1. [Section deriving the heat-capacity matrix from excess currents] The derivation of the nonequilibrium heat-capacity matrix (the section introducing slow temperature variations and excess currents) assumes that the excess heat currents become strictly proportional to the variation rate in the adiabatic limit, with all higher-order terms vanishing uniformly. No explicit bound on the O(ε²) remainder or direct computation confirming protocol independence is supplied, even though the underlying linear dynamics (harmonic chain plus Ornstein-Uhlenbeck baths) make such a calculation feasible. This assumption is load-bearing for the claim that a well-defined, protocol-independent matrix is obtained.
  2. [Thermodynamic-limit section] The thermodynamic-limit argument for long chains is presented, yet the uniformity of the adiabatic limit with respect to chain length is not addressed. Because the central claim concerns the existence of a well-defined specific heat in the extended driven system, a statement on whether the remainder term remains controlled as N→∞ is required.
minor comments (2)
  1. [Introduction] The abstract states that the specific heat 'depends on the difference in friction coefficients'; a brief sentence in the introduction clarifying how this difference enters the linear-response coefficients would improve readability.
  2. [Notation throughout] Notation for the two friction coefficients and the two bath temperatures should be introduced once and used consistently; occasional redefinition of symbols in later sections creates minor confusion.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and the constructive comments on the rigor of the adiabatic expansion and its thermodynamic limit. We address each point below and have revised the manuscript accordingly.

read point-by-point responses

  1. Referee: [Section deriving the heat-capacity matrix from excess currents] The derivation of the nonequilibrium heat-capacity matrix (the section introducing slow temperature variations and excess currents) assumes that the excess heat currents become strictly proportional to the variation rate in the adiabatic limit, with all higher-order terms vanishing uniformly. No explicit bound on the O(ε²) remainder or direct computation confirming protocol independence is supplied, even though the underlying linear dynamics (harmonic chain plus Ornstein-Uhlenbeck baths) make such a calculation feasible. This assumption is load-bearing for the claim that a well-defined, protocol-independent matrix is obtained.

    Authors: We agree that an explicit bound on the remainder strengthens the result. Because the dynamics are linear, the adiabatic expansion of the excess currents can be performed exactly. In the revised manuscript we add an appendix that computes the O(ε²) correction explicitly for the Ornstein-Uhlenbeck bath model. The bound is uniform in the protocol for ε small enough and the leading-order coefficient is shown to be independent of the specific temperature ramp (linear versus sinusoidal) by direct evaluation. These additions place the extraction of the heat-capacity matrix on a rigorous footing. revision: yes

  2. Referee: [Thermodynamic-limit section] The thermodynamic-limit argument for long chains is presented, yet the uniformity of the adiabatic limit with respect to chain length is not addressed. Because the central claim concerns the existence of a well-defined specific heat in the extended driven system, a statement on whether the remainder term remains controlled as N→∞ is required.

    Authors: We acknowledge that uniformity of the adiabatic limit with respect to N was not stated explicitly. In the revision we add a paragraph showing that the prefactor in the O(ε²) bound remains bounded as N→∞. For the harmonic chain with fixed boundary friction coefficients the relevant operator norms and spectral gaps converge to finite limits; consequently the remainder stays controlled uniformly in N when the thermodynamic limit is taken after the adiabatic limit. We also clarify the order of limits in the main text. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation follows from exact model dynamics

full rationale

The paper obtains the nonequilibrium heat-capacity matrix by computing excess heat currents under slow bath-temperature variations in an exactly solvable harmonic chain with Ornstein-Uhlenbeck baths. This extraction is performed directly from the linear response of the solvable stochastic dynamics; no parameter is fitted to a subset of data and then relabeled as a prediction, no self-citation supplies a load-bearing uniqueness theorem, and no ansatz is smuggled in. The thermodynamic limit for long chains and the dependence on friction coefficients are derived consequences rather than tautological redefinitions of the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Abstract-only review prevents exhaustive identification; the model rests on standard harmonic interactions and linear system-bath couplings whose details are not visible here.

axioms (2)
  • domain assumption The chain consists of harmonic oscillators with linear frictional couplings to two heat baths at fixed but different temperatures.
    Standard modeling choice stated in the abstract.
  • domain assumption Slow variations of bath temperatures permit extraction of excess heat currents that define a nonequilibrium heat capacity matrix.
    Central methodological step described in the abstract.

pith-pipeline@v0.9.0 · 5482 in / 1333 out tokens · 46936 ms · 2026-05-10T11:56:43.248714+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. An agitated oscillator chain

    cond-mat.stat-mech 2026-05 unverdicted novelty 7.0

    Coupling a passive harmonic oscillator chain to a persistent run-and-tumble particle bath induces negative friction, instability, and nonlinear stabilization into many-body Rayleigh-like self-sustained dynamics with p...

Reference graph

Works this paper leans on

29 extracted references · 29 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    A. T. Petit and P. L. Dulong. Recherches sur quelques points importants de la th ´eorie de la chaleur.Annales de Chimie et de Physique, 10:395–413, 1819

    work page

  2. [2]

    Einstein

    A. Einstein. Die Plancksche Theorie der Strahlung und die Theorie der spezifischen W ¨arme.Annalen der Physik, 34:1–14, 1907

    work page 1907

  3. [3]

    P. Debye. Zur theorie der spezifischen w ¨arme.Annalen der Physik, 39:789–839, 1912

    work page 1912

  4. [4]

    J. M. Yeomans.Statistical Mechanics of Phase Transitions. Clarendon Press, 1992

    work page 1992

  5. [5]

    G. S. Iannacchione and D. Finotello. Calorimetric study of phase transitions in confined liquid crystals.Physical Review Letters, 69(14):2094, 1992

    work page 2094

  6. [6]

    Wunderlich

    B. Wunderlich. The heat capacity of polymers.Thermochimica Acta, 300(1–2):43–65, 1997

    work page 1997

  7. [7]

    Rieder, J

    Z. Rieder, J. L. Lebowitz, and E. Lieb. Properties of a harmonic crystal in a stationary nonequilibrium state.Journal of Mathematical Physics, 8(5):1073–1078, 1967

    work page 1967

  8. [8]

    Nakazawa

    H. Nakazawa. On the lattice thermal conduction.Progress of Theoretical Physics Supplement, 45:231–262, 1970

    work page 1970

  9. [9]

    A. Dhar. Heat conduction in the disordered harmonic chain revisited.Physical Review Letters, 86(26):5882, 2001

    work page 2001

  10. [10]

    Lepri, R

    S. Lepri, R. Livi, and A. Politi. Thermal conduction in classical low-dimensional lattices.Physics Reports, 377(1):1–80, 2003

    work page 2003

  11. [11]

    Lepri, R

    S. Lepri, R. Livi, and A. Politi. Studies of thermal conductivity in Fermi–Pasta–Ulam-like lattices.Chaos: An Interdisciplinary Journal of Nonlinear Science, 15(1), 2005

    work page 2005

  12. [12]

    K. Aoki, J. Lukkarinen, and H. Spohn. Energy transport in weakly anharmonic chains.Journal of Statistical Physics, 124(5):1105–1129, 2006

    work page 2006

  13. [13]

    A. Dhar. Heat transport in low-dimensional systems.Advances in Physics, 57(5):457–537, 2008

    work page 2008

  14. [14]

    Lukkarinen and H

    J. Lukkarinen and H. Spohn. Anomalous energy transport in the FPU-βchain.Communications on Pure and Applied Mathematics, 61(12):1753–1786, 2008

    work page 2008

  15. [15]

    Mej ´ıa-Monasterio, A

    C. Mej ´ıa-Monasterio, A. Politi, and L. Rondoni. Heat flux in one-dimensional systems.Physical Review E, 100(3):032139, 2019

    work page 2019

  16. [16]

    Boksenbojm, C

    E. Boksenbojm, C. Maes, K. Neto ˇcn´y, and J. Pe ˇsek. Heat capacity in nonequilibrium steady states.Europhysics Letters, 96(4):40001, 2011

    work page 2011

  17. [17]

    Pe ˇsek, E

    J. Pe ˇsek, E. Boksenbojm, and K. Netoˇcn´y. Model study on steady heat capacity in driven stochastic systems.Open Physics, 10(3):692– 701, 2012

    work page 2012

  18. [18]

    Pe ˇsek.Heat Processes in Non-Equilibrium Stochastic Systems

    J. Pe ˇsek.Heat Processes in Non-Equilibrium Stochastic Systems. PhD thesis, Charles University in Prague, 2014

    work page 2014

  19. [19]

    Maes and K

    C. Maes and K. Neto ˇcn´y. Nonequilibrium calorimetry.Journal of Statistical Mechanics: Theory and Experiment, 2019(11):114004, 2019

    work page 2019

  20. [20]

    Khodabandehlou, C

    F. Khodabandehlou, C. Maes, and K. Neto ˇcn´y. A Nernst heat theorem for nonequilibrium jump processes.The Journal of Chemical Physics, 158(20), 2023

    work page 2023

  21. [21]

    Kittel and P

    C. Kittel and P. McEuen.Introduction to Solid State Physics. John Wiley & Sons, 2018

    work page 2018

  22. [22]

    C. Maes, K. Neto ˇcn`y, and M. Verschuere. Heat conduction networks.Journal of Statistical Physics, 111(5):1219–1244, 2003

    work page 2003

  23. [23]

    Quasistatic response for nonequilibrium processes: evaluating the Berry potential and curvature

    A. Beyen, F. Khodabandehlou, and C. Maes. Quasistatic response for nonequilibrium processes: evaluating the Berry potential and curvature.arXiv, 2512.01654, 2025

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  24. [24]

    Khodabandehlou, C

    F. Khodabandehlou, C. Maes, I. Maes, and K. Netoˇcn´y. The vanishing of excess heat for nonequilibrium processes reaching zero ambient temperature.Annales Henri Poincar ´e, 2023

    work page 2023

  25. [25]

    Khodabandehlou, C

    F. Khodabandehlou, C. Maes, and K. Neto ˇcn´y. On the Poisson equation for nonreversible Markov jump processes.Journal of Mathe- matical Physics, 65(4), 2024

    work page 2024

  26. [26]

    Dhar and D

    A. Dhar and D. Roy. Heat transport in harmonic lattices.Journal of Statistical Physics, 125(4):801–820, 2006

    work page 2006

  27. [27]

    Santra and U

    I. Santra and U. Basu. Activity driven transport in harmonic chains.SciPost Physics, 13(2):041, 2022

    work page 2022

  28. [28]

    H. Spohn. Long range correlations for stochastic lattice gases in a nonequilibrium steady state.Journal of Physics A: Mathematical and General, 16(18):4275–4291, 1983

    work page 1983

  29. [29]

    Krekels, C

    S. Krekels, C. Maes, I. Santra, and R. Zhai. Negative differential heat conductivity in a harmonic chain coupled to a particle reservoir. arXiv, 2026

    work page 2026