The on-the-fly surface-hopping program system Newton-X: Application to ab initio simulation of the nonadiabatic photodynamics of benchmark systems

Intersystem Crossing and Triplet Dynamics in an Iron(II) N-Heterocyclic Carbene Photosensitizer

The electronic excited states of the iron(II) complex [Fe^II(tpy)(pyz-NHC)]²⁺ [tpy = 2,2′:6′,2″-terpyridine; pyz-NHC = 1,1′-bis(2,6-diisopropylphenyl)pyrazinyldiimidazolium-2,2′-diylidene] and their relaxation pathways have been theoretically investigated. To this purpose, trajectory surface-hopping simulations within a linear vibronic coupling model including a 244-dimensional potential energy surface (PES) with 20 singlet and 20 triplet coupled states have been used. The simulations show that, after excitation to the lowest-energy absorption band of predominant metal-to-ligand charge-transfer character involving the tpy ligand, almost 80% of the population undergoes intersystem crossing to the triplet manifold in about 50 fs, while the remaining 20% decays through internal conversion to the electronic ground state in about 300 fs. The population transferred to the triplet states is found to deactivate into two different regions of the PESs, one where the static dipole moment is small and shows increased metal-centered character and another with a large static dipole moment, where the electron density is transferred from the tpy to pyz-NHC ligand. Coherent oscillations of 400 fs are observed between these two sets of triplet populations, until the mixture equilibrates to a ratio of 60:40. Finally, the importance of selecting suitable normal modes is highlighted—a choice that can be far from straightforward in transition-metal complexes with hundreds of degrees of freedom.

Trajectory surface-hopping simulations with a linear vibronic coupling model reveal the competition of major intersystem crossing versus minor internal conversion dynamics in an iron(II) N-heterocyclic carbene dye. The triplet population bifurcates into two regions of the potential energy surfaces, characterized by small and large static dipole moments due to different electronic character and showing coherent oscillations of 400 fs until both triplet populations coexist in a mixture of 60:40.

SSAIMS-Stochastic-Selection Ab Initio Multiple Spawning for Efficient Nonadiabatic Molecular Dynamics

Ab initio multiple spawning provides a powerful and accurate way of describing the excited-state dynamics of molecular systems, whose strength resides in the proper description of coherence effects during nonadiabatic processes thanks to the coupling of trajectory basis functions. However, the simultaneous propagation of a large number of trajectory basis functions can be numerically inconvenient. We propose here an elegant and simple solution to this issue, which consists of (i) detecting uncoupled groups of coupled trajectory basis functions and (ii) selecting stochastically one of these groups to continue the ab initio multiple spawning dynamics. We show that this procedure can reproduce the results of full ab initio multiple spawning dynamics in cases where the uncoupled groups of trajectory basis functions stay uncoupled throughout the dynamics (which is often the case in high-dimensional problems). We present and discuss the aforementioned idea in detail and provide simple numerical applications on indole, ethylene, and protonated formaldimine, highlighting the potential of stochastic-selection ab initio multiple spawning.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.

Identification of a receiver triplet state in the ultrafast intersystem crossing of carbonylpyrenes

The near-unity triplet quantum yield of photoexcited carbonyl functionalized pyrenes is theoretically investigated. The estimated energetics of singlet-triplet manifolds and relevant spin-orbit coupling parameters strongly suggest triplet state formation via the S₁→ T₄/T₅ pathway. Quantum wavepacket dynamics of triplet manifolds within the linear vibronic coupling approach reveal that the receiver triplet state would undergo rapid internal conversion decay to the lower triplet state(s), facilitating efficient triplet generation by minimizing the reverse intersystem crossing possibilities. On the basis of these results, a unified mechanism is proposed to describe the ultrafast intersystem crossing process in these molecules.

Competing ultrafast photoinduced electron transfer and intersystem crossing of [Re(CO) 3 (Dmp)(His124)(Trp122)] + in Pseudomonas aeruginosa azurin: a nonadiabatic dynamics study

We present a computational study of sub-picosecond nonadiabatic dynamics in a rhenium complex coupled electronically to a tryptophan (Trp) side chain of Pseudomonas aeruginosa azurin, a prototypical protein used in the study of electron transfer in proteins. To gain a comprehensive understanding of the photoinduced processes in this system, we have carried out vertical excitation calculations at the TDDFT level of theory as well as nonadiabatic dynamics simulations using the surface hopping including arbitrary couplings (SHARC) method coupled to potential energy surfaces represented with a linear vibronic coupling model. The results show that the initial photoexcitation populates both singlet metal-to-ligand charge transfer (MLCT) and singlet charge-separated (CS) states, where in the latter an electron was transferred from the Trp amino acid to the complex. Subsequently, a complex mechanism of simultaneous intersystem crossing and electron transfer leads to the sub-picosecond population of triplet MLCT and triplet CS states. These results confirm the assignment of the sub-ps time constants of previous experimental studies and constitute the first computational evidence for the ultrafast formation of the charge-separated states in Re-sensitized azurin.

Electronic Excited States and UV-Vis Absorption Spectra of the Dihydropyrene/Cyclophanediene Photochromic Couple: a Theoretical Investigation

Dihydropyrene (DHP)/cyclophanediene (CPD) is a fascinating photoswitchable organic system displaying negative photochromism. Upon irradiation in the visible region, the colored DHP can be converted to its open-ring CPD colorless isomer, which can be converted back to DHP by UV light. DHP and CPD thus possess very different absorption spectra whose absorption bands have never been assigned in detail so far. In this work, we characterize the vertical electronic transitions of the first six and seven excited states of DHP and CPD, respectively, aiming for a realistic comparison with experiment. We used state-of-the-art electronic structure methods [e.g., complete active space second-order perturbation theory (CASPT2), n-electron valence-state perturbation theory (NEVPT2), extended multiconfigurational quasi-degenerate perturbation theory (XMCQDPT2), and third-order algebraic diagrammatic construction ADC(3)] capable of describing differential electron correlation. Vertical transition energies were also computed with time-dependent density functional theory (TD-DFT) and compared to these accurate methods. After the reliability of TD-DFT was validated for the main optical transitions, this efficient method was used to simulate the absorption spectra of DHP and CPD in the framework of the Franck-Condon Herzberg-Teller approximation and also using the nuclear ensemble approach. Overall, for both methods, the simulated absorption spectra reproduce nicely the main spectral features of the DHP and CPD isomers, that is, the main four absorption bands of increasing intensity of DHP and the absorption rise below 300 nm for CPD.

Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications

This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.

Polarizable embedding QM/MM: the future gold standard for complex (bio)systems?

We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.

Nonadiabatic Molecular Dynamics on Graphics Processing Units: Performance and Application to Rotary Molecular Motors

Nonadiabatic molecular dynamics (NAMD) simulations of molecular systems require the efficient evaluation of excited-state properties, such as energies, gradients, and nonadiabatic coupling vectors. Here, we investigate the use of graphics processing units (GPUs) in addition to central processing units (CPUs) to efficiently calculate these properties at the time-dependent density functional theory (TDDFT) level of theory. Our implementation in the FermiONs++ program package uses the J-engine and a preselective screening procedure for the calculation of Coulomb and exchange kernels, respectively. We observe good speed-ups for small and large molecular systems (comparable to those observed in ground-state calculations) and reduced (down to sublinear) scaling behavior with respect to the system size (depending on the spatial locality of the investigated excitation). As a first illustrative application, we present efficient NAMD simulations of a series of newly designed light-driven rotary molecular motors and compare their S1 lifetimes. Although all four rotors show different S1 excitation energies, their ability to rotate upon excitation is conserved, making the series an interesting starting point for rotary molecular motors with tunable excitation energies.

Synthesis, Electronic Spectroscopy, and Photochemistry of Methacrolein Oxide: A Four-Carbon Unsaturated Criegee Intermediate from Isoprene Ozonolysis

Ozonolysis of isoprene, one of the most abundant volatile organic compounds in the earth's atmosphere, generates the four-carbon unsaturated methacrolein oxide (MACR-oxide) Criegee intermediate. The first laboratory synthesis and direct detection of MACR-oxide is achieved through reaction of photolytically generated, resonance-stabilized iodoalkene radicals with oxygen. MACR-oxide is characterized on its first π* ← π electronic transition using a ground-state depletion method. MACR-oxide exhibits a broad UV-visible spectrum peaked at 380 nm with weak oscillatory structure at long wavelengths ascribed to vibrational resonances. Complementary theory predicts two strong π* ← π transitions arising from extended conjugation across MACR-oxide with overlapping contributions from its four conformers. Electronic promotion to the 1¹ππ* state agrees well with experiment, and results in nonadiabatic coupling and prompt release of O ¹D products observed as anisotropic velocity-map images. This UV-visible detection scheme will enable study of its unimolecular and bimolecular reactions under thermal conditions of relevance to the atmosphere.

Nonadiabatic Dynamics Simulation Predict Intersystem Crossing in Nitroaromatic Molecules on a Picosecond Time Scale

Previous time-resolved spectroscopic experiments and static quantum-chemical calculations attributed nitronaphthalene derivatives one of the fastest time scales for intersystem crossing within organic molecules, reaching the 100 fs mark. Nonadiabatic dynamics simulations on three nitronaphthalene derivatives challenge this view, showing that the experimentally observed ∼100 fs process corresponds to internal conversion in the singlet manifolds. Intersystem crossing, instead, takes place on a longer time scale of ∼1 ps. The dynamics simulations further reveal that the spin transitions occur via two distinct pathways with different contribution for the three systems, which are determined by electronic factors and the torsion of the nitro group. This study, therefore, indicates that the existence of sub-picosecond intersystem crossing in other nitroaromatic molecules should be questioned.

Excited state dynamics of CH2I2 and CH2BrI studied with UV pump VUV probe photoelectron spectroscopy

We compare the excited state dynamics of diiodomethane (CH₂I₂) and bromoiodomethane (CH₂BrI) using time resolved photoelectron spectroscopy. A 4.65 eV UV pump pulse launches a dissociative wave packet on excited states of both molecules and the ensuing dynamics are probed via photoionization using a 7.75 eV probe pulse. The resulting photoelectrons are measured with the velocity map imaging technique for each pump-probe delay. Our measurements highlight differences in the dynamics for the two molecules, which are interpreted with high-level ab initio molecular dynamics (trajectory surface hopping) calculations. Our analysis allows us to associate features in the photoelectron spectrum with different portions of the excited state wave packet represented by different trajectories. The excited state dynamics in bromoiodomethane are simple and can be described in terms of direct dissociation along the C-I coordinate, whereas the dynamics in diiodomethane involve internal conversion and motion along multiple dimensions.

Performance of trajectory surface hopping method in the treatment of ultrafast intersystem crossing dynamics

We carried out extensive studies to examine the performance of the fewest-switches surface hopping method in the description of the ultrafast intersystem crossing dynamic of various singlet-triplet (S-T) models by comparison with the results of the exact full quantum dynamics. Different implementation details and some derivative approaches were examined. As expected, it is better to perform the trajectory surface hopping calculations in the spin-adiabatic representation or by the local diabatization approach, instead of in the spin-diabatic representation. The surface hopping method provides reasonable results for the short-time dynamics in the S-T model with weak spin-orbital coupling (diabatic coupling), although it does not perform well in the models with strong spin-orbital coupling (diabatic coupling). When the system accesses the S-T potential energy crossing with rather high kinetic energy, the trajectory surface hopping method tends to produce a good description of the nonadiabatic intersystem crossing dynamics. The impact of the decoherence correction on the performance of the trajectory surface hopping is system dependent. It improves the result accuracy in many cases, while its influence may also be minor for other cases.

Computational insight into excited states of the ring-opening radicals from the pyrolysis of furan biofuels

The low-lying valence excited states and Rydberg states of the radical species from the ring-opening reactions in pyrolysis of furan biofuels have been determined by extensive density functional theory and sophisticated wave function theory calculations. The radicals 1-C₄ H₅ O-2, 2-furylCH₂ , and 4-C₆ H₇ O with the delocalized π-type single electron are predicted to be most stable among the reactive species here for furan, 2-methyfuran, and 2,5-dimethylfuran, respectively. Predicted vertical transition energies by TD-CAM-B3LYP show good agreement with those by CASPT2. Some among the electronic excitations to low-lying states can take place in the visible light region, and they may be involved in the combustion process. Further surface hopping dynamics simulations on the excited states of the most stable ring-opening radical 1-C₄ H₅ O-2 of furan as an example reveal that 89.9% sampling trajectories at the initial excited state of 2² A"(π¹ π*² ) decay to the 1² A'(n¹ π*² ) state within an average of 384 fs, and then 81.2% trajectories at the 1² A' state go to the ground state within an average of 114 fs. At the end of the simulation for 1000 fs, 18.8% trajectories still stay on the excited states of 2² A" and 1² A', suggesting that the reactive radicals in the ground state are mainly responsible for the combustion chemistry of furan biofuels. © 2018 Wiley Periodicals, Inc.

Asymmetric substitution changes the UV-induced nonradiative decay pathway and the spectra behaviors of β-diketones

Asymmetric substitution has not been termed as an essential factor in studying photo-induced ultrafast dynamics of molecular system. Asymmetric 4-hydroxybut-3-en-2-one (HEO), together with symmetric malonaldehyde (MA) and acetylacetone (AA), have been provided as target sample to study the nonradiative decay (ND) processes of β-diketones. An effective ND pathway of the three molecules is presented that their excited second (S₂) states transfer to first (S₁) state by nonadiabatic surface hopping, and then transfer to triplet (T₁) state by crossing minimum energy crossing point (MECP), after which decay to ground (S₀) state through MECP. More importantly, the asymmetric substitution of HEO induces the proton transfer in the S₁ state and generates a proton-transferred conformer with lowest energy, which does not occur for MA and AA. This change exploits a new ND pathway that the S₁ state decays to the proton transferred T₁ state and then undergoes reverse proton transfer to S₀ state through the MECPs between the three states. The two pathways of HEO give detailed energy and geometric information on surface hopping of S₂/S₁ and MECPs of S₁/T₁/S₀, and interpret the reason of the ND pathway while not spectra emission. This result is significantly different from the previous reported ND pathway of photoisomerization or conical intersection between different states. This work shows that asymmetric substitution changes the molecular structure and then changes their spectra behaviors.

Electronic spectroscopy of methyl vinyl ketone oxide: A four-carbon unsaturated Criegee intermediate from isoprene ozonolysis

Ozonolysis of isoprene, one of the most abundant volatile organic compounds in the atmosphere, proceeds through methyl vinyl ketone oxide (MVK-oxide), methacrolein oxide, and formaldehyde oxide (CH₂OO) Criegee intermediates. The present study focuses on MVK-oxide, a four-carbon unsaturated carbonyl oxide intermediate, using vacuum ultraviolet photoionization at 118 nm and UV-visible induced depletion of the m/z = 86 mass channel to characterize its first π^* ← π electronic transition. The electronic spectrum is broad and unstructured with its peak at 388 nm (3.2 eV). The MVK-oxide spectrum is shifted to a significantly longer wavelength than CH₂OO and alkyl-substituted Criegee intermediates studied previously due to extended conjugation across the vinyl and carbonyl oxide groups. Electronic excitation results in rapid dissociation at λ ≤ 430 nm to methyl vinyl ketone and O ¹D products, the latter detected by 2 + 1 resonance enhanced multiphoton ionization using velocity map imaging. Complementary electronic structure calculations (CASPT2(12,10)/AVDZ) predict two π^* ← π transitions with significant oscillator strength for each of the four conformers of MVK-oxide with vertical excitation energies (and corresponding wavelengths) in the 3.1-3.6 eV (350-400 nm) and 4.5-5.5 eV (220-280 nm) regions. The computed electronic absorption profile of MVK-oxide, based on a Wigner distribution of ground state configurations and summed over the four conformers, is predicted to peak at 397 nm. UV-visible spectroscopy on the first π^* ← π transition is shown by a combination of experiment and theory to provide a sensitive method for detection of the MVK-oxide Criegee intermediate that will enable further studies of its photochemistry and unimolecular and bimolecular reaction dynamics.

An efficient approximate algorithm for nonadiabatic molecular dynamics

We propose a modification to the nonadiabatic surface hopping calculation method formulated in a paper by Yu et al. [Phys. Chem. Chem. Phys. 16, 25883 (2014)], which is a multidimensional extension of the Zhu-Nakamura theory with a practical diabatic gradient estimation algorithm. In our modification, their diabatic gradient estimation algorithm, which is based on a simple interpolation of the adiabatic potential energy surfaces, is replaced by an algorithm using the numerical derivatives of the adiabatic gradients. We then apply the algorithm to several models of nonadiabatic dynamics, both analytic and ab initio models, to numerically demonstrate that our method indeed widens the applicability and robustness of their method. We also discuss the validity and limitations of our new nonadiabatic surface hopping method while considering in mind potential applications to excited-state dynamics of biomolecules or unconventional nonadiabatic dynamics such as radiation decay processes in ultraintense X-ray fields.

Effects of π-expansion, an additional hydroxyl group, and substitution on the excited state single and double proton transfer of 2-hydroxybenzaldehyde and its relative compounds: TD-DFT static and dynamic study

The effects of π-expansion, an extra hydroxyl group, and substituents on the photophysical properties, the excited state single proton transfer and the double proton transfer of 2-hydroxybenzaldehyde and its relatives have been theoretically investigated using TD-DFT.

Ab initio quantum-chemical computations of the absorption cross sections of HgX2 and HgXY (X, Y = Cl, Br, and I): molecules of interest in the Earth's atmosphere

The electronic-structure properties of the low-lying electronic states and the absorption cross sections of mercury halides have been determined within the UV-vis spectrum range (170 nm ≤ λ_photon ≤ 600 nm).

Finite-temperature Wigner phase-space sampling and temperature effects on the excited-state dynamics of 2-nitronaphthalene

The concept of finite temperature Wigner phase-space sampling allowing the population of vibrationally excited states is introduced and employed to study temperature effects on the absorption spectrum of 2-nitronaphtalene (2NN) and its relaxation dynamics.

Ab Initio Prediction of Fluorescence Lifetimes Involving Solvent Environments by Means of COSMO and Vibrational Broadening

The fluorescence lifetime is a key property of fluorophores that can be utilized for microenvironment probing, analyte sensing, and multiplexing as well as barcoding applications. For the rational design of lifetime probes and barcodes, theoretical methods have been developed to enable the ab initio prediction of this parameter, which depends strongly on interactions with solvent molecules and other chemical species in the emitteŕs immediate environment. In this work, we investigate how a conductor-like screening model (COSMO) can account for variations in fluorescence lifetimes that are caused by such fluorophore-solvent interactions. Therefore, we calculate vibrationally broadened fluorescence spectra using the nuclear ensemble method to obtain distorted molecular geometries to sample the electronic transitions with time-dependent density functional theory (TDDFT). The influence of the solvent on fluorescence lifetimes is accounted for with COSMO. For example, for 4-hydroxythiazole fluorophore containing different heteroatoms and acidic and basic moieties in aprotic and protic solvents of varying polarity, this approach was compared to experimentally determined lifetimes in the same solvents. Our results demonstrate a good correlation between theoretically predicted and experimentally measured fluorescence lifetimes except for the polar solvents ethanol and acetonitrile that can specifically interact with the heteroatoms and the carboxylic acid of the thiazole derivative.

Photoreduction of gaseous oxidized mercury changes global atmospheric mercury speciation, transport and deposition

Anthropogenic mercury (Hg(0)) emissions oxidize to gaseous Hg(II) compounds, before deposition to Earth surface ecosystems. Atmospheric reduction of Hg(II) competes with deposition, thereby modifying the magnitude and pattern of Hg deposition. Global Hg models have postulated that Hg(II) reduction in the atmosphere occurs through aqueous-phase photoreduction that may take place in clouds. Here we report that experimental rainfall Hg(II) photoreduction rates are much slower than modelled rates. We compute absorption cross sections of Hg(II) compounds and show that fast gas-phase Hg(II) photolysis can dominate atmospheric mercury reduction and lead to a substantial increase in the modelled, global atmospheric Hg lifetime by a factor two. Models with Hg(II) photolysis show enhanced Hg(0) deposition to land, which may prolong recovery of aquatic ecosystems long after Hg emissions are lowered, due to the longer residence time of Hg in soils compared with the ocean. Fast Hg(II) photolysis substantially changes atmospheric Hg dynamics and requires further assessment at regional and local scales.

Reduction of gaseous Hg(II) compounds drives atmospheric mercury wet and dry deposition to Earth surface ecosystems. Global Hg models assume this reduction takes place in clouds. Here the authors report a new gas-phase Hg photochemical mechanism that changes atmospheric mercury lifetime and its deposition to the surface.

Simulated and Experimental Time-Resolved Photoelectron Spectra of the Intersystem Crossing Dynamics in 2-Thiouracil

We report time-dependent photoelectron spectra recorded with a single-photon ionization setup and extensive simulations of the same spectra for the excited-state dynamics of 2-thiouracil (2TU) in the gas phase. We find that single-photon ionization produces very similar results as two-photon ionization, showing that the probe process does not have a strong influence on the measured dynamics. The good agreement between the single-photon ionization experiments and the simulations shows that the norms of Dyson orbitals allow for qualitatively describing the ionization probabilities of 2TU. This reasonable performance of Dyson norms is attributed to the particular electronic structure of 2TU, where all important neutral and ionic states involve similar orbital transitions and thus the shape of the Dyson orbitals do not strongly depend on the initial neutral and final ionic state. We argue that similar situations should also occur in other biologically relevant thio-nucleobases, and that the time-resolved photoelectron spectra of these bases could therefore be adequately modeled with the techniques employed here.

Nonadiabatic dynamics: The SHARC approach

We review the Surface Hopping including ARbitrary Couplings (SHARC) approach for excited-state nonadiabatic dynamics simulations. As a generalization of the popular surface hopping method, SHARC allows simulating the full-dimensional dynamics of molecules including any type of coupling terms beyond nonadiabatic couplings. Examples of these arbitrary couplings include spin-orbit couplings or dipole moment-laser field couplings, such that SHARC can describe ultrafast internal conversion, intersystem crossing, and radiative processes. The key step of the SHARC approach consists of a diagonalization of the Hamiltonian including these couplings, such that the nuclear dynamics is carried out on potential energy surfaces including the effects of the couplings-this is critical in any applications considering, for example, transition metal complexes or strong laser fields. We also give an overview over the new SHARC2.0 dynamics software package, released under the GNU General Public License, which implements the SHARC approach and several analysis tools. The review closes with a brief survey of applications where SHARC was employed to study the nonadiabatic dynamics of a wide range of molecular systems. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and KineticsSoftware > Simulation MethodsSoftware > Quantum Chemistry.

Charge Transfer Driven Structural Relaxation in a Push-Pull Azobenzene Dye-Semiconductor Complex

Photoexcited structural dynamics in azo-compounds may differ fundamentally whether the push-pull photochromic azo-compound is isolated or forms a heterogeneous charge transfer complex, due to a sudden oxidation of the chromophore. Herein, we use a quantum-classical self-consistent approach that incorporates nonadiabatic excited-state electronic quantum dynamics into molecular mechanics to study the photoexcited dynamics of the push-pull azo-compound para-Methyl Red in the gas phase and sensitizing the (101) anatase surface of TiO₂. We find that the photoinduced S₂/S₀ trans-to- cis isomerization of para-Methyl Red in the gas phase occurs through a pedal-like torsion around the ϕ_CNNC dihedral angle, without evidence to support the inversion mechanism, likewise in the parent azobenzene molecule. However, the photoexcited structural relaxation of the charge transfer complex para-Methyl Red/TiO₂ contrasts essentially with the isolated azo-compounds. Immediately after photoexcitation, the excited electron flows into the TiO₂ conduction band, with an injection time constant of ≃5 fs, and no indication of isomerization is observed during the 1.5 ps simulations. Instead, a strong vibronic relaxation occurs that excites the NN stretching mode of the azo-group, which is ultimately ascribed to the NA relaxation, and delocalization, of the hole wavepacket.

Trajectory surface hopping study of propane photodissociation dynamics at 157 nm

The photodissociation dynamics of propane molecules has been studied using the quasiclassical trajectory surface hopping (TSH) method in conjunction with Tully's fewest switches algorithm. The trajectories are propagated on potential energy surfaces computed on-the-fly using the multiconfiguration and multireference ab initio method starting in the lowest excited singlet state (HOMO → 3s Rydberg state) of propane at 157 nm with the emphasis on the site specificity of atomic hydrogen elimination, molecular hydrogen elimination, and their product branching ratios. Our dynamics simulation revealed that there are three primary dissociation channels: the atomic hydrogen elimination, the molecular hydrogen elimination, and the C-C bond scission. The trajectories indicate that the H₂ elimination from the internal carbon atom (2,2-H₂ elimination) and terminal carbon atom (1,1-H₂ elimination) is the major process and follows a three centred synchronous concerted mechanism. 1,2-H₂ and 1,3-H₂ eliminations on the other hand are minor processes and exclusively follow the roaming mediated nonadiabatic dynamics. The probability of elimination of the hydrogen atom from two terminal groups (terminal hydrogen elimination) is greater than that from the internal CH₂ group (internal hydrogen elimination). Almost 83% of atomic hydrogen elimination occurs through the asynchronous concerted mechanism from the terminal carbon atom via triple dissociation leading to CH₃ + C₂H₄ + H products. This finding is in good agreement with a recent experimental observation. The present TSH study indicates that approximately one-third of the trajectories those resulted in a triple dissociation channel, CH₃ + C₂H₄ + H completed in the ground singlet state following a nonadiabatic path (hopping from the first excited singlet S₁ to the ground state S₀) via the C-C and C-H dissociation coordinate conical intersection S₁/S₀. The products CH₃(1 ²A₂^″) + C₂H₄(¹A_g) + H, obtained are ground state methyl radicals and ground state ethylene. The trajectories those ended in a triple dissociation channel CH₃ + C₂H₄ + H adiabatically in the S₁ state lead to CH₃(1 ²A₂^″) + C₂H₄ (1 ³B₁) + H, where singlet methyl radicals and triplet ethylene are formed in their corresponding lowest electronic state via a spin conserving route. Two channels, CH₄ + CH₃CH and C₂H₆ + CH_2, are found to have minor contributions. In the case of methane elimination, the trajectories that follow an adiabatic path lead to CH₃CH(1 ¹A″) + CH_4,(1 ¹A₁), where ethylidene is in the excited state and methane is in the ground state. Methane elimination via nonadiabatic path leads to CH₃CH(1¹A^') + CH₄(1 ¹A₁), where both ethylidene and methane are in the ground electronic state. Ethane eliminations follow the adiabatic path leading to C₂H₆(1 ¹A_1g) + CH₂(1 ¹B₁) where ethane is in the ground state and methylene is in the first excited state.

Surface hopping dynamics including intersystem crossing using the algebraic diagrammatic construction method

We report an implementation for employing the algebraic diagrammatic construction to second order [ADC(2)] ab initio electronic structure level of theory in nonadiabatic dynamics simulations in the framework of the SHARC (surface hopping including arbitrary couplings) dynamics method. The implementation is intended to enable computationally efficient, reliable, and easy-to-use nonadiabatic dynamics simulations of intersystem crossing in organic molecules. The methodology is evaluated for the 2-thiouracil molecule. It is shown that ADC(2) yields reliable excited-state energies, wave functions, and spin-orbit coupling terms for this molecule. Dynamics simulations are compared to previously reported results using high-level multi-state complete active space perturbation theory, showing favorable agreement.