[論文レビュー] 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

ひとこと要約

本論文では、多電子原子におけるLS結合スキームにおけるスピン・角運動量積分を効率的に行うためのFORTRAN 77ベースのソフトウェアツール、SAIライブラリを提示する。第二量子化、結合テンソル形式、および軌道、スピン、擬スピン空間における角運動量理論を活用する。二体演算子の行列要素の高速計算を可能とし、特に開いたf殻への対応を拡張し、多配置ハートリー・フォックおよび配置相互作用法における大規模な原子構造計算を向上させる。

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殻構造を含む角運動量手法の拡張。
• 多配置ハートリー・フォックおよび配置相互作用計算における角運動量積分の効率性とスケーラビリティの向上。
• 特に有効演算子または摂動的演算子の行列要素を必要とするコードへの統合を想定したモジュラーかつ再利用可能なコンponentの提供。
• LS結合における複雑な原子およびイオンの正確なエネルギー準位構造の計算、相対論的補正をα²階層まで含む。

提案手法

• 二体演算子の行列要素を、軌道空間およびスピン空間におけるランク(k_l, k_s)を持つ非可約テンソル演算子の形式で表すために、結合テンソル形式における第二量子化を採用。
• 軌道、スピン、擬スピンの三つの結合空間における角運動量理論を用い、電子相関および結合の体系的取り扱いを可能にする。
• 角運動量結合の図式的手法を適用し、縮約行列要素および再結合係数を効率的に計算。
• 標準的な量（縮約係数の分数的親和性、完全に縮約された行列要素など）の評価に、前処理用のSQライブラリと統合。
• 非直交軌道の問題を回避し、数値的安定性を保証するように設計された、角積分の再帰的アルゴリズムを実装。
• 開いたf^N殻を有する構造に対して最適化され、従来の実装よりも高速な実行が達成されている。

実験結果

リサーチクエスチョン

• RQ1 LS結合された多電子原子における二体演算子の角運動量積分を、大規模な原子構造計算においてどのように高速化できるか？
• RQ2 第二量子化および結合テンソル形式を用いて、開いたf殻を含むスピン・角運動量行列要素を効率的に計算する最良の方法は何か？
• RQ3 擬スピン形式と図式的手法の組み合わせにより、複雑な構造における角運動量結合の計算複雑度を低減できるか？
• RQ4 f殻系において、新しいライブラリは従来の実装と比較して、速度およびスケーラビリティの面でどの程度優れているか？
• RQ5 MCHFやATSP_MCHFのような既存の原子構造パッケージへのモジュラー統合は、どのように可能になるか？

主な発見

• SAIライブラリは、LS結合における二体演算子のスピン・角運動量行列要素を効率的に計算可能であり、従来手法と比較して実行速度が顕著に向上している。
• ライブラリは、従来の実装で取り扱いが困難であった開いたf殻（例：f^N構造）への角運動量積分能力を成功裏に拡張している。
• 結合テンソル形式および擬スピン理論の使用により、行列要素の体系的かつコンactな表現が可能となり、アルゴリズムの複雑度が低減されている。
• l = 0, 1, 2, 3およびl ≥ 3のl²に対応し、特にレアアースイオンを含む複雑な原子の正確な取り扱いが可能である。
• ベンチマーク結果では、さまざまな項（例：Ce⁺における4f^5 5d^1）について正確に行列要素が計算されており、既知のデータ（例：テストケースにおける²G状態のエネルギーが1.00000000Rとして計算）と一致している。
• MCHF_ASPおよびATSP_MCHFパッケージと互換性があり、完全な相関および相対論的補正を含む大規模な原子構造計算における能力を強化している。