摘要
本文提出了一种用于原子大规模相关计算的算法，该算法显著降低了随着构建N粒子态所需单粒子函数数量增加而带来的计算复杂度上升。这种优化体现在将哈密顿矩阵与基态的系数或振幅在诸如戴维森方法之类的迭代过程中进行运算的阶段。该算法基于在网格上表示径向单粒子函数，并利用电子-电子相互作用的多极展开，通过对旋量、系数和势能进行一系列变换来实现。此外，它还通过原型法处理角积分，并分别计算构型态函数的相关系数，而非针对成对的构型态函数进行处理。通过双电子原子的计算验证了该算法的效率，并从所需操作次数方面分析了该算法与传统方法在各个阶段的差异。
