[论文解读] Computation of the entropy of polynomials orthogonal on an interval

V S Buyarov, J. S. Dehesa+2 more|ArXiv.org|Oct 16, 2003

Quantum Mechanics and Non-Hermitian Physics被引用 6

一句话总结

本文提出了一种高效算法，仅使用正交多项式三系数递推关系的系数，即可计算实区间上正交多项式的信息熵。该方法利用概率测度之间互能的级数表达式，实现了对勒让德多项式及其在量子系统中球谐函数应用的精确计算。

ABSTRACT

We give an effective method to compute the entropy for polynomials orthogonal on a segment of the real axis that uses as input data only the coefficients of the recurrence relation satisfied by these polynomials. This algorithm is based on a series expression for the mutual energy of two probability measures naturally connected with the polynomials. The particular case of Gegenbauer polynomials is analyzed in detail. These results are applied also to the computation of the entropy of spherical harmonics, important for the study of the entropic uncertainty relations as well as the spatial complexity of physical systems in central potentials.

研究动机与目标

• 开发一种数值稳定且精确的算法，用于计算有限区间上正交多项式的熵。
• 克服因多项式零点处奇点导致的直接数值积分方法的数值不稳定性。
• 实现对由正交多项式控制的系统熵的实用计算，特别是在量子力学中的应用。
• 将该方法扩展至通过球谐函数计算量子系统的空间复杂度，利用勒让德多项式的熵。
• 提供适用于中心势场系统（如类氢原子和双原子分子）的计算框架。

提出的方法

• 以三系数递推关系的系数作为输入，避免对奇异函数进行直接积分。
• 应用从正交多项式导出的两个概率测度之间互能的级数展开。
• 将熵的计算转化为涉及递推关系系数和特殊函数的求和形式。
• 通过截断级数和插值方法，实现对勒让德多项式的算法2计算，适用于非整数参数。
• 通过已知的归一化和角向依赖公式，将球谐函数的熵与勒让德多项式的熵关联起来。
• 在计算基础情况后，利用三次样条插值处理非整数参数值，以在精度与计算成本之间取得平衡。

实验结果

研究问题

• RQ1 是否可以仅通过递推系数而无需直接数值积分，高效计算实区间上正交多项式的熵？
• RQ2 如何利用三系数递推关系的系数计算正交多项式的熵？
• RQ3 勒让德多项式的熵在不同参数值和次数下的行为如何？
• RQ4 球谐函数的熵如何随量子数 $ l $ 和 $ m $ 变化，这又对空间复杂度意味着什么？
• RQ5 插值策略能否可靠地估计勒让德多项式中非整数参数下的熵？

主要发现

• 所提出的算法仅使用递推系数即可实现对多项式熵的精确计算，避免了奇异被积函数带来的数值问题。
• 对于勒让德多项式，熵的计算是可行且稳定的，结果通过级数截断和插值得到验证。
• 球谐函数 $ S_{l,m}[Y] $ 的熵在 $ m = 0 $ 附近最高，并随 $ |m| $ 接近 $ l $ 而单调递减，表明在高 $ |m| $ 时空间局域性更强。
• 球谐函数 $ Y_{l,m} $ 的熵可表示为勒让德多项式 $ E_{l-|m|}^{|m|+1/2} $ 的熵，从而可直接应用该算法。
• 对 $ C_n^ u $ 中非整数 $ u $ 的熵值进行插值，在 $ u $ 较小时会产生较大误差，限制了该方法在该参数区间的可靠性。
• 该算法具有高效性和可扩展性，当使用少量计算得到的基点进行插值时，$ E_{200}^ u $ 的计算时间具有竞争力。