论文部分内容阅读
二硼化镁、锐钛矿、氧化镁、氧化锶和超硬材料立方氮化硼以其卓越的物理化学性能引起了从事实验、理论研究的物理学家和地球物理学家的广泛关注。本论文对二硼化镁、锐钛矿、氧化镁、氧化锶和超硬材料立方氮化硼的弹性和热动力学性质进行了研究。对固体的微观电子结构进行研究,我们可以认识固体各种物理性质的本质。固体是由大量的原子组成,每个原子又有原子核和电子,原则上说,直接求解这个多粒子系统的方程是不现实的,必须采用一些近似和简化。在研究中,对二硼化镁、锐钛矿采用了分子动力学与经验势(Buckingham势)相结合的方法;对氧化镁、氧化锶和超硬材料立方氮化硼研究时,采用了第一性原理。对于氧化镁和氧化锶我们分别采用由Zupan和Causa等人为研究碱土金属氧化物常压室温下的弹性常数所提出的全电子基组;对于立方氮化硼我们采用由Orlando等人于1990年运用CRYSTAL程序的Hartree-Fock方法研究常压室温下的钻石、硅、BN、BP、SiC和AlP所提出的全电子基组。我们在高压下使用这些全电子基组时重新进行了优化。在模拟氧化镁和氧化锶性质过程中,我们应用了原子轨道线性组合的Hartree-Fock自洽场(SCF)方法:在模拟立方氮化硼性质过程中,我们应用了基于Becke和Perdew-Wang非局域交换和关联函数的密度泛函理论(DFT)理论。论文首先对多原子体系研究历史与现状和固体力学的研究历史与现状作了简要的回顾;比较全面地阐述了第一性原理计算方法与分子动力学模拟方法相关理论;详细地叙述了晶体的弹性理论,列出了七大晶系的弹性常数矩阵。最后,论文研究了二硼化镁、锐钛矿、氧化镁、氧化锶和超硬材料立方氮化硼的弹性性质、压缩性、热动力学性质。计算结果与已有的实验和理论结果吻合较好。获得一些国内外未见报道的新的物理结果。这些结果不仅对于理解物质在高压下的物理机制有益,而且为物质在高压下的应用奠定了理论基础。在利用GULP程序研究时,为二硼化镁和锐钛矿建立了常温常压下的Buckingham势。并以此势函数较为详细地研究了常温常压下二硼化镁和锐钛矿的弹性性质、压缩性以及部分热动力学性质。其结果与已有的实验和理论结果吻合得较好,并首次获得了锐钛矿常温常压下的弹性常数。研究发现,由于c轴方向的Mg-B结合键比a或b轴方向的B-B结合键更容易压缩,二硼化镁的c轴方向比a或b轴方向更容易压缩;由于c轴方向的Ti-O和O-O结合键比a或b轴方向的Ti-Ti结合键更容易压缩,锐钛矿的c轴方向比a或b轴方向b也更容易压缩。在锐钛矿的模拟中,我们把每个离子的部分电荷作为拟合参数,这解决了由于电荷的敏感性引起的电荷在局域环境与Ti和O原子之间相互转移带给分子动力学模拟的困难。这种处理方法对于以后处理类似问题有一定的参考意义的。研究发现,由于随着压强增大相互作用的近邻原子增多,二硼化镁和锐钛矿的压缩率随着压强增大非线性变化加强。最后指出了,在12GPa范围内,我们为常温常压下二硼化镁和锐钛矿建立的Buckingham势函数是基本可靠的,但是随着压强的增高此势函数的可靠性必然会降低,我们应该考虑多体势。其原因是用两体势的方式来处理多原子体系的晶体本就是一种简化。随着压强增大,原子间相互作用更为复杂,此时两体势的精度必然要降低。在利用CRSTAL程序研究氧化镁、氧化锶和超硬材料立方氮化硼时,我们首先研究了氧化镁、氧化锶和立方氮化硼常温常压下的平衡晶格常数和弹性常数并与相关的实验和理论值进行了比较,结果表明我们的计算结果与已有的实验值或理论值吻合较好。值得注意的是,我们首次计算了常温常压下的B2型MgO与SrO的弹性常数,并利用弹性常数应满足弹性稳定性判定原则,在理论上我们指出常温常压下的B2型MgO与SrO是不稳定的。在上面工作的基础上,通过三阶Brich-Murnaghan状态方程了获得对应的P-V关系。据此,我们计算了MgO、SrO和立方氮化硼高压下的弹性常数。在我们所讨论的压强范围内B1型MgO与SrO的c44计算值并没有消失,这一点是与Tsuchiya等人的结果不一致。不过它与Hunter等人的相变理论是吻合的。研究发现,由于Mg2+和O2-离子的半径差比Sr2+和O2-的半径差更大并且Mg2+离子的电负性比Sr2+离子的电负性更大,这导致了高压下MgO中的非中心多体力比SrO的更大。通过Cauchy关系,发现立方氮化硼的非中心力随压强增加而增大,并且几乎成线性关系变化。据此,我们认为,在高压下随着非中心多体力场作用的加大,高压下处理立方氮化硼时,我们必须考虑非谐作用的影响,也有必要讨论三阶或四阶弹性常数;在研究立方氮化硼弹性各向异性的0-200GPa范围内中,我们没有发现呈向弹性各向同性的压强点;发现立方氮化硼的c11比c12和c44对压强作用更为敏感;运用立方晶体弹性常数的稳定性判定原则:c11+2c12>0,c44>0,c11-c12>0研究立方氮化硼的力学稳定性时,并没有发现不稳定点。同时,研究中,我们首次详细的讨论了立方氮化硼的弹性对称性;利用弹性与热动力学性质的关系,首次计算了B1型MgO与SrO高压下的德拜温度、纵波和横波弹性波波速;发现,德拜温度、纵波和横波弹性波波速随压强的增大而增大;首次计算了立方氮化硼高压下的德拜温度、纵波和横波弹性波波速、热膨胀系数、热容和泊松比率;发现,纵波和横波弹性波波速、德拜温度和泊松比率随压强的增大而增大,热膨胀系数和热容随压强的增大而减小(这是由于增加压强的效应相当于降低温度)。这些研究无疑对MgO、SrO和立方氮化硼的进一步理论探索提供了丰富的数据。