基于第一性原理与模型构建的二维和三维拓扑半金属特性研究

基于第一性原理与模型构建的二维和三维拓扑半金属特性研究一、引言1.1研究背景与意义拓扑半金属作为凝聚态物理领域的重要研究对象，近年来吸引了广泛的关注。其独特的电子结构和拓扑性质，不仅为基础物理研究带来了革命性的突破，也在材料科学和电子学等领域展现出巨大的应用潜力，成为推动相关领域发展的关键因素之一。在凝聚态物理中，拓扑半金属是一类具有特殊电子结构的材料，其能带结构呈现出独特的线性色散关系，电子在其中表现出类似于无质量粒子的行为。这种特殊的电子结构赋予了拓扑半金属许多新奇的物理性质，如超高的载流子迁移率、独特的光学响应以及奇异的磁输运性质等。这些性质使得拓扑半金属在高速电子学、量子计算、传感器等领域具有潜在的应用价值，可能为下一代信息技术的发展提供关键材料支持。从材料科学的角度来看，拓扑半金属的研究为新型材料的设计和开发提供了新的思路和方法。通过对拓扑半金属的深入研究，科学家们可以探索如何利用其独特的物理性质，开发出具有高性能的电子器件、磁性材料和光学材料等。例如，拓扑半金属的高载流子迁移率使其有望应用于高速电子器件中，提高器件的运行速度和降低能耗；其独特的光学响应则可能为新型光电器件的设计提供基础，实现更高效率的光-电转换。在电子学领域，拓扑半金属的出现为解决当前电子器件面临的一些挑战提供了新的途径。随着集成电路技术的不断发展，器件尺寸越来越小，面临着能耗增加、散热困难等问题。拓扑半金属的低能耗和高载流子迁移率等特性，使其有可能成为下一代电子器件的理想材料，为实现高性能、低能耗的电子器件提供了可能。此外，拓扑半金属在自旋电子学中的潜在应用也备受关注，有望推动自旋电子学器件的发展，实现信息的高效存储和处理。1.2国内外研究现状国内外关于二维和三维拓扑半金属的研究在理论和实验方面都取得了显著进展。在理论研究上，科学家们通过多种理论模型和计算方法，深入探究拓扑半金属的电子结构和拓扑性质。密度泛函理论（DFT）被广泛用于计算材料的电子结构和能量，从而预测拓扑半金属的存在及其性质。通过DFT计算，研究人员成功预言了多种拓扑半金属材料，如Na₃Bi和Cd₃As₂等狄拉克半金属。此外，紧束缚模型、k・p微扰理论等也被用于研究拓扑半金属的低能激发特性和拓扑不变量，为理解其物理性质提供了理论基础。在实验研究方面，随着先进实验技术的不断发展，如角分辨光电子能谱（ARPES）、扫描隧道显微镜（STM）和输运测量等，为拓扑半金属的研究提供了有力手段。ARPES能够直接测量材料的电子能带结构，从而验证理论预言的拓扑半金属态。通过ARPES实验，在Na₃Bi和Cd₃As₂等材料中观测到了狄拉克锥型的能带结构，证实了它们作为狄拉克半金属的特性。STM则可以在实空间对材料表面的原子和电子态进行成像，研究拓扑半金属的表面态和缺陷等性质。输运测量能够研究拓扑半金属的电学和磁学输运性质，如发现外尔半金属中的手性负磁阻效应等。尽管国内外在拓扑半金属研究上取得了诸多成果，但当前研究仍存在一些问题和挑战。在材料制备方面，高质量、大尺寸的拓扑半金属材料的制备仍然面临困难，这限制了对其物理性质的深入研究和实际应用。材料中的杂质和缺陷会对拓扑半金属的电子结构和物理性质产生显著影响，如何精确控制材料的质量和纯度是亟待解决的问题。在理论研究方面，虽然现有的理论模型能够解释部分拓扑半金属的性质，但对于一些复杂的拓扑半金属体系，理论描述还不够完善，需要进一步发展更加精确的理论模型。此外，拓扑半金属与其他材料的复合和集成技术还不够成熟，如何实现拓扑半金属与其他材料的有效复合，以发挥其独特性质，也是未来研究的重要方向。1.3研究内容与方法本研究旨在通过第一性原理和模型研究，深入探讨几种二维和三维拓扑半金属的电子结构、拓扑性质及相关特性。具体研究内容包括：运用第一性原理计算方法，基于密度泛函理论，对选定的二维和三维拓扑半金属材料进行原子结构优化和电子结构计算，精确获取其能带结构、电子态密度等信息，从而分析其电子结构特征。通过计算拓扑不变量，如陈数、Z₂不变量等，确定拓扑半金属的拓扑性质，明确其拓扑分类，并研究拓扑性质与材料结构和电子相互作用的关系。构建合适的理论模型，如紧束缚模型，进一步研究拓扑半金属的低能激发特性和电子输运性质，解释实验中观察到的现象，预测新的物理性质。结合第一性原理计算和模型研究结果，探讨拓扑半金属在电子学、能源等领域的潜在应用，为其实际应用提供理论指导。在研究方法上，以第一性原理计算为核心，利用VASP、QuantumESPRESSO等计算软件进行材料的电子结构计算。通过调整计算参数和方法，确保计算结果的准确性和可靠性。同时，运用理论模型进行辅助分析，将理论模型的结果与第一性原理计算结果相互验证和补充，深入理解拓扑半金属的物理性质。此外，密切关注国内外相关实验研究进展，将理论研究结果与实验数据进行对比和分析，进一步完善理论模型，提高对拓扑半金属的认识。二、理论基础与研究方法2.1第一性原理计算基础2.1.1密度泛函理论（DFT）密度泛函理论（DFT）是第一性原理计算的核心理论框架，在多电子体系的研究中发挥着关键作用。其基本原理基于Hohenberg-Kohn定理，该定理指出，多电子体系的基态能量是电子密度的唯一泛函。这意味着，通过确定电子密度，就能获取体系的所有基态性质，从而将复杂的多体问题转化为相对简单的单电子问题。Hohenberg-Kohn第一定理表明，对于一个处在外部势场V_{ext}(r)中的多电子体系，其基态电子密度\rho(r)唯一地决定了外部势场（除了一个常数）。由于外部势场决定了体系的哈密顿量，也就意味着电子密度唯一地决定了体系的所有性质。Hohenberg-Kohn第二定理则进一步阐述，存在一个以基态电子密度为变量的普适能量泛函E[\rho]，对于给定的外部势场V_{ext}(r)，当电子密度取体系的基态电子密度时，能量泛函E[\rho]达到最小值，且这个最小值就是体系的基态能量。在实际应用中，为了求解多电子体系的能量和波函数，Kohn和Sham提出了著名的Kohn-Sham方程。该方程的推导过程是将多电子体系的总能量E表示为电子动能T、电子-电子相互作用能（包括库仑相互作用能V_{ee}和交换关联能E_{xc}）以及电子与外部势场的相互作用能V_{ext}之和，即E=T+V_{ee}+E_{xc}+V_{ext}。为了将多体问题转化为单体问题，Kohn-Sham引入了一个虚构的无相互作用电子体系，这个体系具有与真实体系相同的电子密度\rho(r)。对于这个虚构的体系，电子的动能可以用单电子波函数\psi_{i}(r)表示为T_{s}=-\frac{\hbar^{2}}{2m}\sum_{i}\int\psi_{i}^{*}(r)\nabla^{2}\psi_{i}(r)dr，库仑相互作用能可以表示为经典的静电能V_{H}=\frac{1}{2}\iint\frac{\rho(r)\rho(r')}{|r-r'|}drdr'。交换关联能E_{xc}则包含了电子之间的交换作用和关联作用，目前还无法精确求解，通常采用近似方法处理，如局域密度近似（LDA）和广义梯度近似（GGA）等。通过变分原理，对能量泛函E[\rho]关于单电子波函数\psi_{i}(r)求变分，得到Kohn-Sham方程：\left[-\frac{\hbar^{2}}{2m}\nabla^{2}+V_{ext}(r)+V_{H}(r)+V_{xc}(r)\right]\psi_{i}(r)=\epsilon_{i}\psi_{i}(r)其中，V_{xc}(r)=\frac{\deltaE_{xc}[\rho]}{\delta\rho(r)}是交换关联势，\epsilon_{i}是单电子本征能量。通过求解Kohn-Sham方程，可以得到单电子波函数\psi_{i}(r)和电子密度\rho(r)=\sum_{i}|\psi_{i}(r)|^{2}，进而计算出体系的总能量和其他物理性质。DFT在处理多电子体系时具有显著优势。它能够有效地简化多体问题的计算复杂度，使得对大规模体系的计算成为可能。在研究材料的电子结构时，DFT可以精确计算出能带结构、电子态密度等关键信息，为理解材料的物理性质提供重要依据。然而，DFT也存在一定的局限性，主要体现在对交换关联能的近似处理上。目前的近似方法虽然在许多情况下能够给出合理的结果，但对于一些强关联体系，如过渡金属氧化物等，这些近似方法往往无法准确描述电子之间的强相互作用，导致计算结果与实验值存在较大偏差。2.1.2赝势方法在第一性原理计算中，赝势方法是一种重要的技术手段，用于简化对电子-离子相互作用的描述。其核心概念是用一个相对简单的赝势来代替真实的离子势，从而在不损失关键物理信息的前提下，降低计算的复杂性。在真实的原子体系中，原子核与内层电子构成离子实，价电子在离子实的势场中运动。价电子波函数在离子实内部的区域变化剧烈，存在很多节点，这是由于需要使价电子与内层电子波函数正交，径向函数乘积积分为零，因而离子实内部出现节点使一部分区域为正，一部分为负。这种剧烈变化使得直接求解包含所有电子的薛定谔方程计算量极大。赝势的作用就是把离子实的内部势能用假想的势能取代真实的势能，但在求解波动方程时，不改变能量本征值和离子实之间区域的波函数。由赝势求出的波函数叫赝波函数，在离子实之间的区域真实的势和赝势给出同样的波函数。通过这种方式，去除了波函数在内壳层的剧烈振荡部分，使其在内壳层区域更加平缓，从而允许只以相对很少的平面波来展开波函数，减少了计算量。常见的赝势类型包括模守恒赝势（Norm-ConservingPseudopotential）和超软赝势（Ultra-SoftPseudopotential）。模守恒赝势要求赝波函数和真实波函数在离子实外的区域具有相同的对数导数，以保证在该区域的散射性质相同，从而守恒波函数的一些重要性质，如电荷密度等。超软赝势则进一步放宽了对波函数的限制，它允许赝波函数在离子实附近有更大的灵活性，通过引入投影算子来保证计算的准确性。超软赝势的优势在于可以使用较少的平面波基组来描述电子波函数，大大提高了计算效率，尤其适用于处理包含大原子或复杂体系的计算。赝势方法在第一性原理计算中有着广泛的应用。在计算材料的晶体结构时，通过使用赝势可以快速准确地得到原子的平衡位置和晶格常数；在研究材料的电子结构和光学性质时，赝势方法能够有效地计算能带结构、态密度等物理量，为理解材料的电学和光学行为提供理论支持。2.1.3平面波基组平面波基组是第一性原理计算中用于展开电子波函数的一种常用基组，具有独特的数学形式和物理意义。在量子力学中，根据波函数的完备性，任何一个函数都可以用一组完备的基函数展开。平面波函数由于其简单的形式和良好的数学性质，成为了展开电子波函数的理想选择。在周期性体系中，如晶体，电子波函数满足布洛赫定理，即\psi_{k}(r)=e^{ik\cdotr}u_{k}(r)，其中k是波矢，u_{k}(r)是具有晶格周期性的函数，u_{k}(r+R)=u_{k}(r)，R是晶格矢量。根据傅里叶变换的性质，具有周期性的函数u_{k}(r)可以展开为平面波的叠加：u_{k}(r)=\sum_{G}c_{k+G}e^{i(k+G)\cdotr}其中，G是倒格矢，c_{k+G}是展开系数。因此，电子波函数\psi_{k}(r)可以表示为平面波的线性组合：\psi_{k}(r)=\sum_{G}c_{k+G}e^{i(k+G)\cdotr}平面波截断能是平面波基组中的一个关键参数，它决定了参与展开的平面波的数量。截断能E_{cut}定义为最大的平面波动能，即E_{cut}=\frac{\hbar^{2}(k+G)_{max}^{2}}{2m}，其中(k+G)_{max}是最大的波矢。当截断能增大时，参与展开的平面波数量增多，波函数的描述更加精确，计算结果也更加准确。但同时，计算量会随着平面波数量的增加而急剧增大，导致计算时间增长和计算资源消耗增加。因此，选择合适的截断能至关重要，需要在计算精度和计算效率之间进行权衡。通常的做法是通过对不同截断能下的计算结果进行收敛性测试，观察体系能量、原子力等物理量随截断能的变化情况，当这些物理量在某一截断能下基本保持不变时，即认为计算结果已收敛，该截断能即为合适的取值。在计算晶体的能带结构时，若截断能设置过低，可能会导致能带结构的不准确，无法正确反映材料的电子性质；而过高的截断能虽然能提高计算精度，但会大大增加计算成本。在实际应用中，对于简单的晶体结构，较低的截断能可能就足以得到较为准确的结果；而对于复杂的体系，如含有过渡金属元素的化合物，可能需要较高的截断能才能保证计算精度。2.2拓扑半金属相关理论2.2.1拓扑不变量与Berry曲率拓扑不变量是拓扑学中的核心概念，在拓扑半金属的研究中，它用于表征材料的拓扑性质，这些性质在连续形变下保持不变，与材料的具体几何形状和尺寸无关。常见的拓扑不变量包括陈数（Chernnumber）和Z_2指标等，它们各自具有独特的物理意义和计算方法，能够为拓扑半金属的分类和性质研究提供重要依据。陈数是一种用于描述二维拓扑材料的拓扑不变量，它与材料的能带结构和Berry曲率密切相关。在二维体系中，Berry曲率\Omega_{n}(k)可以被看作是动量空间中的一种“曲率”，它描述了电子在动量空间中运动时的几何相位变化。对于第n个能带，Berry曲率的定义为：\Omega_{n}(k)=i\sum_{m\neqn}\frac{\langleu_{n}(k)|\nabla_{k}H(k)|u_{m}(k)\rangle\times\langleu_{m}(k)|\nabla_{k}H(k)|u_{n}(k)\rangle}{(E_{n}(k)-E_{m}(k))^{2}}其中，|u_{n}(k)\rangle是第n个能带的布洛赫波函数，H(k)是哈密顿量，E_{n}(k)是第n个能带的能量。陈数C_{n}则是Berry曲率在整个布里渊区上的积分：C_{n}=\frac{1}{2\pi}\iint_{BZ}\Omega_{n}(k)dk_{x}dk_{y}陈数是一个整数，它的非零值表明材料具有非平凡的拓扑性质，这种材料被称为陈绝缘体或量子反常霍尔绝缘体。在量子反常霍尔效应中，材料在零磁场下就能够表现出量子化的霍尔电导，其霍尔电导值与陈数成正比，即\sigma_{xy}=C_{n}\frac{e^{2}}{h}，其中e是电子电荷，h是普朗克常数。Z_2指标是用于描述三维拓扑材料的拓扑不变量，它通过对材料的能带结构进行奇偶性分析来确定。在三维体系中，Z_2指标的计算较为复杂，通常需要考虑材料在不同动量空间方向上的性质。具体来说，Z_2指标可以通过计算时间反演不变动量点（TRIM）处的波函数相位来确定。如果在这些特殊点处的波函数相位满足一定的奇偶性条件，则材料具有非平凡的Z_2拓扑性质，被称为拓扑绝缘体。拓扑绝缘体的表面存在着受拓扑保护的表面态，这些表面态具有独特的电子结构和输运性质，如表面态电子的色散关系呈现线性，且具有无质量狄拉克费米子的特征，这使得拓扑绝缘体在自旋电子学和量子计算等领域具有潜在的应用价值。Berry曲率与拓扑性质之间存在着紧密的联系。Berry曲率在动量空间中的分布决定了拓扑不变量的值，进而决定了材料的拓扑性质。当Berry曲率在布里渊区上的积分不为零时，材料具有非平凡的拓扑性质，存在着受拓扑保护的边界态或表面态。在拓扑半金属中，Berry曲率的奇异性往往与狄拉克点、外尔点等特殊的能带结构特征相关，这些奇异点周围的Berry曲率分布决定了拓扑半金属的拓扑分类和物理性质。2.2.2能带结构与费米面分析能带结构和费米面是研究拓扑半金属电子结构的重要工具，通过对它们的分析可以揭示拓扑半金属的许多关键特征，如狄拉克点、外尔点等，从而深入理解拓扑半金属的物理性质。能带结构描述了电子在晶体中的能量与动量之间的关系，它是研究固体电子性质的基础。在拓扑半金属中，能带结构呈现出独特的特征。狄拉克半金属具有线性色散的狄拉克锥型能带结构，在狄拉克点处，导带和价带相切，且具有线性的色散关系，电子在狄拉克点附近表现出类似于无质量相对论粒子的行为，其能量与动量满足E(k)=\pm\hbarv_{F}|k-k_{D}|，其中v_{F}是费米速度，k_{D}是狄拉克点的波矢。外尔半金属则是在狄拉克半金属的基础上，由于时间反演对称性或空间反演对称性的破缺，狄拉克点分裂为一对具有相反手性的外尔点。外尔点是动量空间中的一种拓扑缺陷，其周围的能带结构也呈现出线性色散关系，但与狄拉克点不同的是，外尔点具有手性，类似于三维空间中的磁单极子，其手性决定了电子的运动方向和输运性质。费米面是指在绝对零度下，电子占据的最高能量状态所对应的等能面。在拓扑半金属中，费米面的形状和性质与能带结构密切相关。通过分析费米面，可以确定材料中电子的分布和运动状态。在狄拉克半金属中，费米面通常是一个围绕狄拉克点的球面，其半径与电子的浓度有关。当费米面与狄拉克点相交时，会出现零能的狄拉克费米子，这些狄拉克费米子具有高迁移率和独特的输运性质，使得狄拉克半金属在高速电子学和量子器件中具有潜在的应用价值。在外尔半金属中，费米面可能呈现出复杂的形状，由于外尔点的存在，费米面会出现一些特殊的特征，如费米面的嵌套和连接不同外尔点的费米弧等。费米弧是外尔半金属表面态的重要特征，它连接着具有相反手性的外尔点在表面上的投影，是外尔半金属拓扑性质的直接体现，在输运测量中，费米弧的存在会导致一些奇异的物理现象，如手性负磁阻效应等。为了确定拓扑半金属的能带结构和费米面，通常采用第一性原理计算方法和实验测量手段相结合的方式。通过第一性原理计算，可以精确地计算出材料的电子结构和能带结构，预测狄拉克点、外尔点等特征的位置和性质。常用的计算软件如VASP、QuantumESPRESSO等，基于密度泛函理论，能够有效地计算材料的电子结构和能量。实验上，则可以利用角分辨光电子能谱（ARPES）、扫描隧道显微镜（STM）等技术直接测量材料的能带结构和费米面。ARPES能够测量材料表面的电子能量和动量分布，从而得到能带结构的信息，通过ARPES实验，已经在许多拓扑半金属材料中观测到了狄拉克锥和外尔点等特征；STM则可以在实空间对材料表面的原子和电子态进行成像，直接观察费米面的形状和电子的局域态密度，为研究拓扑半金属的电子结构提供了直观的证据。2.3模型构建方法2.3.1紧束缚模型紧束缚模型是凝聚态物理中用于描述电子在晶体中运动的一种重要理论模型，其基本假设基于电子在原子附近的行为。该模型认为，电子主要被原子核吸引，并受到相邻原子核的弱相互作用。在晶体中，原子按一定的晶格结构排列，电子在这样的周期性势场中运动。紧束缚模型假设电子在某个原子附近时，主要受到该原子势场的作用，而受到其他原子势场的作用相对较弱，可以看作是一种微扰。基于上述假设，紧束缚模型通过将原子轨道线性组合（LCAO）的方式来构建晶体中的电子波函数。以简单的一维晶格为例，设晶格常数为a，原子轨道为\phi(r-R_i)，其中R_i=ia为第i个原子的位置，则晶体中的电子波函数\psi(r)可以表示为：\psi(r)=\sum_{i}c_{i}\phi(r-R_{i})其中，c_i是系数，表示电子在第i个原子轨道上出现的概率幅。通过将这个波函数代入薛定谔方程，并利用近邻近似等方法进行求解，可以得到电子的能级和波函数。在近邻近似下，只考虑最近邻原子之间的相互作用，此时哈密顿量矩阵元可以表示为：H_{ij}=\langle\phi(r-R_{i})|H|\phi(r-R_{j})\rangle当i=j时，H_{ii}=\epsilon_0，表示电子在单个原子上的能量；当三、二维拓扑半金属的第一性原理研究3.1二维拓扑半金属材料的选择与结构建模3.1.1材料体系介绍在二维拓扑半金属的研究中，选取了具有代表性的材料体系ZrTe₂和Pt₂HgSe₃进行深入探究。ZrTe₂作为一种典型的二维过渡金属硫族化合物，具有独特的晶体结构和物理性质。其晶体结构属于六方晶系，空间群为P6₃/mmc，原子在二维平面上呈周期性排列，形成了层状结构。每一层由Zr原子和Te原子组成，Zr原子位于层的中心位置，被六个Te原子以八面体配位方式包围，层与层之间通过范德华力相互作用结合在一起。这种层状结构赋予了ZrTe₂一些特殊的性质，如较强的各向异性，电子在平面内和平面间的传输性质存在显著差异。在电学性质方面，ZrTe₂表现出金属性，其载流子迁移率较高，这与它的电子结构密切相关。理论研究表明，ZrTe₂在常压下具有非平庸的拓扑性质，其能带结构中存在狄拉克点，电子在狄拉克点附近表现出线性色散关系，类似于无质量的狄拉克费米子，这种特性使得ZrTe₂在拓扑电子学领域具有潜在的应用价值。Pt₂HgSe₃是一种在巴西发现的天然可解理矿物，近期的研究发现它具有独特的拓扑性质。其晶体结构较为复杂，属于三方晶系，空间群为P3̅m1。在晶体结构中，Pt、Hg和Se原子通过共价键和离子键相互作用，形成了特定的原子排列方式。这种结构使得Pt₂HgSe₃在电子结构上表现出独特的特征，理论计算预测该材料单层时是一种带隙高达0.5eV的量子自旋霍尔绝缘体，而体材料是一种罕见的双拓扑半金属。在双拓扑半金属态下，Pt₂HgSe₃的能带结构中存在多个狄拉克点和外尔点，这些特殊的能带结构特征赋予了它丰富的拓扑性质和物理现象，如表面态的存在、拓扑保护的输运性质等，使其成为研究拓扑物理和新型量子材料的重要体系。3.1.2结构优化与稳定性分析利用第一性原理计算方法，基于密度泛函理论，使用VASP软件对ZrTe₂和Pt₂HgSe₃的原子结构进行优化。在计算过程中，采用广义梯度近似（GGA）来处理交换关联能，选用平面波基组来展开电子波函数，并设置合适的平面波截断能和k点网格，以确保计算结果的准确性和收敛性。通过结构优化，得到了材料的平衡晶格常数、原子坐标以及晶胞体积等参数。对于ZrTe₂，优化后的晶格常数a和c与实验值或其他理论计算值具有较好的一致性，验证了计算方法的可靠性。在优化过程中，通过最小化体系的总能量，使原子逐渐调整到最稳定的位置，从而得到了稳定的晶体结构。为了分析结构的稳定性，计算了材料的形成能和结合能。形成能是指从组成元素的稳定单质状态形成化合物时的能量变化，其计算公式为：E_{form}=E_{total}-\sum_{i}n_{i}E_{i}其中，E_{total}是化合物的总能量，n_{i}是组成元素i的原子数，E_{i}是元素i的稳定单质的能量。形成能的负值越大，表明化合物越稳定。对于ZrTe₂，计算得到的形成能为负值，说明其在热力学上是稳定的，能够稳定存在。结合能则是指将晶体中的原子完全分离所需的能量，它反映了原子之间的相互作用强度。结合能越大，晶体结构越稳定。通过计算ZrTe₂和Pt₂HgSe₃的结合能，发现它们的结合能都较大，表明这两种材料的原子之间具有较强的相互作用，晶体结构较为稳定。此外，还通过计算声子谱来进一步验证结构的稳定性。如果声子谱中不存在虚频，说明晶体结构在动力学上是稳定的。计算结果表明，ZrTe₂和Pt₂HgSe₃的声子谱中均无虚频，进一步证实了它们的结构稳定性。在结构优化过程中，观察到体系能量随着原子位置的调整逐渐降低，最终收敛到一个稳定的值，这也直观地表明了结构优化的有效性和结构的稳定性。3.2电子结构与拓扑性质计算3.2.1能带结构与态密度计算通过第一性原理计算，得到了ZrTe₂和Pt₂HgSe₃的能带结构和态密度。对于ZrTe₂，其能带结构呈现出明显的特征。在费米能级附近，存在着线性色散的能带交叉，形成了狄拉克锥。狄拉克锥的存在表明ZrTe₂具有拓扑半金属的特性，电子在狄拉克点附近表现出类似于无质量相对论粒子的行为。通过对能带结构的分析，确定了狄拉克点的位置和能量，以及狄拉克锥的色散关系。同时，计算了ZrTe₂的态密度，发现态密度在狄拉克点处呈现出尖锐的峰，这与能带结构中狄拉克锥的特征相吻合。态密度的分布反映了电子在不同能量状态下的分布情况，在费米能级附近，态密度主要由Zr原子的d轨道和Te原子的p轨道贡献，这表明这些轨道在电子的输运和拓扑性质中起着重要作用。Pt₂HgSe₃的能带结构更为复杂，由于其双拓扑半金属的特性，存在多个狄拉克点和外尔点。在能带结构中，可以观察到多个线性色散的区域，这些区域对应着不同的狄拉克点和外尔点。通过仔细分析能带结构，确定了各个狄拉克点和外尔点的位置、能量以及手性等特征。计算得到的态密度也显示出丰富的结构，在费米能级附近，态密度存在多个峰和谷，这些特征与能带结构中的狄拉克点和外尔点密切相关。通过对态密度的分波分析，明确了Pt、Hg和Se原子的不同轨道对态密度的贡献，发现Pt原子的d轨道和Se原子的p轨道在费米能级附近的态密度贡献较大，这对理解Pt₂HgSe₃的电子结构和拓扑性质具有重要意义。在计算能带结构和态密度时，采用了高精度的计算参数和方法，以确保结果的准确性。通过改变k点网格的密度和平面波截断能等参数进行收敛性测试，发现当k点网格密度达到一定值，平面波截断能足够大时，计算结果基本收敛，能够准确反映材料的电子结构特征。3.2.2Berry曲率与拓扑不变量计算为了确定ZrTe₂和Pt₂HgSe₃的拓扑性质，计算了它们的Berry曲率和拓扑不变量。对于ZrTe₂，首先计算了其在动量空间中的Berry曲率分布。Berry曲率在狄拉克点附近呈现出明显的峰值，表明狄拉克点周围存在着非平凡的拓扑结构。通过对Berry曲率在整个布里渊区上的积分，得到了拓扑不变量陈数。计算结果表明，ZrTe₂的陈数不为零，这证实了它具有非平庸的拓扑性质，属于拓扑半金属。陈数的非零值意味着ZrTe₂存在受拓扑保护的边界态，这些边界态具有独特的电子结构和输运性质，为其在拓扑电子学中的应用提供了基础。对于Pt₂HgSe₃，由于其复杂的能带结构，计算拓扑不变量相对复杂。除了计算陈数外，还需要考虑Z_2指标等其他拓扑不变量。通过对动量空间中不同方向上的Berry曲率进行分析，结合对称性分析，确定了Pt₂HgSe₃的拓扑不变量。计算结果表明，Pt₂HgSe₃具有非平凡的拓扑性质，其拓扑不变量的值与理论预测的双拓扑半金属的特征相符。在计算过程中，利用了对称性的性质来简化计算，通过对晶体结构的对称性分析，确定了一些特殊的动量空间点和路径，在这些点和路径上计算Berry曲率和拓扑不变量，不仅减少了计算量，还提高了计算结果的准确性。通过计算Berry曲率和拓扑不变量，明确了ZrTe₂和Pt₂HgSe₃的拓扑性质，为进一步研究它们的物理性质和应用提供了重要依据。3.3外界因素对二维拓扑半金属性质的影响3.3.1压力效应研究压力作为一种有效的外部调控手段，能够显著改变二维拓扑半金属的晶体结构、电子结构和拓扑性质。对于ZrTe₂，通过第一性原理计算研究了压力对其晶体结构的影响。随着压力的增加，ZrTe₂的晶格常数逐渐减小，层间距也相应减小，这是由于原子间的距离在压力作用下被压缩。在电子结构方面，压力导致ZrTe₂的能带结构发生明显变化。狄拉克点的位置和能量发生移动，狄拉克锥的形状也有所改变，这是因为压力改变了原子间的相互作用，进而影响了电子的能量和波函数。在拓扑性质上，压力会诱导ZrTe₂发生拓扑相变。当压力达到一定值时，拓扑不变量陈数发生变化，表明材料的拓扑性质发生了改变，从一种拓扑态转变为另一种拓扑态。这种拓扑相变伴随着电子结构的显著变化，如能带的重新分布和狄拉克点的消失或出现新的狄拉克点。同时，压力还会影响ZrTe₂的电学性质，如载流子浓度和迁移率等，这些变化与晶体结构和电子结构的改变密切相关。对于Pt₂HgSe₃，压力对其影响同样显著。随着压力的增加，Pt₂HgSe₃发生了晶体结构相变，从常压下的P3̅m1结构转变为另一种结构，这一相变过程伴随着原子位置的重新排列和键长、键角的变化。在电子结构上，压力导致能带结构中狄拉克点和外尔点的位置和能量发生改变，进而影响了材料的拓扑性质。研究发现，在一定压力范围内，Pt₂HgSe₃能够保持其双拓扑半金属的特性，但当压力进一步增大时，可能会发生拓扑相变，拓扑不变量发生变化。压力还对Pt₂HgSe₃的超导性质产生影响，实验和理论计算表明，在特定压力下，Pt₂HgSe₃可能会出现超导电性，这为研究拓扑与超导的关联提供了重要的实验和理论依据。在研究压力效应时，通过构建不同压力下的晶体结构模型，利用第一性原理计算软件进行模拟，得到了压力与晶体结构、电子结构和拓扑性质之间的定量关系，为深入理解压力对二维拓扑半金属的影响提供了详细的信息。3.3.2电场效应研究电场是另一种重要的外部调控因素，能够对二维拓扑半金属的电子结构和输运性质产生显著的调控作用。对于ZrTe₂，在施加电场的情况下，其电子结构发生明显变化。通过第一性原理计算，发现电场能够改变ZrTe₂的能带结构，使狄拉克点的位置和能量发生移动，狄拉克锥的形状也发生畸变。这是因为电场与电子相互作用，改变了电子的势能，从而影响了电子的能量和波函数。在输运性质方面，电场对ZrTe₂的电导率和载流子迁移率产生影响。随着电场强度的增加，电导率和载流子迁移率呈现出非线性的变化趋势，这是由于电场影响了电子的散射过程和输运路径。通过分析电场作用下电子的散射机制，发现电场可以改变电子与杂质、声子等的散射概率，从而影响输运性质。电场还可能导致ZrTe₂出现一些新的物理现象，如谷极化等，这为其在谷电子学器件中的应用提供了潜在的可能性。对于Pt₂HgSe₃，电场对其电子结构和输运性质的调控作用同样显著。施加电场后，Pt₂HgSe₃的能带结构中狄拉克点和外尔点的性质发生改变，导致其拓扑性质也可能发生变化。通过计算电场作用下的拓扑不变量，发现当电场强度达到一定值时，拓扑不变量会发生变化，表明材料的拓扑性质发生了转变。在输运性质上，电场能够调节Pt₂HgSe₃的电阻和霍尔效应等。随着电场强度的变化，电阻呈现出不同的变化规律，霍尔系数也发生改变，这为利用电场调控Pt₂HgSe₃的电学性质提供了理论依据。通过研究电场效应，揭示了电场与二维拓扑半金属电子结构和输运性质之间的内在联系，为开发基于二维拓扑半金属的新型电子器件提供了理论指导，如利用电场调控实现高性能的场效应晶体管等。在研究电场效应时，采用了在晶体结构模型上施加均匀电场的方法，通过调整电场强度和方向，系统地研究了电场对材料性质的影响，为深入理解电场调控机制提供了详细的数据和理论分析。四、二维拓扑半金属的模型研究4.1二维拓扑半金属的紧束缚模型构建4.1.1模型参数确定紧束缚模型在描述二维拓扑半金属的电子行为时，精确确定模型参数至关重要。对于选定的ZrTe₂和Pt₂HgSe₃等二维拓扑半金属材料，模型参数的确定基于对其晶体结构和电子结构的深入分析。以ZrTe₂为例，其晶体结构为六方晶系，空间群为P6₃/mmc。在构建紧束缚模型时，首先考虑其原子间的相互作用。电子在Zr原子和Te原子附近的运动受原子势场的主导，相邻原子势场的作用视为微扰。通过分析原子轨道的重叠情况，确定跳跃积分参数。由于ZrTe₂中Zr-Te键的存在，电子在Zr和Te原子间的跳跃积分是影响电子输运的重要因素。根据量子力学原理，跳跃积分与原子轨道的重叠程度以及原子间的距离相关。通过对ZrTe₂晶体结构的精确测量和理论计算，确定了Zr-Te原子间的距离，进而利用Slater-Koster机制计算不同原子轨道间的跳跃积分。对于s轨道和p轨道之间的跳跃积分，其表达式与原子间连线方向和轨道的取向有关，通过考虑这些几何因素，精确计算出相应的跳跃积分值。在位能参数的确定则与原子的电子云分布和原子的化学环境有关。在ZrTe₂中，Zr原子的d轨道和Te原子的p轨道对电子的束缚能不同，通过分析原子的电子结构和晶体中的电荷分布，确定了每个原子轨道的在位能。考虑到晶体中电子的共有化运动，对在位能进行适当的修正，以反映电子在不同原子轨道间的转移概率。对于Pt₂HgSe₃，其晶体结构更为复杂，属于三方晶系，空间群为P3̅m1。在确定模型参数时，同样需要考虑原子间的相互作用和电子云分布。由于Pt₂HgSe₃中存在多种原子（Pt、Hg、Se），原子间的键合方式和电子云重叠情况更为复杂。通过第一性原理计算和实验测量相结合的方法，确定了不同原子间的跳跃积分和在位能。利用第一性原理计算得到原子间的电荷密度分布，从而推断出原子轨道的重叠程度，进而计算跳跃积分；通过分析晶体中原子的化学环境和电子云分布，确定在位能参数。在计算过程中，考虑到晶体的对称性对参数的影响，利用对称性原理简化计算过程，确保参数的准确性和可靠性。4.1.2模型验证与优化为了验证紧束缚模型的准确性，将模型计算结果与第一性原理计算结果进行对比。对于ZrTe₂，利用紧束缚模型计算其能带结构，并与基于密度泛函理论的第一性原理计算得到的能带结构进行比较。在比较过程中，重点关注狄拉克点的位置、狄拉克锥的形状以及能带的色散关系。通过对比发现，紧束缚模型在一定程度上能够重现第一性原理计算得到的能带特征，但在某些细节上仍存在差异。狄拉克点的能量和位置在两种计算方法中存在一定的偏差，这可能是由于紧束缚模型在处理原子间相互作用时采用了近似方法，忽略了一些高阶相互作用。为了优化紧束缚模型，对模型参数进行调整。通过改变跳跃积分和在位能的数值，观察能带结构的变化，使其与第一性原理计算结果更加接近。采用最小二乘法等优化算法，以第一性原理计算结果为目标函数，对模型参数进行优化，使得紧束缚模型计算得到的能带结构与第一性原理计算结果的误差最小化。在优化过程中，还考虑了模型的物理意义，确保参数的调整符合实际的物理情况。除了能带结构，还对比了两种方法计算得到的态密度、电荷密度等物理量，进一步验证模型的准确性和优化效果。通过优化，紧束缚模型能够更准确地描述ZrTe₂的电子结构，为后续的电子输运性质研究提供了可靠的基础。对于Pt₂HgSe₃，同样进行了模型验证与优化。将紧束缚模型计算得到的能带结构、拓扑不变量等与第一性原理计算结果进行对比，发现模型在描述复杂的能带结构和拓扑性质时存在一定的局限性。为了优化模型，引入了更多的原子轨道和相互作用项，以更全面地考虑电子的行为。通过调整模型参数和改进模型结构，使得紧束缚模型能够更好地拟合第一性原理计算结果，提高了模型对Pt₂HgSe₃电子结构和拓扑性质的描述能力。4.2基于模型的电子输运性质研究4.2.1电导率计算利用紧束缚模型计算二维拓扑半金属的电导率，能够深入理解其电子输运机制。对于ZrTe₂，根据Landauer-Büttiker理论，电导率可以通过计算电子在样品中的传输概率得到。在紧束缚模型框架下，电子的传输概率与能带结构和电子散射过程密切相关。通过计算不同能量下电子在ZrTe₂中的波函数和散射矩阵，确定电子的传输概率。考虑到ZrTe₂中存在的杂质和晶格振动等散射机制，采用自洽Born近似方法来处理电子-杂质散射和电子-声子散射。对于电子-杂质散射，假设杂质为点缺陷，通过计算杂质势与电子波函数的相互作用，得到杂质散射矩阵元，进而确定杂质对电子传输概率的影响。对于电子-声子散射，利用声子的色散关系和电声耦合强度，计算电子-声子散射矩阵元，考虑声子的发射和吸收过程对电子传输的影响。通过求解含时薛定谔方程，得到电子在散射过程中的波函数演化，从而计算出电导率。计算结果表明，ZrTe₂的电导率与电子结构和拓扑性质密切相关。在狄拉克点附近，由于电子具有线性色散关系和高迁移率，电导率呈现出与传统金属不同的特性。随着能量的变化，电导率受到电子散射的影响而发生变化。当考虑杂质散射时，电导率会降低，这是因为杂质散射增加了电子的散射概率，阻碍了电子的传输。而电子-声子散射对电导率的影响则与温度有关，在低温下，声子散射较弱，电导率主要受杂质散射影响；在高温下，声子散射增强，对电导率的影响逐渐增大。对于Pt₂HgSe₃，其电导率的计算同样基于紧束缚模型和Landauer-Büttiker理论。由于Pt₂HgSe₃具有复杂的能带结构和拓扑性质，在计算电导率时需要考虑更多的因素。除了电子-杂质散射和电子-声子散射外，还需要考虑其特殊的拓扑表面态对电子输运的影响。由于拓扑表面态的存在，电子在表面的传输具有独特的性质，可能会导致电导率的各向异性。通过计算不同方向上电子在Pt₂HgSe₃中的传输概率，分析电导率的各向异性特征。结果表明，Pt₂HgSe₃在某些方向上的电导率明显高于其他方向，这与拓扑表面态的电子传输特性有关，为其在电子学器件中的应用提供了重要的理论依据。4.2.2热导率计算计算二维拓扑半金属的热导率，有助于研究其热输运性质，探讨拓扑性质对热导率的影响。在ZrTe₂中，热导率由电子和声子共同贡献。对于电子热导率，根据Wiedemann-Franz定律，电子热导率与电导率和温度成正比，通过前面计算得到的电导率，结合温度和洛伦兹数，可以计算出电子热导率。对于声子热导率，采用声子玻尔兹曼输运方程进行计算。在紧束缚模型基础上，考虑声子的色散关系、声子-声子散射以及声子与电子的相互作用。首先，通过晶格动力学计算得到ZrTe₂的声子色散关系，确定声子的频率和波矢。然后，利用格林函数方法计算声子-声子散射的散射率，考虑三声子散射过程中的能量和动量守恒。对于声子与电子的相互作用，通过电声耦合矩阵元计算电子对声子的散射率，分析电子对声子热输运的影响。通过求解声子玻尔兹曼输运方程，得到声子的分布函数，进而计算出声子热导率。计算结果表明，在ZrTe₂中，电子和声子对热导率的贡献在不同温度下有所不同。在低温下，电子热导率占主导地位，这是因为低温下声子的激发较少，声子热导率较低；随着温度升高，声子热导率逐渐增加，在高温下，声子热导率可能超过电子热导率。拓扑性质对ZrTe₂的热导率也有影响，由于狄拉克点附近电子的特殊性质，电子热导率在狄拉克点附近可能会出现异常变化。此外，声子的散射过程也受到拓扑性质的影响，拓扑保护的表面态可能会影响声子的散射路径，从而对声子热导率产生影响。对于Pt₂HgSe₃，热导率的计算同样考虑电子和声子的贡献。由于其复杂的晶体结构和拓扑性质，声子的色散关系和声子-声子散射过程更为复杂。在计算过程中，利用第一性原理计算得到的晶体结构和电子结构信息，精确计算声子的色散关系和声子-声子散射率。考虑到Pt₂HgSe₃中的双拓扑半金属特性，分析拓扑性质对电子和声子热输运的综合影响。结果发现，Pt₂HgSe₃的热导率在不同温度和方向上呈现出复杂的变化规律，这与它的拓扑性质以及电子和声子的相互作用密切相关，为研究其在热管理等领域的应用提供了理论基础。4.3模型预测与实验验证4.3.1新物性预测基于紧束缚模型，对二维拓扑半金属可能存在的新物性进行预测，为实验研究提供理论指导。对于ZrTe₂，通过对其电子结构和输运性质的深入研究，预测可能存在拓扑磁电效应。在紧束缚模型中，考虑磁场和电场对电子的作用，通过计算电子在电磁场中的哈密顿量和波函数，分析电子的运动轨迹和能量变化。理论分析表明，当ZrTe₂处于外加磁场和电场中时，由于其拓扑性质和电子的相对论性特性，会产生拓扑磁电效应。在这种效应下，材料的电极化强度和磁化强度之间存在耦合关系，即外加电场会引起材料的磁化，而外加磁场则会导致材料的极化。这种拓扑磁电效应的产生源于ZrTe₂中狄拉克点附近电子的特殊性质，电子的运动受到磁场和电场的共同作用，使得电子的波函数发生变化，从而导致电极化和磁化的耦合。预测ZrTe₂在极低温和强磁场条件下可能出现量子振荡现象。在紧束缚模型中，考虑电子在强磁场下的朗道量子化，通过计算电子在朗道能级间的跃迁概率和能量变化，分析量子振荡的特性。当磁场强度变化时，朗道能级会发生移动，电子在不同朗道能级间的跃迁会导致物理量的周期性振荡，如电导率、磁化率等。这种量子振荡现象可以作为研究ZrTe₂拓扑性质和电子结构的重要手段，通过测量量子振荡的频率和幅度，可以获取材料的费米面信息和电子有效质量等参数。对于Pt₂HgSe₃，基于紧束缚模型预测其可能存在的新型超导态。考虑到Pt₂HgSe₃的双拓扑半金属特性和电子-电子相互作用，通过计算电子配对的可能性和超导能隙的大小，分析超导态的形成机制。在紧束缚模型中，引入电子-电子相互作用项，利用BCS理论框架，研究电子在不同动量空间中的配对情况。预测结果表明，在一定条件下，Pt₂HgSe₃中的电子可能会形成库珀对，从而出现超导态，这为实验上探索新型超导材料提供了理论依据。4.3.2与实验结果对比分析将模型预测结果与实验数据进行对比，是检验模型可靠性和局限性的重要手段。对于ZrTe₂，在拓扑磁电效应的预测方面，目前相关实验研究相对较少，但已有一些初步的实验尝试。部分实验通过在ZrTe₂样品上施加电场和磁场，测量其电极化和磁化的变化，但由于实验条件的限制和测量精度的问题，实验结果与理论预测之间存在一定的差异。理论预测的拓扑磁电效应的强度和变化规律在实验中未能完全得到验证，这可能是由于实验样品的质量、杂质含量以及测量方法等因素的影响。未来需要进一步优化实验条件，提高测量精度，以更准确地验证理论预测。在量子振荡现象的研究中，已有一些实验在ZrTe₂中观察到了类似量子振荡的信号。实验通过在低温和强磁场下测量ZrTe₂的电导率和磁化率等物理量，发现随着磁场强度的变化，这些物理量出现了周期性的振荡。与理论预测相比，实验观测到的量子振荡频率和幅度与理论计算结果在一定程度上相符，但也存在一些偏差。理论计算中假设的电子散射机制和样品的理想状态与实际实验情况存在差异，实际样品中存在的杂质和缺陷会影响电子的散射过程，从而导致量子振荡信号的变化。通过进一步研究杂质和缺陷对量子振荡的影响，对理论模型进行修正，能够提高理论与实验的一致性。对于Pt₂HgSe₃的新型超导态预测，目前实验上尚未有明确的观测结果。这可能是由于实现预测的超导条件较为苛刻，实验制备高质量的Pt₂HgSe₃样品以及达到所需的低温和压力等条件存在困难。未来需要进一步改进实验技术，优化样品制备方法，以验证理论预测的新型超导态的存在。通过与实验结果的对比分析，不断完善紧束缚模型，提高其对二维拓扑半金属物理性质的预测能力，为实验研究提供更准确的理论指导。五、三维拓扑半金属的第一性原理研究5.1三维拓扑半金属材料的选择与结构建模5.1.1材料体系介绍本研究选取了MoTe₂和WTe₂作为典型的三维拓扑半金属材料进行深入研究。MoTe₂具有丰富的晶体结构和独特的物理性质，常见的晶体结构包括1T'相和Td相。在1T'相的MoTe₂晶体结构中，其空间群为C2/m，原子排列形成层状结构，每一层由Mo原子和Te原子组成，Mo原子被六个Te原子以八面体配位方式包围。这种结构赋予了MoTe₂一定的各向异性，在电学性质方面，1T'相的MoTe₂表现出金属性，且理论研究表明它具有非平庸的拓扑性质，其能带结构中存在狄拉克点，电子在狄拉克点附近呈现出线性色散关系，具有拓扑半金属的特征。当温度低于250K时，1T'相的MoTe₂会发生晶体结构相变，转变为Td相，Td相的MoTe₂空间群为Pnma，在该相中可以存在第二类外尔费米子，展现出更为丰富的拓扑物理性质。WTe₂同样具有独特的晶体结构和物理特性，其晶体结构属于单斜晶系，空间群为P2₁/c。在WTe₂的晶体结构中，W原子和Te原子通过共价键相互连接，形成了特定的原子排列方式。这种结构使得WTe₂具有较强的各向异性，在电学性质上，WTe₂表现出金属性，并且具有巨大的非饱和磁阻特性，这与它的电子结构密切相关。研究发现，WTe₂是一种第二类外尔半金属，其能带结构中存在外尔点，外尔点周围的电子具有特殊的色散关系和拓扑性质，使得WTe₂在拓扑半金属的研究中具有重要地位。5.1.2结构优化与稳定性分析利用基于密度泛函理论的第一性原理计算方法，借助VASP软件对MoTe₂和WTe₂的原子结构进行细致优化。在计算过程中，选用广义梯度近似（GGA）来处理交换关联能，以平面波基组展开电子波函数，并精心设置合适的平面波截断能和k点网格，从而保障计算结果的准确性与收敛性。通过结构优化，成功获取了材料的平衡晶格常数、原子坐标以及晶胞体积等关键参数。对于MoTe₂的1T'相，优化后的晶格常数a、b、c以及角度β等参数与实验值或其他理论计算值具有良好的一致性，这充分验证了计算方法的可靠性。在优化过程中，通过最小化体系的总能量，促使原子逐步调整至最稳定的位置，进而得到稳定的晶体结构。为了深入分析结构的稳定性，计算了材料的形成能和结合能。形成能是衡量从组成元素的稳定单质状态形成化合物时的能量变化，其计算公式为E_{form}=E_{total}-\sum_{i}n_{i}E_{i}，其中E_{total}是化合物的总能量，n_{i}是组成元素i的原子数，E_{i}是元素i的稳定单质的能量。形成能的负值越大，表明化合物越稳定。对于MoTe₂，计算得到的形成能为负值，这表明其在热力学上是稳定的，能够稳定存在。结合能则是指将晶体中的原子完全分离所需的能量，它反映了原子之间的相互作用强度，结合能越大，晶体结构越稳定。通过计算MoTe₂和WTe₂的结合能，发现它们的结合能都较大，这表明这两种材料的原子之间具有较强的相互作用，晶体结构较为稳定。此外，还通过计算声子谱来进一步验证结构的稳定性。若声子谱中不存在虚频，说明晶体结构在动力学上是稳定的。计算结果表明，MoTe₂和WTe₂的声子谱中均无虚频，进一步证实了它们的结构稳定性。在结构优化过程中，观察到体系能量随着原子位置的调整逐渐降低，最终收敛到一个稳定的值，这也直观地表明了结构优化的有效性和结构的稳定性。5.2电子结构与拓扑性质计算5.2.1能带结构与态密度计算通过第一性原理计算，获得了MoTe₂和WTe₂精确的能带结构和态密度。对于MoTe₂的1T'相，在费米能级附近，能带结构呈现出显著的特征，存在线性色散的能带交叉，进而形成了狄拉克锥。狄拉克锥的出现有力地表明MoTe₂具有拓扑半金属的特性，电子在狄拉克点附近表现出类似于无质量相对论粒子的行为。通过对能带结构的深入分析，准确确定了狄拉克点的位置和能量，以及狄拉克锥的色散关系。同时，计算了MoTe₂的态密度，发现态密度在狄拉克点处呈现出尖锐的峰，这与能带结构中狄拉克锥的特征高度吻合。态密度的分布清晰地反映了电子在不同能量状态下的分布情况，在费米能级附近，态密度主要由Mo原子的d轨道和Te原子的p轨道贡献，这表明这些轨道在电子的输运和拓扑性质中起着至关重要的作用。当MoTe₂从1T'相转变为Td相后，能带结构发生明显变化，出现了外尔点，外尔点周围的能带呈现出独特的线性色散关系，对应态密度也出现新的特征峰，这些变化与Td相中的第二类外尔费米子密切相关。WTe₂的能带结构同样展现出独特的性质，由于其作为第二类外尔半金属的特性，能带结构中存在外尔点。通过仔细分析能带结构，确定了外尔点的位置、能量以及手性等关键特征。计算得到的态密度也显示出丰富的结构，在费米能级附近，态密度存在多个峰和谷，这些特征与能带结构中的外尔点紧密相关。通过对态密度的分波分析，明确了W原子和Te原子的不同轨道对态密度的贡献，发现W原子的d轨道和Te原子的p轨道在费米能级附近的态密度贡献较大，这对深入理解WTe₂的电子结构和拓扑性质具有重要意义。在计算能带结构和态密度时，采用了高精度的计算参数和方法，通过改变k点网格的密度和平面波截断能等参数进行严格的收敛性测试，发现当k点网格密度达到一定值，平面波截断能足够大时，计算结果基本收敛，能够准确反映材料的电子结构特征。5.2.2Berry曲率与拓扑不变量计算为了准确确定MoTe₂和WTe₂的拓扑性质，对它们的Berry曲率和拓扑不变量进行了精确计算。对于MoTe₂的1T'相，首先计算了其在动量空间中的Berry曲率分布。Berry曲率在狄拉克点附近呈现出明显的峰值，这表明狄拉克点周围存在着非平凡的拓扑结构。通过对Berry曲率在整个布里渊区上的积分，得到了拓扑不变量陈数。计算结果表明，MoTe₂的1T'相陈数不为零，这证实了它具有非平庸的拓扑性质，属于拓扑半金属。陈数的非零值意味着MoTe₂存在受拓扑保护的边界态，这些边界态具有独特的电子结构和输运性质，为其在拓扑电子学中的应用提供了坚实基础。当MoTe₂转变为Td相后，通过对动量空间中不同方向上的Berry曲率进行分析，结合对称性分析，确定其拓扑不变量，进一步明确了其作为第二类外尔半金属的拓扑特性。对于WTe₂，由于其复杂的能带结构和外尔点的存在，计算拓扑不变量相对复杂。除了考虑陈数外，还需要深入分析与外尔点相关的拓扑特征。通过对动量空间中不同方向上的Berry曲率进行详细分析，结合对称性分析，确定了WTe₂的拓扑不变量。计算结果表明，WTe₂具有非平凡的拓扑性质，其拓扑不变量的值与理论预测的第二类外尔半金属的特征相符。在计算过程中，充分利用了对称性的性质来简化计算，通过对晶体结构的对称性分析，确定了一些特殊的动量空间点和路径，在这些点和路径上计算Berry曲率和拓扑不变量，不仅减少了计算量，还显著提高了计算结果的准确性。通过计算Berry曲率和拓扑不变量，明确了MoTe₂和WTe₂的拓扑性质，为进一步研究它们的物理性质和应用提供了关键依据。5.3外界因素对三维拓扑半金属性质的影响5.3.1压力效应研究压力作为一种有效的外部调控手段，能够显著改变三维拓扑半金属的晶体结构、电子结构和拓扑性质。对于MoTe₂，通过第一性原理计算深入研究了压力对其晶体结构的影响。随着压力的逐渐增加，MoTe₂的晶格常数逐渐减小，原子间的距离被压缩，层间距也相应减小。在电子结构方面，压力导致MoTe₂的能带结构发生明显变化。狄拉克点的位置和能量发生移动，狄拉克锥的形状也有所改变，这是因为压力改变了原子间的相互作用，进而影响了电子的能量和波函数。在拓扑性质上，压力会诱导MoTe₂发生拓扑相变。当压力达到一定值时，拓扑不变量陈数发生变化，表明材料的拓扑性质发生了改变，从一种拓扑态转变为另一种拓扑态。这种拓扑相变伴随着电子结构的显著变化，如能带的重新分布和狄拉克点的消失或出现新的狄拉克点。同时，压力还会影响MoTe₂的电学性质，如载流子浓度和迁移率等，这些变化与晶体结构和电子结构的改变密切相关。当MoTe₂处于Td相时，压力还可能影响外尔点的性质，导致外尔半金属相的变化。对于WTe₂，压力对其影响同样显著。随着压力的增加，WTe₂的晶体结构发生变化，原子的位置和键长、键角都有所改变。在电子结构上，压力导致能带结构中外尔点的位置和能量发生改变，进而影响了材料的拓扑性质。研究发现，在一定压力范围内，WTe₂能够保持其第二类外尔半金属的特性，但当压力进一步增大时，可能会发生拓扑相变，拓扑不变量发生变化。压力还对WTe₂的超导性质产生影响，实验和理论计算表明，在特定压力下，WTe₂可能会出现超导电性，这为研究拓扑与超导的关联提供了重要的实验和理论依据。在研究压力效应时，通过构建不同压力下的晶体结构模型，利用第一性原理计算软件进行模拟，得到了压力与晶体结构、电子结构和拓扑性质之间的定量关系，为深入理解压力对三维拓扑半金属的影响提供了详细的信息。5.3.2磁场效应研究磁场是另一个重要的外部调控因素，能够对三维拓扑半金属的电子结构和输运性质产生显著的调控作用。对于MoTe₂，在施加磁场的情况下，其电子结构发生明显变化。通过第一性原理计算，发现磁场能够改变MoTe₂的能带结构，使狄拉克点的位置和能量发生移动，狄拉克锥的形状也发生畸变。这是因为磁场与电子相互作用，改变了电子的势能，从而影响了电子的能量和波函数。在输运性质方面，磁场对MoTe₂的电导率和载流子迁移率产生影响。随着磁场强度的增加，电导率和载流子迁移率呈现出非线性的变化趋势，这是由于磁场影响了电子的散射过程和输运路径。通过分析磁场作用下电子的散射机制，发现磁场可以改变电子与杂质、声子等的散射概率，从而影响输运性质。磁场还可能导致MoTe₂出现一些新的物理现象，如反常霍尔效应等，这为其在磁电器件中的应用提供了潜在的可能性。当MoTe₂处于Td相时，磁场对外尔半金属的输运性质影响更为复杂，可能导致手性磁效应等现象的变化。对于WTe₂，磁场对其电子结构和输运性质的调控作用同样显著。施加磁场后，WTe₂的能带结构中外尔点的性质发生改变，导致其拓扑性质也可能发生变化。通过计算磁场作用下的拓扑不变量，发现当磁场强度达到一定值时，拓扑不变量会发生变化，表明材料的拓扑性质发生了转变。在输运性质上，磁场能够调节WTe₂的电阻和霍尔效应等。随着磁场强度的变化，电阻呈现出不同的变化规律，霍尔系数也发生改变，这为利用磁场调控WTe₂的电学性质提供了理论依据。通过研究磁场效应，揭示了磁场与三维拓扑半金属电子结构和输运性质之间的内在联系，为开发基于三维拓扑半金属的新型磁电器件提供了理论指导，如利用磁场调控实现高性能的磁传感器等。在研究磁场效应时，采用了在晶体结构模型上施加均匀磁场的方法，通过调整磁场强度和方向，系统地研究了磁场对材料性质的影响，为深入理解磁场调控机制提供了详细的数据和理论分析。六、三维拓扑半金属的模型研究6.1三维拓扑半金属的紧束缚模型构建6.1.1模型参数确定对于三维拓扑半金属MoTe₂和WTe₂，构建紧束缚模型时，准确确定模型参数至关重要。以MoTe₂为例，其1T'相晶体结构为层状，原子间的相互作用较为复杂。在确定跳跃积分参数时，需考虑Mo原子与Te原子之间不同轨道的相互作用。由于Mo原子的d轨道与Te原子的p轨道存在较强的杂化作用，电子在这些轨道间的跳跃积分对电子结构和输运性质有显著影响。通过分析晶体结构中原子间的距离和轨道取向，利用Slater-Koster方法计算不同原子轨道间的跳跃积分。对于最近邻的Mo-Te原子，s-p轨道间的跳跃积分t_{sp}与原子间距离r_{Mo-Te}以及轨道夹角有关，根据量子力学原理，t_{sp}\propto\frac{1}{r_{Mo-Te}^3}\cos\theta，其中\theta为轨道夹角。通过精确测量MoTe₂晶体结构中原子的坐标和键长，确定r_{Mo-Te}和\theta的值，进而计算出t_{sp}。同时，考虑次近邻原子间的相互作用，对跳跃积分进行修正，以更准确地描述电子的行为。在位能参数的确定与原子的电子云分布和化学环境密切相关。Mo原子在MoTe₂晶体中处于特定的配位环境，其d轨道的在位能受到周围Te原子的影响。通过分析Mo原子的电子结构和晶体中的电荷分布，确定其d轨道的在位能\epsilon_d。考虑到晶体中电子的共有化运动，对在位能进行修正，以反映电子在不同原子轨道间的转移概率。对于Te原子，同样根据其电子结构和在晶体中的位置，确定其p轨道的在位能\epsilon_p。在确定这些参数时，还需考虑温度等因素对原子电子云分布的影响，以提高模型的准确性。对于WTe₂，其晶体结构为单斜晶系，原子排列方式与MoTe₂不同，原子间的相互作用也有所差异。在确定跳跃积分参数时，需要考虑W原子与Te原子间的多种轨道组合的相互作用。由于WTe₂中存在较强的各向异性，不同方向上原子间的跳跃积分不同。通过分析晶体结构在不同方向上的原子排列和键长，分别计算不同方向上的跳跃积分。对于沿a轴方向的W-Te原子间的跳跃积分t_{a}，与沿b轴方向的跳跃积分t_{b}，由于原子间距离和轨道取向的不同，其值也不同。通过精确的结构分析和量子力学计算，确定这些跳跃积分的值。在位能参数的确定同样考虑W原子和Te原子的电子结构和化学环境，确定W原子d轨道和Te原子p轨道的在位能，并根据晶体的对称性和电子的共有化运动进行修正。在确定模型参数过程中，充分利用第一性原理计算结果和实验测量数据，相互验证和补充，以确保参数的准确性和可靠性。6.1.2模型验证与优化为验证紧束缚模型的准确性，将模型计算结果与第一性原理计算结果进行对比。以MoTe₂为例，利用紧束缚模型计算其能带结构，并与基于密度泛函理论的第一性原理计算得到的能带结构进行细致比较。在比较过程中，重点关注狄拉克点和外尔点（1T'相主要关注狄拉克点，Td相关注外尔点）的位置、狄拉克锥和外尔点附近能带的形状以及能带的色散关系。通过对比发现，紧束缚模型在一定程度上能够重现第一性原理计算得到的能带特征，但在某些细节上仍存在差异。狄拉克点的能量和位置在两种计算方法中存在一定的偏差，这可能是由于紧束缚模型在处理原子间相互作用时采用了近似方法，忽略了一些高阶相互作用和电子-电子关联效应。为优化紧束缚模型，对模型参数进行调整。采用最小二乘法等优化算法，以第一性原理计算结果为目标函数，对跳跃积分和在位能等参数进行优化，使得紧束缚模型计算得到的能带结构与第一性原理计算结果的误差最小化。在优化过程中，还考虑了模型的物理意义，确保参数的调整符合实际的物理情况。除了能带结构，还对比了两种方法计算得到的态密度、电荷密度等物理量，进一步验证模型的准确性和优化效果。通过优化，紧束缚模型能够更准确地描述MoTe₂的电子结构，为后续的电子输运性质研究提供了可靠的基础。对于WTe₂，同样进行了模型验证与优化。将紧束缚模型计算得到的能带结构、拓扑不变量等与第一性原理计算结果进行对比，发现模型在描述复杂的能带结构和拓扑性质时存在一定的局限性。为了优化模型，引入了更多的原子轨道和相互作用项，以更全面地考虑电子的行为。考虑到WTe₂中电子-声子相互作用对电子结构的影响，在模型中加入电子-声子相互作用项，通过调整相关参数，使得紧束缚模型能够更好地拟合第一性原理计算结果，提高了模型对WTe₂电子结构和拓扑性质的描述能力。6.2基于模型的电子输运性质研究6.2.1磁电阻计算利用紧束缚模型计算三维拓扑半金属的磁电阻，能够深入探究其电子输运与拓扑性质的关联。对于MoTe₂，根据Landauer-Büttiker理论，磁电阻与电子在磁场中的传输概率密切相关。在紧束缚模型框架下，考虑磁场对电子的洛伦兹力作用，电子在磁场中的哈密顿量发生变化。通过引入Peierls代换，将磁场效应纳入哈密顿量中，从而计算电子在磁场中的波函数和散射矩阵。在计算过程中，考虑MoTe₂中存在的杂质和晶格振动等散射机制。对于杂质散射，采用自洽Born近似方法，计算杂质势与电子波函数的相互作用，得到杂质散射矩阵元，确定杂质对电子传输概率的影响。对于晶格振动引起的电子-声子散射，利用声子的色散关系和电声耦合强度，计算电子-声子散射矩阵元，考虑声子的发射和吸收过程对电子传输的影响。通过求解含时薛定谔方程，得到电子在散射过程中的波函数演化，进而计算出不同磁场强度下的磁电阻。计算结果表明，MoTe₂的磁电阻与电子结构和拓扑性质紧密相关。在狄拉克点附近，由于电子具有线性色散关系和高迁移率，磁电阻呈现出与传统金属不同的特性。当磁场强度变化时，磁电阻出现非线性变化，这是因为磁场改变了电子的运动轨迹和散射概率。在低温下，杂质散射对磁电阻的影响较大；随着温度升高，电子-声子散射逐渐增强，对磁电阻的贡献增大。此外，拓扑性质对磁电阻也有显著影响，由于MoTe₂具有非平凡的拓扑性质，存在受拓扑保护的边界态，这些边界态的电子输运特性会影响磁电阻的大小和变化规律。对于WTe₂，其磁电阻的计算同样基于紧束缚模型和Landauer-Büttiker理论。由于WTe₂是第二类外尔半金属，具有独特的能带结构和拓扑性质，在计算磁电阻时需要考虑更多的因素。除了电子-杂质散射和电子-声子散射外，还需要考虑外尔点附近电子的特殊行为以及拓扑表面态对电子输运的影响。在强磁场下，外尔点附近的电子会发生朗道量子化，形成手性朗道能级，这会导致磁电阻出现量子振荡现象。通过计算电子在不同朗道能级间的跃迁概率和能量变化，分析量子振荡的特性，研究其与外尔点拓扑电荷和手性的关系。结果表明，WTe₂的磁电阻在强磁场下呈现出复杂的变化规律，与外尔半金属的拓扑性质密切相关，为研究外尔半金属的电子输运性质提供了重要的理论依据。6.2.2霍尔效应计算计算三维拓扑半金属的霍尔效应，有助于深入理解其拓扑性质对输运性质的影响，为拓扑量子器件的应用提供理论支持。对于MoTe₂，根据霍尔效应的原理，霍尔电压与电流、磁场以及材料的电子结构有关。在紧束缚模型中，通过计算电子在磁场中的运动轨迹和散射过程，确定霍尔系数。考虑电子在磁场中的洛伦兹力作用，电子在垂直于电流和磁场的方向上发生偏转，从而在样品两侧产生霍尔电压。利用紧束缚模型计算电子的波函数和散射矩阵，结合玻尔兹曼输运方程，计算不同磁场强度和电流密度下的霍尔系数。在计算过程中，考虑MoTe₂中杂质和晶格振动对电子散射的影响，以及电子之间的相互作用。通过自洽Born近似方法处理杂质散射，利用声子的色散关系和电声耦合强度计算电子-声子散射，考虑电子-电子相互作用对电子输运的影响。计算结果表明，MoTe₂的霍尔系数与拓扑性质密切相关。由于MoTe₂具有非平凡的拓扑性质，存在受拓扑保护的边界态，这些边界态的电子对霍尔效应有重要贡献。在狄拉克点附近，霍尔系数呈现出与传统金属不同的行为，这是由于狄拉克费米子的特殊性质导致的。随着磁场强度的变化，霍尔系数发生变化，这与电子在磁场中的运动状态和散射概率的改变有关。此外，温度对霍尔系数也有影响，在低温下，杂质散射主导霍尔系数的变化；随着温度升高，电子-声子散射增强，霍尔系数的温度依赖性发生变化。对于WTe₂，其霍尔效应的计算同样考虑其特殊的能带结构和拓扑性质。由于WTe₂是第二类外尔半金属，外尔点的存在使得霍尔效应更为复杂。在计算霍尔系数时，需要考虑外尔点附近电子的手性和拓扑表面态的影响。外尔点附近的电子具有手性，其运动方向与手性相关，这会导致霍尔效应出现一些特殊的性质。拓扑表面态的电子输运特性也会影响霍尔系数的大小和变化规律。通过计算电子在不同动量空间中的运动轨迹和散射过程，分析外尔点和拓扑表面态对霍尔效应的影响机制。结果表明，WTe₂的霍尔系数在不同磁场和温度下呈现出复杂的变化规律，与外尔半金属的拓扑性质密切相关，为研究基于WTe₂的拓扑量子器件提供了理论基础。6.3模型预测与实验验证6.3.1新物性预测基于紧束缚模型，对三维拓扑半金属可能存在的新物性进行预测，为实验研究提供理论指导。对于MoTe₂，通过对其电子结构和输运性质的深入研究，预测可能存