to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added and two surface hopping algorithms are included to enable nonadiabatic calculations. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree–Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

The study concerns the relaxation of electronic excited states of the DNA nucleoside deoxycytidine (dCyd) and its methylated analogue 5-methyldeoxycytidine (5mdCyd), known to be involved in the formation of UV-induced lesions of the genetic code. Due to the existence of four closely lying and potentially coupled excited states, the deactivation pathways in these systems are particularly complex and have not been assessed so far. Here, we provide a complete mechanistic picture of the excited state relaxation of dCyd/5mdCyd in three solvents-water, acetonitrile, and tetrahydrofuran-by combining femtosecond fluorescence experiments, addressing the effect of solvent proticity on the relaxation dynamics of dCyd and 5mdCyd for the first time, and two complementary quantum mechanical approaches (CASPT2/MM and PCM/TD-CAM-B3LYP). The lowest energy ππ* state is responsible for the sub-picosecond lifetime observed for dCyd in all the solvents. In addition, computed excited state absorption and transient IR spectra allow one, for the first time, to assign the tens of picoseconds time constant, reported previously, to a dark state (nπ*) involving the carbonyl lone pair. A second low-lying dark state, involving the nitrogen lone pair (nπ*), does significantly participate in the excited state dynamics. The 267 nm excitation of dCyd leads to a non-negligible population of the second bright ππ* state, which affects the dynamics, acting mainly as a "doorway" state for the nπ* state. The solvent plays a key role governing the interplay between the different excited states; unexpectedly, water favors population of the dark states. In the case of 5mdCyd, an energy barrier present on the main nonradiative decay route explains the 6-fold lengthening of the excited state lifetime compared to that of dCyd, observed for all the examined solvents. Moreover, C5-methylation destabilizes both nπ* and nπ* dark states, thus preventing them from being populated.