The multiconfigurational time-dependent Hartree (MCTDH) approach facilitates multidimensional quantum dynamics calculations by representing the wavepacket in an optimal set of time-dependent basis functions, called single-particle functions. Choosing these single-particle functions to be themselves multidimensional wavefunctions which are represented using a MCTDH representation, a multilayer MCTDH scheme has been constructed and used for quantum dynamics calculations treating up to 1000 degrees of freedom rigorously [Wang and Thoss, J. Chem. Phys. 199, 1289 (2003)]. The present work gives a practical scheme which facilitates the application of the multilayer MCTDH approach, which previously has only been employed to study systems described by model-type Hamiltonians, to molecular systems described by more complicated Hamiltonians and general potential energy surfaces. A multilayer extension of the correlation discrete variable representation (CDVR) scheme employed in MCTDH calculations studying quantum dynamics on general potential energy surfaces is developed and tested in a simple numerical application. The resulting multilayer MCTDH/CDVR approach might offer a perspective to rigorously describe the quantum dynamics of larger polyatomic systems.

Thermal rate constants are calculated for the H + CH 4 f CH 3 + H 2 reaction employing the potential energy surface of Espinosa-García (Espinosa-García, J. J. Chem. Phys. 2002, 116, 10664). Two theoretical approaches are used. First, we employ the multiconfigurational time-dependent Hartree method combined with flux correlation functions. In this way rate constants in the range 225-400 K are obtained and compared with previous results using the same theoretical method but the potential energy surface of Wu et al. (Wu, T.; Werner, H.-J.; Manthe, U. Science 2004, 306, 2227). It is found that the Espinosa-García surface results in larger rate constants. Second, a harmonic quantum transition state theory (HQTST) implementation of instanton theory is used to obtain rate constants in a temperature interval from 20 K up to the crossover temperature at 296 K. The HQTST estimates are larger than MCTDH ones by a factor of about three in the common temperature range. Comparison is also made with various tunneling corrections to transition state theory and quantum instanton theory.

A full-dimensional quantum dynamics simulation of a hydrogen atom reacting with methane on an accurate ab initio potential energy surface is reported. Based on first-principles theory, thermal rate constants are predicted with an accuracy comparable to (or even exceeding) experimental precision. The theoretical prediction is within the range of the significantly varied experimental rate constants reported by different groups. This level of accuracy has previously been achieved only for smaller, three-or four-atom reactive systems. Comparison with classical transition state theory confirms the importance of quantum mechanical tunneling for the rate constant below 400 kelvin.