Observation of the Sub-100 Femtosecond Population of a Dark State in a Thiobase Mediating Intersystem Crossing

Transient-absorption spectroscopy of dendrimers via nonadiabatic excited-state dynamics simulations

The efficiency of light-harvesting and energy transfer in multi-chromophore ensembles underpins natural photosynthesis. Dendrimers are highly branched synthetic multi-chromophoric conjugated supra-molecules that mimic these natural processes. After photoexcitation, their repeated units participate in a number of intramolecular electronic energy relaxation and redistribution pathways that ultimately funnel to a sink. Here, a model four-branched dendrimer with a pyrene core is theoretically studied using nonadiabatic molecular dynamics simulations. We evaluate excited-state photoinduced dynamics of the dendrimer, and demonstrate on-the-fly simulations of its transient absorption pump-probe (TA-PP) spectra. We show how the evolutions of the simulated TA-PP spectra monitor in real time photoinduced energy relaxation and redistribution, and provide a detailed microscopic picture of the relevant energy-transfer pathways. To the best of our knowledge, this is the first of this kind of on-the-fly atomistic simulation of TA-PP signals reported for a large molecular system.

Thio and Seleno-Psoralens as Efficient Triplet Harvesting Photosensitizers for Photodynamic Therapy

The Psoralen (Pso) molecule finds extensive applications in photo-chemotherapy, courtesy of its triplet state forming ability. Sulfur and selenium replacement of exocyclic carbonyl oxygen of organic chromophores foster efficient triplet harvesting with near unity triplet quantum yield. These triplet-forming photosensitizers are useful in Photodynamic Therapy (PDT) applications for selective apoptosis of cancer cells. In this work, we have critically assessed the effect of the sulfur and selenium substitution at the exocyclic carbonyl (TPso and SePso, respectively) and endocyclic oxygen positions of Psoralen. It resulted in a significant redshifted absorption spectrum to access the PDT therapeutic window with increased oscillator strength. The reduction in singlet-triplet energy gap and enhancement in the spin-orbit coupling values increase the number of intersystem crossing (ISC) pathways to the triplet manifold, which shortens the ISC lifetime from 10-5 s for Pso to 10-8 s for TPso and 10-9 s for SePso. The intramolecular photo-induced electron transfer process, a competitive pathway to ISC, is also considerably curbed by exocyclic functionalizations. In addition, a maximum of 115 GM of two-photon absorption (2PA) with IR absorption (660-1050 nm) confirms that the Psoralen skeleton can be effectively tweaked via single chalcogen atom replacement to design a suitable PDT photosensitizer.

Resolving Photoinduced Femtosecond Three-Dimensional Solute-Solvent Dynamics through Surface Hopping Simulations

Photoinduced dynamics in solution is governed by mutual solute-solvent interactions, which give rise to phenomena like solvatochromism, the Stokes shift, dual fluorescence, or charge transfer. Understanding these phenomena requires simulating the solute's photoinduced dynamics and simultaneously resolving the three-dimensional solvent distribution dynamics. If using trajectory surface hopping (TSH) to this aim, thousands of trajectories are required to adequately sample the time-dependent three-dimensional solvent distribution functions, and thus resolve the solvent dynamics with sub-Ångstrom and femtosecond accuracy and sufficiently low noise levels. Unfortunately, simulating thousands of trajectories with TSH in the framework of hybrid quantum mechanical/molecular mechanical (QM/MM) can be prohibitively expensive when employing ab initio electronic structure methods. To tackle this challenge, we recently introduced a computationally efficient approach that combines efficient linear vibronic coupling models with molecular mechanics (LVC/MM) via electrostatic embedding [Polonius et al., JCTC 2023, 19, 7171-7186]. This method provides solvent-embedded, nonadiabatically coupled potential energy surfaces while scaling similarly to MM force fields. Here, we employ TSH with LVC/MM to unravel the photoinduced dynamics of two small thiocarbonyl compounds solvated in water. We describe how to estimate the number of trajectories required to produce nearly noise-free three-dimensional solvent distribution functions and present an analysis based on approximately 10,000 trajectories propagated for 3 ps. In the electronic ground state, both molecules exhibit in-plane hydrogen bonds to the sulfur atom. Shortly after excitation, these bonds are broken and reform perpendicular to the molecular plane on timescales that differ by an order of magnitude due to steric effects. We also show that the solvent relaxation dynamics is coupled to the electronic dynamics, including intersystem crossing. These findings are relevant to advance the understanding of the coupled solute-solvent dynamics of solvated photoexcited molecules, e.g., biologically relevant thio-nucleobases.

WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

Sub-100-fs energy transfer in coenzyme NADH is a coherent process assisted by a charge-transfer state

Excitation energy transfer (EET) is a key photoinduced process in biological chromophoric assemblies. Here we investigate the factors which can drive EET into efficient ultrafast sub-ps regimes. We demonstrate how a coherent transport of electronic population could facilitate this in water solvated NADH coenzyme and uncover the role of an intermediate dark charge-transfer state. High temporal resolution ultrafast optical spectroscopy gives a 54±11 fs time constant for the EET process. Nonadiabatic quantum dynamical simulations computed through the time-evolution of multidimensional wavepackets suggest that the population transfer is mediated by photoexcited molecular vibrations due to strong coupling between the electronic states. The polar aqueous solvent environment leads to the active participation of a dark charge transfer state, accelerating the vibronically coherent EET process in favorably stacked conformers and solvent cavities. Our work demonstrates how the interplay of structural and environmental factors leads to diverse pathways for the EET process in flexible heterodimers and provides general insights relevant for coherent EET processes in stacked multichromophoric aggregates like DNA strands.

Ultrafast photochemistry and electron-diffraction spectra in n → (3s) Rydberg excited cyclobutanone resolved at the multireference perturbative level

We study the ultrafast time evolution of cyclobutanone excited to the singlet n → Rydberg state through non-adiabatic surface-hopping simulationsperformed at extended multi-state complete active space second-order perturbation (XMS-CASPT2) level of theory. These dynamics predict relaxation to the ground-state with a timescale of 822 ± 45 fs with minimal involvement of the triplets. The major relaxation path to the ground-state involves a three-state degeneracy region and leads to a variety of fragmented photoproducts. We simulate the resulting time-resolved electron-diffraction spectra, which track the relaxation of the excited state and the formation of various photoproducts in the ground state.

Quantum-Chemistry Study of the Photophysical Properties of 4-Thiouracil and Comparisons with 2-Thiouracil

DNA in living beings is constantly damaged by exogenous and endogenous agents. However, in some cases, DNA photodamage can have interesting applications, as it happens in photodynamic therapy. In this work, the current knowledge on the photophysics of 4-thiouracil has been extended by further quantum-chemistry studies to improve the agreement between theory and experiments, to better understand the differences with 2-thiouracil, and, last but not least, to verify its usefulness as a photosensitizer for photodynamic therapy. This study has been carried out by determining the most favorable deactivation paths of UV-vis photoexcited 4-thiouracil by means of the photochemical reaction path approach and an efficient combination of the complete-active-space second-order perturbation theory//complete-active-space self-consistent field (CASPT2//CASSCF), (CASPT2//CASPT2), time-dependent density functional theory (TDDFT), and spin-flip TDDFT (SF-TDDFT) methodologies. By comparing the data computed herein for both 4-thiouracil and 2-thiouracil, a rationale is provided on the relatively higher yields of intersystem crossing, triplet lifetime and singlet oxygen production of 4-thiouracil, and the relatively higher yield of phosphorescence of 2-thiouracil.

Toward a Stochastic Complete Active Space Second-Order Perturbation Theory

In this work, an internally contracted stochastic complete active space second-order perturbation theory, stochastic-CASPT2, is reported. The method relies on stochastically sampled reduced density matrices (RDMs) up to rank four and contractions thereof with the generalized Fock matrix. A new protocol for calculating higher-order RDMs in full configuration interaction quantum Monte Carlo (FCIQMC) has been designed based on (1) restricting sampling of the corresponding excitations to a deterministic subspace, (2) averaging the RDMs from independent dynamics and (3) projecting them onto the closest positive semi-definite matrix. Our protocol avoids previously encountered numerical conditioning problems in the orthogonalization of the perturber overlap matrix stemming from numerical noise. The chromium dimer CASSCF(12,12)/CASPT2 binding curve is computed as a proof of concept.

The impact of the chalcogen-substitution element and initial spectroscopic state on excited-state relaxation pathways in nucleobase photosensitizers: a combination of static and dynamic studies

The substitution of oxygen with chalcogen in carbonyl group(s) of canonical nucleobases gives an impressive triplet generation, enabling their promising applications in medicine and other emerging techniques. The excited-state relaxation S2(ππ*) → S1(nπ*) → T1(ππ*) has been considered the preferred path for triplet generation in these nucleobase derivatives. Here, we demonstrate enhanced quantum efficiency of direct intersystem crossing from S2 to triplet manifold upon substitution with heavier chalcogen elements. The excited-state relaxation dynamics of sulfur/selenium substituted guanines in a vacuum is investigated using a combination of static quantum chemical calculations and on-the-fly excited-state molecular dynamics simulations. We find that in sulfur-substitution the S2 state predominantly decays to the S1 state, while upon selenium-substitution the S2 state deactivation leads to simultaneous population of the S1 and T2,3 states in the same time scale and multi-state quasi-degeneracy region S2/S1/T2,3. Interestingly, the ultrafast deactivation of the spectroscopic S3 state of both studied molecules to the S1 state occurs through a successive S3 → S2 → S1 path involving a multi-state quasi-degeneracy S3/S2/S1. The populated S1 and T2 states will cross the lowest triplet state, and the S1 → T intersystem crossing happens in a multi-state quasi-degeneracy region S1/T2,3/T1 and is accelerated by selenium-substitution. The present study reveals the influence of both the chalcogen substitution element and initial spectroscopic state on the excited-state relaxation mechanism of nucleobase photosensitizers and also highlights the important role of multi-state quasi-degeneracy in mediating the complex relaxation process. These theoretical results provide additional insights into the intrinsic photophysics of nucleobase-based photosensitizers and are helpful for designing novel photo-sensitizers for real applications.

Thieno[3,4-d]pyrimidin-4(3H)-thione: an effective, oxygenation independent, heavy-atom-free photosensitizer for cancer cells

All-organic, heavy-atom-free photosensitizers based on thionation of nucleobases are receiving increased attention because they are easy to make, noncytotoxic, work both in the presence and absence of molecular oxygen, and can be readily incorporated into DNA and RNA. In this contribution, the DNA and RNA fluorescent probe, thieno[3,4-d]pyrimidin-4(1H)-one, has been thionated to develop thieno[3,4-d]pyrimidin-4(3H)-thione, which is nonfluorescent and absorbs near-visible radiation with about 60% higher efficiency. Steady-state absorption and emission spectra are combined with transient absorption spectroscopy and CASPT2 calculations to delineate the electronic relaxation mechanisms of both pyrimidine derivatives in aqueous and acetonitrile solutions. It is demonstrated that thieno[3,4-d]pyrimidin-4(3H)-thione efficiently populates the long-lived and reactive triplet state generating singlet oxygen with a quantum yield of about 80% independent of solvent. It is further shown that thieno[3,4-d]pyrimidin-4(3H)-thione exhibits high photodynamic efficacy against monolayer melanoma cells and cervical cancer cells both under normoxic and hypoxic conditions. Our combined spectroscopic, computational, and in vitro data demonstrate the excellent potential of thieno[3,4-d]pyrimidin-4(1H)-thione as a heavy-atom-free PDT agent and paves the way for further development of photosensitizers based on the thionation of thieno[3,4-d]pyrimidine derivatives. Collectively, the experimental and computational results demonstrate that thieno[3,4-d]pyrimidine-4(3H)-thione stands out as the most promising thiobase photosensitizer developed to this date.

Direct temporal characterization of sub-3-fs deep UV pulses generated by resonant dispersive wave emission

We report on the complete temporal characterization of ultrashort pulses, generated by resonant dispersive wave emission in gas-filled hollow-capillary fibers, with energy in the microjoule range and continuously tunable from the deep-ultraviolet to the ultraviolet. Temporal characterization of such ultrabroad pulses, particularly challenging in this spectral region, was performed using an all-in-vacuum setup for self-diffraction frequency resolved optical gating (SD-FROG). Sub-3-fs pulses were measured, tunable from 250 nm to 350 nm, with a minimum pulse duration of 2.4 ± 0.1 fs.

Probing C-I bond fission in the UV photochemistry of 2-iodothiophene with core-to-valence transient absorption spectroscopy

The UV photochemistry of small heteroaromatic molecules serves as a testbed for understanding fundamental photo-induced chemical transformations in moderately complex compounds, including isomerization, ring-opening, and molecular dissociation. Here, a combined experimental-theoretical study of 268 nm UV light-induced dynamics in 2-iodothiophene (C4H3IS) is performed. The dynamics are experimentally monitored with a femtosecond extreme ultraviolet (XUV) probe that measures iodine N-edge 4d core-to-valence transitions. Experiments are complemented by density functional theory calculations of both the pump-pulse induced valence excitations and the XUV probe-induced core-to-valence transitions. Possible intramolecular relaxation dynamics are investigated by ab initio molecular dynamics simulations. Gradual absorption changes up to ∼0.5 to 1 ps after excitation are observed for both the parent molecular species and emerging iodine fragments, with the latter appearing with a characteristic rise time of 160 ± 30 fs. Comparison of spectral intensities and energies with the calculations identifies an iodine dissociation pathway initiated by a predominant π → π* excitation. In contrast, initial excitation to a nearby n⟂ → σ* state appears unlikely based on a significantly smaller oscillator strength and the absence of any corresponding XUV absorption signatures. Excitation to the π → π* state is followed by contraction of the C-I bond, enabling a nonadiabatic transition to a dissociative π→σC-I* state. For the subsequent fragmentation, a relatively narrow bond-length region along the C-I stretch coordinate between 230 and 280 pm is identified, where the transition between the parent molecule and the thienyl radical + iodine atom products becomes prominent in the XUV spectrum due to rapid localization of two singly occupied molecular orbitals on the two fragments.

Nonlinear Molecular Electronic Spectroscopy via MCTDH Quantum Dynamics: From Exact to Approximate Expressions

We present an accurate and efficient approach to computing the linear and nonlinear optical spectroscopy of a closed quantum system subject to impulsive interactions with an incident electromagnetic field. It incorporates the effect of ultrafast nonadiabatic dynamics by means of explicit numerical propagation of the nuclear wave packet. The fundamental expressions for the evaluation of first- and higher-order response functions are recast in a general form that can be used with any quantum dynamics code capable of computing the overlap of nuclear wave packets evolving in different states. Here we present the evaluation of these expressions with the multiconfiguration time-dependent Hartree (MCTDH) method. Application is made to pyrene, excited to its lowest bright excited state S₂ which exhibits a sub-100-fs nonadiabatic decay to a dark state S₁. The system is described by a linear vibronic coupling Hamiltonian, parametrized with multiconfiguration electronic structure methods. We show that the ultrafast nonadiabatic dynamics can have a remarkable effect on the spectral line shapes that goes beyond simple lifetime broadening. Furthermore, a widely employed approximate expression based on the time scale separation of dephasing and population relaxation is recast in the same theoretical framework. Application to pyrene shows the range of validity of such approximations.

Thiolumazines as Heavy-Atom-Free Photosensitizers for Applications in Daylight Photodynamic Therapy: Insights from Ultrafast Excited-State Dynamics

Finding appropriate photosensitizers (PSs) for daylight photodynamic therapy (dPDT) applications is extremely challenging, even though heavy-atom-free photosensitizers (HAFPSs) such as thiocarbonyl-modified nucleobases have shown a ray of hope. Few attempts have been made to find alternative natural products for dPDT applications. Pteridine heterocycles consisting of a pyrazine ring and a pyrimidine ring, such as lumazine, which exhibit many structural similarities to the alloxazine ring of the flavin molecule, could be an option for HAFPSs. The photophysical and quantum mechanical studies of the thio-modified lumazines revealed that sequential thiomodifications in lumazine result in a bathochromic shift. Additionally, higher tissue penetration depths were observed for thiolumazines. The fluorescence quenching in the case of thiomodified lumazines was explained using triplet state formation, whereas the contribution from the photoinduced electron transfer process cannot be ignored. It was also noticed that a strong one-photon absorption influenced the two-photon absorption (TPA) process, leading to a self-focusing effect in the visible spectral region. The higher tissue penetration and larger TPA cross section are the hallmark characteristics of the thiolumazines to be considered as potential HAFPSs for dPDT applications.

Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water

Rationalizing the photochemistry of nucleobases where an oxygen is replaced by a heavier atom is essential for applications that exploit near-unity triplet quantum yields. Herein, we report on the ultrafast excited-state deactivation mechanism of 6-selenoguanine (6SeGua) in water by combining nonadiabatic trajectory surface-hopping dynamics with an electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. We find that the predominant relaxation mechanism after irradiation starts on the bright singlet S₂ state that converts internally to the dark S₁ state, from which the population is transferred to the triplet T₂ state via intersystem crossing and finally to the lowest T₁ state. This S₂ → S₁ → T₂ → T₁ deactivation pathway is similar to that observed for the lighter 6-thioguanine (6tGua) analogue, but counterintuitively, the T₁ lifetime of the heavier 6SeGua is shorter than that of 6tGua. This fact is explained by the smaller activation barrier to reach the T₁/S₀ crossing point and the larger spin-orbit couplings of 6SeGua compared to 6tGua. From the dynamical simulations, we also calculate transient absorption spectra (TAS), which provide two time constants (τ₁ = 131 fs and τ₂ = 191 fs) that are in excellent agreement with the experimentally reported value (τ_exp = 130 ± 50 fs) (Farrel et al. J. Am. Chem. Soc. 2018, 140, 11214). Intersystem crossing itself is calculated to occur with a time scale of 452 ± 38 fs, highlighting that the TAS is the result of a complex average of signals coming from different nonradiative processes and not intersystem crossing alone.

Attosecond spectroscopy for the investigation of ultrafast dynamics in atomic, molecular and solid-state physics

Since the first demonstration of the generation of attosecond pulses (1 as = 10^-18s) in the extreme-ultraviolet spectral region, several measurement techniques have been introduced, at the beginning for the temporal characterization of the pulses, and immediately after for the investigation of electronic and nuclear ultrafast dynamics in atoms, molecules and solids with unprecedented temporal resolution. The attosecond spectroscopic tools established in the last two decades, together with the development of sophisticated theoretical methods for the interpretation of the experimental outcomes, allowed to unravel and investigate physical processes never observed before, such as the delay in photoemission from atoms and solids, the motion of electrons in molecules after prompt ionization which precede any notable nuclear motion, the temporal evolution of the tunneling process in dielectrics, and many others. This review focused on applications of attosecond techniques to the investigation of ultrafast processes in atoms, molecules and solids. Thanks to the introduction and ongoing developments of new spectroscopic techniques, the attosecond science is rapidly moving towards the investigation, understanding and control of coupled electron-nuclear dynamics in increasingly complex systems, with ever more accurate and complete investigation techniques. Here we will review the most common techniques presenting the latest results in atoms, molecules and solids.

Time-Resolved Optical Pump-Resonant X-ray Probe Spectroscopy of 4-Thiouracil: A Simulation Study

We theoretically monitor the photoinduced ππ* → nπ* internal conversion process in 4-thiouracil (4TU), triggered by an optical pump. The element-sensitive spectroscopic signatures are recorded by a resonant X-ray probe tuned to the sulfur, oxygen, or nitrogen K-edge. We employ high-level electronic structure methods optimized for core-excited electronic structure calculation combined with quantum nuclear wavepacket dynamics computed on two relevant nuclear modes, fully accounting for their quantum nature of nuclear motions. We critically discuss the capabilities and limitations of the resonant technique. For sulfur and nitrogen, we document a pre-edge spectral window free from ground-state background and rich with ππ* and nπ* absorption features. The lowest sulfur K-edge shows strong absorption for both ππ* and nπ*. In the lowest nitrogen K-edge window, we resolve a state-specific fingerprint of the ππ* and an approximate timing of the conical intersection via its depletion. A spectral signature of the nπ* transition, not accessible by UV-vis spectroscopy, is identified. The oxygen K-edge is not sensitive to molecular deformations and gives steady transient absorption features without spectral dynamics. The ππ*/nπ* coherence information is masked by more intense contributions from populations. Altogether, element-specific time-resolved resonant X-ray spectroscopy provides a detailed picture of the electronic excited-state dynamics and therefore a sensitive window into the photophysics of thiobases.

Ultrafast Excited-State Decay Mechanisms of 6-Thioguanine Followed by Sub-20 fs UV Transient Absorption Spectroscopy

Understanding the primary steps following UV photoexcitation in sulphur-substituted DNA bases (thiobases) is fundamental for developing new phototherapeutic drugs. However, the investigation of the excited-state dynamics in sub-100 fs time scales has been elusive until now due to technical challenges. Here, we track the ultrafast decay mechanisms that lead to the electron trapping in the triplet manifold for 6-thioguanine in an aqueous solution, using broadband transient absorption spectroscopy with a sub-20 fs temporal resolution. We obtain experimental evidence of the fast internal conversion from the S₂(ππ*) to the S₁(nπ*) states, which takes place in about 80 fs and demonstrates that the S₁(nπ*) state acts as a doorway to the triplet population in 522 fs. Our results are supported by MS-CASPT2 calculations, predicting a planar S₂(ππ*) pseudo-minimum in agreement with the stimulated emission signal observed in the experiment.

Coupled- and Independent-Trajectory Approaches Based on the Exact Factorization Using the PyUNIxMD Package

We present mixed quantum-classical approaches based on the exact factorization framework. The electron-nuclear correlation term in the exact factorization enables us to deal with quantum coherences by accounting for electronic and nuclear nonadiabatic couplings effectively within classical nuclei approximation. We compare coupled- and independent-trajectory approximations with each other to understand algorithms in description of the bifurcation of nuclear wave packets and the correct spatial distribution of electronic wave functions along with nuclear trajectories. Finally, we show numerical results for comparisons of coupled- and independent-trajectory approaches for the photoisomerization of a protonated Schiff base from excited state molecular dynamics (ESMD) simulations with the recently developed Python-based ESMD code, namely, the PyUNIxMD program.

Following excited-state chemical shifts in molecular ultrafast x-ray photoelectron spectroscopy

The conversion of photon energy into other energetic forms in molecules is accompanied by charge moving on ultrafast timescales. We directly observe the charge motion at a specific site in an electronically excited molecule using time-resolved x-ray photoelectron spectroscopy (TR-XPS). We extend the concept of static chemical shift from conventional XPS by the excited-state chemical shift (ESCS), which is connected to the charge in the framework of a potential model. This allows us to invert TR-XPS spectra to the dynamic charge at a specific atom. We demonstrate the power of TR-XPS by using sulphur 2p-core-electron-emission probing to study the UV-excited dynamics of 2-thiouracil. The method allows us to discover that a major part of the population relaxes to the molecular ground state within 220–250 fs. In addition, a 250-fs oscillation, visible in the kinetic energy of the TR-XPS, reveals a coherent exchange of population among electronic states.

Imaging the charge flow in photoexcited molecules would provide key information on photophysical and photochemical processes. Here the authors demonstrate tracking in real time after photoexcitation the change in charge density at a specific site of 2-thiouracil using time-resolved X-ray photoelectron spectroscopy.

One order of magnitude increase of triplet state lifetime observed in deprotonated form selenium substituted uracil

Selenium nucleic acids possess unique properties and have been demonstrated to have a wide range of applications such as in DNA X-ray crystallography and novel medical therapies. However, as a heavy atom, selenium substitution may easily alter the photophysical properties of a nucleic acid by red-shifting the absorption spectra and introducing effective intersystem crossing to triplet excited states. In present work, the excited state dynamics of a naturally occurring selenium substituted uracil (2-selenuracil, 2SeU) is studied by using femtosecond transient absorption spectroscopy as well as quantum chemistry calculations. Ultrafast intersystem crossing to the lowest triplet state (T₁) and effective non-radiative decay of this state to the ground state (S₀) are demonstrated in the neutral form 2SeU. However, the triplet lifetime of the deprotonated form 2SeU is found to be almost one order of magnitude longer than that in the neutral one. Quantum chemistry calculations indicate that the short triplet lifetime in 2SeU is due to excited state population decay through a crossing point between T₁ and S₀. In the deprotonated form, shortening the N1-C2 bond length makes the structural distortion more difficult and brings a larger energy barrier on the pathway to the T₁/S₀ crossing point, resulting in one order of magnitude increase of the triplet state lifetime. Our study reveals one key factor to regulate the triplet lifetime of 2SeU and sets the stage to further investigate the photophysical and photochemical properties of 2SeU-containing DNA/RNA duplexes.

Coherent vibrational modes promote the ultrafast internal conversion and intersystem crossing in thiobases

Thionated nucleobases are obtained by replacing oxygen with sulphur atoms in the canonical nucleobases. They absorb light efficiently in the near-ultraviolet, populating singlet states which undergo intersystem crossing to the...

Conical Intersection Passages of Molecules Probed by X-ray Diffraction and Stimulated Raman Spectroscopy

Conical intersections (CoIns) play an important role in ultrafast relaxation channels. Their monitoring remains a formidable experimental challenge. We theoretically compare the probing of the S₂ → S₁ CoIn passage in 4-thiouracil by monitoring its vibronic coherences, using off-resonant X-ray-stimulated Raman spectroscopy (TRUECARS) and time-resolved X-ray diffraction (TRXD). The quantum nuclear wavepacket (WP) dynamics provides an accurate picture of the photoinduced dynamics. Upon photoexcitation, the WP oscillates among the Franck-Condon point, the S₂ minimum, and the CoIn with a 70 fs period. A vibronic coherence first emerges at 20 fs and can be observed until the S₂ state is fully depopulated. The distribution of the vibronic frequencies involved in the coherence is recorded by the TRUECARS spectrogram. The TRXD signal provides spatial images of electron densities associated with the CoIn. In combination, the two signals provide a complementary picture of the nonadiabatic passage, which helps in the study of the underlying photophysics in thiobases.

Ab Initio Quasiclassical Simulation of Femtosecond Time-Resolved Two-Dimensional Electronic Spectra of Pyrazine

Two-dimensional (2D) electronic spectroscopy is a powerful nonlinear technique which provides spectroscopic information on two frequency axes as well as dynamical information as a function of the so-called waiting time. Herein, an ab initio theoretical framework for the simulation of electronic 2D spectra has been developed. The method is based on the classical approximation to the doorway-window representation of three-pulse photon-echo signals and the description of nuclear motion by classical trajectories. Nonadiabatic effects are taken into account by a trajectory surface-hopping algorithm. 2D electronic spectra were simulated with ab initio on-the-fly trajectory calculations using the ADC(2) electronic-structure method for the pyrazine molecule, which is a benchmark system for ultrafast radiationless decay through conical intersections. It is demonstrated that 2D spectroscopy with subfemtosecond UV pulses can provide unprecedented detailed information on the ultrafast photodynamics of polyatomic molecules.

In Silico Ultrafast Nonlinear Spectroscopy Meets Experiments: The Case of Perylene Bisimide Dye

Spectroscopy simulations are of paramount importance for the interpretation of experimental electronic spectra, the disentangling of overlapping spectral features, and the tracing of the microscopic origin of the observed signals. Linear and nonlinear simulations are based on the results drawn from electronic structure calculations that provide the necessary parameterization of the molecular systems probed by light. Here, we investigate the applicability of excited-state properties obtained from linear-response time-dependent density functional theory (TDDFT) in the description of nonlinear spectra by employing the pseudowavefunction approach and compare them with benchmarks from highly accurate RASSCF/RASPT2 calculations and with high temporal resolution experimental results. As a test case, we consider the prediction of femtosecond transient absorption and two-dimensional electronic spectroscopy of a perylene bisimide dye in solution. We find that experimental signals are well reproduced by both theoretical approaches, showing that the computationally cheaper TDDFT can be a suitable option for the simulation of nonlinear spectroscopy of molecular systems that are too large to be treated with higher-level RASSCF/RASPT2 methods.

Simulation and Analysis of the Transient Absorption Spectrum of 4-(N,N-Dimethylamino)benzonitrile (DMABN) in Acetonitrile

4-(N,N-Dimethylamino)benzonitrile (DMABN) is a well-known model compound for dual fluorescence-in sufficiently polar solvents, it exhibits two distinct fluorescence emission bands. The interpretation of its transient absorption (TA) spectrum in the visible range is the subject of a long-standing controversy. In the present study, we resolve this issue by calculating the TA spectrum on the basis of nonadiabatic molecular dynamics simulations. An unambiguous assignment of spectral signals to specific excited-state structures is achieved by breaking down the calculated spectrum into contributions from twisted and nontwisted molecular geometries. In particular, the much-discussed excited-state absorption band near 1.7 eV (ca. 700 nm) is attributed to the near-planar locally excited (LE) minimum on the S1 state. On the technical side, our study demonstrates that the second-order approximate coupled cluster singles and doubles (CC2) method can be used successfully to calculate the TA spectra of moderately large organic molecules, provided that the system in question does not approach a crossing between the lowest excited state and the singlet ground state within the time frame of the simulation.

Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations

The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.

Quantum and semiclassical dynamical studies of nonadiabatic processes in solution: achievements and perspectives

We concisely review the main methodological approaches to model nonadiabatic dynamics in isotropic solutions and their applications. Three general classes of models are identified as the most used to include solvent effects in the simulations. The first model describes the solvent as a set of harmonic collective modes coupled to the solute degrees of freedom, and the second as a continuum, while the third explicitly includes solvent molecules in the calculations. The issues related to the use of these models in semiclassical and quantum dynamical simulations are discussed, as well as the main limitations and perspectives of each approach.

iSPECTRON: A simulation interface for linear and nonlinear spectra with ab-initio quantum chemistry software

We introduce iSPECTRON, a program that parses data from common quantum chemistry software (NWChem, OpenMolcas, Gaussian, Cobramm, etc.), produces the input files for the simulation of linear and nonlinear spectroscopy of molecules with the Spectron code, and analyzes the spectra with a broad range of tools. Vibronic spectra are expressed in term of the electronic eigenstates, obtained from quantum chemistry computations, and vibrational/bath effects are incorporated in the framework of the displaced harmonic oscillator model, where all required quantities are computed at the Franck-Condon point. The program capabilities are illustrated by simulating linear absorption, transient absorption and two dimensional electronic spectra of the pyrene molecule. Calculations at two levels of electronic structure theory, time-dependent density functional theory (with NWChem) and RASSCF/RASPT2 (with OpenMolcas) are presented and compared where possible. The iSPECTRON program is available online at https://github.com/ispectrongit/iSPECTRON/ and distributed open source under the terms of the Educational Community License version 2.0 (ECL 2.0).

Ab Initio Surface-Hopping Simulation of Femtosecond Transient-Absorption Pump-Probe Signals of Nonadiabatic Excited-State Dynamics Using the Doorway-Window Representation

An ab initio theoretical framework for the simulation of femtosecond time-resolved transient absorption (TA) pump-probe (PP) spectra with quasi-classical trajectories is presented. The simulations are based on the classical approximation to the doorway-window (DW) representation of third-order four-wave-mixing signals. The DW formula accounts for the finite duration and spectral shape of the pump and probe pulses. In the classical DW formalism, classical trajectories are stochastically sampled from a positive definite doorway distribution, and the signals are evaluated by averaging over a positive definite window distribution. Nonadiabatic excited-state dynamics is described by a stochastic surface-hopping algorithm. The method has been implemented for the pyrazine molecule with the second-order algebraic-diagrammatic construction (ADC(2)) ab initio electronic-structure method. The methodology is illustrated by ab initio simulations of the ground-state bleach, stimulated emission, and excited-state absorption contributions to the TA PP spectrum of gas-phase pyrazine.

New Insights about the Photostability of DNA/RNA Bases: Triplet nπ* State Leads to Effective Intersystem Crossing in Pyrimidinones

The high photostability of DNA/RNA nucleobases is attributed to the effective internal conversions of their bright ¹ππ* states to the ground state through conical intersections. Intersystem crossing (ISC) from singlet to triplet excited states is a minor decay pathway in nucleobases and it is observed with ∼1-2% quantum yields (QYs) in pyrimidine bases. Presumably, ISC in pyrimidines takes place from the dark singlet ¹nπ* state to the lowest triplet ³ππ* state. However, recent studies showed that ISC from the initial populated bright ¹ππ* state to higher energy triplet ³nπ* states indeed occurs in the subpicosecond timescale. Such a mechanism is still poorly understood since direct observation of this pathway is challenging. Herein, excited state dynamics of three pyrimidinones, which share the same skeleton with pyrimidine bases, is investigated in different solvents. Compared to canonical pyrimidine bases, removing the oxygen atom at the C4 position revokes the low-lying dark ¹nπ* state in pyrimidinones, resulting in direct ISC from the S₁ (¹ππ*) state to triplet T₃ (³nπ*) state with much higher QYs. Meanwhile, hydrogen bonding between the carbonyl group in pyrimidinones and protic solvents can accelerate vibrational cooling of the hot S₁ (¹ππ*) state, leading to higher fluorescence QYs and smaller ISC rate constants. These results not only evidence the hypothesis of the direct ¹ππ* → ³nπ* ISC mechanism, but also contribute to a better understanding of triplet formation in pyrimidines.

Electronic Relaxation Pathways in Heavy-Atom-Free Photosensitizers Absorbing Near-Infrared Radiation and Exhibiting High Yields of Singlet Oxygen Generation

Heavy-atom-free photosensitizers (HAF-PSs) based on thionation of carbonyl groups of readily accessible organic compounds are rapidly emerging as a versatile class of molecules. However, their photochemical properties and electronic relaxation mechanisms are currently unknown. Investigating the excited-state dynamics is essential to understand their benefits and limitations and to develop photosensitizers with improved photochemical properties. Herein, the photochemical and electronic-structure properties of two of the most promising HAF-PSs developed to date are revealed. It is shown that excitation of thio-4-(dimethylamino)naphthalamide and thionated Nile Red with near-infrared radiation leads to the efficient population of the triplet manifold through multiple relaxation pathways in hundreds of femtoseconds. The strong singlet-triplet couplings in this family of photosensitizers should enable a broad range of applications, including in photodynamic therapy, photocatalysis, photovoltaics, organic LEDs, and photon up-conversion.

Investigating the mapping of chromophore excitations onto the electron detachment spectrum: photodissociation spectroscopy of iodide ion-thiouracil clusters

Laser photodissociation spectroscopy (3.1-5.7 eV) has been applied to iodide complexes of the non-native nucleobases, 2-thiouracil (2-TU), 4-thiouracil (4-TU) and 2,4-thiouracil (2,4-TU), to probe the excited states and intracluster electron transfer as a function of sulphur atom substitution. Photodepletion is strong for all clusters (I-·2-TU, I-·4-TU and I-·2,4-TU) and is dominated by electron detachment processes. For I-·4-TU and I-·2,4-TU, photodecay is accompanied by formation of the respective molecular anions, 4-TU- and 2,4-TU-, behaviour that is not found for other nucleobases. Notably, the I-·2TU complex does not fragment with formation of its molecular anion. We attribute the novel formation of 4-TU- and 2,4-TU- to the fact that these valence anions are significantly more stable than 2-TU-. We observe further similar behaviour for I-·4-TU and I-·2,4-TU relating to the general profile of their photodepletion spectra, since both strongly resemble the intrinsic absorption spectra of the respective uncomplexed thiouracil molecule. This indicates that the nucleobase chromophore excitations are determining the clusters' spectral profile. In contrast, the I-·2-TU photodepletion spectrum is dominated by the electron detachment profile, with the near-threshold dipole-bound excited state being the only distinct spectral feature. We discuss these observations in the context of differences in the dipole moments of the thionucleobases, and their impact on the coupling of nucleobase-centred transitions onto the electron detachment spectrum.

Quantum dynamics of the ππ*/nπ* decay of the epigenetic nucleobase 1,5-dimethyl-cytosine in the gas phase

We study the ultrafast dynamics of 1,5-dimethyl-cytosine, a model for 5-methyl-cytidine, after photoexcitation to the first two bright ππ* states, focusing on the possible population transfer to dark nπ* states. To that end we propagate the initial wave packets on the coupled potential energy surfaces of the seven lowest energy excited states modelled with a diabatic linear vibronic coupling (LVC) model, considering all the vibrational coordinates. Time-evolution is computed by the multilayer version of the multiconfigurational time dependent Hartree (ML-MCTDH) method. The LVC Hamiltonian is parametrized with time-dependent density functional theory (TD-DFT) calculations adopting PBE0 and CAM-B3LYP functionals, which provide a different energy gap between the lowest energy nπ* states and the spectroscopic ππ* state. Population of the lowest ππ* flows to a dark nπ* state which involves a lone pair (LP) of the carbonyl oxygen (nOπ*), but the extent of such transfer is much larger according to PBE0 than to CAM-B3LYP. Photoexcitation to the second bright state gives rise to much richer dynamics with an ultrafast (50 fs) complete decay to the lowest ππ*, to nOπ* and to another nπ* in which the excited electron comes from the LP of the ring nitrogen. We perform a detailed analysis of the vibronic dynamics both in terms of normal modes and valence coordinates (bond lengths and angles). The comparison with the analogous dynamics in 1-methyl-cytosine, a model for cytidine, provides insights into the effect of methylation at carbon 5 on the electronic and nuclear dynamics.

Role of Multistate Intersections in Photochemistry

It has been generally accepted that the intersection of potential energy surfaces can facilitate nonadiabatic transitions and plays a crucial role in photochemistry. Although most previous studies have focused on the conical intersection of two electronic states, multistate intersections are common in polyatomic molecules, and their key roles in photochemistry have been uncovered by electronic structure calculations and nonadiabatic dynamics simulations. In this Perspective, the algorithms for searching two- or three-state intersections are first examined with an emphasis on the latest development in a general algorithm for location of multistate intersections. Then, we focus on intersystem crossing (ISC) that occurs in the region of multistate intersection, paying more attention to how the state-specific spin-orbit coupling interaction influences nonadiabatic ISC processes. Finally, the interweaving of nonadiabatic dynamics simulation and electronic structure calculation has been recognized as a correct way to ascertain the vital roles of multistate intersections in photochemical reactions.

Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going

Computational spectroscopy is becoming a mandatory tool for the interpretation of the complex, and often congested, spectral maps delivered by modern non-linear multi-pulse techniques. The fields of Electronic Structure Methods, Non-Adiabatic Molecular Dynamics, and Theoretical Spectroscopy represent the three pillars of the virtual ultrafast optical spectrometer, able to deliver transient spectra in silico from first principles. A successful simulation strategy requires a synergistic approach that balances between the three fields, each one having its very own challenges and bottlenecks. The aim of this Perspective is to demonstrate that, despite these challenges, an impressive agreement between theory and experiment is achievable now regarding the modeling of ultrafast photoinduced processes in complex molecular architectures. Beyond that, some key recent developments in the three fields are presented that we believe will have major impacts on spectroscopic simulations in the very near future. Potential directions of development, pending challenges, and rising opportunities are illustrated.

Ultrafast Dynamics of the Two Lowest Bright Excited States of Cytosine and 1-Methylcytosine: A Quantum Dynamical Study

The nonadiabatic quantum dynamics (QD) of cytosine and 1-methylcytosine in the gas phase is simulated for 250 fs after a photoexcitation to one of the first two bright states. The nuclear wavepacket is propagated on the coupled diabatic potential energy surfaces of the lowest seven excited states, including ππ*, nπ*, and Rydberg states along all the vibrational degrees of freedom. We focus in particular on the interplay between the bright and the dark nπ* states, not considering the decay to the ground electronic state. To run these simulations, we implemented an automatic general procedure to parametrize linear vibronic coupling (LVC) models with time-dependent density functional theory (DFT) computations and interfaced it with Gaussian package. The wavepacket was propagated with the multilayer version of the multiconfigurational time dependent Hartree method. Two different density functionals, PBE0 and CAM-B3LYP, which provide a different description of the relative stability of the lowest energy dark states, were used to parametrize the LVC Hamiltonian. Part of the photoexcited population on lowest HOMO-LUMO transition (π_Hπ_L^*) decays within less than 100 fs to a nπ* state which mainly involves a promotion of an electron from the oxygen lone pair to the LUMO (n_Oπ_L^*). The population of the second ππ* state decays almost completely, in <100 fs, not only to π_Hπ_L^* and to n_Oπ_L^* states but also to another nπ_L^* state involving the nitrogen lone pair. The efficiency of the adopted protocol allowed us to check the accuracy of the predictions by repeating the QD simulations with different LVC Hamiltonians parametrized either at the ground-state minimum or at stationary structures of different relevant excited states.

Internal conversion and intersystem crossing dynamics of uracil upon double thionation: a time-resolved photoelectron spectroscopy study in the gas phase

The photophysical properties of 2,4-dithiouracil (2,4-DTU) in the gas phase are studied by time-resolved photoelectron spectroscopy (TRPES) with three different excitation wavelengths in direct extension of previous work on uracil (U), 2-thiouracil (2-TU) and 4-thiouracil (4-TU). Non-radiative deactivation in the canonical nucleobases like uracil mainly occurs via internal conversion (IC) along singlet excited states, although intersystem crossing (ISC) to a long-lived triplet state was confirmed to play a minor role. In thionated uracils, ISC to the triplet state becomes ultrafast and highly efficient with a quantum yield near unity; however, the lifetime of the triplet state is strongly dependent on the position of the sulfur atom. In 2-TU, ISC back to the ground state occurs within a few hundred picoseconds, whereas the population remains trapped in the lowest triplet state in the case of 4-TU. Upon doubling the degree of thionation, ISC remains highly efficient and dominates the photophysics of 2,4-DTU. However, several low-lying excited states contribute to competing IC and ISC pathways and a complex deactivation mechanism, which is evaluated here based on TRPES measurements and discussed in the context of the singly thionated uracils.

Detection of the thietane precursor in the UVA formation of the DNA 6-4 photoadduct

Notwithstanding the central biological role of the (6-4) photoadduct in the induction of skin cancer by sunlight, crucial mechanistic details about its formation have evaded characterization despite efforts spanning more than half a century. 4-Thiothymidine (4tT) has been widely used as an important model system to study its mechanism of formation, but the excited-state precursor, the intermediate species, and the time scale leading to the formation of the (6-4) photoadduct have remained elusive. Herein, steady-state and time-resolved spectroscopic techniques are combined with new and reported quantum-chemical calculations to demonstrate the excited state leading to the formation of the thietane intermediate, its rate, and the formation of the (6-4) photoadduct using the 5’-TT(4tT)T(4tT)TT-3’ DNA oligonucleotide. Efficient, sub-1 ps intersystem crossing leads to the population of a triplet minimum of the thietane intermediate in as short as 3 ps, which intersystem crosses to its ground state and rearranges to form the (6-4) photoadduct.

The mechanisms of formation of the (6-4) photoproducts in DNA damage by sunlight is still debated. Here the authors show, by optical spectroscopies and computations, the details of the formation of a (6-4) photoadduct via the thietane intermediate in a single-stranded DNA oligonucleotide.

Observation of Enhanced Dissociative Photochemistry in the Non-Native Nucleobase 2-Thiouracil

We present the first study to measure the dissociative photochemistry of 2-thiouracil (2-TU), an important nucleobase analogue with applications in molecular biology and pharmacology. Laser photodissociation spectroscopy is applied to the deprotonated and protonated forms of 2-TU, which are produced in the gas-phase using electrospray ionization mass spectrometry. Our results show that the deprotonated form of 2-thiouracil ([2-TU-H]⁻) decays predominantly by electron ejection and hence concomitant production of the [2-TU-H]· free-radical species, following photoexcitation across the UVA-UVC region. Thiocyanate (SCN⁻) and a m/z 93 fragment ion are also observed as photodecay products of [2-TU-H]⁻ but at very low intensities. Photoexcitation of protonated 2-thiouracil ([2-TU·H]⁺) across the same UVA-UVC spectral region produces the m/z 96 cationic fragment as the major photofragment. This ion corresponds to ejection of an HS· radical from the precursor ion and is determined to be a product of direct excited state decay. Fragment ions associated with decay of the hot ground state (i.e., the ions we would expect to observe if 2-thiouracil was behaving like UV-dissipating uracil) are observed as much more minor products. This behaviour is consistent with enhanced intersystem crossing to triplet excited states compared to internal conversion back to the ground state. These are the first experiments to probe the effect of protonation/deprotonation on thionucleobase photochemistry, and hence explore the effect of pH at a molecular level on their photophysical properties.

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

The excited-state relaxation mechanism of potential UVA-activated phototherapeutic molecules: trajectory surface hopping simulations of both 4-thiothymine and 2,4-dithiothymine

Recent experimental investigations of the photochemical properties of a series of sulfur-substituted pyrimidine derivatives provide insights into the phototherapeutic potential of these nucleobase variants. Herein we elucidate the triplet formation mechanism of two prospective UVA-activated phototherapeutic molecules, 4-thiothymine and 2,4-dithiothymine, upon photo-excitation by applying the trajectory surface hopping dynamics at the LR-TDDFT level. Our simulations reasonably reproduce the experimental time constants and demonstrate the preferred triplet formation pathway which starts from the S₁(n_Sπ*) state for both molecules. It is found that deactivation of the first bright state to the S₁(n_Sπ*) state proceeds through a mechanism involving elongation of the C5-C6 and C4-S8 bond-lengths and C2-pyramidalization in 4-thiothymine and involving elongation of the C5-C6 and C2-S7 bond-lengths in 2,4-dithiothymine. The intersystem crossing of 2,4-dithiothymine occurs either at geometries characterized by elongated C5-C6 and C2-S7 bond-lengths or at geometries showing elongated C5-C6 and C4-S8 bond-lengths as seen in 4-thiothymine. The solvents are found to affect the S₂ state decay of 4-thiothymine, leading to a competing pathway between S₂→ S₁ and S₂→ T₃. This study provides a molecular-level understanding of the underlying excited-state relaxation of the two UVA-activated thiopyrimidines, which may be linked to their potential applications in pharmacological science and also prove helpful for designing more effective phototherapeutic agents.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Excited State Lifetimes of Sulfur-Substituted DNA and RNA Monomers Probed Using the Femtosecond Fluorescence Up-Conversion Technique

Sulfur-substituted DNA and RNA nucleobase derivatives (a.k.a., thiobases) are an important family of biomolecules. They are used as prodrugs and as chemotherapeutic agents in medical settings, and as photocrosslinker molecules in structural-biology applications. Remarkably, excitation of thiobases with ultraviolet to near-visible light results in the population of long-lived and reactive triplet states on a time scale of hundreds of femtoseconds and with near-unity yields. This efficient nonradiative decay pathway explains the vanishingly small fluorescence yields reported for the thiobases and the scarcity of fluorescence lifetimes in the literature. In this study, we report fluorescence lifetimes for twelve thiobase derivatives, both in aqueous solution at physiological pH and in acetonitrile. Excitation is performed at 267 and 362 nm, while fluorescence emission is detected at 380, 425, 450, 525, or 532 nm. All the investigated thiobases reveal fluorescence lifetimes that decay in a few hundreds of femtoseconds and with magnitudes that depend and are sensitive to the position and degree of sulfur-atom substitution and on the solvent environment. Interestingly, however, three thiopyrimidine derivatives (i.e., 2-thiocytidine, 2-thiouridine, and 4-thiothymidine) also exhibit a small amplitude fluorescence component of a few picoseconds in aqueous solution. Furthermore, the N-glycosylation of thiobases to form DNA or RNA nucleoside analogues is demonstrated as affecting their fluorescence lifetimes. In aqueous solution, the fluorescence decay signals exciting at 267 nm are equal or slower than those collected exciting at 362 nm. In acetonitrile, however, the fluorescence decay signals recorded upon 267 nm excitation are, in all cases, faster than those measured exciting at 362 nm. A comparison to the literature values show that, while both the DNA and RNA nucleobase and thiobase derivatives exhibit sub-picosecond fluorescence lifetimes, the 1ππ* excited-state population in the nucleobase monomers primarily decay back to the ground state, whereas it predominantly populates long-lived and reactive triplet states in thiobase monomers.

Modeling multidimensional spectral lineshapes from first principles: application to water-solvated adenine

We theoretically describe spectral lineshape from first principles, providing insight into solvent–solute interactions in terms of static and dynamic disorder and how these shape experimental signals in linear and non-linear optical spectroscopies.

Quantum Chemical Modeling of the Photoinduced Activity of Multichromophoric Biosystems

Multichromophoric biosystems represent a broad family with very diverse members, ranging from light-harvesting pigment-protein complexes to nucleic acids. The former are designed to capture, harvest, efficiently transport, and transform energy from sunlight for photosynthesis, while the latter should dissipate the absorbed radiation as quickly as possible to prevent photodamages and corruption of the carried genetic information. Because of the unique electronic and structural characteristics, the modeling of their photoinduced activity is a real challenge. Numerous approaches have been devised building on the theoretical development achieved for single chromophores and on model Hamiltonians that capture the essential features of the system. Still, a question remains: is a general strategy for the accurate modeling of multichromophoric systems possible? By using a quantum chemical point of view, here we review the advancements developed so far highlighting differences and similarities with the single chromophore treatment. Finally, we outline the important limitations and challenges that still need to be tackled to reach a complete and accurate picture of their photoinduced properties and dynamics.