壳模型哈密顿量本征值的美与奇.PDFVIP

第30 卷第3 期 原子 核 物 理评论 Vol. 30, No. 3 2013 年9 月 Nuclear Physics Review Sep. , 2013 文章编号: 1007-4627(2013) 03-0260-08 壳模型哈密顿量本征值的美与奇 赵玉民 ( 上海交通大学物理与天文系，粒子物理与核物理研究所，上海 200240 ) 摘要: 本征值问题是自然科学中基本运算之一，对于超大矩阵的对角化是当今许多科学问题的瓶颈。 在应用原子核壳模型理论研究较重的原子核结构时，因为壳模型组态太大，通常的方法是基于各种物 理考虑做某些组态截断，另一个思路是利用新的算法和飞速发展的计算机资源对这些大矩阵对角化或 者近似对角化。总结了本课题组近年来在壳模型哈密顿量本征值近似方面研究的主要结果，包括最低 本征值半经验公式及多种外推方法、本征值与对角元的相关性等。 关键词: 壳模型；本征值；哈密顿量；组态截断 中图分类号: O571.2 文献标志码: A DOI: 10.11804/NuclPhysRev.30.03.260 [2] [3] 1 引言 论 和角动量投影方法 等。 另外一个思路是近似对角化，即得到近似的本征 现代科学许多问题与矩阵本征值有关，长期以来 值。如果近似程度足够好，人们也可以使用近似的 人们基于不同的条件发展了多种矩阵对角化的方法。 本征值作为理论计算结果，例如东京研究组的外推 [4] 例如Jacobi 算法。该算法的实质是把矩阵经过一系列 法 。这些外推法可以用现有的计算能力处理维数在 的正交变换，逐次消去绝对值最大的非零矩阵元。这 现有基础上再多几个数量级的矩阵。然而，其在实际 种方法适用于维数不太大( 比如小于3 000) 的实对称 计算中的可靠性还有待于更多的检验。 矩阵，给出所有的本征值和本征函数。当维数进一步 人们对本征值问题的兴趣起源于随机相互作用 [5] 增加时，循环次数增加，收敛的速度越来越慢。 的原子核结构理论工作 ，那里需要比较准确地估计 在原子核壳模型计算中，人们感兴趣的往往是比 给定自旋的壳模型哈密顿量的本征值。过去估计本 [6] 较低的能级，所以不需要把所有本征值全部求出来。 征值的方法，比较常用的是误差函数方法 以及通过 [1] ´ [7] 在这种情况下Lanczos 算法 有较大优越性。如果只 矩阵的迹和宽度等的经验方法。Velazquez 和Zuker 、