Role of water in intramolecular proton transfer reactions of formamide and thioformamide

Ultrafast excited state dynamics of pyridine N-oxide derivative in solution; femtosecond fluorescence up-conversion and theoretical calculations

In this study we have investigated 2-ethylamino-4-nitro-6-methyl pyridine N-oxide (2E6M) molecule that belongs to important group of Proton Coupled Electron Transfer (PCET) compounds where both the charge transfer (CT) and proton transfer processes in excited states may proceed. In this case, this is possible due to the donors and acceptors of electrons and protons in this system, as well as due to the presence of intramolecular {N-H… O [2,566(3) Å}, hydrogen bond.Using stationary and time-resolved spectroscopy, as well as quantum chemical calculations on the DFT and TD DFT B3LYP/6-31G (d,p) level of theory, a partial CT nature of the S₀ → S₁ transition in both tautomeric forms (N and T) has been revealed. Additionally, the excited state intramolecular proton transfer (ESIPT) process shown to be more favorable in apolar and weakly polar solvents than in strongly polar acetonitrile (E_N(S1) > E_T(S1). The displacement of charge from the amine group and the ring to the nitro group has been observed on the changing shapes of the HOMO and LUMO orbitals involved in this transition what further quantitatively allowed to realize the increase in the dipole moment of both forms in the electronic excited state. The calculations show that in two solvents with radically different polarity (heptane, acetonitrile), dipole moments of both forms are very similar [in acetonitrile u_N(S₁) and u_T(S₁) are 11.0 D and 11.5 D, respectively]. Hence, in polar media both forms can be stabilized in a comparable manner. This made it difficult for us to assign a single fluorescent band in acetonitrile to one of the tautomeric forms. However, it seems that due to application of time-resolved spectroscopy, this problem has been clarified. The TCSPC decay curve in acetonitrile with an ultrafast lifetime assigned to the (N) form, along with the femtosecond up-conversion signals that demonstrated only an ultrafast decay without any rise-time of a new excited (T) species, allowed us to conclude that in 2E6M in strongly polar solvent the ESIPT does not occur.The unique fluorescence band origins from the (N) form. In protic solvents, the significant kinetic isotopic effects have provided us with conclusive evidence for the presence of the solvent-assisted ESIPT process. Furthermore, it was noticed that the fluorescence lifetime in D₂O (100-120 fs) estimated from the up-conversion signals is about 40 times shorter relative to methanol. This may suggest that the sine qua non for the ESIPT process in 2E6M in protic solvents is the formation of a complex with a solvent molecule in the hydrogen bridge between the proton donor and proton acceptor, respectively.

Mechanistic insight into the tautomerization of histidine initiated by water-catalyzed N-H and C-H cleavages

The N-H and C-H activation is of great significance in organic chemistry and chemical industry fields, especially, in the utilization of petroleum raw materials. High NδH (tautomer of natural histidine) content would increase Alzheimer's disease risk. To inhibit this and improve the activation of N-H and C-H bonds, the isomerization mechanism from NδH to NεH of histidine-containing dipeptide catalyzed by water cluster was explored. The results discovered that water cluster assists this reaction by reducing the activation energies from 68.20 to 9.60 kcal mol^-1, and its size not only affects the reaction rate but also determines the reaction pathway in a degree. Moreover, water cluster, taken as a potential green catalyst, is more effective on the reactions involving N-H and C-H bond cleavages than reported common toxic organometallic compounds and has different catalytic mechanisms. This work also provides some theoretical guidance for the modulation of Alzheimer's disease induced by histidine isomerization.

Improved algorithm for the direct dynamics variational multi-configurational Gaussian method

The Direct Dynamics variational Multi-Configurational Gaussian (DD-vMCG) method provides a fully quantum mechanical solution to the time-dependent Schrödinger equation for the time evolution of nuclei with potential surfaces calculated on-the-fly using a quantum chemistry program. Initial studies have shown its potential for flexible and accurate simulations of non-adiabatic excited-state molecular dynamics. In this paper, we present developments to the DD-vMCG algorithm that improve both its accuracy and efficiency. First, a new, efficient parallel algorithm to control the DD-vMCG database of quantum chemistry points is presented along with improvements to the Shepard interpolation scheme. Second, the use of symmetry in describing the potential surfaces is introduced along with a new phase convention in the propagation diabatization. Benchmark calculations on the allene radical cation including all degrees of freedom then show that the new scheme is able to produce a consistent non-adiabatic coupling vector field. This new DD-vMCG version thus opens the route for effectively and accurately treating complex chemical systems using quantum dynamics simulations.

Proton transfer in the benzimidazolone and benzimidazolthione tautomerism process catalyzed by polar protic solvents

The tautomeric equilibrium of benzimidazolone and benzimidazolthione have been studied by the density functional theory method using the CAM-B3LYP functional together with the 6-311G(d,p) basis set. Two separate mechanisms have been investigated: a direct intramolecular transfer using the polarizable continuum model and an indirect proton transfer assisted by a molecule of solvent (C6H12, H2O, CH3OH, and H2O2). In both cases, the results obtained indicate that ketone and thione are the most stable forms. However, the enhanced height of the activation barrier for the four-center mechanisms describing the tautomerism reaction as a direct intramolecular transfer implicates a relatively disadvantaged process. The participation of a polar protic solvent molecule allows the lowering of the activation energy barrier. Potential energy profiles of keto-enol and thio-enol tautomerism assisted by methanol and water are very different. The former one describes a concerted mechanism but the latter does not because it is associated with asynchronous processes that take place during the thio-enol tautomerism.

Hydrogenation and hydration of carbon dioxide: a detailed characterization of the reaction mechanisms based on the reaction force and reaction electronic flux analyses

A computational DFT study of the reaction mechanism of hydrogenation and hydration of carbon dioxide is presented. It has been found that hydrogenation and hydration are endoenergetic reactions that are carried out in two steps, passing by a stable intermediate that is surrounded by energy barriers of 70 kcal/mol and 10 kcal/mol for hydrogenation and 50 kcal/mol and 10 kcal/mol for hydration. Using the reaction force analysis, we were able to characterize the physical nature of the activation barriers and found that activation energies are mostly due to structural rearrangements. An interesting difference in the reaction mechanisms disclosed by the reaction force and electronic flux analyses is that while in the hydrogenation reaction the mechanisms is conditioned by the H₂ cleavage with a high energy barrier, in the hydration reaction the formation of a transient four member ring structure favoring an attractive local hydrogen bond interaction pushes the reaction toward the product with a considerably lower energy barrier.

Physico-chemical profiles of the wobble ↔ Watson-Crick G*·2AP(w) ↔ G·2AP(WC) and A·2AP(w) ↔ A*·2AP(WC) tautomerisations: a QM/QTAIM comprehensive survey

This study is intended to clarify in detail the tautomeric transformations of the wobble (w) G*·2AP(w) and A·2AP(w) nucleobase mispairs involving 2-aminopurine (2AP) into the Watson-Crick (WC) G·2AP(WC) and A*·2AP(WC) base mispairs (asterisks denote mutagenic tautomers of the DNA bases), respectively, by quantum-mechanical methods and Bader's Quantum Theory of Atoms in Molecules. Our previously reported methodology has been used, which allows the evolution of the physico-chemical parameters to be tracked along the entire internal reaction coordinate (IRC), not exclusively in the stationary states of these reactions. These biologically important G*·2AP(w) ↔ G·2AP(WC) and A·2AP(w) ↔ A*·2AP(WC) w ↔ WC tautomerisations, which are involved in mutagenic tautomerically-conformational pathways, determine the origin of the transitions and transversions induced by 2AP. In addition, it is established that they proceed through planar, highly stable, zwitterionic transition states and they exhibit similar physico-chemical profiles and stages of sequential intrapair proton transfer, followed by spatial rearrangement of the nucleobases relative to each other within the base pairs. These w ↔ WC tautomerisations occur non-dissociatively and are accompanied by a significant alteration in geometry (from wobble to Watson-Crick and vice versa) and redistribution of the specific intermolecular interactions, which can be divided into 10 patterns including AHB H-bonds and loosened A-H-B covalent bridges along the IRC of tautomerisation. Based on the redistribution of the geometrical and electron-topological parameters of the intrapair hydrogen bonds, exactly 9 key points have been allocated to characterize the evolution of these reactions.

Curve crossing in a manifold of coupled electronic states: direct quantum dynamics simulations of formamide

Quantum dynamics simulations are an important tool to evaluate molecular behaviour including the, often key, quantum nature of the system. In this paper we present an algorithm that is able to simulate the time evolution of a molecule after photo-excitation into a manifold of states. The direct dynamics variational multi-configurational Gaussian (DD-vMCG) method circumvents the computational bottleneck problems of traditional grid-based methods by computing the potential energy functions on-the-fly, i.e. only where required. Unlike other commonly used direct dynamics methods, DD-vMCG is fully quantum mechanical. Here, the method is combined with a novel on-the-fly diabatisation scheme to simulate the short-time dynamics of the key molecule formamide and its acid analogue formimidic acid. This is a challenging test system due to the nature and large number of excited states, and eight coupled states are included in the calculations. It is shown that the method is able to provide unbiased information on the product channels open after excitation at different energies and demonstrates the potential to be a practical scheme, limited mainly by the quality of the quantum chemistry used to describe the excited states.

A QM/QTAIM detailed look at the Watson-Crick↔wobble tautomeric transformations of the 2-aminopurine·pyrimidine mispairs

This work is devoted to the careful QM/QTAIM analysis of the evolution of the basic physico-chemical parameters along the intrinsic reaction coordinate (IRC) of the biologically important 2AP·T(WC)↔2AP·T*(w) and 2AP·C*(WC)↔2AP·C(w) Watson-Crick(WC)↔wobble(w) tautomeric transformations obtained at each point of the IRC using original authors' methodology. Established profiles reflect the high similarity between the courses of these processes. Basing on the scrupulous analysis of the profiles of their geometric and electron-topological parameters, it was established that the dipole-active WC↔w tautomerizations of the Watson-Crick-like 2AP·T(WC)/2AP·C*(WC) mispairs, stabilized by the two classical N3H⋯N1, N2H⋯O2 and one weak C6H⋯O4/N4 H-bonds, into the wobble 2AP·T*(w)/2AP·C(w) base pairs, respectively, joined by the two classical N2H⋯N3 and O4/N4H⋯N1 H-bonds, proceed via the concerted stepwise mechanism through the sequential intrapair proton transfer and subsequent large-scale shifting of the bases relative each other, through the planar, highly stable, zwitterionic transition states stabilized by the participation of the four H-bonds - N1⁺H⋯O4^-/N4^-, N1⁺H⋯N3^-, N2⁺H⋯N3^-, and N2⁺H⋯O2^-. Moreover, it was found out that the 2AP·T(WC)↔2AP·T*(w)/2AP·C*(WC)↔2AP·C(w) tautomerization reactions occur non-dissociatively and are accompanied by the consequent replacement of the 10 unique patterns of the specific intermolecular interactions along the IRC. Obtained data are of paramount importance in view of their possible application for the control and management of the proton transfer, e.g. by external electric or laser fields.

A QM/QTAIM research under the magnifying glass of the DPT tautomerisation of the wobble mispairs involving 2-aminopurine

In this study, a comprehensive survey of the changes of the physico-chemical parameters at each point of the IRC for the biologically important T·2AP*(w) ↔ T*·2AP(w) and G·2AP*(w) ↔ G*·2AP(w) DPT tautomerisation reactions involved in the point mutations (transitions and transversions) induced by 2-aminopurine (2AP) in DNA is provided.

Study of antiradical mechanisms with dihydroxybenzenes using reaction force and reaction electronic flux

Phenolic compounds represent an important category of antioxidants, DFT reactivity indexes proved that the substitution of electron donating groups in dihydroxybenzenes, favor the spontaneity of the hydrogen atom transfer mechanism.