Квантовая статистическая механика

Статистическая механика квантово-механических систем

Квантовая статистическая механика — это статистическая механика, применяемая к квантово-механическим системам . В квантовой механике статистический ансамбль ( распределение вероятностей по возможным квантовым состояниям ) описывается оператором плотности S , который является неотрицательным, самосопряженным , следовым оператором следа 1 в гильбертовом пространстве H, описывающем квантовую систему. Это можно показать с помощью различных математических формализмов для квантовой механики .

Ожидание

Из классической теории вероятностей мы знаем, что ожидание случайной величины X определяется ее распределением D X по формуле

E ( X ) = R d λ D X ( λ ) {\displaystyle \mathbb {E} (X)=\int _{\mathbb {R} }d\lambda \operatorname {D} _{X}(\lambda )}

предполагая, конечно, что случайная величина интегрируема или что случайная величина неотрицательна. Аналогично, пусть A будет наблюдаемой квантово-механической системой. A задается плотно определенным самосопряженным оператором на H . Спектральная мера A определяется как

E A ( U ) = U d λ E ( λ ) , {\displaystyle \operatorname {E} _{A}(U)=\int _{U}d\lambda \operatorname {E} (\lambda ),}

однозначно определяет A и наоборот, однозначно определяется A. E Aбулев гомоморфизм из борелевских подмножеств R в решетку Q самосопряженных проекций H. По аналогии с теорией вероятностей, если задано состояние S , мы вводим распределение A относительно S , которое является вероятностной мерой , определенной на борелевских подмножествах R с помощью

D A ( U ) = Tr ( E A ( U ) S ) . {\displaystyle \operatorname {D} _{A}(U)=\operatorname {Tr} (\operatorname {E} _{A}(U)S).}

Аналогично, ожидаемое значение A определяется с точки зрения распределения вероятностей D A следующим образом:

E ( A ) = R d λ D A ( λ ) . {\displaystyle \mathbb {E} (A)=\int _{\mathbb {R} }d\lambda \,\operatorname {D} _{A}(\lambda ).}

Обратите внимание, что это ожидание относится к смешанному состоянию S , которое используется в определении D A.

Замечание . По техническим причинам необходимо рассматривать отдельно положительные и отрицательные части A, определяемые функциональным исчислением Бореля для неограниченных операторов.

Можно легко показать:

E ( A ) = Tr ( A S ) = Tr ( S A ) . {\displaystyle \mathbb {E} (A)=\operatorname {Tr} (AS)=\operatorname {Tr} (SA).}

След оператора А записывается следующим образом:

Tr ( A ) = m m | A | m . {\displaystyle \operatorname {Tr} (A)=\sum _{m}\langle m|A|m\rangle .}

Обратите внимание, что если Sчистое состояние, соответствующее вектору , то: ψ {\displaystyle \psi }

E ( A ) = ψ | A | ψ . {\displaystyle \mathbb {E} (A)=\langle \psi |A|\psi \rangle .}

Энтропия фон Неймана

Особое значение для описания случайности состояния имеет энтропия фон Неймана S, формально определяемая как

H ( S ) = Tr ( S log 2 S ) {\displaystyle \operatorname {H} (S)=-\operatorname {Tr} (S\log _{2}S)} .

На самом деле, оператор S log 2 S не обязательно является следовым классом. Однако, если S — неотрицательный самосопряженный оператор, не являющийся следовым классом, мы определяем Tr( S ) = +∞. Также отметим, что любой оператор плотности S может быть диагонализирован, то есть его можно представить в некотором ортонормированном базисе (возможно, бесконечной) матрицей вида

[ λ 1 0 0 0 λ 2 0 0 0 λ n ] {\displaystyle {\begin{bmatrix}\lambda _{1}&0&\cdots &0&\cdots \\0&\lambda _{2}&\cdots &0&\cdots \\\vdots &\vdots &\ddots &\\0&0&&\lambda _{n}&\\\vdots &\vdots &&&\ddots \end{bmatrix}}}

и мы определяем

H ( S ) = i λ i log 2 λ i . {\displaystyle \operatorname {H} (S)=-\sum _{i}\lambda _{i}\log _{2}\lambda _{i}.}

Соглашение заключается в том, что , поскольку событие с вероятностью ноль не должно вносить вклад в энтропию. Это значение является расширенным действительным числом (то есть находится в [0, ∞]), и это, очевидно, унитарный инвариант S . 0 log 2 0 = 0 {\displaystyle \;0\log _{2}0=0}

Замечание . Действительно возможно, что H( S ) = +∞ для некоторого оператора плотности S . Фактически T — диагональная матрица

T = [ 1 2 ( log 2 2 ) 2 0 0 0 1 3 ( log 2 3 ) 2 0 0 0 1 n ( log 2 n ) 2 ] {\displaystyle T={\begin{bmatrix}{\frac {1}{2(\log _{2}2)^{2}}}&0&\cdots &0&\cdots \\0&{\frac {1}{3(\log _{2}3)^{2}}}&\cdots &0&\cdots \\\vdots &\vdots &\ddots &\\0&0&&{\frac {1}{n(\log _{2}n)^{2}}}&\\\vdots &\vdots &&&\ddots \end{bmatrix}}}

T — неотрицательный следовой класс, и можно показать, что T log 2 T не является следовым классом.

Теорема . Энтропия — унитарный инвариант.

По аналогии с классической энтропией (обратите внимание на сходство в определениях), H( S ) измеряет количество случайности в состоянии S . Чем более разбросаны собственные значения, тем больше энтропия системы. Для системы, в которой пространство H конечномерно, энтропия максимальна для состояний S , которые в диагональной форме имеют представление

[ 1 n 0 0 0 1 n 0 0 0 1 n ] {\displaystyle {\begin{bmatrix}{\frac {1}{n}}&0&\cdots &0\\0&{\frac {1}{n}}&\dots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &{\frac {1}{n}}\end{bmatrix}}}

Для такого S , H( S ) = log 2 n . Состояние S называется максимально смешанным состоянием.

Напомним, что чистое состояние — это одна из форм

S = | ψ ψ | , {\displaystyle S=|\psi \rangle \langle \psi |,}

для ψ — вектор нормы 1.

Теорема . H( S ) = 0 тогда и только тогда, когда S является чистым состоянием.

Ибо S является чистым состоянием тогда и только тогда, когда его диагональная форма имеет ровно один ненулевой элемент, равный 1.

Энтропию можно использовать как меру квантовой запутанности .

Канонический ансамбль Гиббса

Рассмотрим ансамбль систем, описываемых гамильтонианом H со средней энергией E. Если H имеет чисто точечный спектр и собственные значения H стремятся к +∞ достаточно быстро, e r H будет неотрицательным оператором следового класса для каждого положительного r . E n {\displaystyle E_{n}}

Канонический ансамбль Гиббса описывается состоянием

S = e β H Tr ( e β H ) . {\displaystyle S={\frac {\mathrm {e} ^{-\beta H}}{\operatorname {Tr} (\mathrm {e} ^{-\beta H})}}.}

Где β таково, что среднее значение энергии по ансамблю удовлетворяет условию

Tr ( S H ) = E {\displaystyle \operatorname {Tr} (SH)=E}

и

Tr ( e β H ) = n e β E n = Z ( β ) {\displaystyle \operatorname {Tr} (\mathrm {e} ^{-\beta H})=\sum _{n}\mathrm {e} ^{-\beta E_{n}}=Z(\beta )}

Это называется функцией распределения ; это квантово-механическая версия канонической функции распределения классической статистической механики. Вероятность того, что система, выбранная случайным образом из ансамбля, будет находиться в состоянии, соответствующем собственному значению энергии, равна E m {\displaystyle E_{m}}

P ( E m ) = e β E m n e β E n . {\displaystyle {\mathcal {P}}(E_{m})={\frac {\mathrm {e} ^{-\beta E_{m}}}{\sum _{n}\mathrm {e} ^{-\beta E_{n}}}}.}

При определенных условиях канонический ансамбль Гиббса максимизирует энтропию фон Неймана состояния, подчиняющегося требованию сохранения энергии. [ необходимо разъяснение ]

Большой канонический ансамбль

Для открытых систем, где энергия и число частиц могут колебаться, система описывается большим каноническим ансамблем , описываемым матрицей плотности

ρ = e β ( i μ i N i H ) Tr ( e β ( i μ i N i H ) ) . {\displaystyle \rho ={\frac {\mathrm {e} ^{\beta (\sum _{i}\mu _{i}N_{i}-H)}}{\operatorname {Tr} \left(\mathrm {e} ^{\beta (\sum _{i}\mu _{i}N_{i}-H)}\right)}}.}

где N 1 , N 2 , ... — операторы числа частиц для различных видов частиц, которые обмениваются с резервуаром. Обратите внимание, что это матрица плотности, включающая гораздо больше состояний (с переменным N) по сравнению с каноническим ансамблем.

Великая функция распределения — это

Z ( β , μ 1 , μ 2 , ) = Tr ( e β ( i μ i N i H ) ) {\displaystyle {\mathcal {Z}}(\beta ,\mu _{1},\mu _{2},\cdots )=\operatorname {Tr} (\mathrm {e} ^{\beta (\sum _{i}\mu _{i}N_{i}-H)})}

Смотрите также

Ссылки

Retrieved from "https://en.wikipedia.org/w/index.php?title=Quantum_statistical_mechanics&oldid=1251270348"