液态铝基合金微观结构及其性质的分子动力学模拟研究

利用分子动力学研究了液态铝基合金Al80Fe20和Al80Cu20在 800~2200K温度范围内的微观结构和性质的变化。扩散系数的计算结果显示两合金有着不同的动力学行为。温度大约1400K以下液态Al80Cu20的扩散系数随温度的变化比Arrhenius关系所预期的要慢些,所得结果与Brillo研究小组的实验结果基本一致。当液态铝基合金的动力学发生异常变化时,其微观二十面体短程序的含量和作用也发生相应的变化。因此,熔体的奇特动力学行为可能由于其微观局域结构短程序随着温度的变化而异常地改变而引起的。     

液态铝基合金微观结构及其性质的分子动力学模拟研究(英文)

利用分子动力学研究了液态铝基合金Al80Fe20和Al80Cu20在800~2200K温度范围内的微观结构和性质的变化.扩散系数的计算结果显示两合金有着不同的动力学行为.温度大约1400K以下液态Al80Cu20的扩散系数随温度的变化比Arrhenius关系所预期的要慢些,所得结果与Brillo研究小组的实验结果基本一致.当液态铝基合金的动力学发生异常变化时,其微观二十面体短程序的含量和作用也发生相应的变化.因此,熔体的奇特动力学行为可能由于其微观局域结构短程序随着温度的变化而异常地改变而引起的.     

Cu-Ni合金结构的分子动力学模拟

采用EAM作用势对Cu Ni合金的结构特性进行了MD模拟研究 .通过FZ结构因子可发现 ,Cu含量的变化对结构因子的波动影响很小 ,键取向序参数和键对也表现出相似的变化规律 ,这表明液态Cu Ni合金对成份变化不敏感 ,体系中的化学序较弱 .将Cu70 Ni3 0 合金熔体的FZ结构因子与Waseda的实验结果进行对比 ,发现二者吻合得较好 ,表明EAM势可以很好地描绘Cu Ni合金的结构特性 .在快速冷却过程中 ,除了Cu2 0 Ni80 合金外 ,其他合金成份的双体分布函数的第二峰都发生了劈裂 ,标志着体系最终形成了非晶结构 ,而Cu2 0 Ni80 合金的双体分布函数却表现出晶体峰的特征 .通过对键取向序参数、键型指数以及铜镍原子的有效扩散系数的分析表明 ,在快速冷却过程中 ,Cu2 0 Ni80 合金最终形成了hcp晶体结构     

分子动力学模拟组分对Co_(1415-x)Al_x团簇结构特征的影响

本文采用分子动力学模拟方法结合镶嵌原子势,研究了在200 K时二元(Co Al)1415团簇的结构随Co原子浓度的变化情况.利用径向分布函数、对分析技术和键取向序参数方法研究了微观局部结构情况,研究结果表明:(Co Al)1415团簇的组分对最终冷却结构影响较大,Co原子浓度为100%~70%的团簇表现出不完全的六角密排结构特征;Co原子浓度为50%的团簇具有局部的体心立方体结构特征;Co原子浓度为30%~10%时,表现出部分区域的二十面体和缺陷二十面体结构特征.     

分子动力学模拟组分对Co1415-xAlx团簇结构特征的影响

本文采用分子动力学模拟方法结合镶嵌原子势,研究了在200 K时二元（CoAl）1415团簇的结构随Co原子浓度的变化情况。利用径向分布函数、对分析技术和键取向序参数方法研究了微观局部结构情况,研究结果表明： （CoAl）1415团簇的组分对最终冷却结构影响较大,Co原子浓度为100%~70%的团簇表现出不完全的六角密排结构特征;Co原子浓度为50%的团簇具有局部的体心立方体结构特征;Co原子浓度为30%~10%时,表现出部分区域的二十面体和缺陷二十面体结构特征。     

First-principles molecular dynamics simulation of biased electrode/solution interface

First-principles molecular dynamics simulations have been carried out for water in contact with Pt(1 1 1) surface. To apply negative bias potential to the water/Pt interface, excess electrons were added to our slab model using the recently developed computational scheme called “effective screening medium (ESM)”. Water molecules located away from the surface reoriented themselves to screen the electric field, but they responded differently near the surface. Water molecules nearest to the surface, forming a distinct layered structure with the hydrogen atom directed to the surface, increased the density with increasing field. On these bases, we discuss microscopic aspects of the electric double layer.     

Torsional buckling behavior of boron-nitride nanotubes using molecular dynamics simulations

Based on molecular dynamics simulations, the mechanical properties and buckling behavior of boron-nitride nanotubes under the action of torsional load are investigated. According to the results, the torsional properties of a boron-nitride nanotube are higher than those of its carbon counterpart. The effect of geometrical parameters on these parameters is also investigated. It is observed that by increasing the radius of nanotubes of the same length, unlike the critical shear strain, the critical torque considerably increases. The effect of chirality is also found to be negligible in the cases of critical shear strain and buckling mode, unlike the critical torque.     

Monte Carlo molecular dynamics simulations for two-dimensional magnets

A combined Monte Carlo molecular dynamics simulation technique is used to study thedynamic structure factor on a square lattice for isotropic Heisenberg and planar classical ferromagnetic spin Hamiltonians.     

Determination of the activation energy by stochastic analyses of molecular dynamics simulations of dislocation processes

An investigation is reported of the probability and the probability density of thermal activation of stress-driven dislocation processes, as simulated using molecular dynamics (MD). Stochastic analyses of the survival probability are found to lead to simple relationships between the loading history and the distribution of the interaction time and strength. It is shown that the determination of the activation energy associated to a thermally activated event can be achieved by a reduction of the stochastic process to a process obeying the Poisson's distribution, preserving the activation probability at the survival time. The method is applied to the kink-pair mechanism for screw dislocations in iron. Predictions are compared with experimental results and with other methods reported in the literature, which allows the difference in the approximations and in the assumptions considered in these models to be underlined.     

Properties of MgO at high pressures: Shell-model molecular dynamics simulation

Shell-model molecular dynamics (MD) simulation has been performed to investigate the melting of the major Earth-forming mineral: periclase (MgO), at elevated temperatures and high pressures, based on the thermal instability analysis. The interatomic potential is taken to be the sum of pair-wise additive Coulomb, van der Waals attraction, and repulsive interactions. The MD simulation with selected Lewis–Catlow (LC) potential parameters is found to be very successful in describing the melting behavior for MgO, by taking account of the overheating of a crystalline solid at ambient pressure. The thermodynamic melting curve is estimated on the basis of the thermal instability MD simulations and compared with the available experimental data and other theoretical results in the pressure ranges 0–150 GPa. Our simulated melting curve of MgO is consistent with results obtained from Lindemann melting equation and two-phase simulated data at constant pressure by Belonoshko and Dubrovinsky, in the pressure below 20 GPa. The extrapolated melting temperatures in the lower mantle are in good agreement with the results obtained from Wang's empirical model up to 100 GPa. Compared with experimental measurements, our results are substantially higher than that determined by Zerr and Boehler, and the discrepancy between DAC and MD melting temperatures may be well explained with different melting mechanisms. Meanwhile, the radial distribution functions (RDFs) of Mg–Mg, O–Mg, and O–O ion pairs near the melting temperature have been investigated.     

A modified equation of state for Xe at high pressures by molecular dynamics simulation

The exact equation of state （EOS） for the fission gas Xe is necessary for the accurate prediction of the fission gas behavior in uranium dioxide nuclear fuel, However, the comparison with the experimental data indicates that the applicable pressure ranges of existing EOS for Xe published in the literature cannot cover the overpressure of the rim fission gas bubble at the typical UO2 fuel pellet rim structure. Based on the interatomic potential of Xe, the pressure-volume-temperature data are calculated by the molecular dynamics （MD） simulation. The results indicate that the data of MD simulation with Ross and McMahan＇s potential [M. Ross and A. K. McMahan 1980 Phys. Rev. B 21 1658] are in good agreement with the experimental data. A preferable EOS for Xe is proposed based on the MD simulation. The comparison with the MD simulation data shows that the proposed EOS can be applied at pressures up to 550 MPa and 3 GPa and temperatures 900 K and 1373 K respectively. The applicable pressure range of this EOS is wider than those of the other existing EOS for Xe published in the literature.     

Medium-range structural order in liquid Ni20Al80 alloy: Experimental and molecular dynamics studies

In this Letter, structural properties of liquid Ni₂₀Al₈₀ have been studied through molecular dynamics (MD) simulations, and results agree well with X-ray experiments. The prepeak in the structure factor manifests the existence of medium-range structural order (MRSO) in liquid Ni₂₀Al₈₀. Analyses in MRSO reveal the dominance of ideal icosahedra and defective icosahedra which, furthermore, form clusters and share common atoms.     

钴和铁熔体短程序的分子动力学模拟

通过分子动力学模拟研究了金属钴和铁熔体从普通液态到过冷液态普通微观结构的性质.所计算两体分布函数与可获得的实验结果基本一致,从普通液态到过冷液态的局域结构中的原子配位数为11和12的分布几率随着温度的降低而变小,而原子的配位数大于12的分布几率反而增加,角分布函数位于55.有一个明显的峰,位于110°有一个宽展的峰,位于150°有一个肩膀,表明金属钴和铁熔体的微观局域结构要比规则的二十面体团簇的结构复杂得多.我们的模拟结果表明金属熔体中的二十面体短程序随着温度的降低而增加,并在过冷液态中而占优势.     

氦热力学性质的分子动力学研究

 利用经典分子动力学和第一性原理分子动力学,研究了氦在高压下的熔化曲线、状态方程和非金属-金属转变。得到了氦在温度小于4.5 eV、 密度0.3~5.0 g/cm³范围内的状态方程,并把氦的熔化曲线的压强范围拓展到了50 GPa。氦的能隙宽度曲线表明,温度大大降低了氦的金属化密度。     

Investigation of the structural and dynamic basis of kinesin dissociation from microtubule by atomistic molecular dynamics simulations

Kinesin is a molecular motor that can step processively on microtubules via the hydrolysis of ATP molecules. An important factor characterizing the processivity of the kinesin motor is its dissociation from the microtubule. Here, using all-atom molecular dynamics simulations, we studied the dissociation process of the kinesin head in weak-microtubule-binding or ADP state from tubulin on the basis of the available high-resolution structural data for the head and tubulin. By analyzing the simulated snapshots of the structure of the head-tubulin complex we provided detailed structural and dynamic information for the dissociation process. We found that the dissociation of the head along different directions relative to the tubulin exhibits very different dynamic behaviors. Moreover, the potential forms or energy landscapes of the interaction between the head and tubulin along different directions were determined. The studies have important implications for the detailed molecular mechanism of the dissociation of the kinesin motor and thus are critical to the mechanism of its processivity.     

Energy flow in thermostatted non-equilibrium molecular dynamics simulations

Energy transfers between internal kinetic and potential energy reservoirs in a simple liquid are studied by setting the temperature of one energy reservoir to a different value from that of the others and computing the resultant energy flows. In the first set of simulations, the x-directional kinetic temperature was artificially raised above the other five, and in the second, the x-directional configurational temperature was artificially raised above the other five. In both cases, external energy flows balanced, but unexpected energy flows between different directional components of the potential energy were observed. Additional simulations showed that these energy flows occurred regardless of the arrangement of thermostats imposed on the six degrees of freedom and the addition of shear. Heat flow between degrees of freedom that were ostensibly at the same temperature was anomalously observed. It was concluded that a different breakdown of the contributions to the configurational energy that is consistent with the definition of the directional configurational temperatures is required.     

含氢碳膜的生长机制: 分子动力学模拟研究低能量CH基团的作用

探索等离子增强化学气相沉积(PECVD)技术中含氢碳膜的生长机理, 制备出常态超润滑含氢碳膜是表面工程技术领域的目标之一. 基于REBO势函数, 采用分子动力学模拟方法, 通过对比研究CH基团在清洁金刚石和吸氢金刚石表面的沉积行为, 发现低能量CH基团在清洁金刚石(111)面上的吸附效率大于98%, 而在吸氢金刚石(111)面上的吸附效率低于1%. 结果表明PECVD法制备含氢碳膜时, 低能量CH基团对薄膜生长的贡献主要来自于其在表面非饱和C位置的选择性吸附.     

First-principles study of lattice dynamics and thermodynamics ofosmium under pressure

We have performed the first-principles linear response calculations of the lattice dynamics, thermal equation of state and thermodynamical properties of hcp Os metal by using the plane-wave pseudopotential method. The thermodynamical properties are deduced from the calculated Helmholtz free energy by taking into account the electronic contribution and lattice vibrational contribution. The phonon frequencies at Gamma point are consistent with experimental values and the dispersion curves at various pressures have been determined. The calculated volume, bulk modulus and their pressure derivatives as a function of temperature are in excellent agreement with the experimental results. The calculated specific heat indicates that the electronic contribution is important not only at very low temperatures but also at high temperatures due to the electronic thermal excitation. The calculated Debye temperature at a very low temperature is in good agreement with experimental values and drops to a constant until 100~K.     

氦热力学性质的分子动力学研究

利用经典分子动力学和第一性原理分子动力学,研究了氦在高压下的熔化曲线、状态方程和非金属-金属转变。得到了氦在温度小于4.5 eV、 密度0.3~5.0 g/cm3范围内的状态方程,并把氦的熔化曲线的压强范围拓展到了50 GPa。氦的能隙宽度曲线表明,温度大大降低了氦的金属化密度。     

Numerical distortion and effects of thermostat in molecular dynamics simulations of single-walled carbon nanotubes

In this paper, single-walled carbon nanotubes （SWCNTs） are studied through molecular dynamics （MD） simulation. The simulations are performed at temperatures of 1 and 300K separately, with atomic interactions characterized by the second Reactive Empirical Bond Order （REBO） potential, and temperature controlled by a certain thermostat, i.e. by separately using the velocity scaling, the Berendsen scheme, the Nose-Hoover scheme, and the generalized Langevin scheme. Results for a （5,5） SWCNT with a length of 24.5 nm show apparent distortions in nanotube configuration, which can further enter into periodic vibrations, except in simulations using the generalized Langevin thermostat, which is ascribed to periodic boundary conditions used in simulation. The periodic boundary conditions may implicitly be applied in the form of an inconsistent constraint along the axis of the nanotube. The combination of the inconsistent constraint with the cumulative errors in calculation causes the distortions of nanotubes. When the generalized Langevin thermostat is applied, inconsistently distributed errors are dispersed by the random forces, and so the distortions and vibrations disappear. This speculation is confirmed by simulation in the case without periodic boundary conditions, where no apparent distortion and vibration occur. It is also revealed that numerically induced distortions and vibrations occur only in simulation of nanotubes with a small diameter and a large length-to-diameter ratio. When MD simulation is applied to a system with a particular geometry, attention should be paid to avoiding the numerical distortion and the result infidelity.