On the nature of electron correlation in C60

The ground state restricted Hartree Fock (RHF) wave function of C(60) is found to be unstable with respect to spin symmetry breaking, and further minimization leads to a significantly spin contaminated unrestricted Hartree Fock (UHF) solution ( = 7.5, 9.6 for singlet and triplet, respectively). The nature of the symmetry breaking in C(60) relative to the radicaloid fullerene, C(36), is assessed by energy lowering of the UHF solution, , and the unpaired electron number. We conclude that the high value of each of these measures in C(60) is not attributable to strong correlation behavior as is the case for C(36). Instead, their origin is from the collective effect of relatively weak, global correlations present in the π space of both fullerenes. Second order perturbation (MP2) calculations of the singlet triplet gap are significantly more accurate with RHF orbitals than UHF orbitals, while orbital optimized opposite spin second order correlation (O2) performs even better.     

Scaled second-order perturbation corrections to configuration interaction singles: efficient and reliable excitation energy methods

Two modifications of the perturbative doubles correction to configuration interaction with single substitutions (CIS(D)) are suggested, which are excited state analogues of ground state scaled second-order M?ller-Plesset (MP2) methods. The first approach employs two parameters to scale the two spin components of the direct term of CIS(D), starting from the two-parameter spin-component scaled (SCS) MP2 ground state, and is termed SCS-CIS(D). An efficient resolution-of-the-identity (RI) implementation of this approach is described. The second approach employs a single parameter to scale only the opposite-spin direct term of CIS(D), starting from the one-parameter scaled opposite-spin (SOS) MP2 ground state, and is called SOS-CIS(D). By utilizing auxiliary basis expansions and a Laplace transform, a fourth-order algorithm for SOS-CIS(D) is described and implemented. The parameters that describe SCS-CIS(D) and SOS-CIS(D) are optimized based on a training set that includes valence excitations of various organic molecules and Rydberg transitions of water and ammonia, and they significantly improve upon CIS(D) itself. The accuracy of the two methods is found to be comparable. This arises from a strong correlation between the same-spin and the opposite-spin portions of the excitation energy terms. The methods are successfully applied to the zincbacteriochlorin-bacteriochlorin charge-transfer transition, for which time-dependent density functional theory, with presently available exchange-correlation functionals, is known to fail. The methods are also successfully applied to describe various electronic transitions outside of the training set. The efficiency of the SOS-CIS(D) and the auxiliary basis implementation of CIS(D) and SCS-CIS(D) are confirmed with a series of timing tests.     

Symmetry breaking in benzene and larger aromatic molecules within generalized valence bond coupled cluster methods

The origin of symmetry breaking (SB) in benzene in generalized valence bond methods is investigated within a coupled cluster formalism that correlates all valence electrons. Retention of a limited number of pair correlation amplitudes (as in the perfect- and imperfect-pairing models) that incompletely describes interpair correlations leads to symmetry breaking as the orbitals and amplitudes are optimized. Local correlation models that are exact for one, two, and three interacting pairs at the doubles excitation level are compared against the exact pair correlation treatment, which correlates four interacting pairs at once in the connected double substitution operator. For simplicity, this comparison is performed with a second-order model of electron correlation, which is reasonably faithful to the infinite-order result. The significant SB known for the one-pair model (perfect pairing) is not eliminated at the two-pair level, but is virtually eliminated at the three-pair level. Therefore, a tractable hybrid model is proposed, which combines three-pair correlations at the second-order level and infinite-order treatment for the strong imperfect-pairing correlations involving one and two-pair correlations. This model greatly reduces SB in benzene and larger delocalized pi systems such as naphthalene and the phenalenyl cation and anion. The resulting optimized orbitals are localized in the sigma space but exhibit significant delocalization in the pi space. This means that correlation effects associated with different resonance structures are treated in a more balanced way than if the pi orbitals localize, leading to reduced SB.     

Improvement of the coupled-cluster singles and doubles method via scaling same- and opposite-spin components of the double excitation correlation energy

There has been much interest in cost-free improvements to second-order M?ller-Plesset perturbation theory (MP2) via scaling the same- and opposite-spin components of the correlation energy (spin-component scaled MP2). By scaling the same- and opposite-spin components of the double excitation correlation energy from the coupled-cluster of single and double excitations (CCSD) method, similar improvements can be achieved. Optimized for a set of 48 reaction energies, scaling factors were determined to be 1.13 and 1.27 for the same- and opposite-spin components, respectively. Preliminary results suggest that the spin-component scaled CCSD (SCS-CCSD) method will outperform all MP2 type methods considered for describing intermolecular interactions. Potential energy curves computed with the SCS-CCSD method for the sandwich benzene dimer and methane dimer reproduce the benchmark CCSD(T) potential curves with errors of only a few hundredths of 1 kcal mol(-1) for the minima. The performance of the SCS-CCSD method suggests that it is a reliable, lower cost alternative to the CCSD(T) method.     

Scaled opposite spin second order Møller-Plesset theory with improved physical description of long-range dispersion interactions

Separate scaling of the same-spin and opposite spin contributions to the second-order M?ller-Plesset energy can yield statistically improved performance for a variety of chemical problems. If only the opposite spin contribution is scaled, it is also possible to reduce the computational complexity from fifth order to fourth order in system size, with very little degradation of the results. However neither of these scaled MP2 energies recovers the full MP2 result for the dispersion energy of nonoverlapping systems. This deficiency is addressed in this work by using a distance-dependent scaling of the opposite spin correlation energy. The resulting method is compared against the previously proposed scaled MP2 methods on a range of problems involving both short and long-range interactions.     

Scaled opposite-spin second order Møller-Plesset correlation energy: an economical electronic structure method

A simplified approach to treating the electron correlation energy is suggested in which only the alpha-beta component of the second order M?ller-Plesset energy is evaluated, and then scaled by an empirical factor which is suggested to be 1.3. This scaled opposite-spin second order energy (SOS-MP2), where MP2 is M?ller-Plesset theory, yields results for relative energies and derivative properties that are statistically improved over the conventional MP2 method. Furthermore, the SOS-MP2 energy can be evaluated without the fifth order computational steps associated with MP2 theory, even without exploiting any spatial locality. A fourth order algorithm is given for evaluating the opposite spin MP2 energy using auxiliary basis expansions, and a Laplace approach, and timing comparisons are given.     

Assessing spin-component-scaled second-order M?ller-plesset theory using anharmonic frequencies

Four common parametrisations of spin-component-scaled second-order M?ller-Plesset (MP2) theory are benchmarked by calculating the anharmonic vibrational frequencies of a test suite consisting of eighteen diatomic and five small molecules. Of the four methods, the scaled opposite-spin MP2 (SOS-MP2), the variable-scaling opposite-spin MP2 (VOS-MP2) and the spin-component-scaled MP2 (SCS-MP2) methods perform statistically better than standard MP2 theory, while the spin-component scaled for nucleic bases MP2 (SCSN-MP2) performs worse. Vibrations of closed-shell diatomic molecules are slightly more accurately described by the SOS-MP2 method of Head-Gordon (ε(MAD) =51 cm(-1) ) than the SCS-MP2 method of Grimme (ε(MAD) =61 cm(-1)) or the size-consistent parametrisation of VOS-MP2 (ε(MAD) =54 cm(-1)). For open-shell diatomic molecules, the SOS-MP2 (ε(MAD) =83 cm(-1)) and SCS-MP2 (ε(MAD) =81 cm(-1)) methods are of similar accuracy, while VOS-MP2 is slightly better (ε(MAD) =77 cm(-1)). Since the VOS-MP2 and SOS-MP2 methods tend to have smaller deviations from experiment, and they can be made computationally more economical than the SCS-MP2 or MP2 methods, we suggest that they should be the preferred ab initio method for computing vibrational frequencies in large molecules.     

Quantifying the symmetry content of the electronic structure of molecules: molecular orbitals and the wave function

A scheme to quantify the symmetry content of the electronic wave function and molecular orbitals for arbitrary molecules is developed within the formalism of Continuous Symmetry Measures (CSMs). After defining the symmetry operation expectation values (SOEVs) as the key quantity to gauge the symmetry content of molecular wavefunctions, we present the working equations to be implemented in order to carry out real calculations using standard quantum chemistry software. The potentialities of a symmetry analysis using this new method are shown by means of some illustrative examples such as the changes induced in the molecular orbitals of a diatomic molecule by an electronegativity perturbation, the breaking of orbital symmetry along the dissociation path of the H(2) molecule, the changes in the molecular orbitals upon a geometrical distortion of the benzene molecule, and the inversion symmetry content in the different spin states of the [Fe(CH(3))(4)](2-) complex.     

The accuracy of dipole moments from spin-component scaled CC2 in ground and electronically excited states

The accuracy of dipole moments calculated from wave function methods based on second-order perturbation theory is investigated in the ground and electronically excited states. Results from the approximate coupled-cluster singles-and-doubles model, CC2, M?ller-Plesset perturbation theory, MP2, and the algebraic diagrammatic construction through second-order, ADC(2), are discussed together with the spin-component scaled and the scaled opposite-spin variants of these methods. The computed dipole moments show a very good correlation with data from high-resolution spectroscopy. Compared to the unscaled methods, the spin-component scaling increases the accuracy of the results and improves the robustness of the calculations. An accuracy about 0.2 to 0.1 D in the ground state and about 0.3 to 0.2 D in the electronically excited states can be achieved with these approaches.     

The consequences of neglecting permutation symmetry in the description of many-electrons systems

The consequences of neglecting the permutation symmetry of the Hamiltonian of many-electrons system are examined. From the comparison of wave functions based on methods, which take (generalized valence bond [GVB]) and do not take (Hartree-Fock) the permutation symmetry into account, it is shown that neglecting the permutation symmetry leads to false concepts, misinterpretations, and unjustifiable approximations when dealing with many-electrons systems, atoms, and molecules. In particular, it is shown that how the double occupancy of atomic and molecular orbitals, the exchange integral, the correlation energy, and the so-called “nondynamic” correlation energy are related to neglecting the permutation symmetry.     

Semiempirical double-hybrid density functional with improved description of long-range correlation

The recently proposed new family of "double-hybrid" density functionals [Grimme, S. J. Chem. Phys. 2006, 124, 34108] replaces a fraction of the semi-local correlation energy by a non-local correlation energy expression that employs the Kohn-Sham orbitals in second-order many-body perturbation theory. These functionals have provided results of high accuracy over a wide range of properties but fail to accurately describe long-range van der Waals interactions. In this work, a distance-dependent scaling factor for the non-local correlation energy is introduced to address this problem, and two new double-hybrid density functionals are proposed. The new functionals are optimized with the finite cc-pVTZ basis on training sets of atomization energies and intermolecular interaction energies. They are compared against (scaled) second-order M?ller-Plesset perturbation theories and popular density functionals including the hybrid-GGA functional B3-LYP and the first double-hybrid functional (B2-PLYP). Tests are performed on an extensive set including reaction energies, barrier heights, weakly interacting complexes, transition-metal systems, molecular geometries, and harmonic vibrational frequencies. Within the cc-pVTZ atomic orbital basis, we have demonstrated the ability to find a parametrization scheme which is simultaneously able to describe thermochemistry and weakly bound systems with a satisfactory degree of accuracy.     

Real versus artifactual symmetry-breaking effects in Hartree-Fock, density-functional, and coupled-cluster methods

We have examined the relative abilities of Hartree-Fock, density-functional theory (DFT), and coupled-cluster theory in describing second-order (pseudo) Jahn-Teller (SOJT) effects, perhaps the most commonly encountered form of symmetry breaking in polyatomic molecules. As test cases, we have considered two prototypical systems: the 2Sigmau+ states of D( infinity h) BNB and C3+ for which interaction with a low-lying 2Sigmag+ excited state leads to symmetry breaking of the nuclear framework. We find that the Hartree-Fock and B3LYP methods correctly reproduce the pole structure of quadratic force constants expected from exact SOJT theory, but that both methods appear to underestimate the strength of the coupling between the electronic states. Although the Tamm-Dancoff (CIS) approximation gives excitation energies with no relationship to the SOJT interaction, the random-phase-approximation (RPA) approach to Hartree-Fock and time-dependent DFT excitation energies predicts state crossings coinciding nearly perfectly with the positions of the force constant poles. On the other hand, the RPA excited-state energies exhibit unphysical curvature near their crossings with the ground (reference) state, a problem arising directly from the mathematical structure of the RPA equations. Coupled-cluster methods appear to accurately predict the strength of the SOJT interactions between the 2Sigmau+ and 2Sigmag+ states, assuming that the inclusion of full triple excitations provides a suitable approximation to the exact wave function, and are the only methods examined here which predict symmetry breaking in BNB. However, coupled-cluster methods are plagued by artifactual force constant poles arising from the response of the underlying reference molecular orbitals to the geometric perturbation. Furthermore, the structure of the "true" SOJT force constant poles predicted by coupled-cluster methods, although correctly positioned, has the wrong structure.     

Hybrid correlation models based on active-space partitioning: seeking accurate O(N5) ab initio methods for bond breaking

Moller-Plesset second-order (MP2) perturbation theory remains the least expensive standard ab initio method that includes electron correlation, scaling as O(N5) with the number of molecular orbitals N. Unfortunately, when restricted Hartree-Fock orbitals are employed, the potential energy curves calculated with this method are of little use at large interatomic separations because of the divergent behavior of MP2 in these regions. In our previous study [J. Chem. Phys. 122, 234110 (2005)] we combined the MP2 method with the singles and doubles coupled cluster (CCSD) method to produce a hybrid method that retains the computational scaling of MP2 and improves dramatically the shape of the MP2 curves. In this work we expand the hybrid methodology to several other schemes. We investigate a new, improved MP2-CCSD method as well as a few other O(N5) methods related to the Epstein-Nesbet pair correlation theory. Nonparallelity errors across the dissociation curve as well as several spectroscopic constants are computed for BH, HF, H2O, CH+, CH4, and Li2 molecules with the 6-31G* basis set and compared with the corresponding full configuration interaction results. We show that among the O(N5) methods considered, our new hybrid MP2-CCSD method is the most accurate and significantly outperforms MP2 not only at large interatomic separations, but also near equilibrium geometries.     

Constrained-pairing mean-field theory. V. Triplet pairing formalism

Describing strong (also known as static) correlation caused by degenerate or nearly degenerate orbitals near the Fermi level remains a theoretical challenge, particularly in molecular systems. Constrained-pairing mean-field theory has been quite successful, capturing the effects of static correlation in bond formation and breaking in closed-shell molecular systems by using singlet electron entanglement to model static correlation at mean-field computational cost. This work extends the previous formalism to include triplet pairing. Additionally, a spin orbital extension of the "odd-electron" formalism is presented as a method for understanding electron entanglement in molecules.     

Direct energy functional minimization under orthogonality constraints

The direct energy functional minimization problem in electronic structure theory, where the single-particle orbitals are optimized under the constraint of orthogonality, is explored. We present an orbital transformation based on an efficient expansion of the inverse factorization of the overlap matrix that keeps orbitals orthonormal. The orbital transformation maps the orthogonality constrained energy functional to an approximate unconstrained functional, which is correct to some order in a neighborhood of an orthogonal but approximate solution. A conjugate gradient scheme can then be used to find the ground state orbitals from the minimization of a sequence of transformed unconstrained electronic energy functionals. The technique provides an efficient, robust, and numerically stable approach to direct total energy minimization in first principles electronic structure theory based on tight-binding, Hartree-Fock, or density functional theory. For sparse problems, where both the orbitals and the effective single-particle Hamiltonians have sparse matrix representations, the effort scales linearly with the number of basis functions N in each iteration. For problems where only the overlap and Hamiltonian matrices are sparse the computational cost scales as O(M2N), where M is the number of occupied orbitals. We report a single point density functional energy calculation of a DNA decamer hydrated with 4003 water molecules under periodic boundary conditions. The DNA fragment containing a cis-syn thymine dimer is composed of 634 atoms and the whole system contains a total of 12,661 atoms and 103,333 spherical Gaussian basis functions.     

Accuracy and limitations of second-order many-body perturbation theory for predicting vertical detachment energies of solvated-electron clusters

Vertical electron detachment energies (VDEs) are calculated for a variety of (H(2)O)(n)(-) and (HF)(n)(-) isomers, using different electronic structure methodologies but focusing in particular on a comparison between second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster theory with noniterative triples, CCSD(T). For the surface-bound electrons that characterize small (H(2)O)(n)(-) clusters (n< or = 7), the correlation energy associated with the unpaired electron grows linearly as a function of the VDE but is unrelated to the number of monomers, n. In every example considered here, including strongly-bound "cavity" isomers of (H(2)O)(24)(-), the correlation energy associated with the unpaired electron is significantly smaller than that associated with typical valence electrons. As a result, the error in the MP2 detachment energy, as a fraction of the CCSD(T) value, approaches a limit of about -7% for (H(2)O)(n)(-) clusters with VDEs larger than about 0.4 eV. CCSD(T) detachment energies are bounded from below by MP2 values and from above by VDEs calculated using second-order many-body perturbation theory with molecular orbitals obtained from density functional theory. For a variety of both strongly- and weakly-bound isomers of (H(2)O)(20)(-) and (H(2)O)(24)(-), including both surface states and cavity states, these bounds afford typical error bars of +/-0.1 eV. We have found only one case where the Hartree-Fock and density functional orbitals differ qualitatively; in this case the aforementioned bounds lie 0.4 eV apart, and second-order perturbation theory may not be reliable.     

Treating molecules in arbitrary spin states using the parametric two-electron reduced-density-matrix method

Minimizing the electronic energy with respect to a parameterized two-electron reduced density matrix (2-RDM) is known as a parametric variational 2-RDM method. The parametric 2-RDM method with the M 2-RDM parametrization [D. A. Mazziotti, Phys. Rev. Lett. 101, 253002 (2008)] is extended to treat molecules in arbitrary spin states. Like its singlet counterpart, the M parametric 2-RDM method for arbitrary spin states is derived using approximate N-representability conditions, which allow it to capture more correlation energy than coupled cluster with single and double excitations at a lower computational cost. We present energies, optimized bond lengths, potential energy curves, and occupation numbers for a set of molecules in a variety of spin states using the M and K parametric 2-RDM methods as well as several wavefunction methods. We show that the M parametric 2-RDM method can describe bond breaking of open-shell molecules like triplet B(2) and singlet and triplet OH(+) even in the presence of strong correlation. Finally, the computed 2-RDMs are shown to be nearly N-representable at both equilibrium and non-equilibrium geometries.     

Orbital-dependent correlation energy in density-functional theory based on a second-order perturbation approach: success and failure

We have developed a second-order perturbation theory (PT) energy functional within density-functional theory (DFT). Based on PT with the Kohn-Sham (KS) determinant as a reference, this new ab initio exchange-correlation functional includes an exact exchange (EXX) energy in the first order and a correlation energy including all single and double excitations from the KS reference in the second order. The explicit dependence of the exchange and correlation energy on the KS orbitals in the functional fits well into our direct minimization approach for the optimized effective potential, which is a very efficient method to perform fully self-consistent calculations for any orbital-dependent functionals. To investigate the quality of the correlation functional, we have applied the method to selected atoms and molecules. For two-electron atoms and small molecules described with small basis sets, this new method provides excellent results, improving both second-order Moller-Plesset expression and any conventional DFT results significantly. For larger systems, however, it performs poorly, converging to very low unphysical total energies. The failure of PT based energy functionals is analyzed, and its origin is traced back to near degeneracy problems due to the orbital- and eigenvalue-dependent algebraic structure of the correlation functional. The failure emerges in the self-consistent approach but not in perturbative post-EXX calculations, emphasizing the crucial importance of self-consistency in testing new orbital-dependent energy functionals.