研究目的
研究在赝势-密度泛函理论中采用傅里叶变换计算部分芯电荷密度时引入的数值误差,并提出避免这些误差的解决方案。
研究成果
该研究表明,采用傅里叶变换在实空间网格中计算部分芯电荷密度的常用方法,可能导致交换关联势出现高达1电子伏特的显著数值误差。问题源于部分芯电荷在倒易空间中的缓衰减特性。研究建议在截断半径rcore的球体内使用更平滑的部分芯电荷函数,或直接在实空间网格中计算部分芯电荷ρcore partial以避免此问题。
研究不足
使用傅里叶变换和不足的平面波截断所导致的问题仅出现在包含大真空区域的体系计算中(如低维材料和体材料表面),并影响那些在真空区域具有较大成分的未占据态能量。误差大小因体系不同而存在很大差异。
1:实验设计与方法选择:
本研究采用基于平面波的量子 espresso、PARATEC 代码和基于实空间的 PARSEC 代码进行计算。所有计算均在交换关联泛函的局域密度近似(LDA)下进行,使用 Perdew?Zunger 参数化方法,并采用 Troullier?Martin 方法生成的模守恒赝势。
2:样本选择与数据来源:
研究聚焦于通过超胞模拟的低维材料或表面,具体包括单层 MoSe? 和铜 [100] 表面。
3:实验设备与材料清单:
研究使用了美国国家能源研究科学计算中心(NERSC)和德克萨斯高级计算中心(TACC)提供的计算资源。
4:实验步骤与操作流程:
研究比较了通过傅里叶变换和实空间直接计算部分芯电荷密度所得的结果,分析了这些方法对交换关联势和能带结构准确性的影响。
5:数据分析方法:
研究分析了交换关联势中的数值误差及其对材料能带结构的影响,并检验了计算相对于平面波截断值的收敛性。
独家科研数据包,助您复现前沿成果,加速创新突破
获取完整内容