Pathway approach to ultrafast photochemistry: potential surfaces, conical intersections and isomerizations of small polyenes

Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.

From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions

This review discusses how ultrafast organic photochemical reactions are controlled by conical intersections, highlighting that decay to the ground-state at multiple points of the intersection space results in their multi-mode character.

Ultrafast Excited State Dynamics of Spatially Confined Organic Molecules

This Feature Article highlights the role of spatial confinement in controlling the fundamental behavior of molecules. Select examples illustrate the value of using space as a tool to control and understand excited-state dynamics through a combination of ultrafast spectroscopy and conventional steady-state methods. Molecules of interest were confined within a closed molecular capsule, derived from a cavitand known as octa acid (OA), whose internal void space is sufficient to accommodate molecules as long as tetracene and as wide as pyrene. The free space, i.e., the space that is left following the occupation of the guest within the host, is shown to play a significant role in altering the behavior of guest molecules in the excited state. The results reported here suggest that in addition to weak interactions that are commonly emphasized in supramolecular chemistry, the extent of empty space (i.e., the remaining void space within the capsule) is important in controlling the excited-state behavior of confined molecules on ultrafast time scales. For example, the role of free space in controlling the excited-state dynamics of guest molecules is highlighted by probing the cis-trans isomerization of stilbenes and azobenzenes within the OA capsule. Isomerization of both types of molecule are slowed when they are confined within a small space, with encapsulated azobenzenes taking a different reaction pathway compared to that in solution upon excitation to S₂. In addition to steric constraints, confinement of reactive molecules in a small space helps to override the need for diffusion to bring the reactants together, thus enabling the measurement of processes that occur faster than the time scale for diffusion. The advantages of reducing free space and confining reactive molecules are illustrated by recording unprecedented excimer emission from anthracene and by measuring ultrafast electron transfer rates across the organic molecular wall. By monitoring the translational motion of anthracene pairs in a restricted space, it has been possible to document the pathway undertaken by excited anthracene from inception to the formation of the excimer on the excited-state surface. Similarly, ultrafast electron transfer experiments pursued here have established that the process is not hindered by a molecular wall. Apparently, the electron can cross the OA capsule wall provided the donor and acceptor are in close proximity. Measurements on the ultrafast time scale provide crucial insights for each of the examples presented here, emphasizing the value of both "space" and "time" in controlling and understanding the dynamics of excited molecules.

Internal Conversion between Bright (11Bu+) and Dark (21Ag-) States in s-trans-Butadiene and s-trans-Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1¹B_u⁺ and 2¹A_g^-, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N_o}, and DEA-EOMCC(4p2h){N_u} or DEA-EOMCC(3p1h,4p2h){N_u}, consistently predict that the 1¹B_u⁺/2¹A_g^- crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1¹B_u⁺ and 2¹A_g^- potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT2, do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR(4,4), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

A systematic model study quantifying how conical intersection topography modulates photochemical reactions

Despite their important role in photochemistry and expected presence in most polyatomic molecules, conical intersections have been thoroughly characterized in a comparatively small number of systems. Conical intersections can confer molecular photoreactivity or photostability, often with remarkable efficacy, due to their unique structure: at a conical intersection, the adiabatic potential energy surfaces of two or more electronic states are degenerate, enabling ultrafast decay from an excited state without radiative emission, known as nonadiabatic transfer. Furthermore, the precise conical intersection topography determines fundamental properties of photochemical processes, including excited-state decay rate, efficacy, and molecular products that are formed. However, these relationships have yet to be defined comprehensively. In this article, we use an adaptable computational model to investigate a variety of conical intersection topographies, simulate resulting nonadiabatic dynamics, and calculate key photochemical observables. We varied the vibrational mode frequencies to modify conical intersection topography systematically in four primary classes of conical intersections and quantified the resulting rate, total yield, and product yield of nonadiabatic decay. The results reveal that higher vibrational mode frequencies reduce nonadiabatic transfer, but increase the transfer rate and resulting photoproduct formation. These trends can inform progress toward experimental control of photochemical reactions or tuning of molecules' photochemical properties based on conical intersections and their topography.

Selective bond formation triggered by short optical pulses: quantum dynamics of a four-center ring closure

Making bonds with attopulses: quantum dynamics of the ring closure of norbornadiene to quadricyclane.

Ultrafast ring closing of a diarylethene-based photoswitchable nucleoside

Deoxyuridine nucleosides embodied into diarylethenes form an especial class of photoswitchable compounds that are designed to stack and pair with DNA bases. The molecular geometry can be switched between "open" and "closed" isomers by a pericyclic reaction that affects the stability of the surrounding double helix. This potentially enables light-induced control of DNA hybridization at microscopic resolution. Despite its importance for the optimization of DNA photoswitches, the ultrafast photoisomerization mechanism of these diarylethenes is still not well understood. In this work, femtosecond transient absorption spectroscopy is applied to study the ring closing reaction upon UV excitation with 45 fs pulses. Excited-state absorption decays rapidly and gives rise to the UV-Vis difference spectrum of the "closed" form within ≈15 ps. Time constants of 0.09, 0.49 and 6.6 ps characterize the multimodal dynamics, where a swift recurrence in the signal anisotropy indicates transient population of the intermediate 21A-like state.

Ultrafast trans → cis Photoisomerization Dynamics of Alkyl-Substituted Stilbenes in a Supramolecular Capsule

Ultrafast spectroscopy reveals the effects of confinement on the excited-state photoisomerization dynamics for a series of alkyl-substituted trans-stilbenes encapsulated in the hydrophobic cavity of an aqueous supramolecular organic host-guest complex. Compared with the solvated compounds, encapsulated trans-stilbenes have broader excited-state absorption spectra, excited-state lifetimes that are 3-4 times longer, and photoisomerization quantum yields that are 1.7-6.5 times lower in the restricted environment. The organic capsule disrupts the equilibrium structure and restricts torsional rotation around the central C═C double bond in the excited state, which is an important motion for the relaxation of trans-stilbene from S₁ to S₀. The location and identity of alkyl substituents play a significant role in determining the excited-state dynamics and photoisomerization quantum yields by tuning the relative crowding inside the capsule. The results are discussed in terms of distortions of the ground- and excited-state potential energy surfaces, including the topology of the S₁-S₀ conical intersection.

Photocleavage of coumarin dimers studied by femtosecond UV transient absorption spectroscopy

Coumarins are a class of UV absorbing compounds which exhibit fast, photoinduced cyclobutane ring formation and cleavage reactions. The photophysics behind such processes hold significant relevance for biomedical and photoresponsive materials research. In order to further understand the underlying dynamics of the cleavage reaction, and develop strategies for increasing the reaction efficiency, UV transient absorption spectroscopy was applied to three unsubstituted, isomeric coumarin dimers: anti-head-to-head (anti-hh), syn-head-to-head (syn-hh) and syn-head-to-tail (syn-ht). The experiments performed under 280 nm excitation and broadband (300-620 nm) probing revealed that the cleavage reaction of coumarin dimers occurs through non-radiative, short-lived (<200 fs) singlet states. From the data, two branched kinetic models were developed to describe the monomer formation and dimer relaxation dynamics, identify possible intermediate states, and determine the quantum yields of the dimer cleavage. The anti-hh dimer shows the highest cleavage efficiency with a value of about 20%. The differences in the cleavage efficiency for the three isomers are interpreted in terms of differing steric hindrances of the benzene groups attached to the cyclobutane ring and charge delocalisation of the intermediate state.

π-Conjugated Organometallic Isoindigo Oligomer and Polymer Chromophores: Singlet and Triplet Excited State Dynamics and Application in Polymer Solar Cells

An isoindigo based π-conjugated oligomer and polymer that contain cyclometalated platinum(II) "auxochrome" units were subjected to photophysical characterization, and application of the polymer in bulk heterojunction polymer solar cells with PCBM acceptor was examined. The objective of the study was to explore the effect of the heavy metal centers on the excited state properties, in particular, intersystem crossing to a triplet (exciton) state, and further how this would influence the performance of the organometallic polymer in solar cells. The materials were characterized by electrochemistry, ground state absorption, emission, and picosecond-nanosecond transient absorption spectroscopy. Electrochemical measurements indicate that the cyclometalated units have a significant impact on the HOMO energy level of the chromophores, but little effect on the LUMO, which is consistent with localization of the LUMO on the isoindigo acceptor unit. Picosecond-nanosecond transient absorption spectroscopy reveals a transient with ∼100 ns lifetime that is assigned to a triplet excited state that is produced by intersystem crossing from a singlet state on a time scale of ∼130 ps. This is the first time that a triplet state has been observed for isoindigo π-conjugated chromophores. The performance of the polymer in bulk heterojunction solar cells was explored with PC61BM as an acceptor. The performance of the cells was optimum at a relatively high PCBM loading (1:6, polymer:PCBM), but the overall efficiency was relatively low with power conversion efficiency (PCE) of 0.22%. Atomic force microscopy of blend films reveals that the length scale of the phase separation decreases with increasing PCBM content, suggesting a reason for the increase in PCE with acceptor loading. Energetic considerations show that the triplet state in the polymer is too low in energy to undergo charge separation with PCBM. Further, due to the relatively low LUMO energy of the polymer, charge transfer from the singlet to PCBM is only weakly exothermic, which is believed to be the reason that the photocurrent efficiency is relatively low.

Theoretical study of the hydroxyl radical addition to uracil and photochemistry of the formed U6OH• adduct

Hydroxyl radical ((•)OH) is produced in biological systems by external or endogenous agents. It can damage DNA/RNA by attacking pyrimidine nucleobases through the addition to the C5═C6 double bond. The adduct resulting from the attachment at the C5 position prevails in the experimental measurements, although the reasons for this preference remain unclear. The first aim of this work is therefore to shed light on the comprehension of this important process. Thus, the thermal (•)OH addition to the C5═C6 double bond of uracil has been studied theoretically by using DFT, MP2, and the multiconfigurational CASPT2//CASSCF methodologies. The in-vacuo results obtained with the latter protocol plus the analysis of solvent effects support the experimental observation. A significant lower barrier height is predicted for the C5 pathway with respect to that of the C6 route. In contrast to the C5 adduct, the C6 adduct is able to absorb visible light. Hence, the second aim of the work is to study the photochemistry of this species using the CASPT2//CASSCF methodology within the framework of the photochemical reaction path approach (PRPA). The nonradiative decay to the ground state of this compound has been characterized. A photoreactive character is predicted for the C6 adduct in the excited states according to the presence of excited-state minima along the main decay channel. Finally, a new mechanism of photodissociation has been explored, which implies the photoinduced regeneration of the canonical nucleobase by irradiating with visible light, being therefore relevant in RNA protection against damage by reactive oxygen species.

Energetics and dynamics of the low-lying electronic states of constrained polyenes: implications for infinite polyenes

Steady-state and ultrafast transient absorption spectra were obtained for a series of conformationally constrained, isomerically pure polyenes with 5-23 conjugated double bonds (N). These data and fluorescence spectra of the shorter polyenes reveal the N dependence of the energies of six (1)B(u)(+) and two (1)A(g)(-) excited states. The (1)B(u)(+) states converge to a common infinite polyene limit of 15,900 ± 100 cm(-1). The two excited (1)A(g)(-) states, however, exhibit a large (~9000 cm(-1)) energy difference in the infinite polyene limit, in contrast to the common value previously predicted by theory. EOM-CCSD ab initio and MNDO-PSDCI semiempirical MO theories account for the experimental transition energies and intensities. The complex, multistep dynamics of the 1(1)B(u)(+) → 2(1)A(g)(-) → 1(1)A(g)(-) excited state decay pathways as a function of N are compared with kinetic data from several natural and synthetic carotenoids. Distinctive transient absorption signals in the visible region, previously identified with S* states in carotenoids, also are observed for the longer polyenes. Analysis of the lifetimes of the 2(1)A(g)(-) states, using the energy gap law for nonradiative decay, reveals remarkable similarities in the N dependence of the 2(1)A(g)(-) decay kinetics of the carotenoid and polyene systems. These findings are important for understanding the mechanisms by which carotenoids carry out their roles as light-harvesting molecules and photoprotective agents in biological systems.

Controlling the efficiency of an artificial light-harvesting complex

Adaptive femtosecond pulse shaping in an evolutionary learning loop is applied to a bioinspired dyad molecule that closely mimics the early-time photophysics of the light-harvesting complex 2 (LH2) photosynthetic antenna complex. Control over the branching ratio between the two competing pathways for energy flow, internal conversion (IC) and energy transfer (ET), is realized. We show that by pulse shaping it is possible to increase independently the relative yield of both channels, ET and IC. The optimization results are analyzed by using Fourier analysis, which gives direct insight to the mechanism featuring quantum interference of a low-frequency mode. The results from the closed-loop experiments are repeatable and robust and demonstrate the power of coherent control experiments as a spectroscopic tool (i.e., quantum-control spectroscopy) capable of revealing functionally relevant molecular properties that are hidden from conventional techniques.

Investigation of the Influence of Solute−Solvent Interactions on the Vibrational Energy Relaxation Dynamics of Large Molecules in Liquids

The influence of solute-solvent interactions on the vibrational energy relaxation dynamics of perylene and substituted perylenes in the first singlet excited-state upon excitation with moderate (<0.4 eV) excess energy has been investigated by monitoring the early narrowing of their fluorescence spectrum. This narrowing was found to occur on timescales ranging from a few hundreds of femtoseconds to a few picoseconds. Other processes, such as a partial decay of the fluorescence anisotropy and the damping of a low-frequency oscillation due to the propagation of a vibrational wavepacket, were found to take place on a very similar time scale. No significant relationship between the strength of nonspecific solute-solvent interactions and the vibrational energy relaxation dynamics of the solutes could be evidenced. On the other hand, in alcohols the spectral narrowing is faster with a solute having H-bonding sites, indicating that this specific interaction tends to favor vibrational energy relaxation. No relationship between the dynamics of spectral narrowing and macroscopic solvent properties, such as the thermal diffusivity, could be found. On the other hand, a correlation between this narrowing dynamics and the number of low-frequency modes of the solvent molecules was evidenced. All these observations cannot be discussed with a model where vibrational energy relaxation occurs via two consecutive and dynamically well-separated steps, namely ultrafast intramolecular vibrational redistribution followed by slower vibrational cooling. On the contrary, the results indicate that both intra- and intermolecular vibrational energy redistribution processes are closely entangled and occur, at least partially, on similar timescales.

Forward and Backward Pericyclic Photochemical Reactions Have Intermediates in Common, Yet Cyclobutenes Break the Rules

Photochemical pericyclic reactions are believed to proceed via a so-called pericyclic minimum on the lowest excited potential surface (S(1)), which is common to both the forward and backward reactions. Such a common intermediate has never been directly detected. The photointerconversion of 1,3-butadiene and cyclobutene is the prevailing prototype for such reactions, yet only diene ring closure proceeds with the stereospecificity that the Woodward-Hoffmann rules predict. This contrast seems to exclude a common intermediate. Using ultrafast spectroscopy, we show that the excited states of two cyclobutene/diene isomeric pairs are linked by not one, but by two common minima, p* and ct*. Starting from the diene side (cyclohepta-1,3-diene and cycloocta-1,3-diene), electrocyclic ring closure passes via the pericyclic minimum p*, whereas ct* is mainly responsible for cis-trans isomerization. Starting from the corresponding cyclobutenes (bicyclo[3.2.0]heptene-6 and bicyclo[4.2.0]octene-7), the forbidden isomer is formed from ct*. The path branches at the first (S(2)/S(1)) conical intersection towards p* and ct*. The fact that the energetically unfavorable ct* path can compete is ascribed to a dynamic effect: the momentum in C=C twist direction, acquired--such as in other olefins--in the Franck-Condon region of the cyclobutenes.