MP2 energy and density for large molecular systems with internal error control using the Divide-Expand-Consolidate scheme.

Divide-Expand-Consolidate (DEC) is a local correlation method where the inherent locality of the electron correlation problem is used to express the correlated calculation on a large molecular system in terms of small independent fragment calculations employing small subsets of local HF orbitals. A crucial feature of the DEC scheme is that the sizes of the local orbital spaces are determined in a black box manner during the calculation. In this way it is ensured that the correlation energy has been determined to a predefined precision compared to a conventional calculation. In the present work we apply the DEC scheme to calculate the correlation energy as well as the electron density matrix for the insulin molecule using second order Møller-Plesset (MP2) theory. This is the first DEC calculation on a molecular system which is too large to be treated using a conventional MP2 implementation. The fragmentation errors for the insulin DEC calculation are carefully analyzed using internal consistency checks. It is demonstrated that size-intensive properties are determined to the same precision for small and large molecules. For example, the percentage of correlation energy recovered and the error per electron in the correlated density matrix depend only on the predefined precision and not on the molecular size.

[1]  R. Bartlett,et al.  Transferability in the natural linear-scaled coupled-cluster effective Hamiltonian approach: Applications to dynamic polarizabilities and dispersion coefficients. , 2008, The Journal of chemical physics.

[2]  H. Werner,et al.  The effect of local approximations in coupled-cluster wave functions on dipole moments and static dipole polarisabilitiesDedicated to Prof. W. Meyer on the occasion of his 65th birthday. , 2004 .

[3]  S. F. Boys.  Construction of Some Molecular Orbitals to Be Approximately Invariant for Changes from One Molecule to Another , 1960 .

[4]  R. Bartlett,et al.  The full CCSDT model for molecular electronic structure , 1987 .

[5]  Tatiana Korona,et al.  Transition strengths and first-order properties of excited states from local coupled cluster CC2 response theory with density fitting. , 2007, The Journal of chemical physics.

[6]  Benny G. Johnson,et al.  THE CONTINUOUS FAST MULTIPOLE METHOD , 1994 .

[7]  Branislav Jansík,et al.  Linear scaling coupled cluster method with correlation energy based error control. , 2010, The Journal of chemical physics.

[8]  Jing Ma,et al.  Linear scaling local correlation approach for solving the coupled cluster equations of large systems , 2002, J. Comput. Chem..

[9]  Georg Hetzer,et al.  Low-order scaling local electron correlation methods. I. Linear scaling local MP2 , 1999 .

[10]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[11]  T. Korona.  The effect of local approximations on first-order properties from expectation-value coupled cluster theory , 2011 .

[12]  Philippe Y. Ayala,et al.  Linear scaling coupled cluster and perturbation theories in the atomic orbital basis , 1999 .

[13]  Frederick R Manby,et al.  Analytical energy gradients for local second-order Møller-Plesset perturbation theory using density fitting approximations. , 2004, The Journal of chemical physics.

[14]  Christian Ochsenfeld,et al.  A reformulation of the coupled perturbed self-consistent field equations entirely within a local atomic orbital density matrix-based scheme , 1997 .

[15]  Trygve Helgaker,et al.  Implementation of the incremental scheme for one-electron first-order properties in coupled-cluster theory. , 2009, The Journal of chemical physics.

[16]  Kazuo Kitaura,et al.  Coupled-cluster theory based upon the fragment molecular-orbital method. , 2005, The Journal of chemical physics.

[17]  Weitao Yang,et al.  Insights into Current Limitations of Density Functional Theory , 2008, Science.

[18]  H. Nakai,et al.  How does it become possible to treat delocalized and/or open-shell systems in fragmentation-based linear-scaling electronic structure calculations? The case of the divide-and-conquer method. , 2012, Physical chemistry chemical physics : PCCP.

[19]  K. Kitaura,et al.  Analytic energy gradient for second-order Møller-Plesset perturbation theory based on the fragment molecular orbital method. , 2011, The Journal of chemical physics.

[20]  Hans-Joachim Werner,et al.  Comparison of explicitly correlated local coupled-cluster methods with various choices of virtual orbitals. , 2012, Physical chemistry chemical physics : PCCP.

[21]  Peter Pulay,et al.  Localizability of dynamic electron correlation , 1983 .

[22]  T. Daniel Crawford,et al.  Local correlation in coupled cluster calculations of molecular response properties , 2004 .

[23]  Hans-Joachim Werner,et al.  Local treatment of electron correlation in coupled cluster theory , 1996 .

[24]  S. Hirata,et al.  Fragment and localized orbital methods in electronic structure theory. , 2012, Physical chemistry chemical physics : PCCP.

[25]  Philippe Y. Ayala,et al.  Linear scaling second-order Moller–Plesset theory in the atomic orbital basis for large molecular systems , 1999 .

[26]  H. Werner,et al.  Local treatment of electron excitations in the EOM-CCSD method , 2003 .

[27]  Trygve Helgaker,et al.  Robust and Reliable Multilevel Minimization of the Kohn-Sham Energy. , 2009, Journal of chemical theory and computation.

[28]  Guntram Rauhut,et al.  Analytical energy gradients for local second-order Mo/ller–Plesset perturbation theory , 1998 .

[29]  Christian Ochsenfeld,et al.  Linear and sublinear scaling formation of Hartree-Fock-type exchange matrices , 1998 .

[30]  Christian Ochsenfeld,et al.  An atomic orbital-based reformulation of energy gradients in second-order Møller-Plesset perturbation theory. , 2008, The Journal of chemical physics.

[31]  Tatiana Korona,et al.  Local CC2 electronic excitation energies for large molecules with density fitting. , 2006, The Journal of chemical physics.

[32]  Hermann Stoll,et al.  The correlation energy of crystalline silicon , 1992 .

[33]  T. Crawford,et al.  Local correlation domains for coupled cluster theory: optical rotation and magnetic-field perturbations. , 2008, Physical chemistry chemical physics : PCCP.

[34]  Paul G. Mezey,et al.  A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions , 1989 .

[35]  Mikael P. Johansson,et al.  A stepwise atomic, valence-molecular, and full-molecular optimisation of the Hartree-Fock/Kohn-Sham energy. , 2009, Physical chemistry chemical physics : PCCP.

[36]  Michael Dolg,et al.  Fully automated implementation of the incremental scheme: application to CCSD energies for hydrocarbons and transition metal compounds. , 2007, The Journal of chemical physics.

[37]  Dimitrios G Liakos,et al.  Efficient and accurate approximations to the local coupled cluster singles doubles method using a truncated pair natural orbital basis. , 2009, The Journal of chemical physics.

[38]  J. Olsen,et al.  Coupled-cluster theory in a projected atomic orbital basis. , 2006, The Journal of chemical physics.

[39]  F. Neese,et al.  Efficient and accurate local approximations to coupled-electron pair approaches: An attempt to revive the pair natural orbital method. , 2009, The Journal of chemical physics.

[40]  Daniel Kats,et al.  Local CC2 response method for triplet states based on Laplace transform: excitation energies and first-order properties. , 2010, The Journal of chemical physics.

[41]  T. Crawford,et al.  Localized optimized orbitals, coupled cluster theory, and chiroptical response properties. , 2012, Physical chemistry chemical physics : PCCP.

[42]  M. Head‐Gordon,et al.  An accurate local model for triple substitutions in fourth order Møller–Plesset theory and in perturbative corrections to singles and doubles coupled cluster methods , 2000 .

[43]  M. Nooijen,et al.  Dynamically screened local correlation method using enveloping localized orbitals. , 2006, The Journal of chemical physics.

[44]  W. Pangborn,et al.  The structure of T6 human insulin at 1.0 A resolution. , 2003, Acta crystallographica. Section D, Biological crystallography.

[45]  Hans-Joachim Werner,et al.  An efficient local coupled cluster method for accurate thermochemistry of large systems. , 2011, The Journal of chemical physics.

[46]  Takeshi Ishikawa,et al.  Fragment molecular orbital calculation using the RI-MP2 method , 2009 .

[47]  Klaus Ruedenberg,et al.  Localized Atomic and Molecular Orbitals , 1963 .

[48]  T. Daniel Crawford,et al.  Locally correlated equation-of-motion coupled cluster theory for the excited states of large molecules , 2002 .

[49]  Stinne Høst,et al.  Local orbitals by minimizing powers of the orbital variance. , 2011, The Journal of chemical physics.

[50]  M. Urban,et al.  Optimized virtual orbitals for correlated calculations : Towards large scale CCSD(T) calculations of molecular dipole moments and polarizabilities , 2006 .

[51]  Marcin Ziółkowski,et al.  A Locality Analysis of the Divide-Expand-Consolidate Coupled Cluster Amplitude Equations. , 2011, Journal of chemical theory and computation.

[52]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[53]  Frederick R Manby,et al.  Tensor factorizations of local second-order Møller-Plesset theory. , 2010, The Journal of chemical physics.

[54]  P. Čárský,et al.  Efficient evaluation of Coulomb integrals in a mixed Gaussian and plane-wave basis using the density fitting and Cholesky decomposition. , 2012, The Journal of chemical physics.

[55]  Poul Jørgensen,et al.  The divide-expand-consolidate family of coupled cluster methods: numerical illustrations using second order Møller-Plesset perturbation theory. , 2012, The Journal of chemical physics.

[56]  Kazuo Kitaura,et al.  Exploring chemistry with the fragment molecular orbital method. , 2012, Physical chemistry chemical physics : PCCP.

[57]  Christian Ochsenfeld,et al.  Rigorous integral screening for electron correlation methods. , 2005, The Journal of chemical physics.

[58]  Paweł Sałek,et al.  Linear-scaling implementation of molecular response theory in self-consistent field electronic-structure theory. , 2007, The Journal of chemical physics.

[59]  R. Bartlett,et al.  A full coupled‐cluster singles and doubles model: The inclusion of disconnected triples , 1982 .

[60]  Wei Li,et al.  A refined cluster-in-molecule local correlation approach for predicting the relative energies of large systems. , 2012, Physical chemistry chemical physics : PCCP.

[61]  Trygve Helgaker,et al.  The augmented Roothaan-Hall method for optimizing Hartree-Fock and Kohn-Sham density matrices. , 2008, The Journal of chemical physics.

[62]  M. Head‐Gordon,et al.  A fifth-order perturbation comparison of electron correlation theories , 1989 .

[63]  Poul Jørgensen,et al.  Trust Region Minimization of Orbital Localization Functions. , 2012, Journal of chemical theory and computation.

[64]  Masato Kobayashi,et al.  Two‐level hierarchical parallelization of second‐order Møller–plesset perturbation calculations in divide‐and‐conquer method , 2011, J. Comput. Chem..

[65]  Wei Li,et al.  Improved design of orbital domains within the cluster-in-molecule local correlation framework: single-environment cluster-in-molecule ansatz and its application to local coupled-cluster approach with singles and doubles. , 2010, The journal of physical chemistry. A.

[66]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[67]  R. Fletcher.  Practical Methods of Optimization , 1988 .

[68]  Christian Ochsenfeld,et al.  Tighter multipole-based integral estimates and parallel implementation of linear-scaling AO-MP2 theory. , 2008, Physical chemistry chemical physics : PCCP.

[69]  Masato Kobayashi,et al.  Extension of linear-scaling divide-and-conquer-based correlation method to coupled cluster theory with singles and doubles excitations. , 2008, The Journal of chemical physics.

[70]  Guntram Rauhut,et al.  Analytical energy gradients for local coupled-cluster methods , 2001 .

[71]  Rodney J Bartlett,et al.  A natural linear scaling coupled-cluster method. , 2004, The Journal of chemical physics.

[72]  Thomas Kjærgaard,et al.  Molecular gradient for second-order Møller-Plesset perturbation theory using the divide-expand-consolidate (DEC) scheme. , 2012, The Journal of chemical physics.