什麼是布里淵函數?
布里淵函數 \(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 到 \(\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\)),並在 \(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:磁能 vs 熱能
輸入 \(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?兩個 \(\coth\) 項在 \(x = 0\) 時都會發散,但兩者的差具有有限極限 0;本工具回傳的正是此極限值。
\(J\) 可以輸入哪些值?正整數與半整數(\(1/2\)、\(1\)、\(3/2\)、\(2\)、……)。\(J = 0\) 無效,因為會發生除以零的情形;若輸入非半整數的值,計算器也會提出警告。
如何取得朗之萬函數?將 \(J\) 輸入為「inf」(或「infinity」),即可使用 \(L(x) = \coth(x) - \frac{1}{x}\)。