通过第一性原理方法对等离子体纳米粒子的光谱进行建模需要耗费大量的计算资源,因此需要具有高准确度/计算成本比的方法。本文,我们表明,如果在辅助基组中每个原子仅采用一个 s 型函数,并采用适当优化的指数,则可以大大简化时间相关密度泛函理论 (TDDFT) 方法。这种方法(称为 TDDFT-as,代表辅助 s 型)可以预测不同尺寸和形状的银纳米粒子的激发能量,与参考 TDDFT 计算相比,平均误差仅为 12 meV。TDDFT-as 方法类似于线性响应处理的紧束缚近似方案,但适用于原子跃迁电荷,这里精确计算(即没有来自群体分析的近似)。我们发现,原子跃迁电荷的精确计算大大改善了宽能量范围内的吸收光谱。
材料的基本性质由原子核势能、电子质量和相互排斥力下的电子决定。不同材料之间的变量是离子势。计算电子性质的逻辑程序是从势到电子分布。这使得从原子、分子到固体的材料性质的实际计算成为可能。由于许多人的努力,这种方法已经蓬勃发展。这个概念类似于将人类人口分布的预测从丘陵和山谷的景观转变为从人口分布确定景观。在原子系统中,量子的怪癖允许这种切换,但决定它只是量子态断层扫描中的一个切片。作者分享了他从这个切片开发的经验,但接近用人口切换景观的强大概念。
稿件的所有计算资源均可在 Git 存储库 [1] 和相关数据文件 [2] 中找到。其中包括用于生成输入文件、运行计算、处理和分析数据以及生成图形的脚本。文件组织在存储库中的 README.md 文件中描述。所有 DFT 计算均使用 FHI-aims [3] 完成,其使用原子中心基组和数值径向部分。我们使用严格的默认基组和网格设置,这可确保本文研究的范德华 (vdW) 体系的结合能数值收敛到 0.1 kcal/mol。MBD 计算借助于集成到 FHI-aims 中的 Libmbd 库 [4] 执行,并且可使用当前开发版本在 FHI-aims 中直接执行 MBD-NL 计算。我们目前的实现不包括函数导数δαVV'[n]/δn,因此本文在自洽PBE密度上评估MBD-NL,而导数的实现仍在进行中。重要的是,已发现由范德华相互作用引起的电子密度变化对相互作用能和核力的影响可忽略不计[5]。S66、X23和S12L集的PBE、PBE0和VV10能量取自[6],其使用与本文相同的数值设置。对于分子晶体,所有DFT和MBD计算均使用逆空间中密度至少为0.8˚A的k点网格。对于硬固体,我们使用了[7]中的k点密度。所有分子和晶体几何形状均直接取自各自的基准集,未进行任何松弛。表 I 报告了 MBD-NL、MBD@rsSCS 和 VV10 与 PBE 和 PBE0 函数结合对一组有机分子晶体(X23,[11])、一组超分子复合物(S12L,[12])和一组 26 种层状材料(称为“26”,[10])的性能。在标准范德华数据集中,S12L 是唯一一个 MBD-NL 与 PBE 和 PBE0 函数结合时实现不同性能的数据集。这主要是因为 PBE 与 π – π 复合物的结合力比 PBE0 略强。对于大型 π – π 复合物,半局部 DFT 和长程范德华模型之间的适当平衡尚不清楚 [6]。在“26”集中,MBD@rsSCS 哈密顿量对 26 种化合物中的 20 种具有负特征值。然而,为了获得有限能量,我们使用了 Gould 等人提出的特征值重标度 [9]。图 1 比较了由 PBE-NL 计算的混合有机/无机界面的结合能曲线以及 Ruiz 等人的 MBD@rsSCS 和 TS 方法的表面变体 [13]。表 II 列出了 DFT+MBD 的时序示例