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

МИФИ_Вычметоды КФ

.pdf
Скачиваний:
0
Добавлен:
13.02.2026
Размер:
9.98 Mб
Скачать

Следовательно, таким градиентным преобразованием волновой функции удается убрать векторный потенциал из уравнения Шредингера, но при этом появляется фазовый множитель в граничных условиях [18].

Из (6.42) и (6.45) – (6.46) следует, что для ввода в систему поля и тока достаточно на границе системы ввести фазовый множитель

ei χ = ei

e

Φ

2πiΦ

 

2πhc

 

(6.47)

hc

= e Φ0 ; Φ0

=

;

 

e

 

 

 

 

 

 

 

здесь Φ0 – квант потока. Удобно магнитный поток отсчитывать в

единицах этой величины. Таким образом, при вводе магнитного потока методом изменения граничных условий (6.47) вся система в тороидальной геометрии пронизывается магнитным полем. Более того, так как узлы системы эквивалентны, поток может быть введен в любом месте и даже распределен равномерно или неравномерно вдоль замкнутого контура по оси x. При этом результаты расчета не изменятся, если равны циркуляция векторного потенциала и введенный в систему суммарный поток.

В частном случае системы, описываемой единой волновой функцией Ψ , квантово-механическое выражение для плотности тока, переносимого частицами системы массы m и заряда e , имеет следующий вид:

r

 

 

e

(

 

*

)

 

eh

(

*

 

 

 

* )

 

e2

 

 

 

 

2

r

(6.48)

 

 

 

 

 

 

 

 

j =

2m

Ψ PΨ + h.c.

=

2mi

Ψ

Ψ − Ψ

 

Ψ

 

mc

 

Ψ

 

 

A ,

 

 

h

 

e r

 

 

 

 

 

 

 

 

 

 

 

 

 

здесь P =

 

 

 

A

– оператор

обобщенного

импульса. Первая

i

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

часть выражения (6.48) для тока – парамагнитный вклад, приводящий к усилению магнитного поля в среде, вторая – диамагнитный, приводящий к экранировке внешнего поля. Ток, циркулирующий в системе, зависит от волновой функции и векторного потенциала, при этом если градиенты волновой функции малы (мал парамагнитный вклад), то ток и векторный потенциал просто пропорциональны друг другу:

r

 

e2

 

 

2

r

eh

 

 

2

 

 

j

 

 

Ψ

 

A = −

 

χ

Ψ

 

.

(6.49)

mc

 

m

 

 

 

 

 

 

141

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последнее соотношение справедливо, например, в сверхпроводниках, и приводит к объяснению эффекта выталкивания индукции магнитного поля из объема массивного сверхпроводника – эффект Мейсснера – Оксенфельда [19, 20].

Таким образом, для описания системы во внешнем магнитном поле и в присутствии тока следует внести магнитный поток Φ либо на границе, либо распределить его вдоль заданного контура.

Рассмотрим теперь модель Бозе – Хаббарда (6.14) во внешнем магнитном поле. Введем распределенный поток так, чтобы при движении частицы вдоль оси x при каждом перемещении с узла на узел добавлялась определенная часть циркуляции векторного потенциала. Кинетическая часть гамильтониана при градиентном преобразовании волновой функции преобразуется следующим образом:

 

 

 

 

e r 2

 

 

 

 

 

 

 

 

 

 

p −

c

A

 

 

 

 

 

 

 

 

 

*

 

 

 

3

 

 

 

 

 

 

tij (A) = Ψi

 

 

 

 

Ψj d r

 

 

 

 

(6.50)

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i′)*

p2

Ψj′ d3r = d3i*

p

 

e

−iχ(r )+iχ(r

)

= tije

i(χ

−χ )

.

 

 

 

Ψj

i j

 

j

i

2m

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Появляется фазовый множитель с разностью фаз конечного и начального состояния в показателе экспоненты [21]. В общем случае, при перескоке частицы с узла i на узел j следует записать следующее слагаемое:

ai+ aj ai+aj exp

2πi j r

r

(6.51)

A dl .

Φ0 i

 

 

Очевидно, этот фазовый множитель появляется из-за градиентной перенормировки волновой функции. Если движение частицы перпендикулярно контуру обхода (т.е. перпендикулярно оси x), то интеграл в экспоненте равен нулю, т.е. при любых движениях частицы по системе в (6.51) учитывается только проекция этого движения на ось x. Прыжки частицы туда и обратно вдоль оси x будут давать вклады c противоположными по знаку фазами, так что гамильтониан будет эрмитовым, и его кинетическая часть примет следующий вид:

142

+

 

2πi j r

r

(6.52)

 

 

Φ

 

Adl

 

H = −tai aj exp

0

.

 

<ij>

 

 

i

 

 

 

Если положить, что фаза монотонно наращивается при движении

r

 

частицы вдоль оси x, т.е. векторный потенциал A = Φ

ex

постоянен

Lx

 

 

и, следовательно, интеграл в фазовом множителе в (6.52) одинаков для ближайших соседей, то для случая кубической симметрии или тороидальной геометрии имеем:

 

 

 

 

2πiΦ

 

 

2πiΦ

 

 

 

 

 

H = −tai+ai+er

x

eΦ0Lx

+ ai++er

aie

Φ0Lx

− t(ai+ai+er

+ ai++er

y

ai) −

 

 

 

 

x

 

 

 

 

y

 

 

i

 

 

 

 

 

 

 

 

 

i

 

 

(6.53)

 

 

 

 

 

 

 

U ni(ni − 1) + Vninj.

 

− t(ai+ ai+erz + ai++erz ai) +

 

 

 

i

 

 

 

 

 

2

i

<ij>

 

 

 

После формирования гамильтоновой матрицы получаем задачу на собственные значения эрмитовой матрицы. Недиагональные матричные элементы, возникшие из-за движения вдоль направления x, будут иметь фазовый множитель, зависящий от введенной в систему градиентной фазы. Получившиеся в результате диагонализации собственные значения также будут зависеть от Φ .

Можно задать другое распределение векторного потенциала, положив его нулю везде, кроме границы (Lx ,y,z) ↔ (1,y,z) . При

прохождении границы полагаем, что

 

 

 

 

 

 

2πiΦ

 

a+

 

a

a+

 

a e Φ0 ;

 

Lx ,y,z

1,y,z

Lx ,y,z

1,y,z

 

(6.54)

 

 

 

 

 

 

 

2πiΦ

a+

 

 

a+

 

 

 

a

Lx ,y,z

a

 

e

Φ0 .

 

1,y,z

 

1,y,z

Lx ,y,z

 

 

 

Таким образом, всю фазу можно учитывать только на границе, перпендикулярной оси x, и результаты расчета от этого не

изменятся. В частном случае, если Φ = Φ20 , т.е. в систему введена половина кванта потока, что эквивалентно разности фаз χ = π ,

143

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

a+

a

−a+

a

;

Lx ,y,z 1,y,z

Lx ,y,z 1,y,z

(6.55)

a+

 

 

−a+

 

 

a

Lx ,y,z

a

Lx ,y,z

.

1,y,z

 

1,y,z

 

 

Оператор тока в узельном представлении имеет вид [21]:

 

 

 

 

 

∂H

ite

 

 

2πiΦ

 

 

2πiΦ

 

 

 

 

 

 

 

 

Φ

L

 

 

 

Φ

L

 

 

 

 

 

jx = −

 

=

ai+ai+er x e

 

0

 

x

 

− ai++er x aie

0

 

x

,

(6.56)

 

 

∂Ax

 

 

 

 

 

 

 

 

 

 

 

 

 

hc i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для малых Φ получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ite

 

(ai+ai+er

 

 

2πiΦ

(ai+ai+er

 

 

 

 

 

 

 

jx

 

 

− ai++er

ai) +

+ai++er

ai) .

(6.57)

 

 

 

 

 

Φ→0

hc

 

i

 

x

x

 

Φ L

x

 

 

 

x

 

 

 

x

 

 

 

 

 

 

 

 

 

 

0

 

 

 

i

 

 

 

 

 

 

 

В отсутствие фазы ток равен нулю (первое слагаемое в (6.57) при усреднении дает нуль, так как направления вдоль и против оси x эквивалентны), а при наличии магнитного поля напрямую зависит от введенного магнитного потока.

Рассмотрим теперь зависимость энергии системы свободных частиц на решетке от введенного магнитного потока. Возьмем для расчета вариант с равномерно распределенным векторным потенциалом:

 

 

 

 

2πiΦ

 

 

2πiΦ

 

 

H = −tai+ai+er

x

eΦ0Lx

+ ai++er

x

aie

Φ0Lx

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− t(ai+ai+er y + ai++er y ai) − t(ai+ ai+erz

 

+ ai++erz ai); .

i

 

 

 

i

 

 

 

 

 

i x,y,z

Применим фурье-преобразование в следующем виде:

aj =

 

1

akei(kxxj +kyyj +kzzj ).

 

 

 

 

LxLyLz kxkykz

(6.58)

(6.59)

144

Подставив (6.59) в выражение для гамильтониана (6.58), получим:

 

H = −t

1

 

 

 

 

ak+ak (Tx+ + Tx+ Ty+ + Ty+ Tz+ + Tz);

 

 

 

 

 

 

 

 

 

 

 

LxLyLz kx ky kz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k′x k′y k′z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xi yi zi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2πiΦ

 

 

 

 

 

 

 

 

 

 

 

i

k

x

(x

+e

x

)−k′

x

i

+(k

y

−k

)y

i

+(k

−k′ )z

i

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T+

= e

 

 

 

 

 

 

i

 

 

 

 

 

 

x

 

 

 

 

 

 

 

y

 

 

 

 

z

z

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ0Lx

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2πiΦ

 

 

 

 

 

 

 

 

 

 

 

i(k

 

x

 

−k′

 

(x

+e

x

)+(k

y

−k′

)y

+(k

−k

′ )z

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

x

 

i

 

 

x

 

 

i

 

 

 

 

 

 

 

 

 

 

y

 

 

i

 

 

z

 

z

 

 

 

 

 

 

 

 

;

 

 

 

 

 

= e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ0Lx

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T+

= ei((kx −k′x )xi +ky (yi +ey )−k′yyi +(kz −k′z )zi);

 

 

 

 

 

(6.60)

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

= ei((kx −k′x )xi +kyyi −k′y (yi +ey )+(kz −k′z )zi);

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T+

= ei((kx −k′x )xi +(ky −k′y )yi +kz (zi +ez )−k′zzi );

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

= ei((kx −k′x )xi +(ky −k′y )yi +kzzi −k′z (zi +ez )).

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δk

 

 

k′

 

 

=

1

 

 

 

ei(kx −k′x )x ,

 

 

 

 

 

 

 

 

 

 

 

 

(6.61)

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

x

 

 

 

 

 

Lx

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i(k

e

 

+

2πiΦ

)

 

 

 

−i(k

 

e

 

+

2πiΦ

)

 

 

 

H = −t ak+

k

y

k

z

ak

k

y

k

 

e

 

 

x

 

x

 

 

Φ0Lx

 

 

+ e

 

 

 

x

 

x

 

 

Φ0Lx

 

 

kxkykz

x

 

 

 

 

 

 

 

x

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ikyey

+ e

−ikyey

 

 

 

 

(6.62)

 

 

− t akxkykz akxkykz

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kxkykz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

ak

 

 

 

 

 

 

 

 

 

ikzez

 

+ e

−ikzez

 

 

 

 

 

 

 

 

 

 

− t ak

x

k

y

k

z

k

k

 

 

e

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

kxkykz

 

 

 

 

 

 

 

 

 

x

 

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда спектр системы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2πΦ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. (6.63)

ε(k,Φ) = −2t cos(kxex

 

 

 

Φ0Lx

 

 

) + cos(kyey ) + cos(kzez )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

ex , ey , ez

 

 

 

 

 

 

 

 

 

 

модули

 

 

 

 

 

векторов

 

 

 

 

трансляции

по

соответствующим осям, равные периоду решетки в этих направлениях. Таким образом, получаем обычный одночастичный спектр в приближении сильной связи, только по направлению kx

появляется вклад от фазы [21]. Фактически все импульсы частиц

145

системы получили фазовый сдвиг в направлении приложенной фазы или внешнего тока (рис. 6.5).

Рис. 6.5. Сдвиг спектра системы при учете градиентной фазы вдоль оси x, k = 2πΦ / Φ0Lx

Перебирая одночастичные состояния, находим спектр всей системы.

Заметим, что все фазовые зависимости спектра, рассмотренные выше, справедливы и для ферми-систем [21].

Полная энергия системы, рассчитанная с учетом фазы, удовлетворяет условию периодичности по фазе

E0 (Φ + Φ0 ) = E0 (Φ) .

(6.64)

Это свойство можно получить непосредственно из спектральной зависимости одночастичных состояний. Рассмотрим конкретный пример – одномерную систему из Na = Lx = 6 узлов и N = 3 частиц.

В отсутствие фазы разрешенные импульсы в этой системе будут (рис. 6.6)

k0

= 0; k±1

= ±

π

; k±2

= ±

; k3 = π .

(6.65)

3

3

 

 

 

 

 

 

 

В случае конечного значения Φ все импульсы получают сдвиг, равный πΦ /3Φ0 . Так как в основном состоянии все частицы имеют

нулевой импульс, то

 

 

π Φ

 

 

E0

 

 

 

 

 

(6.66)

 

 

 

= −6t cos

3 Φ0

.

 

 

 

 

 

146

 

 

 

 

 

Рис. 6.6. Одночастичные состояния для свободных частиц на цепочке из шести узлов

Выражение (6.66) справедливо в интервале 0 ≤

 

Φ

 

 

1

до тех пор,

 

 

2

 

 

Φ0

 

 

 

 

 

пока сдвиг импульсов не будет равен k =

0.5

 

k−1

− k0

 

=

π

. В

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

случае дальнейшего возрастания фазы минимальная разрешенная одночастичная энергия будет отвечать не импульсу k0 , а импульсу

k−1 (рис. 6.7), все частицы перейдут на этот уровень, и энергия основного состояния станет равной

 

 

π

 

Φ

 

(6.67)

E0

 

 

(1 −

 

 

= −6t cos

3

Φ0

) .

 

 

 

 

 

 

Далее при достижении значения Φ = Φ0 ситуация повторится. Получаем, таким образом, периодическую зависимость E0 (Φ) с периодом Φ0 (рис. 6.8).

Ток в системе также будет периодичен:

 

 

 

 

2πt

 

 

 

πΦ

 

 

 

 

 

Φ0

 

 

 

 

 

 

 

 

 

 

≤ Φ ≤

;

 

∂E

Φ0

sin

0

; 0

2

J = −

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.68)

∂Φ

2πt

π

 

Φ

 

 

 

Φ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ0

sin

3

(

Φ0

1)

;

2

≤ Φ ≤ Φ0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

147

 

 

 

 

 

 

 

 

Рис. 6.7. Сдвиг энергетических уровней при градиентном преобразовании

Рис. 6.8. Зависимость энергии основного состояния от фазы (6.66) – (6.67)

Подобные зависимости имеют место и для других физических величин в различных невзаимодействующих системах, в том числе и с ферми-статистикой.

В термодинамическом пределе, при Na → ∞, Lx → ∞ , для

невзаимодействующего бозе-газа можно получить следующую зависимость энергии основного состояния от фазы:

при 0 ≤ Φ 1 заполнены только состояния с нулевым импульсом,

Φ0 2

поэтому

148

 

 

 

 

 

E0 (Φ) = εk

 

 

 

 

 

 

 

 

 

 

 

2πΦ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nk

k=k0

 

 

 

 

 

 

 

 

 

 

Lx

→∞ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −2tNcos

 

Φ0Lx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.69)

 

 

 

 

 

 

 

 

 

πΦ

 

 

 

 

 

 

(

 

 

 

 

 

)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ Φ0

 

 

 

 

 

 

 

 

 

 

 

= −2tN 1 − 2

 

 

 

 

 

 

 

= E0 (0) + 4tπ

N

 

 

 

 

;

 

 

 

 

 

 

 

Φ0L

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при

1

 

Φ

≤ 1

заполнены только состояния с импульсом

k = k−1 ,

2

 

 

 

 

 

Φ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E0 (Φ) = εk

 

 

 

 

 

 

 

 

 

 

2πΦ

 

 

 

 

 

 

 

 

 

 

 

 

 

nk

 

 

 

 

 

 

 

 

 

 

 

Lx →∞ =

 

 

 

 

 

 

k=k−1

= −2tNcos

Φ

0Lx

 

Lx

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.70)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

(

 

 

 

 

 

)2

 

 

 

 

 

 

 

 

 

π

 

 

 

Φ

 

 

 

 

 

 

2

 

 

 

1 − Φ Φ0

 

 

 

 

 

 

 

= −2tN 1 − 2

 

 

 

 

 

= E0 (0) + 4tπ N

 

 

 

 

 

2

 

 

 

.

 

 

 

 

L

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ0

 

 

 

 

 

 

 

 

 

 

 

Lx

 

 

 

 

 

Последнее слагаемое в (6.69) и (6.70) – вклад фазы или токовое слагаемое, оно зависит от размерности системы. В трехмерной ситуации оно расходится пропорционально линейному размеру системы L с ростом L, так как при NNa → const имеем, что

Na = LxLyLz ~ L3x ; в двумерной ситуации Na = LxLy ~ L2x оно выходит

на постоянное значение, а в одномерном случае спадает пропорционально 1 /L , так как Na = Lx .

Периодичность энергии и любых других характеристик системы

существует даже при наличии взаимодействия между частицами. Покажем это.

Из (6.46) и (6.47) следует, что при сдвиге разности фаз χ на 2π

волновая функция не меняется, т.е. любые решения уравнения Шредингера и любые квантово-механические средние должны быть периодическими функциями Φ с периодом Φ0 .

Далее энергия и термодинамические функции системы являются

четными

функциями Φ . Это следует из (6.46), так как замена

Φ → −Φ

эквивалентна изменению направления оси x. Таким

 

149

образом, вид любой термодинамической функции, в частности логарифма статистической суммы

Z = e−En / T

(6.71)

n

 

как функции Φ , качественно будет такой, как показано на рис. 6.9 (заметим, что зависимость, показанная на рис. 6.9, согласуется с точным решением (6.67) в пределе T → 0 , когда F = −TlnZ → E ).

Рис. 6.9. Качественная зависимость логарифма статистической суммы от фазы

Более того, можно показать, что в идеальной, без диссипации энергии, макроскопической системе реализуются только те значения магнитного потока, которые кратны кванту Φ0 , т.е. имеет

место квантование магнитного потока [18]:

 

Φ = nΦ0 , n− целое .

(6.72)

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

J = −

∂F

=

T ∂ lnZ

.

(6.73)

∂Φ

 

 

 

∂Φ

 

Отсюда следует, что в термодинамическом пределе в равновесии реализуются только экстремальные значения Φ , которые соответствуют условию

∂F

(6.74)

 

 

∂Φ = 0 .

 

150