Negative Differential Heat Conductivity in a Harmonic Chain Coupled to a Particle Reservoir

Abstract: When coupling thermal baths at different temperatures, negative differential thermal conductivity is typically attributed to nonlinear interactions in the connecting medium. In this work, we demonstrate that such an effect can arise purely from the nature of the thermal baths and their coupling with the medium. Specifically, we construct a bath composed of overdamped thermal particles, which is coupled to one end of a harmonic chain, while the other end is connected to a standard Langevin heat bath. By analyzing the steady-state heat current, we observe significant negative differential thermal conductivity. In particular, as the temperature difference between the two baths diverges, the steady-state heat current through the chain vanishes. The effect is thermokinetic: we compute the effective dissipative coefficient and we find that it scales inversely with the square of the temperature of the particle bath in the high-temperature limit, resulting in an asymptotic decoupling between the bath and the chain. Our results highlight that nonequilibrium transport properties can be strongly influenced by the structure of the environment and its coupling to the system, even in otherwise linear systems.

• The paper presents a finding that increasing the left bath temperature reduces the heat current in an otherwise linear harmonic chain.
• It employs a unique thermokinetic coupling methodology that integrates out bath dynamics via a von Mises-type periodic function.
• The study highlights potential for tuning thermal transport properties in mesoscopic systems through strategic reservoir engineering.

Negative Differential Heat Conductivity from Thermokinetic Bath Coupling in Harmonic Chains

Introduction and Background

This work examines energy transport in a driven harmonic chain coupled to two distinct thermal environments: one end is connected to a standard Langevin heat bath, and the other to a reservoir of overdamped thermal particles. Crucially, the chain itself is strictly linear—no nonlinear bulk interactions are present. The study is motivated by the persistent question in nonequilibrium statistical mechanics concerning how the nature and coupling mechanisms of thermal reservoirs affect macroscopic transport, especially in low-dimensional or soft-matter settings where positional and elastic degrees of freedom coexist.

Prior results in the literature demonstrate that negative differential conductivity (NDC)—where increasing the temperature bias leads to diminished current—typically arises in systems with nonlinear interactions in the bulk or with active/kinetically-constrained boundary conditions. The novelty here is the emergence of NDC in a fully harmonic chain, not from nonlinearity within the chain, but from the unique, thermally driven, non-Langevin coupling at the boundary.

Model Formulation and Analysis

The model consists of a one-dimensional chain of $N$ coupled harmonic oscillators, each with positions $\{q_j\}$ and momenta $\{p_j\}$. The rightmost oscillator couples to a conventional Langevin reservoir at temperature $T_\mathrm{R}$. The leftmost oscillator interacts locally, via a von Mises-type periodic function, with a bath composed of $N_a$ overdamped Brownian particles at temperature $T_L$. The Hamiltonian includes standard quadratic terms for the chain and an interaction term capturing the effect of particle "bombardment".

A careful reduction is performed to integrate out the fast bath variables under a regime of time-scale separation. The left oscillator dynamics acquires three contributions from the bath: a (quasi-static) mean force, a fluctuating force exhibiting fluctuation-dissipation properties, and—of central importance—an effective friction term. This friction coefficient $\gamma_\mathrm{eff}$ grows with the square of the particle bath temperature: for high $T_L$, $\gamma_\mathrm{eff} \propto 1/T_L^2$. This scaling implies that, as the left bath temperature increases, the energy transfer from chain to bath becomes less efficient due to an emergent decoupling at high $T_L$. Simulation results for the single oscillator confirm this analytic scaling.

With this effective description, the full chain dynamics is recast in terms of temperature-dependent effective friction and stiffness at the left boundary.

Heat Transport and Emergence of NDC

The steady-state heat current is evaluated using nonequilibrium Green's function methodology. For canonical harmonic chains coupled at both boundaries to Langevin reservoirs, current is size-independent and grows linearly with temperature bias, as given by Rieder-Lieb-Lebowitz theory.

In contrast, here, the key finding is the highly non-monotonic dependence of steady-state current $\{q_j\}$0 on $\{q_j\}$1. Specifically, for fixed $\{q_j\}$2, the stationary current $\{q_j\}$3 initially increases with $\{q_j\}$4, attains a maximum, and then decreases rapidly, vanishing as $\{q_j\}$5 for large $\{q_j\}$6. This regime constitutes clear negative differential heat conductivity—an increase in driving temperature $\{q_j\}$7 reduces current, contrary to standard intuitive expectations from Fourier-type transport laws. The mechanism is purely thermokinetic: it arises from the temperature dependence of the effective friction mediated by the non-Langevin particle bath, absent any nonlinearities in the chain.

Numerically, the authors confirm the predicted scaling and the pronounced non-monotonicity for various parameters. They demonstrate that varying the standard Langevin bath temperature $\{q_j\}$8 yields conventional linear behavior, further isolating the effect to the unique nature and coupling of the left particle bath.

Theoretical and Practical Implications

The demonstration that NDC can arise solely from reservoir engineering in an otherwise linear system revises the understanding of boundary-driven transport far from equilibrium. This challenges the common paradigm that attributes all such phenomena to nonlinearity or strong kinetic constraints in either the bulk or reservoirs. The work underscores the need to model bath-system coupling explicitly, especially in active matter, viscoelastic networks, and synthetic or biological systems where environmental degrees of freedom do not merely serve as featureless noise sources.

Theoretically, this research paves the way for a more general framework for open-system transport phenomena in which the kinetic origin and temperature scaling of bath-induced dissipation must be accounted for. Practically, it raises the possibility of tuning thermal transport properties in mesoscopic systems via bath engineering, potentially offering designs for thermal regulators or rectifiers in nano/biotechnological applications without recourse to nonlinearity or active control.

Future Directions

Potential future research includes extending this analysis to higher dimensions, complex networked chains, multiphysics couplings (e.g., electrothermal analogues), or baths with internal structure and interacting degrees of freedom. It would be particularly instructive to explore quantum analogues, the role of finite system sizes, and experimental realizations in colloidal and metamaterial platforms. Furthermore, investigating feedback-controlled baths or baths with memory (non-Markovian dynamics) could yield even richer nonequilibrium transport phenomenology.

Conclusion

This study establishes that negative differential heat conductivity can emerge in a strictly linear harmonic chain due purely to the thermokinetic structure of the boundary reservoir, specifically when coupled to an overdamped particle bath whose effective dissipative coupling scales inversely with the square of the bath temperature. This finding highlights how boundary conditions and reservoir-system couplings can fundamentally reconfigure transport properties, and motivates a reevaluation of theoretical approaches to nonequilibrium processes in low-dimensional and hybrid systems.