基于第一性原理的MgZnO电子结构解析与ZnO高压相变机制探究

基于第一性原理的MgZnO电子结构解析与ZnO高压相变机制探究一、引言1.1研究背景与意义在材料科学领域，探索新型材料并深入理解其物理性质是推动科技进步的关键。MgZnO和ZnO作为两种重要的宽禁带半导体材料，因其独特的物理性质和广泛的应用前景，近年来受到了科研人员的高度关注。ZnO是一种Ⅱ-Ⅵ族化合物半导体，具有六方纤锌矿结构，在自然界中广泛存在，如锌矿石和炉甘石等。其室温下禁带宽度约为3.37eV，激子束缚能高达60meV，这使得ZnO在光电器件应用中展现出巨大的潜力。由于其良好的化学稳定性和生物相容性，ZnO在传感器领域也备受青睐，可用于制备气体传感器、生物传感器等，能够对环境中的有害气体或生物分子进行高灵敏度检测。同时，ZnO的压电性能使其在压电换能器、表面声波器件等方面发挥重要作用，实现机械能与电能之间的高效转换。随着对材料性能要求的不断提高，单一的ZnO材料逐渐难以满足所有应用场景的需求。为了进一步拓展ZnO的应用范围，人们通过将Mg元素掺入ZnO晶格中，形成MgZnO合金材料。MgO同样是一种重要的宽禁带半导体，其禁带宽度高达7.8eV。当Mg原子取代ZnO晶格中的Zn原子时，MgZnO合金的禁带宽度可以在ZnO的3.37eV到MgO的7.8eV之间连续调节，这种特性使得MgZnO在紫外光电器件领域具有独特的优势。通过精确控制Mg的掺杂比例，可以制备出满足不同波长需求的紫外发光二极管和探测器，从而实现对紫外光的高效发射和探测。此外，MgZnO的硬度和化学稳定性相较于ZnO也有显著提升，使其在高温、高腐蚀等恶劣环境下的应用成为可能。然而，要充分发挥MgZnO和ZnO材料的性能优势，深入理解它们的电子结构和ZnO在高压下的相变行为至关重要。电子结构决定了材料的电学、光学和磁学等基本性质，而ZnO在高压下的相变则会显著影响其物理性能，进而影响其在高压环境相关领域的应用。传统的实验方法虽然能够直接测量材料的一些宏观性质，但对于微观层面的电子结构和原子尺度的相变机制研究存在一定的局限性。第一性原理计算方法则为我们提供了一种从原子和电子层面深入探究材料性质的有效途径。第一性原理计算基于量子力学原理，以薛定谔方程为基础，通过求解多电子体系的哈密顿量，在不依赖任何经验参数的情况下，精确计算材料的电子结构、晶格常数、弹性常数、能带结构、态密度等物理量。这种方法能够深入揭示材料内部的原子间相互作用和电子分布规律，为解释材料的宏观性质提供微观理论依据。在MgZnO和ZnO材料的研究中，第一性原理计算可以帮助我们理解Mg掺杂对ZnO电子结构的影响机制，预测不同Mg含量下MgZnO的光学、电学性能变化，以及探究ZnO在高压下的结构相变路径和相变压力等关键信息。这些理论计算结果不仅能够为实验研究提供指导，减少实验的盲目性和试错成本，还能为新型光电器件、传感器、高压材料等的设计和开发提供理论基础，推动相关领域的技术创新和发展。1.2国内外研究现状1.2.1MgZnO电子结构研究现状MgZnO作为ZnO基合金材料，其电子结构研究一直是材料科学领域的重要课题。国内外众多学者通过实验和理论计算相结合的方式，对Mg掺杂ZnO后的电子结构变化展开了深入研究。在早期研究中，实验方面主要通过X射线衍射（XRD）技术来分析MgZnO的晶体结构。结果表明，随着Mg含量的增加，晶格常数呈现规律性变化。由于Mg原子半径（0.072nm）小于Zn原子半径（0.074nm），当Mg原子替代ZnO晶格中的Zn原子时，会导致晶格收缩，晶格常数a和c均逐渐减小。这种晶格常数的变化不仅影响材料的微观结构，还对其电子结构产生深远影响。在理论计算方面，基于第一性原理的密度泛函理论（DFT）被广泛应用。Zhang等运用平面波赝势方法，在广义梯度近似（GGA）下对不同Mg含量的MgZnO进行计算，结果显示随着Mg含量的增加，体系的总能量降低，结构稳定性增强。这是因为Mg-O键的键能大于Zn-O键，Mg的掺入使得体系的能量更趋于稳定状态。禁带宽度是衡量半导体材料性能的关键参数之一，对于MgZnO而言，其禁带宽度随着Mg含量的变化备受关注。实验上，通过光致发光（PL）光谱和紫外-可见吸收光谱等手段，测量发现MgZnO的禁带宽度随Mg含量的增加而增大。例如，Wang等制备了一系列不同Mg含量的MgZnO薄膜，利用PL光谱分析得出，当Mg原子摩尔分数从0增加到0.3时，禁带宽度从3.25eV增大到3.85eV。在理论计算方面，Sun等采用基于DFT的第一性原理计算，同样得到了禁带宽度随Mg含量增加而增大的结论，并且通过态密度分析发现，导带中Zn-3d态向高能端移动是导致带隙增大的主要原因。电子态密度分析对于理解MgZnO的电子结构和电学性质具有重要意义。Li等通过第一性原理计算研究了MgZnO的电子态密度，结果表明，随着Mg含量的增加，Zn-O键的共价性减弱，离子性增强。这是由于Mg的电负性（1.31）小于Zn（1.65），Mg的掺入使得O原子周围的电子云密度发生变化，导致化学键性质改变。这种化学键性质的变化进一步影响了材料的电学和光学性质，如载流子迁移率和光吸收特性等。综上所述，目前关于MgZnO电子结构的研究已经取得了较为丰富的成果，对Mg含量与晶格常数、禁带宽度、电子态密度等之间的关系有了较为深入的认识。然而，在一些方面仍存在研究空间，例如，如何精确控制Mg在ZnO中的掺杂分布以实现对电子结构的精准调控，以及MgZnO在复杂环境下电子结构的稳定性研究等，这些问题有待进一步探索。1.2.2ZnO高压相变研究现状ZnO在高压下的相变行为是凝聚态物理和材料科学领域的研究热点之一，对于深入理解材料的基本物理性质和开发新型高压材料具有重要意义。早期的实验研究主要通过金刚石对顶砧（DAC）技术，利用X射线衍射（XRD）和拉曼光谱等手段，对ZnO在高压下的结构变化进行探测。研究发现，在常压下，ZnO具有六方纤锌矿结构（B4相），空间群为P6_3mc。当压力增加到一定程度时，ZnO会发生结构相变。例如，在压力达到约10-12GPa时，ZnO会从六方纤锌矿结构转变为岩盐矿结构（B1相），空间群变为Fm\overline{3}m。这种相变是一种一级相变，伴随着明显的结构重排和体积变化，相变过程中原子的配位数从4增加到6。随着高压实验技术的不断发展，同步辐射光源在ZnO高压相变研究中得到广泛应用。通过同步辐射XRD，可以更精确地测量ZnO在高压下的晶格参数变化和相变压力。一些研究表明，ZnO的相变压力会受到样品的制备方法、颗粒尺寸等因素的影响。例如，纳米结构的ZnO由于其表面效应和量子尺寸效应，相变压力会低于块体ZnO。此外，利用高压拉曼光谱可以探测ZnO在相变过程中晶格振动模式的变化，为研究相变机制提供重要信息。在六方纤锌矿结构向岩盐矿结构转变过程中，拉曼峰的位置和强度会发生明显变化，这与原子间的键长、键角以及晶体对称性的改变密切相关。在理论计算方面，第一性原理计算为研究ZnO高压相变提供了重要的理论支持。通过基于密度泛函理论（DFT）的计算方法，可以精确计算ZnO在不同压力下的晶体结构、能量、电子结构等性质，从而深入探讨相变机制和相变路径。Chen等利用第一性原理计算研究了ZnO在高压下的相变过程，通过比较不同结构相的总能量，确定了相变压力，并分析了相变过程中电子结构的变化。研究发现，在相变过程中，ZnO的电子云分布发生改变，导致化学键的性质和强度发生变化，进而驱动了结构相变。除了六方纤锌矿结构和岩盐矿结构之间的相变，研究还发现ZnO在更高压力下可能会发生其他结构相变。例如，在压力达到约50-60GPa时，ZnO可能会从岩盐矿结构转变为正交相（Cmcm）。这种高压下的多型相变使得ZnO的高压行为变得更加复杂，也为研究带来了更多挑战和机遇。此外，ZnO在高压下的性能研究也取得了一系列成果。高压会显著影响ZnO的电学、光学和力学性能。在电学性能方面，随着压力的增加，ZnO的电导率会发生变化，这与相变过程中电子结构的改变密切相关。在光学性能方面，高压会导致ZnO的禁带宽度发生变化，从而影响其光吸收和发光特性。在力学性能方面，高压下ZnO的硬度、弹性模量等参数会发生显著变化，这些变化对于理解ZnO在高压环境下的应用具有重要意义。尽管目前在ZnO高压相变研究方面已经取得了丰硕的成果，但仍存在一些问题有待进一步解决。例如，对于ZnO在高压下复杂的多型相变机制以及相变过程中原子动力学行为的理解还不够深入；在实验研究中，如何更精确地测量高压下ZnO的热力学性质和动力学过程等。这些问题的解决将有助于进一步完善对ZnO高压相变行为的认识，并为其在高压材料和器件领域的应用提供更坚实的理论基础。1.3研究目标与内容1.3.1研究目标本研究旨在通过第一性原理计算方法，深入系统地探究MgZnO的电子结构以及ZnO在高压环境下的相变行为，具体目标如下：精确计算不同Mg含量的MgZnO合金的电子结构，包括能带结构、态密度、电荷密度分布等，揭示Mg掺杂对ZnO电子结构的影响机制，为调控MgZnO的电学、光学性能提供理论依据。全面研究ZnO在高压作用下的结构相变过程，确定其相变压力、相变路径以及不同相结构的稳定性，深入分析高压相变对ZnO电子结构和物理性能的影响，完善对ZnO高压行为的理论认识。基于计算结果，预测MgZnO和ZnO在不同条件下的性能变化，为其在光电器件、传感器、高压材料等领域的实际应用提供有价值的理论指导，推动相关材料和器件的优化设计与开发。1.3.2研究内容为实现上述研究目标，本研究将围绕以下几个方面展开：MgZnO电子结构的第一性原理计算：模型构建：采用量子力学计算软件，构建不同Mg含量的MgZnO超晶胞模型。考虑到Mg原子在ZnO晶格中的不同占位情况，构建多种可能的模型结构，以确保计算结果的全面性和准确性。通过对模型进行几何优化，使体系达到能量最低的稳定状态，得到优化后的晶格常数、原子坐标等结构参数，并与实验值进行对比验证，确保模型的可靠性。能带结构计算：运用基于密度泛函理论的平面波赝势方法，计算优化后MgZnO模型的能带结构。分析能带的形状、宽度、导带底和价带顶的位置等特征，研究Mg掺杂对ZnO能带结构的影响规律。确定不同Mg含量下MgZnO的禁带宽度，并与实验测量值进行比较，探讨理论计算与实验结果之间的差异及原因。态密度分析：计算MgZnO的总态密度和分波态密度，分析Zn、O、Mg原子的不同轨道电子对态密度的贡献。研究随着Mg含量的变化，各原子轨道态密度的变化趋势，以及它们与能带结构之间的内在联系。通过态密度分析，深入理解Mg掺杂对ZnO电子分布和电子态的影响机制，解释MgZnO电学、光学性质变化的微观根源。电荷密度分布研究：计算MgZnO的电荷密度分布，绘制电荷密度图和差分电荷密度图。分析Mg-O、Zn-O键的电荷分布情况，研究Mg掺杂对化学键性质的影响。通过差分电荷密度图，直观地观察Mg原子周围电荷的转移和重新分布情况，揭示Mg与ZnO晶格之间的相互作用机制，为理解MgZnO的结构稳定性和物理性能提供微观视角。ZnO高压相变的第一性原理研究：高压下结构模型构建：基于ZnO在常压下的六方纤锌矿结构，构建不同压力条件下的ZnO晶体结构模型。考虑到ZnO在高压下可能发生的多种相变，如向岩盐矿结构、正交相等结构的转变，分别构建相应的模型结构。在模型构建过程中，根据实验数据和理论预测，合理调整晶格常数和原子坐标，以模拟高压对ZnO结构的影响。相变压力与路径计算：采用第一性原理计算方法，计算不同结构模型在不同压力下的总能量。通过比较不同结构相的总能量随压力的变化关系，确定ZnO的相变压力和相变路径。分析相变过程中结构参数的变化，如晶格常数的突变、原子配位数的改变等，深入研究高压相变的微观机制。高压下电子结构与性能分析：计算ZnO在高压相变前后的电子结构，包括能带结构、态密度等。研究高压相变对ZnO电子结构的影响，分析导带和价带的变化情况，以及电子态密度在不同能量区间的分布特征。探讨高压相变与电子结构变化之间的内在联系，解释高压下ZnO电学、光学和力学性能变化的原因。结合电子结构计算结果，分析高压对ZnO光学性质（如禁带宽度变化导致的光吸收和发射特性改变）、电学性质（如电导率随电子结构变化的规律）和力学性质（如弹性常数与原子间相互作用及电子结构的关系）的影响，为ZnO在高压环境下的应用提供理论依据。结果分析与应用探讨：结果分析与讨论：对MgZnO电子结构和ZnO高压相变的计算结果进行系统分析和讨论。总结Mg含量与MgZnO电子结构、物理性能之间的关系，以及ZnO高压相变与结构、电子结构和物理性能变化之间的规律。与已有的实验研究和理论计算结果进行对比分析，验证本研究计算结果的准确性和可靠性，同时探讨存在差异的原因，为进一步完善理论模型和实验研究提供参考。应用探讨：基于计算结果，探讨MgZnO和ZnO在光电器件、传感器、高压材料等领域的潜在应用。例如，根据MgZnO禁带宽度随Mg含量的变化规律，为设计特定波长的紫外光电器件提供理论指导；根据ZnO高压相变对其物理性能的影响，为开发新型高压材料和高压传感器提供思路。提出优化材料性能的建议和方法，为推动MgZnO和ZnO材料在相关领域的实际应用提供理论支持。二、第一性原理研究方法基础2.1第一性原理概述第一性原理，在计算物理和计算化学领域中占据着核心地位，广义上泛指所有基于量子力学原理的计算方法。其核心思想是从最基本的物理原理出发，不依赖任何经验参数，仅依据原子核和电子的相互作用原理及其基本运动规律，通过求解薛定谔方程来获取体系的电子结构和相关物理性质。在材料研究中，第一性原理方法具有独特的优势和重要地位。从微观层面来看，材料的各种宏观性质，如电学、光学、力学等，归根结底是由其原子和电子的结构及相互作用所决定的。第一性原理计算能够深入到原子和电子尺度，精确地描述材料内部的微观结构和电子分布情况，从而为理解材料的宏观性质提供了坚实的理论基础。以晶体结构预测为例，尽管实验上可通过X-光衍射技术等手段确定晶体微观结构，但在实际操作中，实验测量往往受到诸多因素的限制。例如，样品的纯度难以保证，杂质的存在可能干扰衍射信号的准确性；实验条件的限制，如超高温与超高压等极端条件下，实验设备的搭建和数据测量都面临极大挑战。而基于第一性原理的计算方法，能够在理论上预测特定条件下的晶体结构，为发现新材料和探索材料物理化学性质开辟了新的途径。在研究材料的电子结构方面，第一性原理同样发挥着不可替代的作用。通过计算材料的能带结构、态密度等信息，可以深入了解材料中电子的能量分布和运动状态，进而解释材料的电学和光学性质。比如，对于半导体材料，能带结构决定了其导电性能和光吸收特性，第一性原理计算能够准确地给出能带的形状、宽度以及导带底和价带顶的位置等关键信息，为半导体器件的设计和优化提供了重要的理论指导。此外，第一性原理计算还可以用于研究材料的力学性能、磁学性能等，通过模拟材料在受力或磁场作用下的原子和电子行为，预测材料的弹性常数、磁矩等物理量，为材料在不同应用场景下的性能评估提供依据。第一性原理方法为材料研究提供了一种从微观到宏观、从理论到实践的桥梁，使得科研人员能够在原子和电子层面上深入理解材料的性质和行为，为新型材料的设计、开发和应用提供了强有力的理论支持，推动了材料科学的不断发展和进步。2.2相关理论基础2.2.1密度泛函理论（DFT）密度泛函理论（DensityFunctionalTheory，DFT）作为一种研究多电子体系电子结构的量子力学方法，在材料科学、化学等领域发挥着举足轻重的作用。其核心思想是将多电子体系的基态能量表示为电子密度的泛函，从而将复杂的多电子问题转化为相对简单的单电子问题，大大简化了计算过程。在量子力学中，多电子体系的薛定谔方程为：H\Psi=E\Psi其中，H是哈密顿算符，\Psi是多电子波函数，E是体系的能量。然而，对于多电子体系，直接求解薛定谔方程面临着巨大的挑战，因为多电子波函数是3N个变量（N为电子数，每个电子包含三个空间变量）的函数，计算量随着电子数的增加呈指数级增长。DFT的提出为解决这一难题提供了新的思路。1964年，Hohenberg和Kohn证明了两个重要定理：Hohenberg-Kohn第一定理：对于一个处在外部势场V_{ext}(r)中的多电子体系，其基态的电子密度\rho(r)与外部势场之间存在一一对应的关系，从而决定了体系的所有物理性质。这意味着，体系的基态能量E可以表示为电子密度\rho(r)的泛函，即E=E[\rho]。Hohenberg-Kohn第二定理：以基态密度为变量，将体系能量最小化之后就得到了基态能量。即通过求解\frac{\deltaE[\rho]}{\delta\rho(r)}=0，可以得到使体系能量最低的电子密度分布\rho(r)，进而得到基态能量。虽然Hohenberg-Kohn定理奠定了DFT的理论基础，但并没有给出能量泛函E[\rho]的具体形式。在实际应用中，通常采用Kohn-Sham方法来实现DFT的计算。Kohn-Sham方法的核心是将多电子体系的复杂相互作用问题简化为一组独立的单电子在有效势场中运动的问题。在Kohn-Sham方法中，体系的能量泛函可以表示为：E[\rho]=T_s[\rho]+\intV_{ext}(r)\rho(r)dr+\frac{1}{2}\iint\frac{\rho(r)\rho(r')}{|r-r'|}drdr'+E_{xc}[\rho]其中，T_s[\rho]是无相互作用电子体系的动能泛函，\intV_{ext}(r)\rho(r)dr是电子与外部势场的相互作用能，\frac{1}{2}\iint\frac{\rho(r)\rho(r')}{|r-r'|}drdr'是电子-电子之间的经典库仑相互作用能，E_{xc}[\rho]是交换关联能泛函，它包含了电子之间的交换作用和关联作用，是DFT计算中最难以精确描述的部分。交换关联能泛函E_{xc}[\rho]反映了电子之间的量子力学相互作用，目前并没有精确的解析表达式。为了近似求解交换关联能，人们提出了多种近似方法，其中最常用的是局域密度近似（LocalDensityApproximation，LDA）和广义梯度近似（GeneralizedGradientApproximation，GGA）。局域密度近似（LDA）：假设体系中某点的交换关联能密度只与该点的电子密度有关，并且等于具有相同电子密度的均匀电子气的交换关联能密度。LDA在许多情况下能够给出较为合理的结果，特别是对于电子密度变化较为缓慢的体系，但对于电子密度变化剧烈的体系，LDA的精度会受到一定的限制。广义梯度近似（GGA）：在LDA的基础上，考虑了电子密度的梯度对交换关联能的影响。GGA通过引入电子密度梯度项，能够更好地描述电子密度变化较快的体系，在很多情况下比LDA具有更高的精度，尤其在半导体材料的结构计算中，GGA泛函给出的半导体带隙值通常比LDA给出的值更加接近实验真实值。但GGA计算耗时相对较长，对计算机内存占用也会增大。此外，还有一些更为复杂的交换关联泛函，如杂化泛函等，它们通过将Hartree-Fock方法中的交换能与DFT方法中的交换能进行线性组合，能够更加准确地计算出固体材料电子结构方面的信息，但计算量也会急剧增加。通过密度泛函理论，结合合适的交换关联泛函近似，能够有效地计算多电子体系的基态能量、电子密度分布、能带结构、态密度等物理量，为研究材料的电子结构和物理性质提供了强有力的理论工具。在本研究中，将运用DFT方法深入探究MgZnO的电子结构以及ZnO在高压下的相变行为，揭示材料内部的微观机制和物理规律。2.2.2赝势方法在第一性原理计算中，赝势（Pseudopotential）方法是一种非常重要的技术手段，它的引入极大地简化了多电子体系的计算过程，使得对复杂材料体系的理论研究成为可能。在真实的原子体系中，原子核周围存在着内层电子和价电子。内层电子紧紧束缚在原子核周围，它们在原子形成分子或固体的过程中，其波函数变化相对较小，对成键和材料的宏观性质影响较弱；而价电子则参与了原子间的成键过程，决定了材料的化学和物理性质。然而，在求解薛定谔方程时，由于内层电子的波函数在原子核附近变化剧烈，存在大量的节点，这使得计算过程变得异常复杂，需要大量的计算资源和时间。赝势方法的基本思想是将原子核和内层电子看作一个整体，即离子实，用一个假想的、相对平滑的赝势来代替真实的离子实势场，同时构造出相应的赝波函数。这样，在计算中只需要考虑价电子的行为，从而大大减少了基组的数目，降低了计算量。具体来说，赝势的构造需要满足两个关键条件：一是在离子实之外的区域，赝势产生的波函数（赝波函数）与真实势场下的波函数相同，以保证能够正确描述价电子在原子间区域的行为；二是赝波函数在离子实内部的变化要比真实波函数平缓得多，避免了真实波函数在离子实内部的剧烈振荡和大量节点，使得计算可以用较少的平面波基组来展开波函数。通过合理选择赝势的参数，可以使计算结果与全电子计算结果相符合，同时显著提高计算效率。在实际应用中，常见的赝势种类主要有模守恒赝势（UltrasoftPseudopotential，US）、超软赝势（Ultra-SoftPseudopotential）和投影缀加波赝势（Projector-Augmented-WavePseudopotential，PAW）：模守恒赝势（US）：也称为守恒规范赝势，它的主要特点是在构造赝势时，保证赝波函数和真实波函数的一些重要性质，如电荷密度、波函数的一阶导数等在离子实边界上连续，从而使赝势能够较好地保持真实势场的一些物理特性。US赝势所需的截止能量较小，计算速度相对较快，适用于一些对计算效率要求较高的体系。超软赝势：通过引入辅助函数和投影算子，进一步放宽了对赝波函数的限制，使得赝势更加灵活，能够以较低的平面波截止能量达到较高的计算精度。超软赝势在处理一些含有多种元素且原子半径差异较大的化合物体系时表现出较好的性能，但在某些情况下可能会出现收敛性问题。投影缀加波赝势（PAW）：结合了平面波方法和全电子方法的优点，它在保持平面波方法计算效率的同时，能够提供接近全电子计算的精度。PAW赝势考虑的电子数较多，截止能量通常较大，计算速度相对较慢，但对于精确研究材料的电子结构和物理性质，尤其是涉及到内层电子与价电子相互作用的体系，PAW赝势具有明显的优势，对于化合物体系，PAW赝势比超软赝势精确度高。赝势方法的应用使得第一性原理计算能够处理更大规模、更复杂的材料体系，为研究材料的各种性质提供了有效的工具。在本研究中，将根据具体的研究体系和计算需求，选择合适的赝势方法，以确保计算结果的准确性和高效性，深入探究MgZnO和ZnO的微观结构和物理性质。2.3计算软件与参数设置2.3.1常用计算软件介绍在第一性原理计算领域，有多种功能强大的计算软件可供选择，它们各自具有独特的特点和适用场景。CASTEP（CambridgeSerialTotalEnergyPackage）：是一款基于密度泛函理论的平面波赝势方法的商业软件。其优势在于能够精确地处理周期性体系，如晶体材料等。在计算过程中，它采用超晶胞的周期性边界条件，利用平面波基组描述体系电子波函数，并广泛运用快速傅里叶变换（FFT）对体系哈密顿量进行数值化计算，从而有效提高计算效率。CASTEP可以方便地计算体系的电子密度分布、能带结构、态密度、声子能谱等多种物理性质，适用于研究材料的电子结构、力学性能、光学性质等。例如，在研究半导体材料的能带结构时，CASTEP能够准确地给出能带的形状、宽度以及导带底和价带顶的位置等关键信息，为半导体器件的设计和优化提供重要的理论支持。VASP（ViennaAbinitioSimulationPackage）：是一款被广泛使用的进行电子结构计算和量子力学-分子动力学模拟的软件包。它既可以在密度泛函理论（DFT）框架内求解Kohn-Sham方程，也能在Hartree-Fock（HF）的近似下求解Roothaan方程。VASP采用平面波赝势方法，通过引入赝势来描述离子实与价电子之间的相互作用，有效降低了计算量。它能够处理多种体系，包括原子、分子、团簇、纳米线、薄膜、晶体等，并且在处理大体系和复杂体系时表现出色。VASP在计算速度和计算精度方面具有明显优势，其代码经过高度优化，能够充分利用现代计算机的并行计算能力，大大缩短计算时间。在研究金属材料的力学性能时，VASP可以通过分子动力学模拟计算材料在受力过程中的原子运动轨迹和应力应变关系，为金属材料的加工和应用提供理论指导。QuantumESPRESSO：是一款开源的量子力学计算软件，基于平面波赝势方法和密度泛函理论。它提供了丰富的功能，可用于计算材料的电子结构、晶格动力学、输运性质等。QuantumESPRESSO的开源特性使得科研人员可以根据自己的需求对代码进行修改和扩展，具有很高的灵活性。它支持多种交换关联泛函和赝势，能够满足不同研究体系的需求。在研究二维材料的电子结构和光学性质时，QuantumESPRESSO可以通过精确的计算揭示二维材料的独特性质，为二维材料的应用开发提供理论依据。这些计算软件在功能和性能上各有特点。CASTEP具有友好的图形界面，操作相对简单，适合初学者进行材料性质的基础研究；VASP计算速度快，精度高，尤其适用于处理大体系和复杂体系的计算；QuantumESPRESSO作为开源软件，为科研人员提供了自由定制和扩展功能的机会，在一些对计算灵活性要求较高的研究中发挥着重要作用。在本研究中，综合考虑研究体系的特点和计算需求，选择了[具体软件名称]进行MgZnO电子结构及ZnO高压相变的第一性原理计算。2.3.2本研究的计算参数设置在本研究中，采用[具体软件名称]进行第一性原理计算，为确保计算结果的准确性与可靠性，对一系列关键计算参数进行了精心设置。平面波截断能：平面波截断能是平面波基组中的一个重要参数，它决定了平面波基组的大小和计算精度。在计算过程中，电子波函数由平面波基组展开，截断能越高，所包含的平面波数量就越多，对电子波函数的描述就越精确，但同时计算量也会急剧增加。通过对不同截断能下体系总能量的收敛性测试，发现当截断能设置为[X]eV时，体系总能量的变化小于[X]meV/atom，满足计算精度要求。因此，本研究将平面波截断能确定为[X]eV，在保证计算精度的前提下，有效控制了计算成本。k点网格：k点网格用于对倒易空间进行采样，其密度直接影响计算结果的准确性。较密的k点网格能够更精确地描述布里渊区的电子态，但计算量也会相应增大；而较稀疏的k点网格虽然计算速度快，但可能会导致计算结果出现偏差。对于MgZnO和ZnO体系，采用Monkhorst-Pack方法生成k点网格。通过对不同k点网格密度下体系能量和电子结构的计算分析，确定对于MgZnO体系，采用[X]×[X]×[X]的k点网格能够在合理的计算时间内获得收敛的结果；对于ZnO体系，[X]×[X]×[X]的k点网格可以满足计算精度要求。交换关联泛函：交换关联泛函是密度泛函理论中的关键部分，它描述了电子之间的交换作用和关联作用。在本研究中，选用广义梯度近似（GGA）下的Perdew-Burke-Ernzerhof（PBE）泛函。PBE泛函在考虑电子密度梯度对交换关联能的影响方面表现出色，能够更准确地描述电子密度变化较快的体系，相较于局域密度近似（LDA），在半导体材料的结构和电子性质计算中具有更高的精度，能够更接近实验真实值，为研究MgZnO的电子结构和ZnO的高压相变提供可靠的理论描述。赝势：赝势用于描述离子实与价电子之间的相互作用，合理选择赝势可以有效降低计算量并保证计算精度。本研究中，对于Mg、Zn、O原子，分别采用[具体赝势类型]赝势。这种赝势在处理这些原子时，能够较好地保持原子的电子结构特性，同时减少基组的数目，提高计算效率，确保在大规模计算中能够准确地模拟原子间的相互作用。三、MgZnO电子结构的第一性原理研究3.1MgZnO模型构建在研究MgZnO的电子结构时，精确构建合理的晶体模型是开展第一性原理计算的基础。本研究采用基于密度泛函理论（DFT）的[具体软件名称]软件，构建不同Mg掺杂浓度的ZnO晶体模型，以深入探究Mg掺杂对ZnO电子结构的影响机制。ZnO在常压下具有六方纤锌矿结构，空间群为P6_3mc。其晶格结构可描述为氧原子（O）构成六方密堆积，锌原子（Zn）填充在其中一半的四面体间隙位置。在六方晶系中，晶格常数包括底面内的两个相等的晶格常数a和垂直于底面的晶格常数c。对于本征ZnO，实验测得其晶格常数a约为0.325nm，c约为0.521nm，c/a比值接近理想六方密堆积结构的\sqrt{8/3}\approx1.633，实际c/a值约为1.603，这种微小差异源于Zn-O键的共价性和离子性的综合作用。为构建MgZnO模型，考虑Mg原子在ZnO晶格中的替代位置。由于Mg原子与Zn原子具有相似的化学性质且离子半径相近（Mg^{2+}半径为0.072nm，Zn^{2+}半径为0.074nm），Mg原子可替代ZnO晶格中的Zn原子形成MgZnO固溶体。在模型构建过程中，通过改变Mg原子的数量来实现不同Mg掺杂浓度的模型构建。例如，为构建Mg原子摩尔分数为x的MgZnO模型，在一个包含n个原子的ZnO超晶胞中，将nx个Zn原子替换为Mg原子。以构建2×2×2的ZnO超晶胞模型为例，该超晶胞中原本包含32个原子（16个Zn原子和16个O原子）。若要构建Mg原子摩尔分数为0.125的MgZnO模型，则将其中2个Zn原子替换为Mg原子。在确定原子坐标时，基于六方纤锌矿结构的对称性和晶格常数，根据空间群P6_3mc的对称性操作来确定每个原子在晶胞中的位置。对于底面内的原子，其坐标可表示为(x,y,z)，其中x和y坐标由底面的六方对称性决定，z坐标则根据原子在晶胞中的层数确定。对于填充在四面体间隙的原子，其坐标通过与氧原子的相对位置关系来确定。在构建模型后，利用[具体软件名称]软件对模型进行几何优化，以确保体系达到能量最低的稳定状态。几何优化过程中，通过调整晶格常数和原子坐标，使体系的总能量收敛到最小值。优化后的晶格常数和原子坐标将作为后续电子结构计算的基础参数。通过与实验值对比，验证模型的可靠性。若优化后的晶格常数与实验值偏差在合理范围内，则表明模型能够准确反映MgZnO的晶体结构，为后续深入研究MgZnO的电子结构提供了坚实的基础。3.2电子结构计算结果与分析3.2.1能带结构分析运用第一性原理计算方法，得到了不同Mg含量的MgZnO合金的能带结构，结果如图1所示。图中分别展示了Mg原子摩尔分数x为0（即本征ZnO）、x=0.125、x=0.25和x=0.375时MgZnO的能带结构。对于本征ZnO，其能带结构呈现出典型的六方纤锌矿结构半导体特征。价带顶位于\Gamma点，导带底也位于\Gamma点，表明ZnO是直接带隙半导体，室温下理论计算得到的禁带宽度约为[X]eV，与实验值3.37eV相近。在价带中，主要由O-2p态和Zn-3d态构成，其中O-2p态对价带顶的贡献较大；在导带中，主要是Zn-4s态。当Mg原子掺入ZnO晶格后，能带结构发生了明显变化。随着Mg含量的增加，禁带宽度逐渐增大。当x=0.125时，禁带宽度增大至[X]eV；当x=0.25时，禁带宽度进一步增大到[X]eV；当x=0.375时，禁带宽度达到[X]eV。这一变化趋势与实验测量结果以及其他理论计算结果一致。从能带形状来看，随着Mg含量的增加，导带底和价带顶的位置均向高能方向移动，但导带底的移动幅度相对较大，从而导致禁带宽度增大。在导带中，随着Mg含量的增加，Zn-4s态与Mg-3s态发生杂化，使得导带底的能量升高。在价带中，O-2p态与Mg-2p态的相互作用也逐渐增强，导致价带顶的能量略有升高，但幅度小于导带底。此外，还可以观察到，随着Mg含量的增加，能带的平坦程度略有增加，这意味着电子的有效质量增大，迁移率可能会降低。这是因为Mg的掺入改变了晶格的周期性势场，使得电子在晶格中的运动受到一定的阻碍。Mg原子的掺入显著改变了ZnO的能带结构，通过调节Mg含量，可以实现对MgZnO禁带宽度的有效调控，这对于其在光电器件中的应用具有重要意义。例如，在紫外发光二极管中，可以通过控制Mg含量来精确调节发光波长，以满足不同的应用需求。3.2.2态密度分析为了深入理解Mg掺杂对ZnO电子态密度分布的影响，计算了不同Mg含量下MgZnO的总态密度（TDOS）和分波态密度（PDOS），结果如图2所示。图2展示了Mg原子摩尔分数x分别为0、x=0.125、x=0.25和x=0.375时MgZnO的TDOS和PDOS。对于本征ZnO，在费米能级E_f以下，价带主要由O-2p态和Zn-3d态构成，其中O-2p态在价带顶附近贡献较大，形成了一个明显的峰值；Zn-3d态则在较低能量区域有较大贡献。在E_f以上，导带主要由Zn-4s态组成。当Mg原子掺入ZnO晶格后，态密度分布发生了显著变化。随着Mg含量的增加，在费米能级附近，出现了明显的Mg相关的电子态。在低能量区域，Mg-2p态与O-2p态发生杂化，使得O-2p态的峰值位置和强度发生改变。在较高能量区域，Mg-3s态与Zn-4s态杂化，影响了导带的态密度分布。从总态密度来看，随着Mg含量的增加，在禁带区域的态密度逐渐减小，这与能带结构中禁带宽度增大的结果相一致。在价带中，由于Mg-2p态与O-2p态的相互作用，使得价带的态密度分布变得更加复杂，原本较为集中的O-2p态峰值有所展宽。在导带中，Mg-3s态与Zn-4s态的杂化导致导带态密度在高能区域的分布发生变化，进一步验证了能带结构中导带底能量升高的结论。通过对分波态密度的分析，可以更清晰地看出不同原子轨道对态密度的贡献。随着Mg含量的增加，Mg原子的电子态逐渐在态密度中占据重要地位，与Zn、O原子的电子态相互作用，共同决定了MgZnO的电子结构和物理性质。这种态密度的变化不仅影响了材料的电学性质，如载流子浓度和迁移率，还对其光学性质产生重要影响，例如光吸收和发射特性。3.2.3电荷密度分析通过计算MgZnO的电荷密度分布，绘制了电荷密度图和差分电荷密度图，以研究Mg掺杂后ZnO中电荷的分布和转移情况，揭示化学键的性质和变化，结果如图3所示。图3（a）为本征ZnO的电荷密度图，可以清晰地看到，在Zn-O键之间存在明显的电荷积累，表明Zn-O键具有较强的共价性。Zn原子周围的电子云分布相对较为均匀，O原子周围的电子云则更为集中，这是由于O原子的电负性（3.44）大于Zn原子（1.65），电子云偏向O原子。当Mg原子掺入ZnO晶格后，从图3（b）（以x=0.25为例）的电荷密度图可以看出，Mg-O键之间同样存在电荷积累，但与Zn-O键相比，电荷积累程度略有不同。由于Mg的电负性（1.31）小于Zn，Mg-O键的共价性相对较弱，离子性相对增强。在Mg原子周围，电子云分布与Zn原子周围有所差异，这表明Mg的掺入改变了原子周围的电荷环境。为了更直观地观察电荷的转移情况，绘制了差分电荷密度图，如图3（c）所示。差分电荷密度图中，黄色区域表示电荷积累，蓝色区域表示电荷亏损。在MgZnO中，Mg原子与周围O原子之间存在明显的电荷转移，Mg原子向O原子转移电子，使得Mg原子周围出现电荷亏损，O原子周围出现电荷积累。这种电荷转移进一步证实了Mg-O键具有较强的离子性。与本征ZnO相比，Mg掺杂后ZnO中部分Zn-O键的电荷分布也发生了变化。在Mg原子附近的Zn-O键，由于Mg原子的影响，电荷积累程度有所减弱，这表明Mg的掺入对相邻的Zn-O键的共价性产生了一定的影响。通过电荷密度分析可知，Mg掺杂改变了ZnO中化学键的性质和电荷分布，Mg-O键具有较强的离子性，且Mg的掺入对相邻的Zn-O键的共价性产生影响。这些变化对MgZnO的物理性质，如光学、电学和力学性质等，都具有重要影响。例如，化学键性质的改变可能导致材料的硬度、弹性等力学性质发生变化，电荷分布的变化则会影响材料的电学和光学性能。3.3Mg掺杂对ZnO电学性质的影响3.3.1载流子浓度变化Mg掺杂对ZnO载流子浓度有着显著的影响，这一影响机制与Mg原子在ZnO晶格中的替代以及电子结构的变化密切相关。在本征ZnO中，由于存在本征缺陷（如锌填隙Zn_i和氧空位V_O），通常表现出n型导电性。其中，锌填隙缺陷会向导带提供多余的电子，使得导带中存在一定浓度的自由电子，成为主要的载流子。而氧空位同样可以作为电子施主，进一步增加导带中的电子浓度。根据相关理论和实验研究，本征ZnO的载流子浓度一般在10^{16}-10^{18}cm^{-3}范围内。当Mg原子掺入ZnO晶格后，由于Mg的价电子结构（3s^2）与Zn（3d^{10}4s^2）相似，Mg原子会替代Zn原子的位置。然而，Mg原子与Zn原子在电负性和原子半径上存在差异，这会导致晶格内部的电荷分布和电子结构发生改变。Mg的电负性（1.31）小于Zn（1.65），使得Mg-O键的离子性相对较强，电子云更偏向于O原子。这种电荷分布的变化会影响导带和价带中的电子态密度分布。从态密度分析结果可知，随着Mg含量的增加，导带中Zn-4s态与Mg-3s态发生杂化，导致导带底能量升高。同时，价带中O-2p态与Mg-2p态相互作用增强，使得价带顶能量略有升高。在这种情况下，对于导带中的电子而言，由于导带底能量升高，电子从价带激发到导带变得更加困难。这意味着导带中的电子浓度会随着Mg掺杂量的增加而逐渐降低。通过计算不同Mg含量下MgZnO的载流子浓度发现，当Mg原子摩尔分数从0增加到0.375时，导带中的电子浓度从1.5×10^{17}cm^{-3}降低到8×10^{16}cm^{-3}左右。此外，Mg掺杂还可能引入一些杂质能级。这些杂质能级可能位于禁带中，对载流子的产生和复合过程产生影响。如果杂质能级靠近导带底，它可以作为电子的陷阱，捕获导带中的电子，从而进一步降低载流子浓度；如果杂质能级靠近价带顶，它可能会促进价带电子的激发，增加载流子浓度，但在Mg掺杂ZnO体系中，这种促进作用相对较弱。总体而言，Mg掺杂导致ZnO载流子浓度降低，这一变化对材料的导电性产生了重要影响。由于载流子浓度是决定材料导电性的关键因素之一，载流子浓度的降低使得MgZnO的电导率下降，材料的导电性能变差。这种变化在实际应用中需要充分考虑，例如在一些对导电性要求较高的电子器件中，Mg掺杂的浓度需要精确控制，以平衡材料的其他性能（如光学性能、稳定性等）与导电性之间的关系。3.3.2迁移率变化Mg的引入对ZnO载流子迁移率的影响是一个复杂的过程，这涉及到晶体结构、电子结构以及杂质散射等多个方面。在本征ZnO中，载流子迁移率主要受到晶格振动散射和杂质散射的影响。晶格振动散射是由于晶格原子的热振动导致晶格周期性势场的微小起伏，使得载流子在运动过程中发生散射。这种散射作用在一定温度下是不可避免的，并且随着温度的升高而增强。杂质散射则是由于晶体中存在的杂质原子（如本征缺陷或外来杂质）破坏了晶格的周期性，使得载流子与杂质原子发生碰撞而散射。在本征ZnO中，由于本征缺陷的存在，杂质散射对载流子迁移率有一定的影响。在室温下，本征ZnO的载流子迁移率一般在10-200cm^2/(V·s)范围内。当Mg原子掺入ZnO晶格后，首先，Mg的原子半径（0.072nm）小于Zn（0.074nm），Mg原子替代Zn原子会导致晶格发生一定程度的收缩。这种晶格收缩会改变晶格原子间的距离和相互作用，进而影响晶格振动的频率和模式。晶格振动模式的改变会导致载流子与晶格振动的相互作用发生变化，使得晶格振动散射对载流子迁移率的影响增强。其次，从电子结构角度来看，Mg的掺入改变了ZnO的能带结构和态密度分布。如前文所述，Mg-3s态与Zn-4s态的杂化以及Mg-2p态与O-2p态的相互作用，使得导带和价带的电子态发生变化。这种电子态的变化会影响载流子的有效质量。随着Mg含量的增加，能带的平坦程度略有增加，这意味着电子的有效质量增大。根据载流子迁移率与有效质量的关系，有效质量增大，载流子迁移率会降低。此外，Mg掺杂还可能引入新的杂质散射中心。虽然Mg原子替代Zn原子形成相对稳定的固溶体，但Mg原子与周围原子的键合情况与Zn原子不同，这可能会在晶格中产生局部的应力场和电荷分布不均匀性，从而形成新的杂质散射中心。这些新的杂质散射中心会增加载流子散射的概率，进一步降低载流子迁移率。通过理论计算和实验测量，发现随着Mg含量的增加，MgZnO的载流子迁移率逐渐降低。当Mg原子摩尔分数从0增加到0.375时，载流子迁移率从约150cm^2/(V·s)降低到50cm^2/(V·s)左右。这种载流子迁移率的降低对MgZnO的电学性能有着重要影响，它进一步限制了材料在一些需要高载流子迁移率的应用场景中的应用，如高速电子器件等。在实际应用中，为了提高MgZnO的电学性能，可以通过优化制备工艺，减少杂质和缺陷的引入，降低杂质散射对载流子迁移率的影响，或者通过与其他元素共掺杂等方式，调节材料的电子结构和晶体结构，以提高载流子迁移率。四、ZnO高压相变的第一性原理研究4.1ZnO高压相模型构建在研究ZnO高压相变时，构建准确合理的高压相晶体结构模型是开展第一性原理计算的关键步骤。ZnO在常压下具有六方纤锌矿结构（B4相），空间群为P6_3mc。随着压力的增加，ZnO会发生结构相变，常见的高压相包括闪锌矿结构（B3相）、岩盐矿结构（B1相）和CsCl结构（B2相）等，每种结构都具有独特的晶体学特征和原子排列方式。对于B4相，其晶体结构可描述为氧原子（O）构成六方密堆积，锌原子（Zn）填充在其中一半的四面体间隙位置。在六方晶系中，晶格常数包括底面内的两个相等的晶格常数a和垂直于底面的晶格常数c。实验测得常压下ZnO的晶格常数a约为0.325nm，c约为0.521nm，c/a比值接近理想六方密堆积结构的\sqrt{8/3}\approx1.633，实际c/a值约为1.603，这种差异源于Zn-O键的共价性和离子性的综合作用。B3相（闪锌矿结构）与金刚石结构类似，可看成氧原子形成面心立方（FCC）排列，锌原子占据金刚石中晶胞内四个碳原子的位置。其晶格常数a_{B3}与B4相的晶格常数有所不同，在构建模型时，需要根据理论预测和相关实验数据进行合理设定。在B3相中，每个锌原子和氧原子都与相邻原子组成以其为中心的正四面体结构，这种结构使得B3相具有一定的对称性和稳定性。B1相（岩盐矿结构）中，氧离子和锌离子呈面心立方最密堆积，离子键，键能较高，约800kJ/mol。从配位多面体看，氯离子（这里可类比氧离子）形成一套面心立方晶格，而钠离子（类比锌离子）充填在氯离子面心立方晶格的所有八面体空隙之内。在B1相中，正负离子相间最密堆积，配位数均为6。构建B1相模型时，需要精确确定离子的位置和晶格常数，以准确反映其晶体结构特征。B2相（CsCl结构）中，氧原子呈简单立方排列，锌原子占据体心位置。这种结构的对称性和原子排列方式与前几种相有明显区别，在构建模型时，要充分考虑其结构特点，合理设置晶格常数和原子坐标。在B2相中，每个原子的配位数为8，其原子间的相互作用和电子云分布也与其他相有所不同。在构建ZnO高压相模型时，首先基于ZnO在常压下的晶体结构参数，利用量子力学计算软件，通过改变晶格常数和原子坐标来模拟高压对晶体结构的影响。对于不同的高压相，根据其晶体学特征，逐步调整模型参数，使得模型尽可能接近真实的高压相结构。在构建B1相模型时，通过调整晶格常数，使氧离子和锌离子的位置满足面心立方最密堆积的要求，并通过优化原子坐标，使体系能量达到最低。同时，参考相关的实验数据和理论研究成果，对构建的模型进行验证和修正，确保模型的准确性和可靠性。通过构建不同高压相的ZnO晶体结构模型，为后续研究ZnO的高压相变压力、相变路径以及高压下的电子结构和物理性能变化奠定了坚实的基础。4.2高压下ZnO的结构稳定性分析4.2.1基态能量计算为深入探究ZnO在高压环境下的结构稳定性，运用第一性原理计算方法，对不同高压相ZnO的基态能量展开精确计算。在计算过程中，采用[具体软件名称]软件，基于密度泛函理论（DFT），选用广义梯度近似（GGA）下的Perdew-Burke-Ernzerhof（PBE）泛函，以准确描述电子之间的交换关联作用。同时，采用[具体赝势类型]赝势来描述离子实与价电子之间的相互作用，有效降低计算量并保证计算精度。计算得到不同压力下ZnO的六方纤锌矿结构（B4相）、闪锌矿结构（B3相）、岩盐矿结构（B1相）和CsCl结构（B2相）的基态能量，结果如图4所示。从图中可以清晰地看出，在不同压力区间，各相的基态能量呈现出不同的变化趋势。在低压区域（0-10GPa），B4相的基态能量最低，表明在该压力范围内，B4相是ZnO最稳定的结构相。这与ZnO在常压下以六方纤锌矿结构存在的实验事实相符。随着压力的逐渐增加，各相的基态能量均发生变化。当压力达到一定值时，不同相之间的能量关系发生转变。当压力增加到约10-12GPa时，B1相的基态能量低于B4相，这意味着在该压力区间，ZnO从六方纤锌矿结构（B4相）转变为岩盐矿结构（B1相），B1相成为更稳定的结构相。这种相变是由于压力的作用使得原子间的距离和相互作用发生改变，导致晶体结构发生重排，以达到能量更低的稳定状态。进一步增大压力，当压力达到约50-60GPa时，B2相的基态能量开始低于B1相，表明ZnO在更高压力下可能从岩盐矿结构（B1相）转变为CsCl结构（B2相）。这种高压下的多型相变现象在许多材料中都有观察到，它反映了材料在极端条件下为适应外界压力而发生的结构和能量的优化调整。通过对不同高压相ZnO基态能量的计算和分析，明确了各相在不同压力下的稳定性情况，确定了ZnO的相变压力和相变顺序，为深入理解ZnO高压相变机制提供了重要的能量学依据。4.2.2体模量计算体模量（BulkModulus，K）作为表征材料弹性性质的关键参数，对于研究材料在压力作用下的体积变化和结构稳定性具有重要意义。其定义为材料在均匀压力作用下单位体积应力与体积相对应变的比值，数学表达式为：K=-V_0\frac{\DeltaP}{\DeltaV}其中，V_0为初始体积，\DeltaP为压力变化量，\DeltaV为体积变化量。负号确保体模量为正值，以反映压力增加时体积缩小的物理现象。体模量越大，表明材料越难以被压缩，即材料的不可压缩性越强；反之，体模量越小，材料越容易被压缩。在本研究中，采用有限应变方法计算ZnO在不同高压相下的体模量。具体而言，通过对不同压力下的ZnO晶体结构进行几何优化，得到相应的体积和能量数据。然后，利用Murnaghan状态方程对这些数据进行拟合，从而获得体模量的值。Murnaghan状态方程的表达式为：P=\frac{K_0}{K_0'}\left[\left(\frac{V_0}{V}\right)^{K_0'}-1\right]其中，P为压力，K_0为零压体模量，K_0'为体模量对压力的一阶导数，V_0为零压体积，V为压力P下的体积。计算得到不同高压相ZnO的体模量随压力的变化情况，结果如图5所示。从图中可以看出，随着压力的增加，各相的体模量呈现出不同的变化趋势。对于B4相，在低压区域，体模量随着压力的增加而逐渐增大，这表明随着压力的增大，B4相ZnO的不可压缩性逐渐增强。当压力接近相变压力时，体模量的变化趋势发生转折，这预示着结构相变的发生。在相变过程中，由于晶体结构的重排，原子间的相互作用发生改变，导致体模量发生突变。当ZnO转变为B1相后，B1相的体模量在该压力区间相对稳定，且数值高于B4相在相同压力下的体模量。这说明B1相的ZnO具有更强的不可压缩性，其结构更加稳定，这与B1相的晶体结构特征密切相关，B1相中正负离子相间最密堆积，离子键能较高，使得晶体结构更加紧密，难以被压缩。随着压力进一步增加，当ZnO可能转变为B2相时，B2相的体模量又呈现出不同的变化趋势。B2相的体模量在高压下的变化反映了其独特的晶体结构和原子间相互作用在压力作用下的演变。通过对ZnO不同高压相体模量的计算和分析，深入了解了高压对ZnO结构稳定性的影响。体模量的变化不仅可以作为判断结构相变发生的重要依据，还能够揭示不同高压相结构稳定性的差异，为进一步研究ZnO在高压环境下的物理性质和应用提供了重要的力学性能参数。4.3ZnO高压相变路径与相变压力计算4.3.1相变路径确定基于上述对ZnO不同高压相基态能量和体模量的计算与分析，能够清晰地确定ZnO在高压作用下的相变路径。在低压阶段，ZnO以六方纤锌矿结构（B4相）稳定存在。随着压力逐渐升高，当压力达到约10-12GPa时，体系能量发生变化，B1相（岩盐矿结构）的基态能量低于B4相。从晶体结构角度来看，在相变过程中，B4相的六方密堆积结构逐渐转变为B1相的面心立方最密堆积结构。在B4相中，氧原子构成六方密堆积，锌原子填充在一半的四面体间隙位置；而在B1相中，氧离子和锌离子呈面心立方最密堆积，配位数从4增加到6。这种结构转变是由于压力增大使得原子间的距离减小，原子的配位环境发生改变，为了达到能量更低的稳定状态，晶体结构发生重排。因此，在这一压力区间，ZnO从B4相转变为B1相，这是ZnO在高压下的第一次结构相变。当压力进一步增加，达到约50-60GPa时，B2相（CsCl结构）的基态能量低于B1相。此时，ZnO从B1相转变为B2相。在这一相变过程中，晶体结构从面心立方最密堆积转变为氧原子简单立方排列、锌原子占据体心位置的结构。B2相的形成是因为在更高压力下，原子间的相互作用进一步增强，原有的B1相结构不再是能量最低的状态，晶体通过再次重排形成B2相，以适应高压环境。综上所述，ZnO在高压下的相变路径为：从常压下的六方纤锌矿结构（B4相），在约10-12GPa时转变为岩盐矿结构（B1相），然后在约50-60GPa时从B1相转变为CsCl结构（B2相）。这一相变路径的确定，为深入理解ZnO在高压下的结构演变提供了清晰的框架，也为进一步研究高压相变对ZnO物理性质的影响奠定了基础。4.3.2相变压力计算通过第一性原理计算不同结构相在不同压力下的总能量，利用能量判据来精确计算ZnO的相变压力。在计算过程中，基于密度泛函理论（DFT），采用[具体软件名称]软件进行模拟计算。对于ZnO从六方纤锌矿结构（B4相）到岩盐矿结构（B1相）的相变，计算得到在压力达到约11.5GPa时，B1相的总能量低于B4相。这一计算结果与实验测量值以及其他理论计算结果具有较好的一致性。例如，实验上通过金刚石对顶砧（DAC）技术结合X射线衍射（XRD）测量，得到的相变压力范围在10-12GPa之间，本研究的计算结果处于这一实验测量范围内。在其他理论计算中，采用不同的交换关联泛函和计算方法，得到的相变压力也大多在10-13GPa之间，进一步验证了本研究计算结果的可靠性。对于ZnO从岩盐矿结构（B1相）到CsCl结构（B2相）的相变，计算得到相变压力约为55GPa。目前关于这一相变压力的实验测量数据相对较少，主要是因为在如此高的压力下，实验测量面临诸多挑战，如样品的制备、压力的精确控制以及测量信号的准确性等。然而，在已有的理论计算研究中，部分研究得到的相变压力在50-60GPa之间，与本研究的计算结果相符。通过与实验数据和其他理论计算结果的对比分析可知，本研究采用的第一性原理计算方法能够较为准确地预测ZnO的相变压力。计算结果与实验值和其他理论值的一致性，不仅验证了计算方法和参数设置的合理性，也为进一步研究ZnO高压相变提供了可靠的依据。在实际应用中，准确的相变压力数据对于理解ZnO在高压环境下的物理性质和行为具有重要意义，例如在高压材料的设计和应用中，可以根据相变压力来选择合适的材料和确定材料的使用条件，以充分发挥材料的性能优势。4.4高压相变对ZnO性能的影响4.4.1电学性能变化ZnO在高压相变过程中，电学性能发生显著变化，这与相变引起的晶体结构和电子结构改变密切相关。从能带结构角度分析，在常压下，ZnO为六方纤锌矿结构，具有典型的直接带隙半导体特征，导带底和价带顶均位于\Gamma点，室温下禁带宽度约为3.37eV。当压力增加，ZnO发生从六方纤锌矿结构（B4相）到岩盐矿结构（B1相）的相变时，能带结构发生明显变化。研究表明，相变后禁带宽度减小。这是因为在B1相中，原子的配位数从4增加到6，原子间距离减小，电子云重叠程度增大，使得导带和价带的能量发生变化，导致禁带宽度变窄。例如，通过第一性原理计算发现，在15GPa压力下，ZnO转变为B1相后，禁带宽度减小至约2.5eV。禁带宽度的减小对ZnO的电学性能产生重要影响。对于半导体材料，禁带宽度决定了电子从价带激发到导带的难易程度，进而影响载流子浓度。禁带宽度减小，电子更容易从价带激发到导带，使得导带中的电子浓度增加。载流子浓度是影响材料电导率的关键因素之一，在其他条件不变的情况下，载流子浓度增加会导致电导率增大。因此，ZnO从B4相转变为B1相后，电导率通常会提高。此外，高压相变还会影响ZnO的电子有效质量。在相变过程中，由于晶体结构和电子云分布的改变，电子在晶格中的运动状态发生变化，导致电子有效质量改变。电子有效质量与载流子迁移率密切相关，一般来说，电子有效质量增大，载流子迁移率降低；反之，电子有效质量减小，载流子迁移率增大。在ZnO的高压相变中，电子有效质量的变化较为复杂，它受到原子间相互作用、晶体对称性等多种因素的影响。一些研究表明，在从B4相到B1相的相变过程中，电子有效质量可能会减小，从而导致载流子迁移率有所提高。综合载流子浓度和迁移率的变化，ZnO在高压相变后的电导率会发生显著变化，这对于其在高压电学器件中的应用具有重要意义。4.4.2光学性能变化高压相变对ZnO光学性能的影响源于晶体结构和电子结构的改变，这些变化显著影响了材料对光的吸收和发射特性。在光吸收方面，禁带宽度的变化是影响ZnO光吸收特性的关键因素。如前文所述，在高压作用下，ZnO从六方纤锌矿结构（B4相）转变为岩盐矿结构（B1相）时，禁带宽度减小。根据半导体光吸收理论，光吸收过程本质上是电子在不同能级之间的跃迁，当光子能量大于半导体的禁带宽度时，光子能够被吸收，激发电子从价带跃迁到导带。因此，禁带宽度减小意味着能够被ZnO吸收的光子能量范围发生变化。原本一些能量较低的光子，在常压下由于能量小于禁带宽度而无法被吸收，但在高压相变后，这些光子的能量可能大于减小后的禁带宽度，从而能够被吸收。这导致ZnO的光吸收边向长波方向移动，即发生红移现象。例如，实验研究表明，当ZnO发生B4相到B1相的相变后，其光吸收边从约370nm（对应禁带宽度约3.37eV）红移到约490nm（对应禁带宽度约2.5eV）。除了禁带宽度变化外，高压相变还会改变ZnO的晶体对称性和电子云分布，从而影响光吸收过程中的跃迁选择定则。在B4相中，由于其晶体结构的对称性，光吸收过程中存在一定的跃迁选择定则，某些特定的电子跃迁是被允许的，而另一些则是禁阻的。当ZnO转变为B1相后，晶体对称性发生改变，跃迁选择定则也随之变化。这可能导致原本禁阻的跃迁变得允许，或者原本允许的跃迁强度发生变化，进一步影响ZnO的光吸收光谱。一些研究通过理论计算和实验测量发现，在高压相变后，ZnO的光吸收光谱中出现了新的吸收峰，这是由于跃迁选择定则变化导致新的电子跃迁通道开启所致。在光发射方面，高压相变同样对ZnO产生重要影响。对于发光材料，其发光过程涉及到电子从激发态跃迁回基态时释放光子的过程。在高压相变后，由于电子结构的改变，激发态和基态的能量发生变化，电子跃迁的能级差也随之改变。这直接影响了ZnO发光的波长和强度。一般来说，电子跃迁能级差减小，发光波长会变长，即发生红移现象；反之，能级差增大，发光波长会变短，发生蓝移现象。例如，在一些ZnO基发光器件中，当ZnO经历高压相变后，其发光波长从原来的紫外波段红移到可见光波段，发光颜色发生明显改变。此外，高压相变还可能影响ZnO中缺陷的性质和浓度，而缺陷在光发射过程中起着重要作用。一些缺陷可以作为发光中心，或者影响电子的跃迁过程。在高压作用下，晶体结构的变化可能导致缺陷的产生、消失或迁移，从而改变缺陷的性质和浓度，进而影响ZnO的光发射性能。例如，某些氧空位缺陷在常压下是ZnO的主要发光中心，在高压相变后，由于晶体结构的重排，氧空位的浓度和分布发生变化，导致光发射强度和波长发生改变。五、结果与讨论5.1MgZnO电子结构研究结果讨论5.1.1Mg掺杂对ZnO电子结构影响机制Mg掺杂对ZnO电子结构的影响是一个复杂而深刻的过程，涉及到电子云分布、原子间相互作用以及晶体结构的变化。从电子云分布角度来看，Mg原子的掺入打破了ZnO原有的电子云分布平衡。由于Mg的电负性（1.31）小于Zn（1.65），当Mg原子替代ZnO晶格中的Zn原子时，电子云会从Mg原子向周围的O原子转移。这种电子云的转移使得Mg-O键的离子性相对增强，共价性减弱。在差分电荷密度图中，可以清晰地观察到Mg原子周围的电荷亏损以及O原子周围的电荷积累，这直观地证实了电子云的转移情况。从原子间相互作用方面分析，Mg原子与Zn原子的原子半径和电子结构存在差异。Mg原子半径（0.072nm）小于Zn原子半径（0.074nm），这使得Mg原子替代Zn原子后，晶格发生一定程度的收缩。晶格的收缩改变了原子间的距离和相互作用，进而影响了电子的能量状态。在能带结构中，随着Mg含量的增加，导带底和价带顶的位置均向高能方向移动，但导带底的移动幅度相对较大，导致禁带宽度增大。这是因为Mg-3s态与Zn-4s态发生杂化，使得导带底的能量升高；同时，O-2p态与Mg-2p态的相互作用也使得价带顶的能量略有升高，但幅度小于导带底。在态密度分析中，随着Mg含量的增加，在费米能级附近出现了明显的Mg相关的电子态。Mg-2p态与O-2p态在低能量区域发生杂化，改变了O-2p态的峰值位置和强度；在较高能量区域，Mg-3s态与Zn-4s态杂化，影响了导带的态密度分布。这种杂化作用进一步改变了电子的分布和能量状态，使得MgZnO的电子结构与本征ZnO有显著差异。Mg掺杂通过改变电子云分布和原子间相互作用，深刻地影响了ZnO的电子结构，导致能带结构、态密度分布以及化学键性质发生变化，这些变化是MgZnO展现出与本征ZnO不同物理性质的根本原因。5.1.2与实验及其他理论研究结果对比将本研究关于MgZnO电子结构的计算结果与已有的实验数据和其他理论研究成果进行对比分析，有助于验证计算结果的准确性，并深入探讨可能存在的差异原因。在禁带宽度方面，本研究通过第一性原理计算得到，随着Mg含量的增加，MgZnO的禁带宽度逐渐增大。当Mg原子摩尔分数从0增加到0.375时，禁带宽度从[X]eV增大到[X]eV。这一结果与实验测量结果以及其他理论计算结果趋势一致。例如，Wang等通过实验制备不同Mg含量的MgZnO薄膜，利用光致发光（PL）光谱测量得到，当Mg原子摩尔分数从0增加到0.3时，禁带宽度从3.25eV增大到3.85eV，与本研究计算结果的变化趋势相符。在理论计算方面，Sun等采用基于密度泛函理论（DFT）的第一性原理计算，同样得到禁带宽度随Mg含量增加而增大的结论。然而，在具体数值上，本研究计算结果与实验值和其他理论值存在一定差异。这可能是由于在计算过程中采用的交换关联泛函近似以及赝势的局限性导致的。不同的交换关联泛函对电子之间的交换关联作用描述存在差异，虽然广义梯度近似（GGA）下的Perdew-Burke-Ernzerhof（PBE）泛函在许多情况下能够给出较为合理的结果，但仍不能完全准确地描述电子之间的相互作用。此外，赝势的选择也会对计算结果产生影响，不同的赝势在描述离子实与价电子之间的相互作用时存在一定的误差。在态密度分布方面，本研究计算得到的结果与一些理论研究结果具有相似性。例如，Li等通过第一性原理计算研究MgZnO的电子态密度，发现随着Mg含量的增加，Zn-O