Розподіл Ґіббса

Матеріал з testwiki
Перейти до навігації Перейти до пошуку

Розподіл Ґіббса — розподіл, що визначає кількості частинок в різних квантових станах. Ґрунтується на таких постулатах статистики:

  1. Всі доступні мікростани системи рівноймовірні.
  2. Рівновазі відповідає найімовірніший розподіл (підсистем за станами).
  3. Ймовірність перебування підсистеми в деякому стані визначається лише енергією стану.

Розподіл Ґіббса являє собою найзагальнішу і зручну основу для побудови рівноважної статистичної механіки.

Кількісний розгляд

Статистична сума

G=N!N1!N2!,(0)

як і в термодинаміці, має зміст відносної ймовірність знаходження системи в певному микростанів. І, дивлячись на співвідношення Больцмана S=klnG, легко зрозуміти, що станам з максимальною ентропією відповідає максимальна статистична вага. Потрібно врахувати, що в системі постійні число частинок

iNi=N=const(1)

і повна енергія

iNiεi=E=const.(2)

Факторіал великих чисел (а числа N і Ni великі; тими з них, які малі, можна знехтувати) знаходиться за формулою Стірлінга: N!=2πN(Ne)Nexp(ϑ12N), де 0<ϑ<1. Цю точну формулу можна замінити наближеною

N!=2πN(Ne)N,(3)

так як відносна помилка в обчисленнях за цією формулою не перевершує e112N1112N, вже при n=10 вона менше одного відсотока. Із співвідношень (0), (1) і (3) випливає наступне:

G=N!i2πNiNiNieNi=N!ieNi(i2π)(iNiNiNi)=N!eiNi(2π)0,5NiNiNiNi=N!eiNi(2π)0,5NiNiNi+0,5.

Чисельник тут є функція від N, і можна ввести позначення

C(N)=N!eiNi(2π)0,5N,

що дає

G=C(N)iNiNi+0,5.(4)

Тоді з формули Больцмана S=klnG слідує

S=ki((Ni+0,5)lnNi)+const.

Тут можна знехтувати 0,5 порівняно з Ni. Тоді

S=ki(NilnNi)+const.(5)

Максимум ентропії (5) із урахуванням співвідношень (1) і (2), використовуючи метод невизначених множників, буде при умовах

lnNidNi=0,dNi=0,εidNi=0.

Звідси (lnNi+β+αεi)dNi=0, де α и β — множники Лагранжа, не залежні від змінних Ni. У системі є m змінні і три рівняння - отже, будь-які дві залежать від інших; відповідно можна вважати залежними N1 та N2 і вибрати множники Лагранжа так, щоб коефіцієнти при dN1 и dN2 звернулися в 0. Тоді при інших dNi змінні N3, N4, … можна прийняти за незалежні, і при них коефіцієнти також будуть рівні 0. Так отримано

lnNi+β+αεi=0,

звідси

N¯i=N0eαεi,

де N0=eβ — нова константа.

Для визначення сталої α можна скласти систему в теплопровідні стінки і квазістатично змінювати її температуру. Зміна енергії газ а одно dE=εidN¯i, а зміна ентропії (зі співвідношення (5)) дорівнює dS=klnN¯idN¯i=kαεidN¯i. Так як dE=TdS, то звідси α=1kT, і тому

N¯i=N0eεikT.(6)

Термостат

Отримано найбільш ймовірне розподіл системи. Для довільної макроскопічної системи (системи в термостаті), оточеній протяжної середовищем (термостатом), температура якої підтримується постійною, виконується співвідношення (6) - розподіл Гіббса: їм визначається відносна ймовірність того, що система при термодинамічній рівновазі знаходиться в i-вому квантовому стані.

Див. також

Джерела

  • Базаров И. П., Геворкян Э. В., Николаев П. Н. Термодинамика и статистическая физика. Теория равновесных систем. — М.: МГУ, 1986. — 312 с.
  • Квасников И. А. Термодинамика и статистическая физика. Теория равновесных систем. Статистическая физика. — Том 2. — М.: УРСС, 2002. — 430 с.
  • Кубо Р. Статистическая механика. — М.: Мир, 1967. — 452 c.
  • Сивухин Д. В. Общий курс физики. — В 5 т. — Т. II. Термодинамика и молекулярная физика. — М.: ФИЗМАТЛИТ, 2005.
  • Терлецкий Я. П. Статистическая физика. — 2-е изд. — М.: Высшая школа, 1973. — 277 c.