The C4H7+ cation. A theoretical investigation

The potential energy surface of the C4H7+ cation has been investigated with ab initio quantum chemical theory. Extended basis set calculations, including electronic correlation, show that cyclobutyl and cyclopropylcarbinyl cation are equally stable isomers. The saddle point connecting these isomers lies 0.6 kcal/mol above the minima. The global C4H7+ minimum corresponds to the 1-methylallyl cation, which is 9.0 kcal/mol more stable than the cyclobutyl and the cyclopropylcarbinyl cation and 9.5 kcal/mol below the 2-methylallyl cation. These results are in excellent agreement with experimental data.     

Isotope effects on the unimolecular dissociation of ionized 3-methyl-2-butanol: reactions via a long-lived C-H-C hydrogen-bridged ion-neutral complex

A pronounced isotope effect causes metastable CD3CHOHCH(CH3)2+* ions to expell C3H6D2 in preference to C3H7D in a ratio of approximately 33:1; a number of related compounds show similar effects. High-level ab initio calculations suggest that the reactant alcohol molecular ion possesses an extraordinarily long alpha-carbon-carbon bond and that the reaction proceeds via the formation of an intermediate hydrogen-bridged complex of propane and ionized vinyl alcohol, in which the bridging hydrogen atom is almost midway between the two carbon termini. The isotopic preference reflects the difference between the zero-point vibrational energies of the isotopically different product pairs rather than kinetic isotope effects on the hydrogen atom transfer reactions that precede dissociation.     

Photodissociation dynamics of pyrimidine

Photodissociation of pyrimidine at 193 and 248 nm was investigated separately using vacuum ultraviolet photoionization at 118.4 and 88.6 nm and multimass ion imaging techniques. Six dissociation channels were observed at 193 nm, including C4N2H4 --> C4N2H3 + H and five ring opening dissociation channels, C4N2H4 --> C3NH3 + HCN, C4N2H4 --> 2C2NH2, C4N2H4 --> CH3N + C3NH, C4N2H4 --> C4NH2 + NH2, and C4N2H4 --> CH2N + C3NH2. Only the first four channels were observed at 248 nm. Photofragment translational energy distributions and dissociation rates indicate that dissociation occurs in the ground electronic state after internal conversion at both wavelengths. The dissociation rates were found to be >5 x 10(7) and 1 x 10(6) s(-1) at 193 and 248 nm, respectively. Comparison with the potential energies from ab initio calculations have been made.     

Sensitivity of the H + CH3 → CH4 recombination rate constant to the shape of the CH stretching potential

Using a potential-energy surface obtained in part from ab initio calculations, the H + CH₃ → CH₄ bimolecular rate constant at T = 300 K is determined from a Monte Carlo classical trajectory study. Representing the CH stretching potential with a standard Morse function instead ofthe ab initio curve increases the calculated rate constant by an order of magnitude. The experimental recombination rate constant is intermediate of the rate constants calculated with the Morse and ab initio stretching potentials.Two properties of the H + CH₃ α CH₄ potential-energy surface which significantly affect the recombination rate constant are the shape of the CH stretching potential and the attenuation of the H₃CH bending frequencies. Ab initio calculations with a hierarchy of basis sets and treatment of electron correlation indicate the latter is properly described [13]. The exact shape of the CH stretching potential is not delineated by the ab initio calculations, since the ab initio calculations are not converged for bond lengths of 2.0–3.0 Å [12]. However, the form of this stretching potential deduced from the highest-level ab initio calculations, and fit analytically by eq. (2), is significantly different from a Morse function. The experimental recombination rate constant is intermediate of the rate constants calculated with the Morse and ab initio CH stretching potentials. This indicates that the actual CH potential energy curve lies between the Morse and ab initio curves. This is consistent with the finding that potential energy curves for diatomics are not well described by a Morse function [12].     

Photodissociation dynamics of C6HxF6-x (x = 1-4) at 193 nm

Photodissociation of fluorine-substituted benzenes, including 1,3-difluorobenzene, 1,2,4-trifluorobenzene, 1,2,4,5-tetrafluorobenzene, and pentafluorobenzene, at 193 nm under collision-free conditions has been studied in separate experiments using multimass ion imaging techniques. HF elimination was found to be the major dissociation channel for all of these molecules. Small amounts of photofragments of C(6)H(3)F(2) and C(6)H(2)F(3) from 1,3-difluorobenzene and 1,2,4-trifluorobenzene, respectively, were also observed. They correspond to the minor dissociation channel of hydrogen elimination. Dissociation rates and fragment translational energy distributions obtained from experimental measurements suggest that HF and hydrogen elimination reactions occur in the ground electronic state. The potential energy surface obtained from ab initio calculations indicates that the four-center reaction in the ground electronic state is the major dissociation mechanism for the HF eliminations. A comparison with the RRKM calculation has been made.     

Photodissociation dynamics of pyridine

Photodissociation of pyridine, 2,6-d2-pyridine, and d5-pyridine at 193 and 248 nm was investigated separately using multimass ion imaging techniques. Six dissociation channels were observed at 193 nm, including C5NH5 --> C5NH4 + H (10%) and five ring opening dissociation channels, C5NH5 --> C4H4 + HCN, C5NH5 --> C3H3 + C2NH2, C5NH5 --> C2H4 +C3NH, C5NH5 --> C4NH2 + CH3 (14%), and C5NH5 --> C2H2 + C3NH3. Extensive H and D atom exchanges of 2,6-d2-pyridine prior to dissociation were observed. Photofragment translational energy distributions and dissociation rates indicate that dissociation occurs in the ground electronic state after internal conversion. The dissociation rate of pyridine excited by 248-nm photons was too slow to be measured, and the upper limit of the dissociation rate was estimated to be 2x10(3) s(-1). Comparisons with potential energies obtained from ab initio calculations and dissociation rates obtained from the Rice-Ramsperger-Kassel-Marcus theory have been made.     

Structure, stability and vibrational spectrum of the hydrogen-bonded complex between HNO3 and H2O. Ab initio and DFT studies

The structure, stability and vibrational spectrum of the binary complex between HONO2 and H2O have been investigated using ab initio calculations at SCF and MP2 levels with different basis sets and B3LYP/6-31G(d,p) calculations. Full geometry optimization was made for the complex studied. It was established that the hydrogen-bonded H2O...HONO2 complex has a planar structure. The corrected values of the dissociation energy at the SCF and MP2 levels and B3LYP calculations are indicative of relatively strong OH...O hydrogen-bonded interaction. The changes in the vibrational characteristics (vibrational frequencies and infrared intensities) arising from the hydrogen bonding between HONO2 and H2O have been estimated by using the ab initio calculations at SCF and MP2 levels and B3LYP/6-31G(d,p) calculations. It was established that the most sensitive to the complexation is the stretching O-H vibration from HONO2. In agreement with the experiment, its vibrational frequency in the complex is shifted to lower wavenumbers. The predicted frequency shift with the B3LYP/6-31G(d,p) calculations (-439 cm(-1)) is in the best agreement with the experimentally measured (-498 cm(-1)). The intensity of this vibration increases dramatically upon hydrogen bonding. The ab initio calculations at the SCF level predict an increase up to five times; at the MP2 level up to 10 times and the B3LYP/6-31G(d,p) predicted increase is up to 17 times. The good agreement between the predicted values of the frequency shifts and those experimentally observed show that the structure of the hydrogen-bonded complex H2O...HONO2 is reliable.     

Hydrogen spillover mechanism on a Pd-doped Mg surface as revealed by ab initio density functional calculation

The hydrogenation kinetics of Mg is slow, impeding its application for mobile hydrogen storage. We demonstrate by ab initio density functional theory (DFT) calculations that the reaction path can be greatly modified by adding transition metal catalysts. Contrasting with Ti doping, a Pd dopant will result in a very small activation barrier for both dissociation of molecular hydrogen and diffusion of atomic H on the Mg surface. This new computational finding supports-for the first time by ab initio simulation-the proposed hydrogen spillover mechanism for rationalizing experimentally observed fast hydrogenation kinetics for Pd-capped Mg materials.     

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.     

A vacuum ultraviolet photoionization mass spectrometric study of acetone

The photoionization and dissociative photoionization of acetone have been studied at the photon energy range of 8-20 eV. Photoionization efficiency spectra for ions CH3COCH3+, CH3+, C2H3+, C3H3+, C3H5+, CH(2-)CO+, CH3CO+, C3H4O+, and CH3COCH2+ have been measured. In addition, the energetics of the dissociative photoionization has been examined by ab initio Gaussian-3 (G3) calculations. The computational results are useful in establishing the dissociation channels near the ionization thresholds. With the help of G3 results, the dissociation channels for the formation of the fragment ions CH3CO+, CH2CO+, CH3+, C3H3+, and CH3COCH2+ have been established. The G3 results are in fair to excellent agreement with the experimental data.     

C?H Bond dissociation of acetylene: Local density functional calculations

The C? H bond dissociation energy of acetylene was computed by both ab initio approaches and density functional theory in a local density approximation (DFT–LDA ). Structures and energies for acetylene and its dissociation products (the ethynyl and hydrogen radicals) are presented and compared. Using directly computed HCCH and HCC· energies and the exact H· value, the DFT–LDA calculations are found to yield C? H dissociation energies ranging from 129 to 131 kcal/mol, in good agreement with recent experimental and the highest level theoretical results. The DFT–LDA results show little dependence upon the computational procedure used to obtain geometries.     

Theoretical prediction of the ionization energies of the C4H7 radicals: 1-methylallyl, 2-methylallyl, cyclopropylmethyl, and cyclobutyl radicals

The ionization energies (IEs) for the 1-methylallyl, 2-methylallyl, cyclopropylmethyl, and cyclobutyl radicals have been calculated by the wave function based ab initio CCSD(T)/CBS approach, which involves the approximation to the complete basis set (CBS) limit at the coupled cluster level with single and double excitations plus quasiperturbative triple excitation [CCSD(T)]. The zero-point vibrational energy correction, the core-valence electronic correction, and the scalar relativistic effect correction are included in these calculations. The present CCSD(T)/CBS results are then compared with the IEs determined in the photoelectron experiment by Schultz et al. [J. Am. Chem. Soc. 106, 7336 (1984)] The predicted IE value (7.881 eV) of 2-methylallyl radical is found to compare very favorably with the experimental value of 7.90+/-0.02 eV. Two ionization transitions for cis-1-methylallyl and trans-1-methylallyl radicals have been considered here. The comparison between the predicted IE values and the previous measurements shows that the photoelectron peak observed by Schultz et al. likely corresponds to the adiabatic ionization transition for the trans-1-methylallyl radical to form trans-1-methylallyl cation. Although a precise IE value for the cyclopropylmethyl radical has not been directly determined, the experimental value deduced indirectly using other known energetic data is found to be in good accord with the present CCSD(T)/CBS prediction. We expect that the Franck-Condon factor for ionization transition of c-C4H7-->bicyclobutonium is much less favorable than that for ionization transition of c-C4H7-->planar-C4H7+, and the observed IE in the previous photoelectron experiment is likely due to the ionization transition for c-C4H7-->planar-C4H7+. Based on our CCSD(T)/CBS prediction, the ionization transition of c-C4H7-->bicyclobutonium with an IE value around 6.92 eV should be taken as the adiabatic ionization transition for the cyclobutyl radical. The present study provides support for the conclusion that the CCSD(T)/CBS approach with high-level energetic corrections can be used to provide reliable IE predictions for C4 hydrocarbon radicals with an uncertainty of +/-22 meV. The CCSD(T)/CBS predictions to the heats of formation for the aforementioned radicals and cations are also presented.     

Ab initio multireference configuration-interaction study of hydrogen molecule activation by Cs-promoted Pt clusters

The adsorption of the H2 molecule on CsnPt(5-n) bcc (111) clusters for Cs/Pt rates of 20%, 40%, and 80% is studied using ab initio multiconfigurational self-consistent field plus multireference configuration-interaction variational and perturbative calculations. The H2 interaction with the clusters is studied in ground and excited states with geometry optimization, where the hydrogen adsorption takes place by a Pt atom. These calculations are compared with those of H2 adsorption on Pt4. The most stable configurations of Cs/Pt4 and Cs2Pt3 clusters (Cs/Pt rates of 20% and 40%) are a doublet and a closed-shell singlet, respectively. Both clusters capture and activate the hydrogen molecule and their behaviors resemble Pt4. The H2 capture distances are, respectively, similar and smaller than Pt4 capture distances, while the H-H bond dissociation distances are similar and bigger than those of Pt4; however, none of them presents activation barriers. The most stable Cs4Pt cluster (Cs/Pt rate of 80%) is also a closed-shell singlet; it also captures and activates the hydrogen molecule and shows a different behavior as compared with Cs/Pt4, Cs2Pt3, and Pt4 clusters. The capture distance is quite smaller and is obtained after surmounting an activation barrier. For all clusters studied here, no hydrogen absorption was observed, only the adsorption of H2.     

Photodissociation dynamics of benzyl alcohol at 193 nm

Photodissociation dynamics of benzyl alcohol, C(6)H(5)CH(2)OH and C(6)H(5)CD(2)OH, in a molecular beam was investigated at 193 nm using multimass ion imaging techniques. Four dissociation channels were observed, including OH elimination and H(2)O elimination from the ground electronic state, H atom elimination (from OH functional group), and CH(2)OH elimination from the triplet state. The dissociation rate on the ground state was found to be 7.7 × 10(6) s(-1). Comparison to the potential energy surfaces from ab initio calculations, dissociation rate, and branching ratio from Rice-Ramsperger-Kassel-Marcus calculations were made.     

Theoretical Studies on Structures and Properties of Endohedral Fullerenes Complexes： XHn@C32（X=F,O, N,C; n=l-4）

Theoretical studies on structures and properties of endohedral fullerene complexes formed by encapsulating small molecules of HF, H20, NH3, and CH4 in a C32 fullerene cage, were carried out by ab initio method. Current calculations reveal that these processes to encase them in fullerene are energetically unfavorable because of the small cavity size of C32. The red shift in the F-H stretching frequency indicates the potential existence of hydrogen bonding between the HF molecule and the carbon cage.     

Endgroup-assisted siloxane bond cleavage in the gas phase

Unimolecular dissociation of H(2)N(CH(2))(3)SiOSi(CH(2))(3)NH(3)(+) generates SiC(5)H(16)NO(+) and SiC(5)H(14)N(+). The formation of SiC(5)H(16)NO(+) involves dissociation of a Si[bond]O bond and formation of an O[bond]H bond through rearrangement. The fragmentation mechanism was investigated utilizing ab initio calculations and Fourier transform ion cyclotron resonance (FTICR) mass spectrometry in combination with hydrogen/deuterium (H/D) exchange reactions. Sustained off-resonance irradiation collision-induced dissociation (SORI-CID) studies of the fully deuterated ion D(2)N(CH(2))(3)SiOSi(CH(2))(3)ND(3)(+) provided convincing evidence for a backbiting mechanism which involves hydrogen transfer from the terminal amine group to the oxygen to form a silanol-containing species. Theoretical calculations indicated decomposition of H(2)N(CH(2))(3)SiOSi(CH(2))(3)NH(3)(+) through a backbiting mechanism is the lowest energy decomposition channel, compared with other alternative routes. Two mechanisms were proposed for the fragmentation process which leads to the siloxane bond cleavage and the SORI-CID results of partially deuterated precursor ions suggest both mechanisms should be operative. Rearrangement to yield a silanol-containing product ion requires end groups possessing a labile hydrogen atom. Decomposition of disiloxane ions with end groups lacking labile hydrogen atoms yielded product ions from direct bond cleavages.     

Hydrogen-bonded complexes of nitrous and nitric acids with ethene Ab initio and DFT studies

The complexes formed by ethene with nitric and nitrous (trans and cis) acids have been investigated by ab initio (SCF and MP2) and B3LYP calculations with 6-311++G(d,p) basis set. Full geometry optimisation has been performed for the complexes studied. The most stable structures of the complexes are established. Bearing in mind the corrected values of the dissociation energy the studied hydrogen-bonded complexes can be ordered as follows: C(2)H(4)...HONO(2)>C(2)H(4)...HONO-trans>C(2)H(4)...HONO-cis. In the complexes the acids act as proton donors forming the pi-type of hydrogen bond with ethene. The predicted changes in the vibrational characteristics (vibrational frequencies and infrared intensities) arising from the hydrogen bonding are in good agreement with the experimentally measured. The predicted frequency shift of the stretching OH vibration in the nitric acid is largest (-210 cm(-1)), followed by the shifts in the trans-HONO (-141 cm(-1)) and cis-HONO (-109 cm(-1)). The calculations predict an increase of the IR intensity of the stretching O-H vibration in the complexes from 6 to 10 times.     

A gate to organokrypton chemistry: HKrCCH

An organic molecule containing krypton, HKrCCH, is reported. The preparation of HKrCCH includes 193-nm photolysis of H2C2/Kr solid mixtures at 8 K and subsequent thermal mobilization of hydrogen atoms at >/=30 K. The identification is based on infrared absorption spectroscopy and supported by ab initio calculations which show ionic and covalent contributions to the bonding. We believe that a series of similar organokrypton molecules can be prepared as computationally demonstrated for HKrC4H and HKrC3H3. These results feature a generally novel way for activating chemically the H-CC- group, which can find practical applications of the krypton catalysis.     

Calculating geochemical reaction pathways--exploration of the inner-sphere water exchange mechanism in Al(H2O)6(3+)(aq) + nH2O with ab Initio calculations and molecular dynamics

We have simulated exchange of inner-sphere and bulk water molecules for different sizes of Al3+(aq) clusters, Al(H2O)63+ + nH2O for n = 0, 1, 6, or 12, with ab initio and molecular dynamics simulations, in order to understand how robust the ab initio method is for identifying hydrolytic reaction pathways of particular importance to geochemistry. In contrast to many interfacial reactions, this particular elementary reaction is particularly simple and well-constrained by experiment. Nevertheless, we find that a rich array of parallel reaction pathways depend sensitively on the details of the solvation sphere and structure and that larger clusters are not necessarily better. Inner-sphere water exchange in Al3+(aq) may occur through two Langford-Gray dissociative pathways, one in which the incoming and outgoing waters are cis, the other in which they are trans to one another. A large majority of exchanges in the molecular dynamics simulations occurred via the trans mechanism, in contrast to the predictions of the ab initio method. In Al(H2O)63+ + H2O, the cis mechanism has a transition state of 84.3 kJ/mol, which is in good agreement with previous experimental and ab initio results, while the trans mechanism has only a saddle point with two negative frequencies, not a transition state, at 89.7 kJ/mol. In addition to the exchange mechanisms, dissociation pathways could be identified that were considerably lower in energy than experiment and varied considerably between 60 and 100 kJ/mol, depending on the particular geometry and cluster size, with no clear relation between the two. Ab initio calculations using large clusters with full second coordination spheres (n = 12) were unable to find dissociation or exchange transition states because the network of hydrogen bonds in the second coordination sphere was too rigid to accommodate the outgoing inner-sphere water. Our results indicate that caution should surround ab initio simulation of complicated dynamic processes such as hydrolysis, ion exchange, and interfacial reactions that involve several steps. Dynamic methods of simulation need to accompany static methods such as ab initio calculation, and it is best to consider simulated pathways as hypotheses to be tested experimentally rather than definitive properties of the reaction.     

Communication: theoretical exploration of Au(+)+H2, D2, and HD reactive collisions

A quasi-classical study of the endoergic Au(+)((1)S) + H(2)(X(1)Σ(g)(+)) → AuH(+) ((2)Σ(+)) + H((2)S) reaction, and isotopic variants, is performed to compare with recent experimental results [F. Li, C. S. Hinton, M. Citir, F. Liu, and P. B. Armentrout, J. Chem. Phys. 134, 024310 (2011)]. For this purpose, a new global potential energy surface has been developed based on multi-reference configuration interaction ab initio calculations. The quasi-classical trajectory results show a very good agreement with the experiments, showing the same trends for the different isotopic variants of the hydrogen molecule. It is also found that the total dissociation into three fragments, Au(+)+H+H, is the dominant reaction channel for energies above the H(2) dissociation energy. This results from a well in the entrance channel of the potential energy surface, which enhances the probability of H-Au-H insertion.