RESUMO
WavePacket is an open-source program package for numerical simulations in quantum dynamics. Building on the previous Part I (Schmidt and Lorenz, Comput. Phys. Commun. 2017, 213, 223] and Part II (Schmidt and Hartmann, Comput. Phys. Commun. 2018, 228, 229] which dealt with quantum dynamics of closed and open systems, respectively, the present Part III adds fully classical and mixed quantum-classical propagation techniques to WavePacket. There classical phase-space densities are sampled by trajectories which follow (diabatic or adiabatic) potential energy surfaces. In the vicinity of (genuine or avoided) intersections of those surfaces, trajectories may switch between them. To model these transitions, two classes of stochastic algorithms have been implemented: (1) Tully's fewest switches surface hopping and (2) Landau-Zener-based single switch surface hopping. The latter one offers the advantage of being based on adiabatic energy gaps only, thus not requiring nonadiabatic coupling information any more. The present work describes the MATLAB version of WavePacket 6.1.0, which is essentially an object-oriented rewrite of previous versions, allowing to perform fully classical, quantum-classical and quantum-mechanical simulations on an equal footing, that is, for the same physical system described by the same WavePacket input. The software package is hosted and further developed at the Sourceforge platform, where also extensive Wiki-documentation as well as numerous worked-out demonstration examples with animated graphics are available. © 2019 Wiley Periodicals, Inc.
RESUMO
The excited state dynamics of an unnatural base pair (UBP) d5SICS/dNaM were investigated by accurate ab-initio calculations. Time-dependent density functional and high-level multireference calculations (MS-CASPT2) were performed to elucidate the excitation of this UBP and its excited state relaxation mechanism. After excitation to the bright state S2 (ππ*), it decays to the S1 state and then undergoes efficient intersystem crossing to the triplet manifold. The presence of sulfur atom in d5SICS leads to strong spin-orbit coupling (SOC) and a small energy gap that facilitates intersystem crossing from S1 (ns π*) to T2 (ππ*) followed by internal conversion to T1 state. Similarly in dNaM, the deactivation pathway follows analogous trends. CASPT2 calculations suggest that the S1 (ππ*) state is a dark state below the accessible S2 (ππ*) bright state. During the ultrafast deactivation, it exhibits bond length inversion. From S1 state, significant SOC leads the population transfer to T3 due to a smaller energy gap. Henceforth, fast internal conversion occurs from T3 to T2 followed by T1 . From time-dependent trajectory surface hopping dynamics, it is found that excited state relaxation occurs on a sub-picosecond timescale in d5SICS and dNaM. Our findings strongly suggest that there is enough energy available in triplet state of UBP to generate reactive oxygen species and induce phototoxicity with respect to cellular DNA.