基于第一性原理剖析硼氮共掺杂超小直径碳纳米管的结构与性能

基于第一性原理剖析硼氮共掺杂超小直径碳纳米管的结构与性能一、引言1.1研究背景与意义碳纳米管（CarbonNanotubes，CNTs）自1991年被日本电气公司的饭岛澄男发现以来，凭借其独特的结构和优异的性能，在材料科学、纳米技术等众多领域引发了广泛关注与深入研究。碳纳米管可视为由石墨烯片层卷曲而成的无缝、中空管体，依据石墨烯片层数，分为单壁碳纳米管（Single-WalledCarbonNanotubes，SWCNTs）和多壁碳纳米管（Multi-WalledCarbonNanotubes，MWCNTs）。单壁碳纳米管由一层石墨烯片构成，直径通常在0.75-3nm，长度在1-50μm；多壁碳纳米管含有多层石墨烯片，层数在2-50之间，层间距约0.34±0.01nm，与石墨层间距相近，典型直径和长度分别为2-30nm和0.1-50μm。碳纳米管的结构赋予其一系列卓越特性。在机械性能方面，其强度和刚度极高，单壁碳纳米管强度比钢高出100倍以上，多壁碳纳米管强度也远超其他材料，同时具备超弹性和高韧性，使其在制备高强度材料领域极具应用潜力。电学性能上，碳纳米管电导率可达铜的5000倍以上，且兼具半导体和金属特性，这使其能够制备成各类电子器件，如场效应晶体管、逻辑门和存储器等，这些器件具有体积小、速度快、功耗低等优势，在电子领域应用广泛。热学性能上，碳纳米管热传导性能出色，导热系数比铜高3000倍以上，可用于制备高效散热材料和热电材料。在光学性能方面，碳纳米管对光的吸收和发射表现出独特性质，可应用于光学传感器、显示器等领域。此外，碳纳米管还具有良好的生物学性能，生物相容性佳且具备药物传递性能，可用于生物医学材料的制备，如药物载体、癌症治疗剂和基因传递载体等。正是由于这些优异性能，碳纳米管在众多领域展现出巨大的应用价值。在电子领域，碳纳米管可用于制造高性能电子器件，推动电子设备向小型化、高速化、低功耗方向发展；在能源领域，作为锂离子电池电极材料或导电添加剂，能提高电池充放电速率和循环稳定性，在超级电容器电极材料、氢气储存材料以及太阳能电池透明导电电极等方面也具有重要应用；在复合材料领域，与其他材料复合可形成高性能复合材料，如碳纳米管增强塑料、碳纳米管增强陶瓷等，广泛应用于航空航天、汽车制造等工业领域；在生物医学领域，可用于药物传递、疾病诊断和治疗等，为生物医学的发展提供新的途径和方法。然而，纯碳纳米管在实际应用中仍存在一定局限性。例如，在电子学应用中，其电学性能的可调控性有限，难以满足不同电子器件对特定电学性能的精确需求；在能源存储领域，与某些新型电极材料相比，其能量密度提升空间受限，影响了在高能量密度电池中的应用；在催化领域，其催化活性和选择性有待进一步提高，以实现更高效的化学反应。为了克服这些局限性，拓展碳纳米管的应用范围，对其进行改性研究显得尤为必要。掺杂作为一种有效的改性手段，通过向碳纳米管中引入外来原子，能够改变其电子结构、晶体结构和表面性质，从而显著提升碳纳米管的性能。在众多掺杂元素中，硼（B）和氮（N）由于原子尺寸与碳原子相近，且具有独特的电子特性，成为备受关注的掺杂选择。硼原子具有缺电子特性，氮原子具有富电子特性，将硼、氮原子引入碳纳米管后，会与碳原子形成不同的化学键和电子云分布，进而对碳纳米管的电学、力学、热学等性能产生显著影响。硼氮共掺杂相较于单一元素掺杂，具有独特的优势。一方面，硼和氮的协同作用能够更有效地调节碳纳米管的电子结构，产生一些单一掺杂所不具备的新特性。例如，在某些情况下，硼氮共掺杂可使碳纳米管产生独特的电子态分布，从而增强其电学性能的可调控性。另一方面，硼氮共掺杂还可以改善碳纳米管的表面化学活性，提高其与其他材料的相容性和界面结合力，这在复合材料制备等领域具有重要意义。超小直径碳纳米管（直径小于1nm）由于其量子限域效应更为显著，展现出与常规碳纳米管不同的物理化学性质，在纳米电子学、单分子传感器等前沿领域具有潜在的应用价值。然而，超小直径碳纳米管的性能同样需要进一步优化以满足实际应用需求。硼氮共掺杂为超小直径碳纳米管的性能提升提供了新的途径，通过精确控制硼氮原子的掺杂位置和浓度，有望实现对超小直径碳纳米管性能的精准调控，从而推动其在相关领域的实际应用。本研究运用第一性原理方法，对硼氮共掺杂超小直径碳纳米管展开深入研究。第一性原理计算基于量子力学原理，从原子和电子层面出发，无需借助经验参数，能够准确预测材料的结构和性能，为研究硼氮共掺杂超小直径碳纳米管的微观结构和性能变化提供了有力的理论工具。通过本研究，旨在揭示硼氮共掺杂对超小直径碳纳米管结构、电子性质、力学性能等方面的影响机制，为新型碳纳米管材料的设计与开发提供理论依据，推动碳纳米管在更多领域的实际应用和技术突破。1.2国内外研究现状碳纳米管自被发现以来，其独特的结构和优异性能吸引了全球科研人员的广泛关注，在理论和实验研究方面均取得了丰硕成果。在理论研究上，众多学者运用分子动力学模拟、第一性原理计算等方法，对碳纳米管的原子结构、电子性质、力学性能等进行深入探究。例如，通过分子动力学模拟，揭示了碳纳米管在拉伸、弯曲等力学作用下的原子运动和变形机制，为其在力学应用领域提供了理论依据。第一性原理计算则从电子层面出发，精确分析了碳纳米管的电子态密度、能带结构等，为其电学性能的研究和应用开发奠定了基础。在实验研究方面，科研人员不断探索新的制备方法和工艺，以提高碳纳米管的质量和产量，并对其性能进行实验表征。化学气相沉积法（CVD）、电弧放电法、激光烧蚀法等是目前常用的制备方法。其中，CVD法能够在较低温度下生长碳纳米管，且可通过控制反应条件实现对碳纳米管生长位置、管径和手性的一定程度控制，在大规模制备和应用方面具有优势；电弧放电法可制备出高质量的碳纳米管，但产量较低，制备过程中产生的杂质较多，后续提纯工艺较为复杂；激光烧蚀法制备的碳纳米管质量高，管径分布均匀，但设备昂贵，制备成本高，产量难以满足大规模应用需求。在碳纳米管的应用研究领域，其在电子学、能源、复合材料、生物医学等多个方面展现出巨大的应用潜力。在电子学领域，碳纳米管凭借其优异的电学性能，可用于制造高性能的电子器件。例如，碳纳米管场效应晶体管（CNT-FET）具有高电子迁移率、低功耗等优点，有望成为下一代集成电路的关键元件；在逻辑电路中，基于碳纳米管构建的逻辑门能够实现高速、低功耗的逻辑运算，为实现芯片的小型化和高性能化提供了新途径；在传感器方面，碳纳米管对某些气体分子具有特殊的吸附和电学响应特性，可制备高灵敏度的气体传感器，用于检测环境中的有害气体或生物分子。在能源领域，碳纳米管作为锂离子电池电极材料或导电添加剂，能够显著提高电池的充放电速率和循环稳定性；在超级电容器中，碳纳米管的高比表面积和良好导电性使其成为理想的电极材料，可实现快速充放电和高能量密度存储；此外，碳纳米管在太阳能电池透明导电电极、氢气储存材料等方面也具有重要的研究价值和应用前景。在复合材料领域，碳纳米管与其他材料复合后，能够显著增强复合材料的力学性能、电学性能和热学性能。如碳纳米管增强塑料、碳纳米管增强陶瓷等复合材料，在航空航天、汽车制造等工业领域得到广泛应用，可减轻部件重量，提高其强度和耐久性。在生物医学领域，碳纳米管的生物相容性和药物传递性能使其可用于药物载体、癌症治疗剂和基因传递载体等方面，为疾病的诊断和治疗提供了新的策略和方法。尽管碳纳米管在各个领域取得了显著进展，但在实际应用中仍面临一些挑战。例如，在电子学应用中，碳纳米管的电学性能受制备方法和生长条件影响较大，难以精确控制其电学性质，导致器件性能的一致性和稳定性较差；在能源存储领域，虽然碳纳米管能够提高电池和超级电容器的部分性能，但与理论预期仍存在差距，能量密度和功率密度的进一步提升仍是研究的重点和难点；在复合材料中，碳纳米管与基体材料的界面结合强度有待提高，以充分发挥其增强作用；在生物医学应用中，碳纳米管的生物安全性和长期稳定性还需要深入研究。为了克服这些挑战，拓展碳纳米管的应用范围，对其进行改性研究成为当前的研究热点之一。掺杂作为一种有效的改性手段，通过向碳纳米管中引入外来原子，能够改变其电子结构和晶体结构，从而改善其性能。硼氮共掺杂碳纳米管作为一种新型的掺杂碳纳米管材料，近年来受到了广泛关注。在实验研究方面，众多科研团队通过不同的制备方法成功合成了硼氮共掺杂碳纳米管。例如，北京大学邹如强课题组开发了简单溶液-热解法制备硼氮共掺杂碳纳米管材料（BCN），并研究了形貌控制机制，该材料展现出优越的电催化氧还原性能；他们还通过精确形貌控制，成功制备了系列硼氮共掺杂碳纳米管包覆中空Ni、Co、Mn等金属氧化物材料并用于高性能锂离子电池。这些研究表明，硼氮共掺杂能够显著改善碳纳米管在特定领域的性能，为其实际应用提供了新的可能。在理论研究方面，基于第一性原理的计算方法被广泛应用于研究硼氮共掺杂碳纳米管的结构和性能。韦建卫等人基于第一性原理的密度泛函理论、结合非平衡格林函数方法，计算了硼氮共掺杂情况下单壁碳纳米管的电子结构和输运特性，发现硼氮原子更趋向于形成沿管轴方向的硼氮原子对，且硼氮原子对共掺杂显著提升了半导体性单壁管（10，0）的输运特性，而对于金属型（5，5）管的掺杂使得其在小偏压区间内表现出明显的半导体特性。这些理论研究成果为理解硼氮共掺杂碳纳米管的性能变化提供了微观层面的解释，为实验研究和材料设计提供了重要的理论指导。然而，目前对于硼氮共掺杂碳纳米管的研究仍存在一些不足之处。一方面，实验研究中，制备高质量、高纯度且具有特定结构和性能的硼氮共掺杂碳纳米管仍面临挑战，制备工艺的复杂性和成本较高限制了其大规模应用。另一方面，理论研究主要集中在对常见管径和手性的碳纳米管进行硼氮共掺杂的研究，对于超小直径碳纳米管的硼氮共掺杂研究相对较少。超小直径碳纳米管由于其独特的量子限域效应，展现出与常规碳纳米管不同的物理化学性质，在纳米电子学、单分子传感器等前沿领域具有潜在的应用价值。但目前关于硼氮共掺杂对超小直径碳纳米管结构和性能影响的系统研究还较为缺乏，其掺杂机制和性能调控规律尚未完全明确。本研究将聚焦于硼氮共掺杂超小直径碳纳米管，运用第一性原理方法，深入研究硼氮共掺杂对超小直径碳纳米管结构、电子性质、力学性能等方面的影响。通过系统地改变硼氮原子的掺杂位置、浓度和比例，全面分析其对超小直径碳纳米管性能的调控规律，揭示硼氮共掺杂的作用机制。本研究旨在填补该领域在超小直径碳纳米管硼氮共掺杂研究方面的不足，为新型碳纳米管材料的设计与开发提供理论依据，推动碳纳米管在更多前沿领域的实际应用。1.3研究内容与方法本研究聚焦于硼氮共掺杂超小直径碳纳米管，旨在深入探究其结构、性能及掺杂机制，为新型碳纳米管材料的设计与应用提供坚实的理论依据。具体研究内容涵盖以下几个关键方面：结构特性研究：构建多种不同硼氮原子掺杂位置、浓度和比例的超小直径碳纳米管模型。通过第一性原理计算，精确优化这些模型的几何结构，获取其稳定的原子坐标和键长、键角等结构参数。深入分析硼氮共掺杂对超小直径碳纳米管管径、管长、手性以及原子间相互作用的影响，揭示掺杂导致的结构畸变和重构规律，为理解其性能变化奠定结构基础。电子性质研究：运用第一性原理方法，计算硼氮共掺杂超小直径碳纳米管的电子结构，包括能带结构、电子态密度、电荷密度分布等。从电子层面分析硼氮掺杂对碳纳米管电学性能的影响机制，例如探讨掺杂如何改变碳纳米管的金属-半导体特性、载流子浓度和迁移率等。研究硼氮原子与碳原子之间的电荷转移和轨道杂化情况，阐明其对电子结构和电学性能的调控作用。力学性能研究：基于第一性原理，对硼氮共掺杂超小直径碳纳米管进行力学性能模拟。计算其在拉伸、弯曲、扭转等不同力学载荷下的应力-应变关系、杨氏模量、剪切模量等力学参数。分析硼氮共掺杂对碳纳米管力学性能的影响，揭示掺杂增强或削弱力学性能的微观机制，如掺杂对碳-碳键强度、原子间结合力以及缺陷形成的影响。输运性质研究：结合非平衡格林函数方法与第一性原理，计算硼氮共掺杂超小直径碳纳米管的电子输运特性，包括电子透射系数、电流-电压曲线等。研究不同掺杂情况对碳纳米管电子输运能力的影响，分析其在纳米电子器件应用中的潜力和优势。探讨掺杂引入的杂质能级和缺陷对电子散射的影响，从而深入理解输运性质的变化规律。本研究采用基于量子力学原理的第一性原理计算方法，具体借助MaterialsStudio软件中的CASTEP模块来开展相关研究。第一性原理计算从原子和电子的基本相互作用出发，无需依赖任何经验参数，能够准确地预测材料的结构和性能，为深入探究硼氮共掺杂超小直径碳纳米管的微观特性提供了有力的工具。在计算过程中，采用平面波赝势方法（PWPM）来描述电子与原子核之间的相互作用，通过将电子波函数展开为平面波基组，结合赝势近似，有效地降低了计算量，同时保证了计算精度。电子交换关联能的计算采用广义梯度近似（GGA）下的Perdew-Burke-Ernzerhof（PBE）泛函，该泛函能够较好地描述材料中的电子交换和关联效应，适用于多种材料体系的研究。平面波截断能和k点网格的选取通过收敛性测试来确定，以确保计算结果的准确性和可靠性。在结构优化过程中，采用BFGS算法对原子坐标进行迭代优化，直至体系的能量和受力收敛到设定的阈值以下。通过以上计算方法和参数设置，能够准确地模拟硼氮共掺杂超小直径碳纳米管的结构和性能，为研究工作提供可靠的数据支持。二、理论基础与计算方法2.1第一性原理概述第一性原理，在计算物理与计算化学领域是极为关键的概念，广义上泛指所有基于量子力学原理的计算。其核心思想深深扎根于量子力学的基本原理，将由多个原子构成的体系视作由多个电子与原子核所组成的系统，并且在对问题进行处理时，尽可能减少对经验参数的依赖，最大限度地进行“非经验性”处理。在实际操作中，仅需依据电子质量、光速、质子与中子质量等少数基本的物理常量，便能够开展对体系能量以及电子结构等物理性质的计算。从理论根源来讲，第一性原理的计算主要是围绕着求解薛定谔方程展开。对于一个包含N个电子和M个原子核的多体系统，其哈密顿量可以表示为：\hat{H}=-\frac{\hbar^2}{2m_e}\sum_{i=1}^{N}\nabla_{i}^{2}-\sum_{i=1}^{N}\sum_{A=1}^{M}\frac{Z_Ae^2}{r_{iA}}+\frac{1}{2}\sum_{i\neqj}^{N}\frac{e^2}{r_{ij}}-\frac{\hbar^2}{2}\sum_{A=1}^{M}\frac{1Z_A}{m_A}\nabla_{A}^{2}+\frac{1}{2}\sum_{A\neqB}^{M}\frac{Z_AZ_Be^2}{R_{AB}}其中，第一项表示电子的动能，第二项是电子与原子核之间的库仑吸引能，第三项为电子-电子之间的库仑排斥能，第四项代表原子核的动能，第五项则是原子核-原子核之间的库仑排斥能。\hbar为约化普朗克常数，m_e是电子质量，m_A为原子核质量，Z_A是原子核A的电荷数，r_{iA}是电子i与原子核A之间的距离，r_{ij}是电子i和电子j之间的距离，R_{AB}是原子核A和原子核B之间的距离。薛定谔方程的表达式为\hat{H}\Psi=E\Psi，其中\Psi是系统的波函数，E是系统的能量。然而，由于多体系统中电子之间存在着复杂的相互作用，直接精确求解上述完整的薛定谔方程在计算上是极具挑战性的，甚至对于中等规模的体系都几乎难以实现。为了能够在实际计算中有效处理多体系统，科学家们发展了一系列的近似方法和理论框架。其中，密度泛函理论（DensityFunctionalTheory，DFT）便是在第一性原理计算中被广泛应用的重要理论。DFT的核心在于将多电子系统的总能量表达为电子密度的函数。Hohenberg-Kohn定理为DFT奠定了坚实的理论基础，该定理指出：体系的基态电子密度唯一地决定了体系的外势（例如由原子核产生的势场），进而决定了体系的所有性质。这意味着，只要确定了电子密度，体系的其他性质，如能量、电子结构等都可以通过相应的计算得出。存在一个关于电子密度的泛函，其最小值对应于系统的基态能量。基于Hohenberg-Kohn定理，Kohn和Sham提出了Kohn-Sham方程，通过引入一个虚构的非相互作用电子系统，使得多电子体系的问题可以转化为相对简单的单电子问题进行求解。在Kohn-Sham方法中，多电子体系的总能量可以表示为：E[\rho]=T_s[\rho]+V_{ne}[\rho]+V_{ee}[\rho]+E_{xc}[\rho]其中，T_s[\rho]是虚构的非相互作用电子系统的动能，V_{ne}[\rho]是电子与原子核之间的相互作用能，V_{ee}[\rho]是电子-电子之间的经典库仑相互作用能，E_{xc}[\rho]是交换关联能。交换关联能描述了电子之间的交换作用和关联作用，它包含了多体系统中电子相互作用的复杂性，目前无法精确求解，需要采用近似方法进行处理。常见的近似方法有局域密度近似（LocalDensityApproximation，LDA）和广义梯度近似（GeneralizedGradientApproximation，GGA）等。LDA假设体系的交换关联能可以通过均匀电子气的交换关联能来近似，其表达式为：E_{xc}^{LDA}[\rho]=\intd^3r\rho(r)\epsilon_{xc}(\rho(r))其中，\epsilon_{xc}(\rho(r))是密度为\rho(r)的均匀电子气的交换关联能密度。虽然LDA在许多材料模拟计算中能够给出较为合理的结果，但是对于电子密度空间变化较大的体系，其计算误差会显著增大。GGA则在LDA的基础上进行了改进，它不仅考虑了电子密度的大小，还考虑了电子密度的梯度信息，将交换关联能表示为电子密度及其梯度的泛函，即：E_{xc}^{GGA}[\rho]=\intd^3rf(\rho(r),\nabla\rho(r))其中，f(\rho(r),\nabla\rho(r))是一个与电子密度和电子密度梯度相关的函数。与LDA相比，GGA通常能够提供更为准确的能量和结构描述，在半导体材料电子结构的计算中，GGA泛函给出的半导体带隙值通常更接近实验真实值。除了DFT之外，在第一性原理计算中还有其他一些重要的理论和方法，例如Hartree-Fock方法。Hartree-Fock方法基于单电子近似，将多电子体系的波函数表示为单电子波函数的Slater行列式。通过变分原理求解体系的能量，得到Hartree-Fock方程。虽然Hartree-Fock方法能够精确计算交换能，但是由于它没有考虑电子之间的关联能，对于一些体系的计算结果与实际情况存在偏差。在实际应用中，第一性原理计算具有诸多显著的优势。它能够在原子和电子层面上深入研究材料的微观结构和性质，揭示材料性能的内在物理机制。与实验方法相比，第一性原理计算不受实验条件的限制，可以对各种复杂的材料体系进行理论研究。例如，对于一些在极端条件下（如高温、高压、强磁场等）才能合成或表现出特殊性能的材料，通过第一性原理计算可以在理论上预测其结构和性质，为实验研究提供重要的指导。此外，第一性原理计算还可以在材料设计阶段发挥重要作用。通过对不同原子组成、结构和电子态的理论计算和模拟，可以快速筛选出具有潜在优异性能的材料体系，大大缩短新材料的研发周期，降低研发成本。在材料科学领域，第一性原理计算已经成为研究材料结构、性能和相互作用的重要工具。它可以用于预测材料的晶体结构稳定性、电子结构（如能带结构、电子态密度等）、力学性能（如弹性常数、杨氏模量等）、热学性能（如热膨胀系数、热导率等）以及光学性能（如吸收光谱、发射光谱等）。在化学领域，第一性原理计算可用于研究化学反应机理、分子的电子结构和化学反应活性等。在物理学领域，它对于研究凝聚态物理中的各种现象，如超导、磁性、半导体物理等，提供了深入的理论分析手段。随着计算机技术的飞速发展，第一性原理计算的计算能力和效率不断提高，能够处理的体系规模也越来越大。同时，新的计算方法和理论不断涌现，进一步拓展了第一性原理计算的应用范围和精度。例如，近年来发展起来的杂化泛函方法，将Hartree-Fock方法中的交换能与DFT方法中的交换能进行线性组合，能够更准确地描述材料的电子结构。然而，第一性原理计算也并非完美无缺。尽管通过各种近似方法在一定程度上降低了计算量，但对于一些复杂的多体系统，其计算量仍然非常巨大，需要消耗大量的计算资源和时间。此外，目前的近似方法对于某些特殊体系，如强关联电子体系（如过渡金属氧化物等）的描述还存在一定的局限性，计算结果与实验值可能存在较大偏差。因此，不断改进和发展第一性原理计算方法，提高计算精度和效率，仍然是当前计算材料科学和计算化学领域的重要研究方向。2.2密度泛函理论密度泛函理论（DensityFunctionalTheory，DFT）作为第一性原理计算的重要理论基础，在材料科学、量子化学等领域中发挥着举足轻重的作用。其核心思想是将多电子体系的总能量表达为电子密度的泛函，从而把复杂的多体问题转化为相对简单的单电子问题进行求解。1964年，Hohenberg和Kohn提出的两个定理为DFT奠定了坚实的理论根基。Hohenberg-Kohn第一定理指出，对于一个处在外部势场V_{ext}(r)中的多电子体系，其基态电子密度\rho(r)唯一地决定了该体系的外部势场（除了一个常数项外），进而决定了体系的所有性质。这意味着，体系的所有物理性质，如能量、电子结构、电荷分布等，都可以通过电子密度\rho(r)来确定。数学上可表示为，若有两个不同的外部势场V_{ext1}(r)和V_{ext2}(r)，它们对应的基态电子密度分别为\rho_1(r)和\rho_2(r)，如果\rho_1(r)=\rho_2(r)，那么V_{ext1}(r)和V_{ext2}(r)仅相差一个常数。Hohenberg-Kohn第二定理表明，存在一个关于电子密度\rho(r)的普适泛函E[\rho]，当电子密度取体系基态的电子密度时，该泛函达到最小值，且这个最小值就是体系的基态能量。即体系的基态能量E_0满足E_0=\min_{\rho}E[\rho]，其中\min_{\rho}表示对所有可能的电子密度\rho求最小值。基于Hohenberg-Kohn定理，Kohn和Sham提出了著名的Kohn-Sham方程，这是DFT实际应用中的关键步骤。在Kohn-Sham方法中，将多电子体系的总能量E[\rho]表示为以下几个部分：E[\rho]=T_s[\rho]+V_{ne}[\rho]+V_{ee}[\rho]+E_{xc}[\rho]其中，T_s[\rho]是虚构的非相互作用电子系统的动能，它是一个仅与电子密度\rho(r)有关的泛函。在这个虚构系统中，电子之间没有相互作用，每个电子独立地在一个有效势场中运动。V_{ne}[\rho]是电子与原子核之间的相互作用能，它可以通过电子密度\rho(r)和原子核的位置及电荷来计算。V_{ee}[\rho]是电子-电子之间的经典库仑相互作用能，也称为Hartree能，可表示为V_{ee}[\rho]=\frac{1}{2}\iint\frac{\rho(r)\rho(r')}{|r-r'|}drdr'，它描述了电子之间的静电排斥作用。E_{xc}[\rho]是交换关联能，它包含了电子之间的交换作用和关联作用。交换作用源于电子的费米子特性，使得相同自旋的电子之间存在相互回避的效应；关联作用则描述了电子之间由于库仑相互作用而产生的动态相关性。交换关联能是DFT中最为复杂且难以精确求解的部分，目前主要通过近似方法来处理。在实际计算中，确定交换关联能E_{xc}[\rho]的近似形式至关重要，不同的近似方法会对计算结果的准确性和计算效率产生显著影响。以下介绍几种常用的交换关联泛函：局域密度近似（LocalDensityApproximation，LDA）：LDA是最早提出的一种交换关联能近似方法，其基本假设是体系的交换关联能可以通过均匀电子气的交换关联能来近似。在LDA中，将非均匀电子体系分割成许多微小的区域，在每个小区域内近似认为电子密度是均匀的，然后用均匀电子气的交换关联能密度\epsilon_{xc}(\rho)来计算该区域的交换关联能。整个体系的交换关联能E_{xc}^{LDA}[\rho]可表示为E_{xc}^{LDA}[\rho]=\int\rho(r)\epsilon_{xc}(\rho(r))dr，其中\rho(r)是位置r处的电子密度。LDA在许多材料模拟计算中能够给出较为合理的结果，尤其是对于电子密度变化较为缓慢的体系。这是因为在这些体系中，均匀电子气的假设相对合理，能够较好地描述电子之间的相互作用。然而，对于电子密度空间变化较大的体系，LDA的计算误差会明显增大。例如，在一些具有强局域化电子的体系或表面体系中，电子密度在空间的变化较为剧烈，LDA无法准确描述电子之间的交换关联作用，导致计算结果与实际情况存在较大偏差。广义梯度近似（GeneralizedGradientApproximation，GGA）：为了改进LDA的局限性，GGA在LDA的基础上进行了拓展，它不仅考虑了电子密度的大小，还考虑了电子密度的梯度信息。GGA认为交换关联能不仅与电子密度\rho(r)有关，还与电子密度的梯度\nabla\rho(r)相关。其交换关联能E_{xc}^{GGA}[\rho]的表达式为E_{xc}^{GGA}[\rho]=\intf(\rho(r),\nabla\rho(r))dr，其中f(\rho(r),\nabla\rho(r))是一个与电子密度和电子密度梯度相关的函数。与LDA相比，GGA能够更准确地描述电子密度变化较大的体系。在半导体材料的电子结构计算中，GGA泛函给出的半导体带隙值通常比LDA给出的值更接近实验真实值。这是因为GGA考虑了电子密度梯度信息，能够更好地反映半导体中电子的非均匀分布情况，从而更准确地描述电子之间的交换关联作用。目前常用的GGA泛函有Perdew-Wang(PW91)、Perdew-Burke-Ernzerhof(PBE)等形式。其中，PBE泛函由于其计算精度和计算效率的良好平衡，在材料科学研究中得到了广泛的应用。杂化泛函（HybridFunctional）：杂化泛函的提出是为了进一步提高交换关联能计算的精度。它将Hartree-Fock方法中的精确交换能与DFT方法中的交换能进行线性组合，从而得到计算体系的交换关联泛函。这种方法能够在一定程度上弥补DFT方法中交换能计算的不足。以HSE06杂化泛函为例，它在固体材料电子结构的计算中表现出色，能够准确地计算出固体材料的电子结构信息，为研究固体材料的光学特性以及光电子转移等问题提供了有力的支持。然而，杂化泛函的使用会显著增加计算量，这是因为Hartree-Fock方法中的精确交换能计算涉及到多中心积分，计算过程较为复杂。因此，目前杂化泛函主要适用于对计算精度要求较高的小体系计算。在实际应用中，需要根据具体的研究体系和计算需求来选择合适的交换关联泛函。对于电子密度变化较为平缓的体系，LDA可能已经能够满足计算精度要求，且计算效率较高；对于电子密度变化较大的体系，GGA通常能提供更准确的结果；而对于对交换能计算精度要求极高的小体系，杂化泛函则是更好的选择。同时，随着理论研究的不断深入和计算机技术的发展，新的交换关联泛函也在不断涌现，以满足不同体系和研究目的的需求。2.3计算软件与模型构建本研究借助MaterialsStudio软件中的CASTEP模块开展相关模拟计算工作。MaterialsStudio是一款专为材料科学领域研究者开发的可运行于PC上的模拟软件，其包含多种功能强大的模块。其中，CASTEP模块基于第一性原理，采用平面波赝势方法，能够精确地模拟材料的结构和性质，在材料的晶体结构优化、电子结构计算、力学性能分析等方面具有广泛应用。在使用CASTEP模块进行计算时，其操作流程清晰规范。首先，需要构建材料的初始模型，确定原子的种类、位置和晶胞参数等信息；然后，设置计算参数，如电子交换关联泛函、平面波截断能、k点网格等；接着，进行结构优化，通过迭代计算使体系能量和受力收敛到设定的阈值，得到稳定的原子坐标和结构参数；最后，基于优化后的结构进行各种性质的计算，如电子结构、力学性能等。该模块的优势在于其图形化界面操作方便，建模和结果分析容易，能够直观地展示材料的结构和计算结果，便于研究者理解和分析。同时，CASTEP模块的计算精度较高，能够为研究提供可靠的数据支持。在构建超小直径碳纳米管模型时，依据碳纳米管的结构特点进行构建。碳纳米管可看作是由石墨烯片层卷曲而成，其结构由手性矢量(n,m)来确定，手性矢量定义为\vec{C}_h=n\vec{a}_1+m\vec{a}_2，其中\vec{a}_1和\vec{a}_2是石墨烯的两个基矢。对于超小直径碳纳米管，选取(5,5)和(6,0)等手性类型进行研究。以(5,5)单壁碳纳米管为例，其构建步骤如下：在MaterialsStudio软件中，利用“BuildNanotube”工具，选择“Single-Walled”类型，输入手性指数(5,5)，设置合适的管长（如包含10个重复单元），即可生成初始的(5,5)单壁碳纳米管模型。此时，模型中的碳原子处于理想的晶格位置，但在实际计算前，还需对模型进行结构优化。在CASTEP模块中，设置结构优化参数，采用BFGS算法，平面波截断能设置为500eV（通过收敛性测试确定），k点网格采用Monkhorst-Pack方法生成，设置合适的k点密度（如5×5×1）。经过结构优化后，得到稳定的(5,5)单壁碳纳米管结构，其原子坐标和键长、键角等结构参数发生了微小变化，体系能量达到最低，此时的结构可用于后续的性能计算。同理，按照类似的方法构建(6,0)单壁碳纳米管模型，并进行结构优化。通过对不同手性超小直径碳纳米管模型的构建和优化，为研究硼氮共掺杂对其性能的影响提供了基础。构建硼氮共掺杂模型时，在已优化的超小直径碳纳米管模型基础上进行操作。以(5,5)单壁碳纳米管的硼氮共掺杂模型构建为例，采用替位掺杂的方式，即让硼原子（B）和氮原子（N）分别替代碳纳米管中的碳原子。为了研究不同掺杂位置和浓度对碳纳米管性能的影响，设计多种掺杂方案。例如，当掺杂浓度为5%时，选取碳纳米管上每隔20个碳原子的位置进行硼氮共掺杂，即B原子和N原子相隔20个碳原子进行替位掺杂；当掺杂浓度为10%时，每隔10个碳原子进行硼氮共掺杂。在MaterialsStudio软件中，通过“ReplaceAtom”工具，将选定位置的碳原子分别替换为硼原子和氮原子，从而构建出不同掺杂浓度和位置的硼氮共掺杂超小直径碳纳米管模型。构建完成后，同样在CASTEP模块中对这些模型进行结构优化，优化参数与未掺杂碳纳米管模型优化时保持一致。经过结构优化，硼氮共掺杂模型中的原子位置发生调整，体系达到稳定状态。此时，由于硼氮原子的引入，碳纳米管的键长、键角以及原子间相互作用发生了变化。通过分析这些结构变化，可以初步了解硼氮共掺杂对超小直径碳纳米管结构的影响。同时，这些优化后的硼氮共掺杂模型将作为后续电子性质、力学性能等研究的基础，通过对不同掺杂模型的计算和分析，深入揭示硼氮共掺杂对超小直径碳纳米管性能的影响机制。三、硼氮共掺杂超小直径碳纳米管的结构特征3.1碳纳米管的基本结构碳纳米管是一种具有独特结构的纳米材料，可看作是由石墨烯片层沿着特定方向卷曲而成的无缝、中空管体。其结构主要由碳原子以sp^2杂化方式相互连接形成六边形网格，构成管的侧壁。从碳原子层数角度，碳纳米管可分为单壁碳纳米管（SWCNTs）和多壁碳纳米管（MWCNTs）。单壁碳纳米管仅由一层石墨烯片卷曲而成，管径范围通常在0.75-3nm，长度在1-50μm，因其结构简单，在研究碳纳米管的本征性质时是重要的模型体系。多壁碳纳米管则由多层石墨烯片同轴卷曲而成，层数在2-50之间，相邻层间距约0.34±0.01nm，与石墨层间距相近，其管径一般在2-30nm，长度在0.1-50μm。多壁碳纳米管由于多层结构，在力学性能、电学性能等方面表现出与单壁碳纳米管不同的特性，且在一些应用场景中，如高强度复合材料制备，具有独特的优势。碳纳米管的结构还可依据手性矢量(n,m)进行分类，手性矢量定义为\vec{C}_h=n\vec{a}_1+m\vec{a}_2，其中\vec{a}_1和\vec{a}_2是石墨烯的两个基矢。根据n和m的取值不同，碳纳米管呈现出三种主要类型：扶手椅型、锯齿型和手性型。当n=m时，形成扶手椅型碳纳米管，其结构具有高度对称性，电学性质表现为金属性；当m=0时，得到锯齿型碳纳米管，其电学性质既可能是金属性，也可能是半导体性，取决于n的值；当n\neqm且m\neq0时，为手性型碳纳米管，具有独特的螺旋结构，电学性质通常为半导体性。不同手性的碳纳米管，由于其原子排列方式和电子云分布的差异，在电学、力学、光学等性能上存在显著不同。例如，扶手椅型碳纳米管的金属性使其在电子学领域可用于制备导电线路和电极等；半导体性的手性型碳纳米管则在纳米电子器件，如场效应晶体管的制备中具有重要应用价值。超小直径碳纳米管，通常指直径小于1nm的碳纳米管，由于其极小的管径，量子限域效应更为显著，展现出与常规碳纳米管不同的物理化学性质。在量子限域效应的作用下，超小直径碳纳米管的电子态密度分布发生变化，导致其电学性质对管径和手性的变化更为敏感。例如，一些超小直径碳纳米管可能表现出与宏观材料不同的电子输运特性，电子在其中的运动呈现出量子化特征，这为其在纳米电子学领域，如单电子晶体管、量子比特等器件的应用提供了潜在可能。在力学性能方面，超小直径碳纳米管由于管径小，结构中的缺陷相对较少，理论上具有更高的强度和弹性模量。然而，其极小的尺寸也使得在实际应用中对其操控和与其他材料的集成面临挑战。在光学性能上，超小直径碳纳米管的量子限域效应导致其光吸收和发射特性发生改变，可应用于单分子荧光传感器、纳米光源等领域。3.2硼氮共掺杂的方式与位置在硼氮共掺杂超小直径碳纳米管的研究中，硼氮原子的掺杂方式主要为替位掺杂，即硼（B）和氮（N）原子分别取代碳纳米管中的碳原子。这种掺杂方式会改变碳纳米管的原子排列和电子结构，进而影响其性能。以扶手椅型(5,5)超小直径碳纳米管为例，当进行硼氮共掺杂时，硼原子和氮原子在管上的掺杂位置存在多种可能性。它们可以相邻替位，即硼原子和氮原子直接取代相邻的两个碳原子位置，这种相邻替位会使硼氮原子间形成较强的相互作用，对碳纳米管的局部电子云分布产生显著影响。也可以间隔一定数量的碳原子进行替位，如间隔2个、3个碳原子等。不同的间隔距离会导致硼氮原子对碳纳米管电子结构的影响范围和程度有所不同。当间隔距离较大时，硼氮原子的影响相对独立，主要作用于其周围的局部区域；而间隔距离较小时，硼氮原子的影响区域可能会相互重叠，产生协同效应。为了系统研究不同掺杂位置对结构稳定性的影响，构建一系列具有不同硼氮掺杂位置的(5,5)碳纳米管模型。在模型1中，设置硼氮原子相邻替位，即硼原子和氮原子直接取代相邻的两个碳原子；在模型2中，使硼氮原子间隔2个碳原子替位；在模型3中，让硼氮原子间隔3个碳原子替位。对这些模型进行结构优化后，通过计算体系的总能量来评估结构稳定性，总能量越低，结构越稳定。计算结果表明，硼氮原子相邻替位的模型1总能量相对较高，结构稳定性较差。这是因为相邻替位时，硼氮原子的电子云相互作用较强，导致局部原子间的应力较大，从而使体系能量升高。而间隔2个碳原子替位的模型2和间隔3个碳原子替位的模型3，总能量相对较低，结构稳定性较好。其中，间隔3个碳原子替位的模型3总能量最低，结构最为稳定。这是由于适当的间隔距离能够在一定程度上缓解硼氮原子引入带来的结构应力，使体系达到更稳定的状态。进一步分析不同掺杂位置下碳纳米管的键长和键角变化，也能为结构稳定性提供依据。在硼氮原子相邻替位的模型1中，硼氮原子与周围碳原子形成的键长和键角与纯碳纳米管相比，发生了较大的畸变。硼原子与相邻碳原子的键长变长，氮原子与相邻碳原子的键长变短，同时相关键角也发生明显改变。这种较大的结构畸变增加了体系的能量，降低了结构稳定性。在间隔2个碳原子替位的模型2和间隔3个碳原子替位的模型3中，键长和键角的变化相对较小。随着间隔碳原子数的增加，键长和键角逐渐接近纯碳纳米管的数值，结构畸变程度减小，体系能量降低，结构稳定性增强。通过对不同掺杂位置下碳纳米管的电荷密度分布进行分析，也能深入理解其对结构稳定性的影响。硼氮原子的引入会导致碳纳米管电荷分布发生变化。在相邻替位的模型1中，硼氮原子周围的电荷密度分布出现明显的不均匀性，电荷在硼氮原子之间发生重新分布。这种不均匀的电荷分布会产生较大的静电相互作用，增加体系的能量，不利于结构稳定。而在间隔替位的模型中，电荷密度分布的不均匀性相对较小。随着间隔碳原子数的增加，电荷分布逐渐趋于均匀，静电相互作用减小，体系能量降低，结构稳定性提高。3.3掺杂对碳纳米管结构参数的影响硼氮共掺杂会对超小直径碳纳米管的管径、键长、键角等结构参数产生显著影响。以扶手椅型(5,5)超小直径碳纳米管为例，在进行硼氮共掺杂后，管径会发生变化。当硼氮原子以一定浓度和位置掺杂时，由于硼氮原子与碳原子的原子半径差异，会导致碳纳米管的局部结构发生畸变，进而影响管径大小。通过第一性原理计算得到，未掺杂的(5,5)碳纳米管管径约为0.68nm，而在硼氮共掺杂浓度为5%（硼氮原子间隔3个碳原子替位）时，管径增大至0.71nm。这是因为硼原子半径（0.85Å）和氮原子半径（0.75Å）均略大于碳原子半径（0.77Å），硼氮原子的引入使碳纳米管管壁受到向外的张力，从而导致管径增大。在键长方面，硼氮共掺杂同样会引起明显变化。在碳纳米管中，原本的碳-碳键长约为0.142nm。当硼原子取代碳原子时，硼与相邻碳原子形成的B-C键长会变长，例如在上述5%掺杂浓度的模型中，B-C键长达到0.150nm。这是因为硼原子的缺电子特性，使其与碳原子形成的共价键电子云分布相对较松散，导致键长增加。而氮原子取代碳原子后，N-C键长则略有缩短，约为0.140nm。这是由于氮原子的电负性大于碳原子，对电子的吸引能力更强，使得N-C键的电子云更靠近氮原子，从而键长缩短。键角的变化也是硼氮共掺杂对碳纳米管结构影响的重要方面。在未掺杂的碳纳米管中，碳原子以sp^2杂化形成六边形网格，C-C-C键角接近120°。硼氮共掺杂后，由于硼氮原子与周围碳原子的相互作用改变，键角发生明显变化。在硼氮原子相邻替位的情况下，与硼氮原子相连的碳原子的键角会显著偏离120°，例如在某个模型中，与硼氮原子相邻的C-C-C键角减小至115°。这是因为硼氮原子的引入破坏了原本碳纳米管的对称性和原子间的平衡，使得局部原子间的相互作用力发生改变，从而导致键角畸变。而在间隔替位的情况下，键角的变化相对较小，但仍与未掺杂时存在差异。例如在硼氮原子间隔3个碳原子替位的模型中，远离硼氮原子的区域键角接近120°，但在硼氮原子附近区域，键角会略微减小至118°左右。硼氮共掺杂对超小直径碳纳米管管径、键长和键角的影响，会进一步改变碳纳米管的原子间相互作用和结构稳定性。管径的变化会影响碳纳米管的空间构型和电子云分布，进而影响其电学性能；键长和键角的改变会直接影响原子间的结合力和相互作用能，对碳纳米管的力学性能和化学活性产生影响。这些结构参数的变化是理解硼氮共掺杂超小直径碳纳米管性能变化的重要基础。四、电子性质分析4.1能带结构本征超小直径碳纳米管的能带结构表现出与管径和手性密切相关的特性。以扶手椅型(5,5)超小直径碳纳米管为例，其具有金属性，能带结构中导带和价带在费米能级处相交，不存在明显的带隙。这是由于扶手椅型碳纳米管的特殊原子排列方式，使得其电子态在费米能级附近具有连续的分布，电子能够自由地在导带和价带之间跃迁，从而表现出良好的导电性。而锯齿型(6,0)超小直径碳纳米管则呈现出半导体性，存在一定的固有带隙。这是因为其原子排列方式导致电子态在费米能级附近出现了能隙，电子从价带跃迁到导带需要克服一定的能量障碍。这种能带结构的差异源于不同手性碳纳米管中碳原子的排列方式和电子云分布的不同，进而影响了电子的运动和跃迁特性。硼氮共掺杂对超小直径碳纳米管的能带结构产生了显著影响，进而改变了其带隙和电子跃迁特性。当对(5,5)碳纳米管进行硼氮共掺杂时，由于硼氮原子的引入，打破了原本碳纳米管的电子结构对称性。硼原子的缺电子特性和氮原子的富电子特性，使得电子云在硼氮原子周围发生重新分布。从能带结构上看，在费米能级附近出现了新的杂质能级。这些杂质能级的出现改变了电子的跃迁路径和能量需求。在某些掺杂情况下，原本金属性的(5,5)碳纳米管在小偏压区间内表现出明显的半导体特性。这是因为硼氮共掺杂引入的杂质能级使得电子在费米能级附近的跃迁受到限制，原本连续的导带和价带在小偏压下出现了能隙，电子需要克服一定的能量才能实现跃迁，从而表现出半导体的电学特性。对于(6,0)碳纳米管，硼氮共掺杂同样改变了其能带结构。在未掺杂时，(6,0)碳纳米管具有一定的固有带隙。硼氮共掺杂后，带隙大小发生了变化。这是因为硼氮原子与碳原子之间的电荷转移和轨道杂化，改变了价带和导带的能量位置。当硼氮原子与碳原子形成化学键时，其轨道相互作用使得价带和导带的电子态密度分布发生改变，进而导致带隙的变化。在一些掺杂模型中，带隙增大，这意味着电子从价带跃迁到导带需要更高的能量，材料的电学性能变得更加稳定；而在另一些掺杂模型中，带隙减小，电子跃迁所需能量降低，材料的电学活性增强。这种带隙的变化对超小直径碳纳米管的电学性能有着重要影响，例如在半导体器件应用中，带隙的精确调控对于器件的性能和功能起着关键作用。4.2态密度研究本征和掺杂碳纳米管的总态密度（TDOS）和分波态密度（PDOS），有助于深入理解硼氮原子对电子态分布的影响。本征(5,5)超小直径碳纳米管的总态密度在费米能级（EF）处呈现出明显的特征，由于其具有金属性，费米能级穿过导带和价带，使得在EF处态密度不为零，电子具有较高的占据概率。从分波态密度来看，碳原子的2s和2p轨道对总态密度有重要贡献。在低能量区域，2s轨道的贡献较大，形成了较为平坦的态密度分布；而在较高能量区域，2p轨道的贡献显著增强，特别是在费米能级附近，2p轨道的态密度变化较为剧烈，这与(5,5)碳纳米管的金属性密切相关，表明2p轨道上的电子在电子输运等过程中起着关键作用。硼氮共掺杂后，(5,5)碳纳米管的总态密度在费米能级附近发生了显著变化。由于硼氮原子的引入，在费米能级附近出现了新的态密度峰，这是由于硼氮原子与碳原子之间的电子相互作用导致的。从分波态密度分析，硼原子的2p轨道和氮原子的2p轨道对费米能级附近的态密度产生了重要影响。硼原子的2p轨道在费米能级以下出现了一个明显的态密度峰，这是因为硼原子的缺电子特性，使得其2p轨道与碳原子的2p轨道相互作用后，在该能量区域形成了新的电子态。氮原子的2p轨道则在费米能级以上产生了态密度峰，这是由于氮原子的富电子特性，其2p轨道与碳原子的2p轨道相互作用后，在较高能量区域引入了新的电子态。这种硼氮原子的分波态密度变化，导致了总态密度在费米能级附近的重新分布，进而改变了碳纳米管的电学性质。对于本征(6,0)超小直径碳纳米管，由于其具有半导体性，总态密度在费米能级处存在明显的能隙，态密度为零。分波态密度显示，碳原子的2s和2p轨道同样在不同能量区域对总态密度有贡献。在价带区域，2p轨道的贡献较大，形成了较高的态密度；而在导带区域，2s和2p轨道的态密度相对较低，这与(6,0)碳纳米管的半导体特性相符，表明价带中的电子主要占据2p轨道，而导带中的电子占据2s和2p轨道的概率相对较小。硼氮共掺杂(6,0)碳纳米管后，总态密度的能隙发生了变化。在一些掺杂情况下，能隙减小，这是由于硼氮原子的引入导致了新的电子态在能隙中出现。从分波态密度来看，硼原子和氮原子的2p轨道在能隙中引入了新的态密度峰。硼原子的2p轨道态密度峰靠近价带顶，氮原子的2p轨道态密度峰靠近导带底。这种在能隙中的态密度变化，使得电子从价带跃迁到导带所需的能量降低，从而改变了(6,0)碳纳米管的电学性能。硼氮原子的掺杂还导致了价带和导带中碳原子2p轨道态密度的变化，进一步影响了电子在价带和导带之间的跃迁概率。4.3电荷密度与电荷转移通过计算硼氮共掺杂超小直径碳纳米管的电荷密度分布和电荷转移情况，能够深入了解硼氮原子与碳原子间的电荷相互作用及其对电学性质的影响。以(5,5)超小直径碳纳米管为例，在未掺杂时，碳原子之间的电荷分布相对均匀，电子云呈对称分布在碳-碳键上。而当进行硼氮共掺杂后，电荷密度分布发生了显著变化。由于硼原子的电负性（2.04）小于碳原子（2.55），氮原子的电负性（3.04）大于碳原子，在硼氮原子与碳原子形成化学键的区域，电子云会发生重新分布。在硼原子与碳原子形成的B-C键中，电子云会向碳原子方向偏移，导致硼原子周围的电荷密度降低，呈现出一定的正电荷特征；而在氮原子与碳原子形成的N-C键中，电子云会向氮原子方向偏移，使氮原子周围的电荷密度升高，表现出一定的负电荷特征。通过Bader电荷分析，可以定量地研究硼氮原子与碳原子之间的电荷转移情况。计算结果表明，当硼原子取代碳原子时，硼原子向周围碳原子转移了约0.25个电子，这使得硼原子带正电荷，周围碳原子带负电荷。而氮原子取代碳原子后，氮原子从周围碳原子获得了约0.23个电子，导致氮原子带负电荷，周围碳原子带正电荷。这种电荷转移现象打破了原本碳纳米管中电荷分布的对称性，在硼氮原子附近形成了局部的电荷不均匀区域。电荷密度分布和电荷转移的变化对超小直径碳纳米管的电学性质产生了重要影响。电荷的重新分布改变了碳纳米管中电子的能量状态和运动特性。由于硼氮原子周围电荷密度的变化，使得电子在这些区域的运动受到阻碍或促进。在硼原子附近，由于电荷密度降低，电子的能量升高，电子在该区域的运动变得相对困难；而在氮原子附近，电荷密度升高，电子的能量降低，电子在该区域的运动相对容易。这种电子能量状态和运动特性的改变，会影响碳纳米管的电学性能，如电导率、载流子迁移率等。电荷转移还会导致碳纳米管中出现局部的电场，这些电场会对电子的输运产生影响，进一步改变碳纳米管的电学性质。五、力学性质研究5.1杨氏模量杨氏模量是衡量材料抵抗弹性形变能力的重要力学参数，它反映了材料在弹性限度内应力与应变的比例关系。对于碳纳米管而言，杨氏模量是评估其力学性能的关键指标之一，它不仅与碳纳米管的原子结构、化学键特性密切相关，还受到掺杂等因素的显著影响。为了深入探究硼氮共掺杂对超小直径碳纳米管杨氏模量的影响，本研究运用第一性原理方法，对本征和硼氮共掺杂的超小直径碳纳米管进行了细致的模拟计算。在计算过程中，严格遵循相关的理论和方法，确保计算结果的准确性和可靠性。对于本征超小直径碳纳米管，以扶手椅型(5,5)和锯齿型(6,0)两种常见手性结构为例进行分析。通过精确的计算，得到扶手椅型(5,5)本征超小直径碳纳米管的杨氏模量约为1.25TPa，而锯齿型(6,0)本征超小直径碳纳米管的杨氏模量约为1.18TPa。这一结果与前人的研究成果具有一定的一致性，充分验证了本研究计算方法的可靠性。不同手性的本征碳纳米管杨氏模量存在差异，主要是由于其原子排列方式和电子云分布的不同。扶手椅型碳纳米管具有高度的对称性，原子间的相互作用力较为均匀，使得其在抵抗弹性形变时表现出较高的能力，从而具有较高的杨氏模量；而锯齿型碳纳米管的原子排列方式相对复杂，电子云分布也存在一定的不均匀性，导致其杨氏模量略低于扶手椅型碳纳米管。当对超小直径碳纳米管进行硼氮共掺杂后，杨氏模量发生了明显的变化。在硼氮共掺杂浓度为5%（硼氮原子间隔3个碳原子替位）的情况下，扶手椅型(5,5)碳纳米管的杨氏模量降低至1.10TPa，锯齿型(6,0)碳纳米管的杨氏模量降低至1.05TPa。这表明硼氮共掺杂在一定程度上削弱了超小直径碳纳米管的弹性性能。硼氮共掺杂导致杨氏模量降低的原因主要有以下几个方面：原子半径差异引起的结构畸变：硼原子半径（0.85Å）和氮原子半径（0.75Å）与碳原子半径（0.77Å）存在一定差异，当硼氮原子取代碳原子后，会导致碳纳米管的局部结构发生畸变。这种结构畸变使得原子间的键长和键角发生改变，破坏了原本碳纳米管中原子间的平衡状态，从而降低了原子间的相互作用力，最终导致杨氏模量下降。电子结构变化对原子间相互作用的影响：硼氮原子的引入改变了碳纳米管的电子结构。硼原子的缺电子特性和氮原子的富电子特性，使得电子云在硼氮原子周围发生重新分布，导致原子间的电子云重叠程度发生变化，进而改变了原子间的相互作用力。这种电子结构的变化使得碳纳米管在受到外力作用时，原子间的结合力减弱，更容易发生弹性形变，表现为杨氏模量的降低。杂质原子与碳原子间的键合特性：硼氮原子与碳原子形成的B-C键和N-C键的键能与碳-碳键键能存在差异。B-C键的键长较长，键能相对较低；N-C键虽然键长略短，但由于氮原子的电负性较大，与碳原子之间的电子云分布不均匀，也在一定程度上影响了原子间的结合力。这些键合特性的改变使得碳纳米管的整体力学性能发生变化，杨氏模量降低。5.2拉伸与弯曲特性为了深入研究硼氮共掺杂对超小直径碳纳米管拉伸与弯曲特性的影响，采用基于第一性原理的方法进行模拟分析。在拉伸特性研究中，构建本征和硼氮共掺杂的超小直径碳纳米管模型，对其施加轴向拉伸载荷，通过计算应力-应变曲线来分析拉伸性能。以扶手椅型(5,5)碳纳米管为例，本征(5,5)碳纳米管在拉伸过程中，应力随着应变的增加呈现出线性增长的趋势，直至达到屈服点。其屈服强度约为100GPa，断裂应变约为12%。这表明本征(5,5)碳纳米管具有较高的拉伸强度和一定的延展性。当对(5,5)碳纳米管进行硼氮共掺杂后，拉伸性能发生了显著变化。在硼氮共掺杂浓度为5%（硼氮原子间隔3个碳原子替位）的情况下，屈服强度降低至85GPa，断裂应变减小至10%。这说明硼氮共掺杂削弱了(5,5)碳纳米管的拉伸强度和延展性。这主要是因为硼氮原子的引入导致碳纳米管结构发生畸变，原子间的键长和键角改变，使得原子间的结合力减弱。在拉伸过程中，这些结构缺陷更容易引发应力集中，导致材料更早地发生屈服和断裂。在弯曲特性研究方面，对本征和硼氮共掺杂的超小直径碳纳米管模型施加弯曲载荷，分析其弯曲刚度和弯曲变形行为。本征(6,0)碳纳米管具有较高的弯曲刚度，在较小的弯曲角度下，能够保持较好的结构完整性。随着弯曲角度的增加，碳纳米管开始发生塑性变形，当弯曲角度达到一定程度时，碳纳米管会发生断裂。硼氮共掺杂对(6,0)碳纳米管的弯曲特性同样产生了明显影响。硼氮共掺杂后，(6,0)碳纳米管的弯曲刚度降低，更容易发生弯曲变形。在相同的弯曲载荷下，硼氮共掺杂(6,0)碳纳米管的弯曲角度更大，且在较小的弯曲角度下就可能出现塑性变形和断裂。这是由于硼氮共掺杂破坏了碳纳米管的结构稳定性，降低了原子间的相互作用力，使得碳纳米管在弯曲过程中抵抗变形的能力下降。通过对拉伸和弯曲特性的研究可以看出，硼氮共掺杂在一定程度上降低了超小直径碳纳米管的力学性能。然而，在某些特定的应用场景中，这种性能变化也可能带来新的应用机遇。在一些需要材料具有一定柔韧性的纳米器件中，硼氮共掺杂导致的弯曲刚度降低和延展性变化，可能使碳纳米管更适合作为柔性连接部件或传感器的敏感元件。在复合材料中，硼氮共掺杂碳纳米管与基体材料的结合方式和界面性能也会受到其拉伸和弯曲特性变化的影响，需要进一步研究优化，以充分发挥其在复合材料中的作用。5.3力学性能与结构的关联超小直径碳纳米管的力学性能与结构参数、电子性质之间存在着紧密的内在联系，硼氮共掺杂通过改变这些因素，深刻影响着碳纳米管的力学性能。从结构参数角度来看，管径、键长和键角的变化是导致力学性能改变的重要因素。在超小直径碳纳米管中，管径的大小直接影响着其力学性能。一般来说，管径越小，碳纳米管的结构越紧凑，原子间的相互作用力越强，力学性能越好。然而，硼氮共掺杂后，由于硼氮原子半径与碳原子半径的差异，会使碳纳米管管径发生变化。如前文所述，硼氮原子的引入使碳纳米管管径增大，这会导致原子间的相互作用力减弱，从而降低碳纳米管的力学性能。键长和键角的改变同样对力学性能产生显著影响。在碳纳米管中，碳-碳键的键长和键角相对稳定，保证了其良好的力学性能。硼氮共掺杂后，B-C键和N-C键的键长与碳-碳键不同，且键角也发生畸变。B-C键长变长，键能相对较低，使得原子间的结合力减弱；N-C键虽键长略短，但由于氮原子电负性大，与碳原子间电子云分布不均匀，也影响了原子间的结合力。这些键长和键角的变化，破坏了碳纳米管原本稳定的结构，使其在受力时更容易发生变形和断裂，进而降低了力学性能。电子性质方面，能带结构、态密度以及电荷密度与电荷转移等因素与力学性能密切相关。能带结构的变化直接影响电子的运动状态和能量分布，进而影响碳纳米管的力学性能。硼氮共掺杂使超小直径碳纳米管的能带结构发生改变，出现新的杂质能级。这些杂质能级的存在改变了电子的跃迁方式和能量状态，使得电子在碳纳米管中的分布不再均匀，从而影响了原子间的相互作用力。在金属性的(5,5)碳纳米管中，硼氮共掺杂导致其在小偏压区间内表现出半导体特性，电子的运动受到限制，这会使碳纳米管的力学性能发生变化。态密度的变化反映了电子在不同能量状态下的分布情况，对力学性能也有重要影响。硼氮共掺杂后，在费米能级附近出现了新的态密度峰，这表明电子在该能量区域的分布发生了改变。电子分布的变化会影响原子间的电子云重叠程度，进而改变原子间的相互作用力。如硼原子的2p轨道在费米能级以下出现态密度峰，氮原子的2p轨道在费米能级以上出现态密度峰，这种变化导致了原子间相互作用的改变，从而影响了碳纳米管的力学性能。电荷密度与电荷转移的变化同样对力学性能产生影响。硼氮原子与碳原子间的电荷转移，使得原子周围的电荷密度发生改变，形成局部的电荷不均匀区域。这种电荷分布的不均匀会产生额外的静电相互作用，影响原子间的结合力。在硼原子附近，电荷密度降低，电子能量升高，原子间结合力减弱；在氮原子附近，电荷密度升高，电子能量降低，原子间结合力相对增强。但总体上，硼氮共掺杂导致的电荷分布变化会破坏碳纳米管原本的电荷平衡，降低其力学性能。六、输运性质探讨6.1电子输运特性为深入探究硼氮共掺杂对超小直径碳纳米管电子输运特性的影响，本研究结合非平衡格林函数方法与第一性原理，对本征和硼氮共掺杂的超小直径碳纳米管进行了全面计算与分析。在电子透射谱方面，本征扶手椅型(5,5)超小直径碳纳米管展现出独特的特性。由于其具有金属性，在费米能级附近，电子透射系数较高，表明电子能够较为自由地通过碳纳米管。在费米能级处，电子透射系数接近1，这意味着电子在该能量状态下几乎可以无障碍地传输，体现了其良好的导电性能。而对于本征锯齿型(6,0)超小直径碳纳米管，由于其半导体性，在费米能级附近存在能隙，电子透射系数在能隙范围内几乎为零。这表明在该能量区间，电子难以通过碳纳米管，只有当电子能量跨越能隙进入导带时，电子透射系数才会显著增加。硼氮共掺杂对超小直径碳纳米管的电子透射谱产生了显著影响。以(5,5)碳纳米管为例，在硼氮共掺杂后，费米能级附近的电子透射系数发生了明显变化。由于硼氮原子的引入，打破了原本碳纳米管的电子结构对称性，在费米能级附近出现了新的杂质能级。这些杂质能级导致电子在传输过程中发生散射，使得电子透射系数在某些能量区间降低。在杂质能级对应的能量位置，电子透射系数出现了明显的低谷，表明电子在这些能量状态下的传输受到了阻碍。而在其他能量区间，电子透射系数也有所改变，这是由于硼氮共掺杂引起的电子云重新分布以及原子间相互作用的变化，影响了电子的传输路径和散射概率。对于(6,0)碳纳米管，硼氮共掺杂同样改变了其电子透射谱。硼氮共掺杂后，原本在能隙中几乎为零的电子透射系数在某些能量位置出现了非零值。这是因为硼氮原子与碳原子之间的电荷转移和轨道杂化，在能隙中引入了新的电子态，使得电子在这些能量状态下有了一定的传输概率。随着硼氮掺杂浓度的增加，能隙中的电子透射系数进一步增大，表明硼氮共掺杂在一定程度上拓宽了电子的传输通道，改变了(6,0)碳纳米管的半导体特性。在电导率方面，本征(5,5)碳纳米管由于其金属性，具有较高的电导率。根据计算结果，其室温下的电导率可达10^6S/cm量级。这是由于其电子结构中导带和价带在费米能级处相交，电子能够自由地在导带和价带之间跃迁，形成良好的导电通道。而本征(6,0)碳纳米管由于存在能隙，电导率相对较低。在室温下，其电导率约为10^-2S/cm量级。电子需要克服能隙的能量障碍才能实现跃迁，这限制了电子的传输，从而导致电导率较低。硼氮共掺杂对超小直径碳纳米管的电导率也产生了显著影响。对于(5,5)碳纳米管，硼氮共掺杂后，电导率在一定程度上降低。这是因为硼氮共掺杂引入的杂质能级和电子散射，阻碍了电子的传输。随着硼氮掺杂浓度的增加，电导率进一步下降。当硼氮掺杂浓度从5%增加到10%时，电导率下降了约20%。对于(6,0)碳纳米管，硼氮共掺杂后，电导率在一定程度上增加。这是因为硼氮共掺杂在能隙中引入了新的电子态，拓宽了电子的传输通道。在某些掺杂情况下，电导率可提高约一个数量级。然而，当硼氮掺杂浓度过高时，电导率的增加趋势逐渐减缓，甚至出现下降的情况。这是因为过高的掺杂浓度会导致杂质能级的重叠和电子散射的加剧，反而不利于电子的传输。6.2热输运性质热导率是衡量材料热传导能力的关键参数，在材料的热管理应用中起着至关重要的作用。对于超小直径碳纳米管而言，热导率的研究不仅有助于深入理解其热输运机制，还为其在高性能热管理材料领域的应用提供了理论基础。本征超小直径碳纳米管的热导率展现出与管径和手性紧密相关的特性。以扶手椅型(5,5)超小直径碳纳米管为例，通过基于第一性原理的计算，结合声子输运理论，得到其在室温下的热导率可达5000W/(m・K)左右。这一数值远远高于许多传统的热传导材料，如铜的热导率约为400W/(m・K)。(5,5)碳纳米管具有较高热导率的原因主要在于其原子结构的高度有序性和碳原子之间较强的共价键作用。在这种结构中，声子能够较为自由地传播，散射概率较低，从而使得热传导效率较高。对于锯齿型(6,0)超小直径碳纳米管，其室温下的热导率约为4000W/(m・K)。虽然略低于(5,5)碳纳米管，但仍然表现出优异的热传导性能。不同手性碳纳米管热导率的差异源于其原子排列方式和电子云分布的不同，进而影响了声子的传播特性。硼氮共掺杂对超小直径碳纳米管的热导率产生了显著影响。以(5,5)碳纳米管为例，在硼氮共掺杂浓度为5%（硼氮原子间隔3个碳原子替位）时，热导率降低至3500W/(m・K)左右。这表明硼氮共掺杂在一定程度上削弱了碳纳米管的热传导能力。硼氮共掺杂导致热导率降低的原因主要有以下几个方面：声子散射增强：硼氮原子的引入，由于其原子质量和原子半径与碳原子不同，会在碳纳米管内部产生晶格畸变。这种晶格畸变会破坏声子传播的周期性势场，使得声子在传播过程中更容易发生散射。声子散射概率的增加，导致声子的平均自由程减小，从而降低了热导率。在硼氮原子附近，晶格畸变较为明显，声子在此处的散射尤为显著。电子-声子相互作用改变：硼氮共掺杂改变了碳纳米管的电子结构，使得电子-声子相互作用发生变化。硼原子的缺电子特性和氮原子的富电子特性，导致电子云在硼氮原子周围重新分布。这种电子结构的变化会影响电子与声子之间的耦合强度。在某些情况下，电子-声子相互作用的增强会导致声子的能量损耗增加，从而降低热导率。电子-声子相互作用的改变还会影响声子的色散关系，进一步影响声子的传播和热输运过程。杂质能级的影响：硼氮共掺杂引入了新的杂质能级，这些杂质能级会对声子的传播产生影响。杂质能级的存在使得电子在能级之间的跃迁变得更加复杂，从而增加了电子与声子相互作用的机会。电子与声子的相互作用会导致声子的散射和能量转移，进而降低热导率。杂质能级还可能与声子发生共振散射，进一步增强声子的散射概率，降低热导率。热膨胀系数是描述材料在温度变化时尺寸变化特性的重要参数，它反映了材料内部原子间的相互作用和热振动特性。对于超小直径碳纳米管，热膨胀系数的研究对于其在高温环境下的应用具有重要意义。本征超小直径碳纳米管的热膨胀系数表现出与常规材料不同的特性。以(5,5)超小直径碳纳米管为例，在低温范围内（0-100K），其热膨胀系数呈现出负值，这意味着随着温度的升高，碳纳米管的长度会略微缩短。这种负热膨胀特性源于碳纳米管独特的原子结构和振动模式。在低温下，碳原子的振动主要以面内振动为主，这种振动模式使得碳纳米管在受热时，原子间的距离会略微减小，从而导致长度缩短。随着温度的升高（100-300K），热膨胀系数逐渐变为正值，且数值较小。在这个温度范围内，碳原子的振动模式逐渐变得复杂，面外振动的贡献逐渐增加。面外振动使得碳纳米管在受热时，原子间的距离增大，从而表现出正的热膨胀。但由于碳纳米管中碳原子之间的共价键作用较强，限制了原子的热振动幅度，因此热膨胀系数相对较小。硼氮共掺杂对超小直径碳纳米管的热膨胀系数产生了明显的影响。在硼氮共掺杂浓度为5%（硼氮原子间隔3个碳原子替位）的情况下，(5,5)碳纳米管的热膨胀系数发生了显著变化。在低温范围内，负热膨胀系数的绝对值减小，表明硼氮共掺杂削弱了碳纳米管的负热膨胀特性。这是因为硼氮原子的引入改变了碳纳米管的原子间相互作用和振动模式。硼氮原子与碳原子之间的键合特性与碳-碳键不同，导致原子间的相互作用力发生变化。这种变化使得碳原子的振动模式发生改变，面内振动的贡献相对减小，从而削弱了负热膨胀特性。在高温范围内，热膨胀系数增大，表明硼氮共掺杂使得碳纳米管在