Large-scale shell-model investigation of $2\nu$ECEC in $^{132}$Ba and $^{78}$Kr

本論文は、$^{132}$Baおよび$^{78}$Krにおける二重電子捕獲の、検証済みの有効相互作用に基づいた更新された核行列要素と半減期予測を提供し、将来の実験的取り組みを支援するための、大規模な殻模型による調査を提示するものである。

要約

全体像：目に見えない幽霊を捕まえる

原子核を、小さくて混み合ったダンスフロアだと想像してみてください。通常、ダンサー（陽子と中性子）は非常に安定しており、パートナーを変えることはありません。しかし時折、2人のダンサーが全く同時に場所を入れ替える決断をすることがあります。これは二重電子捕獲と呼ばれる、非常に稀なイベントです。

この特定の「ダンス」では、原子核の中の2つの陽子が、原子の外殻から2つの電子を捕まえ、中性子へと変化します。この現象は極めて稀にしか起こらないため、一度発生するのを確認するには、信じられないほど長い時間（数兆年）がかかります。科学者たちが、この現象が正確にどれくらいの時間を要するか（半減期）を知りたい理由は、ニュートリノ（幽霊のように微小な粒子）の性質といった、宇宙の根本的なルールを理解する助けになるからです。

この論文の著者たちは、いわば建築家でありエンジニアです。彼らはこれらのイベントを捕まえるための新しい機械を作ったのではなく、代わりに、ダンスフロアがどのように振る舞い、どれくらい待つ必要があるのかを予測するための、超詳細なコンピュータ・シミュレーションを構築しました。

主役の二人：132Ba と 78Kr

研究者たちは、この稀なダンスの候補となる2つの特定の原子（原子核）に焦点を当てました。

1. バリウム132 (132Ba): このダンスを行うのではないかと科学者が疑っている重い原子ですが、まだ誰もその現場を捉えていません。彼らは古い地質学的な手がかりから、それが起こり得るということだけを知っています。
2. クリプトン78 (78Kr): 科学者がこのダンスが起きることを最近確認した原子ですが、測定値にはまだ少し曖昧さが残っています。

手法：「レゴ」シミュレーション

何が起こるかを予測するために、科学者たちはラージスケール・シェルモデルと呼ばれる手法を用いました。

• 例え話: 何十億個ものレゴブロックで作られた複雑な構造物が、どのように耐えられるかを予測しようとしている場面を想像してください。単に推測するのではなく、一つ一つのブロックが隣のブロックとどのように接続されているかを正確に知る必要があります。
• ツール: 科学者たちは、陽子と中性子がどのように相互作用するかをコンピュータに教える、巨大なデジタル「レゴセット」（有効相互作用と呼ばれます）を使用しました。
  • 132Baについては、SN100PNというセットを使用しました。
  • 78Krについては、GWBXGというセットを使用しました。

アップグレード: 彼らの以前のクリプトンに関する研究では、レゴの建物の「一階部分」しか見ていませんでした。今回の新しい研究では、モデルを拡張して「上の階」（高エネルギー準位）を含めました。これは、スカイスクレイパーが風でどのように揺れるかを理解するためには、基礎だけでなく、最上階まで見る必要があると気づいたようなものです。

設計図の検証：シミュレーションは機能したか？

この稀なダンスについての予測を信頼する前に、科学者たちはシлоベーションが正確であることを確認しなければなりませんでした。彼らは、関与する原子の「通常の」振る舞いをチェックすることでこれを行いました。

• エネルギー準位: コンピュータが原子の正しい「振動」やエネルギー状態を予測しているかどうかをチェックしました。
• 形状: 原子が球形なのか、あるいは少し潰れた卵のような形（変形）をしているのかをチェックしました。

結果: コンピュータ・シミュレーションは、実世界の実験データとほぼ完璧に一致しました。それは、橋の縮尺模型を作り、それが実際の橋と全く同じように耐えられることを確認するようなものでした。これにより、彼らの稀なダンスに関する予測も信頼できるものであるという自信が得られました。

主な知見：「待ち時間」の予測

1. バリウム132（謎の候補）

バリウム132がこのダンスを行うのをまだ誰も見ていないため、科学者たちは理論的な基準値を提示しました。

• 予測: もし約 7.33 × 10²⁴ 年（7の後に24個のゼロが続く数字です！）待てば、その現象を目撃できるかもしれないと計算しました。
• なぜ重要か: これは将来の実験のための「ターゲット」となります。これは科学者たちに対し、「100年待ってもダメだ。数兆年を待てるほどの検出器を作る必要がある」と伝えているのです。彼らの計算は、現在科学者が設定している最小限の限界よりもはるかに長く、探索はまだ継続中であることを意味しています。

2. クリプトン78（確認された候補）

科学者はすでにクリプトン78がこのダンスを行うのを見ていますが、測定値にはバラつきがあります。

• 予測: 新しく、より詳細なシミュレーションは、待ち時間が 8.78 × 10²² 年 であると予測しています。
• 改善点: 彼らの以前の研究（より小さなレゴセットを用いたもの）では、わずかに異なる時間を予測していました。モデルに「上の階」を加えたことで、新しい予測は最近の実験で実際に観測されている値に非常に近くなりました。これは、ぼやけた写真を高精細な画像にアップグレードしたようなもので、画像がより鮮明で正確になったのです。

「ボリュームノブ」（軸結合）

シミュレーションの難しい部分の一つは、コンピュータは宇宙のあらゆる微細な力をすべて知っているわけではないということです。これを修正するために、科学者たちは有効軸結合定数 ($g_{eff}^A$) と呼ばれる「ボリュームノブ」を使用します。

• 例え話: 曲を録音しているのですが、マイクが高音域を逃している場面を想像してください。そこで、マイクが逃した部分を補うために、ボリューム（ノブ）を上げます。
• 科学者たちは、待ち時間の予測がどのように変化するかを見るために、さまざまな「ボリューム設定」をテストしました。異なる設定を用いても、彼らの結果はこれまでに分かっていることと矛盾しませんでした。

結論：何を学んだのか？

この論文は次のように結論付けています。

1. シミュレーションは堅牢である: 彼らのコンピュータ・モデルは、これらの原子がどのように振る舞うかを記述する上で非常に優れています。
2. バリウム132: 彼らは、この崩壊を見るためにどれくらい待たなければならないかについて、現時点で最高の理論的推測を提供しました。これは、実験家たちが自分たちの検出器がどの程度の感度を持つ必要があるかを知る助けになります。
3. クリプトン78: より大きく複雑なモデルを見ることで、予測を改善し、実世界のデータにより一致させることができました。

要約すると、これらの科学者たちは、原子のダンスフロアのより優れた地図を作り上げたのです。彼らはまだダンサーを捕まえてはいませんが（バリウムの場合）、どこでいつ探すべきかについて、より優れた見通しを得ています。そしてクリプトンの場合、彼らの地図は以前のものよりもずっと正確になっています。

技術概要

技術要約：$^{132}$Baおよび$^{78}$Krにおける2νECECに関する大規模殻模型的研究

問題提起
二次の弱相互作用過程、特に二重電子捕獲（2νECEC）の研究は、基礎的なニュートリノの性質を探索し、核構造モデルを検証するために不可欠である。二重ベータ崩壊（2νββ）は確立されている一方で、2νECECモードは、通常、半減期が長くQ値が低いため、あまり探索されてこなかった。$^{132}$Baのような候補核（$^{132}$Baについては半減期の下限値のみが存在する）や、近年観測された$^{78}$Krのような核に対して、将来の実験的探索を導くためには、核行列要素（NME）および半減期の信頼できる理論的予測が極めて重要である。既存の理論的推定値は顕著なモデル依存性を示しており、これらの稀な崩壊モードを制約するために、独立して検証された堅牢な手法が必要とされている。

手法
著者らは、大規模核殻模型（NSM）計算を用いて、$^{132}$Baおよび$^{78}$Krにおける2νECEC過程を調査している。

• 有効相互作用： 本研究では、$^{132}$BaにはSN100PN有効相互作用を、$^{78}$KrにはGWBXG相互作用を用いている。
• モデル空間：
  • $^{132}$Ba（親核）の場合、モデル空間は、両方のプロトンおよび中性子に対して$0g_{7/2}1d_{2s}0h_{11/2}$軌道を含み、$^{100}$Snを不活性コアとする。計算の実現可能性を確保するために、構成を特定の粒子・空孔限界に制限する打ち切りが適用された。
  • $^{78}$Krについては、中性子のモデル空間を拡張することで、先行研究を拡張している。先行研究は$N=50$殻閉鎖に限定されていたが、本研究では、$N=50$を越えて$0g_{7/2}$、$1d_{5/2}$、$1d_{3/2}$、$2s_{1/2}$軌道への一粒子・一空孔（1p-1h）励起を含めている。
• 計算ツール： 一体型遷移密度（OBTD）にはNUSHELLXコードを使用し、エネルギースペクトルおよび遷移確率はKSHELLコードを用いて計算した。
• 検証： 相互作用の信頼性は、まず、親核、中間核、および娘核の分光学的性質（低励起エネルギースペクトルおよび減少四重極遷移確率 $B(E2)$）を実験データと比較することによって検証された。
• NME計算： ガモフ・テラー（GT）NMEは、標準的な一体型GT演算子を用いて計算された。中間核（$^{132}$Csおよび$^{78}$Br）における$1^+$状態のエネルギーに対するNMEの累積的な寄与を分析し、収束性を確認した。多体相関や二体電流による欠損を補うために、有効軸性ベクトル結合定数（$g^A_{eff}$）が適用された。

主な貢献

1. $^{132}$Baにおける初の殻模型NME： 本研究は、$^{132}$Baにおける2νECEC過程のNMEに関する初の殻模型計算を提示している。
2. $^{78}$Krの改良されたモデル空間： 本研究は、中性子のモデル空間を$N=50$殻閉鎖を超えて拡張することにより、構成制限のある従来の研究の限界に対処し、$^{78}$Krの更新された改良された殻模型推定値を提供している。
3. 分光学的検証： 著者らは、採用された有効相互作用が、崩壊鎖に関与する核の分光学的性質（エネルギー準位および$B(E2)$強度）を実験データと良好に再現することを示し、NME計算のための信頼できる基盤を確立した。
4. 感度分析： 中間核$^{132}$Csにおける未確定の最低$1^+$状態に関する$^{132}$Baの計算について、感度分析が行われた。

結果

• 分光学： 計算結果は、$^{132}$Ba、$^{132}$Cs、$^{132}$Xe、$^{78}$Kr、$^{78}$Br、および$^{78}$Seの実験データと極めて良好な一致を示した。$^{132}$Baについては、結果は中程度のプロレート変形（$\beta_2 \approx +0.18$）を示唆しており、一方で$^{78}$Krは形状混合を伴うオブレート（扁平）な固有形状を示している。
• $^{132}$BaのNMEおよび半減期：
  • 累積NMEは、中間核における6.3 MeV付近で0.0654の値に飽和する。
  • $g^A_{eff} = 0.94$ を用いた場合、予測される半減期は $7.33 \times 10^{24}$ 年である。これは現在の実験的な下限値（$> 2.2 \times 10^{21}$ 年）よりも約3桁長い。
  • この計算は、$^{132}$Csの最低$1^+$状態が100 keV以下であれば堅牢であるが、この状態が500 keVを超えるとNMEの顕著な抑制が起こる。
• $^{78}$KrのNMEおよび半減期：
  • 累積NMEは、飽和後に最終的に 0.2311 に達し、これは制限されたモデル空間を用いた従来の数値である0.1997からの改善である。
  • $g^A_{eff} = 0.8$ を用いた場合、予測される半減期は $8.78 \times 10^{22}$ 年である。
  • この結果は、最近測定された実験値 $[1.9^{+1.3}_{-0.7} \pm 0.3] \times 10^{22}$ 年と合理的な一致を見せている。
  • また、本研究では、競合する陽電子放出チャネル（$2\nu\beta^+\beta^+$ および $2\nu EC\beta^+$）の半減期をそれぞれ $4.68 \times 10^{26}$ 年および $1.19 \times 10^{23}$ 年と推定している。

意義および主張
著者らは、この研究が、この稀な崩壊モードを制約するための将来の実験的努力を支援するための $^{132}$Baのベースラインとなる理論的予測 を提供すると主張している。$^{78}$Krについては、より制限されたモデル空間を用いた以前の研究と比較して、更新され改良された殻模型の推定値 を提供している。結果は、有効軸性結合の使用が、標準的な一体型演算子には明示的に含まれていない多体効果（二体電流など）を部分的に吸収していることを認めつつも、典型的な理論的不確かさの範囲内で妥当であると提示されている。論文は、分光学的合致が波動関数を検証する一方で、崩落演算子は$g^A_{eff}$による繰り込みを必要とすることを強調している。著者らは、彼らの知見が2νECEC NMEを評価するための基礎となる核波動関数の信頼性を支持していると結論付けているが、より大きなモデル空間や改良されたアイソスピン対称性相互作用を用いることで、これらの予測をさらに精緻化できる可能性があるとも述べている。