Multiconfiguration time-dependent Hartree and classical dynamics studies of the photodissociation of HF and HCl molecules adsorbed on ice: extension to three dimensions.

The 3D photodissociation dynamics of HCl and HF molecules adsorbed on ice is studied by quantum and classical simulations. The quantum calculations are carried out with the multiconfiguration time-dependent Hartree (MCTDH) approach. Dynamical observables like angular distributions in the momentum space of the H fragments, absorption cross sections are computed. The results are compared with our previous 2D studies. As expected, less encapsulation of the H atom between the ice surface and the halogen atom is obtained in the 3D study, resulting in less pronounced interference structures in the photoabsorption cross sections and in a decrease of the classical rainbow peaks observed in the 2D scheme. Although the amplitudes of the oscillations corresponding to quantum interferences in the asymptotic angular distribution of the H fragment are different between the 2D and 3D results, the qualitative pattern of the oscillations is similar in the 2D and 3D approaches. In addition, a good agreement is observed for the angular distribution between the classical and the quantum calculations.

Vibronic transitions and quantum dynamics in molecular oligomers: a theoretical analysis with an application to aggregates of perylene bisimides.

Vibronic absorption spectra of molecular aggregates consisting of up to N = 9 monomer units are calculated employing methods of time-dependent quantum mechanics. Taking one vibrational degree of freedom for each monomer into account and treating one-exciton excited electronic states leads to a problem with N vibrations and N electronically coupled states. The demanding quantum propagation is carried out within the multiconfiguration time-depended Hartree method (MCTDH). Spectral features of and population transfer in the aggregates are analyzed as a function of the aggregate size and the strength of the electronic coupling. With a model for oligomers of perylene bisimides, it is shown how measured temperature-dependent absorption spectra correlate with the aggregate size. Furthermore, the exciton localization and dynamics in these aggregates are investigated.

Multistate vibronic interactions in difluorobenzene radical cations. II. Quantum dynamical simulations.

The multistate vibronic dynamics in the X-D electronic states of all three difluorobenzene radical cations are investigated theoretically by an ab initio quantum dynamical approach. The vibronic coupling scheme and the ab initio values of the system parameters are adopted from Paper I [S. Faraji and H. Koppel, J. Chem. Phys. 129, 074310 (2008)]. Extensive calculations by wave-packet propagation have been performed with the aid of the multiconfiguration time-dependent Hartree method. Five coupled electronic potential energy surfaces and 10 (11 in the case of the orthoisomer) vibrational degrees of freedom have been included in these calculations. The nonadiabatic interactions lead to the restructuring of the photoelectron spectral envelopes. Ultrafast internal conversion processes within the electronic manifolds in question demonstrate the strength of the nonadiabatic coupling effects and complement the analogous findings for the electronic spectra. The internal conversion dynamics is characterized by a stepwise transfer of the electronic population to the lowest electronic state on a time scale of femtoseconds to picoseconds. A difference between the three isomers is found to be related to the weaker interaction between the sets of X-A and B-C-D states (with high-energy conical intersections) in the meta isomer, as compared to the other isomers. The implications of these findings for the qualitative understanding of the fluorescence dynamics of fluorinated benzene radical cations are discussed.

Classical and quantum studies of the photodissociation of a HX (X=Cl,F) molecule adsorbed on ice.

The photodissociation dynamics of a HX (X = Cl,F) molecule adsorbed on a hexagonal ice surface at T = 0 K is studied using time-dependent quantum wave packets and quasiclassical trajectories. The relevant potential energy surfaces are calculated using high-level ab initio methods. We present here two dimensional calculations for the dynamics of the hydrogen photofragment for both HCl and HF molecules. The purpose of this paper is to compare the photodissociation dynamics of the two molecules which are adsorbed on the ice surface with different equilibrium geometries. The total photodissociation cross section and the angular distribution are calculated. The comparison with classical trajectory calculations provides evidence for typical quantum effects and reveals rainbow structures.

Multimode Jahn-Teller and pseudo-Jahn-Teller interactions in the cyclopropane radical cation: complex vibronic spectra and nonradiative decay dynamics.

The complex vibronic spectra and the nonradiative decay dynamics of the cyclopropane radical cation (CP+) are simulated theoretically with the aid of a time-dependent wave packet propagation approach using the multireference time-dependent Hartree scheme. The theoretical results are compared with the experimental photoelectron spectrum of cyclopropane. The ground and first excited electronic states of CP+ are of X2E' and A2E'' type, respectively. Each of these degenerate electronic states undergoes Jahn-Teller (JT) splitting when the radical cation is distorted along the degenerate vibrational modes of e' symmetry. The JT split components of these two electronic states can also undergo pseudo-Jahn-Teller (PJT)-type crossings via the vibrational modes of e'', a1'' and a2'' symmetries. These lead to the possibility of multiple multidimensional conical intersections and highly nonadiabatic nuclear motions in these coupled manifolds of electronic states. In a previous publication [J. Phys. Chem. A 2004, 108, 2256], we investigated the JT interactions alone in the X2E' ground electronic manifold of CP+. In the present work, the JT interactions in the A2E'' electronic manifold are treated, and our previous work is extended by considering the coupling between the X2E' and A2E'' electronic states of CP+. The nuclear dynamics in this coupled manifold of two JT split doubly degenerate electronic states is simulated by considering fourteen active and most relevant vibrational degrees of freedom. The vibronic level spectra and the ultrafast nonradiative decay of the excited cationic states are examined and are related to the highly complex entanglement of electronic and nuclear degrees of freedom in this prototypical molecular system.

Multiconfiguration time-dependent Hartree method applied to molecular dissociation on surfaces: H2 + Pt(111).

Four-dimensional quantum dynamics calculations are performed on the dissociative chemisorption of H(2) on Pt(111) using the multiconfiguration time-dependent Hartree method. The aim of this work is to study the performance of the multiconfiguration time-dependent Hartree method for a gas-surface reaction by comparison with the standard time-dependent wave-packet propagation method. The initial-state resolved dissociation probability of H(2) is calculated within two four-dimensional models. The first four-dimensional model treats explicitly the rotational motion of the molecule and the H(2) dissociation is studied above two different sites of the Pt(111). For this model, only a potential-energy surface of general form was available. This potential was refitted to a sum of product form to allow efficient calculations with the multiconfiguration time-dependent Hartree method. The second model focuses on the description of the center-of-mass motion parallel to the surface, the rotational motion of the molecule being frozen. These four-dimensional quantum dynamics calculations yield important insights which can help with performing full six-dimensional calculations on H(2) dissociating on Pt(111). The multiconfiguration time-dependent Hartree method is shown to be particularly efficient for computing initial-state selective dissociation probabilities for the system studied, with a good accuracy and a reduced amount of memory and computational time when compared to the standard time-dependent wave-packet method.

On the unphysical impact of complex absorbing potentials on the Hamiltonian and its remedy.

The introduction of complex absorbing potentials as numerical tools to stabilize or increase the efficiency of calculations based on wave-packet propagation or on eigenvalue problems has the drawback of causing a modification of the Hamilton operator of the problem. In this work the consequences of such a modification are analyzed and the corrections required in order to properly describe the original physical process are derived. As an example, the decay of excited molecular states is considered: it is shown that the standard time-independent expression for the decay spectrum loses its validity when a complex absorbing potential is introduced in the nuclear Hamilton operator of the problem. To remedy the situation, a new, very stable formula is derived and tested on relevant model studies. Numerical examples are discussed.

Wave packet study of the UV photodissociation of the Ar2HBr complex.

The photodissociation dynamics of the Ar2HBr van der Waals molecule is studied using the multiconfiguration time-dependent Hartree method. Standard Jacobian coordinates are used to describe the molecule. Two four-dimensional calculations are carried out where the rotation of the Ar2 molecule and, in addition, either the vibration of Ar2-Br or that of Ar2 are frozen. The time-evolution of the probability density in the different modes and the calculation of the dissociative flux show that the dissociating hydrogen atom preferentially moves out of the plane defined by Ar2 and Br. A comprehensive study of the cage effect in the process is presented.

On the interatomic Coulombic decay in the Ne dimer.

The interatomic Coulombic decay (ICD) in the Ne dimer is discussed in view of the recent experimental results. The ICD electron spectrum and the kinetic energy release of the Ne+ fragments resulting after Coulomb explosion of Ne2 (2+) are computed and compared to the measured ones. A very good agreement is found, confirming the dynamics predicted for this decay mechanism. The effect of the temperature on the electron spectrum is briefly investigated.

Theoretical study of excitations in furan: spectra and molecular dynamics.

The excitation spectra and molecular dynamics of furan associated with its low-lying excited singlet states 1A2(3s), 1B2(V), 1A1(V'), and 1B1(3p) are investigated using an ab initio quantum-dynamical approach. The ab initio results of our previous work [J. Chem. Phys. 119, 737 (2003)] on the potential energy surfaces (PES) of these states indicate that they are vibronically coupled with each other and subject to conical intersections. This should give rise to complex nonadiabatic nuclear dynamics. In the present work the dynamical problem is treated using adequate vibronic coupling models accounting for up to four coupled PES and thirteen vibrational degrees of freedom. The calculations were performed using the multiconfiguration time-dependent Hartree method for wave-packet propagation. It is found that in the low-energy region the nuclear dynamics of furan is governed mainly by vibronic coupling of the 1A2(3s) and 1B2(V) states, involving also the 1A1(V') state. These interactions are responsible for the ultrafast internal conversion from the 1B2(V) state, characterized by a transfer of the electronic population to the 1A2(3s) state on a time scale of approximately 25 fs. The calculated photoabsorption spectrum of furan is in good qualitative agreement with experimental data. Some assignments of the measured spectrum are proposed.

Cumulative isomerization probability studied by various transition state wave packet methods including the MCTDH algorithm. Benchmark: HCN-->CNH isomerization.

The 3D cumulative isomerization probability N(E) for the transfer of a light particle between two atoms is computed by one time-independent and two time-dependent versions of the transition state wave packet (TSWP) method. The time-independent method is based on the direct expansion of the microcanonical projection operator on Chebyshev polynomials. In the time-dependent TSWP methods, the propagations are carried out by the split operator scheme and the multiconfiguration time-dependent Hartree (MCTDH) algorithm. This is the very first implementation of the TSWP method in the Heidelberg MCTDH package [G. W. Worth, M. H. Beck, A. Jackle, and H.-D. Meyer, The MCDTH package, Version 8.2 (2000); H.-D Meyer, Version 8.3 (2002). See http://www.pci.uni-heidelberg.de/tc/usr/mctdh/]. The benchmark is the HCN-->CNH isomerization for zero total angular momentum. Particular insights are given into the tunneling region. In larger systems, the time-dependent version of TSWP making use of the MCTDH algorithm will permit to treat more and more modes quantum mechanically, for very accurate results. Therefore, it was important to calibrate the implementation. Besides, we also assess the efficiency of a reduced dimensionality approach by comparing the new exact 3D calculations of N(E) for the HCN-->CNH isomerization with results obtained via 1D or 2D active subspaces. This suggests that, it should be possible to take directly benefit of the present 3D approaches, adapted for triatomic Jacobi coordinates to compute N(E) for H-transfer in larger systems, via three active coordinates. The prerequisite is then the simplification of the reduced 3D kinetic energy operator with rigid constraint to take the form corresponding to a pseudo triatomic system in Jacobi coordinates with effective masses. This last step is checked in the methoxy radical and malonaldehyde. Finally, different ways to obtain reliable eigenvectors of the flux operator associated with a dividing surface are revisited.