Track synthesis, materials and catalysis as it is published.
Search papers, labs, and topics across Lattice.
The authors formulate a reduced-scaling implementation of the coupled-cluster-based static embedding method (MPCC) by integrating density fitting across all embedding layers, a Laplace-transformed Sylvester equation solver, and the distinguishable cluster approximation. This addresses the prohibitive computational scaling that typically prevents coupled-cluster precision from being applied to extended molecular systems where standard perturbation theory fails. The resulting approach brings the low-level solver down to $\mathcal{O}(N^4)$ and the screened interaction to a non-iterative $N_{\rm frag}\mathcal{O}(N^4)$ scaling, matching the accuracy of canonical MPCC across polyacetylene chains, ozone bond-breaking, and challenging S22 dimers.
Coupled-cluster accuracy no longer requires steep global scaling: combining density fitting with distinguishable cluster approximations slashes quantum embedding costs down to $\mathcal{O}(N^4)$ without sacrificing chemical accuracy.
An improved algorithm to implement our previously proposed static quantum embedding method MPCC (J. Chem. Phys. 161, 164107 (2024)) is presented. It uses the density-fitting (DF) approximation for the electron-repulsion integrals (ERIs) across all layers of the embedding method: the low-level method, the screened interaction, and the high-level method. Additionally, an approximate factorized solution of the Sylvester equation is employed based on the Laplace transformation for the low-level method. To build the screened interaction, we have explored an approximate construction of it based on the distinguishable cluster approximation (DCA) by Kats and Manby (J. Chem. Phys. 141, 061101 (2014)), which amounts to neglecting certain higher-order exchange and ladder diagrams. Combining all these techniques, we have obtained an $\mathcal{O}(N^4)$ algorithm for the low-level method; a non-iterative $N_{\rm frag} \mathcal{O}(N^4)$ algorithm for the screened interaction; and an iterative $\mathcal{O}(N_{\rm frag}^6)$ algorithm for the fragment solver, where $N_{\rm frag}$ is the number of fragment orbitals and $N$ is the total number of orbitals. The resulting low-scaling implementation of the MPCC method is applied to trans polyacetylenes, that is, C$_{2n}$H$_{2n+2}$ molecules for chain lengths $n=1-8$, and the timing and accuracy of the method are compared to the original unfactorized implementation. Moreover, the approximations are further validated by applying the method to the potential energy surface of the \ce{O3} molecule. Finally, we apply the method to the selected molecules in the S22 dataset for which the performance of the MP2 method is known to be poor. The results demonstrate that the DF-based implementation of the MPCC method provides a significant reduction in computational cost while maintaining accuracy comparable to that of the original implementation.