A four parameters modified Gaussian function describing shape of asymmetric alpha lines

技術要約：非対称アルファ線のための4パラメータ修正ガウス関数

問題提起
アルファ線スペクトロメトリーにおいて、スペクトル線は理想的なガウス分布ではなく、頻繁に非対称な形状を示す。この非対称性は、主にアルファ粒子が半導体検出器のデッドレイヤー（入射窓）や試料物質を通過する際のエネルギー減衰に起因する。この減衰の程度は、検出器に対する粒子のランダムな放出角度に依存する。その結果、ピークの低エネルギー側は幅広くなり、指数関数的な減衰に近づく一方で、高エネルギー側は統計的なゆらぎによる修正のみを受けたガウス分布を維持する。

歴史的に、この形状を記述するための単一の、一般的に受け入れられた数学的公式は存在しなかった。実用的な解決策としては、スペクトル線の特定のセクションに対してマルチパラメータの関数セットを適用する方法が一般的であった。様々なソフトウェアパッケージ（例：ALFUN-A、MicroSAMPO、ALFA）が存在するが、広く採用されているALFコンピュータコードで使用されている特定の解析関数は、本研究以前には文献に詳細が記載されていなかった。

手法
本論文は、わずか4つの自由パラメータを用いて非対称なアルファピーク形状を記述するために設計された、独自の解析関数を提示している。この関数は、1993年頃に開発された（以前のガンマ線用コードであるP.I.M.P.から修正された）アルファスペクトル評価用のFortranベースのツールであるALFコンピュータコードに実装されている。

手法の核となるのは、以下の修正ガウス関数 $N(x)$ である。ここで $x$ はチャンネル番号を表す。

$$N(x) = \alpha \cdot \exp\left( \frac{-(x-M)^2}{\beta \{|x-M| - (x-M)\} + \gamma} \right)$$

パラメータ：

• $\alpha$: 振幅（最大値におけるカウント数）。
• $M$: ピーク最大値の位置。
• $\beta$ および $\gamma$: 試料および検出器の幾何学的構造に固有の形状パラメータ（ともに $>0$）。

関数の挙動：

• 高エネルギー側 ($x > M$): 関数は標準的なガウス分布に収束する。
• 低エネルギー側 ($x < M$): 関数はより複雑になり、ピーク最大値から大きく離れた位置（$x \ll M$）では通常の指数関数へと移行する。
• パラメータの役割: $\gamma$ はピークの右側の形状を制御し、$\beta$ と $\gamma$ は共に左側の形状を制御する。

フィッティング手順（デコンボリューション）:
ALFコードは、最大10個の重なり合ったアルファピークをデコンボリューションするために、非線形関数フィッティング手順を利用する。プロセスは以下の2つの明確なステップで実行される。

1. 初期形状の決定: 単一の明確なスペクトル線（または最大値を含む断片）をフィッティングして、4つのパラメータ（$\alpha, M, \beta, \gamma$）を決定する。これには、線形近似（テイラー展開）を用いた10回の反復による標準的な非線形最小二乗法が用いられる。
2. ピークのデコンボリューション: $\beta$ と $\gamma$ が決定されると、これらは個々のピークではなく、試料の品質および試料・検出器の幾何学的構造に依存するため、全スペクトルに対して固定定数として扱われる。第2ステップでは、すべての重なり合ったピークの振幅（$\alpha_i$）と位置（$M_i$）を同時にフィッティングする。これは、関数の和を級数へと展開し、クラコビアン（Cracovian）アルゴリズムを用いて変数ベクトルの増分を解くことで達成される。

クラコビアンアルゴリズムは、パラメータ値の決定に使用されるプロセス内で、適合パラメータの不確かさも同時に計算できるという具体的な利点として強調されている。

主な貢献と結果

• 解析的な簡潔性: 本論文は、複雑なマルチセグメント・フィッティングを必要とせずに、アルファ線の複雑な非対称性を効果的にモデル化するコンパクトな4パラメータ関数を導入している。
• 検証と長期的な有用性: この関数は1993年から日常的に使用されている。過去30年間にわたり、15,000以上のスペクトルに適用されてきた。検証は、計算されたピーク面積と総カウント数の比較によって行われ、すべてのテストが肯定的に合格した。
• デコンボリューション能力: 本手法は、動物の骨から得られたTh分画試料を用いて示されたように、部分的に重なり合ったピークを分離することに成功した。
• 予期せぬ有用性: 著者らは、対数チャンネルスケールを用いた液体シンチレーションスペクトロメータによるベータ線スペクトル（具体的には $^{90}\text{Sr}$ / $^{90}\text{Y}$）の解析においても、この関数が驚くほど有用であることを指摘しており、誤差3%以内の結果を得ている。

意義
本論文の意義は、主にALFソフトウェア内におけるこの特定の4パラメータ関数の長期的な検証成功にある。これまで文書化されていなかったALFの解析的基礎を形式化することで、本研究は、重なり合ったアルファピークをデコンボリューションするための堅牢で数学的に定義されたツールを提供している。本手法は、解析的な厳密さと計算効率のバランスを実現しており、従来のガウスフィッティングが失敗する複雑なスペクトルにおいて、ピーク面積の精密な決定を可能にしている。クラコビアンアルゴリズムの使用により、フィッティングと同時にパラメータの不確かさが厳密に計算されることが保証されている。