Interplay of radiative and nonradiative transitions in surface hopping with radiation-molecule interactions.

We implemented a method for the treatment of field induced transitions in trajectory surface hopping simulations, in the framework of the local diabatization scheme, especially suited for on-the-fly dynamics. The method is applied to a simple one-dimensional model with an avoided crossing and compared with quantum wavepacket dynamics. The results show the importance of introducing a proper decoherence correction to surface hopping, in order to obtain meaningful results. Also the energy conservation policy of standard surface hopping must be revised: in fact, the quantum wavepacket energetics is well reproduced if energy absorption/emission is allowed for in the hops determined by radiation-molecule coupling. To our knowledge, this is the first time the issues of decoherence and energy conservation have been analyzed in depth to devise a mixed quantum-classical method for dynamics with molecule-field interactions.

Mixed quantum-classical dynamics in the adiabatic representation to simulate molecules driven by strong laser pulses.

The dynamics of molecules under strong laser pulses is characterized by large Stark effects that modify and reshape the electronic potentials, known as laser-induced potentials (LIPs). If the time scale of the interaction is slow enough that the nuclear positions can adapt to these externally driven changes, the dynamics proceeds by adiabatic following, where the nuclei gain very little kinetic energy during the process. In this regime we show that the molecular dynamics can be simulated quite accurately by a semiclassical surface-hopping scheme formulated in the adiabatic representation. The nuclear motion is then influenced by the gradients of the laser-modified potentials, and nonadiabatic couplings are seen as transitions between the LIPs. As an example, we simulate the process of adiabatic passage by light induced potentials in Na(2) using the surface-hopping technique both in the diabatic representation based on molecular potentials and in the adiabatic representation based on LIPs, showing how the choice of the representation is crucial in reproducing the results obtained by exact quantum dynamical calculations.