在本教程中，您将对平行板纳米电容器进行原子尺度的模拟，这意味着两个平行金属表面之间的间隙在纳米范围内。这目前不是一个真实的器件，目的只是为了展示一种可用于计算其他结构电容的方法，包括金属/半导体/金属结构。

打开 Builder ，按照以下步骤操作：

• 点击 Add From Database，添加金到 Stash。
• 打开 Builders Cleave plugin，输入Miller 指数为 [111]，创建 Au(111) 电极。
• 利用 Bulk Tools Repeat，将结构在 C 方向上扩至 6 倍。
• 选中并删除中间的两个原子—它们的原子指数为 3 和 14。
• 选中表面的原子（指数为 7 和 8），用工具栏的按钮 将其转变为鬼原子。
• 两表面间的距离还是有点太小。为更改间距，可打开 Coordinate Tools Translate，选中所有右手边表面上的原子（包含离得最近的鬼原子），将这些原子在 Z 方向上移动 + 2 Å。
• 点击 Device Tools Device From Bulk，采用默认设置，将块体结构转变为器件构型。

最终的器件如下图所示。

点击 按钮将器件构型发送到 Script Generator ，按如下所示设置脚本：

• 添加一个 New Calculator ，采用 DFT-LDA 引擎和 SingleZetaPolarized 基组。K 点网格设为 9 × 9 × 100。

• 添加以下分析项目 ：MullikenPopulation, ElectrostaticDifferencePotential 和 ElectronDifferenceDensity。
• 将输出文件名改为 au_vacuumgap.nc。
• 保存脚本为 au_vacuumgap.py。

发送脚本到 Job Manager ，运行计算。在笔记本电脑上需要 15-20 分钟。

为了得到电容，计算 $dq/dV$ 是很有必要的，$dq$ 是表面上的感应电荷，$V$ 是偏压。利用以下脚本执行偏压分别为 0.25 、0.5 、0.75 和 1.0 V 时的计算。

device_configuration = nlread('au_vacuumgap.nc', DeviceConfiguration)[-1]
calculator = device_configuration.calculator()
for voltage in [0.25, 0.5, 0.75, 1.0]:
    # Set new calculator with modified electrode voltages on the configuration
    # use the self consistent state of  the old calculation as starting input.
    device_configuration.setCalculator(
          calculator(electrode_voltages=(0.5*voltage*Volt, -0.5*voltage*Volt)),
                     initial_state=device_configuration)
    device_configuration.update()
    nlsave('au_vacuumgap.nc', device_configuration, object_id="SCF %s" % voltage)
    # -------------------------------------------------------------
    # Mulliken population
    # -------------------------------------------------------------
    mulliken_population = MullikenPopulation(device_configuration)
    nlsave('au_vacuumgap.nc', mulliken_population, object_id="Mulliken %s" % voltage)
    # -------------------------------------------------------------
    # Electrostatic difference potential
    # -------------------------------------------------------------
    electrostatic_difference_potential = ElectrostaticDifferencePotential(device_configuration)
    nlsave('au_vacuumgap.nc', electrostatic_difference_potential, object_id="EDP %s" % voltage)
    # -------------------------------------------------------------
    # Electron difference density
    # -------------------------------------------------------------
    electron_difference_density = ElectronDifferenceDensity(device_configuration)
    nlsave('au_vacuumgap.nc', electron_difference_density, object_id="EDD %s" % voltage)

您可以在此处下载脚本：↓bias_loop.py。用 Job Manager 运行或从终端执行：

$ atkpython bias_loop.py > bias_loop.out

在笔记本电脑上将需要约 1 小时完成计算。

我们将介绍两种不同的方法计算电容：

• 一种基于感应电荷的简便方法；
• 一种基于评估静电能量的更普遍方法。

有限偏压将在两个金表面上形成正/负电荷网络，如下图所示。

现在，原则上电容可以由公式 $C=dq/dV$ 计算得到，$q$ 为偏差函数在一个表面上的总累积电荷。因此，从直线的斜率就可以获得纳米电容器的电容：

上图和基于 Mulliken 群的整个数据分析都可以通过运行脚本：↓mulliken-analysis.py 实现。脚本很长，但这是一个学习 Python 和 QuantumATK 中自定义数据分析的绝佳机会。

脚本显示由数据点拟合的直线得到电容为 5.77*10-22 法拉。为了正确全面地看待这个数字，需要考虑理想平行板电容器，可表达为 $C = \varepsilon A/d$，此处 $vacuum permittivity $ 为真空介电常数，$A$ 为由器件结构计算得到的每个平板的面积。然后您就可以计算平行板间距 $d=A/C$。该步骤也在脚本中完成，并显示结果为 $d$ = 11.05 Å。

将该值与表面金原子之间的距离（13.77 Å）进行比较，所谓图像平面的位置估计高于金表面 1.36 Å。

要使用更为通用和准确的方法计算电容，您需要首先计算静电能量 $E$，其以偏压 $V$ 为变量的函数表达式为：

$$ E(V)=\frac{1}{2}\int \delta v(\mathbf{r},V) \delta n(\mathbf{r}, V)d\mathbf{r} $$

在这里，$\delta v$ 为诱导静电势，$\delta n$ 为诱导密度。

利用脚本 ↓ potential_plot.py 从已存在的 NetCDF 文件中提取 1 V 偏压下的电势和密度，并生成诱导密度和静电势的图。脚本还报告这两个参数：

• 真空间隙中间的电场 $E$ = 0.0905 V/Å；
• 两图像平面间距 $d$ = 11.88 Å。

蓝色圆圈对应于原子位置。主要的电荷累积发生在表面正上方非常狭窄的区域，外表面原子和鬼原子之间。

由此，可以将抛物线按照公式 $E(V) = CV^2/2$ 拟合后提取出电容。用脚本 ↓ v-e-plot.py 计算静电能量并拟合得到电容，如下所示。

脚本报告的电容为 $C$ = 5.45*10^-22 法拉，与采用 Mulliken 电荷分析的简单方法得到的电容值一致。

脚本还可以报告 $d$ = 11.69 Å，对应于金表面上方 1.04 Å 处的图像平面，非常接近于诱导密度峰值的位置。