Основы статистической механики: ансамбли и распределения

Мост от микромира к макромиру

Газ в комнате содержит ~10²⁵ молекул. Записать уравнения движения для каждой — практически невозможно. Но нам и не нужна такая детальная информация: мы хотим знать температуру, давление, теплоёмкость. Статистическая механика — это мост от микроскопических степеней свободы к макроскопическим термодинамическим свойствам через вероятностные методы.

Ключевая идея Больцмана: когда у нас нет информации о микросостоянии, мы считаем все микросостояния с одинаковой энергией равновероятными. Это гипотеза равных вероятностей — фундамент статистической механики.

Микроканонический ансамбль: изолированная система

Условия: Фиксированы E (энергия), V (объём), N (число частиц). Система изолирована.

Число микросостояний: Ω(E, V, N) — количество способов, которыми система может иметь данную энергию.

Энтропия Больцмана:

S = k_B ln Ω

Разбор: k_B = 1.38 × 10⁻²³ Дж/К — постоянная Больцмана, связывает микроскопическое (Ω) с макроскопическим (S). Логарифм: потому что энтропия аддитивна (S_AB = S_A + S_B для несвязанных систем), а Ω мультипликативно (Ω_AB = Ω_A × Ω_B).

Температура из энтропии: 1/T = ∂S/∂E|{V,N}. Давление: p/T = ∂S/∂V|{E,N}. Химический потенциал: μ/T = −∂S/∂N|_{E,V}. Это фундаментальное соотношение термодинамики — всё из одной функции S(E,V,N)!

Канонический ансамбль: система при постоянной температуре

Условия: Фиксированы T, V, N. Система в тепловом контакте с большим резервуаром (термостатом).

Вероятность микросостояния s:

P(s) = e^{−E(s)/k_BT} / Z = e^{−βE(s)} / Z, β = 1/(k_BT)

Статистическая сумма:

Z(T, V, N) = Σₛ e^{−βE(s)}

Физический смысл: Z — нормировочный множитель, но содержит всю термодинамическую информацию.

Связь с термодинамикой: Свободная энергия Гельмгольца F = −k_BT ln Z. Из F легко получить: ⟨E⟩ = −∂(ln Z)/∂β, давление p = −∂F/∂V, энтропия S = −∂F/∂T.

Флуктуации энергии: ⟨(ΔE)²⟩ = k_BT² C_V, где C_V — теплоёмкость. Относительные флуктуации: ΔE/⟨E⟩ ~ 1/√N → 0 при N → ∞. При N ~ 10²⁵ флуктуации энергии составляют 10⁻¹²·⁵ от среднего — практически ноль.

Числовой пример: двухуровневая система (спины)

Система N независимых спинов, каждый имеет два состояния: εᵢ ∈ {0, Δ} (Δ — щель).

Статсумма одного спина: z = 1 + e^{−βΔ}. Всей системы: Z = zᴺ.

Средняя энергия: ⟨E⟩ = −N ∂(ln z)/∂β = NΔ e^{−βΔ}/(1 + e^{−βΔ}) = NΔ/(1 + e^{βΔ})

При T → 0 (β → ∞): ⟨E⟩ → 0 (все в основном состоянии). При T → ∞ (β → 0): ⟨E⟩ → NΔ/2 (равновероятны оба состояния).

Теплоёмкость: C_V = k_B (βΔ)² e^{βΔ}/(1 + e^{βΔ})². Максимум при k_BT ≈ 0.42Δ — «аномалия Шоттки». Это характерно для двухуровневых систем — парамагнетики, туннельные дефекты в стёклах.

Большой канонический ансамбль: переменное число частиц

Условия: Фиксированы T, V, μ (химический потенциал). Система обменивается частицами с резервуаром.

P(s) = e^{−β(E(s) − μN(s))} / Ξ

Большая статсумма: Ξ = Σₛ e^{−β(E−μN)} = Σ_N Z(T,V,N) e^{βμN}

Квантовые газы: Для невзаимодействующих частиц: квантовая статсумма отличается от классической знаком «±» в знаменателе.

Фермионы (полуцелый спин): ⟨n_k⟩ = 1/(e^{β(ε_k−μ)} + 1) — распределение Ферми–Дирака. Каждый уровень заполняется максимум одной частицей (принцип Паули).

Бозоны (целый спин): ⟨n_k⟩ = 1/(e^{β(ε_k−μ)} − 1) — распределение Бозе–Эйнштейна. Бозоны «любят» компанию — они стремятся в уже заполненные состояния.

Реальные применения: Теплоёмкость металлов — Ферми-газ электронов. Чёрное тело (излучение) — Бозе-газ фотонов (μ=0). Сверхтекучий ⁴He — бозоны. Квантовый компьютер — двухуровневые системы (кубиты).

Бозе-Эйнштейновский конденсат (БЭК): При T < T_c все бозоны «падают» в основное квантовое состояние — система ведёт себя как единая макроскопическая квантовая волна. Предсказан в 1924 году, экспериментально реализован в 1995 году (Корнелл, Виман, Кеттерле; Нобелевская премия 2001). Сверхтекучесть жидкого ⁴He — частный случай БЭК. Конденсат применяется для создания атомных лазеров и прецизионных измерений констант тяготения.

Статистическая физика в материаловедении и технологиях

Статистическая механика связывает микроскопические свойства вещества с макроскопическими термодинамическими характеристиками. В полупроводниковой промышленности распределение Ферми–Дирака определяет число носителей в зонах проводимости и валентной: при легировании кремния фосфором уровень Ферми смещается, управляя электрическими свойствами. Дизайн новых батарей для электромобилей требует понимания свободной энергии Гельмгольца как функции концентрации ионов лития — задача статистической механики смешанных систем. В разработке сплавов метод Монте-Карло моделирует конфигурационную энтропию и тепловое расширение кристаллических решёток, предсказывая диаграммы фазового равновесия. В белковой биологии уравнение состояния полимера описывает конформационную энтропию полипептидной цепи: именно статистическая механика объясняет, почему белок сворачивается в уникальную нативную структуру — это минимум свободной энергии. В физике конденсированного состояния теория стекла, магнетизм и квантовые критические явления описываются через разделение функции разбиения — центрального объекта статистической механики. В атмосферной физике распределение Больцмана объясняет барометрическую формулу: давление атмосферы экспоненциально убывает с высотой, задавая профиль плотности атмосферы.

Метод Монте-Карло в статистической физике (алгоритм Метрополиса, 1953) позволяет вычислять термодинамические средние для сложных систем без аналитического вычисления функции разбиения: цепь Маркова с переходными вероятностями Больцмана сходится к каноническому ансамблю. Это применяется в молекулярном моделировании лекарств, разработке новых материалов и моделировании климата — везде, где система слишком сложна для точного аналитического подхода.

Задание: (а) Двухуровневая система: ε ∈ {0, Δ}. Вычислите Z, ⟨E⟩(T), C_V(T). Постройте C_V/k_B в зависимости от k_BT/Δ. При каком T максимум? (б) Идеальный газ: Z_N = (V/λ³)^N/N!, λ = h/√(2πmk_BT) — тепловая длина де Бройля. Найдите F, p, S, ⟨E⟩. Получите уравнение состояния pV = Nk_BT.