近饱和区过热水蒸汽性质的分子动力学模拟研究

本文根据统计热力学理论,应用分子动力学模拟方法,对过热水蒸气的热物性进行了研究.选取实验压力为10MPa,利用TIP4P模型和不同经典方程,模拟计算温度从320℃到800℃时水蒸气的压力.从结果可以看出:近饱和区标准维里方程模拟压力值偏离实验值较远,径向状态方程和段式状态方程模拟结果的误差要小于维里方程的.随温度的增加,三种方程的模拟结果均趋于实验值.另外,本文还对影响模拟结果的一些因素进行了详细的讨论.     

临界点附近H2气体的焓和熵的热力学计算

采用热力学方法和Benedict-Webb-Rubin(BWR)方程及相应的参数,计 算了临界温度附近H2气体的绝对焓和熵值,同时拟合了一组高准确度的定压比热容与温度 关联式参数.通过计算值与实验值的比较发现:在临界温度以上,计算值与实验值符合良好 ,对低于临界温度的未饱和H2气体也能计算出较好了结果.但对临界点和温度低于临界温 度、压强高于临界压强的过饱和H2气体的计算误差比较大,说明此时的H2气体的热力学性质特别,应当专门处理.     

高密度氩气的原子间相互作用与状态方程

本文从高密度气体的原子间相互作用出发进行思考,在修改范德瓦耳方程的基础上,提出了一种新的适用于高密度气体的状态方程,并用以对氩气临界区的等温压缩线进行计算,结果与实验值很符合.     

于串珠分子动力学的Cl+XCl→XCl+Cl(X=H，D，Mu)反应的速率系数和动力学同位素效应研究

本文用串珠分子动力学计算典型的重-轻-重反应Cl+XCl→XCl+Cl(X=H,D,Mu)的热速率系数和动力学同位素效应. 对于Cl+HCl反应,串珠分子动力学方法计算得到的结果与实验值之间的高度一致性为串珠分子动力学理论的准确性提供了有力的证据. 串珠分子动力学结果也与其他理论方法的计算结果相一致,包括正则变分过渡态理论和量子动力学方法. Cl+MuCl反应在过渡态附近存在势能双峰,留下了自由能阱,它来自于势能面峰值处反应系统的模式软化. 同时,如此明显的自由能阱表明存在反应共振. 对于Cl+DCl反应,与实验值相比,串珠分子动力学速率系数再次给出了非常准确的结果. 唯一的例外是在312.5 K的温度下,串珠分子动力学速率系数和其他理论方法计算得到的结果彼此接近,却略低于实验值,表明误差可能源于实验测量或本文使用的势能面不够精确.     

基于串珠分子动力学的Cl+XCl→XCl+Cl(X=H,D,Mu)反应的速率系数和动力学同位素效应研究（英文）

本文用串珠分子动力学计算典型的重-轻-重反应Cl+XCl→XCl+Cl(X=H,D,Mu)的热速率系数和动力学同位素效应.对于Cl+HCl反应,串珠分子动力学方法计算得到的结果与实验值之间的高度一致性为串珠分子动力学理论的准确性提供了有力的证据.串珠分子动力学结果也与其他理论方法的计算结果相一致,包括正则变分过渡态理论和量子动力学方法.Cl+MuCl反应在过渡态附近存在势能双峰,留下了自由能阱,它来自于势能面峰值处反应系统的模式软化.同时,如此明显的自由能阱表明存在反应共振.对于Cl+DCl反应,与实验值相比,串珠分子动力学速率系数再次给出了非常准确的结果.唯一的例外是在312.5K的温度下,串珠分子动力学速率系数和其他理论方法计算得到的结果彼此接近,却略低于实验值,表明误差可能源于实验测量或本文使用的势能面不够精确.     

高温高压下液态水声速的研究——不同状态方程的准确性验证

 研究水的热力学状态方程,对于理解地球及行星科学等起着重要作用,但由于高温高压区域的实验数据较少,该区域的状态方程主要依赖于由低压部分外延或分子动力学模拟计算得到。采用布里渊散射技术测量熔解曲线附近液态水的声速,低温区采用电加热系统,高温区采用激光加热布里渊散射系统,分析比较了由实验测量得到的声速值与用经验状态方程计算的结果之间的差别。结果表明,在温度不超过673 K、压力不超过6.0 GPa的范围内,Abramson方程的计算结果与实验测量结果在误差范围内一致,而Saul 和IAPWS-95的预言值比实验测量值偏高,并且温度越高偏差越大。在压力为21 GPa、温度为890~1 100 K时,实验测量出的水的声速比状态方程预言的结果偏高。     

高压下钙钛矿结构MgSiO3的分子动力学研究

利用分子动力学方法,研究了高温高压下钙钛矿结构MgSiO3的状态方程.研究表明,分子动力学模拟结果很好地再现了广泛温度和压强范围内钙钛矿结构MgSiO3的摩尔体积.温度300 K压强上升到120 GPa模拟的钙钛矿结构MgSiO3状态方程和有效的实验结果基本一致.在更高温度和更高压强下模拟的钙钛矿结构MgSiO3状态方程和他人的计算值吻合的很好.另外,还分别计算了温度300 K,900 K,1500 K和2500 K压强上升到120 GPa时MgSiO3的体积压缩率.     

基于傅里叶变换红外光谱法CO2气体碳同位素比检测研究

介绍了基于傅里叶变换红外技术检测CO₂气体碳同位素比的新方法, 详细介绍了如何从HITRAN红外数据库中提取气体标准吸收截面; 介绍了基于非线性最小二乘法反演CO₂气体碳同位素比和整套实验装置的组成及实验步骤. 从理论和实验分析两方面讨论了温度和气压变化对δ¹³CO₂值的影响规律. 对于同一CO₂标准气体, 采用FTIR和同位素质谱法两种技术进行了δ¹³CO₂值对比检测, 两种测量技术的平均值差异仅为0.25%. 从实验结果可以看出, FTIR技术可以实现对CO₂气体碳同位素比的检测.     

MgSiO3状态方程的分子动力学模拟

利用分子动力学方法,研究了高温高压下钙钛矿结构MgSiO_3的状态方程。研究表明,分子动力学模拟结果精确地再现了广泛温度和压强范围内MgSiO_3的摩尔体积。在300 K压强上升到140 GPa模拟的MgSiO_3状态方程和有效的实验值、他人的拟合值以及基于局域密度近似的第一原理计算结果基本一致。并且更高温度和更高压强下模拟的MgSiO_3状态方程和他人的计算值吻合的很好。另外,还分别计算了300、900、2000和3000 K压强上升到120 GPa时MgSiO_3的体积压缩率。     

邮票上的物理学史——气体动理论和统计力学的建立

气体动理论是热学的一种微观理论 ,用分子的运动来解释物质的宏观性质 .它的基本概念有两个 :一是物质由大量分子和原子组成 ;二是热现象是这些分子无规运动的一种表现形式 .首先看气体的状态方程 .状态方程是描写气体状态的宏观参数 ｐ、Ｖ、Ｔ所满足的方程 ,它本身是一项宏观性质 .但是任何微观理论都应当能够导出状态方程 ,而且真实气体的状态方程即范德瓦尔斯方程的确是由气体动理论导出的 .英国化学家玻意耳于 1 662年 (他的邮票见本专栏“真空和大气压”一节 )由实验发现 ,对于一定量的气体 ,在恒定温度下 ,其体积与压力成反比 ,即 ｐ…     

Molecular Dynamics Study of gases H2,D2 and T2

The classical Molecular dynamics simulation has been used to study the equation of state of gas H2,D2 and T2.It has also been investigated that the isotope mass affects on the accuracy of equation of state.Our calculated Iesults show that the classical effect is principal and the isotope mass effects on the equation of state are obvious for the much light gases.At the same time,some useful theoretical data of equation of state for these gases have been provided.It is found that the classical simulation is still effective to the quantum gas.However,the quantum mechanics simulation and the improvement of intermolecular interaction potential are necessary if more accurate computational results are expected.     

氢同位素转动跃迁拉曼频移的理论计算

激光拉曼光谱作为一种非侵入式的物质结构和成分分析的测试手段特别适合于核工业中含氚放射性同位素气体的定性和定量分析。将经典物理学方法和量子力学理论相结合,采用非刚性转子校正模型,计算获得了H2、HD、HT、D2、DT、T26种氢同位素气体分子的转动跃迁拉曼频移,为氢同位素气体拉曼光谱谱峰的辨识提供理论数据。通过实验测定H2和D2的转动拉曼光谱,得到这2种分子的转动拉曼图谱和拉曼频移,并将测定结果与计算结果进行比较,验证了理论计算方法的可行性和正确性。     

Isotope fractionation in aqua-gas systems: Cl(2)-HCl-Cl(-), Br(2)-HBr-Br(-) and H(2)S-S(2-)

We report calculated values of isotope fractionation factors between chlorine, bromine and sulphide hydrated anions and respective gaseous compounds: hydrogen chloride, hydrogen bromide, molecular chlorine and bromine and hydrogen sulphide. For the calculation of the reduced partition function ratios (β-factors) of hydrated Cl(-), Br(-) and S(2-) anions, we used a model of a cluster composed of the considered ion surrounded by two shells of H(2)O molecules. Only the electrostatic interaction between ion and water molecules treated as electric dipoles was taken into account. The β-factors for the gaseous compounds (HCl, Cl(2), HBr, Br(2) and H(2)S) were calculated from vibrational frequencies reported by Urey and Greiff [Isotopic Exchange Equilibria, J. Am. Chem. Soc. 57, 321 (1935)] and Schauble et al. [Theoretical Estimates Equilibrium Chlorine-Isotope Fractionation, Geochim. Cosmochim. Acta 67, 3267 (2003)]. Low-temperature isotope fractionation between chlorine-hydrated anion and hydrogen chloride attains 1.55-1.68‰ (this work), which is in good agreement with experimental data (1.4-1.8‰) [Z.D. Sharp, J.D. Barnes, T.P. Fischer and M. Halick, An Experimental Determination of Chlorine Isotope Fractionation in Acid Systems and Applications to Volcanic Fumaroles, Geochim. Cosmochim. Acta 74, 264 (2010)]. The predicted isotope fractionations for hydrated bromine and HBr, Br(2) gases are very small, 1000?ln α, do not exceed 0.8‰; thus, the expected variations of bromine isotope composition in aqua-gas systems will require enhanced precision for their detection. In contrast, the sulphur isotope fractionation between H(2)S( gas ) and S(2-) attains 6.0‰ at room temperature and drops nearly linearly to 3.1‰ at 350°C.     

稀有气体在电离区压缩特性研究

采用Saha方程加Debye-Hückel修正简单近似模型,给出了稀有气体He, Ne, Ar, Kr, Xe在电离区的物态方程以及离化度.计算结果与已有实验和理论计算进行比较,验证了模型的可靠性.通过对稀有气体等温和冲击压缩特性理论分析,讨论了其压缩和电离的规律性.论证了在稀有气体He, Ne, Ar, Kr, Xe中,气体Xe具有较好的抗压缩性.     

Inert fluorinated gas T1 calculator

The physics of spin-rotation interaction in roughly spherical perfluorinated gas molecules has been studied extensively. But, it is difficult to calculate a spin-lattice relaxation time constant T1 for any given temperature and pressure using the published literature. We give a unified parameterization that makes use of the Clausius equation of state, Lennard-Jones collision dynamics, and a formulaic temperature dependence for collision cross section for rotational change. The model fits T1s for SF6, CF4, C2F6, and c-C4F8 for temperatures from 180 to 360 K and pressures from 2 to 210 kPa and in mixtures with other common gases to within our limits of measurement. It also fits previous data tabulated according to known number densities. Given a pressure, temperature, and mixture composition, one can now calculate T1s for common laboratory conditions with a known accuracy, typically 0.5%. Given the success of the model's formulaic structure, it is likely to apply to even broader ranges of physical conditions and to other gases that relax by spin-rotation interaction.     

Strong isotope effects on the charge transfer in slow collisions of He2+ with atomic hydrogen, deuterium, and tritium

Probabilities and cross sections for charge transfer by He2+ impact on atomic hydrogen (H), deuterium (D), and tritium (T) at low collision energies are calculated. The results are obtained using an ab initio theory, which solves the time-dependent Schr?dinger equation. For the H target, excellent agreement is achieved between the present and previous results. Differences by orders of magnitude are observed between the cross sections for H, D, and T. A method is introduced to separate the contributions of charge-transfer mechanisms due to radial and rotational coupling. The large differences observed for H, D, and T are attributed to isotope effects in the rotational coupling mechanism.     

Application of a new equation of state to energy carriers

A new equation of state (IR EOS) recently reported for liquids and gases has been utilized to predict the densities of some energy carriers at different temperatures, pressures. The ability of IR EOS is examined by comparing its results with experimental data for some energy carriers in homogeneous gas, homogeneous liquid and gas–liquid transition region from low to very high pressures. The IR EOS gives excellent results in homogenous gas and homogeneous liquid region while its predictions in gas–liquid transition have more deviations. The average absolute deviation between calculated and experimental densities for 968 data points of 12 energy carriers is 0.33% over the entire range of data with a maximum pressure of 1000 MPa.     

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.     

自由电子激光振荡器的二维数值模拟

在柱二维坐标系中,用数值计算方法求解了摇摆器内电子模拟方程组、光场方程组和无源区(不含摇摆器区)经非自适应方法变换后的无源旁轴波动方程并设计了相应的程序R_2D。用该程序计算的数值结果表明,在有源条件下二维数值模拟结果与高斯光束近似下的数值模拟结果符合较好;取完全相同的计算条件R-2D程序与洛斯阿拉莫斯实验室(LANL)的FELEX程序的计算结果符合较好;初步认为LANL实验数据的物理图象合理,数值结果基本可信;最后给出北京自由电子激光器(BFEL)的一些物理数据。