二维硫化亚锗（GeS）凭借独特褶皱晶体结构与优异光电特性，在光电器件领域备受关注，但 ABACUS+PyATB 的光学性质第一性原理计算流程繁琐，手动配置环境、串联自洽计算与后处理步骤门槛高。本文借助 AsterFire Go 平台 ABACUS 智算 Agent，仅用一句自然语言指令，自动完成全套 DFT 计算，高效解析单层 GeS 光学各向异性，清晰揭示其面内、面外光学响应差异，为二维半导体光电性能快速模拟提供全新对话式科研方案。
AsterFire Go 是一个用聊天就能做科学计算的平台。药物设计、蛋白模拟、分子动力学、量子化学、DFT 计算……打字告诉它你想算什么，最快速度给你结果。平台致力于让大家轻松调用各类算力，不用琢磨环境与集群配置、纠结方法、调试参数，把精力留给科研本身。现阶段封装了 SPONGE（对标 GROMACS）、ABACUS（对标 VASP）、ORCA（学术用途）、ProteinMPNN、DSDP（对标 Autodock Vina）、Bindcraft、RFdiffusion 等软件，同时还有文献检索总结、pdf 高保真翻译等百余种工具，均可通过自然语言调用或串入工作流。
AsterFire Go 用户交流群
产品内测期间加入官方用户群可得 1 个月专业版会员，含 4 个流动星币 + 50 GB 存储空间 + 百余种高级模型 + 全算力消费 9 折优惠。
对初次接触第一性原理计算的用户而言，操作难点不仅在于某一个参数的设置，更在于把各阶段正确衔接起来。对于光学计算，自洽算完还要把 HR / rR / SR 等矩阵元文件（H(R)：实空间哈密顿量；S(R)：交叠矩阵；r(R)：位置矩阵）正确交付 PyATB，才能获得最终结果。本例展示的重点，就是用简短的自然语言实现整个目标，让任务自动连续推进，直接展示核心指标与最终结果。
| 传统 DFT 光学性质计算 | AsterFire Go 平台 | |
|---|---|---|
| 操作方式 | 手动准备结构、计算输入与后处理输入，命令行操作 | 以自然语言启动计算，自动规划，对话式操作 |
| 环境部署 | 手动部署 ABACUS 与 PyATB 环境：依赖、并行、求解器、后处理库等等 | 依托平台强大的科学计算综合能力，无需本地环境配置 |
| 文件准备 | 手动准备 STRU、INPUT、KPT、赝势与数值原子轨道文件，还要另写 PyATB 的输入 | 自动完成结构识别，以及输入文件的准备与校验 |
| 参数配置 | 自行测试收敛性，决定基组、截断能、K 点、占据带数、展宽 η 与能量步长等参数 | Agent 自动完成可靠且精确的参数配置 |
| 任务调度 | 自洽计算与光学后处理之间手动搬运 HR / rR / SR 矩阵文件，手动串联各阶段 | 自动任务调度，由 Agent 组织任务、传递文件，衔接各阶段结果 |
| 结果整理 | 手动从输出文件提取 ε1、ε2，再换算出 α、n、k、R 等并绘图 | 自动完成结果整理 |
| 以本案例为例 | 需分步手动实施计算 | 一句话完成全套计算 |
ABACUS（Atomic-orbital Based Ab-initio Computation at USTC，中文名：原子算筹）是一款国产开源的以 DFT 计算为主的第一性原理软件，由中国科学技术大学何力新教授、中科院物理研究所任新国研究员和北京大学陈默涵研究员主导开发，拥有完全自主知识产权。其主要开发人员来自中国科学技术大学、中国科学院物理研究所、北京大学、北京科学智能研究院和合肥综合性国家科学中心人工智能研究院等单位。ABACUS 支持平面波和数值原子轨道两种基组，对各种体系与性质的计算支持良好。
星使智算结合 ABACUS 开发者的专家知识，以及官方手册，各领域知识，和大量用户的实践反馈，打造了 ABACUS 智能体。通过它，对话即可完成 ABACUS 计算，整个过程中用户无需手动输入关键参数。
进入 ABACUS 智算 Agent，上传已经优化好的 GeS 结构文件，并输入本次使用的指令：
计算优化后的 GeS 四原子二维结构的光学性质
这条指令同时说明了研究对象（GeS 四原子二维结构）、前置状态（结构已经优化过）和预期输出（光学性质）。进行其他计算时，也应说清这几样。若对晶型、真空层厚度等已有明确、特殊的要求，应一并写入。
Agent 识别、校验完结构后，会列出计算设置、平台资源探测结果、作业流水线规划等（此处为 结构输入 → ABACUS 计算 → 后处理）。用户确认上述内容后，Agent 开始执行计算。本次任务的要点包括：不重复做结构弛豫；调用 ABACUS（LCAO 基组）+ PyATB，PBE 泛函（默认），非自旋极化、无 SOC；输出为光吸收系数谱（cm-1，0-30 eV）与介电函数实部 / 虚部，光谱展宽 η = 0.1 eV。用户审核 Agent 的计划时，若有其他计算要求，也可随时补充说明。
对于不熟 ABACUS 的用户，核对计划时，只需重点看「算什么」和「用什么结构算」，而不必逐项复述软件参数，智能体可以给出很好的计算设置。
图 1：描述 GeS 光学性质计算需求并确认任务计划。
确认后，Agent 开始准备输入文件并调用平台资源执行计算。本次任务，Agent 选择使用 1 张 V100 GPU 和 cuSolver 求解器。Agent 会根据任务，选用合适的高性能 GPU 计算资源，以及相应的求解器，以高效完成计算。
在整个执行过程中，Agent 会自动进行电子步监测、节点状态核查、输出合理性等一系列必要检查。例如本例中，它通过平台的状态核查判定 abacus_optical 工具处于「计算步骤运行中」，并以真实的 Slurm 作业号作为启动证据；同时确认：等待输入的节点为空、没有失败节点。如遇电子步不收敛等情况，Agent 会自动调整参数并重新计算（失败节点按调试循环，最多重跑 3 轮）——在运维层面上，绝大多数问题同样不需要用户手动解决。与此同时，用户可以在平台上预览各文件，无需下载，随时查看计算真实的运行状态。
本例中，ABACUS 自洽计算 20 步收敛，随后的 PyATB 光学后处理扫描 8000 个 k 点。整个工作流共耗时 17 分 23 秒。
任务结束后，平台返回「GeS 二维单层光学性质计算（ABACUS）· 执行报告」，内附吸收系数谱、介电函数实部 / 虚部、光学张量三分量图谱与汇总的原始数据。本次执行的关键指标为：光子能量采样区间 0-30 eV、共 3001 个采样点，吸收峰能量 6.85 eV，峰值吸收系数 1.02×106 cm-1，费米能级 -1.1187 eV，占据能带 20 条。同时，任务面板把本次执行的全部产物整理为 14 项关键文件，可一次下载。
图 2：执行计算、查看报告与结果。
自洽计算中，电子步 20 步收敛，总能为 -4453.3654 eV，费米能级为 -1.1187 eV；单点计算给出的最大残余力约 0.14 eV/Å、最大应力约 1.8 kbar（总压 1.14 kbar）。这说明，该二维单层构型合理，可以直接作为光学性质计算的起点。
光学性质对结构的敏感程度高于总能：如果用户误将未优化结构上传给智能体，吸收峰的位置与强度都会随之漂移，因此，在自洽结束后，Agent 额外核验了结构的可靠性（无用户提醒）。
图 3：GeS 单层的晶体结构、光吸收系数谱、介电函数实部，以及折射率与消光系数。
吸收谱给出沿三个方向的吸收系数。最高的吸收峰位于 6.85 eV，峰值吸收系数约 1.02×106 cm-1，这个峰主要来自面外（zz）分量；面内两个分量的主峰则分别落在 5.21 eV（约 7.26×105 cm-1，xx）与 4.47 eV（约 6.58×105 cm-1，yy），三分量的平均吸收峰在 6.40 eV（约 7.35×105 cm-1）。整体来看，吸收主要集中在 2-8 eV 这一紫外 - 可见光区间，在这一波段，GeS 的透过性能较差。
吸收系数本身直接由介电函数计算得到，这两者是对于线性光学响应的两种不同方面的表述。介电函数实部 ε1 描述介质对电场的极化响应：低频极限则给出光学介电常数（ε∞）的三个分量，对于 GeS，ε∞ 分别为 3.43（xx）、3.55（yy）与 2.23（zz），面内方向明显大于面外方向。随光子能量升高，ε1 依次在 5.07 eV（xx）、4.47 eV（yy）与 6.73 eV（zz）附近越过零转为负值——这三个零点恰好落在各自的吸收峰上。介电函数虚部 ε2 描述能量的耗散，其振幅在 3.57 eV（xx）、2.88 eV（yy）与 6.30 eV（zz）附近出现极大值，与吸收峰一一呼应。折射率 n 与消光系数 k 同理：面内的 n 在 2.6–2.8 eV 之间取到最大值，面外分量则要到 6 eV 前后才取到峰值；k 的极大值出现在各自吸收峰附近——面内约 1.38（xx）与 1.45（yy），面外约 1.52（zz）。
另外，还需提及，此处光学性质仅包含电子贡献，因此不能给出正确的中、远红外光学性质。对于这一波段，须额外考虑声子对光学性质的修正（相关计算非常耗时）。
把三个分量放在一起对比，二维 GeS 的光学各向异性一目了然：面内的 xx 与 yy 分量在 3-5 eV 形成第一个吸收波段，面外的 zz 分量则在这一区间响应较弱（4 eV 处约为面内的七分之一），而在 6.5-7 eV 出现整条谱上最强的吸收峰。反射率 R 也有类似趋势——低能区面内两个方向明显高于面外，而到 8 eV 附近，面外方向升至 0.42，面内则落回至 0.2 以下。
这种「面内低能活跃、面外高能活跃」的差别，根源在于单层褶皱结构，面内与面外的成键方式不同，因此使跃迁矩阵元对偏振方向敏感。对器件应用而言，这意味着可以用入射光的偏振方向来选择性激发某一类跃迁。
图 4：GeS 单层的面内与面外吸收对比、平均吸收系数、介电函数虚部与反射率。
回看整个任务，我们仅用一句话完成了整套计算。在操作上，用户主要参与需求描述和计划确认；而输入文件的准备、自洽计算、后处理、两者的衔接、结果汇总等等，全部由 Agent 处理。对刚开始接触二维材料光学性质计算的研究者来说，这种方式有助于先建立对 ABACUS 与 PyATB 联合流程，整体的、系统的认识。之后若有兴趣，研究者可以再深入研究相关理论，逐步了解计算方法，与各种参数的具体影响。
（小提示：也可以与 ABACUS Agent 对话学习相关知识哦~密度泛函的相关理论，已经以知识库的形式封装在 Agent 中啦）
任务结束后，Agent 会推荐下一步的计算方向，如用 HSE 得到更准确的带隙与吸收边、计算激子效应（BSE）后的光学吸收谱、计算声子谱与电声耦合，等等。用户如有需要，可再用简单几个字开启更为深入、复杂的计算分析。