Abstract:

Vibrational energy flow and conformational transitions following excitation of the OH stretching mode of the most stable conformer of glycine are studied by classical trajectories. ``On the fly'' simulations with the PM3 semiempirical electronic structure method for the potential surface are used. Initial conditions are selected to correspond to the v = 1 excitation of the OH stretch. The main findings are: (1) An an equilibrium-like ratio is established between the populations of the 3 lowest-lying conformers after about 10 picoseconds. (2) There is a high probability throughout the 150 ps of the simulations for finding the molecule in geometries far from the equilibrium structures of the lowest-energy conformers. (3) Energy from the initial excited OH (v = 1) stretch flows preferentially to 5 other vibrational modes, including the bending motion of the H atom. (4) RRK theory yields conformational transition rates that deviate substantially from the classical trajectory results. Possible implication of these results for vibrational energy flow and conformational transitions in small biological molecules are discussed.