On the spectroscopic and thermochemical properties of ClO, BrO, IO, and their anions

A coupled cluster composite approach has been used to accurately determine the spectroscopic constants, bond dissociation energies, and heats of formation for the X1(2)II(3/2) states of the halogen oxides ClO, BrO, and IO, as well as their negative ions ClO-, BrO-, and IO-. After determining the frozen core, complete basis set (CBS) limit CCSD(T) values, corrections were added for core-valence correlation, relativistic effects (scalar and spin-orbit), the pseudopotential approximation (BrO and IO), iterative connected triple excitations (CCSDT), and iterative quadruples (CCSDTQ). The final ab initio equilibrium bond lengths and harmonic frequencies for ClO and BrO differ from their accurate experimental values by an average of just 0.0005 A and 0.8 cm-1, respectively. The bond length of IO is overestimated by 0.0047 A, presumably due to an underestimation of molecular spin-orbit coupling effects. Spectroscopic constants for the spin-orbit excited X2(2)III(1/2) states are also reported for each species. The predicted bond lengths and harmonic frequencies for the closed-shell anions are expected to be accurate to within about 0.001 A and 2 cm-1, respectively. The dissociation energies of the radicals have been determined by both direct calculation and through use of negative ion thermochemical cycles, which made use of a small amount of accurate experimental data. The resulting values of D0, 63.5, 55.8, and 54.2 kcal/mol for ClO, BrO, and IO, respectively, are the most accurate ab initio values to date, and those for ClO and BrO differ from their experimental values by just 0.1 kcal/mol. These dissociation energies lead to heats of formation, DeltaH(f) (298 K), of 24.2 +/- 0.3, 29.6 +/- 0.4, and 29.9 +/- 0.6 kcal/mol for ClO, BrO, and IO, respectively. Also, the final calculated electron affinities are all within 0.2 kcal/mol of their experimental values. Improved pseudopotential parameters for the iodine atom are also reported, together with revised correlation consistent basis sets for this atom.     

Ab initio characterization of XH3 (X=N, P). I. Ammonia, phosphine and their related ions and radicals: structure and thermochemistry

Ammonia, phosphine and their related cations, anions and radicals have been investigated at a high level of accuracy. The singles and doubles coupled cluster method including a perturbational correction for connected triple excitations, CCSD(T), in conjunction with correlation consistent basis sets ranging in size from triple to sextuple zeta have been employed. Extrapolation to the complete basis set limit has been used with accurate treatments of core–valence correlation effects in order to accurately predict structures, ionization potentials, electron affinities as well as N–H and P–H bond dissociation energies. For all the species studied, harmonic vibrational frequencies have also been evaluated in order to obtain zero-point corrections.     

Anharmonic frequencies of CX2Y2 (X, Y = O, N, F, H, D) isomers and related systems obtained from vibrational multiconfiguration self-consistent field theory

Accurate anharmonic frequencies are provided for molecules of current research, i.e., diazirines, diazomethane, the corresponding fluorinated and deuterated compounds, their dioxygen analogs, and others. Vibrational-state energies were obtained from state-specific vibrational multiconfiguration self-consistent field theory (VMCSCF) based on multilevel potential energy surfaces (PES) generated from explicitly correlated coupled cluster, CCSD(T)-F12a, and double-hybrid density functional calculations, B2PLYP. To accelerate the vibrational structure calculations, a configuration selection scheme as well as a polynomial representation of the PES have been exploited. Because experimental data are scarce for these systems, many calculated frequencies of this study are predictions and may guide experiments to come.     

The low-lying electronic excited states of NiCO

Highly correlated coupled cluster methods with single and double excitations (CSSD) and CCSD with perturbative triple excitations were used to predict molecular structures and harmonic vibrational frequencies for the electronic ground state X 1Sigma+, and for the 3Delta, 3Sigma+, 3Phi, 1 3Pi, 2 3Pi, 1Sigma+, 1Delta, and 1Pi excited states of NiCO. The X 1Sigma+ ground state's geometry is for the first time compared with the recently determined experimental structure. The adiabatic excitation energies, vertical excitation energies, and dissociation energies of these excited states are predicted. The importance of pi and sigma bonding for the Ni-C bond is discussed based on the structures of excited states.     

Calculation of the O-H stretching vibrational overtone spectrum of the water dimer

The O-H stretching vibrational overtone spectrum of the water dimer has been calculated with the dimer modeled as two individually vibrating monomer units. Vibrational term values and absorption intensities have been obtained variationally with a computed dipole moment surface and an internal coordinate Hamiltonian, which consists of exact kinetic energy operators within the Born-Oppenheimer approximation of the monomer units. Three-dimensional ab initio potential energy and dipole moment surfaces have been calculated using the internal coordinates of the monomer units using the coupled cluster method including single, double, and perturbative triple excitations [CCSD(T)] with the augmented correlation consistent valence triple zeta basis set (aug-cc-pVTZ). The augmented correlation consistent valence quadruple zeta basis set (aug-cc-pVQZ), counterpoise correction, basis set extrapolation to the complete basis set limit, relativistic corrections, and core and valence electron correlations effects have been included in one-dimensional potential energy surface cuts. The aim is both to investigate the level of ab initio and vibrational calculations necessary to produce accurate results when compared with experiment and to aid the detection of the water dimer under atmospheric conditions.     

A high‐accuracy theoretical study of the CHnP Systems n = 1–3

We have performed high‐level electronic structure computations on the most important species of the CH_nP systems n = 1–3 to characterize them and provide reliable information about the equilibrium and vibrationally averaged molecular structures, rotational constants, vibrational frequencies (harmonic and anharmonic), formation enthalpies, and vertical excitation energies. Those chemical systems are intermediates for several important reactions and also prototypical phosphorus‐carbon compounds; however, they are often elusive to experimental detection. The present results significantly complement their knowledge and can be used as an assessment of the experimental information when available. The explicitly correlated coupled‐cluster RCCSD(T)‐F12 method has been used for geometry optimizations and vibrational frequency calculations. Vibrational configuration interaction theory has been used to account for anharmonicity effects. Basis‐set limit extrapolations have been carried out to determine accurate thermochemical quantities. Electronic excited states have been calculated with coupled‐cluster approaches and also by means of the multireference configuration interaction method. © 2013 Wiley Periodicals, Inc.     

Accurate Quantum‐Chemical Prediction of Enthalpies of Formation of Small Molecules in the Gas Phase

The coupled‐cluster approach, including single and double excitations and perturbative corrections for triple excitations, is capable of predicting molecular electronic energies and enthalpies of formation of small molecules in the gas phase with very high accuracy (specifically, with error bars less than 5 kJ mol^?1), provided that the electronic wavefunction is dominated by the Hartree–Fock configuration. This capability is illustrated by calculations on molecules containing O–H and O–F bonds, namely OH, FO, H₂O, HOF, and F₂O. To achieve this very high accuracy, it is imperative to account for electron‐correlation effects in a quantitative manner, either by using explicitly correlated two‐particle basis functions (R12 functions) or by extrapolating to the limit of a complete basis. Besides taking into account harmonic zero‐point vibrational energies, it is also necessary to account for anharmonic corrections to the zero‐point vibrational energies, to include the core orbitals into the coupled‐cluster calculations, and to account for spin–orbit corrections and scalar relativistic effects. These additional corrections constitute small but significant contributions in the range of 1–4 kJ mol^?1 to the enthalpies of formation of the aforementioned molecules. The highly accurate coupled‐cluster results, obtained by employing R12 functions and by including various corrections, are compared with standard Kohn–Sham density‐functional calculations as well as with the Gaussian‐2 and complete‐basis‐set model chemistries.     

Structure and bonding of the MCN molecules, M=Cu,Ag,Au,Rg

High-precision calculations are reported for the title series with M=Cu, Ag, Au, using CCSD(T) with the latest pseudopotentials and basis sets up to cc-pVQZ. The bond lengths for M=Cu, Ag, Au agree with experiment within better than 1 pm. The role of deep-core excitations is studied. The second-order spin-orbit effects are evaluated at the density functional theory level, including M=Rg. A qualitative bonding analysis suggests multiple M-C bonding. The calculated vibrational frequencies are expected to be more accurate than the available experimental estimates. The M-C bond lengths obey Cu相似文献   

Rotational spectrum, potential energy surface, and bound states of the weakly bound complex He-N2O

Pure rotational transitions of the weakly bound complex He-N(2)O and three minor isotopomers (He-(14)N(15)NO, He-(15)N(14)NO, and He-(15)N(15)NO) were measured in the frequency region from 6 to 20 GHz. Predictions for the microwave transition frequencies were based on the infrared work by Tang and McKellar [J. Chem. Phys. 117, 2586 (2002)]. In the case of (14)N containing isotopomers, nuclear quadrupole hyperfine structure of the rotational transitions was observed and analyzed. The resulting spectroscopic parameters were used to determine geometrical and dynamical information about the complex. An ab initio potential energy surface was calculated at the coupled cluster level of theory with single and double excitations and perturbative inclusion of triple excitations. This surface was constructed using the augmented correlation consistent polarized valence triple zeta basis set for all atoms with the inclusion of bond functions for the van der Waals bond. Bound state calculations were done to determine the energies of low-lying rovibrational levels that are supported by the potential energy surface. The resulting transition energies agree with the experimental values to 1% or better.     

Combustion chemistry: important features of the C3H5 potential energy surface, including allyl radical, propargyl + H2, allene + H, and eight transition states

The C(3)H(5) potential energy surface (PES) encompasses molecules of great significance to hydrocarbon combustion, including the resonantly stabilized free radicals propargyl (plus H(2)) and allyl. In this work, we investigate the interconversions that take place on this PES using high level coupled cluster methodology. Accurate geometries are obtained using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)] combined with Dunning's correlation consistent quadruple-ζ basis set cc-pVQZ. The energies for these stationary points are then refined by a systematic series of computations, within the focal point scheme, using the cc-pVXZ (X = D, T, Q, 5, 6) basis sets and correlation treatments as extensive as coupled cluster with full single, double, and triple excitation and perturbative quadruple excitations [CCSDT(Q)]. Our benchmarks provide a zero-point vibrational energy (ZPVE) corrected barrier of 10.0 kcal mol(-1) for conversion of allene + H to propargyl + H(2). We also find that the barrier for H addition to a terminal carbon atom in allene leading to propenyl is 1.8 kcal mol(-1) lower than that for the addition to a central atom to form the allyl radical.     

The structures of m-benzyne and tetrafluoro-m-benzyne

The structures of m-benzyne and its fluorinated derivative, tetrafluoro-m-benzyne, were investigated using coupled cluster methods including triple excitations [CCSD(T) and CCSDT], different reference wave functions (spin-restricted Hartree-Fock, spin-unrestricted Hartree-Fock, and Brueckner), and different basis sets [6-31G(d,p) and correlation-consistent valence triple-zeta (cc-pVTZ)]. The inclusion of triple excitations in conjunction with d- and f-type polarization functions is paramount to correctly describe through-bond delocalization of the monocyclic form. At the highest level of theory, the C1-C3 distance of the minimum energy form of m-benzyne is 2.0 A and the profile of the potential energy surface along the C1-C3 distance is that of an asymmetric, single well, in agreement with previous density-functional theory and coupled cluster studies. In addition, the calculated CCSD(T) fundamental frequencies are in excellent agreement with the measured infrared frequencies, thus confirming the monocyclic form of m-benzyne. For tetrafluoro-m-benzyne, however, the increased eclipsing strain between the ring-external C-X bonds stabilizes the bicyclo[3.1.0]hexatriene form: the C1-C3 distance is calculated at the CCSD(T)/cc-pVTZ level to be approximately 1.75 A, which is in the range of elongated CC bonds. Computed harmonic vibrational frequencies compare reasonably well with the experimental neon-matrix difference spectrum and provide further evidence for the existence of a bicyclic form.     

Mixed uranium chloride fluorides UF6-nCln and methoxyuranium fluorides UF6-n(OCH3)n: a theoretical study of equilibrium geometries,vibrational frequencies,and the role of the f orbitals

The title compounds, the uranium (VI) fluoride chlorides (UF6-nCln, n = 0-6) and methoxyuranium (VI) fluorides [UF6-n(OCH3)n, n = 0-5], have been studied using relativistic density functional theory. Applying the B3LYP hybrid functional and an effective core potential on uranium, equilibrium geometries have been calculated for these molecules. In addition, harmonic vibrational frequencies have been computed for the chloride fluorides. Calculated frequencies have been compared to experiment where possible. All experimentally observed bands have been assigned, based on these calculations. The average deviation between theoretical and experimental frequencies is 15.6 cm-1 for 23 experimental modes. Theory always underestimates the experimental frequencies. This can be explained by the calculated bond lengths that are somewhat too long. The electronic structure of the uranium (VI) chloride fluorides has been investigated using scalar relativistic calculations and the PW91 functional. Periodic trends in the role and bonding contribution of the uranium 5f orbitals are discussed.     

Tailored coupled cluster singles and doubles method applied to calculations on molecular structure and harmonic vibrational frequencies of ozone

To assess the separation of dynamic and nondynamic correlations and orbital choice, we calculate the molecular structure and harmonic vibrational frequencies of ozone with the recently developed tailored coupled cluster singles and doubles method (TCCSD). We employ the Hartree-Fock and complete active space (CAS) self-consistent field (SCF) orbitals to perform TCCSD calculations. When using the Hartree-Fock orbitals, it is difficult to reproduce the experimental vibrational frequency of the asymmetric stretching mode. On the other hand, the TCCSD based on the CASSCF orbitals in a correlation consistent polarized valence triple zeta basis yields excellent results with the two symmetric vibrations differing from the experimental harmonic values by 2 cm(-1) and the asymmetric vibration differing by 9 cm(-1).     

Ab initio characterization of (CH3IO3) isomers and the CH3O2 + IO reaction pathways

The geometries, harmonic vibrational frequencies, relative energetics, and enthalpies of formation of (CH(3)IO(3)) isomers and the reaction CH(3)O(2) + IO have been investigated using quantum mechanical methods. Optimization has been performed at the MP2 level of theory, using all electron and effective core potential, ECP, computational techniques. The relative energetics has been studied by single-point calculations at the CCSD(T) level. Methyl iodate, CH(3)OIO(2), is found to be the lowest-energy isomer showing particular stabilization. The two nascent association minima, CH(3)OOOI and CH(3)OOIO, show similar stabilities, and they are considerably higher located than CH(3)OIO(2). Interisomerization barriers have been determined, along with the transition states involved in various pathways of the reaction CH(3)O(2) + IO.     

High-spin versus broken symmetry-Effect of DFT spin density representation on the geometries of three diiron (III) model compounds

Unrestricted density functional theory calculations have been conducted on three diiron(III) synthetic model compounds containing antiferromagnetically coupled high-spin (HS) irons for which crystallographic structures and Raman spectral data are available. Three density functionals have been employed: BPW91, PWC, and BOP. The study compares the effects on optimized geometries and harmonic vibrational frequencies of spin-paired (SP) low-spin, HS, and broken symmetry antiferromagnetically coupled singlet representations of the spin density distribution. The geometries around the diiron centers in the HS and broken symmetry (BS) representations are found to be similar, both markedly different from those arising from the SP representation. Small differences between the HS and BS results are seen in bond lengths, angles, Raman frequencies, and spin densities associated with oxo and peroxo bridges between the irons.     

The energetics and structural properties of diazomethyl (HCNN) and cyanomidyl (HNCN) radicals and their related cations and anions from ab initio calculations

The molecular structures and energetics of diazomethyl (HCNN) and cyanomidyl (HNCN) radicals and their related cations (HCNN(+),HNCN(+)) and anions (HCNN(-),HNCN(-)) are reported at a high level of accuracy. The singles and doubles coupled-cluster method including a perturbational correction for connected triple excitations with systematic sequences of correlation consistent basis sets have been employed. Extrapolation to the complete basis set limit has been used with accurate treatments of core-valence correlation effects in order to accurately predict molecular properties, ionization potentials, electron affinities as well as C-H and N-H bond dissociation energies. For all the species studied, harmonic vibrational frequencies have also been evaluated in order to obtain zero-point corrections to ionization potentials, electron affinities, and dissociation energies.     

Accurate ab initio determination of spectroscopic and thermochemical properties of mono- and dichlorocarbenes

The best technically feasible values for the triplet-singlet energy gap and the enthalpies of formation of the HCCl and CCl2 radicals have been determined through the focal-point approach. The electronic structure computations were based on high-level coupled cluster (CC) methods, up to quadruple excitations (CCSDTQ), and large-size correlation-consistent basis sets, ranging from aug-cc-pVDZ to aug-cc-pV6Z, followed by extrapolation to the complete basis set limit. Small corrections due to core correlation, relativistic effects, diagonal Born-Oppenheimer correction, as well as harmonic and anharmonic zero-point vibrational energy corrections have been taken into account. The final estimates for the triplet-singlet energy gap, T0(?), are 2170+/-40 cm-1 for HCCl and 7045+/-60 cm-1 for CCl2, favoring the singlet states in both cases. Complete quartic force fields in internal coordinates have been computed for both the X and ? states of both radicals at the frozen-core CCSD(T)/aug-cc-pVQZ level. Using these force fields vibrational energy levels of {HCCl, DCCl, CCl2} up to {6000, 5000, 7000} cm-1 were calculated both by second-order vibrational perturbation theory (VPT2) and variationally. These results, especially the variational ones, show excellent agreement with the experimentally determined energy levels. The enthalpies of formation of HCCl (X1A') and CCl2(X1A1), at 0 K, are 76.28+/-0.20 and 54.54+/-0.20 kcal mol-1, respectively.