An efficient internally contracted multiconfiguration–reference configuration interaction method
01 Nov 1988-Journal of Chemical Physics (American Institute of Physics)-Vol. 89, Iss: 9, pp 5803-5814
TL;DR: In this article, a new internally contracted direct multiconfiguration-reference configuration interaction (MRCI) method is described which allows the use of much larger reference spaces than any previous MRCI method.
Abstract: A new internally contracted direct multiconfiguration–reference configuration interaction (MRCI) method is described which allows the use of much larger reference spaces than any previous MRCI method. The configurations with two electrons in the external orbital space are generated by applying pair excitation operators to the reference wave function as a whole, while the singly external and internal configurations are standard uncontracted spin eigenfunctions. A new efficient and simple method for the calculation of the coupling coefficients is used, which is well suited for vector machines, and allows the recalculation of all coupling coefficients each time they are needed. The vector H⋅c is computed partly in a nonorthogonal configuration basis. In order to test the accuracy of the internally contracted wave functions, benchmark calculations have been performed for F−, H2O, NH2, CH2, CH3, OH, NO, N2, and O2 at various geometries. The deviations of the energies obtained with internally contracted and uncontracted MRCI wave functions are mostly smaller than 1 mH and typically 3–5 times smaller than the deviations between the uncontracted MRCI and the full CI. Dipole moments, electric dipole polarizabilities, and electronic dipole transition moments calculated with uncontracted and contracted MRCI wave functions also are found to be in close agreement. The efficiency of the method is demonstrated in large scale calculations for the CN, NH3, CO2, and Cr2 molecules. In these calculations up to 3088 reference configurations and up to 154 orbitals were employed. The biggest calculation is equivalent to an uncontracted MRCI with more than 78 million configurations.
Citations
More filters
TL;DR: Molpro (available at http://www.molpro.net) is a general-purpose quantum chemical program as discussed by the authors, which uses local approximations combined with explicit correlation treatments, highly accurate coupled-cluster calculations are now possible for molecules with up to approximately 100 atoms.
Abstract: Molpro (available at http://www.molpro.net) is a general-purpose quantum chemical program. The original focus was on high-accuracy wave function calculations for small molecules, but using local approximations combined with explicit correlation treatments, highly accurate coupled-cluster calculations are now possible for molecules with up to approximately 100 atoms. Recently, multireference correlation treatments were also made applicable to larger molecules. Furthermore, an efficient implementation of density functional theory is available.
2,999 citations
TL;DR: In this article, the correlationconsistent polarized core-valence basis sets (cc•pVXZ) for the atoms boron through neon have been extended to treat core and core•valence correlation effects, and the exponents of the core functions were determined by minimizing the difference between allelectron and valence only correlation energies obtained from HF+1+2 calculations on the ground states of the atoms.
Abstract: The correlation‐consistent polarized valence basis sets (cc‐pVXZ) for the atoms boron through neon have been extended to treat core and core‐valence correlation effects. Basis functions were added to the existing cc‐pVXZ sets to form correlation‐consistent polarized core‐valence sets (cc‐pCVXZ) in the usual pattern: Double zeta added (1s1p), triple zeta added (2s2p1d), quadruple zeta added (3s3p2d1f), and quintuple zeta added (4s4p3d2f1g). The exponents of the core functions were determined by minimizing the difference between all‐electron and valence‐only correlation energies obtained from HF+1+2 calculations on the ground states of the atoms. With the cc‐pCVXZ sets, core, core‐valence, and valence correlation energies all converge exponentially toward apparent complete basis set (CBS) limits, as do the corresponding all‐electron singles and doubles CI energies. Several test applications of the new sets are presented: The first two ionization potentials of boron, the 3P–5S separation in carbon, and the X...
2,666 citations
TL;DR: In this paper, correlation consistent and augmented correlation consistent basis sets for the third row main group atoms gallium through krypton were determined for the gallium atom, and the results showed good convergence to an apparent complete basis set limit.
Abstract: Valence correlation consistent and augmented correlation consistent basis sets have been determined for the third row, main group atoms gallium through krypton. The methodology, originally developed for the first row atoms, was first applied to the selenium atom, resulting in the expected natural groupings of correlation functions (although higher angular momentum functions tend to be relatively more important for the third row atoms as they were for the second row atoms). After testing the generality of the conclusions for the gallium atom, the procedure was used to generate correlation consistent basis sets for all of the atoms gallium through krypton. The correlation consistent basis sets for the third row main group atoms are as follows: cc-pVDZ: (14s11p6d)/[5s4p2d]; cc-pVTZ: (20s13p9d1f )/[6s5p3d1f]; cc-pVQZ: (21s16p12d2 f1g)/[7s6p4d2 f1g]; cc-pV5Z: (26s17p13d3f2g1h)/[8s7p5d3f2g1h]. Augmented sets were obtained by adding diffuse functions to the above sets (one for each angular momentum present in the set), with the exponents of the additional functions optimized in calculations on the atomic anions. Test calculations on the atoms as well as selected molecules with the new basis sets show good convergence to an apparent complete basis set limit.
2,164 citations
TL;DR: In this paper, the coupled cluster method restricted to single and double excitations (CCSD) is considered for the case of a spin restricted Hartree-Fock open shell reference determinant.
Abstract: The coupled cluster method restricted to single and double excitations (CCSD) is considered for the case of a spin restricted Hartree–Fock open shell reference determinant. A spin–orbital based formulation, in which the cluster operator spans exactly the minimal first order interacting space, is presented, and computationally optimal working equations are given. In the limit of a large number of closed shell orbitals, the cost is identical to that of an optimum treatment of an equivalent closed shell problem, which is obtained as a special case of the formulation presented. The theory is applied to the calculation of a number of diatomic potential energy functions and compared with spin‐unrestricted theory.
1,822 citations
TL;DR: In this paper, the authors used correlation consistent Gaussian basis sets from double to sextuple zeta quality to estimate the barrier height of the H+H2 exchange reaction.
Abstract: Using systematic sequences of correlation consistent Gaussian basis sets from double to sextuple zeta quality, the classical barrier height of the H+H2 exchange reaction has been calculated by multireference configuration interaction (MRCI) methods. The MRCI calculations for collinear H3 have also been calibrated against large basis set full CI (FCI) results, which demonstrate that the MRCI treatment leads to energies less than 1 μhartree (≤0.001 kcal/mol) above the FCI energies. The dependence of both the H2 and H3 total energies on the basis set is found to be very regular, and this behavior has been used to extrapolate to the complete basis set (CBS) limits. The resulting estimate of the H–H–H CBS limit yields a classical barrier height, relative to exact H+H2, of 9.60±0.02 kcal/mol; the best directly calculated value for the barrier is equal to 9.62 kcal/mol. These results are in excellent agreement with recent quantum Monte Carlo calculations.
1,614 citations
References
More filters
TL;DR: In this article, a new method for evaluating one-particle coupling coefficients in general configuration interaction calculation is presented, through repeated application and use of resolutions of the identity, two-, three-and four-body coupling coefficients and density matrices may be built in a simple and efficient way.
2,455 citations
TL;DR: In this paper, a second-order optimisation procedure for general complete active space (CAS) wavefunctions is described. But this method is restricted to very long complete active spaces.
2,365 citations
TL;DR: In this paper, a configuration selection method for CI calculations is discussed and applied in which the energy lowering produced in a secular equation by the addition of a given test species to a series of dominant configurations is used as an ordering parameter.
Abstract: A configuration selection method for CI calculations is discussed and applied in which the energy lowering produced in a secular equation by the addition of a given test species to a series of dominant configurations is used as an ordering parameter. Configurations with energy lowerings below a given energy cut-off value are not included in the final secular equations but instead a method of estimating the combined effect of the neglected species on the corresponding non-selected CI results is developed. The influence of the choice of main configurations used in the selection process is given close examination as well as the importance of the MO basis employed in the treatment as a whole; in the latter case a non-iterative procedure for obtaining approximate natural orbitals for such calculations is suggested. The resulting configuration selection procedure is equally applicable to all types of electronic states in any nuclear geometry and the results of the associated CI calculations are seen to be essentially equivalent to a complete treatment in which all single- and double-excitation species with respect to aseries of dominant configurations in a given state are included.
1,126 citations
TL;DR: In this article, a multi-reference double-excitation CI (MRD-CI) method is discussed and its results are compared with those of related techniques, which employs a configuration selection procedure to order the various generated species according to their energy-lowering capability and then uses an energy extrapolation procedure based on perturbation theory to obtain suitably accurate estimates of the eigenvalues of the entire MRD- CI space.
Abstract: Implementation of a multi-reference double-excitation CI (MRD-CI) method is discussed and its results are compared with those of related techniques. This approach employs a configuration selection procedure to order the various generated species according to their energy-lowering capability and then uses an energy extrapolation procedure based on perturbation theory to obtain suitably accurate estimates of the eigenvalues of the entire MRD-CI space. By employing this selection procedure it is possible to test from 2000 to 4000 symmetry-adapted functions (SAF's) per second of CPU time on an IBM 370–168 system, thereby allowing one to apply the energy extrapolation quite conveniently to CI spaces consisting of several hundred thousand species. By systematically increasing the number of reference configurations in the MRD-CI it is clear that the limit of a full CI can be approached and as a result such a computational procedure appears to be generally valid for any type of electronic state and for any nuclea...
852 citations