1. 引言aiida-qeq 是一个基于 AiiDA 工作流框架的 Python 包用于自动化调用 QEQ电荷平衡算法为分子和材料体系分配原子电荷。它封装了 QEQ 程序的输入输出将电荷分配任务无缝集成到 AiiDA 的数据流和计算图管理中适合需要高通量计算原子电荷的研究人员。本文将从功能、安装、语法、参数、实际案例以及常见错误等方面系统介绍 aiida-qeq 的使用方法。2. 核心功能aiida-qeq 的核心功能包括QEQ 计算封装将 QEQ 可执行程序封装为 AiiDA 计算插件支持自动生成输入文件、解析输出文件。电荷分配基于电负性均衡原理为体系中的每个原子分配部分电荷。工作流集成提供完整的工作流WorkChain可串联结构优化、电荷计算和结果存储。数据节点管理自动将输入结构、输出电荷、日志等保存为 AiiDA 数据节点便于追溯和复用。多力场支持支持多种电负性参数集适配不同分子力场。3. 安装方法aiida-qeq 依赖 AiiDA 核心库和 QEQ 可执行程序。推荐使用 conda 或 pip 安装。3.1 安装 AiiDA 核心conda create -n aiida-qeq python3.9 conda activate aiida-qeq pip install aiida-core3.2 安装 aiida-qeqpip install aiida-qeq3.3 安装 QEQ 可执行程序QEQ 程序通常随分子模拟软件如 LAMMPS一起分发。需要将可执行文件路径配置到 AiiDA 中reentry scan verdi quicksetup verdi code create core.code.installed --label qeq --description QEQ executable --computer localhost --filepath /path/to/qeq4. 基本语法与参数aiida-qeq 的核心接口是QeqCalculation计算类和QeqWorkChain工作流类。4.1 计算类参数参数类型说明structureStructureData输入原子结构必填parametersDictQEQ 运行参数必填codeCodeQEQ 可执行程序必填qeq_parametersDict电负性参数集可选4.2 常用运行参数参数键默认值说明cutoff10.0近程相互作用截断距离埃tol1e-6电荷收敛容差maxiter200最大迭代次数qeq_scale1.0电荷缩放因子4.3 工作流调用示例from aiida.orm import load_code, Dict, StructureData from aiida_qeq.workchains import QeqWorkChain structure StructureData() ... 从文件或数据库构建结构 params Dict(dict{ cutoff: 10.0, tol: 1e-6, maxiter: 200, }) code load_code(qeqlocalhost) builder QeqWorkChain.get_builder() builder.structure structure builder.parameters params builder.code code from aiida.engine import run result run(builder) print(电荷结果节点:, result[charge])5. 9 个实际应用案例5.1 案例一水分子电荷分配对单个水分子进行 QEQ 电荷计算验证氧原子和氢原子的电荷分布是否符合预期。from aiida.orm import StructureData, Dict, load_code from aiida.engine import run from aiida_qeq.calculations import QeqCalculation 构建水分子结构 structure StructureData() structure.append_atom((O, (0.0, 0.0, 0.0))) structure.append_atom((H, (0.957, 0.0, 0.0))) structure.append_atom((H, (-0.240, 0.927, 0.0))) params Dict(dict{cutoff: 8.0, tol: 1e-6}) code load_code(qeqlocalhost) builder QeqCalculation.get_builder() builder.structure structure builder.parameters params builder.code code result run(builder) print(原子电荷:, result[charge].get_array(charges))5.2 案例二甲烷分子电荷计算计算甲烷分子中碳原子和氢原子的 QEQ 电荷用于力场参数验证。from ase.build import molecule from aiida.orm import StructureData, Dict, load_code from aiida.engine import run from aiida_qeq.calculations import QeqCalculation ase_struct molecule(CH4) structure StructureData(asease_struct) params Dict(dict{cutoff: 10.0, tol: 1e-7, maxiter: 300}) code load_code(qeqlocalhost) builder QeqCalculation.get_builder() builder.structure structure builder.parameters params builder.code code result run(builder) charges result[charge].get_array(charges) print(甲烷电荷:, charges)5.3 案例三乙醇分子电荷分配对乙醇分子进行电荷计算分析羟基氧和氢的电荷差异。from ase.build import molecule from aiida.orm import StructureData, Dict, load_code from aiida.engine import run from aiida_qeq.calculations import QeqCalculation ase_struct molecule(C2H5OH) structure StructureData(asease_struct) params Dict(dict{cutoff: 12.0, tol: 1e-6}) code load_code(qeqlocalhost) builder QeqCalculation.get_builder() builder.structure structure builder.parameters params builder.code code result run(builder) charges result[charge].get_array(charges) print(乙醇电荷:, charges)5.4 案例四二氧化硅晶体电荷计算对 α-石英晶胞进行 QEQ 电荷计算验证周期性体系的电荷分布。from aiida.orm import StructureData, Dict, load_code from aiida.engine import run from aiida_qeq.calculations import QeqCalculation 从 CIF 文件读取结构 from aiida.orm import CifData cif CifData(file/path/to/quartz.cif) structure cif.get_structure() params Dict(dict{cutoff: 10.0, tol: 1e-6, maxiter: 500}) code load_code(qeqlocalhost) builder QeqCalculation.get_builder() builder.structure structure builder.parameters params builder.code code result run(builder) charges result[charge].get_array(charges) print(SiO2 电荷:, charges)5.5 案例五金属有机框架MOF节点电荷计算对 MOF 材料的次级构筑单元进行电荷分配用于后续吸附模拟。from aiida.orm import StructureData, Dict, load_code from aiida.engine import run from aiida_qeq.calculations import QeqCalculation 假设已从数据库加载 MOF 结构 structure load_node(mof_structure_uuid) params Dict(dict{cutoff: 15.0, tol: 1e-7, maxiter: 1000}) code load_code(qeqlocalhost) builder QeqCalculation.get_builder() builder.structure structure builder.parameters params builder.code code result run(builder) charges result[charge].get_array(charges) print(MOF 节点电荷:, charges)5.6 案例六聚合物链段电荷分配对聚乙烯链段进行 QEQ 电荷计算验证碳氢电荷分布规律。from ase.build import molecule from aiida.orm import StructureData, Dict, load_code from aiida.engine import run from aiida_qeq.calculations import QeqCalculation 构建聚乙烯链段C10H22 from ase import Atoms atoms Atoms(C10H22) 使用 ASE 的几何构建工具生成合理构象 from ase.build import fcc111, add_adsorbate 此处简化实际应使用分子构建工具 structure StructureData(aseatoms) params Dict(dict{cutoff: 10.0, tol: 1e-6}) code load_code(qeqlocalhost) builder QeqCalculation.get_builder() builder.structure structure builder.parameters params builder.code code result run(builder) charges result[charge].get_array(charges) print(聚乙烯链段电荷:, charges)5.7 案例七离子液体阴阳离子电荷计算对 1-丁基-3-甲基咪唑六氟磷酸盐离子对进行电荷分配。from aiida.orm import StructureData, Dict, load_code from aiida.engine import run from aiida_qeq.calculations import QeqCalculation 从 XYZ 文件读取离子对结构 from aiida.orm import SinglefileData xyz SinglefileData(file/path/to/ion_pair.xyz) structure StructureData() 解析 XYZ 并构建结构此处省略解析代码 params Dict(dict{cutoff: 12.0, tol: 1e-7, maxiter: 800}) code load_code(qeqlocalhost) builder QeqCalculation.get_builder() builder.structure structure builder.parameters params builder.code code result run(builder) charges result[charge].get_array(charges) print(离子对电荷:, charges)5.8 案例八表面吸附体系电荷重分配计算小分子在金属表面吸附后的电荷重新分布用于催化机理研究。from ase.build import fcc111, add_adsorbate from ase import Atoms from aiida.orm import StructureData, Dict, load_code from aiida.engine import run from aiida_qeq.calculations import QeqCalculation 构建 Pt(111) 表面并吸附 CO slab fcc111(Pt, size(2, 2, 4), vacuum10.0) co Atoms(CO, positions[(0.0, 0.0, 2.0), (0.0, 0.0, 2.8)]) add_adsorbate(slab, co, 2.0, ontop) structure StructureData(aseslab) params Dict(dict{cutoff: 10.0, tol: 1e-6, maxiter: 600}) code load_code(qeqlocalhost) builder QeqCalculation.get_builder() builder.structure structure builder.parameters params builder.code code result run(builder) charges result[charge].get_array(charges) print(吸附体系电荷:, charges)5.9 案例九高通量电荷计算工作流批量计算多个分子结构的电荷并将结果汇总存储。from aiida.orm import Dict, load_code, load_group from aiida.engine import run from aiida_qeq.workchains import QeqWorkChain code load_code(qeqlocalhost) group load_group(molecules_set) results [] for node in group.nodes: builder QeqWorkChain.get_builder() builder.structure node builder.parameters Dict(dict{cutoff: 10.0, tol: 1e-6}) builder.code code result run(builder) results.append((node.uuid, result[charge].get_array(charges))) for uuid, charges in results: print(f结构 {uuid}: {charges})6. 常见错误与使用注意事项6.1 常见错误错误现象可能原因解决方法找不到 QEQ 可执行文件Code 未正确配置使用verdi code show检查路径电荷不守恒参数qeq_scale设置不当检查缩放因子确保总电荷为零计算不收敛maxiter过小或tol过严增大迭代次数或放宽容差结构解析失败结构文件格式不支持转换为 CIF 或 XYZ 格式周期性边界错误晶胞参数缺失确保 StructureData 包含正确的晶胞6.2 使用注意事项参数合理性cutoff应大于体系最大原子间距否则会遗漏近程相互作用。电荷单位QEQ 输出的电荷单位为元电荷 e注意与力场单位换算。结构完整性输入结构必须包含所有原子坐标和元素种类周期性体系需提供晶胞向量。参数集匹配不同力场需要匹配对应的电负性参数集混用会导致电荷异常。版本兼容aiida-qeq 对 AiiDA 核心版本有要求升级前请查阅官方文档。结果验证建议对已知体系如水、甲烷进行测试确认电荷分布符合化学直觉。数据追溯所有计算结果自动存入 AiiDA 数据库可通过verdi node show查看完整 provenance。7. 总结aiida-qeq 将 QEQ 电荷分配算法与 AiiDA 工作流引擎深度集成为分子和材料模拟提供了可追溯、可复现的电荷计算方案。通过本文介绍的 9 个案例读者可以快速上手从单分子到周期性体系、从单次计算到高通量扫描的各类场景。建议在实际项目中先验证参数合理性再批量运行并善用 AiiDA 的数据管理能力追踪每一步计算。《动手学PyTorch建模与应用:从深度学习到大模型》是一本从零基础上手深度学习和大模型的PyTorch实战指南。全书共11章前6章涵盖深度学习基础包括张量运算、神经网络原理、数据预处理及卷积神经网络等后5章进阶探讨图像、文本、音频建模技术并结合Transformer架构解析大语言模型的开发实践。书中通过房价预测、图像分类等案例讲解模型构建方法每章附有动手练习题帮助读者巩固实战能力。内容兼顾数学原理与工程实现适配PyTorch框架最新技术发展趋势。