Time evolution of ML-MCTDH wavefunctions. II. Application of the projector splitting integrator

The multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach can suffer from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of regularization of the EOMs. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here, we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to 10⁶ bath modes and find that for these problems, the PSI requires roughly 3-4 orders of magnitude fewer Hamiltonian evaluations and 2-3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH and 2-3/1-2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to ∼1.3×10⁹ parameters to obtain accurate dynamics.