LDHs主体层板与卤素阴离子超分子作用的理论研究

构建了LDHs主客体作用模型, 采用混合密度泛函B3LYP方法, 在6-31G(d)水平上进行结构优化和频率分析, 然后分别用6-31G(d)和6-311++G(d, p)基组计算主客体相互作用能, 从几何参数、电荷布居、前线轨道、能量以及热力学参数等角度探讨LDHs主体层板与卤素阴离子(F?, Cl?)间的超分子作用. 计算结果表明, LDHs主体层板复合卤素阴离子是一个自发过程. LDHs主客体间存在着较强的超分子作用, 主要包括静电和氢键作用, 相互作用能分别为?592.45和?444.01 kJ·mol?1. LDHs主体层板与卤素阴离子的前线轨道发生作用, 电子容易从卤素阴离子的HOMO向层板的LUMO转移, 形成的组装产物Mg6Al(OH)14+?F?比Mg6Al(OH)14+?Cl?稳定.     

类水滑石主体层板与客体CO2-3、H2O间的超分子作用

通过构建类水滑石双层计算模型, 采用混合密度泛函B3LYP/6-31G(d)//B3LYP/3-21G方法计算类水滑石(LDHs-CO3-yH2O)的结构与能量, 探讨LDHs限域空间中客体阴离子及水分子的分布形态以及主客体超分子作用. 计算结果表明, 客体阴离子与水分子以平行层板的方式存在于水滑石层间. 主客体发生作用时, CO2-3的HOMO轨道向层板的LUMO轨道转移电子. 所形成的LDHs-CO3主客体作用要强于LDHs-F以及LDHs-Cl, 与其离子交换性能相一致. 水滑石去水结构(LDHs-CO3)水合过程, 氢键作用较静电作用更占优势, 并且layer-water型氢键要强于anion-water 型氢键. 此外, 水合能计算表明LDHs水合具有一定的饱和量.     

类水滑石主体层板与客体CO2-3、H2O间的超分子作用

通过构建类水滑石双层计算模型,采用混合密度泛函B3LYP/6-31G(d)//B3LYP/3-21G方法计算类水滑石(LDHs-CO3-yH2O)的结构与能量,探讨LDHs限域空间中客体阴离子及水分子的分布形态以及主客体超分子作用.计算结果表明,客体阴离子与水分子以平行层板的方式存在于水滑石层间.主客体发生作用时,CO2-3的HOMO轨道向层板的LUMO轨道转移电子.所形成的LDHs-CO3主客体作用要强于LDHs-F以及LDHs-Cl.与其离子交换性能相一致.水滑石去水结构(LDHs-CO3)水合过程,氢键作用较静电作用更占优势,并且layer-water型氢键要强于anion-water型氢键.此外,水合能计算表明LDHs水合具有一定的饱和量.     

密度泛函理论研究间硝基苯基吡咯酰胺识别卤素阴离子

采用混合密度泛函B3LYP方法, 在LANL2DZ水平上优化间硝基苯基吡咯酰胺(NPC)及其卤素阴离子复合物的几何构型, 从几何结构参数、电荷布居、前线轨道、结合能以及热力学参数等角度探讨复合物形成过程中主体分子的构象变化以及主客体间的超分子作用. 研究结果表明, 吡咯酰胺主体分子识别卤素阴离子是一个自发过程. 识别过程包括阴离子诱导下主体分子构象转变和“活性”构象与阴离子形成氢键两个步骤. 计算表明主体分子复合阴离子过程中, 两者的前线轨道发生作用, 电子容易从卤素阴离子的HOMO向主体分子的LUMO转移. 在所形成的四种复合物NPC-X^－ (X^－＝F^－, Cl^－, Br^－, I^－)中, 以NPC-F^－最稳定, 与文献报道吡咯酰胺主体分子识别卤素阴离子能力的实验结果相一致.     

一种天然产物Wangzaozin A的细胞毒活性

从总序香茶菜Isodon racemosa (Hemsl) Hara植物中分离得到一个对人类肿瘤细胞Bel-7402和HD-8910具有毒活性的对映-贝壳衫烷型二萜Wangzaozin A化合物(1). 应用密度泛函理论(DFT)B3LYP方法, 对该分子的几何构型进行优化, 结果表明用B3LYP/6-31G(d)优化的几何参数与它的X射线衍射结构参数基本一致. 在优化的几何构型基础上, 采用规范不变原子轨道(GIAO)法, 在B3LYP理论水平分别用6-31G(d), 6-31G(d,p), 6-31+G(d,p)和6-31++G(d,p)基组进行核磁共振(NMR)化学位移值计算, 预测的1H和13CNMR化学位移值与实验值吻合; 统计误差分析表明, 用B3LYP/6-31G(d)优化的分子构型接近实际的分子构型. 因此, DFT方法适用这一类型化合物的构型和NMR参数进行预测. 在几何优化的基础上, 在B3LYP/6-31G(d)水平上, 对Wangzaozin A分子的静电位(MEP)进行理论计算. MEP三维图表明, 在Wangzaozin A分子中α-亚甲基环戊酮的羰基和羟基附近出现富电子区域(负电位), 起着供电子作用, 与受体的正电子区域结合. 这些结果从理论上支持了α-亚甲基环戊酮结构是一种抗肿瘤活性中心的看法.     

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

摘要利用从头计算法分别在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能量和分子轨道能级三个方面讨论了外电场对分子导线的影响, 给出了聚乙炔分子导线性质与外电场变化的定量关系.     

吡唑啉-噁二唑类有机分子的电子光谱和二阶非线性光学性质

使用密度泛函理论在B3LYP/6-31＋G＊水平上对吡唑啉-噁二唑类（A1～A10）10个分子进行几何构型优化,在此基础上进行自然界轨道电荷分析,采用TD-DFT（TDB3LYP/6-31＋G＊）计算电子吸收光谱,有限场FF方法（B3LYP/6-31＋＋G＊＊）计算二阶非线性光学性质？0.计算结果表明,此类分子的？0值...     

密度泛函理论研究DMSO溶剂中邻-二(吡咯-2-甲酰胺基)亚苯与卤素阴离子间的超分子作用

采用混合密度泛函B3LYP方法, 在LANL2DZ水平上优化DMSO溶剂中邻-二(吡咯-2-甲酰胺基)亚苯(PFP)及其卤素阴离子复合物的几何构型, 从几何结构参数、电荷布居、前线轨道、结合能以及热力学参数等角度探讨复合物形成过程中主体分子的构象变化以及主客体间的超分子作用. 研究结果表明, 邻-二(吡咯-2-甲酰胺基)亚苯识别卤素阴离子是一个自发过程. 复合阴离子过程中, 两者的前线轨道发生作用, 电子容易从卤素阴离子的HOMO向主体分子的LUMO转移. 在所形成的四种复合物PFP:X^－(X^－＝F^－, Cl^－, Br^－, I^－)中, 以PFP:F^－最稳定, 与文献报道邻-二(吡咯-2-甲酰胺基)亚苯识别卤素阴离子能力的实验结果相一致. 设计合成含多个吡咯酰胺基团、并能形成合适“杯口”大小的主体化合物, 将增强对卤素阴离子的识别能力.     

已烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p),6-31+G(d,p),6-311G(d,p)及6- 311+G(d,p)基组,分别对1-C_6H_(12)~+,2-C_6H_(12)~+和3-C_6H_(12)~+的各种构 象进行了几何构型优化,并在B3LYP/6-311G(d,p)水平上进行了频率分析计算,在 各优化构型上,使用B3LYP和MP2(full)方法进行了超精细结构的计算。计算的3- C_6H_(12)~+的超精细偶合常数比以往的计算结果更好;1-C_6H_(12)~+和2-C_6H_ (12)~+的超精细偶合常数目前尚无实验数据报道,本计算预言了它们的超精细偶合 常数和最稳定构型。     

一种室温离子液体与水相互作用的密度泛函理论研究

采用量子化学计算中的密度泛函理论(DFT),在B3LYP/6-31+G(d,p)计算水平上研究了离子液体1-乙基-3-甲基咪唑四氟硼酸盐([Emim]BF4)及其与水分子形成的复合物的稳定构型和相互作用能;经振动频率分析得到了[Emim]BF4及其与水的复合物的红外光谱.计算结果表明,相对于水分子与阳离子的作用而言,水分子与阴离子的作用对离子液体结构的影响更大.与此同时,实验测得的[Emim]BF4的红外光谱与计算结果吻合.     

Theoretical study of thiourea derivatives as chemosensors for fluoride and acetate anions

The different interactions between a chemosensor,1-(naphthalen-2-yl)-3-(6-nitrobenzothiazol-2-yl)-thiourea(1),and F-,acetate(AcO-),Cl-,and Br-anions have been investigated theoretically at the B3LYP/6-31+G(d,p) level with the basis set superposition error(BSSE) correction.It was found that the high selectivity of compound 1 for F-can be ascribed to the ability of the anion to deprotonate the N-H fragment of the host sensor,while the chemosensor also has a strong affinity for AcO-by virtue of the formation of a hydrogen-bonded complex.Intramolecular charge transfer(ICT) causes the colorimetric signaling of compound 1 after interaction with F-/AcO-.A study of substituent effects suggested that the O/NH-and O/S-substituted derivatives are also expected to be promising candidates for chromogenic F-/AcO-chemosensors.     

Microhydration of the guanine-cytosine (GC) base pair in the neutral and anionic radical states: a density functional study

A density functional study of the effects of microhydration on the guanine-cytosine (GC) base pair and its anion radical is presented. Geometries of the GC base pair in the presence of 6 and 11 water molecules were fully optimized in the neutral (GC-nH2O) and anion radical [(GC-nH2O)*-] (n = 6 and 11) states using the B3LYP method and the 6-31+G** basis set. Further, vibrational frequency analysis at the same level of theory (B3LYP/6-31+G**) was also performed to ensure the existence of local minima in these hydrated structures. It was found that water molecules surrounding the GC base pair have significant effects on the geometry of the GC base pair and promote nonplanarity in the GC base pair. The calculated structures were found to be in good agreement with those observed experimentally and obtained in molecular dynamics (MD) simulation studies. The water molecules in neutral GC-nH2O complexes lie near the ring plane of the GC base pair where they undergo hydrogen bonding with both GC and each other. However, in the GC anion radical complexes (GC-nH2O, n = 6, 11), the water molecules are displaced substantially from the GC ring plane. For GC-11H2O*-, a water molecule is hydrogen-bonded with the C6 atom of the cytosine base. We found that the hydration shell initially destabilizes the GC base pair toward electron capture as a transient anion. Energetically unstable diffuse states in the hydration shell are suggested to provide an intermediate state for the excess electron before molecular reorganization of the water molecules and the base pair results in a stable anion formation. The singly occupied molecular orbital (SOMO) in the anion radical complexes clearly shows that an excess electron localizes into a pi orbital of cytosine. The zero-point-energy (ZPE-) corrected adiabatic electron affinities (AEAs) of the GC-6H2O and GC-11H2O complexes, at the B3LYP/6-31+G** level of theory, were found to be 0.74 and 0.95 eV, respectively. However, the incorporation of bulk water as a solvent using the polarized continuum model (PCM) increases the EAs of these complexes to 1.77 eV.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

Comparison of Some Computational Methods for Geometry Optimisation of Phosphorus Acid Derivatives

Several economical methods for geometry optimization, that should be applicable to larger molecules, have been evaluated for 19 phosphorus acid derivatives. MP2/cc-pVDZ geometry optimizations are used as reference points and the geometries obtained from the other methods are evaluated with respect to deviations in bond lengths and angles, from the reference geometries. The geometry optimization methods are also compared to the much used B3LYP/6-31G(d) method. Single point energies obtained by subsequent EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) calculations on the respective equilibrium geometries are also reported relative to the energies obtained from the reference geometries. The geometries from HF/MIDI! optimizations were closer to those of the references than the geometries of the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! optimizations. The EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) single point energies obtained from the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! geometries gave a mean absolute deviation (MAD) from that of the reference geometries of 1.4-3.9 kcal mol m 1 . The HF/MIDI! geometries, however, gave EDF1/6-31+G(d) and B3LYP/6-31+G(d,p) energies with a MAD of only about 0.5 and 0.55 kcal mol m 1 respectively from the energies obtained with the reference geometries. Thus, use of HF/MIDI! for geometry optimization of phosphorus acids is a method that gives geometries of near-MP2 quality, resulting in a fair accuracy of energies in subsequent single point calculations, at a much lower computational cost other methods that give similar accuracies.     

Conformational study of the structure of free 18-crown-6

A conformational search was performed for 18-crown-6 using the CONLEX method at the MM3 level. To have a more accurate energy order of the predicted conformations, the predicted conformations were geometry optimized at the HF/STO-3G level and the 198 lowest energy conformations, according to the HF/STO-3G energy order, were geometry optimized at the HF/6-31+G level. In addition, the 47 nonredundant lowest energy conformations, according to the MP2/6-31+G energy order at the HF/6-31+G optimized geometry, hereafter the MP2/6-31+G//HF/6-31+G energy order, were geometry optimized at the B3LYP/6-31+G level. According to the MP2/6-31+G//B3LYP/6-31+G energy order, three conformations had energies lower than the experimentally known Ci conformation of 18c6. At the MP2/6-31+G//B3LYP/6-31+G level, the S6 lowest energy conformation is more stable by 1.96 kcal/mol than this Ci conformation. This was confirmed by results at the MP2/6-31+G level with an energy difference of 1.84 kcal/mol. Comparison between the structure of the S6 conformation of 18c6 and the S4 lowest energy conformation of 12-crown-4, as well as other important conformations of both molecules, is made. It is concluded that the correlation energy is necessary to have an accurate energy order of the predicted conformations. A rationalization of the conformational energy order in terms of the hydrogen bonding and conformational dihedral angles is given. It is also suggested that to have a better energy order of the predicted conformations at the MM3 level, better empirical force fields corresponding to the hydrogen bond interactions are needed.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

A comparative effect of ring annelation on cation-benzene and anion-benzene

DFT/B3LYP calculations were carried out on several π-complexes formed by cations and anions with annelated benzene, respectively. The binding energies obtained with standard method were corrected by basis set superposition error (BSSE) and zero-point energy (ZPE) during the geometry optimization for all complexes at the same levels of theory, respectively. Some different aspects of the π–cation have been compared to those of π–anion, involving in binding energy changes in effect of ring annelation, the aromaticity of the ring upon complexation, Mulliken and NBO charge-transfer. The effect of BSSE correction during the optimization is very important in some π–anion complexes whether or not using diffuse functions in basis set, and results with at least one set of diffuse functions 6-31+G(d) basis set is a little better than results obtained by 6-31G(d, p) basis set for some π–anion especially for F⁻ complexes.     

Heat of hydrogenation of 1,5-dehydroquadricyclane. A computational and experimental study of a highly pyramidalized alkene

The radical anion of the highly pyramidalized alkene 1,5-dehydroquadricyclane (1) was generated in the gas phase from the Squires reaction of 1,5-bis(trimethylsilyl)quadricyclane with F-/F2. The electron binding energy and proton affinity of 1*- were determined by bracketing experiments to be 0.6 +/- 0.1 eV and 386 +/- 5 kcal/mol, respectively. These values are in good agreement with values predicted by density functional theory (B3LYP/6-31+G*) and ab initio (CASPT2/6-31+G*) calculations. The experimental heat of hydrogenation of 1, obtained from a thermochemical cycle, was found to be 91 +/- 9 kcal/mol. This value of deltaH(H2) leads to values of 67 +/- 9 kcal/mol for the olefin strain energy (OSE) of 1, 172 +/- 9 kcal/mol for its heat of formation, and 23 +/- 9 kcal/mol for its pi bond dissociation enthalpy. Since the retro-Diels-Alder reaction of neutral 1 is computed to be highly exothermic, the finding that 1*- apparently does not undergo a retro-Diels-Alder reaction is of particular interest. The B3LYP/6-31+G* optimized geometry of 1 suggests that the bonding in this alkene is partially delocalized, presumably because the highly pyramidalized double bond in 1 interacts with the distal cyclopropane bonds in a manner that eventually leads to a retro-Diels-Alder reaction. The good agreement of the B3LYP and (2/2)CASPT2 values for the heat of hydrogenation and OSE of 1 with the experimentally derived values provides indirect evidence for the correctness of the B3LYP prediction that the equilibrium geometry of 1 lies part way along the reaction coordinate to the transition structure for the retro-Diels-Alder reaction.