Главная страница
Навигация по странице:

  • 4.10. Термодинамика кристаллической решетки. Теория Дебая

  • Физика. Конспекты лекций и задачи казань 2015 удк 536. 7(07) ббк 22. 317 А62


    Скачать 1.13 Mb.
    НазваниеКонспекты лекций и задачи казань 2015 удк 536. 7(07) ббк 22. 317 А62
    АнкорФизика
    Дата24.05.2022
    Размер1.13 Mb.
    Формат файлаpdf
    Имя файлаThermodynamics_and_statistical_physics.pdf
    ТипКонспект
    #546850
    страница11 из 22
    1   ...   7   8   9   10   11   12   13   14   ...   22

    k
    ) жидкости; в соответствии с законами
    Рис. 4.3

    k
    N
    C
    V
    3/2


    0
    Рис. 4.2

    86 сохранения энергии и импульса
    M
    0
    V
    2
    /2 = M
    0
    V′
    2
    /2 +

    k
    ,
    M
    0
    V
    = M
    0
    V
    '+

    k
    , или

    (k
    V) =
    k
    +

    2
    k
    2
    /2M
    0
    . Второе слагаемое для макроскопических тел с большой массой M
    0
    ничтожно мало, поэтому минимальная величина скорости V, при которой возникает тормозящее возбуждение, достигается при V
    k и равна V =
    k
    /

    k
    . При спектре вида

    k
    =

    2
    k
    2
    /2M возможно возникновение тормозящих возбуждений при сколь угодно малой скорости V. Если же спектр фононного типа,

    k
    =

    vk
    (v
     скорость распространения возбуждений), то при скоростях V < v тормозящие возбуждения не возникают, и движение тела происходит без трения. На рис. 4.3 показана схема реального спектра возбуждений в HeII. Теоретически такая схема может быть получена при учете взаимодействий между атомами гелия. Минимальная скорость, начиная с которой возникает сопротивление движению, определяется, очевидно, значением k, соответствующим минимуму на кривой дисперсии
    (k).
    Пунктиром на рисунке отмечена линейная связь между
     и k, пунктиром с точками
     спектр свободных частиц.
    Отметим, что бозе-эйнштейновская конденсация при низких температурах и малых концентрациях наблюдалась недавно на металлах
    23
    Na (3s
    1
    ) и
    87
    Rb (5s
    1
    ). Объём конденсата в рубидии составлял доли кубического миллиметра, число атомов N

    10 3
    –10 7
    , температура конденсации была порядка 10
    −8
    К. Авторы этих работ были удостоены нобелевской премии 2001 года.
    4.9. Черное излучение
    Электромагнитное поле в полости
    (кубическом ящике) посредством разложения векторного потенциала в ряд Фурье можно представить в виде суперпозиции плоских волн:
    ,
    )
    (
    ],
    *
    [
    0
    t
    i
    i
    i
    e
    t
    e
    t
    e
    t
    t
    k
    k
    k
    r
    k
    k
    r
    k
    k
    k
    A
    a
    )
    (
    a
    )
    (
    a
    )
    A(r,







    (4.33) причем

    k
    = ck, (
    a
    k
    k) = 0 (условие поперечности волн). При выводе этих соотношений используются волновое уравнение, калибровка
     = 0, условие Лоренца div
    A
    = 0. Накладывая периодические граничные условия, вновь получаем возможные значения волнового вектора

    87
    k = (2
    /L)(n
    x
    , n
    y
    , n
    z
    ), n
    i
     целые числа. Энергия поля в полости
    (E = (
    1/c)A/t, B = rotA):
    H
    =







    *
    2
    )
    (
    8 1
    2 2
    2 2
    k
    k
    k
    k
    a
    a
    B
    E
    c
    V
    dV
    Переходя к вещественным переменным, характеризующим поле:
    *)
    (
    4
    /
    *),
    (
    4
    /
    s s
    2
    s s
    s s
    2
    s
    k
    k
    k
    k
    k
    k
    k
    k
    a
    a
    c
    V
    i
    Q
    P
    a
    a
    c
    V
    Q









    (s
     индекс поляризации, различает две независимые амплитуды волны в плоскости, перпендикулярной направлению распространения
    k), имеем
    H
    ).
    (
    s
    2
    s
    2 2
    s
    2 1




    k
    k
    k
    k
    Q
    P
    (4.34)
    Детали вычислений, приводящих к этому гамильтониану поля, см., например, в книге Ландау и Лифшиц, Теория поля, 1960. Таким образом, с каждой волной
    ks соотносится осциллятор с частотой

    k
    , причем Q
    k
    s
    , P
    k
    s
    оказываются каноническими координатами и импульсами для поля. Число осцилляторов поля бесконечно, в силу чего в классической теории полная энергия и теплоемкость поля бесконечны.
    Квантование поля осуществляется переходом к операторам поля, причем операторы координат и импульсов подчиняются каноническим перестановочным соотношениям
    ]
    ˆ
    ,
    ˆ
    [
    s s
    ss
    i
    Q
    P




    kk
    k
    k

    Возможные энергии поля складываются из энергий независимых осцилляторов: E(...n
    k
    s
    ...) =


    s s
    k
    k
    k
    n

    (считая энергию основного состояния равной нулю), а соответствующие стационарные состояния определяются набором квантовых чисел возбуждения n
    ks
    осцилляторов
    (...n
    k
    s
    ...)
    ...n
    k
    s

    . Квант возбуждения осциллятора удобно считать квазичастицей поля
     фотоном  с энергий


    k
    и импульсом
    k
     Тогда энергия и импульс поля складываются из энергий и импульсов всех фотонов, а их средние значения определяются средним числом фотонов. Переходы между соседними стационарными состояниями осцилляторов можно рассматривать как порождение или уничтожение фотонов; операторы рождения (b
    ks+
    ) и уничтожения (b
    ks
    ) фотонов естественно выражаются через координаты и импульсы осцилляторов:
    ).
    (
    2
    /
    ˆ
    ),
    (
    2
    /
    ˆ
    s s
    s s
    s s
    k
    k
    k
    k
    k
    k
    k
    k
    b
    b
    i
    P
    b
    b
    Q











    88
    При таком определении
    1
    ,
    1 1
    s s
    s s
    s s
    s s






    k
    k
    k
    k
    k
    k
    k
    k
    n
    n
    n
    b
    n
    n
    n
    b
    Средняя энергия осциллятора при температуре окружения

    (ср. (2.30)):
    ]
    1
    )
    /
    [exp(
    ,
    1
    s s
    s














    k
    k
    k
    k
    k


    n
    n
    (4.35)
    Как видно, среднее число фотонов дается распределением Бозе –
    Эйнштейна, причем химический потенциал
     = 0. Последнее обстоятельство можно связать с тем, что общее число фотонов N не фиксировано и подлежит определению из условия минимума свободной энергии: (
    F/N)
    ,V
    (= ) = 0.
    Переход от суммирования по
    k к интегрированию осуществляется так же, как переход от суммы по импульсам к интегралу для систем частиц:





    s
    3
    )
    ,
    (
    )
    2
    (
    )
    ,
    (
    k
    k
    k
    k
    s
    d
    s
    f
    V
    s
    f
    Для функций, зависящих только от частоты,








    s
    2 3
    2
    ,
    )
    (
    )
    (
    k
    d
    f
    c
    V
    f
    k
    то есть плотность распределения осцилляторов по частотам
    D
    (
    ) = V
    2
    /

    2
    c
    3
    . Поэтому энергия излучения, приходящаяся на интервал частот (
    ,  + d) (формула Планка, рис. 4.4):
    3 2 3
    ( )
    exp(
    / ) 1
    V
    d
    dE
    c
     
     

      


    (4.36)

    d
    dE

    1

    2
    >

    1

    m1

    m2

    Рис. 4.4

    89
    Распределение энергии по длинам волн выглядит следующим образом:
    2 5
    8
    ( )
    exp(2
    ) 1
    V c
    d
    dE
    c


     


     


    В предельных случаях низких и высоких частот получаются формулы
    Рэлея – Джинса:
    dE =
    (V

    2
    /
    2
    c
    3
    )d
     (

     << ) и Вина
    dE =
    (V

    /

    2
    c
    3
    )

    3
    exp(


    /)d (

     >> ).
    Максимум спектральной плотности энергии излучения dE/d

    достигается на частоте


    m
    = 2,822
    , линейно возрастающей с температурой. (При исследовании распределения энергии излучения по длинам волн максимум приходится на
    2 4.965 .
    m
    c
      


    ) Эта формула позволяет оценить поверхностную температуру звезд на основе спектра их излучения:
    T
    (K) =


    m
    /2,822k
    B
    Полная энергия излучения в полости:
    E =
    ,
    4
    )
    (
    15 4
    3 4
    2
    VT
    c
    c
    V
    SB





    где постоянная Стефана-Больцмана

    SB
    =

    2
    k
    B4
    /60

    3
    c
    2
    = 5,67.10
    5
    г.с
    3
    .K
    4
    Теплоемкость C
    V
    =
    4E/
    , энтропия  = C
    V
    /3, свободная энергия
    F =
    E/3. Давление излучения p = E/3V = (4
    SB
    /3c)T
    4
    .
    Оценки при
    Т
    = 300 К дают E
    rad
     5.10
    5
    (эрг/см
    3
    ) против E
     7.10 5
    для атомарного идеального газа с N/V = 10 19
    ; p
    rad
     2,5.10
    15
    T
    4
    против p
    gas
     1,4.10 3
    T
    , что объясняет возможность пренебрежения излучением при обычных условиях.
    Исследуем энергию излучения из полости через небольшое отверстие. Энергия, перенесенная за 1 сек через единицу поверхности полости в направлении
    nd
     (рис. 4.5), определяется фотонами с kn, находящимися в цилиндре длиной c с осью, параллельной
    n (объем цилиндра ccos
    )
    dI
    (
    )=(E/V)(c/4)cosd
    (в силу изотропности излучения в полости на направление d
     приxодится доля d
    /4 полной энергии излучения).

    90
    Полная интенсивность излучения на единицу поверх-ности I =

    SB
    T
    4
    (закон
    Стефана – Больцмана).
    Излучение вне полости уже не равновесно, однако спектральный состав его тот же, что и в полости.
    4.10. Термодинамика кристаллической решетки. Теория Дебая
    Твердое тело можно представить как совокупность N атомов, занимающих определенное (в среднем) положение в пространстве вследствие сильной связи с остальными атомами и совершающих малые колебания относительно этих равновесных положений. Из механики известно, что такую систему в гармоническом приближении
    (при учете лишь членов второго порядка в разложении энергии системы в ряд по смещениям атомов из равновесных положений) можно представить как набор 3N
    6 (
    3N) независимых осцилляторов.
    Кванты соответствующих упругих колебаний
     фононы  образуют систему невзаимодействующих бозонов.
    При таком подходе учитывается лишь движение атомов как целых
     решеточная часть всей системы. Электронные или спиновые степени свободы при необходимости учитываются отдельно. Статсумма решетки в гармоническом приближении:
    Z =






    l
    l
    1
    )]
    /
    exp(
    1
    [

    (произведение по всем осцилляторам решетки); свободная энергия

    F = F
    0
    +






    l
    l
    )],
    /
    exp(
    1
    ln[

    (4.37) энергия E = E
    0
    +





    l
    l
    l
    ,
    1
    )
    /
    exp(

    где F
    0
    = E
    0
    = N

    0
     потенциальная энергия системы атомов в равновесном положении,

    0
     энергия связи, отнесенная к одному атому. К F
    0
    относится и энергия нулевых колебаний.
    Рис. 4.5

    91
    Представление о характере решеточ- ных колебаний можно составить, рассмат- ривая либо колебания кристалла, облада- ющего идеальной периодической структу- рой, либо колебания сплошной однородной упругой среды. В изотропной упругой среде объемом
    L
    3
    распространяются продольные и поперечные звуковые волны с волновыми векторами
    k = (2
    /L)(n
    x
    , n
    y
    , n
    z
    ) и частотами, соответственно,

    l
    =v
    l
    k
    ,

    t
    =v
    t
    k
    ,(4.38) где v
    l
    , v
    t
     скорости продольного и поперечного звука, определяемые упругими константами и плотностью среды; для каждого
    k имеются две независимые поперечные волны. С каждой волной можно связать осциллятор, как и в случае электромагнитного излучения, но, поскольку теперь общее число осцилляторов ограничено, величины
    k должны быть ограничены сверху (k
    0

    /a, a  размеры кристаллической ячейки). Более реалистичная атомная модель кристалла приводит и к значительно более сложной картине колебаний. В кристаллах со сложной структурой элементарной ячейки (в ячейке содержится более чем один атом) помимо акустических ветвей колебаний, приближенно напоминающих звуковые волны сплошной среды, имеются оптические колебания с

    opt
    (k
    =
    0)
     0. Схема кривых дисперсии волн в кристалле приведена на рис. 4.6.
    Модель Дебая колебаний твердого тела допускает, что все колебания сводятся к акустическим с законом дисперсии (4.38); имеется предельная частота колебаний

    D
    (дебаевская частота), определяемая из условия, что общее число осцилляторов равно 3N.
    Плотность распределения акустических колебаний с данной поляризацией s (= l, t) по частотам определяется как и в случае световых волн, D
    s
    (
    ) = V
    2
    /2

    2
    vs
    3
    . Поэтому общее число колебаний в интервале от 0 до

    D
    :
    3N =
     


















    s
    D
    D
    t
    l
    s
    D
    v
    V
    v
    v
    V
    d
    D
    0 3
    2 3
    3 3
    1 3
    3 2
    2 2
    1 2
    )
    (
    Отсюда


    D
    v
    l
    v
    t
    k
    k
    0
    Рис. 4.6

    92 3
    /
    1 2
    )
    /
    6
    (
    V
    N
    v
    D



    , (4.39) а температура Дебая
     =


    D
    (/k
    B
    ). Энергия твердого тела:
    1 2
    3 0
    /
    3 3
    2 0










    D
    e
    d
    v
    V
    E
    E


    (4.40)
    Свободная энергия в дебаевской модели получается из (4.37) интегрированием по частям:
    F = F
    0
    +


























    /
    0 3
    /
    3 3
    }.
    1
    )
    1
    ln(
    {
    3
    x
    e
    dx
    x
    e
    N
    (4.41)
    При низких температурах,
     << , энергия, свободная энергия, энтропия и теплоемкость равны:
    E
    = N

    0 3
    4 5
    3











    N
    , F = F
    0
    3 4
    5











    N
    ,
     =
    3 4
    5 4









    N
    ,
    C =
    3 4
    5 12









    N
    (4.42)
    (для твердых тел C
    p
    C
    V
    , см. задачу 4.40). При высоких температурах,
     >> ,
    E = N

    0
    + 3N
    , F = F
    0
    3Nln(e
    1/3
    /),  = 3Nln(e
    4/3
    /), C = 3N (4.43)
    (закон Дюлонга и Пти).
    Свободная энергия (4.41) при любой температуре имеет вид
    F = F
    0
    +
    f(/), так что теплоемкость C = 
    2
    F
    /
    
    2
    = (/)f'' (штрихи означают производные по
    /). Таким образом, теплоемкости дебаевских тел, находящихся в «соответственных состояниях» (при одинаковых приведенных температурах
    /), совпадают (закон соответственных состояний).
    Уравнение состояния Дебая:
    p
    =
    F
    0
    /
    V

    (
    /V)[f(/)  f
    '
    (
    /)(/)] = E
    0
    /
    V + E
    D
    /V,(4.44) где E
    D
    =
    f  f
    '
     дебаевская часть энергии решетки,  = (V/)(/V) =
    =
    dln/dlnV  постоянная Грюнайзена (можно пренебречь ее зависимостью от температуры).
    Формула Грюнайзена связи коэффициента линейного расширения изотропных тел с теплоемкостью:
     =
    3 3
    1 1
    3 1
    V
    C
    p
    p
    V
    V
    l
    l
    V
    V
    V
    p









    

    


    














    




    93
    Модель Дебая довольно грубая, но она хорошо описывает термодинамические свойства твердого тела при очень низких температурах, когда основной вклад в эти свойства вносят низкочастотные акустические колебания, а другие колебания
    «заморожены». При высоких же температурах все осцилляторы вносят одинаковый классический вклад в энергию и теплоемкость, и детали моделей не играют роли. В качестве усовершенствования модели оптические колебания сложных кристаллов иногда рассматривают отдельно, полагая, что все колебания одной моды обладают одинаковой частотой (комбинированная модель Дебая и Эйнштейна).
    1   ...   7   8   9   10   11   12   13   14   ...   22


    написать администратору сайта