Excitation function measurement of $^{144}$Sm($α$,n) reaction at sub-Coulomb energies and detailed covariance analysis

本論文では、分子堆積法により作製された高純度アルミ箔上の濃縮$^{144}$Sm$_2$O$_3$標的を用いた積層箔照射法により、Coulomb 障壁以下のエネルギー領域（14〜21 MeV）における$^{144}$Sm($\alpha$,n)反応の断面積を測定し、ビームエネルギーの不確かさのシミュレーションや共分散・相関行列を用いた詳細な誤差解析を行った上で、既存の実験データおよびハウザー・ファッシュバック統計モデルによる理論予測と比較検討した。

要約

1. 研究の目的：宇宙のレシピと、新しい薬の材料

この実験は、大きく分けて 2 つの重要な目的を持っています。

• 宇宙のレシピ（天体物理学）：
  宇宙には「p 核（ピー核）」と呼ばれる、非常に珍しい元素があります。これらは通常の星の進化では作られず、超新星爆発のような激しい出来事でしか作られない「宇宙の希少食材」です。
  この研究では、「サマリウム（Sm）」という元素にアルファ粒子をぶつけて、どうやって別の元素（ガドリニウム）に変えるかを調べました。これは、宇宙でこれらの珍しい元素が「どうやって作られたか」というレシピ（反応の確率）を解明するための重要な一歩です。

• 新しい薬の材料（医療）：
  実験で作られた「ガドリニウム -147」という物質は、がんの診断に使われる**SPET（単一光子放射断層撮影）という検査に使える可能性があります。
  ちょうど「カメラのフィルム」のように、体内の特定の場所を鮮明に写し出すのに最適な半減期（38 時間）と、強い光（ガンマ線）を出すため、「魔法の薬（マグネトファーマシューティカル）」**として期待されているのです。

2. 実験の方法：「重ねたパン」に弾を撃つ

研究者たちは、以下のような方法で実験を行いました。

• 積み重ねたフォイル（紙）の山：
  薄いアルミの板の上に、サマリウムという物質を非常に薄く（髪の毛の太さの数千分の一程度）塗った「ターゲット」を、何枚も積み重ねました。これを**「積み重ねたパン」**と想像してください。
• エネルギーを調整する：
  加速器から発射されたアルファ粒子（弾丸）は、最初とても速い（28 MeV）ですが、それでは「パン」の奥まで届きすぎてしまいます。そこで、アルミの板を挟んで**「減速装置」**を通すことで、弾丸のスピードを徐々に落とし、14〜21 MeV という「ちょうどいい速度」で各層に当てるようにしました。
• 反応を待つ：
  弾が当たると、サマリウムがアルファ粒子を飲み込み、中性子を吐き出して「ガドリニウム」に変わります。これを**「活性化」**と呼びます。

3. データの解析：「料理の味」を正確に測る

実験が終わると、できたガドリニウムがどんな光（ガンマ線）を出しているかを測りました。しかし、ここにはいくつかの「落とし穴」がありました。

• エネルギーのブレ（ジャストインタイムの難しさ）：
  弾丸が「パン」の層を通過するたびに、少しずつスピードが落ち、バラつき（エネルギーの広がり）が生まれます。
  研究者たちは、**「GEANT4」というコンピューターシミュレーション（仮想空間での再現）を使って、弾丸がどの層で、どのくらいの速度で止まったかを精密に計算しました。まるで、「風が吹く中で、何枚も傘を重ねて通るボールの軌跡を、スーパーコンピュータで追跡する」**ような作業です。
• 誤差の計算（ covariance 分析）：
  通常、実験結果には「誤差（±）」がつきます。しかし、この研究では**「誤差同士の関係性（共分散）」**まで詳しく分析しました。
  • 例え話： 「料理の味」を測る時、塩の量、火の強さ、材料の鮮度がそれぞれ誤差を持っています。これらは独立ではなく、**「塩を多く入れれば火の強さの調整も変わる」**ように、互いに影響し合っています。
  • この研究では、**「どの誤差が、どの誤差と手を取り合っているか」を行列（表）にまとめて、結果の信頼性を極限まで高めました。これはこの反応の測定において「初めて」**行われた詳細な分析です。

4. 理論との比較：「レシピ本」でチェック

実験で得られた結果（どのくらいの確率で反応が起きたか）を、「TALYS」という理論計算ソフトで予測された値と比較しました。

• 理論にはいくつかの「レシピ（モデル）」がありますが、どれが正しいか試行錯誤しました。
• その結果、「アルファ粒子が原子核にどう近づくか（光学ポテンシャル）」を決めるレシピが、結果に最も大きな影響を与えることが分かりました。
• 実験データは、いくつかの理論モデルとよく一致しましたが、特に低いエネルギー領域では、理論が少し過小評価する傾向があることも分かりました。

5. 結論：なぜこれがすごいのか？

この研究の最大の功績は、「単に数値を測っただけ」ではなく、「その数値がどれくらい信頼できるか、そして誤差同士がどう関係しているか」まで、徹底的に解明した点にあります。

• 宇宙論への貢献： 宇宙で元素が作られるプロセスを、より正確にシミュレーションできるようになりました。
• 医療への貢献： がん診断に使われる新しい放射性同位体の製造方法を確立し、その安全性や効率を高めるデータを提供しました。

つまり、この論文は**「原子核という小さな世界で、正確無比な『ものさし』と『計算機』を使って、宇宙の秘密と未来の医療の鍵を解き明かした」**という物語なのです。

技術概要

以下は、提示された論文「Excitation function measurement of 144Sm(α,n) reaction at sub-Coulomb energies and detailed covariance analysis（サブクーロンエネルギー領域における 144Sm(α,n) 反応の励起関数測定および詳細な共分散解析）」の技術的サマリーです。

1. 研究の背景と課題 (Problem)

• 核天体物理学における重要性: 144Sm(α,n)147Gd 反応は、鉄より重い元素における p 過程（p-process）の理解、特に p 核（中性子不足核）の生成メカニズムを解明する上で極めて重要です。144Sm は魔法数を持つ p 核であり、太陽系における存在量が比較的多い（3.08%）ため、その反応断面積の正確な測定は天体物理モデルの精度向上に不可欠です。
• 医療応用: 反応生成物である 147Gd は、半減期が約 38 時間であり、229.32 keV の強いガンマ線（60.7%）を放出するため、単一光子放出断層撮影（SPET）用の潜在的な磁気医薬品として注目されています。
• 既存データの不足と不確実性: 従来の測定データは限られており、特にクーロン障壁（約 21.8 MeV）以下のエネルギー領域（14–21 MeV）における詳細な励起関数データが不足していました。また、従来の測定では、ビーム強度、標的厚さ、検出器効率など、複数の要因に由来する誤差の相関（共分散）を考慮した包括的な不確実性評価が十分に行われていませんでした。

2. 実験手法と方法論 (Methodology)

• 実験装置と条件:
  • 場所：Variable Energy Cyclotron Centre (VECC)、カルカッタ（インド）。
  • 手法：積層箔活性化法（Stacked foil activation technique）とオフラインガンマ線分光法。
  • ビーム：28 MeV のアルファ粒子ビームをアルミデグレーダー箔で減速し、標的箔に対して 14–21 MeV のエネルギー範囲で照射。
  • 標的：同位体濃縮（67%）された 144Sm2O3 粉末を、高純度アルミ箔上に分子堆積法で堆積（厚さ 280–350 μg/cm2）。
• エネルギー精度の向上（GEANT4 シミュレーション）:
  • ビームがデグレーダー箔や標的箔を通過する際のエネルギー損失とエネルギー広がり（エネルギー・ストラグル）を正確に評価するため、GEANT4 フレームワークを用いた詳細なモンテカルロシミュレーションを実施。
  • これにより、各標的位置における実効的な照射エネルギーとその不確実性（1σ）を定量化しました。
• ガンマ線測定と効率較正:
  • HPGe 検出器（40% 相対効率、1.33 MeV でのエネルギー分解能 1.8 keV）を用いて、147Gd の特徴的なガンマ線（229.3 keV）を測定。
  • 検出器効率の較正には 152Eu 標準線源を使用。
  • 重要: 検出器に近い配置（12.5 mm）での測定における「同時計数効果（coincidence summing）」と「有限の標的幾何形状」による補正を、EFFTRAN コードを用いて詳細に行い、効率曲線と共分散行列を算出。
• 共分散解析（Covariance Analysis）:
  • 測定された断面積の不確実性を、統計的誤差だけでなく、ビーム電流、標的厚さ、ガンマ線強度、検出器効率などの入力パラメータ間の相関を含めて評価。
  • ヤコビアン行列を用いて、入力パラメータの共分散から測定断面積の共分散行列および相関行列を構築しました。これは本反応において初めて行われた詳細な解析です。

3. 主要な貢献と成果 (Key Contributions & Results)

• 新たな励起関数データの取得:
  • クーロン障壁以下の 5 つのエネルギー点（14.09, 16.03, 17.68, 19.34, 20.90 MeV）において、144Sm(α,n)147Gd 反応の断面積を初めて高精度で測定しました。
  • 得られた断面積は、14 MeV で約 0.79 mb、21 MeV で約 595 mb であり、エネルギーとともに急激に増加する傾向を示しました。
• 包括的な不確実性評価:
  • 各測定点における相対的不確実性は約 25% でしたが、これには標的厚さの不確実性（15–25%）が主要な要因として寄与していました。
  • 測定値間の相関行列を初めて提示し、測定値同士が約 7–8% のレベルで相関していることを明らかにしました。これにより、反応速度などの物理量を推定する際の誤差評価の精度が向上します。
• 理論モデルとの比較（TALYS 2.0）:
  • 432 通りのパラメータ組み合わせ（8 種類のα光学ポテンシャル、6 種類の核レベル密度、9 種類のガンマ線強度関数）を用いた Hauser-Feshbach 統計モデル計算（TALYS 2.0）と比較しました。
  • 結果: 断面積の予測値は、主にα光学ポテンシャル（AOMP）の選択に敏感であることが判明しました。
    • AOMP3（Demetriou-Goriely ポテンシャル）と AOMP6（Avrigeanu らのポテンシャル）は、低エネルギー領域のデータをよく再現しましたが、高エネルギー側では若干過小評価する傾向がありました。
    • AOMP7 や AOMP8 は高エネルギー側を良く再現しますが、低エネルギー側で大きく過大評価しました。
  • 全パラメータの組み合わせによる理論予測の範囲（グレーバンド）は、実験データと概ね一致しましたが、特に低エネルギー領域での理論モデルの精度向上の必要性が示唆されました。

4. 意義と結論 (Significance & Conclusion)

• 天体物理モデルへの寄与: 本研究で得られたサブクーロンエネルギー領域の高精度な断面積データと、詳細な共分散情報は、p 過程における 144Sm の燃焼シミュレーションの精度を高め、宇宙における元素合成の理解を深める上で重要です。
• 医療用アイソトープ生産: 147Gd の生成断面積が明確化されたことで、医療用放射性同位体の生産効率の最適化や、必要な照射条件の設計に役立ちます。
• 手法論的革新: 本論文は、核反応断面積の測定において、単なる値の提示にとどまらず、測定値間の相関を含む「共分散行列」を体系的に評価・報告した点において画期的です。このアプローチは、将来の核データ評価（Nuclear Data Evaluation）や、反応速度論的な計算において、誤差の過小評価や過大評価を防ぐための重要な基盤となります。

総じて、本論文は実験的測定、高度なシミュレーション、および統計的誤差解析を統合し、144Sm(α,n) 反応に関する最も包括的なデータセットと評価を提供した重要な研究です。