索末菲-库默尔函数的渐近展开(大参数情况)
1. 背景与问题引入
索末菲-库默尔函数是微分方程
\[\frac{d^2 w}{dz^2} + \left( \frac{1}{z} \right) \frac{dw}{dz} + \left( 1 - \frac{\nu^2}{z^2} \right) w = 0 \]
的解,常用 \(F_\nu(z)\) 表示。在物理问题中(如波传播、衍射分析),常需研究 \(|z| \gg 1\) 或 \(|z| \gg |\nu|\) 时的渐近行为。大参数渐近展开旨在用简单的解析近似代替复杂函数,便于计算和分析。
2. 渐近展开的核心思想
渐近展开不要求级数收敛,而是要求截断误差在参数趋于无穷时迅速衰减。对于大 \(|z|\),索末菲-库默尔函数的渐近形式通常包含指数项、代数衰减项和振荡项,其结构为:
\[F_\nu(z) \sim e^{\pm i z} z^{\alpha} \sum_{k=0}^\infty \frac{a_k(\nu)}{z^k} \quad (|z| \to \infty). \]
其中符号 \(\sim\) 表示渐近等价,即截断到 \(z^{-N}\) 的误差为 \(o(|z|^{-N})\)。
3. 推导方法:鞍点法与积分表示
索末菲-库默尔函数有积分表示:
\[F_\nu(z) = \frac{1}{\sqrt{\pi}} \int_C e^{i z \cosh t - \nu t} \, dt, \]
路径 \(C\) 取决于边界条件。大 \(z\) 时,指数因子 \(e^{i z \cosh t}\) 剧烈振荡,主要贡献来自鞍点(驻相点)。求解 \(\frac{d}{dt}(i z \cosh t) = 0\) 得鞍点 \(t = 0\)。通过鞍点法展开被积函数,可得渐近级数:
- 在鞍点附近展开 \(\cosh t \approx 1 + \frac{t^2}{2}\);
- 高斯积分给出主导项 \(e^{i z} / \sqrt{z}\);
- 高阶修正通过展开到 \(t^4, t^6, \ldots\) 并逐项积分得到。
4. 具体渐近公式
索末菲-库默尔函数的大参数渐近展开分为两类(对应线性无关解):
\[F_\nu(z) \sim \frac{e^{i z}}{\sqrt{z}} \sum_{k=0}^\infty \frac{a_k(\nu)}{(i z)^k} + \frac{e^{-i z}}{\sqrt{z}} \sum_{k=0}^\infty \frac{b_k(\nu)}{(-i z)^k}, \]
其中系数 \(a_k, b_k\) 由递归关系确定。例如,首项为:
\[F_\nu(z) \approx \sqrt{\frac{2}{\pi z}} \cos\left( z - \frac{\nu\pi}{2} - \frac{\pi}{4} \right) + O(|z|^{-3/2}). \]
该公式在 \(|\arg z| < \pi\) 时成立,且需注意斯托克斯现象(不同区域相位跳变)。
5. 误差分析与有效性
- 截断误差:若取前 \(N\) 项,误差量级为 \(|z|^{-N-1/2}\)。
- 有效性条件:要求 \(|z| \gg |\nu|\) 且 \(|z| \gg 1\),否则需改用小参数展开或其他方法。
- 复数域扩展:通过解析延拓,公式可推广到复平面,但需考虑分支切割(通常沿负实轴)。
6. 物理应用示例
在波导或散射问题中,索末菲-库默尔函数描述远离源点的场分布。例如,当 \(z = kr\)(\(k\) 为波数,\(r\) 为距离),大 \(r\) 渐近展开给出辐射场的球面波衰减形式:
\[\psi(r) \propto \frac{e^{i k r}}{r} \left( 1 + \frac{c}{r} + \cdots \right), \]
其中 \(c\) 由 \(a_k(\nu)\) 决定,反映角动量 \(\nu\) 对远场相位的影响。
总结
索末菲-库默尔函数的大参数渐近展开通过鞍点法系统导出,结合指数振荡与代数衰减,为波传播、量子散射等物理问题提供了高效计算工具。其精度随 \(|z|\) 增大而提高,但需注意适用条件和复平面上的解析行为。