3-硝基-1,2,4-三唑-5-酮与NH3及H2O分子间相互作用的理论研究

在DFT—B3LYP／6—311 G^**水平上,求得3-硝基-1,2,4-三唑-5-酮(NTO)／NH3和NTO／H2O两种超分子体系势能面上5种全优化构型．经基组叠加误差(BSSE)和零点能(ZPE)校正,求得NTO与NH3和H2O的分子间最大相互作用能依次为-37．58和-30．14kJ／mol,表明NTO与NH3的分子问相互作用强于与H2O的作用．超分子体系中电子均由NH3或H2O向NTO转移,相互作用能主要由强氢键所贡献,由自然键轨道分析揭示了相瓦作用的本质．对优化构型进行振动分析,并基于统计热力学求得200．0～800．0K温度范围从单体形成超分子的热力学性质变化．发现由NTO和NH3形成超分子Ⅱ和Ⅲ在常温下可自发进行;而NTO和H2O只在低温下才能自发形成Ⅳ,Ⅴ和Ⅵ超分子．     

NH_(2-3)~(0-1 )离解能等的高级ab initio计算与评价

FeiQi等人[1]最近对H2N－H、H2N －H的离解能以及NH3、NH2的电离能，用真空紫外光电离实验进行了重新测定，并同时得到了H2N－H 的离解能．他们还在QCISD（T）／6－311＋G（3df，ZP）／／MPZ（full）／6刁IG（d）水平（GZ理论的参考水平）上对这些数据及其它相关分子的某些性质     

类锗烯H2C=GeLiCl与RH(R=F,OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了不饱和类锗烯H2C=GeLiCl与RH(R=F, OH, NH2)的插入反应. 在B3LYP/6-311+G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C=GeLiCl与HF、H2O 或NH3发生插入反应的机理相同. QCISD/6-311++G(d,p)//B3LYP/6-311+G(d,p)计算的三个反应的势垒分别为173.53、194.48和209.05 kJ·mol-1, 反应热分别为60.18、72.93和75.34 kJ·mol-1. 相同条件下发生插入反应时, 反应活性顺序都是H—F>H—OH>H—NH2.     

三重态类硅烯HB＝SiLiF的结构及其与R—H(R=F,OH,NH_2)的插入反应

用密度泛函理论(DFT)和二次组态相互作用(QCISD)方法研究了三重态类硅烯HB=SiLiF的结构及其与RH(R=F,OH,NH2)的插入反应.计算结果表明,类硅烯HB=SiLiF有三种平衡构型,其中四元环构型能量最低,是其存在的主要构型.HB=SiLiF与HF,H2O和NH3发生插入反应的机理相同.QCISD/6-311++G(d,p)//B3LYP/6-311+G(d,p)计算的三个反应的势垒分别为124.85,140.67和148.16kJ·mol-1,反应热分别为-2.22,20.08和23.22kJ·mol-1.相同条件下发生插入反应时,反应活性都是H—FH—OHH—NH2.     

类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应. 在B3LYP/6- 311＋G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C＝SiLiBr与HF, H2O或NH3发生插入反应的机理相同. QCISD/6-311＋＋G(d,p)//B3LYP/6-311＋G(d,p)计算的三个反应的势垒分别为148.62, 164.42和165.07 kJ•mol－1, 反应热分别为－69.63, －43.02和－28.27 kJ•mol－1. 相同条件下发生插入反应时, 反应活性都是H—F＞H—OH＞H—NH2.     

1,2,4-三氮杂苯-(H2O)n复合物氢键相互作用的密度泛函理论研究

用密度泛函理论方法在B3LYP／6—31 G**水平上对1，2，4-三氮杂苯-(H2O)n(n=1，2，3)氢键复合物的基态进行了结构优化和能量计算，结果表明复合物之间存在较强的氢键作用，所有稳定复合物结构中形成一个N…H--O氢键并终止于弱O…H—C氢键的氢键水链的构型最稳定．同时，用含时密度泛函理论方法(TD—DPT)在TD—B3LYP／6—31 G**水平上计算了1,2，4-三氮杂苯单体及其氢键复合物的单重态第-1(n，π*)垂直激发能．     

2002年全国高中学生化学竞赛理论试题答案

第1题1－1(NH4)3H4PMo12O42·H2O＋26NaOH＝12Na2MoO4＋Na2HPO4＋3NH3＋17H2ONH3＋HCI＝NH4 Cl计量关系 P2O5－(NH4)3H4PMo12O42·H2O－NaOH1－2－46P2O5的摩尔质量:142g·mol-1P2O5%:(1．026mol·dm－3×0．04000dm-3－0．1022mol·dm－30．01520dm-3×142g·mol-1(46×0．385g)＝0．317＝31．7%(计量关系正确得3分,计算合理结果错扣1分,有效数字错扣0．5分)     

C6H5—H…X分子间氢键的理论计算

利用密度泛函理论B3LYP方法, 在6-311+G(3df,2p)水平上对C6H5—H…X型分子间氢键进行了几何构型优化、氢键相互作用能、电子密度分布等计算. 其中C6H6为质子供体, HCOH、H2O、NH3、CH2NH和HCN为质子受体. 从电荷布居分析、自然键轨道等角度详细地讨论了C6H5—H…X 体系中, 共轭π键、O和N的不同键型结构对氢键形成的影响以及孤电子对与C—H 反键轨道之间的相互作用(n→σ*)等.     

催化动力学分光光度法测定微量铅的研究

本文报道在NH3-NH4C1介质中，利用铅催化H2O2氧化胭脂红腿色这一反应，建立催化动力学光度法测定微量铅的新方法。方法的检出限为0．091μg／25mL，线性范围0．10-0.90μg／25mL，用于分析汽油中的铅，结果满意。     

甲醛与甲酸相互作用的从头算研究

在MP2／6—31G(d)和MP2(FC)／6—311 G(d，p)水平上，对H2CO和HCOOH以及设计的4种构型H2CO…HCOOH复合物等进行几何全优化计算，经振动频率分析，确认它们为势能面上的稳定驻点．然后在MP2／6—311 G(2df，2p)水平上进行单点能计算和基组重叠误差(BSSE)校正以获得相互作用能，并利用自然键轨道理论探讨H2CO和H(X)OH相互作用的本质。     

NH0-1+2-3离解能等的高级ab initio计算与评价

A recent experimental determination[1] of the dissociation energies (D0) for H2N-H, H2N+-H and H2N-H+, the ionization energies for NH3 and NH2 resulted in large deviations when compared with those of the earlier values and the QCISD(T)/6-311+G(3df,2p) ab initio calculations. We have performed some higher level ab initio calculations on these data by utilizing the Gaussian 92/DFT and Gaussian 94 pakages of programs and have assessed the available experimental values. Our calculations were carried out at the QCISD (TQ)/aug-cc-pVDZ, G2(QCI), QCISD(T)/6-311 ++G(3df,3pd) and QCISD(T)/aug-cc-pVTZ levels of theory. Geometries were optimized at both of the MP2(full)/6-31G(d) and the MP2(full)/6-31(d,p) levels, and were compared with those of the experiments if available. The MP2(full)/6-31G(d,p) tight-optimized geometries for the neutrals are closer to those of the experiments than those of the MP2 (full)/6-31G(d), and are in excellent agreement with the experimental results as shown in Table 1. In this case, we assumed that the optimized geometries for the cations would be better if p polarization functions are added to the hydrogen atoms. We firstly noted that the symmetry of the NH3+ cation was D3h, other than Cs. as reported in ref.[1]. All of the zero-point energies and the final geometries are calculated at the MP2(full)/6-31G(d,p) level of theory. We have also repeated the QCISD(T )/6-311 + G(3df,2p) calculations of ref. [1], because we could not identify their level of goemetry optimization. It is found that the total energy, -55.244 19 Hartrees, for NH2+(1A1 ) in ref.[1] might be in error. Our result is -55.336 29 Hartrees at the same level of theory. At our highest level [QCISD(T)/aug-cc-pVTZ] of calculations as shown in Table 3, the D0 (temperature at zero Kelvin) values of H2N-H, H2N+-H(3B1for NH2+ ) and H2N- H+ are 4.51, 5.49 and 8.00 eV, respectively. These data reported in re f.[1] were 4.97, 5.59 and 8.41 eV, respectively. Our result on D0(H2N-H) supports the work of ref.[2,3,5,6]. The ionization energies (IE) for NH3 and NH2 (3B1 for NH2+) at our highest level are 10.11 and 11.09 eV while in ref.[1] were 10.16 and 10.78 eV, respectively. For the latter, our result supports the experiment of ref.[3]. Our predicted D0 for HN2+-H and IE for NH2 (1A1 for each NH2+) are 6.80 and 12.39 eV, respectively. These values differ greatly from the predicted values (9.29 and 14.88 eV) of ref.[1] where the total energy of NH2+(1A1) might be in error. The D0 value for HN-H has not been found in ref.[1]. Our result supports the work of ref.[3]. We have also derived all of these values at the temperature of 298K and under the pressure of 101kPa at several levels of thoery as shown in Table 3. On examining the experiment of ref.[1] in detail, it is easy to find that all of the larger deviations might be from a too high value of the appearance potential of proton AP(H+). Indeed, ref.[1] has mentioned that the determintion of AP(H+), due to kinetic shift, would lead to a hihger value for the dissociation energy as has been pointed out by Berkowitz and Ruscic. In this work, we concluded that, besides some mistakes in the theoretical calculations of ref.[1], the dissociation energies for H2N-H and H2N-H+,the IE for NH2 (3B1 for NH2+) might also be unreliable and need to be re-examined.
     

Modern valence-bond description of chemical reaction mechanisms: the 1,3-dipolar addition of diazomethane to ethene

The electronic mechanism for the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane (CH2N2) to ethene (C2H4) is described through spin-coupled (SC) calculations at a sequence of geometries along the intrinsic reaction coordinate obtained at the MP2/6-31G(d) level of theory. It is shown that the bonding rearrangements occurring during the course of this reaction follow a heterolytic pattern, characterized by the movement of three well-identifiable orbital pairs, which are initially responsible for the pi bond in ethene and the C-N pi bond and one of the N-N pi bonds in diazomethane and are retained throughout the entire reaction path from reactants to product. Taken together with our previous SC study of the electronic mechanism of the 1,3-dipolar cycloaddition of fulminic acid (HCNO) to ethyne (C2H2) (Theor. Chim. Acc. 1998, 100, 222), the results of the present work suggest strongly that most gas-phase concerted 1,3-dipolar cycloaddition reactions can be expected to follow a heterolytic mechanism of this type, which does not involve an aromatic transition state. The more conventional aspects of the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane to ethene, including optimized transition structure geometry, electronic activation energy, activation barrier corrected for zero-point energies, standard enthalpy, entropy and Gibbs free energy of activation, have been calculated at the HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), MP2/6-31G(d,p), QCISD/6-31G(d) and CCD/6-31G(d) levels of theory. We also report the CCD/6-311++G(2d, 2p)//CCD/6-31G(d), MP4(SDTQ)/6-311++G(2d,2p)//CCD/6-31G(d) and CCSD(T)/6-311++G(2d, 2p)//CCD/6-31G(d) electronic activation energies.     

Theoretical conformational analysis for neurotransmitters in the gas phase and in aqueous solution. Norepinephrine

The natural neurotransmitter (R)-norepinephrine takes the monocationic form in 93% abundance at the physiological tissue pH of 7.4. Ab initio and DFT/B3LYP calculations were performed for 12 protonated conformers of (R)-norepinephrine in the gas phase with geometry optimizations up to the MP2/6-311++G level, and with single-point calculations up to the QCISD(T) level at the HF/6-31G-optimized geometries. Four monohydrates were studied at the MP2/6-31G//HF/6-31G level. In the gas phase, the G1 conformer is the most stable with phenyl.NH(3)(+) gauche and HO(alc).NH(3)(+) gauche arrangements. A strained intramolecular hydrogen bond was found for conformers (G1 and T) with close NH(3)(+) and OH groups. Upon rotation of the NH(3)(+) group as a whole unit about the C(beta)-C(alpha) axis, a 3-fold potential was calculated with free energies for barriers of 3-12 kcal/mol at the HF/6-31G level. Only small deviations were found in MP2/6-311++G single-point calculations. A 2-fold potential was calculated for the phenyl rotation with free energies of 11-13 kcal/mol for the barriers at T = 310 K and p = 1 atm. A molecular mechanics docking study of (R)-norepinephrine in a model binding pocket of the beta-adrenergic receptor shows that the ligand takes a conformation close to the T(3) arrangement. The effect of aqueous solvation was considered by the free energy perturbation method implemented in Monte Carlo simulations. There are 4-5 strongly bound water molecules in hydrogen bonds to the conformers. Although hydration stabilizes mostly the G2 form with gauche phenyl.NH(3)(+) arrangement and a water-exposed NH(3)(+) group, the conformer population becomes T > G1 > G2, in agreement with the PMR spectroscopy measurements by Solmajer et al. (Z. Naturforsch. 1983, 38c, 758). Solvent effects reduce the free energies for barriers to 3-6 and 9-12 kcal/mol for rotations about the C(beta)-C(alpha) and the C(1)(ring)-C(beta) axes, respectively.     

CH3 CH2 +N(4S)反应机理的理论研究

The reaction for CH3CH2+N(4S) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single point calculations for all the stationary points were carried out at the QCISD(T)/ 6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2CH2+3NH and H2CN＋CH3, and the minor products are the CH3CHN+H in the reaction. The majority of the products CH2CH2+3NH are formed via a direct hydrogen abstraction channel. The products H2CN＋CH3 are produced via an addition/dissociation channel. The products CH3CHN+H are produced via an addition/dissociation channel.     

Molecular structure of the octatetranyl anion, C8H-: a computational study

The equilibrium molecular structure of the octatetranyl anion, C8H(-), which has been recently detected in two astronomical environments, is investigated with the aid of both ab initio post-Hartree-Fock and density functional theory (DFT) calculations. The model chemistry adopted in this study was selected after a series of benchmark calculations performed on molecular acetylene for which accurate gas-phase structural data are available. Geometry optimizations performed at the CCSD/6-311+G(2d,p), QCISD/6-311+G(2d,p), and MP4(SDQ)/6-311+G(2d,p) levels of theory yield for C8H(-) an interesting polyyne-type structure that defies the chemical formula displaying a simple alternation of triple and single carbon-carbon bonds, [:C[triple bond]C-C[triple bond]C-C[triple bond]C-C[triple bond]CH](1-). In the optimized geometry of C8H(-), as one proceeds from the naked carbon atom on one side of the chain to the CH unit on the opposite side of the chain, the short (formally triple) carbon-carbon bonds decrease in length from 1.255 to 1.213 A whereas the long (formally single) carbon-carbon bonds increase (albeit only slightly) in length from 1.362 to 1.378 A (CCSD results). In striking contrast, both MP2 and DFT (B3LYP and PBE0) calculations fail in reproducing the pattern of the carbon-carbon bond lengths obtained with the CCSD, QCISD, and MP4 methods. The structures of three shorter n-even chains, C(n)H(-) (n = 2, 4, and 6), along with those of four n-odd compounds (n = 3, 5, 7, and 9) are also investigated at the CCSD/6-311+G(2d,p) level of theory.     

Weakly bound complexes of N2O: an ab initio theoretical analysis toward the design of N2O receptors

Ab initio calculations at MP2/6-311++G(2d,2p) and MP2/6-311++G(3df,3pd) computational levels have been used to analyze the interactions between nitrous oxide and a series of small and large molecules that act simultaneously as hydrogen bond donors and electron donors. The basis set superposition error (BSSE) and zero point energy (ZPE) corrected binding energies of small N2O complexes (H2O, NH3, HOOH, HOO*, HONH2, HCO2H, H2CO, HCONH2, H2CNH, HC(NH)NH2, SH2, H2CS, HCSOH, HCSNH2) vary between -0.93 and -2.90 kcal/mol at MP2/6-311++G(3df,3pd) level, and for eight large complexes of N2O they vary between -2.98 and -3.37 kcal/mol at the MP2/6-311++G(2d,2p) level. The most strongly bound among small N2O complexes (HCSNH2-N2O) contains a NH..N bond, along with S-->N interactions, and the most unstable (H2S-N2O) contains just S-->N interactions. The electron density properties have been analyzed within the atoms in molecules (AIM) methodology. Results of the present study open a window into the nature of the interactions between N2O with other molecular moieties and open the possibility to design N2O abiotic receptors.     

Proton affinity correlations between hydrogen and dihydrogen bond acceptors

Several series of hydrogen- and dihydrogen-bonded complexes with HCN, C2H2, HF, H2O, CH3CONH2, and CH3COOH as donors and H2O, MeOH, EtOH, MeOMe, NH3, NH2Me, NHMe2, NMe3, NEtMe2, and BH3-NMe3 as acceptors were investigated using the MP2/6-311++G(d,p) level of theory. The total lowering of the X-H stretching frequencies in the hydrogen-bonded complexes were linearly correlated with the proton affinities of the accepting bases. From comparison of hydrogen- and dihydrogen-bonded complexes, a scaling factor to estimate the exact proton affinity of a dihydrogen bond acceptor was developed. Further, the scaling factor involving linear donors (1.204) is marginally higher than that involving nonlinear donor molecules (1.162). Finally, it was found that, given identical conditions, a hydrogen bond will be about 16-20% stronger than a corresponding dihydrogen bond.     

Ab initio study of some peroxides. HOOH, CH3OOH and, CH3OOCH3

The geometries of HOOH, CH₃OOH, and CH₃OOCH₃, were optimized with different basis sets (3-21G, 6-31G^*(*) and D95^**) at different levels of theory (HF, MP2, MP4, and CI). HF/3-21G optimizations result in planar trans conformations for all three peroxides. HF/6-31G^** calculations predict skew conformations for HOOH and CH₃OOH, but a planar trans struture for CH₃OOCH₃. For the larger basis set the calculated bond lengths, especially the O-O bonds, are too short. Optimizations for HOOH including electron correlation at the MP2, MP3, MP4, CI, and CCD level improve the agreement for bond lengths and the OOH angle, but result in dihedral angles Which are too large by 3– 8°. In the case of CH₃OOCH₃, similar calculations at the MP2 and CI level predict planar trans structures instead of the experimentally observed skew conformation. On the other hand, MP4 single point calculations at MP2 optimized parameters result in a correct skew structure. For all three peroxides a computationally “economic” method, i.e., single point calculations at MP2 or MP4 level with HF/3-21G optimized parameters, result in close agreement between calculated and experimental structures.     

Theoretical calculations of glycine and alanine gas-phase acidities

The gas-phase acidities of glycine and alanine were determined by using a variety of high level theoretical methods to establish which of these would give the best results with accessible computational efforts. MP2, MP4, QCISD, G2 ab initio procedures, hybrid Becke3-LYP (B3LYP) and gradient corrected Becke-Perdew (BP) and Perdew-Wang and Perdew (PWP) nonlocal density functionals were used for the calculations. A maximum deviation of approximately 13 and 18 kJ/mol from experimental data was observed for the computed delta Hacid and delta Gacid values, respectively. The best result was obtained at G2 level, but comparable reliability was reached when the considerably less time consuming B3LYP, BP, and PWP density functional approaches were employed.