基于原子间相互作用的低温硅单晶负热膨胀机制的研究

运用晶格热膨胀的微扰理论,推导了硅晶体的热膨胀系数与原子间两体作用和三体作用的三阶力常数之间的关系公式,在此基础上对热膨胀系数的实验数据进行拟合,计算了硅晶体内的原子间两体和三体的三阶力常数,发现三体作用的三阶力常数为正数,是硅晶体在低温下具有负热膨胀性质的根本原因。计算与分析的结果表明,运用Stillinger-Weber模型得到的三阶力常数为负数,据此不可能计算得到低温下的负热膨胀系数,因此应该对该模型进行修正。     

硅纳米晶体薄膜热膨胀性质的分子动力学研究

基于Stillinger-Weber势对硅纳米晶体薄膜的热膨胀性质进行了分子动力学模拟.研究表明,硅纳米晶体薄膜表面层原子的二聚现象引起薄膜收缩,而原子之间的非和谐势能引起薄膜膨胀;在约400K以下的低温段,由于硅纳米晶体薄膜表面层原子发生二聚的原子数目随温度的升高而明显增多,而原子间非和谐势能很小,故此时二聚主导热膨胀性质,热膨胀系数为负;在高温段(约400K以上),由于发生二聚的原子数目随温度升高不再显著地增加并渐趋于稳定,而原子间非和谐势能逐渐显著并主导热膨胀性质,故热膨胀系数为正.     

硅纳米晶体薄膜热膨胀性质的分子动力学研究

基于Stillinger-Weber势对硅纳米晶体薄膜的热膨胀性质进行了分子动力学模拟。研究表明,硅纳米晶体薄膜表面层原子的二聚现象引起薄膜收缩,而原子之间的非和谐势能引起薄膜膨胀;在约400K以下的低温段,由于硅纳米晶体薄膜表面层原子发生二聚的原子数目随温度的升高而明显增多,而原子间非和谐势能很小,故此时二聚主导热膨胀性质,热膨胀系数为负;在高温段（约400K以上）,由于发生二聚的原子数目随温度升高不再显著地增加并渐趋于稳定,而原子间非和谐势能逐渐显著并主导热膨胀性质,故热膨胀系数为正。     

由热学性质获取氩晶体原子间各阶力常数

本文基于晶格动力学和量子力学微扰理论推导了氩晶体的热膨胀系数和比热与原子间相互作用的各阶力常数之间的关系公式,在此基础上根据热膨胀系数和比热的数据计算了氩晶体内的原子间相互作用的各阶力常数,并根据这些力常数绘制了原子间相互作用势能曲线,经比对发现该势能曲线与Morse势能曲线能较好吻合,这表明,本文提出的从热膨胀系数和比热获取各阶力常数的方法是正确的。     

由热学性质获取氩晶体原子间各阶力常数

本文基于晶格动力学和量子力学微扰理论推导了氩晶体的热膨胀系数和比热与原子间相互作用的各阶力常数之间的关系公式,在此基础上根据热膨胀系数和比热的数据计算了氩晶体内的原子间相互作用的各阶力常数,并根据这些力常数绘制了原子间相互作用势能曲线,经比对发现该势能曲线与Morse势能曲线能较好吻合,这表明,本文提出的从热膨胀系数和比热获取各阶力常数的方法是正确的。     

纳米晶体线的热膨胀性质

将原子间相互作用势的非谐项作为微扰,运用声子数表象中的晶格原子振动位移和晶格振动哈密顿公式,推导了纳米晶体线的热膨胀系数公式,并进行了数值计算.研究结果表明,纳米晶体线越细,其热膨胀系数越大.     

硅晶体中点缺陷结合过程的分子动力学研究

采用分子动力学方法模拟研究了硅晶体中的空位和间隙原子的结合过程. 研究中采用了Stilliger-Weber三体经验势描述原子间的相互作用, 系统分别在低温300K和高温1400K进行弛豫. 计算中发现空位和间隙原子倾向于通过<111>方向结合，而<110>方向上存在着势垒. 通过势垒值的计算, 对Tang和Zawadzki势垒计算值的差异进行了解释. 关键词： 分子动力学 空位与间隙原子 扩散     

HMX不同晶型热膨胀特性及相变的ReaxFF分子动力学模拟

HMX（octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine）晶体不同晶型对温度的响应对于深入认识含能材料在外加载荷下的感度、稳定性、相变等具有重要意义．采用ReaxFF—lg势函数和等温等压分子动力学（NPT-MD）方法,研究了β-,δ-和α-HMX三种固态晶型在T=303—503K内的晶体结构和分子结构．结果表明,ReaxFF—lg势函数能够合理地描述HMX晶体的热膨胀行为,计算得到的晶体结构和热膨胀系数与实验数据比较符合．三种晶型的线膨胀系数表明,β-HMX晶体具有明显的热膨胀各向异性;δ-HMX晶体C轴的热膨胀特性与a,b轴略有不同;α-HMX晶体三个方向的热膨胀特性基本相同．三种晶型的体膨胀系数为8-HMX〉α-HMX〉β-HMX,表明δ-HMX晶体对温度最敏感,这可能是δ-HMX晶体具有更高感度的原因．在T=443K时,β-HMX晶体的晶胞参数突变,且部分分子构型由“椅式”转变为“船式-椅式”．采用两相热力学模型（2PT）得到的亥姆赫兹自由能表明,β—HMX晶体在T=423—443K时发生相变．β-,δ-和α-HMX晶体分别在T=303—423K,T=443—503K和T=363-423K内稳定存在．     

低温相H_2S晶体振动光谱计算

本文介绍了计算分子晶体基频振动的力常数模型并应用于Ｈ２Ｓ晶体低温相的拉曼，红外光谱分析。使用１１个可调参数（包括３个分子内部力常数和８个分子间作用力参数）计算得到６９个光学模式的振动频率及其对应的本征矢量。计算结果与实验值相符甚好。对已观测到的谱线作了全面指认。并给出了未观测到谱线的位置。计算表明，晶体中氢键为最重要的分子间相互作用。它使Ｈ２Ｓ分子Ｓ—Ｈ键力在晶体中减弱７％。     

����������ѧ���ʵ��ܶȷ��������о�

运用密度泛函理论方法,采用MS中CASTEP模块计算了铀晶体的能量随体积变化的关系,并绘制了E-V曲线。运用准谐振子德拜模型计算了压强从0~100GPa、温度从0~2000K下的晶体铀的相对体积,发现相对体积随着温度和压强的增加而减小。并成功得到了热膨胀系数、热容、以及体弹模量随温度和压强的变化关系。预测了铀的热力学性质。为铀晶体在高温高压下的实验研究提供了理论依据。     

A continuum thermal stress theory for crystals based on interatomic potentials

This paper presents a new continuum thermal stress theory for crystals based on interatomic potentials.The effect of finite temperature is taken into account via a harmonic model.An EAM potential for copper is adopted in this paper and verified by computing the effect of the temperature on the specific heat,coefficient of thermal expansion and lattice constant.Then we calculate the elastic constants of copper at finite temperature.The calculation results are in good agreement with experimental data.The thermal stress theory is applied to an anisotropic crystal graphite,in which the Brenner potential is employed.Temperature dependence of the thermodynamic properties,lattice constants and thermal strains for graphite is calculated.The calculation results are also in good agreement with experimental data.     

General qualifications for the sign of thermal expansion coefficients at low temperatures

It is shown by arguments of plausibility that when interatomic forces include mainly central two body interactions, then the sign of the thermal expansion coefficient for isotropic solids at low temperature can be negative only for solids that lack symmetry of inversion. It is deduced that for RbBr and RbI, angular three-body forces play an important role in the nearest neighbor interactions.     

Negative thermal expansion in crystals with the zincblende structure: an EXAFS study of CdTe

The extended x-ray absorption fine structure (EXAFS) has been measured at both the K edges of cadmium and tellurium in CdTe, from liquid helium to room temperature, in order to investigate the local thermodynamic behaviour. The temperature dependences of the structural parameters obtained from the separate analysis of the two edges are perfectly consistent. The positive contribution to the thermal expansion due to the bond stretching and the negative contribution due to the tension effects are disentangled and quantified in terms of the bond thermal expansion and the perpendicular mean square relative displacement. The comparison with previous EXAFS results for Ge and CuCl shows that relevant correlations can be established between a number of local parameters measured by means of EXAFS and the properties of the lattice negative thermal expansion of tetrahedrally bonded semiconductors. The effective force constants derived from the EXAFS are compared with the force constants of a valence force field model.     

Analysis of anharmonic and thermoelastic properties of non-centrosymmetric crystals

In the present paper an analytical potential form is used for overlap repulsive energy, derived by Harrison from quantum mechanical considerations, along with the composite effect of three-body forces and intersublattice displacement. The short-range overlap parameters in Harrison's potential form have direct correlation with the valence state energies for outermost electrons. The potential model is applied to calculate the third and fourth order elastic constants, first and second pressure derivatives of second order elastic constants, Grüneisen parameter and its volume dependence, Anderson parameter, and thermal expansion coefficient for three non-centrosymmetric crystals, viz. CaF₂, SrF₂ and BaF₂. The calculated values of various physical quantities are found to be in good agreement with experimental data.The authors are grateful to Dr. Mansour Khalef, the Head of Physics Department, TNRC, Tajura (Tripoli) for the facilities and encouragements.     

A general Mobius inversion transform formula for hcp lattices and its application

In this paper,we present a general Mobius inversion transform formula for hcp lattices.This formula can be applied to hcp lattices with a non-ideal c/a value and to obtain the pair potential between atoms in these lattices from the cohesive energy.Also,the three-body interaction among atoms in the lattices can be taken into account in the method.This method gives a useful means to obtain interatomic interactions in the interatomic force model.The method has been applied to zinc,and the pair potential obtained is used to calculate the phonon dispersion relations for some high-symmetry directions.It is found that,by properly considering a three-body interaction,one can acquire satisfactory results.     

An approximate multi-property empirical isotropic interatomic pair potential for the krypton dimer

An approximate empirical isotropic interatomic potential for krypton interaction is developed by simultaneously fitting the Morse-Morse-Morse-Spline-van der Waals potential form to the pressure second virial coefficient, viscosity, thermal conductivity and depolarized interaction-induced light scattering data. Absolute zeroth and second moments of the two-and three-body spectra, the pressure third virial coefficient and isotopic thermal diffusion factor have been measured and compared with theoretical calculations using various models for the interatomic potential. The results show that it is the most accurate potential yet reported for this system. The use of the new potential in lattice sum calculations yields good results for several properties of solid krypton.