[論文レビュー] Analytical properties of the R^(1/m) luminosity law

Luca Ciotti, G. Bertin|arXiv (Cornell University)|Nov 5, 1999

Galaxies: Formation, Evolution, Phenomena参考文献 5被引用数 24

ひとこと要約

この論文は、Sersic $ R^{1/m} $則における次元なしスケール因子 $ b(m) $ に対する非常に高精度な漸近展開を導出しており、$ m \geq 1 $ の全範囲で正確な解析的計算を可能にする。導出された式 $ b(m) = 2m - \frac{1}{3} + \frac{4}{405m} + \frac{46}{25515m^2} $ は、$ m=1 $ でさえも相対誤差が $ 10^{-6} $ 未満に抑えられ、従来の近似式を上回り、銀河の光度プロファイルにおける理論的・観測的応用を統合する。

ABSTRACT

In this paper we describe some analytical properties of the R^{1/m} law proposed by Sersic (1968) to categorize the photometric profiles of elliptical galaxies. In particular, we present the full asymptotic expansion for the dimensionless scale factor b(m) that is introduced when referring the profile to the standard effective radius. Surprisingly, our asymptotic analysis turns out to be useful even for values of m as low as unity, thus providing a unified analytical tool for observational and theoretical investigations based on the R^{1/m} law for the entire range of interesting photometric profiles, from spiral to elliptical galaxies.

研究の動機と目的

• Sersic $ R^{1/m} $則における次元なしスケール因子 $ b(m) $ の閉形式漸近展開を導出すること。これは、一般に数値的にしか解けない。
• $ m=1 $（指数的）から $ m=10 $（de Vaucouleurs型に近い）までの全範囲で有効な単一の解析的ツールを提供することで、理論的および観測的取り扱いを統合すること。
• C89、C91、PS97のような既存の補間式を、特に低 $ m $ における精度を高めるために上回ること。
• 大 $ m $ におけるSersicプロファイルの漸近的収束が $ R^{-2} $ のべき乗則に近づくこと、およびその近似が成り立つ半径範囲を明確にすること。

提案手法

• ガンマ関数のスターリング近似を用いて、$ \gamma(\alpha, x) = \Gamma(\alpha)/2 $ の解 $ x $ に対する漸近展開を導出する。ここで $ \alpha = 2m $ である。
• 漸近近似を段階的に改善するための列 $ x_n = \alpha + \sum_{k=0}^{n-1} \frac{c_k}{\alpha^k} $ を構築する。
• 不完全ガンマ関数 $ \gamma(\alpha, x) $ の漸近展開を用いて係数 $ c_k $ を導出し、$ m $ の逆数の累乗の級数を得る。
• 実用的な公式を得るため、展開を4項で切り捨て、$ b(m) = 2m - \frac{1}{3} + \frac{4}{405m} + \frac{46}{25515m^2} $ を得る。
• 数値的解 $ \gamma(2m, b) = \Gamma(2m)/2 $ との比較により、$ 1 \leq m \leq 10 $ の範囲で公式の妥当性を検証する。
• 大 $ m $ におけるSersicプロファイルの半径的挙動を、拡張座標 $ \xi = \ln \eta $ を導入することで分析し、$ R^{-2} $ べき乗則と比較する。

実験結果

リサーチクエスチョン

• RQ1 Sersic $ R^{1/m} $則における $ b(m) $ の閉形式漸近展開を導出し、$ m=1 $ のような小さな $ m $ に対しても精度が保たれるようにできるか。
• RQ2 大 $ \alpha $ における不完全ガンマ関数 $ \gamma(\alpha, x) $ の漸近的性質は、$ \gamma(2m, b) = \Gamma(2m)/2 $ の解をどのように示唆するか。
• RQ3 Sersicプロファイル $ I(R) \propto e^{-b \eta^{1/m}} $ が $ R^{-2} $ べき乗則に近づく半径範囲は何か。また、その誤差は $ m $ にどのように依存するか。
• RQ4 既存の $ b(m) $ の補間式（特に $ m=1 $ の場合）と比較して、新しい漸近展開はどの程度精度が優れているか。
• RQ5 切り捨てられた漸近展開の最大相対誤差は何か。また、$ m \geq 1 $ の全範囲で $ 10^{-6} $ 未満に保たれるか。

主な発見

• 4項に切り捨てた導出された漸近展開 $ b(m) = 2m - \frac{1}{3} + \frac{4}{405m} + \frac{46}{25515m^2} $ は、$ m \geq 1 $ の全範囲で相対誤差が $ 10^{-6} $ 未満である。$ m=1 $ に対しても同様に成立する。
• $ m=1 $ では、切り捨て展開の相対誤差は7.3%であり、$ m=4 $ では1.8%である。これは、$ b \simeq 2m - 1/3 $ のような従来の近似式よりも優れた精度を示している。
• 通常の漸近級数の期待とは反して、この展開は小さな $ m $ に対しても有効であることが示され、$ m=1 $ から $ m=10 $ までのすべてのSersicプロファイルを統一的に扱える解析的ツールとなる。
• Sersicプロファイルは大 $ m $ において $ R^{-2} $ べき乗則に漸近的に近づく。この近似は、$ \eta \lesssim e^{1/3} \approx 1.39 $ の半径範囲で最も正確であり、相対誤差が最小化される。
• $ R^{-2} $ 近似は、$ 1 < \eta < e^{1/3} $ の領域でSersicプロファイルを低く見積もる。最大相対誤差は $ \sim \frac{1}{36m} $ のオーダーであり、$ m $ が大きい場合には小さい。
• 相対誤差 $ \epsilon $ 以内で $ R^{-2} $ 近似が成り立つ半径範囲は、$ (1-3\epsilon)^m \lesssim \eta \lesssim (1+3\epsilon)^m e^{1/3} $ で境界づけられ、非対称な適用範囲を示している。