Здесь уже совсем ничего непонятно.
Извиняюсь за долгий ответ, забыл про Ваше сообщение.
Попробую описать эту простую модель более подробно.
У Вас имеется куча молекул (

штуков), и каждая из них находится в состоянии

, где

-- амплитуда и фаза соответствующего вклада, причём амплитуды (

) для всех молекул одинаковы, а вот фазы -- могут быть разными. Такое организовать очень просто: например, у Вас много молекул в длинной трубке, Вы светите на молекулы лазером с длиной волны существенно меньше, чем размер трубки, в результате чего молекулы в разных кусках трубы будут "чувствовать" разную фазу электрического поля, что приведёт к разным фазам молекул в разных кусках трубки, при этом молекулы будут продолжать лететь в своих направлениях, в результате чего через некоторое время молекулы с разными фазами окажутся в одних и тех же кусках трубы.
Чтобы посчитать наблюдаемую величину

, нам произвести следующую операцию:

,
где

-- это матрица плотности чистого состояния молекулы номер
i с элементами

.
То бишь, матрица плотности ансамбля оказывается суммой матриц плотностей чистых состояний индивидуальных молекул.
Если у нас молекул очень много, мы можем заменить суммирование на интегрирование, и в ансамбле у нас найдётся молекула с любой разностью фаз

, причём каждая из этих разностей фаз будет равновероятна. В результате мы получаем, что диагональные элементы матрицы плотности будут

(заселённости даны амплитудами состояний), а вот недиагональные окажутся

.