Molecular simulation study of the glass transition in a soft primitive model for ionic liquids

In this paper, we present a molecular dynamics study of the glass transition for a soft-core primitive model for ionic liquids, in which cations are fully flexible chains of tangent soft spherical monomers, being the positively charged monomer at one of the ends of the chain, and anions as charged soft spheres. We have monitored transport coefficients such as the self-diffusion coefficients and the shear viscosity, as well as correlation functions such as the mean-square displacement, the self-intermediate scattering function, and probes of heterogeneous dynamics such as the van Hove distribution function and the four-points susceptibility. The analysis of these properties indicates that, for a given pressure, the glass transition shows a weak temperature dependence on the cation length, occurring first for short-chain than for long-chain ionic liquids.     

Freezing and melting of binary mixtures confined in a nanopore

This paper reports on a Grand Canonical Monte Carlo study of the freezing and melting of Lennard–Jones Ar/Kr mixtures confined in a slit pore composed of two strongly attractive structureless walls. For all molar compositions and temperatures, the pore, which has a width of 1.44?nm, accommodates two contact layers and one inner layer. Different wall/fluid interactions are considered, corresponding to pore walls that have a larger affinity for either Ar or Kr. The solid/liquid phase diagram of the confined mixture is determined and results compared with data for the bulk mixture. The structure of the confined mixture is studied using 2D order parameters and both positional g(r) and bond orientational G₆(r) pair correlation functions. It is found that in the confined solid phase, both the contact and inner layers have a hexagonal crystal structure. It is shown that the freezing temperature of the Ar/Kr confined mixture is higher than the bulk freezing point for all molar compositions. Also, it is found that the freezing temperature becomes larger as the ratio α of the wall/fluid to the fluid/fluid interactions increases, in agreement with previous simulation studies on pure substances confined in nanopores. In the case of pore walls having a stronger affinity for Kr atoms (ε _Ar/W<ε _Kr/W), it is observed that both the contact and inner layers of the confined mixture undergo, at the same temperature, a transition from the liquid phase to the crystal phase. The freezing of Ar/Kr mixtures confined between the walls having a stronger affinity for Ar (ε _Ar/W?>?ε _Kr/W) is more complex: for Kr molar concentration lower than 0.35, we observe the presence of an intermediate state between all layers being 2D hexagonal crystals and all the layers being liquid. This intermediate state consists of a crystalline contact layer and a liquid-like inner layer. It is also shown that the qualitative variations of the increase of freezing temperature with the molar composition depend on the affinity of the pore wall for the different components. These results confirm that, in addition to the parameter α the ratio of the wall/fluid interactions for the two species, η=?_Ar/W/?_Kr/W, is a key variable in determining the freezing and melting behaviour of the confined mixture.     

Molecular dynamics of tert-butyl chloride confined in CPG studied by NMR

The molecular dynamics of tert-butyl chloride (TBC) confined to Controlled Pore Glass matrices of 25 and 7.5 nm were investigated by measuring NMR linewidths, lineshapes, and 1H and 2H spin--lattice relaxation times. The behaviour of confined TBC can be explained assuming that the guest molecules form two distinct phases; the surface-affected phase, composed of molecules located at the pore surface, and the bulk-like phase located at the centre of the pores. The bulk-like component of confined TBC, at the temperatures corresponding to the phase III, is characterized by two dynamically different subphases.     

Molecular dynamics studies of liquids using a Beowulf computer

Molecular dynamics is a popular methodology for investigating the properties of liquids. In this article, the historical development of the subject and its current status will be briefly reviewed. The different parallelisation strategies that are commonly used are discussed, highlighting their relative strengths and weaknesses. Particular attention is given to the software and hardware aspects of implementing these algorithms on the 'Beowulf class' of parallel computers. Finally, three different examples of parallel molecular dynamics studies on a Beowulf computer will be discussed, that are indicative of the range of potential applications.     

Dynamics of layering transitions in confined liquids

Multiple beam interferometry and video microscopy were used to investigate the layering transition of thin liquid films of 1-undecanol confined between atomically smooth mica surfaces. The expulsion of a molecularly thin lubricant layer was followed directly in two dimensions. Overall, the dynamics of the transition follows theoretical predictions based on two-dimensional hydrodynamics. Frequently, pockets of liquid remain trapped inside the contact area at the end of the transition. The trapped pockets undergo shape transformations to minimize elastic and interfacial energy.     

Molecular insight into structures of monocationic and dicationic ionic liquids in mica slits

Structures of ionic liquids (ILs) 1-decyl-3-methylimidazolium bis(trifluoromethanesulfonyl)azanide ([C₁₀mim][TFSA]) and 1-decyl-dimethylimidazolium bis(trifluoromethanesulfonyl)azanide ([C₁₀(mim)₂](TFSA)₂) in different-sized mica slits have been investigated using molecular dynamics simulations. Ion density and angular distributions for monocationic IL [C₁₀mim][TFSA] were analysed to elucidate the IL structures under different surface charges and especially their changes in the direction perpendicular to the surfaces. [C₁₀mim][TFSA] formes in bilayers, compatible with existing models of ILs with long alkyl chains. For dicationic IL [C₁₀(mim)₂](TFSA)₂, cations adjacent to the mica surface tend to stay parallel to the surface with both positively charged rings absorbed. While near the centre of the slit, dications show the weak tendency of orientation distribution, more random than [C₁₀mim]⁺ ions. Structures of [C₁₀(mim)₂](TFSA)₂ cannot be described by bilayer models. Additionally, the in-plane arrangement of [C₁₀mim][TFSA] is more ordered when K⁺ ions completely neutralise the negative charge of the mica surface, and [C₁₀mim]⁺ ions tend to be located in hexagonal mica lattices with two aluminium atoms in replacement of silicon atoms. [TFSA]^? ions are constrained by the neighbouring K⁺ ions absorbed onto mica lattices.     

Aggregates in mixtures of ionic liquids

Surfactant-like ionic liquids (C₁₆MIMCl, C₁₆MIMBF₄) aggregate in mixtures with another ionic liquid (EAN). Critical aggregate concentrations and estimations of object sizes are given and compared to aqueous systems and other relevant literature data. The investigated mixtures are stable up to more than 200 °C and can probably be used to extend the limited temperature range of water-based colloids.     

Molecular simulation study of triangle-well fluids confined in slit pores

Grand canonical Monte Carlo simulations are used to study the behaviour of triangle-well (TW) fluids with variable well widths confined inside slit pores. The effect of individual factors influencing the properties of confined fluids such as fluid–fluid interactions, pore size and pore wall–fluid interactions are obtained using simulations as it is difficult to experimentally determine the same. An interesting observation of this study is that inside the narrow pore of slit height h* = 5 at the high-pressure condition of P* = 0.8, for the TW fluid with long-range attraction or for the fluid at a low temperature for even a short-range attraction, the density profiles show layering such that there is a sticking tendency of the particles at centre, while there is a depletion of particles near the wall (as the layers at the centre have higher density peak heights than near the walls).     

Structural properties of ionic liquids

Two types of ionic liquids are considered: molten salts and aqueous solutions of strong electrolytes. It is shown that in both cases, modern theoretical methods are capable of making definite predictions about the spatial correlations between the ions themselves and between the ions and water molecules. It is further shown that the technique of neutron diffraction as applied to isotopically enriched samples allows, for the first time, detailed comparisons to be made between theory and experiment. The results of these comparisons are discussed and areas in which the theory is in urgent need of revision are identified.     

嵌段共聚物受限于软孔内的自组装

利用自洽场方法研究两嵌段共聚物受限于接枝均聚物链(聚合物刷)圆孔中的自组装相形貌.研究表明,当圆孔内径一定时,嵌段比f和聚合物刷C的体积分数φ_C是调控嵌段共聚物相形貌的主要因素,聚合物刷的弹性熵也起着重要作用.当f＝0.7时,在聚合物刷的浸润下,贴近刷表面处AB嵌段共聚物构成环层状结构,随着φ_C的减小这种结构会周而复始地出现.当f处于0. 关键词： 嵌段共聚物 圆孔受限 聚合物刷 自洽场     

Au诱导形成有序Si纳米孔阵列及其应用

以自组装聚苯乙烯小球(PS)单层膜为掩膜,利用Au对Si表面的催化氧化作用以及KOH溶液对单晶Si的各向异性腐蚀特性,在Si(100)面上制备了一系列尺寸小于100 nm有序可控的Si纳米孔阵列.扫描电镜(SEM)和原子力显微镜(AFM)等的测试结果显示:当PS小球溶液与甲醇溶液的体积比为9:11时,可形成大面积无缺陷的单层膜;但当体积比过大时,会导致类似双层膜结构的形成;而当体积比过小时,会诱导形成点缺陷和线缺陷.对PS小球及溅射Au处理过的Si晶片进行KOH溶液腐蚀,随着腐蚀时间变长,纳米孔的横向尺寸和深度增大,其形貌由圆形逐渐变为倒金字塔型,当腐蚀时间超过10 min,纳米孔阵列的有序性遭到破坏.采用离子束溅射技术在倒金字塔型纳米孔衬底上获得了有序Ge/Si纳米岛,而在圆形纳米孔衬底上获得了有序Ge/Si纳米环.进一步对有序Ge/Si纳米岛及纳米环的形成机理进行了解释.