Методичка по осциллятору, квантованию, матрице плотности (1129351), страница 9
Текст из файла (страница 9)
Поскольку λ не зависит от состояния, энтропия с полученными вероятностямимаксимальна.Вспомним теперь основные свойства матрицы плотности, рассмотренные в параграфе (3.2), а именно, свойство74(7.7) и (7.8). Следовательно можно выразить энтропию согласно определению (4.1) через матрицу плотности в диагональном представлении:X(4.4)S=−ρν lnρν .νОбратим внимание, что выражение (4.4) можно переписать, используя выражение (1.6):XS=−hν 0 |(ρν |νihν|)|ν 00 ihν 00 |(lnρν |νihν|)|ν 0 i = −T rρ̂ν lnρ̂ν .ν,ν 0 ,ν 00Здесь ρ̂ν обозначает матрицу плотности, записанную в диагональном представлении.Поскольку от диагонального представления всегда можно перейти к произвольному, запишем теперь определениеэнтропии через оператор ρ̂ в общем виде:S = −T rρ̂lnρ̂.(4.5)В дальнейшем будем использовать термины как матрицаплотности, так и статистический оператор.Определим теперь вид статистического оператора в рассмотренном выше случае, когда все состояния ансамбля систем равновероятны.
Проварьируем определение (4.5) приусловии равенства единице следа оператора:T r(1 + lnρ̂ + λ)δ ρ̂ = 0,или1− c − число.(4.6)1+λЗададим теперь дополнительные сведения об ансамблесистем. А именно, пусть ансамбль характеризуется энергией, которая согласно свойствам матрицы плотности поопределению естьρ̂ =bhEi = E = T rρ̂H,75(4.7)b – гамильтониан систем ансамбля.где HВновь потребуем максимума энтропии, но теперь ещепри одном дополнительном условии (4.7):b ρ̂ = 0.T r(1 + lnρ̂ + λ + β H)δПоскольку все вариации произвольны, получаемилиblnρ̂ = −1 − λ − β H,bρ̂ = e−(1+λ) e−β H .(4.8)Первая экспонента в формуле (4.8) может быть выражена через статистическую сумму из условия нормировкиматрицы плотности T rρ̂ = 1 :be1+λ T re−β H ≡ Z.(4.9)Можно теперь переписать выражение (4.8), используя определение (4.9):bρ̂ = Z −1 e−β H .(4.10)Оставшийся неопределенный параметр β находится из условия (4.7):bb −β H = −hEi = Z −1 T rHe∂lnZ.∂β(4.11)Таким образом, параметр β определяется величиной средней энергии.
В термодинамическом пределеβ=1.T(4.12)Здесь и далее мы измеряем температуру в единицах энергии ( или полагаем постоянную Больцмана k = 1.)764.5Уравнение для матрицы плотностив координатном представленииЗапишем уравнение Лиувилля (3.3) для матрицы плотности в координатном представлении, воспользовавшись методами теории представлений, изложенными в главе 11 :0b ρ̂(t)|x0 i − ihx|ρ̂(t)H|xb 0 i = 0.ˆhx|ρ̇(t)|xi + ihx|H(5.1)Обозначим для краткости матричные элементы в координатном представлении как:b 0 i ≡ Axx0 .hx|A|x(5.2)В обозначениях (7.5) уравнение (5.1) принимает вид:XX∂ρxx0 (t) + iHxy ρyx0 − iρxy Hyx0 = 0∂tyy(5.3)Для непрерывного базиса все суммы имеют смысл интегральных и поэтому заменяются интегралами, а матричные элементы представляются ядрами этих интегральныхоператоров.
Следовательно, можно записать в координатном представлении:ZXbH(x)ρ(x, y) ≡ H(x, x0 )ρ(x0 , y)dx0 =Hxx0 ρx0 y . (5.4)x0Иными словами, можно суммирование обозначить какXbHxx0 ρx0 y = H(x)ρ(x,y).(5.5)x01В этом параграфе мы полагаем все константы единицами: ~ =1, m = 1 и т.д.77Совершенно аналогично записывается это соотношение в“обратном” порядке:Xb tr (x0 )ρ(x, x0 , t).ρxy Hyx0 = H(5.6)yТаким образом, уравнение Лиувилля в координатном представлении (5.3) принимает вид:∂bb tr (x0 )ρ(x, x0 , t) = 0 (5.7)ρ(x, x0 , t) + iH(x)ρ(x,x0 , t) − iH∂tГамильтониан частицы в координатном представлении имеет дифференциальную форму:p̂21 ∂2bH(x)=+ U (x) = −+ U (x)22 ∂x2(5.8)и в силу эрмитовостиb tr (x) = H(x),bHпоэтому уравнение (5.7) можно переписать в форме дифференциального уравнения:µ 2¶∂∂i∂20ρ(x, x0 , t)+ρ(x, x , t)−−∂t2 ∂x0 2 ∂x2¡¢+ U (x) − U (x0 ) ρ(x, x0 , t) = 0.(5.9)78Глава 4Представления матрицыплотности4.1Функция ВигнераВ предыдущих главах мы видели, что между классическими квантовым описанием систем существует принципиальное отличие: вместо физических величин, задаваемых числовыми функциями, определены операторы, подчиняющиеся определенной алгебре.
Физические величины определяются средними значениями соответствующих операторов.Принципиальное отличие заключается в том, что операторы между собой не всегда коммутируют, что означаетневозможность одновременной измеримости соответствующих физических величин. В классической механике такихпроблем нет. Однако в 1932 году Е. Вигнер открыл такоепредставление в квантовой механике, которое оказываетсянаиболее близким к классическим представлениям.
ПодходВигнера затем в 1949 году был более полно развит Дж.Мойалем. Такое представление задается функцией Вигнера W (q, p), где q и p – обобщенные координаты и импульс79частицы (физические величины).Запишем (одночастичную) матрицу плотности в координатном представлении, считая спектр дискретным:Xρ(r, r0 ) =wn ϕn (r)ϕ∗n (r0 ).(1.1)nКак помним, диагональный элемент матрицы плотности вкоординатном представлении определяет вероятность обнаружить частицу в точке с координатой r).Совершенно аналогично в импульсном представлениидиагональный элемент ρ(p, p) определяет вероятность обнаружить частицу с значением импульса p.
Установим связьмежду этими величинами.Запишем матрицу плотности в импульсном представлении:Xρ(p, p0 ) =wn ϕn (p)ϕ∗n (p0 ) =n=Zdrdr0 ρ(r, r0 ) expµ¶i(pr − p0 r0 ) .~(1.2)Матрица плотности в координатном представлении удобна для вычисления средних значений функций координат:Rdrρ(r, r)V (r)T rVb ρ.(1.3)hV (r)i == RT rρdrρ(r, r)Импульсное представление удобно для вычисления средних значений функций импульсов (например, кинетическойэнергии):Rcρ(dp/(2π~)3 )ρ(p, p)W (p)T rWRhW (p)i ==.(1.4)T rρ(dp/(2π~)3 )ρ(p, p)В классической механике такая задача решалась бы с помощью одной функции распределения в фазовом пространстве F (p, r), которая имеет смысл плотности вероятности.80Причем распределение вероятности в координатном пространстве определяется просто частичным интегрированием по импульсному пространству:Zdpf (r) = F (p, r),(1.5)(2π~)3а распределение вероятности в импульсном пространстве,соответственно, интегрированием по координатной частифазового пространства:Zf (p) = F (p, r)dr.(1.6)Возникает естественное желание построить и в квантовой механике такую функцию распределения, котораябы удовлетворяла условиям (1.5) и (7.2).
Для этого рассмотрим диагональные компоненты матрицы плотности вимпульсном представлении (1.2), которые определяют вероятности и должны соответствовать функции распределения. Поскольку искомая функция должна одновременно зависеть как от импульсов, так и от координат, проведем “частичное” преобразование Фурье в формуле (1.2),рассматривая формально матрицу плотности в координатном представлении как функцию переменных (x + x0 )/2 и(x − x0 )/2. Переходя затем к новым обозначениям x/2 →r, x0 → u, и оставляя преобразование только по переменной u, получаем функцию Вигнера:Z³uu ´ −i~−1 pue.(1.7)W (p, r) = duρ r + , r −22Формула (1.7) обратима, и для матрицы плотности можнозаписать выражение через функцию Вигнера:¶µZdpr + r0−100ρ(r, r ) =ei~ p(r−r ) .(1.8)W p,3(2π~)281Упражнение.
Получить обратное преобразование (1.8).Убедимся, что функция Вигнера удовлетворяет необходимым условиям. Прежде всего проинтегрируем по импульсу:µZ¶ZZ³dpuu´du −i~−1 puW (p, r) = duρ r + , r −e=(2π~)322(2π~)3Z³u´uδ(u) = ρ(r, r) = f (r).= duρ r + , r −22Таким образом, получили, что первое условие (1.5) выполняется. Для проверки выполнения второго условия (7.2)проделаем выкладки в “обратном” порядке:µ¶Zif (p) =ρ(p, p) = drdr0 ρ(r, r0 ) exp − p(r − r0 ) =~ZZ³´uu −i~−1 pu= dxduρ x + , x −e= dxW (p, x).22При преобразовании интеграла мы сделали замену переменных:uur = x + , r0 = x − .22Легко видеть, что введенная таким образом функция Вигнера W (p, r) позволяет получать средние значения величин, зависящих только либо от координат, либо от импульсов.
Нельзя получить правильный результат для функций,зависящих одновременно и от импульсов, и от координат.В этом нет ничего удивительного, поскольку координатыи импульс в квантовой механике связаны соотношениемнеопределенностей. Рассмотрим этот вопрос более подробно.824.2Некоторые свойства преобразованияФурье. Уравнение Мойала для функции ВигнераДальнейшее изложение материала в той или иной степенибудет связано с преобразованием Фурье матрицы плотности, а также соответствующих операторов.
Поэтому удобновспомнить некоторые полезные соотношения. Для простоты будем рассматривать одномерное преобразование. Намногомерный случай все результаты можно легко распространить.Пусть дано преобразование Фурье некоторой функцииf (x), т.е. 1Z1f (x) = √f˜(k)eikx dk(2.1)2πи, соответственно,Z1˜f (x) = √f (x)e−ikx dx.(2.2)2πПродифференцируем выражение (7.3) по координате и получим:Z hi1∂f (x) = √ik f˜(k) eikx dk.∂x2πИными словами, мы можем формально установить соответствие:∂f (x) −→ ik f˜(k).(2.3)∂xУмножим теперь выражение (7.3) на x :Z1f˜(k)xeikx dk.(2.4)xf (x) = √2π1Мы здесь применяем так называемое симметричное преобразование Фурье, хотя в физике пре переходе от координатного к импульсному представлению чаще всего используется несимметричноепреобразование.83Поскольку можно представитьxeikx = −i∂ ³ ikx ´,e∂kвыражение (2.4) можно переписать ка굶Z∂1˜f (k) −ixf (x) = √eikx dk.∂k2π(2.5)Беря интеграл по частям и считая значение подынтегральной функции на бесконечно удаленных пределах равнымнулю, получаем:¶Z µ∂1i f˜(k) eikx dk.(2.6)xf (x) = √∂k2πСледовательно, имеет место соответствие:xf (x) −→ i∂ ˜f (k).∂k(2.7)Можно сказать, если имеется некоторое выражени嵶∂f (x),(2.8)Φ x,∂xего образ Фурье дается выражение춵∂Φ i , ik f˜(k).∂k(2.9)Применим теперь полученные соотношения для матрицы плотности в координатном представлении и функцииВигнера.
Действительно, как следует из формул (1.7) и(1.8), рассматриваемые функции связаны между собой преобразованием Фурье. Запишем указанные соотношения в84одномерном случае в безразмерной форме, как и все остальное в этом параграфе, в момент времени t = 0 :Z ³u ´ −ipuuedu(2.10)W (q, p) = ρ q + , q −22и обратное преобразование:¶µZx + x000, p eip(x−x ) dp.ρ(x, x ) = W2(2.11)Умножим соотношение (2.11)на x и занесем эту переменную под знак интеграла:µ¶Zx + x000xρ(x, x ) = W, p xeip(x−x ) dp.(2.12)2Возникает соблазн поступить так же, как и в предыдущемслучае: продифференцировать по p а затем проинтегрировать по частям.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
