МИФИ_Вычметоды КФ
.pdf
силами, т.е. в поле упругих волн – звуковых колебаний, а затем перейти в непрерывный предел.
Рассмотрим одномерную периодическую систему с периодом a из одинаковых атомов массы M, соединенных пружинами с одинаковой жесткостью k (рис. 4.1).
Напишем уравнение движения для такой системы. Для этого пронумеруем все атомы и рассмотрим уравнение Ньютона для n-го атома. На него действуют упругие силы со стороны (n − 1) -го атома,
направленные влево, и со стороны (n + 1) -го атома, направленные вправо.
Рис. 4.1. Периодическая система атомов массы M, соединенных одинаковыми пружинами с жесткостью k
Если обозначить смещение атома из положения равновесия как Un ,
то
MU&&n = −k(Un − Un−1) + k(Un+1 − Un ) = k(Un+1 + Un−1 − 2Un). (4.26)
Перейдем теперь к непрерывному пределу, устремляя период a к нулю, и разложим смещение в ряд до второго порядка по a:
Un+1 |
= U(xn + a) = U(xn) + U′(xn )a + |
1 U′′(xn)a2 |
; |
|
|
2 |
(4.27) |
|
|
|
|
Un−1 |
= U(xn − a) = U(xn) − U′(xn)a + |
1 U′′(xn)a2. |
|
|
|
2 |
|
Тогда уравнение движения приобретает привычный вид волнового уравнения:
61
&& |
2 |
∂2U |
|
|
MU(x) = ka |
|
|
. |
(4.28) |
|
∂x2 |
|||
|
|
|
|
Введем линейную плотность ρ = Ma и модуль упругости G = ka , тогда
&& |
∂2U |
|
|
|
ρU(x) = G |
|
. |
(4.29) |
|
∂x2 |
||||
|
|
|
Квадрат скорости звука при этом определяется как c2 = Gρ . Имеем,
таким образом, одномерное волновое уравнение, описывающее собственные колебания струны с линейной плотностью ρ и
модулем упругости G.
Напишем гамильтониан системы – сумму кинетической и потенциальной энергий атомов:
MU& 2 |
|
|
(U |
|
− U |
|
)2 |
|
|
|
|
U& 2 |
|
|
(∂U |
n |
/∂x)2 |
|
|
|||||||
H = ∑ |
n |
+ k |
|
|
n |
|
|
|
n−1 |
|
|
|
= a∑ ρ |
n + G |
|
|
|
= |
||||||||
2 |
|
|
|
|
2 |
|
|
|
|
2 |
||||||||||||||||
n |
|
|
|
|
|
|
|
|
|
|
|
n |
2 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
& |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
(4.30) |
||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
(∂U/∂x) |
|
|
|
|
|
|||||||
|
|
|
∫ |
|
|
|
(U(x)) |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
= |
dx |
ρ |
|
|
|
|
|
|
|
|
+ G |
|
|
|
. |
|
|
|
|
|||||
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Введем фурье-компоненты смещений по координате: |
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
U(x,t) = |
1 |
|
|
∑U(k,t)eikx , |
|
|
|
|
|
|
(4.31) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
тогда гамильтониан можно переписать следующим образом:
|
|
& |
& |
* |
(k,t) |
|
* |
(k,t) |
|
|
U(k,t)U |
|
U(k,t)U |
(4.32) |
|||||
H = a∑ ρ |
|
|
|
|
+ Gk2 |
|
. |
||
|
2 |
|
|
2 |
|||||
k |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
В последнем выражении учтено, что для того, чтобы смещение было действительной величиной, необходимо потребовать, чтобы
U* (−k,t) = U(k,t) .
62
Пользуясь уравнением на собственные колебания, можно получить зависимость от времени фурье-компонент смещений. Действительно,
&& |
2 |
U(k,t) U(k,t) = U1 (k)e |
iωkt |
+ U2 |
(k)e |
−iωkt |
; |
||||||||||||||||
U(k,t) = −c |
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
ωk = c |
|
k |
|
. |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Далее воспользуемся условием U* (k,t) = U(−k,t) . Имеем: |
|
||||||||||||||||||||||
|
* |
|
|
|
* |
(k)e |
−iωkt |
|
+ U |
* |
(k)e |
iωkt |
; |
|
|
|
|
||||||
U (k,t) = U |
|
|
|
|
|
2 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
U* (−k,t) = U* (−k)e−iωkt |
+ U* |
(−k)eiωkt |
U* (−k) = U |
2 |
(k). |
||||||||||||||||||
|
1 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
||
Окончательно находим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
U(k,t) = U(k)e−iωkt + U* (−k)eiωkt ; |
|
|
|
|
|||||||||||||||||||
U(x,t) = |
1 |
∑(U(k)eikx−iωkt |
+ U* (k)e−ikx+iωkt ). |
|
|
|
|||||||||||||||||
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
L |
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Введем обобщенные координаты и импульсы: Q(k,t) = U(k,t);
P(k,t) = ρQ& (k,t).
Перепишем оператор энергии в обобщенных координатах:
P P |
|
ρω2 |
|
|
||
H = ∑ |
k |
−k |
|
k |
|
|
2ρ |
+ |
2 |
QkQ−k , |
|||
k |
|
|
|
|||
|
|
|
|
|
||
(4.33)
(4.34)
(4.35)
(4.36)
(4.37)
здесь Pk ,Q−k – сопряженные обобщенные координаты и импульсы.
Дифференцируя оператор энергии по импульсу, координате, а также используя уравнения движения, можно получить уравнения Гамильтона
|
& |
∂H |
|
|
|
Q−k = |
|
; |
|
||
∂Pk |
|
||||
|
|
|
(4.38) |
||
|
|
∂H |
|||
|
& |
|
|||
Pk = − |
|
|
. |
||
∂Q−k |
|||||
|
|
|
|||
63
Для выражения гамильтониана через операторы в представлении чисел заполнения введем операторы ak (t), ak+ (t) следующим образом:
|
|
|
|
|
Pk (t) |
|
|
|
|
ak (t) = |
ρωk Qk (t) + i |
|
|
; |
|||||
|
|
|
|
||||||
|
2h |
2ρωkh |
|||||||
|
|
|
|
|
(4.39) |
||||
|
ρωk |
|
|
|
|||||
ak+ (t) = |
Q−k (t) − i |
|
P−k (t) |
|
. |
||||
|
|
|
|
|
|||||
|
|
||||||||
|
2h |
|
2ρωkh |
||||||
Введенные таким образом операторы рождения и уничтожения обладают всеми атрибутами операторов, фигурирующих ранее в задаче для одного осциллятора. В этом несложно убедиться, вводя правила коммутации для обобщенных координат и импульсов: [Pk ,Q−k ] = −ih .
Задача 4.3. Проверить коммутационные соотношения для операторов ak (t), ak+ (t) .
Окончательно гамильтониан системы представляется как сумма энергий гармонических операторов:
|
|
1 |
|
(4.40) |
H = ∑hωk ak+ ak + |
2 |
, |
||
k |
|
|
|
|
при этом через эти операторы можно выразить и смещение:
U(k,t) = Q(k,t) = |
|
|
h |
|
(ak (t) + a−+k (t)) ; |
|
||
|
|
|
|
|||||
|
|
|
|
|
2ρωk |
(4.41) |
||
|
|
|
|
|
|
|
|
|
U(x,t) = ∑ |
h |
|
|
(ak (t)eikx + ak+ (t)e−ikx ) . |
|
|||
2ρLω |
|
|||||||
k |
|
|
|
|
||||
|
k |
|
|
|
|
|||
Зная решение задачи для осциллятора, далее можно рассчитать и термодинамические свойства данной системы.
Теперь рассмотрим формализм вторичного квантования (или представление чисел заполнения) в общем случае.
64
4.3.Формализм вторичного квантования
Волновые функции в шредингеровском представлении неудобны для описания ансамбля тождественных частиц. Действительно, если рассмотреть функцию двух тождественных частиц Ψ(a,b) ,
то функция Ψ(b,a) описывает ту же амплитуду вероятности того,
что одна частица находится в точке a, а вторая – в точке b. Следовательно, значения этих функций не могут задаваться произвольно, и должны различаться только фазовым множителем, так как модули квадратов этих волновых функций описывают одну и ту же физическую величину и должны быть равны:
Ψ(a,b) 2 = Ψ(b,a) 2 . Эту проблему устраняют, потребовав, чтобы
волновая функция двух тождественных частиц была либо симметрична, либо антисимметрична относительно перестановки аргументов. Для системы, состоящей из многих частиц, такая проблема еще более усугубляется, так как следует потребовать симметрию или антисимметрию любой пары частиц. В случае антисимметрии результат представим в виде так называемого детерминанта Слэтера [3], работать с которым чрезвычайно неудобно. Иногда предлагается также альтернативный метод – схема Юнга [1] – для учета тождественности, но он также неудобен для практического применения.
Предлагаемый ниже аппарат вторичного квантования (или формализм чисел заполнения) позволяет избежать вышеуказанных трудностей, удобен для работы (в том числе, для формализации задач на компьютере), и специально направлен на описание систем тождественных (неразличимых) частиц. Будет показано, как естественным образом все операторы физических величин можно представить через операторы рождения и уничтожения фермиевского или бозевского вида и как практически работать с ними.
65
4.3.1.Одночастичный базис
Введем вакуумную волновую функцию 0
, обозначающую
состояние рассматриваемой системы, не содержащее ни одной частицы. Затем добавим в систему одну частицу, и рассмотрим полный набор таких одночастичных состояний. В качестве одночастичных состояний можно выбрать, например, плоские
волны, образующие полный набор: k
= eikx . В формализме чисел
заполнения такие состояния будут представлены следующим образом:
k ≡ ak+ |
0 . |
(4.42) |
Это тождество означает, что введенный оператор рождения ak+ по
определению рождает частицу в одночастичном состоянии, описываемом плоской волной с волновым вектором k.
Итак, одночастичное состояние определяется как результат действия оператора рождения на вакуумное состояние. Волновая функция, сопряженная с одночастичной, представляется в виде
0 ak , где ak – сопряженный с ak+ оператор (принято говорить, что
кет-вектору |
k |
соответствует бра-вектор |
k |
). Из-за |
|
|
|
|
|
ортонормированности одночастичного базиса справедливо следующее очевидное соотношение:
0 |
akak+′ |
0 = δkk' . |
(4.43) |
|
|
|
|
Для случая одной частицы пока не получено никаких новых результатов, кроме переобозначений волновых функций. Формализм вторичного квантования становится полезным при рассмотрении систем двух и более частиц.
Подчеркнем, что полный набор – необязательно плоские волны, а квантовое число k – необязательно волновой вектор. Более того, функции любого другого базиса всегда можно записать в виде линейных комбинаций плоских волн. Значит, и операторы
66
рождения, соответствующие новому базису, также можно выразить через линейную комбинацию операторов ak+ .
Пусть имеется новый набор базисных функций n
. Тогда
n = ∑ |
|
k k |
|
n . |
(4.44) |
|
|
||||
k |
|
|
|
|
|
Это соотношение следует из полноты базиса и тождественного соотношения 1 ≡ ∑ k
k . Соответствующий набор операторов
k
рождения дается следующими линейными комбинациями операторов ak+ :
an+ = ∑ak+ k |
|
n . |
(4.45) |
|
|||
k |
|
|
|
Оператор an+ рождает частицу в состоянии n
.
В частном случае, если ввести координатные волновые функции, являющиеся фурье-образом от плоских волн
|
x = ∑e−ikx |
|
k , |
(4.46) |
|
|
|||
|
k |
|
|
|
то соответствующие операторы рождения |
|
|||
a+ (x) ≡ Ψ + (x) = ∑e−ikxak+ |
(4.47) |
|||
|
k |
|
||
рождают частицу в пространственной точке x (их называют полевыми операторами – операторами поля частиц).
4.3.2. Двухчастичный и многочастичный базис. Коммутационные соотношения
Рассмотрим теперь полный набор состояний двух частиц. В обычном (шредингеровском) представлении можно взять, например, произведение плоских волн
k′,k′′ = eik′xeik′′x , |
(4.48) |
но такой выбор неоднозначен – заданным значениям векторов k′,k′′ соответствуют два разных состояния: в одном первая
67
частица имеет импульс k′ , а вторая – k′′ , в другом – наоборот. Положим, что частицы неразличимы, т.е. осуществляется только одно физическое состояние, а именно состояние, в котором одна частица обладает волновым вектором k′ , а другая – волновым вектором k′′ . Получается, что в шредингеровском представлении такое вырождение по симметрии мы должны учесть, хотя никакой дополнительной физической информации этот факт не вносит в задачу.
Посмотрим, как записать это состояние в формализме вторичного квантования:
k′,k′′ → a+ a+ |
0 . |
(4.49) |
k′ k′′ |
|
|
В результате действия двух операторов рождения на вакуумное состояние рождаются две частицы, одна частица с волновым вектором k′ , а другая – c волновым вектором k′′ . Очень важно, что номера частиц при этом не вводятся, однако следует помнить, что одно и то же вторично квантованное двухчастичное состояние можно записать двумя способами. Так, для заданных векторов k′ , k′′ можно записать еще одно состояние
k′′,k′ → a+ |
a+ |
0 . |
(4.50) |
k′′ |
k′ |
|
|
Это то же самое физическое состояние, поэтому отличаться эти функции могут только множителем, т.е.
a+ |
a+ |
= Fa+ a+ |
, |
(4.51) |
||
k′′ |
k′ |
|
k′ |
k′′ |
|
|
a+ a+ |
= |
1 a+ |
a+ . |
|
||
k′ |
k′′ |
|
F k′′ |
k′ |
|
|
Последние соотношения (коммутационные) должны выполняться в любом базисе и должны быть инвариантными относительно переходов от одного базиса к другому. Ясно, что множитель F не должен зависеть от k′ , k′′ и их взаимного расположения в коммутационном соотношении. Единственные числа, которые удовлетворяют этим условиям, это F = ±1 . Таким образом, возможны только два вида коммутационных соотношений:
ak+′′ak+′ = +ak+′ak+′′ , |
(4.52) |
или
68
ak+′′ak+′ = −ak+′ak+′′ . |
(4.53) |
Итак, возможна либо коммутация, либо антикоммутация операторов рождения. Частицы, операторы которых коммутируют, называются бозонами, частицы с антикоммутационными соотношениями – фермионами. Фермионами являются электроны, протоны и нейтроны, слагающие атомы, а бозонами являются, например, фотоны – кванты света, фононы – кванты звуковых колебаний, нейтральные атомы гелия и др. Все элементарные частицы
разделены на эти два основных класса.
Проанализируем важнейшие свойства квантовых систем, обладающих бозеили ферми-симметрией. Пока речь не идет о термодинамике, обсуждаются лишь квантовые состояния систем.
Посмотрим, что произойдет, если волновые вектора совпадают: k′ = k′′ . Для случая бозонов получаем тривиальное тождество, а для фермионов находим
ak+′ak+′ = −ak+′ak+′ ≡ 0 . |
(4.54) |
Получаем важнейшее свойство фермионов – в одном и том же квантовом состоянии два фермиона находиться не могут – знаменитый принцип (или запрет) Паули.
Еще одно свойство – связь между коммутацией (перестановкой) операторов рождения и перестановкой частиц. Рассмотрим движение двух частиц в системе их центра масс. Тогда соотношения коммутации примут следующий вид
a+ |
a+ |
= ±a+ |
a+ |
+q |
, |
(4.55) |
K |
+q K−q |
K |
−q K |
|
|
причем 2K – импульс центра масс, 2q – импульс относительного движения.
Пусть центр масс движется равномерно и описывается плоской волной, а нас интересует функция относительного движения. Тогда
волновая |
функция |
двух |
частиц |
представима |
в |
виде |
|
r |
r |
r |
r |
rr |
r |
|
|
Ψ(r1 |
,r2) = Ψ(r,K) = Φ(K)ϕ(r) = e2iKr ∑g(q)ϕKr |
(q) , где r,K |
– радиусы- |
||||
|
|
|
|
q |
|
|
|
векторы |
относительного движения и |
движения центра |
масс; |
||||
69
функция Φ(K) отвечает за движение системы как целого, а ϕ(r) – за внутренние степени свободы. Последнюю функцию разложим в
ряд |
Фурье, |
|
компоненты |
|
которого |
представимы |
|
|
в |
числах |
|||
заполнения: ϕr |
|
r |
a+ |
|
0 . Учтем уравнения коммутации: |
||||||||
(q) = a+ |
|
||||||||||||
|
K |
|
K |
+q K−q |
|
|
|
|
|
|
|
|
|
∑g(q)aK+ +qaK+ −q |
|
0 = ±∑g(q)aK+ −qaK+ +q |
|
0 = ±∑g(−q)aK+ +qaK+ −q |
|
0 . |
(4.56) |
||||||
|
|
|
|||||||||||
q |
|
|
q |
|
|
|
|
q |
|
|
|
|
|
Здесь |
g(q) – |
|
коэффициенты Фурье. |
Последнее |
соотношение |
||||||||
получается заменой q → −q под знаком суммы. Видно, что функции g(q) должны быть четными для бозонов и нечетными для
фермионов (для простоты мы пренебрегли здесь учетом других степеней свободы, например спинов частиц, что, однако, не изменит общности результатов).
Четность функций g(q) напрямую связана со значением
орбитального момента количества движения, так что орбитальный момент количества движения должен быть четным для бозонов и нечетным для фермионов. Но процедура изменения знака импульса относительного движения эквивалентна процедуре перестановки между собой частиц. Более того, функции g(q) – это
шредингеровские волновые функции относительного движения, только в импульсном представлении. Поэтому отсюда следует, что бозонная волновая функция должна быть симметричной относительно перестановки частиц, а фермионная –
антисимметричной.
Эрмитово сопряженная двухчастичная волновая функция (бравектор) вводится следующим очевидным образом:
k′,k′′ |
|
→ 0 |
|
ak′′ak′. |
(4.57) |
|
|
Нормировка этой волновой функции
0 |
a |
k′′ |
a |
k′ |
a+ a+ |
0 = N |
k′k′′ |
(4.58) |
|
|
|
k′ k′′ |
|
|
|||
|
|
|
|
|
|
|
|
|
должна быть выбрана из соображений удобства. Более общее соотношение ортонормированности имеет следующий вид:
0 |
ak′′′ak′ak+′ak+′′ |
0 = Nk′k′′δk′′′k′′ . |
(4.59) |
70
