固液相变的分子动力学研究Ⅰ─—(KI)_(108)离子簇的熔化和凝固(英)

通过分子动力学模拟，考查并分析了（KI）108离子簇的结构、能量和相变的动力学行为．在加热和冷却过程中，离子簇再现了熔化和凝固现象，而且熔化起始于立方体的其一顶点，熔化的离子簇不是球形的，说明了离子簇的非湿特征．根据结晶的成核速率，讨论了电子衍射实验中观察KI凝固的可能性     

Adsorption of water molecules in slit pores

We report molecular dynamics (MD) simulations on the adsorption of water in attractive and repulsive slit pores, where the slit and a bulk region are in contact with each other. Water structure, surface force and adsorption behavior are investigated as a function of the overall density in the bulk region. The gas–liquid transition in both types of pores occurs at similar densities of the bulk region.     

Hierarchical Structure‐Within‐Structure Morphologies in A2‐star‐(B‐alt‐C) Molecules

We employ dissipative particle dynamics (DPD) to examine the self‐assembly behavior of A₂‐star‐(B‐alt‐C) molecules. We successfully observe various types of hierarchical structure‐within‐structures, such as A‐formed spheres in the matrix formed by B and C alternating layers, hexagonally packed A‐formed cylinders in the matrix with B and C segregated layers, B and C alternating layers‐within‐lamellae, coaxial B and C alternating domains within hexagonally packed BC‐formed cylinders in the A‐matrix, and co‐centric BC‐alternating domains within BC‐formed spheres in the A‐matrix, by increasing the A composition. Generally speaking, the small length‐scale B and C segregated domains are in parallel to the large length‐scale structures. This hierarchical periodicity along the same axis as well as the various characteristic structures, that the A₂‐star‐(B‐alt‐C) copolymers display, are quite different from those in A‐block‐(B‐graft‐C) coil‐comb copolymers. Moreover, it is interesting to find that when the copolymer chain length increases, though the hierarchical structure type is maintained, the number of small length‐scale lamellae that can form within the large length‐scale structure increases. These hierarchical structures under various compositions are reported theoretically for the first time in the copolymer systems consisting of the alternating blocks, and are in good agreement with the most recent experimental work by Matsushita and co‐workers (Macromolecules 2007 , 40, 4023).     

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.     

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.