Search

Mode-dependent H atom tunneling dynamics of the O–H bond predissociation of the S1 phenol has been theoretically analyzed. As the tunneling is governed by the complicated multi-dimensional potential energy surfaces that are dynamically shaped by the upper-lying S1(ππ*)/S2(πσ*) conical intersection, the mode-specific tunneling dynamics of phenol (S1) has been quite formidable to be understood. Herein, we have examined the topography of the potential energy surface along the particular S1 vibrational mode of interest at the nuclear configurations of the S1 minimum and S1/S2 conical intersection. The effective adiabatic tunneling barrier experienced by the reactive flux at the particular S1 vibrational mode excitation is then uniquely determined by the topographic shape of the potential energy surface extended along the conical intersection seam coordinate associated with the particular vibrational mode. The resultant multi-dimensional coupling of the specific vibrational mode to the tunneling coordinate is then reflected in the mode-dependent tunneling rate as well as nonadiabatic transition probability. Remarkably, the mode-specific experimental result of the S1 phenol tunneling reaction [K. C. Woo and S. K. Kim, J. Phys. Chem. A 123, 1529–1537 (2019)] (in terms of the qualitative and relative mode-dependent dynamic behavior) could be well rationalized by semi-classical calculations based on the mode-specific topography of the effective tunneling barrier, providing the clear conceptual insight that the skewed potential energy surfaces along the conical intersection seam (strongly or weakly coupled to the tunneling reaction coordinate) may dictate the tunneling dynamics in the proximity of the conical intersection. [ABSTRACT FROM AUTHOR]

The πσ*‐mediated photochemistry of heteroaromatic molecules has provoked the investigation of the conical intersection dynamics. The Born–Oppenheimer approximation fails at the conical intersection where the S1 (ππ*) and S2 (πσ*) states cross. The nonadiabatic transitions are much influenced by the nuclear configuration of the reactive flux particularly in the curve‐crossing region encountered along the reaction pathway. In this article, we focus on the tunneling dynamics of phenols and thiophenols. The O (S)H bond cleavage occurs via tunneling through the barrier which is dynamically shaped by the upper‐lying S1/S2 conical intersection in terms of the couplings at the individual branching planes as well as along the (3N‐8) dimensional seam coordinates. State‐specific tunneling rates and their interpretation are given for phenol, substituted phenols, thiophenol, ortho‐substituted thiophenols, and benzenediols including their 1:1 water clusters. The completely orthogonal modes to the tunneling coordinate are very critical in the dynamic shaping of the reaction barrier. [ABSTRACT FROM AUTHOR]