MCTDH
Content
 Brief Description of MCTDH
 MCTDH Applications
 The Heidelberg MCTDH Package
 Literature Downloads
 Benchmark Downloads
 Miscellanies
 List of MCTDH Publications
 MCTDH documentation (V8.5), (V8.4), (V8.3). (Not all links of the documentation will work.
If you are an MCTDH user, please use the documentation which comes with the package).
Important Note
There shall be a "High dimensional quantum dynamics: challenges and opportunities" workshop in Rostock, second half of 2016. The workshop will be organized by Oliver Kühn.
There has been a "High dimensional quantum dynamics: challenges and opportunities" workshop in Mittelwihr near Strasbourg, September 2  6, 2014.
See: http://quantique.ustrasbg.fr/HDQD2014/
There has been a "High dimensional quantum dynamics: challenges and opportunities" workshop in Birmingham, April 1214, 2012, organized by Graham Worth.
See: http://www.stchem.bham.ac.uk/~worthgrp/hdqd2012/index.html
There has been a MCTDH workshop in Montpellier, February 2528, 2008, organized by Fabien Gatti.
http://www.lsd.univmontp2.fr/Workshop_quantum_dynamics/
April 15th, 2009 there has appeared the book:
Multidimensional Quantum Dynamics:
MCTDH Theory and Applications
H.D. Meyer, F. Gatti, and G. A. Worth, editors
For a preview (cover page, list of contents, and introduction) see here (PDF).
For an announcement and order form click here (PDF).
You can view the book at books.google.com .
For further reading go to Literature Downloads. There you will find, among others, the MCTDH review and my lecture notes on MCTDH (intro_MCTDH).
Brief Description of MCTDH
MCTDH stands for Multi Configuration Time Dependent Hartree. MCTDH is a general algorithm to solve the timedependent Schrödinger equation for multidimensional dynamical systems consisting of distinguishable particles. MCTDH can thus determine the quantal motion of the nuclei of a molecular system evolving on one or several coupled electronic potential energy surfaces. MCTDH by its very nature is an approximate method. However, it can be made as accurate as any competing method, but its numerical efficiency deteriorates with growing accuracy.
MCTDH is designed for treating multidimensional problems, in particular problems that are difficult or even impossible to attack in a conventional way. There is no or only little gain when treating systems with less than three degrees of freedom by MCTDH. However, for convenience − not for numerical speed − be have used the MCTDH package even for onedimensional problems. MCTDH will in general be best suited for systems with 4 to 12 degrees of freedom. Because of hardware limitations it may in general not be possible to treat much larger systems. For a certain class of problems, however, one can go much further. Already some time ago (1998/1999) we have performed converged MCTDH calculations on the absorption spectrum of pyrazine accounting for the motion of all 24 (!) internal degrees of freedom [30,35,36]. Similar calculations have been done for the ionisation spectra of Allene [56], Butatrien [57] and Benzene [66]. (The numbers is square brackets refer to the List of MCTDH Publications ). The spinboson model was treated by Haobin Wang including 80 (!!) vibrational modes [49], and we have recently used MCTDH to study the multidimensional HenonHeiles Hamiltonian including up to 32 degrees of freedom [70] and a Morse oscillator coupled to a harmonic bath including up to 61 degrees of freedom [76].
The MCTDH program package has been generalised to enable the propagation of density operators [39,45,46,99], and recently it has been generalized to compute eigenstates ("improved relaxation") [75,128].
MCTDH has the following advantages
 The MCTDH algorithm can be very fast. In fact it may be several orders of magnitude faster than a conventional treatment. Whether MCTDH is fast or not depends on the system under consideration (see below).
 The MCTDH algorithm is small, i.e it requires an amazingly small amount of central memory (RAM). The MCTDH representation of the wave function is also very compact, making it possible to store a large number of wave functions for later analysis.
MCTDH has the following disadvantages
 The MCTDH algorithm requires the Hamiltonian to be expanded as a sum of products of oneparticle operators. The kinetic energy is usually already in this MCTDH form. The potential is in MCTDH form only when dealing with model problems. In general one has to transform [19,31] the potential (obtained e.g. by quantum chemistry calculations) to MCTDH form. There exist an efficient algorithm to accomplish this transformation. However, the MCTDH computation time grows linearly with the number of Hamiltonian terms and MCTDH will become slow, if an accurate representation of the potential requires very many terms.
 The MCTDH algorithm is fast only if the time dependent wavepacket can  at each instant of time  be expanded into a small (but optimised) product basis set. This makes it difficult for MCTDH to propagate a wavepacket when the system is highly chaotic. It is also difficult for MCTDH to propagate a wavepacket for very long times (hundreds of vibrational periods, say). However, we recently [70] have propagated a wavepacket on a 4D HenonHeiles potential over 300 time units (i.e. about 50 vibrational periods). Thus even a chaotic system may be studied for rather long times. The 14D HenonHeiles potential was also investigated. Here the propagation was stopped after 25 time units, as this was enough to reproduce the spectrum.
MCTDH Applications
Problems treated with MCTDH are:
 Photodissociation: absorption spectra and product distributions. (NOCl, NO_{2}, CH_{3}I, ICN, ArHBr, Ar_{2}HBr, HCl/ice) [2,3,4,9,13,71,104,106].
 Computation of absorption (or ionisation) spectra of polyatomic molecules. Vibronically coupled systems. Pyrazine [20,30,35,36]; Allene [56]; Butatriene [57]; Benzene [66], Furan [93], Pentatetraene [114], Cyclopropane [135], .Phenide [139], Difluorobenzene [171].
 Computation of resonance Raman spectra of polyatomic molecules. (CH_{3}I) [25].
 Vibrational predissociation. Cl_{2}Ne [12] ; I_{2}...Ne_{2} cluster [60].
 Moleculesurface scattering : rotational excitation and diffractive scattering of diatomic molecules colliding with a hard but corrugated surface. (H_{2}/LiF, N_{2}/LiF) [17,23,50]. (See also [42,43]). Dissociation from a surface. (CH_{4}/Ni, CH_{3}I/MgO) [11,16,28]. Dissociative adsorption of H_{2} on Cu(100) [92], N2 on stepped Ru(0001) [105], and H_{2} on Pt(111) [115].
 Inelastic scattering scattering, cross−sections: H_{2} + H_{2} [112,148,157,187].
 Reactive scattering, cross−sections: H + H_{2} > H_{2} + H , H + D_{2} > HD + D [15,29,34,54,65]. H_{2} + OH > H_{2}O + H [195].
 Reactive scattering, thermal reaction rates: [24,27,33,37,47,48,62,67,74,124,140,155,193].
 Computation of eigenenergies by combining MCTDH with the filtering method. CO_{2} [53], HO_{2} [55], toluene [88], and HFCO [118].
 SpinBoson model [49,58,59,77].
 Nuclear dynamics during electronscattering processes: CO_{2} + e^{− } [62,70], H_{2}O + e^{− } [86,127,133].
 Proton transfer reactions [63,80,98,100,102].
 Isomerisation of HONO [83,84] and IVR of CF_{3}H [85], Toluene [88], HFCO [118] and DFCO [134].
 Highdimensional HenonHeiles potentials [70].
 Morse oscillator coupled to a bath of 60 oscillators [76].
 Small bosonic systems [129,130,141,142,158,161,162,178,197]
 Optimal control [126,161].
 Computation of eigenstates with improved relaxation [128,174].
 Density operator propagation on model systems (up to 4D) [39,45,46] and on vibrational relaxation of CO on Cu (6D) [99].
 Dynamics and infrared spectra of the Zundel cation H_{5}O_{2}^{+} [143,152,153,169,176].
There are more applications. See the List of MCTDH Publications for a comprehensive overview. (The citations [...] refer to this list.)
The Heidelberg MCTDH Package
The Heidelberg MCTDH package is a set of programs for multidimensional quantum dynamics, and can do much more than wavefunction (or density operator) propagation using the MCTDH algorithm. For example, numerically exact propagations are also possible using a shortiterative Lanczos integrator. As a byproduct of the improved relaxation method, it is also possible to generate a desired eigenfunction of an operator, and for smaller systems the spectrum of a Hamiltonian may be obtained by Lanczos diagonalisation of a Hamiltonian in a DVR basis. The package consists of about 60 programs, the largest and most important of which is called simply mctdh. Using keyword input read from text files, it is able to set up a system using a range of DVRs or FFT for the primitive basis, form an initial wavefunction (density operator), and propagate this wavefunction (density operator) in time. Perhaps one of the most powerful features of the program is that it uses a text input to generate the operator. In fact, if the operator has a simple analytic form, it is often possible to implement it without having to touch the code. Routines coding for more complicated functions can also be linked to the program. If potential functions are not in MCTDH form, there is the potfit program to make the transformation. It is also possible to use the CDVR method, or even to use the potential as it is, of course with resulting inefficient propagation. Timedependent Hamiltonians can also be used. Other important programs of the package are a set of analysis tools. These include filter, which performs a filter analysis of the autocorrelation function, and flux, which does the flux analysis. Various programs plot one and twodimensional graphs of the wavefunction and the potential energy surface, and simple movies can be made. Other programs can be used to check the convergence of a calculation, generate a spectrum from the autocorrelation function, etc. All plotting uses the freely available Gnuplot program, often driven using interactive menus. The package consists of more than 250,000 lines of code (mainly FORTRAN 77, some C, and including documentation and examples). The documentation consists of both online documentation (in HTML) and the MCTDH User's Guide (in Latex). The installation is performed by convenient installation scripts. We have run MCTDH on DEC−alpha, IBM−RS6000, Cray, Sun, Silicon Graphics and HP computers, and in particular on Linux−PC's.
The MCTDH program package is distributed on request to interested researchers. If you want to work with MCTDH, please send an email to
HansDieter.Meyer at pci.uniheidelberg.de
Literature Downloads
The first description of the MCTDH algorithm was given in [1] (see the list of MCTDH publications). A reader interested in the MCTDH algorithm should start with reading [2]. The ref. [26] discusses all technical developments up to mid 1997. A comprehensive overview on MCTDH can be found in [40]. This review article, published in Physics Reports, is the most detailed description of MCTDH available. It does not only describe MCTDH as such, but also some of the analyse techniques, as well as some technical aspects (e.g. DVR's). The preprint version (LaTeX/postscript) of this review can be downloaded (see below).
More recently we wrote a feature article on MCTDH [72] which appeared in Theoretical Chemistry Accounts (TCA). This article covers all the new developments (e.g. density operator propagation) of MCTDH.
Some small errors have been found in the review published as Physics Reports 324, 1 (2000). A correction list can be downloaded from below.
(The download version of the review is already corrected).
A brief description of the MCTDH method was prepared for an CCP6 workshop (March 2001). This writeup can be downloaded from below.
In 2008 there appeared two review articles, one on multimode nonadiabatic dynamics (Worth et al) and one one computation of vibrational energies (Bowman et al), which you may find interesting. The preprint form of these articles can be downloaded from below.
In May 2009 there will appear a book on MCTDH:
Multidimensional Quantum Dynamics: MCTDH Theory and Applications The first few pages of the book can be downloaded from below.
The notes of a 2010 lecture "Introduction to MCTDH" can also be found below. These notes have recently been corrected (June 2013).
To get an impression of the MCTDH project see the MCTDH documentation. (Not all links will work! If you are a MCTDH user, please use the HTMLdocumentation which comes with the package.)
There is also the MCTDH Users Guide. In contrast to the review, which describes the MCTDH algorithm, the guide describes the use of the Heidelberg MCTDH program package.
Some sets of slides from talks at the MCTDH workshop July 2001 are available: The MCTDH Program: Structure and Development and Operators and the MCTDH Program
Finally, a MCTDH bibtex file and the results of some benchmark calculations can be downloaded from this site.
The following files are for download  

gzipped Postscript  Description  
CCP6, 116kB  CCP6, 110kB  Short description of MCTDH, prepared for an CCP6 workshop (2001) 
Review, 460kB  Review, 708kB  MCTDH Review, preprint version, published in: Physics Reports 324 (2000), 1. 
corrections, 44kB  corrections, 28kB  List of corrections for the review. (The download version is already corrected). 
guide85, 484kB  guide85, 964kB  MCTDH Users Guide. Description of the use of the Heidelberg MCTDH Package (8.5 branch). 
guide84, 484kB  guide84, 964kB  MCTDH Users Guide. Description of the use of the Heidelberg MCTDH Package (8.4 branch). 
guide83, 956kB  guide83, 480kB  MCTDH Users Guide. Description of the use of the Heidelberg MCTDH Package (8.3 branch). 
feature, 200kB  feature, 252kB  MCTDH Feature Article, preprint version, published in: Theor. Chem. Acc. 109 (2003), 251 
conical intersection, 1.3 MB  conical intersection, 2.4 MB  Multidimensional dynamics involving a conical intersection: Wavepacket calculations using the MCTDH method Worth, Meyer, and Cederbaum, in: Conical Intersections, W. Domcke, D. R. Yarkony, and H. Köppel, Eds., World Scientific Co. (2004) ISBN 9812386726 
MCTDHbook  Multidimensional Quantum Dynamics: MCTDH Theory and Applications H.D. Meyer, G. A. Worth, and F. Gatti, editors WileyVCH, Weinheim, April 2009. First few pages. 

bowman review  Variational Quantum Approaches for Computing Vibrational Energies of Polyatomic Molecules J. M. Bowman, T. Carrington, and H.D. Meyer Molecular Physics 106, 21452182 (2008). 

IRPCreview  Using the MCTDH wavepacket ptopagation method to describe multimode nonadiabatic dynamics G. A. Worth, H.D. Meyer, H. Köppel, L.S. Cederbaum and I. Burghardt International Reviews in Physical Chemistry, Vol. 27, No. 3, (2008) 569606. (Taylor and Francis). 

program, 94kB  program, 127kB  Slides (MCTDH workshop): The MCTDH Program: Structure and Development. 
operators, 128kB  operators, 183kB  Slides (MCTDH workshop): Operators and the MCTDH Program. 
intro_MCTDH, 1.7MB  Introduction to MCTDH, lecture notes (HDM 2010), written by DPR (2011). Corrected (2013).  
meyer_rev_2011, 1.1MB 
Studying molecular quantum dynamics with the multiconfiguration timedependent Hartree (MCTDH) method . H.D. Meyer, 2012. 

mctdh.bib.gz, 12kB  mctdh.bib, 32kB  Bibtex file "mctdh.bib" containing references to MCTDH articles. 
MCTDH benchmarks  On this directory one finds the data of some MCTDH benchmark calculations. (targz) 
Miscellanies
The MCTDH program package is distributed on request to interested researchers. If you want to work with MCTDH, please send an email to
HansDieter.Meyer at pci.uniheidelberg.de
We have recently investigated the multidimensional HenonHeiles system with MCTDH (J.Chem.Phys. 117 (2002), 10499). As the results of these benchmark calculations may be of interest to others (e.g. for checking their method against converged fully quantal results), the data may be downloaded from MCTDH benchmarks. There you will find also the results of the Morse oscillator coupled to a bath calculation (J.Chem.Phys. 119 (2003), 24) and of the two pyrazine calculations.
If you are writing an article, the MCTDH bibtex file may be of some help. The bibtex file can be found under Literature Downloads
Prof. Dr. H.D. Meyer, Universität Heidelberg, Im Neuenheimer Feld 229
email: HansDieter.Meyer at pci.uniheidelberg.de
Phone: +49 − 6221 − 54 52 10