[论文解读] The library of subroutines for calculation of matrix elements of two-particle operators for many-electron atoms

Gediminas Gaigalas|ArXiv.org|May 14, 2004

Atomic and Molecular Physics参考文献 15被引用 7

一句话总结

本论文介绍了 SAI 库，这是一个基于 FORTRAN 77 的软件工具，用于在 LS 耦合方案下对多电子原子实现高效的自旋-角动量积分，该工具利用了二次量子化、耦合张量形式以及轨道、自旋和赝自旋空间中的角动量理论。它能够加速双粒子算符矩阵元的计算，特别是扩展了对开放 f 亚层的支持，并提升了多组态 Hartree-Fock 和组态相互作用方法中的大规模原子结构计算效率。

ABSTRACT

In this paper a library for spin--angular integration in LS-coupling for many-electron atoms is presented. The software is an implementation of a methodology based on the second quantization in coupled tensorial form, the angular momentum theory in 3 spaces (orbital, spin and quasispin), and the graphical technique of angular momentum. This implementation extends applications of the relevant methodology to open f- shells and leads to faster execution of angular integration codes. The possibility of using some library routines for solving various angular momentum problems in atomic physics is also discussed.

研究动机与目标

• 开发一个高性能软件库，用于计算 LS 耦合多电子原子中双粒子算符的自旋-角动量矩阵元。
• 将现有的角动量方法扩展至包含此前计算困难的开放 f 亚层构型。
• 提高多组态 Hartree-Fock 和组态相互作用计算中角动量积分的效率与可扩展性。
• 提供一个模块化、可重用的组件，便于集成到原子结构代码中，特别是那些需要有效算符或微扰算符矩阵元的代码。
• 在 LS 耦合框架下，支持复杂原子和离子的精确能级谱计算，包括至 α² 阶的相对论修正。

提出的方法

• 采用耦合张量形式的二次量子化，将双粒子算符的矩阵元表示为轨道空间和自旋空间中秩为 (k_l, k_s) 的不可约张量算符。
• 在三个耦合空间（轨道、自旋和赝自旋）中应用角动量理论，实现对电子关联和耦合的系统处理。
• 应用角动量耦合的图解技术，高效计算约化矩阵元和再耦合系数。
• 与先前的 SQ 库集成，用于评估标准量，如约化分数父系数和完全约化矩阵元。
• 实现一种递归算法来处理角积分，避免非正交轨道带来的复杂性，并确保数值稳定性。
• 针对具有开放 f^N 亚层的构型进行优化，执行速度显著优于以往实现。

实验结果

研究问题

• RQ1 如何加速 LS 耦合多电子原子中双粒子算符的角动量积分，以支持大规模原子结构计算？
• RQ2 在二次量子化和耦合张量形式下，计算涉及开放 f 亚层的自旋-角动量矩阵元的最高效方法是什么？
• RQ3 结合赝自旋形式和图解技术，是否能降低复杂构型中角动量耦合的计算复杂度？
• RQ4 与以往实现相比，该新库在 f 亚层体系中的速度和可扩展性提升程度如何？
• RQ5 该库如何以模块化方式集成到现有原子结构软件包（如 MCHF 或 ATSP_MCHF）中，以实现更广泛的应用？

主要发现

• SAI 库实现了 LS 耦合下双粒子算符自旋-角动量矩阵元的高效计算，执行速度显著优于以往方法。
• 该库成功将角动量积分能力扩展至开放 f 亚层（例如 f^N 构型），这些构型在以往标准实现中难以处理。
• 采用耦合张量形式和赝自旋理论，使得矩阵元能够以系统化且紧凑的方式表示，降低了算法复杂度。
• 该库支持 l = 0, 1, 2, 3 及 l ≥ 3 时的 l²，能够准确处理包括稀土离子在内的复杂原子。
• 基准测试结果表明，该库能正确计算各种项（如 Ce⁺ 中的 4f^5 5d^1）的矩阵元，其结果与已知数据一致（例如，测试案例中 ²G 态能量计算为 1.00000000R）。
• 该库与 MCHF_ASP 和 ATSP_MCHF 软件包兼容，显著增强了其在包含全相关和相对论修正的大规模原子结构计算中的能力。