Chemical potentials and phase equilibria of Lennard-Jones chain fluids

Computer simulations (molecular dynamics) were performed for ensembles of flexible tangent Lennard-Jones chains consisting of n sites (n = 1, 2, 4, 8, and 16). From these simulations, the orthobaric liquid and vapour densities were calculated not only with the traditional method of simulating a liquid film in coexistence with vapour, but also using the rigorous thermodynamic condition of satisfying the chemical potential equality between the phases in equilibrium. The agreement with literature data, as far as such exist, is excellent.     

Synthesis of Cu nanoparticles by condensation from the gas phase

A series of computer experiments was conducted by the raw material evaporation and subsequent condensation to determine the most efficient regimes of copper nanoparticles synthesis. By variation of the cooling rate of the metal vapor the formation of Cu clusters were studied. The computer experiments showed the influence of different final temperatures on the shape of the resultant particles. This allowed to determine the conditions for a predominant formation of worm-like clusters of a spherical shape in the gas phase synthesis.     

固-液界面热整流现象的气体运动论分析

本文采用非平衡态分子动力学模拟方法,研究了固-液界面体系中的热整流现象。系统的两端为固体,而流体处于中间夹层,两端固-液界面的势能参数不同。模拟计算结果表明,此系统实现热整流的原因是左右两端固体壁面对流体分子吸附特性的不对称性。基于气体运动论的分析表明,高低温固体壁面所吸附流体分子数、温度协调系数和系统压强在正反向传热过程中均存在一定差异,从而形成了热整流效应;并且两端固-液界面差异性越大,热整流现象越明显。     

金属纳米薄膜在石墨基底表面的动力学演化

采用分子动力学方法研究了金属Au和Pt纳米薄膜在石墨(烯)基底表面的动力学演化过程,探讨了金属薄膜和石墨(烯)基底间的相互作用对金属纳米薄膜在固态基底表面的去湿以及脱附的动力学演化的影响.研究结果表明,在高温下,相同层数的Au和Pt纳米薄膜在单层石墨基底表面上存在不同的去湿现象,主要表现为厚度较小的Pt纳米薄膜在去湿过程中有纳米空洞形成,而同样厚度的Au薄膜在去湿过程中没有形成空洞.Au和Pt两种金属薄膜在高温下都去湿形成纳米液滴,这些液滴最终都以一定的速度脱离基底.在模拟的薄膜厚度范围内(0.2—2.3 nm),Au和Pt纳米液滴脱离基底的速度随厚度增加表现出不同的变化规律.Pt纳米液滴的脱离速度随薄膜初始厚度的增加先增加后减少,而Au脱离速度随厚度的增加先减少,达到一个临界厚度后脱离速度突然迅速增加.利用薄膜与基底间相互作用的不同导致去湿过程中的黏滞耗散不同,定性分析了这种变化规律的原因.此外,进一步研究还发现金属液滴的脱离时间与薄膜厚度和模拟温度的依赖关系,发现脱离时间随薄膜厚度的增加而增加,随模拟温度的升高而减小.这些研究结果可以为金属镀膜、浮选、表面清洁、器件表面去湿等工业生产过程提供理论指导.     

Can a difference in molecular weights cause an eruption in a driven flow of self-organizing immiscible system?

Driven flow of a non-equilibrium non-conservative (NENC) system with a mixture of immiscible particles (A,B of molecular weight M_A, M_B) exhibits self-organizing patterns (segregation, phase-separation, etc.) in steady-state. The flow response (v) of mass flux density (j) to bias (H), in steady-state is found to be sensitive to molecular weight ratio (α = M_B/M_A). While the flux density (j) responds linearly to bias for both components (A, B) at α = 1, onset of eruptive response occurs at extreme bias (H ↦ 1) at α > 1 where v ↦∞ for heavier (B) and v ↦- ∞ for lighter (A) constituents. Difference in molecular weights (M_A, M_B) is not only critical to eruptive flow but also in controlling the flow response prior to this crossover.     

A molecular dynamics study of the structural change differences between Au225 and Au369 clusters on MgO surfaces at low temperature

The differences in structural change between Au 225 and Au 369 clusters with their(111) facets supported on MgO(100) surfaces at 5 K are studied by using molecular-dynamics simulations with the atomic interchange potentials of the Au/MgO interface.The parameters are obtained from the ab initio energies using the Chen-Mo¨bius inversion method.Analyses of the pair distribution functions show that the two Au clusters use different deformation processes to adjust the distances between the interface atoms,owing to the misfit between the atom distances among the clusters and the substrates.The local structural changes are identified by atomic density profiles.     

Friction and normal forces of model friction modifier additives in simulations of boundary lubrication

ABSTRACT

Interaction forces between solid surfaces are often mitigated by adsorbed molecules that control normal and friction forces at nanoscale separations. Molecular dynamics simulations were conducted of opposing semi-ordered monolayers of united-atom chains on sliding surfaces to relate friction and normal forces to imposed sliding velocity and inter-surface separation. Practical examples include adsorbed friction-modifier molecules in automatic transmission fluids. Friction scenarios in the simulations had zero, one, or two fluid layers trapped between adsorbed monolayers. Sliding friction forces increased with sliding velocity at each stable separation. Lower normal forces were obtained than in most previous nanotribology molecular simulations and were relatively independent of sliding speed. Distinguishing average frictional force from its fluctuations showed the importance of system size. Uniform velocities were obtained in the sliding direction across each adsorbed film, with a gradient across the gap containing trapped fluid. The calculated friction stress was consistent with measurements reported using a surface forces apparatus, indicating that drag between an adsorbed layer and trapped fluid can account sufficiently for sliding friction in friction modifier systems. An example is shown in which changes in molecular organisation parallel to the surface led to a large change in normal force but no change in friction force.     

A molecular dynamics study of the structural change differences between Au225 and Au369 clusters on MgO surfaces at low temperature

The differences in structural change between Au₂₂₅ and Au₃₆₉ clusters with their (111) facets supported on MgO(100) surfaces at 5 K are studied by using molecular-dynamics simulations with the atomic interchange potentials of the Au/MgO interface. The parameters are obtained from the ab initio energies using the Chen-Möbius inversion method. Analyses of the pair distribution functions show that the two Au clusters use different deformation processes to adjust the distances between the interface atoms, owing to the misfit between the atom distances among the clusters and the substrates. The local structural changes are identified by atomic density profiles.     

载能氢原子与石墨(001)面碰撞过程中的能量传递行为的分子动力学研究

本文采用分子动力学方法研究了单一载能氢原子与石墨碰撞时氢原子被石墨反射、 吸附和石墨被氢原子穿透的发生系数以及碰撞中的能量传递机理. 研究发现: 与单层石墨相比, 多层石墨之间的长程相互作用增加了氢原子发生反射的能量范围, 尤其当入射能量大于20.0 eV时, 对反射过程的影响很明显; 当氢原子的入射能量大于25.0 eV时, 有一定的概率穿透四层石墨; 当氢原子入射能量高于28.0 eV时, 载能氢原子的能量传递给第二层石墨烯的比传递给第一层石墨烯的多. 这些结果对理解聚变反应中, 碳基材料的化学腐蚀及氚滞留有重要意义.     

棱形石墨烯纳米孔道内离子传输性能的分子动力学模拟

纳米通道的尺寸、结构和表面化学对其内部溶液的分布结构和输运性质有着重大影响.本文研究了一种全新的菱形石墨烯纳米通道.这种理想的通道与最近被广泛研究的金属有机框架材料(MOF)的内部结构类似,有着与传统的碳纳米管截然不同的内部结构.本文使用分子动力学模拟的方法研究在不同尺寸的菱形石墨烯纳米通道内的KCl溶液的性质,并将其与同尺寸的单壁碳纳米管进行了比较.研究结果表明在小孔道内(<20?)其内部的溶液结构呈现若干个密度极高的聚集区域,即出现了结晶化的趋势．这一研究结果,将为MOF的结构设计提供思路,从而有望实现类似于生物离子通道的高选择性和高传输能力.     

Structure and phase transition of a two—dimensional dusty plasma

The structure and phase transition of a two-dimensional (2D) dusty plasma have been investigated in detail bymolecular dynamics simulation. Pair correlation function, static structure factor, mean square displacement, and bondangle correlation function have been calculated to characterize the structural properties. The variation of internalenergy, shear modulus, particle trajectories and structural properties with temperature has been monitored to studythe phase transition of the 2D dusty plasma system. The simulation results are in favour of a two-step continuoustransition for this kind of plasma.     

单空位缺陷对载能氢原子与石墨层间碰撞的能量交换的影响的分子动力学研究

本文采用分子动力学方法研究空位缺陷对石墨层中碳氢粒子碰撞的影响.将氢原子以不同能量分别向单空位缺陷边缘的两个碳原子轰击,分析了入射氢原子的能量损失、发生吸附反应的能量范围和靶原子的能量传递过程.研究发现,单空位缺陷边缘的碳氢粒子更易发生吸附反应;在碳氢粒子正碰过程中,氢原子随入射能量变化出现了双反射区域;碳氢粒子在空位缺陷边缘吸附后,形成了高结合能的sp²结构,并出现悬挂键,其临近的碳碳键能未降低;单空位缺陷边缘的碳原子吸附氢原子能量的能力强而传递能量的能力弱.这些结果对理解聚变反应 关键词： 面向等离子体材料 分子动力学方法 单空位缺陷     

Tethered polyelectrolytes under the action of an electrical field

For a polyelectrolyte undergoing electrophoretic motion, it is predicted (D. Long, J.L. Viovy, A. Ajdari, Phys. Rev. Lett. 76, 3858 (1996); D. Long, A. Ajdari, Electrophoresis 17, 1161 (1996)) that the mechanical force necessary to stall the molecule is substantially smaller than the sum of electrical forces applied on all monomers. In fact, it should be proportional to its hydrodynamic friction coefficient and therefore to the size of its conformation. In our work we examine this prediction using coarse-grained molecular-dynamics simulations in which we explicitly include the polymer, the solvent, the counterions and salt. The electrophoretic mobility of polyelectrolytes is evaluated, the mechanical force necessary to keep the molecules tethered is measured and the resulting anisotropic polymer conformations are observed and quantified. Our results corroborate Long et al.'s prediction.     

小尺寸铜团簇冷却与并合过程中结构变化的原子尺度研究

研究Cu_N（N=57,58,59）熔融铜团簇在冷却过程以及300 K时两个具有二十面体结构Cu₅₅团簇在并合过程中的结构变化.对这些小尺寸团簇的结构变化采用基于嵌入原子方法的正则系综分子动力学进行计算机模拟.通过对模拟结果的分析表明,小团簇的冷却和并合过程存在阶段变化的特点.降温过程中Cu_N（N=57,58,59）团簇的原子运动及其微观结构变化表现出较大差异,由此导致这三类团簇内原子排布的不同,其中Cu₅₉团簇结构的有序程度最低.在两个Cu₅₅团簇并合早期阶段,这两个团簇相接触后发生变形导致原子位置出现较大改变,在随后的并合过程中,原子扩散引起原子局部位置调整导致所并合体系的结构发生变化.远离两个团簇接触区的原子仍保持其并合前的结构. 关键词： 团簇 分子动力学 计算机模拟 表面     

Si表面间水平碳纳米管束的分子动力学模拟研究

本文采用分子动力学模拟方法研究了Si表面间单壁水平碳纳米管束SWCNT (10,10)的变形和摩擦特性.系统在弛豫平衡后,首先对碳纳米管束施加压力至碳纳米管或Si表面结构破坏.之后在无压力和高压力两种情况下使上表面沿水平方向做剪切运动以研究碳纳米管束的摩擦特性.结果表明,由于碳纳米管的柔韧性,碳纳米管束在加载过程中出现明显变形,但直至3.8 GPa高压下并无结构破坏.系统无压力时SWCNT (10,10)在原地轻微随机滚动,压力为3.8 GPa时,碳纳米管束出现了整体的轻微滑动,同时伴随无规律的轻微滚动, 关键词： 碳纳米管束 摩擦 分子动力学模拟     

二元Lennard-Jones液体的相分离过程及其扩散性质的分子动力学研究

采用分子动力学(MD)模拟方法,研究了二元体系中相分离过程、粒子的扩散系数以及相分离域尺寸大小随温度的变化规律.发现,相分离域随温度的生长过程可以分为两个阶段,分别是温度比较高的快速生长阶段和低温时的稳定生长阶段;相分离体系中系统的扩散激活能不是常数,而是一个随温度变化的函数,并且当温度高于60 K时,满足关系式E(T)=a+bT^c.讨论了组元尺寸的变化对相分离过程的影响.结果表明,随两组元中某一组元 关键词： 相分离 扩散 分子动力学模拟     

载能氢同位素原子与石墨(001)面碰撞的分子动力学研究

采用分子动力学方法研究了载能H同位素原子与石墨晶体碰撞的同位素效应. 碳氢系统的强共价键作用和石墨层间的弱van der Waals力分别用REBO和Ito半经验势函数来描述. 研究发现: 随着入射原子质量的增加, 上表面吸附几率和反射几率的峰值都会向高能区移动; 相比于H, ²H入射原子, ³H入射原子具有较高的吸附几率——包括上表面吸附和内部吸附; 穿透石墨晶体, ²H, ³H原子所需的能量较高; 原子质量和原子入射能量都会影响入射粒子与不同石墨层之间的能量传递过程. 这些结果对理解碳基材料的³H滞留机制有重要意义.     

Molecular dynamical simulation of the canonical ensemble

We apply a technique to simulate the canonical ensemble, mixing molecular dynamics and Monte Carlo techniques, in which particles suffer virtual hard shocks. In the limit of infinite time the system approaches a Boltzmann distribution. A good approximation to the Boltzmann distribution is achieved in computationally accessible time for some model systems including the one-dimensional jellium.     

应力诱发NiAl单晶马氏体相变的分子动力学模拟

利用嵌入原子势（EAM）,对NiAl单晶在外应力作用下的动态拉伸过程进行了分子动力学模拟.应力-应变曲线分析以及原子构型分析表明外应力诱发NiAl合金发生了马氏体相变,原子结构由B2相转变为L1₀相.通过研究原子构型的演化过程,发现马氏体相变是通过多个{110}孪晶面的扩展和湮灭作用来完成的.同时探讨了马氏体相变的微观机理. 关键词： 马氏体相变 NiAl 分子动力学模拟     

Molecular dynamics simulation of the melting of a twist Σ=5 grain boundary

The existence of a possible grain boundary disordering transition of the melting type in a =5 (001) twist boundary of aluminium bicrystal below the melting temperature was investigated using a constant pressure molecular dynamics simulation. The calculated melting temperature T _cm of the bulk Al is about 960 K. The total internal energy, the structure factor, and the pair distribution function were calculated at different layers across the grain boundary. The mean atomic volume, the grain boundary energy, and the thermal expansion coefficients were also calculated using the same simulation method. This simulation also allows us to image the grain boundary structure at different temperatures. The equilibrium grain boundary structure at 300 K retains the periodicity of the coincident site lattice, so that the lowest energy structure corresponds to the coincident site arrangement of the two ideal crystals. With increasing temperature, the total internal energy of the atoms for both the perfect crystal and the grain boundary increases, as do the number of layers in the grain boundary. The grain boundary core exists and the perfect crystal structure still exists outside the grain boundary at 0.9375 T _cm. However, two atomic layers of the equilibrium grain boundary structure at 0.9375 T _cm lose the coincident site lattice periodicity and attain a structure with liquid-like disorder. Therefore, partial melting of the grain boundary has occurred at the temperature above 0.9375 T _cm which is in agreement with the experimental results.