In this paper effects of higher order Jahn-Teller coupling terms on the nonadiabatic dynamics are studied. Of particular interest is the case when the potential energy surfaces of the degenerate state show pronounced anharmonicity. In order to demonstrate the effects a two-dimensional E multiply sign in circle e Jahn-Teller model system is treated which is based on the e(') stretching vibration of the photoactive (2)E(') state of NO(3) as a realistic example. The sixth order E multiply sign in circle e Jahn-Teller Hamiltonian is derived in the diabatic representation which is valid for any system with a C(3) rotation axis. This diabatization scheme is compared to lower-order Jahn-Teller Hamiltonians and to symmetry adapted as well as ad hoc approximations. Lower-order potentials result in pronounced quantitative and qualitative differences in the dynamics, including differences in the evolution of mean values, the autocorrelation functions (and thus the corresponding spectra), and the electronic population evolution. In the particular example treated, the results of fourth and fifth order potentials are very similar to the sixth order reference system. In contrast, the approximate sixth order Hamiltonians, though the corresponding adiabatic surfaces seem to be nearly identical, results in pronounced differences. The possible consequences for the dynamics of realistic systems with higher dimensionality are briefly discussed.

An analytical anharmonic six-dimensional three-sheeted potential energy surface for the ground and first excited states of the ammonia cation has been developed which is tailored to model the ultrafast photoinduced dynamics. Selected ab initio cuts, obtained by multireference configuration interaction calculations, have been used to determine the parameters of a diabatic representation for this Jahn-Teller and pseudo-Jahn-Teller system. The model includes higher-order coupling terms both for the Jahn-Teller and for the pseudo-Jahn-Teller matrix elements. The relaxation to the ground state is possible via dynamical pseudo-Jahn-Teller couplings involving the asymmetric bending and stretching coordinates. The photoelectron spectrum of NH3 and the internal conversion dynamics of NH3+ have been determined by wave packet propagation calculations employing the multiconfigurational time-dependent Hartree method. Three different time scales are found in the dynamics calculations for the second absorption band. The ultrafast Jahn-Teller dynamics of the two excited states occurs on a 5 fs time scale. The major part of the internal conversion to the ground state takes place within a short time scale of 20 fs. This fast internal conversion is, however, incomplete and the remaining excited state population does not decay completely even within 100 fs.

Benchmark calculations of the tunneling splitting in malonaldehyde using the full dimensional potential proposed by Yagi et al. [J. Chem. Phys. 115, 10647 (2001)] are reported. Two exact quantum dynamics methods are used: the multiconfigurational time-dependent Hartree (MCTDH) approach and the diffusion Monte Carlo based projection operator imaginary time spectral evolution (POITSE) method. A ground state tunneling splitting of 25.7±0.3 cm−1 is calculated using POITSE. The MCTDH computation yields 25 cm−1 converged to about 10% accuracy. These rigorous results are used to evaluate the accuracy of approximate dynamical approaches, e.g., the instanton theory.