Troubleshooting Time-Dependent Density-Functional Theory for Photochemical Applications: Oxirane

The development of analytic-gradient methodology for excited states within conventional time-dependent density-functional theory (TDDFT) would seem to offer a relatively inexpensive alternative to better established quantum-chemical approaches for the modeling of photochemical reactions. However, even though TDDFT is formally exact, practical calculations involve the use of approximate functionals, in particular the TDDFT adiabatic approximation, whose use in photochemical applications must be further validated. Here, we investigate the prototypical case of the symmetric CC ring opening of oxirane. We demonstrate by direct comparison with the results of high-quality quantum Monte Carlo calculations that, far from being an approximation on TDDFT, the Tamm-Dancoff approximation (TDA) is a practical necessity for avoiding triplet instabilities and singlet near instabilities, thus helping maintain energetically reasonable excited-state potential energy surfaces during bond breaking. Other difficulties one would encounter in modeling oxirane photodynamics are pointed out but none of these is likely to prevent a qualitatively correct TDDFT/TDA description of photochemistry in this prototypical molecule.Comment: 19 pages, 17 figures, submitted to the Journal of Chemical Physic

Theoretical Investigation of the Infrared Spectra of the H 5

Reduced dimensional quantum dynamics calculations of the infrared spectrum of the H5 + and D5 + clusters are reported in both low, 300-2200 cm-1, and high, 2400-4500 cm -1, energy regions. The proposed four-dimensional quantum model describes the motion of the proton between the two vibrating hydrogen molecules. The simulations are performed using time-dependent and time-independent approaches within the multiconfiguration time-dependent Hartree method. Propagation of the wavepackets includes an absorbing scheme to deal with vibrational dissociating states, and to assign the different spectral lines, block improved relaxation computations are performed for both bound and predissociative vibrational states of the systems. The reported computations make use of an analytical ab initio-based potential energy, and >on the fly> DFT dipole moment surfaces. The predominant features in the spectra are assigned to the excitations of the shared-proton stretch mode, and above dissociation the symmetric and antisymmetric stretching of the two H2 and the breathing mode of H3 + are also involved. The computed infrared absorption spectra for both cations are in very good agreement with the recent experimental measurements available from multiple-photon dissociation and mass-selected single-photon photodissociation spectroscopy techniques. Comparison of the present results with previous theoretical calculations on these systems is also presented. Such comparisons between different theoretical approaches and experimental measurements can serve to evaluate the approximations employed, and to guide higher-order computations. © 2013 American Chemical Society.This work has been supported by the MICINN grants FIS2010- 18132 and FIS2011-29596-C02-01, Consolider-Ingenio 2010 Programme CSD2009-00038 (MICINN), and COST Action CM1002 (CODECS).Peer Reviewe

Suitable coordinates for quantum dynamics: Applications using the multiconfiguration time-dependent Hartree (MCTDH) algorithm

International audienc

Different Flavors of Nonadiabatic Molecular Dynamics

The Born‐Oppenheimer approximation constitutes a cornerstone of our understanding of molecules and their reactivity, partly because it introduces a somewhat simplified representation of the molecular wavefunction. However, when a molecule absorbs light containing enough energy to trigger an electronic transition, the simplistic nature of the molecular wavefunction offered by the Born‐Oppenheimer approximation breaks down as a result of the now non‐negligible coupling between nuclear and electronic motion, often coined nonadiabatic couplings. Hence, the description of nonadiabatic processes implies a change in our representation of the molecular wavefunction, leading eventually to the design of new theoretical tools to describe the fate of an electronically‐excited molecule. This Overview focuses on this quantity—the total molecular wavefunction—and the different approaches proposed to describe theoretically this complicated object in non‐Born‐Oppenheimer conditions, namely the Born‐Huang and Exact‐Factorization representations. The way each representation depicts the appearance of nonadiabatic effects is then revealed by using a model of a coupled proton–electron transfer reaction. Applying approximations to the formally exact equations of motion obtained within each representation leads to the derivation, or proposition, of different strategies to simulate the nonadiabatic dynamics of molecules. Approaches like quantum dynamics with fixed and time‐dependent grids, traveling basis functions, or mixed quantum/classical like surface hopping, Ehrenfest dynamics, or coupled‐trajectory schemes are described in this Overview