Angelo Giussani

MOLCAS 8: New Capabilities for Multiconfigurational Quantum Chemical Calculations across the Periodic Table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas–Kroll–Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC‐PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large‐scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method 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, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Excitation of Nucleobases from a Computational Perspective I: Reaction Paths

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.

Photoinduced Formation Mechanism of the Thymine–Thymine (6–4) Adduct

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.

Essential on the Photophysics and Photochemistry of the Indole Chromophore by Using a Totally Unconstrained Theoretical Approach.

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.

Modeling the high-energy electronic state manifold of adenine: Calibration for nonlinear electronic spectroscopy

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 Lb/La transitions with the pump pulse pair and probing in the Vis to near ultraviolet spectral window.

Ultraviolet vision: photophysical properties of the unprotonated retinyl Schiff base in the Siberian hamster cone pigment

Abstract The Siberian hamster ultraviolet (SHUV) visual pigment has an unprotonated Schiff-base (SB) retinyl chromophore in the dark state, which becomes protonated after photoexcitation during the early stages of the photobleaching cycle. While the photochemical relaxation processes of the SHUV remain poorly understood, they are expected to show significant differences when compared to those of the protonated SB (PSB) chromophore in visual rhodopsin. Here, we report a study of the photophysical properties of the SHUV unprotonated SB (SHUV-USB), based on multiconfigurational and multireference perturbative methods within a hybrid quantum mechanics/molecular mechanics scheme. Comparisons of multireference and time-dependent density functional theory results indicate that both methodologies predict an ionic excited state (S1), similar to the PSB of rhodopsin, although its minimum has even bond-lengths in the central region of the retinyl polyene chain. The analysis of excited-state manifolds at the Franck–Condon region and S1 minimum configuration indicates that the skeletal relaxation initiated in the S1 surface is likely to involve S1/S2 surface crossing. These results provide valuable insights for future studies of the SHUV-USB photoisomerization mechanism.

Multiple Decay Mechanisms and 2D-UV Spectroscopic Fingerprints of Singlet Excited Solvated Adenine-Uracil Monophosphate.

Abstract The decay channels of singlet excited adenine uracil monophosphate (ApU) in water are studied with CASPT2//CASSCF:MM potential energy calculations and simulation of the 2D‐UV spectroscopic fingerprints with the aim of elucidating the role of the different electronic states of the stacked conformer in the excited state dynamics. The adenine 1La state can decay without a barrier to a conical intersection with the ground state. In contrast, the adenine 1Lb and uracil S(U) states have minima that are separated from the intersections by sizeable barriers. Depending on the backbone conformation, the CT state can undergo inter‐base hydrogen transfer and decay to the ground state through a conical intersection, or it can yield a long‐lived minimum stabilized by a hydrogen bond between the two ribose rings. This suggests that the 1Lb, S(U) and CT states of the stacked conformer may all contribute to the experimental lifetimes of 18 and 240 ps. We have also simulated the time evolution of the 2D‐UV spectra and provide the specific fingerprint of each species in a recommended probe window between 25 000 and 38 000 cm−1 in which decongested, clearly distinguishable spectra can be obtained. This is expected to allow the mechanistic scenarios to be discerned in the near future with the help of the corresponding experiments. Our results reveal the complexity of the photophysics of the relatively small ApU system, and the potential of 2D‐UV spectroscopy to disentangle the photophysics of multichromophoric systems.

Spectroscopic fingerprints of DNA/RNA pyrimidine nucleobases in third-order nonlinear electronic spectra

Accurate ab initio modeling of spectroscopic signals in nonlinear electronic spectra, such as bidimensional (2D) spectra, requires the computation of the electronic transitions induced by the incoming pump/probe pulses, resulting in a challenging calculation of many electronic excited states. A protocol is thus required to evaluate the variations of spectral properties, like transition energies and dipole moments, with the computational level, and to estimate the sensitivity of the spectra to these variations. Such a protocol is presented here within the framework of complete and restricted active space self-consistent field (CASSCF/RASSCF) theory and its second-order perturbation theory extensions (CASPT2/RASPT2). The electronic excited-state manifolds of pyrimidine nucleobases (thymine, uracil, and cytosine) are carefully characterized in vacuo employing high-level RAS(0,0|10,8|2,12)//SS-RASPT2 calculations. The results provide a reference data set that can be used for optimizing computational efforts and costs, as required for studying computationally more demanding multichromophoric systems (e.g., di- and oligonucleotides). The spectroscopic signatures of the 2D electronic spectrum of a perfectly stacked uracil–cytosine dimer model are characterized, and experimental setups are proposed that can resolve non-covalent interchromophoric interactions in canonical pyrimidine nucleobase-stacked dimers.

Probing deactivation pathways of DNA nucleobases by two-dimensional electronic spectroscopy: first principles simulations

The SOS//QM/MM [Rivalta et al., Int. J. Quant. Chem., 2014, 114, 85] method consists of an arsenal of computational tools allowing accurate simulation of one-dimensional (1D) and bi-dimensional (2D) electronic spectra of monomeric and dimeric systems with unprecedented details and accuracy. Prominent features like doubly excited local and excimer states, accessible in multi-photon processes, as well as charge-transfer states arise naturally through the fully quantum-mechanical description of the aggregates. In this contribution the SOS//QM/MM approach is extended to simulate time-resolved 2D spectra that can be used to characterize ultrafast excited state relaxation dynamics with atomistic details. We demonstrate how critical structures on the excited state potential energy surface, obtained through state-of-the-art quantum chemical computations, can be used as snapshots of the excited state relaxation dynamics to generate spectral fingerprints for different de-excitation channels. The approach is based on high-level multi-configurational wavefunction methods combined with non-linear response theory and incorporates the effects of the solvent/environment through hybrid quantum mechanics/molecular mechanics (QM/MM) techniques. Specifically, the protocol makes use of the second-order Perturbation Theory (CASPT2) on top of Complete Active Space Self Consistent Field (CASSCF) strategy to compute the high-lying excited states that can be accessed in different 2D experimental setups. As an example, the photophysics of the stacked adenine-adenine dimer in a double-stranded DNA is modeled through 2D near-ultraviolet (NUV) spectroscopy.

Photophysics of 1-Aminonaphthalene: A Theoretical and Time-Resolved Experimental Study

The photophysics of 1-aminonaphthalene (1-napthylamine, AMN) has been investigated on the basis of a constructive experimental-theoretical interplay derived from time-resolved measurements and high-level quantum-chemical ab initio CASPT2//CASSCF calculations. Transient ionization signals at femtosecond resolution were collected for AMN cold isolated molecules following excitation from the vibrationless ground level to a number of vibrational states (within the pump resolution) in the lowest accessible excited state and further multiphoton ionization probing at 500, 800, and 1300 nm. Theory predicts two pipi* states, (1)L(b) and (1)L(a), as the lowest singlet electronic excitations, with adiabatic transitions from S(0) at 3.50 and 3.69 eV, respectively. Since the associated oscillator strength for the lowest transition is exceedingly small, the (1)L(b) state is not expected to become populated significantly and the (1)L(a) state appears as the main protagonist of the AMN photophysics. Though calculations foresee a surface crossing between (1)L(a) and the lower (1)L(b) states, no dynamical signature of it is observed in the time-dependent measurements. In the relaxation of (1)L(a), the radiant emission competes with the intersystem crossing and internal conversion channels. The rates of these mechanisms have been determined at different excitation energies. The internal conversion is mediated by a (1)L(a)/S(0) conical intersection located 0.7 eV above the (1)L(a) minimum. The relaxation of a higher-lying singlet excited state, observed above 40 000 cm(-1) (4.96 eV) and calculated at 5.18 eV, has been also explored.