基于第一性原理探究MgB₂掺杂系统的电子与光学特性

基于第一性原理探究MgB₂掺杂系统的电子与光学特性一、引言1.1MgB₂材料概述MgB₂作为一种备受瞩目的超导材料，自2001年被发现具有39K的超导转变温度以来，便在超导领域掀起了研究热潮。它的出现为超导材料的发展注入了新的活力，其独特的物理性质和潜在的应用价值吸引了众多科研人员的关注。MgB₂是一种金属间化合物，同时也具有离子化合物的特征，其镁层与硼层交替排列。在晶体结构方面，MgB₂属于六方晶系，空间群为P6/mmm，这种结构中，硼原子呈蜂巢状的平面排列，每一个硼原子都与周围三个硼原子相连，形成了二维的六角网状结构；而镁原子则位于硼原子层的上下两侧，呈六方密堆积排列，填充在硼原子层之间的空隙中，这种有序的原子排列方式赋予了MgB₂独特的物理性质。从电子结构角度来看，MgB₂的价电子主要由镁的3s电子和硼的2s、2p电子构成。硼原子的2p电子在形成共价键时，会与周围的硼原子相互作用，形成具有一定能量分布的能带结构。其中，硼原子的2p电子与相邻硼原子的2p电子之间的相互作用，对其超导特性有着至关重要的影响。通过对MgB₂电子结构的研究发现，其电子态密度在费米面附近存在明显的特征，这与超导电子配对机制密切相关。超导特性是MgB₂材料的核心优势之一，其超导转变温度（Tc）高达39K，这一温度相较于传统的低温超导材料有了显著提升，使其能够在液氮温区（77K）以下实现超导态。与高温超导体相比，MgB₂虽然转变温度相对较低，但其各向异性小，相干长度大，这使得它在制备约瑟夫森结器件等超导电子学器件时具有独特的优势。例如，在超导量子干涉器件（SQUID）中，MgB₂的这些特性有助于提高器件的灵敏度和稳定性，从而在微弱磁场探测等领域发挥重要作用。在超导高频器件方面，MgB₂的低损耗特性也为其在通信领域的应用提供了可能，有望实现高速、低功耗的信号传输。与低温超导体如NbTi和Nb₃Sn相比，MgB₂超导器件能够在液氮温区运行，这大大降低了制冷成本。液氮作为制冷剂，来源广泛且价格相对低廉，使得基于MgB₂的超导应用更加经济可行。从成本效益的角度来看，MgB₂在大规模应用中具有显著的优势，有望在电力传输、储能等领域实现商业化应用。在长距离输电线路中，利用MgB₂超导线材可以大幅降低输电损耗，提高能源传输效率，减少对环境的影响。在储能领域，MgB₂超导储能系统能够实现快速的能量存储和释放，为智能电网的稳定运行提供有力支持。MgB₂在超导领域具有重要的地位和广阔的应用潜力。它不仅为超导物理的研究提供了新的模型体系，有助于深入理解超导现象的本质和机制，推动超导理论的发展；而且在实际应用中，MgB₂在超导电子学、电力传输、医疗设备、交通运输等诸多领域都展现出了巨大的应用价值，有望成为下一代超导技术的核心材料，为解决能源、信息、医疗等领域的关键问题提供创新的解决方案。1.2MgB₂掺杂研究现状自2001年MgB₂的超导特性被发现以来，对其掺杂的研究便成为超导领域的一个重要方向。早期的研究主要聚焦于探索掺杂元素对MgB₂超导性能的影响，随着研究的深入，逐渐扩展到对其电子结构、晶体结构以及光学性质等多方面的研究。常见的掺杂元素种类繁多，可大致分为金属元素和非金属元素。金属元素如Al、Ti、Zr、Hf等，非金属元素如C、N、Si等。不同的掺杂元素因其原子结构和化学性质的差异，在MgB₂晶格中会产生不同的作用机制。根据掺杂元素在MgB₂晶格中的位置和作用，可将掺杂分为替代式掺杂和间隙式掺杂。替代式掺杂是指掺杂元素取代MgB₂晶格中的Mg或B原子，如C替代B原子；间隙式掺杂则是掺杂元素填充在MgB₂晶格的间隙位置，如某些金属原子填充在晶格间隙中。在电子结构影响方面，理论计算和实验研究表明，掺杂会改变MgB₂的电子态密度分布。当C掺杂时，C原子的外层电子结构与B原子不同，C替代B后会使B原子周围的电子云分布发生变化，进而影响费米面附近的电子态密度。这种变化会改变电子之间的相互作用，对超导电子配对产生影响，最终影响超导转变温度。一些金属元素的掺杂也会导致电子结构的变化，如Al掺杂会引入额外的电子态，改变费米面的形状和电子的散射机制。在超导电性方面，掺杂对MgB₂的超导转变温度（Tc）、临界电流密度（Jc）和上临界场（Hc2）等关键超导性能参数有显著影响。许多研究表明，适量的C掺杂可以提高MgB₂在高磁场下的临界电流密度，这是因为C掺杂引入了有效的磁通钉扎中心，增强了对磁通线的钉扎作用，从而抑制了磁通线的运动，提高了超导载流能力。然而，C掺杂通常会导致超导转变温度有所下降，这是由于C替代B后，改变了B-B键的强度和电子结构，影响了超导电子配对机制。不同的掺杂元素对超导性能的影响存在差异，一些金属元素的掺杂可能会同时提高临界电流密度和上临界场，但对超导转变温度的影响则因元素而异。在光学性质方面，MgB₂掺杂系统的研究相对较少，但近年来也逐渐受到关注。掺杂会改变MgB₂的能带结构，进而影响其光学吸收、反射和发射等性质。理论计算预测，某些掺杂元素可能会在MgB₂的能带结构中引入新的能级，这些能级会参与光的吸收和发射过程，导致光学吸收边的移动和新的发光峰的出现。在实验研究中，通过光谱分析等手段也证实了掺杂对MgB₂光学性质的影响。例如，通过测量掺杂MgB₂样品的紫外-可见吸收光谱，发现掺杂后吸收光谱的形状和吸收峰位置发生了变化，这与理论计算中能带结构的变化相呼应。对掺杂MgB₂的红外光学性质研究也发现，掺杂会改变其晶格振动模式，从而影响红外吸收特性。尽管目前对MgB₂掺杂系统的研究取得了一定的进展，但仍存在许多问题和挑战。不同掺杂元素之间的协同作用机制尚不完全清楚，如何通过多元素共掺实现对MgB₂性能的优化还需要进一步探索。在实际应用中，如何制备高质量、均匀掺杂的MgB₂材料，以及如何解决掺杂带来的材料加工性能下降等问题，也是亟待解决的关键问题。未来，随着研究的不断深入和技术的不断进步，有望进一步揭示MgB₂掺杂系统的内在物理机制，为其在超导电子学、能源等领域的应用提供更坚实的理论基础和技术支持。1.3研究目的与意义本研究旨在深入探究MgB₂掺杂系统的电子结构和光学性质，通过系统的理论计算和分析，揭示掺杂对MgB₂微观结构和宏观性能的影响机制，为MgB₂超导材料的性能优化和应用拓展提供坚实的理论基础。从科学研究的角度来看，深入理解MgB₂掺杂系统的电子结构和光学性质，有助于揭示超导现象的微观物理机制。超导现象的微观机制一直是凝聚态物理领域的研究热点和难点，虽然BCS理论能够解释常规超导体的超导机制，但对于MgB₂这种具有独特晶体结构和电子特性的超导材料，其超导机制仍存在许多未解之谜。通过研究掺杂对MgB₂电子结构的影响，如电子态密度分布、能带结构变化以及电子-声子相互作用的改变等，可以进一步揭示其超导电子配对的微观过程，为完善超导理论提供重要的实验和理论依据。对掺杂MgB₂光学性质的研究，能够为理解材料的光与物质相互作用提供新的视角。光学性质是材料电子结构的宏观体现，通过研究掺杂导致的光学吸收、反射、发射等性质的变化，可以反推材料电子结构的改变，为研究电子结构提供一种间接但有效的手段。这种研究方法有助于深入理解光在超导材料中的传播和相互作用规律，丰富和拓展光与物质相互作用的理论体系。在应用层面，研究MgB₂掺杂系统具有至关重要的意义。MgB₂作为一种有潜力的超导材料，其在超导电子学、能源、医疗等领域展现出广阔的应用前景。然而，目前MgB₂材料的实际应用仍面临一些挑战，如临界电流密度较低、上临界场不够高等问题，这些问题限制了其在强磁场和高电流密度环境下的应用。通过掺杂来优化MgB₂的性能是解决这些问题的关键途径之一。深入了解掺杂对电子结构和光学性质的影响，能够为开发新型超导材料提供理论指导。例如，通过合理选择掺杂元素和控制掺杂浓度，可以调控MgB₂的电子结构，增强电子-声子相互作用，从而提高超导转变温度和临界电流密度；也可以引入有效的磁通钉扎中心，改善材料在磁场中的性能。在光学性质方面，对掺杂MgB₂光学性质的研究成果，为其在光电器件中的应用提供了可能。例如，利用掺杂MgB₂独特的光学性质，可以开发新型的超导光探测器、发光二极管等光电器件，拓展其在光通信、光学传感等领域的应用。二、理论基础与研究方法2.1第一性原理第一性原理，又被称作从头算，是基于量子力学原理发展而来的重要理论计算方法。其核心思想在于将由多个原子构成的体系视为由多个电子和原子核组成的复杂系统，从微观层面出发，依据量子力学的基本原理对该系统进行处理，并且在处理过程中尽量减少对经验参数的依赖，最大程度地实现“非经验性”计算。在量子力学的框架下，微观粒子的运动遵循薛定谔方程。对于一个包含多个电子和原子核的体系，其薛定谔方程可以表示为：\hat{H}\Psi=E\Psi其中，\hat{H}是哈密顿算符，它包含了体系中所有粒子的动能以及粒子之间的相互作用势能；\Psi是体系的波函数，描述了体系中粒子的量子状态；E则是体系的能量本征值。然而，直接求解多电子体系的薛定谔方程是极其困难的，因为电子之间存在着复杂的相互作用，这使得方程的求解面临着巨大的计算量和数学上的复杂性。为了简化计算，第一性原理引入了一系列近似处理方法。其中，最常用的近似是Born-Oppenheimer近似，该近似基于原子核质量远大于电子质量这一事实，将电子的运动和原子核的运动分开考虑。在Born-Oppenheimer近似下，原子核可以看作是固定在空间中的点电荷，电子则在这些固定的原子核产生的势场中运动。这样，就可以将多电子体系的薛定谔方程简化为只包含电子的薛定谔方程，大大降低了计算的难度。密度泛函理论（DensityFunctionalTheory，DFT）是第一性原理计算中广泛应用的一种理论方法。其核心概念是将体系的基态能量表示为电子密度的泛函，即体系的所有性质都可以通过电子密度来描述。在DFT中，引入了交换关联泛函来描述电子之间的交换作用和关联作用，从而将多电子问题转化为单电子问题进行求解。Kohn-Sham方程是DFT的核心方程，通过求解该方程可以得到体系的电子密度和能量。在材料研究领域，第一性原理具有诸多显著的优势。第一性原理计算能够从原子和电子层面深入揭示材料的微观结构与宏观性能之间的内在联系。通过计算材料的电子结构，如能带结构、态密度等，可以了解电子在材料中的分布和运动状态，进而解释材料的电学、光学、磁学等性质。在研究半导体材料时，通过第一性原理计算能带结构，可以准确确定材料的禁带宽度，这对于理解半导体的导电特性以及设计半导体器件具有重要的指导意义。在研究磁性材料时，计算电子的自旋态和磁矩分布，可以深入探究材料的磁性起源和磁相互作用机制。第一性原理计算还能够预测材料的新性质和新现象，为新材料的研发提供理论依据。在探索新型超导材料时，通过第一性原理计算不同元素组合和晶体结构下材料的电子结构和超导特性，可以筛选出具有潜在超导性能的材料体系，大大缩短了新材料研发的周期和成本。在设计新型催化剂时，计算材料表面与反应物分子之间的相互作用能和反应路径，可以预测材料的催化活性和选择性，指导新型高效催化剂的设计。第一性原理计算在材料研究中的应用范围极为广泛，涵盖了材料科学的各个领域。在晶体结构预测方面，第一性原理可以根据材料的化学成分和原子间相互作用，预测材料可能存在的晶体结构，并通过计算能量来确定最稳定的结构。在研究金属合金时，通过第一性原理计算不同原子排列方式下合金的能量和结构稳定性，可以预测合金的相图和析出相，为合金的成分设计和热处理工艺优化提供指导。在研究陶瓷材料时，计算晶体结构和化学键特性，可以深入了解陶瓷材料的硬度、脆性等力学性能的微观本质。在研究半导体材料时，计算能带结构和载流子迁移率等电学性质，可以为半导体器件的设计和性能优化提供理论支持。在研究超导材料时，计算电子-声子相互作用和超导能隙等超导特性，可以深入理解超导机制，为新型超导材料的研发提供理论指导。在研究磁性材料时，计算磁矩和磁相互作用等磁性性质，可以为磁性存储和磁性传感器等器件的设计提供依据。在研究光学材料时，计算介电函数、吸收系数等光学性质，可以为光电器件的设计和光学材料的应用提供理论支持。在研究纳米材料时，由于纳米材料具有尺寸效应和表面效应等特殊性质，第一性原理计算可以从原子尺度研究纳米材料的结构和性能，为纳米材料的制备和应用提供理论指导。在研究复合材料时，计算不同相之间的界面相互作用和应力分布等力学性能，可以为复合材料的设计和优化提供理论依据。在研究生物材料时，计算生物分子与材料表面的相互作用和生物相容性等性质，可以为生物材料的设计和应用提供理论支持。2.2密度泛函理论密度泛函理论（DensityFunctionalTheory，DFT）是一种用于研究多电子体系电子结构的量子力学方法，在凝聚态物理、计算材料学和计算化学等领域有着广泛的应用。其核心在于将体系的基态能量表示为电子密度的泛函，即体系的所有性质都可通过电子密度来描述，这一理论的提出极大地简化了多电子体系的计算问题。Hohenberg-Kohn定理是密度泛函理论的重要基石，它包含两个重要内容。第一定理指出，对于一个处在外部势场V_{ext}(\vec{r})中的多电子体系，其基态电子密度\rho(\vec{r})是唯一确定的，并且体系的基态能量是电子密度的泛函。这意味着，只要确定了基态电子密度，体系的所有基态性质都可以通过该电子密度的泛函来确定。从数学角度来看，对于给定的外部势场V_{ext}(\vec{r})，体系的哈密顿量\hat{H}可以表示为：\hat{H}=\sum_{i=1}^{N}\left(-\frac{\hbar^{2}}{2m}\nabla_{i}^{2}+V_{ext}(\vec{r}_{i})\right)+\frac{1}{2}\sum_{i\neqj}\frac{e^{2}}{\vert\vec{r}_{i}-\vec{r}_{j}\vert}其中，N为电子数，\vec{r}_{i}为第i个电子的位置，m为电子质量，\hbar为约化普朗克常数，e为电子电荷。根据第一定理，存在一个泛函E[\rho]，使得体系的基态能量E_{0}满足：E_{0}=\min_{\rho}E[\rho]并且，当\rho取基态电子密度\rho_{0}(\vec{r})时，E[\rho]取得最小值E_{0}。第二定理表明，可以通过变分原理来求解体系的基态电子密度和基态能量。即对于一个试探电子密度\rho(\vec{r})，通过计算能量泛函E[\rho]，并不断调整\rho(\vec{r})，使其满足能量最小化条件，最终得到的电子密度就是体系的基态电子密度，对应的能量就是基态能量。为了具体求解能量泛函，Kohn和Sham提出了Kohn-Sham方程。Kohn-Sham方法的核心思想是将多电子体系的复杂相互作用问题转化为一组单电子在有效势场中运动的问题。在Kohn-Sham框架下，体系的能量泛函E[\rho]可以表示为：E[\rho]=T_{s}[\rho]+E_{H}[\rho]+E_{xc}[\rho]+\intV_{ext}(\vec{r})\rho(\vec{r})d\vec{r}其中，T_{s}[\rho]是无相互作用电子体系的动能泛函，E_{H}[\rho]是电子-电子库仑相互作用能，即Hartree能，E_{xc}[\rho]是交换关联能泛函，\intV_{ext}(\vec{r})\rho(\vec{r})d\vec{r}是电子与外部势场的相互作用能。无相互作用电子体系的动能泛函T_{s}[\rho]可通过求解单电子的Kohn-Sham轨道\psi_{i}(\vec{r})得到：T_{s}[\rho]=-\frac{\hbar^{2}}{2m}\sum_{i=1}^{N}\int\vert\nabla\psi_{i}(\vec{r})\vert^{2}d\vec{r}Hartree能E_{H}[\rho]描述了电子之间的经典库仑相互作用，其表达式为：E_{H}[\rho]=\frac{1}{2}\iint\frac{\rho(\vec{r})\rho(\vec{r}')}{\vert\vec{r}-\vec{r}'\vert}d\vec{r}d\vec{r}'交换关联能泛函E_{xc}[\rho]则包含了电子之间的交换作用和关联作用，由于交换关联能的精确形式难以确定，通常采用近似方法来处理。在实际计算中，常用的近似方法有局域密度近似（LocalDensityApproximation，LDA）和广义梯度近似（GeneralizedGradientApproximation，GGA）。局域密度近似假设体系中某点的交换关联能密度只与该点的电子密度有关，其表达式为：E_{xc}^{LDA}[\rho]=\int\rho(\vec{r})\epsilon_{xc}^{LDA}(\rho(\vec{r}))d\vec{r}其中，\epsilon_{xc}^{LDA}(\rho)是均匀电子气的交换关联能密度，可通过理论计算或实验数据拟合得到。LDA在处理一些简单体系时取得了较好的结果，但对于电子密度变化较大的体系，其精度会受到一定限制。广义梯度近似则考虑了电子密度的梯度对交换关联能的影响，其交换关联能泛函不仅依赖于电子密度\rho(\vec{r})，还依赖于电子密度的梯度\nabla\rho(\vec{r})，表达式为：E_{xc}^{GGA}[\rho]=\int\rho(\vec{r})\epsilon_{xc}^{GGA}(\rho(\vec{r}),\nabla\rho(\vec{r}))d\vec{r}GGA在一定程度上改进了LDA的不足，对于许多体系的计算结果更加准确，尤其是对于表面、界面等电子密度变化较大的体系。GGA有多种具体的形式，如Perdew-Burke-Ernzerhof（PBE）泛函等，不同的GGA泛函在不同的体系中表现出不同的性能。2.3计算程序与方法本研究使用MaterialsStudio软件中的CASTEP模块进行第一性原理计算。CASTEP是一款基于密度泛函理论的量子力学程序包，它能够对晶体、表面和分子等体系进行精确的电子结构计算和性质预测。该程序包具有强大的功能和友好的操作界面，能够方便地实现模型搭建、参数设置以及结果分析等一系列计算流程。在模型构建方面，首先需要获取MgB₂的晶体结构信息。MgB₂属于六方晶系，空间群为P6/mmm，晶格常数a=3.086\mathring{A}，c=3.524\mathring{A}。利用MaterialsStudio软件中的晶体结构构建工具，按照MgB₂的晶体结构参数搭建其原胞模型。对于掺杂体系，根据不同的掺杂类型（替代式掺杂或间隙式掺杂），在原胞模型中相应位置替换或添加掺杂原子，构建出掺杂MgB₂的超晶胞模型。为了确保计算结果的准确性，超晶胞的大小需要合理选择，既要包含足够多的原子以体现掺杂的影响，又要控制计算量在可承受范围内。经过测试，选择了2\times2\times2的超晶胞，其中包含32个原子（16个Mg原子和16个B原子），对于掺杂体系，则根据掺杂浓度相应调整原子组成。在参数设置方面，采用广义梯度近似（GGA）中的Perdew-Burke-Ernzerhof（PBE）泛函来描述电子之间的交换关联作用。PBE泛函在处理固体材料的电子结构和性质时表现出较好的准确性和可靠性，能够较为准确地描述MgB₂及其掺杂体系的电子相互作用。赝势选用超软赝势（USPP），超软赝势能够在保证计算精度的前提下，有效降低计算量，提高计算效率。它通过对离子实与价电子之间的相互作用进行近似处理，使得在平面波基组下的计算更加高效。平面波截断能设置为400eV，这个截断能值是通过对不同截断能下的计算结果进行收敛性测试确定的。在收敛性测试中，逐步增加平面波截断能，观察体系总能量、原子受力等物理量的变化情况，当截断能达到400eV时，这些物理量的变化趋于稳定，表明此时的计算结果已经收敛，能够满足计算精度要求。布里渊区积分采用Monkhorst-Pack方法，k点网格设置为6\times6\times4。k点网格的选取直接影响到计算结果的精度和计算量，通过对不同k点网格下的计算结果进行对比分析，发现6\times6\times4的k点网格能够在保证计算精度的同时，兼顾计算效率。在计算过程中，使用BFGS算法进行几何优化，BFGS算法是一种拟牛顿算法，它通过迭代更新海森矩阵的近似逆来搜索能量极小值，具有收敛速度快、稳定性好的优点。在几何优化过程中，设置原子受力收敛标准为0.01eV/Å，能量收敛标准为1\times10^{-6}eV/atom，最大位移收敛标准为5\times10^{-4}Å。当体系的原子受力、能量和最大位移满足这些收敛标准时，认为几何优化达到收敛，得到的结构为体系的稳定结构。三、Al、Li掺杂MgB₂系统研究3.1理论模型与计算设置为深入探究Al、Li掺杂对MgB₂系统电子结构和光学性质的影响，构建合理准确的理论模型是首要任务。在本研究中，基于MgB₂的六方晶系结构（空间群为P6/mmm），采用超晶胞模型来模拟掺杂体系。对于Al、Li掺杂MgB₂的超晶胞构建，选用了2\times2\times2的超晶胞，此超晶胞中包含32个原子，其中16个Mg原子和16个B原子。在构建Al掺杂模型时，考虑到Al原子与B原子的原子半径和化学性质的相似性，将Al原子替代B原子的位置进行掺杂。具体而言，在超晶胞中随机选取一定比例的B原子，用Al原子进行替换，以模拟不同掺杂浓度下的体系。在研究5%Al掺杂时，在超晶胞中选取1个B原子，用Al原子替换，从而实现对该掺杂浓度体系的模拟研究。在Li掺杂模型构建中，由于Li原子半径较小，更倾向于填充在晶格间隙位置，因此将Li原子放置在MgB₂晶格的间隙位置。通过精确计算和分析，确定Li原子在晶格间隙中的稳定位置，以确保模型的准确性和可靠性。经过结构优化计算，发现Li原子填充在硼原子层之间的间隙位置时，体系能量最低，结构最为稳定。在确定原子位置后，对晶胞参数进行细致的优化计算。通过第一性原理计算中的几何优化步骤，不断调整晶胞的棱长a和c以及棱间夹角，使得体系的总能量达到最低，原子受力收敛到极小值，从而得到稳定的晶胞参数。在未掺杂的MgB₂超晶胞中，经过优化计算得到的晶格常数a=3.090\mathring{A}，c=3.528\mathring{A}，与实验值和其他理论计算结果相符。对于Al、Li掺杂体系，随着掺杂原子的引入，晶胞参数会发生相应的变化。在5%Al掺杂体系中，优化后的晶格常数a略微增大至3.095\mathring{A}，c变化较小为3.529\mathring{A}；在3%Li掺杂体系中，a减小为3.084\mathring{A}，c增大为3.532\mathring{A}。这些晶胞参数的变化反映了掺杂原子与MgB₂晶格之间的相互作用，对后续的电子结构和光学性质研究具有重要影响。本研究采用基于密度泛函理论（DFT）的CASTEP模块进行计算。在计算过程中，交换关联泛函选择广义梯度近似（GGA）中的Perdew-Burke-Ernzerhof（PBE）泛函，该泛函在处理固体材料的电子结构和性质时能够提供较为准确的结果。赝势选用超软赝势（USPP），它可以有效降低计算量，提高计算效率，同时保证计算精度。平面波截断能设置为400eV，此截断能值是通过对不同截断能下的计算结果进行收敛性测试确定的。当截断能达到400eV时，体系的总能量、原子受力等物理量的变化趋于稳定，满足计算精度要求。布里渊区积分采用Monkhorst-Pack方法，k点网格设置为6\times6\times4。k点网格的选取对计算结果的精度和计算量有显著影响，通过对不同k点网格下的计算结果进行对比分析，发现6\times6\times4的k点网格能够在保证计算精度的同时，兼顾计算效率。在几何优化过程中，使用BFGS算法进行结构优化，设置原子受力收敛标准为0.01eV/Å，能量收敛标准为1\times10^{-6}eV/atom，最大位移收敛标准为5\times10^{-4}Å。当体系的原子受力、能量和最大位移满足这些收敛标准时，认为几何优化达到收敛，得到的结构为体系的稳定结构。在电子结构计算中，采用自洽场（SCF）迭代方法，通过不断迭代计算，使体系的电子密度和能量达到自洽，从而得到准确的电子结构信息。3.2计算结果与分析3.2.1几何优化结果经过几何优化后，得到了未掺杂MgB₂以及Al、Li掺杂MgB₂体系的稳定晶胞结构，其晶胞结构如图1所示，（a）为未掺杂MgB₂，（b）为5%Al掺杂MgB₂，（c）为3%Li掺杂MgB₂。在未掺杂的MgB₂中，硼原子呈蜂巢状的平面排列，每一个硼原子都与周围三个硼原子相连，形成二维的六角网状结构；镁原子位于硼原子层的上下两侧，呈六方密堆积排列，填充在硼原子层之间的空隙中，这种结构与MgB₂的理论晶体结构一致。图1MgB₂及掺杂体系晶胞结构对于5%Al掺杂的MgB₂体系，Al原子替代了部分B原子的位置。从优化后的结构可以看出，Al原子的掺入导致周围B原子的位置发生了一定程度的调整，以适应新的原子间相互作用。这种调整主要表现为B原子与Al原子之间的距离缩短，而与其他B原子的距离略有变化。Al-B键长为1.61Å，相较于未掺杂时B-B键长1.70Å有所缩短，这是由于Al原子的电负性与B原子不同，Al-B之间的化学键具有更强的极性，导致原子间的吸引力增强，从而使键长缩短。Al原子的掺入也使得晶胞参数发生了变化，晶格常数a略微增大至3.095\mathring{A}，c变化较小为3.529\mathring{A}。a的增大可能是由于Al原子的原子半径略大于B原子，替代B原子后占据了更大的空间，从而导致晶胞在a方向上的尺寸增加；而c方向上的变化较小，说明Al原子的掺入对该方向上的原子堆积方式影响相对较小。在3%Li掺杂的MgB₂体系中，Li原子填充在硼原子层之间的间隙位置。优化后的结构显示，Li原子的存在使得周围的Mg原子和B原子的位置也发生了一定的位移。Li原子与周围B原子之间形成了一定的相互作用，Li-B距离为2.05Å。这种相互作用导致B原子向Li原子靠近，从而改变了B原子层的平面结构，使其出现了一定程度的扭曲。Li原子的掺入对晶胞参数的影响与Al掺杂不同，晶格常数a减小为3.084\mathring{A}，c增大为3.532\mathring{A}。a的减小可能是因为Li原子填充在间隙位置后，对a方向上的原子间作用力产生了影响，使得原子间的距离缩短；而c的增大则可能是由于Li原子与硼原子层之间的相互作用，使得硼原子层之间的距离增大，进而导致c方向上的尺寸增加。通过分析掺杂前后晶胞参数的变化以及原子间距离的改变，可以探讨掺杂对MgB₂晶体结构稳定性的影响。一般来说，晶体结构的稳定性与体系的能量密切相关，能量越低，结构越稳定。在Al掺杂体系中，虽然晶胞参数发生了变化，但体系的总能量相较于未掺杂时略有降低，这表明Al掺杂在一定程度上增强了晶体结构的稳定性。这可能是因为Al-B键的形成增强了原子间的相互作用，使得晶体结构更加紧密。而在Li掺杂体系中，体系的总能量略有升高，说明Li掺杂可能会降低晶体结构的稳定性。这可能是由于Li原子填充在间隙位置后，破坏了原有的原子堆积方式，导致原子间的相互作用发生了变化，从而使体系的能量升高。3.2.2掺杂与键布居分析键布居分析是研究原子间化学键性质的重要方法，它通过计算原子间的电子分布情况，来揭示化学键的强度和性质。在本研究中，对Al、Li掺杂后MgB₂系统中各原子间的键布居进行了计算，计算结果如表1所示，其中列出了未掺杂MgB₂以及5%Al掺杂、3%Li掺杂MgB₂体系中B-B、B-Mg等键的布居数。表1MgB₂及掺杂体系中键布居数体系B-B键布居数B-Mg键布居数Al-B键布居数（Al掺杂）Li-B键布居数（Li掺杂）未掺杂MgB₂0.350.20--5%Al掺杂MgB₂0.320.180.38-3%Li掺杂MgB₂0.330.19-0.25在未掺杂的MgB₂中，B-B键的布居数为0.35，B-Mg键的布居数为0.20。B-B键具有较高的布居数，表明B-B之间存在较强的共价键相互作用，这与硼原子的成键特性相符，硼原子通过共享电子形成稳定的共价键网络，对MgB₂的结构和性质起着重要的支撑作用。B-Mg键的布居数相对较低，说明B-Mg之间的化学键具有一定的离子性，镁原子向硼原子提供电子，形成离子键相互作用，这种离子键与共价键的共存，赋予了MgB₂独特的物理性质。当MgB₂中掺入5%的Al后，B-B键的布居数降低至0.32，B-Mg键的布居数降低至0.18，而新形成的Al-B键布居数为0.38。B-B键布居数的降低表明Al掺杂削弱了B-B键的强度，这是因为Al原子替代B原子后，改变了周围电子云的分布，使得B-B之间的电子共享程度降低。B-Mg键布居数的降低则说明Al掺杂也对B-Mg键产生了影响，可能是由于Al原子的引入改变了晶体中的电荷分布，从而减弱了B-Mg之间的离子键相互作用。Al-B键具有较高的布居数，说明Al-B之间形成了较强的共价键，这是由于Al和B的电负性差异较小，易于形成共价键，且Al-B键的形成在一定程度上补偿了B-B键和B-Mg键强度的减弱。在3%Li掺杂的MgB₂体系中，B-B键布居数为0.33，B-Mg键布居数为0.19，Li-B键布居数为0.25。与未掺杂相比，B-B键和B-Mg键的布居数变化较小，说明Li掺杂对B-B键和B-Mg键的影响相对较弱。Li-B键的形成表明Li原子与B原子之间存在一定的相互作用，Li-B键布居数为0.25，说明Li-B键具有一定的强度，这种相互作用可能会对MgB₂的电子结构和物理性质产生影响。由于Li原子填充在间隙位置，其对B-B键和B-Mg键的影响主要是通过改变晶体中的电荷分布和原子间的相互作用力来实现的。通过对比不同掺杂下B-B、B-Mg等键的变化，可以发现Al掺杂对化学键的影响更为显著，它主要通过改变B-B键和B-Mg键的强度来影响MgB₂的结构和性质；而Li掺杂对化学键的影响相对较小，主要是通过与B原子形成新的Li-B键来改变晶体的电子结构和性质。这些结果为深入理解掺杂对MgB₂性能的影响提供了微观层面的依据。3.2.3能带结构与态密度分析能带结构和态密度是描述材料电子结构的重要物理量，它们能够反映电子在材料中的能量分布和状态密度分布情况，对于理解材料的电学、光学等性质具有重要意义。本研究通过第一性原理计算，绘制了未掺杂MgB₂以及Al、Li掺杂MgB₂体系的能带图和态密度图，如图2所示，（a）为未掺杂MgB₂的能带图，（b）为5%Al掺杂MgB₂的能带图，（c）为3%Li掺杂MgB₂的能带图，（d）为未掺杂MgB₂的态密度图，（e）为5%Al掺杂MgB₂的态密度图，（f）为3%Li掺杂MgB₂的态密度图。图2MgB₂及掺杂体系能带图和态密度图在未掺杂MgB₂的能带图中，可以清晰地看到能带的分布情况。费米能级E_F位于能带结构的中间位置，价带顶和导带底之间存在一定的能隙，表明MgB₂是一种具有一定带隙的半导体材料。价带主要由B原子的2p电子和Mg原子的3s电子贡献，导带则主要由B原子的2p电子激发形成。在费米能级附近，能带较为平坦，说明电子的有效质量较大，电子的移动性较差。从态密度图中可以看出，在费米能级处，态密度不为零，这表明在费米能级附近存在一定数量的电子态，这些电子态对MgB₂的电学性质有着重要的影响。当MgB₂中掺入5%的Al后，能带结构发生了明显的变化。能带展宽现象较为显著，尤其是在费米能级附近，能带的斜率增大，这意味着电子的有效质量减小，电子的移动性增强。Al原子的掺入引入了新的电子态，这些电子态与原有的B原子和Mg原子的电子态相互作用，导致能带结构的改变。在态密度图中，费米能级处的态密度发生了明显的变化，相较于未掺杂时，态密度有所增加，这表明Al掺杂增加了费米能级附近的电子态密度。Al原子的3s和3p电子与B原子的2p电子发生杂化，形成了新的电子态，这些新的电子态在费米能级附近贡献了更多的态密度。在3%Li掺杂的MgB₂体系中，能带结构也发生了变化，但变化程度相对较小。能带在费米能级附近的斜率略有增加，说明Li掺杂对电子的有效质量和移动性有一定的影响，但影响不如Al掺杂显著。从态密度图中可以看出，费米能级处的态密度也有所增加，但增加的幅度小于Al掺杂体系。Li原子的2s电子与B原子的2p电子发生相互作用，虽然也引入了新的电子态，但由于Li原子的电子云分布和掺杂浓度的影响，对费米能级附近态密度的贡献相对较小。费米能级附近的态密度与超导性密切相关。根据BCS理论，超导电子配对发生在费米能级附近，态密度越大，电子配对的可能性就越大，超导转变温度也就越高。在本研究中，Al、Li掺杂后费米能级附近态密度的增加，可能会对MgB₂的超导性能产生影响。Al掺杂体系中费米能级附近态密度的显著增加，预示着Al掺杂可能会提高MgB₂的超导转变温度；而Li掺杂体系中态密度的增加幅度较小，对超导转变温度的影响可能相对较弱。但超导性能还受到其他因素的影响，如电子-声子相互作用等，需要进一步综合分析。3.2.4掺杂与超导电性研究超导电性是MgB₂材料的重要特性之一，而掺杂对其超导电性有着显著的影响。结合BCS理论，根据前面的计算结果，深入分析Al、Li掺杂对超导转变温度的影响及其物理机制。BCS理论认为，超导现象是由于电子之间通过晶格振动（声子）的介导形成库珀对，这些库珀对能够在晶格中无电阻地移动，从而导致超导电性。超导转变温度T_c与电子-声子相互作用强度、费米能级处的态密度等因素密切相关，其关系可以用以下公式表示：T_c=\frac{\theta_D}{1.14}\exp\left(-\frac{1}{N(0)V_{ep}}\right)其中，\theta_D是德拜温度，N(0)是费米能级处的态密度，V_{ep}是电子-声子相互作用势。在未掺杂的MgB₂中，电子-声子相互作用主要发生在B原子的2p电子与晶格振动之间。从前面的计算结果可知，费米能级处的态密度为一定值，此时的超导转变温度为T_{c0}。当MgB₂中掺入5%的Al后，从能带结构和态密度分析可知，费米能级附近的态密度N(0)显著增加。这是因为Al原子的掺入引入了新的电子态，这些电子态与原有的B原子和Mg原子的电子态相互作用，使得费米能级附近的电子态密度增大。根据BCS理论公式，N(0)的增大将导致\frac{1}{N(0)V_{ep}}的值减小，从而使得\exp\left(-\frac{1}{N(0)V_{ep}}\right)的值增大，最终导致超导转变温度T_c升高。Al-B键的形成也可能改变了晶格的振动模式，进而影响电子-声子相互作用强度V_{ep}，但具体的影响还需要进一步的研究和分析。在3%Li掺杂的MgB₂体系中，虽然费米能级附近的态密度N(0)也有所增加，但增加的幅度小于Al掺杂体系。Li原子的2s电子与B原子的2p电子相互作用，引入了一定数量的新电子态，但由于Li原子的电子云分布和掺杂浓度的限制，对费米能级附近态密度的贡献相对较小。根据BCS理论公式，N(0)增加幅度较小，使得\frac{1}{N(0)V_{ep}}的减小幅度相对较小，\exp\left(-\frac{1}{N(0)V_{ep}}\right)的增大也相对较小，因此超导转变温度T_c的升高幅度不如Al掺杂明显。Li原子填充在晶格间隙位置，可能对晶格的对称性和振动模式产生一定的影响，从而间接影响电子-声子相互作用，但这种影响相较于Al掺杂更为复杂，需要进一步深入研究。Al、Li掺杂对MgB₂的超导转变温度有着不同程度的影响，其物理机制主要与掺杂引起的费米能级附近态密度的变化以及对电子-声子相互作用的影响有关。通过对掺杂与超导电性的研究，为进一步优化MgB₂超导材料的性能提供了理论指导。3.3本章小结本章基于第一性原理，通过精心构建理论模型并进行精确计算，深入研究了Al、Li掺杂MgB₂系统，全面分析了掺杂对其几何结构、电子结构以及超导电性的影响。在几何结构方面，对未掺杂MgB₂以及Al、Li掺杂MgB₂体系进行几何优化后，发现掺杂导致晶胞参数和原子间距离发生明显变化。5%Al掺杂使晶格常数a增大，B-B键长缩短，这是由于Al原子半径略大于B原子且Al-B键的形成增强了原子间相互作用；3%Li掺杂使晶格常数a减小、c增大，B原子层出现扭曲，Li-B键长为2.05Å，表明Li原子填充在间隙位置对晶体结构产生了显著影响。从晶体结构稳定性来看，Al掺杂使体系总能量略有降低，增强了结构稳定性；Li掺杂使体系总能量略有升高，降低了结构稳定性。键布居分析显示，未掺杂MgB₂中B-B键布居数为0.35，B-Mg键布居数为0.20。5%Al掺杂后，B-B键布居数降至0.32，B-Mg键布居数降至0.18，Al-B键布居数为0.38，说明Al掺杂削弱了B-B键和B-Mg键强度，同时形成了较强的Al-B键；3%Li掺杂后，B-B键布居数为0.33，B-Mg键布居数为0.19，Li-B键布居数为0.25，表明Li掺杂对B-B键和B-Mg键影响较小，主要通过形成Li-B键改变晶体性质。能带结构与态密度分析表明，未掺杂MgB₂具有一定带隙，价带主要由B原子的2p电子和Mg原子的3s电子贡献，导带主要由B原子的2p电子激发形成。5%Al掺杂后，能带展宽，费米能级附近态密度显著增加，这是由于Al原子的3s和3p电子与B原子的2p电子杂化，引入了新的电子态；3%Li掺杂后，能带在费米能级附近斜率略有增加，态密度也有所增加，但幅度小于Al掺杂体系，Li原子的2s电子与B原子的2p电子相互作用引入的新电子态对费米能级附近态密度贡献相对较小。结合BCS理论，探讨了掺杂对超导电性的影响。Al、Li掺杂均使费米能级附近态密度增加，可能影响超导性能。5%Al掺杂体系态密度显著增加，可能提高超导转变温度；3%Li掺杂体系态密度增加幅度较小，对超导转变温度影响相对较弱。超导性能还受电子-声子相互作用等因素影响，需要进一步综合研究。Al、Li掺杂对MgB₂的电子结构和超导性能产生了不同程度的影响，这些研究结果为进一步优化MgB₂超导材料的性能提供了重要的理论依据，也为探索新型超导材料的掺杂策略提供了有益的参考。四、Be、C共掺MgB₂系统研究4.1模型构建与计算方法为深入研究Be、C共掺对MgB₂系统电子结构和光学性质的影响，构建合理的理论模型并选择合适的计算方法至关重要。本研究基于MgB₂的六方晶系结构（空间群为P6/mmm），采用超晶胞模型来模拟Be、C共掺体系。在构建Be、C共掺MgB₂的超晶胞时，选用2\times2\times2的超晶胞，该超晶胞包含32个原子，其中16个Mg原子和16个B原子。考虑到Be原子半径与Mg原子相近，化学性质也有一定的相似性，将Be原子替代Mg原子的位置进行掺杂；而C原子半径与B原子相近，C原子替代B原子的位置进行掺杂。通过精确计算和分析，确定了不同掺杂浓度下Be、C原子在超晶胞中的位置。在研究3%Be和5%C共掺时，在超晶胞中选取1个Mg原子用Be原子替换，选取1个B原子用C原子替换，从而实现对该共掺浓度体系的模拟研究。在确定原子位置后，对晶胞参数进行优化计算。通过第一性原理计算中的几何优化步骤，不断调整晶胞的棱长a和c以及棱间夹角，使得体系的总能量达到最低，原子受力收敛到极小值，从而得到稳定的晶胞参数。在未掺杂的MgB₂超晶胞中，经过优化计算得到的晶格常数a=3.090\mathring{A}，c=3.528\mathring{A}，与实验值和其他理论计算结果相符。对于Be、C共掺体系，随着共掺原子的引入，晶胞参数会发生相应的变化。在3%Be和5%C共掺体系中，优化后的晶格常数a增大至3.098\mathring{A}，c增大为3.535\mathring{A}。这些晶胞参数的变化反映了共掺原子与MgB₂晶格之间的相互作用，对后续的电子结构和光学性质研究具有重要影响。本研究采用基于密度泛函理论（DFT）的CASTEP模块进行计算。在计算过程中，交换关联泛函选择广义梯度近似（GGA）中的Perdew-Burke-Ernzerhof（PBE）泛函，该泛函在处理固体材料的电子结构和性质时能够提供较为准确的结果。赝势选用超软赝势（USPP），它可以有效降低计算量，提高计算效率，同时保证计算精度。平面波截断能设置为400eV，此截断能值是通过对不同截断能下的计算结果进行收敛性测试确定的。当截断能达到400eV时，体系的总能量、原子受力等物理量的变化趋于稳定，满足计算精度要求。布里渊区积分采用Monkhorst-Pack方法，k点网格设置为6\times6\times4。k点网格的选取对计算结果的精度和计算量有显著影响，通过对不同k点网格下的计算结果进行对比分析，发现6\times6\times4的k点网格能够在保证计算精度的同时，兼顾计算效率。在几何优化过程中，使用BFGS算法进行结构优化，设置原子受力收敛标准为0.01eV/Å，能量收敛标准为1\times10^{-6}eV/atom，最大位移收敛标准为5\times10^{-4}Å。当体系的原子受力、能量和最大位移满足这些收敛标准时，认为几何优化达到收敛，得到的结构为体系的稳定结构。在电子结构计算中，采用自洽场（SCF）迭代方法，通过不断迭代计算，使体系的电子密度和能量达到自洽，从而得到准确的电子结构信息。4.2计算结果与讨论4.2.1几何结构与键布居经过几何优化后，得到了未掺杂MgB₂以及Be、C共掺MgB₂体系的稳定晶胞结构，其晶胞结构如图3所示，（a）为未掺杂MgB₂，（b）为3%Be和5%C共掺MgB₂。在未掺杂的MgB₂中，硼原子呈蜂巢状的平面排列，每一个硼原子都与周围三个硼原子相连，形成二维的六角网状结构；镁原子位于硼原子层的上下两侧，呈六方密堆积排列，填充在硼原子层之间的空隙中，这种结构与MgB₂的理论晶体结构一致。图3MgB₂及共掺体系晶胞结构对于3%Be和5%C共掺的MgB₂体系，Be原子替代了部分Mg原子的位置，C原子替代了部分B原子的位置。从优化后的结构可以看出，Be、C共掺导致周围原子的位置发生了明显的调整，以适应新的原子间相互作用。Be原子的掺入使得周围Mg原子的位置发生了一定的位移，C原子的掺入则使得周围B原子的位置发生了变化。Be-B键长为1.65Å，C-B键长为1.58Å。Be-B键长相较于未掺杂时Mg-B键长略有缩短，这是由于Be原子的电负性与Mg原子不同，Be-B之间的化学键具有更强的极性，导致原子间的吸引力增强，从而使键长缩短。C-B键长也相较于未掺杂时B-B键长有所缩短，这是因为C原子的电负性大于B原子，C-B之间形成了较强的共价键，使得原子间距离减小。Be、C共掺也使得晶胞参数发生了变化，晶格常数a增大至3.098\mathring{A}，c增大为3.535\mathring{A}。a和c的增大可能是由于Be、C原子的原子半径与被替代的Mg、B原子半径存在差异，共掺后原子占据了更大的空间，从而导致晶胞在a和c方向上的尺寸增加。通过分析共掺前后晶胞参数的变化以及原子间距离的改变，可以探讨共掺对MgB₂晶体结构稳定性的影响。一般来说，晶体结构的稳定性与体系的能量密切相关，能量越低，结构越稳定。在Be、C共掺体系中，体系的总能量相较于未掺杂时略有升高，这表明Be、C共掺可能会降低晶体结构的稳定性。这可能是因为Be、C原子的掺入破坏了原有的原子堆积方式，导致原子间的相互作用发生了变化，从而使体系的能量升高。键布居分析是研究原子间化学键性质的重要方法，它通过计算原子间的电子分布情况，来揭示化学键的强度和性质。在本研究中，对Be、C共掺后MgB₂系统中各原子间的键布居进行了计算，计算结果如表2所示，其中列出了未掺杂MgB₂以及3%Be和5%C共掺MgB₂体系中B-B、B-Mg等键的布居数。表2MgB₂及共掺体系中键布居数体系B-B键布居数B-Mg键布居数Be-B键布居数（共掺）C-B键布居数（共掺）未掺杂MgB₂0.350.20--3%Be和5%C共掺MgB₂0.300.170.320.36在未掺杂的MgB₂中，B-B键的布居数为0.35，B-Mg键的布居数为0.20。B-B键具有较高的布居数，表明B-B之间存在较强的共价键相互作用，这与硼原子的成键特性相符，硼原子通过共享电子形成稳定的共价键网络，对MgB₂的结构和性质起着重要的支撑作用。B-Mg键的布居数相对较低，说明B-Mg之间的化学键具有一定的离子性，镁原子向硼原子提供电子，形成离子键相互作用，这种离子键与共价键的共存，赋予了MgB₂独特的物理性质。当MgB₂中掺入3%Be和5%C后，B-B键的布居数降低至0.30，B-Mg键的布居数降低至0.17，而新形成的Be-B键布居数为0.32，C-B键布居数为0.36。B-B键布居数的降低表明Be、C共掺削弱了B-B键的强度，这是因为C原子替代B原子后，改变了周围电子云的分布，使得B-B之间的电子共享程度降低。B-Mg键布居数的降低则说明Be、C共掺也对B-Mg键产生了影响，可能是由于Be原子替代Mg原子以及C原子替代B原子后，改变了晶体中的电荷分布，从而减弱了B-Mg之间的离子键相互作用。Be-B键和C-B键具有较高的布居数，说明Be-B和C-B之间形成了较强的共价键，这是由于Be、C与B的电负性差异较小，易于形成共价键，且Be-B键和C-B键的形成在一定程度上补偿了B-B键和B-Mg键强度的减弱。通过对比不同掺杂下B-B、B-Mg等键的变化，可以发现Be、C共掺对化学键的影响较为显著，它主要通过改变B-B键和B-Mg键的强度，以及形成新的Be-B键和C-B键来影响MgB₂的结构和性质。这些结果为深入理解共掺对MgB₂性能的影响提供了微观层面的依据。4.2.2能带与态密度分析能带结构和态密度是描述材料电子结构的重要物理量，它们能够反映电子在材料中的能量分布和状态密度分布情况，对于理解材料的电学、光学等性质具有重要意义。本研究通过第一性原理计算，绘制了未掺杂MgB₂以及Be、C共掺MgB₂体系的能带图和态密度图，如图4所示，（a）为未掺杂MgB₂的能带图，（b）为3%Be和5%C共掺MgB₂的能带图，（c）为未掺杂MgB₂的态密度图，（d）为3%Be和5%C共掺MgB₂的态密度图。图4MgB₂及共掺体系能带图和态密度图在未掺杂MgB₂的能带图中，可以清晰地看到能带的分布情况。费米能级E_F位于能带结构的中间位置，价带顶和导带底之间存在一定的能隙，表明MgB₂是一种具有一定带隙的半导体材料。价带主要由B原子的2p电子和Mg原子的3s电子贡献，导带则主要由B原子的2p电子激发形成。在费米能级附近，能带较为平坦，说明电子的有效质量较大，电子的移动性较差。从态密度图中可以看出，在费米能级处，态密度不为零，这表明在费米能级附近存在一定数量的电子态，这些电子态对MgB₂的电学性质有着重要的影响。当MgB₂中掺入3%Be和5%C后，能带结构发生了明显的变化。能带展宽现象较为显著，尤其是在费米能级附近，能带的斜率增大，这意味着电子的有效质量减小，电子的移动性增强。Be、C共掺引入了新的电子态，这些电子态与原有的B原子和Mg原子的电子态相互作用，导致能带结构的改变。在态密度图中，费米能级处的态密度发生了明显的变化，相较于未掺杂时，态密度有所增加，这表明Be、C共掺增加了费米能级附近的电子态密度。Be原子的2s和2p电子与B原子的2p电子发生杂化，C原子的2s和2p电子也与B原子的2p电子发生杂化，形成了新的电子态，这些新的电子态在费米能级附近贡献了更多的态密度。费米能级附近的态密度与超导性密切相关。根据BCS理论，超导电子配对发生在费米能级附近，态密度越大，电子配对的可能性就越大，超导转变温度也就越高。在本研究中，Be、C共掺后费米能级附近态密度的增加，可能会对MgB₂的超导性能产生影响。Be、C共掺体系中费米能级附近态密度的显著增加，预示着Be、C共掺可能会提高MgB₂的超导转变温度。但超导性能还受到其他因素的影响，如电子-声子相互作用等，需要进一步综合分析。4.2.3超导电性分析超导电性是MgB₂材料的重要特性之一，而掺杂对其超导电性有着显著的影响。结合BCS理论，根据前面的计算结果，深入分析Be、C共掺对超导转变温度的影响及其物理机制。BCS理论认为，超导现象是由于电子之间通过晶格振动（声子）的介导形成库珀对，这些库珀对能够在晶格中无电阻地移动，从而导致超导电性。超导转变温度T_c与电子-声子相互作用强度、费米能级处的态密度等因素密切相关，其关系可以用以下公式表示：T_c=\frac{\theta_D}{1.14}\exp\left(-\frac{1}{N(0)V_{ep}}\right)其中，\theta_D是德拜温度，N(0)是费米能级处的态密度，V_{ep}是电子-声子相互作用势。在未掺杂的MgB₂中，电子-声子相互作用主要发生在B原子的2p电子与晶格振动之间。从前面的计算结果可知，费米能级处的态密度为一定值，此时的超导转变温度为T_{c0}。当MgB₂中掺入3%Be和5%C后，从能带结构和态密度分析可知，费米能级附近的态密度N(0)显著增加。这是因为Be、C共掺引入了新的电子态，这些电子态与原有的B原子和Mg原子的电子态相互作用，使得费米能级附近的电子态密度增大。根据BCS理论公式，N(0)的增大将导致\frac{1}{N(0)V_{ep}}的值减小，从而使得\exp\left(-\frac{1}{N(0)V_{ep}}\right)的值增大，最终导致超导转变温度T_c升高。Be-B键和C-B键的形成也可能改变了晶格的振动模式，进而影响电子-声子相互作用强度V_{ep}，但具体的影响还需要进一步的研究和分析。Be、C共掺对MgB₂的超导转变温度有着显著的影响，其物理机制主要与共掺引起的费米能级附近态密度的变化以及对电子-声子相互作用的影响有关。通过对共掺与超导电性的研究，为进一步优化MgB₂超导材料的性能提供了理论指导。4.2.4光学性质研究光学性质是材料的重要物理性质之一，它反映了材料与光相互作用的特性。通过第一性原理计算，可以得到材料的光电导谱、介电函数等光学性质，从而深入了解材料的光学行为。本研究计算了未掺杂MgB₂以及Be、C共掺MgB₂体系的光电导谱和介电函数，结果如图5所示，（a）为未掺杂MgB₂的光电导谱，（b）为3%Be和5%C共掺MgB₂的光电导谱，（c）为未掺杂MgB₂的介电函数实部，（d）为3%Be和5%C共掺MgB₂的介电函数实部。图5MgB₂及共掺体系光电导谱和介电函数在光电导谱中，未掺杂MgB₂在低能区（0-5eV）有一个较弱的光电导峰，这主要是由于价带中的电子跃迁到导带引起的。随着能量的增加，光电导逐渐减小。当MgB₂中掺入3%Be和5%C后，光电导谱发生了明显的变化。在低能区，光电导峰的强度增强，且峰位发生了移动。这是因为Be、C共掺改变了材料的电子结构，使得电子跃迁的概率和能量发生了变化。新形成的Be-B键和C-B键引入了新的电子态，这些电子态参与了光的吸收和发射过程，导致光电导峰的变化。在高能区，共掺体系的光电导也有所增加，这可能是由于共掺导致能带结构的变化，使得更多的电子跃迁通道被打开。介电函数是描述材料在电场作用下电极化性质的物理量，其实部\varepsilon_1(\omega)反映了材料对光的折射和吸收特性。在未掺杂MgB₂中，介电函数实部在低频区（0-2eV）有一个较大的值，随着频率的增加，逐渐减小。这表明在低频区，MgB₂对光的折射能力较强。当掺入Be、C后，介电函数实部在低频区的值略有减小，且在高频区的变化趋势也发生了改变。这是因为共掺改变了材料的电子云分布和化学键性质，从而影响了材料的电极化响应。Be、C共掺引入的新电子态和化学键变化，使得材料在不同频率下的极化能力发生了变化，进而导致介电函数实部的改变。通过对比不同掺杂下的光学性质差异，可以揭示Be、C共掺对MgB₂光学性质的影响机制。Be、C共掺主要通过改变材料的电子结构和化学键性质，影响电子跃迁和电极化响应，从而改变材料的光电导谱和介电函数。这些结果为MgB₂材料在光电器件等领域的应用提供了理论依据，有助于进一步探索其在光学方面的潜在应用价值。4.3本章小结本章运用第一性原理，通过构建精确的理论模型并实施严谨的计算流程，深入剖析了Be、C共掺MgB₂系统，全面阐释了共掺对其几何结构、电子结构、超导电性以及光学性质的影响。在几何结构层面，对未掺杂MgB₂以及Be、C共掺MgB₂体系进行几何优化后，发现共掺致使晶胞参数和原子间距离产生显著变化。3%Be和5%C共掺使晶格常数a增大至3.098\mathring{A}，c增大为3.535\mathring{A}，Be-B键长为1.65Å，C-B键长为1.58Å。这是由于Be、C原子的原子半径与被替代的Mg、B原子半径存在差异，且Be-B、C-B键的形成改变了原子间的相互作用。从晶体结构稳定性考量，共掺使体系总能量略有升高，可能降低了结构稳定性。键布居分析表明，未掺杂MgB₂中B-B键布居数为0.35，B-Mg键布居数为0.20。3%Be和5%C共掺后，B-B键布居数降至0.30，B-Mg键布居数降至0.17，Be-B键布居数为0.32，C-B键布居数为0.36。这表明共掺削弱了B-B键和B-Mg键强度，同时形成了较强的Be-B键和C-B键，新键的形成在一定程度上补偿了原有键强度的减弱，从微观层面揭示了共掺对MgB₂结构和性质的影响。能带结构与态密度分析显示，未掺杂MgB₂具有一定带隙，价带主要由B原子的2p电子和Mg原子的3s电子贡献，导带主要由B原子的2p电子激发形成。3%Be和5%C共掺后，能带展宽，费米能级附近态密度显著增加，这是因为Be、C共掺引入了新的电子态，Be原子的2s和2p电子、C原子的2s和2p电子与B原子的2p电子杂化，在费米能级附近贡献了更多态密度。结合BCS理论探讨超导电性可知，Be、C共掺使费米能级附近态密度增加，可能提高超导转变温度。根据BCS理论公式，态密度N(0)增大导致\frac{1}{N(0)V_{ep}}减小，\exp\left(-\frac{1}{N(0)V_{ep}}\right)增大，从而使超导转变温度T_c升高，但超导性能还受电子-声子相互作用等因素影响，需进一步综合研究。在光学性质方面，计算了未掺杂MgB₂以及Be、C共掺MgB₂体系的光电导谱和介电函数。共掺后光电导谱在低能区峰强度增强且峰位移动，高能区光电导增加，这是因为共掺改变了电子结构，引入的新电子态参与光的吸收和发射过程。介电函数实部在低频区的值略有减小，高频区变化趋势改变，这是由于共掺改变了电子云分布和化学键性质，影响了电极化响应。Be、C共掺对MgB₂的电子结构、超导电性和光学性质产生了显著影响，共掺补偿作用在各性质中均有体现，这些研究成果为进一步优化MgB₂超导材料的性能以及探索其在光电器件等领域的应用提供了重要的理论依据。五、Ag掺杂MgB₂系统研究5.1研究模型与计算过程为深入探究Ag掺杂对MgB₂系统电子结构和光学性质的影响，构建准确合理的研究模型并采用科学严谨的计算过程至关重要。本研究基于MgB₂的六方晶系结构（空间群为P6/mmm），利用MaterialsStudio软件中的CASTEP模块开展第一性原理计算。在构建Ag掺杂MgB₂的超晶胞模型时，选用2\times2\times2的超晶胞，该超晶胞包含32个原子，其中16个Mg原子和16个B原子。考虑到Ag原子半径与Mg原子半径有一定差异，且化学性质也有所不同，本研究将Ag原子替代Mg原子的位置进行掺杂。通过精确计算和分析，确定了不同掺杂浓度下Ag原子在超晶胞中的位置。在研究5%Ag掺杂时，在超晶胞中选取1个Mg原子用Ag原子替换，从而实现对该掺杂浓度体系的模拟研究。在确定原子位置后，对晶胞参数进行细致的优化计算。通过第一性原理计算中的几何优化步骤，不断调整晶胞的棱长a和c以及棱间夹角，使得体系的总能量达到最低，原子受力收敛到极小值，从而得到稳定的晶胞参数。在未掺杂的MgB₂超晶胞中，经过优化计算得到的晶格常数a=3.090\mathring{A}，c=3.528\mathring{A}，与实验值和其他理论计算结果相符。对于5%Ag掺杂体系，随着Ag原子的引入，晶胞参数会发生相应的变化。优化后的晶格常数a减小至3.082\mathring{A}，c减小为3.522\mathring{A}。这些晶胞参数的变化反映了Ag原子与MgB₂晶格之间的相互作用，对后续的电子结构和光学性质研究具有重要影响。在计算过程中，交换关联泛函选择广义梯度近似（GGA）中的Perdew-Burke-Ernzerhof（PBE）泛函，该泛函在处理固体材料的电子结构和性质时能够提供较为准确的结果。赝势选用超软赝势（USPP），它可以有效降低计算量，提高计算效率，同时保证计算精度。平面波截断能设置为400eV，此截断能值是通过对不同截断能下的计算结果进行收敛性测试确定的。当截断能达到400eV时，体系的总能量、原子受力等物理量的变化趋于稳定，满足计算精度要求。布里渊区积分采用Monkhorst-Pack方法，k点网格设置为6\times6\times4。k点网格的选取对计算结果的精度和计算量有显著影响，通过对不同k点网格下的计算结果进行对比分析，发现6\times6\times4的k点网格能够在保证计算精度的同时，兼顾计算效率。在几何优化过程中，使用BFGS算法进行结构优化，设置原子受力收敛标准为0.01eV/Å，能量收敛标准为1\times10^{-6}eV/atom，最大位移收敛标准为5\times10^{-4}Å。当体系的原子受力、能量和最大位移满足这些收敛标准时，认为几何优化达到收敛，得到的结构为体系的稳定结构。在电子结构计算中，采用自洽场（SCF）迭代方法，通过不断迭代计算，使体系的电子密度和能量达到自洽，从而得到准确的电子结构信息。5.2结果分析5.2.1几何优化与键布居经过几何优化后，获得了未掺杂MgB₂以及5%Ag掺杂MgB₂体系的稳定晶胞结构，如图6所示，（a）为未掺杂MgB₂，（b）为5%Ag掺杂MgB₂。在未掺杂的MgB₂中，硼原子呈蜂巢状的平面排列，每一个硼原子都与周围三个硼原子相连，形成二维的六角网状结构；镁原子位于硼原子层的上下两侧，呈六方密堆积排列，填充在硼原子层之间的空隙中，这种结构与MgB₂的理论晶体结构一致。图6MgB₂及Ag掺杂体系晶胞结构对于5%Ag掺杂的MgB₂体系，Ag原子替代了部分Mg原子的位置。从优化后的结构可以看出，Ag原子的掺入导致周围原子的位置发生了明显的调整，以适应新的原子间相互作用。Ag原子的掺入使得周围B原子的位置发生了一定的位移，B原子与Ag原子之间的距离有所改变。Ag-B键长为2.15Å，相较于未掺杂时Mg-B键长2.25Å略有缩短，这是由于Ag原子的电负性与Mg原子不同，Ag-B之间的化学键具有更强的极性，导致原子间的吸引力增强，从而使键长缩短。Ag掺杂也使得晶胞参数发生了变化，晶格常数a减小至3.082\mathring{A}，c减小为3.522\mathring{A}。a和c的减小可能是由于Ag原子的原子半径小于Mg原子，替代Mg原子后占据的空间减小，从而导致晶胞在a和c方向上的尺寸收缩。通过分析掺杂前后晶胞参数的变化以及原子间距离的改变，可以探讨Ag掺杂对MgB₂晶体结构稳定性的影响。一般来说，晶体结构的稳定性与体系的能量密切相关，能量越低，结构越稳定。在Ag掺杂体系中，体系的总能量相较于未掺杂时略有升高，这表明Ag掺杂可能会降低晶体结构的稳定性。这可能是因为Ag原子的掺入破坏了原有