We analytically compute the thermal conductivity of anharmonic crystals with self-consistent stochastic reservoirs. We develop an integral representation for the heat current, assume the approximation of discrete times, and in a perturbative analysis that is rigorously supported by a convergent cluster expansion compute the thermal conductivity for a chain with quartic anharmonic on-site potential. In the high anharmonicity regime, the result for the dependence of the conductivity on temperature is the same as for the system without inner reservoirs. The presented formalism is quite general and is also valid for inhomogeneous systems in any space dimension.

• Received 18 January 2013

DOI:https://doi.org/10.1103/PhysRevE.87.032158