Photo-oxidation of methanol in complexes with pyrido[2,3-b]pyrazine: a nonadiabatic molecular dynamics study

Excited-state Proton-Coupled Electron Transfer (PCET) constitutes a key step in the photo-oxidation of small, electron-rich systems possessing acidic hydrogen atoms, such as water or alcohols, which can play a vital role in green hydrogen production. In this contribution, we employ ab initio quantum-chemical methods and on-the-fly nonadiabatic molecular dynamics simulations to study the mechanism and the photodynamics of PCET in 1 : 1 pyrido[2,3-b]pyrazine complexes with methanol. We find the process to be ultrafast and efficient when the intramolecular hydrogen bond is formed with one of the β-positioned nitrogen atoms. The complex exhibiting a hydrogen bond with an isolated nitrogen site, on the contrary, shows much lower reactivity. We explain this effect with the stabilization of the reactive ππ* charge-transfer electronic state in the former case.

Surface hopping modeling of charge and energy transfer in active environments

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

Two-Sided Impact of Water on the Relaxation of Inner-Valence Vacancies of Biologically Relevant Molecules

After ionization of an inner-valence electron of molecules, the resulting cation-radicals store substantial internal energy which, if sufficient, can trigger ejection of an additional electron in an Auger decay usually followed by molecule fragmentation. In the environment, intermolecular Coulombic decay (ICD) and electron-transfer mediated decay (ETMD) are also operative, resulting in one or two electrons being ejected from a neighbor, thus preventing the fragmentation of the initially ionized molecule. These relaxation processes are investigated theoretically for prototypical heterocycle-water complexes of imidazole, pyrrole, and pyridine. It is found that the hydrogen-bonding site of the water molecule critically influences the nature and energetics of the electronic states involved, opening or closing certain relaxation processes of the inner-valence ionized system. Our results indicate that the relaxation mechanisms of biologically relevant systems with inner-valence vacancies on their carbon atoms can strongly depend on the presence of the electron-density donating or accepting neighbor, either water or another biomolecule.

Factors contributing to halogen bond strength and stretch or contraction of internal covalent bond

The halogen bond formed by a series of Lewis acids TF₃X (T = C, Si, Ge, Sn, Pb; X = Cl, Br, I) with NH₃ is studied by quantum chemical calculations. The interaction energy is closely mimicked by the depth of the σ-hole on the X atom as well as the full electrostatic energy. There is a first trend by which the hole is deepened if the T atom to which X is attached becomes more electron-withdrawing: C > Si > Ge > Sn > Pb. On the other hand, larger more polarizable T atoms are better able to transmit the electron-withdrawing power of the F substituents. The combination of these two opposing factors leaves PbF₃X forming the strongest XBs, followed by CF₃X, with SiF₃X engaging in the weakest bonds. The charge transfer from the NH₃ lone pair into the σ*(TX) antibonding orbital tends to elongate the covalent TX bond, and this force is largest for the heavier X and T atoms. On the other hand, the contraction of this bond deepens the σ-hole at the X atom, which would enhance both the electrostatic component and the full interaction energy. This bond-shortening effect is greatest for the lighter X atoms. The combination of these two opposing forces leaves the T-X bond contracting for X = Cl and Br, but lengthening for I.

Trajectory Surface Hopping for a Polarizable Embedding QM/MM Formulation

We present the implementation of trajectory surface-hopping nonadiabatic dynamics for a polarizable embedding QM/MM formulation. Time-dependent density functional theory was used at the quantum mechanical level of theory, whereas the molecular mechanics description involved the polarizable AMOEBA force field. This implementation has been obtained by integrating the surface-hopping program Newton-X NS with an interface between the Gaussian 16 and the Tinker suites of codes to calculate QM/AMOEBA energies and forces. The implementation has been tested on a photoinduced electron-driven proton-transfer reaction involving pyrimidine and a hydrogen-bonded water surrounded by a small cluster of water molecules and within a large water droplet.

Water Oxidation and Hydrogen Evolution with Organic Photooxidants: A Theoretical Perspective

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 Ångströ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-H₂O 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 ¹ππ* 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.

Suppression of Charge Recombination by Auxiliary Atoms in Photoinduced Charge Separation Dynamics with Mn Oxides: A Theoretical Study

Charge separation is one of the most crucial processes in photochemical dynamics of energy conversion, widely observed ranging from water splitting in photosystem II (PSII) of plants to photoinduced oxidation reduction processes. Several basic principles, with respect to charge separation, are known, each of which suffers inherent charge recombination channels that suppress the separation efficiency. We found a charge separation mechanism in the photoinduced excited-state proton transfer dynamics from Mn oxides to organic acceptors. This mechanism is referred to as coupled proton and electron wave-packet transfer (CPEWT), which is essentially a synchronous transfer of electron wave-packets and protons through mutually different spatial channels to separated destinations passing through nonadiabatic regions, such as conical intersections, and avoided crossings. CPEWT also applies to collision-induced ground-state water splitting dynamics catalyzed by Mn₄CaO₅ cluster. For the present photoinduced charge separation dynamics by Mn oxides, we identified a dynamical mechanism of charge recombination. It takes place by passing across nonadiabatic regions, which are different from those for charge separations and lead to the excited states of the initial state before photoabsorption. This article is an overview of our work on photoinduced charge separation and associated charge recombination with an additional study. After reviewing the basic mechanisms of charge separation and recombination, we herein studied substituent effects on the suppression of such charge recombination by doping auxiliary atoms. Our illustrative systems are X–Mn(OH)₂ tied to N-methylformamidine, with X=OH, Be(OH)₃, Mg(OH)₃, Ca(OH)₃, Sr(OH)₃ along with Al(OH)₄ and Zn(OH)₃. We found that the competence of suppression of charge recombination depends significantly on the substituents. The present study should serve as a useful guiding principle in designing the relevant photocatalysts.

Enhancement of tetrel bond involving tetrazole-TtR3 (Tt = C, Si; R = H, F). Promotion of SiR3 transfer by a triel bond

The combination of a CR3 (R = H, F) with a tetrazole can result in a moderate carbon bond, which can be further strengthened by a triel bond. On the other hand, SiR3 group is half transferred between the two N atoms in these conditions.

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-H₂O 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 HzH₂ in an exothermic reaction with a low reaction barrier. The time scales and branching ratios of competing channels in H-atom photodetachment from the HzH₂ 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.

Barrier-Lowering Effects of Baird Antiaromaticity in Photoinduced Proton-Coupled Electron Transfer (PCET) Reactions

Many popular organic chromophores that catalyze photoinduced proton-coupled electron transfer (PCET) reactions are aromatic in the ground state but become excited-state antiaromatic in the lowest ππ* state. We show that excited-state antiaromaticity makes electron transfer easier. Two representative photoinduced electron transfer processes are investigated: (1) the photolysis of phenol and (2) solar water splitting of a pyridine-water complex. In the selected reactions, the directions of electron transfer are opposite, but the net result is proton transfer following the direction of electron transfer. Nucleus-independent chemical shifts (NICS), ionization energies, electron affinities, and PCET energy profiles of selected [4n] and [4n + 2] π-systems are presented, and important mechanistic implications are discussed.

Unraveling the Mechanism for Tuning the Fluorescence of Fluorescein Derivatives: The Role of the Conical Intersection and nπ* State

Although a large number of fluorescein derivatives have been developed and applied in many different fields, the general mechanisms for tuning the fluorescence of fluorescein derivatives still remain uncovered. Herein, we found that the fluorescence quenching of neutral form of fluorescein derivatives in acidic medium resulted from a dark nπ* state, whereas the fluorescence of the anionic form of fluorescein derivatives in the gas phase and alkaline solutions was tuned by minimal energy conical intersection (MECI). The formation of MECI involved significant rotation of benzene ring and flip-flop motion of xanthene moiety, which would be restricted by intermolecular hydrogen bonding and lowering temperature. The energy barrier for reaching MECI depended on the substituents in the benzene moiety in accordance with experimentally observed substituent effects. These unprecedented mechanisms would lead to a recognition of fluorescein derivatives and could provide a correct and instructive design strategy for further developing new fluorescein derivatives.

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.

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 S₁ 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.

Performance of Mixed Quantum-Classical Approaches on Modeling the Crossover from Hopping to Bandlike Charge Transport in Organic Semiconductors

In the present study, several mixed quantum-classical (MQC) methods are applied to on-the-fly nonadiabatic molecular dynamics simulations of hole transport in molecular organic semiconductors (OSCs). The tested MQC methods contain the mean-field Ehrenfest (MFE), trajectory surface hopping (TSH) approaches based on Tully's fewest switches surface hopping (FSSH) and the global flux surface hopping (GFSH), the latter in the diabatic/adiabatic representation, and a Landau-Zener type trajectory surface hopping (LZSH). We also tested several correction schemes which were proposed to identify trivial crossings and to remove unphysical long-range charge transfers due to decoherence corrections. In addition, several cost-effective approaches for the nuclear velocity adjustment after an energy-allowed/energy-forbidden hop are investigated with respect to detailed balance and internal consistency conditions. To model a broad spectrum of OSCs with different charge transport characteristics, we derived from the anthracene structural model the construction of two additional models by uniformly scaling down the electronic couplings by the factors of 0.1 and 0.5. Anthracene shows a bandlike charge transport mechanism, characterized by slightly delocalized charge carriers 'diffusing' through the crystal. For smaller couplings, the mechanism changes to a hopping type, characteristically differing in the charge delocalization and temperature dependence. The MFE and corrected adiabatic TSH approaches are able to quantitatively reproduce the expected behavior, while the diabatic LZSH method fails for large couplings, as do approaches which are based on the hopping of localized charge between neighboring sites. Moreover, we find that while the hole mobility of the anthracene crystal simulated using the celebrated Marcus theory is in good agreement with the experimental value, its agreement has to be regarded as an accident due to the overestimation of the prefactor in the Marcus rate equation.

Charge separation and successive reconfigurations of electronic and protonic states in a water-splitting catalytic cycle with the Mn4CaO5 cluster. On the mechanism of water splitting in PSII

Charge separation, reloading of electrons and protons, and O₂ generation in a catalytic cycle for water splitting with Mn₄CaO₅ in PSII.

Photoinduced Water-Heptazine Electron-Driven Proton Transfer: Perspective for Water Splitting with g-C3N4

Heptazine-assembled polymeric carbon nitride (CN) materials have fascinated the research community as a photocatalyst for hydrogen evolution while less attention has been devoted to the mechanistic features of the host materials. Using excited-state nonadiabatic dynamics simulations, the molecular-level picture of the decomposition of heptazine hydrogen bonded to water molecule(s) (heptazine-water complex) into heptazinyl and hydroxyl biradical products is revealed. Dynamics simulations show that hydrogen detachment from the water molecule to the heptazine occurs within tens of femtoseconds and suggest that excited-state deactivation via N-H······O-H electron-driven proton transfer (EDPT) is the dominant and most relevant excited-state deactivation process in heptazine-water complexes leading to conical intersection. The observation of photorelaxation-induced water splitting by heptazine is proof of the water-splitting reaction principle, which presents further challenges for computational and experimental investigations of the deactivation of heptazinyl and OH biradical products for efficient hydrogen evolution.

Mechanism of Photoinduced Dihydroazulene Ring-Opening Reaction

The photoinduced ring-opening reaction is a key process in the functioning of dihydroazulene/vinylheptafulvene (DHA/VHF) photoswitches. Over the years, the mechanism of this reaction has been extensively debated. Herein, by means of nonadiabatic trajectory dynamics simulations and quantum chemistry calculations, we present the first detailed and comprehensive investigation on the mechanism of the photoinduced ring-opening reaction of DHA. The results show the crucial role of the excited-state ring planarization process for the bond breaking. Our dynamics simulations show that the DHA ring opening is an ultrafast reaction that does not follow exponential kinetics but exhibits ballistic dynamics. Upon photoexcitation, the planarization occurs within 300-500 fs. This leads to the ring-opening reaction and concurrent decay of the molecule to the ground state within 100 fs through an S₁ → S₀ internal conversion process toward forming the VHF isomer. These results are consistent with previous ultrafast time-resolved experiments and lead to a thorough understanding of the DHA/VHF photoconversion.