ブリルアン関数とは
ブリルアン関数 \(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 の値（x の初期値）、点の間隔（増分）、生成する行数を指定します。x の値は \(x_i = \text{startX} + i\cdot\text{stepX}\)（i = 0 〜 count−1）の式で生成されます。本ツールはすべての \((x, B_J(x))\) の組を表示し、曲線を描画します。
計算式の解説
有限の J に対して、関数は2つの双曲線余接（coth）の組み合わせで表されます：
$$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\)）、\(x \to \pm\infty\) で \(\pm 1\) に飽和します。coth は原点で特異点をもつため、x が原点（またはごく近傍）にある場合は、解析的な極限値である 0 を返します。
計算例
J = 1/2（すなわち 2J = 1）で x = 1 のときを考えます。このとき \((2J+1)/2J = 2\)、\(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 因子（分光的スプリット因子）、角運動量に対する磁気モーメントを関連付ける無次元数。純スピンの場合 \(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})\)；これはブリルアン関数に 2 回現れ、大きい \(u\) に対して 1 に近づきます。 |
| ランジュヴァン関数 \(L(x)\) | \(J \to \infty\) のときのブリルアン関数の古典的極限：\(L(x) = \coth x - 1/x\)。自由に回転する古典的磁気双極子（配向の量子化なし）を記述します。 |
よくある質問
なぜ \(B_J(0)\) は 0 と表示されるのですか？ x = 0 では両方の coth 項が発散しますが、その差は有限の極限値 0 をもちます。本ツールはその極限値を表示します。
J に有効な値は何ですか？ 正の整数および半整数（1/2、1、3/2、2、…）です。J = 0 はゼロ除算となるため無効で、半整数以外の入力に対しては警告が表示されます。
ランジュバン関数を求めるには？ J に「inf」（または「infinity」）と入力すると、\(L(x) = \coth(x) - \frac{1}{x}\) が使用されます。
