Surface hopping investigation of the relaxation dynamics in radical cations

Exact-Factorization-Based Surface Hopping without Velocity Adjustment

While surface hopping has emerged as a powerful method for simulating non-adiabatic dynamics in large molecules, the ad hoc nature of the necessary velocity adjustments and decoherence corrections in the algorithm somewhat reduces its reliability. Here we propose a new scheme that eliminates these aspects by combining the nuclear equation from the quantum-trajectory surface-hopping approach with the electronic equation derived from the exact-factorization approach. The resulting method, denoted QTSH-XF, yields a surface-hopping method on firmer ground than previous and is shown to successfully capture dynamics in Tully models and in a linear vibronic coupling model of the photoexcited uracil cation.

Significance of Energy Conservation in Coupled-Trajectory Approaches to Nonadiabatic Dynamics

Through approximating electron-nuclear correlation terms in the exact factorization approach, trajectory-based methods have been derived and successfully applied to the dynamics of a variety of light-induced molecular processes, capturing quantum (de)coherence effects rigorously. These terms account for the coupling among the trajectories, recovering the nonlocal nature of quantum nuclear dynamics that is completely overlooked in traditional independent-trajectory algorithms. Nevertheless, some of the approximations introduced in the derivation of some of these methods do not conserve the total energy. We analyze energy conservation in the coupled-trajectory mixed quantum-classical (CTMQC) algorithm and explore the performance of a modified algorithm, CTMQC-E, where some of the terms are redefined to restore energy conservation. A set of molecular models is used as a test, namely, 2-cis-penta-2,4-dienimium cation, bis(methylene) adamantyl radical cation, butatriene cation, uracil radical cation, and neutral pyrazine.

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

The exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron-ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron-nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

Excited-state dynamics of o-nitrophenol studied with UV pump-VUV probe time-resolved photoelectron and photoion spectroscopy

Time-resolved photoionization measurements were performed on o-nitrophenol pumped with UV laser pulses at a central wavelength of 255 nm (4.9 eV) and probed with vacuum ultraviolet (VUV) pulses at 153 nm (8.1 eV). The photoelectron spectrum and time of flight mass spectrum for ions were recorded at each pump-probe delay. The measurements are interpreted with the aid of electronic structure calculations for both the neutral and ionic states. Evidence is found for the formation of a bicyclic intermediate followed by NO dissociation through a process of internal conversion and intersystem crossing. The combination of photoelectron and photoion spectroscopy, together with computational results, provides strong evidence of intersystem crossing that is difficult to establish with only a single technique.

Nonadiabatic Excited State Dynamics of Organic Chromophores: Take-Home Messages

Nonadiabatic excited state dynamics are important in a variety of processes. Theoretical and experimental developments have allowed for a great progress in this area, while combining the two is often necessary and the best approach to obtain insight into the photophysical behavior of molecules. In this Feature Article we use examples of our recent work combining time-resolved photoelectron spectroscopy with theoretical nonadiabatic dynamics to highlight important lessons we learned. We compare the nonadiabatic excited state dynamics of three different organic molecules with the aim of elucidating connections between structure and dynamics. Calculations and measurements are compared for uracil, 1,3-cyclooctadiene, and 1,3-cyclohexadiene. The comparison highlights the role of rigidity in influencing the dynamics and the difficulty of capturing the dynamics accurately with calculations.

Accurate Non-Adiabatic Couplings from Optimally-Tuned Range-Separated Hybrid Functionals

Precise theoretical calculations of non-adiabatic couplings, which describe the interaction between two Born-Oppenheimer surfaces, are important for the modeling of radiationless decay mechanisms in photochemical processes. Here, we demonstrate that accurate non-adiabatic couplings can be calculated in the framework of linear-response time-dependent density functional theory by using non-empirical, optimally-tuned range-separated hybrid (OT-RSH) functionals. We focus on molecular radicals, in which ultrafast non-radiative decay plays a crucial role, to find that the OT-RSH functional compares well to wave-function based reference data and competes with the accuracy of semi-empirical CAM-B3LYP calculations. Our findings show that the OT-RSH approach provides very accurate non-adiabatic couplings and, therefore, provides a computationally efficient alternative to wave-function based techniques.

Stable excited dication: trapping on the S1 state of formaldehyde dication after strong field ionization

Combined theoretical and experimental work examines the dynamics of dication formaldehyde produced by strong field ionization. Trajectory surface hopping dynamics on the first several singlet electronic states of the formaldehyde dication are used to examine the relaxation pathways and dissociation channels, while kinetic energy distributions after strong field ionization of formaldehyde and deuterated formaldehyde are used to confirm the theoretical predictions. We find that the first excited state of the formaldehyde dication is stable, neither decays to the ground state nor dissociates, even though the ground state and higher lying states are directly dissociative. The stability of the first excited state is explained by its symmetry which does not allow for radiative or nonradiative transitions to the ground state and by large barriers to dissociate on the excited state surface.

Conformer-Specific Dissociation Dynamics in Dimethyl Methylphosphonate Radical Cation

The dynamics of the dimethyl methylphosphonate (DMMP) radical cation after production by strong field adiabatic ionization have been investigated. Pump-probe experiments using strong field 1300 nm pulses to adiabatically ionize DMMP and a 800 nm non-ionizing probe induce coherent oscillations of the parent ion yield with a period of about 45 fs. The yields of two fragments, PO2C2H7+ and PO2CH4+, oscillate approximately out of phase with the parent ion, but with a slight phase shift relative to each other. We use electronic structure theory and nonadiabatic surface hopping dynamics to understand the underlying dynamics. The results show that while the cation oscillates on the ground state along the P=O bond stretch coordinate, the probe excites population to higher electronic states that can lead to fragments PO2C2H7+ and PO2CH4+. The computational results combined with the experimental observations indicate that the two conformers of DMMP that are populated under experimental conditions exhibit different dynamics after being excited to the higher electronic states of the cation leading to different dissociation products. These results highlight the potential usefulness of these pump-probe measurements as a tool to study conformer-specific dynamics in molecules of biological interest.

Electronic coherences in nonadiabatic molecular photophysics revealed by time-resolved photoelectron spectroscopy

Time-resolved photoelectron spectroscopy (TRPES) is a promising technique for the study of ultrafast molecular processes, such as the nonadiabatic dynamics taking place at conical intersections. Directly accessing the evolution of the coherences generated at the conical intersection should provide most valuable dynamical information. However, the signals are dominated by background contributions due to state populations, and most theoretical treatments completely neglect the role of the coherences. Here we show that distinguishable signatures of molecular coherences appear in TRPES. These can be recorded using currently available ultrashort pulses and unambiguously extracted at the postprocessing stage. The technique thus provides direct access to nonadiabatic coherence dynamics.

Exact-Factorization-Based Surface Hopping for Multistate Dynamics

A surface-hopping algorithm recently derived from the exact factorization approach, SHXF [Ha et al. J. Phys. Chem. Lett. 2018, 9, 1097], introduces an additional term in the electronic equation of surface hopping that couples electronic states through the quantum momentum. This term not only provides a first-principles description of decoherence, but here we show it is crucial to accurately capture nonadiabatic dynamics when more than two states are occupied at any given time. Using a vibronic coupling model of the uracil cation, we show that the lack of this term in traditional surface-hopping methods, including those with decoherence corrections, leads to failure to predict the dynamics through a three-state intersection, while SHXF performs similarly to the multiconfiguration time-dependent Hartree quantum dynamics benchmark.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

Time Resolved Photoelectron Spectroscopy as a Test of Electronic Structure and Nonadiabatic Dynamics

We compare different levels of theory for simulating excited state molecular dynamics and use time-resolved photoelectron spectroscopy measurements to benchmark the theory. We perform trajectory surface hopping simulations for uracil excited to the first bright state (ππ*) using three different levels of theory (CASSCF, MRCIS, and XMS-CASPT2) in order to understand the role of dynamical correlation in determining the excited state dynamics, with a focus on the coupling between different electronic states and internal conversion back to the ground state. These dynamics calculations are used to simulate the time-resolved photoelectron spectra. The comparison of the calculated and measured spectra allows us to draw conclusions regarding the relative insights and quantitative accuracy of the calculations at the three different levels of theory, demonstrating that detailed quantitative comparisons of time-resolved photoelectron spectra can be used to benchmark methodology.

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.

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.

Ultrafast and radiationless electronic excited state decay of uracil and thymine cations: computing the effects of dynamic electron correlation

In this article we characterise the radiationless decay of the first few electronic excited states of the cations of DNA/RNA nucleobases uracil and thymine, including the effects of dynamic electron correlation on energies and geometries (optimised with XMS-CASPT2). In both systems, we find that one state of ²n and another two of ²π⁺ character can be populated following photoionisation, and their different minima and interstate crossings are located. We find strong similarities between uracil and thymine cations: with accessible conical intersections suggesting that depopulation of their electronic excited states takes place on ultrafast timescales in both systems, suggesting that they are photostable in agreement with previous theoretical (uracil⁺) evidence. We find that dynamic electron correlation separates the energy levels of the "3-state" conical intersection (D₂/D₁/D₀)_CI previously located with CASSCF for uracil⁺, which will therefore have a different geometry and higher energy. Simulating the electronic and vibrational absorptions allows us to characterise spectral fingerprints that could be used to monitor these cation photo-processes experimentally.

A molecular beam apparatus for performing single photon initiated dissociative rearrangement reactions (SPIDRR) with transition metal cation bound organic clusters

The study of gas-phase ion-molecule reactions has been influential in the investigation of transition metal mediated bond activation and catalysis. We have furthered this field by developing a new technique capable of measuring the microcanonical kinetics for reactions between transition metal cations and neutral organic molecules. This novel method has been designated as single photon initiated dissociative rearrangement reaction (SPIDRR) technique and provides a nearly direct measurement of microcanonical reaction rate constants. For this reason, SPIDRR offers unique insight into reaction mechanisms and dynamics by assessing the energy dependence of the microcanonical rate constant, as well as measuring product branching fractions and kinetic isotope effects. The following paper provides a detailed overview of SPIDRR and its advantages in the field of gas-phase catalysis research.

Coexistence of Different Electron-Transfer Mechanisms in the DNA Repair Process by Photolyase

DNA photolyase has been the topic of extensive studies due to its important role of repairing photodamaged DNA, and its unique feature of using light as an energy source. A crucial step in the repair by DNA photolyase is the forward electron transfer from its cofactor (FADH(-) ) to the damaged DNA, and the detailed mechanism of this process has been controversial. In the present study, we examine the forward electron transfer in DNA photolyase by carrying out high-level ab initio calculations in combination with a quantum mechanical/molecular mechanical (QM/MM) approach, and by measuring fluorescence emission spectra at low temperature. On the basis of these computational and experimental results, we demonstrate that multiple decay pathways exist in DNA photolyase depending on the wavelength at excitation and the subsequent transition. This implies that the forward electron transfer in DNA photolyase occurs not only by superexchange mechanism but also by sequential electron transfer.

Molecular Double Ionization Using Strong Field Few-Cycle Laser Pulses

We study strong field double ionization of a series of organic molecules by making use of coincidence detection of fragment ions. We measure the double ionization yield as a function of pulse duration, intensity, polarization, and molecular conjugation. For conjugated molecules we find strong enhancement in the double ionization rate over what one would expect on the basis of tunneling or multiphoton ionization rates. Calculations reveal a correlation between the electronic structure of the different molecules and the observed double ionization yields, highlighting the removal of electrons from inner orbitals.