Non-Gaussian wave packet dynamics in anharmonic potential: Cumulant expansion treatment

Abstract

This manuscript utilizes cumulant expansion as an alternative algebraic approach to evaluating integrals and solving a system of nonlinear differential equations for probing anharmonic dynamics in condensed phase systems using Morse oscillator. These integrals and differential equations become harder to solve as the anharmonicity of the system goes beyond that of Morse oscillator description. This algebraic approach becomes critically important in case of Morse oscillator as it tends to exhibit divergent dynamics and numerical uncertainties at low temperatures. The autocorrelation function is calculated algebraically and compared to the exact one for they match perfectly. It is also compared to the approximate autocorrelation function using the differential equations technique reported in Toutounji (2014) for weak and strong electron–phonon coupling cases. It is found that the present cumulant method is more efficient, and easier to use, than the exact expression. Deviation between the approximate autocorrelation function and the exact autocorrelation function starts to arise as the electron–phonon coupling strength increases. The autocorrelation function obtained using cumulants identically matches the exact autocorrelation function, thereby surpassing the approach presented in Toutounji (2014). The advantage of the present methodology is its applicability to various types of electron–phonon coupling cases. Additionally, the herein approach only uses algebraic techniques, thereby avoiding both the divergence integral and solving a set of linear first- and second-order partial differential equations as was done in previous work. Model calculations are presented to demonstrate the accuracy of the herein work.