Ab initio many-body investigation of structure and stability of two-fold rings in silicates

In this paper we present ab initio many-body calculations on the strain energy of W silica, taken as a model system for edge-sharing tetrahedral SiO(2) systems with respect to corner-sharing ones as in alpha quartz. The mean-field results were obtained using the restricted Hartree-Fock approach, while the many-body effects were taken into account by the second-order M?ller-Plesset perturbation theory and the coupled-cluster approach. Correlation contributions are found to play an important role to determine the stability of edge-sharing units. The most sophisticated method used in our calculation, i.e., the coupled-cluster approach with single and double excitations, yields a strain energy of 0.0427 a.u. per Si(2)O(4) unit with respect to alpha quartz, which is even smaller than the value obtained by a previous density functional theory calculation.     

Importance of structural motifs in liquid hydrogen fluoride

The investigation of liquid phases by means of accurate electronic structure methods is a demanding task due to the high computational effort. We applied second-order M?ller-Plesset perturbation theory and high-level quantum chemical calculations using the coupled-cluster method with single, double and perturbative triple excitations in combination with Dunnings correlation-consistent basis sets up to quintuple ζ quality. Based on these calculations, we extrapolated the correlation energy to the basis set limit in order to improve the results even further. For comparison to the correlated electronic structure methods, density functional calculations employing different functionals are presented as well. The investigated species are a cyclic pentamer as well as a set of branched structures. The quantum cluster equilibrium method is employed for the investigation of the liquid-phase structure of hydrogen fluoride. The pentamer is found to be present to a high extent and in the case of the MP2/QZVP data, its presence improves the results significantly. Accounting for branched structures slightly improves results, so that they are found to be present but not to dominate in liquid hydrogen fluoride. Concerning both the interaction energy and the result of the quantum cluster equilibrium calculation the basis set has a major influence, whereas the difference between M?ller-Plesset perturbation theory and coupled-cluster calculations is less pronounced.     

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.     

Intermolecular potential energy surface and second virial coefficients for the water-CO2 dimer

A five-dimensional potential energy surface is calculated for the interaction of water and CO(2), using second-order M?ller-Plesset perturbation theory and coupled-cluster theory with single, double, and perturbative triple excitations. The correlation energy component of the potential energy surface is corrected for basis set incompleteness. In agreement with previous studies, the most negative interaction energy is calculated for a structure with C(2v) symmetry, where the oxygen atom of water is close to the carbon atom of CO(2). Second virial coefficients for the water-CO(2) pair are calculated for a range of temperatures, and their uncertainties are estimated. The virial coefficients are shown to be in close agreement with the available experimental data.     

Conformations of allyl amine: theory vs experiment

The relative stabilities of the five conformers of allyl amine, a medium-size aliphatic molecule, were estimated by applying ab initio quantum mechanical methods at several levels of theory. The second-order M?ller-Plesset perturbation method (MP2), quadratic configuration interaction including single and double excitations (QCISD), coupled-cluster with single and double excitations (CCSD) and CCSD plus perturbative triple excitations [CCSD(T)] were applied. The Dunning correlation consistent basis sets (through aug-cc-pVQZ and cc-pV5Z) were employed. The MP2 energies relative to the energy of the cis-trans conformer reported here appear to approach the basis set limit. The predicted allyl amine conformer energies approaching the Hartree-Fock basis set limit are 158 cm-1 (cis-gauche), -5 cm-1 (gauche-trans), and -146 cm-1 (gauche-gauche). The same three relative energies near the MP2 basis set limit are 135, 103, and 50 cm-1, respectively. The analogous energies deduced from experiment are 173 +/- 12, 92 +/- 8, and 122 +/- 5 cm-1. The theoretical results obtained in the present study suggest that satisfactory predictions of the conformer energetics of allyl amine may be achieved only by theoretical methods that incorporate consideration of correlation effects in conjunction with large basis sets. Evaluation of the zero-point vibrational energy corrections is critical, due to the very small classical energy differences between the five conformers of allyl amine. Agreement between theory and experiment for the gauche-gauche conformational energy remains problematical.     

Accurate ab initio structure, dissociation energy, and vibrational spectroscopy of the F(-)-CH4 anion complex

Accurate equilibrium structure, dissociation energy, global potential energy surface (PES), dipole moment surface (DMS), and the infrared vibrational spectrum in the 0-3000 cm(-1) range of the F(-)-CH4 anion complex have been obtained. The equilibrium electronic structure calculations employed second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster (CC) method up to single, double, triple, and perturbative quadruple excitations using the aug-cc-p(C)VXZ [X = 2(D), 3(T), 4(Q), and 5] correlation-consistent basis sets. The best equilibrium geometry has been obtained at the all-electron CCSD(T)/aug-cc-pCVQZ level of theory. The dissociation energy has been determined based on basis set extrapolation techniques within the focal-point analysis (FPA) approach considering (a) electron correlation beyond the all-electron CCSD(T) level, (b) relativistic effects, (c) diagonal Born-Oppenheimer corrections (DBOC), and (d) variationally computed zero-point vibrational energies. The final D(e) and D0 values are 2398 +/- 12 and 2280 +/- 20 cm(-1), respectively. The global PES and DMS have been computed at the frozen-core CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ levels of theory, respectively. Variational vibrational calculations have been performed for CH4 and F(-)-CH4 employing the vibrational configuration interaction (VCI) method as implemented in Multimode.     

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.     

Interaction energy contributions of H-bonded and stacked structures of the AT and GC DNA base pairs from the combined density functional theory and intermolecular perturbation theory approach

Stacked and Watson-Crick structures of DNA base pairs are investigated with the DFT-SAPT variant of intermolecular perturbation theory, yielding a rigorous decomposition of the interaction energy into electrostatic, induction, dispersion, and exchange contributions. Their interplay in the various structures is analyzed. Total interaction energies extrapolated to the complete basis set limit are compared with corresponding second-order M?ller-Plesset and estimated coupled-cluster theory results.     

An ab initio study of van der Waals potential energy parameters for silver clusters

We employ ab initio calculations of van der Waals complexes to study the potential energy parameters (C(6) coefficients) of van der Waals interactions for modeling of the adsorption of silver clusters on the graphite surface. Electronic structure calculations of the (Ag(2))(2), Ag(2)-H(2), and Ag(2)-C(6)H(6) complexes are performed using a coupled-cluster approach that includes single, double, and perturbative triple excitations (CCSD(T)), M?ller-Plesset second-order perturbation theory (MP2), and spin-component-scaled MP2 (SCS-MP2) methods. Using the atom pair approximation, the C(6) coefficients for silver-silver, silver-hydrogen, and silver-carbon atom systems are obtained after subtracting the energies of quadrupole-quadrupole interactions from the total electronic energy.     

Structure and stability of the organo-noble gas molecules XNgCCX and XNgCCNgX (Ng = Kr, Ar; X = F, Cl)

Density functional theory (B3LYP) and ab initio theory [second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster theory including single, double, and quasiperturbative triple excitations (CCSD(T))] have been used in combination with the standard and augmented correlation consistent basis sets (cc-pVnZ and aug-cc-pVnZ, where n = D, T, and Q) to investigate potential new noble gas compounds. Two classes of molecules were studied: XNgCCNgX and XNgCCX, where Ng = Kr and Ar and X = F and Cl. These molecules were characterized by finding the ground-state structures and calculating the relative energies, charge distributions, and vibrational frequencies. In addition, transition-state structures were also determined and decomposition pathways were identified through intrinsic reaction coordinate calculations.     

Is spin-component scaled second-order Møller-Plesset perturbation theory an appropriate method for the study of noncovalent interactions in molecules?

Testing of the spin-component scaled second-order M?ller-Plesset (SCS-MP2) method for the computation of noncovalent interaction energies is done with a database of 165 biologically relevant complexes. The effects of the spin-scaling procedure (i.e., MP2 vs SCS-MP2), the basis set size, and the corrections for basis set superposition error (BSSE) are systematically examined. When using two-point basis set extrapolations for the correlation energy, augmentation of the atomic orbital basis with computationally costly diffuse functions is found to be obsolete. In general, SCS-MP2 also improves results for noncovalent interactions statistically on MP2, and significant outliers are removed. Moreover, it is shown that effects of BSSE and one-particle basis set incompleteness almost cancel each other in the case of triple-zeta sets (SCS-MP2/TZVPP or SCS-MP2/cc-pVTZ without counterpoise correction), which opens a practical route to efficient computations for large systems. We recommend SCS-MP2 as the preferred quantum chemical wave function based method for the noncovalent interactions in large biologically relevant systems when reasonable coupled-cluster with single and double and perturbative triple excitations (CCSD(T)) calculations cannot be performed anymore. A comparison to MP2 and CCSD(T) interaction energies for n-alkane dimers, however, indicates (and this also holds to a lesser extent for hydrogen-bonded systems) limitations of SCS-MP2 when treating chemically "saturated" interactions. The different behavior of second-order perturbation theory for saturated and for stacked pi-systems is discussed.     

Benchmarking the performance of spin-component scaled CC2 in ground and electronically excited states

A generalization of the spin-component scaling and scaled opposite-spin modifications of second-order M?ller-Plesset perturbation theory to the approximate coupled-cluster singles-and-doubles model CC2 (termed SCS-CC2 and SOS-CC2) is discussed and a preliminary implementation of ground and excited state energies and analytic gradients is reported. The computational results for bond distances, harmonic frequencies, adiabatic and 0-0 excitation energies are compared with experimental results to benchmark their performance. It is found that both variants of the spin-scaling increase the robustness of CC2 against strong correlation effects and lead for this method even to somewhat larger improvements than those observed for second-order M?ller-Plesset perturbation theory. The spin-component scaling also enhances systematically the accuracy of CC2 for 0-0 excitation energies for pi --> pi* and n --> pi* transitions, if geometries are determined at the same level.     

Theoretical studies on hydroquinone-benzene clusters

High-level ab initio calculations were carried out to evaluate the interaction between the hydroquinone and benzene molecules. The intermolecular interaction energy was calculated using the M?ller-Plesset second-order perturbation theory at the complete basis set limit and also at the coupled cluster theory with single, double, and perturbatively triple excitations. The calculated binding energy is larger than the benzene dimer interaction energy. The T-shaped cluster (T-a) and the parallel conformation (P-a) are calculated to be nearly isoenergetic. Owing to the large energy gain in the attraction by electron correlation, the dispersion interaction is important for the attraction.     

How accurate is the density functional theory combined with symmetry-adapted perturbation theory approach for CH-pi and pi-pi interactions? A comparison to supermolecular calculations for the acetylene-benzene dimer

Five different orientations of the acetylene-benzene dimer including the T-shaped global minimum structure are used to assess the accuracy of the density functional theory combined with symmetry adapted perturbation theory (DFT-SAPT) approach in its density-fitting implementation (DF-DFT-SAPT) for the study of CH-pi and pi-pi interactions. The results are compared with the outcome of counterpoise corrected supermolecular calculations employing second-order M?ller-Plesset (MP2), spin-component scaled MP2 (SCS-MP2) and single and double excitation coupled cluster theory including perturbative triple excitations (CCSD(T)). For all considered orientations MP2 predicts much deeper potential energy curves with considerably shifted minima compared to CCSD(T) and DFT-SAPT. In spite of being an improvement over the results of MP2, SCS-MP2 tends to underestimate the well depth while DFT-SAPT, employing an asymptotically corrected hybrid exchange-correlation potential in conjunction with the adiabatic local density approximation for the exchange-correlation kernel, is found to be in excellent agreement with CCSD(T). Furthermore, DFT-SAPT provides a detailed understanding of the importance of the electrostatic, induction and dispersion contributions to the total interaction energy and their repulsive exchange corrections.     

Ab initio structural and vibrational investigation of sulfuric acid monohydrate

We employ ab initio methods to find stable geometries and to calculate potential energy surfaces and vibrational wavenumbers for sulfuric acid monohydrate. Geometry optimizations are carried out with the explicitly correlated coupled-cluster approach that includes single, double, and perturbative triple excitations (CCSD(T)-F12a) with a valence double-ζ basis set (VDZ-F12). Four different stable geometries are found, and the two lowest are within 0.41 kJ mol(-1) (or 34 cm(-1)) of each other. Vibrational harmonic wavenumbers are calculated at both the density-fitted local spin component scaled second-order M?ller-Plesset perturbation theory (DF-SCS-LMP2) with the aug-cc-pV(T+d)Z basis set and the CCSD-F12/VDZ-F12 level. Water O-H stretching vibrations and two highly anharmonic large-amplitude motions connecting the three lowest potential energy minima are considered by limiting the dimensionality of the corresponding potential energy surfaces to small two- or three-dimensional subspaces that contain only strongly coupled vibrational degrees of freedom. In these anharmonic domains, the vibrational problem is solved variationally using potential energy surfaces calculated at the CCSD(T)-F12a/VDZ-F12 level.     

Dissolution nature of the lithium hydroxide by water molecules

The structures, stabilities, thermodynamic quantities, dissociation energies, infrared spectra, and electronic properties of LiOH hydrated by up to seven water molecules are investigated by using the density-functional theory and the M?ller-Plesset second-order perturbation theory (MP2). Further accurate analysis based on the coupled-cluster theory with singles, doubles, and perturbative triples excitations agrees with the MP2 results. The Li-OH stretch mode significantly shifts with the increase of water molecules, and it eventually disappears upon dissociation. It is revealed that seven water molecules are needed for the stable dissociation of LiOH (as a completely dissociated conformation), in contrast to the cases of RbOH and CsOH which require four and three water molecules, respectively.     

An ab initio theoretical prediction: an antiaromatic ring pi-dihydrogen bond accompanied by two secondary interactions in a "wheel with a pair of pedals" shaped complex FH . . . C4H4 . . . HF

By the counterpoise-correlated potential energy surface method (interaction energy optimization), the structure of the pi H-bond complex FH cdots, three dots, centered FH . . . C4H4 . . . HF has been obtained at the second-order M?ller-Plesset perturbation theory (MP2/aug-cc-pVDZ) level. Intermolecular interaction energy of the complex is calculated to be -7.8 kcal/mol at the coupled-cluster theory with single, double substitutions and perturbatively linked triple excitations CCSD (T)/aug-cc-pVDZ level. The optimized structure is a "wheel with a pair of pedals" shaped (1mid R:1) structure in which both HF molecules almost lie on either vertical line passing through the middle-point of the C[Double Bond]C bond on either side of the horizontal plane of the C4 ring for cyclobutadiene. In the structure, an antiaromatic ring pi-dihydrogen bond is found, in which the proton acceptor is antiaromatic 4 electron and 4 center pi bond and the donors are both acidic H atoms of HF molecules. In accompanying with the pi-dihydrogen bond, two secondary interactions are exposed. The first is a repulsive interaction between an H atom of HF and a near pair of H atoms of C4H4 ring. The second is the double pi-type H bond between two lone pairs on a F atom and a far pair of H atoms.     

Extending the molecular size in accurate quantum-chemical calculations: the equilibrium structure and spectroscopic properties of uracil

The equilibrium structure of uracil has been investigated using both theoretical and experimental data. With respect to the former, quantum-chemical calculations at the coupled-cluster level in conjunction with a triple-zeta basis set have been carried out. Extrapolation to the basis set limit, performed employing the second-order M?ller-Plesset perturbation theory, and inclusion of core-correlation and diffuse-function corrections have also been considered. Based on the available rotational constants for various isotopic species together with corresponding computed vibrational corrections, the semi-experimental equilibrium structure of uracil has been determined for the first time. Theoretical and semi-experimental structures have been found in remarkably good agreement, thus pointing out the limitations of previous experimental determinations. Molecular and spectroscopic properties of uracil have then been studied by means of the composite computational approach introduced for the molecular structure evaluation. Among the results achieved, we mention the revision of the dipole moment. On the whole, it has been proved that the computational procedure presented is able to provide parameters with the proper accuracy to support experimental investigations of large molecules of biological interest.     

Model of peptide bond-aromatic ring interaction: correlated ab initio quantum chemical study

Aromatic ring-peptide bond interactions (modeled as benzene and formamide, N-methylformamide and N-methylacetamide) are studied by means of advanced computational chemistry methods: second-order M?ller-Plesset (MP2), coupled-cluster single and double excitation model [CCSD(T)], and density functional theory with dispersion (DFT-D). The geometrical preferences of these interactions as well as their interaction energy content, in both parallel and T-shaped arrangements, are investigated. The stabilization energy reaches a value of over 5 kcal/mol for the N-methylformamide-benzene complex at the CCSD(T)/complete basis set (CBS) level. Decomposition of interaction energy by the DFT-symmetry-adapted perturbation treatment (SAPT) technique shows that the parallel and T-shaped arrangements, although similar in their total interaction energies, differ significantly in the proportion of electrostatic and dispersion terms.     

Why benchmark-quality computations are needed to reproduce 1-adamantyl cation NMR chemical shifts accurately

While the experimental (1)H NMR chemical shiftsof the 1-adamantyl cation can be computed within reasonably small error bounds, the usual Hartree-Fock and density functional quantum-chemical computations, as well as those based on rather elaborate second-order M?ller-Plesset perturbation theory, fail to reproduce its experimental (13)C NMR chemical shifts satisfactorily. This also is true even if the NMR shielding calculations treat electron correlation adequately by the coupled-cluster singles and doubles model augmented by a perturbative correction for triple excitations (i.e., at the CCSD(T) level) with quadruple-ζ basis sets. We demonstrate that good agreement can be achieved if highly accurate 1-adamantyl cation equilibrium geometries based on parallel computations of CCSD(T) gradients are employed for the NMR shielding computations.