Search

The inversion of the energies of the lowest singlet (S1) and triplet (T1) excited states in violation of Hund’s multiplicity rule is a rare phenomenon in stable organic molecules. S1-T1 inversion has significant consequences for the photophysics and photochemistry of organic chromophores. In this work, ab initio computational methods were employed to explore the possibility of S1-T1 inversion in hexagonal polycyclic aromatic and heteroaromatic compounds. Although the singlet-triplet energy gap ΔST = ES1 – ET1 decreases with increasing size of hexagonal polycyclic aromatics, it remains positive up to kekulene (19 rings). However, symmetric substitution of C-C pairs by B-N groups in the interior, keeping the conjugation of the outer rim intact, results in compounds with robustly negative ΔST. The non-overlapping pattern of the densities of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) is identified as the decisive criterion for S1-T1 inversion. These findings establish the existence of a new family of boron carbon nitrides with inverted singlet-triplet gaps.

Ab initio surface-hopping dynamics simulations provide temporally and structurally resolved insight into the photoinduced proton-coupled electron transfer reaction in hydrogen-bonded complexes of trianisoleheptazine with phenol and methoxyphenol.

In this perspective, we discuss a novel water-splitting scenario, namely the direct oxidation of water molecules by organic photooxidants in hydrogen-bonded chromophore-water complexes. In comparison with the established scenario of semiconductor-based water splitting, the distance of electron transfer processes is thereby reduced from mesoscopic scales to the {\AA}ngstr{\o}m scale and the time scale is reduced from milliseconds to femtoseconds, which suppresses competing loss processes. The concept is illustrated by computational studies for the heptazine-H2O complex. The excited-state landscape of this complex has been characterized with ab initio electronic-structure methods and the proton-coupled electron-transfer dynamics has been explored with nonadiabatic dynamics simulations. A unique feature of the heptazine chromophore is the existence of a low-lying and exceptionally long-lived 1{\pi}{\pi}* state in which a substantial part of the photon energy can be stored for hundreds of nanoseconds and is available for the oxidation of water molecules. The calculations reveal that the absorption spectra and the photochemical functionalities of heptazine chromophores can be systematically tailored by chemical substitution. The options of harvesting hydrogen and the problems posed by the high reactivity of OH radicals are discussed., Comment: 24 pages, 2 tables, 8 figures

Recently, a derivative of the heptazine (tris-triazine) molecule, trianisole-heptazine (TAHz), was synthesized and was shown to catalyze the oxidation of water to hydroxyl radicals under 365 nm LED light in a homogeneous reaction (E. J. Rabe et al., J. Phys. Chem. Lett. 2018, 9, 6257-6261). The possibility of water photo-oxidation with a precisely defined molecular catalyst in neat solvents opens new perspectives for clarifying the fundamental reaction mechanisms involved in water oxidation photocatalysis. In the present work, the effects of chemical substituents on the three CH positions of Hz on the photocatalytic reactivity were explored with wave function-based ab initio electronic-structure calculations for hydrogen-bonded complexes of Hz and three selected Hz derivatives (TAHz, trichloro-Hz, and tricyano-Hz) with a water molecule. While anisole is an electron-donating substituent, Cl is a weakly electron-withdrawing substituent and CN is a strongly electron-withdrawing substituent. It is shown that the barrier for the photoinduced abstraction of an H atom from the water molecule is raised (lowered) by electron-donating (electron-withdrawing) substituents. The highly mobile and reactive hydroxyl radicals generated by water oxidation can recombine with the reduced chromophore radicals to yield photohydrates. The effect of substituents on the thermodynamics of the photohydration reaction was computed. Among the four chromophores studied, TAHz stands out on account of the metastability of its photohydrate, which suggests self-healing of the photocatalyst after oxidation of TAHzH radicals by OH radicals. In addition, the effect of substituents on the H atom photodetachment reaction from the reduced chromophores, which closes the catalytic cycle, has been investigated. The energy of the repulsive 2πσ* state, which drives the photodetachment reaction is lowered (raised) by electron-donating (electron withdrawing) substituents. All four chromophores exhibit inverted S1/T1 gaps. This feature eliminates long-lived triplet states and thus avoids the activation of molecular oxygen to highly reactive singlet oxygen.

We demonstrate chemical tuning and laser-driven control of intermolecular H atom abstraction from protic solvent molecules. Using multipulse ultrafast pump-push-probe transient absorption (TA) spec...

UV irradiation of RNA leads to the formation of intra- and inter-strand crosslinks of cyclobutane type. Despite the importance of this reaction, relatively little is known about how the mutual orientation of the two bases affects the outcome of the reaction. Here we report a comparative nonadiabatic molecular dynamics study of face-to-back (F2B) and face-to-face (F2F) stacked uracil-water clusters. The computations were performed using the second-order algebraic-diagrammatic-construction (ADC(2)) method. We found that F2B stacked uracil-water clusters either relax non-reactively to the ground state by an ethylenic twist around the CC bond or remain in the lowest nπ* state in which the two bases gradually move away from each other. This finding is consistent with the low propensity for the formation of intra-strand cyclobutane dimers between adjacent RNA bases. On the contrary, in F2F stacked uracil-water clusters, in addition to non-reactive deactivation, we found a pro-reactive deactivation pathway, which may lead to the formation of cyclobutane uracil dimers in the electronic ground state. On a qualitative level, the observed photodynamics of F2F stacked uracil-water clusters explains the greater propensity of RNA to form inter-strand cyclobutane-type crosslinks.

Two-dimensional (2D) electronic spectroscopy is a powerful nonlinear technique which provides spectroscopic information on two frequency axes as well as dynamical information as a function of the so-called waiting time. Herein, an ab initio theoretical framework for the simulation of electronic 2D spectra has been developed. The method is based on the classical approximation to the doorway-window representation of three-pulse photon-echo signals and the description of nuclear motion by classical trajectories. Nonadiabatic effects are taken into account by a trajectory surface-hopping algorithm. 2D electronic spectra were simulated with ab initio on-the-fly trajectory calculations using the ADC(2) electronic-structure method for the pyrazine molecule, which is a benchmark system for ultrafast radiationless decay through conical intersections. It is demonstrated that 2D spectroscopy with subfemtosecond UV pulses can provide unprecedented detailed information on the ultrafast photodynamics of polyatomic molecules.

In the past decade, polymeric carbon nitrides consisting of heptazine (Hz) building blocks emerged as highly promising materials for photocatalytic hydrogen evolution from water or sacrificial electron donors with near-ultraviolet light. However, the complexity of these materials and their poor characterization at the atomic level are detrimental to the unraveling of the detailed mechanisms of the reactions leading to hydrogen evolution. Recently, it has been shown that a derivative of the Hz molecule, trianisole-heptazine (TAHz), is a potent photobase, which readily oxidizes various derivatives of phenol and even water by an excited-state proton-coupled electron-transfer (PCET) reaction. Energy profiles along minimum-energy reaction paths and relaxed PCET potential-energy surfaces, which previously were computed with ab initio electronic-structure methods for complexes of Hz and TAHz with protic substrates, led to qualitative insights. To obtain more quantitative insight, reaction dynamics simulations are required. In the present work, the time scales of the electron and proton transfer processes and the branching ratios of competing channels were explored with ab initio on-the-fly quasiclassical surface-hopping trajectory simulations for the hydrogen-bonded Hz-H2O complex. By the analysis of representative trajectories, detailed insight into the interplay of various nonadiabatic electronic transitions, electron transfer, proton transfer, and vibrational energy relaxation is obtained. The HzH radicals which are formed by the photoreduction of Hz can disproportionate to Hz and HzH2 in an exothermic reaction with a low reaction barrier. The time scales and branching ratios of competing channels in H-atom photodetachment from the HzH2 molecule also were explored with ab initio surface-hopping simulations. These results delineate for the first time a quantitatively supported scenario of water oxidation and hydrogen evolution with a molecular carbon nitride photocatalyst.

A computational protocol for simulating time-resolved photoelectron signals of medium-sized molecules is presented. The procedure is based on a trajectory surface-hopping description of the excited-state dynamics and a combined Dyson orbital and multicenter B-spline approach for the computation of cross sections and asymmetry parameters. The accuracy of the procedure has been illustrated for the case of ultrafast internal conversion of gas-phase pyrazine excited to the 1B2u(ππ*) state. The simulated spectra and the asymmetry map are compared to the experimental data, and a very good agreement was obtained without applying any energy-dependent rescaling or broadening. An interesting side result of this work is the finding that the signature of the 1Au(nπ*) state is indistinguishable from that of the 1B3u(nπ*) state in the time-resolved photoelectron spectrum. By locating four symmetrically equivalent minima on the lowest-excited (S1) adiabatic potential energy surface of pyrazine, we revealed the strong vibronic coupling of the 1Au(nπ*) and 1B3u(nπ*) states near the S1 ← S0 band origin.

To inform prospective design rules for controlling aza-arene photochemistry, we studied hydrogen-bonded complexes of 2,5,8-tris(4-methoxyphenyl)-1,3,4,6,7,9,9b-heptaazaphenalene (TAHz), a molecular...

Monitoring of wave packet dynamics at conical intersections by time- and frequency-resolved fluorescence spectroscopy has been investigated theoretically for a three-state two-mode model of a conical intersection coupled to a dissipative environment. The ideal and the actually measurable time- and frequency-gated fluorescence spectra are accurately and efficiently simulated by combining the hierarchy equations-of-motion method for dissipative quantum dynamics with the methodology of the equation-of-motion phase-matching approach for the calculation of spectroscopic signals. It is shown that time- and frequency-resolved fluorescence spectra reveal essential aspects of the wave packet dynamics at conical intersections and the effects of environment-induced dissipation. The results of the present work indicate that fluorescence up-conversion spectroscopy with femtosecond time resolution is an efficient tool for the characterization of ultrafast dynamics at conical intersections. Accepted version L.P.C. acknowledges support by a postdoctoral fellowship of the Alexander von Humboldt-Foundation. M.F.G. and W.D. acknowledge support from the Deutsche Forschungsgemeinschaft through a research grant and through the DFG Cluster of Excellence Munich-Centre for Advanced Photonics (http:// www.munich-photonics.de).

The mechanisms of photoinduced reactions of adenine with water molecules in hydrogen-bonded adenine-water complexes were investigated with ab initio wave-function-based electronic-structure calculations. Two excited-state electron/proton transfer reaction mechanisms have been characterized: H-atom abstraction from water by photoexcited adenine as well as H-atom transfer from photoexcited adenine or the (adenine+H) radical to water. In the water-to-adenine H-atom transfer reaction, an electron from one of the p orbitals of the water molecule fills the hole in the n (π) orbital of the nπ* (ππ*) excited state of adenine, resulting in a charge-separated electronic state. The electronic charge separation is neutralized by the transfer of a proton from the water molecule to adenine, resulting in the (adenine+H)OH biradical in the electronic ground state. In the adenine-to-water H-atom transfer reaction, πσ* states localized at the acidic sites of adenine provide the mechanism for the photoejection of an electron from adenine, which is followed by proton transfer to the hydrogen-bonded water molecule, resulting in the (adenine-H)H3O biradical. The energy profiles of the photoreactions have been computed as relaxed scans with the ADC(2) electronic-structure method. These reactions, which involve the reactivity of adenine with hydrogen-bonded water molecules, compete with the well-established intrinsic excited-state deactivation mechanisms of adenine via ring-puckering or ring-opening conical intersections. By providing additional decay channels, the electron/proton exchange reactions with water can account for the observed significantly shortened excited-state lifetime of adenine in aqueous environments. These findings indicate that adenine possibly was not only a photostabilizer at the beginning of life, but also a primordial photocatalyst for water splitting.

An ab initio theoretical framework for the simulation of femtosecond time-resolved transient absorption (TA) pump-probe (PP) spectra with quasi-classical trajectories is presented. The simulations are based on the classical approximation to the doorway-window (DW) representation of third-order four-wave-mixing signals. The DW formula accounts for the finite duration and spectral shape of the pump and probe pulses. In the classical DW formalism, classical trajectories are stochastically sampled from a positive definite doorway distribution, and the signals are evaluated by averaging over a positive definite window distribution. Nonadiabatic excited-state dynamics is described by a stochastic surface-hopping algorithm. The method has been implemented for the pyrazine molecule with the second-order algebraic-diagrammatic construction (ADC(2)) ab initio electronic-structure method. The methodology is illustrated by ab initio simulations of the ground-state bleach, stimulated emission, and excited-state absorption contributions to the TA PP spectrum of gas-phase pyrazine.

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.

We studied spectroscopic signatures of the nonadiabatic dynamics at conical intersections formed by the lowest excited singlet states in pyrazine. We considered two ab initio models of conical intersections in the excited states of pyrazine developed by Sala et al. [Phys. Chem. Chem. Phys. 16, 15957 (2014)]: a two-state (B2u and B3u), five-mode model and a three-state (B2u, B3u, and Au), nine-mode model. We simulated the signals of three widely used techniques: time- and frequency-resolved fluorescence spectroscopy, transient absorption pump–probe spectroscopy, and electronic two-dimensional spectroscopy. The signals were calculated through third-order response functions, which, in turn, were evaluated with the numerically accurate multiple Davydov ansatz. We establish spectroscopic signatures of the optically dark Au state and demonstrate that the key features of the photoinduced dynamics, such as electronic/nuclear populations, electronic/nuclear coherences, and electronic/nuclear energy transfer processes, are imprinted in the spectroscopic signals. We show that a fairly complete picture of the nonadiabatic dynamics at conical intersections can be obtained when several spectroscopic techniques are combined. Provided that the time resolution is sufficient, time- and frequency-resolved fluorescence may provide the best visualization of the nonadiabatic dynamics near conical intersections.

The effects of high pulse intensity and chirp on two-dimensional electronic spectroscopy signals are experimentally investigated in the highly non-perturbative regime using atomic rubidium vapor as clean model system. Data analysis is performed based on higher-order Feynman diagrams and non-perturbative numerical simulations of the system response. It is shown that higher-order contributions may lead to a fundamental change of the static appearance and beating-maps of the 2D spectra and that chirped pulses enhance or suppress distinct higher-order pathways. We further give an estimate of the threshold intensity beyond which the high-intensity effects become visible for the system under consideration.

International audience; The photocatalytic oxidation of water with molecular or polymeric N-heterocyclic chromophores is a topic of high current interest in the context of artificial photosynthesis, that is, the conversion of solar energy to clean fuels. Hydrogen-bonded clusters of N-heterocycles with water molecules in a molecular beam are simple model systems for which the basic mechanisms of photochemical water oxidation can be studied under well-defined conditions. In this work, we explored the photoinduced H-atom transfer reaction in pyrimidine-water clusters yielding pyrimidinyl and hydroxyl radicals with laser spectroscopy, mass spectrometry and trajectory-based ab initio molecular dynamics simulations. The oxidation of water by photoexcited pyrimidine is unequivocally confirmed by the detection of the pyrimidinyl radical. The dynamics simulations provide information on the time scales and branching ratios of the reaction. While relaxation to local minima of the S1 potential-energy surface is the dominant reaction channel, the H-atom transfer reaction occurs on ultrafast time scales (faster than about 100 fs) with a branching ratio of a few percent. From the relaxed population in the S1 state, H-atom transfer yielding pyrimidinyl and hydroxyl radicals can occur on much longer time scales (picoseconds to nanoseconds) by H-atom tunneling.

Electron diffraction correlates the excited-state decay of pyridine with its ring distortion

We present a conspectus of recent joint spectroscopic and computational studies which provided novel insight into the photochemistry of hydrogen-bonded complexes of the heptazine (Hz) chromophore with hydroxylic substrate molecules (water and phenol). It was found that a functionalized derivative of Hz, tri-anisole-heptazine (TAHz), can photooxidize water and phenol in a homogeneous photochemical reaction. This allows the exploration of the basic mechanisms of the proton-coupled electron-transfer (PCET) process involved in the water photooxidation reaction in well-defined complexes of chemically tunable molecular chromophores with chemically tunable substrate molecules. The unique properties of the excited electronic states of the Hz molecule and derivatives thereof are highlighted. The potential energy landscape relevant for the PCET reaction has been characterized by judicious computational studies. These data provided the basis for the demonstration of rational laser control of PCET reactions in TAHz-phenol complexes by pump-push-probe spectroscopy, which sheds light on the branching mechanisms occurring by the interaction of nonreactive locally excited states of the chromophore with reactive intermolecular charge-transfer states. Extrapolating from these results, we propose a general scenario which unravels the complex photoinduced water-splitting reaction into simple sequential light-driven one-electron redox reactions followed by simple dark radical-radical recombination reactions., Comment: 33 pages, 5 tables, 17 figures

Recent experiments have shown that photoexcited pyridine can abstract a hydrogen atom from a water molecule in pyridine-water clusters containing at least four water molecules. To explain these findings, we explored the electron-driven proton-transfer reaction from water to pyridine in pyridine-(H2O)n, n = 1–4, complexes with ab initio methods. It is shown that the photoinduced electron/proton transfer reaction is energetically possible for all clusters. The calculations reveal that the hydrogen bond between pyridine and the adjacent water molecule is weakened (strengthened) in the 1nπ∗ (1ππ∗) excited state, which is unfavorable (favorable) for the H-atom transfer reaction. For pyridine-(H2O)n clusters with n = 1–3, the steepest descent path leads to a local minimum of 1nπ∗ character, while for the pyridine-(H2O)4 cluster, this path leads to a local minimum of 1ππ∗ character. The transition state calculations show the presence of this 1ππ∗ minimum substantially reduces the barrier height for the H-atom transfer reaction. These results provide a tentative explanation of the experimental observations.

To gain mechanistic understanding of heptazine-based photochemistry, we synthesized and studied 2,5,8-tris(4-methoxyphenyl)-1,3,4,6,7,9,9b-heptaazaphenalene (TAHz), a model molecular photocatalyst chemically related to carbon nitride. On the basis of time-resolved photoluminescence (TR-PL) spectroscopy, we kinetically reveal a new feature that emerges in aqueous dispersions of TAHz. Using global target analysis, we spectrally and kinetically resolve the new emission feature to be blue shifted from the steady-state luminescence, and observe a fast decay component exhibiting a kinetic isotope effect (KIE) of 2.9 in H2O versus D2O, not observed in the steady-state PL. From ab initio electronic-structure calculations, we attribute this new PL peak to the fluorescence of an upper excited state of mixed nπ*/ππ* character. In water, the KIE suggests the excited state is quenched by proton-coupled electron transfer, liberating hydroxyl radicals that we detect using terephthalic acid. Our findings are consistent w...

The ultrafast transient dynamics of a free-base tetraphenylporphyrin (H2TPP) surface-mounted metal organic framework (SURMOF) is spectroscopically characterized and compared with that of H2TPP in e...

Polymeric carbon-nitride materials consisting of triazine or heptazine units have recently attracted vast interest as photocatalysts for water splitting with visible light. Adopting the hydrogen-bonded triazine-water complex as a model system, we explored the photochemical reaction mechanisms involved in the water splitting reaction in this system, using wavefunction-based ab initio electronic-structure methods. It is shown that photoexcited triazine can abstract a hydrogen atom from the water molecule by the sequential transfer of an electron and a proton from water to triazine, resulting in the triazinyl-hydroxyl biradical in the electronic ground state. It is furthermore shown that the excess hydrogen atom of the triazinyl radical can be photodetached by a second photon, which regenerates the triazine molecule. The hydrogen-bonded water molecule is thus decomposed into hydrogen and hydroxyl radicals in a biphotonic photochemical reaction. These results shed light on the molecular mechanisms of the water-oxidation reaction catalyzed by triazine-based organic polymers.

We present a theoretical description of double-pump femtosecond single-molecule signals with fluorescence detection. We simulate these signals in the weak-field regime for a model mimicking a chromophore with a Franck-Condon-active vibrational mode. We establish several signatures of these signals which are characteristic for the weak-field regime. The signatures include the quenching of vibrational beatings by electronic dephasing and a pronounced tilt of the phase-time profiles in the two-dimensional (2D) maps. We study how environment-induced slow modulations of the electronic dephasing and relevant chromophore parameters (electronic energy, orientation, vibrational frequency and relative shift of the potential energy surfaces) affect the signals.

The Renner vibronic-coupling problem in

The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.

According to Hund’s rule, the lowest triplet state (T1) is lower in energy than the lowest excited singlet state (S1) in closed-shell molecules. The exchange integral lowers the energy of the tripl...

We study the effect of a dissipative environment on the orientational relaxation of a three-dimensional quantum linear rotor. We provide a derivation of the Hamiltonian of a linear rotor coupled to a harmonic bath from first principles, confirming earlier conjectures. The dynamics generated by this Hamiltonian is investigated by the hierarchical equations-of-motion method assuming a Drude spectral density of the bath. We perform numerically accurate simulations and analyze the behavior of orientational correlation functions and the rotational structures of infrared absorption and Raman scattering spectra. We explore the features of orientational correlation functions and their spectra for a wide range of system-bath couplings, bath memory times, and temperatures. We discuss the signatures of the orientational relaxation in the underdamped regime, the strongly damped regime, and the librational regime. We show that the behavior of orientational correlation functions and their spectra can conveniently be analyzed in terms of three characteristic times, which are explicitly expressed in terms of the parameters of the Hamiltonian.

We introduce and analyze the concept of the "ideal" time and frequency resolved transient-absorption pump-probe (PP) signal. The ideal signal provides the most direct link between the "real" (measurable) PP signal and the material system dynamics. The simulation of PP signals involves two steps. (i) The ideal signal, which exhibits perfect time and frequency resolution, is calculated. For this purpose, the probe pulse is replaced by an auxiliary continuous-wave pulse. (ii) The real signal is obtained by the convolution of the ideal signal with the appropriate time- and frequency-gate function, which depends on the envelope of the actual probe pulse. This concept has been used to simulate integral and dispersed PP signals for a model system exhibiting nonadiabatic and dissipative dynamics. The ideal signal is computed with the two-pulse equation-of-motion phase-matching approach which has been extended to take excited-state absorption into account. We demonstrate how the ideal signal, an object exhibiting the features of moving wave packets as well as stationary spectra, is related to real signals detected with short (good temporal resolution) or long (good frequency resolution) probe pulses.

The reduction of carbon dioxide to fuels or chemical feedstocks with solar energy is one of the grand goals of current chemistry. In recent years, electrochemical, photoelectrochemical, and photochemical experiments provided hints of an unexpected catalytic role of the pyridine molecule in the reduction of carbon dioxide to formic acid or methanol. In particular, it has been suggested that the 1-pyridinyl radical (PyH) may be able to reduce carbon dioxide to the hydroxy-formyl radical. However, extensive theoretical studies of the thermodynamics and kinetics of the reaction called this interpretation of the experimental observations into question. Using ab initio computational methods, we investigated the photochemistry of the hydrogen-bonded PyH···CO2 complex. Our results reveal that carbon dioxide can be reduced to the hydroxy-formyl radical by a proton-coupled electron-transfer (PCET) reaction in excited states of the PyH···CO2 complex. In contrast to the ground-state PCET reaction, which exhibits a su...

It has been found in recent molecular beam experiments that pyridine molecules photoexcited at 255 nm can abstract hydrogen atoms from hydrogen-bonded water molecules in pyridine–water clusters, resulting in pyridinyl-hydroxyl radical pairs. The reaction could only be detected for clusters containing at least four water molecules. To provide insight into the mechanisms of this reaction, we performed ab initio excited-state trajectory surface-hopping dynamics simulations for two pyridine–water complexes, containing one and four water molecules, respectively, using the second-order algebraic-diagrammatic-construction (ADC(2)) electronic-structure method. A computationally efficient surface-hopping algorithm based on the Landau–Zener formula has been used to evaluate the transition probability between electronic states. The formation of the pyridinyl radical via an electron-driven proton transfer (EDPT) process from water to pyridine is confirmed by the simulations. The analysis of the competing excited-state reaction mechanisms up to 500 fs reveals that adiabatic relaxation to local minima of the S1(nπ*) potential-energy surface is the dominant channel in both clusters, followed by internal conversion to the electronic ground state via so-called ring-puckering conical intersections. The efficiency of the latter contribution is weakly dependent of the size of the clusters. The EDPT reaction occurs on the fastest time scales (faster than 200 fs) with a branching ratio of several percent. It is found to be four times more efficient in the pyridine–(H2O)4 cluster than in the pyridine–H2O cluster, which is qualitatively consistent with the experimental observations. A detailed understanding of the photoinduced reaction mechanisms in complexes of N-heterocyclic chromophores with water molecules is of relevance for future systematic knowledge-based developments of optimized materials for photocatalytic water splitting with sunlight.

The ultrafast nonadiabatic dynamics of a two-electronic-state four-vibrational-mode conical intersection coupled to a finite bath with up to 20 harmonic oscillators has been investigated by employing the multiple Davydov D

Unraveling the many facets of coherent and incoherent exciton motion in an ensemble of chromophores is an inherently complex quantum mechanical problem that has triggered a vivid debate on the role of quantum effects in molecular materials and biophysical systems. Here the dynamics of a statistical ensemble of molecular aggregates consisting of identical chromophores is investigated within a new theoretical framework. Taking account of intrinsic properties of the system, the Hamiltonian of the aggregate is partitioned into two mutually commuting vibrational and vibronic operators. This representation paves the way for an analysis that reveals the role of static disorder in ensembles of aggregates. Using analytical methods, it is demonstrated that after a critical time τD ≃ 2π/σ (σ being the dispersion of the disorder) any dynamic variable of the aggregate exhibits purely vibrational dynamics. This result is confirmed by exact numerical calculations of the time-dependent site populations of the aggregate. These findings may be useful for the interpretation of recent femtosecond spectroscopic experiments on molecular aggregates.

We present a first-principles study of the static and dynamic aspects of the strong Jahn–Teller (JT) and pseudo-JT (PJT) effects in niobium tetrafluoride, NbF4, in the manifold of its electronic ground state, 2E, and its first excited state, 2T2. The complex topography of the full-dimensional multi-sheeted adiabatic JT/PJT surfaces is analyzed computationally at the complete-active-space self-consistent-field (CASSCF) and multireference second-order perturbation levels of electronic structure theory, providing a detailed characterization of minima, saddle points, and minimum-energy conical intersection points. The calculations reveal that the tetrahedral (Td) configuration of NbF4 undergoes strong JT distortions along the bending mode of e symmetry, yielding tetragonal molecular structures of D2d symmetry with Td → D2d stabilization energies of about 2000 cm−1 in the X2E state and about 6400 cm−1 in the A2T2 state. In addition, there exists strong X2E−A2T2 PJT coupling via the bending mode of t2 symmetry, which becomes important near the crossing seam of the X2E and A2T2 potential energy surfaces. A five-state five-mode JT/PJT vibronic-coupling Hamiltonian is constructed in terms of symmetry-invariant polynomial expansions of the X2E and A2T2 diabatic potential energy surfaces in the e and t2 bending coordinates. The parameters of the Hamiltonian are determined by a least-squares fit of its eigenvalues to the CASSCF ab initio data. The vibronic spectra and the time evolution of adiabatic electronic population probabilities are computed with the multi-configuration time-dependent Hartree method. The complexity of the spectra reflects the effects of the exceptionally strong E × e and T2 × e JT couplings and (E + T2) × (e + t2) PJT coupling. The time evolution of the populations of the adiabatic electronic states after the initial preparation of the A2T2 state reveals the femtosecond nonadiabatic dynamics through a multidimensional seam of conical intersection. These results represent the first study of the static and dynamical JT/PJT effects in the X2E and A2T2 electronic states of NbF4.

Eumelanin is a naturally occurring skin pigment which is responsible for developing a suntan. The complex structure of eumelanin consists of π-stacked oligomers of various indole derivatives, such as the monomeric building block 5,6-dihydroxyindole (DHI). In this work, we present an ab initio wave-function study of the absorption behavior of DHI oligomers and of doubly and triply π-stacked species of these oligomers. We have simulated the onset of the electronic absorption spectra by employing the MP2 and the linear-response CC2 methods. Our results demonstrate the effect of an increasing degree of oligomerization of DHI and of an increasing degree of π-stacking of DHI oligomers on the onset of the absorption spectra and on the degree of red-shift toward the visible region of the spectrum. We find that π-stacking of DHI and its oligomers substantially red-shifts the onset of the absorption spectra. Our results also suggest that the optical properties of biological eumelanin cannot be simulated by considering the DHI building blocks alone, but instead the building blocks indole-semiquinone and indole-quinone have to be considered as well. This study contributes to advancing the understanding of the complex photophysics of the eumelanin biopolymer.

The reactivity of photoexcited 9H-adenine with hydrogen-bonded water molecules in the 9H-adenine-(H2 O)5 cluster is investigated by using ab initio electronic structure methods, focusing on the photoreactivity of the three basic sites of 9H-adenine. The energy profiles of excited-state reaction paths for electron/proton transfer from water to adenine are computed. For two of the three sites, a barrierless or nearly barrierless reaction path towards a low-lying S1 -S0 conical intersection is found. This reaction mechanism, which is specific for adenine in an aqueous environment, can explain the substantially shortened excited-state lifetime of 9H-adenine in water. Depending on the branching ratio of the nonadiabatic dynamics at the S1 -S0 conical intersection, the electron/proton transfer process can enhance the photostability of 9H-adenine in water or can lead to the generation of adenine-H(⋅) and OH(⋅) free radicals. Although the branching ratio is yet unknown, these findings indicate that adenine might have served as a catalyst for energy harvesting by water splitting in the early stages of the evolution of life.

International audience