Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

лекции по статам

.pdf
Скачиваний:
201
Добавлен:
10.05.2015
Размер:
5.3 Mб
Скачать

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

Наиболее важные операторы динамических переменных для многочастичной системы имеют вид

 

 

 

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