Using molecular similarity to construct accurate semiempirical electron structure theories

Benjamin G. Janesko,David Yaron
DOI: https://doi.org/10.1063/1.1785771
2004-09-15
Abstract:Ab initio electronic structure methods give accurate results for small systems, but do not scale well to large systems. Chemical insight tells us that molecular functional groups will behave approximately the same way in all molecules, large or small. This molecular similarity is exploited in semiempirical methods, which couple simple electronic structure theories with parameters for the transferable characteristics of functional groups. We propse that high-level calculations on small molecules provide a rich source of parametrization data. In principle, we can select a functional group, generate a large amount of ab initio data on the group in various small-molecule environments, and "mine" this data to build a sophisticated model for the group's behavior in large molecules. This work details such a model for electron correlation: a semiempirical, subsystem-based correlation functional that predicts a subsystem's two-electron density as a functional of its one-electron density. This model is demonstrated on two small systems: chains of linear, minimal-basis (H-H)5, treated as a sum of four overlapping (H-H)2 subsystems; and the aldehyde group of a set of HOC-R molecules. The results provide an initial demonstration of the feasibility of this approach.
Chemical Physics
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop an efficient and accurate semi - empirical electronic structure theory to overcome the computational complexity problem of ab - initio methods in dealing with large - molecule systems. Specifically, the author proposes a semi - empirical correlation function based on subsystems. By using high - precision ab - initio data of small molecules to parameterize the model, it is able to predict the electron - correlation effects in large - molecule systems. This method aims to combine the principles of molecular similarity and nearsightedness and utilize the transfer characteristics of functional groups in different molecular environments to construct a semi - empirical model capable of dealing with large - molecule systems. ### Main contributions of the paper 1. **Proposed a new semi - empirical model**: - This model is based on the electron - correlation function of subsystems and can predict the two - electron density matrix of subsystems as a function of its one - electron density matrix. - The model is parameterized by high - precision ab - initio data, ensuring the accuracy and reliability of the model. 2. **Verified the effectiveness of the model**: - The author verified the feasibility of the model through tests on two small systems. These systems include: - A linear minimum - basis - set hydrogen - atom chain (H - H)₅, which is decomposed into four overlapping (H - H)₂ subsystems. - The aldehyde groups in a set of HOC - R molecules. - The results show that this model performs well in predicting the electron - correlation energy of these systems and is superior to the traditional MP2 method. 3. **Explored the sources of error**: - The author analyzed in detail the three main approximate assumptions in the model: nearsightedness, principal - component dimension reduction, and the prediction from the one - electron density matrix to the two - electron density matrix. - By comparing different approximate methods, the author evaluated the influence of each assumption on the final result. 4. **Extended the application range of the model**: - The author discussed how to apply this model to density - functional theory (DFT), especially in combination with exact exchange and semi - empirical correlation functions. - Preliminary results show that this model can also provide good prediction effects in DFT calculations. ### Main formulas - **One - electron density matrix**: \[ 1D(a,b)=\langle\Phi|a^{\dagger}ab|\Phi\rangle \] - **Two - electron density matrix**: \[ 2D(ac,bd)=\frac{1}{2}\langle\Phi|a^{\dagger}a^{\dagger}cabd|\Phi\rangle \] - **Cumulative expansion of the two - electron density matrix**: \[ 2D(ac,bd)=\frac{1}{2}1D(a,b)1D(c,d)-\frac{1}{2}1D(a,d)1D(b,c)+2\Delta(ac,bd) \] - **Correlation energy**: \[ E_{\text{corr}}=\sum_{abcd}\langle ac|bd\rangle2\Delta(ac,bd) \] - **Prediction of the subsystem two - electron density matrix**: \[ 2\Delta[1D](ac,bd)=2\Delta_{\text{avg}}(ac,bd)+\sum_{j}2\Delta_{j}(ac,bd)\left(\alpha_{j}+\sum_{i}\left(\gamma_{ij}(1D|1D_{i})+\sigma_{ij}(1D|1D_{i})^{2}\right)\right) \] ### Conclusion This paper successfully solves the computational complexity problem of ab - initio methods in dealing with large - molecule systems by developing and verifying a new semi - empirical electron - correlation function. By using high - precision ab - initio data of small molecules, the author constructs a model that can accurately predict the electron - correlation effects in large - molecule systems. This method provides new tools and ideas for future research on large - molecule systems.