有机物对蝴蝶幼虫麻醉活性与结构参数的关系

应用量子化学HF/3-21G从头算法得到了20种有机分子的优势构象,利用HF/3-21G和分子图形学技术获得其电子结构参数和几何结构参数,并将这些参数与有机物对蝴蝶幼虫的麻醉活性相关联.结果表明,有机物麻醉活性与分子的范德华体积、分子最低空轨道能和分子偶极矩之间存在良好的多元线性相关性,成功地建立有机物对蝴蝶幼虫麻醉活性的构效关系式.     

脂肪醇盐离子质子亲和能与量子化学参数的关系

采用半经验量子化学方法（MNDO，AM1和PM3）、从头算（采用STO-3G，3-21G，6-31G，6-31G^*,^6-31G^**,6-311G,6-311G^*和6-311G^**基组）HF方法计算羟基和6种脂肪醇盐离子以获得电子结构数据。通过相关分析发现脂肪醇盐离子中的氧原子净电荷，HOMO能级与质子亲和能之间具有良好的线性的相关性。     

原子类型电拓扑指数在脂肪醇类化合物毒性预测研究中的应用

采用原子类型电拓扑指数(ETSIAT)为结构描述子, 对梨形四膜虫、番茄、红蜘蛛、发光细菌荧光素酶、呆鲦鱼以及鼠的脂肪醇毒性进行了定量构效关系(QSAR)研究. 在偏最小二乘法建立QSAR模型的基础上, 采用内部和外部两种验证方法对模型进行验证. 研究结果表明: ETSIAT可以很好地表征脂肪醇与多种生物毒性相关的结构信息, 所建模型具有良好的稳定性和外部预测能力. 偏最小二乘模型分析结果显示: 疏水性可能是影响脂肪醇毒性大小的主要因素, 其毒性随C链长度的增加而增强, 且相同C原子数的直链脂肪醇毒性明显高于支链脂肪醇; 由于空间位阻效应, 长链脂肪醇与短链脂肪醇可能存在不同的毒性作用机理.     

预测氯代芳烃对戈卑鱼毒性的理论研究

为建立氯代芳烃对戈卑鱼毒性相关的QSAR模型,分析了22种氯代芳烃的结构特征,计算了各个分子的分子连接性指数.通过多元线性逐步回归分析研究,筛选了其中4种分子连接性指数~0X,~2X,~3X和~4X_(PC),建立了这4种指数与氯代芳烃对戈卑鱼毒性的定量结构-活性相关的QSAR模型,方程的相关系数R=0.979,判定系数R_(adj)~2=0.959,标准误差S=0.134,经检验该模型具有良好的稳定性和预测能力.     

甲氧基多溴联苯醚色谱相对保留时间的QSRR研究

基于有机化合物定量结构与色谱相对保留时间的关系,对26个甲氧基多溴联苯醚化合物进行了B3LYP/6-31G*水平上的结构优化,并在优化结构基础上进行了分子连接性指数和量子化学结构参数的提取。应用SPSS统计软件对甲氧基多溴联苯醚化合物的色谱相对保留时间与分子连接性指数进行了多元线性回归(MLR)。结果表明,三阶类分子连接性指数可以很好地表达甲氧基多溴联苯醚化合物色谱相对保留时间与结构描述符之间的定量关系。用外部测试集的方法对模型预测能力进行了验证,结果表明:所建立的QSRR模型具有较强的稳健性和预测能力。     

脂肪醇气相色谱保留指数与结构的相关性研究

在分子图的邻接矩阵基础上,构建了一个化合物均价连接性指数mL,mL=∑(Ai·Aj·Ak…)-0.5,其中一阶指数1L及定位基参数β与25种脂肪醇在6种固定相(SE-30,OV-3,OV-7,OV-11,OV-17和OV-25)上的气相色谱保留指数I显著相关,相关系数均大于0.98。所建定量结 构-保留关系(QSRR)模型具有良好的稳定性和预测能力,较好地揭示了脂肪醇在不同固定相上气相色谱保留指数的变化规律。     

一个新的连接性指数对脂肪醇的QSPR/QSAR研究

在分子拓扑理论的基础上 ,提出了一个新的连接性指数nχz,并用 nχz 研究了脂肪醇化合物在水中的溶解度 (Sw)、正辛醇 /水分配系数 (Kow)及对水生生物急性毒性 (-LC50 ) ,给出了相关方程 ,研究结果表明 ,相关性良好 ,新方法计算方便 ,物理意义明确 ,预测值与相应的实验值较吻合     

分子连接性指数与溶液吸附自由能的相关性

测定了硅胶自四氯化碳中吸附脂肪醇及活性炭自水中吸附芳香化合物的吸附等温线, 等温线均可用Langmuir方程描述。利用Langmuir参数计算了各体系的吸附标准自由能变化。计算了各种吸附质的分子连接性指数。各种分子的分子连接性指数与吸附标准自由能变化间有相关性, 并给出了相关方程, 同时作了初步的讨论。     

拓扑指数与烷烃色谱保留指数的定量相关研究

基于分子图论, 计算了64个链烷烃的一阶连接性指数(¹χ)和奇偶指数(OEI), 并对气相色谱保留指数进行定量结构-性质相关(QSPR)研究, 得到多元线性回归模型RI＝55.7631＋25.7317OEI＋111.4508¹χ (R²＝0.9960, RMS＝9.36, ARD＝1.02%). 交叉验证和预测结果表明, 所建立的QSPR模型具有良好的稳定性和预测能力. 与文献结果比较, 本文所用的参数少, 且计算简便.     

非金属二元氢化物pka与量子化学参数的关系

在HF／3－21G水平下，对13种非金属二元氢化物分子（HF，HCl,HBr,HI,H2O,H2S,H2Se,H2Te,NH3,PH3,AsH3,CH4和SiH4)进行了几何构型全优化和电子结构计算，将获得的电子结构数据与这些分子的pka值相关联，进行逐步回归分析，结果显示，pka与非氢原子的键价，分子的最低空轨道能间存在显的二元线性相关性，相关系数R＝0．994。pka计算值与实验值符合得较好，表明非金属二元氢化物中非氢原子的键价，分子的最低空轨道能在决定其酸性上起重要作用。     

在小硅化物SIX_4~m（X＝H，He，Li，Be，O，F，Ne，Na，Mg，S，Cl，Ar；m＝0，±4）中分子成键倾向的从头算研究

用从头算方法，在ＨＦ／ＳＴＯ－３Ｇ、ＨＦ／３－２１Ｇ和ＨＦ／６－３１Ｇ水平上研究了小硅化物的成键倾向性。计算结果表明，所研究的分子势能曲线均有稳定的极小值（ＳｉＬｉ_４除外），与已知的稳定分子ＳｉＨ_４、ＳｉＦ_４和ＳｉＣｌ_４比较，含惰性元素的未知分子ＳｉＨｅ、ＳｉＮｅ和ＳｉＡｒ比含碱金属和碱土金属的未知分子ＳｉＬｉ_４、ＳｉＮａ_４、ＳｉＢｅ和ＳｉＭｇ有较明显的成键可能性，而未知分子ＳｉＳ介于其间，所以ＳｉＨｅ、ＳｉＮｅ和ＳｉＡｒ有可能成为稳定的新分子的候选对象。     

A comparative molecular field analysis (CoMFA) study using semiempirical, density functional, ab initio methods and pharmacophore derivation using DISCOtech on sigma 1 ligands

The Comparative Molecular Field Analysis (CoMFA) was developed to investigate a three-dimensional quantitative structure activity relationship (3D-QSAR) model of ligands for the sigma 1 receptor. The starting geometry of sigma-1 receptor ligands was obtained from the Tripos force field minimizations and conformations were decided from DISCOtech using the SYBYL 6.8. program. The structures of 48 molecules were fully optimized at the ab initio HF/3-21G* and semiempirical AM1 calculations using GAUSSIAN 98. The electrostatic charges were calculated using several methods such as semiempirical AM1, density functional B3LYP/3-21G*, and ab initio HF/3-21G*, MP2/3-21G* calculations within GAUSSIAN 98. Using the optimized geometries, the CoMFA results derived from the HF/3-21G method were better than those from AM1. The best CoMFA was obtained from HF/3-21G* optimized geometry and charges (R2 = 0.977). Using the optimized geometries, the CoMFA results derived from the HF/3-21G methods were better than those from AM1 calculations. The training set of 43 molecules gave higher R2 (0.989-0.977) from HF/3-21G* optimized geometries than R2 (0.966-0.911) values from AM1 optimized geometries. The test set of five molecules also suggested that HF/3-21G* optimized geometries produced good CoMFA models to predict bioactivity of sigma 1 receptor ligands but AM1 optimized geometries failed to predict reasonable bioactivity of sigma 1 receptor ligands using different calculations for atomic charges.     

Synthesis,Crystal Structure and Theoretical Calculation of [2-（2,6-Dioxacyclohexyl）-5-methoxylphenols]2

The title compound [2-（2,6-dioxacyclohexyl）-5-methoxylphenols]2 was synthesized by the reaction of 2-hydroxyl-3-methoxylbenzaldehyde and 1,3-propanediol in the presence of DMF-DMS adduct and characterized by IR spectrum, UV-Vis spectrum and X-ray crystallography. The title compound belongs to monoclinic, space group P21 with a = 9.8967（10）, b = 8.2240（9）, c = 13.3654（14） A^°, β= 90.016（2）°, C11H14O4, Mr = 210.22, V = 1087.8（2）A^°3, Z = 4, Dc = 1.284 Mg/m^3, F（000） = 448,μ = 0.098 mm^-1, the final R = 0.0300 and wR = 0.0761 for 2070 observed reflections with I 〉 2σ（I）. The molecules are connected via intermolecular O-H…O hydrogen bonds into a 2D network structure. Moreover, theoretical investigations of the title compound with HF/6-31G（d） method were performed, and its stability, frontier molecular orbital composition and Mulliken charge distribution were also discussed. The compound is a bis-molecule. The two molecules stay together and could not be separated. Two crystallographically independent molecules exist in an asymmetric unit. The bond lengths and bond angles of the two molecules are slightly different with each other.     

Mixing parameters for geometry optimization using the Hamiltonian algorithm

We study the mixing parameters for the search of an optimal geometry using the Hamiltonian algorithm (HA) combined with ab initio molecular orbital calculations. We choose the C?CC?CC?CC dihedral angle of the butane molecule as an example. HF/3-21G level calculations are employed as the molecular orbital calculations. The distributions of the eigenvalues of mixing coefficients are fitted with the linear, quadratic, and quartic functions. Analyses of HA calculations both up to 2,000 and 60,000 iterative calculations show a possibility that the mixing process reduces the number of iterations. The low energy HF/3-21G, B3LYP/6-31G**, and PCM B3LYP/aug-cc-pVDZ optimized structures of the N-acetyl l-histidine N??-methyl amide and four water molecule supermolecule were also determined using the HA optimization method and compared to the recently determined thought to be global minimum energy structure.     

Vibrational spectra and ab initio analysis of tert-butyl, trimethylsilyl, trimethylgermyl, and trimethylstannyl derivatives of 3,3-dimethylcyclopropene VIII. 3,3-Dimethyl-1,2-bis-(trimethylstannyl)cyclopropene

The quantum mechanical force fields of 3,3-dimethyl-1,2-bis-(tert-butyl)cyclopropene (I), 3,3-dimethyl-1,2-bis-(trimethylsilyl)cyclopropene (II), 3,3-dimethyl-1,2-bis-(trimethylgermyl)cyclopropene (III), and 3,3-dimethyl-1,2-bis-(trimethylstannyl)cyclopropene (IV) were calculated at the HF/3-21G*//HF/3-21G* level. The scale factors which were optimized previously for the HF/3-21G*//HF/3-21G* quantum mechanical force field of 3,3-dimethyl-1-(trimethylsilyl)cyclopropene were used for correction of the force fields of these molecules. Good agreement between the frequencies calculated from these scaled force fields and the well-analyzed and assigned experimental frequencies of II and III suggests the transferability of these scale factors and the possibility of the spectroscopically accurate prediction of the vibrational spectrum of IV. Some regularities in the changes of the vibrational frequencies were found for this molecular series.     

Molecular Dipole Moment Computed with Ab Initio MKS Charges

1 INTRODUCTION Dipole moment (DM) can be measured by a variety of experimental methods[1] or computed with atomic charge distribution directly derived from molecular orbital calculations[2~10]. The average errors in DM computed with AM1, PM3 and MNDO methods for 125 organic compounds are 0.35 D, 0.45 D and 0.45 D, respectively[2]. Feller et al. reported that the average errors in the computed DM at the HF level are 0.65 D, 0.49 D and 0.30 D with the basis set STO-3G, 3-21G and…     

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

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

Ab Initio Study of Ring Flipping of the Overcrowded Peri-Substituted Naphthalenes

Ab inintio molecular orbital and density functional theory method were used to investigate the structural and dynamic behavior of 1,8-di-tert-butyl naphthalene (1), 1,8-bis(trimethylsilyl)naphthalene (2), 1,8-bis(trimethylgermyl)naphthalene (3), and 1,8-bis(trimethylstannyl)naphthalene (4). HF/3-21G//HF/3-21G results revealed that the ring flipping barrier height of compound 1–4 is 92.59, 32.13, 26.76, and 15.46 kJ mol^?1 respectively. The obtained results show that the transition state structure for ring flipping of the bulky-groups is in a planar form with naphthalene ring. Contrary to compound 1, the ring flipping of compounds 2–4 occurred easily at room temperature. Also, MP2/3-21G//HF/3-21G energy calculation, show that the enantiomerization energy of compounds 1–4 are 97.99, 33.24, 26.80, and 15.38 kJ·mol^?1 respectively. The required energy for ring inversion of compounds 1–4 are 85.09, 27.26, 21.54, and 10.21 kJ mol^?1 respectively, as calculated by B3LYP/3-21G//HF/3-21G method. It can be concluded that the lower energy barrier of the ring flipping of compounds 2–4 is related to the increasing of the bond lengths of Si—C, Ge—C, and Sn—C, in contrast to C—C bond.     

某些含氧酸中原子净电荷与pK_a的关系

在量子化学从头算 (分别采用HF/STO 3G ,HF/STO 3G ,HF/ 3 2 1G和MP2 /STO 3G)水平下 ,对 11种含氧酸分子 (H2 O ,H2 O2 ,HNO2 ,HNO3 ,HClO ,HClO2 ,HClO3 ,HCO2 H ,CH3 CO2 H ,HBrO和HIO)进行几何构型全优化和电子结构计算 ,将获得的电子结构数据与 pKa 实验值进行逐步回归分析 ,结果发现含氧酸中羟基氢原子净电荷、含氧酸中与羟基氧相连非氢原子净电荷两项与 pKa 呈良好的二元线性关系 ,表明含氧酸中羟基氧所连原子的净电荷在决定含氧酸的酸性上起重要作用