C_（60）分子间相互作用势及其晶体压缩系数的计算

根据Ｃ（６０）晶体升华焓的实验值求得了Ｃ（６０）分子间相互作用势的解析表达式。用该表达式计算了Ｃ（６０）分别为面心立方和六方晶格时晶体的结合能，讨论了不同晶格的稳定性问题，解释了Ｃ（６０）晶体一般是面心立方的原因。本文还根据Ｃ（６０）分子间相互作用势计算了Ｃ（６０）晶体的压缩系数，并与文献的实验值进行了比较，所得结果与实验值相符合。     

C60固体中的相变

本文用（ｅｘｐ－６－１）势对固体Ｃ６０中分子间的相互作用进行计算，给出常温以下固体Ｃ６０中分子存在的状态，计算表明，常温下Ｃ６０高速旋转，随着温度降低，旋转速度逐渐减小，并且旋转轴向也逐渐趋于单一，当降到－２６０Ｋ时发生一级相变，此时ＦＣＣ晶胞中（０，０，０），１／２，１／２，０），（１，１／２，１／２）及（１／２，０，１／２）四套格点上的分子各只绕［１１１］三重轴旋转，Ｃ６０晶体从ＦＣＣ结构变为     

A2BC60和A3C60晶体中的短程相互作用

研究了Ａ３Ｃ６０和Ａ２ＢＣ６０晶体的结构和稳定性。计算了它们的Ｍａｄｅｌｕｎｇ常数、结合能、晶格常数和体积弹性模量,并讨论了晶体中的短程相互作用。结果显示,短程相互作用结构和八面体空隙碱金属的稳定性具有极大的影响。ＣｓＫ２Ｃ６０、ＲｂＫ２Ｃ６０和Ｋ３Ｃ６０的晶格常数计算结果与实验测量值非常符合。     

碳六十内嵌复合物中钠及氯与碳六十之间的相互作用

本文用量化学及分子力学方法计算了（Ｎａ＠Ｃ６０）及（Ｃｌ＠Ｃ６０）中Ｎａ及Ｃｌ与Ｃ６０之间的相互作用，总的相互作用是吸引。在静电相互作用、色散相互作用及电子云重排斥相互作用三项中前一项占９０％以上，有着决定性的影响，色散相互作用约占５％。对电子云重叠排斥作用来说，由于Ｃｌ＾－的半径比Ｎａ＾＋的大很多，所以在（Ｃｌ＠Ｃ６０）中这种作用要比在（Ｎａ＠Ｃ６０）中大很多。     

有机金属化合物中二级化学键的理论研究（V）——二级化学键生成的分子力学

用分子力学方法对一些金属汞化合物进行了探讨和研究，通过与对考体系对比以及内旋转势垒计算的方式，说明了在诸如cis-ClHgCHCHCl(I),o-C6H4HgCL2(IV),Cl-Hg-CH2-CH2-CN(VIIa),Cl-Hg-CH2-CH2-C6H5(XIa)分子中，非相邻的Cl,Hg原子间存在着弱相互作用。     

β-HMX晶体结构及其性质的高水平计算研究

对β-HMX晶体进行DFT-B3LYP/6-31G计算研究，求得其能带结构和电子结构，探讨了结构-性能关系。从带隙约为5.7eV推知HMX晶体为绝缘体，计算所得晶格间相互作用能为149.29kJ.mol^-^1；计算的升华热(154.08kJ.mol^-^1)与实验值相符。晶体中相邻分子的硝基氧之间点电荷静电势要高，表明该处晶格易于断裂、并可能成为冲击引爆中心。从原子间距和原子间重叠布居分析，发现β-HMX分子之间存在微弱的氢键；环外以硝胺键较强，环内以N-C键较弱，预示该二部位易于引发分解和起爆。     

无定形1,2-聚丁二烯的分子链结构与分子链的相互作用

本工作用聚合物玻璃化转变温度,内聚能密度和分子链内旋转势垒及异构化能表征分子链的相互作用;提出了在混合溶剂中用特性粘数法测定内聚能密度,和从聚合物玻璃化转变时的自由体积及稀溶液性质计算分子链内旋转参数等方法。     

金属富勒烯的价π轨道成键,杜瓦结构和无π轨道自由基结构

对富勒烯的价π轨道成键和无π轨道自由基进行了理论研究，得到内接和外接一种和多种金属富勒烯的各种可能对称性和结构，除了Ｃ６０（Ｉｈ），Ｃ３４（Ｃ３ｖ）和Ｃ２８（Ｔｄ）外，我们还研究了Ｃ３３（Ｃ３），Ｃ５６（Ｔｄ）和Ｃ６８（Ｔｄ）。这样的价键观点将为分子轨道计算以及分子构件配置的相反观点提供新的见解。     

ClnAlNHn和HnAlNHn(n=1～3)的振动光谱与化学键性质的研究

用从头计算分子轨道法和密度泛函理论，在ＨＦ／６－３１Ｇ＊和Ｂ３ＬＹＰ／６－３１Ｇ＊水平上对ＣｌｎＡｌＮＨｎ和ＨａＡｌＮＨｎ（ｎ＝１～３）及其碎片分子的几何构型、电子结构、振动光谱和化学热力学性质进行了理论研究。结果表明，优化几何参数与实验值相。ＣｌＡｌＮＨ２和Ｈ２ＡｌＮＨ２分子中，Ａｌ－Ｎ键为由一个σ键和一个π键组成的双重键，旋转势垒分别为３４．１０和５４．３５ＫＪ．ｍｏｌ＾－１，而Ｃｌ３ＡｌＮＨ     

分子内旋转势垒的计算

乙烷类型分子中C—C单键的旋转,既不是完全自由的,也不是完全受阻的。这是K.S.Pitzer等人研究了乙烷分子的熵和热焓后得出的结论。他所建立的乙烷分子配分函数,只有考虑到在分子内存在有旋转势垒,才能使这些量的计算结果与实验值相符合。由此他得到的乙烷旋转势垒大约为3Kcal/mol。     

Concomitant Polymorphism in an Organic Solid: Molecular and Crystal Structure and Intra- and Intermolecular Potential Contributions to tert-Butyl and Methyl Group Rotation

We investigate the relationship between structure (crystal and molecular) and tert-butyl and methyl group dynamics in 2-(tert-butyl)-9-(4-(tert-butyl)phenyl)anthracene. Powder and single-crystal X-ray diffraction, taken together, show that different polycrystalline samples recrystallized from different solvents have different amounts of at least four polymorphs (crystallites having different crystal structures), of which we have identified three by single crystal X-ray diffraction. The molecules in the asymmetric units of the different crystal structures differ by the dihedral angle the tert-butylphenyl group makes with the anthracene moiety. Ab initio electronic structure calculations on the isolated molecule show that very little intramolecular energy is required to change this angle over a range of about 60° which is probably the origin of the concomitant polymorphism (crystals of more than one polymorph in a polycrystalline sample). Solid state ¹H nuclear magnetic resonance (NMR) spin-lattice relaxation experiments support the powder and single-crystal X-ray results and provide average NMR activation energies (closely related to rotational barriers) for the rotation of the tert-butyl groups and their constituent methyl groups. These barriers have both an intramolecular and an intermolecular component. The latter is sensitive to the crystal structure. The intramolecular components of the rotational barriers of the two tert-butyl groups in the isolated molecule are investigated with ab initio electronic structure calculations.     

CF3 rotation in 3-(trifluoromethyl)phenanthrene. X-ray diffraction and ab initio electronic structure calculations

The molecular and crystal structure of 3-(trifluoromethyl)phenanthrene has been determined by X-ray diffraction. The structure of the isolated molecule has been calculated using electronic structure methods at the HF/3-21G, HF/6-31G, MP2/6-31G and B3LYP/6-31G levels. The potential energy surfaces for the rotation of the CF3 group in both the isolated molecule and cluster models for the crystal were computed using electronic structure methods. The barrier height for CF3 rotation in the isolated molecule was calculated to be 0.40 kcal mol(-1) at B3LYP/6-311+G//B3LYP/6-311+G. The B3LYP/6-31G calculated CF3 rotational barrier in a 13-molecule cluster based on the X-ray data was found to be 2.6 kcal mol(-1). The latter is in excellent agreement with experimental results from the NMR relaxation experiments reported in the companion paper (Beckmann, P. A.; Rosenberg, J.; Nordstrom, K.; Mallory, C. W.; Mallory, F. B. J. Phys. Chem. A 2006, 110, 3947). The computational results on the models for the solid state suggest that the intermolecular interaction between nearest neighbor pairs of CF3 groups in the crystal accounts for roughly 75% of the barrier to rotation in the solid state. This pair is found to undergo cooperative reorientation. We attribute the CF3 reorientational disorder in the crystal as observed by X-ray diffraction to the presence of a pair of minima on the potential energy surface and the effects of librational motion.     

Modeling fluid diffusion using the lattice density functional theory approach: counterdiffusion in an external field

The dependence of the diffusivity on temperature, pressure, and composition is not understood well; consequently, data is preferred significantly over correlations in most practical situations. Even in dilute gases, the contributions of attractions and repulsions to the diffusivity are difficult to understand on a molecular level without performing simulations. We have developed a Lattice Density Functional Theory (LDFT) approach for modeling diffusion to supplement existing methods that are very rigorous but computationally demanding. The LDFT approach is analogous to the van der Waals equation in how it accounts for molecular interactions in that it has first-order corrections to ideal behavior; it is an extension of the Equilibrium LDFT for adsorption and phase behavior. In this work, the LDFT approach is presented and demonstrated by modeling the problem of color counterdiffusion in an externally-applied potential field. This potential field, in combination with the intermolecular potential function, creates a diffusion regime in which repulsions cause oscillations in the density profile. Using the LDFT approach, the oscillations were described and attributed to nearest-neighbor and next nearest-neighbor interactions. The LDFT approach gives qualitative and quantitative agreement with dual control-volume Grand Canonical Molecular Dynamics simulations.     

A complete nearest-neighbor force field model for C60

A force field model is developed for C(60) that features 13 force constants representing all interactions between nearest-neighboring atoms. The model is compared with, and tested against, other force field models in the literature. Force constants for C(60) are then deduced by fitting the model to the 14 known optically accessible vibrational frequencies of the molecule. Finally, the model is fitted to two existing theoretical calculations of the complete vibrational spectrum of C(60). Fair agreement is obtained with the theoretical calculations, implying that interactions with atoms other than nearest neighbors are small.     

Pnictogen Interactions with Nitrogen Acceptors

Stabilizing nitrogen pnictogen bond interactions were measured using molecular rotors. Intramolecular C=O⋅⋅⋅N interactions were formed in the bond rotation transition states which lowered the rotational barriers and increased the rates of rotation, as measured by EXSY NMR. The pnictogen interaction energies show a very strong correlation with the positive electrostatic potential on nitrogen, which was consistent with a strong electrostatic component. In contrast, the NBO perturbation and pyramidalization analyses show no correlation, suggesting that the orbital-orbital component is minor. The strongest C=O⋅⋅⋅N pnictogen interactions were comparable to C=O⋅⋅⋅C=O interactions and were stronger than C=O⋅⋅⋅Ph interactions, when measured using the same N-phenylimide rotor system. The ability of the nitrogen pnictogen interactions to stabilize transition states and enhance kinetic processes demonstrates their potential in catalysis and reaction design.     

立方面心C60晶体晶格振动的研究

采用C60分子之间相互作用势的Kihra形式,研究了立方面心C60晶体的晶格振动问题,得到了质心振动沿[111]、[110]及[100]方向的声子散射圆频率分布曲线及C60晶格振动频率的态密度分布.采用所得到的C60晶格振动频率的态密度分布,计算了晶体C60在298 K时的等压热容,所得数值与实验值相符.     

The quenching of isopropyl group rotation in van der Waals molecular solids

X-ray diffraction experiments are employed to determine the molecular and crystal structure of 3-isopropylchrysene. Based on this structure, electronic structure calculations are employed to calculate methyl group and isopropyl group rotational barriers in a central molecule of a ten-molecule cluster. The two slightly inequivalent methyl group barriers are found to be 12 and 15 kJ mol(-1) and the isopropyl group barrier is found to be about 240 kJ mol(-1), meaning that isopropyl group rotation is completely quenched in the solid state. For comparison, electronic structure calculations are also performed in the isolated molecule, determining both the structure and the rotational barriers, which are determined to be 15 kJ mol(-1) for both the isopropyl group and the two equivalent methyl groups. These calculations are compared with, and are consistent with, previously published NMR (1)H spin-lattice relaxation experiments where it was found that the barrier for methyl group rotation was 11+/-1 kJ mol(-1) and that the barrier for isopropyl group rotation was infinite on the solid state NMR time scale.     

Effect of packing and tilt on the rotational barriers of an amino, nitro-substituted phenylene ethynylene trimer

Rotational potentials are computed for heptamers and nonamers of an amino, nitro-substituted phenylene ethynylene trimer molecule. A herringbone and a parallel-slipped packing arrangement are considered. The effect of tilting the molecules with respect to the surface as well as the effect of the gold support are also taken into account. The herringbone structure with the molecules perpendicular to the surface has a low rotational barrier (2 kcal/mol). Tilting the molecules by 30 degrees increases the rotational barriers significantly (16 kcal/mol). The parallel-slipped structure has rotational barriers of 7 kcal/mol. Including the effect of the gold support increases the rotational barriers for tilted molecules but has little effect on perpendicular molecules.