Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals.

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Excited-state singlet-triplet inversion in hexagonal aromatic and heteroaromatic compounds.

The inversion of the energies of the lowest singlet (S1) and lowest 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, wave-function based ab initio computational methods were employed to explore the possibility of S1-T1 inversion in hexagonal polycyclic aromatic and heteroaromatic compounds. In these molecules, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) are two-fold degenerate. The HOMO-LUMO transition gives rise to three singlet and three triplet excited states. While the singlet-triplet energy gap ΔST, defined as the energy difference between the S1 state and the T1 state, is clearly positive for benzene, it is predicted to be close to zero for borazine, the boron nitride analogue of benzene. Although ΔST decreases with increasing size of hexagonal polycyclic aromatics, it remains positive up to circumcoronene (19 rings). However, symmetry-preserving 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. These findings establish the existence of a new family of boron carbon nitrides with inverted singlet-triplet gaps.

Reduction of CO2 with Hydrated Electrons: Ab Initio Computational Studies for Finite-Size Cluster Models.

In the present work, the mechanisms of the reduction of the CO2 molecule with hydrated electrons to the hydroxyl-formyl (HOCO) radical were studied with ab initio computational methods. Hydrated hydronium radicals, H3O(H2O)n (n = 0,3,6), are considered as finite-size models of the hydrated electron in liquid water. The investigation of cluster models allows the application of high-accuracy electronic-structure methods, which are not computationally feasible in condensed-phase simulations. Reaction paths and potential-energy (PE) profiles of the proton-coupled electron-transfer reaction from hydrated H3O radicals to the CO2 molecule were explored on the ground-state PE surface. The computationally efficient unrestricted second-order Møller-Plesset method is employed, and its accuracy has been carefully benchmarked in comparison with complete-active-space self-consistent-field and multi-reference second-order perturbation calculations. The results provide insights into the interplay of electron transfer from the diffuse Rydberg-type unpaired electron of H3O to the CO2 molecule, the contraction of the electron cloud by the re-hybridization of the carbon atom of CO2, and proton transfer from the nearest water molecule to the CO2- anion, followed by Grotthus-type proton rearrangements to form stable clusters. Starting from local energy minima of hydrogen-bonded CO2-H3O(H2O)n complexes, the reaction to form HOCO-(H2O)n+1 complexes is exothermic by about 1.3 eV (125 kJ/mol). The reaction is barrier controlled with a barrier of the order of a few tenths of an electron volt, depending on size and conformation of the water cluster. This barrier is at least an order of magnitude lower than the barrier of the reaction of CO2 with any closed-shell partner molecule. The HOCO radicals can recombine by H-atom transfer (disproportionation), resulting in formic acid or a dihydroxycarbene product, as well as by the formation of a C-C bond, resulting in oxalic acid. The strong exothermicity of these radical-radical recombination reactions likely results in the fragmentation of the closed-shell products formic acid and oxalic acid, which explains the strong specificity for CO formation observed in recent experiments of Hamers and co-workers.

Ab initio trajectory surface-hopping dynamics studies of excited-state proton-coupled electron transfer reactions in trianisoleheptazine-phenol complexes.

The excited-state proton-coupled electron-transfer (PCET) reaction in hydrogen-bonded complexes of trianisoleheptazine (TAHz), a chromophore related to polymeric carbon nitrides widely used in hydrogen-evolution photocatalysis, with several phenol derivatives were recently studied by Schlenker and coworkers with time-resolved photoluminescence quenching and pump-probe experiments. A pronounced dependence of the PCET reactivity on the electron-donating/electron-withdrawing character of the substituents on phenol was found, with indications of a barrierless or nearly barrierless PCET reaction for the most strongly electron-donating substituent, methoxy. In the present work, the excited-state PCET dynamics was explored with first-principles nonadiabatic dynamics simulations using the TDDFT/ωB97X-D electronic-structure model for two selected complexes, TAHz-phenol and TAHz-methoxyphenol. The qualitative reliability of the TDDFT/ωB97X-D electronic-structure model was assessed by extensive benchmarking of excitation energies and potential-energy profiles against a wave-function-based ab initio method, the algebraic-diagrammatic construction of second order (ADC(2)). The nonadiabatic dynamics simulations provide temporally and structurally resolved insights into paradigmatic PCET reactions in TAHz-phenol complexes. The radiationless relaxation of the photoexcited bright 1ππ* state to the long-lived dark S1 state of TAHz occurs in less than 100 fs. The ensuing PCET reaction on the adiabatic S1 surface is faster in TAHz-methoxyphenol complexes than in TAHz-phenol complexes due to a lower H-atom-transfer barrier, as observed in the experiments. The relaxation of the complexes to the electronic ground state is found to occur exclusively via PCET within the 250 fs time window covered by the present simulations, confirming the essential role of the hydrogen bond for the fluorescence quenching process. The absolute values of the computed PCET time constants are significantly shorter than those extracted from time-resolved photoluminescence measurements for mixtures of TAHz with phenolic substrates in toluene. The possible origins of this discrepancy are discussed.

On the propensity of formation of cyclobutane dimers in face-to-face and face-to-back uracil stacks in solution.

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.

Ab Initio Electronic Structure Study of the Photoinduced Reduction of Carbon Dioxide with the Heptazinyl Radical.

The photocatalytic conversion of carbon dioxide to liquid fuels with electrons taken from water with solar photons is one of the grand goals of renewable energy research. Polymeric carbon nitrides recently emerged as metal-free materials with promising functionalities for hydrogen evolution from water as well as the activation of carbon dioxide. Molecular heptazine (Hz), the building block of polymeric carbon nitrides, is one the strongest known organic photo-oxidants and has been shown to be able to photo-oxidize water with near-visible light, resulting in reduced (hydrogenated) heptazine (HzH) and OH radicals. In the present work, we explored with ab initio computational methods whether the HzH chromophore is able to reduce carbon dioxide to the hydroxy-formyl (HOCO) radical in hydrogen-bonded HzH-CO2 complexes by the absorption of a photon. In remarkable contrast to the high barrier for carbon dioxide activation in the electronic ground state, the excited-state proton-coupled electron transfer (PCET) reaction is nearly barrierless, but requires the diabatic passage of three conical intersections. The possibility of barrierless carbon dioxide activation by excited-state PCET has so far not been taken into consideration in the interpretation of photocatalytic carbon dioxide reduction on carbon nitride materials.

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.

Triangular boron carbon nitrides: an unexplored family of chromophores with unique properties for photocatalysis and optoelectronics.

It has recently been shown that cycl[3.3.3]azine and heptazine (1,3,4,6,7,9,9b-heptaazaphenalene) as well as related azaphenalenes exhibit inverted singlet and triplet states, that is, the energy of the lowest singlet excited state (S1) is below the energy of the lowest triplet excited state (T1). This feature is unique among all known aromatic chromophores and is of outstanding relevance for applications in photocatalysis and organic optoelectronics. Heptazine is the building block of the polymeric material graphitic carbon nitride which is an extensively explored photocatalyst in hydrogen evolution photocatalysis. Derivatives of heptazine have also been identified as efficient emitters in organic light emitting diodes (OLEDs). In both areas, the inverted singlet-triplet gap of heptazine is a highly beneficial feature. In photocatalysis, the absence of a long-lived triplet state eliminates the activation of atmospheric oxygen, which is favourable for long-term operational stability. In optoelectronics, singlet-triplet inversion implies the possibility of 100% fluorescence efficiency of electron-hole recombination. However, the absorption and luminescence wavelengths of heptazine and the S1-S0 transition dipole moment are difficult to tune for optimal functionality. In this work, we employed high-level ab initio electronic structure theory to devise and characterize a large family of novel heteroaromatic chromophores, the triangular boron carbon nitrides. These novel heterocycles inherit essential spectroscopic features from heptazine, in particular the inverted singlet-triplet gap, while their absorption and luminescence spectra and transition dipole moments are widely tuneable. For applications in photocatalysis, the wavelength of the absorption maximum can be tuned to improve the overlap with the solar spectrum at the surface of earth. For applications in OLEDs, the colour of emission can be adjusted and the fluorescence yield can be enhanced.

Ab Initio Nonadiabatic Surface-Hopping Trajectory Simulations of Photocatalytic Water Oxidation and Hydrogen Evolution with the Heptazine Chromophore.

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.

Strong static and dynamic Jahn-Teller and pseudo-Jahn-Teller effects in niobium tetrafluoride.

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 X̃2E state and about 6400 cm-1 in the Ã2T2 state. In addition, there exists strong X̃2E-Ã2T2 PJT coupling via the bending mode of t2 symmetry, which becomes important near the crossing seam of the X̃2E and Ã2T2 potential energy surfaces. A five-state five-mode JT/PJT vibronic-coupling Hamiltonian is constructed in terms of symmetry-invariant polynomial expansions of the X̃2E and Ã2T2 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 Ã2T2 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 X̃2E and Ã2T2 electronic states of NbF4.

Effects of high pulse intensity and chirp in two-dimensional electronic spectroscopy of an atomic vapor.

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.

Photoinduced water oxidation in pyrimidine-water clusters: a combined experimental and theoretical study.

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.

Molecular Design of Heptazine-Based Photocatalysts: Effect of Substituents on Photocatalytic Efficiency and Photostability.

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.

Multi-faceted spectroscopic mapping of ultrafast nonadiabatic dynamics near conical intersections: A computational study.

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.

Photooxidation of water with heptazine-based molecular photocatalysts: Insights from spectroscopy and computational chemistry.

We present a conspectus of recent joint spectroscopic and computational studies that 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 that unravels the complex photoinduced water-splitting reaction into simple sequential light-driven one-electron redox reactions followed by simple dark radical-radical recombination reactions.

UV absorption spectra of DNA bases in the 350-190 nm range: assignment and state specific analysis of solvation effects.

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.

Photoinduced electron-driven proton transfer from water to an N-heterocyclic chromophore: nonadiabatic dynamics studies for pyridine-water clusters.

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.

Mechanisms of photoreactivity in hydrogen-bonded adenine-H2O complexes.

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.

Role of the Pyridinyl Radical in the Light-Driven Reduction of Carbon Dioxide: A First-Principles Study.

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 substantial barrier, the excited-state PCET reaction is barrierless but requires the passage through two conical intersections. Our results provide a tentative explanation of the catalytic role of the PyH radical in the reduction of CO2 with the qualification that the absorption of a photon by PyH is necessary.

Singlet-Triplet Inversion in Heptazine and in Polymeric Carbon Nitrides.

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 triplet state and raises the energy of the singlet state of the same orbital character, leading to a positive singlet-triplet energy gap (ΔST). Exceptions are known for biradicals and charge-transfer excited states of large molecules in which the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) are spatially separated, resulting in a small exchange integral. In the present work, we discovered with ADC(2), CC2, EOM-CCSD, and CASPT2 calculations that heptazine (1,3,4,6,7,9,9b-heptaazaphenalene or tri-s-triazine) exhibits an inverted S1/T1 energy gap (ΔST ≈ -0.25 eV). This appears to be the first example of a stable closed-shell organic molecule exhibiting S1/T1 inversion at its equilibrium geometry. The origins of this phenomenon are the nearly pure HOMO-LUMO excitation character of the S1 and T1 states and the lack of spatial overlap of HOMO and LUMO due to a unique structure of these orbitals of heptazine. The S1/T1 inversion is found to be extremely robust, being affected neither by substitution of heptazine nor by oligomerization of heptazine units. Using time-resolved photoluminescence and transient absorption spectroscopy, we investigated the excited-state dynamics of 2,5,8-tris(4-methoxyphenyl)-1,3,4,6,7,9,9b-heptaazaphenalene (TAHz), a chemically stable heptazine derivative, in the presence of external heavy atom sources as well as triplet-quenching oxygen. These spectroscopic data are consistent with TAHz singlet excited state decay in the absence of a low-energy triplet loss channel. The absence of intersystem crossing and an exceptionally low radiative rate result in unusually long S1 lifetimes (of the order of hundreds of nanoseconds in nonaqueous solvents). These features of the heptazine chromophore have profound implications for organic optoelectronics as well as for water-splitting photocatalysis with heptazine-based polymers (e.g., graphitic carbon nitride) which have yet to be systematically explored and exploited.