基于第一性原理探究石墨烯氮化硼及Ⅵ族硫化物异质结热输运特性

基于第一性原理探究石墨烯氮化硼及Ⅵ族硫化物异质结热输运特性一、绪论1.1研究背景与意义在当今信息化时代，电子芯片作为各类电子设备的核心部件，其性能的提升对于推动信息技术的发展至关重要。随着半导体微电子工业的不断进步，芯片的集成度持续提高，尺寸不断缩小，这使得单位面积上的热通量急剧增加，散热问题成为制约芯片性能提升的关键瓶颈。按照Arrhenius方程，在半导体芯片的工作温度范围内，芯片温度每上升10℃，其寿命就会降低50%。例如，在高性能计算机中，大量的晶体管密集工作，产生的热量若不能及时散发，不仅会导致芯片性能下降，出现运行速度变慢、数据处理错误等问题，严重时甚至会烧毁芯片，使设备无法正常运行。为了解决芯片的散热问题，深入研究材料的热输运性质显得尤为重要。低维材料，如石墨烯、氮化硼及Ⅵ族硫化物等，由于其独特的原子结构和电子特性，展现出了优异的热学性能，成为了研究的热点。其中，石墨烯具有超高的热导率，理论值可达5000W/(m・K)以上，这使得它在热管理领域具有巨大的应用潜力，有望作为高效的散热材料用于芯片散热。氮化硼则具有良好的绝缘性和较高的热导率，能够在保证电气隔离的同时有效地传导热量，可用于制作芯片的绝缘散热层。Ⅵ族硫化物如二硫化钼、二硫化钨等，也因其独特的热学和电学性质，在电子器件中具有潜在的应用价值。当这些低维材料形成异质结时，由于不同材料之间的界面相互作用，其热输运性质会发生显著变化。异质结中的界面热阻、声子散射等因素都会对热输运过程产生重要影响，进而影响整个器件的热性能。因此，研究石墨烯、氮化硼及Ⅵ族硫化物异质结的热输运性质，对于理解低维材料体系中的热传导机制，开发新型的高性能热管理材料和器件具有重要的理论意义和实际应用价值。从理论层面来看，研究异质结的热输运性质有助于深入揭示低维材料体系中声子的传播、散射以及界面热阻的微观机制。通过第一性原理计算等方法，可以从原子和电子层面详细分析异质结中不同材料之间的相互作用，以及这些作用对热输运过程的影响。这不仅能够丰富和完善低维材料热学理论，还为进一步优化材料的热性能提供了理论基础。在实际应用方面，对异质结热输运性质的深入了解，能够为新型热管理材料和器件的设计提供指导。例如，在芯片散热领域，可以根据异质结的热输运特性，设计出更高效的散热结构，提高芯片的散热效率，降低芯片温度，从而提升芯片的性能和可靠性。此外，在其他需要高效热管理的领域，如新能源汽车、航空航天等，这些研究成果也具有重要的应用前景。通过开发基于异质结的高性能热管理材料和器件，可以有效地解决这些领域中的散热问题，推动相关技术的发展和进步。1.2二维材料概述二维材料是指电子仅可在两个维度的平面内自由运动的材料，其原子层厚度在纳米量级，通常只有几个原子层甚至单原子层。这种独特的原子结构赋予了二维材料许多优异的物理性质，使其在众多领域展现出巨大的应用潜力。以下将对本研究中涉及的石墨烯、六方氮化硼、过渡性金属硫化物（如二硫化钼、二硫化钨）等二维材料的结构和基本特性进行详细介绍。1.2.1石墨烯石墨烯是由碳原子以sp²杂化轨道组成六角型呈蜂巢晶格的二维碳纳米材料，其结构如图1.1(a)所示。它是一种由单个原子层构成的平面薄膜，具有极高的晶体质量和原子级平整度。石墨烯中的碳原子之间通过共价键相互连接，形成了稳定的六角形晶格结构，这种结构赋予了石墨烯许多优异的物理性质。在力学性能方面，石墨烯具有出色的强度和柔韧性。理论计算表明，石墨烯的拉伸强度高达130GPa，是钢铁的数百倍，这使得它在柔性电子器件和高强度复合材料等领域具有潜在的应用价值。例如，在可穿戴电子设备中，石墨烯可以作为柔性基板，为其他电子元件提供支撑，同时还能承受一定程度的弯曲和拉伸，保证设备的正常运行。石墨烯的电学性能也十分优异，它是一种零带隙的半金属材料，具有极高的载流子迁移率，室温下的迁移率可达到2×10⁵cm²/(V・s)，这使得石墨烯在高速电子器件和高频电路等方面具有广阔的应用前景。例如，利用石墨烯制作的场效应晶体管，其开关速度比传统的硅基晶体管更快，能够大大提高集成电路的运行速度。此外，石墨烯还具有良好的导电性，其电导率可达10⁶S/m，可用于制作高性能的电极材料，如在超级电容器中，石墨烯电极能够提高电容器的充放电效率和循环稳定性。在热学性能上，石墨烯拥有超高的热导率，理论值可达5000W/(m・K)以上，这使其成为热管理领域的理想材料。在电子芯片中，石墨烯可以作为散热材料，将芯片产生的热量快速传导出去，有效降低芯片温度，提高芯片的性能和可靠性。例如，将石墨烯与其他材料复合制备成散热薄膜，应用于智能手机、电脑等电子设备的散热模块中，能够显著提高设备的散热效率，延长设备的使用寿命。【配图1张：(a)石墨烯结构示意图；(b)六方氮化硼结构示意图；(c)二硫化钼结构示意图；(d)二硫化钨结构示意图】1.2.2六方氮化硼六方氮化硼（h-BN）的结构与石墨烯类似，同样具有层状结构，如图1.1(b)所示。它由硼原子和氮原子交替排列组成六角形的蜂窝状晶格，层间通过较弱的范德华力相互作用结合在一起。这种结构使得六方氮化硼在保留了类似于石墨烯的一些特性的同时，还具有自身独特的性质。六方氮化硼是一种宽带隙半导体材料，其带隙约为5.9eV，这一特性使其具有良好的绝缘性能，可用于制作电子器件中的绝缘层，有效隔离不同的导电部分，防止电流泄漏，确保电子器件的安全稳定运行。例如，在集成电路中，六方氮化硼可以作为绝缘介质，将不同的晶体管隔开，提高电路的集成度和性能。在热学性能方面，六方氮化硼具有较高的热导率，室温下其热导率可达300-400W/(m・K)，同时还具有较低的热膨胀系数，这使得它在高温环境下能够保持稳定的性能，不易因温度变化而发生变形或损坏。因此，六方氮化硼常被用作高温结构材料和散热材料，在电子封装、航空航天等领域发挥着重要作用。例如，在航空发动机的高温部件中，使用六方氮化硼材料可以提高部件的耐高温性能和热稳定性，保证发动机在恶劣的工作环境下正常运行。此外，由于其良好的绝缘性和高导热性，六方氮化硼还可以用于制作电子设备的散热基板，将芯片产生的热量快速传导出去，提高设备的散热效率。六方氮化硼还具有优异的化学稳定性和润滑性。它在高温、强酸碱等恶劣环境下都能保持稳定的化学性质，不易被腐蚀。同时，其层状结构使得六方氮化硼具有类似于石墨的润滑性能，可作为高温润滑剂使用，在高温、高压等特殊工况下，能够有效降低摩擦系数，减少机械部件的磨损，提高机械设备的运行效率和寿命。1.2.3过渡性金属硫化物（以二硫化钼、二硫化钨为例）过渡性金属硫化物（TMDs）是一类具有代表性的二维材料，其中二硫化钼（MoS₂）和二硫化钨（WS₂）是研究较为广泛的两种材料。它们的结构由过渡金属原子（Mo或W）与硫原子组成，具有类似三明治的层状结构，如图1.1(c)和(d)所示。在每一层中，过渡金属原子位于中间，两侧是硫原子，通过共价键相互连接形成六边形的晶格结构。层与层之间则通过较弱的范德华力相互作用。二硫化钼是一种直接带隙半导体材料，其块体材料的带隙约为1.2eV，而单层二硫化钼的带隙可增加到1.8eV左右，这种带隙的变化使其在半导体器件领域具有重要的应用价值。例如，单层二硫化钼可用于制作场效应晶体管，由于其直接带隙的特性，能够实现更高的开关速度和更低的功耗，有望应用于下一代高性能集成电路中。此外，二硫化钼还具有良好的光学性质，能够吸收和发射特定波长的光，可用于制备光电探测器、发光二极管等光电器件。在光电探测器中，二硫化钼能够将光信号转化为电信号，对微弱的光信号具有较高的响应灵敏度，可应用于光通信、生物医学检测等领域。二硫化钨同样是一种具有层状结构的过渡金属硫化物，其块体为间接带隙半导体，带隙约为1.3eV，单层二硫化钨则转变为直接带隙半导体，带隙约为2.0eV。与二硫化钼类似，二硫化钨也具有独特的电学、光学和力学性能。在电学方面，二硫化钨的载流子迁移率较高，可用于制作高性能的电子器件。在光学方面，二硫化钨对光的吸收和发射特性使其在光电器件中具有潜在的应用，如可用于制备光电传感器、激光发射器等。此外，二硫化钨还具有较好的机械性能，能够在一定程度上承受外力的作用，这为其在柔性电子器件中的应用提供了可能。例如，将二硫化钨制成柔性薄膜，应用于可穿戴电子设备中，能够实现设备的轻量化和柔性化，提高用户的佩戴舒适度。1.3二维异质结研究现状二维异质结是由两种或多种不同的二维材料通过范德华力相互堆叠而成的新型材料体系。由于不同二维材料之间的协同效应和界面相互作用，二维异质结展现出了许多独特的物理性质，在高速电子器件、高效光电器件、高性能传感器等领域具有广阔的应用前景，因此受到了科学界和产业界的广泛关注。在实验制备方面，随着材料制备技术的不断发展，目前已经实现了多种二维异质结的高质量制备。化学气相沉积（CVD）技术是制备二维异质结的常用方法之一。通过精确控制反应气体的流量、温度和生长时间等参数，能够在衬底上逐层生长不同的二维材料，从而实现二维异质结的可控合成。例如，利用CVD技术成功制备了石墨烯/六方氮化硼异质结，这种异质结结合了石墨烯的高导电性和六方氮化硼的良好绝缘性，在电子器件中具有潜在的应用价值，可用于制作高性能的场效应晶体管。机械剥离法也是制备二维异质结的重要手段。该方法通过使用胶带等工具从体材料上剥离出单层或少数层的二维材料，然后将不同的二维材料逐层转移并堆叠在一起，形成二维异质结。这种方法虽然制备效率相对较低，但能够获得高质量的二维材料和异质结，有利于研究其本征物理性质。例如，通过机械剥离法制备的二硫化钼/二硫化钨异质结，在光电器件中表现出了独特的光学和电学性质，可用于制备高性能的光电探测器。此外，分子束外延（MBE）技术也可用于制备高质量的二维异质结，该技术能够在原子尺度上精确控制材料的生长，制备出具有精确原子结构和界面的异质结，为研究异质结的微观结构和物理性质提供了有力的手段。在理论研究方面，第一性原理计算、分子动力学模拟等方法被广泛应用于二维异质结的研究中。第一性原理计算基于量子力学原理，能够从原子和电子层面出发，精确计算二维异质结的电子结构、力学性能、光学性质等。通过第一性原理计算，研究人员深入分析了不同二维材料之间的界面相互作用，揭示了异质结中电荷转移、能带结构变化等微观机制。例如，对于石墨烯/过渡金属硫化物异质结，第一性原理计算表明，界面处存在电荷转移现象，这导致了异质结的能带结构发生变化，进而影响其电学和光学性质。分子动力学模拟则主要用于研究二维异质结的热学性能和动态行为。通过模拟原子的运动轨迹和相互作用，能够得到二维异质结的热导率、热膨胀系数等热学参数，以及在外界扰动下的结构稳定性和动态响应。例如，利用分子动力学模拟研究了二维异质结在不同温度下的热输运过程，发现界面处的声子散射对热导率有显著影响。此外，理论研究还涉及二维异质结的输运性质、磁性等方面，为深入理解其物理性质和开发新型应用提供了理论支持。然而，目前关于二维异质结热输运的理论研究仍存在一些不足。在理论模型方面，虽然现有的第一性原理计算和分子动力学模拟方法能够提供一定的信息，但这些模型往往存在一定的局限性。例如，第一性原理计算在处理大规模体系时计算量巨大，难以模拟长时间尺度的热输运过程；分子动力学模拟中使用的力场参数往往是基于经验或半经验方法确定的，对于复杂的二维异质结体系，力场的准确性可能受到质疑，这会影响模拟结果的可靠性。在实验验证方面，由于二维异质结的热输运性质测量较为困难，实验技术还不够成熟，目前相关的实验数据相对较少，这使得理论研究缺乏足够的实验验证和支撑，导致理论与实验之间的对比和验证存在一定的困难，限制了对二维异质结热输运性质的深入理解和准确把握。1.4研究内容与方法本研究旨在运用第一性原理计算方法，深入探究石墨烯、氮化硼及Ⅵ族硫化物异质结的热输运性质，具体研究内容与方法如下：1.4.1研究内容异质结体系构建：构建石墨烯/氮化硼、石墨烯/二硫化钼、石墨烯/二硫化钨、氮化硼/二硫化钼、氮化硼/二硫化钨、二硫化钼/二硫化钨等多种异质结体系。通过对不同二维材料的合理组合，考虑不同的堆叠方式和界面间距，以获得具有代表性的异质结模型。例如，对于石墨烯/氮化硼异质结，研究AB堆叠和AA堆叠两种方式下异质结的结构和性质差异。在AB堆叠中，石墨烯的碳原子位于氮化硼六角形晶格的中心上方，而在AA堆叠中，石墨烯的碳原子与氮化硼的硼原子或氮原子直接对齐，分析这两种堆叠方式对异质结热输运性质的影响。热输运性质计算：运用基于密度泛函理论的第一性原理计算方法，借助VASP、QuantumEspresso等计算软件，计算异质结的热导率、界面热阻等热输运性质。通过计算声子色散关系、声子态密度和声子散射率等参数，深入分析声子在异质结中的传播和散射机制。例如，在计算石墨烯/二硫化钼异质结的热导率时，通过分析声子色散曲线，确定不同频率声子的传播速度和模式，进而研究声子在界面处的散射情况对热导率的影响。影响因素分析：系统研究界面结构、原子间相互作用、温度等因素对异质结热输运性质的影响。通过改变异质结的界面结构，如引入缺陷、杂质或改变界面粗糙度，分析其对声子散射和热输运的影响。研究不同原子间相互作用强度下异质结的热输运性质变化，以及温度升高时声子散射机制的变化对热导率和界面热阻的影响。例如，在研究界面缺陷对氮化硼/二硫化钨异质结热输运性质的影响时，通过在界面处引入硼原子空位或硫原子空位，计算热输运性质的变化，分析缺陷对声子散射路径和散射概率的影响。1.4.2研究方法第一性原理计算：基于密度泛函理论，将电子的动能、电子-电子相互作用能以及电子与原子核的相互作用能等用电子密度的泛函来表示。通过求解Kohn-Sham方程，得到体系的电子结构和能量，进而计算异质结的各种物理性质。在计算过程中，采用平面波赝势方法，将离子实对电子的作用用赝势来代替，以降低计算量。同时，选择合适的交换关联泛函，如广义梯度近似（GGA）中的PBE泛函，来描述电子之间的交换关联作用，以保证计算结果的准确性。声子计算方法：采用有限位移法计算声子色散关系。通过给原子施加微小的位移，计算体系的能量变化，进而得到力常数矩阵，通过对角化力常数矩阵得到声子的频率和本征矢，从而获得声子色散关系。利用密度泛函微扰理论（DFPT）计算声子态密度，该方法能够在倒空间中高效地计算声子频率和态密度，为分析声子的性质提供了有力的工具。数据处理与分析：对计算得到的数据进行深入分析，运用Origin、Matlab等软件绘制图表，直观展示异质结的热输运性质与各影响因素之间的关系。通过拟合数据，建立数学模型，进一步揭示热输运性质的变化规律。例如，通过对不同温度下异质结热导率数据的拟合，建立热导率与温度的函数关系，分析热导率随温度的变化趋势和机制。二、理论基础与计算方法2.1第一性原理基本理论2.1.1多粒子体系的薛定谔方程在量子力学中，多粒子体系的运动状态由波函数\Psi(r_1,r_2,\cdots,r_N,t)来描述，其中r_i表示第i个粒子的坐标，t为时间。多粒子体系的薛定谔方程是描述其量子行为的核心方程，其形式为：i\hbar\frac{\partial\Psi(r_1,r_2,\cdots,r_N,t)}{\partialt}=\hat{H}\Psi(r_1,r_2,\cdots,r_N,t)其中i是虚数单位，\hbar是约化普朗克常数，\hat{H}为哈密顿算符，它包含了体系中所有粒子的动能以及粒子间的相互作用势能。对于由N个电子和M个原子核组成的体系，哈密顿算符可表示为：\hat{H}=-\sum_{i=1}^{N}\frac{\hbar^2}{2m_e}\nabla_{i}^{2}-\sum_{j=1}^{M}\frac{\hbar^2}{2M_j}\nabla_{j}^{2}+\sum_{i=1}^{N}\sum_{j=1}^{M}\frac{Z_je^2}{|r_i-R_j|}+\frac{1}{2}\sum_{i\neqj}^{N}\frac{e^2}{|r_i-r_j|}+\frac{1}{2}\sum_{k\neql}^{M}\frac{Z_kZ_le^2}{|R_k-R_l|}其中，第一项表示电子的动能，m_e为电子质量，\nabla_{i}^{2}是对第i个电子坐标的拉普拉斯算符；第二项表示原子核的动能，M_j为第j个原子核的质量，\nabla_{j}^{2}是对第j个原子核坐标的拉普拉斯算符；第三项表示电子与原子核之间的库仑吸引能，Z_j是第j个原子核的电荷数，e为电子电荷量，r_i和R_j分别为电子和原子核的位置矢量；第四项表示电子之间的库仑排斥能；第五项表示原子核之间的库仑排斥能。多粒子体系的薛定谔方程从本质上描述了微观粒子的波动性和粒子性，它是量子力学的基本方程之一，在描述微观体系中占据着核心地位。通过求解薛定谔方程，可以得到体系的波函数，进而获得体系的各种物理性质，如能量、电子密度分布等。然而，由于多粒子体系中粒子间相互作用的复杂性，精确求解薛定谔方程是极其困难的，通常需要采用一系列的近似方法来简化计算。2.1.2绝热近似绝热近似，又称Born-Oppenheimer近似，是由MaxBorn和J.RobertOppenheimer在1927年提出的一种重要的近似方法，在简化多粒子体系计算中发挥着关键作用。其基本原理基于原子核和电子质量的巨大差异，原子核的质量远大于电子的质量，大约是电子质量的1836倍。这种质量上的悬殊导致原子核的运动速度比电子慢得多，当电子在原子核形成的势场中快速运动时，原子核几乎可以看作是静止的。因此，可以将电子的运动和原子核的运动分开考虑，把多粒子体系的薛定谔方程分解为电子和原子核两部分的薛定谔方程。具体来说，对于一个包含N个电子和M个原子核的体系，在绝热近似下，先固定原子核的位置，求解电子的薛定谔方程：\left[-\sum_{i=1}^{N}\frac{\hbar^2}{2m_e}\nabla_{i}^{2}+\sum_{i=1}^{N}\sum_{j=1}^{M}\frac{Z_je^2}{|r_i-R_j|}+\frac{1}{2}\sum_{i\neqj}^{N}\frac{e^2}{|r_i-r_j|}\right]\Psi_{e}(r_1,r_2,\cdots,r_N;R_1,R_2,\cdots,R_M)=E_{e}(R_1,R_2,\cdots,R_M)\Psi_{e}(r_1,r_2,\cdots,r_N;R_1,R_2,\cdots,R_M)得到电子的波函数\Psi_{e}和能量E_{e}，这里的能量E_{e}是原子核位置(R_1,R_2,\cdots,R_M)的函数。然后，将电子的能量E_{e}作为原子核运动的势能，求解原子核的薛定谔方程：\left[-\sum_{j=1}^{M}\frac{\hbar^2}{2M_j}\nabla_{j}^{2}+E_{e}(R_1,R_2,\cdots,R_M)\right]\Psi_{n}(R_1,R_2,\cdots,R_M)=E\Psi_{n}(R_1,R_2,\cdots,R_M)从而得到原子核的波函数\Psi_{n}和体系的总能量E。绝热近似的适用范围主要取决于体系中电子和原子核运动的耦合程度。在大多数分子和固体体系中，电子和原子核运动的时间尺度差异明显，绝热近似能够很好地成立。然而，在一些特殊情况下，如分子的光解离过程、固体中的电子-声子强耦合体系等，电子和原子核的运动耦合较强，绝热近似可能不再适用，需要考虑非绝热效应。在研究分子的光激发过程时，电子激发态和基态之间的耦合可能导致电子和原子核的运动不能简单地分离，此时就需要采用非绝热理论来描述体系的动力学过程。2.1.3Hartree-Fock近似Hartree-Fock近似是一种用于处理多电子体系中电子间相互作用的重要方法。在多电子体系中，电子之间存在着复杂的库仑相互作用，这使得精确求解多电子体系的薛定谔方程变得极为困难。Hartree-Fock近似的基本思路是将多电子问题简化为单电子问题，通过自洽场方法来处理电子间的相互作用。该方法假设每个电子都在其他电子的平均势场中独立运动，体系的总波函数可以表示为单电子波函数的Slater行列式，以保证波函数满足费米子的反对称性。对于一个包含N个电子的体系，其总波函数\Psi(r_1,r_2,\cdots,r_N)可写为：\Psi(r_1,r_2,\cdots,r_N)=\frac{1}{\sqrt{N!}}\begin{vmatrix}\psi_1(r_1)&\psi_1(r_2)&\cdots&\psi_1(r_N)\\\psi_2(r_1)&\psi_2(r_2)&\cdots&\psi_2(r_N)\\\vdots&\vdots&\ddots&\vdots\\\psi_N(r_1)&\psi_N(r_2)&\cdots&\psi_N(r_N)\end{vmatrix}其中\psi_i(r_j)表示第i个单电子波函数在第j个电子坐标处的值。基于上述假设，Hartree-Fock方程可以通过变分原理推导得出，其形式为：\left[-\frac{\hbar^2}{2m_e}\nabla_{i}^{2}+V_{ext}(r_i)+V_{H}(r_i)+V_{x}(r_i)\right]\psi_i(r_i)=\epsilon_i\psi_i(r_i)其中，第一项是电子的动能项；V_{ext}(r_i)是外部势场，主要来自原子核的吸引势；V_{H}(r_i)是Hartree势，描述了其他电子对第i个电子的平均库仑排斥作用，可表示为：V_{H}(r_i)=\sum_{j\neqi}^{N}\int\frac{e^2|\psi_j(r_j)|^2}{|r_i-r_j|}dr_jV_{x}(r_i)是交换势，考虑了电子的反对称性导致的交换能，它体现了由于泡利不相容原理，自旋相同的电子之间存在的一种特殊的相互作用。交换势的计算较为复杂，在Hartree-Fock近似中通常采用一些近似方法来处理。\epsilon_i是第i个单电子的能量。Hartree-Fock近似在处理原子数较少的系统时是一种很方便的近似方法，能够给出较为准确的结果。然而，它也存在一些局限性。Hartree-Fock近似完全忽略了电子关联作用，即电子之间除了平均库仑相互作用和交换作用之外的瞬时相互作用。这种忽略使得Hartree-Fock方法在处理一些强关联体系时，如过渡金属化合物、高温超导体等，计算结果与实验值存在较大偏差。随着电子数的增多，Hartree-Fock方法的计算量呈指数增加，对计算机的内存和CPU运算速度要求极高，这使得它在处理大规模多电子体系时面临巨大的挑战。2.2密度泛函理论2.2.1Hobenberg-Kohn定理密度泛函理论（DFT）是一种研究多电子体系电子结构的量子力学方法，在凝聚态物理和计算化学领域有着广泛的应用。Hobenberg-Kohn定理为密度泛函理论提供了坚实的理论基础。Hobenberg-Kohn第一定理指出，对于一个共同的外部势v(r)，相互作用的多粒子系统的所有基态性质都由（非简併）基态的电子密度分布n(r)唯一地决定。这意味着电子密度分布n(r)包含了体系基态的所有信息，是描述多粒子系统基态性质的基本变量。例如，对于一个原子体系，其基态的能量、电子云分布等性质都可以通过基态的电子密度分布n(r)来确定。从物理意义上讲，电子密度分布反映了电子在空间中的概率分布情况，而多粒子系统的各种性质正是由电子的分布和相互作用所决定的。Hobenberg-Kohn第二定理表明，如果n(r)是体系正确的密度分布，则E[n(r)]是最低的能量，即体系的基态能量。这一定理为通过变分原理求解体系的基态能量提供了依据。在实际计算中，可以通过构造合适的能量泛函E[n(r)]，并对其进行变分，找到使能量最小的电子密度分布n(r)，从而得到体系的基态能量。例如，在计算分子的基态能量时，可以将分子的电子密度分布作为变量，通过优化能量泛函来确定分子的最稳定结构和基态能量。Hobenberg-Kohn定理的证明采用了反证法。假设有另一个外部势v’(r)，其基态\Psi’也会产生相同的n(r)。由于v(r)â‰

v’(r)，所以\Psiâ‰

\Psi’（除非v’(r)-v(r)=const）。\Psi与\Psi’满足不同的Schrödinger方程：H\Psi=E\Psi，H’\Psi’=E’Ψ’。利用基态能量最小原理，有E=\langle\Psi|H|\Psi\rangle\lt\langle\Psi'|H|\Psi'\rangle，E'=\langle\Psi'|H'|\Psi'\rangle\lt\langle\Psi|H'|\Psi\rangle。将带撇的与不带撇的交换可得E'\lt\langle\Psi|H'|\Psi\rangle，E\lt\langle\Psi'|H|\Psi'\rangle，这两个式子相互矛盾，表明v’(r)不可能产生同样的n(r)，所以v(r)是n(r)的唯一泛函。由于v(r)决定整个H，即系统的基态能量是n(r)的唯一泛函。同理，动能T和电子-电子相互作用能U也是n(r)的唯一泛函。可定义普适函数F[n]，使得系统的基态能量泛函可写为E[n]=F[n]+\intv(r)n(r)dr。Hobenberg-Kohn定理的重要意义在于，它打破了传统量子力学中以波函数作为基本变量的局限，将多体波函数问题转化为了泛函最小化问题，大大简化了多电子体系问题的复杂度。传统的波函数方法需要处理一个包含多个变量的积分方程，计算量巨大，而密度泛函理论以电子密度作为基本变量，电子密度仅是三个变量的函数，在概念和计算上都更加方便。这使得密度泛函理论在处理大规模多电子体系时具有明显的优势，为研究分子和凝聚态的性质提供了一种有效的手段。2.2.2Kohn-Sham方程Kohn-Sham方程是密度泛函理论的核心，它基于Hobenberg-Kohn定理，将复杂的多体问题简化为单电子问题，从而使得多电子体系的计算成为可能。在密度泛函理论中，系统的基态能量泛函可表示为：E[n]=T_s[n]+E_H[n]+E_{xc}[n]+\intv_{ext}(r)n(r)dr其中，T_s[n]是无相互作用电子体系的动能泛函，E_H[n]是Hartree能，描述电子之间的经典库仑相互作用，E_{xc}[n]是交换关联能泛函，包含了电子之间的交换能和关联能，这部分是多体相互作用的复杂部分，目前还无法精确计算，\intv_{ext}(r)n(r)dr是外部势场与电子的相互作用能。为了求解能量泛函的最小值，Kohn和Sham引入了一个假设的无相互作用的参考系统，该系统的电子在一个有效势场V_{eff}(r)中运动。有效势场V_{eff}(r)包含了外部势v_{ext}(r)、Hartree势V_H(r)和交换关联势V_{xc}(r)，即：V_{eff}(r)=v_{ext}(r)+V_H(r)+V_{xc}(r)其中，Hartree势V_H(r)可表示为：V_H(r)=\int\frac{n(r')}{|r-r'|}dr'交换关联势V_{xc}(r)定义为交换关联能泛函E_{xc}[n]对电子密度n(r)的泛函导数：V_{xc}(r)=\frac{\deltaE_{xc}[n]}{\deltan(r)}对于这个无相互作用的参考系统，单电子的Kohn-Sham方程为：\left[-\frac{\hbar^2}{2m}\nabla^2+V_{eff}(r)\right]\phi_i(r)=\epsilon_i\phi_i(r)其中，-\frac{\hbar^2}{2m}\nabla^2是电子的动能算符，\phi_i(r)是第i个单电子的波函数，\epsilon_i是相应的本征能量。体系的电子密度n(r)可通过单电子波函数\phi_i(r)表示为：n(r)=\sum_{i=1}^{N}|\phi_i(r)|^2其中N是电子总数。Kohn-Sham方程的求解过程是一个自洽迭代的过程。首先，猜测一个初始的电子密度分布n^{(0)}(r)，根据上述公式计算出有效势场V_{eff}^{(0)}(r)。然后，将V_{eff}^{(0)}(r)代入Kohn-Sham方程，求解得到一组新的单电子波函数\phi_i^{(1)}(r)。根据这组新的波函数计算出新的电子密度分布n^{(1)}(r)。接着，用n^{(1)}(r)重新计算有效势场V_{eff}^{(1)}(r)，再次代入Kohn-Sham方程求解，如此反复迭代，直到前后两次迭代得到的电子密度分布或能量的差异小于某个设定的收敛阈值，此时认为计算达到收敛，得到的结果即为体系的基态电子结构和能量。例如，在计算晶体的电子结构时，通过不断迭代求解Kohn-Sham方程，最终可以得到晶体中电子的能量本征值和波函数，进而计算出晶体的各种物理性质，如能带结构、态密度等。2.2.3局域密度近似（LDA）和广义梯度近似（GGA）在密度泛函理论的实际应用中，交换关联能泛函E_{xc}[n]的精确计算是一个关键问题，但由于其复杂性，目前无法得到精确的表达式，因此通常采用近似方法来处理。局域密度近似（LDA）和广义梯度近似（GGA）是两种常用的近似方法。局域密度近似（LDA）假设电子密度在空间中缓慢变化，对于缓变的n(r)或高电子密度情况，交换关联能泛函E_{xc}[n]可以近似表示为：E_{xc}^{LDA}[n]=\intn(r)\epsilon_{xc}(n(r))dr其中，\epsilon_{xc}(n(r))是均匀电子气的交换关联能密度，它是电子密度n(r)的函数。\epsilon_{xc}(n(r))可以从均匀自由电子气的理论结果得到，对于不同的r，有不同的n(r)，相应的有不同的\epsilon_{xc}(n(r))。一种计算\epsilon_{xc}(n(r))的近似公式为（在Hartree单位下）：\epsilon_{xc}(n)=\frac{-0.458}{r_s}-\frac{0.0311}{r_s}\lnr_s+\frac{0.0155}{r_s}\ln(1+1.709r_s^{\frac{4}{3}})其中r_s是自由电子气的电子“半径”，定义为r_s=(\frac{3}{4\pin})^{\frac{1}{3}}。LDA的优点是计算相对简单，计算量较小，在一些电子密度变化较为平缓的体系中，能够给出较为合理的结果。例如，对于金属体系，LDA能够较好地描述电子的行为，计算得到的金属的电子结构和一些物理性质与实验结果有较好的一致性。然而，LDA也存在明显的局限性。它假设电子密度在空间中是均匀的，忽略了电子密度的梯度变化对交换关联能的影响。在电子密度变化较大的区域，如分子的边界、固体表面等，LDA的计算结果会出现较大偏差。对于分子体系，LDA往往会高估分子的结合能，导致计算得到的分子结构与实际情况存在一定差异。广义梯度近似（GGA）在LDA的基础上进行了改进，考虑了电子密度的空间梯度对交换关联能的影响。GGA的交换关联能泛函一般形式为：E_{xc}^{GGA}[n]=\intn(r)\epsilon_{xc}(n(r),\nablan(r))dr其中\epsilon_{xc}(n(r),\nablan(r))不仅是电子密度n(r)的函数，还与电子密度的梯度\nablan(r)有关。GGA通过引入电子密度梯度项，能够更好地描述电子密度变化较快的区域，从而提高了计算的准确性。在计算半导体材料的能带结构时，GGA考虑了电子密度在原子附近的快速变化，计算得到的能带间隙比LDA更接近实验值。常见的GGA泛函有PW91、PBE等。PBE泛函是目前应用较为广泛的一种GGA泛函，它在计算精度和计算效率之间取得了较好的平衡。PBE泛函的交换关联能密度\epsilon_{xc}^{PBE}(n,\nablan)的表达式较为复杂，包含了多个与电子密度和密度梯度相关的项。与LDA相比，GGA在处理电子密度变化较大的体系时具有明显的优势，但GGA的计算量相对较大，计算时间较长。在计算大分子体系时，GGA的计算成本会显著增加，对计算资源的要求更高。2.3声子玻尔兹曼输运方程在研究材料的热输运性质时，声子玻尔兹曼输运方程（PBTE）是一个重要的理论工具，它能够有效地描述声子在材料中的输运过程。声子是晶体中原子集体振动的量子化激发，类似于光子是电磁场的量子化激发。在晶体中，原子的热振动可以看作是一系列简谐振动的叠加，这些简谐振动的能量量子就是声子。声子的运动和相互作用决定了材料的热导率等热输运性质。声子玻尔兹曼输运方程的一般形式为：\frac{\partialn_{\mathbf{q}\nu}}{\partialt}+\dot{\mathbf{q}}\cdot\nabla_{\mathbf{r}}n_{\mathbf{q}\nu}+\dot{\mathbf{r}}\cdot\nabla_{\mathbf{q}}n_{\mathbf{q}\nu}=\left(\frac{\partialn_{\mathbf{q}\nu}}{\partialt}\right)_{scatt}其中，n_{\mathbf{q}\nu}是波矢为\mathbf{q}、频率为\nu的声子分布函数，表示在单位相空间体积内声子的数量。\frac{\partialn_{\mathbf{q}\nu}}{\partialt}是声子分布函数随时间的变化率；\dot{\mathbf{q}}\cdot\nabla_{\mathbf{r}}n_{\mathbf{q}\nu}描述了声子在实空间中的扩散对分布函数的影响，其中\dot{\mathbf{q}}是波矢\mathbf{q}的时间变化率，\nabla_{\mathbf{r}}是实空间的梯度算符；\dot{\mathbf{r}}\cdot\nabla_{\mathbf{q}}n_{\mathbf{q}\nu}表示声子在倒空间中的漂移对分布函数的影响，\dot{\mathbf{r}}是位置矢量\mathbf{r}的时间变化率，\nabla_{\mathbf{q}}是倒空间的梯度算符；\left(\frac{\partialn_{\mathbf{q}\nu}}{\partialt}\right)_{scatt}是由于声子散射导致的分布函数的变化率，声子散射是声子与晶体中的缺陷、杂质、其他声子等相互作用，从而改变声子的运动状态的过程。方程左边第一项表示声子分布函数随时间的直接变化，在稳态情况下，该项为零，即系统达到稳定状态，声子分布不再随时间变化。第二项考虑了声子在实空间中的扩散，当声子在材料中传播时，由于温度梯度等因素，声子会从高温区域向低温区域扩散，从而改变声子的分布。第三项则描述了声子在倒空间中的漂移，这与声子的色散关系密切相关，不同波矢的声子具有不同的能量和群速度，在外部作用下，声子的波矢会发生变化，导致声子在倒空间中的漂移。方程右边的散射项是声子玻尔兹曼输运方程的关键部分，它包含了声子与各种散射源的相互作用，如声子-声子散射、声子-杂质散射、声子-缺陷散射等。这些散射过程会使声子的能量、动量发生改变，从而影响声子的分布函数。在实际计算中，通常采用弛豫时间近似来简化散射项，假设声子的散射过程可以用一个弛豫时间\tau_{\mathbf{q}\nu}来描述，即\left(\frac{\partialn_{\mathbf{q}\nu}}{\partialt}\right)_{scatt}=-\frac{n_{\mathbf{q}\nu}-n_{\mathbf{q}\nu}^0}{\tau_{\mathbf{q}\nu}}，其中n_{\mathbf{q}\nu}^0是声子的平衡分布函数，通常满足玻色-爱因斯坦分布。声子玻尔兹曼输运方程的物理意义在于，它从微观层面描述了声子的输运过程，将声子的运动、散射与材料的热输运性质联系起来。通过求解声子玻尔兹曼输运方程，可以得到声子的分布函数，进而计算出材料的热导率等热输运参数。在研究石墨烯/氮化硼异质结的热输运性质时，利用声子玻尔兹曼输运方程，结合第一性原理计算得到的声子色散关系和声子散射率等参数，能够深入分析声子在异质结界面处的散射情况，以及不同波矢、频率的声子对热导率的贡献，从而揭示异质结热输运的微观机制。2.4计算方法与软件工具在本研究中，运用了多种先进的计算方法和专业软件工具，以深入探究石墨烯、氮化硼及Ⅵ族硫化物异质结的热输运性质。这些计算方法和软件工具相互配合，为研究提供了强大的技术支持，使得我们能够从原子和电子层面精确地分析异质结的结构和性质。在电子结构计算方面，采用基于密度泛函理论（DFT）的VASP（ViennaAbinitioSimulationPackage）软件进行第一性原理计算。VASP利用平面波基组和赝势或投影增强波方法（PAW）来描述电子状态，能够精确计算材料的电子结构。在计算过程中，选择广义梯度近似（GGA）中的PBE泛函来描述电子之间的交换关联作用，以提高计算的准确性。在构建石墨烯/氮化硼异质结模型后，通过VASP软件计算得到异质结的电子密度分布、能带结构等信息，从而深入分析异质结中电子的行为和相互作用。VASP软件还能进行结构优化，通过调整原子的位置和晶格参数，使异质结的总能量达到最小，得到最稳定的结构，为后续的热输运性质计算奠定基础。MaterialsStudio软件在模型构建和可视化分析中发挥了重要作用。利用其强大的建模功能，能够快速、准确地构建石墨烯、氮化硼及Ⅵ族硫化物异质结的初始结构模型。通过直观的图形界面，可以对原子的排列方式、晶格常数等进行精确设置和调整，确保模型的准确性和合理性。在构建石墨烯/二硫化钼异质结模型时，能够清晰地展示石墨烯和二硫化钼的堆叠方式、原子间的相对位置等信息，为后续的计算和分析提供了直观的参考。MaterialsStudio软件还可以对计算结果进行可视化分析，如绘制电子云分布、电荷密度图等，帮助我们更直观地理解异质结的电子结构和物理性质。在声子计算方面，借助Phonopy软件计算声子色散关系和声子态密度。Phonopy采用有限位移法，通过给原子施加微小的位移，计算体系的能量变化，进而得到力常数矩阵，通过对角化力常数矩阵得到声子的频率和本征矢，从而获得声子色散关系。在计算氮化硼/二硫化钨异质结的声子性质时，利用Phonopy软件得到的声子色散曲线，可以清晰地看到不同频率声子的传播模式和速度，分析声子在异质结中的传播特性。结合声子态密度的计算结果，能够进一步了解不同频率声子的分布情况，为研究声子对热输运的贡献提供了重要依据。Thirdorder软件用于计算三阶力常数，进而得到声子散射率。通过计算声子之间的相互作用，确定声子散射的概率和方式，从而深入了解声子在异质结中的散射机制。在研究石墨烯/二硫化钨异质结的热输运过程中，利用Thirdorder软件计算得到的声子散射率，分析声子在界面处的散射情况，以及不同声子散射过程对热导率的影响。将计算得到的声子色散关系、声子态密度和声子散射率等参数输入到ShengBTE软件中，求解声子玻尔兹曼输运方程，从而得到异质结的热导率和界面热阻等热输运性质。ShengBTE软件能够综合考虑各种因素对声子输运的影响，精确计算热输运参数。在计算石墨烯/氮化硼异质结的热导率时，ShengBTE软件根据输入的参数，考虑了声子-声子散射、声子-界面散射等因素，得到了准确的热导率数值，为分析异质结的热输运性能提供了关键数据。三、MoS_2/WS_2异质结热传输性质研究3.1引言在二维材料异质结的研究领域中，MoS_2/WS_2异质结因其独特的结构和性质，成为了热输运研究的重要对象。MoS_2和WS_2均属于过渡金属硫化物，具有相似的晶体结构，均由过渡金属原子（Mo或W）与硫原子组成类似三明治的层状结构，层内原子通过共价键相互连接，层间则通过较弱的范德华力相互作用。这种相似性使得MoS_2与WS_2能够形成稳定的异质结，且二者在原子尺度上的匹配度较高，晶格失配度较小，为研究界面特性对热输运的影响提供了理想的模型体系。从应用前景来看，MoS_2/WS_2异质结在电子器件领域展现出了巨大的潜力。在纳米电子学中，该异质结可用于构建高性能的晶体管、逻辑电路等。其独特的能带结构和电学性质，有望实现更高的电子迁移率和开关速度，降低器件的功耗，从而推动纳米电子器件向更小尺寸、更高性能的方向发展。在光电子学领域，MoS_2/WS_2异质结可用于制备光电探测器、发光二极管等光电器件。由于不同材料之间的界面效应，该异质结能够实现对光的高效吸收和发射，提高光电器件的响应速度和灵敏度，在光通信、生物医学检测等领域具有重要的应用价值。而热输运性质作为影响器件性能和稳定性的关键因素，直接关系到这些应用的实际效果。例如，在高功率电子器件中，如果热输运不畅，会导致器件温度升高，进而影响器件的电学性能，甚至引发器件的热失效。因此，深入研究MoS_2/WS_2异质结的热输运性质，对于优化器件设计、提高器件性能和可靠性具有重要的指导意义。在热输运研究方面，MoS_2/WS_2异质结具有独特的价值。与单一的二维材料相比，异质结中的界面成为了热输运过程中的关键因素。界面处原子的排列方式、电子云分布以及原子间的相互作用等都与体相材料不同，这些差异会导致声子在界面处发生散射，从而影响热导率和界面热阻等热输运性质。研究MoS_2/WS_2异质结的热输运性质，有助于深入理解声子在不同材料界面处的散射机制，揭示界面热阻的微观起源。这不仅对于丰富和完善低维材料热学理论具有重要意义，还能够为设计和制备具有低界面热阻、高导热性能的新型异质结材料提供理论依据。此外，通过研究温度、界面结构等因素对MoS_2/WS_2异质结热输运性质的影响，可以进一步拓展对热输运规律的认识，为在不同工作条件下优化材料的热性能提供指导。3.2计算模型与方法为了深入研究MoS_2/WS_2异质结的热输运性质，构建了合理的计算模型，并采用了一系列精确的计算方法和参数设置。在模型构建方面，使用MaterialsStudio软件构建MoS_2/WS_2异质结模型。首先，分别构建单层MoS_2和单层WS_2的原胞模型。单层MoS_2由一个Mo原子和两个S原子组成类似三明治的结构，Mo原子位于中间层，两侧是S原子，原子间通过共价键相互连接。同理，单层WS_2由一个W原子和两个S原子组成类似结构。将单层MoS_2和单层WS_2沿c轴方向垂直堆叠，形成MoS_2/WS_2异质结模型。在堆叠过程中，考虑了不同的原子相对位置，以模拟不同的界面结构。设置底部的S原子和W原子分别与顶部的Mo原子和S原子完全重叠的结构，这种原子排列方式形成了一种特定的界面结构，有助于研究该界面结构对热输运性质的影响。对构建好的异质结模型进行结构优化，使其达到能量最低的稳定状态。在优化过程中，充分考虑了原子间的相互作用和范德华力的影响。范德华力虽然相对较弱，但在二维材料的层间相互作用中起着重要作用，它会影响异质结的层间距和原子的相对位置，进而对热输运性质产生影响。通过优化，得到了稳定的MoS_2/WS_2异质结结构，其晶格参数、原子坐标等信息为后续的计算提供了准确的基础。【配图1张：MoS_2/WS_2异质结结构示意图，展示原子的排列和堆叠方式】在计算过程中，采用基于密度泛函理论（DFT）的VASP软件进行第一性原理计算。在参数设置方面，选择广义梯度近似（GGA）中的PBE泛函来描述电子之间的交换关联作用。PBE泛函在处理电子密度变化较为复杂的体系时，能够较好地考虑电子之间的交换和关联效应，从而提高计算结果的准确性。平面波截断能设置为500eV，这一能量值能够保证平面波基组对电子波函数的精确描述，在该截断能下，计算结果能够较好地收敛，同时又不会过度增加计算量。k点网格采用Monkhorst-Pack方法生成，设置为15×15×1，这种k点网格设置能够在保证计算精度的前提下，合理控制计算成本。在计算过程中，对力和能量的收敛标准进行了严格设置，力的收敛标准设为10⁻⁴eV/Å，能量的收敛标准设为10⁻⁶eV，以确保计算结果的准确性和可靠性。在结构优化过程中，当原子受力小于10⁻⁴eV/Å，体系总能量变化小于10⁻⁶eV时，认为结构达到稳定状态，此时得到的原子坐标和晶格参数即为优化后的结果。利用密度泛函微扰理论（DFPT）结合Phonopy软件计算声子色散关系和声子态密度。在计算过程中，采用3×3×1均匀q点网格对电子态采样，以准确描述声子的色散特性。通过给原子施加微小的位移，计算体系的能量变化，进而得到力常数矩阵，通过对角化力常数矩阵得到声子的频率和本征矢，从而获得声子色散关系。结合声子态密度的计算结果，可以深入了解不同频率声子的分布情况，以及声子在异质结中的传播特性。采用Thirdorder软件计算三阶力常数，进而得到声子散射率。在计算过程中，充分考虑了声子之间的相互作用，通过精确计算声子之间的散射概率和方式，深入了解声子在异质结中的散射机制。声子散射是影响热输运性质的重要因素，不同的声子散射过程会导致声子的能量和动量发生变化，从而影响热导率等热输运参数。将计算得到的声子色散关系、声子态密度和声子散射率等参数输入到ShengBTE软件中，求解声子玻尔兹曼输运方程，从而得到异质结的热导率和界面热阻等热输运性质。ShengBTE软件能够综合考虑各种因素对声子输运的影响，精确计算热输运参数。在计算过程中，考虑了声子-声子散射、声子-界面散射等因素，通过迭代求解声子玻尔兹曼输运方程，得到准确的热导率和界面热阻数值，为分析MoS_2/WS_2异质结的热输运性能提供关键数据。3.3结果与讨论3.3.1晶格热导率通过求解声子玻尔兹曼输运方程，计算得到了MoS_2/WS_2异质结在不同温度下的晶格热导率，结果如图3.1所示。从图中可以明显看出，随着温度的升高，MoS_2/WS_2异质结晶格热导率呈现出逐渐下降的趋势。在低温区域（如50K），晶格热导率相对较高，这主要是因为在低温下，声子的散射机制主要以声子-边界散射为主，声子平均自由程较长，能够较为自由地传播，因此热导率较高。随着温度的升高，声子-声子散射逐渐增强，特别是Umklapp散射过程变得更加频繁。Umklapp散射会导致声子的动量发生较大变化，使得声子的平均自由程减小，从而阻碍了声子的传播，导致热导率下降。在高温区域（如300K及以上），声子-声子散射占据主导地位，热导率随温度的升高而显著降低。【配图1张：MoS_2/WS_2异质结晶格热导率随温度变化曲线】与单层MoS_2和单层WS_2的晶格热导率相比，MoS_2/WS_2异质结的热导率介于两者之间。这一现象主要归因于异质结的界面特性。界面处原子的排列方式与体相材料不同，原子间的相互作用也存在差异，这会导致声子在界面处发生散射。声子在从MoS_2层传播到WS_2层时，由于界面处原子势场的变化，声子的传播方向和能量会发生改变，部分声子会被散射回原层，从而降低了声子的传输效率，导致热导率下降。但由于MoS_2和WS_2具有相似的晶体结构和原子排列方式，界面处的晶格失配度较小，这在一定程度上减小了声子散射的强度，使得异质结的热导率没有远低于单层材料，而是处于两者之间。3.3.2晶体结构优化后的MoS_2/WS_2异质结晶体结构如图3.2所示。异质结形成了稳定的六方结构，晶格常数a=b=3.19\mathring{A}，层间距为3.34\mathring{A}。在结构中，MoS_2层和WS_2层通过范德华力相互作用结合在一起。Mo-S键长为2.423\mathring{A}，W-S键长为2.425\mathring{A}，相比于单层MoS_2和WS_2中的键长（分别为2.381\mathring{A}和2.384\mathring{A}）有所增加，这是由于层间范德华力的作用使得原子间的距离增大。【配图1张：优化后的MoS_2/WS_2异质结晶体结构示意图，清晰展示原子位置和键长】这种晶体结构对热输运有着重要的影响。范德华力虽然相对较弱，但它决定了异质结的层间距和原子的相对位置，从而影响了声子在层间的传输。由于范德华力的存在，声子在层间传播时会受到一定的阻碍，这是导致异质结热导率低于单层材料的一个重要原因。此外，界面处原子的排列方式和键长的变化，会改变声子的色散关系和声子态密度，进而影响声子的传播和散射。界面处的原子排列可能会引入一些局域化的振动模式，这些模式会与声子发生相互作用，导致声子散射增强，进一步降低热导率。3.3.3色散关系和态密度MoS_2/WS_2异质结的声子色散关系和声子态密度如图3.3所示。在声子色散关系图中，可以观察到不同频率的声子模式。在低频区域，主要存在声学声子模式，包括纵向声学（LA）声子、横向声学（TA）声子和弯曲声学（ZA）声子。这些声学声子模式的频率随着波矢的变化而呈现出一定的规律，它们在热输运中起着重要的作用，主要负责热量的传输。在高频区域，主要是光学声子模式，光学声子的频率相对较高，能量较大，但由于其与其他声子的相互作用较强，散射概率较大，对热导率的贡献相对较小。【配图1张：MoS_2/WS_2异质结声子色散关系和声子态密度图】声子态密度反映了不同频率声子的分布情况。从图中可以看出，在某些频率范围内，声子态密度存在峰值，这表明在这些频率下，声子的数量较多。这些峰值的位置和强度与晶体结构密切相关，晶体结构的周期性和原子间的相互作用决定了声子的振动模式和能量分布，从而影响声子态密度。在MoS_2/WS_2异质结中，由于界面的存在，声子态密度在界面附近会发生明显的变化，这是因为界面处原子的排列和相互作用与体相不同，导致了声子振动模式的改变。声子色散关系和态密度与热导率之间存在着密切的内在联系。声子的群速度与色散关系密切相关，群速度越大，声子携带热量的能力越强。在热导率的计算公式中，声子的群速度和声子态密度都是重要的参数。声子态密度较大的频率范围内，参与热输运的声子数量较多，对热导率的贡献也较大。而声子色散关系的变化会导致声子群速度的改变，进而影响热导率。如果声子色散曲线发生变化，使得某些频率的声子群速度减小，那么这些声子对热导率的贡献也会相应减小。3.3.4格林戴森参数格林戴森参数是描述声子-声子相互作用强度的重要参数，通过计算得到了MoS_2/WS_2异质结的格林戴森参数。在异质结中，格林戴森参数的值与声子的频率和波矢有关。在低频声学声子区域，格林戴森参数相对较小，这表明低频声学声子之间的相互作用较弱，声子能够较为自由地传播，有利于热输运。随着声子频率的增加，格林戴森参数逐渐增大，特别是在高频光学声子区域，格林戴森参数较大，说明高频光学声子之间的相互作用较强，声子容易发生散射，这会阻碍热输运。【配图1张：MoS_2/WS_2异质结格林戴森参数随声子频率变化图】格林戴森参数在描述声子-声子相互作用对热输运影响中起着关键作用。声子-声子相互作用是导致声子散射的重要原因之一，而格林戴森参数直接反映了这种相互作用的强度。当格林戴森参数较小时，声子-声子散射较弱，声子的平均自由程较长，热导率较高；反之，当格林戴森参数较大时，声子-声子散射较强，声子的平均自由程减小，热导率降低。在MoS_2/WS_2异质结中，由于界面的存在，界面处原子的振动模式与体相不同，这会导致声子-声子相互作用的变化，从而影响格林戴森参数。界面处的原子振动可能会激发一些新的声子模式，这些声子模式与原有的声子模式相互作用，使得格林戴森参数在界面附近发生变化，进而影响热输运过程。3.3.5声子群速度和声子寿命计算得到的MoS_2/WS_2异质结声子群速度和声子寿命随声子频率的变化如图3.4所示。声子群速度反映了声子携带能量的传播速度，在低频声学声子区域，声子群速度较大，这意味着低频声学声子能够快速地传播热量，对热导率有较大的贡献。随着声子频率的增加，声子群速度逐渐减小，在高频光学声子区域，声子群速度较小，这使得高频光学声子在热输运中的作用相对较弱。【配图1张：MoS_2/WS_2异质结声子群速度和声子寿命随声子频率变化图】声子寿命则表示声子在传播过程中保持其状态的平均时间。在异质结中，声子寿命与声子-声子散射、声子-界面散射等因素密切相关。在低频区域，声子寿命相对较长，这是因为低频声子的散射概率较小，能够在较长时间内保持其传播状态。随着频率的升高，声子-声子散射增强，声子寿命逐渐减小。特别是在高频区域，声子-声子散射和界面散射都较为强烈，导致声子寿命显著缩短。声子群速度和声子寿命对MoS_2/WS_2异质结热输运有着重要的影响。根据热导率的计算公式，热导率与声子群速度和声子寿命成正比。声子群速度越大，声子在单位时间内传播的距离越远，携带的热量就越多；声子寿命越长，声子在传播过程中被散射的概率越小，能够更有效地传输热量。因此，在低频区域，由于声子群速度较大且声子寿命较长，热导率较高；而在高频区域，声子群速度较小且声子寿命较短，热导率较低。界面的存在会改变声子的散射情况，从而影响声子群速度和声子寿命。界面处的原子排列和相互作用与体相不同，会导致声子在界面处发生散射，使得声子群速度减小，声子寿命缩短，进而降低热导率。3.4本章小结本章通过第一性原理计算和求解声子玻尔兹曼输运方程，对MoS_2/WS_2异质结的热输运性质进行了深入研究，取得了一系列有意义的结果。在晶格热导率方面，明确了其随温度升高而下降的变化规律，揭示了低温下声子-边界散射主导，高温下声子-声子散射主导的机制，同时指出异质结热导率介于单层MoS_2和WS_2之间，这与界面处原子排列和相互作用导致的声子散射密切相关。在晶体结构研究中，确定了优化后异质结的稳定六方结构，以及晶格常数、层间距和键长等参数，阐明了范德华力对键长和热输运的影响，以及界面原子排列对声子色散关系和声子态密度的改变。对声子色散关系和声子态密度的分析，清晰地展示了不同频率声子模式的分布和特性，以及它们与晶体结构和热导率的内在联系。格林戴森参数的计算，定量地描述了声子-声子相互作用强度对热输运的影响，为理解声子散射机制提供了关键依据。通过研究声子群速度和声子寿命，明确了它们对热导率的重要作用，以及界面散射对其产生的影响。本研究的创新点在于，采用多种先进的计算方法和软件工具，从多个角度深入分析了MoS_2/WS_2异质结的热输运性质，全面揭示了声子在异质结中的传播、散射机制以及与热输运性质的关联。然而，研究也存在一些不足之处。在计算过程中，虽然采用了较为精确的近似方法，但交换关联能泛函的近似处理仍可能导致一定的误差，对计算结果的准确性产生影响。在模型构建方面，仅考虑了特定的原子堆叠方式和界面结构，实际的异质结可能存在更多复杂的结构和缺陷，未来研究可进一步拓展模型，考虑更多实际因素，以更全面地研究异质结的热输运性质。四、Graphene/h-BN异质结的热输运研究4.1引言在二维材料异质结的研究体系中，Graphene/h-BN异质结凭借其独特的结构与性质，成为了热输运领域的重要研究对象。石墨烯（Graphene）作为一种由碳原子以sp²杂化轨道组成六角型呈蜂巢晶格的二维碳纳米材料，具备诸多优异性能。其力学性能卓越，拉伸强度高达130GPa，是钢铁的数百倍，这使其在柔性电子器件和高强度复合材料等领域展现出潜在应用价值，如在可穿戴电子设备中作为柔性基板，为其他电子元件提供稳定支撑，且能承受一定程度的弯曲和拉伸，保障设备正常运行。电学性能上，石墨烯是零带隙的半金属材料，载流子迁移率室温下可达2×10⁵cm²/(V・s)，电导率可达10⁶S/m，在高速电子器件和高频电路方面应用前景广阔，例如其制作的场效应晶体管开关速度比传统硅基晶体管更快，能显著提升集成电路的运行速度；在超级电容器中作为电极材料，可提高电容器的充放电效率和循环稳定性。在热学性能方面，石墨烯拥有超高的热导率，理论值可达5000W/(m・K)以上，这使其成为热管理领域的理想材料，在电子芯片散热等方面发挥着关键作用。六方氮化硼（h-BN）的结构与石墨烯类似，同样具有层状结构，由硼原子和氮原子交替排列组成六角形的蜂窝状晶格，层间通过较弱的范德华力相互作用结合在一起。h-BN是宽带隙半导体材料，带隙约为5.9eV，具有良好的绝缘性能，可用于制作电子器件中的绝缘层，有效隔离不同导电部分，防止电流泄漏，确保电子器件安全稳定运行，如在集成电路中作为绝缘介质，将不同晶体管隔开，提高电路集成度和性能。在热学性能方面，h-BN具有较高的热导率，室温下可达300-400W/(m・K)，同时热膨胀系数较低，在高温环境下性能稳定，不易因温度变化而变形或损坏，常被用作高温结构材料和散热材料，在电子封装、航空航天等领域发挥重要作用，例如在航空发动机高温部件中使用，可提高部件耐高温性能和热稳定性；也可用于制作电子设备散热基板，快速传导芯片产生的热量，提高设备散热效率。此外，h-BN还具有优异的化学稳定性和润滑性，在高温、强酸碱等恶劣环境下化学性质稳定，不易被腐蚀，其层状结构赋予的润滑性能使其可作为高温润滑剂使用，在高温、高压等特殊工况下，有效降低摩擦系数，减少机械部件磨损，提高机械设备运行效率和寿命。当石墨烯与六方氮化硼形成异质结时，由于二者之间的协同效应和界面相互作用，Graphene/h-BN异质结展现出许多独特的物理性质，在多个领域具有广阔的应用前景。在电子器件领域，该异质结结合了石墨烯的高导电性和六方氮化硼的良好绝缘性，可用于制作高性能的场效应晶体管。在石墨烯/六方氮化硼异质结场效应晶体管中，六方氮化硼作为绝缘衬底，能够有效减少漏电现象，提高晶体管的开关性能和稳定性；而石墨烯作为导电沟道，凭借其高载流子迁移率，可实现高速电子传输，从而提高晶体管的工作频率和运算速度，有望推动集成电路向更小尺寸、更高性能方向发展。在传感器领域，Graphene/h-BN异质结对某些气体分子具有特殊的吸附和电学响应特性，可用于制备高灵敏度的气体传感器。当特定气体分子吸附在异质结表面时，会引起石墨烯电学性能的变化，通过检测这种变化可实现对气体的高灵敏度检测，在环境监测、生物医学检测等领域具有重要应用价值。在光电器件领域，Graphene/h-BN异质结可用于制备发光二极管、光电探测器等。由于异质结界面处的能带结构和电子态密度发生变化，能够实现对光的高效吸收和发射，提高光电器件的响应速度和灵敏度，在光通信、生物医学成像等领域具有潜在应用。热输运性质作为影响Graphene/h-BN异质结在上述应用中性能和稳定性的关键因素，直接关系到这些应用的实际效果。在高功率电子器件中，若热输运不畅，会导致器件温度升高，进而影响器件电学性能，甚至引发热失效。在高频电子器件中，热量的积累会导致载流子迁移率下降，器件性能恶化。因此，深入研究Graphene/h-BN异质结的热输运性质，对于优化器件设计、提高器件性能和可靠性具有重要的指导意义。从理论层面来看，研究Graphene/h-BN异质结的热输运性质，有助于深入理解声子在不同材料界面处的散射机制，揭示界面热阻的微观起源，丰富和完善低维材料热学理论。从实际应用角度出发，通过研究温度、界面结构等因素对Graphene/h-BN异质结热输运性质的影响，可以为设计和制备具有低界面热阻、高导热性能的新型异质结材料提供理论依据，推动其在电子、能源、航空航天等领域的广泛应用。4.2模型构建与计算方法为深入探究Graphene/h-BN异质结的热输运性质，构建了精准的原子模型，并采用一系列先进的计算方法和合理的参数设置。运用MaterialsStudio软件构建Graphene/h-BN异质结的原子模型。首先，分别构建单层石墨烯和单层六方氮化硼的原胞模型。单层石墨烯由碳原子以sp²杂化轨道组成六角型呈蜂巢晶格结构，每个碳原子与相邻的三个碳原子通过共价键相连，键长约为1.42Å。单层六方氮化硼由硼原子和氮原子交替排列组成六角形的蜂窝状晶格，B-N键长约为1.45Å，由于硼原子和氮原子的电负性差异，B-N键具有一定的极性。将单层石墨烯和单层六方氮化硼沿c轴方向垂直堆叠，形成Graphene/h-BN异质结模型。在堆叠过程中，考虑了不同的原子相对位置，以模拟不同的界面结构。设置了AB堆叠方式，即石墨烯的碳原子位于氮化硼六角形晶格的中心上方，这种堆叠方式下，石墨烯与六方氮化硼之间的原子相互作用较为稳定，界面处的原子排列呈现出一定的规律性，有助于研究该界面结构对热输运性质的影响。同时，也考虑了AA堆叠方式，此时石墨烯的碳原子与氮化硼的硼原子或氮原子直接对齐，这种堆叠方式下，界面处的原子排列和相互作用与AB堆叠有所不同，通过对比两种堆叠方式下的热输运性质，能够更全面地了解界面结构对异质结热输运的影响。对构建好的异质结模型进行结构优化，使其达到能量最低的稳定状态。在优化过程中，充分考虑了原子间的相互作用和范德华力的影响。范德华力虽然相对较弱，但在二维材料的层间相互作用中起着重要作用，它会影响异质结的层间距和原子的相对位置，进而对热输运性质产生影响。通过优化，得到了稳定的Graphene/h-BN异质结结构，其晶格参数、原子坐标等信息为后续的计算提供了准确的基础。【配图1张：Graphene/h-BN异质结结构示意图，展示原子的排列和堆叠方式】在计算过程中，采用基于密度泛函理论（DFT）的VASP软件进行第一性原理计算。在参数设置方面，选择广义梯度近似（GGA）中的PBE泛函来描述电子之间的交换关联作用。PBE泛函在处理电子密度变化较为复杂的体系时，能够较好地考虑电子之间的交换和关联效应，从而提高计算结果的准确性。平面波截断能设置为500eV，这一能量值能够保证平面波基组对电子波函数的精确描述，在该截断能下，计算结果能够较好地收敛，同时又不会过度增加计算量。k点网格采用Monkhorst-Pack方法生成，设置为15×15×1，这种k点网格设置能够在保证计算精度的前提下，合理控制计算成本。在计算过程中，对力和能量的收敛标准进行了严格设置，力的收敛标准设为10⁻⁴eV/Å，能量的收敛标准设为10⁻⁶eV，以确保计算结果的准确性和可靠性。在结构优化过程中，当原子受力小于10⁻⁴eV/Å，体系总能量变化小于10⁻⁶eV时，认为结构达到稳定状态，此时得到的原子坐标和晶格参数即为优化后的结果。利用密度泛函微扰理论（DFPT）结合Phonopy软件计算声子色散关系和声子态密度。在计算过程中，采用3×3×1均匀q点网格对电子态采样，以准确描述声子的色散特性。通过给原子施加微小的位移，计算体系的能量变化，进而得到力常数矩阵，通