Surface hopping dynamics methods using the coupled cluster to approximated second order (CC2), the algebraic diagrammatic construction scheme to second order (ADC(2)), and the time-dependent density functional theory (TDDFT) were developed and implemented into the program system Newton-X. These procedures are especially well-suited to simulate nonadiabatic processes involving various excited states of the same multiplicity and the dynamics in the first excited state toward an energetic minimum or up to the region where a crossing with the ground state is found. 9H-adenine in the gas phase was selected as the test case. The results showed that dynamics with ADC(2) is very stable, whereas CC2 dynamics fails within 100 fs, because of numerical instabilities present in the case of quasi-degenerate excited states. ADC(2) dynamics correctly predicts the ultrafast character of the deactivation process. It predicts that C2-puckered conical intersections should be the preferential pathway for internal conversion for low-energy excitation. C6-puckered conical intersection also contributes appreciably to internal conversion, becoming as important as C2-puckered for high-energy excitations. In any case, H-elimination plays only a minor role. TDDFT based on a long-range corrected functional fails to predict the ultrafast deactivation. In the comparison with several other methods previously used for dynamics simulations of adenine, ADC(2) has the best performance, providing the most consistent results so far.

No abstract

Ab initio nonadiabatic dynamics simulations of cis-to-trans isomerization of azobenzene upon S(1) (n-π*) excitation are carried out employing the fewest-switches surface hopping method. Azobenzene photoisomerization occurs purely as a rotational motion of the central CNNC moiety. Two nonequivalent rotational pathways corresponding to clockwise or counterclockwise rotation are available. The course of the rotational motion is strongly dependent on the initial conditions. The internal conversion occurs via an S(0)/S(1) crossing seam located near the midpoint of both of these rotational pathways. Based on statistical analysis, it is shown that the occurrence of one or other pathway can be completely controlled by selecting adequate initial conditions.

Newton-X is an open-source computational platform to perform nonadiabatic molecular dynamics based on surface hopping and spectrum simulations using the nuclear ensemble approach. Both are among the most common methodologies in computational chemistry for photophysical and photochemical investigations. This paper describes the main features of these methods and how they are implemented in Newton-X. It emphasizes the newest developments, including zero-point-energy leakage correction, dynamics on complex-valued potential energy surfaces, dynamics induced by incoherent light, dynamics based on machine-learning potentials, exciton dynamics of multiple chromophores, and supervised and unsupervised machine learning techniques. Newton-X is interfaced with several third-party quantum-chemistry programs, spanning a broad spectrum of electronic structure methods.

We present surface hopping dynamics on potential energy surfaces resulting from the spin-orbit splitting, i.e., surfaces corresponding to the eigenstates of the total electronic Hamiltonian including the spin-orbit coupling. In this approach, difficulties arise because of random phases of degenerate eigenvectors and possibility of crossings of the resulting mixed states. Our implementation solves these problems and allows propagation of the coefficients both in the representation of the spin free Hamiltonian and directly in the "diagonal representation" of the mixed states. We also provide a detailed discussion of the state crossing and point out several peculiarities that were not mentioned in the previous literature. We also incorporate the effect of the environment via the quantum mechanics/molecular mechanics approach. As a test case, we apply our methodology to deactivation of thiophene and selenophene in the gas phase, ethanol solution, and bulk liquid phase. First, 100 trajectories without spin-orbit coupling have been calculated for thiophene starting both in S and S states. A subset of 32 initial conditions starting in the S state was then used for gas phase simulations with spin-orbit coupling utilizing the 3-step integrator of SHARC, our implementation of the 3-step propagator in Newton-X and two new "one-step" approaches. Subsequently, we carried out simulations in ethanol solution and bulk liquid phase for both thiophene and selenophene. For both molecules, the deactivation of the S state proceeds via the ring opening pathway. The total population of triplet states reaches around 15% and 40% after 80 fs for thiophene and selenophene, respectively. However, it only begins growing after the ring opening is initiated; hence, the triplet states do not directly contribute to the deactivation mechanism. For thiophene, the resulting deactivation lifetime of the S state was 68 fs in the gas phase, 76 fs in ethanol solution, and 78 fs in the liquid phase, in a good agreement with the experimental value of 80 fs (liquid phase). For selenophene, the obtained S lifetime was 60 fs in the gas phase and 62 fs for both ethanol solution and liquid phase. The higher rate of intersystem crossing to the triplet states in selenophene is likely the reason for the lower fluorescence observed in selenium containing polymer compounds.