Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary.
ELKH-BME Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary.
J Chem Theory Comput. 2023 May 23;19(10):2850-2862. doi: 10.1021/acs.jctc.3c00101. Epub 2023 May 3.
Our reduced-cost scheme [ , , 094111] based on the frozen virtual natural orbital and natural auxiliary function approaches is extended to core excitations. The efficiency of the approximation is presented for the second-order algebraic-diagrammatic construction [ADC(2)] method invoking the core-valence separation (CVS) and density fitting approaches. The errors introduced by the present scheme are comprehensively analyzed for more than 200 excitation energies and 80 oscillator strengths, including C, N, and O K-edge excitations, as well as 1 → π* and Rydberg transitions. Our results show that significant savings can be gained in computational requirements at the expense of a moderate error. That is, the mean absolute error for the excitation energies, being lower than 0.20 eV, is an order of magnitude smaller than the intrinsic error of CVS-ADC(2), while the mean relative error for the oscillator strengths is between 0.06 and 0.08, which is still acceptable. As significant differences for different types of excitations cannot be observed, the robustness of the approximation is also demonstrated. The improvements in the computational requirements are measured for extended molecules. In this case, an overall 7-fold speedup is obtained in the wall-clock times, while dramatic reductions in the memory requirements are also achieved. In addition, it is also proved that the new approach enables us to perform CVS-ADC(2) calculations within reasonable runtime for systems of 100 atoms using reliable basis sets.
我们基于冻结虚拟自然轨道和自然辅助函数方法的低成本方案[,, 094111]被扩展到核心激发。我们提出了核心价层分离(CVS)和密度拟合方法的二阶代数图论构造[ADC(2)]方法中该近似的效率。通过对超过 200 个激发能和 80 个振子强度的分析,包括 C、N 和 O K 边激发以及 1→π*和里德堡跃迁,全面分析了本方案引入的误差。我们的结果表明,以适度的误差为代价,可以在计算要求方面节省大量的资源。也就是说,激发能的平均绝对误差低于 0.20 eV,比 CVS-ADC(2)的固有误差小一个数量级,而振子强度的平均相对误差在 0.06 到 0.08 之间,仍然是可以接受的。由于不同类型的激发没有明显的差异,因此也证明了该近似的稳健性。在扩展分子中测量了计算要求的改进。在这种情况下,在壁钟时间上获得了 7 倍的总体加速,同时也显著减少了内存需求。此外,还证明了新方法使我们能够在合理的运行时间内使用可靠的基组对 100 个原子的系统进行 CVS-ADC(2)计算。