The ab initio limit quartic force field of BH3

The complete quartic force field of BH(3) has been converged to the ab initio limit by extrapolation of core-valence correlation-consistent basis set series (cc-pCVXZ, X = T, Q, 5) of all-electron CCSD(T) (coupled-cluster singles and doubles with perturbative triples) energy points. Additional computations including full coupled-cluster treatments through quadruple excitations (CCSDTQ), scalar relativistic effects, and diagonal Born-Oppenheimer corrections (DBOC) were concurrently executed. Within second-order vibrational perturbation theory (VPT2) our quartic force field yields the fundamental frequencies nu(1) = 2502.3 cm(-1), nu(2) = 1147.2 cm(-1), nu(3) = 2602.1 cm(-1), and nu(4) = 1196.5 cm(-1), in excellent agreement with observed gas-phase fundamentals, displaying a mean absolute error of only 0.3 cm(-1). Our converged prediction for the equilibrium bond length of BH(3) is r(e) = 1.1867 A.     

Anharmonic force field of cis- and trans-formic acid from high-level ab initio calculations, and analysis of resonance polyads

The quadratic, cubic, and semidiagonal quartic force fields of cis- and trans-formic acid have been calculated using three different levels of theory. They all give satisfactory results, including the one at the lowest level of theory which is the MP2 method employing a basis set of triple-zeta quality. The results are used to theoretically analyze resonance polyads, including the one involving the 4(1), 5(1), 6(1), 8(1), 7(1)9(1), and 9(2) vibrational states. A semiexperimental equilibrium structure is derived from experimental ground state rotational constants and rovibrational interaction parameters calculated from the ab initio force field. The ab initio structure calculated at the CCSD(T) level of theory using a basis set of quintuple-zeta quality is in excellent agreement with the semiexperimental structure.     

SiH2Cl2: ab initio anharmonic force field, dipole moments, and infrared vibrational transitions

The vibrational spectra of SiH2Cl2 have been recorded in the 1000-13,000 cm(-1) region, utilizing the Fourier-transform spectroscopy and Fourier-transform intracavity laser absorption spectroscopy. Totally 61 band centers and intensities are derived from the infrared spectra. An ab initio quartic force field is obtained by applying the second-order Moller-Plesset perturbation theory and correlation-consistent polarized valence triplet-zeta basis sets [J. Chem. Phys. 90, 1007 (1989); 98, 1358 (1993)]. Most observed bands are assigned by the vibration analysis based on the second-order perturbation theory. Reduced-dimensional ab initio dipole moment functions (two dimensional and three dimensional) have also been calculated to investigate the absolute band intensities of the SiH2 chromophore. The calculated values agree reasonably with the observed ones.     

Ab initio characterization of C5

In this paper, the structure and spectroscopic parameters of the C5 cluster are determined using multiconfigurational quantum chemical methods as implemented in the MOLCAS software. A number of spectroscopic properties (band center positions, l-doubling parameters, and rotational constants) have been characterized. From the new results, the assignments of previous astrophysical observations [J. Goicoechea et al., Astrophys. J. 609, 225 (2004)] are discussed. A detailed exploration of the global potential energy surface confirms that C5 has a X1Sigmag+ linear isomer of prominent stability and, at least, three minimum energy structures showing singlet electronic ground states. Two of them are cyclic and one has a nonplanar geometry. Vertical and adiabatic electronic transitions and vibrational spectroscopic parameters are determined for the most stable linear isomer using multiconfigurational second order perturbation theory (CASPT2) using an active space containing 12 valence orbitals with 12 active electrons and extended ANO-type basis sets. The infrared spectrum has been analyzed from an anharmonic force field derived form the local surface, determined from the energies of a grid of 1350 geometries. The force field includes four coupling terms. The CASPT2 band center position of the nu7(piu) anharmonic fundamental has been calculated to be at 102 cm(-1), which validates the assignment to C5 of the pattern of bands centered at 102 cm(-1) observed with the ISO telescope.     

Study of correlation effects on stretching force constants of the H3N⋯LiF lithium-bonded and H3N⋯HF hydrogen-bonded complexes

All quadratic, cubic and quartic force constants associated with high and low vibrational modes of the H₃N⋯HF hydrogen-bonded and H₃N⋯LiF lithium-bonded complexes have been calculated employing the Møller—Plesset perturbation theory to the second order (MP2) with the 4-31G** basis set.     

Ab-initio molecular geometry and normal coordinate analysis of pyrrolidine molecule.

The molecular geometry of pyrrolidine was quantum mechanically calculated using the split valence 6-31G** basis set. Electron correlation energy has been computed employing MP2 method. The molecule showed an envelope form puckered structure with inter-plane angle of 36.4 degrees and has a total energy of -132976.80 kcal mol(-1) of which a -464.86 kcal mol(-1) electron correlation energy. The twist form of the molecule showed a twist angle of 10.2 degrees from planarity and has a total energy of -132976.05 kcal mol(-1) involving -464.097 kcal mol(-1) electron correlation energy. The normal coordinates of the molecule were theoretically analyzed on the basis of the Cs point symmetry of the envelope form. Using initial set of force constants obtained from the ab-initio calculations the fundamental vibrational frequencies were computed. The IR and laser Raman spectra of Pyrrolidine molecules were measured. All the observed vibrational bands including combination bands and overtones were assigned to normal modes with the aid of the potential energy distribution values obtained from normal coordinate calculations. The molecular force field was obtained by refining the initial set of force constants using the least square fit method. The molecular force field was determined by refining the initial set of force constants using the least square fit method instead of using the less accurate scaling factor methods. The determined molecular force field has produced simulated frequencies best match to the observed values. The low frequency molecular out-of-plane deformation modes were observed in both infrared and Raman spectra at 298 and 163 cm(-1). The barrier of ring twisting estimated from the observed ring out-of-plane vibrational mode at 163 cm(-1) was found 3.1 kcal mol(-1).     

Anharmonic vibrational frequencies and vibrationally averaged structures and nuclear magnetic resonance parameters of FHF-

The anharmonic vibrational frequencies of FHF(-) were computed by the vibrational self-consistent-field, configuration-interaction, and second-order perturbation methods with a multiresolution composite potential energy surface generated by the electronic coupled-cluster method with various basis sets. Anharmonic vibrational averaging was performed for the bond length and nuclear magnetic resonance indirect spin-spin coupling constants, where the latter computed by the equation-of-motion coupled-cluster method. The calculations placed the vibrational frequencies at 580 (nu(1)), 1292 (nu(2)), 1313 (nu(3)), 1837 (nu(1) + nu(3)), and 1864 cm(-1) (nu(1) + nu(2)), the zero-point H-F bond length (r(0)) at 1.1539 A, the zero-point one-bond spin-spin coupling constant [(1)J(0)(HF)] at 124 Hz, and the bond dissociation energy (D(0)) at 43.3 kcal/mol. They agreed excellently with the corresponding experimental values: nu(1) = 583 cm(-1), nu(2) = 1286 cm(-1), nu(3) = 1331 cm(-1), nu(1) + nu(3) = 1849 cm(-1), nu(1) + nu(2) = 1858 cm(-1), r(0) = 1.1522 A, (1)J(0)(HF) = 124+/-3 Hz, and D(0) = 44.4+/-1.6 kcal/mol. The vibrationally averaged bond lengths matched closely the experimental values of five excited vibrational states, furnishing a highly dependable basis for correct band assignments. An adiabatic separation of high- (nu(3)) and low-frequency (nu(1)) stretching modes was examined and found to explain semiquantitatively the appearance of a nu(1) progression on nu(3). Our calculations predicted a value of 186 Hz for experimentally inaccessible (2)J(0)(FF).     

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).     

Trans-1-chloro-2-fluoroethylene: microwave spectra and anharmonic force field

For the first time the millimeter-wave spectra of the trans-35ClHC=CHF and trans-37ClHC=CHF isotopomers have been observed in natural abundance. Many DeltaJ=0, +/-1 DeltaK(-1)=+1 transitions for 35ClHC=CHF and DeltaJ=0 DeltaK(-1)=+1 transitions for 37ClHC=CHF have been detected and assigned. This allowed us to accurately determine the vibrational ground-state rotational constants, quartic and some sextic centrifugal distortion constants, and nuclear quadrupole coupling constants for both 35Cl and 37Cl. The experimental investigation has been supported by highly accurate theoretical predictions. As far as ab initio computations are concerned, the complete set of cubic and quartic force constants have been evaluated by numerical differentiation of the analytic second-order M?ller-Plesset many-body perturbation theory/correlation consistent polarized valence triple zeta second derivatives. The anharmonic part of the force field completes the theoretical study on the equilibrium structure, dipole moment, chlorine quadrupolar tensor, and harmonic force field previously carried out by the same authors.     

Infrared overtone spectroscopy and vibrational analysis of a Fermi resonance in nitric acid: Experiment and theory

High resolution infrared spectra of nitric acid have been recorded in the first OH overtone region under jet-cooled conditions using a sequential IR-UV excitation method. Vibrational bands observed at 6933.39(3), 6938.75(4), and 6951.985(3) cm(-1) (origins) with relative intensities of 0.42(1), 0.38(1), and 0.20(1) are attributed to strongly mixed states involved in a Fermi resonance. A vibrational deperturbation analysis suggests that the optically bright OH overtone stretch (2nu1) at 6939.2(1) cm(-1) is coupled directly to the nu1 + 2nu2 state at 6946.4(1) cm(-1) and indirectly to the 3nu2 + nu3 + nu7 state at 6938.5(1) cm(-1). Both the identity of the zero-order states and the indirect coupling scheme are deduced from complementary CCSD(T) calculations in conjunction with second-order vibrational perturbation theory. The deperturbation analysis also yields the experimental coupling between 2nu1 and nu1 + 2nu2 of -6.9(1) cm(-1), and that between the two dark states of +5.0(1) cm(-1). The calculated vibrational energies and couplings are in near quantitative agreement with experimentally derived values except for a predicted twofold stronger coupling of 2nu1 to nu1 + 2nu2. Weaker coupling of the strongly mixed states to a dense background of vibrational states via intramolecular vibrational energy redistribution is evident from the experimental linewidths of 0.08 and 0.25 cm(-1) for the higher energy and two overlapping lower energy bands, respectively. A comprehensive rotational analysis of the higher energy band yields spectroscopic parameters and the direction of the OH overtone transition dipole moment.     

Potential-energy surface, van der Waals energy spectrum, and vibronic transitions in s-tetrazine-argon complex

The van der Waals vibrational states and the structure of the vibronic spectrum of s-tetrazine-argon complex have been studied by the ab initio methods. The potential-energy surface of the ground S(0) electronic state of the complex has been constructed by fitting the analytical many-body expansion to a large set of the interaction energy values computed using the second-order M?ller-Plesset perturbation theory combined with the standard aug-cc-pVDZ basis set. The equilibrium structure of the complex found is that with argon located above the tetrazine ring at a distance of 3.394 A. The calculated dissociation energy of 354 cm(-1) is compatible with the experiment. The van der Waals energy spectrum calculated from the potential-energy surface is explained analyzing a correlation with a simpler energy spectrum of benzene-argon. A new assignment of the S(0)-S(1) vibronic spectrum is proposed on the basis of the rigorous selection rules, vibrational energy levels in S(0) and S(1) states and vibronic transition intensities calculated from the electronic transition dipole moment surfaces.     

Vibrational anharmonicity and harmonic force fields for dichloromethane from quantum-chemical calculations

Anharmonic and related constants have been calculated for CH2Cl2, CD2Cl2, and CHDCl2 by using the program Gaussian03 and B3LYP and MP2 models. Bases used were 6-311++G** and cc-pVTZ. The size of grid used in the B3LYP/6-311++G** model had a noticeable effect on resulting data. Features of the MP2/6-311++G** calculations suggested a deleterious effect of the absence of f functions in this basis set. The need for the replacement of second-order terms in the perturbation theory formulas for the vibrational anharmonic constants x ij in the presence of Fermi resonance was explored, and minor resonances were found associated with the cubic constants varphi 122 and varphi 299 (d 0 isotopomer), phi122 and phi849 (d2), and phi278 (d1). Computed xij values for nuCH and nuCD motions agree quite well with earlier experimental data. Observed anharmonic frequencies, nu obsd, were corrected to "observed" harmonic frequencies, omega obsd, by using computed differences Delta = omegaQC-nuQC. These differences Delta are larger for the antisymmetric nuasCH2 mode than for symmetric nusCH2 motion. This fact made it necessary to use differing scale factors for the two kinds of CH stretching force constants in a subsequent scaling of the harmonic force field to nuobsd. Force field scaling was also carried out by refining to omega obsd. In both approaches, the B3LYP models required differing scale factors for symmetric and antisymmetric CCl stretching force constants, indicating a failure to compute an accurate C-Cl stretch-stretch interaction force constant. The MP2/cc-pVTZ force field was preferred. Both scaled and unscaled harmonic force fields were used to calculate centrifugal distortion constants (CDCs) and contributions to the vibrational dependence of the rotational constants (alphas). Variations in the CDCs can, in part, be explained by the magnitudes of the frequencies used in the scaling process.     

Infrared spectra of permanganate anions in potassium perchlorate matrices: vibrational anharmonicity,effective symmetry and vibrational mode mixing effects

Fourier transform infrared spectra of MnO4- anions isomorphously isolated in potassium perchlorate matrices were recorded at room and low temperature (LT, approximately 100 K). On the basis of the detected second-order vibrational transitions involving the dopant species nu3 mode components, anharmonicity constants and harmonic eigenvalues for these modes were calculated. Despite the fact that, rigorously speaking, the appearance of the spectra of dopant permanganate anions may be explained in terms of a Cs site symmetry, the nu3 stretching region resembles an approximate A1+E splitting (characteristic for a local C3nu or even higher symmetry), which is not expected, even within the latent symmetry approach. We explain such spectral patterns on the basis of vibrational mode mixing (a 'Fermi-like' resonance) of the MnO4- nu1 mode with the nu3, site-group component. With the results of degenerate case stationary perturbation theory, we show that in the present case the Fermi-like resonance is predominantly responsible for the observed spectral features. The appearance of the region of second-order vibrational transitions in the spectra of dopant permanganate anions may be better explained in terms of the (rigorous) crystallographic Cs site group (corresponding to the crystallographic Pnma space group), instead of the 'latent' (effective) symmetry site group C2nu (corresponding to the latent symmetry space group Imma).     

Vibrational spectra of stannane: Harmonic force field, Raman and IR intensities from ab initio correlated wavefunctions

The vibrational spectra and the harmonic force field for the stannane molecule has been calculated from ab initio SCF calculations using an effective core potential and two double-zeta basis sets for the valence electrons. Polarized functions and electron correlation effects on the calculated force constants and frequencies are analyzed carefully. At HF SCF level the calculated frequencies are higher than the experimental ones by about 8.15% while at second-order Møller-Plesset perturbation level, the mean absolute percentage deviation of the frequencies is found to be 3.75%. This mean absolute percentage deviation is notably improved at SDCI (3.2%) and SDQCI (2.2%) levels of theory. The integrated molar absorption coefficients and the scattering activities are calculated and their values agree reasonably with experiment.     

Ab-initio molecular geometry and normal coordinate analysis of tetrahydrothiophene molecule.

The molecular geometry of tetrahydrothiophene (THT) was quantum mechanically calculated using the split valence 6-31G** basis set. Electron correlation energy has been computed employing MP2 method. The molecule showed a twist form puckered structure with a twist torsion angle of 13 degrees and has a total energy of -347,877.514 kcal/mol of which a 436.715 kcal/mol electron correlation energy. The envelope form of the molecule showed an inter-plane angle of 22 degrees and has a total energy of -347,874.430 kcal/mol involving -436.558 kcal/mol electron correlation energy. The normal coordinates of the molecule were theoretically analyzed and the fundamental vibrational frequencies were calculated. The IR and laser Raman spectra of THT molecule was measured. All the observed vibrational bands including combination bands and overtones were assigned to normal modes with the aid of the potential energy distribution values obtained from normal coordinate calculations. The molecular force field was determined by refining the initial set of force constants using the least square fit method instead of using the less accurate scaling factor methods. The determined molecular force field has produced simulated frequencies which best match the observed values. The lowest-energy modes of vibration were two molecular out-of-plane deformations, observed at 114 and 166 cm(-1). The barrier of ring twisting estimated from the observed ring out-of-plane vibrational mode at 114 cm(-1) was estimated.     

Vibrational anharmonicity and vibrational Stark effect of sulfate ions trapped in potassium, rubidium and cesium chromate

Room and low temperature (approximately 100 K) FT-IR and Raman spectra of the sulfate doped K2CrO4, Rb2CrO4 and Cs2CrO4 were recorded. The positions of the nu1, nu3, and nu4 fundamental mode components of the dopant anions were measured. Nine (out of possible ten) second-order stretch-stretch vibrational transitions of the dopant anions were detected. On the basis of these data, the anharmonicity constants and the corresponding harmonic eigenvalues were calculated for several vibrational transitions using second-order perturbation theory expressions. The anharmonicity of the studied second-order transitions of the type nu1 + nu3i increases in the order (SO4/K2CrO4) < (SO4/Rb2CrO4) < (SO4/Cs2CrO4), while for those of the type nu3j + nu3i it basically follows the trend: (SO4/K2CrO4) > (SO4/Rb2CrO4) < (SO4/Cs2CrO4). The measured relative Stark splittings of the nu3 and nu4 mode components of the dopant SO4(2-) anions, as well as the average X13i/3i3j values decrease in the order (SO4/K2CrO4) > (SO4/Rb2CrO4) > (SO4/Cs2CrO4). In all cases, the splitting is larger for nu3 than for nu4 modes, indicating a smaller angular than the bond length distortion. The theory of vibrational Stark effect suggests that the observed frequency shifts of the nu1, nu3 and nu4 mode components may be attributed to the increase of the field strength at the doped anion site going from K2CrO4 to Cs2CrO4. The Stark splitting of these modes, on the other hand, implies that the internal crystalline field vector is almost parallel to the (hypothetical) C2 axis of the slightly distorted dopant tetrahedral anions.     

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.     

Anharmonic rovibrational analysis for disilacyclopropenylidene (Si2CH2)

The global minimum on the Si(2)CH(2) electronic singlet potential energy surface has been theoretically predicted to be a peculiar hydrogen bridged (Si···H···Si) disilacyclopropenylidene structure (Si(2)CH(2)). An accurate quartic force field for Si(2)CH(2) has been determined employing ab initio coupled-cluster theory with single and double excitations and a perturbative treatment for triple excitations [CCSD(T)], in combination with the correlation consistent core-valence quadruple zeta (cc-pCVQZ) basis set. The vibration-rotation coupling constants, equilibrium and zero-point vibration corrected rotational constants, centrifugal distortion constants, and harmonic and fundamental vibrational frequencies for six isotopologues of Si(2)CH(2) are predicted using vibrational second-order perturbation theory (VPT2). The anharmonic corrections for the vibrational motions involving the H bridged bonds are found to be more than 5% with respect to the corresponding harmonic vibrational frequencies. In this light, an experimental detection and characterization of disilacyclopropenylidene (Si(2)CH(2)) is highly desired.     

Equilibrium structure of sulfuric acid

The equilibrium structure of the more stable conformer of H2SO4, of C2 symmetry, has been calculated ab initio using the CCSD(T) method and taking into account the core correlation correction. The accuracy of this structure has been checked by comparing it to that of similar molecules and by estimating the effects of basis set enlargement and of diffuse functions. Furthermore, the quadratic, cubic, and quartic force fields have been calculated at the MP2 level of theory using a basis set of triple-zeta quality. The spectroscopic constants derived from the force field are in satisfactory agreement with the experimental ones. The resulting band origins are compared to literature infrared values, including those for overtone and combination bands. Normal modes of vibration are pictured. Using this force field, semiexperimental equilibrium rotational constants are determined which allows us to check the accuracy of the ab initio structure and to refine it using a mixed regression method.     

Quartic force field predictions of the fundamental vibrational frequencies and spectroscopic constants of the cations HOCO+ and DOCO+

Only one fundamental vibrational frequency of protonated carbon dioxide (HOCO(+)) has been experimentally observed in the gas phase: the ν(1) O-H stretch. Utilizing quartic force fields defined from CCSD(T)/aug-cc-pVXZ (X = T,Q,5) complete basis set limit extrapolated energies modified to include corrections for core correlation and scalar relativistic effects coupled to vibrational perturbation theory and vibrational configuration interaction computations, we are predicting the full set of gas phase fundamental vibrational frequencies of HOCO(+). Our prediction of ν(1) is within less than 1 cm(-1) of the experimental value. Our computations also include predictions of the gas phase fundamental vibrational frequencies of the deuterated form of the cation, DOCO(+). Additionally, other spectroscopic constants for both systems are reported as part of this study, and a search for a cis-HOCO(+) minimum found no such stationary point on the potential surface indicating that only the trans isomer is stable.