Molecular dynamics of MgSiO3 perovskite melting

The melting curve of MgSiO分子动力学 MgSiO3钙钛矿 熔化温度 高压melting temperature, molecular dynamics, high pressureProject supported by the National Natural Science Foundation of China (Grant Nos 10274055 and 10376021),the Natural Science Foundation of Gansu Province, China (Grant No 3ZS051-A25-027) and the Scientific Research Foundation of Education Bureau of Gansu Province, China (Grant No 0410-01).2005-01-125/8/2005 12:00:00 AMThe melting curve of MgSiO3 perovskite is simulated using molecular dynamics simulations method at high pressure. It is shown that the simulated equation of state of MgSiO3 perovskite is very successful in reproducing accurately the experimental data. The pressure dependence of the simulated melting temperature of MgSiO3 perovskite reproduces the stability of the orthorhombic perovskite phase up to high pressure of 130GPa at ambient temperature, consistent with the theoretical data of the other calculations. It is shown that its transformation to the cubic phase and melting at high pressure and high temperature are in agreement with recent experiments.     

4种等温物态方程适用性的比较研究

 通过分析现有文献，认为4种两参数的物态方程在高压物理研究中较为流行，它们是Vinet、Baonza、Morse和Born-Meyer（BM）方程。将这4种方程用于拟合50种材料的实验压缩数据，得出了零压下压缩模量及其一阶压力导数的优化取值，并计算了4种方程的平均压力误差。结果表明，Morse方程的精度最高，对50种材料的平均拟合误差为0.557 6%；BM方程次之，平均误差为0.615 1%；Vinet和Baonza方程的误差大一些，分别为0.788 9%和0.833 3%。对8种典型材料计算了压力误差随压强变化的曲线，所得结果与平均误差的趋势一致，也是Morse方程的精度最高。在宽广压力范围的高压物态方程研究中，推荐使用Morse方程。     

高温高压下Zr2AlC的结构及热学性质的第一性原理

利用密度泛函理论研究了高温高压下Zr₂AlC的结构和热力学性质,计算得到Zr₂AlC的晶格参数与实验值符合较好．研究了Zr₂AlC的弹性常数、体模量、剪切模量和杨氏模量等力学性质随压力变化的趋势．同时研究了维氏硬度随压力的变化趋势．通过计算得到的杨氏模量预测了Zr₂AlC的弹性各向异性．最后,基于准简谐德拜模型,成功预测了Zr₂AlC的德拜温度、热容、热膨胀系数和Grüneisen参数随着压强和温度的变化关系.     

低孔隙度疏松铝的高压声速与冲击熔化

对含微孔洞疏松度m=1.04的疏松铝进行了冲击加载-卸载实验,利用DISAR(distance interferometer system for any reflector)测得了53至99 GPa五个冲击压力下疏松铝/LiF界面粒子速度波剖面,获得了各压力下的纵波声速和其中三个压力点的体波声速,确定出疏松铝的冲击熔化压力约为81 GPa,确定出高压下冲击熔化前的泊松比约为0.372.通过分析,微孔洞明显降低了冲击熔化压力,引起的非谐振效应明显,状态方程计算中考虑非谐效应,非谐因子l 关键词： 低孔隙度 疏松铝 声速 冲击熔化     

液氘Hugoniot曲线的理论计算

 用液体变分微扰理论加量子力学修正和指数6（EXP-6）相互作用势，计算了液氘的一次和二次冲击压力为0～80 GPa范围Hugoniot曲线。计算结果表明：冲击压力与实验符合较好。但冲击温度明显高于实验值，对其偏高的原因进行了探讨。     

Parallel Molecular Dynamics Simulations of Ejection from the Metal Cu and Al Under Shock Loading

Large-scale non-equilibrium molecular dynamics simulations are used to investigate the ejection of the metal under a shock loading. The present work focus on the dynamic process of ejection from the metal Cu and Al surface groove under shock loading, using parallel MD implementation and the Morse potential. The ejected mass coefficient and the size distribution of ejected particles （cluster for atoms） are investigated with changes of the half-angle or the depth of groove and shock strength.     

Pressure dissociation of dense hydrogen

The self-consistent fluid variational model （SFVM） has been used to describe the pressure dissociation of dense hydrogen at high temperatures. This paper focuses on a mixture of hydrogen atoms and molecules and is devoted to the study of the phenomenon of pressure dissociation at finite temperatures. The equation of state and dissociation degree have been calculated from the free energy functions in the range of temperature 2000-10,000K and density 0.02-1.0g/cm^3, which can be compared with other approaches and experiments. The pressure dissociation is found to occur in higher density range, while temperature dissociation is a more gradual effect.     

高密度流体He+H_2混合物物态方程研究

选择高密度流体He +H2 混合物作为研究对象 ,用Ree混合规则和vanderWaals单组分流体变分微扰理论加量子力学一级修正模型编制计算程序。作为对计算模型及其程序的检验 ,首先用已有的α exp 6优化势参数 ,计算了T =30 0K的He +H2 混合物的等温相平衡线 ,得到了与实验值和分子动力学 (MD)数值模拟一致的结果 ,然后进一步计算了 0～ 60GPa和 5 0～ 70 0 0K压力温度范围内的流体He+H2 混合物 (He∶H2 分别为 1∶1、1∶3、3∶1摩尔比 )的高压物态方程。与Monte Carlo模拟数据所进行的比较表明 ,在低温下 ,量子力学修正对热力学量的计算是重要的     

多体相互作用对高压固氦状态方程的影响

采用量子化学自洽场方法及原子团簇理论对固氦晶格原子势能进行多体展开，定量确定了最近邻原子间距R在0.26—0.175nm之间，短程的两体到五体相互作用对晶体结合能及压缩特性的贡献，并对不同间距下多体展开式的截断位置进行了讨论. 结果表明：当最近邻原子间距R在0.26—0.175 nm之间取值，在具有hcp相结构的晶格中，氦原子势能的多体展开式是一收敛的交错型级数，该级数中两体项、四体项为正值，而三体项、五体项为负值. 当R在0.26—0.23 nm之间取值时, 仅保留级数中两体和三体项就能很好地描述晶格原子势能；当R在0.23—0.20 nm之间取值，多体展式需要展开到四体相互作用项；当R在0.20—0.175 nm之间取值必须再加入五体项贡献. 考虑到五体相互作用后，理论计算已能圆满地解释目前固氦等温压缩数据，其最大压力达 60GPa. 关键词： 固氦 状态方程 多体相互作用     

H_2+He流体混合物在部分离解区的物态方程

H2+He流体混合物在高温高压下由于氢的离解化学反应形成由H2,H,He三种粒子构成的混合体系,此时粒子间的相互作用较为复杂,离解能也会由于粒子间的这种复杂相互作用而降低.本文利用自洽流体变分理论来研究部分离解区H2+He流体混合物的高温高压物态方程,模型考虑了各种粒子间的相互作用及由温致和压致效应引起的离解能降低的自洽变分修正,并通过自洽流体变分过程对非理想的离解平衡方程求解得到粒子数密度分布,进而对自由能求导获得体系的热力学状态参量.计算结果与已有的冲击波实验数据、蒙特卡罗模拟及其他理论计算进行了比较.