在计算化学和物理化学领域，电离势（Ionization Potential, IP）和电子亲和势（Electron Affinity, EA）是两个核心概念，它们描述了原子、分子或离子在电子转移过程中的能量变化。这些性质不仅在理解物质的化学反应性、氧化还原行为、光谱特性等方面至关重要，还广泛应用于材料科学、有机合成和环境化学等领域。

1. 电离势（IP）的定义与物理意义

电离势是指从一个中性原子、分子或离子中移除一个电子所需的最小能量。这个过程本质上是氧化反应，因为失去电子会使系统变成阳离子。简单来说，它衡量了一个系统抓住电子的难易程度：IP值越高，意味着电子越难被移除，系统越稳定；反之则越容易发生电离。

从物理角度看，电离势源于量子力学中的电子轨道能级。根据 Koopmans 定理（一个近似理论），IP 可以粗略近似为最高占据分子轨道（HOMO）的负能量值，即 IP ≈ -ε_HOMO。这是因为 HOMO 是最外层电子所在的轨道，移除它所需的能量最低。但在实际中，这个近似有局限性，因为它忽略了电子重排和弛豫效应。

电离势可以进一步分为两种主要类型，以考虑几何结构的变化：

• 绝热电离势（Adiabatic Ionization Potential, AIP）：这是最常见的定义，指的是中性分子（或原子）从其优化几何结构（ground state）电离到阳离子优化几何结构所需的能量。计算公式为：AIP = E(阳离子, 优化几何) - E(中性分子, 优化几何)。这里的“绝热”意味着过程允许系统在电离后缓慢调整结构，以达到最低能量状态。这更接近实验中的热力学值，比如在气相中通过光电子谱（Photoelectron Spectroscopy）测得的电离能。
• 垂直电离势（Vertical Ionization Potential, VIP）：这是一种“瞬时”定义，指的是在中性分子的几何结构不变的情况下，直接移除电子所需的能量。公式为：VIP = E(阳离子 @ 中性几何) - E(中性分子, 优化几何)。这里的“垂直”源于 Franck-Condon 原理，类似于光谱中的垂直跃迁，不允许结构弛豫。VIP 通常比 AIP 高，因为阳离子在不优化的几何下能量更高。

在实际应用中，IP 值通常以电子伏特（eV）为单位报告。例如，对于水分子（H₂O），实验 AIP 约为 12.6 eV，这意味着需要这个能量才能从水中移除一个电子，形成 H₂O⁺。高 IP 的分子（如惰性气体）不易氧化，而低 IP 的（如碱金属）则易于失去电子。

2. 电子亲和势（EA）的定义与物理意义

电子亲和势是指一个中性原子、分子或离子捕获一个电子，形成阴离子时释放出的能量。它本质上是还原反应，因为获得电子会使系统带负电。EA 值越高（正值越大），意味着系统越容易接受电子，越具有氧化性；如果 EA 为负值，则表示捕获电子需要额外能量，过程不自发。

物理上，EA 可以近似为最低未占据分子轨道（LUMO）的负能量值，即 EA ≈ -ε_LUMO（同样基于 Koopmans 定理）。但这个近似对 EA 的误差更大，因为阴离子往往更不稳定，电子相关效应更显著。

类似于 IP，EA 也分为两种类型：

• 绝热电子亲和势（Adiabatic Electron Affinity, AEA）：定义为中性分子从其优化几何捕获电子，到阴离子优化几何时释放的能量。公式：AEA = E(中性分子, 优化几何) - E(阴离子, 优化几何)。这考虑了结构弛豫，是大多数文献和实验（如电子捕获解离谱）报告的值。正 AEA 表示过程放热，自发发生。
• 垂直电子亲和势（Vertical Electron Affinity, VEA）：指在中性分子的几何不变下，捕获电子释放的能量。公式：VEA = E(中性分子, 优化几何) - E(阴离子 @ 中性几何)。VEA 通常比 AEA 小（绝对值更小），因为阴离子在不优化的几何下能量更高，不如弛豫后稳定。

举例来说，氯原子（Cl）的 EA 约为 3.6 eV，这意味着 Cl 捕获电子形成 Cl⁻ 时释放大量能量，因此氯气是强氧化剂。相比之下，许多有机分子（如苯）的 EA 为负值，表示它们不易形成稳定阴离子。

IP 和 EA 的关系密切：一个系统的 IP 往往与另一个系统的 EA 相关联，形成氧化还原对。例如，在电化学中，IP 和 EA 可以帮助估算半电池电位或带隙能量（Band Gap = IP - EA，在半导体中）。

3. 在 Gaussian 软件中的计算原理与注意事项

前置准备工作（强烈建议先做好这些）

1. 准备好分子的初始几何结构（xyz 或 pdb 格式，用 GaussView 画或从数据库下载）。
2. 选择计算级别：

• 快筛/中等精度：ωB97X-D 或 ωB97M-V / def2-TZVP（带 GD3BJ 色散）或 6-311++G(d,p)
• 高精度：ωB97M-V / def2-QZVPPD 或 maug-cc-pVTZ（单点）
• 基准：DLPNO-CCSD(T)/CBS（用 ORCA 或其他软件更方便，Gaussian 也支持但慢）
  3. 阴离子计算特别容易不收敛，所以基组必须带 diffuse 函数（aug-、maug-、++G、def2-QZVPPD 等）。
  4. 所有计算都建议加上 int=ultrafine（网格更细，SCF 更稳）和 empiricaldispersion=gd3bj（色散校正）。
  5. 用同一套方法/基组计算所有物种，避免基组叠加误差（BSSE 在这里不明显，但能量一致性重要）。

一、计算绝热电离势（Adiabatic IP）和绝热电子亲和势（Adiabatic EA）——最推荐、最常见的做法

步骤 1：优化并频率计算中性分子（Neutral, N 电子）

这是起点，所有后续计算都依赖它的几何。

输入文件示例（以苯为例，ωB97X-D/def2-TZVP 级别）：

text

%chk=neutral.chk

%mem=16GB

%nprocshared=8

# opt freq ωB97X-D/def2-TZVP empiricaldispersion=gd3bj int=ultrafine

 

Neutral molecule optimization + freq

 

0 1

C   0.000   1.391   0.000

C   1.204   0.696   0.000

... （其余坐标）

• 运行后，确认无虚频（所有频率 >0），得到优化几何和零点能（ZPE）。
• 输出中找：Sum of electronic and zero-point Energies（简称 E_ZPE）或 Sum of electronic and thermal Enthalpies（带热校正的焓，298K 标准）。

步骤 2：优化并频率计算阳离子（Cation, N-1 电子，用于 IP）

用中性分子的优化几何作为起点（避免从头猜结构），最关键：用中性分子的波函数作为初猜（guess=read），大大提高收敛率。

输入文件：

text

%chk=cation.chk

%mem=16GB

%nprocshared=8

# guess=read geom=checkpoint opt freq ωB97X-D/def2-TZVP empiricaldispersion=gd3bj int=ultrafine SCF=xqc

 

Cation (IP calculation)

 

+1 2   ← 电荷 +1，自旋多重度 2（大多数情况为双重态）

• geom=checkpoint 自动读 chk 的最后优化几何。
• guess=read 读 chk 的波函数。
• 如果 SCF 不收敛，改成 SCF=qc 或 SCF=xqc（二次收敛）。
• 运行后确认无虚频，得到阳离子的 E_ZPE 或热焓。

步骤 3：优化并频率计算阴离子（Anion, N+1 电子，用于 EA）

同上，但电荷 -1，自旋通常 2。

输入文件：

text

%chk=anion.chk

%mem=16GB

%nprocshared=8

# guess=read geom=checkpoint opt freq ωB97X-D/def2-TZVP empiricaldispersion=gd3bj int=ultrafine SCF=xqc

 

Anion (EA calculation)

 

-1 2

• 阴离子最容易出问题：虚频、SCF 不收敛、自旋污染。
• 解决办法：用更大的 diffuse 基组、ultrafine 网格、xqc 算法；有时试试不同的泛函（如 TPSSh 对阴离子友好）。
• 确认无虚频，得到阴离子的 E_ZPE 或热焓。

步骤 4：提取能量并计算 IP/EA（单位 eV）

从三个 .log 文件的末尾提取：

• E_neutral = 中性分子的 Sum of electronic and zero-point Energies（或 thermal Enthalpies）
• E_cation = 阳离子的对应值
• E_anion = 阴离子的对应值

计算公式（Hartree 单位）：

• Adiabatic IP (AIP) = E_cation - E_neutral
• Adiabatic EA (AEA) = E_neutral - E_anion

单位转换：1 Hartree = 27.211386 eV，所以乘以 27.211386 得到 eV 值。

如果用了热焓（Sum of electronic and thermal Enthalpies），得到的是 298 K 下的焓差，更接近实验热值。

可选步骤 5：更高精度单点能

几何用 def2-TZVP 优化好后，用更大基组单点：

text

%chk=neutral_big.chk

# geom=checkpoint guess=read ωB97M-V/def2-QZVPPD empiricaldispersion=gd3bj int=ultrafine

 

Single point for IP/EA

对阳离子和阴离子也做同样单点。然后用这些单点电子能量 + 之前 TZVP 级别的 ZPE 校正，得到更准的结果。

二、计算垂直电离势（Vertical IP）和垂直电子亲和势（Vertical EA）

更简单、更快，只需优化中性分子。

步骤：

1. 只做步骤 1：优化中性分子 + 频率，得到几何和 ZPE。
2. 用中性几何，对阳离子做单点（不要 opt）：

text

%chk=cation_vert.chk

# guess=read geom=checkpoint ωB97X-D/def2-TZVP empiricaldispersion=gd3bj int=ultrafine

 

Vertical cation single point

 

+1 2

3. 同理，对阴离子做单点：

text

%chk=anion_vert.chk

# guess=read geom=checkpoint ωB97X-D/def2-TZVP empiricaldispersion=gd3bj int=ultrafine

 

Vertical anion single point

 

-1 2

4. 计算：

• VIP = E(cation 单点) - E(neutral)
• VEA = E(neutral) - E(anion 单点)

垂直值通常比绝热值大 0.1–0.5 eV（阳离子弛豫后更稳定）。

三、常见问题

• 阴离子不收敛 → 加大 diffuse（用 aug-cc-pVDZ 先试试）、int=ultrafine、SCF=(xqc,maxcycle=200)、或换泛函。
• 自旋污染 → 检查输出 ⟨S²⟩，正常双重态应 ≈0.75，若远大于此，说明污染严重，可能需用限制开壳或投影方法。
• 能量取哪个？ → 文献高精度用 Zero-point corrected；热力学性质用 thermal Enthalpies；快筛直接用电子能量也行。
• 单位换算：Gaussian 输出是 Hartree，乘 211386 得 eV；或乘 627.5095 得 kcal/mol。
• 批量计算 → 用脚本提取能量（Multiwfn 或自写 python 脚本）。