Mo高压熔化的分子动力学模拟

 采用第一原理方法计算了钼在零温下的结构,表明钼在500 GPa以下一直保持bcc结构（常温）,与实验一致。在零压附近计算了E-V关系,利用Murnaghan物态方程拟合得到了零压体积及其模量,与实验结果符合得很好。采用第一原理分子动力学模拟了钼的高压熔化性质。采用NVT系综计算了128个原子的系统,初始构形为bcc结构,体积分别为0.015 48、0.012 19、0.010 98、0.009 84、0.009 10 nm³/atom,计算了几个温度点,拟合得到了熔化曲线,熔化温度明显高于金刚石压砧（DAC）实验结果;将初始构形改变为fcc结构,模拟其熔化特性,得到的熔化温度明显下降,与激光加载DAC实验结果一致,认为可能的原因是钼熔化后形成的液体结构类似于fcc结构,而不是常态时的bcc结构。     

量子分子动力学模拟液体钚的输运性质

采用第一原理分子动力学（QMD）方法模拟液体钚的输运性质.计算的粘性和扩散系数在较低温度时与文献有明显差异,在实验测量范围内,模拟结果与实验一致,温度升高时数值模拟结果趋于一致.利用QMD的模拟结果计算了应力自相关函数和速度自相关函数,结果表明：在温度较低时,液体钚呈现明显的强关联特性.对于具有强关联特性的液体,利用较短时间的QMD模拟结果,通过简单e指数拟合外推到t→∞得到的扩散系数和粘性具有较大偏差,这是造成本文模拟结果与文献结果出现差异的主要原因.通过增加QMD模拟时间步数,获得了更为准确的输运性质.     

Isentropic Compression of Deuterium by Quantum Molecular Dynamics

In this work we present our calculations of the compression isentrope of deuterium. We use an ab initio molecular dynamics approach to compute pressure and internal energy of dense deuterium in the range of temperatures 293–25000 K and densities 0.9–4.3 g/cm³. We then restore the isentrope and compare the results with experimental data and other theories. The position of the calculated isentrope slightly depends on the initial point and agrees with the low‐pressure part and the highest‐pressure point of the measurements. However, we do not observe the density jump registered experimentally at P ≈ 130 GPa (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

纳米尺度下毛细流动的分子动力学模拟

采用分子动力学模拟研究纳米尺度下壁面润湿性对毛细流动的影响,主要考虑纳米通道两侧壁面润湿性相同与不同两种情况。研究表明：两侧壁面润湿性相同条件下,毛细流动随着壁面润湿性增强而加快, 毛细高度随时间的变化早期偏离Lucas-Washburn理论,但后期与其预测符合。在纳米通道两侧壁面润湿性不同的情况下,液面会发生振荡,两侧壁面毛细高度也不相等,且液面振荡的幅度和两侧壁面毛细高度差都随着两侧壁面润湿性差异的增大而增大。基于能量转化分析,提出两侧壁面湿润性不同情况下纳米通道中毛细流动发生的条件以及毛细流动快慢的判别依据。研究结果加深了对纳米尺度下毛细流动机理的认识,并为相关工程应用提供理论参考。     

液态氧化铝结构的分子动力学模拟

利用分子动力学模拟液态氧化铝从2663 K到过冷温度2223 K的结构变化,通过总双体分布函数得到的液态氧化铝的原子第一、第二近邻距离与Ansell研究小组的实验结果基本一致.通过配位数和偏角分布函数的分析,液态氧化铝的局域结构主要是由氧的三重铝原子配位和铝原子四重(AlO₄)和五重(AlO₅)氧原子配位组成.没有观察到可能因铝原子高、低配位数氧原子的团簇分离而出现的液-液相变.     

Molecular Dynamics Simulations of Warm Dense Carbon

We present classical and DFT‐based molecular dynamics (MD) simulations of carbon in the warm dense matter regime (? = 3.7 g/cc, 0.86 eV < T < 8.62 eV [T < 100 eV for classical MD]). Two different classical interatomic potentials are used: 1. LCBOP, designed to simulate condensed (e.g. solid) phases of C, and 2. linearly screened Coulomb (Yukawa) potentials. It is shown that LCBOP over‐predicts minima and maxima in the pair correlation functions of liquid‐C in this regime when compared to the DFT‐MD results. The screened Coulomb model, while under‐correlating at low‐T, seems to produce the correct qualitative features in the static ionic pair distributions at the highest‐T. However, both approaches predict the decay in the ionic contribution of the specific heat as T → ∞ to be much slower than that predicted by a model based on DFT‐MD. These differences in the MD‐derived equations of state in warm dense regimes could have important consequences when using classical inter‐ionic forces such as these in large‐scale MD simulations aimed at studying, for instance, processes of relevance to inertial confinement fusion when C is used as an ablator material. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Study of Friction between Liquid Crystals and Crystalline Surfaces by Molecular Dynamic Simulations

The lubrication characteristics of liquid crystal (LC) molecules sheared between two crystalline surfaces obtained from molecular dynamics (MD) simulations are reported in this article. We consider a coarse-grained rigid bead-necklace model of the LC molecules confined between two atomic surfaces subject to different shearing velocities. A systematic study shows that the slip length of LC lubrication changes significantly as a function of the LC-surface interaction energy, which can be well described though a theoretical curve. The slip length increases as shear rate increases at high LC-surface interaction energy. However, this trend can not be observed for low interaction energy. The orientation of the LC molecules near the surface is found to be guided by the atomics surfaces. The influence of temperature on the lubrication characteristics is also discussed in this article.     

液态合金NiAl凝固过程中微观结构转变的分子动力学模拟

采用分子动力学模拟方法对液态NiAl凝固过程进行了研究 ,考察了不同冷却速度下液态NiAl结构变化特点 ,原子间相互作用势采用F S多体势 ,结构分析采用键取向序和对分析技术 .计算结果表明 ,冷却速度对液态NiAl结构转变有重要影响 ,在不同的冷却速度下 ,NiAl凝固过程出现了明显不同 ,冷速为 4× 10 13 和4× 10 12 K/s时 ,NiAl快速凝固为无序的非晶体结构 ;而在较慢的 8× 10 11K/s冷速下 ,NiAl凝固为晶态结构 .给出了不同冷却速度下液态NiAl结构转变的微观信息 .     

高压下MgO熔化特性的分子动力学模拟

利用分子动力学方法,模拟研究了高压下MgO的熔化特性．通过晶体的现代熔化理论,对MgO的分子动力学模拟熔化温度进行了修正,得到了高温高压下MgO的熔化温度．计算得到的MgO熔化曲线和已有的实验及其它理论结果在0-135 GPa进行了比较,发现修正得到的MgO熔化温度和由Lindemann熔化方程及两相方法得到的结果在压力低于15 GPa时符合很好．同时,MgO熔化模拟有效解释了一阶相变分子动力学过程中出现的过热熔化现象．     

Shear Viscosity of Liquid‐Phase Yukawa Plasmas from Molecular Dynamics Simulations on Graphics Processing Units

We have carried out million‐particle equilibrium molecular dynamics simulations of 3‐dimensional Yukawa liquids in order to determine the shear viscosity coefficient. The computations have been executed on Graphics Processing Unit (GPU) architectures with our largely parallelized code. The results cover the strongly coupled liquid phase, with Γ up to the vicinity of the freezing transition, for the 1 ≤ κ ≤ 3 domain of the screening parameter of the Yukawa potential. The good agreement of the present results with those obtained from earlier simulations of significantly smaller systems (consisting of several hundred to several thousand particles) verifies that the viscosity data derived in these smaller scale simulations are also acceptable (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamics of the phase-change material GeTe across the structural phase transition

Results of inelastic neutron scattering experiments and ab initio molecular dynamics simulations for GeTe – the parent compound of phase-change materials are reported. The inelastic neutron spectra of GeTe powder samples have been determined in the temperature range extending from 300 to 700 K. The phonon peaks undergo thermal shifts resulting from anharmonic effects being weaker for acoustic than optic modes. A small concentration of free charge carries arising from the presence of Ge-vacancies was found not to affect significantly the neutron weighted phonon densities of states of GeTe. The spectral pattern changes qualitatively across the structural phase transition, but the local structure of GeTe remains hardly affected, as confirmed by the analysis of temperature dependence of the pairdistribution function obtained from ab initio molecular dynamics investigations. The present theoretical studies support in a wide extent our experimental observations and also those provided by local probe methods.     

超高压水的分子动力学模拟

 采用平衡分子动力学（EMD）方法,模拟研究了温度范围为243~348 K、压强范围为0.1~400 MPa条件下水的热力学性质、结构和动力学性质,模拟结果与实验值吻合较好。模拟结果表明,随着压强的增大,水分子间的氢键作用增强,扩散系数减小;随着温度的升高,水分子间的氢键作用减弱,有序程度下降,扩散系数增大。但在过冷水中,扩散系数随压强的增大有增加的趋势。     

纳米流体分散稳定性的分子动力学研究初探

纳米流体的制备和分散稳定性是纳米流体研究工作的前提。本文利用平衡分子动力学模拟方法,尝试从微观角度研究纳米流体的分散稳定性。以纳米颗粒Cu在液态CO_2中的分散稳定为算例,通过跟踪每个固体分子和液体分子的位置,分别观察了颗粒周围的液体吸附现象和颗粒的团聚现象。通过改变固液分子间的引力和斥力,观察其对颗粒周围液体分子的密...     

增强采样分子动力学模拟在生物分子研究中的进展

分子动力学模拟能够描述蛋白质分子在行使生物学功能过程中涉及的构象变化,已发展成为中物学研究中重要的计算工具.由于生物分子的构象分布存在崎岖的自由能面,在较为复杂的生物体系的模拟中,传统的分子动力学模拟的构象采样能力受到极大限制,模拟的时间尺度与真实的生物学过程之间仍存在差距.增强采样是解决这一问题的有效手段.本文综述了两类增强采样方法即约束型和无约束型增强采样算法的理论基础、最新进展及其在生物分子中的典型应用,同时也简要总结了组合型增强采样算法近些年的发展.     

Configurational Entropy Relaxation of Silica Glass—Molecular Dynamics Simulations

Vitreous silica was modelled using molecular dynamics (MD). The glass structure was transferred into an undirected graph and decomposed into disjoint structural units that were ideally mixed to calculate the configurational entropy. The Debye relaxation model was suggested to simulate the evolution of entropy during the cooling of the system. It was found that the relaxation of the configurational entropy of MD corresponds to the effective cooling rate of 6.3 × 10⁶ Ks⁻¹ and its extrapolation to 0.33 Ks⁻¹ mimics the glass transition with T_g; close to the experimental value. Debye relaxation correctly describes the observed MD evolution of configurational entropy and explains the existence of freezing-in temperature and the shape of the curve in the transition region.     

高压下TiN的结构转变、弹性和热力学性质的第一性原理计算

 利用基于密度泛函的第一性原理,计算了高压下TiN的结构转变、弹性和热力学性质。计算结果表明：在压力作用下,TiN经历了从NaCl型结构到CsCl型结构的转变,转变压力为348 GPa;TiN的弹性系数随着压力的增加呈线性增加规律。此外,还给出了德拜温度和热容量这两个重要热力学量与温度和（或）压力的依赖关系。     

含过渡金属铜离子的西佛碱型液晶化合物的结构分析

通过红外光谱和拉曼光谱的研究表明,在金属配合物ＶＯＨ１２－Ｃｕ中,尾链以全反式构象存在,且流动性增加,同时由于羰基于上的氧与另一分子中的配位中心离子铜之间存在相互作用,氧向铜提供电子,使分子呈现交错层状排列,既分子Ａ中羰基的位置与分子Ｂ中铜的位置上下对应,分子Ｂ中羰基的位置与分子Ｃ中铜的位置上下对应。     

Molecular dynamics of reactions between (4,0) zigzag carbon nanotube and hydrogen peroxide under extreme conditions

ABSTRACT

Single-wall carbon nanotubes (CNTs) have been suggested as potential materials for use in next-generation gas sensors. The sidewall functionalisation of CNTs facilitates gas molecule adsorption. In this study, density functional theory (DFT)-based ab initio molecular dynamics simulations are performed for a periodic zigzag single-wall (4,0) CNT surrounded by a monolayer of hydrogen peroxide molecules in an attempt to find conditions that favour sidewall functionalisation. The dependency of dynamics on charge states of the system is examined. It is found negative charges favour reactions that result in the functionalisation of the CNT. First principles molecular dynamics of defect formation yields chemically reasonable structure of stable defects, which can be reproduced in CNTs of any diameter and chirality. The explored hydroxyl and hydroperoxyl defects increase conductivity in a large diameter (10,0) CNT, while decrease conductivities in a small diameter (4,0) CNT.     

Pressure and Electrical Resistivity Measurements on Hot Expanded Metals: Comparisons with Quantum Molecular Dynamics Simulations and Average‐Atom Approaches

We present experimental results on pressures and resistivities of expanded nickel and titanium at respective densities of 0.1 g/cm³ and 0.2 g/cm³, and in a range of temperature of 1‐3 eV that corresponds to the warm dense matter (WDM) regime. These data are used to benchmark different theoretical approaches. A comparison is presented between fully 3‐dimensional quantum molecular dynamics (QMD) methods, based on density functional theory, with average‐atom (AA) methods, that are essentially one dimensional. AA methods are used to identify interband transitions and photoionization thresholds. In this regime the evaluation of the thermodynamic properties as well as electrical properties is difficult due to the concurrence of density and thermal effects which directly drive the metal‐non‐metal transition. QMD simulations are also helpful to give a precise estimation of the temperature of experiments which is not directly accessible [1] (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)