Conformers and intramolecular hydrogen bonding of the oxalic acid monomer and its anions

Density functional theory, B3LYP/6‐31G** and B3LYP/6‐311+G(2d,p), and ab initio MP2/6‐31G** calculations have been carried out to investigate the conformers, transition states, and energy barriers of the conformational processes of oxalic acid and its anions. QCISD/6‐31G** geometrical optimization is also performed in the stable forms. Its calculated energy differences between the two most stable conformers are very near to the related observed value at 7.0 kJ/mol. It is found that the structures and relative energies of oxalic acid conformers predicted by these methods show similar results, and that the conformer L1 (C_2h) with the double‐interfunctional‐groups hydrogen bonds is the most stable conformer. The magnitude of hydrogen bond energies depends on the energy differences of various optimized structures. The hydrogen bond energies will be about 32 kJ/mol for interfunctional groups, 17 kJ/mol for weak interfunctional groups, 24 kJ/mol for intra‐COOH in (COOH)₂, and 60 kJ/mol for interfunctional groups in (COOH)COO⁻¹ ion if calculated using the B3LYP/6‐311+G(2d,p) method. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 541–551, 2000     

二氯亚甲基锗烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet dichloromethylene germylene and ethylene has been investigated with B3LYP/6-31G＊ method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies for the involved conformations were calculated by CCSD（T）//B3LYP/6-31G＊ method. On the basis of the surface energy profile obtained with CCSD（T）// B3LYP/6-31G＊ method for the cycloaddition reaction between singlet dichloromethylene germylene and ethylene, it can be predicted that the dominant reaction pathway is that an intermediate INT1 is firstly formed between the two reactants through a barrier-free exothermic reaction of 61.7 kJ/mol, and the intermediate INT1 then isomerizes to an active four-membered ring product P2.1 via a transition state TS2, an intermediate INT2 and a transition state TS2.1, in which energy barriers are 57.7 and 42.2 kJ/mol, respectively.     

Broken symmetry density functional study of a mixed‐valence unsymmetrical dinuclear iron complex

Density functional theory (DFT) calculations with different exchange‐correlation functionals were performed for a mixed valence Fe(II)/Fe(III) binuclear complex with μ‐methoxo and two μ‐carboxylate bridging ligands, (1) with geometry optimizations being performed for all possible spin multiplicities (M_S = 2, 4, 6, 8, and 10). Within the exchange‐correlation functionals studied, only the hybrid GGA functionals B3P and B3LYP and also the pure GGA functional RPBE, predicts the geometry with high spin (S = 9/2) to be more stable than the geometry with low spin state (S = 1/2) by 20 kcal/mol, in agreement with the experimental findings. These functionals also predict the same stability order for the different spin states, being M_S = 10>8>6>2>4. The meta‐GGA functionals TPSS and TPSSh and also the pure GGA functionals BLYP and BP86 predict different stability orders. The computed average EPR g‐tensor, g_av, of 2.03, at the B3LYP level, is in good agreement with the experimental findings. Heisenberg exchange coupling constants, J, were calculated within the broken‐symmetry formalism, at the B3LYP level, showing that the two iron centers are antiferromagnetic coupling, with a very weak coupling constant of about ?7 cm^?1, in good agreement with the experimental value. Additionally, the effect of using different multiplicities of the reference geometries on the computed J value is discussed. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

尿酸质子转移反应的理论研究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难以进行。     

Bioinorganic chemistry modeled with the TPSSh density functional

In this work, the TPSSh density functional has been benchmarked against a test set of experimental structures and bond energies for 80 transition-metal-containing diatomics. It is found that the TPSSh functional gives structures of the same quality as other commonly used hybrid and nonhybrid functionals such as B3LYP and BP86. TPSSh gives a slope of 0.99 upon linear fitting to experimental bond energies, whereas B3LYP and BP86, representing 20% and 0% exact exchange, respectively, give linear fits with slopes of 0.91 and 1.07. Thus, TPSSh eliminates the large systematic component of the error in other functionals, reducing rms errors from 46-57 to 34 kJ/mol. The nonhybrid version of the functional, TPSS, gives a slope of 1.08, similar to BP86, implying that using 10% exact exchange is the main reason for the success of TPSSh. Typical bioinorganic reactions were then investigated, including spin inversion and electron affinity in iron-sulfur clusters, and breaking or formation of bonds in iron proteins and cobalamins. The results show that differences in reaction energies due to exact exchange can be much larger than the usually cited approximately 20 kJ/mol, sometimes exceeding 100 kJ/mol. The TPSSh functional provides energies approximately halfway between nonhybrids BP86 and TPSS, and 20% exact exchange hybrid B3LYP: Thus, a linear correlation between the amount of exact exchange and the numeric value of the reaction energy is observed in all these cases. For these reasons, TPSSh stands out as a most promising density functional for use and further development within the field of bioinorganic chemistry.     

Effect of trimethylsilyl substitution on the structure and properties of phthalocyanine: A density functional theory study

The geometries of four isomers of the trimethylsilyl substituted phthalocyanine (Pc)— I , II , III , and IV —have been optimized at the B3LYP/3‐21G level of density functional theory. Normal‐mode vibrational analyses have been performed and their standard thermodynamic functions, molar fractions, and electronic absorption spectra calculated. Single‐point energies have been calculated at the B3LYP/6‐311G* level for all isomers to evaluate the heats of formation from an isodesmic reaction. It is found that substitution has little influence on the geometry and electronic structures of the Pc framework. The corresponding geometric parameters in various isomers are close. According to the B3LYP/6‐311G*//B3LYP/3‐21G results, substitution at the peripheral position of the isoindole with an inner hydrogen is most favorable. The energies increase in the order of IV < II < III < I , and the energy difference between IV and I is 5.75 kJ/mol. The molar fractions of IV , II , III , and I are 0.80, 0.17, 0.02, and 0.02 and the heats of formation are 2009.96, 2010.10, 2015.85, and 2016.52 kJ/mol, respectively. This indicates that nonperipheral substituted Pcs have higher energy and little production because they are not stable under the considered conditions. The electronic spectra of the substituted Pcs calculated using the ZINDO method have two strong Q absorption bands around 700 nm and one B band around 300 nm that are slightly shifted compared with those in Pc. The ratios of the oscillator strength of the B band to the Q bands are much lowered by trimethylsilyl substitution. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

亚甲基硅烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet methylidenesilene and ethylene has been investigated with MP2/6-31G^＊ and B3LYP/6-31G^＊ methods, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies of the involved conformers were calculated by CCSD（T）//MP2/6-31G＊ and CCSD（T）//B3LYP/6-31 G＊ methods, respectively. The results show that the dominant reaction pathway of the cycloaddition reaction is that a complex intermediate is firstly formed between the two reactants through a barrier-free exothermic reaction of 13.3 kJ/mol, and the complex is then isomefized to a four-membered ring product P2,1 via a transition state TS2.1 with a barrier of 32.0 kJ/mol.     

亚烷基锡烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet alkylidenestannylene and ethylene has been investigated with MP2/3-21 G^* and B3LYP/3-21 G* methods, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies for the involved conformations were calculated by CCSD(T)//MP2/3-2 IG^* and CCSD(T)//B3LYP/3-21G^* methods, respectively. The results show that the dominant reaction pathway of the cycloaddition is that an intermediate (INT) is firstly formed between the two reactants through a barrier-free exothermic reaction of 39.7 kJ/mol, and the intermediate then isomerizes to a four-membered ring product (P2.1) via a transition state TS2.1 with a barrier of 66.8 kJ/mol.     

DFT Study of Hydrogen-Bonded 1,3,5-Triazine-Water Complexes

The 1,3,5-triazine-water hydrogen bonding interactions have been investigated using the density functional theory B3LYP method and 6-31 ＋＋G^＊＊ basis, obtaining one, two and seven energy minima of the ground states for the 1,3,5-triazine-water, 1,3,5-triazine-（water）2 and 1,3,5-triazine-（water）3 complexes respectively. The fully optimized geometries and binding energies were reported for the various stationary points. The global minima of 1,3,5-triazine-（water）2 and 1,3,5-triazine-（water）3 complexes have a hydrogen bond N…H-O and a chain of water molecules, terminated by a hydrogen bond O…H-C. The binding energies are 13.38, 39.52 and 67.79 kJ/mol for the most stable 1,3,5-triazine-water, 1,3,5-triazine-（water）2 and 1,3,5-triazine-（water）3 complexes respectively, after the basis set superposition error and zero point energy corrections. The H-O symmetric stretching modes of water in the complexes are red-shifted relative to those of the monomer water. In addition, the NBO analysis indicates that inter-molecule charge transfer is 0.02145 e, 0.02501 e and 0.02777 e for the most stable 1 ： 1, 1 ： 2 and 1 ： 3 complexes between 1,3,5-triazine and water, respectively.     

铁原子与氮分子间的相互作用——单端位构型

用杂化密度泛函B3LYP方法在6-311+G(d)基组水平上研究了Fe 原子与N2分子相互作用的单端位构型的直线形和弯曲形两种结构的平衡几何结构、电子结构、轨道布局及红外光谱等性质. 计算结果表明, 由于强的σ-σ电子对互斥作用, 基组态4s²3d⁶的Fe原子不能与N2分子发生化学作用; 当Fe 原子呈现可与N2之间发生σ-π授予反馈作用的激发组态时, Fe 与N2分子之间可形成稳定的结构; 在得到的多个电子态中, 能量最低的是直线形的1³∑^-, 比Fe(a⁵D)和N2(¹∑⁺g )能量高21.6 kJ·mol-1, 同时存在几个能量相近的电子态, 如1³∏、1³Φ; 弯曲形都是不稳定态, 可能是连接直线形和单侧双配位构型的过渡态; 单端位构型产物相对于基态的反应物均是热力学不稳定的; 单端位构型中Fe对N2的活化作用很小, N—N 键长增加不超过7 pm.     

Density Functional Theory Study on Hydrogen Bonding Interaction of Catechin-(H2O)n

Density functional theory B3LYP method with 6‐31G* basis set has been used to optimize the geometries of the catechin, water and catechin‐(H₂O)_n complexes. The vibrational frequencies have been studied at the same level to analyze these complexes. Six and eleven stable structures for the catechin‐H₂O and catechin‐(H₂O)₂ have been found, respectively. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) have been utilized to investigate the hydrogen bonds involved in all the systems. The interaction energies of all the complexes corrected by basis set superposition error, are from ?13.27 to ?83.56 kJ/mol. All calculations also indicate that there are strong hydrogen‐bonding interactions in catechin‐water complexes. The strong hydrogen‐bonding contributes to the interaction energies dominantly. The O–H stretching motions in all the complexes are red‐shifted relative to that of the monomer.     

B20N20团簇的结构与稳定性

B3LYP/6-31G* density functional theory calculations have been carried out on the structure and stability of ten B20N20 clusters. It was found that two new proposed isomers with two octagons, twelve hexagons, eight squares in Cab and C2 symmetry were more stable than the isomer with sixteen hexagons and six squares in C2 symmetry which was previously deemed to the most stable by 79.5 and 13.8 kJ/mol respectively. The isomer with two decagons in S10 symmetry is much higher in energy than the most stable structure in C4h symmetry by 637.2 kJ/mol.     

Comparative studies on the structures,infrared spectrum,and thermodynamic properties of phthalocyanine using ab initio Hartree–Fock and density functional theory methods

The molecular structure, vibrational spectrum, standard thermodynamic functions, and enthalpy of formation of free base phthalocyanine (Pc) have been studied using the density functional theory B3LYP procedure, as well as the ab initio Hartree–Fock method. Various basis sets 3‐21G, 6‐31G*, and LANL2DZ have been employed. The results obtained at various levels are discussed and compared with each other and with the available experimental data. It is shown that calculations performed at the Hartree–Fock level cannot produce a reliable geometry and related properties such as the dipole moment of Pc and similar porphyrin‐based systems. Electron correlation must be included in the calculations. The basis set has comparatively less effect on the calculated results. The results derived at the B3LYP level using the smaller 3‐21G and LANL2DZ basis sets are very close to those produced using the medium 6‐31G* basis set. The geometry of Pc obtained at the B3LYP level has D_2h symmetry and the diameter of the central macrocycle is about 4 Å. The enthalpy of formation of Pc in the gas phase has been predicted to be 1518.50 kJ/mol at the B3LYP/6‐311G(2d,2p)//B3LYP/6‐31G* level via an isodesmic reaction. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Theoretical Study on the Intermolecular Interactions between 1,1‐Diamino‐2,2‐Dinitroethylene and H2O

The density functional method was applied to the study of 1,1‐diamino‐2,2‐dinitroethylene (Fox‐7)/H₂O dimer. All the possible dimers ( 1, 2 and 3 ), as well as the monomers, were fully optimized with the DFT method at the B3LYP/6‐311++G^** level. The basis set superposition errors (BSSE) are 4.62, 4.07 and 3.45 kJ/mol, and the zero point energy (ZPE) corrections for the interaction energies are 7.94, 5.66 and 6.40 kJ/mol for 1, 2 and 3 , respectively. Dimer 1 is the most stable, judged by binding energy. After BSSE and ZPE corrections, the greatest corrected intermolecular interaction energy of dimer 1 was predicted to be ?29.36 kJ/mol. The charge redistribution mainly occurs on the adjacent N–H··· O atoms and N–O··· H atoms between submolecules. The oxygen in the nitro group acts as a moderate hydrogen acceptor as compared to water oxygen. Based on the statistical thermodynamic method, the standard thermodynamic functions, heat capacities (C⁰_P), entropies (S⁰_T) and thermal corrections to enthalpy (H⁰_T), and the changes of thermodynamic properties on going from monomer to dimer over the temperature range 200.00‐700.00 K were predicted. It is energetically or thermodynamically favorable for Fox‐7 to bind with H₂O and to form dimer 1 at room temperature.     

Kinetic Study and NBO Analysis of the Dehydrogenation Mechanism of Five‐membered Ring Heterocyclic 2,5‐Dihydro‐[furan,thiophene, selenophene

The theoretical study of the dehydrogenation of 2,5‐dihydro‐[furan ( 1 ), thiophene ( 2 ), and selenophene ( 3 )] was carried out using ab initio molecular orbital (MO) and density functional theory (DFT) methods at the B3LYP/6‐311G**//B3LYP/6‐311G** and MP2/6‐311G**//B3LYP/6‐311G** levels of theory. Among the used methods in this study, the obtained results show that B3LYP/6‐311G** method is in good agreement with the available experimental values. Based on the optimized ground state geometries using B3LYP/6‐311G** method, the natural bond orbital (NBO) analysis of donor‐acceptor (bond‐antibond) interactions revealed that the stabilization energies associated with the electronic delocalization from non‐bonding lone‐pair orbitals [LP(e)_X3] to δ*_{C(1)  H(2)} antibonding orbital, decrease from compounds 1 to 3 . The LP(e)_X3→δ*_{C(1)  H(2)} resonance energies for compounds 1 – 3 are 23.37, 16.05 and 12.46 kJ/mol, respectively. Also, the LP(e)_X3→δ*_{C(1)  H(2)} delocalizations could fairly explain the decrease of occupancies of LP(e)_X3 non‐bonding orbitals in ring of compounds 1 – 3 ( 3 > 2 > 1 ). The electronic delocalization from LP(e)_X3 non‐bonding orbitals to δ*_{C(1)  H(2)} antibonding orbital increases the ground state structure stability, Therefore, the decrease of LP(e)_X3→δ*_{C(1)  H(2)} delocalizations could fairly explain the kinetic of the dehydrogenation reactions of compounds 1 – 3 (k ₁ >k ₂ >k ₃ ). Also, the donor‐acceptor interactions, as obtained from NBO analysis, revealed that the (_C(4)C(7)→δ*_{C(1)  H(2)} resonance energies decrease from compounds 1 to 3 . Further, the results showed that the energy gaps between (_C(4)C(7) bonding and δ*_{C(1)  H(2)} antibonding orbitals decrease from compounds 1 to 3 . The results suggest also that in compounds 1 – 3 , the hydrogen eliminations are controlled by LP(e)→δ* resonance energies. Analysis of bond order, natural bond orbital charges, bond indexes, synchronicity parameters, and IRC calculations indicate that these reactions are occurring through a concerted and synchronous six‐membered cyclic transition state type of mechanism.     

Theoretical study on the isomerization and tautomerism in barbituric acid

Isomerization and tautomerism of 16 isomers of barbituric acid (BA) were studied at the MP2 and B3LYP levels of theory. Activation energies (E _a), imaginary frequencies (υ), and Gibbs free energies (ΔG ^#) of the amine-imine and keto-enol tautomerisms and O–H internal rotations were calculated. The activation energies of amine-imine tautomerisms were in the range of 110–200 kJ/mol and for keto-enol tautomerisms were larger than 200 kJ/mol. The calculated activation energies of internal O–H rotations were smaller than 60 kJ/mol. Effect of micro-hydration on the transition state structures and activation energies of the tautomerisms were also investigated. Water molecule catalyzed the tautomerisms and decreased the activation energies of both the amine-imine and keto-enol tautomerisms about 100–120 kJ/mol.     

A Theoretical Study on the Equilibrium Structures and Relative Stabilities of H_2SO

1INTRODUCTIONTheconformationandvariousphysicalprope-rtiesofdichalcogencompoundsR-X-Y-R?(X,Y=O,S,Se,andTe)havebeenatopicofseveralexperi-mentalandtheoreticalinvestigations[1].CompoundscontainingS=SandS=Obondshavelongbeenpro-posedasintermediatesinorganicsynthesisand,onoccasion,asstableentities[2～18].Despitethecom-monlyusedrepresentationofS–Obondinsulfoxi-desandothersulfinylderivativesasS=O,sulfoxidesareinmanywaysbestdescribedasylideswithhighlypolarizedS–Oσ-bondbecauseofelectro-sta…     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.