高岭石吸附乙烯和苯的Delft分子力学研究

运用Ｄｅｌｆｔ分子力学（ＤＭＭ）程序及其粘土和共轭烯烃力场，计算研究了高岭石对乙烯和苯的吸附作用，探讨了吸附对粘土晶体和有机分子的结构、电荷分布和能量的影响，求得了高岭石吸附乙烯和苯的吸附热等重要物理量．     

适用于TATB,RDX,HMX含能材料的全原子力场的建立与验证

报道一个适用于三种常见的含能材料分子三硝基三氨基苯(TATB),环三亚甲基三硝胺(RDX),环四亚甲基四硝胺(HMX)的全原子力场.力场采用广泛使用的力场函数形式,其中键参数通过拟合量子化学密度泛函计算的数据获得,电荷参数和范德华参数通过拟合相应的分子晶体的物理性质(密度和升华焓)优化得到.通过计算分子和分子晶体的性质显示该力场可以用来准确地预测分子结构、分子振动频率和分子晶体的晶胞参数、密度和升华焓.进一步的验证显示该力场可用来较为准确地预测分子晶体的状态方程和机械模量.     

硝胺及其甲基衍生物的从头计算研究 3: 甲硝胺及其同位素衍生物的谐性力场和振动光谱

用TEXAS从头计算程序,取STO-4-21G基组,计算了甲硝胺的谐性力场和振动光谱.直接理论计算的谐性力场经由其他分子转移来的经验校正因子校正后,提供了甲硝胺振动基频的预测.预测值和甲硝胺分子在气相中的振动光谱实验值之间的平均偏差为31cm^-1.为了获得更合适的气相甲硝胺振动力场和预测它的同位素衍生物的振动光谱,我们优化了一组新的校正因子,使理论值和实验值的平均偏差减为8.9cm^-1.用这组校正因子得到的力场预测了三个同位素衍生物的振动光谱,其同位素位移的理论预测值和实验值符合良好.     

硝胺及其甲基衍生物的从头计算研究 Ⅲ.甲硝胺及其同位素衍生物的谐性力场和振动光谱

用TEXAS从头计算程序,取STO-4-21G基组,计算了甲硝胺的谐性力场和振动光谱.直接理论计算的谐性力场经由其他分子转移来的经验校正因子校正后,提供了甲硝胺振动基频的预测.预测值和甲硝胺分子在气相中的振动光谱实验值之间的平均偏差为31cm~(-1).为了获得更合适的气相甲硝胺振动力场和预测它的同位素衍生物的振动光谱,我们优化了一组新的校正因子,使理论值和实验值的平均偏差减为8.9cm~(-1).用这组校正因子得到的力场预测了三个同位素衍生物的振动光谱,其同位素位移的理论预测值和实验值符合良好.     

硝胺及其甲基衍生物的从头计算研究 3: 甲硝胺及其同位素衍生物的谐性力场和振动光谱

用TEXAS从头计算程序,取STO-4-21G基组,计算了甲硝胺的谐性力场和振动光谱.直接理论计算的谐性力场经由其他分子转移来的经验校正因子校正后,提供了甲硝胺振动基频的预测.预测值和甲硝胺分子在气相中的振动光谱实验值之间的平均偏差为31cm^-1.为了获得更合适的气相甲硝胺振动力场和预测它的同位素衍生物的振动光谱,我们优化了一组新的校正因子,使理论值和实验值的平均偏差减为8.9cm^-1.用这组校正因子得到的力场预测了三个同位素衍生物的振动光谱,其同位素位移的理论预测值和实验值符合良好.     

基于能量转移机理的水溶性荧光共轭聚合物体系的设计制备及应用

水溶性共轭聚合物由于具有优异的光学性能,如强大的光捕获能力和独特的分子线效应,在化学、生物和医疗领域都有良好的应用前景。共轭聚合物受光激发后,所产生的激子可以沿着长程共轭的π骨架自由迁移,在分子骨架的任何位点都可以转移给能量受体,从而实现高效的能量传递。因此,水溶性共轭聚合物适合于构建各种能量传递体系,从而实现荧光信号放大、多色发光、活性氧产生效率提高等光学性能的优化。本文简述了基于水溶性共轭聚合物构建能量转移体系的原理和方法,总结了其在生物传感、生物成像、光动力杀菌和光动力抗癌等领域的应用,最后对水溶性共轭聚合物存在的主要问题以及未来的发展方向进行了分析和展望。     

粘土的分子力学力场优化及结构、性能研究

以DELPHI分子力学程序优化了粘土力场。计算研究了高岭石、地开石、珍珠陶土、叶蜡石1Tc和叶蜡石2M等五种粘土的结构、电荷分布和红外光谱等性质，所得结果与实验结果比较相符。     

XSO_2NCO(X=F,Cl)分子振动光谱的从头算研究——谐性力场和频率分布

采用从头算方法以 6 31G 基组对XSO2 NCO(X =F ,Cl)分子的振动谐性力场和红外光谱进行了研究 .由最小二乘法通过拟合FSO2 NCO分子的理论和实验频率优化了一套力常数校正因子 ,然后将其迁移到ClSO2 NCO分子中 ,以产生一个纯理论预测的振动基频 .利用校正后的力场计算得到FSO2 NCO和ClSO2 NCO的振动频率 ,与实验值比较平均偏差分别为 3cm-1和 5cm-1.根据用校正后力场进行简正坐标分析得到的势能分布 (PED)和从头算红外光谱强度值对此二分子的振动基频的指认进行了讨论 .     

TATB基PBX结合能和力学性能的理论研究

以SCF-MO-AM1方法和MM-COMPASS力场, 对TATB (1,3,5-三氨基-2,4,6-三硝基苯)与系列高聚物组成的PBX(高聚物粘结炸药)尺寸匹配原子簇, 分别进行全优化几何构型计算, 发现两种方法求得的结合能存在良好的线性关系. 对TATB (3×3×4)超晶胞及其与系列氟聚物组成的双组分PBX, 实施COMPASS力场下的分子动力学(MD)周期性模拟计算, 首次求得其弹性系数、模量和泊松比, 发现添加少量高聚物即能有效改善炸药的力学性能.     

共轭烯烃的分子几何构型、电子结构和生成热的DMM和AM1、PM3的比较研究

运用Ｄｅｌｆｔ分子力学（ＤＭＭ）力场和程序以及半经验分子轨道ＡＭ１和ＰＭ３方法计算研究了丁二烯、苯、甲苯、联苯、苯乙烯、富烯、、环辛四烯、［２，２］对环烷和菲等１０个共轭烯烃分子的几何构型、电子结构和生成热．ＤＭＭ计算的几何构型和生成热与实验结果相吻合，电荷分布结果与从头计算结果较接近．ＡＭ１和ＰＭ３计算的几何构型较好，但计算的生成热与实验结果偏差较大．ＰＭ３计算值比ＡＭ１的稍好．     

Structures and stabilities of diacetylene-expanded polyhedranes by quantum mechanics and molecular mechanics

The structures, heats of formation, and strain energies of diacetylene (buta-1,3-diynediyl) expanded molecules have been computed with ab initio and molecular mechanics calculations. Expanded cubane, prismane, tetrahedrane, and expanded monocyclics and bicyclics were optimized at the HF/6-31G(d) and B3LYP/6-31G(d) levels. The heats of formation of these systems were obtained from isodesmic equations at the HF/6-31G(d) level. Heats of formation were also calculated from Benson group equivalents. The strain energies of these expanded molecules were estimated by several independent methods. An adapted MM3 molecular mechanics force field, specifically parametrized to treat conjugated acetylene units, was employed for one measure of strain energy and as an additional method for structural analysis. Expanded dodecahedrane and icosahedrane were calculated by this method. Expanded molecules were considered structurally in the context of their potential material applications.     

Molecular mechanics calculations of conjugated amides and an ab initio investigation of acrylamide and its β-amino derivative: Conformational analysis and rotational barriers

Conformations and rotational barriers in a series of conjugated primary and tertiary amides have been analyzed by a modified MM2(91) force field, which treats the amide nitrogen as part of the conjugated system by redefining the atom type for the nitrogen. Ab initio molecular orbital calculations at the MP2/6-31G* level have been performed on the stable conformers and transition structures of acrylamide and β-trans-aminoacrylamide. The results have been used, with published experimental and computational data, to generate parameters for the MM2 force field. The force field has been applied to various conjugated amides, such as reduced nicotinamide adenine dinucleotide (NADH) and NAD⁺ analogues, nicotinamide, urea, vinylogous urea derivatives, and nucleic acid bases. The fundamental difference between primary and tertiary conjugated amides with respect to both conformation and barrier is highlighted. © 1996 by John Wiley & Sons, Inc.     

MM3 force field prediction of the enantioselective preference in the asymmetric synthesis of a chiral 2-cyclohexen-1-ol using a chiral lithium amide reagent

Enantioselective preference in the asymmetric synthesis where cyclohexene oxide is transformed enantioselectively to chiral (S)- or (R)-2-cyclohexen-1-ol by the reaction with the appropriate chiral lithium amide reagent has been evaluated theoretically using the MM3 force field. The plausible possible structures for each precursor (reaction intermediate complex) leading to a (S)- or (R)-2-cyclohexen-1-ol have been optimized with the extended MM3 force field applicable to the lithium amide functional group, and the populations of their (S)- or (R)-reaction intermediate complexes at an ambient temperature (298 K) were calculated. The initial structure for evaluating the reaction intermediates of this asymmetric synthesis was constructed on the basis of the optimized ab initio transition state structure (MP2/6-31+G^∗) comprising lithium amide LiNH₂ and propene oxide. To the thus obtained transition state structure composed of LiNH₂ and propene oxide, the other remaining Cartesian coordinates for the actual reaction intermediates composed of the chiral lithium amides and cyclohexene oxide were added to make the reaction intermediate structure. The conformational search for the reaction intermediate has been carried out by using the Stochastic search Algorithm, and the optimized geometries and their conformational energies (steric energies) have been calculated by the MM3 force field. The populations calculated from the conformational energies of the reaction intermediate leading to the (S)- or (R)-2-cyclohexen-1-ol were shown to be linearly well correlated with the experimentally reported enantiomer excess (% ee) values. The critical factors to control the enantioselectivity were investigated on the basis of the optimized structures of the reaction intermediate complexes. The MM3 force field approach was shown to be applicable to the theoretical evaluation of the enantioselectivity and be useful for designing a new functional chiral lithium amide reagent for the asymmetric synthesis.     

分子力场进展

分子力学（简称ＭＭ）是近年来化学家常用的一种计算方法。与量子力学从头计算和半经验方法相比，用分子力学处理大分子可以大大节省计算时间，而且，在大多数情况下，用分子力学方法计算得到的分子几何构型参数与实验值之间的差值可在实验误差范围之内。所以，分子力学是研究生物化学体系的有效和可行的手段。分子力学的核心是分子力场。本文介绍了分子力场的量子力学背景、分子力场和光谱力场之间的关系。分子力场的一般形式、分力     

Excited States of Conjugated Hydrocarbons Using the Molecular Mechanics - Valence Bond (MMVB) Method: Conical Intersections and Dynamics

We developed the molecular mechanics—valence bond (MMVB) method following an original suggestion of Jean-Paul Malrieu and coworkers. By coupling a parameterized Heisenberg Hamiltonian to a standard classical force field (MM2), reliable ground and excited state geometries of conjugated hydrocarbons can be rapidly optimized. The MMVB method was central to our development of algorithms for locating conical intersections and calculating their associated decay dynamics. Here, we briefly review the chemical applications of MMVB to date, and present two new studies using the photostability of pyracylene and the excited state decay dynamics of the photochromic dihydroazulene/vinylheptafulvene (DHA/VHF) reaction.     

Statistical mechanically averaged molecular properties of liquid water calculated using the combined coupled cluster/molecular dynamics method

Liquid water is investigated theoretically using combined molecular dynamics (MD) simulations and accurate electronic structure methods. The statistical mechanically averaged molecular properties of liquid water are calculated using the combined coupled cluster/molecular mechanics (CC/MM) method for a large number of configurations generated from MD simulations. The method includes electron correlation effects at the coupled cluster singles and doubles level and the use of a large correlation consistent basis set. A polarizable force field has been used for the molecular dynamics part in both the CC/MM method and in the MD simulation. We describe how the methodology can be optimized with respect to computational costs while maintaining the quality of the results. Using the optimized method we study the energetic properties including the heat of vaporization and electronic excitation energies as well as electric dipole and quadrupole moments, the frequency dependent electric (dipole) polarizability, and electric-field-induced second harmonic generation first and second hyperpolarizabilities. Comparisons with experiments are performed where reliable data are available. Furthermore, we discuss the important issue on how to compare the calculated microscopic nonlocal properties to the experimental macroscopic measurements.     

Force field calculations (MM3) on glyoxal,quinones, and related compounds

A general force field type of calculation has been devised in connection with MM3 to treat 1,2- and 1,4-diketones, both when they are not conjugated (as in derivatives of glyoxal) and when they are conjugated (as in derivatives of ortho- and para-benzoquinone). The molecular structures, moments of inertia, dipole moments, and vibrational spectra have been examined for about 15 compounds, some in several conformations. Ab initio calculations (6-31G*) have been used to determine quantities that have not been previously defined by experiment. In general, the force field permits the calculation of the structures with high accuracy, and the spectroscopic and conformational energy data with fair accuracy. © 1994 by John Wiley & Sons, Inc.     

Optimization of a molecular mechanics force field for type-II polyoxometalates focussing on electrostatic interactions: a case study

In this study, we have focussed on type-II polyanions such as [M(7)O(24)](6-), and we have developed and validated optimized force fields that include electrostatic and van der Waals interactions. These contributions to the total steric energy are described by the nonbonded term, which encompasses all interactions between atoms that are not transmitted through the bonds. A first validation of a stochastic technique based on genetic algorithms was previously made for the optimization of force fields dedicated to type-I polyoxometalates. To describe the new nonbonded term added in the functional, a fixed-charged model was chosen. Therefore, one of the main issues was to analyze that which partial atomic charges could be reliably used to describe these interactions in such inorganic compounds. Based on several computational strategies, molecular mechanics (MM) force field parameters were optimized using different types of atomic charges. Moreover, the influence of the electrostatic and van der Waals buffering constants and 1,4-interactions scaling factors used in the force field were also tested, either being optimized as well or fixed with respect to the values of CHARMM force field. Results show that some atomic charges are not well adapted to CHARMM parameters and lead to unrealistic MM-optimized structures or a MM divergence. As a result, a new scaling factor has been optimized for Quantum Theory of Atoms in Molecules charges and charges derived from the electrostatic potential such as ChelpG. The force fields optimized can be mixed with the CHARMM force field, without changing it, to study for the first time hepta-anions interacting with organic molecules.