什么是布里渊函数？
布里渊函数 \(B_J(x)\) 用来描述由总角动量量子数为 \(\text{J}\) 的原子组成的顺磁体的磁化强度。在统计力学中，它的无量纲自变量为 \(x = g\cdot\mu_B\cdot\text{J}\cdot B / (k_B\cdot T)\)，也就是磁能与热能之比。本计算器是一个纯数学的特殊函数求值工具：你直接输入 \(x\)（不涉及任何物理单位），它就会返回 \(B_J(x)\) 的数值表和图像。当 \(\text{J}\) 取值非常大时，该函数会趋近于经典的朗之万函数 \(L(x)\)，只要输入“inf”即可选用这一极限。
如何使用本计算器
\(\text{J}\) 可以填整数、半整数（如 1/2、3/2），也可以填小数（如 0.5）。若要取朗之万极限，请输入“inf”。接着设定第一个 \(x\) 值（\(x\) 的初始值）、相邻取值点之间的间隔（步长），以及要生成多少行。各个 \(x\) 值按 \(x_i = \text{startX} + i\cdot\text{stepX}\) 生成，其中 \(i\) 从 0 取到 \(\text{count}-1\)。工具会逐一列出每一对 \((x, B_J(x))\)，并绘制出对应曲线。
公式解析
当 \(\text{J}\) 为有限值时，函数由两个双曲余切相减构成：
$$B_J(x) = \frac{2\text{J}+1}{2\text{J}}\coth\!\left(\frac{2\text{J}+1}{2\text{J}}\,x\right) - \frac{1}{2\text{J}}\coth\!\left(\frac{x}{2\text{J}}\right)$$它是奇函数，因此 \(B_J(-x) = -B_J(x)\)；曲线过原点（\(B_J(0)=0\)），并在 \(x \to \pm\infty\) 时饱和趋向 \(\pm 1\)。由于 \(\coth\) 在零点存在奇点，计算器在 \(x\) 等于（或极其接近）原点时返回 0，这正是该函数正确的解析极限。
实例演算
取 \(\text{J} = 1/2\)（即 \(2\text{J} = 1\)），在 \(x = 1\) 处计算。此时 \((2\text{J}+1)/2\text{J} = 2\)、\(1/2\text{J} = 1\)，于是 $$B_{1/2}(1) = 2\cdot\coth(2) - \coth(1) = 2(1.037314) - 1.313035 = 0.761594$$可作验证：当 \(\text{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因子（谱线裂分因子），是一个无量纲数，将磁矩与角动量联系起来。对于纯自旋\(g \approx 2\)；对于轨道和自旋角动量的组合，由朗德公式给出。 |
| \(\mu_B\) | 玻尔磁子，原子磁矩的自然单位，\(\mu_B = e\hbar/(2m_e) \approx 9.274\times10^{-24}\ \text{J/T}\)。 |
| \(k_B\) | 玻尔兹曼常数，\(k_B \approx 1.381\times10^{-23}\ \text{J/K}\)，将温度转换为热能\(k_B T\)。 |
| \(B\) | 磁感应强度（磁场），单位为特斯拉（T）。较大的\(B\)增加\(x\)并将系统驱向饱和。 |
| \(T\) | 绝对温度，单位为开尔文（K）。较高的\(T\)增加热随机化，减小\(x\)和磁化。 |
| \(\coth\) | 双曲余切，\(\coth(u) = \cosh(u)/\sinh(u) = (e^{u}+e^{-u})/(e^{u}-e^{-u})\)；它在布里渊函数中出现两次，对于大的\(u\)趋向于1。 |
| 朗之万函数\(L(x)\) | 当\(J \to \infty\)时布里渊函数的经典极限：\(L(x) = \coth x - 1/x\)。它描述自由旋转的经典磁偶极子（方向的量子化）。 |
常见问题
为什么 \(B_J(0)\) 显示为 0？在 \(x = 0\) 处两个 \(\coth\) 项都会发散，但它们的差值具有有限极限 0，工具返回的正是这个极限值。
\(\text{J}\) 可以取哪些值？正整数和半整数（1/2、1、3/2、2、……）。\(\text{J} = 0\) 无效，因为会导致除以零；若输入的不是半整数，计算器会给出提示。
怎样得到朗之万函数？将 \(\text{J}\) 设为“inf”（或“infinity”），即可使用 $$L(x) = \coth(x) - \frac{1}{x}$$
