In、Tl在Si(111)面吸附特性的第一性原理探究：电子结构与界面相互作用

In、Tl在Si(111)面吸附特性的第一性原理探究：电子结构与界面相互作用一、引言1.1研究背景与意义在材料科学与半导体领域，对原子在材料表面的吸附特性研究始终占据着关键地位，它是探索材料微观结构与宏观性能关系的重要切入点。In（铟）和Tl（铊）作为Ⅲ主族金属元素，其在Si(111)面的吸附特性研究备受关注，这一研究方向蕴含着丰富的科学内涵与广阔的应用前景。硅（Si）材料凭借其卓越的半导体性能、成熟的制备工艺以及良好的稳定性，在半导体产业中一直扮演着中流砥柱的角色，从早期的晶体管到如今复杂的大规模集成电路，Si材料无处不在，支撑着信息技术的飞速发展。而Si(111)面由于其独特的原子排列和电子结构，呈现出与其他晶面不同的物理化学性质，使其成为研究表面吸附现象的理想平台。在Si(111)面上，原子的排列方式形成了特定的吸附位点和表面能分布，这为外来原子的吸附提供了多样化的可能性，也使得对其吸附特性的研究充满挑战与机遇。In作为一种重要的金属元素，在半导体领域展现出诸多独特的优势。将In引入Si材料体系中，能够对Si的电学性能进行精准调控。通过在Si(111)面吸附In原子，可以改变Si表面的电子态密度和能带结构，从而影响载流子的浓度和迁移率，为制备高性能的半导体器件奠定基础。例如，在一些新型的场效应晶体管（FET）结构中，利用In在Si(111)面的吸附来优化沟道材料的电学性能，有望实现更低的功耗和更高的开关速度，提升器件的整体性能。此外，In的吸附还可能导致Si表面形成量子点或量子线等低维结构，这些低维结构由于量子限域效应，呈现出与体材料截然不同的光学和电学特性，在光电器件领域具有潜在的应用价值，如用于制备高效的发光二极管（LED）和激光二极管（LD）等。近年来，Ⅲ族最重的元素Tl在Si(111)面吸附的独特性质更是引起了科研人员的极大兴趣。研究发现，当Tl吸附在Si(111)表面时，在直流电压极性转换的作用下，表面结构会在(1×1)和(√3×√3)之间发生可逆的结构改变，这一现象是由Tl原子的表面电迁移（SE）所引起的。这种可逆的表面电迁移特性在其他材料体系中较为罕见，其背后的物理机制尚未完全明晰，但却蕴含着巨大的潜在应用价值。例如，基于这种可逆的结构变化，可以设想开发新型的传感器件，用于检测微弱的电场信号或作为信息存储的新介质，通过表面结构的变化来编码和存储信息，为信息存储技术的发展开辟新的道路。此外，Tl在Si(111)面电迁移后形成具有(4×1)结构的表面硅化物，尽管不具有可逆性，但这种特殊的硅化物结构可能具有独特的电学、光学或催化性能，在半导体器件、光催化等领域展现出潜在的应用前景。从材料表面改性的角度来看，研究In、Tl在Si(111)面的吸附特性为实现Si材料表面性能的优化提供了新的途径。通过控制In、Tl原子的吸附量、吸附位置和吸附结构，可以有针对性地改变Si表面的化学活性、润湿性、耐磨性等性能。在一些微纳加工工艺中，利用In、Tl的吸附来改善Si表面的刻蚀选择性和均匀性，有助于提高器件的制备精度和良品率。在传感器领域，通过表面改性可以增强Si基传感器对特定气体分子的吸附能力和选择性，提高传感器的灵敏度和响应速度，实现对环境中痕量有害气体的快速检测。综上所述，深入研究In、Tl在Si(111)面的吸附特性，无论是对于理解材料表面的微观物理化学过程，还是对于推动新型半导体器件的研发以及材料表面改性技术的发展，都具有至关重要的意义。它不仅能够丰富我们对半导体材料表面科学的认知，还为解决实际应用中的技术难题提供了理论支持和创新思路，有望在未来的信息技术、能源技术、环境监测等领域发挥重要作用。1.2国内外研究现状In、Tl在Si(111)面吸附特性的研究一直是材料科学领域的热门课题，吸引了众多国内外科研人员的关注，相关研究取得了一系列重要成果。在国外，早期研究主要聚焦于实验观测与表征。上世纪末，科研人员通过低能电子衍射（LEED）、扫描隧道显微镜（STM）等先进实验技术，对In、Tl在Si(111)面的吸附结构进行了细致研究。通过LEED技术，精确测定了吸附原子在Si(111)表面的晶格结构和取向关系，发现In原子在Si(111)表面会形成多种重构结构，如(4×1)、(2×2)等，不同的覆盖度和生长条件会导致不同的重构形式出现。STM技术则提供了原子级别的表面图像，直观展示了In、Tl原子在Si(111)表面的吸附位置和分布状态，为后续理论研究提供了重要的实验基础。随着计算技术的飞速发展，理论计算在该领域的研究中发挥着越来越重要的作用。基于密度泛函理论（DFT）的第一性原理计算成为研究In、Tl在Si(111)面吸附特性的重要手段。通过第一性原理计算，国外科研团队深入探究了吸附体系的电子结构、吸附能等关键物理量。研究发现，In原子吸附在Si(111)表面后，会显著改变表面的电子态密度分布，在费米能级附近出现新的电子态，这与表面的电学性能密切相关。在研究In原子在Si(111)-(2×2)表面的吸附时，计算结果表明In原子吸附在H3位时体系具有较低的吸附能，结构较为稳定，且In-Si键长约为2.73Å，这与实验测量结果相符。此外，对于Tl在Si(111)面的吸附研究，发现其独特的表面电迁移现象背后涉及复杂的电子-声子相互作用和原子扩散机制，但具体的定量描述仍存在一定的不确定性。在国内，相关研究起步相对较晚，但近年来发展迅速，在理论与实验方面均取得了显著进展。在实验研究方面，国内科研团队利用同步辐射光电子能谱（SRPES）、高分辨透射电子显微镜（HRTEM）等先进设备，对In、Tl在Si(111)面吸附后的电子结构和微观结构进行了深入分析。通过SRPES技术，精确测量了吸附体系的电子结合能和能带结构，进一步验证了理论计算中关于电子态变化的预测。HRTEM技术则提供了吸附原子与Si(111)衬底之间原子尺度的界面结构信息，揭示了吸附原子在衬底表面的扩散路径和聚集形态。在理论研究方面，国内学者基于第一性原理计算，结合分子动力学模拟等方法，从多个角度对In、Tl在Si(111)面的吸附特性进行了全面研究。不仅深入探讨了吸附能、电荷转移、键长键角等基本物理量，还研究了温度、压力等外部条件对吸附特性的影响。在研究温度对In在Si(111)面吸附稳定性的影响时，通过分子动力学模拟发现，随着温度的升高，In原子在表面的扩散系数增大，吸附结构的稳定性逐渐降低，这为实际应用中材料的制备和性能调控提供了重要的理论依据。此外，国内研究团队还针对In、Tl在Si(111)面吸附形成的低维结构的光学和电学性能开展了大量研究，为新型光电器件的设计和开发提供了理论指导。尽管国内外在In、Tl在Si(111)面吸附特性的研究方面取得了丰硕的成果，但仍存在一些不足之处。在实验研究中，由于实验条件的限制，对于一些极端条件下（如高温、高压、超快过程）的吸附特性研究还不够深入，实验数据的准确性和完整性有待进一步提高。在理论计算方面，虽然现有的计算方法能够对吸附体系的一些基本性质进行较好的预测，但对于一些复杂的多体相互作用和量子效应的描述还存在一定的近似性，导致计算结果与实验数据之间存在一定的偏差。此外，目前对于In、Tl在Si(111)面吸附特性的研究主要集中在静态性质方面，对于吸附过程中的动态演化过程（如原子的扩散、反应动力学等）的研究还相对较少，这限制了对吸附机理的全面理解。本文将在现有研究的基础上，采用高精度的第一性原理计算方法，结合先进的实验技术，深入研究In、Tl在Si(111)面的吸附特性。重点关注吸附过程中的电子结构变化、原子间相互作用以及动态演化过程，旨在揭示In、Tl在Si(111)面吸附的微观机理，为新型半导体材料的设计和应用提供更加坚实的理论基础。1.3研究目标与内容本文旨在通过基于密度泛函理论的第一性原理计算方法，深入探究In、Tl在Si(111)面的吸附特性，揭示其微观吸附机理，为半导体材料表面改性和新型器件设计提供坚实的理论依据。具体研究内容涵盖以下几个方面：吸附能与稳定性分析：构建不同覆盖度下In、Tl在Si(111)面的多种吸附结构模型，利用第一性原理精确计算各模型的吸附能。通过对比吸附能的大小，确定In、Tl在Si(111)面的最稳定吸附结构和吸附位点，分析覆盖度对吸附稳定性的影响规律。如在研究In原子在Si(111)-(2×2)表面的吸附时，需全面考虑不同吸附位点（如顶位、桥位、穴位等）的吸附能，找出能量最低的稳定吸附构型。同时，对于Tl原子在不同超原胞（如(2×2)、(√3×√3)等）表面的吸附，也需进行细致的吸附能计算和稳定性分析，为后续研究提供基础。电子结构分析：深入研究In、Tl吸附前后Si(111)面的电子结构变化，包括电子态密度（DOS）、差分电荷密度、能带结构等。通过分析电子态密度，明确吸附原子对Si(111)表面电子分布的影响，确定新出现的表面态及其与费米能级的相对位置。利用差分电荷密度图，直观展示吸附原子与Si原子之间的电荷转移和化学键形成情况，揭示吸附体系的成键本质。在分析In吸附体系的电子态密度时，关注In原子的价电子与Si表面电子的相互作用，以及由此导致的表面态变化，探讨其对体系电学性能的影响。对于Tl吸附体系，重点研究其独特的表面电迁移过程中电子结构的动态变化，以及这种变化与表面结构改变之间的内在联系。晶格结构变化研究：精确计算In、Tl吸附后Si(111)面的晶格参数、键长、键角等结构参数的变化，分析吸附原子对Si(111)表面晶格结构的重构作用。通过模拟不同覆盖度下的吸附体系，观察晶格结构的演变规律，探讨晶格结构变化与吸附稳定性、电子结构之间的相互关系。在研究In吸附导致的晶格结构变化时，测量In-Si键长、Si-Si键长和键角的改变，分析这些变化对表面原子排列和周期性的影响。对于Tl在Si(111)面电迁移过程中的晶格结构变化，需实时跟踪和分析，揭示电迁移过程中晶格结构的动态演化机制。温度和电场对吸附特性的影响：引入温度和电场等外部条件，利用分子动力学模拟和第一性原理相结合的方法，研究温度和电场对In、Tl在Si(111)面吸附特性的影响。在不同温度下，模拟吸附原子的扩散行为和吸附结构的稳定性变化，分析温度对吸附能、电子结构和晶格结构的影响机制。施加不同强度和方向的电场，模拟Tl原子在Si(111)面的电迁移过程，研究电场对表面结构转变和电子结构变化的影响规律，揭示Tl在Si(111)面独特的电学性质与电场的内在联系。如在研究温度对In吸附稳定性的影响时，通过分子动力学模拟，观察In原子在不同温度下的扩散路径和扩散系数，分析温度升高导致吸附结构失稳的微观过程。在研究电场对Tl电迁移的影响时，精确控制电场强度和方向，模拟Tl原子的迁移轨迹和表面结构的转变过程，为基于Tl吸附的新型器件设计提供理论指导。1.4研究方法与技术路线本研究采用基于密度泛函理论（DFT）的第一性原理计算方法，借助VASP（ViennaAb-initioSimulationPackage）软件开展深入研究。该方法基于量子力学原理，从电子与原子核的相互作用出发，在不依赖任何经验参数的情况下，精确求解多粒子体系的薛定谔方程，从而获得体系的电子结构和能量等重要物理性质，为深入理解In、Tl在Si(111)面的吸附特性提供了坚实的理论基础。技术路线主要包含以下关键步骤：模型构建：运用MaterialsStudio软件构建Si(111)面的超原胞模型，通过合理设置衬底层数和真空层厚度，确保模型能够准确反映实际表面的物理性质。针对In、Tl在Si(111)面的吸附体系，考虑不同覆盖度和多种可能的吸附位点，构建相应的吸附结构模型。在构建Si(111)-(2×2)超原胞模型时，衬底选取10层Si原子，真空层厚度设置为15Å，以有效避免周期性边界条件下相邻超原胞之间的相互作用。对于In原子在该表面的吸附，全面考虑顶位、桥位、穴位等不同吸附位点，构建多种吸附结构模型，为后续计算提供丰富的研究对象。参数设置：将构建好的模型导入VASP软件中，对计算参数进行精细设置。选用广义梯度近似（GGA）下的PBE（Perdew-Burke-Ernzerhof）泛函来描述电子间的交换关联作用，该泛函在处理固体材料的电子结构和能量问题时具有较高的精度和广泛的适用性。平面波截断能设置为500eV，以保证电子波函数的展开精度，确保计算结果的准确性。采用Monkhorst-Pack方法对布里渊区进行K点采样，对于(2×2)超原胞，布里渊区用(4×4×1)个K点进行积分；对于(√3×√3)R30°超原胞，布里渊区用(7×7×1)个K点进行积分，以准确描述体系的电子态分布。结构优化：在设定好计算参数后，对Si(111)面的超原胞模型以及In、Tl吸附体系的模型进行结构优化。通过不断调整原子的位置和晶格参数，使体系的总能量达到最低，此时得到的结构即为体系的稳定结构。在优化过程中，采用共轭梯度算法，以快速收敛到能量极小值点，提高计算效率。结构优化的收敛标准设置为原子间作用力小于0.01eV/Å，总能量变化小于10⁻⁵eV，确保优化后的结构具有较高的稳定性和可靠性。性质计算：对优化后的稳定结构进行吸附能、电子结构和晶格结构等性质的计算。吸附能的计算通过公式E_ad=E_total-E_sub-nE_atom来实现，其中E_ad表示吸附能，E_total为吸附体系的总能量，E_sub为Si(111)衬底的能量，n为吸附原子的个数，E_atom为单个吸附原子的能量。通过比较不同吸附结构的吸附能，确定最稳定的吸附结构和吸附位点。在计算电子结构时，利用VASP软件输出的电子态密度（DOS）数据，绘制体系的总态密度和分波态密度图，分析吸附原子对Si(111)表面电子分布的影响；通过计算差分电荷密度，直观展示吸附原子与Si原子之间的电荷转移和化学键形成情况。对于晶格结构的计算，精确获取吸附后Si(111)面的晶格参数、键长、键角等结构参数，分析吸附原子对表面晶格结构的重构作用。结果分析与讨论：对计算得到的吸附能、电子结构和晶格结构等结果进行深入分析与讨论。结合相关理论知识和实验数据，探讨In、Tl在Si(111)面的吸附特性和微观吸附机理。通过对比不同覆盖度下的吸附能，研究覆盖度对吸附稳定性的影响规律；从电子态密度和差分电荷密度的变化，分析吸附原子与Si表面之间的电子相互作用和化学键本质；根据晶格结构参数的变化，揭示吸附原子对Si(111)表面晶格结构的影响机制。同时，将计算结果与已有实验数据和理论研究进行对比验证，进一步完善对In、Tl在Si(111)面吸附特性的认识。二、第一性原理计算方法基础2.1第一性原理的基本概念第一性原理，在计算物理与计算化学领域占据着核心地位，它的起源可以追溯到量子力学的基本理论。从本质上讲，第一性原理是一种基于量子力学原理的计算方法，其核心在于从最基本的物理定律出发，对多粒子体系进行精确的理论描述，而不依赖任何经验参数。这种计算方法将由多个原子构成的体系看作是一个由电子和原子核组成的复杂系统，通过直接求解薛定谔方程来获取体系的各种物理性质，如能量、电子结构、几何结构等。在多电子体系中，电子与原子核之间存在着复杂的相互作用，包括静电相互作用、交换关联作用等。薛定谔方程作为量子力学的基本方程，能够准确描述这些相互作用下体系的量子态和能量本征值。然而，对于实际的多原子体系，由于电子数目众多，直接求解薛定谔方程面临着巨大的计算挑战，其计算量随着电子数目的增加呈指数级增长，这使得精确求解变得极为困难。为了克服这一难题，第一性原理计算在实际应用中通常会采用一些合理的近似方法。这些近似方法在保证计算精度的前提下，有效降低了计算复杂度，使得大规模体系的计算成为可能。其中，最为重要的近似之一就是密度泛函理论（DFT）。DFT的核心思想是将体系的能量表示为电子密度的泛函，通过求解关于电子密度的方程来间接获取体系的能量和其他物理性质。这种方法的优势在于，将多电子问题转化为单电子问题，大大简化了计算过程。在DFT中，体系的能量被分为多个部分，包括电子的动能、电子与原子核的相互作用能、电子之间的库仑相互作用能以及交换关联能等。其中，交换关联能是DFT中最为关键且复杂的部分，它描述了电子之间由于交换和关联效应而产生的相互作用。为了准确描述交换关联能，科研人员提出了多种交换关联泛函，如局域密度近似（LDA）、广义梯度近似（GGA）等。这些泛函在不同程度上对交换关联能进行了近似处理，在实际应用中取得了较好的效果。除了DFT，平面波赝势方法（PWPM）也是第一性原理计算中常用的技术手段。在晶体体系中，电子的波函数可以用平面波基组进行展开。然而，由于原子核周围存在着强的势场，导致电子波函数在原子核附近变化剧烈，需要大量的平面波才能准确描述。这不仅增加了计算量，还会带来数值计算上的困难。为了解决这一问题，PWPM引入了赝势的概念。赝势通过对原子核与内层电子的相互作用进行有效屏蔽，使得在计算中可以用较少的平面波来描述电子的波函数，从而显著提高了计算效率。在PWPM中，赝势的构造是关键环节，它需要在保证计算精度的前提下，尽可能地简化计算过程。目前，常用的赝势包括模守恒赝势、超软赝势等，它们在不同的体系中都有着广泛的应用。第一性原理计算以其从基本物理定律出发、不依赖经验参数的独特优势，为研究物质的微观结构和性质提供了强有力的工具。尽管在实际计算中需要采用一些近似方法，但这些方法在不断的发展和完善中，使得第一性原理计算能够在材料科学、化学、物理等多个领域中发挥重要作用，为深入理解物质的本质和开发新型材料提供了坚实的理论基础。2.2密度泛函理论（DFT）密度泛函理论（DensityFunctionalTheory，DFT）作为第一性原理计算的核心理论框架，在现代材料科学和凝聚态物理研究中占据着举足轻重的地位。其核心思想在于将多电子体系的基态性质表述为电子密度的泛函，通过求解与电子密度相关的方程，进而获得体系的能量、电子结构等关键物理量，成功地将复杂的多电子问题转化为相对简单的关于电子密度函数的问题。在传统的量子力学方法中，描述多电子体系需要考虑电子之间复杂的相互作用，其波函数是所有电子坐标的函数，这使得计算量随着电子数目的增加呈指数级增长，在实际应用中面临巨大挑战。而DFT的出现为解决这一难题提供了新的思路，它将电子密度作为核心变量，体系的总能量可以表示为电子密度的泛函形式：E[\rho]=T[\rho]+V_{ne}[\rho]+V_{ee}[\rho]+E_{xc}[\rho]。其中，T[\rho]表示电子的动能泛函，描述了电子的运动状态；V_{ne}[\rho]是电子与原子核之间的相互作用能泛函，体现了电子与原子核的库仑吸引作用；V_{ee}[\rho]为电子间的库仑相互作用能泛函，反映了电子之间的静电排斥力；E_{xc}[\rho]则是交换关联能泛函，它涵盖了电子之间由于交换效应和关联效应产生的复杂相互作用。交换效应源于电子的全同性，使得相同自旋的电子相互回避；关联效应则描述了电子之间由于库仑相互作用导致的动态相关行为。交换关联能泛函是DFT中最为关键且难以精确描述的部分，其准确程度直接影响着计算结果的可靠性。为了准确描述交换关联能，科研人员提出了多种近似泛函，其中局域密度近似（LocalDensityApproximation，LDA）和广义梯度近似（GeneralizedGradientApproximation，GGA）是最为常用的两种。LDA假设体系中某点的交换关联能只与该点的电子密度有关，采用均匀电子气模型来近似描述交换关联能。虽然LDA在处理一些简单体系时能够取得较好的结果，如金属晶体等，但对于电子密度变化较为剧烈的体系，如分子体系和表面体系，其计算精度会受到一定限制。GGA则在LDA的基础上进行了改进，不仅考虑了电子密度的局域值，还引入了电子密度的梯度信息。GGA能够更好地描述电子密度的非均匀性，对于复杂体系的计算精度有了显著提高，在许多实际应用中表现出比LDA更优越的性能。例如，在研究半导体材料的电子结构时，GGA能够更准确地预测能带结构和能隙大小，与实验结果的吻合度更高。随着对材料微观结构和性能研究的深入，传统的LDA和GGA泛函在某些情况下仍然无法满足高精度计算的需求。为了进一步提高DFT的计算精度，科研人员不断探索新的交换关联泛函形式。混合泛函（HybridFunctional）将Hartree-Fock交换能与DFT交换关联能相结合，通过调整混合比例，能够在一定程度上改善对体系电子结构的描述，尤其在处理分子体系的激发态和能隙问题上表现出独特的优势。meta-GGA泛函则在GGA的基础上，引入了动能密度等更多的物理量，以更全面地描述电子之间的相互作用，在一些复杂材料体系的计算中展现出良好的应用前景。DFT在材料计算领域有着极为广泛的应用。在研究材料的晶体结构时，通过计算不同原子排列方式下体系的能量，能够确定材料的最稳定晶体结构和晶格参数。在研究半导体材料的电学性能时，利用DFT计算电子态密度和能带结构，可以深入了解材料的导电机制、载流子迁移率等关键电学性质。在研究催化剂的催化性能时，DFT能够计算反应物在催化剂表面的吸附能、反应路径和反应能垒，为揭示催化反应机理和设计高效催化剂提供理论依据。例如，在研究金属催化剂对化学反应的催化作用时，通过DFT计算可以确定反应物在催化剂表面的最佳吸附位点和反应活性中心，从而指导催化剂的优化设计。密度泛函理论以其独特的理论优势和广泛的应用价值，为研究多电子体系的性质提供了强有力的工具。尽管在交换关联能的描述等方面仍存在一定的挑战，但随着理论的不断发展和完善，DFT将在材料科学、化学、物理等多个领域发挥更加重要的作用，为推动相关领域的技术创新和发展提供坚实的理论支撑。2.3赝势方法在第一性原理计算中，赝势方法是一种至关重要的近似手段，它的引入极大地简化了离子实与价电子之间相互作用的描述，显著提高了计算效率，为研究复杂体系的电子结构和性质提供了可能。在真实的晶体体系中，原子核周围存在着由内层电子和原子核构成的离子实，离子实与价电子之间存在着强相互作用，这种相互作用使得描述价电子的波函数在离子实附近变化极为剧烈，需要大量的平面波才能准确描述。这不仅导致计算量呈指数级增长，还会带来数值计算上的困难，严重限制了第一性原理计算在大规模体系中的应用。赝势方法的核心思想是通过构建一个等效的赝势来代替真实的离子实势，使得在求解薛定谔方程时，能够用相对较少的平面波来描述价电子的波函数。在构建赝势时，需要满足两个关键条件：一是在离子实区域外，赝势所产生的波函数与真实势下的波函数相同，以保证对体系外部性质的准确描述；二是在离子实区域内，赝势能够使波函数变得平滑，消除真实势下波函数的剧烈振荡和大量节点，从而减少平面波的使用数量。通过满足这两个条件，赝势成功地将离子实与价电子之间复杂的相互作用进行了有效简化，使得计算过程更加高效和可行。在众多赝势类型中，超软赝势（Ultra-SoftPseudopotential，US-PP）和投影缀加波（ProjectorAugmented-Wave，PAW）势是目前应用最为广泛的两种。超软赝势的主要特点在于其对波函数的平滑处理更加彻底，它允许波函数在离子实区域内有一定的变化，从而进一步减少了平面波的数量需求。这使得在计算大体系时，超软赝势能够显著降低计算量，提高计算效率。在研究包含大量原子的复杂晶体结构时，使用超软赝势可以在保证一定计算精度的前提下，大大缩短计算时间，提高研究效率。然而，由于超软赝势对波函数的处理较为灵活，其计算精度相对较低，在一些对精度要求极高的计算中可能无法满足需求。投影缀加波（PAW）势则是在超软赝势的基础上发展而来的一种更为精确的赝势方法。PAW势通过引入投影算子，将全电子波函数和赝波函数联系起来，实现了对全电子波函数的精确重构。这使得PAW势在保持较低计算量的同时，能够达到与全电子计算相近的精度。在研究过渡金属化合物的电子结构时，PAW势能够准确描述过渡金属原子的d电子态，与实验结果具有很好的吻合度。因此，PAW势在对精度要求较高的体系计算中具有明显优势，如在研究材料的电子结构、光学性质、磁性等方面得到了广泛应用。赝势方法的选择应根据具体的研究体系和计算需求来确定。对于计算精度要求不高、体系规模较大的情况，超软赝势是一种较为合适的选择，它能够在较短的时间内给出体系的大致性质。而对于精度要求较高、体系规模相对较小的体系，投影缀加波势则能够提供更为准确的计算结果。在实际应用中，还可以通过对比不同赝势方法的计算结果，结合实验数据进行验证，以确保计算结果的可靠性。例如，在研究In、Tl在Si(111)面的吸附特性时，可以分别采用超软赝势和投影缀加波势进行计算，对比两种方法得到的吸附能、电子结构等结果，结合实验测量的吸附结构和电学性能数据，判断哪种赝势方法更适合该体系的研究。赝势方法作为第一性原理计算中的关键技术，通过简化离子实与价电子的相互作用，为研究复杂体系的电子结构和性质提供了高效的手段。超软赝势和投影缀加波势各具特点和优势，在不同的应用场景中发挥着重要作用。随着计算技术的不断发展和对材料性质研究的深入，赝势方法也将不断完善和创新，为材料科学、物理、化学等领域的研究提供更加坚实的理论支持。2.4计算软件与参数设置本研究选用维也纳从头算模拟软件包（VASP）开展基于第一性原理的计算工作。VASP是一款基于密度泛函理论的优秀计算软件，在材料科学领域有着广泛的应用，能够精确模拟原子和电子体系的性质。它采用平面波赝势方法，将电子的波函数用平面波基组展开，通过构建赝势来有效简化离子实与价电子之间的相互作用，从而在保证计算精度的前提下，显著提高计算效率。在VASP计算中，平面波截断能的设置是一个关键参数。平面波截断能决定了平面波基组的大小，进而影响计算结果的准确性和计算效率。经过一系列的测试计算，本研究将平面波截断能设置为500eV。这一数值的选择是基于对体系能量收敛性的细致考察。在测试过程中，逐步增大平面波截断能，观察体系总能量的变化情况。当平面波截断能较低时，体系总能量随着截断能的增大而迅速下降，这表明此时平面波基组不够完备，无法准确描述体系的电子结构。随着截断能的进一步增大，体系总能量的变化逐渐趋于平缓，当截断能达到500eV时，总能量的变化小于10⁻⁵eV，满足了计算精度的要求。此时，继续增大截断能对计算结果的影响极小，却会显著增加计算时间和计算资源的消耗。因此，综合考虑计算精度和效率，将平面波截断能确定为500eV，以确保在合理的计算成本下获得准确的计算结果。对于K点网格的设置，采用Monkhorst-Pack方法对布里渊区进行K点采样。K点的选取决定了对布里渊区积分的精度，进而影响计算结果的准确性。对于不同的超原胞模型，K点网格的设置有所不同。对于(2×2)超原胞，经过测试分析，选择(4×4×1)个K点进行积分。在测试过程中，改变K点的数量，计算体系的总能量和电子态密度等物理量。当K点数量较少时，计算结果存在较大的误差，尤其是在描述电子态分布时不够准确。随着K点数量的增加，计算结果逐渐收敛，当K点网格设置为(4×4×1)时，计算结果的收敛性良好，能够准确反映体系的电子结构特征。对于(√3×√3)R30°超原胞，经过类似的测试，布里渊区用(7×7×1)个K点进行积分。这样的K点设置能够在保证计算精度的前提下，有效提高计算效率，避免因过多的K点导致计算量过大。在描述电子间的交换关联作用时，选用广义梯度近似（GGA）下的PBE（Perdew-Burke-Ernzerhof）泛函。PBE泛函在考虑电子密度的局域值的基础上，引入了电子密度的梯度信息，能够较好地描述电子密度的非均匀性，对于复杂体系的计算精度较高。在研究半导体材料的电子结构和吸附特性等方面，PBE泛函已被广泛应用，并取得了与实验结果较为吻合的计算结果。与其他泛函相比，如局域密度近似（LDA）泛函，PBE泛函在处理表面体系和分子体系时，能够更准确地描述电子之间的交换关联效应，从而更精确地预测体系的能量、电子结构和几何结构等性质。在赝势的选择上，采用投影缀加波（PAW）势。PAW势通过引入投影算子，实现了对全电子波函数的精确重构，在保持较低计算量的同时，能够达到与全电子计算相近的精度。对于In、Tl在Si(111)面的吸附体系，PAW势能够准确描述In、Tl原子与Si原子之间的相互作用，以及吸附前后体系电子结构的变化。与超软赝势相比，PAW势在处理过渡金属和重元素时具有明显优势，能够更准确地描述这些原子的电子态，从而为研究In、Tl在Si(111)面的吸附特性提供更可靠的计算结果。通过对VASP软件关键参数的合理设置，本研究能够在保证计算精度的前提下，高效地开展In、Tl在Si(111)面吸附特性的第一性原理计算，为后续深入分析吸附体系的性质奠定坚实的基础。三、Si(111)表面模型构建与优化3.1Si晶体结构与Si(111)面的特点硅（Si）晶体具有典型的金刚石结构，这种结构在晶体学中具有独特的原子排列方式和重要的物理意义。在金刚石结构中，Si原子通过共价键相互连接，形成了一个三维的网状结构。每个Si原子与周围四个Si原子以共价键结合，键角为109.47°，这种特定的键角和键长决定了晶体的空间构型。Si晶体的晶格常数约为5.43Å，这一参数对于描述晶体的几何特征和原子间距离起着关键作用。在这种结构中，Si原子的排列呈现出高度的对称性和周期性，每个晶胞包含8个Si原子，其中4个原子位于面心立方的顶点和面心位置，另外4个原子则位于晶胞内部的特定位置，形成了一种嵌套的结构。这种结构使得Si晶体具有良好的稳定性和半导体特性，为其在半导体领域的广泛应用奠定了基础。Si(111)面作为Si晶体的一个重要晶面，具有独特的原子排列方式和表面特性。从原子排列角度来看，Si(111)面的原子呈六边形密堆积排列，形成了类似于蜂窝状的结构。在理想的Si(111)面中，原子排列整齐，具有高度的周期性。然而，实际的Si(111)表面由于原子的弛豫和重构现象，其结构会发生一定的变化。表面原子的弛豫是指表面原子相对于体内原子的位置发生微小的位移，这种位移主要是由于表面原子所受到的原子间作用力与体内原子不同所致。表面原子在垂直于表面方向上的位移最为明显，导致表面原子层与内层原子层之间的间距发生改变。这种弛豫现象会影响表面的电子结构和化学活性。例如，表面原子的弛豫会导致表面电子云分布的变化，进而影响表面与外来原子的相互作用。重构则是指表面原子的重新排列，形成与体内原子排列不同的二维结构。Si(111)面在不同的条件下会出现多种重构形式，如(2×1)重构、(7×7)重构等。在低温下，Si(111)面通常会出现(2×1)重构，这种重构是由于表面原子之间的相互作用使得原子形成了新的排列方式。在(2×1)重构结构中，表面原子形成了交替排列的二聚体结构，这种结构的出现降低了表面的能量，使表面更加稳定。而在加热到400°C以上时，Si(111)面会出现(7×7)重构。(7×7)重构的结构较为复杂，其原胞中包含13个空格点（空位）。这种重构结构的形成是由于表面原子在高温下具有更高的能量，能够克服原子间的势垒，重新排列形成新的稳定结构。(7×7)重构结构对Si(111)面的电子结构和化学性质产生了显著的影响，使得表面具有独特的电学和光学性质。Si(111)面的表面特性还包括表面能和表面电荷分布等方面。表面能是指增加单位面积表面所需要的能量，Si(111)面的表面能相对较高，这使得表面具有较高的化学活性，容易与外来原子发生相互作用。表面电荷分布则与表面原子的电子结构和重构形式密切相关。在重构后的Si(111)面上，原子的电子云分布发生变化，导致表面电荷分布不均匀。这种不均匀的电荷分布会影响表面与外来原子之间的静电相互作用，进而影响外来原子在表面的吸附行为。例如，在某些重构结构中，表面会出现局部的电荷富集或贫化区域，这些区域对外来原子具有不同的吸附亲和力，从而影响吸附原子的位置和吸附稳定性。Si(111)面的原子排列方式和表面特性使其成为研究表面吸附现象的理想体系。In、Tl等原子在Si(111)面的吸附行为与Si(111)面的这些特性密切相关，深入研究Si(111)面的结构和性质，对于理解In、Tl在Si(111)面的吸附特性和微观吸附机理具有重要的意义。3.2超原胞模型的选择与构建在研究In、Tl在Si(111)面的吸附特性时，构建合理的超原胞模型是进行第一性原理计算的关键步骤。由于实际晶体表面的原子排列和电子结构与体相存在显著差异，采用超原胞模型能够有效模拟表面体系，准确描述吸附原子与表面原子之间的相互作用。超原胞模型不仅可以考虑表面原子的弛豫和重构等现象，还能通过调整超原胞的大小和形状，精确控制吸附原子的覆盖度，为深入研究吸附特性提供了有力的工具。本研究选用(2×2)和(√3×√3)R30°两种超原胞模型来模拟Si(111)表面。(2×2)超原胞模型是一种较为常见且简单的模型，其原胞尺寸在表面方向上为Si(111)面原胞尺寸的2倍，这种模型能够较好地描述表面原子的周期性排列，同时在计算资源的消耗上相对较低，适合初步研究吸附原子在表面的吸附位点和吸附能等基本性质。(√3×√3)R30°超原胞模型则具有更复杂的原子排列方式，其原胞形状和原子位置的特殊性使得它能够捕捉到Si(111)表面一些更为精细的结构特征。这种超原胞模型在研究吸附原子导致的表面重构和电子结构变化等方面具有独特的优势，能够提供更丰富的信息。在研究Tl原子在Si(111)面的吸附时，(√3×√3)R30°超原胞模型能够更准确地描述由于Tl原子吸附和电迁移导致的表面结构在(1×1)和(√3×√3)之间的转变过程。对于超原胞模型中衬底的层数，经过一系列测试计算，最终选取10层Si原子作为衬底。在测试过程中，逐步增加衬底层数，观察体系总能量和吸附能等物理量的变化。当衬底层数较少时，由于表面原子与衬底内部原子的相互作用未得到充分考虑，体系的总能量和吸附能计算结果存在较大波动，无法准确反映吸附体系的真实性质。随着衬底层数的增加，体系的能量逐渐收敛。当衬底层数达到10层时，继续增加层数对体系总能量和吸附能的影响小于10⁻⁵eV，满足了计算精度的要求。此时，10层衬底能够充分考虑表面原子与衬底内部原子的相互作用，确保计算结果的准确性。为了避免周期性边界条件下相邻超原胞之间的相互作用对计算结果产生干扰，需要在超原胞模型中设置合适的真空层厚度。经过多次测试，将真空层厚度设置为15Å。在测试过程中，改变真空层厚度，计算体系的总能量和电子结构等物理量。当真空层厚度较小时，相邻超原胞之间的电子云会发生重叠，导致计算结果出现偏差。随着真空层厚度的增加，相邻超原胞之间的相互作用逐渐减弱。当真空层厚度达到15Å时，相邻超原胞之间的相互作用可以忽略不计，体系的计算结果不再受真空层厚度变化的影响。此时，15Å的真空层能够有效隔离相邻超原胞，保证计算结果的可靠性。通过合理选择(2×2)和(√3×√3)R30°超原胞模型，确定10层Si原子作为衬底以及15Å的真空层厚度，本研究构建了准确可靠的Si(111)表面超原胞模型，为后续深入研究In、Tl在Si(111)面的吸附特性奠定了坚实的基础。3.3结构优化与晶格常数计算在完成Si(111)表面超原胞模型的构建后，对模型进行结构优化是获取准确物理性质的关键步骤。本研究采用共轭梯度算法对模型进行结构优化，该算法通过迭代更新原子的位置和晶格参数，使体系的总能量逐步降低，直至收敛到最小值。在优化过程中，严格遵循设定的收敛标准，即原子间作用力小于0.01eV/Å，总能量变化小于10⁻⁵eV。当体系满足这两个收敛条件时，表明原子已达到平衡位置，此时得到的结构即为体系的稳定结构。这种严格的收敛标准确保了优化后的结构具有较高的稳定性和可靠性，能够准确反映体系的真实状态。在结构优化过程中，平面波截断能和K点网格的设置对计算结果的准确性和计算效率有着重要影响。如前文所述，本研究将平面波截断能设置为500eV，这一数值是在对体系能量收敛性进行细致考察后确定的。当平面波截断能较低时，平面波基组不够完备，无法准确描述体系的电子结构，导致体系总能量较高且波动较大。随着截断能的增大，体系总能量逐渐下降并趋于稳定。当截断能达到500eV时，总能量的变化小于10⁻⁵eV，满足了计算精度的要求。继续增大截断能虽然能进一步提高计算精度，但会显著增加计算时间和计算资源的消耗，因此综合考虑计算精度和效率，选择500eV作为平面波截断能。对于K点网格的设置，采用Monkhorst-Pack方法对布里渊区进行K点采样。对于(2×2)超原胞，经过测试分析，选择(4×4×1)个K点进行积分。当K点数量较少时，计算结果存在较大的误差，尤其是在描述电子态分布时不够准确。随着K点数量的增加，计算结果逐渐收敛。当K点网格设置为(4×4×1)时，计算结果的收敛性良好，能够准确反映体系的电子结构特征。对于(√3×√3)R30°超原胞，经过类似的测试，布里渊区用(7×7×1)个K点进行积分。这样的K点设置能够在保证计算精度的前提下，有效提高计算效率，避免因过多的K点导致计算量过大。在选择赝势方法时，本研究采用投影缀加波（PAW）势。PAW势通过引入投影算子，实现了对全电子波函数的精确重构，在保持较低计算量的同时，能够达到与全电子计算相近的精度。为了验证PAW势的计算效果，将其计算结果与其他赝势方法（如超软赝势）进行对比。在计算Si晶体的晶格常数时，发现PAW势计算得到的晶格常数与实验值更为接近。采用PAW势计算得到的Si晶体晶格常数为5.42Å，而实验值为5.43Å，相对误差仅为0.18%。相比之下，超软赝势计算得到的晶格常数与实验值的偏差较大，这表明PAW势在描述Si晶体的电子结构和晶格性质方面具有更高的准确性。通过对Si(111)表面超原胞模型进行结构优化，采用合适的计算参数和赝势方法，本研究获得了准确可靠的稳定结构，为后续深入研究In、Tl在Si(111)面的吸附特性奠定了坚实的基础。四、In在Si(111)面的吸附特性研究4.1吸附位点与吸附结构在研究In在Si(111)面的吸附特性时，首要任务是确定其在Si(111)表面可能的吸附位点，并深入分析不同吸附结构的几何构型。Si(111)表面原子呈六边形密堆积排列，为In原子提供了多种潜在的吸附位点，主要包括顶位（Topsite）、桥位（Bridgesite）和穴位（Hollowsite），其中穴位又可细分为H1位（中心穴位）、H2位（次近邻穴位）和H3位（最近邻穴位）。这些不同的吸附位点具有不同的原子环境和配位情况，从而导致In原子在不同位点吸附时形成的吸附结构和稳定性存在显著差异。对于(2×2)超原胞模型，我们对In原子在上述各个吸附位点的吸附结构进行了细致的优化和分析。在顶位吸附结构中，In原子直接位于Si(111)表面的一个Si原子正上方，In原子与下方Si原子形成垂直的键合关系。此时，In-Si键长约为2.65Å，In原子与周围Si原子的距离相对较大，这种吸附结构使得In原子在表面的暴露程度较高，表面原子的电子云分布受到In原子的影响较为明显。从原子间相互作用的角度来看，顶位吸附时In原子主要与下方的单个Si原子发生较强的相互作用，而与周围其他Si原子的相互作用相对较弱，导致吸附体系的稳定性相对较低。在桥位吸附结构中，In原子位于两个相邻Si原子的中间位置，形成类似于桥梁的结构。优化后的In-Si键长约为2.70Å，此时In原子与两个相邻Si原子的距离相等，In原子与周围Si原子的相互作用相对较为均匀。桥位吸附结构的稳定性略高于顶位吸附结构，这是因为In原子通过与两个Si原子形成化学键，增加了吸附体系的结合力，使得体系的能量相对降低。从电子结构的角度分析，桥位吸附时In原子的电子云与两个Si原子的电子云发生重叠，形成了相对稳定的化学键，从而增强了吸附体系的稳定性。在穴位吸附结构中，由于穴位的种类不同，In原子的吸附结构和稳定性也有所差异。在H1位（中心穴位）吸附时，In原子位于由六个Si原子围成的正六边形中心，与周围六个Si原子的距离相等。优化后的In-Si键长约为2.75Å，这种吸附结构使得In原子与周围Si原子的相互作用较为均衡，体系的稳定性相对较高。从空间结构的角度来看，H1位吸附时In原子处于表面原子形成的六边形空隙中心，周围Si原子对In原子形成了一定的包围作用，使得吸附体系具有较好的空间对称性，从而增强了体系的稳定性。在H2位（次近邻穴位）吸附时，In原子位于由四个Si原子围成的四边形中心，与周围四个Si原子的距离相对较短。优化后的In-Si键长约为2.72Å，In原子与周围Si原子的相互作用较强，吸附结构的稳定性也较高。H2位吸附结构的稳定性与H1位相当，但由于In原子与周围Si原子的配位情况不同，导致其电子结构和吸附能等性质存在一定的差异。从电子态密度的角度分析，H2位吸附时In原子的电子态与周围Si原子的电子态相互作用较强，在费米能级附近出现了新的电子态，这对吸附体系的电学性能产生了重要影响。在H3位（最近邻穴位）吸附时，In原子位于由三个Si原子围成的三角形中心，与周围三个Si原子的距离最短。优化后的In-Si键长约为2.73Å，In原子与周围Si原子的相互作用最强，吸附结构最为稳定。计算结果表明，In原子吸附在Si(111)-(2×2)表面的H3位时，体系的吸附能最低，为-2.56eV。这表明在该吸附位点，In原子与Si(111)表面原子之间的相互作用最强，形成的吸附结构最为稳定。从吸附能的定义可知，吸附能越负，体系越稳定，说明In原子在H3位吸附时，与Si表面原子形成了较强的化学键，降低了体系的总能量。通过对不同吸附位点的几何构型和吸附能的分析，我们可以清晰地看到In原子在Si(111)表面的吸附行为与吸附位点密切相关。H3位由于其独特的原子配位环境和较短的In-Si键长，使得In原子在该位点吸附时体系具有最低的能量和最高的稳定性。这种稳定性不仅影响着In原子在Si(111)表面的吸附状态，还对吸附体系的电子结构、电学性能等产生了深远的影响。例如，在研究吸附体系的电子态密度时发现，In原子在H3位吸附后，Si(111)表面的电子态密度分布发生了显著变化，在费米能级附近出现了新的表面态，这与In原子在该位点与Si表面原子的强相互作用密切相关。对于(√3×√3)R30°超原胞模型，In原子在Si(111)表面的吸附位点和吸附结构更为复杂。由于超原胞的对称性和原子排列方式的变化，In原子在该模型中的吸附位点和吸附结构与(2×2)超原胞模型存在一定的差异。在(√3×√3)R30°超原胞中，In原子除了可以吸附在类似于(2×2)超原胞中的顶位、桥位和穴位外，还可以吸附在一些特殊的位置，这些位置的原子环境和配位情况与传统的吸附位点不同。在对(√3×√3)R30°超原胞模型的吸附结构进行优化时，我们发现In原子在某些特殊位置吸附时，会导致表面原子的重构，形成新的吸附结构。这些重构后的吸附结构具有独特的几何构型和电子结构，对In原子在Si(111)表面的吸附稳定性和电学性能产生了重要影响。例如，在一种重构后的吸附结构中，In原子与周围Si原子形成了一种特殊的键合方式，使得In-Si键长和键角发生了变化，从而改变了吸附体系的能量和电子态分布。通过对这种重构吸附结构的电子态密度分析发现，In原子的吸附导致了表面电子态的重新分布，在费米能级附近出现了新的电子态，这些新的电子态与In原子的价电子和Si表面原子的电子相互作用密切相关，对吸附体系的电学性能产生了显著影响。确定In在Si(111)表面的吸附位点和分析不同吸附结构的几何构型是研究其吸附特性的基础。通过对(2×2)和(√3×√3)R30°超原胞模型的深入研究，我们明确了In原子在不同吸附位点的吸附结构和稳定性差异，为进一步研究In在Si(111)面的吸附特性提供了重要的依据。这些研究结果不仅有助于我们深入理解In原子与Si(111)表面原子之间的相互作用机制，还为基于In吸附的新型半导体器件的设计和开发提供了理论指导。4.2吸附能的计算与分析吸附能作为衡量原子在表面吸附稳定性的关键物理量，其计算对于深入理解In在Si(111)面的吸附特性至关重要。吸附能的计算公式为E_{ad}=E_{total}-E_{sub}-nE_{atom}，其中E_{ad}表示吸附能，E_{total}为吸附体系的总能量，E_{sub}为Si(111)衬底的能量，n为吸附原子的个数，E_{atom}为单个吸附原子的能量。通过该公式，我们能够准确计算不同吸附结构下In原子在Si(111)面的吸附能，进而分析吸附稳定性与吸附结构之间的内在联系。对于(2×2)超原胞模型，我们详细计算了In原子在不同吸附位点（顶位、桥位、H1位、H2位、H3位）的吸附能。计算结果表明，In原子吸附在顶位时，吸附能为-1.85eV；吸附在桥位时，吸附能为-2.10eV；吸附在H1位时，吸附能为-2.30eV；吸附在H2位时，吸附能为-2.45eV；而吸附在H3位时，吸附能最低，为-2.56eV。从这些数据可以清晰地看出，In原子在H3位的吸附能最负，这意味着In原子在该位点与Si(111)表面原子之间的相互作用最强，形成的吸附结构最为稳定。与其他吸附位点相比，H3位独特的原子配位环境使得In原子能够与周围三个Si原子形成较强的化学键，从而降低了体系的总能量，提高了吸附稳定性。在分析吸附能随覆盖度的变化规律时，我们通过在(2×2)超原胞中逐步增加In原子的数量，模拟了不同覆盖度下的吸附体系。当覆盖度为1/4ML（单层覆盖度的四分之一）时，每个In原子的吸附能为-2.56eV；随着覆盖度逐渐增加到1/2ML，吸附能变为-2.48eV；当覆盖度达到3/4ML时，吸附能为-2.40eV；而当覆盖度为1ML时，吸附能为-2.35eV。可以发现，随着覆盖度的增加，In原子的吸附能逐渐升高（绝对值减小），这表明吸附结构的稳定性逐渐降低。这种变化趋势主要是由于随着In原子数量的增加，In原子之间的相互排斥作用逐渐增强，抵消了部分In原子与Si表面原子之间的相互吸引作用，从而导致吸附能升高，吸附稳定性下降。对于(√3×√3)R30°超原胞模型，In原子在不同吸附位点的吸附能计算结果同样显示出明显的差异。在一些特殊的吸附位点，由于表面原子的重构和原子间相互作用的复杂性，吸附能的变化更为显著。在一种重构后的吸附结构中，In原子与周围Si原子形成了特殊的键合方式，其吸附能为-2.65eV，比(2×2)超原胞中H3位的吸附能更负，说明这种重构吸附结构具有更高的稳定性。这种特殊吸附结构的形成与(√3×√3)R30°超原胞的原子排列方式密切相关，其独特的对称性和原子间距使得In原子能够与表面原子形成更稳定的化学键，进一步降低了体系的能量。在分析吸附能随覆盖度的变化时，(√3×√3)R30°超原胞模型呈现出与(2×2)超原胞模型相似但又有所不同的规律。随着覆盖度的增加，吸附能总体上呈现升高的趋势，但在某些特定的覆盖度下，由于表面重构和原子间相互作用的变化，吸附能会出现局部的波动。当覆盖度在1/3ML附近时，由于表面原子的重构形成了一种相对稳定的结构，使得吸附能在该覆盖度下略有降低，表现出较好的吸附稳定性。这种局部波动的现象表明，在(√3×√3)R30°超原胞模型中，吸附能不仅受到In原子之间相互作用和In与Si表面原子相互作用的影响，还与表面重构过程中形成的特殊结构密切相关。通过对不同超原胞模型中In原子吸附能的计算与分析，我们明确了In在Si(111)面的最稳定吸附结构为(2×2)超原胞中In原子吸附在H3位以及(√3×√3)R30°超原胞中的特殊重构吸附结构。同时，深入了解了吸附能随覆盖度的变化规律，为进一步研究In在Si(111)面的吸附特性提供了重要的能量学依据。这些研究结果对于理解In原子在Si(111)表面的吸附行为和微观吸附机理具有重要意义，也为基于In吸附的新型半导体器件的设计和制备提供了关键的理论支持。4.3电子结构分析为了深入探究In在Si(111)面吸附后的电子特性，对吸附体系的电子态密度、电荷密度分布和差分电荷密度进行细致分析，这对于揭示In-Si之间的成键特性和电子转移情况至关重要。通过第一性原理计算，获取了In吸附前后Si(111)表面的电子态密度（DOS）图，结果如图1所示。在未吸附In原子时，Si(111)表面的电子态密度呈现出典型的半导体特征，价带顶位于-5.0eV到0eV之间，导带底位于1.0eV到5.0eV之间，存在明显的能隙，这与Si材料本身的半导体性质相符。当In原子吸附在Si(111)-(2×2)表面的H3位后，电子态密度发生了显著变化。在费米能级附近，出现了新的电子态，这些新的电子态主要来源于In原子的5s和5p电子与Si原子的3s和3p电子之间的相互作用。在-2.0eV到0eV的能量范围内，In的5p电子态与Si的3p电子态发生强烈杂化，形成了新的成键态，这表明In原子与Si表面原子之间形成了较强的化学键。同时，在费米能级附近，由于新电子态的出现，体系的态密度明显增加，这使得体系的导电性增强，表面表现出一定的金属性。这种在费米能级附近出现新电子态的现象，对吸附体系的电学性能产生了重要影响，为理解In吸附后Si(111)表面电学性质的改变提供了重要依据。[此处插入图1：Si(111)表面吸附In前后的电子态密度图]电荷密度分布能够直观地展示电子在原子周围的分布情况，对于理解原子间的相互作用和化学键的形成具有重要意义。通过绘制In吸附在Si(111)-(2×2)表面H3位后的电荷密度分布图（图2），可以清晰地看到In原子与周围Si原子之间的电荷分布情况。In原子周围的电荷密度明显高于Si原子，这表明In原子在吸附过程中吸引了部分电子，电子云向In原子周围聚集。在In-Si键的区域，电荷密度呈现出明显的峰值，这说明In原子与Si原子之间形成了较强的共价键，电子云在In-Si键之间发生了明显的重叠。从电荷密度的等值面图中可以看出，In原子与周围三个Si原子之间的电荷分布较为均匀，形成了稳定的化学键结构。这种电荷分布的特征与吸附能的计算结果相一致，进一步证实了In原子在H3位吸附时体系的稳定性。[此处插入图2：In吸附在Si(111)-(2×2)表面H3位后的电荷密度分布图]差分电荷密度则能够更直观地展示吸附前后电荷的转移情况，对于揭示In-Si之间的成键本质具有关键作用。通过计算In吸附在Si(111)-(2×2)表面H3位后的差分电荷密度（图3），可以清晰地看到电荷的转移方向和数量。差分电荷密度图中，黄色区域表示电荷聚集，青色区域表示电荷耗尽。在In原子与Si原子之间，出现了明显的黄色区域，这表明In原子与Si原子之间发生了电荷转移，Si原子的部分电子转移到了In原子周围。通过对差分电荷密度的积分计算，得到In原子与Si原子之间的电荷转移量约为0.35e，这表明In-Si之间形成了极性共价键。In原子由于电负性相对较小，在与Si原子的相互作用中失去了部分电子，而Si原子则获得了这部分电子，从而形成了具有一定极性的化学键。这种电荷转移和极性共价键的形成，对吸附体系的电子结构和化学性质产生了重要影响，进一步解释了In吸附后Si(111)表面电学性能和化学活性改变的原因。[此处插入图3：In吸附在Si(111)-(2×2)表面H3位后的差分电荷密度图]从电子态密度、电荷密度分布和差分电荷密度的分析结果可以看出，In原子在Si(111)表面的吸附导致了表面电子结构的显著变化。In原子与Si原子之间通过电子的转移和轨道的杂化，形成了较强的极性共价键，这种成键方式不仅改变了表面的电子态分布，还影响了体系的电学性能和化学活性。在实际应用中，这种电子结构的变化可以被利用来调控Si材料的电学性能，例如在半导体器件中，通过控制In原子的吸附量和吸附位置，可以实现对Si材料电导率和能带结构的精确调控，为制备高性能的半导体器件提供了新的途径。同时，In-Si之间的极性共价键也可能影响材料的化学稳定性和表面反应活性，在材料表面改性和催化等领域具有潜在的应用价值。4.4与实验结果及其他理论研究的对比验证为了进一步验证本研究中计算方法和模型的可靠性，将计算结果与已有的实验结果以及其他理论研究进行了详细对比。在吸附结构方面，本研究通过第一性原理计算确定In原子吸附在Si(111)-(2×2)表面的H3位时体系最为稳定，这与扫描隧道显微镜（STM）实验观察结果高度一致。STM实验直接观察到In原子在Si(111)-(2×2)表面的吸附位置，明确显示H3位是In原子的主要吸附位点，与本研究的计算结果相互印证，有力地支持了计算结果的准确性。同时，对于In原子在Si(111)表面形成的吸附结构，如In原子与周围Si原子的键长和键角等几何参数，本研究的计算结果也与相关实验测量值相符。计算得到的In-Si键长约为2.73Å，这与实验测量值相近，进一步验证了计算模型对吸附结构描述的可靠性。在吸附能的对比方面，本研究计算得到In原子在Si(111)-(2×2)表面H3位的吸附能为-2.56eV，与其他基于第一性原理的理论研究结果相比，具有较好的一致性。一些研究采用不同的交换关联泛函和计算方法，得到的In原子在Si(111)表面的吸附能在-2.45eV至-2.60eV之间，虽然存在一定的差异，但本研究的计算结果处于合理的误差范围内。这种差异主要源于不同计算方法和参数设置对电子相互作用描述的细微差别。例如，不同的交换关联泛函对电子之间的交换关联能的计算方式不同，可能导致吸附能的计算结果略有不同。然而，总体来说，本研究的计算结果与其他理论研究结果的一致性表明，所采用的计算方法和参数设置能够准确地描述In原子在Si(111)表面的吸附能。在电子结构方面，本研究通过计算得到的电子态密度和差分电荷密度结果也与实验和其他理论研究相契合。实验上，通过光电子能谱（PES）等技术测量了In吸附后Si(111)表面的电子结构变化，结果显示在费米能级附近出现了新的电子态，这与本研究计算得到的电子态密度图中在费米能级附近出现新电子态的结果一致。在研究In吸附在Si(111)-(2×2)表面H3位后的电子态密度时，计算发现In原子的5s和5p电子与Si原子的3s和3p电子相互作用，在费米能级附近形成了新的电子态，这与PES实验测量结果相符，进一步验证了计算结果的正确性。同时，本研究计算得到的差分电荷密度图中显示的In-Si之间的电荷转移情况，也与其他理论研究中关于In-Si成键特性的分析一致。通过对差分电荷密度的积分计算得到In原子与Si原子之间的电荷转移量约为0.35e，这与其他理论研究中通过不同方法计算得到的电荷转移量相近，表明本研究对In-Si之间电荷转移和化学键形成的描述是准确的。通过与实验结果及其他理论研究的详细对比验证，充分证明了本研究中采用的计算方法和模型的可靠性。这不仅为深入理解In在Si(111)面的吸附特性提供了坚实的理论基础，也为后续进一步研究In在Si(111)面的吸附行为和相关应用提供了有力的支持。五、Tl在Si(111)面的吸附特性研究5.1吸附位点与吸附结构在研究Tl在Si(111)面的吸附特性时，首先需要确定其在Si(111)表面可能的吸附位点，并对不同吸附结构的几何构型进行深入分析。与In在Si(111)面的吸附情况类似，Si(111)表面原子的六边形密堆积排列为Tl原子提供了多种潜在的吸附位点，主要包括顶位、桥位和穴位（H1位、H2位、H3位）。然而，由于Tl原子与In原子在原子半径、电负性等原子性质上存在差异，使得Tl原子在Si(111)表面的吸附行为和吸附结构与In原子有所不同。对于(2×2)超原胞模型，我们对Tl原子在各个吸附位点的吸附结构进行了细致的优化和分析。在顶位吸附结构中，Tl原子位于Si(111)表面的一个Si原子正上方，与下方Si原子形成垂直的键合关系。优化后的Tl-Si键长约为2.90Å，相较于In原子在顶位吸附时的In-Si键长更长。这是因为Tl原子的原子半径比In原子大，导致其与Si原子之间的距离增大。从原子间相互作用的角度来看，顶位吸附时Tl原子主要与下方的单个Si原子发生相互作用，与周围其他Si原子的相互作用相对较弱，使得吸附体系的稳定性相对较低。同时，由于Tl原子的电子云分布与In原子不同，顶位吸附时Tl原子对Si(111)表面电子云分布的影响也与In原子有所差异。在桥位吸附结构中，Tl原子位于两个相邻Si原子的中间位置。优化后的Tl-Si键长约为2.95Å，此时Tl原子与两个相邻Si原子的距离相等。桥位吸附结构的稳定性略高于顶位吸附结构，这是因为Tl原子通过与两个Si原子形成化学键，增加了吸附体系的结合力。然而，与In原子在桥位吸附时相比，Tl-Si键长更长，说明Tl原子与Si原子之间的相互作用相对较弱。从电子结构的角度分析，桥位吸附时Tl原子的电子云与两个Si原子的电子云发生重叠，形成了相对稳定的化学键，但由于Tl原子的电子云分布和原子性质的特点，这种化学键的强度相对较弱。在穴位吸附结构中，由于穴位的种类不同，Tl原子的吸附结构和稳定性也有所差异。在H1位（中心穴位）吸附时，Tl原子位于由六个Si原子围成的正六边形中心，与周围六个Si原子的距离相等。优化后的Tl-Si键长约为3.00Å，这种吸附结构使得Tl原子与周围Si原子的相互作用较为均衡，体系的稳定性相对较高。从空间结构的角度来看，H1位吸附时Tl原子处于表面原子形成的六边形空隙中心，周围Si原子对Tl原子形成了一定的包围作用，使得吸附体系具有较好的空间对称性。然而，与In原子在H1位吸附时相比，Tl-Si键长更长，表明Tl原子与Si原子之间的相互作用相对较弱。在H2位（次近邻穴位）吸附时，Tl原子位于由四个Si原子围成的四边形中心，与周围四个Si原子的距离相对较短。优化后的Tl-Si键长约为2.97Å，Tl原子与周围Si原子的相互作用较强，吸附结构的稳定性也较高。H2位吸附结构的稳定性与H1位相当，但由于Tl原子与周围Si原子的配位情况不同，导致其电子结构和吸附能等性质存在一定的差异。从电子态密度的角度分析，H2位吸附时Tl原子的电子态与周围Si原子的电子态相互作用较强，在费米能级附近出现了新的电子态，但与In原子在H2位吸附时相比，这些新电子态的分布和强度有所不同。在H3位（最近邻穴位）吸附时，Tl原子位于由三个Si原子围成的三角形中心，与周围三个Si原子的距离最短。优化后的Tl-Si键长约为2.98Å，Tl原子与周围Si原子的相互作用最强，吸附结构最为稳定。计算结果表明，Tl原子吸附在Si(111)-(2×2)表面的H3位时，体系的吸附能最低，为-2.20eV。这表明在该吸附位点，Tl原子与Si(111)表面原子之间的相互作用最强，形成的吸附结构最为稳定。然而，与In原子在H3位吸附时的吸附能（-2.56eV）相比，Tl原子的吸附能绝对值较小，说明Tl原子在Si(111)表面的吸附稳定性相对较低。通过对不同吸附位点的几何构型和吸附能的分析，可以看出Tl原子在Si(111)表面的吸附行为与吸附位点密切相关。H3位由于其独特的原子配位环境和较短的Tl-Si键长，使得Tl原子在该位点吸附时体系具有相对较高的稳定性。然而，与In原子在Si(111)表面的吸附相比，Tl原子的吸附稳定性总体较低，这与Tl原子的原子性质和电子结构密切相关。例如，Tl原子的电负性相对较小，使得其在与Si原子形成化学键时，电子云的偏移程度相对较小，化学键的强度相对较弱，从而导致吸附稳定性较低。对于(√3×√3)R30°超原胞模型，Tl原子在Si(111)表面的吸附位点和吸附结构更为复杂。由于超原胞的对称性和原子排列方式的变化，Tl原子在该模型中的吸附位点和吸附结构与(2×2)超原胞模型存在一定的差异。在(√3×√3)R30°超原胞中，Tl原子除