[论文解读] On evaluation of integrals involving Bessel functions

D. Babusci, G. Dattoli|arXiv (Cornell University)|Nov 3, 2011

Algebraic and Geometric Analysis参考文献 1被引用 3

一句话总结

本文提出了一种符号算子方法，用于计算涉及贝塞尔函数、指数函数、对数函数及其乘积的定积分——这些在量子场论中至关重要。通过使用上下文演算将贝塞尔函数视为形式指数函数，该方法能够对复杂积分进行闭式求值，包括费曼图振幅中出现的积分，结果以伽马函数和 digamma 函数表示。

ABSTRACT

We introduce a symbolic method for the evaluation of definite integrals containing combinations of various functions, including exponentials, logarithm and products of Bessel functions of different types. The method we develop is naturally suited for the evaluation of integrals associated with specific Feynman diagrams.

研究动机与目标

• 开发一种系统性的符号方法，用于评估包含贝塞尔函数、指数函数和对数函数的积分。
• 通过减去同阶无穷大来处理发散积分的挑战，从而获得有限结果。
• 将该方法扩展至不同类型贝塞尔函数的乘积，以及修正贝塞尔函数（Macdonald 函数）和诺伊曼函数（Neumann functions）。
• 实现对量子场论中费曼图振幅相关积分的高效计算。
• 为使用计算代数系统实现自动化符号计算奠定基础。

提出的方法

• 利用由 $\hat{c}^\mu \varphi(0) = \varphi(\mu) = \frac{1}{\Gamma(\mu+1)}$ 定义的上下文算子 $\hat{c}$，将贝塞尔函数表示为形式指数函数。
• 应用恒等式 $J_\nu(x) = \left(\frac{x}{2}\hat{c}\right)^\nu \exp\left\{-\hat{c}\left(\frac{x}{2}\right)^2\right\} \varphi(0)$，以代数方式表达贝塞尔函数。
• 通过发散积分的正则化程序 $\lim_{\nu \to -1} \int x^\nu dx = \ln x$ 处理由对数项引起的发散。
• 通过 $K_n(x) = \frac{\pi}{2} \lim_{\nu \to n} \frac{I_{-\nu}(x) - I_\nu(x)}{\sin(\nu\pi)}$ 将方法扩展至修正贝塞尔函数，并使用上下文算子重新表述。
• 通过引入广义算子 $\hat{d}$ 表示 $J_\mu(ax)J_\nu(bx)$，将符号方法应用于贝塞尔函数乘积，从而导出级数表示。
• 利用广义超几何型函数，推导出如 $\int_0^\infty x^\alpha J_\mu(ax)J_\nu(bx) dx$ 等积分的闭式表达式。

实验结果

研究问题

• RQ1 如何对涉及 $x^\nu$ 的发散积分进行正则化，以获得有限结果（如 $\ln x$）？
• RQ2 能否通过符号算子将贝塞尔函数视为初等函数，以简化积分求值？
• RQ3 积分 $\int_0^\infty x^\mu J_\lambda(px) \ln(bx) dx$ 的闭式表达式是什么？
• RQ4 如何对不同类型贝塞尔函数的乘积进行符号积分？
• RQ5 该方法能否扩展至量子场论中涉及修正贝塞尔函数和诺伊曼函数的积分？

主要发现

• 通过正则化方法推导出积分 $\int_0^\infty x^\mu e^{-a x^2} \ln(bx) dx = \frac{1}{4\sqrt{a^{\mu+1}}} \Gamma\left(\frac{\mu+1}{2}\right) \left[\ln\left(\frac{b^2}{a}\right) + \psi\left(\frac{\mu+1}{2}\right)\right]$。
• 通过符号方法获得表达式 $A_\lambda(p,b) = \int_0^\infty J_\lambda(px) \ln(bx) dx = \frac{1}{p} \left[\ln\left(\frac{2b}{p}\right) + \psi\left(\frac{\lambda+1}{2}\right)\right]$。
• 乘积 $J_\mu(ax)J_\nu(bx)$ 被表示为 $\left(\frac{a x}{2}\right)^\mu \left(\frac{b x}{2}\right)^\nu \exp\left\{-\hat{d} \left(\frac{x}{2}\right)^2\right\} \varphi_{\mu,\nu}(0;a,b)$，从而实现系统的积分计算。
• 积分 $\Omega_{\alpha,\mu,\nu}(a,b) = \int_0^\infty x^\alpha J_\mu(ax)J_\nu(bx) dx$ 被计算为 $2^\alpha a^\mu b^\nu \Gamma\left(\frac{\alpha+\mu+\nu+1}{2}\right) \varphi_{\mu,\nu}\left(-\frac{\alpha+\mu+\nu+1}{2};a,b\right)$。
• 该方法成功计算出 $\Xi_{\alpha,\mu,\nu}(a,b) = \int_0^\infty x^\alpha J_\mu(ax)Y_\nu(bx) dx$，使用公式 $\Xi = \frac{1}{\sin(\nu\pi)} \left[ \cos(\nu\pi)\Omega_{\alpha,\mu,\nu} - \Omega_{\alpha,\mu,-\nu} \right]$。
• 利用 $\hat{c}^\nu e^{-\hat{c}x} \varphi(0)$ 表达了三类贝塞尔函数 $C_\nu(x)$，其积分为 $\int_0^\infty C_\nu(x) dx = \frac{1}{\Gamma(\nu)}$ 和 $\int_0^\infty x^\mu C_\nu(x^2) dx = \frac{1}{2} \frac{\Gamma\left(\frac{\mu+1}{2}\right)}{\Gamma\left(\nu - \frac{\mu-1}{2}\right)}$。