The curvature of the conical intersection seam: an approximate second-order analysis

We present a method for analyzing the curvature (second derivatives) of the conical intersection hyperline at an optimized critical point. Our method uses the projected Hessians of the degenerate states after elimination of the two branching space coordinates, and is equivalent to a frequency calculation on a single Born-Oppenheimer potential-energy surface. Based on the projected Hessians, we develop an equation for the energy as a function of a set of curvilinear coordinates where the degeneracy is preserved to second order (i.e., the conical intersection hyperline). The curvature of the potential-energy surface in these coordinates is the curvature of the conical intersection hyperline itself, and thus determines whether one has a minimum or saddle point on the hyperline. The equation used to classify optimized conical intersection points depends in a simple way on the first- and second-order degeneracy splittings calculated at these points. As an example, for fulvene, we show that the two optimized conical intersection points of C2v symmetry are saddle points on the intersection hyperline. Accordingly, there are further intersection points of lower energy, and one of C2 symmetry - presented here for the first time - is found to be the global minimum in the intersection space

Photoisomerization for a model protonated Schiff base in solution: Sloped/peaked conical intersection perspective

The topographical character of conical intersections (CIs)-either sloped or peaked-has played a fundamental and important role in the discussion of the efficiency of CIs as photochemical "funnels." Here this perspective is employed in connection with a recent study of a model protonated Schiff base (PSB) cis to trans photoisomerization in solution [Malhado et al., J. Phys. Chem. A 115, 3720 (2011)]. In that study, the calculated reduced photochemical quantum yield for the successful production of trans product versus cis reactant in acetonitrile solvent compared to water was interpreted in terms of a dynamical solvent effect related to the dominance, for the acetonitrile case, of S-1 to S-0 nonadiabatic transitions prior to the reaching the seam of CIs. The solvent influence on the quantum yield is here re-examined in the sloped/peaked CI topographical perspective via conversion of the model's two PSB internal coordinates and a nonequilibrium solvent coordinate into an effective branching space description, which is then used to re-analyze the generalized Langevin equation/surface hopping results. The present study supports the original interpretation and enriches it in terms of topographical detail. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4754505]

Theoretical analysis of photoinduced H-atom elimination in thiophenol

The photoinduced hydrogen elimination reaction in thiophenol via the conical intersections of the dissociative (1)pi sigma* excited state with the bound (1)pi pi* excited state and the electronic ground state has been investigated with ab initio electronic-structure calculations and time-dependent quantum wave-packet calculations. A screening of the coupling constants of the symmetry-allowed coupling modes at the (1)pi pi*-(1)pi sigma* and (1)pi sigma*-S-0 conical intersection shows that the SH torsional mode is by far the most important coupling mode at both conical intersections. A model including three intersecting potential-energy surfaces (S-0, (1)pi pi*, (1)pi sigma*) and two nuclear degrees of freedom (SH stretch and SH torsion) has been constructed on the basis of ab initio complete-active-space self-consistent field and multireference second-order perturbation theory calculations. The nonadiabatic quantum wave-packet dynamics initiated by optical excitation of the (1)pi pi* and (1)pi sigma* states has been explored for this three-state two-coordinate model. The photodissociation dynamics is characterized in terms of snapshots of time-dependent wave packets, time-dependent electronic population probabilities, and the branching ratio of the (2)sigma/(2)pi electronic states of the thiophenoxyl radical. The dependence of the timescale of the photodissociation process and the branching ratio on the initial excitation of the SH stretching and SH torsional vibrations has been analyzed. It is shown that the node structure, which is imposed on the nuclear wave packets by the initial vibrational preparation as well as by the transitions through the conical intersections, has a profound effect on the photodissociation dynamics. The effect of additional weak coupling modes of CC twist (nu(16a)) and ring-distortion (nu(16b)) character has been investigated with three-dimensional and four-dimensional time-dependent wave-packet calculations, and has been found to be minor. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4709608]

A velocity map ion-imaging study on ketene photodissociation at 208 and 213 nm: Rotational dependence of product angular anisotropy

Photodissociation dynamics of ketene following excitation at 208.59 and 213.24 nm have been investigated using the velocity map ion-imaging method. Both the angular distribution and translational energy distribution of the CO products at different rotational and vibrational states have been obtained. No significant difference in the translational energy distributions for different CO rotational state products has been observed at both excitation wavelengths. The anisotropy parameter beta is, however, noticeably different for different CO rotational state products at both excitation wavelengths. For lower rotational states of the CO product, beta is smaller than zero, while beta is larger than zero for CO at higher rotational states. The observed rotational dependence of angular anisotropy is interpreted as the dynamical influence of a peculiar conical intersection between the B-1(1) excited state and (1)A(2) state along the C-S-I coordinate.

Automatic generation of active coordinates for quantum dynamics calculations: application to the dynamics of benzene photochemistry

A new practical method to generate a subspace of active coordinates for quantum dynamics calculations is presented. These reduced coordinates are obtained as the normal modes of an analytical quadratic representation of the energy difference between excited and ground states within the complete active space self-consistent field method. At the Franck-Condon point, the largest negative eigenvalues of this Hessian correspond to the photoactive modes: those that reduce the energy difference and lead to the conical intersection; eigenvalues close to 0 correspond to bath modes, while modes with large positive eigenvalues are photoinactive vibrations, which increase the energy difference. The efficacy of quantum dynamics run in the subspace of the photoactive modes is illustrated with the photochemistry of benzene, where theoretical simulations are designed to assist optimal control experiments

2-Aminopurine non-radiative decay and emission in aqueous solution: A theoretical study

The minimum energy path along the lowest-lying pi pi* excited state of 2-aminopurine was calculated to elucidate the mechanisms of radiationless decay and emission in water. The sequential Monte Carlo quantum mechanics approach with a multiconfigurational and perturbative description of the wave function was employed to compute the minimum, transition state, and conical intersection. It was found that the barrier in the potential energy surface to access the conical intersection funnel increases in aqueous environment, making the system prone to enlarge the emission yield. These results rationalize the observed enhancement of emission in 2-aminopurine upon increasing of the solvent polarity. (c) 2008 Elsevier B.V. All rights reserved.

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-zeta 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)(pi pi*) state can be populated. The first one takes place via a conical intersection involving the bright (1)(pi pi*) and the lowest (1)(n pi*) states, ((1)pi pi*/(1)n pi*)(CI), from which a low energy singlet-triplet crossing, ((1)n pi*/(3)pi pi*)(STC), connecting the (1)(n pi*) state to the lowest (3)(pi pi*) triplet state is accessible. The second mechanism arises via a singlet-triplet crossing, ((1)pi pi*/(3)n pi*)(STC), leading to a conical intersection in the triplet manifold, ((3)n pi*/(3)pi pi*)(CI), evolving to the lowest (3)(pi pi*) state. Further radiationless decay to the ground state is possible through a (gs/(3)pi pi*)(STC).

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.

On the Mechanisms of Triplet Excited State Population in 8-Azaadenine

The photophysics of 8-azaadenine (8-AA) has been studied with the CASPT2//CASSCF protocol and ANO-L double-zeta basis sets. Stationary equilibrium structures, surface crossings, minimum energy paths, and linear interpolations have been used to study possible mechanisms to populate the lowest triplet state, T-1 (3)(pi pi*), capable of sensitizing molecular oxygen. Our results show that two main mechanisms can occur after photoexcitation to the S-2 (1)(pi pi*) state. The first one is through the S-2/S-1 conical intersection (((1)pi pi*/(1)n pi*)(Cl)), leading to the S-1 ((1)n pi*) state minimum, (S-1 ((1)n pi*))(min), where a singlet-triplet crossing, ((1)n pi*/(3)pi pi*)(STC), is accessible. The second one starts with the ((1)pi pi*/(3)n pi*)(STC) at the (S-2((1)pi pi*))(min), from which the system can evolve to the (T-2 ((3)n pi*))(min), with subsequent population of the T-1 excited electronic state, due to the ((3)n pi*/(3)pi pi*)(Cl) conical intersection.

Computational chemistry studies of UV induced processes in human skin

This thesis presents and uses the techniques of computational chemistry to explore two different processes induced in human skin by ultraviolet light. The first is the transformation of urocanic acid into a immunosuppressing agent, and the other is the enzymatic action of the 8-oxoguanine glycosylase enzyme. The photochemistry of urocanic acid is investigated by time-dependent density functional theory. Vertical absorption spectra of the molecule in different forms and environments is assigned and candidate states for the photochemistry at different wavelengths are identified. Molecular dynamics simulations of urocanic acid in gas phase and aqueous solution reveals considerable flexibility under experimental conditions, particularly for for the cis isomer where competition between intra- and inter-molecular interactions increases flexibility. A model to explain the observed gas phase photochemistry of urocanic acid is developed and it is shown that a reinterpretation in terms of a mixture between isomers significantly enhances the agreement between theory and experiment , and resolves several peculiarities in the spectrum. A model for the photochemistry in the aqueous phase of urocanic acid is then developed, in which two excited states governs the efficiency of photoisomerization. The point of entrance into a conical intersection seam is shown to explain the wavelength dependence of photoisomerization quantum yield. Finally some mechanistic aspects of the DNA repair enzyme 8-oxoguanine glycosylase is investigated with density functional theory. It is found that the critical amino acid of the active site can provide catalytic power in several different manners, and that a recent proposal involving a SN1 type of mechanism seems the most efficient one.

Excited-State Structure, Vibrations, and Nonradiative Relaxation of Jet-Cooled 5-Fluorocytosine

The S0 → S1 vibronic spectrum and S1 state nonradiative relaxation of jet-cooled keto-amino 5-fluorocytosine (5FCyt) are investigated by two-color resonant two-photon ionization spectroscopy at 0.3 and 0.05 cm–1 resolution. The 000 rotational band contour is polarized in-plane, implying that the electronic transition is 1ππ*. The electronic transition dipole moment orientation and the changes of rotational constants agree closely with the SCS-CC2 calculated values for the 1ππ* (S1) transition of 5FCyt. The spectral region from 0 to 300 cm–1 is dominated by overtone and combination bands of the out-of-plane ν1′ (boat), ν2′ (butterfly), and ν3′ (HN–C6H twist) vibrations, implying that the pyrimidinone frame is distorted out-of-plane by the 1ππ* excitation, in agreement with SCS-CC2 calculations. The number of vibronic bands rises strongly around +350 cm–1; this is attributed to the 1ππ* state barrier to planarity that corresponds to the central maximum of the double-minimum out-of-plane vibrational potentials along the ν1′, ν2′, and ν3′ coordinates, which gives rise to a high density of vibronic excitations. At +1200 cm–1, rapid nonradiative relaxation (knr ≥ 1012 s–1) sets in, which we interpret as the height of the 1ππ* state barrier in front of the lowest S1/S0 conical intersection. This barrier in 5FCyt is 3 times higher than that in cytosine. The lifetimes of the ν′ = 0, 2ν1′, 2ν2′, 2ν1′ + 2ν2′, 4ν2′, and 2ν1′ + 4ν2′ levels are determined from Lorentzian widths fitted to the rotational band contours and are τ ≥ 75 ps for ν′ = 0, decreasing to τ ≥ 55 ps at the 2ν1′ + 4ν2′ level at +234 cm–1. These gas-phase lifetimes are twice those of S1 state cytosine and 10–100 times those of the other canonical nucleobases in the gas phase. On the other hand, the 5FCyt gas-phase lifetime is close to the 73 ps lifetime in room-temperature solvents. This lack of dependence on temperature and on the surrounding medium implies that the 5FCyt nonradiative relaxation from its S1 (1ππ*) state is essentially controlled by the same ∼1200 cm–1 barrier and conical intersection both in the gas phase and in solution.

Intersystem crossing rates of S1 state keto-amino cytosine at low excess energy

The amino-keto tautomer of supersonic jet-cooled cytosine undergoes intersystem crossing (ISC) from the v = 0 and low-lying vibronic levels of its S1(¹ππ*) state. We investigate these ISC rates experimentally and theoretically as a function of S1 state vibrational excess energy Eexc. The S1 vibronic levels are pumped with a ~5 ns UV laser, the S1 and triplet state ion signals are separated by prompt or delayed ionization with a second UV laser pulse. After correcting the raw ISC yields for the relative S1 and T1ionization cross sections, we obtain energy dependent ISC quantum yields Q corr ISC =1% –5%. These are combined with previously measured vibronic state-specific decay rates, giving ISC rates kISC = 0.4–1.5 ⋅ 10⁹ s⁻¹, the corresponding S1⇝S0internal conversion (IC) rates are 30–100 times larger. Theoretical ISC rates are computed using SCS-CC2 methods, which predict rapid ISC from the S1; v = 0 state with kISC = 3 ⋅ 10⁹ s⁻¹ to the T1(³ππ*) triplet state. The surprisingly high rate of this El Sayed-forbidden transition is caused by a substantial admixture of ¹nOπ* character into the S1(¹ππ*) wave function at its non-planar minimum geometry. The combination of experiment and theory implies that (1) below Eexc = 550 cm⁻¹ in the S1 state, S1⇝S0internal conversion dominates the nonradiative decay with kIC ≥ 2 ⋅ 10¹⁰ s⁻¹, (2) the calculated S1⇝T1 (¹ππ*⇝³ππ*) ISC rate is in good agreement with experiment, (3) being El-Sayed forbidden, the S1⇝T1 ISC is moderately fast (kISC = 3 ⋅ 10⁹ s⁻¹), and not ultrafast, as claimed by other calculations, and (4) at Eexc ~ 550 cm⁻¹ the IC rate increases by ~50 times, probably by accessing the lowest conical intersection (the C5-twist CI) and thereby effectively switching off the ISC decay channels.

The relaxation dynamics of the excited electronic states of retinal in bacteriorhodopsin by two-pump-probe femtosecond studies

We present the results of two-pump and probe femtosecond experiments designed to follow the relaxation dynamics of the lowest excited state (S1) populated by different modes. In the first mode, a direct (S0 → S1) radiative excitation of the ground state is used. In the second mode, an indirect excitation is used where the S1 state is populated by the use of two femtosecond laser pulses with different colors and delay times between them. The first pulse excites the S0 → S1 transition whereas the second pulse excites the S1 → Sn transition. The nonradiative relaxation from the Sn state populates the lowest excited state. Our results suggest that the S1 state relaxes faster when populated nonradiatively from the Sn state than when pumped directly by the S0 → S1 excitation. Additionally, the Sn → S1 nonradiative relaxation time is found to change by varying the delay time between the two pump pulses. The observed dependence of the lowest excited state population as well as its dependence on the delay between the two pump pulses are found to fit a kinetic model in which the Sn state populates a different surface (called S′1) than the one being directly excited (S1). The possible involvement of the Ag type states, the J intermediate, and the conical intersection leading to the S0 or to the isomerization product (K intermediate) are discussed in the framework of the proposed model.

Spin boson models for quantum decoherence of electronic excitations of biomolecules and quantum dots in a solvent

We give a theoretical treatment of the interaction of electronic excitations (excitions) in biomolecules and quantum dots with the surrounding polar solvent. Significant quantum decoherence occurs due to the interaction of the electric dipole moment of the solute with the fluctuating electric dipole moments of the individual molecules in the solvent. We introduce spin boson models which could be used to describe the effects. of decoherence on the quantum dynamics of biomolecules which undergo light-induced conformational change and on biomolecules or quantum dots which are coupled by Forster resonant energy transfer.

Ground and Electronic Excited States from Pairing Matrix Fluctuation and Particle-Particle Random Phase Approximation

The accurate description of ground and electronic excited states is an important and challenging topic in quantum chemistry. The pairing matrix fluctuation, as a counterpart of the density fluctuation, is applied to this topic. From the pairing matrix fluctuation, the exact electron correlation energy as well as two electron addition/removal energies can be extracted. Therefore, both ground state and excited states energies can be obtained and they are in principle exact with a complete knowledge of the pairing matrix fluctuation. In practice, considering the exact pairing matrix fluctuation is unknown, we adopt its simple approximation --- the particle-particle random phase approximation (pp-RPA) --- for ground and excited states calculations. The algorithms for accelerating the pp-RPA calculation, including spin separation, spin adaptation, as well as an iterative Davidson method, are developed. For ground states correlation descriptions, the results obtained from pp-RPA are usually comparable to and can be more accurate than those from traditional particle-hole random phase approximation (ph-RPA). For excited states, the pp-RPA is able to describe double, Rydberg, and charge transfer excitations, which are challenging for conventional time-dependent density functional theory (TDDFT). Although the pp-RPA intrinsically cannot describe those excitations excited from the orbitals below the highest occupied molecular orbital (HOMO), its performances on those single excitations that can be captured are comparable to TDDFT. The pp-RPA for excitation calculation is further applied to challenging diradical problems and is used to unveil the nature of the ground and electronic excited states of higher acenes. The pp-RPA and the corresponding Tamm-Dancoff approximation (pp-TDA) are also applied to conical intersections, an important concept in nonadiabatic dynamics. Their good description of the double-cone feature of conical intersections is in sharp contrast to the failure of TDDFT. All in all, the pairing matrix fluctuation opens up new channel of thinking for quantum chemistry, and the pp-RPA is a promising method in describing ground and electronic excited states.