Sebastian Mai

A general method to describe intersystem crossing dynamics in trajectory surface hopping

Intersystem crossing is a radiationless process that can take place in a molecule irradiated by UV-Vis light, thereby playing an important role in many environmental, biological and technological processes. This paper reviews different methods to describe intersystem crossing dynamics, paying attention to semiclassical trajectory theories, which are especially interesting because they can be applied to large systems with many degrees of freedom. In particular, a general trajectory surface hopping methodology recently developed by the authors, which is able to include nonadiabatic and spin-orbit couplings in excited-state dynamics simulations, is explained in detail. This method, termed SHARC, can in principle include any arbitrary coupling, what makes it generally applicable to photophysical and photochemical problems, also those including explicit laser fields. A step-by-step derivation of the main equations of motion employed in surface hopping based on the fewest-switches method of Tully, adapted for the inclusion of spin-orbit interactions, is provided. Special emphasis is put on describing the different possible choices of the electronic bases in which spin-orbit can be included in surface hopping, highlighting the advantages and inconsistencies of the different approaches.

Ruthenium(II) photosensitizers of tridentate click-derived cyclometalating ligands: a joint experimental and computational study.

A systematic series of heteroleptic bis(tridentate)ruthenium(II) complexes of click-derived 1,3-bis(1,2,3-triazol-4-yl)benzene N^C^N-coordinating ligands was synthesized, analyzed by single crystal X-ray diffraction, investigated photophysically and electrochemically, and studied by computational methods. The presented comprehensive characterization allows a more detailed understanding of the radiationless deactivation mechanisms. Furthermore, we provide a fully optimized synthesis and systematic variations towards redox-matched, broadly and intensely absorbing, cyclometalated ruthenium(II) complexes. Most of them show a weak room-temperature emission and a prolonged excited-state lifetime. They display a broad absorption up to 700 nm and high molar extinction coefficients up to 20 000 M(-1)cm(-1) of the metal-to-ligand charge transfer bands, resulting in a black color. Thus, the complexes reveal great potential for dye-sensitized solar-cell applications.

Singlet and triplet excited-state dynamics study of the keto and enol tautomers of cytosine.

The photoinduced excited-state dynamics of the keto and enol forms of cytosine have been investigated by using ab initio surface-hopping to gain an understanding of the outcome of molecular beam femtosecond pump-probe photoionisation spectroscopy experiments. Both singlet and triplet states were included in the dynamics. The results show that triplet states play a significant role in the relaxation of the keto tautomer, whereas they are less important in the enol tautomer. In both forms, the T1 state minimum was found to be too low in energy to be detected in standard photoionisation spectroscopy experiments and therefore experimental decay times should arise from simultaneous relaxation to the ground state and additional intersystem crossing followed by internal conversion to the T1 state. In agreement with available experimental lifetimes, we observed three decay constants of 7, 270 and 1900 fs, the first two coming from the keto tautomer and the third from the enol form. Deactivation of the enol tautomer is due to internal conversion to the ground state through two ethylenic-type S1/S0 conical intersections.

Intersystem Crossing Pathways in the Noncanonical Nucleobase 2-Thiouracil: A Time-Dependent Picture

The deactivation mechanism after ultraviolet irradiation of 2-thiouracil has been investigated using nonadiabatic dynamics simulations at the MS-CASPT2 level of theory. It is found that after excitation the S2 quickly relaxes to S1, and from there intersystem crossing takes place to both T2 and T1 with a time constant of 400 fs and a triplet yield above 80%, in very good agreement with recent femtosecond experiments in solution. Both indirect S1 → T2 → T1 and direct S1 → T1 pathways contribute to intersystem crossing, with the former being predominant. The results contribute to the understanding of how some noncanonical nucleobases respond to harmful ultraviolet light, which could be relevant for prospective photochemotherapeutic applications.

The origin of efficient triplet state population in sulfur-substituted nucleobases

Elucidating the photophysical mechanisms in sulfur-substituted nucleobases (thiobases) is essential for designing prospective drugs for photo- and chemotherapeutic applications. Although it has long been established that the phototherapeutic activity of thiobases is intimately linked to efficient intersystem crossing into reactive triplet states, the molecular factors underlying this efficiency are poorly understood. Herein we combine femtosecond transient absorption experiments with quantum chemistry and nonadiabatic dynamics simulations to investigate 2-thiocytosine as a necessary step to unravel the electronic and structural elements that lead to ultrafast and near-unity triplet-state population in thiobases in general. We show that different parts of the potential energy surfaces are stabilized to different extents via thionation, quenching the intrinsic photostability of canonical DNA and RNA nucleobases. These findings satisfactorily explain why thiobases exhibit the fastest intersystem crossing lifetimes measured to date among bio-organic molecules and have near-unity triplet yields, whereas the triplet yields of canonical nucleobases are nearly zero.

A Static Picture of the Relaxation and Intersystem Crossing Mechanisms of Photoexcited 2-Thiouracil.

Accurate excited-state quantum chemical calculations on 2-thiouracil, employing large active spaces and up to quadruple-ζ quality basis sets in multistate complete active space perturbation theory calculations, are reported. The results suggest that the main relaxation path for 2-thiouracil after photoexcitation should be S2 → S1 → T2 → T1, and that this relaxation occurs on a subpicosecond time scale. There are two deactivation pathways from the initially excited bright S2 state to S1, one of which is nearly barrierless and should promote ultrafast internal conversion. After relaxation to the S1 minimum, small singlet–triplet energy gaps and spin–orbit couplings of about 130 cm–1 are expected to facilitate intersystem crossing to T2, from where very fast internal conversion to T1 occurs. An important finding is that 2-thiouracil shows strong pyramidalization at the carbon atom of the thiocarbonyl group in several excited states.

Electronic and Structural Elements That Regulate the Excited-State Dynamics in Purine Nucleobase Derivatives

The excited-state dynamics of the purine free base and 9-methylpurine are investigated using experimental and theoretical methods. Femtosecond broadband transient absorption experiments reveal that excitation of these purine derivatives in aqueous solution at 266 nm results primarily in ultrafast conversion of the S2(ππ*) state to the vibrationally excited 1nπ* state. Following vibrational and conformational relaxation, the 1nπ* state acts as a doorway state in the efficient population of the triplet manifold with an intersystem crossing lifetime of hundreds of picoseconds. Experiments show an almost 2-fold increase in the intersystem crossing rate on going from polar aprotic to nonpolar solvents, suggesting that a solvent-dependent energy barrier must be surmounted to access the singlet-to-triplet crossing region. Ab initio static and surface-hopping dynamics simulations lend strong support to the proposed relaxation mechanism. Collectively, the experimental and computational results demonstrate that the accessibility of the nπ* states and the topology of the potential energy surfaces in the vicinity of conical intersections are key elements in controlling the excited-state dynamics of the purine derivatives. From a structural perspective, it is shown that the purine chromophore is not responsible for the ultrafast internal conversion in the adenine and guanine monomers. Instead, C6 functionalization plays an important role in regulating the rates of radiative and nonradiative relaxation. C6 functionalization inhibits access to the 1nπ* state while simultaneously facilitating access to the 1ππ*(La)/S0 conical intersection, such that population of the 1nπ* state cannot compete with the relaxation pathways to the ground state involving ring puckering at the C2 position.

Efficient and Flexible Computation of Many-Electron Wave Function Overlaps

A new algorithm for the computation of the overlap between many-electron wave functions is described. This algorithm allows for the extensive use of recurring intermediates and thus provides high computational efficiency. Because of the general formalism employed, overlaps can be computed for varying wave function types, molecular orbitals, basis sets, and molecular geometries. This paves the way for efficiently computing nonadiabatic interaction terms for dynamics simulations. In addition, other application areas can be envisaged, such as the comparison of wave functions constructed at different levels of theory. Aside from explaining the algorithm and evaluating the performance, a detailed analysis of the numerical stability of wave function overlaps is carried out, and strategies for overcoming potential severe pitfalls due to displaced atoms and truncated wave functions are presented.

Benzophenone Ultrafast Triplet Population : Revisiting the Kinetic Model by Surface-Hopping Dynamics

The photochemistry of benzophenone, a paradigmatic organic molecule for photosensitization, was investigated by means of surface-hopping ab initio molecular dynamics. Different mechanisms were found to be relevant within the first 600 fs after excitation; the long-debated direct (S1 → T1) and indirect (S1 → T2 → T1) mechanisms for population of the low-lying triplet state are both possible, with the latter being prevalent. Moreover, we established the existence of a kinetic equilibrium between the two triplet states, never observed before. This fact implies that a significant fraction of the overall population resides in T2, eventually allowing one to revisit the usual spectroscopic assignment proposed by transient absorption spectroscopy. This finding is of particular interest for photocatalysis as well as for DNA damages studies because both T1 and T2 channels are, in principle, available for benzophenone-mediated photoinduced energy transfer toward DNA.

Excitation of Nucleobases from a Computational Perspective II: Dynamics

This chapter is devoted to unravel the relaxation processes taking place after photoexcitation of isolated DNA/RNA nucleobases in gas phase from a time-dependent perspective. To this aim, several methods are at hand, ranging from full quantum dynamics to various flavours of semiclassical or ab initio molecular dynamics, each with its advantages and its limitations. As this contribution shows, the most common approach employed up to date to learn about the deactivation of nucleobases in gas phase is a combination of the Tully surface hopping algorithm with on-the-fly CASSCF calculations. Different dynamics methods or, even more dramatically, different electronic structure methods can provide different dynamics. A comprehensive review of the different mechanisms suggested for each nucleobase is provided and compared to available experimental time scales. The results are discussed in a general context involving the effects of the different applied electronic structure and dynamics methods. Mechanistic similarities and differences between the two groups of nucleobases - the purine derivatives (adenine and guanine) and the pyrimidine derivatives (thymine, uracil, and cytosine) - are elucidated. Finally, a perspective on the future of dynamics simulations in the context of nucleobase relaxation is given.