Dynamics Studies on Molecular Diffusion in Zeolites

1 INTRODUCTION Zeolites have attracted much attention in both scientific and industrial areas for their special characteristics and multiple uses. Zeolites are composed of TO4 (T = Si, Al, Ga, P, etc.) tetrahedra which are connected with each other by sharing oxygen atoms to produce a complex and repetitive three-dimensional atom network with regular molecular dimension cavities joined by channels. Shape selectivity is the most important property of zeolites. Combined with the possib…     

Molecular Dynamics Simulation of Diffusion of Vitamin C in Water Solution

Under different temperatures and concentrations, the diffusion of Vitamin C (VC) in water solution was examined by molecular dynamics simulation. The diffusion coefficients were calculated based on the Einstein equation. The influences of temperature, concentration, and simulation time on the diffusion coefficient were discussed. The results showed that at higher temperature and lower concentration the normal diffusions appear relatively late, but the linear range of mean square displacement curves continues longer than that at lower temperature and higher concentration. At the same temperature, the normal diffusion time increases and the diffusion coefficient decreases as the simulation concentration increases. These simulation results are in good agreement with experiments. Analyses of the pair correlation functions of the simulation systems showed that hydrogen bonds are mainly formed between the hydrogen atoms of VC molecules and oxygen atoms of H₂O molecules, rather than between the O atoms of VC molecules and H atoms of H₂O molecules. The diffusion coefficient is higher as the interaction between water molecules and VC molecules is stronger when VC concentration is lower. The water in the model systems affects the diffusion of VC molecules by the short‐range repulsion of O(H₂O)‐O(H₂O) pairs and the non‐bond interaction of H(H₂O)‐H(H₂O) pairs. The short‐range repulsion of O(H₂O)‐O(H₂O) pairs is greater when VC concentration is higher, the diffusion of VC is weaker. The greater the non‐bond interaction of H(H₂O)‐H(H₂O) pairs is, the higher the VC diffusion is. It is expected that this study can provide a theoretical direction for the experiments on the mass transfer of VC in water solution.     

Molecular Dynamics Simulations of Ethene Diffusion in Orthorhombic and Monoclinic MFI

运用分子动力学（ＭＤ）模拟研究了３００Ｋ下乙烯在正交晶系和单斜晶系ＭＦＩ中的吸附、极化和扩散行为，获得了温度、能量、分子扩散速度、分子与沸石骨架间的相互作用、分子在沸石孔道中的吸附能以及分子的偶极矩等计算结果；同时发现，分子在直孔道和Ｚ型孔道的交叉部位的边界处极化最大，且在交叉部位的中心区域处出现的几率最大．     

氧气在聚丙烯内吸附和扩散的分子模拟

采用巨正则Monte Carlo和分子动力学模拟相结合的方法研究了氧气在不同聚合度的聚丙烯内的吸附和扩散. 模拟结果表明, 随聚丙烯聚合度的增加, 聚丙烯对氧气的吸附量逐渐增加, 而氧气在聚丙烯内的扩散系数减小; 当聚合度增大到一定程度时, 吸附量和扩散系数都趋于一稳定值. 随温度的升高, 氧气在聚丙烯内的吸附量减少, 而扩散系数增大. 本文还应用自由体积理论探讨了氧气在聚合物内扩散的机理, 发现氧气在聚丙烯内以空穴形式存在的自由体积之间扩散, 即氧气先在一个空穴内不停振动, 然后通过聚丙烯链段运动形成的通道跳跃到下一个空穴来完成扩散. 结果表明, 较高聚合度的聚合物材料在常温及低温下使用对于其在食品包装材料中的应用是有利的, 这为食品包装材料行业相关产品的应用开发提供了一定的指导和依据.     

Unveiling Carbon Dioxide and Ethanol Diffusion in Carbonated Water-Ethanol Mixtures by Molecular Dynamics Simulations

The diffusion of carbon dioxide (CO2) and ethanol (EtOH) is a fundamental transport process behind the formation and growth of CO2 bubbles in sparkling beverages and the release of organoleptic compounds at the liquid free surface. In the present study, CO2 and EtOH diffusion coefficients are computed from molecular dynamics (MD) simulations and compared with experimental values derived from the Stokes-Einstein (SE) relation on the basis of viscometry experiments and hydrodynamic radii deduced from former nuclear magnetic resonance (NMR) measurements. These diffusion coefficients steadily increase with temperature and decrease as the concentration of ethanol rises. The agreement between theory and experiment is suitable for CO2. Theoretical EtOH diffusion coefficients tend to overestimate slightly experimental values, although the agreement can be improved by changing the hydrodynamic radius used to evaluate experimental diffusion coefficients. This apparent disagreement should not rely on limitations of the MD simulations nor on the approximations made to evaluate theoretical diffusion coefficients. Improvement of the molecular models, as well as additional NMR measurements on sparkling beverages at several temperatures and ethanol concentrations, would help solve this issue.     

ITQ-3分子筛中二元混合分子扩散性质的模拟

采用平衡分子动力学方法先是模拟研究了纯组分Ar、SF_6以及CF4在ITQ-3分子筛中的扩散行为,结果表明,在窄孔道中Ar的扩散系数随负载的增加呈现先增加再减小的变化趋势,在宽孔道中则随负载逐渐减小。纯组分的大分子SF_6以及CF4只在y方向扩散,扩散系数与小分子在z方向的变化规律一致。随后,又模拟了二元混合组分SF_6和CF4在ITQ-3分子筛中的扩散行为,模拟结果与各自纯组分SF_6和CF4进行了对比,发现二元混合的SF_6和CF4整体变化趋势与其单组分保持一致,在y方向的扩散系数都是随着负载的增加而先增大后减小,但是两者之间还存在一种相互作用,使得在整个负载范围内扩散较快的CF4的扩散系数比纯组分的扩散系数相对小一些,而纯组分时扩散相对缓慢的SF_6在二元混合中的扩散系数则相对增大。     

柔性金属—有机骨架材料中已烷的扩散研究

本工作将MM3力场进行了扩展,使其可用于描述其它IRMOF材料的柔韧性;在此基础上,我们采用分子动力学模拟研究了柔性IRMOF-1和-16材料中已烷的扩散。 本文重点研究了温度和分子数对己烷自扩散系数、扩散机理,以及骨架柔性的影响。结果表明,分子数是影响扩散路径的重要因素。其次,IRMOF-16的柔性强于IRMOF-1。 工作的结论有助于进一步研究链状分子在柔性MOF材料中的扩散。     

Microscopic Diffusion Dynamics of Silver Complex‐Based Room‐Temperature Ionic Liquids Probed by Quasielastic Neutron Scattering

Quasielastic neutron scattering is used to probe the microscopic diffusion dynamics of the hydrogen‐bearing cations of two different silver complex‐derived room‐temperature ionic liquids, [Ag(propylamine)₂⁺][Tf₂N^?] (Tf=trifluoromethanesulfonyl) and [Ag(1‐pentene)⁺][Tf₂N^?]. In the temperature range from 300 to 340 K, analysis of the scattering momentum transfer dependence of the data provides evidence for three distinct diffusion components. The slowest component describes the long‐range cationic translational diffusion. A possible link between the microscopic diffusion parameters and the structural features of the cations comprising these two ionic liquids is discussed.     

石墨烯在甘油/尿素剥离液中的稳定行为的分子动力学模拟研究

现有的实验方法很难实时观测到石墨烯在液相剥离溶剂中的结构演变,尤其是石墨烯稳定的微观机理尚不明确.本工作通过分子动力学方法,模拟了多层石墨烯和U型石墨烯在不同的物质的量比下的甘油/尿素溶剂中的结构变化,研究剥离液对石墨烯稳定性的影响.结果 表明,多层石墨烯在不同溶剂体系中的稳定性差异不显著;而U型石墨烯在各溶剂体系的稳...     

水化锂蒙脱石层间结构和动力特征的分子模拟

利用分子力场和分子动力学(MD)的方法研究了Li-蒙脱石的结构构型, 层间阳离子的水化行为、水分子的结构特征以及它们的扩散性质. 分子力场构型优化结果表明: Li-蒙脱石的层间距、体积和密度与层间水含量有关; MD模拟的动画轨迹显示Li-蒙脱石层间Li^＋的位置与层间电荷位置有关. 均方根位移和自扩散系数的计算结果表明: 层间阳离子、水分子在Li-蒙脱石一、二层水合物中的扩散受到上下粘土片表面的限制, 在三层水合物中开始离开粘土层面向其它方向快速扩散. 径向分布函数及其结构因子的分析结果表明Li^＋在一、二、三层水合物中有不同的水合层; 层间水分子的结构特征说明其在蒙脱石层间有水合水分子和自由水分子之分, 且它们的比值在一、二和三层水合物中有所不同.     

受限于纳米碳管中的乙醇分子的结构和扩散的分子模拟研究

利用分子模拟研究了常温常压下受限于(8,8) (管径1.081 nm)和(15,15) (管径 2.035 nm)单壁纳米碳管中的乙醇分子. 对受限分子的径向密度分布和氢键等静态性质以及扩散性质进行了分析. 结果显示在管内乙醇分子的平均氢键数目和主体相一致. 乙醇分子在(8,8)碳管内具有高度有序的结构, 而在(15,15)碳管内由于空间的增大导致结构有序度的降低, 其中分子取向已呈随机分布. 进一步对扩散系数的分析发现, 在管内乙醇分子的轴向扩散系数低于主体相, 特别在(8,8)碳管内乙醇分子几乎丧失了轴向扩散能力.     

Molecular Dynamics Simulation of the Fracture in Polymer‐Exfoliated Layered Silicate Nanocomposites

Summary: The MD technique was used to investigate the fracture behavior in fully exfoliated layered silicate (nanoplatelet)‐polymer nanocomposites. MD results reveal that the addition of the nanoplatelets can improve the fracture strength of polymers. The interactions between the surface of the nanoplatelets and the segments of the polymer, and the relaxation time of polymer chains have significant influences on the fracture strength of the polymer. For polymers with T_g below room temperature, such as polyurethane, or close to room temperature, such as nylon, the nanoplatelets are always working for the enhancement of the mechanical properties. However, for polymers with T_g above room temperature, such as epoxy and polystyrene, the addition of the nanoplatelets is not working well for toughening these polymers. If the nanoplatelets are to enhance the mechanical properties of these polymers, it is necessary to build up a stress relaxation interface between the polymer and the nanoplatelet in order to reduce the effect of the difference between the relaxation time of nanofillers and that of polymers.

Force per area versus distance curves as a function of the difference of the relaxation times of the nanoplatelets and polymer chains.     

Modeling of Structural,Energetic, and Dynamic Properties of Few‐Atom Silver Clusters Embedded in Polynucleotide Strands by Using Molecular Dynamics

This work concerns the study of the structural, energetic, and dynamic properties of fluorescent systems composed of silver clusters stabilized by polynucleotide strands. To do so, classical interaction potentials relative to silver, neutral and cationic, were introduced in the AMBER force field. Molecular dynamics simulations allowed analysis of the nature and force of the interactions between the various parts of the nucleic oligomers and the silver clusters. Conformational analyses were necessary to explore the flexibility of the supramolecular assemblies, specifically by radial distribution functions and Ramachandran‐type maps.     

锂-,钠-,钾水化蒙脱石层间结构的分子动力学模拟

运用分子动力学(molecular dynamics, MD)方法分别研究了含有32, 64和96个水分子的Li-, Na-, K-蒙脱石层间阳离子与水分子的位置和结构. 计算结果表明蒙脱石层间阳离子位置与四面体和八面体电荷位置及离子的大小有关. 一层水合物中可以观察到三种阳离子都能和四面体电荷与八面体电荷位置分别形成内、外配位作用. 二层水合物中, 仍然可以观察到Li^＋和Na^＋与电荷位置的配位作用, 但是已经开始向层中其他方向扩散, 而K^＋仍然在粘土的表面附近. 三层水合物中, Li^＋, Na^＋开始从电荷位置和表面分离, K^＋也开始向层间其他方向扩散. 水分子在所有三种水合物中都分散于层间各个方向. 径向分布函数的分析结果表明层间三种阳离子组织水分子的能力不同, 水化作用随着阳离子半径的增大而减弱; 此外层中水分子的聚合程度随着水分子的增加而加强, 水分子的结构也不同于模拟的液体水分子的结构; 说明蒙脱石层间阳离子的溶剂化作用对水分子的组织起着重要的作用.     

ZIF-8材料中CH4/H2吸附与扩散的分子模拟

采用平衡分子动力学和巨正则系综蒙特卡洛模拟方法对多孔沸石咪唑酯骨架材料(ZIF-8)中CH4和H2分子的吸附与扩散特性进行了比较研究.结果表明,采用柔性力场能够很好地复制ZIF-8在不同压力或温度条件下的晶体结构,也能准确地计算不同温度下CH4和H2分子在ZIF-8中的扩散系数,特别是高温下CH4分子因能够摆脱ZIF-8骨架笼口的空间限制而使其扩散系数出现大幅提升.同时,该力场也能粗略地模拟CH4和H2分子在ZIF-8中的等温吸附曲线,通过自编程序得到吸附和扩散平衡时CH4和H2分子在ZIF-8单元晶胞内的几率密度分布数据,并利用VMD软件可视化.结果显示CH4和H2分子在ZIF-8中的优先吸附位置均在大孔中心靠近咪唑环的区域,但CH4分子的优先吸附位置有两个不同层次,而H2分子的优先吸附位置只有一个层次,说明CH4和H2分子在ZIF-8中存在着不同的吸附机理.     

Structure of Polymer Brushes in Cylindrical Tubes: A Molecular Dynamics Simulation

Summary: Molecular dynamics simulations of a coarse‐grained bead‐spring model of flexible macromolecules tethered with one end to the surface of a cylindrical pore are presented. Chain length N and grafting density σ are varied over a wide range and the crossover from “mushroom” to “brush” behavior is studied for three pore diameters. The monomer density profile and the distribution of the free chain ends are computed and compared to the corresponding model of polymer brushes at flat substrates. It is found that there exists a regime of N and σ for large enough pore diameter where the brush height in the pore exceeds the brush height on the flat substrate, while for large enough N and σ (and small enough pore diameters) the opposite behavior occurs, i.e. the brush is compressed by confinement. These findings are used to discuss the corresponding theories on polymer brushes at concave substrates.

Snapshot picture of a brush grafted inside of a cylinder, for , , displaying different chains in distinct colors in order to be able to distinguish them. Top shows a side view of the cylinder, and the lower part a view of the cross‐section. Note that the particles forming the cylindrical wall are not displayed.     

超临界CO2中甲醇和乙醇无限稀释扩散系数的分子动力学模拟与实验测定

采用分子动力学模拟(MD)方法对甲醇和乙醇在超临界二氧化碳中的无限稀释扩散系数进行了模拟计算, 并应用泰勒分散理论, 采用超临界色谱仪对模拟结果进行了实验验证. 模拟计算值与实验值较吻合, 且变化规律基本一致, 表明采用这种新方法可以准确有效地预测超临界体系的扩散性质, 能够方便地应用于工程设计.     

诺氟沙星-DNA复合物的分子动力学模拟

采用分子模建的方法构建了诺氟沙星-DNA复合物的初始结构, 通过2 ns的分子动力学(MD)模拟研究表明: 诺氟沙星能够和双螺旋d[ATATCGATAT]₂形成稳定的复合物, 药物分子可紧密结合在DNA的小沟区域, 并且能够与DNA的鸟嘌呤碱基形成两个稳定的氢键. 在分子水平上提供了诺氟沙星直接与双螺旋DNA相互作用的结构及复合物的动态变化情况.     

Dynamics of water and sodium in gels under salt‐induced phase transition

Synthetic and biological gels undergo a sharp volume phase transition when subjected to a variety of environmental changes. Water and ion dynamics within swollen and compact phases are critical for understanding fundamental concepts in cellular (specifically neuronal) biophysics, for models of bound, free, or ordered water in complex environments; and for practical applications such as the design of gels for drug release, biomimetics, sensors, or actuators. In this work, we find, for the first time, basic physical parameters that shed light on the interaction of gels with water and electrolytes, across a volume phase transition. Water within a gel can be separated into bound and free populations with high exchange rate. We show that free water dynamics in compact gels are the same as those in pure water. Bound water was found to comprise a single layer around the polymers in both phases, with a correlation time three orders of magnitude higher than that of free water. Most importantly, salt‐induced phase transition was found to be different from a standard coil‐globule transition (e.g., temperature‐induced), with no rejection of bound water as the gel compacts. © 2015 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53, 1620–1628     

Sorption and Diffusion of Propylene and Ethylene in Heterophasic Polypropylene Copolymers

Summary: Sorption experiments of ethylene and propylene in different polypropylene powder samples, both homopolymer and heterophasic copolymers with different rubber content, have been carried out in a high-pressure magnetic suspension balance at 10 bars pressure and 70 °C. The gross solubilities measured can be well correlated with the rubber content of the polymer samples. Solubility of ethylene and propylene in the rubber phase differ from solubility in the amorphous fraction of the homopolymer, especially the concentration ratio of propylene to ethylene differs significantly between rubber phase and amorphous fraction of the homopolymer. From the slope of monomer uptake, information on kinetics of mass-transfer can be gained. No significant differences were observed in terms of mass-transfer for ethylene and propylene. With increasing rubber content, effective diffusion coefficients increased slightly. By combined sorption studies with powder samples and compressed films, information about both effective diffusion coefficients and the effective length scale of diffusion could be gained. It could be shown, that the particle radius is not the characteristic length of diffusion in the studied powder samples. Mass transfer of nearly all samples could be described by a constant diffusion length of 120 to 130 µm, independently on particle size. This indicates that the effective scale of diffusion in polymer particles is in between microparticle and macroparticle scale used in classical particle modeling.