用平均场理论研究盘状液晶分子间氢键对清亮点的影响

我们前期的合成和实验发现, 对于一系列含酰胺基支链的苯并菲分子, 其生成的液晶的清亮点的高低与分子间氢键的数目有关. 采用Gay-Berne势模拟液晶分子的双体势, 用量子化学和分子力学相结合的ONIOM (HF/STO-3G: UFF)方法优化含酰胺基支链的苯并菲分子的单体和双体的构型和能量, 从而获得Gay-Berne势的参数. 在此基础上, 根据平均场理论, 解一个关于序参量和温度的积分方程, 求出液晶的清亮点. 计算结果与实验一致, 可以解释清亮点和分子间氢键的关系.     

含乙烯氧基柔链的苯并菲盘状液晶清亮点的平均场理论研究

我们前期的合成和实验发现, 对于一系列含端烯氧基醚链的苯并菲分子, 它们没有分子间氢键, 其生成的液晶的清亮点的高低与端烯氧基醚链的数目有关. 本工作采用Gay-Berne势模拟液晶分子的双体势, 用量子化学和分子力学相结合的ONIOM(MP2/6-31G*(0.25):UFF)方法优化含端烯氧基醚链的苯并菲分子的单体和双体的构型和能量, 从而获得Gay-Berne势的参数. 在此基础上, 根据平均场理论, 解一个关于序参量和温度的积分方程, 求出液晶的清亮点. 计算结果与实验趋势相同, 可以解释清亮点和分子端烯氧基醚链数目的关系.     

紫杉醇的核磁共振谱及其分子几何构型的从头算研究

应用规范不变原子轨道GIAO法在HF/6-31G水平计算了MM3、PM3和HF/STO-3G三种最优分子构型下紫杉醇的δ1H-NMR和δ13C-NMR.对NMR化学位移计算值与实验值的对比和误差的统计和分析表明,采用HF/STO-3G优化的几何构型计算得到的化学位移误差最小,因此这一构型与实际分子最为接近.从理论构型得到4-OAc与3′-H的距离为0.2663nm,4-OAc与3′-Ph的距离为0.3104nm,4-OAc与2-OCOPh的距离为0.3707nm,支持了Williams等关于紫杉醇分子内存在NOE(nuclearoverhausereffect)效应及Velde等关于在极性溶剂中4-OAc与3′-Ph、2-OCOPh形成分子内疏水簇的观点.     

含苯并五元环方酸衍生物的二阶非线性光学性质的理论研究

利用从头算HF/6-31G*方法对用苯并五元环取代的方酸衍生物体系SQ1～SQ15进行几何构型优化和电子结构计算. 以优化后的构型为基础, 应用ZINDO方法计算电子光谱. 同时应用从头算CPHF/6-31G*方法和半经验FF/AM1, FF/PM3, FF/MNDO等有限场方法计算了分子的二阶非线性光学系数β_μ. 研究几何结构、电子结构和前线分子轨道能与β_μ之间的关系, 为设计性能优良的有机非线性光学材料提供理论指导, 并对这四种计算方法的结果进行了比较.     

笼状氮化铝团簇的结构特征

建立了一种设计笼状氮化铝团簇的原则与方法,以此设计了（AlN)_n (n = 2 ～ 4)团簇的几何结构。采用量子化学AM1,HF/STO-3G和HF/LANL2DZ方法,对设计 的所有分子的同分异构体进行了优化,筛选出了低能量构型。完成了最低能量构型 结合能和能量二次差分值的计算,找到了（AlN)_n团簇的稳定性规律。用密度泛函 理论（DFT）的B3LYP/6-31G方法,对（AlN)_n (n = 8 ～ 15)等8种团簇进行了验 证性的计算,所得结果与HF/LANL2DZ方法一致。为下一步在实验上合成较大尺寸（ AlN)_n团簇提供了理论依据。     

9-羟基苯并萘酮分子内氢传递过程的理论计算

借助MNDO、AM 1、HF/ 3 2 1G、HF/ 6 31G 量子化学理论计算方法 ,对 9 羟基苯并萘酮 (9 HPO)的分子内氢传递过程进行了理论探讨 ,并与X射线衍射的结构数据作了比较。发现 :1 .A构型的 9 HPO比B构型的更稳定 ,两个等同的A可以通过分子内氢传递相互转化 ,反应的过渡态就是B .2 .B是一种严格对称的平面构型 ,它的H电荷和偶极矩都较A增大 ,所以分子内氢传递速率在极性溶剂中将加快。3.共轭在 9 HPO分子内氢传递反应中起了重要的作用。 4.对于类似分子IHT的计算 ,先用半经验方法或小基组的从头计算方法优化结构 ,再用大基组的从头计算方法计算单点能可以得到较好的结果。     

甲氧基苯基偶氮衍生物电子结构与二阶NLO性质的HF/FF方法研究

采用量子化学HF方法在6-31G水平上优化6个甲氧基苯基偶氮衍生物分子的几何构型，利用HF／6-31G。方法计算它们的偶极矩、电荷分布、前线分子轨道能级并结合有限场(FF)方法计算二阶非线性光学系数．结果表明，偶氮苯衍生物分子具有很好的共轭性，在给吸电子基团作用下，电荷转移明显，展现示出较强的极性．偶氮苯衍生物分子与苯乙烯、Schiff碱类衍生物相似，也具有很好二阶非线性光学活性，同时六元杂环取代的偶氮苯衍生物分子二阶非线性光学系数比未取代的大，五元杂环取代结果相反．     

甲醛与乙醛,甲醚,硝基甲烷相互作用的从头算研究

用6-31G、全构型优化,研究了甲醛与乙醛、甲醚、硝基甲烷的相互作用。结果表明所有超分子稳定构型都包含2个C—H—O氢键的平面环状结构。H—O距离为0.228～0.264 nm,作用能为—19～—24 kJ/mol,与二聚水的作用能接近。稳定性主要取决于甲基上取代基Y吸电子能力以及环状结构中氢键张力。STO-3G不很适用于研究这类分子的相互作用。     

用平均场近似研究液晶分子结构与清亮点的关系

液晶的平均场理论假设每个液晶分子都是处在其它周围所有分子的平均作用势场中, 而不考虑该分子与周围各个分子之间的单独关联, 一个分子所感受到的分子场就是所有其它分子给予这个分子的作用场对位置和对取向的平均值. 当周围分子不是完全有序排列时, 其混乱程度可以用序参量S表示(S≤1), 中央一个分子受的势能为－Sg(θ). 其中g(θ)是在排列整齐的平均场中该分子的势能函数的负值. 在此基础上, 可以解一个关于序参量和温度的积分方程, 求出清亮点. 本文采用Gay-Berne势模拟液晶分子的双体势, 用分子力学UFF方法优化液晶分子(二环辛烷和硼笼的衍生物)的单体和双体的构型和能量, 从而获得Gay-Berne势的参数, 得到液晶的清亮点. 计算结果与实验基本一致.     

电场作用下分子导线的理论研究

摘要利用从头计算法分别在HF/6-31G, HF/6-31G*, HF/6-31G**, HF/6-31+G, HF/6-31++G, HF/6-31+G*, HF/6-31+G**, HF/D95+*, B3LYP/6-31G*和B3LYP/6-31+G*水平上计算了5个单体的聚乙炔分子导线, 从几何构型、 SCF能量和分子轨道能级三个方面讨论了外电场对分子导线的影响, 给出了聚乙炔分子导线性质与外电场变化的定量关系.     

Quantum chemical calculations based on ONIOM and the DFT methods in the inclusion complex: doxycycline/2-O-Me-β-cyclodextrin

A computational study of inclusion complexes of 2-methyl-βCD with Doxycycline tautomeric (enol and keto form) has been performed with several combinations of ONIOM hybrid calculations. The reliability of the ONIOM2 calculations at the integrated level, ONIOM2 (M05-2X/6-31G(d): M05-2X/3-21G*), ONIOM2 (M05-2X/6-31G(d):HF/3-21G*), ONIOM2 (B3LYP/6-31G(d):HF/3-21G*), ONIOM2 (B3LYP/6-31G(d):B3LYP/3-21G*) and ONIOM2 (B3PW91/6-31G(d):B3PW91/3-21G*) was examined. Their complexation, binding, deformation and stabilization energies, and geometrical data were compared with those of the target geometry structure optimized at the M05-2X/6-31G(d) level of theory. Mixed combinations ONIOM2 (M05-2X 6-31G(d):HF 3-21G*) and ONIOM2 (B3LYP 6-31G(d):HF 3-21G*) reproduces nearly the target geometry structure and provides realistic energetic results at a relatively low computational cost.     

Searching for conical intersections of potential energy surfaces with the ONIOM method: application to previtamin D

We demonstrate that the ONIOM method can be used to optimize a conical intersection between the ground and first excited-state potential energy surfaces of previtamin D (precalciferol), with excitation localized in a small part of the molecule: the hexatriene chromophore. These calculations were up to 100 times faster with little loss of accuracy compared to a full non-ONIOM Target calculation. The most accurate ONIOM method combination was CASSCF/4-31G//ROHF/STO-3G(Triplet): in comparison to the Target (CASSCF/4-31G), bond lengths and angles in the hexatriene model region were calculated to within 0.02 A and 0.7 degrees , respectively, and the energy difference between the conical intersection and nearest associated S 1 minimum to within 0.5 kcal x mol (-1). All of the low-level methods selected produced accurate geometries, including the UFF molecular mechanics and AM1 semiempirical methods, suggesting a cheap and efficient way of initially optimizing conical intersections geometries. Furthermore, ONIOM allows for an assessment of the localization of excited states, providing some fundamental insight into the physical processes involved.     

Imidazo[1,2-a]pyrazine-3,6-diones derived from alpha-amino acids: a theoretical mechanistic study of their formation via pyrolysis and silica-catalyzed process

Imidazo[1,2-a]pyrazine-3,6-diones are unusual compounds composed of three alpha-amino acid fragments. These bicyclic amidines (BCAs) form under high temperatures or with the use of strong dehydrating reagents. We gave insight into the mechanisms of BCA formation via gas-phase pyrolytic and silica-catalyzed reactions of glycine (Gly) and alpha-aminoisobutyric acid (AIB) with related diketopiperazines (DKPs), using quantum chemical calculations. The entire process requires four steps: (1) O-acylation of DKP with free or silica-bonded amino acid, (2) acyl transfer from the oxygen to the nitrogen atom, (3) intramolecular condensation of the N-acyl DKP into a cyclol, and (4) elimination of water. To study step (1) at silica surface (modeled by H7Si8O12-OH cluster), we employed two-level ONIOM calculations (AM1:UFF, B3LYP/3-21G:UFF and B3LYP/6-31G(d):UFF); all gas-phase reactions were studied at the AM1, B3LYP/3-21G and B3LYP/6-31G(d) levels. The catalytic effect of silica was observed for both Gly and AIB: the activation energy in the O-acylation at the surface was lower by more than 9 kcal mol(-1) as compared to the gas-phase process. Contrary to the exothermic O-acylation, the gas-phase transfer reaction (step 2) was exothermic in both cases, but more favorable for Gly. The cyclocondensation of N-acylated DKPs into BCAs (steps 3 and 4) is endothermic for Gly and exothermic for AIB.     

Experimental and theoretical investigation of the first overtone spectrum of 1,3,5-trinitrotoluene

Vibrational overtone spectroscopy is a powerful tool for studying intramolecular and intermolecular interactions. We report on a combined experimental and modeling study of the C-H stretch first overtone of bulk 1,3,5-trinitrotoluene (TNT) and TNT on fumed-silica powder. We recorded the overtone spectra by laser photoacoustic spectroscopy and compared them with those predicted with the harmonically coupled anharmonic oscillator model in the 5600-6600 cm(-1) region. The model systems included single molecules and hybrid quantum and molecular mechanical (QM:MM) clusters to account for the effects of intermolecular interactions on the observed spectra. We performed the hybrid QM:MM calculations at the HF/6-31+G(d,p), B3LYP/6-31+G(d,p), and MP2/6-31+G(d,p) levels of theory and with the universal force field (UFF) to account for van der Waals and electrostatic effects from surrounding molecules. Overtone spectra calculated from the MP2 level of theory, using a HF/3-21+G* calculation to assign molecular charges in the MM layer, and the Merz-Singh-Kollman population analysis for assigning partial charge in the QM layer and determining the transition dipole moment agreed best with the experimental data.     

A QM/MM-based computational investigation on the catalytic mechanism of saccharopine reductase

Saccharopine reductase from Magnaporthe grisea, an NADPH-containing enzyme in the α-aminoadipate pathway, catalyses the formation of saccharopine, a precursor to L-lysine, from the substrates glutamate and α-aminoadipate-δ-semialdehyde. Its catalytic mechanism has been investigated using quantum mechanics/molecular mechanics (QM/MM) ONIOM-based approaches. In particular, the overall catalytic pathway has been elucidated and the effects of electron correlation and the anisotropic polar protein environment have been examined via the use of the ONIOM(HF/6-31G(d):AMBER94) and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) methods within the mechanical embedding formulism and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) and ONIOM(MP2/6-311G(d,p)//HF/6-31G(d):AMBER94) within the electronic embedding formulism. The results of the present study suggest that saccharopine reductase utilises a substrate-assisted catalytic pathway in which acid/base groups within the cosubstrates themselves facilitate the mechanistically required proton transfers. Thus, the enzyme appears to act most likely by binding the three required reactant molecules glutamate, α-aminoadipate-δ-semialdehyde and NADPH in a manner and polar environment conducive to reaction.     

A theoretical study on the structure and properties of phenothiazine derivatives and their radical cations

Semiempirical CNDO, AM1, PM3 and ab initio HF/STO-3G, HF/3-21G(d), and HF/6-31(d) methods were employed in the geometry optimization of the phenothiazine and the corresponding radical cation. The results obtained from the PM3 performances were as good as those from the ab initio calculations in the structure optimization of both phenothiazine and phenothiazine radical cation. The PM3 method was used to optimize the structures of a series of N-substituted phenothiazine derivatives and their radical cations. The PM3-optimized results were then analyzed with the ab initio calculation at the 6-311G(d,p) level, which yielded the total energy, frontier molecular orbitals, dipole moments, and charge and spin density distributions of the phenothiazine derivatives and their radical cations.     

Intermolecular Interaction of HMX： an Application of ONIOM Methodology

IntroductionExtensive studies on intermolecular interac-tions have been made in the past decades due totheir importance in a wide range of physical,chem-ical and biological fields.Researches on the weakintermolecular interactions began with hydrogenbonds.Scheiner summarized the ab initio investiga-tions on hydrogen bonding in detail[1] .With the abinitio method supermolecular structures and bind-ing energies can be predicted notonly for H- bondedsystems but also for other systems[2 _ 5] when …     

Theoretical study of vitamin properties from combined QM-MM methods: Comparison of chemical shifts and energy

The combination of quantum mechanics (QM) and molecular mechanics (MM) methods has become an alternative tool for many applications for which pure QM and MM are not suitable. The QM-MM method has been used for different types of problems, for example, structural biology, surface phenomena, and the liquid phase. In this paper, we have implemented these methods for vitamins, an important kind of biological molecule, and then compared results. The calculations were done by the full ab initio method (HF/3–21 g and HF/6–31 g) and QM-MM (ONIOM) method with HF(3–21 g)/AM1/UFF; then, we found that the geometry obtained by the QM-MM method is very accurate and this rapid method can be used in place of time consuming ab initio methods for large molecules. A comparison of energy values in the QM-MM and QM methods is given. We compare chemical shifts and conclude that the QM-MM method is a perturbed full QM method. The text was submitted by the authors in English.     

Bisaziridine tetracarboxylates as building blocks in the stereoselective synthesis of C60-fullerene diads and dumbbell-like bis-C60-fullerene triads

The synthesis of alkoxycarbonyl-substituted bisaziridines with the two aziridine units connected by conjugated p-phenylene, partly conjugated 1,1'-biphenyl-4,4'-diyl, and nonconjugated 4,4'-methylenediphenyl linkers was developed. The reaction of azomethine ylides derived from the bisaziridines with fullerene C(60) was optimized and used for the stereoselective preparation of both the monoadducts (C(60)-linker-aziridine dicarboxylate), and the dumbbell bisadducts (C(60)-linker-C(60)). The reasons for the observed selectivity of the azomethine ylide formation and cycloaddition were theoretically studied at the DFT B3LYP/6-31G(d) level or at the ONIOM B3LYP/6-31G(d):B3LYP/STO-3G level for fullerene-containing molecules.     

Adsorption of NH(3) and H(2)O in acidic chabazite. Comparison of ONIOM approach with periodic calculations

The adsorption of NH(3) and H(2)O in acidic chabazite has been studied with the B3LYP method within the cluster approach (5T, 48T clusters) and the periodic approach adopting a Si/Al = 11/1 chabazite and a basis set of polarized double-zeta quality. The 5T cluster has been treated fully ab initio at the B3LYP level whereas the 48T cluster has been treated with the ONIOM2 scheme using B3LYP as the high level of theory and the MNDO, AM1, and HF/3-21G methods as low levels of theory. Periodic calculations show that the adsorption of NH(3) in acidic chabazite takes place through an ion pair (NH(4)(+)-CHA(-)) structure, the computed adsorption energy being -32 kcal/mol. The adsorption of H(2)O leads to a hydrogen bonded (H(2)O-HCHA) complex with the computed adsorption energy of -20 kcal/mol. All ONIOM combinations provide similar structures to those obtained with periodic calculations. Adsorption energies, however, are sensitive to the low level used, especially for NH(3). The ONIOM B3LYP:HF/3-21G method is the one that provides more satisfactory results. Present results show that, for larger zeolites, the ONIOM scheme can be successfully applied while drastically reducing the cost of a fully ab initio treatment.