Что такое функция Бриллюэна?
Функция Бриллюэна \(B_J(x)\) описывает намагниченность парамагнетика, состоящего из атомов с квантовым числом полного момента импульса J. В статистической механике безразмерный аргумент задаётся как \(x = g\cdot\mu_B\cdot J\cdot B / (k_B\cdot T)\) — это отношение магнитной энергии к тепловой. Данный калькулятор представляет собой чисто математический вычислитель специальной функции: вы напрямую задаёте \(x\) (без физических единиц), а он возвращает таблицу и график \(B_J(x)\). При очень больших J функция стремится к классической функции Ланжевена \(L(x)\), которую можно выбрать, введя «inf».
Как пользоваться калькулятором
Введите J в виде целого числа, полуцелой дроби (например, 1/2 или 3/2) либо десятичного значения вроде 0.5. Для предела Ланжевена укажите «inf». Затем задайте начальное значение \(x\) (Initial value of x), шаг между точками (Increment) и количество строк, которые нужно сгенерировать. Значения \(x\) вычисляются по формуле \(x_i = \text{startX} + i\cdot\text{stepX}\) при \(i\) от 0 до \(\text{count}-1\). Инструмент выводит каждую пару \((x, B_J(x))\) и строит кривую.
Разбор формулы
Для конечного J функция представляет собой комбинацию двух гиперболических котангенсов: $$B_J(x) = \frac{2J+1}{2J}\coth\!\left(\frac{2J+1}{2J}\,x\right) - \frac{1}{2J}\coth\!\left(\frac{x}{2J}\right)$$ Это нечётная функция, поэтому \(B_J(-x) = -B_J(x)\); она проходит через начало координат \((B_J(0)=0)\) и насыщается до \(\pm 1\) при \(x \to \pm\infty\). Поскольку \(\coth\) имеет особенность в нуле, калькулятор возвращает 0 для \(x\) в начале координат (или предельно близко к нему), что соответствует правильному аналитическому пределу.
Разобранный пример
Возьмём \(J = 1/2\) (тогда \(2J = 1\)) при \(x = 1\). Тогда \(\frac{2J+1}{2J} = 2\) и \(\frac{1}{2J} = 1\), что даёт $$B_{1/2}(1) = 2\cdot\coth(2) - \coth(1) = 2(1.037314) - 1.313035 = 0.761594.$$ Для проверки: при \(J = 1/2\) функция Бриллюэна совпадает с \(\tanh(x)\), а \(\tanh(1) = 0.761594\). Для случая Ланжевена при \(x = 2\): $$L(2) = \coth(2) - \frac{1}{2} = 1.037314 - 0.5 = 0.537314.$$
Интерпретация результата функции Бриллюэна
Значение, возвращаемое калькулятором, \(B_J(x)\), является безразмерным и ограничено между 0 и 1. Физически оно равно относительной намагниченности парамагнетика — отношению фактической намагниченности к намагниченности насыщения:
$$B_J(x) = \frac{M}{M_\text{sat}}, \qquad 0 \le B_J(x) \le 1.$$Результат 0 означает отсутствие чистого выравнивания магнитных моментов (нулевое поле или бесконечная температура), а результат, приближающийся к 1, означает, что каждый момент полностью выравнен с приложенным полем (полное насыщение).
Аргумент x: магнитная энергия и тепловая энергия
Входная переменная \(x\) — это отношение магнитной энергии (энергии Зеемана) момента к доступной тепловой энергии:
$$x = \frac{g\,\mu_B\,J\,B}{k_B\,T}.$$Когда \(x\) мало, случайная тепловая агитация \(k_B T\) доминирует над выравнивающей магнитной энергией, поэтому моменты почти рандомизированы; когда \(x\) велико, магнитная энергия берет верх и моменты фиксируются в выравнивании.
Режим малых x (режим Кюри)
Для \(x \ll 1\) функция Бриллюэна линейна по \(x\):
$$B_J(x) \approx \frac{J+1}{3J}\,x.$$Подставив \(x = g\mu_B J B /(k_B T)\), получим намагниченность, пропорциональную \(B/T\), что в точности соответствует закону Кюри: восприимчивость падает как \(1/T\). Это режим, применимый к обычным парамагнетикам в лабораторных полях при комнатной температуре, где \(B_J(x)\) обычно намного меньше 1.
Режим высоких x и насыщения
Для \(x \gg 1\) оба гиперболических котангенса стремятся к 1 и функция насыщается:
$$B_J(x) \to 1.$$Это соответствует сильным полям и/или очень низким температурам, где практически все магнитные моменты указывают вдоль поля и намагниченность больше не может увеличиваться. На графике это отображается как плато, приближающееся к горизонтальной линии \(B_J=1\). При \(J \to \infty\) кривая приближается к классической функции Ланжевена \(L(x)=\coth x - 1/x\).
Ключевые термины и переменные
| Символ / Термин | Значение |
|---|---|
| \(J\) | Квантовое число полного углового момента магнитного иона (объединяет орбитальный и спиновый вклады). Может быть целым или полуцелым числом (например, 1/2, 1, 3/2, 2). Оно определяет форму кривой и количество доступных состояний \(m_J\), равное \(2J+1\). |
| \(x\) | Безразмерный аргумент функции, \(x = g\mu_B J B/(k_B T)\) — отношение магнитной энергии (энергии Зеемана) к тепловой энергии. Это горизонтальная ось таблицы и графика. |
| \(g\) | Фактор Ланде (спектроскопический множитель расщепления), безразмерное число, связывающее магнитный момент с угловым моментом. Для чистого спина \(g \approx 2\); для объединенного орбитального и спинового углового момента оно определяется формулой Ланде. |
| \(\mu_B\) | Магнетон Бора, естественная единица атомного магнитного момента, \(\mu_B = e\hbar/(2m_e) \approx 9.274\times10^{-24}\ \text{Дж/Тл}\). |
| \(k_B\) | Постоянная Больцмана, \(k_B \approx 1.381\times10^{-23}\ \text{Дж/К}\), преобразующая температуру в тепловую энергию \(k_B T\). |
| \(B\) | Магнитная индукция (магнитное поле), измеряемая в тесла (Тл). Большее значение \(B\) увеличивает \(x\) и подталкивает систему к насыщению. |
| \(T\) | Абсолютная температура в кельвинах (К). Более высокая \(T\) увеличивает тепловую рандомизацию, уменьшая \(x\) и намагниченность. |
| \(\coth\) | Гиперболический котангенс, \(\coth(u) = \cosh(u)/\sinh(u) = (e^{u}+e^{-u})/(e^{u}-e^{-u})\); он появляется дважды в функции Бриллюэна и стремится к 1 при большом \(u\). |
| Функция Ланжевена \(L(x)\) | Классический предел функции Бриллюэна при \(J \to \infty\): \(L(x) = \coth x - 1/x\). Она описывает свободно вращающиеся классические магнитные диполи (без квантования ориентации). |
Часто задаваемые вопросы
Почему \(B_J(0)\) отображается как 0? Оба слагаемых с \(\coth\) расходятся при \(x = 0\), но их разность имеет конечный предел, равный 0; именно этот предел и показывает калькулятор.
Какие значения J допустимы? Положительные целые и полуцелые числа (1/2, 1, 3/2, 2, ...). Значение \(J = 0\) недопустимо, так как приводит к делению на ноль, а при вводе значений, не являющихся полуцелыми, калькулятор выдаёт предупреждение.
Как получить функцию Ланжевена? Введите «inf» (или «infinity») в поле J, чтобы использовать \(L(x) = \coth(x) - \frac{1}{x}\).
