纳米铝热剂Al/SiO_2层状结构铝热反应的分子动力学模拟

利用分子动力学模拟方法和反应力场势函数研究了Al/SiO_2层状纳米体系的铝热反应,模拟了在不同初始温度下(600,700,800,900,1000和1100 K)绝热反应的结构变化和能量性质,发现Al/SiO_2体系的铝热反应是自加热的氧化还原反应,当初始温度为900和1000 K时,Al经历了熔化前的一个临界状态,与SiO_2的铝热反应比较活跃,系统温度随着反应时间的增加不断升高,当初始温度为600,700,800和1100 K时,初始温度越高,在Al和SiO_2界面形成的Al-O层越薄,系统发生铝热反应达到的最终绝热温度越高,所用的时间(有效反应时间τ)越短,即界面扩散阻挡层的厚度对铝热反应的自加热速率产生了影响,初始温度为600,700,800,1100K时的自加热速率分别为3.4,3.5,4.7和5.4K/p8.A1/SiO_2体系的铝热反应析出了Si单质,与实验结果相符合。     

Dissipative particle dynamics study of velocity autocorrelation function and self-diffusion coefficient in terms of interaction potential strength

This research focuses on numerically investigating the self-diffusion coefficient and velocity autocorrelation function (VACF) of a dissipative particle dynamics (DPD) fluid as a function of the conservative interaction strength. Analytic solutions to VACF and self-diffusion coefficients in DPD were obtained by many researchers in some restricted cases including ideal gases, without the account of conservative force. As departure from the ideal gas conditions are accentuated with increasing the relative proportion of conservative force, it is anticipated that the VACF should gradually deviate from its normally expected exponentially decay. This trend is confirmed through numerical simulations and an expression in terms of the conservative force parameter, density and temperature is proposed for the self-diffusion coefficient. As it concerned the VACF, the equivalent Langevin equation describing Brownian motion of particles with a harmonic potential is adapted to the problem and reveals an exponentially decaying oscillatory pattern influenced by the conservative force parameter, dissipative parameter and temperature. Although the proposed model for obtaining the self-diffusion coefficient with consideration of the conservative force could not be verified due to computational complexities, nonetheless the Arrhenius dependency of the self-diffusion coefficient to temperature and pressure permits to certify our model over a definite range of DPD parameters.     

Vapour-liquid equilibria of the two- and three-dimensional monoatomic classical fluids interacting via double Yukawa potential

We have carried out Monte Carlo simulations in Gibbs ensemble for two-and three-dimensional double Yukawa fluid. We have compared liquid-vapour equilibrium curve with that of Lennard-Jones, when parameters occurring in double Yukawa potential are chosen to fit Lennard-Jones potential. The results are in good agreement. The role of repulsive and attractive contributions for the potential on the liquid-vapour coexistence region as well as on critical temperature and critical density has been studied. The critical temperature is found to be more sensitive than the critical density to the variation in repulsive and attractive parts of the potential. Also, the range of the attractive interaction directly influences range of the liquid vapour coexistence region. It has been found that smaller the values of the attractive parameter, larger is the coexistence region.     

Diffusion and thermite reaction process of film-honeycomb Al/NiO nanothermite:Molecular dynamics simulations using ReaxFF reactive force field

The diffusion and thermite reaction process of Al/NiO nanothermite composed of Al nanofilm and NiO nano honeycomb are investigated by molecular dynamics simulations in combination with the Reax FF. The diffusion and thermite reaction are characterized by measuring energy release, adiabatic reaction temperature, and activation energy. Based on time evolution of atomic configuration and mean square displacement, the initialization of the thermite reaction process of Al/NiO nanothermite results from the diffusion of Al atoms. Under the microcanonical ensemble, it is found that the adiabatic reaction temperature of the thermite reaction process of Al/NiO nanothermite reaches over 5500 K, and activation energy is 8.43 k J/mol. The release energy of the thermite reaction process of Al/NiO nanothermite is 2.2 k J/g, which is in accordance with the available experimental value. With the same initial temperature, the adiabatic reaction temperature of the thermite reaction process of Al/NiO nanothermite has a tendency to decrease dramatically as the equivalence ratio increases. On the basis of chemical bond analysis, the initial temperature and equivalence ratio have great effects on the thermite reaction process, but do not significantly affect the average length of Al–Ni nor Al–O bond. Overall, the thermite reaction of film-honeycomb Al/NiO nanothermite is a complicated process instead of a theoretical equation.     

用FS多体势模型模拟金属铜的冷却过程

采用常温、常压分子动力学模拟方法模拟了在周期性边界条件下由500个原子构成的液态Cu模型系统的凝固过程,考察了不同降温速率下液态Cu的凝固行为.模拟结果很好地重现了实验值.模拟中原子间作用势采用FS势,结构分析采用双体分布函数、对分析技术、内能、均方位移等方法,得到了原子体系微观结构组态变化的重要信息,并利用能量分析的方法对体系微观结构的变化进行了说明.     

分子动力学模拟压水反应堆中氢气对水的影响

通过分子动力学方法模拟了在常温常压下(1 atm, 298 K)和在压水堆环境下(155 atm, 626 K), 水分子数为256, 氢分子数为0, 25, 50, 75和100等不同数目时, 粒子系统的动力学性质和微观结构, 分析了不同氢气对水中溶解氧的影响. 从模拟结果可知, 在常温常压和压水堆环境下, 当氢粒子数分别为0, 25, 50, 75和100时, 粒子系统的均方位移会随氢分子数增加而增加, 并且常温常压下的增长幅度远小于压水堆环境下的增长幅度, 如压水堆环境下氢分子数为75时系统的均方位移约是常温常压下氢分子数为75时系统的均方位移的6.02倍, 比压水堆环境下氢分子数0时系统的均方位移增加了131.88%. 此外, 粒子系统的微观结构, 从径向分布函数看, 在常温常压下随着氢分子数目的增加而小幅度增加, 这与常温常压下因氢气溶解在水中增大了氧离子周围的粒子密度相符合. 而在压水堆环境下, 氢分子数为75, 50, 25与为0时的水比较, 其径向分布均不会有太大的变化, 而分子数为100时会出现明显增加, 与为0时的水比较其径向分布增加了22.00%. 模拟结果表明, 往压水堆中的水加入氢气能明显地抑制水中的溶解氧.     

分子动力学方法研究金属Ti中He小团簇的迁移

采用分子动力学方法模拟了不同温度下He原子及He团簇在金属Ti中的迁移特性,并计算了扩散前系数和激活能. 研究发现这种扩散的各向异性非常显著,具体表现在He原子及He团簇在不同方向上扩散系数的前系数完全不同,但它们的激活能却相同. 研究表明:在预测金属中He的扩散行为时,必须采用动态模拟方法才能得到准确的前系数,仅仅考虑势垒的静态模拟方法是不行的. 另外,还发现一个不同于直觉的现象,即较低温度下He二聚物的扩散系数比单个He原子的扩散系数大;此外,在所模拟的温度范围内Arrhenius方程能够很好地描述它们的扩散. 这说明动力学模拟对预测金属中He的扩散行为具有重要的意义. 关键词： 分子动力学 扩散系数 各向异性 氦     

Self-diffusion of particles interacting through a square-well or square-shoulder potential

The diffusion coefficient and velocity autocorrelation function for a fluid of particles interacting through a square-well or square-shoulder potential are calculated from a kinetic theory similar to the Davis-Rice-Sengers theory and the results are compared to those of computer simulations. At low densities the theory yields too low estimates due to the neglect of correlations between subsequent partial collisions of identical pairs; in particular, the neglect of boundstate effects appears important. At intermediate densities the theory makes reasonable predictions and at high densities it produces too high values, due to the neglect of ring terms and other correlated collision events. The results for the square-shoulder potential generally exhibit better agreement between theory and simulations than do those for the square-well potential.     

Diffusion in a periodic potential with a local perturbation

We consider the diffusion of a particle at X_t in a drift field derived from a smooth potential of the formV+B, whereV is periodic andB is a bump of compact support. With no bump,B=0, the mean squared displacementE(t) E |X _t – X₀|² =D(V)t +C +O(e ^–t ),>0, in any dimension. WhenB0, we establish in one dimension the asymptotic expansion , 0, ast. Our analysis relies on the Nash estimates developed in previous work for the transition density of the process and their consequences for the analytic structure,of the Laplace transform ofE(t).     

Molecular dynamics results for the radial distribution functions of highly asymmetric hard sphere mixtures

Molecular dynamics (MD) results for the radial distribution functions of mixtures of large and small hard spheres are reported for size ratios whose (large/small) values are 1, 2.5, 5, 7.5, and 10 in the region where the concentration of the large spheres is very small. The MD contact values of these functions are compared with formulae due to Boublik, Mansoori, Carnahan, Starling, Leland, Grundke, and Henderson, Viduna and Smith, Henderson, Trokhymchuk, Woodcock, and Chan, as well as new formulae that are considered here. The new formulae give good agreement for the large–small contact values and reasonably good agreement for the large–large contact values. The Viduna–Smith formula is satisfactory for the small–small contact value and quite reasonable for the small–large contact value. Undoubtedly, further improvements are possible. These results give insight into what may be called the colloidal limit, where the size ratio is exceedingly large while the concentration of the large spheres is exceedingly small, and into the passage to this limit.     

Molecular dynamics calculations of shear viscosity time-correlation functions for hard spheres

The time-correlation function for shear viscosity is evaluated for hard spheres at volumes of 1.6 and 3 times the close-packed volume by a Monte Carlomolecular dynamics technique. At both densities, the kinetic part of the timecorrelation function is consistent, within its rather large statistical uncertainty, with the long-timet ^–3/2 tail predicted by the mode-coupling theory. However, at the higher density, the time-correlation function is dominated by the cross and potential terms out to 25 mean free times, whereas the mode-coupling theory predicts that these are asymptotically negligible compared to the kinetic part. The total time-correlation function decays roughly ast ^–3/2, with much larger than the mode-coupling value, similar to the recent observations by Evans in his nonequilibrium simulations of argon and methane. The exact value of the exponent is, however, not very precisely determined. By analogy with the case of the velocity autocorrelation function, for which results are also presented at these densities, it is argued that it is quite possible that at high density the asymptotic behavior is not established until times substantially longer than those attainable in the present work. At the lower density, the cross and potential terms are of the same magnitude as the kinetic part, and all are consistent with the mode-coupling predictions within the relatively large statistical uncertainties.Work performed under the auspices of U.S. Department of Energy.     

N-丁基吡啶四氟硼酸盐/水二元体系的分子动力学模拟研究

利用分子动力学模拟方法研究了离子液体N-丁基吡啶四氟硼酸盐([BPy]BF₄)与Tip4p模型水分子的二元体系的微观结构.比较了各组分间的径向分布函数,结果表明随着离子液体比例的增加,水与阴离子、水与阳离子头部吡啶环、阳离子头部与阴离子上相关原子间的径向分布函数峰值都呈现递增的趋势,而阳离子上丁基链末端碳原子间的径向分布函数没有明显变化;空间分布函数则直观地反映出阴离子主要分布在阳离子的吡啶环周围,水分子在阴离子周围近似呈均等分布,且几率随离子液体比例增大而增加;另外还探讨了不同离子液体比例下的二元体系中氢键的数目和寿命,结果均呈现一定规律的变化.     

FexO-SiO2-CaO-MgO-“NiO”系镍渣势函数及分子动力学模拟

从非经验参数角度构建基于两原子模型的Fe_xO-SiO₂-CaO-MgO-"NiO"渣的势函数,利用分子动力学模拟揭示调控镍渣组分与结构及物理化学性能之间的关系。结果表明： BMH（Born-Mayer-Huggins）势函数能够较好表征Fe_xO-SiO₂-CaO-MgO-"NiO"系镍渣的势能。当镍渣中CaO含量为15 wt.%时,Si⁴⁺-Si⁴⁺间配位数最小,此时熔渣的聚合度最低,有利于扩散。当Fe元素由Fe²⁺转换为Fe³⁺后,和O^2-结合能力更强,会导致熔渣更难以扩散,因而镍渣的黏度迅速升高,造成冶炼条件恶化,因此在镍闪速熔炼时要严格控制Fe²⁺/Fe³⁺的比例。模拟计算的黏度与实测值吻合较好,表明构建的势函数能够较好地反映镍渣的物化性能。     

绝缘油在水和酸共同作用下的分子动力学模拟

为探讨水分和温度对变压器油中酸扩散行为的影响,在分子模拟软件内建立了水分子、矿物油分子、甲酸分子的复合模型,水的体积分数分别为1%,3%,5%,甲酸的质量分数为3%.通过分子动力学方法模拟计算了甲酸分子在不同含水量油模型中的扩散系数、径向分布函数以及复合模型的总动能;此外,选取了水体积分数为1%的复合模型为研究对象,以温度为变量进行动力学模拟.结果表明:随着水含量的增多,甲酸的扩散系数并没有逐渐增大,甲酸分子与水分子之间的氢键数目呈逐渐减少的趋势.而随着温度的升高,油模型具有的动能越大,酸分子的扩散逐渐增强.油模型中,含水量越多,水和酸在热运动下碰撞油分子的能量越大;温度越高,相互之间的碰撞越激烈.模拟结果说明水分增加或温度增高,均对油分子的结构稳定性造成一定程度的破坏.     

Molecular dynamics study of self-diffusion in Zr

We employed a recently developed semi-empirical Zr potential to determine the diffusivities in hcp and bcc Zr via molecular dynamics simulation. The point defect concentration was determined directly from molecular dynamics (MD) simulation rather than from theoretical methods using T = 0 calculations. Our MD simulation indicates that the diffusion proceeds via the interstitial mechanism in hcp Zr, and both vacancy and interstitial mechanisms contribute to diffusivity in bcc Zr. The agreement with the experimental data is excellent for hcp Zr and rather good for bcc Zr at high temperatures, but there is considerable disagreement at low temperatures.     

Size effect in the melting and freezing behaviors of Al/Ti core-shell nanoparticles using molecular dynamics simulations

The thermal stability of Ti@Al core/shell nanoparticles with different sizes and components during continuous heating and cooling processes is examined by a molecular dynamics simulation with embedded atom method. The thermodynamic properties and structure evolution during continuous heating and cooling processes are investigated through the characterization of the potential energy, specific heat distribution, and radial distribution function(RDF). Our study shows that, for fixed Ti core size, the melting temperature decreases with Al shell thickness, while the crystallizing temperature and glass formation temperature increase with Al shell thickness. Diverse melting mechanisms have been discovered for different Ti core sized with fixed Al shell thickness nanoparticles. The melting temperature increases with the Ti core radius. The trend agrees well with the theoretical phase diagram of bimetallic nanoparticles. In addition, the glass phase formation of Al–Ti nanoparticles for the fast cooling rate of 12 K/ps, and the crystal phase formation for the low cooling rate of 0.15 K/ps. The icosahedron structure is formed in the frozen 4366 Al–Ti atoms for the low cooling rate.     

导电石墨烯孔道内双电层结构的研究

为有效开发和利用新能源,人们迫切需要高性能的超级电容器提供能量的存储和转换.在超级电容器中双电层结构扮演着关键性的角色.本文利用分子动力学方法通过建立开放的石墨烯纳米孔道(1～2 nm),研究了KCl溶液在纳米孔道内的双电层结构,同时也比较了恒电量模拟(Q)和恒电势模拟法(U)下双电层结构的异同.结果表明在恒电势模拟法考虑了导电石墨烯壁的镜像作用使结果更符合实验中的材料系统.而石墨烯壁的镜像作用能额外吸附离子从而增强孔道内部的阴阳离子,这可能有助于电极电容的提升.通过对不同孔道高度的研究,本文发现水分子作为介电材料在水基超级电容器中发挥着决定性的作用.它能在很大程度上抵消不同离子和不同孔道高度下双电层的变化,从而在不同情况下获得了相似的电容.     

Origin of the log-periodic oscillations in the quantum dynamics of electrons in quasiperiodic systems

Recently, the occurrence of log-periodic oscillations in the quantum dynamics of electrons was reported for the one-dimensional Fibonacci quasicrystal by Lifshitz and Even-Dar Mandel. We apply a real-space renormalization group approach to show that these log-periodic oscillations are related to the underlying quasiperiodic structure of the Fibonacci quasicrystal. We find that they originate from the superposition of bonding and antibonding states associated to strongly coupled atoms in the chains, and that they show a hierarchical structure closely related to the atomic configurations.     

超临界Lennard-Jones流体结构特性分子动力学研究

研究超临界流体在不同压力和温度的结构特征有助于深刻理解并有效利用超临界流体.本文采用分子动力学方法模拟超临界压力、拟临界温度附近流体的结构及密度波动曲线的排列熵,分析状态参数变化的影响.结果表明,定压下,径向分布函数随温度升高,第一峰值位置逐渐向右移动,但右移幅度随着压力偏离临界点距离的增大而减弱,近临界压力时,出现峰值最高点的工况和等温压缩系数的极值点位置一致,压力增大,该现象消失.低压力拟临界点时易出现面积大、相对集中且分布稳定的高/低密度区,无明显嵌套现象.静态结构因子存在一定发散行为,发散的最大值和等温压缩系数极值点所处工况符合.低压力时密度时间序列的波动幅度最大,类周期现象较明显.在分子间势能、等温压缩系数和热运动效应的共同作用下,当压力(P)为1.1倍的临界压力(Pc)时,排列熵在0.99倍的拟临界温度(Tpc)达到最小值,P = 1.3Pc和1.5Pc时,最小排列熵与等温压缩系数的最大值工况点保持一致,压力继续增大,各模拟工况密度和排列熵的波动减弱,流体均匀性增强.     

Molecular dynamics study of structural evolution of aluminum during rapid quenching under different pressures

We present a molecular dynamics simulation for liquid Al during the rapid quenching under different pressures. The pair analysis technique and the probabilities of bond energy distribution of inherent structures have been employed to reveal the structural characteristics of liquid and glassy Al. During the liquid-glass quenching process, the bond pairs representing the degree of icosahedral short-range ordering are largely enhanced, whereas the bond pairs being related to fcc and hcp crystalline order increase at first then decrease. The pressure effect on various bond pairs for liquids is larger than for glasses. Two kinds of bond pairs, which exist in large proportion in the amorphous, are demonstrated for the transformation from 1431 to 1541 bond pairs when decreasing temperature or increasing pressure below glass transformation temperature (T_g). Although the sum of these two pairs keeps unchanged below T_g, the role of them is like a bridge which links another two different kinds of bonded pairs, icosahedral ordering and crystalline ordering pairs.