Nonequilibrium stationary state of a harmonic crystal with alternating masses

We analyze the nonequilibrium steady states (NESS) of a one-dimensional harmonic chain of N atoms with alternating masses connected to heat reservoirs at unequal temperatures. We find that the temperature profile defined through the local kinetic energy T(j)≡〈p_j²〉/m_j oscillates with period two in the bulk of the system. Depending on boundary conditions, either the heavier or the lighter particles in the bulk are hotter. We obtain explicit integral expressions for the bulk temperature profile and steady state current in the limit N→∞. These depend on whether N is odd or even. We also study similar temperature oscillations in the NESS of systems with noise in the dynamics. These die out as N→∞