An automated method to find reaction mechanisms and solve the kinetics in organometallic catalysis.

A novel computational method is proposed in this work for use in discovering reaction mechanisms and solving the kinetics of transition metal-catalyzed reactions. The method does not rely on either chemical intuition or assumed a priori mechanisms, and it works in a fully automated fashion. Its core is a procedure, recently developed by one of the authors, that combines accelerated direct dynamics with an efficient geometry-based post-processing algorithm to find transition states (Martinez-Nunez, E., J. Comput. Chem.2015, 36, 222-234). In the present work, several auxiliary tools have been added to deal with the specific features of transition metal catalytic reactions. As a test case, we chose the cobalt-catalyzed hydroformylation of ethylene because of its well-established mechanism, and the fact that it has already been used in previous automated computational studies. Besides the generally accepted mechanism of Heck and Breslow, several side reactions, such as hydrogenation of the alkene, emerged from our calculations. Additionally, the calculated rate law for the hydroformylation reaction agrees reasonably well with those obtained in previous experimental and theoretical studies.

Trajectory dynamics study of collision-induced dissociation of the Ar + CH4 reaction at hyperthermal conditions: vibrational excitation and isotope substitution.

We investigate the role of vibrational energy excitation of methane and two deuterated species (CD(4) and CH(2)D(2)) in the collision-induced dissociation (CID) process with argon at hyperthermal energies. The quasi-classical trajectory method has been applied, and the reactive Ar + CH(4) system has been modeled by using a modified version of the CH(4) potential energy surface of Duchovic et al. (J. Phys. Chem. 1984, 88, 1339) and the Ar-CH(4) intermolecular potential function obtained by Troya (J. Phys. Chem. A 2005, 109, 5814). This study clearly shows that CID is markedly enhanced with vibrational excitation and, to a lesser degree, with collision energy. In general, CID increases by exciting stretch vibrational modes of the reactant molecule. For the direct dissociation of CH(4), however, the CID cross sections appear to be essentially independent of which vibrational mode is initially excited. In all situations studied, the CID cross sections are always greater for the Ar + CD(4) reaction than for the Ar + CH(4) one, the Ar + CH(2)D(2) being an intermediate situation. A detailed analysis of the energy transfer processes, including their relation with CID, is also presented.

On the conformational memory in the photodissociation of formic acid.

The photodissociation of formic acid at 248 and 193 nm was investigated by classical trajectory and RRKM calculations using an interpolated potential energy surface, iteratively constructed using the B3LYP/aug-cc-pVDZ level of calculation. Several sampling schemes in the ground electronic state were employed to explore the possibility of conformational memory in formic acid. The CO/CO2 branching ratios obtained from trajectories initiated at the cis and at the trans conformers are almost identical to each other and in very good accordance with the RRKM results. In addition, when a specific initial excitation that simulates more rigorously the internal conversion process is used, the calculated branching ratio does not vary with respect to those obtained from cis and trans initializations. This result is at odds with the idea of conformational memory in the ground state proposed recently for the interpretation of the experimental results. It was also found that the calculated CO vibrational distributions after dissociation of the parent molecule at 248 nm are in agreement with the experimental available data.

Trajectory dynamics study of the Ar + CH4 dissociation reaction at high temperatures: the importance of zero-point-energy effects.

Large-scale classical trajectory calculations have been performed to study the reaction Ar + CH4--> CH3 +H + Ar in the temperature range 2500 < or = T/K < or = 4500. The potential energy surface used for ArCH4 is the sum of the nonbonding pairwise potentials of Hase and collaborators (J. Chem. Phys. 2001, 114, 535) that models the intermolecular interaction and the CH4 intramolecular potential of Duchovic et al. (J. Phys. Chem. 1984, 88, 1339), which has been modified to account for the H-H repulsion at small bending angles. The thermal rate coefficient has been calculated, and the zero-point energy (ZPE) of the CH3 product molecule has been taken into account in the analysis of the results; also, two approaches have been applied for discarding predissociative trajectories. In both cases, good agreement is observed between the experimental and trajectory results after imposing the ZPE of CH3. The energy-transfer parameters have also been obtained from trajectory calculations and compared with available values estimated from experiment using the master equation formalism; in general, the agreement is good.