National Research Council of Canada, Ottawa, Ontario K1A 0R6, Canada.
J Chem Phys. 2012 Dec 14;137(22):224105. doi: 10.1063/1.4769198.
A new instanton approach is reported to tunneling at zero-temperature in multidimensional (MD) systems in which a "light particle" is transferred between two equivalent "heavy" sites. The method is based on two concepts. The first is that an adequate MD potential energy surface can be generated from input of the stationary configurations only, by choosing as a basis the normal modes of the transition state. It takes the form of a double-minimum potential along the mode with imaginary frequency and coupling terms to the remaining (harmonic) oscillators. Standard integrating out of the oscillators gives rise to an effective 1D instanton problem for the adiabatic potential, but requires evaluation of a nonlocal term in the Euclidean action, governed by exponential (memory) kernels. The second concept is that this nonlocal action can be treated as a "perturbation," for which a new approximate instanton solution is derived, termed the "rainbow" solution. Key to the approach is avoidance of approximations to the exponential kernels, which is made possible by a remarkable conversion property of the rainbow solution. This leads to a new approximation scheme for direct evaluation of the Euclidean action, which avoids the time-consuming search of the exact instanton trajectory. This "rainbow approximation" can handle coupling to modes that cover a wide range of frequencies and bridge the gap between the adiabatic and sudden approximations. It suffers far fewer restrictions than these conventional approximations and is proving particularly effective for systems with strong coupling, such as proton transfer in hydrogen bonds. Comparison with the known exact instanton action in two-dimensional models and application to zero-level tunneling splittings in two isotopomers of malonaldehyde are presented to show the accuracy and efficiency of the approach.
报道了一种新的瞬子方法,用于在多维(MD)系统中零温度下的隧穿,其中“轻粒子”在两个等效的“重”位置之间转移。该方法基于两个概念。第一个是,可以仅通过输入静态配置来生成适当的 MD 势能表面,选择过渡态的本征模作为基础。它具有沿虚频模态的双势阱形式和与剩余(谐波)振荡器的耦合项。标准的振荡器积分会导致绝热势的有效 1D 瞬子问题,但需要在欧几里得作用中评估非局部项,由指数(记忆)核控制。第二个概念是,这个非局部作用可以被视为“微扰”,对于这个微扰,我们推导出了一个新的近似瞬子解,称为“彩虹”解。该方法的关键是避免对指数核进行近似,这是彩虹解的一个显著转换性质使得成为可能。这导致了一种新的用于直接评估欧几里得作用的近似方案,避免了对精确瞬子轨迹的耗时搜索。这种“彩虹近似”可以处理覆盖广泛频率范围的模式的耦合,并弥合绝热和突发近似之间的差距。它比这些传统的近似受到的限制要少得多,并且对于具有强耦合的系统,例如氢键中的质子转移,证明特别有效。与二维模型中已知的精确瞬子作用进行比较,并将其应用于丙二醛的两种同位素的零能级隧穿分裂,以显示该方法的准确性和效率。