AI Scientist Platform 实战:苯乙基间苯二酚的分子轨道与静电势
以苯乙基间苯二酚的分子轨道与静电势分析任务为例,了解 Mira 如何结合项目、Skills、任务与文件,完成构象采样、DFT、HOMO/LUMO 和 ESP 工作流。
这次 AI Scientist Platform 演示记录 Mira 如何把苯乙基间苯二酚的分子轨道与分子表面静电势任务组织为可审阅的研究记录:项目保留问题与上下文,Skills 完成专业步骤,任务执行计算,文件保留构象、DFT、前线分子轨道和分子表面静电势结果。它展示平台工作流,不针对产品功效或安全性作出结论。
量子化学任务常常需要在结构、计算工具、输出文件和下一步研究计划之间切换。Mira 将这些步骤组织为可保存的工作流:确认输入和计算边界,得到苯乙基间苯二酚的轨道与静电势结果,再由研发人员决定哪些假设值得进入下一轮计算或实验。

截图展示了任务问题、项目文件、Skills 调用、任务状态和输出产物如何保留在同一个工作台中。
作者:Mira · 发布于:2026-07-14
本文依据:一次已完成的 Mira 任务记录及其输出文件编写。计算在气相模型下完成,结果用于展示量子化学工作流;不构成对特定商业原料、产品功效或安全性的结论。
苯乙基间苯二酚计算案例
本次任务要求 Mira 分析苯乙基间苯二酚的分子轨道信息和静电势。该任务采样 20 个构象,以 wb97x-d4/def2-svpd 优化选定结构,并输出 HOMO/LUMO 与 ESP 结果。本文以这次具体计算说明平台工作流。
平台如何承载一次计算?
这次任务使用了四项平台能力:
- 项目:保留问题、假设和任务历史。
- Skills:完成输入准备、构象采样、计算和分析等专业步骤。
- 任务:以一次可追踪的执行串起结构与结果。
- 文件:保留坐标、日志、波函数、CSV 和图像,供团队复查。
研究人员可以复查输入、方法、输出和结论边界,而不是只接收一段没有证据的回答。
演示输入与计算链路
本次演示使用了一个间苯二酚衍生物模型:SMILES 为 Oc1cc(O)cc(CCc2ccccc2)c1,分子式为 C₁₄H₁₄O₂,体系为中性单重态。这里展示输入,是为了说明任务实际计算的对象;它不用于标识任何商业原料。输入结构、电荷和自旋决定了后续构象采样、DFT 计算和结果解释的边界。
任务按以下步骤完成:
- 确认结构:解析 SMILES,并生成 2D 结构图供核对。
- 构象采样:使用 CREST 和 GFN2-xTB 得到 20 个构象,选择最低能量构象进入 DFT。
- 几何优化与计算:GPU4PySCF 以 wb97x-d4/def2-svpd 完成 DFT 几何优化和单点计算,并导出 Molden 波函数文件。
- 轨道与静电势分析:用 Multiwfn 生成 HOMO/LUMO、ESP 和 cube 文件,再由 VMD 渲染轨道与 ESP 图像。
- 文件归档:保留构象坐标、优化日志、Molden 波函数、Multiwfn 数据表、cube 文件、VMD 图片和渲染脚本。
CREST、GPU4PySCF、Multiwfn 和 VMD 分别完成科学计算与可视化。Mira 负责提交计算节点、衔接中间文件、调用批处理脚本、跟踪长任务并整理结果文件。研究人员可以复查计算记录,或自行调整可视化参数。
本次演示的计算输出
下表只汇总模型分子的计算输出,不将它们外推为产品表现或生物活性。
| 指标 | 本次结果 | 用途 |
|---|---|---|
| 构象采样 | 20 个构象 | 选择可说明的 DFT 起点 |
| 最低 CREST 构象能量 | -45.24634655 Hartree | 后续 DFT 的输入 |
| DFT 优化方法 | wb97x-d4/def2-svpd | 记录计算层级,便于复查 |
| 优化后总能量 | -691.808182 Hartree | 当前气相模型下的计算记录 |
| HOMO 能量 | -0.3082520 Eh,约 -8.39 eV | 同一模型内的电子结构比较 |
| LUMO 能量 | 0.0611860 Eh,约 1.67 eV | 与 HOMO 一起阅读 |
| HOMO-LUMO 能隙 | 约 10.05 eV | 同一模型内的结构描述符 |
| ESP 表面范围 | -26.0 至 50.0 kcal/mol | 分子表面相对静电环境 |
两个羟基氧原子附近是本次结果中较明显的富电子区域,平均表面静电势约为 -18.6 和 -18.5 kcal/mol。整个表面中,正电势区域约占 58%,负电势区域约占 42%。这些数据可用于选择后续的对照模型或相互作用问题。
如何阅读 HOMO、LUMO 和 ESP?
HOMO 可以帮助观察相对容易参与给电子过程的电子分布区域,LUMO 则提供潜在受电子区域的一个视角。两者的能量和空间分布适合在相同计算层级下比较相关模型分子。
分子表面静电势(ESP) 把不同区域的静电环境映射到分子表面。羟基氧附近的负电势区域值得被标记下来,但它不等于已经证明了某个靶点相互作用或真实世界结果。溶剂、浓度、pH 和实验条件都会影响真实体系。
平台如何支持复用?
单次计算的价值,取决于团队能否复查和复用。一个稳妥的流程通常包括:
- 写清问题和计算边界。
- 记录 SMILES、电荷和自旋。
- 保留构象数、方法、日志、波函数和图表。
- 对相关模型分子做同一层级的平行计算。
- 由研究人员决定是否加入溶剂模型、不同方法或实验验证。
PySCF 等电子结构工具让可编排的量子化学流程成为可能。Mira 在项目中保存输入与输出,调用专业 Skills,通过任务执行较长计算,再将文件和限制条件交给下一位审阅者。
如果你也关心模型选择如何影响 Agent 工作流,可阅读 GLM-5.2 免费体验指南。
研究人员仍需做什么?
一次气相计算不能替代实验验证、安全性评价或领域特定模型。它也不能确认具体靶点相互作用、真实世界表现或法规合规。
研究人员仍要判断输入是否正确、方法是否适合问题,以及结果是否支持下一步。准备好用自己的研究问题验证这条链路时,可以创建 Mira 项目。
常见问题
这次计算能直接证明模型分子的真实世界性质吗?
不能。它提供构象、轨道和静电势等分子层面的信息,帮助提出可检验的假设;任何真实世界结论仍需使用合适的模型和实验验证。
为什么要先做构象采样,再做 DFT?
柔性分子可能存在多个合理构象。先采样可以减少把任意初始结构直接送进高精度计算的风险。本次演示从 20 个构象中选择最低能量构象作为 DFT 起点。
HOMO-LUMO 能隙越小,模型分子就一定越好吗?
不是。能隙只是在既定模型和计算条件下描述电子结构的一个指标。它必须与分子结构、目标问题、环境条件以及实验结果一起解读。
Mira 如何帮助团队复用这类计算?
Mira 可以把项目问题、输入结构、Skills、任务过程、输出文件和审阅结论保留在同一处。团队随后可沿用相同流程分析相关模型分子,并在同一上下文中比较结果。

