Theoretical study of intermolecular proton transfer reaction in isolated 5-hydroxyisoxazole–water complexes

A systematic investigation in isolated 5-hydroxyisoxazole–water complexes (5-HIO · (H₂O)_nn = 1–3) is performed at the DFT level, employing B3LYP/6-31G(d, p) basis set. Single-point energy calculations are also performed at the MP2 level using B3LYP/6-31G(d, p) optimized geometries and the 6-311++G(d, p) basis set. The computational results show that the keto tautomer K₂ is the most stable isomer in the gas phase, and the tautomer K₁ to be the next most stable tautomer. Hydrogen bonding between HIO and the water molecule(s) will dramatically lower the barrier by a concerted multiple proton transfer mechanism. The proton transfer process of 3WE_cis ↔ 3WK₁ and 2WE_trans ↔ 2WK₂ is found to be more efficient in two tautomerization, and the barrier heights are 7.03 and 14.15 kcal/mol at B3LYP/6-31G(d, p) level, respectively. However, the proton transfer reaction between E_cis and K₁ cannot happen without solvent-assisted.     

Theoretical study on molecular structures, intramolecular proton transfer reaction, and solvent effects of 1-phenyl-3-methyl-4-(6-hydro-4-amino-5-sulfo-2,3-pyrazine)-pyrazole-5-one

The molecular structures and intramolecular proton transfer reaction of 1-phenyl-3-methyl-4-(6-hydro-4-amino-5-sulfo-2,3-pyrazine)-pyrazole-5-one have been investigated with both ab initio method and the density functional theory. The solvent effects are simulated using the self-consistent reaction field (SCRF) method within the framework of the polarizable continuum model (PCM). The results show that the computed geometrical parameters at the B3LYP levels are in better agreement with experimental values than those at the RHF levels, and the choice of functional in DFT plays an important role in describing the title compound. It is found that strong hydrogen bonds (O–H···N and O···H–N) exist in the title compound, and in the proton transfer process, the O–H bond is broken while the N–H bond is formed. In addition, the order of stability of the isomers remains the same in different solvents, while the barrier height of the proton transfer reaction and dipole moments for the title compound grow with the increase of the solvent polarity. Eventually, the NBO analysis shows that the strength of the hydrogen bond reduces with the increase of the solvent polarity.     

水催化2个酯分子相互转化反应的理论研究

从印楝植物内生真菌Phomopsis sp.培养液中分离得到的4-acetoxymultiplolide(1)和1-acetoxymultiplo-lide(2)在室温及水存在下能够相互转化. 提出二者相互转化最可能的4个途径(机理A~D). 在B3LYP/6-311+G(d,p)水平进行气相条件的优化, 结果表明, 无水催化的机理A中TS1和TS2的活化能均显著大于120 kJ/mol, 2个分子水催化的机理D中TS1和TS2的活化能则显著降低. 计算结果显示水的溶剂化效应能进一步降低机理D中TS1和TS2的活化能. 在MP2/6-311++G(2d,2p)//B3LYP/6-311+G(d,p)水平计算了单点能, 得到在水相时机理D中TS1和TS2的活化能分别为106.24和107.37 kJ/mol. 因此, 机理D是化合物1 和2在室温下及水存在时相互转化最可能的途径, 该途径是一种特殊的水催化分子内酯的醇解反应, 也是一种经典的亲核加成反应, 通过一种新的叔醇中间体实现.     

A CASSCF-CASPT2 study of the excited-state intramolecular proton transfer reaction in 1-amino-3-propenal using different active spaces

In this work we analyze how the choice of the active space in the CASSCF (the complete-active-space multiconfiguration self-consistent-field method) and CASPT2 (the second-order perturbation theory based on the CASSCF reference wave function) calculations affects the computed potential energy curves (PECs) for the intramolecular proton transfer reaction in the ground state and the two lowest lying singlet excited states of 1-amino-3-propenal. As anticipated, the results revealed that, qualitatively, the proton transfer in the different states can be correctly described even by minimal active spaces, which include the orbitals involved in the electronic excitation of the considered state and the antibonding sigma orbital corresponding to the bond formed by the molecule with the migrating hydrogen atom. However, quantitatively, the relative energies of the two tautomers and the energy barriers computed at the CASSCF level change when the active space is increased, indicating importance of the dynamic electron correlation. Introducing the dynamic correlation effects via CASPT2 makes the calculated energy parameters more uniform among the different active spaces. The analysis suggested certain optimal active spaces for studying proton transfer reactions in systems similar to 1-amino-3-propenal. The PEC calculations for excited states showed that the results are sensitive to the molecular geometries used in the calculations, particularly near the transition point. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1422–1431 (1999)     

AMT异构体互变机理的理论研究

用密度泛函理论(DFT)的B3LYP/6-311G(d, p)和Müller-Plesset微扰理论的MP2/6-31G(d)方法，优化了AMT(2-氨基-5-巯基-1,3,4-噻二唑)各种异构体和过渡态结构的几何构型，并对它们的电子结构、振动光谱和化学键性质进行了研究.还研究了AMT异构体的互变机理，提出了AMT异构体abcda的循环式互变途径.进一步完成了对AMT异构体成键方式的自然键轨道(NBO)分析.     

Theoretical study of the reaction mechanism of proton transfer in glycinamide

To investigate the tautomerism of glycinamide that is induced by proton transfer, we present detailed theoretical studies on the reaction mechanism of both the isolated gas phase and H₂O‐assisted proton transfer process of glycinamide, using density functional theory calculations by means of the B3LYP hybrid functional. Twenty‐six geometries, including 10 significant transition states, were optimized, and these geometrical parameters are discussed in detail. The relative order of the activation energy for hydrogen atom transfer of all the conformers has been systematically explored in this essay. For the amido hydrogen atom transfer process, the relative order of the activation energy is: IV < II < III < I, while in the carbonic hydrogen atom transfer process, the relative order is IV > II > III > I. Meanwhile, the most favorable structure for both the amido hydrogen atom transfer and the carbonic hydrogen atom transfer has been found. The involvement of the water molecule not only can stabilize the transition states and the ground states, but can also reduce the activation energy greatly. The superior catalytic effect of H₂O has been discussed in detail. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Theoretical study on proton tautomerism of chryosphanol and its analogues in the ground and excited states

Proton tautomerism of 1,8-dihydroxy-3-methyl-anthraquinone and its analogues were studied using HF and CIS methods with 6-31g(d,p) basis set for the ground and singlet excited states. The calculations indicate that the compound exists two strong intramolecular hydrogen bonds (IHB), and shows similar characters in its proton transfer processes considering the geometries and Mulliken charge population. Calculation results further show that intramolecular proton transfer (IPT) is not favored in view of the energy trend for chryosphanol, which has two hydrogens of hydroxyl groups bond with a common oxygen of carbonyl group and exists two IHBs in the peri region. However, it exhibits normal intramolecular proton transfer for the derivatives of chryosphanol, which have only one pair of adjacent hydroxyl group and carbonyl group existing in the peri region. Hereby, it can be conjectured from a theoretical point of view that IPT is absent in the similar structure such as hypericin’s peri region. Calculation results on the photophysical process show that the isomerization process is competitive with the intersystem crossing process, which facilitates the increase of triplet state quantum efficiency and photosensitive activity.     

尿酸质子转移反应的理论研究I——分子内质子转移研究

使用量子化学中的Hartree-Fock方法和密度泛函理论中的B3LYP方法，分别在3－21G^*和6－31G（d)水平上，计算了尿酸分子从三羰基异构体向三羟基异构体的转化。结果表明，转化过程经历了单羟基和双羟基异构体2种中间物和3种过渡态时的分子内质子转移（IPT），转移中的H原邻近的N，O和C原子形成了具有四元环结构的过渡态。随着IPT的进行，N－H键逐渐被削弱和断裂，O－H键则逐渐生成。3个反应的活化能分别为190.3kJ/mol,181.4kJ/mol和249.9kJ/mol(B3LYP/6-31G(d))。较高的活化能表明在室温下，无催化剂的IPT难以进行。     

Theoretical studies of conjugation and substituent effect on intramolecular proton transfer in the ground and excited states

The ground- and excited-state intramolecular proton transfer (GSIPT and ESIPT) for 8-hydroxy-4H-naphthalen-1-one (HNA), 5-hydroxynaphthoquinone (HNQ), 1-hydroxy-anthraquione (HAQ), 7-hydroxy-1-indenone (7HIN), 5,8-dihydroxynaphthoquinone (DHNQ) and 4,9-dihydroxyperylene-3,10-quinone (DHP) are studied at B3LYP/6-31G(d,p) and TD B3LYP/6-31G(d,p) level. The calculated results show that the PES of GSIPT for HNA, HNQ and HAQ exhibit a single minimum in the enol zone, while for 7-HIN, DHNQ and DHP exhibit a double minimum and a high barrier between the two minima. The barrierless ESIPT for HNA is predicted, however, the PES of ESIPT for HNQ, HAQ, 7HIN, DHNQ and DHP exhibit a high barrier in the S₁ tautomerism.     

Theoretical Study on the Structure and Vibrational Spectra for4-methyl-3-penten-2-one

The4-methyl-3-penten-2-oneisanimportanta,0-unsaturedketonemoleculeandanimportantligandoftheorganometalliccompounds.Untilnow,somepropertiesofphoto-chemistryandexcitedstateshavebeenexperimentallystudied'-#andconformationofthemoleculewerealsostudiedinexperiments"'.Standardinfraredgratingspectrumwasalsoobtainedin1970'.However,therearenodensityfunctionaltheory(DFT)calculationsofthismoleculeintheliterature.Recently,densityfunctionaltheoryhasbeenacceptedbytheahinitioquantumchemistrycommunityasacost…     

含不同质子供体的2-苯基苯并三唑衍生物激发态质子转移的理论研究

用密度泛函理论(DFT)和二级微扰理论(MP2)研究了带不同质子供体的2-苯基苯并三唑衍生物: 2-(2-羟苯基)苯并三唑(H-TIN), 2-(2-氨苯基)苯并三唑(APyBT)和2-(2-巯苯基)苯并三唑(MPyBT)的激发态分子内质子转移(ESIPT)性质以及它们作为紫外光吸收剂的光物理机制. 结果表明, 在基态时三个化合物的最稳定异构体是均存在分子内氢键的正常构型N, 而互变异构体T和其扭曲构型T_twisted都是不稳定的. 激发态势能曲线表明H-TIN和APyBT的ESIPT分别需要克服约7.06和20.7 kJ/mol的能垒, 而MPyBT的ESIPT无需能垒|同时结合分子轨道, 电荷差分密度三维立体图的分析结果表明三个化合物都能发生ESIPT, 并且伴随有扭曲分子内电荷转移, 这些原因均表明它们都具有好的紫外光稳定作用.     

5-氟胞嘧啶气相及水助质子转移异构化的理论研究

采用密度泛函B3LYP/6-311G**方法,对6种5-氟胞嘧啶异构体孤立分子的稳定性及质子转移引起的酮式-烯醇式、氨基式-亚胺式互变异构反应机理进行了计算研究,获得了零点能、吉布斯自由能及质子转移过程的反应焓、活化能、活化吉布斯自由能和速率常数等参数.计算结果表明,气相中烯醇-氨基式FC4是最稳定的异构体.分子内质子转移设计了FC1→FC2和FC1→FC6两条通道,分别标记为P(1)和P(2),各通道速控步骤的活化能和速率常数分别为155.9 kJ·mol-1,4.70×10-15 s-1和173.1 kJ·mol-1,1.41×10-18 s-1.水助催化时,相应通道P(3) 和P(4) 速控步骤的活化能和速率常数分别为51.0 kJ·mol-1,1.41×103 s-1和88.2 kJ·mol-1,4.53×10-3 s-1.可见,水分子的加入极大地降低了质子转移的活化能垒.另外发现,水分子参与形成协同的双质子转移机理比水助单质子转移机理更利于降低活化能垒.     

Theoretical study of the intramolecular CH/π interaction effect on rotation energy barriers in 1-pentene, 2,2′-diisopropyl biphenyl and some amino and nitro derivatives

The 2,2′-diisopropyl biphenyl conformers, and their amino and nitro para-disubstituted derivatives present two typical characteristics: In the ground state, CH/π interactions may induce local structures by positioning H atoms above some C atoms of the unsaturated cycles, and next the main skeleton is of biphenyl type. From ab initio theoretical calculations, we analyse first these characteristics separately by considering smaller molecules, i.e. 1-pentene for CH/π interactions and the biphenyl molecule itself. Sophisticated methods can be used for 1-pentene. We point out that the CH/π interaction present in the syn-conformer is not sufficiently stabilizing to compensate steric repulsions and the anti-conformer is found as the ground state. In the case of the biphenyl molecule, like many authors did before, and experiment compared, we were not able to improve significantly the calculated rotation energy barrier between the ground state and the conformation with a coplanar arrangement of the π-cycles. MP2 and B3LYP calculations, with basis sets of double-ζ plus polarization quality, on 2,2′-diisopropyl biphenyl conformers and their amino and nitro para-disubstituted derivatives, emphasize the difficulties found with 1-pentene and biphenyl, but damped down. The electron-donating and electron-withdrawing group effects of the amino and nitro substituents are analysed in term of σ and π contributions. This is mainly a π effect, which imposes its behaviour to the total electronic population only when considering also the atoms bonded to the amino and nitro groups. An increase of the electronic population on the atoms of the CH/π bond, mainly located on the C atom of the π system, is observed and is rather a σ effect. Moreover we show that CH/π interactions in the ground state only arise between H atoms of CH in iPr groups with C_ipso, not with C_ortho, and fit the experimental substituent effect on the basicity of the aryl group.     

Theoretical study of the reactivity of (CH3)2AlCH2I promoted cyclopropanation reactions

Density functional theory calculations are reported for the cyclopropanation reactions of (CH₃)₂AlCH₂I with ethylene for two reaction channels: methylene transfer and carbometalation. These computational results suggest that the methylene transfer process is favored and the competition from the carbometalation pathway is negligible.     

Coupling characteristics and proton transfer mechanisms of guanine–Na monohydrate

The coupling characteristics and the proton transfer mechanisms of guanine–Na⁺ monohydrate are determined in this investigation after the implementation of the geometry optimization and the harmonic vibrational frequency calculations. There are two elementary coupling modes: the interaction of monohydrated sodium ion with two heteroatoms which form a ringed coupling, and hydrogen-bond involved coupling mode. Two potential reaction pathways, coupling mode and hydration have been taken into account, and the accurate values of binding energy are corrected for basis set superposition error (BSSE) and zero-point vibrational energy (ZPVE). Relative energies of the hydrated guanine–sodium ion complexes indicate that the ringed-coupling complexes are predominant geometries with much lower energies. Monohydrated sodium ion coupling with O6 and N7 generates the most stable geometry with a five-member cycle. Sodium ion plays an important role in the tautomerization for guanine–sodium ion complexes. This investigation indicates that the stable cation-π complexes cannot be optimized for guanine–sodium ion monohydrate. Amino-involved coupling often gives rise to a twisted four-membered cycle with unrealistic distribution of positive charge and higher energies. The rotation of amino group is likely to lead to the redistribution of the base pair hydration bonding. Effective distribution of the positive charge is an important factor in the stabilization of biological systems and binding energies for the monohydrated guanine–sodium ion complexes. The enolic coupling complex has the higher energy than the keto type due to the hindrance for the positive charge.     

2-(benzimidazol-2-yl)-3-hydroxychromone derivatives: spectroscopic properties and a possibile alternative intramolecular proton phototransfer

Three 2-(benzimidazol-2-yl)-3-hydroxychromone derivatives were synthesized. Their spectroscopic and fluorescent properties, due to excited state intramolecular proton transfer (ESIPT) from OH to carbonyl, were studied. Theoretical possibility of an alternative intramolecular H-bonding and experimental evidence for such behavior are discussed.      

Theoretical Studies on the Structure and Spectrum of Imidazole-Chloranil Charge Transfer Complex

UV-Vis absorption spectra of the molecular complex formed by imidazole (Im) and chloranil (CA) were measured in chloroform. The stoichiometry of the imidazole-chloranil (Im-CA) complex was determined as 1:1 by applying Benesi-Hildebrand's equation and Job's continuous variation method. Density function theory (DFT) and MP2 calculations were performed to study the structures and the binding energies of the Im-CA complex. The calculations located four conformations (denoted as S1-S4) for the Im-CA complex, two edge(Im)-to-face(CA) linked and two edge(Im)-to-edge(CA) linked. It was found that the edge-to-face conformers are more stable than the edge-to-edge ones. The bonding characteristics of these conformers were investigated with natural population analysis (NPA), topological analysis of electron density, and natural bond orbital (NBO) analysis. It was revealed that the edge-to-face conformers are charge-transfer (CT) complexes whereas the edge-to-edge conformers are the hydrogen bond complexes. For the most stable conformation of the Im-CA complex (S1), the charge transfer interaction of the imidazole n(N15) lone pair orbital with the chloranil ∏*(C1=O7) orbital plays a crucial role in the Im-CA binding, and the binding is further strengthened by the O7… H20 hydrogen bond. The electronic excitation energies of the complex (S1) were calculated with time-dependent DFT (TDDFT), and the observed UV-Visible spectrum of the complex was analyzed based on the computed results.     

Theoretical study on the structure and UV–visible spectrum of pyridine–chloranil charge-transfer complex

The ground-state structure of the charge-transfer complex formed by pyridine (Py) as electron donor and chloranil (CA) as acceptor has been studied by full geometry optimization at the MP2 and DFT levels of theory. Binding energies were calculated and counterpoise corrections were used to correct the BSSE. Both MP2 and DFT indicate that the pyridine binds with chloranil to form an inclined T-shape structure, with the pyridine plane perpendicular to the chloranil. The CP and ZPE corrected binding energies were calculated to be 14.21 kJ/mol by PBEPBE/6-31G(d) and 23.21 kJ/mol by MP2/6-31G(d). The charge distribution of the ground state Py–CA complex was evaluated with the natural population analysis, showing a net charge transfer from Py to CA. Analysis of the frontier molecular orbitals reveals a σ–π interaction between CA and Py, and the binding is reinforced by the attraction of the O₇ atom of CA with the H₂₃ atom of Py. TD-DFT calculations have been performed to analyze the UV–visible spectrum of Py–CA complex, revealing both the charge transfer transitions and the weak symmetry-relieved chloranil π–π^* transition in the UV–visible region.