Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of Chemical Theory and Computation

Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials

Kentaro Hino, Yuki Kurashige

Published May 23, 2022• 11 citations
Publisher

Abstract

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.

DOI
10.1021/acs.jctc.2c00243

Related papers

Comparison of Matrix Product State and Multiconfiguration Time-Dependent Hartree Methods for Nonadiabatic Dynamics of Exciton Dissociation

Maximilian F. X. Dorfner, Dominik Brey, Irène Burghardt +1 more

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.

Journal of Chemical Theory and Computation•Oct 4, 2024•10 citations

Matrix product state formulation of the multiconfiguration time-dependent Hartree theory

Yuki Kurashige

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.

arXiv:1808.01110•The Journal of Chemical Physics•Nov 21, 2018•39 citations

Time-dependent density matrix renormalization group quantum dynamics for realistic chemical systems

Xiaoyu Xie, Yuyang Liu, Yao Yao +3 more

Electronic and/or vibronic coherence has been found by recent ultrafast spectroscopy experiments in many chemical, biological, and material systems. This indicates that there are strong and complicated interactions between electronic states and vibration modes in realistic chemical systems. Therefore, simulations of quantum dynamics with a large number of electronic and vibrational degrees of freedom are highly desirable. Due to the efficient compression and localized representation of quantum states in the matrix-product state (MPS) formulation, time-evolution methods based on the MPS framework, which we summarily refer to as tDMRG (time-dependent density-matrix renormalization group) methods, are considered to be promising candidates to study the quantum dynamics of realistic chemical systems. In this work, we benchmark the performances of four different tDMRG methods, including global Taylor, global Krylov, and local one-site and two-site time-dependent variational principles (1TDVP and 2TDVP), with a comparison to multiconfiguration time-dependent Hartree and experimental results. Two typical chemical systems of internal conversion and singlet fission are investigated: one containing strong and high-order local and nonlocal electron-vibration couplings and the other exhibiting a continuous phonon bath. The comparison shows that the tDMRG methods (particularly, the 2TDVP method) can describe the full quantum dynamics in large chemical systems accurately and efficiently. Several key parameters in the tDMRG calculation including the truncation error threshold, time interval, and ordering of local sites were also investigated to strike the balance between efficiency and accuracy of results.

arXiv:1908.10588•The Journal of Chemical Physics•Dec 9, 2019•75 citations

Computing vibrational eigenstates with tree tensor network states (TTNS)

Henrik R. Larsson

We present how to compute vibrational eigenstates with tree tensor network states (TTNSs), the underlying ansatz behind the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. The eigenstates are computed with an algorithm that is based on the density matrix renormalization group (DMRG). We apply this to compute the vibrational spectrum of acetonitrile (CH3CN) to high accuracy and compare TTNSs with matrix product states (MPSs), the ansatz behind the DMRG. The presented optimization scheme converges much faster than ML-MCTDH-based optimization. For this particular system, we found no major advantage of the more general TTNS over MPS. We highlight that for both TTNS and MPS, the usage of an adaptive bond dimension significantly reduces the amount of required parameters. We furthermore propose a procedure to find good trees.

arXiv:1909.13831•The Journal of Chemical Physics•Nov 25, 2019•77 citations

n -Mode Quantized Anharmonic Vibronic Hamiltonians for Matrix Product State Dynamics

Valentin Barandun, Nina Glaser, Markus Reiher

High Resolution Image Download MS PowerPoint Slide Theoretical predictions of photochemical processes are essential for interpreting and understanding spectral features. Reliable quantum dynamics calculations of vibronic systems require precise modeling of anharmonic effects in the potential energy surfaces and off-diagonal nonadiabatic coupling terms. In this work, we present the n -mode quantization of all vibronic Hamiltonian terms comprised of general high-dimensional model representations. We expand the existing vibrational DMRG formalism by applying the n -mode quantization to all potential energy surfaces entering the Hamiltonian as well as to all off-diagonal coupling terms. This is accompanied by the introduction of a novel matrix product state architecture employing tailored local site operators, which allow an effective encoding of the vibronic wave function in a tensor-train format. This results in a second-quantized framework for accurate vibronic calculations employing the density matrix renormalization group algorithm. We demonstrate the accuracy and reliability of this approach by calculating the excited-state quantum dynamics of maleimide. We analyze convergence and the choice of parameters of the underlying time-dependent density matrix renormalization group algorithm for the n -mode vibronic Hamiltonian, demonstrating that it enables accurate calculations of complex photochemical dynamics.

arXiv:2511.03936•Journal of Chemical Theory and Computation•Feb 12, 2026•0 citations

Thouless theorem for matrix product states and subsequent post density matrix renormalization group methods

Sebastian Wouters, Naoki Nakatani, Dimitri Van Neck +1 more

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.

arXiv:1305.1761•Physical Review B•Aug 12, 2013•46 citations
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.
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.
Electronic and/or vibronic coherence has been found by recent ultrafast spectroscopy experiments in many chemical, biological, and material systems. This indicates that there are strong and complicated interactions between electronic states and vibration modes in realistic chemical systems. Therefore, simulations of quantum dynamics with a large number of electronic and vibrational degrees of freedom are highly desirable. Due to the efficient compression and localized representation of quantum states in the matrix-product state (MPS) formulation, time-evolution methods based on the MPS framework, which we summarily refer to as tDMRG (time-dependent density-matrix renormalization group) methods, are considered to be promising candidates to study the quantum dynamics of realistic chemical systems. In this work, we benchmark the performances of four different tDMRG methods, including global Taylor, global Krylov, and local one-site and two-site time-dependent variational principles (1TDVP and 2TDVP), with a comparison to multiconfiguration time-dependent Hartree and experimental results. Two typical chemical systems of internal conversion and singlet fission are investigated: one containing strong and high-order local and nonlocal electron-vibration couplings and the other exhibiting a continuous phonon bath. The comparison shows that the tDMRG methods (particularly, the 2TDVP method) can describe the full quantum dynamics in large chemical systems accurately and efficiently. Several key parameters in the tDMRG calculation including the truncation error threshold, time interval, and ordering of local sites were also investigated to strike the balance between efficiency and accuracy of results.
We present how to compute vibrational eigenstates with tree tensor network states (TTNSs), the underlying ansatz behind the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. The eigenstates are computed with an algorithm that is based on the density matrix renormalization group (DMRG). We apply this to compute the vibrational spectrum of acetonitrile (CH3CN) to high accuracy and compare TTNSs with matrix product states (MPSs), the ansatz behind the DMRG. The presented optimization scheme converges much faster than ML-MCTDH-based optimization. For this particular system, we found no major advantage of the more general TTNS over MPS. We highlight that for both TTNS and MPS, the usage of an adaptive bond dimension significantly reduces the amount of required parameters. We furthermore propose a procedure to find good trees.
High Resolution Image Download MS PowerPoint Slide Theoretical predictions of photochemical processes are essential for interpreting and understanding spectral features. Reliable quantum dynamics calculations of vibronic systems require precise modeling of anharmonic effects in the potential energy surfaces and off-diagonal nonadiabatic coupling terms. In this work, we present the n -mode quantization of all vibronic Hamiltonian terms comprised of general high-dimensional model representations. We expand the existing vibrational DMRG formalism by applying the n -mode quantization to all potential energy surfaces entering the Hamiltonian as well as to all off-diagonal coupling terms. This is accompanied by the introduction of a novel matrix product state architecture employing tailored local site operators, which allow an effective encoding of the vibronic wave function in a tensor-train format. This results in a second-quantized framework for accurate vibronic calculations employing the density matrix renormalization group algorithm. We demonstrate the accuracy and reliability of this approach by calculating the excited-state quantum dynamics of maleimide. We analyze convergence and the choice of parameters of the underlying time-dependent density matrix renormalization group algorithm for the n -mode vibronic Hamiltonian, demonstrating that it enables accurate calculations of complex photochemical dynamics.
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.