In this paper, we present benchmark results for dissipative dynamics of a harmonic oscillator coupled to an anharmonic bath of Morse oscillators. The microscopic Hamiltonian has been chosen so that the anharmonicity can be adjusted as a free parameter, and its effect can be isolated. This leads to a temperature dependent spectral density of the bath, which is studied for ohmic and lorentzian cases. Also, we compare numerically exact multiconfiguration time-dependent Hartree results with approximate solutions using continuous configuration time-dependent self-consistent field and local coherent state approximation.

The nonadiabatic photodissociation dynamics of (HI)2 is simulated by applying a wave packet approach which starts from the I*-HI complex (where I* denotes the I(2P1/2) excited electronic state) produced after the photodissociation of the first HI moiety within (HI)2. In the model, two excited electronic potential surfaces corresponding to I*-HI(A 1Pi1) and I-HI(A 1Pi1), which interact through spin-rotation coupling, are considered. The simulations show that upon photodissociation of HI within I*-HI, the dissociating H fragment undergoes intracluster collisions with the I* atom. Some of these collisional events induce an electronically nonadiabatic transition which causes the deactivation of I* to the I ground electronic state. The probability of such nonadiabatic process is found to be 0.37%. Most of the photodissociation process takes place in the upper excited electronic surface [that of the I*-HI(A 1Pi1) complex], where H dissociation is found to be mainly direct or involving weak H/I* intracluster collisions. These weak collisions with high collisional angular momentum, and therefore high collisional impact parameters associated, are responsible for most of the probability of nonadiabatic transitions found. The type of H/I* collisions leading to nonadiabatic transitions appears to be closely related to the nature of the spin-rotation coupling between the two excited electronic states involved.

A potential energy surface for the ground electronic state of the Ar-HI van der Waals complex is calculated at the coupled-cluster with single and double excitations and a noniterative perturbation treatment of triple excitations [CCSD(T)] level of theory. Calculations are performed using for the iodine atom a correlation consistent triple-zeta valence basis set in conjunction with large-core Stuttgart-Dresden-Bonn relativistic pseudopotential, whereas specific augmented correlation consistent basis sets are employed for the H and Ar atoms supplemented with an additional set of bond functions. In agreement with previous studies, the equilibrium structure is found to be linear Ar-I-H, with a well depth of 205.38 cm(-1). Another two secondary minima are also predicted at a linear and bent Ar-H-I configurations with well depths of 153.57 and 151.57 cm(-1), respectively. The parametrized CCSD(T) potential is used to calculate rovibrational bound states of Ar-HI/Ar-DI complexes, and the vibrationally averaged structures of the different isomers are determined. Spectroscopic constants are also computed from the CCSD(T) surface and their comparison with available experimental data demonstrates the quality of the present surface in the corresponding configuration regions.

Typical (sub)system-bath quantum dynamical problems are often investigated by means of (approximate) reduced equations of motion. Wavepacket approaches to the dynamics of the whole system have gained momentum in recent years and there is hope that properly designed approximations to the wavefunction will allow one to correctly describe the subsystem evolution. The continuous-configuration time-dependent self-consistent field (CC-TDSCF) and local coherent-state approximation (LCSA) methods, for instance, use a simple Hartree product of bath single-particle-functions for each discrete variable representation (DVR) state introduced in the Hilbert space of the subsystem. Here we focus on the above two methods and replace the DVR states with the eigenstates of the subsystem Hamiltonian, i.e., we adopt an energy-local representation for the subsystem. We find that stable and semiquantitative results are obtained for a number of dissipative problems, at the same (small) computational cost of the original methods. Furthermore, we find that both methods give very similar results, thus suggesting that coherent-states are well suited to describe (local) bath states. As a whole, present results highlight the importance of the system basis-set in the selected-multiconfiguration expansion of the wavefunction. They suggest that accurate and yet computationally cheap methods may be simply obtained from CC-TDSCF/LCSA by letting the subsystem states be variationally optimized.

The sequential photodissociation dynamics of (HI)2 is studied by means of a nonadiabatic wave packet treatment starting from the I*-HI complex. The model reproduces the main experimental findings for photolysis with 266 nm radiation. The results confirm that some of the H atoms dissociated from the I*-HI complex deactivate the I* atom through a HI* intracluster collision which induces an I*-->I electronically nonadiabatic transition. As a consequence, these H fragments become very fast by acquiring nearly all the I* excitation energy, equivalent to the I*I spin-orbit splitting. A most interesting result is the high production of bound I2 fragments in highly excited rovibrational states in the photolysis, indicating that the H dissociation is mainly direct.