Analytic solution for heat flow through a general harmonic network

Abstract

We present an analytic expression for the heat current through a general harmonic network coupled with Ohmic reservoirs. We use a method that enables us to express the stationary state of the network in terms of the eigenvectors and eigenvalues of a generalized cubic eigenvalue problem. In this way, we obtain exact formulas for the heat current and the local temperature inside the network. Our method does not rely on the usual assumptions of weak coupling to the environments or on the existence of an infinite cutoff in the environmental spectral densities. We use this method to study nonequilibrium processes without the weak coupling and Markovian approximations. As a first application of our method, we revisit the problem of heat conduction in two- and three-dimensional crystals with binary mass disorder. We complement previous results showing that for small systems the scaling of the heat current with the system size greatly depends on the strength of the interaction between system and reservoirs. This somewhat counterintuitive result seems not to have been noticed before.