第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

引言

第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
在上一章《第八章:杂化泛函!| 2026新版VASP基础教程》中,华算科技朱老师详细介绍了密度泛函理论基础泛函与近似。VASP通过赝势与平面波基组方法,运用密度泛函理论求解材料性质。本章将正式引入这一理论核心,从物理图像出发,赝势与基组及其在描述不同材料体系时的关键作用。
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

PAW赝势原理与背景

第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
赝势(Pseudopotential)是电子结构计算中为了解决原子核附近电子波函数振荡剧烈、平面波基组难以收敛的问题而提出的近似方法。传统赝势方法(Pseudopotential Approach)虽然极大简化了计算,但由于在计算中不可避免地丢失了核心区域电子的完整信息(如奇异点),导致其在处理第一族元素、过渡金属及稀土元素时精度受限。
为了弥补这一缺陷,投影增强波(Projector Augmented Wave,PAW)方法应运而生。PAW不仅保留了赝势方法的计算效率优势,还通过数学上的巧妙变换,完全恢复了全电子波函数(All-Electron Wavefunction)的精度。它可以被看作是赝势方法与全电子原子轨道方法(All-Electron Methods)的有机结合。
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

物理与数学原理

第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
PAW方法的核心思想是可逆的线性变换。它引入了一个辅助的平滑波函数(通常记作第九章:VASP中的赝势与基组!| 2026新版VASP基础教程),并通过一个称为“增强算符”(augmentation operator)的变换,将其映射为真实的全电子波函数Ψ
具体而言,PAW方法将空间划分为两个区域:
增强球体(Augmentation Sphere):围绕每个原子核的球形区域,记作ΩII为原子编号)。在此区域内,电子波函数波动剧烈,无法直接使用平面波描述。
间隙区域(Interstitial Region):球体之外的空间,电子行为更为平滑。
在间隙区域,由于波函数足够平滑,PAW直接使用平面波展开:
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
而在增强球体内部,PAW引入了局部原子基函数(Partial Waves)第九章:VASP中的赝势与基组!| 2026新版VASP基础教程和其对应的赝基函数第九章:VASP中的赝势与基组!| 2026新版VASP基础教程,以及用于投影的投影函数(Projectors)第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
这种划分利用了化学键的局域性:键合主要发生在球体之间的间隙区域(这里平面波足够),而球体内部主要涉及核与核心电子(这里需要局部化的原子轨道来描述)。
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
3. 核心方程与变换
PAW方法的数学描述如下:
全电子波函数的分解:
全电子波函数Ψ可以表示为辅助波函数第九章:VASP中的赝势与基组!| 2026新版VASP基础教程与在每个增强球体内的局部校正之和:
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
这表明Ψ第九章:VASP中的赝势与基组!| 2026新版VASP基础教程在球体外完全相同,仅在球体内部存在差异(通过局部基函数校正)。
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

电子密度的构造

第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
由于密度是能量泛函的关键变量,PAW必须精确构造总电子密度n(r)。其表达式为:
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
其中:
ñ(r)是由辅助波函数第九章:VASP中的赝势与基组!| 2026新版VASP基础教程直接计算得到的平滑密度。
nI(r)ñI(r)分别是全电子密度和赝密度在球体内的局部贡献,后者用于校正ñ在球体内部的误差。
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

与传统赝势的关系

第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
PAW方法实际上蕴含了传统赝势方法。如果在PAW的能量泛函中,对密度n1和 ñ1进行泰勒展开并截断到一阶,得到的近似形式正是传统赝势的表达式。这意味着:
PAW是更一般的理论:在特定近似下退化为传统赝势。
优越性:PAW通过保留了全电子波函数的完整信息(虽然在计算中只显式处理平滑波函数),避免了赝势中可能出现的“幽灵态”(Ghost State)问题,特别适合处理过渡金属、稀土元素和高压体系。
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

实际实现与应用

第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
在实际的第一性原理计算软件(如VASP、Quantum ESPRESSO、ABINIT)中,PAW方法已成为默认的高精度计算方案。其实现通常包括以下步骤:
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

生成PAW库

预先计算全电子原子结构,生成局部基函数、赝基函数和投影函数。
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

设置计算参数

用户选择平面波能量截断(ENCUT),该参数决定了辅助波函数的精度。
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

计算与收敛

由于第九章:VASP中的赝势与基组!| 2026新版VASP基础教程平滑且仅在球体外准确,PAW通常比传统赝势收敛更快,但在极端体系下仍然保留了全电子精度。
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

本章要点总结

第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

PAW赝势原理知识点梳理

物理与数学原理:介绍了投影函数、局部原子基函数
核心方程与变换:介绍了全电子波函数与电子密度
与传统赝势的关系:PAW赝势具有全电子波函数信息
第九章:VASP中的赝势与基组!| 2026新版VASP基础教程

下一步学习建议

下一章将正式引入本次教程的核心—KS方程。我们将从理论背景与意义、数学形式、求解流程,详细介绍KS方程,以及它在处理不同材料体系时的独特优势的应用,敬请期待!

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
(0)
上一篇 1天前
下一篇 7小时前

相关推荐