Pair Transfer and Reaction Dynamics in $^{40,48}$Ca + $^{96}$Zr Collisions Below the Coulomb Barrier

技術要約：クーロン障壁以下における $^{40,48}\text{Ca} + ^{96}\text{Zr}$ 衝突の対転移および反応ダイナミクス

問題提起
サブバリア融合反応は、核構造と相関、特にペアリングの役割を調べるための敏感なプローブとなる。実験的観測は、二中性子転移断面積の増強（例：$^{40}\text{Ca} + ^{96}\text{Zr}$ および $^{60}\text{Ni} + ^{116}\text{Sn}$ において）を裏付けているが、時間依存フレームワーク内でBCS近似を用いた従来の理論的取り扱いは、実験的に観測された一中性子と二中性子の転移確率の比（$P_{2n}/P_{1n}^2$）の増強因子を再現することに失敗してきた。ペアリングを含まない理論モデルは0.25から0.5の間の比を予測するが、実験値はそれよりも大幅に大きい（$^{40}\text{Ca} + ^{96}\text{Zr}$ では約3）。さらに、ペアリング相関、動的な核の変形能、および転移に関与する特定の角運動量チャネルの間の相互作用は、微視的に未解明のままである。

手法
本研究では、融合閾値以下の $^{40,48}\text{Ca} + ^{96}\text{Zr}$ の正面衝突を調査するために、時間依存密度汎関数理論（TDDFT）を超流動系へと拡張した時間依存超流動局所密度近似（TDSLDA）を用いる。計算には、核エネルギー密度汎関数（NEDF）であるSeaLL1を使用する。ペアリングの効果を孤立させるため、著者らは参照状態として、同じコードと汎関数を用いた、ペアリングを含まない標準的な時間依存ハートリー・フォック（TDHF）理論を用いた並行シミュレーションを行う。

主要な技術的要素は以下の通りである：

• 射影演算子： 転移確率を微視的に抽出するために、反応後に中性子セクターに対して粒子数射影演算子を適用する。この手法は、SimenelによってTDHF用に提案され、超流動系へと一般化されたものであり、特定の数の核子が放出される断片の空間領域内に存在する状態である確率を決定する。
• 角運動量射影： 転移チャネルにおける中性子の $K$ 分布を特徴付けるために、角運動量 $K$ に対する追加の射影器が用いられ、逐次転移と相関した対転移を区別する。
• エネルギー基準の設定： 融合閾値はNEDFのパラメータ調整に対して敏感であるため（SeaLL1は実験値より約4 MeV高い閾値を予測する）、実験データとの一貫した比較を確実にするために、すべての観測量は計算された閾値エネルギーに対して記録される。

主要な結果

1. 動的変形能と平均転移： ペアリング相関は、原子核の動的な応答を著しく変化させる。TDSLDAでは、$^{96}\text{Zr}$ はクーロン場によって誘起された大きな分離距離において扁平（oblate）な四重極変形を発達させ、$\beta_{20} \approx -0.15$ から $-0.2$ で安定する。対照的に、TDHFは$^{96}\text{Zr}$ が接触まで球状を維持し、接触後に長球（prolate）変形を発達することを予測する。このTDHFにおける長球形状は、より大きな密度重なりと長い相互作用時間を可能にし、結果としてTDSLDAと比較してより大きな平均中性子転移をもたらす。これは、ペアリングが動的変形能（断熱性）を高める一方で、これらの特定の系においては平均転移の大きさを抑制することを示している。
2. 転移確率および増強比：
   • 一中性子転移： ペアリングが含まれる場合、平均転移の減少と一致して、一中性子転移の確率は小さくなる。
   • 二中性子転移： 二中性子転移の確率は、特に $^{48}\text{Ca} + ^{96}\text{Zr}$ 反応において、ペアリングによって著しく増強される。
   • 増強比 ($P_{2n}/P_{1n}^2$)： TDSLDA計算は、実験的に観測された増強因子（$^{40}\text{Ca} + ^{96}\text{Zr}$ では約3）を再現するが、従来のBCSベースのアプローチは、この比を過小評価していた。本研究は、$^{48}\text{Ca} + ^{96}\text{Zr}$ においてさらに大きな増強比を予測しており、これはその系における最終的なペアリングギャップの大きさに関連している。
3. $K$-角運動量分布： 角運動量射影の解析により、一中性子転移は主に $|K|=1/2$ の状態（衝突軸に沿って広がる波動関数に関連する）を伴う一方で、二中性子転移は $K=0$ チャネルによって支配されていることが明らかになった。この $K=0$ の支配は、相関したクーパー対の転移を意味しており、増強比に対する微視的な説明を与える。ここでの対の破壊メカニズムは、核分裂におけるものよりも、サブバリア衝突における励起エネルギーが低いため、それほど重要ではない。

意義と主張
本論文は、完全なTDSLDAフレームワークを用いて、低エネルギー重イオン衝突における一中性子および二中性子転移確率に関する初の微視的研究を提供すると主張している。主な意義は、従来の近似（TDHF+BCSなど）が正確に捉えることができなかった二中性子転移の実験的な増強因子を、TDSLDAが正常に再現できることを示した点にある。

著者らは、ペアリング相関は原子核の動的変形能と強く結びついており、それが平均転移を決定する一方で、対転移の増強は $K=0$ チャネルの占有によって駆動されると結論付けている。これらの結果は、TDSLDAが核の超流動性と反応ダイナミクスの相互作用を記述するための信頼できる微視的フレームワークであることを裏付けている。増強比 $P_{2n}/P_{1n}^2$ は、ペアリング相関の強度と理論的取り扱いの両方を制約するための重要な観測量として特定されている。