[论文解读] BSEPACK User's Guide

Meiyue Shao, Chao Yang|arXiv (Cornell University)|Dec 23, 2016

Zeolite Catalysis and Synthesis参考文献 3被引用 3

一句话总结

BSEPACK 是一个用于在分布式内存高性能计算（HPC）系统上求解贝特-萨佩特（Bethe–Salpeter）特征值问题的并行 ScaLAPACK 风格 Fortran 90 库，可高效计算材料科学中的电子吸收光谱。该库支持完整的 BSE 和塔姆-丹科夫近似（Tamm–Dancoff approximation, TDA），并提供多种求解器，包括直接对角化和迭代 Lanczos 方法，专为扩展体系的光学性质大规模从头计算而优化。

ABSTRACT

This is the user manual for the software package BSEPACK (Bethe--Salpeter Eigenvalue Solver Package).

研究动机与目标

• 开发一个可扩展的并行库，用于在大规模电子结构计算中求解贝特-萨佩特特征值问题。
• 实现通过完整 BSE 和塔姆-丹科夫近似（TDA）框架准确高效地计算材料吸收光谱。
• 支持多种求解策略，包括直接对角化和迭代 Lanczos 方法，以提升在高性能计算（HPC）平台上的性能。
• 提供一个可移植、模块化的 Fortran 90 库，依赖 MPI 和 ScaLAPACK，适用于集成到从头电子结构代码中。
• 通过全面测试和与 ScaLAPACK 兼容的标准二维块循环数据布局，确保正确性和性能。

提出的方法

• 该库以结构化形式实现贝特-萨佩特哈密顿矩阵：H = [[A, B], [-B, -A]]，其中 A 和 B 为复矩阵，且要求 H 为厄米特矩阵并正定。
• 采用 ScaLAPACK 约定的二维块循环方式在处理器网格上分布矩阵，以在分布式内存系统上实现高效的负载均衡和通信模式。
• 核心求解器 PDBSEIG（实数）和 PZBSEIG（复数）计算 H 的谱分解，返回特征值 Λ 和特征向量 X1、X2，通过特定酉条件进行归一化。
• 对于吸收光谱计算，该库使用带多种展宽函数（高斯型、洛伦兹型）的 Lanczos 算法，并支持直接和迭代求解路径。
• 与优化的 BLAS、LAPACK、ScaLAPACK 和 BLACS 库集成，提供多种求解器组合（如 FULLBSE+DIRECT、FULLBSE+LANCZOS）以实现灵活性和高性能。
• 库中包含示例程序和测试套件（如 runall.sh），用于验证正确性，输入数据以标准化文本格式组织，用于 A、B、d 和 σ。

实验结果

研究问题

• RQ1 如何在大规模分布式内存 HPC 系统上高效求解贝特-萨佩特特征值问题，同时最小化通信开销？
• RQ2 在完整 BSE 和 TDA 框架下，计算吸收光谱时，迭代求解器与直接求解器的最佳组合是什么？
• RQ3 如何有效利用二维块循环数据布局，在数千个处理器上扩展大规模贝特-萨佩特矩阵的求解？
• RQ4 在吸收光谱计算中，不同展宽方法（高斯型、洛伦兹型）的性能特征和数值稳定性如何？
• RQ5 如何设计一个模块化、可移植且可维护的库，以支持电子结构计算中的实数与复数运算？

主要发现

• 该库成功通过 TESTING 目录中全部 8192 项测试，测试通过率为 100%，表明实现正确且鲁棒。
• 该库支持多种求解器组合，包括通过 LAPACK HEEV/HEEVD/HEEVR/HEEVX 实现的直接对角化，以及使用高斯型和洛伦兹型展宽的迭代 Lanczos 方法。
• 通过 Lanczos 算法计算的吸收光谱在 ω 的 1024 个采样点上实现了准确结果，如示例程序 absorption_real.f 所示。
• 输入格式已标准化：矩阵 A、B 和 d 按行主序读取，实部与虚部分开；σ 作为单个实数读取。
• 该库设计用于与现有从头电子结构代码互操作，接口清晰，遵循 ScaLAPACK 约定，并支持实数与复数运算。
• 采用修改后的 BSD 许可证，允许开放使用和分发，学术使用需引用。