Miquel Huix-Rotllant

Shape of Multireference, Equation-of-Motion Coupled-Cluster, and Density Functional Theory Potential Energy Surfaces at a Conical Intersection.

We report and characterize ground-state and excited-state potential energy profiles using a variety of electronic structure methods along a loop lying on the branching plane associated with a conical intersection (CI) of a reduced retinal model, the penta-2,4-dieniminium cation (PSB3). Whereas the performance of the equation-of-motion coupled-cluster, density functional theory, and multireference methods had been tested along the excited- and ground-state paths of PSB3 in our earlier work, the ability of these methods to correctly describe the potential energy surface shape along a CI branching plane has not yet been investigated. This is the focus of the present contribution. We find, in agreement with earlier studies by others, that standard time-dependent DFT (TDDFT) does not yield the correct two-dimensional (i.e., conical) crossing along the branching plane but rather a one-dimensional (i.e., linear) crossing along the same plane. The same type of behavior is found for SS-CASPT2(IPEA=0), SS-CASPT2(IPEA=0.25), spin-projected SF-TDDFT, EOM-SF-CCSD, and, finally, for the reference MRCISD+Q method. In contrast, we found that MRCISD, CASSCF, MS-CASPT2(IPEA=0), MS-CASPT2(IPEA=0.25), XMCQDPT2, QD-NEVPT2, non-spin-projected SF-TDDFT, and SI-SA-REKS yield the expected conical crossing. To assess the effect of the different crossing topologies (i.e., linear or conical) on the PSB3 photoisomerization efficiency, we discuss the results of 100 semiclassical trajectories computed by CASSCF and SS-CASPT2(IPEA=0.25) for a PSB3 derivative. We show that for the same initial conditions, the two methods yield similar dynamics leading to isomerization quantum yields that differ by only a few percent.

Assessment of Density Functional Theory for Describing the Correlation Effects on the Ground and Excited State Potential Energy Surfaces of a Retinal Chromophore Model

In the quest for a cost-effective level of theory able to describe a large portion of the ground and excited potential energy surfaces of large chromophores, promising approaches are rooted in various approximations to the exact density functional theory (DFT). In the present work, we investigate how generalized Kohn-Sham DFT (GKS-DFT), time-dependent DFT (TDDFT), and spin-restricted ensemble-DFT (REKS) methods perform along three important paths characterizing a model retinal chromophore (the penta-2,4-dieniminium cation) in a region of near-degeneracy (close to a conical intersection) with respect to reference high-level multiconfigurational wave function methods. If GKS-DFT correctly describes the closed-shell charge transfer state, only TDDFT and REKS approaches give access to the open-shell diradical, one which sometimes corresponds to the electronic ground state. It is demonstrated that the main drawback of the usual DFT-based methods lies in the absence of interactions between the charge transfer and the diradicaloid configurations. Hence, we test a new computational scheme based on the State-averaged REKS (SA-REKS) approach, which explicitly includes these interactions into account. The State-Interaction SA-REKS (SI-SA-REKS) method significantly improves on the REKS and the SA-REKS results for the target system. The similarities and differences between DFT and wave function-based approaches are analyzed according to (1) the active space dimensions of the wave function-based methods and (2) the relative electronegativities of the allyl and protonated Schiff base moieties.

Concurrent Effects of Delocalization and Internal Conversion Tune Charge Separation at Regioregular Polythiophene-Fullerene Heterojunctions

Quantum-dynamical simulations are used to investigate the interplay of exciton delocalization and vibronically induced internal conversion processes in the elementary charge separation steps at regioregular donor-acceptor heterojunctions. Ultrafast internal conversion leads to efficient deexcitation within the excitonic and charge transfer manifolds, thus modifying the charge separation dynamics. We address a model donor-acceptor junction representative of regioregular P3HT-PCBM, using high-dimensional quantum dynamics simulations by multiconfigurational methods. While partial trapping into an interfacial charge separated state occurs, long-range charge-separated states are accessed as previously demonstrated in the work of Tamura and Burghardt [J. Am. Chem. Soc. 2013, 135, 16364]. For an H-aggregate type, stacked donor species, the initial bright state undergoes ultrafast internal conversion within the excitonic manifold, creating multiple charge transfer pathways before reaching the lowest-energy dark exciton, which is uncoupled from the charge transfer manifold. This process profoundly affects the charge separation mechanism and efficiency. For small energetic offsets between the interfacial excitonic and charge transfer states, a delocalized initial bright state proves less prone to electron-hole capture by the interfacial trap than a localized, vibronic wavepacket close to the interface. For both delocalized and localized initial states, a comparable yield of free carriers is obtained, which is found to be optimal for energetic offsets of the order of the Coulomb barrier to charge separation. Interfacial trapping is significantly reduced as the barrier height decreases with fullerene aggregation. Despite the high-dimensional nature of the system, charge separation is an ultrafast coherent quantum process exhibiting oscillatory features as observed in recent experiments.

Assessment of density functional theory based ΔSCF (self-consistent field) and linear response methods for longest wavelength excited states of extended π-conjugated molecular systems

Computational investigation of the longest wavelength excitations in a series of cyanines and linear n-acenes is undertaken with the use of standard spin-conserving linear response time-dependent density functional theory (TD-DFT) as well as its spin-flip variant and a ΔSCF method based on the ensemble DFT. The spin-conserving linear response TD-DFT fails to accurately reproduce the lowest excitation energy in these π-conjugated systems by strongly overestimating the excitation energies of cyanines and underestimating the excitation energies of n-acenes. The spin-flip TD-DFT is capable of correcting the underestimation of excitation energies of n-acenes by bringing in the non-dynamic electron correlation into the ground state; however, it does not fully correct for the overestimation of the excitation energies of cyanines, for which the non-dynamic correlation does not seem to play a role. The ensemble DFT method employed in this work is capable of correcting for the effect of missing non-dynamic correlation in the ground state of n-acenes and for the deficient description of differential correlation effects between the ground and excited states of cyanines and yields the excitation energies of both types of extended π-conjugated systems with the accuracy matching high-level ab initio multireference calculations.

Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules.

Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems.

Theoretical study of the photochemical generation of triplet acetophenone

Acetophenone has a rich photochemistry, which strongly depends on the absorbing state. For example, the excitation to the lowest singlet excited state (S1) leads to a triplet population with a phosphorescence quantum yield of one, while the excitation to S2 leads to photocleavage reactions. Here, we rationalize the photochemistry of acetophenone after being absorbed into the S1, S2 and S3 states by performing a systematic study of all the singlet and triplet minimum energy structures and state crossings between the relevant electronic states. We calculate these structures at the complete-active space self-consistent field (CASSCF) level of theory and at the correlated extended second-order quasi-degenerate multi-reference perturbation theory (XMCQDPT2), emphasizing the importance of correlation effects in the determination of structures.

Stabilization of Mixed Frenkel-Charge Transfer Excitons Extended Across Both Strands of Guanine–Cytosine DNA Duplexes

The photoreactive pathways that may lead to DNA damage depend crucially upon the nature of the excited electronic states. The study of alternating guanine-cytosine duplexes by fluorescence spectroscopy and quantum mechanical calculations identifies a novel type of excited states that can be populated following UVB excitation. These states, denoted High-energy Emitting Long-lived Mixed (HELM) states, extend across both strands and arise from mixing between cytosine Frenkel excitons and guanine-to-cytosine charge transfer states. They emit at energies higher than ππ* states localized on single bases, survive for several nanoseconds, are sensitive to the ionic strength of the solution, and are strongly affected by the structural transition from the B form to the Z form. Their impact on the formation of lesions of the genetic code needs to be assessed.

Density-functional methods for excited states

Current status and recent developments in linear response TD-DFT.- Density matrix functional theory (DMFT) and linear response time-dependent DMFT (TD-DMFT) for excited states.- Ensemble DFT for strongly correlated molecules and excited states.- Dynamical processes in open quantum systems from a TD-DFT perspective.- Computational Molecular Spectroscopy with TD-DFT.- Non-linear response and multi-dimensional spectroscopy using time-dependent DFT.- Absorption spectroscopy, emissive properties and ultra-fast intersystem crossing processes in transition metal complexes: TD-DFT and Spin-Orbit Coupling.- Excitons and extended systems in TD-DFT.- Conical intersections with DFT methods.- Multi-scale modeling of photochemistry with TD-DFT.- Semi-classical surface-hopping dynamics using potential energy surfaces from time-dependent DFT calculations.

Description of Conical Intersections with Density Functional Methods

Conical intersections are perhaps the most significant mechanistic features of chemical reactions occurring through excited states. By providing funnels for efficient non-adiabatic population transfer, conical intersections govern the branching ratio of products of such reactions, similar to what the transition states do for ground-state reactivity. In this regard, intersections between the ground and the lowest excited states play a special role, and the correct description of the potential energy surfaces in their vicinity is crucial for understanding the mechanism and dynamics of excited-state reactions. The methods of density functional theory, such as time-dependent density functional theory, are widely used to describe the excited states of large molecules. However, are these methods suitable for describing the conical intersections or do they lead to artifacts and, consequently, to erroneous description of reaction dynamics? Here we address the first part of this question and analyze the ability of several density functional approaches, including the linear-response time-dependent approach as well as the spin-flip and ensemble formalisms, to provide the correct description of conical intersections and the potential energy surfaces in their vicinity. It is demonstrated that the commonly used linear-response time-dependent theory does not yield a proper description of these features and that one should instead use alternative computational approaches.

Photophysics of Acetophenone Interacting with DNA: Why the Road to Photosensitization is Open

Deoxyribonucleic acid photosensitization, i.e. the photoinduced electron‐ or energy‐transfer of chromophores interacting with DNA, is a crucial phenomenon that triggers important DNA lesions such as pyrimidine dimerization, even upon absorption of relatively low‐energy radiation. Oxidative lesions may also be produced via the photoinduced production of reactive oxygen species. Aromatic ketones, and acetophenone in particular, are well known for their sensitization effects. In this contribution we model the structural and dynamical properties of the acetophenone/DNA aggregates as well as their spectroscopic and photophysical properties using high‐level hybrid quantum mechanics/molecular mechanics methods. We show that the key steps of the photochemistry of acetophenone in gas phase are conserved in the macromolecular environment and thus an ultrafast singlet–triplet conversion of acetophenone is expected prior to the transfer to DNA.