《Physical Chemistry Chemical Physics》：Perspectives on theoretical tools for studying molecular anionsOpen Access
编辑推荐：
本文围绕研究人员所在研究组参与发展的四类理论化学工具进行思考，这些工具已被证明可用于分子阴离子研究。四类工具包括：(i) 直接计算强度量电子亲和能（EA）的运动方程（EOM）方法，而非由广延量的中性分子能量减去广延量的阴离子能量；(ii) 为研究偶极束缚阴离子
本文围绕研究人员所在研究组参与发展的四类理论化学工具进行思考，这些工具已被证明可用于分子阴离子研究。四类工具包括：(i) 直接计算强度量电子亲和能（EA）的运动方程（EOM）方法，而非由广延量的中性分子能量减去广延量的阴离子能量；(ii) 为研究偶极束缚阴离子而系统构建额外弥散原子轨道基函数的方法；(iii) 由振动/转动能量转移导致电子发射的速率及相关角分布的解析表达式；(iv) 利用常规量子化学工具确定电子亚稳阴离子的能量和寿命的方法。本文还对这些研究领域未来值得开展的改进提出建议，期望能启发更年轻的理论化学家从事相关工作。
研究背景方面，分子阴离子（MAs）的电子结合能通常较小，EA远小于电离势（IP），早期电子结构方法计算EA的误差常达±0.5 eV，难以满足要求。其原因包括：阴离子电子云更弥散，需要额外径向弥散原子轨道基函数；阴离子—中性相关能差必须高精度处理；传统方法分别计算广延量的中性总能量与阴离子总能量再相减，计算成本高且存在规模依赖。实验上，光电脱附可一步直接测量EA，这促使研究人员发展可与实验测量方式相对应的理论工具。研究人员回顾了50余年来研究组参与发展的MAs理论工具，认为这些工具推动了弱束缚电子体系研究，但仍需在多参考处理、基组构建、亚稳态外推可靠性等方面改进，其意义在于为MAs提供直接、平衡且可与实验对接的量子化学框架。
关键技术方法如下（不超过250字）：采用EOM（equations of motion，运动方程）、GF（Green's function，格林函数）及CC（coupled-cluster，耦合簇）参考的EOM-EA方法直接求EA；以RSPT（Rayleigh–Schr?dinger perturbation theory，Rayleigh–Schr?dinger微扰理论）和CC波函数构造Ω算符并求解矩阵本征值或迭代方程。偶极束缚阴离子采用aug-cc-pVNZ加额外弥散s、p（必要时d）函数，偶合指数约3.0~5.0。亚稳态采用稳定化图与GPA（generalized Padé approximate，广义Padé近似）、CAP（complex absorbing potential，复吸收势）、Coulomb/非Coulomb稳定化及ACCC（analytic continuation of the coupling constant，耦合常数解析延拓）/RAC（regularized analytic continuation，正则化解析延拓）。无样本队列来源。
研究结果如下。II. 直接EA计算方法：II.A 动机表明，早期EA计算受弥散基函数和相关能差限制，研究人员希望像光电脱附实验一样直接计算强度量EA。II.B EOM理论用于EA的起点表明，通过选取RSPT形式的中性参考波函数ψ
neutral、以二次量子化算符构造Ω并导出反对易子矩阵方程，可迭代求得EA；零阶对应Koopmans定理，高阶项包含轨道弛豫和差分电子相关。II.C killer条件与分区步骤表明，RSPT参考函数不严格满足Ω
+ψ
neutral=0，但满足killer条件的参考函数与RSPT系数差异小，仍可保证EA在三阶精度；分区步骤假定少数c系数主导，对某些多主导构型并不可靠。II.D 与GF理论的联系表明，EOM-EA方程可在多体GF框架下推出，等价于三阶电子传播子；早期小线性阴离子计算用于检验额外基函数和三阶贡献，Poul J?rgensen和Ken Jordan的加入推动了GF联系及偶极束缚（DB）阴离子研究。II.E 重大改进与遗留工作表明，Marcel Nooijen、Rodney Bartlett、John Stanton、So Hirata、Piotr Piecuch和Anna Krylov等发展的CC-EOM方法更强大且保持大小广延性，并扩展到IP和激发能；Vincent Ortiz和So Hirata推进了GF理论；多参考及双主导构型问题仍需改进。III. 偶极束缚阴离子的特殊工具：1996年与Ken Jordan、Alexander Boldyrev、Maciej Gutowski和Piotr Skurski等表明，DB阴离子外层弥散轨道可极化，色散样相关贡献重要，Koopmans定理（KT）级不可靠；2000年与Piotr Skurski、Maciej Gutowski等提出以aug-cc-pVNZ为起点系统添加额外弥散s、p（有时d）函数的策略，并指出ghost原子用于几何导数时需谨慎。IV. 亚稳分子阴离子的工具：IV.A 振动—转动诱导电子发射表明，1981年推导的速率表达式由电子耦合项和振动波函数导数耦合项决定，Condon样近似下主要给出δv=±1跃迁；振动坐标Q的对称性约束电子跃迁矩阵元和角分布。IV.B 本征亚稳分子阴离子电子态表明，稳定化图可提取能量和线宽，简单线性双diabatic模型不足，GPA以高阶多项式/有理分式改进；Spiridoula Matsika等将其用于尿嘧啶阴离子π
*形状共振（shape resonance），多个避免交叉需高阶或多态GPA；CAP、Coulomb/非Coulomb稳定化及ACCC/RAC可外推得到复数EA与半宽Γ/2，但中间离散连续态（DC）可能影响可靠性。V. 简要总结与讨论表明，理论工具应聚焦于与实验测量方式相连的量；DB阴离子被认为可能存在于星际空间并参与界面电子捕获，相关贡献不可忽略；MAs电子结合能与振动/转动能级间隔相近，需理解能量和角动量耦合；本征亚稳电子态已有常规电子结构码确定能量和寿命的方法，但仍需改进。
研究结论翻译：第一，将理论工具聚焦于与实验测量方式相连的量是有用的；EA实验可直接测量阴离子—中性能量差，理论也应如此，这是响应方法的核心，EOM和GF方法重要但仍需改进。第二，偶极束缚阴离子目前被认为可能存在于星际空间，也可能在极性分子液—气或固—气界面捕获电子时起作用；其量子化学研究需要特殊基组，且不能忽略相关贡献。第三，MAs的电子结合能常与振动能级间隔相近，偶极束缚阴离子还可有相近的转动能级间隔；因此需要理解振动/转动运动与弱束缚电子运动之间的能量和角动量耦合，包括电子脱附的倾向、速率和角分布。第四，对于本征亚稳的电子态，已有多种量子化学方法可利用常规电子结构程序确定其能量和寿命；研究人员强调了常用方法并提出仍需改进的方向。
