Efficient Implementation of NOCI-MP2 Using the Resolution of the Identity Approximation with Application to Charged Dimers and Long C-C Bonds in Ethane Derivatives.

Autor: Yost SR; Department of Chemistry and Biochemistry , Texas State University , San Marcos , Texas 78666 , United States.; Department of Chemistry , University of California , Berkeley , California 94720 , United States.; Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States., Head-Gordon M; Department of Chemistry , University of California , Berkeley , California 94720 , United States.; Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.
Jazyk: angličtina
Zdroj: Journal of chemical theory and computation [J Chem Theory Comput] 2018 Sep 11; Vol. 14 (9), pp. 4791-4805. Date of Electronic Publication: 2018 Aug 28.
DOI: 10.1021/acs.jctc.8b00697
Abstrakt: An efficient implementation of the perturb-then-diagonalize nonorthogonal configuration interaction method with second-order Møller-Plesset perturbation theory (NOCI-MP2) is presented. Relative to other low scaling multireference perturbation theories, NOCI-MP2 often requires a much smaller active space because of the use of nonorthogonal reference configurations. Reworking the NOCI-MP2 equations with the resolution of the identity (RI) approximation enables the method to have the same memory requirements and computational scaling as single reference RI-MP2. The working equations are extended to include single substitutions as required when the reference determinants do not satisfy the Hartree-Fock equations. A detailed computational algorithm is presented along with timings to establish the performance of the implementation. NOCI-MP2 is applied to the binding energy and charge resonance energy in dication and monocation π dimers, as well as didiamantane ethane, and hexaphenylethane. A well-defined set of nonorthogonal determinants are obtained using absolutely localized molecular orbitals (ALMOs), as solutions to the self-consistent field for molecular interactions (SCF-MI) equations corresponding to covalent and ionic determinants. Agreement with experimental information where available, and other multireference methods, is satisfactory, with the use of an 0.3 au level shift to guard against large MP2 amplitudes. For didiamantane ethane and hexaphenylethane, large dispersion forces help stabilize the molecules despite the steric repulsion. By contrast, in the case of hexaphenylethane, the energy penalty from the geometric distortion of the fragments significantly weakens the bond.
Databáze: MEDLINE