A matrix product state formulation of the multiconfiguration time-dependent Hartree (MPS-MCTDH) theory is presented. The Hilbert space that is spanned by the direct products of the phonon degree of freedoms, which is linearly parameterized in the MCTDH ansatz and thus results in an exponential increase in the computational cost, is parametrized by the MPS form. Equations of motion based on the Dirac-Frenkel time-dependent variational principle is derived by using the tangent space projection and the projector-splitting technique for the MPS, which have been recently developed. The mean-field operators, which appear in the equation of motion of the MCTDH single particle functions, are written in terms of the MPS form and efficiently evaluated by a sweep algorithm that is similar to the density-matrix renormalized group sweep. The efficiency and convergence of the MPS approximation to the MCTDH are demonstrated by quantum dynamics simulations of extended excitonic molecular systems.
A matrix product state formulation of the multiconfiguration time-dependent Hartree (MPS-MCTDH) theory is presented. The Hilbert space that is spanned by the direct products of the phonon degree of freedoms, which is linearly parameterized in the MCTDH ansatz and thus results in an exponential increase in the computational cost, is parametrized by the MPS form. Equations of motion based on the Dirac-Frenkel time-dependent variational principle is derived by using the tangent space projection and the projector-splitting technique for the MPS, which have been recently developed. The mean-field operators, which appear in the equation of motion of the MCTDH single particle functions, are written in terms of the MPS form and efficiently evaluated by a sweep algorithm that is similar to the density-matrix renormalized group sweep. The efficiency and convergence of the MPS approximation to the MCTDH are demonstrated by quantum dynamics simulations of extended excitonic molecular systems.
A numerical integrator is proposed for solving the multiconfiguration time-dependent Hartree (MCTDH) equations of motion, which are widely used in computations of molecular quantum dynamics. In contrast to existing integrators, the proposed algorithm does not require inverses of ill-conditioned density matrices and obviates the need for their regularization, allowing for large stepsizes also in the case of near-singular density matrices. The nonlinear MCTDH equations are split into a chain of linear differential equations that can all be efficiently solved by Lanczos approximations to the action of Hermitian-matrix exponentials, alternating with orthogonal matrix decompositions. The integrator is an extension to the Tucker tensor format of recently proposed projector-splitting integrators for the dynamical low-rank approximation by matrices and tensor trains (or matrix product states). The integrator is time-reversible and preserves both the norm and the total energy.
Multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) theory is a rigorous and powerful method to simulate quantum dynamics in complex many-body systems. This approach extends the original MCTDH theory of Meyer, Manthe, and Cederbaum to include dynamically contracted layers in a recursive way, within which the equations of motion are determined from the Dirac-Frenkel variational principle. This paper presents the general derivation of the theory and analyzes the important features that make the ML-MCTDH method numerically efficient. Furthermore, we discuss the generalization of the theory to treat many-body identical particles (fermions or bosons) as well as calculating energy eigenstates via the improved relaxation method.
The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n -mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n -body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH 2 O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C 2 n H 2 n +2 . The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.
The excited-state dynamics of organic molecules, molecular aggregates, and donor-acceptor clusters is typically governed by the interplay of electronic excitations and, due to their flexibility and soft bonding, by the interaction with their vibrations. This interaction in these systems can be characterized by a few relevant electronic states that are coupled to numerous vibrational normal modes, encompassing a vast configurational space of the molecules. The full quantum simulation of these type of systems has been long dominated by the multiconfiguration time-dependent Hartree (MCTDH) approach and its multilayer variants, which are considered the gold standard in the presence of electron-vibration coupling with a large number of modes. Recently, also the matrix product state ansatz (MPS) with appropriate time-evolution schemes has been applied to these types of Hamiltonians. In this article, we provide a numerical comparison of excited-state dynamics between the MCTDH and MPS approaches for two electron-vibration coupled systems. Notably, we consider two models for exciton dissociation at a P3HT:PCBM heterojunction, comprising two electronic states and 100 vibrational modes, and 26 electronic states and 113 vibrational modes, respectively. While both methods agree very well for the first model, more pronounced deviations are found for the second model. We trace back the divergence between the methods to the different way entanglement is treated.
The similarities between Hartree-Fock (HF) theory and the density matrix renormalization group (DMRG) are explored. Both methods can be formulated as the variational optimization of a wave-function Ansatz. Linearization of the time-dependent variational principle near a variational minimum allows to derive the random phase approximation (RPA). We show that the nonredundant parameterization of the matrix product state (MPS) tangent space [J. Haegeman, J. I. Cirac, T. J. Osborne, I. Pi\ifmmode \check{z}\else \v{z}\fi{}orn, H. Verschelde, and F. Verstraete, Phys. Rev. Lett. 107, 070601 (2011)] leads to the Thouless theorem for MPS, i.e., an explicit nonredundant parameterization of the entire MPS manifold, starting from a specific MPS reference. Excitation operators are identified, which extends the analogy between HF and DMRG to the Tamm-Dancoff approximation (TDA), the configuration interaction (CI) expansion, and coupled cluster theory. For a small one-dimensional Hubbard chain, we use a CI-MPS Ansatz with single and double excitations to improve on the ground state and to calculate low-lying excitation energies. For a symmetry-broken ground state of this model, we show that RPA-MPS allows to retrieve the Goldstone mode. We also discuss calculations of the RPA-MPS correlation energy. With the long-range quantum chemical Pariser-Parr-Pople Hamiltonian, low-lying TDA-MPS and RPA-MPS excitation energies for polyenes are obtained.
Dynamical electronic- and vibrational-structure theories have received a growing interest in the past few years due to their ability to simulate spectra recorded with ultrafast experimental techniques. The exact time evolution of a molecular system can, in principle, be obtained from the time-dependent version of full configuration interaction. Such an approach is, however, limited to few-atom systems due to the exponential increase of its cost with the system dimension. In the present work, we overcome this unfavorable scaling by employing the time-dependent density matrix renormalization group (TD-DMRG) which parametrizes the time-dependent wave function as a matrix product state. The time-dependent Schrödinger equation is then integrated with a sweep-based algorithm, as in standard time-independent DMRG. Unlike other TD-DMRG approaches, the one presented here leads to a set of coupled equations that can be integrated exactly. The resulting theory enables us to study real- and imaginary-time evolutions of Hamiltonians comprising more than 20 degrees of freedom that are challenging for current state-of-the-art quantum dynamics algorithms. We apply our algorithm to the simulation of quantum dynamics of models of increasing complexity, ranging from simple excitonic Hamiltonians to more complex ab initio vibronic ones.