Semispectroscopic and quantitative structure-property relationship estimates of the equilibrium and vibrationally averaged structure and dipole moment of 1-buten-3-yne

Systematic quantum chemical calculations have been performed to obtain precise estimates of the equilibrium and vibrationally averaged molecular structure and electric dipole moment of vinylacetylene (VA, 1-buten-3-yne). Anharmonic (cubic and semi-diagonal quartic) MP2/cc-pVTZ force fields in normal coordinates were computed to account for anharmonic vibrational effects, including zero-point contributions to the rotational constants and the electric dipole moment. A simultaneous weighted least-squares structural refinement was performed, resulting in the best semispectroscopic estimate of the re structure of VA. The refinement was based on experimentally measured ground-state rotational constants of two isotopologs of VA corrected to equilibrium values using MP2/cc-pVTZ vibration-rotation interaction constants and all-electron CCSD(T)/aug-cc-pVTZ structural constraints. The semispectroscopic re structure of VA agrees excellently with the high-level CCSD(T)/aug-cc-pVTZ ab initio structure. The most dependable, CCSD(T)/cc-pVQZ//CCSD(T)/aug-cc-pVTZ equilibrium electric dipole moment of VA, in D, is mua= 0.4088, mub= 0.0004, and muc= 0. The vibrationally corrected a-component of 0.4214 D is in excellent agreement with one of the available experimental values. The present analysis shows that mub is negligible even after vibrational correction. A simple quantitative structure-property relationship (QSPR) model resulted in a highly similar estimate, 0.45 D, for the electric dipole moment of VA.     

Ab initio potential energy and dipole moment surfaces of (H2O)2

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for the water dimer, (H2O)2. The CCSD(T)-PES is a very precise fit to 19,805 ab initio energies obtained with the coupled-cluster (CCSD(T)) method, using an aug-cc-pVTZ basis. The standard counterpoise correction was applied to approximately eliminate basis set superposition errors. The fit is based on an approach that incorporates the permutational symmetry of identical atoms [Huang, X.; Braams, B.; Bowman, J. M. J. Chem.Phys. 2005, 122, 044308]. The DMS is a fit to the dipole moment obtained with M?ller-Plesset (MP2) theory, using an aug-cc-pVTZ basis. The PES has an RMS fitting error of 31 cm(-1) for energies below 20,000 cm(-1), relative to the global minimum. This surface can describe various internal floppy motions, including various monomer inversions, and isomerization pathways. Ten characteristic stationary points have been located on the surface, four of which are transition structures and the rest are higher order saddle points. Their geometrical and vibrational properties are presented and compared with available previous theoretical work. The CCSD(T)-PES and MP2-DMS dissociate correctly (and symmetrically) to two H2O monomers, with D(e) = 1665.7 cm(-1) (19.93 kJ/mol). Accurate quantum calculations of the zero-point energy of the dimer (using diffusion Monte Carlo) and the monomers (using a vibrational configuration interaction approach) are reported, and these together with D(e) give a value of D0 of 1042 cm(-1) (12.44 kJ/mol). A best estimated value is 1130 cm(-1) (13.5 kJ/mol).     

Flexible, ab initio potential, and dipole moment surfaces for water. I. Tests and applications for clusters up to the 22-mer

We report full-dimensional, ab initio potential energy and dipole moment surfaces, denoted PES and DMS, respectively, for arbitrary numbers of water monomers. The PES is a sum of 1-, 2-, and 3-body potentials which can also be augmented by semiempirical long-range higher-body interactions. The 1-body potential is a spectroscopically accurate monomer potential, and the 2- and 3-body potentials are permutationally invariant fits to tens of thousands of CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ electronic energies, respectively. The DMS is a sum of 1- and 2-body DMS, which are covariant fits to tens of thousands MP2/aug-cc-pVTZ dipole moment data. We present the details of these new 2- and 3-body potentials and then extensive applications and tests of this PES are made to the structures, classical binding energies, and harmonic frequencies of water clusters up to the 22-mer. In addition, we report the dipole moment for these clusters at various minima and compare the results against available and new ab initio calculations.     

Ab initio potential energy and dipole moment surfaces for H5O2 +

Full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for H(5)O(2) (+). Tens of thousands of coupled-cluster [CCSD(T)] and second-order Moller-Plesset (MP2) calculations of electronic energies, using aug-cc-pVTZ basis, were done. The energies were fit very precisely in terms of all the internuclear distances, using standard least-square procedures, however, with a fitting basis that satisfies permutational symmetry with respect to like atoms. The H(5)O(2) (+) PES is a fit to 48 189 CCSD(T) energies, containing 7962 polynomial coefficients. The PES has a rms fitting error of 34.9 cm(-1) for the entire data set up to 110 000 cm(-1). This surface can describe various internal floppy motions, including the H atom exchanges, monomer inversions, and monomer torsions. First- and higher-order saddle points have been located on the surface and compared with available previous theoretical work. In addition, the PES dissociates correctly (and symmetrically) to H(2)O+H(3)O(+), with D(e)=11 923.8 cm(-1). Geometrical and vibrational properties of the monomer fragments are presented. The corresponding global DMS fit (MP2 based) involves 3844 polynomial coefficients and also dissociates correctly.     

A quantum chemistry study of the Cl atom reaction with formaldehyde

The elementary vapor-phase reaction between Cl atoms and HCHO has been studied by ab initio methods. Calculations at the MP2, MP3, MP4(SDTQ), CCSD, CCSD(T), and MRD-CI levels of theory show that the reaction is characterized by a low electronic barrier; excluding the effects of spin-orbit splitting in Cl, our best estimate at the MRD-CI/aug-cc-pVTZ//RHF-RCCSD(T)/aug-cc-pVTZ level of theory predicts a Born-Oppenheimer barrier height of 0.7 kJ mol-1. The energies of the lowest six electronic states as resulting from MRD-CI calculations are presented at discrete points along the reaction path, and two avoided crossings are found in the transition state region. The spin-orbit splitting in Cl is also calculated along the reaction path; it is not negligible in the transition state region and is found to increase the barrier by only 1.4 kJ mol-1 at the RCCSD(T)/aug-cc-pVTZ transition state geometry. The minimum energy path of the reaction connects an energetically weakly stabilized adduct on the flat potential surface on the reactant side and an energetically strongly stabilized postreaction adduct. The reaction rate coefficient and the kinetic isotope effects were calculated using improved canonical variational theory with small curvature tunneling (ICVT/SCT), and the results were compared to experimental data. The experimental reaction rate coefficient is reproduced within its uncertainty limits by variational transition state theory with interpolated single-point energy corrections (ISPE) at the MP4(SDTQ) level of theory and by conventional transition state theory with interpolated optimized energies (IOE) at the MRD-CI//RCCSD(T) level of theory and interpolated optimized geometries at the RCCSD(T) level of theory on an MP2/aug-cc-pVTZ potential energy surface when employing scaled vibrational frequencies.     

State of the art theoretical study and comparison to experiment for the phenol...argon complex

The phenol...argon complex was studied by means of various high level ab initio quantum mechanics methods and high resolution threshold ionization spectroscopy. The structure and stabilization energy of different conformers were determined. Stabilization energy of van der Waals bonded and H-bonded PhOH...Ar complex determined at CCSD(T) complete basis set (CBS) level for CP-RI-MP2/cc-pVTZ/Ar aug-cc-pVTZ geometries amount to 434 and 285 cm(-1). The CCSD(T)/CBS were constructed either as a sum of MP2/CBS interaction energy and CCSD(T) correction term [difference between CCSD(T) and MP2 correlation energies determined with medium basis set] or directly from CCSD(T)/aug-cc-pVDZ and aug-cc-pVTZ energies. Both schemes provide very similar values. Harmonic vibrational analysis revealed that the H-bonded structure does not represent energy minimum but first order transition structure. The respective imaginary vibrational mode (16 cm(-1)) connects two possible argon locations -- above and below the phenol aromatic ring. Including the DeltaZPVE, we obtained stabilization enthalpy at 0 K of 389 cm(-1). This value is marginally higher (25-35 cm(-1), 0.07-0.10 kcal/mol) than the experimental value. The determination of DeltaZPVE constitutes the most significant error and possible improvements should come from more accurate evaluation of the (nonharmonic) vibrational frequencies.     

The calculated infrared spectrum of Cl- H2O using a full dimensional ab initio potential surface and dipole moment surface

We report a full dimensional ab initio-based global potential energy surface (PES) and dipole moment surface (DMS) for Cl(-)H(2)O. Both surfaces are symmetric with respect to interchange of the H atoms. The PES is a fit to thousands of electronic energies calculated using the coupled-cluster method (CCSD(T)) with a moderately large basis (aug-cc-pVTZ). The infrared spectrum and vibrational dynamics are reported and compared to experiment. These results are analyzed by examination of wave function and the dipole surface.     

Structures and stability of medium silicon clusters. II. Ab initio molecular orbital calculations of Si12-Si20

Ab initio all-electron molecular-orbital calculations are carried out to study the structures and relative stability of low-energy silicon clusters (Si(n),n = 12-20). Selected geometric isomers include those predicted by Ho et al. [Nature (London) 392, 582 (1998)] based on an unbiased search with tight-binding/genetic algorithm, as well as those found by Rata et al. [Phys. Rev. Lett. 85, 546 (2000)] based on density-functional tight-binding/single-parent evolution algorithm. These geometric isomers are optimized at the M?ller-Plesset (MP2) MP2/6-31G(d) level. The single-point energy at the coupled-cluster single and double substitutions (including triple excitations) [CCSD(T)] CCSD(T)/6-31G(d) level for several low-lying isomers are further computed. Harmonic vibrational frequency analysis at the MP2/6-31G(d) level of theory is also undertaken to assure that the optimized geometries are stable. For Si12-Si17 and Si19 the isomer with the lowest-energy at the CCSD(T)/6-31G(d) level is the same as that predicted by Ho et al., whereas for Si18 and Si20, the same as predicted by Rata et al. However, for Si14 and Si15, the vibrational frequency analysis indicates that the isomer with the lowest CCSD(T)/6-31G(d) single-point energy gives rise to imaginary frequencies. Small structural perturbation onto the Si14 and Si15 isomers can remove the imaginary frequencies and results in new isomers with slightly lower MP2/6-31G(d) energy; however the new isomers have a higher single-point energy at the CCSD(T)/6-31G(d) level. For most Si(n) (n = 12-18,20) the low-lying isomers are prolate in shape, whereas for Si19 a spherical-like isomer is slightly lower in energy at the CCSD(T)/6-31G(d) level than low-lying prolate isomers.     

Electron affinity of Al13: a correlated electronic structure study

Neutral and anionic 13-atom aluminum clusters are studied with high-level, fully ab initio methods: second-order perturbation theory (MP2) and coupled cluster theory with singles, doubles, and perturbative triples (CCSD(T)). Energies and vibrational frequencies are reported for icosahedral and decahedral isomers, and are compared with density functional theory results. At the MP2 level of theory, with all of the basis sets employed, the icosahedral structure is energetically favored over the decahedral structure for both the neutral and anionic Al(13) clusters. Hessian calculations imply that only the icosahedral structures are potential energy minima. The CCSD(T)/aug-cc-pVTZ adiabatic electron affinity of Al(13) is found to be 3.57 eV, in excellent agreement with experiment.     

In search of definitive signatures of the elusive NCCO radical

Previous experimental assignments of the fundamental vibrational frequencies of NCCO have been brought into question by subsequent unsuccessful attempts to observe IR signatures of this radical at these frequencies. Here we compute the fundamental vibrational frequencies by applying second-order vibrational perturbation theory to the complete quartic force field computed at the all-electron (AE) coupled cluster singles, doubles, and perturbative triples level [CCSD(T)] with the correlation-consistent, polarized core-valence quadruple-zeta (cc-pCVQZ) basis set, which has tight functions to correctly describe core correlation. The AE-CCSD(T)/cc-pCVQZ geometric parameters are r(e)(N-C)=1.1623 A, r(e)(C-C)=1.4370 A, r(e)(C-O)=1.1758 A, theta(e)(N-C-C)=168.55 degrees , and theta(e)(C-C-O)=132.22 degrees . Our CCSD(T)/cc-pCVQZ values of the characteristic stretching frequencies nu(1) and nu(2) are 2171 and 1898 cm(-1), respectively, in stark contrast to the experimentally derived values of 2093 and 1774 cm(-1). Finally, focal-point extrapolations using correlation-consistent basis sets cc-pVXZ (X=D,T,Q,5,6) and electron correlation treatments as extensive as full coupled cluster singles, doubles, and triples (CCSDT) with perturbative accounting of quadruple excitations [CCSDT(Q)] determine the vibrationless barrier to linearity of NCCO and the dissociation energy (D(0)) of NCCO-->NC+CO to be 8.4 and 26.5 kcal mol(-1), respectively. Using our precisely determined dissociation energy, we recommend a new 0 K enthalpy of formation for NCCO of 50.9+/-0.3 kcal mol(-1).     

Ab initio global potential-energy surface for H5(+) --> H3(+) + H2

An accurate global potential-energy surface (PES) is reported for H5(+) based on more than 100,000 CCSD(T)/aug-cc-pVTZ ab initio energies. This PES has full permutational symmetry with respect to interchange of H atoms and dissociates to H3(+) and H2. Ten known stationary points of H5(+) are characterized and compared to previous ab initio calculations. Quantum diffusion Monte Carlo calculations are performed on the PES to obtain the zero-point energy of H5(+) and the anharmonic dissociation energy (D0) of H5(+) --> H3(+) + H2. The rigorous zero-point state of H4D+ is also calculated and discussed within the context of a strictly classical approach to obtain the branching ratio of the reaction H4D+ --> H3(+) + HD and H2D+ + H2. Such an approach is taken using the PES and critiqued based on the properties of the quantum zero-point state. Finally, a simple procedure for adding the long range-interaction energy is described.     

Aliphatic C-H/pi interactions: Methane-benzene, methane-phenol, and methane-indole complexes

Noncovalent C-H/pi interactions are prevalent in biochemistry and are important in molecular recognition. In this work, we present potential energy curves for methane-benzene, methane-phenol, and methane-indole complexes as prototypes for interactions between C-H bonds and the aromatic components of phenylalanine, tyrosine, and tryptophan. Second-order perturbation theory (MP2) is used in conjunction with the aug-cc-pVDZ and aug-cc-pVTZ basis sets to determine the counterpoise-corrected interaction energy for selected complex configurations. Using corrections for higher-order electron correlation determined with coupled-cluster theory through perturbative triples [CCSD(T)] in the aug-cc-pVDZ basis set, we estimate, through an additive approximation, results at the very accurate CCSD(T)/aug-cc-pVTZ level of theory. Symmetry-adapted perturbation theory (SAPT) is employed to determine the physically significant components of the total interaction energy for each complex.     

An ab initio based global potential energy surface describing CH5+ --> CH3+ + H2

A full-dimensional, ab initio based potential energy surface (PES) for CH(5)(+), which can describe dissociation is reported. The PES is a precise fit to 36173 coupled-cluster [CCSD(T)] calculations of electronic energies done using an aug-cc-pVTZ basis. The fit uses a polynomial basis that is invariant with respect to permutation of the five H atoms, and thus describes all 120 equivalent minima. The rms fitting error is 78.1 cm(-1) for the entire data set of energies up to 30,000 cm(-1) and a normal-mode analysis of CH(5)(+) also verifies the accuracy of the fit. Two saddle points have been located on the surface as well and compared with previous theoretical work. The PES dissociates correctly to the fragments CH(3)(+) + H(2) and the equilibrium geometry and normal-mode analyses of these fragments are also presented. Diffusion Monte Carlo calculations are done for the zero-point energies of CH(5)(+) (and some isotopologs) as well as for the separated fragments of CH(5)(+), CH(3)(+) + H(2) and those of CH(4)D(+), CH(3)(+) + HD and CH(2)D(+) + H(2). Values of D(0) are reported for these dissociations. A molecular dynamics calculation of CH(4)D(+) dissociation at one total energy is also performed to both validate the applicability of the PES for dynamics studies as well as to test a simple classical statistical prediction of the branching ratio of the dissociation products.     

Ab initio and analytic intermolecular potentials for Ar-CH(3)OH

Ab initio calculations at the CCSD(T)/aug-cc-pVTZ level of theory were used to characterize the Ar-CH(3)OH intermolecular potential energy surface (PES). Potential energy curves were calculated for four different Ar + CH(3)OH orientations and used to derive an analytic function for the intermolecular PES. A sum of Ar-C, Ar-O, Ar-H(C), and Ar-H(O) two-body potentials gives an excellent fit to these potential energy curves up to 100 kcal mol(-1), and adding an additional r(-n) term to the Buckingham two-body potential results in only a minor improvement in the fit. Three Ar-CH(3)OH van der Waals minima were found from the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ calculations. The structure of the global minimum is in overall good agreement with experiment (X.-C. Tan, L. Sun and R. L. Kuczkowski, J. Mol. Spectrosc., 1995, 171, 248). It is T-shaped with the hydroxyl H-atom syn with respect to Ar. Extrapolated to the complete basis set (CBS) limit, the global minimum has a well depth of 0.72 kcal mol(-1) with basis set superposition error (BSSE) correction. The aug-cc-pVTZ basis set gives a well depth only 0.10 kcal mol(-1) smaller than this value. The well depths of the other two minima are within 0.16 kcal mol(-1) of the global minimum. The analytic Ar-CH(3)OH intermolecular potential also identifies these three minima as the only van der Waals minima and the structures predicted by the analytic potential are similar to the ab initio structures. The analytic potential identifies the same global minimum and the predicted well depths for the minima are within 0.05 kcal mol(-1) of the ab initio values. Combining this Ar-CH(3)OH intermolecular potential with a potential for a OH-terminated alkylthiolate self-assembled monolayer surface (i.e., HO-SAM) provides a potential to model Ar + HO-SAM collisions.     

Kinetic isotope effects for Cl + CH4 ⇌ HCl + CH3 calculated using ab initio semiclassical transition state theory

Calculations were carried out for 25 isotopologues of the title reaction for various combinations of (35)Cl, (37)Cl, (12)C, (13)C, (14)C, H, and D. The computed rate constants are based on harmonic vibrational frequencies calculated at the CCSD(T)/aug-cc-pVTZ level of theory and X(ij) vibrational anharmonicity coefficients calculated at the CCSD(T) /aug-cc-pVDZ level of theory. For some reactions, anharmonicity coefficients were also computed at the CCSD(T)/aug-cc-pVTZ level of theory. The classical reaction barrier was taken from Eskola et al. [J. Phys. Chem. A 2008, 112, 7391-7401], who extrapolated CCSD(T) calculations to the complete basis set limit. Rate constants were calculated for temperatures from ～100 to ～2000 K. The computed ab initio rate constant for the normal isotopologue is in good agreement with experiments over the entire temperature range (～10% lower than the recommended experimental value at 298 K). The ab initio H/D kinetic isotope effects (KIEs) for CH(3)D, CH(2)D(2), CHD(3), and CD(4) are in very good agreement with literature experimental data. The ab initio (12)C/(13)C KIE is in error by ～2% at 298 K for calculations using X(ij) coefficients computed with the aug-cc-pVDZ basis set, but the error is reduced to ～1% when X(ij) coefficients computed with the larger aug-cc-pVTZ basis set are used. Systematic improvements appear to be possible. The present SCTST results are found to be more accurate than those from other theoretical calculations. Overall, this is a very promising method for computing ab initio kinetic isotope effects.     

Highly accurate potential-energy and dipole moment surfaces for vibrational state calculations of methane

Full-dimensional ab initio potential-energy surface (PES) and dipole moment surface are constructed for a methane molecule at the CCSD(T)/cc-pVTZ and MP2/cc-pVTZ levels of theory, respectively, by the modified Shepard interpolation method based on the fourth-order Taylor expansion [MSI(4th)]. The reference points for the interpolation have been set in the coupling region of CH symmetric and antisymmetric stretching modes so as to reproduce the vibrational energy levels related to CH stretching vibrations. The vibrational configuration-interaction calculations have been performed to obtain the energy levels and the absorption intensities up to 9000 cm(-1) with the use of MSI(4th)-PES. The calculated fundamental frequencies and low-lying vibrational energy levels show that MSI(4th) is superior to the widely employed quartic force field, giving a better agreement with the experimental values. The absorption bands of overtones as well as combination bands, which are caused by purely anharmonic effects, have been obtained up to 9000 cm(-1). Strongly coupled states with visible intensity have been found in the 6500-9000 cm(-1) region where the experimental data are still lacking.     

Ab initio based potential energy surfaces and Franck-Condon analysis of ionization thresholds of cyclic-C3H and linear-C3H

We report a Franck-Condon analysis in reduced dimensionality of the ionization thresholds of linear(l)-C3H and cyclic(c)-C3H using MP2-based potential energy surfaces and CCSD(T)/aug-cc-pVTZ calculations of electronic energies at selected geometries. The potential energy surfaces are fits to tens of thousands of MP2/aug-cc-pVTZ energies for the neutral and cation systems. These fits properly describe the invariance of the potential with respect to all permutations of the three C atoms. The realism of the potential surfaces is assessed by comparing stationary-point structures, energies, and normal-mode frequencies with previous high-level ab initio calculations. Several key vibrational modes in this ionization process are located at saddle points and so a numerical approach to obtain the Franck-Condon factors for those modes is done. On the basis of this analysis combined with a simple harmonic treatment of the energies of the remaining modes and key electronic energy differences obtained with CCSD(T)/aug-cc-pVTZ calculations, we find the threshold ionization energy of l-C3H to be 9.06 eV and for c-C3H we estimate the threshold to be in the range 9.70-9.76 eV. We estimate these values are accurate to within +/-0.05 eV.     

Dihydrogen bonding vs metal-σ interaction in complexes between H2 and metal hydride

The complexes formed by hydrogen with metal hydrides (LiH, NaH, BeH(2), MgH(2), BH(3), AlH(3), Li(2)H(2), Na(2)H(2), Be(2)H(4), and Mg(2)H(4)) have been theoretically studied at the MP2/aug-cc-pVTZ, MP2/aug-cc-pVQZ and CCSD(T)/aug-cc-pVTZ//CCSD/aug-cc-pVTZ levels of theory. The hydrogen molecule can act as a Lewis acid or base. In the first case, a dihydrogen bonded complex is obtained and in the second an interaction between the σ-bond of the hydrogen molecule and an empty orbital of the metal atoms is found. Quantum theory of atoms in molecules and natural bond orbitals methods have been applied to analyze the intermolecular interactions. Additionally, the cooperativity effects are analyzed for selected complexes with two H(2) molecules where both kinds of interactions exist simultaneously.     

Large-amplitude dynamics in vinyl radical: the role of quantum tunneling as an isomerization mechanism

We report tunneling splittings associated with the large amplitude 1,2 H-atom migration to the global minima in the vinyl radical. These are obtained using a recent full-dimensional ab initio potential energy surface (PES) [A. R. Sharma, B. J. Braams, S. Carter, B. C. Shepler, and J. M. Bowman, J. Chem. Phys. 130(17), 174301 (2009)] and independently, directly calculated "reaction paths." The PES is a multidimensional fit to coupled cluster single and double and perturbative treatment of triple excitations coupled-cluster single double triple (CCSD(T)) with the augmented correlation consistent triple zeta basis set (aug-cc-pVTZ). The reaction path potentials are obtained from a series of CCSD(T)/aug-cc-pVnTZ calculations extrapolated to the complete basis set limit. Approximate 1D calculations of the tunneling splitting for these 1,2-H atom migrations are obtained using each of these potentials as well as quite different 1D Hamiltonians. The splittings are calculated over a large energy ranges, with results from the two sets of calculations in excellent agreement. Though negligibly slow (>1 s) for the vibrational ground state, this work predicts tunneling-promoted 1,2 hydride shift dynamics in vinyl to exhibit exponential growth with internal vibrational excitation, specifically achieving rates on the sub-μs time scale at energies above E ≈ 7500 cm(-1). Most importantly, these results begin to elucidate the possible role of quantum isomerization through barriers without dissociation, in competition with the more conventional picture of classical roaming permitted over a much narrower window of energies immediately below the bond dissociation limit. Furthermore, when integrated over a Boltzmann distribution of thermal energies, these microcanonical tunneling rates are consistent with sub-μs time scales for 1,2 hydride shift dynamics at T > 1400 K. These results have potential relevance for combustion modeling of low-pressure flames, as well as recent observations of nuclear spin statistical mixing from high-resolution IR/microwave spectroscopy on vinyl radical.