Tunneling of coupled methyl quantum rotors in 4-methylpyridine: Single rotor potential versus coupling interaction.

Abstract

We study the influence of rotational coupling between a pair of methyl rotators on the tunneling spectrum in condensed phase. Two interacting adjacent methyl groups are simulated within a coupled-pair model composed of static rotational potential created by the chemical environment and the interaction potential between two methyl groups. We solve the two-dimensional time-independent Schrödinger equation analytically by expanding the wave functions on the basis set of two independent free-rotor functions. We investigate three scenarios which differ with respect to the relative strength of single-rotor and coupling potential. For each scenario, we illustrate the dependence of the energy level scheme on the coupling strength. It is found that the main determinant of splitting energy levels tends to be a function of the ratio of strengths of coupling and single-rotor potential. The tunnel splitting caused by coupling is maximized for the coupled rotors in which their total hindering potential is relatively shallow. Such a weakly hindered methyl rotational potential is predicted for 4-methylpyridine at low temperature. The experimental observation of multiple tunneling peaks arising from a single type of methyl group in 4-methylpyridine in the inelastic neutron scattering spectrum is widely attributed to the rotor-rotor coupling. In this regard, using a set of first-principles calculations combined with the nudged elastic band method, we investigate the rotational potential energy surface (PES) of the coaxial pairs of rotors in 4-methylpyridine. A Numerov-type method is used to numerically solve the two-dimensional time-independent Schrödinger equation for the calculated 2D-density functional theory profile. Our computed energy levels reproduce the observed tunneling transitions well. Moreover, the calculated density distribution of the three methyl protons resembles the experimental nuclear densities obtained from the Fourier difference method. By mapping the calculated first-principles PES on the model, it is confirmed that the hindering potential in 4-methylpyridine consists of proportionally shallow single-rotor potential to coupling interaction.