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

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

端羟基聚丁二烯/增塑剂共混物相容性的分子动力学模拟

固体推进剂和炸药的力学性能在很大程度上依赖于配方中高分子粘结剂与增塑剂的相容性. 本文对相容和非相容两种体系进行了分子动力学(MD)模拟, 以考察分子模拟方法的实用性. 为预测固体推进剂中端羟基聚丁二烯(HTPB)与增塑剂癸二酸二辛酯(DOS)、硝化甘油(NG)的相容性, 采用MD模拟方法在COMPASS力场下, 对HTPB、DOS、NG和共混物HTPB/DOS、HTPB/NG的密度、内聚能密度及溶度参数等进行了模拟计算. 通过比较溶度参数差值(△δ)的大小、分子间径向分布函数和模拟前后体系密度变化情况均可以预测HTPB/DOS属于相容体系,而HTPB/NG属于不相容体系, 与实验结果一致. 径向分布函数分析同时揭示了HTPB/增塑剂组分之间的相互作用及本质. 本文的模拟方法可以作为预测聚合物与增塑剂相容性的有利工具, 也可以为固体推进剂和炸药的配方设计提供理论指导.     

分子动力学模拟无定形BCN体系结构特征

为研究无定形BCN材料的原子结构特征, 采用分子动力学方法模拟了高温下无定形BCN体系的动力学行为, 计算并分析了引入C原子对体系的径向分布函数(RDF)和配位数的影响. 模拟结果表明, C原子对BCN体系的RDF曲线的峰位置及形状影响很小, 但使其偏径向分布函数(PRDF)的一些峰向右偏移; 引入的C原子部分占据了原来N原子的位置, 使得B-N的配位数降低. 在模拟条件下, 未观察到BN相与C相分离.     

Mechanical Properties of Interlocked-ring Polymers: A Molecular Dynamics Simulation Study

Interlocked-ring polymers, also known as polycatenanes, possess an interesting molecular architecture. These polymers are composed of many interlocked rings in a linear chain. The topological constrain between neighboring rings distinguishes the interlockedring polymer from its linear counterpart. Here we present extensive molecular dynamic simulations on the interlocked-ring polymers and analyze the static properties of the polymer. By applying external forces to the polymer, we also study the force-extension curves of the polymer, which provides rich information about the mechanical properties of the interlockedring polymers.     

Dynamics of a Spreading Nanodroplet: A Molecular Dynamic Simulation

The spreading of polymer nanodroplets upon a sudden change from partial to complete wetting on an ideally flat and structureless solid substrate has been studied by molecular dynamic simulations using a coarse‐grained bead‐spring model of flexible macromolecules. Tanner's law for the growth of the lateral droplet radius {R(t) ∝ t^0.1} is found to hold as long as the droplet does not disintegrate into individually moving chains. The data for the contact angle θ following from Tanner's law correspond to a dependence on time {θ(t) ∝ t^−0.3}. Our analysis of the mean square displacements of the polymer centers of mass reveals several dynamic regimes during the process of spreading. PACS numbers: 68.10.Gw, 05.70.Ln, 61.20.Ja, 8.45.Gd.

Molecular dynamics results for the average mean square displacement of all polymer chains plotted vs. time for a broad range of values for ε_wall.     

Polymer Structures and Glass Transition: A Molecular Dynamics Simulation Study

Full atomistic molecular dynamics (MD) simulations on five polymers with different chain backbone (C—C, Si—O, and C—O) and different side groups (—H, one —CH₃, and two —CH₃) are performed to study the effects of chain flexibility and side groups on the glass transition of polymers. Molecular dynamics simulations of NPT (constant pressure and constant temperature) dynamics are carried out to obtain specific volume as a function of temperature for polyethylene (PE), poly(propylene) (PP), polyisobutylene (PIB), poly(oxymethylene) (POM), and poly(dimethylsiloxane) (PDMS). The volumetric glass transition temperature has been determined as the temperature marking the discontinuity in slope of the plots of V–T simulation data. Various energy components at different temperatures of the polymers are investigated and their roles played in the glass transition process are analyzed. In order to understand the polymer chain conformations above and below the glass transition temperature, dihedral angle distributions of polymer chains at various temperatures are also studied.     

Self-Assembly of NaOL-DDA Mixtures in Aqueous Solution: A Molecular Dynamics Simulation Study

The self-assembly behaviors of sodium oleate (NaOL), dodecylamine (DDA), and their mixtures in aqueous solution were systematically investigated by large-scale molecular dynamics simulations, respectively. The interaction mechanisms between the surfactants, as well as the surfactants and solvent, were revealed via the radial distribution function (RDF), cluster size, solvent-accessible surface area (SASA), hydrogen bond, and non-bond interaction energy. Results showed that the molecules more easily formed aggregates in mixed systems compared to pure systems, indicating higher surface activity. The SASA values of DDA and NaOL decreased significantly after mixing, indicating a tighter aggregation of the mixed surfactants. The RDF results indicated that DDA and NaOL strongly interacted with each other, especially in the mixed system with a 1:1 molar ratio. Compared to van der Waals interactions, electrostatic interactions between the surfactant molecules were the main contributors to the improved aggregation in the mixed systems. Besides, hydrogen bonds were found between NaOL and DDA in the mixed systems. Therefore, the aggregates in the mixed systems were much more compact in comparison with pure systems, which contributed to the reduction of the repulsive force between same molecules. These findings indicated that the mixed NaOL/DDA surfactants had a great potential in application of mineral flotation.     

RbCl熔解的分子动力学模拟研究

采用等压分子动力学模拟方法，研究了从晶相到液相不同温度下ＲｂＣｌ体系的结构和性质，等压模拟咄等容怀致的平衡性质和结构特征，键序参数计算，不计及导热过程情形下，体系在振荡驰豫时间量级内完成熔化过程，由于不需要密度数据，等压模拟有望发展成为材料设计中的一种手段。     

液体和非晶态NiAl3合金结构的从头算分子动力学模拟

应用从头算分子动力学方法模拟了液体以及淬冷形成的NiAl₃合金体系,得到了它们的对相关函数、结构因子、键对分析信息.结果分析表明,在淬冷条件下得到的体系呈现非晶态性质,且非晶态结构类似于液态NiAl₃合金的结构,可以用液体结构近似描述非晶态性质.还进行了电子结构分析,得到液体NiAl₃合金的电子态密度和电荷分布.在液体镍铝合金中,镍为电子受体,部分电子由铝向镍转移,支持了Candy等人的XPS实验结果.镍铝间强烈作用,形成带有弱共价键性质的金属键.镍在合金中相当分散,这能部分解释由淬冷形成的NiAl₃合金制得的骨架镍催化剂活性增强的原因.     

液态金属结构变化的分子动力学模拟研究

随着计算机技术的飞速发展，已有可能将分子动力学方法应用于液态金属微观结构组态瞬时变化的定量模拟研究，使这一极为复杂的微观过程呈现出一幅十分清晰的物理图象，并已取得了许多重要的成果［1－4］．本文在作者原有工作的基础上［2－4］，对液态金属Al的结构变化进行了分子动力学计算机模拟研究，得到一幅非常清晰的关于金属熔融后继续升温过程中结构组态如何逐步变化的图景：随着温度的升高，其与高有序度相关的键型数目越来越少，而与无序度相关的键型数目却不断增加，即系统的总趋势是无序度增加．但却发现，即使在达到1773K（1．9…     

核糖核酸酶Sa表面水和尿素分子的分布和动力学行为的分子动力学模拟

利用分子动力学模拟研究了在不同尿素浓度下,核糖核酸酶Sa(RNase Sa)表面水和尿素分子的分布和动力学行为。 结果表明,尿素分子可与RNase Sa酶形成较强的相互作用,并取代其表面的水分子而富集在蛋白质表面。 尿素分子更倾向与RNase Sa酶的疏水残基作用,与RNase Sa酶主链形成氢键的能力更强。 尿素分子的平动和转动远远慢于水分子的平动和转动。 RNase Sa酶表面水分子的平动和转动随着尿素浓度增加而逐渐变慢,但RNase Sa酶表面尿素分子的动力学并不依赖于尿素浓度变化。 本研究中明晰的RNase Sa酶表面水和尿素分子分布和动力学有助于理解水和尿素分子对蛋白质稳定性的影响。     

纳米器件的分子动力学模拟

分子动力学广泛应用于分子体系, 但受限于计算能力, 一般难以用于纳米器件的研究. 本文采用自主开发的超大规模分子动力学仿真程序(NanoMD), 构建了原子级的纳米齿轮模型, 并以此为代表实现了对具有高速转动特性的纳米器件的模拟. 通过位错和缺陷分析, 确定了纳米器件在高速转动过程中的应力分布以及失效机制, 并明确了以极限弹性转速为依据的材料强度衡量模式. 研究发现纳米器件在极限转速方面存在明显的尺寸效应: 随着器件直径的减小而单调增大, 随着轴径的缩小而先增大后减小.     

A Study About Nanocluster Deposition of Thin-film Formation by Molecular Dynamics Simulation

The thin-film growth has been confirmed to be assembled by an enormous number of clusters in ICBD method. In sequence of clusters’ depositions proceeds to form the thin-film to understand quantitatively the interaction mechanisms between the cluster atoms and the substrate atoms, we use molecular dynamics simulation with EAM potential. The quantitative of flatness of deposition and percent of disordered atoms were proposed to evaluate the property of thin-film. In this simulation, three different Co cluster sizes of 55, 70, and 100 atoms with different velocities (100 up to 800 m/s) were deposited on a Al(0 0 1) substrate whose temperatures were set between 300 and 500 K. The simulations begin at specific equilibrium temperature of clusters and the substrate. The simulations are performed at different temperatures of the clusters and substrate and for different sizes of clusters. We showed that the percent of disordered atoms of substrate are affected by the cluster size and velocity of the clusters. Temperature dependence of the number of disordered atoms for different cluster’s velocity was observed. We investigated the effect of cluster size and initial velocity of cluster on the value of flatness.     

氧气在碱溶液中结构与扩散系数的分子动力学模拟

采用分子动力学方法对不同温度(25-120 ℃)及碱浓度(1:100-1:5, 摩尔比)下NaOH和KOH溶液中的氧气进行了模拟. 本文考察了NaOH及KOH溶液中溶剂-溶剂、氧气-溶剂及氧气-溶质的径向分布函数, 并采用爱因斯坦方程从均方位移曲线中计算得到了氧气及溶质离子的扩散系数. 结果显示随着碱浓度的升高, 氧气扩散系数逐渐减少; 在相同条件下, 氧气在NaOH溶液中扩散系数小于在KOH溶液中的扩散系数; 溶质离子扩散系数的变化规律与氧气一致. 通过与现有实验结果对比, 发现了分子动力学方法的可靠性及用于研究实验受限领域的优越性.     

热辅助存储磁盘硅掺杂非晶碳薄膜氧化的ReaxFF反应力场分子动力学模拟

Heat-assisted magnetic recording (HAMR) is one of the promising ways to extend the magnetic recording area density to 1 Tb·in^-2 in hard disk drives (HDDs).High temperature induced by laser heating can cause carbon overcoat (COC) oxidation.Reactive molecular dynamics (MD) simulations are performed to investigate the oxidation process of silicon-doped amorphous carbon (a-C:Si) films for HAMR application.The atomic details of the structure evolution and oxidation process are investigated, and, the oxidation mechanism of the a-C:Si film is clarified.The effect of the duration of laser irradiation on the oxidation of the a-C:Si film is investigated.The oxidation occurs during heating and the beginning of cooling process.Both volume expansion during heating process and cluster of carbon atoms during cooling process increase the rate of sp² carbon.Because of the decrease in the amount of unsaturated silicon atoms and low diffusion coefficient of atomic oxygen, the oxidation rate of the a-C:Si film decreases with laser irradiation cycles.The molecular oxygen is the oxidant due to surface defect of a-C:Si film.The atomic strains break the O-O bonds in Si-O-O-Si linkages and rearrange the surface oxide layers, and process the oxidation of the a-C:Si film.     

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.     

Order and Phase Behavior of a Cylinder Forming Diblock Copolymers and Nano‐Particles Mixture in Confinement: A Molecular Dynamics Study

We study a coarse grained model of cylinder forming diblock copolymers and nano‐particles (NPs) mixture confined between Lennard–Jones hard walls. Two models for non‐selective interactions between monomers and NPs are applied. In the case of purely repulsive interactions between NPs and monomers (athermal case) strong segregation of NPs at the film surfaces and the formation of droplets of particles inside the copolymer film can be observed. For weakly attractive interactions between NPs and monomers (thermal case) formation of droplets of particles disappears and segregation on the film surfaces depend on temperature. The uptake of NPs by the copolymer film in the thermal case displays a non‐monotonic dependence on temperature which can be qualitatively explained by a mean‐field model. In both cases of non‐selective interactions NPs are preferentially localized at the interface between the microphase domains.

     

Study of Imidazolium Ionic Liquids:Temperature-dependent Fluorescence and Molecular Dynamics Simulation

The steady-state fluorescence spectra and molecular dynamics simulations were explored to investigate the temperature dependent organization in some imidazolium ionic liquids:1-butyl-3-methylimidazolium hexafluo-rophosphate([bmim][PF6]),1-ethyl-3-methylimidazolium ethylsulfate([emim][EtSO4]) and 1-butyl-3-methylimida-zolium tetrafluoroborate([bmim][BF4]).The pure room temperature ionic liquids(ILs) exhibit a large red shift at more than an excitation wavelength of around 340 nm,which demonstrates the hetero...     

Molecular Dynamics Simulation Study of ClF in Water: Halogen Bonding Interaction in Liquid

Does the halogen bonding interaction co-exist in liquid when it competes with the hydrogen bonding interaction? The classical molecular dynamics simulations for the solvation properties of ClF molecule in water are performed with the Lennard-Jones plus Coulomb electrostatic potential parameters that are optimized with ab initio interaction energy calculations for the pre-reactive H₂O…ClF complex. We find that the halogen bonding interactions occur between O and Cl atoms and have the comparable strength and population with respect to the hydrogen bonding interactions of Cl…H.     

异型碳纳米管储氢性能的分子动力学模拟研究

采用分子动力学(MD)方法对三种理想的Y型碳纳米管[记为Y(4,4), Y(6,6), Y(10,0)]和三种L型碳纳米管[记为L(9,0), L(6,6), L(10,0)]之储氢性能进行了模拟研究, 并与相应的直线型碳纳米管的储氢能力进行了比较, 同时考察了温度、碳纳米管的直径和螺旋性以及缺陷的位置和大小对异型碳纳米管储氢性能的影响. 结果表明, 在室温和低温条件下, 异型碳纳米管的储氢量高于直线型碳纳米管的储氢量, 且其储氢量大小随温度的降低和碳管直径的增大而增加, 椅式碳纳米管的储氢性能优于齿式碳纳米管, 而缺陷的位置和大小对异型碳管之储氢性能的影响则因碳管的形貌和直径的大小不同而存在差异.