[论文解读] Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method

Julio César Cruz, Jorge Garza+2 more|arXiv (Cornell University)|Mar 10, 2022

Nuclear physics research studies参考文献 101被引用 2

一句话总结

本文提出了一种随机的、四组分的相对论性二阶多体微扰理论（MP2）方法，通过在实时间和虚时间中采用蒙特卡罗积分，绕过了计算成本高昂的四索引积分变换。该方法在系统尺寸上实现了潜在的二次方缩放（O(n²)），显著优于传统方法的O(n⁵)缩放，同时保持了高精度，并可在4096个处理器上实现92%的强可扩展性。

ABSTRACT

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-CDF (cumulative distribution function) method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with demonstrated 92% strong scalability going from 64 to 4096 processors for a fixed job size.

研究动机与目标

• 开发一种适用于强电子相关性和相对论效应的重元素化合物的可扩展电子相关方法。
• 通过消除四索引积分变换，克服传统四组分相对论性MP2方法中O(n⁵)的计算瓶颈。
• 证明随机积分可在空间紧凑的重系统中实现低于立方的缩放（潜在O(n²)），而这些系统中局部相关方法会失效。
• 利用高性能计算实现对含大量电子和基函数的系统进行大规模相对论性相关计算。

提出的方法

• 通过蒙特卡罗（MC）积分对涉及四组分旋量和库仑势的13维积分进行随机评估。
• 使用马尔可夫链蒙特卡罗（Metropolis）采样进行实空间电子坐标积分，使用逆累积分布函数（inverse-CDF）采样进行虚时间积分。
• 以随机采样替代确定性的四索引积分变换，避免了对大型积分的存储和变换。
• 利用独立蒙特卡罗采样固有的并行性，可在数千个处理器上高效分布。
• 采用重要性采样与权重函数，以降低统计误差并提高收敛性。
• 与狄拉克-哈特ree-福克参考方法集成，从一开始就包含相对论效应。

实验结果

研究问题

• RQ1 随机积分能否将四组分相对论性MP2的缩放从O(n⁵)降低至低于立方？
• RQ2 在空间紧凑的重分子中，该随机方法是否在显著降低计算成本的同时保持了高精度？
• RQ3 该方法能否在大规模并行架构（如数千个CPU或GPU）上实现高效缩放？
• RQ4 在不同系统尺寸和基组尺寸下，统计误差是否可控且可预测？
• RQ5 该方法是否适用于局部相关近似失效的系统，如镧系和锕系元素？

主要发现

• 计算成本随电子数或基函数数的增加，最坏情况为立方缩放，潜在可低至O(n²)，相比传统O(n⁵) MP2有显著改进。
• 在64到4096个处理器之间扩展时，该方法实现了92%的强可扩展性，表明并行效率极佳。
• 统计误差可控，且随采样量增加而减小，可实现高精度结果。
• 该随机方法成功避免了存储和变换四索引积分的需要，消除了相对论相关计算中的主要瓶颈。
• 该方法适用于空间紧凑的重系统，如锕系和镧系元素，而传统局部相关方法在这些系统中会失效。
• 数值证据表明，与确定性方法相比，该方法可在显著降低计算成本的同时实现高精度。