
лекции по статам
.pdf
физической величины для системы из одинаковых бозонов можно выразить через операторы рождения и уничтожения.
Наиболее важные операторы динамических переменных для многочастичной системы имеют вид
|
|
|
N |
|
|
|
|
X |
|
|
|
^(1) |
X |
|
^(2) |
|
1 |
|
|
||
|
^(1) |
|
|
^(2) |
|
(4.7) |
||||
|
A |
= |
A |
(qi); |
A |
= |
2 |
A |
(qi; qj); |
|
|
|
|
i=1 |
|
|
|
|
i6=j |
|
|
|
|
|
|
|
|
|
|
|
|
|
^(1) |
|
|
|
|
|
|
|
|
^(2) |
(qi; qj) |
где A |
(qi) оператор, действующий на координаты i-ой частицы, A |
оператор, действующий на координаты частиц с номерами i и j. Динамическая пе-
^
ременная, оператор которой A(1) есть сумма операторов для отдельных частиц, называется аддитивной динамической переменной; динамические переменные
^
с операторами A(2) обычно называются динамическими переменными бинарного типа. В некоторых задачах квантовой механики встречаются операторы более сложной конструкции, но мы не будем здесь ими заниматься.
В качестве иллюстрации приведем простые примеры динамических переменных вида (??). Аддитивной динамической переменной является кинетическая энергия
^ |
2 |
=2m, а примером динамической переменной бинарного |
|||||
частиц системы T = |
i p^i |
||||||
|
|
^ |
|
U( ~r |
|
~r |
). Рекомен- |
энергия взаимодействия U = (1=2) |
i6=j |
|
|||||
типа может служить P |
|
|
j i |
jj |
|
P
дуем читателю самому вспомнить, какие из других ранее встречавшихся динамических переменных относятся к аддитивным переменным и к переменным бинарного типа.
^(1) |
^(2) |
|
Обычно операторы динамических переменных A |
и A |
легче всего построить |
в координатном представлении. Чтобы найти их матричные элементы в представлении чисел заполнения, нужно вычислить интегралы
0 |
^(1) |
jfnlgi = |
hfnlgjA |
||
0 |
^(2) |
jfnlgi = |
hfnlgjA |
|
|
i |
Z |
fnl0 g A (qi) fnlg dq1 dqN ; |
|||
X |
|
(s) |
^(1) |
(s) |
|||
|
|
|
|
|
|||
|
2 i=j |
Z |
|
|
(4.8) |
||
|
fnl0 g A (qi; qj) fnlg dq1 dqN |
||||||
1 |
X |
(s) |
^(2) |
(s) |
|||
|
|
|
6 |
|
|
|
|
с симметризованными базисными волновыми функциями (??). На первый взгляд задача кажется безнадежной из-за огромного числа переменных при больших значениях числа частиц N. Впрочем, есть и упрощающие обстоятельства. Так как
^ ^
частицы одинаковы, то вид операторов A(1)(qi) и A(2)(qi; qj) одинаков для любых
номеров частиц. Далее, многочастичные базисные волновые функции (fsn) g есть
l
произведения одночастичных ортонормированных волновых функций 'l(q). Поэтому удается получить более или менее простые выражения для матричных элементов. Мы не будем приводить соответствующие выкладки2. Оказывается, что
^(1) |
и любых операторов |
матричные элементы (??) любых аддитивных операторов A |
2Интересующийся читатель может обратиться к учебникам по квантовой механике (см., например, [?]).
3
^
бинарного типа A(2) совпадают с матричными элементами следующих операторов, записанных через операторы рождения и уничтожения:
A^(1) |
= Xl; l0 |
hl0jA^(1)jli a^ly0 a^l; |
|
|
|
||
^(2) |
|
1 |
l; |
X0 0 |
y |
y |
|
|
|
0 0 ^(2) |
|
||||
A |
= |
2 |
|
hl m jA |
jlmi a^m0 |
a^l0 |
a^l a^m: |
m; l ; m
Здесь h 0j ^(1)j i и h 0 0j ^(2)j i матричные элементы операторов ^(1) l A l l m A lm A
одночастичным волновым функциям 'l(q):
hl0jA^(1)jli = Z |
'l0 (q) A^(1)(q)'l(q) dq; |
|
hl0m0jA^(2)jlmi = Z |
'l0 (q1)'m0 (q2) A^(2)(q1; q2) 'l(q1)'m(q2) dq1 dq2: |
(4.9)
(4.10)
^
и A(2) по
(4.11)
(4.12)
Ясно, что матричные элементы (??) и (??) вычислить намного проще, чем матричные элементы (??) с многочастичными волновыми функциями.
Главное достоинство формул (??) и (??) состоит в том, что вычисление средних значений и матричных элементов операторов физических величин сводится теперь к вычислению средних значений и матричных элементов операторов, построенных из операторов рождения и уничтожения, которые довольно просто действуют на базисные векторы состояния системы jn1; n2; : : : ; nl; : : :i и удовлетворяют простым коммутационным соотношениям (??).
4

Статистическая физика
5.Представление чисел заполнения для фермионов (вторичное квантование)
5.1.Базисные квантовые состояния
В качестве базисных квантовых состояний для системы тождественных фермионов возьмем состояния, которые описываются антисимметричными волновыми функ-
циями (fan) g(q1; : : : ; qN ). Эти состояния полностью характеризуются набором чисел
l
заполнения nl. При этом, для каждого l, nl = 0; 1. Обозначим их
jfnlgi = jn1; n2; ; nl; i
Как известно, волновые функции этих состояний ортонормированы. В координатном представлении
(a) |
(q1; : : : ; qN ) = hq1; : : : ; qN jfnlgi: |
(5.1) |
fnlg |
Любой вектор состояния системы фермионов j (t)i может быть представлен в виде разложения по базисным векторам
Xf lg |
; : : : i; |
(5.2) |
j (t)i = C (fnlg; t) j : : : ; nl |
n
которое формально совпадает с формулой для бозонов, но в данном случае аргументом амплитуды вероятности C(fnlg; t) является последовательность нулей и единиц. Совокупность амплитуд вероятности волновая функция системы фермионов в представлении чисел заполнения.
5.2.Операторы рождения и уничтожения фермионов
Операторы рождения и уничтожения строятся так, чтобы формулы для операторов физических величин имели точно такой же вид, как и для бозонов (см. лекцию 4). Это удалось сделать П. Йордану и Е. Вигнеру в 1928 г.
Выпишем правила действия операторов рождения и уничтожения на базисные векторы состояния системы фермионов:
a^l j : : : ; nl; : : :i = ( 1) l pnl j : : : ; 1 nl; : : :i;
(5.3)
a^yl j : : : ; nl; : : :i = ( 1) l p1 nl j : : : ; 1 nl; : : :i:
l 1
X
Здесь l = nl число заполненных одночастичных состояний, предшеству-
k=1
ющих состоянию jli. Правила оказались более сложными, чем для бозонов. Впрочем, для практических вычислений сами эти правила используются очень
1

редко. Обычно достаточно знать основные свойства операторов рождения и уничтожения, которые мы теперь рассмотрим.
1. При действии операторов рождения и уничтожения на квантовые состояния системы не возникает нефизических состояний с числами заполнения больше единицы и меньше нуля.
Док-во: Согласно правилам действия операторов рождения и уничтожения,
a^l j : : : ; 0; : : :i = 0; |
a^l j : : : ; 1; : : :i = ( 1) l j : : : ; 0; : : :i; |
a^ly j : : : ; 0; : : :i = ( 1) l j : : : ; 1; : : :i; |
(5.4) |
a^ly j : : : ; 1; : : :i = 0: |
2. Квадраты операторов рождения и уничтожения равны нулю.
Док-во: Это проверяется с помощью предыдущих формул. Таким образом,
(^al)2 a^la^l = 0; (^aly)2 a^lya^ly = 0: |
(5.5) |
3. Оператор числа частиц в одночастичном состоянии jli имеет вид n^l = a^yl a^l (как и для бозонов).
Док-во: (Нужно помнить, что nl = 0; 1)
a^yl a^l j : : : ; nl; : : :i = ( 1) l pnl a^yl j : : : ; 1 nl; : : :i = nl j : : : ; nl; : : :i:
4. Коммутационные соотношения для операторов рождения и уничтожения.
Введем понятие антикоммутатора операторов
^ ^ ^ ^ ^ ^ |
антикоммутатор. |
(5.6) |
fA; Bg AB + BA |
Тогда из правил (??) следует, что
fa^l; a^l0 g = fa^ly; a^ly0 g = 0; |
fa^l; a^ly0 g = ll0 : |
(5.7) |
Они заменяют коммутационные соотношения для бозонов. В квантовой механике операторы рождения и уничтожения, удовлетворяющие соотношениям (??), принято называть ферми-операторами. С помощью ферми-операторов описывается, например, система электронов в кристаллах.
Док-во последнего соотношения в (??).
Пусть l 6= l0 и l < l0 (случай, когда l > l0 самим). Нужно показать, что
a^la^yl0 + a^yl0 a^ j ; nl; nl0 ; i = 0:
2

Рассмотрим
a^la^yl0 j ; nl; nl0 ; i = a^l( 1) l0 p1 nl0 j ; nl; ; 1 nl0 ; i
Снова применяя основное правило, получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
a^ a^y |
|
|
|
|
|
|
|
= ( 1) l0 |
|
|
1) l p |
|
|
p |
|
|
|
|
|
|
|
|
|
|
|
|
|
(5.8) |
|
; n ; |
|
n |
|
|
; |
i |
( |
|
|
|
1 |
|
nl |
|
j |
; 1 |
|
nl; |
|
; 1 |
|
nl ; |
i |
||||||
|
0 |
n |
|
|
|||||||||||||||||||||||||
l |
l0 j |
l |
|
l |
|
|
|
|
|
l |
|
0 |
|
|
|
0 |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Действуя теперь операторами в обратном порядке, получим
a^yl0 a^lj ; nl; ; nl0 ; i = a^yl0 ( 1) l pnl j ; 1 nl; ; nl0 ; i
Если nl = 0, то результат действия будет 0 как и в (??). Пусть nl = 1. Тогда число заполненных одночастичных состояний, предшествующих состоянию jli стало на единицу меньше. Поэтому
a^yl0 a^lj ; nl; ; nl0 ; i = ( 1) l0 ( 1) l pnl p1 nl0 j ; 1 nl; ; 1 nl0 ; i
(5.9) Сравнивая это равенство с (??), видим, что результаты действия на любое базисное состояние системы a^yl0 a^l и a^la^yl0 отличаются только знаком.
Пусть теперь l = l0. Тогда |
|
|
|
|
a^la^ly j : : : ; nl; : : :i = ( 1) l |
p |
|
a^lj : : : ; 1 nl; : : :i = (1 nl)j : : : ; nl; : : :i: |
|
1 nl |
||||
Итак, |
|
|
|
|
a^la^ly j : : : ; nl; : : :i = (1 nl) j : : : ; nl; : : :i: |
(5.10) |
Сравнивая это с
a^yl a^l j : : : ; nl; : : :i = nl j : : : ; nl; : : :i;
видим, что результат действия a^la^yl + a^yl a^l на любое базисное состояние системы эквивалентно действию единичного оператора.
Резюме:
Напомним, что операторы рождения и уничтожения для фермионов построены так, что формулы для операторов физических величин в представлении чисел заполнения имеют точно такой же вид, что и для бозонов. Таким образом, различие между статистиками Бозе-Эйнштейна и Ферми-Дирака описывается лишь коммутационными соотношениями между операторами рождения и уничтожения.
Пример из алгебры операторов рождения и уничтожения.
Вычислим коммутатор [^al; n^l0 ] сразу для двух статистик.
Пишем
[^al; n^l0 ] = a^ln^l0 n^l0 a^l = a^la^yl0 a^l0 a^yl0 a^l0 a^l
3
Далее используется такой прием: в каждом слагаемом все операторы рождения нужно расположить слева от всех операторов уничтожения, используя коммутационные соотношения
a^la^l0 a^l0 a^l = 0; |
a^lya^ly0 a^ly0 a^ly = 0; |
a^la^ly0 a^ly0 a^l = ll0 ; |
где верхний знак относится к бозонам, а нижний к фермионам. Согласно по-
следнему равенству,
a^la^yl0 = ll0 a^yl0 a^l
Возвращаясь к нашему примеру, получим
[^al; n^l0 ] = ll0 a^yl0 a^l a^l0 a^yl0 a^l0 a^l = ll0 a^l a^yl0 a^la^l0 a^yl0 a^l0 a^l
Два последних члена, где операторы расставлены в нужном порядке, точно сокращаются (при перестановке операторов уничтожения для фермионов знак еще раз изменится). Поэтому для обеих статистик
[^al; n^l0 ] = ll0 a^l
Самим так же показать, что
[^ayl ; n^l0 ] = ll0 a^yl
Ответ, конечно, можно получить быстрее, если вспомнить, что операторы рождения и уничтожения являются эрмитово сопряженными друг к другу.
4
Статистическая физика
6.Идеальные квантовые газы
Идеальный квантовый газ модель, в которой частицы не взаимодействуют друг с другом. В зависимости от типа статистики, которой подчиняются частицы, существуют две модели бозе-газ и ферми-газ.
Примеры: газообразный He4 бозе-газ, газообразный He3 ферми-газ, газ свободных электронов ферми-газ.
6.1.Гамильтониан идеального квантового газа
^
Обозначение: hi гамильтониан одной частицы с номером i. Гамильтониан всего газа
|
N |
^ |
Xi |
^ |
|
H = |
hi |
|
=1 |
^
Если частицы одинаковы, то все hi имеют одинаковый вид. Простейший пример: частицы без структуры
^ |
|
p^i2 |
hi |
= |
|
|
||
|
2m |
(6.1)
(6.2)
Более сложная модель молекулярный идеальный газ. В этом случае
^ гамильтониан молекулы, включающий операторы кинетической энергии ядер hi
и электронов, а также операторы их взаимодействия (главное взаимодействие кулоновское).
В представлении чисел заполнения |
|
|
H^ = Xl; l0 |
hl0jh^jli a^ly0 a^l |
(6.3) |
где jli базисные квантовые состояния одной частицы.
Часто бывает удобно выбрать в качестве jli стационарные состояния частицы, удовлетворяющие уравнению
^ |
jli |
(6.4) |
h jli = "l |
где l набор всех квантовых чисел, определяющих стационарное состояние частицы, а "l уровни энергии частицы. Тогда
H^ = X"l a^lya^l |
= X"l n^l |
(6.5) |
l |
l |
|
При таком выборе одночастичных состояний квантовые состояния jfnkgi являются собственными состояниями гамильтониана системы (т.е. стационарными состояниями). Проверка элементарна.
1

Уровни энергии системы
XX
Efnlg = "l nl; |
nl = N: |
(6.6) |
l |
l |
|
Полученный результат физически очевиден.
6.2.Основное состояние квантового идеального газа
Бозе-газ
Eосн = "min N |
(6.7) |
Ферми-газ
В основном состоянии частицы занимают одночастичные состояния jli с минимально возможной энергией. До некоторого значения энергии "F все состояния будут заполнены, а состояния с энергией "l > "F будут свободны. Максимальное значение энергии одночастичных занятых состояний "F называется энергией Ферми.
Таким образом, энергия основного состояния идеального ферми-газа
XX
Eосн = |
"l; |
1 = N |
(6.8) |
"l< "F "l< "F
Воставшейся части лекции вычислим энергию основного состояния ферми-газа
иэнергию Ферми для частиц без внутренней структуры (наиболее интересный пример газ свободных электронов).
6.3.Спектр импульса частиц в конечном объеме
Для частиц без внутренней структуры одночастичный гамильтониан имеет вид
^ |
|
|
p^i2 |
|
hi |
= |
|
|
(6.9) |
|
||||
|
|
2m |
Естественно взять в качестве одночастичных состояний jp;~ msi. Волновые функции таких состояний имеют вид
1 |
ei~p ~r ;ms |
(6.10) |
'p;m~s (~r; ) = pV |
Проблема с координатной волновой функцией волновые функции с различными p~ не ортогональны друг к другу.
Замечание. Волновые функции можно нормировать на -функцию; они будут ортогональны, но тогда нужно считать, что частица может двигаться во всем пространстве. Это нефизическое предположение, так как газ всегда занимает конечный объем.
Рассмотрим сначала одномерное движение
2

1 |
eipxx=~: |
(6.11) |
'p (x) = p |
xLx
Потребуем, чтобы все 'px (x) были периодическими функциями с периодом Lx, т. е. чтобы выполнялось условие
'p |
(0) = 'p |
(Lx): |
(6.12) |
|
x |
x |
|
Ясно, что теперь px не может быть любым действительным числом. Найдем его возможные значения. Подставляя в (??) явное выражение (??) для собственной функции, приходим к условию
eipxLx=~ = 1: |
(6.13) |
Это условие выполняется, если pxLx=~ кратно 2 . Отсюда следует, что
px = |
2 ~ |
nx; nx = 0; 1; 2; : : : |
(6.14) |
Lx |
Отметим, что для px получился дискретный спектр значений. Теперь, вычисляя скалярное произведение h 'px j'px0 i легко проверить, что
|
h 'px j'px0 i = px; px0 ; |
|
|
(6.15) |
|||
где введен символ Кронекера |
|
|
|
|
|
|
|
px; px0 = ( |
0; |
px |
= px0 |
: |
(6.16) |
||
|
|
1; |
px |
= px0 |
; |
|
|
|
|
|
|
6 |
|
|
|
Соотношение (??) показывает, что собственные функции (??) образуют ортонормированный набор, если спектр значений импульса определяется формулой (??).
Переходу к неограниченной области движения частицы соответствует предел Lx ! 1. При этом спектр значений импульса (??) становится непрерывным. В собственных функциях (??) невозможно выполнить предельный переход, так как множитель перед экспонентой стремится к нулю1. Напомним, однако, что сама по себе волновая функция не является наблюдаемой, поэтому нужно определить правило для предельного перехода в выражениях, которые встречаются при вычислении наблюдаемых величин. Обычно приходится иметь дело с суммами типа
px |
f(px) nx |
f |
|
2Lx~ nx |
; |
(6.17) |
X |
X |
|
|
|
|
|
где f(px) гладкая функция импульса. Покажем, что в пределе Lx ! 1 сумма преобразуется в интеграл по определенному правилу.
1Хотя при Lx ! 1 сама волновая функция всюду стремится к нулю, размер области движения стремится к бесконечности, поэтому при любом Lx остается справедливым соотношение (??).
3

Изобразим дискретные собственные значения импульса (??) точками на оси px. При большом Lx соседние точки расположены очень близко друг к другу, т. е. спектр импульса “почти непрерывный”. Возьмем интервал px, который значительно превышает расстояние 2 ~=Lx между соседними собственными значениями импульса. Число точек, попавших в интервал px равно pxLx=2 ~. Разделим теперь всю ось px на интервалы px, пронумеруем их индексом i и запишем приближенное значение суммы (??) в виде
Xpx |
|
L |
|
|
f(px) |
2 x~ Xi |
f(pxi) px; |
(6.18) |
где pxi произвольная точка в интервале с номером i. Например, это может быть середина интервала. Пусть Lx растет. Одновременно можно уменьшать px. При этом сумма в правой части (??) будет стремиться к интегралу от функции f(px) по всем px. Мы получаем правило перехода от дискретного спектра импульса к непрерывному:
Lx |
px |
f(px) ! |
1 |
f(px) dpx |
при Lx ! 1: |
(6.19) |
||
2 ~ Z |
||||||||
1 |
X |
|
1 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
Обобщим все сказанное выше на трехмерный случай2. Дискретный спектр импульса соответствует предположению, что частица находится в параллелепипеде со сторонами Lx, Ly Lz и объемом V = LxLyLz. Нормированные на единицу и ортогональные друг к другу собственные функции импульса имеют вид
1 |
ei~p ~r=~; |
(6.20) |
'p~ (~r ) = pV |
а дискретный спектр проекций импульса дается формулами
p |
|
= |
2 ~ |
n ; |
p |
|
= |
2 ~ |
n |
; |
p |
|
= |
2 ~ |
n |
; |
(6.21) |
x |
|
y |
|
z |
|
||||||||||||
|
|
|
x |
|
|
|
y |
|
|
|
|
z |
|
|
|||
|
|
|
Lx |
|
|
|
Ly |
|
|
|
|
Lz |
|
|
где nx, ny, nz независимо принимают значения 0; 1; 2; : : : Правило перехода к непрерывному спектру формулируется следующим образом. Если f(p~ ) f(px; py; pz) гладкая функция, то при V ! 1 имеем
|
V |
px; py; pz f(p~ ) ! |
(2 ~)3 |
1 |
dpx |
1 |
dpy |
1 |
dpz f(p~ ) : |
(6.22) |
|
Z |
Z |
Z |
|||||||
|
1 |
X |
1 |
|
|
|
|
|
|
|
|
|
|
1 |
1 |
|
1 |
|
|
Обычно для простоты в формулах (??) полагают Lx = Ly = Lz = L. Тогда объем области движения V = L3.
2Обоснование приводимых утверждений оставляем читателю в качестве упражнения.
4