On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2)

Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.

Ultrafast decay dynamics of electronically excited 2-ethylpyrrole

The excited-state decay dynamics of 2-ethylpyrrole following UV excitation in the wavelength range of 254.8-218.0 nm is investigated in detail using the femtosecond time-resolved photoelectron imaging method. The time-resolved photoelectron spectra and photoelectron angular distributions at all pump wavelengths are carefully analysed and the following picture is derived: at the longest pump wavelengths (254.8, 248.3 and 246.1 nm), 2-ethylpyrrole is excited to the S₁(¹πσ*) state having a lifetime of about 50 fs. At 248.3, 246.1 and 237.4 nm, another excited state of Rydberg character is excited. The lifetime of this state is ∼570 fs at 237.4 nm and becomes slightly longer at other two pump wavelengths. At the shortest pump wavelengths (230.8 and 218.0 nm), 2-ethylpyrrole is excited to a state which is tentatively assigned to the 1¹ππ* state, having a lifetime of 75 ± 15 and 48 ± 10 fs for the longer and shorter pump wavelengths, respectively. Internal conversion to the S₁(¹πσ*) state might be one of the decay mechanisms of the 1¹ππ* state.

Caveat when using ADC(2) for studying the photochemistry of carbonyl-containing molecules

Several electronic-structure methods are available to study the photochemistry and photophysics of organic molecules. Among them, ADC(2) stands as a sweet spot between computational efficiency and accuracy. As a result, ADC(2) has recently seen its number of applications booming, in particular to unravel the deactivation pathways and photodynamics of organic molecules. Despite this growing success, we demonstrate here that care has to be taken when studying the nonradiative pathways of carbonyl-containing molecules, as ADC(2) appears to suffer from a systematic flaw.

Simulation of UV absorption spectra and relaxation dynamics of uracil and uracil-water clusters

Despite many studies, the mechanisms of nonradiative relaxation of uracil in the gas phase and in aqueous solution are still not fully resolved. Here we combine theoretical UV absorption spectroscopy with nonadiabatic dynamics simulations to identify the photophysical mechanisms that can give rise to experimentally observed decay time constants. We first compute and theoretically assign the electronic spectra of uracil using the second-order algebraic-diagrammatic-construction (ADC(2)) method. The obtained electronic states, their energy differences and state-specific solvation effects are the prerequisites for understanding the photodynamics. We then use nonadiabatic trajectory-surface-hopping dynamics simulations to investigate the photoinduced dynamics of uracil and uracil-water clusters. In contrast to previous studies, we found that a single mechanism - the ethylenic twist around the C[double bond, length as m-dash]C bond - is responsible for the ultrafast component of the nonradiative decay, both in the gas phase and in solution. Very good agreement with the experimentally determined ultrashort decay time constants is obtained.

Time-Resolved Photoelectron Spectroscopy Studies of Isoxazole and Oxazole

The excited state relaxation pathways of isoxazole and oxazole upon excitation with UV-light were investigated by nonadiabatic ab initio dynamics simulations and time-resolved photoelectron spectroscopy. Excitation of the bright ππ*-state of isoxazole predominantly leads to ring-opening dynamics. Both the initially excited ππ*-state and the dissociative πσ*-state offer a combined barrier-free reaction pathway, such that ring-opening, defined as a distance of more than 2 Å between two neighboring atoms, occurs within 45 fs. For oxazole, in contrast, the excited state dynamics is about twice as slow (85 fs) and the quantum yield for ring-opening is lower. This is caused by a small barrier between the ππ*-state and the πσ*-state along the reaction path, which suppresses direct ring-opening. Theoretical findings are consistent with the measured time-resolved photoelectron spectra, confirming the timescales and the quantum yields for the ring-opening channel. The results indicate that a combination of time-resolved photoelectron spectroscopy and excited state dynamics simulations can explain the dominant reaction pathways for this class of molecules. As a general rule, we suggest that the antibonding σ*-orbital located between the oxygen atom and a neighboring atom of a five-membered heterocyclic system provides a driving force for ring-opening reactions, which is modified by the presence and position of additional nitrogen atoms.

Excited-state intramolecular proton transfer driven by conical intersection in hydroxychromones

Nonradiative decay pathways associated with vibronically coupled S₁ (ππ*)-S₂ (nπ*) potential energy surfaces of 3- and 5-hydroxychromones are investigated by employing the linear vibronic coupling approach. The presence of a conical intersection close to the Franck-Condon point is identified based on the critical examination of computed energetics and structural parameters of stationary points. We show that very minimal displacements of relevant atoms of intramolecular proton transfer geometry are adequate to drive the molecule toward the conical intersection nuclear configuration. The evolving wavepacket on S₁ (ππ*) bifurcates at the conical intersection: a part of the wavepacket moves to S₂ (nπ*) within a few femtoseconds while the other decays to S₁ minimum. Our findings indicate the possibility of forming the proton transfer tautomer product via S₂ (nπ*), competing with the traditional pathway occurring on S₁ (ππ*).

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.

Improved Ehrenfest Approach to Model Correlated Electron-Nuclear Dynamics

Mixed quantum-classical mechanical descriptions are critical to modeling coupled electron-nuclear dynamics, i.e., nonadiabatic molecular dynamics, relevant to photochemical and photophysical processes. We introduce an efficient description of such dynamics in terms of an effective Hamiltonian that not only properly captures electron-nuclear correlation effects but also helps develop an efficient computational method. In particular, we introduce a coupled Gaussian wavepacket parametrization of the nuclear wave function, which generalizes the Ehrenfest approach to account for electron-nuclei correlations. We test this new approach, Ehrenfest-Plus, on a suite of model problems that probe electron-nuclear correlation in nonadiabatic transitions. The high accuracy of our approach, combined with mixed quantum-classical efficiency, opens a path for improved simulation of nonadiabatic molecular dynamics in realistic molecular systems.

Multistate hybrid time-dependent density functional theory with surface hopping accurately captures ultrafast thymine photodeactivation

We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born–Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions.

Ultrafast photodissociation dynamics of 2-ethylpyrrole: adding insight to experiment with ab initio multiple cloning

Fully quantum nonadiabatic dynamics calculation of 2-ethylpyrrole photodissociation shows the experimentally obtained ultrafast time-constant has two components.

Acetylacetone photodynamics at a seeded free-electron laser

The first steps in photochemical processes, such as photosynthesis or animal vision, involve changes in electronic and geometric structure on extremely short time scales. Time-resolved photoelectron spectroscopy is a natural way to measure such changes, but has been hindered hitherto by limitations of available pulsed light sources in the vacuum-ultraviolet and soft X-ray spectral region, which have insufficient resolution in time and energy simultaneously. The unique combination of intensity, energy resolution, and femtosecond pulse duration of the FERMI-seeded free-electron laser can now provide exceptionally detailed information on photoexcitation–deexcitation and fragmentation in pump-probe experiments on the 50-femtosecond time scale. For the prototypical system acetylacetone we report here electron spectra measured as a function of time delay with enough spectral and time resolution to follow several photoexcited species through well-characterized individual steps, interpreted using state-of-the-art static and dynamics calculations. These results open the way for investigations of photochemical processes in unprecedented detail.

The first steps in photochemical processes involve changes in electronic and geometric structure on extremely short timescales. Here, the authors report femtosecond dynamics in prototypical acetylacetone, by pump-probe photoexcitation-photoemission experiments and static and dynamics calculations.

Hydrogen Bonds and Life in the Universe

The scientific community is allocating more and more resources to space missions and astronomical observations dedicated to the search for life beyond Earth. This experimental endeavor needs to be backed by a theoretical framework aimed at defining universal criteria for the existence of life. With this aim in mind, we have explored which chemical and physical properties should be expected for life possibly different from the terrestrial one, but similarly sustained by genetic and catalytic molecules. We show that functional molecules performing genetic and catalytic tasks must feature a hierarchy of chemical interactions operating in distinct energy bands. Of all known chemical bonds and forces, only hydrogen bonds are able to mediate the directional interactions of lower energy that are needed for the operation of genetic and catalytic tasks. For this reason and because of the unique quantum properties of hydrogen bonding, the functional molecules involved in life processes are predicted to have extensive hydrogen-bonding capabilities. A molecular medium generating a hydrogen-bond network is probably essential to support the activity of the functional molecules. These hydrogen-bond requirements constrain the viability of hypothetical biochemistries alternative to the terrestrial one, provide thermal limits to life molecular processes, and offer a conceptual framework to define a transition from a “covalent-bond stage” to a “hydrogen-bond stage” in prebiotic chemistry.

Nonradiative Relaxation Mechanisms of UV Excited Phenylalanine Residues: A Comparative Computational Study

The present work is directed toward understanding the mechanisms of excited state deactivation in three neutral model peptides containing the phenylalanine residue. The excited state dynamics of theγL(g+)folded form of N-acetylphenylalaninylamide (NAPA B) and its amide-N-methylated derivative (NAPMA B) is reviewed and compared to the dynamics of the monohydrated structure of NAPA (NAPAH). The goal is to unravel how the environment, and in particular solvation, impacts the photodynamics of peptides. The systems are investigated using reaction path calculations and surface hopping nonadiabatic dynamics based on the coupled cluster doubles (CC2) method and time-dependent density functional theory. The work emphasizes the role that excitation transfer from the phenylππ*to amidenπ*state plays in the deactivation of the three systems and shows how the ease of out-of-plane distortions of the amide group determines the rate of population transfer between the two electronic states. The subsequent dynamics on thenπ*state is barrierless along several pathways and leads to fast deactivation to the ground electronic state.

Ultrafast excited-state dynamics of 2,4-dimethylpyrrole

Ultrafast excited-state dynamics of 2,4-dimethylpyrrole are studied in detail following deep UV excitation.

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.

How does tetraphenylethylene relax from its excited states?

Photocyclization play a key role in the deactivation mechanism of tetraphenylethylene.

Assessment of Approximate Coupled-Cluster and Algebraic-Diagrammatic-Construction Methods for Ground- and Excited-State Reaction Paths and the Conical-Intersection Seam of a Retinal-Chromophore Model

As a minimal model of the chromophore of rhodopsin proteins, the penta-2,4-dieniminium cation (PSB3) poses a challenging test system for the assessment of electronic-structure methods for the exploration of ground- and excited-state potential-energy surfaces, the topography of conical intersections, and the dimensionality (topology) of the branching space. Herein, we report on the performance of the approximate linear-response coupled-cluster method of second order (CC2) and the algebraic-diagrammatic-construction scheme of the polarization propagator of second and third orders (ADC(2) and ADC(3)). For the ADC(2) method, we considered both the strict and extended variants (ADC(2)-s and ADC(2)-x). For both CC2 and ADC methods, we also tested the spin-component-scaled (SCS) and spin-opposite-scaled (SOS) variants. We have explored several ground- and excited-state reaction paths, a circular path centered around the S1/S0 surface crossing, and a 2D scan of the potential-energy surfaces along the branching space. We find that the CC2 and ADC methods yield a different dimensionality of the intersection space. While the ADC methods yield a linear intersection topology, we find a conical intersection topology for the CC2 method. We present computational evidence showing that the linear-response CC2 method yields a surface crossing between the reference state and the first response state featuring characteristics that are expected for a true conical intersection. Finally, we test the performance of these methods for the approximate geometry optimization of the S1/S0 minimum-energy conical intersection and compare the geometries with available data from multireference methods. The present study provides new insight into the performance of linear-response CC2 and polarization-propagator ADC methods for molecular electronic spectroscopy and applications in computational photochemistry.