从头计算研究乙酰胆碱构象和分子静电势

在Hartree－Fock和二级Moller－Plesset微扰理论MP2水平，用6－31G对乙酰胆碱进行了从头计算研究，发现了旁式和反式两种稳定构象．经零点振动能校正后，在Hartree－Fock级别旁式较反式稳定0．24kJ·mol~(-1)；而在MP2级别，反式较旁式稳定2．07kJ·mol~(-1).谐振动频率在HF／6－31G水平计算.还给出了两种构象的MP2分子静电势图.     

分子间相互作用的Ab Initio计算研究——氢键复合物体系(HF)_2中的非键相互作用

本文以6-311 G~(**)为基组对氢键复合物(HF)_2进行了Ab Initio计算,给出了2个HF在距离及取向平衡点附近的势能面。用势诱导最小二乘拟合法计算了体系中各原子的原子电荷。用(exp-6-1)势函数分析法分解了2个HF间的氢键相互作用能,表明其中主要的相互作用是电荷转移而不是静电性的。2个HF分子间的势能曲线呈Morse函数性质。     

密度泛函理论研究镥二聚体低电子能态的相对稳定性、结构以及成键性质

用密度泛函理论方法研究了镥二聚体（Lu2）低能量电子态的性质,计算了电子态相对能量、平衡键长、振动频率以及基态解离能,考察了密度泛函性质、相对论有效势种类以及Hartree—Fock交换作用大小对计算结果的影响．结果表明,无论采用何种密度泛函和相对论有效势,体系的基态都为三重态,与其他一些基于分子轨道理论的从头计算方法得到的结论是一致的．另外,与分子轨道从头计算结果以及实验结果比较发现,采用杂化密度泛函理论和Stuttgart小核有效势计算得到的结果总体吻合最好．最后,特别分析研究了B3LYP计算中Hartree—Fock交换作用大小对基态键长和基态解离能的影响,发现随着交换作用的增大,键长增长,解离能减小,这是由于5d轨道杂化导致的共价成键作用减弱造成的．     

氮原子在Ni平坦和缺陷表面上的吸附和振动

 构造了氮－镍相互作用的5－参数Morse势，研究了氮原子在Ni（\r\n100），Ni（110）和Ni（111）平坦表面的吸附和振动，获得了氮原子\r\n在三个低指数表面的吸附位、吸附构型、结合能和本征振动等数据，计\r\n算结果与实验结果非常吻合．同时，与Ni（100）表面对比，系统研究\r\n了氮原子在Ni（510）台阶面的吸附和扩散．计算结果表明，氮原子在\r\n台阶下部形成最稳定的吸附态，台阶对下台面上扩散的氮原子形成捕获\r\n势，对上台面上扩散的氮原子形成反射势．     

石墨炉原子吸收法测定甲壳素中的砷

采用石墨炉原子吸收法直接测定甲壳素中的砷。以HF－HClO4溶解消化试样,用Mg(NO3)2作基体改进剂,灰化温度为1000℃,原子化温度为2300℃。该法相对标准偏差为1．02％,回收率在97．2％－105．8％之间,方法快速、简便,结果准确。     

Prediction of Crystal Structure and Density for Benzotrifuroxan(BTF) by AAPM

因苯并三氧化呋咱（ＢＴＦ）是具有广泛应用的很重要的化合物，且其晶体结构和性能的信息缺乏，故本文基于原子原子势方法（ＡＡＰＭ）预示了它的晶体结构和密度．分析了静电相互作用对ＢＴＦ晶体结构的影响．ＢＴＦ的非共面结构不能简单地不考虑电荷模型而加以预测．运用ａｂｉｎｉｔｉｏ计算求得的净原子电荷计算了它的结构参数和密度．结果发现，ＭＫＳ电荷能提供结构和密度预示与实测之间的最佳对应．由整体探索求得的分子堆积给出了合理成功的ＢＴＦ的结构和密度预测．密度预测的非精确度不超过０．０２ｇ·ｃｍ－３．     

用非线性统计理论方程计算气-液平衡

用非线性统计理论方法计算气－ 液平衡组成,参数有明确的物理意义,并可用纯组分的物性参数进行估算．用此方法计算了气－ 液平衡计算中最困难的二元含水体系的平衡组成．结果表明该方法简单,计算精度高,有理论和实际意义．     

热处理对Nylon—11结构的影响

对不同温度热处理过的Nylon－11样品，采用广角X射线衍射（WAXD），根据X射线衍射强度理论，用多重峰图解分峰方法研究了Nylon－11结构的变化，并导出了计算Nylon－11结晶度的公式．所得结果与密度法、DSC法进行了比较．同时也给出了Nylon－11的有关热力学参数，     

最弱受约束电子势模型理论下使用双广义拉盖尔多项式计算氦原子基态能量

计算了在最弱受约束电子势模型理论下使用双广义拉盖尔多项式的氦原子基态能量. 给出了氦原子基态能量期望值的表达式. 通过搜索它的极小, 找到了氦原子基态能量的极小值. 将我们的结果与双灼(zeta)函数的Hartree-Fock(HF)方法的结果进行比较后发现, 我们的结果略优于双灼(zeta)函数的Hartree-Fock(HF)方法的结果. 如此表明线性组合技术可以应用在最弱受约束电子势模型理论中.     

量子化学赝势法半从头计算程序MQMMP的设计

作者设计了赝势法半从头计算程序MQMMP,并计算了第二周期原子和N_2、HCI分子,其结果相当好地逼近用全电子从头计算法得到的结果。这表明赝势法半从头计算为大分子、尤其是含重原子的大分子的非经验量子化学计算提供了较为简便的方法。     

Quantum chemical simulation of silicon tetrachloride hydrogenation

The quantum chemical calculation of the activation parameters of the reduction of a SiCl₄ molecule with molecular hydrogen, atomic hydrogen, and atomic chlorine was performed. The energy parameters were determined within the scope of the density functional theory (DFT) with the complete geometry optimization by the unrestricted UB3LYP/6-311+G(d) method. The calculated activation energies allow one to exclude the participation of molecular hydrogen in processes of dehalogenation of chlorosilanes.     

TiC、TiC1-x、(Ti1-xNbx)C电子结构的计算

采用离散变分x_α法（DV-X_α法）对TiC理想晶体、空位和掺杂缺陷结构中的电子结构进行了计算．通过选取分子簇模型,模拟了理想晶体、空位和掺杂缺陷情况．采用多重散射离散变分X_α法,通过自治迭代来求解局域密度泛函方程,得到了各个分子簇模型的电子结构．分析计算结果发现,在理想TiC结构的态密度图中,费米能级位于两峰之间．但在费米能级处的电子态密度不为零,这提供了TiC导电性的来源．在空位模型中,发现电子态密度在费米能级处有较大的值,说明空位的存在有利于提高TiC的导电能力．对于Nb掺杂后的电子结构,在费米能级处存在一个电子态密度峰,因而也有利于提高其导电性．在计算的过程中考虑到了分子簇模型边界条件带来的电行转移效应对电子结构的影响,通过提供适当的环境势,得到了较精确的计算结果．与已有的计算结果进行了对比,有较好的一致性．     

New approach to the semiempirical calculation of atomic charges for polypeptides and large molecular systems

Starting from the bond polarization theory (BPT), a new semiempirical method for the calculation of net atomic charges is developed. The bond polarization theory establishes a linear dependence of atomic charges from the bond polarization energy. This energy is calculated from the hybrid orbitals forming a bond and the point charges within the neighborhood. Empirical parameters are introduced for the polarity of an unpolarized bond and for the change of the atomic charge with σ- and π-bond polarization. Because these parameters are linear, they can be calibrated directly using net atomic charges from ab initio calculations. This procedure was performed using the charges from STO3G calculations on a set of 18 amino acids. Using the two parameters for CH, OH, σ-CO, and NH bonds and the three parameters for CC, CO, and CN bonds, the 350 ab initio charges can be reproduced with high accuracy by solving sets of linear equations for the charges. The calculation of charges for large molecular systems including all inter- and intramolecular mutual polarizations requires only a few seconds (up to 100 atoms) or minutes (700 atoms) on a PC. This procedure is well suited for the application in molecular mechanics or molecular dynamics programs to overcome the limitations of most force fields used up to now. One of the weakest points in these programs is the use of fixed or topological charges to define the electrostatic potential. As an application of the new method, we calculated the interaction energy of an ion with valinomycin. This ring molecule forms octahedral oxygen cages around ions like potassium and acts thereby as selective ion carrier. To accomplish this function, valinomycin has to strip off the hydratization spheres of the ions, and therefore its preference for certain types of ions could be deduced from the interaction energies. © 1994 by John Wiley & Sons, Inc.     

Effect of the twisted alignment on the liquid crystal wave‐front corrector

The effect of twisted alignment on the phase modulation of a liquid crystal wave‐front corrector was investigated. First, the effect of twisted alignment is discussed in terms of the modulation principle of the liquid crystal molecule. Only partial incident light can be modulated because of the effect of the twisted alignment. Other unmodulated light will affect the correction accuracy and the resolution of the image. The blazed grating method is proposed to solve this problem. Adaptive correction was performed without the blazed grating method and the correction results are poor. A similar adaptive correction experiment was performed with the blazed grating method and a better correction result is obtained. The residual averaged wave‐front errors are PV = 0.101λ and RMS = 0.015λ and a resolvable image is obtained.     

Density-difference maps in quantum chemistry

Electron density-difference maps, used to study the changes that occur when a molecule changes its state or when the nuclei of a molecule change their relative positions, are generally useful only if the atomic densities cancel when one molecular density distribution is subtracted from the other. When, as in the case of the nonrigid internal rotation in ethane, such a cancellation of atomic densities is not possible the method of simple subtraction is no longer appropriate. It is shown that useful density-difference maps can nevertheless be obtained, when the changes in geometric parameters are small, by the calculation of two generalized density-difference functions: a point difference function which allows a comparison of the densities at corresponding points in the two systems, and a volume difference function to compare the amounts of charge in corresponding regions. The method is illustrated by consideration of a change in bond length of the nitrogen molecule and by the nonrigid internal rotation in ethane.Presented at the XIVth Quantum Theory Conference in Canterbury, England, September 1981.     

Building molecular charge distributions from fragments: Application to HIV-1 protease inhibitors

Interaction energies are a function of the molecular charge distribution. In previous work, we found that the set of atomic partial charges giving the best agreement with experimental vacuum dipole moments were from density functional theory calculations using an extended basis set. Extension of such computations to larger molecules requires an atomic partial charge calculation beyond present computational resources. A solution to this problem is the calculation of atomic partial charges for segments of the molecule and reassociation of such fragments to yield partial charges for the entire molecule. Various partitions and reassociation methods for five molecules relevant to HIV-1 protease inhibitors are examined. A useful method of reassociation is introduced in which atomic partial charges for a large molecule are computed by fitting to the combined electrostatic potential calculated from the fragment partial charges. As expected, the best sites for partitions are shown to be carbon—carbon rather than carbon—nitrogen bonds. © 1997 by John Wiley & Sons, Inc.     

Prediction of molecular solvation free energy based on the optimization of atomic solvation parameters with genetic algorithm

We propose an improved solvent contact model to estimate the solvation free energy of an organic molecule from individual atomic contributions. The modification of the solvation model involves the optimization of three kinds of parameters in the solvation free energy function: atomic fragmental volume, maximum atomic occupancy, and atomic solvation parameters. All of these atomic parameters for 24 atom types are developed by the operation of a standard genetic algorithm in such a way as to minimize the difference between experimental and calculated solvation free energies. The data set for experimental solvation free energies is divided into a training set of 131 compounds and a test set of 24 compounds. Linear regressions with the optimized atomic parameters yield fits with the squared correlation coefficients (r2) of 0.89 and 0.86 for the training set and for the test set, respectively. Overall, the results indicate that the improved solvent contact model with the newly developed atomic parameters would be a useful tool for rapid calculation of molecular solvation free energies in aqueous solution.     

A study of amino-protecting groups using the polarizable continuum model (PCM)

In order to better understand the performance of 1,2-dimethyl-5-acetyl barbituric acid (DMB) as an amino protecting group relative to 5,5-dimethylcyclohexane-1,3-dione (DMD), ab initio calculations were performed. pK_a calculations using the PCM model indicated that both molecules are more acidic in the enol form. Therefore, the protecting reaction of these molecules should involve the anions formed from the loss of a proton from the enol compounds. Contrary to what would be expected, the larger efficiency exhibited by the DMB molecule cannot be attributed to an extension of the electronic conjugation effect. In the absence of any other noticeable effect that could be responsible for the greater efficiency of the DMB molecule, we are inclined to believe that the difference could be accounted for by the presence of two independent centers of conjugation.This paper is dedicated to Jacopo Tomasi in recognition of his outstanding contribution to the field of computational chemistry in solution. The authors are honored to contribute to this volume; especially so for two of them (COS and MACN) who have the privilege of his friendship.Acknowledgements The authors would like to thank the Brazilian research agencies CNPq, CAPES and FAPERJ for the financial support. C. O. da Silva thanks the Dipartimento di Chimica e Chimica Industriale, University of Pisa, where the MCSCF calculations were performed.     

Outer shell ionization potentials for ethane by a many-body Green's function method

A calculation based upon the many-body Green's function method is employed to obtain the outer shell vertical ionization potentials of the ethane molecule. An extended basis set is employed to represent the approximate optical potential, derived by the functional derivative approach, as well as the one-electron Green's function. The results obtained confirm a²E_g state for the ion arising from the first ionization process. Supported by a fellowship of the Scuola Normale Superiore.