1 基本概念
1.1 定义
平面波基组是由一组具有不同波矢的平面波函数构成的基函数集合，常用于展开波函数、电子态或其他满足波动描述的物理量。它的核心特点是形式统一、表达简洁，并且与周期性边界条件高度契合，因此在晶体物理和相关数值计算中应用广泛。
1.2 数学表达
平面波基组通常以指数函数的形式写出，其波动特征由波矢决定。对于周期性体系，波函数可以表示为若干个平面波的线性组合，从而将连续问题转化为离散的系数求解问题。
1.2.1 平面波函数形式
在三维空间中，平面波通常写作 \[ \phi_{\mathbf{G}}(\mathbf{r}) = e^{i(\mathbf{k}+\mathbf{G})\cdot\mathbf{r}} \] 其中，\(\mathbf{k}\) 为晶体动量，\(\mathbf{G}\) 为倒格矢。该形式沿空间各方向保持恒定振幅，仅相位随位置变化，是描述传播性波动的基本函数之一。
1.2.2 波矢与频率关系
对于满足波动方程的体系，波矢的大小与能量或频率之间通常存在直接联系。在线性、自由粒子或近似自由粒子的情形下，波矢越大，对应的空间振荡越快，所需描述的细节也越多，因此在数值展开中往往需要更高阶的平面波参与。
1.3 适用对象
平面波基组最适合处理具有规则周期结构的体系，但在适当处理后，也可用于某些非严格周期的对象。
1.3.1 周期性体系
晶体、层状材料、超晶胞和周期势场系统是平面波基组的典型应用对象。由于这些体系满足周期边界条件，平面波可以自然地在晶胞内展开，计算流程也较为统一。
1.3.2 准周期与有限体系的处理
对于分子、团簇或表面吸附等有限体系，通常会借助足够大的真空层构建超胞，使其在数值上呈现周期复制结构。这样虽然引入了人为周期性，但在适当设置单元尺寸后，可较好近似孤立体系的性质。
2 数学性质
2.1 正交性与完备性
平面波基组的重要优势之一在于其良好的正交结构与系统展开能力，这使其在数值线性代数和变分计算中十分方便。
2.1.1 正交归一关系
在周期盒或晶胞定义下，不同波矢对应的平面波之间通常满足正交关系。经过适当归一化后，基函数可形成正交归一组，从而简化矩阵元计算并提高数值稳定性。
2.1.2 完备展开
在无限多基函数的极限下，平面波集合可以构成完备基，从而任意足够平滑的周期函数都可被展开。这一性质为其系统性精度提升提供了理论基础：只需不断增加基函数数量，结果便可逐步逼近精确解。
2.2 周期边界条件
平面波与周期边界条件天然兼容，因此常作为周期体系的标准展开工具。
2.2.1 晶胞中的表示
在晶胞内部，波函数通常满足平移后的相位关系。平面波形式能够直接反映这一特征，使得晶体中的电子态可按布洛赫定理组织为“平面波乘以周期函数”的形式。
2.2.2 倒空间表述
平面波展开将问题从实空间转移到倒空间，体系中的空间变化由不同倒空间分量来表征。这样不仅便于处理周期结构，也适合借助傅里叶变换实现势能、密度和波函数之间的快速转换。
2.3 截断与收敛
实际计算中不可能使用无限多个平面波，因此需要在有限能量范围内进行截断。
2.3.1 截断能量
常见做法是设定平面波动能上限，只有满足条件的波矢分量才被纳入基组。这个上限通常称为截断能量，是控制计算精度与规模的关键参数。
2.3.2 基组收敛性
随着截断能量提高，基组规模增大，计算结果一般会逐渐收敛。对不同体系而言，收敛速度会有所差别；结构更平滑的势场通常更容易收敛，而具有强局域变化的对象则往往需要更多平面波。
3 构建方法
3.1 倒空间离散化
平面波基组的实际构建通常建立在倒空间离散化之上，即将连续波矢转化为一组离散的倒格空间点。
3.1.1 晶格与倒格矢
实空间晶格决定倒空间的几何结构，倒格矢则给出了允许的平面波分量。通过晶格参数可以生成对应的倒格空间，从而系统枚举可用的波矢集合。
3.1.2 k 点采样
在周期体系中，不同布里渊区波矢处的性质可能不同，因此通常需要对若干代表性 k 点进行采样。合理的 k 点网格可显著提升积分近似的准确度，是平面波计算中不可缺少的步骤。
3.2 基组截断策略
截断策略决定了哪些平面波被保留，以及如何控制计算成本。
3.2.1 球形截断
最常见的方案是以倒空间中的球面作为界限，仅保留动能低于阈值的分量。该方法具有几何上简单、各向同性较好的优点，适合多数三维体系。
3.2.2 方向性截断
在低维材料或各向异性很强的系统中，有时会采用方向性更强的截断方式，例如在某些方向上保留更多波矢分量，以更有效地描述结构特征和边界效应。
3.3 规范化处理
为了使不同基函数在数值上具有统一尺度，平面波需要进行规范化处理。
3.3.1 归一化系数
归一化系数用于确保基函数满足预定的积分条件，便于构造正交矩阵和求解本征值问题。不同程序或理论约定下，系数写法可能略有差异，但目的都在于保持一致的数值尺度。
3.3.2 相位约定
平面波的相位因子通常可有不同书写方式，例如正负号约定或指数前因子的差别。只要全体系内部保持一致，相位选择不会改变物理结果，但会影响公式书写与程序实现方式。
4 典型应用
4.1 固体物理计算
平面波基组是固体电子结构计算中的基础工具之一，尤其适用于晶体、半导体和金属体系。
4.1.1 能带结构
通过平面波展开并求解本征值问题，可以获得电子能带随 k 点变化的关系。能带结构能够反映导体、绝缘体和半导体的差异，也常用于分析有效质量、带隙和电子输运性质。
4.1.2 电子态密度
在获得大量波矢点和能级后，可以进一步统计电子态密度。态密度有助于理解费米能级附近的电子分布，并为热学、电学与光学性质分析提供基础。
4.2 量子化学模拟
在量子化学中，平面波基组并不总是首选，但在某些场景下具有独特优势。
4.2.1 分子与团簇的周期近似
对于分子或小团簇，可将其放入较大的周期超胞中进行模拟。这种方法便于使用统一的周期性框架，但需要足够大的真空区以减少镜像相互作用。
4.2.2 与赝势结合的计算
平面波方法常与赝势结合，以减少对核附近快速振荡波函数的直接描述需求。这样既能降低基组规模，也能提升对较重元素体系的计算可行性。
4.3 波动方程求解
除电子结构外，平面波展开也可用于多种波动型方程的数值求解。
4.3.1 薛定谔方程
在量子力学中，平面波可作为求解薛定谔方程的展开基底。通过将波函数写成平面波线性组合，可以将微分方程转化为矩阵本征值问题，便于数值处理。
4.3.2 泊松方程与场展开
电势、静电场和某些标量场也可在平面波基下展开。对于周期边界条件下的泊松方程，倒空间求解往往比实空间离散更高效，尤其适合大规模周期系统。
5 计算实现
5.1 快速傅里叶变换
平面波方法的高效实现与快速傅里叶变换密切相关。
5.1.1 实空间与倒空间转换
许多物理量在实空间中较容易处理，而平面波系数则在倒空间中更简洁。快速傅里叶变换可在两者之间高效转换，从而避免直接计算大量积分。
5.1.2 计算复杂度优化
借助快速算法，传统上高成本的卷积和变换过程可显著加速。这也是平面波方法在大规模计算中仍具竞争力的重要原因之一。
5.2 存储与并行化
当基组规模增大时，数据组织和并行策略对性能影响很大。
5.2.1 数据布局
平面波系数通常需要按波矢、k 点、能带或自旋等维度组织。合理的数据布局有助于提升缓存命中率，减少通信开销，并改善整体效率。
5.2.2 多线程与分布式计算
现代平面波程序一般支持多线程和分布式并行。前者适合单节点内加速，后者则用于跨节点处理大规模体系，两者配合可显著扩大可计算问题的尺寸。
5.3 常见软件中的实现
平面波基组已成为多类电子结构程序的标准配置之一。
5.3.1 平面波程序框架
典型程序通常包含结构输入、基组生成、哈密顿量构造、自洽迭代和结果后处理等模块。不同软件在细节上各有差别，但整体流程较为相似。
5.3.2 输入参数与控制项
常见控制项包括截断能量、k 点网格、赝势类型、收敛阈值以及并行设置等。这些参数直接影响计算精度、速度和稳定性，需要根据体系特征加以调整。
6 优缺点
6.1 优势
平面波基组因其规范统一和理论清晰而广受欢迎。
6.1.1 系统性强
通过提高截断能量即可逐步改善精度，且收敛行为通常较容易监控。这种可控性使其在高精度计算中很有吸引力。
6.1.2 适合周期体系
由于天然满足周期边界条件，平面波在晶体和其他周期结构中具有明显优势，公式处理简洁，数值实现也较成熟。
6.2 局限性
尽管优点突出，平面波基组也存在一些固有不足。
6.2.1 对局域特征描述效率较低
当体系中存在强局域电子态、尖锐势阱或核附近快速变化时，平面波需要大量高频分量才能准确表示，因此效率可能不如局域基组。
6.2.2 计算成本较高
与某些更紧凑的基组相比，平面波方法往往需要更大的存储空间和更高的运算量，特别是在大体系或高截断能量下，这一问题更为明显。
6.3 与其他基组比较
平面波常与其他数值基组方法对照使用，以便选择更合适的计算方案。
6.3.1 高斯型基组
高斯型基组在分子量子化学中非常常见，优点是对局域轨道表达高效，且积分计算有成熟算法。相比之下，平面波更适合周期体系，但在孤立分子问题中有时不如高斯基紧凑。
6.3.2 实空间网格方法
实空间网格方法直接在空间点上离散方程，直观且易处理局域边界条件。平面波方法则在周期问题上更自然，且倒空间处理常带来较高效率。两者在不同问题类型中各有优势。
7 相关扩展
7.1 赝势与超软赝势
赝势技术常与平面波基组配套使用，用于简化原子核附近的复杂行为。
7.1.1 核电子处理
核电子通常分布在原子核附近，变化剧烈且对化学性质影响相对间接。赝势方法通过用等效势替代真实核势，减少对这些电子的显式求解需求。
7.1.2 价电子近似
在多数材料和化学问题中，真正决定键合与响应性质的是价电子。赝势近似保留价电子的主要物理效应，同时显著降低计算难度。
7.2 增广平面波方法
为了弥补纯平面波在核附近描述不足的问题，发展出多种增广方案。
7.2.1 LAPW 方法
线性增广平面波方法在原子球内部使用更适合局域行为的函数，在球外仍采用平面波展开。它兼顾了局域精细结构与周期体系处理能力。
7.2.2 PAW 方法
投影缀加波方法通过构造变换，将平面波的平滑表示与真实波函数的局域信息结合起来。该方法在精度和效率之间取得了较好平衡，应用十分广泛。
7.3 多尺度与混合基组
为了适应更复杂的问题，平面波常与其他基组形成混合策略。
7.3.1 平面波-局域基组合
在某些计算中，周期部分用平面波处理，而局域强相关区域则用局域基函数补充。这种组合可在保持整体便利性的同时，提高对特殊区域的描述能力。
7.3.2 自适应基组策略
自适应方案会根据体系局部特征动态调整基函数分布，使资源更多投入到变化剧烈的区域。此类方法常用于追求更高效率或更强灵活性的数值框架。