Fourier's law for a harmonic crystal with self-consistent stochastic reservoirs

We consider a d-dimensional harmonic crystal in contact with a stochastic Langevin type heat bath at each site. The temperatures of the "exterior" left and right heat baths are at specified values T-L and T-R, respectively, while the temperatures of the "interior" baths are chosen self-consistently so that there is no average flux of energy between them and the system in the steady state. We prove that this requirement uniquely fixes the temperatures and the self consistent system has a unique steady state. For the infinite system this state is one of local thermal equilibrium. The corresponding heat current satisfies Fourier's law with a finite positive thermal conductivity which can also be computed using the Green-Kubo formula. For the harmonic chain (d = 1) the conductivity agrees with the expression obtained by Bolsterli, Rich, and Visscher in 1970 who first studied this model. In the other limit, d >> 1, the stationary infinite volume heat conductivity behaves as (l(d)d)(-1) where l(d) is the coupling to the intermediate reservoirs. We also analyze the effect of having a non-uniform distribution of the heat bath couplings. These results are proven rigorously by controlling the behavior of the correlations in the thermodynamic limit.