Hierarchical Wavepacket Propagation Framework via ML-MCTDH for Molecular Reaction Dynamics

TL;DR

This paper introduces a hierarchical wavepacket propagation framework using ML-MCTDH for efficient simulation of molecular reaction dynamics, combining mode partitioning, SOP representations, and tensor network techniques.

Contribution

The work develops a novel hierarchical ML-MCTDH framework with tensor network reformulation, enabling scalable and accurate wavepacket propagation for complex molecular systems.

Findings

Hierarchical mode partitioning improves computational efficiency.

SOP representation of KEO and PES facilitates calculations.

Tensor network reformulation offers alternative computational strategies.

Abstract

This work presents a computational framework for studying reaction dynamics via wavepacket propagation, employing the multiconfiguration time-dependent Hartree (MCTDH) method and its multilayer extension (ML-MCTDH) as the core methodologies. The core idea centers on the concept of modes that combine several coordinates along with their hierarchical separations because the degrees of freedom are too numerous to be efficiently treated as a single mode. First, the system is partitioned into several fragments within the same layer, and these fragments are further decomposed. Repeating this process, a hierarchical separation of modes emerges, until modes of a manageable size are achieved. Accordingly, the coordinates frame can be designed hierarchically. Second, the kinetic energy operator (KEO) is derived as a sum-of-products (SOP) of single-particle differential operators through polyspherical approach, while the potential energy surface (PES) is expressed in a similar SOP form of single-particle potentials (SPPs) through (1) reconstruction approaches using an existing PES or (2) direct approaches based on a computed database. Third, the nuclear wave function is expressed in a multi-layer expansion form, where each term is a product of single-particle functions (SPFs) that are further expanded by the SPFs in deeper layer. This expansion form is also adopted by variational eigensolver for electronic wave function. Finally, the Dirac-Frenkel variational principle leads to a set of working equations whose solutions reproduce reaction dynamics. In addition, the hierarchical framework can be rearranged by the mathematical language of tensor network (TN) or tree tensor network (TTN). In this work, we compare the methods represented by function with those in the form of TN or TTN. We also discuss the limitations of the present framework and propose solutions, providing further perspectives.