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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

The role of adenine excimers in the photophysics of oligonucleotides.

Energies and structures of different arrangements of the stacked adenine homodimer have been computed at the ab initio CASPT2 level of theory in isolation and in an aqueous environment. Adenine dimers are shown to form excimer singlet states with different degrees of stacking and interaction. A model for a 2-fold decay dynamics of adenine oligomers can be supported in which, after initial excitation in the middle UV range, unstacked or slightly stacked pairs of nucleobases will relax by an ultrafast internal conversion to the ground state, localizing the excitation in the monomer and through the corresponding conical intersection with the ground state. On the other hand, long-lifetime intrastrand stacked excimer singlet states will be formed in different conformations, including neutral and charge transfer dimers, which originate the red-shifted emission observed in the oligonucleotide chains and that will evolve toward the same monomer decay channel after surmounting an energy barrier. By computing the transient absorption spectra for the different structures considered and their relative stability in vacuo and in water, it is concluded that in the adenine homodimers the maximum-overlap face-to-face orientations are the most stable excimer conformations observed in experiment.

Electrostatic control of the photoisomerization efficiency and optical properties in visual pigments: on the role of counterion quenching.

Hybrid QM(CASPT2//CASSCF/6-31G*)/MM(Amber) computations have been used to map the photoisomerization path of the retinal chromophore in Rhodopsin and explore the reasons behind the photoactivity efficiency and spectral control in the visual pigments. It is shown that while the electrostatic environment plays a central role in properly tuning the optical properties of the chromophore, it is also critical in biasing the ultrafast photochemical event: it controls the slope of the photoisomerization channel as well as the accessibility of the S(1)/S(0) crossing space triggering the ultrafast decay. The roles of the E113 counterion, the E181 residue, and the other amino acids of the protein pocket are explicitly analyzed: it appears that counterion quenching by the protein environment plays a key role in setting up the chromophore's optical properties and its photochemical efficiency. A unified scenario is presented that discloses the relationship between spectroscopic and mechanistic properties in rhodopsins and allows us to draw a solid mechanism for spectral tuning in color vision pigments: a tunable counterion shielding appears as the elective mechanism for L<-->M spectral modulation, while a retinal conformational control must dictate S absorption. Finally, it is suggested that this model may contribute to shed new light into mutations-related vision deficiencies that opens innovative perspectives for experimental biomolecular investigations in this field.

Determination of the electron-detachment energies of 2'-deoxyguanosine 5'-monophosphate anion: influence of the conformation.

The vertical electron-detachment energies (VDEs) of the singly charged 2'-deoxyguanosine 5'-monophosphate anion (dGMP-) are determined by using the multiconfigurational second-order perturbation CASPT2 method at the MP2 ground-state equilibrium geometry of relevant conformers. The origin of the unique low-energy band in the gas phase photoelectron spectrum of dGMP-, with maximum at around 5.05 eV, is unambiguously assigned to electron detachment from the highest occupied molecular orbital of pi-character belonging to guanine fragment of a syn conformation. The presence of a short H-bond linking the 2-amino and phosphate groups, the guanine moiety acting as proton donor, is precisely responsible for the pronounced decrease of the computed VDE with respect to that obtained in other conformations. As a whole, the present research supports the nucleobase as the site with the lowest ionization potential in negatively charged (deprotonated) nucleotides at the most stable conformations as well as for B-DNA-like type arrangements, in agreement with experimental evidence.

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.

Do fluorescence and transient absorption probe the same intramolecular charge transfer state of 4-(dimethylamino)benzonitrile?

We present here the results of time-resolved absorption and emission experiments for 4-(dimethylamino)benzonitrile in solution, which suggest that the fluorescent intramolecular charge transfer (ICT) state may differ from the twisted ICT (TICT) state observed in transient absorption.

Linear and nonlinear optical properties of a series of Ni-dithiolene derivatives.

Some linear and nonlinear optical (NLO) properties of Ni(SCH)(4) and several of its derivatives have been computed by employing a series of basis sets and a hierarchy of methods (e.g., HF, DFT, coupled cluster, and multiconfigurational techniques). The electronic structure of Ni(SCH)(4) has been also analyzed by using CASSCF/CASPT2, ab initio valence bond, and DFT methods. In particular we discuss how the diradicaloid character (DC) of Ni(SCH)(4) significantly affects its NLO properties. The quasidegeneracy of the two lowest-energy singlet states 1 (1)A(g) and 1 (1)B(1u), the clear DC nature of the former, and the very large number of low-lying states enhance the NLO properties values. These particular features are used to interpret the NLO properties of Ni(SCH)(4). The DC of the considered derivatives has been estimated and correlated with the NLO properties. CASVB computations have shown that the structures with Ni(II) are the dominant ones, while those with Ni(0) and Ni(IV) have negligible weight. The weights of the four diradical structures were discussed in connection with the weight of the structures, where all the electrons are paired. Comparative discussion of the properties of Ni(SCH)(4) with those of tetrathia fulvalene demonstrates the very large effect of Ni on the properties of the Ni-dithiolene derivatives. A similar remarkable effect on the NLO properties is produced by one or two methyl or C(3)S groups. The considered Ni-dithiolene derivatives have exceptionally large NLO properties. This feature in connection with their other physical properties makes them ideal candidates for photonic applications.

A three-state model for the photophysics of guanine.

The nonadiabatic photochemistry of the guanine molecule (2-amino-6-oxopurine) and some of its tautomers has been studied by means of the high-level theoretical ab initio quantum chemistry methods CASSCF and CASPT2. Accurate computations, based by the first time on minimum energy reaction paths, states minima, transition states, reaction barriers, and conical intersections on the potential energy hypersurfaces of the molecules lead to interpret the photochemistry of guanine and derivatives within a three-state model. As in the other purine DNA nucleobase, adenine, the ultrafast subpicosecond fluorescence decay measured in guanine is attributed to the barrierless character of the path leading from the initially populated 1(pi pi* L(a)) spectroscopic state of the molecule toward the low-lying methanamine-like conical intersection (gs/pi pi* L(a))CI. On the contrary, other tautomers are shown to have a reaction energy barrier along the main relaxation profile. A second, slower decay is attributed to a path involving switches toward two other states, 1(pi pi* L(b)) and, in particular, 1(n(O) pi*), ultimately leading to conical intersections with the ground state. A common framework for the ultrafast relaxation of the natural nucleobases is obtained in which the predominant role of a pi pi*-type state is confirmed.

Molecular basis of DNA photodimerization: intrinsic production of cyclobutane cytosine dimers.

Based on CASPT2 results, the present contribution establishes for the first time that cytosine photodimer formation (C< >C) is mediated along the triplet and singlet manifold by a singlet-triplet crossing, (T1/S0)X, and by a conical intersection, (S1/S0)CI, respectively. The former can be accessed in a barrierless way from a great variety of photochemical avenues and exhibits a covalent single bond between the ethene C6-C6' carbon atoms of each monomer. The efficiency of the stepwise triplet mechanism, however, would be modulated by the effectiveness of the intersystem crossing mechanism. The results provide the grounds for the understanding of the potential photogenotoxicity of endogenous and exogenous compounds via triplet-triplet sensitization, with a lower bound for cytosine oligonucleotides predicted to be 2.70 eV, and give support to the traditional view of the primary role of triplet excited states in the photochemistry of DNA, a well-known source of photoproducts in solution under triplet photosensitization conditions. The function played by singlet excimers (excited dimers) to explain both the red-shifted fluorescence and photoreaction is highlighted. A rationale on the pronounced wavelength dependence of the observed fluorescence is offered. Geometrical arrangements at the time of light irradiation close to, but energetically above, (S1/S0)CI are suggested as reactive orientations that become prone to produce C< >C directly, with no energy barrier. Because of the outstanding intrinsic ability of cytosine to form stable relaxed excimers, the system located near the bound relaxed excimer has to accumulate enough vibrational energy to surmount a small barrier of 0.2 eV to reach (S1/S0)CI, making the overall process to proceed at a slower relative rate as compared to other compounds such as thymine, which is not susceptible of forming so stable excimers.