摘要：aiida-vibroscopy 是 AiiDA 生态中专注于振动光谱计算的插件，基于 phonopy 与 Quantum ESPRESSO 实现红外、拉曼、非谐声子与介电张量计算，并通过 AiiDA 的 provenance 机制保证全流程可复现。本文系统梳理其核心功能、安装配置、命令行与 Python API 语法、关键参数含义，并给出 16 个可直接参考的实际应用案例，最后总结常见报错与使用注意事项。

一、aiida-vibroscopy 是什么

aiida-vibroscopy 是一个开源的 AiiDA 插件包，用于自动化计算固体材料的振动光谱性质。它把第一性原理计算中繁琐的“结构优化 → 力常数计算 → 声子谱 → 介电张量 → 红外/拉曼强度”流程封装为可复现的 AiiDA 工作流，用户只需提供结构，即可得到完整的光谱数据。

其核心依赖包括：

• AiiDA Core：负责工作流调度、数据溯源与结果存储。
• aiida-quantumespresso：提供 pw.x、ph.x 等计算插件。
• phonopy：处理有限位移与力常数、声子谱后处理。
• spglib：对称性分析，用于减少位移计算量。

二、核心功能概览

2.1 红外与拉曼光谱

通过计算 Born 有效电荷与介电张量，结合声子本征矢，得到红外吸收强度与拉曼散射张量，输出可直接绘图的频率-强度数据。

2.2 非谐声子与声子寿命

支持三阶力常数计算，进而评估非谐效应、声子寿命与晶格热导相关量。

2.3 介电张量与 Born 电荷

基于密度泛函微扰理论（DFPT）计算电子介电张量与原子 Born 有效电荷张量。

2.4 全流程可复现

每一步计算都作为 AiiDA 节点存储，可追溯输入结构、赝势、参数与输出，便于批量对比与论文复现。

三、安装与环境配置

3.1 前置条件

• Python 3.9 及以上版本
• 已安装并配置好的 AiiDA 环境
• 已配置 Quantum ESPRESSO 代码（pw.x、ph.x）
• phonopy、spglib 依赖

3.2 安装命令

# 建议在独立虚拟环境中安装 python -m venv aiiDA-vib-env source aiiDA-vib-env/bin/activate 安装插件 pip install aiida-vibroscopy 验证安装 verdi plugin list aiida.workflows | grep vibroscopy

3.3 配置计算代码

verdi code create core.code.installed \ --label qe-pw \ --computer localhost \ --remote-abs-path /opt/qe/bin/pw.x \ --default-calc-job-plugin quantumespresso.pw verdi code create core.code.installed --label qe-ph --computer localhost --remote-abs-path /opt/qe/bin/ph.x --default-calc-job-plugin quantumespresso.ph

四、语法与关键参数

4.1 命令行入口

aiida-vibroscopy --help aiida-vibroscopy run --structure Si.cif --protocol moderate

4.2 Python API 基本用法

from aiida import load_profile from aiida.orm import StructureData, Dict, load_code from aiida_vibroscopy.workflows.phonons.ir_raman import IRamanSpectraWorkChain load_profile() structure = StructureData.from_cif("Si.cif") pw_code = load_code("qe-pw") ph_code = load_code("qe-ph") builder = IRamanSpectraWorkChain.get_builder() builder.structure = structure builder.pw_code = pw_code builder.ph_code = ph_code builder.protocol = "moderate" builder.settings = Dict(dict={"displacement_generator": {"distance": 0.01}}) from aiida.engine import run results = run(IRamanSpectraWorkChain, **builder)

4.3 主要参数说明

五、16 个实际应用案例

案例 1：硅（Si）的拉曼光谱计算

金刚石结构硅是验证流程的标准体系。使用 moderate 协议计算 Γ 点光学声子，可得到约 520 cm⁻¹ 的拉曼峰，用于校验安装与参数是否正确。

builder.structure = StructureData.from_cif("Si.cif") builder.protocol = "moderate" builder.ir_raman = True

案例 2：金刚石的红外与拉曼联合计算

金刚石为高对称结构，仅 Γ 点有光学支。通过同时开启 ir_raman 与 dielectric，可验证 Born 电荷为零、拉曼张量各向同性的物理预期。

案例 3：NaCl 离子晶体的红外活性分析

NaCl 中 Born 有效电荷显著，红外吸收强度大。该案例适合演示如何从输出中提取 TO/LO 分裂与红外强度。

案例 4：TiO₂ 金红石相的声子谱

金红石相 TiO₂ 具有多个拉曼活性模式，是研究光催化材料振动性质的典型体系，可用于对比实验拉曼峰位。

案例 5：GaAs 极性半导体的 LO-TO 分裂

通过介电张量与 Born 电荷计算，可复现 GaAs 在 Γ 点的 LO-TO 分裂，验证非解析项修正是否正确。

案例 6：MgO 的介电张量计算

MgO 为高对称离子晶体，介电张量各向同性。该案例用于检查 DFPT 介电计算流程与收敛性。

案例 7：石墨烯的拉曼 G 峰与 2D 峰

二维材料需注意真空层厚度。通过设置足够大的 c 轴真空层，可计算石墨烯的 G 峰与 2D 峰位置。

案例 8：MoS₂ 单层的振动模式分析

单层 MoS₂ 的 E′ 与 A₁′ 模式对层数敏感，可用于判断层数与应变状态，是二维材料表征的经典案例。

案例 9：BaTiO₃ 铁电体的软模计算

铁电体软模与相变密切相关。该案例演示如何计算软模频率随温度或应变的演化趋势。

案例 10：ZnO 的拉曼光谱与缺陷峰

ZnO 的 E₂(high) 峰对晶格应力敏感，可用于分析薄膜应力状态与缺陷引入的额外峰。

案例 11：AlN 压电材料的 Born 电荷

AlN 具有强压电性，Born 有效电荷是理解压电响应的关键。该案例演示如何提取 Born 电荷张量。

案例 12：CaF₂ 萤石结构的红外模式

萤石结构具有典型的三重简并红外活性模式，适合作为教学案例演示红外强度计算。

案例 13：PbTe 热电材料的非谐声子

PbTe 的低热导率源于强非谐性。开启 anharmonic 参数计算三阶力常数，可评估声子寿命。

案例 14：SrTiO₃ 量子顺电体的介电响应

SrTiO₃ 的介电常数对温度敏感，该案例演示如何计算静态介电张量并与实验对比。

案例 15：二维 h-BN 的红外活性模式

h-BN 单层的 E₁ᵤ 模式为红外活性，需注意真空层与偶极修正设置，避免周期性镜像相互作用。

案例 16：批量筛选多种氧化物的拉曼峰位

结合 AiiDA 的 Group 与批量提交功能，可对一组氧化物结构批量运行工作流，并汇总拉曼峰位用于材料筛选。

from aiida.orm import Group, QueryBuilder qb = QueryBuilder() qb.append(StructureData, filters={"attributes.kinds": {"like": "%O%"}}) structures = [row[0] for row in qb.all()] for s in structures: b = IRamanSpectraWorkChain.get_builder() b.structure = s b.protocol = "fast" b.pw_code = pw_code b.ph_code = ph_code submit(IRamanSpectraWorkChain, **b)

六、常见错误与排查

6.1 赝势家族未安装

报错信息通常为Pseudo family not found。解决方式是执行aiida-pseudo install family安装对应赝势家族，并在 builder 中正确指定名称。

6.2 对称性导致位移数量异常

若位移数量远少于预期，通常是 spglib 识别出更高对称性。可通过 settings 中的对称性开关关闭自动对称化进行排查。

6.3 ph.x 计算不收敛

常见原因是截断能过低或 k 点过稀。建议先做收敛性测试，逐步提高 ecutwfc 与 k 点密度。

6.4 二维材料真空层不足

真空层小于 15 Å 时，周期性镜像会导致虚假相互作用，表现为声子频率异常或出现虚频。

6.5 内存或磁盘不足

三阶力常数计算量随原子数急剧增长，建议先用小体系验证流程，再扩展到大规模计算。

七、使用注意事项

• 先做收敛性测试：截断能与 k 点直接影响声子频率精度，务必先收敛再计算光谱。
• 结构预处理：输入结构应先充分弛豫，残余应力会导致虚频。
• 协议选择：fast 用于流程验证，moderate 用于常规研究，precise 用于发表级数据。
• 版本兼容：aiida-vibroscopy 与 aiida-quantumespresso 版本需匹配，升级前查阅发布说明。
• 结果校验：建议用已知体系（如 Si、NaCl）验证流程，再推广到新材料。
• 数据管理：定期备份 AiiDA profile 数据库，避免计算记录丢失。

八、总结

aiida-vibroscopy 把复杂的振动光谱计算流程标准化、可复现化，显著降低了第一性原理光谱计算的门槛。掌握其安装配置、核心参数与典型工作流后，可以高效完成从结构到红外/拉曼光谱的完整分析。建议从硅、NaCl 等标准体系入手，逐步扩展到二维材料、铁电体与热电材料等复杂体系。

《AI提示工程必知必会》为读者提供了丰富的AI提示工程知识与实战技能，主要包括各类提示词的应用，如问答式、指令式、状态类、建议式、安全类和感谢类提示词，以及如何通过实战演练掌握提示词的使用技巧；使用提示词进行文本摘要、改写重述、语法纠错、机器翻译等语言处理任务，以及在数据挖掘、程序开发等领域的应用；AI在绘画创作上的应用，百度文心一言和阿里通义大模型这两大智能平台的特性与功能，以及市场调研中提示词的实战应用。通过阅读《AI提示工程必知必会》，读者可掌握如何有效利用AI提示工程提升工作效率，创新工作流程，并在职场中脱颖而出。