[論文レビュー] 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つの確率測度間の相互エネルギーの級数展開を適用する。
• エントロピー計算を再帰関係の係数と特殊関数を含む和に変換する。
• 非整数パラメータに対し、切り捨て級数と補間を用いてギーゲンバウサー多項式用のアルゴリズム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} $ のエントロピーを用いて表現可能であり、アルゴリズムの直接的適用が可能になる。
• パラメータ $ u $ が小さい範囲では、$ C_n^ u $ の非整数 $ u $ におけるエントロピー値の補間は大きな誤差を生じ、その範囲での信頼性が制限される。
• アルゴリズムは効率的かつスケーラブルであり、少数の基本点を事前に計算した上で補間を用いることで、$ E_{200}^ u $ の計算時間も競争力を持つ。