[论文解读] Electromagnetic Strength Distributions from the Ab Initio No-Core Shell Model

Christina Stumpf, Tobias Wolfgruber+1 more|arXiv (Cornell University)|Sep 20, 2017

Electromagnetic Scattering and Analysis参考文献 5被引用 3

一句话总结

本论文提出一种从头计算方法，结合无核壳模型（NCSM）与Lanczos强度函数方法，计算从16O到24O的氧同素体中电磁跃迁强度分布（E0、E1和E2）。采用规范两体和三体核子相互作用，实现从低激发态到巨共振的全能量范围覆盖，揭示了碎片化、精细结构及亚基态模式，同时暴露了RPA类近似方法的局限性，并为核相互作用提供了定量约束。

ABSTRACT

We present an ab initio approach for the description of collective excitations and transition strength distributions of arbitrary nuclei up into the sd-shell that based on the No-Core Shell Model in combination with the Lanczos strength-function method. Starting from two- and three-nucleon interactions from chiral effective field theory, we investigate the electric monopole, dipole, and quadrupole response of the even oxygen isotopes from 16-O to 24-O. The method describes the full energy range from low-lying excitations to the giant resonance region and beyond in a unified and consistent framework, including a complete description of fragmentation and fine-structure. This opens unique opportunities for understanding dynamic properties of nuclei from first principles and to further constrain nuclear interactions. We demonstrate the computational efficiency and the robust model-space convergence of our approach and compare to established approximate methods, such as the Random Phase Approximation, shedding new light on their deficiencies.

研究动机与目标

• 开发一种基于第一性原理的ab initio框架，用于计算sd壳区的电磁跃迁强度分布，涵盖低激发态和巨共振。
• 通过提供从第一性原理出发的完全一致且收敛的描述，克服RPA和SRPA等近似方法的局限性。
• 通过将计算得到的强度分布与实验数据对比，检验现代规范核相互作用（NN+3N）的预测能力。
• 研究中子富集的氧同素体中亚基态偶极和四极共振的出现机制。
• 评估基态关联和非局部性在决定超越平均场预期的共振能量中的作用。

提出的方法

• 采用无核壳模型（NCSM）结合规范两体和三体核子相互作用求解核多体问题，使用谐振子基矢截断至$N_{\text{max}}$。
• 应用重要性截断以控制基矢维度的阶乘增长并改善收敛性。
• 将Lanczos强度函数方法应用于基态本征态向量，以高效且高收敛性地计算跃迁强度分布。
• 该方法可直接计算全能量谱范围内的强度函数，包括碎片化和精细结构，而无需依赖连续区近似。
• 采用同位旋分解以识别贡献于E2强度的中子主导的低激发$2^+$态，表明其具有亚基态四极特征。
• 通过增加$N_{\text{max}}$的收敛性测试对方法进行验证，并与RPA和SRPA近似进行比较。

实验结果

研究问题

• RQ1 ab initio计算的氧同素体电磁强度分布与实验数据相比如何，特别是对于标量巨单极共振？
• RQ2 现代规范相互作用（NN+3N）在多大程度上再现了氧同素体链中E1和E2强度分布的观测系统性？
• RQ3 尽管相互作用再现了体积模量和基态半径，为何预测的E0、E1和E2模式的共振能量仍偏高？
• RQ4 低激发、中子主导的$2^+$态在E2强度分布中起什么作用，其与亚基态四极共振有何关联？
• RQ5 RPA和SRPA的缺陷在ab initio基准计算背景下如何体现？

主要发现

• NNLO$_{\text{sat}}$相互作用预测16O中标量E0巨单极共振的中心能量为19.8 MeV，高于实验值13.5 MeV，尽管其再现了核物质的体积模量。
• 同位旋对E1强度分布随20O向24O变化而展宽，低能强度增强，与亚基态偶极共振的出现一致。
• 在NNLO$_{\text{sat}}$相互作用下，偶极极化率从16O的0.48 fm³增加到24O的1.09 fm³，但所有值仍低于16O的实验值0.58 fm³。
• 18O及更重同素体中低激发$2^+$态对标量E2强度有显著贡献，表明存在亚基态四极共振，同位旋分解证实其具有强中子主导性。
• NN+3N(400)相互作用使基态半径低估约10%，导致其强度分布相对NNLO$_{\text{sat}}$发生偏移，后者被调制以再现实验半径。
• Lanczos方法在$N_{\text{max}}$增加时表现出稳健收敛性，实现了全能量范围内无需连续区近似的可靠强度函数计算。