《Journal of Chemical Theory and Computation》：User-Defined Electrostatic Potentials in DFT Supercell Calculations: Implementation and Application to Electrified Interfaces
编辑推荐：
本文介绍了一种基于维也纳从头算模拟软件包（VASP）-Python接口的创新性方法，用于在密度泛函理论（DFT）超晶格计算中灵活地施加用户自定义的静电势。该方法解决了在模拟电化学过程、界面现象及外偏压作用下材料行为时，传统DFT代码应用外部电场不便的难题。通过在标准VASP软件内部直接操控静电势，并辅以必需的能量和力修正，本工作实现了对任意外加电场的精确控制。研究通过多个案例（包括带电表面分子吸附、场离子显微镜、电化学固-水界面及隐式溶剂模型等）展示了该方法的强大应用潜力，为从原子尺度深入理解电场调控下的物理与化学过程提供了高度灵活且稳健的计算框架。
在模拟分子和材料对外部静电势和电场的响应方面，密度泛函理论（DFT）凭借其在计算成本与精度间的良好平衡，已成为不可或缺的工具。然而，在实际DFT计算中施加用户定义的静电势并非易事，通常需要直接修改特定DFT程序的源代码。这不仅要求用户深入理解庞大复杂的代码库，还容易引入错误，且难以在不同软件版本间维护和复现。本研究提出并详细阐述了一种在DFT超晶格计算中实现任意静电势控制的方法，其核心在于利用最新发布的VASP 6.5.0版本中的Python接口，该接口允许通过Python脚本在源代码内部直接修改局域势等关键物理量。
理论计算细节
VASP-Python接口的工作流程
该Python接口为DFT计算提供了前所未有的灵活性。它能够在计算初始化阶段以及能量最小化或分子动力学（AIMD）运行的离子弛豫循环中，动态访问并修改局域势、核电荷、力及总能量等关键量。为实现外部静电势Vext的施加，只需在电子循环中将Vext添加到总势中。然而，由于接口设计的特性，Vext与原子核之间的相互作用并未被自动计入。因此，为了获得正确的能量和力，必须通过Python插件额外加入相应的修正项：对第I个原子核的能量修正ΔEI= -ZIVext(RI)和力修正ΔFI= -ZI?ext(RI)，其中ZI是核电荷，?ext是对应于Vext的外加电场。研究以真空中的中性氢原子为模型系统进行了基准测试，结果清晰地验证了未修正的能量和力会呈现与理论预测一致的线性偏差，而应用上述修正后，系统的物理行为得以正确恢复。
模拟带电表面与电势控制
在电化学等相关领域，模拟带电表面和界面是DFT计算的重要应用。这通常通过设置一个补偿带电模型表面电荷的对电极来实现，从而在超晶格系统中引入电场。本文介绍了两种对电极方案：氖计算对电极（Ne-CCE）和电荷密度对电极（CDCE）。Ne-CCE通过调节氖原子的核核心电荷来精确控制对电极电荷。而CDCE则更进一步，它不采用显式原子，而是引入一个具有高斯分布的外部对电极电荷密度ρext，并通过求解泊松方程计算其对应的静电势Vext，再通过Python接口加入计算。CDCE方法的优势在于突破了Ne能带对最大电场强度的限制，拓宽了可模拟的电场范围。
为了在AIMD模拟中更真实地再现恒电位实验条件，研究采用了热电位计技术，并通过VASP-Python接口重新实现了该方法。该技术的核心是在电极电荷nelectrode中引入涨落项，以在小尺度DFT元胞内模拟宏观恒电位动力学。在每个时间步，根据目标电位Φ0和瞬时电位Φ(t)的差异，结合系统本征电容和热涨落，动态更新电极电荷，从而在平均意义上将电极电位维持在设定值。
VASP中的实现
在实际操作层面，用户需要创建一个名为vasp_plugin.py的Python输入文件，其中定义了三个核心回调函数：local_potential、force_and_stress和occupancies。在local_potential函数中，将计算好的Vext添加到VASP的内部静电势上，并同时加入核电荷的能量修正。当体系带电时，还需手动实现偶极修正以处理长程库仑相互作用。force_and_stress函数则负责添加上文推导的原子核受力修正。occupancies函数在热电位计模拟中至关重要，它在每个离子步结束时被调用，用于根据热电位计公式更新电极电荷nelectrode，并相应地调整体系总电子数以保持电荷中性。通过设置INCAR文件中的相应标签，即可激活这些Python插件功能。
案例研究
带电表面吸附
以Au(111)表面为模型体系，展示了CDCE方法在模拟带电表面的能力。通过施加不同的电极电荷nelectrode，可以在金表面诱导出相应的静电偏压（电压）。研究计算了单个H原子在Au(111)表面不同吸附位点（顶位、hcp空心位、fcc空心位）的吸附能随外加电场的变化。结果表明，在零场下fcc空心位最稳定，但在电场达到约0.05 V/?时，hcp空心位在能量上变得更为有利。通过计算电荷密度差，可以可视化氢原子吸附引起的电荷转移，以及在施加外电位下Au-H键的进一步极化。此外，AIMD模拟显示，在Au(111)表面上的乙醛分子其构型会随表面电荷的改变而发生显著变化：表面带正电时，电负性的氧原子与表面强烈相互作用，距离较近；表面带负电时，氧原子被排斥，导致分子发生旋转。
场离子显微镜与原子探针层析成像
在原子探针层析成像（APT）等实验中，强电场会导致表面原子的场致脱附、电离和蒸发。研究模拟了具有台阶的Li(952)表面在300K下的动力学行为。在无外场时，台阶拐角处的原子保持束缚状态。然而，当施加1.25 V/?的均匀电场后，AIMD模拟清楚地捕捉到该拐角原子在约4.5皮秒内自发脱离台阶并在表面迁移的过程。这从原子尺度直接证实了强电场可以 destabilize 台阶/拐角原子，促进其以吸附原子形式释放，为理解APT实验机制提供了重要见解。
电化学固-液界面
将CDCE方法与热电位计技术相结合，可以实现对固-液界面的恒电位AIMD模拟。研究以Au(111)/水界面为例，在模拟过程中，Python插件根据设定目标电位动态调整金表面的电荷。模拟结果显示，电极电位和电荷均围绕目标值波动。施加的电位对界面水结构有显著影响：在负电位下，氧原子被金表面排斥，第一密度峰向界面外移动；在正电位下，表面电荷吸引氧原子，使水分子密度分布的主峰更靠近金属表面。这些结构变化伴随着界面振动光谱、电容和水层动力学性质的改变，突显了在电化学界面模拟中进行显式电位控制的重要性。
量子力学/分子力学隐式溶剂模型
利用修改外部势Vext的能力，可以方便地引入各种溶剂化模型。研究展示了一种基于VASP-Python接口的QM/平均场-MM隐式溶剂化工作流程。以水溶液中的Mg2+离子为例，首先对水环境进行经典分子动力学模拟，并对溶剂电荷密度进行平均，得到平均的溶剂电荷密度ρsolventMM。随后，将此电荷密度转化为静电势VsolventMM，并作为外部势Vext加入到Kohn-Sham哈密顿量中。结果显示，与真空中的孤立Mg2+离子产生的长程抛物线型电势不同，在隐式溶剂模型下，镁离子的+2电荷被周围水分子屏蔽，电势分布呈现出反映离子周围形成 distinct 溶剂化层的振荡特征。该方法可轻松扩展到表面和界面等更复杂的几何结构，用户可定义任意的QM区域，并将其与灵活的MM环境描述耦合，从而高度适应各种电化学和界面问题。
结论
本工作介绍了一种在DFT超晶格计算中施加任意静电势的灵活方案。该方案通过VASP-Python接口实现，允许用户提供任何所需的外部势Vext，并将其直接纳入Kohn-Sham哈密顿量。通过推导并应用必需的能量和力修正项，确保了总能量和原子力的物理正确性。通过表面吸附、场离子显微镜、电化学界面以及隐式溶剂模型等多个案例研究，充分展示了该方法在探究电场诱导现象方面的强大能力与广泛应用前景。该方法将外部势作为用户可控的输入，在保留VASP标准工作流程的同时，提供了对静电环境的完全控制，为在成熟的DFT框架内研究从催化、腐蚀到纳米器件物理等广泛领域中的场效应现象，提供了一个通用、易扩展的计算框架。
