С формальной точки зрения это безусловно не так. Пусть у вас имеется система невзаимодействующих магнитных моментов, двух типов (с разными g-факторами). Можно подобрать резонансную частоту так, чтобы фотоны взаимодейстовали с одной подсистемой и не взаимодействовали со второй. Тогда мы пробуем только одну подсистему, но вторая никоим образом термостатом для нее не является. Взаимодействия-то между подсистемами нет. Это не тот случай, который релевантен для статфизики, и не тот, что мы рассматриваем. Вот он и исключается требованием некоммутативности операторов Лиувилля для подсистемы и для всей системы.