Nuclear $\gamma$-Ray Cascades as Markov Processes

原子の内部を、エネルギーが住まう、活気にあふれた多層階のスカイスクレイパー（超高層ビル）として想像してみてください。原子が励起されると（例えば、星の中での衝突によって）、高いフロアへと跳ね上がります。そこには永遠に留まることはできず、最も安定した状態である地上階へと戻りたがります。しかし、必ずしもエレベーターで直行するわけではありません。時には、ガンマ線と呼ばれる小さなエネルギーの跳躍による、曲がりくねった階段を一段ずつ降りていかなければならないこともあります。新しいフロアに着陸するたびに、次にどの階段を進むべきかを決めなければなりません。星々の混沌とした高エネルギーのキッチンの中で、これらの決定が、ある特別な長寿命元素がどれほど生成されるかを決定づけるのです。科学者たちは何十年もの間、これらの階段をマッピングしようとしてきましたが、地図はしばしば食い違い、その計算も非常に複雑です。なぜなら、すべてのステップが前のステップに依存しているからです。

この論文は、「マルコフ連鎖」という数学的ツールを用いて、この迷路を解くための巧妙な新しい方法を紹介しています。これは、システムが状態から状態へと移動する際の確率を追跡するための、洗練された手法です。著者は、励起された原子のエネルギー準位を、通過していく「一時的な（transient）」部屋として扱い、最終的な安定した基底状態を、一度入ると出られない「吸収的な（absorbing）」部屋として扱っています。この崩壊プロセス全体を巨大な確率ゲームへと変えることで、原子が基底状態に辿り着く確率と、短寿命の「アイソマー（同位体）」状態に辿り着く確率を正確に算出しています。鍵となる発見は、アルミニウム26の特定の重要なエネルギー準位（6398 keV）において、基底状態に辿り着く確率が 0.68 ± 0.06 (1σ) ± 0.13 (2σ) であるということです。この数値は極めて重要です。なぜなら、これが私たちの銀河系に存在する、実際に宇宙で光り輝いている「長寿命な種類」のアルミニウム26がどれくらいあるのかを教えてくれるからです。

原子迷路の物語

アルミニウム26の世界に飛び込んでみましょう。これは、宇宙の灯台として機能する放射性同位体です。それが崩壊するとき、非常に特定のガンマ線（高エネルギーの光線）を 1.809 MeV で放出します。天文学者たちは、この光を銀河中に発見しており、それは今まさに星の中で新しい元素が鍛えられている証拠となっています。しかし、どれほどのアルミニウム26が作られているのかを理解するには、その誕生に関する特定の詳細を知る必要があります。すなわち、マグネシウム25の原子核に陽子が衝突すると、それらが融合して、励起されたアルミニウム26が生成されるのです。この新しい原子は落ち着きがなく、落ち着こうと躍起になっています。

問題は、この落ち着きのない原子には、落ち着くための主に2つの方法があることです。一つは、非常に長い時間（約 7.2 × 10⁵ 年）生存することのできる「基底状態」まで一気に降りること。もう一つは、228 keV にある「捕らえられた」状態（アイソマー）に落ちることですが、これだとわずか 6.35 秒 で崩壊してしまいます。もし基底状態に辿り着けば、それは長い間生き残り、私たちの望遠鏡に見える存在となります。もし短寿命の状態に辿り着けば、すぐに消えてしまいます。天体物理学者が抱く大きな疑問は、「原子が長い道のりを選ぶ確率（これを f₀ と呼びましょう）はいくらなのか？」ということです。

長年、科学者たちは分岐比（ある経路を取る確率と別の経路を取る確率）を測定することで f₀ を計算しようとしてきました。しかし、データはバラバラでした。実験ごとに異なる数値が得られ、崩壊が連鎖（カスケード）として起こるため、一つのステップの測定における小さな誤差が、最終的な答えに非線形で混乱を招くような影響を与えてしまうのです。従来の数学は、これらの誤差が単純で独立していると仮定しがちですが、実際にはそうではありません。最初のステップの確率を間違えると、二番目のステップの予測も間違ってしまうため、誤差が積み重なっていくのです。

新しい「ゲームボード」のアプローチ

著者であるアタナシオス・プサルティス氏は、これらのステップを個別の推測として扱うのではなく、崩壊全体を一つの繋がったゲームとして扱うことにしました。彼は、確率論の概念である**吸収マルコフ連鎖（Absorbing Markov Chain）**を用いました。

スカイスクレイパーの高層階からスタートするボードゲームを想像してみてください。サイコロを振り、どの廊下を進むかを決めます。いくつかの廊下は他のフロアへと繋がり（一時的な状態）、いくつかの廊下は出口のドアへと繋がっています（吸収的な状態）。一度出口のドアを通り抜けると、ゲームは終了し、上に戻ることはできません。この原子のゲームにおいて、「出口のドア」は基底状態と短寿命のアイソマーです。「廊下」はガンマ線の遷移です。

この論文では、あらゆる可能な経路をマッピングする巨大な行列（数字のグリッド）を構築しています。

• Q は、「廊下」（励起状態）間の移動確率を表します。
• R は、「出口のドア」（基底状態またはアイソマー状態への着地）へ歩いて出る確率を表します。
• I は、出口のドアそのものを表します。

特定の数学的操作（行列の逆行列計算）を行うことで、著者は、経路がどれほど紆余曲折していようとも、特定の出口に辿り着く正確な確率を計算することができます。この方法は正確であり、曖昧な近似に頼ることはありません。

「ディリクレ」のダイスで不確実性を制御する

ここからが、この論文の真に巧妙な点です。実験データは完璧ではありません。あらゆる測定には「ゆとり（不確実性）」があります。もし各ステップに対して、そのゆとりの中でランダムに数値を選んでしまうと、誤って物理法則を破ってしまう可能性があります。例えば、「左に行く確率が60%、右に行く確率が60%」という経路を選んでしまうかもしれません。合計すると120%になってしまい、これは不可能です。確率は常に100%にならなければなりません。

これを解決するために、著者はディリクレ分布と呼ばれる特別な確率分布を使用しました。これは、魔法のように連動している「重み付きのダイス」のようなものです。もし一つの経路で高い数字が出たら、他のダイスは合計が常に100%になるように自動的に調整されます。これにより、シミュレーションされるすべての「ゲームの実行」が、物理法則を遵守するように設計されています。

著者は、6398 keV の特定のエネルギー準位（反応における 92 keV の共鳴エネルギーに対応）に対して、このシミュレーションを 5,000回 実行しました。彼は4つの主要な実験（Champagne et al., Kankainen et al., Lotay et al., Zhang et al.）からのデータを取得し、それぞれのデータセットに対してディリクレのダイスをマルコフ連鎖モデルの中で転がしました。

結果：より鮮明な姿

計算が終わったとき、この論文は 6398 keV 状態の基底状態への供給確率（f₀）について、非常に具体的な答えを導き出しました。4つの実験データを「階層的ベイズ法」（すべての実験を信頼しつつ、それらが隠れた相違を持っている可能性も認める高度な手法）を用いて組み合わせた結果は以下の通りです。

f₀ = 0.68 ± 0.06 (1σ) ± 0.13 (2σ)

これは、原子が基底状態に辿り着く確率が 68% であり、誤差の範囲が 0.06 であることを意味します。より確信を持ちたい場合（95%の信頼区間）には、その範囲は 0.13 まで広がります。

また、この論文はユニークなことも行いました。それは、階段のどの特定のステップが最も混乱の原因となっているかを突き止めたことです。数学的な分析を通じて、著者は 6398 keV から 2070 keV への遷移 が、不確実性の「ボトルネック」であることを特定しました。4つのデータセットのうち3つにおいて、この単一のジャンプが誤差の大部分（場合によっては 77% から 90% 以上）を占めていたのです。

なぜこれが重要なのか

これは単なる数学パズルの解決ではありません。f₀ の値は、星の中で生成されるアルミニウム26の量を計算する方法に直接影響を与えます。もしこの数値を間違えれば、銀河がどのように進化し、この元素がどれほど存在するのかという私たちのモデルが狂ってしまう可能性があるのです。

この論文は、これまでのより単純なカスケード計算の方法が、ステップ間の相関を正しく扱えなかったために的外れであったことを裏付けています。新しいマルコフ連鎖のアプローチはより高速で正確であり、科学者に明確な「今後のToDoリスト」を提示しています。もし、銀河に対する理解の不確実性を減らしたいのであれば、すべてのステップを測定するのではなく、6398 keV から 2070 keV へのあの特定のジャンプに最高の望遠鏡を集中させるべきである、というリストです。

要するに、著者は原子のスカイスクレイパーのためのより優れた地図を作り上げました。どこで階段が滑りやすく、どこで道が明快であるかを示し、私たちが宇宙の灯台をより高い精度で数えられるようにしてくれたのです。