An ab initio study of the electronic structure of BCl3(2+) and its decomposition pathways

Ab initio methods have been used to characterise the lowest energy potential energy surfaces of 1BCl3(2+) and 3BCl3(2+). The methods used are MP2, CCSD, CCSD(T) and MCSCF with 6-311G(d), cc-pVTZ and aug-cc-pVTZ basis sets. While the singlet surface is relatively straight-forward, the triplet surface is very complicated, with many stationary points in close energetic proximity. The singlet surface can fragment to the following products (1BCl + 1Cl+ + 1Cl+), (1Cl+ + 1B+ + 1Cl2), (2BCl+ + 2Cl), while the triplet surface can fragment to (1BCl2+ + 3Cl+) and (2BCl2(2+) + 2Cl). 2BCl2(2+) can further fragment to (1Cl+ + 2BCl+). These results are in good agreement with previous experimental data from coincidence mass spectroscopy. [symbol: see text]1 diagnostic values have been calculated for all of the stationary points of BCl3(2+), using the method of Lee and Taylor. These data, together with CCSD/CCSD(T) energy differences and MCSCF calculations, are used to conclude that most of the stationary points on the singlet surface are well represented using single reference methods. The stationary points of the triplet system have [symbol: see text]1 diagnostic values greater than those for the singlet system, as expected when using the closed-shell [symbol: see text]1 diagnostic method of Lee and Taylor. However, all of the structures have acceptable single reference wavefunctions if the open-shell upper limit of Rienstra-Kiracofe et al. (0.045) is used, a conclusion fully supported by CCSD/CCSD(T) energy differences. CCSD(T) energies determined for the fragmentation asymptotes have been compared with experimental data collated from the NIST Theoretical and Scientific Data website, and the generally very good agreement between theory and experiment reinforces the reliability of the CCSD(T) method.     

The importance of dihydrogen complexes HnGe(H2)+ (n=0,1) to the chemistry of cationic germanium hydrides: advanced theoretical and mass spectrometric analysis

Investigations of [Ge,Hn]-/0/- (n = 2,3) have been performed using a four-sector mass spectrometer. The results reveal that the complexes HnGe(H2)+ (n = 0,1) play an important role in the unimolecular dissociation of the metastable cations. Theoretical calculations support the experimental observations in most instances, and the established view that the global minimum of [Ge,H2]+ is an inserted structure may need reexamination; CCSD(T,full)/cc-pVTZ//CCSD(T)/6-311 ++ G(d,p) and B3LYP/cc-pVTZ studies of three low-lying cation states (2A1 HGeH+, 2B2 Ge(H2)+ and 2B1 Ge(H2)+) indicate a very small energy difference (ca. 4 kcal mol(-1)) between 2A1 HGeH+ and 2B2 Ge(H2)+; B3LYP favours the ion-molecule complex, whereas coupled-cluster calculations favour the inserted structure for the global minimum. Single-point multireference (MR) averaged coupled-pair functional and MR-configuration interaction calculations give conflicting results regarding the global minimum. We also present theoretical evidence indicating that the orbital-crossing point implicated in the spin-allowed metastable dissociation HGeH+* --> Ge(H2)+* --> Ge+ + H2 lies above the H-loss asymptote. Thus, a quantum-mechanical tunneling mechanism is invoked to explain the preponderance of the H2-loss signal for the metastable ion.     

Computational studies on the ground and excited states of BrOOBr

In this work, theoretical computations for the ground and excited states of BrOOBr have been performed at high-level ab initio molecular orbital theories. The ground-state geometries of BrOOBr in different forms (trans, cis, and twist form) have been optimized at the couple-cluster CCSD(T) level of theory with cc-pVTZ and aug-cc-pVTZ basis sets, which indicates that at CCSD(T)/cc-pVTZ level of theory, the twist form is 4.96 kcal/mol more stable than the trans form and 10.67 kcal/mol more stable than the cis form; at the CCSD(T)/aug-cc-pVTZ basis set the twist form is 4.33 kcal/mol more stable than the trans form and 9.54 kcal/mol more stable than the cis form. The vertical excitation energies and potential-energy curves for the singlet and triplet low-lying excited states of BrOOBr were calculated at both the complete active space self-consistent-field (CASSCF) level of theory and the multireference internally contracted configuration interaction (MRCI) level of theory. The differences of potential-energy curves at CASSCF and MRCI levels of theory are found for the BrOOBr excited states. At CASSCF level of theory, none of the BrOOBr excited states are bound. However, at MRCI level of theory, all the BrOOBr states studied in this work are bound or slightly bound at the Frank-Condon region. In addition, the scalar relativistic effect and the spin-orbital coupling effect on the vertical excitation energies of the electronic states of BrOOBr were estimated.     

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.     

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.     

Tetrahedrane--dossier of an unknown

To probe whether tetrahedrane should be isolable the thermodynamics and kinetics of C4H4 singlet and triplet structures were studied extensively at the CCSD(T)/cc-pVTZ//CCSD(T)/cc-pVDZ, CCSD(T)/cc-pVDZ, CCSD(T)/cc-pVDZ//B3 LYP/6-311G**, and B3 LYP/6-311G** levels of theory. The reaction of cyclopropene with atomic carbon, which was previously suggested to involve tetrahedrane as a reactive intermediate, was re-examined experimentally with low-temperature matrix-isolation techniques. While experimental and theoretical results exclude the intermediacy of tetrahedrane in the above reaction, it is predicted to be an isolable molecule. Among the many C4H4 species, we pay special attention to the electronic effects on the ground state multiplicity of the respective carbenes.     

Theoretical study on the potential energy surface of NC<Subscript>3</Subscript>P isomers

Density functional theory (DFT) calculations have been used to study the isomerization process in the NC₃P system. At the DFT/B3LYP/6-311G(d) level, 28 triplet and 28 singlet minima were obtained on their respective potential energy surfaces. The linear triplet ³NCCCP is the lowest-energy structure among the isomers. On the triplet PES, only linear isomers ³NCCCP, ³CNCCP, ³CCCNP, and ³CCNCP possess great kinetic and thermodynamic stabilities to exist under low-temperature conditions (such as in the dense interstellar clouds). At the same time, one chain-like and four three-membered-ring isomers on the singlet PES have been located with high kinetic and thermodynamic stabilities. Further CCSD(T)/6-311G(2df)//QCISD/6-311G(d), CCSD(T)/cc-pVTZ//DFT/B3LYP/cc-pVTZ, and CASPT2(14,12)/cc-pVQZ//CASSCF(14,12)/cc-p VQZ calculations are performed on the structures, frequencies, and energies of the relevant species. The bonding natures were analyzed and the results were compared with the analogous NC₃N and NC₂P molecules so as to aid their future experimental or astrophysical detection.     

Gibbs energy of activation for thermal isomerization of (1Z)‐acetaldehyde hydrazone and (1Z)‐acetaldehyde N,N‐dimethylhydrazone by Gaussian‐4 theory and CCSD(T)/CBS computations

In this article, we examined the Gibbs energy of activation for the Z/E thermal isomerization reaction of (1Z)‐acetaldehyde hydrazone and (1Z)‐acetaldehyde N,N‐dimethylhydrazone, at 298.15 K in the solvent of cyclohexane. We carried out computations employing both the Gaussian‐4 (G4) theory and the coupled cluster method using both single and double substitutions and triple excitations noniteratively, CCSD(T). The CCSD(T) energy is extrapolated to the complete basis set (CBS). We compared the calculated results to the available experimental observation. It appeared that both G4 and CCSD(T)/CBS computations overestimated the experimental value by as much as about 6 and 12 kcal/mol in the present two cases. We discussed possible sources of error and proposed the experimental kinetic data could be questionable. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

An ab initio study of the Cl(P)+C2H6→C2H5+HCl abstraction reaction

Electronic energies, geometries, and harmonic vibration frequencies for the reactants, products, and transition state for the Cl(³P)+C₂H₆→C₂H₅+HCl abstraction reaction were evaluated at the HF and MP2 levels using several correlation consistent polarized-valence basis sets. Single-point calculations at PMP2, MP4, QCISD(T), and CCSD(T) levels were also carried out. The values of the forward activation energies obtained at the MP4/cc-pVTZ, QCISD(T)/cc-pVTZ, and CCSD(T)/cc-pVTZ levels using the MP2/cc-pVTZ structures are equal to −0.1, −0.4, and −0.3 kcal/mol, respectively. The experimental value is equal to 0.3±0.2 kcal/mol. We found that the MP2/aug-cc-pVTZ adiabatic vibration energy for the reaction (−2.4 kcal/mol) agrees well with the experimental value −(2.2–2.6) kcal/mol. Rate constants calculated with the zeroth-order interpolated variational transition state (IVTST-0) method are in good agreement with experiment. In general, the theoretical rate constants differ from experiment by, at most, a factor of 2.6.     

Ab initio calculation of bowl, cage, and ring isomers of C20 and C20-

High-level ab initio calculations have been carried out to reexamine relative stability of bowl, cage, and ring isomers of C(20) and C(20)(-). The total electronic energies of the three isomers show different energy orderings, strongly depending on the hybrid functionals selected. It is found that among three popular hybrid density-functional (DF) methods B3LYP, B3PW91, PBE1PBE, and a new hybrid-meta-DF method TPSSKCIS, only the PBE1PBE method (with cc-pVTZ basis set) gives qualitatively correct energy ordering as that predicted from ab initio CCSD(T)/cc-pVDZ [CCSD(T)-coupled-cluster method including singles, doubles, and noniterative perturbative triples; cc-pVDZ-correlation consistent polarized valence double zeta] as well as from MP4(SDQ)/cc-pVTZ [MP4-fourth-order Moller-Plesset; cc-pVTZ-correlation consistent polarized valence triple zeta] calculations. Both CCSD(T) and MP4 calculations indicate that the bowl is most likely the global minimum of neutral C(20) isomers, followed by the fullerene cage and ring. For the anionic counterparts, the PBE1PBE calculation also agrees with MP4/cc-pVTZ calculation, both predicting that the bowl is still the lowest-energy structure of C(20)(-) at T=0 K, followed by the ring and the cage. In contrast, both B3LYP/cc-pVTZ and B3PW91/cc-pVTZ calculations predict that the ring is the lowest-energy structure of C(20)(-). Apparently, this good reliability in predicting the energy ordering renders the hybrid PBE method a leading choice for predicting relative stability among large-sized carbon clusters and other carbon nanostructures (e.g., finite-size carbon nanotubes, nano-onions, or nanohorns). The relative stabilities derived from total energy with Gibbs free-energy corrections demonstrate a changing ordering in which ring becomes more favorable for both C(20) and C(20)(-) at high temperatures. Finally, photoelectron spectra (PES) for the anionic C(20)(-) isomers have been computed. With binding energies up to 7 eV, the simulated PES show ample spectral features to distinguish the three competitive C(20)(-) isomers.     

H-C-SiH3: direct generation and spectroscopic identification of ethylidene's cousin

Ground-state triplet silaethylidene, generated directly by the reaction of 3P carbon atoms with silane under matrix isolation conditions in solid Ar (10-12 K), has been thoroughly characterized by the EPR and IR spectra of both the parent and perdeuterated isotopologs. A theoretical anharmonic vibrational analysis based on a CCSD(T)/cc-pVTZ complete quartic force field gave remarkable agreement with the experimental IR fundamentals, generally within 10 cm-1 and without any empirical scaling of the ab initio frequencies. Silaethylidene exhibits a Cs minimum with a H-C-Si angle near 153 degrees , but the barrier to H-C-Si linearity (C3v symmetry) is only 0.24 kcal mol-1. This minuscule barrier can be surmounted by zero-point vibrations, as evident from the EPR data. The triplet stabilizing effect of the electropositive SiH3 group amounts to about 15 kcal mol-1.     

Ion-molecule rate constants and branching ratios for the reaction of N(3) (+)+O(2) from 120 to 1400 K

The kinetics of the reaction of N(3) (+) with O(2) has been studied from 120 to 1400 K using both a selected ion flow tube and high-temperature flowing afterglow. The rate constant decreases from 120 K to approximately 1200 K and then increases slightly up to the maximum temperature studied, 1400 K. The rate constant compares well to most of the previous measurements in the overlapping temperature range. Comparing the results to drift tube data shows that there is not a large difference between increasing the translational energy available for reaction and increasing the internal energy of the reactants over much of the range, i.e., all types of energies drive the reactivity equally. The reaction produces both NO(+) and NO(2) (+), the latter of which is shown to be the higher energy NOO(+) linear isomer. The ratio of NOO(+) to NO(+) decreases from a value of over 2 at 120 K to less than 0.01 at 1400 K because of dissociation of NOO(+) at the higher temperatures. This ratio decreases exponentially with increasing temperature. High-level theoretical calculations have also been performed to compliment the data. Calculations using multi-reference configuration interaction theory at the MRCISD(Q)/cc-pVTZ level of theory show that singlet NOO(+) is linear and is 4.5 eV higher in energy than ONO(+). A barrier of 0.9 eV prevents dissociation into NO(+) and O((1)D); however, a crossing to a triplet surface connects to NO(+) and O((3)P) products. A singlet and a triplet potential energy surface leading to products have been determined using coupled cluster theory at the CCSD(T)/aug-cc-pVQZ level on structures optimized at the Becke3-Lee, Yang, and Parr (B3LYP)/aug-cc-pVTZ level of theory. The experimental results and reaction mechanism are evaluated using these surfaces.     

A theoretical study of the H-abstraction reactions from HOI by moist air radiolytic products (H, OH, and O (3P)) and iodine atoms (2P(3/2))

The rate constants of the reactions of HOI molecules with H, OH, O ((3)P), and I ((2)P(3/2)) atoms have been estimated over the temperature range 300-2500 K using four different levels of theory. Geometry optimizations and vibrational frequency calculations are performed using MP2 methods combined with two basis sets (cc-pVTZ and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVTZ, cc-pVQZ, 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Reaction enthalpies at 0 K were calculated at the CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory and compared to the experimental values taken from the literature. Canonical transition-state theory with an Eckart tunneling correction is used to predict the rate constants as a function of temperature. The computational procedure has been used to predict rate constants for H-abstraction elementary reactions because there are actually no literature data to which the calculated rate constants can be directly compared. The final objective is to implement kinetics of gaseous reactions in the ASTEC (accident source term evaluation code) program to improve speciation of fission products, which can be transported along the reactor coolant system (RCS) of a pressurized water reactor (PWR) in the case of a severe accident.     

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.     

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.     

A coupled-cluster study of the structure and vibrational spectra of pyrazole and imidazole

Harmonic force fields were calculated at the corresponding optimized geometries for pyrazole and imidazole at the HF, B3LYP, MP2, CCSD and CCSD(T) levels using the 6-31G* basis set and at the HF and B3LYP levels using the cc-pVTZ basis set. The agreement between the calculated and experimental geometries by the CCSD and CCSD(T) methods was generally similar to that obtained with the B3LYP and MP2 methods. The force fields were scaled using one-scale-factor (1SF), 3SF and 7SF scaling schemes. The scale factors were varied with respect to the experimental frequencies. Using 7SF scaling, the root-mean-square (RMS) deviation of the calculated frequencies from the experimental frequencies by the HF, B3LYP, MP2, CCSD and CCSD(T) methods and the 6-31G* basis set was 16, 7, 13, 11 and 11 cm(-1), respectively. This shows that the B3LYP method is preferred for force field calculations over the perturbative MP2, CCSD and CCSD(T) methods. Using 1SF scaling, the CCSD(T) scale factor was 0.931, the highest among the five methods used but close to that obtained with the B3LYP method and the cc-pVTZ basis set with lower RMS deviation.     

S3O2--an unusual structure with pi*-pi* interaction

The structures and relative stabilities of 15 S3O2 isomers have been investigated by G3X(MP2), CCSD(T)/aug-cc-pVTZ and MRCI/CASSCF calculations. The global energy minimum is a three-membered sulfur ring with two adjacent sulfoxide groups in a trans conformation, i.e. a vic-disulfoxide of C2 symmetry. The SS bond lengths are 2.136 (2x) and 2.354 angstroms at the CCSD(T)/cc-pVTZ level of theory. There is a strong interaction between the pi* orbitals of the two S=O moieties both in the trans and in the almost degenerate cis conformer. The corresponding chain-like singlet and triplet isomers of connectivity OSSSO lie close in energy (ca. 67 kJ mol(-1)) while five-membered and branched four-membered rings are significantly less stable. The structure of S3O2 is in contrast to that of the isoelectronic analogue S5, which exists as a five-membered twisted heterocycle.     

Ab initio study of hydrazinyl radical: toward a DMBE potential energy surface

A series of stationary structures of the hydrazinyl radical have been characterized by optimization at the CCSD(T)/cc-pVTZ level of theory. CCSD(T)/aug-cc-pVXZ single-point calculations have also been carried out at the optimized geometries with basis sets of different cardinal numbers (X = T, Q), which were used to obtain accurate energies via extrapolation to the complete basis set limit. A discussion on the analytical modeling of the potential energy surface of hydrazinyl is also presented.     

Theoretical study on the mechanism of the gas-phase reaction of diborane(3) anion with carbon disulfide

The complex potential energy surface of the gas-phase reaction of HB(H)BH- with CS2 to give three low-lying products [B2H3S]- + CS, [BH2CS]- + HBS, and [BH3CS] + BS-, involving nine [B2H3CS2]- isomers and 12 transition states, has been investigated at the CCSD(T)/6-311++G(d,p)/B3LYP/6-311++G(d,p) level. Our calculations are in harmony with the recent experimental and theoretical results, and reveal some new bonding and kinetic features of this reaction system. Our theoretical results may help the further identification of the products [BH2CS]- + HBS and [BH3CS] + BS- and may provide useful information on the chemical behaviors of other electron-deficient boron hydride anions.     

Prototypical triplet alkyl phosphonatocarbenes

The current case study focuses on the generation, identification, and characterization of two representative mono- and disubstituted alkyl phosphonatocarbenes by means of matrix isolation techniques in conjunction with density functional theory [B3LYP/6-311++G(d,p)] and coupled cluster [CCSD(T)/cc-pVXZ, X = D, T] computations. The EPR measurements identify both carbenes as triplet ground-state species with D values of 0.660 and 0.623 cm(-1), respectively, exhibiting persistency toward intramolecular reactions (the EPR signal observable in perfluoromethylcyclohexane up to around 70 K for the disubstituted molecule). While the reaction of the carbene center of the conformationally rich tetramethyl bisphosphonatocarbene with the CH bonds of the methyl groups leads to phosphaoxetane at room temperature, its fragmentation via a Wittig-type reaction during high vacuum flash pyrolysis (HVFP) results in dimethyl vinylphosphonate and methyl metaphosphate. The latter has been observed for the first time as an isolated entity.