Philipp Marquetand

SHARC: ab Initio Molecular Dynamics with Surface Hopping in the Adiabatic Representation Including Arbitrary Couplings.

We present a semiclassical surface-hopping method which is able to treat arbitrary couplings in molecular systems including all degrees of freedom. A reformulation of the standard surface-hopping scheme in terms of a unitary transformation matrix allows for the description of interactions like spin-orbit coupling or transitions induced by laser fields. The accuracy of our method is demonstrated in two systems. The first one, consisting of two model electronic states, validates the semiclassical approach in the presence of an electric field. In the second one, the dynamics in the IBr molecule in the presence of spin-orbit coupling after laser excitation is investigated. Due to an avoided crossing that originates from spin-orbit coupling, IBr dissociates into two channels: I + Br((2)P3/2) and I + Br*((2)P1/2). In both systems, the obtained results are in very good agreement with those calculated from exact quantum dynamical simulations.

Femtosecond Intersystem Crossing in the DNA Nucleobase Cytosine.

Ab initio molecular dynamics including nonadiabatic and spin-orbit couplings on equal footing is used to unravel the deactivation of cytosine after UV light absorption. Intersystem crossing (ISC) is found to compete directly with internal conversion in tens of femtoseconds, thus making cytosine the organic compound with the fastest triplet population calculated so far. It is found that close degeneracy between singlet and triplet states can more than compensate for very small spin-orbit couplings, leading to efficient ISC. The femtosecond nature of the ISC process highlights its importance in photochemistry and challenges the conventional view that large singlet-triplet couplings are required for an efficient population flow into triplet states. These findings are important to understand DNA photostability and the photochemistry and dynamics of organic molecules in general.

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.

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.

Communication: GAIMS—Generalized Ab Initio Multiple Spawning for both internal conversion and intersystem crossing processes

Full multiple spawning is a formally exact method to describe the excited-state dynamics of molecular systems beyond the Born-Oppenheimer approximation. However, it has been limited until now to the description of radiationless transitions taking place between electronic states with the same spin multiplicity. This Communication presents a generalization of the full and ab initio multiple spawning methods to both internal conversion (mediated by nonadiabatic coupling terms) and intersystem crossing events (triggered by spin-orbit coupling matrix elements) based on a spin-diabatic representation. The results of two numerical applications, a model system and the deactivation of thioformaldehyde, validate the presented formalism and its implementation.

Nonadiabatic ab initio molecular dynamics including spin–orbit coupling and laser fields

Nonadiabatic ab initio molecular dynamics (MD) including spin-orbit coupling (SOC) and laser fields is investigated as a general tool for studies of excited-state processes. Up to now, SOCs are not included in standard ab initio MD packages. Therefore, transitions to triplet states cannot be treated in a straightforward way. Nevertheless, triplet states play an important role in a large variety of systems and can now be treated within the given framework. The laser interaction is treated on a non-perturbative level that allows nonlinear effects like strong Stark shifts to be considered. As MD allows for the handling of many atoms, the interplay between triplet and singlet states of large molecular systems will be accessible. In order to test the method, IBr is taken as a model system, where SOC plays a crucial role for the shape of the potential curves and thus the dynamics. Moreover, the influence of the nonresonant dynamic Stark effect is considered. The latter is capable of controlling reaction barriers by electric fields in time-reversible conditions, and thus a control laser using this effect acts like a photonic catalyst. In the IBr molecule, the branching ratio at an avoided crossing, which arises from SOC, can be influenced.