Investigation of ionic diffusion on molten Li-KCI interface by molecular dynamics method

采用分子动力学(MD)方法研究了熔融Li(电极)-KCl(电解质)界面上离子的扩散行为.熔融界面上离子的扩散动力学通过离子质心的均方位移(MSD)和速度自相关函数(VACF)进行研究,扩散系数由MSD(t)函数线性区间的斜率和VACF(t)函数积分得到.模拟结果表明,在熔融的Li-KCl界面上,Li+离子在浓度梯度的驱动下穿过界面发生定向迁移,导致双电层的形成和外电路上电流的输出.Li+离子的扩散系数比K+和Cl-离子的大7～8倍,说明在界面上Li+离子是主要的载荷子,热电池的电荷传输机制主要与Li+离子的扩散运动有关.由Nemst-Einstein公式对电导率进行估算,由扩散到KCl层中的Li+离子产生的电导率约为0.4 S·cm-2,对应的电流密度估算值为3.27×105 A-cm-2.     

石墨烯/聚乙烯界面导热性能的分子动力学模拟

采用非平衡分子动力学模拟(NEMD)方法研究了石墨烯/聚乙烯纳米复合材料的界面导热性能,主要考察了石墨烯层数、尺寸对界面热阻的影响.研究结果表明:当石墨烯层数为一层时,界面热导为46.79 MW/(m~2K),随着石墨烯层数的增加,界面热导下降;但石墨烯层数超过四层后,界面热导趋于恒定接近39.00 MW/(m~2K);随着石墨烯尺寸的增大,石墨烯中较长波长声子被引发并对界面热传导起到主要的作用,最终导致界面热导逐渐增大.     

非离子表面活性剂在气液界面的分子动力学模拟

采用分子动力学方法研究了十二烷基醇聚氧乙烯醚系列非离子表面活性剂单分子层在气液界面的微观结构,并通过表面张力的计算考察了表面活性剂分子结构与性能的关系.研究结果表明,随着表面活性剂分子乙氧基基团个数的增加,模拟所得的表面张力的变化趋势与实验一致,所有分子的计算误差在5 mN·m-1以内.同时,随着乙氧基基团数目的增加,...     

多元体系油水界面上常见表面活性剂行为的分子动力学模拟

采用分子动力学模拟研究了以十二烷基苯磺酸钠(SDBS)为代表的阴离子型表面活性剂,以十二烷基三甲基溴化铵(DTAB)为代表的阳离子型表面活性剂,以壬基酚聚氧乙烯醚(NPE)为代表的非离子型表面活性剂,以十二烷基二甲基甜菜碱(Betaine)为代表的两性表面活性剂及空白实验.模拟了表面活性剂在油水界面上的行为,考察了表面活性剂分子与石油分子之间的径向分布函数(RDF)、石油分子在竖直方向的均方位移(MSD)、油水界面张力(IFT)、石油层与岩石层之间的相互作用能、石油层的相对浓度在竖直方向的分布及石油分子质心位置随模拟时间的变化关系等,讨论了不同表面活性剂的洗油性能.结果表明:(1)SDBS,NPE和Betaine分子初始状态下呈近似的规律排列,非极性端部分插入油相中,极性端延伸进入水相中;随后表面活性剂的极性端表现出聚集趋势,逐渐形成一个外部亲油内部亲水的一个胶束状粒子,粒子随模拟的进行逐渐融入到油层当中;DTAB从开始的近似规则排列逐渐变为无规排列,但是始终保持亲油端插入到油相中,亲水端位于油水界面上.(2)表面活性剂分子与石油分子之间的相互作用强弱顺序为Betaine≈DTABSDBSNPE.(3)由质心高度和动力过程中的图像截图分析,表面活性剂洗油效果的顺序为BetaineSDBSNPEDTABNone.模拟结果与实际的驱油结果一致,从分子层面上解释了不同表面活性剂洗油的规律.     

电解质溶液界面结构的分子动力学模拟研究

电解质溶液界面结构的研究不仅具有重要的理论意义, 而且具有一定的实用价值. 采用分子动力学模拟研究了LiCl, LiBr, LiI, NaI, KI, CsI水溶液中阴阳离子在1×105 Pa和300 K下的气液界面分布情况, 探讨离子水化与电解质溶液界面结构的关系, 并分析阳离子水化能力的强弱对共存阴离子在界面富集分布的影响. 通过对模拟结果的分析发现, 离子的水化能力越强, 就越能形成稳定的水化结构而处于本体相中, 水化能力越弱, 则越易在界面富集. 该机理合理地解释了离子在界面的分布现象, 阳离子水化能力一般较其共存阴离子强而处于本体相, 阴离子则趋向在界面处富集; 不同阴离子在界面的密度分布也与阴离子的水化能力相关, 阴离子水化能力越弱, 其在界面富集程度越高, 不同阴离子在界面的富集趋势为Cl－＜Br－＜I－; 阳离子水化能力的强弱也影响其共存阴离子在界面的富集程度, 阳离子的水化能力越弱, 其共存阴离子在界面的富集程度就越低.     

超临界L-J流体混合物互扩散性质的分子动力学模拟

采用分子动力学方法计算了超临界Lennard-Jones(L-J)流体混合物的扩散性质, 分析了超临界条件下二元L-J混合流体(Ar-Kr体系)中各组分的自扩散系数及Maxwell-Stefan互扩散系数随组分的变化情况. 结果表明, 与Darken公式能很好地应用于常规条件下的Ar-Kr体系不同, 超临界条件下的Maxwell-Stefan扩散系数明显大于Darken公式的预测值, 不同原子间的速度互相关函数不可忽略且为正相关.     

气液界面上阴离子表面活性剂单层膜的分子动力学模拟

用分子动力学方法研究了阴离子表面活性剂十二烷基硫酸钠(SDS)在气液界面上的结构和动力学性质. 选择单分子占有面积分别为0.45和0.68 nm²的两个模拟体系, 通过径向分布函数表征了单层膜的厚度, 并根据疏水链中碳原子与极性头中硫原子之间组成的矢量分布和取向函数, 对比了不同界面单层膜的有序排列情况. 结果表明在分子占有面积较小达到饱和吸附的情况下, 界面上的SDS具有较好的有序性. 通过计算气液界面附近水分子的扩散系数发现: 由于氢键和静电作用的影响, 界面区域内的水分子较本体溶液中的水分子有较弱的迁移能力.     

腐蚀介质在缓蚀剂膜中扩散行为的分子动力学模拟

采用分子动力学模拟方法,从缓蚀剂膜阻碍腐蚀介质粒子(H2O、H3O+和HCO3-)向金属表面扩散的角度,研究了4种1-R1-2-十一烷基-咪唑啉缓蚀剂(R1:羧甲基(A),羟乙基(B),氨乙基(C),氢(D))抑制碳钢CO2腐蚀的缓蚀机理,并对其缓蚀性能进行了理论评价.腐蚀介质粒子在不同缓蚀剂膜中的扩散系数、粒子与膜的相互作用能以及膜的自扩散性能的计算结果表明:4种缓蚀剂均可形成稳定的缓蚀剂膜,能有效阻碍腐蚀介质粒子向金属表面的扩散,达到抑制或延缓腐蚀的目的;随亲水支链(R1)极性的增加,缓蚀剂膜对腐蚀介质粒子扩散行为的抑制能力逐渐增强;同种缓蚀剂膜对正负离子H3O+和HCO3-比对中性的H2O分子具有更强的扩散抑制能力.综合计算及分析结果,4种缓蚀剂缓蚀性能的理论评价结果为ABCD,与文献实验结果吻合.     

气体在无定型聚异戊二烯中扩散的分子动力学模拟

采用分子动力学模拟方法研究了多个温度下氧气、氮气及甲烷在无定型顺式1,4-聚异戊二烯中的扩散系数。在模拟过程中,使用COMPASS力场作为分子力场。应用COMPASS力场的势能函数,聚合物的密度及玻璃化转变温度的计算结果与实验值有较好吻合。在278-378 K的温度范围内,通过3或1.5 ns时长的正则系综动力学模拟,计算了不同温度下氧气、氮气及甲烷的扩散系数。结果表明,根据爱因斯坦关系式计算得到的扩散系数与实验结果比较接近。对气体扩散系数与温度的关系进一步研究,发现在278-378 K温度范围内,甲烷的扩散系数随温度变化的半对数曲线图是非线性的,而氧气和氮气的扩散系数随温度变化的半对数曲线图是线性的。本文研究结果有助于理解温度对气体扩散的影响机制,并为高温下气体在天然橡胶中扩散系数的测定及天然橡胶热氧老化建模分析提供依据。     

阴离子表面活性剂在油水界面聚集的分子动力学模拟

用分子动力学方法模拟了油、水和阴离子表面活性剂组成的混合溶液从初始“均相”到“油水两相”分离的动力学过程, 研究了十二烷基苯磺酸钠(SDBS)在界面分离过程中的作用. 模拟发现, 油水两相能够在短时间内分离达到平衡, 形成一个明显的油水界面; 在SDBS存在情况下, 油水界面的分离时间随着SDBS浓度的增加逐渐增加, 达到平衡时SDBS会在界面处形成一个明显的界面膜, 并对油水界面处的水分子有限制作用. 模拟表明, 分子动力学方法可以作为实验的一种补充, 为实验提供必要的微观分子结构信息.     

Molecular dynamics simulation of ionic transport on molten Li–KCl interface

A molecular dynamics study is performed to determine the dynamics and transport properties of the ions on the molten interface between anode metal Li and electrolyte KCl. Radial distribution function of the ionic pair and the behavior of the mean‐square displacement (MSD) as a function of time (t) indicate that KCl and metal Li are in the molten state at 2,200 K in the canonical ensemble. The dynamics of the ionic transport are characterized by studying MSD for the centers of mass of the ions at different temperatures. Diffusion coefficient is evaluated from the linear slope of the MSD (t) function in the range of 0–500 ps. The MSD and diffusion coefficient of the Li⁺ ions are much larger than those of the Cl^? and K⁺ ions due to the difference in ionic characteristic. The transport process has been dominated by the Li⁺ ions on the molten interface and the Li⁺ ions are main charge carriers. The energy barrier of the Li⁺ ions transporting into the molten KCl is fitted to be 5.28 kcal/mol in the light of the activation model. The electrical conductivity of the Li⁺ ions transporting into the molten KCl are calculated from the Nernst–Einstein formula to be in the range of 0.2–0.3 S cm^?1. The current density resulted from the Li⁺ ions through the interface are estimated to be an order of 10⁶ A cm^?2, which may be the value corresponding to a larger concentration gradient of the Li⁺ ions. Simulated results at different temperatures show that the diffusion coefficient, conductivity and current density have increased with the temperature. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Application of molecular dynamics simulations in molecular property prediction II: diffusion coefficient

In this work, we have evaluated how well the general assisted model building with energy refinement (AMBER) force field performs in studying the dynamic properties of liquids. Diffusion coefficients (D) have been predicted for 17 solvents, five organic compounds in aqueous solutions, four proteins in aqueous solutions, and nine organic compounds in nonaqueous solutions. An efficient sampling strategy has been proposed and tested in the calculation of the diffusion coefficients of solutes in solutions. There are two major findings of this study. First of all, the diffusion coefficients of organic solutes in aqueous solution can be well predicted: the average unsigned errors and the root mean square errors are 0.137 and 0.171 × 10(-5) cm(-2) s(-1), respectively. Second, although the absolute values of D cannot be predicted, good correlations have been achieved for eight organic solvents with experimental data (R(2) = 0.784), four proteins in aqueous solutions (R(2) = 0.996), and nine organic compounds in nonaqueous solutions (R(2) = 0.834). The temperature dependent behaviors of three solvents, namely, TIP3P water, dimethyl sulfoxide, and cyclohexane have been studied. The major molecular dynamics (MD) settings, such as the sizes of simulation boxes and with/without wrapping the coordinates of MD snapshots into the primary simulation boxes have been explored. We have concluded that our sampling strategy that averaging the mean square displacement collected in multiple short-MD simulations is efficient in predicting diffusion coefficients of solutes at infinite dilution.     

Thermal diffusion and molecular weight calibration of poly(ethylene-co-vinyl acetate)s by thermal field-flow fractionation

Copolymer characterization is accomplished with respect to measurement of thermal diffusion coefficient (D_T) and molecular weight determination by thermal field-flow fractionation. The examined copolymers are the eight poly(ethylene-co-vinyl acetate)s [P(E-V)] having different compositions of vinyl acetate ranging from 25 to 70% and the molecular weight from 110,000 to 285,000, and three polyvinyl acetate standards as component homopolymer. The carrier solvents are tetrahydrofuran, toluene, and chlorobenzene which have different viscosities and thermal conductivities. Measured D_T values vary from 1.36 × 10^?8 to 5.97 × 10_?8 cm²/(s . K) which are dependent on the composition of copolymers and types of carriers. These values increase linearly with the increase of weight percent of vinyl acetate. It is possible to estimate D_T values of polyethylene from the extrapolated intercept in the plots of D_T vs. vinyl acetate wt % of copolymer. Tetrahydrofuran is found to be the appropriate carrier solvent for the separation of P(E-V) copolymers since D_T varies greatly with the increase of wt % in THF. Attempts are made to correlate the measured retention data with molecular sizes of copolymers for the construction of the molecular weight calibration curve. Good correlations (r² ≥ 0.931) are found in which D/D_T values of polymers vary inversely with the product of hydrodynamic volume by weight ratio of vinyl acetate. Based on this relationship, the unknown molecular weight of copolymer sample can be determined from component homopolymers for which standards are readily available. © 1995 John Wiley & Sons, Inc.     

The absorption and diffusion of polyethylene chains on the carbon nanotube: The molecular dynamics study

Classical molecular dynamics simulations have been used to investigate the absorption and diffusion behavior of polyethylene (PE) chains on the surface of the side‐wall of the carbon nanotube (CNT). Different degrees of polymerization from 50 to 80 at separate temperatures of 300, 400, 500, and 600 K are considered. Through the simulation, it is examined that the PE chains are absorbed on the surface of CNT and form stable composites with the nanotube as capsules. It is found that the most probable distance between the CNT and the C atoms in backbone of PE molecules only attribute to the temperature, and at T = 300 K, this distance is about 3.8 Å. Furthermore, the pattern of the composites mainly depends on the temperature and the length matching of the chains and the CNT. In particular, the PE chains keep approximately linear conformation, and extend along the axis of the CNT at the room temperature. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 46: 272–280, 2008     

Investigation of chiral recognition by molecular micelles with molecular dynamics simulations

Molecular dynamics simulations were used to characterize the binding of the chiral drugs chlorthalidone and lorazepam to the molecular micelle poly-(sodium undecyl-(L)-leucine-valine). The project’s goal was to characterize the nature of chiral recognition in capillary electrophoresis separations that use molecular micelles as the chiral selector. The shapes and charge distributions of the chiral molecules investigated, their orientations within the molecular micelle chiral binding pockets, and the formation of stereoselective intermolecular hydrogen bonds with the molecular micelle were all found to play key roles in determining where and how lorazepam and chlorthalidone enantiomers interacted with the molecular micelle.     

甲烷水合物导热机理的分子动力学模拟

甲烷水合物导热系数是甲烷水合物勘探、开采、储运以及其他应用过程中一个十分重要的物理参数.我们采用平衡分子动力学（EMD）方法Green-Kubo理论计算温度203.15～263.15K、压力范围3～100MPa、晶穴占有率为0～1的sI甲烷水合物的导热系数,采用的水分子模型包括TIP4P、TIP4P-Ew、TIP4P-FQ、TIP4P/2005、TIP4P/Ice.研究了主客体分子、外界温压条件等对甲烷水合物导热性能的影响.研究结果显示甲烷水合物的低导热性能由主体分子构建的sI笼型结构决定,而客体分子进入笼型结构后,使得笼型结构导热性能增强,同时进入笼型结构的客体分子越多,甲烷水合物导热性能越强.研究结果还显示在高温区域（T〉TDebye/3）内不同温度作用下,所有sI水合物具有相似的导热规律.压力对导热系数有一定影响,尤其是在较高压力条件下,压力越高,导热系数越大.而在不同温度和不同压力作用过程中,密度的改变对导热系数的增大或减小几乎没有影响.     

Investigation of the change in the electronic properties of FeF3 by the introduction of oxygen using a molecular orbital method

FeF₃ has attracted interest as a conversion‐reaction‐based positive electrode material in applications to lithium ion batteries. However, slow reaction kinetics is a major drawback due to its poor electrical conductivity. The electronic features of FeF₃ were examined using the DV‐Xα molecular orbital method. This article reports the effects of oxygen doping on the bonding characteristics and electrical conductivity. An analysis of the bond overlap population and spatial distribution of electrons showed that the Fe? O bond has a more covalent nature than the Fe? F bond. New energy levels were generated in the original band gap region through an interaction between the Fe3d and O2p orbitals with the introduction of oxygen. The electrical conductivity of FeF₃ is expected to be increased by the partial substitution of oxygen for fluorine due to the higher covalent character and the formation of new energy levels. © 2013 Wiley Periodicals, Inc.     

The effect of corrugation of pore wall on the thermal diffusion in nanopores by molecular dynamics simulations

In this work we have studied the effect of corrugation on the thermal diffusion (soret effect) in isotopic and non-isotopic fluid mixtures confined in a slit pore. We used a boundary driven non-equilibrium molecular dynamics to simulate thermal diffusion in Lennard–Jones (LJ) binary mixtures confined in structureless Steele 10-4-3 and atomistic Lennard–Jones pore walls. The results showed that for the isotopic mixture thermal diffusion factor for both wall types agrees and the corrugation of the LJ wall has no effect in isotopic mixture. However, for non-isotopic mixture confined in atomistic LJ pore the component with stronger attraction adsorbs more to the wall than the structureless Steele wall. The effect of corrugation of pore wall on the thermal diffusion is noticeable in narrow slit pore and mixture with large difference in molecular attraction parameter of components.     

Molecular dynamics study of the water/n-alkane interface

Molecular dynamics simulations on the interface between liquid water and liquid n-alkane (including octane, nonane, decane, undecane and dodecane) have been performed with the purpose to study the interfacial properties: (Ⅰ) density profile; (Ⅱ) molecular orientation; (Ⅲ) interfacial tension and the temperature effect on the interfacial tension. Simulation results show that at the interface the structures of both water and n-alkane are different from those in the bulk. Water has an orientational preference ...