Prediction Challenge: Simulating Rydberg photoexcited cyclobutanone with surface hopping dynamics based on different electronic structure methods.

This research examines the nonadiabatic dynamics of cyclobutanone after excitation into the n → 3s Rydberg S2 state. It stems from our contribution to the Special Topic of the Journal of Chemical Physics to test the predictive capability of computational chemistry against unseen experimental data. Decoherence-corrected fewest-switches surface hopping was used to simulate nonadiabatic dynamics with full and approximated nonadiabatic couplings. Several simulation sets were computed with different electronic structure methods, including a multiconfigurational wavefunction [multiconfigurational self-consistent field (MCSCF)] specially built to describe dissociative channels, multireference semiempirical approach, time-dependent density functional theory, algebraic diagrammatic construction, and coupled cluster. MCSCF dynamics predicts a slow deactivation of the S2 state (10 ps), followed by an ultrafast population transfer from S1 to S0 (<100 fs). CO elimination (C3 channel) dominates over C2H4 formation (C2 channel). These findings radically differ from the other methods, which predicted S2 lifetimes 10-250 times shorter and C2 channel predominance. These results suggest that routine electronic structure methods may hold low predictive power for the outcome of nonadiabatic dynamics.

Temperature effects on the internal conversion of excited adenine and adenosine.

This work aims to elucidate the dependence of the excited-state lifetime of adenine and adenosine on temperature. So far, it has been experimentally shown that while adenine's lifetime is unaffected by temperature, adenosine's lifetime strongly depends on it. However, the non-Arrhenius temperature dependence has posed a challenge in explaining this phenomenon. We used surface hopping to simulate the dynamics of adenine and adenosine in the gas phase at 0 and 400 K. The temperature effects were observed under the initial conditions via Wigner sampling with thermal corrections. Our results confirm that adenine's excited-state lifetime does not depend on temperature, while adenosine's lifetime does. Adenosine's dependency is due to intramolecular vibrational energy transfer from adenine to the ribose group. At 0 K, this transfer reduced the mean kinetic energy of adenine's moiety so much that internal conversion is inhibited, and the lifetime elongated by a factor of 2.3 compared to that at 400 K. The modeling also definitively ruled out the influence of viscosity, which was proposed as an alternative explanation previously.

Surface hopping modeling of charge and energy transfer in active environments.

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

Using diketopyrrolopyrroles to stabilize double excitation and control internal conversion.

Diketopyrrolopyrrole (DPP) is a pivotal functional group to tune the physicochemical properties of novel organic photoelectronic materials. Among multiple uses, DPP-thiophene derivatives forming a dimer through a vinyl linker were recently shown to quench the fluorescence observed in their isolated monomers. Here, we explain this fluorescence quenching using computational chemistry. The DPP-thiophene dimer has a low-lying doubly excited state that is not energetically accessible for the monomer. This state delays the fluorescence allowing internal conversion to occur first. We characterize the doubly excited state wavefunction by systematically changing the derivatives to tune the π-scaffold size and the acceptor and donor characters. The origin of this state's stabilization is related to the increase in the π-system and not to the charge-transfer features. This analysis delivers core conceptual information on the electronic properties of organic chromophores arranged symmetrically around a vinyl linker, opening new ways to control the balance between luminescence and internal conversion.

Probing the electronic structure and photophysics of thiophene-diketopyrrolopyrrole derivatives in solution.

Diketopyrrolopyrroles are a popular class of electron-withdrawing unit in optoelectronic materials. When combined with electron donating side-chain functional groups such as thiophenes, they form a very broad class of donor-acceptor molecules: thiophene-diketopyrrolopyrroles (TDPPs). Despite their widescale use in biosensors and photovoltaic materials, studies have yet to establish the important link between the electronic structure of the specific TDPP and the critical optical properties. To bridge this gap, ultrafast transient absorption with 22 fs time resolution has been used to explore the photophysics of three prototypical TDPP molecules: a monomer, dimer and polymer in solution. Interpretation of experimental data was assisted by a recent high-level theoretical study, and additional density functional theory calculations. These studies show that the photophysics of these molecular prototypes under visible photoexcitation are determined by just two excited electronic states, having very different electronic characters (one is optically bright, the other dark), their relative energetic ordering and the timescales for internal conversion from one to the other and/or to the ground state. The underlying difference in electronic structure alters the branching between these excited states and their associated dynamics. In turn, these factors dictate the fluorescence quantum yields, which are shown to vary by ∼1-2 orders of magnitude across the TDPP prototypes investigated here. The fast non-radiative transfer of molecules from the bright to dark states is mediated by conical intersections. Remarkably, wavepacket signals in the measured transient absorption data carry signatures of the nuclear motions that enable mixing of the electronic-nuclear wavefunction and facilitate non-adiabatic coupling between the bright and dark states.

Mechanistic Aspects of the Photophysics of UVA Filters Based on Meldrum Derivatives.

Skin photoprotection against UVA radiation is crucial, but it is hindered by the sparsity of approved commercial UVA filters. Sinapoyl malate (SM) derivatives are promising candidates for a new class of UVA filters. They have been previously identified as an efficient photoprotective sunscreen in plants due to their fast nonradiative energy dissipation. Combining experimental and computational results, in our previous letter (J. Phys. Chem. Lett. 2021, 12, 337-344) we showed that coumaryl Meldrum (CMe) and sinapoyl Meldrum (SMe) are outstanding candidates for UVA filters in sunscreen formulations. Here, we deliver a comprehensive computational characterization of the excited-state dynamics of these molecules. Using reaction pathways and excited-state dynamics simulations, we could elucidate the photodeactivation mechanism of these molecules. Upon photoexcitation, they follow a two-step logistic decay. First, an ultrafast and efficient relaxation stabilizes the excited state alongside a 90° twisting around the allylic double bond, giving rise to a minimum with a twisted intramolecular excited-state (TICT) character. From this minimum, internal conversion to the ground state occurs after overcoming a 0.2 eV barrier. Minor differences in the nonradiative decay and fluorescence of CMe and SMe are associated with an additional minimum present only in the latter.

A three-state model for the photo-Fries rearrangement.

A three-state model for the photo-Fries rearrangement (PFR) is proposed based on multiconfigurational calculations. It provides a comprehensive mechanistic picture of all steps of the reaction, from the photoabsorption to the final tautomerization. The three states participating in the PFR are an aromatic 1ππ*, which absorbs the radiation; a pre-dissociative 1nπ*, which transfers the energy to the dissociative region; and a 1πσ*, along which dissociation occurs. The transfer from 1ππ* to 1nπ* involves pyramidalization of the carbonyl carbon, while transfer from 1nπ* to 1πσ* takes place through CO stretching. Different products are available after a conical intersection with the ground state. Among them is a recombined radical intermediate, which can yield ortho-PFR products after an intramolecular 1,3-H tunneling. The three-state model is developed for phenyl acetate, the basic prototype for the PFR, and it reconciles the theory with a series of observations from time-resolved spectroscopy. It also delivers a rational way to optimize PFR yields, since, as shown for four different systems, diverse substituents can change the energetic order of the 1ππ* and 1nπ* states, preventing or enhancing the PFR.

Recommendations for Velocity Adjustment in Surface Hopping.

This study investigates velocity adjustment directions after hopping in surface hopping dynamics. Using fulvene and a protonated Schiff base (PSB4) as case studies, we investigate the population decay and reaction yields of different sets of dynamics with the velocity adjusted in either the nonadiabatic coupling, gradient difference, or momentum directions. For the latter, in addition to the conventional algorithm, we investigated the performance of a reduced kinetic energy reservoir approach recently proposed. Our evaluation also considered velocity adjustment in the directions of approximate nonadiabatic coupling vectors. While results for fulvene are susceptible to the adjustment approach, PSB4 is not. We correlated this dependence to the topography near the conical intersections. When nonadiabatic coupling vectors are unavailable, the gradient difference direction is the best adjustment option. If the gradient difference is also unavailable, a semiempirical vector direction or the momentum direction with a reduced kinetic energy reservoir becomes an excellent option to prevent an artificial excess of back hoppings. The precise velocity adjustment direction is less crucial for describing the nonadiabatic dynamics than the kinetic energy reservoir's size.

Classification of doubly excited molecular electronic states.

Electronic states with partial or complete doubly excited character play a crucial role in many areas, such as singlet fission and non-linear optical spectroscopy. Although doubly excited states have been studied in polyenes and related systems for many years, the assignment as singly vs. doubly excited, even in the simplest case of butadiene, has sparked controversies. So far, no well-defined framework for classifying doubly excited states has been developed, and even more, there is not even a well-accepted definition of doubly excited character as such. Here, we present a solution: a physically motivated definition of doubly excited character based on operator expectation values and density matrices, which works independently of the underlying orbital representation, avoiding ambiguities that have plagued earlier studies. Furthermore, we propose a classification scheme to differentiate three cases: (i) two single excitations occurring within two independent pairs of orbitals leaving four open shells (DOS), (ii) the promotion of both electrons to the same orbital, producing a closed-shell determinant (DCS), and (iii) a mixture of singly and doubly excited configurations not aligning with either one of the previous cases (Dmix). We highlight their differences in underlying energy terms and explain their signatures in practical computations. The three cases are illustrated through various high-level computational methods using dimers for DOS, polyenes for Dmix, and cyclobutane and tetrazine for DCS. The conversion between DOS and DCS is investigated using a well-known photochemical reaction, the photodimerization of ethylene. This work provides a deeper understanding of doubly excited states and may guide more rigorous discussions toward improving their computational description while also giving insight into their fundamental photophysics.

Understanding the Impact of Symmetrical Substitution on the Photodynamics of Sinapate Esters Using Gas-Phase Ultrafast Spectroscopy.

Two model biomimetic systems, ethyl sinapate (ES) and its symmetrical analogue, diethyl 2-(4-hydroxy-3,5-dimethoxybenzylidene)malonate (or diethyl sinapate, DES), are stripped to their core fundamentals through gas-phase spectroscopy to understand the underlying photophysics of photothermal materials. Following photoexcitation to the optically bright S1(ππ*) state, DES is found to repopulate the electronic ground state over 3 orders of magnitude quicker than its nonsymmetrical counterpart, ES. Our XMS-CASPT2 calculations shed light on the experimental results, revealing crucial differences in the potential energy surfaces and conical intersection topography between ES and DES. From this work, a peaked conical intersection, seen for DES, shows vital importance for the nonradiative ground-state recovery of photothermal materials. This fundamental comparative study highlights the potential impact that symmetrical substitution can have on the photodynamics of sinapate esters, providing a blueprint for future advancement in photothermal technology.

Role of the Hydrogen Bond on the Internal Conversion of Photoexcited Adenosine.

Experiments and theory have revealed that hydrogen bonds modify the excited-state lifetimes of nucleosides compared to nucleobases. Nevertheless, how these bonds impact the internal conversion is still unsettled. This work simulates the non-adiabatic dynamics of adenosine conformers in the gas phase with and without hydrogen bonds between the sugar and adenine moieties. The isomer containing the hydrogen bond (syn) exhibits a significantly shorter excited-state lifetime than the isomer without it (anti). However, internal conversion through electron-driven proton transfer between sugar and adenine plays only a minor (although non-negligible) role in the photophysics of adenosine. Either with or without hydrogen bonds, photodeactivation preferentially occurs following the ring-puckering pathways. The role of the hydrogen bond is to avoid the sugar rotation relative to adenine, shortening the distance to the ring-puckering internal conversion.

Direct structural observation of ultrafast photoisomerization dynamics in sinapate esters.

Sinapate esters have been extensively studied for their potential application in 'nature-inspired' photoprotection. There is general consensus that the relaxation mechanism of sinapate esters following photoexcitation with ultraviolet radiation is mediated by geometric isomerization. This has been largely inferred through indirect studies involving transient electronic absorption spectroscopy in conjunction with steady-state spectroscopies. However, to-date, there is no direct experimental evidence tracking the formation of the photoisomer in real-time. Using transient vibrational absorption spectroscopy, we report on the direct structural changes that occur upon photoexcitation, resulting in the photoisomer formation. Our mechanistic analysis predicts that, from the photoprepared ππ* state, internal conversion takes place through a conical intersection (CI) near the geometry of the initial isomer. Our calculations suggest that different CI topographies at relevant points on the seam of intersection may influence the isomerization yield. Altogether, we provide compelling evidence suggesting that a sinapate ester's geometric isomerization can be a more complex dynamical process than originally thought.

Fewest switches surface hopping with Baeck-An couplings.

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

New Generation UV-A Filters: Understanding Their Photodynamics on a Human Skin Mimic.

The sparsity of efficient commercial ultraviolet-A (UV-A) filters is a major challenge toward developing effective broadband sunscreens with minimal human- and eco-toxicity. To combat this, we have designed a new class of Meldrum-based phenolic UV-A filters. We explore the ultrafast photodynamics of coumaryl Meldrum, CMe, and sinapyl Meldrum (SMe), both in an industry-standard emollient and on a synthetic skin mimic, using femtosecond transient electronic and vibrational absorption spectroscopies and computational simulations. Upon photoexcitation to the lowest excited singlet state (S1), these Meldrum-based phenolics undergo fast and efficient nonradiative decay to repopulate the electronic ground state (S0). We propose an initial ultrafast twisted intramolecular charge-transfer mechanism as these systems evolve out of the Franck-Condon region toward an S1/S0 conical intersection, followed by internal conversion to S0 and subsequent vibrational cooling. Importantly, we correlate these findings to their long-term photostability upon irradiation with a solar simulator and conclude that these molecules surpass the basic requirements of an industry-standard UV filter.

Towards developing novel and sustainable molecular light-to-heat converters.

Light-to-heat conversion materials generate great interest due to their widespread applications, notable exemplars being solar energy harvesting and photoprotection. Another more recently identified potential application for such materials is in molecular heaters for agriculture, whose function is to protect crops from extreme cold weather and extend both the growing season and the geographic areas capable of supporting growth, all of which could help reduce food security challenges. To address this demand, a new series of phenolic-based barbituric absorbers of ultraviolet (UV) radiation has been designed and synthesised in a sustainable manner. The photophysics of these molecules has been studied in solution using femtosecond transient electronic and vibrational absorption spectroscopies, allied with computational simulations and their potential toxicity assessed by in silico studies. Following photoexcitation to the lowest singlet excited state, these barbituric absorbers repopulate the electronic ground state with high fidelity on an ultrafast time scale (within a few picoseconds). The energy relaxation pathway includes a twisted intramolecular charge-transfer state as the system evolves out of the Franck-Condon region, internal conversion to the ground electronic state, and subsequent vibrational cooling. These barbituric absorbers display promising light-to-heat conversion capabilities, are predicted to be non-toxic, and demand further study within neighbouring application-based fields.

On the origin of the shift between vertical excitation and band maximum in molecular photoabsorption.

The analysis of the photoabsorption spectra of molecules shows that the band maximum is usually redshifted in comparison to the vertical excitation. We conducted a throughout analysis of this shift based on low-dimensional analytical and numerical model systems, showing that its origin is rooted in the frequency change between the ground and the excited states in multidimensional systems. Moreover, we deliver a benchmark of ab initio results for the shift based on a comparison of vertical excitations and band maxima calculated with the nuclear ensemble approach for the 28 organic molecules in the Mülheim molecular dataset. The mean value of the shift calculated over 60 transitions is 0.11 ± 0.08 eV. The mean value of the band width is 0.32 ± 0.14 eV. Graphical abstract .

Correction: Dipolar vinyl sulfur fluorescent dyes. Synthesis and photophysics of sulfide, sulfoxide and sulfone based D-π-A compounds.

[This corrects the article DOI: 10.1039/C6RA27989A.].