Combined theoretical and experimental study of the photophysics of asulam.

The photophysics of the neutral molecular form of the herbicide asulam has been described in a joint experimental and theoretical, at the CASPT2 level, study. The unique π → π* aromatic electronic transition (f, ca. 0.5) shows a weak red-shift as the polarity of the solvent is increased, whereas the fluorescence band undergoes larger red-shifts. Solvatochromic data point to higher dipole moment in the excited state than in the ground state (µ(g) < µ(e)). The observed increase in pKa in the excited state (pKa* - pKa, ca. 3) is consistent with the results of the Kamlet-Abboud-Taft and Catalán et al. multiparametric approaches. Fluorescence quantum yield varies with the solvent, higher in water (ϕ(f) = 0.16) and lower in methanol and 1-propanol (approx. 0.02). Room temperature fluorescence lifetime in aqueous solution is (1.0 ± 0.2) ns, whereas the phosphorescence lifetime in glassy EtOH at 77 K and the corresponding quantum yield are (1.1 ± 0.1) s and 0.36, respectively. The lack of mirror image symmetry between modified absorption and fluorescence spectra reflects different nuclear configurations in the absorbing and emitting states. The low value measured for the fluorescence quantum yield is justified by an efficient nonradiative decay channel, related with the presence of an easily accessible conical intersection between the initially populated singlet bright (1)(L(a) ππ*) state and the ground state (gs/ππ*)(CI). Along the main decay path of the (1)(L(a) ππ*) state the system undergoes an internal conversion process that switches part of the population from the bright (1)(L(a) ππ*) to the dark (1)(L(b) ππ*) state, which is responsible for the fluorescence. Additionally, singlet-triplet crossing regions have been found, a fact that can explain the phosphorescent emission detected. An intersystem crossing region between the phosphorescent state (3)(L(a) ππ*) and the ground state has been characterized, which contributes to the nonradiative deactivation of the excitation energy.

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.

Proton/Hydrogen Transfer Mechanisms in the Guanine-Cytosine Base Pair: Photostability and Tautomerism.

Proton/hydrogen-transfer processes have been broadly studied in the past 50 years to explain the photostability and the spontaneous tautomerism in the DNA base pairs. In the present study, the CASSCF/CASPT2 methodology is used to map the two-dimensional potential energy surfaces along the stretched NH reaction coordinates of the guanine-cytosine (GC) base pair. Concerted and stepwise pathways are explored initially in vacuo, and three mechanisms are studied: the stepwise double proton transfer, the stepwise double hydrogen transfer, and the concerted double proton transfer. The results are consistent with previous findings related to the photostability of the GC base pair, and a new contribution to tautomerism is provided. The C-based imino-oxo and imino-enol GC tautomers, which can be generated during the UV irradiation of the Watson-Crick base pair, have analogous radiationless energy-decay channels to those of the canonical base pair. In addition, the C-based imino-enol GC tautomer is thermally less stable. A study of the GC base pair is carried out subsequently taking into account the DNA surroundings in the biological environment. The most important stationary points are computed using the quantum mechanics/molecular mechanics (QM/MM) approach, suggesting a similar scenario for the proton/hydrogen-transfer phenomena in vacuo and in DNA. Finally, the static model is complemented by ab initio dynamic simulations, which show that vibrations at the hydrogen bonds can indeed originate hydrogen-transfer processes in the GC base pair. The relevance of the present findings for the rationalization of the preservation of the genetic code and mutagenesis is discussed.

On the deactivation mechanisms of adenine-thymine base pair.

In this contribution, the multiconfigurational second-order perturbation theory method based on a complete active space reference wave function (CASSCF/CASPT2) is applied to study all possible single and double proton/hydrogen transfers between the nucleobases in the adenine-thymine (AT) base pair, analyzing the role of excited states with different nature [localized (LE) and charge transfer (CT)], and considering concerted as well as step-wise mechanisms. According to the findings, once the lowest excited states, localized in adenine, are populated during UV irradiation of the Watson-Crick base pair, the proton transfer in the N-O bridge does not require high energy in order to populate a CT state. The latter state will immediately relax toward a crossing with the ground state, which will funnel the system to either the canonical structure or the imino-enol tautomer. The base pair is also capable of repairing itself easily since the imino-enol species is unstable to thermal conversion.

Distinct photophysics of the isomers of B18H22 explained.

The photophysics of the two isomers of octadecaborane(22), anti- and syn-B(18)H(22), have been studied by UV-vis spectroscopic techniques and theoretical computational methods. In air-saturated hexane, anti-B(18)H(22) shows fluorescence with a high quantum yield, Φ(F) = 0.97, and singlet oxygen O(2)((1)Δ(g)) production (Φ(Δ) ∼ 0.008). Conversely, isomer syn-B(18)H(22) shows no measurable fluorescence, instead displaying much faster, picosecond nonradiative decay of excited singlet states. Computed potential energy hypersurfaces (PEHs) for both isomers rationalize these data, pointing to a deep S(1) minimum for anti-B(18)H(22) and a conical intersection (CI) between its S(0) and S(1) states that lies 0.51 eV higher in energy. Such an energy barrier to nonradiative relaxation is not present in the PEH of syn-B(18)H(22), and the system therefore has sufficient initial energy on excitation to reach the (S(0)/S(1)) CI and to then decay to the ground state without fluorescence. The computational analysis of the geometries at stationary points along the PEH of both isomers shows that the determining factor for the dissimilar photophysics of anti- and syn-B(18)H(22) may be due to the significant differences in the geometrical rearrangements at their respective conical intersections. Thus, the syn isomer shows one very large, B-B elongation of 1.2 Å from 1.8 Å in the ground state to 3.0 Å at the CI, whereas the anti isomer shows smaller elongations (below 1 Å) in several B-B connectivities at its (S(0)/S(1))(CI). The absorbed energy in S(1) for the anti-B(18)H(22) is therefore redistributed vibrationally into several regions of the molecule rather than almost completely into a single vibrational mode as in the case for the syn isomer. The consequent prolonged S(1) lifetime for the anti isomer allows for relaxation via fluorescence.

Progress and challenges in the calculation of electronic excited states.

A detailed understanding of the properties of electronic excited states and the reaction mechanisms that molecules undergo after light irradiation is a fundamental ingredient for following light-driven natural processes and for designing novel photonic materials. The aim of this review is to present an overview of the ab initio quantum chemical and time-dependent density functional theory methods that can be used to model spectroscopy and photochemistry in molecular systems. The applicability and limitations of the different methods as well as the main frontiers are discussed. To illustrate the progress achieved by excited-state chemistry in the recent years as well as the main challenges facing computational chemistry, three main applications that reflect the authors' experience are addressed: the UV/Vis spectroscopy of organic molecules, the assignment of absorption and emission bands of organometallic complexes, and finally, the obtainment of non-adiabatic photoinduced pathways mediated by conical intersections. In the latter case, special emphasis is put on the photochemistry of DNA. These applications show that the description of electronically excited states is a rewarding but challenging area of research.

On the N1-H and N3-H Bond Dissociation in Uracil by Low Energy Electrons: A CASSCF/CASPT2 Study.

The dissociative electron-attachment (DEA) phenomena at the N1-H and N3-H bonds observed experimentally at low energies (<3 eV) in uracil are studied with the CASSCF/CASPT2 methodology. Two valence-bound π(-) and two dissociative σ(-) states of the uracil anionic species, together with the ground state of the neutral molecule, are proven to contribute to the shapes appearing in the experimental DEA cross sections. Conical intersections (CI) between the π(-) and σ(-) are established as the structures which activate the DEA processes. The N1-H and N3-H DEA mechanisms in uracil are described, and experimental observations are interpreted on the basis of two factors: (1) the relative energy of the (U-H)(-) + H fragments obtained after DEA with respect to the ground-state equilibrium structure (S0) of the neutral molecule (threshold for DEA) and (2) the relative energy of the CIs also with respect to S0 (band maxima). The π1(-) state is found to be mainly responsible for the N1-H bond breaking, whereas the π2(-) state is proved to be involved in the cleavage of the N3-H bond.

A distance-dependent parameterization of the extended Hubbard model for conjugated and aromatic hydrocarbons derived from stretched ethene.

The Hubbard model, which is widely used in physics but is mostly unfamiliar to chemists, provides an attractive yet simple model for chemistry beyond the self consistent field molecular orbital approximation. The Hubbard model adds an effective electron-electron repulsion when two electrons occupy the same atomic orbital to the familiar Hückel Hamiltonian. Thus it breaks the degeneracy between excited singlet and triplet states and allows an explicit treatment of electron correlation. We show how to evaluate the parameters of the model from high-level ab initio calculations on two-atom fragments and then to transfer the parameters to large molecules and polymers where accurate ab initio calculations are difficult or impossible. The recently developed MS-RASPT2 method is used to generate accurate potential energy curves for ethene as a function of carbon-carbon bond length, which are used to parameterize the model for conjugated hydrocarbons. Test applications to several conjugated/aromatic molecules show that even though the model is very simple, it is capable of reasonably accurate predictions for bond lengths, and predicts molecular excitation energies in reasonable agreement with those from the MS-RASPT2 method.

Intramolecular charge transfer and dual fluorescence of 4-(dimethylamino)benzonitrile: ultrafast branching followed by a two-fold decay mechanism.

In this contribution we present new experimental and theoretical results for the intramolecular charge transfer (ICT) reaction underlying the dual fluorescence of 4-(dimethylamino)benzonitrile (DMABN), which indicate that the fully twisted ICT (TICT) state is responsible for the time-resolved transient absorption spectrum while a distinct partially twisted ICT (pTICT) structure is suggested for the fluorescent ICT state.

Ultrafast ring-opening/closing and deactivation channels for a model spiropyran-merocyanine system.

The photochemistry of a model merocyanine-spiropyran system was analyzed theoretically at the MS-CASPT2//SA-CASSCF(14,12) level. Several excited singlet states were studied in both the closed spiropyran and open merocyanine forms, and the paths to the different S(1)/S(0) conical intersections found were analyzed. After absorption of UV light from the spiropyran form, there are two possible ultrafast routes to efficient conversion to the ground state; one involves the rupture of the C(spiro)-O bond leading to the open form and the other involves the lengthening of the C(spiro)-N bond with no photoreaction. From the merocyanine side the excited state can reach a very broad S(1)/S(0) conical intersection region that leads the system to the closed form after rotation of the central methine bond. Alternatively, rotation of the other methine bonds connects the system through different S(1)/S(0) conical intersections to several merocyanine isomers. The present work provides a theoretical framework for the recent experimental results (Buback , J. J. Am. Chem. Soc. 2010, 132, 1610-1619) and sheds light on the complex photochemistry of these kinds of compounds.

On the relaxation mechanisms of 6-azauracil.

The nonadiabatic photochemistry of 6-azauracil has been studied by means of the CASPT2//CASSCF protocol and double-ζ plus polarization ANO basis sets. Minimum energy states, transition states, minimum energy paths, and surface intersections have been computed in order to obtain an accurate description of several potential energy hypersurfaces. It is concluded that, after absorption of ultraviolet radiation (248 nm), two main relaxation mechanisms may occur, via which the lowest (3)(ππ*) state can be populated. The first one takes place via a conical intersection involving the bright (1)(ππ*) and the lowest (1)(nπ*) states, ((1)ππ*/(1)nπ*)(CI), from which a low-energy singlet-triplet crossing, ((1)nπ*/(3)ππ*)(STC), connecting the (1)(nπ*) state to the lowest (3)(ππ*) triplet state is accessible. The second mechanism arises via a singlet-triplet crossing, ((1)ππ*/(3)nπ*)(STC), leading to a conical intersection in the triplet manifold, ((3)nπ*/(3)ππ*)(CI), evolving to the lowest (3)(ππ*) state. Further radiationless decay to the ground state is possible through a (gs/(3)ππ*)(STC).

Multiconfigurational Second-Order Perturbation Theory Restricted Active Space (RASPT2) Method for Electronic Excited States: A Benchmark Study.

The recently developed second-order perturbation theory restricted active space (RASPT2) method has been benchmarked versus the well-established complete active space (CASPT2) approach. Vertical excitation energies for valence and Rydberg excited states of different groups of organic (polyenes, acenes, heterocycles, azabenzenes, nucleobases, and free base porphin) and inorganic (nickel atom and copper tetrachloride dianion) molecules have been computed at the RASPT2 and multistate (MS) RASPT2 levels using different reference spaces and compared with CASPT2, CCSD, and experimental data in order to set the accuracy of the approach, which extends the applicability of multiconfigurational perturbation theory to much larger and complex systems than previously. Relevant aspects in multiconfigurational excited state quantum chemistry such as the valence-Rydberg mixing problem in organic molecules or the double d-shell effect for first-row transition metals have also been addressed.

Theoretical investigations of the IR spectroscopy of Ni(C(2)S(2)H(2))(2). A case study of the P_VMWCI(2) algorithm including anharmonic effects.

The near infrared (NIR) spectra of bis(ethylene-1,2-dithiolato)nickel, Ni(C(2)S(2)H(2))(2) are fully interpreted here by applying a method developed for efficient automatic computation of both the infrared wave numbers and the intensities. The employed procedure uses parallel variational multiple window configuration interaction wave functions, the so-named P_VMWCI(2) algorithm, which incorporates both the mechanical and the electric anharmonic effects. It is shown that inclusion of anharmonicities is crucial for correctly assigning the fundamental, combination, and overtone vibrational frequencies in the infrared spectrum of the target system, for which conflicting assignments are found in literature.

Asymmetry and non-adiabaticity in fragmentation of disulfide bonds upon electron capture.

Although it has been generally assumed that electron attachment to disulfide derivatives leads to a systematic and significant activation of the S-S bond, we show, by using [CH(3)SSX] (X = CH(3), NH(2), OH, F) derivatives as model compounds, that this is the case only when the X substituents have low electronegativity. Through the use of MP2, QCI and CASPT2 molecular orbital (MO) methods, we elucidate, for the first time, the mechanisms that lead to unimolecular fragmentation of disulfide derivatives after electron attachment. Our theoretical scrutiny indicates that these mechanisms are more intricate than assumed in previous studies. The most stable products, from a thermodynamic viewpoint, correspond to the release of neutral molecules; CH(4), NH(3), H(2)O, and HF. However, the barriers to reach these products depend strongly on the electronegativity of the X substituents. Only for very electronegative substituents, such as OH or F, the loss of H(2)O or HF is the most favorable process, and likely the only one observed. This is possible because of two concomitant factors, 1) the extra electron for [CH(3)SSX](-) (X = OH, F) occupies a sigma*(S-X) MO, which favors the cleavage of the S-X bond, and 2) the activation barriers associated with the hydrogen transfer process to produce H(2)O and HF are rather low. Only when the substituents are less electronegative (X = H, CH(3), NH(2)) the extra electron is located in a sigma*(S-S) orbital and the cleavage of the disulfide bridge becomes the most favorable process. The intimate mechanism associated with the S-S bond dissociation process also depends strongly on the nature of the substituent. For X = H or CH(3) the process is strictly adiabatic, while for X = NH(2) it proceeds through a conical intersection (CI) associated with the charge reorganization necessary to obtain, from a molecular anion with the extra electron delocalized in a sigma*(S-S) antibonding orbital, two fragments with the proper charge localization.

Thioxanthone: on the shape of the first absorption band.

The equilibrium ground state geometry of thioxanthone (TX) has been investigated and its effect on the vertical excitation energies and photophysical behaviour has been explained. In line with this purpose, the first absorption band of TX has been simulated and analysed in detail. The calculations show that TX is planar, C(2v) symmetric in its ground state. The energy of the low-lying excited states seems to be rather insensitive along the butterfly motion coordinate. The shoulder in the first absorption band (at around 3.43 eV) is shown to be caused by vibrational progression of various in-plane modes and does not justify the hypothesis that two photophysically distinct conformers of TX exist. The calculated vertical excitation spectrum (in vacuum) has been compared with the experimental absorption bands in non-polar solvents.

Electron capture activation of the disulfide bond. The role of the asymmetry and electronegativity.

The effects of electron capture on the structure of XSSX' disulfide derivatives in which the substituents attached to the sulfur atoms have different electronegativites have been investigated at different levels of theory, namely DFT, MP2, QCISD and CASSCF/CASPT2. Although it has been generally assumed that electron attachment to disulfide derivatives leads to a systematic and significant activation of the S-S bond, our results show that this is the case only when the substituents X or X' have low electronegativity. Otherwise, the S-S bond in the anion remains practically unperturbed and only the S-X bond is largely activated or even broken, because the extra electron occupies the sigma*(S-X) rather than the sigma*(S-S) antibonding orbital. Our results also show that S-S activation yields a system with a unique anion, whereas when the S-X activation is significant, two stable anionic species, stretched and bent, are formed.

Ultrafast decay of the excited singlet states of thioxanthone by internal conversion and intersystem crossing.

The experimental ultrafast photophysics of thioxanthone in several aprotic organic solvents at room temperature is presented, measured using femtosecond transient absorption together with high-level ab initio CASPT2 calculations of the singlet- and triplet-state manifolds in the gas phase, including computed state minima and conical intersections, transition energies, oscillator strengths, and spin-orbit coupling terms. The initially populated singlet pi pi* state is shown to decay through internal conversion and intersystem crossing processes via intermediate n pi* singlet and triplet states, respectively. Two easily accessible conical intersections explain the favorable internal conversion rates and low fluorescence quantum yields in nonpolar media. The presence of a singlet-triplet crossing near the singlet pi pi* minimum and the large spin-orbit coupling terms also rationalize the high intersystem crossing rates. A phenomenological kinetic scheme is proposed that accounts for the decrease in internal conversion and intersystem crossing (i.e. the very large experimental crescendo of the fluorescence quantum yield) with the increase of solvent polarity.

Calibration of Cholesky Auxiliary Basis Sets for Multiconfigurational Perturbation Theory Calculations of Excitation Energies.

The accuracy of auxiliary basis sets derived from Cholesky decomposition of two-electron integrals is assessed for excitation energies calculated at the state-average complete active space self-consistent field (CASSCF) and multiconfigurational second order perturbation theory (CASPT2) levels of theory using segmented as well as generally contracted atomic orbital basis sets. Based on 196 valence excitations in 26 organic molecules and 72 Rydberg excitations in 3 organic molecules, the results show that Cholesky auxiliary basis sets can be used without compromising the accuracy of the multiconfigurational methods. Specifically, with a decomposition threshold of 10(-4) au, the mean error due to the Cholesky auxiliary basis set is 0.001 eV, or smaller, decreasing with increasing atomic orbital basis set quality.

Singlet-Triplet States Interaction Regions in DNA/RNA Nucleobase Hypersurfaces.

The present study provides new insight into the intrinsic mechanisms for the population of the triplet manifold in DNA nucleobases by determining, at the multiconfigurational CASSCF/CASPT2 level, the singlet-triplet states crossing regions and the main decay paths for their lowest singlet and triplet states after near-UV irradiation. The studied singlet-triplet interacting regions are accessible along the minimum energy path of the initially populated singlet bright (1)ππ* state. In particular, all five natural DNA/RNA nucleobases have, at the end of the main minimum energy path and near a conical intersection of the ground and (1)ππ* states, a low-energy, easily accessible, singlet-triplet crossing region directly connecting the lowest singlet and triplet ππ* excited states. Adenine, thymine, and uracil display additional higher-energy crossing regions related to the presence of low-lying singlet and a triplet nπ* state. These funnels are absent in guanine and cytosine, which have the bright (1)ππ* state lower in energy and less accessible nπ* states. Knowledge of the location and accessibility of these regions, in which the singlet-triplet interaction is related to large spin-orbit coupling elements, may help to understand experimental evidence such as the wavelength dependence measured for the triplet formation quantum yield in nucleobases and the prevalence of adenine and thymine components in the phosphorescence spectra of DNA.

MOLCAS 7: the next generation.

Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented in this report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two-electron integrals and in the generation of so-called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self-consistent field, density functional theory, 2nd order perturbation theory, complete-active space self-consistent field multiconfigurational reference 2nd order perturbation theory, and coupled-cluster methods. The report further elaborates on the implementation of a restricted-active space self-consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas-Kroll-Hess transformation for one-component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so-called picture-change-free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies.