[論文レビュー] User guide for the hfbcs-qrpa(v1) code

G. Colò, X. Roca-Maza|arXiv (Cornell University)|Feb 12, 2021

Nuclear physics research studies被引用数 4

ひとこと要約

hfbcs-qrpa(v1)コードは、スカイリムRPAフレームワークを拡張し、開殻球核におけるハートリー・フォック・BCS（HF-BCS）平均場計算と準粒子ランダム位相近似（QRPA）を組み込むことで、対相関と集団励起状態の自己無撞着な記述を可能にする。主な結果として、収束した対ギャップ（例：Δn = 1.45 MeV）、電荷半径（Rch = 4.64 fm）、遷移強度が得られ、エネルギー加重和則および数値的検証により検証された。

ABSTRACT

We briefly explain the structure of the Hartree-Fock-Bardeen-Cooper-Schrieffer (HFBCS) and Quasi-particle Random Phase Approximation (QRPA) code hfbcs-qrpa, highlighting only the differences with the code skyrme_rpa that has been previously published [G. Colò, L. Cao, N. Van Giai, and L. Capelli, Comp. Phys. Comm. 184, 142 (2013)]. The code deals with open shell spherical systems and non-spin flip, non-charge-exchange and natural parity $(-1)^l$ excitations.

研究の動機と目的

• スカイリムRPAコードを、開殻球核におけるHF-BCSとQRPAを用いた対相関を含めるように拡張すること。
• 電荷交換を伴わない自然パリティ遷移における準粒子励起状態と集団応答を計算する自己無撞着なフレームワークを提供すること。
• 密度依存性ゼロレンジ対力に基づき、手動または自動によるペアリング窓選択を柔軟に可能とすること。
• 電荷半径、対ギャップ、遷移確率、和則などの観測量を数値的収束チェックとともに計算すること。
• 出力ファイルとして強度関数、遷移確率、エネルギー加重和則を含む、使いやすくドキュメント化されたコードを提供すること。

提案手法

• 球対称性を仮定した座標空間で、スカイリム相互作用（例：SkM*、SLy4）とゼロレンジ密度依存性対力に基づくHF-BCS方程式を解く。
• 全エネルギーの2通りの評価方法（直接機能積分と準粒子エネルギーの総和に再配置項を加えたもの）を実装し、収束チェックを可能にする。
• 体積または表面型対力（パラメータV₀とx（x=0,1,0.5）で調整可能）を用い、自動モードでは120SnにおけるΔn=1.245 MeVを再現するようにV₀を固定する。
• 有限サイズの陽子および中性子電荷分布、異常磁気モーメント、スピン軌道寄与を用いて電荷半径を計算する。
• スカラーおよびアイソスカラー励起モードにおけるQRPA線形応答を計算し、遷移強度、和則（m(0), m(-1), m(1), m(3)）、重心エネルギーを算出する。
• 出力ファイルを生成：hfbcs-qrpa.out（遷移強度）、strength_is.outおよびstrength_iv.out（畳み込み強度関数）、可視化用にローレンツ型畳み込みを適用する。

実験結果

リサーチクエスチョン

• RQ1 HF-BCSフレームワークにおけるBCS対相関の組み込みが、開殻核における核の基底状態および対ギャップの記述に与える影響は何か？
• RQ2 体積型、表面型、混合型の異なる対力タイプが、計算された対ギャップおよび遷移強度に与える影響は何か？
• RQ3 計算されたエネルギー加重和則および重心エネルギーが、スカラーおよびアイソスカラーモードにおいて期待される値とどの程度一致するか？
• RQ4 コードの数値的実装が、関数的評価と準粒子エネルギー総和の両方のエネルギー評価方法間でどの程度一貫性を保っているか？
• RQ5 コードが120Snのような核において、実験的観測量（電荷半径、遷移確率など）をどの程度正確に再現できるか？

主な発見

• コードは例題系で中性子対ギャップΔn = 1.45 MeVを正確に計算し、実験的傾向と整合的である。
• 電荷半径はRch = 4.6418 fmとして計算され、陽子サイズ（R²p = 0.64 fm²）、中性子サイズ（R²n = -0.11 fm²）、スピン軌道寄与の寄与が含まれる。
• スカラー型エネルギー加重和則（EWSR）はm(1)(IS) = 2.0875×10⁶ MeVであり、97.5％の消耗度を示し、良好な収束性を示している。
• アイソスカラー型EWSRはm(1)(IV) = 2.3841×10⁶ MeVであり、D.C.極限はm(1)[D.C.]IS = 2.1407×10⁶ MeVであり、強い集団性を示している。
• スカラー型遷移の重心エネルギーはE_central(IS) = 7.9021 MeV、アイソスカラー型はE_central(IV) = 21.187 MeVであり、異なる集団性のパターンを反映している。
• 数値的チェックによりエネルギーの一貫性が確認され、|E_int/E_HF - 1| = 3.68×10⁻⁴と高く、収束精度が優れている。