液H2+D2混合冲击压缩特性的理论计算

在液体变分微扰理论加量子力学修正基础上,用指数6有效作用势对液H2及其同位素液D2的混合物态方程进行了计算,首次给出了液H2+D2混合冲击压缩Hugoniot曲线和冲击温度.     

固氢固氦等温压缩线的理论计算

用原子团簇理论和完全量子力学计算,研究了固氦和固氢的等温压缩线.在固氢的分子间相互作用势计算中,采用了球形势近似和等效氦原子模型.计算结果与100GPa以下的实验等温线相比较,发现两者有很好的符合. 关键词：     

高密度氦状态方程的分子动力学模拟研究

 用分子动力学模拟和与固氦实验等温线比较的方法,确定了高密度氦指数6势函数中劲度系数α的优选值为12.7。然后我们用这个新的α值计算了fcc固氦的状态方程及其径向分布函数,发现当把ρ限定为1.60 g/cm³时,其径向分布函数的第二个峰将在2 000 K到3 040 K区间消失,表明此时固氦的长程有序性降低,或者说发生了固-液相变过程。     

采用从头计算方法研究液态氦原子间等效对势

 提出了一种具有体心立方排列的原子团簇He₉模型,用于模拟液氦中邻近原子的近程平均分布特征,并运用量子从头计算方法和原子团簇理论计算技术,首次从理论上计算了液态氦原子间等效作用对势。计算给出的等效对势能较好地再现液氦的等温压缩线及冲击压缩线,所提出的液体等效结构的原子团簇模型能比较合理地描述液氦近程结构特征以及原子间多体相互作用规律。还将理论对势函数与经验Exp-6势进行了比较。     

Fe/FeO/FeS混合物的Hugoniot线研究

 在50～210 GPa的压力范围内,用二级轻气炮和电探针技术对平均密度ρ₀＝(6.69±0.06) g/cm³的三组元Fe/FeO/FeS（质量分数分别为58.96%、35.83%、5.21%）混合物的Hugoniot线进行了实验测量,所得的Hugoniot参数为：C₀＝(3.97±0.07) km/s,λ＝1.58±0.03。该混合物Hugoniot线的实测结果与用体积可加性原理计算得到的相同组分混合物的Hugoniot线的符合性很好;根据实验数据还计算了混合物的0 K等温压缩线,发现它与体积可加性原理对单质Fe、FeO和FeS的0 K等温压缩线的计算结果相一致,证明了实验结果的合理性与体积可加性原理的适用性,也表明了此混合物在冲击压缩过程中没有发生过可察觉的化学反应。研究结果亦为今后对外地核各种浓度比Fe-O-S体系候选组分高温高压物态方程及物性的进一步研究奠定了基础。     

微喷射粒子团簇探测算法及其应用

在团簇多重标号技术基础上,提出了团簇探测的新算法,并以此作为大尺度二维和三维分子动力学模拟冲击加载金属材料微喷射粒子团簇分析的子程序.以二维为例描述了算法的逻辑结构,测试分析了算法的运行效率,并将其应用到二维分子动力学模拟冲击加载金属Cu和Al微喷射粒子的团簇统计分析中,对微喷射粒子团簇尺寸分布得到一些规律性认识.     

面心立方金属铈的高压电子结构研究

本文简要介绍了我们开发的基于密度泛函理论和动力学平均场理论（DFT+DMFT）的全自洽的的计算程序,并运用该程序研究了面心立方金属铈的等结构相变. 动力学平均场采用的是目前公认的最为准确和高效的连续时间蒙特卡罗方法作为杂质求解器. 我们首先分别计算了 和 相的电子态密度,该结果可以正确的描述两个相的电子结构特征. 然后计算了不同温度下的相变前后的电子结构,从中可以得到电子结构在不同温度下的相变过程,与之前的实验结果一致. 这一工作为以后进一步的金属铈的等结构相变研究打下了基础.     

Extension of Compressibility-Route Cubic Equations of State and the Radial Distribution Functions at Contact to Multi-Component Hard-Sphere Mixtures

The equation of state(EOS) for hard-sphere fluid derived from compressibility routes of Percus-Yevick theory(PYC) is extended. The two parameters are determined by fitting well-known virial coefficients of pure fluid.The extended cubic EOS can be directly extended to multi-component mixtures, merely demanding the EOS of mixtures also is cubic and combining two physical conditions for the radial distribution functions at contact(RDFC) of mixtures.The calculated virial coefficients of pure fluid and predicted compressibility factors and RDFC for both pure fluid and mixtures are excellent as compared with the simulation data. The values of RDFC for mixtures with extremely large size ratio 10 are far better than the BGHLL expressions in literature.     

稠密氦气辐射不透明度的实验研究

用二级氢气炮作为冲击压缩加载工具和多通道瞬态辐射高温计作为主要测量系统,测量了冲击压缩氦气等离子体的光辐亮度历史(初始温度293 K,初始压力为0.6 MPa和1.2 MPa两种).根据实测记录信号波形的有关特征量,计算得到了氦等离子体的光谱吸收系数κ(λ)和冲击压缩铝基板表面的光反射率R.结果发现:受冲击铝基板表面的光反射率～0.4(比其初始反射率0.8约降低一半),与Erskine的数据相同;对于光谱吸收系数的实测数据,本文目前未能给出合理的物理解释.     

Simulated Melting Curve of NaC1 up to 200 kbar

The melting curve of NaCI is studied up to 200 kbar by means of the shell-model molecular dynamics method, usiug massive shell core interaction potentials. The model for the interatomic interaction is shown to produce reasonable results at a wide range of pressures in bulk transitions. The pressure dependence of the melting curve of NaCl was calculated and the result was modified on the assumption of overheating due to the small system size and small time scale simulation. The final result is in good agreement with the corrected experimental values, accounting for melting mechanisms such as surface heating or superheating. Therefore, it is believed that bulk transition simulation at constant pressure indeed provides a useful tool for studying the melting transition.