Efficient time-dependent vibrational coupled cluster computations with time-dependent basis sets at the two-mode coupling level: full and hybrid TDMVCC[2]

TL;DR

This paper presents an efficient implementation of time-dependent vibrational coupled cluster theory with basis functions that adapt over time, enabling accurate and scalable simulations of molecular vibrational dynamics.

Contribution

It introduces a computationally efficient TDMVCC[2] method with hybrid approaches, reducing scaling and improving convergence for large molecular systems.

Findings

Cubic scaling with degrees of freedom achieved for wave function propagation.

Hybrid TDMVCC/TDH approach lowers computational cost.

Faster convergence demonstrated on benzoic acid vibrational redistribution.

Abstract

The computation of the nuclear quantum dynamics of molecules is challenging, requiring both accuracy and efficiency to be applicable to systems of interest. Recently, theories have been developed for employing time-dependent basis functions (denoted modals) with vibrational coupled cluster theory (TDMVCC). The TDMVCC method was introduced along with a pilot implementation, which illustrated good accuracy in benchmark computations. In this paper we report an efficient implementation of TDMVCC covering the case where the wave function and Hamiltonian contain up to two-mode couplings. After a careful regrouping of terms, the wave function can be propagated with a cubic computational scaling with respect to the number of degrees of freedom. We discuss the use of a restricted set of active one-mode basis functions for each mode, as well as two interesting limits: i) the use of a full active basis where the variational modal determination amounts essentially to the variational determination of a time-dependent reference state for the cluster expansion; and ii) the use of a single function as active basis for some degrees of freedom. The latter case defines a hybrid TDMVCC/TDH approach which can obtain an even lower computational scaling. The resulting computational scaling for hybrid and full TDMVCC[2] are illustrated for polyaromatic hydrocarbons (PAHs) with up to 264 modes. Finally, computations on the internal vibrational redistribution on benzoic acid (39 modes) are used to show the faster convergence of TDMVCC/TDH hybrid computations towards TDMVCC compared to simple neglect of some degrees of freedom.