Full CI benchmark calculations for molecular properties

Full CI calculations of first- and second-order properties are presented to provide benchmark results for comparisons with other methods, such as multireference CI(MRCI). The full CI(FCI) polarizability of F^– is computed using a double zeta plus polarization plus diffuse basis set. These FCI results are compared to those obtained at other levels of theory; the CASSCF/MRCI with Davidson correction results are in excellent agreement with the FCI. Differences between the polarizability results computed as a (numerical) second derivative of the energy or as an induced dipole moment are also discussed. FCI calculations are presented for the dipole moment and polarizability of HF, CH₂ and SiH₂ using a DZP basis set. Again, the CASSCF/MRCI values are in excellent agreement with the FCI results, whereas SDCI values, whether computed as an expectation value or as an energy derivative, are much worse. The results obtained using the CPF approach are in considerably better agreement with the FCI results than SDCI, and are similar in quality to the SDCI energy derivative results with the inclusion of Davidson's correction.     

Heisenberg behavior of some carbon‐beryllium compounds: How well truncated‐CI approaches work

This works tries to establish the performance of truncated CI calculations on the evaluation of magnetic coupling parameters with respect to available FCI estimates on a set of carbon‐beryllium clusters. First‐, second‐ and third‐neighbor magnetic coupling constants have been evaluated and many body effective parameters as the cyclic terms. They result from the fitting of the low‐lying states to the eigenvalues of an extended Heisenberg Hamiltonian, involving not only two‐body isotropic terms but also cyclic terms. SDCI and DDCI calculations have been carried out and their performance compared with FCI ones. The impact of the basis set choice and size‐consistency errors have been explored. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Modeling the adiabatic connection in H2

Full configuration interaction (FCI) data are used to quantify the accuracy of approximate adiabatic connection (AC) forms in describing the ground state potential energy curve of H2, within spin-restricted density functional theory (DFT). For each internuclear separation R, accurate properties of the AC are determined from large basis set FCI calculations. The parameters in the approximate AC form are then determined so as to reproduce these FCI values exactly, yielding an exchange-correlation energy expressed entirely in terms of FCI-derived quantities. This is combined with other FCI-derived energy components to give the total electronic energy; comparison with the FCI energy quantifies the accuracy of the AC form. Initial calculations focus on a [1/1]-Padé-based form. The potential energy curve determined using the procedure is a notable improvement over those from existing DFT functionals. The accuracy near equilibrium is quantified by calculating the bond length and vibrational wave numbers; errors in the latter are below 0.5%. The molecule dissociates correctly, which can be traced to the use of virtual orbital eigenvalues in the slope in the noninteracting limit, capturing static correlation. At intermediate R, the potential energy curve exhibits an unphysical barrier, similar to that noted previously using the random phase approximation. Alternative forms of the AC are also considered, paying attention to size extensivity and the behavior in the strong-interaction limit; none provide an accurate potential energy curve for all R, although good accuracy can be achieved near equilibrium. The study demonstrates how data from correlated ab initio calculations can provide valuable information about AC forms and highlight areas where further theoretical progress is required.     

Application of the iterative difference-dedicated configuration interaction method to the determination of excitation energies in some benchmark systems: Be, CH+, BH and CH2

The iterative difference-dedicated CI method (IDDCI) has been applied to determine excitation energies in small systems for which benchmark FCI and other high-level calculations exist. Transitions to excited singlet and triplet states in Be and vertical transitions in CH⁺, BH and CH₂ are reported. The deviations from FCI results are lower than 0.1 eV and compare advantageously with SDCI including size-consistency corrections, (SC)²SDCI, and with coupled cluster calculations including the effect of triples, especially for the states which have a predominant double excitation character. The IDDCI procedure has been speeded up by using smaller subspaces for optimizing the molecular orbitals. Received: 17 January 1997 / Accepted: 31 July 1997     

A sparse matrix based full-configuration interaction algorithm

We present an algorithm related to the full-configuration interaction (FCI) method that makes complete use of the sparse nature of the coefficient vector representing the many-electron wave function in a determinantal basis. Main achievements of the presented sparse FCI (SFCI) algorithm are (i) development of an iteration procedure that avoids the storage of FCI size vectors; (ii) development of an efficient algorithm to evaluate the effect of the Hamiltonian when both the initial and the product vectors are sparse. As a result of point (i) large disk operations can be skipped which otherwise may be a bottleneck of the procedure. At point (ii) we progress by adopting the implementation of the linear transformation by Olsen et al. [J. Chem Phys. 89, 2185 (1988)] for the sparse case, getting the algorithm applicable to larger systems and faster at the same time. The error of a SFCI calculation depends only on the dropout thresholds for the sparse vectors, and can be tuned by controlling the amount of system memory passed to the procedure. The algorithm permits to perform FCI calculations on single node workstations for systems previously accessible only by supercomputers.     

Electron correlation effects in cobalt fluorides CoFn

The molecular cobalt fluorides CoF₂, CoF₃ and CoF₄ are studied and compared by employing different basis sets as well as Quantum Information Theory (QIT) to investigate their correlation effects. These prototypical monomers may be systematically extended in size yielding a novel quasi 1-dimensional, strongly correlated model system consisting of cobalt atoms bridged by oxygen atoms and fluorine termination on both ends. Accurate correlation energies are obtained using Full Configuration Interaction (FCI) and Full Configuration Interaction Quantum Monte Carlo (FCIQMC) calculations and the results are compared to Coupled Cluster and Density Matrix Renormalization Group (DMRG) energies. The analysis indicates the cobalt atom requires a larger number of one-electron basis functions than fluorine and the use of localized molecular orbitals may facilitate calculations for the extended systems.     

A full-configuration interaction "nuclear orbital" method to study doped 3HeN clusters (N< or =4)

An efficient full configuration interaction (FCI) treatment, based on the Jacobi-Davidson algorithm, is developed in order to study small doped (3)He(N) clusters. The state of each He atom in a given cluster is described by a set of wave-functions which by extention of the quantum-chemistry notation are caller here "nuclear orbitals". The FCI treatment is applied to the calculation of binding energies and helium natural orbitals of (3)He(N)...Br(2)(X) complexes. In agreement with our previous calculations using a Hartree-Fock approach [Phys. Rev. Lett. 93, 053401 (2004)], in which the He-He interaction is modified at small distances to account for short-range correlation effects, the lowest-energy states of each multiplet are found to be very close in energy. The natural orbital analysis, in turn, indicates the adequacy of the "nuclear orbital" approach in these systems.     

Potential energy curves for the X1∑ g, B1△g and B''''1∑ g states of C2 using MRCI and approximate CI methods

The potential energy curves for the X1∑ g, B1△g and B′1∑ g states of C2 have been studied by using MRCI and approximate CI methods, and are benchmarked against the calculations of full configuration interaction (FCI). The results obtained by MRCI method agree with the FCI very well, and even are accurate enough to compare other approximate methods as benchmark, when the calculations of FCI are not feasible. The approximate CI methods mentioned in this paper are reliable for treating chemical problems.     

Full configuration interaction benchmark calculations for transition moments

Full configuration-interaction (FCI) calculations have been performed for the ã ¹A₁–b¹B₁and ã ¹A₁–(2)¹A₁transitions in CH₂ and for selected dipole and quadrupole transitions in BeO. The FCI transition moments are compared to those obtained from correlation treatments that truncate the n-particle expansion. The state-averaged MCSCF/SOCI and FCI results agree well, even for BeO, where the CASSCF level nonorthogonal transition moment differs from the state-averaged CASSCF transition moment.     

Interactions in diatomic dimers involving closed-shell metals

Interaction energies of dimers containing alkaline earth (Be, Mg, and Ca) metals have been investigated using symmetry-adapted perturbation theory (SAPT) and supermolecular (SM) methods. Also, to enable broader comparisons, some calculations have been performed on the Zn dimer and on the He-Mg dimer. Although all of the investigated metallic atoms have closed electronic shells, the quasidegeneracy of the ground states of these atoms with the lowest-lying excited states leads to convergence problems in theories based on a single-determinant reference state. The main goal of the present work was to establish how the quality of the interaction energies computed using various electronic-structure methods changes across the range of atoms. We show that although the convergence problems become somewhat less severe with the increase of the atomic number, single-determinant-based methods do not provide reliable interaction energies for any of the investigated metallic dimers even at the level of the coupled-cluster method with single, double, and noniterative triple excitations [CCSD(T)]. However, interaction energies accurate to within a few percent can be obtained if CCSD(T) calculations in large basis sets are extrapolated to the complete basis set limit and followed by full configuration interaction (FCI) calculations with a frozen-core (FC) approximation. Since the systems considered contain only two valence electrons, FCI/FC calculations have been feasible for all of them except for Zn2, providing the best theoretical estimates of the binding energies to date. We found that a large part of the error of the SAPT results originates from limiting some exchange components to terms proportional to the squares of the intermonomer orbital overlap integrals. When the neglected terms were approximately accounted for, the accuracy improved significantly and became comparable to that of CCSD(T), allowing us to obtain for the first time a physical interpretation of the interaction energies in metallic dimers.     

Accurate pair interaction energies for helium from supermolecular Gaussian geminal calculations

Nonrelativistic clamped-nuclei pair interaction energy for ground-state helium atoms has been computed for 12 interatomic separations ranging from 3.0 to 9.0 bohr. The calculations applied the supermolecular approach. The major part of the interaction energy was obtained using the Gaussian geminal implementation of the coupled-cluster theory with double excitations (CCD). Relatively small contributions from single, triple, and quadruple excitations were subsequently included employing the conventional orbital coupled-cluster method with single, double, and noniterative triple excitations [CCSD(T)] and the full configuration interaction (FCI) method. For three distances, the single-excitation contribution was taken from literature Gaussian-geminal calculations at the CCSD level. The orbital CCSD(T) and FCI calculations used very large basis sets, up to doubly augmented septuple- and sextuple-zeta size, respectively, and were followed by extrapolations to the complete basis set limits. The accuracy of the total interaction energies has been estimated to be about 3 mK or 0.03% at the minimum of the potential well. For the attractive part of the well, the relative errors remain consistently smaller than 0.03%. In the repulsive part, the accuracy is even better, except, of course, for the region where the potential goes through zero. For interatomic separations smaller than 4.0 bohr, the relative errors do not exceed 0.01%. Such uncertainties are significantly smaller than the expected values of the relativistic and diagonal Born-Oppenheimer contributions to the potential.     

A graphical symmetric group approach for a spin adapted full configuration interaction: partitioning of a configuration graph into sets of closed-shell and open-shell graphs

We developed a spin adapted full configuration interaction (FCI) method which was expected to be effective for parallel processing. The graphical symmetric group approach (GSGA) was employed, where a configuration graph was partitioned into several sets of closed-shell and open-shell graphs. The configuration state functions (CSFs) bearing the same number of closed-shells and open-shells were assembled in a configuration group. The graphical approach provided a number to identify each CSF in a sequential order within the group. Combination of this partitioning and an intermediate configuration-driven algorithm in calculating the so-called σ vectors allowed us to use symbolic coupling constants. Furthermore, this combination made it easy to implement an efficient algorithm suitable to task-distributed parallel procedure for evaluating σ vectors. A program was written and some test calculations were carried out with high parallel efficiency. The largest size of FCI used 10 million CSFs (20 million determinants).     

Full configuration interaction study of the ground state of closed-shell cyclic PPP polyenes

Full configuration interaction (FCI ) calculations are reported for the closed-shell cyclic polyenes C_NH_N (N = 6, 10, 14, 18) in the Pariser–Parr–Pople (PPP ) approximation, as a function of the hopping parameter β. A wide range of values of β is considered, from a highly correlated situation, β = 0, to a very weakly correlated limit, β = ?10 eV. An estimate of the role of higher than two-body connected cluster components was done, through a partial (i.e., limited to two-body terms in the exponent) cluster analysis performed on the FCI wave function. The comparison with the approximate coupled pair theory that accounts for quadruply excited clusters [see P. Piecuch and J. Paldus, Theor. Chim. Acta 78 , 65 (1990)] shows a good agreement in the whole range of the hopping parameter, particularly when the contribution of connected triples excitations is also taken into account. © 1994 John Wiley & Sons, Inc.     

Full configuration interaction calculation of singlet excited states of Be3

The full configuration interaction (FCI) study of the singlets vertical spectrum of the neutral beryllium trimer has been performed using atomic natural orbitals [3s2p1d] basis set. The FCI triangular equilibrium structure of the ground state has been used to calculate the FCI vertical excitation energies up to 4.8 eV. The FCI vertical ionization potential for the same geometry and basis set amounts to 7.6292 eV. The FCI dipole and quadrupole transition moments from the ground state are reported as well. The FCI electric quadrupole moment of the X (3)A(1) (') ground state has been also calculated with the same basis set (Theta(zz)=-2.6461 a.u., Theta(xx)=Theta(yy)=-1/2Theta(zz)). Twelve of the 19 calculated excited singlets are doubly excited states. Most of the states have large multiconfigurational character. These results provide benchmark values for electronic correlation multireference methods. (4ex6MO)CAS-SDCI values for the same energies and properties are also reported.     

Comparison of low-order multireference many-body perturbation theories

Tests have been made to benchmark and assess the relative accuracies of low-order multireference perturbation theories as compared to coupled cluster (CC) and full configuration interaction (FCI) methods. Test calculations include the ground and some excited states of the Be, H(2), BeH(2), CH(2), and SiH(2) systems. Comparisons with FCI and CC calculations show that in most cases the effective valence shell Hamiltonian (H(v)) method is more accurate than other low-order multireference perturbation theories, although none of the perturbative methods is as accurate as the CC approximations. We also briefly discuss some of the basic differences among the multireference perturbation theories considered in this work.     

CI calculations of long-range C6 dispersion coefficients for BH–BH

Large scale CI calculations of frequency-dependent dipole polarizabilities (FDPs) of ground state BH at the experimental distance of R = 2.329a₀ have been performed using two sets of contracted GTO functions involving up to 316 × 10⁶ symmetry-adapted Slater determinants, comparing recently developed interpolation and variational techniques. Final calculation of BH–BH dispersion coefficients was done using the larger basis set of 268 AOs in the frozen-core FCI approximation and the lower bounds resulting from the variational approach.