Abstract:We demonstrate how using semi-stochastic Density Fitting (ss-DF) can accelerate self-consistent Density Matrix Embedding Theory (DMET) calculations by reducing the number of auxiliary orbitals in the three-indexed DF integrals. This results in time savings in building the Hartree Fock (HF) Coulomb and the Exchange Matrix and in integral transformation from the Atomic Orbital (AO) basis to the Embedding orbital (EO) basis. We then apply ss-DF to a range of hydrogen bonded clusters. First, we examine how the amount of deterministic space impacts the quality of the calculation in a $\left(\mathrm{H}_{2} \mathrm{O}\right)_{10}$ cluster. Next, we test the computational efficiency of ss-DF when compared to deterministic DF (d-DF) in water clusters containing 6 to 30 molecules of water using a triple-$\zeta$ basis set. Finally, we perform numerical structural optimisations on water and hydrogen fluoride clusters, revealing that DMET can recover weak interactions when using a back-transformed energy formula. This work therefore demonstrates the potential of using stochastic resolution of identity in quantum embedding theories and their potential to recover weak interactions.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to accelerate the calculations of the self - consistent Density Matrix Embedding Theory (DMET) by using the semi - stochastic Density Fitting (ss - DF) method. Specifically, the paper aims to save the time for constructing the Hartree - Fock (HF) Coulomb matrix and exchange matrix, as well as the integral transformation time from the Atomic Orbital (AO) basis to the Embedding Orbital (EO) basis by reducing the number of auxiliary orbits in the three - index integrals. In addition, the paper also explores the effect of applying the ss - DF method in hydrogen - bonded clusters, especially its performance in computational efficiency and recovering weak interactions. ### Main research contents: 1. **Method introduction**: - **Semi - stochastic Density Fitting (ss - DF)**: By setting a threshold, the auxiliary orbits are divided into deterministic and stochastic parts, thereby reducing the computational cost while maintaining the computational accuracy. - **Density Matrix Embedding Theory (DMET)**: A quantum embedding method used to deal with electron - correlation problems in large - molecule systems. 2. **Experimental design**: - **Testing of hydrogen - bonded clusters**: First, test the influence of different deterministic spaces on the calculation quality in a hydrogen - bonded cluster containing 10 water molecules. - **Comparison of computational efficiency**: In water clusters containing 6 to 30 water molecules, compare the computational efficiency of deterministic Density Fitting (d - DF) and semi - stochastic Density Fitting (ss - DF). - **Structure optimization**: Perform numerical structure optimization on water and hydrofluoric acid clusters to verify the performance of DMET in recovering weak interactions (such as hydrogen bonds). 3. **Result analysis**: - **Influence of threshold on accuracy**: As the random threshold increases, the computational precision and accuracy gradually decrease, especially when the threshold is greater than 0.2. - **Computational efficiency**: In large water clusters, ss - DF - DMET significantly reduces the computational time compared to d - DF - DMET, especially when the number of basis functions exceeds 1,000. - **Recovery of weak interactions**: Through structure optimization, the effectiveness of DMET in recovering weak interactions such as hydrogen bonds is verified, with a recovery rate of about 97%. ### Conclusion: The paper successfully shows that using semi - stochastic Density Fitting (ss - DF) can significantly accelerate DMET calculations and exhibits good computational efficiency and the ability to recover weak interactions when dealing with large hydrogen - bonded clusters. This method provides an effective solution for dealing with the electronic - structure problems of complex - molecule systems.

Accelerating Density Matrix Embedding with Stochastic Density Fitting Theory – An Application to Hydrogen Bonded Clusters

Accelerating Density Matrix Embedding with Stochastic Density Fitting Theory: An Application to Hydrogen Bonded Clusters

Density Fitted Fragment Embedding – Principles and Applications

Extending Density Functional Embedding Theory for Covalently Bonded Systems.

A Reduced Scaling Approach to Obtaining Embedding Orbital Coefficients: Benchmarking DMET-CCSD(T) Relative Energies in Water Clusters

Enhancing robustness and efficiency of density matrix embedding theory via semidefinite programming and local correlation potential fitting

Subspace Density Matrix Functional Embedding Theory: Theory, Implementation, and Applications to Molecular Systems

Frozen-Density Embedding for Including Environmental Effects in the Dirac-Kohn-Sham Theory: An Implementation Based on Density Fitting and Prototyping Techniques

Static Embedding with Pair Coupled Cluster Doubles Based Methods

Can Density Matrix Embedding Theory with the Complete Activate Space Self-Consistent Field Solver Describe Single and Double Bond Breaking in Molecular Systems?

Embedded density functional theory for covalently bonded and strongly interacting subsystems

Overlapped Embedded Fragment Stochastic Density Functional Theory for Covalently Bonded Materials

High-Throughput Condensed-Phase Hybrid Density Functional Theory for Large-Scale Finite-Gap Systems: The SeA Approach

How does HF-DFT achieve chemical accuracy for water clusters?

Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance

A Practical Guide to Density Matrix Embedding Theory in Quantum Chemistry

Density Functionals with Quantum Chemical Accuracy: From Machine Learning to Molecular Dynamics

Embedding theory contributions to average atom models for warm dense matter

Implementation of frozen density embedding in CP2K and OpenMolcas: CASSCF wavefunctions embedded in a Gaussian and plane wave DFT environment