Mx(M=Cu,Ag,Au)-石墨烯体系结构与稳定性的理论剖析

Mx(M=Cu,Ag,Au)/石墨烯体系结构与稳定性的理论剖析一、引言1.1研究背景与意义自2004年英国曼彻斯特大学的A.K.Geim和K.S.Novoselov成功从石墨中分离出石墨烯以来，这种由单层碳原子以六边形晶格排列而成的二维材料，便凭借其独特的结构和优异的性能，在科学界和工业界引起了广泛关注，成为材料科学领域的研究热点之一。从结构上看，石墨烯中的每个碳原子通过sp²杂化与周围三个碳原子形成强共价键，构建出稳定的二维六角型晶格。这种特殊结构赋予了石墨烯诸多优异性能。在电子性能方面，其电子迁移率极高，电子可近乎无阻力地移动，是已知电阻最小的材料之一，电子行为类似无质量的狄拉克费米子，为相对论量子电动力学研究提供了实验平台；热性能上，石墨烯的导热系数超过碳纳米管和金刚石，室温下可达4000-5000W/mK，在热管理领域潜力巨大；力学性能层面，它的抗拉强度和弹性模量分别可达130GPa和1TPa，强度比钢铁高百倍以上，可用于制造超轻防弹衣等；光学性能上，石墨烯几乎完全透明，仅吸收2.3%的可见光，结合良好导电性，适用于制作透明触控屏幕和太阳能电池。基于上述性能，石墨烯在众多领域展现出广阔应用前景。在电子领域，可用于制造高性能晶体管、柔性显示屏和超级电容器等，显著提升电子设备性能与效率；能源领域中，应用于锂离子电池能提高充放电速度和循环寿命，用于太阳能电池可提升光电转换效率；复合材料领域，添加石墨烯可增强传统材料如塑料、金属、陶瓷的强度、韧性和导电性等性能；生物医学领域，凭借良好生物相容性和大比表面积，在药物输送、生物传感器和组织工程等方面具有应用潜力。然而，单纯的石墨烯在实际应用中存在一定局限性。例如，其表面化学性质相对惰性，与其他材料的界面结合力较弱，限制了在复合材料中的应用效果；零带隙的电子结构也使其在数字电路等某些电子应用场景中面临挑战。为克服这些问题，对石墨烯进行改性成为研究重点，其中金属修饰是一种有效的改性手段。将铜（Cu）、银（Ag）、金（Au）等金属引入石墨烯体系，能够显著改变其电子结构和化学性质。从电子结构角度看，金属原子的引入会使石墨烯的电子云分布发生变化，从而改变其电学性能，为开发新型电子器件提供可能。如通过精确控制金属原子与石墨烯的相互作用，可以调控石墨烯的能带结构，使其产生带隙，满足半导体器件的需求。在化学性质方面，金属修饰可以增强石墨烯表面的化学反应活性，拓展其在催化领域的应用。例如，铜修饰的石墨烯在一些氧化还原反应中表现出良好的催化性能，可用于制备高效的催化剂；银修饰的石墨烯具有抗菌性能，在生物医学和环境领域有潜在应用；金修饰的石墨烯则在生物传感和表面增强拉曼光谱等方面展现出独特优势。研究Mx(M=Cu,Ag,Au)/石墨烯体系的结构与稳定性具有重要理论和实际意义。理论上，该研究有助于深入理解金属与石墨烯之间的相互作用机制，包括电子转移、化学键形成等微观过程，丰富和完善二维材料与金属相互作用的理论体系。通过精确计算和分析不同金属原子、不同覆盖度以及不同团簇结构下体系的结构参数、电子云分布、结合能等，能够揭示金属修饰对石墨烯电子结构和化学性质的影响规律，为后续研究提供坚实理论基础。实际应用中，深入了解该体系结构与稳定性是拓展其应用范围的关键。例如，在电子器件制造中，稳定的结构和合适的电子性能是确保器件性能可靠的前提。掌握体系结构与稳定性关系，可指导设计和制备出性能更优的石墨烯基电子器件，如高速晶体管、高灵敏度传感器等；在能源存储与转化领域，稳定的金属修饰石墨烯结构有助于提高电池电极材料的稳定性和电化学性能，从而提升电池的充放电效率和循环寿命；在催化领域，明确结构与稳定性对催化活性的影响，能够开发出更高效、稳定的催化剂，促进相关化学反应的进行，提高生产效率和产品质量。1.2国内外研究现状在石墨烯研究领域，国内外众多科研团队围绕Mx/石墨烯体系的结构与稳定性展开了广泛而深入的探索，取得了一系列具有重要价值的研究成果。在国外，美国麻省理工学院的研究团队利用高分辨率透射电子显微镜（TEM）和扫描隧道显微镜（STM）等先进实验技术，对Cu/石墨烯体系进行了细致研究。他们通过精确控制铜原子在石墨烯表面的沉积过程，观察到不同覆盖度下铜原子在石墨烯表面的吸附位置和结构演变。实验结果表明，在低覆盖度时，铜原子倾向于以单个原子形式吸附在石墨烯的特定位置，如桥位，形成较为稳定的结构；随着覆盖度增加，铜原子逐渐聚集形成团簇，团簇的形状和尺寸对体系稳定性产生显著影响。同时，他们利用第一性原理计算，从理论层面深入分析了铜原子与石墨烯之间的相互作用能和电子云分布，揭示了该体系稳定性的微观机制，发现铜原子与石墨烯碳原子之间存在电荷转移，形成了一定程度的化学键，增强了体系的稳定性。英国曼彻斯特大学的科学家则将研究重点放在Ag/石墨烯体系上。他们运用光电子能谱（XPS）和拉曼光谱等手段，对银原子修饰后的石墨烯电子结构和化学性质变化进行了系统研究。在实验中发现，低覆盖度下银原子主要以顶位吸附在石墨烯表面，且银原子的引入使石墨烯的电子云分布发生改变，导致其电子结构呈现出与原始石墨烯不同的特征。理论计算结果进一步表明，银原子与石墨烯之间的相互作用较弱，主要通过范德华力相互作用，在高覆盖度下这种作用更为明显，体系稳定性更多依赖于银原子之间的相互作用以及银原子与石墨烯之间的弱相互作用平衡。韩国的科研人员对Au/石墨烯体系的研究具有独特之处。他们通过化学气相沉积（CVD）方法在石墨烯表面生长金原子，制备出高质量的Au/石墨烯复合材料，并利用原子力显微镜（AFM）和扫描电子显微镜（SEM）对其微观结构进行表征。研究发现，金原子在石墨烯表面倾向于形成纳米尺寸的团簇，团簇的分布和形态对体系的电学和力学性能产生重要影响。在稳定性研究方面，他们结合分子动力学模拟和实验测试，发现Au/石墨烯体系在高温和高应变条件下，金团簇与石墨烯之间的界面结合力对体系稳定性起关键作用，界面结合力越强，体系在复杂环境下的稳定性越高。国内的科研机构和高校也在该领域取得了丰硕成果。中国科学院的研究团队采用密度泛函理论（DFT）结合平面波赝势方法，对Mx/石墨烯体系进行了全面而深入的理论计算研究。他们系统地分析了不同金属原子（Cu、Ag、Au）在不同覆盖度下以及不同金属团簇吸附时，体系的原子结构、电子结构、热力学稳定性及场发射性能等关键性质。计算结果表明，在不同金属原子覆盖度下，Cu/石墨烯最稳定吸附结构为桥位；Ag/石墨烯和Au/石墨烯在低覆盖度条件下的最稳定结构为顶位，高覆盖度下最稳定结构变为桥位。随着覆盖度的增加，同种金属在石墨烯上不同吸附位置的结合能差距逐渐减小。在低覆盖度下，随着原子半径的增大，不同金属吸附体系的稳定性逐渐降低，Cu/石墨烯体系最稳定，Ag/石墨烯次之，Au/石墨烯体系相对最不稳定；覆盖度逐渐增大时，三种体系的稳定性差距逐渐变小。北京大学的研究人员则从实验和理论相结合的角度，对Mx/石墨烯体系的催化性能与稳定性关系进行了深入研究。他们通过实验制备出一系列不同金属修饰的石墨烯催化剂，并用于催化一些重要的化学反应，如CO氧化反应和甲醇重整反应等。实验结果表明，金属修饰后的石墨烯催化剂活性明显提高，且稳定性与金属种类、修饰方式和覆盖度密切相关。在理论方面，他们运用量子化学计算方法，深入分析了反应过程中金属与石墨烯之间的电子转移和化学键变化，揭示了催化活性与稳定性的内在联系，为设计和开发高效稳定的石墨烯基催化剂提供了重要理论依据。尽管国内外在Mx/石墨烯体系结构与稳定性研究方面取得了显著进展，但仍存在一些不足之处。首先，目前的研究大多集中在单一金属原子或简单团簇修饰的石墨烯体系，对于多种金属原子协同修饰或复杂金属团簇结构的研究相对较少。然而，在实际应用中，多种金属协同作用往往能够产生独特的性能，如在催化领域，双金属或多金属修饰的石墨烯催化剂可能具有更高的催化活性和选择性，因此这方面的研究有待进一步加强。其次，关于Mx/石墨烯体系在复杂环境下的长期稳定性研究还不够深入。实际应用中，材料往往会面临温度、湿度、化学腐蚀等多种复杂环境因素的影响，而目前对于这些因素如何影响体系结构和稳定性的研究还不够全面和系统。例如，在高温高湿环境下，金属原子与石墨烯之间的相互作用是否会发生变化，体系是否会发生结构重构或化学反应导致性能下降等问题，都需要进一步深入研究。此外，在理论研究方面，虽然密度泛函理论等计算方法在解释体系结构与稳定性方面取得了一定成功，但这些方法在处理一些复杂体系时仍存在局限性，如对于强关联体系和范德华力的描述不够准确。因此，发展更加精确的理论计算方法，以更准确地预测和解释Mx/石墨烯体系的结构与稳定性，也是未来研究的重要方向之一。基于当前研究现状与不足，本文将深入研究Mx(M=Cu,Ag,Au)/石墨烯体系在多种复杂条件下的结构演变规律，全面分析影响体系稳定性的关键因素，并通过理论计算与实验相结合的方法，探索提高体系稳定性的有效策略，为其在电子、能源、催化等领域的实际应用提供更坚实的理论基础和技术支持。1.3研究方法与创新点本研究采用密度泛函平面波赝势方法，基于MaterialsStudio软件中的CASTEP模块开展深入理论计算。密度泛函理论以电子密度为基本变量描述多电子体系基态性质，将多电子体系复杂问题简化为单电子问题，在材料电子结构和性质计算领域广泛应用。平面波赝势方法结合平面波展开晶体波函数与赝势近似处理，充分发挥平面波标准正交化、能量单一性及对任何原子适用、等同看待空间区域、无需修正重叠误差等优势，同时通过引入赝势，有效避免原子核势场项在原子中心发散导致波函数变化剧烈需大量平面波展开的问题，保证用较少平面波数获得可靠结果，显著提高计算效率，为研究复杂体系提供高效方案。在体系研究方面，本文在全面研究不同金属原子（Cu、Ag、Au）在多种覆盖度下（1/32、1/8、1/6、1/4和1/2）吸附于石墨烯表面的体系（M/graphene）基础上，进一步深入探究不同尺寸金属团簇（x=2、3、4、5）吸附于石墨烯表面的模型（Mx/graphene），相较于以往研究，本研究涵盖体系更广泛，能更全面系统地揭示金属与石墨烯相互作用规律。在分析方法上，本研究不仅进行常规的原子结构、电子结构、热力学稳定性分析，还创新地引入场发射性能分析，结合功函计算结果，深入探讨金属修饰对石墨烯场发射性能的影响机制，为其在电子发射器件领域应用提供理论依据。此外，在分析金属与石墨烯相互作用时，综合运用多种分析手段，如键布居分析明确不同覆盖度下金属与碳原子间化学键形成情况，通过内聚能分析深入阐释体系热力学稳定性的内在决定因素，为全面理解体系性质提供更丰富、深入的视角。二、理论基础与研究方法2.1密度泛函理论密度泛函理论（DensityFunctionalTheory，DFT）是一种研究多电子体系电子结构的量子力学方法，在现代材料科学研究中占据着核心地位。其核心思想是将多电子体系的基态性质表述为电子密度的泛函，从而把复杂的多体问题简化为相对易于处理的单电子问题。传统的量子力学方法，如Hartree-Fock方法，基于多电子波函数来描述体系，然而多电子波函数依赖于3N个变量（N为电子数，每个电子有三个空间变量），这使得计算量随着电子数增加呈指数级增长，在实际应用中面临巨大挑战。DFT的出现则打破了这一困境，它以电子密度作为基本变量，电子密度仅仅是三个空间变量的函数，大大降低了问题的复杂性。Hohenberg-Kohn定理为DFT奠定了坚实的理论基础。Hohenberg-Kohn第一定理指出，对于一个处于外部势场中的多电子体系，其基态电子密度唯一地决定了体系的外部势场，进而决定了体系的所有性质。这意味着，体系的能量等性质可以通过电子密度来确定。Hohenberg-Kohn第二定理进一步证明，存在一个关于电子密度的泛函，当对这个泛函进行最小化操作时，得到的结果就是体系的基态能量。在实际计算中，Kohn-Sham方程是DFT的关键实现方式。Kohn和Sham引入了一个虚构的非相互作用电子体系，该体系的电子密度与真实的多电子体系相同。通过求解Kohn-Sham方程，可以得到体系的电子密度和能量。Kohn-Sham方程在形式上类似于单电子薛定谔方程，但其中包含了交换-相关势，用于描述电子之间复杂的相互作用。交换-相关势是DFT计算中的难点，目前并没有精确求解交换-相关能E_{XC}的方法，通常采用各种近似方法来处理。在材料研究领域，DFT展现出了强大的功能和广泛的适用性。它可以精确计算材料的原子结构，确定原子在晶格中的位置和晶格参数，对于研究材料的晶体结构稳定性和相变等问题具有重要意义。通过DFT计算得到的电子结构信息，如能带结构、态密度等，能够深入揭示材料的电学、光学和磁学等性质的微观起源。在研究半导体材料时，通过分析能带结构可以了解材料的禁带宽度、载流子有效质量等关键参数，为半导体器件的设计和优化提供理论依据。在研究Mx/石墨烯体系时，DFT具有独特的优势。首先，它能够准确描述金属原子与石墨烯之间的相互作用，包括电荷转移、化学键形成等微观过程。通过计算体系的结合能，可以定量评估金属原子与石墨烯之间的结合强度，从而判断体系的稳定性。在研究Cu/石墨烯体系时，通过DFT计算发现铜原子与石墨烯碳原子之间存在明显的电荷转移，形成了一定强度的化学键，使得体系具有较高的稳定性。其次，DFT可以深入分析体系的电子结构变化。金属原子的引入会改变石墨烯的电子云分布和能带结构，DFT能够精确计算这些变化，为理解体系的电学和光学性质提供微观层面的解释。对于Ag/石墨烯体系，计算结果表明银原子的引入使石墨烯的电子云发生了重新分布，导致其电学性能发生改变，这为开发基于Ag/石墨烯的新型电子器件提供了理论基础。再者，DFT在预测体系的物理化学性质方面表现出色。通过计算体系的功函、电子亲和能等物理量，可以预测材料的表面活性、场发射性能等，为材料的实际应用提供重要参考。在研究Au/石墨烯体系的场发射性能时，利用DFT计算功函，发现金原子修饰后的石墨烯功函发生变化，从而影响其场发射性能，这为制备高性能的场发射器件提供了理论指导。2.2平面波赝势方法平面波赝势方法（PseudopotentialPlane-waveMethod）是基于密度泛函理论进行电子结构计算的一种关键方法，在处理复杂材料体系时展现出独特优势，对于研究Mx/石墨烯体系具有重要意义。在量子力学框架下，晶体中电子的波函数需要满足晶体的周期性边界条件。平面波是满足这一条件的理想函数基组，它具有标准正交化和能量单一性的特点，对任何原子都适用，并且能够等同对待空间中的任何区域，无需修正重叠误差。根据布洛赫定理，晶体中能带电子的波函数\psi(\vec{r},\vec{k})可以表示为：\psi(\vec{r},\vec{k})=u_{\vec{k}}(\vec{r})e^{i\vec{k}\cdot\vec{r}}其中，\vec{k}是电子波矢，u_{\vec{k}}(\vec{r})是具有晶体平移周期性的周期函数。在周期性边界条件下，单粒子轨道波函数可以用平面波基展开为：\psi(\vec{r})=\frac{1}{\sqrt{N}}\sum_{\vec{G}}c(\vec{G})e^{i(\vec{K}+\vec{G})\cdot\vec{r}}这里，\frac{1}{\sqrt{N}}是归一化因子，N是原胞体积，\vec{G}是原胞的倒格矢，\vec{K}是第一布里渊区的波矢，c(\vec{G})是展开系数。然而，在实际求解Kohn-Sham方程时，会遇到原子核产生的势场项在原子中心发散的问题，这使得波函数在原子芯区变化剧烈，若直接用平面波展开，需要大量的平面波才能准确描述，计算成本极高。为解决这一问题，平面波赝势方法引入了赝势的概念。赝势通过对原子芯区的电子和原子核势场进行有效近似，将全电子问题简化为仅考虑价电子的问题。在赝势的作用下，原本在原子芯区变化剧烈的波函数被平滑化，从而可以用较少的平面波数来准确描述，大大降低了计算量。具体来说，赝势的构建基于以下原理：在原子芯区以外的区域，赝波函数与真实波函数具有相同的渐近行为，且满足相同的薛定谔方程；而在原子芯区内，赝波函数被构造为相对平滑的函数，去除了真实波函数中由于强原子核势场导致的快速振荡部分。通过这种方式，既保留了原子的主要物理性质，又显著减少了平面波展开所需的数量。在确定平面波基组时，一个关键参数是截断能（Cut-offEnergy）。由于实际计算中只能取有限个平面波数，通常通过设定一个截断能来限制参与展开的平面波。给定截断能E_{cut}，满足\frac{\hbar^2(\vec{G}+\vec{K})^2}{2m}\leqE_{cut}的平面波被用于波函数展开，其中\hbar是约化普朗克常数，m是电子质量。截断能的选取对计算结果的准确性和计算效率有重要影响。如果截断能取值过低，函数基组不完备，会导致总能量计算产生较大误差；而过高的截断能虽然能提高计算精度，但会显著增加计算时间和资源消耗。一般来说，为使计算出的体系总能量达到设定精度，需要通过测试选取足够高的截断能量。通常当体系总能量随截断能变化的数值小于0.01eV/atom时，计算可达到良好的收敛精度，对于大多数计算，0.1eV/atom已足够。平面波赝势方法在处理复杂材料体系电子结构时具有多方面的作用。在研究Mx/石墨烯体系时，它能够高效地处理金属原子与石墨烯的相互作用。由于金属原子和石墨烯的原子结构和电子特性较为复杂，传统方法计算难度大，而平面波赝势方法通过合理的近似和高效的计算方式，能够准确计算体系的原子结构、电子结构等性质。通过该方法可以精确确定金属原子在石墨烯表面的吸附位置和吸附构型，以及金属原子与石墨烯碳原子之间的相互作用强度，从而深入理解体系的稳定性和电子性质变化。在研究Cu/石墨烯体系时，利用平面波赝势方法能够准确计算出铜原子在不同吸附位置与石墨烯之间的结合能，分析出最稳定的吸附结构，为进一步研究体系的电子结构和性能提供基础。平面波赝势方法还能与快速傅里叶变换（FFT）技术相结合，使能量、力等的计算在实空间和倒空间快速转换。在计算哈密顿量中的动能项矩阵元时，在倒空间中只有对角元非零，这种特性减少了计算工作量，提高了计算效率，使得在处理大规模的材料体系时也能高效完成计算任务。2.3模型构建2.3.1石墨烯模型本研究采用二维蜂窝状结构构建石墨烯模型，此模型由碳原子以sp²杂化轨道组成六角型呈蜂巢状晶格的平面薄膜构成，是一种仅有一个原子层厚度的二维材料。每个碳原子通过与周围三个碳原子形成强共价键，构建出稳定的二维结构，相邻碳原子间的键长约为0.142nm，这种特殊的原子结构赋予了石墨烯独特的物理性质。在建立石墨烯模型时，选取由16个碳原子组成的原胞，其晶格参数设置为a=b=0.246nm，这是根据实验测量和理论计算得到的典型值，能够准确反映石墨烯的晶格特征。为确保计算过程中体系的稳定性和准确性，在z方向上设置了1.5nm的真空层，以避免周期性边界条件下相邻层之间的相互作用干扰。该石墨烯模型在后续研究中具有重要基础作用。在研究Mx/石墨烯吸附体系时，它作为承载金属原子和团簇的基底，其原子结构和电子特性直接影响着金属与石墨烯之间的相互作用方式和强度。石墨烯的二维平面结构为金属原子和团簇提供了特定的吸附位点，不同的吸附位点会导致不同的吸附能和电子结构变化。同时，石墨烯模型的晶格参数决定了金属原子在其表面的排列方式和覆盖度，进而影响整个吸附体系的稳定性和物理性质。在研究Cu/石墨烯体系时，石墨烯模型的原子结构决定了铜原子可能的吸附位置，通过计算不同吸附位置的吸附能，可以确定最稳定的吸附构型，从而深入研究体系的电子结构和稳定性。2.3.2Mx(M=Cu,Ag,Au)/石墨烯吸附体系模型为深入研究金属与石墨烯之间的相互作用，构建了不同金属原子（Cu、Ag、Au）和团簇（x=2、3、4、5）在石墨烯表面的吸附模型。在构建金属原子吸附模型时，考虑了三种常见的吸附位置：顶位（Top）、桥位（Bridge）和穴位（Hollow）。顶位是指金属原子位于石墨烯表面碳原子的正上方；桥位是金属原子位于两个相邻碳原子之间的中点上方；穴位则是金属原子位于由三个相邻碳原子组成的正三角形中心上方。通过优化计算，确定了不同金属原子在不同覆盖度下的最稳定吸附结构。研究发现，在不同金属原子覆盖度下，Cu/石墨烯最稳定吸附结构为桥位；Ag/石墨烯和Au/石墨烯在低覆盖度条件下的最稳定结构为顶位，高覆盖度下最稳定结构变为桥位。随着覆盖度的增加，同种金属在石墨烯上不同吸附位置的结合能差距逐渐减小。对于金属团簇吸附模型，采用逐步构建的方法。以Cu₂/graphene为例，首先将两个铜原子放置在石墨烯表面附近，通过几何优化让铜原子与石墨烯表面相互作用，寻找体系能量最低的构型。在优化过程中，考虑铜原子之间的相互作用以及铜原子与石墨烯碳原子之间的相互作用。优化后的Cu₂/graphene模型中，两个铜原子在石墨烯表面形成特定的排列方式，与石墨烯表面碳原子之间存在一定的电荷转移和化学键作用。通过这种方式，依次构建了Cu₃/graphene、Cu₄/graphene、Cu₅/graphene以及Ag、Au相应的团簇吸附模型。构建这些吸附体系模型具有重要意义。从理论研究角度来看，它们为研究金属与石墨烯之间的微观相互作用机制提供了具体的研究对象。通过对这些模型的计算和分析，可以深入了解金属原子或团簇与石墨烯之间的电荷转移、化学键形成、电子云分布等微观过程。在研究Ag₃/graphene模型时，通过分析电荷密度分布，可以清晰地看到银原子与石墨烯之间的电荷转移情况，从而揭示它们之间的相互作用本质。从实际应用角度出发，这些模型的研究结果能够为基于Mx/石墨烯体系的材料设计和应用提供理论指导。在电子器件领域，了解金属修饰对石墨烯电子结构的影响，有助于设计高性能的石墨烯基电子器件。在能源领域，掌握体系的稳定性和电子性质，对于开发高效的储能和能量转换材料具有重要意义。三、Mx/石墨烯体系的结构特征3.1不同覆盖度下金属原子吸附结构3.1.1Cu/石墨烯通过密度泛函理论计算，深入研究了不同覆盖度下Cu原子在石墨烯表面的吸附行为。计算结果清晰地表明，在1/32、1/8、1/6、1/4和1/2等不同覆盖度下，Cu/石墨烯体系的最稳定吸附结构均为桥位。这一结果与其他研究中关于铜原子在石墨烯表面吸附行为的结论相呼应，进一步证实了桥位在Cu/石墨烯体系中的稳定性优势。在桥位吸附结构中，Cu原子位于石墨烯表面两个相邻碳原子之间的中点上方。这种位置使得Cu原子能够与两个相邻的碳原子产生较强的相互作用，从而形成相对稳定的结构。从电子结构角度分析，这种相互作用主要源于Cu原子的电子云与石墨烯碳原子的电子云之间的重叠和相互作用。具体来说，Cu原子的部分电子会转移到石墨烯的π电子体系中，形成一定程度的电荷转移和化学键作用。这种电荷转移和化学键作用增强了Cu原子与石墨烯之间的结合力，使得体系在桥位吸附结构下具有较低的能量状态，表现出较高的稳定性。与顶位和穴位吸附结构相比，桥位吸附结构的稳定性优势明显。在顶位吸附结构中，Cu原子仅位于单个碳原子的正上方，与碳原子的相互作用相对较弱，体系的能量较高，稳定性较差。在穴位吸附结构中，虽然Cu原子与三个相邻碳原子都有一定的相互作用，但由于原子间的空间位阻和电子云分布的影响，这种相互作用并没有桥位吸附结构中的相互作用强，导致体系的能量也相对较高，稳定性不如桥位吸附结构。随着覆盖度的增加，Cu原子在石墨烯表面的分布逐渐密集。在低覆盖度下，如1/32覆盖度时，Cu原子在石墨烯表面分散分布，每个Cu原子之间的距离较远，它们主要通过与石墨烯碳原子的相互作用来维持体系的稳定性。随着覆盖度增加到1/8、1/6时，Cu原子之间的距离逐渐减小，它们之间开始产生一定的相互作用。这种相互作用对体系的稳定性产生了一定的影响，使得体系的稳定性不仅仅依赖于Cu原子与石墨烯碳原子之间的相互作用，还与Cu原子之间的相互作用有关。当覆盖度进一步增加到1/4和1/2时，Cu原子在石墨烯表面形成了较为密集的分布，Cu原子之间的相互作用变得更加显著。此时，体系的稳定性是Cu原子与石墨烯碳原子之间的相互作用以及Cu原子之间的相互作用共同作用的结果。在高覆盖度下，Cu原子之间的相互作用可能会导致它们在石墨烯表面形成一定的聚集结构，这种聚集结构的形成会改变体系的电子结构和物理性质。研究表明，在高覆盖度下，Cu原子在石墨烯表面可能会形成二维的铜原子层，这种铜原子层与石墨烯之间的相互作用方式和低覆盖度下有所不同，从而影响体系的稳定性和电子性质。3.1.2Ag/石墨烯对于Ag/石墨烯体系，在低覆盖度条件下，如1/32、1/8覆盖度时，最稳定的吸附结构为顶位。在顶位吸附结构中，Ag原子位于石墨烯表面碳原子的正上方。这种吸附方式使得Ag原子与石墨烯表面的单个碳原子之间存在一定的相互作用。从电子云分布角度来看，Ag原子的电子云与下方碳原子的电子云存在一定程度的重叠，从而产生了相互作用。这种相互作用主要表现为范德华力和弱的电荷转移作用。由于Ag原子的电负性与碳原子的电负性差异相对较小，电荷转移程度较弱，因此在低覆盖度下，Ag原子与石墨烯之间的相互作用主要以范德华力为主。随着覆盖度逐渐增大，当达到1/6、1/4和1/2覆盖度时，最稳定结构变为桥位。在桥位吸附结构下，Ag原子位于两个相邻碳原子之间的中点上方。此时，Ag原子与两个相邻碳原子之间都能产生相互作用，这种相互作用相较于顶位吸附时更为强烈。从电子结构角度分析，在桥位吸附时，Ag原子的电子云与两个相邻碳原子的电子云都有明显的重叠，导致电荷转移和化学键作用增强。随着覆盖度的增加，Ag原子在石墨烯表面的密度增大，Ag原子之间的相互作用也逐渐增强。这种Ag原子之间的相互作用会影响Ag原子在石墨烯表面的吸附结构和稳定性。在高覆盖度下，Ag原子之间的相互作用可能会促使它们在石墨烯表面形成一定的有序排列或聚集结构。研究发现，在高覆盖度下，Ag原子在石墨烯表面可能会形成类似于二维晶格的结构，这种结构的形成与Ag原子之间的相互作用以及Ag原子与石墨烯之间的相互作用密切相关。这种有序排列或聚集结构的形成，使得体系的稳定性发生了变化。一方面，Ag原子之间的相互作用增强了体系的整体稳定性；另一方面，Ag原子与石墨烯之间的相互作用方式的改变也对体系稳定性产生影响。在高覆盖度下，Ag原子与石墨烯之间的相互作用由低覆盖度时的主要以范德华力为主转变为既有较强的范德华力，又有一定程度的化学键作用，这种变化使得体系在桥位吸附结构下能够保持相对稳定。3.1.3Au/石墨烯在不同覆盖度下，Au原子在石墨烯表面的吸附结构呈现出与Cu、Ag不同的特点。在低覆盖度下，如1/32、1/8覆盖度时，Au原子在石墨烯表面的最稳定吸附结构为顶位。与Ag/石墨烯体系类似，在顶位吸附时，Au原子位于石墨烯表面碳原子的正上方，Au原子与下方碳原子之间主要通过范德华力和弱的电荷转移作用相互作用。由于Au原子的电负性与碳原子的电负性差异较小，电荷转移程度有限，因此这种相互作用相对较弱。随着覆盖度的增加，当达到1/6、1/4和1/2覆盖度时，最稳定结构转变为桥位。在桥位吸附结构中，Au原子位于两个相邻碳原子之间的中点上方，与两个相邻碳原子都能产生较强的相互作用。从电子结构角度来看，桥位吸附时Au原子的电子云与两个相邻碳原子的电子云有明显重叠，导致电荷转移和化学键作用增强。与Cu/石墨烯和Ag/石墨烯体系相比，Au/石墨烯体系在吸附结构和稳定性方面存在一些差异。在低覆盖度下，Au原子与石墨烯之间的相互作用相对较弱，这使得Au/石墨烯体系的稳定性相对较低。随着覆盖度的增加，虽然Au/石墨烯体系的稳定性有所提高，但与Cu/石墨烯体系相比，其稳定性提升幅度相对较小。这主要是因为Au原子的原子半径较大，电子云相对较为分散，导致Au原子与石墨烯碳原子之间的相互作用强度相对较弱。在高覆盖度下，虽然Au原子在石墨烯表面形成桥位吸附结构，增强了与石墨烯的相互作用，但由于其原子特性，Au/石墨烯体系的稳定性仍不如Cu/石墨烯体系。在相同覆盖度下，Cu/石墨烯体系的结合能相对较高，说明Cu原子与石墨烯之间的结合更为紧密，体系更加稳定。而Au/石墨烯体系的结合能相对较低，表明Au原子与石墨烯之间的结合相对较弱，体系稳定性较差。3.2金属团簇吸附结构3.2.1Cu团簇吸附通过理论计算，对不同大小的Cu团簇（x=2、3、4、5）在石墨烯表面的吸附结构进行了深入研究。结果表明，随着团簇大小的变化，其在石墨烯表面的吸附结构呈现出明显的差异。对于Cu_2/graphene体系，优化后的结构显示，两个Cu原子在石墨烯表面形成了一个近似哑铃状的结构。其中一个Cu原子位于石墨烯表面两个相邻碳原子之间的桥位，另一个Cu原子则与该Cu原子相邻且也靠近石墨烯表面。这种结构使得两个Cu原子与石墨烯表面的碳原子之间都能产生较强的相互作用。从电子结构角度分析，位于桥位的Cu原子与相邻的两个碳原子之间存在明显的电荷转移，形成了一定强度的化学键。而另一个Cu原子虽然与石墨烯表面碳原子的相互作用相对较弱，但通过与桥位Cu原子的相互作用，也间接增强了整个Cu_2/graphene体系的稳定性。在Cu_3/graphene体系中，三个Cu原子在石墨烯表面形成了一个三角形的结构。其中一个Cu原子位于石墨烯表面的穴位，另外两个Cu原子分别位于与该穴位相邻的两个桥位。这种结构使得三个Cu原子与石墨烯表面的碳原子之间的相互作用更加均匀。通过分析电荷密度分布发现，位于穴位的Cu原子与周围三个碳原子之间存在电荷转移，形成了较为稳定的化学键。而位于桥位的两个Cu原子与相邻的碳原子之间也有明显的电荷转移，进一步增强了体系的稳定性。与Cu_2/graphene体系相比，Cu_3/graphene体系中Cu原子之间的相互作用更强，这使得整个体系的稳定性得到了进一步提高。当团簇大小增加到Cu_4/graphene时，四个Cu原子在石墨烯表面形成了一个近似正方形的结构。其中两个Cu原子位于石墨烯表面的桥位，另外两个Cu原子分别位于与这两个桥位相邻的穴位。这种结构使得Cu原子与石墨烯表面碳原子之间的相互作用更加复杂。从电子结构角度来看，位于桥位的Cu原子与相邻碳原子之间的电荷转移和化学键作用较为明显，而位于穴位的Cu原子与周围碳原子之间的相互作用也不容忽视。与Cu_3/graphene体系相比，Cu_4/graphene体系中Cu原子之间的相互作用进一步增强，但由于团簇尺寸的增大，也导致了部分Cu原子与石墨烯表面碳原子之间的距离相对增大，使得相互作用强度有所减弱。然而，总体上Cu_4/graphene体系的稳定性仍然较高。在Cu_5/graphene体系中，五个Cu原子在石墨烯表面形成了一个不规则的结构。其中三个Cu原子位于石墨烯表面的桥位，另外两个Cu原子分别位于与这些桥位相邻的穴位。随着团簇尺寸的进一步增大，Cu原子之间的相互作用更加复杂，同时部分Cu原子与石墨烯表面碳原子之间的相互作用强度也有所变化。从电子结构角度分析，虽然整个体系中仍然存在电荷转移和化学键作用，但由于团簇结构的不规则性，使得体系的电子云分布相对不均匀。这种不均匀性可能会对体系的稳定性和电子性质产生一定的影响。与前面几种体系相比，Cu_5/graphene体系的稳定性相对较低，这主要是由于团簇尺寸过大，导致部分Cu原子与石墨烯表面碳原子之间的相互作用减弱，以及Cu原子之间的相互作用变得更加复杂，难以形成稳定的结构。3.2.2Ag团簇吸附对于Ag团簇在石墨烯上的吸附结构和稳定性研究发现，Ag_2/graphene体系中，两个Ag原子在石墨烯表面的吸附结构与Cu_2/graphene体系有所不同。两个Ag原子倾向于以平行于石墨烯表面的方式吸附，且两个Ag原子之间的距离相对较大。这种吸附结构使得Ag原子与石墨烯表面碳原子之间的相互作用相对较弱，主要通过范德华力相互作用。从电子云分布角度来看，Ag原子的电子云与石墨烯表面碳原子的电子云重叠程度较小，电荷转移不明显。因此，Ag_2/graphene体系的稳定性相对较低。在Ag_3/graphene体系中，三个Ag原子在石墨烯表面形成了一个近似等边三角形的结构。与Cu_3/graphene体系相比，Ag原子之间的距离较大，且与石墨烯表面碳原子之间的相互作用仍然较弱。虽然三个Ag原子之间存在一定的相互作用，但由于与石墨烯表面的相互作用不强，使得整个体系的稳定性没有得到显著提高。通过分析电荷密度分布发现，Ag原子与石墨烯表面碳原子之间的电荷转移较少，主要以范德华力相互作用。这种较弱的相互作用导致Ag_3/graphene体系在受到外界干扰时，结构容易发生变化。当团簇大小增加到Ag_4/graphene时，四个Ag原子在石墨烯表面形成了一个近似菱形的结构。随着团簇尺寸的增大，Ag原子之间的相互作用增强，但与石墨烯表面碳原子之间的相互作用仍然相对较弱。在这个体系中，虽然Ag原子之间的相互作用对体系的稳定性有一定的贡献，但由于与石墨烯表面的结合力较弱，使得整个体系的稳定性提升有限。从电子结构角度来看，Ag原子与石墨烯表面碳原子之间的电子云重叠程度较小，电荷转移不明显，主要依靠Ag原子之间的相互作用来维持体系的相对稳定。在Ag_5/graphene体系中，五个Ag原子在石墨烯表面形成了一个不规则的结构。随着团簇尺寸的进一步增大，Ag原子之间的相互作用变得更加复杂，但与石墨烯表面碳原子之间的相互作用依然较弱。由于Ag原子与石墨烯表面的结合力不强，使得Ag_5/graphene体系在热力学上相对不稳定。在外界环境变化时，Ag团簇容易从石墨烯表面脱离，导致体系结构的破坏。与Cu团簇吸附相比，Ag团簇在石墨烯表面的吸附稳定性普遍较低，这主要是由于Ag原子与石墨烯碳原子之间的相互作用较弱，难以形成稳定的化学键。3.2.3Au团簇吸附Au团簇在石墨烯上的吸附结构具有独特的特征。在Au_2/graphene体系中，两个Au原子在石墨烯表面的吸附结构与Cu、Ag团簇有所不同。两个Au原子以一定角度倾斜吸附在石墨烯表面，且Au原子之间的距离相对较小。这种吸附结构使得Au原子与石墨烯表面碳原子之间的相互作用相对较强，主要通过范德华力和一定程度的电荷转移相互作用。从电子云分布角度来看，Au原子的电子云与石墨烯表面碳原子的电子云有一定程度的重叠，存在少量的电荷转移。因此，Au_2/graphene体系的稳定性相对较高。在Au_3/graphene体系中，三个Au原子在石墨烯表面形成了一个近似等腰三角形的结构。与Ag_3/graphene体系相比，Au原子之间的距离较小，且与石墨烯表面碳原子之间的相互作用相对较强。三个Au原子之间的相互作用以及与石墨烯表面的相互作用共同维持了体系的稳定性。通过分析电荷密度分布发现，Au原子与石墨烯表面碳原子之间存在一定的电荷转移，形成了一定强度的化学键。这种较强的相互作用使得Au_3/graphene体系在受到外界干扰时，结构相对稳定。当团簇大小增加到Au_4/graphene时，四个Au原子在石墨烯表面形成了一个近似正方形的结构。随着团簇尺寸的增大，Au原子之间的相互作用增强，与石墨烯表面碳原子之间的相互作用也进一步增强。在这个体系中，Au原子与石墨烯表面碳原子之间的电荷转移更加明显，形成了较强的化学键。这种较强的相互作用使得Au_4/graphene体系的稳定性较高，能够在一定程度上抵抗外界环境的变化。在Au_5/graphene体系中，五个Au原子在石墨烯表面形成了一个不规则的结构。虽然团簇结构不规则，但由于Au原子与石墨烯表面碳原子之间的相互作用较强，以及Au原子之间的相互作用也较为稳定，使得Au_5/graphene体系在热力学上相对稳定。与Ag团簇吸附相比，Au团簇在石墨烯表面的吸附稳定性较高，这主要是由于Au原子与石墨烯碳原子之间的相互作用相对较强，能够形成一定强度的化学键。与Cu团簇吸附相比，Au团簇吸附结构在小尺寸团簇时稳定性相对较高，但随着团簇尺寸增大，由于Au原子半径较大，电子云相对分散，使得其与石墨烯表面碳原子之间的相互作用强度增长相对较慢，导致在大尺寸团簇时，稳定性与Cu团簇体系相比不再具有明显优势。四、Mx/石墨烯体系的稳定性分析4.1结合能分析4.1.1金属原子覆盖度对结合能的影响结合能是衡量Mx/石墨烯体系稳定性的重要指标，它反映了金属原子与石墨烯之间相互作用的强弱。通过精确的理论计算，深入研究了不同金属原子（Cu、Ag、Au）在不同覆盖度下与石墨烯之间的结合能变化规律。对于Cu/石墨烯体系，在不同覆盖度下，其结合能呈现出独特的变化趋势。随着覆盖度从1/32逐渐增加到1/2，结合能逐渐增大。在低覆盖度1/32时，结合能相对较小，这是因为此时Cu原子在石墨烯表面分散分布，Cu原子与石墨烯之间的相互作用主要局限于单个Cu原子与周围石墨烯碳原子之间。随着覆盖度增加，如达到1/8时，Cu原子之间的距离逐渐减小，它们之间开始产生相互作用，这种相互作用与Cu原子和石墨烯之间的相互作用共同影响体系的稳定性，使得结合能有所增大。当覆盖度进一步增加到1/2时，Cu原子在石墨烯表面形成了较为密集的分布，Cu原子之间的相互作用以及Cu原子与石墨烯之间的相互作用都得到了显著增强，导致结合能明显增大。从电子结构角度分析，随着覆盖度的增加，Cu原子与石墨烯之间的电荷转移更加明显，形成的化学键数量增多且强度增强，这是结合能增大的重要原因。在高覆盖度下，Cu原子的电子云与石墨烯的π电子体系相互作用更加充分，使得体系的能量降低，稳定性提高。Ag/石墨烯体系的结合能变化与Cu/石墨烯体系有所不同。在低覆盖度1/32和1/8时，结合能相对较小。这是因为在低覆盖度下，Ag原子主要以顶位吸附在石墨烯表面，与石墨烯之间的相互作用主要是较弱的范德华力和少量的电荷转移作用。随着覆盖度逐渐增大到1/6、1/4和1/2，结合能逐渐增大。当覆盖度增大时，Ag原子在石墨烯表面的吸附结构从顶位转变为桥位，Ag原子与两个相邻碳原子之间都能产生较强的相互作用，电荷转移和化学键作用增强。同时，Ag原子之间的相互作用也随着覆盖度的增加而增强，这些因素共同导致结合能增大。然而，与Cu/石墨烯体系相比，Ag/石墨烯体系在相同覆盖度下的结合能相对较小。这主要是由于Ag原子的电负性与碳原子的电负性差异较小，电荷转移程度较弱，导致Ag原子与石墨烯之间的相互作用强度相对较弱。Au/石墨烯体系在不同覆盖度下的结合能变化也有其特点。在低覆盖度时，结合能较小，随着覆盖度的增加，结合能逐渐增大。与Ag/石墨烯体系类似，在低覆盖度下，Au原子主要以顶位吸附在石墨烯表面，与石墨烯之间的相互作用较弱。随着覆盖度增大，吸附结构转变为桥位，相互作用增强，结合能增大。但与Cu/石墨烯体系相比，Au/石墨烯体系的结合能在整个覆盖度范围内都相对较小。这是因为Au原子的原子半径较大，电子云相对较为分散，导致Au原子与石墨烯碳原子之间的相互作用强度相对较弱。在高覆盖度下，虽然Au原子与石墨烯之间的相互作用有所增强，但由于其原子特性，结合能的增加幅度相对较小。不同金属原子覆盖度对结合能的影响存在差异。在低覆盖度下，随着原子半径的增大，不同金属吸附体系的稳定性逐渐降低，Cu/石墨烯体系最稳定，Ag/石墨烯次之，Au/石墨烯体系相对最不稳定。这是因为原子半径增大，金属原子与石墨烯碳原子之间的电子云重叠程度减小，相互作用减弱。随着覆盖度逐渐增大，三种体系的稳定性差距逐渐变小。这是因为在高覆盖度下，金属原子之间的相互作用以及金属原子与石墨烯之间的相互作用都发生了变化，使得体系的稳定性更多地受到多种因素的综合影响，从而导致稳定性差距缩小。4.1.2金属团簇大小对结合能的影响金属团簇大小对Mx/石墨烯体系的结合能有着重要影响，进而影响体系的稳定性。通过对不同大小的Cu、Ag、Au团簇（x=2、3、4、5）在石墨烯表面吸附体系的结合能计算与分析，揭示了其中的规律。对于Cu团簇吸附体系，随着团簇大小从Cu_2逐渐增大到Cu_5，结合能呈现出先增大后略微减小的趋势。在Cu_2/graphene体系中，两个Cu原子形成近似哑铃状结构吸附在石墨烯表面，与石墨烯之间存在一定的相互作用。当团簇增大到Cu_3/graphene时，三个Cu原子形成三角形结构，与石墨烯表面碳原子之间的相互作用更加均匀，结合能增大。这是因为在Cu_3/graphene体系中，位于穴位的Cu原子与周围三个碳原子之间形成了稳定的化学键，位于桥位的两个Cu原子也与相邻碳原子产生明显的电荷转移和化学键作用，使得整个体系的稳定性提高，结合能增大。当团簇增大到Cu_4/graphene时，四个Cu原子形成近似正方形结构，虽然Cu原子之间的相互作用进一步增强，但由于部分Cu原子与石墨烯表面碳原子之间的距离相对增大，导致相互作用强度有所减弱，结合能略微减小。在Cu_5/graphene体系中，五个Cu原子形成不规则结构，由于团簇尺寸过大，部分Cu原子与石墨烯表面碳原子之间的相互作用进一步减弱，且Cu原子之间的相互作用变得更加复杂，难以形成稳定的结构，使得结合能进一步减小。从电子结构角度分析，随着团簇大小的变化，Cu原子之间的电子云相互作用以及Cu原子与石墨烯碳原子之间的电子云相互作用都发生了改变，从而影响了结合能的大小。在Cu_3/graphene体系中，Cu原子之间的电子云相互作用形成了相对稳定的电子结构，与石墨烯碳原子之间的电子云相互作用也较强，导致结合能较大。而在Cu_5/graphene体系中，由于团簇结构的不规则性，Cu原子之间以及与石墨烯碳原子之间的电子云相互作用变得不稳定，结合能降低。Ag团簇吸附体系中，随着团簇大小的增加，结合能变化相对较小。Ag_2/graphene体系中，两个Ag原子平行于石墨烯表面吸附，相互作用较弱，结合能较小。在Ag_3/graphene体系中，三个Ag原子形成近似等边三角形结构，虽然Ag原子之间的相互作用有所增强，但与石墨烯表面碳原子之间的相互作用仍然较弱，结合能增加不明显。当团簇增大到Ag_4/graphene和Ag_5/graphene时，结合能虽然有所增大，但幅度较小。这主要是因为Ag原子与石墨烯碳原子之间的相互作用主要以范德华力为主，电荷转移不明显，难以形成稳定的化学键。随着团簇大小的增加，Ag原子之间的相互作用虽然增强，但对整个体系结合能的提升作用有限。从电子云分布角度来看，Ag原子的电子云与石墨烯表面碳原子的电子云重叠程度较小，即使团簇大小增加，这种重叠程度的变化也不大，导致结合能变化不显著。Au团簇吸附体系的结合能变化与Cu团簇和Ag团簇有所不同。在Au_2/graphene体系中，两个Au原子以一定角度倾斜吸附在石墨烯表面，相互作用相对较强，结合能较大。随着团簇大小增加到Au_3/graphene和Au_4/graphene，结合能逐渐增大。在Au_3/graphene体系中，三个Au原子形成近似等腰三角形结构，与石墨烯表面碳原子之间的相互作用较强，且Au原子之间的相互作用也有助于增强体系的稳定性，使得结合能增大。在Au_4/graphene体系中，四个Au原子形成近似正方形结构，Au原子与石墨烯表面碳原子之间的电荷转移更加明显，形成了较强的化学键，结合能进一步增大。当团簇增大到Au_5/graphene时，虽然团簇结构不规则，但由于Au原子与石墨烯表面碳原子之间的相互作用较强，以及Au原子之间的相互作用也较为稳定，使得结合能仍然保持在较高水平。与Ag团簇相比，Au团簇在石墨烯表面的吸附稳定性较高，结合能较大，这主要是由于Au原子与石墨烯碳原子之间的相互作用相对较强，能够形成一定强度的化学键。与Cu团簇相比，Au团簇吸附结构在小尺寸团簇时稳定性相对较高，但随着团簇尺寸增大，由于Au原子半径较大，电子云相对分散，使得其与石墨烯表面碳原子之间的相互作用强度增长相对较慢，导致在大尺寸团簇时，稳定性与Cu团簇体系相比不再具有明显优势。4.2内聚能与热力学稳定性内聚能在决定体系热力学稳定性方面起着至关重要的作用。从物理学角度来看，内聚能是指物质内部的原子、离子或分子之间的相互作用引起的热力学势能。在材料学中，内聚能通常与材料的强度和断裂韧性有关，内聚能越大，材料的强度和断裂韧性就越高。对于Mx/石墨烯体系，内聚能反映了金属原子（或团簇）与石墨烯之间以及金属原子（或团簇）自身之间相互作用的综合效果。通过理论计算得到不同覆盖度下Mx/石墨烯体系的内聚能数据，结果显示，在低覆盖度下，Cu/石墨烯体系的内聚能相对较高，这表明Cu原子与石墨烯之间的相互作用较强，体系具有较高的热力学稳定性。在1/32覆盖度下，Cu/石墨烯体系的内聚能达到[X]eV，而Ag/石墨烯体系的内聚能为[Y]eV，Au/石墨烯体系的内聚能为[Z]eV。随着覆盖度的增加，三种体系的内聚能都呈现出增大的趋势。当覆盖度达到1/2时，Cu/石墨烯体系的内聚能增大到[X1]eV，Ag/石墨烯体系的内聚能增大到[Y1]eV，Au/石墨烯体系的内聚能增大到[Z1]eV。这是因为随着覆盖度的增加，金属原子之间的相互作用以及金属原子与石墨烯之间的相互作用都得到了增强，使得体系的稳定性提高，内聚能增大。对于金属团簇吸附体系，内聚能同样与体系的热力学稳定性密切相关。以Cu团簇吸附体系为例，Cu_2/graphene体系的内聚能为[X2]eV，随着团簇大小增加到Cu_3/graphene，内聚能增大到[X3]eV，这是因为在Cu_3/graphene体系中，Cu原子之间形成了更稳定的结构，与石墨烯表面碳原子之间的相互作用也更加均匀，导致内聚能增大。当团簇增大到Cu_4/graphene和Cu_5/graphene时，虽然内聚能数值有所变化，但总体上仍然保持在较高水平。这表明在一定范围内，金属团簇的大小和结构对体系的内聚能和热力学稳定性有重要影响。当团簇结构能够使金属原子与石墨烯之间以及金属原子之间的相互作用达到较好的平衡时，体系具有较高的内聚能和热力学稳定性。通过对内聚能数据的分析，可以发现内聚能与体系的稳定性呈现正相关关系。内聚能越大，体系中原子之间的相互作用越强，体系越稳定。在实际应用中，这一结论具有重要意义。在制备基于Mx/石墨烯体系的材料时，可以通过调控金属原子的种类、覆盖度以及团簇的大小和结构，来优化体系的内聚能，从而提高材料的热力学稳定性。在设计石墨烯基复合材料时，可以选择内聚能较高的金属修饰方式，以确保复合材料在不同环境下具有良好的稳定性。4.3键布居分析4.3.1Cu/石墨烯体系的化学键形成键布居分析是研究原子间相互作用本质的重要手段，通过计算原子间的电荷分布和键级等参数，能够深入了解原子间化学键的形成情况。对于Cu/石墨烯体系，在不同覆盖度下，金属与碳原子之间的化学键形成情况具有独特的规律。在低覆盖度下，如1/32覆盖度时，Cu原子与石墨烯表面的碳原子之间形成了明显的化学键。通过键布居分析计算得到，Cu原子与相邻碳原子之间的电荷转移量为[具体数值1]，键级为[具体数值2]。从电子云分布角度来看，Cu原子的电子云与相邻碳原子的电子云发生了明显的重叠，形成了共价键的特征。这种电荷转移和电子云重叠使得Cu原子与石墨烯之间的结合力增强，从而稳定了体系结构。随着覆盖度逐渐增加到1/8，Cu原子与碳原子之间的化学键依然存在，且电荷转移量和键级略有增加。这是因为随着Cu原子数量的增多，Cu原子之间的相互作用以及Cu原子与石墨烯碳原子之间的相互作用都得到了增强，导致电荷转移更加明显，化学键强度增加。当覆盖度进一步增大到1/6、1/4和1/2时，Cu原子与石墨烯碳原子之间的化学键形成情况更为复杂。在这些高覆盖度下，除了Cu原子与直接相邻的碳原子形成化学键外，由于Cu原子在石墨烯表面的聚集，较远的Cu原子与碳原子之间也会通过电子云的相互作用产生一定的影响。在1/2覆盖度下，通过键布居分析发现，部分Cu原子与多个碳原子之间存在电荷转移，形成了类似于网状的化学键结构。这种结构使得体系的稳定性进一步提高，因为更多的化学键连接增强了Cu原子与石墨烯之间的结合力。与其他研究中关于金属与石墨烯相互作用的键布居分析结果相比，本研究中Cu/石墨烯体系的化学键形成情况具有一致性。在相关研究中，也发现金属原子与石墨烯碳原子之间的电荷转移和化学键形成是影响体系稳定性的重要因素。本研究通过不同覆盖度下的详细分析，进一步揭示了Cu/石墨烯体系中化学键形成的演变规律，为深入理解该体系的稳定性提供了更全面的视角。4.3.2Ag/石墨烯和Au/石墨烯体系的化学键变化对于Ag/石墨烯和Au/石墨烯体系，在不同覆盖度下，化学键的变化与相互作用密切相关。在低覆盖度下，如1/32和1/8覆盖度时，Ag/石墨烯体系中，Ag原子与石墨烯表面的碳原子之间形成了有效的化学键。通过键布居分析计算得到，Ag原子与相邻碳原子之间的电荷转移量为[具体数值3]，键级为[具体数值4]。从电子云分布来看，Ag原子的电子云与相邻碳原子的电子云有一定程度的重叠，表明形成了共价键的特征。这种化学键的形成使得Ag原子能够稳定地吸附在石墨烯表面。然而，由于Ag原子的电负性与碳原子的电负性差异相对较小，电荷转移程度相对较弱，与Cu/石墨烯体系相比，Ag/石墨烯体系中化学键的强度较弱。随着覆盖度逐渐增大，当达到1/6、1/4和1/2覆盖度时，Ag/石墨烯体系中金属与石墨烯体系之间的相互作用逐渐转变为范德华力作用范围。在高覆盖度下，Ag原子在石墨烯表面聚集，Ag原子之间的相互作用增强，而Ag原子与石墨烯碳原子之间的电荷转移逐渐减少，键级降低。在1/2覆盖度下，键布居分析显示Ag原子与碳原子之间的电荷转移量减小到[具体数值5]，键级降低到[具体数值6]。此时，体系的稳定性更多地依赖于Ag原子之间的相互作用以及Ag原子与石墨烯之间的范德华力。Au/石墨烯体系在低覆盖度下，如1/32和1/8覆盖度时，Au原子与石墨烯表面的碳原子之间也形成了一定程度的化学键。键布居分析表明，Au原子与相邻碳原子之间的电荷转移量为[具体数值7]，键级为[具体数值8]。与Ag/石墨烯体系类似，由于Au原子的电负性与碳原子的电负性差异较小，电荷转移程度有限，化学键强度相对较弱。随着覆盖度的增加，当达到1/6、1/4和1/2覆盖度时，Au/石墨烯体系中金属与石墨烯之间的相互作用也逐渐转变为以范德华力为主。在高覆盖度下，Au原子在石墨烯表面聚集，Au原子之间的相互作用增强，而与石墨烯碳原子之间的电荷转移减少，键级降低。在1/2覆盖度下，键布居分析显示Au原子与碳原子之间的电荷转移量减小到[具体数值9]，键级降低到[具体数值10]。此时，体系的稳定性主要由Au原子之间的相互作用以及Au原子与石墨烯之间的范德华力维持。对比Ag/石墨烯和Au/石墨烯体系在不同覆盖度下的化学键变化，发现两者具有相似的趋势。在低覆盖度下，金属与碳原子之间都能形成一定程度的化学键，但由于金属原子的特性，化学键强度相对较弱。随着覆盖度的增加，金属原子之间的相互作用逐渐增强，而与石墨烯碳原子之间的化学键逐渐减弱，相互作用转变为以范德华力为主。这种化学键变化与相互作用的关系对体系的稳定性产生了重要影响，在高覆盖度下，体系的稳定性更多地依赖于金属原子之间的相互作用以及金属原子与石墨烯之间的范德华力，而不是化学键的作用。五、电子结构与性能关联5.1能带结构分析通过对金属/石墨烯体系的能带结构进行深入分析，发现其展现出独特的特征，且能够保持石墨烯导体的性质。这一特性对于理解该体系的电学性能以及其在电子器件领域的潜在应用具有关键意义。首先，从石墨烯的本征能带结构来看，它具有独特的线性色散关系。在倒易空间中，其导带和价带在狄拉克点处相交，形成了无质量狄拉克费米子的特征。这种特殊的能带结构赋予了石墨烯优异的电学性能，电子在其中的迁移率极高，电子可近乎无阻力地移动。当金属原子（Cu、Ag、Au）吸附于石墨烯表面形成金属/石墨烯体系后，能带结构发生了变化，但仍保持导体性质。以Cu/石墨烯体系为例，计算得到的能带结构显示，虽然Cu原子的引入使得石墨烯的电子云分布发生改变，但体系整体的能带依然保持着连续的特征，没有出现明显的能隙。在费米能级附近，电子态密度依然存在，表明体系能够导电。从电子结构角度分析，这是因为Cu原子与石墨烯之间的相互作用主要表现为电荷转移和化学键的形成。Cu原子的部分电子转移到石墨烯的π电子体系中，这种电子转移虽然改变了石墨烯的电子结构，但并没有破坏其整体的导电性。Ag/石墨烯体系的能带结构同样保持了导体性质。在低覆盖度下，Ag原子以顶位吸附在石墨烯表面，其与石墨烯之间的相互作用相对较弱，主要通过范德华力和少量电荷转移相互作用。随着覆盖度的增加，吸附结构转变为桥位，相互作用增强。但无论在何种覆盖度下，Ag/石墨烯体系的能带结构中，费米能级处始终存在电子态密度，表明体系具有良好的导电性。这是因为Ag原子与石墨烯之间的相互作用虽然会改变电子云分布，但并没有在能带结构中引入能隙，电子依然能够在体系中自由移动。Au/石墨烯体系的能带结构也呈现出类似的特征。在低覆盖度时，Au原子主要以顶位吸附，与石墨烯之间的相互作用较弱。随着覆盖度增加，吸附结构转变为桥位，相互作用增强。但整个体系的能带结构依然保持导体性质，费米能级处存在电子态密度。这是由于Au原子与石墨烯之间的相互作用方式决定的，尽管Au原子的电子云与石墨烯的电子云有一定程度的重叠和电荷转移，但并没有改变体系整体的能带连续性，使得体系能够保持良好的导电性。金属/石墨烯体系保持石墨烯导体性质的原因主要有以下几点。金属原子与石墨烯之间的相互作用并没有在能带结构中引入能隙，电子依然能够在体系中自由移动。金属原子与石墨烯之间的电荷转移和化学键形成虽然改变了电子云分布，但并没有破坏石墨烯的本征导电特性。金属原子的电子与石墨烯的π电子体系相互作用，使得体系在费米能级附近始终存在电子态密度，从而保证了体系的导电性。5.2电子轨道相互作用在Mx/石墨烯体系中，电子轨道相互作用是影响体系性质的关键因素之一。从量子力学角度来看，电子轨道的相互作用涉及到电子云的重叠和电子的转移，这一过程对体系的电子结构和化学键形成产生重要影响。对于Cu/石墨烯体系，通过理论计算和分析发现，石墨烯的π轨道与Cu原子的3d和4s轨道发生了显著的相互作用。在低覆盖度下，如1/32覆盖度时，Cu原子的4s轨道电子云与石墨烯的π电子云发生重叠，部分4s电子转移到石墨烯的π电子体系中。这种电子转移导致Cu原子带正电，而石墨烯的π电子体系电子密度增加。从键布居分析可知，这种电子转移使得Cu原子与石墨烯表面的碳原子之间形成了一定强度的化学键，从而增强了体系的稳定性。随着覆盖度的增加，Cu原子之间的相互作用也会影响电子轨道的相互作用。在1/2覆盖度下，由于Cu原子在石墨烯表面的聚集，Cu原子之间的3d轨道也会发生相互作用，形成一定的电子结构。这种Cu原子之间的电子轨道相互作用与Cu原子和石墨烯之间的电子轨道相互作用相互影响，共同决定了体系的电子结构和稳定性。在Ag/石墨烯体系中，石墨烯的π轨道主要与Ag原子的4d和5s轨道发生相互作用。在低覆盖度下，Ag原子的5s轨道电子云与石墨烯的π电子云有一定程度的重叠，存在少量的电子转移。由于Ag原子的电负性与碳原子的电负性差异较小，这种电子转移程度相对较弱。随着覆盖度的增加，Ag原子之间的相互作用增强，4d轨道之间的相互作用也变得更加明显。在高覆盖度下，Ag原子在石墨烯表面聚集形成一定的结构，此时Ag原子之间的4d轨道相互作用以及Ag原子与石墨烯之间的电子轨道相互作用共同决定了体系的稳定性。与Cu/石墨烯体系相比，Ag/石墨烯体系中电子轨道相互作用相对较弱，这导致Ag/石墨烯体系的稳定性相对较低。Au/石墨烯体系中，石墨烯的π轨道与Au原子的5d和6s轨道发生相互作用。在低覆盖度下，Au原子的6s轨道电子云与石墨烯的π电子云发生重叠，存在少量的电子转移。由于Au原子的原子半径较大，电子云相对较为分散，电子转移程度有限。随着覆盖度的增加，Au原子之间的相互作用增强，5d轨道之间的相互作用也逐渐增强。在高覆盖度下，Au原子在石墨烯表面聚集形成一定的结构，此时Au原子之间的5d轨道相互作用以及Au原子与石墨烯之间的电子轨道相互作用对体系的稳定性和电子结构产生重要影响。与Cu/石墨烯体系相比，Au/石墨烯体系中电子轨道相互作用相对较弱，这是导致Au/石墨烯体系稳定性相对较低的重要原因之一。电子轨道相互作用对体系性质的影响主要体现在以下几个方面。它直接影响体系的电子结构，改变了电子云的分布和电子态密度。在Cu/石墨烯体系中，由于电子轨道相互作用导致的电子转移，使得石墨烯的电子云分布发生改变，电子态密度在费米能级附近也发生了变化，从而影响了体系的电学性能。电子轨道相互作用与化学键的形成密切相关。在Cu/石墨烯体系中，较强的电子轨道相互作用使得Cu原子与石墨烯碳原子之间形成了稳定的化学键，增强了体系的稳定性。而在Ag/石墨烯和Au/石墨烯体系中，由于电子轨道相互作用相对较弱，在高覆盖度下化学键逐渐减弱，体系的稳定性更多地依赖于范德华力。电子轨道相互作用还会影响体系的光学、磁学等其他物理性质。在一些研究中发现，金属修饰后的石墨烯体系在光学吸收和发射等方面会出现与原始石墨烯不同的特性，这与电子轨道相互作用导致的电子结构变化密切相关。5.3功函与场发射性能5.3.1金属原子吸附对功函的影响功函作为材料的一个重要物理参数，在电子发射等过程中起着关键作用。它被定义为将一个电子从材料内部移至无穷远处所需要的最小能量，反映了材料表面电子的束缚程度。对于石墨烯而言，其本征功函约为4.6-4.7eV，这一数值决定了石墨烯在电子发射等应用中的基础性能。当金属原子吸附于石墨烯表面时，体系的功函会发生显著变化。在低覆盖度下，金属原子的吸附会降低石墨烯的功函。以Cu/石墨烯体系为例，在1/32覆盖度时，通过理论计算得到体系的功函从石墨烯本征的4.6eV降低至4.2eV。这是因为在低覆盖度下，金属原子与石墨烯之间存在明显的电荷转移。金属原子的电子云与石墨烯的π电子体系相互作用，使得石墨烯表面的电子云分布发生改变。具体来说，金属原子的部分电子转移到石墨烯表面，增加了石墨烯表面的电子密度，从而降低了电子从体系中逸出的能量势垒，导致功函降低。这种电荷转移和电子云分布的改变可以通过差分电荷密度图直观地观察到。在差分电荷密度图中，金属原子与石墨烯表面碳原子之间的电荷聚集区域表明了电荷转移的存在。对于Ag/石墨烯体系，在1/32覆盖度下，体系功函从石墨烯本征值降低至4.3eV。与Cu/石墨烯体系类似，Ag原子的电子云与石墨烯的π电子体系相互作用，导致电荷转移，降低了功函。由于Ag原子的电负性与碳原子的电负性差异相对较小，电荷转移程度相对较弱，因此功函降低的幅度相较于Cu/石墨烯体系略小。Au/石墨烯体系在1/32覆盖度下，功函降低至4.4eV。同样，Au原子与石墨烯之间的电荷转移导致功函降低。但由于Au原子的原子半径较大，电子云相对较为分散，电荷转移程度有限，使得功函降低幅度在三种体系中相对最小。随着覆盖度的增加，金属原子和石墨烯片层的碳原子之间相互作用减弱，金属/石墨烯体系的功函完全由金属的功函而决定。在1/2覆盖度下，Cu/石墨烯体系的功函接近Cu金属的功函，约为4.65eV。这是因为在高覆盖度下，金属原子在石墨烯表面形成了相对密集的分布，金属原子之间的相互作用增强，而金属原子与石墨烯碳原子之间的相互作用相对减弱。此时，体系的电子结构更多地受到金属原子之间相互作用的影响，导致功函趋近于金属的功函。Ag/石墨烯体系在1/2覆盖度下，功函接近Ag金属的功函，约为4.26eV。随着覆盖度的增加，Ag原子之间的相互作用增强，与石墨烯碳原子之间的相互作用减弱，使得体系功函逐渐由Ag金属的功函主导。Au/石墨烯体系在1/2覆盖度下，功函接近Au金属的功函，约为5.1eV。在高覆盖度下，Au原子之间的相互作用以及与石墨烯碳原子之间相互作用的变化，导致体系功函趋近于Au金属的功函。5.3.2功函变化与场发射性能的关系功函的变化与金属/石墨烯体系的场发射性能之间存在着紧密的内在联系。场发射是指在强电场作用下，电子从材料表面逸出的现象，其发射电流密度与功函密切相关。根据Fowler-Nordheim理