[论文解读] Variational approach for the electronic structure calculation on the second-order reduced density matrices and the $N$-representability problem

Maho Nakata, Mituhiro Fukuda+1 more|arXiv (Cornell University)|Oct 20, 2010

Advanced Chemical Physics Studies参考文献 54被引用 5

一句话总结

本文提出了一种基于二阶约化密度矩阵（2-RDM）的变分电子结构计算方法，通过将问题表述为半定规划（SDP），避免了传统方法的指数级缩放问题。该方法通过施加近似N-可表示性条件（P、Q、G、T1、T2′），实现了与当前量子化学领域金标准CCSD(T)方法相当的精度，对于O₂⁺等具有挑战性的开壳层体系，相关能可达全配置相互作用（full CI）的132%。

ABSTRACT

The reduced-density-matrix method is an promising candidate for the next generation electronic structure calculation method; it is equivalent to solve the Schrödinger equation for the ground state. The number of variables is the same as a four electron system and constant regardless of the electrons in the system. Thus many researchers have been dreaming of a much simpler method for quantum mechanics. In this chapter, we give a overview of the reduced-density matrix method; details of the theories, methods, history, and some new computational results. Typically, the results are comparable to the CCSD(T) which is a sophisticated traditional approach in quantum chemistry.

研究动机与目标

• 开发一种基于二阶约化密度矩阵（2-RDM）的变分电子结构方法，避免全配置相互作用方法的指数级缩放问题。
• 通过应用近似必要条件（P、Q、G、T1、T2′）解决2-RDM理论中的N-可表示性问题，以确保物理一致性。
• 证明2-RDM方法能够实现与CCSD(T)相当的相关能，后者是当前量子化学领域的金标准。
• 利用超级计算机和SDPARA半定规划求解器，求解大规模2-RDM问题，突破现有计算可行性的极限。

提出的方法

• 将电子结构问题表述为在2-RDM空间上的能量泛函变分最小化问题，受物理约束条件限制。
• 以2-RDM作为基本变量，使变量数量恒定（四个），与体系大小无关，优于波函数基方法。
• 施加必要的N-可表示性条件——P、Q、G、T1和T2′，以确保2-RDM对应于有效的N电子密度矩阵。
• 将问题重新表述为半定规划（SDP），在正半定矩阵上最小化线性泛函。
• 在超级计算机上使用SDPARA求解器处理大规模2-RDM问题，包括NH₃、C₂、CH₃和O₂⁺等体系。
• 通过与57个原子和分子体系的全配置相互作用（full CI）和CCSD(T)基准值比较，验证了结果的准确性。

实验结果

研究问题

• RQ1 变分2-RDM方法是否能在多种原子和分子体系中实现与CCSD(T)相当的精度？
• RQ2 近似N-可表示性条件（P、Q、G、T1、T2′）在恢复电子相关能方面的有效性如何？
• RQ3 该2-RDM方法在O₂⁺等具有挑战性的开壳层体系和强关联体系中的表现如何？
• RQ4 在现代SDP求解器和超级计算资源支持下，2-RDM方法在多大程度上可扩展至大体系？

主要发现

• 采用P、Q、G、T1和T2′条件的2-RDM方法，其基态能量与当前量子化学领域金标准CCSD(T)相当。
• 对于NH₂⁻、CH₂、NH₃和CH₃，该方法分别恢复了全配置相互作用相关能的100.4%、100.4%、100.4%和100.3%。
• 对于C₂，该方法实现了全配置相互作用相关能的101.2%，表明在多参考区域也具有高精度。
• 对于O₂⁺，由于计算限制仅使用P、Q和G条件，该方法仍达到了全配置相互作用相关能的132%，凸显其在开壳层体系中的鲁棒性。
• 该方法显著优于HF方法，57个测试体系中与全配置相互作用的能差始终低于HF方法。
• 迄今为止求解的最大规模2-RDM问题包括NH₃、C₂、CH₃和O₂⁺，均使用超级计算机上的SDPARA求解器求解，证明了其可扩展性。