过渡金属用于二极管的第一性原理模拟：密度泛函理论计算中不准确性的挑战与克服方法

《Nanoscale Advances》：First principles simulations of transition metals for diodes: challenges and approaches to overcome inaccuracies in calculations

　　

研究人员采用密度泛函理论（Density Functional Theory，DFT）研究了4B族（Cr、Mo、W）、8B族（Ni、Pd、Pt）和9B族（Cu、Ag、Au）过渡金属的电子、光学及力学性质。使用三种交换关联（exchange–correlation，xc）泛函——广义梯度近似Perdew–Burke–Ernzerhof（PBE）、meta广义梯度近似R2SCAN及杂化泛函Heyd–Scuseria–Ernzerhof（HSE06）——进行比较，以评估其对不同物理性质预测能力的差异。态密度（Density of States，DOS）计算结果表明PBE与R2SCAN吻合较好，而HSE06因引入精确Hartree–Fock交换导致价带与导带区域出现显著偏移。光学性质分析显示，对于d带填满或近填满的金属，PBE和R2SCAN能更高精度地复现实验曲线；而对于具更局域化d电子态的金属，HSE06表现更优。力学性质以杨氏模量（Young's modulus）表征，结果显示未发s→d电子激发（electronic d promotion）的金属具有较高弹性模量。本工作阐明了电子结构预测、键杂化（bond hybridization）及电子激发如何决定所考察金属元素的性质，指出了三种计算方法各自的优缺点，并为纯金属体系的可靠模拟提供指导原则。

论文解读：过渡金属用于二极管的第一性原理模拟——DFT计算不准确性的挑战与克服方法

该文由Marina Petrossi、Elaheh Mohebbi、Eleonora Pavoni、Pierluigi Stipa、Davide Mencarelli、Emiliano Laudadio及Luca Pierantoni完成，发表于《Nanoscale Advances》。

研究背景与目的

原子尺度模拟是材料科学研究金属元素本征性质的重要工具，其中密度泛函理论（Density Functional Theory，DFT）被广泛应用。然而DFT中交换关联（exchange–correlation，xc）泛函存在自相互作用误差（self-interaction error）与离域化误差（delocalization error），不同泛函对不同过渡金属的适用性尚无统一共识。过渡金属因可能存在ns→(n?1)d电子激发（electronic promotion），电子组态复杂，加之文献中计算方法可靠性常相互矛盾，纯金属状态下其电子、光学及力学性质的系统性对比研究较为缺乏。二极管等器件中金属电极/接触层的性能与金属的功函数、光学常数及力学稳定性密切相关，因此需在器件制备前通过可靠模拟筛选最优金属构型。本研究旨在通过对九种典型过渡金属（4B族Cr、Mo、W；8B族Ni、Pd、Pt；9B族Cu、Ag、Au）分别采用PBE、R²SCAN及HSE06三种xc泛函开展DFT计算，并结合分子动力学（Molecular Dynamics，MD）模拟力学性质，系统评估各泛函对不同金属不同性质预测的优劣，为二极管用金属的第一性原理模拟提供泛函选择依据。

主要关键技术方法

研究人员选取周期表中4B族（Cr、Mo、W，体心立方bcc，空间群229）、8B族（Ni、Pd、Pt，面心立方fcc，空间群225）及9B族（Cu、Ag、Au，fcc，空间群225）共九种纯过渡金属为研究对象。电子、光学及态密度（Density of States，DOS）计算均采用Quantum ATK软件包进行DFT模拟，分别使用GGA-PBE、meta-GGA R²SCAN及杂化HSE06泛函，配以Norm-conserving PseudoDojo赝势（Pseudopotential，PP）及LCAO Ultra基组，截断能150 hartree，k点密度20 ?^?3（光学计算为28 ?^?3），考虑各金属特定的等离子体频率与逆寿命时间以纳入带内贡献。力学性质通过GROMACS 2023.3进行经典MD模拟，采用修正CHARMM27力场描述金属–金属Lennard-Jones相互作用，构建10×10×100 nm³盒子并施加周期性边界条件，经能量最小化和298 K平衡后运行100 ns生产步，沿z方向施加单轴拉/压应变由应力–应变曲线初始斜率求取面内杨氏模量（Young's modulus）。

研究结果

Electronic properties（电子性质）

通过对比三种泛函得到的DOS曲线发现：4B族Cr、Mo、W中PBE与R²SCAN峰形及位置高度相似，HSE06使价带相对费米能级下移且导带态密度降低，反映精确交换对激发态的描述作用；8B族Ni、Pd、Pt的PBE与R²SCAN结果几乎重合，HSE06在费米能级以下一致性好但在6 eV以上出现明显峰位偏移及强度变化；9B族Cu、Ag、Au同样PBE与R²SCAN接近，HSE06引起更显著的费米能级远离及d带偏移。结论：PBE与R²SCAN对金属DOS描述相近，R²SCAN引入的动能密度项未显著改变金属元素描述；HSE06因含精确HF交换使d电子态发生能移，对巡游（itinerant）d电子金属易过移d子带，对局域化d电子更有利。

Optical properties（光学性质）

在1–4 eV能量范围内对比折射率n与消光系数k与实验值：4B族金属中Cr的n在1.3–2.3 eV内HSE06较接近实验值，但k全体方法均高估；Mo三种泛函n略有高估但2 eV后吻合，k R²SCAN稍优；W的n从约1.4 eV起R²SCAN最佳，k全体略高估低能区。8B族Ni、Pd、Pt的n与k以PBE和R²SCAN与实验符合最好，HSE06在低能区n略高估。9B族Cu的n以HSE06较接近实验，k以PBE/R²SCAN为佳；Ag的n与k HSE06与实验几乎完全一致；Au的n HSE06佳，k在>2.5 eV低估而PBE/R²SCAN更好。结论：d带填满或近填满具巡游d电子的金属（8B、9B中Cu、Ag、Au的光学低频区）PBE与R²SCAN再现实验精度更高；具更局域化d电子态者（如Cr的部分光学区、Ag全波段、Au的n）HSE06表现更优。物理机制为d带越满越金属性、电子越巡游、屏蔽越强，合适精确交换比例应很小，半局域泛函更适宜；反之局域化d态受自相互作用误差影响大，HSE06可部分纠正此误差。

Mechanical properties（力学性质）

通过MD计算杨氏模量并与实验值比对：bcc结构Cr（252 GPa）、Mo（315 GPa）、W（360 GPa）高于fcc结构金属；fcc中Ni（204 GPa，无s→d激发）＞Pt（165 GPa，部分激发）＞Pd（119 GPa，双s电子激发填d带），Cu（120 GPa，激发填d带）＞Au（80 GPa）＞Ag（72 GPa）。未发生ns→(n?1)d电子激发的金属保持较高杨氏模量，电子激发削弱化学键强致模量下降；bcc相较fcc具更高硬度但与弹性模量无简单键长相关性。结论：晶体相与键杂化（特别是有无电子激发）是决定过渡金属杨氏模量的关键因素，未激发金属弹性更强。

总结与结论（Conclusions）翻译浓缩

本研究用PBE、R²SCAN及HSE06三种xc泛函考察了4B、8B、9B族九种过渡金属的电子、光学及力学性质。DOS显示PBE与R²SCAN给出相似轮廓，R²SCAN引入动能密度项未显著改变金属描述，HSE06因含精确交换引起d带偏移。partially filled d带金属（Cr、Mo、W）对泛函敏感；8B族在费米能级下三者吻合好；9B族随价电子增多差异增大。光学谱表明对具巡游d电子金属PBE与R²SCAN更一致，对某些具局域化d态金属HSE06更接近实验。力学结果确认bcc金属较fcc具更高杨氏模量，未发s→d电子激发的金属弹性模量更高。本研究为二极管设计用纯金属体系DFT模拟的泛函选择提供指导，后续将考察功函数、电阻率和热稳定性等参数。