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

книги из ГПНТБ / Беккер Р. Теория теплоты

.pdf
Скачиваний:
15
Добавлен:
25.10.2023
Размер:
15.46 Mб
Скачать

ность в точке сот . Интеграл по спектральному распреде­ лению дает общее число /V собственных колебаний.

Рисунок 86 показывает несколько типичных форм ко­ лебаний цепи, соответствующих определенным значени­ ям k. При граничном колебании (k — п/а; eihan=(— l ) n ) соседние узлы решетки имеют одинаковую и противопо­ ложно направленную амплитуду. Если в линейной решет­ ке выделить две части с шагом в каждой из них, то обе

Рис. 86. Форма колебаний цепи при частных зна­ чениях к.

эти части решетки колеблются с одинаковой амплитудой в противоположных направлениях. При таких колебаниях центры отдельных пружин неподвижны. Каждый атом движется таким образом, как если бы обе соседние пру­ жины были закреплены в их центрах. Но коэффициент жесткости для половины длины пружины равен 2\, поэто­ му с о 2 = 4 / / т , так как при колебании каждый раз нужно учитывать две пружины. Другой типичный случай возни­

кает при k

=

n/2a. Если записать

колебание при таком

значении k

в действительной форме:

 

qn

= cos kan cos a>t =

cos ~ - cos (at,

то можно заметить, что в этом случае колеблется только одна часть решетки (п = 0, 4-2, + 4 , ...), в то время как другая остается неподвижной. Следовательно, каждая отдельная колеблющаяся материальная точка движется

340

так, как если бы она колебалась между двумя соседни­ ми неподвижными точками. Поэтому (a2 — 2flm.

При очень небольших значениях k(ka<^l) ампли­ туды соседних атомов почти равны. В этом случае сме­ щения могут описываться с помощью медленно изменя­

ющейся непрерывной функции от

координат

 

 

q(X)

= e i ( k X - a t ) ,

 

(63.12)

где функция

q(na)=qn

определяет смещение

атома

номер п. Это

предельный случай

упругих волн. С по­

мощью уравнения колебаний (63.12) описываются упру­ гие волны в цепи, причем a/k—c, очевидно, представля­ ет собой скорость звука, т. е. скорость, с которой распро­ страняются волны. В упругой области связь между час­

тотой и «вектором распространения» k

звуковых волн

согласно уравнению

(63.8)

определяется из выражения

 

 

®=®J-~-

 

(63.13)

Частота пропорциональна к , скорость звука в соот­

ветствии с этим

равна:

 

 

 

£

= ш

« = а т / - £ - = с .

(63.14)

| k

I

2

j

т

 

На основании этого можно легко уяснить себе, как, зная

скорость звука, рассчитать

начальный

участок спектра,

не

рассматривая

подробнее

структуру

решетки

цепи.

Так

как известно,

что собственные колебания

имеют

форму звуковых волн по уравнению (63.12) и, кроме того,

известна

связь

между со и k, то можно снова выдвинуть

условие

периодичности

на длине L=Na

и опять полу­

чить уравнение

(63.7)

для возможных

значений k1. Мо­

гут использоваться также другие уравнения, а в уравне­

ние (63.11) вместо dco/dk можно

подставить

скорость

звука. Для упругого спектра это дает 2 :

 

г у п п (&) ^ю ~ ~ ^ ю =

~ ~

(63.15)

пс

ясо„.

 

1 Однако значения k не ограничиваются интервалом 2я/а. Это ограничение проявляется только вследствие того, что с целью уста­ новления числа собственных колебаний равным числу атомов берется лишь N самых низких значений.

2 То, что упругое спектральное распределение здесь строго по­ стоянно, обусловливается особенностями линейной схемы.

341

Если сравнить это выражение с действительным спектром решетки по уравнению (63.11), то можно за­ метить, что 2у П р точно соответствует постоянному на­ чальному участку спектра решетки

(со < co,n; г =

—).

Следовательно, теория упругости заменяет действи­ тельную функцию со(&) прямой ы — ck. Предельная час­ тота спектра соответствует k = n/a. Отсюда видно, что

I1 \ \

/

|

— ,

X i

/1*

-я/ а

о

л / а

 

 

 

 

 

 

 

 

 

Рис. 87.

Спектр

решетки

 

 

и упругий спектр цепи.

 

 

Упругое

 

спектральное

рас­

 

 

пределение

совпадает

со

 

 

спектром

 

 

решетки

 

при

ма­

 

 

лых частотах. Если

упругий

 

 

спектр

должен

 

содержать

 

 

такое же

 

число

частот,

как

 

 

и спектр

 

решетки,

то

он

 

 

должен

 

быть

оборван

на

 

 

частоте

со

 

 

 

 

 

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

ших длин волн. Учет структуры решетки в общем

случае

смещает

рассчитанные

по теории

упругости

частоты

в область меньших значений частот. Предельная

частота

по Дебаю и наибольшая частота

решетки

различаются

на множитель я/2, так

что упругий

спектр

простирается

в область

значительно

более высоких частот (рис. 87).

Особенно поучительно рассмотрение линейной цепи с двумя ато­

мами различной массы М и

т(М^т)

(рис. 88). Пусть положения

атомов при ненапряженных пружинах снова заданы с помощью расстояний па. В точках 2па находятся атомы с массой М, в точках ( 2 « + 1 ) а — атомы с массой т. Таким образом, рассматриваем такую

342

же схему, как и выше, только теперь атомы, чередуясь, имеют массы Мит. Уравнения движения.

 

 

= — / {2<72« — 2 ? 2 „ + i — q2n-i};

,

(63.16)

 

 

—~

f {2q2n+i

Чт —

Qin+i)

 

 

 

можно снова решить с помощью

оправдавших

себя

условий

 

 

qtn

= А е Ш а п - ш ;

 

)

 

 

 

 

 

 

 

 

}

 

(63 17)

 

 

? M + i = B e ' * ( 2 " + I , < H " " B ' - 1

 

 

- 5 а

-to.

0

а

 

За

 

-9

О

 

 

О

^

 

 

Т

 

 

ЯР

TI

 

ЧЗ

%/,

Рис. 88. Линейная цепь с двумя различными массами.

При этом амплитуды для различных видов атомов (Л, В) при­ нимаются различными. Диапазон значений К в данном случае огра­ ничен интервалом

ял

— < k < — , 2а 2а

так как решетка имеет период 2а. Если в качестве граничного ус­ ловия примем периодичность решения с периодом 2aN, то возмож­ ными значениями k будут:

я

I N

\

- 1 , 0 , 1

IN

\

N

k =

/; / = — — — 1

— - 1 ) , — .

Na

\ 2

/

'

\.2

 

) 2 .

 

 

 

 

 

 

(63.18;

Число

этих значений

равно N,

число атомов

каждого

вида в ин­

тервале периодичности 2aN также равно N. После подстановки урав­ нений (63.17) в уравнения (63.18) получаем систему уравнений для определения со2

 

МаРА =

2]А — 2/ cos kaB,

\

 

 

mco2 £ =

2/3 — 2f cos kaA,

j\

(63.19)

детерминант которой должен быть равен нулю.

 

 

(УИсо2 — 2/) (тсо2 — 2/) — 4/2 cos2

ka = 0.

(63.20)

Это квадратное

по со2 уравнение для каждого из значений k да­

ет две собственные

частоты:

 

 

 

 

со" = •Mm

\M + m — VM2

+ m* + 2mMcos2ka\ ;

(63.21

со2, = — -

IЛ1.+ т + Ум2

+ ш 2 + 2тМ cos 2ka\ .

(63.22)

Mm

 

 

 

 

 

Всего получаем 2N собственных частот, т. е. ровно столько же частот, сколько атомов в интервале периодичности 2aN.

.343

Из выражения (63.19) получаем для отношения амплитуд:

А_

cos ka

со2

— 2//от

М_

 

В

со2 — 2 / / М '

 

(63.23)

cos ka

 

 

 

 

 

 

т

 

Связь между со и k распадается на две ветви (рис. 89). К «аку­ стической» ветви (со) относятся меньшие частоты, здесь соседние ато­ мы колеблются в одинаковых фазах ( Л / В > 0 ) ; малые частоты соот-

Рис. 89. Спектр и форма колебаний линейной цепи из двух разновид­ ностей атомов с отношением масс М : от= 4 : 1.

В оптической ветви ш2 >2f/m, тогда согласно (63.23) Л/В<0. В акустической вет-

т

2 <2 //т, тогда согласно (63.23) Л/В>0. Форма колебаний (а — Ь) для трех частных значений k объяснена в тексте.

ветствуют упругим (акустическим) колебаниям цепи. В «оптической» ветви 1 (со+) соседние атомы колеблются с противоположно направ­ ленными амплитудами ( Л / В < 0 ) ; ее частоты лежат выше, чем для акустической ветви.

Особенно характерными

здесь

являются случаи (рис. 89):

 

2

М +

от А

от

 

1. k =

0, со

I = 2/

, — =— — . Обе части

решетки

 

f

Mm

В

М

 

колеблются

относительно друг друга как целые. Это точно

соответст­

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

только

амплитуды

соседних

атомов

имеют различное значение,

 

л

2

2/

 

 

2.

& = ——, со I

= —

, 4 = 0.

Часть решетки с более тяжелы-

 

^

от

 

 

ми атомами не движется, соседние атомы части решетки от колеблют­

ся с

противоположно направленной амплитудой. Частота со2 =2//от

1

Колебания, относящиеся к оптической ветви, особенно в окрест^

постях k = 0, лежат в основе оптических явлений.

344

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

3. k= — , со2

=

— ,

 

 

 

 

 

~

 

М

 

 

 

В—О. При таких

колебаниях

 

 

 

неподвижна часть

решетки

 

1

 

т, в то время как соседние

М= 1

 

атомы

с

массой

М колеб­

 

j Sу

Спектр

лются в противофазе.

 

Решетка

Дебая

На

рис. 90, а приведено

 

 

 

 

1

изменение

со (k)

для

одина­

 

 

ковых масс:

 

 

 

М-

 

|

и + (*) =

(4//М)

cos2 te/2;

 

а)

«'-уж '"JD

to 2 _ W =

(4//М)

sin 2 ta/2 .

 

 

 

 

 

 

 

 

 

 

(63.24)

 

 

 

Это изменение идентич­ но тому, которое было полу­ чено для линейной цепи с одинаковыми атомами, здесь только благодаря различию обеих частей решетки ис­ кусственно изменено деле­ ние на интервалы к. Значе­ ния k по уравнению (63.17) лежат в два раза плотнее, чем по уравнению (63.7), так как в этом случае мы рассматривали 2N атомов с массой М. На рис. 90, б при­ ведена зависимость при

Рис.

 

90.

Спектр

решетки

 

 

 

и спектр

Дебая.

 

 

 

 

 

 

а — при

одинаковых

массах

 

 

 

(М=т);

б — при М : m=9 : 1;

 

 

 

в

М : т=100 : 1.

При

боль­

MJm=t0O:is

 

 

шом

отношении

масс

спектр

 

 

 

Дебая

не

дает

хорошего

Спектр

 

 

приближения.

Частоты

оп­

Дебая

 

 

тической

ветви

практически

 

 

могут

заменяться

одной

час­

.

 

си

тотой

и „ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\

(',)

шс

-\Щ

 

 

 

 

 

 

 

 

т т

 

 

 

 

 

 

 

 

 

 

 

345

очень большой разнице масс (М^>т). Тогда корень в уравнениях (63.21), (66.22) можно разложить в ряд ':

 

«1 » (2//Ж) sin2 ка\

ш _ = V

2f/M

\ sin ka

|;

 

 

 

 

 

со2 . =

(2//УИ) ( l +

 

cos2

to};

ш +

=

 

 

(63.25)

 

Обе ветви проходят далеко друг от друга.

 

 

 

 

 

 

Спектральное

распределение (рис.

90,6)

получают

так

же, как

и в случае цепи с однородными атомами. При

сильно

отличающих­

ся

массах

также

получают спектр

из

двух

лежащих

далеко

друг

от

друга

функций распределения

(рис. 90, s).

Акустическая

ветвь

дает такой же спектр, как и в

случае равных масс.

Напротив,

оп­

тическая ветвь почти монохроматична. Частоты оптической

ветви

расположены симметрично по отношению

 

 

 

 

 

 

а

диапазон

частот

имеет относительную

ширину

т/2М

в соответст­

вии с уравнением

(63.25). Следовательно, в этом случае оптическую

часть спектра можно описать с помощью выражения

(62.5) для

мо­

дели Эйнштейна с

частотой сов. Соотношения

для пространственной

решетки похожи.

Д л я решеток,

компоненты

которых

имеют

очень

различные массы, по методу Дебая можно приближенно описывать лишь акустическую ветвь; для описания оптической ветви нужно привлекать метод Эйнштейна (например, для йодистого лития с от­ ношением масс 127:7). Описание полного спектрального распреде­ ления с помощью только приближенного метода Дебая в этом слу­ чая очень плохое.

64. КОЛЕБАНИЯ ПРОСТРАНСТВЕННОЙ РЕШЕТКИ

При рассмотрении пространственной решетки ограни­ чимся простейшей кубической решеткой с постоянной решетки а (рис. 91,а—в). Положение атомов этой ре­ шетки в равновесном состоянии определяется выраже­ нием

m

am,

(64.1)

 

где m = ( m b т2, Ы3) —вектор, имеющий лишь целочис­ ленные компоненты. Для обсуждения колебаний пред-

1 Акустическая ветвь имеет ту же форму, что и для цепи с од­

нородными атомами, в которой массы М связаны

пружинами

дли­

ной 2а. При таких колебаниях малые массы не

играют

никакой

роли. Таким образом, можно вернуться к формуле

(63.4а),

заменив

в ней а на и f на f/2,

так как вдвое более длинная пружина

имеет

половинный коэффициент

жесткости.

 

 

 

346

ставим, что ближайшие

соседние

атомы

этой

решетки

(6 соседних атомов,

расположенных

на

расстоянии а)

связаны пружинами /, в то время как

12

следующих со­

седних атомов (находящихся на расстоянии

|^2а) сое­

динены пружинами с жесткостью fСмещение

 

атома

m из состояния покоя

обозначим

вектором

q m

с тремя

 

 

 

 

 

 

 

m

компонентами этого

вектора по координатным

осям qi

( t = l , 2, 3, или i=x,

у, z). При таких

обозначениях вы-

Рис. 91. Элементарная ячейка.

а — положение атомов в простейшей кубической решетке а, б — атомы решетки связаны пружинами / с ближайшими соседями и пружинами /' со следующими по удалению соседями; в к определению смещений q.

ражение для потенциальной

энергии,

представленное

в виде квадратичного

ряда, имеет вид:

 

 

 

 

 

m , n

 

 

 

а уравнения движения

i I

 

 

 

 

 

 

 

 

Mqf

=—

ф™ qy .

 

(64.3)

 

 

 

n,l

 

 

 

Коэффициенты

ряда

потенциальной

энергии

имеют

очень наглядный

смысл,

который лучше

всего

выявля­

ется из уравнений движения.

O f f представляет собой

силу, действующую на атом

номер m в направлении i,

если атом номер п сдвинут в направлении / на отрезок s, в то время как все другие атомы остаются в состоянии

1 Пружинная связь лишь с ближайшими соседями не дала бы стабильной решетки.

347

покоя. На этом основании эти коэффициенты легко опре­

делить

исходя

из

принятых

пружинных

соединений

(рис. 92). Если рассматривать два атома m и п и сдви­

нуть атом номер п в направлении

/ на

расстояние s,

то

сила, действующая на атом m при небольшом

смещении

s, имеет направление связи обоих атомов. Величина этой

силы равна проекции смещения на соединительную ли­

нию,

умноженной

на

соответствующий

 

коэффициент

 

 

 

 

 

 

 

 

 

 

Рис. 92. К расчету сил, дей­

 

 

 

 

 

 

 

 

 

 

 

ствующих на атом п в на­

 

 

 

 

 

 

 

 

 

 

 

правлении Rn , если атом т

 

 

 

 

 

 

 

 

 

 

 

сдвинут на s в направлении

 

 

 

 

 

 

 

 

 

 

 

R m .

Если

смещение

s

имеет

 

 

 

 

 

 

 

 

 

 

 

компоненту

в

направлении

 

 

 

 

 

 

 

 

 

 

 

соединительной

оси,

возни­

 

 

 

 

 

 

 

 

 

 

 

кает

сила.

 

 

 

 

 

 

жесткости

пружины

f | m

- n |

, который зависит

только

от

значения

| т — п | , т.

е. от

расстояния

между

 

атомами:

 

 

 

 

ф т п =

_

 

Я т - п |

 

( n i — п ) г ( т — п ) /

 

 

 

, g 4 4

 

 

 

 

 

"

 

 

 

 

1 т — п | | т — п |

 

 

 

 

'

при

тфп.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Э Т О

определение

справедливо только при m=^=n. Ко­

эффициенты

OftN

определяют,

используя

силу,

действу­

ющую на атом т , когда только этот атом выведен из со­

стояния

 

покоя,

а

остальная

решетка

неподвижна.

Коэффициент

Ф ™ Т

непосредственно

связан

с

 

коэффи­

циентами, определяемыми из выражения (64.4). Ибо

если все атомы

получат

одно

и то же

смещение

q" =

s,

независимо от

п,

то

в

уравнение

движения

не

могут

войти никакие силы. Следовательно, должно выполнять­

ся условие

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ц ф т ' "

=

0

 

 

 

 

 

 

(64.5)

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

или же

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ф

Г

=

-

£

ФТ.

 

 

 

 

 

 

(64.5а)

 

 

 

 

 

 

 

 

п(фт)

 

 

 

 

 

 

 

 

 

 

Коэффициенты

Ф™"

зависят

только

от т — п . Поэ­

тому в уравнении (64.5а) безразлично, производится ли

348

суммирование по всем n ( ^ m ) или по всем п—т(Ф<д):

 

 

 

 

Ф5Г=-Е Ф П = Ф П .

 

 

 

(64.56)

 

 

 

 

 

 

п ( ^ 0 )

 

 

 

 

 

 

 

 

Естественно, что

Ф ^ т

не

зависит

от т . При пру­

жинном

соединении

согласно уравнению

( 6 4 . 4 )

 

 

 

 

 

ф п т . = 2

 

 

 

 

 

 

 

(64.5в)

 

 

 

 

 

 

n(=*0)

 

 

 

 

 

 

 

 

Собственные частоты снова

могут

определяться по

уравнению волны:

 

 

 

 

 

 

 

 

 

 

 

 

q m

= AS mm-*t).

 

q?

=

A . ei

« * « " - « , < ) •

 

( 6 4 . 6 )

 

Амплитуда

А = ( Л Ь

A2,

Аъ)

и вектор

распростране­

ния

k =

(ku

k2, ko) для случая

пространственной

решет­

ки

являются

векторами.

Если

ввести

эти

уравнения

в уравнение

движения,

то

получим линейную

систему

уравнений:

 

 

 

 

 

 

 

 

 

 

 

 

 

М О У М ,

= £ ф ™ У к

( R n

-

R m ) Al

= £

Ф°" g ' k

R 4 , ( 6 4 . 7 )

 

 

 

 

п , /

 

 

 

 

 

п/

 

 

 

 

 

Л>М, = £

(к) Л ;

/„ (к) =

£

 

 

 

(64.7а)

с симметричной трехмерной

матрицей

к.

которая сама

зависит

от

вектора

распространения

Для

 

каждого

вектора

к получаем

три частоты ^(k)

 

и три

вектора

А~( ", которые определяют форму колебания. Так как

значения к,

которые в каждой из компонент

отличаются

в

раз,

согласно

уравнению ( 6 4 . 6 ) дают

совпадаю­

щее

колебание, то

ограничения,

оправдавшие

себя

для значений к в линейной решетке,

справедливы

для

всех трех компонент к:

 

 

 

л

 

-—а< к < а

( 6 4 . 8 )

 

_ _ £ L < A . < J L >

а" а

В соответствии с этим векторы к ограничены кубиком с ребром 2л/а, построенным исходя из нулевой точки.

349

Соседние файлы в папке книги из ГПНТБ