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.

In this article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational

Nitrobenzene is the simplest nitroaromatic compound and yet is characterized by a challenging and rich photophysics and photochemistry. In the present contribution, the main decay paths undertaken by the system after UV absorption from both the brightest (Lππ*) and the lowest (nπ*) singlet excited states have been characterized by means of CASPT2//CASSCF computations. The obtained results match with the main photophysical properties experimentally reported: the lack of fluorescence and phosphorescence emission is justified by the presence of accessible conical intersections and intersystem crossing regions between, respectively, the (nπ*) and (nπ*) states and the ground state, while the high triplet quantum yield is attributable to the strong coupling between the (nπ*) and (ππ*) states along the main decay path of the former. Two not previously reported singlet-triplet crossing regions, termed (T1/S0) and (T1/S0), have been here documented, from which the ground state can decay toward NO and phenoxy radical production and toward the formation of an epoxide ring structure, respectively. A possible mechanism leading to the photoisomerization of the nitro into the nitrite group, believed to be a key step in the photodegradation of nitrobenzene, has been proposed, based on the geometrical deformation recorded along the decay path leading from the (nπ*) state back to the original ground state through a conical intersection characterized by a significant shortening of the carbon-nitrogen bond.

The main intrinsic photochemical events in nucleobases can be described on theoretical grounds within the realm of non-adiabatic computational photochemistry. From a static standpoint, the photochemical reaction path approach (PRPA), through the computation of the respective minimum energy path (MEP), can be regarded as the most suitable strategy in order to explore the electronically excited isolated nucleobases. Unfortunately, the PRPA does not appear widely in the studies reported in the last decade. The main ultrafast decay observed experimentally for the gas-phase excited nucleobases is related to the computed barrierless MEPs from the bright excited state connecting the initial Franck-Condon region and a conical intersection involving the ground state. At the highest level of theory currently available (CASPT2//CASPT2), the lowest excited (1)(ππ*) hypersurface for cytosine has a shallow minimum along the MEP deactivation pathway. In any case, the internal conversion processes in all the natural nucleobases are attained by means of interstate crossings, a self-protection mechanism that prevents the occurrence of photoinduced damage of nucleobases by ultraviolet radiation. Many alternative and secondary paths have been proposed in the literature, which ultimately provide a rich and constructive interplay between experimentally and theoretically oriented research.

Pump-probe electronic spectroscopy using femtosecond laser pulses has evolved into a standard tool for tracking ultrafast excited state dynamics. Its two-dimensional (2D) counterpart is becoming an increasingly available and promising technique for resolving many of the limitations of pump-probe caused by spectral congestion. The ability to simulate pump-probe and 2D spectra from ab initio computations would allow one to link mechanistic observables like molecular motions and the making/breaking of chemical bonds to experimental observables like excited state lifetimes and quantum yields. From a theoretical standpoint, the characterization of the electronic transitions in the visible (Vis)/ultraviolet (UV), which are excited via the interaction of a molecular system with the incoming pump/probe pulses, translates into the determination of a computationally challenging number of excited states (going over 100) even for small/medium sized systems. A protocol is therefore required to evaluate the fluctuations of spectral properties like transition energies and dipole moments as a function of the computational parameters and to estimate the effect of these fluctuations on the transient spectral appearance. In the present contribution such a protocol is presented within the framework of complete and restricted active space self-consistent field theory and its second-order perturbation theory extensions. The electronic excited states of adenine have been carefully characterized through a previously presented computational recipe [Nenov et al., Comput. Theor. Chem. 1040-1041, 295-303 (2014]. A wise reduction of the level of theory has then been performed in order to obtain a computationally less demanding approach that is still able to reproduce the characteristic features of the reference data. Foreseeing the potentiality of 2D electronic spectroscopy to track polynucleotide ground and excited state dynamics, and in particular its expected ability to provide conformational dependent fingerprints in dimeric systems, the performances of the selected reduced level of calculations have been tested in the construction of 2D electronic spectra for the in vacuo adenine monomer and the unstacked adenine homodimer, thereby exciting the L b /L a transitions with the pump pulse pair and probing in the Vis to near ultraviolet spectral window. C 2015 AIP Publishing LLC. [http://dx

The photoinduced mechanism leading to the formation of the thymine-thymine (6-4) photolesion has been studied by using the CASPT2//CASSCF approach over a dinucleotide model in vacuo. Following light absorption, localization of the excitation on a single thymine leads to fast singlet-triplet crossing that populates the triplet (3)(nπ*) state of thymine. This state, displaying an elongated C(4)═O bond, triggers (6-4) dimer formation by reaction with the C(5)═C(6) double bond of the adjacent thymine, followed by a second intersystem crossing, which acts as a gate between the excited state of the reactant and the ground state of the photoproduct. The requirement of localized excitation on just one thymine, whose main decay channel (by radiationless repopulation of its ground state) is nonphotochemical, can rationalize the experimentally observed low quantum yield of formation for the thymine-thymine (6-4) adduct.

The photophysics and deactivation pathways of the noncanonical 5-azacytosine nucleobase were studied using the CASPT2//CASSCF protocol. One of the most significant differences with respect to the parent molecule cytosine is the presence of a dark (1)(nNπ*) excited state placed energetically below the bright excited state (1)(ππ*) at the Franck-Condon region. The main photoresponse of the system is a presumably efficient radiationless decay back to the original ground state, mediated by two accessible conical intersections involving a population transfer from the (1)(ππ*) and the (1)(nNπ*) states to the ground state. Therefore, a minor contribution of the triplet states in the photophysics of the system is expected, despite the presence of a deactivation path leading to the lowest (3)(ππ*) triplet state. The global scenario on the photophysics and photochemistry of the 5-azacytosine system gathered on theoretical grounds is consistent with the available experimental data, taking especially into account the low values of the singlet-triplet intersystem crossing and fluorescence quantum yields observed.

Indole is a chromophore present in many different molecules of biological interest, such as the essential amino acid tryptophan and the neurotransmitter serotonin. On the basis of CASPT2//CASSCF quantum chemical calculations, the photophysical properties of the system after UV irradiation have been studied through the exploration of the potential energy hypersurfaces of the singlet and triplet low-lying valence excited states. In contrast to previous studies, the present work has been carried out without imposing any restriction to the geometry of the molecule (C1 symmetry) and by performing minimum energy path calculations, which is the only instrument able to provide the lowest-energy evolution of the system. Relevant findings to the photophysics of bare indole have been obtained, which compete with the currently accepted mechanism for the energy decay in the molecule. The results show the presence of a conical intersection (CI) between the initially populated (1)(La ππ*) and the (1)(Lb ππ*) state, easily accessible through a barrierless pathway from the Franck-Condon region. At this CI region, part of the population is switched from the bright (1)(La ππ*) state to the (1)(Lb ππ*) state, and the system evolves toward a minimum structure from which the expected fluorescence takes place. The reported low values of the fluorescence quantum yield are explained by means of a new nonradiative mechanism specific for the (1)(Lb ππ*) state, in which the presence of an ethene-like CI between the (Lb ππ*) and ground states is the main feature.