The correlation discrete variable representation revisited

TL;DR

This paper revisits the correlation discrete variable representation (CDVR) to improve quantum dynamics calculations, introducing a revised approach that enhances efficiency and accuracy without increasing computational costs.

Contribution

The authors present a revised non-hierarchical CDVR method that avoids explicit projection, scales favorably with system size, and uses artificial SPFs to systematically improve quadrature accuracy.

Findings

Revised CDVR achieves n^4 scaling with the number of SPFs.

Demonstrates high accuracy in photodissociation and vibrational state calculations.

No increased CPU time for 24D pyrazine system compared to sum of products methods.

Abstract

The correlation discrete variable representation (CDVR) enables efficient quantum dynamics calculation with the multi-layer multi-configurational time-dependent Hartree (MCTDH) approach on general potential energy surfaces. It employs a time-dependent quadrature to compute potential energy matrix elements, thereby eliminating the need to refit the potential to a sum of products form. The non-hierarchical CDVR conserves the inherent symmetry properties of tree-shaped wavefunction representations and drastically reduces the number of grid points compared to the original hierarchical CDVR. However, it requires projection on the space spanned by the single-hole functions (SHFs) at each node of the tree, which can introduce unphysical couplings for unconverged basis sets. In this work, the non-hierarchical CDVR is revisited and a revised approach that avoids explicit projection on the single-hole space is introduced. The computational costs of the revised approach scale favorably with the number of single-particle functions (SPFs): for a tree with three edges at each node and $n$ SPFs at each edge, a n^4 scaling is achieved. Furthermore, a revised scheme that uses artificial SPFs to systematically increase the accuracy of the CDVR quadrature is presented. Computations studying the photodissociation of NOCl, the vibrational states of methyl, and the non-adiabatic quantum dynamics of photoexcited pyrazine demonstrate the accuracy and efficiency of the revised non-hierarchical CDVR. Notably, for the 24-dimensional pyrazine system the use of the CDVR does not increase the required CPU time compared to calculations utilizing the sum of products form of the vibronic coupling model.