恢复蛋白Recoverin正则模式分析的理论研究

采用旋转平移块方法对Ca2+/豆蔻酰基开关进行了正则模式分析(NMA). 研究结果表明, 恢复蛋白(Recoverin)的T态的N-末端与C-末端易于发生刚体逆向旋转, 一旦结合Ca2+, 很容易发生构象变化, 形成具有双向构象转变特征的I态. I态是一个中间结构, 既可以发生构象回转到T态, 又可以继续相对旋转到R态, 使豆蔻酰基完全暴露, 从而行使其信号传导生物功能. 从低频振动模式分析可以看出, 恢复蛋白具有构象转变这一本质属性.     

稀土化合物Ho2Co14Fe3的磁性分析

在分析稀土磁性常用的哈密顿中，引入了自旋－自旋偶极作用，用参数拟合的方法对Ｈｏ２Ｃｏ１４Ｆｅ３的磁化曲线进行了分析，计算结果与实验数据符合得很好，并且提示了Ｈｏ２Ｃｏ１４Ｆｅ３磁矩的突变现象。     

Cis–trans isomerization and spin multiplicity dependences on the static first hyperpolarizability for the two-alkali-metal-doped saddle[4]pyrrole compounds

Doping two alkali-metal atoms (Li and Na) into the saddle-shaped saddle[4]pyrrole forms four new two-alkali-metal-doped compounds with alkalide or electride characteristic. They are cis-LiNa(saddle[4]pyrrole) isomers 1(singlet) and 2(triplet), and trans-Li(saddle[4]pyrrole)Na isomers 3(singlet) and 4(triplet). The four structures with all-real frequencies are obtained at the density functional theory (DFT) B3LYP/6-311+G(d) level. All calculations of electric properties have been carried out at the second order Møller–Plesset perturbation theory (MP2) level. The order of the β ₀ values is 3.54 × 10³ for trans-4(triplet) <1.51 × 10⁴ for cis-1(singlet) <3.57 × 10⁴ for cis-2(triplet) <2.34 × 10⁵ a.u. for trans-3(singlet). The static first hyperpolarizability (β ₀) depends on the cis–trans isomerization and spin multiplicity. The result demonstrates that the cis–trans isomerization and spin multiplicity controls of the second-order NLO response are possible.     

Theoretical study on structures and stability of C4P isomers

The structures, energetics, spectroscopies, and stabilities of doublet C(4)P isomeric species are explored at the DFT/B3LYP, QCISD, and CCSD(T) (singlet-point) levels. A total of 12 minimum isomers and 27 interconversion transition states are located. At the CCSD(T)/6-311G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest-lying isomer is a floppy CCCCP 1 (0.0 kcal/mol) mainly featuring a cumulenic structure |C=C=C=C=P*|, which differs much from the analogous C4N radical (|*C-C[triple bond]C-C[triple bond]N|). The quasi-linearity and the low bending mode of 1 are in contrast to the previous prediction. The second energetically followed isomer PC-cCCC 3 (14.9 kcal/mol) possesses a CCC ring-bonded to CP. The two low-lying isomers are separated by a high-energy ring-closure/open transition state (26.5 kcal/mol) and thus are very promising candidates for future laboratory and astrophysical detection. Furthermore, four high-energy isomers, that is, two bent isomers CCPCC 2 (68.4 kcal/mol) and CCPCC 2' (68.5 kcal/mol) and two cagelike species 10 (56.0 kcal/mol) and 11 (67.9 kcal/mol), are also stabilized by considerable barriers. The present work is the first detailed potential energy survey of CnP clusters and can provide useful information for the investigation of larger CnP radicals and for understanding the isomerism of P-doped C vaporization processes.     

Theoretical elucidation of the rhodium-catalyzed [4 + 2] annulation reactions

The reaction mechanism of the Rh-catalyzed [4 + 2] annulation of 4-alkynals with isocyanates is unraveled using density functional calculations. The reaction mechanisms of the model system and the real substituted system have been investigated and the results are compared. From our theoretical results based on the model and real substituted system, it is shown that (a) the rate-determining step is the Rh-H addition to the alkyne, (b) the formation of the cyclopentenone G and glutarimide K represents a severe competition, and (c) the product selectivity should be controlled by the amount of the isocyanates. In addition, it is demonstrated that there exist steric effects in the real substituted system, but missed in model system. Our calculations also show that although the results obtained on the model system could explain the mechanism in principle, the real substituted system could reflect the mechanism more exactly and make the reaction proceed with regioselectivity.     

NC3O分子体系的异构化及其结构和性能的理论研究

采用 DFT, QCISD 及 CCSD(T)方法, 对 NC3O 分子体系的异构化进行了系统的研究, 得到了5个稳定的异构体, 其中包括4个链状异构体和1个带有支链的平面结构异构体. 这5个稳定异构体都具有较高的动力学和热力学稳定性.     

预测[C,O,S]体系的稳定异构体

采用 DFT, QCISD及CCSD(T)方法对单重态和三重态的COS分子体系势能面进行了理论计算, 在 QCISD/6-311G(d) 水平上得到4个过渡态连接的6个稳定构型. 经动力学及热力学分析发现只有一个单重态线性的分子11(O—C—S) 能够稳定存在.     

聚乙烯在羟基化β-石英(100)表面上的有序吸附

利用分子动力学模拟方法研究了聚乙烯链在羟基化β-石英(100)表面上的吸附. 结果表明, 吸附基底上的规则图案起到模板的作用, 无论是真空还是溶液环境, 都会引导聚乙烯链在表面形成二维沿着[110]方向取向的折叠构型. 为了与聚乙烯的吸附相比较, 进一步研究了聚氧化乙烯链在相同表面上的吸附情况. 结果表明, 极性链与非极性链在极性表面上的吸附情况完全不同.     

P14原子簇的理论研究

以Ｐ_２、Ｐ_４（Ｔｄ）、Ｐ_８（Ｃ_（２ｖ））及Ｐ_（１０）（Ｃ_（２ｖ））为结构单元设计了Ｐ_２＋Ｐ_４（Ｔ_ｄ）＋Ｐ_８（Ｃ_（２ｖ））→Ｐ_（１４）（Ｃ_ｓ）（Ａ）、Ｐｓ（Ｃ_（２ｖ）＋３Ｐ_２→Ｐ_（１４）（Ｃ_（２ｖ））（Ｉ）（Ｂ）、Ｐ_（１０）（Ｃ_（２ｖ）＋２Ｐ_２→Ｐ１４（Ｃ_（２ｖ））（Ⅱ）（Ｃ），利用从头算Ｇａｕｓｓｉａｎ－９４程序，选择６－３１Ｇ基组，对４种结构单元及Ｐ_（１４）原子簇的３种几何构型Ｃｓ、Ｃ_（２ｖ）（Ⅰ）、Ｃ_（２ｖ）（Ⅱ）进行全优化。相对能量的计算结果表明，Ｐ_（１４）（Ｃ_ｓ）与Ｐ_（１４）（Ｃ_（２ｖ））（Ⅰ）构型稳定。进一步设计（７／２）Ｐ_４（Ｔ_ｄ）→Ｐ_（１４）（Ｄ）及７Ｐ_２→Ｐ_（１４）（Ｅ），其相对能量表明，Ｐ_（１４）（Ｃ_ｓ）构型是稳定的．磷与磷连结单键键长范围为０．２２０～０．２２８ｎｍ，双键键长范围为０．２００～０．２０２ｎｍ，与实验结果大体相当。     

Direct ab initio dynamics study of radical C4H (X̃2Σ+) + CH4 reaction

The methane (CH(4)) hydrogen abstraction reaction by linear butadiynyl radical C(4)H (CCCCH) has been investigated by direct ab initio dynamics over a wide temperature range of 100-3000 K, theoretically. The potential energy surfaces (PESs) have been constructed at the CCSD(T)/aug-cc-pVTZ//BB1K/6-311G(d,p) levels of theory. Two different hydrogen abstraction channels by C(1) and C(4) of C(4)H (C(1)C(2)C(3)C(4)H) have been considered. The results indicate that the C(1) position of C(4)H is a more reactive site. The electron transfer behaviors of two possible channels are also analyzed by quasi-restricted orbital (QRO) in detail. The rate constants calculated by canonical variational transition-state theory (CVT) with the small-curvature tunneling correction (SCT) are in excellent agreement with available experimental values. The normal and three-parameter expressions of Arrhenius rate constants are also provided within 100-3000 K. It is expected to be helpful for further studies on the reaction dynamics behaviors over a wide temperature range where no experimental data is available so far.