Coupled-cluster dynamic polarizabilities including triple excitations

Dynamic polarizabilities for open- and closed-shell molecules were obtained by using coupled-cluster (CC) linear response theory with full treatment of singles, doubles, and triples (CCSDT-LR) with large basis sets utilizing the NWChem software suite. By using four approximate CC methods in conjunction with augmented cc-pVNZ basis sets, we are able to evaluate the convergence in both many-electron and one-electron spaces. For systems with primarily dynamic correlation, the results for CC3 and CCSDT are almost indistinguishable. For systems with significant static correlation, the CC3 tends to overestimate the triples contribution, while the PS(T) approximation [J. Chem. Phys. 127, 164105 (2007)] produces mixed results that are heavily dependent on the accuracies provided by noniterative approaches used to correct the equation-of-motion CCSD excitation energies. Our results for open-shell systems show that the choice of reference (restricted open-shell Hartree-Fock versus unrestricted Hartree-Fock) can have a significant impact on the accuracy of polarizabilities. A simple extrapolation based on pentuple-zeta CCSD calculations and triple-zeta CCSDT calculations reproduces experimental results with good precision in most cases.     

Active-space coupled-cluster methods through connected quadruple excitations

Coupled-cluster methods that include just a subset of all connected triple, quadruple, or both excitation amplitudes, according to the ansatz of and Adamowicz co-workers [Int. Rev. Phys. Chem. 12, 339 (1993); J. Chem. Phys. 99, 1875 (1993); 100, 5792 (1994)] and Piecuch et al. [J. Chem. Phys. 110, 6103 (1999)], have been implemented into parallel execution programs. They are applicable to closed- and open-shell species and they take advantage of real Abelian point-group symmetry. A symbol manipulation program has been invoked to automate the implementation. These methods have been applied to the singlet-triplet separations of five triatomic hydrides (CH2, NH2+, SiH2, PH2+, and AsH2+) with consideration of scalar relativistic effects. They have been shown to be remarkably effective with errors arising from the use of a very small subset of higher-order excitations being no more than a few tenths of 1 kcal/mol.     

General-model-space state-universal coupled-cluster methods for excited states: diagonal noniterative triple corrections

The recently developed multireference, general-model-space, state-universal coupled-cluster approach considering singles and doubles (GMS SU CCSD) has been extended to account perturbatively for triples, similar to the ubiquitous single-reference CCSD(T) method. The effectiveness of this extension in handling of excited states and its ability to account for the static and nondynamic correlation effects when considering spin- and/or space-symmetry degenerate levels within the spin-orbital formalism is examined on the example of low-lying excitation energies of the C2, N2, and CO molecules and a comparison is made with the (N,N)-CCSD method used for the same purpose. It is shown that while the triple corrections are very effective in improving the absolute energies, they have only a modest effect on the corresponding excitation energies, which may be even detrimental if both the ground- and excited-state levels cannot be given a balanced treatment. While the triple corrections help to avoid the symmetry-breaking effects arising due to the use of the spin-orbital formalism, they are much less effective in this regard than the (N,N)-CCSD approach.     

Coupled cluster methods including triple excitations for excited states of radicals

We report an extension of the coupled cluster iterative-triples model, CC3, to excited states of open-shell molecules, including radicals. We define the method for both spin-unrestricted Hartree-Fock (UHF) and spin-restricted open-shell Hartree-Fock (ROHF) reference determinants and discuss its efficient implementation in the PSI3 program package. The program is streamlined to use at most O(N(7)) computational steps and avoids storage of the triple-excitation amplitudes for both the ground- and excited-state calculations. The excitation-energy program makes use of a Lowdin projection formalism (comparable to that of earlier implementations) that allows computational reduction of the Davidson algorithm to only the single- and double-excitation space, but limits the calculation to only one excited state at a time. However, a root-following algorithm may be used to compute energies for multiple states of the same symmetry. Benchmark applications of the new methods to the lowest valence (2)B(1) state of the allyl radical, low-lying states of the CH and CO(+) diatomics, and the nitromethyl radical show substantial improvement over ROHF- and UHF-based CCSD excitation energies for states with strong double-excitation character or cases suffering from significant spin contamination. For the allyl radical, CC3 adiabatic excitation energies differ from experiment by less than 0.02 eV, while for the (2)Sigma(+) state of CH, significant errors of more than 0.4 eV remain.     

A note on the treatment of quadruple excitations in configuration interaction

Perturbation theory is used to justify the approximation in configuration interaction of the coefficients of quadruple excitations as products of coefficients of double excitations. Corresponding energy expressions are presented and tested in a simple application to the beryllium atom. The ideas are likely to be of most use in connection with large configuration interaction calculations.     

Semi-empirical molecular orbital methods including dispersion corrections for the accurate prediction of the full range of intermolecular interactions in biomolecules

Semi-empirical calculations including an empirical dispersive correction are used to calculate intermolecular interaction energies and structures for a large database containing 156 biologically relevant molecules (hydrogen-bonded DNA base pairs, interstrand base pairs, stacked base pairs and amino acid base pairs) for which MP2 and CCSD(T) complete basis set (CBS) limit estimates of the interaction energies are available. The dispersion corrected semi-empirical methods are parameterised against a small training set of 22 complexes having a range of biologically important non-covalent interactions. For the full molecule set (156 complexes), compared to the high-level ab initio database, the mean unsigned errors of the interaction energies at the corrected semi-empirical level are 1.1 (AM1-D) and 1.2 (PM3-D) kcal mol(-1), being a significant improvement over existing AM1 and PM3 methods (8.6 and 8.2 kcal mol(-1)). Importantly, the new semi-empirical methods are capable of describing the diverse range of biological interactions, most notably stacking interactions, which are poorly described by both current AM1 and PM3 methods and by many DFT functionals. The new methods require no more computer time than existing semi-empirical methods and therefore represent an important advance in the study of important biological interactions.     

A comparative assessment of the perturbative and renormalized coupled cluster theories with a noniterative treatment of triple excitations for thermochemical kinetics, including a study of basis set and core correlation effects

The CCSD, CCSD(T), and CR-CC(2,3) coupled cluster methods, combined with five triple-zeta basis sets, namely, MG3S, aug-cc-pVTZ, aug-cc-pV(T+d)Z, aug-cc-pCVTZ, and aug-cc-pCV(T+d)Z, are tested against the DBH24 database of diverse reaction barrier heights. The calculations confirm that the inclusion of connected triple excitations is essential to achieving high accuracy for thermochemical kinetics. They show that various noniterative ways of incorporating connected triple excitations in coupled cluster theory, including the CCSD(T) approach, the full CR-CC(2,3) method, and approximate variants of CR-CC(2,3) similar to the triples corrections of the CCSD(2) approaches, are all about equally accurate for describing the effects of connected triply excited clusters in studies of activation barriers. The effect of freezing core electrons on the results of the CCSD, CCSD(T), and CR-CC(2,3) calculations for barrier heights is also examined. It is demonstrated that to include core correlation most reliably, a basis set including functions that correlate the core and that can treat core-valence correlation is required. On the other hand, the frozen-core approximation using valence-optimized basis sets that lead to relatively small computational costs of CCSD(T) and CR-CC(2,3) calculations can achieve almost as high accuracy as the analogous fully correlated calculations.     

Coupled-cluster theory with simplified linear-r(12) corrections: the CCSD(R12) model

A simplified singles-and-doubles linear-r(12) corrected coupled-cluster model, denoted CCSD(R12), is proposed and compared with the complete singles-and-doubles linear-r(12) coupled-cluster method CCSD-R12. An orthonormal auxiliary basis set is used for the resolution-of-the-identity approximation to calculate three-electron integrals needed in the linear-r(12) Ansatz. Basis-set convergence is investigated for a selected set of atoms and small molecules. In a large basis, the CCSD(R12) model provides an excellent approximation to the full linear-r(12) energy contribution, whereas the magnitude of this contribution is significantly overestimated at the level of second-order perturbation theory.     

Coupled-cluster singles and doubles for extended systems

Coupled-cluster theory with connected single and double excitation operators (CCSD) and related approximations, such as linearized CCSD, quadratic configuration interaction with single and double excitation operators, coupled-cluster with connected double excitation operator (CCD), linearized CCD, approximate CCD, and second- and third-order many-body perturbation theories, are formulated and implemented for infinitely extended one-dimensional systems (polymers), on the basis of the periodic boundary conditions and distance-based screening of integrals, density matrix elements, and excitation amplitudes. The variation of correlation energies with the truncation radii of short- and long-range lattice sums and with the number of wave vector sampling points in the first Brillouin zone is examined for polyethylene, polyacetylene, and polyyne, and is shown to be a function of the degree of pi-electron conjugation or the fundamental band gaps. The t2 and t1 amplitudes in the atomic orbital (AO) basis are obtained by first computing the t amplitudes in the Bloch-orbital basis and subsequently back-transforming them into the AO basis. The plot of these AO-based t amplitudes as a function of unit cells also indicates that the t2 amplitudes of polyacetylene and polyyne exhibit appreciably slower decay than those of polyethylene, although the asymptotic decay behavior is invariably 1/r3. The AO-based t1 amplitudes appear to correlate strongly with the electronic structure, and they decay seemingly exponentially for polyethylene whereas they stay at a constant magnitude across the seventh nearest neighbors of polyacetylene and polyyne, which attests to far reaching effects of nondynamical electron correlation mediated by orbital rotation. Nonetheless, the unit cell contributions to the correlation energies taper below 10(-6) hartree after 15 A for all three polymers. The basis set dependence of the decay behavior of t2 amplitudes is also examined for linear hydrogen fluoride polymer (HF)infinity and linear beryllium polymer (Be)infinity employing the STO-3G, 6-31G, and 6-31G* basis sets, and proves to be rather small.     

The non-covalent functionalisation of carbon nanotubes studied by density functional and semi-empirical molecular orbital methods including dispersion corrections

Density functional theory (DFT-D) and semi-empirical (PM3-D) methods having an added empirical dispersion correction have been used to study the binding of a series of small molecules and planar aromatic molecules to single-walled carbon nanotubes (CNTs). For the small molecule set, the PM3-D method gives a mean unsigned error (MUE) in the binding energies of 1.2 kcal mol(-1) when judged against experimental reference data for graphitic carbon. This value is close to the MUE for this method compared to high-level ab initio data for biological complexes. The PM3-D and DFT-D calculations describing the adsorption of the planar organic molecules (benzene, bibenzene, naphthalene, anthracene, TCNQ and DDQ) on the outer-walls of both semi-conducting and metallic CNTs give similar binding energies for benzene and DDQ, but do not display a stronger adsorption on [6,6] compared to [10,0] structures shown by another DFT study.     

Fast noniterative orbital localization for large molecules

We use Cholesky decomposition of the density matrix in atomic orbital basis to define a new set of occupied molecular orbital coefficients. Analysis of the resulting orbitals ("Cholesky molecular orbitals") demonstrates their localized character inherited from the sparsity of the density matrix. Comparison with the results of traditional iterative localization schemes shows minor differences with respect to a number of suitable measures of locality, particularly the scaling with system size of orbital pair domains used in local correlation methods. The Cholesky procedure for generating orthonormal localized orbitals is noniterative and may be made linear scaling. Although our present implementation scales cubically, the algorithm is significantly faster than any of the conventional localization schemes. In addition, since this approach does not require starting orbitals, it will be useful in local correlation treatments on top of diagonalization-free Hartree-Fock optimization algorithms.     

Universal state-selective corrections to multi-reference coupled-cluster theories with single and double excitations

The recently proposed universal state-selective (USS) corrections [K. Kowalski, J. Chem. Phys. 134, 194107 (2011)] to approximate multi-reference coupled-cluster (MRCC) energies can be commonly applied to any type of MRCC theory based on the Jeziorski-Monkhorst [B. Jeziorski and H. J. Monkhorst, Phys. Rev. A 24, 1668 (1981)] exponential ansatz. In this paper we report on the performance of a simple USS correction to the Brillouin-Wigner and Mukherjee's MRCC approaches employing single and double excitations (USS-BW-MRCCSD and USS-Mk-MRCCSD). It is shown that the USS-BW-MRCCSD correction, which employs the manifold of single and double excitations, can be related to a posteriori corrections utilized in routine BW-MRCCSD calculations. In several benchmark calculations we compare the USS-BW-MRCCSD and USS-Mk-MRCCSD results with the results obtained with the full configuration interaction method.     

Coupled-cluster methods with perturbative inclusion of explicitly correlated terms: a preliminary investigation

We propose to account for the large basis-set error of a conventional coupled-cluster energy and wave function by a simple perturbative correction. The perturbation expansion is constructed by L?wdin partitioning of the similarity-transformed Hamiltonian in a space that includes explicitly correlated basis functions. To test this idea, we investigate the second-order explicitly correlated correction to the coupled-cluster singles and doubles (CCSD) energy, denoted here as the CCSD(2)(R12) method. The proposed perturbation expansion presents a systematic and easy-to-interpret picture of the "interference" between the basis-set and correlation hierarchies in the many-body electronic-structure theory. The leading-order term in the energy correction is the amplitude-independent R12 correction from the standard second-order M?ller-Plesset R12 method. The cluster amplitudes appear in the higher-order terms and their effect is to decrease the basis-set correction, in accordance with the usual experience. In addition to the use of the standard R12 technology which simplifies all matrix elements to at most two-electron integrals, we propose several optional approximations to select only the most important terms in the energy correction. For a limited test set, the valence CCSD energies computed with the approximate method, termed , are on average precise to (1.9, 1.4, 0.5 and 0.1%) when computed with Dunning's aug-cc-pVXZ basis sets [X = (D, T, Q, 5)] accompanied by a single Slater-type correlation factor. This precision is a roughly an order of magnitude improvement over the standard CCSD method, whose respective average basis-set errors are (28.2, 10.6, 4.4 and 2.1%). Performance of the method is almost identical to that of the more complex iterative counterpart, CCSD(R12). The proposed approach to explicitly correlated coupled-cluster methods is technically appealing since no modification of the coupled-cluster equations is necessary and the standard M?ller-Plesset R12 machinery can be reused.     

Many-body perturbation theory,coupled-pair many-electron theory,and the importance of quadruple excitations for the correlation problem

Many-body (diagrammatic) perturbation theory (MBPT ), coupled-pair many-electron theory (CPMET ), and configuration interaction (CI ) are investigated with particular emphasis on the importance of quadruple excitations in correlation theories. These different methods are used to obtain single, double, and quadruple excitation contributions to the correlation energy for a series of molecules including CO₂, HCN, N₂, CO, BH₃, and NH₃. It is demonstrated that the sum of double and quadruple excitation diagrams through fourth-order perturbation theory is usually quite close to the CPMET result for these molecules at equilibrium geometries. The superior reliability of the CPMET model as a function of internuclear separation is illustrated by studying the ¹∑ potential curve of Be₂. This molecule violates the assumption common to nondegenerate perturbation theory that only a single reference function is important and this causes improper behavior of the potential curve as a function of R. This is resolved once the quadruple excitation terms are fully included by CPMET .     

The lengths of molybdenum to molybdenum quadruple bonds: correlations, explanations, and corrections

A systematic search of the Cambridge Crystallographic Database has yielded structures of 467 molecules containing, or reported to contain, a molybdenum to molybdenum quadruple bond, Mo[quadruple bond]Mo. We have arranged these data as a histogram, in order to determine the "normal" range of Mo[quadruple bond]Mo lengths, 2.06-2.17 A, as well as to examine molecules in which the Mo[quadruple bond]Mo length deviates from the norm. We discuss 24 molecules with exceptionally long Mo[quadruple bond]Mo lengths. Another molecule, Mo(2)(pyNC(O)CH(3))(4), was previously reported to have an unusually short Mo[quadruple bond]Mo bond. However, the structure is highly disordered (probably twinned); thus we have prepared and determined the structure of an analogous molecule, Mo(2)(pyNC(O)CH(2)CH(3))(4), which displays a completely normal Mo[quadruple bond]Mo bond length of over 2.08 A.     

An activation analysis system for short-lived radioisotopes including automatic dead-time corrections with a microcomputer

A system based on an IBM-PC microcomputer coupled to a Canberra Series 80 multichannel analyser has been developed for activation analysis with short-lived radioisotopes. The data transfer program can store up to 77 gamma-ray spectra on a floppy disk. A spectrum analysis program, DVC, has been written to determine peak areas interactively, to correct the counting losses, and to calculate elemental concentrations.     

Simple methods for dead-time and pile-up corrections in analytical gamma-ray spectrometry

In the paper the mathematical methods for dead-time and pile-up corrections are discussed. The dead-time correction formulae for a system of two short-lived isotopes and constant component (background) are reported which mathematical problem has not been solved so far. The new electronic circuit for simultaneous measurements of clock- and live- (or dead-) times is described. It is shown that using this circuit one can correct the counting loss for both effects simultaneously. Finally, the advantages and disadvantages of mathematical methods for dead-time and pile-up corrections are discussed based on the authors' several-year laboratory practice.