单位体系与能量-力一致性新手完全指南)
【免费下载链接】pyraimdNew PYRAIMD updated from 2021 version developed from Nanoacademic Technology Inc.项目地址https://gitcode.com/gh_mirrors/py/pyraimd点击查看免费下载做分子动力学模拟能量差一个 Hartree、步长差一个数量级整条轨迹就会全盘报废。Pyramid包名pyraimd2是一个 Python 写的从头算分子动力学框架参考方法负责算能量和力机器学习势负责推进原子运动框架再决定何时检查、校验和更新模型。它的核心亮点之一是energetic force-error 方法——把力算得准不准直接和原子怎么动关联起来。对新手来说理解它的单位体系和**能量-力一致性energy-force consistency**检查是少走弯路的第一步。 30 秒认识 Pyramid (pyraimd2)一句话把第一性原理引擎PySCF、Quantum ESPRESSO、机器学习原子间势MLIP如 MACE和 MD 循环装进同一个可扩展的 Python 框架里。Python import 名是pyraimd2命令行工具是pyramid别名pyraimd2核心安装只需要 NumPy 和 ASE不需要 torch、PySCF 或外部程序新手路径见 docs/user_guide.md五分钟上手流程在 README.md一个完全离线、几秒跑完的玩具示例在 examples/harmonic_adaptive/建议先跑它再碰单位。 Pyramid 单位体系五个固定单位写在字段名里Pyramid 最省心的设计是单位直接编码在配置字段名里而且全项目固定不变不允许你另外指定物理量固定单位典型配置字段位置 / 位移Å埃, angstrom结构文件坐标能量 / 残余功eVenergy、max_energy_drift_eV力 / 力误差预算eV/Åforce_budget_eV_A、fmax_eV_A时间fs飞秒timestep_fs温度K开尔文temperature_K应力可选eV/ųASE Voigt 序后端返回这个约定在 docs/configuration.md 中明确写出Units are part of the field names and fixed everywhere代码层面同样是全局契约见 src/pyraimd2/config.py。好处很直接你看到timestep_fs 1.0它就是 1 飞秒不用再去查某个引擎的默认单位是秒还是飞秒。这也是 API 层的统一承诺——docs/api.md 开头就声明公开接口使用 Å、eV、eV/Å、fs 和 K。新手最该记住的几个带单位字段timestep_fsMD 步长fsMTSrespa里它是内层步长temperature_KNVT 浴的目标温度Kfriction_per_fsNVT 必需的浴耦合系数注意单位是每飞秒fmax_eV_A几何优化收敛阈值eV/Å默认 0.05force_budget_eV_A自适应 MD 的力误差接受预算eV/Åmax_wall_hours唯一的小时字段是墙钟预算而非物理时间完整字段表在 docs/configuration.md 的 Configuration reference 一节。 引擎边界上的单位换算Hartree、Bohr、Rydberg 都去哪了底层引擎各说各话PySCF 内部用Hartree/BohrQuantum ESPRESSO 用Rydberg。Pyramid 的做法是在引擎边界处一次性换算成 ASE 单位之后全链路只流通 eV 和 ÅPySCF 引擎位置除以units.Bohr送进去能量乘以units.Hartree、力乘以units.Hartree / units.Bohr换回来见 src/pyraimd2/engines/pyscf_engine.pyQE 引擎RY_EV units.Hartree / 2.0力从 Ry/Bohr 换算到 eV/Å见 src/pyraimd2/engines/qe_engine.py后端边界统一约定写在 docs/architecture.md 的 Units and configurations 一节给自写后端的启示实现引擎协议时convert units at that boundary and return energies consistent with the forces——换算发生在你的边界上返回的能量必须和力自洽架构文档原话见 docs/architecture.md。✅ 能量-力一致性如何检查力真的是能量的负梯度力误差预算再小前提是力来自一个保守势F -∇E。如果某个模型的力和能量对不上能量学检查锚定能量、端点功就没有意义。Pyramid 用两种机制兜底1. 能力声明capabilities。每个后端必须声明force_consistent、forces_conservative等属性。内置的谐振子玩具后端就诚实声明了两项全真src/pyraimd2/backends/harmonic.py跨能量类型如参考端报free_energy、代理端报energy的组合要求双方都声明force_consistentTrue且forces_conservativeTrue否则拒绝运行。2. 有限差分一致性检查。在线更新代理模型时GuardedUpdater会对候选模型做扰动探针沿两个内部方向相对坐标呼吸模 可复现伪随机方向做 ±step 有限差分比较数值导数与−F·u取最坏偏差有限差分步长energy_force_fd_step_A默认 1e-3 Å绝对容差energy_force_atol_eV默认 1e-3 eV相对容差energy_force_rtol默认 0.05超差即以energy_force_inconsistent拒绝候选并回滚父模型代码见 src/pyraimd2/loop/online.py 与默认参数定义 src/pyraimd2/loop/online.py。有个精巧的细节值得新手品味探针刻意不用整体平移方向。因为平移不变势的力总和为零只查平移方向会漏掉整体缩放力但能量不变的作弊模型——文档注释里写得很直白src/pyraimd2/loop/online.py。⚠️ 新手避坑清单不要猜单位看字段名。Pyramid 的字段名自带单位写timestep_fs 1就是 1 fs如果你从别的软件抄配置步长常以 ps 给出先换算。NVT 必须给friction_per_fs。它是每飞秒的耦合系数漏写会在 validate 阶段直接报错。calibrate-scale的配对数据要带force_unit eV/angstrom即力的单位必须是 eV/Åexamples/calibrate_scale/流程见 docs/user_guide.md。比较轨迹时时间必须精确匹配pyramid compare按物理时间逐点对齐容差是 1e-9 fs 级绝对容差且只读已完成轨迹docs/api.md。未知字段是错误不是静默忽略。配置里写错字段名比如timestep_fs2会在加载时立刻报出点号字段名和修复建议这是 src/pyraimd2/config.py 的设计原则。 继续深入五分钟离线跑通examples/harmonic_adaptive/配置字段与 CLI 语义docs/configuration.md架构与单位/边界约定docs/architecture.md公开 API 与单位声明docs/api.md能量学方法力误差 → 原子运动docs/energetic_force_error.md掌握五个固定单位 字段名即单位 力必须是能量负梯度这三件事你在 Pyramid 里就不会再被单位坑到了。赞分享【免费下载链接】pyraimdNew PYRAIMD updated from 2021 version developed from Nanoacademic Technology Inc.项目地址https://gitcode.com/gh_mirrors/py/pyraimd点击查看免费下载相关推荐10分钟装好Pyramid(pyraimd2)第一性原理分子动力学框架pip安装与可选后端避坑完全教程10分钟装好Pyramid pyraimd2 第一性原理分子动力学框架pip安装与可选后端避坑完全教程 PyramidPython 包名 pyraimd2NVT系综下的自适应分子动力学与模型在线更新Pyramid(pyraimd2)GuardedUpdater守护式更新完全指南NVT系综下的自适应分子动力学与模型在线更新Pyramid pyraimd2 GuardedUpdater守护式更新完全指南 Pyramid包名 pyraiPyramid(pyraimd2)架构完全拆解让能量力计算与动力学决策解耦的Engine与Surrogate协议Pyramid pyraimd2 架构完全拆解让能量力计算与动力学决策解耦的Engine与Surrogate协议 Pyramid 是一个 Python 从头算创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考
