Effect of Initial Conditions Sampling on Surface Hopping Simulations in the Ultrashort and Picosecond Time Range. Azomethane Photodissociation as a Case Study

: We tested the effect of different ways of sampling the initial conditions in surface hopping simulations, with a focus on the initial energy distributions and on the treatment of the zero point energy (ZPE). As a test case, we chose the gas phase photodynamics of azomethane, which features different processes occurring in overlapping time scales: geometry relaxation in the excited state, internal conversion, photoisomerization, and fast and slow dissociation. The simulations, based on a semiempirical method, had a sufficiently long duration (10 ps) to encompass all of the above processes. We tested several variants of methods based on the quantum mechanical (QM) distributions of the nuclear coordinates q and momenta p, which yield, at least on the average over a large sampling set, the correct QM energy, namely the ZPE when starting from the ground vibrational state. We compared the QM samplings with the classical Boltzmann (CB) distribution obtained by a thermostated trajectory, whereby thermal effects are taken into account, but the ZPE is utterly ignored. We found that most QM and CB approaches yield similar results as to short time dynamics and decay lifetimes, whereas the rate of the ground state dissociation reaction CH3NNCH3 → CH3NN + CH3 is sharply affected by the sampling method. With QM samplings a large fraction of trajectories dissociate promply (<1 ps) after decay to the ground state and with rates of the order of 10-1 ps-1 after the first ps. Instead, the CB samplings yield a much smaller fraction of prompt dissociations and much lower rates at long times. We provided evidence that the ZPE "leaks" from high frequency modes to the reactive ones (N-C bond elongations), therefore unphysically increasing the dissociation rates with QM samplings. We show that an effective way to take into account the ZPE and to avoid the "leaking" problem is to add the ZPE to the potential energy surfaces as a function of the most relevant internal coordinates. Then, Boltzmann sampling can be done as usual, so this approach is suitable also for condensed state dynamics. In the tests we present here, the ZPE correction method yields dissociation rates intermediate between QM and uncorrected Boltzmann samplings.