Post-Transition-State Direct Dynamics Simulations on the Ozonolysis of Catechol

On-the-fly dynamics simulations are performed for the reaction of catechol + O3. The post transition state (TS) dynamics is studied at temperatures of 400 and 500 K. The PM7 semiempirical method is employed for calculating the potential energy gradient needed for integrating Hamilton’s equations of motion. This semiempirical method provides excellent agreement in terms of energy and geometry of the TSs as well as minimum energy states of the system with respect to B3LYP/6-311+G (2df, 2p) calculated results. In the dynamics, first, a peroxyacid is formed, which further dissociates to different fragments. Four major channels forming CO, CO2, H2O, and small carboxylic acid (SCA) fragments are seen in this reaction. Rates of each of the channels and the overall unimolecular reaction are calculated at both temperatures. Branching ratios of all these product channels are calculated and compared with experiment. The minimum energy profile of CO2, CO, and H2O channels are calculated. A qualitative estimate of activation energies for all the channels are obtained and compared with the explicit TS energies of three product channels, which ultimately correlate with the reaction probabilities.