[论文解读] Large-scale ab initio simulations based on systematically improvable atomic basis

Pengfei Li, Xiaohui Liu+6 more|arXiv (Cornell University)|Feb 28, 2015

Electron and X-Ray Spectroscopy Techniques被引用 4

一句话总结

本文提出 ABACUS，一种基于陈-郭-何方案导出的系统可改进数值原子基组的从头算 DFT 代码。该方法可实现对分子、固体、表面和缺陷的大规模、高精度模拟，通过优化的基组与线性标度算法实现高效的力计算。

ABSTRACT

We present a first-principles computer code package (ABACUS) that is based on density functional theory and numerical atomic basis sets. Theoretical foundations and numerical techniques used in the code are described, with focus on the accuracy and transferability of the hierarchical atomic basis sets as generated using a scheme proposed by Chen, Guo and He [J. Phys.:Condens. Matter extbf{22}, 445501 (2010)]. Benchmark results are presented for a variety of systems include molecules, solids, surfaces, and defects. All results show that the ABACUS package with its associated atomic basis sets is an efficient and reliable tool for simulating both small and large-scale materials.

研究动机与目标

• 开发一种可扩展、高精度且可靠的从头算 DFT 代码，用于大规模材料模拟。
• 实现一种系统可改进的原子基组生成方案，确保在各类材料中具备高可转移性与高精度。
• 通过数值原子轨道实现高效的力计算，包括费曼-海勒曼力、布卢伊力与非正交性力。
• 支持原子轨道与平面波基组，用于交叉验证与一致性检查。
• 集成先进泛函，如 LDA、GGA（PBE）以及用于范德华相互作用校正的 DFT-D2。

提出的方法

• 采用密度泛函理论（DFT）与 Kohn-Sham 形式，以数值原子轨道作为基函数。
• 利用陈-郭-何方案生成具有高可转移性的分层、系统可改进的原子基组。
• 采用极点展开与选择性逆运算（PEXSI）方法，实现大体系下低于立方标度（≤O(N²)）的计算效率。
• 通过双中心积分技术计算力的费曼-海勒曼力、布卢伊力与非正交性力贡献。
• 利用径向网格插值与实球谐函数计算原子轨道对原子坐标的导数。
• 支持统一赝势格式（UPF）中的赝势，并包含 LDA、GGA 与 DFT-D2 泛函。

实验结果

研究问题

• RQ1 基于系统可改进原子基组的从头算代码能否在各类材料中实现高精度与高可转移性？
• RQ2 此类代码在采用线性标度假设下，能否高效完成大规模 DFT 模拟？
• RQ3 在非正交原子轨道体系中，力计算的精度如何，包括布卢伊力与非正交性力的贡献？
• RQ4 双基组能力（原子轨道与平面波）能否确保基准模拟中的一致性与可靠性？
• RQ5 优化后的基组在表面、缺陷与扩展固体等复杂体系中的表现如何？

主要发现

• ABACUS 软件包成功实现了对小体系与大体系的从头算大规模模拟，展现出高效率与高可靠性。
• 对分子、固体、表面与缺陷的基准测试结果证实了陈-郭-何方案生成的基组具有高精度与高可转移性。
• PEXSI 方法的应用使 Kohn-Sham DFT 计算的标度低至 O(N^α)，其中 α ≤ 2，适用于绝缘体与金属体系。
• 力计算准确且高效，费曼-海勒曼力、布卢伊力与非正交性力通过优化的数值与解析技术实现。
• 双基组能力（原子轨道与平面波）支持交叉验证与一致性检查，显著增强了模拟结果的可信度。
• LDA、GGA（PBE）与 DFT-D2 泛函的实现可准确处理交换-关联效应，包括范德华相互作用。