The intermolecular potential function of Smith–Thakkar type for C60

The intermolecular potential function of Smith–Thakkar type for C₆₀ has been proposed, and its expression is as follows

The unit of u(r) is J/mol, r is the distance between two C₆₀ molecules center and the unit is nm. Some properties of C₆₀ in the gas and crystal have been studied using the interaction potential of Smith–Thakkar type, such as stability of C₆₀ crystals, virial coefficient and lattice dynamics.     

Theoretical study of the effect of frit quality on chromatography using computational fluid dynamics

Summary Frits at both ends of a chromatographic column, especially for a preparative column, have significant influence on the flow distribution within the column and thus the column efficiency. However, frits have received little attention from chromatographers in the past. Here a theoretical study was conducted with the aid of CFD software FLUENT to investigate the effect of frits on the performance of homogeneous and heterogeneous chromatographic columns. A dimensionless number,FQ, was applied to characterize frit quality. This study visualized how frit quality affects the flow distribution and the concentration band, the shape of eluted pulse at the colum exit and column efficiency. Simulation results show that the development length of the flow distribution is related toFQ but has nothing to do with the packing heterogeneity. The curvature of the concentration band in a column depends onFQ and packing quality. This study shows column efficiency can be improved significantly by increasingFQ and/or frit permeability.     

高分子玻璃化转变过程的分子动力学模拟本科教学实验

玻璃化转变是高分子专业教学中的重要内容,转变过程中自由体积的变化和分子链段的运动是理解玻璃化转变的难点。在计算机上使用分子模拟的方法得到三维的可视化图形,可以更直观地观察到自由体积和分子链的变化。提炼近年分子模拟技术在玻璃化转变研究中的最新科研成果,得到既有理论基础又适合本科教学的课程素材。分子模拟可以实现了仪器测试无法达到的超快速升降温,并与仪器测试得到的实验结果相对照,验证了普通实验很难验证的理论观点,从而拓展了本科实验的范围,达到了很好的教学效果。     

熔体快速冷凝过程的微观结构演化

采用分子动力学模拟技术研究了金属间化合物AuCu3熔体的双体分布函数、键对、多面体、配位数等在快速凝固条件下随温度的变化情况，详细考察了AuCu3中微观组团随温度的演化特点.结果表明，AuCu3熔体降温至700 K时双体分布函数的第二峰已发生劈裂，液态金属中已经产生了非晶态；同时液体中的键对数及多面体数与温度的关系都表明，在上述向非晶转变的过程中，AuCu3熔体的确发生了微观结构组态的变化，其中以液体中的缺陷多面体随温度变化最为剧烈.     

Theoretical study and rate constant calculation for the reactions of SH (SD) with Cl2, Br2, and BrCl

The mechanisms of the SH (SD) radicals with Cl2 (R1), Br2 (R2), and BrCl (R3) are investigated theoretically, and the rate constants are calculated using a dual-level direct dynamics method. The optimized geometries and frequencies of the stationary points are calculated at the MP2/6-311G(d,p) and MPW1K/6-311G(d,p) levels. Higher-level energies are obtained at the approximate QCISD(T)/6-311++G(3df, 2pd) level using the MP2 geometries as well as by the multicoefficient correlation method based on QCISD (MC-QCISD) using the MPW1K geometries. Complexes with energies less than those of the reactants or products are located at the entrance or the exit channels of these reactions, which indicate that the reactions may proceed via an indirect mechanism. The enthalpies of formation for the species XSH/XSD (X = Cl and Br) are evaluated using hydrogenation working reactions method. By canonical variational transition-state theory (CVT), the rate constants of SH and SD radicals with Cl2, Br2, and BrCl are calculated over a wide temperature range of 200-2000 K at the a-QCISD(T)/6-311++G(3df, 2pd)//MP2/6-311G(d, p) level. Good agreement between the calculated and experimental rate constants is obtained in the measured temperature range. Our calculations show that for SH (SD) + BrCl reaction bromine abstraction (R3a or R3a') leading to the formation of BrSH (BrSD) + Cl in a barrierless process dominants the reaction with the branching ratios for channels 3a and 3a' of 99% at 298 K, which is quite different from the experimental result of k3a'/k3' = 54 +/- 10%. Negative activation energies are found at the higher level for the SH + Br2 and SH + BrCl (Br-abstraction) reactions; as a result, the rate constants show a slightly negative temperature dependence, which is consistent with the determination in the literature. The kinetic isotope effects for the three reactions are "inverse". The values of kH/kD are 0.88, 0.91, and 0.69 at room temperature, respectively, and they increase as the temperature increases.     

The stability of poly(m-carborane-siloxane) elastomers exposed to heat and gamma radiation

Poly(m-carboranyl-siloxane) elastomers containing a mixture of di-methyl- and methylphenyl-silyl units were synthesised using the ferric chloride catalysed condensation reaction between di-chloro-diorganosilane and bis(di-methylmethoxysilyl)-m-carborane. These elastomeric materials were originally developed to have greater stability to extreme thermal environments and retain tailorable physical and chemical properties relative to comparable non-carborane containing elastomers. Prepared samples were aged either by heating in air at elevated temperatures or by gamma irradiation from a ⁶⁰Co source. Multinuclear (¹H, ¹³C and ¹¹B) solid and solution state nuclear magnetic resonance (NMR) was used to assess degradation. This included measurements of segmental chain dynamics using a solid-echo pulse sequence reflecting changes in crosslink density and assessing changes to the carborane fragment by ¹¹B and ¹H Magic Angle Spinning (MAS) methods. Thermogravimetric measurements were also performed to assess thermal stability. Gamma radiation (to a dose of 1 MGy) was found to induce only a small degree of elastomer hardening as evidenced by a reduction in segmental chain dynamics. The carborane cage however, remained intact at these dose levels. Thermal degradation was observed to lead to oxidative crosslinking, the degree of which is dependent on temperature. At temperatures below 350 °C, only small changes in segmental dynamics were observed commensurate with only minor weight loss at this temperature. At temperatures above 350 °C, the degradation of the elastomer increased dramatically with decreased segmental dynamics and presumed partial oxidation of the carborane cage. The integrity of the m-carborane cage and the segmental dynamics were found to be significantly reduced at temperatures above 580 °C, in line with the known cage rearrangement temperature for icosahedral carboranes.     

Quantized Hamilton Dynamics

The paper describes the quantized Hamilton dynamics (QHD) approach that extends classical Hamiltonian dynamics and captures quantum effects, such as zero point energy, tunneling, decoherence, branching, and state-specific dynamics. The approximations are made by closures of the hierarchy of Heisenberg equations for quantum observables with the higher order observables decomposed into products of the lower order ones. The technique is applied to the vibrational energy exchange in a water molecule, the tunneling escape from a metastable state, the double-slit interference, the population transfer, dephasing and vibrational coherence transfer in a two-level system coupled to a phonon, and the scattering of a light particle off a surface phonon, where QHD is coupled to quantum mechanics in the Schrödinger representation. Generation of thermal ensembles in the extended space of QHD variables is discussed. QHD reduces to classical mechanics at the first order, closely resembles classical mechanics at the higher orders, and requires little computational effort, providing an efficient tool for treatment of the quantum effects in large systems.     

Multipeak distributions of first passage times in bistable dynamics in a model of a thermochemical system.

A master equation is used to study transitions between the stable limit cycle and stable focus in the two-variable bistable system. The distribution function of the mean first passage time between these attractors and the relative dispersion of the mean first return time from the stable focus to itself as a function of the intensity of fluctuations are calculated and discussed. A coherence resonance is observed for the return time from the focus to itself.     

Quantum chemical study of alternating N/B and N/C π bonds in (un)charged even‐membered 4n and 4n+2 cyclic systems and their related open systems

With semi‐empirical MO and ab initio calculations at different levels, we investigated the π conjugation of alternating X? Y bonds with X? Y for N/B and N/C combinations in an open and cyclic arrangement. Although the intrinsic symmetry is lost for the acyclic even‐membered compounds, the alternation is still reflected in its geometry and electron‐density transfer. For the cyclic π compounds, we focused our attention on borazine N₃B₃H₆ (D_3h symmetry), which is isoelectronic with benzene (D_6h symmetry). Specific attention is given to the electrophilic behavior of borazine with respect to CH and SiMe. The dynamics based on the results of FT‐ICR mass spectrometry was studied in more detail. In addition, the results of the cyclic systems with 4n and 4n+2 π electrons concerning their geometries are compared with the corresponding carbon compounds. Attention is also given to the dication of borazine, because of the corresponding triplet ground state of the benzene dication. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

超离子导体CaF2中的Ca2+亚晶格和F-亚晶格

CaF_2在熔化以前为超离子导体相。一些实验和理论的研究表明,在超离子导体相中,Ca~(2+)仍维持原有的面心立方骨架,而F~-则在Ca~(2+)骨架中运动。早期的分子动力学模拟结果表明Ca~(2+)的均方根位移仅约0.3A,而F~-的扩散系数可达2.6×10~(-5)cm·s~(-1),已是熔盐扩散系数量级。近年来的中子散射实验表明在扩散离子和近邻离子间存在着某种动力学相关。为解释这些事实,新近Gillan的分子动力学模拟表明扩散离子伴随着F~-亚晶格变形,而Kaneko和Ueda的分子动力学模拟则表明扩散离子伴随着近邻离子在同一运动方向的相关运动。进一步的研究尚待进行。八十年代初,Nelson等人提出了描述晶体、非晶和液态中键取向的键球谐函数方法。