什么是朗之万函数?
朗之万函数定义为 \(L(x) = \coth(x) - \frac{1}{x}\),其中 \(\coth(x)\) 是双曲余切函数。它最早出现在保罗·朗之万(Paul Langevin)的经典顺磁性理论中,用来描述在外磁场作用下、一群可自由转动的磁偶极子的平均磁化强度。同样的函数也出现在介电极化的朗之万-德拜理论中,以及自由连接高分子链的统计力学里——在那里它把链的伸长量与施加的张力联系起来。从数学上看,它是布里渊函数(Brillouin function)在 \(J\) 趋于无穷大时的极限。
如何使用本计算器
在自变量 \(x\) 处输入任意实数——正数、负数或零均可,计算器即可返回 \(L(x)\) 的值。该函数为奇函数,因此 \(L(-x) = -L(x)\),其输出值始终严格落在 -1 与 1 之间。当 \(x\) 取较大值时,结果会逐渐饱和趋近于 +1(\(x\) 为较大负数时则趋近于 -1),这反映了偶极子的完全取向排列。
公式详解
$$L(x) = \coth(x) - \frac{1}{x}$$其中 \(\coth(x) = \frac{\cosh(x)}{\sinh(x)} = \frac{e^x + e^{-x}}{e^x - e^{-x}}\)。在 \(x = 0\) 处两项各自发散,但它们的差存在一个可去奇点,其值等于 0。由于在零点附近直接相减两个很大的数会引发灾难性抵消(catastrophic cancellation),因此当 \(|x|\) 非常小时,本工具会改用泰勒级数 \(L(x) \approx \frac{x}{3} - \frac{x^3}{45} + \frac{2x^5}{945}\) 来计算,由此给出原点处众所周知的斜率 \(\frac{1}{3}\)。
计算示例
当 \(x = 1\) 时:\(\cosh(1) = 1.5430806348\),\(\sinh(1) = 1.1752011936\),故 \(\coth(1) = 1.3130352855\)。于是 $$L(1) = 1.3130352855 - 1 = 0.3130352855$$对于较小的自变量 \(x = 0.1\),级数计算给出 \(\frac{0.1}{3} - \frac{0.001}{45} \approx 0.0333111\),与直接计算的结果一致。
常见问题
为什么 \(L(0) = 0\)? \(\coth(x)\) 和 \(\frac{1}{x}\) 在 \(x = 0\) 处都会趋于无穷大,但二者恰好相互抵消;它们差值的极限正好是 0,并且函数从原点起以斜率 \(\frac{1}{3}\) 线性上升。
\(L(x)\) 的取值范围是多少? 朗之万函数单调递增,值域为 \((-1, 1)\);它无限趋近于 \(\pm 1\) 这两条渐近线,但永远无法真正达到。
它与布里渊函数有什么关系? 当 \(J\) 趋于无穷大(即经典连续自旋的情形)时,布里渊函数 \(B_J(x)\) 退化为朗之万函数。