We have investigated the bifurcation dynamics of the H2CO− + CH3Cl →CH3CH2O + Cl− (C-substitution channel)/H2CO + CH3 + Cl− (electron transfer channel) reaction using three different molecular dynamics (MD) methods: quasi-classical trajectory (QCT), classical MD, and ring-polymer MD simulations. All MD calculations were performed using the on-the-fly approach at the ab initio molecular orbital theory level. We found that the C-substitution channel is dominant in the classical and ring-polymer MD calculations, while the electron transfer channel is dominant in the QCT calculations. The different bifurcation mechanisms in the types of MD simulations were investigated using supervised machine-learning classification. This difference is attributed to the varying efficiency of the intramolecular vibrational redistribution process among the different MD methods. Thus, in the bifurcation mechanism of the H2CO− + CH3Cl reaction, the bifurcation fractions are primarily determined by the vibrational dynamics occurring in the later stages of the reaction.
Read full abstract