It would seem to me that applying your thermostat to all atoms is at odds with your desire to heat a specific region.
To paraphrase Sytemantics: Complex systems almost never work. And those that do are built from simple systems that work.
Have you successfully built a simulation where you have defined and heated a specific static region?