Multi-configuration time-dependent Hartree (MCTDH) is an approach to quantum molecular dynamics, an algorithm to solve the time-dependent Schrödinger equation for multidimensional dynamical systems consisting of distinguishable particles. The nuclei of molecules is one example of such particles and their vibrational motion is a form of time-dependence. The method uses an overall wavefunction composed of products of single-particle wavefunctions as first proposed by Douglas Hartree in 1927. The "multiconfiguration" part of the method refers to combining multiple such products.
MCTDH can predict the motion of the nuclei of a molecular system evolving on one or several coupled electronic potential energy surfaces. It is an approximate method whose numerical efficiency decreases with growing accuracy.
MCTDH is suited for multi-dimensional problems, in particular for problems that are difficult or even impossible to solve in conventional ways.
Where the number of configurations is given by the product . The single particle functions (SPFs), , are expressed in a time-independent basis set:
Where is a primitive basis function, in general a Discrete Variable Representation (DVR) that is dependent on coordinate . If , one returns to the Time Dependent Hartree (TDH) approach. In MCTDH, both the coefficients and the basis function are time-dependent and optimized using the variational principle.
Where:
Which is subject to the boundary conditions . After integration, one obtains:
Where only the time derivative is to be varied. We can rewrite this norm squared term as a scalar product, and vary the bra and ket side of the product:
If each variation of is an allowed variation, then both the Lagrangian and the McLanchlan Variational Principle turn into the Dirac-Frenkel Variational Principle:
Which simplest and thus preferred method of deriving the equations of motion.
The original ansatz of MCTDH generates a single layer tensor tree; however, there is a limit to the size and complexity this single layer can handle. This prompted the development of a multilayer (ML)-MCTDH ansatz by Manthe which was then generalized by Vendrell and Meyer.
Multiple layers are generated through the creation of a tensor tree of nodes linking the modes (DOFs). Solving the tree layout is an NP-hard problem, but strides have been taken to automate this process through mode correlations by Mendive-Tapia.
The generalized ML expansion of Meyer can be written as follows:
Where the coordinates are combined as
Where the equations of motion are now represented as follows:
The SPF EOMs are formally defined the same for all layers:
Where is a Hermitian gauge operator defined as follows:
The first verification of the MCTDH method was with the NOCl molecule. Its size and asymmetry makes it a perfect test bed for MCTDH: it is small and simple enough for its numerics to be manually verified, yet complicated enough for it to already squeeze advantages against conventional product-basis methods.
The solvation of the hydronium ion is a topic of continued research. Researchers have been able to successfully use MCTDH to model the Zundel and Eigen ions in close agreement with experiment.
For a typical input in ML-MCTDH to be run, a node tree, potential energy surface, and equations of motion must be generated by the user. These prerequisitesâÂÂalong with total compute timeâÂÂsoft-cap the size of systems able to be studied with ML-MCTDH; however, advances in neural networks have been shown to address the difficulty of the generation of potential energy surfaces. These issues can also by circumvented by using the spin-boson or other similar bath models that do not pose the same assignment challenges.