Gaussian-based multiconfiguration time-dependent Hartree: A two-layer approach. II. Application to vibrational energy transport in a molecular chain.

We report on first applications of the Two-Layer Gaussian-based Multi-Configuration Time-Dependent Hartree (2L-GMCTDH) method [Römer et al., J. Chem. Phys. 138, 064106 (2013)] for high-dimensional quantum propagation using variational Gaussian basis sets. This method circumvents the limitations of conventional variational Gaussian wavepacket (GWP) methods by introducing a hierarchical wavefunction representation with a fully flexible first layer composed of orthogonal single-particle functions, which are in turn expressed as superpositions of GWPs of fixed width. The method is applied to a model Hamiltonian describing vibrational energy transport through a molecular chain. The model combines bilinear site-to-site couplings with site-local couplings induced by cubic anharmonicities. We report on simulation results for realizations comprising 5 sites with 35 vibrational modes and 18 sites with 90 vibrational modes, which are shown to be in excellent agreement with reference calculations by the Multi-Layer MCTDH method.