1 基本概念

1.1 模型定义

紧束缚模型是一类用于描述晶体中电子运动的近似理论。其基本出发点是:电子并非完全自由,而是主要受原子势场束缚,电子波函数在空间上以原子位置附近为中心,并可通过相邻原子间的跃迁形成延展的晶体态。该模型常被用来构建能带结构,并研究电子在固体中的分布、传播与相互作用特征。

1.2 核心假设

紧束缚模型通常包含几个关键假设。其一,电子在单个原子附近具有较强的局域性,原子轨道是合适的基底;其二,晶体中各原子之间存在有限的轨道重叠,电子可以在原子间发生跃迁;其三,主要相互作用可通过少数参数描述,如在位能和跃迁积分。由于这些假设,模型既保留了原子尺度信息,又能反映周期晶体中的能带形成。

1.3 适用物理图景

紧束缚模型适合描述原子间距较小、轨道重叠明确但电子仍带有局域特征的体系。它在共价晶体、某些金属、低维材料以及分子轨道问题中表现尤为有效。与强调连续自由传播的模型相比,紧束缚框架更适合从“原子如何拼接成晶体”这一角度理解材料电子结构。

1.3.1 电子局域化与离域化

该模型的物理图景介于局域化和离域化之间。电子虽然以原子附近为中心,但由于跃迁作用,它们可以在整个晶体中传播,形成具有一定带宽的能带。局域化强弱与跃迁强度密切相关:跃迁弱时电子更接近原子态,跃迁强时则更趋向离域态。

1.3.2 原子轨道重叠

原子轨道重叠是紧束缚模型的基础。相邻原子轨道之间的重叠越明显,跃迁概率通常越大,能带也更宽。不同轨道的对称性、方向性和空间分布会影响重叠形式,因此 s、p、d 等轨道在模型中往往需要分别处理。

1.4 与其他模型的关系

紧束缚模型与其他电子结构模型并非彼此孤立,而是处于同一理论谱系中。它常与自由电子模型、近自由电子模型、劈裂原子轨道模型等并列使用,用于从不同极限理解固体电子行为。

1.4.1 自由电子模型

自由电子模型把电子近似看作在均匀势场中运动的粒子,强调动能主导而忽略晶格的细致结构。相比之下,紧束缚模型强调原子势的束缚作用,更适合描述由局域轨道形成的能带。两者分别对应“电子近似自由”和“电子近似束缚”的不同极限。

1.4.2 劈裂原子轨道模型

劈裂原子轨道模型可视为紧束缚思想的直观表述:孤立原子的离散能级在晶体中因相互作用而分裂成一组密集能级,最终形成能带。该观点有助于理解为何原子数增多后,原先的原子能级会演化为连续或近连续的带状结构。

2 理论基础

2.1 晶体中的电子态

晶体电子态通常以周期势场中的波函数来描述。紧束缚模型不直接从连续波动方程的复杂解析解出发,而是借助原子轨道的叠加来构造近似本征态,从而将电子态表示为一组离散基函数的线性组合。

2.1.1 波函数展开

在紧束缚处理中,晶体波函数通常写成各原子轨道的加权和。系数反映电子在不同原子或不同轨道上的分布概率,而相位因子则体现周期性与晶格动量信息。这样的展开方式使问题从连续空间转化为离散矩阵问题。

2.1.2 Bloch 定理背景

晶体具有平移对称性,因此电子态满足 Bloch 定理,即波函数可写为平面波因子乘以周期函数。紧束缚模型与 Bloch 定理并不矛盾,实际上常将原子轨道按 Bloch 和的形式组合,以满足晶体平移对称性并构造能带本征态。

2.2 哈密顿量构造

紧束缚模型的核心是哈密顿量的参数化表示。哈密顿量由在位能、相邻原子间跃迁项以及可能的相互作用项组成,通过这些项可以近似刻画电子在晶格中的能量与运动。

2.2.1 在位能项

在位能项描述电子停留在某个原子轨道上时的能量,通常对应原子轨道本身的能级。不同元素、不同轨道的在位能往往不同,因此该项决定了能带的基准位置与多个轨道之间的相对能量顺序。

2.2.2 跃迁积分

跃迁积分用于描述电子从一个原子轨道跳到另一个轨道的幅度。它是紧束缚模型中最重要的参数之一,通常与原子间距离、轨道对称性以及重叠程度相关。跃迁积分越大,电子越容易在晶格中传播,能带也越容易展宽。

2.2.3 相互作用近似

在多数基础应用中,紧束缚模型只保留最主要的相互作用,例如最近邻跃迁或少量次近邻项。这样做可显著简化计算,同时保留主要物理特征。对于复杂体系,也可逐步加入自旋轨道耦合、晶场效应或电子关联修正

2.3 近邻近似

近邻近似是假设电子主要在有限范围内发生跃迁。它反映了原子轨道重叠随距离快速衰减的事实,因此能以较少参数获得较好的近似效果。

2.3.1 最近邻跃迁

最近邻跃迁只考虑相邻原子之间的电子跃迁,常用于一维链、蜂窝晶格等基础模型。该近似简洁且物理透明,是建立解析结果和直观图像的重要起点。

2.3.2 次近邻跃迁

次近邻跃迁考虑更远一层原子间的耦合。虽然其幅度通常小于最近邻项,但在某些材料中会显著影响能带非对称性、态密度以及费米能附近的精细结构,因此常用于提高模型精度

2.3.3 截断策略

截断策略是指在有限范围内保留跃迁项,其余项忽略。截断半径的选择取决于目标材料的电子结构特征与所需精度。通常,低维体系和局域轨道体系更适合采用较短范围的截断。

3 数学形式

3.1 线性组合原子轨道方法

线性组合原子轨道方法是紧束缚模型的常见数学表述。它把晶体电子态写成一组原子轨道的线性叠加,既便于保留局域物理,又能自然引入晶体周期性。

3.1.1 基函数选择

基函数通常选取与材料中主要价电子相关的原子轨道,如 s、p、d 轨道。基函数的选择直接影响模型的复杂度与适用范围。若体系以某类轨道主导,则采用对应的少轨道基底即可得到较好的近似。

3.1.2 系数求解

展开系数通过求解本征值问题获得。对于给定哈密顿量,系数决定不同原子轨道在某个能态中的权重与相对相位,从而给出能带形状、波函数分布和对称性信息。

3.2 能带色散关系

能带色散关系描述能量随晶体动量变化的规律,是紧束缚模型最常见的输出之一。不同晶格结构和跃迁范围会导致不同的色散形式。

3.2.1 一维链模型

在一维链中,最近邻紧束缚模型通常可导出简洁的色散关系,能量随波矢呈周期性变化。该模型常用于说明带宽、带边和有效质量等概念,是最经典的入门实例。

3.2.2 二维晶格模型

二维晶格中,晶格对称性更加丰富,色散关系会反映不同方向上的各向异性。例如方格晶格与蜂窝晶格会呈现不同的能带结构,进而影响态密度与输运行为。

3.2.3 三维晶格模型

三维体系的紧束缚色散更接近真实固体。由于原子排列更复杂,常需考虑多个轨道和多条跃迁路径。三维模型广泛用于半导体、金属和过渡金属材料的电子结构分析

3.3 本征值问题

紧束缚模型最终会归结为一个矩阵本征值问题。通过求解本征值和本征矢,可获得能带和电子态分布。

3.3.1 矩阵表示

在选定基底后,哈密顿量可写成矩阵形式。矩阵元素对应在位能和跃迁项,维数由基函数数量决定。对于多原子、多轨道体系,矩阵通常较大,但结构清晰,便于数值处理。

3.3.2 广义本征方程

当基函数之间并不正交时,需要求解广义本征方程,而不仅是普通本征方程。此时还要引入重叠矩阵,以正确反映不同原子轨道之间的非正交性。

3.3.3 数值求解方法

数值上常通过矩阵对角化、迭代算法或稀疏线性代数方法求解能带。对于大规模体系,可结合周期性简化或超胞技术,以获得可计算的近似结果。

4 典型模型

4.1 一维原子链

一维原子链是紧束缚理论中最基础的模型之一,能够直接展示能带如何由原子轨道拼接而成。

4.1.1 单原子链

单原子链通常假定每个晶格点只有一个等价原子轨道。该模型能解析给出能带色散,并清晰说明跃迁强度如何影响带宽。它是理解紧束缚思想的标准教材模型。

4.1.2 双原子链

双原子链包含两个不同或等价但非单原子重复的基元,因而能表现出更复杂的能带分裂。该模型常用于说明多个原子基元如何导致多个能带分支,并引入带隙或带交叉现象。

4.1.3 能隙形成

在双原子链等模型中,若晶格中存在不等价原子或不同在位能,原本可能连续连接的能级会在某些动量处打开能隙。能隙的出现通常与对称性破缺、轨道混合或周期势差异有关。

4.2 石墨烯紧束缚模型

石墨烯是紧束缚模型应用最具代表性的二维体系之一。其蜂窝晶格结构与 p 轨道主导的电子态,使得简单模型就能捕捉到重要的低能物理。

4.2.1 蜂窝晶格

石墨烯的晶格由两个互穿的子晶格构成。最简单的紧束缚模型只考虑最近邻 π 轨道跃迁,即可得到典型的锥形能带结构。这一模型成功解释了其特殊的电子色散特征。

4.2.2 Dirac 锥近似

在布里渊区特定点附近,石墨烯的能带可近似为线性色散,表现得类似无质量 Dirac 粒子。Dirac 锥近似使石墨烯在低能区呈现出与普通半导体不同的输运性质。

4.2.3 次近邻修正

加入次近邻跃迁后,石墨烯模型可更准确地描述电子-空穴对称性的偏离以及能带的细微修正。这类修正常用于提高理论与实验结果的一致性。

4.3 过渡金属体系

过渡金属材料的电子结构通常涉及多个 d 轨道,紧束缚模型在其中具有较高的实用价值。

4.3.1 d 轨道贡献

d 轨道空间延展性较强,且方向性明显,因此对能带结构、磁性和导电性影响显著。模型中若忽略 d 轨道,往往难以正确描述过渡金属的主要电子性质。

4.3.2 晶场分裂

在晶体环境中,d 轨道会因周围配位场而发生能级分裂。晶场分裂改变不同轨道的相对能量,进而影响电子占据顺序、磁矩大小与光谱特征。

4.3.3 多轨道耦合

过渡金属体系常需同时考虑多个轨道之间的耦合。多轨道紧束缚模型能够描述轨道混合、带交叉和局域磁性等现象,是研究复杂材料的重要工具。

5 物理性质分析

5.1 能带宽度

能带宽度是衡量电子在晶体中离域程度的重要指标,通常与跃迁强度直接相关。

5.1.1 跃迁强度影响

跃迁积分越大,电子更容易在晶格中传播,能带通常更宽;反之,跃迁较弱时能带较窄,电子更接近局域态。带宽变化会进一步影响导电性和有效质量。

5.1.2 轨道重叠程度

轨道重叠越强,带宽一般越大。由于不同轨道的方向性差异,s 轨道、p 轨道和 d 轨道对应的带宽特征常有明显区别,这也是材料电子性质差异的重要来源。

5.2 态密度

态密度描述单位能量区间内可供电子占据的态数。紧束缚模型通过能带结构可直接导出态密度,并分析其峰值与奇点。

5.2.1 van Hove 奇点

在某些能量处,态密度会出现尖峰或非解析行为,这类特征称为 van Hove 奇点。它们通常对应能带极值或鞍点,常对光学响应和电输运有显著影响。

5.2.2 维度效应

不同维度的体系其态密度形状差异明显。一般而言,一维体系更易出现尖锐奇点,二维体系常呈阶跃或峰结构,三维体系则通常更平滑。这一差异是低维材料物性的重要基础。

5.3 电学性质

紧束缚模型能够从能带占据和带隙结构出发,解释材料是导体、半导体还是表现出更复杂的绝缘行为。

5.3.1 导体

若费米能位于部分填充的能带中,电子可在外场作用下容易响应,材料通常表现为导体。紧束缚模型可通过能带交叠或带内占据来描述这类情况。

5.3.2 半导体

若价带与导带之间存在有限带隙,且费米能位于带隙内,材料通常表现为半导体。紧束缚模型常被用来分析带隙大小及其与晶格结构、杂化方式之间的关系。

5.3.3 金属-绝缘体行为

在某些条件下,体系的电学行为会在金属与绝缘体之间切换。紧束缚模型可以从能带填充、对称性破缺、局域化增强等角度提供解释,尽管更精确的描述往往需要额外相互作用修正。

5.4 电子局域化

电子局域化指电子波函数在空间上被限制于有限区域内。紧束缚模型本身强调局域轨道,因此特别适合分析局域化及其偏离情形。

5.4.1 Anderson 局域化

当体系存在随机无序时,电子波可能因多次散射而局域化,无法自由传播。这一现象称为 Anderson 局域化,常被用于研究无序介质中的输运抑制。

5.4.2 缺陷诱导态

晶体缺陷会在禁带或带边附近产生局域态,改变局部电子分布与电学响应。紧束缚模型能够较方便地将缺陷作为局部参数修正,从而分析缺陷态的出现与性质。

6 参数化与计算

6.1 参数来源

紧束缚模型的实用性很大程度上依赖参数的合理选取。常见参数来源包括实验拟合和第一性原理计算结果。

6.1.1 实验拟合

通过将模型结果与光谱、输运或其他实验数据匹配,可以反推出跃迁积分和在位能等参数。该方法直观有效,尤其适合已有较完整实验数据的体系。

6.1.2 第一性原理导出

也可从更基本的电子结构计算中提取紧束缚参数,用于构建有效模型。这种方式能够提高参数的物理一致性,并减少纯经验拟合带来的不确定性。

6.2 Slater-Koster 参数

Slater-Koster 方法为不同轨道间的跃迁积分提供了一套系统化表示,广泛用于多轨道紧束缚模型。

6.2.1 σ 键与 π 键

σ 键和 π 键对应不同的轨道重叠方式。前者通常沿键轴方向重叠较强,后者则与侧向重叠有关。两类键的区分有助于建立方向依赖的跃迁参数。

6.2.2 轨道方向余弦

轨道方向余弦用于将键方向与轨道对称性联系起来。通过方向余弦,可以把复杂的三维几何关系转化为标准化参数组合,便于在晶体中系统计算跃迁矩阵元。

6.2.3 参数表构建

参数表通常列出不同轨道对、不同键类型和不同距离对应的跃迁积分。构建规范的参数表有助于保持模型可复用性,也便于在不同材料间进行比较。

6.3 数值实现

紧束缚模型在实际应用中常依赖数值实现,以处理复杂晶格、大规模超胞或多轨道体系。

6.3.1 矩阵对角化

最常见的数值步骤是对哈密顿矩阵进行对角化,从而获得能量本征值和本征态。对于较小系统,可直接使用常规线性代数算法;对于大规模问题,则多采用稀疏矩阵技术。

6.3.2 周期边界条件

周期边界条件可模拟无限晶体,避免边界效应干扰。它使动量空间分析更加自然,也便于获得标准能带结构。

6.3.3 超胞方法

超胞方法通过把多个原胞组合成更大的计算单元,以描述缺陷、杂质、无序或畸变结构。该方法在研究真实材料中的局部扰动时非常常用。

7 扩展与改进

7.1 多轨道紧束缚模型

多轨道模型在单一轨道近似不足时尤为重要。它通过引入多个原子轨道,能够更完整地表示材料的电子结构。

7.1.1 s-p 模型

s-p 模型考虑 s 轨道与 p 轨道之间的耦合,常用于半导体和简单共价材料的基础描述。该模型在精度与简洁性之间具有较好平衡。

7.1.2 sp3 模型

sp3 模型进一步引入四个杂化轨道,适合描述具有强共价键特征的晶体结构。它在传统半导体和许多四面体配位材料中应用广泛。

7.1.3 sp3d5s* 模型

sp3d5s* 模型包含更多轨道自由度,能更准确地拟合复杂能带,尤其适用于需要同时考虑导带、价带和高能带的材料。

7.2 自旋相关扩展

在考虑电子自旋后,紧束缚模型可以描述更丰富的磁性和自旋输运现象。

7.2.1 自旋轨道耦合

自旋轨道耦合将电子自旋与轨道运动联系起来,常导致能带劈裂或拓扑相关效应。引入这一项后,模型可更好地处理重元素体系和低对称结构。

7.2.2 磁性有序

在磁性材料中,电子自旋可能呈现铁磁或反铁磁排列。紧束缚模型可通过自旋依赖的在位项或交换项刻画这种有序状态。

7.2.3 Zeeman 项

当外磁场存在时,Zeeman 项描述磁场与电子自旋或轨道磁矩的相互作用。它常用于研究磁场下的能级移动、简并劈裂与磁响应。

7.3 相互作用修正

基础紧束缚模型主要处理单电子图像,而实际材料中电子间相互作用也可能不可忽略,因此需要更高阶修正。

7.3.1 Hubbard 模型联系

Hubbard 模型可看作在紧束缚框架上加入局域电子-电子排斥的简化形式。它常用于研究磁性、关联效应和局域电子态的形成。

7.3.2 电子-电子关联

电子-电子关联会改变能带占据、准粒子寿命和激发谱结构。对于强关联体系,单纯的紧束缚近似往往不足,需要额外理论方法补充。

7.3.3 电子-声子耦合

电子与晶格振动的耦合会影响电阻、超导和局域畸变等现象。将电子-声子耦合纳入紧束缚模型后,可更深入研究温度效应与结构响应。

8 应用领域

8.1 材料电子结构

紧束缚模型是材料电子结构分析中的常用工具,尤其适合需要兼顾物理直观与计算效率的场景。

8.1.1 半导体能带分析

在半导体中,紧束缚模型可用于估算带隙、有效质量和能带边位置,并帮助理解掺杂、应变或晶格改变对电子结构的影响。

8.1.2 金属与合金

对于金属和合金体系,该模型可用于分析能带交叠、轨道杂化以及成分变化引起的电子态调整。其参数化形式也便于处理不同元素混合后的有效行为。

8.1.3 二维材料

二维材料中原子层数少、轨道耦合明确,紧束缚模型往往能以较低成本刻画其能带、边缘态和低维输运特征,因此应用非常广泛。

8.2 分子轨道理论

紧束缚思想不仅适用于周期晶体,也可推广到分子体系,用于分析共轭结构和分子轨道分布。

8.2.1 共轭体系

在共轭分子中,π 电子可沿多个原子离域。紧束缚模型常用于描述这类离域电子的轨道分裂与稳定化趋势。

8.2.2 芳香性分析

芳香性与分子中电子的环状离域有关。通过紧束缚模型可以定性判断环状共轭体系的能级填充与稳定性,为芳香性分析提供直观图像。

8.2.3 大分子近似

对于较大的有机分子或聚合物,紧束缚模型可作为有效近似,用较少参数捕捉主要轨道结构,避免完整量子化学计算带来的高成本。

8.3 器件模拟

紧束缚模型在纳米电子学中很常见,能够用于器件级别的输运和局域态分析。

8.3.1 纳米线

纳米线体系具有明显的一维或准一维特征,紧束缚模型可用于研究其子带结构、边界效应和尺寸量子化现象。

8.3.2 量子点

量子点中的电子被限制在有限区域,能级离散化明显。紧束缚框架适合描述局域势阱中的电子态、耦合与隧穿过程。

8.3.3 低维输运

在低维器件中,电子输运往往受几何约束、界面和局域散射影响。紧束缚模型能与输运理论结合,用于分析电导、隧穿与干涉效应。

9 局限性与评价

9.1 近似边界

紧束缚模型并非适用于所有体系,其准确性依赖于电子局域程度、轨道选择和参数质量。

9.1.1 强关联体系限制

对于电子强烈相互作用的材料,单电子紧束缚图像可能无法完整反映真实物理。此时需要引入更复杂的多体理论或关联修正。

9.1.2 强无序体系限制

当无序程度很高时,简单的有限跃迁近似可能不足以描述真实局域化与散射行为。尤其在强随机扰动下,模型参数本身也可能变得空间依赖且难以准确确定。

9.2 优点

尽管存在局限,紧束缚模型仍因其简洁和高效而在理论与计算中占据重要地位。

9.2.1 物理直观

该模型直接从原子轨道出发,便于理解能带从何而来,也容易把计算结果与化学键、轨道重叠等概念联系起来。

9.2.2 计算效率高

相较于更高精度的第一性原理方法,紧束缚模型通常参数更少、矩阵更稀疏,因此可在较低计算成本下处理较大体系。

9.2.3 可解释性强

由于模型参数具有明确物理意义,研究者可以通过调整在位能、跃迁积分等项,清楚地追踪各因素对能带与物性的影响。

9.3 发展现状

现代紧束缚研究已不再局限于传统解析模型,而是与更精确的计算方法和数据驱动技术结合得越来越紧密。

9.3.1 与第一性原理结合

当前常见做法是先通过第一性原理计算获得参考能带,再构建紧束缚有效模型进行简化描述。这样既保持物理精度,也提高了后续分析的灵活性。

9.3.2 与机器学习参数化结合

近年来,机器学习也被用于辅助提取或优化紧束缚参数。通过数据驱动方法,可以更快地建立适用于复杂材料体系的有效模型,并提升参数拟合的自动化程度。