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.     

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.     

Application of the three body force shell model to the lattice dynamics of calcium oxide

The lattice dynamics of GaO has been studied on the basis of the three body force shell model, which takes into account the effect of many body interactions in the lattice potential. The dispersion curves obtained by plottingω vsq agree fairly well with the experiments. It is concluded that the value of the molecular electronic polarizability of the solid must be must small than that determined experimentally which suggests that the interaction system in the solid may have a substantial covalent character.     

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

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

壳层模型分子动力学在CaO研究中的应用

利用壳层分子动力学方法结合有效的对势,研究了高压条件下CaO的熔化曲线。研究表明,分子动力学模拟结果精确地再现了广泛压强范围内CaO的状态方程。研究中考虑了分子动力学模拟熔化存在的过热现象,通过晶体的现代熔化理论,对CaO的分子动力学模拟熔化温度进行了修正,获得了高温高压下CaO正确的熔化温度。因此,常压下引入壳层模型的分子动力学为研究物质熔化提供了一个很好的方法,这种方法可进一步推广到其它物质的高压熔化研究中。     

非晶态石英的变电荷分子动力学模拟

基于迭代变电荷方法,用分子动力学模拟了非晶态石英的结构与振动特征.首先利用熔化-猝火方法得到了非晶石英的平衡结构.在此基础上获得了体系不同原子对之间的对关联函数、键角分布函数和振动频谱等,结果与实验数据均符合较好.变电荷方法的计算结果表明,非晶石英体系内粒子的电荷与石英晶体内粒子电荷显著不同,并且出现了较大的涨落. 关键词： 分子动力学 变电荷 非晶石英     

Relaxation dynamics and power spectra of liquid water: a molecular dynamics simulation study

We present molecular dynamics simulations of liquid water at normal and supercooled conditions. Autocorrelation functions (ACFs) of several structural quantities and their fourier transforms are obtained and analysed. Structural correlations and relaxation times increase linearly with degree of supercooling. Power spectra of ACFs show increase in librational motion of liquid water with cooling. These modes intensify with supercooling because of structuring and ordering of water molecules. Overall, liquid water structure is homogenous over the temperatures and pressures studied and undergoes fluctuation–dissipation in its local-density variations [English and Tse, Phys. Rev. Lett. 106, 037801 (2011)].     

非绝热分子动力学的量子路径模拟

基于近期发展的经典-量子混合模拟非绝热分子动力学的量子路径方案,本文对5个典型势能面模型进行了模拟,包括单交叉模型、双交叉模型、拓展耦合模型、哑铃模型以及双弓模型.由于难以在严格意义上得到退相干速率,数值模拟中,我们比较了三个不同的退相干速率公式,包括冻结高斯波包近似退相干速率、能量分辨速率以及力分辨速率.在模拟过程中,我们恰当地处理了势能面跳跃时的能量守恒和力的反向问题.通过与全量子动力学模拟的精确结果进行对比发现,对于结构较简单的势能面模型,三种退相干速率都能得到较好的结果;然而对于较复杂的势能面模型,由于复杂量子干涉的原因,与其他混合经典-量子动力学方案类似,量子路径方案仍然难以得到较准确的结果.如何发展更加有效的混合经典-量子模拟方案,是未来研究的重要课题.     

基于分子动力学模拟的主链型液晶聚合物的新模型

在传统的Gay-Berne (GB)/Lennard-Jones (LJ)模型的基础上，发展了一种用于模拟半刚性主链型液晶聚合物(LCP)的分子级模型，命名为Solo-LJ-SP-GB 模型.单一的LJ联合体和非线性弹簧被用于描述LCP分子中的间隔体.用分子动力学模拟半刚性主链型LCP系统(该系统由169条分子链组成，每两个刚性体之间的间隔体个数为6)时，该模型所需的计算时间不到传统的GB/LJ 模型所需时间的十分之一，大大地提高了计算效率.通过采用该模型模拟半刚性主链型LCP的相变问题，观察到了与半刚性主链型LCP分子中间隔体个数相关的热力学的奇偶效应以及从等方相到向列相的相转变过程.这些模拟结果与当前的试验结果相当符合，从而表明了该模型可以较为准确地描述出半刚性主链型LCP的结构特性. 关键词： Solo-LJ-SP-GB模型 液晶聚合物 分子动力学模拟     

分子动力学模拟方法在非线性光学中的应用

有机分子非线性光学材料在频率转换、电光调制和双光子吸收等方面具有重要的应用.介绍了近年来分子动力学模拟方法在有机分子非线性光学性质理论研究中的主要应用,包括电场极化效应、局域场因子、非线性极化率和双光子吸收等.此外,结合最新的科研工作,介绍了分子动力学模拟方法在溶剂效应和聚集效应研究中发挥的重要作用. 关键词： 分子动力学模拟 溶剂效应 聚集效应 双光子吸收